Distribution of spontaneous mutants and inferences about the replication mode of the RNA bacteriophage phi6

Microbe Profile: Bacteriophage ϕ6: a model for segmented RNA viruses and the evolutionary consequences of viral 'sex'

Bacteriophage ϕ6 is a segmented dsRNA virus with a lipid envelope, which are unusual traits in bacterial viruses but common in eukaryotic viruses. This uniqueness allowed ϕ6 and its Pseudomonad hosts to serve as a molecular model for RNA genetics, mutation, replication, packaging, and reassortment in both bacterial and eukaryotic viruses. However, an additional uniqueness of ϕ6, created by its high mutation rate, was its use as an experimental system to study key questions such as the evolution of sex (segment reassortment), host-pathogen interactions, mutational load, rates of adaptation, genetic and phenotypic complexity, and game theory.

Developing Phage Therapy That Overcomes the Evolution of Bacterial Resistance

The global rise of antibiotic resistance in bacterial pathogens and the waning efficacy of antibiotics urge consideration of alternative antimicrobial strategies. Phage therapy is a classic approach where bacteriophages (bacteria-specific viruses) are used against bacterial infections, with many recent successes in personalized medicine treatment of intractable infections. However, a perpetual challenge for developing generalized phage therapy is the expectation that viruses will exert selection for target bacteria to deploy defenses against virus attack, causing evolution of phage resistance during patient treatment. Here we review the two main complementary strategies for mitigating bacterial resistance in phage therapy: minimizing the ability for bacterial populations to evolve phage resistance and driving (steering) evolution of phage-resistant bacteria toward clinically favorable outcomes. We discuss future research directions that might further address the phage-resistance problem, to foster widespread development and deployment of therapeutic phage strategies that outsmart evolved bacterial resistance in clinical settings.

Single-cell mutation rate of turnip crinkle virus (-)-strand replication intermediates

Viruses with single-stranded, positive-sense (+) RNA genomes incur high numbers of errors during replication, thereby creating diversified genome populations from which new, better adapted viral variants can emerge. However, a definitive error rate is known for a relatively few (+) RNA plant viruses, due to challenges to account for perturbations caused by natural selection and/or experimental set-ups. To address these challenges, we developed a new approach that exclusively profiled errors in the (-)-strand replication intermediates of turnip crinkle virus (TCV), in singly infected cells. A series of controls and safeguards were devised to ensure errors inherent to the experimental process were accounted for. This approach permitted the estimation of a TCV error rate of 8.47 X 10-5 substitution per nucleotide site per cell infection. Importantly, the characteristic error distribution pattern among the 50 copies of 2,363-base-pair cDNA fragments predicted that nearly all TCV (-) strands were products of one replication cycle per cell. Furthermore, some of the errors probably elevated error frequencies by lowering the fidelity of TCV RNA-dependent RNA polymerase, and/or permitting occasional re-replication of progeny genomes. In summary, by profiling errors in TCV (-)-strand intermediates incurred during replication in single cells, this study provided strong support for a stamping machine mode of replication employed by a (+) RNA virus.

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.

Plant Virus Adaptation to New Hosts: A Multi-scale Approach

Viruses are studied at each level of biological complexity: from within-cells to ecosystems. The same basic evolutionary forces and principles operate at each level: mutation and recombination, selection, genetic drift, migration, and adaptive trade-offs. Great efforts have been put into understanding each level in great detail, hoping to predict the dynamics of viral population, prevent virus emergence, and manage their spread and virulence. Unfortunately, we are still far from this. To achieve these ambitious goals, we advocate for an integrative perspective of virus evolution. Focusing in plant viruses, we illustrate the pervasiveness of the above-mentioned principles. Beginning at the within-cell level, we describe replication modes, infection bottlenecks, and cellular contagion rates. Next, we move up to the colonization of distal tissues, discussing the fundamental role of random events. Then, we jump beyond the individual host and discuss the link between transmission mode and virulence. Finally, at the community level, we discuss properties of virus-plant infection networks. To close this review we propose the multilayer network theory, in which elements at different layers are connected and submit to their own dynamics that feed across layers, resulting in new emerging properties, as a way to integrate information from the different levels.

Emergence of Southern Rice Black-Streaked Dwarf Virus in the Centuries-Old Chinese Yuanyang Agrosystem of Rice Landraces

Southern rice black-streaked dwarf virus (SRBSDV), which causes severe disease symptoms in rice (Oriza sativa L.) has been emerging in the last decade throughout northern Vietnam, southern Japan and southern, central and eastern China. Here we attempt to quantify the prevalence of SRBSDV in the Honghe Hani rice terraces system (HHRTS)-a Chinese 1300-year-old traditional rice production system. We first confirm that genetically diverse rice varieties are still being cultivated in the HHRTS and categorize these varieties into three main genetic clusters, including the modern hybrid varieties group (MH), the Hongyang improved modern variety group (HY) and the traditional indica landraces group (TIL). We also show over a 2-year period that SRBSDV remains prevalent in the HHRTS (20.1% prevalence) and that both the TIL (17.9% prevalence) and the MH varieties (5.1% prevalence) were less affected by SRBSDV than were the HY varieties (30.2% prevalence). Collectively we suggest that SRBSDV isolates are freely moving within the HHRTS and that TIL, HY and MH rice genetic clusters are not being preferentially infected by particular SRBSDV lineages. Given that SRBSDV can cause 30-50% rice yield losses, our study emphasizes both the need to better monitor the disease in the HHRTS, and the need to start considering ways to reduce its burden on rice production.

Noise-induced bistability in the quasi-neutral coexistence of viral RNAs under different replication modes

Evolutionary and dynamical investigations into real viral populations indicate that RNA replication can range between the two extremes represented by so-called 'stamping machine replication' (SMR) and 'geometric replication' (GR). The impact of asymmetries in replication for single-stranded (+) sense RNA viruses has been mainly studied with deterministic models. However, viral replication should be better described by including stochasticity, as the cell infection process is typically initiated with a very small number of RNA macromolecules, and thus largely influenced by intrinsic noise. Under appropriate conditions, deterministic theoretical descriptions of viral RNA replication predict a quasi-neutral coexistence scenario, with a line of fixed points involving different strands' equilibrium ratios depending on the initial conditions. Recent research into the quasi-neutral coexistence in two competing populations reveals that stochastic fluctuations fundamentally alter the mean-field scenario, and one of the two species outcompetes the other. In this article, we study this phenomenon for viral RNA replication modes by means of stochastic simulations and a diffusion approximation. Our results reveal that noise has a strong impact on the amplification of viral RNAs, also causing the emergence of noise-induced bistability. We provide analytical criteria for the dominance of (+) sense strands depending on the initial populations on the line of equilibria, which are in agreement with direct stochastic simulation results. The biological implications of this noise-driven mechanism are discussed within the framework of the evolutionary dynamics of RNA viruses with different modes of replication.

STABILITY OF A STOCHASTICALLY PERTURBED MODEL OF INTRACELLULAR SINGLE-STRANDED RNA VIRUS REPLICATION

Compared to the replication of double-stranded RNA and DNA viruses, the replication of single-stranded viruses requires the production of a number of intermediate strands that serve as templates for the synthesis of genomic-sense strands. Two theoretical extreme mechanisms for replication for such single-stranded viruses have been proposed; one extreme being represented by the so-called linear stamping machine and the opposite extreme by the exponential growth. Of course, real systems are more complex and examples have been described in which a combination of such extreme mechanisms can also occur: a fraction of the produced progeny resulting from a stamping-machine type of replication that uses the parental genome as template, whereas other fraction of the progeny results from the replication of other progeny genomes. Martínez et al. 1 , Sardanyés et al. 2 and Fornés et al. 3 suggested and analyzed a deterministic model of single-stranded RNA (ssRNA) virus intracellular replication that incorporated variability in the replication mechanisms.

To explore how stochasticity can affect this mixed-model principal properties, in this paper, we consider the stability of a stochastically perturbed model of ssRNA virus replication within a cell. Using the direct Lyapunov method, we found sufficient conditions for the stability in probability of equilibrium states for this model. This result confirms that this heterogeneous model of single-stranded RNA virus replication is stable with respect to stochastic perturbations of the environment.

Viral replication modes in single-peak fitness landscapes: A dynamical systems analysis

Positive-sense, single-stranded RNA viruses are important pathogens infecting almost all types of organisms. Experimental evidence from distributions of mutations and from viral RNA amplification suggest that these pathogens may follow different RNA replication modes, ranging from the stamping machine replication (SMR) to the geometric replication (GR) mode. Although previous theoretical work has focused on the evolutionary dynamics of RNA viruses amplifying their genomes with different strategies, little is known in terms of the bifurcations and transitions involving the so-called error threshold (mutation-induced dominance of mutants) and lethal mutagenesis (extinction of all sequences due to mutation accumulation and demographic stochasticity). Here we analyze a dynamical system describing the intracellular amplification of viral RNA genomes evolving on a single-peak fitness landscape focusing on three cases considering neutral, deleterious, and lethal mutants. We analytically derive the critical mutation rates causing lethal mutagenesis and error threshold, governed by transcritical bifurcations that depend on parameters α (parameter introducing the mode of replication), replicative fitness of mutants (k₁), and on the spontaneous degradation rates of the sequences (ϵ). Our results relate the error catastrophe with lethal mutagenesis in a model with continuous populations of viral genomes. The former case involves dominance of the mutant sequences, while the latter, a deterministic extinction of the viral RNAs during replication due to increased mutation. For the lethal case the critical mutation rate involving lethal mutagenesis is μ_c=1-ɛ/α. Here, the SMR involves lower critical mutation rates, being the system more robust to lethal mutagenesis replicating closer to the GR mode. This result is also found for the neutral and deleterious cases, but for these later cases lethal mutagenesis can shift to the error threshold once the replication mode surpasses a threshold given by α=ϵ/k₁.

Modeling the probability distribution of the bacterial burst size via a game-theoretic approach

Based on previous studies, empirical distribution of the bacterial burst size varies even in a population of isogenic bacteria. Since bacteriophage progenies increase linearly with time, it is the lysis time variation that results in the bacterial burst size variations. Here, the burst size variation is computationally modeled by considering the lysis time decisions as a game. Each player in the game is a bacteriophage that has initially infected and lysed its host bacterium. Also, the payoff of each burst size strategy is the average number of bacteria that are solely infected by the bacteriophage progenies after lysis. For calculating the payoffs, a new version of ball and bin model with time dependent occupation probabilities (TDOP) is proposed. We show that Nash equilibrium occurs for a range of mixed burst size strategies that are chosen and played by bacteriophages, stochastically. Moreover, it is concluded that the burst size variations arise from choosing mixed lysis strategies by each player. By choosing the lysis time and also the burst size stochastically, the released bacteriophage progenies infect a portion of host bacteria in environment and avoid extinction. The probability distribution of the mixed burst size strategies is also identified.

Tobamoviruses as Models for the Study of Virus Evolution

The study of tobacco mosaic virus and other tobamovirus species has greatly contributed to the development of all areas of virology, including virus evolution. Research with tobamoviruses has been pioneer, or particularly significant, in all major areas of research in this field, including: the characterization of the genetic diversity of virus populations, the mechanisms and rates of generation of genetic diversity, the analysis of the genetic structure of virus populations and of the factors that shape it, the adaptation of viruses to hosts and the evolution of host range, and the evolution of virus taxa and of virus-host interactions. Many of these continue to be hot topics in evolutionary biology, or have been identified recently as such, including (i) host-range evolution, (ii) predicting the overcoming of resistance in crops, (iii) trade-offs between virus life-history traits in virus evolution, and (iv) the codivergence of viruses and hosts at different taxonomical and spatial scales. Tobamoviruses may be particularly appropriate to address these topics with plant viruses, as they provide convenient experimental systems, and as the detailed knowledge on their molecular and structural biology allows the analysis of the mechanisms behind evolutionary processes. Also, the extensive information on parameters related to infection dynamics and population structure may facilitate the development of realistic models to predict virus evolution. Certainly, tobamoviruses will continue to be favorite system for the study of virus evolution.

Cystoviral RNA-directed RNA polymerases: Regulation of RNA synthesis on multiple time and length scales

•

Role of the RNA polymerase in the cystoviral life-cycle.

•

Spatio-temporal regulation of RNA synthesis in cystoviruses.

•

Emerging role of conformational dynamics in polymerase function.

P2, an RNA-directed RNA polymerase (RdRP), is encoded on the largest of the three segments of the double-stranded RNA genome of cystoviruses. P2 performs the dual tasks of replication and transcription de novo on single-stranded RNA templates, and plays a critical role in the viral life-cycle. Work over the last few decades has yielded a wealth of biochemical and structural information on the functional regulation of P2, on its role in the spatiotemporal regulation of RNA synthesis and its variability across the Cystoviridae family. These range from atomic resolution snapshots of P2 trapped in functionally significant states, in complex with catalytic/structural metal ions, polynucleotide templates and substrate nucleoside triphosphates, to P2 in the context of viral capsids providing structural insight into the assembly of supramolecular complexes and regulatory interactions therein. They include in vitro biochemical studies using P2 purified to homogeneity and in vivo studies utilizing infectious core particles. Recent advances in experimental techniques have also allowed access to the temporal dimension and enabled the characterization of dynamics of P2 on the sub-nanosecond to millisecond timescale through measurements of nuclear spin relaxation in solution and single molecule studies of transcription from seconds to minutes. Below we summarize the most significant results that provide critical insight into the role of P2 in regulating RNA synthesis in cystoviruses.

Stability of A Coevolving Host-parasite System Peaks at Intermediate Productivity

Habitat productivity may affect the stability of consumer-resource systems, through both ecological and evolutionary mechanisms. We hypothesize that coevolving consumer-resource systems show more stable dynamics at intermediate resource availability, while very low-level resource supply cannot support sufficiently large populations of resource and consumer species to avoid stochastic extinction, and extremely resource-rich environments may promote escalatory arms-race-like coevolution that can cause strong fluctuations in species abundance and even extinction of one or both trophic levels. We tested these ideas by carrying out an experimental evolution study with a model bacterium-phage system (Pseudomonas fluorescens SBW25 and its phage SBW25Φ2). Consistent with our hypothesis, this system was most stable at intermediate resource supply (fewer extinction events and smaller magnitude of population fluctuation). In our experiment, the rate of coevolution between bacterial resistance and phage infectivity was correlated with the magnitude of population fluctuation, which may explain the different in stability between levels of resource supply. Crucially, our results are consistent with a suggestion that, among the two major modes of antagonistic coevolution, arms race is more likely than fluctuation selection dynamics to cause extinction events in consumer-resource systems. This study suggests an important role of environment-dependent coevolutionary dynamics for the stability of consumer-resource species systems, therefore highlights the importance to consider contemporaneous evolutionary dynamics when studying the stability of ecosystems, particularly those under environmental changes.

Mechanisms of viral mutation

The remarkable capacity of some viruses to adapt to new hosts and environments is highly dependent on their ability to generate de novo diversity in a short period of time. Rates of spontaneous mutation vary amply among viruses. RNA viruses mutate faster than DNA viruses, single-stranded viruses mutate faster than double-strand virus, and genome size appears to correlate negatively with mutation rate. Viral mutation rates are modulated at different levels, including polymerase fidelity, sequence context, template secondary structure, cellular microenvironment, replication mechanisms, proofreading, and access to post-replicative repair. Additionally, massive numbers of mutations can be introduced by some virus-encoded diversity-generating elements, as well as by host-encoded cytidine/adenine deaminases. Our current knowledge of viral mutation rates indicates that viral genetic diversity is determined by multiple virus- and host-dependent processes, and that viral mutation rates can evolve in response to specific selective pressures.

Single-Cell Analysis of RNA Virus Infection Identifies Multiple Genetically Diverse Viral Genomes within Single Infectious Units

Genetic diversity enables a virus to colonize novel hosts, evade immunity, and evolve drug resistance. However, viral diversity is typically assessed at the population level. Given the existence of cell-to-cell variation, it is critical to understand viral genetic structure at the single-cell level. By combining single-cell isolation with ultra-deep sequencing, we characterized the genetic structure and diversity of a RNA virus shortly after single-cell bottlenecks. Full-length sequences from 881 viral plaques derived from 90 individual cells reveal that sequence variants pre-existing in different viral genomes can be co-transmitted within the same infectious unit to individual cells. Further, the rate of spontaneous virus mutation varies across individual cells, and early production of diversity depends on the viral yield of the very first infected cell. These results unravel genetic and structural features of a virus at the single-cell level, with implications for viral diversity and evolution.

•

Genetic diversity is found within viral units infecting an individual cell

•

Viral particles within the same infectious unit interact genetically and functionally

•

The production of viral genetic diversity is marked by widespread cell-to-cell variation

•

Single-cell virus yields determine the early production viral genetic diversity

Genetically Determined Variation in Lysis Time Variance in the Bacteriophage φX174

Researchers in evolutionary genetics recently have recognized an exciting opportunity in decomposing beneficial mutations into their proximal, mechanistic determinants. The application of methods and concepts from molecular biology and life history theory to studies of lytic bacteriophages (phages) has allowed them to understand how natural selection sees mutations influencing life history. This work motivated the research presented here, in which we explored whether, under consistent experimental conditions, small differences in the genome of bacteriophage φX174 could lead to altered life history phenotypes among a panel of eight genetically distinct clones. We assessed the clones' phenotypes by applying a novel statistical framework to the results of a serially sampled parallel infection assay, in which we simultaneously inoculated each of a large number of replicate host volumes with ∼1 phage particle. We sequentially plated the volumes over the course of infection and counted the plaques that formed after incubation. These counts served as a proxy for the number of phage particles in a single volume as a function of time. From repeated assays, we inferred significant, genetically determined heterogeneity in lysis time and burst size, including lysis time variance. These findings are interesting in light of the genetic and phenotypic constraints on the single-protein lysis mechanism of φX174. We speculate briefly on the mechanisms underlying our results, and we discuss the potential importance of lysis time variance in viral evolution.

Experimentally guided models reveal replication principles that shape the mutation distribution of RNA viruses

Life history theory posits that the sequence and timing of events in an organism's lifespan are fine-tuned by evolution to maximize the production of viable offspring. In a virus, a life history strategy is largely manifested in its replication mode. Here, we develop a stochastic mathematical model to infer the replication mode shaping the structure and mutation distribution of a poliovirus population in an intact single infected cell. We measure production of RNA and poliovirus particles through the infection cycle, and use these data to infer the parameters of our model. We find that on average the viral progeny produced from each cell are approximately five generations removed from the infecting virus. Multiple generations within a single cell infection provide opportunities for significant accumulation of mutations per viral genome and for intracellular selection.

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

Viruses with genetic information made up of molecules of RNA can multiply quickly, but not very accurately. This means that many errors, or mutations, occur when the RNA is copied to create new viruses. The advantage of this rapid, but mistake-filled, RNA replication process is that some of the mutations will be beneficial to the virus. This allows viruses to rapidly evolve, for example, to develop resistance against drugs.

The poliovirus is an RNA virus that can cause paralysis and death in humans. To prevent such infections, scientists have extensively studied the poliovirus and have developed effective vaccines against it that have eliminated the virus from all but a few countries. Because so much is known about the poliovirus and because it has a very simple structure, scientists continue to use the poliovirus as a model to study virus behavior.

One unknown aspect of the poliovirus' behavior is how it replicates after invading a cell. Are all of its RNA copies made from the original viral RNA that first infected the cell, in what is known as a ‘stamping machine’ model? Or do the new copies of the RNA also get copied themselves in a ‘geometric replication mode’ that increases the likelihood of mutations and enables the virus to evolve more rapidly?

Viral RNA molecules are copied by one of the virus's own proteins and so before the viral RNA can be replicated, it must first be translated to form viral proteins. When and where replication begins depends on the concentration of translated proteins around the RNA and so replication tends to begin in particular areas of the cell at different times. Schulte, Draghi et al. used mathematical modeling to create computer simulations of the number of polioviruses in a cell that take into account these time and space constraints. By including random elements in the model, the simulated behavior more accurately follows experimentally recorded data than previously used models.

The results of the model led Schulte, Draghi et al. to conclude that the poliovirus replicates by the ‘geometric mode’; as new copies of the poliovirus RNA are made, each copy goes on to make more copies. This means that in a single infected cell there are multiple generations of RNA, and each generation may undergo distinct mutations that are passed on to the next set of RNA copies. In fact, Schulte, Draghi et al. found that the average virus released from an infected cell is the great-great-great-granddaughter of the original virus that infected the cell. With so many different generations of virus coexisting in a cell, there are a lot of opportunities for new genetic combinations to occur and for viruses to evolve new abilities.

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

New enveloped dsRNA phage from freshwater habitat

Cystoviridae is a family of bacteriophages with a tri-segmented dsRNA genome enclosed in a tri-layered virion structure. Here, we present a new putative member of the Cystoviridae family, bacteriophage ϕNN. ϕNN was isolated from a Finnish lake in contrast to the previously identified cystoviruses, which originate from various legume samples collected in the USA. The nucleotide sequence of the virus reveals a strong genetic similarity (~80 % for the L-segments, ~55 % for the M-segments and ~84 % for the S-segments) to Pseudomonas phage ϕ6, the type member of the virus family. However, the relationship between ϕNN and other cystoviruses is more distant. In general, proteins located in the internal parts of the virion were more conserved than those exposed on the virion surface, a phenomenon previously reported among eukaryotic dsRNA viruses. Structural models of several putative ϕNN proteins propose that cystoviral structures are highly conserved.

Frequency and fitness consequences of bacteriophage φ6 host range mutations

Viruses readily mutate and gain the ability to infect novel hosts, but few data are available regarding the number of possible host range-expanding mutations allowing infection of any given novel host, and the fitness consequences of these mutations on original and novel hosts. To gain insight into the process of host range expansion, we isolated and sequenced 69 independent mutants of the dsRNA bacteriophage Φ6 able to infect the novel host, Pseudomonas pseudoalcaligenes. In total, we found at least 17 unique suites of mutations among these 69 mutants. We assayed fitness for 13 of 17 mutant genotypes on P. pseudoalcaligenes and the standard laboratory host, P. phaseolicola. Mutants exhibited significantly lower fitnesses on P. pseudoalcaligenes compared to P. phaseolicola. Furthermore, 12 of the 13 assayed mutants showed reduced fitness on P. phaseolicola compared to wildtype Φ6, confirming the prevalence of antagonistic pleiotropy during host range expansion. Further experiments revealed that the mechanistic basis of these fitness differences was likely variation in host attachment ability. In addition, using computational protein modeling, we show that host-range expanding mutations occurred in hotspots on the surface of the phage's host attachment protein opposite a putative hydrophobic anchoring domain.

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.

Correlation between mutation rate and genome size in riboviruses: mutation rate of bacteriophage Qβ

Genome sizes and mutation rates covary across all domains of life. In unicellular organisms and DNA viruses, they show an inverse relationship known as Drake’s rule. However, it is still unclear whether a similar relationship exists between genome sizes and mutation rates in RNA genomes. Coronaviruses, the RNA viruses with the largest genomes (∼30 kb), encode a proofreading 3′ exonuclease that allows them to increase replication fidelity. However, it is unknown whether, conversely, the RNA viruses with the smallest genomes tend to show particularly high mutation rates. To test this, we measured the mutation rate of bacteriophage Qβ, a 4.2-kb levivirus. Amber reversion-based Luria–Delbrück fluctuation tests combined with mutant sequencing gave an estimate of 1.4 × 10⁻⁴ substitutions per nucleotide per round of copying, the highest mutation rate reported for any virus using this method. This estimate was confirmed using a direct plaque sequencing approach and after reanalysis of previously published estimates for this phage. Comparison with other riboviruses (all RNA viruses except retroviruses) provided statistical support for a negative correlation between mutation rates and genome sizes. We suggest that the mutation rates of RNA viruses might be optimized for maximal adaptability and that the value of this optimum may in turn depend inversely on genome size.

Intergenerational phenotypic mixing in viral evolution

Viral particles (virions) are made of genomic material packaged with proteins, drawn from the pool of proteins in the parent cell. It is well known that when virion concentrations are high, cells can be coinfected with multiple viral strains that can complement each other. Viral genomes can then interact with proteins derived from different strains, in a phenomenon known as phenotypic mixing. But phenotypic mixing is actually far more common: viruses mutate very often, and each time a mutation occurs, the parent cell contains different types of viral genomes. Due to phenotypic mixing, changes in viral phenotypes can be shifted by a generation from the mutations that cause them. In the regime of evolutionary invasion and escape, when mutations are crucial for the virus to survive, this timing can have a large influence on the probability of emergence of an adapted strain. Modeling the dynamics of viral evolution in these contexts thus requires attention to the mutational mechanism and the determinants of fitness.

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.

Quantifying the diversification of hepatitis C virus (HCV) during primary infection: estimates of the in vivo mutation rate

Hepatitis C virus (HCV) is present in the host with multiple variants generated by its error prone RNA-dependent RNA polymerase. Little is known about the initial viral diversification and the viral life cycle processes that influence diversity. We studied the diversification of HCV during acute infection in 17 plasma donors, with frequent sampling early in infection. To analyze these data, we developed a new stochastic model of the HCV life cycle. We found that the accumulation of mutations is surprisingly slow: at 30 days, the viral population on average is still 46% identical to its transmitted viral genome. Fitting the model to the sequence data, we estimate the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. To confirm these results we analyzed the frequency of stop codons (N = 10) among all possible non-sense mutation targets (M = 898,335), and found a mutation rate of 2.8–3.2×10⁻⁵, consistent with the estimate from the dynamical model. The slow accumulation of mutations is consistent with slow turnover of infected cells and replication complexes within infected cells. This slow turnover is also inferred from the viral load kinetics. Our estimated mutation rate, which is similar to that of other RNA viruses (e.g., HIV and influenza), is also compatible with the accumulation of substitutions seen in HCV at the population level. Our model identifies the relevant processes (long-lived cells and slow turnover of replication complexes) and parameters involved in determining the rate of HCV diversification.

Hepatitis C virus (HCV) is a RNA virus that infects over 170 million people across the world. It leads to a chronic infection in the majority of people who are infected (>70%). Most people only discover that they are infected long after initial infection. Thus, it is difficult to study the very early events in infection. Here we study 17 individuals during the earliest possible stages of infection, from before the virus is detectable in the plasma to around 35 days post-infection. We focus on understanding the viral kinetics and the diversification of HCV during this acute phase of infection. During chronic infection HCV is present in the host as a swarm of multiple variants generated by its error prone copying. We studied the early diversification of HCV during acute infection using a new mathematical model of HCV replication. We found that after a phase of fast increase in viral load, accompanied by viral diversification, there is a stabilization of viral load and diversity levels. Using our model, we were able to estimate for the first time the HCV mutation rate during acute infection. We estimated the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. We also used a different approach, based on results of classical genetics, to calculate HCV's mutation rate and obtained consistent results (2.8–3.2×10⁻⁵).

Evolutionary genomics of host-use in bifurcating demes of RNA virus phi-6

BACKGROUND

Viruses are exceedingly diverse in their evolved strategies to manipulate hosts for viral replication. However, despite these differences, most virus populations will occasionally experience two commonly-encountered challenges: growth in variable host environments, and growth under fluctuating population sizes. We used the segmented RNA bacteriophage ϕ6 as a model for studying the evolutionary genomics of virus adaptation in the face of host switches and parametrically varying population sizes. To do so, we created a bifurcating deme structure that reflected lineage splitting in natural populations, allowing us to test whether phylogenetic algorithms could accurately resolve this 'known phylogeny'. The resulting tree yielded 32 clones at the tips and internal nodes; these strains were fully sequenced and measured for phenotypic changes in selected traits (fitness on original and novel hosts).

RESULTS

We observed that RNA segment size was negatively correlated with the extent of molecular change in the imposed treatments; molecular substitutions tended to cluster on the Small and Medium RNA chromosomes of the virus, and not on the Large segment. Our study yielded a very large molecular and phenotypic dataset, fostering possible inferences on genotype-phenotype associations. Using further experimental evolution, we confirmed an inference on the unanticipated role of an allelic switch in a viral assembly protein, which governed viral performance across host environments.

CONCLUSIONS

Our study demonstrated that varying complexities can be simultaneously incorporated into experimental evolution, to examine the combined effects of population size, and adaptation in novel environments. The imposed bifurcating structure revealed that some methods for phylogenetic reconstruction failed to resolve the true phylogeny, owing to a paucity of molecular substitutions separating the RNA viruses that evolved in our study.

The three faces of riboviral spontaneous mutation: spectrum, mode of genome replication, and mutation rate

Riboviruses (RNA viruses without DNA replication intermediates) are the most abundant pathogens infecting animals and plants. Only a few riboviral infections can be controlled with antiviral drugs, mainly because of the rapid appearance of resistance mutations. Little reliable information is available concerning i) kinds and relative frequencies of mutations (the mutational spectrum), ii) mode of genome replication and mutation accumulation, and iii) rates of spontaneous mutation. To illuminate these issues, we developed a model in vivo system based on phage Qß infecting its natural host, Escherichia coli. The Qß RT gene encoding the Read-Through protein was used as a mutation reporter. To reduce uncertainties in mutation frequencies due to selection, the experimental Qß populations were established after a single cycle of infection and selection against RT⁻ mutants during phage growth was ameliorated by plasmid-based RT complementation in trans. The dynamics of Qß genome replication were confirmed to reflect the linear process of iterative copying (the stamping-machine mode). A total of 32 RT mutants were detected among 7,517 Qß isolates. Sequencing analysis of 45 RT mutations revealed a spectrum dominated by 39 transitions, plus 4 transversions and 2 indels. A clear template•primer mismatch bias was observed: A•C>C•A>U•G>G•U> transversion mismatches. The average mutation rate per base replication was ≈9.1×10⁻⁶ for base substitutions and ≈2.3×10⁻⁷ for indels. The estimated mutation rate per genome replication, μ_g, was ≈0.04 (or, per phage generation, ≈0.08), although secondary RT mutations arose during the growth of some RT mutants at a rate about 7-fold higher, signaling the possible impact of transitory bouts of hypermutation. These results are contrasted with those previously reported for other riboviruses to depict the current state of the art in riboviral mutagenesis.

Viral disease is a subject of major concern in public health. Diseases produced by riboviruses (RNA viruses sensu stricto) represent a special urgency, because these viruses display an exceptional capability to generate resistance mutations against antiviral drugs. Unfortunately, little is known about the rate and nature of spontaneous mutation in riboviruses. Thus, characterization of their mutation process may be helpful in the development of improved ways to counteract riboviral diseases. In this study, we investigated the mutation process in vivo of a model ribovirus, the bacteriophage Qß, focusing on three key aspects: i) the kinds and relative frequencies of mutations, ii) the mode of genome replication, and iii) the rate of spontaneous mutation. Our results, combined with other information about riboviral mutagenesis, depict a ribovirus mutation spectrum largely dominated by transitions, a predominantly linear mode of genome replication, and a mutation rate per genome replication on the order of 0.04 for bacteriophages and plant viruses but perhaps an order of magnitude higher for mammalian riboviruses.

Biomedical implications of viral mutation and evolution

Mutation rates vary hugely across viruses and strongly determine their evolution. In addition, viral mutation and evolution are biomedically relevant because they can determine pathogenesis, vaccine efficacy and antiviral resistance. We review experimental methods for estimating viral mutation rates and how these estimates vary across viral groups, paying special attention to the more general trends. Recent advances positing a direct association between viral mutation rates and virulence, or the use of high-fidelity variants as attenuated vaccines, are also discussed. Finally, we review the implications of viral mutation and evolution for the design of rational antiviral therapies and for efficient epidemiological surveillance.

Quasispecies spatial models for RNA viruses with different replication modes and infection strategies

Empirical observations and theoretical studies suggest that viruses may use different replication strategies to amplify their genomes, which impact the dynamics of mutation accumulation in viral populations and therefore, their fitness and virulence. Similarly, during natural infections, viruses replicate and infect cells that are rarely in suspension but spatially organized. Surprisingly, most quasispecies models of virus replication have ignored these two phenomena. In order to study these two viral characteristics, we have developed stochastic cellular automata models that simulate two different modes of replication (geometric vs stamping machine) for quasispecies replicating and spreading on a two-dimensional space. Furthermore, we explored these two replication models considering epistatic fitness landscapes (antagonistic vs synergistic) and different scenarios for cell-to-cell spread, one with free superinfection and another with superinfection inhibition. We found that the master sequences for populations replicating geometrically and with antagonistic fitness effects vanished at low critical mutation rates. By contrast, the highest critical mutation rate was observed for populations replicating geometrically but with a synergistic fitness landscape. Our simulations also showed that for stamping machine replication and antagonistic epistasis, a combination that appears to be common among plant viruses, populations further increased their robustness by inhibiting superinfection. We have also shown that the mode of replication strongly influenced the linkage between viral loci, which rapidly reached linkage equilibrium at increasing mutations for geometric replication. We also found that the strategy that minimized the time required to spread over the whole space was the stamping machine with antagonistic epistasis among mutations. Finally, our simulations revealed that the multiplicity of infection fluctuated but generically increased along time.

Dynamics of alternative modes of RNA replication for positive-sense RNA viruses

We propose and study nonlinear mathematical models describing the intracellular time dynamics of viral RNA accumulation for positive-sense single-stranded RNA viruses. Our models consider different replication modes ranging between two extremes represented by the geometric replication (GR) and the linear stamping machine replication (SMR). We first analyse a model that quantitatively reproduced experimental data for the accumulation dynamics of both polarities of turnip mosaic potyvirus RNAs. We identify a non-degenerate transcritical bifurcation governing the extinction of both strands depending on three key parameters: the mode of replication (α), the replication rate (r) and the degradation rate (δ) of viral strands. Our results indicate that the bifurcation associated with α generically takes place when the replication mode is closer to the SMR, thus suggesting that GR may provide viral strands with an increased robustness against degradation. This transcritical bifurcation, which is responsible for the switching from an active to an absorbing regime, suggests a smooth (i.e. second-order), absorbing-state phase transition. Finally, we also analyse a simplified model that only incorporates asymmetry in replication tied to differential replication modes.

Dynamics of a plant RNA virus intracellular accumulation: stamping machine vs. geometric replication

The tremendous evolutionary potential of RNA viruses allows them to thrive despite host defense mechanisms and endows them with properties such as emergence, host switching, and virulence. The frequency of mutant viruses after an infectious process results from the interplay between the error rate of the viral replicase, from purifying mechanisms acting after transcription on aberrant RNAs, and from the amplification dynamics of virus RNA positive (+) and negative (-) strands. Two extreme scenarios describing viral RNA amplification are the geometric growth, in which each RNA strand serves as template for the synthesis of complementary strands with the same efficiency, and the stamping machine, where a strand is reiteratively used as template to synthesize multiple copies of the complementary. The resulting mutation frequencies are completely different, being geometric growth largely more mutagenic than stamping machine. In this work we evaluate the contribution of geometric growth and stamping machine to the overall genome amplification of the plant (+)-strand RNA virus turnip mosaic virus (TuMV). By means of transfection experiments of Nicotiana benthamiana protoplasts with a TuMV cDNA infectious clone and by using strand-specific quantitative real-time PCR, we determined the amplification dynamics of viral (+) and (-) RNA during a single-cell infectious process. A mathematical model describing the amplification of each viral strand was fitted to the data. Analyses of the model parameters showed that TuMV (+) and (-) RNA amplification occurs through a mixed strategy with ∼93% of genomes produced via stamping machine and only ∼7% resulting from geometric growth.

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.

Viral mutation rates

Accurate estimates of virus mutation rates are important to understand the evolution of the viruses and to combat them. However, methods of estimation are varied and often complex. Here, we critically review over 40 original studies and establish criteria to facilitate comparative analyses. The mutation rates of 23 viruses are presented as substitutions per nucleotide per cell infection (s/n/c) and corrected for selection bias where necessary, using a new statistical method. The resulting rates range from 10(-8) to 10(-6) s/n/c for DNA viruses and from 10(-6) to 10(-4) s/n/c for RNA viruses. Similar to what has been shown previously for DNA viruses, there appears to be a negative correlation between mutation rate and genome size among RNA viruses, but this result requires further experimental testing. Contrary to some suggestions, the mutation rate of retroviruses is not lower than that of other RNA viruses. We also show that nucleotide substitutions are on average four times more common than insertions/deletions (indels). Finally, we provide estimates of the mutation rate per nucleotide per strand copying, which tends to be lower than that per cell infection because some viruses undergo several rounds of copying per cell, particularly double-stranded DNA viruses. A regularly updated virus mutation rate data set will be available at www.uv.es/rsanjuan/virmut.

The rate and spectrum of spontaneous mutations in a plant RNA virus

Knowing mutation rates and the molecular spectrum of spontaneous mutations is important to understanding how the genetic composition of viral populations evolves. Previous studies have shown that the rate of spontaneous mutations for RNA viruses widely varies between 0.01 and 2 mutations per genome and generation, with plant RNA viruses always occupying the lower side of this range. However, this peculiarity of plant RNA viruses is based on a very limited number of studies. Here we analyze the spontaneous mutational spectrum and the mutation rate of Tobacco etch potyvirus, a model system of positive sense RNA viruses. Our experimental setup minimizes the action of purifying selection on the mutational spectrum, thus giving a picture of what types of mutations are produced by the viral replicase. As expected for a neutral target, we found that transitions and nonsynonymous (including a few stop codons and small deletions) mutations were the most abundant type. This spectrum was notably different from the one previously described for another plant virus. We have estimated that the spontaneous mutation rate for this virus was in the range 10(-6)-10(-5) mutations per site and generation. Our estimates are in the same biological ballpark that previous values reported for plant RNA viruses. This finding gives further support to the idea that plant RNA viruses may have lower mutation rates than their animal counterparts.

Does mutational robustness inhibit extinction by lethal mutagenesis in viral populations?

Lethal mutagenesis is a promising new antiviral therapy that kills a virus by raising its mutation rate. One potential shortcoming of lethal mutagenesis is that viruses may resist the treatment by evolving genomes with increased robustness to mutations. Here, we investigate to what extent mutational robustness can inhibit extinction by lethal mutagenesis in viruses, using both simple toy models and more biophysically realistic models based on RNA secondary-structure folding. We show that although the evolution of greater robustness may be promoted by increasing the mutation rate of a viral population, such evolution is unlikely to greatly increase the mutation rate required for certain extinction. Using an analytic multi-type branching process model, we investigate whether the evolution of robustness can be relevant on the time scales on which extinction takes place. We find that the evolution of robustness matters only when initial viral population sizes are small and deleterious mutation rates are only slightly above the level at which extinction can occur. The stochastic calculations are in good agreement with simulations of self-replicating RNA sequences that have to fold into a specific secondary structure to reproduce. We conclude that the evolution of mutational robustness is in most cases unlikely to prevent the extinction of viruses by lethal mutagenesis.

The high mutation rate of RNA viruses, such as HIV, allows them to rapidly evolve resistance to host defenses and antiviral drugs. A new approach to treating these viruses—lethal mutagenesis—turns the mutation rate of these viruses against them. It uses mutagens to increase the viruses' mutation rates so much that the accumulation of harmful mutations drives viral populations to extinction. Is there any way that a virus could adapt to a drug that increases its mutation rate? One way is that the virus could evolve so that mutations tend to be less harmful. In previous experimental work, there have been reports that virus populations can differ in robustness. Yet, the evolution of mutational robustness did not seem to inhibit extinction by lethal mutagenesis. In this work, we model viral populations under lethal mutagenesis in order to see when viruses might escape extinction by evolving robustness to mutations. We find that viruses can benefit from robustness only at relatively low mutation rates because the extent to which robustness increases fitness is rapidly drowned out by the extent to which higher mutation rates decrease fitness. The implication is that the evolution of mutational robustness is not a fundamental impediment to lethal mutagenesis therapy.

Simple genomes, complex interactions: epistasis in RNA virus

Owed to their reduced size and low number of proteins encoded, RNA viruses and other subviral pathogens are often considered as being genetically too simple. However, this structural simplicity also creates the necessity for viral RNA sequences to encode for more than one protein and for proteins to carry out multiple functions, all together resulting in complex patterns of genetic interactions. In this work we will first review the experimental studies revealing that the architecture of viral genomes is dominated by antagonistic interactions among loci. Second, we will also review mathematical models and provide a description of computational tools for the study of RNA virus dynamics and evolution. As an application of these tools, we will finish this review article by analyzing a stochastic bit-string model of in silico virus replication. This model analyzes the interplay between epistasis and the mode of replication on determining the population load of deleterious mutations. The model suggests that, for a given mutation rate, the deleterious mutational load is always larger when epistasis is predominantly antagonistic than when synergism is the rule. However, the magnitude of this effect is larger if replication occurs geometrically than if it proceeds linearly.

The relationship between mutation frequency and replication strategy in positive-sense single-stranded RNA viruses

For positive-sense single-stranded RNA virus genomes, there is a trade-off between the mutually exclusive tasks of transcription, translation and encapsidation. The replication strategy that maximizes the intracellular growth rate of the virus requires iterative genome transcription from positive to negative, and back to positive sense. However, RNA viruses experience high mutation rates, and the proportion of genomes with lethal mutations increases with the number of replication cycles. Thus, intracellular mutant frequency will depend on the replication strategy. Introducing apparently realistic mutation rates into a model of viral replication demonstrates that strategies that maximize viral growth rate could result in an average of 26 mutations per genome by the time plausible numbers of positive strands have been generated, and that virus viability could be as low as 0.1 per cent. At high mutation rates or when a high proportion of mutations are deleterious, the optimal strategy shifts towards synthesizing more negative strands per positive strand, and in extremis towards a 'stamping-machine' replication mode where all the encapsidated genomes come from only two transcriptional steps. We conclude that if viral mutation rates are as high as current estimates suggest, either mutation frequency must be considerably higher than generally anticipated and the proportion of viable viruses produced extremely small, or replication strategies cannot be optimized to maximize viral growth rate. Mechanistic models linking mutation frequency to replication mechanisms coupled with data generated through new deep-sequencing technologies could play an important role in improving the estimates of viral mutation rate.

Replication mode and landscape topology differentially affect RNA virus mutational load and robustness

Regardless of genome polarity, intermediaries of complementary sense must be synthesized and used as templates for the production of new genomic strands. Depending on whether these new genomic strands become themselves templates for producing extra antigenomic ones, thus giving rise to geometric growth, or only the firstly synthesized antigenomic strands can be used to this end, thus following Luria's stamping machine model, the abundances and distributions of mutant genomes will be different. Here we propose mathematical and bit string models that allow distinguishing between stamping machine and geometric replication. We have observed that, regardless the topology of the fitness landscape, the critical mutation rate at which the master sequence disappears increases as the mechanism of replication switches from purely geometric to stamping machine. We also found that, for a wide range of mutation rates, large-effect mutations do not accumulate regardless the scheme of replication. However, mild mutations accumulate more in the geometric model. Furthermore, at high mutation rates, geometric growth leads to a population collapse for intermediate values of mutational effects at which the stamping machine still produces master genomes. We observed that the critical mutation rate was weakly dependent on the strength of antagonistic epistasis but strongly dependent on synergistic epistasis. In conclusion, we have shown that RNA viruses may increase their robustness against the accumulation of deleterious mutations by replicating as stamping machines and that the magnitude of this benefit depends on the topology of the fitness landscape assumed.

Point mutation rate of bacteriophage PhiX174

The point mutation rate of phage PhiX174 was determined using the fluctuation test. After identifying the genetic changes associated with the selected phenotype, we obtained an estimate of 1.0x10(-6) substitutions per base per round of copying, which is consistent with Drake's rule (0.003 mutations per genome per round of copying in DNA-based microorganisms).

RNA conformational changes in the life cycles of RNA viruses, viroids, and virus-associated RNAs

The rugged nature of the RNA structural free energy landscape allows cellular RNAs to respond to environmental conditions or fluctuating levels of effector molecules by undergoing dynamic conformational changes that switch on or off activities such as catalysis, transcription or translation. Infectious RNAs must also temporally control incompatible activities and rapidly complete their life cycle before being targeted by cellular defenses. Viral genomic RNAs must switch between translation and replication, and untranslated subviral RNAs must control other activities such as RNA editing or self-cleavage. Unlike well characterized riboswitches in cellular RNAs, the control of infectious RNA activities by altering the configuration of functional RNA domains has only recently been recognized. In this review, we will present some of these molecular rearrangements found in RNA viruses, viroids and virus-associated RNAs, relating how these dynamic regions were discovered, the activities that might be regulated, and what factors or conditions might cause a switch between conformations.

Robustness promotes evolvability of thermotolerance in an RNA virus

BACKGROUND

The ability for an evolving population to adapt to a novel environment is achieved through a balance of robustness and evolvability. Robustness is the invariance of phenotype in the face of perturbation and evolvability is the capacity to adapt in response to selection. Genetic robustness has been posited, depending on the underlying mechanism, to either decrease the efficacy of selection, or increase the possibility of future adaptation. However, the true effect of genetic robustness on evolvability in biological systems remains uncertain.

RESULTS

Here we demonstrate that genetic robustness increases evolvability of thermotolerance in laboratory populations of the RNA virus phi6. We observed that populations founded by robust clones evolved greater resistance to heat shock, relative to populations founded by brittle (less-robust) clones. Thus, we provide empirical evidence for the idea that robustness can promote evolvability in this environment, and further suggest that evolvability can arise indirectly via selection for robustness, rather than through direct selective action.

CONCLUSION

Our data imply that greater tolerance of mutational change is associated with virus adaptability in a new niche, a finding generally relevant to evolutionary biology, and informative for elucidating how viruses might evolve to emerge in new habitats and/or overcome novel therapies.

Structural domains within the 3' untranslated region of Turnip crinkle virus

The genomes of positive-strand RNA viruses undergo conformational shifts that complicate efforts to equate structures with function. We have initiated a detailed analysis of secondary and tertiary elements within the 3' end of Turnip crinkle virus (TCV) that are required for viral accumulation in vivo. MPGAfold, a massively parallel genetic algorithm, suggested the presence of five hairpins (H4a, H4b, and previously identified hairpins H4, H5, and Pr) and one H-type pseudoknot (Psi(3)) within the 3'-terminal 194 nucleotides (nt). In vivo compensatory mutagenesis analyses confirmed the existence of H4a, H4b, Psi(3) and a second pseudoknot (Psi(2)) previously identified in a TCV satellite RNA. In-line structure probing of the 194-nt fragment supported the coexistence of H4, H4a, H4b, Psi(3) and a pseudoknot that connects H5 and the 3' end (Psi(1)). Stepwise replacements of TCV elements with the comparable elements from Cardamine chlorotic fleck virus indicated that the complete 142-nt 3' end, and subsets containing Psi(3), H4a, and H4b or Psi(3), H4a, H4b, H5, and Psi(2), form functional domains for virus accumulation in vivo. A new 3-D molecular modeling protocol (RNA2D3D) predicted that H4a, H4b, H5, Psi(3), and Psi(2) are capable of simultaneous existence and bears some resemblance to a tRNA. The related Japanese iris necrotic ring virus does not have comparable domains. These results provide a framework for determining how interconnected elements participate in processes that require 3' untranslated region sequences such as translation and replication.

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.

Host density impacts relative fitness of bacteriophage Phi6 genotypes in structured habitats

Spatially structured environments may impact evolution by restricting population sizes, limiting opportunities for genetic mixis, or weakening selection against deleterious genotypes. When habitat structure impedes dispersal, low-productivity (less virulent) infectious parasites may benefit from their prudent exploitation of local hosts. Here we explored the combined ability for habitat structure and host density to dictate the relative reproductive success of differentially productive parasites. To do so, we allowed two RNA bacteriophage Phi6 genotypes to compete in structured and unstructured (semi-solid versus liquid) habitats while manipulating the density of Pseudomonas hosts. In the unstructured habitats, the more-productive phage strain experienced a relatively constant fitness advantage regardless of starting host density. By contrast, in structured habitats, restricted phage dispersal may have magnified the importance of local productivity, thus allowing the relative fitness of the less-productive virus to improve as host density increased. Further data suggested that latent period (duration of cellular infection) and especially burst size (viral progeny produced per cell) were the phage "life-history" traits most responsible for our results. We discuss the relevance of our findings for selection occurring in natural phage populations and for the general evolutionary epidemiology of infectious parasites.