- Research article
- Open Access
- Open Peer Review
Envelope protein gene based molecular characterization of Japanese encephalitis virus clinical isolates from West Bengal, India: a comparative approach with respect to SA14-14-2 live attenuated vaccine strain
BMC Infectious Diseases volume 13, Article number: 368 (2013)
Increasing virulence of Japanese encephalitis virus (JEV), a mosquito-borne zoonotic pathogen is of grave concern because it causes a neurotrophic killer disease Japanese Encephalitis (JE) which, in turn, is responsible globally for viral acute encephalitis syndrome (AES). Despite the availability of vaccine, JE/AES cases and deaths have become regular features in the different rural districts of West Bengal (WB) state, India, indicating either the partial coverage of vaccine or the emergence of new strain of JEV. Therefore, a study was undertaken to characterize and compare the complete envelope (E) protein gene based molecular changes/patterns of JEVs circulating in WB.
Total of 98 AES case-patients’ samples were tested to detect the presence of JEV specific immunoglobulin M (IgM) antibody by Mac-ELISA method. Only JEV IgM negative samples with a history of ≤3 days’ illness were screened for virus isolation and RT-PCR. E gene sequences of JEV isolates were subjected to molecular phylogeny and immunoinformatics analysis.
Present study confirmed JEV etiology in 39.7% and 29.1% of patients presenting ≤15 days’ febrile illness, as determined by Mac-ELISA and RT-PCR respectively. Phylogenetic analysis based on complete E gene sequences of JEV isolates showed the co-circulation of JEV genotype I (GI) with genotype III (GIII). This study also demonstrated that isolate-specific crucial amino acid substitutions were closely related to neurovirulence/neuroinvasiveness of JE. On the basis of immunoinformatics analysis, some substitutions were predicted to disrupt T-cell epitope immunogenicity/antigenicity that might largely influence the outcome of vaccine derived from JEV GIII SA14-14-2 strain and this has been observed in a previously vaccinated boy with mild JE/AES due to JEV GI infection.
Based on molecular evolutionary and bioinformatic approaches, we report evolution of JEV at a local level. Such naturally occurring evolution is likely to affect the disease profile and the vaccine efficacy to protect against JEV GI may demand careful evaluation.
The flaviviruses persist as causative agents for a wide range of human related infectious diseases worldwide. However, the most common causes of Acute Encephalitis Syndrome (AES) are associated with the Japanese encephalitis virus (JEV) serocomplex that includes Murray Valley encephalitis virus (MVEV), West Nile virus (WNV), St. Louis encephalitis virus (SLEV), and prototypical member JEV within/belonging to the Flavivirus genus under the family Flaviviridae[1, 2]. The mosquito-borne JEV is the sole etiologic agent of Japanese encephalitis (JE); a neurotropic killer disease being one of the major causes of viral acute encephalitis in human. Since the isolation of this virus in Japan in 1935 , it has been detected worldwide becoming a major public health problem. Worldwide case-fatality rate of JE was recorded to be 30% approximately with 30-50% of survivors having irreversible neuropsychiatric sequelae .
The natural transmission cycle of JEV involves an enzoonotic (sylvatic) mosquito-bird-mosquito and/or mosquito-pig-mosquito cycle; primarily involving Culex spp. mosquitoes as primary vectors , wading birds as reservoir host , pigs as amplifying host  and Humans are the accidental “dead end” hosts .
Like other flaviviruses, JEV, an enveloped positive-sense single stranded RNA (~ 11 kb in length) virus contains single open reading frame (ORF) encoding a polyprotein that is processed into three structural (C, M, and E) and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) proteins, flanked by 5′- and 3′-non-translated regions (NTRs) [9–11]. Among the three structural proteins, the envelope (E) protein is considered as the most antigenic part of the viral genome and is found to be involved in the majority of the biological properties of the virus, such as binding to the cell receptor, inducing immunological responses (neutralization, passive protection and antibody dependent enhancement), virion assembly and fusion activity at low pH [12, 13]. In addition, the amino acid substitutions in E protein have a major role in determining the neuorovirulence or neuroinvasiveness . The nucleotide sequence of full-length E gene of JEV is an established/reliable phylogenetic marker because this region has got no selective pressure supporting obscure long-term evolutionary relationship. Based on the nucleotide sequence of E gene, JEV can be divided into five distinct genotypes . Mostly genotype III (GIII) is circulated in the Southeast Asian countries, including Japan, South Korea, China, Taiwan, Vietnam, Philippines, and India . However, it has been recently documented that GIII is replaced by genotype I (GI) in South Korea, Thailand and China .
In India, the existence of JEV was first reported serologically in 1954 . However, the disease was first recognized in India at Vellore (in the state of Tamil Nadu) in 1955 . Since then, epidemics of JE in different states have been recorded [19, 20]. It was mentioned that genotype III is predominant in India, but recently genotype I has been introduced in this country [16, 21].
The state of West Bengal (WB) is situated at the eastern part (23°00′N, 87°00′E) of India with an area and population of 88,752 km2 and 91,347,736 respectively bounded on the north by Sikkim and Bhutan, on the east by Assam and Bangladesh, on the south by the Bay of Bengal and on the west by Orissa, Bihar, Jharkhand and Nepal . In 1973 JE outbreak was first recorded in the rural districts of Burdwan and Bankura in the state of WB where 700 cases and 300 deaths were reported [23, 24]. Thereafter, several JE outbreaks took place in the state [25–27]. As per published literature, the State Health Department, Govt. of WB has conducted vaccination programme against JE in different rural districts of WB . But still sporadic JE/AES cases and deaths are being reported every year from the state . The people of the state dependent on cultivation work in the water-logged paddy fields which serve as congenial home for mosquito breeding . Moreover, to raise their economic status they usually take up piggery and mini-poultry in their own hut bringing animals in their close association. In addition, environmental factors of this state also favor JEV transmission . Moreover, the reports of JE incidences or endemicity of JE in the state might be the indications of partial vaccination although the emergence of mutated/new strain of JEV could not be excluded.
The present study was carried out to characterize the genetic variation in the E protein of JEV circulating in the WB state and to observe the difference and/or crucial amino acid changes in E protein of the WB isolates in comparison with SA14-14-2 live attenuated vaccine strain in order to enable us to understand the variations in antigenic and neutralizing properties of different JEV WB isolates, followed by analyzing the molecular impact of amino acid substitutions on envelope protein structure of those isolates in relation to disease severity.
Case enrollment and sample collection
Clinically diagnosed 98 AES cases (28 from Malda, 15 from Hooghly, 25 from Midnapore, 11 from Birbhum, 10 from Howrah and 9 from South 24 Parganas) of all age groups having high grade fever for ≤15 days with the symptoms viz. headache, vomiting, unconsciousness/coma, convulsion/seizure, abnormal movements, delirium, altered sensorium, neck rigidity, presence of kernig’s sign etc. were admitted in abovenoted 6 district government hospitals in the state of WB from August to December both in 2011 and 2012. The study was duly approved by the joint ethical committee of ICMR (Indian Council of Medical Research) virus unit and NICED (National Institute of Cholera and Enteric Diseases), Kolkata, India. Informed consents were obtained in prescribed proforma from the patients or legal guardians or relatives of the patients before the collection of the samples. A total of 98 [62 serum and 36 cerebrospinal fluid (CSF)] samples were collected and/or referred from the above said clinically diagnosed AES case-patients to ICMR Virus Unit, Kolkata, maintaining the cold chain, for the detection of JEV infection. Cerebral malaria and bacteriological etiology were ruled out on the basis of clinical observations by the physicians of the hospitals concerned. All sera and CSF samples were stored in aliquots at −80°C until use.
All the samples were tested to detect the presence of JEV specific immunoglobulin M (IgM) antibody using JEV IgM-capture ELISA test kit (supplied by National Institute of Virology, Pune, India) based on the manufacturers’ protocol. Optical density (O.D) was measured at 450 nm using an ELISA reader (Titertek Multiskan Plus, Lab systems Finland, Type-314).
Isolation of virus from samples
Only 24 JEV IgM negative samples with a history of ≤3 days’ illness were screened and 200 μl of each of them was inoculated to C6/36 cell line for virus isolation as previously described . The tissue culture fluids were collected from the samples producing prominent cytopathic effect (CPE), and centrifuged at 1000 × g for 5 minutes and the supernatants of tissue culture fluids (STF) were kept in aliquots at −80°c till the isolation/extraction of viral RNA, followed by E gene amplification through RT-PCR test.
RNA extraction, E gene amplification and nucleotide sequencing
Viral RNA was extracted from 140 μl of STF utilizing QIAamp RNA viral kit (Qiagen, GmbH, Hilden, Germany), following the manufacturer’s instructions. Extracted/purified RNA was applied as template for cDNA synthesis using avian myeloblastosis virus (AMV) reverse transcriptase (RT) and the JEV E gene specific reverse primer. In brief, the target RNA was converted to cDNA in 25 μl mixture volumes containing the following components: 11.76 μl of dH2O (MiliQ grade), 5 μl of 5× AMV RT buffer, 1.6 μl of 10 mM deoxynucleotide triphosphates (dNTPS) (Invitrogen, USA), 0.1 μl of 0.1 M dithiothreitol (DTT) (Promega, USA), 1.5 μl of 10 μM JEV E gene specific reverse primer , 0.04 μl of 10 U/μl AMV RT (Promega, USA) and 5 μl of RNA (50 pg to 1 μg). The reaction mixture was allowed to cDNA synthesis for 1 h at 42°c. The cDNA was subsequently used for polymerase chain reaction (PCR) amplification in the total of 20 μl reaction mixture that was prepared using 9.3 μl of dH2O (MiliQ grade), 2 μl of 10× DreamTaq DNA polymerase buffer, 0.5 μl of 10 mM deoxynucleotide triphosphates (dNTPS) (Invitrogen, USA), 1.5 μl of 10 μM JEV E gene specific primer pairs , 0.2 μl of 5U/μl DreamTaq DNA polymerase (Fermentas Inc., USA) and 5 μl of cDNA. PCR reaction conditions were as follows: initial denaturation at 94°c for 5 min, followed by 35 cycles of denaturation (at 94°c for 30 s), primer annealing (at 67°c for30 s), and primer extension (at 72°c for 1 min 30 s). Then a final extension step was carried out at 72°c for 5 min. The PCR products were separated by electrophoresis on 1% agarose gel, stained with ethidium bromide.
RT-PCR amplicons were purified with the help of the Qiagen gel extraction kit (Qiagen, GmbH, Hilden, Germany), according to the manufacturer’s protocol and subjected to direct sequencing using the BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, CA, USA), in accordance with manufacturer’s specifications. The products were purified by ethanol precipitation and were analyzed by an automated DNA sequencer, 3130XL Genetic Analyzer (PE Applied Biosystems, Foster city, CA, USA). The freely available Finch TV software (http://www.geospiza.com) was used to edit and correct the 1,500 nucleotides which generated complete E gene sequences of JEV isolates. Those sequences were subjected to BLAST search.
Multiple sequence alignment (MSA) and phylogenetic analysis
The E gene sequences of JEV strains/isolates used in MSA and phylogenetic analysis in this study were listed in Additional file 1. MSA and phylogenetic analysis were performed by CLUSTALW (http://www.ebi.ac.uk/Tools/clustalw2/index.html) and MEGA version 5.0 software (http://www.megasofteware.net). The phylogenetic tree was constructed by the neighbor-joining method, tested with Kimura 2-parameter model and evaluated by 1000 bootstrap pseudo replicates. The strain MVEV-1-51 was used as an out group for generating the rooted tree. Moreover, the structure-based MSA of JEV isolates with respect to SA14-14-2 vaccine strain were analyzed by an online server named Porter (http://www.distill.ucd.ie/porter/) to confirm the alteration of E protein secondary structure regarding disease severity.
Mapping of amino acid substitutions
Amino acid sequences were deduced from nucleotide sequences of the JEV isolates by Transeq (http://www.ebi.ac.uk/Tools/emboss/transeq/). The comparisons were made to identify the JEV isolates-specific amino acid substitutions in relation to JEV SA14-14-2 vaccine strain and were subsequently mapped onto the predicted 3-dimensional structures of isolate-specific E protein, as modelled based on predicted similarities to the homologous (95.8% - 97.8% identity) E protein of JEV SA14-14-2 vaccine strain [PDB:3P54] using ESyPred3D web server (http://www.fundp.ac.be/sciences/biologie/urbm/bioinfo/esypred/). Finally, the predicted JEV isolate-specific E protein structures were viewed and manipulated with the help of PyMol software (http://www.pymol.org/).
T-cell epitope prediction
T-cell epitopes of the JEV isolate-specific structural E proteins were predicted using the EpiJen online server (http://www.ddg-pharmfac.net/epijen/) with respect to JEV SA14-14-2 vaccine strain. The isolate-specific amino acid sequences of E protein were put to use as input to EpiJen including JEV SA14-14-2 vaccine strain and the server was run in the default mode. The appropriate proteasomal and tap cut-off values were selected during epitope predictions for the 18 different HLA alleles.
Serology, virus isolation and RT-PCR
Out of 98 samples, only 39 (39.7%) were reactive to JEV specific IgM antibody, of which 24 (61.5%) and 15 (38.4%) samples were CSF and serum respectively. Of the remaining 59, only 24 JEV IgM negative samples having the history of ≤3 days of febrile illness were screened and subsequently subjected to tissue culture for virus isolation. This resulted in 11 samples producing prominent CPE (characterized by the roundening of cells, increased granularity and vacuolation, followed by cell death and disruption of the monolayer by detachment of the dead cells) of which only 7 (29.1%) samples consisting of 5 (71.4%) from CSF and 2 (28.5%) from serum were determined to be JE with positive findings of viral RNA by RT-PCR test.
The basic characteristics of the patients from whom the samples were collected and selected/subjected to tissue culture for virus isolation, were presented in Table 1. The patients from the districts of Midnapore, Malda, Birbhum, Howrah and South 24 Parganas in WB state were all pediatric-adolescents, with an average age of 8–13 years whereas the patients belonging to Hooghly district were all adults, with an average age of 47 years. The summarized data of the patients from whom the viruses were isolated in the study was depicted in Table 2.
However, we have a total of 46 (39 IgM + 7 RT-PCR positive) JE cases representing 46.9% of 98 cases. Out of 46 cases, 33 (27 IgM + 6 RT-PCR positive) were pediatric-adolescents (71.7%) and remaining 13 (12 IgM + 1 RT-PCR positive) were found to be adult cases (28.2%). Moreover, the occurrence of JEV infection was recorded for the months of August to October in 2011 and 2012, with the maximum number of cases (71.4%) observed in September in 2 consecutive years.
Molecular phylogeny and sequence analysis of isolates
The Figure 1 represented the phylogenetic tree derived from 7 E gene sequences of JEV isolates along with 40 wild-type JEV strains, including 11 from India and 29 from worldwide (Additional file 1). Phylogram showed 6 E gene sequences of the isolates [GenBank:KC526869-KC526871 and KC802020-KC802022] belonging to GIII and comprising 98%-100% nucleotide similarity with each other and 93%-98% nucleotide similarity with other Indian GIII strains, having the highest similarity (97%-98%) with Indian P20778 strain [GenBank:Z34096] and those GIII E gene sequences also amounted to 95%-96% nucleotide similarity with GIII SA14-14-2 vaccine strain [GenBank:D90195]. On the other hand, 1 E gene sequence of another isolate [GenBank:KC526872] belonging to GI that was most similar (96%) with Japanese GI strain Ishikawa [GenBank:AB051292], followed by 94% nucleotide similarity with other Indian GI isolate JEV-GKP-0945054 [GenBank:HM156572] and also having 87% nucleotide similarity with above said vaccine strain. Above all, the E gene nucleotide sequences of GIII isolates comprised 89%-90% nucleotide identity with GI isolate.
The isolates IND/11/WB/JEV48 [GenBank:KC802020] and IND/11/WB/JEV49 [GenBank:KC802021] were identical (having 100% nucleotide similarity) with isolate IND/11/WB/JEV46 [GenBank:KC526869] and all of them showed same clinical outcome in the patient 4, 5 and 2 (Table 2). Moreover, 100% nucleotide similarity was also observed in between the isolate IND/12/WB/JEV51 [GenBank:KC802022] and IND/12/WB/JEV50 [GenBank:KC526871] and both of them were associated with mild JE in the patient 7 and 6 (Table 2). Since, the isolates IND/11/WB/JEV48, IND/11/WB/JEV49 and IND/12/WB/JEV51 were 100% identical with isolates IND/11/WB/JEV46 and IND/12/WB/JEV50, molecular modelling and T-cell epitope prediction data of the isolates IND/11/WB/JEV48, IND/11/WB/JEV49 and IND/12/WB/JEV51 were not shown in Figure 2 and/or in Tables 3 and 4 respectively except the other isolates.
Among several nucleotide changes, some were non-synonymous (data not shown) resulting in amino acid substitutions in JEV isolate-specific E protein in comparison with above said vaccine strain (Table 3). Moreover, the 8 amino acid substitutions (N2H, F107L, K138E, V176I, A177T, G244E, G261S and M279K) were common to all isolates whereas isolate-specific unique amino acid substitutions were also shown by a structure-based MSA of JEV isolates and SA14-14-2 vaccine strain [Additional file 2 (b)] which might confirm the alteration of E protein secondary structure with respect to disease severity.
Molecular modelling of isolates-specific E protein
Homology modelling of JEV isolates-specific 3-dimensional structures of E protein mapped with amino acid substitutions in relation to SA14-14-2 vaccine strain have been shown in Figure 2 (a)–(d). The amino acid substitutions were pointed out to areas with major secondary structures likely to affect the protein conformation regarding disease severity.
T-cell epitope prediction
Amino acid substitutions in E proteins of JEV isolates with potential T-cell epitopes were predicted in Table 4 and the substitutions in peptide-epitope sequences of the isolates with greater/lower IC50 values/scores were also given in Additional file 3 and Additional file 4 (actual Amino acid residues in peptide-epitope sequences of vaccine strain SA14-14-2 with IC50 values/scores were not shown). As a result, it is interesting to note that majority of the identified mutations falling within predicted T-cell peptide-epitopes might influence/alter immunogenicity/antigenicity in relation to said vaccine strain. While ranking the epitopes for mutation mapping, epitopes with calculated IC50 values/scores above 100 were not included since such peptide-epitopes were unlikely to yield significant immunogenicity.
The virulence of JEV etiology is of grave public health concern in the state of WB. The present study reveals that 39/98 (39.7%) and 7/24 (29.1%) samples were positive to JE by IgM-capture ELISA and RT-PCR method respectively. This observation proves that JE occurs recently and the ELISA negative acute samples with febrile illness for ≤ 3 days should be subjected to RT-PCR test to confirm the total number of JE cases.
Pediatric-adolescents (71.7%) were higher with JE cases than adults (28.2%) as because pediatrics were infected possibly due to lack of immunity and adolescents were directly exposed to the mosquito vector (Culex sp.) bite, as they usually took active part in cultivation in water-logged crop fields where vectors usually breed. In contrast, low numbers of JE cases in adults were possibly due to the development of immunity by sub-clinical infections in natural time scale. In the present study, the maximum number of JE cases was found to occur in the monsoon period i.e. in the month of September when the Culex mosquitoes breed in the paddy fields covered with stagnant rain water.
However, we proposed that 52 [(59–24) + (24–7)] samples with a history of 1–15 days’ illness were true JE negative possibly due to either mishandling of samples which damaged the IgM antibody/the viral titre or the presence of another etiology responsible for AES.
In molecular phylogeny, JEV GIII isolates i.e. IND/11/WB/JEV46 [GenBank:KC526869], IND/11/WB/JEV47 [GenBank:KC526870], IND/11/WB/JEV48 [GenBank:KC802020], IND/11/WB/JEV49 [GenBank:KC802021], IND/12/WB/JEV50 [GenBank:KC526871] and IND/12/WB/JEV51 [GenBank:KC802022] showed closer clustering, even though all of 6 isolates were from a geographically distant location or different districts of the state (Table 2, Figure 1). This indicated similarities among the isolates circulating in the 2 consecutive years 2011, 2012. Moreover, based on the present study, we stated that there had been simultaneous circulation of both JEV GI [GenBank:KC526872] and GIII [GenBank:KC526869] in the district of Midnapore, WB in 2011 corroborating with our earlier findings . It would be interesting to determine whether GI has become the predominant genotype or co-circulate with GIII in WB, India; additional JEV isolates from human clinical samples will be required to confirm the findings. Our study also reveals that the JEV GI [GenBank:KC526872] was isolated from a 9-year old boy-patient who had already been immunized with live attenuated JE vaccine derived from GIII strain SA14-14-2 and the boy was found to have clinically developed moderate JE with high fever including unconsciousness, neck rigidity and convulsion. In this connection, the efficacy of the vaccine to protect against GI isolate of JEV needs careful evaluation. In addition, we have identified the samples from 3 adult patients (average age of 47 years) with AES belonging to district Hooghly, WB were subjected to tissue culture for JEV isolation (Table 1) and out of these 3 samples, only 1 was confirmed as JE positive by RT-PCR (Table 2). This observation suggested that adults in this district were at higher risk for JE. Therefore, the state health department, government of WB should take initiative about JE vaccination to this population to overcome the disease burden.
The 500 amino acids long JEV E protein is encoded by E gene of 1500 nucleotides in length. The E protein of flaviviruses, including JEV, plays an important role in immunogenicity, tissue tropism, cell fusion and infection, and virus maturation . The predicted 3-dimensional structures of isolate-specific E protein were modelled onto the crystal structure of JEV SA14-14-2 vaccine strain E protein  and composed of separate structural domains [designated as domain I (DI), domain II (DII) and domain III (DIII)] including stem (ST) and transmembrane (TM) region like other flaviviruses [32, 33]. DI, which is referred to as the central domain, consists of 127 residues (1–51, 135–193 and 283–299); DII, which is called as dimerization domain, consists of 172 residues (52–134 and 194–282); the most antigenic part of E protein is DIII, a continuous stretch of 100 residues (300–399) that has been implicated in receptor binding; a stretch of 52 residues (400–451) is considered as ST region of E protein not simply to support the protein, but rather, it plays a key role in viral entry and maintenance of the protein complexes on the virion surface followed by TM, a 49 residues (452–500) long region playing a major role in the processing, sub-cellular localization and assembly of E protein [Additional file 2 (a)].
In the present study, in order to identify the possible correlation of clinical severity (i.e. neurovirulence/neuroinvasiveness) and altered immunogenicity of JEV infection in WB with isolate-specific amino acid changes of E protein, we analyzed E protein amino acid changes of the viral isolates with respect to SA14-14-2 vaccine strain. Our study showed greater amino acid changes in the JEV isolates circulating in WB than that of SA14-14-2 (Table 3). It was found that 8 critical amino acid changes at E107 (F→L), E138 (K→E), E176 (V→I), E177 (A→T), E264 (H→Q), E279 (M→K), E315 (V→A) and E439 (R→K) (Table 3) of 3 isolates i.e. IND/11/WB/JEV45-IND/11/WB/JEV47 (except the isolate IND/12/WB/JEV50) display neurovirulence in the patient 1–3 respectively (Table 2). That type of severity due to JEV infection worldwide or in earlier reports [14, 34–37] has been linked to the amino acid alterations in the virus. Moreover, the above said eight critical amino acid changes including E123 (S→R) and E306 (E→G) in the isolate IND/11/WB/JEV47 and only E306 (E→G) in the isolate IND/11/WB/JEV46 (Table 3) might introduce much more neurovirulence/neuroinvasiveness corroborating with previous findings [38, 39] and the increased virulence of the isolates is likely to cause death and severe JE/AES in the patient 3 and 2 respectively (Table 2). In contrast, although the isolate IND/12/WB/JEV50 showed critical amino acid changes, it did not include H→Q at E264, V→A at E315 and R→K at E439 which were found to be closely associated with neurovirulence (Table 3). Therefore, it might result in low neurovirulence in the patient 4 who developed mild JE/AES (Table 2). Moreover, regarding the low neurovirulence property of the isolate IND/12/WB/JEV50, it was speculated that a novel conservative amino acid substitution i.e. Q158P was mapped onto an α-helix of glycosylated E0F0 loop [31, 40] in DI of that isolate’s E protein [Figure 2 (a)] and was found to have potential to break (as because the proline is a helix breaker) the secondary structural elements surrounding that position and might also affect the N154 glycosylation site  of E protein of the isolate [Additional file 2 (b)], being less severe for the patient 4.
The crystal structure of E protein of WNV, a Flavivirus closely related to JEV , illustrated that the finger like projection of a stretch of 13 residues (98–110) long fusion loop (FL) is also present in the DII of JEV E protein, essential for viral infectivity, membrane fusion and is recognized as neutralizing epitope [31, 43]. In the study, Hopp and Woods hydrophilicity prediction  revealed that the non-conservative substitution of basic lysine (K) for polar asparagine (N) resulted in decreased (IC50 value/score) hydrophilicity of 98DRGWGKGCGLFKG110 peptide sequence of the isolate IND/11/WB/JEV45 E protein [Table 4, Additional file 2(b) and Additional file 3]. This in turn, indicates altered antigenicity of the isolate in comparison with SA14-14-2 vaccine strain. Moreover, we found that another non-conservative amino acid change i.e. E107 (F→L) at the tip of FL region of E protein of all JEV isolates circulating in WB might lead to escape from antibody neutralization or neutralizing epitope with decreased/increased IC50 value/score in relation to said vaccine strain [Table 4, Additional file 2(b), Additional file 3 and Additional file 4]. This finding was supported by previous studies [31, 45].
It is noteworthy that the amino acid substitutions at E64 (S→T), E244 (G→E), E306 (E→G) and E311 (A→R) in the DII and DIII of JEV isolate-specific E protein with respect to said vaccine strain [Table 4, Additional file 2(b), Additional file 3 and Additional file 4] might largely change immunogenicity (with decreased/increased IC50 value/score) of the isolates corroborating with earlier reports [33, 40, 46–49]. According to earlier workers [32, 40, 50], the surface exposed BC loop (E229-E333), DE loop (E365-368) and the putative receptor binding RGD motif (E387-E389), adjacent to FG loop (E389-E391) region were recognized as most antigenic part of DIII region of WNV E protein (closely related to JEV). In the present study, some amino acid substitutions (i.e. S327Q/T, S329T, S331K, A366S, G388E/K and D389E) were observed within and/or adjacent to BC loop, DE loop and RGD associated FG loop region of the JEV isolates’ (IND/11/WB/JEV45-IND/11/WB/JEV47) E protein (Table 3). In Hopp and Woods hydrophilicity prediction , both conservative and non-conservative substitutions of S327Q/T, S329T, S331K (polar to charged), A366S (non-polar to polar) and G388E/K (non-polar to charged) resulting in increased hydrophilicity of the 327 QYTGKDG333 / 327 TYTGSDG333 / 327 TYSGSDG333, 365SSNV368 and 387REDKQ391 / 387RKDKQ391 peptide sequences of above mentioned JEV isolates in relation to SA14-14-2 vaccine strain might lead to escape from antibody neutralization or neutralizing epitope due to increased/decreased IC50 value/score [Table 4, Additional file 2(b), Additional file 3 and Additional file 4] whereas another conservative amino acid change at E389 (D→E) in the BC loop region of E protein of the isolate IND/11/WB/JEV45 retained the same hydrophilicity of the 387RGEKQ391 peptide sequence in comparison with said vaccine strain.
Several mutations in the structural E protein gene region were seen to fall within the predicted T-cell epitopes (Table 4) as a result of analysis of the changes that might affect the immunogenicity of these isolates. Involvement of epitopes presented by major HLA alleles, such as HLA-A*02 and HLA-A*11 were shown by the present study. Moreover, HLA-A*11, one of the main HLA types, is found to be present in population across the world particularly with its high prevalence in South East Asia . When the variations in these epitopes take place, they may have the potential to affect the cell-mediated immune response against JEV infection.
It is worthy to mention that substitutions of non-polar to charged amino acid residues at E432 (G→R), E463 (L→R), E467 (L→E) and E481 (A→D) in the isolate-specific ST and TM region of E protein (Table 3) might facilitate the processing, sub-cellular localization and assembly of E protein of the isolates, which was proved by earlier study . Interestingly, the present JEV GI isolate IND/11/WB/JEV45 had a non-conservative amino acid substitution at E222 (A→S) corresponding to DII region of its E protein [Table 3, Figure 2 (d), Additional file 2 (b)] which was also observed in the previous study . Such type of mutation might influence/facilitate vector competency and viral infectivity. This was similar to the recent findings . All WB isolates showed amino acid mutations at E2 (Asn to His) and E261 (Gly to Ser) and the isolate-specific mutations i.e. N8H of IND/11/WB/JEV47; R416K of IND/11/WB/JEV47 and IND/12/WB/JEV50 and V372L, M374I and W396R of IND/11/WB/JEV45 were also observed in the study (Table 3), but the biological significance of those mutations were unknown.
In summary, our study reveals that JEV GI and GIII co-circulate with genetic changes in the state of WB, India during 2011–12. The JEV GI isolate IND/11/WB/JEV45 [GenBank:KC526872] from a 9-year old boy-patient who had been vaccinated with the live attenuated JE vaccine derived from JEV GIII SA14-14-2 strain but developed moderate JE/AES as because the GI isolate differed from the JEV GIII vaccine strain by means of some amino acid substitutions which were responsible for both neurovirulence and alteration of immunogenicity/antigenicity. Therefore, the vaccine efficacy to protect against JEV GI may demand careful evaluation. However, using a bioinformatics approach we point out that this study is a localized example of molecular evolution of E protein gene of JEV isolates in nature, which is likely to affect the disease profile. Further, functional studies on these mutant isolates would help us to understand the correlation of the genetic changes with respect to JEV virulence and pathogenesis.
The study was duly approved by the joint ethical committee of ICMR (Indian Council of Medical Research) virus unit and NICED (National Institute of Cholera and Enteric Diseases), Kolkata, India.
Lindenbach BD, Thiel HJ, Rice CM: Flaviviridae: the viruses and their replication. Fields virology. Volume 1. Edited by: Knipe DM, Howley PM. 2007, Philadelphia, PA: Lippincott Williams & Wilkins, 1101-1152. 4
Mackenzie JS, Gubler DJ, Petersen LR: Emerging flaviviruses: the spread and resurgence of Japanese encephalitis, West Nile and dengue viruses. Nat Med. 2004, 10 (suppl 12): S98-S109.
Tanaka M, Aira Y, Igarashi A: Comparative nucleotide and amino acid sequences of five Japanese encephalitis virus strains isolated in Japan and China. Trop Med. 1991, 33: 15-21.
World Health Organization: Japanese encephalitis vaccines-WHO position paper. Wkly Epidemiol Rec. 2006, 81: 331-340.
Geevarghese G, Kanogia PC, Mishra AC: Japanese encephalitis-Vector Biology. 2004, NIV Pune Year Book, Himayatnagar: Orient Longman Publication
The Center for Food Security and Public Health, Iowa State University: Japanese Encephalitis: Japanese B Encephalitis, Arbovirus B, Mosquito–borne Encephalitis Virus. http://www.cfsph.iastate.edu/Factsheets/pdfs/japanese_encephalitis.pdf,
Reuben R, Gajanana A: Japanese encephalitis in India. Indian J Pediatr. 1997, 64: 243-251. 10.1007/BF02752458.
Diagana M, Preux PM, Dumas M: Japanese encephalitis revisited. J Neurol Sci. 2007, 262: 165-170. 10.1016/j.jns.2007.06.041.
Chambers TJ, Hahn CS, Galler R, Rice CM: Flavivirus genome organisation, expression and replication. Annu Rev Microbio. 1990, 44: 649-688. 10.1146/annurev.mi.44.100190.003245.
Fields BN, Knipe DM, Howley PM: Fields Virology. 2007, Philadelphia, PA: Lippincott Williams & Wilkins
Sumiyoshi H, Mori C, Fuke I, Morita K, Kuhara S, Kondou J, Kikuchi Y, Nagamtu H, Igarashi A: Complete nucleotide sequence of the Japanese encephalitis virus genome RNA. Virology. 1987, 161: 497-510. 10.1016/0042-6822(87)90144-9.
Monath TP: Flaviviruses. Virology. Edited by: Fields BN. 1990, New York: Raven Press, 763-814.
Rice CM, Strauss EG, Strauss JH: Structure of the Flavivirus genome. The Togaviridae and Flaviviridae. Edited by: Schlesinger S, Schlesinger MJ. 1986, New York: Plenum Press, 279-326.
McMinn PC: The molecular basis of virulence of the encephalitogenic flaviviruses. J Gen Virol. 1997, 78: 2711-2722.
Solomon T, Ni H, Beasley DW, Ekkelenkamp M, Cardosa MJ, Barrett AD: Origin and evolution of Japanese encephalitis virus in Southeast Asia. J Virol. 2003, 77: 3091-3098. 10.1128/JVI.77.5.3091-3098.2003.
Fulmali PV, Sapkal GN, Athawale S, Gore MM, Mishra AC, Bondre VP: Introduction of Japanese encephalitis virus genotype I, India. Emerg Infect Dis. 2011, 17: 319-321. 10.3201/eid1702.100815.
Smithburn KC, Kerr JA, Gatne PB: Neutralising antibodies against certain viruses in the sera of residents of India. J Immunol. 1954, 72: 248-251.
Namachivayam V, Umayal K: Proceedings of National Conference. Japanese Encephalitis. 1982, New Delhi: Indian Council of Medical Research, 30-33.
Dhillon GP, Raina VK: Epidemiology of Japanese encephalitis in context with Indian scenario. Indian J Med Assoc. 2008, 106: 660-663.
Kabilan L, Rajendran R, Arunachalam N, Ramesh S, Srinivasan S, Samuel PP, Dash AP: Japanese encephalitis in India: an overview. Indian J Pediatr. 2004, 71: 609-615. 10.1007/BF02724120.
Sarkar A, Taraphdar D, Mukhopadhyay SK, Chakrabarti S, Chatterjee S: Molecular evidence for the occurrence of Japanese encephalitis virus genotype I and III infection associated with acute encephalitis in patients of West Bengal, India, 2010. Virol J. 2012, 9: 271-10.1186/1743-422X-9-271.
Sarkar A, Taraphdar D, Mukhopadhyay BB, Kumar M, Mukhopadhyay SK, Chatterjee S: Influence of socio-economic status and environmental factors on serologically diagnosed Japanese encephalitis cases in the state of West Bengal, India during 2005–2010. Health. 2012, 4: 6-12. 10.4236/health.2012.41002.
Ghosh SN, Rodrigues FM, Seth GP, Tongaonkar SS, Padbidri VS, Gupta NP: Investigations on the outbreak of Japanese encephalitis in Burdwan district, West Bengal. Part II. Serological survey of human population. Indian J Med Res. 1975, 63: 1472-1477.
Rodrigues FM, Ghosh SN, Banerjee K, Chatterjee AK, Gupta NP: A post-epidemic serological survey of humans in Bankura district, West Bengal, following the epidemic of Japanese encephalitis in 1973. Indian J Med Res. 1975, 63: 1478-1485.
Rajagopalan PK, Panicker KN: A note on the 1976 epidemic of Japanese encephalitis in Burdwan district West Bengal. Indian J Med Res. 1978, 68: 3938-
Banerjee K, Sengupta SN, Dandawate CN, Tongaonkar SS, Gupta NP: Virological and serological investigations of an epidemic of encephalitis which occurred at Bankura district, West Bengal. Indian J Med Res. 1976, 64: 121-130.
Mukhopadhyay BB, Mukherjee B, Bagchi SB, Chakraborty M, Mukherjee KK: An epidemiological investigation of Japanese encephalitis outbreak in Burdwan, District of west Bengal during 1987–1988. Indian J Public Health. 1990, 34: 107-116.
Sarkar A, Taraphdar D, Mukhopadhyay SK, Chakrabarti S, Chatterjee S: Serological and molecular diagnosis of Japanese encephalitis reveals an increasing public health problem in the state of West Bengal, India. Trans R Soc Trop Med Hyg. 2012, 106: 15-19. 10.1016/j.trstmh.2011.08.011.
Programme Implementation Plan 2011-12/National Rural Health Mission/Health and Family Welfare Department, West Bengal: Disease control programme: scheme/program under national rural health mission. http://pipnrhm-mohfw.nic.in/index_files/non_high_focus_large/westbengal/4d.pdf,
Westaway EG: Flavivirus replication strategy. Adv Virus Res. 1987, 33: 45-90.
Luca VC, Abimansour J, Nelson CA, Fremont DH: Crystal structure of the Japanese encephalitis virus envelope protein. J Virol. 2012, 86: 2337-2346. 10.1128/JVI.06072-11.
Rey FA, Heinz FX, Mandle C, Kunz C, Harrison SC: The envelope glycoprotein from tick-borne encephalitis virus at 2 Aº resolution. Nature. 1995, 375: 291-298. 10.1038/375291a0.
Kolaskar AS, Kulkarni-Kale U: Prediction of three-dimensional structure and mapping of conformational epitopes of envelope glycoprotein of Japanese encephalitis virus. Virology. 1999, 261: 31-42. 10.1006/viro.1999.9859.
Qi SC, Xiang ML, Qiao YZ, Dan DW, Huan CC, Ping Q: Isolation and molecular characterization of genotype 1 Japanese encephalitis virus, SX09S-01, from pigs in China. Virol J. 2011, 8: 472-10.1186/1743-422X-8-472.
Zhang JS, Zhao QM, Guo XF, Zuo SQ, Cheng JX, Jia N, Wu C, Dai PF, Zhao JY: Isolation and genetic characteristics of human genotype 1 Japanese encephalitis virus, China, 2009. PLoS ONE. 2011, 6: e16418-10.1371/journal.pone.0016418.
Zhao Z, Date T, Li Y, Kato T, Miyamoto M, Yasui K, Wakita T: Characterization of the E-138 (Glu/Lys) mutation in Japanese encephalitis virus by using a stable, full-length, infectious cDNA clone. J Gen Virol. 2005, 86: 2209-2220. 10.1099/vir.0.80638-0.
Thomas PM, Arroyo J, Levenbook I, Zhang ZX, Catalan J, Draper K, Guirakhoo F: Single mutation in the flavivirus envelope protein hinge region increases neurovirulence for mice and monkeys but decreases viscerotropism for monkeys: relevance to development and safety testing of live, attenuated vaccines. J Virol. 2002, 76: 1932-1943. 10.1128/JVI.76.4.1932-1943.2002.
Tajima S, Nerome R, Nukui Y, Kato F, Takasaki T, Kurane I: A single mutation in the Japanese encephalitis virus E protein (S123R) increases its growth rate in mouse neuroblastoma cells and its pathogenicity in mice. Virology. 2010, 396: 298-304. 10.1016/j.virol.2009.10.035.
Ni H, Barrett ADT: Molecular differences between wild-type Japanese encephalitis virus strains of high and low mouse neuroinvasiveness. J Gen Virol. 1996, 77: 1449-1455. 10.1099/0022-1317-77-7-1449.
Ryuta K, Kar K, Anthony K, Gould LH, Ledizet M, Fikrig E, Marasco WA, Koski RA, Yorgo M: Crystal structure of West Nile virus envelope glycoprotein reveals viral surface epitopes. J Virol. 2006, 80: 11000-11008. 10.1128/JVI.01735-06.
Beasley DW, Whiteman MC, Zhang S, Huang CY, Schneider BS, Smith DR, Gromowski GD, Higgs S, Kinney RM, Barrett AD: Envelope protein glycosylation status influences mouse neuroinvasion phenotype of genetic lineage 1 West Nile virus strains. J Virol. 2005, 79: 8339-8347. 10.1128/JVI.79.13.8339-8347.2005.
Grant EN, Nelson CA, Chen BR, Diamond MS, Fremont DH: Crystal structure of the West Nile virus envelope glycoprotein. J Virol. 2006, 80: 11467-11474. 10.1128/JVI.01125-06.
Sultana H, Foellmer HG, Neelakanta G, Oliphant T, Engle M, Ledizet M, Krishnan MN, Bonafé N, Anthony KG, Marasco WA, Kaplan P, Montgomery RR, Diamond MS, Koski RA, Fikrig E: Fusion loop peptide of the West Nile virus envelope protein is essential for pathogenesis and is recognized by a therapeutic cross-reactive human monoclonal antibody. J Immunol. 2009, 183: 650-660. 10.4049/jimmunol.0900093.
Hopp TP, Woods KR: Prediction of protein antigenic determinants from amino acid sequences. Proc Natl Acad Sci. 1981, 78: 3824-3828. 10.1073/pnas.78.6.3824.
Cherrier MV, Kaufmann B, Nybakken GE, Lok SM, Warren JT, Chen BR, Nelson CA, Kostyuchenko VA, Holdaway HA, Chipman PR, Kuhn RJ, Diamond MS, Rossmann MG, Fremont DH: Structural basis for the preferential recognition of immature flaviviruses by a fusion-loop antibody. EMBO J. 2009, 28: 3269-3276. 10.1038/emboj.2009.245.
Takada K, Masaki H, Konishi E, Takahashi M, Kurane I: Definition of an epitope on Japanese encephalitis virus (JEV) envelope protein recognized by JEV-specific murine CD8+ cytotoxic T lymphocytes. Arch Virol. 2000, 145: 523-534. 10.1007/s007050050043.
Ni H, Barrett AD: Nucleotide and deduced amino acid sequence of the structural protein genes of Japanese encephalitis viruses from different geographical locations. J Gen Virol. 1995, 76: 401-407. 10.1099/0022-1317-76-2-401.
Lin CW, Wu SC: A functional epitope determinant on domain III of the Japanese encephalitis virus envelope protein interacted with neutralizing-antibody combining sites. J Virol. 2003, 77: 2600-2606. 10.1128/JVI.77.4.2600-2606.2003.
Schuh AJ, Li L, Tesh RB, Innis BL, Barrett AD: Genetic characterization of early isolates of Japanese encephalitis virus: genotype II has been circulating since at least 1951. J Gen Virol. 2010, 91: 95-102. 10.1099/vir.0.013631-0.
Zhang S, Bovshik EI, Maillard R, Gromowski GD, Volk DE, Schein CH, Huang CY, Gorenstein DG, Lee JC, Barrett AD, Beasley DW: Role of BC loop residues in structure, function and antigenicity of the West Nile virus envelope protein receptor-binding domain III. Virology. 2010, 403: 85-91. 10.1016/j.virol.2010.03.038.
Sidney J, Grey HM, Southwood S, Celis E, Wentworth PA, del Guercio MF, Kubo RT, Chesnut RW, Sette A: Definition of an HLA-A3-like supermotif demonstrates the overlapping peptide-binding repertoires of common HLA molecules. Hum Immunol. 1996, 45: 79-93. 10.1016/0198-8859(95)00173-5.
Cocquerel L, Wychowski C, Minner F, Penin F, Dubuisson J: Charged residues in the transmembrane domains of hepatitis C virus glycoproteins play a major role in the processing, subcellular localization, and assembly of these envelope proteins. J Virol. 2000, 74: 3623-3633. 10.1128/JVI.74.8.3623-3633.2000.
Tang WF, Ogawa M, Eshita Y, Aono H, Makino Y: Molecular evolution of Japanese encephalitis virus isolates from swine in Oita, Japan during 1980–2009. Infect Genet Evol. 2010, 10: 329-336. 10.1016/j.meegid.2009.12.005.
Pujhari SK, Prabhakar S, Ratho RK, Modi M, Sharma M, Mishra B: A novel mutation (S227T) in domain II of the envelope gene of Japanese encephalitis virus circulating in North India. Epidemiol Infect. 2011, 139: 849-856. 10.1017/S0950268810001937.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/13/368/prepub
This work was supported by the grant from The Department of Science and Technology (DST), Government of West Bengal, India [No. 396(Sanc.)/ST/P/S&T/9G-27/2007]. The authors are indebted to Dr. Sekhar Chakrabarti, the Officer-In-Charge, ICMR Virus Unit, for allowing them to carry out the work. We also gratefully acknowledge the help we received from NIV, Pune for providing us with the ELISA kit for the detection of IgM antibody to Japanese encephalitis virus. The enthusiastic help obtained from the physicians of the district hospitals, for providing us the samples from clinically suspected JE/AES for this study, is gratefully acknowledged.
The authors declare that they have no competing interests.
AS conceived the study, the design, and drafted the manuscript. AS carried out serology and molecular work. AS, AB, BKP participated in immunoinformatics analysis. AS, AB, BKP, SKM and SC contributed to the data analysis and data interpretation. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Sources of the JEV strains/isolates used in the phylogenetic analysis in this study. All strains/isolates are JEV, except Murray Valley encephalitis virus strain (MVE-1-51) used as an out group for phylogenetic analysis in this study. */** Isolates were 100% identical with each other. (PDF 70 KB)
Additional file 2: (a) Schematic structural representation of E protein of JEV isolates. Here, DI, DII (FL), DIII, ST and TM were symbolized for domain I (red), domain II with orange colored fusion loop (yellow), domain III (blue), stem (magenta) and transmembrane (gray) region respectively including the span/stretch of respective domain/region with above shown numbering, according to previous reports [32, 33]. (b) A comparison of structure-based multiple sequence alignment of JEV isolates-specific E protein with respect to SA14-14-2 vaccine strain. Domains/regions were colored as Figure 2 (a). Secondary structural elements (described in terms of E-extended strand/β-sheet, H-α-helix and C-coil/turn/bend) were shown below each sequence alignment. Underlined bold faces amino acids represent the JEV isolates-specific amino acid substitutions in E protein. Amino acid substitutions leading to escape from antibody neutralization or neutralizing epitopes or related to neuorovirulence/neuroinvasiveness or both found in the JEV isolates with respect to SA14-14-2 vaccine strain were represented by shaded dot, dot and bold asterisk (*) respectively. One amino acid substitution at E158 (Q→P) in the isolate IND/12/WB/JEV50 resulting in secondary structural changes marked as green shaded bold faces. (PDF 79 KB)
Additional file 3: Amino acid substitutions in E protein epitopes of JEV WB isolates associated with HLA-A alleles as predicted by EpiJen serve. The values given within the bracket indicates the 50% inhibitory concentration (IC50) of the peptide, a measure of the binding affinity. An IC50 value < 50 is considered as a good affinity. Amino acid substitutions in the predicted epitopes of E protein are marked as bold. (PDF 63 KB)
Additional file 4: Amino acid substitutions in E protein epitopes of JEV WB isolates associated with HLA-B alleles as predicted by EpiJen server. The values given within the bracket indicates the 50% inhibitory concentration (IC50) of the peptide, a measure of the binding affinity. An IC50 value < 50 is considered as a good affinity. Amino acid substitutions in the predicted epitopes of E protein are marked as bold. (PDF 199 KB)
About this article
Cite this article
Sarkar, A., Banik, A., Pathak, B.K. et al. Envelope protein gene based molecular characterization of Japanese encephalitis virus clinical isolates from West Bengal, India: a comparative approach with respect to SA14-14-2 live attenuated vaccine strain. BMC Infect Dis 13, 368 (2013) doi:10.1186/1471-2334-13-368
- Japanese encephalitis virus
- Envelope protein gene
- Molecular phylogeny
- Genotype III
- Genotype I
- Homology modeling
- T-cell epitope
- West Bengal