Determining the probability of fixation of beneficial mutations is critically important for building predictive models of adaptive evolution. Despite considerable theoretical work, models of fixation probability have stood untested for nearly a century. However, recent advances in experimental and theoretical techniques permit the development of models with testable predictions. We developed a new model for the probability of surviving genetic drift, a major component of fixation probability, for novel beneficial mutations in the fungus Aspergillus nidulans, based on the life-history characteristics of its colony growth on a solid surface. We tested the model by measuring the probability of surviving drift in 11 adapted strains introduced into wild-type populations of different densities. We found that the probability of surviving drift increased with mutant invasion fitness, and decreased with wild-type density, as expected. The model accurately predicted the survival probability for the majority of mutants, yielding one of the first direct tests of the extinction probability of beneficial mutations.
A fundamental theme of theoretical population genetics that has received little empirical attention is the fixation probability of a beneficial mutation (), the likelihood that the mutation will ultimately be present in every individual of the population. First addressed by Haldane, Fisher and Wright nearly 90 years ago, the best-known theoretical result is Haldane's approximation for , , for populations of constant size undergoing weak selection [1–3]. Despite a sustained interest in theoretical models, strict model assumptions and an inability to track individual mutations have prevented empirical tests of their predictions. However, recent advances in theoretical and experimental techniques permit us now to develop models with predictions that can be tested directly .
The first component of is determining the probability that novel beneficial mutations survive genetic drift, the random loss of mutants while rare. The processes driving genetic drift are specific to the organism and its growth environment, which must be incorporated into fixation models through a thorough understanding of life-history traits. Although most models (see ) use relative mutant fecundity as a predictor for surviving drift, recent work has shown that other measures of fitness may be more applicable (such as survival or generation time) [5,6]. We developed and tested a model for the probability of surviving drift for mutations arising in the fungus Aspergillus nidulans based on the life history of a radially growing colony , with mycelial growth rate  and initial frequency as model parameters. To test the model, we introduced individual mutant spores with a known fitness advantage into populations of wild-type spores of a different colour, and then visually determined the fate of the mutant. This approach let us separate loss due to drift into components with different biological bases. Our model accurately predicted the empirical probability of surviving drift; this is one of the first direct tests of any model of fixation probability.
(a) Experimental protocol
For growth media, we used minimal medium (MM) with pH 5.8, consisting of: NaNO3 6.0 g l−1; KH2PO4; 1.5 g l−1; MgSO4 · 7H2O; 0.5 g l−1; NaCl 0.5 g l−1; agar 12 g l−1 and 1 mg l−1 each of FeSO4, ZnSO4, MnCl2 and CuSO4, and complete medium (CM), which is MM plus: tryptone 10 g l−1; yeast extract 5 g l−1 and vitamin solution (after autoclaving) . After autoclaving, glucose was added to MM at 0.1 g l−1 (0.01%) or to CM at 4 g l−1 (0.4%).
We used 11 mutant strains of A. nidulans with beneficial mutations resulting from a long-term evolution experiment in which spore samples were transferred weekly to MM after approximately 50 nuclear generations. These were adapted from two ancestral strains with either green or yellow spores. The growth rate of each strain in millimetre per hour was determined by measuring colony diameter after 5 days of growth on CM; these were used to obtain invasion fitness (see §2b for definition).
We estimated the probabilities of escaping the inoculum (pe) and forming a sector () by inoculating an average of one mutant spore with either 100 or 1000 wild-type spores in a 5 µl droplet on CM. The number of mutant spores was Poisson distributed (λ = 1), which we took into account in tests of the model (see the electronic supplementary material, §3). Forty replicate populations were used per strain-treatment combination, totalling 880 populations. Colonies were photographed (Panasonic DMC-FZ8) after 3 days of growth at 37°C and storage at 4°C for between 1 and 7 days. We categorized colonies into four groups by the presence of mutant mycelium: ‘escaped inoculum (with sector)’, ‘escaped inoculum (without sector)’, ‘blocked in inoculum’ or ‘undetectable’ (figure 1). Proportion of colony area occupied by mutant mycelium was measured with ImageJ (v. 1.44) using ‘Area Fraction’ . Data are available through Dryad (http://dx.doi.org/10.5061/dryad.0rc67).
Extinction probabilities are extremely sensitive to life-history and experimental details . To match theory to experiment, we consider a wild-type fungus growing radially in a circular colony with radial growth rate w. The colony begins as an inoculum in the centre of the plate, containing wild-type spores. The inoculum also contains one mutant spore with growth rate m > w.
Since is large, wild-type ‘sub-colonies’ growing radially from individual spores are likely to physically trap the mutant sub-population within the inoculum. In the electronic supplementary material, we compute the probability that the wild-type hyphae will ultimately block the mutant sub-colony from reaching the inoculum boundary. Assuming that the mutant hyphae will take the ‘shortest path’ to the boundary, we derive a lower bound for the probability that the mutant escapes the inoculum: where denotes the invasion fitness, a measure arising from the geometry of colony growth (see the electronic supplementary figure S1 and §S1):
Note that for large-effect mutations (), is very close to the relative fitness , while for small mutational effects .
This approximation for pe gives a lower bound because it assumes the mutant sub-colony escapes along the shortest path. In the electronic supplementary material, we relax this assumption, yielding a less elegant but more accurate numerical estimate of pe.
We also developed a simulation approach in which a circular plate is modelled as a grid of small pixels. Initially, wild-type spores are distributed randomly across the inoculum area, and a single-mutant spore is added at a random position within the inoculum. Each spore grows radially as a sub-colony until it is surrounded by occupied pixels, at which point it is occluded and ceases growing.
In §3, we demonstrate that, when Ni = 1000, most mutant spores that escape the inoculum successfully form sectors. However, it is unclear whether forming a sector guarantees fixation in a serial transfer regime. Letting A denote the fractional area of the beneficial sector, we also model the serial transfer process in the electronic supplementary material and predict that whenever the product is sufficiently large (i.e. ), after serial transfer is well-approximated by pe.
For the 11 mutants, the probability of sector formation () increased with invasion fitness () and decreased with wild-type frequency () (logistic regression: , for both parameters: , p < 0.001; figure 2a). The probability of escaping the inoculum (pe) also increased with and decreased with (logistic regression: , for both parameters: , p < 0.0001; figure 2c,d), as did average sector size (ANCOVA: model F2,118 = 10.2, adjusted R2 = 0.13, p < 0.0001; figure 2b).
Logically, was smaller than pe for both (one-tailed paired t-test: t10 = 4.96, p < 0.0005) and (one-tailed paired t-test: t10 = 2.28, p < 0.05). However, for , seven of 11 mutant strains had identical and pe, and thus the difference was driven by four mutant strains which had large proportions of occluded sectors (see §4). The conditional probability of forming a sector, given that an escape had occurred (), increased with , but the effect of was conditional on : for , the relationship between and was positive; for , it was negative (logistic regression: , main effects and , , ; interaction effect z =−1.85, p = 0.06).
Figure 2 also demonstrates that, for spores/inoculum, pe was well-approximated only by the simulation results (dotted line), as expected given the model assumptions that is large and the mutant colonizes a straight path to the inoculum boundary. For , the analytical and numerical models (dashed and solid lines, respectively) provided good estimates of pe, with four exceptions. Grey lines show the classical prediction for fixation probability in a haploid population of constant size , corrected for the possibility of multiple mutant spores per inoculum (see the electronic supplementary material).
The product of initial wild-type frequency and average sector area, , ranged from 9.22 to 13.3 for , and from 66.8 to 238 for . Thus, when , suggesting that pe is a good estimator of the ultimate fixation probability after serial transfer for these populations.
Long-standing models for the probability of escaping genetic drift have remained untested owing to a mismatch between appropriate empirical data and model assumptions. We developed and tested a new model for the probability of surviving genetic drift for beneficial mutations in the fungus A. nidulans based on a thorough understanding of its life-history characteristics during colony expansion on a surface. Our model incorporated initial frequency and invasion fitness and accurately predicted the probability of surviving genetic drift for 11 beneficial mutants.
While genetic drift is typically characterized by the random sampling of alleles , in our system drift arises primarily from the random position of spores within the inoculum. If we consider inocula as ‘generations’ in a serial transfer regime, the random sampling of alleles from one generation to the next occurs largely through the random placement of spores in the inoculum, since spores that are close to the boundary have a greater chance of contributing to the next generation. Thus, our model captures the physical mechanism underpinning genetic drift. We also identified a second component of drift in our system, which blocked mutants after they had escaped the inoculum. Experimentally, we were able to observe the effect of and on these components independently, which led to the discovery of a surprising dynamic between and drift after escaping the inoculum. Although the overall probability of escaping drift decreased with as expected, the conditional probability of sector formation () increased with , meaning the chance of being blocked outside the inoculum was marginally greater for than . This results from an increased number of ‘circuitous escapes’ permitted by the lower wild-type density, where the mutant does not take the shortest path to escape the inoculum. These mutants tend to form narrowed, irregularly shaped growing fronts (as in figure 1b) that presumably result from single hyphae and are thus more easily occluded by the wild-type. The probability of a circuitous escape is also expected to increase with , but this effect likely depends on , i.e. there is a minimum that permits a circuitous escape for each . This is a possible explanation for the increased pe we observed in four of the populations, and for the negative relationship between and overall. These results suggest that the relationships between drift and fitness or initial frequency are not always straightforward—the spatial patterns of growth must also be considered.
Determining the probability of escaping genetic drift is the first component of a model for the ultimate fixation of beneficial mutations, . How much additional loss should be expected due to population bottlenecks if populations would be serially transferred until the mutant fixes or goes extinct? Our theoretical results suggest that pe is a good estimator of ultimate fixation probability for populations with large (where A is sector area). Furthermore, intuition from serial transfer experiments with this system suggests that visible sectors are unlikely to be lost due to bottlenecks, so that pe is likely to predict when . Populations from the satisfy both these requirements, suggesting that pe is a good estimate of for these populations. Further work will measure fixation probability directly and test its predicted dependence on pe and . Quantifying the relationship between the probability of fixation of genotypes and competitive fitness when common will aid our understanding of the variation in life-history traits among different organisms.
The classical prediction for  yielded a surprisingly good fit to our data when ; however, this prediction changes negligibly for and fits the data very poorly in this case. In contrast, our predictions vary sensitively with wild-type density, as do the experimental data. Although our model is specific to our experimental system, this work nevertheless demonstrates the predictive power of models of fixation probabilities that are tailored to the experimental system, and the power of current experimental techniques to examine the most fundamental assumptions of population genetics.
We thank Marijke Slakhorst for strains and technical assistance. This work was supported by NSERC grants to L.M.W. and D.R.G.
One contribution to a Special Feature on ‘Experimental evolution’ organized by Paul Sniegowski, Thomas Bataillon and Paul Joyce.
- Received April 5, 2012.
- Accepted June 6, 2012.
- This journal is © 2012 The Royal Society