Model and test in a fungus of the probability that beneficial mutations survive drift

Danna R. Gifford, J. Arjan G. M. de Visser, Lindi M. Wahl

Abstract

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.

1. Introduction

A fundamental theme of theoretical population genetics that has received little empirical attention is the fixation probability of a beneficial mutation (Embedded Image), 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 Embedded Image, Embedded Image, for populations of constant size undergoing weak selection [13]. 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 [4].

The first component of Embedded Image 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 [4]) 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 [7], with mycelial growth rate [8] 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.

2. Methods

(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) [9]. 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 (Embedded Image) 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’ [10]. Data are available through Dryad (http://dx.doi.org/10.5061/dryad.0rc67).

Figure 1.

Colonies containing mutant and wild-type spores after 3 days of growth. (a) Escaped inoculum (with sector)—mutant reaches colony boundary, (b) escaped inoculum (no sector)—mutant escapes inoculum but is occluded by wild-type (c) blocked in inoculum—mutant germinates but does not escape inoculum, and (d) undetectable—mutant is absent, does not germinate or fails to produce spores.

(b) Modelling

Extinction probabilities are extremely sensitive to life-history and experimental details [4]. 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 Embedded Image wild-type spores. The inoculum also contains one mutant spore with growth rate m > w.

Since Embedded Image 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: Embedded Image where Embedded Image denotes the invasion fitness, a measure arising from the geometry of colony growth (see the electronic supplementary figure S1 and §S1): Embedded Image

Note that for large-effect mutations (Embedded Image), Embedded Image is very close to the relative fitness Embedded Image, while for small mutational effects Embedded Image.

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 Embedded Image is sufficiently large (i.e. Embedded Image), Embedded Image after serial transfer is well-approximated by pe.

3. Results

For the 11 mutants, the probability of sector formation (Embedded Image) increased with invasion fitness (Embedded Image) and decreased with wild-type frequency (Embedded Image) (logistic regression: Embedded Image, for both parameters: Embedded Image, p < 0.001; figure 2a). The probability of escaping the inoculum (pe) also increased with Embedded Image and decreased with Embedded Image (logistic regression: Embedded Image, for both parameters: Embedded Image, 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).

Figure 2.

Both (a) the proportion of colonies with a mutant sector, Embedded Image, and (b) the average proportion of colony area made up by mutant mycelium, A, increase with mutant invasion fitness, Embedded Image and decrease with Embedded Image (black symbols, Embedded Image; grey symbols, Embedded Image). (c) For Embedded Image, the proportion of mutants escaping the inoculum droplet, pe (circles), closely matches simulation results and exceeds Embedded Image (triangles, from (a)). (d) For Embedded Image, observed pe (circles) closely matches modelled pe and observed Embedded Image (triangles, from (b)), barring four exceptional data points (see §4). Error bars are Embedded Images.e.; invasion fitness has a s.e. ≤ 0.029 for all points (not shown). Dashed line, analytical model; solid line, numerical model; dotted line, simulation results; grey line, classical prediction [11].

Logically, Embedded Image was smaller than pe for both Embedded Image (one-tailed paired t-test: t10 = 4.96, p < 0.0005) and Embedded Image (one-tailed paired t-test: t10 = 2.28, p < 0.05). However, for Embedded Image, seven of 11 mutant strains had identical Embedded Image 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 (Embedded Image), increased with Embedded Image, but the effect of Embedded Image was conditional on Embedded Image: for Embedded Image, the relationship between Embedded Image and Embedded Image was positive; for Embedded Image, it was negative (logistic regression: Embedded Image, main effects Embedded Image and Embedded Image, Embedded Image, Embedded Image; interaction effect z =−1.85, p = 0.06).

Figure 2 also demonstrates that, for Embedded Image spores/inoculum, pe was well-approximated only by the simulation results (dotted line), as expected given the model assumptions that Embedded Image is large and the mutant colonizes a straight path to the inoculum boundary. For Embedded Image, 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 Embedded Image [11], 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, Embedded Image, ranged from 9.22 to 13.3 for Embedded Image, and from 66.8 to 238 for Embedded Image. Thus, Embedded Image when Embedded Image, suggesting that pe is a good estimator of the ultimate fixation probability after serial transfer for these populations.

4. Discussion

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 [4], 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 Embedded Image and Embedded Image on these components independently, which led to the discovery of a surprising dynamic between Embedded Image and drift after escaping the inoculum. Although the overall probability of escaping drift decreased with Embedded Image as expected, the conditional probability of sector formation (Embedded Image) increased with Embedded Image, meaning the chance of being blocked outside the inoculum was marginally greater for Embedded Image than Embedded Image. 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 Embedded Image, but this effect likely depends on Embedded Image, i.e. there is a minimum Embedded Image that permits a circuitous escape for each Embedded Image. This is a possible explanation for the increased pe we observed in four of the Embedded Image populations, and for the negative relationship between Embedded Image and Embedded Image 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, Embedded Image. 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 Embedded Image (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 Embedded Image when Embedded Image. Populations from the Embedded Image satisfy both these requirements, suggesting that pe is a good estimate of Embedded Image for these populations. Further work will measure fixation probability directly and test its predicted dependence on pe and Embedded Image. 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 Embedded Image [11] yielded a surprisingly good fit to our data when Embedded Image; however, this prediction changes negligibly for Embedded Image 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.

Acknowledgements

We thank Marijke Slakhorst for strains and technical assistance. This work was supported by NSERC grants to L.M.W. and D.R.G.

Footnotes

  • 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.

References

View Abstract