Mapping dynamic QTL for plant height in triticale
- Tobias Würschum†1Email author,
- Wenxin Liu†2,
- Lucas Busemeyer3,
- Matthew R Tucker4,
- Jochen C Reif5,
- Elmar A Weissmann6,
- Volker Hahn1,
- Arno Ruckelshausen3 and
- Hans Peter Maurer1
© Würschum et al.; licensee BioMed Central Ltd. 2014
Received: 14 April 2014
Accepted: 8 May 2014
Published: 19 May 2014
Plant height is a prime example of a dynamic trait that changes constantly throughout adult development. In this study we utilised a large triticale mapping population, comprising 647 doubled haploid lines derived from 4 families, to phenotype for plant height by a precision phenotyping platform at multiple time points.
Using multiple-line cross QTL mapping we identified main effect and epistatic QTL for plant height for each of the time points. Interestingly, some QTL were detected at all time points whereas others were specific to particular developmental stages. Furthermore, the contribution of the QTL to the genotypic variance of plant height also varied with time as exemplified by a major QTL identified on chromosome 6A.
Taken together, our results in the small grain cereal triticale reveal the importance of considering temporal genetic patterns in the regulation of complex traits such as plant height.
Plant height in small grain cereals is an important agronomic trait affecting crop performance, particularly lodging and consequently grain yield and grain quality. The reduction of crop height has therefore been an important breeding goal for many decades . The identification of variants that reduce height without adversely affecting yield potential is still an important breeding target. Today, crop height in small grain cereals can be strongly modified by plant breeders through major dwarfing or semi-dwarfing genes. The introduction of the so-called Reduced height (Rht) genes into bread wheat was a key component of the ‘green revolution’. Alleles of Rht-B1 and Rht-D1 are nowadays the major sources of semi-dwarfism in wheat and actually increase grain yield in most environments [2, 3]. The Rht-1 loci encode DELLA proteins that integrate hormonal (gibberellin) and environmental signals to affect plant growth . Rht-1 is encoded by three homoeoloci located on group 4 chromosomes of the A, B and D genome of wheat . Two further important Rht genes, Rht7 and Rht8, were found to be located on chromosomes 2A and 2D, respectively, and may represent a homoeologous series . In rye, the major dwarfing gene Ddw1 is located on chromosome 5R and may represent a homeolog of Rht12[5–7]. Ddw1 has recently been shown to segregate in elite triticale germplasm and to exhibit a strong effect on plant height [8, 9]. In addition, the photoperiod insensitive alleles of the major photoperiod regulator Ppd-1, located on group 2 chromosomes, can also have pleiotropic effects on plant height .
The genetic control underlying plant height is usually studied by assessing the trait once the plants have reached their final height. Plant height, however, is a dynamic trait  that shows strong phenotypic changes during the adult plants’ development. Yan et al.  performed QTL mapping in a rice doubled haploid population phenotyped for plant height at multiple time points. This analysis revealed some QTL that were detected at all examined time points but also many QTL that could only be identified at one or some time points. Further support for the temporal changes of the genetic control underlying dynamic traits came from a recent study on biomass in triticale . This study revealed that the entire genetic architecture of biomass accumulation is under temporal dynamic control.
In this study we employed triticale (×Triticosecale Wittmack; 2n = 6x = 42; AABBRR) as a model species for small grain cereals and assessed plant height at three different time points by the precision phenotyping platform ‘BreedVision’ . A large doubled haploid population with 647 individuals derived from 4 families formed the basis for multiple-line cross QTL mapping  to identify plant height QTL at all three time points. In particular our objectives were to identify main and epistatic QTL, assess the contribution of these QTL to the genotypic variance of plant height, and to unravel the temporal genetic patterns underlying the phenotypic development of plant height.
Summary statistics for plant height (cm) at the three developmental stages (PH1-PH3)
Results of QTL mapping at three developmental stages (PH1-PH3) and fivefold cross-validation
Many traits of agronomic or biological importance undergo dynamic phenotypic changes during development. Plant height in cereals is a prime example for such a dynamic trait as it changes constantly throughout most of the plants’ adult life. Nevertheless, studies on the genetic control of plant height have so far mainly focussed on the final height, thus neglecting the trait’s developmental dynamics. In this study we therefore investigated plant height at three different time points to unravel the genetic dynamics underlying plant height development.
Phenotypic development of plant height
A prerequisite for dynamic QTL mapping is the assessment of the trait at multiple time points. In this study, we used the recently described ‘BreedVision’ precision phenotyping platform  to assess plant height at three different time points, approximately corresponding to the stages when awns become visible, late flowering and very early dough development . Based on an elaborate calibration experiment and employing triticale as model species, this precision phenotyping platform has been shown to accurately predict plant height . In this study, we observed high heritabilities above 0.9 for predicted plant height at all three time points, based on a large mapping population evaluated at 4 environments. These high heritabilities in combination with the high prediction accuracy demonstrate the great potential of precision phenotyping, especially for a temporal assessment of traits as required to dissect their underlying genetic dynamics.
The R 2 values between plant height at the three time points were highest for the contiguous time points and decreased slightly with increasing time between measurements (Additional file 1: Figure S1). This is comparable to the temporal development of another dynamic trait, biomass yield, in triticale . These results indicate that although the measurements at the different time points are not independent, there is variation that is unaccounted for. These differences between genotypes in phenotypic development thus indicate differences in the underlying genetics. The trait distributions approximately followed a normal distribution except for family EAW78 (Figure 1). This family has recently been shown to segregate for a major plant height QTL on chromosome 5R  which corresponds to the rye dwarfing gene Ddw1. This major QTL explained 40 - 50% of the genotypic variance at all three time points and is thus the most likely reason for the observed bimodal tendency in this family. Taken together, these results illustrate the phenotypic plasticity of plant height during development which underlines the need for a temporal assessment of dynamic traits.
The genetic architecture of plant height
The QTL with the highest contribution to the genotypic variance were identified on chromosomes 5R and 6A (Additional file 1: Table S1) and the QTL on chromosome 5R segregating in family EAW78 has recently been identified as Ddw1. This illustrates that the rye derived Ddw1 also exhibits a strong effect in triticale in the presence of the wheat A and B genomes. The QTL on chromosome 6A does not appear to represent one of the known major Rht loci and may be a novel QTL. No QTL was detected on chromosome 4B where Rht-B1 is located in wheat .
The QTL frequency distributions further supported the majority of the QTL as they were identified in a high number of runs. The cross-validated proportion of genotypic variance explained by the detected QTL amounted to a considerable 55.3 to 58.2% (Table 2). This can to some extent be explained by the major QTL on chromosome 5R but must also be attributed to other QTL spread over the entire genome. Our results thus indicate that plant height in triticale is a complex trait controlled by some major or medium effect QTL in addition to many small effect QTL.
Epistasis refers to interactions between the alleles at two or more genetic loci in the genome . Epistatic QTL contribute to the genetic architecture of many complex traits and have recently been reported for different crops including wheat, maize, and rapeseed [18–22]. The orthogonal contrasts between family means and the means of the respective parents can indicate the presence of epistasis. Non-significant orthogonal contrasts as observed here, however, do not exclude the presence of epistasis and the full 2-dimensional epistasis scans did indeed detect epistatic QTL for each of the three time points. The contribution of these epistatic interactions to the genotypic variance was rather small. It must be noted however that the power to detect epistatic QTL more strongly depends on population size than the power to detect main effect QTL and consequently many epistatic QTL may have remained undetected. Thus, while their individual contribution to the genotypic variance of plant height may be small, their combined contribution may be substantial assuming a higher number of epistatic interactions. Collectively, these results show that both main effect and epistatic QTL contribute to the genetic architecture of plant height in triticale.
Temporal genetic dynamics underlying plant height development
Dynamic traits are characterized by temporal phenotypic changes suggesting that the genetic control underlying such traits may also vary with time. QTL mapping for plant height in rice at multiple time points revealed many loci that were specific for one or a few of the time points . Employing triticale as model species for small grain cereals, Busemeyer et al.  and has recently reported dynamic QTL mapping for biomass and also reported developmental stage specific QTL.
In this study we assessed plant height of triticale at three different time points. Our results corroborate those from rice , as we also identified QTL that appear to contribute to the genetic control underlying plant height at all investigated time points, whereas other QTL were only identified at one or two of the time points. This illustrates that QTL underlying plant height are often not static and the trait not controlled by a fixed set of loci throughout crop height development. Rather, the genetic control of plant height undergoes rapid temporal changes. This is also illustrated by the major QTL on chromosome 6A. While this QTL was identified at all three time points, its contribution to the genotypic variance of plant height increased from 3.5% at PH1 to 10.7% at PH2 and finally 11.5% at PH3 (Figure 4). This suggests that the effects of this QTL become more pronounced after flowering of the plants. While at PH1 many QTL with small effect were identified, the least number of QTL was detected at PH3 for which the majority of the genotypic variance was contributed by the two major QTL on chromosomes 5R and 6A. This illustrates that assessing only the final plant height as represented by PH3 is not fully representative for the genetic architecture underlying plant height and would miss a number of QTL active during earlier development.
Interestingly, we found a similar dynamic genetic pattern for the epistatic interactions which also changed with development (Figure 5). Collectively, our results revealed that the entire genetic architecture underlying plant height shows dynamic temporal changes during crop development.
In this study, we employed a precision phenotyping platform to assess plant height of triticale at three time points. We show that both main effect and epistatic QTL are not static but rather are mainly detected at only one or two of the examined developmental stages. In addition, we observed variable contributions of QTL to the genotypic variance of plant height, as exemplified by the major QTL detected on chromosome 6A. Taken together, our results illustrate the temporal dynamics of the genetic control underlying plant height which emphasizes the need for multiple assessments of such dynamic traits.
Plant material, field trials and phenotypic data
The plant material and the field trials used in this study have been described in Busemeyer et al. . Phenotypic data for plant height were obtained by non-invasive prediction based on the ‘BreedVision’ precision phenotyping platform . Plant height was predicted at the three developmental stages: PH1 = BBCH stage 49 (awns visible), PH2 = BBCH 69 (late flowering), and PH3 = BBCH 81 (very early dough development) , in a mapping population consisting of 647 doubled haploid (DH)  triticale lines. The DH lines are derived from 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 alpha 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.7 for PH1, 4.6 for PH2 and 4.7 for PH3 were used corresponding to an experiment-wise type I error of P < 0.10, based on 2,000 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 . Plant height 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 previously [35, 36]. QTL frequency distributions were assessed by a 5 cM sliding window.
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 .
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 Kim Möller, Katharina Alheit, Agnes Rölfing-Finze, Hans Häge, Jacek Till and Justus von Kittlitz for outstanding work in the field and generation of the phenotypic data.
- Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Snape J: Meta-QTL analysis of the genetic control of crop height in elite European winter wheat germplasm. Mol Breed. 2012, 29: 159-171. 10.1007/s11032-010-9534-x.View ArticleGoogle Scholar
- Flintham JE, Börner A, Worland AJ, Gale MD: Optimizing wheat grain yield: effects of Rht (gibberellin-insensitive) dwarfing genes. J Agric Sci. 1997, 128: 11-25. 10.1017/S0021859696003942.View ArticleGoogle Scholar
- Wilhelm EP, Mackay IJ, Saville RJ, Korolev AV, Balfourier F, Greenland AJ, Boulton MI, Powell W: Haplotype dictionary for the Rht-1 loci in wheat. Theor Appl Genet. 2012, 126: 1733-1747.View ArticleGoogle Scholar
- Archard P, Cheng H, De Grauwe L, Decat J, Schoutteten H, Moritz T, Van Der Straeten D, Peng JR, Harberd NP: Integration of plant responses to environmentally activated phytohormonal signals. Science. 2006, 311: 91-94. 10.1126/science.1118642.View ArticleGoogle Scholar
- Börner A, Plaschke J, Korzun V, Worland AJ: The relationships between the dwarfing genes of wheat and rye. Euphytica. 1996, 89: 69-75. 10.1007/BF00015721.View ArticleGoogle Scholar
- Börner A, Korzun V, Voylokov AV, Weber WE: Detection of quantitative trait loci on chromosome 5R of rye (Secale cereale L.). Theor Appl Genet. 1999, 98: 1087-1090. 10.1007/s001220051171.View ArticleGoogle Scholar
- Korzun V, Melz G, Börner A: RFLP mapping of the dwarfing (Ddw1) and hairy peduncle (Hp) genes on chromosome 5 of rye (Secale cereale L.). Theor Appl Genet. 1996, 92: 1073-1077. 10.1007/BF00224051.PubMedView 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. 2014, 127: 251-260. 10.1007/s00122-013-2214-6.PubMedView ArticleGoogle Scholar
- Kalih R, Maurer HP, Hackauf B, Miedaner T: Effect of a rye dwarfing gene on plant height, heading stage, and Fusarium head blight in triticale (×Triticosecale Wittmack). Theor Appl Genet. 2014, In pressGoogle Scholar
- Börner A, Worland AJ, Plaschke J, Schumann E, Law CN: Pleiotropic effects of genes for reduced height (Rht) and day-length insensitivity (Ppd) on yield and its components for wheat grown in middle Europe. Plant Breed. 1993, 111: 204-216. 10.1111/j.1439-0523.1993.tb00631.x.View ArticleGoogle Scholar
- Wu R, Lin M: Functional mapping – how to map and study the genetic architecture of dynamic complex 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.PubMedPubMed CentralGoogle Scholar
- Busemeyer L, Ruckelshausen A, Möller K, Melchinger AE, Alheit KV, Maurer HP, 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-PubMedPubMed CentralView 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
- 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
- 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. Ann Appl Biol. 1991, 119: 561-601. 10.1111/j.1744-7348.1991.tb04895.x.View 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
- 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
- Steinhoff J, Liu W, Reif JC, Porta GD, 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, 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
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Alheit KV, Reif JC, Maurer HP, Hahn V, Weissmann EA, Miedaner T, Würschum T: Detection of segregation distortion loci in triticale (×Triticosecale Wittmack) based on a high-density DArT marker consensus genetic linkage map. BMC Genomics. 2011, 12: 380-10.1186/1471-2164-12-380.PubMedPubMed CentralView 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
- 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.PubMedPubMed CentralGoogle Scholar
- Gilmour AR, Gogel BG, Cullis BR, Thompson R: ASReml User Guide Release 3.0. 2009, Hemel Hempstead, HP1 1ES, UK: VSN International LtdGoogle 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
- Steinhoff J, Liu W, Maurer HP, Würschum T, Longin CFH, 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
- SAS Institute: SAS/STAT 9.2 User's Guide. 2008, Cary, NC, USAGoogle 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.PubMedPubMed CentralGoogle Scholar
- Lander ES, Botstein S: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.PubMedPubMed CentralGoogle 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.PubMedPubMed CentralGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Würschum T, Kraft T: Cross-validation in association mapping and its relevance for the estimation of QTL parameters of complex traits. Heredity. 2014, 112: 463-468. 10.1038/hdy.2013.126.PubMedPubMed CentralView 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.PubMedPubMed CentralView 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/4.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.