Diminishing returns of population size in the rate of RNA virus adaptation

Whenever an asexual viral population evolves by adapting to new environmental conditions, beneficial mutations, the ultimate cause of adaptation, are randomly produced and then fixed in the population. The larger the population size and the higher the mutation rate, the more beneficial mutations can be produced per unit time. With the usually high mutation rate of RNA viruses and in a large enough population, several beneficial mutations could arise at the same time but in different genetic backgrounds, and if the virus is asexual, they will never be brought together through recombination. Thus, the best of these genotypes must outcompete each other on their way to fixation. This competition among beneficial mutations has the effect of slowing the overall rate of adaptation. This phenomenon is known as clonal interference. Clonal interference predicts a speed limit for adaptation as the population size increases. In the present report, by varying the size of evolving vesicular stomatitis virus populations, we found evidence clearly demonstrating this speed limit and thus indicating that clonal interference might be an important factor modulating the rate of adaptation to an in vitro cell system. Several evolutionary and epidemiological implications of the clonal interference model applied to RNA viruses are discussed.     

Rate of Adaptation in Sexuals and Asexuals: A Solvable Model of the Fisher–Muller Effect

The adaptation of large asexual populations is hampered by the competition between independently arising beneficial mutations in different individuals, which is known as clonal interference. In classic work, Fisher and Muller proposed that recombination provides an evolutionary advantage in large populations by alleviating this competition. Based on recent progress in quantifying the speed of adaptation in asexual populations undergoing clonal interference, we present a detailed analysis of the Fisher–Muller mechanism for a model genome consisting of two loci with an infinite number of beneficial alleles each and multiplicative (nonepistatic) fitness effects. We solve the deterministic, infinite population dynamics exactly and show that, for a particular, natural mutation scheme, the speed of adaptation in sexuals is twice as large as in asexuals. This result is argued to hold for any nonzero value of the rate of recombination. Guided by the infinite population result and by previous work on asexual adaptation, we postulate an expression for the speed of adaptation in finite sexual populations that agrees with numerical simulations over a wide range of population sizes and recombination rates. The ratio of the sexual to asexual adaptation speed is a function of population size that increases in the clonal interference regime and approaches 2 for extremely large populations. The simulations also show that the imbalance between the numbers of accumulated mutations at the two loci is strongly suppressed even by a small amount of recombination. The generalization of the model to an arbitrary number L of loci is briefly discussed. If each offspring samples the alleles at each locus from the gene pool of the whole population rather than from two parents, the ratio of the sexual to asexual adaptation speed is approximately equal to L in large populations. A possible realization of this scenario is the reassortment of genetic material in RNA viruses with L genomic segments.     

Contrasting dynamics of a mutator allele in asexual populations of differing size

Mutators have been shown to hitchhike in asexual populations when the anticipated beneficial mutation supply rate of the mutator subpopulation, NU(b) (for subpopulation of size N and beneficial mutation rate U(b)) exceeds that of the wild-type subpopulation. Here, we examine the effect of total population size on mutator dynamics in asexual experimental populations of Saccharomyces cerevisiae. Although mutators quickly hitchhike to fixation in smaller populations, mutator fixation requires more and more time as population size increases; this observed delay in mutator hitchhiking is consistent with the expected effect of clonal interference. Interestingly, despite their higher beneficial mutation supply rate, mutators are supplanted by the wild type in very large populations. We postulate that this striking reversal in mutator dynamics is caused by an interaction between clonal interference, the fitness cost of the mutator allele, and infrequent large-effect beneficial mutations in our experimental populations. Our work thus identifies a potential set of circumstances under which mutator hitchhiking can be inhibited in natural asexual populations, despite recent theoretical predictions that such populations should have a net tendency to evolve ever-higher genomic mutation rates.     

Interfering waves of adaptation promote spatial mixing

A fundamental problem of asexual adaptation is that beneficial substitutions are not efficiently accumulated in large populations: Beneficial mutations often go extinct because they compete with one another in going to fixation. It has been argued that such clonal interference may have led to the evolution of sex and recombination in well-mixed populations. Here, we study clonal interference, and mechanisms of its mitigation, in an evolutionary model of spatially structured populations with uniform selection pressure. Clonal interference is much more prevalent with spatial structure than without, due to the slow wave-like spread of beneficial mutations through space. We find that the adaptation speed of asexuals saturates when the linear habitat size exceeds a characteristic interference length, which becomes shorter with smaller migration and larger mutation rate. The limiting speed is proportional to μ(1/2) and μ(1/3) in linear and planar habitats, respectively, where the mutational supply μ is the product of mutation rate and local population density. This scaling and the existence of a speed limit should be amenable to experimental tests as they fall far below predicted adaptation speeds for well-mixed populations (that scale as the logarithm of population size). Finally, we show that not only recombination, but also long-range migration is a highly efficient mechanism of relaxing clonal competition in structured populations. Our conservative estimates of the interference length predict prevalent clonal interference in microbial colonies and biofilms, so clonal competition should be a strong driver of both genetic and spatial mixing in those contexts.     

The two-mutant problem: clonal interference in evolutionary graph theory

In large asexual populations, clonal interference, whereby different beneficial mutations compete to fix in the population simultaneously, may be the norm. Results extrapolated from the spread of individual mutations in homogeneous backgrounds are found to be misleading in such situations: clonal interference severely inhibits the spread of beneficial mutations. In contrast with results gained in systems with just one mutation striving for fixation at any one time, the spatial structure of the population is found to be an important factor in determining the fixation probability when there are two beneficial mutations.     

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.     

Fitness effects of fixed beneficial mutations in microbial populations

Beneficial mutations are intuitively relevant to understanding adaptation, yet not all beneficial mutations are of consequence to the long-term evolutionary outcome of adaptation. Many beneficial mutations-mostly those of small effect-are lost due either to (1) genetic drift or to (2) competition among clones carrying different beneficial mutations, a phenomenon called the "Hill-Robertson effect" for sexual populations and "clonal interference" for asexual populations. Competition among clones becomes more prevalent with increasing genetic linkage and increasing population size, and it is thus generally characteristic of microbial populations. Together, these two phenomena suggest that only those beneficial mutations of large fitness effect should achieve fixation, despite the fact that most beneficial mutations produced are predicted to have very small fitness effects. Here, we confirm this prediction-both empirically and theoretically-by showing that fitness effects of fixed beneficial mutations follow a distribution whose mode is positive.     

THE ADAPTATION RATE OF ASEXUALS: DELETERIOUS MUTATIONS,CLONAL INTERFERENCE AND POPULATION BOTTLENECKS

The rate at which a population adapts to its environment is a cornerstone of evolutionary theory, and recent experimental advances in microbial populations have renewed interest in predicting and testing this rate. Efforts to understand the adaptation rate theoretically are complicated by high mutation rates, to both beneficial and deleterious mutations, and by the fact that beneficial mutations compete with each other in asexual populations (clonal interference). Testable predictions must also include the effects of population bottlenecks, repeated reductions in population size imposed by the experimental protocol. In this contribution, we integrate previous work that addresses each of these issues, developing an overall prediction for the adaptation rate that includes: beneficial mutations with probabilistically distributed effects, deleterious mutations of arbitrary effect, population bottlenecks, and clonal interference.     

THE EFFECTS OF POPULATION BOTTLENECKS ON CLONAL INTERFERENCE, AND THE ADAPTATION EFFECTIVE POPULATION SIZE

Clonal interference refers to the competition that arises in asexual populations when multiple beneficial mutations segregate simultaneously. A large body of theoretical and experimental work now addresses this issue. Although much of the experimental work is performed in populations that grow exponentially between periodic population bottlenecks, the theoretical work to date has addressed only populations of a constant size. We derive an analytical approximation for the rate of adaptation in the presence of both clonal interference and bottlenecks, and compare this prediction to the results of an individual-based simulation, showing excellent agreement in the parameter regime in which clonal interference prevails. We also derive an appropriate definition for the effective population size for adaptive evolution experiments in the presence of population bottlenecks. This "adaptation effective population size" allows for a good approximation of the expected rate of adaptation, either in the strong-selection weak-mutation regime, or when clonal interference comes into play. In the multiple mutation regime, when the product of the population size and mutation rate is extremely large, these results no longer hold.     

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.     

Harnessing recombination to speed adaptive evolution in Escherichia coli

Evolutionary engineering typically involves asexual propagation of a strain to improve a desired phenotype. However, asexual populations suffer from extensive clonal interference, a phenomenon where distinct lineages of beneficial clones compete and are often lost from the population given sufficient time. Improved adaptive mutants can likely be generated by genetic exchange between lineages, thereby reducing clonal interference. We present a system that allows continuous in situ recombination by using an Esherichia coli F-based conjugation system lacking surface exclusion. Evolution experiments revealed that Hfr-mediated recombination significantly speeds adaptation in certain circumstances. These results show that our system is stable, effective, and suitable for use in evolutionary engineering applications.     

The influence of deleterious mutations on adaptation in asexual populations

We study the dynamics of adaptation in asexual populations that undergo both beneficial and deleterious mutations. In particular, how the deleterious mutations affect the fixation of beneficial mutations was investigated. Using extensive Monte Carlo simulations, we find that in the "strong-selection weak mutation (SSWM)" regime or in the "clonal interference (CI)" regime, deleterious mutations rarely influence the distribution of "selection coefficients of the fixed mutations (SCFM)"; while in the "multiple mutations" regime, the accumulation of deleterious mutations would lead to a decrease in fitness significantly. We conclude that the effects of deleterious mutations on adaptation depend largely on the supply of beneficial mutations. And interestingly, the lowest adaptation rate occurs for a moderate value of selection coefficient of deleterious mutations.     

The speed of evolution and maintenance of variation in asexual populations

BACKGROUND: The rate at which beneficial mutations accumulate determines how fast asexual populations evolve, but this is only partially understood. Some recent clonal-interference models suggest that evolution in large asexual populations is limited because smaller beneficial mutations are outcompeted by larger beneficial mutations that occur in different lineages within the same population. This analysis assumes that the important mutations fix one at a time; it ignores multiple beneficial mutations that occur in the lineage of an earlier beneficial mutation, before the first mutation in the series can fix. We focus on the effects of such multiple mutations. RESULTS: Our analysis predicts that the variation in fitness maintained by a continuously evolving population increases as the logarithm of the population size and logarithm of the mutation rate and thus yields a similar logarithmic increase in the speed of evolution. To test these predictions, we evolved asexual budding yeast in glucose-limited media at a range of population sizes and mutation rates. CONCLUSIONS: We find that their evolution is dominated by the accumulation of multiple mutations of moderate effect. Our results agree with our theoretical predictions and are inconsistent with the one-by-one fixation of mutants assumed by recent clonal-interference analysis.     

Clonal interference, multiple mutations and adaptation in large asexual populations

Two important problems affect the ability of asexual populations to accumulate beneficial mutations and hence to adapt. First, clonal interference causes some beneficial mutations to be outcompeted by more-fit mutations that occur in the same genetic background. Second, multiple mutations occur in some individuals, so even mutations of large effect can be outcompeted unless they occur in a good genetic background that contains other beneficial mutations. In this article, we use a Monte Carlo simulation to study how these two factors influence the adaptation of asexual populations. We find that the results depend qualitatively on the shape of the distribution of the fitness effects of possible beneficial mutations. When this distribution falls off slower than exponentially, clonal interference alone reasonably describes which mutations dominate the adaptation, although it gives a misleading picture of the evolutionary dynamics. When the distribution falls off faster than exponentially, an analysis based on multiple mutations is more appropriate. Using our simulations, we are able to explore the limits of validity of both of these approaches, and we explore the complex dynamics in the regimes where neither one is fully applicable.     

Competition between high- and higher-mutating strains of Escherichia coli

Experimental studies have shown that a mutator allele can readily hitchhike to fixation with beneficial mutations in an asexual population having a low, wild-type mutation rate. Here, we show that a genotype bearing two mutator alleles can supplant a population already fixed for one mutator allele. Our results provide experimental support for recent theory predicting that mutator alleles will tend to accumulate in asexual populations by hitchhiking with beneficial mutations, causing an ever-higher genomic mutation rate.     

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

A comparison of directed evolution approaches using the beta-glucuronidase model system

Protein engineers can alter the properties of enzymes by directing their evolution in vitro. Many methods to generate molecular diversity and to identify improved clones have been developed, but experimental evolution remains as much an art as a science. We previously used DNA shuffling (sexual recombination) and a histochemical screen to direct the evolution of Escherichia coli beta-glucuronidase (GUS) variants with improved beta-galactosidase (BGAL) activity. Here, we employ the same model evolutionary system to test the efficiencies of several other techniques: recursive random mutagenesis (asexual), combinatorial cassette mutagenesis (high-frequency recombination) and a versatile high-throughput microplate screen. GUS variants with altered specificity evolved in each trial, but different combinations of mutagenesis and screening techniques effected the fixation of different beneficial mutations. The new microplate screen identified a broader set of mutations than the previously employed X-gal colony screen. Recursive random mutagenesis produced essentially asexual populations, within which beneficial mutations drove each other into extinction (clonal interference); DNA shuffling and combinatorial cassette mutagenesis led instead to the accumulation of beneficial mutations within a single allele. These results explain why recombinational approaches generally increase the efficiency of laboratory evolution.     

The evolution of mutation rates: separating causes from consequences

Natural selection can adjust the rate of mutation in a population by acting on allelic variation affecting processes of DNA replication and repair. Because mutation is the ultimate source of the genetic variation required for adaptation, it can be appealing to suppose that the genomic mutation rate is adjusted to a level that best promotes adaptation. Most mutations with phenotypic effects are harmful, however, and thus there is relentless selection within populations for lower genomic mutation rates. Selection on beneficial mutations can counter this effect by favoring alleles that raise the mutation rate, but the effect of beneficial mutations on the genomic mutation rate is extremely sensitive to recombination and is unlikely to be important in sexual populations. In contrast, high genomic mutation rates can evolve in asexual populations under the influence of beneficial mutations, but this phenomenon is probably of limited adaptive significance and represents, at best, a temporary reprieve from the continual selection pressure to reduce mutation. The physiological cost of reducing mutation below the low level observed in most populations may be the most important factor in setting the genomic mutation rate in sexual and asexual systems, regardless of the benefits of mutation in producing new adaptive variation. Maintenance of mutation rates higher than the minimum set by this "cost of fidelity" is likely only under special circumstances.     

A Ruby in the Rubbish: Beneficial Mutations, Deleterious Mutations and the Evolution of Sex

This study presents a mathematical model in which a single beneficial mutation arises in a very large population that is subject to frequent deleterious mutations. The results suggest that, if the population is sexual, then the deleterious mutations will have little effect on the ultimate fate of the beneficial mutation. However, if most offspring are produced asexually, then the probability that the beneficial mutation will be lost from the population may be greatly enhanced by the deleterious mutations. Thus, sexual populations may adapt much more quickly than populations where most reproduction is asexual. Some of the results were produced using computer simulation methods, and a technique was developed that allows treatment of arbitrarily large numbers of individuals in a reasonable amount of computer time. This technique may be of prove useful for the analysis of a wide variety of models, though there are some constraints on its applicability. For example, the technique requires that reproduction can be described by Poisson processes.     

Beneficial mutation selection balance and the effect of linkage on positive selection

When beneficial mutations are rare, they accumulate by a series of selective sweeps. But when they are common, many beneficial mutations will occur before any can fix, so there will be many different mutant lineages in the population concurrently. In an asexual population, these different mutant lineages interfere and not all can fix simultaneously. In addition, further beneficial mutations can accumulate in mutant lineages while these are still a minority of the population. In this article, we analyze the dynamics of such multiple mutations and the interplay between multiple mutations and interference between clones. These result in substantial variation in fitness accumulating within a single asexual population. The amount of variation is determined by a balance between selection, which destroys variation, and beneficial mutations, which create more. The behavior depends in a subtle way on the population parameters: the population size, the beneficial mutation rate, and the distribution of the fitness increments of the potential beneficial mutations. The mutation-selection balance leads to a continually evolving population with a steady-state fitness variation. This variation increases logarithmically with both population size and mutation rate and sets the rate at which the population accumulates beneficial mutations, which thus also grows only logarithmically with population size and mutation rate. These results imply that mutator phenotypes are less effective in larger asexual populations. They also have consequences for the advantages (or disadvantages) of sex via the Fisher-Muller effect; these are discussed briefly.