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.     

A test for epistasis among induced mutations in Caenorhabditis elegans

Synergistic epistasis, in which deleterious mutations tend to magnify each other's effects, is a necessary component of the mutational deterministic hypothesis for the maintenance of sexual production. We tested for epistasis for life-history traits in the soil nematode Caenorhabditis elegans by inducing mutations in two genetic backgrounds: a wild-type strain and a set of genetically loaded lines that contain large numbers of independent mildly detrimental mutations. There was no significant difference between the effect of new mutations on the wild-type background and the genetically loaded background for four out of five fitness correlates. In these four cases, the maximum level of epistasis compatible with the data was very low. The fifth trait, late productivity, is not likely to be an important component of fitness. This suggests either that specific environmental conditions are required to cause epistasis or that synergistic epistasis is not a general phenomenon. We also suggest a new mechanism by which deleterious mutations may provide an advantage to sexual reproduction under low selection coefficients.     

Epistasis in a model of molecular signal transduction

Biological functions typically involve complex interacting molecular networks, with numerous feedback and regulation loops. How the properties of the system are affected when one, or several of its parts are modified is a question of fundamental interest, with numerous implications for the way we study and understand biological processes and treat diseases. This question can be rephrased in terms of relating genotypes to phenotypes: to what extent does the effect of a genetic variation at one locus depend on genetic variation at all other loci? Systematic quantitative measurements of epistasis – the deviation from additivity in the effect of alleles at different loci – on a given quantitative trait remain a major challenge. Here, we take a complementary approach of studying theoretically the effect of varying multiple parameters in a validated model of molecular signal transduction. To connect with the genotype/phenotype mapping we interpret parameters of the model as different loci with discrete choices of these parameters as alleles, which allows us to systematically examine the dependence of the signaling output – a quantitative trait – on the set of possible allelic combinations. We show quite generally that quantitative traits behave approximately additively (weak epistasis) when alleles correspond to small changes of parameters; epistasis appears as a result of large differences between alleles. When epistasis is relatively strong, it is concentrated in a sparse subset of loci and in low order (e.g. pair-wise) interactions. We find that focusing on interaction between loci that exhibit strong additive effects is an efficient way of identifying most of the epistasis. Our model study defines a theoretical framework for interpretation of experimental data and provides statistical predictions for the structure of genetic interaction expected for moderately complex biological circuits.     

Intraspecific trait variation influences physiological performance and fitness in the South Africa shrub genus Protea (Proteaceae)

Background and AimsGlobal plant trait datasets commonly identify trait relationships that are interpreted to reflect fundamental trade-offs associated with plant strategies, but often these trait relationships are not identified when evaluating them at smaller taxonomic and spatial scales. In this study we evaluate trait relationships measured on individual plants for five widespread Protea species in South Africa to determine whether broad-scale patterns of structural trait (e.g. leaf area) and physiological trait (e.g. photosynthetic rates) relationships can be detected within natural populations, and if these traits are themselves related to plant fitness.MethodsWe evaluated the variance structure (i.e. the proportional intraspecific trait variation relative to among-species variation) for nine structural traits and six physiological traits measured in wild populations. We used a multivariate path model to evaluate the relationships between structural traits and physiological traits, and the relationship between these traits and plant size and reproductive effort.Key ResultsWhile intraspecific trait variation is relatively low for structural traits, it accounts for between 50 and 100 % of the variation in physiological traits. Furthermore, we identified few trait associations between any one structural trait and physiological trait, but multivariate regressions revealed clear associations between combinations of structural traits and physiological performance (R² = 0.37–0.64), and almost all traits had detectable associations with plant fitness.ConclusionsIntraspecific variation in structural traits leads to predictable differences in individual-level physiological performance in a multivariate framework, even though the relationship of any particular structural trait to physiological performance may be weak or undetectable. Furthermore, intraspecific variation in both structural and physiological traits leads to differences in plant size and fitness. These results demonstrate the importance of considering measurements of multivariate phenotypes on individual plants when evaluating trait relationships and how trait variation influences predictions of ecological and evolutionary outcomes.     

The capture of heritable variation for genetic quality through social competition

In theory, females of many species choose mates based on traits that are indicators of male genetic quality. A fundamental question in evolutionary biology is why genetic variation for such indicator traits persists despite strong persistent selection imposed by female preference, which is known as the lek paradox. One potential solution to the lek paradox suggests that the traits that are targets of mate choice should evolve condition-dependent expression and that condition should have a large genetic variance. Condition is expected to exhibit high genetic variance because it is affected by a large number of physiological processes and hence, condition-dependent traits should 'capture' variation contributed by a large number of loci. We suggest that a potentially important cause of variation in condition is competition for limited resources. Here, we discuss a pair of models to analyze the evolutionary genetics of traits affected by success in social competition for resources. We show that competition can contribute to genetic variation of 'competition-dependent' traits that have fundamentally different evolutionary properties than other sources of variation. Competition dependence can make traits honest indicators of genetic quality by revealing the relative competitive ability of males, can provide a component of heritable variation that does not contribute to trait evolution, and can help maintain heritable variation under directional selection. Here we provide a general introduction to the concept of competition dependence and briefly introduce two models to demonstrate the potential evolutionary consequences of competition-dependent trait expression.     

Fisher’s Geometrical Model Emerges as a Property of Complex Integrated Phenotypic Networks

Models relating phenotype space to fitness (phenotype–fitness landscapes) have seen important developments recently. They can roughly be divided into mechanistic models (e.g., metabolic networks) and more heuristic models like Fisher’s geometrical model. Each has its own drawbacks, but both yield testable predictions on how the context (genomic background or environment) affects the distribution of mutation effects on fitness and thus adaptation. Both have received some empirical validation. This article aims at bridging the gap between these approaches. A derivation of the Fisher model “from first principles” is proposed, where the basic assumptions emerge from a more general model, inspired by mechanistic networks. I start from a general phenotypic network relating unspecified phenotypic traits and fitness. A limited set of qualitative assumptions is then imposed, mostly corresponding to known features of phenotypic networks: a large set of traits is pleiotropically affected by mutations and determines a much smaller set of traits under optimizing selection. Otherwise, the model remains fairly general regarding the phenotypic processes involved or the distribution of mutation effects affecting the network. A statistical treatment and a local approximation close to a fitness optimum yield a landscape that is effectively the isotropic Fisher model or its extension with a single dominant phenotypic direction. The fit of the resulting alternative distributions is illustrated in an empirical data set. These results bear implications on the validity of Fisher’s model’s assumptions and on which features of mutation fitness effects may vary (or not) across genomic or environmental contexts.     

A Developmental Systems Perspective on Epistasis: Computational Exploration of Mutational Interactions in Model Developmental Regulatory Networks

The way in which the information contained in genotypes is translated into complex phenotypic traits (i.e. embryonic expression patterns) depends on its decoding by a multilayered hierarchy of biomolecular systems (regulatory networks). Each layer of this hierarchy displays its own regulatory schemes (i.e. operational rules such as +/− feedback) and associated control parameters, resulting in characteristic variational constraints. This process can be conceptualized as a mapping issue, and in the context of highly-dimensional genotype-phenotype mappings (GPMs) epistatic events have been shown to be ubiquitous, manifested in non-linear correspondences between changes in the genotype and their phenotypic effects. In this study I concentrate on epistatic phenomena pervading levels of biological organization above the genetic material, more specifically the realm of molecular networks. At this level, systems approaches to studying GPMs are specially suitable to shed light on the mechanistic basis of epistatic phenomena. To this aim, I constructed and analyzed ensembles of highly-modular (fully interconnected) networks with distinctive topologies, each displaying dynamic behaviors that were categorized as either arbitrary or functional according to early patterning processes in the Drosophila embryo. Spatio-temporal expression trajectories in virtual syncytial embryos were simulated via reaction-diffusion models. My in silico mutational experiments show that: 1) the average fitness decay tendency to successively accumulated mutations in ensembles of functional networks indicates the prevalence of positive epistasis, whereas in ensembles of arbitrary networks negative epistasis is the dominant tendency; and 2) the evaluation of epistatic coefficients of diverse interaction orders indicates that, both positive and negative epistasis are more prevalent in functional networks than in arbitrary ones. Overall, I conclude that the phenotypic and fitness effects of multiple perturbations are strongly conditioned by both the regulatory architecture (i.e. pattern of coupled feedback structures) and the dynamic nature of the spatio-temporal expression trajectories displayed by the simulated networks.     

A population genetics model for multiple quantitative traits exhibiting pleiotropy and epistasis

We study a population genetics model of an organism with a genome of L(tot)loci that determine the values of T quantitative traits. Each trait is controlled by a subset of L loci assigned randomly from the genome. There is an optimum value for each trait, and stabilizing selection acts on the phenotype as a whole to maintain actual trait values close to their optima. The model contains pleiotropic effects (loci can affect more than one trait) and epistasis in fitness. We use adaptive walk simulations to find high-fitness genotypes and to study the way these genotypes are distributed in sequence space. We then simulate the evolution of haploid and diploid populations on these fitness landscapes and show that the genotypes of populations are able to drift through sequence space despite stabilizing selection on the phenotype. We study the way the rate of drift and the extent of the accessible region of sequence space is affected by mutation rate, selection strength, population size, recombination rate, and the parameters L and T that control the landscape shape. There are three regimes of the model. If LTL(tot), there are many small peaks that can be spread over a wide region of sequence space. Compensatory neutral mutations are important in the population dynamics in this case.     

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.     

The effect of induced mutations on quantitative traits in Arabidopsis thaliana: Natural versus artificial conditions

Mutations are the ultimate source of all genetic variations. New mutations are expected to affect quantitative traits differently depending on the extent to which traits contribute to fitness and the environment in which they are tested. The dogma is that the preponderance of mutations affecting fitness will be skewed toward deleterious while their effects on nonfitness traits will be bidirectionally distributed. There are mixed views on the role of stress in modulating these effects. We quantify mutation effects by inducing mutations in Arabidopsis thaliana (Columbia accession) using the chemical ethylmethane sulfonate. We measured the effects of new mutations relative to a premutation founder for fitness components under both natural (field) and artificial (growth room) conditions. Additionally, we measured three other quantitative traits, not expected to contribute directly to fitness, under artificial conditions. We found that induced mutations were equally as likely to increase as decrease a trait when that trait was not closely related to fitness (traits that were neither survivorship nor reproduction). We also found that new mutations were more likely to decrease fitness or fitness‐related traits under more stressful field conditions than under relatively benign artificial conditions. In the benign condition, the effect of new mutations on fitness components was similar to traits not as closely related to fitness. These results highlight the importance of measuring the effects of new mutations on fitness and other traits under a range of conditions.     

Condition dependence of sexually dimorphic colouration and longevity in the ambush bug Phymata americana

Sexually selected traits that are costly are predicted to be more condition dependent than nonsexually selected traits. Assuming resource limitation, increased allocation to a sexually selected trait may also come at a cost to other fitness components. To test these predictions, we varied adult food ration to manipulate condition in the colour dimorphic bug, Phymata americana. We compared the degree of condition dependence in a sexually selected trait expressed in males to a nonsexually selected trait expressed in males and females. We also evaluated the effects of condition on longevity of both sexes. We found that the expression of these colour pattern traits was strongly influenced by both diet and age. As expected, the strength of condition dependence was much more pronounced in the sexually selected, male-limited trait but the nonsexual trait also exhibited significant condition dependence in both sexes. The sexually selected male trait also exhibited a higher coefficient of phenotypic variation than the nonsexually selected trait in males and females. Diet had contrasting effects on male and female longevity; increased food availability had positive effects on female lifespan but these effects were not detected in males, suggesting that males allocated limited resources preferentially to sexually selected traits. These results are consistent with the expectation that optimal allocation to various fitness components differs between the sexes.     

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.     

Patterns and Causes of Signed Linkage Disequilibria in Flies and Plants

Most empirical studies of linkage disequilibrium (LD) study its magnitude, ignoring its sign. Here, we examine patterns of signed LD in two population genomic data sets, one from Capsella grandiflora and one from Drosophila melanogaster. We consider how processes such as drift, admixture, Hill–Robertson interference, and epistasis may contribute to these patterns. We report that most types of mutations exhibit positive LD, particularly, if they are predicted to be less deleterious. We show with simulations that this pattern arises easily in a model of admixture or distance-biased mating, and that genome-wide differences across site types are generally expected due to differences in the strength of purifying selection even in the absence of epistasis. We further explore how signed LD decays on a finer scale, showing that loss of function mutations exhibit particularly positive LD across short distances, a pattern consistent with intragenic antagonistic epistasis. Controlling for genomic distance, signed LD in C. grandiflora decays faster within genes, compared with between genes, likely a by-product of frequent recombination in gene promoters known to occur in plant genomes. Finally, we use information from published biological networks to explore whether there is evidence for negative synergistic epistasis between interacting radical missense mutations. In D. melanogaster networks, we find a modest but significant enrichment of negative LD, consistent with the possibility of intranetwork negative synergistic epistasis.     

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.     

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.     

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.     

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.     

Genetic complexity and quantitative trait loci mapping of yeast morphological traits

Functional genomics relies on two essential parameters: the sensitivity of phenotypic measures and the power to detect genomic perturbations that cause phenotypic variations. In model organisms, two types of perturbations are widely used. Artificial mutations can be introduced in virtually any gene and allow the systematic analysis of gene function via mutants fitness. Alternatively, natural genetic variations can be associated to particular phenotypes via genetic mapping. However, the access to genome manipulation and breeding provided by model organisms is sometimes counterbalanced by phenotyping limitations. Here we investigated the natural genetic diversity of Saccharomyces cerevisiae cellular morphology using a very sensitive high-throughput imaging platform. We quantified 501 morphological parameters in over 50,000 yeast cells from a cross between two wild-type divergent backgrounds. Extensive morphological differences were found between these backgrounds. The genetic architecture of the traits was complex, with evidence of both epistasis and transgressive segregation. We mapped quantitative trait loci (QTL) for 67 traits and discovered 364 correlations between traits segregation and inheritance of gene expression levels. We validated one QTL by the replacement of a single base in the genome. This study illustrates the natural diversity and complexity of cellular traits among natural yeast strains and provides an ideal framework for a genetical genomics dissection of multiple traits. Our results did not overlap with results previously obtained from systematic deletion strains, showing that both approaches are necessary for the functional exploration of genomes.     

An experimental assessment of artificial within-family selection for fitness in conservation programs

Equalizing familiar contributions is the simplest recommended strategy to maintain genetic diversity in conservation programs. However, this method implies a relaxation of natural selection and the possibility of accumulation of deleterious mutations. Computer simulations have shown that performing selection within families for fitness traits in a conservation program can be useful to alleviate such problems. We thus carried out an experiment with the model species Drosophila melanogaster in order to assess whether or not selection for fitness traits can be useful. We considered a fitness trait (pupa productivity) that was first checked to perform as a typical fitness component. The trait showed an inbreeding depression of 1.2 per 1 % increase in inbreeding and an asymmetrical response to selection with average realized heritabilities of about 0.04 in the upward direction and an order of magnitude larger (0.36) in the downward direction. The management experiment indicated that artificial within-family selection for fitness had only a marginal success for two reasons. First, there was not an appreciable decline in fitness across the experiment despite the low population sizes assumed (N = 10 or 20), even in the population not subjected to selection. This result is compatible with fitness models which imply the segregation of few deleterious mutations of large effect. Second, artificial selection within families had a limited impact on the trait, as one expects for a typical fitness component with very low heritability.     

Developmental associations between traits: covariance and beyond

Statistical associations between phenotypic traits often result from shared developmental processes and include both covariation between the trait values and more complex associations between higher moments of the joint distribution of traits. In this article, an analytical technique for calculating the covariance between traits is presented on the basis of (1). the distribution of underlying genetic and environmental variation that jointly influences the traits and (2). the mechanics of how these underlying factors influence the development of each trait. It is shown that epistasis can produce patterns of covariation between traits that are not seen in additive models. Applying this approach to a trait in parents and the same trait in their offspring allows us to study the consequences of epistasis for the evolution of additive genetic variance and heritability. This analysis is then extended to the study of more complicated associations between traits. It is shown that even traits that are not correlated may exhibit developmental associations that influence their joint evolution.