key: cord-030504-442pnl5s
authors: Sun, Haiyan; Gao, Min; Cui, Dawei
title: Molecular characteristics of the VP1 region of enterovirus 71 strains in China
date: 2020-08-14
journal: Gut Pathog
DOI: 10.1186/s13099-020-00377-2
sha: 
doc_id: 30504
cord_uid: 442pnl5s

BACKGROUND: 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. METHODS: 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. RESULTS: 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%). CONCLUSION: 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.

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 Gen-Bank 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. (H07, FY7VP5, and FY-23 K-B) and 27 known genotypes of EV71 strains were considered reference sequences, as presented in Table 1 .

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.

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.

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) .

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 Table 3 and Fig. 5 .

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 ). 

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 ).

It has been reported that since 2004, the major EV71 strains in mainland China basically belong to the C4 genotype [18, 33] . 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

H, histidine; Q, glutamine; R, arginine; N, asparagine; D, aspartic acid; S, serine; G, glycine; E, glutamic acid; K, lysine; A, alanine; L, leucine; V, valine; I, isoleucine; F, phenylalanine; T, threonine; P, proline 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.

Enterovirus 71 outbreak in the People's Republic of China in 2008

Hand, foot, and mouth disease in China: patterns of spread and transmissibility

Hand, foot, and mouth disease in China, 2008-12: an epidemiological study

Virology, epidemiology, pathogenesis, and control of enterovirus 71

EV71 infection induces neurodegeneration via activating TLR7 signaling and IL-6 production

Severe paediatric conditions linked with EV-A71 and EV-D68, France

Clinical features, diagnosis, and management of enterovirus 71

Therapeutic and prevention strategies against human enterovirus 71 infection

An epidemic of enterovirus 71 infection in Taiwan. Taiwan Enterovirus Epidemic Working Group

Genetic analysis of enterovirus 71 isolated from fatal and non-fatal cases of hand, foot and mouth disease during an epidemic in Taiwan

Epidemic hand, foot and mouth disease caused by human enterovirus 71

Molecular epidemiology of enterovirus 71 in peninsular Malaysia

Population dynamics and genetic diversity of C4 strains of human enterovirus 71 in Mainland China

Retrospective seroepidemiology indicated that human enterovirus 71 and coxsackievirus A16 circulated wildly in central and southern China before large-scale outbreaks from

Change of major genotype of enterovirus 71 in outbreaks of hand-foot-and-mouth disease in Taiwan between 1998 and

Genetic analysis of the VP1 region of enterovirus 71 reveals the emergence of genotype A in central China

Evolutionary trajectory of the VP1 gene of human enterovirus 71 genogroup B and C viruses

An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China

Enterovirus 71 infection and neurological complications

The persistent circulation of enterovirus 71 in People's Republic of China: causing emerging nationwide epidemics since

Identification of a C2 subgenogroup strain of enterovirus 71 in a retrospective study in Shandong Province, China, from 1990 to 2010

Molecular epidemiology of human enterovirus associated with aseptic meningitis in

Human enterovirus 71 epidemics: what's next?

Epidemiology of enterovirus 71 infections in Taiwan

Mutations in the non-structural protein region contribute to intra-genotypic evolution of enterovirus 71

Review of enterovirus 71 vaccines

):e79783. • fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year • At BMC

Molecular epidemiology of human enterovirus 71 in the United Kingdom from 1998 to

Appearance of intratypic recombination of enterovirus 71 in Taiwan from

Detection of recombination breakpoints in the genomes of human enterovirus 71 strains isolated in the Netherlands in epidemic and non-epidemic years

The association of recombination events in the founding and emergence of subgenogroup evolutionary lineages of human enterovirus 71

Complete genome characterization of the 2017 dengue outbreak in Xishuangbanna, a Border City of China

Complete genome analysis of the C4 subgenotype strains of enterovirus 71: predominant recombination C4 viruses persistently circulating in China for 14 years

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

Epidemiology and etiology of hand, foot, and mouth disease in Fujian province

A novel finding for enterovirus virulence from the capsid protein VP1 of EV71 circulating in mainland China

Patterns of polymorphism and divergence in the VP1 gene of enterovirus 71 circulating in the Asia-Pacific region between 1994 and 2013

Enterovirus A71 VP1 variation A289T decreases the central nervous system infectivity via attenuation of interactions between VP1 and vimentin in vitro and in vivo

Enterovirus 71 binding to PSGL-1 on leukocytes: VP1-145 acts as a molecular switch to control receptor interaction

Antiviral activity of pyridyl imidazolidinones against enterovirus 71 variants

Molecular epidemiology analysis of enterovirus 71 strains

The variations of VP1 protein might be associated with nervous system symptoms caused by enterovirus 71 infection

Amino acid variation at VP1-145 of Enterovirus 71 determines attachment receptor usage and neurovirulence in human scavenger receptor B2 transgenic mice

Cooperative effect of the VP1 amino acids 98E, 145A and 169F in the productive infection of mouse cell lines by enterovirus 71 (BS strain)

A novel murine model expressing a chimeric mSCARB2/hSCARB2 receptor is highly susceptible to oral infection with clinical isolates of Enterovirus 71

An infectious clone of enterovirus 71(EV71) that is capable of infecting neonatal immune competent mice without adaptive mutations

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

EV71 vaccine, a new tool to control outbreaks of hand, foot and mouth disease (HFMD)

Epidemiology of hand, foot and mouth disease in China

An inactivated enterovirus 71 vaccine in healthy children

Efficacy, safety, and immunogenicity of an enterovirus 71 vaccine in China

Immunogenicity, safety, and lot consistency of a novel inactivated enterovirus 71 vaccine in Chinese children aged 6 to 59 months

Molecular comparison and evolutionary analyses of VP1 nucleotide sequences of new African human enterovirus 71 isolates reveal a wide genetic diversity

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.Authors' contributions 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.