THE ACTION OF STABILIZING SELECTION,MUTATION, AND DRIFT ON EPISTATIC QUANTITATIVE TRAITS

For a quantitative trait under stabilizing selection, the effect of epistasis on its genetic architecture and on the changes of genetic variance caused by bottlenecking were investigated using theory and simulation. Assuming empirical estimates of the rate and effects of mutations and the intensity of selection, we assessed the impact of two‐locus epistasis (synergistic/antagonistic) among linked or unlinked loci on the distribution of effects and frequencies of segregating loci in populations at the mutation‐selection‐drift balance. Strong pervasive epistasis did not modify substantially the genetic properties of the trait and, therefore, the most likely explanation for the low amount of variation usually accounted by the loci detected in genome‐wide association analyses is that many causal loci will pass undetected. We investigated the impact of epistasis on the changes in genetic variance components when large populations were subjected to successive bottlenecks of different sizes, considering the action of genetic drift, operating singly (D), or jointly with mutation (MD) and selection (MSD). An initial increase of the different components of the genetic variance, as well as a dramatic acceleration of the between‐line divergence, were always associated with synergistic epistasis but were strongly constrained by selection.     

Epistasis can increase multivariate trait diversity in haploid non-recombining populations

We evaluate the effect of epistasis on genetically-based multivariate trait variation in haploid non-recombining populations. In a univariate setting, past work has shown that epistasis reduces genetic variance (additive plus epistatic) in a population experiencing stabilizing selection. Here we show that in a multivariate setting, epistasis also reduces total genetic variation across the entire multivariate trait in a population experiencing stabilizing selection. But, we also show that the pattern of variation across the multivariate trait can be more even when epistasis occurs compared to when epistasis is absent, such that some character combinations will have more genetic variance when epistasis occurs compared to when epistasis is absent. In fact, a measure of generalized multivariate trait variation can be increased by epistasis under weak to moderate stabilizing selection conditions, as well as neutral conditions. Likewise, a measure of conditional evolvability can be increased by epistasis under weak to moderate stabilizing selection and neutral conditions. We investigate the nature of epistasis assuming a multivariate-normal model genetic effects and investigate the nature of epistasis underlying the biophysical properties of RNA. Increased multivariate diversity occurs for populations that are infinite in size, as well as populations that are finite in size. Our model of finite populations is explicitly genealogical and we link our findings about the evenness of eigenvalues with epistasis to prior work on the genealogical mapping of epistatic effects.     

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.     

EPISTASIS AS A SOURCE OF INCREASED ADDITIVE GENETIC VARIANCE AT POPULATION BOTTLENECKS

The role of epistasis in evolution and speciation has remained controversial. We use a new parameterization of physiological epistasis to examine the effects of epistasis on levels of additive genetic variance during a population bottleneck. We found that all forms of epistasis increase average additive genetic variance in finite populations derived from initial populations with intermediate allele frequencies. Average additive variance continues to increase over many generations, especially at larger population sizes (N = 32 to 64). Additive-by-additive epistasis is the most potent source of additive genetic variance in this situation, whereas dominance-by-dominance epistasis contributes smaller amounts of additive genetic variance. With additive-by-dominance epistasis, additive genetic variance decreases at a relatively high rate immediately after a population bottleneck, rebounding to higher levels after several generations. Empirical examples of epistasis for murine adult body weight based on measured genotypes are provided illustrating the varying effects of epistasis on additive genetic variance during population bottlenecks.     

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions

Abstract We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.     

THE EVOLUTION OF STRONG REPRODUCTIVE ISOLATION

Felsenstein distinguished two ways by which selection can directly strengthen isolation. First, a modifier that strengthens prezygotic isolation can be favored everywhere. This fits with the traditional view of reinforcement as an adaptation to reduce deleterious hybridization by strengthening assortative mating. Second, selection can favor association between different incompatibilities, despite recombination. We generalize this "two allele" model to follow associations among any number of incompatibilities, which may include both assortment and hybrid inviability. Our key argument is that this process, of coupling between incompatibilities, may be quite different from the usual view of reinforcement: strong isolation can evolve through the coupling of any kind of incompatibility, whether prezygotic or postzygotic. Single locus incompatibilities become coupled because associations between them increase the variance in compatibility, which in turn increases mean fitness if there is positive epistasis. Multiple incompatibilities, each maintained by epistasis, can become coupled in the same way. In contrast, a single-locus incompatibility can become coupled with loci that reduce the viability of haploid hybrids because this reduces harmful recombination. We obtain simple approximations for the limits of tight linkage, and strong assortment, and show how assortment alleles can invade through associations with other components of reproductive isolation.     

Mapping quantitative trait loci controlling silking date in a diallel cross among four lines of maize

 We describe and apply an interval mapping method for quantitative trait locus (QTL) detection using F₃ and testcross progenies derived from F₂ populations obtained from a diallel cross among four elite lines of maize. Linear model-based procedures were used for the test and estimation of putative QTL effects together with genetic interactions including epistasis. We mapped QTL associated with silking date and explored their genetic effects. Ten QTL were detected, and these explained more than 40% of the phenotypic variance. Most of these QTL had consistent and stable effects among genetic backgrounds and did not show significant epistasis. QTL-by-environment interaction was important for four QTL and was essentially due to changes in magnitude of allelic effects. These results show the efficiency of our method in several genetic situations as well as the power of the diallel design in detecting QTL simultaneously over several populations. Received: 2 September 1996 / Accepted: 20 December 1996     

Learning genetic epistasis using Bayesian network scoring criteria

Background  

Gene-gene epistatic interactions likely play an important role in the genetic basis of many common diseases. Recently, machine-learning and data mining methods have been developed for learning epistatic relationships from data. A well-known combinatorial method that has been successfully applied for detecting epistasis is Multifactor Dimensionality Reduction (MDR). Jiang et al. created a combinatorial epistasis learning method called BNMBL to learn Bayesian network (BN) epistatic models. They compared BNMBL to MDR using simulated data sets. Each of these data sets was generated from a model that associates two SNPs with a disease and includes 18 unrelated SNPs. For each data set, BNMBL and MDR were used to score all 2-SNP models, and BNMBL learned significantly more correct models. In real data sets, we ordinarily do not know the number of SNPs that influence phenotype. BNMBL may not perform as well if we also scored models containing more than two SNPs. Furthermore, a number of other BN scoring criteria have been developed. They may detect epistatic interactions even better than BNMBL.     

Association Mapping for Epistasis and Environmental Interaction of Yield Traits in 323 Cotton Cultivars under 9 Different Environments

Improving yield is a major objective for cotton breeding schemes, and lint yield and its three component traits (boll number, boll weight and lint percentage) are complex traits controlled by multiple genes and various environments. Association mapping was performed to detect markers associated with these four traits using 651 simple sequence repeats (SSRs). A mixed linear model including epistasis and environmental interaction was used to screen the loci associated with these four yield traits by 323 accessions of Gossypium hirsutum L. evaluated in nine different environments. 251 significant loci were detected to be associated with lint yield and its three components, including 69 loci with individual effects and all involved in epistasis interactions. These significant loci explain ∼ 62.05% of the phenotypic variance (ranging from 49.06% ∼ 72.29% for these four traits). It was indicated by high contribution of environmental interaction to the phenotypic variance for lint yield and boll numbers, that genetic effects of SSR loci were susceptible to environment factors. Shared loci were also observed among these four traits, which may be used for simultaneous improvement in cotton breeding for yield traits. Furthermore, consistent and elite loci were screened with −Log₁₀ (P-value) >8.0 based on predicted effects of loci detected in different environments. There was one locus and 6 pairs of epistasis for lint yield, 4 loci and 10 epistasis for boll number, 15 loci and 2 epistasis for boll weight, and 2 loci and 5 epistasis for lint percentage, respectively. These results provided insights into the genetic basis of lint yield and its components and may be useful for marker-assisted breeding to improve cotton production.     

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.     

Efficiency of triple test cross for detecting epistasis with marker information

The triple test cross (TTC) is an experimental design for detecting epistasis and estimating the components of genetic variance for quantitative traits. In this paper, we extend the analysis to include molecular information. The statistical power of the mating design was assessed under a model assuming that a finite number of loci affect the trait in question. Formulae are developed for the analysis with or without marker information relating to the recombination fraction between loci, the genetical properties of quantitative trait controlled by the quantitative trait loci (QTL), the linkage phases of the parents and population size. Application of these formulae showed that the recombination fraction between genes and the magnitude and the types of epistasis have important interactions in their effects on power. The results demonstrate that the TTC may have increased power to detect epistasis when marker information is present. However, the simulation experiments show that the standard deviation of the estimated expected mean square was higher with one marker than that with two, whereas the corresponding value without marker information was the lowest. In addition, we demonstrate that the relative position of QTL and markers and the number of markers can both affect the power of epistatic detection.     

EPISTATIC INTERACTIONS DETERMINE THE MUTATIONAL PATHWAYS AND COEXISTENCE OF LINEAGES IN CLONAL ESCHERICHIA COLI POPULATIONS

Understanding how diversity emerges in a single niche is not fully understood. Rugged fitness landscapes and epistasis between beneficial mutations could explain coexistence among emerging lineages. To provide an experimental test of this notion, we investigated epistasis among four pleiotropic mutations in rpoS, mglD, malT, and hfq present in two coexisting lineages that repeatedly fixed in experimental populations of Escherichia coli. The mutations were transferred into the ancestral background individually or in combination of double or triple alleles. The combined competitive fitness of two or three beneficial mutations from the same lineage was consistently lower than the sum of the competitive fitness of single mutants—a clear indication of negative epistasis within lineages. We also found sign epistasis (i.e., the combined fitness of two beneficial mutations lower than the ancestor), not only from two different lineages (i.e., hfq and rpoS) but also from the same lineage (i.e., mglD and malT). The sign epistasis between loci of different lineages indeed indicated a rugged fitness landscape, providing an epistatic explanation for the coexistence of distinct rpoS and hfq lineages in evolving populations. The negative and sign epistasis between beneficial mutations within the same lineage can further explain the order of mutation acquisition.     

Incomplete loss of a conserved trait: function,latitudinal cline,and genetic constraints

Retention of nonfunctional traits over evolutionary time is puzzling, because the cost of trait production should drive loss. Indeed, several studies have found nonfunctional traits are rapidly eliminated by selection. However, theory suggests that complex genetic interactions and a lack of genetic variance can constrain evolution, including trait loss. In the mustard family Brassicaceae the conserved floral condition includes four long and two short stamens, but we show that short stamens in the highly self‐pollinating mustard Arabidopsis thaliana do not significantly increase selfed seed set, suggesting that the trait has lost most or all of its function after the transition to selfing. We find that short stamen loss is common in native populations. Loss is incomplete and decreases with increasing latitude, a cline unexplained by correlations with flowering time or ovule count (which also vary with latitude). Using recombinant inbred lines derived from a cross between plants at the latitudinal extremes of the native range, we found three QTLs affecting short stamen number, with epistasis among them constraining stamen loss. Constraints on stamen loss from both epistasis and low genetic variance may be augmented by high selfing rates, suggesting that these kinds of constraints may be common in inbred species.     

Two-stage two-locus models in genome-wide association

Studies in model organisms suggest that epistasis may play an important role in the etiology of complex diseases and traits in humans. With the era of large-scale genome-wide association studies fast approaching, it is important to quantify whether it will be possible to detect interacting loci using realistic sample sizes in humans and to what extent undetected epistasis will adversely affect power to detect association when single-locus approaches are employed. We therefore investigated the power to detect association for an extensive range of two-locus quantitative trait models that incorporated varying degrees of epistasis. We compared the power to detect association using a single-locus model that ignored interaction effects, a full two-locus model that allowed for interactions, and, most important, two two-stage strategies whereby a subset of loci initially identified using single-locus tests were analyzed using the full two-locus model. Despite the penalty introduced by multiple testing, fitting the full two-locus model performed better than single-locus tests for many of the situations considered, particularly when compared with attempts to detect both individual loci. Using a two-stage strategy reduced the computational burden associated with performing an exhaustive two-locus search across the genome but was not as powerful as the exhaustive search when loci interacted. Two-stage approaches also increased the risk of missing interacting loci that contributed little effect at the margins. Based on our extensive simulations, our results suggest that an exhaustive search involving all pairwise combinations of markers across the genome might provide a useful complement to single-locus scans in identifying interacting loci that contribute to moderate proportions of the phenotypic variance.     

WILL POPULATION BOTTLENECKS AND MULTILOCUS EPISTASIS INCREASE ADDITIVE GENETIC VARIANCE?

We apply new analytical methods to understand the consequences of population bottlenecks for expected additive genetic variance. We analyze essentially all models for multilocus epistasis that have been numerically simulated to demonstrate increased additive variance. We conclude that for biologically plausible models, large increases in expected additive variance--attributable to epistasis rather than dominance--are unlikely. Naciri-Graven and Goudet (2003) found that as the number of epistatically interacting loci increases, additive variance tends to be inflated more after a bottleneck. We argue that this result reflects biologically unrealistic aspects of their models. Specifically, as the number of loci increases, higher-order epistatic interactions become increasingly important in these models, with an increasing fraction of the genetic variance becoming nonadditive, contrary to empirical observations. As shown by Barton and Turelli (2004), without dominance, conversion of nonadditive to additive variance depends only on the variance components and not on the number of loci per se. Numerical results indicating that more inbreeding is needed to produce maximal release of additive variance with more loci follow directly from our analytical results, which show that high levels of inbreeding (F > 0.5) are needed for significant conversion of higher-order components. We discuss alternative approaches to modeling multilocus epistasis and understanding its consequences.     

The Evolution of Epistasis and Its Links With Genetic Robustness, Complexity and Drift in a Phenotypic Model of Adaptation

The epistatic interactions among mutations have a large effect on the evolution of populations. In this article we provide a formalism under which epistatic interactions among pairs of mutations have a distribution whose mean can be modulated. We find that the mean epistasis is correlated to the effect of mutations or genetic robustness, which suggests that such formalism is in good agreement with most in silico models of evolution where the same pattern is observed. We further show that the evolution of epistasis is highly dependant on the intensity of drift and of how complex the organisms are, and that either positive or negative epistasis could be selected for, depending on the balance between the efficiency of selection and the intensity of drift.     

Diallel and generation means analyses for the components of resistance to Cercospora arachidicola in peanut

Summary The inheritance of the components of partial resistance to Cercospora arachidicola Hori in peanut (Arachis hypogaea L.) was examined in two five-parent diallels and in the six generations of two single crosses in greenhouse tests. The Griffing (1956) analysis indicated general combining ability (GCA) to be of most importance, yet large ratios of SCA/GCA sum of squares suggested nonadditive genetic variance as well. Reciprocal effects were found for lesion area and lesion number/10 cm² leaf area. The importance of nonadditive genetic variance was substantiated by the lack of fit for the additive-dominance model in the Hayman's analysis (1954 a, b). Further evidence from the Hayman's analysis indicated that epistasis may be important in determining the inheritance of some of the components of resistance. Additive gene effects alone accounted for the genetic variability observed among the generation means from two single crosses for all components of resistance except latent period. There was evidence that epistasis was an important mode of gene action for the inheritance of latent period.Paper No. 10172 of the Journal Series of the North Carolina Agricultural Research Service, Raleigh, NC 27601, USA     

EPISTASIS AND THE EVOLUTION OF ADDITIVE GENETIC VARIANCE IN POPULATIONS THAT PASS THROUGH A BOTTLENECK

Traditional models of genetic drift predict a linear decrease in additive genetic variance for populations passing through a bottleneck. This perceived lack of heritable variance limits the scope of founder-effect models of speciation. We produced 55 replicate bottleneck populations maintained at two male-female pairs through four generations of inbreeding (average F = 0.39). These populations were formed from an F₂ intercross of the LG/J and SM/J inbred mouse strains. Two contemporaneous control strains maintained with more than 60 mating pairs per generation were formed from this same source population. The average level of within-strain additive genetic variance for adult body weight was compared between the control and experimental lines. Additive genetic variance for adult body weight within experimental bottleneck strains was significantly higher than expected under an additive genetic model This enhancement of additive genetic variance under inbreeding is likely to be due to epistasis, which retards or reverses the loss of additive genetic variance under inbreeding for adult body weight in this population. Therefore, founder-effect speciation processes may not be constrained by a loss of heritable variance due to population bottlenecks.     

ON THE EFFECT OF FOUNDER EVENTS ON EPISTATIC GENETIC VARIANCE

Mayr (1963) proposed that small isolated propagules from a large panmictic population would occasionally undergo a genetic revolution due to loss of genetic variability. More recently Templeton (1980a) has suggested that founder events may be much more important in systems that have strong epistasis. Because of the work of these and other authors it becomes an interesting theoretical problem to study the distribution of epistatic variance in a population following a founder event. In the model presented here measures of coancestry (Cockerham, 1967, 1984; Cockerham and Weir, 1973; Weir and Cockerham, 1973, 1977; Tachida and Cockerham, unpubl.) are used to examine the effect of founder events on additive-by-additive epistasis. Using this approach, the coancestries, or intraclass correlations, within individuals and within demes, together with the genetic variance components in the ancestral population are used to obtain the variance within and among demes following a founder event. Examples are analyzed for single founder events of 1–25 individuals and multiple founder events of two individuals. Following a single founder event, the contribution of the additive variance to the variance within demes relative to the additive variance in the ancestral population is always less than one. However, the contribution of epistatic variance to the variance within demes relative to the epistatic variance in the ancestral population is always greater than one. Thus, while a founder event decreases the contribution of additive variance to the variance within demes, it increases the contribution of epistatic variance to the variance within demes. The contribution of epistatic variance to the variance among demes following a single founder event is not qualitatively different from the contribution of additive variance to the variance among demes. These results indicate that epistatic variance is less likely than additive variance to cause a genetic revolution following a single founder event. When populations undergo multiple founder events the situation changes considerably. Epistatic variance may contribute as much as four times its original value to the variance among demes, while additive variance can contribute maximally twice its original value to the variance among demes. Thus, epistasis, which is relatively unimportant following a single founder event, may have major evolutionary implications if drift is allowed to continue for several generations.