Genetic dynamics underlying phenotypic development of biomass yield in triticale
© Liu et al.; licensee BioMed Central Ltd. 2014
Received: 25 November 2013
Accepted: 6 June 2014
Published: 10 June 2014
The nature of dynamic traits with their phenotypic plasticity suggests that they are under the control of a dynamic genetic regulation. We employed a precision phenotyping platform to non-invasively assess biomass yield in a large mapping population of triticale at three developmental stages.
Using multiple-line cross QTL mapping we identified QTL for each of these developmental stages which explained a considerable proportion of the genotypic variance. Some QTL were identified at each developmental stage and thus contribute to biomass yield throughout the studied developmental phases. Interestingly, we also observed QTL that were only identified for one or two of the developmental stages illustrating a temporal contribution of these QTL to the trait. In addition, epistatic QTL were detected and the epistatic interaction landscape was shown to dynamically change with developmental progression.
In summary, our results reveal the temporal dynamics of the genetic architecture underlying biomass accumulation in triticale and emphasize the need for a temporal assessment of dynamic traits.
Quantitative trait locus (QTL) mapping approaches are popular genomic tools to dissect the genetic architecture underlying complex traits and to identify QTL . The basis for these approaches are mapping populations, which are commonly phenotyped at only a single time point. It is well known, however, that many traits of agronomic importance are under the control of complex dynamic regulation . Consequently, the traditional static examination completely neglects the developmental dynamics underlying trait formation. For example, biomass changes with time and thus, two genotypes can be phenotypically identical at a certain developmental stage while the temporal patterns of genetic control may vary between them. Until now, only little is known about the genetic dynamics of complex traits. Yan et al.  evaluated plant height in rice at different developmental stages and observed some QTL that could be detected at all stages whereas others were only detectable at one or some of them. Busemeyer et al.  recently reported a dynamic mapping for biomass in triticale and also observed developmental stage specific QTL.
In order to assess the temporal changes in the genetic control of trait formation, the phenome of the plants must be assessed at several time points. A key component for monitoring the phenotypic changes of plants is the development of appropriate phenotyping technologies which for crops must enable phenotyping under field conditions [5, 6]. Busemeyer et al.  have recently described the development of the ‘BreedVision’ precision phenotyping platform for non-invasive, high-throughput and high-dimensional phenotyping of small grain cereals under field conditions. This platform incorporates light curtains, laser distance sensors, a 3D-Time-of-Flight camera and hyperspectral imaging and traits are predicted based on sensor fusion, i.e., the combination of sensors and their information. This platform is of particular interest as it enables the evaluation of traits not amenable to traditional phenotyping and in addition, can replace destructive measurements thereby enabling multiple measurements of the same plants. Biomass, for example, is traditionally evaluated by harvesting yield plots with a field chopper. The major disadvantage of this approach is, that it is destructive, thus prohibiting the assessment of other traits as well as the biomass development of the plots over time. The precision phenotyping platform has been calibrated for biomass yield of triticale which yielded a high prediction accuracy and heritability for the predicted biomass .
Triticale (× Triticosecale Wittmack L.) shows a broad variation for biomass yield  and is therefore well suited to study the genetics underlying this trait in small grain cereals. A powerful approach for QTL mapping is to employ multiple segregating families [9, 10]. Busemeyer et al.  used biomass data from four families predicted based on measurements of the precision phenotyping platform at three developmental stages in combination with an association mapping approach to detect QTL. Association mapping is an identity-by-state approach and the employed biometric model has been shown to be rather conservative . Consistenly, the analysis of part of these data with an identity-by-descent approach identified substantially more QTL . The aim of this study was therefore to reanalyse the triticale biomass data  employing an identity-by-descent approach to assess the temporal dynamics of QTL contributions to biomass yield. Here, we describe the detection of main effect QTL for biomass yield in triticale at three different developmental stages, the variation in the contribution of these QTL to the genotypic variance over time, and the dynamics of epistatic QTL contributing to the genetic architecture underlying phenotypic development of biomass yield.
Results of QTL mapping and fivefold cross-validation
Many traits of biological or agronomic importance undergo dynamic changes with time and the developmental progression of the individuals. It therefore appears obvious that these changes will at least in part be reflected by similar changes in the underlying genetics. However, a temporal assessment of the genetic architecture of complex traits has thus far largely been neglected. In this study we therefore performed a dynamic mapping of biomass yield QTL in triticale to decipher the genetic dynamics underlying phenotypic development of this trait.
Phenotypic development of biomass yield
A prerequiste for efficient QTL detection in multiple families is the precise estimation of phenotypic values . In addition, for traits which are traditionally assessed by destructive measurements, alternative non-invasive approaches are required to enable monitoring dynamic changes of the trait over time. In this study we focused on biomass yield as the accumulation of biomass is central to agricultural productivity, employing triticale as model crop for small grain cereals. The earliest developmental stage (BM1) was chosen at BBCH stage 49 as this represents the earliest time point when triticale would be harvested for biomass. Biomass was predicted based on sensor measurements of the precision phenotyping platform and we observed high heritabilities for all three developmental stages at which biomass yield was assessed. This in combination with the high prediction accuracies underlines the great potential of precision phenotyping, especially for a temporal assessment of dynamic traits.
Biomass yield was characterized by a progression with time and the development of the plants (Figure 1). The correlations of biomass yield at the different time points were highest for the directly successive time points (BM1-BM2, BM2-BM3) and decreased with increasing temporal and thus developmental distance between time points (BM1-BM3). The variation in phenotypic values increased with development of the plants (Figure 1, Additional file 1: Table S1). This does not appear to be caused by an increased influence of the environment as the contribution of the genotype-by-environment interaction variance was strongest at BM1. Rather, the effects of the genetic factors contributing to biomass yield become more pronounced the longer they act on the trait. This is exemplified by the development of biomass in family EAW78. Gowda et al.  have recently shown that grain yield, heading time, spikes per square meter, 1000-kernel weight, and early plant height are key contributors to early biomass yield. Family EAW78 segregates for a major plant height QTL which likely corresponds to the rye dwarfing gene Ddw1 and which causes the observed bimodal distribution . However, the distribution of the biomass yield values in family EAW78 shows a strong temporal plasticity. At BM1 the biomass values are normally distributed, then turn into the bimodel distribution at BM2 which at BM3 becomes less pronounced. Thus, Ddw1 appears to have the strongest contribution to biomass yield at BM2 after which its influence on biomass declines or that of other factors becomes more prominent. This illustrates the plasticity of the trait and the temporal contribution of genetic factors to its expression and highlights the need for a temporal assessment of the underlying genetics.
Detection of biomass yield QTL
Busemeyer et al.  used an association mapping approach to analyze the data but identified only the major QTL on chromosomes 5A and 5R. By contrast, the linkage mapping approach applied here appears more powerful with regard to QTL detection in this data set as we identified nine or ten QTL for biomass yield at all three developmental stages. These QTL were supported by the QTL frequency distributions (Table 1, Figure 2, Additional file 1: Figure S1). Consequently, the proportion of genotypic variance (p G ) explained by the QTL detected here was approximately twice as high as that described in the previous study. For all three developmental stages the p G was around 60 percent and even cross-validated still amounted to a considerable 30 to 39 percent (Table 1). This is in accordance with the quantitative nature of the trait which implies many QTL with effects too small to be detected in QTL mapping given a reasonable population size.
Contribution of epistasis to biomass yield
Epistasis refers to interactions between the alleles at two or more genetic loci . It contributes to the genetic architecture of many complex traits and has recently been reported for different crops including maize, wheat and rapeseed [10, 13, 15–17]. The orthogonal contrasts between family means and the means of the respective parents were only significant for BM3 in family EAW78 which indicates the presence of epistasis. It must be noted however, that non-significant orthogonal contrasts do not imply abscence of epistasis. Our analysis revealed epistatic QTL for all three developmental stages (Figure 5). The detection power for epistatic QTL, however, more strongly depends on the population size than that for the detection of main effect QTL. Thus, despite the relatively large mapping population employed here, many epistatic QTL may have remained undetected due to insufficient QTL detection power. In contrast to the main effect QTL for which similar numbers were detected for each developmental stage, we identified substantially more epistatic QTL for BM1 as compared to BM2 and BM3. This corroborates previous findings from Busemeyer et al.  and suggests a varying contribution of epistasis to the genetic architecture of biomass yield over time. The proportion of genotypic variance explained by these epistatic QTL was rather small and averaged 2.4, 1.2 and 1.5 percent for BM1, BM2 and BM3, respectively. Nevertheless, this result illustrates the contribution of epistasis to the genetic architecture of biomass yield, especially assuming that a number of epistatic QTL had remained undetected.
Temporal genetic patterns of biomass regulation
The nature of dynamic traits with their temporal changes suggests at least some plasticity in the underlying genetic control. We observed an almost equal number of QTL for biomass yield for each of the three developmental stages but the Venn diagram indicated both specific and overlapping QTL (Figures 2 and 3). The detection of a QTL at all three time points indicates that it stably contributes to biomass yield throughout the studied development, ranging from the stage where the awns are just visible to very early dough development. By contrast, some QTL were only detected at one or two of these developmental stages. Two QTL were identified at two developmental stages which in both cases were successive stages, i.e., BM1 and BM2 or BM2 and BM3. For each developmental stage we identified four or five QTL that were only detected for that particular stage. As illustrated by the LOD profiles (Figure 2) this does not appear to be caused by peaks close to the significance threshold, being by chance slightly above or below it for the different time points. By contrast, it appears to reflect the dynamic genetics underlying biomass accumulation in triticale. This is supported by the variable contribution of the loci detected for any of the three developmental stages to the genotypic variance at BM1, BM2 and BM3, respectively (Figure 4). Interestingly, we observed a similar plasticity for the epistatic interactions suggesting that the entire genetic architecture of biomass yield undergoes dynamic changes during the developmental progression of the plants.
In this study, we employed phenomics data of biomass yield generated at three developmental stages by precision phenotyping of a large mapping population of triticale. We show the phenotypic plasticity of this trait and demonstrate that this is reflected in a similar plasticity of the underlying genetics. Thus, the genetics of dynamic traits should best be assessed in a temporal manner to capture all the genetic factors that contribute to the trait during development.
Plant material, field trials and phenotypic data
The plant material, field trials and the collection of phenotypic data used in this study have been described by Busemeyer et al. . In brief, phenotypic data for biomass yield were obtained by non-invasive prediction based on a precision phenotyping platform . A calibration experiment was performed based on 25 diverse triticale genotypes at three developmental stages: BM1 = BBCH stage 49 (awns visible), BM2 = BBCH 69 (late flowering), and BM3 = BBCH 81 (very early dough development) . The plants were harvested with a field chopper to determine reference fresh weight based on which the calibration models were established. With these calibration models biomass yield was predicted at the three developmental stages in a mapping population consisting of 647 doubled haploid (DH) [19, 20] triticale lines. The results for BM3 have in part been reported by Alheit et al. . The population consisted of four families designated DH06 (131), DH07 (120), EAW74 (200), and EAW78 (196) which have been described by Alheit et al. . The DH lines were grown in partially replicated designs  including common checks with 960 plots per location, at two locations in two years. Phenotypic data were analyzed by ordinary lattice analysis of variance . Variance components were determined by the restricted maximum likelihood (REML) method assuming a full random model and heritability (h 2 ) on an entry-mean basis was estimated from the variance components as the ratio of genotypic to phenotypic variance . Best linear unbiased estimates (BLUEs) were estimated across environments assuming fixed effects for the genotype. All statistical analyses were performed using ASReml 3.0 .
Multiple-line cross QTL mapping
where Y was a N × 1 column vector of the BLUE values of phenotypic data of N progenies coming from P families. J was a N × P matrix whose elements were 1 or 0 according to whether or not individual i belonged to family p and M was a P × 1 vector of family specific means. Xq (Xc) a N × P matrix containing the expected number (ranging from 0 to 2) of allele k for each individual in family p at QTL q (cofactor c), and Bq (Bc) was a P × 1 vector of the expected allele substitution effects of QTL q (cofactor c) in family p. ϵ was the vector of the residuals.
Cofactor selection was performed using PROC GLMSELECT implemented in the statistical software SAS . The presence of a putative QTL in an interval was tested using a likelihood-ratio test with the statistical software R . LOD-thresholds of 4.8 for BM1, 4.6 for BM2 and 4.6 for BM3 were used corresponding to an experiment-wise type I error of P < 0.10, based on 2000 permutations . Cofactors were excluded within a distance to the marker interval under consideration smaller than 10 cM and the support interval of a QTL was defined as a LOD fall-off of 1.0 expressed as position on the chromosome in centimorgans (cM) . The proportion of genotypic variance explained by the detected QTL was estimated as R 2 adj /h 2 . Biomass yield QTL were declared as overlapping between the three developmental stages if they fell within an arbitrarily defined 10 cM interval surrounding the QTL. Fivefold cross-validation was done as described by Liu et al. .
The epistasis scan for pairwise interactions was done with the model described above which was extended by the term X q’ B q’ for the second locus and the interaction term between the two loci q and q’ X qq’ B qq’ . We used an α-level of 0.05 and followed the suggestion of Holland et al.  dividing the α-level by the number of possible independent pairwise interactions between chromosome regions, assuming two separate regions per chromosome (P < 5.3e-5). The circular plots illustrating the epistatic interactions were created with Circos . The detected epistatic QTL were illustrated by these circular plots (Figure 5a) showing the interactions among different chromosome regions. The temporal development of one QTL detected at BM2 involving chromosomes 6A and 5B was illustrated by the 3-dimensional interaction landscape between these two chromosomes. In Figure 5b the –log10(P values) of all tested pairwise interactions were plotted and consequently the higher the peak, the stronger the association of the epistatic interaction with the trait. Significant interactions are marked in red.
This research was funded by the German Federal Ministry of Education and Research (BMBF) under the promotional reference 0315414. This publication only reflects the views of the authors. We acknowledge the handling of the funding by the Project Management Organisation Jülich (PtJ). We thank Lucas Busemeyer, Kim Möller and Katharina Alheit for their high dedication while recording and generating the phenotypic data, and Agnes Rölfing-Finze, Hans Häge, Jacek Till and Justus von Kittlitz for outstanding work in the field.
- Würschum T: Mapping QTL for agronomic traits in breeding populations. Theor Appl Genet. 2012, 125: 201-210. 10.1007/s00122-012-1887-6.PubMedView ArticleGoogle Scholar
- Wu R, Lin M: Functional mapping – how to map and study the genetic architecture of dynamic camplex traits. Nat Rev Genet. 2006, 7: 229-237.PubMedView ArticleGoogle Scholar
- Yan J, Zhu J, He C, Benmoussa M, Wu P: Molecular dissection of developmental behaviour of plant height in rice (Oryza sativa L.). Genetics. 1998, 150: 1257-1265.PubMed CentralPubMedGoogle Scholar
- Busemeyer L, Ruckelshausen A, Möller K, Melchinger AE, Alheit KV, Maurer HP, Hahn V, Weissmann EA, Reif JC, Würschum T: Precision phenotyping of biomass accumulation in triticale reveals temporal genetic patterns of regulation. Sci Rep. 2013, 3: 2442-PubMed CentralPubMedView ArticleGoogle Scholar
- Montes JM, Melchinger AE, Reif JR: Novel throughput phenotyping platforms in plant genetic studies. Trends Plant Sci. 2007, 12: 433-436. 10.1016/j.tplants.2007.08.006.PubMedView ArticleGoogle Scholar
- White JW, Andrade-Sanchez P, Gore MA, Bronson KF, Coffelt TA, Conley MM, Feldmann KA, French AN, Heun JT, Hunsaker DJ, Jenks MA, Kimball BA, Roth RL, Strand RJ, Thorp KR, Wall GW, Wang G: Field-based phenomics for plant genetics research. Field Crops Res. 2012, 133: 101-112.View ArticleGoogle Scholar
- Busemeyer L, Mentrup D, Möller K, Wunder E, Alheit K, Hahn V, Maurer HP, Reif JC, Würschum T, Müller J, Rahe F, Ruckelshausen A: Breedvision - A multi-sensor platform for non-destructive field-based phenotyping in plant breeding. Sensors (Switzerland). 2013, 13: 2830-2847. 10.3390/s130302830.View ArticleGoogle Scholar
- Gowda M, Hahn V, Reif JC, Longin CFH, Alheit KV, Maurer HP: Potential for simultaneous improvement of grain and biomass yield in Central European winter triticale germplasm. Field Crops Res. 2011, 121: 153-157. 10.1016/j.fcr.2010.12.003.View ArticleGoogle Scholar
- Steinhoff J, Liu W, Maurer HP, Würschum T, Longin FH, Ranc N, Reif JC: Multiple-line cross quantitative trait locus mapping in european elite maize. Crop Sci. 2011, 51: 2505-2516. 10.2135/cropsci2011.03.0181.View ArticleGoogle Scholar
- Steinhoff J, Liu W, Reif JC, Della Porta G, Ranc N, Würschum T: Detection of QTL for flowering time in multiple families of elite maize. Theor Appl Genet. 2012, 125: 1539-1551. 10.1007/s00122-012-1933-4.PubMedView ArticleGoogle Scholar
- Würschum T, Liu W, Gowda M, Maurer HP, Fischer S, Schechert A, Reif JC: Comparison of biometrical models for joint linkage association mapping. Heredity. 2012, 108: 332-340. 10.1038/hdy.2011.78.PubMed CentralPubMedView ArticleGoogle Scholar
- Alheit KV, Busemeyer L, Liu W, Maurer HP, Gowda M, Hahn V, Weissmann S, Ruckelshausen A, Reif JC, Würschum T: Multiple-line cross QTL mapping for biomass yield and plant height in triticale (x Triticosecale Wittmack). Theor Appl Genet. 2013, 127: 251-260.PubMedView ArticleGoogle Scholar
- Liu W, Reif JC, Ranc N, Porta GD, Würschum T: Comparison of biometrical approaches for QTL detection in multiple segregating families. Theor Appl Genet. 2012, 125: 987-998. 10.1007/s00122-012-1889-4.PubMedView ArticleGoogle Scholar
- Carlborg Ö, Haley CS: Epistasis: Too often neglected in complex trait studies?. Nat Rev Genet. 2004, 5: 618-625. 10.1038/nrg1407.PubMedView ArticleGoogle Scholar
- Buckler ES, Holland JB, Bradbury PJ, Acharya CB, Brown PJ, Browne C, Ersoz E, Flint-Garcia S, Garcia A, Glaubitz JC, Goodman MM, Harjes C, Guill K, Kroon DE, Larsson S, Lepak NK, Li H, Mitchell SE, Pressoir G, Peiffer JA, Rosas MO, Rocheford TR, Romay MC, Romero S, Salvo S, Villeda HS, Da Silva HS, Sun Q, Tian F, Upadyayula N, et al: The genetic architecture of maize flowering time. Science. 2009, 325: 714-718. 10.1126/science.1174276.PubMedView ArticleGoogle Scholar
- Reif JC, Maurer HP, Korzun V, Ebmeyer E, Miedaner T, Würschum T: Mapping QTLs with main and epistatic effects underlying grain yield and heading time in soft winter wheat. Theor Appl Genet. 2011, 123: 283-292. 10.1007/s00122-011-1583-y.PubMedView ArticleGoogle Scholar
- Würschum T, Maurer HP, Dreyer F, Reif JC: Effect of inter- and intragenic epistasis on the heritability of oil content in rapeseed (Brassica napus L.). Theor Appl Genet. 2013, 126: 435-441. 10.1007/s00122-012-1991-7.PubMedView ArticleGoogle Scholar
- Lancashire PD, Bleiholder H, van Boom TD, Langelüddeke P, Stauss R, Weber E, Witzenberger A: A uniform decimal code for growth stages of crops and weeds. Annals of Applied Biology. 1991, 119: 561-601. 10.1111/j.1744-7348.1991.tb04895.x.View ArticleGoogle Scholar
- Würschum T, Tucker MR, Reif JC, Maurer HP: Improved efficiency of doubled haploid generation in hexaploid triticale by in vitro chromosome doubling. BMC Plant Biol. 2012, 12: 109-10.1186/1471-2229-12-109.PubMed CentralPubMedView ArticleGoogle Scholar
- Würschum T, Tucker MR, Maurer HP: Stress treatments influence efficiency of microspore embryogenesis and green plant regeneration in hexaploid triticale (×Triticosecale Wittmack L.). In Vitro Cell Dev Biol Plant. 2014, 50: 143-148. 10.1007/s11627-013-9539-3.View ArticleGoogle Scholar
- Alheit KV, Reif JC, Maurer HP, Hahn V, Weissmann EA, Miedaner T, Würschum T: Detection of segregation distortion loci in triticale (x Triticosecale Wittmack) based on a high-density DArT marker consensus genetic linkage map. BMC Genomics. 2011, 12: 380-10.1186/1471-2164-12-380.PubMed CentralPubMedView ArticleGoogle Scholar
- Williams E, Piepho H-P, Whitaker D: Augmented p-rep designs. Biom J. 2011, 53: 19-27. 10.1002/bimj.201000102.PubMedView ArticleGoogle Scholar
- Cochran WG, Cox GM: Experimental Designs. 1957, New York: WileyGoogle Scholar
- Melchinger AE, Utz HF, Schön CC: Quantitative trait locus (QTL) mapping using different testers and independent population samples in maize reveals low power of QTL detection and larger bias in estimates of QTL effects. Genetics. 1998, 149: 383-403.PubMed CentralPubMedGoogle Scholar
- Gilmour AR, Gogel BG, Cullis BR, Thompson R: ASReml user Guide Release 3.0. VSN International Ltd, Hemel Hempstead. 2009, UK: HP1 1ESGoogle Scholar
- Blanc G, Charcosset A, Mangin B, Gallais A, Moreau L: Connected populations for detecting quantitative trait loci and testing for epistasis: An application in maize. Theor Appl Genet. 2006, 113: 206-224. 10.1007/s00122-006-0287-1.PubMedView ArticleGoogle Scholar
- SAS Institute: SAS/STAT 9.2 User's guide. 2008, Cary NCGoogle Scholar
- R Development Core Team: R: a language and environment for statistical computing. 2010, R foundation for statistical computing, http://www.R-project.org,Google Scholar
- Doerge RW, Churchill GA: Permutation tests for multiple loci affecting a quantitative character. Genetics. 1996, 142: 285-294.PubMed CentralPubMedGoogle Scholar
- Lander ES, Botstein S: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.PubMed CentralPubMedGoogle Scholar
- Utz HF, Melchinger AE, Schön CC: Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics. 2000, 154: 1839-1849.PubMed CentralPubMedGoogle Scholar
- Liu W, Maurer HP, Reif JC, Melchinger AE, Utz HF, Tucker MR, Ranc N, Della Porta G, Würschum T: Optimum design of family structure and allocation of resources in association mapping with lines from multiple crosses. Heredity. 2013, 110: 71-79. 10.1038/hdy.2012.63.PubMed CentralPubMedView ArticleGoogle Scholar
- Holland JB, Portyanko VA, Hoffmann DL, Lee M: Genomic regions controlling vernalization and photoperiod responses in oat. Theor Appl Genet. 2002, 105: 113-126. 10.1007/s00122-001-0845-5.PubMedView ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645. 10.1101/gr.092759.109.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to 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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.