Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Curr Biol. Author manuscript; available in PMC 2010 November 18.
Published in final edited form as:
PMCID: PMC2987722

The Speed of Evolution and Maintenance of Variation in Asexual Populations



The rate at which beneficial mutations accumulate determines how fast asexual populations evolve, but 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 e ects of such multiple mutations.


Our analysis predicts that the variation in fitness that a continuously evolving population maintains increases as the logarithm of the population size and logarithm of the mutation rate, yielding 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.


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.

How do the mutation rate, population size, and the magnitude of beneficial mutations determine the rate at which asexual populations evolve? This question is important for comparing different experimental populations with each other and with natural populations, as well as for understanding the effects of factors that could increase the rate of evolution such as sex [110] or mutator phenotypes [6, 1115]. Lenski and others have found that, above a certain size, laboratory asexual populations do not evolve much faster at large population sizes and mutation rates than at small ones [4, 1619]. Their work is consistent with recent “clonal interference” theories of asexual evolution, in which the speed of evolution is limited in large populations because most beneficial mutations are outcompeted by larger beneficial mutations that occur in other lineages in the same population [2026]. However, these experiments did not monitor the distributions of fitnesses within populations as they evolved, which are important in discriminating between different theories.

We have measured the speed of evolution and the distribution of fitnesses in evolving asexual populations of budding yeast at a range of population sizes and two mutation rates. Our results are inconsistent with recent clonal interference analyses. Rather, they agree with theoretical predictions that the speed of evolution in large populations is dominated by multiple mutations that occur in the same lineage before the first mutation in the series has fixed. These mutations piggyback on each other in the sense that the presence of one beneficial mutation in a lineage helps another mutation in that lineage to outcompete other beneficial mutations of the same size that occur in other lineages that have fewer or smaller beneficial mutations. The balance between clonal interference and multiple mutations sets the typical size of the mutations that accumulate on top of each other.


Theories for the evolution of asexual populations come in three flavors. The first applies to small populations. If the effective size of a population is N and the beneficial mutation rate is Ub, new beneficial mutations arise at a rate NUb. Most of these mutations are lost by genetic drift before they become abundant enough to reliably prosper. But for a mutation with a selective advantage s, there is a probability s that the mutant lineage will survive drift and grow common enough for its selective advantage to take over [27]. We call this establishment of the mutant population. If nothing else interferes, mutations that establish will eventually outcompete and eliminate the genotype that gave rise to them. The time it takes a mutant to get from establishment to being half of the population is about 1sln[Ns], while the time between the establishment of successive mutations is 1NUbs. Thus when NUb1ln[Ns], mutations fix much more rapidly than they are established, so evolution is mutation-limited. The rate of evolution ν — defined as the rate of change of the mean fitness (more precisely, the mean log fitness) of the population — is νNUb[left angle bracket]s2[right angle bracket], where [left angle bracket]s2[right angle bracket] is the mean square s. This is known as the one locus regime (sometimes called the strong selection weak mutation regime); to make clear its nature, we refer to it as the successional fixation regime.

When mutations establish faster than they can fix, different mutations occur and spread through the population concurrently, and can interfere with each other (the “Hill-Robertson effect” [28, 29]). This concurrent mutations situation is more complicated than the successional fixation regime, and various analyses have considered different aspects of the dynamics. One approach (Fig. 1a) focuses on competition between mutations that have different fitness effects [2026]. This has been called “clonal interference” by its creators; we refer to it as “one-by-one clonal interference,” since it assumes that mutations fix one at a time. This analysis considers a mutation B with fitness advantage sB that becomes established in a population in which a different lineage with mutation A (with fitness advantage sA) is already spreading. If sA > sB, the lineage that carries B will be eliminated, whereas if sA < sB, lineage B can overtake and eliminate A. This process “wastes” some beneficial mutations and thus slows down the speed of evolution. Because more mutations are lost in larger populations, this analysis predicts that ν increases slowly as NUb rises; the details depend on the distribution of the strengths of beneficial mutations, which we call ρ(s).

Figure 1
Schematics of one-by-one clonal interference and multiple-mutations analyses. (a) One-by-one clonal interference. A population (the original wild-type (WT)) gives rise to several different beneficial mutations, which must compete. Here mutation C will ...

Although one-by-one clonal interference is one important aspect of the large-NUb dynamics, there is another crucial feature that affects the evolution of such populations. Even if a more-fit mutation B occurs before an earlier but less-fit mutation A fixes, A is not always doomed: an individual with mutation A can get an additional mutation C. If the combined fitness of A and C exceeds that of B, mutation A (along with C) can fix after all. One-by-one clonal interference neglects these complications. It assumes that mutations only occur in the majority population (“wild-type”), and that the most-fit such mutant outcompetes all others, becomes the new wild-type, and the process repeats. This is the one-by-one assumption. Yet the creation of multiple mutants is not a small effect: in populations large enough that one-by-one clonal interference is important, double mutants will routinely appear. Thus considering one-by-one clonal interference alone is incomplete. In typical laboratory yeast and microbial populations, N is often large enough for even triple and quadruple beneficial mutations to regularly occur before the first mutation in the series fixes [30]. Recent simulation studies [31, 32], as well as some experiments [33], have also indirectly hinted at the importance of these multiple mutation effects.

When multiple mutations are common, a different picture of the evolutionary dynamics is necessary. We must consider the competition between mutants that arise in lineages that already have other beneficial mutations. The beneficial mutations that matter most are those that occur in individuals that already have many others. Less-fit individuals which get an additional mutation will usually still be less fit than the most-fit individuals in the population, and hence doomed to eventual extinction. This effect also “wastes” beneficial mutations, and thus also causes the speed of evolution to increase only slowly as the total beneficial mutation rate NUb rises.

One-by-one clonal interference analyses focus on the competition between beneficial mutations arising from their different fitness effects while ignoring the competition between mutations based on the fitness of the individual they occurred in. We take here the opposite approach, focusing on the effects of multiple mutations accumulating in the same lineage (Fig. 1b). Specifically, our “multiple-mutations” analysis considers the possibility that a second beneficial mutation arises in the lineage of an earlier mutation A well before that lineage dominates the population. This creates a new more-fit lineage A*. A further mutation can convert A* into A**, but mutations in individuals other than A* are wasted. This is true even though mutations in A* individuals are much rarer than those in other individuals (since A* individuals are rare), because the mutations that confer a particular benefit in A* are more strongly selected for than mutations that confer the same benefit but occur in less fit lineages. Mutations in A** then create a still more-fit lineage A***, and so on. This establishment of mutations on top of existing unfixed mutations increases the width of the population fitness distribution (i.e. it increases the variation in fitness). But this tendency is countered by the selection against the remainder of the population, which is less fit than these multiply mutant individuals; this selection increases the mean fitness of the population and reduces the variation in fitness. Eventually these competing forces balance, creating a steady state shape of the evolving fitness distribution, with the broadening produced by the continuing establishment of mutations balanced by the narrowing by selection. When there is a large supply of possible beneficial mutations, this distribution moves continuously towards higher fitness while maintaining the shape set by the beneficial mutation-selection balance.

As NUb increases, the fitness distribution gets broader both because multiple mutations happen faster and because larger less-fit subpopulations take longer to eliminate. Concomitantly, the evolution gets faster because the speed of evolution is roughly equal to the fitness variance of the population. Our key result is that the speed of evolution and the variance in fitness both increase logarithmically with N and logarithmically with Ub (Box 1) — but not with the combination NUb. Our analysis is described in much more detail elsewhere [30].

Box 1

In the simple model of multiple mutations of the same effect s, there are two factors that determine the speed of evolution. The first is the dynamics of already established populations, which is dominated by selection. We define the lead of the distribution, qs, as the difference between the fitness of the most-fit individual and the mean fitness of the population; the fittest individuals have q more beneficial mutations than the mean. Once it is established, the fittest population grows exponentially, first at rate qs but more slowly as selection increases the mean fitness. Growing from its establishment upon reaching about 1qs individuals (which size it reaches rapidly if it does so at all [9]) until it reaches a large fraction of N will thus take time ln(Nqs)(qs2) (since qs2 is the average growth rate of the mutant relative to the mean growth rate during the period between establishment and fixation), and in this time the mean fitness will increase by qs. Therefore ν ≈ (qs)2/[2 ln(Nqs)].

The other factor is the dynamics of the most-fit subpopulations (the “nose” of the fitness distribution), where new mutations are essential. A more-fit mutant that moves the nose forward by s will be established a time τq after the previous most-fit mutant. Thus the nose advances at a speed ν = s/[left angle bracket]τq[right angle bracket], where [left angle bracket]τq[right angle bracket] is the average τq. After it is established, the fittest population nq will grow exponentially at rate qs and produce mutants at a rate UbnqUb1qseqst. Many new mutants will establish soon after the time τ at which Ubqs0τnq(t)dt=1, so the time it takes a new mutant to establish is τq1qsln(sUb). This means the nose advances at rate ν = s/[left angle bracket]τq[right angle bracket] ~ qs2/ ln(s/Ub). Yet we argued above that the bulk of the population fixes the speed of the mean via the selection pressure: ν ≈ (qs)2/[2 ln(Nqs)]. In steady state, the speed of the mean must equal the speed of the nose — the mutation-selection balance. This implies that q2ln[Ns]ln[sUb] and v2s2ln[Ns]ln2[sUb].

This crude argument neglects some important details of the stochastic process at the nose, but the basic qualitative behavior follows from this intuitive reasoning. We present a more detailed and careful analysis (and simulations to test the theory) elsewhere [30], and find



These are the two key predictions we test in this work. In comparing with experiments, we must also account for transient effects, which lead to slower adaptation before the steady state mutation-selection balance is reached. These transient effects are greater in larger populations; details of this effect are presented elsewhere [30].

The results in Box 1 reflect the simplest multiple-mutations model (Fig. 1b), which assumes that all mutations have the same effect, ŝ, so that one-by-one clonal interference is absent by definition. Others have recently studied similar models in regimes relevant to other situations [34, 35], and earlier work [13] took the first steps in the analysis described above but did not correctly account for stochastic effects that control the timing between the establishment of successive mutations, or for the balance between mutation and selection.

One-by-one clonal interference and the simplest multiple-mutation analyses are both incomplete. Each neglects the important effect that dominates the other. Neglecting multiple mutations never gives a complete picture of the dynamics, because whenever one-by-one clonal interference is important, so are multiple mutations. However, our multiple-mutation analysis can partially account for the effects of one-by-one clonal interference, because in many situations mutations with a small range of fitness effects around some value s dominate the evolution (Fig. 1c) [30]. Mutations much smaller than s occur frequently, but grow too slowly to interfere with mutations of size s (i.e. they are wasted due to clonal interference with the mutations of size s). And, unless the distribution of mutational effects ρ(s) falls off very slowly with s (slower than 1/s3), mutations much more beneficial than s happen rarely enough to have little overall impact. The multiple-mutation analysis cannot predict s, as this depends in a subtle way [30] on the unknown distribution of mutational effects and the resulting clonal interference processes. However, if we take s as a parameter to fit from experiments, the multiple-mutation theory implicitly accounts for one-by-one clonal interference effects, provided we redefine Ub to be the mutation rate towards beneficial mutations of roughly this size. However, if ρ(s) falls off slower than exponentially with s, s will depend significantly on N and Ub, and the behavior is more complicated.

Experimental Results and Discussion

Our multiple-mutation analysis predicts that the speed of evolution and the variation in fitness within a population both increase logarithmically with N and with Ub. It also shows that neglecting multiple mutations is a serious flaw of one-by-one clonal interference analyses. We set out to test these predictions by evolving asexual populations of diploid budding yeast in glucose-limited media for 500 generations at three different effective population sizes ranging from 1400 to 3.5 × 106, each with two different mutation rates: “nonmutator” populations and msh2Δ “mutator” populations with Ub estimated to be 10 times higher [36] (based on the elevation of mutation rate at two particular loci, and hence only a rough estimate; see Methods). The generation time in this media was initially about 130 minutes, compared to 90 minutes in rich media.

We periodically measured the fitness of each entire population by mixing a sample of it with a derivative of the ancestral strain that had been labelled with yellow fluorescent protein, growing the mixed population for 20 generations, and determining the ratio of the two strains at the beginning and end of the assay by using flow cytometry to distinguish labelled from unlabelled cells (see Methods). The total fitness changes over the 500 generations were used to obtain the average speed of evolution (Fig. 2a). We also measured the distribution of fitnesses within some of the evolved populations by isolating 96 individuals from each and then measuring their individual fitnesses (Fig. 3).

Figure 2
(a) Average speed of adaptation (ν) of our experimental populations, in percent fitness increase per generation, versus ln[NUb] (scaled to the smallest nonmutator population). N is the effective population size (which takes on three values, N ...
Figure 3
Fitness distributions of evolved populations after 500 generations, for three independently evolved lines of the: (a) largest mutator populations, (b) smallest mutator populations, (c) largest nonmutator populations, (d) smallest nonmutator populations. ...

Our data are clearly inconsistent with the simple successional fixation prediction, ν linear in NUb (p < 0.001). Other simple interpretations are ruled out by the observed time dependence of the mean fitness of our populations (Fig. 2b). The rate of fitness increase is roughly constant, in particular showing no evidence of slowing down as the experiment progresses (if anything, a slight speeding up is seen). This indicates that neither antagonistic epistasis (i.e. the combined effect of two beneficial mutations being less than the sum of their separate effects) nor a limited supply of beneficial mutations (i.e. “running out” of beneficial mutations) can be responsible for the observed weak dependence of ν on N. Note that the batch culture environment remains the same throughout our experiment, and the populations are in exponential phase throughout, so environmental changes cannot explain these results.

Our data are also inconsistent with one-by-one clonal interference analyses, because of their assumption that mutations fix singly in succession. This inconsistency is most apparent for our largest populations. Our large nonmutator populations increased in fitness by about 4 to 7 percent in 500 generations. This is not enough time for two or more mutations adding up to 4 to 7 percent to fix one by one (i.e. successionally). For example, two 3.5 percent effect mutations would take a minimum of 1000 generations to fix successionally; all other combinations adding to between 4 and 7 percent would take even longer. A similar argument applies to our large mutator populations.

Thus if single beneficial mutations fix successionally, one large such mutation must be responsible for almost the entire observed fitness increase. However, this is also inconsistent with the data. Fig. 2b shows that the mean fitness of our populations increases smoothly, and the individual profiles are similar to each other. Both features imply that the evolution is not dominated by single large mutations. If it were, the mean fitness would remain constant for a time and then rapidly increase by the amount of the large mutation. For a 7 percent mutation, for example, most of the increase in fitness would occur in just 30 generations (fixation times are much longer because mutants spend a long time while rare). The fitnesses of different populations would also show a wide range of kinetics depending on whether their large effect mutations occurred early or late (Fig. 4a). Yet this is not at all what we see. Instead, the gradual increase in fitness and similar kinetics between lines strongly suggest that many smaller mutations are steadily accumulating (Fig. 4b). This cannot happen unless multiple mutations sweep together: successional sweeps of small effect mutations would take far too long (Fig. 4c). A more detailed discussion, including other inconsistencies with one-by-one clonal interference and special circumstances in which one-by-one clonal interference could produce the observed results, is presented in the supplemental material.

Figure 4
Computer simulations showing the kinetics of the increase in mean fitness in 10 simulated populations of size N = 3.5 × 106, the same as our large experimental populations. (a) Assuming a single 7% effect mutation is responsible for the evolution ...

The above arguments suggest that the multiple-mutation analysis is the correct explanation for our results. A key qualitative prediction of this analysis is that the width of the fitness distributions in large populations should be greater than in small populations. In contrast, one-by-one clonal interference predicts that fitness distributions will show pronounced fluctuations over time for any population size: narrow and dominated by a single clone most of the time or, if measured during a selective sweep, clearly bimodal. Similar behavior would arise from a simple successional fixation (one-locus) analysis, except that the rate of sweeps would increase dramatically in large populations. In actuality, for both mutators and nonmutators, we find that the fitness distributions of large populations are broader than of small ones (Fig. 3a-d).

These predictions can be made quantitative: the expected widths of the fitness distributions and the speeds of evolution are given by the formulas in Box 1. These predictions depend on just two unknown parameters: the typical size of the beneficial mutations responsible for the fitness increase, s, and the rate at which these beneficial mutations occur, Ub. One might worry that for any experimental data, there would be a combination of Ub and s that would produce a good fit. This is not so. In the smallest populations, NUb is so small that they can only be in the successional fixation regime. These populations tightly constrain Ub and s in a way that is independent of the multiple-mutation theory, ruling out arbitrary Ub and s which might have yielded good fits to the other data. Within these constraints, we fit Ub and s from the data, yielding values Ub = 2.4 × 10−4 for mutator populations (hence Ub = 2.4 × 10−5 for nonmutators), and s = 2 percent. Details of the theory-independent constraints and the fit to data are described in the Methods. The resulting comparison between theory and experiment is summarized in Figs. Figs.2a2a and and3h.3h. The predicted increases in mean fitness (which give the speeds of evolution shown in Fig. 2a) and widths of the fitness distributions are each within a single fitness increment s of the experiments — as accurately as theory could possibly predict. There are, however, small systematic discrepancies: the theory overestimates the mean fitness increases for mutator populations and underestimates their width, while making the opposite errors in nonmutators. This is likely due to deleterious mutations, which we now consider.

Deleterious mutations complicate the shapes of the fitness distributions. However, their effects are most pronounced on the less-fit side of the distributions: on the more-fit side, all the clones are depleted similarly by deleterious mutations and the modifications of the shape of the distribution are small. Thus in the analysis described above we only use the more-fit side, above the median. But deleterious mutations may indirectly a ect the more-fit side of the fitness distributions, for example by decreasing the median fitness. They will also cause a reduction in the mean fitness of the population and hence could reduce the observed speed of evolution.

Unfortunately, the effects of deleterious mutations depend on the unknown distribution of their fitness decrements, so precise predictions are impossible. We can, however, estimate their maximum impact by looking at the small population fitness distributions. The small nonmutator fitness distributions are no wider to the right than expected from the measured experimental errors (Fig. 3d, e), which implies that deleterious mutations do not significantly reduce the mean fitness in nonmutators, nor do they a ect the above-median fitness distribution width. In other words, the contribution of deleterious mutations in nonmutators is minor and limited to a slight increase in the width of the less-fit tail. In mutators, on the other hand, the width to the right of the median in the small populations could be entirely due to reduction of the median by deleterious mutations, entirely due to beneficial mutations, or due to some combination of the two. This means that in all the mutator populations, deleterious mutations may decrease the mean fitness by at most two percent and broaden their fitness distributions by convolving them (defined in [37]) with a distribution of standard deviation 1.4 percent. These shifts lead to slight changes in the best-fit Ub and s which lead to the opposite shifts in the predicted results for nonmutators. These corrections roughly account for the systematic discrepancies between experiments and the multiple-mutation theory.

Although the multiple-mutation picture better explains our experiments, clonal interference must nevertheless also occur. Mutations of very small effect are certainly being regularly wasted, and this process partially determines the typical size, s, of the mutations that dominate the evolution. However, our data indicates that the effect omitted in one-by-one clonal interference analyses, the accumulation of multiple mutations, is crucial. After fitting s from data to implicitly account for clonal interference effects, the simple multiple-mutation theory is consistent with our experiments, especially once we consider the additional effects of deleterious mutations.

Several other recent experimental studies have also found that, as in our experiments, the speed of adaptation increases less than linearly with population size and mutation rate [4, 1619]. This has been taken as support for one-by-one clonal interference. But our multiple-mutations analysis also predicts a specific form of this less than linear dependence on N and Ub, albeit for different reasons. This earlier work is not suffciently detailed to distinguish between one-by-one clonal interference and our multiple mutations model.

If the beneficial mutations of size s ≈ 2% are point mutations, combining the estimate of Ub with the per base pair mutation rate of order 10−9 per generation [38] suggests that the target size for beneficial mutations in our experiments is a few thousand base pairs. This is substantially higher than the beneficial mutation rates in several earlier studies done in different environments [39, 40], but closer to recent estimates by Joseph and Hall [41]. It is possible that there are several targets of roughly a hundred base pairs, such as genes where inactivating one of the two copies in a diploid conveys an advantage, or a number of much smaller mutational hot spots (as found by [42]), whose mutation rate is much higher than the average per base pair mutation rate — perhaps having evolved to allow rapid mutational switches between different metabolic environments encountered in the natural history of budding yeast.

Finally, we note that the logarithmic increase in the speed of evolution with N and Ub in the large-NUb multiple-mutations regime is dramatically slower than the linear successional mutations regime. The difference has many implications. For example, the potential advantage of sex in combining mutations from different lineages becomes more pronounced in large populations, while mutator phenotypes become less useful as population sizes increase.

Experimental Procedures

Experimental Evolution Protocol

The experimental lines were established from a single W303 diploid with a heterozygous deletion at both ime1 (required for sporulation) and msh2 (this deletion elevates the mutation rate). The ime1 deletion was created by amplifying from plasmid pFA6-kanMX4 using primers F1 (5′-GAAAAA AATAAT AAAAGA AAAGCT TTTCTA TTCCTC TCCCCA CAAACA AAGGTC GACGGA TCCCCG GGTT-3′) and R1 (5′-AATGGA TATATT TTGAGG GAAGGG GGAAGA TTGTAG TACTTT TCGAGA ATCGAT GAATTC GAGCTC GTT-3′), creating a deletion disruption cassette used to delete ime1 and replace it with KANmx6, which confers resistance to G418 [43]. The msh2 deletion cassette was amplified from plasmid pAG25 using primers F1 (5′-AAAAAT CTCTTT ATCTGC TGACCT AACATC AAAATC CTCAGA TTAAAA GTGGTC GACGGA TCCCCG GGTT-3′) and R1 (5′-ATCTAT ATATTA TCTATC GATTCT CACTTA AGATGT CGTTGT AATATT AATCGA TGAATT CGAGCT CGTT-3′), creating a deletion disruption cassette used to delete msh2 and replace it with NAT 1, which confers resistance to Clonat [44]. This resulted in a genotype matαmataime1Δ::KANmx6IME1msh2Δ::NAT1MSH2URA3ura3CAN1can1ade2ade2. This was sporulated, and the spore clones mated to create two independent homozygous ime1Δ msh2Δ “mutators,” and two independent homozygous ime1Δ “nonmutators.” We refer to all these clones as asexuals, as they are unable to sporulate. Each genotype was used to found 9 lines, 3 lines at each of 3 different effective population sizes, from two independent diploids. All initial lines were approximately clonal. The elevation of Ub in the mutators is impractical to measure, as it depends only on the unknown loci where beneficial mutations can occur. We thus assume that mutators have Ub of order 10 times that of nonmutators, as indicated by previous measurements of msh2Δ cells at CAN1 and URA3 [36]. This is only a rough estimate, but our theoretical predictions are fairly insensitive to this.

Each initial population was grown overnight in YEPEG (YEP+2.5% ethanol, 2% glycerol, where YEP is 2% yeast extract, 2% peptone, 0.0025% adenine and tryptophan)[45] to eliminate possible petite phenotypes, and then placed in the selective media at the appropriate bottleneck population size to begin the experiment. Selection was performed in YEP+0.05% dextrose, ampicillin at 100 μg/mL and tetracycline at 25 μg/mL (LG). Cultures were grown in 15 mL of LG in 50 mL test tubes, continuously mixed in roller drums at 30 °C until they reached a density of about 1.5×107 cells/mL, half the saturation density in this media. Each was then counted with a Coulter counter and an appropriate amount of media transferred to fresh LG. The amounts transferred were chosen to achieve bottleneck population sizes of Nb = 100, Nb = 7100, and Nb = 5.5 × 105.

The differing bottleneck sizes in the same volume of media required different numbers of generations G between transfers, G ≈ 21, G ≈ 14.5, and G ≈ 8.5 for small, medium, and large populations respectively. This gives effective population sizes separated by factors of 50: N = 1.4 × 103, N = 6.9 × 104, and N = 3.5 × 106 (see below).

The serial dilution process was carried out for approximately 500 total generations of all of the lines. Beginning from the same original starting lines, the entire experiment was then repeated. Aliquots from the evolving cultures were frozen periodically throughout the experiment in LG+7.5% glycerol.

This experiment was originally intended to compare sexual and asexual populations. Accordingly, all lines were periodically put through sporulation cycles. The asexual lines we discuss in this paper are unable to sporulate, but for consistency experienced the same sporulation conditions as sexual lines. In each sporulation cycle, the cells were transferred at 1.5 × 107 cells/mL to YEPA (YEP+2% potassium acetate)[45]. They were grown in YEPA for 6 hours, then in SPM (2% potassium acetate, 0.02% ra nose, 0.0005% adenine, ampicillin at 100 μg/mL and tetracycline at 25 μg/mL)[45] for 48 hours, placed on YPD (YEP+2% dextrose)[45] plates to mate overnight and then returned to LG. Total growth during this process was about 2 generations. In the first 500-generation experiment, sporulation cycles were carried out every 100 generations. In the second, independent experiment, they were carried out every 60 generations. The fitness distributions shown are all from the latter populations.

Fitness Measurements

To measure fitness, we competed the evolved strains against reference strains. These reference strains were constructed from spores from the original starting line, transformed to create Hygromycin-resistant strains (for colony-counting assays) and a Hygromycin-resistant strain expressing yellow fluorescent protein (YFP, for fluorescence activated cell sorting (FACS) assays). This reference strain was constructed using a pVENUs plasmid containing YFP linked to a HIS3 marker, generously provided by Kurt Thorn (Bauer Center for Genomics Research, Harvard University). A 1433-bp fragment containing a Hygromycin resistance cassette was obtained by digesting pAG32 with EcoRV and Mlu1. The HIS3 marker was then removed from pVENUs by digesting this plasmid with EcoRV and Mlu1 to remove the 1055-bp fragment containing the HIS3 coding region. The Hygromycin resistance cassette was ligated with the resulting 3671-bp vector fragment containing a YFP coding region to generate pVEN2. The Hygromycin-resistant YFP strain was created by amplifying from pVEN2 using primers F1 (5′-CCGCTG TCGGTA TGGGTG CCGGTG CTCTAG CTGCTG CTGCTA TGTTGT TAGGTC GACGGA TCCCCG GGTT-3′) and R1 (5′-CGAAAA TTTTGA AAAAAG CCATAT AGATAT TATAAA AAATCA GAGATT TCTCGA TGAATT CGAGCT CGTT-3′), creating a cassette used to fuse the YFP, linked to Hygromycin resistance, to cwp2 [44]. Although Cwp2 is a cell wall protein, the Cwp2-YFP fusion is cytoplasmic.

To measure the mean fitness of an evolved population, approximately 10 μL of a reference strain was unfrozen onto a YPD plate and left to grow overnight. These cells were then placed in LG overnight to reacclimate. At the same time, about 10 μL of the evolved strain was unfrozen and placed in LG overnight. These two lines were then mixed at a density of order 1 × 105 cells/mL, and the ratio ri of the number of reference to evolved individuals was measured. The cells were grown to about 1.5 × 107 cells/mL, and then 10 μL of the culture was transferred to new media. They were again allowed to grow and diluted. After growing a third time to about 1.5 × 107 cells/mL (a total of about t = 19 generations of growth) the ratio of evolved to reference rf was measured. The fitness difference was defined to be Δs=1tln[rirf]. The fitnesses of the starting lines were simultaneously measured to determine the initial fitness advantage or disadvantage of the reference strain. These were subtracted from the results to yield comparisons between the evolved and starting fitness.

To measure the ratio of reference to evolved populations by colony counting (used for measuring the mean fitness in the first of the 500-generation experiments), we diluted and plated the mixture of the two strains on YPD plates. After individual colonies were visible on the plates, we replica-plated these plates to YPD and Hygromycin plates (YPD+Hygromycin B at 300 μg/mL [44]) and counted the number of hygromycin-resistant reference individuals and hygromycin-sensitive evolved individuals. These counts were done in replicate on 4 plates, to give a total of 500 to 800 colonies counted. To measure the ratio by FACS (used for all other measurements), an aliquot of the mixture of the two strains was washed out of LG, placed in PBS, sonicated, and then analyzed with a MoFlo FACS machine (DakoCytomation, Carpinteria, CA). 5 × 103 to 1 × 105 cells per mixture were analyzed. The cells clearly divide into YFP-labelled and unlabelled individuals, yielding the ratio of reference to evolved cells.

To measure the fitness distribution within a population, we unfroze that population and diluted and plated it to single colonies on a YPD plate. 96 randomly chosen individual colonies were picked from this plate and placed in LG overnight to reacclimate. These were then mixed with a reference strain (also plated and placed in LG to reacclimate overnight) and the fitness measurements proceeded as above.

Fits to Theoretical Results

The less-fit sides of the fitness distributions are heavily influenced by deleterious mutations. These bias the shapes, particularly of this side of the distributions. In addition, we occasionally observe a single individual 8 to 10 percent less fit than the mean, presumably because of a large-effect deleterious mutation. This means that the standard deviation, or any other statistic that depends on the whole distribution, can be heavily influenced by deleterious mutations. Therefore in making statistical comparisons between fitness distributions and with theory, we use only the more-fit half of each distribution. We define the center of the distribution to be the median, xm, (which is less influenced by outlying deleterious mutations than the mean) and use the second moment to the right of the median as a measure of the width, σ>2=2MΣi(xixm)2, where xi are the fitnesses of the M/2 measurements above the median. To test the theory, we compare values of σ> between large and small populations at the same mutation rate. Other measures of the width give similar results.

In the parameter ranges relevant for this experiment, the theory predicts an approximately Gaussian fitness distribution (except in the tails where the number of individuals is too small to be sampled anyway). This means we can relate the theoretical q (see Box 1) to the observed σ>2 by


in fitting to the multiple-mutation theory. The lead of the distribution is much larger in a large population than a small one given the same σ>, because one percent of a large population represents many more individuals than one percent of a small population. Thus the difference of the leads between large and small populations is much greater than the difference in σ>.

To fit the theory to the data, we calculated the mean speeds of adaptation and leads of fitness distributions from our experiments, and estimated the standard error in these values from the inter-line variation. We defined the best-fit values of s and Ub to be those that minimized the squared deviations of the theoretical predictions from the experimental values, measured in units of the standard errors. This fit was done within the constraints imposed by the smallest NUb populations. These constraints depend only on bounding arguments, and are independent of the multiple-mutation theory. The one-locus theory result is an upper bound on the speed of evolution at given parameters. Thus if we assume that the smallest mutator populations are in the one-locus regime, for these populations to have evolved as much as observed yields an estimated lower bound on Ub and s (we cannot use the smallest nonmutator populations because these are consistent with not having acquired any beneficial mutations at all). This implies that in order to achieve the mean fitness gains seen in the small mutator populations, and to explain the between-line variations, Ub in mutators must be at least of order 10−4.5 and s at least about 0.01. At the same time, the smallest NUb populations evolved by at most a few percent. This sets an upper bound on s of a few percent, and given the lower bound on s, it also implies that Ub in mutators is at most of order 10−3.5. Since the multiple-mutation theory depends only logarithmically on Ub, this factor-of-ten range of Ub is a tight constraint.

There is no single effective population size relevant for comparison to theory. Rather, the population size enters in two ways: the time it takes selection to change the mean fitness, and the rate at which new mutations establish. The harmonic mean of the population size is relevant for the former, and Wahl and Gerrish [46] found that the effective size N = NbG ln 2 determines the latter. Only the latter is relevant for single-locus dynamics, because here fixation times are negligible. The two definitions are very similar for all comparisons to our theory made in this paper (which depend only logarithmically on N), so we use the latter formula throughout and neglect the difference. For the discussion of the inconsistency with clonal interference, we use the harmonic mean effective population size where appropriate (i.e. when discussing fixation times).

Supplementary Material



We are grateful to Dawn Thompson and Sheri Simmons for many useful discussions and help maintaining the lines, to John Wakeley for stimulating discussions, and to Larry Shumway, George Kenty, and Doug Melton for help with FACS. We thank members of the Murray lab for comments on the manuscript. This work was supported by grant GM 68763 from the National Institutes of Health.


Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.


[1] Crow JF, Kimura M. Evolution in sexual and asexual populations. American Naturalist. 1965;909:439.
[2] Smith J Maynard. Evolution in sexual and asexual populations. American Naturalist. 1968;102(927):469.
[3] Smith J Maynard. What use is sex? J. Theor. Biol. 1971;30(2):319. [PubMed]
[4] Colegrave N. Sex releases the speed limit on evolution. Nature. 2002;420(6916):664–666. [PubMed]
[5] Otto SP, Lenormand T. Resolving the paradox of sex and recombination. Nature Reviews Genetics. 2002;3(4):252–261. [PubMed]
[6] Tenaillon O, Le Nagard H, Godelle B, Taddei F. Mutators and sex in bacteria: Conflict between adaptive strategies. PNAS. 2000;97(19):10465–10470. [PubMed]
[7] Barton NH, Otto Sarah P. Evolution of recombination due to random drift. Genetics. 2005;169(4):2353–2370. [PubMed]
[8] Barton NH. Linkage and the limits to natural-selection. Genetics. 1995;140(2):821–841. [PubMed]
[9] Otto SP, Barton NH. The evolution of recombination: Removing the limits to natural selection. Genetics. 1997;147(2):879–906. [PubMed]
[10] Otto SP, Barton NH. Selection for recombination in small populations. Evolution. 2001;55(10):1921–1931. [PubMed]
[11] Cox EC, Gibson TC. Selection for high mutation rates in chemostats. Genetics. 1974;77(2):169–184. [PubMed]
[12] LeClerc JE, Li BG, Payne WL, Cebula TA. High mutation frequencies among escherichia coli and salmonella pathogens. Science. 1996;274(5290):1208–1211. [PubMed]
[13] Sniegowski PD, Gerrish PJ, Lenski RE. Evolution of high mutation rates in experimental populations of e-coli. Nature. 1997;387(6634):703–705. [PubMed]
[14] Taddei F, Radman M, MaynardSmith J, Toupance B, Gouyon PH, Godelle B. Role of mutator alleles in adaptive evolution. Nature. 1997;387(6634):700–702. [PubMed]
[15] Giraud A, Matic I, Tenaillon O, Clara A, Radman M, Fons M, Taddei F. Costs and benefits of high mutation rates: Adaptive evolution of bacteria in the mouse gut. Science. 2001;291(5513):2606–2608. [PubMed]
[16] de Visser Jagm, Zeyl CW, Gerrish PJ, Blanchard JL, Lenski RE. Diminishing returns from mutation supply rate in asexual populations. Science. 1999;283(5400):404–406. [PubMed]
[17] Miralles R, Moya A, Elena SF. Diminishing returns of population size in the rate of rna virus adaptation. J. Virology. 2000;74(8):3566–3571. [PMC free article] [PubMed]
[18] Miralles R, Gerrish PJ, Moya A, Elena SF. Clonal interference and the evolution of rna viruses. Science. 1999;285(5434):1745–1747. [PubMed]
[19] de Visser JAGM, Rozen DE. Limits to adaptation in asexual populations. Journal of Evolutionary Biology. 2005;18(4):779–788. [PubMed]
[20] Gerrish P. The rhythm of microbial adaptation. Nature. 2001;413(6853):299–302. [PubMed]
[21] Johnson T, Barton NH. The effect of deleterious alleles on adaptation in asexual populations. Genetics. 2002;162(1):395–411. [PubMed]
[22] Gerrish Philip, Lenski Richard. The fate of competing beneficial mutations in an asexual population. Genetica. 1998;102/103:127–144. [PubMed]
[23] Wilke CO. The speed of adapation in large asexual populations. Genetics. 2004;167:2045–53. [PubMed]
[24] Kim Y, Stephan W. Selective sweeps in the presence of interference among partially linked loci. Genetics. 2003;164(1):389–398. [PubMed]
[25] Campos PRA, De Oliveira VM. Mutational effects on the clonal interference phenomenon. Evolution. 2004;58(5):932–937. [PubMed]
[26] Orr HA. The rate of adaptation in asexuals. Genetics. 2000;155(2):961–968. [PubMed]
[27] More generally, there is a probability As that a mutant lineage will survive drift, where A is a constant of order 1 which depends on the specific stochastic model of the dynamics. Our theoretical results are based on a continuous-time branching process model where A = 1, but the value of A is unimportant in the analysis presented here See Johnson T, Gerrish PJ. The fixation probability of a beneficial allele in a population dividing by binary fission. Genetica. 2002;115(3):283–287. [PubMed]
[28] Hill WG, Robertson A. Effect of linkage on limits to artificial selection. Genetical Research. 1966;8(3):269. [PubMed]
[29] Felsenstein Joseph. The evolutionary advantage of recombination. Genetics. 1974;78:737–756. [PubMed]
[30] Desai Michael M., Fisher Daniel S. Beneficial mutation-selection balance and the effect of linkage on positive selection. Submitted to Genetics. 2006 [PubMed]
[31] Yedid Gabriel, Bell Graham. Microevolution in an electronic microcosm. The American Naturalist. 2001;157(5):465–487. [PubMed]
[32] Bachtrog Doris, Gordo Isabel. Adaptive evolution of asexual populations under muller′s ratchet. Evolution. 2004;58(7):1403–1413. [PubMed]
[33] Shaver AC, Dombrowski PG, Sweeney JY, Treis T, Zappala RM, Sniegowski PD. Fitness evolution and the rise of mutator alleles in experimental escherichia coli populations. Genetics. 2002;162(2):557–566. [PubMed]
[34] Rouzine Igor, Wakeley John, Coffin JM. The solitary wave of asexual evolution. PNAS. 2003;100(2):587–592. [PubMed]
[35] Kim Yuseob, Orr H. Allen. Adaptation in sexuals vs. asexuals: clonal interference and the fisher-muller model. Genetics. 2005;171:1377–1386. [PubMed]
[36] Lang Greg, Murray Andrew W. unpublished results. 2005
[37] Feller William. An Introduction to Probability Theory and Its Applications. Wiley; United States: 1968.
[38] Drake JW. A constant rate of spontaneous mutation in dna-based microbes. PNAS. 1991;88(16):7160–7164. [PubMed]
[39] Lenski RE, Rose MR, Simpson SC, Tadler SC. Long-term experimental evolution in escherichia-coli .1. adaptation and divergence during 2,000 generations. American Naturalist. 1991;138(6):1315–1341.
[40] Zeyl C, Vanderford T, Carter M. An evolutionary advantage of haploidy in large yeast populations. Science. 2003;299(5606):555–558. [PubMed]
[41] Joseph Sarah B., Hall David W. Spontaneous mutations in diploid saccharomyces cerevisiae: More beneficial than expected. Genetics. 2004;168(4):1817–1825. [PubMed]
[42] Dunham Maitreya J., Badrane Hassan, Ferea Tracy, Adams Julian, Brown Patrick O., Rosenzweig Frank, Botstein David. Characteristic genome rearrangements in experimental evolution of saccharomyces cerevisiae. Proceedings of the National Academy of Sciences. 2002;99(25):16144–16149. [PubMed]
[43] Wach A, Brachat A, Pohlmann R, Philippsen P. New heterologous modules for classical or pcr-based gene disruptions in saccharomyces cerevisiae. Yeast. 1994;10(3):1793–1808. [PubMed]
[44] Goldstein AL, McCusker JH. Three new dominant drug resistance cassettes for gene disruption in saccharomyces cerevisiae. Yeast. 1999;15(14):1541–53. [PubMed]
[45] Sherman F, Fink G, Lawrence C. Methods in Yeast Genetics. Cold Spring Harbor Laboratory Press; Cold Spring Harbor, NY: 1974.
[46] Wahl LM, Gerrish PJ. The probability that beneficial mutations are lost in populations with periodic bottlenecks. Evolution. 2001;55(12):2606–2610. [PubMed]