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.     

Evolutionary design on a budget: robustness and optimality of bacteriophage T7

Exploring how biological systems have been 'designed' by evolution to achieve robust behaviours is now a subject of increasing research effort. Yet, it still remains unclear how environmental factors may contribute to this process. This issue is addressed by employing a detailed computer model for the intracellular growth of phage T7. More than 150 000 in silico T7 mutants were generated and the rates and efficiencies of their growth in two host environments, namely, a realistic environment that offered finite host resources for the synthesis of phage functions and a hypothetical environment where the phage was supplied infinite host resources, were evaluated. Results revealed two key properties of phage T7. First, T7 growth was overall robust with respect to perturbations in its parameters, but fragile with respect to changes in the ordering of its genetic elements. Secondly, the wild-type T7 had close to optimal fitness in the finite environment. Furthermore, a strong correlation was found between fitness and growth efficiency in the finite environment. The results underscore the potential importance of the environment in shaping robust design of a biological system. In particular, the strong correlation between fitness and growth efficiency suggests that T7 may have evolved to maximise its growth rate by minimising waste of finite resources.     

Evolution can favor antagonistic epistasis

The accumulation of deleterious mutations plays a major role in evolution, and key to this are the interactions between their fitness effects, known as epistasis. Whether mutations tend to interact synergistically (with multiple mutations being more deleterious than would be expected from their individual fitness effects) or antagonistically is important for a variety of evolutionary questions, particularly the evolution of sex. Unfortunately, the experimental evidence on the prevalence and strength of epistasis is mixed and inconclusive. Here we study theoretically whether synergistic or antagonistic epistasis is likely to be favored by evolution and by how much. We find that in the presence of recombination, evolution favors less synergistic or more antagonistic epistasis whenever mutations that change the epistasis in this direction are possible. This is because evolution favors increased buffering against the effects of deleterious mutations. This suggests that we should not expect synergistic epistasis to be widespread in nature and hence that the mutational deterministic hypothesis for the advantage of sex may not apply widely.     

An Experimental Test for Synergistic Epistasis and Its Application in Chlamydomonas

Theoretically, one of the most general benefits of sex is given by its function in facilitating selection against deleterious mutations. This advantage of sex may be deterministic if deleterious mutations affect the fitness of an individual in a synergistic way, i.e., if mutations increase each others' negative fitness effect. We present a new test for synergistic epistasis that considers the skewness of the log fitness distribution of offspring from a cross. We applied this test to data of the unicellular alga Chlamydomonas moewussii. For this purpose, two crosses were made: one between two strains that are presumed to have accumulated slightly deleterious mutations, the other between two strains without a history of mutation accumulation. Fitness was measured by estimating the two parameters of logistic growth in batch culture, the maximum growth rate (r) and the carrying capacity (K). The finding of a negatively skewed distribution for K in the accumulation cross suggests synergism between mutations affecting the carrying capacity, while the absence of skewness for r in both crosses is consistent with independent effects of mutations affecting this parameter. The results suggest a possible alternative explanation for the general observation that sex is related to constant environments, where selection on K predominates, while asexual reproduction is found in more variable environments, where selection on r is more important.     

The Environment Affects Epistatic Interactions to Alter the Topology of an Empirical Fitness Landscape

The fitness effect of mutations can be influenced by their interactions with the environment, other mutations, or both. Previously, we constructed 32 ( = 2⁵) genotypes that comprise all possible combinations of the first five beneficial mutations to fix in a laboratory-evolved population of Escherichia coli. We found that (i) all five mutations were beneficial for the background on which they occurred; (ii) interactions between mutations drove a diminishing returns type epistasis, whereby epistasis became increasingly antagonistic as the expected fitness of a genotype increased; and (iii) the adaptive landscape revealed by the mutation combinations was smooth, having a single global fitness peak. Here we examine how the environment influences epistasis by determining the interactions between the same mutations in two alternative environments, selected from among 1,920 screened environments, that produced the largest increase or decrease in fitness of the most derived genotype. Some general features of the interactions were consistent: mutations tended to remain beneficial and the overall pattern of epistasis was of diminishing returns. Other features depended on the environment; in particular, several mutations were deleterious when added to specific genotypes, indicating the presence of antagonistic interactions that were absent in the original selection environment. Antagonism was not caused by consistent pleiotropic effects of individual mutations but rather by changing interactions between mutations. Our results demonstrate that understanding adaptation in changing environments will require consideration of the combined effect of epistasis and pleiotropy across environments.     

Population Growth Inflates the Per-Individual Number of Deleterious Mutations and Reduces Their Mean Effect

This study addresses the question of how purifying selection operates during recent rapid population growth such as has been experienced by human populations. This is not a straightforward problem because the human population is not at equilibrium: population genetics predicts that, on the one hand, the efficacy of natural selection increases as population size increases, eliminating ever more weakly deleterious variants; on the other hand, a larger number of deleterious mutations will be introduced into the population and will be more likely to increase in their number of copies as the population grows. To understand how patterns of human genetic variation have been shaped by the interaction of natural selection and population growth, we examined the trajectories of mutations with varying selection coefficients, using computer simulations. We observed that while population growth dramatically increases the number of deleterious segregating sites in the population, it only mildly increases the number carried by each individual. Our simulations also show an increased efficacy of natural selection, reflected in a higher fraction of deleterious mutations eliminated at each generation and a more efficient elimination of the most deleterious ones. As a consequence, while each individual carries a larger number of deleterious alleles than expected in the absence of growth, the average selection coefficient of each segregating allele is less deleterious. Combined, our results suggest that the genetic risk of complex diseases in growing populations might be distributed across a larger number of more weakly deleterious rare variants.     

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs’ deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.     

Parasites and deleterious mutations: interactions influencing the evolutionary maintenance of sex

The restrictive assumptions associated with purely genetic and purely ecological mechanisms suggest that neither of the two forces, in isolation, can offer a general explanation for the evolutionary maintenance of sex. Consequently, attention has turned to pluralistic models (i.e. models that apply both ecological and genetic mechanisms). Existing research has shown that combining mutation accumulation and parasitism allows restrictive assumptions about genetic and parasite parameter values to be relaxed while still predicting the maintenance of sex. However, several empirical studies have shown that deleterious mutations and parasitism can reduce fitness to a greater extent than would be expected if the two acted independently. We show how interactions between these genetic and ecological forces can completely reverse predictions about the evolution of reproductive modes. Moreover, we demonstrate that synergistic interactions between infection and deleterious mutations can render sex evolutionarily stable even when there is antagonistic epistasis among deleterious mutations, thereby widening the conditions for the evolutionary maintenance of sex.     

TEST OF INTERACTION BETWEEN GENETIC MARKERS THAT AFFECT FITNESS IN ASPERGILLUS NIGER

In this paper we study whether and how a number of arbitrarily chosen marker mutations interact in their effect on fitness, which is relevant for our understanding of the evolution of sex. If epistasis is synergistic, the main function of sex may be to facilitate selection against deleterious mutations. We use strains of the filamentous fungus Aspergillus niger with variable combinations of marker mutations that have been obtained by isolating segregants from a diploid between a wild-type strain and a related strain carrying a marker mutation on each of its eight chromosomes. The marker mutations include five auxotrophic and two resistance mutations. As a measure of fitness the mycelium growth rate on supplemented medium has been used. The results suggest that the marker mutations have independent effects on fitness, and hence they do not support the deterministic mutation hypothesis of the evolution of sex. The apparent linear relationship between mutation number and log fitness is the result of interactions of opposite type (i.e., synergistic and antagonistic) that cancel each other's effect. However, due to an isolation bias caused by the fact that not all possible strains with many mutations could be isolated, the results may be relatively biased towards an antagonistic relationship between mutation number and log fitness.     

The population genetics of X-autosome synthetic lethals and steriles

Epistatic interactions are widespread, and many of these interactions involve combinations of alleles at different loci that are deleterious when present in the same individual. The average genetic environment of sex-linked genes differs from that of autosomal genes, suggesting that the population genetics of interacting X-linked and autosomal alleles may be complex. Using both analytical theory and computer simulations, we analyzed the evolutionary trajectories and mutation-selection balance conditions for X-autosome synthetic lethals and steriles. Allele frequencies follow a set of fundamental trajectories, and incompatible alleles are able to segregate at much higher frequencies than single-locus expectations. Equilibria exist, and they can involve fixation of either autosomal or X-linked alleles. The exact equilibrium depends on whether synthetic alleles are dominant or recessive and whether fitness effects are seen in males, females, or both sexes. When single-locus fitness effects and synthetic incompatibilities are both present, population dynamics depend on the dominance of alleles and historical contingency (i.e., whether X-linked or autosomal mutations occur first). Recessive synthetic lethality can result in high-frequency X-linked alleles, and dominant synthetic lethality can result in high-frequency autosomal alleles. Many X-autosome incompatibilities in natural populations may be cryptic, appearing to be single-locus effects because one locus is fixed. We also discuss the implications of these findings with respect to standing genetic variation and the origins of Haldane's rule.     

THE CONTRIBUTION OF NEW MUTATIONS TO GENOTYPE-ENVIRONMENT INTERACTION FOR FITNESS IN DROSOPHILA MELANOGASTER

Many studies have documented the existence of genotype-environment interaction (GEI) for traits closely related to fitness in natural populations. A type of GEI that is commonly observed is changes in the fitness ranking of genetic groups (families, clones, or inbred lines) in different environments. We refer to such changes in ranking as crossing of reaction norms for fitness. A common interpretation of crossing of reaction norms for fitness is that selection favors different alleles in the different environments (i.e., that “trade-offs” exist). If this is the case, selection could maintain genetic variation, and even lead to reproductive isolation between subpopulations using different environments. Even if the same alleles are favored in every environment, however, deleterious mutations that vary in the magnitude of their effect depending on environment could cause reaction norms for fitness to cross. If deleterious mutations with environment-dependent effects are responsible for maintaining much of the variation leading to crossing of reaction norms for fitness in natural populations, it should be possible to observe crossing of reaction norms for fitness among otherwise genetically identical lines bearing newly arisen spontaneous mutations. We examined the contribution of new mutations to GEI for fitness in Drosophila melanogaster. Eighteen lines were derived from a common, highly inbred base stock, and maintained at a population size of 10 pairs for over 200 generations, to allow them to accumulate spontaneous mutations. Because of the small population size of the lines, selection against mildly deleterious mutations should have been relatively ineffective. The lines were tested for productivity (number of surviving adult progeny from a standard number of parents) in five different environmental treatments, comprising different food media, temperatures, and levels of competition. The lines showed highly significant GEI for productivity, owing largely to considerable changes in ranking in the different environments. We conclude that mutations that are deleterious on average, but whose quantitative effects depend on environment, could be responsible for maintaining much of the variation leading to crossing of reaction norms for fitness that has been observed in samples of D. melanogaster from the wild.     

Factors affecting the genetic load in Drosophila: synergistic epistasis and correlations among fitness components

Two factors that can affect genetic load, synergistic epistasis and sexual selection, were investigated in Drosophila melanogaster. A set of five chromosomal regions containing visible recessive mutations were put together in all combinations to create a full set of 32 homozygous lines fixed for different numbers of known mutations. Two measures of fitness were made for each line: productivity (a combined measure of fecundity and egg-to-adult survivorship) and competitive male mating success. Productivity, but not male mating success, showed a pattern of strong average synergistic epistasis, such that the log fitness declined nonlinearly with increasing numbers of mutations. Synergistic epistasis is known to reduce the mutation load. Both fitness components show some positive and some negative interactions between specific sets of mutations. Furthermore, alleles with deleterious effects on productivity tend to also diminish male mating success. Given that male mating success can affect relative fitness without changing the mean productivity of a population, these additional effects would lead to lower frequencies and lower fixation rates of deleterious alleles without higher costs to the mean fitness of the population.     

Pleiotropic Effects of a Mitochondrial–Nuclear Incompatibility Depend upon the Accelerating Effect of Temperature in Drosophila

Interactions between mitochondrial and nuclear gene products that underlie eukaryotic energy metabolism can cause the fitness effects of mutations in one genome to be conditional on variation in the other genome. In ectotherms, the effects of these interactions are likely to depend upon the thermal environment, because increasing temperature accelerates molecular rates. We find that temperature strongly modifies the pleiotropic phenotypic effects of an incompatible interaction between a Drosophila melanogaster polymorphism in the nuclear-encoded, mitochondrial tyrosyl-transfer (t)RNA synthetase and a D. simulans polymorphism in the mitochondrially encoded tRNA^Tyr. The incompatible mitochondrial–nuclear genotype extends development time, decreases larval survivorship, and reduces pupation height, indicative of decreased energetic performance. These deleterious effects are ameliorated when larvae develop at 16° and exacerbated at warmer temperatures, leading to complete sterility in both sexes at 28°. The incompatible genotype has a normal metabolic rate at 16° but a significantly elevated rate at 25°, consistent with the hypothesis that inefficient energy metabolism extends development in this genotype at warmer temperatures. Furthermore, the incompatibility decreases metabolic plasticity of larvae developed at 16°, indicating that cooler development temperatures do not completely mitigate the deleterious effects of this genetic interaction. Our results suggest that the epistatic fitness effects of metabolic mutations may generally be conditional on the thermal environment. The expression of epistatic interactions in some environments, but not others, weakens the efficacy of selection in removing deleterious epistatic variants from populations and may promote the accumulation of incompatibilities whose fitness effects will depend upon the environment in which hybrids occur.     

Abiotic stress does not magnify the deleterious effects of spontaneous mutations

Although the effects of deleterious alleles often are predicted to be greater in stressful environments, there is no theoretical basis for this prediction and the empirical evidence is mixed. Here we characterized the effects of three types of abiotic stress (thermal, oxidative and hyperosmotic) on two sets of nematode (Caenorhabditis elegans) mutation accumulation (MA) lines that differ by threefold in fitness. We compared the survival and egg-to-adult viability between environments (benign and stressful) and between fitness categories (high-fitness MA, low-fitness MA). If the environment and mutation load have synergistic effects on trait means, then the difference between the high and low-fitness MA lines should be larger in stressful environments. Although the stress treatments consistently decreased survival and/or viability, we did not detect significant interactions between fitness categories and environment types. In contrast, we did find consistent evidence for synergistic effects on (micro)environmental variation. The lack of signal in trait means likely reflects the very low starting fitness of some low-fitness MA lines, the potential for cross-stress responses and the context dependence of mutational effects. In addition, the large increases in the environmental variance in the stressful environments may have masked small changes in trait means. These results do not provide evidence for synergism between mutation and stress.     

Estimation of Deleterious-Mutation Parameters in Natural Populations

The rate and average effects of spontaneous deleterious mutations are important determinants of the evolution of breeding systems and of the vulnerability of small populations to extinction. Nevertheless, few attempts have been made to estimate the properties of such mutations, and those studies that have been performed have been extremely labor intensive, relying on long-term, laboratory mutation-accumulation experiments. We present an alternative to the latter approach. For populations in which the genetic variance for fitness is a consequence of selection-mutation balance, the mean fitness and genetic variance of fitness in outbred and inbred generations can be expressed as simple functions of the genomic mutation rate, average homozygous effect and average dominance coefficient of new mutations. Using empirical estimates for the mean and genetic variance of fitness, these expressions can then be solved to obtain joint estimates of the deleterious-mutation parameters. We employ computer simulations to evaluate the degree of bias of the estimators and present some general recommendations on the application of the technique. Our procedures provide some hope for obtaining estimates of the properties of deleterious mutations from a wide phylogenetic range of species as well as a mechanism for testing the validity of alternative models for the maintenance of genetic variance for fitness.     

The fitness effect of mutations across environments: Fisher's geometrical model with multiple optima

When are mutations beneficial in one environment and deleterious in another? More generally, what is the relationship between mutation effects across environments? These questions are crucial to predict adaptation in heterogeneous conditions in a broad sense. Empirical evidence documents various patterns of fitness effects across environments but we still lack a framework to analyze these multivariate data. In this article, we extend Fisher's geometrical model to multiple environments determining distinct peaks. We derive the fitness distribution, in one environment, among mutants with a given fitness in another and the bivariate distribution of random mutants’ fitnesses across two or more environments. The geometry of the phenotype‐fitness landscape is naturally interpreted in terms of fitness trade‐offs between environments. These results may be used to fit/predict empirical distributions or to predict the pattern of adaptation across heterogeneous conditions. As an example, we derive the genomic rate of substitution and of adaptation in a metapopulation divided into two distinct habitats in a high migration regime and show that they depend critically on the geometry of the phenotype‐fitness landscape.     

The effects of parasitism and inbreeding on the competitive ability in Daphnia magna: evidence for synergistic epistasis

Synergistic epistasis for fitness is often assumed in models of how selection acts on the frequency and distribution of deleterious mutations. Evidence for synergistic epistasis would exist if the logarithm of fitness declines more quickly with number of deleterious mutations, than predicted by a linear decline. This can be studied indirectly by quantifying the effect of different levels of inbreeding on fitness. Here, six sets (different genetic backgrounds) of three increasingly inbred Daphnia magna clones were used to assess their relative fitness according to changes in frequency in a competition experiment against a tester clone. A novelty of the mating procedure was that the inbreeding coefficients (F) of the three clones belonging to each set increased in steps of 0.25 independent of the (unknown) inbreeding coefficient of the common ancestor. The equal increase of the inbreeding coefficients is important, because deviations influence the quantification of inbreeding depression, its variance and the detection of epistasis. In a simple mathematical model we show that when working with a partially inbred population inbreeding depression is underestimated, the variance of fitness is increased, and the detection of epistasis more difficult. Further, to examine whether an interaction between inbreeding and parasitism exists, each inbred clone was tested with and without a microsporidium infection (Octosporea bayeri). We found a nonlinear decrease of the logarithm of fitness across the three levels of inbreeding, indicating synergistic epistasis. The interaction term between parasitism and inbreeding was not significant. Our results suggest that deleterious mutations may be purged effectively once the level of inbreeding is high, but that parasitism seems not to influence this effect.     

Environment changes epistasis to alter trade‐offs along alternative evolutionary paths

The fitness effect of a mutation can depend on both its genetic background, known as epistasis, and the prevailing external environment. Many examples of these dependencies are known, but few studies consider both aspects in combination, especially as they affect mutations that have been selected together. We examine interactions between five coevolved mutations in eight diverse environments. We find that mutations are, on average, beneficial across environments, but that there is high variation in their fitness effects, including many examples of mutations conferring a cost in some, but not other, genetic background‐environment combinations. Indeed, even when global interaction trends are accounted for, specific local mutation interactions are common and differed across environments. One consequence of this dependence is that the range of trade‐offs in genotype fitness across selected and alternative environments are contingent on the particular evolutionary path followed over the mutation landscape. Finally, although specific interactions were common, there was a consistent pattern of diminishing returns epistasis whereby mutation effects were less beneficial when added to genotypes of higher fitness. Our results underline that specific mutation effects are highly dependent on the combination of genetic and external environments, and support a general relationship between a genotype's current fitness and its potential to increase in fitness.     

Emergence of long‐term balanced polymorphism under cyclic selection of spatially variable magnitude

A fundamental question in evolutionary biology is what promotes genetic variation at nonneutral loci, a major precursor to adaptation in changing environments. In particular, balanced polymorphism under realistic evolutionary models of temporally varying environments in finite natural populations remains to be demonstrated. Here, we propose a novel mechanism of balancing selection under temporally varying fitnesses. Using forward‐in‐time computer simulations and mathematical analysis, we show that cyclic selection that spatially varies in magnitude, such as along an environmental gradient, can lead to elevated levels of nonneutral genetic polymorphism in finite populations. Balanced polymorphism is more likely with an increase in gene flow, magnitude and period of fitness oscillations, and spatial heterogeneity. This polymorphism‐promoting effect is robust to small systematic fitness differences between competing alleles or to random environmental perturbation. Furthermore, we demonstrate analytically that protected polymorphism arises as spatially heterogeneous cyclic fitness oscillations generate a type of storage effect that leads to negative frequency dependent selection. Our findings imply that spatially variable cyclic environments can promote elevated levels of nonneutral genetic variation in natural populations.     

Epistasis among Deleterious Mutations in the HIV-1 Protease

A central goal in molecular evolution is to understand how genetic interactions between protein mutations shape protein function and fitness. While intergenic epistasis has been extensively explored in eukaryotes, bacteria, and viruses, intragenic epistatic interactions have been insufficiently studied. Here, we employ a model system in which lambda phage fitness correlates with the enzymatic activity of human immunodeficiency virus type 1 (HIV-1) protease to systematically determine the epistatic interactions between intragenic pairs of deleterious protein substitutions. We generated 114 genotypes of the HIV-1 protease, each carrying pairs of nucleotide substitution mutations whose separated and combined deleterious effects on fitness were then determined. A high proportion (39%) of pairs displayed lethality. Several pairs exhibited significant interactions for fitness, including positive and negative epistasis. Significant negative epistatic interactions predominated (15%) over positive interactions (2%). However, the average ± SD epistatic effect, ē = 0.0025 ± 0.1334, was not significantly different from zero (p = 0.8368). Notably, epistatic interactions, regardless of epistatic direction, tend to be more frequent in the context of less deleterious mutations. In the present study, the high frequencies of lethality and negative epistasis indicate that the HIV-1 protease is highly sensitive to the effects of deleterious mutations. Therefore, proteins may not be as robust to mutational change as is usually expected.