Evolution at increased error rate leads to the coexistence of multiple adaptive pathways in an RNA virus

Evolutionary Adaptation of an RNA Bacteriophage to Repeated Freezing and Thawing Cycles

Bacteriophage fitness is determined by factors influencing both their replication within bacteria and their ability to maintain infectivity between infections. The latter becomes particularly crucial under adverse environmental conditions or when host density is low. In such scenarios, the damage experienced by viral particles could lead to the loss of infectivity, which might be mitigated if the virus undergoes evolutionary optimization through replication. In this study, we conducted an evolution experiment involving bacteriophage Qβ, wherein it underwent 30 serial transfers, each involving a cycle of freezing and thawing followed by replication of the surviving viruses. Our findings show that Qβ was capable of enhancing its resistance to this selective pressure through various adaptive pathways that did not impair the virus replicative capacity. Notably, these adaptations predominantly involved mutations located within genes encoding capsid proteins. The adapted populations exhibited higher resistance levels than individual viruses isolated from them, and the latter surpassed those observed in single mutants generated via site-directed mutagenesis. This suggests potential interactions among mutants and mutations. In conclusion, our study highlights the significant role of extracellular selective pressures in driving the evolution of phages, influencing both the genetic composition of their populations and their phenotypic properties.

Standing Genetic Diversity and Transmission Bottleneck Size Drive Adaptation in Bacteriophage Qβ

A critical issue to understanding how populations adapt to new selective pressures is the relative contribution of the initial standing genetic diversity versus that generated de novo. RNA viruses are an excellent model to study this question, as they form highly heterogeneous populations whose genetic diversity can be modulated by factors such as the number of generations, the size of population bottlenecks, or exposure to new environment conditions. In this work, we propagated at nonoptimal temperature (43 °C) two bacteriophage Qβ populations differing in their degree of heterogeneity. Deep sequencing analysis showed that, prior to the temperature change, the most heterogeneous population contained some low-frequency mutations that had previously been detected in the consensus sequences of other Qβ populations adapted to 43 °C. Evolved populations with origin in this ancestor reached similar growth rates, but the adaptive pathways depended on the frequency of these standing mutations and the transmission bottleneck size. In contrast, the growth rate achieved by populations with origin in the less heterogeneous ancestor did depend on the transmission bottleneck size. The conclusion is that viral diversification in a particular environment may lead to the emergence of mutants capable of accelerating adaptation when the environment changes.

Senecavirus A Enhances Its Adaptive Evolution via Synonymous Codon Bias Evolution

Synonymous codon bias in the viral genome affects protein translation and gene expression, suggesting that the synonymous codon mutant plays an essential role in influencing virulence and evolution. However, how the recessive mutant form contributes to virus evolvability remains elusive. In this paper, we characterize how the Senecavirus A (SVA), a picornavirus, utilizes synonymous codon mutations to influence its evolution, resulting in the adaptive evolution of the virus to adverse environments. The phylogenetic tree and Median-joining (MJ)-Network of these SVA lineages worldwide were constructed to reveal SVA three-stage genetic development clusters. Furthermore, we analyzed the codon bias of the SVA genome of selected strains and found that SVA could increase the GC content of the third base of some amino acid synonymous codons to enhance the viral RNA adaptive evolution. Our results highlight the impact of recessive mutation of virus codon bias on the evolution of the SVA and uncover a previously underappreciated evolutionary strategy for SVA. They also underline the importance of understanding the genetic evolution of SVA and how SVA adapts to the adverse effects of external stress.

Intra-Population Competition during Adaptation to Increased Temperature in an RNA Bacteriophage

Evolution of RNA bacteriophages of the family Leviviridae is governed by the high error rates of their RNA-dependent RNA polymerases. This fact, together with their large population sizes, leads to the generation of highly heterogeneous populations that adapt rapidly to most changes in the environment. Throughout adaptation, the different mutants that make up a viral population compete with each other in a non-trivial process in which their selective values change over time due to the generation of new mutations. In this work we have characterised the intra-population dynamics of a well-studied levivirus, Qβ, when it is propagated at a higher-than-optimal temperature. Our results show that adapting populations experienced rapid changes that involved the ascent of particular genotypes and the loss of some beneficial mutations of early generation. Artificially reconstructed populations, containing a fraction of the diversity present in actual populations, fixed mutations more rapidly, illustrating how population bottlenecks may guide the adaptive pathways. The conclusion is that, when the availability of beneficial mutations under a particular selective condition is elevated, the final outcome of adaptation depends more on the occasional occurrence of population bottlenecks and how mutations combine in genomes than on the selective value of particular mutations.

Newly emerged resistance-breaking variants of cucumber mosaic virus represent ongoing host-interactive evolution of an RNA virus

Understanding the evolutionary history of a virus and the mechanisms influencing the direction of its evolution is essential for the development of more durable strategies to control the virus in crop fields. While the deployment of host resistance in crops is the most efficient means to control various viruses, host resistance itself can act as strong selective pressure and thus play a critical role in the evolution of virus virulence. Cucumber mosaic virus (CMV), a plant RNA virus with high evolutionary capacity, has caused endemic disease in various crops worldwide, including pepper (Capsicum annuum L.), because of frequent emergence of resistance-breaking variants. In this study, we examined the molecular and evolutionary characteristics of recently emerged, resistance-breaking CMV variants infecting pepper. Our population genetics analysis revealed that the high divergence capacity of CMV RNA1 might have played an essential role in the host-interactive evolution of CMV and in shaping the CMV population structure in pepper. We also demonstrated that nonsynonymous mutations in RNA1 encoding the 1a protein enabled CMV to overcome the deployed resistance in pepper. Our findings suggest that resistance-driven selective pressures on RNA1 might have contributed in shaping the unique evolutionary pattern of CMV in pepper. Therefore, deployment of a single resistance gene may reduce resistance durability against CMV and more integrated approaches are warranted for successful control of CMV in pepper.

Evolutionary Dynamics in the RNA Bacteriophage Qβ Depends on the Pattern of Change in Selective Pressures

The rate of change in selective pressures is one of the main factors that determines the likelihood that populations can adapt to stress conditions. Generally, the reduction in the population size that accompanies abrupt environmental changes makes it difficult to generate and select adaptive mutations. However, in systems with high genetic diversity, as happens in RNA viruses, mutations with beneficial effects under new conditions can already be present in the population, facilitating adaptation. In this work, we have propagated an RNA bacteriophage (Qβ) at temperatures higher than the optimum, following different patterns of change. We have determined the fitness values and the consensus sequences of all lineages throughout the evolutionary process in order to establish correspondences between fitness variations and adaptive pathways. Our results show that populations subjected to a sudden temperature change gain fitness and fix mutations faster than those subjected to gradual changes, differing also in the particular selected mutations. The life-history of populations prior to the environmental change has great importance in the dynamics of adaptation. The conclusion is that in the bacteriophage Qβ, the standing genetic diversity together with the rate of temperature change determine both the rapidity of adaptation and the followed evolutionary pathways.

Evolutionary adaptation of an RNA bacteriophage to the simultaneous increase in the within-host and extracellular temperatures

Bacteriophages are the most numerous biological entities on Earth. They are on the basis of most ecosystems, regulating the diversity and abundance of bacterial populations and contributing to the nutrient and energy cycles. Bacteriophages have two well differentiated phases in their life cycle, one extracellular, in which they behave as inert particles, and other one inside their hosts, where they replicate to give rise to a progeny. In both phases they are exposed to environmental conditions that often act as selective pressures that limit both their survival in the environment and their ability to replicate, two fitness traits that frequently cannot be optimised simultaneously. In this study we have analysed the evolutionary ability of an RNA bacteriophage, the bacteriophage Qβ, when it is confronted with a temperature increase that affects both the extracellular and the intracellular media. Our results show that Qβ can optimise its survivability when exposed to short-term high temperature extracellular heat shocks, as well as its replicative ability at higher-than-optimal temperature. Mutations responsible for simultaneous adaptation were the same as those selected when adaptation to each condition proceeded separately, showing the absence of important trade-offs between survival and reproduction in this virus.

Differences in adaptive dynamics determine the success of virus variants that propagate together

Virus fitness is a complex parameter that results from the interaction of virus-specific characters (e.g. intracellular growth rate, adsorption rate, virion extracellular stability, and tolerance to mutations) with others that depend on the underlying fitness landscape and the internal structure of the whole population. Individual mutants usually have lower fitness values than the complex population from which they come from. When they are propagated and allowed to attain large population sizes for a sufficiently long time, they approach mutation-selection equilibrium with the concomitant fitness gains. The optimization process follows dynamics that vary among viruses, likely due to differences in any of the parameters that determine fitness values. As a consequence, when different mutants spread together, the number of generations experienced by each of them prior to co-propagation may determine its particular fate. In this work we attempt a clarification of the effect of different levels of population diversity in the outcome of competition dynamics. To this end, we analyze the behavior of two mutants of the RNA bacteriophage Qβ that co-propagate with the wild-type virus. When both competitor viruses are clonal, the mutants rapidly outcompete the wild type. However, the outcome in competitions performed with partially optimized virus populations depends on the distance of the competitors to their clonal origin. We also implement a theoretical population dynamics model that describes the evolution of a heterogeneous population of individuals, each characterized by a fitness value, subjected to subsequent cycles of replication and mutation. The experimental results are explained in the framework of our theoretical model under two non-excluding, likely complementary assumptions: (1) The relative advantage of both competitors changes as populations approach mutation-selection equilibrium, as a consequence of differences in their growth rates and (2) one of the competitors is more robust to mutations than the other. The main conclusion is that the nearness of an RNA virus population to mutation-selection equilibrium is a key factor determining the fate of particular mutants arising during replication.

Adaptive multiscapes: an up-to-date metaphor to visualize molecular adaptation

Background

Wright’s metaphor of the fitness landscape has shaped and conditioned our view of the adaptation of populations for almost a century. Since its inception, and including criticism raised by Wright himself, the concept has been surrounded by controversy. Among others, the debate stems from the intrinsic difficulty to capture important features of the space of genotypes, such as its high dimensionality or the existence of abundant ridges, in a visually appealing two-dimensional picture. Two additional currently widespread observations come to further constrain the applicability of the original metaphor: the very skewed distribution of phenotype sizes (which may actively prevent, due to entropic effects, the achievement of fitness maxima), and functional promiscuity (i.e. the existence of secondary functions which entail partial adaptation to environments never encountered before by the population).

Results

Here we revise some of the shortcomings of the fitness landscape metaphor and propose a new “scape” formed by interconnected layers, each layer containing the phenotypes viable in a given environment. Different phenotypes within a layer are accessible through mutations with selective value, while neutral mutations cause displacements of populations within a phenotype. A different environment is represented as a separated layer, where phenotypes may have new fitness values, other phenotypes may be viable, and the same genotype may yield a different phenotype, representing genotypic promiscuity. This scenario explicitly includes the many-to-many structure of the genotype-to-phenotype map. A number of empirical observations regarding the adaptation of populations in the light of adaptive multiscapes are reviewed.

Conclusions

Several shortcomings of Wright’s visualization of fitness landscapes can be overcome through adaptive multiscapes. Relevant aspects of population adaptation, such as neutral drift, functional promiscuity or environment-dependent fitness, as well as entropic trapping and the concomitant impossibility to reach fitness peaks are visualized at once. Adaptive multiscapes should aid in the qualitative understanding of the multiple pathways involved in evolutionary dynamics.

Reviewers

This article was reviewed by Eugene Koonin and Ricard Solé.

Impact of increased mutagenesis on adaptation to high temperature in bacteriophage Qβ

RNA viruses replicate with very high error rates, which makes them more sensitive to additional increases in this parameter. This fact has inspired an antiviral strategy named lethal mutagenesis, which is based on the artificial increase of the error rate above a threshold incompatible with virus infectivity. A relevant issue concerning lethal mutagenesis is whether incomplete treatments might enhance the adaptive possibilities of viruses. We have addressed this question by subjecting an RNA virus, the bacteriophage Qβ, to different transmission regimes in the presence or the absence of sublethal concentrations of the mutagenic nucleoside analogue 5-azacytidine (AZC). Populations obtained were subsequently exposed to a non-optimal temperature and analyzed to determine their consensus sequences. Our results show that previously mutagenized populations rapidly fixed a specific set of mutations upon propagation at the new temperature, suggesting that the expansion of the mutant spectrum caused by AZC has an influence on later evolutionary behavior.

Artifact-Free Quantification and Sequencing of Rare Recombinant Viruses by Using Drop-Based Microfluidics

Recombination is an important driver in the evolution of viruses and thus is key to understanding viral epidemics and improving strategies to prevent future outbreaks. Characterization of rare recombinant subpopulations remains technically challenging because of artifacts such as artificial recombinants, known as chimeras, and amplification bias. To overcome this, we have developed a high-throughput microfluidic technique with a second verification step in order to amplify and sequence single recombinant viruses with high fidelity in picoliter drops. We obtained the first artifact-free estimate of in vitro recombination rate between murine norovirus strains MNV-1 and WU20 co-infecting a cell (P(rec) = 3.3 × 10(-4) ± 2 × 10(-5) ) for a 1205 nt region. Our approach represents a time- and cost-effective improvement over current methods, and can be adapted for genomic studies requiring artifact- and bias-free selective amplification, such as microbial pathogens, or rare cancer cells.

Fitness Landscapes of Functional RNAs

The notion of fitness landscapes, a map between genotype and fitness, was proposed more than 80 years ago. For most of this time data was only available for a few alleles, and thus we had only a restricted view of the whole fitness landscape. Recently, advances in genetics and molecular biology allow a more detailed view of them. Here we review experimental and theoretical studies of fitness landscapes of functional RNAs, especially aptamers and ribozymes. We find that RNA structures can be divided into critical structures, connecting structures, neutral structures and forbidden structures. Such characterisation, coupled with theoretical sequence-to-structure predictions, allows us to construct the whole fitness landscape. Fitness landscapes then can be used to study evolution, and in our case the development of the RNA world.

Temporal dynamics of intrahost molecular evolution for a plant RNA virus

Populations of plant RNA viruses are highly polymorphic in infected plants, which may allow rapid within-host evolution. To understand tobacco etch potyvirus (TEV) evolution, longitudinal samples from experimentally evolved populations in the natural host tobacco and from the alternative host pepper were phenotypically characterized and genetically analyzed. Temporal and compartmental variabilities of TEV populations were quantified using high throughput Illumina sequencing and population genetic approaches. Of the two viral phenotypic traits measured, virulence increased in the novel host but decreased in the original one, and viral load decreased in both hosts, though to a lesser extent in the novel one. Dynamics of population genetic diversity were also markedly different among hosts. Population heterozygosity increased in the ancestral host, with a dominance of synonymous mutations fixed, whereas it did not change or even decreased in the new host, with an excess of nonsynonymous mutations. All together, these observations suggest that directional selection is the dominant evolutionary force in TEV populations evolving in a novel host whereas either diversifying selection or random genetic drift may play a fundamental role in the natural host. To better understand these evolutionary dynamics, we developed a computer simulation model that incorporates the effects of mutation, selection, and drift. Upon parameterization with empirical data from previous studies, model predictions matched the observed patterns, thus reinforcing our idea that the empirical patterns of mutation accumulation represent adaptive evolution.

Emerging viruses: why they are not jacks of all trades?

In order to limit the impact of the recent pandemics ignited by viral host jumps, it is necessary to better understand the ecological and evolutionary factors influencing the early steps of emergence and the interactions between them. Antagonistic pleiotropy, that is, the negative fitness effect in the primary host of mutations allowing the infection of and the multiplication in a new host, has long been thought to be the main limitation to the evolution of generalist viruses and thus to emergence. However, the accumulation of experimental examples contradicting the hypothesis of antagonistic pleiotropy has highlighted the importance of other factors such as the epistasis between mutations increasing the adaptation to a new host. Epistasis is pervasive in viruses, affects the shape of the adaptive landscape and consequently the accessibility of evolutionary pathways. Finally, recent studies have gone steps further in the complexity of viral fitness determinism and stressed the potential importance of the epistatic pleiotropy and of the impact of host living conditions.

Changes in protein domains outside the catalytic site of the bacteriophage Qβ replicase reduce the mutagenic effect of 5-azacytidine

The high genetic heterogeneity and great adaptability of RNA viruses are ultimately caused by the low replication fidelity of their polymerases. However, single amino acid substitutions that modify replication fidelity can evolve in response to mutagenic treatments with nucleoside analogues. Here, we investigated how two independent mutants of the bacteriophage Qβ replicase (Thr210Ala and Tyr410His) reduce sensitivity to the nucleoside analogue 5-azacytidine (AZC). Despite being located outside the catalytic site, both mutants reduced the mutation frequency in the presence of the drug. However, they did not modify the type of AZC-induced substitutions, which was mediated mainly by ambiguous base pairing of the analogue with purines. Furthermore, the Thr210Ala and Tyr410His substitutions had little or no effect on replication fidelity in untreated viruses. Also, both substitutions were costly in the absence of AZC or when the action of the drug was suppressed by adding an excess of natural pyrimidines (uridine or cytosine). Overall, the phenotypic properties of these two mutants were highly convergent, despite the mutations being located in different domains of the Qβ replicase. This suggests that treatment with a given nucleoside analogue tends to select for a unique functional response in the viral replicase.

IMPORTANCE

In the last years, artificial increase of the replication error rate has been proposed as an antiviral therapy. In this study, we investigated the mechanisms by which two substitutions in the Qβ replicase confer partial resistance to the mutagenic nucleoside analogue AZC. As opposed to previous work with animal viruses, where different mutations selected sequentially conferred nucleoside analogue resistance through different mechanisms, our results suggest that there are few or no alternative AZC resistance phenotypes in Qβ. Also, despite resistance mutations being highly costly in the absence of the drug, there was no sequential fixation of secondary mutations. Bacteriophage Qβ is the virus with the highest reported mutation rate, which should make it particularly sensitive to nucleoside analogue treatments, probably favoring resistance mutations even if they incur high costs. The results are also relevant for understanding the possible pathways by which fidelity of the replication machinery can be modified.

Experimental evolution of an emerging plant virus in host genotypes that differ in their susceptibility to infection

This study evaluates the extent to which genetic differences among host individuals from the same species condition the evolution of a plant RNA virus. We performed a threefold replicated evolution experiment in which Tobacco etch potyvirus isolate At17b (TEV-At17b), adapted to Arabidopsis thaliana ecotype Ler-0, was serially passaged in five genetically heterogeneous ecotypes of A. thaliana. After 15 passages we found that evolved viruses improved their fitness, showed higher infectivity and stronger virulence in their local host ecotypes. The genome of evolved lineages was sequenced and putative adaptive mutations identified. Host-driven convergent mutations have been identified. Evidences supported selection for increased translational efficiency. Next, we sought for the specificity of virus adaptation by infecting all five ecotypes with all 15 evolved virus populations. We found that some ecotypes were more permissive to infection than others, and that some evolved virus isolates were more specialist/generalist than others. The bipartite network linking ecotypes with evolved viruses was significantly nested but not modular, suggesting that hard-to-infect ecotypes were infected by generalist viruses whereas easy-to-infect ecotypes were infected by all viruses, as predicted by a gene-for-gene model of infection.

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.

Adaptation to fluctuating temperatures in an RNA virus is driven by the most stringent selective pressure

The frequency of change in the selective pressures is one of the main factors driving evolution. It is generally accepted that constant environments select specialist organisms whereas changing environments favour generalists. The particular outcome achieved in either case also depends on the relative strength of the selective pressures and on the fitness costs of mutations across environments. RNA viruses are characterized by their high genetic diversity, which provides fast adaptation to environmental changes and helps them evade most antiviral treatments. Therefore, the study of the adaptive possibilities of RNA viruses is highly relevant for both basic and applied research. In this study we have evolved an RNA virus, the bacteriophage Qβ, under three different temperatures that either were kept constant or alternated periodically. The populations obtained were analyzed at the phenotypic and the genotypic level to characterize the evolutionary process followed by the virus in each case and the amount of convergent genetic changes attained. Finally, we also investigated the influence of the pre-existent genetic diversity on adaptation to high temperature. The main conclusions that arise from our results are: i) under periodically changing temperature conditions, evolution of bacteriophage Qβ is driven by the most stringent selective pressure, ii) there is a high degree of evolutionary convergence between replicated populations and also among populations evolved at different temperatures, iii) there are mutations specific of a particular condition, and iv) adaptation to high temperatures in populations differing in their pre-existent genetic diversity takes place through the selection of a common set of mutations.

Evolution of tertiary structure of viral RNA dependent polymerases

Viral RNA dependent polymerases (vRdPs) are present in all RNA viruses; unfortunately, their sequence similarity is too low for phylogenetic studies. Nevertheless, vRdP protein structures are remarkably conserved. In this study, we used the structural similarity of vRdPs to reconstruct their evolutionary history. The major strength of this work is in unifying sequence and structural data into a single quantitative phylogenetic analysis, using powerful a Bayesian approach.

The resulting phylogram of vRdPs demonstrates that RNA-dependent DNA polymerases (RdDPs) of viruses within Retroviridae family cluster in a clearly separated group of vRdPs, while RNA-dependent RNA polymerases (RdRPs) of dsRNA and +ssRNA viruses are mixed together. This evidence supports the hypothesis that RdRPs replicating +ssRNA viruses evolved multiple times from RdRPs replicating +dsRNA viruses, and vice versa. Moreover, our phylogram may be presented as a scheme for RNA virus evolution. The results are in concordance with the actual concept of RNA virus evolution. Finally, the methods used in our work provide a new direction for studying ancient virus evolution.

Strong epistatic interactions within a single protein

A central question in protein molecular evolution is whether an amino acid that occurs at a given site makes an independent contribution to fitness or whether its contribution depends on other amino acid sites in the protein sequence, a phenomenon known as intragenic epistasis. In the absence of intragenic epistasis, natural selection acts on a protein mutation independent of its genetic background, and the experimentally determined fitness for a mutation should be the same across all genetic backgrounds. We tested this hypothesis by using site-directed mutagenesis to introduce a well-characterized deleterious single amino acid substitution in 56 different hepatitis C virus NS3 protease variants. The catalytic efficiency of the mutated proteases was determined and compared with the corresponding wild-type protein. Fitness effects ranged from lethality to significantly beneficial. Although primarily deleterious and lethal effects were observed (41 and 5 out of 56 tested variants, respectively), deleterious effects ranged from near neutral (-26.7% reduction of fitness) to near lethal (-98.5%). Our findings demonstrate that the introduced amino acid substitution has different fitness effects in different protein variants and provide independent support for the relevant role of intragenic epistasis in protein evolution.