Environment determines fidelity for an RNA virus replicase

Evolutionary deletions within the SARS-CoV-2 genome as signature trends for virus fitness and adaptation

Coronaviruses are large RNA viruses that can infect and spread among humans and animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 2019, has evolved since its first detection in December 2019. Deletions are a common occurrence in SARS-CoV-2 evolution, particularly in specific genomic sites, and may be associated with the emergence of highly competent lineages. While deletions typically have a negative impact on viral fitness, some persist and become fixed in viral populations, indicating that they may confer advantageous benefits for the virus's adaptive evolution. This work presents a literature review and data analysis on structural losses in the SARS-CoV-2 genome and the potential relevance of specific signatures for enhanced viral fitness and spread.

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

L-shaped distribution of the relative substitution rate (c/μ) observed for SARS-COV-2's genome, inconsistent with the selectionist theory, the neutral theory and the nearly neutral theory but a near-neutral balanced selection theory: Implication on "neutralist-selectionist" debate

The genomic substitution rate (GSR) of SARS-CoV-2 exhibits a molecular clock feature and does not change under fluctuating environmental factors such as the infected human population (10°-107), vaccination etc. The molecular clock feature is believed to be inconsistent with the selectionist theory (ST). The GSR shows lack of dependence on the effective population size, suggesting Ohta's nearly neutral theory (ONNT) is not applicable to this virus. Big variation of the substitution rate within its genome is also inconsistent with Kimura's neutral theory (KNT). Thus, all three existing evolution theories fail to explain the evolutionary nature of this virus. In this paper, we proposed a Segment Substitution Rate Model (SSRM) under non-neutral selections and pointed out that a balanced mechanism between negative and positive selection of some segments that could also lead to the molecular clock feature. We named this hybrid mechanism as near-neutral balanced selection theory (NNBST) and examined if it was followed by SARS-CoV-2 using the three independent sets of SARS-CoV-2 genomes selected by the Nextstrain team. Intriguingly, the relative substitution rate of this virus exhibited an L-shaped probability distribution consisting with NNBST rather than Poisson distribution predicted by KNT or an asymmetric distribution predicted by ONNT in which nearly neutral sites are believed to be slightly deleterious only, or the distribution that is lack of nearly neutral sites predicted by ST. The time-dependence of the substitution rates for some segments and their correlation with the vaccination were observed, supporting NNBST. Our relative substitution rate method provides a tool to resolve the long standing "neutralist-selectionist" controversy. Implications of NNBST in resolving Lewontin's Paradox is also discussed.

Influence of viral genome properties on polymerase fidelity

The first step of viral evolution takes place during genome replication via the error-prone viral polymerase. Among the mutants that arise through this process, only a few well-adapted variants will be selected by natural selection, renewing the viral genome population. Viral polymerase-mediated errors are thought to occur stochastically. However, accumulating evidence suggests that viral polymerase-mediated mutations are heterogeneously distributed throughout the viral genome. Here, we review work that supports this concept and provides mechanistic insights into how specific features of the viral genome could modulate viral polymerase-mediated errors. A predisposition to accumulate viral polymerase-mediated errors at specific loci in the viral genome may guide evolution to specific pathways, thus opening new directions of research to better understand viral evolutionary dynamics.

The Ups and Downs of an Out-of-the-Box Scientist with a Curious Mind

My early life was challenging, and not conducive to the study of science, but my first introduction to viruses was an epiphany for me. I spent the whole of my career dedicated to understanding viruses, driven largely by curiosity. This led me down many different avenues of study, and to work with many wonderful colleagues, most of whom remain friends. Some highlights of my career include the discovery of a mutualistic three-way symbiosis involving a virus, a fungus, and a plant; genetic mapping of a pathogenicity gene in tomato; uncovering a virus in 1,000-year-old corncobs; exploring virus biodiversity in wild plants; and establishing a system to use a fungal virus to understand the epidemiology of its host. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Effects of Host and Pathogenicity on Mutation Rates in Avian Influenza A Viruses

Mutation is the primary determinant of genetic diversity in influenza viruses. The rate of mutation, measured in an absolute time-scale, is likely to be dependent on the rate of errors in copying RNA sequences per replication and the number of replications per unit time. Conditions for viral replication are probably different among host taxa, potentially generating the host-specificity of the viral mutation rate, and possibly between highly and low pathogenic viruses. This study investigated whether mutation rates per year in avian influenza A viruses depend on host taxa and pathogenicity. We inferred mutation rates from the rates of synonymous substitutions, which are assumed to be neutral and thus equal to mutation rates, at four segments that code internal viral proteins (PB2, PB1, PA, NP). On the phylogeny of all avian viral sequences for each segment, multiple distinct subtrees (clades) were identified that represent viral subpopulations, which are likely to have evolved within particular host taxa. Using simple regression analysis, we found that mutation rates were significantly higher in viruses infecting chickens than domestic ducks, and in those infecting wild shorebirds than wild ducks. Host-dependency of the substitution rate was also confirmed by Bayesian phylogenetic analysis. However, we did not find evidence that the mutation rate is higher in highly pathogenic than in low pathogenic viruses. We discuss these results considering viral replication rate as the major determinant of mutation rate per unit time.

Phylodynamic Study of the Conserved RNA Structure Encompassing the Hemagglutinin Cleavage Site Encoding Region of H5 and H7 Low Pathogenic Avian Influenza Viruses

Highly Pathogenic Avian Influenza Viruses (HPAIV) evolve from Low Pathogenic Avian Influenza Viruses (LPAIV) of the H5 and H7 subtypes. This evolution is characterized by the acquisition of a multi-basic cleavage site (MBCS) motif in the hemagglutinin (HA) that leads to an extended viral tropism and severe disease in poultry. One key unanswered question is whether the risk of transition to HPAIV is similar for all LPAIV H5 or H7 strains, or whether specific determinants in the HA sequence of some H5 or H7 LPAIV strains correlate with a higher risk of transition to HPAIV. Here we determined if specific features of the conserved RNA stem loop located at the hemagglutinin cleavage site-encoding region could be detected along the LPAIV to HPAIV evolutionary pathway. Analysis of the thermodynamic stability of the predicted RNA structures showed no specific patterns common to HA sequences leading to HPAIV and distinct from those remaining LPAIV. However, RNA structure clustering analysis revealed that most of the American lineage ancestors leading to H7 emergences via recombination shared the same vRNA structure topology at the HA1/HA2 boundary region. Our study thus identified predicted secondary RNA structures present in the HA of H7 viruses, which could promote genetic recombination and acquisition of a MBCS.

Rapid protein sequence evolution via compensatory frameshift is widespread in RNA virus genomes

Background

RNA viruses possess remarkable evolutionary versatility driven by the high mutability of their genomes. Frameshifting nucleotide insertions or deletions (indels), which cause the premature termination of proteins, are frequently observed in the coding sequences of various viral genomes. When a secondary indel occurs near the primary indel site, the open reading frame can be restored to produce functional proteins, a phenomenon known as the compensatory frameshift.

Results

In this study, we systematically analyzed publicly available viral genome sequences and identified compensatory frameshift events in hundreds of viral protein-coding sequences. Compensatory frameshift events resulted in large-scale amino acid differences between the compensatory frameshift form and the wild type even though their nucleotide sequences were almost identical. Phylogenetic analyses revealed that the evolutionary distance between proteins with and without a compensatory frameshift were significantly overestimated because amino acid mismatches caused by compensatory frameshifts were counted as substitutions. Further, this could cause compensatory frameshift forms to branch in different locations in the protein and nucleotide trees, which may obscure the correct interpretation of phylogenetic relationships between variant viruses.

Conclusions

Our results imply that the compensatory frameshift is one of the mechanisms driving the rapid protein evolution of RNA viruses and potentially assisting their host-range expansion and adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04182-9.

Evolutionary dynamics of the SARS-CoV-2 ORF8 accessory gene

The new SARS-CoV-2 poses a significant threat to human health but many aspects of its basic biology remain unknown. Its genome encodes accessory genes that differ significantly within coronaviruses and contribute to the virus pathogenicity. Among accessory genes, open reading frame 8 (ORF8) stands out by being highly variable and showing structural changes suspected to be related with the virus ability to spread. However, the function of ORF8 remains to be elucidated, making it less studied than other SARS-CoV-2 genes. Here I show that ORF8 is poorly conserved among related coronaviruses. The ORF8 phylogeny built using 11,113 SARS-CoV-2 sequences revealed traces of a typical expanding population with a small number of highly frequent lineages. Interestingly, I detected several nonsense mutations and three main deletions in the ORF8 gene that either remove or significantly change the ORF8 protein. These findings suggest that SARS-CoV-2 can persist without a functional ORF8 protein. Deletion breakpoints were found located in predicted hairpins suggesting a possible involvement of these elements in the rearrangement process. Although the function of ORF8 remains to be elucidated, its structural plasticity and high diversity suggest an important role in SARS-CoV-2 pathogenicity.

Cassava brown streak virus Ham1 protein hydrolyses mutagenic nucleotides and is a necrosis determinant

Cassava brown streak disease (CBSD) is a leading cause of cassava losses in East and Central Africa, and is currently having a severe impact on food security. The disease is caused by two viruses within the Potyviridae family: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), which both encode atypical Ham1 proteins with highly conserved inosine triphosphate (ITP) pyrophosphohydrolase (ITPase) domains. ITPase proteins are widely encoded by plant, animal, and archaea. They selectively hydrolyse mutagenic nucleotide triphosphates to prevent their incorporation into nucleic acid and thereby function to reduce mutation rates. It has previously been hypothesized that U/CBSVs encode Ham1 proteins with ITPase activity to reduce viral mutation rates during infection. In this study, we investigate the potential roles of U/CBSV Ham1 proteins. We show that both CBSV and UCBSV Ham1 proteins have ITPase activities through in vitro enzyme assays. Deep-sequencing experiments found no evidence of the U/CBSV Ham1 proteins providing mutagenic protection during infections of Nicotiana hosts. Manipulations of the CBSV_Tanza infectious clone were performed, including a Ham1 deletion, ITPase point mutations, and UCBSV Ham1 chimera. Unlike severely necrotic wild-type CBSV_Tanza infections, infections of Nicotiana benthamiana with the manipulated CBSV infectious clones do not develop necrosis, indicating that that the CBSV Ham1 is a necrosis determinant. We propose that the presence of U/CBSV Ham1 proteins with highly conserved ITPase motifs indicates that they serve highly selectable functions during infections of cassava and may represent a euphorbia host adaptation that could be targeted in antiviral strategies.

Naturally occurring mutations in PB1 affect influenza A virus replication fidelity, virulence, and adaptability

Background

Mutations in the PB1 subunit of RNA-dependent RNA polymerase (RdRp) of influenza A virus can affect replication fidelity. Before the influenza A/H1N1 pandemic in 2009, most human influenza A/H1N1 viruses contained the avian-associated residue, serine, at position 216 in PB1. However, near the onset of the 2009 pandemic, human viruses began to acquire the mammalian-associated residue, glycine, at PB1–216, and PB1–216G became predominant in human viruses thereafter.

Methods

Using entropy-based analysis algorithm, we have previously identified several host-specific amino-acid signatures that separated avian and swine viruses from human influenza viruses. The presence of these host-specific signatures in human influenza A/H1N1 viruses suggested that these mutations were the result of adaptive genetic evolution that enabled these influenza viruses to circumvent host barriers, which resulted in cross-species transmission. We investigated the biological impact of this natural avian-to-mammalian signature substitution at PB1–216 in human influenza A/H1N1 viruses.

Results

We found that PB1–216G viruses had greater mutation potential, and were more sensitive to ribavirin than PB1–216S viruses. In oseltamivir-treated HEK293 cells, PB1–216G viruses generated mutations in viral neuraminidase at a higher rate than PB1–216S viruses. By contrast, PB1–216S viruses were more virulent in mice than PB1–216G viruses. These results suggest that the PB1-S216G substitution enhances viral epidemiological fitness by increasing the frequency of adaptive mutations in human influenza A/H1N1 viruses.

Conclusions

Our results thus suggest that the increased adaptability and epidemiological fitness of naturally arising human PB1–216G viruses, which have a canonical low-fidelity replicase, were the biological mechanisms underlying the replacement of PB1–216S viruses with a high-fidelity replicase following the emergence of pdmH1N1. We think that continued surveillance of such naturally occurring PB1–216 variants among others is warranted to assess the potential impact of changes in RdRp fidelity on the adaptability and epidemiological fitness of human A/H1N1 influenza viruses.

Electronic supplementary material

The online version of this article (10.1186/s12929-019-0547-4) contains supplementary material, which is available to authorized users.

Impact of Cultivated Hosts on the Recombination of Cucumber Mosaic Virus

Cucumber mosaic virus (CMV) is one of the most successful viruses known, infecting over 1,200 species of plants. Like other single-stranded RNA viruses, CMV is known to have a high potential for population diversity due to error-prone replication and short generation times. Recombination is also a mechanism that allows viruses to adapt to new hosts. Host genes have been identified that impact the recombination of RNA viruses by using single-cell yeast systems. To determine the impact that the natural plant host has on virus recombination, we used a high-recombination-frequency strain of CMV, LS-CMV, which belongs to subgroup II, in three different cultivated hosts: Capsicum annuum cv. Marengo (pepper), Nicotiana tabacum cv. Xanthi nc (tobacco), and Cucurbita pepo cv. Black Beauty (zucchini). The recombination frequency was calculated by using an RNA 3 reporter carrying restriction enzyme sites created by introducing silent mutations. Our results show that the recombination frequency of LS-CMV is correlated with the infected host. The recombination events in pepper were 1.8-fold higher than those in tobacco and 5-fold higher than those in zucchini. Furthermore, we observed the generation of defective RNAs in inoculated pepper plants, but not in tobacco or zucchini. These results indicate that the host is involved in both intra- and intermolecular recombination events and that hosts like pepper could foster more rapid evolution of the virus. In addition, we report for the first time the production of defective RNAs in a CMV subgroup II isolate.IMPORTANCE Recombination is an important mechanism used by viruses for their diversification and to adapt to diverse hosts. Understanding the host role in the mechanisms of evolution is important for virus disease management and controlling the emergence of new strains. This study shows the impact that cultivated hosts are playing in the evolution of CMV. Furthermore, our results and previous studies show how some specific hosts could be an ideal environment for the emergence of new viral strains.

Chikungunya Virus Fidelity Variants Exhibit Differential Attenuation and Population Diversity in Cell Culture and Adult Mice

Chikungunya virus (CHIKV) is a reemerging global health threat that produces debilitating arthritis in people. Like other RNA viruses with high mutation rates, CHIKV produces populations of genetically diverse genomes within a host. While several known CHIKV mutations influence disease severity in vertebrates and transmission by mosquitoes, the role of intrahost diversity in chikungunya arthritic disease has not been studied. In this study, high- and low-fidelity CHIKV variants, previously characterized by altered in vitro population mutation frequencies, were used to evaluate how intrahost diversity influences clinical disease, CHIKV replication, and antibody neutralization in immunocompetent adult mice inoculated in the rear footpads. Both high- and low-fidelity mutations were hypothesized to attenuate CHIKV arthritic disease, replication, and neutralizing antibody levels compared to wild-type (WT) CHIKV. Unexpectedly, high-fidelity mutants elicited more severe arthritic disease than the WT despite comparable CHIKV replication, whereas a low-fidelity mutant produced attenuated disease and replication. Serum antibody developed against both high- and low-fidelity CHIKV exhibited reduced neutralization of WT CHIKV. Using next-generation sequencing (NGS), the high-fidelity mutations were demonstrated to be genetically stable but produced more genetically diverse populations than WT CHIKV in mice. This enhanced diversification was subsequently reproduced after serial in vitro passage. The NGS results contrast with previously reported population diversities for fidelity variants, which focused mainly on part of the E1 gene, and highlight the need for direct measurements of mutation rates to clarify CHIKV fidelity phenotypes.IMPORTANCE CHIKV is a reemerging global health threat that elicits debilitating arthritis in humans. There are currently no commercially available CHIKV vaccines. Like other RNA viruses, CHIKV has a high mutation rate and is capable of rapid intrahost diversification during an infection. In other RNA viruses, virus population diversity associates with disease progression; however, potential impacts of intrahost viral diversity on CHIKV arthritic disease have not been studied. Using previously characterized CHIKV fidelity variants, we addressed whether CHIKV population diversity influences the severity of arthritis and host antibody response in an arthritic mouse model. Our findings show that CHIKV populations with greater genetic diversity can cause more severe disease and stimulate antibody responses with reduced neutralization of low-diversity virus populations in vitro The discordant high-fidelity phenotypes in this study highlight the complexity of inferring replication fidelity indirectly from population diversity.

Analysis of quasispecies variation in single and mixed viral infection

Population diversity was examined in individual and natural mixed infections of Cowpea chlorotic mottle virus (CCMV) and Cucumber mosaic virus (CMV) isolates in two systemic hosts, cowpea and Nicotiana benthamiana. Isolates of CCMV and CMV obtained from a cowpea field in Arkansas were separated biologically in cowpea and tobacco plants, respectively. After separation, individual and mixed cultures of both viruses were serially passaged ten times by mechanical inoculation in cowpea and N. benthamiana. High-fidelity reverse transcriptase-polymerase chain reaction (HiFi RT-PCR) of RNA 3, followed by cDNA cloning and sequence analysis was used to assess the quasispecies cloud size of CCMV and CMV populations in passages zero and ten of each host species. The levels of population variation were generally higher in individual infections of CCMV-Car1 and-Car2 isolates, and the CMV-Car2 isolate compared with mixed infections, in both host species, although the significance of the differences varied depending on how mutations were counted. There were no significant differences in the levels of population variation in individual and mixed infections of the CMV-Car1 isolate. Partially fixed mutations were observed in both individual and mixed infections of the CCMV-Car2 isolate in N. benthamiana and CMV-Car1 and-Car2 isolates in both cowpea and N. benthamiana.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

Deep sequencing for discovery and evolutionary analysis of plant viruses

The advent of next generation sequencing (NGS), or deep sequencing, has allowed great advances to be made in discovery, diagnostics, and evolutionary studies in plant viruses. Various methods have been used for enrichment for virus-specific nucleic acids, each of which have some drawbacks. Many novel viruses have been discovered in plants by NGS technologies, and there is a good deal of promise for more comprehensive studies in virus evolution. However, each aspect of using NGS has its caveats that need to be considered, and there is still a need for better tools of analysis, as well as method for validation of sequence variation.

Local adaptation of plant viruses: lessons from experimental evolution

For multihost pathogens, adaptation to multiple hosts has important implications for both applied and basic research. At the applied level, it is one of the main factors determining the probability and severity of emerging disease outbreaks. At the basic level, it is thought to be a key mechanism for the maintenance of genetic diversity both in host and pathogen species. In recent years, a number of evolution experiments have assessed the fate of plant virus populations replicating within and adapting to one single or to multiple hosts species. A first group of these experiments tackled the existence of trade-offs in fitness and virulence for viruses evolving either within a single hosts species or alternating between two different host species. A second set of experiments explored the role of genetic variability in susceptibility and resistance to infection among individuals from the same host species in the extent of virus local adaptation and of virulence. In general, when a single host species or genotype is available, these experiments show that local adaptation takes place, often but not always associated with a fitness trade-off. However, alternating between different host species or infecting resistant host genotypes may select for generalist viruses that experience no fitness cost. Therefore, the expected cost of generalism, arising from antagonistic pleiotropy and other genetic mechanisms generating fitness trade-offs between hosts, could not be generalized and strongly depend on the characteristics of each particular pathosystem. At the genomic level, these studies show pervasive convergent molecular evolution, suggesting that the number of accessible molecular pathways leading to adaptation to novel hosts is limited.

Genetic Drift, Purifying Selection and Vector Genotype Shape Dengue Virus Intra-host Genetic Diversity in Mosquitoes

Due to their error-prone replication, RNA viruses typically exist as a diverse population of closely related genomes, which is considered critical for their fitness and adaptive potential. Intra-host demographic fluctuations that stochastically reduce the effective size of viral populations are a challenge to maintaining genetic diversity during systemic host infection. Arthropod-borne viruses (arboviruses) traverse several anatomical barriers during infection of their arthropod vectors that are believed to impose population bottlenecks. These anatomical barriers have been associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. Whether these patterns result from stochastic sampling (genetic drift) rather than natural selection, and/or from the influence of vector genetic heterogeneity has not been elucidated. Here, we used deep sequencing of full-length viral genomes to monitor the intra-host evolution of a wild-type dengue virus isolate during infection of several mosquito genetic backgrounds. We estimated a bottleneck size ranging from 5 to 42 founding viral genomes at initial midgut infection, irrespective of mosquito genotype, resulting in stochastic reshuffling of the variant repertoire. The observed level of genetic diversity increased following initial midgut infection but significantly differed between mosquito genetic backgrounds despite a similar initial bottleneck size. Natural selection was predominantly negative (purifying) during viral population expansion. Taken together, our results indicate that dengue virus intra-host genetic diversity in the mosquito vector is shaped by genetic drift and purifying selection, and point to a novel role for vector genetic factors in the genetic breadth of virus populations during infection. Identifying the evolutionary forces acting on arboviral populations within their arthropod vector provides novel insights into arbovirus evolution.

During infection of their arthropod vectors, arthropod-borne viruses (arboviruses) such as dengue viruses traverse several anatomical barriers that are believed to cause dramatic reductions in population size. Such population bottlenecks challenge the maintenance of viral genetic diversity, which is considered critical for fitness and adaptability of arboviruses. Anatomical barriers in the vector were previously associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. However, the relative role of random processes and natural selection, and the influence of vector genetic heterogeneity have not been elucidated. In this study, we used high-throughput sequencing to monitor dengue virus genetic diversity during infection of several genetic backgrounds of their mosquito vector. Our results show that initial infection of the vector is randomly founded by only a few tens of individual virus genomes. The overall level of viral genetic diversity generated during infection was predominantly under purifying selection but differed significantly between mosquito genetic backgrounds. Thus, in addition to random evolutionary forces and the purging of deleterious mutations that shape dengue virus genetic diversity during vector infection, our results also point to a novel role for vector genetic factors in the genetic breadth of virus populations.

Plant Virus Diversity and Evolution

Historically, the majority of plant virology focused on agricultural systems. Recent efforts have expanded our knowledge of the true diversity of plant viruses by studying those viruses that infect wild, undomesticated plants. Those efforts have provided answers to basic ecological questions regarding viruses in the wild, and insights into evolutionary questions, regarding the origins of viruses. While much work has been done, we have merely scratched the surface of the diversity that is estimated to exist. In this chapter we discuss the state of our knowledge of virus diversity, both in agricultural systems as well as in native wild systems, the border between these two systems and how viruses adapt and move across this border into an artificial, domesticated environment. We look at how this diversity has affected our outlook on viruses as a whole, shifting our past view of viruses as purely antagonistic entities of destruction to one where viruses are in a mutually beneficial relationship with their hosts. Additionally, we discuss the current work that plant virology has put forth regarding the evolutionary mechanisms, the life histories, and the deep evolution of viruses.

Mutation and recombination frequencies reveal a biological contrast within strains of Cucumber mosaic virus

Recent in planta studies have shown that strains Fny and LS of Cucumber mosaic virus (CMV) display differential genetic diversities, Fny and LS having higher and lower mutation frequencies, respectively (J. S. Pita and M. J. Roossinck, J Virol 87:790–797, 2012 http://dx.doi.org/10.1128/JVI.01891-12). In this article, we show that these virus strains have differential recombination frequencies as well. However, the high-diversity Fny strain is a low-recombination virus, whereas the very-low-diversity LS strain is instead a high-recombination virus. Unlike the mutation frequency that was determined by both RNAs 1 and 2, the control elements of recombination frequency reside predominantly within RNA 2, specifically within the 2a gene.

IMPORTANCE

Recombination is an important mechanism in virus evolution that can lead to increased or decreased variation and is a major player in virus speciation events that can lead to emerging viruses. Although viral genomes show very frequent evidence of recombination, details of the mechanism involved in these events are still poorly understood. We show here that the reciprocal effects of high mutation frequency and low recombination frequency (and vice versa) involve the RNA-dependent RNA polymerase of the virus, and we speculate that these evolutionary events are related to differences in processivity for two strains of the same virus.

Sequence-Specific Fidelity Alterations Associated with West Nile Virus Attenuation in Mosquitoes

High rates of error-prone replication result in the rapid accumulation of genetic diversity of RNA viruses. Recent studies suggest that mutation rates are selected for optimal viral fitness and that modest variations in replicase fidelity may be associated with viral attenuation. Arthropod-borne viruses (arboviruses) are unique in their requirement for host cycling and may necessitate substantial genetic and phenotypic plasticity. In order to more thoroughly investigate the correlates, mechanisms and consequences of arbovirus fidelity, we selected fidelity variants of West Nile virus (WNV; Flaviviridae, Flavivirus) utilizing selection in the presence of a mutagen. We identified two mutations in the WNV RNA-dependent RNA polymerase associated with increased fidelity, V793I and G806R, and a single mutation in the WNV methyltransferase, T248I, associated with decreased fidelity. Both deep-sequencing and in vitro biochemical assays confirmed strain-specific differences in both fidelity and mutational bias. WNV fidelity variants demonstrated host-specific alterations to replicative fitness in vitro, with modest attenuation in mosquito but not vertebrate cell culture. Experimental infections of colonized and field populations of Cx. quinquefaciatus demonstrated that WNV fidelity alterations are associated with a significantly impaired capacity to establish viable infections in mosquitoes. Taken together, these studies (i) demonstrate the importance of allosteric interactions in regulating mutation rates, (ii) establish that mutational spectra can be both sequence and strain-dependent, and (iii) display the profound phenotypic consequences associated with altered replication complex function of flaviviruses.

West Nile virus (WNV) is the most geographically widespread arthropod-borne virus (arbovirus) in the world. Like most arboviruses, WNV is a RNA virus which is highly mutable and exists in nature as genetically diverse mutant swarms. Although many recent studies have investigated the relationship between virus mutation rate and viral fitness, this had not previously been determined for WNV or other flaviviruses. We identified WNV mutations associated with variation in mutation rate using cell culture passage in the presence of a mutagen and engineered these mutations into an infectious WNV clone in order to investigate the causes and consequences of altered fidelity. Our results demonstrate that interactions among proteins which comprise the WNV replication complex can significantly alter both the extent and types of mutations that occur. In addition, we show that both increasing and decreasing WNV fidelity has host-specific effects on replication in cell culture and is associated with nearly complete ablation of WNV infection in mosquito vectors. These results have significant implications for our understanding of arbovirus evolution, replication complex function and arboviral fitness in mosquitoes, and identify important targets to study the determinants and mechanisms of vector competence and arbovirus fidelity.

DEVELOPMENTALLY REGULATED PLASMA MEMBRANE PROTEIN of Nicotiana benthamiana contributes to potyvirus movement and transports to plasmodesmata via the early secretory pathway and the actomyosin system

The intercellular movement of plant viruses requires both viral and host proteins. Previous studies have demonstrated that the frame-shift protein P3N-PIPO (for the protein encoded by the open reading frame [ORF] containing 5'-terminus of P3 and a +2 frame-shift ORF called Pretty Interesting Potyviridae ORF and embedded in the P3) and CYLINDRICAL INCLUSION (CI) proteins were required for potyvirus cell-to-cell movement. Here, we provide genetic evidence showing that a Tobacco vein banding mosaic virus (TVBMV; genus Potyvirus) mutant carrying a truncated PIPO domain of 58 amino acid residues could move between cells and induce systemic infection in Nicotiana benthamiana plants; mutants carrying a PIPO domain of seven, 20, or 43 amino acid residues failed to move between cells and cause systemic infection in this host plant. Interestingly, the movement-defective mutants produced progeny that eliminated the previously introduced stop codons and thus restored their systemic movement ability. We also present evidence showing that a developmentally regulated plasma membrane protein of N. benthamiana (referred to as NbDREPP) interacted with both P3N-PIPO and CI of the movement-competent TVBMV. The knockdown of NbDREPP gene expression in N. benthamiana impeded the cell-to-cell movement of TVBMV. NbDREPP was shown to colocalize with TVBMV P3N-PIPO and CI at plasmodesmata (PD) and traffic to PD via the early secretory pathway and the actomyosin motility system. We also show that myosin XI-2 is specially required for transporting NbDREPP to PD. In conclusion, NbDREPP is a key host protein within the early secretory pathway and the actomyosin motility system that interacts with two movement proteins and influences virus movement.

Getting to Know Viral Evolutionary Strategies: Towards the Next Generation of Quasispecies Models

Viral populations are formed by complex ensembles of genomes with broad phenotypic diversity. The adaptive strategies deployed by these ensembles are multiple and often cannot be predicted a priori. Our understanding of viral dynamics is mostly based on two kinds of empirical approaches: one directed towards characterizing molecular changes underlying fitness changes and another focused on population-level responses. Simultaneously, theoretical efforts are directed towards developing a formal picture of viral evolution by means of more realistic fitness landscapes and reliable population dynamics models. New technologies, chiefly the use of next-generation sequencing and related tools, are opening avenues connecting the molecular and the population levels. In the near future, we hope to be witnesses of an integration of these still decoupled approaches, leading into more accurate and realistic quasispecies models able to capture robust generalities and endowed with a satisfactory predictive power.

Genetic diversity and evolution of satellite RNAs associated with the bamboo mosaic virus

Satellite RNAs (satRNAs) are subviral agents that depend on cognate helper viruses for genome replication and encapsidation. Their negative impacts on helper viruses have been exploited to control plant viral diseases. SatBaMV is a commonly found satRNA associated with Bamboo mosaic virus (BaMV) that infects diverse bamboo species in the field. To investigate the genetic diversity and evolution of satRNAs, we examined seven satBaMV populations derived from five bamboo species and cultivars from Taiwan, China, and India and one from the greenhouse. We found 3 distinct clades among the seven populations. Clade I is consisted of all satBaMV isolates, except for those from Dendrocalamus latiflorus in Taiwan and Bambusa vulgaris in India, which belong to Clades II and III, respectively. Interestingly, nucleotide diversity was lower for Clade I than II and III. However, the nucleotide diversity did not seem to depend on bamboo species or geographic location. Our population genetic analyses revealed the presence of excessive low-frequency polymorphic sites, which suggests that the satBaMV population was under purifying selection and/or population expansion. Further analysis of P20, the only satBaMV gene that encodes a non-structural protein involved in the long-distance movement of satBaMV, showed evidence of purifying selection. Taken together, our results suggest that purifying selection against defective P20 protein is responsible at least in part for the evolution of the satBaMV genome.

RNA virus evolution at variable error rate

ABSTRACT: 

RNA viruses replicate their genomes with very high error rates, which leads to the generation of a large genetic diversity that makes them highly adaptable to most environmental pressures, including antiviral drugs and immune responses. However, since most mutations are deleterious, an excess of errors can be very negative for RNA viruses, entailing that error rates must be finely regulated. Currently, the manipulation of the error rate is emerging as a promising antiviral therapy that could minimize the problem of virus adaptation to classical treatments. This review provides a detailed analysis of the different outcomes that can result from the variation of the error rate in RNA viruses, on the basis of the more relevant findings obtained in experimental studies.

Variability in the mutation rates of RNA viruses

ABSTRACT: 

It is well established that RNA viruses show extremely high mutation rates, but less attention has been paid to the fact that their mutation rates also vary strongly, from 10-6 to 10-4 substitutions per nucleotide per cell infection. The causes explaining this variability are still poorly understood, but candidate factors are the viral genome size and polarity, host-specific gene expression patterns, or the intracellular environment. Differences between animal and plant viruses, or between arthropod-borne and directly transmitted viruses have also been postulated. Finally, RNA viruses may be able to regulate the rate at which new mutations spread in the population by modifying features of the viral infection cycle, such as lysis time.

Variation in RNA virus mutation rates across host cells

It is well established that RNA viruses exhibit higher rates of spontaneous mutation than DNA viruses and microorganisms. However, their mutation rates vary amply, from 10⁻⁶ to 10⁻⁴ substitutions per nucleotide per round of copying (s/n/r) and the causes of this variability remain poorly understood. In addition to differences in intrinsic fidelity or error correction capability, viral mutation rates may be dependent on host factors. Here, we assessed the effect of the cellular environment on the rate of spontaneous mutation of the vesicular stomatitis virus (VSV), which has a broad host range and cell tropism. Luria-Delbrück fluctuation tests and sequencing showed that VSV mutated similarly in baby hamster kidney, murine embryonic fibroblasts, colon cancer, and neuroblastoma cells (approx. 10⁻⁵ s/n/r). Cell immortalization through p53 inactivation and oxygen levels (1–21%) did not have a significant impact on viral replication fidelity. This shows that previously published mutation rates can be considered reliable despite being based on a narrow and artificial set of laboratory conditions. Interestingly, we also found that VSV mutated approximately four times more slowly in various insect cells compared with mammalian cells. This may contribute to explaining the relatively slow evolution of VSV and other arthropod-borne viruses in nature.

RNA viruses show high rates of spontaneous mutation, a feature that profoundly influences viral evolution, disease emergence, the appearance of drug resistances, and vaccine efficacy. However, RNA virus mutation rates vary substantially and the factors determining this variability remain poorly understood. Here, we investigated the effects of host factors on viral replication fidelity by measuring the viral mutation rate in different cell types and under various culturing conditions. To carry out these experiments we chose the vesicular stomatitis virus (VSV), an insect-transmitted mammalian RNA virus with an extremely wide cellular and host tropism. We found that the VSV replication machinery was robust to changes in cellular physiology driven by cell immortalization or shifts in temperature and oxygen levels. In contrast, VSV mutated significantly more slowly in insect cells than in mammalian cells, a finding may help us to understand why arthropod-borne viruses tend to evolve more slowly than directly transmitted viruses in nature.

Fixation of emerging interviral recombinants in cucumber mosaic virus populations

Interstrain recombinants were observed in the progenies of the Cucumber mosaic virus (CMV) reassortant L(1)L(2)F(3) containing RNAs 1 and 2 from LS-CMV and RNA 3 from Fny-CMV. We characterized these recombinants, and we found that their fixation was controlled by the nature of the replicating RNAs 1 and 2. We demonstrate that the 2b gene partially affects this fixation process, but only in the context of homologous RNAs 1 and 2.

Adaptive evolution by recombination is not associated with increased mutation rates in Maize streak virus

BACKGROUND

Single-stranded (ss) DNA viruses in the family Geminiviridae are proving to be very useful in real-time evolution studies. The high mutation rate of geminiviruses and other ssDNA viruses is somewhat mysterious in that their DNA genomes are replicated in host nuclei by high fidelity host polymerases. Although strand specific mutation biases observed in virus species from the geminivirus genus Mastrevirus indicate that the high mutation rates in viruses in this genus may be due to mutational processes that operate specifically on ssDNA, it is currently unknown whether viruses from other genera display similar strand specific mutation biases. Also, geminivirus genomes frequently recombine with one another and an alternative cause of their high mutation rates could be that the recombination process is either directly mutagenic or produces a selective environment in which the survival of mutants is favoured. To investigate whether there is an association between recombination and increased basal mutation rates or increased degrees of selection favoring the survival of mutations, we compared the mutation dynamics of the MSV-MatA and MSV-VW field isolates of Maize streak virus (MSV; Mastrevirus), with both a laboratory constructed MSV recombinant, and MSV recombinants closely resembling MSV-MatA. To determine whether strand specific mutation biases are a general characteristic of geminivirus evolution we compared mutation spectra arising during these MSV experiments with those arising during similar experiments involving the geminivirus Tomato yellow leaf curl virus (Begomovirus genus).

RESULTS

Although both the genomic distribution of mutations and the occurrence of various convergent mutations at specific genomic sites indicated that either mutation hotspots or selection for adaptive mutations might elevate observed mutation rates in MSV, we found no association between recombination and mutation rates. Importantly, when comparing the mutation spectra of MSV and TYLCV we observed similar strand specific mutation biases arising predominantly from imbalances in the complementary mutations G → T: C → A.

CONCLUSIONS

While our results suggest that recombination does not strongly influence mutation rates in MSV, they indicate that high geminivirus mutation rates are at least partially attributable to increased susceptibility of all geminivirus genomes to oxidative damage while in a single stranded state.

Mapping viral functional domains for genetic diversity in plants

Cucumber mosaic virus (CMV) comprises numerous isolates with various levels of in-host diversity. Subgroup-distinctive features of the Fny and LS strains provided us with a platform to genetically map the viral control elements for genetic variation in planta. We found that both RNAs 1 and 2 controlled levels of genetic diversity, and further fine mapping revealed that the control elements of mutation frequency reside within the first 596 amino acids (aa) of RNA 1. The 2a/2b overlapping region of the 2a protein also contributed to control of viral genetic variation. Furthermore, the 3' nontranslated region (NTR) of RNA 3 constituted a hot spot of polymorphism, where the majority of fixed mutations found in the population were clustered. The 2b gene of CMV, a viral suppressor of gene silencing, controls the abundance of the fixed mutants in the viral population via a host-dependent mechanism.

The evolutionary genetics of emerging plant RNA viruses

Over the years, agriculture across the world has been compromised by a succession of devastating epidemics caused by new viruses that spilled over from reservoir species or by new variants of classic viruses that acquired new virulence factors or changed their epidemiological patterns. Viral emergence is usually associated with ecological change or with agronomical practices bringing together reservoirs and crop species. The complete picture is, however, much more complex, and results from an evolutionary process in which the main players are ecological factors, viruses' genetic plasticity, and host factors required for virus replication, all mixed with a good measure of stochasticity. The present review puts emergence of plant RNA viruses into the framework of evolutionary genetics, stressing that viral emergence begins with a stochastic process that involves the transmission of a preexisting viral strain into a new host species, followed by adaptation to the new host.

Changes in population dynamics in mutualistic versus pathogenic viruses

Although generally regarded as pathogens, viruses can also be mutualists. A number of examples of extreme mutualism (i.e., symbiogenesis) have been well studied. Other examples of mutualism are less common, but this is likely because viruses have rarely been thought of as having any beneficial effects on their hosts. The effect of mutualism on the population dynamics of viruses is a topic that has not been addressed experimentally. However, the potential for understanding mutualism and how a virus might become a mutualist may be elucidated by understanding these dynamics.

Mutant spectra in virus behavior

RNA viruses replicate as complex mutant spectra, also termed ‘mutant clouds’, known as viral quasispecies. While this is a widely observed viral population structure, it is less known that a number of biologically relevant features of this important group of viral pathogens depend on (or are strongly influenced by) the complexity and composition of mutant spectra. Among them, fitness increase or decrease depending on intrapopulation complementation or interference, selection triggered by memory genomes, pathogenic potential of viruses, disease evolution and the response to antiviral treatments. Quasispecies represent the recognition of complex behavior in viruses, and it is an oversimplification to equate such a population structure with the classic polymorphism of population biology. Darwinian principles acting on genome collectivities that replicate with high error rates provide a unique population structure prone to flexible and largely unpredictable behavior.

Genetic bottlenecks during systemic movement of Cucumber mosaic virus vary in different host plants

Genetic bottlenecks are stochastic events that narrow variation in a population. We compared bottlenecks during the systemic infection of Cucumber mosaic virus (CMV) in four host plants. We mechanically inoculated an artificial population of twelve CMV mutants to young leaves of tomato, pepper, Nicotiana benthamiana, and squash. The inoculated leaves and primary and secondary systemically infected leaves were sampled at 2, 10, and 15 days post-inoculation. All twelve mutants were detected in all of the inoculated leaves. The number of mutants recovered from the systemically infected leaves of all host species was reduced significantly, indicating bottlenecks in systemic movement. The recovery frequencies of a few of the mutants were significantly different in each host probably due to host-specific selective forces. These results have implications for the differences in virus population variation that is seen in different host plants.

Tomato bushy stunt virus recombination guided by introduced microRNA target sequences

Previously we described Tomato bushy stunt virus (TBSV) vectors, which retained their capsid protein gene and were engineered with magnesium chelatase (ChlH) and phytoene desaturase (PDS) gene sequences from Nicotiana benthamiana. Upon plant infection, these vectors eventually lost the inserted sequences, presumably as a result of recombination. Here, we modified the same vectors to also contain the plant miR171 or miR159 target sequences immediately 3' of the silencing inserts. We inoculated N. benthamiana plants and sequenced recombinant RNAs recovered from noninoculated upper leaves. We found that while some of the recombinant RNAs retained the microRNA (miRNA) target sites, most retained only the 3' 10 and 13 nucleotides of the two original plant miRNA target sequences, indicating in planta miRNA-guided RNA-induced silencing complex cleavage of the recombinant TBSV RNAs. In addition, recovered RNAs also contained various fragments of the original sequence (ChlH and PDS) upstream of the miRNA cleavage site, suggesting that the 3' portion of the miRNA-cleaved TBSV RNAs served as a template for negative-strand RNA synthesis by the TBSV RNA-dependent RNA polymerase (RdRp), followed by template switching by the RdRp and continued RNA synthesis resulting in loss of nonessential nucleotides.

Effect of citrus hosts on the generation, maintenance and evolutionary fate of genetic variability of citrus exocortis viroid

Citrus exocortis viroid (CEVd) populations are composed of closely related haplotypes whose frequencies in the population result from the equilibrium between mutation, selection and genetic drift. The genetic diversity of CEVd populations infecting different citrus hosts was studied by comparing populations recovered from infected trifoliate orange and sour orange seedling trees after 10 years of evolution, with the ancestral population maintained for the same period in the original host, Etrog citron. Furthermore, populations isolated from these trifoliate orange and sour orange trees were transmitted back to Etrog citron plants and the evolution of their mutant spectra was studied. The results indicate that (i) the amount and composition of the within-plant genetic diversity generated varies between these two hosts and is markedly different from that which is characteristic of the original Etrog citron host and (ii) the genetic diversity found after transmitting back to Etrog citron is indistinguishable from that which is characteristic of the ancestral Etrog citron population, regardless of the citrus plant from which the evolved populations were isolated. The relationship between the CEVd populations from Etrog citron and trifoliate orange, both sensitive hosts, and those from sour orange, which is a tolerant host, is discussed.

Upper-limit mutation rate estimation for a plant RNA virus

It is generally accepted that mutation rates of RNA viruses are inherently high due to the lack of proofreading mechanisms. However, direct estimates of mutation rate are surprisingly scarce, in particular for plant viruses. Here, based on the analysis of in vivo mutation frequencies in tobacco etch virus, we calculate an upper-bound mutation rate estimation of 3x10(-5) per site and per round of replication; a value which turns out to be undistinguishable from the methodological error. Nonetheless, the value is barely on the lower side of the range accepted for RNA viruses, although in good agreement with the only direct estimate obtained for other plant viruses. These observations suggest that, perhaps, differences in the selective pressures operating during plant virus evolution may have driven their mutation rates towards values lower than those characteristic of other RNA viruses infecting bacteria or animals.

Phylogenetic analysis reveals the emergence, evolution and dispersal of carnivore parvoviruses

Canine parvovirus (CPV), first recognized as an emerging virus of dogs in 1978, resulted from a successful cross-species transmission. CPV emerged from the endemic feline panleukopenia virus (FPV), or from a closely related parvovirus of another host. Here we refine our current understanding of the evolution and population dynamics of FPV and CPV. By analysing nearly full-length viral sequences we show that the majority of substitutions distinguishing CPV from FPV are located in the capsid protein gene, and that this gene is under positive selection in CPV, resulting in a significantly elevated rate of molecular evolution. This provides strong phylogenetic evidence for a prominent role of the viral capsid in host adaptation. In addition, an analysis of the population dynamics of more recent CPV reveals, on a global scale, a strongly spatially subdivided CPV population with little viral movement among countries and a relatively constant population size. Such limited viral migration contrasts with the global spread of the virus observed during the early phase of the CPV pandemic, but corresponds to the more endemic nature of current CPV infections.

Changes in the intraisolate genetic structure of Beet necrotic yellow vein virus populations associated with plant resistance breakdown

The causal agent of rhizomania disease, Beet necrotic yellow vein virus (BNYVV), typically produces asymptomatic root-limited infections in sugar beets (Beta vulgaris) carrying the Rz1-allele. Unfortunately, this dominant resistance has been recently overcome. Multiple cDNA clones of the viral pathogenic determinant p25, derived from populations infecting susceptible or resistant plants, were sequenced to identify host effects on the viral population structure. Populations isolated from compatible plant-virus interactions (susceptible plant-wild type virus and resistant plant-resistant breaking variants) were large and relatively homogeneous, whereas those from the incompatible interaction (resistant plant-avirulent type virus) were small and highly heterogeneous. All populations from susceptible plants had the same dominant haplotype, whereas those from resistant cultivars had a different haplotype surrounded by a spectrum of mutants. Selection and diversification analyses suggest an evolutionary trajectory of BNYVV with positive selection for changes required to overcome resistance, followed by elimination of hitchhiking mutations through purifying selection.

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.