The contribution of epistatic pleiotropy to the genetic architecture of covariation among polygenic traits in mice

The contribution that pleiotropic effects of individual loci make to covariation among traits is well understood theoretically and is becoming well documented empirically. However, little is known about the role of epistasis in determining patterns of covariation among traits. To address this problem we combine a quantitative trait locus (QTL) analysis with a two-locus model to assess the contribution of epistasis to the genetic architecture of variation and covariation of organ weights and limb bone lengths in a backcross population of mice created from the M16i and CAST/Ei strains. Significant epistasis was exhibited by 14 pairwise combinations of QTL for organ weights and 10 combinations of QTL for limb bone lengths, which contributed, on average, about 5% of the variation in organ weights and 8% in limb bone lengths beyond that of single-locus QTL effects. Epistatic pleiotropy was much more common in the limb bones (seven of 10 epistatic combinations affecting limb bone lengths were pleiotropic) than the organs (three of the 14 epistatic combinations affecting organ weights were pleiotropic). In both cases, epistatic pleiotropy was less common than single-locus pleiotropy. Epistatic pleiotropy accounted for an average of 6% of covariation among organ weights and 21% of covariation among limb bone lengths, which represented an average of one-fifth (for organ weights) and one-third (for limb bone lengths) of the total genetic covariance between traits. Thus, although epistatic pleiotropy made a smaller contribution than single-locus pleiotropy, it clearly made a significant contribution to the genetic architecture of variation/covariation.     

Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight

Pleiotropy is an aspect of genetic architecture underlying the phenotypic covariance structure. The presence of genetic variation in pleiotropy is necessary for natural selection to shape patterns of covariation between traits. We examined the contribution of differential epistasis to variation in the intertrait relationship and the nature of this variation. Genetic variation in pleiotropy was revealed by mapping quantitative trait loci (QTLs) affecting the allometry of mouse limb and tail length relative to body weight in the mouse-inbred strain LG/J by SM/J intercross. These relationship QTLs (rQTLs) modify relationships between the traits affected by a common pleiotropic locus. We detected 11 rQTLs, mostly affecting allometry of multiple bones. We further identified epistatic interactions responsible for the observed allometric variation. Forty loci that interact epistatically with the detected rQTLs were identified. We demonstrate how these epistatic interactions differentially affect the body size variance and the covariance of traits with body size. We conclude that epistasis, by differentially affecting both the canalization and mean values of the traits of a pleiotropic domain, causes variation in the covariance structure. Variation in pleiotropy maintains evolvability of the genetic architecture, in particular the evolvability of its modular organization.     

The genetic architecture of domestication in the chicken: effects of pleiotropy and linkage

The extent of pleiotropy and epistasis in quantitative traits remains equivocal. In the case of pleiotropy, multiple quantitative trait loci are often taken to be pleiotropic if their confidence intervals overlap, without formal statistical tests being used to ascertain if these overlapping loci are statistically significantly pleiotropic. Additionally, the degree to which the genetic correlations between phenotypic traits are reflected in these pleiotropic quantitative trait loci is often variable, especially in the case of antagonistic pleiotropy. Similarly, the extent of epistasis in various morphological, behavioural and life-history traits is also debated, with a general problem being the sample sizes required to detect such effects. Domestication involves a large number of trade-offs, which are reflected in numerous behavioural, morphological and life-history traits which have evolved as a consequence of adaptation to selective pressures exerted by humans and captivity. The comparison between wild and domestic animals allows the genetic analysis of the traits that differ between these population types, as well as being a general model of evolution. Using a large F(2) intercross between wild and domesticated chickens, in combination with a dense SNP and microsatellite marker map, both pleiotropy and epistasis were analysed. The majority of traits were found to segregate in 11 tight 'blocks' and reflected the trade-offs associated with domestication. These blocks were shown to have a pleiotropic 'core' surrounded by more loosely linked loci. In contrast, epistatic interactions were almost entirely absent, with only six pairs identified over all traits analysed. These results give insights both into the extent of such blocks in evolution and the development of domestication itself.     

Pleiotropic effects on mandibular morphology II: differential epistasis and genetic variation in morphological integration

The evolution of morphological modularity through the sequestration of pleiotropy to sets of functionally and developmentally related traits requires genetic variation in the relationships between traits. Genetic variation in relationships between traits can result from differential epistasis, where epistatic relationships for pairs of loci are different for different traits. This study maps relationship quantitative trait loci (QTLs), specifically QTLs that affect the relationship between individual mandibular traits and mandible length, across the genome in an F2 intercross of the LG/J and SM/J inbred mouse strains (N = 1045). We discovered 23 relationship QTLs scattered throughout the genome. All mandibular traits were involved in one or more relationship QTL. When multiple traits were affected at a relationship QTL, the traits tended to come from a developmentally restricted region of the mandible, either the muscular processes or the alveolus. About one-third of the relationship QTLs correspond to previously located trait QTLs affecting the same traits. These results comprise examples of genetic variation necessary for an evolutionary response to selection on the range of pleiotropic effects.     

An epistatic genetic basis for physical activity traits in mice

We recently identified several (4-8) quantitative trait loci (QTL) for 3 physical activity traits (daily distance, duration, and speed voluntarily run) in an F(2) population of mice derived from an original intercross of 2 strains that exhibited large differences in activity. These QTL cumulatively explained from 11% to 34% of the variation in these traits, but this was considerably less than their total genetic variability estimated from differences among inbred strains. We therefore decided to test whether epistatic interactions might account for additional genetic variation in these traits in this same population of mice. We conducted a full genome epistasis scan for all possible interactions of QTL between each pair of 20 chromosomes. The results of this scan revealed an abundance of epistasis, with QTL throughout the genome being involved in significant interactions. Overall, epistatic effects contributed an average of 26% of the total variation among the 3 activity traits. These results suggest that epistatic interactions of genes may play as important a role in the genetic architecture of physical activity traits as single-locus effects and need to be considered in future candidate gene identification studies.     

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.     

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.     

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.     

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.     

Epistasis plays an important role as genetic basis of heterosis in rice

Ｔｈｅｇｅｎｅｔｉｃｂａｓｉｓｏｆｈｅｔｅｒｏｓｉｓｉｓｓｔｉｌｌａｄｅｂａｔｉｎｇｉｓｓｕｅ.Ｔｗｏｈｙｐｏｔｈｅｓｅｓ,ｔｈｅｄｏｍｉｎａｎｃｅｈｙｐｏｔｈｅｓｉｓａｎｄｔｈｅｏｖｅｒｄｏｍｉｎａｎｃｅｈｙｐｏｔｈｅｓｉｓ,ｂｏｔｈｐｒｏｐｏｓｅｄｉｎ1908[1—3],ｈａｖｅｃｏｍｐｅｔｅｄｆｏｒｍｏｓｔｐａｒｔｏｆｔｈｉｓｃｅｎｔｕｒｙ.Ａｌｔｈｏｕｇｈｍａｎｙｒｅｓｅａｒｃｈｅｒｓｐｒｅｆｅｒｏｎｅｈｙｐｏｔｈｅｓｉｓｔｏｔｈｅｏｔｈｅｒ,ｅｘｐｅｒｉｍｅｎｔａｌｄａｔａａｌｌｏｗｉｎｇｆｏｒｃｒ…     

Simultaneous mapping of epistatic QTL in DU6i × DBA/2 mice

We have mapped epistatic quantitative trait loci (QTL) in an F₂ cross between DU6i × DBA/2 mice. By including these epistatic QTL and their interaction parameters in the genetic model, we were able to increase the genetic variance explained substantially (8.8%–128.3%) for several growth and body composition traits. We used an analysis method based on a simultaneous search for epistatic QTL pairs without assuming that the QTL had any effect individually. We were able to detect several QTL that could not be detected in a search for marginal QTL effects because the epistasis cancelled out the individual effects of the QTL. In total, 23 genomic regions were found to contain QTL affecting one or several of the traits and eight of these QTL did not have significant individual effects. We identified 44 QTL pairs with significant effects on the traits, and, for 28 of the pairs, an epistatic QTL model fit the data significantly better than a model without interactions. The epistatic pairs were classified by the significance of the epistatic parameters in the genetic model, and visual inspection of the two-locus genotype means identified six types of related genotype–phenotype patterns among the pairs. Five of these patterns resembled previously published patterns of QTL interactions.     

Genetic variation underlying pod size and color traits of common bean depends on quantitative trait loci with epistatic effects

Common bean is an important vegetable legume in many regions of the world. Size and color of fresh pods are the key factors for deciding the commercial acceptance of bean as a fresh vegetable. The genetic basis of important horticultural traits of common bean is still poorly understood, which hinders DNA marker-assisted breeding in this crop. Here we report the identification of single-locus and epistatic quantitative trait loci (QTLs), as well as their environment interaction effects for six pod traits, namely width, thickness, length, size index, beak length and color, using an Andean intra-gene pool recombinant inbred line population from a cross between a cultivated common bean and an exotic nuña bean. The QTL analyses performed detected a total of 23 QTLs (single-locus QTLs and epistatic QTLs): five with only individual additive effects and six with only epistatic effects, while the remaining twelve showed both effects. These QTLs were distributed across linkage groups (LGs) 1, 2, 4, 6, 7, 8, 9, 10 and 11; particularly noteworthy are the QTLs for pod size co-located on LGs 1 and 4, indicative of tight linkage or genes with pleiotropic effects governing these traits. Overall, the results obtained showed that additive and epistatic effects are the major genetic basis of pod size and color traits. The mapping of QTLs including epistatic loci for the six pod traits evaluated provides support for implementing marker-assisted selection toward genetic improvement of common bean.     

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.     

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.     

Indirect genetic effects from ecological interactions in Arabidopsis thaliana

Indirect genetic effects arise when genes expressed in one individual affect the expression of traits in other individuals. The importance of indirect genetic effects has been recognized for a diversity of evolutionary processes including kin selection, sexual selection, community structure and multilevel selection, but data regarding their genetic architecture and prevalence throughout the genome remain scarce, especially for interactions between unrelated individuals. Using a set of 411 Bay-0 x Shahdara Arabidopsis recombinant inbred lines grown with Landsberg neighbours, we examined quantitative trait loci (QTL) having direct and indirect effects on size, developmental, and fitness related traits. Using an interval mapping approach, we identified 15 QTL with direct effects and found that 13 of these QTL had significant indirect effects on trait expression in neighbouring plants. These results suggest widespread pleiotropy, as nearly all direct effect QTL have associated pleiotropic indirect effects. Paradoxically, most indirect effects were of the same sign as direct effects, creating a pattern of nearly universal positive pleiotropy that makes most covariances between direct and indirect effects positive. These results are consistent with a complex genetic basis for intraspecific interactions, but suggest that interactions between neighbouring plants are largely positive, rather than negative as would be expected for competition. In addition to their evolutionary and ecological importance, these pleiotropic relationships between DGE and IGE loci have implications for quantitative genetic studies of natural populations as well as experimental design considerations. Additionally, studies that ignore IGEs may over- or underestimate quantitative genetic parameters, as well as the effect of and variance contributed by QTL.     

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.     

Dynamic QTL and epistasis analysis on seedling root traits in upland cotton

Roots are involved in acquisition of water and nutrients, as well as in providing structural support to plant. The root system provides a dynamic model for developmental analysis. Here, we investigated quantitative trait loci (QTL), dynamic conditional QTL and epistatic interactions for seedling root traits using an upland cotton F₂ population and a constructed genetic map. Totally, 37 QTLs for root traits, 35 dynamic conditional QTLs based on the net increased amount of root traits (root tips, forks, length, surface area and volume) (i) after transplanting 10 days compared to 5 days, and (ii) after transplanting 15 days to 10 days were detected. Obvious dynamic characteristic of QTL and dynamic conditional QTL existed at different developmental stages of root because QTL and dynamic conditional QTL had not been detected simultaneously. We further confirmed that additive and dominance effects of QTL qRSA-chr1-1 in interval time 5 to 10 DAT (days after transplant) offset the effects in 10 to 15 DAT. Lots of two-locus interactions for root traits were identified unconditionally or dynamically, and a few epistatic interactions were only detected simultaneously in interval time of 5–10 DAT and 10–15 DAT, suggesting different interactive genetic mechanisms on root development at different stages. Dynamic conditional QTL and epistasis effects provide new attempts to understand the dynamics of roots and provide clues for root architecture selection in upland cotton.     

QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat

In bread wheat, single-locus and two-locus QTL analyses were conducted for seven yield and yield contributing traits using two different mapping populations (P I and P II). Single-locus QTL analyses involved composite interval mapping (CIM) for individual traits and multiple-trait composite interval mapping (MCIM) for correlated yield traits to detect the pleiotropic QTLs. Two-locus analyses were conducted to detect main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL × environment interactions (QE and QQE). Only a solitary QTL for spikelets per spike was common between the above two populations. HomoeoQTLs were also detected, suggesting the presence of triplicate QTLs in bread wheat. Relatively fewer QTLs were detected in P I than in P II. This may be partly due to low density of marker loci on P I framework map (173) than in P II (521) and partly due to more divergent parents used for developing P II. Six QTLs were important which were pleiotropic/coincident involving more than one trait and were also consistent over environments. These QTLs could be utilized efficiently for marker assisted selection (MAS).     

A framework for detecting and characterizing genetic background-dependent imprinting effects

Genomic imprinting, where the effects of alleles depend on their parent-of-origin, can be an important component of the genetic architecture of complex traits. Although there has been a rapidly increasing number of studies of genetic architecture that have examined imprinting effects, none have examined whether imprinting effects depend on genetic background. Such effects are critical for the evolution of genomic imprinting because they allow the imprinting state of a locus to evolve as a function of genetic background. Here we develop a two-locus model of epistasis that includes epistatic interactions involving imprinting effects and apply this model to scan the mouse genome for loci that modulate the imprinting effects of quantitative trait loci (QTL). The inclusion of imprinting leads to nine orthogonal forms of epistasis, five of which do not appear in the usual two-locus decomposition of epistasis. Each form represents a change in the imprinting status of one locus across different classes of genotypes at the other locus. Our genome scan identified two different locus pairs that show complex patterns of epistasis, where the imprinting effect at one locus changes across genetic backgrounds at the other locus. Thus, our model provides a framework for the detection of genetic background-dependent imprinting effects that should provide insights into the background dependence and evolution of genomic imprinting. Our application of the model to a genome scan supports this assertion by identifying pairs of loci that show reciprocal changes in their imprinting status as the background provided by the other locus changes.     

Shared quantitative trait loci underlying the genetic correlation between continuous traits

We review genetic correlations among quantitative traits in light of their underlying quantitative trait loci (QTL). We derive an expectation of genetic correlation from the effects of underlying loci and test whether published genetic correlations can be explained by the QTL underlying the traits. While genetically correlated traits shared more QTL (33%) on average than uncorrelated traits (11%), the actual number of shared QTL shared was small. QTL usually predicted the sign of the correlation with good accuracy, but the quantitative prediction was poor. Approximately 25% of trait pairs in the data set had at least one QTL with antagonistic effects. Yet a significant minority (20%) of such trait pairs have net positive genetic correlations due to such antagonistic QTL 'hidden' within positive genetic correlations. We review the evidence on whether shared QTL represent single pleiotropic loci or closely linked monotropic genes, and argue that strict pleiotropy can be viewed as one end of a continuum of recombination rates where r=0. QTL studies of genetic correlation will likely be insufficient to predict evolutionary trajectories over long time spans in large panmictic populations, but will provide important insights into the trade-offs involved in population and species divergence.