Skip to main content

Genetic dynamics underlying phenotypic development of biomass yield in triticale



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 [1]. 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 [2]. 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. [3] 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. [4] 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. [7] 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 [4].

Triticale (× Triticosecale Wittmack L.) shows a broad variation for biomass yield [8] 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. [4] 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 [11]. Consistenly, the analysis of part of these data with an identity-by-descent approach identified substantially more QTL [12]. The aim of this study was therefore to reanalyse the triticale biomass data [4] 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.


In the entire mapping population with 647 DH lines derived from four families, we observed significant (P < 0.01) genotypic variances σ G 2 and genotype-by-environment interaction variances σ G × E 2 for biomass yield at all three developmental stages (BM1 – BM3) (Additional file 1: Table S1). The ratio of σ G 2 to σ G × E 2 ranged from approximately 3:1 for BM1 to 7:1 for BM2. The heritabilities were high and ranged from 0.81 for BM1 to 0.91 for BM2. Phenotypic correlations between biomass yield at the three developmental stages were 0.84 for BM1-BM2, 0.89 for BM2-BM3 and 0.75 for BM1-BM3 (all significant at P < 0.01). The phenotypic values of the parents differed to varying degrees ranging from Δ0.4 to 1.2 for BM1, from Δ0.0 to 1.7 for BM2 and from Δ0.1 to 1.9 for BM3 (Figure 1). Orthogonal contrasts of the means of the families and their respective parents were not significant except for family EAW78 for BM3 (P < 0.05). The trait distributions approximately followed a normal distribution except for family EAW78 for which the distribution at BM2 and BM3 became bimodal. The variation within each family increased with increasing developmental stage (Additional file 1: Table S1, Figure 1) and in each family DH lines transgressed the respective parents for all three developmental stages. Taken together, the data set captures the progression of biomass yield over time and is therefore well suited to study the underlying genetic dynamics.

Figure 1
figure 1

Phenotypic development of biomass yield. Histograms of the phenotypic values of biomass yield at three developmental stages (BM1-BM3) for the entire population and for each of the four families (DH06, DH07, EAW74, EAW78).

Multiple-line cross QTL mapping revealed 10 QTL for BM1 and BM2 and 9 QTL for BM3 (Table 1, Additional file 1: Table S2, Figure 2). Interestingly, we observed biomass QTL that were detected at all three developmental stages as well as QTL that were only detected at one or two of the developmental stages (Figure 3). Combined, the QTL explained between 62.2 (BM1) and 64.7 (BM2) percent of the genotypic variance. The proportion of genotypic variance explained by single QTL ranged from 2.5 to 28.0 for BM1, from 1.7 to 34.0 for BM2 and from 2.2 to 35.6 percent for BM3 (Additional file 1: Table S2). The major QTL detected at all three developmental stages was located on chromosome 5R. Another QTL which consistently contributed more than 5 percent to the genotypic variance was identified on chromosome 5A and on chromosome 6A a QTL contributing a similar proportion to BM2 and BM3 was found. We used fivefold cross-validation to assess the quality of the obtained QTL results (Table 1). This revealed a bias in the estimation of the proportion of genotypic variance which cross-validated still ranged between 34.0 (BM1) and 44.8 (BM2) percent. The QTL frequency distributions revealed some QTL that were detected in a high number of runs whereas other QTL identified with the full data set were only detected in few runs (Additional file 1: Figure S1). By contrast, few QTL were detected in a considerable number of runs which have not been identified in the full data set (e.g., BM2 QTL on chromosome 7A).We next assessed the contribution to the genotypic variance at BM1, BM2 and BM3 for all loci detected as QTL for any of the three developmental stages. This substantiated the presence of chromosomal regions which stably contribute to biomass yield as opposed to regions which only contribute to biomass yield in a temporally restricted manner (Figure 4).The full 2-dimensional epistasis scan revealed epistatic QTL for all three developmental stages (Figure 5a). Eight epistatic QTL were identified for BM1, three for BM2 and four for BM3. The contribution of these epistatic QTL to the genotypic variance was small and ranged between 1.0 and 3.5 percent. The epistatic interaction landscape undergoes temporal changes with development as illustrated in Figure 5b for the QTL detected at BM2 involving chromosomes 6A and 5B.

Table 1 Results of QTL mapping and fivefold cross-validation
Figure 2
figure 2

Main effect QTL. QTL for biomass yield detected at three developmental stages (BM1-BM3).

Figure 3
figure 3

Specific and overlapping QTL. Venn diagram for biomass yield QTL detected at three developmental stages (BM1-BM3).

Figure 4
figure 4

Temporal dynamics of QTL. Temporal development of the contribution of the QTL detected for any of the three developmental stages (BM1-BM3) to the proportion of explained genotypic variance for biomass yield.

Figure 5
figure 5

Epistatic QTL. (a) Epistatic QTL for biomass yield at three developmental stages (BM1-BM3). (b) Temporal development of the epistatic interaction landscape for the QTL detected at BM2 involving chromosomes 6A and 5B.


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 [13]. 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. [8] 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 [12]. 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. [4] 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 [14]. 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, 1517]. 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. [4] 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. [4]. In brief, phenotypic data for biomass yield were obtained by non-invasive prediction based on a precision phenotyping platform [7]. 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) [18]. 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. [21]. The population consisted of four families designated DH06 (131), DH07 (120), EAW74 (200), and EAW78 (196) which have been described by Alheit et al. [12]. The DH lines were grown in partially replicated designs [22] including common checks with 960 plots per location, at two locations in two years. Phenotypic data were analyzed by ordinary lattice analysis of variance [23]. Variance components were determined by the restricted maximum likelihood (REML) method assuming a full random model and heritability (h2) on an entry-mean basis was estimated from the variance components as the ratio of genotypic to phenotypic variance [24]. Best linear unbiased estimates (BLUEs) were estimated across environments assuming fixed effects for the genotype. All statistical analyses were performed using ASReml 3.0 [25].

Multiple-line cross QTL mapping

The DH lines were genotyped with DArT markers and QTL mapping was done based on the integrated consensus linkage map described by Alheit et al. [12]. For QTL mapping, an additive genetic model was chosen and the joint analysis across families was performed with a model assuming specific QTL effects for every family (disconnected model) [26] as described in detail by Steinhoff et al. [9]. In brief, the multiple-line cross QTL mapping model was:

Y = JM + X q B q + c q X c B c + ϵ

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 [27]. The presence of a putative QTL in an interval was tested using a likelihood-ratio test with the statistical software R [28]. 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 [29]. 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) [30]. The proportion of genotypic variance explained by the detected QTL was estimated as R2 adj /h2 [31]. 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. [32].

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. [33] 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 [34]. 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.


  1. Würschum T: Mapping QTL for agronomic traits in breeding populations. Theor Appl Genet. 2012, 125: 201-210. 10.1007/s00122-012-1887-6.

    Article  PubMed  Google Scholar 

  2. 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.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    CAS  PubMed Central  PubMed  Google Scholar 

  4. 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-

    Article  PubMed Central  PubMed  Google Scholar 

  5. 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.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. 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.

    Article  Google Scholar 

  8. 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.

    Article  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. 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.

    Article  PubMed  Google Scholar 

  11. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  12. 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.

    Article  PubMed  Google Scholar 

  13. 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.

    Article  PubMed  Google Scholar 

  14. Carlborg Ö, Haley CS: Epistasis: Too often neglected in complex trait studies?. Nat Rev Genet. 2004, 5: 618-625. 10.1038/nrg1407.

    Article  CAS  PubMed  Google Scholar 

  15. 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.

    Article  CAS  PubMed  Google Scholar 

  16. 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.

    Article  PubMed  Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. 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.

    Article  Google Scholar 

  19. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  20. 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.

    Article  Google Scholar 

  21. 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.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Williams E, Piepho H-P, Whitaker D: Augmented p-rep designs. Biom J. 2011, 53: 19-27. 10.1002/bimj.201000102.

    Article  PubMed  Google Scholar 

  23. Cochran WG, Cox GM: Experimental Designs. 1957, New York: Wiley

    Google Scholar 

  24. 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.

    CAS  PubMed Central  PubMed  Google Scholar 

  25. Gilmour AR, Gogel BG, Cullis BR, Thompson R: ASReml user Guide Release 3.0. VSN International Ltd, Hemel Hempstead. 2009, UK: HP1 1ES

    Google Scholar 

  26. 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.

    Article  CAS  PubMed  Google Scholar 

  27. SAS Institute: SAS/STAT 9.2 User's guide. 2008, Cary NC

    Google Scholar 

  28. R Development Core Team: R: a language and environment for statistical computing. 2010, R foundation for statistical computing,,

    Google Scholar 

  29. Doerge RW, Churchill GA: Permutation tests for multiple loci affecting a quantitative character. Genetics. 1996, 142: 285-294.

    CAS  PubMed Central  PubMed  Google Scholar 

  30. Lander ES, Botstein S: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.

    CAS  PubMed Central  PubMed  Google Scholar 

  31. 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 Central  PubMed  Google Scholar 

  32. 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.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. 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.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references


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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Tobias Würschum.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TW, HPM, VH, JCR, AR designed the study. HPM, VH, EAW collected phenotypes. WL, MG, TW carried out analyses. TW wrote the paper. All authors discussed the results and commented on the manuscript.

Electronic supplementary material


Additional file 1: Table S1: Summary statistics for biomass yield [t ha-1] at the three developmental stages. Table S2. QTL detected for biomass yield at three developmental stages (BM1-BM3). Figure S1. QTL frequency distributions. (PDF 50 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, W., Gowda, M., Reif, J.C. et al. Genetic dynamics underlying phenotypic development of biomass yield in triticale. BMC Genomics 15, 458 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: