Temporal variation in selection accelerates mutational decay by Muller's ratchet

Asexual species accumulate deleterious mutations through an irreversible process known as Muller's ratchet. Attempts to quantify the rate of the ratchet have ignored the role of temporal environmental heterogeneity even though it is common in nature and has the potential to affect overall ratchet rate. Here we examine Muller's ratchet in the context of conditional neutrality (i.e., mutations that are deleterious in some environmental conditions but neutral in others) as well as more subtle changes in the strength (but not sign) of selection. We find that temporal variation increases the rate of the ratchet (mutation accumulation) and the rate of fitness decline over that of populations experiencing constant selection of equivalent average strength. Temporal autocorrelation magnifies the effects of temporal heterogeneity and can allow the ratchet to operate at large population sizes in which it would be halted under constant selection. Classic studies of Muller's ratchet show that the rate of fitness decline is maximized when selection is of a low but intermediate strength. This relationship changes quantitatively with all forms of temporal heterogeneity studied and changes qualitatively when there is temporal autocorrelation in selection. In particular, the rate of fitness decline can increase indefinitely with the strength of selection with some forms of temporal heterogeneity. Our finding that temporal autocorrelation in selection dramatically increases ratchet rate and rate of fitness decline may help to explain the paucity of asexual taxa.     

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.     

Dynamic mutation-selection balance as an evolutionary attractor

The vast majority of mutations are deleterious and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet can lead to a rapid degradation of population fitness. Evidently, the long-term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Hence any stable evolutionary state of a population in a static environment must involve a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that such a state can exist for any population size N and mutation rate U and calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low ε suffices to achieve stability, even in small populations in the face of high mutation rates and weak selection, maintaining a well-adapted population in spite of Muller's ratchet. This may explain the maintenance of mitochondria and other asexual genomes.     

Muller's ratchet and the pattern of variation at a neutral locus

The levels and patterns of variation at a neutral locus are analyzed in a haploid asexual population undergoing accumulation of deleterious mutations due to Muller's ratchet. We find that the movement of Muller's ratchet can be associated with a considerable reduction in genetic diversity below classical neutral expectation. The extent to which variability is reduced is a function of the deleterious mutation rate, the fitness effects of the mutations, and the population size. Approximate analytical expressions for the expected genetic diversity are compared with simulation results under two different models of deleterious mutations: a model where all deleterious mutations have equal effects and a model where there are two classes of deleterious mutations. We also find that Muller's ratchet can produce a considerable distortion in the neutral frequency spectrum toward an excess of rare variants.     

Effect of deleterious mutation-accumulation on the fitness of RNA bacteriophage MS2

Abstract.— RNA viruses show the highest mutation rate in nautre. It has been extensively demonstrated that, in the absence of purifying selection, RNA viruses accumulate deleterious mutations at a high rate. However, the parameters describing this accumulation are, in general, poorly understood. The present study reports evidences for fitness declines by the accumulation of deleterious mutations in the bacteriophage MS2. We estimated the rate of fitness decline to be as high as 16% per bottleneck transfer. In addition, our results agree with an additive model of fitness effects.     

Multiple infection dynamics has pronounced effects on the fitness of RNA viruses

Several factors play a role during the replication and transmission of RNA viruses. First, as a consequence of their enormous mutation rate, complex mixtures of genomes are generated immediately after infection of a new host. Secondly, differences in growth and competition rates drive the selection of certain genetic variants within an infected host. Thirdly, but not less important, a random sampling occurs at the moment of viral infectious passage from an infected to a healthy host. In addition, the availability of hosts also influences the fate of a given viral genotype. When new hosts are scarce, different viral genotypes might infect the same host, adding an extra complexity to the competition among genetic variants. We have employed a two‐fold approach to analyse the role played by each of these factors in the evolution of RNA viruses. First, we have derived a model that takes into account all the preceding factors. This model employs the classic Lotka‐Volterra competition equations but it also incorporates the effect of mutation during RNA replication, the effect of the stochastic sampling at the moment of infectious passage among hosts and, the effect of the type of infection (single, coinfection or superinfection). Secondly, the predictions of the model have been tested in an in vitro evolution experiment. Both theoretical and experimental results show that in infection passages with coinfection viral fitness increased more than in single infections. In contrast, infection passages with superinfection did not differ from the single infection. The coinfection frequency also affected the outcome: the larger the proportion of viruses coinfecting a host, the larger increase in fitness observed.     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

The evolution of low mutation rates in experimental mutator populations of Saccharomyces cerevisiae

Mutation is the source of both beneficial adaptive variation and deleterious genetic load, fueling the opposing selective forces than shape mutation rate evolution. This dichotomy is well illustrated by the evolution of the mutator phenotype, a genome-wide 10- to 100-fold increase in mutation rate. This phenotype has often been observed in clonally expanding populations exposed to novel or frequently changing conditions. Although studies of both experimental and natural populations have shed light on the evolutionary forces that lead to the spread of the mutator allele through a population, significant gaps in our understanding of mutator evolution remain. Here we use an experimental evolution approach to investigate the conditions required for the evolution of a reduction in mutation rate and the mechanisms by which populations tolerate the accumulation of deleterious mutations. We find that after ～6,700 generations, four out of eight experimental mutator lines had evolved a decreased mutation rate. We provide evidence that the accumulation of deleterious mutations leads to selection for reduced mutation rate clones in populations of mutators. Finally, we test the long-term consequences of the mutator phenotype, finding that mutator lines follow different evolutionary trajectories, some of which lead to drug resistance.     

Mutation accumulation and the formation of range limits

The dynamics of range formation are important for understanding and predicting species distributions. Here, we focus on a process that has thus far been overlooked in the context of range formation; the accumulation of mutation load. We find that mutation accumulation severely reduces the extent of a range across an environmental gradient, especially when dispersal is limited, growth rate is low and mutations are of intermediate deleterious effect. Our results illustrate the important role deleterious mutations can play in range formation. We highlight this as a necessary focus for further work, noting particularly the potentially conflicting effects dispersal may have in reducing mutation load and simultaneously increasing migration load in marginal populations.     

Epistasis and its relationship to canalization in the RNA virus phi 6

Although deleterious mutations are believed to play a critical role in evolution, assessing their realized effect has been difficult. A key parameter governing the effect of deleterious mutations is the nature of epistasis, the interaction between the mutations. RNA viruses should provide one of the best systems for investigating the nature of epistasis because the high mutation rate allows a thorough investigation of mutational effects and interactions. Nonetheless, previous investigations of RNA viruses by S. Crotty and co-workers and by S. F. Elena have been unable to detect a significant effect of epistasis. Here we provide evidence that positive epistasis is characteristic of deleterious mutations in the RNA bacteriophage phi 6. We estimated the effects of deleterious mutations by performing mutation-accumulation experiments on five viral genotypes of decreasing fitness. We inferred positive epistasis because viral genotypes with low fitness were found to be less sensitive to deleterious mutations. We further examined environmental sensitivity in these genotypes and found that low-fitness genotypes were also less sensitive to environmental perturbations. Our results suggest that even random mutations impact the degree of canalization, the buffering of a phenotype against genetic and environmental perturbations. In addition, our results suggest that genetic and environmental canalization have the same developmental basis and finally that an understanding of the nature of epistasis may first require an understanding of the nature of canalization.     

Applying the genetic theories of ageing to the cytoplasm: cytoplasmic genetic covariation for fitness and lifespan

Two genetic models exist to explain the evolution of ageing – mutation accumulation (MA) and antagonistic pleiotropy (AP). Under MA, a reduced intensity of selection with age results in accumulation of late‐acting deleterious mutations. Under AP, late‐acting deleterious mutations accumulate because they confer beneficial effects early in life. Recent studies suggest that the mitochondrial genome is a major player in ageing. It therefore seems plausible that the MA and AP models will be relevant to genomes within the cytoplasm. This possibility has not been considered previously. We explore whether patterns of covariation between fitness and ageing across 25 cytoplasmic lines, sampled from a population of Drosophila melanogaster, are consistent with the genetic associations predicted under MA or AP. We find negative covariation for fitness and the rate of ageing, and positive covariation for fitness and lifespan. Notably, the direction of these associations is opposite to that typically predicted under AP.     

Human Cytomegalovirus Protein pUL117 Targets the Mini-Chromosome Maintenance Complex and Suppresses Cellular DNA Synthesis

Modulation of host DNA synthesis is essential for many viruses to establish productive infections and contributes to viral diseases. Human cytomegalovirus (HCMV), a large DNA virus, blocks host DNA synthesis and deregulates cell cycle progression. We report that pUL117, a viral protein that we recently identified, is required for HCMV to block host DNA synthesis. Mutant viruses in which pUL117 was disrupted, either by frame-shift mutation or by a protein destabilization-based approach, failed to block host DNA synthesis at times after 24 hours post infection in human foreskin fibroblasts. Furthermore, pUL117-deficient virus stimulated quiescent fibroblasts to enter S-phase, demonstrating the intrinsic ability of HCMV to promote host DNA synthesis, which was suppressed by pUL117. We examined key proteins known to be involved in inhibition of host DNA synthesis in HCMV infection, and found that many were unlikely involved in the inhibitory activity of pUL117, including geminin, cyclin A, and viral protein IE2, based on their expression patterns. However, the ability of HCMV to delay the accumulation of the mini-chromosome maintenance (MCM) complex proteins, represented by MCM2 and MCM4, and prevent their loading onto chromatin, was compromised in the absence of pUL117. When expressed alone, pUL117 slowed cell proliferation, delayed DNA synthesis, and inhibited MCM accumulation. Knockdown of MCM proteins by siRNA restored the ability of pUL117-deficient virus to block cellular DNA synthesis. Thus, targeting MCM complex is one mechanism pUL117 employs to help block cellular DNA synthesis during HCMV infection. Our finding substantiates an emerging picture that deregulation of MCM is a conserved strategy for many viruses to prevent host DNA synthesis and helps to elucidate the complex strategy used by a large DNA virus to modulate cellular processes to promote infection and pathogenesis.     

Increased burst size in multiply infected cells can alter basic virus dynamics

ABSTRACT: BACKGROUND: The dynamics of viral infections have been studied extensively in a variety of settings, both experimentally and with mathematical models. The majori-ty of mathematical models assumes that only one virus can infect a given cell at a time. It is, however, clear that especially in the context of high viral load, cells can become infected with multiple copies of a virus, a process called coinfection. This has been best demonstrated experimentally for human immunodeficiency virus (HIV), although it is thought to be equally relevant for a number of other viral infections. In a previously explored mathematical model, the viral output from an infected cell does not depend on the number of viruses that reside in the cell, i.e. viral replication is limited by cellular rather than viral factors. In this case, basic virus dynamics properties are not altered by coinfection. Results: Here, we explore the alternative assumption that multiply infected cells are characterized by an increased burst size and find that this can fundamentally alter model predictions. Under this scenario, establishment of infection may not be solely determined by the basic reproductive ratio of the virus, but can depend on the initial virus load. Upon infection, the virus population need not follow straight exponential growth. Instead, the exponential rate of growth can increase over time as virus load becomes larger. Moreover, the model suggests that the ability of anti-viral drugs to suppress the virus population can depend on the virus load upon initiation of therapy. This is because more coinfected cells, which produce more virus, are present at higher virus loads. Hence, the degree of drug resistance is not only determined by the viral genotype, but also by the prevalence of coinfected cells. Conclusions: Our work shows how an increased burst size in multiply infected cells can alter basic infection dynamics. This forms the basis for future experimental testing of model assumptions and predictions that can distinguish between the different scenarios.     

Recombination in HIV and the evolution of drug resistance: for better or for worse?

The rapid evolution of drug resistance remains a major obstacle for HIV therapy. The capacity of the virus for recombination is widely believed to facilitate the evolution of drug resistance. Here, we challenge this intuitive view. We develop a population genetic model of HIV replication that incorporates the processes of mutation, cellular superinfection, and recombination. We show that cellular superinfection increases the abundance of low fitness viruses at the expense of the fittest strains due to the mixing of viral proteins during virion assembly. Moreover, we argue that whether recombination facilitates the evolution of drug resistance depends critically on how resistance mutations interact to determine viral fitness. Contrary to the commonly held belief, we find that, under the most plausible biological assumptions, recombination is expected to slow down the rate of evolution of multi-drug-resistant virus during therapy.     

The extinction time under mutational meltdown driven by high mutation rates

Mutational meltdown describes an eco‐evolutionary process in which the accumulation of deleterious mutations causes a fitness decline that eventually leads to the extinction of a population. Possible applications of this concept include medical treatment of RNA virus infections based on mutagenic drugs that increase the mutation rate of the pathogen. To determine the usefulness and expected success of such an antiviral treatment, estimates of the expected time to mutational meltdown are necessary. Here, we compute the extinction time of a population under high mutation rates, using both analytical approaches and stochastic simulations. Extinction is the result of three consecutive processes: (a) initial accumulation of deleterious mutations due to the increased mutation pressure; (b) consecutive loss of the fittest haplotype due to Muller''s ratchet; (c) rapid population decline toward extinction. We find accurate analytical results for the mean extinction time, which show that the deleterious mutation rate has the strongest effect on the extinction time. We confirm that intermediate‐sized deleterious selection coefficients minimize the extinction time. Finally, our simulations show that the variation in extinction time, given a set of parameters, is surprisingly small.     

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.     

Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus

Considerable uncertainty surrounds the evolutionary rates of and selection pressures acting on arthropod-borne RNA viruses (arboviruses). In particular, it is unclear why arboviruses such as dengue virus show substantial genetic variation within individual humans and mosquitoes yet low long-term rates of amino acid substitution. To address this question, I compared patterns of nonsynonymous variation in populations of dengue virus sampled at different levels of evolutionary divergence. Although nonsynonymous variation was abundant in viral populations within individual humans, there was a marked reduction in the frequency of nonsynonymous mutations in interhost comparisons. Moreover, intrahost genetic variation corresponded to a random pattern of mutation, and most of the sites that exhibited nonsynonymous variation within hosts were invariant at deeper phylogenetic levels. This loss of long-term nonsynonymous variation is the signature of extensive purifying selection such that more than 90% of all nonsynonymous mutations are deleterious. Consequently, although arboviruses are able to successfully adapt to diverse cell types, they are characterized by a high rate of deleterious mutation.     

RNA interference as a tool for exploring HIV-1 robustness

Short interfering RNAs (siRNAs) that target viral genes can efficiently inhibit human immunodeficiency virus type 1 (HIV-1) replication. Nevertheless, there is the potential for viral escape, particularly with a highly mutable target such as HIV-1. We present a novel strategy for anticipating and preventing viral escape using second-generation siRNAs. The evolutionary capacity of HIV-1 was tested by exerting strong selective pressure on a highly conserved sequence in the HIV-1 genome. We assayed the antiviral efficacy of five overlapping siRNAs directed against an essential region of the HIV-1 protease. Serial viral transfers in U87-CD4-CXCR4 cells were performed using four of the siRNAs. This procedure was repeated until virus breakthrough was detected. After several serial culture passages, resistant virus with a single point mutation in the targeted region was detected in the culture supernatants. The emergence of resistant virus was confirmed by molecular cloning and DNA sequencing of viral RNA. The most common escape route was the D30N mutation. Importantly, the addition of a second-generation siRNA that matched the D30N mutation restored viral inhibition and delayed development of escape variants. Passages performed with both siRNAs prevented the emergence of the D30N escape mutant and forced the virus to develop new escape routes. Thus, second-generation siRNAs can be used to block escape from RNA interference (RNAi) and to search for new RNAi escape routes. The protocol described here may be useful for exploring the sequence space available for HIV-1 evolution and for producing attenuated or deleterious viruses.     

Elevated levels of expression associated with regions of the Drosophila genome that lack crossing over

The recombinational environment influences patterns of molecular evolution through the effects of Hill-Robertson interference. Here, we examine genome-wide patterns of gene expression with respect to recombinational environment in Drosophila melanogaster. We find that regions of the genome lacking crossing over exhibit elevated levels of expression, and this is most pronounced for genes on the entirely non-crossing over fourth chromosome. We find no evidence for differences in the patterns of gene expression between regions of high, intermediate and low crossover frequencies. These results suggest that, in the absence of crossing over, selection to maintain control of expression may be compromised, perhaps due to the accumulation of deleterious mutations in regulatory regions. Alternatively, higher gene expression may be evolving to compensate for defective protein products or reduced translational efficiency.