Journal of ISSN: 2373-6453JHVRV

Human Virology & Retrovirology
Review Article
Volume 3 Issue 2 - 2016
Molecular Analysis of Enteroviruses
Fei Zhou*, Matthew VN O’Sullivan, Jonathan R Iredell, Vitali Sintchenko, Gwendolyn L Gilbert and Dominic E Dwyer
Centre for Infectious Diseases and Microbiology (CIDM), Institute of Clinical Pathology and Medical Research (ICPMR), Westmead Hospital, The University of Sydney, Australia
Received:February 16, 2016 | Published: February 26, 2016
*Corresponding author: Fei Zhou, Centre for Infectious Diseases and Microbiology (CIDM), Institute of Clinical Pathology and Medical Research (ICPMR), Westmead Hospital, Darcy Road, Westmead, New South Wales, Australia, 2145, Tel: (612) 9845 6255; Fax: (612) 9893 8659; Email:
Citation: Zhou F, O’Sullivan MVN, Iredell JR, Sintchenko V, Gilbert GL, Dwyer DE (2016) Molecular Analysis of Enteroviruses. J Hum Virol Retrovirol 3(2): 00083. DOI: 10.15406/jhvrv.2016.03.00083


Enteroviruses (EVs) are common human pathogens which usually infect the gastrointestinal and respiratory tracts and can spread to other organs or systems. They are characterized by high genomic plasticity, primarily due to high mutation and recombination rates. Improved molecular diagnostic methods and genetic sequence analyses are beginning to discover the complex characteristics of individual serotypes and genotypes. Understanding the tempo and pattern of molecular diversity and evolution is of great importance in the pathogenesis of EVs, information which will assist in disease prevention and control. In this review, we will focus on molecular analysis of EVs, including current diagnosis, epidemiology and evolution.

Keywords: Enterovirus; Molecular diagnosis; Molecular epidemiology; Molecular diversity; Evolution; EVs


PCR: Polymerase Chain Reaction; IRES: Internal Ribosome Entry Site; UTRs: Untranslated Regions; RT: Reverse Transcription; EVs: Enterovirus


The order Picornavirales consists of the families Picornaviridae, Dicistroviridae, Iflaviridae, Marnaviridae and Secoviridae and two unassigned genera, Bacillarnavirus and Labyrnavirus. The family Picornaviridae currently consists of 50 species grouped into 29 genera ( Enteroviruses (EVs) (genus Enterovirus, family Picornaviridae) are common human pathogens associated with a wide spectrum of conditions, ranging from asymptomatic infection to serious illness, such as aseptic meningitis, meningoencephalitis, myocarditis and acute flaccid paralysis, especially in children [1] and immunocompromised patients [2]. The original classification of EVs consisted of polioviruses (PVs), coxsackieviruses A (CV-As) or B (CV-Bs) and echoviruses (Es), based on biological activity and disease. Currently, the genus EV is classified to 12 species, EV A to H, J and Rhinovirus A to C ( Seven of the species, EV A to D (formerly named Human enterovirusA to D) and Rhinovirus A to C (formerly named Human rhinovirus A to C) are known to infect humans [3].

The EV genome is a positive single-stranded RNA molecule of approximately 7,500 nucleotides, comprising a single open reading frame flanked 5’ and 3’ by untranslated regions (UTRs) (Figure 1). The cloverleaf structure of domain (stem-loop) I in the 5’ UTR is important for viral replication [4]; while domains II to VI encompass the internal ribosome entry site (IRES) which directs translation of the mRNA by internal ribosome binding [5]. The coding region, divided into three subregions (P1, P2 and P3), consists of a single open reading frame encoding a polyprotein. The P1 region encodes four structural capsid proteins (VP4, VP2, VP3 and VP1). VP1, VP2 and VP3 are located at the surface of the viral capsid and are exposed to immune pressure, whereas VP4 is located inside the capsid.

Figure 1: Structure of the enterovirus RNA genome, with the genome-linked protein VPg at the 5’ end, the 5’ UTR, the protein coding region (P1, P2 and P3), the 3’ UTR and the poly (A) tail. Coding regions for the viral proteins are indicated.

The VP1 capsid protein is the most external and immunodominant of the picornavirus capsid proteins [6], and contains most neutralization epitopes. The P1 coding region for the capsid proteins provides reliable correlation between sequence relatedness and the previous definition of serotype [7]. This also appears to be true for the various individual capsid protein regions, except for VP4, whose sequence does not always correlate with serotype [8,9]. Recent molecular studies suggest that the VP1 nucleotide sequence correlates well with antigenic typing by neutralization and can be used for virus identification and evolutionary studies [7,10,11]. VP1 sequencing has supported and clarified early serologic data, and assisted in proposals of new EV types [7,12,13]. In addition, the non-structural proteins are encoded in the P2 (2A, 2B and 2C) and P3 (3A, 3B, 3C and 3D polymerase) regions.

Although the EV antigenic properties are very important, the introduction of molecular typing methods and a reassessment of the limitations of the original classification (PV, CV-A or CV-B and echoviruses) have resulted in the development of the current classification scheme. Currently, the genus EV is classified to 12 species, EV A to H, J and Rhinovirus A to C ( A picornavirus species is a class of phylogenetically related serotypes or strains which would normally be expected to share

  1. A limited range of hosts and cellular receptors,
  2. A significant degree of compatibility in proteolytic processing, replication, encapsidation and genetic recombination,
  3. Essentially identical genome maps ( EV-A, -B, -C and -D species consist of 25, 63, 23 and five (sero) types, respectively (Table 1)


In this review, we will focus on molecular analysis of EVs, including current diagnosis, epidemiology and evolution.


(Sero) Types


CV-A2 to -A8, -A10, -A12, -A14, -A16, EV-A71, 9 new types (EV-A76, -A89, -A90, -A91, -A92, -A114, -A119, -A120 and -A121) and the simian EVs  SV19, SV43, SV46 and baboon EV A13 (BA13)


CV-B1 to -B5 (including swine vesicular disease virus [SVDV]), -B6, -A9, 28 echovirus serotypes (E-1, -2, -3, -4, -5, -6, -7, -9, -11, -12, -13, -14, -15, -16, -17, -18, -19, -20, -21, -24, -25, -26, -27, -29, -30, -31, -32 and -33), EV-B69, 26 new types (EV-B73,  -B74, -B75, -B77, -B78, -B79, -B80, -B81, -B82, -B83, -B84, -B85, -B86, -B87, -B88, -B93, -B97, -B98, -B100, -B101, -B106, -B107, -B110, -B111, -B112 and -B113) and the simian enterovirus SA5


PV-1 to -3, CV-A1, -A11, -A13, -A17, -A19, A20, -A21, -A22, -A24 and 11 new types (EV-C95, -C96, -C99, -C102, -C104, -C105, -C109, -C113, -C116, -C117 and -C118)


EV-D68, -D70, -D94 and 2 new types (EV-D111 and -D120)

Table 1: EV-A, B, C and D species and their corresponding (sero) types a.


Molecular diagnosis of enteroviral infections

The traditional approaches to detecting and characterizing EVs are based on the time-consuming and labour-intensive procedures of viral isolation in cell culture and neutralization with mixed hyperimmune equine serum pools and specific monovalent polyclonal antisera for confirmation [14]. Due to distinct advantages in speed and convenience, molecular biology techniques are supplanting the traditional approaches of EV detection and characterization [10,14]. Continued development and optimization of these methods will bring about further application of these techniques for EV diagnostics. As with virus isolation and serotyping, molecular approaches can detect the presence of EV in a specimen or confirm an isolate, and with some procedures, to further characterize the detected virus [15,16].

The first category of nucleic acid tests is based on the polymerase chain reaction (PCR). This is applied primarily to detect EV genomes in clinical specimens, cell cultures and biopsy or autopsy tissues. The design of EV-specific primers (and probes) in the 5’ UTR to target all known EV members is probably practical and reasonable [14,16-19]. Therefore, 5’ UTR EV-specific primers (and probes) can be used as an effective tool to confirm whether an isolate is an EV, and the design of generic assays lowers the risk of missing newly emerging or still unrecognized variants [20,21]. Currently, PCR is commonly used to directly detect viruses in clinical specimens from individuals suspected of EV infection. The major advantage of the pan-EV PCR is that rapid EV detection is possible, even with very small amounts of clinical material such as cerebrospinal fluid [22]. However, the sensitivity may be lower in certain specimens (e.g., stool) [23]. Furthermore, by changing the target, PCR can be used to characterize a particular EV. Certain genomic sequences within the capsid-coding region (especially VP1) correspond to the conserved antigenic property of the virus, since the antigenic property that defines serotype is a property of the viral capsid proteins [7]. Investigations of the PV capsid regions have led to primer design that can selectively amplify isolates from a single serotype, but not from other heterologous serotypes [24]. These principles have been extended to other EV serotypes or receptor groups, providing efficient methods for rapid EV characterization, such as type-specific reverse transcription (RT)-qPCR for sensitive detection of E30 in cerebrospinal fluid specimens [25] and a real-time RT-PCR assay for the detection of EV-D68 in clinical specimens [26].

The second approach utilizes nucleic acid probes, and was the first application of molecular techniques to detect EV genomes. Early probes were made by synthesizing cDNA from CV-B viral RNA. Then, the cDNA was cloned into bacteria, amplified, and radiochemically labelled to provide a reagent for detecting EV genomes by nucleic acid hybridization [27]. Further refinement of this technique included the use of RNA probes to improve sensitivity, the introduction of nonradiochemical probe labelling, and the application of these techniques to in situ hybridization in myocardial tissue [28,29]. Application of bioinformatics in probe design enables detection of EVs from formalin-fixed tissue samples on different taxonomic levels by advanced in situ hybridization technology [30]. Nucleic acid probe hybridization has also been broadly used in the WHO PV eradication program to determine whether an isolate is a vaccine-derived or a wild-type PV [31,32]. Microarray techniques for the detection and characterization of EVs have been developed and applied in some studies, but they can only identify limited serotypes [33,34]. In addition, a modified RT-PCR-based reverse line blot hybridization assay for synchronous detection of 5’ UTR and VP1 region has been developed as a quick, accurate and efficient approach to characterizing common and nonserotypeable EVs [11,14].

The ease and low cost of this strategy will improve the diagnosis and epidemiologic investigation of EV infections and outbreaks where other typing methods are unavailable. The third group of molecular methods is based on EV genomic sequencing for characterization at the highest levels of specificity and accuracy [10,16,35]. The information from the nucleotide sequence of a given virus represents its fundamental characterization [36,37]. It is practical to use sequence information to assign an EV isolate to a particular (sero) type [37]. The molecular typing system is based on RT-PCR and nucleotide sequencing of the 3’ half or the entire VP1 genomic region [11,35], and these can function as an excellent surrogate for serotyping by neutralization. The (sero)type of an unknown isolate is inferred by comparison of the partial or complete VP1 sequence with a database containing VP1 sequences for the prototype and variant strains of all known EV (sero)types. The principles of EV assignation include:

  1. A partial or complete VP1 nucleotide sequence identity of ≥75% (>85% or 88% amino acid sequence identity) between a clinical EV isolate and a serotype prototype strain may be used to establish the serotype of the isolate, on the provision that the second highest score is <70%;
  2. A best-match nucleotide sequence identity of <70% may indicate that the isolate represents an unknown or new serotype;
  3. A sequence identity between 70% and 75% indicates that further characterization is required before the isolate can be identified firmly [7,38]. Using these criteria, strains of homologous (sero) types can be easily distinguished from heterologous (sero) types. This approach can greatly reduce the time for typing EV isolates, and can be used to type isolates that are difficult or impossible to type using conventional serotyping methods. Molecular identification and analysis of nonserotypeable EV isolates is relevant to epidemiological surveillance, to determine relationships between types and diseases, and to assist in discovery of new types [14]. Furthermore, whole genome sequencing technology has been developed recently, which can enhance the diagnosis of EVs and reveal epidemiological and evolutionary dynamics of EVs [1,39].

Other molecular methods that assist serotype-specific identification include combined multiplex RT-PCR and microarray to detect and differentiate EV-A71 and CV-A16 [33], molecular characterization of PV strains by RT-PCR-restriction fragment length polymorphisms [40], an integrated micro-RT-PCR system for automatic detection of microorganisms including EV-A71 [41], and an RT-loop-mediated isothermal amplification system for rapid and highly sensitive detection of EVs in stool samples [42].

Molecular epidemiology of enteroviruses

Methods of molecular epidemiology: EVs can be transmitted by both faecal-oral and respiratory routes. Studies of the variations of enteroviral proteins or nucleic acids have provided significant epidemiologic information [43,44]. Molecular epidemiologic investigations have aided in our understanding of EVs by:

  1. Providing the opportunity for specific strain identification [45];
  2. Providing insights into EV classification and taxonomy [46];
  3. Clarifying outbreak origins [47];
  4. Identifying distinct strains transmitted between outbreaks [48];
  5. Determining the evolutionary genetics of specific (sero) types [49].

Molecular typing approaches (including PCR and genomic sequencing) allow rapid and accurate identification of EVs. The major method applied to obtain epidemiologic information is direct analysis of genomic variations of multiple sequences (sometimes local sequences compared with GenBank database sequences) using different types of software (such as MEGA and BioNumerics) following nucleic acid sequencing in different regions of EVs [18,50]. Phylogenetic analysis of sequences in different genes (such as 5’ UTR and VP1) can be used to discriminate between variants within a (sero) type, to enhance the epidemiologic linkage of EV strains, and to study the evolution of a given (sero) type or among different (sero) types [17-19]. Analysis of whole genome sequences has been used to reveal epidemiologic and evolutionary dynamics of EVs through time and space [1,39,51]. The combination of molecular typing and phylogenetic sequence analysis will assist with both individual patient diagnosis and public health measures. Moreover, some novel methods using bioinformatic technology have appeared. For example, the combination of multiple sequence alignment with the linear mapping hash function is a computationally efficient way of gene sequence clustering and can be a valuable tool for the assessment of similarity, clustering of different microbial genomes (such as VP1 region of EV-A71), identifying reference sequences, and for the study of evolution of bacteria and viruses [53].

Molecular epidemiology of polioviruses

The application of nucleic acid sequencing to the studies of wild-type PV isolates from different areas has remarkably expanded the power of epidemiologic studies [53,54]. Importantly, phylogenetic trees or lineage maps can be used for demonstration of sequence relationships [17,55,57]. Closely related viruses can be easily detected by analysing random mutations in different PV genomes; while more distantly related viruses can be clustered into distinct geographic groups of endemic circulation. From 2001, the WHO polio network laboratories routinely sequence the complete VP1 gene of all wild-type and vaccine-derived PV isolates for comparison. Using VP1 sequence analysis, the genetic diversity of PV strains can be exploited in molecular surveillance [55].

Moreover, the rapid evolution of PVs permits the use of comparative VP1 nucleotide sequencing to resolve the fine structure of a poliomyelitis outbreak, and to elucidate individual chains of transmission [58]. Based on a nucleic acid sequence database of PV strains worldwide, it has been possible to develop rapid approaches to tracking wild-type PV strains [59]. Recently, in response to the observation that the large majority of vaccine-derived PV isolates are type 2, WHO has developed a plan for coordinated worldwide replacement of trivalent oral PV vaccines (tOPV) with bivalent OPV (bOPV; types 1 and 3) in April 2016, preceded by the introduction of at least one dose of inac­tivated PV vaccine into routine immuniza­tion schedules worldwide [60].

Molecular epidemiology of non-polio enteroviruses

Research on the molecular epidemiology of non-polio EVs has concentrated on the evolutionary features following the comparison of viral isolates within a (sero)type over time [19,61], the comparison of isolates from different (sero)types [18,62], and between different genera within the Picornaviridae[63].

EV-A71 and CV-A16 are members of EV-A species, and are major and independent etiologic agents of hand-foot-mouth disease [64,65]. To understand the molecular epidemiology and evolution of EV-A71, whole-genome sequences [1] and several distinct regions of the EV-A71 genome have been employed for analysis, such as 5’ UTR [19,66], 3D polymerase region [67], VP4 region [68], and VP1 region [19,49,61,69,70]. However, the VP1 gene is considered to be the most informative and robust region for evolutionary study due to a high degree of diversity and lack of involvement in recombination [71]. EV-A71 strains isolated worldwide are classified into four genogroups: A-D [72].

Genogroups B and C have been differentiated into subgenogroups B0-B5 and C1-C5 [72]. Recent studies have proposed three new genogroups, including E (in sub-Saharan Africa), F (in Madagascar) [73] and G (in India) [74]. Other new subgenogroups proposed recently include C4a, C4b, C6, C7 and B6 [75,76]. The analysis of EV-A71 based on VP1 region, in combination with 5’ UTR [19] or 3D polymerase region [67], could provide more objective information on molecular epidemiology and evolution of EV-A71. For CV-A16, the VP1 and VP4 genes are mostly used for classifying the virus [19,77-79]. It has been classified on VP1 analysis into geno groups A and B, with geno group B further divided into subgenogroups B1 (including B1a, B1b and B1c), B2 [64,80] and a possible B3 [19]. Molecular analysis of VP1, combined with 5’ UTR, may also enhance the reliability of molecular surveillance and the understanding of genetic evolution of CVA16 [19].

In EV-B species, several CV-B outbreaks or epidemics have been recorded. CV-B1 infection accounted for 23% of reported cases with severe neonatal morbidity and multiple deaths between 2007 and 2008 in the USA [81]. CV-B5 has caused widespread epidemics in the USA, particularly in 1967, 1972 and 1983 [82]. The sequence analyses of these CV-B viruses have been based on VP1 region, 5’ UTR and 3D [18,81-84]. Four main genotypes (GI-GIV) and four small clusters (C1-C4) of CV-B1 have been described based on the VP1 region [84]. E11 causes a full spectrum of enteroviral diseases, including meningitis [85], but it has been particularly associated with cases of severe neonatal disease [86]. A few studies have analysed the genetic diversity among E11 VP1 (-2A), 3D polymerase and 5’ UTR sequences [18,85-87], which divided E11 strains into different subgroups. The E11-B cluster has been shown to represent a phylogenetically distinct subgroup (with altered virulence) within the E11 serotype [87]. Since its characterization in 1958, E30 has been increasingly associated with meningitis outbreaks worldwide [88]. Several studies have shown different clustering of E30 variants based on VP1 and/or 5’ UTR sequences [18,89-91].

Investigations have shown that EV-C (sero) types (including PVs) cluster into three subgroups in the VP1 region, designated as A to C [92,93]. CV-A11 groups together with the strains of subgroups B or C depending on the method of phylogenetic tree construction [93]. The major worldwide CV-A24 epidemics have been attributed to four genotypes of CV-A24 variant strains (I to IV) on the basis of phylogenetic analysis of the 3C protease or VP1 region [94,95].

After the initial identification of EV-D68 (belonging to EV-D species) in 1962, detection of this virus was rarely reported until the early 2000s. EV-D68 has caused widespread respiratory illness recently. Neurological syndromes including acute flaccid paralysis following EV-D68 infection have also been reported in a small number of cases [39,96]. Whole-genome sequences [39] and several distinct regions of the EV-D68 genome have been used for analysis of molecular epidemiology, such as 5’ UTR, VP1 and VP4/ VP2 regions [97,98]. The EV-D68 strains circulating in recent years are classified into three genetic groups [99,100]. For the VP1 region, strains belonging to lineages 1 and 2 are likely to form clusters containing viruses from the same geographical origins on the phylogenetic tree. In contrast, strains belonging to lineage 3 are less likely to show geographical clustering patterns. These data suggest that transmission at the global level might be more common for lineage 3 than for lineages 1 and 2 [101].

Molecular diversity and evolution of enteroviruses

Most RNA virus genomes have much higher evolution rates than DNA viruses [102]. EVs are well known for high genomic plasticity due to high mutation and recombination rates. Analysis of whole genome sequences has been used to reveal epidemiologic and evolutionary dynamics of EVs through time and space [1,39,51]. Studies of EVs suggested that sequence variations in the UTRs, the VP1 region and the 3D polymerase could affect virus infection and replication capability in vitro and in vivo [103,104].

Mutation of enteroviruses

Base misincorporation during chain elongation and lack of 3’ to 5’ exonuclease proofreading ability in RNA polymerases can cause very high error rates in EVs; there is a rate of spontaneous mutation of approximately one mutation per genome per replication [105]. High RNA virus error rates are necessary to enable survival of the virus population under selective pressure [106]. The molecular diversity and evolution of EVs occurs through genetic drift and, over much longer periods, antigenic diversification in the structural gene region encoding the virus capsid, including VP1 [61,107]. Sequence changes occur rapidly in the capsid region [61]. Studies have also proved that mutations in the 5’ UTR markedly decrease multiplication efficiency [108], alter cell tropism [109], and attenuate virulence [110].

In EV-A species, the overall rapid evolutionary rate of the EV-A71 VP1 gene has been estimated to be 4.2 × 103 and 3.4 × 10-3 [61] or 4.5 -4.6 × 10-3 and 4.2 × 10-3 [49] all substitutions per nucleotide per year for genogroups B and C, respectively. The Ala to Val substitution at VP1 position 170 appears to be associated with increased neurovirulence [69]. Additionally, the nucleotide mutation (A→G) at position 485 in the 5’ UTR may be one of the molecular determinants of EV-A71 attenuation of neurovirulence in cynomolgus monkeys [111]. The average evolutionary rate calculated for CV-A16 B1a and B2 was 0.91 × 10-2 synonymous substitutions per nucleotide per year based on VP1 gene [64]. The evolutionary rate of CV-A16 is relatively slow given that the prototype CV-A16 strain was first identified in South Africa in 1951, and that all other CV-A16 strains formed a single genogroup B following approximately 60 years of evolution [64]. Comparing two regions (5’ UTR and VP1) of EV-A71 and CV-A16 viruses provides evidence of epidemiologic linkage of EV-A strains, and mutation in the two regions plays a vital role in the evolution of these viruses [19].

In EV-B species, the broad spectrum of CV-B associated disease reflects the existence of multiple strains within a single serotype with various degrees of virulence [112]. Individual nucleotide substitutions in the noncoding and coding regions of the viral genome determine virulence [113]. Major pathogenic determinants of CV infections have been localized to the 5’ UTR and the capsid protein coding region [114]. The estimated rate of CV-B1 evolution in the VP1 gene was 7.73 × 10-3 substitutions per site per year [84]. Strong negative selection shaped the evolution of VP1 and 3CD loci in CVB5 strains, but compelling phylogenetic data suggested that immune selection pressure resulted in the emergence of the two genogroups (A and B) with opposed evolutionary pathways [115]. The genogroups A and B also differed in the temporal occurrence of the amino acid changes. Studies of E30 meningitis outbreaks have shown that point mutations give rise to substantial genetic diversity in VP1. Substitutions in VP1 region occurred predominantly at synonymous sites, with VP1 showing a rapid substitution rate of 8.3/8.4 × 10-3 substitutions per site per year, and recombination frequency was tightly correlated with VP1 divergence [116,117]. The predicted rates of evolution in the VP1 region for E9 and E11 were 5.8 × 10-3 and 4.8 × 10-3 substitutions per site per year, respectively [118]. Moreover, the study of Spanish E6 strains revealed high rates of sequence change of the VP1 region (around 1.1× 10-2 substitutions per site per year) and there was a direct relationship between recombination frequency and VP1 sequence divergence for E6 [119]. The relatively high sequence diversity among the 5’ UTR and VP1 regions of serotypes in EV-B species suggests that nucleotide substitution is probably the dominant evolutionary mechanism in these regions of the genome [18].

In EV-C species, PVs are among the most rapidly evolving viruses known, with most evolution appearing to be random genetic drift since >80% of nucleotide substitutions within the coding region generate synonymous codons. Nucleotide substitutions accumulate at an overall rate of approximately 1% per year at all sites, and at approximately 3% per year at synonymous sites, calculated primarily on the VP1 region or P1/capsid region [106,120,121]. Evolution rates appear to be similar across serotypes and between wild-type and vaccine-derived PVs [56,106]. A common feature of the Sabin oral PV vaccine strains is the presence of nucleotide substitutions in the IRES of 5’ UTR, which in serotypes 1 and 3 are critical attenuating mutations. Additional mutations encoding amino acid substitutions in the capsid region (including VP1) contribute to and stabilize the attenuated phenotype [122]. In CV-A24 variant strains, the estimated evolutionary rates were 1.0 - 4.1 × 10-3 or 6.27 - 6.67 × 10-4 per nucleotide per year based on the 3C protease gene [94,123-125], 4.2 × 10-3 per nucleotide per year based on the VP4 gene [123], and 1.83-2.17 × 10-3 per nucleotide per year based on the 3’ VP1 gene [125].

The mean nucleotide substitution rates of EV-C subgroup B (ranging from 1.170 to 3.625 × 10-3 with combined high-probability distribution interval of 1.170 to 3.625 for partial VP1) are slightly lower than the corresponding estimates for other EV serotypes [93]. Intraserotypic genetic change of VP1 is dominated by silent mutations accompanied by amino acid polymorphism occurring dominantly at immunogenic sites, while interserotypic differences of VP1 include permanent fixation and insertion/deletion of distinct ‘signature’-amino acids that could be a result of larger scale changes in the capsid structure [93]. The relatively high sequence diversity among the 5’ UTR and VP1 regions of EV-C viruses also suggests that nucleotide substitution is probably the dominant evolutionary mechanism in these regions of the genome [17].

In EV-D species, the mean rate of nucleotide substitution for VP1 gene of EV-D68 was 6.2 × 10-3 substitutions per site per year [97]. Mutations have accumulated in two specific regions of the EV-D68 VP1 sequences: the BC and DE loops [100]. Selection analyses conducted on the VP1 amino acid sequences in several studies have revealed that most of the positively selected codons are found in the BC and DE loop regions [126,127]. These findings suggest that antigenic epitopes might be located in the BC and DE loop regions. Unique amino acid mutations in the VP2 and VP3 regions that can differentiate recently detected strains from the Fermon strain, as well as recent strains of one genetic lineage from another, have been identified; however, such mutations were found less frequently than those in the VP1 region [127]. Strains of lineage 1 were likely to have deletions at both positions 681-704 and 721-731 in the 5′UTR. In contrast, lineage 3 strains had deletions only at positions 624-704, whereas those of lineage 2 had no deletions in the 5′UTR [101].

Recombination of enteroviruses

The evolutionary diversity of EVs is also due to intraserotypic / interserotypic (or interspecies) recombination, the majority of which occurs within P2, P3 non-structural coding regions and 5’ UTR, with little or no recombination within the P1 capsid coding region [85,128-130]; however, the survival of offspring with recombination within the P1 capsid coding region is rare [8,128]. In general, only the P1 capsid region (mainly VP2-VP3-VP1) is inherited as a single unit, suggesting that it is the primary determinant of EV identity [129-131]. Due to recombination, no general correlation is found between sequence relatedness and specific serotype in genome segments outside the capsid region.

RNA recombination is believed to be coupled with the process of genome RNA replication. It occurs by template switching during negative strand synthesis, as first demonstrated in PV-infected cells [132] and subsequently in cell-free extracts [133]. Recombination can occur between PVs and other EV-C viruses within the 5’ UTR and non-structural coding regions; this shuffling of different non-structural coding regions may lead to serotypes with selective advantages [3,46,120]. Recombination within the 5’ UTR and the non-structural coding regions was also found among EV-A and EV-B members (within the same species) [129;134]. The frequency of recombination differs between EV species, being higher in EV-B species than in EV-A species [8] and very low in EV-D species [135,136].

In EV-A species, EV-A71 intraserotypic and interserotypic recombination have been reported [19,67,137]. Recombination would be facilitated by the presence of several subgenogroups circulating during the same period in one country or region, such as the co-circulation of four subgenogroups B3, B4, C1 and C2 in Malaysia between 1997 and 2000 [138]. Recent studies revealed that some CV-A16 isolates might have recombination with the EV-A71 isolates in the 5’ UTR or 2A-2B region, and such recombination events might change the epidemic trend [80,139].

In EV-B species, although the major neurotropism determinant in the CV-B3 genome was in the capsid region, viruses containing the CV-A9 capsid were also able to initiate infection in the central nervous system provided they contained the CV-B3 5’ UTR [114]. The presence of the 5’ UTR of CV-A9 clearly enhanced muscle tissue tropism. The viability of chimeric viruses constructed between CV-A9 and CV-B3 suggested that recombination could take place in nature, partly explaining the genetic diversity and pathogenic variation among CVs [114]. Studies have demonstrated that recombination can occur between the 5’ UTR and capsid, the capsid and 3D, or both in some CV-B cases [18,84,140].

Recombination could potentially lead to new strains that might cause outbreaks of disease. Furthermore, studies of echoviruses suggest that recombination has played a significant role in the evolution of these echoviruses. For example, E30 can recombine with any EV-B strain in the non-structural region or 5’ UTR [18.90,141]. In the non-structural coding region of E9, E11 and E30, sequence change has been shown by frequent recombination events, leading to the acquisition of novel 3D polymerase sequences interspersed with sequences of other EV-B isolates [118]. The first evidence for natural intertypic recombination for E24 was described recently, which showed that this E24 strain recombined with a CV-B1 stain at the 2C region [51]. In fact, interserotypic recombination may play a major role in the evolution of some EV-B viruses [18]. Due to the need of gene regulation, restricted genetic variation has been noted mainly in the non-structural coding region, and in some cases recombination occurs in the 5’ UTR.

In EV-C species, natural interserotypic recombination within the capsid region of PVs appears to be rare, probably because structural incompatibilities restrict the fitness of most interserotypic recombinants [142]. In contrast, interserotypic recombinants excreted by recent vaccinees exposed to tOPV and recombinants among circulating wild-type PVs and vaccine-derived PVs usually have crossovers restricted to the noncapsid region [142]. PVs recombine actively with each other and with other closely related EV-C members, a process which might offset the effects of accumulation of deleterious mutations arising in some lineages [46,121,142,143]. In fact, the recombination between PVs and other EV-C viruses is a normal feature of PV evolution in nature. Recombination may lead to a lack of linkage between noncapsid sequences and the neutralisation serotype of a virus within a species [46].

In addition, the prototype EV-C strains are recombinant relative to each other [92], and circulating EVs commonly recombine with wild-type and vaccine PVs [144]. Studies have shown that, compared with VP1, the 5’ UTR of CV-A (sero)types and PVs has frequently been subject to recombination [17,145], indicating that interserotypic recombination may play a major role in the evolution of some of these viruses. The phylogenetic pattern of the EV-C strains in a recent study suggests non-random recombination between strains of the same species [3]. The frequency of recombination appears to differ between species, types and intraspecies groupings, which is possibly due to cell tropism, sequence similarity and the phylogenetic position of a given virus strain [3].

Recombination also occurs in new EV (sero) types. For instance, the 5’ UTR of an EV-C109 strain in a study is likely the product of an interspecies recombination event between ancestral members of the EV-A and EV-C groups [146]. Another analysis of the P3 region of an EV-C104 strain indicates that this virus has undergone a rearrangement within the last decade, likely due to recombination with an EV-C117-like virus [130].

Analyses of EV prototype strains and clinical isolates suggest that interserotypic recombination is a frequent event and that it generally occurs among viruses of the same species. The exception is in the 5’ UTR, where only a single genetic group can be identified within EV-A and -B species and a second within EV-C and-D species [129,134,137]. This reflects a different pressure for evolutionary change exerted on the 5’ UTR and the rest of the genome, largely shaped by the functional significance of the specific genomic region and its interaction with the stable cellular environment [147,148]. Likewise, recombination may delink a given capsid sequence (serotype) from its original 5’ UTR, explaining why 5’ UTR is generally unsuitable as a molecular alternative to serotype identification [129]. In general, the EV genome could be recognized as a stable symbiosis of genes, and EV species consist of a finite set of capsid genes responsible for different serotypes and a continuum of non-structural protein genes that seem to evolve in a relatively independent manner [149].


EVs are characterized by high genomic plasticity, primarily due to high mutation and recombination rates. Improved molecular diagnostic methods and genetic sequence analyses are beginning to discover the complex characteristics of individual serotypes and genotypes. Understanding the tempo and pattern of molecular diversity and evolution is of great importance in the pathogenesis of EVs, information which will assist in disease prevention and control [18].

The vast majority of EVs infect the gastrointestinal and respiratory tracts, and can spread to other organs or systems. Improvements in surveillance, serological surveys and detailed genetic and antigenic characterization of viral populations may help to elucidate biological mechanisms [150]. Better characterization of these human pathogens may help to develop vaccines or antiviral treatments and to monitor the emergence of new strains [151].


  1. Geoghegan JL, Tan le V, Kuhnert D, Halpin RA, Lin X, et al. (2015) Phylodynamics of Enterovirus A71-Associated Hand, Foot, and Mouth Disease in Viet Nam. J Virol 89(17): 8871-8879.
  2. Li L, Ivanova O, Driss N, Tiongco-Recto M, da SR, et al. (2014) Poliovirus excretion among persons with primary immune deficiency disorders: summary of a seven-country study series. J Infect Dis 210(Suppl 1): S368-S372.
  3. Smura T, Blomqvist S, Vuorinen T, Ivanova O, Samoilovich E, et al. (2014) Recombination in the evolution of enterovirus C species sub-group that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99. PLoS One 9(4): e94579.
  4. Bailey JM, Tapprich WE (2007) Structure of the 5' nontranslated region of the coxsackievirus b3 genome: Chemical modification and comparative sequence analysis. J Virol 81(2): 650-668.
  5. Fernandez-Miragall O, Lopez dQ, Martinez-Salas E (2009) Relevance of RNA structure for the activity of picornavirus IRES elements. Virus Res 139(2): 172-182.
  6. Rossmann MG, Arnold E, Erickson JW, Frankenberger EA, Griffith JP, et al. (1985) Structure of a human common cold virus and functional relationship to other picornaviruses. Nature 317(6033): 145-153.
  7. Oberste MS, Michele SM, Maher K, Schnurr D, Cisterna D, N et al. (2014) Molecular identification and characterization of two proposed new enterovirus serotypes, EV74 and EV75. J Gen Virol 85(Pt 11): 3205-3212.
  8. Simmonds P, Welch J (2006) Frequency and dynamics of recombination within different species of human enteroviruses. J Virol 80(1): 483-493.
  9. Perera D, Shimizu H, Yoshida H, Tu PV, Ishiko H, et al. (2010) A comparison of the VP1, VP2, and VP4 regions for molecular typing of human enteroviruses. J Med Virol 82(4): 649-657.
  10. Oberste MS, Maher K, Kilpatrick DR, Pallansch MA (1999) Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to picornavirus classification. J Virol 73(3): 1941-1948.
  11. Zhou F, Kong F, McPhie K, Ratnamohan M, Donovan L, et al. (2009) Identification of 20 common human enterovirus serotypes by use of a reverse transcription-PCR-based reverse line blot hybridization assay. J Clin Microbiol 47(9): 2737-2743.
  12. Oberste MS, Maher K, Michele SM, Belliot G, Uddin M, et al. (2005) Enteroviruses 76, 89, 90 and 91 represent a novel group within the species Human Enterovirus A. J Gen Virol 86(Pt 2): 445-451.
  13. Oberste MS, Maher K, Nix WA, Michele SM, Uddin M, et al. (2007) Molecular identification of 13 new enterovirus types: EV79-88, EV97 and EV100-101, members of the species Human Enterovirus B. Virus Res 128(1-2): 34-42.
  14. Zhou F, Kong F, McPhie K, Ratnamohan M, Gilbert GL, et al. (2010) Molecular identification and analysis of nonserotypeable human enteroviruses. J Clin Microbiol 48(4): 1276-1282.
  15. Chen CH, Hsu BM, Wan MT (2008) Molecular detection and prevalence of enterovirus within environmental water in Taiwan. J Appl Microbiol 104(3): 817-823.
  16. Iturriza-Gomara M, Megson B, Gray J (2006) Molecular detection and characterization of human enteroviruses directly from clinical samples using RT-PCR and DNA sequencing. J Med Virol 78(2): 243-253.
  17. Zhou F, O’Sullivan MVN, Iredell JR, Sintchenko V, Gilbert GL, et al. (2005) Molecular analysis of enterovirus C species using the 5' untranslated region and VP1 region. J Hum Virol Retrovirol 2(2): 00038.
  18. Zhou F, Wang Q, Sintchenko V, Gilbert GL, O'Sullivan MV, et al. (2014) Use of the 5' untranslated region and VP1 region to examine the molecular diversity in enterovirus B species. J Med Microbiol 63(Pt 10): 1339-1355.
  19. Zhou F, Kong F, Wang B, McPhie K, Gilbert GL, et al. (2011) Molecular characterization of enterovirus 71 and coxsackievirus A16 using the 5' untranslated region and VP1 region. J Med Microbiol 60(Pt 3): 349-358.
  20. Tapparel C, Cordey S, Van BS, Turin L, Lee WM, et al. (2009) New molecular detection tools adapted to emerging rhinoviruses and enteroviruses. J Clin Microbiol 47(6): 1742-1749.
  21. Thoelen I, Lemey P, Van DD, I, Beuselinck K, Lindberg AM, et al. (2003) Molecular typing and epidemiology of enteroviruses identified from an outbreak of aseptic meningitis in Belgium during the summer of 2000. J Med Virol 70(3): 420-429.
  22. Tavakoli NP, Wang H, Nattanmai S, Dupuis M, Fusco H, et al. (2008) Detection and typing of enteroviruses from CSF specimens from patients diagnosed with meningitis/encephalitis. J Clin Virol 43(2): 207-211.
  23. Shieh YS, Wait D, Tai L, Sobsey MD (1995) Methods to remove inhibitors in sewage and other fecal wastes for enterovirus detection by the polymerase chain reaction. J Virol Methods 54(1): 51-66.
  24. Kilpatrick DR, Nottay B, Yang CF, Yang SJ, Da Silva E, et al. (1998) Serotype-specific identification of polioviruses by PCR using primers containing mixed-base or deoxyinosine residues at positions of codon degeneracy. J Clin Microbiol 36(2): 352-357.
  25. Osterback R, Kalliokoski T, Lahdesmaki T, Peltola V, Ruuskanen O, et al. (2015) Echovirus 30 meningitis epidemic followed by an outbreak-specific RT-qPCR. J Clin Virol 69: 7-11.
  26. Wylie TN, Wylie KM, Buller RS, Cannella M, Storch GA (2015) Development and Evaluation of an Enterovirus D68 Real-Time Reverse Transcriptase PCR Assay. J Clin Microbiol 53(8): 2641-2647.
  27. Tracy S (1985) Comparison of genomic homologies in the coxsackievirus B group by use of cDNA: RNA dot-blot hybridization. J Clin Microbiol 21(3): 371-374.
  28. Hohenadl C, Klingel K, Mertsching J, Hofschneider PH, Kandolf R (1991) Strand-specific detection of enteroviral RNA in myocardial tissue by in situ hybridization. Mol Cell Probes 5(1): 11-20.
  29. Kandolf R, Ameis D, Kirschner P, Canu A, Hofschneider PH (1987) In situ detection of enteroviral genomes in myocardial cells by nucleic acid hybridization: an approach to the diagnosis of viral heart disease. Proc Natl Acad Sci U S A 84(17): 6272-6276.
  30. Laiho JE, Oikarinen S, Oikarinen M, Larsson PG, Stone VM, et al. (2015) Application of bioinformatics in probe design enables detection of enteroviruses on different taxonomic levels by advanced in situ hybridization technology. J Clin Virol 69: 165-171.
  31. De L, Yang CF, Da Silva E, Boshell J, Caceres P, et al. (1997) Genotype-specific RNA probes for direct identification of wild polioviruses by blot hybridization. J Clin Microbiol 35(11): 2834-2840.
  32. De L, Nottay B, Yang CF, Holloway BP, Pallansch M, et al. (1995) Identification of vaccine-related polioviruses by hybridization with specific RNA probes. J Clin Microbiol 33(3): 562-571.
  33. Chen TC, Chen GW, Hsiung CA, Yang JY, Shih SR, et al. (2006) Combining multiplex reverse transcription-PCR and a diagnostic microarray to detect and differentiate enterovirus 71 and coxsackievirus A16. J Clin Microbiol 44(6): 2212-2219.
  34. Susi P, Hattara L, Waris M, Luoma-Aho T, Siitari H, et al. (2009) Typing of enteroviruses by use of microwell oligonucleotide arrays. J Clin Microbiol 47(6): 1863-1870.
  35. Oberste MS, Maher K, Williams AJ, Dybdahl-Sissoko N, Brown BA, et al. (2006) Species-specific RT-PCR amplification of human enteroviruses: a tool for rapid species identification of uncharacterized enteroviruses. J Gen Virol 87(Pt 1): 119-128.
  36. Witso E, Palacios G, Cinek O, Stene LC, Grinde B, et al. (2006) High prevalence of human enterovirus A infections in natural circulation of human enteroviruses. J Clin Microbiol 44(11): 4095-4100.
  37. Oberste MS, Maher K, Kilpatrick DR, Flemister MR, Brown BA, et al. (2006) Typing of human enteroviruses by partial sequencing of VP1. J Clin Microbiol 37(5): 1288-1293.
  38. Brown BA, Maher K, Flemister MR, Naraghi-Arani P, Uddin M, et al. (2009) Resolving ambiguities in genetic typing of human enterovirus species C clinical isolates and identification of enterovirus 96, 99 and 102. J Gen Virol 90(Pt 7): 1713-1723.
  39. Greninger AL, Naccache SN, Messacar K, Clayton A, Yu G, et al. (2015) A novel outbreak enterovirus D68 strain associated with acute flaccid myelitis cases in the USA (2012-14): a retrospective cohort study. Lancet Infect Dis 15(6): 671-682.
  40. Baicus A, Combiescu M, Persu A, Oprisan G, Aubert-Combiescu A, et al. (2006) The molecular characterization of poliovirus strains by the RT-PCR-RFLP assay and its use in the active surveillance for acute flaccid paralysis cases in Romania between 2001-2006. Roum Arch Microbiol Immunol 65(3-4): 120-130.
  41. Lien KY, Lee WC, Lei HY, Lee GB (2007) Integrated reverse transcription polymerase chain reaction systems for virus detection. Biosens Bioelectron 22(8): 1739-1748.
  42. Arita M, Ling H, Yan D, Nishimura Y, Yoshida H, et al. (2009) Development of a reverse transcription-loop-mediated isothermal amplification (RT-LAMP) system for a highly sensitive detection of enterovirus in the stool samples of acute flaccid paralysis cases. BMC Infect Dis 9: 208.
  43. Trallero G, Avellon A, Otero A, De MT, Perez C, et al. (2010) Enteroviruses in Spain over the decade 1998-2007: virological and epidemiological studies. J Clin Virol 47(2): 170-176.
  44. Simonen ML, Roivainen M, Iber J, Burns C, Hovi T (2010) Outbreak of poliomyelitis in Finland in 1984-85 - Re-analysis of viral sequences using the current standard approach. Virus Res 147(1): 91-97.
  45. Wang JR, Tsai HP, Huang SW, Kuo PH, Kiang D, et al. (2002) Laboratory diagnosis and genetic analysis of an echovirus 30-associated outbreak of aseptic meningitis in Taiwan in 2001. J Clin Microbiol 40(12): 4439-4444.
  46. Brown B, Oberste MS, Maher K, Pallansch MA (2003) Complete genomic sequencing shows that polioviruses and members of human enterovirus species C are closely related in the noncapsid coding region. J Virol 77(16): 8973-8984.
  47. Afif H, Sutter RW, Kew OM, Fontaine RE, Pallansch MA, et al. (1997) Outbreak of poliomyelitis in Gizan, Saudi Arabia: cocirculation of wild type 1 polioviruses from three separate origins. J Infect Dis 175 (Suppl 1): S71-S75.
  48. Savolainen C, Hovi T, Mulders MN (2001) Molecular epidemiology of echovirus 30 in Europe: succession of dominant sublineages within a single major genotype. Arch Virol 146(3): 521-537.
  49. Tee KK, Lam TT, Chan YF, Bible JM, Kamarulzaman A, et al. Evolutionary genetics of human enterovirus 71: origin, population dynamics, natural selection, and seasonal periodicity of the VP1 gene. J Virol 84(7): 3339-3350.
  50. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30(12): 2725-2729.
  51. Bian LL, Yao X, Mao QY, Gao F, Wang YP, et al. (2015) Complete Genome Sequence Analysis of Echovirus 24 Associated with Hand-Foot-and-Mouth Disease in China in 2012. Genome Announc 3(1): e01456-1e01514.
  52. Helal M, Kong F, Chen SC, Zhou F, Dwyer DE, et al. (2012) Linear normalised hash function for clustering gene sequences and identifying reference sequences from multiple sequence alignments. Microb Inform Exp 2(1): 2.
  53. Rico-Hesse R, Pallansch MA, Nottay BK, Kew OM (1987) Geographic distribution of wild poliovirus type 1 genotypes. Virology 160(2): 311-322.
  54. Manor Y, Blomqvist S, Sofer D, Alfandari J, Halmut T, et al. (2007) Advanced environmental surveillance and molecular analyses indicate separate importations rather than endemic circulation of wild type 1 poliovirus in Gaza district in 2002. Appl Environ Microbiol 73(18): 5954-5958.
  55. Blomqvist S, Savolainen C, Laine P, Hirttio P, Lamminsalo E, et al. (2004) Characterization of a highly evolved vaccine-derived poliovirus type 3 isolated from sewage in Estonia. J Virol 78(9): 4876-4883.
  56. Kew O, Morris-Glasgow V, Landaverde M, Burns C, Shaw J, et al. (2002) Outbreak of poliomyelitis in Hispaniola associated with circulating type 1 vaccine-derived poliovirus. Science 296(5566): 356-359.
  57. Yang CF, Naguib T, Yang SJ, Nasr E, Jorba J, et al. (2003) Circulation of endemic type 2 vaccine-derived poliovirus in Egypt from 1983 to 1993. J Virol 77(15): 8366-8377.
  58. Shulman LM, Handsher R, Yang CF, Yang SJ, Manor J, et al. (2000) Resolution of the pathways of poliovirus type 1 transmission during an outbreak. J Clin Microbiol 38(3): 945-952.
  59. Lipskaya GY, Chervonskaya EA, Belova GI, Maslova SV, Kutateladze TN, et al. (1995) Geographical genotypes (genotypes) of poliovirus case isolates from the former Soviet Union: relatedness to other known poliovirus genotypes. J Gen Virol 76 (Pt 7): 1687-1699.
  60. Diop OM, Burns CC, Sutter RW, Wassilak SG, Kew OM (2015) Update on vaccine-derived polioviruses - Worldwide, January 2014-March 2015. MMWR Morb Mortal Wkly Rep 64(23): 640-646.
  61. Brown BA, Oberste MS, Alexander JP, Kennett ML, Pallansch MA (1999) Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998. J Virol 73(12): 9969-9975.
  62. Kottaridi C, Bolanaki E, Mamuris Z, Stathopoulos C, Markoulatos P (2006) Molecular phylogeny of VP1, 2A, and 2B genes of echovirus isolates: epidemiological linkage and observations on genetic variation. Arch Virol 151(6): 1117-1132.
  63. Victoria JG, Kapoor A, Li L, Blinkova O, Slikas B, et al. (2009) Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis. J Virol 83(9): 4642-4651.
  64. Zhang Y, Wang D, Yan D, Zhu S, Liu J, et al. (2010) Molecular evidence of persistent epidemic and evolution of subgenotype B1 coxsackievirus A16-associated hand, foot, and mouth disease in China. J Clin Microbiol 48(2): 619-622.
  65. McMinn P, Stratov I, Nagarajan L, Davis S (2001) Neurological manifestations of enterovirus 71 infection in children during an outbreak of hand, foot, and mouth disease in Western Australia. Clin Infect Dis 32(2): 236-242.
  66. AbuBakar S, Chee HY, Al-Kobaisi MF, Xiaoshan J, Chua KB, et al. (1999) Identification of enterovirus 71 isolates from an outbreak of hand, foot and mouth disease (HFMD) with fatal cases of encephalomyelitis in Malaysia. Virus Res 61(1): 1-9.
  67. Bible JM, Iturriza-Gomara M, Megson B, Brown D, Pantelidis P, et al. (2008                       ) Molecular epidemiology of human enterovirus 71 in the United Kingdom from 1998 to 2006. J Clin Microbiol 46(10): 3192-3200.
  68. Chu PY, Lin KH, Hwang KP, Chou LC, Wang CF, et al. (2001) Molecular epidemiology of enterovirus 71 in Taiwan. Arch Virol 146(3): 589-600.
  69. McMinn P, Lindsay K, Perera D, Chan HM, Chan KP, et al. (2001) Phylogenetic analysis of enterovirus 71 strains isolated during linked epidemics in Malaysia, Singapore, and Western Australia. J Virol 75(16): 7732-7738.
  70. Sanders SA, Herrero LJ, McPhie K, Chow SS, Craig ME, et al. (2006) Molecular epidemiology of enterovirus 71 over two decades in an Australian urban community. Arch Virol 151(5): 1003-1013.
  71. Bible JM, Pantelidis P, Chan PK, Tong CY (2007) Genetic evolution of enterovirus 71: epidemiological and pathological implications. Rev Med Virol 17(6): 371-379.
  72. McMinn PC (2012) Recent advances in the molecular epidemiology and control of human enterovirus 71 infection. Curr Opin Virol 2(2): 199-205.
  73. Bessaud M, Razafindratsimandresy R, Nougairede A, Joffret ML, Deshpande JM, et al. (2014) Molecular comparison and evolutionary analyses of VP1 nucleotide sequences of new African human enterovirus 71 isolates reveal a wide genetic diversity. PLoS One 9(3): e90624.
  74. Saxena VK, Sane S, Nadkarni SS, Sharma DK, Deshpande JM (2015) Genetic diversity of enterovirus A71, India. Emerg Infect Dis 21(1): 123-126.
  75. Rao CD, Yergolkar P, Shankarappa KS (2012) Antigenic diversity of enteroviruses associated with nonpolio acute flaccid paralysis, India, 2007-2009. Emerg Infect Dis 18(11):1833-1840.
  76. Yip CC, Lau SK, Lo JY, Chan KH, Woo PC, et al. (2013) Genetic characterization of EV71 isolates from 2004 to 2010 reveals predominance and persistent circulation of the newly proposed genotype D and recent emergence of a distinct lineage of subgenotype C2 in Hong Kong. Virol J 10: 222.
  77. Perera D, Yusof MA, Podin Y, Ooi MH, Thao NT, et al. (2007) Molecular phylogeny of modern coxsackievirus A16. Arch Virol 152(6): 1201-1208.
  78. Li L, He Y, Yang H, Zhu J, Xu X, et al. (2005) Genetic characteristics of human enterovirus 71 and coxsackievirus A16 circulating from 1999 to 2004 in Shenzhen, People's Republic of China. J Clin Microbiol 43(8): 3835-3839.
  79. Iwai M, Masaki A, Hasegawa S, Obara M, Horimoto E, et al. (2009) Genetic changes of coxsackievirus A16 and enterovirus 71 isolated from hand, foot, and mouth disease patients in Toyama, Japan between 1981 and 2007. Jpn J Infect Dis 62(4): 254-259.
  80. Zhang C, Zhu R, Yang Y, Chi Y, Yin J, et al. (2015) Phylogenetic analysis of the major causative agents of hand, foot and mouth disease in Suzhou city, Jiangsu province, China, in 2012-2013. Emerg Microbes Infect 4(2): e12.
  81. Wikswo ME, Khetsuriani N, Fowlkes AL, Zheng X, Penaranda S, et al. (2009) Increased activity of Coxsackievirus B1 strains associated with severe disease among young infants in the United States, 2007-2008. Clin Infect Dis 49(5): e44-e51.
  82. Kopecka H, Brown B, Pallansch M (1995) Genotypic variation in coxsackievirus B5 isolates from three different outbreaks in the United States. Virus Res 38(2-3): 125-136.
  83. Rezig D, Ben YA, Ben AH, Bahri O, Triki H (2004) Molecular characterization of coxsackievirus B5 isolates. J Med Virol 72(2): 268-274.
  84. Chu PY, Tyan YC, Chen YS, Chen HL, Lu PL, et al. (2015) Transmission and Demographic Dynamics of Coxsackievirus B1. PLoS One 10(6): e0129272.
  85. Chevaliez S, Szendroi A, Caro V, Balanant J, Guillot S, et al. (2004) Molecular comparison of echovirus 11 strains circulating in Europe during an epidemic of multisystem hemorrhagic disease of infants indicates that evolution generally occurs by recombination. Virology 325(1): 56-70.
  86. Oberste MS, Nix WA, Kilpatrick DR, Flemister MR, Pallansch MA (2003) Molecular epidemiology and type-specific detection of echovirus 11 isolates from the Americas, Europe, Africa, Australia, southern Asia and the Middle East. Virus Res 91(2): 241-248.
  87. Lukashev AN, Lashkevich VA, Koroleva GA, Ilonen J, Karganova GG, et al. (2003) Molecular epidemiology of enteroviruses causing uveitis and multisystem hemorrhagic disease of infants. Virology 307(1): 45-53.
  88. Croker C, Civen R, Keough K, Ngo V, Marutani A, et al. (2015) Aseptic meningitis outbreak associated with echovirus 30 among high school football players--Los Angeles County, California, 2014. MMWR Morb Mortal Wkly Rep 63(51): 1228.
  89. Oberste MS, Maher K, Kennett ML, Campbell JJ, Carpenter MS, et al. (1999) Molecular epidemiology and genetic diversity of echovirus type 30 (E30): genotypes correlate with temporal dynamics of E30 isolation. J Clin Microbiol 37(12): 3928-3933.
  90. Lukashev AN, Ivanova OE, Eremeeva TP, Gmyl LV (2008) Analysis of echovirus 30 isolates from Russia and new independent states revealing frequent recombination and reemergence of ancient lineages. J Clin Microbiol 46(2): 665-670.
  91. Palacios G, Casas I, Cisterna D, Trallero G, Tenorio A, et al. (2002) Molecular epidemiology of echovirus 30: temporal circulation and prevalence of single lineages. J Virol 76(10): 4940-4949.
  92. Brown B, Oberste MS, Maher K, Pallansch MA (2003) Complete genomic sequencing shows that polioviruses and members of human enterovirus species C are closely related in the noncapsid coding region. J Virol 77(16): 8973-8984.
  93. Smura T, Blomqvist S, Vuorinen T, Ivanova O, Samoilovich E, et al. (2014) The evolution of Vp1 gene in enterovirus C species sub-group that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99. PLoS One 9(4): e93737.
  94. Lin KH, Chern CL, Chu PY, Chang CH, Wang HL, et al. (2001) Genetic analysis of recent Taiwanese isolates of a variant of coxsackievirus A24. J Med Virol 64(3): 269-274.
  95. Wu D, Ke CW, Mo YL, Sun LM, Li H, et al. (2033) Multiple outbreaks of acute hemorrhagic conjunctivitis due to a variant of coxsackievirus A24: Guangdong, China, 2007. J Med Virol 80(10): 1762-1768.
  96. Levy A, Roberts J, Lang J, Tempone S, Kesson A, et al. (2015) Enterovirus D68 disease and molecular epidemiology in Australia. J Clin Virol 69: 117-121.
  97. Tokarz R, Firth C, Madhi SA, Howie SR, Wu W, et al. (2012) Worldwide emergence of multiple clades of enterovirus 68. J Gen Virol 93(Pt 9): 1952-1958.
  98. Imamura T, Suzuki A, Lupisan S, Okamoto M, Aniceto R, et al. (2013) Molecular evolution of enterovirus 68 detected in the Philippines. PLoS One 8(9): e74221.
  99. Ikeda T, Mizuta K, Abiko C, Aoki Y, Itagaki T, et al. (2012) Acute respiratory infections due to enterovirus 68 in Yamagata, Japan between 2005 and 2010. Microbiol Immunol 56(2): 139-143.
  100. Meijer A, van der Sanden S, Snijders BE, Jaramillo-Gutierrez G, Bont L, et al. (2012) Emergence and epidemic occurrence of enterovirus 68 respiratory infections in The Netherlands in 2010. Virology 423(1): 49-57.
  101. Imamura T, Oshitani H (2015) Global reemergence of enterovirus D68 as an important pathogen for acute respiratory infections. Rev Med Virol 25(2): 102-114.
  102. Domingo E, Holland JJ (1997) RNA virus mutations and fitness for survival. Annu Rev Microbiol 51: 151-178.
  103. Cordey S, Petty TJ, Schibler M, Martinez Y, Gerlach D, et al. (2012) Identification of site-specific adaptations conferring increased neural cell tropism during human enterovirus 71 infection. PLoS Pathog 8(7): e1002826.
  104. Kung CM, King CC, Lee CN, Huang LM, Lee PI, et al. (2007) Differences in replication capacity between enterovirus 71 isolates obtained from patients with encephalitis and those obtained from patients with herpangina in Taiwan. J Med Virol 79(1): 60-68.
  105. Drake JW, Holland JJ (1999) Mutation rates among RNA viruses. Proc Natl Acad Sci U S A 96(24): 13910-13913.
  106. Gavrilin GV, Cherkasova EA, Lipskaya GY, Kew OM, Agol VI (2000) Evolution of circulating wild poliovirus and of vaccine-derived poliovirus in an immunodeficient patient: a unifying model. J Virol 74(16): 7381-7390.
  107. Martin J, Dunn G, Hull R, Patel V, Minor PD (200) Evolution of the Sabin strain of type 3 poliovirus in an immunodeficient patient during the entire 637-day period of virus excretion. J Virol 74(7): 3001-3010.
  108. Willian S, Tracy S, Chapman N, Leser S, Romero J, et al. (2000) Mutations in a conserved enteroviral RNA oligonucleotide sequence affect positive strand viral RNA synthesis. Arch Virol 145(10): 2061-2086.
  109. Shiroki K, Ishii T, Aoki T, Ota Y, Yang WX, et al. (1997) Host range phenotype induced by mutations in the internal ribosomal entry site of poliovirus RNA. J Virol 71(1): 1-8.
  110. De JN, Franco D, Paul A, Wimmer E, Cello J (2005) Mutation of a single conserved nucleotide between the cloverleaf and internal ribosome entry site attenuates poliovirus neurovirulence. J Virol 79(22): 14235-14243.
  111. Arita M, Shimizu H, Nagata N, Ami Y, Suzaki Y, et al. Temperature-sensitive mutants of enterovirus 71 show attenuation in cynomolgus monkeys. J Gen Virol 86(Pt 5): 1391-1401.
  112. Yin H, Berg AK, Westman J, Hellerstrom C, Frisk G (2002) Complete nucleotide sequence of a Coxsackievirus B-4 strain capable of establishing persistent infection in human pancreatic islet cells: effects on insulin release, proinsulin synthesis, and cell morphology. J Med Virol 68(4): 544-557.
  113. Chapman NM, Tu Z, Tracy S, Gauntt CJ (1994) An infectious cDNA copy of the genome of a non-cardiovirulent coxsackievirus B3 strain: its complete sequence analysis and comparison to the genomes of cardiovirulent coxsackieviruses. Arch Virol 135(1-2): 115-130.
  114. Harvala H, Kalimo H, Dahllund L, Santti J, Hughes P, et al. (2002) Mapping of tissue tropism determinants in coxsackievirus genomes. J Gen Virol 83(Pt 7): 1697-1706.
  115. Henquell C, Mirand A, Richter J, Schuffenecker I, Bottiger B, et al. (2013) Phylogenetic patterns of human coxsackievirus B5 arise from population dynamics between two genogroups and reveal evolutionary factors of molecular adaptation and transmission. J Virol 87(22): 12249-12259.
  116. McWilliam Leitch EC, Bendig J, Cabrerizo M, Cardosa J, et al. (2009) Transmission networks and population turnover of echovirus 30. J Virol 83(5): 2109-2118.
  117. Bailly JL, Chambon M, Henquell C, Icart J, Peigue-Lafeuille H (2000) Genomic variations in echovirus 30 persistent isolates recovered from a chronically infected immunodeficient child and comparison with the reference strain. J Clin Microbiol 38(2): 552-557.
  118. McWilliam Leitch EC, Cabrerizo M, Cardosa J, Harvala H, Ivanova OE, et al. (2010) Evolutionary dynamics and temporal/geographical correlates of recombination in the human enterovirus echovirus types 9, 11, and 30. J Virol 84(18): 9292-9300.
  119. Cabrerizo M, Trallero G, Simmonds P (2014) Recombination and evolutionary dynamics of human echovirus 6. J Med Virol 86(5): 857-864.
  120. Liu HM, Zheng DP, Zhang LB, Oberste MS, Pallansch MA, et al. (2000) Molecular evolution of a type 1 wild-vaccine poliovirus recombinant during widespread circulation in China. J Virol 74(23): 11153-11161.
  121. Yang CF, Chen HY, Jorba J, Sun HC, Yang SJ, et al. (2005) Intratypic recombination among lineages of type 1 vaccine-derived poliovirus emerging during chronic infection of an immunodeficient patient. J Virol 79(20): 12623-12634.
  122. Kew OM, Sutter RW, de Gourville EM, Dowdle WR, Pallansch MA (2005) Vaccine-derived polioviruses and the endgame strategy for global polio eradication. Annu Rev Microbiol 59: 587-635.
  123. Yeo DS, Seah SG, Chew JS, Lim EA, Liaw JC, et al. (2007) Molecular identification of coxsackievirus A24 variant, isolated from an outbreak of acute hemorrhagic conjunctivitis in Singapore in 2005. Arch Virol 152(11): 2005-2016.
  124. Miyamura K, Takeda N, Tanimura M, Ogino T, Yamazaki S, et al. (1990) Evolutionary study on the Coxsackievirus A 24 variant causing acute hemorrhagic conjunctivitis by oligonucleotide mapping analysis of RNA genome. Arch Virol 114(1-2): 37-51.
  125. Chu PY, Ke GM, Chang CH, Lin JC, Sun CY, et al. (2009) Molecular epidemiology of coxsackie A type 24 variant in Taiwan, 2000-2007. J Clin Virol 45(4): 285-291.
  126. Linsuwanon P, Puenpa J, Suwannakarn K, Auksornkitti V, Vichiwattana P, et al. (2012) Molecular epidemiology and evolution of human enterovirus serotype 68 in Thailand, 2006-2011. PLoS One 7(5): e35190.
  127. Imamura T, Okamoto M, Nakakita S, Suzuki A, Saito M, et al. (2014) Antigenic and receptor binding properties of enterovirus 68. J Virol 88(5): 2374-2384.
  128. Andersson P, Edman K, Lindberg AM (2002) Molecular analysis of the echovirus 18 prototype: evidence of interserotypic recombination with echovirus 9. Virus Res 85(1): 71-83.
  129. Oberste MS, Maher K, Pallansch MA (2004) Evidence for frequent recombination within species human enterovirus B based on complete genomic sequences of all thirty-seven serotypes. J Virol 78(2): 855-867.
  130. Tokarz R, Haq S, Sameroff S, Howie SR, Lipkin WI (2013) Genomic analysis of coxsackieviruses A1, A19, A22, enteroviruses 113 and 104: viruses representing two clades with distinct tropism within enterovirus C. J Gen Virol 94(Pt 9): 1995-2004.
  131. Lukashev AN, Lashkevich VA, Ivanova OE, Koroleva GA, Hinkkanen AE, et al. (2005) Recombination in circulating Human enterovirus B: independent evolution of structural and non-structural genome regions. J Gen Virol 86(Pt 12): 3281-3290.
  132. Kirkegaard K, Baltimore D (1986) The mechanism of RNA recombination in poliovirus. Cell 47(3): 433-443.
  133. Tang RS, Barton DJ, Flanegan JB, Kirkegaard K (1997) Poliovirus RNA recombination in cell-free extracts. RNA 3(6): 624-633.
  134. Oberste MS, Penaranda S, Maher K, Pallansch MA (2004) Complete genome sequences of all members of the species Human enterovirus A. J Gen Virol 85(Pt 6): 1597-1607.
  135. Smura TP, Junttila N, Blomqvist S, Norder H, Kaijalainen S, et al. (2007) Enterovirus 94, a proposed new serotype in human enterovirus species D. J Gen Virol 88(Pt 3): 849-858.
  136. Harvala H, Sharp CP, Ngole EM, Delaporte E, Peeters M, et al. (2011) Detection and genetic characterization of enteroviruses circulating among wild populations of chimpanzees in Cameroon: relationship with human and simian enteroviruses. J Virol 85(9): 4480-4486.
  137. Chen X, Zhang Q, Li J, Cao W, Zhang JX, et al. (2010) Analysis of recombination and natural selection in human enterovirus 71. Virology 398(2): 251-261.
  138. Herrero LJ, Lee CS, Hurrelbrink RJ, Chua BH, Chua KB, et al. (2003) Molecular epidemiology of enterovirus 71 in peninsular Malaysia, 1997-2000. Arch Virol 148(7): 1369-1385.
  139. Liu W, Wu S, Xiong Y, Li T, Wen Z, et al. (2014) Co-circulation and genomic recombination of coxsackievirus A16 and enterovirus 71 during a large outbreak of hand, foot, and mouth disease in Central China. PLoS One 9(4): e96051.
  140. Oberste MS, Penaranda S, Pallansch MA (2004) RNA recombination plays a major role in genomic change during circulation of coxsackie B viruses. J Virol 78(6): 2948-2955.
  141. Chen GW, Huang JH, Lo YL, Tsao KC (2007) Chang SC Mosaic genome structure of echovirus type 30 that circulated in Taiwan in 2001. Arch Virol 152(10): 1807-1817.
  142. Liu HM, Zheng DP, Zhang LB, Oberste MS, Kew OM, et al. (2003) Serial recombination during circulation of type 1 wild-vaccine recombinant polioviruses in China. J Virol 77(20): 10994-11005.
  143. Cuervo NS, Guillot S, Romanenkova N, Combiescu M, Aubert-Combiescu A, et al. (2001) Genomic features of intertypic recombinant sabin poliovirus strains excreted by primary vaccinees. J Virol 75(13): 5740-5751.
  144. Liu HM, Zheng DP, Zhang LB, Oberste MS, Kew OM, et al. (2013) Serial recombination during circulation of type 1 wild-vaccine recombinant polioviruses in China. J Virol 77(20): 10994-11005.
  145. Rakoto-Andrianarivelo M, Guillot S, Iber J, Balanant J, Blondel B, et al. (2007) Co-circulation and evolution of polioviruses and species C enteroviruses in a district of Madagascar. PLoS Pathog 3(12): e191.
  146. Yozwiak NL, Skewes-Cox P, Gordon A, Saborio S, Kuan G, et al. (2010) Human enterovirus 109: a novel interspecies recombinant enterovirus isolated from a case of acute pediatric respiratory illness in Nicaragua. J Virol 84(18): 9047-9058.
  147. Lukashev AN, Lashkevich VA, Ivanova OE, Koroleva GA, Hinkkanen AE, et al. (2003) Recombination in circulating enteroviruses. J Virol 77(19): 10423-10431.
  148. Poyry T, Kinnunen L, Hyypia T, Brown B, Horsnell C, et al. (1996) Genetic and phylogenetic clustering of enteroviruses. J Gen Virol 77 (Pt 8): 1699-1717.
  149. Siafakas N, Papaventsis D, Levidiotou-Stefanou S, Vamvakopoulos NC, Markoulatos P (2005) Classification and structure of echovirus 5'-UTR sequences. Virus Genes 31(3): 293-306.
  150. Pons-Salort M, Parker EP, Grassly NC (2015) The epidemiology of non-polio enteroviruses: recent advances and outstanding questions. Curr Opin Infect Dis 28: 479-487.
  151. Tapparel C, Siegrist F, Petty TJ, Kaiser L (2013) Picornavirus and enterovirus diversity with associated human diseases. Infect Genet Evol 14: 282-293.
© 2014-2016 MedCrave Group, All rights reserved. No part of this content may be reproduced or transmitted in any form or by any means as per the standard guidelines of fair use.
Creative Commons License Open Access by MedCrave Group is licensed under a Creative Commons Attribution 4.0 International License.
Based on a work at
Best viewed in Mozilla Firefox | Google Chrome | Above IE 7.0 version | Opera |Privacy Policy