Recombination is a well-known phenomenon for enteroviruses. However, the actual extent of recombination in circulating nonpoliovirus enteroviruses is not known. We have analyzed the phylogenetic relationships in four genome regions, VP1, 2A, 3D, and the 5′ nontranslated region (NTR), of 40 enterovirus B strains (coxsackie B viruses and echoviruses) representing 11 serotypes and isolated in 1981 to 2002 in the former Soviet Union states. In the VP1 region, strains of the same serotype expectedly grouped with their prototype strain. However, as early as the 2A region, phylogenetic grouping differed significantly from that in the VP1 region and indicated recombination within the 2A region. Moreover, in the 5′ NTR and 3D region, only 1 strain of 40 grouped with its prototype strain. Instead, we observed a major group in both the 5′ NTR and the 3D region that united most (in the 5′ NTR) or all (in the 3D region) of the strains studied, regardless of the serotype. Subdivision within that major group in the 3D region correlated with the time of virus isolation but not with the serotype. Therefore, we conclude that a majority, if not all, circulating enterovirus B strains are recombinants relative to the prototype strains, isolated mostly in the 1950s. Moreover, the ubiquitous recombination has allowed different regions of the enterovirus genome to evolve independently. Thus, a novel model of enterovirus genetics is proposed: the enterovirus genome is a stable symbiosis of genes, and enterovirus species consist of a finite set of capsid genes responsible for different serotypes and a continuum of nonstructural protein genes that seem to evolve in a relatively independent manner.

The complete nucleotide sequences of eight Human enterovirus B (HEV-B) strains were determined, representing five serotypes, E6, E7, E11, CVB3 and CVB5, which were isolated in the former Soviet Union between 1998 and 2002. All strains were mosaic recombinants and only the VP2–VP3–VP1 genome region was similar to that of the corresponding prototype HEV-B strains. In seven of the eight strains studied, the 2C–3D genome region was most similar to the prototype E30, EV74 and EV75 strains, whilst the remaining strain was most similar to the prototype E1 and E9 strains in the non-structural protein genome region. Most viruses also bore marks of additional recombination events in this part of the genome. In the 5′ non-translated region, all strains were more similar to the prototype E9 than to other enteroviruses. In most cases, recombination mapped to the VP4 and 2ABC genome regions. This, together with the star-like topology of the phylogenetic trees for these genome regions, identified these genome parts as recombination hot spots. These findings further support the concept of independent evolution of enterovirus genome fragments and indicate a requirement for more advanced typing approaches. A range of available phylogenetic methods was also compared for efficient detection of recombination in enteroviruses.

Enterovirus 71 (EV71) is responsible for frequent large-scale outbreaks of hand, foot, and mouth disease worldwide and represent a major etiological agent of severe, sometimes fatal neurological disease. EV71 variants have been classified into three genogroups (GgA, GgB, and GgC), and the latter two are further subdivided into subgenogroups B1 to B5 and C1 to C5. To investigate the dual roles of recombination and evolution in the epidemiology and transmission of EV71 worldwide, we performed a large-scale genetic analysis of isolates ( n = 308) collected from 19 countries worldwide over a 40-year period. A series of recombination events occurred over this period, which have been identified through incongruities in sequence grouping between the VP1 and 3Dpol regions. Eleven 3Dpol clades were identified, each specific to EV71 and associated with specific subgenogroups but interspersed phylogenetically with clades of coxsackievirus A16 and other EV species A serotypes. The likelihood of recombination increased with VP1 sequence divergence; mean half-lives for EV71 recombinant forms (RFs) of 6 and 9 years for GgB and GgC overlapped with those observed for the EV-B serotypes, echovirus 9 (E9), E30, and E11, respectively (1.3 to 9.8 years). Furthermore, within genogroups, sporadic recombination events occurred, such as the linkage of two B4 variants to RF-W instead of RF-A and of two C4 variants to RF-H. Intriguingly, recombination events occurred as a founding event of most subgenogroups immediately preceding their lineage expansion and global emergence. The possibility that recombination contributed to their subsequent spread through improved fitness requires further biological and immunological characterization.

Genome stability is a prerequisite for the production and use of adenoviruses for therapy of genetic diseases and cancer. To test the premise that the adenoviral genome is stable, the phylogenetic relationships of 16 adenovirus C (AdC) field isolates were studied in four genome regions: hexon, fiber, polymerase and E1A. The phylogenetic relationships in the fiber gene concurred with those in the hexon region. In contrast, the non-structural regions had marks of frequent recombination, to the point that an isolate of one serotype could contain non-structural proteins that were identical to the genes from a different serotype. Our results suggest that recombination among circulating adenoviruses is very frequent and plays an important role in shaping the phylogenetic relationships of adenovirus genomes. Analysis of the available complete genome sequences of AdB, AdC and AdD species showed that recombination shuffles genome fragments within a species, but not between species. One of the AdC field isolates possessed the fiber gene of AdC type 6, but a hexon gene that was distinct from all AdC serotypes. This strain could not be typed unambiguously in a neutralization test and might represent a novel serotype of AdC. Comparison of the right end (nt 18838–33452) of this isolate with that of the ATCC Ad6 strain showed clear evidence of multiple recombination events.

The relationship between virus evolution and recombination in species B human enteroviruses was investigated through large-scale genetic analysis of echovirus type 9 (E9) and E11 isolates ( n = 85 and 116) from 16 European, African, and Asian countries between 1995 and 2008. Cluster 1 E9 isolates and genotype D5 and A E11 isolates showed evidence of frequent recombination between the VP1 and 3Dpol regions, the latter falling into 23 (E9) and 43 (E11) clades interspersed phylogenetically with 46 3Dpol clades of E30 and with those of other species B serotypes. Remarkably, only 2 of the 112 3Dpol clades were shared by more than one serotype (E11 and E30), demonstrating an extremely large and genetically heterogeneous recombination pool of species B nonstructural-region variants. The likelihood of recombination increased with geographical separation and time, and both were correlated with VP1 divergence, whose substitution rates allowed recombination half-lives of 1.3, 9.8, and 3.1 years, respectively, for E9, E11, and E30 to be calculated. These marked differences in recombination dynamics matched epidemiological patterns of periodic epidemic cycles of 2 to 3 (E9) and 5 to 6 (E30) years and the longer-term endemic pattern of E11 infections. Phylotemporal analysis using a Bayesian Markov chain Monte Carlo method, which placed recombination events within the evolutionary reconstruction of VP1, showed a close relationship with VP1 lineage expansion, with defined recombination events that correlated with their epidemiological periodicity. Whether recombination events contribute directly to changes in transmissibility that drive epidemic behavior or occur stochastically during periodic population bottlenecks is an unresolved issue vital to future understanding of enterovirus molecular epidemiology and pathogenesis.

Globally, echovirus 30 (E30) is one of the most frequently identified enteroviruses and a major cause of meningitis. Despite its wide distribution, little is known about its transmission networks or the dynamics of its recombination and geographical spread. To address this, we have conducted an extensive molecular epidemiology and evolutionary study of E30 isolates collected over 8 years from a geographically wide sample base (11 European countries, Asia, and Australia). 3Dpol sequences fell into several distinct phylogenetic groups, interspersed with other species B serotypes, enabling E30 isolates to be classified into 38 recombinant forms (RFs). Substitutions in VP1 and 3Dpol regions occurred predominantly at synonymous sites (ratio of nonsynonymous to synonymous substitutions, 0.05) with VP1 showing a rapid substitution rate of 8.3 × 10⁻³substitutions per site per year. Recombination frequency was tightly correlated with VP1 divergence; viruses differing by evolutionary distances of >0.1 (or 6 years divergent evolution) almost invariably (>97%) had different 3Dpol groups. Frequencies of shared 3Dpol groups additionally correlated with geographical distances, with Europe and South Asia showing turnover of entirely distinct virus populations. Population turnover of E30 was characterized by repeated cycles of emergence, dominance, and disappearance of individual RFs over periods of 3 to 5 years, although the existence and nature of evolutionary selection underlying these population replacements remain unclear. The occurrence of frequent “sporadic” recombinants embedded within VP1 groupings of other RFs and the much greater number of 3Dpol groups than separately identifiable VP1 lineages suggest frequent recombination with an external diverse reservoir of non-E30 viruses.