Inferring nonneutral evolution from contrasting patterns of polymorphisms and divergences in different protein coding regions of enterovirus 71 circulating in Taiwan during 1998-2003

Background  

Enterovirus (EV) 71 is one of the common causative agents for hand, foot, and, mouth disease (HFMD). In recent years, the virus caused several outbreaks with high numbers of deaths and severe neurological complications. Despite the importance of these epidemics, several aspects of the evolutionary and epidemiological dynamics, including viral nucleotide variations within and between different outbreaks, rates of change in immune-related structural regions vs. non-structural regions, and forces driving the evolution of EV71, are still not clear.     

Incidence rates of enterovirus 71 infections in young children during a nationwide epidemic in Taiwan, 2008-09

Objective

Enterovirus 71 (EV71) is causing life-threatening outbreaks in tropical Asia. In Taiwan and other tropical Asian countries, although nationwide EV71 epidemics occur cyclically, age-specific incidence rates of EV71 infections that are critical to estimate disease burden and design vaccine trials are not clear. A nationwide EV71 epidemic occurred in 2008–09 in Taiwan, which provided a unique opportunity to estimate age-specific incidence rates of EV71 infections.

Study Design

We prospectively recruited 749 healthy neonates and conducted follow-ups from June 2006 to December 2009. Sera were obtained from participants at 0, 6, 12, 24, and 36 months of age for measuring EV71 neutralizing antibody titers. If the participants developed suspected enterovirus illnesses, throat swabs were collected for virus isolation.

Results

We detected 28 EV71 infections including 20 symptomatic and 8 asymptomatic infections. Age-specific incidence rates of EV71 infection increased from 1.71 per 100 person-years at 0–6 months of age to 4.09, 5.74, and 4.97 per 100 person-years at 7–12, 13–24, and 25–36 months of age, respectively. Cumulative incidence rate was 15.15 per 100 persons by 36 months of age, respectively.

Conclusions

Risk of EV71 infections in Taiwan increased after 6 months of age during EV71 epidemics. The cumulative incidence rate was 15% by 36 months of age, and 29% of EV71 infections were asymptomatic in young children.     

The complete genome of Enterovirus 71 China strain

Five overlapping clones covering the full genome of Enterovirus 71 China strain SHZH98 were obtained and then the sequences were determined by the chain termination method. It showed that the full length of EV71 SHZH98 genome (not including Poly A tail) is 7408 bp. There are some diversities on the lengths and sequences of 5' UTR and 3' UTR between SHZH98 and the other EV71 strains. In P1 capsid region, which is closely associated with viral immunogenicity, EV71 strain SHZH98 shares the highest homology with Taiwan strains; but in P2 and P3 non-structural gene regions there are higher identities with Coxsakievirus A16 and EV71 strains MS, BrCr than with Taiwan strains. Phylogenetic tree constructed by structural gene region indicates that China strain SHZH98 has a closer relationship with Taiwan strains, however, in the non-coding region it has a closer relationship with Coxsakievirus A16, EV71 strains MS and BrCr. EV71 China strain was analyzed at the molecular level. The results will contribute to the     

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.     

The complete genome of Enterovirus 71 China strain

Five overlapping clones covering the full genome of Enterovirus 71 China strain SHZH98 were obtained and then the sequences were determined by the chain termination method. It showed that the full length of EV71 SHZH98 genome (not including Poly A tail) is 7408 bp. There are some diversities on the lengths and sequences of 5′ UTR and 3′ UTR between SHZH98 and the other EV71 strains. In P1 capsid region, which is closely associated with viral immunogenicity, EV71 strain SHZH98 shares the highest homology with Taiwan strains; but in P2 and P3 non-structural gene regions there are higher identities with Coxsakievirus A16 and EV71 strains MS, BrCr than with Taiwan strains. Phylogenetic tree constructed by structural gene region indicates that China strain SHZH98 has a closer relationship with Taiwan strains, however, in the non-coding region it has a closer relationship with Coxsakievirus A16, EV71 strains MS and BrCr. EV71 China strain was analyzed at the molecular level. The results will contribute to the basic study on enteroviruses and the EV71 prevention in China.     

The Association between Enterovirus 71 Infections and Meteorological Parameters in Taiwan

Background

Enterovirus 71 (EV71) infections are a significant cause of neurological disorder and death in children worldwide. Seasonal variations in EV71 infections have been recognized, but the mechanisms responsible for this phenomenon remain unknown. The purpose of this study was to examine the relationship between meteorological parameters and EV71 infection.

Methods and Findings

We analyzed the number of EV71 infections and daily climate data collected in Taiwan between 1998 and 2008 and used Poisson regression analysis and case-crossover methodology to evaluate the association between weather variability and the incidence of EV71 infection. A total of 1,914 EV71-infected patients were reported between 1998 and 2008. The incidence of EV71 infections reflected significant summertime seasonality (for oscillation, p<0.001). The incidence of EV71 infections began to rise at temperatures above 13°C (r² = 0.76, p<0.001); at temperatures higher than approximately 26°C (r² = 0.94, p<0.05), the incidence began to decline, producing an inverted V-shaped relationship. The increase in the incidence with increasing relative humidity was positive and linear (r² = 0.68, p<0.05). EV71 infection was most highly correlated with temperature and relative humidity in the period that likely preceded the infection.

Conclusion

Our study provides quantitative evidence that the rate of EV71 infection increased significantly with increasing mean temperature and relative humidity in Taiwan.     

Epidemiology and seroepidemiology of human enterovirus 71 among Thai populations

Background

Human enterovirus 71 (EV71) is an important pathogen caused large outbreaks in Asian-Pacific region with severe neurological complications and may lead to death in young children. Understanding of the etiological spectrum and epidemic changes of enterovirus and population’s immunity against EV71 are crucial for the implementation of future therapeutic and prophylactic intervention.

Results

A total of 1,182 patients who presented with the symptoms of hand foot and mouth disease (67.3%) or herpangina (HA) (16.7%) and admitted to the hospitals during 2008-2013 were tested for enterovirus using pan-enterovirus PCR targeting 5′-untranslated region and specific PCR for viral capsid protein 1 gene. Overall, 59.7% were pan-enterovirus positive comprising 9.1% EV71 and 31.2% coxsackievirus species A (CV-A) including 70.5% CV-A6, 27.6% CV-A16, 1.1% CV-A10, and 0.8% CV-A5. HFMD and HA occurred endemically during 2008-2011. The number of cases increased dramatically in June 2012 with the percentage of the recently emerged CV-A6 significantly rose to 28.4%. Co-circulation between different EV71 genotypes was observed during the outbreak. Total of 161 sera obtained from healthy individuals were tested for neutralizing antibodies (NAb) against EV71 subgenotype B5 (EV71-B5) using microneutralization assay. The seropositive rate of EV71-B5 was 65.8%. The age-adjusted seroprevalence for individuals was found to be lowest in children aged >6 months to 2 years (42.5%). The seropositive rate remained relatively low in preschool children aged > 2 years to 6 years (48.3%) and thereafter increased sharply to more than 80% in individuals aged > 6 years.

Conclusions

This study describes longitudinal data reflecting changing patterns of enterovirus prevalence over 6 years and demonstrates high seroprevalences of EV71-B5 NAb among Thai individuals. The rate of EV71 seropositive increased with age but without gender-specific significant difference. We identified that relative lower EV71 seropositive rate in early 2012 may demonstrate widely presented of EV71-B5 in the population before account for a large outbreak scale epidemic occurred in 2012 with due to a relatively high susceptibility of the younger population.     

Complete sequence analyses of enterovirus 71 strains from fatal and non-fatal cases of the hand,foot and mouth disease outbreak in Singapore (2000)

Enterovirus 71 (EV71) is a major aetiological agent of hand, foot and mouth disease (HFMD). In recent years, several outbreaks in East Asia were associated with neurological complications and numerous deaths. An outbreak in Singapore in October 2000 afflicted thousands of children, resulting in four fatal cases from three of whom EV71 was isolated. The genomes of two representative EV71 strains isolated from a fatal case and a surviving patient were completely sequenced, and their nucleotide and amino acid sequences compared with known EV71 strains. The two outbreak strains were classified under genogroup B, together with those previously isolated in Singapore, Malaysia and Japan. Comparative sequence analysis of the two Singapore strains revealed 99% nucleotide similarity, while their deduced amino acid sequences were almost identical except for residue 1506 in the 3A non-structural region. Given that the outbreak involved closely related genetic variants of EV71, the broad spectrum of disease severity may be attributed to critical factors such as varying viral inoculation doses or differing host immune responses following infection, but is less likely to be due to the emergence of EV71 strains with heightened virulence.     

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

Emerging epidemics of hand-foot-and-mouth disease (HFMD) associated with enterovirus 71 (EV71) has become a serious concern in mainland China. It caused 126 and 353 fatalities in 2008 and 2009, respectively. The epidemiologic and pathogenic data of the outbreak collected from national laboratory network and notifiable disease surveillance system. To understand the virological evolution of this emerging outbreak, 326 VP1 gene sequences of EV71 detected in China from 1987 to 2009 were collected for genetic analyses. Evidence from both traditional and molecular epidemiology confirmed that the recent HFMD outbreak was an emerging one caused by EV71 of subgenotype C4. This emerging HFMD outbreak is associated with EV71 of subgenotype C4, circulating persistently in mainland China since 1998, but not attributed to the importation of new genotype. Originating from 1992, subgenotype C4 has been the predominant genotype since 1998 in mainland China, with an evolutionary rate of 4.6∼4.8×10⁻³ nucleotide substitutions/site/year. The phylogenetic analysis revealed that the majority of the virus during this epidemic was the most recent descendant of subgenotype C4 (clade C4a). It suggests that the evolution might be one of the potential reasons for this native virus to cause the emerging outbreak in China. However, strong negative selective pressure on VP1 protein of EV71 suggested that immune escape might not be the evolving strategy of EV71, predicting a light future for vaccine development. Nonetheless, long-term antigenic and genetic surveillance is still necessary for further understanding.     

Characterization of an outbreak of hand, foot, and mouth disease in Nanchang, China in 2010

Recent outbreaks of human enterovirus 71 (EV71) infection and EV71-associated hand, foot, and mouth disease (HFMD) in China have affected millions and potentially lead to life-threatening complications in newborns. Furthermore, these outbreaks represent a significant global public health issue in the world. Understanding the epidemiology of HFMD and EV71 infection and their transmission patterns in China is essential for controlling outbreaks. However, no studies on the outbreaks of HFMD and EV71 infection in China during 2010 have been reported. In this report, we carried out an epidemiological analysis to study an outbreak of HFMD and EV71 infection in 2010 in the city of Nanchang in the Jiangxi province of People's Republic of China. From April 7 to May 11, 2010, a total of 109 HFMD cases were reported, and in this report the HFMD cases were studied by both epidemiological and laboratory analyses. The epidemiological study indicates that children aged younger than 8 years old represented more than 90% of the reported cases, with the age group of 1-3 years containing the highest number of cases. Laboratory studies detected a high prevalence of EV71 amongst the cases in our study, suggesting EV71 as a common enterovirus found in HFMD cases in Nanchang. Phylogenetic analysis of the sequence of the VP1 region of four EV71 isolates indicated that the Nanchang strains belong to the C4 subgenotype commonly found in China during outbreaks in 2008 but contain distinct variations from these strains. Our study for the first time characterizes the epidemiology of HFMD and EV71 infection in China in 2010 and furthermore, provides the first direct evidence of the genotype of EV71 circulating in Nanchang, China. Our study should facilitate the development of public health measures for the control and prevention of HFMD and EV71 infection in at-risk individuals in China.     

Phylogenetic Analysis of Enterovirus 71 Circulating in Beijing,China from 2007 to 2009

The major pathogens of hand, foot and mouth disease (HFMD) in Beijing, China from 2007 to 2009 were identified in this study. A total of 186 HFMD cases were included, and 136 cases (73%) were positive for enterovirus (EV). In 2007, 75% (27/36) were Coxsackievirus A16 (CA16) positive and 19% (7/36) were Enterovirus 71 (EV71) positive cases. However, EV71 was the predominant virus in 2008, when 56% (31/55) of the cases were positive for EV71 and 22% (12/55) were positive for CA16. In 2009, EV71 and CA16, with positive rates of 36% (16/45) and 29% (13/45), respectively, were still the major pathogens of HFMD. Phylogenetic analysis revealed that the dominant genotype of EV71 was C4, with co-circulation of genotype A in 2009. The prevalent cluster of the EV71 subgenotype C4 changed over time. A proposed new sublineage of EV71, C4a-2, was the predominant virus associated with the Beijing and nationwide HFMD outbreaks since 2008 and amino acid substitution, which possibly link to the central nervous system tropism of EV71, was found in genotype A viruses. Persistent surveillance of HFMD-associated pathogens is required for predicting potential emerging viruses and related disease outbreaks.     

Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998

Enterovirus 71 (EV71) (genus Enterovirus, family Picornaviridae), a common cause of hand, foot, and mouth disease (HFMD), may also cause severe neurological diseases, such as encephalitis and poliomyelitis-like paralysis. To examine the genetic diversity and rate of evolution of EV71, we have determined and analyzed complete VP1 sequences (891 nucleotides) for 113 EV71 strains isolated in the United States and five other countries from 1970 to 1998. Nucleotide sequence comparisons demonstrated three distinct EV71 genotypes, designated A, B, and C. The genetic variation within genotypes (12% or fewer nucleotide differences) was less than the variation between genotypes (16.5 to 19.7%). Strains of all three genotypes were at least 94% identical to one another in deduced amino acid sequence. The EV71 prototype strain, BrCr-CA-70, isolated in California in 1970, is the sole member of genotype A. Strains isolated in the United States and Australia during the period from 1972 to 1988, a 1994 Colombian isolate, and isolates from a large HFMD outbreak in Malaysia in 1997 are all members of genotype B. Although strains of genotype B continue to circulate in other parts of the world, none have been isolated in the United States since 1988. Genotype C contains strains isolated in 1985 or later in the United States, Canada, Australia, and the Republic of China. The annual rate of evolution within both the B and C genotypes was estimated to be approximately 1.35 x 10(-2) substitutions per nucleotide and is similar to the rate observed for poliovirus. The results indicate that EV71 is a genetically diverse, rapidly evolving virus. Its worldwide circulation and potential to cause severe disease underscore the need for additional surveillance and improved methods to identify EV71 in human disease.     

Characterization of full-length enterovirus 71 strains from severe and mild disease patients in northeastern China

Human enterovirus 71 (EV71)-associated hand, foot, and mouth disease (HFMD) has been a leading cause of childhood infection in China since 2008. Epidemic and molecular characteristics of HFMD have been examined in many areas of China, including the central and southern regions. However, clinical and genetic characterization of EV71 in the northeastern region of China is scarce. In this study, a series of analyses were performed on seven full-length EV71 sequences from HFMD patients who had either severe or mild disease. We have determined that these seven circulating EV71 viruses from Changchun, China are actually complex recombinant viruses involving multiple type A human enterovirus (HEV). Classified as EV71 subtype C4 (EV71 C4), these Changchun EV71 viruses contain genetic recombination events between the CA4, CA5, EV71B4 and EV71C1 strains. Most of the structural protein region (P1) of these viruses resembled that of the prototype EV71 C1 strains. The non-structural protein domains (P2 and P3) showed a high degree of similarity with CA4, CA5 and EV71 B4 in different regions. The 5'UTR had unclassified recombination,while partial 3D region of these viruses showed a high degree of similarity to CA16. Phylogenetic analysis of full-length or partial sequences of isolates from severe or mild disease patients in Changchun always formed a single cluster in various phylogenetic analyses of different genomic regions, suggesting that all seven strains originated from one single common ancestor. There was no correlation between viral genomic sequence and virulence. Thus, we found that circulating recombinant forms of EV71 are prevalent among HFMD patients in Northeastern China. The existence of a unique cluster of EV71 related viruses in Northeast China has important implications for vaccine development that would address the increasing prevalence of HFMD.     

An attenuated strain of enterovirus 71 belonging to genotype a showed a broad spectrum of antigenicity with attenuated neurovirulence in cynomolgus monkeys

Enterovirus 71 (EV71) is a causative agent of hand, foot, and mouth disease and is also sometimes associated with serious neurological disorders. In this study, we characterized the antigenicity and tissue specificity of an attenuated strain of EV71 [EV71(S1-3')], which belongs to genotype A, in a monkey infection model. Three cynomolgus monkeys were inoculated with EV71(S1-3'), followed by lethal challenge with the parental virulent strain EV71(BrCr-TR) via an intravenous route on day 45 postinoculation of EV71(S1-3'). Monkeys inoculated with EV71(S1-3') showed a mild neurological symptom (tremor) but survived lethal challenge by virulent EV71(BrCr-TR) without exacerbation of the symptom. The immunized monkey sera showed a broad spectrum of neutralizing activity against different genotypes of EV71, including genotypes A, B1, B4, C2, and C4. For the strains examined, the sera showed the highest neutralization activity against the homotype (genotype A) and the lowest neutralization activity against genotype C2. The order of decreasing neutralization activity of sera was as follows: A > B1 > C4 > B4 > C2. To examine the tissue specificity of EV71(S1-3'), two monkeys were intravenously inoculated with EV71(S1-3'), followed by examination of virus distribution in the central nervous system (CNS) and extraneural tissues. In the CNS, EV71(S1-3') was isolated only from the spinal cord. These results indicate that EV71(S1-3') acts as an effective antigen, although this attenuated strain was still neurotropic when inoculated via the intravenous route.     

Statistical tests for natural selection on regulatory regions based on the strength of transcription factor binding sites

Background  

Although cis -regulatory changes play an important role in evolution, it remains difficult to establish the contribution of natural selection to regulatory differences between species. For protein coding regions, powerful tests of natural selection have been developed based on comparisons of synonymous and non-synonymous substitutions, and analogous tests for regulatory regions would be of great utility.     

Seroepidemiology and Molecular Epidemiology of Enterovirus 71 in Russia

Enterovirus 71 (EV71) is an emerging human pathogen causing massive epidemics of hand, foot and mouth disease with severe neurological complications in Asia. EV71 also circulates in Europe, however it does not cause large outbreaks. The reason for distinct epidemiological patterns of EV71 infection in Europe and Asia and the risk of EV71 epidemic in Europe and Russia remain unknown. Seroepidemiology of EV71 and molecular epidemiology of occasional EV71 isolates were studied to explore circulation of EV71 in Russia. In six regions of Russian Federation, seroprevalence of EV71 in sera collected in 2008 ranged from 5% to 20% in children aged 1–2 years and from 19% to 83% in children aged 3–5 years. The seroprevalence among elder children was significantly higher (41–83% vs. 19–27%) in Asian regions of Russia. EV71 strains identified in Russia in 2001–2011 belonged to subtypes C1 and C2, while genotype C4 that was causing epidemics in Asia since 1998 emerged in 2009 and became dominant in 2013.     

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.