JIA Qing Jun, CHEN Xin Yu, LI De Zhou, XU Juan Juan, XU Zhi Gang, DUAN Zhi Liang, and WEN Jin Sheng,#
Comparative Genom ic Analysis of Enterovirus 71 Revealed Six New Potential Neurovirulence-associated Sites*
JIA Qing Jun1,?, CHEN Xin Yu2,?, LI De Zhou3, XU Juan Juan1, XU Zhi Gang1, DUAN Zhi Liang2,#, and WEN Jin Sheng1,#
In the present study, the complete genomes of four common (4/EV71/Wenzhou/CHN/2014, 15/ EV71/Wenzhou/CHN/2014, 116/EV71/Wenzhou/ CHN/2014, and 120/EV71/Wenzhou/CHN/2014) and two virulent (11/EV71/Wenzhou/CHN/2014 and 109/EV71/Wenzhou/CHN/2014) enterovirus 71 (EV71) isolates were sequenced and described. They are 7405 bp in length and belong to EV71 sub-genotype C4 (C4a cluster). Nucleotide sequence alignment revealed six nucleotide variations (GP151→TP151, GP199→AP199, GP261→TP261, AP328→CP328, GP422→AP422, and GP437→TP437) in the two virulent isolates w ithin the 5'UTR of the IRES element. RNA secondary structure predictions of IRES and FCE indicated that the common isolates shared sim ilar structures, which were different from those of the virulent isolates. Moreover, the GP114→CP114and GP151→TP151mutations in the virulent isolates contributed to the formation of the unique RNA secondary structures in SL II. Furthermore, nucleotide/am ino acid sequence alignments of 82 EV71 isolates indicated that six sites (TP488and CP577in the 5'UTR; AsnP57in 2A; IleP56in 3C; CP10and AP47in the 3'UTR) are potentially associated w ith the neurovirulence of EV71. Finally, the 3D structures of 2A were analogous, whereas the structures of VP1 and 3C were variable.
Enterovirus 71 (EV71) is a non-enveloped virus containing a positive-sense single-stranded RNA genome (approximately 7.4 kb) and belongs to the Enterovirus genus of the fam ily Picornaviridae. Structurally, six stem loop (SL) structures, SL I to SL VI, are contained in the EV71 5'UTR. Functionally, the 5'UTR is divided into three regions: a 5' term inal cloverleaf (1 to 89-101 bp), an internal ribosome entry site (IRES) element (123-126 bp to 602-605 bp, ranging from SL II to SL VI), and a hypervariable region (the last 120 bp)[1]. This type of positive-sense single-stranded RNA virus utilizes a cap-independent translation initiation mechanism to regulate the expression of its viral genome[2]. Jerrey et al. found that the functional cis-acting element (FCE, spanning SL IV to SL VI w ithin the 5'UTR) enabled Poliovirus mRNA to translate in a cap-independent fashion and that the major determ inant of this FCE is located in a domain(nts320-631)[3].Unlikeother hand-foot-mouth disease (HFMD) pathogens [such as Coxsackie virus A (CVA)], EV71 infection is more likely to cause severe neurological complications, includingbrainstemencephalitis,aseptic encephalitis, meningitis, poliomyelitis-like paralysis, neurogenic pulmonary edema, and even death[4]. Some studies have reported that the nucleotide mutations w ithin the 5'UTR (such as the stem loop II, SL II, nts 105-179) and some am ino acid variations are associated w ith EV71 neurovirulence[5-6]. Generally, common EV71 isolates cause m ild and self-lim ited herpes outbreaks on the skin and mucosae, whereas virulent EV71 isolates (‘virulent' particularly refers to neurovirulent) lead to severe neurological complications and even death[7].
In the present study, patients diagnosed with either m ild HFMD (with herpes only on the skin and mucosa) or severe HFMD (w ith both herpes and meningitis) were recruited from the Second Affiliated Hospital of Wenzhou Medical University. Throat swab samples or cerebrospinal fluid sampleswere filtered through a 0.22-μm sterile filter, and the filtered suspension was added to either RD cells or Vero cells, which were maintained in RPMI-1640 medium supplemented with 10% FBS, 100 U/m L penicillin, and 100 μg/m L streptomycin. After an infection period of 1.5 h in a 37 °C/5% CO2incubator, the cell supernatant was discarded, and fresh RPMI-1640 medium w ithout FBS was added to the cells. After a 4-d incubation, the supernatant was blind inoculated into freshly prepared cells several times until the cells exhibited a cytopathic effect (CPE), as scored by ‘+++' (the majority of cells had shrunk and become dislodged from the bottom of the culture bottles). Based on EV71-specific primers (TableS1inthewebsiteofBES, www.besjournal.com),RT-PCRassayswere conducted to amplify six overlapping gene segments (638 bp, 1575 bp, 1433 bp, 1624 bp, 2286 bp, and 786 bp), which span the complete genome of EV71, and the PCR products were sequenced. The complete genome sequences and deduced am ino acid sequences of six EV71 isolates were deposited in the GenBank database. The GenBank Accession Numbers of four common isolates (4/EV71/ Wenzhou/CHN/2014, 15/EV71/Wenzhou/CHN/2014, 116/EV71/Wenzhou/CHN/2014, and 120/EV71/ Wenzhou/CHN/2014) were KT008669, KT345960, KT008672, and KT345959, respectively. Two virulent isolates(11/EV71/Wenzhou/CHN/2014and 109/EV71/Wenzhou/CHN/2014) were given the following accession numbers: KT008670 and KT008671, respectively. In this study, six nucleotide variations (GP151→TP151, GP199→AP199, GP261→TP261, AP328→CP328, GP422→AP422, and GP437→TP437) within the IRES region were observed simultaneously in the two virulent EV71 isolates when compared w ith the four common isolates (Figure 1A). RNA secondary structure predictions for IRES, which includes domain I to domain V, were carried out. The results showed that the four common isolates shared sim ilar RNA structures, whereas the two virulent isolates had a large circular stem loop structure between domain II and domain III. This structural difference between the common and virulent isolates may have an impact on RNA translation and neurovirulence (Figure 2). In order to learn more about the subtle RNA structure variations, we also performed RNA secondary structure predictions for two specific functional regions: SL II and FCE. The results indicated that the variation of GP114→CP114in virulentisolate11/EV71/Wenzhou/CHN/2014 contributed to the formation of an additional stem loop structure and that the mutation of GP151→TP151in both virulent isolates produced an enhanced stem loop structure (w ith two more nucleotides) (Figure 3A). Moreover, the variations of GP114→CP114and GP151→TP151resulted in higher m inimum free energy (-12.4 kcal/mol and -17 kcal/mol, respectively) and contributed to the formation of relatively stable secondary structures (Figure 3A). The predicted secondary structures of the FCE of the four common isolates were sim ilar, whereas the structures of two virulent isolates were variable (Figure 3B). Wen et al. demonstrated that the 5'UTR RNA secondary structures of virulent EV71 isolates differed from those of common isolates[8]. However, they did not perform a detailed analysis of additional subtle variations in key regions w ithin the 5'UTR. In this study, instead of repeating previous studies on predictions of the RNA secondary structure of the whole 5'UTR, we selected the IRES for RNA secondary structures analysis and specifically studied two special regions (SL II and FCE) that were previously reported to be associated w ith neurovirulence and viral mRNA translational efficiency, respectively. In agreement w ith previous studies on the structure of the whole 5'UTR, the secondary structures of SL II and FCE were sim ilar among the four common isolates but differed from those of the two virulent isolates. Furthermore, sequence alignment of the 3'UTR (84 nts) suggested that the nucleotides in the 3'UTR were relatively conserved and that only four nucleotide positions were variable. The variation of TP10→CP10or AP10was found in two virulent isolates when compared w ith four common isolates (Figure 1B).
Next, 76 EV71 isolates along w ith the previous six EV71 isolates were divided into ‘Virulent' and‘Common' groups (Table S2 in the website of BES, www.besjournal.com). The nucleotides/am ino acids that exhibited large differences between the‘Virulent' and ‘Common' groups were picked and subjected to Chi-square tests. For example, in the 488th position of the 5'UTR, there were two ‘T's in the ‘Common' group but seven ‘T's in the ‘Virulent' group. The results of the Chi-square tests showed that this difference was statistically significant, and it was concluded that this difference was associated with the neurovirulence of EV71 (χ2=4.204, P= 0.04, α=0.05) (Table 1). As a result, two alternativenucleotides (TP488and CP577) in the 5'UTR, two alternative nucleotides (CP10and AP47) in the 3'UTR, and three alternative am ino acids (GlnP22in VP1; AsnP57in 2A; IleP56in 3C) were considered to be relevant to neurovirulence. To expound the impact of these neurovirulence-related am ino acids on the structures of the corresponding proteins, 3D structurepredictionsoftwareforproteins (https://sw issmodel.expasy.org/interactive/xTem Lr/ models/) was used to predict the structures of VP1, 2A, and 3C. We wanted to determ ine whether, like the predicted RNA structures of the 5'UTR, the 3D structures of these isolates differed between the common and virulent isolates. Unfortunately, the 3D structures were not noticeably different between the common and virulent isolates. 2A had almost the same structure for all six isolates, whereas both VP1 and 3C were variable among the six isolates. The on-line software (ProMod3 Version 1.0.0.) used for the analysis is the most commonly used software for the 3D structure analysis of proteins, and all sequence identities between our am ino acid sequences and the models built on-line were above 90%, w ith several sequence identities reaching 100% (seq identities >50% were conclusive; for more information, please refer to the SWISS-MODEL Homology Modeling Reports in the supplementary materials) (Figure S2 in the website of BES, www.besjournal.com). As a result, these 3D results are probably attributed to the following factors: (1) the lim ited number of samples on 3D structure predictions; (2) a single am ino acid change in a pivotal site may contribute more to EV71 neurovirulencethanthestructureofthe corresponding protein; and (3) a crucial am ino acid may determ ine EV71 vitality when binding to host cells. In summary, further studies are needed to determ ine how these variations affect the neurovirulence of EV71.
Researchers have explored the association of certainnucleotides/am inoacidswiththe neurovirulence of EV71. For example, Liu et al. reported four neurovirulence-related am ino acid substitutions,HisP22→GlnP22,GluP98→LysP98, GluP145→GlyP145, and AlaP289→ThrP289, in VP1[9]. ValP263→IleP263in 3D of virulent isolates was demonstrated to be a unique am ino acid variation that distinguished the virulent isolates from common isolates[10]. Among the seven nucleotides/am ino acids identified in the present study, GlnP22in VP1 has already been described[9]. However, neither the four common isolates nor the two virulent isolates had this neurovirulence-associated am ino acid. Of the six new ly identified neurovirulence-associated nucleotides/am ino acids, five variations were unique to the two virulent isolates. It should be emphasized that we searched and analyzed the whole genome of almost all virulent EV71 isolates (C4 subgenotype) that had been described in either published papers or the GenBank database. Moreover, all analyzed EV71 isolates were carefully chosen. All common isolates were from patients w ith herpes only on the skin and mucosae and w ithout any neurological symptoms. Virulent isolates were isolated from patients w ith not only herpes but also at least one of the follow ing severe neurological complications: brain stem encephalitis, aseptic encephalitis, meningitis,poliomyelitis-likeparalysis,and neurogenic pulmonary edema. Most importantly, the number of EV71 isolates (either common or virulent isolates) utilized in the present study was higher than in any previous EV71 genome sequence analysis of C4. Therefore, all of these factor make our findings more convincing.
Taken together, we demonstrated the genetic characteristics of six native EV71 isolates and identified six novel nucleotide/am ino acid sites (TP488and CP577in the 5'UTR; AsnP57in 2A; IleP56in 3C; CP10and AP47in the 3'UTR) that were concluded to be associated with the neurovirulence of EV71. Furthermore, the RNA structure predictions of IRES, SL II, and FCE suggested that there may be a relationship between EV71 neurovirulence and the RNA secondary structures of IRES, SL II, and FCE in the 5'UTR. The results of the present study add to the know ledge on EV71 origin, evolution, and neurovirulence and will aid in developing safe and effective vaccines. In the future, further studies should be conducted to determ ine whether and how these am ino acid/nucleotide variations, as well as changes to the corresponding RNA and protein structures, affect neurovirulence using cell culture systems, animal models, and molecular biology techniques.
We would like to thank the participants for their cooperation in this study.
?These authors contributed equally to this work.
#Correspondence should be addressed to WEN Jin Sheng, E-mail: w js78@wmu.edu.cn, Tel: 86-577-86689910; DUAN Zhi Liang, E-mail: dzl3032@163.com, Tel: 86-577-88816276.
Accepted: October 1, 2016
REFERENCES
1. Skinner MA, Racaniello VR, Dunn G, et al. New model for the secondary structure of the 5' non-coding RNA of poliovirus is supported by biochem ical and genetic data that also show that RNA secondary structure is important in neurovirulence. J Mol Biol, 1989; 207, 379-92.
2. Bedard KM, Sem ler BL. Regulation of picornavirus gene expression. M icrobes Infect, 2004; 6, 702-13.
3. Pelletier J, Kaplan G, Racaniello VR, et al. Cap-independent translation of poliovirus mRNA is conferred by sequence elements w ithin the 5' noncoding region. Mol Cell Biol, 1988; 8, 1103-12.
4. Mou J, Dawes M, Li Y, et al. Severe hand, foot and mouth disease in Shenzhen, South China: What matters most? Epidem iol Infect, 2014; 142, 776-88.
5. Yeh MT, Wang SW, Yu CK, et al. A single nucleotide in stem loop II of 5'-untranslated region contributes to virulence of enterovirus 71 in m ice. PLoS ONE, 2011; 6, e27082.
6. Kataoka C, Suzuki T, Kotani O, et al. The role of VP1 am ino acid residue 145 of enterovirus 71 in viral fitness and pathogenesis in a cynomolgus monkey model. PLoS Pathog, 2015; 11, e1005033.
7. Wang Y, Zhu QQ, Zeng M, et al. Complete genome sequence of a human enterovirus 71 strain isolated from a fatal case in Shanghai, China, in 2012. Genome Announc, 2014; 2. pii: e00457-14.
8. Wen HL, Si LY, Yuan XJ, et al. Complete genome sequencing and analysis of six enterovirus 71 strains w ith different clinical phenotypes. Virol J, 2013; 10, 115.
9. Liu YJ, Fu C, Wu SY, et al. A novel finding for enterovirus virulence from the capsid protein VP1 of EV71 circulating in mainland China. Virus Genes, 2014; 48, 260-72.
10.Chang GH, Lin L, Luo YJ, et al. Sequence analysis of six enterovirus 71 strains w ith different virulence in humans. Virus Res, 2010; 151, 66-73.
d
*oi: 10.3967/bes2016.103 This study was funded by Natural Science Foundation of Zhejiang (LQ14C010006); National Natural Science Foundation of China (81501363); and Planned Science and Technology Project of Zhejiang (2014C33261).
1. Institute of Arboviruses, School of Basic Medical Sciences, Wenzhou Medical University, Wenzhou 325000, Zhejiang, China; 2. Department of Clinical Laboratory, The Second Affiliated Hospital & Yuying Children's Hospital of Wenzhou Medical University, Wenzhou 325000, Zhejiang, China; 3. Department of Liver, The Second Hospital of Ningbo, Ningbo 315010, Zhejiang, China
Biographical notes of the s: JIA Qing Jun, male, born in 1990, Master's degree candidate, majoring in pathogen biology; CHEN Xin Yu, male, born in 1976, associate chief technician, Master's degree candidate, majoring in biomedical engineering.
May 13, 2016;
Biomedical and Environmental Sciences2016年10期