Genomic Characteristics and E Protein Bioinformatics Analysis of JEV Isolates from South China from 2011 to 2018

Japanese encephalitis is a mosquito-borne zoonotic epidemic caused by the Japanese encephalitis virus (JEV). JEV is not only the leading cause of Asian viral encephalitis, but also one of the leading causes of viral encephalitis worldwide. To understand the genetic evolution and E protein characteristics of JEV, 263 suspected porcine JE samples collected from South China from 2011 to 2018 were inspected. It was found that 78 aborted porcine fetuses were JEV-nucleic-acid-positive, with a positive rate of 29.7%. Furthermore, four JEV variants were isolated from JEV-nucleic-acid-positive materials, namely, CH/GD2011/2011, CH/GD2014/2014, CH/GD2015/2015, and CH/GD2018/2018. The cell culture and virus titer determination of four JEV isolates showed that four JEV isolates could proliferate stably in Vero cells, and the virus titer was as high as 108.5 TCID 50/mL. The whole-genome sequences of four JEV isolates were sequenced. Based on the phylogenetic analysis of the JEV E gene and whole genome, it was found that CH/GD2011/2011 and CH/GD2015/2015 belonged to the GIII type, while CH/GD2014/2014 and CH/GD2018/2018 belonged to the GI type, which was significantly different from that of the JEV classical strain CH/BJ-1/1995. Bioinformatics tools were used to analyze the E protein phosphorylation site, glycosylation site, B cell antigen epitope, and modeled 3D structures of E protein in four JEV isolates. The analysis of the prevalence of JEV and the biological function of E protein can provide a theoretical basis for the prevention and control of JEV and the design of antiviral drugs.


Introduction
Japanese encephalitis (JE) is an acute zoonotic disease caused by the Japanese encephalitis virus (Japanese encephalitis virus, JEV) [1,2]. JEV is a member of the genus Favivirus of the arboviridae (Faviviridae). Pigs are the primary hosts and source of infection of JEV, and other domestic animals such as sheep, cattle, horses, mules, and dogs are other primary sources of infection or act as storage hosts. JEV causes significant economic losses to the pig industry, but also causes severe viral encephalitis in human beings [3]. It is reported that the fatality rate of JEV infection can reach 30%, and nearly half of the survivors are reported to have permanent neuropsychiatric sequelae [4,5]. JEV has long been considered the leading cause of Asian viral encephalitis, of which about 50 percent occurs in China. The occurrence of JE is closely related to the environmental climate (temperature, humidity) and mosquito density, and it is most common in rainy areas [6]. Compared with the north, this is also an important reason why the distribution of Japanese encephalitis is concentrated in the south of China [7].

Viruses and Cells
The vaccine strain CH/SA14-14-2 MSV/2018 (Genbank: MH258849.1), Vero cells and JEV polyclonal antibody used in the study were preserved by the Department of Microbiology and Immunology, School of Veterinary Medicine, South China Agricultural University.

Identification of Clinical Materials and Virus Isolation
From 2011 to 2018, 263 samples of pigs suspected of JE were collected from several pig farms in Guangdong, Guangxi and Hainan provinces. The diseased materials were placed in a precooled milk bowl. Next, 3.0 mL serum-free MEM cell culture medium was added, ground and mixed repeatedly, then frozen and thawed three times. The supernatant was centrifuged and filtered with a 0.22 µM filter to remove bacteria and stored at −80 • C. According to the instructions of the Takara virus DNA/RNA extraction kit, virus RNA was extracted, and JEV cDNA synthesis was carried out with the instructions  (Table 1) were used for the amplification of JEV genome. The DNA reactants of JEV were purified according to the instructions of the JEV gel purification kit made by the OMEGA Company. The purified product was sent to Shanghai Bioengineering Co., Ltd. for sequencing. The positive samples were filtered to remove bacteria, inoculated to grow into monolayer Vero cells, and incubated at 37 • C for one h. The cell supernatant was discarded; then, the samples were supplemented with 2% FBS MEM at 37 • C to continue the culture, and were observed day by day. When the cells showed pathological changes, the cell cultures were harvested in time and identified by RT-PCR. If there were no pathological changes, the cell cultures were gathered on the 5th day of culture, frozen and thawed repeatedly three times and stored at −80 • C. Table 1. Primers for JEV gene amplification.

Serological Identification, Subculture and Determination of Virus Content
The cytopathic change culture identified by PCR was mixed with JEV polyclonal antibody and incubated at 37 • C for one h; then, Vero cells were inoculated and cultured at 37 • C 5% CO 2 . A control group was set up and observed day by day. After plaque purification, the JEV-positive isolates were diluted with a serum-free MEM culture medium ten times. The diluted virus solution was inoculated into 96-well plate Vero cells. After incubation for 2 h, the supernatant was discarded. After cleaning with MEM cell culture medium, 100 µL of MEM culture medium containing 2% FBS was added to each well, and a cell blank control group was set up. Then, 5% CO2 was cultured at 37 • C for seven days and observed daily. The number of cytopathic pores was recorded, and the virus titer of porcine JEV was calculated by the Reed-Muench method [23].

Analysis of Genomic Characteristics and Genetic Variation of Amino Acids
To understand the genetic evolution characteristics of the whole genome of JEV isolates, the sequences of 40 domestic and foreign JEV strains and classical strains were collected from NCBI (https://www.ncbi.nlm.nih.gov/) (accessed on 20 April 2021) (Supplementary Table S1). The obtained gene sequences were spliced with the SeqMan program on DNAstar software (Headquarters: Madison, WI, USA), and the spliced sequences and reference strains were compared with MAFFT software (Research Institute for Microbial Diseases, Osaka University, Suita, Japan) [24]. Then, the genomes and E genes were analyzed by SDTv1.2. We used ClustalW v2.0 (Kyoto University Bioinformatics Center, Kyoto, Japan) to compare nucleotides and deduce amino acid sequences, and edited them manually using MEGA X (Department of Biological Sciences, Tokyo Metropolitan University, Hachioji, Japan) [25]. Using the maximum compound likelihood nucleotide substitution model realized by the adjacent join (NJ) method in MEGAX, the phylogenetic analysis of the complete coding region of nucleotide sequences (nt) and each protein gene data set was carried out. One thousand bootstrap copies were used to evaluate the robustness of the system diagram.

Bioinformatics Analysis
The N-glycosylation sites of JEV isolates were predicted by NetNGlyc4.0Server (https: //www.cbs.dtu.dk/services/NetNGlyc/) (accessed on 20 February 2022), and the phosphorylation sites of the proteins of JEV isolates and reference strains were indicated by NetPhos 3.1Server (https://www.cbs.dtu.dk/services/NetPhos/) (accessed on 20 February 2022). The T cell epitopes of JEV isolates and vaccine reference strains were predicted by the IEDB website (http://www.iedb.org/) (accessed on 22 February 2022) using SYFPEITHI, the database of MHC ligands and peptide motifs, and the antigenicity and B cell epitopes of JEV isolates were analyzed by Protean tool on DNAstar software and ABCpred prediction (http://crdd.Osdd.Net/raghava/abcpred/ABC_submission.Html) (accessed on 22 February 2022). The amino acid sequences of E protein of four JEV isolates and the vaccine strain CH/SA14-14-2 MSV/2018 were compared by a blast to select the optimal template PDB number (more than 60% homology). Then, Modeller software (version 10.3, California Institute for Quantitative Biomedical Research, USA)was used for multi-template modeling [26], and finally, this was visually displayed by PyMol software (Company: DeLano Scientific LLC, San Carlos, CA, USA) [27].

Clinical Sample Detection and Virus Isolation and Identification
From 2011 to 2018, 263 samples were collected from 14 pig farms in the Guangdong and Guangxi provinces for JEV testing (Figure 1). Among the 263 samples, 78 were positive for the porcine JEV nucleic acid test, with a positive rate of 29.7% (78/263). After the positive samples were treated and used for virus isolation on Vero cells, four JEV isolates were obtained, namely, CH/GD2011/2011, CH/GD2014/2014, CH/GD2015/2015, and CH/GD2018/2018. Four strains of cDNA were used as templates for PCR identification with JEV-specific primers, and the results were amplified bands of about 437bp ( Figure 2). Variants CH/GD2011/2011, CH/GD2014/2014, CH/GD2015/2015, and CH/GD2018/2018 were identified by porcine Japanese-encephalitis-specific positive serum. The results showed that CPE was not found in the neutralization group, serum control group, or the cell control group of the four JEV isolates, while CPE was found in the virus control group. The pathological changes were characterized by aggregation, enlargement, separation, and syncytial formation ( Figure 3). The results showed that all four isolates could be neutralized by the specific serum of porcine Japanese encephalitis, and all four isolates were porcine JEV. Four porcine JEV isolates were continuously cultured on Vero cells for four generations, and the F4 generation of the isolated strain was cloned and purified three times. The virus content of the purified strain was continuously passaged on Vero cells for the third time (Supplementary Table S2). The results showed that when all the isolated strains were transferred to F10 generation, the virus propagated stably, and the virus content was 10 7.5~1 0 8.7 TCID 50/mL. Using four JEV isolates of cDNA as the template and ddH 2 O as the control, PCR amplification was carried out using JEV whole-gene amplification primers (Supplementary Figure S1). The corresponding bands of four isolates were amplified by JEV complete gene amplification primers, consistent with the expected size. All the amplified products were purified, sent to Shanghai Bioengineering Technology Service Co., Ltd. for gene sequencing, and then spliced with SeqMan.

Nucleotide and Amino Acid Similarity and Genetic Evolution Analysis of JEV Isolates
The genome and E gene sequences of JEV isolates and different genotype reference strains were analyzed by the MAFFT algorithm in MAFFT software, and the genomic phylogeny was analyzed by MEGA X software. Finally, the NJ phylogenetic tree based on the JEV genome and E gene was obtained ( Figure 4). The results showed that the four isolates and reference strains were divided into five clusters regarding genomic and genetic evolution: GI, GII, GIII, GIV and GV. Variants CH/GD2011/2011 and CH/GD2015/2015 were located in the GIII cluster, and their genetic relationships were close to the CH/SC2016/2016 strain and the CH/SA14-14-2 MSV/2018 strain. Strains CH/GD2014/2014 and CH/GD2018/2018 were located in the GI cluster. Their genetic relationship was closer to the CH/SD0810/2008 strain and the JPN/Mo-Kagawa/2020 strain. Four JEV isolates were more closely related to the GII, GIV and GV JEV classical strains (Figure 4a). The NJ phylogenetic tree based on the E gene of four JEV isolates and reference strains showed that CH/GD2011/2011 and CH/GD2015/2015 variants were also located in the GIII cluster. The genetic relationship was closer to the CH/Hubei/2015 strain and CH/SA14-14-2 MSV/2018 strain, while CH/GD2014/2014 and CH/GD2018/2018 variants were also located in the GI cluster. Their genetic relationship was closer to the CH/SD0810/2008 strain and CH/anheal/2017 strain (Figure 4b).

Nucleotide and Amino Acid Similarity and Genetic Evolution Analysis of JEV Isolates
The genome and E gene sequences of JEV isolates and different genotype reference strains were analyzed by the MAFFT algorithm in MAFFT software, and the genomic phylogeny was analyzed by MEGA X software. Finally, the NJ phylogenetic tree based on the JEV genome and E gene was obtained ( Figure 4). The results showed that the four isolates and reference strains were divided into five clusters regarding genomic and genetic evolution: GI, GII, GIII, GIV and GV. Variants CH/GD2011/2011 and CH/GD2015/2015 were located in the GIII cluster, and their genetic relationships were close to the CH/SC2016/2016 strain and the CH/SA14-14-2 MSV/2018 strain. Strains CH/GD2014/2014 and CH/GD2018/2018 were located in the GI cluster. Their genetic relationship was closer to the CH/SD0810/2008 strain and the JPN/Mo-Kagawa/2020 strain. Four JEV isolates were more closely related to the GII, GIV and GV JEV classical strains ( Figure 4a). The NJ phylogenetic tree based on the E gene of four JEV isolates and reference strains showed that CH/GD2011/2011 and CH/GD2015/2015 variants were also located in the GIII cluster. The genetic relationship was closer to the CH/Hubei/2015 strain and CH/SA14-14-2 MSV/2018 strain, while CH/GD2014/2014 and CH/GD2018/2018 variants were also located in the GI cluster. Their genetic relationship was closer to the CH/SD0810/2008 strain and CH/anheal/2017 strain ( Figure 4b).  The genomic nucleotides and deduced amino acid sequences of four JEV isolates and GI, GII, GIII, GIV and GV JEV classical strains were analyzed. The results showed that the nucleotide similarity and deduced amino acid sequence similarity of the E gene of the four JEV isolates and reference strains were 76.2~99.9% and 89.6~100.0%, respectively. In contrast, the genomic nucleotide sequence similarity was 77.9~99.2%. (Figure 5). At the The genomic nucleotides and deduced amino acid sequences of four JEV isolates and GI, GII, GIII, GIV and GV JEV classical strains were analyzed. The results showed that the nucleotide similarity and deduced amino acid sequence similarity of the E gene of the four JEV isolates and reference strains were 76.2~99.9% and 89.6~100.0%, respectively. In contrast, the genomic nucleotide sequence similarity was 77.9~99.2%. (Figure 5). At the same time, the nucleotides and deduced amino acid sequences of C, PrM, NS1, NS2a, NS2b, NS3, NS4a, NS4b and NS5 genes of the isolate were 75.8%~100.0% and 78.9%~100.0%, respectively (Supplementary Tables S3 and S4). The genomic nucleotides and deduced amino acid sequences of four JEV isolates and GI, GII, GIII, GIV and GV JEV classical strains were analyzed. The results showed that the nucleotide similarity and deduced amino acid sequence similarity of the E gene of the four JEV isolates and reference strains were 76.2~99.9% and 89.6~100.0%, respectively. In contrast, the genomic nucleotide sequence similarity was 77.9~99.2%. (Figure 5). At the same time, the nucleotides and deduced amino acid sequences of C, PrM, NS1, NS2a, NS2b, NS3, NS4a, NS4b and NS5 genes of the isolate were 75.8%~100.0% and 78.9%~100.0%, respectively (Supplement Tables S3 and S4). The protein encoded by the JEV E gene is the main structural protein and antigen protein of JEV which contains the core antigen determinant and is closely related to the virus' adsorption, entry and pathogenicity [28]. The amino acid sequence analysis of E protein of four JEV isolates and reference strains showed that both CH/GD2011/2011 and CH/GD2015/2015 variants had the same amino acid mutations at positions 222, 327 and 366, which were S (serine) to A (alanine), T (threonine) to S (serine), and S (serine) to A (alanine), respectively. The mutation characteristics were similar to those of classical GIII strains such as CH/SA14-14-2/2018, indicating that the above amino acid sites were high-frequency mutations in the E protein of JEV GIII strains. Compared with the reference strain, the CH/GD2014/2014 variant and the CH/GD2018/2018 variant had no amino acid mutation, which indicated that the amino acid sequence of E protein of the JEV GI strain was higher than that of the JEV GI strain. The amino acid sequence of the E protein of the GIII strain was prone to mutation. The 107th amino acid of the E protein of the CH/GD2011/2011 variant mutated from L (leucine) to F (phenylalanine), the 176th amino acid from I (isoleucine) to V (valine), the 244th amino acid from E (glutamic acid) to G (glycine), and the 279th amino acid from K (lysine) to M (methionine). It is suggested that neurotoxicity may be weakened by the mutation of the above amino acid sites ( Figure 6).

Glycosylation Site Analysis
Glycosylation modification is a significant post-translational modification. Glycosylation modification affects the spatial conformation, activity, transport and localization of proteins and plays a vital role in signal transduction, molecular recognition and immunity. At the same time, glycosylation plays an essential role in virus replication, infection, immunity and virulence [29][30][31]. The N-glycosylation sites of the proteins encoded by the E gene of four JEV isolates were analyzed by NetNGlyc1.0Server. All four JEV E proteins had a potential N-glycosylation site at N154, and the confidence level reached 9/9. The glycosylation sites of E protein of the four JEV isolates were the same as the potential glycosylation sites of the CH/SA14-14-2MSV/2018 protein of the CH/BJ-1/1995 strain and the vaccine strain of JEV. The potential glycosylation sites of JEVE protein did not change (the results are not shown), suggesting that the glycosylation sites of the JEVE protein may be conserved in different genotypes and strains.
protein of four JEV isolates and reference strains showed that both CH/GD2011/2011 and CH/GD2015/2015 variants had the same amino acid mutations at positions 222, 327 and 366, which were S (serine) to A (alanine), T (threonine) to S (serine), and S (serine) to A (alanine), respectively. The mutation characteristics were similar to those of classical GIII strains such as CH/SA14-14-2/2018, indicating that the above amino acid sites were highfrequency mutations in the E protein of JEV GIII strains. Compared with the reference strain, the CH/GD2014/2014 variant and the CH/GD2018/2018 variant had no amino acid mutation, which indicated that the amino acid sequence of E protein of the JEV GI strain was higher than that of the JEV GI strain. The amino acid sequence of the E protein of the GIII strain was prone to mutation. The 107th amino acid of the E protein of the CH/GD2011/2011 variant mutated from L (leucine) to F (phenylalanine), the 176th amino acid from I (isoleucine) to V (valine), the 244th amino acid from E (glutamic acid) to G (glycine), and the 279th amino acid from K (lysine) to M (methionine). It is suggested that neurotoxicity may be weakened by the mutation of the above amino acid sites ( Figure 6).

Glycosylation Site Analysis
Glycosylation modification is a significant post-translational modification. Glycosylation modification affects the spatial conformation, activity, transport and localization of proteins and plays a vital role in signal transduction, molecular recognition

Phosphorylation Site Analysis
Protein phosphorylation is a type of post-translational modification of proteins. Phosphorylation is mainly concentrated on tyrosine, threonine, and serine residues, and an additional charge is added to the protein after phosphorylation. At the same time, the unique size and charge characteristics of covalently linked phosphates also allow phosphorylated proteins to specifically recognize phosphorylated proteins through phosphorylationspecific binding domains in other proteins, thus promoting inducible protein-protein interactions [32]. Biological software deduced the E gene sequences of four JEV isolates into amino acid sequences. The phosphorylation sites of E protein of four JEV isolates were predicted by NetPhos3.1Server. It was found that E protein had 72~79 potential phosphorylation  Tables S5 and S6). These phosphorylation sites suggest that the E protein of the JEV isolate may play an essential role in virus replication and immune regulation.  Table S7). The results showed that the E protein of four JEV isolates had good antigenicity and could stimulate the body to induce a humoral immune response (Figure 7). T cell epitopes refer to antigenic epitopes that can be recognized by T cell receptors and induce cellular immune responses. Using the existing database of MHC ligands and polypeptide motif SYFPEITHI, the critical value of the software prediction score was set to 25 to predict the T cell epitopes of the virus proteins of the four JEV isolates and vaccine strains, respectively. Prediction results: JEV CH/GD2011/2011 variant E protein may contain 11 potential candidate T cell epitopes; JEV CH/GD2015/2015 variant E protein may contain 12 potential candidate T cell epitopes; the JEV CH/GD2014/2014 and JEV CH/GD2018/2018 variants may contain 10 potential candidate T cell epitopes and 10 repeats with the JEV CH/GD2011/2011 variant; and the CH/GD2015/2015 variant has one more "GSQEGGLHLALAGAI" T cell epitope (Table 2).  T cell epitopes refer to antigenic epitopes that can be recognized by T cell receptors and induce cellular immune responses. Using the existing database of MHC ligands and polypeptide motif SYFPEITHI, the critical value of the software prediction score was set to 25 to predict the T cell epitopes of the virus proteins of the four JEV isolates and vaccine strains, respectively. Prediction results: JEV CH/GD2011/2011 variant E protein may contain 11 potential candidate T cell epitopes; JEV CH/GD2015/2015 variant E protein may contain 12 potential candidate T cell epitopes; the JEV CH/GD2014/2014 and JEV CH/GD2018/2018 variants may contain 10 potential candidate T cell epitopes and 10 repeats with the JEV CH/GD2011/2011 variant; and the CH/GD2015/2015 variant has one more "GSQEGGLHLALAGAI" T cell epitope (Table 2). , respectively. Moreover, 3p54 is the crystal structure of the E protein of the JEV vaccine strain obtained by X-ray diffraction. Finally, the subunit structure of the E protein of JEV CH/SA14-14-2MSV/2018 and four JEV isolates in their natural state was simulated and displayed by Modeller software. The tertiary structures of the above proteins were compared. The results showed only one amino acid difference in the tertiary structure between JEV CH/GD2011/2011 and CH/SA14-14-2MSV/2018. This variation site was at position 312, changing from K (lysine) to R (arginine). There was no significant difference in the overall structure between the two strains ( Figure 8). of the E protein of the JEV vaccine strain obtained by X-ray diffraction. Finally, the subunit structure of the E protein of JEV CH/SA14-14-2MSV/2018 and four JEV isolates in their natural state was simulated and displayed by Modeller software. The tertiary structures of the above proteins were compared. The results showed only one amino acid difference in the tertiary structure between JEV CH/GD2011/2011 and CH/SA14-14-2MSV/2018. This variation site was at position 312, changing from K (lysine) to R (arginine). There was no significant difference in the overall structure between the two strains ( Figure 8).  There were 11 amino acid differences between JEV CH/GD2015/2015 and CH/SA14-14-2MSV/2018 in the modeled 3D structures (Figure 9) and the variation sites were at the 107th, 129th, 138th, 176th~177th, 232, 244, 264, 279, 315 and 439 amino acids ( Figure 10): F (phenylalanine) to L (leucine), T (threonine) to S (serine), K (lysine) to E (glutamic acid), V (valine) to I (isoleucine), A (alanine) to T (threonine), A (alanine) to T (threonine), G (glycine) to E (glutamic acid), H (histidine) mutated to L (leucine), M (methionine) mutated to K (lysine), V (valine) mutated to A (alanine), and R (arginine) mutated to K (lysine).

Discussion
JE is an acute mosquito-borne infectious disease considered the leading cause of Asian viral encephalitis, of which about 50% occurs in China. Although the number of cases of JEV infection has decreased significantly since the launch of the national vaccination program in the 1970s, the number of reported cases remains high [6,33]. The prevalence of JE is related to climatic factors such as annual precipitation, and it is most common in southwestern and eastern China [7]. JEV can be replicated in mosquitoes and transmitted to pigs by hematophagy so that pigs can become the temporary host of JEV and the source of human infection [34]. In recent years, intensive pig farming has been accelerated, which may increase the probability of pigs being infected with JEV. In this study, 78 porcine JEV nucleic acids were detected in 263 samples collected from 14 different pig farms in southern China from 2011 to 2018, with a positive rate of 29.7%. Considering the limited number of pig farms collected in this survey and the limitations of collecting serum, brain tissue, and aborted stillbirth samples, this can still reflect that there is a certain degree of JEV natural infection in pigs in Guangdong, Guangxi, and other regions. Guangdong and Guangxi are located in the tropics and subtropics, with a warm and humid climate, high mosquito density, large number of domestic pigs, and natural conditions for JEV transmission and reproduction. Although there is no large-scale epidemic of Japanese encephalitis, its potential harm cannot be ignored.
Four strains of viruses were isolated from Vero cells to explore the biological characteristics of porcine JEV which were positive for JEV infection. Different cell culture and virus titer determination of four JEV isolates showed that four JEV isolates could proliferate stably in Vero cells, and the virus titer was as high as 10 8.5 TCID 50/mL.
To understand the characteristics of the genetic evolution of JEV isolates, we determined the whole genome sequences of four JEV isolates in this study. In this study, we performed the phylogenetic analysis and genetic development of JEV by sequencing the entire genome nucleotide sequences of JEV and the nucleotide and amino acid sequences of the E protein encoded by the E gene. JEV is classified into GI, GII, GIII, GIV, and GV types. The JEV variants isolated in this study, CH/GD2011/2011 and CH/GD2015/2015, were classified as GIII at the whole genome nucleotide sequence, E gene nucleotide and amino acid sequence levels and the CH/GD2014/2014 and CH/GD2018/2018 were classified as GI. This is compatible with other studies stating that the predominant circulating strains were GIII and GI [3]. The insertion versus deletion analysis results of the amino acid sequence of each protein of JEV also corroborated the above conclusions. E protein variations, as an essential virulence protein of JEV, and related amino acid sites significantly impact the neurotoxicity and invasiveness of JEV. The current study shows that the amino acid sites involved in the virulence of the E protein include amino acids 107, 138, 123, 176, 177, 244, 279, 315 and 439; however, only the CH/GD2011/2011 harbored mutations that attenuated its neurotoxicity at amino acids 107, 176, 244 and 279. Specific amino acid mutations in the E protein are critical for viral virulence, including E107 and E138, which are JEV virulence regulatory genes; residue E107 is located in a highly conserved hairpin motif in the extended domain II. In contrast, this motif region contains a fusion peptide, and the mutation of E-107 can alter the fusion properties of E protein in cell culture [35]. Zheng et al. [28] demonstrated that the acidity/basicity of E138 residue affects JEV neurovirulence, and the mutation of essential residue in E138 contributes to the attenuation of neurovirulence. In contrast, the mutation of acidic residues enhanced the neurovirulence of the virus. Yang et al. [15] investigated the role of five amino acid residues (E-107, E-138, E-176, E-177, E-279) in the E protein in weakly virulent strains compared with the virulent parental strain, and the findings showed that amino acids at positions E-107 and E-138 played a critical role in the attenuation of neurovirulence in the CH/SA14-14-2MSV/2018 strain, within which E-107 (L→F) and E-138 (E→K) attenuated the neurotoxicity of CH/SA14-14-2MSV/2018. Notably, in the present study, the CH/GD2011/2011 conformed to the above mutations of the E protein at residues E-107 and E-138, providing a reference for subsequent studies on the pathogenicity of JEV isolates. Interestingly, the CH/GD2015/2015 consisted of methionine to threonine mutation at E-129, commonly found in GI or GIII type I to serine, as the cause of the mutation, and whether new mutations occurred in JEV awaits further confirmation. E-176 and E-177 are located in the main domain I of the E protein, a region rich in conformational epitopes sensitive to low pH. The area around the E-176 and E-177 clusters are also considered a virulence locus [36]. Zhou et al. [37] by obtaining the SCYA 2010901 strain after passaging the SCYA 201201 strain 120 times, compared the amino acid sequence of SCYA 2010901 with those of other JEV strains and found that although a single E-176 (I→R) mutation did not affect the viral growth in BHK-21 cells, it significantly reduced the neurotoxicity of JEV in vivo, and the four JEV isolates in this study, namely, E-176 (I→V) and E-177 (T→A), had no mutations in the remaining three. Except for E-176 and E-177, residue E-244 (E→G) and residue E-279 (K→M) also had the same mutation pattern in CH/GD2011/2011 vs. CH/SA14-14-2MSV/2018. E-244 is located in extended domain II and is an essential site outside the critical surveillance site for JE vaccine safety, and E-279 is located in a hinge region in the extended domain II β chain, which may play a regulatory role in E protein function. The present study's analysis of amino acid genetic variations in the E protein revealed that the CH/GD2011/2011 strain, which shares most of the same amino acid mutation patterns as the vaccine strain, may become a successful vaccine strain for the prevention of JEV infection.
The N-glycosylation of the protein plays an essential role in protein folding and can obtain the correct spatial conformation of the protein. Zhang [38] confirmed that the mutation of the N-glycosylation site N154 in the JEV E protein significantly enhanced the induced humoral immune response. Compared with the classical strains, the E protein of the four JEV isolates in this study was more conservative. Dechtawewat et al. [39] used immunoprecipitation with specific antibodies and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to study the potential phosphorylation sites on the NS1 protein of dengue virus (DENV), revealing that DENV NS1 phosphorylation has a functional effect on viral amplification during DENV infection. These findings emphasize the importance of phosphorylation and can also be used as experimental evidence to verify the predicted phosphorylation sites in this study, which may pave the way for the future design of targeted specific antiviral drugs. The E protein is a vital antigen protein, and the results show that its phosphorylation site has high homology among different strains. It is suggested that the E protein of the four JEV isolates can perform its biological function typically. The changes in some E protein phosphorylation sites or N-glycosylation sites may be similar to the mutation of their amino acid sites, which can enhance or reduce the neurotoxicity of the virus strains.
Identifying antigenic epitopes provides a scientific basis for understanding disease aetiology, immune surveillance, developing diagnostic kits, and designing epitope-based subunit vaccines [40]. However, epitope identification is expensive and time-consuming because it requires the experimental screening of many potential epitope candidates. It is necessary to predict epitopes through computer prediction methods and tools for B and T cells. The software prediction found that the four JEV isolates had potential B cell epitopes at positions 24,39,183,198, 332, 433 and 455, respectively. Compared with JEV CH/SA14-14-2MSV/2018, both CH/GD2014/2014 and CH/GD2018/2018 had a potential B cell epitope "LSLPWTSPSSTAWRNR" at 221-236, which provided a new insight for the screening of the JEV antigen epitope vaccine. The E protein of the four isolates all had dozens of T cell epitopes. It is worth noting that the confidence values of "GNYSAQVGASQAAKF" and "FLATGGVLVFLATNV" of the epitopes starting at 153,484 were above 30. The results show that the epitopes can be used to develop diagnostic kits and subunit vaccines based on epitopes in JEV.
Genes directly control biological traits by exploiting the structure of proteins. The molecular structure of a protein is very complex and can be divided into primary, secondary, tertiary and quaternary structures. Its secondary system can be used to further predict the tertiary and quaternary structure and to speculate the function of proteins [41,42]. Protein structural homology modeling provides a structural model for life science research when no experimental structure is available. The actual structural model should not only correctly represent the overall folding of a single protein chain but also correctly represent the atomic details of the interaction with essential cofactors and ligands [43]. In this study, the existing JEVE protein template was selected to simulate the E protein structure of the isolated JEV strain to show its three-dimensional structure and mutation site. The interaction between fusion loop (DII) and fusion loop pocket (DI-DIII) is particularly important for maintaining the JEV E dimer. The confirmed structure reveals the interaction between flavivirus' fusion loop and fusion loop pocket, including the two amino acids in the fusion ring located in W101 and F108, respectively [44]. Interestingly, the W101 site of JEV E interacts with K312, N313, A315, and V323, while F108 interacts with A315 and D316 on the relative subunits [45]. In the amino acid mutation and tertiary structure of binding E protein, it was not difficult to see that the four JEV isolates did not mutate at E-101 or E-108. CH/GD2014/2014, CH/GD2015/2015, and CH/GD2018/2018 did not show mutations in E-312, E-313, E-315, E-316 or E-323. However, like CH/SA14-14-2MSV/2018, CH/GD2011/2011 mutated from A (alanine) to V (valine) in E-315. It is worth noting that CH/GD2011/2011 mutated from K (lysine) to R (arginine) in E-312.The protein modeled 3D structures of four JEV isolates was compared with that of the CH/SA14-14-2MSV/2018. Except that there is only one amino acid difference in modeled 3D structures between JEV CH/GD2011/2011 and CH/SA14-14-2MSV/2018: this variation site is at position 312, changing from K (lysine) to R (arginine). There are 11-13 amino acid differences between CH/GD2015/2015, CH/GD2014/2014, CH/GD2018/2018, and CH/SA14-14-2MSV/2018 with the evolution of the JEV virus. We do not know whether the existing effective JEV vaccine can effectively contain the spread of JEV in the future.
Through the prediction and analysis of glycosylation sites, phosphorylation sites, B cell epitopes, and T cell epitopes of the JEV E protein, as well as the study of the modeled 3D structures, we not only have a deeper understanding of the structure of JEV and E protein, but have also laid a foundation for the future study of the function and mechanisms of structural and nonstructural proteins. This can also lay the foundation for the research and design of more accurate and effective immunogenicity and antigenicity of the JE vaccine.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/vaccines10081303/s1. Table S1: Referenced JEV strains used in this experiment; Figure S1: Genomic amplification of four JEV isolates; Table S2: Virus content assay of pig JEV isolates; Table S3: Nucleotide sequence similarity analysis of the genes of JEV isolates; Table S4: Amino acid sequence similarity analysis of the gene of JEV isolates; Table S5: Common phosphorylation sites of E proteins of JEV isolates and their corresponding kinases; Table S6: Different phosphorylation sites of E proteins of JEV isolates and their corresponding kinases; Table S7: B-cell antigenic epitopes of JEV isolate E protein.