Antigenic and Genetic Diversity of Human Enterovirus 71 from 2009 to 2012, Taiwan

Different subgenogroups of enterovirus 71 (EV-71) have caused numerous outbreaks of hand, foot, and mouth disease worldwide, especially in the Asia-Pacific region. During the development of a vaccine against EV-71, the genetic and antigenic diversities of EV-71 isolates from Taiwan were analyzed by phylogenetic analyses and neutralization tests. The results showed that the dominant genogroups had changed twice, from B to C and from C to B, between 2009 and 2012. The subgenogroup B5 (B5b cluster) was dominant in 2008-2009 but was replaced by subgenogroup C4 in 2010-2011. From the end of 2011 to 2012, the re-emerging subgenogroup B5 (B5c cluster) was identified as the dominant subgenogroup of EV-71 outbreaks, and subgenogroups C2 and C4 were detected in sporadic cases. Interestingly, the amino acid substitution at position 145 in the VP1 gene was observed in some strains isolated from patients with acute flaccid paralysis. Furthermore, thirty-five strains and their corresponding serum samples were used to analyze the cross-protections and antigenic diversities among different subgenogroups (C4a, C5, B4, B5b, B5c, and C2-like) of EV-71. Evident antigenic diversity existed only for the C2-like subgenogroup, which was not effectively neutralized by other serum samples. In contrast, the anti-C2-like serum sample showed broad cross-reactivity against all other subgenogroups. Therefore, these results may provide valuable information for the selection of EV-71 vaccine candidates and the evolution of EV-71 subgenogroups in Taiwan from 2009 to 2012.     

Evolutionary Dynamics and Temporal/Geographical Correlates of Recombination in the Human Enterovirus Echovirus Types 9, 11, and 30

The relationship between virus evolution and recombination in species B human enteroviruses was investigated through large-scale genetic analysis of echovirus type 9 (E9) and E11 isolates (n = 85 and 116) from 16 European, African, and Asian countries between 1995 and 2008. Cluster 1 E9 isolates and genotype D5 and A E11 isolates showed evidence of frequent recombination between the VP1 and 3Dpol regions, the latter falling into 23 (E9) and 43 (E11) clades interspersed phylogenetically with 46 3Dpol clades of E30 and with those of other species B serotypes. Remarkably, only 2 of the 112 3Dpol clades were shared by more than one serotype (E11 and E30), demonstrating an extremely large and genetically heterogeneous recombination pool of species B nonstructural-region variants. The likelihood of recombination increased with geographical separation and time, and both were correlated with VP1 divergence, whose substitution rates allowed recombination half-lives of 1.3, 9.8, and 3.1 years, respectively, for E9, E11, and E30 to be calculated. These marked differences in recombination dynamics matched epidemiological patterns of periodic epidemic cycles of 2 to 3 (E9) and 5 to 6 (E30) years and the longer-term endemic pattern of E11 infections. Phylotemporal analysis using a Bayesian Markov chain Monte Carlo method, which placed recombination events within the evolutionary reconstruction of VP1, showed a close relationship with VP1 lineage expansion, with defined recombination events that correlated with their epidemiological periodicity. Whether recombination events contribute directly to changes in transmissibility that drive epidemic behavior or occur stochastically during periodic population bottlenecks is an unresolved issue vital to future understanding of enterovirus molecular epidemiology and pathogenesis.Human enteroviruses (HEV) are single-stranded, positive-sense RNA viruses in the virus family Picornaviridae. Infections with these viruses are enterically transmitted and normally cause subclinical or mild symptoms, but they are also capable of causing a wide array of often severe disease presentations, including aseptic meningitis, encephalitis, and acute flaccid paralysis. Echovirus type 9 (E9) and E11 are serotypes of species B enteroviruses and are among the most common etiological agents of aseptic meningitis (14). The first epidemiological study of E9 investigated isolates obtained from sick children in 1995 (24). E9 is chiefly associated with mild infections affecting children over the age of 5 years and teenagers, although it can cause more severe syndromes in neonates and immunosuppressed patients (14). E11 infections, which predominantly affect infants less than 1 year old, have been associated with large outbreaks of uveitis in infants (15, 21), sepsis, and other neonatal systemic illnesses with high mortality rates (14, 23).Circulating strains and genotypes of E9 (2, 8, 11, 14) and E11 (4, 7, 21) have been characterized in many different countries through analysis of the structural genome regions, principally VP1. In common with other mammalian RNA viruses, both E9 and E11 show rapid accumulation of nucleotide substitutions over time, but their epidemiologies are distinct. E9 is associated with widespread, large-scale seasonal outbreaks and displays a regular epidemic pattern of outbreaks occurring approximately every 3 years (8, 14). Outbreaks of E11 are less common, occur irregularly, and often last for several years (14, 25).In contrast to time-related diversification observed in the capsid-encoding regions of the HEV genome, nonstructural region genes are subject to frequent recombination (7, 18, 19, 26, 30), leading to the concept of separate, modular evolution of the structural and nonstructural regions of the HEV genome (20). Recombination in the nonstructural gene region of HEV is limited to members of the same species (27). In a recent large-scale investigation of echovirus type 30 (E30) molecular epidemiology, phylogenetic analysis of the 3Dpol region revealed the existence of a number of discrete, bootstrap-supported clades, allowing circulating E30 variants to be classified into a series of distinct recombinant forms (RFs) (22). Individual RFs became very rapidly disseminated over large geographical distances through repeated cycles of emergence, dominance, and disappearance in a 3- to 5-year cycle.In the current study, we carried out parallel investigations of the molecular epidemiology and evolution of clinically presenting isolates of E9 and E11 collected in 16 countries in Europe, central Asia, Southeast Asia, and West Africa over a 14-year period. By comprehensive sequencing of these isolates in the structural (VP1) and nonstructural (3Dpol) genome regions, it was possible to explore the dynamics of sequence diversification and recombination. The turnover of individual recombinant groups and the phylogeographical correlates of recombination were determined and compared to those of E30, providing new insights into the nature of the global spread of these important viral pathogens and the role of recombination in enteroviral evolution.     

Frequency and dynamics of recombination within different species of human enteroviruses

Enteroviruses are members of the family Picornaviridae that cause widespread infections in human and other mammalian populations. Enteroviruses are genetically and antigenically highly variable, and recombination within and between serotypes contributes to their genetic diversity. To investigate the dynamics of the recombination process, sequence phylogenies between three regions of the genome (VP4, VP1, and 3Dpol) were compared among species A and B enterovirus variants detected in a human population-based survey in Scotland between 2000 and 2001, along with contemporary virus isolates collected in the same geographical region. This analysis used novel bioinformatic methods to quantify phylogenetic compatibility and correlations with serotype assignments of evolutionary trees constructed for different regions of the enterovirus genome. Species B enteroviruses showed much more frequent, time-correlated recombination events than those found for species A, despite the equivalence in population sampling, concordant with a linkage analysis of previously characterized enterovirus sequences obtained over longer collection periods. An analysis of recombination among complete genome sequences by computation of a phylogenetic compatibility matrix (PCM) demonstrated sharply defined boundaries between the VP2/VP3/VP1 block and sequences to either side in phylogenetic compatibility. The PCM also revealed equivalent or frequently greater degrees of incompatibility between different parts within the nonstructural region (2A-3D), indicating the occurrence of extensive recombination events in the past evolution of this part of the genome. Together, these findings provide new insights into the dynamics of species A and B enterovirus recombination and evolution and into the contribution of structured sampling to documenting reservoirs, emergence, and spread of novel recombinant forms in human populations.     

Monitoring Antigenic Variations of Enterovirus 71: Implications for Virus Surveillance and Vaccine Development

Enterovirus 71 (EV71) causes life-threatening epidemics in Asia and can be phylogenetically classified into three major genogroups (A∼C) including 11 genotypes (A, B1∼B5, and C1∼C5). Recently, EV71 epidemics occurred cyclically in Taiwan with different genotypes. In recent years, human studies using post-infection sera obtained from children have detected antigenic variations among different EV71 strains. Therefore, surveillance of enterovirus 71 should include phylogenetic and antigenic analysis. Due to limitation of sera available from children with EV71 primary infection, suitable animal models should be developed to generate a panel of antisera for monitoring EV71 antigenic variations. Twelve reference strains representing the 11 EV71 genotypes were grown in rhabdomyosarcoma cells. Infectious EV71 particles were purified and collected to immunize rabbits. The rabbit antisera were then employed to measure neutralizing antibody titers against the 12 reference strains and 5 recent strains. Rabbits immunized with genogroup B and C viruses consistently have a lower neutralizing antibody titers against genogroup A (≧8-fold difference) and antigenic variations between genogroup B and C viruses can be detected but did not have a clear pattern, which are consistent with previous human studies. Comparison between human and rabbit neutralizing antibody profiles, the results showed that ≧8-fold difference in rabbit cross-reactive antibody ratios could be used to screen EV71 isolates for identifying potential antigenic variants. In conclusion, a rabbit model was developed to monitor antigenic variations of EV71, which are critical to select vaccine strains and predict epidemics.     

A Novel Universal Neutralizing Monoclonal Antibody against Enterovirus 71 That Targets the Highly Conserved “Knob” Region of VP3 Protein

Hand, foot and mouth disease caused by enterovirus 71(EV71) leads to the majority of neurological complications and death in young children. While putative inactivated vaccines are only now undergoing clinical trials, no specific treatment options exist yet. Ideally, EV71 specific intravenous immunoglobulins could be developed for targeted treatment of severe cases. To date, only a single universally neutralizing monoclonal antibody against a conserved linear epitope of VP1 has been identified. Other enteroviruses have been shown to possess major conformational neutralizing epitopes on both the VP2 and VP3 capsid proteins. Hence, we attempted to isolate such neutralizing antibodies against conformational epitopes for their potential in the treatment of infection as well as differential diagnosis and vaccine optimization. Here we describe a universal neutralizing monoclonal antibody that recognizes a conserved conformational epitope of EV71 which was mapped using escape mutants. Eight escape mutants from different subgenogroups (A, B2, B4, C2, C4) were rescued; they harbored three essential mutations either at amino acid positions 59, 62 or 67 of the VP3 protein which are all situated in the “knob” region. The escape mutant phenotype could be mimicked by incorporating these mutations into reverse genetically engineered viruses showing that P59L, A62D, A62P and E67D abolish both monoclonal antibody binding and neutralization activity. This is the first conformational neutralization epitope mapped on VP3 for EV71.     

The Role of VP1 Amino Acid Residue 145 of Enterovirus 71 in Viral Fitness and Pathogenesis in a Cynomolgus Monkey Model

Enterovirus 71 (EV71), a major causative agent of hand, foot, and mouth disease, occasionally causes severe neurological symptoms. We identified P-selectin glycoprotein ligand-1 (PSGL-1) as an EV71 receptor and found that an amino acid residue 145 in the capsid protein VP1 (VP1-145) defined PSGL-1-binding (PB) and PSGL-1-nonbinding (non-PB) phenotypes of EV71. However, the role of PSGL-1-dependent EV71 replication in neuropathogenesis remains poorly understood. In this study, we investigated viral replication, genetic stability, and the pathogenicity of PB and non-PB strains of EV71 in a cynomolgus monkey model. Monkeys were intravenously inoculated with cDNA-derived PB and non-PB strains of EV71, EV71-02363-EG and EV71-02363-KE strains, respectively, with two amino acid differences at VP1-98 and VP1-145. Mild neurological symptoms, transient lymphocytopenia, and inflammatory cytokine responses, were found predominantly in the 02363-KE-inoculated monkeys. During the early stage of infection, viruses were frequently detected in clinical samples from 02363-KE-inoculated monkeys but rarely in samples from 02363-EG-inoculated monkeys. Histopathological analysis of central nervous system (CNS) tissues at 10 days postinfection revealed that 02363-KE induced neuropathogenesis more efficiently than that induced by 02363-EG. After inoculation with 02363-EG, almost all EV71 variants detected in clinical samples, CNS, and non-CNS tissues, possessed a G to E amino acid substitution at VP1-145, suggesting a strong in vivo selection of VP1-145E variants and CNS spread presumably in a PSGL-1-independent manner. EV71 variants with VP1-145G were identified only in peripheral blood mononuclear cells in two out of four 02363-EG-inoculated monkeys. Thus, VP1-145E variants are mainly responsible for the development of viremia and neuropathogenesis in a non-human primate model, further suggesting the in vivo involvement of amino acid polymorphism at VP1-145 in cell-specific viral replication, in vivo fitness, and pathogenesis in EV71-infected individuals.     

肠道病毒71型安徽、河南株的分离与VP1序列进化分析

旨在研究手足口病患者中肠道病毒71型分离株的病毒基因型特征。采集手足口病患者的粪便标本,进行病毒分离和逆转录-聚合酶链式反应(RT-PCR)特异性扩增进行鉴定,同时选取其中9株EV71分离株,对其抗原决定簇部位VP1区进行核酸序列测定,并参考EV71 A、B、C各基因型的参考株和以往中国EV71的分离株进行同源分析和构建系统发生树。结果显示,所分析的9株病毒株均为C4亚型,3株安徽株H7、H8和H9的VP1序列相似度很高(≥98.8%,其中H7、H9的相似度为100.0%),4株河南株H3、H4、H5和H6相似度较高(≥98.4%,H3、H4和H5≥99.6%,其中H3、H4的相似度为100.0%),它们同河南株H1、H5的相似度也较高(≥97.2%),河南株H2虽然与其他河南株具有较高的序列相似度,但进化分析表明,其与安徽株同源性较高。结果表明,安徽株H7、H8和H9株变异速率明显加快,这可能导致了手足口病在安徽省的率先爆发与大流行,河南株H2最初可能由安徽传入河南。     

Identification of Site-Specific Adaptations Conferring Increased Neural Cell Tropism during Human Enterovirus 71 Infection

Enterovirus 71 (EV71) is one of the most virulent enteroviruses, but the specific molecular features that enhance its ability to disseminate in humans remain unknown. We analyzed the genomic features of EV71 in an immunocompromised host with disseminated disease according to the different sites of infection. Comparison of five full-length genomes sequenced directly from respiratory, gastrointestinal, nervous system, and blood specimens revealed three nucleotide changes that occurred within a five-day period: a non-conservative amino acid change in VP1 located within the BC loop (L97R), a region considered as an immunogenic site and possibly important in poliovirus host adaptation; a conservative amino acid substitution in protein 2B (A38V); and a silent mutation in protein 3D (L175). Infectious clones were constructed using both BrCr (lineage A) and the clinical strain (lineage C) backgrounds containing either one or both non-synonymous mutations. In vitro cell tropism and competition assays revealed that the VP1₉₇ Leu to Arg substitution within the BC loop conferred a replicative advantage in SH-SY5Y cells of neuroblastoma origin. Interestingly, this mutation was frequently associated in vitro with a second non-conservative mutation (E167G or E167A) in the VP1 EF loop in neuroblastoma cells. Comparative models of these EV71 VP1 variants were built to determine how the substitutions might affect VP1 structure and/or interactions with host cells and suggest that, while no significant structural changes were observed, the substitutions may alter interactions with host cell receptors. Taken together, our results show that the VP1 BC loop region of EV71 plays a critical role in cell tropism independent of EV71 lineage and, thus, may have contributed to dissemination and neurotropism in the immunocompromised patient.     

Enterovirus Migration Patterns between France and Tunisia

The enterovirus (EV) types echovirus (E-) 5, E-9, and E-18, and coxsackievirus (CV-) A9 are infrequently reported in human diseases and their epidemiologic features are poorly defined. Virus transmission patterns between countries have been estimated with phylogenetic data derived from the 1D/VP1 and 3CD gene sequences of a sample of 74 strains obtained in France (2000–2012) and Tunisia (2011–2013) and from the publicly available sequences. The EV types (E-5, E-9, and E-18) exhibited a lower worldwide genetic diversity (respective number of genogroups: 4, 5, and 3) in comparison to CV-A9 (n = 10). The phylogenetic trees estimated with both 1D/VP1 and 3CD sequence data showed variations in the number of co-circulating lineages over the last 20 years among the four EV types. Despite the low number of genogroups in E-18, the virus exhibited the highest number of recombinant 3CD lineages (n = 10) versus 4 (E-5) to 8 (E-9). The phylogenies provided evidence of multiple transportation events between France and Tunisia involving E-5, E-9, E-18, and CV-A9 strains. Virus spread events between France and 17 other countries in five continents had high probabilities of occurrence as those between Tunisia and two European countries other than France. All transportation events were supported by BF values > 10. Inferring the source of virus transmission from phylogenetic data may provide insights into the patterns of sporadic and epidemic diseases caused by EVs.     

肠道病毒71型外壳蛋白VP1全长在大肠杆菌中的高效表达及初步活性测定

目的:重组表达肠道病毒71型(EV71)外壳蛋白VP1全长,用于研制血清学检测试剂和疫苗研发。方法:在获得EV71全长基因并测序正确的基础上,将外壳蛋白VP1全长基因克隆到表达载体pET28a(+)上,构建重组表达质粒pET28a(+)/VP1,转化大肠杆菌BL21,IPTG诱导表达,利用Ni2+亲和层析柱对重组蛋白进行纯化,采用双抗原夹心检测技术评价重组抗原与27份EV71抗体阳性血清和18份阴性血清的反应情况。结果:重组EV71-VP1蛋白在大肠杆菌中诱导6 h后可获得高效表达,能与27份EV71抗体阳性血清中的21份发生阳性反应,EV71双抗原夹心检测与中和血清测试结果具有很好的一致性(P0.05)。结论:实现了肠道病毒71型外壳蛋白VP1的高效表达,为肠道病毒71型诊断试剂和疫苗的研究奠定了基础。     

2008～2010年河南省人肠道病毒71型基因特征和重组分析

对河南省2008～2010年河南省人肠道病毒71型进行基因特征及重组特点研究。对河南省2008～2010年分离的5株肠道病毒EV71型构建VP1序列系统进化树并分析其全基因组序列的重组特点。VP1序列系统进化分析显示2008～2010年河南株均属于C4基因型的C4a亚群,Bootscan分析和5’NCR、P1、P2、P3区的进化树证实C4基因型在2A～2B处存在EV71的B基因型和C基因型的型内重组及在3B～3C处存在EV71的B基因型和CA16/G-10间的型间重组。2008～2010年河南EV71分离株为C4基因型的C4a亚群,与2004年以来的中国大陆优势株流行趋势完全一致,EV71C4基因型存在基因型内和型间双重组现象。     

A retrospective overview of enterovirus infection diagnosis and molecular epidemiology in the public hospitals of Marseille, France (1985-2005)

Human enteroviruses (HEV) are frequent human pathogens and, associated in particular with large outbreaks of aseptic meningitis. Here, we have compiled a database of clinical HEV isolates from the Public Hospitals of Marseille, from 1985 to 2005. Amongst 654 isolates that could be characterized by complete sequencing of the VP1 gene, 98% belonged to species HEV-B; the most frequently isolated serotypes were Echovirus E30, E11, E7, E6 and E4. The high incidence of E30 and the recent emergence of E13 are consistent with reports worldwide and peak HEV isolation occurred mostly in the late spring and summer months. The proportion of echoviruses has decreased across the years, while that of coxsackieviruses has increased. Stool (the most frequent sample type) allowed detection of all identified serotypes. MRC5 (Human lung fibroblasts) cell line was the most conducive cell line for HEV isolation (84.9% of 10 most common serotype isolates, 96.3% in association with BGM (Buffalo green monkey kidney cells)). Previous seroneutralization-based serotype identification demonstrated 55.4% accuracy when compared with molecular VP1 analysis. Our analysis of a large number of clinical strains over 20 years reinforced the validity of VP1 serotyping and showed that comparative p-distance scores can be coupled with phylogenetic analysis to provide non-ambiguous serotype identification. Phylogenetic analysis in the VP1, 2C and 3D regions also provided evidence for recombination events amongst clinical isolates. In particular, it identified isolates with dissimilar VP1 but almost identical nonstructural regions.     

Population Dynamics and Genetic Diversity of C4 Strains of Human Enterovirus 71 in Mainland China, 1998–2010

Background

Since 1997, several countries within the Asian Pacific region have been affected by one or more massive outbreaks of Hand Foot and Mouth Disease (HFMD). Virus typing experiments revealed that these outbreaks were caused by strains of human enterovirus 71 (EV71) belonging to several different, recently emerged subgenogroups. In mainland China, a different situation was observed. The first outbreak, localized in Shangdong Province, was reported in 2007, and was followed by a wide-spread outbreak in mainland China in 2008. Since then, numbers of reported HFMD cases have been persistently high.

Methodology/Principal Findings

To gain insight in the epidemiological behavior of EV71 in China, we studied genetic diversity and EV71 population dynamics to address whether the increase in number of reported EV71 infections reflects a real increase in viral spread or is just the result of increased awareness and surveillance. We used systematically collected VP1 gene sequences of 257 EV71 strains collected in Guangdong province from 2008 to 2010 as part of HFMD surveillance activities, and supplemented them with 305 GenBank EV71 reference stains collected in China from 1998 to 2010. All isolates from Guangdong Province belonged to subgenogroup C4. Viral population dynamics indicated that the increased reporting of HFMD in China since 2007 reflects a real increase in viral spread and continued replacement of viral lineages through time. Amino acid sequence comparisons revealed substitution of amino acid in residues 22, 145 and 289 through time regularly with the VP1 gene of EV71 strains isolated in mainland China from 1998 to 2010.

Conclusions

EV71 strains isolated in mainland China mainly belonged to subgenogroup C4. There was exponential growth of the EV71 virus population in 2007 and 2008. There was amino acid substitution through time regularly with the VP1 gene which possibly increased viral spread and/or ability of the virus to circulate persistently among the Chinese population.     

Molecular Evolution of the Human Enteroviruses: Correlation of Serotype with VP1 Sequence and Application to Picornavirus Classification

Sixty-six human enterovirus serotypes have been identified by serum neutralization, but the molecular determinants of the serotypes are unknown. Since the picornavirus VP1 protein contains a number of neutralization domains, we hypothesized that the VP1 sequence should correspond with neutralization (serotype) and, hence, with phylogenetic lineage. To test this hypothesis and to analyze the phylogenetic relationships among the human enteroviruses, we determined the complete VP1 sequences of the prototype strains of 47 human enterovirus serotypes and 10 antigenic variants. Our sequences, together with those available from GenBank, comprise a database of complete VP1 sequences for all 66 human enterovirus serotypes plus additional strains of seven serotypes. Phylogenetic trees constructed from complete VP1 sequences produced the same four major clusters as published trees based on partial VP2 sequences; in contrast to the VP2 trees, however, in the VP1 trees strains of the same serotype were always monophyletic. In pairwise comparisons of complete VP1 sequences, enteroviruses of the same serotype were clearly distinguished from those of heterologous serotypes, and the limits of intraserotypic divergence appeared to be about 25% nucleotide sequence difference or 12% amino acid sequence difference. Pairwise comparisons suggested that coxsackie A11 and A15 viruses should be classified as strains of the same serotype, as should coxsackie A13 and A18 viruses. Pairwise identity scores also distinguished between enteroviruses of different clusters and enteroviruses from picornaviruses of different genera. The data suggest that VP1 sequence comparisons may be valuable in enterovirus typing and in picornavirus taxonomy by assisting in the genus assignment of unclassified picornaviruses.Human enteroviruses (family Picornaviridae) infect millions of people worldwide each year, resulting in a wide range of clinical outcomes ranging from inapparent infection to mild respiratory illness (common cold), hand-foot-and-mouth disease, acute hemorrhagic conjunctivitis, aseptic meningitis, myocarditis, severe neonatal sepsis-like disease, and acute flaccid paralysis (reviewed in references 43 and 45). In the United States, enteroviruses are responsible for 30,000 to 50,000 meningitis hospitalizations per year as a result of 30 million to 50 million infections. Serologic studies have distinguished 66 human enterovirus serotypes on the basis of an antibody neutralization test (43), and additional antigenic variants have been defined within several of the serotypes on the basis of reduced or nonreciprocal cross-neutralization between prototype and variant strains (6, 8, 68, 71, 72). On the basis of their pathogenesis in humans and experimental animals, the enteroviruses were originally classified into four groups, polioviruses, coxsackie A viruses (CA), coxsackie B viruses (CB), and echoviruses, but it was quickly realized that there were significant overlaps in the biological properties of viruses in the different groups (8). The more recently isolated enteroviruses have been named with a system of consecutive numbers: EV68, EV69, EV70, and EV71 (42).A comparison of nucleotide and deduced amino acid sequences at the 5′ end of VP2 has identified four major phylogenetic groups within the Enterovirus genus: CA16-like viruses (cluster A), a CB-like group containing all CB and echoviruses as well as CA9 and EV69 (cluster B), poliovirus-like viruses (cluster C), and EV68 and EV70 (cluster D) (23, 24, 49, 53, 54, 73). However, pairwise alignments and phylogenetic analyses within these groups demonstrated that the VP2 sequence does not fully correlate with serotype, as viruses known to belong to the same serotype often failed to cluster together (2, 49). (E22 and E23 are genetically distinct from enteroviruses [24], and their reclassification into a separate genus has been proposed [45]).VP1 is the most external and immunodominant of the picornavirus capsid proteins (58). A number of major neutralization sites reside in the VP1 proteins of many picornaviruses (reviewed in references 40 and 44), but the specific epitopes responsible for serotype specificity and intratypic variation have not been identified. Similarly, the genetic correlates of serotype identity remain unknown. If the important serotype-specific neutralization sites reside in VP1, then the VP1 sequence or some portion thereof would be predicted to correlate with serotype. Studies on the three serotypes of poliovirus have shown that a partial VP1 sequence correlates well with serotype (32). In addition, genetic lineages based on the VP1 sequence can be used to define poliovirus reservoirs and chains of transmission (reviewed in reference 30). To test whether the VP1 sequence might be applied to the classification of nonpolio enteroviruses and to the analysis of the phylogenetic relationships among the human enteroviruses, we determined the complete VP1 nucleotide sequences for 47 human enterovirus prototypes and 10 well-characterized antigenic variants. These data, together with previously available sequences, comprise a database of complete VP1 sequences for all known human enterovirus serotypes and 12 natural antigenic variants. This database will be useful for molecular epidemiologic studies of enteroviral disease outbreaks, to obtain a better understanding of the genetic correlates of serotype, and for the development of enteroviral molecular diagnostic reagents.     

Genetic Diversity of Human Enterovirus 68 Strains Isolated in Kenya Using the Hypervariable 3′- End of VP1 Gene

Reports of increasing worldwide circulation of human enterovirus-68 (EV68) are well documented. Despite health concerns posed by resurgence of these viruses, little is known about EV68 strains circulating in Kenya. In this study, we characterized 13 EV68 strains isolated in Kenya between 2008 and 2011 based on the Hypervariable 3′- end of the VP1 gene. Viral RNA was extracted from the isolates and partial VP1 gene amplified by RT-PCR, followed by nucleotide sequencing. Alignment of deduced amino acid sequences revealed substitutions in Kenyan EV68 isolates absent in the prototype reference strain (Fermon). The majority of these changes were present in the BC and DE-loop regions, which are associated with viral antigenicity and virulence. The Kenyan strains exhibited high sequence homology with respect to those from other countries. Natural selection analysis based on the VP1 region showed that the Kenyan EV68 isolates were under purifying selection. Phylogenetic analysis revealed that majority (84.6%) of the Kenyan strains belonged to clade A, while a minority belonged to clades B and C. Overall, our results illustrate that although EV68 strains isolated in Kenya were genetically and antigenically divergent from the prototype strain (Fermon), they were closely related to those circulating in other countries, suggesting worldwide transmissibility. Further, the presence of shared mutations by Kenyan EV68 strains and those isolated in other countries, indicates evolution in the VP1 region may be contributing to increased worldwide detection of the viruses. This is the first study to document circulation of EV68 in Kenya.     

Co-Circulation and Genomic Recombination of Coxsackievirus A16 and Enterovirus 71 during a Large Outbreak of Hand,Foot, and Mouth Disease in Central China

A total of 1844 patients with hand, foot, and mouth disease (HFMD), most of them were children of age 1–3-year-old, in Central China were hospitalized from 2011 to 2012. Among them, 422 were infected with coxsackievirus A16 (CVA16), 334 were infected with enterovirus 71 (EV71), 38 were co-infected with EV71 and CVA16, and 35 were infected with other enteroviruses. Molecular epidemiology analysis revealed that EV71 and CVA16 were detected year-round, but EV71 circulated mainly in July and CVA16 circulated predominantly in November, and incidence of HFMD was reduced in January and February and increased in March. Clinical data showed that hyperglycemia and neurologic complications were significantly higher in EV71-infected patients, while upper respiratory tract infection and C-reactive protein were significantly higher in CVA16-associated patients. 124 EV71 and 80 CVA16 strains were isolated, among them 56 and 68 EV71 strains were C4a and C4b, while 25 and 55 CVA16 strains were B1a and B1b, respectively. Similarity plots and bootscan analyses based on entire genomic sequences revealed that the three C4a sub-genotype EV71 strains were recombinant with C4b sub-genotype EV71 in 2B–2C region, and the three CVA16 strains were recombinant with EV71 in 2A–2B region. Thus, CVA16 and EV71 were the major causative agents in a large HFMD outbreak in Central China. HFMD incidence was high for children among household contact and was detected year-round, but outbreak was seasonal dependent. CVA16 B1b and EV71 C4b reemerged and caused a large epidemic in China after a quiet period of many years. Moreover, EV71 and CVA16 were co-circulated during the outbreak, which may have contributed to the genomic recombination between the pathogens. It should gain more attention as there may be an upward trend in co-circulation of the two pathogens globally and the new role recombination plays in the emergence of new enterovirus variants.     

深圳地区2001～2004年肠道病毒71型部分VP1区基因分析

肠道病毒71型（Enterovirus type71,EV71）自1974年Sehmidt等人首次报道从美国加利福尼亚暴发的表现为中枢神经系统症状的患者标本中分离到后,世界上许多国家相继报道了EV71在不同地区的流行。近年来EV71在亚洲地区的流行呈上升趋势。目前已知EV71的感染可以导致手足口病、无菌性脑膜炎、脑炎和脊髓灰质炎样的麻痹性疾病等多种与神经系统相关的疾病。EV71已经取代脊髓灰质炎病毒成为肠道病毒中最受瞩目的成员。根据病毒衣壳蛋白VP1核苷酸序列的差异,可将EV71分为A、B、C3个基因型。     

Phylogenetic Analysis of Enterovirus 71 Strains Isolated during Linked Epidemics in Malaysia, Singapore, and Western Australia

Enterovirus 71 (EV71) is a frequent cause of hand, foot, and mouth disease (HFMD) epidemics associated with severe neurological sequelae in a small proportion of cases. There has been a significant increase in EV71 epidemic activity throughout the Asia-Pacific region since 1997. Recent HFMD epidemics in this region have been associated with a severe form of brainstem encephalitis associated with pulmonary edema and high case fatality rates. In this study, we show that four genetic lineages of EV71 have been prevalent in the Asia-Pacific region since 1997, including two previously undescribed genogroups (B3 and B4). Furthermore, we show that viruses belonging to genogroups B3 and B4 have circulated endemically in Southeast Asia during this period and have been the primary cause of several large HFMD or encephalitis epidemics in Malaysia, Singapore, and Western Australia.     

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

Background

Clinical manifestations of enterovirus 71 (EV71) range from herpangina, hand-foot-and-mouth disease (HFMD), to severe neurological complications. Unlike the situation of switching genotypes seen in EV71 outbreaks during 1998–2008 in Taiwan, genotype B5 was responsible for two large outbreaks in 2008 and 2012, respectively. In China, by contrast, EV71 often persists as a single genotype in the population and causes frequent outbreaks. To investigate genetic changes in viral evolution, complete EV71 genome sequences were used to analyze the intra-genotypic evolution pattern in Taiwan, China, and the Netherlands.

Results

Genotype B5 was predominant in Taiwan’s 2008 outbreak and was re-emergent in 2012. EV71 strains from both outbreaks were phylogenetically segregated into two lineages containing fourteen non-synonymous substitutions predominantly in the non-structural protein coding region. In China, genotype C4 was first seen in 1998 and caused the latest large outbreak in 2008. Unlike shifting genotypes in Taiwan, genotype C4 persisted with progressive drift through time. A majority of non-synonymous mutations occurred in residues located in the non-structural coding region, showing annual increases. Interestingly, genotype B1/B2 in the Netherlands showed another stepwise evolution with dramatic EV71 activity increase in 1986. Phylogeny of the VP1 coding region in 1971–1986 exhibited similar lineage turnover with genotype C4 in China; however, phylogeny of the 3D-encoding region indicated separate lineage appearing after 1983, suggesting that the 3D-encoding region of genotype B2 was derived from an unidentified ancestor that contributed to intra-genotypic evolution in the Netherlands.

Conclusions

Unlike VP1 coding sequences long used for phylogenetic study of enteroviruses due to expected host immune escape, our study emphasizes a dominant role of non-synonymous mutations in non-structural protein regions that contribute to (re-)emergent genotypes in continuous stepwise evolution. Dozens of amino acid substitutions, especially in non-structural proteins, were identified via genetic changes driven through intra-genotypic evolution worldwide. These identified substitutions appeared to increase viral fitness in the population, affording valuable insights not only for viral evolution but also for prevention, control, and vaccine against EV71 infection.     

Hand foot and mouth disease due to enterovirus 71 in Malaysia

Hand foot and mouth disease is a febrile sickness complex characterized by cutaneous eruption (exanthem) on the palms and soles with simultaneous occurrence of muco-cutanous vesiculo-ulcerative lesions (enanthem) affecting the mouth.The illness is caused by a number of enteroviruses with coxsackievirus A16 and enterovirus 71 as the main causative agents.Human enterovirus 71 (EV71) belongs to the species Human enterovirus A under the genus Enterovirus within the family Picornaviridae.EV71 has been associated with an array of clinical diseases including hand foot and mouth disease (HFMD),aseptic meningitis,encephalitis and poliomyelitis-like acute flaccid paralysis.A large outbreak of HFMD due to highly neurovirulent EV71 emerged in Malaysia in 1997,and caused 41deaths amongst young children.In late 2000,a recurrence of an outbreak of HFMD occurred in Malaysia with S fatalities in peninsular Malaysia.Outbreak of HFMD due to EV71 recurred in 2003 with an unknown number of cases and mortalities.A similar outbreak of HFMD with 2 recorded deaths in young children occurred in peninsular Malaysia in late 2005 and this was followed by a larger outbreak in Sarawak (Malaysian Borneo) with 6 reported fatalities in the early part of 2006.The current on-going outbreak of HFMD started in peninsular Malaysia in epidemiological week 12 of 2010.As with other HFMD outbreaks in Malaysia,both EV71 and CA16 were the main aetiological viruses isolated.In similarity with the HFMD outbreak in 2005,the isolation of CA16 preceded the appearance of EV71.Based on the VP 1 gene nucleotide sequences,4 sub-genogroups of EV71 (C1,C2,B3 and B4) co-circulated and caused the outbreak of hand,foot and mouth disease in peninsular Malaysia in 1997.Two sub-genogroups (C1 and B4) were noted to cause the outbreak in 2000 in both peninsular Malaysia and Sarawak.EV71 of sub-genogroup B5 with smaller contribution from sub-genogroup C1 caused the outbreak in 2003.In the 2005 outbreak,besides the EV71 strains of sub-genogroup C1,EV71 strains belonging to sub-genogroup B5 were isolated but formed a cluster which was distinct from the EV71 strains from the sub-genogroup B5 isolated in 2003.The four EV71 strains isolated from clinical specimens of patients with hand,foot and mouth disease in the Sarawak outbreak in early 2006 also belonged to sub-genogroup B5.Phylogenetic analysis of the VP1 gene suggests that the EV71 strains causing the outbreak in Sarawak could have originated from peninsular Malaysia.Epidemiological and molecular data since 1997 show the recurrence of HFMD due to EV71 in Malaysia every 2 to 4 years.In each of the past outbreaks,more than one sub-genogroup of the virus co-circulate.