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.     

Epistatic effect between ACACA and FABP2 gene on abdominal fat traits in broilers

Epistasis is generally defined as the interaction between two or more genes or their mRNA or protein products to influence a single trait. Experimental evidence suggested that epistasis could be important in the determination of the genetic architecture of complex traits in domestic animals. Acetyl-coenzyme A carboxylase alpha (ACACA) and fatty acid binding protein 2 (FABP2) are both key factors of lipogenesis and transport. They may play a crucial role in the weight variability of abdominal adipose tissue in the growing chicken. In this study, the polymorphisms of c.2292GA in ACACA and c.-561AC in FABP2 were detected among individuals from two broiler lines which were divergently selected for abdominal fat content. Epistasis between the two SNPs on abdominal fat weight (AFW) and abdominal fat percentage (AFP) was analyzed. The additive × additive epistatic components between these two SNPs were found significant or suggestively significant on both AFW and AFP in lean lines of the 9th and 10th generation; whereas, it was not significantly associated with either AFW or AFP in fat lines. At the same time, there were not any other significant epistatic components found in both generations or in both lines. Significant epistatic effects between these two SNPs found only in the lean lines could partly be due to the fact that the abdominal fat traits in these two experimental lines have been greatly modified by strong artificial selection. The results suggested that the epistasis mode may be different between the lean and fat chicken lines. Our results could be helpful in further understanding the genetic interaction between candidate genes contributing to phenotypic variation of abdominal fat content in broilers.     

High throughput analyses of epistasis for swine body dimensions and organ weights

High throughput analyses were performed to detect epistatic QTL in 17 body dimension and organ weight traits from a large F₂ pig population derived from a White Duroc and Erhualian intercross. The analyses used a nested test framework to handle multiple tests and a combined search algorithm to map epistatic QTL with empirical genome‐wide thresholds derived via prior permutation. Alternative statistical models (e.g. including vs. excluding carcass weight as a covariate) were tested to develop an in‐depth understanding of the role of epistasis in these kinds of traits. Epistasis signals were detected in only two or three traits under each statistical model studied. The interaction component of each pair of epistatic QTL explained a small proportion (0.7 to 2.1%) of the phenotypic variance in general. About half of the detected epistatic QTL pairs involved one of the two major QTL on porcine chromosomes 7 and 4. In those traits, the Erhualian allele consistently increased the phenotypes for the chromosome 7 QTL but decreased them for the chromosome 4 QTL. Models including carcass weight as covariate detected epistasis in body dimension traits whereas those excluding carcass weight found epistasis in organ weight traits. In addition, the epistasis results suggested that a QTL on chromosome 14 could be important for a number of organ weight traits. Using the high‐throughput analysis tool to examine different statistical models was essential for the generation of a complete picture of epistasis in a whole category of traits.     

Genetic dissection of chromosome substitution lines of cotton to discover novel Gossypium barbadense L. alleles for improvement of agronomic traits

We recently released a set of 17 chromosome substitution (CS-B) lines (2n = 52) that contain Gossypium barbadense L. doubled-haploid line ‘3-79’ germplasm systematically introgressed into the Upland inbred ‘TM-1’ of G. hirsutum (L.). TM-1 yields much more than 3-79, but cotton from the latter has superior fiber properties. To explore the use of these quasi-isogenic lines in studying gene interactions, we created a partial diallel among six CS-B lines and the inbred TM-1, and characterized their descendents for lint percentage, boll weight, seedcotton yield and lint yield across four environments. Phenotypic data on the traits were analyzed according to the ADAA genetic model to detect significant additive, dominance, and additive-by-additive epistasis effects at the chromosome and chromosome-by-chromosome levels of CS-B lines. For example, line 3-79 had the lowest boll weight, seedcotton yield and lint yield, but CS-B22Lo homozygous dominance genetic effects on seedcotton and lint yield were nearly four times those of TM-1, and its hybrids with TM-1 had the highest additive-by-additive epistatic effects on seedcotton and lint yield. CS-B14sh, 17, 22Lo and 25 produced positive homozygous dominance effects on lint yield, whereas doubly heterozygous combinations of CS-B14sh with CS-B17, 22Lo and 25 produced negative dominance effects, suggesting that epistatic effects between genes in these chromosomes strongly affect lint yield. The results underscore the opportunities to systematically identify genomic regions harboring genes that impart agronomically significant effects via epistatic interactions. The chromosome-by-chromosome approach significantly complements other strategies to detect and quantify epistatic interaction effects, and the quasi-isogenic nature of families and lines from CS-B intermatings will facilitate high-resolution localization, development of markers for selection and map-assisted identification of genes involved in strong epistatic effects.     

Unraveling epistasis with triple testcross progenies of near-isogenic lines

Libraries of near-isogenic lines (NILs) are a powerful plant genetic resource to map quantitative trait loci (QTL). Nevertheless, QTL mapping with NILs is mostly restricted to genetic main effects. Here we propose a two-step procedure to map additive-by-additive digenic epistasis with NILs. In the first step, a generation means analysis of parents, their F₁ hybrid, and one-segment NILs and their triple testcross (TTC) progenies is used to identify in a one-dimensional scan loci exhibiting QTL-by-background interactions. In a second step, one-segment NILs with significant additive-by-additive background interactions are used to produce particular two-segment NILs to test for digenic epistatic interactions between these segments. We evaluated our approach by analyzing a random subset of a genomewide Arabidopsis thaliana NIL library for growth-related traits. The results of our experimental study illustrated the potential of the presented two-step procedure to map additive-by-additive digenic epistasis with NILs. Furthermore, our findings suggested that additive main effects as well as additive-by-additive digenic epistasis strongly influence the genetic architecture underlying growth-related traits of A. thaliana.     

Statistical power for detecting epistasis QTL effects under the F-2 design

Epistasis refers to gene interaction effect involving two or more genes. Statistical methods for mapping quantitative trait loci (QTL) with epistasis effects have become available recently. However, little is known about the statistical power and sample size requirements for mapping epistatic QTL using genetic markers. In this study, we developed analytical formulae to calculate the statistical power and sample requirement for detecting each epistasis effect under the F-2 design based on crossing inbred lines. Assuming two unlinked interactive QTL and the same absolute value for all epistasis effects, the heritability of additive × additive (a × a) effect is twice as large as that of additive × dominance (a × d) or dominance × additive (d × a) effect, and is four times as large as that of dominance × dominance (d × d) effect. Consequently, among the four types of epistasis effects involving two loci, ''a × a'' effect is the easiest to detect whereas ''d × d'' effect is the most difficult to detect. The statistical power for detecting ''a × a'' effect is similar to that for detecting dominance effect of a single QTL. The sample size requirements for detecting ''a × d'', ''d × a'' and ''d × d'' are highly sensitive to increased distance between the markers and the interacting QTLs. Therefore, using dense marker coverage is critical to detecting those effects.     

Both additivity and epistasis control the genetic variation for fruit quality traits in tomato

The effect of a gene involved in the variation of a quantitative trait may change due to epistatic interactions with the overall genetic background or with other genes through digenic interactions. The classical populations used to map quantitative trait loci (QTL) are poorly efficient to detect epistasis. To assess the importance of epistasis in the genetic control of fruit quality traits, we compared 13 tomato lines having the same genetic background except for one to five chromosome fragments introgressed from a distant line. Six traits were assessed: fruit soluble solid content, sugar content and titratable acidity, fruit weight, locule number and fruit firmness. Except for firmness, a large part of the variation of the six traits was under additive control, but interactions between QTL leading to epistasis effects were common. In the lines cumulating several QTL regions, all the significant epistatic interactions had a sign opposite to the additive effects, suggesting less than additive epistasis. Finally the re-examination of the segregating population initially used to map the QTL confirmed the extent of epistasis, which frequently involved a region where main effect QTL have been detected in this progeny or in other studies.     

Environment Determines Epistatic Patterns for a ssDNA Virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype–fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype–fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype–phenotype relationship.     

Epistasis and the temporal change in the additive variance-covariance matrix induced by drift

The effect of population bottlenecks on the components of the genetic covariance generated by two neutral independent epistatic loci has been studied theoretically (additive, covA; dominance, covD; additive-by-additive, covAA; additive-by-dominance, covAD; and dominance-by-dominance, covDD). The additive-by-additive model and a more general model covering all possible types of marginal gene action at the single-locus level (additive/dominance epistatic model) were considered. The covariance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium over replicates randomly derived from the base population, after t consecutive bottlenecks of equal size N (derived components). Formulae were obtained in terms of the allele frequencies and effects at each locus, the corresponding epistatic effects and the inbreeding coefficient Ft. These expressions show that the contribution of nonadditive loci to the derived additive covariance (covAt) does not linearly decrease with inbreeding, as in the pure additive case, and may initially increase or even change sign in specific situations. Numerical examples were also analyzed, restricted for simplicity to the case of all covariance components being positive. For additive-by-additive epistasis, the condition covAt > covA only holds for high frequencies of the allele decreasing the metric traits at each locus (negative allele) if epistasis is weak, or for intermediate allele frequencies if it is strong. For the additive/dominance epistatic model, however, covAt > covA applies for low frequencies of the negative alleles at one or both loci and mild epistasis, but this result can be progressively extended to intermediate frequencies as epistasis becomes stronger. Without epistasis the same qualitative results were found, indicating that marginal dominance induced by epistasis can be considered as the primary cause of an increase of the additive covariance after bottlenecks. For all models, the magnitude of the ratio covAt/covA was inversely related to N and t.     

Dynamic genetic interactions determine odor-guided behavior in Drosophila melanogaster

Understanding the genetic architecture of complex traits requires identification of the underlying genes and characterization of gene-by-gene and genotype-by-environment interactions. Behaviors that mediate interactions between organisms and their environment are complex traits expected to be especially sensitive to environmental conditions. Previous studies on the olfactory avoidance response of Drosophila melanogaster showed that the genetic architecture of this model behavior depends on epistatic networks of pleiotropic genes. We performed a screen of 1339 co-isogenic p[GT1]-element insertion lines to identify novel genes that contribute to odor-guided behavior and identified 55 candidate genes with known p[GT1]-element insertion sites. Characterization of the expression profiles of 10 p[GT1]-element insertion lines showed that the effects of the transposon insertions are often dependent on developmental stage and that hypomorphic mutations in developmental genes can elicit profound adult behavioral deficits. We assessed epistasis among these genes by constructing all possible double heterozygotes and measuring avoidance responses under two stimulus conditions. We observed enhancer and suppressor effects among subsets of these P-element-tagged genes, and surprisingly, epistatic interactions shifted with changes in the concentration of the olfactory stimulus. Our results show that the manifestation of epistatic networks dynamically changes with alterations in the environment.     

Less-than-Additive Epistatic Interactions of Quantitative Trait Loci in Tomato

Epistasis plays a role in determining the phenotype, yet quantitative trait loci (QTL) mapping has uncovered little evidence for it. To address this apparent contradiction, we analyzed interactions between individual Lycopersicon pennellii chromosome segments introgressed into an otherwise homogeneous genetic background of L. esculentum (cv. M82). Ten different homozygous introgression lines, each containing from 4 to 58 cM of introgressed DNA, were crossed in a half diallele scheme. The 45 derived double heterozygotes were evaluated in the field for four yield-associated traits, along with the 10 single heterozygotes and M82. Of 180 (45 X 4) tested interactions, 28% were epistatic (P < 0.05) on both linear and geometric scales. The detected epistasis was predominately less-than-additive, i.e., the effect of the double heterozygotes was smaller than the sum of the effects of the corresponding single heterozygotes. Epistasis was also found for homozygous linked QTL affecting fruit mass and total soluble solids. Although the frequency of epistasis was high, additivity was the major component in the interaction of pairs of QTL. We propose that the diminishing additivity of QTL effects is amplified when more loci are involved; this mode of epistasis may be an important factor in phenotype canalization and in breeding.     

Epistatic interaction is an important genetic basis of grain yield and its components in maize

A population of 294 recombinant inbred lines (RIL) derived from Yuyu22, an elite maize hybrid extending broadly in China, has been constructed to investigate the genetic basis of grain yield, and associated yield components in maize. The main-effect quantitative trait loci (QTL), digenic epistatic interactions, and their interactions with the environment for grain yield and its three components were identified by using the mixed linear model approach. Thirty-two main-effect QTL and forty-four pairs of digenic epistatic interactions were detected for the four measured traits in four environments. Our results suggest that both additive effects and epistasis (additive × additive) effects are important genetic bases of grain yield and its components in the RIL population. Only 30.4% of main-effect QTL for ear length were involved in epistatic interactions. This implies that many loci in epistatic interactions may not have significant effects for traits alone but may affect trait expression by epistatic interaction with the other loci.     

Epistatic contributions to quantitative traits in Tribolium castaneum

Summary Triple-testcross experiments were used to analyze epistatic contributions to % hatchability of eggs, age of pupation, number of eggs laid in 24-hour period, and survival from hatching to day 35. Seven diverse inbred lines and the F₁ produced by crossing the two tester lines were examined for the presence of epistasis. There was evidence of epistasis for each of the 4 traits in at least one of the 8 lines tested. Epistasis was a major source of variation in survival in all of the lines tested.     

The effect of epistasis between linked genes on quantitative trait locus analysis

The effect of epistasis between linked genes on quantitative trait locus (QTL) analysis was studied as a function of their contribution to the phenotypic variance and their genetic distance by simulation of F2 (at least 200 individuals) and recombinant inbred line (RIL) populations. Data sets were replicated 100 times. For F2 populations, the presence of epistasis improves the detection of QTLs having effects in opposite directions. Epistasis between linked QTLs (26.5 cM) was poorly detected even when its contribution was relatively high compared to the main effects, and was null for heritabilities lower than 0.10. The detection of false-positive main effects is strongly affected by the distance between epistatic QTLs. The closer they are (≤11.5 cM), the higher the probability of detecting false-positive main-effect QTLs and the lower the probability of detecting epistatic effects. In this case, the presence of main-effect QTLs is due to the deviation of the heterozygote from the homozygotes at each linked interacting QTL and is algebraically explained by the joint effect of the linkage and the additive-by-additive interaction, resulting in a heterosis at a single genomic region in the absence of simulated dominant genetic effects. The number of false-positive main effects only reached nominal levels at about 100 cM. For RIL populations, the number of false positives or the detection of existing epistasis does not depend on the distance, and the power to detect epistatic QTLs is much higher even with small sample sizes and low contributions to the trait. RIL populations are highly recommended to detect epistatic QTLs and to better infer the genetic architecture of a quantitative trait.     

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.     

Dynamic Epistasis under Varying Environmental Perturbations

Epistasis describes the phenomenon that mutations at different loci do not have independent effects with regard to certain phenotypes. Understanding the global epistatic landscape is vital for many genetic and evolutionary theories. Current knowledge for epistatic dynamics under multiple conditions is limited by the technological difficulties in experimentally screening epistatic relations among genes. We explored this issue by applying flux balance analysis to simulate epistatic landscapes under various environmental perturbations. Specifically, we looked at gene-gene epistatic interactions, where the mutations were assumed to occur in different genes. We predicted that epistasis tends to become more positive from glucose-abundant to nutrient-limiting conditions, indicating that selection might be less effective in removing deleterious mutations in the latter. We also observed a stable core of epistatic interactions in all tested conditions, as well as many epistatic interactions unique to each condition. Interestingly, genes in the stable epistatic interaction network are directly linked to most other genes whereas genes with condition-specific epistasis form a scale-free network. Furthermore, genes with stable epistasis tend to have similar evolutionary rates, whereas this co-evolving relationship does not hold for genes with condition-specific epistasis. Our findings provide a novel genome-wide picture about epistatic dynamics under environmental perturbations.     

Influence of Gene Interaction on Complex Trait Variation with Multilocus Models

Although research effort is being expended into determining the importance of epistasis and epistatic variance for complex traits, there is considerable controversy about their importance. Here we undertake an analysis for quantitative traits utilizing a range of multilocus quantitative genetic models and gene frequency distributions, focusing on the potential magnitude of the epistatic variance. All the epistatic terms involving a particular locus appear in its average effect, with the number of two-locus interaction terms increasing in proportion to the square of the number of loci and that of third order as the cube and so on. Hence multilocus epistasis makes substantial contributions to the additive variance and does not, per se, lead to large increases in the nonadditive part of the genotypic variance. Even though this proportion can be high where epistasis is antagonistic to direct effects, it reduces with multiple loci. As the magnitude of the epistatic variance depends critically on the heterozygosity, for models where frequencies are widely dispersed, such as for selectively neutral mutations, contributions of epistatic variance are always small. Epistasis may be important in understanding the genetic architecture, for example, of function or human disease, but that does not imply that loci exhibiting it will contribute much genetic variance. Overall we conclude that theoretical predictions and experimental observations of low amounts of epistatic variance in outbred populations are concordant. It is not a likely source of missing heritability, for example, or major influence on predictions of rates of evolution.     

Dissecting genetic networks underlying complex phenotypes: the theoretical framework

Great progress has been made in genetic dissection of quantitative trait variation during the past two decades, but many studies still reveal only a small fraction of quantitative trait loci (QTLs), and epistasis remains elusive. We integrate contemporary knowledge of signal transduction pathways with principles of quantitative and population genetics to characterize genetic networks underlying complex traits, using a model founded upon one-way functional dependency of downstream genes on upstream regulators (the principle of hierarchy) and mutual functional dependency among related genes (functional genetic units, FGU). Both simulated and real data suggest that complementary epistasis contributes greatly to quantitative trait variation, and obscures the phenotypic effects of many 'downstream' loci in pathways. The mathematical relationships between the main effects and epistatic effects of genes acting at different levels of signaling pathways were established using the quantitative and population genetic parameters. Both loss of function and "co-adapted" gene complexes formed by multiple alleles with differentiated functions (effects) are predicted to be frequent types of allelic diversity at loci that contribute to the genetic variation of complex traits in populations. Downstream FGUs appear to be more vulnerable to loss of function than their upstream regulators, but this vulnerability is apparently compensated by different FGUs of similar functions. Other predictions from the model may account for puzzling results regarding responses to selection, genotype by environment interaction, and the genetic basis of heterosis.     

Epistasis and Its Contribution to Genetic Variance Components

We present a new parameterization of physiological epistasis that allows the measurement of epistasis separate from its effects on the interaction (epistatic) genetic variance component. Epistasis is the deviation of two-locus genotypic values from the sum of the contributing single-locus genotypic values. This parameterization leads to statistical tests for epistasis given estimates of two-locus genotypic values such as can be obtained from quantitative trait locus studies. The contributions of epistasis to the additive, dominance and interaction genetic variances are specified. Epistasis can make substantial contributions to each of these variance components. This parameterization of epistasis allows general consideration of the role of epistasis in evolution by defining its contribution to the additive genetic variance.     

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.