Genotype V Japanese Encephalitis Virus Is Emerging

1Japanese encephalitis (JE), which is caused by JE virus (JEV), is one of the most important viral encephalitis in the world –. It is prevalent mostly in Asia including eastern Asia –, southern Asia and southeast Asia,. JE has extended its geographic range to the south Pacific region, including Australia,. An estimated 3 billion persons live in countries where JE is endemic –.

2JEV has a zoonotic transmission cycle between mosquitoes (principally of the genus Culex) and vertebrate hosts such as bats, water birds and pigs,. Human beings contract JEV when bitten by infected mosquitoes. Around 35,000–50,000 JE cases are reported each year, of which 10,000–15,000 are fatal –,. Approximately 50% of JE patients present severe neurological and mental sequelae such as motor deficits, and cognitive and language impairment –.

3JEV is a member of the genus Flavivirus, family Flaviviridae. Like other flaviviruses, the JEV genome is a single-stranded positive-sense RNA of approximately 11 kb in length. It is capped at its 5′ end and has a single open reading frame (ORF) that encodes a polyprotein. The ORF is flanked by 5′ and 3′ untranslated regions (UTRs). The viral structural proteins are encoded by the 5′ one-third of the ORF and consist of the capsid (C), membrane (M; formed by proteolytic cleavage of its precursor protein PrM) and envelope (E) proteins. The remaining 3′ region encodes non-structural proteins (NS1 to NS5),.

4JEVs have been divided into five genotypes (genotype I, II, III, IV, V), based on nucleotide sequence of E gene. Genotypes I-IV have been isolated from many vectors –, bats, and patients,,,, in Asia (including eastern, southern and southeast Asia) and Australia. To date, the Muar strain, which was isolated from specimens of brain tissues of patients with viral encephalitis in Malaya in 1952, is the only example of genotype V JEV,. Since that time, no genotype V JEV has been detected. In this study, genotype V JEV was isolated from Culex tritaeniorhynchus collected in China in 2009. This suggests that genotype V JEV is re-emerging in Asian country after a 57 year hiatus.

5C6/36 (Aedes albopictus) cell line was grown in minimal essential medium (HyClone) with Hanks' salt solution supplemented with 10% fetal bovine serum (FBS, HyClone), 2 mM glutamine, 0.12% NaHCO3, and 100 U ml−1 penicillin and streptomycin. Cells were propagated and maintained at 28°C. BHK-21 cells were grown in minimal essential medium (HyClone) with Earl's balanced salt solution supplemented with 10% FBS, 2 mM glutamine, 0.12% NaHCO3, and 100 U ml−1 penicillin and streptomycin. BHK-21 cells were propagated and maintained at 37°C under a 5% CO2 atmosphere.

6An arbovirus survey was conducted in Tibet in the summer of 2009. Mosquitoes were collected in Mainling County (altitude 2900 m) and Medog County (altitude 1000 m) in the Nyingchi area of Tibet. Mosquito samples were collected using mosquito-trapping lamps (Wuhan Lucky Star Environmental Protection Tech Co., Ltd., Hubei, China) in the evening. Collection locations were proximal to sites of frequent human activity. Collection nets containing mosquitoes were frozen for 30 min at −20°C and transferred onto an ice plate for determination of mosquito species (blood-fed and male mosquitoes were discarded). Female mosquitoes were identified to species level by morphologic characteristics and sorted into pools of 100 specimens according to species. The pools were put into collection tubes individually and stored in liquid nitrogen,.

7Mosquito pools were added to 1.5 ml minimal essential medium (HyClone), supplemented with 2 mM glutamine, 0.12% NaHCO3, and 100 U ml−1 penicillin and streptomycin, followed by grinding in a pre-cooled sterile plastic grinding tube using a TissueLyser (GIAGEN, Germany). Homogenized samples were centrifuged at 17,000× g in a microcentrifuge for 20 min at 4°C, and the clarified supernatants were used to inoculate monolayers of BHK-21 and C6/36 cells and incubated at 37°C and 28°C, respectively. The cells were observed daily to check for development of cytopathic effects (CPE). A sample was regarded as virus-positive if it caused CPE in successive cell passages,.

8Viral supernatants were applied to six-well plates (Corning, USA) of confluent BHK-21 cells and incubated for one hour. Plates were first overlaid with medium containing 75% agarose and then with medium containing neutral red vital stain after three days incubation at 37°C in a 5% CO2 incubator. Plaques of different sizes and shape were shattered in 500 ul MEM medium after being picked out using a sterile pipette tip. As described previously, this process was repeated until a single plaque-shaped virus was obtained.

9Viral RNA was extracted from 140 ul supernatant from virus-infected BHK-21 cell cultures using a Viral RNA Mini Kit (QIAGEN, Germany) according to the manufacturer's instructions. cDNA was synthesized using a Ready-to-Go You-Prime First-Strand Beads Kit (GE healthcare, UK) and random hexanucleotide primers. PCR amplification using universal primers specific for flaviviruses, alphaviruses and bunyaviruses was conducted for identification of virus isolates. Primers (Table 1) were designed for full-length genome amplification and sequencing of JEV using the PREMIER Primer 5 software package. Thermal cycling parameters were as follows: one cycle of denaturation (94°C, 5 min) followed by 35 cycles of 94°C denaturation for 30 s, 55°C annealing for 30 s, and 72°C extension for 1 min. The programme ended with an extension step at 72°C for 10 min. Amplified products were examined by agarose gel electrophoresis (1%), purified using a QIAquick Gel Extraction kit (QIAGEN, Germany), and then sequenced directly. Sequencing of the 5′ UTR and 3′ UTR were determined by using 5′ RACE and 3′ RACE system for Rapid Amplification of cDNA Ends (Invitrogen) respectively. 5′ RACE was performed according to standard protocols (Invitrogen 5′ RACE kit). 3′ RACE was performed by first adding a polyA tail using PolyA polymerase (New England Biolabs) and then conduct RT-PCR with gene specific primers and an oligo-dT-adapter primer.

10The full-length genome of the XZ0934 strain was compiled using SeqMan in the Lasergene software package (DNASTAR). Nucleotide and amino acid sequence alignments were generated by ClustalX version 2.0.9,. Analysis of nucleotide and deduced amino acid sequence identities was performed using GeneDoc and MegAlign in the Lasergene software package (DNASTAR).

11Full-length nucleotide sequences of 32 selected JEV strains of varying genotype isolated from different locations and sources, and across a number of years, were downloaded from GenBank (Table 2). The reported structural gene sequence of genotype V JEV (Muar strain) was used to build phylogenetic trees. Neighbor-joining phylogenetic trees based on nucleotide sequences were constructed using MEGA version 4.0.2,. The robustness of phylogenetic constructions was evaluated by bootstrapping using 1000 replicates.

12To better understand the phylogenetic relationship between genotype V JEV and other flaviviruses, full-length nucleotide sequences of previously published JEV strains and other flaviviruses were downloaded from GenBank, including sequences from Murray Valley encephalitis virus (MVEV), West Nile virus (WNV), Kunjin virus (KUNV), St. Louis encephalitis virus (SLEV), Dengue virus 1 (DENV1, Dengue virus 2 (DENV2), Dengue virus 3 (DENV3), Dengue virus 4 (DENV4), Yellow fever virus (YFV), Powassan virus (POWV), Langat virus (LANV), Louping ill virus (LIV), Tick-borne encephalitis virus (TBEV) and Culex flavivirus (Table 3).

13After homogenized supernatants were inoculated onto monolayers of BHK-21 and C6/36 cells, a single pool containing 100 specimens of Culex tritaeniorhynchus yielded a virus isolate designated XZ0934. The supernatant of pool XZ0934 caused cytopathic effects (CPE) in BHK-21 and C6/36 cells in successive cell passages. The C6/36 cells became aggregate, and showed fusion and shedding while the BHK-21 cells became aggregate and began shedding by 72 h post-infection. All plaques in BHK-21 cell monolayers were of identical size (mean 1.5 mm, n = 10). Two plaques were picked from them and subjected to a second round of plaque purification. The resultant data were consistent with the former.

14Viral RNA was extracted and amplified by PCR using primers specific for flaviviruses, alphaviruses and bunyaviruses. XZ0934 was positive when primers specific for flaviviruses (FU1/cFD2) were used, and nucleotide sequencing confirmed that XZ0934 was a JEV. To ensure the consistency of different viral plaques, six purified plaques were picked and amplified using flavivirus-specific primers (FU1/cFD2). The nucleotide and amino acid sequence identities of the six purified plaques were 100%. This indicates that each of the six purified plaques was generated by an identical JEV strain. Therefore, one plaque was selected for full-length genome amplification and sequencing.

15Recent reports have suggested that JEVs currently circulating in China belong to genotypes I and III –, –. Thus, 32 primers were designed using the complete sequences of genotype I JEV Ishikawa (GenBank accession number AB051292) and 48 from the sequence of genotype III JEV Beijing-1 (GenBank accession number L48961). These were used for amplification of the entire XZ0934 genome. PCRs were positive with 4 genotype I and 10 genotype III primers. Based on obtained nucleotide sequences, primers were designed to close the majority of gaps between assembled contigs by PCR amplification in order to determine the whole genome of XZ0934. A further 24 primers (Table 1) were designed and used to verify the accuracy of sequencing. The complete genome (10,983 nt) of XZ0934 was sequenced (GenBank accession no. JF915894) and found to possess one open reading frame (ORF). When the complete genome sequence of isolate XZ0934 was compared with those of 62 known JEV isolates (genotypes I–IV) in Genbank, sequence identities varied from 78.6% (KV1899, K94P05) to 79.7% (CC27-L1) and amino acid sequence identity from 90.0% (KV1899, K94P05) to 91.6% (K87P39). Thus, these data reveal low similarity between XZ0934 and genotype I–IV JEVs.

16Because the structural gene sequence of genotype V (Muar) has been reported, an identity analysis of JEV structural genes (C, PrM, M, E) of XZ0934, Muar and other selected genotype I–IV JEV strains was conducted (Table 4). C gene sequence homology varied from 78.2% (G IV, JKT 6468) to 88.5% (G V, Muar) for nucleotides and 72.4% (G IV, JKT 6468) to 85.8% (G V, Muar) for amino acids. That of the PrM gene varied from 71.7% (G IV, JKT 6468) to 84.1% (G V, Muar) for nucleotides and 81.5% (G IV, JKT 6468) to 90.2% (G V, Muar) for amino acids. M gene sequence homology varied from 80.0% (G IV, JKT 6468) to 95.6% (G V, Muar) for nucleotides and 85.3% (G IV, JKT 6468) to 100.0% (G V, Muar) for amino acids. E gene sequence homology varied from 77.0% (G I, Ishikawa) to 86.0% (G V, Muar) for nucleotides and 89.4% (G I, Ishikawa) to 93.2% (G V, Muar) for amino acids. These data demonstrate that the structural gene sequence homology of XZ0934 was higher with genotype V JEV (Muar) than with other genotype I–IV JEV strains.

17To establish the phylogenetic relationship between XZ0934 and other JEV strains, a phylogenetic tree was constructed using the complete genome sequences of XZ0934 and 32 selected JEV strains (genotypes I–IV). Murray Valley encephalitis virus (MVEV) was used as an outgroup. Five distinct phylogenetic groups were identified. The XZ0934 strain, which was isolated from China, formed a branch divergent from other genotype I–IV JEV strains (Figure 1A). Therefore, XZ0934 should be regarded as a novel, non-genotype I–IV, JEV isolate.

18To study their phylogenetic relationship, a phylogenetic tree was constructed using the reported structural gene nucleotide sequences of Muar, XZ0934, and other JEV strains (genotype I–IV). No matter which structural gene was used to construct the phylogenetic tree, the topology was similar. Five distinct phylogenetic groups were evident in each tree. XZ0934 and Muar fell into the same group when the tree was constructed using the C (Figure 1B), PrM (Figure 1C), M (Figure 1D) or E (Figure 1E) genes. This result suggested that XZ0934 was a novel genotype V JEV isolate.

19A phylogenetic tree was also constructed using genomic nucleotide sequences in order to understand the phylogenetic relationship between XZ0934 and other flaviviruses. Data indicated that XZ0934 was indeed a JEV rather than any of the other 14 flaviviruses (Figure 2).

20In recent years, the sequence of the JEV viral envelope (E) gene has been used by various authors to perform phylogenetic analyses,,,. Based on the resultant data, JEV strains have been divided into five genotypes (genotypes I-V). Genotypes I and III are distributed widely in Asia, including Japan, Korea, China, India, Vietnam and Philippines. Genotype II includes isolates from southern Thailand, Malaysia, Indonesia, and northern Australia. Genotype IV has been isolated only in Indonesia. The Muar strain, isolated in Malaya in 1952, is regarded as the only genotype V JEV isolate,,. In this study, phylogenetic analysis of structural genes and whole genome sequences also suggested the existence of five JEV genotypes. Thus, the E gene is confirmed to be a useful phylogenetic marker for JEV.

21Primers designed for JEV genotypes I and III were used for full-length amplification of XZ0934. Of these, only a few (4/32 genotype I and 10/48 genotype III) resulted in successful amplification. This suggests a low whole genome sequence homology between XZ0934 and genotype I and I JEV isolates. In order to further understand the differences between XZ0934 and other JEV strains (genotype I–IV), an identity analysis was conducted using the full-length nucleotide sequences of XZ0934 and 62 known JEV isolates (genotypes I–IV) in Genbank. Data suggested that XZ0934 and the genotype I–IV JEV strains were dissimilar. The nucleotide sequence identity varied from 78.6% to 79.7% and amino acid sequence identity from 90.0% to 91.6%. Indeed, the sequence divergence ranged from 20.3% to 21.4% (nt) and 8.4%–10.0% (aa). It has been suggested that the nucleotide sequence divergence between different JEV genotypes is ∼10%. The sequence divergence (20.3%–21.4%) between XZ0934 and the genotype I–IV JEVs was greater than 10%, suggesting that XZ0934 is not a member of JEV genotypes I–IV. To confirm that XZ0934 was a JEV and not some other flavivirus, 14 flavivirus strains, including mosquito-borne and tick-borne flaviviruses, were used to build a phylogenetic tree. The data indicated that XZ0934 was indeed a JEV (Fig. 2).

22Four viral encephalitis cases were reported in Malaya (n = 1) and Singapore (n = 3) in the summer of 1952. All patients exhibited high fever, vomiting, headache, disturbance of consciousness, stiff neck and deep coma with rapid progression to death by respiratory failure. Four virus strains were isolated from brain tissue specimens and identified as JEV by neutralization test using the Japanese Nakayama JEV strain. Of these, the Muar strain, isolated from a 19-year-old male patient in Malaya in 1952, has been assigned to genotype V based on the E gene sequence,,,. During the following 57 years (1952–2009), no genotype V JEV has been reported. In this study, XZ0934, isolated from Culex tritaeniorhynchus collected in China, has been identified as a genotype V JEV, based on phylogenetic analysis using both full-length genome and structural gene nucleotide sequences. This represents only the second instance of isolation of genotype V JEV worldwide since 1952. Thus, genotype V JEV is not limited to southeast Asia and has begun to be emerge in the world.

23Many factors may contribute to spread of JEV,, such as changed agricultural practices (which provide new breeding sites for mosquitoes), animal husbandry (which provides host animals for transmission), migrating birds and even wind-blown mosquitoes,. Each of the five known JEV genotypes originated in the Indonesia-Malaysia region, so why has genotype V JEV not been detected for 57 years? How did it spread to China from southeast Asia, a distance of thousands of kilometers? Does this virus exist somewhere along the path from Malaysia to China? All these issues are worthy of further study. Moreover, genotype V JEV was first isolated from human specimens, suggesting a high pathogenicity and the possibility of viral encephalitis. Therefore, increased surveillance and more effective diagnosis of viral encephalitis caused by genotype V JEV is an issue of great concern to nations in which JEV is endemic.