Long-term adaptation of epistatic genetic networks

Gene networks are likely to govern most traits in nature. Mutations at these genes often show functional epistatic interactions that lead to complex genetic architectures and variable fitness effects in different genetic backgrounds. Understanding how epistatic genetic systems evolve in nature remains one of the great challenges in evolutionary biology. Here we combine an analytical framework with individual-based simulations to generate novel predictions about long-term adaptation of epistatic networks. We find that relative to traits governed by independently evolving genes, adaptation with epistatic gene networks is often characterized by longer waiting times to selective sweeps, lower standing genetic variation, and larger fitness effects of adaptive mutations. This may cause epistatic networks to either adapt more slowly or more quickly relative to a nonepistatic system. Interestingly, epistatic networks may adapt faster even when epistatic effects of mutations are on average deleterious. Further, we study the evolution of epistatic properties of adaptive mutations in gene networks. Our results show that adaptive mutations with small fitness effects typically evolve positive synergistic interactions, whereas adaptive mutations with large fitness effects evolve positive synergistic and negative antagonistic interactions at approximately equal frequencies. These results provide testable predictions for adaptation of traits governed by epistatic networks and the evolution of epistasis within networks.     

SELECTION BIASES THE PREVALENCE AND TYPE OF EPISTASIS ALONG ADAPTIVE TRAJECTORIES

The contribution to an organism's phenotype from one genetic locus may depend upon the status of other loci. Such epistatic interactions among loci are now recognized as fundamental to shaping the process of adaptation in evolving populations. Although little is known about the structure of epistasis in most organisms, recent experiments with bacterial populations have concluded that antagonistic interactions abound and tend to deaccelerate the pace of adaptation over time. Here, we use the NK model of fitness landscapes to examine how natural selection biases the mutations that substitute during evolution based on their epistatic interactions. We find that, even when beneficial mutations are rare, these biases are strong and change substantially throughout the course of adaptation. In particular, epistasis is less prevalent than the neutral expectation early in adaptation and much more prevalent later, with a concomitant shift from predominantly antagonistic interactions early in adaptation to synergistic and sign epistasis later in adaptation. We observe the same patterns when reanalyzing data from a recent microbial evolution experiment. These results show that when the order of substitutions is not known, standard methods of analysis may suggest that epistasis retards adaptation when in fact it accelerates it.     

The Consistency of Beneficial Fitness Effects of Mutations across Diverse Genetic Backgrounds

Parallel and convergent evolution have been remarkably common observations in molecular adaptation but primarily in the context of the same genotype adapting to the same conditions. These phenomena therefore tell us about the stochasticity and limitations of adaptation. The limited data on convergence and parallelism in the adaptation of different genotypes conflict as to the importance of such events. If the effects of beneficial mutations are highly context dependent (i.e., if they are epistatic), different genotypes should adapt through different mutations. Epistasis for beneficial mutations has been investigated but mainly through measurement of interactions between individually beneficial mutations for the same genotype. We examine epistasis for beneficial mutations at a broader genetic scale by measuring the fitness effects of two mutations beneficial for the ssDNA bacteriophage ID11 in eight different, related genotypes showing 0.3-3.7% nucleotide divergence from ID11. We found no evidence for sign epistasis, but the mutations tended to have much smaller or no effects on fitness in the new genotypes. We found evidence for diminishing-returns epistasis; the effects were more beneficial for lower-fitness genotypes. The patterns of epistasis were not determined by phylogenetic relationships to the original genotype. To improve our understanding of the patterns of epistasis, we fit the data to a model in which each mutation had a constant, nonepistatic phenotypic effect across genotypes and the phenotype-fitness map had a single optimum. This model fit the data well, suggesting that epistasis for these mutations was due to nonlinearity in the phenotype-fitness mapping and that the likelihood of parallel evolution depends more on phenotype than on genotype.     

Recombination Accelerates Adaptation on a Large-Scale Empirical Fitness Landscape in HIV-1

Recombination has the potential to facilitate adaptation. In spite of the substantial body of theory on the impact of recombination on the evolutionary dynamics of adapting populations, empirical evidence to test these theories is still scarce. We examined the effect of recombination on adaptation on a large-scale empirical fitness landscape in HIV-1 based on in vitro fitness measurements. Our results indicate that recombination substantially increases the rate of adaptation under a wide range of parameter values for population size, mutation rate and recombination rate. The accelerating effect of recombination is stronger for intermediate mutation rates but increases in a monotonic way with the recombination rates and population sizes that we examined. We also found that both fitness effects of individual mutations and epistatic fitness interactions cause recombination to accelerate adaptation. The estimated epistasis in the adapting populations is significantly negative. Our results highlight the importance of recombination in the evolution of HIV-I.     

Environment Determines Epistatic Patterns for a ssDNA Virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype–fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype–fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype–phenotype relationship.     

Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors.     

Exploring the complexity of the HIV-1 fitness landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.     

Pervasive,yet idiosyncratic,epistatic pleiotropy during adaptation in a behaviourally complex microbe

Understanding how multiple mutations interact to jointly impact multiple ecologically important traits is critical for creating a robust picture of organismal fitness and the process of adaptation. However, this is complicated by both environmental heterogeneity and the complexity of genotype‐to‐phenotype relationships generated by pleiotropy and epistasis. Moreover, little is known about how pleiotropic and epistatic relationships themselves change over evolutionary time. The soil bacterium Myxococcus xanthus employs several distinct social traits across a range of environments. Here, we use an experimental lineage of M. xanthus that evolved a novel form of social motility to address how interactions between epistasis and pleiotropy evolve. Specifically, we test how mutations accumulated during selection on soft agar pleiotropically affect several other social traits (hard agar motility, predation and spore production). Relationships between changes in swarming rate in the selective environment and the four other traits varied greatly over time in both direction and magnitude, both across timescales of the entire evolutionary lineage and individual evolutionary time steps. We also tested how a previously defined epistatic interaction is pleiotropically expressed across these traits. We found that phenotypic effects of this epistatic interaction were highly correlated between soft and hard agar motility, but were uncorrelated between soft agar motility and predation, and inversely correlated between soft agar motility and spore production. Our results show that ‘epistatic pleiotropy’ varied greatly in magnitude, and often even in sign, across traits and over time, highlighting the necessity of simultaneously considering the interacting complexities of pleiotropy and epistasis when studying the process of adaptation.     

The impact of high‐order epistasis in the within‐host fitness of a positive‐sense plant RNA virus

RNA viruses are the main source of emerging infectious diseases because of the evolutionary potential bestowed by their fast replication, large population sizes and high mutation and recombination rates. However, an equally important property, which is usually neglected, is the topography of the fitness landscape. How many fitness maxima exist and how well they are connected is especially interesting, as this determines the number of accessible evolutionary pathways. To address this question, we have reconstructed a region of the fitness landscape of tobacco etch potyvirus constituted by mutations observed during the experimental adaptation of the virus to the novel host Arabidopsis thaliana. Fitness was measured for many genotypes and showed the existence of multiple peaks and holes in the landscape. We found prevailing epistatic effects between mutations, with cases of reciprocal sign epistasis being common among pairs of mutations. We also found that high‐order epistasis was as important as pairwise epistasis in their contribution to fitness. Therefore, results suggest that the landscape was rugged due to the existence of holes caused by lethal genotypes, that a very limited number of potential neutral paths exist and that it contained a single adaptive peak.     

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41-215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.     

Population Size Affects Adaptation in Complex Ways: Simulations on Empirical Adaptive Landscapes

Do large populations always outcompete smaller ones? Does increasing the mutation rate have a similar effect to increasing the population size, with respect to the adaptation of a population? How important are substitutions in determining the adaptation rate? In this study, we ask how population size and mutation rate interact to affect adaptation on empirical adaptive landscapes. Using such landscapes, we do not need to make many ad hoc assumption about landscape topography, such as about epistatic interactions among mutations or about the distribution of fitness effects. Moreover, we have a better understanding of all the mutations that occur in a population and their effects on the average fitness of the population than we can know in experimental studies. Our results show that the evolutionary dynamics of a population cannot be fully explained by the population mutation rate \(N\mu\); even at constant \(N\mu\), there can be dramatic differences in the adaptation of populations of different sizes. Moreover, the substitution rate of mutations is not always equivalent to the adaptation rate, because we observed populations adapting to high adaptive peaks without fixing any mutations. Finally, in contrast to some theoretical predictions, even on the most rugged landscapes we study, small population size is never an advantage over larger population size. These result show that complex interactions among multiple factors can affect the evolutionary dynamics of populations, and simple models should be taken with caution.     

The strength of genetic interactions scales weakly with mutational effects

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.     

Perspective: Sign epistasis and genetic constraint on evolutionary trajectories

Epistasis for fitness means that the selective effect of a mutation is conditional on the genetic background in which it appears. Although epistasis is widely observed in nature, our understanding of its consequences for evolution by natural selection remains incomplete. In particular, much attention focuses only on its influence on the instantaneous rate of changes in frequency of selected alleles via epistatic contribution to the additive genetic variance for fitness. Thus, in this framework epistasis only has evolutionary importance if the interacting loci are simultaneously segregating in the population. However, the selective accessibility of mutational trajectories to high fitness genotypes may depend on the genetic background in which novel mutations appear, and this effect is independent of population polymorphism at other loci. Here we explore this second influence of epistasis on evolution by natural selection. We show that it is the consequence of a particular form of epistasis, which we designate sign epistasis. Sign epistasis means that the sign of the fitness effect of a mutation is under epistatic control; thus, such a mutation is beneficial on some genetic backgrounds and deleterious on others. Recent experimental innovations in microbial systems now permit assessment of the fitness effects of individual mutations on multiple genetic backgrounds. We review this literature and identify many examples of sign epistasis, and we suggest that the implications of these results may generalize to other organisms. These theoretical and empirical considerations imply that strong genetic constraint on the selective accessibility of trajectories to high fitness genotypes may exist and suggest specific areas of investigation for future research.     

Epistasis buffers the fitness effects of rifampicin- resistance mutations in Pseudomonas aeruginosa

Epistatic interactions between resistance mutations in antibiotic-free environments potentially play a crucial role in the spread of resistance in pathogen populations by determining the fitness cost associated with resistance. We used an experimental evolution approach to test for epistatic interactions between 14 different pairs of rifampicin mutations in the pathogenic bacterium Pseudomonas aeruginosa in 42 different rifampicin-free environments. First, we show that epistasis between rifampicin-resistance mutations tends to be antagonistic: the fitness effect of having two mutations is generally smaller than that predicted from the effects of individual mutations on the wild-type. Second, we show that sign epistasis between resistance mutations is both common and strong; most notably, pairs of deleterious resistance mutations often partially or completely compensate for each others' costs, revealing a novel mechanism for compensatory adaptation. These results suggest that antagonistic epistasis between intragenic resistance mutations may be a key determinant of the cost of antibiotic resistance and compensatory adaptation in pathogen populations.     

Detecting epistasis from an ensemble of adapting populations

The role that epistasis plays during adaptation remains an outstanding problem, which has received considerable attention in recent years. Most of the recent empirical studies are based on ensembles of replicate populations that adapt in a fixed, laboratory controlled condition. Researchers often seek to infer the presence and form of epistasis in the fitness landscape from the time evolution of various statistics averaged across the ensemble of populations. Here, we provide a rigorous analysis of what quantities, drawn from time series of such ensembles, can be used to infer epistasis for populations evolving under weak mutation on finite‐site fitness landscapes. First, we analyze the mean fitness trajectory—that is, the time course of the ensemble average fitness. We show that for any epistatic fitness landscape and starting genotype, there always exists a non‐epistatic fitness landscape that produces the exact same mean fitness trajectory. Thus, the presence of epistasis is not identifiable from the mean fitness trajectory. By contrast, we show that two other ensemble statistics—the time evolution of the fitness variance across populations, and the time evolution of the mean number of substitutions—can detect certain forms of epistasis in the underlying fitness landscape.     

Variable epistatic effects between mutations at host recognition sites in phiX174 bacteriophage

Epistatic interactions between mutations are widespread. Theoretical investigations have shown that variability in epistatic effects influences fundamental evolutionary processes, yet few empirical studies have identified causes or the extent of this variation. We examined variation in epistatic effects of mutations at two host recognition sites in phiX174 bacteriophage. We calculated epistatic effects from the sum of fitness effects (log scale) of two single mutants and their corresponding double mutant for five combinations of mutations in six conditions. We found that epistatic effects differed in sign, degree, and variability across conditions. The data highlight that even between single mutations at the same two sites the sign and variability of epistatic effects are affected by environment. We discuss these findings in the context of studying the role of epistasis in evolution.     

Impact of epistasis and pleiotropy on evolutionary adaptation

Evolutionary adaptation is often likened to climbing a hill or peak. While this process is simple for fitness landscapes where mutations are independent, the interaction between mutations (epistasis) as well as mutations at loci that affect more than one trait (pleiotropy) are crucial in complex and realistic fitness landscapes. We investigate the impact of epistasis and pleiotropy on adaptive evolution by studying the evolution of a population of asexual haploid organisms (haplotypes) in a model of N interacting loci, where each locus interacts with K other loci. We use a quantitative measure of the magnitude of epistatic interactions between substitutions, and find that it is an increasing function of K. When haplotypes adapt at high mutation rates, more epistatic pairs of substitutions are observed on the line of descent than expected. The highest fitness is attained in landscapes with an intermediate amount of ruggedness that balance the higher fitness potential of interacting genes with their concomitant decreased evolvability. Our findings imply that the synergism between loci that interact epistatically is crucial for evolving genetic modules with high fitness, while too much ruggedness stalls the adaptive process.     

Epistasis between beneficial mutations and the phenotype-to-fitness Map for a ssDNA virus

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic-the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.     

EPISTASIS BETWEEN NEW MUTATIONS AND GENETIC BACKGROUND AND A TEST OF GENETIC CANALIZATION

Abstract The importance for fitness of epistatic interactions among mutations is poorly known, yet epistasis can exert important effects on the dynamics of evolving populations. We showed previously that epistatic interactions are common between pairs of random insertion mutations in the bacterium Escherichia coli . In this paper, we examine interactions between these mutations and other mutations by transducing each of twelve insertion mutations into two genetic backgrounds, one ancestral and the other having evolved in, and adapted to, a defined laboratory environment for 10,000 generations. To assess the effect of the mutation on fitness, we allowed each mutant to compete against its unmutated counterpart in that same environment. Overall, there was a strong positive correlation between the mutational effects on the two genetic backgrounds. Nonetheless, three of the twelve mutations had significantly different effects on the two backgrounds, indicating epistasis. There was no significant tendency for the mutations to be less harmful on the derived background. Thus, there is no evidence supporting the hypothesis that the derived bacteria had adapted, in part, by becoming buffered against the harmful effects of mutations.     

Selection history and epistatic interactions impact dynamics of adaptation to novel environmental stresses

In rapidly changing environments, selection history may impact the dynamics of adaptation. Mutations selected in one environment may result in pleiotropic fitness trade-offs in subsequent novel environments, slowing the rates of adaptation. Epistatic interactions between mutations selected in sequential stressful environments may slow or accelerate subsequent rates of adaptation, depending on the nature of that interaction. We explored the dynamics of adaptation during sequential exposure to herbicides with different modes of action in Chlamydomonas reinhardtii. Evolution of resistance to two of the herbicides was largely independent of selection history. For carbetamide, previous adaptation to other herbicide modes of action positively impacted the likelihood of adaptation to this herbicide. Furthermore, while adaptation to all individual herbicides was associated with pleiotropic fitness costs in stress-free environments, we observed that accumulation of resistance mechanisms was accompanied by a reduction in overall fitness costs. We suggest that antagonistic epistasis may be a driving mechanism that enables populations to more readily adapt in novel environments. These findings highlight the potential for sequences of xenobiotics to facilitate the rapid evolution of multiple-drug and -pesticide resistance, as well as the potential for epistatic interactions between adaptive mutations to facilitate evolutionary rescue in rapidly changing environments.