- Research article
- Open Access
- Open Peer Review
The variations of VP1 protein might be associated with nervous system symptoms caused by enterovirus 71 infection
BMC Infectious Diseases volume 14, Article number: 243 (2014)
The VP1 protein of enterovirus 71 (EV71) is an important immunodominant protein which is responsible for host-receptor binding. Nevertheless, the relationship between VP1 and neurovirulence is still poorly understood. In this study, we investigated the relationship between mutation of VP1 and neurovirulent phenotype of EV71 infection.
One hundred and eighty-seven strains from Genbank were included, with a clear clinical background. They were divided into two groups, one with nervous system symptoms and one with no nervous system symptoms. After alignment, the significance of amino acid variation was determined by using the χ2 test and a phylogenetic tree was constructed with MEGA software (version 5.1).
We showed no significant difference in neurovirulence between genotype B and C. Interestingly, we found that variations of E145G/Q, E164D/K and T292N/K were associated with nervous system infection in genotype B. In the case of genotype C, the N31D mutation increased the risk for nervous complications, whereas I262V mutation decreased the risk of nervous complications. We used a 3D model of VP1 to demonstrate the potential molecular basis for EV71 nervous system tropism.
Distinct variations are shown to be associated with neurovirulent phenotype in the different genotype. Detection of variation in genotypes and subtypes may be important for the prediction of clinical outcomes.
Enterovirus 71 (EV71) is a member of human enterovirus species, which is one of the most common causative agents for herpangina and hand, foot and mouth disease (HFMD) [1, 2]. EV71 infection can also result in acute cardiogenic shock, pulmonary edema, and severe neurological diseases such as encephalitis and meningitis, the major causes for mortality among infants and young children [3, 4], as well as in adults . Recently, outbreaks of EV71 infection have been frequently reported throughout the world, including China , South Korea , Taiwan , Malaysia  and Singapore .
Although molecular epidemiological studies have been extensively conducted by different groups [6, 9, 11–16], little has been done to reveal the correlation of specific viral genotypes or residue variations and nervous system tropism[17, 18]. The mechanism of the neurological complications of enterovirus infection is poorly understood. Neurological virulence is one of the most severe complications responsible for death . It has been reported that patients with central nervous system symptoms accounted for 52.9% (92/174) of EV71 infection . Cordey et al. have demonstrated that an L97R alteration of VP1 protein enhances the neuronal tropism of EV71. Some adaptations of VP1, 5′ NCR, protease 2A determine viral virulence . These results suggested that the sequence variations may contribute to neural infection and neurological complications.
VP1 is one of structural proteins of EV71 virus, which is involved in forming the pentameric icosahedral structure and important for virus to bind receptors of host cells. In this study, we analyzed the correlation between VP1 variations and the neurological virulence caused by EV71 infection.
Nucleotide sequences screening
We took the advantage of BLASTN tool to search VP1 gene of EV71, and cross-validated in the database of EXPASY. The sequences were doubly checked by going through 324 papers one by one. The complete VP1 sequences of 187 EV71 strains (Additional file 1: Table S1) were downloaded from Genbank as of Sep 30, 2012. Seventy-six strains isolated from patients with nervous system symptoms (e.g., encephalitis, aseptic meningitis, and poliomyelitis-like paralysis), were referred to as nervous system EV71 strains (NS-EV71), and the other 111 strains from only HFMD and mild HFMD were considered as non-nervous system EV71 strains (nNS-EV71).
Alignment and phylogenetic tree construction
Sequence alignment and phylogenetic tree was constructed with MEGA software (version 5.1) using the neighbor-joining method. The robustness of the constructed phylogeny was estimated by bootstrap analysis with 1,000 pseudo-replicate data sets.
3D reconstruction of VP1 protein
The crystal structure of the VP1-VP4 protein complex of EV 71 (PDB ID: 3VBS) was obtained from NCBI database (http://www.ncbi.nlm.nih.gov/Structure/mmdb/mmdbsrv.cgi?uid=97658). The 3D structure of VP1 protein was modeled by using free software Cn3D (version 4.3.1, http://www.ncbi.nlm.nih.gov/Structure/CN3D/cn3d.shtml).
The significance of the amino acid variations associated with neurovirulent phenotype was detected using the χ2 test by SPSS 13.0 software. A p value less than 0.05 was considered as statistically significance.
Phylogenetic analysis of VP1
One hundred and eighty-seven strains were carefully screened from Genbank and Pubmed, and divided into two groups. NS-EV71 group contained 76 strains, and nNS-EV71 had 111 strains. A phylogenetic tree was constructed by using 891 nucleotides of VP1 sequences of 187 EV71 strains (Figure 1 and Additional file 1: Figure S1A/B for genotype B, C). These strains were divided into three groups: A, B (B1–B5) and C (C1–C4) genotypes. The strain AF135899 was included in B genotype, but did not belong to B1-B5. Except subtypes B1 and B2, there were no concentrated clusters in NS-EV71 strains in other subtypes of genotype B and C. More importantly, some EV71 subtypes are geographical and time specific, such as, C4 in China and Taiwan, and B4-5 mainly in Malaysia (Additional file 1: Figure S1A/B).
Genotype/subtype analysis in EV71
The genotype or subtype numbers of NS-EV71 and nNS-EV71 cases and their distribution worldwide are listed in Additional file 1: Table S2. Based on SPSS analysis, genotype B or C was not associated with nervous system symptoms (χ2 = 2.021, p = 0.103). However, when further analyzed subtypes, we found that the subtypes B1–B3 were positively correlated with nervous system complications (χ2 = 16.565, p = 0.002). In the case of subtypes in genotype C, no obvious correlation between subtypes and neurovirulence was detected (χ2 = 5.396, p = 0.145).
Relationship between amino acid variations and nervous system symptoms
The amino acid sequences were aligned to find the correlation between amino acid variations and nervous system symptoms. Consensus changes were observed to correlate with specific genotype at position 43, 58, 184 and 240 (Table 1). In genotype B, all subtypes had the consensus amino acids 43E, 58 T, 184 T and 240S; and seven other residues displayed high frequency of substitution in both NS-EV71 and nNS-EV71 cases (Table 2). The amino acid changes of E145G/Q (OR, 3.5; 95% CI, 1.228–9.978), E164D/K (OR, 6.652; 95% CI, 1.315–33.647), and T292N/K (OR, 8; 95% CI, 2.294–27.9) were significantly correlated with nervous system symptoms. The other four amino residues were not significantly associated with neurovirulence.
Genotype C was distinct from genotype B for the consensus sequence with 43 K, 58A, 164D, 184S and 240 T. Three residue substitutions (N31D, V170A and V262I) were associated with nervous system symptoms in C genotype (Table 2). The OR and 95% CI of N31D were 3.044 and 1.002–9.252, indicating that N31D variant was significantly associated with EV71 neurovirulence. More interestingly, we noted that C4 subtype had a consensus residue 262I. To our surprise, the risk of V262I substitution in subtype C1–3 was negatively associated with neurovirulence (χ2 = 6.262, p = 0.012; OR, 0.239; 95% CI, 0.075–0.756). As the occurrence rate of A170V was too low, the role of A170V in clinical presentation needs to be verified in larger samples. There was no amino acid position linked with nervous system symptoms in the C4 subtype (data not shown).
Locations of specific amino acid variation of VP1 protein
Aligned with the structure of VP1 (PDB 3VBS), the consensus sequences of genotypes B and C, 184 T/S and 240S/T positioned at the loop of FG and HI regions, protruding the external surface of virus body and the neighbor canyon that is important for the receptor binding; while 43E/K and 58 T/A were located at the N terminus of VP1, buried in the inner of virus body (Figure 2A).
The residue variations of E145G/Q and E164D/K in genotype B were sited at the DE and EF loop, respectively. They located at the protrusion of the canyon surface. However, T292N placed at the C terminus, which was exposed on the virion surface. The alterations of N31D and I262V in genotype C were situated at the N terminus and β-strand I, respectively (Figure 2B).
EV71 is a picornavirus that causes mild HFMD in most cases. However, it is also frequently reported to cause severe nervous system complications. In this study, we found that no significant relationship was exhibited between neurovirulence and genotype or subtype, except subtypes B1 and B2, which is consistent with the previous study , although Ooi MH observed that children with B4 genotype were less likely to have CNS infection than those with other genotypes .
It has been reported that D164E variation was associated with severe cases of EV71 infection . In the severe cases of EV71 infection patients, the nervous system infection was most frequent. Nervous system tropism for the neural infection might be responsible for the neurovirulence of EV71. The attachment to target cell receptors mediated by VP1 protein was the first and key step for EV71 infection . The residue variations of VP1 protein, which binds receptors on the surface of host cells (e.g., neural cells), may change the specificity of VP1-receptor interaction and therefore play an important role in neurological complications.
So far, four receptors, including P-selectin glycoprotein ligand-1(PSGL-1) , scavenger receptor class B member 2 (SCARB2) [25, 26], sialylated glycans  and heparan sulfate glycosaminoglycan , were identified. In our study, we analyzed the complete VP1 protein sequences of 187 strains with detailed clinic information screened from the Genbank database. The results showed that the consensus sequences of genotype B was distinct from genotype C, which consisted of 43E/K, 58 T/A, 184 T/S and 240S/T (Table 1). While 3D structure presented that residue 184 and 240 lay at the protrusion of the pocket that was responsible for recognition and binding to receptors (Figure 2A) [29, 30], as indicated that different genotypes of EV71 may bind different receptors. The variation of E145G/Q and E164D/K in genotype B was correlated with neurovirulent phenotype as reported by a previous study , but not in genotype C; meanwhile, 164D was one of the consensus sequences in genotype C. However, Zaini et al. showed that Q145E mutant of a C4 strain generated a mouse-virulent phenotype . The substitution of T292N/K has a higher risk for nervous system complications. At this site, N292 is the consensus sequence for subtype of B1 and B2, which did not present in other genotypes. This might be the sequence basis of B1 and B2 subtypes for the elevated risk for neurovirulence than those with other subtypes or genotype. It may increase VP1-receptor binding affinity by forming an ionic bond. The residue 292 might be also responsible for receptor binding, similar to K149 of VP2 protein which interacts with one of the receptors for PSGL-1 .
In the case of genotype C strains, the variation N31D clearly increased the risk of nervous system infection, while I262 V decreased the risk. The N31D substitution took place in the inner of the virion structure. The replacement of I262 V mainly occurred in subtype C1–C3 and was situated at the interaction center of VP1, VP2 and VP3, and at the bottom of the canyon. This conversion may influence the hydrophobic interactions among related proteins, resulting in the loose of canyon structure. No significant change was found to be associated with neurovirulent phenotype in the C4 genotype, which mostly causes outbreaks in China, Taiwan and Japan. Cordey S et al. has identified that L97R substitution (Genbank NO.EU414333 which was included in subtype C1), could enhance the neural cell tropism in the EV71 infection . In all 187 strains, there were only two clones with L97R substitution in the database, the other was Genbank NO.AF009539 which belonged into the subtype B2. Therefore, the L97R variation might not be frequent during the human EV71 infection.
Due to the limited isolates (about 180) and their geographical distribution of genotypes in this study, our findings may not represent the full picture of VP1 variations that are associated with neurological disorders. For example, we only identified 4 B1 and 8 B2 isolates that were associated with neurovirulent in Austria and US, respectively (Table 2). Other variations, which might be associated with neurovirulent phenotype, could be missed because they had not be isolated and sequenced. Other factors, for example, the fitness of virus in host would also play crucial role in the disease process. Tee et al. showed that EV71 subgenogroups have been circulating cryptically in human populations for years before causing large-scale HFMD outbreaks . This would make things more complicated. On the other hand, potential biases may also arise in this study. The residue variants identified in this study needed to be further verified by large sample size with more even geographic distributions, and combined with biological functional studies in both neuron cell and animal models.
In conclusion, some variations of EV71 VP1 protein associated with neurovirulent phenotype were mapped by molecular epidemiology in this study. Detection of these variations coupled with the analysis of the genotype or subtype might contribute to the prediction of clinical presentation. Further experimental data will be needed to show the significance of the results demonstrated in the present study.
Hand, foot and mouth disease
Nervous system EV71 strains
Non-nervous system EV71 strains
Huang HI, Weng KF, Shih SR: Viral and host factors that contribute to pathogenicity of enterovirus 71. Future Microbiol. 2012, 7 (4): 467-479. 10.2217/fmb.12.22.
Lan YC, Lin TH, Tsai JD, Yang YC, Peng CT, Shih MC, Lin YJ, Lin CW: Molecular epidemiology of the 2005 enterovirus 71 outbreak in central Taiwan. Scand J Infect Dis. 2011, 43 (5): 354-359. 10.3109/00365548.2010.545995.
Bek EJ, McMinn PC: The pathogenesis and prevention of encephalitis due to human enterovirus 71. Curr Infect Dis Rep. 2012, 14 (4): 397-407. 10.1007/s11908-012-0267-3.
Wang SM, Lei HY, Yu CK, Wang JR, Su IJ, Liu CC: Acute chemokine response in the blood and cerebrospinal fluid of children with enterovirus 71-associated brainstem encephalitis. J Infect Dis. 2008, 198 (7): 1002-1006. 10.1086/591462.
Hamaguchi T, Fujisawa H, Sakai K, Okino S, Kurosaki N, Nishimura Y, Shimizu H, Yamada M: Acute encephalitis caused by intrafamilial transmission of enterovirus 71 in adult. Emerg Infect Dis. 2008, 14 (5): 828-830. 10.3201/eid1405.071121.
Fan X, Jiang J, Liu Y, Huang X, Wang P, Liu L, Wang J, Chen W, Wu W, Xu B: Detection of human enterovirus 71 and Coxsackievirus A16 in an outbreak of hand, foot, and mouth disease in Henan Province, China in 2009. Virus Genes. 2013, 46 (1): 1-9. 10.1007/s11262-012-0814-x.
Ryu WS, Kang B, Hong J, Hwang S, Kim J, Cheon DS: Clinical and etiological characteristics of enterovirus 71-related diseases during a recent 2-year period in Korea. J Clin Microbiol. 2010, 48 (7): 2490-2494. 10.1128/JCM.02369-09.
Huang YP, Lin TL, Hsu LC, Chen YJ, Tseng YH, Hsu CC, Fan WB, Yang JY, Chang FY, Wu HS: Genetic diversity and C2-like subgenogroup strains of enterovirus 71, Taiwan, 2008. Virol J. 2010, 7: 277-10.1186/1743-422X-7-277.
Chua KB, Kasri AR: Hand foot and mouth disease due to enterovirus 71 in Malaysia. Virol Sin. 2011, 26 (4): 221-228. 10.1007/s12250-011-3195-8.
Wu Y, Yeo A, Phoon MC, Tan EL, Poh CL, Quak SH, Chow VT: The largest outbreak of hand; foot and mouth disease in Singapore in 2008: the role of enterovirus 71 and coxsackievirus A strains. Int J Infect Dis. 2010, 14 (12): e1076-e1081. 10.1016/j.ijid.2010.07.006.
Honkanen H, Oikarinen S, Pakkanen O, Ruokoranta T, Pulkki MM, Laitinen OH, Tauriainen S, Korpela S, Lappalainen M, Vuorinen T, Haapala AM, Veijola R, Simell O, Ilonen J, Knip M, Hyöty H: Human enterovirus 71 strains in the background population and in hospital patients in Finland. J Clin Virol. 2013, 56 (4): 348-353. 10.1016/j.jcv.2012.11.018.
Zhang Y, Tan X, Cui A, Mao N, Xu S, Zhu Z, Zhou J, Shi J, Zhao Y, Wang X, Huang X, Zhu S, Zhang Y, Tang W, Ling H, Xu W: Complete genome analysis of the c4 subgenotype strains of enterovirus 71: predominant recombination c4 viruses persistently circulating in china for 14 years. PLoS One. 2013, 8 (2): e56341-10.1371/journal.pone.0056341.
Lee MS, Tseng FC, Wang JR, Chi CY, Chong P, Su IJ: Challenges to licensure of enterovirus 71 vaccines. PLoS Negl Trop Dis. 2012, 6 (8): e1737-10.1371/journal.pntd.0001737.
Chan YF, Wee KL, Chiam CW, Khor CS, Chan SY, Amalina WM, Sam IC: Comparative genetic analysis of VP4, VP1 and 3D gene regions of enterovirus 71 and coxsackievirus A16 circulating in Malaysia between 1997–2008. Trop Biomed. 2012, 29 (3): 451-466.
Mirand A, Schuffenecker I, Henquell C, Billaud G, Jugie G, Falcon D, Mahul A, Archimbaud C, Terletskaia-Ladwig E, Diedrich S, Huemer HP, Enders M, Lina B, Peigue-Lafeuille H, Bailly JL: Phylogenetic evidence for a recent spread of two populations of human enterovirus 71 in European countries. J Gen Virol. 2010, 91 (Pt 9): 2263-2277.
Chan YF, Sam IC, AbuBakar S: Phylogenetic designation of enterovirus 71 genotypes and subgenotypes using complete genome sequences. Infect Genet Evol. 2010, 10 (3): 404-412. 10.1016/j.meegid.2009.05.010.
Cordey S, Petty TJ, Schibler M, Martinez Y, Gerlach D, van Belle S, Turin L, Zdobnov E, Kaiser L, Tapparel C: Identification of site-specific adaptations conferring increased neural cell tropism during human enterovirus 71 infection. PLoS Pathog. 2012, 8 (7): e1002826-10.1371/journal.ppat.1002826.
Li R, Zou Q, Chen L, Zhang H, Wang Y: Molecular analysis of virulent determinants of enterovirus 71. PLoS One. 2011, 6 (10): e26237-10.1371/journal.pone.0026237.
Perez-Velez CM, Anderson MS, Robinson CC, McFarland EJ, Nix WA, Pallansch MA, Oberste MS, Glode MP: Outbreak of neurologic enterovirus type 71 disease: a diagnostic challenge. Clin Infect Dis. 2007, 45 (8): 950-957. 10.1086/521895.
Ryu WS, Kang B, Hong J, Hwang S, Kim A, Kim J, Cheon DS: Enterovirus 71 infection with central nervous system involvement, South Korea. Emerg Infect Dis. 2010, 16 (11): 1764-1766. 10.3201/eid1611.100104.
Herrero LJ, Lee CS, Hurrelbrink RJ, Chua BH, Chua KB, McMinn PC: Molecular epidemiology of enterovirus 71 in peninsular Malaysia, 1997–2000. Arch Virol. 2003, 148 (7): 1369-1385. 10.1007/s00705-003-0100-2.
Ooi MH, Wong SC, Podin Y, Akin W, Del SS, Mohan A, Chieng CH, Perera D, Clear D, Wong D, Blake E, Cardosa J, Solomon T: Human enterovirus 71 disease in Sarawak, Malaysia: a prospective clinical, virological, and molecular epidemiological study. Clin Infect Dis. 2007, 44 (5): 646-656. 10.1086/511073.
Lu J, He YQ, Yi LN, Zan H, Kung HF, He ML: Viral kinetics of enterovirus 71 in human abdomyosarcoma cells. World J Gastroenterol. 2011, 17 (36): 4135-4142. 10.3748/wjg.v17.i36.4135.
Miyamura K, Nishimura Y, Abo M, Wakita T, Shimizu H: Adaptive mutations in the genomes of enterovirus 71 strains following infection of mouse cells expressing human P-selectin glycoprotein ligand-1. J Gen Virol. 2011, 92 (Pt 2): 287-291.
Yamayoshi S, Yamashita Y, Li J, Hanagata N, Minowa T, Takemura T, Koike S: Scavenger receptor B2 is a cellular receptor for enterovirus 71. Nat Med. 2009, 15 (7): 798-801. 10.1038/nm.1992.
Su PY, Liu YT, Chang HY, Huang SW, Wang YF, Yu CK, Wang JR, Chang CF: Cell surface sialylation affects binding of enterovirus 71 to rhabdomyosarcoma and neuroblastoma cells. BMC Microbiol. 2012, 12: 162-10.1186/1471-2180-12-162.
Yang B, Chuang H, Yang KD: Sialylated glycans as receptor and inhibitor of enterovirus 71 infection to DLD-1 intestinal cells. Virol J. 2009, 6: 141-10.1186/1743-422X-6-141.
Tan CW, Poh CL, Sam IC, Chan YF: Enterovirus 71 uses cell surface heparan sulfate glycosaminoglycan as an attachment receptor. J Virol. 2013, 87 (1): 611-620. 10.1128/JVI.02226-12.
Wang X, Zhu JP, Zhang Q, Xu ZG, Zhang F, Zhao ZH, Zheng WZ, Zheng LS: Detection of enterovirus 71 using reverse transcription loop-mediated isothermal amplification (RT-LAMP). J Virol Methods. 2012, 179 (2): 330-334. 10.1016/j.jviromet.2011.11.019.
Plevka P, Perera R, Cardosa J, Kuhn RJ, Rossmann MG: Crystal structure of human enterovirus 71. Science. 2012, 336 (6086): 1274-10.1126/science.1218713.
Zaini Z, McMinn P: A single mutation in capsid protein VP1 (Q145E) of a genogroup C4 strain of human enterovirus 71 generates a mouse-virulent phenotype. J Gen Virol. 2012, 93 (Pt 9): 1935-1940.
Tee KK, Lam TT, Chan YF, Bible JM, Kamarulzaman A, Tong CY, Takebe Y, Pybus OG: Evolutionary genetics of human enterovirus 71: origin, population dynamics, natural selection, and seasonal periodicity of the VP1 gene. J Virol. 2010, 84 (7): 3339-3350. 10.1128/JVI.01019-09.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/243/prepub
This study was partly supported by National Basic Science Research Program project (No. 2013CB911302) from the Ministry of Science and Technology of China, RGC-GRF (No.464512), the Natural Science Foundation of China (No. 30600733 and No.81272483), Guangdong Provincial Science and Technology (NO. 2010A040302003 and No. 2011B031800163), and Shenzhen Strategic Innovative Development Program (JCYJ20120615085419420) and JCYJ20130401172046448.
The authors declare that they have no competing interests.
BZ and LL drafted the manuscript. XW and CW analyzed the data. KH and LZ download sequences from Genbank database. MLH refined the paper. WZ carried out the design. All authors read and approved the final manuscript.
Bao Zhang, Xianbo Wu, Keyong Huang, Ling Li contributed equally to this work.