Reasearch Awards nomination

Email updates

Keep up to date with the latest news and content from Genome Biology and BioMed Central.

Open Access Research

The strength of genetic interactions scales weakly with mutational effects

Andrea Velenich* and Jeff Gore

Author Affiliations

Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

For all author emails, please log on.

Genome Biology 2013, 14:R76  doi:10.1186/gb-2013-14-7-r76


The electronic version of this article is the complete one and can be found online at: http://genomebiology.com/2013/14/7/R76


Received:11 February 2013
Revisions received:1 July 2013
Accepted:26 July 2013
Published:26 July 2013

© 2013 Velenich et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Genetic interactions pervade every aspect of biology, from evolutionary theory, where they determine the accessibility of evolutionary paths, to medicine, where they can contribute to complex genetic diseases. Until very recently, studies on epistatic interactions have been based on a handful of mutations, providing at best anecdotal evidence about the frequency and the typical strength of genetic interactions. In this study, we analyze a publicly available dataset that contains the growth rates of over five million double knockout mutants of the yeast Saccharomyces cerevisiae.

Results

We discuss a geometric definition of epistasis that reveals a simple and surprisingly weak scaling law for the characteristic strength of genetic interactions as a function of the effects of the mutations being combined. We then utilized this scaling to quantify the roughness of naturally occurring fitness landscapes. Finally, we show how the observed roughness differs from what is predicted by Fisher's geometric model of epistasis, and discuss the consequences for evolutionary dynamics.

Conclusions

Although epistatic interactions between specific genes remain largely unpredictable, the statistical properties of an ensemble of interactions can display conspicuous regularities and be described by simple mathematical laws. By exploiting the amount of data produced by modern high-throughput techniques, it is now possible to thoroughly test the predictions of theoretical models of genetic interactions and to build informed computational models of evolution on realistic fitness landscapes.

Keywords:
Epistasis; Evolution; Fitness landscapes; Genetic interactions; Yeast

Background

Genetic interactions [1] have shaped the evolutionary history of life on earth. They have been found to limit the accessibility of evolutionary paths [2], to confine populations to suboptimal evolutionary states and, on larger time scales, to control the rate of speciation [3]. Epistatic interactions can also be relevant to the development of complex human diseases such as diabetes [4]. Complex traits and diseases are determined by a multiplicity of genomic loci [5], whose independent effects and interactions [6] are often necessary to understand the phenotype of interest. Despite the broad implications of epistatic interactions, a quantitative characterization of their typical strength is still lacking. In this study, we consider growth rate in yeast as an example of a complex trait modulated by genetic interactions.

Previous studies [7-10] on the relation between the growth effects of a mutation and its epistatic interactions have often been based on a handful of mutations, and only in recent years has anecdotal evidence started being replaced by robust statements based on large data sets. Perhaps the most impressive of these datasets is the one made publicly available with the publication of the article entitled 'The genetic landscape of a cell' by Costanzo et al. [11]. The genome of the budding yeast Saccharomyce cerevisiae includes approximately 6,000 genes, about 1,000 of which are essential. Viable mutants can be constructed by knocking out any of the approximately 5,000 non-essential genes, by reducing the expression of the essential genes, or by partially compromising the functionality of the gene products. The dataset (see Additional file 1, Figure S1) has been compiled with the growth rates of about 5.4 million double knockout mutants, a sizable fraction of all possible double knockout mutants in yeast. Supported by the Costanzo et al. dataset, we consider the fundamental question of whether mutations with larger effects have stronger genetic interactions.

Additional file 1. Supplementary Figures and Text.

Format: PDF Size: 3.5MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Results and discussion

An unbiased definition of genetic interactions

A basic approach to study genetic interactions is to consider two mutations with known effects on a quantitative trait, and to measure their combined effect in the double mutant [12]. Given [11,13] the growth rates of a wild type S. cerevisiae strain (g00 = 1) and of two single knockout mutants (g01 and g10), the growth rate of the double knockout mutant (g11) is adequately predicted by a multiplicative null model:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M1">View MathML</a>

Equivalently, defining 'log growth' as the logarithm of the relative growth rate,

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M2">View MathML</a>

the log growth of the double knockout mutant is predicted by an additive null model (Figure 1a):

thumbnailFigure 1. The log growth rates of two mutations combine additively. (a) The average effect of a double knockout (G11) as a function of the effects of the single knockouts (G01 and G10) is G11 = G01 + G10. Experimental mean +/- standard deviation (blue line and blue shaded area) and prediction of the additive null model (red line). (b) Given two mutations, there are four possible mutants with their corresponding log growth rates (black dots). If three of the four log growth rates are known, the fourth one can be predicted by a linear extrapolation (red plane), and epistasis can be defined as the linear deviation from such prediction (red arrow). The magnitude of the deviation is the same regardless of which three of four mutants are chosen.

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M3">View MathML</a>

Epistatic interactions are identified as deviations from the null model, but several non-equivalent alternatives exist for quantifying these deviations [14]. The most common definition of epistasis considers the difference between the measured and the predicted growth rates for the double knockout mutant [11]:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M4">View MathML</a>

Importantly, this definition of e subtly constrains the possible values of epistasis. In fact, when combining very deleterious mutations, e cannot be large and negative even when the double knockout mutant is a synthetic lethal mutant:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M5">View MathML</a>

In order to avoid a priori constraints on the intensity of epistasis, genetic interactions can be defined as the ratio between the measured and predicted relative growth rates, leading to:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M6">View MathML</a>

As an example, E = +1 indicates a double mutant whose growth rate is twice as large as would be expected based upon the multiplicative null model, whereas E = -1 indicates a double mutant whose growth rate is half as large as predicted. This definition of epistasis as fold deviation in the multiplicative model for growth rates is equivalent to a natural definition of epistasis as linear deviation in the additive model for log growth rates (Figure 1b):

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M7">View MathML</a>

A second bias of the common definition of epistasis is that e depends on the choice of which genotype is labeled as 'wild type' or '00', a choice which is always arbitrary, but more obviously so when studying engineered organisms or populations evolving in alternating environments [15]. By contrast,

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M8">View MathML</a>

depends only on which pair of genes is considered, being a geometric measure for the 'curvature' of the fitness landscape (Figure 1b).

The definition of E has found some favor in the theoretical literature [7,16], but it is not routinely used to analyze experimental data apart from rare exceptions [8,17]. Its main drawback is that synthetic lethals have a log growth rate of -∞, and require a separate although simpler analysis in which lethal interactions can simply be counted. The definition of E proves instead to be extremely valuable when quantifying the strength of non-lethal genetic interactions.

Epistatic interactions scale weakly with mutational effects

With the appropriate definition of epistasis, a simple relation between the growth rate effects of two mutations and the expected strength of their interaction emerges.

Let us consider two groups of mutations; in the first group, all mutations have log growth effect G01, and in the second group, all mutations have log growth effect G10. We can then build all possible double mutants obtained by combining one mutation from each group. In the absence of epistasis, all the double mutants have a log growth rate

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M9">View MathML</a>

and the distribution of genetic interactions is sharply peaked at E = 0. When epistasis is present, the distribution of genetic interactions has, in general, non-zero mean and standard deviation. Experimentally, however, the mean of genetic interactions is close to zero (this is why the null model remains approximately valid) (Figure 1a; Figure 2d). Even when the mean interaction is vanishing, the difference between the experimental dataset and the ideal case without interactions can be quantified by the finite value of the experimental standard deviation σ(G01, G10), which provides a numerical estimate for the characteristic strength of epistatic interactions.

thumbnailFigure 2. The strength of epistatic interactions scales with the log growth effects of the interacting knockouts. (a) Each dot represents the variance of several thousand epistatic interactions binned according to the log growth effects of the two single knockouts, G01 and G10. The blue surface is the phenomenological fit:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M27">View MathML</a>

(b) Slices of the plot in (a) for G01 = constant. The dots are the same as in (a), and the solid lines represent the corresponding slice of the one-parameter fitting surface. (c) Diagonal slice of the plot in (a) with finer bins (G01 = G10 within 20%, G = mean(G01, G10)). The blue shaded area is the 25 to 75% confidence interval computed by bootstrap; the red line (var(G, G) = 0.079 G) is computed from the phenomenological model, and the dashed gray line, for which var(G, G) is proportional to G2, represents the lower bound to the slope predicted by the Fisher's geometric model. (c, inset) The epistatic interactions between beneficial mutations are vanishingly small, independently of the effect of the combined mutations. (d) Probability density functions p(E') for the strength of genetic interactions between two deleterious knockouts with similar log growth effects. Different colors correspond to knockouts with different effects: the growth rates effects of the single knockouts being combined are close to -38% (red), -22% (yellow), -12% (green), -6% (blue), and -3% (purple). Each curve has been rescaled so that all distributions have a standard deviation = 1. The left tail of the distributions displays a fat tail, describing the occurrence of strong negative genetic interactions (for comparison, the dashed-dotted black line is a normal distribution).

In order to produce reliable numerical results, thousands of growth rates are necessary to characterize the probability distribution of epistasis. We analyzed the Costanzo et al. dataset by binning pairs of mutations according to the log growth effects of their single knockouts G01 and G10, using the method described above to outline the probability distribution of epistasis. We chose bin sizes that grow exponentially with G in order to ensure an approximately constant number of data points in each bin (see Materials and Methods; see Additional file 1, Figure S2). Most bins contain from thousands to tens of thousands of data points. For each bin, we computed

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M10">View MathML</a>

that is, the variance of the random variable E relative to the bin labeled by growth rates G01 and G10. In the rest of the paper we will refer to such variance as var(G01, G10), emphasizing that the variance in the strength of epistatic interactions is, eventually, a function of G01 and G10 (Figure 2a). The square root of the variance, σ(G01, G10), then represents the expected strength of epistasis as a function of the independently varying effects of the two single knockouts. A natural expectation for the dependence of epistasis on the effect of the combined mutations comes from rescaling Figure 1a; if all the log growth effects of single and double knockouts increase by a factor of two, then the strength of epistasis should also increase by a factor of two. Unexpectedly, however, when combining deleterious mutations, the strength of epistatic interactions does grow with the effects of the mutations that are combined, but the dependence is much weaker; when the effect of both single knockouts is doubled, the strength of epistasis increases only by a factor of √2 (Figure 2).

In more detail, we observed that if the effect of the first knockout (G01) is held constant, the dependence of the variance of epistasis on the effect of the second knockout (G10) is well approximated by a Michaelis-Menten law (Figure 2b):

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M11">View MathML</a>

When the effects of both knockouts are free to vary, the requirement that the variance is a symmetric function of its two variables, G01 and G10, implies that K = |G01| and that v is proportional to G01. A one-parameter function which fits the seen variance over the whole range of deleterious fitness effects (Figure 2a) is then:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M12">View MathML</a>

This functional form can also be obtained from a simple model based on diffusion in fitness space (see Additional file 1, Supplementary text 1). An even simpler phenomenological fit, although slightly less accurate, is:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M13">View MathML</a>

(see Additional file 1, Figure S3). Importantly, these functions capture two major features of the data; first, epistasis vanishes when G01 or G10 = 0; second, when the effects of the two knockouts are similar (G01 = G10 = G along the diagonal of the surface in Figure 2a), the variance of epistasis is approximately proportional to G (Figure 2c):

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M14">View MathML</a>

The scaling described above is seen only for deleterious knockouts. When combining the beneficial knockouts in the dataset instead, the strength of epistasis is close to zero (Figure 2c, inset). This might be because the slightly beneficial knockouts are not adaptive mutations, but simply remove genes that are not needed in the conditions chosen for the experiment, so that their interactions are likely to be negligible. However, in apparent contrast to this observation, recent studies [8,18] on adaptive mutations in Escherichia coli suggest that genetic interactions between adaptive mutations are mostly negative. In fact, during adaptation, the prevalence of negative interactions is likely to be caused by biased sampling, because the mutations that fix in the population are likely to be the ones that solve environmental or biological challenges for an organism. Diminishing returns arise because the appearance of multiple 'solutions' to the same challenge is not necessarily preferable over the presence of a single solution. Rather than focusing on mutations that fix during a bout of adaptation, the Costanzo et al. dataset includes a large fraction of all possible pairs of genes in the yeast genome. Because for most pairs the two genes are involved in unrelated biological processes, interactions are often vanishingly small. We did observe, however, that the distribution of epistatic interactions is asymmetric, with a heavy tail of deleterious interactions (Figure 2d).

Experimental uncertainty generates spurious epistatic interactions

When inferring genetic interactions from experimental data, it is important to take into account that each measured growth rate is affected by some uncertainty, and that measurement errors in the growth rates could erroneously be interpreted as genetic interactions. Importantly, for each single and double mutant, the Costanzo et al. dataset provides the mean growth rate together with its estimated experimental uncertainty (the growth rate of each mutant being measured at least four times).

In order to quantify the effect of the experimental uncertainty on the inferred epistatic interactions, we constructed a number of mock datasets, assuming that the null model without epistatic interactions described biology exactly. In these datasets, each single knockout had the same growth rate as in the original dataset, and each double knockout had a growth rate equal to the product of the relative growth rates of the corresponding single knockouts. We then randomized the mock datasets by shifting each growth rate by a random amount sampled from a Student's t-distribution, with width depending on the corresponding experimental uncertainty reported in the original dataset (see Additional file 1, Supplementary text 3). As expected, analysis of these 'noisy' datasets revealed some epistasis, clearly caused by our addition of experimental noise rather than by any biological mechanism. We found that for pairs involving beneficial or neutral mutations, the variance computed in the mock datasets was comparable to or even greater than the variance observed in the original dataset (Figure 3a, black curves; Figure 3b, blue regions). This fact provides an important internal control, suggesting that the experimental noise has not been underestimated. In spite of this, for pairs of knockouts with substantially deleterious effects, experimental noise accounted for less than half of the total observed variance, with the rest representing genuine biological interactions (Figure 3a, red curves; Figure 3b, red regions).

thumbnailFigure 3. Experimental noise does not account for all of the observed variance of epistasis. (a) Comparison of experimentally measured variance (solid lines; shaded areas: 25 to 75% confidence intervals) and variance caused by experimental noise (dashed lines). If one of the two mutations is neutral, noise accounts for all of the observed variance (black). When deleterious mutations are combined, noise accounts for less than half of the observed variance (red, G01 ≈ -0.7). (b) Ratio between total observed variance and noise-generated variance as a function of the log growth of the knockouts being combined. For deleterious knockouts, the ratio can be significantly greater than 1.

We then decomposed the variance observed in the original dataset into a contribution produced by experimental uncertainty and a contribution of biological origin; the strength of epistatic interactions was finally computed as the square root of the biological part of the variance. For deleterious knockouts, the relative difference between epistasis computed from the raw data and from the data after subtracting the experimental noise was less than 30%, emphasizing the significant but not overwhelming contribution of experimental noise to the observed variability. Figure 2(a-c) represents the 'biological' part of the observed epistasis; before subtracting the contribution of the experimental uncertainty, the plots are qualitatively similar, but quantitatively slightly different (see Additional file 1, Figure S4). Importantly, because variances are additive, the estimated contribution of the experimental uncertainty to epistasis is largely independent of the choice of the statistical distribution used to model experimental uncertainty. In two instances, however, the unknown details of the full distribution of experimental noise are important; when outlining the distribution of epistatic interactions (Figure 2d) and when describing the probability to observe sign epistasis (Figure 4b). In those two figures, we plotted the raw data, and did not attempt to deconvolve the contribution of experimental uncertainty.

thumbnailFigure 4. Sign epistasis is less likely to occur between mutations with large effects. (a) Examples of a smooth landscape with paths of monotonically increasing fitness (left) and a rugged landscape characterized by reciprocal sign epistasis (right). (b) Experimentally measured probability of observing sign epistasis as a function of the log growth of two single knockouts with similar effects (G01 = G10 within 20%, G = mean(G01, G10)). The blue shaded area is the standard error of the mean computed by bootstrap.

Comparison between theory and experiment

The scaling of epistasis observed in the Costanzo et al. dataset (Figure 2) is in sharp contrast to the predictions of Fisher's geometric model [19], a popular model of epistasis in which genetic interactions emerge from geometry. As we saw, when the effects of the two knockouts are similar (G01 = G10 = G), the variance of epistasis is approximately proportional to G. By contrast, in the Fisher's model, the variance var(G, G) would grow faster than G2 (Figure 2c; see Additional file 1, Supplementary text 2), a much stronger dependence than the linear dependence observed experimentally.

A concrete numerical example can highlight the importance of the weaker-than-expected scaling of epistasis described in this study. Let us consider two gene knockouts, each of which reduces the relative growth rate by 5%, from 1.0 to 0.95. According to the multiplicative null model, the growth rate of the double knockout will be approximately 0.952, or approximately 0.90. The questions now are: What kind of deviations could be expected around 0.90? Would a growth rate of 0.85 be surprising? What about a growth rate of 0.50?. Let us use the analytic fit discussed in the previous section

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M15">View MathML</a>

Then

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M16">View MathML</a>

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M17">View MathML</a>

and

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M18">View MathML</a>

A +/- one standard deviation interval for the growth rate of the double knockout is then

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M19">View MathML</a>

Notice that it is not unlikely that epistasis will cancel the effect of the second mutation, so that the growth rate of the double knockout mutant is greater than 0.95, that is, greater than the growth rate of either of the single knockout mutants.

Let us now consider two gene knockouts with stronger effects, each of which reduces the growth rate from 1.0 to 0.60. Then

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M20">View MathML</a>

about 10 times as large as the log growth of the single mutants in the previous example. The Fisher's model would predict a σ(G, G) at least 10 times larger than in the previous example (σ(G, G)≥0.76), and an interval of likely growth rates for the double knockout mutants at least as large as

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M21">View MathML</a>

Notice how, once again, it is not unlikely that owing to genetic interactions, the growth rate of the double knockout mutant is greater than 0.60, the growth rate of either of the two single knockout mutants. The analytic model derived from the experimental data leads to a strikingly different conclusion:

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M22">View MathML</a>

and the +/- one standard deviation interval for the growth rate of the double knockout becomes

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M23">View MathML</a>

In this case, a deviation from the null model that is greater than three standard deviations would be needed for the double knockout mutant to have a growth rate greater than that of the single knockout (0.60), making the event extremely unlikely.

Epistasis constrains the evolutionary dynamics

The previous section provided two examples of reciprocal sign epistasis, realized when two deleterious mutations produce a double mutant that is fitter than either of the two single mutants (Figure 4a). In those cases, a fitness valley limits the evolutionary accessibility of the fitter double mutant, and only on longer time scales may the simultaneous appearance of two mutations [20,21] drive a population to the new local fitness maximum. In this context, the scaling behavior of epistasis is of great importance, because it determines the number and the topology of the evolutionarily accessible paths [2,22,23], ultimately affecting the possible outcomes of the evolutionary process.

In order to describe how epistasis shapes the naturally occurring fitness landscapes, let us consider S(G, G), the probability to observe sign epistasis when combining two mutations with similar growth rate effects, G. Here, S(G, G) depends on the typical interaction strength,

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M24">View MathML</a>

In particular, if σ(G, G) is proportional to G, then the probability of observing sign epistasis is independent of G. The Fisher's model implies a super-linear dependence of σ(G, G) on G, thus predicting a greater probability of observing sign epistasis among mutations with strong effects. Instead, if the scaling of σ(G, G) is proportional to √G (Figure 2), then sign epistasis is more likely to occur among mutations with small effects (Figure 3b). When the relative growth rate effects of the single knockouts are small (<2 to 3%), experimental uncertainty prevents us from pinpointing which pairs of genes are epistatic. This does not mean, however, that mutations with small effects do not interact. Assuming that the scaling of epistasis we measured directly for mutations with intermediate and large effects extends to mutations with small effects, a consequence of the observed scaling of epistasis is the roughening of the local fitness landscape in the proximity of an evolutionary optimum; when the fitness effects of available mutations become small [24], epistatic interactions become increasingly relevant [25,26], reducing the accessibility of evolutionary paths and further slowing down the rate of adaptation [27,28]. The evolutionary dynamics on correlated fitness landscapes [10,29] with the realistic correlations described here certainly deserves further experimental and theoretical investigation.

The scaling of genetic interactions may be generic

To date, our analysis has been limited to interactions between entire gene knockouts. Although mutations with extreme effects on gene regulation and horizontal gene transfer are biologically relevant mechanisms for the removal or acquisition of whole genes at once, organisms explore possible genetic variants largely through the accumulation of single point mutations. The Costanzo et al. data et contains thousands of double mutants for which the first mutation is a gene knockout and the second mutation consists of one or more point mutations in a different gene, causing the gene product to misfold in a temperature-sensitive way. Although the distribution of growth rate effects for point mutations is different than for single gene knockouts (see Additional file 1, Figure S2), the statistics of genetic interactions are remarkably similar when combining two single knockouts and when combining a single knockout with a point mutation (Figure 5). A similar scaling is also seen for the epistatic interactions between single gene knockouts and decreased abundance by mRNA perturbation [30] (DamP) perturbations of a second gene (see Additional file 1 Figure S5). The analysis of these hybrid double mutants suggests that the statistics of the interactions between any two genetic perturbations are determined only by their growth rate effects [31], and not by their biological origin in terms of point mutations or gene knockouts.

thumbnailFigure 5. Point mutations have similar epistatic interactions to those of entire gene knockouts. (a) Comparison between the variance observed in double gene knockout mutants (rainbow dots, same as in Figure 2a) and the variance observed in mixed double mutants generated by combining a gene knockout with point mutations in a different gene (black dots). (b) The red curve is the diagonal slice of the plot in (a) (G01 = G10 within 20%, G = mean(G01, G10)), and the red shaded area is the 25 to 75% confidence interval for the mixed double mutant variance. For comparison, the blue curves describe the variance for double gene knockouts as in Figure 2c. As in Figure 2c, the red line has equation var(G, G) = 0.079 G.

A comparison between different definitions of epistasis

Importantly, any quantitative result on epistasis is a consequence of how epistasis is defined. Of particular interest is how strong an epistatic interaction is deemed to be, based upon its ranking when compared with that of other pairs of mutations. Although the 'traditional' definition

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M25">View MathML</a>

and the 'geometric' definition

<a onClick="popup('http://genomebiology.com/2013/14/7/R76/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://genomebiology.com/2013/14/7/R76/mathml/M26">View MathML</a>

agree about the sets of positive and negative interactions, they assign different strengths and, more importantly, different rankings to the same pair of interacting mutations. As an example, if the Costanzo et al. dataset is analyzed using the 'traditional' definition of genetic interactions, then the linear dependence of var(G, G) on G in Figure 2c is replaced by an oddly non-monotonic dependence, displaying weaker interactions for pairs of genes with either very small or very large fitness effects (Figure 6a). As mentioned previously, this decrease in the inferred strength of epistatic interactions for very deleterious mutations is a mathematical consequence of the traditional definition of epistasis, rather than a property of genetic interactions. The same bias would lead us to conclude that genes with strong effects on growth are almost non-interacting (Figure 6b, red line). However, because previous studies have determined that essential genes partake in more interactions than do non-essential genes [32], it is also reasonable to expect that non-lethal genes with large growth effects are involved in more interactions than genes with small growth effects. Indeed, according to the 'geometric' definition of epistasis, the fraction of genes with which a gene interacts steadily increases with the growth rate effect of the gene (Figure 6b, blue line). By contrast, the traditional definition of epistasis, consistently assigns low rankings to interactions between genes with large growth rate defects, as confirmed by a further analysis comparing the two definitions of epistasis against interactions inferred from the Gene Ontology (GO) database [33] (see Additional file 1, Figure S6). According to the geometric definition of epistasis, genetic networks [34] are denser than expected not only among essential gene [32], but also among genes with large growth effects.

thumbnailFigure 6. Comparison between the traditional and the geometric definitions of epistasis (e and E, respectively). (a) Figure equivalent to Figure 2c, using the traditional definition of epistasis. (b) The fraction of genes interacting with a specific gene is a function of the growth rate effect of such gene. Only the 10,000 most interacting pairs the geometric definition (blue) and the traditional definition (red) are considered to be interactions.

Finally, it is important to emphasize that the traditional definition of epistasis remains slightly more successful at discovering the functional relations between genes, as cataloged in the GO database (see Additional file 1, Figure S6). Part of the reason for this could be that some of those functional characterizations were suggested by the traditional definition of epistasis in the first place. It is certainly true, however, that many of the top-ranking interactions according to the geometric definition of epistasis involve single and double mutants with small growth rates; for those mutants, experimental noise is relatively large, and this may cause a few weakly interacting pairs to be incorrectly ranked as strongly interacting. It is likely that the experimental protocols could be easily adjusted to reduce the relative uncertainty on the growth rate of especially slow-growing mutants to avoid this issue (for example, by allowing for a much longer time for growth or by measuring the growth rates of additional replicates).

Conclusions

We analyzed the growth rates of about five million double mutants in the dataset associated with the work by Costanzo et al. We characterized how the strength of genetic interactions depends on the growth effects of the mutations being combined, and found a weaker dependence than that predicted by current theoretical models. Although the results were obtained mainly from entire gene knockouts, there is some evidence that the observed scaling might extend to the interactions between single point mutations. The scaling of epistasis might or might not be generic [35,36]; important drivers could be the harshness of the environment [37], details about the evolutionary history [38-40], sexual versus asexual reproduction [41] and, perhaps most importantly, metabolic [42-45] and genetic complexity [46,47]. In general, the experimentally observed scaling suggests a previously unexplored class of correlated fitness landscapes with tunable roughness, in which epistasis depends explicitly on the effects of the mutations being combined.

A clear limitation of our discussion is that only pair interactions were considered. Although high-throughput experiments will provide data on higher-order interactions, a solid understanding of pair interactions remains necessary before addressing n-mutation interactions. A genuine three-mutation interaction, for instance, should be defined as the unexplained deviation from what can be computed by combining the effects of all relevant mutations and their pair interactions [10,48], perhaps using linear fits within the additive null model for log growth rates.

The results we present here were based on a geometric definition of epistasis. We compared this definition with a more standard definition, highlighting the desirable mathematical properties of the geometric definition and the simple phenomenological relations it produces.

In conclusion, although each epistatic interaction between specific genes depends on biological details and remains largely unpredictable from first principles, we have shown that the statistical properties of an ensemble of interactions can display conspicuous regularities, and can be described by simple mathematical laws.

Materials and methods

The Costanzo et al. dataset is publicly available [49]. The file http://sgadata_costanzo2009_rawdata_101120.txt.gz webcite was downloaded on August 17, 2010 and analyzed with Mathematica (code available at the Gore laboratory website [50]). We restricted our analysis to double knockout mutants whose growth rates were positive numerical values and for which the growth rates of both single mutants were numerical values (see Additional file 1, Figure S1). Some genes appear in the dataset both as query and array genes; care was taken to avoid double counting.

The exponentially growing intervals used for the binning of the log growth rate effects were defined as [-2n, -2n-1] for an appropriate range of integer n's. Owing to the rarity of extremely deleterious mutations, bins for positive n's contained only a few data points, while bins with large negative n's were extremely small. In the figures we reported only bins for n = -7 to 0, containing log growth rate effects ranging from -20 = -1 to -2-8 = -0.0039 or, alternatively, relative growth rate effects ranging from 2-1 = 0.5 to 2-0.0039 = 0.997. Different choices for the binning sizes and positions did not significantly alter the results of the analysis.

In order to quantify the contribution of experimental uncertainty to epistasis, we generated nine randomized mock datasets. The mean level of noise-generated epistasis in these nine datasets is reported in Figure 4 (dashed lines), and we provide an extensive discussion of the choice of Student's t-distributions to generate the mock datasets from the original dataset (see Additional file 1, Supplementary text 3).

The GO database go_201207-assocdb-tables.tar.gz was downloaded from the GO site [51] on July 19, 2012. The MySQL database was queried with Python and analyzed Mathematica (code available upon request).

Abbreviations

DamP: Decreased abundance by mRNA perturbation; GO: Gene Ontology.

Conflict of interest

The authors declare that they have no conflict of interest.

Authors' contributions

AV and JG designed research; AV performed research and analyzed data; AV and JG wrote the paper. Both authors have read and approved the final manuscript.

Acknowledgements

We are grateful to Mingjie Dai for collaboration during the early stages of the study. We thank Kirill Korolev, Pankaj Mehta, and the members of the Gore laboratory for providing comments and advice on the manuscript. This research was funded by an NIH Pathways to Independence Award, NSF CAREER Award, Pew Biomedical Scholars Program, and Alfred P. Sloan Foundation Fellowship.

References

  1. Phillips PC: Epistasis - The essential role of gene interactions in the structure and evolution of genetic systems.

    Nat Rev Genet 2008, 9:855-867. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Weinreich DM, Delaney NF, DePristo MA, Hartl DL: Darwinian evolution can follow only very few mutational paths to fitter proteins.

    Science 2006, 312:111-114. PubMed Abstract | Publisher Full Text OpenURL

  3. Dettman JR, Sirjusingh C, Kohn LM, Anderson JB: Incipient speciation by divergent adaptation and antagonistic epistasis in yeast.

    Nature 2007, 447:585-588. PubMed Abstract | Publisher Full Text OpenURL

  4. Hoh J, Ott J: Mathematical multi-locus approaches to localizing complex human trait genes.

    Nat Rev Genet 2003, 4:701-709. PubMed Abstract | Publisher Full Text OpenURL

  5. Mackay TFC, Stone EA, Ayroles JF: The genetics of quantitative traits: challenges and prospects.

    Nat Rev Genet 2009, 10:565-577. PubMed Abstract | Publisher Full Text OpenURL

  6. Jansen RC: Studying complex biological systems using multifactorial perturbation.

    Nat Rev Genet 2003, 4:145-151. PubMed Abstract | Publisher Full Text OpenURL

  7. Gros PA, Le Nagard H, Tenaillon O: The evolution of epistasis and its links with genetic robustness, complexity and drift in a phenotypic model of adaptation.

    Genetics 2009, 182:277-293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Khan AI, Dinh DM, Schneider D, Lenski RE, Cooper TF: Negative epistasis between mutations in an evolving bacterial population.

    Science 2011, 332:1193-1196. PubMed Abstract | Publisher Full Text OpenURL

  9. Wilke CO, Adami C: Interaction between directional epistasis and average mutational effects.

    Proc R Soc Lond B Biol Sci 2001, 268:1469-1474. Publisher Full Text OpenURL

  10. Beerenwinkel N, Pachter L, Sturmfels B, Elena SF, Lenski R: Analysis of epistatic interactions and fitness landscapes using a new geometric approach.

    BMC Evol Biol 2007, 7:60-73. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. Costanzo M, Baryshnikova A, Bellay J, Kim Y, Spear ED, Sevier CS, Ding H, Koh JLY, Toufighi K, Mostafavi S, Prinz J, St Onge RP, VanderSluis B, Makhnevych T, Vizeacoumar FJ, Alizadeh S, Bahr S, Brost RL, Chen Y, Cokol M, Deshpande R, Li Z, Lin Z-Y, Liang W, Marback M, Paw J, San Luis B-J, Shuteriqi E, Tong AHY, van Dyk N, et al.: The genetic landscape of a cell.

    Science 2010, 327:425-431. PubMed Abstract | Publisher Full Text OpenURL

  12. Dixon SJ, Costanzo M, Baryshnikova A, Andrews B, Boone C: Systematic mapping of genetic interaction networks.

    Annu Rev Genet 2009, 43:601-625. PubMed Abstract | Publisher Full Text OpenURL

  13. Baryshnikova A, Costanzo M, Kim Y, Ding H, Koh J, Toufighi K, Youn J-Y, Ou J, San Luis B-J, Bandyopadhyay S, Hibbs M, Hess D, Gingras A-C, Bader GD, Troyanskaya OG, Brown GW, Andrews B, Boone C, Myers CL: Quantitative analysis of fitness and genetic interactions in yeast on a genome scale.

    Nat Methods 2010, 7:1017-1024. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Mani R, StOnge RP, Hartman JL, Giaever G, Roth FP: Defining genetic interaction.

    Proc Natl Acad Sci USA 2008, 105:3461-3466. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Tan L, Gore J: Slowly switching between environments facilitates reverse evolution in small populations.

    Evolution 2012, 66:3144-3154. PubMed Abstract | Publisher Full Text OpenURL

  16. Peters AD, Lively CM: Epistasis and the maintenance of sex. In Epistasis and the evolutionary process. Edited by Wolf JB, Brodie ED, Wade MJ. Oxford: Oxford University Press; 2000:99-112. OpenURL

  17. Martin G, Elena SF, Lenormand T: Distributions of epistasis in microbes fit predictions from a fitness landscape model.

    Nat Genet 2007, 39:555-560. PubMed Abstract | Publisher Full Text OpenURL

  18. Chou HH, Chiu HC, Delaney NF, Segre D, Marx CJ: Diminishing return epistasis among beneficial mutations decelerates adaptation.

    Science 2011, 332:1190-1192. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Fisher RA: The Genetical Theory of Natural Selection. Oxford: Clarendon Press; 1930.

  20. Weinreich DM, Watson RA, Chao L: Sign epistasis and genetic constraint on evolutionary trajectories.

    Evolution 2005, 59:1165-1174. PubMed Abstract OpenURL

  21. Weissman DB, Desai MM, Fisher DS, Feldman MW: The rate at which asexual populations cross fitness valleys.

    Theor Pop Biol 2009, 75:286-300. Publisher Full Text OpenURL

  22. Poelwijk FJ, Kiviet DJ, Weinreich DM, Tans ST: Empirical fitness landscapes reveal accessible evolutionry paths.

    Nature 2007, 445:383-386. PubMed Abstract | Publisher Full Text OpenURL

  23. Velenich A, Gore J: Synthetic approaches to understanding biological constrains.

    Curr Opin Chem Biol 2012, 16:323-328. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Eyre-Walker A, Keightley PD: The distribution of fitness effects of new mutations.

    Nat Rev Genet 2007, 8:610-618. PubMed Abstract | Publisher Full Text OpenURL

  25. Tan L, Serene S, Chao HX, Gore J: Hidden randomness between fitness landscapes limits reverse evolution.

    Phys Rev Lett 2011, 106:198102. PubMed Abstract | Publisher Full Text OpenURL

  26. Woods RJ, Barrick JE, Cooper TF, Shrestha U, Kauth MR, Lenski RE: Second-order selection for evolvability in a large Escherichia coli population.

    Science 2011, 331:1433-1436. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Orr HA: The population genetics of adaptation the distribution of factors fixed during adaptive evolution.

    Evolution 1998, 52:935-949. Publisher Full Text OpenURL

  28. Orr HA: The genetic theory of adaptation a brief history.

    Nat Rev Genet 2005, 6:119-127. PubMed Abstract | Publisher Full Text OpenURL

  29. Kryazhimskiy S, Tkčik G, Plotkin JB: The dynamics of adaptation on correlated fitness landscapes.

    Proc Natl Acad Sci USA 2009, 106:18638-18643. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Breslow DK, Cameron DM, Collins SR, Schuldiner M, Stewart-Ornstein J, Newman HW, Braun S, Madhani HD, Krogan NJ, Weissman JS: A comprehensive strategy enabling high-resolution functional analysis of the yeast genome.

    Nat Methods 2008, 5:711-718. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Xu L, Barker B, Gu Z: Dynamic epistasis for different alleles of the same gene.

    Proc Natl Acad Sci USA 2012, 109:10420-10425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Davierwala AP, Haynes J, Li Z, Brost RL, Robinson MD, Yu L, Mnaimneh S, Ding H, Zhu H, Chen Y, Cheng X, Brown GW, Boone C, Andrews BJ, Hughes TR: The synthetic genetic interaction spectrum of essential genes.

    Nat Genet 2005, 37:1147-1152. PubMed Abstract | Publisher Full Text OpenURL

  33. Barabási AL, Oltvai ZN: Network biology: understanding the cell's functional organization.

    Nat Rev Genet 2004, 5:101-113. PubMed Abstract | Publisher Full Text OpenURL

  34. The Gene Ontology Consortium: Gene ontology: tool for the unification of biology.

    Nat Genet 2000, 25:25-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Dixon SJ, Fedyshyn Y, Koh JLY, Keshava Prasad TS, Chahwan C, Chua G, Toufighi K, Baryshnikova A, Hayles J, Hoe K-L, Kim D-U, Park H-O, Myers CL, Pandey A, Durocher D, Andrews BJ, Boone C: Significant conservation of synthetic lethal genetic interaction networks between distantly related eukaryotes.

    Proc Natl Acad Sci USA 2008, 105:16653-16658. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Tischler J, Lehner B, Fraser AG: Evolutionary plasticity of genetic interaction networks.

    Nat Genet 2008, 40:390-391. PubMed Abstract | Publisher Full Text OpenURL

  37. Harrison R, Papp B, Pál C, Oliver SG, Delneri D: Plasticity of genetic interactions in metabolic networks of yeast.

    Proc Natl Acad Sci USA 2007, 104:2307-2312. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Wagner A: Gene duplications robustness and evolutionary innovations.

    BioEssays 2008, 30:367-373. PubMed Abstract | Publisher Full Text OpenURL

  39. Wagner A: Distributed robustness versus redundancy as causes of mutational robustness.

    BioEssays 2005, 27:176-188. PubMed Abstract | Publisher Full Text OpenURL

  40. Roguev A, Bandyopadhyay S, Zofall M, Zhang K, Fischer T, Collins SR, Qu H, Shales M, Park H-O, Hayles J, Hoe K-L, Kim D-U, Ideker T, Grewal SI, Weissman JS, Krogan NJ: Conservation and rewiring of functional modules revealed by an epistasis map in fission yeast.

    Science 2008, 322:405-410. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Azevedo RBR, Lohaus R, Srinivasan S, Dang KK, Burch CL: Sexual reproduction selects for robustness and negative epistasis in artificial gene networks.

    Nature 2006, 440:87-90. PubMed Abstract | Publisher Full Text OpenURL

  42. Segrè D, DeLuna A, Church GM, Kishony R: Modular epistasis in yeast metabolism.

    Nat Genet 2005, 37:77-83. PubMed Abstract | Publisher Full Text OpenURL

  43. He X, Qian W, Wang Z, Li Y, Zhang J: Prevalent positive epistasis in Escherichia coli and Saccharomyces cerevisiae metabolic networks.

    Nat Genet 2010, 42:272-276. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Szappanos B, Kovács K, Szamecz B, Honti F, Costanzo M, Baryshnikova A, Gelius-Dietrich G, Lercher MJ, Jelasity M, Myers CL, Andrews BJ, Boone C, Oliver SG, Pál C, Papp B: An integrated approach to characterize genetic interaction networks in yeast metabolism.

    Nat Genet 2011, 43:656-662. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Almaas E, Kovács B, Vicsek T, Oltvai ZN, Barabási AL: Global organization of metabolic fluxes in the bacterium Escherichia coli.

    Nature 2004, 427:839-843. PubMed Abstract | Publisher Full Text OpenURL

  46. Sanjuán R, Elena SF: Epistasis correlates to genomic complexity.

    Proc Natl Acad Sci USA 2006, 103:14402-14405. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Sanjuán R, Nebot MR: A network model for the correlation between epistasis and genomic complexity.

    PLoS One 2008, 3:e2663. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Wood K, Nishida S, Sontag ED, Cluzel P: Mechanism-independent method for predicting response to multidrug combinations in bacteria.

    Proc Natl Acad Sci USA 2012, 109:12254-12259. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. The Genetic Landscape of the Cell. [http://drygin.ccbr.utoronto.ca/~costanzo2009] webcite

  50. Gore Laboratory. [http://www.gorelab.org/software.html] webcite

  51. Gene Ontology. GO Database Downloads. [http://www.geneontology.org/GO.downloads.database.shtml] webcite