Email updates

Keep up to date with the latest news and content from BMC Plant Biology and BioMed Central.

Open Access Research article

QTL architecture of reproductive fitness characters in Brassica rapa

Jennifer M Dechaine1*, Marcus T Brock2 and Cynthia Weinig2

Author Affiliations

1 Department of Biological Sciences, Central Washington University, Ellensburg, WA 98926, USA

2 Department of Botany, University of Wyoming, Laramie, WY 82071, USA

For all author emails, please log on.

BMC Plant Biology 2014, 14:66  doi:10.1186/1471-2229-14-66


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/14/66


Received:16 July 2013
Accepted:3 March 2014
Published:18 March 2014

© 2014 Dechaine et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

Background

Reproductive output is critical to both agronomists seeking to increase seed yield and to evolutionary biologists interested in understanding natural selection. We examine the genetic architecture of diverse reproductive fitness traits in recombinant inbred lines (RILs) developed from a crop (seed oil) × wild-like (rapid cycling) genotype of Brassica rapa in field and greenhouse environments.

Results

Several fitness traits showed strong correlations and QTL-colocalization across environments (days to bolting, fruit length and seed color). Total fruit number was uncorrelated across environments and most QTL affecting this trait were correspondingly environment-specific. Most fitness components were positively correlated, consistent with life-history theory that genotypic variation in resource acquisition masks tradeoffs. Finally, we detected evidence of transgenerational pleiotropy, that is, maternal days to bolting was negatively correlated with days to offspring germination. A QTL for this transgenerational correlation was mapped to a genomic region harboring one copy of FLOWERING LOCUS C, a genetic locus known to affect both days to flowering as well as germination phenotypes.

Conclusions

This study characterizes the genetic structure of important fitness/yield traits within and between generations in B. rapa. Several identified QTL are suitable candidates for fine-mapping for the improvement of yield in crop Brassicas. Specifically, brFLC1, warrants further investigation as a potential regulator of phenology between generations.

Keywords:
Fitness components; Life-history traits; Phenotypic plasticity; Transgenerational effects; Yield; Brassica rapa

Background

Female reproductive fitness is an important measure of overall crop yield for many cultivated species, such as grain crops and oilseeds. In an evolutionary context, reproductive fitness is commonly used to estimate selection, as in studies investigating crop allele persistence in wild populations [1-3]. Fitness and yield are determined by complex interactions among traits expressed throughout plant development, such as flowering time [4,5], fruit shape [6], seed size [7,8] and seed number [9]. Recent quantitative trait loci (QTL) mapping studies in several agricultural systems have suggested that QTL affecting fitness traits are often clustered in the genome, suggesting tight physical linkage or pleiotropic effects on several components of fitness [9-14], but the ubiquity of this result across experimental populations and environments is not well understood.

Elucidating the genetic architecture underlying fitness traits is critical to predicting their response to selection. For example, although it may be desirable to increase both seed size and seed number in oilseed cultivars, such as Brassica rapa, selection for these traits may be constrained if they are controlled by antagonistic pleiotropy of a single locus or by loci in repulsion phase and in such close physical proximity that recombination is limited. A QTL mapping study in Arabidopsis thaliana suggested that antagonistic pleiotropy of a gene region controlling both seed mass and ovule number may lead to the trade-offs between these traits that have been observed in numerous cultivated and natural species [7]. If this region represents structured pleiotropy, in that allelic substitutions at one locus affect both traits, they will be resistant to change even under strong selection [15,16]. In this situation, it may be prudent to search for natural populations in which these associations are weaker or focus selection efforts on loci that are not pleiotropically regulated.

If covariances among desirable traits are primarily due to extensive linkage disequilibrium, fine-mapping through marker-assisted selection could help to disrupt these associations [17]. QTL regions often encompass regions too large to conclusively discriminate between pleiotropy and linkage, but they provide a starting point for fine-mapping and marker-assisted selection of specific loci [18,19]. Recent studies have advanced from QTL screens to map-based cloning and association mapping to successfully identify candidate genes for important fitness traits in seed oil crops, such as Brassica species [6,20-23].

It is also of interest to understand how QTL expression for fitness characters may vary across simulated versus field/natural settings, as plant breeders attempt to extrapolate to the field from controlled experiments. Significant QTL × environment interactions have been detected for fitness traits in A. thaliana plants grown in the field versus greenhouse or growth chamber environments, suggesting significant differences in QTL expression between controlled and natural settings [24-26]. On the other hand, strong genetic correlations and QTL colocalization between the greenhouse and field have been observed for traits with high heritabilities [26]. Extrapolating QTL data from controlled to field environments may thus be appropriate for some traits.

Seed germination is an important but frequently neglected component of plant fitness. Establishment is a more important determinant of the composition of natural plant populations than is seed availability in many species [27,28]. Germination timing influences fitness by altering the expression of flowering time and other important life-history characters in the annual plant species, Arabidopsis thaliana and Campanulastrum americanum[29-31]. In A. thaliana, seed dormancy QTL colocalize with QTL affecting total lifetime fitness, and epistatic loci affect both traits under field conditions [32]. Clearly, germination phenology influences plant fitness within a generation, but it is equally sensitive to transgenerational influences, which are less well understood, especially in regard to their genetic basis.

Germination phenology is highly sensitive to the maternal genotype, phenotype, and environment [30,33-41]. Maternal flowering day and subsequent timing of seed dispersal can alter germination season in the next generation [34]. Seed size may also affect germination. Studies have found a positive association between seed size and germination percentage or timing in numerous plant species [42], although faster emergence times in smaller seeds have also been demonstrated [43-45]. Despite extensive ecological correlations between seed characters and germination, few studies have mapped QTL for these traits within the same study, and those that have found no QTL affecting both seed size and germination percentage or rate [25,46]. In addition, we are unaware of any study mapping QTL for maternal flowering and progeny germination phenology.

In this study, we examine the quantitative-genetic structure and QTL architecture of bolting time and reproductive fitness characters in recombinant inbred lines of Brassica rapa developed from a cross between a seed oil cultivar and wild-like genotype. QTL mapping was conducted in similar greenhouse and field experiments in order to assess how accurately QTL studies in simulated environments can be extrapolated to the field. In the greenhouse environment, we also examine the relationship between germination phenology and the maternal characters. Brassica rapa is an internationally cultivated seed oil and vegetable crop and a parent of the hybrid-origin mustard crops, B. napus (canola) and B. juncea (Indian mustard) [47]. Numerous genetic resources have been developed for this species [22,23,48,49], and the first reference sequence for B. rapa was released in 2011 [50]. Naturalized populations of B. rapa are distributed throughout the United States and often occur in disturbed habitats or the margins of agricultural fields [51,52]. Wild B. rapa are highly phenotypically variable, ranging from smaller, rapidly-developing spring annuals to larger biennials [52]. Although QTL affecting fruit and seed traits have been identified in a few B. rapa cultivars [53], no study mapping reproductive fitness QTL in a cross between a seed oil cultivar and a wild-like population has been reported.

Methods

Study system

The recombinant inbred lines (RILs) used in this study have been previously described in detail [48,54,55]. Briefly, the RILs were derived from a cross between two inbred B. rapa lines, the seed oil cultivar, yellow sarson (R500), and a rapid-cycling genotype (IMB211). Yellow sarson is a highly inbred, annual B. rapa cultivated in India for over 3,000 years [56]. IMB211 is derived from the Wisconsin Fast Plant™ (WFP) population produced by selection for short generation time [57]. The parental genotypes differ greatly for several reproductive fitness traits. R500 produces the large, heavy, yellow seeds typical of B. rapa cultivars, and IMB211 produces smaller, lighter, brown seeds more similar to naturalized B. rapa populations. IMB211 is phenotypically similar to rapid-developing B. rapa in the United States, and its average flowering day and height are within the range of wild-collected B. rapa[52]. The RILs were developed from self-fertilization of a F1 individual resulting from the R500 × IMB211 cross, and the resulting F2s were propagated by single-seed descent to the S6 generation [48].

Experimental environments

Seeds of 147 RILs and the parental lines were planted in an agricultural field at the University of Minnesota, Saint Paul, MN, April 29–31, 2004. In each of three subplots, four plants of each RIL were planted at 20 cm spacing from each other and the next RIL, resulting in three sets of four adjacent replicates within each RIL. This experiment was part of a larger one examining multiple competitive environments, and further details of the design have been previously described [55]. Bolting date was scored as the date when buds first differentiated from the apical meristem. All flowers were naturally open-pollinated in the field. Plants began to naturally senesce (leaves and stem began to turn brown and dry) in the first week of August and were harvested once 50% of the rosette leaves had senesced on all individuals. At harvest, fruits (siliques) were counted on each plant, and four mature (fully developed) fruits were collected at random from 3–5 replicates of a RIL. Seeds were counted in each fruit, and the length of the longest fruit was measured. Seeds were then pooled within a plant, and seeds-per-fruit, seed mass, seed area, and seed color were determined. Seeds-per-fruit is the mean number of seeds in the 2 longest collected fruit. Seed mass was calculated as the average mass of 20 randomly chosen seeds from the collected fruit. Ten seeds from each plant were then photographed, and their area and color were determined using ImagePro software, which quantifies pixels in digital photographs of the sample (ImagePro Plus V4.5 2001). The trait, seed color, was obtained as the score of the first principal component identified in principal components analysis, using varimax rotation, on pixel densities of three colors: red, green, and blue (PROC FACTOR; SAS, 2001). The first principal component, which was weighted most heavily in red and green, explained greater than 67% of the variation. A lower value for seed color indicates darker, more wild-like seeds.

Fitness traits were also recorded in a greenhouse experiment. Due to space limitations, the greenhouse experiment was conducted over two temporally separated trials, spaced a month apart. In the first trial, on March 14th, 2007, we planted 8 replicates of each of the parental lines and 137 of the 160 B. rapa RILs into 3 in3 pots, filled with Metromix 200 soil, in greenhouses at the University of Minnesota, Saint Paul, MN. Bolting date was recorded for each individual using the same method as in the field, and these plants were then harvested.

Seeds were planted for the second experimental trial on April 17th, 2007 using the same RILs and in the manner described above. Six replicates per line were planted. At 5 days after bolting, plants began flowering, and the number of open flowers was recorded daily. The majority of flowers were allowed to self-pollinate naturally, except that flowers at the 6th and 7th nodes were self-fertilized by hand in order to more closely simulate the effect of pollinators in the field. In early June, the plants began to senesce. At senescence, fruits were counted on each plant, and fruits 6 and 7 (or those fruits arising from hand-pollination) were collected. Fruit length was measured on the longest of the fruits, and the seeds in each collected fruit were counted. Seeds were pooled within a plant, and seed traits were measured as in the field.

Germination characteristics of the seeds collected in the greenhouse experiment were also examined. Seeds collected from the 2 parental lines and RILs with at least 3 replicate plants that set seed were pooled and stored in darkness at room temperature. While excluding RILs with fewer than 3 replicate plants artificially eliminates lines with low viability, this sampling approach helps reduce microenvironmental bias. To further limit microenvironmental maternal effects, the same number of seeds were pooled from each replicate within a RIL combination if possible. Six replicates of 5–10 seeds per RIL were germinated in 3 Conviron E7/2 growth chambers (Controlled Environments Ltd., Winnipeg, Canada) in 2 temporal blocks from February 1–12, 2007 and from March 20–30, 2007. One replicate of each RIL was located in each growth chamber and temporal block, and replicates were randomized within a growth chamber. All seeds were planted onto 0.5% agar in 4.5 cm diameter plastic Petri dishes under a green safe light and then allowed to imbibe water for 24 h in darkness before being moved into the growth chambers. Each growth chamber was fitted with 8 white florescent lights (Silvania Octron, F038/871, 32 W), 2 infrared lights (Industrial Infrared, F32T8/IR-750), and 4 incandescent lights (Sylvania Double Life, Soft White, 120 V/75 W) that generated a R:FR of 1.1 and PAR of 130 – 160 μmol m-2 s-1. Germination trials used 12 L:12D lighting and temperatures of 24D:20 N °C. Germination was checked daily for 10 days and scored when the cotyledons appeared green and open to ~90°. Germination percent and day were calculated by dividing the number of seeds germinated by the seed total and by averaging the germination days for all seeds in a replicate Petri dish, respectively.

Quantitative genetic analysis

For all traits, a restricted maximum likelihood (REML) approach (PROC MIXED, SAS 2001) was used to partition variation attributable to RIL (VL), spatial blocking (VB), and residual error (VR) within each environment. The blocking term was subplot for all traits in the field and greenhouse, and temporal block and chamber for the germination traits measured in the growth chamber. Seed color and fruit length met the assumption of normality in both the field and the greenhouse environments, but all remaining traits were transformed using a Box-Cox procedure [58]. All transformations greatly improved normality and homoscedasticity, and results from the transformed data are reported.

Broad-sense heritabilities were calculated as VL/VP, where, again,VL equals the among-RIL variance component and VP equals the total phenotypic variance for a trait. Least-square means, 95% confidence limits, and best linear unbiased predictor (BLUP) deviations were generated for each RIL using the transformed versions of the traits. BLUP deviations and least-square means were summed, and the resulting values were back-transformed for use in the QTL analysis. The parental genotypes, plants that died before bolting, and RILs with fewer than 3 replicates for a trait were removed from the analyses. Final sample sizes ranged from 730 to 1160 (out of a possible 3600) in the field and 910 to 1400 (out of a possible 2466) in the greenhouse. BLUPs were also used to estimate genetic correlations (rG) within each environment (PROC CORR, SAS 2001). Differences in bivariate matrix structure across environments was tested using Flury’s hierarchical common principal components (CPC) analysis, using the "jump-up approach" and Fisher’s Z-test [59-61]. Cross-environment correlations (rGE) for each trait were calculated as cov1,2/√ (V1,1 × V1,2); cov12 is the covariance of a trait across two environments and V1,1 and V1,2 are the among-RIL variances within each environment [62,63].

QTL-analysis. The R500 × IMB211 RILs were genotyped at 227 RFLP and SSR markers across 10 linkage groups, representing the 10 chromosomes in B. rapa. Marker order for the linkage map was estimated from recombination frequencies observed in the entire population of 160 RILs, i.e., from the most complete dataset, in earlier experiments [48]. Centimorgan distances were then re-estimated in the RILs included in the field and greenhouse experiments presented here (those that survived to bolting) to account for slight differences that might exist due to unsampled recombination events. The map distances for each marker locus were calculated from the estimated recombination frequencies using the Kosambi mapping function in RQTL (R Development Core Team 2010). Twelve markers were uninformative, leaving 215 markers at unique positions in the final re-estimated linkage map.

QTL were mapped using the composite interval mapping (CIM) procedure in QTL Cartographer [64]. To control for effects of variation segregating elsewhere in the genome, we identified co-factors using forward-backward regression, and a 5-cM window; a maximum of 5 co-factors was selected for inclusion in the mapping model. We scanned for QTL at 2-cM intervals across the B. rapa genome. The significance threshold of the likelihood-ratio test statistic (LR) was determined for each trait in each treatment independently by randomly permuting the BLUPs 1000 times. QTL with LRs significant at α = 0.05 are reported [65] with 2-LOD support limits [66]. Additive effects and the genetic variance explained by a QTL were calculated in QTL Cartographer, and confirmed using a GLM model with all QTL detected in the genome-wide screen as main effects [67]. For a given locus, a positive additive effect indicates that the IMB211 parent allele conferred a higher value for the trait. Analysis of variance (ANOVA) was used to test for QTL × environment interactions across the field and greenhouse environments; all significant main-effect QTL for the trait, environment, and QTL × environment effects on each genotypic trait mean were included in the model (PROC GLM, SAS 2001).

Results

Descriptive statistics

RILs varied significantly (P < 0.001) in phenotypic expression of all traits in both environments, indicating underlying genetic variation (Additional file 1: Table S1; Figure  1). Broad-sense heritabilities (H2) were smaller in the field (range, 0.08-0.79) than the greenhouse (range, 0.26-0.86) for all traits except fruit total (Table  1). We detected 7 QTL affecting 6 traits in the field and 12 QTL affecting 8 traits in the greenhouse (Table  2A and B; Figure  2). For any trait, a single QTL explained 5–56% (greenhouse) or 5–59% (field) of the variance (Table  2). Alleles derived from IMB211 had both positive and negative effects on the majority of traits with > 1 QTL in one or both environments, with the exception that the IMB211 allele always decreased the trait value for days to bolting, seed mass, and seed color.

Additional file 1: Table S1. Tests of significance in field versus the greenhouse environments. Asterisks indicate P <0.001.

Format: XLSX Size: 44KB Download fileOpen Data

thumbnailFigure 1. Frequency histograms of best-linear unbiased predictors for reproductive fitness traits in the Brassica rapa R500 × IMB211 recombinant inbred line population field and greenhouse environments. Seed color is the first principal component; a lessor value indicates darker, more wild-like seeds. Germination traits were not measured in the field.

Table 1. Quantitative genetic partitioning of variance components and trait means in the field (F) and greenhouse (GH) environments

Table 2. QTL mapping in the field (A) and greenhouse (B)

thumbnailFigure 2. QTL in the field (filled) and greenhouse (dashed) environments with 1-LOD (bars) and 2-LOD (tails) support limits. Positive additive effects of the IMB211 (wild-like) parent allele are indicated with italics.

Genetic architecture between environments

Single trait genotypic correlations between environments were positive and significant for days to bolting, fruit length, seed mass, and seed color (Table  3). Consistent with the correlations, the majority of QTL for these traits did not exhibit significant QTL × environment interactions (Table  2). Seed color had the highest heritability and cross-environment correlation, and we detected only one large-effect QTL for this trait that was expressed in both environments. These results suggest that the QTL architecture for the majority of fitness traits examined in this study was at least partially conserved between the greenhouse and field.

Table 3. Broad-sense genotypic correlations among traits in the field (above diagonal) and greenhouse (below diagonal) environments

By contrast, genotypic correlations were non-significant for fruit total between environments (Table  3). Three of the five QTL affecting fruit total, including all QTL detected in the field for this trait, were environment-specific, in that the QTL × environment effects were significant (Table  2). QTL × environment effects were also significant for one QTL affecting each of days to bolting, fruit length, seeds-per-fruit, and seed mass, indicating that at least some of the genetic loci affecting each reproductive character, except seed color, were environment-specific or that functional differences between alleles at these QTL were manifested only in one environment.

Relationship among components of fitness

A number of bivariate correlations were conserved between environments. RILs that produced more fruit also produced more seeds-per-fruit, and longer fruits generally had heavier seeds and more seeds (although the fruit length/seeds-per-fruit correlation was significantly stronger in the greenhouse) (Table  3). Plants that produced more seeds-per-fruit also had darker seeds (negative value for seed color) in both environments.

Environmental differences in bivariate correlations were also observed. Plants that produced more fruit and more seeds-per-fruit also had heavier seeds in the field, but these correlations were not significant in the greenhouse (although the difference between environments was not large enough to be statistically significant for the seeds-per-fruit/seed mass correlation). These results suggest that genotypic differences in resource acquisition existed among plants grown in the field, that is, genotypes with high acquisition had higher output for all reproductive components than those with lower acquisition. In the greenhouse, plants did not generally produce both more seeds and larger seeds, suggesting that the genotypic differences in acquisition were smaller in this environment. The relationship between seed mass and color also differed between environments. Although the bivariate correlation between seed mass and color was positive in the greenhouse, the scatter plot showed a sharp increase in seed mass in the darkest colored seeds (Additional file 2: Figure S1), indicating that the lightest and darkest seeds had the largest mass. This relationship was not conserved in the field. Interestingly, significant Z-test differences indicated that the relationship between seed color and seed mass or fruit total reversed between environments, although the relationship between these traits was too weak in each environment for the majority of bivariate correlations to be significant.

Additional file 2: Figure S1. Scatterplot showing best-linear unbiased predictor values (BLUPs) for seed mass and seed color in the field (A) and greenhouse (B). A linear trendline (green) and lowess lines (red) are shown for each figure.

Format: DOCX Size: 260KB Download fileOpen Data

The field-specific positive relationship between fruit total, seeds-per-fruit, and seed mass was consistent with FQTL1-1, which decreased the value for all three traits in the field and had a significant QTL × environment interaction effect for fruit total (Table  2). Moreover, GHQTL9-2 had opposing effects on fruit total versus seed mass, and the QTL × environment effects were significant for these traits (Table  2; Figure  2). These results indicate that GHQTL9-2 conferred environment-specific differences in resource allocation to reproductive fitness traits.

Transgenerational fitness

In the greenhouse environment, we examined the relationships between reproductive characters of the maternal parent and germination patterns in the progeny. Interestingly, later bolting maternal plants produced faster germinating seeds, suggesting a phenological relationship across generations (Table  3). This relationship was consistent with GHQTL10-1, which increased time to germination but decreased time to bolting (Table  2; Figure  2), providing evidence of transgenerational antagonistic pleiotropy at this locus. We also observed a positive relationship between both seed mass and seeds-per-fruit and germination percentage, suggesting that more vigorous maternal genotypes partition more resources to offspring and produce many large seeds that also germinate to a higher percentage.

Discussion

In the current study, we examined the genetic architecture of reproductive fitness components, a topic of interest to both agronomists seeking to increase yield and evolutionary biologists interested in predicting response to selection in quantitative traits. QTL that are consistently expressed among environments are favorable targets for marker-assisted selection in the development of widely distributed agricultural varieties [9]. We detected QTL affecting seeds-per-fruit, seed mass, and fruit length on A1 in both the greenhouse and field environments. These results are consistent with a fruit length QTL previously mapped near the center of A1 in B. rapa[53,68], as well as QTL affecting seed mass [6] and multiple reproductive fitness traits on A1 in B. napus[9,69]. Shi et al. (2009) detected several linked consensus QTL each affecting multiple seed or fruit traits on A1 in B. napus, indicating significant pleiotropy and physical linkage among fitness loci on this chromosome. Given the large influence of sections of A1 on reproductive fitness and yield, this chromosome should be a primary target for fine-mapping and candidate gene screening using association analyses in Brassica species.

Aside from the number and size of seeds, seed oil is a common target of selection among B. rapa varieties, and selection for seed oil phenotypes leads to a differentiation in seed color between cultivated and wild genotypes [70,71]. In the current study, seed color was highly conserved across environments. Heritability and the cross-environment correlation for this trait were both over 0.8, and we detected one large-effect QTL affecting only this trait on A9 in both environments. This supports previous studies mapping seed color to one or a few genes on A9 in B. rapa[23,49]. The candidate gene BrTT8 has been identified in this region and shown to control the accumulation of proanthocyanidins in the seed coat [71,72]. These results suggest that QTL mapping data from greenhouse studies can be extrapolated to field work for loci with a large additive genetic component.

Despite some shared genetic control, at least one QTL for all traits other than seed color was detected in only one environment, suggesting environment-specific effects on fitness traits. Fruit total was largely environment-specific; the genotypic correlation across environments was non-significant for fruit total and seeds-per-fruit, and we detected significant QTL × environment effects for 3 of the 5 QTL affecting fruit total. These results indicate that the genetic architecture of overall reproductive output (total fruit production) is more environmentally sensitive than other seed and fruit traits. A recent meta-analysis of QTL in B. napus that found that 63% of QTL affecting total seed yield were microenvironment-specific versus 47% of QTL for all fitness traits measured [9]. Moreover, 85% of seed yield QTL colocalized with QTL for at least one other life history or reproductive fitness trait [9]. In the current study, all but 1 fruit total QTL affected multiple traits, indicating that different aspects of reproductive output are determined by multiple closely linked or pleiotropic loci.

Fine-mapping studies have found clustered, pleiotropic QTL affecting various life history and fitness characters in rice [5], wheat [73], pea [74], and B. napus[9]. Selection for any one of these QTL is likely to have far reaching consequences on plant fitness. If, for example, QTL clusters correspond to genes regulating coordinated steps in seed development, then selection may have desirable outcomes on several traits, such as an increase in seed size and protein content [74]. Alternatively, antagonistic associations among QTL, as has been found for oil traits in B. juncea[75], may result in negative, nontarget consequences on desirable traits. In our study, a parental allele generally had the same direction of effect on all traits at a QTL, thus supporting coordinated gene regulation of reproductive fitness traits. The only exception was GHQTL9-2, which had opposing effects on fruit total and seed mass in the greenhouse environment, leading to conditional neutrality in the relationship between these traits that differed between the greenhouse and field [76]. Moreover, the opposing effects of GHQTL9-2 on fruit total and seed mass suggests that antagonistic pleiotropy may control allocation of resources to seed size versus seed production in B. rapa, similar to what has been described in A. thaliana[7].

In natural systems, seed germination and establishment are no less important indicators of reproductive fitness than seed traits. Our study provides QTL support for transgenerational co-regulation of life history phenology. Later bolting maternal plants produced seeds that germinated earlier, consistent with opposing effects of SQTL10-1 on these traits. Our results support previous work demonstrating an adaptive association between maternal flowering time and germination season in Campanulastrum americanum[30,34]. Early flowering C. americanum plants produce seeds that disperse earlier and thus germinate earlier because of environmental cues, and strong negative genetic correlations between maternal flowering day and offspring germination timing have been observed in this species [35]. Our results provide QTL support for negative across-generation co-regulation of flowering and germination timing in B. rapa, similar to what has been observed in C. americanum.

The observed negative relationship between days to bolting and germination timing in our study was consistent with opposing effects of SQTL10-1 on these traits, suggesting that in this region there is an important regulatory locus or cluster of loci influencing phenology between generations. SQTL10-1 was flanked by BrFLC1, which is collinear to FLC in A. thaliana[77]. A. thaliana FLC has been shown to act through the flowering time pathway to maternally-regulate germination [78] and is subject to transgenerational epigenetic regulation [78-80]. In B. rapa, a study examining flowering time in a backcross population between biennial and annual B. rapa suggested that most of the variation in flowering time between these life histories was due to FLC loci, including FLC1[77]. Here, flowering time was variable, and it is unknown if SQTL10-1 would have affected germination time if flowering had been synchronous among RILs. FLC1 is a likely candidate gene for transgenerational regulation of evolutionarily important life history phenology in multiple plant species and warrants further investigation of its direct and indirect (through flowering time) effects on germination under natural conditions.

Conclusions

This is the first study to have characterized the correlation structure and QTL architecture of important fitness/yield traits in a crop × wild-like population of B. rapa. Days to bolting, fruit length, and seed color, were highly correlated and displayed QTL-colocalization across environments. Genetic characterization of these traits in controlled settings may readily translate to improvement in field environments. Total fruit number was uncorrelated across environments and most QTL affecting this trait were correspondingly environment-specific, suggesting that the genetic architecture of overall reproductive output is more environmentally sensitive than other measured fitness traits. Finally, BrFLC1 was identified as a potential candidate gene affecting transgenerational regulation of flowering time and seed germination.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JMD, MTB, and CW designed research; JMD and MTB collected and analyzed data; JMD and CW wrote and prepared the manuscript. All authors read and approved the final manuscript.

Acknowledgments

We are grateful to T.C. Osborn and F.L. Iniguez-Luy for development and characterization of the B. rapa RILs used in this study. We also thank J. Johnston, L. Demink, Z. German, C. Willis, A. Hansen, and B. Meyer for their contributions to experimental management and data collection. This work was supported by a grant from the National Science Foundation (20091702) to CW.

References

  1. Allainguillaume J, Alexander M, Bullock JM, Saunders M, Allender CJ, King G, Ford CS, Wilkinson MJ: Fitness of hybrids between rapeseed (Brassica napus) and wild Brassica rapa in natural habitats.

    Mol Ecol 2006, 15:1175-1184. PubMed Abstract | Publisher Full Text OpenURL

  2. Cummings CL, Alexander HM, Snow AA, Rieseberg LH, Kim MJ, Culley TM: Fecundity selection in a sunflower crop-wild study: can ecological data predict crop allele changes ?

    Ecological Applications 2009, 12:1661-1671. OpenURL

  3. Dechaine JM, Burger JC, Chapman M, Seiler GJ, Brunick R, Knapp SJ, Burke JM: Fitness effects and genetic architecture of plant-herbivore interactions in sunflower crop-wild hybrids.

    The New Phytologist 2009, 184:828-841. PubMed Abstract | Publisher Full Text OpenURL

  4. Huang X, Zhao Y, Xinghua W, Li C, Wang A, Zhao Q, Li W, Guo Y, Deng L, Zhu C, Fan D, Lu Y, Weng Q, Liu K, Zhou T, Jing Y, Si L, Dong G, Huang T, Lu T, Feng Q, Qian Q, Li J, Han B: Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm.

    Nature Genetics 2012, 44:32-41. OpenURL

  5. Xue W, Xing Y, Weng X, Zhao Y, Tang W, Wang L, Zhou H, Yu S, Zu C, Li X, Zhang Q: Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice.

    Nature Genetics 2008, 40:761-767. PubMed Abstract | Publisher Full Text OpenURL

  6. Zhang L, Li S, Chen L, Yang G: Identification and mapping of a major dominant quantitative trait locus controlling seeds per silique as a single Mendelian factor in Brassica napus L.

    Theor Appl Genet 2012, 125:695-705. PubMed Abstract | Publisher Full Text OpenURL

  7. Alonso-Blanco C, Blankestijn-De Vries H, Hanhart C, Koornneef M: Natural allelic variation at seed size loci in relation to other life history traits of Arabidopsis thaliana.

    Proc Natl Acad Sci 1999, 96:4710-4717. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Xing Y, Zhang Q: Genetic and molecular basis of rice yield.

    Annu Rev Plant Biol 2010, 61:421-442. PubMed Abstract | Publisher Full Text OpenURL

  9. Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, Bancroft I, Zhao J, Meng J: Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus.

    Genetics 2009, 182:851-861. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Maccaferri M, Sanguineti MC, Corneti S, Ortega JLA, Salem MB, Bort J, DeAmbrogio E, del Moral LFG, Demontis A, El-Ahmed A, Maalouf F, Machlab H, Martos V, Moragues M, Motawaj J, Nachit M, Nserallah N, Ouabbau H, Royo C, Slama A, Tuberosa R: Quantitative trait loci for grain yield and adaptation of Durum Wheat (Triticum durum Desf.) across a wide range of water availability.

    Genetics 2008, 178:489-511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Xie X, Jin F, Song M-H, Suh J-P, Hwang H-G, Kim Y-G, McCouch S, Ahn S-N: Fine mapping of a yield-enhancing QTL cluster associated with transgressive variation in an Oryza sativa x O. rufipogon cross.

    Theor Appl Genet 2008, 116:613-622. PubMed Abstract | Publisher Full Text OpenURL

  12. Khavkin E, Coe EH: Mapped genomic locations for developmental functions and QTLs reflect concerted groups in maize (Zea mays L.).

    Theor Appl Genet 1997, 95:343-352. Publisher Full Text OpenURL

  13. Cai H-Y, Diao S, He Y-G, Zhang L-P, Liu S-J, Zhu Y-G, Zhang Z-H: Genetic and physical mapping of qHY-8, a pleiotropic QTL for heading date and yield-related traits in rice.

    Euphytica 2012, 184:109-118. Publisher Full Text OpenURL

  14. Cai H, Morishima H: QTL clusters reflect character associations in wild and cultivated rice.

    Theor Appl Genet 2002, 104(8):1217-1228. PubMed Abstract | Publisher Full Text OpenURL

  15. Roff DA, Fairbairn DJ: The evolution of trade-offs: where are we?

    J Evol Biol 2007, 20:433-447. PubMed Abstract | Publisher Full Text OpenURL

  16. De Jong G: Quantitative genetics of reaction norms.

    J Evol Biol 1990, 3:447-468. Publisher Full Text OpenURL

  17. Collard BCY, Mackill D: Marker-assisted selection: an approach for precision plant breeding in the twenty-first century.

    Philos Trans R Soc Lond B Biol Sci 2008, 363:557-572. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Asins MJ: Present and future of quantitative trait locus analysis in plant breeding.

    Plant Breeding 2002, 121:281-291. Publisher Full Text OpenURL

  19. Cooper M, van Eeuwijk FA, Hammer GL, Podlich DW, Messina C: Modeling QTL for complex traits: detection and context for plant breeding.

    Curr Opin Plant Biol 2009, 12:231-240. PubMed Abstract | Publisher Full Text OpenURL

  20. Bisht NC, Gupta V, Ramchiary N, Sodhi YS, Mukhopadhyay A, Arumugam N, Pental D, Pradhan AK: Fine mapping of loci involved with glucosinolate biosynthesis in oilseed mustard (Brassica juncea) using genomic information from allied species.

    Theor Appl Genet 2009, 118:413-421. PubMed Abstract | Publisher Full Text OpenURL

  21. Huang Z, Ban YY, Yang L, Zhang Y, Li HQ, Xiao ES, Xu AX, Zhang DH: Fine mapping of the yellow seed locus in Brassica juncea L.

    Genome 2012, 55:8-14. PubMed Abstract | Publisher Full Text OpenURL

  22. Zhang JF, Lu Y, Yuan YX, Zhang XW, Geng JF, Chen Y, Cloutier S, McVetty PBE, Li GY: Map-based cloning and characterization of a gene controlling hairiness and seed coat color traits in Brassica rapa.

    Plant MolBiol 2009, 69:553-563. OpenURL

  23. Xiao L, Zhao Z, Du DZ, Yao YM, Xu L, Tang GY: Genetic characterization and fine mapping of a yellow-seeded gene in Dahuang (a Brassica rapa landrace).

    Theor Appl Genet 2012, 124:903-909. PubMed Abstract | Publisher Full Text OpenURL

  24. Weinig C, Ungerer MC, Dorn LA, Kane NC, Toyonaga Y, Halldorsdottir SS, Mackay TFC, Purugganan MD, Schmitt J: Novel loci control variation in reproductive timing in Arabidosis thaliana in natural environments.

    Genetics 2002, 162:1875-1884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Malmberg RL, Held S, Waits A, Mauricio R: Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse.

    Genetics 2005, 171:2015-2027. OpenURL

  26. Brock MT, Dechaine JM, Iniguez-Luy FL, Maloof JN, Stinchcombe JR, Weinig C: Floral genetic architecture: an examination of QTL architecture underlying floral (co)variation across environments.

    Genetics 2010, 186:1451-1465. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Clark CJ, Poulsen JR, Levey DJ, Osenberg CW: Are plant populations seed limited? A critique and meta-analysis of seed addition experiments.

    The American Naturalist 2007, 170:128-142. PubMed Abstract | Publisher Full Text OpenURL

  28. Turnbull LA, Crawley MJ, Rees M: Are Plant Populations Seed-Limited? A Review of Seed Sowing Experiments M INI- REVIE of seed sowing Are plant populations experiments.

    Oikos 2000, 88:225-238. Publisher Full Text OpenURL

  29. Donohue K: Seeds and seasons: interpreting germination timing in the field.

    Seed Science Research 2005, 15:175-187. Publisher Full Text OpenURL

  30. Galloway L, Etterson J: Transgenerational plasticity is adaptive in the wild.

    Science 2007, 318:1134-1136. PubMed Abstract | Publisher Full Text OpenURL

  31. Donohue K: Germination timing influences natural selection on life-history characters in Arabidopsis Thaliana.

    Ecology 2002, 83:1006-1016. Publisher Full Text OpenURL

  32. Huang X, Schmitt J, Dorn L, Griffith C, Effgen S, Takao S, Koornneef M, Donohue K: The earliest stages of adaptation in an experimental plant population: strong selection on QTLS for seed dormancy.

    Molecular Ecology 2010, 19:1335-1351. PubMed Abstract | Publisher Full Text OpenURL

  33. Donohue K, Schmitt J: Maternal environmental effects in plants: adaptive plasticity? In Maternal Effects as Adaptations. Edited by Mousseau T, Fox C. New York, NY: Oxford University Press; 1998:137-158. OpenURL

  34. Galloway L, Burgess K: Manipulation of flowering time: phenological integration and maternal effects.

    Ecology 2009, 90:2139-2148. PubMed Abstract | Publisher Full Text OpenURL

  35. Galloway LF, Etterson JR, McGlothlin JW: Contribution of direct and maternal genetic effects to life-history evolution.

    The New Phytologist 2009, 183:826-838. PubMed Abstract | Publisher Full Text OpenURL

  36. Dechaine JM, Gardner G, Weinig C: Phytochromes differentially regulate seed germination responses to light quality and temperature cues during seed maturation.

    Plant, Cell & Environment 2009, 32:1297-1309. Publisher Full Text OpenURL

  37. Donohue K: Completing the cycle: maternal effects as the missing link in plant life histories.

    Philos Trans R Soc Lond B Biol Sci 2009, 364:1059-1074. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Herman JJ, Sultan SE, Horgan-Kobelski T, Riggs C: Adaptive transgenerational plasticity in an annual plant: grandparental and parental drought stress enhance performance of seedlings in dry soil.

    Integr Comp Bio 2012, 52:77-88. Publisher Full Text OpenURL

  39. Contreras S, Bennett MA, Metzger JD, Tay D: Maternal light environment during seed development affects lettuce seed weight, germability, and storability.

    HortScience 2008, 43:845-852. OpenURL

  40. Gutterman Y: Maternal effects on seeds during development. In SEEDS: The Ecology of Regeneration in Plant Communities. Edited by Fenner M. Melksham, UK: Redwood Press Ltd; 2000:59-84. OpenURL

  41. Roach DA, Wulff RD: Maternal effects in plants.

    Annual Review of Ecology and Systematics 1987, 18:209-235. Publisher Full Text OpenURL

  42. van Molken T, Jorritsma-Wienk LD, van Hoek PHW, de Kroon H: Only size matters for germination in different populations of the dimorphic Tragopogon pratensis subsp. pratensis (Asteraceae).

    Am J Bot 2005, 2004:432-437. OpenURL

  43. Stamp NE: Production and effect of seeds size in a grassland annual (Erodium brachycarpum, Geraniaceae).

    Am J Bot 1990, 77:874-882. Publisher Full Text OpenURL

  44. Susko DJ, Lovett-Doust L: Patterns of seed mass variation and their effects on seedling traits in Alliaria petiolata (Brassicaceae).

    Canadian J Bot 2000, 87:56-66. Publisher Full Text OpenURL

  45. Zhang J: Seed dimorphism in relation to germination and growth of Cakile edentula.

    Canadian Journal of Botany 1993, 71:1231-1235. Publisher Full Text OpenURL

  46. Bettey M, Finch-Savage WE, King GJ, Lynn JR: Quantitative genetic analysis of seed vigour and pre-emergence seedling growth traits in Brassica oleracea.

    New Phytologist 2000, 148:277-286. Publisher Full Text OpenURL

  47. Prakash S, Hinata K: Taxonomy, cytogenetics and origin of crop Brassicas.

    Opera Botanica 1980, 55:180-188. OpenURL

  48. Iniguez-Luy FL, Lukens L, Farnham MW, Amasino RM, Osborn TC: Development of public immortal mapping populations, molecular markers and linkage maps for rapid cycling Brassica rapa and B. oleracea.

    Theor Appl Genet 2009, 120:31-43. PubMed Abstract | Publisher Full Text OpenURL

  49. Kebede B, Bheema K, Greenshields D, Li C, Selvaraj G, Rahman H: Construction of a genetic linkage map and mapping of QTL for seed color in Brassica rapa.

    Genome 2012, 55:813-823. PubMed Abstract | Publisher Full Text OpenURL

  50. Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, Bai Y, Mun J, Bancroft I, Cheng F: The genome of the mesopolyploid crop species Brassica rapa.

    Nat Genet 2011, 43:1035-1039. PubMed Abstract | Publisher Full Text OpenURL

  51. Dorn L, Mitchell-Olds T: Genetics of Brassica campestris. 1. Genetic constraints on evolution of life-history characters.

    Evolution 1991, 45:371-379. Publisher Full Text OpenURL

  52. Mitchell-Olds T: Pleiotropy causes long-term genetic constraints on life-history evolution in Brassica rapa.

    Evolution 1996, 50:1849-1858. Publisher Full Text OpenURL

  53. Lou P, Zhao J, Kim JS, Shen S, Del Carpio DP, Song X, Jin M, Vreugdenhil D, Wang X, Koornneef M, Bonnema G: Quantitative trait loci for flowering time and morphological traits in multiple populations of Brassica rapa.

    J Exp Bot 2007, 58:4005-4016. PubMed Abstract | Publisher Full Text OpenURL

  54. Lou P, Xie Q, Xu X, Edwards CE, Brock MT, Weinig C, McClung CR: Genetic architecture of the circadian clock and flowering time in Brassica rapa.

    Theor Appl Genet 2011, 123:397-409. PubMed Abstract | Publisher Full Text OpenURL

  55. Dechaine JM, Johnston JA, Brock MT, Weinig C: Constraints on the evolution of adaptive plasticity: costs of plasticity to density are expressed in segregating progenies.

    New Phytologist 2007, 176:874-882. PubMed Abstract | Publisher Full Text OpenURL

  56. Hinata K, Prakash S: Ethnobotany and evolutionary origin of Indian oleiferous Brassicae.

    Indian Journal of Genetics and Plant Breeding 1984, 44:102-112. OpenURL

  57. Williams PH, Hill CB: Rapid-cycling populations of Brassica.

    Science 1986, 232:1385-1389. PubMed Abstract | Publisher Full Text OpenURL

  58. Box GEP, Cox DR: An analysis of transformations.

    Journal of the Royal Statistical Society, Series B 1964, 26:211-243. OpenURL

  59. Phillips PC, Arnold SJ: Hierarchical comparison of genetic variance-covariance matrices. I. Using the flury hierarchy.

    Evolution 1999, 53:1506-1515. Publisher Full Text OpenURL

  60. Flury B: Common Principal Components and Related Multivariate Models. New York: John Wiley & Sons; 1998. OpenURL

  61. Mezey JG, Houle D: Comparing G matrices: Are common principal components informative?

    Genetics 2003, 165:411-425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Gurganus MC, Fry JD, Nuzhdin SV, Pasyukova EG, Lyman RF, Mackay TFC: Genotype-environment interaction at quantitative trait loci affecting sensory bristle number in Drosophila melanogaster.

    Genetics 1998, 149:1883-1898. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Robertson A: The sampling variance of genetic correlations.

    Biometrics 1959, 15:469-485. Publisher Full Text OpenURL

  64. Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer Version 2.5. NCSU Statistical Genetics: Raleigh; 2007. OpenURL

  65. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping.

    Genetics 1994, 138:963-971. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. van Ooijen JW: Accuracy of mapping quantitative trait loci in autogamous species.

    Theor Appl Genet 1992, 84:803-811. PubMed Abstract OpenURL

  67. Lynch M, Walsh JB: Genetics and Analysis of Quantitative Traits. Sunderland: Sinauer Associates Inc.; 1998. OpenURL

  68. Mahmood T, Rahman MH, Stringam GR, Yeh F, Good AG: Identification of quantitative trait loci (QTL) for oil and protein contents and their relationships with other seed quality traits in Brassica juncea.

    Theor Appl Genet 2006, 113:1211-1220. PubMed Abstract | Publisher Full Text OpenURL

  69. Yang P, Shu C, Chen L, Xu J, Wu J, Liu K: Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.

    Theor Appl Genet 2012, 125:285-296. PubMed Abstract | Publisher Full Text OpenURL

  70. Chen BY, Heneen WK: Inheritance of seed colour in Brassica campestris L. and breeding for yellow-seeded B. napus L.

    Euphytica 1992, 59:157-163. Publisher Full Text OpenURL

  71. Shirzadegan M, Röbbelen G: Influence of Seed Color and Hull Proportion on Quality Properties of Seeds in Brassica napus L.

    Fette, Seifen, Anstrichmittel 1985, 87:235-237. Publisher Full Text OpenURL

  72. Chen L, Hong M, Zhang Y, Zu F, Wen J, Yi B, Ma C, Shen J, Tu J, Fu T: A Large Insertion in bHLH Transcription Factor BrTT8 Resulting in Yellow Seed Coat in Brassica rapa.

    Plos One 2012, 7:e44145.

    doi:44110.41371/journal.pone.0044145

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  73. Quarrie S, Pekic Quarrie S, Radosevic R, Rancic D, Kaminska A, Barnes J, Leverington M, Ceoloni C, Dodig D: Dissecting a wheat QTL for yield present in a range of environments: from the QTL to candidate genes.

    J Exp Bot 2006, 57:2627-2637. PubMed Abstract | Publisher Full Text OpenURL

  74. Burstin J, Marget P, Huart M, Moessner A, Mangin B, Duchene C, Desprez B, Munier-Jolain N, Duc G: Developmental genes have pleiotropic effects on plant morphology and source capacity, eventually impacting on seed protein content and productivity in pea.

    Plant Physiol 2007, 144:768-781. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Ramchiary N, Padmaja KL, Sharma S, Gubta V, Sodhi YS, Mukhopadhyay A, Arumugam N, Pental D, Pradhan AK: Mapping for yield influencing QTL in Brassica juncea: implications for breeding of a major oilseed crop of dryland areas.

    Theor Appl Genet 2007, 115:807-817. PubMed Abstract | Publisher Full Text OpenURL

  76. Fry JD, Nuzhdin SV, Pasyukova EG, Mackay TF: QTL mapping of genotype-environment interaction for fitness in Drosophila melanogaster.

    Genet Res 1998, 71:133-141. PubMed Abstract | Publisher Full Text OpenURL

  77. Schranz ME, Quijada P, Sung SB, Lukens L, Amasino RM, Osborn TC: Characterization and effects of the replicated flowering time gene FLC in Brassica rapa.

    Genetics 2002, 162:1457-1468. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Chiang GCK, Barua D, Kramer EM, Amasino RM, Donohue K: Major flowering time gene, flowering locus C, regulates seed germination in Arabidopsis thaliana.

    Proc Natl Acad Sci U S A 2009, 106:11661-11666. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Bastow R, Mylne J, Lister C, Lippman Z, Martienssen R, Dean C: Vernalization requires epigenetic silencing of FLC by histone methylation.

    Nature 2004, 427:164-167. PubMed Abstract | Publisher Full Text OpenURL

  80. Sheldon C, Hills M, Lister C, Dean C, Dennis E, Peacock W: Resetting of FLOWERING LOCUS C expression after epigenetic repression by vernalization.

    Proc Natl Acad Sci 2008, 105:2214-2219. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL