Most RNA viruses infecting mammals and other vertebrates show profound suppression of CpG and UpA dinucleotide frequencies. To investigate this functionally, mutants of the picornavirus, echovirus 7 (E7), were constructed with altered CpG and UpA compositions in two 1.1–1.3 Kbase regions. Those with increased frequencies of CpG and UpA showed impaired replication kinetics and higher RNA/infectivity ratios compared with wild-type virus. Remarkably, mutants with CpGs and UpAs removed showed enhanced replication, larger plaques and rapidly outcompeted wild-type virus on co-infections. Luciferase-expressing E7 sub-genomic replicons with CpGs and UpAs removed from the reporter gene showed 100-fold greater luminescence. E7 and mutants were equivalently sensitive to exogenously added interferon-β, showed no evidence for differential recognition by ADAR1 or pattern recognition receptors RIG-I, MDA5 or PKR. However, kinase inhibitors roscovitine and C16 partially or entirely reversed the attenuated phenotype of high CpG and UpA mutants, potentially through inhibition of currently uncharacterized pattern recognition receptors that respond to RNA composition. Generating viruses with enhanced replication kinetics has applications in vaccine production and reporter gene construction. More fundamentally, the findings introduce a new evolutionary paradigm where dinucleotide composition of viral genomes is subjected to selection pressures independently of coding capacity and profoundly influences host–pathogen interactions.

Discrete RNA secondary and higher-order structures, typically local in extent, play a fundamental role in RNA virus replication. Using new bioinformatics analysis methods, we have identified genome-scale ordered RNA structure (GORS) in many genera and families of positive-strand animal and plant RNA viruses. There was remarkably variability between genera that possess this characteristic; for example, hepaciviruses in the family Flaviviridae show evidence for extensive internal base-pairing throughout their coding sequences that was absent in both the related pestivirus and flavivirus genera. Similar genus-associated variability was observed in the Picornaviridae, the Caliciviridae, and many plant virus families. The similarity in replication strategies between genera in each of these families rules out a role for GORS in a fundamentally conserved aspect of this aspect of the virus life cycle. However, in the Picornaviridae, Flaviviridae, and Caliciviridae, the existence of GORS correlated strongly with the ability of each genus to persist in their natural hosts. This raises the intriguing possibility of a role for GORS in the modulation of innate intracellular defense mechanisms (and secondarily, the acquired immune system) triggered by double-stranded RNA, analogous in function to the expression of structured RNA transcripts by large DNA viruses. Irrespective of function, the observed evolutionary conservation of GORS in many viruses imposes a considerable constraint on genome plasticity and the consequent narrowing of sequence space in which neutral drift can occur. These findings potentially reconcile the rapid evolution of RNA viruses over short periods with the documented examples of extreme conservatism evident from their intimate coevolution with their hosts.

We have used high-throughput Illumina sequencing to identify novel recombinants between Deformed wing virus (DWV) and Varroa destructor virus-1 (VDV-1), which accumulate to higher levels than DWV in both honeybees and Varroa destructor mites. The recombinants, VDV-1VVD and VDV-1DVD, exhibit crossovers between the 5′-UTR and the regions encoding the structural (capsid) and non-structural viral proteins. This implies that the genomes are modular and that each region may evolve independently, as demonstrated in human enteroviruses. Individual honeybee pupae were infected with a mixture of observed recombinants and DWV. A strong correlation was observed between VDV-1DVD levels in honeybee pupae and associated mites, suggesting that this recombinant, with a DWV-derived 5′-UTR and non-structural protein region flanking a VDV-1-derived capsid-encoding region, is better adapted to transmission between V. destructor and honeybees than the parental DWV or a recombinant bearing the VDV-1-derived 5′-UTR (VDV-1VVD).

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.

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.