Identifying quantitative trait locus by genetic background interactions in association studies

Association studies are designed to identify main effects of alleles across a potentially wide range of genetic backgrounds. To control for spurious associations, effects of the genetic background itself are often incorporated into the linear model, either in the form of subpopulation effects in the case of structure or in the form of genetic relationship matrices in the case of complex pedigrees. In this context epistatic interactions between loci can be captured as an interaction effect between the associated locus and the genetic background. In this study I developed genetic and statistical models to tie the locus by genetic background interaction idea back to more standard concepts of epistasis when genetic background is modeled using an additive relationship matrix. I also simulated epistatic interactions in four-generation randomly mating pedigrees and evaluated the ability of the statistical models to identify when a biallelic associated locus was epistatic to other loci. Under additive-by-additive epistasis, when interaction effects of the associated locus were quite large (explaining 20% of the phenotypic variance), epistasis was detected in 79% of pedigrees containing 320 individuals. The epistatic model also predicted the genotypic value of progeny better than a standard additive model in 78% of simulations. When interaction effects were smaller (although still fairly large, explaining 5% of the phenotypic variance), epistasis was detected in only 9% of pedigrees containing 320 individuals and the epistatic and additive models were equally effective at predicting the genotypic values of progeny. Epistasis was detected with the same power whether the overall epistatic effect was the result of a single pairwise interaction or the sum of nine pairwise interactions, each generating one ninth of the epistatic variance. The power to detect epistasis was highest (94%) at low QTL minor allele frequency, fell to a minimum (60%) at minor allele frequency of about 0.2, and then plateaued at about 80% as alleles reached intermediate frequencies. The power to detect epistasis declined when the linkage disequilibrium between the DNA marker and the functional polymorphism was not complete.     

Little epistasis for anxiety-related measures in the DeFries strains of laboratory mice

Recent advances in methodologies for testing epistatic interactions, combined with several successes in demonstrating genetic interaction effects in animal and human genetics, have rekindled interest in the role of epistatic influences on complex traits. It has even been suggested that the unacknowledged presence of epistasis vitiates the genetic dissection of human and animal behavior. Here we report a genome-wide interaction analysis of 1636 F₂ mice to show that epistasis is of minimal importance in an animal model of anxiety. By using a sufficiently large sample of F₂ animals, we provide evidence that interaction effects between any two loci contribute less than 5% to the total phenotypic variance in multiple tests of anxiety. We conclude that interactions between loci do not necessarily vitiate the genetic analysis of behavior in at least one animal model of anxiety.     

Genetic expectations of quantitative trait loci main and interaction effects obtained with the triple testcross design and their relevance for the analysis of heterosis

Interpretation of experimental results from quantitative trait loci (QTL) mapping studies on the predominant type of gene action can be severely affected by the choice of statistical model, experimental design, and provision of epistasis. In this study, we derive quantitative genetic expectations of (i) QTL effects obtained from one-dimensional genome scans with the triple testcross (TTC) design and (ii) pairwise interactions between marker loci using two-way analyses of variance (ANOVA) under the F(2)- and the F(infinity)-metric model. The theoretical results show that genetic expectations of QTL effects estimated with the TTC design are complex, comprising both main and epistatic effects, and that genetic expectations of two-way marker interactions are not straightforward extensions of effects estimated in one-dimensional scans. We also demonstrate that the TTC design can partially overcome the limitations of the design III in separating QTL main effects and their epistatic interactions in the analysis of heterosis and that dominance x additive epistatic interactions of individual QTL with the genetic background can be estimated with a one-dimensional genome scan. Furthermore, we present genetic expectations of variance components for the analysis of TTC progeny tested in a split-plot design, assuming digenic epistasis and arbitrary linkage.     

Replication and explorations of high-order epistasis using a large advanced intercross line pedigree

Dissection of the genetic architecture of complex traits persists as a major challenge in biology; despite considerable efforts, much remains unclear including the role and importance of genetic interactions. This study provides empirical evidence for a strong and persistent contribution of both second- and third-order epistatic interactions to long-term selection response for body weight in two divergently selected chicken lines. We earlier reported a network of interacting loci with large effects on body weight in an F(2) intercross between these high- and low-body weight lines. Here, most pair-wise interactions in the network are replicated in an independent eight-generation advanced intercross line (AIL). The original report showed an important contribution of capacitating epistasis to growth, meaning that the genotype at a hub in the network releases the effects of one or several peripheral loci. After fine-mapping of the loci in the AIL, we show that these interactions were persistent over time. The replication of five of six originally reported epistatic loci, as well as the capacitating epistasis, provides strong empirical evidence that the originally observed epistasis is of biological importance and is a contributor in the genetic architecture of this population. The stability of genetic interaction mechanisms over time indicates a non-transient role of epistasis on phenotypic change. Third-order epistasis was for the first time examined in this study and was shown to make an important contribution to growth, which suggests that the genetic architecture of growth is more complex than can be explained by two-locus interactions only. Our results illustrate the importance of designing studies that facilitate exploration of epistasis in populations for obtaining a comprehensive understanding of the genetics underlying a complex trait.     

A penalized maximum likelihood method for estimating epistatic effects of QTL

Although epistasis is an important phenomenon in the genetics and evolution of complex traits, epistatic effects are hard to estimate. The main problem is due to the overparameterized epistatic genetic models. An epistatic genetic model should include potential pair-wise interaction effects of all loci. However, the model is saturated quickly as the number of loci increases. Therefore, a variable selection technique is usually considered to exclude those interactions with negligible effects. With such techniques, we may run a high risk of missing some important interaction effects by not fully exploring the extremely large parameter space of models. We develop a penalized maximum likelihood method. The method developed here adopts a penalty that depends on the values of the parameters. The penalized likelihood method allows spurious QTL effects to be shrunk towards zero, while QTL with large effects are estimated with virtually no shrinkage. A simulation study shows that the new method can handle a model with a number of effects 15 times larger than the sample size. Simulation studies also show that results of the penalized likelihood method are comparable to the Bayesian shrinkage analysis, but the computational speed of the penalized method is orders of magnitude faster.     

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.     

Epistasis and the mutation load: a measurement-theoretical approach

An approximate solution for the mean fitness in mutation-selection balance with arbitrary order of epistatic interaction is derived. The solution is based on the assumptions of coupling equilibrium and that the interaction effects are multilinear. We find that the effect of m-order epistatic interactions (i.e., interactions among groups of m loci) on the load is dependent on the total genomic mutation rate, U, to the mth power. Thus, higher-order gene interactions are potentially important if U is large and the interaction density among loci is not too low. The solution suggests that synergistic epistasis will decrease the mutation load and that variation in epistatic effects will elevate the load. Both of these results, however, are strictly true only if they refer to epistatic interaction strengths measured in the optimal genotype. If gene interactions are measured at mutation-selection equilibrium, only synergistic interactions among even numbers of genes will reduce the load. Odd-ordered synergistic interactions will then elevate the load. There is no systematic relationship between variation in epistasis and load at equilibrium. We argue that empirical estimates of gene interaction must pay attention to the genetic background in which the effects are measured and that it may be advantageous to refer to average interaction intensities as measured in mutation-selection equilibrium. We derive a simple criterion for the strength of epistasis that is necessary to overcome the twofold disadvantage of sex.     

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.     

Mapping the epistatic network underlying murine reproductive fatpad variation

Genome-wide mapping analyses are now commonplace in many species and several networks of interacting loci have been reported. However, relatively few details regarding epistatic interactions and their contribution to complex trait variation in multicellular organisms are available and the identification of positional candidate loci for epistatic QTL (epiQTL) is hampered, especially in mammals, by the limited genetic resolution inherent in most study designs. Here we further investigate the genetic architecture of reproductive fatpad weight in mice using the F(10) generation of the LG,SM advanced intercross (AI) line. We apply multiple mapping techniques including a single-locus model, locus-specific composite interval mapping (CIM), and tests for multiple QTL per chromosome to the 12 chromosomes known to harbor single-locus QTL (slQTL) affecting obesity in this cross. We also perform a genome-wide scan for pairwise epistasis. Using this combination of approaches we detect 199 peaks spread over all 19 autosomes, which potentially contribute to trait variation including all eight original F(2) loci (Adip1-8), novel slQTL peaks on chromosomes 7 and 9, and several novel epistatic loci. Extensive epistasis is confirmed involving both slQTL confidence intervals (C.I.) as well as regions that show no significant additive or dominance effects. These results provide important new insights into mapping complex genetic architectures and the role of epistasis in complex trait variation.     

Evolution of genetic architecture under directional selection

We investigate the multilinear epistatic model under mutation-limited directional selection. We confirm previous results that only directional epistasis, in which genes on average reinforce or diminish each other's effects, contribute to the initial evolution of mutational effects. Thus, either canalization or decanalization can occur under directional selection, depending on whether positive or negative epistasis is prevalent. We then focus on the evolution of the epistatic coefficients themselves. In the absence of higher-order epistasis, positive pairwise epistasis will tend to weaken relative to additive effects, while negative pairwise epistasis will tend to become strengthened. Positive third-order epistasis will counteract these effects, while negative third-order epistasis will reinforce them. More generally, gene interactions of all orders have an inherent tendency for negative changes under directional selection, which can only be modified by higher-order directional epistasis. We identify three types of nonadditive quasi-equilibrium architectures that, although not strictly stable, can be maintained for an extended time: (1) nondirectional epistatic architectures; (2) canalized architectures with strong epistasis; and (3) near-additive architectures in which additive effects keep increasing relative to epistasis.     

Signals of demographic expansion in Drosophila virilis

Background

The antagonistic co-evolution of hosts and their parasites is considered to be a potential driving force in maintaining host genetic variation including sexual reproduction and recombination. The examination of this hypothesis calls for information about the genetic basis of host-parasite interactions – such as how many genes are involved, how big an effect these genes have and whether there is epistasis between loci. We here examine the genetic architecture of quantitative resistance in animal and plant hosts by concatenating published studies that have identified quantitative trait loci (QTL) for host resistance in animals and plants.

Results

Collectively, these studies show that host resistance is affected by few loci. We particularly show that additional epistatic interactions, especially between loci on different chromosomes, explain a majority of the effects. Furthermore, we find that when experiments are repeated using different host or parasite genotypes under otherwise identical conditions, the underlying genetic architecture of host resistance can vary dramatically – that is, involves different QTLs and epistatic interactions. QTLs and epistatic loci vary much less when host and parasite types remain the same but experiments are repeated in different environments.

Conclusion

This pattern of variability of the genetic architecture is predicted by strong interactions between genotypes and corroborates the prevalence of varying host-parasite combinations over varying environmental conditions. Moreover, epistasis is a major determinant of phenotypic variance for host resistance. Because epistasis seems to occur predominantly between, rather than within, chromosomes, segregation and chromosome number rather than recombination via cross-over should be the major elements affecting adaptive change in host resistance.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Epistatic interaction maps relative to multiple metabolic phenotypes

An epistatic interaction between two genes occurs when the phenotypic impact of one gene depends on another gene, often exposing a functional association between them. Due to experimental scalability and to evolutionary significance, abundant work has been focused on studying how epistasis affects cellular growth rate, most notably in yeast. However, epistasis likely influences many different phenotypes, affecting our capacity to understand cellular functions, biochemical networks adaptation, and genetic diseases. Despite its broad significance, the extent and nature of epistasis relative to different phenotypes remain fundamentally unexplored. Here we use genome-scale metabolic network modeling to investigate the extent and properties of epistatic interactions relative to multiple phenotypes. Specifically, using an experimentally refined stoichiometric model for Saccharomyces cerevisiae, we computed a three-dimensional matrix of epistatic interactions between any two enzyme gene deletions, with respect to all metabolic flux phenotypes. We found that the total number of epistatic interactions between enzymes increases rapidly as phenotypes are added, plateauing at approximately 80 phenotypes, to an overall connectivity that is roughly 8-fold larger than the one observed relative to growth alone. Looking at interactions across all phenotypes, we found that gene pairs interact incoherently relative to different phenotypes, i.e. antagonistically relative to some phenotypes and synergistically relative to others. Specific deletion-deletion-phenotype triplets can be explained metabolically, suggesting a highly informative role of multi-phenotype epistasis in mapping cellular functions. Finally, we found that genes involved in many interactions across multiple phenotypes are more highly expressed, evolve slower, and tend to be associated with diseases, indicating that the importance of genes is hidden in their total phenotypic impact. Our predictions indicate a pervasiveness of nonlinear effects in how genetic perturbations affect multiple metabolic phenotypes. The approaches and results reported could influence future efforts in understanding metabolic diseases and the role of biochemical regulation in the cell.     

Epistasis interaction of QTL effects as a genetic parameter influencing estimation of the genetic additive effect

Epistasis, an additive-by-additive interaction between quantitative trait loci, has been defined as a deviation from the sum of independent effects of individual genes. Epistasis between QTLs assayed in populations segregating for an entire genome has been found at a frequency close to that expected by chance alone. Recently, epistatic effects have been considered by many researchers as important for complex traits. In order to understand the genetic control of complex traits, it is necessary to clarify additive-by-additive interactions among genes. Herein we compare estimates of a parameter connected with the additive gene action calculated on the basis of two models: a model excluding epistasis and a model with additive-by-additive interaction effects. In this paper two data sets were analysed: 1) 150 barley doubled haploid lines derived from the Steptoe × Morex cross, and 2) 145 DH lines of barley obtained from the Harrington × TR306 cross. The results showed that in cases when the effect of epistasis was different from zero, the coefficient of determination was larger for the model with epistasis than for the one excluding epistasis. These results indicate that epistatic interaction plays an important role in controlling the expression of complex traits.     

Environmental Dependence of Genetic Constraint

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data.     

Effect of varying epistasis on the evolution of recombination

Whether recombination decelerates or accelerates a population's response to selection depends, at least in part, on how fitness-determining loci interact. Realistically, all genomes likely contain fitness interactions both with positive and with negative epistasis. Therefore, it is crucial to determine the conditions under which the potential beneficial effects of recombination with negative epistasis prevail over the detrimental effects of recombination with positive epistasis. Here, we examine the simultaneous effects of diverse epistatic interactions with different strengths and signs in a simplified model system with independent pairs of interacting loci and selection acting only on the haploid phase. We find that the average form of epistasis does not predict the average amount of linkage disequilibrium generated or the impact on a recombination modifier when compared to results using the entire distribution of epistatic effects and associated single-mutant effects. Moreover, we show that epistatic interactions of a given strength can produce very different effects, having the greatest impact when selection is weak. In summary, we observe that the evolution of recombination at mutation-selection balance might be driven by a small number of interactions with weak selection rather than by the average epistasis of all interactions. We illustrate this effect with an analysis of published data of Saccharomyces cerevisiae. Thus to draw conclusions on the evolution of recombination from experimental data, it is necessary to consider the distribution of epistatic interactions together with the associated selection coefficients.     

水稻株高上位性效应和ＱＥ互作效应的ＱＴＬ遗传研究

利用基因混合模型的ＱＴＬ定位方法研究了由籼稻品种ＩＲ６４和粳稻品种Ａzucena杂交衍生的ＤＨ群体在４个环境中的ＱＴＬ上位性效应和环境互作效应,结果表明,上位性是数量性状的重要遗传基础,并揭示了上位性的几个重要特点,所有的ＱＴＬ都参与了上位性效应的形成,６４％的ＱＴＬ还具有本身的加性效应,因此传统方法对ＱＴＬ加性效应的估算会由于上位性的影响而有偏,其他３６％的ＱＴＬ没有本身的加性效应,却参与了４８％的上位性互作用,这些位点可能通过诱发和修饰其他位点而起作用,上位性的特点还包括,经常发现了一个ＱＴＬ与多个ＱＴＬ发生互作;大效应的ＱＴＬ也参与上位性互作;上位性互作受环境影响,ＱＴＬ与环境的互效应比ＱＴＬ的主效应更多地被检测到,表明数量性状基因的表达易受环境影响。     

Mapping epistatic quantitative trait loci with one-dimensional genome searches

The discovery of epistatically interacting QTL is hampered by the intractability and low power to detect QTL in multidimensional genome searches. We describe a new method that maps epistatic QTL by identifying loci of high QTL by genetic background interaction. This approach allows detection of QTL involved not only in pairwise but also higher-order interaction, and does so with one-dimensional genome searches. The approach requires large populations derived from multiple related inbred-line crosses as is more typically available for plants. Using maximum likelihood, the method contrasts models in which QTL allelic values are either nested within, or fixed over, populations. We apply the method to simulated doubled-haploid populations derived from a diallel among three inbred parents and illustrate the power of the method to detect QTL of different effect size and different levels of QTL by genetic background interaction. Further, we show how the method can be used in conjunction with standard two-locus QTL detection models that use two-dimensional genome searches and find that the method may double the power to detect first-order epistasis.     

A General Model for Multilocus Epistatic Interactions in Case-Control Studies

Background

Epistasis, i.e., the interaction of alleles at different loci, is thought to play a central role in the formation and progression of complex diseases. The complexity of disease expression should arise from a complex network of epistatic interactions involving multiple genes.

Methodology

We develop a general model for testing high-order epistatic interactions for a complex disease in a case-control study. We incorporate the quantitative genetic theory of high-order epistasis into the setting of cases and controls sampled from a natural population. The new model allows the identification and testing of epistasis and its various genetic components.

Conclusions

Simulation studies were used to examine the power and false positive rates of the model under different sampling strategies. The model was used to detect epistasis in a case-control study of inflammatory bowel disease, in which five SNPs at a candidate gene were typed, leading to the identification of a significant three-locus epistasis.     

Comparison of strategies to detect epistasis from eQTL data

Genome-wide association studies have been instrumental in identifying genetic variants associated with complex traits such as human disease or gene expression phenotypes. It has been proposed that extending existing analysis methods by considering interactions between pairs of loci may uncover additional genetic effects. However, the large number of possible two-marker tests presents significant computational and statistical challenges. Although several strategies to detect epistasis effects have been proposed and tested for specific phenotypes, so far there has been no systematic attempt to compare their performance using real data. We made use of thousands of gene expression traits from linkage and eQTL studies, to compare the performance of different strategies. We found that using information from marginal associations between markers and phenotypes to detect epistatic effects yielded a lower false discovery rate (FDR) than a strategy solely using biological annotation in yeast, whereas results from human data were inconclusive. For future studies whose aim is to discover epistatic effects, we recommend incorporating information about marginal associations between SNPs and phenotypes instead of relying solely on biological annotation. Improved methods to discover epistatic effects will result in a more complete understanding of complex genetic effects.