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

The balance between fitness advantages and costs drives adaptation of bacteriophage Qβ to changes in host density at different temperatures

Introduction

Host density is one of the main factors affecting the infective capacity of viruses. When host density is low, it is more difficult for the virus to find a susceptible cell, which increases its probability of being damaged by the physicochemical agents of the environment. Nevertheless, viruses can adapt to variations in host density through different strategies that depend on the particular characteristics of the life cycle of each virus. In a previous work, using the bacteriophage Qβ as an experimental model, we found that when bacterial density was lower than optimal the virus increased its capacity to penetrate into the bacteria through a mutation in the minor capsid protein (A1) that is not described to interact with the cell receptor.

Results

Here we show that the adaptive pathway followed by Qβ in the face of similar variations in host density depends on environmental temperature. When the value for this parameter is lower than optimal (30°C), the mutation selected is the same as at the optimal temperature (37°C). However, when temperature increases to 43°C, the mutation selected is located in a different protein (A2), which is involved both in the interaction with the cell receptor and in the process of viral progeny release. The new mutation increases the entry of the phage into the bacteria at the three temperatures assayed. However, it also considerably increases the latent period at 30 and 37°C, which is probably the reason why it is not selected at these temperatures.

Conclusion

The conclusion is that the adaptive strategies followed by bacteriophage Qβ, and probably other viruses, in the face of variations in host density depend not only on their advantages at this selective pressure, but also on the fitness costs that particular mutations may present in function of the rest of environmental parameters that influence viral replication and stability.

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.

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.

The Single-Stranded RNA Bacteriophage Qβ Adapts Rapidly to High Temperatures: An Evolution Experiment

Single-stranded (ss)RNA viruses are thought to evolve rapidly due to an inherently high mutation rate. However, it remains unclear how ssRNA viruses adapt to novel environments and/or how many and what types of substitutions are needed to facilitate this evolution. In this study, we followed the adaptation of the ssRNA bacteriophage Qβ using thermally adapted Escherichia coli as a host, which can efficiently grow at temperatures between 37.2 and 45.3 °C. This made it possible to evaluate Qβ adaptation to the highest known temperature that supports growth, 45.3 °C. We found that Qβ was capable of replication at this temperature; within 114 days (~1260 generations), we detected more than 34 novel point mutations in the genome of the thermally adapted Qβ population, representing 0.8% of the total Qβ genome. In addition, we returned the 45.3 °C-adapted Qβ populations to 37.2 °C and passaged them for 8 days (~124 generations). We found that the reverse-adapted Qβ population showed little to no decrease in fitness. These results indicate that Qβ can evolve in response to increasing temperatures in a short period of time with the accumulation of point mutations.

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.

Existing Host Range Mutations Constrain Further Emergence of RNA Viruses

RNA viruses are capable of rapid host shifting, typically due to a point mutation that confers expanded host range. As additional point mutations are necessary for further expansions, epistasis among host range mutations can potentially affect the mutational neighborhood and frequency of niche expansion. We mapped the mutational neighborhood of host range expansion using three genotypes of the double-stranded RNA (dsRNA) bacteriophage φ6 (wild type and two isogenic host range mutants) on the novel host Pseudomonas syringae pv. atrofaciens. Both Sanger sequencing of 50 P. syringae pv. atrofaciens mutant clones for each genotype and population Illumina sequencing revealed the same high-frequency mutations allowing infection of P. syringae pv. atrofaciens. Wild-type φ6 had at least nine different ways of mutating to enter the novel host, eight of which are in p3 (host attachment protein gene), and 13/50 clones had unchanged p3 genes. However, the two isogenic mutants had dramatically restricted neighborhoods: only one or two mutations, all in p3. Deep sequencing revealed that wild-type clones without mutations in p3 likely had changes in p12 (morphogenic protein), a region that was not polymorphic for the two isogenic host range mutants. Sanger sequencing confirmed that 10/13 of the wild-type φ6 clones had nonsynonymous mutations in p12, and 2 others had point mutations in p9 and p5. None of these genes had previously been associated with host range expansion in φ6. We demonstrate, for the first time, epistatic constraint in an RNA virus due to host range mutations themselves, which has implications for models of serial host range expansion.IMPORTANCE RNA viruses mutate rapidly and frequently expand their host ranges to infect novel hosts, leading to serial host shifts. Using an RNA bacteriophage model system (Pseudomonas phage φ6), we studied the impact of preexisting host range mutations on another host range expansion. Results from both clonal Sanger and Illumina sequencing show that extant host range mutations dramatically narrow the neighborhood of potential host range mutations compared to that of wild-type φ6. This research suggests that serial host-shifting viruses may follow a small number of molecular paths to enter additional novel hosts. We also identified new genes involved in φ6 host range expansion, expanding our knowledge of this important model system in experimental evolution.

Rare haplotype load as marker for lethal mutagenesis

RNA viruses replicate with a template-copying fidelity, which lies close to an extinction threshold. Increases of mutation rate by nucleotide analogues can drive viruses towards extinction. This transition is the basis of an antiviral strategy termed lethal mutagenesis. We have introduced a new diversity index, the rare haplotype load (RHL), to describe NS5B (polymerase) mutant spectra of hepatitis C virus (HCV) populations passaged in absence or presence of the mutagenic agents favipiravir or ribavirin. The increase in RHL is more prominent in mutant spectra whose expansions were due to nucleotide analogues than to multiple passages in absence of mutagens. Statistical tests for paired mutagenized versus non-mutagenized samples with 14 diversity indices show that RHL provides consistently the highest standardized effect of mutagenic treatment difference for ribavirin and favipiravir. The results indicate that the enrichment of viral quasispecies in very low frequency minority genomes can serve as a robust marker for lethal mutagenesis. The diagnostic value of RHL from deep sequencing data is relevant to experimental studies on enhanced mutagenesis of viruses, and to pharmacological evaluations of inhibitors suspected to have a mutagenic activity.

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.