help button home button Genetics J Exp Med
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Genetics Published Articles Ahead of Print on February 1, 2008.

Genetics, Vol. 178, 1571-1578, March 2008, Copyright © 2008
doi:10.1534/genetics.107.080853

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
genetics.107.080853v1
178/3/1571    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in Genetics
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Google Scholar
Right arrow Articles by Dickinson, W. J.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Dickinson, W. J.

Synergistic Fitness Interactions and a High Frequency of Beneficial Changes Among Mutations Accumulated Under Relaxed Selection in Saccharomyces cerevisiae

W. Joseph Dickinson1

Department of Biology, University of Utah, Salt Lake City, Utah 84112

1 Address for correspondence: Department of Biology, University of Utah, 257 South 1400 East, Salt Lake City, UT 84112.
E-mail: joebiohorn{at}msn.com

Manuscript received August 20, 2007. Accepted for publication December 31, 2007.


    ABSTRACT
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Spontaneous mutations were accumulated for ~4800 generations in 48 lines of yeast protected from effective selection by frequent passage through single-cell bottlenecks. Changes in fitness were evaluated by direct competition with matched parental stocks differing only at a selectively neutral marker locus. Average fitness declined by ~5% over the course of the experiment. The rate of change increased sharply in later generations, strongly suggesting synergistic epistasis. Divergence among lines increased rapidly relative to the change in average fitness and also at an accelerating pace. Both results are well matched by a model assuming that fitness cost increases exponentially (approximately second order) with the number of accumulated mutations. This result is consistent with fitness loss due primarily to interactions between specific pairs of gene products. I also estimate that ~25% of the mutations with detectable fitness effects were beneficial. This result can be explained by the fact that the effects of most mutations are small relative to the distance from a local fitness optimum.


THERE is substantial interest in understanding the cumulative fitness costs of mutations that individually have small effects. Possible interactions have important implications for a number of issues in evolutionary biology (WAGNER et al. 1998; WOLF et al. 2000; KOUYOS et al. 2007). For example, synergistic epistasis (cumulative fitness loss greater than predicted from independent effects) could help account for the maintenance of sexual reproduction and recombination (KONDRASHOV 1988; WEST et al. 1998; DE VISSER and ELENA 2007). Even with asexual reproduction, synergistic interactions might provide an effective limit to Muller's ratchet by disproportionately eliminating those genomes with the largest number of mutations (KONDRASHOV 1994) . Epistasis also has implications for the survival of small populations and is believed to be an important factor in speciation (ORR 1995; FENSTER et al. 1997; WAGNER et al. 1998; WILLETT 2006). Despite the importance of these issues, relevant experimental results are limited, inconsistent, and controversial (KIBOTA and LYNCH 1996; SHABALINA et al. 1997; KONDRASHOV 1998; FRY et al. 1999; GARCIA-DORADO et al. 1999; KEIGHTLEY and EYRE-WALKER 2000; FRY 2004; KOUYOS et al. 2007). Many attempts to detect synergistic epistasis depend on indirect methods, are sensitive to assumptions in the underlying models used to evaluate the data, or are too labor intensive to allow sample sizes that can give unambiguous results (DENG and FU 1998). We need additional experiments with the power and sensitivity to detect directly the fitness consequences of interactions among sets of spontaneous mutations.

Assessment of fitness effects in lines that have accumulated multiple spontaneous mutations is time and labor intensive but is the most direct approach (MUKAI 1969). Alternative hypotheses lead to straightforward testable predictions (KOUYOS et al. 2007). Independent action, the most useful null hypothesis, should produce multiplicative effects in combination. For multiple mutations with small effects that have not been measured individually, a plot of log fitness against the number of mutations (or the time of accumulation) should yield a straight line. Deviations from linearity indicate epistatic interaction and provide information about the nature of those interactions; negative curvature indicates synergistic epistasis while positive curvature reveals antagonistic or "diminishing returns" interactions. I report here the results of a mutation accumulation experiment with yeast using sensitive competition experiments (THATCHER et al. 1998) to directly measure realized fitness relative to a marked standard in a defined environment.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Yeast stocks:
I accumulated mutations in 48 independent lines (MA lines), eight in each of six series, A–F. Lines within each series were initiated from the same cloned parent. Those six progenitors were, in turn, derived from a single laboratory stock, FY10 (leu2{Delta}1, ura3-52, MAT{alpha}). Random insertions into FY10 of LacZ in a construct that also carried LEU2 as a selectable marker were recovered and tested for fitness effects in previous experiments (THATCHER et al. 1998). Preliminary experiments showed that leu2 was not selectively neutral even in rich medium, so leu2 in FY10 was repaired by homologous recombination to produce FY10(LEU+). For these experiments I chose one of the LacZ insertions (TD50) that, within the sensitivity of our methods, is selectively neutral in comparison to FY10(LEU+). Lines A, C, and E were derived from FY10(LEU+) while B, D, and F are from TD50 (carrying the LacZ marker used to score competition experiments). A and B were converted to MATa (HERSKOWITZ and JENSEN 1991), C and D are MAT{alpha}, and E and F are diploids homozygous at all loci other than MAT, produced by mating A to C and B to D, respectively. The manipulations at the MAT locus were intended to facilitate future genetic analyses of the MA lines and are not directly relevant to this experiment. In competition experiments to measure fitness, the parental A line serves as a competitor for all B accumulation lines while the B parent serves reciprocally for the A lines. Similarly, C pairs with D and E with F. Reference stocks of the parental lines were stored frozen and periodic control experiments confirmed that the matched parents continue to be selectively neutral with respect to each other. Thus, each fitness determination is, in essence, a direct comparison to zero time fitness.

Mutation accumulation:
The 48 MA lines were passed through single-cell bottlenecks on a 3- to 4-day cycle for 200 cycles. Each cycle was initiated by picking a random colony of each line and streaking on rich medium (YPD). Random selection was assured by overlaying each plate on a template with a marked target and selecting the colony closest to that target. Each cycle represents ~24 generations (growing to a colony of ~2 x 107 cells). Because the single-cell bottlenecks predominate, the effective population size is ~10 (HARTL and CLARK 1997). Thus, deleterious mutations with a fitness effect < ~10% should be "invisible" to selection and accumulate at a rate that simply reflects mutation rates. Conversely, beneficial mutations of similar magnitude will not be favored by selection.

Preliminary experiments showed that petite mutations, almost all defective in the mitochondrial genome, are fixed in ~5% of the lines in each cycle. These have relatively large fitness effects that would mask the mutations of small effect that I wanted to study, so I eliminated petites as follows: In addition to streaking on YPD, each selected colony was patched to medium with glycerol as the major carbon source (YPG), on which petites do not grow. When the next cycle was initiated, the YPG patch plate was inspected to identify petites, and the corresponding streaks were discarded. To maintain those lines, I returned to the plates from the previous cycle (stored refrigerated for that purpose) and selected a new random colony. The appearance of petites was a minor annoyance, but provides direct confirmation that the bottlenecking procedure was working as expected to allow fixation even of mutations with rather significant deleterious effects. All stocks also were checked regularly for LacZ and MAT phenotypes, and the rare deviations (mostly loss of LacZ in B and D lines and switching of diploids to MATa or MAT{alpha}) were corrected by returning to the most recent available refrigerated or frozen sample with the correct phenotype. The combination of eliminating petites and, to a lesser extent, other phenotypic deviations means that actual elapsed generations is, on average, ~6% less than the nominal value.

Fitness measurements:
My method for determining relative fitness is slightly modified from that previously described (THATCHER et al. 1998). Briefly, an overnight culture of each MA line was mixed 1:1 with a similar culture of the matched parent differing at the LacZ marker and then immediately diluted 1:256 into fresh medium. Overnight cultures go to saturation, so the relative frequencies of the two genotypes in the initial mixture is close to 1:1. The mixed cultures were maintained at 30° in 1 ml of rich medium (YPD) in 16- x 100-mm tubes angled at 25° from the horizontal and rotated at 40 rpm. They were back diluted 1:256 daily. Under these conditions, cultures return to saturation each day, completing eight generations per dilution cycle. Appropriately diluted samples were plated periodically and scored for frequency of LacZ colonies. A selection coefficient is estimated from the slope of a natural log plot of the ratio of the two colony types as a function of generations of competition (LENSKI 1991). This procedure detects fitness differences as small as ~0.5% (see RESULTS).

LacZ detection:
Yeast colonies positive for LacZ do not stain readily on X-gal plates so the following procedure was adopted. Test plates contained 100 mg/liter of X-gal in YPD with 2% agar. Samples from the competition cultures were plated at a dilution calculated to give a few hundred colonies per plate and incubated at 30° for 3 days to give easily visible colonies. Plates were then exposed to chloroform vapor for 30 min and incubated at 37° to develop the blue color. The chloroform kills enough of the cells to release easily detected amounts of β-galactosidase. Blue and white colonies were hand counted.

Adaptation lines:
To estimate the fitness of the parental lines relative to a local fitness optimum, a series of 24 lines (4 each from parents A–F) was allowed to adapt for 160 cycles (1280 generations) to the conditions used for the competition experiments. Effective population sizes were very large, allowing effective selection. Improvements in fitness were then estimated from standard competition experiments in comparison to the matched frozen parental line. The average improvement across all lines was 8% and the maximum was 12%.

An exponential model of synergistic interactions:
Results of the accumulation experiments are evaluated in comparison to a simple model that assumes a linear increase in the average number of mutations, a Poisson distribution of mutation frequency, and a fitness cost that is an exponential function of the number of mutations. The following five-step algorithm was executed in a standard spreadsheet program.

  1. Generate hypothetical sets of MA lines for each accumulation cycle examined in the real experiment, each set having a Poisson distribution around a given mean number of new mutations. The mean is determined from a specified genomewide mutation rate (U) and was varied across trials.
  2. Assign a fitness cost to each MA line proportional to an exponential function of the number of mutations. The exponent was varied across trials.
  3. Compute the mean fitness loss and the variance for each set of MA lines.
  4. Normalize all results to the observed mean and variance at the end of the real experiment (cycle 200).
  5. Convert mean fitness loss to log mean fitness and plot that value and, separately, the fitness variance as a function of cycle number.

Once templates for all steps were generated, a series of trials with different values of input variables explored the response across a plausible parameter space. Thus, two variables are considered (genomewide mutation rate and the exponent assigned to fitness costs) and it is the shapes of the curves that are of interest.


    RESULTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Fitness histories of all 48 MA lines are presented in Figure 1 together with changes in the mean. Three notable features are evident.

  1. There is a marked tendency toward accelerating fitness loss in later generations. For example, the mean change per generation in the third scored time interval is more than five times that in the first interval, and six of the nine most precipitous losses in fitness occur in that interval (with one in the first interval and two in the second).
  2. Lines diverge dramatically from one another and also apparently at an accelerating pace.
  3. Despite the general trend toward decreased fitness and the presumed absence of effective selection, interspersed increases in fitness (individual line segments with positive slopes) are fairly common at all stages of the experiment. I will return to this last point in Beneficial mutations.


Figure 1
View larger version (23K):
[in this window]
[in a new window]
[Download PPT slide]
 
FIGURE 1.— Fitness changes in 48 MA lines. The dots are individual fitness measurements and the lines connect successive measurements on the same line. The open circles and heavier line represent the means. Each accumulation cycle corresponds to ~24 generations, so the data cover ~4800 generations of mutation accumulation under relaxed selection. Cycle zero assays are on replicates of the starting stocks, so the spread among those points reflects the precision of the competition assays. The markedly greater spreads at later times reflect cumulative fitness differences. For stylistic consistency, zero time points are connected arbitrarily (by numerical order of the replicate samples) to later line histories. Since many lines remain clustered near fitness 1.0 throughout the experiment, the scale is chosen to emphasize differences in the range 0.9–1.04. As a result, some points drop off the bottom of the graph, but their positions can be inferred from the slopes of the lines intersecting the 0.9-fitness threshold.

 
With the fitness of the parental lines defined as 1.0, the mean fitness at cycle 45 is 0.996 ± 0.00014. The reported uncertainty is the 95% confidence interval based on standard error of the mean. Even though the fitness distribution across strains does not appear to be normal, it is likely that estimates of means based on a large number of measurements will be so distributed (BAILEY 1995). The change by cycle 45 is not quite significant (P {cong} 0.1). At cycles 130 and 200, the corresponding values are 0.981 ± 0.00045 and 0.948 ± 0.00093, both highly significant (P < 0.01 and P < 0.0005, respectively). Even by cycle 45 there was a 100-fold increase in the variance relative to zero time controls (a change in the standard deviation ~3.7 times the change in the mean), suggesting that there have been real changes but offsetting gains and losses. This is relevant to the evidence for frequent beneficial mutations discussed below. I also evaluated separately the haploid and diploid lines. Surprisingly, the quantitative differences are modest and not significant. For example, at cycle 200, the haploid and diploid means are 0.941 ± 0.0011 and 0.962 ± 0.0017 and the difference is not significant (P {cong} 0.5). It is of interest in this context that many of the initially haploid lines apparently accumulate duplications to approach a diploid chromosome number (M. LYNCH, personal communication). This change probably is not driven by immediate selective advantage since standard competition experiments with parental haploids against matched diploids (A or C against F and B or D against E) reveal no significant differences. In any case, the qualitative trends (negative curvature of a log mean fitness plot and accelerating increase in the variance) that are the focus of the following analysis are the same for nominally haploid and diploid lines, so all are treated together.

Comparisons to an exponential model:
Changes in log mean fitness and in the fitness variance are summarized in Figure 2 in comparison to results of the exponential model of fitness costs described above. The parameter space evaluated in the model covers U over a 30-fold range from 0.00037 through 0.011 and cost exponents from 1 through 4. The range of mutation rates was based on several reasonably careful of estimates of U in yeast (DRAKE 1991; WLOCH et al. 2001b; ZEYL and DEVISSER 2001), ranging from 0.00024 to 0.0031. Perhaps the most direct determination is by WLOCH et al. (2001b) who scored phenotypes segregating in a large number of tetrads. Their estimate is 0.0011. That is almost certainly low for my purposes, perhaps substantially so, since only mutations with detectable phenotypes were counted. Cryptic mutations that could have significant synergistic effects may be much more common than mutations with detectable effects (DAVIES et al. 1999). Thus, I took the above estimates as defining the low end of the range that should be evaluated and went as high as 10-fold above the value given by WLOCH et al. (2001b).


Figure 2
View larger version (16K):
[in this window]
[in a new window]
[Download PPT slide]
 
FIGURE 2.— Changes in log average fitness and in fitness variance compared to predictions from an exponential model. In A–D, the experimental data are shown as solid circles connected with solid lines. Epistasis is revealed as deviations from linearity in the plot of log mean fitness (A and C). Changes in fitness variance (B and D) reflect divergence among lines. In A and B, the experimental data are compared to a model (see text) with fitness costs assumed to be first-order (squares), second-order (open circles) or third-order (triangles) functions of mutation number. The genomewide mutation rate is held constant at 0.0033. In C and D, the same data are compared to only the second-order model but with the genomewide mutation rate set at 0.00037 (squares) or 0.011 (open circles).

 
Theoretical treatments (and simple probability arguments) predict that the exponent assigned to the fitness cost function of synergistic mutations should equal the number of components directly involved in each interaction (LYNCH 1994; ORR 1995). Thus, for example, pairwise interactions should produce a fitness cost that is a second-order function of mutation number (because the probability that both targets in a given pair have been mutated will increase as the square of the total number of hits). The expectation that independent fitness effects should lead to a linear plot of log mean fitness corresponds, of course, to the first-order case. Consistent with synergistic epistasis, the plot of log mean fitness as a function of accumulation cycles has a markedly negative curvature (Figure 2). Specifically, both the accelerating fitness loss and the rapid increase in variance in the experimental results are well matched by the model assuming second-order costs (Figure 2, A and B). A fourth-order trial (not shown) diverges further in the direction established by the third order. This pattern holds across the range of mutation rates examined (Figure 2, C and D). In fact, the model is not very sensitive to mutation rate, but converges toward the data at the higher end of the tested range.

Beneficial mutations:
It is generally believed that the vast majority of new mutations are deleterious, so I examined the frequency and significance of the apparent improvements seen in Figure 1 in more detail. I assume that the tight clustering of the zero time controls (parent vs. parent competitions) is an indication of the precision and reproducibility of the fitness assays. It seems reasonable that variation among such replicates approximates a normal distribution (even if the effects of mutations do not) and it is experimental error that is at issue in identifying "real" changes, so I use the standard deviation of the control experiments (±0.0024) as the point of reference for estimating the probability that a given change is significant. In the context of estimating the fraction of changes that are improvements, it is not entirely clear whether one should worry more about rejecting real changes (e.g., with stringent criteria that incorporate correction for multiple comparisons) or accepting spurious ones. Fortunately, the results in this case are not very sensitive to the chosen criterion. At P < 0.05 (with no correction for multiple comparisons), we "accept" 65 changes throughout the experiment, 16 of which (24.6%) are improvements. Similarly at P < 0.01, 14 of 57 changes (24.6%) are improvements, as are 11 of 48 (22.9%) at P < 0.001 and 10 of 45 (22.2%) at P < 0.00035. The last threshold corresponds to the Bonferroni correction for P < 0.05 for 144 comparisons (48 lines over three intervals). The largest positive changes are ~2% by cycle 45 and ~4% for later comparisons, well below the level at which selection should be effective (~10%), so I conclude that almost a quarter of the significant fitness changes fixed by drift were improvements.

Independence of late fitness changes:
One other feature of the line histories is less evident by inspection; early fitness loss does not predispose to increased fitness loss in later accumulation cycles. Kendall's rank correlation coefficient was used to test for correlation between fitness at each intermediate cycle and the change in fitness over each subsequent interval. None of three comparisons (w45:{Delta}w45->130, w45:{Delta}w45->200 and w130:{Delta}w130->200) was significant (P > 0.2 in all cases), and the tendency in the second comparison was even of reversed sign. A standard correlation analysis of the quantitative data gives substantially similar results. One comparison (w45::{Delta}w45->130) has a correlation coefficient of ~0.6, which is significant, given the large number of data points. However, that result depends entirely on one data point (stock D1, which had the largest loss by cycle 45). With that line dropped from the analysis, the magnitude of the correlation is no longer significant and the sign is reversed. Thus, in general, present fitness is a poor predictor of vulnerability to subsequent fitness loss. One explanation, examined below, is that cumulative fitness is affected by specific interactions that have little to do with the independent fitness effects of individual mutations.


    DISCUSSION
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Second-order fitness costs:
On our central question, the fitness costs of mutations that individually have small (or negligable) effects appear to be synergistic, at least for the system and conditions examined here. Specifically, the data match closely a model that assumes that costs are approximately a second-order function of the number of mutations. A straightforward explanation is that fitness effects of an important class of mutations depend on specific interactions between pairs of gene products. Higher-order effects (interactions between multiple sites) are certainly possible but second-order effects might predominate because many interactions are inherently pairwise or because, by simple probability arguments, appropriately mutated multi-member sets are rare.

What interactions between pairs of gene products are relevant? Mutations in a single pathway might interact (RICE 1998), but according to metabolic control theory (KACSER and BURNS 1981; SZATHMARY 1993) these effects should be antagonistic. Synergistic interactions seem more likely for molecules that make contact with and recognize specific partners. At least 50% of yeast proteins seem to be members of direct-interaction networks (HO et al. 2002; GAVIN et al. 2006). Examples include components of the same molecular complex, cis- and trans-acting factors involved in regulation of the same gene, and partners in a signal transduction event. Proteins in large complexes may interact with multiple partners (KIM et al. 2006) but typically do so through different surfaces or in a mutually exclusive way (hence, still pairwise). In all of these cases, alterations in one component alone may have little effect but may be incompatible with specific alterations in a partner. Second-site noncomplementation, used in genetic screens to identify interaction partners, depends on just such specific incompatibilities. (STEARNS and BOTSTEIN 1988; HAWLEY and WALKER 2003; KOMILI and ROTH 2007).

Mechanisms that enhance robustness in gene networks, whether evolved under stabilizing selection (WAGNER 2000) or emerging as intrinsic properties of complex networks (SIEGAL and BERGMAN 2002), also involve pairwise interactions. Of two principal mechanisms identified by WAGNER (2000), interactions between genes of unrelated function seem more prevalent than duplicate or overlapping functions, but either seems compatible with mutations in specific pairs of sites being more deleterious than single hits. Also, the lack of correlation between early and late fitness changes is entirely consistent with specific pairwise interactions since the effect of a given mutation is contingent on others present in the same background.

Results of a few other MA experiments also suggest synergistic interactions (MUKAI 1969; FRY 2004; AVILA et al. 2006). In the most recent of those, new Drosophila MA lines were established from end-products of an earlier experiment; declines in viability were substantially accelerated relative to that previous experiment. The authors interpret that result in terms of an elevated mutation rate in the later accumulation cycles. That is certainly a formal possibility (also for my data), but synergistic interaction seems a more parsimonious explanation; the mutation rate hypothesis requires that many lines independently acquire mutators.

Other studies have failed to detect synergistic epistatsis or have reported variable results (KOUYOS et al. 2007). However, the designs of many experiments are not well suited to detection of specific pairwise interactions. For example, combining a limited number of known mutations (DE VISSER et al. 1997; ELENA and LENSKI 1997; WHITLOCK and BOURGUET 2000; SZAFRANIEC et al. 2003; SANJUAN and ELENA 2006) should reveal epistasis only sporadically (when appropriate pairs happen to be included). In addition, knockout mutations may not be suitable since there is no opportunity for interaction with a partner (YOOK et al. 2001). In fact, knockouts of two essential components in the same complex might show antagonistic epistasis since either alone could inactivate the complex. Some studies have used insertional mutations that are likely to block the formation of any functional gene product (ELENA and LENSKI 1997, 2001). Others have used a mutator to accelerate the process in MA experients (WLOCH et al. 2001A), but those mutators are known to produce a large proportion of frameshifts that, again, most often will be knockouts (KORONA 1999).

Beneficial mutations and "Fisher's microscope":
It is widely believed that the vast majority of new mutations are harmful (CHARLESWORTH and CHARLESWORTH 1998; KEIGHTLEY and LYNCH 2003). However, the available data compel only the conclusion that the average effect is deleterious (SHAW et al. 2003), and beneficial mutations have been reported to be a significant fraction of total detected mutations in a few MA studies broadly similar to this one. For example, the beneficial fraction is ~0.05 in another study with yeast (JOSEPH and HALL 2004) and almost half in Arabidopsis (SHAW et al. 2000). Other studies have detected small spontaneous improvements or an increase in variance without significant decrease in the mean without explicitly analyzing those features (KEIGHTLEY and CABALLERO 1997; ZEYL and DEVISSER 2001). Interestingly, selection experiments with Escherichia coli detect beneficial mutations of small effect at an elevated frequency when effective population sizes are small (PERFEITO et al. 2007). At Ne = 2 x 104, ~10% of mutations with detectable fitness effects were beneficial (~1000-fold higher than obtained with Ne = 107), and the mean selection coefficient was only 0.013.

R. A. Fisher famously used a microscope analogy to argue that adaptation depends primarily on mutations of small effect (FISHER 1930). Central to his argument was the expectation that the probability of improvement increases as the magnitude of the effect decreases, approaching 50% as a limit. It follows that if most mutations have very small effects (THATCHER et al. 1998), almost half of all mutations will be improvements. How small must the effects be to see frequent improvement? The key is magnitude relative to distance (in phenotypic space) from a nearby optimum (ORR 2001). For variation in one dimension, Fisher's geometrical model gives the frequency of beneficial mutations as Formula where r is the magnitude of the change and d is the distance from an optimum. Mutations affect phenotypes in a multidimensional space, but measuring distance from a local fitness optimum in a defined environment essentially collapses all of those effects to one dimension. When allowed to adapt at large effective population sizes to the exact conditions of the fitness assays (the defined environment), the lines used in this experiment improved in fitness by an avareage of 0.08 and by a maximum of 0.12. The latter value presumably represents a better estimate of the distance from a local optimum. The average selection coefficient for 65 significant changes (P < 0.05) observed in my accumulation lines is 0.055. If we use d = 0.012 and r = 0.0055 in Fisher's expression, the expected frequency of improvements is 0.27—in remarkable agreement with the 0.22–0.25 actually observed.

There have, of course, been refinements and alternatives to Fisher's model in the intervening years (WAXMAN and WELCH 2005; MARTIN and LENORMAND 2006; ORR 2006), but these do not fundamentally challenge the expectation that very small changes have a reasonable chance of being improvements. The present results provide empirical support for that insight. Direct physical interactions between specific partners, as suggested for synergistic fitness loss, also may be relevant to stochastic improvements; second-site suppression of a mutant phenotype is another well-known phenomenon used in genetic screens to identify interaction partners (HAWLEY and WALKER 2003). It is to be expected that the frequency of spontaneous improvement will vary from experiment to experiment since that frequency should depend on how far the starting stocks were from an optimum for the conditions used to measure fitness. Indeed, this leads to a testable prediction; stochastic improvements should largely disappear if strains are adapted to the test environment for many generations at large effective population sizes prior to use in MA experiments.

So the idea that most mutations are harmful might more correctly be applied to the relatively small class of mutations with fitness effects large enough to be detected by conventional methods. And there is no fundamental conflict between the finding that beneficial mutations are relatively common in specific MA experiments and molecular evidence, suggesting that most nonsynonomous mutations are eliminated by purifying selection. The former result relates to specific stocks that may, by chance, sit some distance from the optimum for the chosen experimental conditions; the latter presumably integrates results over the fluctuating environments experienced by natural populations. The present results nevertheless are relevant to natural populations; they suggest that a population displaced from an adaptive optimum or shifting to a new one will not have long to wait for mutations that are beneficial in the new environment. This fits nicely with recent evidence for frequent adaptive substitutions in both coding and noncoding DNA sequences (SMITH and EYRE-WALKER 2002; ANDOLFATTO 2005).


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
I thank my colleagues in the Department of Biology, particularly John S. Parkinson, for permitting me to keep research space following retirement and for providing funds (in lieu of salary for an extra teaching assignment) used to support this research. Jon Seger, Fred Adler, Mike Lynch, and Bernardo Lemos provided encouragement and useful suggestions, and comments from two anonymous reviewers prompted significant improvements.


    LITERATURE CITED
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 

ANDOLFATTO, P., 2005 Adaptive evolution of non-coding DNA in Drosophila. Nature 437: 1149–1152.[CrossRef][Medline]

AVILA, V., D. CHAVARRIAS, E. SANCHEZ, A. MANRIQUE, C. LOPEZ-FANJUL et al., 2006 Increase of the spontaneous mutation rate in a long-term experiment with Drosophila melanogaster. Genetics 173: 267–277.[Abstract/Free Full Text]

BAILEY, T. J., 1995 Statistical Methods in Biology. Cambridge University Press, Cambridge, UK.

CHARLESWORTH, B., and D. CHARLESWORTH, 1998 Some evolutionary consequences of deleterious mutations. Genetica 103: 3–19.[CrossRef]

DAVIES, E. K., A. D. PETERS and P. D. KEIGHTLEY, 1999 High frequency of cryptic deleterious mutations in Caenorhabditis elegans. Science 285: 1748–1751.[Abstract/Free Full Text]

DENG, H. W., and Y. X. FU, 1998 On the three methods for estimating deleterious genomic mutation parameters. Genet. Res. 71: 223–236.[CrossRef][Medline]

DE VISSER, J. A., and S. F. ELENA, 2007 The evolution of sex: empirical insights into the roles of epistasis and drift. Nat. Rev. Genet. 8: 139–149.[Medline]

DE VISSER, J. A., R. F. HOEKSTRA and H. VAN DEN ENDE, 1997 An experimental test for synergistic epistasis and its application in Chlamydomonas. Genetics 145: 815–819.[Abstract]

DRAKE, J. W., 1991 A constant rate of spontaneous mutation in DNA-based microbes. Proc. Natl. Acad. Sci. USA 88: 7160–7164.[Abstract/Free Full Text]

ELENA, S. F., and R. E. LENSKI, 1997 Test of synergistic interactions among deleterious mutations in bacteria. Nature 390: 395–398.[CrossRef][Medline]

ELENA, S. F., and R. E. LENSKI, 2001 Epistasis between new mutations and genetic background and a test of genetic canalization. Evolution 55: 1746–1752.[CrossRef][Medline]

FENSTER, C. B., L. F. GALLOWAY and L. CHAO, 1997 Epistasis and its consequences for the evolution of natural populations. Trends Ecol. Evol. 12: 282–286.[CrossRef]

FISHER, R. A., 1930 The Genetical Theory of Natural Selection. Clarendon Press, Oxford.

FRY, J. D., 2004 On the rate and linearity of viability declines in Drosophila mutation-accumulation experiments: genomic mutation rates and synergistic epistasis revisited. Genetics 166: 797–806.[Abstract/Free Full Text]

FRY, J. D., P. D. KEIGHTLEY, S. L. HEINSOHN and S. V. NUZHDIN, 1999 New estimates of the rates and effects of mildly deleterious mutation in Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 96: 574–579.[Abstract/Free Full Text]

GARCIA-DORADO, A., C. LOPEZ-FANJUL and A. CABALLERO, 1999 Properties of spontaneous mutations affecting quantitative traits. Genet. Res. 74: 341–350.[CrossRef][Medline]

GAVIN, A. C., P. ALOY, P. GRANDI, R. KRAUSE, M. BOESCHE et al., 2006 Proteome survey reveals modularity of the yeast cell machinery. Nature 440: 631–636.[CrossRef][Medline]

HARTL, D. L., and A. G. CLARK, 1997 Principles of Population Genetics. Sinauer Associates, Sunderland, MA.

HAWLEY, R. S., and M. WALKER, 2003 Advanced Genetic Analysis: Finding Meaning in a Genome. Blackwell, Oxford.

HERSKOWITZ, I., and R. E. JENSEN, 1991 Putting the HO gene to work: practical uses for mating-type switching. Methods Enzymol. 194: 132–146.[Medline]

HO, Y., A. GRUHLER, A. HEILBUT, G. D. BADER, L. MOORE et al., 2002 Systematic identification of protein complexes in Saccharomyces cerevisiae by mass spectrometry. Nature 415: 180–183.[CrossRef][Medline]

JOSEPH, S. B., and D. W. HALL, 2004 Spontaneous mutations in diploid Saccharomyces cerevisiae: more beneficial than expected. Genetics 168: 1817–1825.[Abstract/Free Full Text]

KACSER, H., and J. A. BURNS, 1981 The molecular basis of dominance. Genetics 97: 639–666.[Abstract/Free Full Text]

KEIGHTLEY, P. D., and A. CABALLERO, 1997 Genomic mutation rates for lifetime reproductive output and lifespan in Caenorhabditis elegans. Proc. Natl. Acad. Sci. USA 94: 3823–3827.[Abstract/Free Full Text]

KEIGHTLEY, P. D., and A. EYRE-WALKER, 2000 Deleterious mutations and the evolution of sex. Science 290: 331–333.[Abstract/Free Full Text]

KEIGHTLEY, P. D., and M. LYNCH, 2003 Toward a realistic model of mutations affecting fitness. Evolution 57: 683–685.[CrossRef][Medline]

KIBOTA, T. T., and M. LYNCH, 1996 Estimate of the genomic mutation rate deleterious to overall fitness in E. coli. Nature 381: 694–696.[CrossRef][Medline]

KIM, P. M., L. J. LU, Y. XIA and M. B. GERSTEIN, 2006 Relating three-dimensional structures to protein networks provides evolutionary insights. Science 314: 1938–1941.[Abstract/Free Full Text]

KOMILI, S., and F. P. ROTH, 2007 Genetic interaction screens advance in reverse. Genes Dev. 21: 137–142.[Free Full Text]

KONDRASHOV, A. S., 1988 Deleterious mutations and the evolution of sexual reproduction. Nature 336: 435–440.[CrossRef][Medline]

KONDRASHOV, A. S., 1994 Muller's ratchet under epistatic selection. Genetics 136: 1469–1473.[Abstract]

KONDRASHOV, A. S., 1998 Measuring spontaneous deleterious mutation process. Genetica 103: 183–197.[CrossRef]

KORONA, R., 1999 Unpredictable fitness transitions between haploid and diploid strains of the genetically loaded yeast Saccharomyces cerevisiae. Genetics 151: 77–85.[Abstract/Free Full Text]

KOUYOS, R. D., O. K. SILANDER and S. BONHOEFFER, 2007 Epistasis between deleterious mutations and the evolution of recombination. Trends Ecol. Evol. 22: 308–315.[CrossRef][Medline]

LENSKI, R. E., 1991 Quantifying fitness and gene stability in microorganisms. Biotechnology 15: 173–192.[Medline]

LYNCH, M., 1994 Neutral models of phenotypic evolution, pp. 86–108 in Ecological Genetics, edited by L. A. REAL. Princeton University Press, Princeton, NJ.

MARTIN, G., and T. LENORMAND, 2006 A general multivariate extension of Fisher's geometrical model and the distribution of mutation fitness effects across species. Evolution 60: 893–907.[CrossRef][Medline]

MUKAI, T., 1969 The genetic structure of natural populations of Drosophila melanogaster. VII. Synergistic interactions of spontaneous mutant polygenes controlling viability. Genetics 61: 749–761.[Free Full Text]

ORR, H. A., 1995 The population genetics of speciation: the evolution of hybrid incompatibilities. Genetics 139: 1805–1813.[Abstract]

ORR, H. A., 2001 The "sizes" of mutations fixed in phenotypic evolution: a response to Clarke and Arthur. Evol. Dev. 3: 121–123.[CrossRef][Medline]

ORR, H. A., 2006 The distribution of fitness effects among beneficial mutations in Fisher's geometric model of adaptation. J. Theor. Biol. 238: 279–285.[Medline]

PERFEITO, L., L. FERNANDES, C. MOTA and I. GORDO, 2007 Adaptive mutations in bacteria: high rate and small effects. Science 317: 813–815.[Abstract/Free Full Text]

RICE, W. R., 1998 Requisite mutational load, pathway epistasis and deterministic mutation accumulation in sexual versus asexual populations. Genetica 103: 71–81.[CrossRef]

SANJUAN, R., and S. F. ELENA, 2006 Epistasis correlates to genomic complexity. Proc. Natl. Acad. Sci. USA 103: 14402–14405.[Abstract/Free Full Text]

SHABALINA, S. A., L. YAMPOLSKY and A. S. KONDRASHOV, 1997 Rapid decline of fitness in panmictic populations of Drosophila melanogaster maintained under relaxed natural selection. Proc. Natl. Acad. Sci. USA 94: 13034–13039.[Abstract/Free Full Text]

SHAW, R. G., D. L. BYERS and E. DARMO, 2000 Spontaneous mutational effects on reproductive traits of Arabidopsis thaliana. Genetics 155: 369–378.[Abstract/Free Full Text]

SHAW, R., F. SHAW and C. GEYER, 2003 What fraction of mutations reduces fitness? A reply to Keightley and Lynch. Evolution 57: 686–689.[CrossRef]

SIEGAL, M. L., and A. BERGMAN, 2002 Waddington's canalization revisited: developmental stability and evolution. Proc. Natl. Acad. Sci. USA 99: 10528–10532.[Abstract/Free Full Text]

SMITH, N. G., and A. EYRE-WALKER, 2002 Adaptive protein evolution in Drosophila. Nature 415: 1022–1024.[CrossRef][Medline]

STEARNS, T., and D. BOTSTEIN, 1988 Unlinked noncomplementation: isolation of new conditional-lethal mutations in each of the tubulin genes of Saccharomyces cerevisiae. Genetics 119: 249–260.[Abstract/Free Full Text]

SZAFRANIEC, K., D. M. WLOCH, P. SLIWA, R. H. BORTS and R. KORONA, 2003 Small fitness effects and weak genetic interactions between deleterious mutations in heterozygous loci of the yeast Saccharomyces cerevisiae. Genet. Res. 82: 19–31.[CrossRef][Medline]

SZATHMARY, E., 1993 Do deleterious mutations act synergistically? Metabolic control theory provides a partial answer. Genetics 133: 127–132.[Abstract]

THATCHER, J. W., J. M. SHAW and W. J. DICKINSON, 1998 Marginal fitness contributions of nonessential genes in yeast. Proc. Natl. Acad. Sci. USA 95: 253–257.[Abstract/Free Full Text]

WAGNER, A., 2000 Robustness against mutations in genetic networks of yeast. Nat. Genet. 24: 355–361.[CrossRef][Medline]

WAGNER, G. P., M. D. LAUBICHLER and H. BAGHERI-CHAICHIAN, 1998 Genetic measurement of theory of epistatic effects. Genetica 102–103: 569–580.[CrossRef][Medline]

WAXMAN, D., and J. J. WELCH, 2005 Fisher's microscope and Haldane's ellipse. Am. Nat. 166: 447–457.[CrossRef][Medline]

WEST, S. A., A. D. PETERS and N. H. BARTON, 1998 Testing for epistasis between deleterious mutations. Genetics 149: 435–444.[Abstract/Free Full Text]

WHITLOCK, M. C., and D. BOURGUET, 2000 Factors affecting the genetic load in Drosophila: synergistic epistasis and correlations among fitness components. Evolution 54: 1654–1660.[CrossRef][Medline]

WILLETT, C. S., 2006 Deleterious epistatic interactions between electron transport system protein-coding loci in the copepod Tigriopus californicus. Genetics 173: 1465–1477.[Abstract/Free Full Text]

WLOCH, D., R. BORTS and R. KORONA, 2001a Epistatic interactions of spontaneous mutations in haploid strains of the yeast Saccharomyces cerevisiae. J. Evol. Biol. 14: 310–316.[CrossRef]

WLOCH, D. M., K. SZAFRANIEC, R. H. BORTS and R. KORONA, 2001b Direct estimate of the mutation rate and the distribution of fitness effects in the yeast Saccharomyces cerevisiae. Genetics 159: 441–452.[Abstract/Free Full Text]

WOLF, W. B., E. D. I. BRODIE and M. J. WADE, 2000 Epistasis and the Evolutionary Process. Oxford University Press, Oxford.

YOOK, K. J., S. R. PROULX and E. M. JORGENSEN, 2001 Rules of nonallelic noncomplementation at the synapse in Caenorhabditis elegans. Genetics 158: 209–220.[Abstract/Free Full Text]

ZEYL, C., and J. A. DEVISSER, 2001 Estimates of the rate and distribution of fitness effects of spontaneous mutation in Saccharomyces cerevisiae. Genetics 157: 53–61.[Abstract/Free Full Text]

Communicating editor: J. LAWRENCE


Related articles in Genetics:

ISSUE HIGHLIGHTS

Genetics 2008 178: NP. [Full Text]  




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
genetics.107.080853v1
178/3/1571    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in Genetics
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Google Scholar
Right arrow Articles by Dickinson, W. J.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Dickinson, W. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2008 by the Genetics Society of America.