Use of randomization testing to detect multiple epistatic QTLs

Here, we describe a randomization testing strategy for mapping interacting quantitative trait loci (QTLs). In a forward selection strategy, non-interacting QTLs and simultaneously mapped interacting QTL pairs are added to a total genetic model. Simultaneous mapping of epistatic QTLs increases the power of the mapping strategy by allowing detection of interacting QTL pairs where none of the QTL can be detected by their marginal additive and dominance effects. Randomization testing is used to derive empirical significance thresholds for every model selection step in the procedure. A simulation study was used to evaluate the statistical properties of the proposed randomization tests and for which types of epistasis simultaneous mapping of epistatic QTLs adds power. Least squares regression was used for QTL parameter estimation but any other QTL mapping method can be used. A genetic algorithm was used to search for interacting QTL pairs, which makes the proposed strategy feasible for single processor computers. We believe that this method will facilitate the evaluation of the importance at epistatic interaction among QTLs controlling multifactorial traits and disorders.     

Bayesian analyses of multiple epistatic QTL models for body weight and body composition in mice

To comprehensively investigate the genetic architecture of growth and obesity, we performed Bayesian analyses of multiple epistatic quantitative trait locus (QTL) models for body weights at five ages (12 days, 3, 6, 9 and 12 weeks) and body composition traits (weights of two fat pads and five organs) in mice produced from a cross of the F1 between M16i (selected for rapid growth rate) and CAST/Ei (wild-derived strain of small and lean mice) back to M16i. Bayesian model selection revealed a temporally regulated network of multiple QTL for body weight, involving both strong main effects and epistatic effects. No QTL had strong support for both early and late growth, although overlapping combinations of main and epistatic effects were observed at adjacent ages. Most main effects and epistatic interactions had an opposite effect on early and late growth. The contribution of epistasis was more pronounced for body weights at older ages. Body composition traits were also influenced by an interacting network of multiple QTLs. Several main and epistatic effects were shared by the body composition and body weight traits, suggesting that pleiotropy plays an important role in growth and obesity.     

Bayesian analysis for genetic architecture of dynamic traits

The dissection of the genetic architecture of quantitative traits, including the number and locations of quantitative trait loci (QTL) and their main and epistatic effects, has been an important topic in current QTL mapping. We extend the Bayesian model selection framework for mapping multiple epistatic QTL affecting continuous traits to dynamic traits in experimental crosses. The extension inherits the efficiency of Bayesian model selection and the flexibility of the Legendre polynomial model fitting to the change in genetic and environmental effects with time. We illustrate the proposed method by simultaneously detecting the main and epistatic QTLs for the growth of leaf age in a doubled-haploid population of rice. The behavior and performance of the method are also shown by computer simulation experiments. The results show that our method can more quickly identify interacting QTLs for dynamic traits in the models with many numbers of genetic effects, enhancing our understanding of genetic architecture for dynamic traits. Our proposed method can be treated as a general form of mapping QTL for continuous quantitative traits, being easier to extend to multiple traits and to a single trait with repeat records.     

玉米开花期相关性状的QTL分析

利用玉米强优势组合(Mo17×黄早四)自交衍生的191个F_2单株构建了由SSR和AFLP标记组成的分子连锁图谱,用F2进一步自交产生的184个F_(2:3)家系调查散粉期、吐丝期和开花-吐丝间隔期(ASI)的表型值,采用基于混合线性模型的复合区间作图法和相应的作图软件QTLmapper/V2.0,在两个生长环境下定位了与散粉期、吐丝期和ASI相关的QTL数目分别为13、7和5个,检测到3对控制散粉期、17对控制吐丝期和5对控制ASI的上位性效应位点;同时发现了与环境存在显著互作的3个散粉期、3个吐丝期和2个ASI单位点标记区域以及1对散粉期、3对吐丝期和2对ASI上位性效应区域.对玉米散粉期、吐丝期和ASI遗传基础中遗传因素相对作用大小分析表明,加性效应、部分显性效应和上位性效应是玉米开花期相关性状的重要遗传基础.     

水稻株高构成因素的QTL剖析

利用水稻籼粳杂交 (圭 6 30× 0 2 42 8) F1 的花药离体培养建立的一个含 81个 DH家系的作图群体 ,对水稻株高构成因素 (穗长、第 1节间长、……、第 5节间长 )进行基因定位。DH群体中株高构成因素均呈正态分布。相邻的构成因素间呈极显著的正相关 ,而相距较远的构成因素间的相关较弱或不显著。采用 QTL(Quantitative trait lo-cus)分析 ,定位了影响株高构成因素的 6个 QTL:qtl7同时影响穗长和第 1、2、3节间长 ,qtl1 和 qtl2 同时影响第 4和第 5节间长 ,qtl1 0 a和 qtl1 0 b仅影响第 1节间长 ,qtl3 仅影响第 3节间长。采用 QTL 互作分析 ,检测到 19对显著的互作 ,每个构成因素受 2个或 2个以上的 QTL 互作对的影响。并且还发现 ,同一个 QTL 互作对可能影响不同的性状 ,以及一个 QTL 可以分别与不同的 QTL 产生互作而影响同一个性状或影响不同的性状 ,但总的看来 ,加性效应是主要的。这些结果揭示了株高构成因素间相关的遗传基础 ,在水稻育种中运用这些 QTL 将有助于对株高 ,以及对穗长和上部节间长度进行精细的遗传调控。     

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.     

Mapping QTLs with epistatic effects and QTL×environment interactions for plant height using a doubled haploid population in cultivated wheat

Quantitative trait loci (QTLs) for plant height in wheat (Triticum aestivum L.) were studied using a set of 168 doubled haploid (DH) lines, which were derived from the cross Huapei 3/Yumai 57. A genetic linkage map was constructed using 283 SSR and 22 EST-SSR markers. The DH population and the parents were evaluated for wheat plant height in 2005 and 2006 in Tai’an and 2006 in Suzhou. QTL analyses were performed using the software of QTLNetwork version 2.0 based on the mixed linear model. Four additive QTLs and five pairs of epistatic effects were detected, which were distributed on chromosomes 3A, 4B, 4D, 5A, 6A, 7B, and 7D. Among them, three additive QTLs and three pairs of epistatic QTLs showed QTL×environment interactions (QEs). Two major QTLs, Qph4B and Qph4D, which accounted for 14.51% and 20.22% of the phenotypic variation, were located similar to the reported locations of the dwarfing genes Rht1 and Rht2, respectively. The Qph3A-2 with additive effect was not reported in previous linkage mapping studies. The total QTL effects detected for the plant height explained 85.04% of the phenotypic variation, with additive effects 46.07%, epistatic effects 19.89%, and QEs 19.09%. The results showed that both additive effects and epistatic effects were important genetic bases of wheat plant height, which were subjected to environmental modifications, and caused dramatic changes in phenotypic effects. The information obtained in this study will be useful for manipulating the QTLs for wheat plant height by molecular marker-assisted selection (MAS).     

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 x environment interactions in rice. I. heading date and plant height

One hundred twenty six doubled-haploid (DH) rice lines were evaluated in nine diverse Asian environments to reveal the genetic basis of genotype × environment interactions (GEI) for plant height (PH) and heading date (HD). A subset of lines was also evaluated in four water-limited environments, where the environmental basis of G × E could be more precisely defined. Responses to the environments were resolved into individual QTL × environment interactions using replicated phenotyping and the mixed linear-model approach. A total of 37 main-effect QTLs and 29 epistatic QTLs were identified. On average, these QTLs were detectable in 56% of the environments. When detected in multiple environments, the main effects of most QTLs were consistent in direction but varied considerably in magnitude across environments. Some QTLs had opposite effects in different environments, particularly in water-limited environments, indicating that they responded to the environments differently. Inconsistent QTL detection across environments was due primarily to non- or weak-expression of the QTL, and in part to significant QTL × environment interaction effects in the opposite direction to QTL main effects, and to pronounced epistasis. QTL × environment interactions were trait- and gene-specific. The greater GEI for HD than for PH in rice were reflected by more environment-specific QTLs, greater frequency and magnitude of QTL × environment interaction effects, and more pronounced epistasis for HD than for PH. Our results demonstrated that QTL × environment interaction is an important property of many QTLs, even for highly heritable traits such as height and maturity. Information about QTL × environment interaction is essential if marker-assisted selection is to be applied to the manipulation of quantitative traits.Communicated by G. Wenzel     

The main effects, epistatic effects and environmental interactions of QTLs on the cooking and eating quality of rice in a doubled-haploid line population

Amylose content (AC), gel consistency (GC) and gelatinazation temperature (GT) are three important traits that influence the cooking and eating quality of rice. The objective of this study was to characterize the genetic components, including main-effect quantitative trait loci (QTLs), epistatic QTLs and QTL-by-environment interactions (QEs), that are involved in the control of these three traits. A population of doubled haploid (DH) lines derived from a cross between two indica varieties Zhenshan 97 and H94 was used, and data were collected from a field experiment conducted in two different environments. A genetic linkage map consisting of 218 simple sequence repeat (SSR) loci was constructed, and QTL analysis performed using qtlmapper 1.6 resolved the genetic components into main-effect QTLs, epistatic QTLs and QEs. The analysis detected a total of 12 main-effect QTLs for the three traits, with a QTL corresponding to the Wx locus showing a major effect on AC and GC, and a QTL corresponding to the Alk locus having a major effect on GT. Ten digenic interactions involving 19 loci were detected for the three traits, and six main-effect QTLs and two pairs of epistatic QTLs were involved in QEs. While the main-effect QTLs, especially the ones corresponding to known major loci, apparently played predominant roles in the genetic basis of the traits, under certain conditions epistatic effects and QEs also played important roles in controlling the traits. The implications of the findings for rice quality improvement are discussed.     

Unconditional and conditional QTL analysis of kernel weight related traits in wheat (Triticum aestivum L.) in multiple genetic backgrounds

Wheat thousand kernel weight (TKW) is a complex trait, and is largely controlled by several kernel traits, including kernel length (KL) and kernel width (KW). In order to reveal the genetic relationship between TKW and these kernel traits (KW and KL) as accurate as possible, we applied both unconditional and conditional mapping analyses to three distinct genetic populations, one DH population and two RIL populations. This report describes the identifications of 36 unconditional and conditional additive QTLs and 30 pairs of unconditional and conditional epistatic QTLs, all of which are closely associated with TKW. While the conditional additive locus Qtkw1B, detected in the RIL2 population, exhibited the largest contribution, explaining 14.12 % of TKW variance, the unconditional epistatic QTLs Qtkw3A-2/Qtkw5B.1, detected in the DH population, accounted for 11.95 % of phenotypic variance. This study also showed that, compared with unconditional mapping, conditional mapping resulted in very different numbers and different extent of effects of additive and epistatic QTLs that were associated with TKW when TKW was conditioned on kernel traits (KW and KL). These data strongly suggest that KW and KL indeed play a significant role in determining TKW. Furthermore, we demonstrated that the effects of the 25 additive QTLs for TKW were either entirely or largely determined by KW, while the effects of the other 25 additive QTLs for TKW were either entirely or largely affected by KL. We conclude that the conditional mapping can be useful for a better understanding of the interrelationship between the yield contributing traits at the QTL level.     

Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield

Oil content in rapeseed (Brassica napus L.) is generally regarded as a character with high heritability that is negatively correlated with protein content and influenced by plant developmental and yield related traits. To evaluate possible genetic interrelationships between these traits and oil content, QTL for oil content were mapped using data on oil content and on oil content conditioned on the putatively interrelated traits. Phenotypic data were evaluated in a segregating doubled haploid population of 282 lines derived from the F₁ of a cross between the old German cultivar Sollux and the Chinese cultivar Gaoyou. The material was tested at four locations, two each in Germany and in China. QTLMapper version 1.0 was used for mapping unconditional and conditional QTL with additive (a) and locus pairs with additive × additive epistatic (aa) effects. Clear evidence was found for a strong genetic relationship between oil and protein content. Six QTL and nine epistatic locus pairs were found, which had pleiotropic effects on both traits. Nevertheless, two QTL were also identified, which control oil content independent from protein content and which could be used in practical breeding programs to increase oil content without affecting seed protein content. In addition, six additional QTL with small effects were only identified in the conditional mapping. Some evidence was apparent for a genetic interrelationship between oil content and the number of seeds per silique but no evidence was found for a genetic relationship between oil content and flowering time, grain filling period or single seed weight. The results indicate that for closely correlated traits conditional QTL mapping can be used to dissect the genetic interrelationship between two traits at the level of individual QTL. Furthermore, conditional QTL mapping can reveal additional QTL with small effects that are undetectable in unconditional mapping.     

Influence of dent corn genetic backgrounds on QTL detection for plant-height traits and their relationships in high-oil maize

QTL mapping for plant-height traits has not been hitherto reported in high-oil maize. A high-oil maize inbred ‘GY220’ was crossed with two dent maize inbreds (‘8984’ and ‘8622’) to generate two connected F_2:3 populations. Four plant-height traits were evaluated in 284 and 265 F_2:3 families. Single-trait QTL mapping and multiple-trait joint QTL mapping was used to detect QTLs for the traits and the genetic relationship between plant height (PH) and two other plant-height traits. A total of 28 QTLs and 12 pairs of digenic interactions among detected QTLs for four traits were detected in the two F_2:3 families. Only one marker was shared between the two populations. Joint analysis of PH with ear height (EH) and PH with top height (TH) detected 32 additional QTLs. Our results showed that QTL detection for PH was dependent on the genetic background of dent corn inbreds. Multiple-trait joint QTL analysis could increase the number of detected QTLs.     

家蚕茧质性状的QTL定位研究

采用QTLMapper 2．0 QTL作图软件,对F2群体的家蚕全茧量、茧层量、茧层率和蛹体重等性状进行了QTL定位分析,分别检测出7个、6个、2个、8个有显著效应分量的QTLs,分布于7个、5个、2个、7个不同的连锁群。控制全茧量、茧层量的QTLs一般存在复杂的上位性效应。对全茧量性状,有3对QTLs存在显著的加加上位性效应,其中1对还存在加显、显显互作;共有3个QTLs存在显著的显性效应,1个存在显著的加性效应。对茧层量QTLs,发现1对QTLs存在极显著的各项遗传效应,包括上位性效应;1对QTLs被检测到显著的显显互作,1个QTL具有显著的显性效应,并与另一个QTL存在显著的加加互作。茧层率、蛹体重主要受加性或显性的QTLs作用,没有发现茧层率QTLs的上位性效应,蛹体重的有效QTL大都呈现显著的负向显性效应,只有一对QTLs存在显著的加加上位性效应。第2、3、4、11、13、24、34、37、40连锁群是两个或多个性状QTLs分布的共同连锁群。全茧量和茧层量存在共同的QTL或染色体区域,育种上可通过适当选配,利用基因的互作效应,同步改良这两个性状。     

QTL mapping of root traits in a doubled haploid population from a cross between upland and lowland japonica rice in three environments

To genetically dissect drought resistance associated with japonica upland rice, we evaluated a doubled haploid (DH) population from a cross between two japonica cultivars for seven root traits under three different growing conditions (upland, lowland and upland in PVC pipe). The traits included basal root thickness (BRT), total root number (RN), maximum root length (MRL), root fresh weight (RFW), root dry weight (RDW), ratio of root fresh weight to shoot fresh weight (RFW/SFW) and ratio of root dry weight to shoot dry weight (RDW/SDW). The BRT was significantly correlated with the index of drought resistance, which was defined as the ratio of yield under the stress of the upland condition to that under the normal lowland condition. A complete genetic linkage map with 165 molecular markers covering 1,535 cM was constructed. Seven additive quantitative trait loci (QTLs) and 15 pairs of epistatic loci for BRT and RN were identified under upland and lowland conditions, and 12 additive QTLs and 17 pairs of epistatic QTLs for BRT, RN, MRL, RFW, RFW/SFW and RDW/SDW were identified under the PVC pipe condition. Four additive QTLs and one pair of epistatic QTLs controlling IDR were also found. These QTLs individually explained up to 25.6% of the phenotypic variance. QTL × environment (Q × E) interactions were detected for all root traits, and the contributions of these interactions ranged from 1.1% to 19.9%. Five co-localized QTLs controlling RFW and RDW, RFW/SFW, RDW/SDW and IDR, BRT and RN, RN, MRL and IDR were found. Four types of QTLs governing BRT and RN were classified by their detection in the upland and lowland conditions. Some common QTLs for root traits across different backgrounds were also revealed. These co-localized QTLs and common QTLs will facilitate marker-assisted selection for root traits in rice breeding programs.     

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.     

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.     

Characterization of the main effects,epistatic effects and their environmental interactions of QTLs on the genetic basis of yield traits in rice

Main effects, epistatic effects and their environmental interactions of QTLs are all important genetic components of quantitative traits. In this study, we analyzed the main effects, epistatic effects of the QTLs, and QTL by environment interactions (QEs) underlying four yield traits, using a population of 240 recombinant inbred lines from a cross between two rice varieties tested in replicated field trials. A genetic linkage map with 220 DNA marker loci was constructed. A mixed linear model approach was used to detect QTLs with main effects, QTLs involved in digenic interactions and QEs. In total, 29 QTLs of main effects, and 35 digenic interactions involving 58 loci were detected for the four traits. Thirteen QTLs with main effects showed QEs; no QE was detected for the QTLs involved in epistatic interactions. The amount of variations explained by the QTLs of main effect were larger than the QTLs involved in epistatic interactions, which in turn were larger than QEs for all four traits. This study illustrates the ability of the analysis to assess the genetic components underlying the quantitative traits, and demonstrates the relative importance of the various components as the genetic basis of yield traits in this population.     

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.