Next Article in Journal
Comparative Study on Immune Function of the Head and Trunk Kidney in Rainbow Trout Responding to IHNV Infection
Next Article in Special Issue
Growth, Pathogenesis, and Serological Characteristics of the Japanese Encephalitis Virus Genotype IV Recent Strain 19CxBa-83-Cv
Previous Article in Journal
Rabies Virus Populations in Humans and Mice Show Minor Inter-Host Variability within Various Central Nervous System Regions and Peripheral Tissues
Previous Article in Special Issue
Japanese Encephalitis Virus: The Emergence of Genotype IV in Australia and Its Potential Endemicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore

1
Environmental Health Institute, National Environment Agency, Singapore 228231, Singapore
2
School of Biological Sciences, Nanyang Technological University, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(12), 2662; https://doi.org/10.3390/v14122662
Submission received: 18 October 2022 / Revised: 23 November 2022 / Accepted: 24 November 2022 / Published: 28 November 2022
(This article belongs to the Special Issue Japanese Encephalitis Virus)

Abstract

:
Japanese encephalitis virus (JEV) is an important arbovirus in Asia that can cause serious neurological disease. JEV is transmitted by mosquitoes in an enzootic cycle involving porcine and avian reservoirs, in which humans are accidental, dead-end hosts. JEV is currently not endemic in Singapore, after pig farming was abolished in 1992; the last known human case was reported in 2005. However, due to its location along the East-Asian Australasian Flyway (EAAF), Singapore is vulnerable to JEV re-introduction from the endemic regions. Serological and genetic evidence in the last decade suggests JEV’s presence in the local fauna. In the present study, we report the genetic characterization and the first isolation of JEV from 3214 mosquito pools consisting of 41,843 Culex mosquitoes, which were trapped from April 2014 to May 2021. The findings demonstrated the presence of genotype I of JEV (n = 10), in contrast to the previous reports of the presence of genotype II of JEV in Singapore. The genetic analyses also suggested that JEV has entered Singapore on several occasions and has potentially established an enzootic cycle in the local fauna. These observations have important implications in the risk assessment and the control of Japanese encephalitis in non-endemic countries, such as Singapore, that are at risk for JEV transmission.

1. Introduction

Japanese encephalitis (JE) is the most significant vaccine-preventable neurological disease in Asia and the Western Pacific region [1,2]. The first confirmed cases of JE were detected in Japan in 1924, and it was subsequently reported in the rest of Asia, as well as in parts of Russia and Australia. At least 24 countries are currently considered to be at risk for JE transmission [1,3,4]. Based on 2011 estimates, the global annual caseload of JE was approximately 68,000 cases, with at least 13,000 deaths [5]. The global burden of JE is likely to increase as regions that were formerly less favorable to JE experience more hospitable conditions for its spread due to climate change [6,7,8]. This is most evident in the recent JE outbreak in Australia, with 37 locally-acquired cases reported; the outbreak occurred 20 years after the last known human case [9].
Japanese encephalitis virus (JEV), the causative agent of JE, is a positive-sense, single-stranded RNA virus of the Flaviviridae genus. JEV has five known genotypes (GI–GV) [10,11,12,13,14], including two sub-genotypes of GI—GI-a and GI-b [15]. The JEV genome is approximately 11 kilobases long and encodes ten proteins in an open reading frame flanked by two non-coding regions (NCRs) at the 5′ and 3′ ends. The ten proteins that are encoded are three structural proteins (capsid, membrane, and envelope) and seven non-structural (NS) proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) [16,17].
JEV is transmitted by Culex mosquitoes; Culex tritaeniorhynchus is the primary vector [18,19], and other Culex mosquitoes are secondary vectors. Other genera of mosquitoes, including Aedes, Anopheles, and Mansonia, may act as regional or potential secondary vectors (reviewed by Pearce, et al. [20]). The main reservoir hosts of JEV are birds, primarily of the family Ardeidae, while swine act as amplifying hosts [1,21]. Humans are incidental, dead-end hosts [2,22]. Before 2000, GIII was the most prominent JEV genotype in all endemic areas, but it is currently being supplanted by GI in Asia [23,24,25].
Prior to the abolishment of pig farming in 1992, Singapore had an average of 14 JE cases annually, but has had almost no cases since then [26,27]. However, the country has been identified as one of the 24 countries with a JE transmission risk [3,4]. A highly-urbanized island city-state, Singapore nevertheless retains much of its greenery and biodiversity through Nature Reserves (NRs), forest parks and urban parks, and mangroves [28]. Its tropical climate, combined with its varied forest and urban habitats, provide a conducive environment for at least 180 species of mosquitoes, including Cx. tritaeniorhynchus [29]. The country is also home to wild boars [30,31,32] and local wild birds. Additionally, migratory birds stopover in Singapore, which is located along the East-Asian Australasian Flyway (EAAF) [33]. Ardeid birds are present in both the local (e.g., Grey Herons and Yellow Bitterns) and the migratory bird populations (e.g., Little Egrets, Chinese Pond Herons, and Black Bitterns) [34].
A serological study that was carried out from 2010 to 2018, as well as an investigation of serum samples that were collected between 1996 and 1999, indicated JEV exposure in the wild boar and the wild bird populations in Singapore [26,35], suggesting that JEV may be enzootic in Singapore. Entomology-based surveillance that was initiated in 2011 detected five JEV-positive mosquito pools at a local Avian Sanctuary, further supporting this notion. The JEV sequences that were derived from those mosquito pools clustered with JEV GII isolates detected in Indonesia [36]. The seroconversion of sentinel chickens during surveillance carried out from 2013 to 2014 further suggested active JEV transmission in the rural areas of Singapore [35]. Although JEV had previously been detected in mosquito vectors in Singapore [36,37], virus isolation had not been successful [36]. In the present study, we report the first isolation and detection of JEV GI from Culex mosquitoes in Singapore. This study also describes the genetic characteristics of the newly detected JEV strains. The findings suggest that JEV has potentially established enzootic transmission locally, which has important implications on the risk assessment and vector control policies in Singapore.

2. Materials and Methods

2.1. Entomological and Virological Surveillance

Entomological surveillance and virological surveillance of JEV was carried out near-monthly; this occurred in two phases—a single-site phase from April 2014 to August 2016, and a multiple-site phase from September 2019 to May 2021. During the single-site phase, one surveillance site was set up at Park 1, which was in the north-eastern region of Singapore. During the multiple-site phase, four surveillance sites were set up at NR 1, NR 2-Park 2, Park 3, and Park 4 (Figure 1). NR 2 and Park 2 were considered as one site due to their proximity to each other. These study sites were selected based on their proximity to slow-moving or still water bodies and the available information on the presence of JEV hosts, such as ardeid birds and wild boars, and the presence of JEV vectors. These selection criteria were adopted to maximize the possibility of capturing JEV-positive mosquitoes because of the lack of evidence on active transmission of JEV in Singapore.
Mosquitoes were trapped over two consecutive nights during the single-site phase, and for 18–20 h during the multiple-site phase. Each trapping session was conducted overnight, including crepuscular periods, using Centers for Disease Control (CDC) light traps baited with dry ice and incandescent light. Four to nine locations were selected at each surveillance site, and traps were positioned approximately 1.5–2.0 m above ground in shaded areas adjacent to water bodies where possible. Traps were set between 15:00 and 17:00 h and were collected the next day between 09:00 and 11:00 h.
Additionally, Culex genus mosquitoes that were trapped from November 2020 to May 2021 as part of the malaria surveillance activities of the National Environment Agency (NEA) were also included in this study. The Malaria Surveillance Section (MSS) conducts on-going catches using the Night Catcher (a modified form of the CDC light trap [38]) in 45-day cycles at 116 different locations in Singapore. Each MSS trapping session was carried out over a 12-h period overnight, including crepuscular periods. The two MSS trapping locations where JEV-positive pools were observed are indicated as MSS 1 and MSS 2 (Figure 1).

2.2. Mosquito Identification

Mosquito specimens were transported on dry ice to the Environmental Health Institute (EHI), where they were sorted using mosquito identification keys [39,40,41,42,43,44,45]. For the single-site phase, mosquito specimens were identified to species level where possible. For the multiple-site phase, mosquito specimens were initially identified to genus level, until the first detection of JEV, whereupon they were identified to species level where possible. In the case of morphologically similar species, the specimens were identified to the species group. Damaged specimens were identified to the genus level. For the start of the multiple-site phase, all female Culex mosquitoes were pooled by trapping location and trapping month in pools of up to 100 specimens. Otherwise, females from six Culex vector species—Cx. bitaeniorhynchus, Cx. gelidus, Cx. quinquefasciatus, Cx. sitiens, Cx. tritaeniorhynchus, and Culex spp. vishnui subgroup (excluding Cx. tritaeniorhynchus)—were pooled according to the species, trapping location, and trapping month, in pools of up to 100 specimens. All specimen pools were stored at −80 °C until further analysis.

2.3. RNA Extraction from Mosquito Specimens

Pooled mosquito samples were homogenized in 500 μL of universal transport media (Copan Diagnostics, Murrieta, CA, USA) using either the Mixer Mill MM 400 (Retsch Technology GmbH, Haan, Germany) or the 1600 MiniG homogeniser (Spex SamplePrep, Metuchen, NJ, USA), and centrifuged at 5000× g for 5 min. Viral RNA was extracted from the mosquitoes collected during the single-site phase using 140 μL of the homogenate and the QIAamp Viral RNA mini kit (Qiagen Group, Hilden, Germany), according to the manufacturer’s instructions. For mosquitoes collected during the multi-site phase, as well as the mosquitoes obtained from MSS, viral RNA was extracted from 200 µL of mosquito homogenate using the IndiMag Pathogen kit (Indical Bioscience GmbH, Leipzig, Germany) and the KingFisher Flex system (Thermo Fisher Scientific Inc., Waltham, MA, USA), as per manufacturer’s instructions. Eluted RNA and the remaining homogenates were stored at −80 °C until further use.

2.4. Detection of JEV by Reverse Transcriptase-PCR (RT-PCR)

Extracted RNA was screened for JEV, as described in Yap, Mailepessov, Lim, Chan, How, Humaidi, Yeo, Chong, Lam-Phua, Lee, Okumura, Vythilingam and Ng [36]. Briefly, positive RNA controls of JEV (Nakayama strain) were extracted from Vero cell (ATCC CCL-81; ATCC, Manassas, VA, USA) supernatants using a QIAmp Viral RNA Mini kit (Qiagen), as described above. JEV was detected by using a modified RT-PCR assay adopted from Scherret, et al. [46] and Santhosh, et al. [47]. All reactions were performed using the QuantiTect® SYBR® Green RT-PCR kit (Qiagen) in a final volume of 20 µL, with 5 µL of extracted RNA and 10 µM of each primer. The cycling conditions were as follows: reverse transcription at 50 °C for 20 min and inactivation at 95 °C for 15 min; 40 cycles at 94 °C for 15 s, 58 °C for 30 s, and 72 °C for 30 s; melting curve analysis was conducted from 65 to 95 °C, with a slope of 0.1 °C/s, for validation of the amplified products; cycling conditions concluded with a cooling step at 37 °C for 20 s.

2.5. Sequencing and Analysis of the Envelope Gene of JEV

The envelope (E) gene of JEV was amplified using two pairs of primers, as described in Schuh, et al. [48]. Complementary DNA was first synthesized using the Maxima H Minus First Strand cDNA synthesis kit (Thermo Fisher Scientific Inc.), and subsequent amplifications were performed using Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific Inc.), as per the manufacturer’s recommendations. The amplified products (1700 and 800 bp) were purified using the FavorPrep GEL/PCR Purification Kit (Favorgen, Ping Tung, Taiwan) before sequencing at a commercial sequencing facility using the BigDye terminator Cycle Sequencing Kit protocol (Applied Biosystems, Thermo Fisher Scientific Inc., Waltham, MA, USA).
Raw sequence data were analyzed as described in Yap, Mailepessov, Lim, Chan, How, Humaidi, Yeo, Chong, Lam-Phua, Lee, Okumura, Vythilingam and Ng [36]. A phylogenetic tree was constructed using 129 truncated (701 bp) JEV E-gene sequences; this included the newly detected JEV sequences in the present study (n = 10) and 119 JEV sequences that were retrieved from the NCBI GenBank database (Table S1). The sequences were aligned using MUSCLE in MEGA version X 10.2.5 [49], and a phylogenetic tree was constructed using the maximum likelihood (ML) method. The best fitting model, which was the Tamura-Nei model [50] with gamma-distributed substitution sites (TN93+G), was determined via the model selection tool in MEGA version X 10.2.5. The tree was constructed using TN93+G with 1000 bootstrap reiterations and rooted with the truncated (701 bp) E-gene from Murray Valley Encephalitis virus (MVEV; NCBI Reference Sequence no. NC_000943.1). Amino acid sequences were also obtained via MEGA version X 10.2.5. Sequence identity matrices were generated with the MUSCLE tool in EMBL-EBI with default parameters [51,52]. The newly generated E-gene sequences from this study were deposited in the GenBank nucleotide database (GenBank Accession Nos. KT346354, MW802631, MW802633 to MW802637, and ON804794 to ON804796).

2.6. JEV Isolation

Mosquito homogenates that were positive for JEV were filtered (Minisart 0.45 microns; Sartorius, Göttingen, Germany) before inoculation onto C6/36 cells (ATCC CRL-1660; ATCC) at 90% confluency. The C6/36 cells were incubated in Leibovitz L-15 media (Thermo Fisher Scientific Inc.) with 2% fetal bovine serum (FBS) (Thermo Fisher Scientific Inc.) at 33 °C and monitored for one week. Approximately 1/5 of cell culture suspension was then transferred to freshly prepared 90% confluent C6/36 cells, topped up with fresh 2% FBS Leibovitz L-15 media, and incubated again at 33 °C for one week. This step was repeated one more time for a maximum of three passages. If cytopathic effects or syncytia were observed, the cells were harvested and kept at −80 °C.

2.7. Sequencing and Phylogenetic Analysis of the Whole Genomes of JEV

Viral RNA was extracted from cell culture suspension using the QIAmp Viral RNA Mini kit (Qiagen), according to manufacturer’s instructions. Multiple amplicons covering the whole genome were then synthesized using the primers in Table S3 and a protocol similar to that of E-gene amplification, except with an annealing temperature of 60 °C. The amplification products were purified with the FavorPrep GEL/PCR purification kit (Favorgen) and sequenced at a commercial sequencing facility using the BigDye terminator Cycle Sequencing Kit protocol (Applied Biosystems, Thermo Fisher Scientific Inc.).
Raw sequence data were trimmed and assembled using SeqMan Pro 15.3.0 from the Lasergene software suite (DNAStar Inc., Madison, WI, USA) and BioEdit 7.2.5 to obtain consensus sequences. Where a consensus could not be reached for a sequence call in the contig, the position was represented by the appropriate IUPAC ambiguous base symbol.
A whole-genome phylogenetic tree was constructed using 50 JEV sequences; 3 sequences were generated during the present study and 47 sequences were obtained from the NCBI GenBank database (Table S4). The sequences were aligned in MEGA version X 10.2.5; as above, the best fit model was determined to be the Tamura-Nei model, but with invariant and gamma-distributed substitution sites (TN93+G+I). The tree was constructed using TN93+G+I with 1000 bootstrap iterations and rooted with MVEV (NCBI Reference Sequence no. NC_000943.1). The amino acid translation was carried out via the Translate tools from Expasy [53] and the Sequence Manipulation Suite [54]. Sequence identity matrices of whole genome data sets were also generated as above. The whole genome sequences that were generated during the present study were deposited in the GenBank nucleotide database (GenBank Accession Nos. ON804797 to ON804799).

2.8. Mutation Analyses

Amino acid mutations that were unique to the study sequences were determined by comparing them with all of the JEV GI complete coding sequences available in the NCBI database (n = 158) as of 03 November 2022. Genotype I sequences were selected by constructing a phylogenetic tree consisting of all complete coding sequences of JEV (n = 404) available in the NCBI database. A consensus sequence of the GI complete coding sequences (n = 158) was generated using the PAM40 matrix with the cons module from the EMBOSS suite [55].

3. Results

3.1. Detection of JEV in Field-Caught Mosquitoes

A total of 54,308 Culex mosquitoes collected with CDC light traps and Night Catcher traps were screened in 3261 pools (Table 1). JEV was detected by RT-PCR in five Culex spp. pools and five Cx. tritaeniorhynchus pools, and virus isolation was successful in three pools. JEV was not found in the other five Culex species/species group that were screened.
Of the seven JEV-positive pools that were collected with the CDC light traps, one was collected from the north-eastern Park 1 in 2014, and six were collected from the north-western NR 1 in 2019. Of the three positive pools that were collected with the Night Catcher traps, two were collected in the north-eastern MSS 1 in 2020, while one was collected in the northern MSS 2 in 2021 (Table 2).

3.2. Sequencing, Phylogenetic Analysis, and Genetic Characterization of Detected JEVs

Partial E-gene sequences of JEV (701 bp) were obtained from all ten of the positive pools. Direct PCR amplification of complete genome fragments using the RNA extracted from the JEV-positive mosquito pools was not successful. However, complete genome sequences were generated from three pools where JEV was successfully isolated. All newly-generated JEV E-gene sequences (n = 10) belonged to GI (Figure 2 and Figure S2). They shared 91.58–100% nucleotide and 95.71–100% amino acid similarity among them (Table 3). Nine E-gene sequences belonged to two distinct clades (clusters A and B) in sub-genotype GI-a, while the remaining sequence (SG/EHI-CT1372) belonged to GI-b (Figure 2 and Figure S1). SG/EHI-CT1372 clustered with ZJ-09-52 (GenBank Accession No. JN216865), which was reported from China in 2009. SG/EHI-CX135, -193, SG/EHI-CX 202, -209, and -214 formed a distinct cluster (cluster A) that was closely related to the isolate JE_CP_49 (GenBank Accession No. DQ087974) from Thailand in 2003. These five local sequences in cluster A demonstrated high (99.43–100%) nucleotide sequence identity and 100% amino acid sequence identity. SG/EHI-CT1710, -MS_CT310, -MS_CT261, and -M_CT21_45 formed another cluster (cluster B) that was most closely related to a cluster of sequences (639A37Cx-tri, C081, and JEV/sw/Thailand/185/2017; GenBank Accession Nos. KY927815-816, and LC461958) from Cambodia and Thailand in the 2010s. Cluster B sequences showed high nucleotide and amino acid similarity (99–100% and 100%, respectively) between one another, except for SG/EHI-M_CT21_45, which showed an amino acid similarity of 99.57%.
Similar to the E-gene analysis, all three of the whole genomes that have been reported in the present study (SG/EHI-CT1372, -CX135, and -MS_CT261) belonged to GI. They shared 92.99–93.33% nucleotide and 94.39–96.75% amino acid similarity among them, demonstrating a clear genetic distance between one another (Table 4). SG/EHI-CT1372 was the most distant among the three isolates and belonged to GI-b (Figure 3). This isolate was most closely related to an isolate that was reported from China in 1983 (GenBank Accession No. JF706282). The remaining two sequences belonged to GI-a and clustered with JEV isolates that were reported from China (1977 & 1982), Thailand (1982, 1985, & 2014), Australia (2000–2004), and Cambodia (2014 & 2015) (GenBank Accession Nos. KT957422-423, GQ902058-061, KY927815-816, LC461958, and MT253732-734).
Among the 5′-NCRs of the three newly-detected whole genome sequences, there was only one mutation in SG/EHI-CT1372 (T > G) that was located in one of the loops of the first 5′-NCR stem loop structure (SLA). On the other hand, the 3′-NCRs had 28 mutations among the three isolates; 23 of these were in RNA secondary structures (Table S5). A total of 13 of the 23 mutations were positioned in stems of various 3′-NCR structures, while the remaining 10 mutations were in loop structures. A total of 5 of the 23 mutations belonged to significant regions in these structures, as described by Liu, et al. [56] (Table 5).

3.3. Amino Acid Differences between the New Singaporean JEV Sequences and the GI Consensus Sequence Derived from Published GI Sequences

A total of 30 unique amino acid differences were observed in the study sequences (n = 10) as compared to the JEV GI consensus sequence (Table 6 and Table 7). These mutations were novel among all of the JEV GI sequences (n = 158) that were reported in the NCBI database as of November 2022. Of these, 21 mutations were found in the non-structural proteins (Table 7). Of the remaining 9 mutations in the structural polyprotein, 7 substitutions were found in the E protein (Table 6). The mutation profiles of three complete genomes (SG/EHI-CT1372, SG/EHI-CX135, and SG/EHI-MS_CT261) were clearly distinguishable, indicating a clear genetic distinction among them. Notably, the unique mutation profile of SG/EHI-MS_CT261 consisted of more mutations in the non-structural polyprotein than the remaining two isolates. Among the 30 novel mutations, non-conservative substitutions occurred at 12 positions. Their distribution was proportionately higher in the structural polyprotein, suggesting that there is a greater impact due to novel genetic changes on the makeup of the structural proteins, especially the E protein, than the non-structural proteins of these three isolates.
Several other non-unique, but uncommon, non-conservative mutations were also noteworthy. They included Capsid-T90I (polyprotein residue 90; SG/EHI-CT1372), Membrane-D40H (polyprotein residue 167; SG/EHI-MS_CT261), Envelope-I126T (residue 420 in the polyprotein; SG/EHI-MS_CT261), and NS1-Q51L (residue 845 on the polyprotein, SG/EHI-CX135 and SG/EHI-MS_CT261). Capsid-T90I was only detected in a few sequences that were reported previously from Australia (GenBank Accession Nos. AF217620, MT253735-737). Similarly, Membrane-D40H (M-D40H) was detected in a few sequences from China (GenBank Accession Nos. KT229574, KU351668, and MN544780). Envelope-I126T (E-I126T) showed a wide geographical presence in South-east Asia (Malaysia—HQ223286, Cambodia—KY927815-816, and Thailand—LC461958), while NS1-Q51L was detected in three sequences previously reported from Cambodia (GenBank Accession Nos. KY927815-816) and Thailand (GenBank Accession No. LC461958). Besides their potential phenotypic implications, these mutations also provided clues on their likely origin. In agreement with the phylogenetic analyses, these uncommon mutations suggested that SG/EHI-CT1372 was distinct and more likely to have originated from a different geographical area than SG/EHI-MS_CT261 and SG/EHI-CX135.

4. Discussion

All of the JEV sequences (n = 10) generated during the present study from the mosquito pools collected from 2014 to 2021 belonged to GI. This is in contrast to the detection of JEV belonging to GII locally in 2011 [36]. This observation provides further empirical evidence for the establishment of genotype I as the predominant JEV genotype in Asia [23,24,25,59]. Nine of the ten new JEV E-gene sequences belonged to sub-genotype GI-a, which is strongly associated with the tropics, and is often found in South-east Asia [15]. The remaining sequence (SG/EHI-CT1372) belonged to sub-genotype GI-b, which is also known as a temperate sub-genotype, and is typically found in East Asia [15].
The nine GI-a sequences formed two clusters (cluster A and cluster B, Figure 2). Both of these clusters were closely related to JEV sequences from South-east Asian countries located to the north of Singapore. This suggested that local JEV strains were likely to be introduced by migrating ardeid birds via the EAAF, which stretches from Russia and Alaska to Australia and New Zealand through East Asia and South-east Asia.
Cluster A consisted of five sequences detected in Culex mosquitoes within the same area during the same period (Table 2). However, no subsequent instances of JEV belonging to cluster A were observed in the following years, despite continued surveillance in the same location. This suggested that cluster A could not sustain transmission in the local ecological setting and became extinct, similar to what was observed with JEV strains “B’/B/Ishikawa type” and “C’/C/Ishikawa + Kagawa type” in Toyama, Japan by Obara, et al. [65]; this extinction could either have been due to stochastic reasons or non-favorable virus characteristics for its survival in a new environment. In contrast, the detection of cluster B four times from 2014 to 2021 may suggest its persistent presence through an enzootic transmission cycle involving the fauna in the northern and north-eastern parts of Singapore (Figure 1 and Figure 2, Table 2). Although there have been no cases of JE in Singapore for almost two decades [27], previous studies, including serological evidence of JEV infection in the local and migratory birds, and the wild boars, as well as seroconversion in sentinel chickens, have shown that JEV continues to be present in the animal populations in Singapore [26,35,36] and suggest that it may currently be enzootic in Singapore. Alternatively, cluster B JEV strains could have been imported on repeated occasions by migratory birds. These assumptions appear to support the observations that were made by Kuwata, et al. [66] on the importance of vector and/or host movement in transferring JEV strains between different regions and the maintenance of JEV strains in the local fauna within the same region over a period of time.
The detection of different variants of JEV suggested that multiple incursions of JEV have occurred in Singapore, not only over the course of this study (2014–2016 and 2019–2021), but also during the period that was studied previously [36]. Besides the phylogenetic analyses, the mutation profiling also supported this notion. For example, uncommon mutation profiles suggested that SG/EHI-CT1372 was distinct and more likely to have originated from a different geographical area than SG/EHI-MS_CT261 and SG/EHI-CX135.
Although JEV has only been detected in the rural forested areas of Singapore thus far, increasing urbanization and the push for ecological conservation [67,68] result in close interactions between human dwellings and the habitats that are favored by culicine vectors and JEV hosts. A survey of bloodmeals of Culex mosquitoes in Singapore showed that JEV vectors opportunistically feed on humans [69], while a study of changing JEV epidemiological characteristics in metropolitan areas showed an apparent correlation between the presence of herons and the number of JE cases in the absence of pig farming [70]. It is thus necessary to maintain entomological and virological surveillance to keep track of possible enzootic JEV variants, thereby allowing for more leeway in containing any vector–human transmission events.
Besides gaining insights into its molecular epidemiology, the genetic characterization of JEV is also important for understanding virus fitness and virulence. Both the coding and non-coding regions of JEV are important for its replication and virulence [16,17]. While the structural and non-structural proteins provide the backbone of virion structure and replication machinery, 5′- and 3′-NCRs are also involved in genome circularization and host immune response attenuation. Genome circularization, which involves the stem loop structure SLA, is crucial for flavivirus replication, and exoribonuclease-resistant RNA (xrRNA) structures are able to mitigate the degradation of flavivirus RNA [56,71,72,73].
Among the substitutions that were detected in the study sequences (Table 6 and Table 7, and Table S5), two non-conservative substitutions E-T226A (T520A) and E-T231A (T525A) – detected in SG/EHI-CT1372 – reside within domain I, which is involved in E protein stability [74,75,76]. An E-T226R substitution was also shown to increase JEV neurovirulence and neuroinvasiveness in mice [77]; therefore, it may be useful to look into any effects of E-T226A on JEV virulence. Two E protein substitutions were also found at the neutralizing epitope position E-126 (residue 420 in the polyprotein) [76] across the study sequences in cluster A (E-I126V; Figure 2) and across the study sequences in cluster B (E-I126T; Figure 2) – these substitutions may allow for immune escape. However, no substitutions were detected in the study sequences at the E protein positions commonly found to be important for JEV neurovirulence and/or neuroinvasion (Table 6). Another important substitution was M-D40H (D167H), which was detected in SG/EHI-MS_CT261; a 40-residue region in flavivirus membrane (M) protein induces host cell apoptosis [78], and an M-I36F substitution within that region impeded JEV virion assembly in mammalian cells – likely due to the impairment of M protein apoptosis activity [79]. SG/EHI-MS_CT261 also carried an NS4b-A11T (A2283T) substitution that falls within the LIG-WD40 domain of the NS4b protein, and a mutant involving the LIG-WD40 domain has been shown to induce apoptosis and encephalitis in mouse brains [80]. Positive selection at the serine residue 24 of NS4b has been reported in genotype III of JEV that was isolated from pigs [81]. NS4b-P24S (polyprotein residue 2296) in SG/EHI-CX135 therefore warrants further investigations into its potential adaptation to swine—in particular, in the wild boar (Sus scrofa) as is found in Singapore. The substitution NS5-K68R (K2595R) in SG/EHI-MS_CT261 resides in the N-terminal 83-residue region of NS5, which is required for blocking IFN-α signaling [82]. However, NS5-K68R is a conservative change, and is less likely to affect the protein structure. Lastly, the most significant 3′-NCR mutations in the three whole genome sequences generated during this study were in pseudoknot 1 (PK1) of xrRNA1 (Table 5). PKs have been shown to form a mechanical block to exoribonuclease 1 by wrapping around the 5′ end of the xrRNA [83,84,85,86], hence mitigating host immune response [71,72,87]. Further research is needed to determine whether these mutations are important for adaptations to the local ecological niche.
The present study had several limitations. Firstly, the pooling of mosquitoes for screening might have confounded the detection of JEV in multiple mosquitoes of the same pool. Unfortunately, given the large number of specimens, it was not manpower- or cost-effective to screen individual mosquitoes. Secondly, the present study may not reflect the true prevalence of JE in Singapore due to the limited number of study sites, as the primary focus of this study was not to estimate the prevalence of JE in Singapore, but to determine whether there was any evidence of enzootic transmission of JEV in the country. Thirdly, although all sampling sites were selected based on the presence of ardeid birds, wild boars, water bodies, and vectors in their respective localities, vector control activities were strongly enforced in Park 2, 3, and 4 study sites (Figure 1) due to the proximity to human activity. This could affect the trapping of the mosquitoes at those sites. Lastly, the use of CDC light traps may have biased specimen collection. Lord, et al. [88] found that Cx. tritaeniorhynchus was oversampled with light traps in contrast to Cx. pseudovishnui and other potential JEV vectors and posited that this was due to the differences in the biting behaviors and host-seeking strategies rather than being indicative of the actual population density.

5. Conclusions

To the best of our knowledge, the present study reports the first isolation of JEV in Singapore since the abolishment of pig farms in 1992. This study is also the first record of whole genome sequences of JEV obtained from local mosquito pools. The phylogenetic analyses revealed the first description of genotype I of JEV in Singapore, as well as the possible re-establishment of sustained transmission of JEV in local fauna after its elimination in Singapore. We believe that these findings have important implications, not only for Asia and the Western Pacific, as per the recent Australian outbreak [9], but also for places outside of these regions—in particular, places where JEV is not yet enzootic/endemic but are inhabited by JEV vectors and their hosts. These areas do not necessarily need to be along the EAAF, as shown by the detection of JEV in Italy and Angola [89,90,91]. The unfolding effects of climate change tend to impact the size and the distribution of JEV vector populations, as well as vector biological properties [6,7,8], thus increasing the risk of JEV establishing transmission in such regions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14122662/s1, Figure S1: Maps of locations where (A) ardeid birds and (B) wild boars have been sighted. The maps were created with data from Wang and Hails [34] and Lamperty, Choik, Khoo, Amir, Baker, Chua, Chung, Chua, Koh, Lee, Lum, Pereira Mendes, Ngiam, O’Dempsey, Png, Sophie, Tan, Teo, Thomas, Li, Loo, Wardle and Luskin [32] respectively; Figure S2: Full JEV E-gene phylogenetic tree of 129 truncated JEV E-gene sequences; Table S1: JEV E-gene sequences used for the construction of the E-gene phylogeny tree; Table S2: Distribution of field-caught Culex mosquito populations; Table S3: Primers used for whole genome sequencing of JEV isolates; Table S4: JEV whole genome sequences used for the construction of the whole genome phylogeny tree; Table S5: Mutations in the present study JEV 3′-NCR secondary structures; Table S6: Full list of substitutions in JEV sequences.

Author Contributions

Conceptualization, D.M. and L.C.N.; methodology, D.M., M.K. and S.P.Y.; validation, D.M., H.C.H. and L.C.N.; formal analysis, M.J.L.; investigation, M.J.L., Z.Y.L., H.L.Y., H.Y.P., M.H.L., M.H., C.C., J.G. and M.K.; resources, D.M., H.C.H. and L.C.N.; data curation, M.J.L. and M.H.L.; writing—original draft preparation, M.J.L.; writing—review and editing, H.C.H., L.C.N. and D.M.; visualization, M.J.L.; supervision, D.M. and L.C.N.; project administration, D.M. and L.C.N.; funding acquisition, L.C.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data and datasets supporting the conclusions of this article are provided within the article or in the Supplementary Files.

Acknowledgments

We wish to thank the National Parks Board for their support in sample collection and NEA’s MSS for sample and knowledge provision. We would also like to thank our colleagues at the EHI for technical assistance and discussions, particularly Majhalia Torno, Germaine Ng, Jerald Yam, Carmen Koo, Hui Leng Quek, Siew Hwa Tan, Wei-Ping Tien, Joel Aik, Cheong Huat Tan, and Shuzhen Sim.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burke, D.S.; Leake, C.J. Japanese Encephalitis. In The Arboviruses: Epidemiology and Ecology, 1st ed.; Monath, T.P., Ed.; CRC Press: Boca Raton, FL, USA, 1988; Volume III, pp. 63–92. [Google Scholar]
  2. Yun, S.-I.; Lee, Y.-M. Japanese encephalitis. Hum. Vaccines Immunother. 2014, 10, 263–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, H.; Liang, G. Epidemiology of Japanese encephalitis: Past, present, and future prospects. Ther. Clin. Risk Manag. 2015, 11, 435–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Heffelfinger, J.D.; Li, X.; Batmunkh, N.; Grabovac, V.; Diorditsa, S.; Liyanage, J.B.; Pattamadilok, S.; Bahl, S.; Vannice, K.S.; Hyde, T.B.; et al. Japanese Encephalitis Surveillance and Immunization—Asia and Western Pacific Regions 2016. Morb. Mortal. Wkly. Rep. 2017, 66, 579–583. [Google Scholar] [CrossRef] [PubMed]
  5. Campbell, G.L.; Hills, S.L.; Fischer, M.; Jacobson, J.A.; Hoke, C.H.; Hombach, J.M.; Marfin, A.A.; Solomon, T.; Tsai, T.F.; Tsu, V.D.; et al. Estimated global incidence of Japanese encephalitis: A systematic review. Bull. World Health Organ. 2011, 89, 766–774E. [Google Scholar] [CrossRef]
  6. Hsu, S.M.; Yen, A.M.F.; Chen, T.H.H. The impact of climate on Japanese encephalitis. Epidemiol. Infect. 2008, 136, 980–987. [Google Scholar] [CrossRef]
  7. Bai, L.; Morton, L.C.; Liu, Q.-Y. Climate change and mosquito-borne diseases in China: A review. Glob. Health 2013, 9, 10. [Google Scholar] [CrossRef] [Green Version]
  8. Baylis, M.; Barker, C.M.; Caminade, C.; Joshid, B.R.; Pante, G.R.; Rayamajhi, A.; Reisen, W.K.; Impoinvil, D.E. Emergence or improved detection of Japanese encephalitis virus in the Himalayan highlands? Trans. R. Soc. Trop. Med. Hyg. 2016, 110, 209–211. [Google Scholar] [CrossRef] [Green Version]
  9. Japanese Encephalitis—Australia. Available online: https://www.who.int/emergencies/disease-outbreak-news/item/2022-DON365 (accessed on 13 June 2022).
  10. Chen, W.-R.; Tesh, R.B.; Rico-Hesse, R. Genetic variation of Japanese encephalitis virus in nature. J. Gen. Virol. 1990, 71, 2915–2922. [Google Scholar] [CrossRef]
  11. Chen, W.-R.; Rico-Hesse, R.; Tesh, R.B. A new genotype of Japanese encephalitis virus from Indonesia. Am. J. Trop. Med. Hyg. 1992, 47, 61–69. [Google Scholar] [CrossRef]
  12. Uchil, P.D.; Satchidanandam, V. Phylogenetic analysis of Japanese encephalitis virus: Envelope gene based analysis reveals a fifth genotype, geographic clustering, and multiple introductions of the virus into the Indian subcontinent. Am. J. Trop. Med. Hyg. 2001, 65, 242–251. [Google Scholar] [CrossRef]
  13. Li, M.-H.; Fu, S.-H.; Chen, W.-X.; Wang, H.-Y.W.; Guo, Y.-H.; Liu, Q.-Y.; Li, Y.-X.; Luo, H.-M.; Da, W.; Duo Ji, D.Z.; et al. Genotype V Japanese Encephalitis Virus Is Emerging. PLoS Negl. Trop. Dis. 2011, 5, e1231. [Google Scholar] [CrossRef]
  14. Woo, J.H.; Jeong, Y.E.; Jo, J.E.; Shim, S.-M.; Ryou, J.; Kim, K.-C.; Lee, W.J.; Lee, J.-Y. Genetic Characterization of Japanese Encephalitis Virus Genotype 5 Isolated from Patient, South Korea, 2015. Emerg. Infect. Dis. 2020, 26, 1002–1006. [Google Scholar] [CrossRef]
  15. Schuh, A.J.; Ward, M.J.; Leigh Brown, A.J.; Barrett, A.D.T. Phylogeography of Japanese Encephalitis Virus: Genotype Is Associated with Climate. PLoS Negl. Trop. Dis. 2013, 7, e2411. [Google Scholar] [CrossRef] [Green Version]
  16. Chambers, T.J.; Hahn, C.S.; Galler, R.; Rice, C.M. Flavivirus Genome Organization, Expression, and Replication. Annu. Rev. Microbiol. 1990, 44, 649–688. [Google Scholar] [CrossRef]
  17. Markoff, L. 5′- and 3′-noncoding regions in flavivirus RNA. Adv. Virus Res. 2003, 59, 177–228. [Google Scholar] [CrossRef]
  18. Seo, H.-J.; Kim, H.C.; Klein, T.A.; Ramey, A.M.; Lee, J.-H.; Kyung, S.-G.; Park, J.-Y.; Cho, Y.S.; Cho, I.-S.; Yeh, J.-Y. Molecular Detection and Genotyping of Japanese Encephalitis Virus in Mosquitoes during a 2010 Outbreak in the Republic of Korea. PLoS ONE 2013, 8, e55165. [Google Scholar] [CrossRef]
  19. Su, C.-L.; Yang, C.-F.; Teng, H.-J.; Lu, L.-C.; Lin, C.; Tsai, K.-H.; Chen, Y.-Y.; Chen, L.-Y.; Chang, S.-F.; Shu, P.-Y. Molecular Epidemiology of Japanese Encephalitis Virus in Mosquitoes in Taiwan during 2005–2012. PLoS Negl. Trop. Dis. 2014, 8, e3122. [Google Scholar] [CrossRef] [Green Version]
  20. Pearce, J.C.; Learoyd, T.P.; Langendorf, B.J.; Logan, J.G. Japanese encephalitis: The vectors, ecology and potential for expansion. J. Travel Med. 2018, 25, S16–S26. [Google Scholar] [CrossRef] [Green Version]
  21. Lindahl, J.F.; Ståhl, K.; Chirico, J.; Boqvist, S.; Thu, H.T.V.; Magnusson, U. Circulation of Japanese Encephalitis Virus in Pigs and Mosquito Vectors within Can Tho City, Vietnam. PLoS Negl. Trop. Dis. 2013, 4, e2153. [Google Scholar] [CrossRef]
  22. Vaughn, D.W.; Hoke, C.H., Jr. The epidemiology of Japanese encephalitis: Prospects for prevention. Epidemiol. Rev. 1992, 14, 197–221. [Google Scholar] [CrossRef]
  23. Pan, X.-L.; Liu, H.; Wang, H.-Y.; Fu, S.-H.; Liu, H.-Z.; Zhang, H.-L.; Li, M.-H.; Gao, X.-Y.; Wang, J.-L.; Sun, X.-H.; et al. Emergence of Genotype I of Japanese Encephalitis Virus as the Dominant Genotype in Asia. J. Virol. 2011, 85, 9847–9853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Schuh, A.J.; Ward, M.J.; Leigh Brown, A.J.; Barrett, A.D.T. Dynamics of the Emergence and Establishment of a Newly Dominant Genotype of Japanese Encephalitis Virus throughout Asia. J. Virol. 2014, 88, 4522–4532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Gao, X.-Y.; Liu, H.; Li, X.; Fu, S.-H.; Cao, L.; Shao, N.; Zhang, W.; Wang, Q.; Lu, Z.; Lei, W.; et al. Changing Geographic Distribution of Japanese Encephalitis Virus Genotypes, 1935–2017. Vector-Borne Zoonotic Dis. 2019, 19, 35–44. [Google Scholar] [CrossRef] [PubMed]
  26. Ting, S.H.L.; Tan, H.C.; Wong, W.K.; Ng, M.L.; Chan, S.H.; Ooi, E.E. Seroepidemiology of neutralizing antibodies to Japanese encephalitis virus in Singapore: Continued transmission despite abolishment of pig farming? Acta Trop. 2004, 92, 187–191. [Google Scholar] [CrossRef] [PubMed]
  27. Koh, Y.-L.; Tan, B.-H.; Loh, J.J.-P.; Ooi, E.E.; Su, S.-Y.; Hsu, L.-Y. Japanese Encephalitis, Singapore. Emerg. Infect. Dis. 2006, 12, 525–526. [Google Scholar] [CrossRef] [PubMed]
  28. Yee, A.T.K.; Corlett, R.T.; Liew, S.C.; Tan, H.T.W. The vegetation of Singapore—An updated map. Gard. Bull. Singap. 2011, 63, 205–212. [Google Scholar]
  29. Lam-Phua, S.-G.; Yeo, H.; Lee, R.M.-L.; Chong, C.-S.; Png, A.-B.; Foo, S.-Y.; Liew, C.; Ng, L.C.; Tang, C.-S.; Rueda, L.M.; et al. Mosquitoes (Diptera: Culicidae) of Singapore: Updated Checklist and New Records. J. Med. Entomol. 2019, 56, 103–119. [Google Scholar] [CrossRef] [PubMed]
  30. Yong, D.L.; Lee, B.P.Y.H.; Ang, A.L.; Tan, K.-H. The Status on Singapore Island of the Eurasian Wild Pig Sus Scrofa (Mammalia: Suidae). Nat. Singap. 2010, 3, 227–237. [Google Scholar]
  31. Koh, J.J.-M.; Rheindt, F.E.; Ng, Y.X.E.; Webb, E.L. Evidence of genetic connectivity between fragmented pig populations in a tropical urban city-state. Raffles Bull. Zool. 2019, 67, 14–31. [Google Scholar] [CrossRef]
  32. Lamperty, T.; Choik, W.; Khoo, M.; Amir, Z.; Baker, N.; Chua, M.; Chung, Y.; Chua, Y.; Koh, J.J.-M.; Lee, B.; et al. Rewilding by Large Terrestrial Vertebrates in Singapore; Manuscript submitted for publication; University of Washington: Washington, WA, USA, 2022. [Google Scholar]
  33. Lim, H.C.; Posa, M.R.C. Distribution and prey of migratory shorebirds on the northern coastline of Singapore. Raffles Bull. Zool. 2014, 62, 701–717. [Google Scholar]
  34. Wang, L.K.; Hails, C.J. An Annotated Checklist of the Birds of Singapore. Raffles Bull. Zool. 2007, 15, 1–179. [Google Scholar]
  35. Yap, G.; Lim, X.F.; Chan, S.; How, C.B.; Humaidi, M.; Yeo, G.; Mailepessov, D.; Kong, M.; Lai, Y.L.; Okumura, C.; et al. Serological evidence of continued Japanese encephalitis virus transmission in Singapore nearly three decades after end of pig farming. Parasites Vectors 2019, 12, 1–7. [Google Scholar] [CrossRef] [Green Version]
  36. Yap, G.; Mailepessov, D.; Lim, X.F.; Chan, S.; How, C.B.; Humaidi, M.; Yeo, G.; Chong, C.S.; Lam-Phua, S.G.; Lee, R.; et al. Detection of Japanese Encephalitis Virus in Culex Mosquitoes in Singapore. Am. J. Trop. Med. Hyg. 2020, 103, 1234–1240. [Google Scholar] [CrossRef]
  37. Yin-Coggrave, M.; Pong, W.S. A Study of the Mosquito Vectors of Viruses in Singapore. Ann. Trop. Med. Parasitol. 1964, 58, 359–366. [Google Scholar] [CrossRef]
  38. Pang, S.-C.; Andolina, C.; Malleret, B.; Christensen, P.R.; Lam-Phua, S.G.; Bin Abdul Razak, M.A.; Chong, C.S.; Li, D.; Chu, C.S.; Russell, B.; et al. Singapore’s Anopheles sinensis Form A is susceptible to Plasmodium vivax isolates from the western Thailand–Myanmar border. Malar. J. 2017, 16, 1–12. [Google Scholar] [CrossRef] [Green Version]
  39. Rattanarithikul, R.; Harrison, B.A.; Panthusiri, P.; Peyton, E.L.; Coleman, R.E. Illustrated keys to the mosquitoes of Thailand III. Genera Aedeomyia, Ficalbia, Mimomyia, Hodgesia, Coquillettidia, Mansonia, and Uranotaenia. Southeast Asian J. Trop. Med. Public Health 2006, 37 (Suppl. S1), 1–85. [Google Scholar]
  40. Rattanarithikul, R.; Harrison, B.A.; Panthusiri, P.; Coleman, R.E. Illustrated keys to the mosquitoes of Thailand I. Background; geographic distribution; lists of genera, subgenera, and species; and a key to the genera. Southeast Asian J. Trop. Med. Public Health 2005, 36 (Suppl. S1), 1–80. [Google Scholar]
  41. Rattanarithikul, R.; Harrison, B.A.; Harbach, R.E.; Panthusiri, P.; Coleman, R.E. Illustrated keys to the mosquitoes of Thailand IV. Anopheles. Southeast Asian J. Trop. Med. Public Health 2006, 37 (Suppl. S2), 1–128. [Google Scholar]
  42. Rattanarithikul, R.; Harbach, R.E.; Harrison, B.A.; Panthusiri, P.; Jones, J.W.; Coleman, R.E. Illustrated keys to the mosquitoes of Thailand. II. Genera Culex and Lutzia. Southeast Asian J. Trop. Med. Public Health 2005, 36 (Suppl. S2), 1–97. [Google Scholar]
  43. Rattanarithikul, R.; Harbach, R.E.; Harrison, B.A.; Panthusiri, P.; Coleman, R.E.; Richardson, J.H. Illustrated keys to the mosquitoes of Thailand VI. Tribe aedini. Southeast Asian J. Trop. Med. Public Health 2010, 41 (Suppl. S1), 1–225. [Google Scholar]
  44. Rattanarithikul, R.; Harbach, R.E.; Harrison, B.A.; Panthusiri, P.; Coleman, R.E. Illustrated keys to themosquitoes of Thailand V.Genera orthopodomyia, kimia, Malaya, topomyia, tripteroides, and Toxorhynchites. Southeast Asian J. Trop. Med. Public Health 2007, 38 ((Suppl. S2), 1–65. [Google Scholar] [PubMed]
  45. Yeo, H.; Tan, C.-H.; Chong, C.S.; Lam-Phua, S.G. Identification key to the genera of adult female mosquitoes (Diptera: Culicidae) of Singapore. Zootaxa 2019, 4624, 407–423. [Google Scholar] [CrossRef] [PubMed]
  46. Scherret, J.H.; Poidinger, M.; Mackenzie, J.S.; Broom, A.K.; Deubel, V.; Lipkin, W.I.; Briese, T.; Gould, E.A.; Hall, R.A. The Relationships between West Nile and Kunjin Viruses. Emerg. Infect. Dis. 2001, 7, 697–705. [Google Scholar] [CrossRef] [PubMed]
  47. Santhosh, S.R.; Paridaa, M.M.; Dash, P.K.; Pateriya, A.; Pattnaik, B.; Pradhan, H.K.; Tripathi, N.K.; Ambuj, S.; Gupta, N.; Saxena, P.; et al. Development and evaluation of SYBR Green I-based one-step real-time RT-PCR assay for detection and quantitation of Japanese encephalitis virus. J. Virol. Methods 2007, 143, 73–80. [Google Scholar] [CrossRef] [PubMed]
  48. Schuh, A.J.; Li, L.; Tesh, R.B.; Innis, B.L.; Barrett, A.D.T. 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. [Google Scholar] [CrossRef] [PubMed]
  49. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  50. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar] [CrossRef] [Green Version]
  51. Madeira, F.; Pearce, M.; Tivey, A.R.N.; Basutkar, P.; Lee, J.; Edbali, O.; Madhusoodanan, N.; Kolesnikov, A.; Lopez, R. Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 2022; online ahead of print. [Google Scholar] [CrossRef]
  52. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef]
  53. Duvaud, S.; Gabella, C.; Lisacek, F.; Stockinger, H.; Ioannidis, V.; Durinx, C. Expasy, the Swiss Bioinformatics Resource Portal, as designed by its users. Nucleic Acids Res. 2021, 49, W216–W227. [Google Scholar] [CrossRef]
  54. Stothard, P. The Sequence Manipulation Suite: JavaScript Programs for Analyzing and Formatting Protein and DNA Sequences. BioTechniques 2000, 28, 1102–1104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef] [PubMed]
  56. Liu, H.; Zhang, J.; Niu, Y.; Liang, G. The 5′ and 3′ Untranslated Regions of the Japanese Encephalitis Virus (JEV): Molecular Genetics and Higher Order Structures. Front. Microbiol. 2021, 12, 730045. [Google Scholar] [CrossRef]
  57. Hu, T.; Wu, Z.; Wu, S.; Chen, S.; Cheng, A. The key amino acids of E protein involved in early flavivirus infection: Viral entry. Virol. J. 2021, 18, 1–12. [Google Scholar] [CrossRef] [PubMed]
  58. Liu, X.; Zhao, X.; Na, R.; Li, L.; Warkentin, E.; Witt, J.; Lu, X.; Yu, Y.; Wei, Y.; Peng, G.; et al. The structure differences of Japanese encephalitis virus SA14 and SA14-14-2 E proteins elucidate the virulence attenuation mechanism. Protein Cell 2019, 10, 149–153. [Google Scholar] [CrossRef] [Green Version]
  59. Yun, S.-M.; Cho, J.E.; Ju, Y.-R.; Kim, S.Y.; Ryou, J.; Han, M.G.; Choi, W.-Y.; Jeong, Y.E. Molecular epidemiology of Japanese encephalitis virus circulating in South Korea, 1983–2005. Virol. J. 2010, 7, 127. [Google Scholar] [CrossRef] [Green Version]
  60. Gromowski, G.D.; Firestone, C.-Y.; Whitehead, S.S. Genetic Determinants of Japanese Encephalitis Virus Vaccine Strain SA14-14-2 That Govern Attenuation of Virulence in Mice. J. Virol. 2015, 89, 6328–6337. [Google Scholar] [CrossRef] [Green Version]
  61. Yamaguchi, Y.; Nukui, Y.; Kotaki, A.; Sawabe, K.; Saijo, M.; Watanabe, H.; Kurane, I.; Takasaki, T.; Tajima, S. Characterization of a serine-to-asparagine substitution at position 123 in the Japanese encephalitis virus E protein. J. General Virol. 2013, 94, 90–96. [Google Scholar] [CrossRef]
  62. 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. [Google Scholar] [CrossRef]
  63. 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. General Virol. 2005, 86, 2209–2220. [Google Scholar] [CrossRef]
  64. McMinn, P.C. The molecular basis of virulence of the encephalitogenic flaviviruses. J. General Virol. 1997, 78, 2711–2722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Obara, M.; Yamauchi, T.; Watanabe, M.; Hasegawa, S.; Ueda, Y.; Matsuno, K.; Iwai, M.; Horimoto, E.; Kurata, T.; Takizawa, T.; et al. Continuity and Change of Japanese Encephalitis Virus in Toyama Prefecture, Japan. Am. J. Trop. Med. Hyg. 2011, 84, 695–708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kuwata, R.; Torii, S.; Shimoda, H.; Supriyono; Phichitraslip, T.; Prasertsincharoen, N.; Takemae, H.; Bautista, R.C.J.T.; Ebora, V.D.B.M.; Abella, J.A.C.; et al. Distribution of Japanese Encephalitis Virus, Japan and Southeast Asia, 2016–2018. Emerg. Infect. Dis. 2020, 26, 125–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Wong, C.M.L. The Developmental State in Ecological Modernization and the Politics of Environmental Framings: The Case of Singapore and Implications for East Asia. Nat. Cult. 2012, 7, 95–119. [Google Scholar] [CrossRef]
  68. Richards, D.R.; Law, A.; Tan, C.S.Y.; Shaikh, S.F.E.A.; Carrasco, L.R.; Jaung, W.; Oh, R.R.Y. Rapid urbanisation in Singapore causes a shift from local provisioning and regulating to cultural Ecosyst. Serv. use. Ecosyst. Serv. 2020, 46, 101193. [Google Scholar] [CrossRef]
  69. Yeo, G.; Chan, S.; How, C.B.; Humaidi, M.; Lim, X.F.; Mailepessov, D.; Chong, C.S.; Lam-Phua, S.G.; Lee, R.; Hapuarachchi, H.C.; et al. Molecular Analysis of the Bloodmeals of Culex spp. Mosquitoes at Natural Habitats in Singapore to Investigate the Potential Risk of Japanese Encephalitis Virus and West Nile Virus Transmission. Vector-Borne Zoonotic Dis. 2020, 20, 703–714. [Google Scholar] [CrossRef]
  70. Bae, W.; Kim, J.H.; Kim, J.; Lee, J.; Hwang, E.S. Changes of Epidemiological Characteristics of Japanese Encephalitis Viral Infection and Birds as a Potential Viral Transmitter in Korea. J. Korean Med. Sci. 2018, 33, e70. [Google Scholar] [CrossRef]
  71. Clarke, B.D.; Roby, J.A.; Slonchak, A.; Khromykh, A.A. Functional non-coding RNAs derived from the flavivirus 3′ untranslated region. Virus Res. 2015, 206, 53–61. [Google Scholar] [CrossRef]
  72. Charley, P.A.; Wilusz, J. Standing your ground to exoribonucleases: Function of Flavivirus long non-coding RNAs. Virus Res. 2016, 212, 70–77. [Google Scholar] [CrossRef]
  73. Sanford, T.J.; Mears, H.V.; Fajardo, T., Jr.; Locker, N.; Sweeney, T.R. Circularization of flavivirus genomic RNA inhibits de novo translation initiation. Nucleic Acids Res. 2019, 47, 9789–9802. [Google Scholar] [CrossRef] [Green Version]
  74. Wu, K.-P.; Wu, C.-W.; Tsao, Y.-P.; Kuo, T.-W.; Lou, Y.-C.; Lin, C.-W.; Wu, S.-C.; Cheng, J.-W. Structural basis of a flavivirus recognized by its neutralizing antibody: Solution structure of the domain III of the Japanese encephalitis virus envelope protein. J. Biol. Chem. 2003, 278, 46007–46013. [Google Scholar] [CrossRef] [Green Version]
  75. Zhang, Y.; Zhang, W.; Ogata, S.; Clements, D.; Strauss, J.H.; Baker, T.S.; Kuhn, R.J.; Rossmann, M.G. Conformational Changes of the Flavivirus E Glycoprotein. Structure 2004, 12, 1607–1618. [Google Scholar] [CrossRef] [Green Version]
  76. Luca, V.C.; AbiMansour, J.; Nelson, C.A.; Fremont, D.H. Crystal Structure of the Japanese Encephalitis Virus Envelope Protein. J. Virol. 2012, 86, 2337–2346. [Google Scholar] [CrossRef] [Green Version]
  77. Yang, D.; Li, X.-F.; Ye, Q.; Wang, H.-J.; Deng, Y.-Q.; Zhu, S.-Y.; Zhang, Y.; Li, S.-H.; Qin, C.-F. Characterization of live-attenuated Japanese encephalitis vaccine virus SA14-14-2. Vaccine 2014, 32, 2675–2681. [Google Scholar] [CrossRef]
  78. Catteau, A.; Kalinina, O.; Wagner, M.-C.; Deubel, V.; Courageot, M.-P.; Desprès, P. Dengue virus M protein contains a proapoptotic sequence referred to as ApoptoM. J. Gen. Virol. 2003, 84, 2781–2793. [Google Scholar] [CrossRef]
  79. De Wispelaere, M.; Khou, C.; Frenkiel, M.-P.; Desprès, P.; Pardigon, N. A Single Amino Acid Substitution in the M Protein Attenuates Japanese Encephalitis Virus in Mammalian Hosts. J. Virol. 2016, 90, 2676–2689. [Google Scholar] [CrossRef] [Green Version]
  80. Wang, Q.; Xin, X.; Wang, T.; Wan, J.; Ou, Y.; Yang, Z.; Yu, Q.; Zhu, L.; Guo, Y.; Wu, Y.; et al. Japanese Encephalitis Virus Induces Apoptosis and Encephalitis by Activating the PERK Pathway. J. Virol. 2019, 73, e00887-19. [Google Scholar] [CrossRef] [Green Version]
  81. Liu, W.-J.; Zhu, M.; Pei, J.-J.; Dong, X.-Y.; Liu, W.; Zhao, M.-Q.; Wang, J.-Y.; Gou, H.-C.; Luo, Y.-W.; Chen, J.-D. Molecular phylogenetic and positive selection analysis of Japanese encephalitis virus strains isolated from pigs in China. Virus Res. 2013, 178, 547–552. [Google Scholar] [CrossRef]
  82. Lin, R.-J.; Chang, B.-L.; Yu, H.-P.; Liao, C.-L.; Lin, Y.-L. Blocking of Interferon-Induced Jak-Stat Signaling by Japanese Encephalitis Virus NS5 through a Protein Tyrosine Phosphatase-Mediated Mechanism. J. Virol. 2006, 80, 5908–5918. [Google Scholar] [CrossRef]
  83. Chapman, E.G.; Costantino, D.A.; Rabe, J.L.; Moon, S.L.; Wilusz, J.; Nix, J.C.; Kieft, J.S. The Structural Basis of Pathogenic Subgenomic Flavivirus RNA (sfRNA) Production. Science 2014, 344, 307–310. [Google Scholar] [CrossRef] [Green Version]
  84. Kieft, J.S.; Rabe, J.L.; Chapman, E.G. New hypotheses derived from the structure of a flaviviral Xrn1-resistant RNA: Conservation, folding, and host adaptation. RNA Biol. 2015, 12, 1169–1177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Akiyama, B.M.; Laurence, H.M.; Massey, A.R.; Costantino, D.A.; Xie, X.; Yang, Y.; Shi, P.-Y.; Nix, J.C.; Beckham, J.D.; Kieft, J.S. Zika virus produces noncoding RNAs using a multi-pseudoknot structure that confounds a cellular exonuclease. Science 2016, 354, 1148–1152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. MacFadden, A.; O’Donoghue, Z.; Silva, P.A.G.C.; Chapman, E.G.; Olsthoorn, R.C.; Sterken, M.G.; Pijlman, G.; Bredenbeek, P.P.J.; Kieft, J.S. Mechanism and structural diversity of exoribonuclease-resistant RNA structures in flaviviral RNAs. Nat. Commun. 2018, 9, 119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Niu, X.; Sun, R.; Chen, Z.; Yao, Y.; Zuo, X.; Chen, C.; Fang, X. Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs. Nat. Commun. 2021, 12, 6417. [Google Scholar] [CrossRef]
  88. Lord, J.S.; Al-Amin, H.M.; Chakma, S.; Alam, M.S.; Gurley, E.S.; Pulliam, J.R.C. Sampling Design Influences the Observed Dominance of Culex tritaeniorhynchus: Considerations for Future Studies of Japanese Encephalitis Virus Transmission. PLoS Negl. Trop. Dis. 2016, 10, e0004249. [Google Scholar] [CrossRef] [Green Version]
  89. Platonov, A.E.; Rossi, G.; Karan, L.S.; Mironov, K.O.; Busani, L.; Rezza, G. Does the Japanese encephalitis virus (JEV) represent a threat for human health in Europe? Detection of JEV RNA sequences in birds collected in Italy. Eurosurveillance 2012, 17, 20241. [Google Scholar] [CrossRef] [Green Version]
  90. Ravanini, P.; Huhtamo, E.; Ilaria, V.; Crobu, M.G.; Nicosia, A.M.; Servino, L.; Rivasi, F.; Allegrini, S.; Miglio, U.; Magri, A.; et al. Japanese encephalitis virus RNA detected in Culex pipiens mosquitoes in Italy. Eurosurveillance 2012, 17, 20221. [Google Scholar] [CrossRef] [Green Version]
  91. Simon-Loriere, E.; Faye, O.; Prot, M.; Casademont, I.; Fall, G.; Fernandez-Garcia, M.D.; Diagne, M.M.; Kipela, J.-M.; Fall, I.S.; Holmes, E.C.; et al. Autochthonous Japanese Encephalitis with Yellow Fever Coinfection in Africa. N. Engl. J. Med. 2017, 376, 1483–1485. [Google Scholar] [CrossRef]
Figure 1. Locations of the five sites used in JEV entomological surveillance in this study, as well as the two locations where JEV was detected in mosquitoes that were obtained from MSS. The single-site study site is in purple and the multi-site study sites are in light blue. MSS locations are in yellow. The solid shapes indicate areas where JEV-positive pools were observed, while the hatched shapes indicate areas where JEV-positive pools were not observed. See Figure S1 for the approximate locations where ardeid birds and wild boars were sighted in Singapore.
Figure 1. Locations of the five sites used in JEV entomological surveillance in this study, as well as the two locations where JEV was detected in mosquitoes that were obtained from MSS. The single-site study site is in purple and the multi-site study sites are in light blue. MSS locations are in yellow. The solid shapes indicate areas where JEV-positive pools were observed, while the hatched shapes indicate areas where JEV-positive pools were not observed. See Figure S1 for the approximate locations where ardeid birds and wild boars were sighted in Singapore.
Viruses 14 02662 g001
Figure 2. Maximum likelihood phylogenetic analysis of 129 truncated (701 bp) JEV E-gene sequences. This tree was constructed using TN93+G with 1000 bootstrap iterations; the outgroup is MVEV (NCBI Reference Sequence no. NC_000943.1). Only parts of the tree are presented for easier viewing (see Figure S1 for the full tree): (A) Sub-genotype I-a; (B) Partial sub-genotype I-b. The JEV E-gene sequences presented in this study are indicated with cluster A (blue) and cluster B (orange), as well as a star for SG/EHI-CT1372. Each node is labelled with the GenBank Accession No. (where available) followed by the sequence name.
Figure 2. Maximum likelihood phylogenetic analysis of 129 truncated (701 bp) JEV E-gene sequences. This tree was constructed using TN93+G with 1000 bootstrap iterations; the outgroup is MVEV (NCBI Reference Sequence no. NC_000943.1). Only parts of the tree are presented for easier viewing (see Figure S1 for the full tree): (A) Sub-genotype I-a; (B) Partial sub-genotype I-b. The JEV E-gene sequences presented in this study are indicated with cluster A (blue) and cluster B (orange), as well as a star for SG/EHI-CT1372. Each node is labelled with the GenBank Accession No. (where available) followed by the sequence name.
Viruses 14 02662 g002aViruses 14 02662 g002b
Figure 3. Maximum likelihood phylogenetic analysis of 50 JEV whole genome isolates. This tree was constructed using TN93+G+I with 1000 bootstrap iterations; the outgroup is MVEV (NCBI Reference Sequence no. NC_000943.1). The JEV isolates presented in this study are indicated with a star: SG/EHI-CT1372; blue triangle: SG/EHI-CX135; orange circle: SG/EHI-MS_CT261. The colors correspond to the clusters indicated in Figure 2. Genotypes are as indicated. Each node is labelled with the GenBank Accession No. (where available) followed by the isolate name.
Figure 3. Maximum likelihood phylogenetic analysis of 50 JEV whole genome isolates. This tree was constructed using TN93+G+I with 1000 bootstrap iterations; the outgroup is MVEV (NCBI Reference Sequence no. NC_000943.1). The JEV isolates presented in this study are indicated with a star: SG/EHI-CT1372; blue triangle: SG/EHI-CX135; orange circle: SG/EHI-MS_CT261. The colors correspond to the clusters indicated in Figure 2. Genotypes are as indicated. Each node is labelled with the GenBank Accession No. (where available) followed by the isolate name.
Viruses 14 02662 g003
Table 1. Field-caught Culex mosquitoes in the present study. See Table S2 for more details on the distribution of mosquitoes in the present study.
Table 1. Field-caught Culex mosquitoes in the present study. See Table S2 for more details on the distribution of mosquitoes in the present study.
ClassificationCDC Light TrapsNight Catcher Traps
No. of Mosquitoes CollectedNo. of Mosquito PoolsNo. of JEV-Positive PoolsDetection Percentage 1 (%)No. of Mosquitoes CollectedNo. of Mosquito PoolsNo. of JEV-Positive PoolsDetection Percentage 1 (%)
Culex spp. 213,84521352.35----
Cx. bitaeniorhynchus33170032400
Cx. gelidus11467002041000
Cx. quinquefasciatus288293001100
Cx. sitiens1497535002491200
Cx. tritaeniorhynchus14,25892620.2210,22329131.03
Culex spp. vishnui subgroup 3964610210013247100
1 JEV-positive pools/total pools. 2 From September to October 2019 only. 3 Excluding Cx. tritaeniorhynchus.
Table 2. Spatial and temporal information of the JEV-positive pools.
Table 2. Spatial and temporal information of the JEV-positive pools.
JEV Sequence IDLocation CollectedMonth & Year CollectedTrap Used for Collection
SG/EHI-M_CT21_45MSS 2May 2021Night Catcher trap
SG/EHI-MS_CT310MSS 1December 2020Night Catcher trap
SG/EHI-MS_CT261MSS 1November 2020Night Catcher trap
SG/EHI-CT1372NR 1November 2019CDC light trap
SG/EHI-CX135NR 1October 2019CDC light trap
SG/EHI-CX193NR 1October 2019CDC light trap
SG/EHI-CX 202NR 1October 2019CDC light trap
SG/EHI-CX 209NR 1October 2019CDC light trap
SG/EHI-CX 214NR 1October 2019CDC light trap
SG/EHI-CT1710Park 1June 2014CDC light trap
Table 3. Percentage nucleotide and amino acid sequence similarity of the envelope gene sequences obtained from ten mosquito pools. Nucleotide sequence identities (NT ID) are in the upper right triangular area above the shaded cells and amino acid sequence identities (AA ID) are in the lower left triangular area.
Table 3. Percentage nucleotide and amino acid sequence similarity of the envelope gene sequences obtained from ten mosquito pools. Nucleotide sequence identities (NT ID) are in the upper right triangular area above the shaded cells and amino acid sequence identities (AA ID) are in the lower left triangular area.
NT ID (%)SG/EHI-CT1372SG/EHI-CX 202SG/EHI-CX 209SG/EHI-CX 214SG/EHI-CX193SG/EHI-CX135SG/EHI-CT1710SG/EHI-MS_CT261SG/EHI-MS_CT310SG/EHI-M_CT21_45
AA ID (%)
SG/EHI-CT1372-91.8791.8791.7391.7391.5892.1592.7292.7292.58
SG/EHI-CX 20296.57-10099.5799.5799.4394.1594.4494.4494.29
SG/EHI-CX 20996.57100-99.5799.5799.4394.1594.4494.4494.29
SG/EHI-CX 21496.57100100-10099.8694.2994.5894.5894.44
SG/EHI-CX19396.57100100100-99.8694.2994.5894.5894.44
SG/EHI-CX13596.57100100100100-94.1594.4494.4494.29
SG/EHI-CT171096.5796.1496.1496.1496.1496.14-99.1499.1499
SG/EHI-MS_CT26196.5796.1496.1496.1496.1496.14100-10099.86
SG/EHI-MS_CT31096.5796.1496.1496.1496.1496.14100100-99.86
SG/EHI-M_CT21_4596.1495.7195.7195.7195.7195.7199.5799.5799.57-
Table 4. Percentage nucleotide and amino acid sequence similarity of the polyprotein sequences of three study isolates. Nucleotide sequence identities (NT ID) are in the upper right triangular area above the shaded cells and amino acid sequence identities (AA ID) are in the lower left triangular area.
Table 4. Percentage nucleotide and amino acid sequence similarity of the polyprotein sequences of three study isolates. Nucleotide sequence identities (NT ID) are in the upper right triangular area above the shaded cells and amino acid sequence identities (AA ID) are in the lower left triangular area.
NT ID (%)SG/EHI-CT1372SG/EHI-CX135SG/EHI-MS_CT261
AA ID (%)
SG/EHI-CT1372-93.0192.99
SG/EHI-CX13595.47-93.33
SG/EHI-MS_CT26194.3996.75-
Table 5. Mutations involved in significant regions of 3′-NCR secondary structures belonging to the JEV isolates in this study.
Table 5. Mutations involved in significant regions of 3′-NCR secondary structures belonging to the JEV isolates in this study.
IsolateMutationRegion Involved
SG/EHI-CT1372G > APseudoknot 1 (PK1), found in exoribonuclease resistant RNA 1 (xrRNA1)
SG/EHI-CX135, SG/EHI-MS_CT261 T > C
SG/EHI-CX135C > T3′-NCR variable region conserved sequences 1 (3′ vrCS1), found in the 3′-NCR variable region stem loop (3′vrSL)
SG/EHI-CX135G > A
SG/EHI-CX135, SG/EHI-MS_CT261T > C3′-NCR dumbbell structure conserved sequences 2 (3′dbCS2), found in the dumbbell structure 2 (DB2)
Table 6. Novel amino acid substitutions in the structural polyprotein of JEV genotype I. Substitutions that are only detected in the study sequences are shown. The comparison included all of the JEV GI complete coding sequences (n = 158) available in NCBI as of November 2022. Amino acid positions in the E protein that are commonly associated with JEV neurovirulence and/or neuroinvasiveness are also highlighted. Amino acid residues are numbered according to the complete polyprotein of NCBI Reference Sequence NP_059434.1. Dots indicate the residues that are identical to the consensus sequence. See Table S6 for the full list of substitutions in JEV sequences included in the table.
Table 6. Novel amino acid substitutions in the structural polyprotein of JEV genotype I. Substitutions that are only detected in the study sequences are shown. The comparison included all of the JEV GI complete coding sequences (n = 158) available in NCBI as of November 2022. Amino acid positions in the E protein that are commonly associated with JEV neurovirulence and/or neuroinvasiveness are also highlighted. Amino acid residues are numbered according to the complete polyprotein of NCBI Reference Sequence NP_059434.1. Dots indicate the residues that are identical to the consensus sequence. See Table S6 for the full list of substitutions in JEV sequences included in the table.
ProteinCapsidprM/MembraneEnvelope
Amino Acid Position in the Protein36476396107 123 138 176 177 226231264 279 360447492
Polyprotein Amino Acid Position3647190390401417432470471520525558573654741786
JEV GI consensusSVEFLSEITTTQKFGV
SG/EHI-CT1372 L A *A * L
SG/EHI-CX135G *A Y *S *
SG/EHI-CX 202NDNDND NDND
SG/EHI-CX 209NDNDND NDND
SG/EHI-CX 214NDNDND NDND
SG/EHI-CX193NDNDND NDND
SG/EHI-MS_CT261 K I
SG/EHI-CT1710NDNDND NDND
SG/EHI-MS_CT310NDNDND NDND
SG/EHI-M_CT21_45NDNDND NDND
ND = Not determined because only partial E-gene sequences (701 bp) were analyzed for respective sequences. * Non-conservative (polar/non-polar) mutation; Conservative mutation with charge differences. Amino acid positions in the E protein commonly associated with JEV neurovirulence and/or neuroinvasiveness [57,58,59,60,61,62,63,64].
Table 7. Novel amino acid substitutions in the non-structural polyprotein of JEV genotype I. Substitutions only detected in the study sequences are shown. The comparison included all of the JEV GI complete coding sequences (n = 158) available in NCBI as of November 2022. Amino acid residues are numbered according to the complete polyprotein of NCBI Reference Sequence NP_059434.1. Dots indicate the residues that are identical to the consensus sequence. See Table S6 for the full list of substitutions in isolates included in the table and the undetermined amino acid residues of SG/EHI-CX135.
Table 7. Novel amino acid substitutions in the non-structural polyprotein of JEV genotype I. Substitutions only detected in the study sequences are shown. The comparison included all of the JEV GI complete coding sequences (n = 158) available in NCBI as of November 2022. Amino acid residues are numbered according to the complete polyprotein of NCBI Reference Sequence NP_059434.1. Dots indicate the residues that are identical to the consensus sequence. See Table S6 for the full list of substitutions in isolates included in the table and the undetermined amino acid residues of SG/EHI-CX135.
ProteinNS1NS2ANS3NS4ANS4BNS5
Amino Acid Position in the Protein917720821062778433635840254117811918368168280426805830
Polyprotein Amino Acid Position8039711002135615661581158818401862190621772283235023912455259526952807295333323357
JEV GI consensus TDDSEIRLKDIASTVKKREKM
SG/EHI-CT1372I * A *D I *L R
SG/EHI-CX135 ENDND E NDND D A
SG/EHI-MS_CT261 E LKMR VT *A *M * RRM * T *
ND = Not determined because of a poor quality sequence. * Non-conservative (polar/non-polar) mutations.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lim, M.J.; Loh, Z.Y.; Yeo, H.L.; Yenamandra, S.P.; Kong, M.; Pang, H.Y.; Lee, M.H.; Humaidi, M.; Chua, C.; Griffiths, J.; et al. Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore. Viruses 2022, 14, 2662. https://doi.org/10.3390/v14122662

AMA Style

Lim MJ, Loh ZY, Yeo HL, Yenamandra SP, Kong M, Pang HY, Lee MH, Humaidi M, Chua C, Griffiths J, et al. Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore. Viruses. 2022; 14(12):2662. https://doi.org/10.3390/v14122662

Chicago/Turabian Style

Lim, Ming Jie, Zhi Yang Loh, Hui Ling Yeo, Surya Pavan Yenamandra, Marcella Kong, Hao Yang Pang, Meng Han Lee, Mahathir Humaidi, Cliff Chua, Jane Griffiths, and et al. 2022. "Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore" Viruses 14, no. 12: 2662. https://doi.org/10.3390/v14122662

APA Style

Lim, M. J., Loh, Z. Y., Yeo, H. L., Yenamandra, S. P., Kong, M., Pang, H. Y., Lee, M. H., Humaidi, M., Chua, C., Griffiths, J., Ng, L. C., Hapuarachchi, H. C., & Mailepessov, D. (2022). Isolation and Genetic Characterization of Japanese Encephalitis Virus Two Decades after Its Elimination in Singapore. Viruses, 14(12), 2662. https://doi.org/10.3390/v14122662

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop