Mutation-selection balance and the evolutionary advantage of sex and recombination

Mutation-selection balance in a multi-locus system is investigated theoretically, using a modification of Bulmer's infinitesimal model of selection on a normally-distributed quantitative character, taking the number of mutations per individual (n) to represent the character value. The logarithm of the fitness of an individual with n mutations is assumed to be a quadratic, decreasing function of n. The equilibrium properties of infinitely large asexual populations, random-mating populations lacking genetic recombination, and random-mating populations with arbitrary recombination frequencies are investigated. With 'synergistic' epistasis on the scale of log fitness, such that log fitness declines more steeply as n increases, it is shown that equilibrium mean fitness is least for asexual populations. In sexual populations, mean fitness increases with the number of chromosomes and with the map length per chromosome. With 'diminishing returns' epistasis, such that log fitness declines less steeply as n increases, mean fitness behaves in the opposite way. Selection on asexual variants and genes affecting the rate of genetic recombination in random-mating populations was also studied. With synergistic epistasis, zero recombination always appears to be disfavoured, but free recombination is disfavoured when the mutation rate per genome is sufficiently small, leading to evolutionary stability of maps of intermediate length. With synergistic epistasis, an asexual mutant is unlikely to invade a sexual population if the mutation rate per diploid genome greatly exceeds unity. Recombination is selectively disadvantageous when there is diminishing returns epistasis. These results are compared with the results of previous theoretical studies of this problem, and with experimental data.     

Muller''s Ratchet, Epistasis and Mutation Effects

In this study, computer simulation is used to show that despite synergistic epistasis for fitness, Muller's ratchet can lead to lethal fitness loss in a population of asexuals through the accumulation of deleterious mutations. This result contradicts previous work that indicated that epistasis will halt the ratchet. The present results show that epistasis will not halt the ratchet provided that rather than a single deleterious mutation effect, there is a distribution of deleterious mutation effects with sufficient density near zero. In addition to epistasis and mutation distribution, the ability of Muller's ratchet to lead to the extinction of an asexual population under epistasis for fitness depends strongly on the expected number of offspring that survive to reproductive age. This strong dependence is not present in the nonepistatic model and suggests that interpreting the population growth parameter as fecundity is inadequate. Because a continuous distribution of mutation effects is used in this model, an emphasis is placed on the dynamics of the mutation effect distribution rather than on the dynamics of the number of least mutation loaded individuals. This perspective suggests that current models of gene interaction are too simple to apply directly to long-term prediction for populations undergoing the ratchet.     

Epistasis from functional dependence of fitness on underlying traits

Epistasis between mutations in two genes is thought to reflect an interdependence of their functions. While sometimes epistasis is predictable using mechanistic models, its roots seem, in general, hidden in the complex architecture of biological networks. Here, we ask how epistasis can be quantified based on the mathematical dependence of a system-level trait (e.g. fitness) on lower-level traits (e.g. molecular or cellular properties). We first focus on a model in which fitness is the difference between a benefit and a cost trait, both pleiotropically affected by mutations. We show that despite its simplicity, this model can be used to analytically predict certain properties of the ensuing distribution of epistasis, such as a global negative bias, resulting in antagonism between beneficial mutations, and synergism between deleterious ones. We next extend these ideas to derive a general expression for epistasis given an arbitrary functional dependence of fitness on other traits. This expression demonstrates how epistasis relative to fitness can emerge despite the absence of epistasis relative to lower level traits, leading to a formalization of the concept of independence between biological processes. Our results suggest that epistasis may be largely shaped by the pervasiveness of pleiotropic effects and modular organization in biological networks.     

The effects of multilocus balancing selection on neutral variability

We studied the effect of multilocus balancing selection on neutral nucleotide variability at linked sites by simulating a model where diallelic polymorphisms are maintained at an arbitrary number of selected loci by means of symmetric overdominance. Different combinations of alleles define different genetic backgrounds that subdivide the population and strongly affect variability. Several multilocus fitness regimes with different degrees of epistasis and gametic disequilibrium are allowed. Analytical results based on a multilocus extension of the structured coalescent predict that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. Our simulation results show that although variability increases with the number of genetic backgrounds that are maintained in the population, it is reduced by random fluctuations in the frequencies of those backgrounds and does not reach high levels even in very large populations. We also show that previous results on balancing selection in single-locus systems do not extend to the multilocus scenario in a straightforward way. Different patterns of linkage disequilibrium and of the frequency spectrum of neutral mutations are expected under different degrees of epistasis. Interestingly, the power to detect balancing selection using deviations from a neutral distribution of allele frequencies seems to be diminished under the fitness regime that leads to the largest increase of variability over the neutral case. This and other results are discussed in the light of data from the Mhc.     

On the evolution of epistasis III: the haploid case with mutation

Whether interaction between genes is better represented by synergistic or antagonistic epistasis has been a focus of experimental research in bacterial population genetics. Our previous research on evolution of modifiers of epistasis in diploid systems has indicated that the strength of positive or negative epistasis should increase provided linkage disequilibrium is maintained. Here we study a modifier of epistasis in fitness between two loci in a haploid system. Epistasis is modified in the neighborhood of a mutation-selection balance. We show that when linkage in the three-locus system is tight, an increase in the frequency of a modifier allele that induces either more negative or more positive epistasis is possible. Epistasis here can be measured on either an additive or multiplicative scale.     

The decline in fitness with inbreeding: evidence for negative dominance‐by‐dominance epistasis in Drosophila melanogaster

Genetic interactions can play an important role in the evolution of reproductive strategies. In particular, negative dominance‐by‐dominance epistasis for fitness can theoretically favour sex and recombination. This form of epistasis can be detected statistically because it generates nonlinearity in the relationship between fitness and inbreeding coefficient. Measures of fitness in progressively inbred lines tend to show limited evidence for epistasis. However, tests of this kind can be biased against detecting an accelerating decline due to line losses at higher inbreeding levels. We tested for dominance‐by‐dominance epistasis in Drosophila melanogaster by examining viability at five inbreeding levels that were generated simultaneously, avoiding the bias against detecting nonlinearity that has affected previous studies. We find an accelerating rate of fitness decline with inbreeding, indicating that dominance‐by‐dominance epistasis is negative on average, which should favour sex and recombination.     

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.     

Surfing on the seascape: Adaptation in a changing environment

The environment changes constantly at various time scales and, in order to survive, species need to keep adapting. Whether these species succeed in avoiding extinction is a major evolutionary question. Using a multilocus evolutionary model of a mutation‐limited population adapting under strong selection, we investigate the effects of the frequency of environmental fluctuations on adaptation. Our results rely on an “adaptive‐walk” approximation and use mathematical methods from evolutionary computation theory to investigate the interplay between fluctuation frequency, the similarity of environments, and the number of loci contributing to adaptation. First, we assume a linear additive fitness function, but later generalize our results to include several types of epistasis. We show that frequent environmental changes prevent populations from reaching a fitness peak, but they may also prevent the large fitness loss that occurs after a single environmental change. Thus, the population can survive, although not thrive, in a wide range of conditions. Furthermore, we show that in a frequently changing environment, the similarity of threats that a population faces affects the level of adaptation that it is able to achieve. We check and supplement our analytical results with simulations.     

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.     

What lies beneath: natural products from marine organisms as nuclear receptor modulators

The consequences for organism fitness of mutations in a given protein are often thought to be determined to a significant extent by epistasis, that is, by the fact that the effect of a mutation may be strongly dependent on the previous mutational background. Actually, a given mutation could be deleterious or beneficial depending on the background, a situation known as 'sign epistasis'. Under pervasive sign epistasis, many mutational trajectories towards a 'fitter protein' will show a 'dip' in fitness and, it has been previously suggested, only a few trajectories will be available to Darwinian selection. In this issue of the Biochemical Journal, Zhang et al. explore how this simple picture needs to be modified when two rather general and important features are taken into account, namely that many proteins are promiscuous and that living organisms must survive and thrive in environments that change continuously. The multidimensional nature of epistasis for a protein involved in several tasks, together with the fact that different tasks may become critical for organism survival as environmental conditions change, is shown by Zhang et al. to contribute to eliminating fitness dead-ends in protein sequence space. Consequently, many alternative mutational trajectories should allow protein optimization for enhanced organism fitness under changing environmental conditions.     

A simple model of co-evolutionary dynamics caused by epistatic selection

Epistasis is the dependency of the effect of a mutation on the genetic background in which it occurs. Epistasis has been widely documented and implicated in the evolution of species barriers and the evolution of genetic architecture. Here we propose a simple model to formalize the idea that epistasis can also lead to co-evolutionary patterns in molecular evolution of interacting genes. This model epistasis is represented by the influence of one gene substitution on the fitness rank of the resident allele at another locus. We assume that increasing or decreasing fitness rank occur equally likely. In simulations we show that this form of epistasis leads to co-evolution in the sense that the length of an adaptive walk between interacting loci is highly correlated. This effect is caused by episodes of elevated rate of evolution in both loci simultaneously. We find that the influence of epistasis on these measures of co-evolutionary dynamics is relatively robust to the details of the model. The main factor influencing the correlation in evolutionary rates is the probability that a substitution will have an epistatic effect, but the strength of epistasis or the asymmetry of the initial fitness ranks of the alleles have only a minor effect. We suggest that covariance in rates of evolution among loci could be used to detect epistasis among loci.     

The Impact of Macroscopic Epistasis on Long-Term Evolutionary Dynamics

Genetic interactions can strongly influence the fitness effects of individual mutations, yet the impact of these epistatic interactions on evolutionary dynamics remains poorly understood. Here we investigate the evolutionary role of epistasis over 50,000 generations in a well-studied laboratory evolution experiment in Escherichia coli. The extensive duration of this experiment provides a unique window into the effects of epistasis during long-term adaptation to a constant environment. Guided by analytical results in the weak-mutation limit, we develop a computational framework to assess the compatibility of a given epistatic model with the observed patterns of fitness gain and mutation accumulation through time. We find that a decelerating fitness trajectory alone provides little power to distinguish between competing models, including those that lack any direct epistatic interactions between mutations. However, when combined with the mutation trajectory, these observables place strong constraints on the set of possible models of epistasis, ruling out many existing explanations of the data. Instead, we find that the data are consistent with a “two-epoch” model of adaptation, in which an initial burst of diminishing-returns epistasis is followed by a steady accumulation of mutations under a constant distribution of fitness effects. Our results highlight the need for additional DNA sequencing of these populations, as well as for more sophisticated models of epistasis that are compatible with all of the experimental data.     

A sharp minimum on the mean number of steps taken in adaptive walks

It was recently conjectured by H.A. Orr that from a random initial point on a random fitness landscape of alphabetic sequences with one-mutation adjacency, chosen from a larger class of landscapes, no adaptive algorithm can arrive at a local optimum in fewer than on average e-1 steps. Here, using an example in which the mean number of steps to a local optimum equals (A-1)/A, where A is the number of distinct "letters" in the "alphabet" from which sequences are constructed, it is shown that as originally stated, the conjecture does not hold. It is also demonstrated that (A-1)/A is a sharp minimum on the mean number of steps taken in adaptive walks on fitness landscapes of alphabetic sequences with one-mutation adjacency. As the example that achieves the new lower bound has properties that are not often considered as potential attributes for fitness landscapes-non-identically distributed fitnesses and negative fitness correlations for adjacent points-a weaker set of conditions characteristic of more commonly studied fitness landscapes is proposed under which the lower bound on the mean length of adaptive walks is conjectured to equal e-1.     

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.     

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.     

Epistasis and its relationship to canalization in the RNA virus phi 6

Although deleterious mutations are believed to play a critical role in evolution, assessing their realized effect has been difficult. A key parameter governing the effect of deleterious mutations is the nature of epistasis, the interaction between the mutations. RNA viruses should provide one of the best systems for investigating the nature of epistasis because the high mutation rate allows a thorough investigation of mutational effects and interactions. Nonetheless, previous investigations of RNA viruses by S. Crotty and co-workers and by S. F. Elena have been unable to detect a significant effect of epistasis. Here we provide evidence that positive epistasis is characteristic of deleterious mutations in the RNA bacteriophage phi 6. We estimated the effects of deleterious mutations by performing mutation-accumulation experiments on five viral genotypes of decreasing fitness. We inferred positive epistasis because viral genotypes with low fitness were found to be less sensitive to deleterious mutations. We further examined environmental sensitivity in these genotypes and found that low-fitness genotypes were also less sensitive to environmental perturbations. Our results suggest that even random mutations impact the degree of canalization, the buffering of a phenotype against genetic and environmental perturbations. In addition, our results suggest that genetic and environmental canalization have the same developmental basis and finally that an understanding of the nature of epistasis may first require an understanding of the nature of canalization.     

RATES OF FITNESS DECLINE AND REBOUND SUGGEST PERVASIVE EPISTASIS

Unraveling the factors that determine the rate of adaptation is a major question in evolutionary biology. One key parameter is the effect of a new mutation on fitness, which invariably depends on the environment and genetic background. The fate of a mutation also depends on population size, which determines the amount of drift it will experience. Here, we manipulate both population size and genotype composition and follow adaptation of 23 distinct Escherichia coli genotypes. These have previously accumulated mutations under intense genetic drift and encompass a substantial fitness variation. A simple rule is uncovered: the net fitness change is negatively correlated with the fitness of the genotype in which new mutations appear—a signature of epistasis. We find that Fisher's geometrical model can account for the observed patterns of fitness change and infer the parameters of this model that best fit the data, using Approximate Bayesian Computation. We estimate a genomic mutation rate of 0.01 per generation for fitness altering mutations, albeit with a large confidence interval, a mean fitness effect of mutations of ?0.01, and an effective number of traits nine in mutS^? E. coli. This framework can be extended to confront a broader range of models with data and test different classes of fitness landscape models.     

Epistasis,inbreeding depression,and the evolution of self-fertilization

Inbreeding depression resulting from partially recessive deleterious alleles is thought to be the main genetic factor preventing self-fertilizing mutants from spreading in outcrossing hermaphroditic populations. However, deleterious alleles may also generate an advantage to selfers in terms of more efficient purging, while the effects of epistasis among those alleles on inbreeding depression and mating system evolution remain little explored. In this article, we use a general model of selection to disentangle the effects of different forms of epistasis (additive-by-additive, additive-by-dominance, and dominance-by-dominance) on inbreeding depression and on the strength of selection for selfing. Models with fixed epistasis across loci, and models of stabilizing selection acting on quantitative traits (generating distributions of epistasis) are considered as special cases. Besides its effects on inbreeding depression, epistasis may increase the purging advantage associated with selfing (when it is negative on average), while the variance in epistasis favors selfing through the generation of linkage disequilibria that increase mean fitness. Approximations for the strengths of these effects are derived, and compared with individual-based simulation results.     

Gene interactions from maternal effects

Theoretical analyses have demonstrated a potential role for epistasis in many of the most important processes in evolution. These analyses generally assume that an individual's genes map directly to its phenotype and epistasis results from interactions among loci that contribute to the same biochemical or developmental pathways (termed physiological, or within-genotype, epistasis). For many characters, particularly those expressed early in life, an individual's phenotype may also be affected by genes expressed by its parents. The presence of these parental effects allows for interactions between the genes present in the parental and offspring genomes. When the phenotypic effect of a locus in the offspring depends on the alleles possessed by its parents, genotype-by-genotype, or among-genotype, epistasis occurs. The among-genotype epistasis resulting from parental effects may contribute to ruggedness of adaptive landscapes because early mortality often accounts for much of the variance in fitness in populations. To demonstrate how parent-offspring interactions can result in among-genotype epistasis, I use a two-locus model, with one maternal effect locus and one direct effect locus, each with two alleles. Dynamical equations are presented for the two-locus model and are directly contrasted with the dynamical equations derived for a model for physiological epistasis. The relationship between the evolutionary dynamics resulting from these two forms of epistasis is discussed. Three scenarios are presented to illustrate systems in which maternal-offspring, genotype-by-genotype epistasis may occur. The implications of maternal-offspring epistasis for quantitative-trait-loci studies are also discussed.