Evolution of coadaptation in a two-locus epistatic system

Although recent advances in genome biology have dramatically increased our understanding of the contribution of gene interactions to the development of complex phenotypes, we still lack general agreement on the process and mechanisms responsible for the evolution of epistatic systems. Even if genes in a species are indeed integrated into coadapted complexes of interacting components, simple additive evolution may eventually result in epistatic differentiation of populations. Consequently, the prevalence of epistatic gene action does not tell us anything about the role of epistatic selection in the history of population divergence. To elucidate the contribution of epistatic selection in the evolution of coadaptation, we investigate the fixation process of two mutations that interact synergistically to enhance fitness. We show by diffusion analysis and simulations that epistatic selection on cosegregating variants does not by itself promote the evolution of epistatic systems; rather, accumulation of neutral mutations may play a crucial role, creating an appropriate genetic milieu for adaptive evolution in the future generations.     

Haldane's sieve and adaptation from the standing genetic variation

We consider populations that adapt to a sudden environmental change by fixing alleles found at mutation-selection balance. In particular, we calculate probabilities of fixation for previously deleterious alleles, ignoring the input of new mutations. We find that "Haldane's sieve"--the bias against the establishment of recessive beneficial mutations--does not hold under these conditions. Instead probabilities of fixation are generally independent of dominance. We show that this result is robust to patterns of sex expression for both X-linked and autosomal loci. We further show that adaptive evolution is invariably slower at X-linked than autosomal loci when evolution begins from mutation-selection balance. This result differs from that obtained when adaptation uses new mutations, a finding that may have some bearing on recent attempts to distinguish between hitchhiking and background selection by contrasting the molecular population genetics of X-linked vs. autosomal loci. Last, we suggest a test to determine whether adaptation used new mutations or previously deleterious alleles from the standing genetic variation.     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

SELECTION BIASES THE PREVALENCE AND TYPE OF EPISTASIS ALONG ADAPTIVE TRAJECTORIES

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

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.     

Adaptation from standing genetic variation

Populations adapt to novel environments in two distinct ways: selection on pre-existing genetic variation and selection on new mutations. These alternative sources of beneficial alleles can result in different evolutionary dynamics and distinct genetic outcomes. Compared with new mutations, adaptation from standing genetic variation is likely to lead to faster evolution, the fixation of more alleles of small effect and the spread of more recessive alleles. There is potential to distinguish between adaptation from standing variation and that from new mutations by differences in the genomic signature of selection. Here we review these approaches and possible examples of adaptation from standing variation in natural populations. Understanding how the source of genetic variation affects adaptation will be integral for predicting how populations will respond to changing environments.     

Origin, fate, and architecture of ecologically relevant genetic variation

Recent advances in molecular genetics combined with field manipulations are yielding new insight into the origin, evolutionary fate, and genetic architecture of phenotypic variation in natural plant populations, with two surprising implications for the evolution of plant genomes. First, genetic loci exhibiting antagonistic pleiotropy across natural environments appear rare relative to loci that are adaptive in one or more environments and neutral elsewhere. These 'conditionally neutral' alleles should sweep to fixation when they arise, yet genome comparisons find little evidence for such selective sweeps. Second, genes under biotic selection tend to be of larger effect than genes under abiotic selection. Recent theory suggests this may be a consequence of high gene flow among populations under selection for local adaptation.     

Molecular evolution near a two-locus balanced polymorphism

Balancing selection at one locus can increase the amount of selectively neutral variation within neighboring genomic regions. Discrete phenotypic polymorphisms studied in natural populations are frequently determined by sets of interacting genes instead of alternative alleles at single loci. We extend coalescent theory to investigate balancing selection on combinations of linked genes. We find that variation at neutral sites is increased across a much larger genomic region relative to the single-locus models: the entire region lying between the two loci in balanced combination is affected to some degree. Epistatic selection maintains these high levels of neutral variation because it directly opposes the homogenizing effect of recombination. The results of the theory are discussed in relation to published gene sequence data, primarily from Drosophila.     

How does adaptation sweep through the genome? Insights from long-term selection experiments

A major goal in evolutionary biology is to understand the origins and fates of adaptive mutations. Natural selection may act to increase the frequency of de novo beneficial mutations, or those already present in the population as standing genetic variation. These beneficial mutations may ultimately reach fixation in a population, or they may stop increasing in frequency once a particular phenotypic state has been achieved. It is not yet well understood how different features of population biology, and/or different environmental circumstances affect these adaptive processes. Experimental evolution is a promising technique for studying the dynamics of beneficial alleles, as populations evolving in the laboratory experience natural selection in a replicated, controlled manner. Whole-genome sequencing, regularly obtained over the course of sustained laboratory selection, could potentially reveal insights into the mutational dynamics that most likely occur in natural populations under similar circumstances. To date, only a few evolution experiments for which whole-genome data are available exist. This review describes results from these resequenced laboratory-selected populations, in systems with and without sexual recombination. In asexual systems, adaptation from new mutations can be studied, and results to date suggest that the complete, unimpeded fixation of these mutations is not always observed. In sexual systems, adaptation from standing genetic variation can be studied, and in the admittedly few examples we have, the complete fixation of standing variants is not always observed. To date, the relative frequency of adaptation from new mutations versus standing variation has not been tested using a single experimental system, but recent studies using Caenorhabditis elegans and Saccharomyces cerevisiae suggest that this a realistic future goal.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.     

The coalescent process in models with selection, recombination and geographic subdivision

A population genetic model with a single locus at which balancing selection acts and many linked loci at which neutral mutations can occur is analysed using the coalescent approach. The model incorporates geographic subdivision with migration, as well as mutation, recombination, and genetic drift of neutral variation. It is found that geographic subdivision can affect genetic variation even with high rates of migration, providing that selection is strong enough to maintain different allele frequencies at the selected locus. Published sequence data from the alcohol dehydrogenase locus of Drosophila melanogaster are found to fit the proposed model slightly better than a similar model without subdivision.     

Compensatory Evolution and the Origins of Innovations

Cryptic genetic sequences have attenuated effects on phenotypes. In the classic view, relaxed selection allows cryptic genetic diversity to build up across individuals in a population, providing alleles that may later contribute to adaptation when co-opted—e.g., following a mutation increasing expression from a low, attenuated baseline. This view is described, for example, by the metaphor of the spread of a population across a neutral network in genotype space. As an alternative view, consider the fact that most phenotypic traits are affected by multiple sequences, including cryptic ones. Even in a strictly clonal population, the co-option of cryptic sequences at different loci may have different phenotypic effects and offer the population multiple adaptive possibilities. Here, we model the evolution of quantitative phenotypic characters encoded by cryptic sequences and compare the relative contributions of genetic diversity and of variation across sites to the phenotypic potential of a population. We show that most of the phenotypic variation accessible through co-option would exist even in populations with no polymorphism. This is made possible by a history of compensatory evolution, whereby the phenotypic effect of a cryptic mutation at one site was balanced by mutations elsewhere in the genome, leading to a diversity of cryptic effect sizes across sites rather than across individuals. Cryptic sequences might accelerate adaptation and facilitate large phenotypic changes even in the absence of genetic diversity, as traditionally defined in terms of alternative alleles.     

Spatial pattern of nucleotide polymorphism indicates molecular adaptation in the bryophyte Sphagnum fimbriatum

In organisms with haploid-dominant life cycles, natural selection is expected to be especially effective because genetic variation is exposed directly to selection. However, in spore-producing plants with high dispersal abilities, among-population migration may counteract local adaptation by continuously redistributing genetic variability. In this study, we tested for adaptation at the molecular level by comparing nucleotide polymorphism in two genes (GapC and Rpb2) in 10 European populations of the peatmoss species, Sphagnum fimbriatum with variability at nine microsatellite loci assumed to be selectively neutral. In line with previous results, the GapC and Rpb2 genes showed strikingly different patterns of nucleotide polymorphism. Neutrality tests and comparison of population differentiation based on the GapC and Rpb2 genes with neutrally evolving microsatellites using coalescent simulations supported non-neutral evolution in GapC, but neutral evolution in the Rpb2 gene. These observations and the positions of the replacement mutations in the GAPDH enzyme (coded by GapC) indicate a significant impact of replacement mutations on enzyme function. Furthermore, the geographic distribution of alternate GapC alleles and/or linked genomic regions suggests that they have had differential success in the recolonization of Europe following the Last Glacial Maximum.     

Perspective: Sign epistasis and genetic constraint on evolutionary trajectories

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

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.     

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.     

Emergent neutrality in adaptive asexual evolution

In nonrecombining genomes, genetic linkage can be an important evolutionary force. Linkage generates interference interactions, by which simultaneously occurring mutations affect each other's chance of fixation. Here, we develop a comprehensive model of adaptive evolution in linked genomes, which integrates interference interactions between multiple beneficial and deleterious mutations into a unified framework. By an approximate analytical solution, we predict the fixation rates of these mutations, as well as the probabilities of beneficial and deleterious alleles at fixed genomic sites. We find that interference interactions generate a regime of emergent neutrality: all genomic sites with selection coefficients smaller in magnitude than a characteristic threshold have nearly random fixed alleles, and both beneficial and deleterious mutations at these sites have nearly neutral fixation rates. We show that this dynamic limits not only the speed of adaptation, but also a population's degree of adaptation in its current environment. We apply the model to different scenarios: stationary adaptation in a time-dependent environment and approach to equilibrium in a fixed environment. In both cases, the analytical predictions are in good agreement with numerical simulations. Our results suggest that interference can severely compromise biological functions in an adapting population, which sets viability limits on adaptive evolution under linkage.     

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, ‘antagonistic’ pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits ‘adaptive’ pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a ‘drought escape’ strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.     

Evidence for contrasting modes of selection at interacting globin genes in the European rabbit (Oryctolagus cuniculus)

In hybrid zones between genetically differentiated populations, variation in locus-specific rates of introgression may reflect adaptation to different environments or adaptation to different genetic backgrounds. The European rabbit, Oryctolagus cuniculus, is well-suited to studies of such hybrid zone dynamics because it is composed of two genetically divergent subspecies that hybridize in a zone of secondary contact in central Iberia. A species-wide survey of allozyme variation revealed a broad range of locus-specific divergence levels (F(ST) ranged from 0 to 0.54, mean F(ST)=0.16). Interestingly, the two loci that fell at opposite ends of the distribution of F(ST) values, haemoglobin alpha-chain (HBA) and haemoglobin beta-chain (HBB), encode interacting subunits of the haemoglobin protein. The contrasting patterns of spatial variation at these two loci could not be reconciled under a neutral model of population structure. The HBA gene exhibited higher-than-expected levels of population differentiation, consistent with a history of spatially varying selection. The HBB gene exhibited lower-than-expected levels of population differentiation, consistent with some form of spatially uniform selection. Patterns of linkage disequilibrium and allele frequency variation do not appear to fit any simple model of two-locus epistatic selection.