Epistasis from functional dependence of fitness on underlying traits

Mutational robustness changes during long-term adaptation in laboratory budding yeast populations

As an adapting population traverses the fitness landscape, its local neighborhood (i.e., the collection of fitness effects of single-step mutations) can change shape because of interactions with mutations acquired during evolution. These changes to the distribution of fitness effects can affect both the rate of adaptation and the accumulation of deleterious mutations. However, while numerous models of fitness landscapes have been proposed in the literature, empirical data on how this distribution changes during evolution remains limited. In this study, we directly measure how the fitness landscape neighborhood changes during laboratory adaptation. Using a barcode-based mutagenesis system, we measure the fitness effects of 91 specific gene disruption mutations in genetic backgrounds spanning 8000–10,000 generations of evolution in two constant environments. We find that the mean of the distribution of fitness effects decreases in one environment, indicating a reduction in mutational robustness, but does not change in the other. We show that these distribution-level patterns result from differences in the relative frequency of certain patterns of epistasis at the level of individual mutations, including fitness-correlated and idiosyncratic epistasis.

Emergence and propagation of epistasis in metabolic networks

Epistasis is often used to probe functional relationships between genes, and it plays an important role in evolution. However, we lack theory to understand how functional relationships at the molecular level translate into epistasis at the level of whole-organism phenotypes, such as fitness. Here, I derive two rules for how epistasis between mutations with small effects propagates from lower- to higher-level phenotypes in a hierarchical metabolic network with first-order kinetics and how such epistasis depends on topology. Most importantly, weak epistasis at a lower level may be distorted as it propagates to higher levels. Computational analyses show that epistasis in more realistic models likely follows similar, albeit more complex, patterns. These results suggest that pairwise inter-gene epistasis should be common, and it should generically depend on the genetic background and environment. Furthermore, the epistasis coefficients measured for high-level phenotypes may not be sufficient to fully infer the underlying functional relationships.

Accelerated inbreeding depression suggests synergistic epistasis for deleterious mutations in Drosophila melanogaster

Epistasis may have important consequences for a number of issues in quantitative genetics and evolutionary biology. In particular, synergistic epistasis for deleterious alleles is relevant to the mutation load paradox and the evolution of sex and recombination. Some studies have shown evidence of synergistic epistasis for spontaneous or induced deleterious mutations appearing in mutation-accumulation experiments. However, many newly arising mutations may not actually be segregating in natural populations because of the erasing action of natural selection. A demonstration of synergistic epistasis for naturally segregating alleles can be achieved by means of inbreeding depression studies, as deleterious recessive allelic effects are exposed in inbred lines. Nevertheless, evidence of epistasis from these studies is scarce and controversial. In this paper, we report the results of two independent inbreeding experiments carried out with two different populations of Drosophila melanogaster. The results show a consistent accelerated inbreeding depression for fitness, suggesting synergistic epistasis among deleterious alleles. We also performed computer simulations assuming different possible models of epistasis and mutational parameters for fitness, finding some of them to be compatible with the results observed. Our results suggest that synergistic epistasis for deleterious mutations not only occurs among newly arisen spontaneous or induced mutations, but also among segregating alleles in natural populations.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

Power law fitness landscapes and their ability to predict fitness

Whether or not evolution by natural selection is predictable depends on the existence of general patterns shaping the way mutations interact with the genetic background. This interaction, also known as epistasis, has been observed during adaptation (macroscopic epistasis) and in individual mutations (microscopic epistasis). Interestingly, a consistent negative correlation between the fitness effect of beneficial mutations and background fitness (known as diminishing returns epistasis) has been observed across different species and conditions. We tested whether the adaptation pattern of an additional species, Schizosaccharomyces pombe, followed the same trend. We used strains that differed by the presence of large karyotype differences and observed the same pattern of fitness convergence. Using these data along with published datasets, we measured the ability of different models to describe adaptation rates. We found that a phenotype-fitness landscape shaped like a power law is able to correctly predict adaptation dynamics in a variety of species and conditions. Furthermore we show that this model can provide a link between the observed macroscopic and microscopic epistasis. It may be very useful in the development of algorithms able to predict the adaptation of microorganisms from measures of the current phenotypes. Overall, our results suggest that even though adaptation quickly slows down, populations adapting to lab conditions may be quite far from a fitness peak.

Negative selection in humans and fruit flies involves synergistic epistasis

Negative selection against deleterious alleles produced by mutation influences within-population variation as the most pervasive form of natural selection. However, it is not known whether deleterious alleles affect fitness independently, so that cumulative fitness loss depends exponentially on the number of deleterious alleles, or synergistically, so that each additional deleterious allele results in a larger decrease in relative fitness. Negative selection with synergistic epistasis should produce negative linkage disequilibrium between deleterious alleles and, therefore, an underdispersed distribution of the number of deleterious alleles in the genome. Indeed, we detected underdispersion of the number of rare loss-of-function alleles in eight independent data sets from human and fly populations. Thus, selection against rare protein-disrupting alleles is characterized by synergistic epistasis, which may explain how human and fly populations persist despite high genomic mutation rates.

Additive Phenotypes Underlie Epistasis of Fitness Effects

Gene interactions, or epistasis, play a large role in determining evolutionary outcomes. The ruggedness of fitness landscapes, and thus the predictability of evolution and the accessibility of high-fitness genotypes, is determined largely by the pervasiveness of epistasis and the degree of correlation between similar genotypes. We created all possible pairings of three sets of five beneficial first-step mutations fixed during adaptive walks under three different regimes: selection on growth rate alone, on growth rate and thermal stability, and on growth rate and pH stability. All 30 double-mutants displayed negative, antagonistic epistasis with regard to growth rate and fitness, but positive epistasis and additivity were common for the stability phenotypes. This suggested that biophysically simple phenotypes, such as capsid stability, may on average behave more additively than complex phenotypes like viral growth rate. Growth rate epistasis was also smaller in magnitude when the individual effects of single mutations were smaller. Significant sign epistasis, such that the effect of a mutation that is beneficial in the wild-type background is deleterious in combination with a second mutation, emerged more frequently in intragenic mutational pairings than in intergenic pairs, and was evident in nearly half of the double-mutants, indicating that the fitness landscape is moderately uncorrelated and of intermediate ruggedness. Together, our results indicated that mutations may interact additively with regard to phenotype when considered at a basic, biophysical level, but that epistasis arises as a result of pleiotropic interactions between the individual components of complex phenotypes and diminishing returns arising from intermediate phenotypic optima.

What can we learn from fitness landscapes?

A combinatorially complete data set consists of studies of all possible combinations of a set of mutant sites in a gene or mutant alleles in a genome. Among the most robust conclusions from these studies is that epistasis between beneficial mutations often shows a pattern of diminishing returns, in which favorable mutations are less fit when combined than would be expected. Another robust inference is that the number of adaptive evolutionary paths is often limited to a relatively small fraction of the theoretical possibilities, owing largely to sign epistasis requiring evolutionary steps that would entail a decrease in fitness. Here we summarize these and other results while also examining issues that remain unresolved and future directions that seem promising.

Evolutionary constraints in variable environments, from proteins to networks

Environmental changes can not only trigger a regulatory response, but also impose evolutionary pressures that can modify the underlying regulatory network. Here, we review recent approaches that are beginning to disentangle this complex interplay between regulatory and evolutionary responses. Systematic genetic reconstructions have shown how evolutionary constraints arise from epistatic interactions between mutations in fixed environments. This approach is now being extended to more complex environments and systems. The first results suggest that epistasis is affected dramatically by environmental changes and, hence, can profoundly affect the course of evolution. Thus, external environments not only define the selection of favored phenotypes, but also affect the internal constraints that can limit the evolution of these phenotypes. These findings also raise new questions relating to the conditions for evolutionary transitions and the evolutionary potential of regulatory networks.

Mapping the fitness landscape of gene expression uncovers the cause of antagonism and sign epistasis between adaptive mutations

How do adapting populations navigate the tensions between the costs of gene expression and the benefits of gene products to optimize the levels of many genes at once? Here we combined independently-arising beneficial mutations that altered enzyme levels in the central metabolism of Methylobacterium extorquens to uncover the fitness landscape defined by gene expression levels. We found strong antagonism and sign epistasis between these beneficial mutations. Mutations with the largest individual benefit interacted the most antagonistically with other mutations, a trend we also uncovered through analyses of datasets from other model systems. However, these beneficial mutations interacted multiplicatively (i.e., no epistasis) at the level of enzyme expression. By generating a model that predicts fitness from enzyme levels we could explain the observed sign epistasis as a result of overshooting the optimum defined by a balance between enzyme catalysis benefits and fitness costs. Knowledge of the phenotypic landscape also illuminated that, although the fitness peak was phenotypically far from the ancestral state, it was not genetically distant. Single beneficial mutations jumped straight toward the global optimum rather than being constrained to change the expression phenotypes in the correlated fashion expected by the genetic architecture. Given that adaptation in nature often results from optimizing gene expression, these conclusions can be widely applicable to other organisms and selective conditions. Poor interactions between individually beneficial alleles affecting gene expression may thus compromise the benefit of sex during adaptation and promote genetic differentiation.

The pace and outcome of a series of adaptive steps in an evolving lineage depends upon how well different beneficial mutations stack on top of each other. We found that independent beneficial mutations that affected gene expression for a metabolic pathway did not work well together, and were often jointly deleterious. The most beneficial mutations interacted the most poorly with others, which was a trend we found common in other biological systems. Through generating a model that accounted for enzymatic benefits and expression costs, we uncovered that this antagonism was caused by a phenotype to fitness mapping that had an intermediate peak. This allowed us to predict the fitness effect of double mutants and to uncover that the single winning mutations tended to move straight to the peak in a single step. These findings demonstrate the importance of considering the phenotypic changes that cause nonlinear interactions between mutations upon fitness, and thus influence how populations evolve.

Invariance and optimality in the regulation of an enzyme

BACKGROUND

The Michaelis-Menten equation, proposed a century ago, describes the kinetics of enzyme-catalyzed biochemical reactions. Since then, this equation has been used in countless, increasingly complex models of cellular metabolism, often including time-dependent enzyme levels. However, even for a single reaction, there remains a fundamental disconnect between our understanding of the reaction kinetics, and the regulation of that reaction through changes in the abundance of active enzyme.

RESULTS

We revisit the Michaelis-Menten equation under the assumption of a time-dependent enzyme concentration. We show that all temporal enzyme profiles with the same average enzyme level yield identical substrate degradation- a simple analytical conclusion that can be thought of as an invariance principle, and which we validate experimentally using a β-galactosidase assay. The ensemble of all time-dependent enzyme trajectories with the same average concentration constitutes a space of functions. We develop a simple model of biological fitness which assigns a cost to each of these trajectories (in the form of a function of functions, i.e. a functional). We then show how one can use variational calculus to analytically infer temporal enzyme profiles that minimize the overall enzyme cost. In particular, by separately treating the static costs of amino acid sequestration and the dynamic costs of protein production, we identify a fundamental cellular tradeoff.

CONCLUSIONS

The overall metabolic outcome of a reaction described by Michaelis-Menten kinetics is ultimately determined by the average concentration of the enzyme during a given time interval. This invariance in analogy to path-independent phenomena in physics, suggests a new way in which variational calculus can be employed to address biological questions. Together, our results point to possible avenues for a unified approach to studying metabolism and its regulation.

REVIEWERS

This article was reviewed by Sergei Maslov, William Hlavacek and Daniel Kahn.