help button home button Genetics J Neurosci
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Genetics, Vol. 178, 2265-2274, April 2008, Copyright © 2008
doi:10.1534/genetics.107.084871

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Google Scholar
Right arrow Articles by Melchinger, A. E.
Right arrow Articles by Schön, C. C.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Melchinger, A. E.
Right arrow Articles by Schön, C. C.

Genetic Expectations of Quantitative Trait Loci Main and Interaction Effects Obtained With the Triple Testcross Design and Their Relevance for the Analysis of Heterosis

A. E. Melchinger*,1, H. F. Utz* and C. C. Schön{dagger}

* Institute of Plant Breeding, Seed Science, and Population Genetics, University of Hohenheim, 70599 Stuttgart, Germany and {dagger} Plant Breeding, Center of Life and Food Sciences, Technical University of Munich, 85350 Freising, Germany

1 Corresponding author: Institute of Plant Breeding, Seed Science, and Population Genetics, University of Hohenheim, Fruwirthstrasse 21, 70599 Stuttgart, Germany.
E-mail: melchinger{at}uni-hohenheim.de

Manuscript received November 22, 2007. Accepted for publication January 30, 2008.


    ABSTRACT
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
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 F2- and the F{infty}-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.


ESTIMATION of the type of gene action at loci underlying quantitative traits has been a major research focus of quantitative genetics. When defining gene action at single loci we generally distinguish between additive gene action and deviations from additivity due to intralocus allelic interactions, such as dominance and overdominance. With more than one locus, the genotype can also be affected by interlocus interactions, i.e., epistasis. New insights from molecular studies have demonstrated the importance of epistatic interactions in the inheritance of complex traits for a broad spectrum of organisms (SCHADT et al. 2003; BREM et al. 2005). However, the majority of epistatic effects are assumed to be small while available genetic models and experimental designs suffer from limited power of detection. Therefore, new strategies need to be devised for the detection and estimation of epistatic quantitative trait loci (QTL) effects.

Many models have been developed to distinguish between different types of gene action and to estimate the magnitude of genetic effects. ANDERSON and KEMPTHORNE (1954) and GAMBLE (1962) proposed the estimation of additive, dominance, and epistatic effects from first-moment statistics, i.e., generation means. These parameters reflect sums of gene effects over all loci and, consequently, positive and negative effects at individual loci may cancel each other. COCKERHAM (1954) proposed a general model for estimating the type of gene action from second-moment statistics. He developed a set of orthogonal contrasts to partition the genetic variance into additive, dominance, and epistatic components. For populations derived from a cross between two inbred lines, simpler models such as the F2- and the F{infty}-metric model have been described (see VAN DER VEEN 1959). In the presence of epistasis, the choice of metric becomes crucial and depends on the genetic material under study and the experimental design employed for estimating genetic effects. A detailed comparison of the statistical properties and genetic expectations of the F2- and the F{infty}-metric model has been given by KAO and ZENG (2002) and YANG (2004). With an F2 population, the F2 metric is to be preferred because genetic effects under the F2-metric model are orthogonal and thus, in contrast to the F{infty}-metric model, unbiased estimates of genetic effects can be obtained irrespective of the presence of epistasis (ALVAREZ-CASTRO and CALBORG 2007). Furthermore, the genetic variance can be partitioned into eight independent components and genetic covariances are absent, if only first-order interactions are present. MELCHINGER et al. (2007) gave further arguments for the superiority of the F2 metric for the analysis of heterosis with the North Carolina Experiment III (design III), originally devised by COMSTOCK and ROBINSON (1952).

In addition to different genetic models, different experimental designs have been proposed for estimating the type of gene action from second-moment statistics. However, as pointed out by KEARSEY and JINKS (1968), many designs have serious limitations with respect to unbiased estimation of genetic effects. Frequently, absence of epistasis is assumed but a valid test of this assumption is not provided. Furthermore, the different components of variance are estimated with varying precision and power and linkage disequilibrium between loci is not accounted for. The design III partially overcomes these limitations. The experimental units are produced from backcross matings of F2 plants to the two parental lines from which the F2 was derived. Additive and dominance components of variance can be estimated with nearly equal precision under the assumption of diploidy, biallelic and equal gene frequencies, and absence of linkage and epistasis. Thus, an estimate of the average degree of dominance can be obtained from the ratio of the dominance and additive variance components. COCKERHAM and ZENG (1996) extended Comstock and Robinson's analysis of variance (ANOVA) to include linkage and two-locus epistasis for F2 and F3 progenies and developed orthogonal contrasts for marker-aided mapping of QTL using single-marker ANOVA. MELCHINGER et al. (2007) demonstrated the exceptional features of the design III for the identification of QTL contributing to heterosis. They defined a new type of heterotic gene effect, denoted as augmented dominance effect di*, which equals the net contribution of QTL i to midparent heterosis (MPH). It comprises the dominance effect d minus half the sum of additive x additive (aa) epistatic interactions with the genetic background. The novelty of their approach is that QTL significantly contributing to MPH are identified and both dominance and epistasis are accounted for.

An elegant extension of the design III proposed by KEARSEY and JINKS (1968), named the triple testcross (TTC) design, provides a test of significance for the presence of epistasis. In the TTC design, testcrosses are produced not only with the two parental lines but also with the F1 derived from them. For every progeny from a segregating population, e.g., F2 plant or recombinant inbred line (RIL), three sets of data can be created: (i) the average parental testcross performance, (ii) the difference between the parental testcross performances, and (iii) the deviation of testcross progenies with the F1 from the mean of the parental testcrosses. Making use of the advances of marker technology, KEARSEY et al. (2003) and FRASCAROLI et al. (2007) presented experimental results from QTL analyses based on the TTC design with data from Arabidopsis and maize, respectively, but genetic expectations of QTL effects estimated with their respective models were not given and the connection to the analysis of heterosis was not made.

In this study, we give genetic expectations of QTL effects estimated with the TTC design in the presence of epistasis. We show that with the TTC design dominance x additive epistatic interactions of individual QTL with the genetic background can be estimated with one-dimensional genome scans. We also demonstrate that the limitation of the design III in the analysis of heterosis to separate QTL main effects and their epistatic interactions with all other QTL can partially be overcome with the TTC design. Objectives of this study were to (1) extend the theory given by MELCHINGER et al. (2007) to derive quantitative genetic expectations of QTL effects obtained from one-dimensional genome scans with the TTC design using composite-interval mapping (CIM), (2) give quantitative genetic expectations of pairwise interactions between marker loci using two-way ANOVA, (3) derive genetic expectations of variance components of the ANOVA for the TTC progeny tested in a split-plot design for digenic epistasis and arbitrary linkage, and (4) relate our results to the analysis of MPH. Application of our theory to experimental data has been published by KUSTERER et al. (2007).


    THEORY
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Experimental design:
Let us assume a random population of RILs derived from the cross between two homozygous lines P1 and P2. Further, we assume that the RILs are backcrossed to their parental lines and the F1 derived from them, yielding testcross progenies Ht of RIL p (p = 1, ..., n) with testers P1 (t = 1), P2 (t = 2), and F1 (t = 3). The parental line exhibiting superior average testcross performance is denoted as P2.

To obtain maximum precision of progeny means in subsequent QTL analyses we suggest that the testcross progeny Ht are evaluated in a split-plot design with n main plots, each main plot comprising all three testcrosses of the pth RIL. The model for the phenotypic trait values Ytpk can be written as

Formula
with rk being the effect of the kth replication (k = 1, ..., r), gp the genetic effect of the pth RIL, (rg)kp the main plot error term, ht the effect of the tth tester, (hg)tp the interaction between tester t and RIL p, and etpk the subplot error term. While testers are considered fixed, all other effects are assumed random. Following KEARSEY and JINKS (1968), three linear transformations Zs (s = 1, 2, 3) on the performance data are generated with Z1pk = (Y1pk + Y2pk)/2 and Z2pk = Y1pk – Y2pk, and Z3pk = (Y1pk + Y2pk – 2Y3pk). Thus, Zspk denotes the phenotypic value of transformation Zs for RIL p grown in the kth block and Zsp the mean over replicates of Zs for RIL p. Note that Z1 and Z3 are not orthogonal. Thus, partitioning of the genetic variance in the ANOVA is given for the linear transformation Formula instead of Z1pk.

Following MELCHINGER et al. (2007) the genetic constitution of parameters in the expected mean squares of the ANOVA of the TTC design (Table 1) is derived below. Let P1 and P2 differ at loci set Q = {1, ..., q} affecting the quantitative trait of interest and let vi be an indicator variable for the genotype at locus i taking values 0 for homozygous P1 or 2 for homozygous P2. We define the additive effect ai and the dominance effect di at QTL i in accordance with the definition of FALCONER and MACKAY (1996, p. 109), except that ai is negative when the trait-increasing allele is contributed by P1. Epistatic effects between loci i and j are denoted aaij for additive x additive, adij for additive at locus i and dominance at locus j, daij for dominance at locus i and additive at locus j, and ddij for dominance x dominance. The sum of additive x additive epistatic effects over all pairs of QTL is denoted [aa] and aa epistatic interactions of QTL i with the entire genetic background are denoted [aa]. The same notation is followed for ad, da, and dd epistatic effects. To allow extensions to multiple loci we express linkage between loci i and j with the linkage value Formula (SCHNELL 1961), which can be calculated from the recombination frequency rij as Formula. Linkage disequilibrium between loci i and j in the gametic array of progeny derived from cross P1 x P2 is given by the linkage disequilibrium parameter Dij (FALCONER and MACKAY 1996, p. 18). For RILs derived without random mating prior to selfing, Dij can be calculated as Formula. General derivations of expectations and variances for the TTC design based on the theory presented by MELCHINGER et al. (2007) are given in the APPENDIX.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Analysis of variance on linear transformations Formula, Z2, and Z3 calculated from testcross performance of recombinant inbred lines (RIL) evaluated in a split-plot design with n main plots each comprising testcrosses of the pth RIL with testers P1 (t = 1) and P2 (t = 2) and their F1 (t = 3)

 
Assuming digenic epistasis and with Qi denoting the loci set Q excluding element i we obtain

Formula

Formula
and

Formula
As follows from the ANOVA (Table 1), the presence of epistasis can be tested with the linear transformation Z3, the null hypothesis being H0: Formula = 0. In the absence of linkage (Formula = Dij = 0), Formula simplifies to Formula, the sum of additive x additive epistatic effects over all pairs of QTL. Our results on Formula are in agreement with the results of KEARSEY and JINKS (1968), who developed a test of significance for the net contribution of aa effects in the TTC design performed with testcrosses of F2 plants. As mentioned by these authors, if aaij effects have different signs and cancel in the composite effect Formula, then the null hypothesis will be accepted even though strong epistasis may be present.

Assuming digenic epistasis but arbitrary linkage, genetic expectations of the variance components for the TTC are as follows:

Formula

Formula

Formula

We note that with the TTC design partitioning of the genetic variance on the basis of linear transformations Formula, Z2, and Z3 differs from the design III. The progeny variance among RILs in the TTC design is reduced compared with the design III (Formula < Formula). Both variance components comprise the same genetic effects, but the influence of da and dd epistatic effects is decreased by the factor Formula, i.e., by two-thirds for unlinked loci. The variance component from interaction of RILs with the two parental testers (Formula) is identical for the TTC design and the design III. The genetic variance component arising from Z3 is a complex function of different epistatic effects. In the absence of linkage, rejection of the null hypothesis H0: Formula = 0 provides evidence for epistasis of type ad and/or da.

As can be seen from Table 1, the proposed split-plot design in which each of n main plots comprises three subplots, i.e., the three testcrosses of the pth RIL, is advantageous compared with a randomized complete block design (RCB). The standard error of progeny means Z2 and Z3 is calculated from the subplot error, which is expected to be smaller than the error variance of the RCB. Competition effects between the three different testcross progenies are not expected due to equal inbreeding coefficients (F = 0.5). Compared with the design employed by FRASCAROLI et al. (2007), who assigned testers to t main plots and progenies from the same tester to subplots, this design has the advantage that the test for epistasis summed over all loci (Formula) becomes more powerful due to more degrees of freedom of the main plot error ((n 1)(r – 1) > (t – 1)(r – 1)). In addition, the subplot error variance should be decreased due to the small number of subplots per main plot resulting in higher precision of progeny means Z2 and Z3.

QTL analysis with the TTC design:
MELCHINGER et al. (2007) derived quantitative genetic expectations of QTL effects obtained with the design III and RILs. Using CIM and assuming digenic epistasis they demonstrated that in one-dimensional genome scans on Z1 and Z2 the contrast of the two (unobservable) homozygous genotype classes at QTL i equals the augmented additive (ai*) and dominance (di*) effects:

Formula

Formula
They concluded that the sum of QTL effects equals genotypic expectations for the parental difference (PD) and MPH. Thus,

Formula
and

Formula

For the identification of QTL affecting the PD or MPH it is favorable that in ai* and di* main effects and epistatic interactions of QTL i with the genetic background are confounded. However, the dissection of augmented QTL effects into their components is desirable when the relative contribution of the individual effects is of interest, i.e., the additive effect ai at QTL i and [da] epistasis contributing to the PD and the dominance effect di at QTL i and [aa] epistasis contributing to MPH.

The contribution of the additive effect ai at QTL i and its [da] epistatic interactions can be estimated from one-dimensional genome scans with H3 and the linear transformation Z3. Following MELCHINGER et al. (2007), with CIM, i.e., estimating the QTL position and including cofactors in the model, we obtain the following quantitative genetic expectations of QTL effects at QTL i:

Formula

Formula
Thus, with the TTC design, genome scans on H3 and Z3 can be adopted for estimating to what extent the individual effects ai and [da] contribute to the augmented additive effect ai* in the absence of linkage.

On the basis of one-dimensional genome scans, the TTC design does not provide a solution to the dissection of the augmented dominance effect di* into its components, i.e., the dominance effect of QTL i and [aa]. However, as commonly practiced in QTL analyses, digenic epistasis can be estimated by two-way ANOVA on the basis of interactions of two-locus combinations of marker genotypes. In the following, we derive quantitative genetic expectations of contrasts for two-locus marker genotypes with linear transformations Zs of the TTC design. We assume two QTL i and j and two marker loci m1 and m2, each with genotype classes v (vi = 0, 2; vj = 0, 2; Formula; Formula), and define the vector Formula with Formula referring to the conditional probability of the QTL genotype vivj (vivj = 22, 20, 02, 00) given marker genotype Formula (Formula = 22, 20, 02, 00). Following the parameterization of gamete frequencies given by SCHNELL (1961), the four-locus genotype frequencies of RILs Formula can be expressed by using six two-locus Formula and one four-locus linkage disequilibrium parameter Formula as

Formula

Following MELCHINGER et al. (2007), conditional expectations of linear functions Zs (Zs = Z1, Z1*, Z2, Z3) and testcross progenies Ht (t = 1, 2, 3) are obtained by

Formula
and

Formula
with the matrices Formula and Formula denoting the coefficients of genetic effects given the genotype of the parental RIL, the tester t, and for Formula the linear transformation s (for details see the APPENDIX). Formula denotes the vector of genetic effects defined according to the F2 metric. From this, we obtain the expectations of the interaction between markers m1 and m2 for linear functions Zs and testcross progenies Ht as

Formula
and

Formula

For arbitrary linkage between all four loci (i, j, m1, m2) calculations of conditional QTL genotype frequencies become rather unwieldy. Therefore, we exemplify our derivations for the special case in which marker loci m1 and m2 are unlinked (Formula) and QTL i is linked to marker m1 and QTL j to marker m2 (Formula). Conditional probabilities of QTL genotypes at loci i and j are given in Table 2. Summation over QTL i and j yields genotypic expectations for interactions between marker pairs with linear transformations Zs calculated from the TTC design,

Formula

Formula

Formula

Formula
where Q(m1) denotes all loci in set Q in linkage disequilibrium with marker m1.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Frequencies (Formula) of the four possible QTL genotypes of recombinant inbred line (RIL) parents at QTL i and j conditional on the marker genotype v (v = 0, 2) at marker loci m1 and m2 calculated using the two-locus linkage disequilibrium parameters Formula and Formula

 
For the dissection of heterotic-effect di* into its components, estimates of aa epistatic interactions of QTL i with other QTL in the genome are of particular interest and can be obtained from two-way ANOVAs of marker interactions with H3. The genotypic expectation for interactions between marker pairs with testcross progenies H3 is given by

Formula

With cofactors in the model and assuming (i) that QTL linked to markers m1 and m2 interact only with each other and not with other QTL and (ii) Formula, then the genotypic expectations of interactions between markers simplify to

Formula

Formula

Formula

Formula

Formula

While augmented additive and dominance effects are estimated in one-dimensional scans of Z1 and Z2, respectively, estimates of two-way marker interactions on Z1 and Z1* yield a confounded estimate of aa and dd interactions, and Z2 captures ad and da epistasis. A first estimate of genetic background interactions contributing to di* at QTL i can be obtained with two-way marker interactions on H3. Unbiased estimates of dd epistasis can be obtained with two-way ANOVAs on Z3.

It becomes obvious that genetic expectations of two-way marker interactions obtained with the TTC design are not straightforward extensions of effects estimated with the same linear transformation in one-dimensional scans. Note that this is also true under the F{infty} model employed by KEARSEY et al. (2003) and FRASCAROLI et al. (2007) and for the separate analysis of backcross progenies as performed by STUBER et al. (1992). A summary of genetic expectations of QTL effects estimated from one-dimensional genome scans and two-way ANOVAs of marker interactions on (i) the three possible testcrosses (H1, H2, H3), (ii) the RIL lines (H4), and (iii) the linear transformations Zs are given in Table 3 for both the F2- and the F{infty}-metric models. Genetic expectations of QTL effects comprising epistatic effects but no main effects are identical for the two models. Genetic expectations of QTL effects obtained with one-dimensional genome scans on Z1 and Z2 differ, because under the F{infty} metric estimates of QTL main effects are confounded with epistasis. Following YANG (2004), genetic expectations for the F{infty}-metric model are obtained by substituting ai and di of the F2 model with Formula and Formula, respectively. As evident from Table 3, additional complexity is introduced by the use of the F{infty} metric. In the general case of populations with arbitrary gene frequencies the NOIA model devised by ALVAREZ-CASTRO and CALBORG (2007) could be used for transforming the QTL effects determined in such a population to the genetic effects defined under the F2 or the F{infty} metric.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Genetic expectations of QTL effects obtained with the triple testcross (TTC) design and estimated from one-dimensional genome scans Formula and two-way ANOVAs of marker interactions Formula between unlinked marker loci m1 and m2 (Formula) on the three possible testcrosses (H1, H2, H3), the RIL lines (H4), and the linear transformations Zs under the F2- and F{infty}-metric models

 

    DISCUSSION
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Genetic expectations of QTL effects:
In this study, genetic expectations of QTL effects are estimated with the TTC design. Accounting for all types of digenic epistasis and arbitrary linkage, genetic expectations are given for one-dimensional genome scans and two-way marker ANOVAs under both the F2- and the F{infty}-metric models. These theoretical results contribute significantly to the interpretation of QTL mapping experiments estimating the type of gene action with the TTC design and the design III. The advantages of these two designs have been widely recognized but, to date, genetic expectations of QTL effects have been given for the design III using only single-marker ANOVA (COCKERHAM and ZENG 1996) and one-dimensional genome scans with CIM (MELCHINGER et al. 2007).

When making inferences on the predominant type of gene action, profound knowledge on the genetic expectations of QTL effects is crucial, as was demonstrated by COCKERHAM and ZENG (1996) with a reanalysis of data from STUBER et al. (1992). STUBER et al. (1992) estimated the type of gene action in a marker-aided design III experiment with F3 lines of maize. They performed the QTL analysis separately for the testcrosses with each parent (i.e., on H1 and H2) and detected that overdominance but not epistasis plays a major role in the inheritance of grain yield. Accounting for epistasis and performing a joint analysis of both testcrosses, COCKERHAM and ZENG (1996) found mostly QTL with dominant and epistatic gene action. Similarly divergent results were found for experimental studies on rice grain yield. XIAO et al. (1995) found dominance to be the most important type of gene action. When accounting for epistasis in the genetic model, LI et al. (2001) and LUO et al. (2001) detected strong evidence for overdominance and epistasis. Thus, interpretation of experimental results from QTL mapping studies on the predominant type of gene action can be severely affected by the choice of statistical model and the provision of epistasis. As evident from Table 3, genetic expectations of QTL effects estimated with the TTC design and the design III are complex and comprise both main and epistatic effects. This fact has been widely neglected in the literature. The theoretical results presented here assist in the interpretation of experimental results obtained with two of the major designs employed in the analysis of gene action and heterosis. We therefore believe that reanalysis of previously collected data sets with the statistical methods presented here as well as the joint analysis of similar data sets with a special focus on epistasis will be rewarding. The general quantitative genetic theory given in the APPENDIX allows derivation of genetic expectations of QTL effects for experiments where production of testcross progenies was performed with double-haploid lines or F2 or F3 populations.

Analysis of heterosis with the TTC design:
MELCHINGER et al. (2007) demonstrated that in the analysis of the genetic causes of heterosis we need to identify genomic regions that harbor augmented dominance effects di* instead of identifying QTL with maximum dominance di and dd interactions that control F1 performance. With the design III, the confounding of QTL effects with epistatic background variation is desirable for the identification of heterotic QTL. However, it is also a limitation because partitioning of augmented QTL effects ai* and di* into their main and epistatic components is not possible. By adding the F1 as a third tester, this limitation can partially be overcome. The contribution of [da] epistatic interactions of QTL i with the genetic background to the augmented additive effect ai* can be estimated with genome scans on Z3 and two-way marker interactions with H3 can be used to estimate aa interactions of individual QTL with the genetic background. Because the search for interactions may be restricted to those QTL with significant augmented dominance effects (di*), the problem of multiple testing is alleviated. However, two-way marker interactions will provide only a rough estimate of the contributions of [aa] epistasis to di* due to limited power of detection. Furthermore, QTL with a significant positive dominance effect di and positive [aa] epistasis may remain undetected in a genomewide scan with Z2 if the two effects cancel each other, resulting in nonsignificant di* effects. We are currently in the process of developing new experimental designs that allow separate estimation of di and [aa] components contributing to di* and, consequently, MPH in one-dimensional genome scans.

Detection of epistatic interactions:
The development of statistical tools and powerful experimental designs for an efficient identification of genetic interactions is a major challenge in the analysis of quantitative traits. As a result of limited statistical power, detection of significant epistatic QTL interactions has proved difficult in marker-aided studies on complex traits, such as yield, even with dense marker coverage and large populations (e.g., SCHÖN et al. 2004; MIHALJEVIC et al. 2005). On the contrary, the presence of significant epistasis has been demonstrated when clearly defined genes were investigated and efficient molecular tools were at hand. DOEBLEY et al. (1995) demonstrated dependency of QTL effects on genetic background for plant and inflorescence architecture in maize and teosinte. Epistatic interactions of QTL and expression QTL (eQTL) involved in regulation of flowering in Arabidopsis were reported by KEURENTJES et al. (2007). KROYMANN and MITCHELL-OLDS (2005) cloned two QTL for growth rate in Arabidopsis exhibiting significant epistasis with the genetic background. The authors pointed out that the two QTL would not have been detected with classical QTL analysis approaches and that we are likely to introduce an ascertainment bias because QTL with significant epistatic interaction effects might not be representative of the majority of QTL with small effects contributing to gene networks.

In this study, we have developed a one-dimensional genome scan for epistatic interactions of type dominance x additive. QTL detected with CIM on the linear transformation Z3 exhibit significant [da] epistasis with the genetic background. With this method, statistical power of detection is increased compared with statistical tests for epistasis based on interactions of all possible marker pairs, because the number of significance tests is greatly reduced and thus safeguarding against a high false discovery rate becomes less rigorous. Employing diallel crosses of three homozygous parents, JANNINK and JANSEN (2001) proposed a one-dimensional search for significant background interactions of QTL. With simulated data, they reported a twofold increase in power with the proposed one-dimensional search compared with standard two-dimensional searches. BLANC et al. (2006) employed a similar method in an experimental study with multiparental crosses of maize and found substantial evidence for QTL x genetic-background interactions, especially for grain yield. With both methods the partitioning of epistasis into its components is not feasible. Nevertheless, with the increasing availability of QTL mapping populations derived from multiline crosses these methods are a valuable contribution to the identification of QTL with significant interactions with the genetic background.

In addition to genetic interactions of type dominance x additive [da], interactions of type additive x dominance [ad] can be detected in the absence of linkage with a one-dimensional genome scan on the linear transformation Z4 = 2H3H4 if, in addition to the testcrosses of the RILs, their line per se performance (H4) is tested, as was done by FRASCAROLI et al. (2007) and KUSTERER et al. (2007) (see the APPENDIX and Table 3). Epistatic interactions of type [ad] or [da] of QTL i may not play a major role in elite breeding material, because when summed over the entire genetic background, they are likely to cancel each other due to opposite signs of individual interactions. However, so far we have not been able to verify this hypothesis in QTL analyses. Many experimental studies demonstrated only minor importance of significant ad or da epistasis on the basis of two-way ANOVAs of marker interactions (e.g., HUA et al. 2003) but it has been difficult to distinguish between true and false negatives due to the limitations of statistical tests. With the variance component Formula, the TTC design provides a significance test for the presence of ad and da epistasis. With the one-dimensional genome scan on Z3 and Z4 (if data on per se performance of RILs are available) and reasonable sample sizes we can achieve sufficient power to estimate the magnitude of [ad] or [da] interactions of individual QTL as well as the hypothesis [ad] = [da]. Once QTL are identified to interact with the genetic background, two-way marker ANOVAs on Z2 can be performed to obtain individual estimates of QTL interactions between these QTL and other QTL in the genome Formula.

In conclusion, we are still at the beginning of understanding the complex interactions of individual genes and gene networks even with extensive genomic tools at hand. Knowledge about genetic expectations of QTL effects in the presence of epistasis will facilitate the assessment of gene action and function and will help elucidate the quantitative genetic basis of heterosis. As pointed out by JANNINK and JANSEN (2001), marker-assisted transfer of single genes affecting quantitative traits may be a fruitless endeavor if alleles show strong epistasis and fail to interact with the target genome in the same way as with the donor genome. On the other hand, with a more profound understanding of gene interactions, breeders may be empowered to utilize new alleles from nonadapted genetic resources or genetic engineering that exhibit favorable epistasis with the genetic background. In combination with newly developed statistical methods, such as Bayesian approaches (XU and JIA 2007), multiple-interval mapping (KAO et al. 1999), or two-stage analyses (BREM et al. 2005), powerful experimental designs can significantly increase the efficiency of experiments analyzing phenotypic data on agronomic traits such as yield. Furthermore, the same experimental designs can be used for molecular studies on the quantitative genetics of transcription, protein, or metabolite data. We believe that the analysis of gene interactions will be of increasing importance in future molecular and quantitative genetics research and that the theoretical results from this study provide improved analytical tools for the interpretation of a wide range of experimental data.


View this table:
[in this window]
[in a new window]

 
TABLE A1 Coefficients in 4 x 9 matrix Ht (t = 1, 2, 3, 4) for the conditional expectation of the genotypic value of testcross progeny Ht of a RIL with genotype vivj (vivj = 22, 20, 02, 00) at QTL i and j for testcross performance (t = 1, 2, 3) with tester P1 (t = 1), P2 (t = 2), or F1 (t = 3) or per se performance (t = 4)

 

    APPENDIX: GENERAL DERIVATION OF EXPECTATIONS, VARIANCES, AND COVARIANCES OF LINEAR TRANSFORMATIONS ZS FOR THE TTC DESIGN
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Following MELCHINGER et al.'s (2007) expectations, variances and covariances of linear transformations Zs (Zs = Z1, Z1*, Z2, Z3, Z4) are given by

Formula

Formula
and

Formula

Assuming digenic epistasis Formula denotes the frequencies of the four possible genotypes at two QTL i and j, F denotes a diagonal matrix with these frequencies on the diagonal, and E denotes the vector of genetic effects FormulaAD; i.e., Formula.

Elements of the matrix Ht denote the coefficients of genetic effects FormulaAD in the conditional genotypic expectation of testcross progeny Ht of a RIL with genotype vivj (vivj = 22, 20, 02, 00) at QTL i and j for testcross performance (t = 1, 2, 3) with tester P1 (t = 1), P2 (t = 2), or F1 (t = 3) or per se performance (t = 4). Assuming digenic epistasis, elements of Ht are given in Table A1. The matrices Ks are obtained for the TTC design by calculating

Formula


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
This project was supported by the Deutsche Forschungsgemeinschaft (German Research Foundation) under the priority research program "Heterosis in Plants" (research grants ME931/4-1 and ME931/4-2).


    LITERATURE CITED
 TOP
 ABSTRACT
 THEORY
 DISCUSSION
 APPENDIX: GENERAL DERIVATION OF...
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 

ALVAREZ-CASTRO, J. M., and Ö. CALBORG, 2007 A unified model for functional and statistical epistasis and its application in quantitative trait loci analysis. Genetics 176: 1151–1167.[Abstract/Free Full Text]

ANDERSON, V. L., and O. KEMPTHORNE, 1954 A model for the study of quantitative inheritance. Genetics 39: 883–898.[Free Full Text]

BLANC, G., A. CHARCOSSET, B. MANQIN, A. GALLAIS and L. MOREAU, 2006 Connected populations for detecting quantitative trait loci and testing for epistasis: an application in maize. Theor. Appl. Genet. 113: 206–224.[CrossRef][Medline]

BREM, R. B., J. D. STOREY, J. WHITTLE and L. KRUGLYAK, 2005 Genetic interactions between polymorhpisms that affect gene expression in yeast. Nature 436: 701–703.[CrossRef][Medline]

COCKERHAM, C. C., 1954 An extension of the concept of partitioning hereditary variance for analysis of covariances among relatives when epistasis is present. Genetics 39: 859–882.[Free Full Text]

COCKERHAM, C. C., and Z. B. ZENG, 1996 Design III with marker loci. Genetics 143: 1437–1456.[Abstract]

COMSTOCK, R. E., and H. F. ROBINSON, 1952 Estimation of average dominance of genes, pp. 494–516 in Heterosis, edited by J. W. GOWEN. Iowa State College Press, Ames, IA.

DOEBLEY, J., A. STEC and C. GUSTUS, 1995 Teosinte branched1 and the origin of maize: evidence for epistasis and the evolution of dominance. Genetics 141: 333–346.[Abstract]

FALCONER, D. S., and T. F. C. MACKAY, 1996 Introduction to Quantitative Genetics. Longman, Harlow, UK.

FRASCAROLI, E., M. A. CANÉ, P. LANDI, G. PEA, L. GIANFRANCESCHI et al., 2007 Classical genetic and quantitative trait loci analyses of heterosis in a maize hybrid between two elite inbred lines. Genetics 176: 625–644.[Abstract/Free Full Text]

GAMBLE, E. E., 1962 Gene effects in corn (Zea mays L.). I. Separation and relative importance of gene effects for yield. Can. J. Plant Sci. 42: 339–348.

HUA, J., Y. XING, W. WU, C. XU, X. SUN et al., 2003 Single-locus heterotic effects and dominance by dominance interactions can adequately explain the genetic basis of heterosis in an elite rice hybrid. Proc. Natl. Acad. Sci. USA 100: 2574–2579.[Abstract/Free Full Text]

JANNINK, J.-L., and R. JANSEN, 2001 Mapping epistatic quantitative trait loci with one-dimensional genome searches. Genetics 157: 445–454.[Abstract/Free Full Text]

KAO, C. H., and Z. B. ZENG, 2002 Modeling epistasis of quantitative trait loci using Cockerham's model. Genetics 160: 1243–1261.[Abstract/Free Full Text]

KAO, C. H., Z. B. ZENG and R. D. TEASDALE, 1999 Multiple interval mapping for quantitative trait loci. Genetics 152: 1203–1216.[Abstract/Free Full Text]

KEARSEY, M. J., and J. L. JINKS, 1968 A general method of detecting additive, dominance and epistatic variation for metrical traits. I. Theory. Heredity 23: 403–409.

KEARSEY, M. J., H. S. POONI and N. H. SYED, 2003 Genetics of quantitative traits in Arabidopsis thaliana. Heredity 91: 456–464.[CrossRef][Medline]

KEURENTJES, J. J. B., J. FU, I. R. TERPSTRA, J. M. GARCIA, G. VAN DEN ACKERVEKEN et al., 2007 Regulatory network construction in Arabidopsis by using genome-wide gene expression quantitative trait loci. Proc. Natl. Acad. Sci. USA 104: 1708–1713.[Abstract/Free Full Text]

KROYMANN, J., and T. MITCHELL-OLDS, 2005 Epistasis and balanced polymorphism influencing complex trait variation. Nature 435: 95–98.[CrossRef][Medline]

KUSTERER, B., H. P. PIEPHO, H. F. UTZ, C. C. SCHÖN, J. MUMINOVIC et al., 2007 Heterosis for biomass-related traits in Arabidopsis investigated by a novel QTL analysis of the triple testcross design with recombinant inbred lines. Genetics 177: 1839–1850.[Abstract/Free Full Text]

LI, Z. K, L. J. LUO, H. W. MEI, D. L. WANG, Q. Y. SHU et al., 2001 Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield. Genetics 158: 1737–1753.[Abstract/Free Full Text]

LUO, L. J., Z. K. LI, H. W. MEI, Q. Y. SHU, R. TABIEN et al., 2001 Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. II. Grain yield components. Genetics 158: 1755–1771.[Abstract/Free Full Text]

MELCHINGER, A. E., H. F. UTZ, H. P. PIEPHO, Z. B. ZENG and C. C. SCHÖN, 2007 The role of epistasis in the manifestation of heterosis—a systems-oriented approach. Genetics 177: 1815–1825.[Abstract/Free Full Text]

MIHALJEVIC, R., H. F. UTZ and A. E. MELCHINGER, 2005 No evidence for epistasis in hybrid and per se performance of elite European flint maize inbreds from generation means and QTL analyses. Crop Sci. 45: 2605–2613.[Abstract/Free Full Text]

SCHADT, E. E., S. A. MONKS, T. A. DRAKE, A. J. LUSIS, N. CHE et al., 2003 Genetics of gene expression surveyed in maize, mouse and man. Nature 422: 297–302.[CrossRef][Medline]

SCHNELL, F. W., 1961 Some general formulations of linkage effects in inbreeding. Genetics 46: 947–957.[Free Full Text]

SCHÖN, C. C., F. H. UTZ, S. GROH, B. TRUBERG, S. OPENSHAW et al., 2004 Quantitative trait locus mapping based on resampling in a vast maize testcross experiment and its relevance to quantitative genetics for complex traits. Genetics 167: 485–498.[Abstract/Free Full Text]

STUBER, C. W., S. E. LINCOLN, D. W. WOLFF, T. HELENTJARIS and E. S. LANDER, 1992 Identification of genetic factors contributing to heterosis in a hybrid from two elite maize inbred lines using molecular markers. Genetics 132: 823–839.[Abstract]

VAN DER VEEN, J. H., 1959 Tests of non-allelic interaction and linkage for quantitative characters in generations derived from two diploid pure lines. Genetica 30: 201–232.[CrossRef][Medline]

XIAO, J., J. LI, L. YUAN and S. D. TANKSLEY, 1995 Dominance is the major genetic basis of heterosis in rice as revealed by QTL analysis using molecular markers. Genetics 140: 745–754.[Abstract]

XU, S., and Z. JIA, 2007 Genomewide analysis of epistatic effects for quantitative traits in barley. Genetics 175: 1955–1963.[Abstract/Free Full Text]

YANG, R. C., 2004 Epistasis of quantitative trait loci under different gene action models. Genetics 167: 1493–1505.[Abstract/Free Full Text]

Communicating editor: R. W. DOERGE





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Google Scholar
Right arrow Articles by Melchinger, A. E.
Right arrow Articles by Schön, C. C.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Melchinger, A. E.
Right arrow Articles by Schön, C. C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2008 by the Genetics Society of America.