Next Article in Journal
Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem
Next Article in Special Issue
Quantitative Trait Loci and Candidate Genes Associated with Cold-Acclimation and Microdochium nivale Tolerance/Susceptibility in Winter Triticale (x Triticosecale)
Previous Article in Journal
Quinoa for the Brazilian Cerrado: Agronomic Characteristics of Elite Genotypes under Different Water Regimes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Fine-Mapping of Quantitative Trait Loci Controlling Plant Height in Central European Winter Triticale (×Triticosecale Wittmack)

1
State Plant Breeding Institute, University of Hohenheim, 70599 Stuttgart, Germany
2
Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Plants 2021, 10(8), 1592; https://doi.org/10.3390/plants10081592
Submission received: 5 July 2021 / Revised: 26 July 2021 / Accepted: 29 July 2021 / Published: 2 August 2021
(This article belongs to the Special Issue Genetics, Profiling and Breeding of Triticale)

Abstract

:
The quantitatively inherited trait plant height is routinely evaluated in triticale breeding programs as it substantially influences lodging and disease susceptibility, is a main contributor to biomass yield, and is required to improve hybrid seed production by fine-tuning plant height in the female and male parental pools in hybrid breeding programs. In this study, we evaluated a panel of 846 diverse Central European triticale genotypes to dissect the genetic architecture underlying plant height by genome-wide association mapping. This revealed three medium- to large-effect QTL on chromosomes 5A, 4B, and 5R. Genetic and physical fine-mapping of the putative QTL revealed that the QTL on chromosome 5R most likely corresponds to Ddw1 and that the QTL on chromosome 5A is likely to be Rht12. Furthermore, we observed a temporal trend in registered cultivars with a decreasing plant height during the past decades, accompanied by an increasing use of the height-reducing alleles at the identified QTL. In summary, our results shed new light on the genetic control of plant height in triticale and open new avenues for future improvement by breeding.

1. Introduction

Hexaploid triticale (×Triticosecale Wittmack) is a man-made cereal that combines the genomes of tetraploid durum wheat (Triticum durum L.) and diploid rye (Secale cereal L.) [1]. The first triticale cultivars were released in the 1970s [1] and since then the triticale acreage has grown to approximately four million hectares worldwide as of 2019 [2]. Today’s grain use is mainly restricted to animal fodder and bioethanol production [1,3,4]. Its biomass is used as fresh fodder [4] and more recently for biogas production due to a promotional policy in the European Union [5,6,7,8].
Knowledge about the genetic control of plant height is of interest in small-grain cereal breeding programs. Tall plants usually are more susceptible to lodging resulting in substantial grain yield losses, reduced quality, and high drying costs during the harvesting process [9,10,11]. They also have a less favorable harvest index and the tendency consequently was to breed for shorter types. However, taller genotypes are also associated with a reduced Fusarium head blight susceptibility [12,13,14] and generally show an increased biomass yield, for which plant height is an important component of biomass yield [15,16,17]. Biomass yield has gained increasing interest as a breeding goal in triticale breeding programs in recent years due to the high potential of triticale as a bioenergy crop [8,18]. In addition, hybrid breeding has shown potential in triticale, which requires the control of the plant height of the male and female parental components, as the male lines should be taller compared to their female mating partners in order to increase the efficiency of hybrid seed production [19]. Thus, adjustment of plant height is required in triticale breeding programs, but the direction depends on the breeding goals.
Plant height in small-grain cereals is considered a quantitatively inherited trait controlled by many minor as well as by a few major quantitative trait loci (QTL) [17,20,21,22]. However, some height-reducing loci, such as the Rht-1 homoeoloci, show adverse pleiotropic effects on other developmental processes, such as root elongation, early seedling vigor, or coleoptile length, which has been extensively studied in wheat [23,24,25,26]. In wheat and rye, the founding species of triticale, a large number of height-reducing loci are known [21,27,28,29,30]. Of the major QTL located on the A or B genome that are used in wheat breeding and are therefore potentially suitable to be used in triticale breeding, Rht-B1 located on chromosome 4B has previously been reported in spring triticale [31]. QTL mapping in bi-parental triticale families has revealed several QTL on different chromosomes, but most of them with only small effects [14,17,32]. The only major QTL identified in triticale so far is Ddw1, located on chromosome 5 of the rye genome [14,17].
Consequently, the aims of this study were to investigate the phenotypic variation, the genetic architecture, and long-term genetic trends of plant height in Central European winter triticale. In particular, our objectives were to (i) perform genome-wide association mapping in a diversity panel of 846 Central European triticale genotypes, (ii) fine-map the genomic regions showing significant associations with plant height and determine their co-location with known height-reducing loci from wheat and rye, and (iii) evaluate long-term genetic trends and the allele frequencies of the detected QTL in registered cultivars.

2. Results

This study was based on 846 diverse triticale lines comprising 129 registered cultivars and 717 advanced breeding lines. All genotypes were evaluated in multi-location field trials for plant height and in addition for their developmental stage (BBCH scale) at a time point when the majority of the lines were heading, as this trait is often associated with plant height. We observed significant variation for both traits and high heritabilities of 0.81 for plant height and 0.80 for the developmental stage (Table 1). Plant height ranged from 87.3 cm to 126.6 cm and the BBCH stage ranged from 46.3 (opening of the flag leaf sheath) to 59.8 (completion of ear emergence). The correlation between the two traits was 0.24 (p < 0.001).
We used the 129 registered cultivars to investigate long-term trends resulting from breeding in the period from 1982 when the first cultivars of this panel were released until today. Plant height declined in the period from 1982 to 2010 from on average 110.4 cm to 102.4 cm and showed an increasing proportion of taller genotypes from 2011 on (Figure 1). By contrast, the developmental stage did not change over time, except that there may be a slight tendency towards earlier heading in the more recent cultivars, i.e., a higher and thus more advanced BBCH stage.
We performed genome-wide association mapping using 31,823 markers to investigate the genetic architecture underlying the traits plant height and developmental stage (Figure 2). We identified 27 significantly associated markers for plant height and 17 for developmental stage (Supplementary Table S1). Jointly, these markers accounted for 42.16% of the total genotypic variance of plant height and 29.31% of that of developmental stage (Supplementary Table S1). The strongest association was found for the markers on chromosome 5R (Table 2). In a single fit, the QTL on 5A, 4B, and 5R explained more than 10% of the genotypic variance. In a joint fit, the QTL on chromosome 5R was fitted first and retained its high proportion of explained genotypic variance for plant height, whereas that of the 5A and 4B QTL was reduced to 0.40 and 2.34%, respectively. On chromosome 5R, three markers explained more than 1% of the genotypic variance even in a joint fit. This might indicate three separate QTL or that none of the markers alone is able to capture the full variance explained by a single QTL in this region. The fact that all three markers are located in close proximity, at 990.4, 997.0, and 1019.4 cm, suggests that they all identify the same QTL. Thus, the significantly associated markers might correspond to four QTL for plant height, located on chromosomes 5A, 4B, 4R, and 5R. The absolute values of the QTL allele substitution effects ranged between 1.03 and 5.30 cm (Table 2). Notably, the QTL effects varied substantially with the application rate of growth regulators and approximately doubled when only one instead of two rates were applied (Supplementary Figure S1). For the QTL on chromosome 5R, for example, the difference between the two allelic classes was 9.51 cm at the location Hohenheim when growth regulators were applied twice, which increased to 19.19 cm in the observation plots with only one application. For the developmental stage, three QTL were identified, the same QTL on chromosomes 5A, 4B, and 5R were also identified for plant height (Table 2). For these QTL, the allele that increased plant height also advanced the developmental stage, i.e., resulted in an earlier heading.
Next, we genetically and physically fine-mapped the plant height QTL on chromosomes 5A, 4B, and 5R (Figure 3a,b). The marker sequences were BLASTed against the respective reference genomes of wheat and rye to obtain their physical positions and to compare the QTL regions with those of known height-reducing loci. The significantly associated markers of the chromosome 5A QTL are located between 692 and 699 Mbp. For chromosome 4B, the significantly associated markers are located in the region between 657 and 666 Mbp, whereas Rht-B1 is located at the front end of this chromosome. Physical map positions of all significant markers on chromosome 5R are located within the region from 842 to 875 Mbp. The dominant plant height locus Ddw1 has recently been reported to be located at 862 Mbp [33]. Analysis of linkage disequilibrium (LD) among the significantly associated markers of these three QTL regions revealed a possible LD between some markers in the 4B and 5R regions (Figure 3c).
To better understand the effect of the detected QTL on plant height and their utilization in triticale breeding, we used the marker with the strongest association to represent the QTL on chromosomes 5A, 4B, and 5R (Figure 4, significantly associated markers and their allele status for the 129 registered cultivars can be found in the Supplementary Tables S2–S5). We found that in the 129 registered cultivars all three QTL occurred with both alleles and the difference between the two homozygous classes ranged between 3.5 and 6.7 cm. Notably, however, our analysis revealed that the height-reducing alleles often occur in combination (Figure 5). Twenty-four of the cultivars carry none of the height-reducing alleles at these three QTL, 97 carry the height-reducing allele at the 5A QTL, 59 at the 4B QTL, and 20 at the 5R QTL. The 20 cultivars carrying the height-reducing allele at the 5R QTL all also carry the height-reducing alleles at the 5A and 4B QTL. Likewise, the 59 cultivars carrying the height-reducing 4B QTL allele also carry the height-reducing allele at the 5A QTL, except for two. Besides the latter two cultivars, this leaves 40 that only carry the height-reducing allele at the 5A QTL, so at only one of the three QTL. The analysis of the advanced breeding lines and the complete panel showed similar results (Supplementary Figures S2–S4). The only difference compared to the subset of registered cultivars was the higher frequency in the whole panel of genotypes carrying all three QTL and the lower frequency of genotypes carrying only the QTL on chromosome 5A (Figure 5, Supplementary Figure S4).
Last, we investigated the development of the frequencies of the height-reducing alleles at these three QTL in registered cultivars over time. This revealed an increase for all three QTL, as well as an increasing frequency of their combinations (Figure 6). For the 5A QTL, the frequency of the height-reducing allele was already around 0.5 in the cultivars registered before 1991, which then increased to above 0.9 in the cultivars registered in the past decade. For the 4B QTL, the same frequency started at around 0.2 and then increased to around 0.6, and for the QTL on chromosome 5R the frequency of the height-reducing allele was low in the cultivars registered until the year 2000 and only then increased to around 0.5 in the cultivars registered since 2011.

3. Discussion

Plant height is of high importance in triticale breeding programs and in this study, we evaluated a diverse panel of 846 Central European triticale lines to investigate the genetic architecture underlying plant height, fine-map, and characterize identified QTL and draw conclusions for triticale breeding.

3.1. Genome-Wide Association Mapping and Characterization of QTL for Plant Height

Genome-wide association mapping revealed four putative QTL for plant height on chromosomes 5A, 4B, 4R, and 5R, of which the 5A, 4B, and 5R QTL can be considered as medium- to large-effect QTL. The QTL on chromosomes 5A, 4B, and 5R were also identified for developmental stage and thus have pleiotropic effects on both traits. The effects of these QTL on plant height may appear smaller than expected, but were shown to strongly depend on the application rate of growth regulators (Supplementary Figure S1). Notably, owing to the large variation present in the diversity panel, growth regulators were applied in these trials to prevent lodging and to mirror agronomic practice. Thus, the QTL effects can be expected to be larger in non-treated plants and reflect the substantial effect of the QTL even in combination with growth regulators.
In a single fit model with only one QTL, the 5A and 4B QTL explained 13.06 and 18.97% of the genotypic variance, respectively, which reduced to 0.40 and 2.34% in a joint fit of all QTL. Notably, the 5R QTL is fitted first in the joint fit model and thus retains its explained variance. The finding that the other two QTL do not, can have different reasons. It might be due to linkage disequilibrium among the QTL regions. However, at least for the 5A QTL, this LD with the 5R QTL was rather low and therefore LD among these QTL appears as a less likely explanation. A more likely reason is the joint occurrence of the height-reducing alleles of the three QTL in triticale lines observed in the registered cultivars but also in the entire panel (Figure 5). As the height-reducing allele of the 5R QTL occurs together with those of the 5A and 4B QTL, the genotypic classes are always confounded. Regarding the QTL on chromosomes 4B and 5R, the available genetic map positions as well as physical positions linked the significantly associated markers to regions on these two chromosomes. In an unrelated bi-parental population, however, we observed tight linkage among all these markers that segregated in that population (data not shown), indicating that they may actually identify only one QTL. Thus, further research using bi-parental families is required to disentangle the identities and effects of these QTL.
For the major QTL on chromosome 5R, three markers remained which in a joint fit still explained 1.79, 1.91, and 29.38% of the genotypic variance. While this could mean that there are three separate QTL, we concluded that it is more likely that they all identify the same QTL, as they are located in close proximity. Probably even the most strongly associated marker is not in perfect LD with the causal variant, such that also other markers can capture a small part of the variance explained by the QTL. Korzun et al. [34] mapped the dominant dwarfing gene Ddw1 on the long arm of chromosome 5R. Braun et al. [33] recently identified a gibberellin 2-oxidase gene (ScGA2ox12) as a candidate for Ddw1 in rye, which is located at 862 Mbp on chromosome 5R. As Ddw1 has been reported in spring and winter triticale [14,17,31,35] and as the QTL region coincides with the position of Ddw1, we conclude that the 5R QTL corresponds to Ddw1. Ddw1 is present in 20 of the 129 registered cultivars analyzed in this study and appears to have been increasingly used in recent years. However, previous work also illustrated adverse effects of Ddw1, for example on Fusarium head blight resistance [14]. The utilization of this height-reducing locus in triticale, therefore, requires further research on its advantages and limitations.
The QTL on chromosome 5A was located in the region between 692 and 699 Mbp (Figure 3b) and recently the dominant height-reducing Rht12 locus has been reported to be located in the region between 700 and 710 Mbp of chromosome 5A in wheat [36]. It thus appears likely that the QTL on chromosome 5A corresponds to Rht12. To our knowledge, Rht12 has not been reported in triticale germplasm so far, but its effect on plant height has been extensively investigated in wheat [23,37,38]. The height-reducing Rht12 is reported to delay ear emergence in wheat [37,38] but this also depends on photoperiodic sensitivity loci [37]. We observed a pleiotropic effect of the QTL on chromosome 5A, with the allele that reduces plant height delaying development and resulting in later heading. Thus, further research is necessary to investigate the effects of the putative Rht12 locus in triticale on other agronomic traits such as grain yield, ear emergence, disease resistance, early seedling vigor, or coleoptile length.
The height-reducing Rht-B1 locus located on chromosome 4B has been investigated extensively in wheat [12,20,39,40,41], but much less so in triticale [31]. In this study, we identified a QTL on chromosome 4B with pleiotropic effects on plant height and developmental stage (Figure 2; Supplementary Table S1). However, Rht-B1 is located on the short arm of chromosome 4B at approximately 30.8 Mbp, whereas the 4B QTL identified here is located at the opposite end of this chromosome (Figure 3b). Thus, Rht-B1 is no candidate for this QTL, which is in line with previous findings in triticale, where Rht-B1 did not contribute to the genetic control of plant height [32]. This QTL also showed a substantial plant height reduction and a high penetrance in the registered cultivars used in this study and is therefore of interest for future research.

3.2. Long-Term Trends of Plant Height in Triticale

The analysis of the development of plant height from the early 1980s—when the first cultivars of this panel were registered—until today, revealed a constant reduction in plant height (Figure 1) that is in line with previous findings [9]. In recent years, however, the demand for triticale genotypes with a greater biomass yield increased [15,42]. As plant height is one of the main contributors to biomass yield [15,16,17], taller genotypes were favored in breeding programs focusing on biomass yield.
In line with the phenotypic observation, we observed an increasing frequency of the height-reducing alleles at the three investigated QTL (Figure 6). The QTL on chromosome 5A, assumed to be Rht12, was constantly present at a medium to high frequency throughout the evaluated period. By contrast, especially the frequency of the QTL on chromosome 5R, likely Ddw1, increased only more recently since the turn of the millennium.

3.3. Conclusions

In this study, we identified three medium- to large-effect QTL and fine-mapping indicated that two of them correspond to Rht12 and Ddw1. While Ddw1 has been reported in triticale previously [14,17,31,35], Rht12 has not yet been reported in triticale.
The widely used height-reducing locus Rht-B1 has been one of the main contributors to improving grain yield and yield stability in wheat [20,39,43], but can negatively affect grain yield under drought conditions [43,44,45]. With the consequences of climate change [46], breeders will have to improve drought resistance also in triticale breeding material [47]. In wheat, the effect of height-reducing loci on seedling establishment and associated traits such as coleoptile length and early seedling vigor has been widely evaluated under dry conditions [25,26,48]. Genotypes carrying the dwarfing genes Rht8 and Rht12 have been reported to react less pronounced to heat and drought stress [23,37,49]. The effect of height-reducing loci on drought tolerance in triticale lacks behind what is known already in wheat and further research is required to evaluate the putative Rht12 as well as other loci for their effects on drought tolerance and thus their potential for future breeding of climate-resilient triticale.
In conclusion, the diversity panel showed significant genotypic variation for plant height, which can be exploited to tailor plant height in breeding programs, either by phenotypic or marker-assisted selection. In addition, height-reducing loci from wheat, such as Rht-B1 or Rht24 [20,39], do not appear to be exploited in Central European triticale so far and could be introgressed from wheat.

4. Materials and Methods

4.1. Phenotypic Data

This study was based on a total of 846 triticale (×Triticosecale Wittmack) genotypes consisting of a diverse collection of officially registered Central European cultivars (n = 129) and lines in advanced breeding status (n = 717, for further details see Neuweiler et al. [50]). Phenotypic data of the diversity panel were taken from a large field trial (n = 1280), for which genotypes were grown in two separate sets comprising 800 (set A) and 500 (set B) individuals including 20 common checks. All genotypes were planted in yield plots (Y) with a plot size ranging from 5 to 10.5 m² and in addition, in double row observation plots (O) in the growing seasons 2014 and 2015. All field trials were designed as partially replicated α-lattice designs [51]. The yield plot trials had an average replication number of 1.3 and 1.2 for set A and B, respectively, observation plots were replicated twice. The field trials were grown in Eckartsweier (EWE, Y, 48°31′18′′ N, 7°52′18′′ E, 140 m above sea level, masl), Franconia (FRA, Y, 2014: 49°39′58′′ N, 9°47′30′′ E, 310 masl, 2015: 49°49′22′′ N, 10°6′19′′ E, 270 masl), Hohenheim (HOH, O, Y, 48°28′49′′ N, 9°11′16′′ E, 400 masl), Ihinger Hof (IHO, O, Y, 48°44′40′′ N, 8°55′25′′ E, 480 masl), Moosburg (MSB, Y, 48°26′36′′ N, 11°45′22′′ E, 420 masl), and Oberer Lindenhof (OLI, O, 48°28′49′′ N, 9°18′56′′ E, 700 masl) in 2014 and 2015 as follows: in 2014 only set A was evaluated as yield plots at the locations EWE, FRA, HOH, IHO, and MSB. In 2015, sets A and B were grown at HOH, IHO, and OLI as observation plots and in EWE, FRA, HOH, IHO, and MSB as yield plots. Genotypes grown in yield plots were treated with growth regulators twice and in observation plots once to prevent lodging. Plant height (cm) was assessed after flowering from the ground to the tip of the ears, excluding awns. In addition, the developmental stage was scored on a BBCH scale according to Zadoks et al. [52] when the ears of the majority of the genotypes were emerging.

4.2. Statistical Analysis

The linear mixed models used in this study followed the syntax as outlined by Piepho et al. [53], where dot operators specify crossed effects and fixed and random effects are separated by a colon, introducing fixed effects first. We used the following model to obtain best linear unbiased estimates (BLUEs) and least significant differences (LSD) as well as a full random model to determine variance components and heritability estimates:
G:E + E·G + E·T + E·T·R·B + E·T·R·B
where E, G, T, R, and B denotes environments, genotypes, trial effects as a combination between genotype sets A and B and plot types (yield plots, double row observation plots), replications within each environment-trial combination, and incomplete blocks within replications of each environment-trial combination. Trial main effects were nested within environments. Variance components for the diversity panel used in this study were calculated using dummy variables in the above model. Variance components and heritabilities were estimated with the full random model using the restricted maximum likelihood method implemented in the software package ASReml-R 3.0 [54]. We assumed heterogeneous error variances at an environment-trial level. Variance components were tested for significance (p < 0.05, 0.01, 0.001) using a likelihood ratio test [55]. LSDs (p < 0.05) were calculated as an approximation using the twofold of the average standard error or a difference. Broad-sense heritability (H²) was estimated according to Cullis et al. [56] using the mean variance of a difference ( v ¯ B L U P ) between two best linear unbiased predictors and the genotypic variance component ( σ G 2 ) as:
H 2 = 1 v ¯ B L U P 2     σ G 2

4.3. Molecular Data Analysis

Marker data were available for all of the 846 genotypes of the diversity panel and obtained by a genotyping-by-sequencing approach from Diversity Arrays Technology, Canberra, Australia (www.diversityarrays.com, accessed on 14 May 2021). Markers showing more than 20% missing values and a minor allele frequency of 5% or less were discarded. Unmapped significant markers were assigned to their most likely chromosomal position by their linkage disequilibrium with mapped markers. These measures yielded 31,823 markers in total that were used for the genome-wide association study. Positions for 10,192 markers on the A genome and 14,747 markers on the B genome were known [57]. Map positions for markers from the R genome were assigned as described in Neuweiler et al. [50] resulting in 6884 markers with a known map position. In total, 25,721 dominant silico-DArT and 6102 SNP markers were used. To differentiate between the silico-DArT and SNP markers, we assigned a “D” and an “S” prefix for silico-DArT and SNP markers, respectively. To determine the physical positions of markers in the QTL regions, the sequences of the markers were BLASTed against the reference genomes of wheat and rye (IWGSC RefSeq v1.0 and Secale cereale Lo7).

4.4. Genome-Wide Association Study

We performed genome-wide association mapping using a mixed linear model incorporating a kinship matrix [58] implemented in the GAPIT R package [59]. We used a Bonferroni-corrected significance threshold of p < 0.05. The genotypic data, filtered as stated above, were imputed using the k-nearest neighbor genotype imputation technique LD-kNNi implemented in LinkImpute [60]. Furthermore, we calculated the total proportion of the genotypic variance (pG) explained by all detected QTL as pG = R2adj/H2, where R2adj was the adjusted R2 from the linear model, and H2 as the heritability of the trait [61]. To correct for collinearity, we ordered the QTL in the order of their strength of association and calculated their individual proportion of explained genotypic variance accordingly. The pG values for individual QTL were derived by estimating their sums of squares in the linear model. The allele substitution (α) effect was derived as α = a(1 + k(p1p2)) with a as genotypic value of the corresponding locus, k the degree of dominance, as well as p1 and p2 as the allele frequencies. The α-effect corresponds to one-half of the difference of the genotypic values of the corresponding genotypic classes of a QTL when inbred lines are considered. The α-effects were acquired by fitting a linear model for the marker of interest and the corresponding trait. Then the α-effect is represented by the determined regression coefficient [62].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10081592/s1, Figure S1: Effects of the QTL qPH.5A, qPH.4B, and qPH.5R.C shown for observation and yield plots at Hohen-heim and Ihinger Hof. Observation plots were treated once, yield plots twice with growth regulators. The height-reducing QTL allele is shown in blue. Figure S2: Effects of the QTL detected on chromosomes 5A, 4B, and 5R as well as their combinations on plant height, assessed in the 846 registered cultivars and advanced breeding lines. ‘+’ and ‘−’ indicate presence or absence of the height-reducing QTL alleles. If no presence/absence labeling is shown, the allelic state of the respective QTL was not considered. Figure S3: Effects of the QTL detected on chromosomes 5A, 4B, and 5R as well as their combinations on plant height, assessed in the 717 advanced breeding lines. ‘+’ and ‘−’ indicate presence or absence of the height-reducing QTL alleles. If no presence/absence labeling is shown, the allelic state of the respective QTL was not considered. Figure S4: Number of genotypes carrying one or a combination of plant height reducing alleles at the detected QTL for the whole population including registered cultivars and advanced breeding lines (N = 846, left) as well as for the ad-vanced breeding lines alone (N = 717, right). Table S1: Significant markers detected for plant height and development stage by the genome wide association study.

Author Contributions

H.P.M. designed the experiments. J.T. and J.E.N. analyzed the data. J.T., H.P.M. and T.W. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Federal Ministry of Food and Agriculture (BMEL) through its project management body Fachagentur für Nachwachsende Rohstoffe e.V. (FNR) within the grant numbers 22024515, 22406112, 22406212, 22406312, and 22406412.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Mergoum, M.; Singh, P.K.; Peña, R.J.; Lozano-del Río, A.J.; Cooper, K.V.; Salmon, D.F.; Gómez Macpherson, H. Triticale: A “New” Crop with Old Challenges. In Cereals; Carena, M.J., Ed.; Springer US: New York, NY, USA, 2009; pp. 267–287. ISBN 978-0-387-72294-8. [Google Scholar]
  2. FAOSTAT. Statistical Databases and Datasets of the Food and Agriculture Organization of the United Nations. Available online: http://www.fao.org/faostat (accessed on 4 May 2021).
  3. McGoverin, C.M.; Snyders, F.; Muller, N.; Botes, W.; Fox, G.; Manley, M. A review of triticale uses and the effect of growth environment on grain quality. J. Sci. Food Agric. 2011, 91, 1155–1165. [Google Scholar] [CrossRef] [PubMed]
  4. Myer, A.; Del Lozano Río, A.J. Triticale as animal feed. In Triticale Improvement and Production; Mergoum, M., Gómez Macpherson, H., Eds.; Food and Agriculture Organization Of The United Nations: Rome, Italy, 2004; pp. 49–58. [Google Scholar]
  5. EU. Directive (EU) 2018/2001 of the European Parliament and of the Council of 11 December 2018 on the Promotion of the Use of Energy from Renewable Sources (recast). Off. J. Eur. Union 2018, 61, 82–209. [Google Scholar]
  6. EU. Directive 2009/28/EC of the European Parliament and of the Council of 23 April 2009 on the Promotion of the Use of Energy from Renewable Sources and Amending and Subsequently Repealing Directives 2001/77/EC and 2003/30/EC. Off. J. Eur. Union 2009, 52, 16–62. [Google Scholar]
  7. EU. Directive 2001/77/EC of the European Parliament and of the Council of 27 September 2001 on the Promotion of Electricity Produced from Renewable Energy Sources in the Internal Electricity Market. Off. J. Eur. Union 2001, 44, 33–40. [Google Scholar]
  8. Weiland, P. Biogas production: Current state and perspectives. Appl. Microbiol. Biotechnol. 2010, 85, 849–860. [Google Scholar] [CrossRef]
  9. Losert, D.; Maurer, H.P.; Marulanda, J.J.; Würschum, T. Phenotypic and genotypic analyses of diversity and breeding progress in European triticale (×Triticosecale Wittmack). Plant Breed. 2017, 136, 18–27. [Google Scholar] [CrossRef]
  10. Liu, W.; Leiser, W.L.; Maurer, H.P.; Li, J.; Weissmann, S.; Hahn, V.; Würschum, T. Evaluation of genomic approaches for marker-based improvement of lodging tolerance in triticale. Plant Breed. 2015, 134, 416–422. [Google Scholar] [CrossRef]
  11. Rajkumara, S. Lodging in cereals—A review. Agric. Rev. 2008, 29, 55–60. [Google Scholar]
  12. Miedaner, T.; Voss, H.-H. Effect of Dwarfing Rht Genes on Fusarium Head Blight Resistance in Two Sets of Near-Isogenic Lines of Wheat and Check Cultivars. Crop. Sci. 2008, 48, 2115–2122. [Google Scholar] [CrossRef]
  13. Yan, W.; Li, H.B.; Cai, S.B.; Ma, H.X.; Rebetzke, G.J.; Liu, C.J. Effects of plant height on type I and type II resistance to fusarium head blight in wheat. Plant Pathol. 2011, 60, 506–512. [Google Scholar] [CrossRef]
  14. Kalih, R.; Maurer, H.P.; 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, 127, 1527–1536. [Google Scholar] [CrossRef]
  15. Gowda, M.; Hahn, V.; Reif, J.C.; Longin, C.F.H.; Alheit, K.; Maurer, H.P. Potential for simultaneous improvement of grain and biomass yield in Central European winter triticale germplasm. Field Crop. Res. 2011, 121, 153–157. [Google Scholar] [CrossRef]
  16. Losert, D.; Maurer, H.P.; Weissmann, S.; Würschum, T. Hybrid breeding for biomass yield in winter triticale: I. Hybrid performance, trait correlations and heterosis. Plant Breed. 2016, 135, 560–566. [Google Scholar] [CrossRef]
  17. Alheit, K.V.; Busemeyer, L.; Liu, W.; Maurer, H.P.; Gowda, M.; Hahn, V.; Weissmann, S.; Ruckelshausen, A.; Reif, J.C.; Würschum, T. Multiple-line cross QTL mapping for biomass yield and plant height in triticale (×Triticosecale Wittmack). Theor. Appl. Genet. 2014, 127, 251–260. [Google Scholar] [CrossRef]
  18. LTZ. Energiepflanzen für Biogasanlagen. Available online: https://pflanzen.fnr.de/ (accessed on 4 April 2021).
  19. Longin, C.F.H.; Mühleisen, J.; Maurer, H.P.; Zhang, H.; Gowda, M.; Reif, J.C. Hybrid breeding in autogamous cereals. Theor. Appl. Genet. 2012, 125, 1087–1096. [Google Scholar] [CrossRef] [PubMed]
  20. Würschum, T.; Langer, S.M.; Longin, C.F.H. Genetic control of plant height in European winter wheat cultivars. Theor. Appl. Genet. 2015, 128, 865–874. [Google Scholar] [CrossRef]
  21. Börner, A.; Plaschke, J.; Korzun, V.; Worland, A.J. The relationships between the dwarfing genes of wheat and rye. Euphytica 1996, 89, 69–75. [Google Scholar] [CrossRef]
  22. 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. [Google Scholar] [CrossRef]
  23. Ellis, M.H.; Rebetzke, G.J.; Chandler, P.; Bonnett, D.; Spielmeyer, W.; Richards, R.A. The effect of different height reducing genes on the early growth of wheat. Funct. Plant Biol. 2004, 31, 583–589. [Google Scholar] [CrossRef]
  24. Botwright, T.; Rebetzke, G.; Condon, T.; Richards, R. The effect of rht genotype and temperature on coleoptile growth and dry matter partitioning in young wheat seedlings. Funct. Plant Biol. 2001, 28, 417–423. [Google Scholar] [CrossRef]
  25. Rebetzke, G.J.; Appels, R.; Morrison, A.D.; Richards, R.A.; McDonald, G.; Ellis, M.H.; Spielmeyer, W.; Bonnett, D.G. Quantitative trait loci on chromosome 4B for coleoptile length and early vigour in wheat (Triticum aestivum L.). Aust. J. Agric. Res. 2001, 52, 1221–1224. [Google Scholar] [CrossRef]
  26. Bai, C.; Liang, Y.; Hawkesford, M.J. Identification of QTLs associated with seedling root traits and their correlation with plant height in wheat. J. Exp. Bot. 2013, 64, 1745–1753. [Google Scholar] [CrossRef] [Green Version]
  27. Kantarek, Z.; Masojć, P.; Bienias, A.; Milczarski, P. Identification of a novel, dominant dwarfing gene (Ddw4) and its effect on morphological traits of rye. PLoS ONE 2018, 13, e0199335. [Google Scholar] [CrossRef] [PubMed]
  28. Stojałowski, S.; Myśków, B.; Hanek, M. Phenotypic effect and chromosomal localization of Ddw3, the dominant dwarfing gene in rye (Secale cereale L.). Euphytica 2015, 201, 43–52. [Google Scholar] [CrossRef] [Green Version]
  29. McIntosh, R.A.; Dubcovsky, J.; Rogers, W.J.; Morris, C.; Xia, X.C. Catalogue of Gene Symbols for Wheat: 2017 Supplement. Available online: https://shigen.nig.ac.jp/wheat/komugi/genes/symbolClassList.jsp (accessed on 14 May 2021).
  30. McIntosh, R.A.; Yamazaki, Y.; Dubcovsky, Y.J.; Rogers, J.; Morris, C.; Apples, R.; Xia, X.C. Catalogue of Gene Symbols for Wheat. Available online: https://shigen.nig.ac.jp/wheat/komugi/genes/symbolClassList.jsp (accessed on 14 May 2021).
  31. Chernook, A.; Kroupin, P.; Karlov, G.; Soloviev, A.; Korshunova, A.; Rubets, V.; Igonin, V.; Divashuk, M. Effects of Rht-B1b and Ddw1 Dwarfing Genes in Two Connecting Populations of Spring Triticale under Greenhouse Experiment Conditions. Agriculture 2019, 9, 119. [Google Scholar] [CrossRef] [Green Version]
  32. Würschum, T.; Liu, W.; Busemeyer, L.; Tucker, M.R.; Reif, J.C.; Weissmann, E.A.; Hahn, V.; Ruckelshausen, A.; Maurer, H.P. Mapping dynamic QTL for plant height in triticale. BMC Genet. 2014, 15, 59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Braun, E.M.; Tsvetkova, N.; Rotter, B.; Siekmann, D.; Schwefel, K.; Krezdorn, N.; Plieske, J.; Winter, P.; Melz, G.; Voylokov, A.V.; et al. Gene Expression Profiling and Fine Mapping Identifies a Gibberellin 2-Oxidase Gene Co-segregating With the Dominant Dwarfing Gene Ddw1 in Rye (Secale cereale L.). Front. Plant. Sci. 2019, 10, 857. [Google Scholar] [CrossRef]
  34. Korzun, V.; Börner, A.; Melz, G. RFLP mapping of the dwarfing (Ddw1) and hairy peduncle (Hp) genes on chromosome of rye (Secale cereale L.). Theor. Appl. Genet. 1996, 92, 1073–1077. [Google Scholar] [CrossRef] [PubMed]
  35. Kroupin, P.; Chernook, A.; Karlov, G.; Soloviev, A.; Divashuk, M. Effect of Dwarfing Gene Ddw1 on Height and Agronomic Traits in Spring Triticale in Greenhouse and Field Experiments in a Non-Black Earth Region of Russia. Plants 2019, 8, 131. [Google Scholar] [CrossRef] [Green Version]
  36. Sun, L.; Yang, W.; Li, Y.; Shan, Q.; Ye, X.; Wang, D.; Yu, K.; Lu, W.; Xin, P.; Pei, Z.; et al. A wheat dominant dwarfing line with Rht12, which reduces stem cell length and affects gibberellic acid synthesis, is a 5AL terminal deletion line. Plant. J. 2019, 97, 887–900. [Google Scholar] [CrossRef]
  37. Chen, L.; Phillips, A.L.; Condon, A.G.; Parry, M.A.J.; Hu, Y.-G. GA-responsive dwarfing gene Rht12 affects the developmental and agronomic traits in common bread wheat. PLoS ONE 2013, 8, e62285. [Google Scholar] [CrossRef] [Green Version]
  38. Worland, A.J.; Sayers, E.J.; Borner, A. The Genetics and Breeding Potential of Rht12, a Dominant Dwarfing Gene in Wheat. Plant. Breed. 1994, 113, 187–196. [Google Scholar] [CrossRef]
  39. Würschum, T.; Langer, S.M.; Longin, C.F.H.; Tucker, M.R.; Leiser, W.L. A modern Green Revolution gene for reduced height in wheat. Plant. J. 2017, 92, 892–903. [Google Scholar] [CrossRef] [Green Version]
  40. Flintham, J.E.; Börner, A.; Worland, A.J.; Gale, M.D. Optimizing wheat grain yield: Effects of Rht (gibberellin-insensitive) dwarfing genes. J. Agric. Sci. 1997, 128, 11–25. [Google Scholar] [CrossRef]
  41. Wilhelm, E.P.; Mackay, I.J.; Saville, R.J.; Korolev, A.V.; Balfourier, F.; Greenland, A.J.; Boulton, M.I.; Powell, W. Haplotype dictionary for the Rht-1 loci in wheat. Theor. Appl. Genet. 2013, 126, 1733–1747. [Google Scholar] [CrossRef] [PubMed]
  42. Ericsson, K.; Rosenqvist, H.; Nilsson, L.J. Energy crop production costs in the EU. Biomass Bioenergy 2009, 33, 1577–1586. [Google Scholar] [CrossRef] [Green Version]
  43. Worland, A.J.; Korzun, V.; Röder, M.S.; Ganal, M.W.; Law, C.N. Genetic analysis of the dwarfing gene Rht8 in wheat. Part II. The distribution and adaptive significance of allelic variants at the Rht8 locus of wheat as revealed by microsatellite screening. Theor. Appl. Genet. 1998, 96, 1110–1120. [Google Scholar] [CrossRef]
  44. Mathews, K.L.; Chapman, S.C.; Trethowan, R.; Singh, R.P.; Crossa, J.; Pfeiffer, W.; Ginkel, M.; DeLacy, I. Global Adaptation of Spring Bread and Durum Wheat Lines Near-Isogenic for Major Reduced Height Genes. Crop. Sci. 2006, 46, 603–613. [Google Scholar] [CrossRef]
  45. Butler, J.D.; Byrne, P.F.; Mohammadi, V.; Chapman, P.L.; Haley, S.D. Agronomic Performance of Rht Alleles in a Spring Wheat Population across a Range of Moisture Levels. Crop. Sci. 2005, 45, 939–947. [Google Scholar] [CrossRef]
  46. Fahad, S.; Bajwa, A.A.; Nazir, U.; Anjum, S.A.; Farooq, A.; Zohaib, A.; Sadia, S.; Nasim, W.; Adkins, S.; Saud, S.; et al. Crop Production under Drought and Heat Stress: Plant Responses and Management Options. Front. Plant. Sci. 2017, 8, 1147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Ceccarelli, S.; Grando, S.; Maatougui, M.; Michael, M.; Slash, M.; Haghparast, R.; Rahmanian, M.; Taheri, A.; Al-Yassin, A.; Benbelkacem, A.; et al. Plant breeding and climate changes. J. Agric. Sci. 2010, 148, 627–637. [Google Scholar] [CrossRef]
  48. Singh, K.; Khanna-Chopra, R. Physiology and QTL anaylsis of coleoptile length, a trait for drought tolerance in wheat. J. Plant. Biol. 2010, 37, 1–9. [Google Scholar]
  49. Botwright, T.L.; Rebetzke, G.J.; Condon, A.G.; Richards, R.A. Influence of variety, seed position and seed source on screening for coleoptile length in bread wheat (Triticum aestivum L.). Euphytica 2001, 119, 349–356. [Google Scholar] [CrossRef]
  50. Neuweiler, J.E.; Maurer, H.P.; Würschum, T. Long-term trends and genetic architecture of seed characteristics, grain yield and correlated agronomic traits in triticale (×Triticosecale Wittmack). Plant. Breed. 2020, 139, 717–729. [Google Scholar] [CrossRef] [Green Version]
  51. Williams, E.; Piepho, H.P.; Whitaker, D. Augmented p-rep designs. Biom. J. 2011, 53, 19–27. [Google Scholar] [CrossRef] [PubMed]
  52. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  53. Piepho, H.P.; Büchse, A.; Emrich, K. A hitchhiker’s guide to mixed models for randomized experiments. J. Agron. Crop. Sci. 2003, 189, 310–322. [Google Scholar] [CrossRef]
  54. Butler, D.G.; Cullis, B.R.; Gilmour, A.R.; Gogel, B.J. ASReml-R Reference Manual; The State of Queensland, Department of Primary Industries and Fisherie: Queensland, Australia, 2009.
  55. Stram, D.O.; Lee, J.W. Variance components testing in the longitudinal mixed effects model. Biometrics 1994, 50, 1171. [Google Scholar] [CrossRef]
  56. Cullis, B.R.; Smith, A.B.; Coombes, N.E. On the design of early generation variety trials with correlated data. J. Agric. Biol. Environ. Stat. 2006, 11, 381–393. [Google Scholar] [CrossRef]
  57. Li, H.; Vikram, P.; Singh, R.P.; Kilian, A.; Carling, J.; Song, J.; Burgueno-Ferreira, J.A.; Bhavani, S.; Huerta-Espino, J.; Payne, T.; et al. A high density GBS map of bread wheat and its application for dissecting complex disease resistance traits. BMC Genom. 2015, 16, 216. [Google Scholar] [CrossRef] [Green Version]
  58. Yu, J.; Pressoir, G.; Briggs, W.H.; Vroh Bi, I.; Yamasaki, M.; Doebley, J.F.; McMullen, M.D.; Gaut, B.S.; Nielsen, D.M.; Holland, J.B.; et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 2006, 203–208. [Google Scholar] [CrossRef] [PubMed]
  59. Lipka, A.E.; Tian, F.; Wang, Q.; Peiffer, J.; Li, M.; Bradbury, P.J.; Gore, M.A.; Buckler, E.S.; Zhang, Z. GAPIT: Genome association and prediction integrated tool. Bioinformatics 2012, 2397–2399. [Google Scholar] [CrossRef] [Green Version]
  60. Money, D.; Gardner, K.; Migicovsky, Z.; Schwaninger, H.; Zhong, G.-Y.; Myles, S. LinkImpute: Fast and Accurate Genotype Imputation for Nonmodel Organisms. G3 Genes Genomes Genet. 2015, 5, 2383–2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Utz, H.F.; Melchinger, A.E.; Schön, C.C. 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. [Google Scholar] [CrossRef] [PubMed]
  62. Lynch, M.; Walsh, B. Genetics and Analysis of Quantitative Traits, 2nd ed.; Sinauer: Sunderland, MA, USA, 1998; ISBN 978-0878934812. [Google Scholar]
Figure 1. Development of plant height and developmental stage in the 129 registered cultivars dependent on their year of release.
Figure 1. Development of plant height and developmental stage in the 129 registered cultivars dependent on their year of release.
Plants 10 01592 g001
Figure 2. Circular Manhattan plot for plant height (PH) and developmental stage (DS). The red dashed lines indicate the Bonferroni-corrected significance threshold (p < 0.05) and the shaded area the QTL with pleiotropic effects on both traits.
Figure 2. Circular Manhattan plot for plant height (PH) and developmental stage (DS). The red dashed lines indicate the Bonferroni-corrected significance threshold (p < 0.05) and the shaded area the QTL with pleiotropic effects on both traits.
Plants 10 01592 g002
Figure 3. Fine-mapping of the identified QTL. (a) Genetic and (b) physical fine-mapping of the QTL on chromosomes 5A, 4B, and 5R. Red dots and lines indicate significantly associated markers after a Bonferroni multiple-test correction. (c) Linkage disequilibrium (r2) among the markers from the three QTL regions.
Figure 3. Fine-mapping of the identified QTL. (a) Genetic and (b) physical fine-mapping of the QTL on chromosomes 5A, 4B, and 5R. Red dots and lines indicate significantly associated markers after a Bonferroni multiple-test correction. (c) Linkage disequilibrium (r2) among the markers from the three QTL regions.
Plants 10 01592 g003
Figure 4. Effects of the QTL detected on chromosomes 5A, 4B, and 5R as well as their combinations on plant height, assessed in the 129 registered cultivars. ‘+‘ and ‘−’ indicate presence or absence of the height-reducing QTL alleles. If no presence/absence labeling is shown, the allelic state of the respective QTL was not considered.
Figure 4. Effects of the QTL detected on chromosomes 5A, 4B, and 5R as well as their combinations on plant height, assessed in the 129 registered cultivars. ‘+‘ and ‘−’ indicate presence or absence of the height-reducing QTL alleles. If no presence/absence labeling is shown, the allelic state of the respective QTL was not considered.
Plants 10 01592 g004
Figure 5. Number of cultivars with one or a combination of the plant height-reducing alleles at the detected QTL.
Figure 5. Number of cultivars with one or a combination of the plant height-reducing alleles at the detected QTL.
Plants 10 01592 g005
Figure 6. Temporal trends of the utilization of the identified QTL. Frequency of the height-reducing allele at the detected QTL qPH.5A, qPH.4B, and qPH.5R.C as well as their combinations in the registered cultivars dependent on their year of release.
Figure 6. Temporal trends of the utilization of the identified QTL. Frequency of the height-reducing allele at the detected QTL qPH.5A, qPH.4B, and qPH.5R.C as well as their combinations in the registered cultivars dependent on their year of release.
Plants 10 01592 g006
Table 1. Means and ranges of the best linear unbiased estimators as well as genotypic and genotype-by-environment interaction variance components of 846 diverse breeding lines and registered cultivars.
Table 1. Means and ranges of the best linear unbiased estimators as well as genotypic and genotype-by-environment interaction variance components of 846 diverse breeding lines and registered cultivars.
Plant HeightDevelopmental Stage
Min80.243.0
Mean100.853.4
Max126.659.8
LSD0.054.01.3
σ G 2 77.98 ***6.09 ***
σ G × E 2 7.38 ***1.21 ***
σ e 2 18.321.27
H20.810.80
***
Significantly different from zero at the 0.001 probability level.
Mean residual error variance across environment-trial combinations.
Table 2. Results of the genome-wide association mapping for plant height and developmental stage.
Table 2. Results of the genome-wide association mapping for plant height and developmental stage.
QTLMarkerChr.Pos. (cm)p-ValuepG Joint apG Single bα-Effect (Single Fit)pc
Plant height (42.16% pG total)
qPH.5AD105068725A150.62.3 × 10−180.4013.063.720.25
qPH.4BD4371530 d4B70.61.2 × 10−232.3418.97−4.020.36
qPH.4RD105237484R173.07.7 × 10−81.991.351.030.48
qPH.5R.AD105197775R990.47.3 × 10−71.794.562.260.23
qPH.5R.BD43484285R997.05.5 × 10−311.9119.12−3.940.41
qPH.5R.CS4341499 d5R1019.42.0 × 10−4829.3829.385.300.63
Developmental stage (29.31% pG total)
qEC.5AD105068725A150.63.3 × 10−103.2616.381.150.25
qEC.4BD43720074B70.62.7 × 10−150.0312.83−0.880.51
qEC.5R.BD43484285R997.02.4 × 10−1616.7416.74−1.020.37
qEC.5R.CS4341499 d5R1019.44.1 × 10−165.3519.431.190.63
a 
pG values obtained by a joint fit of all significant markers for the respective trait in a linear model; markers were ordered according to their p-value (lowest first).
b 
pG values obtained when each significant marker was fitted in a linear model for the respective trait.
c 
Frequency of trait-increasing allele, i.e., development stages (for earlier genotypes), plant height (for taller genotypes).
d 
Unmapped marker that was assigned to its most probable position based on its LD with mapped markers.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Trini, J.; Maurer, H.P.; Neuweiler, J.E.; Würschum, T. Identification and Fine-Mapping of Quantitative Trait Loci Controlling Plant Height in Central European Winter Triticale (×Triticosecale Wittmack). Plants 2021, 10, 1592. https://doi.org/10.3390/plants10081592

AMA Style

Trini J, Maurer HP, Neuweiler JE, Würschum T. Identification and Fine-Mapping of Quantitative Trait Loci Controlling Plant Height in Central European Winter Triticale (×Triticosecale Wittmack). Plants. 2021; 10(8):1592. https://doi.org/10.3390/plants10081592

Chicago/Turabian Style

Trini, Johannes, Hans Peter Maurer, Jan Eric Neuweiler, and Tobias Würschum. 2021. "Identification and Fine-Mapping of Quantitative Trait Loci Controlling Plant Height in Central European Winter Triticale (×Triticosecale Wittmack)" Plants 10, no. 8: 1592. https://doi.org/10.3390/plants10081592

APA Style

Trini, J., Maurer, H. P., Neuweiler, J. E., & Würschum, T. (2021). Identification and Fine-Mapping of Quantitative Trait Loci Controlling Plant Height in Central European Winter Triticale (×Triticosecale Wittmack). Plants, 10(8), 1592. https://doi.org/10.3390/plants10081592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop