Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.     

Unraveling epistasis with triple testcross progenies of near-isogenic lines

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

Genetics of quantitative traits in Arabidopsis thaliana

The genetic control of 22 quantitative traits, including developmental rates and sizes, was examined in generations of Arabidopsis thaliana derived from the cross between the ecotypes, Columbia (Col) and Landsberg erecta (Ler). The data were obtained from three sets of families raised in the same trial: the 16 basic generations, that is, parents, F(1)'s, F(2)'s, backcrosses, recombinant inbred lines (RILs) and a triple test cross (TTC), the latter produced by crossing the RILs to Col, Ler and their F(1). The data were analysed by two approaches. The first (approach A) involved traditional generation mean and variance component analysis and the second (B), based around the RILs and TTC families, involved marker-based QTL analysis.From (A), genetic differences between Col and Ler were detected for all traits with moderate heritabilities. Height at flowering was the only trait to show heterosis. Dominance was partial to complete for all height traits, and there was no overdominance but there was strong evidence for directional dominance. For most other traits, dominance was ambidirectional and incomplete, with average dominance ratios of around 80%. Epistasis, particularly of the duplicate type that opposes dominance, was a common feature of all traits. The presence of epistasis must imply multiple QTL for all traits.The QTL analysis located 38 significant effects in four regions of chromosomes I, II, IV and V, but not III. QTL affecting rosette size and leaf number were identified in all four regions, with days to maturity on chromosomes IV and V. The only QTL for height was located at the expected position of the erecta gene (chromosome II; 50 cM), but the additive and dominance effects of this single QTL did not adequately explain the generation means. The possible involvement of other interacting height QTL is discussed.     

The role of epistasis in the manifestation of heterosis: a systems-oriented approach

Heterosis is widely used in breeding, but the genetic basis of this biological phenomenon has not been elucidated. We postulate that additive and dominance genetic effects as well as two-locus interactions estimated in classical QTL analyses are not sufficient for quantifying the contributions of QTL to heterosis. A general theoretical framework for determining the contributions of different types of genetic effects to heterosis was developed. Additive x additive epistatic interactions of individual loci with the entire genetic background were identified as a major component of midparent heterosis. On the basis of these findings we defined a new type of heterotic effect denoted as augmented dominance effect di* that comprises the dominance effect at each QTL minus half the sum of additive x additive interactions with all other QTL. We demonstrate that genotypic expectations of QTL effects obtained from analyses with the design III using testcrosses of recombinant inbred lines and composite-interval mapping precisely equal genotypic expectations of midparent heterosis, thus identifying genomic regions relevant for expression of heterosis. The theory for QTL mapping of multiple traits is extended to the simultaneous mapping of newly defined genetic effects to improve the power of QTL detection and distinguish between dominance and overdominance.     

On the genetic control of heterosis for fruit shape in melon (Cucumis melo L.)

The objective of the present work is to study the genetic basis of heterosis for fruit shape (FS) in melon observed in a cross between the Spanish cultivar "Piel de Sapo" (PS) and the Korean accession PI 161375 (Songwang Charmi [SC]) using a set of near-isogenic lines (NILs) with contrasting phenotypes for FS, each carrying a single chromosomal introgression from SC within the genetic background of PS. We investigated the FS of homozygous NILs, hybrids NIL x PS, and all 2-way crosses between NILs to test the main heterosis hypotheses (dominance, overdominance, and epistatic interactions). Gene action of alleles of quantitative trait loci inducing fruit enlargement was dominance, whereas those inducing rounder fruit were additive or recessive. Only minor epistatic interactions were found. Therefore, the most plausible explanation for FS heterosis in this cross is in agreement with the dominance complementation hypothesis. Over 70% of the hybrid heterosis could be achieved by combining just 2 loci, indicating that the genetic control of FS heterosis in this cross is relatively simple. FS is proposed as a reproductive trait in melon because of the high correlation to the number of seeds produced along the fruit longitudinal axis.     

Identification of heterotic metabolite QTL in Arabidopsis thaliana RIL and IL populations

Two mapping populations of a cross between the Arabidopsis thaliana accessions Col-0 and C24 were cultivated and analyzed with respect to the levels of 181 metabolites to elucidate the biological phenomenon of heterosis at the metabolic level. The relative mid-parent heterosis in the F₁ hybrids was <20% for most metabolic traits. The first mapping population consisting of 369 recombinant inbred lines (RILs) and their test cross progeny with both parents allowed us to determine the position and effect of 147 quantitative trait loci (QTL) for metabolite absolute mid-parent heterosis (aMPH). Furthermore, we identified 153 and 83 QTL for augmented additive (Z₁) and dominance effects (Z₂), respectively. We identified putative candidate genes for these QTL using the aracyc database ( http://www.arabidopsis.org/biocyc ), and calculated the average degree of dominance, which was within the dominance and over-dominance range for most metabolites. Analyzing a second population of 41 introgression lines (ILs) and their test crosses with the recurrent parent, we identified 634 significant differences in metabolite levels. Nine per cent of these effects were classified as over-dominant, according to the mode of inheritance. A comparison of both approaches suggested epistasis as a major contributor to metabolite heterosis in Arabidopsis. A linear combination of metabolite levels was shown to significantly correlate with biomass heterosis ( r  = 0.62).     

Additive and over-dominant effects resulting from epistatic loci are the primary genetic basis of heterosis in rice

A set of 148 F9 recombinant inbred lines (RILs) was developed from the cross of an indica cultivar 93-11 and japonica cultivar DTT13,showing strong F1 heterosis.Subsequently,two backcross F1 (BCF1) populations were constructed by backcrossing these 148 RILs to two parents,93-11 and DT713.These three related populations (281BCF1 lines,148 RILs) were phenotyped for six yield-related traits in two locations.Significant inbreeding depression was detected in the population of RILS and a high level of heterosis was observed in the two BCF1 populations.A total of 42 main-effect quantitative trait loci (M-QTLs) and 109 epistatic effect QTL pairs (E-QTLs) were detected in the three related populations using the mixed model approach.By comparing the genetic effects of these QTLs detected in the RILs,BCF1 performance and mid-parental heterosis (HMp),we found that,in both BCF1 populations,the QTLs detected could be classified into two predominant types:additive and over-domlnant loci,which indicated that the additive and over-dominant effect were more important than complete or partially dominance for M-QTLs and E-QTLs.Further,we found that the E-QTLs detected collectively explained a larger portion of the total phenotypic variation than the M-QTLs in both RILs and BCF1 populations.All of these results suggest that additive and over-dominance resulting from epistatic loci might be the primary genetic basis of heterosis in rice.     

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

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

QTL and epistatic analyses of heterosis for seed yield and three yield component traits using molecular markers in rapeseed (Brassica napus L.)

Aiming to explore the basis of heterosis in rapeseed, QTLs for yield and three yield component traits were mapped and the digenic interactions were detected in an F₂ population derived from a cross between two elite rapeseed lines, SI-1300 and Eagle, in this study. Twenty-eight QTLs were detected for the four yield traits, with only two of them detected simultaneously in the Wuhan and Jingmen environments. Additive, partial dominance, dominance, and overdominance effects were all identified for the investigated traits. Dominance (including partial dominance) was shown by 55% of the QTLs, which suggests that dominance is a major genetic basis of heterosis in rapeseed. At the P ?? 0.01 level with 1000 random permutations, 108 and 104 significant digenic interactions were detected in Wuhan and Jingmen, respectively, for the four yield-related traits using all possible locus pairs of molecular markers. Digenic interactions, including additive by additive, additive by dominance, and dominance by dominance, were frequent and widespread in this population. In most cases (78.3%), the interactions occurred among marker loci for which significant effects were not detected by single-locus analysis. Some QTLs (57.1%) detected by single-locus analysis were involved in epistatic interactions. It was concluded that epistasis, along with dominance (including partial dominance), is responsible for the expression of heterosis in rapeseed.     

Something Old and Something New: Wedding Recombinant Inbred Lines with Traditional Line Cross Analysis Increases Power to Describe Gene Interactions

In this paper we present a novel approach to quantifying genetic architecture that combines recombinant inbred lines (RIL) with line cross analysis (LCA). LCA is a method of quantifying directional genetic effects (i.e. summed effects of all loci) that differentiate two parental lines. Directional genetic effects are thought to be critical components of genetic architecture for the long term response to selection and as a cause of inbreeding depression. LCA typically begins with two inbred parental lines that are crossed to produce several generations such as F1, F2, and backcrosses to each parent. When a RIL population (founded from the same P1 and P2 as was used to found the line cross population) is added to the LCA, the sampling variance of several nonadditive genetic effect estimates is greatly reduced. Specifically, estimates of directional dominance, additive x additive, and dominance x dominance epistatic effects are reduced by 92%, 94%, and 56% respectively. The RIL population can be simultaneously used for QTL identification, thus uncovering the effects of specific loci or genomic regions as elements of genetic architecture. LCA and QTL mapping with RIL provide two qualitatively different measures of genetic architecture with the potential to overcome weaknesses of each approach alone. This approach provides cross-validation of the estimates of additive and additive x additive effects, much smaller confidence intervals on dominance, additive x additive and dominance x dominance estimates, qualitatively different measures of genetic architecture, and the potential when used together to balance the weaknesses of LCA or RIL QTL analyses when used alone.     

Gene actions at loci underlying several quantitative traits in two elite rice hybrids

To understand the gene activities controlling nine important agronomic quantitative traits in rice, we applied a North Carolina design 3 (NC III design) analysis to recombinant inbred lines (RILs) in highly heterotic inter- (IJ) and intra-subspecific (II) hybrids by performing the following tasks: (1) investigating the relative contribution of additive, dominant, and epistatic effects for performance traits by generation means analysis and variance component estimates; (2) detecting the number, genomic positions, and genetic effects of QTL for phenotypic traits; and (3) characterizing their mode of gene action. Under an F∞-metric, generation means analysis and variance components estimates revealed that epistatic effects prevailed for the majority of traits in the two hybrids. QTL analysis identified 48 and 66 main-effect QTL (M-QTL) for nine traits in IJ and II hybrids, respectively. In IJ hybrids, 20 QTL (41.7%) showed an additive effect of gene actions, 20 (41.7%) showed partial-to-complete dominance, and 8 (16.7%) showed overdominance. In II hybrids, 34 QTL (51.5%) exhibited additive effects, 14 (21.2%) partial-to-complete dominance, and 18 (27.3%) overdominance. There were 153 digenic interactions (E-QTL) in the IJ hybrid and 252 in the II hybrid. These results suggest that additive effects, dominance, overdominance, and particularly epistasis attribute to the genetic basis of the expression of traits in the two hybrids. Additionally, we determined that the genetic causes of phenotypic traits and their heterosis are different. In the plants we studied, the phenotypic traits investigated and their heterosis were conditioned by different M-QTL and E-QTL, respectively, and were mainly due to non-allelic interactions (epistasis).     

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

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

The genetic basis of heterosis: multiparental quantitative trait loci mapping reveals contrasted levels of apparent overdominance among traits of agronomical interest in maize (Zea mays L.)

Understanding the genetic bases underlying heterosis is a major issue in maize (Zea mays L.). We extended the North Carolina design III (NCIII) by using three populations of recombinant inbred lines derived from three parental lines belonging to different heterotic pools, crossed with each parental line to obtain nine families of hybrids. A total of 1253 hybrids were evaluated for grain moisture, silking date, plant height, and grain yield. Quantitative trait loci (QTL) mapping was carried out on the six families obtained from crosses to parental lines following the "classical" NCIII method and with a multiparental connected model on the global design, adding the three families obtained from crosses to the nonparental line. Results of the QTL detection highlighted that most of the QTL detected for grain yield displayed apparent overdominance effects and limited differences between heterozygous genotypes, whereas for grain moisture predominance of additive effects was observed. For plant height and silking date results were intermediate. Except for grain yield, most of the QTL identified showed significant additive-by-additive epistatic interactions. High correlation observed between heterosis and the heterozygosity of hybrids at markers confirms the complex genetic basis and the role of dominance in heterosis. An important proportion of QTL detected were located close to the centromeres. We hypothesized that the lower recombination in these regions favors the detection of (i) linked QTL in repulsion phase, leading to apparent overdominance for heterotic traits and (ii) linked QTL in coupling phase, reinforcing apparent additive effects of linked QTL for the other traits.     

Quantitative trait loci mapping and the genetic basis of heterosis in maize and rice

Despite its importance to agriculture, the genetic basis of heterosis is still not well understood. The main competing hypotheses include dominance, overdominance, and epistasis. NC design III is an experimental design that has been used for estimating the average degree of dominance of quantitative trait loci (QTL) and also for studying heterosis. In this study, we first develop a multiple-interval mapping (MIM) model for design III that provides a platform to estimate the number, genomic positions, augmented additive and dominance effects, and epistatic interactions of QTL. The model can be used for parents with any generation of selfing. We apply the method to two data sets, one for maize and one for rice. Our results show that heterosis in maize is mainly due to dominant gene action, although overdominance of individual QTL could not completely be ruled out due to the mapping resolution and limitations of NC design III. For rice, the estimated QTL dominant effects could not explain the observed heterosis. There is evidence that additive × additive epistatic effects of QTL could be the main cause for the heterosis in rice. The difference in the genetic basis of heterosis seems to be related to open or self pollination of the two species. The MIM model for NC design III is implemented in Windows QTL Cartographer, a freely distributed software.     

Population-based resequencing analysis of improved wheat germplasm at wheat leaf rust resistance locus Lr21

Little is known about the genetic control of heterosis in the complex polyploid crop species oilseed rape (Brassica napus L.). In this study, two large doubled-haploid (DH) mapping populations and two corresponding sets of backcrossed test hybrids (THs) were analysed in controlled greenhouse experiments and extensive field trials for seedling biomass and yield performance traits, respectively. Genetic maps from the two populations, aligned with the help of common simple sequence repeat markers, were used to localise and compare quantitative trait loci (QTL) related to the expression of heterosis for seedling developmental traits, plant height at flowering, thousand seed mass, seeds per silique, siliques per unit area and seed yield. QTL were mapped using data from the respective DH populations, their corresponding TH populations and from mid-parent heterosis (MPH) data, allowing additive and dominance effects along with digenic epistatic interactions to be estimated. A number of genome regions containing numerous heterosis-related QTL involved in different traits and at different developmental stages were identified at corresponding map positions in the two populations. The co-localisation of per se QTL from the DH population datasets with heterosis-related QTL from the MPH data could indicate regulatory loci that may also contribute to fixed heterosis in the highly duplicated B. napus genome. Given the key role of epistatic interactions in the expression of heterosis in oilseed rape, these QTL hotspots might harbour genes involved in regulation of heterosis (including fixed heterosis) for different traits throughout the plant life cycle, including a significant overall influence on heterosis for seed yield.     

Analysis on additive effects and additive-by-additive epistatic effects of QTLs for yield traits in a recombinant inbred line population of rice

A linkage map consisting of 158 DNA markers were constructed by using a recombinant inbred line (RIL) population derived from the indica-indica rice cross Zhenshan 97B 2 Milyang 46. Quantitative trait loci (QTLs) conditioning grain yield and five yield component traits were determined at the one-locus and two-locus levels, and genotype-by-environment (GE) interactions were analyzed. Thirty-one QTLs were detected to have significant additive effects for yield traits, of which 12 also exhibited significant epistatic effects. Sixteen significant additive-by-additive (AA) interactions were detected, of which nine occurred between QTLs with own additive effects (M_epQTLs), four occurred between QTLs showing epistatic effects only (epQTLs), and three occurred between M_epQTLs and epQTLs. Significant GE interactions were found for six QTLs with additive effects and one AA interaction. Generally, the contributions to the phenotypic variation were higher due to QTL main effects than to epistatic effects. The detection of additive effects and AA effects of a QTL interfered with each other, indicating that the detection of QTLs with main effects, as well as the magnitude and directions of the additive effects, might vary depending on their interactions with other loci.     

A complete solution for dissecting pure main and epistatic effects of QTL in triple testcross design

Epistasis plays an important role in genetics, evolution and crop breeding. To detect the epistasis, triple test cross (TTC) design had been developed several decades ago. Classical procedures for the TTC design use only linear transformations Z(1), Z(2) and Z(3), calculated from the TTC family means of quantitative trait, to infer the nature of the collective additive, dominance and epistatic effects of all the genes. Although several quantitative trait loci (QTL) mapping approaches in the TTC design have been developed, these approaches do not provide a complete solution for dissecting pure main and epistatic effects. In this study, therefore, we developed a two-step approach to estimate all pure main and epistatic effects in the F(2)-based TTC design under the F(2) and F(∞) metric models. In the first step, with Z(1) and Z(2) the augmented main and epistatic effects in the full genetic model that simultaneously considered all putative QTL on the whole genome were estimated using empirical Bayes approach, and with Z(3) three pure epistatic effects were obtained using two-dimensional genome scans. In the second step, the three pure epistatic effects obtained in the first step were integrated with the augmented epistatic and main effects for the further estimation of all other pure effects. A series of Monte Carlo simulation experiments has been carried out to confirm the proposed method. The results from simulation experiments show that: 1) the newly defined genetic parameters could be rightly identified with satisfactory statistical power and precision; 2) the F(2)-based TTC design was superior to the F(2) and F(2:3) designs; 3) with Z(1) and Z(2) the statistical powers for the detection of augmented epistatic effects were substantively affected by the signs of pure epistatic effects; and 4) with Z(3) the estimation of pure epistatic effects required large sample size and family replication number. The extension of the proposed method in this study to other base populations was further discussed.     

家蚕茧质性状的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或染色体区域,育种上可通过适当选配,利用基因的互作效应,同步改良这两个性状。     

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.