Skip to main content

Molecular characteristics of the VP1 region of enterovirus 71 strains in China



Enterovirus 71 (EV71) is the most commonly implicated causative agent of severe outbreaks of paediatric hand, foot, and mouth disease (HFMD).VP1 protein, a capsid protein of EV71, is responsible for the genotype of the virus and is essential for vaccine development and effectiveness. However, the genotypes of EV71 isolates in China are still not completely clear.


The VP1 gene sequences of 3712 EV71 virus strains from China, excluding repetitive sequences and 30 known EV71 genotypes as reference strains, between 1986 and 2019 were obtained from GenBank. Phylogenetic tree, amino acid homology, genetic variation and genotype analyses of the EV71VP1 protein were performed with MEGA 6.0 software.


The amino acid identity was found to be 88.33%–100% among the 3712 EV71 strains, 93.47%–100% compared with vaccine strain H07, and 93.04%–100% compared with vaccine strains FY7VP5 or FY-23 K-B. Since 2000, the prevalent strains of EV71 were mainly of the C4 genotype. Among these, the C4a subgenotype was predominant, followed by the C4b subgenotype; other subgenotypes appeared sporadically between 2005 and 2018 in mainland China. The B4 genotype was the main genotype in Taiwan, and the epidemic strains were constantly changing. Some amino acid variations in VP1 of EV71 occurred with high frequencies, including A289T (20.99%), H22Q (16.49%), A293S (15.95%), S283T (15.11%), V249I (7.76%), N31D (7.25%), and E98K (6.65%).


The C4 genotype of EV71 in China matches the vaccine and should effectively control EV71. However, the efficacy of the vaccine is partially affected by the continuous change in epidemic strains in Taiwan. These results suggest that the genetic characteristics of the EV71-VP1 region should be continuously monitored, which is critical for epidemic control and vaccine design to prevent EV71 infection in children.


Enterovirus 71 (EV71) is a common pathogen of hand, foot, and mouth disease (HFMD) in children, and it is also the most important risk factor for severe cases and deaths [1,2,3,4].Some children with HFMD can have severe neurological effects, such as aseptic meningitis, encephalitis and acute delayed paralysis, develop serious brainstem encephalitis, neurogenic pulmonary oedema and even die [5, 6]. Children with severe neurological diseases who survive often have irreversible sequelae, seriously threatening their health [7, 8].

EV71 was first discovered in 1969; the virus was distributed mainly in the Americas, Europe and other countries in a sporadic form, with outbreaks in some European countries in the 1970s and 1980s [9, 10]. After 1997, EV71 began to emerge and spread in Asia, and the Asia–Pacific region is the most prevalent area for EV71. Indeed, there are reports of EV71 outbreaks in China, Singapore, and Malaysia, among others [11, 12]. Diseases caused by EV71 infection have been widely prevalent in China since 2007 [13]. For example, HFMD pandemics in Linyi city, Shandong Province, and Guangdong, Anhui Province, in 2008 resulted in tens of thousands of childhood infections and death among dozens of children [14, 15].

EV71 belongs to the Enterovirus genus of the RNA virus family. EV71 can be clustered into three genotypes according to nucleotide differences in the VP1 region, including genotype A (BrCr) with only one member, B and C. In contrast, genotypes B and C are divided into five subgenotypes, B1–B5 and C1–C5, respectively, and genotype C4 is further subdivided into C4a and C4b [16,17,18]. In China, genotypeC4 is the main epidemic strain; C4b was the predominant epidemic genotype from 1998 to 2004 and the C4a subgenotype after 2004 in mainland China, though Taiwan strains continue to circulate, including C2, B4 and B5 genotypes [19,20,21,22,23,24,25]. Phase III clinical trials of vaccines from three companies in China have been completed, and the genotypes of their vaccine strains are all C4a subgenotypes [26]. The Vaccine Research and Development Center of National Institutes of Health in Taiwan has also developed an FI-EV71 vaccine based on the B4 subtype (EV71vac), which can cause a robust cross-neutralizing antibody reaction against different EV71 gene subtypes, such as B4, B1, B5 and C4a [26, 27].

Nonetheless, there are many reports on recombination between different genotypes of EV71 [28,29,30,31], suggesting that EV71 has high variability and recombination ability, which may lead to the production of new pathogenic strains. Therefore, genome monitoring of EV71 epidemic strains is of great significance for the prevention and control of EV71 epidemics and can guide the application of the EV71 vaccine to a certain extent. In this study, the VP1 sequences of all EV71 viruses registered in GenBank in China from 1996 to 2019 were collected, and the molecular characteristics of genes were analysed using bioinformatics software to provide a scientific basis for the prevention and control of HFMD epidemics.


Acquisition of the EV71VP1 gene sequence

The complete 891-bp VP1 gene sequence of EV71 strains isolated from children with HFMD from 1986 to December 31, 2019, in China with known collection dates and isolation regions, including mainland China, Hong Kong, Macao and Taiwan, was obtained from the GenBank public database at the National Center for Biotechnology Information (NCBI) PubMed website ( A total of 8340 EV71 strains from China were collected from GenBank, consisting of 6572 strains from mainland China, 26 from Hong Kong, none from Macao, and 1742 from Taiwan. If EV71 strains were isolated the same year and from the same region and had 100% nucleotide homology, only one strain was included. Therefore, a total of 3470 strains from mainland China, 2 from Hong Kong, and 1156 from Taiwan were removed. Ultimately, 3712 strains, consisting of 3102 from mainland China, 24 from Hong Kong, and 586 from Taiwan, were retained for the study. The time and geographical distribution of all EV71 isolates are shown in Figs. 1 and 2. The nucleotide sequence of the complete VP1 gene of all virus strains is 891 bp, and Chinese vaccine strains (H07, FY7VP5, and FY-23 K-B) and 27 known genotypes of EV71 strains were considered reference sequences, as presented in Table 1.

Fig. 1

The number of EV71 strains isolated in different years in China. The 3712 EV71 strains with complete VP1 gene sequences collected in China were obtained from the GenBank public database at the National Center for Biotechnology Information (NCBI) PubMed website ( between 1986 and December 31, 2019. However, no EV71 strains were available for 2019

Fig. 2

The number of EV71 strains isolated in different provinces of China. A total of 3712 EV71 strains were collected in China, including different provinces of mainland China (3102), Hong Kong (24), Macao (0) and Taiwan (586), between 1986 and December 31, 2019

Table 1 Data of reference strains of EV71 genotypes

Construction of an EV71VP1 phylogenetic tree and gene sequence analysis

The complete VP1 sequence of EV71 strains was compared by Molecular Evolutionary Genetics Analysis (MEGA) version 6.0 [32]. In brief, a phylogenetic tree of the EV71VP1 gene was constructed by the adjacency method (neighbour-joining, N-J) with 1000 bootstrap replications. Homology and variation of EV71 strain VP1 gene sequences were analysed by MEGA 6.0.


Phylogeny and homology analysis of the EV71 VP1 region

A phylogenetic tree of the VP1 amino acid sequence was constructed with 3712 EV71 isolates from China and 30 reference strains (Fig. 3). Amino acid identity among the 3712 EV71 strains is 88.33%–100%, 93.47%–100% compared with vaccine strain H07 and 93.04%–100% compared with vaccine strains FY7VP5 or FY-23 K-B. Among the EV71 strains, the C4 genotype accounted for most of the EV71 strains, and the C4a subgenotype was the most common. Moreover, B4, C4b and C2 were found to be important genotypes of EV71; other subgenotypes appeared sporadically.

Fig. 3

Phylogenetic tree of the entire VP1 gene of EV71. The 30 reference strains are indicated in different coloured shapes, as shown below: :HQ328793(H07) :JX025561(FY7VP5) :EU812515(FY-23 K-B) :HQ828086(C4a) :EU753365(C4a) :KU936124(C4a) :KU936125(C4a) :KU936130(C4a) :JQ742002(C4b) :JQ742001(C4b) :AY207625(C3) :AF376081(C3) :AF304457(C2) :AF376110 (C2) :AY125969(C1) :AY125976 (C1) :JN874558 (C5) :KU888174(C5) :U22521(A) :GU434678(A) :AB059814(B1) :AF135886 (B1) :AF009534(B2) :AF009540(B2) :AF376119(B3) :AF376073(B3) :AF376084(B4) :AJ586873 (B4) :AB177815 (B5) :AB177816 (B5)

Genotypic distribution of EV71 strains in different years in China

Among the 3712 strains, the first 6 strains isolated from China in 1986were all of the B3 genotype, after which different genotypes/subgenotypes were detected in China. From 1986 to 2007, there was a small epidemic peak of C2-genotype EV71 in 1998. Between 1999 and 2003, the B4 genotype was the predominant strain of EV71. Most EV71 strains between 2004 and 2005 were of the C3 genotype, and from 2006 to 2018, the C4a subgenotype was the most common epidemic strain. B4 and C4b genotypes also accounted for a certain proportion (Table 2 and Fig. 4).

Table 2 Number of EV71 genotypes for every year in China
Fig. 4

Subgenotype percentages of EV71 for every year in China. In this study, 12 genotypes were detected among 3712 strains. The B4 genotype was the predominant strain between 1999 and 2003 and C3 genotypes from 2004 to 2005; the C4a subgenotype was the most common epidemic strain from 2006 to 2018

Next, the EV71 strains isolated from China, including mainland China, Hong Kong, and Taiwan, were further analysed to explore their genotypic distribution in different years. The results indicated that from 2005 to 2018, the C4a subgenotype was absolutely the predominant strain of EV71, followed by the C4b subgenotype; other subgenotypes appeared sporadically in mainland China. Moreover, among the 24 strains isolated from Hong Kong, there were 10 C4a-subgenotype strains (1 strain in 2008, 8 strains in 2010, and 1 strain in 2012) and 14 C4b-subgenotype strains (1 strain in 2008 and 13 strains in 2010). The results for mainland China are shown in Table 3 and Fig. 5.

Table 3 Number of EV71 genotypes for every year in mainland China
Fig. 5

The genotypic distribution of EV71 every year in mainland China. The C4a subgenotype was absolutely the predominant strain of EV71 from 2005 to 2018, followed by theC4b subgenotype; other subgenotypes appeared sporadically in mainland China

However, the epidemic strains of EV71 in Taiwan from 1986 to 2018 are quite different from those in mainland China. Overall, epidemic strains are not predominantly represented by a single subgenotype but are constantly changing. The first 6 strains isolated in Taiwan in 1986 were all of the B3 genotype. In 1998, C2 was the predominant EV71 genotype in Taiwan. From 1999 to 2003, 2008 to 2009, and 2011 to 2012, the epidemic strains in Taiwan evolved mainly into B4-genotype EV71 strains. Additionally, C4a was an important genotype between 2010 and 2011 (Table 4 and Fig. 6).

Table 4 Number of EV71 genotypes for every year in Taiwan, China
Fig. 6

The subgenotype percentages of EV71 for every year in Taiwan. Epidemic strains of EV71 in Taiwan in recent years are constantly changing. The B3 genotype first emerged in 1986. From 1999 to 2003, 2008 to 2009, and 2011 to 2012, the B4 genotype was the predominant strain. Between 2010 and 2011, C4a was the main genotype

Analysis of amino acid variation of EV71VP1

Analysis of common amino acid variation sites in the EV71VP1 region showed that of the 3712 strains, residue 289 was the major common site, with a total variation rate of 23.55%; A289T was the most common variant, accounting for 20.99%. Moreover, 10.24% of strains from mainland China, including Hong Kong, exhibited this variation; 10.75% of strains from Taiwan also exhibited this variation. There was a total mutation rate of 21.85% at amino acid position 22 in the VP1 region. The frequency of the H22Q variation was detected to be 16.49%, of which mainland China, including Hong Kong, exhibited 5.50%; Taiwan exhibited 10.99%. Furthermore, 23% of the EV71 strains harboured an H22R mutation. The mutation rates of A293S, S283T, V249I, N31D and E98K were 15.95%, 15.11%, 7.76%, 7.76% and 6.65%, respectively. Among the common sites of amino acid variation in the EV71VP1 region, the most common was residue 293, with five types of amino acid variation, among which the A293S mutation was the most common (Table 5).

Table 5 Analysis of common amino acid variation sites in VP1 of EV71


It has been reported that since 2004, the major EV71 strains in mainland China basically belong to the C4 genotype [18, 33]. In this study, it was found that the epidemic strains of EV71 in mainland China have been dominated by the C4 genotype since 2000, mainly the C4b subgenotype from 2000 to 2004, and the proportion of the C4a subgenotype has increased significantly since 2005.From 2005 to 2018, the C4a subtype was dominant; the C4b subtype was the second most common, and other subtypes appeared sporadically; this pattern was similar to previous reports [20, 21, 30, 33]. From 2004 to 2018, most of the EV71 isolates in mainland China were genotype C, and most of them were subtype C4 (especially subtype C4a); other genotypes occurred sporadically (9 strains of genotype B in 2011, 4 strains of genotype A in 2008, 2 strains of genotype A in 2009, 1 strain of genotype An in 2011 and 1 strain of genotype A in 2012). The immune effectiveness of inactivated vaccines mostly depends on the antigenic correlation between epidemic strains and vaccine strains, which are often best for preventing infection of the same subtype virus but are inferior against different subtypes [30]. Recent studies have also shown that the EV71 vaccine (especially in children who receive 2 doses) can effectively prevent and control childhood EV71-associated HFMD but has no protective effect against coxsackievirus (CV) A6 (CVA-6) or CVA16, and there is no explanation for the effectiveness of other subtypes of EV71 (excluding C4a subtypes) [34]. These studies show that vaccine research and development for EV71 combined with CVA6 and CVA16 and other multivalent vaccines might better prevent EV71 infection.

Interestingly, this study found that the EV71 epidemic strains in Taiwan were mainly of the B4 genotype, which was different from those in mainland China; EV71 epidemic strains are also constantly changing, which is consistent with a previous report [26]. According to a human phase 1 clinical trial on adults in 2010, the FI-EV71 vaccine (EV71vac) based on the B4 genotype from Taiwan is safe and induces a high titre of neutralizing antibodies against EV71; it was also highly effective against B1, B5, and C4a strains. However, the titres of neutralizing antibodies against C4b and CVA16 were low in 20% of volunteers, and virus-neutralizing antibodies against the C2 genotype were not detected in 90% of vaccine recipients [26, 27]. These studies indicate that it is necessary to strengthen the monitoring of EV71 genotypes; new multivalent and effective vaccines that can cover local strains should be designed and applied according to the genotypes of the local predominant EV71 epidemic strains to ensure that the vaccine is more accurate in controlling HFMD epidemics.

Some studies have shown that the H22Q mutation in the VP1 protein of EV71 can lead to a decrease in the adsorption capacity of the C4 genotype to host cells [35,36,37]. The amino acid at position 22 of 78.15% of the 3712 strains isolated in China is H (histidine), which suggests that most of the viruses have strong adsorption capacity to host cells. Furthermore, H22Q was detected in 10.99% of all EV71 strains in Taiwan, significantly more prevalent than that in mainland China and Hong Kong (5.50%), suggesting that the adsorption capacity of some strains in Taiwan to host cells is weak in comparison with that of strains in mainland China.

Studies have shown that the A289T EV71VP1 variant is closely related to the occurrence of severe HFMD and that the neurological symptoms caused by EV71 infection are significantly increased when the amino acid at position 289 of VP1 is A (alanine); in contrast, there is low neurotoxicity when the amino acid is T (threonine) [36, 38]. In this study, 76.45% (2838 strains) of the virus strains were found to contain an A (alanine), suggesting that most of these EV71 viruses have high neurotoxicity. Moreover, 10.24% (380 strains) of the strains in mainland China (including Hong Kong) and 10.75% (399 strains) of those in Taiwan contain a T (threonine), suggesting low neurotoxicity. It remains to be further studied whether new mutations such as A289V (valine)/D (aspartic acid)/I (isoleucine) mutations will cause the emergence of severe HFMD.

EV71 can infect human lymphocytes by binding to its receptor molecule P-selectin glycoprotein ligand-1 (PSGL-1). When E (glutamic acid) at position 145 in VP1 is mutated to G (glycine) or Q (glutamine), the virus binds PSGL-1 more readily, whereas its PSGL-1-binding ability is weakened or lost when E is present [39]. In this study, the amino acid at position 145 in most strains was found to be E, with an E145G/Q mutation rate of 6.63%, suggesting that the emergence of this mutation may result in a virus that is more likely to infect human lymphocytes.

It has been reported that the E98K mutation may increase the hydrophobicity of VP1, making it easy for large compounds to enter and interfere with receptor binding, suggesting that E98K mutant viruses are sensitive to larger compounds [40]. Other studies have shown that E145G and N31D mutations are associated with increased virulence of EV71 and may increase the risk of neurological complications but that I262V mutations reduce the risk of neurological complications [41,42,43]. In this study, the E98K, E145G, N31D and I262V mutation rates were 6.65%, 3.13%, 7.25% and 2.77%, respectively. These findings indicate that these mutations may play an important role in the pathogenicity of mild and severe EV71-associated HFMD.

Humans are the only natural host and source of EV71. Indeed, EV71 cannot infect rodents, which is due mainly to the incompatibility between the virus and rodent cells, and the different expression of its scavenger receptor in humans and rodents [44,45,46]. However, some studies have found that simultaneous substitution of K98E, E145A and L169F in VP1 of EV71 can result in infection in mice [44]. Our study showed that among 3712 strains, the mutation frequencies of K98E, E145A and L169F were 93.24%, 0.30% and 0.03%, respectively; however, no strain with all three mutations was found. These findings indicate that humans are still the only host of EV71 in China; nevertheless, the existence of individual mutations does not rule out the emergence over time of strains that can infect other mammals. Therefore, it is important to closely monitor mutation of the key sites of the EV71VP1 protein.

EV71 is the most important pathogen causing severe HFMD in children, which can lead to irreversible sequelae or death, and it is a serious threat to their health [4, 7]. At present, there is no specific treatment for EV71 infection. The development and marketing of an inactivated EV71 vaccine in China is crucial for the prevention of HFMD caused by EV71 infection [26, 47,48,49]. Phase III clinical trials of the EV71 inactivated vaccine approved in China in 2015 have shown protective effectiveness against EV71-associated HFMD of more than 90% [47, 50,51,52]. However, according to molecular epidemiological studies of EV71, EV71 gene mutations occur frequently, leading to genetic diversity [28,29,30,31, 53]. These studies suggest that there is still a need for strengthening surveillance of EV71 genotypes and the development of new EV71 vaccines.

This study had a retrospective design, and there are some limitations. First, the EV71VP1 gene sequences from China analysed in this study were downloaded from the GenBank database but were not tested by us. Due to time constraints, only the VP1 region was analysed and studied. In future studies, we will conduct research on the complete genome sequence of EV71 in China. Second, the Chinese EV71VP1 strains registered in GenBank do not cover all provinces in the country, and the data for some years are missing; thus, some isolates of other genotypes may have been unavailable. Third, it is not clear whether some variations in the amino acid residues found in the study are related to the severity of disease or the route of transmission.


In summary, the prevalent strains of EV71 belong mainly to the C4 genotype. The C4a subgenotype was predominant, the C4b subgenotype was the second most prevalent, whereas other subgenotypes appeared sporadically in mainland China. The B4 genotype was the dominant genotype in Taiwan, and the epidemic strain is continuously changing. Moreover, variation in key positions of the EV71VP1 region is very important for the development of severe HFMD. Taken together, the findings indicate that the genetic characteristics of the EV71VP1 region should be continuously monitored, which is essential for the prevention and control of EV71-associated HFMD in children and EV71 vaccine design.

Availability of data and materials

The available data used and/or analysed during the current study are all included in the manuscript.


  1. 1.

    Yang F, Ren L, Xiong Z, et al. Enterovirus 71 outbreak in the People’s Republic of China in 2008. J ClinMicrobiol. 2009;47(7):2351–2.

    Google Scholar 

  2. 2.

    Wang Y, Feng Z, Yang Y, et al. Hand, foot, and mouth disease in China: patterns of spread and transmissibility. Epidemiology. 2011;22(6):781–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Xing W, Liao Q, Viboud C, et al. Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Solomon T, Lewthwaite P, Perera D, et al. Virology, epidemiology, pathogenesis, and control of enterovirus 71. Lancet Infect Dis. 2010;10(11):778–90.

    PubMed  Google Scholar 

  5. 5.

    Luo Z, Su R, Wang W, et al. EV71 infection induces neurodegeneration via activating TLR7 signaling and IL-6 production. PLoSPathog. 2019;15(11):e1008142.

    CAS  Google Scholar 

  6. 6.

    Antona D, Kossorotoff M, Schuffenecker I, et al. Severe paediatric conditions linked with EV-A71 and EV-D68, France, May to October 2016. Euro Surveill. 2016;21(46):30402.

    PubMed Central  Google Scholar 

  7. 7.

    Ooi MH, Wong SC, Lewthwaite P, et al. Clinical features, diagnosis, and management of enterovirus 71. Lancet Neurol. 2010;9(11):1097–105.

    PubMed  Google Scholar 

  8. 8.

    Kok CC. Therapeutic and prevention strategies against human enterovirus 71 infection. World J Virol. 2015;4(2):78–95.

    PubMed  PubMed Central  Google Scholar 

  9. 9.

    Ho M, Chen ER, Hsu KH, et al. An epidemic of enterovirus 71 infection in Taiwan. Taiwan Enterovirus Epidemic Working Group. N Engl J Med. 1999;341(13):929–35.

    CAS  PubMed  Google Scholar 

  10. 10.

    Shih SR, Ho MS, Lin KH, et al. Genetic analysis of enterovirus 71 isolated from fatal and non-fatal cases of hand, foot and mouth disease during an epidemic in Taiwan, 1998. Virus Res. 2000;68(2):127–36.

    CAS  PubMed  Google Scholar 

  11. 11.

    Chan KP, Goh KT, Chong CY, et al. Epidemic hand, foot and mouth disease caused by human enterovirus 71, Singapore. Emerg Infect Dis. 2003;9(1):78–85.

    PubMed  PubMed Central  Google Scholar 

  12. 12.

    Herrero LJ, Lee CS, Hurrelbrink RJ, et al. Molecular epidemiology of enterovirus 71 in peninsular Malaysia, 1997–2000. Arch Virol. 2003;148(7):1369–85.

    CAS  PubMed  Google Scholar 

  13. 13.

    Guan D, van der Sanden S, Zeng H, et al. Population dynamics and genetic diversity of C4 strains of human enterovirus 71 in Mainland China, 1998–2010. PLoS ONE. 2012;7(9):e44386.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Zhu Z, Zhu S, Guo X, et al. Retrospective seroepidemiology indicated that human enterovirus 71 and coxsackievirus A16 circulated wildly in central and southern China before large-scale outbreaks from 2008. Virol J. 2010;7:300.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Wang JR, Tuan YC, Tsai HP, et al. Change of major genotype of enterovirus 71 in outbreaks of hand-foot-and-mouth disease in Taiwan between 1998 and 2000. J Clin Microbiol. 2002;40(1):10–5.

    PubMed  PubMed Central  Google Scholar 

  16. 16.

    Yu H, Chen W, Chang H, et al. Genetic analysis of the VP1 region of enterovirus 71 reveals the emergence of genotype A in central China in 2008. Virus Genes. 2010;41(1):1–4.

    CAS  PubMed  Google Scholar 

  17. 17.

    van der Sanden S, van der Avoort H, Lemey P, et al. Evolutionary trajectory of the VP1 gene of human enterovirus 71 genogroup B and C viruses. J GenVirol. 2010;91(Pt 8):1949–58.

    Google Scholar 

  18. 18.

    Zhang Y, Zhu Z, Yang W, et al. An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol J. 2010;7:94.

    PubMed  PubMed Central  Google Scholar 

  19. 19.

    Lee KY. Enterovirus 71 infection and neurological complications. Korean J Pediatr. 2016;59(10):395–401.

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    Tan X, Huang X, Zhu S, et al. The persistent circulation of enterovirus 71 in People’s Republic of China: causing emerging nationwide epidemics since 2008. PLoS ONE. 2011;6(9):e25662.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Tao Z, Wang H, Xu A. Identification of a C2 subgenogroup strain of enterovirus 71 in a retrospective study in Shandong Province, China, from 1990 to 2010. J Clin Microbiol. 2012;50:1823–4.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Tao Z, Wang H, Li Y, et al. Molecular epidemiology of human enterovirus associated with aseptic meningitis in Shandong Province, China, 2006–2012. PLoS ONE. 2014;9(2):e89766.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Yip CC, Lau SK, Woo PC, et al. Human enterovirus 71 epidemics: what’s next? EmergHealth Threats J. 2013;6:19780.

    Google Scholar 

  24. 24.

    Chia MY, Chiang PS, Chung WY, et al. Epidemiology of enterovirus 71 infections in Taiwan. PediatrNeonatol. 2014;55(4):243–9.

    Google Scholar 

  25. 25.

    Huang SW, Cheng HL, Hsieh HY, et al. Mutations in the non-structural protein region contribute to intra-genotypic evolution of enterovirus 71. J Biomed Sci. 2014;21(1):33.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Chong P, Liu CC, Chow YH, et al. Review of enterovirus 71 vaccines. Clin Infect Dis. 2015;60(5):797–803.

    CAS  PubMed  Google Scholar 

  27. 27.

    Chou AH, Liu CC, Chang JY, et al. Formalin-inactivated EV71 vaccine candidate induced cross-neutralizing antibody against subgenotypes B1, B4, B5 and C4A in adult volunteers. PLoS ONE. 2013;8(11):e79783.

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    Bible JM, Iturriza-Gomara M, Megson B, et al. Molecular epidemiology of human enterovirus 71 in the United Kingdom from 1998 to 2006. J Clin Mocrobiol. 2008;46(10):3192–200.

    CAS  Google Scholar 

  29. 29.

    Huang SC, Hsu YW, Wang HC, et al. Appearance of intratypic recombination of enterovirus 71 in Taiwan from 2002 to 2005. Virus Res. 2008;131(2):250–9.

    CAS  PubMed  Google Scholar 

  30. 30.

    van der Sanden S, van Eek J, Martin DP, et al. Detection of recombination breakpoints in the genomes of human enterovirus 71 strains isolated in the Netherlands in epidemic and non-epidemic years, 1963–2010. Infect Genet Evol. 2011;11(5):886–94.

    PubMed  Google Scholar 

  31. 31.

    McWilliam LE, Cabrerizo M, Cardosa J, et al. The association of recombination events in the founding and emergence of subgenogroup evolutionary lineages of human enterovirus 71. J Virol. 2012;86(5):2676–85.

    Google Scholar 

  32. 32.

    Wen S, Ma D, Lin Y, et al. Complete genome characterization of the 2017 dengue outbreak in Xishuangbanna, a Border City of China, Burma and Laos. Front Cell Infect Microbiol. 2018;8:148.

    PubMed  PubMed Central  Google Scholar 

  33. 33.

    Zhang Y, Tan X, Cui A, et al. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Li Y, Zhou Y, Cheng Y, et al. Effectiveness of EV-A71 vaccination in prevention of paediatric hand, foot, and mouth disease associated with EV-A71 virus infection requiring hospitalisation in Henan, China, 2017–18: a test-negative case-control study. Lancet Child Adolesc Health. 2019;3(10):697–704.

    PubMed  Google Scholar 

  35. 35.

    Weng Y, Chen W, Huang M, et al. Epidemiology and etiology of hand, foot, and mouth disease in Fujian province, 2008–2014. Arch Virol. 2017;162(2):535–42.

    CAS  PubMed  Google Scholar 

  36. 36.

    Liu Y, Fu C, Wu S, et al. A novel finding for enterovirus virulence from the capsid protein VP1 of EV71 circulating in mainland China. Virus Genes. 2014;48(2):260–72.

    CAS  PubMed  Google Scholar 

  37. 37.

    Wu JS, Zhao N, Pan H, et al. Patterns of polymorphism and divergence in the VP1 gene of enterovirus 71 circulating in the Asia-Pacific region between 1994 and 2013. J Virol Methods. 2013;193(2):713–28.

    CAS  PubMed  Google Scholar 

  38. 38.

    Zhu H, Cao Y, Su W, et al. Enterovirus A71 VP1 variation A289T decreases the central nervous system infectivity via attenuation of interactions between VP1 and vimentin in vitro and in vivo. Viruses. 2019;11(5):467.

    CAS  PubMed Central  Google Scholar 

  39. 39.

    Nishimura Y, Lee H, Hafenstein S, et al. Enterovirus 71 binding to PSGL-1 on leukocytes: VP1-145 acts as a molecular switch to control receptor interaction. PLoS Pathog. 2013;9(7):e1003511.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Chen TC, Liu SC, Huang PN, et al. Antiviral activity of pyridyl imidazolidinones against enterovirus 71 variants. J Biome Sci. 2008;15(3):291–300.

    CAS  Google Scholar 

  41. 41.

    Le TV, Nguyen V, Nguyen QH, et al. Molecular epidemiology analysis of enterovirus 71 strains isolated in Dak Lak, Vietnam, 2011–2016. J Med Virol. 2019;91(1):56–64.

    PubMed  Google Scholar 

  42. 42.

    Zhang B, Wu X, Huang K, et al. The variations of VP1 protein might be associated with nervous system symptoms caused by enterovirus 71 infection. BMC Infect Dis. 2014;14:243.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    Kobayashi K, Sudaka Y, Takashino A, et al. Amino acid variation at VP1-145 of Enterovirus 71 determines attachment receptor usage and neurovirulence in human scavenger receptor B2 transgenic mice. J Virol. 2018;92(15):e00681–18.

    PubMed  PubMed Central  Google Scholar 

  44. 44.

    Victorio CB, Xu Y, Ng Q, et al. Cooperative effect of the VP1 amino acids 98E, 145A and 169F in the productive infection of mouse cell lines by enterovirus 71 (BS strain). Emerg Microbes Infect. 2016;5(6):e60.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Yang CH, Liang CT, Jiang ST, et al. A novel murine model expressing a chimeric mSCARB2/hSCARB2 receptor is highly susceptible to oral infection with clinical isolates of Enterovirus 71. J Virol. 2019;93(11):e00183–19.

    PubMed  PubMed Central  Google Scholar 

  46. 46.

    Zhang H, Song Z, Zou J, et al. An infectious clone of enterovirus 71(EV71) that is capable of infecting neonatal immune competent mice without adaptive mutations. Emerg Microbes Infect. 2020;9(1):427–38.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Zhu FC, Meng FY, Li JX, et al. Efficacy, safety, and immunology of an inactivated alum-adjuvant enterovirus 71 vaccine in children in China: a multicentre, randomised, double-blind, placebo-controlled, phase 3 trial. Lancet. 2013;381(9882):2024–32.

    CAS  PubMed  Google Scholar 

  48. 48.

    Mao QY, Wang Y, Bian L, et al. EV71 vaccine, a new tool to control outbreaks of hand, foot and mouth disease (HFMD). Expert Rev Vaccines. 2016;15(5):599–606.

    CAS  PubMed  Google Scholar 

  49. 49.

    Yang B, Liu F, Liao Q, et al. Epidemiology of hand, foot and mouth disease in China, 2008 to 2015 prior to theintroduction of EV-A71 vaccine. Euro Surveill. 2017;22(50):16–00824.

    PubMed Central  Google Scholar 

  50. 50.

    Li R, Liu L, Mo Z, et al. An inactivated enterovirus 71 vaccine in healthy children. N Engl J Med. 2014;370(9):829–37.

    CAS  PubMed  Google Scholar 

  51. 51.

    Zhu F, Xu W, Xia J, et al. Efficacy, safety, and immunogenicity of an enterovirus 71 vaccine in China. N Engl J Med. 2014;370(9):818–28.

    CAS  PubMed  Google Scholar 

  52. 52.

    Hu YM, Wang X, Wang JZ, et al. Immunogenicity, safety, and lot consistency of a novel inactivated enterovirus 71 vaccine in Chinese children aged 6 to 59 months. Clin Vaccine Immunol. 2013;20(12):1805–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Bessaud M, Razafindratsimandresy R, Nougairede A, et al. Molecular comparison and evolutionary analyses of VP1 nucleotide sequences of new African human enterovirus 71 isolates reveal a wide genetic diversity. PLoS ONE. 2014;9(3):e90624.

    PubMed  PubMed Central  Google Scholar 

Download references


The authors acknowledge all of the clinical staff throughout China who provided the VP1 protein-encoding region gene sequences and information concerning the EV71 strains used in this study. We gratefully acknowledge American Journal Experts (AJE) (Verification code: 053B-2FE3-F3E3-887C-4651) and (Verification code: 0CA5-4784-DED2-FED4-0F6A) for providing English language editing services for our manuscript twice.


This work was supported by grants from the National Natural Science Foundation of China (No. 81871709), the National Key Programs for Infectious Diseases of China (No. 2017ZX10103008-005), the Natural Science Foundation of Zhejiang Province, China (No. LY16H200001), the Scientific Research Projects of Education of Zhejiang Province, China (No. Y201534117), and the Scientific Technology Projects of Health and the Family Planning Commission of Zhejiang Province (No. 2015KYB149).

Author information




DC participated in the design of the study. DC, HS and MG performed the data collection and statistical analysis, and DC and HS drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dawei Cui.

Ethics declarations

Ethics approval and consent to participate

This study was in compliance with ethical standards.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, H., Gao, M. & Cui, D. Molecular characteristics of the VP1 region of enterovirus 71 strains in China. Gut Pathog 12, 38 (2020).

Download citation


  • Molecular characteristics
  • Enterovirus 71
  • VP1
  • Mutation
  • Genotype