Analysis of epistatic interactions and fitness landscapes using a new geometric approach
- Niko Beerenwinkel^{1}Email author,
- Lior Pachter^{2},
- Bernd Sturmfels^{2},
- Santiago F Elena^{3} and
- Richard E Lenski^{4}Email author
DOI: 10.1186/1471-2148-7-60
© Beerenwinkel et al; licensee BioMed Central Ltd. 2007
Received: 17 October 2006
Accepted: 13 April 2007
Published: 13 April 2007
Abstract
Background
Understanding interactions between mutations and how they affect fitness is a central problem in evolutionary biology that bears on such fundamental issues as the structure of fitness landscapes and the evolution of sex. To date, analyses of fitness landscapes have focused either on the overall directional curvature of the fitness landscape or on the distribution of pairwise interactions. In this paper, we propose and employ a new mathematical approach that allows a more complete description of multi-way interactions and provides new insights into the structure of fitness landscapes.
Results
We apply the mathematical theory of gene interactions developed by Beerenwinkel et al. to a fitness landscape for Escherichia coli obtained by Elena and Lenski. The genotypes were constructed by introducing nine mutations into a wild-type strain and constructing a restricted set of 27 double mutants. Despite the absence of mutants higher than second order, our analysis of this genotypic space points to previously unappreciated gene interactions, in addition to the standard pairwise epistasis. Our analysis confirms Elena and Lenski's inference that the fitness landscape is complex, so that an overall measure of curvature obscures a diversity of interaction types. We also demonstrate that some mutations contribute disproportionately to this complexity. In particular, some mutations are systematically better than others at mixing with other mutations. We also find a strong correlation between epistasis and the average fitness loss caused by deleterious mutations. In particular, the epistatic deviations from multiplicative expectations tend toward more positive values in the context of more deleterious mutations, emphasizing that pairwise epistasis is a local property of the fitness landscape. Finally, we determine the geometry of the fitness landscape, which reflects many of these biologically interesting features.
Conclusion
A full description of complex fitness landscapes requires more information than the average curvature or the distribution of independent pairwise interactions. We have proposed a mathematical approach that, in principle, allows a complete description and, in practice, can suggest new insights into the structure of real fitness landscapes. Our analysis emphasizes the value of non-independent genotypes for these inferences.
Background
Understanding the structures of fitness landscapes is central to evolutionary biology. The image of populations evolving on fitness landscapes traces to Sewall Wright's seminal work in the thirties [3]. Since then, several types of fitness landscapes have been discussed in the literature [3–8], resulting in some confusion. The surface of a landscape may represent the relative fitness of individual types, or the average fitness of a population. In the former case, the underlying coordinates describe either the genotypic or phenotypic state of an individual; in the latter case, the coordinates describe either gene frequencies or average phenotypes in a population. Our paper concerns the mathematical analysis and interpretation of fitness landscapes where the height of the surface represents the relative fitness of individuals and the coordinates are different genotypes. In this evolutionary context, fitness measures the expected reproductive success of an individual having a specific genotype in some particular environment. Thus, a fitness landscape is given by assigning to each genotype g its fitness w_{ g }.
If all mutations were strictly additive or multiplicative in their effects on fitness, then it would be rather easy to describe the structure of fitness landscapes and understand the resulting dynamics of adaptation by natural selection. However, many mutations interact with one another in complex ways. For example, two or more mutations may interact such that their combined effect on fitness is much greater or much less than predicted from their individual effects; their combined effect may even be opposite in sign to the expectation based on their individual effects. These deviations from simple expectations are called epistasis [8, 9], and they determine the shape of fitness landscapes, including curvature and ruggedness. Understanding the prevalence and mathematical forms of these epistatic interactions is important for many issues in evolutionary biology including the dynamics of adaptation and divergence [3, 8, 10–12], reproductive isolation and speciation [8, 12–16], the evolution of sexual reproduction [17–20], the robustness of organisms to developmental, environmental, and mutational perturbations [21–24], the persistence of drug-resistant pathogens [25, 26], and more.
Over the last decade, several studies have sought to examine the form and prevalence of epistatic interactions by measuring the fitness effects of numerous mutations alone and in combination in viruses, bacteria, fungi, and animals [2, 27–31]. To date, these analyses have focused on the overall directional curvature of fitness as a function of the number of mutations, on the distribution of pairwise interactions, or on both. In this paper, we introduce a new mathematical approach that allows a more complete description of multi-locus interactions in a fitness landscape. The benefits of this approach become evident from our analysis of a fitness landscape of E. coli that was obtained and first analyzed by Elena and Lenski [2]. They measured relative fitness values for 37 genotypes at 9 loci including the non-mutated parental strain or "wild-type", 9 single mutants, and 27 double mutants that each had a different pair of mutations. Although there are 36 (= 9·8/2) possible double mutants, only 27 were constructed due to limitations of the markers used for strain construction. No triple mutants or higher-order combinations were constructed in this experiment.
Genotype space.
w | a | b | c | r | s | t | x | y | z |
a | ar | as | at | ax | ay | az | |||
b | br | bs | bt | bx | by | bz | |||
c | cr | cs | ct | cx | cy | cz | |||
r | ar | br | cr | rx | ry | rz | |||
s | as | bs | cs | sx | sy | sz | |||
t | at | bt | ct | tx | ty | tz | |||
x | ax | bx | cx | rx | sx | tx | |||
y | ay | by | cy | ry | sy | ty | |||
z | az | bz | cz | rz | sz | tz |
Fitness data.
1.000 | 0.976 | 0.708 | 0.975 | 0.981 | 0.984 | 0.995 | 0.978 | 0.564 | 0.593 |
0.976 | 0.990 | 0.973 | 0.990 | 0.982 | 0.718 | 0.500 | |||
0.708 | 0.964 | 0.684 | 0.694 | 0.782 | 0.664 | 0.510 | |||
0.975 | 0.983 | 0.975 | 0.974 | 0.977 | 0.650 | 0.482 | |||
0.981 | 0.990 | 0.964 | 0.983 | 0.718 | 0.988 | 0.524 | |||
0.984 | 0.973 | 0.684 | 0.975 | 0.724 | 0.986 | 0.490 | |||
0.995 | 0.990 | 0.694 | 0.974 | 0.982 | 0.679 | 0.508 | |||
0.978 | 0.982 | 0.782 | 0.977 | 0.718 | 0.724 | 0.982 | |||
0.564 | 0.718 | 0.664 | 0.650 | 0.988 | 0.986 | 0.679 | |||
0.593 | 0.500 | 0.510 | 0.482 | 0.524 | 0.490 | 0.508 |
The goal of our analysis is to describe the geometry of the E. coli fitness landscape obtained by Elena and Lenski [2]. The fitness of a genotype was measured as the rate of population growth, expressed relative to the growth rate for the wild-type. These authors calculated pairwise interactions by considering the 27 equations of the following sort: w·ar - a·r, where w is the fitness of the wild-type, ar is the fitness of a double mutant, and a and r are fitness values for each associated single mutant. They found many significant deviations from zero, including several cases of both positive and negative epistasis. The key messages of our paper are that additional types of gene interactions exist even within this fitness landscape, and that geometric methods can be used to describe and analyze the system more exhaustively.
Our analysis of this landscape is based on the approach developed by Beerenwinkel et al. [1]. We identify a comprehensive set of 243 gene interactions that includes the 27 standard pairwise tests used by Elena and Lenski. We then compare the epistatic deviations calculated using the new tests to those obtained from the standard tests, and we use the new tests to extract previously unnoticed features of the fitness landscape. Specifically, we investigate how epistasis depends on the fitness loss associated with deleterious mutations. We also consider tests that provide a new perspective on the relative "mixing ability" of different mutations. Here, the mixing ability of any given mutation specifies whether its epistatic interactions with a set of other mutations tend to be positive or negative. Finally, we describe the geometry of the overall fitness landscape by focusing on the three-locus subsystems whose shapes correspond to the triangulations of the four genotopes shown in Figure 1, panels b-d. We also discuss how this geometric formulation reflects the underlying set of gene interactions.
Results
Markov basis of the interaction space
Our first point is that the genotype space in Table 1 allows many more tests of epistasis than the 27 standard tests performed by Elena and Lenski [2]. Notice that the standard test w·ar - a·r compares two genotypes having zero and two mutations with two others each having one mutation. The minimal Markov basis (see Methods for a mathematical definition) of the space of tables of the type displayed in Table 1 reveals an additional 216 non-standard tests of the following two sorts. First, there are 108 "double-double" tests that compare two double mutants with two other double mutants, holding the distribution of the mutant alleles constant, for example, ar·bs - as·br. Second, there are another 108 orthogonal "single-double" tests that compare one single mutant and one double mutant with another single mutant and another double mutant, again holding constant the allele frequencies, for example, a·br - b·ar. For both types of non-standard tests, the two genotypes on the right-hand side can be regarded as the products of recombination between the two genotypes on the left-hand side. Notice that the same relationship also holds true for the standard tests.
Experimental biologists (including two authors of this paper) are likely to raise three concerns about these non-standard tests. What biological insights can non-standard tests provide beyond those obtained using standard tests? Are these additional tests independent of the standard tests? What computational tools are available to perform such tests on other datasets?
As we show in the sections that follow, the non-standard tests are potentially useful in at least three respects. First, they allow one to focus attention on features of epistasis that are not quantifiable by the standard tests. For example, we perform non-standard tests of the "single-double" type to explore whether some mutations are better overall mixers than others. Second, non-standard tests span greater genetic distances than do pairwise tests, allowing more powerful analyses of the structure of fitness landscapes. For example, we use the "double-double" tests to test curvature at genetic distances of four, whereas standard tests allow curvature to be examined only at distances of two. Third, non-standard tests are an integral part of the complete geometric description of a fitness landscape. While this high-dimensional geometry is abstract and even foreign, we describe how biological features of gene interactions, such as mixing ability, are embedded in the geometric shape of the landscape.
The non-standard tests are not independent, in a statistical sense, of the pairwise tests or of one another because all the tests are calculated from the same underlying data. Nonetheless, this complication can be addressed by employing appropriate statistical methods (Tukey's jackknife, Bonferroni correction, etc.) to ensure that significance levels reflect the data structure. Regarding the availability of computational tools, we provide references to programs that automate the calculations of the Markov basis and perform the triangulations necessary to describe the geometry of landscapes, and these tools can be applied to other datasets. In supplementary material, we illustrate the use of these computational tools and show their output (see Additional files 1-5).
Comparing standard and non-standard epistasis terms overall
An important feature of the standard tests is that the sign of epistasis, either positive or negative, is always expressed in reference to the same wild-type strain. The key result reported by Elena and Lenski, based on the standard tests, was that there were many significant epistatic deviations in both directions, in contrast to one hypothesis that predicted negative epistasis should be the general rule [17]. At first glance, the non-standard tests described above would not seem to allow this prediction to be tested, because the sign of the difference is arbitrary; that is, the labels that put one pair of mutants on one side of the equation, and not on the other side, are arbitrary. However, the biological interest in epistasis is often framed in terms of deleterious mutations, and Elena and Lenski ensured that they had a high-fitness wild-type strain by using one that had evolved for 10,000 generations in the exact same environment employed for measuring the relative fitness of all the genotypes in their study. In the same vein, we can anchor all 216 non-standard equations simply by ensuring that the genotype with the highest fitness (of the four used in any given calculation) appears on the left-hand side of the equation with positive sign.
To test this prediction for the fitness peak at the wild-type, we calculated for each standard test of epistasis (such as w·ar - a·r) the average fitness decrement Δ = 1 - (a + r)/2 of the two single mutants relative to the wild-type. We then plot the epistasis values as a function of Δ for all 27 equations. In this way, we can test the hypothesis that epistatic deviations tend more toward positive values when mutations in the wild-type background are more deleterious (large Δ) than when mutations are less deleterious (small Δ). The scope of this analysis can be extended substantially by also including the non-standard "double-double" epistatic equations (such as ar·bs - as·br). We anchored each double-double test with the fittest genotype (say ar) on the left-hand side, and calculated the average fitness decrement, Δ = ar - (as + br)/2, of the two genotypes that appear on the right-hand side of the equation. These tests ask more generally for the relationship between epistasis and the average fitness loss relative to any local fitness peak rather than only the one centred on the wild-type. Moreover, these non-standard tests may have greater power because their "reach" extends over mutational distances of four, rather than the two allowed by standard tests. (We have excluded the non-standard single-double tests from this analysis, because in those equations each genotype has different Hamming distances to the two genotypes on the opposite side of the difference equation. This heterogeneity prohibits direct comparisons with the other tests).
If we want to make the same type of inference about the complete genotype space, rather than the specific subset sampled by Elena and Lenski, we can apply a similar, but not identical, test. More precisely, we want to investigate the correlation between average fitness decrease and epistasis among any single and double mutants of the E. coli genome. In that case, the unit of independent observation becomes the nine single mutations, and the resulting inference concerns the landscape of all possible double mutations derived from the same parental strain. We again applied Tukey's jackknife, this time eliminating in turn all tests containing a particular mutation in any of the four double mutants, rather than eliminating a particular genotype as before. In this case, the association is only marginally significant (mean r = 0.572, t_{ s }= 1.363, 8 d.f., one-tailed p = 0.105 for the standard tests, and mean r = 0.457, t_{ s }= 1.780, 8 d.f., one-tailed p = 0.056 for the non-standard tests), but still supports the trend for increasingly positive epistasis with increasing average fitness loss. Given that these tests aim to infer a complex genome-wide relationship between average mutational effects and the form of epistasis from a small sample of mutations, it is encouraging to find such a strong trend.
We present these alternative analyses to emphasize the subtly different hypotheses that can be addressed by using our mathematical approach. Comparing the last two analyses suggests that individual mutations might have pervasive effects on the shape of the local landscape. While pervasive effects of certain mutations can make it more difficult to test broader generalizations, the precise nature of these pervasive effects is of biological interest. In the next section, we follow the same general approach, but focusing on a different set of epistatic interactions, to examine differences between individual mutations in greater detail.
Some non-standard tests reveal differences in mixing ability
In this section, we use non-standard tests of the "single-double" type to explore a particular aspect of the fitness landscape, specifically whether certain mutations are better mixers than others. The mixing ability of any particular mutation indicates whether its epistatic interactions with other mutations tend to be positive or negative. We can then measure the relative mixing ability of two mutations by holding constant the identity of other mutations with which the two of interest are mixed. Consider the polynomial a·br - b·ar. This test asks, in effect, whether mutation a or b mixes better with a third mutation r. An individual test of this form might be interesting when one has specific knowledge about the identity of the three mutated genes and the position of their products in a metabolic pathway, for example. By contrast, Elena and Lenski emphasized the statistical properties of epistatic interactions. In this context, we can examine related sets of these single-double equations to ask whether one mutation is a better mixer than another in the context of the sample of mutations with which they were each tested.
Any pair of mutations belonging to the same set of three ({a, b, c}, {r, s, t}, or {x, y, z}) was tested with the exact same set of six mutations belonging to the other two sets. For example, the following six equations examine the relative mixing ability of "focal mutations" a and b with respect to "tester mutations" r, s, t, x, y, and z: a·br - b·ar; a·bs - b·as; a·bt - b·at; a·bx - b·ax; a·by - b·ay; and a·bz - b·az. All in all, there are nine groups of six equations each that compare the general mixing abilities of two focal mutations from the same set: a versus b, a versus c, b versus c, r versus s, r versus t, s versus t, x versus y, x versus z, and y versus z. There are another 18 groups of three equations each that compare the mixing abilities of two focal mutations from different sets. (Nine more groups of this type are redundant, in the sense that their values are determined by the first 18 groups.)
Geometry of the fitness landscape
Our final set of results is concerned with the geometric shape of the fitness landscape. Since fitness landscapes are high-dimensional and complicated objects, it is desirable to classify them into a finite set of distinct shapes; with the general idea that fitness landscapes with the same shape are likely to share biological properties. This approach generalizes the classification of bi-allelic two-locus landscapes into those with positive epistasis versus those with negative epistasis. This appealing binary classification has been linked, for example, to the advantage of sex, but it does not extend to higher-dimensional genotype spaces. We present here a notion of the shape of a fitness landscapes for any genotypic space. This concept is intimately related to the interaction tests discussed so far, because the shape is determined by a certain subset of the gene interactions that includes the Markov basis. Thus, the proposed classification of landscapes into shapes can be regarded as a formal summary of all the various standard and non-standard tests.
The fitness landscape studied in this paper consists of the 37 E. coli genotypes and their fitness values as shown in Table 2. The genotope is the set of all possible allele frequencies that can be realized by any population on these genotypes. It is a nine-dimensional figure with 37 vertices, and it contains all the three-dimensional genotopes shown in Figure 1. By the shape of the fitness landscape we mean the triangulation of the genotope that is induced by the fitness values (see Methods for details). Since most of us have trouble visualizing and interpreting nine-dimensional objects, we therefore study the shape by analyzing its restrictions to genotypes on two and three of the nine loci. In so doing, we aim to illustrate how certain features of the fitness landscape - in particular, differences in mixing ability revealed by the non-standard tests - are reflected in the geometry of the fitness landscape.
Consider the bi-allelic two-locus system with genotypes w, a, r, and ar. Its genotope is the unit square in (a, r)-space, and a generic fitness landscape has exactly one of two possible shapes corresponding to either negative or positive epistasis, i.e., to w ar - a r being either negative or positive. Negative epistasis induces the triangulation of the square consisting of the two triangles {w, a, r} and {a, r, ar}, whereas positive epistasis results in the other possible triangulation with {w, a, ar} and {w, r, ar}.
For larger genetic systems, the role of the triangles is played by simplices. The shape of the present fitness landscape on 37 genotypes is a triangulation of the genotope into 362 nine-dimensional simplices (see Additional file 5). The following analysis of the geometry of this space is based on the general framework developed elsewhere [1], as applied to the specific experimental design that generated the E. coli fitness landscape investigated in this paper. The general idea is to investigate and describe the interaction space. This space has finite dimension, but infinitely many elements, and each element of this space represents one interaction. There are many different ways of extracting potentially interesting subsets of the interaction space. We suggest looking at the circuits, which generate the space but are not linearly independent. The signs of the circuits determine the triangulation of the genotope and thus the shape of the fitness landscape. However, the number of circuits grows fast with the size of the fitness landscape; the present E. coli landscape has 772,731 circuits. Therefore, we restrict our attention to the much smaller subset provided by the Markov basis. The circuits and the Markov basis provide a natural generalization of the concept of pairwise epistasis.
The analogous graph of possible shapes for the three-locus subsystem corresponding to the triangular prism (Figure 1c) is a hexagon. In fact, that hexagon occurs three times as a sub-graph in the graph for type (b) illustrated in Figure 5, because the genotope (c) is contained within the genotope (b) in three different orientations. Intuitively, any three-locus genotype space lacking a double and the triple mutant can be regarded as a subspace of the space that only lacks the triple mutant. The six triangulations of genotope (c) are represented in Figure 5 by the three hexagons with vertices {35, 36, 52, 21, 20, 50} (upper left), {36, 37, 56, 25, 23, 52} (upper right), and {37, 35, 50, 22, 24, 56} (bottom), respectively, that are related by symmetry.
The shape of fitness landscapes on type (c) spaces, such as {w, a, b, r, ar, br}, is determined by the signs of the three tests w·ar - a·r, w·br - b·r, and b·ar - a·br. Therefore, the shape summarizes the information about the standard pairwise interactions between mutations a and r and between b and r, as well as the relative mixing ability of a and b with respect to r. For example, five of the six subsystems of type (c) involving mutations y and z have shape 21, while the subspace {w, y, z, b, by, bz} has shape 52. Shape 21 is defined by a negative sign for the second of the three tests above and a positive sign for the other two, whereas shape 52 is defined by positive signs for all three tests. Hence, these shapes reflect positive epistasis involving y (6/6 tests), positive epistasis involving z (5/6 tests), and superior mixing ability of y over z (6/6 tests). As such, the geometric shapes provide more details about the form of mutational interactions than the average values of the tests analyzed in the previous section. On the other hand, these shapes summarize the data by classifying continuous fitness values into discrete shape classes which reflect the sign pattern of the interaction tests.
Discussion and Conclusion
Epistasis occurs whenever mutations interact non-linearly with one another, and it represents a major challenge in describing the mathematical structure of real fitness landscapes. With epistatic interactions, the combined effect of two or more mutations on fitness may be greater than, less than, or opposite in sign to expectations obtained by combining their separate effects. A growing body of empirical research indicates that epistasis is very common in nature [2, 11, 21, 23, 26–34], [39, 40]. However, a complete mathematical description of epistatic interactions has not been forthcoming for any system because the forms of epistasis appear to be diverse, idiosyncratic, and hence complex.
To date, two different aspects of epistasis have served as summary statistics of these interactions. First, studies have used the overall directional curvature of mean fitness as a function of the number of random mutations introduced into the genome of some wild-type organism [2, 27–31]. An older variant of this approach uses the time during which mutations have accumulated in a population subjected to severe bottlenecks as a proxy for estimating the number of mutations [41–44]. In any case, the absence of overall directional curvature does not distinguish between two biological scenarios: (i) most mutations have independent effects such that there is very little epistasis; and (ii) epistasis is common but interactions between mutations are diverse in their directional effects, thereby obscuring overall average curvature [2]. These two scenarios make different predictions about the evolution of a population on a fitness landscape that may confound efforts to understand, for example, the evolution of sexual reproduction [18, 20]. The second summary of epistasis considers the statistical distribution of a particular class of epistatic interactions, typically pairwise [45]. For example, Elena and Lenski [2] found that average fitness as a function of mutation number did not deviate significantly from log-linearity, which might suggest that epistasis is rare. But in a second experiment, they showed that significant pairwise interactions were common, although some were positive and others negative, such that there was no clear trend with respect to overall directional curvature.
The objective of this paper is to introduce biologists to a new mathematical framework for characterizing epistatic interactions between mutations, which goes beyond both overall directional curvature and pairwise interactions by providing a complete geometrical description of the epistatic interactions that define an empirically determined fitness landscape. To that end, we have re-analyzed the dataset from the pairwise experimental design performed by Elena and Lenski [2] using this new approach. In addition to providing an overall geometric description, various biologically motivated tests about the forms of epistasis are embedded in this framework. In particular, the geometric framework allows not only tests of the standard pairwise interactions but also non-pairwise tests that gave new insights into (i) the relationship between the form of epistasis and the individual mutational effects, and (ii) variation between mutations in their mixing ability with other mutations.
The fitness landscape that we analyzed comprises 37 genotypes of E. coli that were constructed by introducing nine mutations into a wild-type strain and constructing a restricted set of 27 double mutants. Despite the absence of any triple or other higher-order mutants in the dataset, our analysis reveals complex new epistatic interactions, beyond the pairwise interactions reported previously. First, our analysis confirms and extends Elena and Lenski's inference that the fitness landscape is complex, such that an overall measure of curvature obscures a complex admixture of interaction types, some with positive and others with negative effects on fitness (Figure 2). Second, we calculated the set of non-standard interactions that contrast two double mutants with two other double mutants, while holding all of the component alleles constant. In doing so, we found a strong correlation between the average fitness decrement associated with mutations and the resulting form of epistasis, such that epistatic deviations tend toward more positive fitness effects when the component mutations are more deleterious (Figure 3). This finding also led us to re-examine interactions based on the standard pairwise tests for evidence of this relationship, and the same trend was evident. This correlation is consistent with previous studies showing that compensatory mutations, which contribute to positive epistasis, become more important as one moves farther away from a local fitness peak [33–36]. More generally, this association emphasizes that any particular epistatic interaction is a local feature in the fitness landscape, and this association identifies one source of variation among local features.
Third, we show that individual mutations contribute in different ways to the complex admixture of epistatic interactions. In particular, we found that some mutations are better mixers than other mutations (Figure 4). That is, double mutants that include mutations that are relatively "good mixers" tend to be more fit than double mutants that harbour "bad mixers", even when the identity of their partner mutations is held constant. Although our analyses have not specifically addressed the evolution of sexual recombination, we suggest that appropriately designed tests of mixing ability may provide valuable insights into this area. A key difference between asexual and sexual reproduction is that the fate of a particular mutation is tied to its fitness effect in the genetic background where it arose in the case of asexual reproduction, while a mutation's fate in a sexual population depends on its effects over many backgrounds [46, 47]. Finally, we determine the overall geometric shape of the fitness landscape, which summarizes all the biologically interesting features described above (Figure 5).
In closing, we would like to raise an issue related to experimental design, one that requires attention when planning studies that might employ these new approaches to testing epistatic interactions and describing fitness landscapes. In their paper, Elena and Lenski [2] viewed it as problematic that the same mutations were used in multiple genotypes, because this compromised the statistical independence of some of their observations. They addressed this problem by applying Bonferroni corrections to their statistical tests, but they could instead have made 27 double-mutant genotypes with no overlap in the constituent mutations. Fortunately, they did not do so because, if they had, it would not have been possible to apply the mathematical approaches we have used to analyze the epistatic interactions and the shape of the underlying fitness landscape. Thus, minimizing the mutational overlap between genetic constructs may simplify statistical analyses, but it also constrains the analysis of more complex forms of epistasis. In particular, the existence of shared mutations allowed us to examine genotypes that spanned mutational distances of 3 (single-double tests) and 4 (double-double tests), which were essential for achieving the new inferences outlined in this study. Therefore, we recommend that future studies of epistatic interactions include many genotypes that share mutations in order to explore the geometry of the fitness landscape more fully. Of course, the inclusion of triple mutants and other higher-order genotypes can extend the reach of our geometric approach but, even then, the reach will be greater still if two or more sets of higher-order genotypes share some mutations.
Methods
Bacterial experiments
Our analyses use the data obtained from the second of two experiments reported by Elena and Lenski [2]. Details of genotype construction and competition assays are presented there. Briefly, they constructed 9 genotypes, all from the same initial strain, that each had a single mutation caused by the insertion of one of three Tn10 mini-transposons carrying different antibiotic-resistance markers. Other work indicated that the site of an insertion is largely responsible for a mutation's effect on fitness [48]. They then constructed 27 genotypes each carrying two of these mutations. Nine of the 36 possible double mutants were not constructed owing to restrictions based on having only three resistance markers, and no triple mutants were constructed in that experiment. The fitness of each of these 9 single mutants and 27 double mutants was measured in competition with the wild-type strain. Fitness measures the growth rate of a genotype realized during competition with the wild-type, and it is also expressed relative to the realized growth rate of the wild-type. Each pair-wise competition was replicated 10-fold and, in our analyses, we have used the median of the 10 estimates as the fitness value for each genotype. The fitness of the wild-type strain was set to 1.
Mathematical framework
Our analysis is based on the mathematical framework presented by Beerenwinkel et al. [1]. The dataset analyzed here induces the genotype space G ⊆ {0,1}^{9} of size 37 consisting of the wild-type strain w = 000000000, the nine single mutants a = 100000000, ..., z = 000000001, and the subset of 27 double mutants ar = 100100000, ..., tz = 000001001 (Table 1). We denote by Δ_{ G }⊂ R^{37} the set of all populations (genotype frequencies) on G. The genotope Π_{ G }⊂ R^{9} is the set of all allele frequencies that can be realized by populations on G. The genotope is a nine-dimensional polytope, defined as the convex hull of the 37 genotypes, i.e., as the set of all convex combinations $\sum _{g\in G}{\lambda}_{g}\cdot g$, with λ_{ g }≥ 0 and ${\sum}_{g\in G}{\lambda}_{g}}=1$). Figure 1 shows the genotopes of all three-locus subsystems induced by G, all of which are contained in Π_{ G }.
There is a natural mapping ρ : Δ_{ G }→ Π_{ G }that assigns to each population its allele frequency spectrum. The kernel of this map defines the interaction space (see Section 3 in [1]). The elements of the interaction space measure the gene interactions that underlie the given fitness landscape on G. They include both the standard and non-standard tests used in our analyses.
The specific experimental design that was used to generate the present fitness dataset allows us to represent the genotypes in the matrix displayed in Table 1. The entries of this table are interpreted as variables. If we regard Table 1, hypothetically, as a contingency table of genotype counts, then taking row sums (or, equivalently, column sums) corresponds precisely to evaluating the mapping ρ. The kernel of ρ is the space of integer tables of the same format with zero margins (row sums). A canonical set of generators for this kernel is the minimal Markov basis [49], which consists of polynomials in the unknown variables, one for each genotype. These polynomials correspond to the minimal generators of the toric ideal associated with ρ (see [51] for an introduction to toric ideals). The Markov basis represents the fundamental gene interactions that underlie the given fitness landscape. It consists of the 27 standard tests and the 216 non-standard tests that are described in the Results section and listed in supplementary material ( see Additional file 4).
A simplex is an n-dimensional analogue of a triangle. Formally, a simplex can be defined as the convex hull of a set of (n + 1) affinely independent points in R^{ n }[50]. This means that the (n + 1) points span an affine space of dimension n. For example, Figure 1d shows a 3-simplex (known as a tetrahedron), which is the convex hull of the points w, a, b, and c. By a triangulation of a polytope we mean its decomposition into a set of simplices. A fitness landscape on G gives rise to a triangulation of the genotope Π_{ G }(see Section 4 in [1]). Thus, the shapes of fitness landscapes on a genotype space G are defined as the induced triangulations of its genotope. These triangulations are determined by the initial ideal of the toric ideal associated with ρ. The derivation of the triangulation from the initial ideal is explained in [[51], Chapter 8]. The specific triangulation induced by the data in Table 2 consists of 362 simplices and is provided as supplementary material (see Additional file 5). In Figure 5, all possible shapes are displayed for the three-locus system that corresponds to the genotope shown in Figure 1b. In fact, the shapes correspond to the vertices of the three-dimensional secondary polytope of this genotope.
Computational methods
The Markov basis can always be derived from the genotype space (the set of measured genotypes) by algebraic computations. However, there is no simple recipe for writing out the Markov basis. We computed the Markov basis independently with the computer algebra systems Macaulay 2 [52] and Singular [53]. The triangulation was computed independently with Macaulay 2 and 4ti2 [54]. All statistical computations were performed in R [55]. In our supplementary materials, we illustrate the use of the Macaulay 2 program for calculating the Markov basis and triangulation of a dataset [see Additional files 1-3].
Declarations
Acknowledgements
This work was supported by the "FunBio" grant from DARPA to Simon Levin (Princeton University). We thank Simon Levin and Ben Mann (DARPA) for facilitating this math-bio collaboration, Peter Bates and Charles Ofria for helpful discussions, Peter Malkin for an improved program for computing circuits, Mike Stillman for writing and improving the Macaulay 2 code, and three anonymous reviewers for suggestions. Collection of the dataset used in this study was funded by a fellowship from the Spanish MEC to S.F.E. and a grant from the NSF to R.E.L. N.B. was funded by a grant from the Bill & Melinda Gates foundation through the Grand Challenges in Global Health Initiative.
Authors’ Affiliations
References
- Beerenwinkel N, Pachter L, Sturmfels B: Epistasis and shapes of fitness landscapes. Statistica Sinica.Google Scholar
- Elena SF, Lenski RE: Test of synergistic interactions among deleterious mutations in bacteria. Nature. 1997, 390: 395-398. 10.1038/37108.View ArticlePubMedGoogle Scholar
- Wright S: The roles of mutation, inbreeding, crossbreeding, and selection in evolution. Proc Sixth Intl Cong Genetics. 1932, 1: 356-366.Google Scholar
- Fisher RA: The Genetical Theory of Natural Selection. 1930, Oxford University PressView ArticleGoogle Scholar
- Simpson GG: The Major Features of Evolution. 1953, Columbia University PressGoogle Scholar
- Provine WB: Sewall Wright and Evolutionary Biology. 1986, University of Chicago PressGoogle Scholar
- Kauffman SA, Levin S: Towards a general theory of adaptive walks on rugged landscapes. J Theor Biol. 1987, 128: 11-45.View ArticlePubMedGoogle Scholar
- Gavrilets S: Fitness Landscapes and the Origin of Species. 2004, Princeton University PressGoogle Scholar
- Wolf JB, Brodie ED, Wade MJ MJ, editors: Epistasis and the Evolutionary Process. 2000, Oxford University PressGoogle Scholar
- Lenski RE, Rose MR, Simpson SC, Tadler SC: Long-term experimental evolution in Escherichia coli I. Adaptation and divergence during 2,000 generations. Am Nat. 1991, 138: 1315-1341. 10.1086/285289.View ArticleGoogle Scholar
- Cooper TF, Rozen DE, Lenski RE: Parallel changes in gene expression after 20,000 generations of evolution in E. coli. Proc Natl Acad Sci USA. 2003, 100: 1072-1077. 10.1073/pnas.0334340100.PubMed CentralView ArticlePubMedGoogle Scholar
- Wade MJ: Epistasis as a genetic constraint within populations and an accelerant of adaptive divergence among them. Epistasis and the Evolutionary Process. Edited by: Wolf JB, Brodie ED, Wade MJ. 2000, Oxford University Press, 213-231.Google Scholar
- Dobzhansky T: Genetics and the Origin of Species. 1937, Columbia University PressGoogle Scholar
- Orr HA: The population genetics of speciation: the evolution of hybrid incompatibilities. Genetics. 1995, 139: 1805-1813.PubMed CentralPubMedGoogle Scholar
- Turelli M, Orr HA: Dominance, epistasis and the genetics of postzygotic isolation. Genetics. 2000, 154: 1663-1679.PubMed CentralPubMedGoogle Scholar
- Demuth JP, Wade MJ: On the theoretical and empirical framework for studying genetic interactions within and among species. Am Nat. 2005, 165: 524-536. 10.1086/429276.View ArticlePubMedGoogle Scholar
- Kondrashov AS: Classification of hypotheses on the advantage of amphimixis. J Heredity. 1993, 84: 372-387.Google Scholar
- Otto SP, Feldman MW: Deleterious mutations, variable epistatic interactions, and the evolution of recombination. Theor Pop Biol. 1997, 51: 134-147. 10.1006/tpbi.1997.1301.View ArticleGoogle Scholar
- Peters AD, Lively CM: Epistasis and the maintenance of sex. Epistasis and the Evolutionary Process. Edited by: Wolf JB, Brodie ED, Wade MJ. 2000, Oxford University Press, 99-112.Google Scholar
- Kouyos RD, Otto SP, Bonhoeffer S: Effect of varying epistasis on the evolution of recombination. Genetics. 2006, 173: 589-597. 10.1534/genetics.105.053108.PubMed CentralView ArticlePubMedGoogle Scholar
- Elena SF, Lenski RE: Epistasis between new mutations and genetic background and a test of genetic canalization. Evolution. 2001, 55: 1746-1752. 10.1554/0014-3820(2001)055[1746:EBNMAG]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Wilke CO, Wang J, Ofria C, Lenski RE, Adami C: Evolution of digital organisms at high mutation rate leads to survival of the flattest. Nature. 2001, 412: 331-333. 10.1038/35085569.View ArticlePubMedGoogle Scholar
- Burch CL, Chao L: Epistasis and its relationship to canalization in the RNA virus ϕ6. Genetics. 2004, 167: 559-567. 10.1534/genetics.103.021196.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice SH: The evolution of developmental interactions: epistasis, canalization, and integration. Epistasis and the Evolutionary Process. Edited by: Wolf JB, Brodie ED, Wade MJ. 2000, Oxford University Press, 82-98.Google Scholar
- Lenski RE: The cost of antibiotic resistance – from the perspective of a bacterium. Ciba Foundation Symposia. 1997, 207: 131-140.PubMedGoogle Scholar
- Bonhoeffer S, Chappey C, Parkin NT, Whitcomb JM, Petropoulos CJ: Evidence for positive epistasis in HIV-1. Science. 2004, 306: 1547-1550. 10.1126/science.1101786.View ArticlePubMedGoogle Scholar
- de Visser JAGM, Hoekstra RF, van den Ende H: Test of interaction between genetic markers that affect fitness in Aspergillus niger. Evolution. 1997, 51: 1499-1505. 10.2307/2411202.View ArticleGoogle Scholar
- Whitlock MC, Bourguet D: Factors affecting the genetic load in Drosophila : synergistic epistasis and correlations among fitness components. Evolution. 2000, 54: 1654-1660. 10.1554/0014-3820(2000)054[1654:FATGLI]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Szafraniec K, Wloch DM, Sliwa P, Borts RH, Korona R: Small fitness effects and weak genetic interactions between deleterious mutations in heterozygous loci of the yeast Saccharomyces cerevisiae. Genet Res. 2003, 82: 19-31. 10.1017/S001667230300630X.View ArticlePubMedGoogle Scholar
- Sanjuán R, Moya A, Elena SF: The contribution of epistasis to the architecture of fitness in an RNA virus. Proc Natl Acad Sci USA. 2004, 101: 15376-15379. 10.1073/pnas.0404125101.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanjuán R, Elena SF: Epistasis correlates to genomic complexity. Proc Natl Acad Sci USA. 2006, 103: 14402-14405. 10.1073/pnas.0604543103.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore FBG, Rozen DE, Lenski RE: Pervasive compensatory adaptation in Escherichia coli. Proc Roy Soc, London B. 2006Google Scholar
- Wilke CO, Lenski RE, Adami C: Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding. BMC Evol Biol. 2003, 3: 3-10.1186/1471-2148-3-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Poon A, Chao L: The rate of compensatory mutation in the DNA bacteriophage φX174. Genetics. 2005, 170: 989-999. 10.1534/genetics.104.039438.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilke CO, Adami C: Interaction between directional epistasis and average mutational effects. Proc Roy Soc, London B. 2001, 268: 1469-1474. 10.1098/rspb.2001.1690.View ArticleGoogle Scholar
- Sanjuán R, Forment J, Elena SF: In silico predicted robustness of viroid RNA structure. II. Interaction between mutation pairs. Mol Biol Evol. 2006, 23: 2123-2130. 10.1093/molbev/msl083.View ArticlePubMedGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry. 1981, New York, FreemanGoogle Scholar
- Lunzer M, Miller SP, Felsheim R, Dean AM: The biochemical architecture of an ancient adaptive landscape. Science. 2005, 310: 499-501. 10.1126/science.1115649.View ArticlePubMedGoogle Scholar
- 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. 10.1126/science.1123539.View ArticlePubMedGoogle Scholar
- Poelwijk FJ, Kiviet DJ, Tans SJ: Evolutionary potential of a duplicated repressor-operator pair: simulating pathways using mutation data. PLoS Comput Biol. 2006, 2: e58-10.1371/journal.pcbi.0020058.PubMed CentralView ArticlePubMedGoogle Scholar
- Mukai T: Genetic structure of natural populations of Drosophila melanogaster. VII. Synergistic interaction of spontaneous mutant polygenes controlling viability. Genetics. 1969, 61: 749-761.PubMed CentralPubMedGoogle Scholar
- Chao L: Evolution of sex in RNA viruses. J Theor Biol. 1988, 133: 99-112. 10.1016/S0022-5193(88)80027-4.View ArticlePubMedGoogle Scholar
- Kibota TT, Lynch M: Estimate of the genomic mutation rate deleterious to overall fitness in E. coli . Nature. 1996, 381: 694-696. 10.1038/381694a0.View ArticlePubMedGoogle Scholar
- Shaw FH, Geyer CJ, Shaw RG: A comprehensive model of mutations affecting fitness and inferences for Arabidopsis thaliana. Evolution. 2002, 56: 453-463. 10.1554/0014-3820(2002)056[0453:ACMOMA]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Phillips PC, Otto SP, Whitlock MC: Beyond the average: The evolutionary importance of epistasis and the variability of epistatic effects. Epistasis and the Evolutionary Process. Edited by: Oxford University Press. 2000, Wolf JB, Brodie ED, Wade MJ, 20-38.Google Scholar
- Malmberg RL: Evolution of epistasis and advantage of recombination in populations of bacteriophage T4. Genetics. 1977, 86: 607-621.PubMed CentralPubMedGoogle Scholar
- Misevic D, Ofria C, Lenski RE: Sexual reproduction reshapes the genetic architecture of digital organisms. Proc Roy Soc, London B. 2006, 273: 457-464. 10.1098/rspb.2005.3338.View ArticleGoogle Scholar
- Elena SF, Ekunwe L, Hajela N, Oden SA, Lenski RE: Distribution of fitness effects caused by random insertion mutations in Escherichia coli. Genetica. 1998, 102/103: 349-358. 10.1023/A:1017031008316.View ArticleGoogle Scholar
- Diaconis P, Sturmfels B: Algebraic algorithms for sampling from conditional distributions. Ann Stat. 1998, 26: 363-397. 10.1214/aos/1030563990.View ArticleGoogle Scholar
- Ziegler G: Lectures on polytopes. 1995, New York, SpringerView ArticleGoogle Scholar
- Sturmfels B: Gröbner bases and convex polytopes. Providence, RI, American Mathematical Society. 1996Google Scholar
- Grayson DR, Stillman ME: Macaulay 2, a software system for research in algebraic geometry. [http://www.math.uiuc.edu/Macaulay2/]
- Greuel G-M, Pfister G, Schönemann H: Singular 3.0. A computer algebra system for polynomial computations. [http://www.singular.uni-kl.de]
- Hemmecke R, Hemmecke R, Malkin P: 4ti2 Version 1.2–Computation of Hilbert bases, Graver bases, toric Gröbner bases, and more. [http://www.4ti2.de]
- R Development Core Team: R: a language and environment for statistical computing. Vienna, Austria, R Foundation for Statistical Computing. 2006, [http://www.R-project.org]Google Scholar
Copyright
This article is published under license to 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.