Skip to main content

Identification of QTL hot spots for malting quality in two elite breeding lines with distinct tolerance to abiotic stress

Abstract

Background

Barley (Hordeum vulgare) is an important crop cultivated across the world. Drought is a major abiotic factor compromising barley yield worldwide, therefore in modern spring barley cultivars superior seed and malting quality characteristics should be combined with reasonable level of drought tolerance. Previously we have identified a number of barley lines demonstrating the superior yield performance under drought conditions. The aim of this work was to perform a QTL analysis of malting quality traits in a doubled haploid (DH) mapping population of two elite barley lines that differ in their reaction pattern to drought stress.

Results

A population of DH lines was developed by crossing two drought-tolerant elite breeding lines, Victoriana and Sofiara, exploiting distinct mechanism of drought tolerance, sustaining assimilation vs remobilization. The mapping population was assayed under field conditions at four distinct locations that differed in precipitation rate. DH lines were genotyped with the Illumina 9 K iSelect assay, and linkage map including 1782 polymorphic markers and covering a total map length of 1140 cM was constructed. The result of quantitative trait loci (QTL) analysis showed that majority of the traits were affected by several main effect QTL and/or QTL x environment (QE) interactions. In total, 57, 41, and 5 QTL were associated with yield-related traits, malting quality traits and seed quality traits, respectively. 11 and 29 of mapped QTL explained more than 10 and 5% of phenotypic variation, respectively. In several chromosomal regions co-localization between QTL for various traits were observed. The largest clusters were detected on chromosomes 3H and 4H.

Conclusions

Our QTL mapping results revealed several novel consistent genomic regions controlling malting quality which could be exploited in marker assisted selection. In this context, the complex QTL region on chromosome 3H seems of particular interest, as it harbors several large effect QTL.

Background

Barley (Hordeum vulgare) is an important crop cultivated throughout the world with Russia, France, Germany, Australia and Canada being the major producers. It is used as human food, livestock feed and for malt products. About 30% of the global barley production (133.5 million tons) is processed as malt by the brewing or the distilling industry [1, 2]. Therefore, one of the primary breeding goals is to improve the malting quality.

Malting quality is the result of complex interactions of numerous trait components, each of which is under control of multiple genes [3, 4]. Therefore, phenotypic selection of malting quality is labor-intensive and costly [5]. Marker-assisted selection is a promising alternative to improve malting quality [6] but necessitates detailed knowledge of its genetic architecture. Several linkage and association mapping studies have been conducted and QTL controlling components of barley malting quality, such as alpha- and beta-amylase activity, malt extract, malt soluble nitrogen, diastatic power, Kolbach index, wort β-glucans, wort viscosity, friability and apparent final attenuation have been identified [7,8,9,10,11,12]. Malting quality QTL were mapped to all seven barley chromosomes, however, the distribution was uneven and in some cases the formation of QTL clusters simultaneously affecting several malting traits have been observed [13,14,15]. Moreover, detected QTL were frequently dependent on genetic background and environmental conditions [16,17,18].

A few QTL explaining a large proportion of the genotypic variance have been used in marker-assisted selection (MAS). Igartua et al. [19] applied MAS for two QTL regions on chromosome 5H affecting several malting quality traits in the Harrington/TR306 population. Selection for the Harrington allele at target regions produced DH lines with superior malting quality characteristics such as low grain protein and β-glucan content, high diastatic power and malt extract. Ayoub et al. [20] introgressed a QTL from chromosome 5H affecting alpha-amylase activity originating from the cultivar Morex into the feed barley cultivar Labelle and observed an increase in alpha-amylase activity. Laido et al. [21] reported successful marker-assisted introgression of major QTL regions for grain protein content, malt extract, friability and viscosity on chromosome 1H in a population consisting of DH lines derived from winter and spring barley cross (Nure/Tremois). MAS for malting quality QTL may greatly benefit the introduction of novel diversity into elite germplasm. This is of particular importance if complex traits such as tolerance to biotic and abiotic stresses are to be introgressed from unadapted germplasm sources.

Previously, we have identified a number of barley lines showing senescence or stay-green phenotype and demonstrating the superior yield performance under drought conditions. The stay-green and senescing lines differed in their assimilation performance under drought stress. Moreover, it was revealed that senescing lines synthesized greater levels of ABA than the stay-green lines under short-term stress, and it continued to maintain a high level of ABA and ABA catabolites under long-term stress. This suggests that a greater flux in ABA metabolism in the senescing lines negatively affected water use efficiency and assimilation [22]. The objective of the current study was to perform a QTL analysis of malting quality traits in a biparental mapping population of two elite barley lines that differ in their reaction pattern to drought stress in order to identify genetic regions associated with twenty two malting-, grain-quality and yield-related traits.

Results

Phenotypic analysis

The BLUEs of the DH population were widely distributed for all examined traits approximating in most cases a normal distribution (Additional file 1: Figure S1; Table 1). For all traits we observed transgressive variation with several DH lines possessing phenotypic values lower or higher than their parental lines. The genotypic variance components\( \Big({\sigma}_G^2 \)) were for all traits except seed nitrogen content (SNC) significantly (P < 0.05) larger than zero (Table 2). For 14 out of the 22 measured traits, significant (P < 0.05) variance components were detected for genotype-by-environment interaction effects \( \Big({\sigma}_{G\times E}^2 \)). The ratio of \( {\sigma}_G^2:{\sigma}_{G\times E}^2 \) was minimum for grain yield (YLD) amounting to 0.27. The broad-sense heritability estimates (h2) ranged from 0.34 for SNC to 0.98 for seed length (SL).

Table 1 Summary statistics of the evaluated quantitative traits for the parents and the Sofiara/Victoriana DH population across environments
Table 2 Estimates of variance components for measured traits across environments

Correlation estimates

Pearson correlation coefficients were estimated between all pairs of traits (Table 3). Grain yield was negatively correlated (P < 0.05) with soluble nitrogen. Correlation coefficients between malting quality traits ranged from − 0.92 to 0.86. Seed quality traits such as SNC and seed starch content (SSC) were moderately negatively correlated with fine-grind malt extract (FGE) (r = − 0.49) and rootlet losses (RTL) (r = − 0.28), respectively.

Table 3 Pearson correlation coefficients between measured traits

QTL analysis

The 1782 mapped SNP markers detected a total of 380 independent loci. Hence, for QTL mapping 380 polymorphic high-quality SNP markers were selected, i.e. one marker per locus. The number of polymorphic markers ranged from 106 for chromosome 1H to 374 for chromosome 2H (Additional file 2: Table S1). Despite the high marker density, we observed some regions such as on the chromosome 1H exhibiting a low number of polymorphic markers. It is likely that these regions are identical by descent in both parents.

Genome-wide mapping resulted in 103 QTL with additive main effects and/or additive x environment interaction effects (QE) (Tables 4, 6, Fig. 1). No QTL were detected for seed carbon content. Twenty-nine QTL explained more the 5% and 11 QTL even more than 10% of the phenotypic variation. Variation explained by main effects (QTL) ranged from 0.41 to 68.77%, which was substantially higher than that of QTL x Environment interaction effects.

Table 4 Main-effect QTL for measured traits detected in the Sofiara/Victoriana DH population in six environments
Fig. 1
figure 1

Location of main-effect QTL on chromosomes. Triangle and cross represent favorable QTL alleles contributed by Victoriana and Sofiara, respectively

QTL and environment interactions for yield-related traits

The QTL analysis identified a total of 50 QTL for five yield-related traits. Seed area (SA) yielded the highest number of QTL (14) with phenotypic variation in the range of 0.17 to 15.87%. QTL of highest LOD score for SA were mapped on chromosome 2H and 5H explaining 15.87 and 4% of phenotypic variation, respectively. One significant environmental interaction effect was found for QSA-1-2. Twelve main effect QTL were observed for SL. The phenotypic variation explained by a single QTL ranged from 0.29 to 18.13%. A major QTL with a LOD value of 41.3 was found on chromosome 2H near marker SCRI_RS_209622. None of the QTL showed significant QE interaction. Thousand grain weight (TGW) and seed breadth (SB) were also among the traits with the highest number of QTL, 9 and 8 respectively. Loci affecting TGW were found to be located on chromosomes 2H, 4H, 5H, 6H and 7H and phenotypic variance explained by individual QTL ranged from 0.90 to 5.16%. Two QTL, QTGW-3-1 and QTGW-3-2 displayed significant QE interactions. QTL for SB were located on each chromosome except 1H and explained from 0.83 to 8.58% of the phenotypic variation. QSB-3 had only significant QE interaction effect. Four significant QTL affecting YLD were detected on three chromosomes (3H, 4H and 6H) and collectively explained 1.27% of phenotypic variation (Additional file 3: Table S2). QYLD-3-2 and QYLD-4 showed QE interaction only whereas QYLD-3-1 and QYLD-6 showed main effects only. The alleles for higher yield of QYLD-3-1 and QYLD-6 were contributed by Victoriana and Sofiara, respectively.

QTL for malting quality-related traits

Overall, 41 significant QTL were identified for twelve malting quality traits, ranging from 1 to 11 QTL per trait. There were five QTL that explained more than 10% of phenotypic variance for such traits as wort viscosity (VIS), soluble nitrogen (SNI), Kolbach Index (KOL), friability (FRI) and wort β-glucan content (BGL). Among all traits, the highest numbers of main effect QTL were detected for KOL (5), which were mapped to four chromosomes. The highest number of QE interaction effects were detected for SNI (9), located to three chromosomes.

Fourteen QTL with phenotypic variation in the range of 0.62 to 27.31% were detected for SNI. They were mapped to four chromosomes (2H, 3H, 4H and 5H). There were two main-effect QTL and nine QE interaction effects which did not show significant main effect. The proportion of collective phenotypic variation explained by main-effects and QE interaction effects were 33.44 and 13.06% (Additional file 3: Table S2). This suggests that the environmental interaction effects were important for the expression of SNI. The phenotypic effects of QE demonstrated large variation across the environments both in the magnitude and direction. The largest main effect QTL (QSNI-3-1) was mapped on chromosome 3H at 60.96 cM near the marker BK_08. The superior allele was contributed by Sofiara and explained 27.31% of the phenotypic variation with a LOD of 25.79.

For KOL six QTL were identified, explaining 0.85 to 13.74% of the phenotypic variance. At five loci the Sofiara alleles had increasing effects on the trait and only at one (QKOL-5-1) caused the reduction. Only single QTL located on chromosome 3H (QKOL-3) with significant main effect showed QE interactions.

A total of 6 QTL located on four chromosomal regions 2H, 3H, 4H and 5H were associated with BGL. Collectively, these four main effect QTL and QE interactions explained 50.80 and 11.99% of the phenotypic variation, respectively. Alleles for the QBGL-3, QBGL-4-1 and QBGL-4-2 for increased BGL content were contributed by Victoriana and at QBGL-5 the Sofiara allele increased the trait value. The major QTL, explaining up to 37.42% of phenotypic variation in BGL was located on 3H chromosome at 60.06 cM.

Four loci controlling FRI were identified on three chromosomes. The phenotypic variance explained by each QTL ranged from 3.93 to 33.81%. In all cases the alleles responsible for increasing the trait value were contributed by Sofiara. Collectively, the four main effect QTL explained 50.78% of the phenotypic variation. The largest QTL (QFRI-3) with LOD value 36.80 was identified on chromosome 3H at position 59.08 cM and explained 33.81% of the phenotypic variance. QFRI-3 also showed QE effect that explained 12.41% of the phenotypic variation.

Two main effect QTL were identified for malt protein content (PRO). The QTL were located on chromosomes 4H and 6H and explained 5.1 and 4.5% of the phenotypic variation, respectively. The Victoriana allele at each locus increased the protein content.

Four QTL for VIS were identified on three chromosomes. The phenotypic variance explained by each QTL ranged from 4.57 to 13.82%. Collectively, three main effect QTL and one QE interaction effect explained 29.25 and 1.67% of the phenotypic variation, respectively. In most cases the Victoriana alleles were responsible for increased trait value. The largest QTL (QVIS-3) with LOD value 13.91 was located on chromosome 3H at position 63.66 cM and explained 13.82% of the phenotypic variance.

For FGE, a total of three QTL were detected that mapped to chromosomes 1H, 3H and 7H. QFGE-3 and QFGE-7, which were main effect QTL collectively explained 21.37% of the phenotypic variation and positives alleles at these loci were contributed by Sofiara. One QE interaction effect (QFGE-1) was detected on chromosome 1H and explained 0.99% of phenotypic variance.

Three main effect QTL were detected for MMC. One was located on chromosome 3H at position 74.64 cM and explained 2.66% of the phenotypic variation. The other two were mapped on chromosome 4H at position 47.94 and 111.57 cM and explained 1.9 and 3.67% of the phenotypic variance, respectively. At the locus QMMC-3 and QMMC-4-2 positive effects were associated with Victoriana alleles and at QMMC-4-1 the Sofiara allele increased the trait value.

Three main effect QTL affecting malting losses were identified. Two QTL for RTL were located on chromosome 2H and 4H and one QTL for respiration losses (REL) on chromosome 2H. The percentage of phenotypic variation explained by individual QTL varied between 0.57 and 2.1%. The positive alleles for the QREL-2 and QRTL-2 were contributed by Sofiara and for QRTL-4 by Victoriana.

Only one QTL was detected for limit of attenuation (FLA) with LOD value 5.26. It was located on chromosome 3H near the marker BK_08 and explained 6.85% of phenotypic variance. The Sofiara allele contributed to the increased trait value.

QTL for seed-quality traits

A total of five significant QTL located on four chromosomes (2H, 3H, 4H and 7H) were detected for three seed-quality traits (Tables 4, 6, Fig. 1). The percentage of phenotypic variation explained by these QTL varied from 1.69 to 2.77%. For SSC, two main effect QTL (QSSC-2 and QSSC-3) were detected. The first QTL was located at the position 126.5 cM on 2H, explaining of the 2.38% of phenotypic variation. The second one was found on 3H at the position 66.4 cM and explained the 1.69% of the phenotypic variance. No significant QE interaction was detected for SSC. Two main effect QTL influencing SNC were observed on the chromosome 4H and 7H at the position 87.4 and 112.7 cM, respectively. These QTL accounted for 2.17–2.77% of the phenotypic variation in the trait. QTL on chromosome 4H (QSNC-4-2) showed no significant main effect while QE interaction effect was detected, explaining 2.01% of the phenotypic variance.

QTL for plant height and flowering time

Three QTL for plant height (PH) were found. QPH-5 explained 0.37% of variance. QPH-7 was located on 7H at position 114.27 cM near the marker SCRI_RS_194841 explaining 0.5% of phenotypic variation. QE interaction effect (QPH-2) was identified on chromosome 2H and explained 1.09% of phenotypic variation.

Four main effect QTL for flowering time (FT) with an LOD > 3 that collectively accounted for 46.66% of phenotypic variance were mapped. Victoriana alleles increased the time to flowering at the QFT-1 and QFT-3 and reduced it at the QFT-4 and QFT-6.

Co-localizations of QTL

For several traits hot spots of co-locating QTL were observed on all chromosomes except for 1H. It is well known that such clustering is a consequence of either close linkage or pleiotropy. The highest level of clustering occurred in two regions: in the centromeric region of chromosome 3H and at the bottom of chromosome 4H (Additional file 4: Figure S2). The cluster on chromosome 3H mainly comprises malting quality QTL for BGL, FRI, KOL, VIS and SNI (R2 > 10%). The cluster on 4H contained QTL for BGL, FRI and PRO together with several QTL for seed morphological traits. QTL for BGL were found to be most often co-localized with loci affecting FRI and VIS. These traits were highly correlated to each other (Table 3). A higher content of one component (BGL/VIS) corresponded to a lower content of the other (FRI). This was corroborated by the observed allelic effects, with Victoriana alleles increasing BGL content and decreasing FRI. In two chromosomal regions (3H, 4H) QTL for BGL co-localized with loci affecting FRI and VIS. In two other regions on chromosome 4H and 5H BGL QTL coincided with QTL for FRI or VIS, respectively. QBGL-3 co-localized not only with QTL for VIS and FRI but also with loci affecting FLA, SSC and SNI. The region involved in BGL variation in the distal region of chromosome 4H also accounted for variation in the protein content. Here, the Victoriana allele caused an increase in BGL and PRO content (Table. 5). Protein and β-glucan metabolism are independent biological processes therefore the co-localisation of QTL in this case is likely due to linkage between different genes than due to pleiotropy. Clusters of QTL affecting yield-related traits such as TGW, SL, SB and SA were observed on chromosomes 2H, 4H, 5H, 6H and 7H. The QTL cluster on chromosome 2H was also associated with SSC and SNI, while the cluster on the chromosome 5H involved QTL influencing BGL and VIS. We also observed that two regions on chromosomes 3H and 6H harbor coinciding QTL for YLD and FT, reflecting the interaction between these traits. The positive alleles for both traits at the cluster on chromosome 3H and 4H were derived from Victoriana and Sofiara, respectively.

Table 5 QTL x Environment interaction effect for measured traits detected in the Sofiara/Victoriana DH population in six environments

Discussion

In the present study comprehensive phenotyping, genotyping and QTL analysis of a population derived from two elite barley lines, that differ for drought tolerance behavior, has led to the identification of numerous QTL for yield and yield components as well as for seed and malting quality. On all chromosomes except 1H the hot spots carrying QTL for three or more traits where identified.

Main effect QTL

Our analysis resulted in identification of 85 main effect QTL for 21 traits. Victoriana alleles increased the trait value in 36 instances, while in 49 cases they were responsible for decreasing of trait expression. For about half of the traits parental alleles of both increasing and decreasing effects were observed. This underscores on the one hand the complex inheritance and on the other hand the transgressive potential of the cross. However, for two traits (PRO and SNI) the increased trait performance was associated with Victoriana alleles; while only Sofiara alleles were responsible for increasing of the trait value for FGE, FRI, FLA, total malting losses (TML), REL and PH. We will focus on QTL which are supposed important to understand the genetic architecture of a given trait and could be useful for further improvement of malting quality.

TGW is important yield component, and many investigations to determine genomic regions controlling this trait have been performed. QTL for TGW were mapped across all seven chromosomes in various bi-parental populations and diverse panels for association mapping studies [23,24,25,26]. In our study, a total of nine QTL were detected for TGW. As expected, high positive correlations with r values in the range of 0.64 to 0.84 were observed between TGW and SL, SB and SA confirming that TGW is mainly a function of seed volume. Indeed, co-localization of each QTL for TGW and at least one of the parameters for seed volume were observed. With regard to potential improving of TGW, the two QTL localized on long arm of 5H and short arm of 7H were of special interest, because they explained about 5% of phenotypic variation. The alleles for increasing TGW at QTGW-5-1 and QTGW-7 were contributed by Sofiara and Victoriana, respectively. The QTGW-5-1 and QTGW-7 were found to correspond to QTL detected previously by Rode et al. [25] and Schnaithmann and Pillen [26], respectively; whereas QTGW-5-2, corresponding to the marker BOPA1_1394–1222 at the peak position 169.4 cM is a newly identified locus. QTGW-6-1 and QTGW-6-2, found in our study, were both derived from Sofiara. They coincided with previously reported QTL for TGW [26], suggesting that these genetic factor(s) influencing TGW could operate in a different genetic background. QTGW-6-2 could be considered of special interest, as it was previously reported that loci at this position were able to increase TGW under nitrogen deficiency [26].

A total of 34 main effect QTL were determined for seed morphological traits. A half (50%) of QTL was contributed by Sofiara alleles while a second half by Victoriana. QSL-1 was coincident with QSA-1-1 and found to correspond to the locus controlling seed length detected by Kalladan et al. [27]. In the same genomic region the gene encoding drought-responsive MYB transcription factor was mapped [28]. The largest QTL, explaining more than 5% of phenotypic variation, were mapped on chromosomes 2H, 4H and 5H. QSL-2-1, mapped to the distal part of chromosome 2H at the marker SCRI_RS_209622 co-localized with QSA-2 and QTGW-2. For all these QTL the alleles increasing the trait value were contributed by Victoriana. That could suggest the presence of the loci with pleiotropic effect in this genomic region. Of interest to note that no QTL for SL were found in this region on chromosome 2H using the set of barley introgression lines S42ILs, produced from the cross Scarlett/ISR42–8. QTL affecting SA was identified on chromosome 2H, however, it was associated with different genomic region [26]. Kalladan et al. [27] recently reported mapping QTL for SL and TGW on chromosome 2H in the DH population Brenda/HS584. Due to absence of common markers, it is difficult to conclude whether they match to the ones identified in our study. However, a previously identified QTL for SL mapped before [27] was assigned to bin 8 according to the BinMap2005 (http://wheat.pw.usda.gov/cgi-bin/cmap), while QSL-2-1 detected in our study belongs to the bin 11. Therefore it could be considered as a newly identified locus. In addition, we have identified several QTL affecting such morphological seed traits as SL, SB and SA on chromosome 4H. QTL for SL mapped on chromosome 4H were not identified in previous studies, either in the S42IL or Brenda/HS584 populations. Therefore they could be considered as newly identified loci. QTL for SB and SA on chromosome 4H have been identified in the QTL study of wild barley introgression lines S42ILs. Most likely two loci (QSB-4 and QSA-4-1) identified in our study are identical to the ones detected by Schnaithmann and Pillen [26]. However, the additional loci controlling SA (QSA-4-2, QSA-4-3 and QSA-4-4) that all were contributed by Victoriana and have been mapped to the distal part of chromosome 4H in our investigation, represent novel QTL. On chromosome 5H, five loci influencing variation of all three examined seeds morphological traits were mapped. Positive alleles were exclusively contributed by Sofiara. Two of them, QSA5–2 and QSL-5, were mapped to identical position. Recently a locus controlling variation in SL was mapped on chromosome 5H [27], however, its position differs from the one detected in our study.

Our analysis revealed 30 main effect malting quality QTL, explaining 0.6–37.4% of the phenotypic variation. Seventeen favorable QTL alleles were contributed by Sofiara and 13 by Victoriana. 56.7% of those QTL were mapped to chromosomes 3H (8 loci) and 4H (9 loci). Two clusters comprising the largest effect QTL associated with 10 different malting quality traits have been identified on these chromosomes. This finding indicates the importance of these chromosomal regions for developing superior malting phenotype.

Malt extract and fermentability

In previous mapping studies, QTL for malt extract have been mapped on all seven chromosomes [15, 21, 29, 30]. FGE in the present population was controlled by two main effect QTL positioned on the short arm of 3H and 7H with the favorable alleles of Sofiara, accounting for 7% of the phenotypic variation. Malt extract QTL on the chromosome 3H was reported by Mohammadi et al. [11]. However, that locus was assigned to the long arm of chromosome 3H (bin 16) and could not be matched with QFGE-3 detected in our study. Position of QFGE-7 did not coincide with the localization of the loci controlling malt extract on chromosome7H reported earlier either [1, 13]. However, QFGE-7 was situated in proximity (about 15 cM) to QTL detected by Wang et al. [15] in the TX9425/Naso Nijo population. Comparison based on marker bins (BinMap2005) revealed that QFGE-7 was assigned to the bin 1 while the one detected by Wang et al. [15] to the bin 2 on chromosome 7H. This finding suggests that the top of chromosome 7 contains genetic factors affecting malt extract which operate in different genetic backgrounds.

Regarding the malt wort fermentability, a single QTL for FLA has been localized on the short arm of chromosome 3H at the marker BK_08 (bin 6). Loci influencing fermentability were reported on chromosomes 4H and 6H [31]. In addition, Swanston et al. [12] mapped QTL for fermentability on 3H at the sdw1 locus (bin11) which was suggested to be determined by ß-amylase activity. QFLA-3 detected in our study was contributed by Sofiara and coincided with fermentability locus detected in the Triumph/Morex population by Elía et al. [10]. Our correlation analyses (Table 3) suggested that a number of traits like PRO, FGE, VIS, FRI and BGL influenced FLA, with VIS having the strongest effect. Moreover, QTL for FRI, BGL, SNI, and VIS co-localized to the same genomic region, providing evidence that extent of endosperm modification and wort viscosity are major factors influencing fermentability.

Viskosity, friability and ß-glucan content

Such traits as VIS, FRI and BGL showed a high correlation to each other. Co-localizations QTL for at least two of these traits were observed on 3H, 4H and 5H. The favorable alleles for FRI QTL were exclusively contributed by Sofiara that was opposed to the loci controlling VIS and BGL which in their majority have been provided by Victoriana. Observed co-localization of loci affecting FRI and VIS could be explained by strong dependence of both traits on the process of degradation of ß-glucans that represent a major component of endosperm cell walls in barley [32]. The BGL QTL with the largest effect has been mapped near the centromeric region of chromosome 3H at the marker BK_08. QBGL-3 matched with QTL for malt ß-glucan reported by Han et al. [33] and was located in close proximity to the region affecting ß-glucan in barley grains [34]. QBGL-4-1 was found to be located in the same region as a QTL for malt ß-glucan reported earlier [33] while the location of QBGL-4-2 and QBGL-5 did not match with positions of the loci controlling BGL detected either in linkage or association mapping [34, 35]. None of mapped BGL QTL coincide with the position of known genes Glb1 and Glb2 on chromosome 1H and 7H, encoding (1,3;1,4)-β-glucanases. This could be explained by lack of genetic differences between the parents at these particular loci. Two loci QVIS-3, QVIS-5, controlling more than 5% of variability in viscosity, were found to match the QTL reported before [15, 21, 31]. In addition we have identified QTL for VIS at the short arm of the chromosome 4H linked to the marker BOPA2_12_30540. Previously, a locus controlling viscosity positioned on chromosome 4H was reported in the VB9524/ND11231_12 population [9], however, its precise location could not be clarified due to lack of common markers in our studies.

A total of four QTL were detected for FRI which were all derived from Sofiara, the parent possessing the favorable alleles (higher trait value). Two of them, QFRI-4-1 and QFRI-7, could be matched to the ones reported before [31, 36, 37]. However, the FRI in our population was mainly influenced by the locus on chromosome 3H which explained 33.8% of genetic variance. QFRI-3 belongs to the bin 6 and was found to be located in close proximity to FRI locus mapped on 3H (bin7) in the population Nure/Tremois [21]. Another interesting novel QTL for increased FRI maps at a distal position of chromosome 4 near marker BOPA1_ABC08009–1–2-304 and accounted for 5% of the phenotypic variation. QFRI-4-2 coincided with QBGL-4-2 and QPRO-4, suggesting that terminal region of chromosome 4H contains closely linked or pleiotropic genes affecting a range of important malting quality traits.

Protein content, soluble nitrogen, Kohlbach index

Nine QTL were found to be associated with PRO, SNI and KOL, traits that are used for characterization of nitrogen fractions of malt. A major QTL for SNI (QSNI-3-1) was mapped on chromosome 3H near the marker BK_08 in the same chromosomal region where locus controlling variation of SNI in Harrington/Mikamo Golden population was reported by Zhou et al. [38]. A novel QTL for SNI was detected on 2H near the marker BOPA1_11591–265. This QTL explained 1% of the phenotypic variation and was found to be co-located with loci affecting TGW, SA and SL. Genomic regions affecting malt protein content are known on all seven chromosomes [15, 38, 39]. However, two QTL QPRO-4 and QPRO-6 detected in the present study could not be matched with any previous findings. A total of five QTL were detected for KOL, and at four loci the increased KOL was contributed by Sofiara. Locations of QKOL-5-1, QKOL-5-2 and QKOL-7 coincide with position for the loci reported previously [15, 16]. The major QTL, QKOL-3, was mapped on 3H with closest marker SCRI_RS_115045 (bin 6). Locus controlling variation in KOL on chromosome 3H was reported in population 88Ab536/Morex [29]. This QTL was linked with the marker 1_1328 associated with bin 14 and therefore could not be matched with the one identified in current study. In addition, a minor QTL for KOL was mapped on 1H near the marker SCRI_RS_205669 (bin 5). In a similar region, a QTL affecting KOL was reported in wild barley introgression lines S42ILs [14]. In particular, line S24IL-102 bearing an introgression of about 85 cM in size showed a strong reduction of KOL up to 11.5%. However, it has to be noted that the large size of introgressed region covering nine bins obviates a precise comparison of QTL.

Seed starch and seed nitrogen content

QSNC-4-1 and QSNC-7 were found to match the QTL for grain protein content or N% reported before in the Vlamingh/Buloke and Brenda/HS584 populations, respectively [27, 30]. Locus controlling seed starch content on chromosome 3H (QSSC-3) was co-localized with QTL for a range of important malting quality traits like KOL, FRI, SNI and VIS. It coincides with QTL controlling SSC in the same region reported before by Kalladan et al. [27]. That suggests that the short arm of chromosome 3H contains genetic factors influencing a number of malting quality characters which stay active in different genetic backgrounds and are of primary importance for breeders. In addition, we have mapped QTL for SSC on 2H chromosome contributed by Victoriana. It was linked to the marker BOPA2_12_10735 (bin 11) and could not be matched to any QTL detected in any of above mentioned studies. The same region harbors a locus for SNI at which the Victoriana allele causes an increase in SNI. High SNI content would be considered a disadvantage, because it finally results in decreasing of malt extract. Moreover, QSSC-2 co-localized with QTL affecting yield-related parameters such as SA, SL and TGW. For all three QTL alleles increasing trait value were contributed exclusively by Victoriana.

QE interaction effects

From the practical point of view it is important to know how the identified loci performed in different environmental conditions. Our analysis revealed that 5 out of 86 main effect QTL showed QE interactions, in particularly QBGL-3, QVIS-3, QFRI-3, QKOL-3 and QSB-3. Phenotypic variation explained by QE interactions amounted to 12.7, 4.6, 11.4, 5.6 and 1.7%, respectively (Table 5). Loci for malting quality traits showed differences only in magnitude, while for QSB-3 differences were observed in both the magnitudes and directions of effects. QBGL-3 showed the largest variation among environments, with QE values of − 81.5, − 46.6 and − 190.7 in E1, E2 and E4, respectively. This suggests that expression of QBGL-3 has similar effect under all tested environments. The QE interaction effects of these QTL were substantially smaller than the corresponding main effects (Table 4). This suggests that the genetic effects of these loci were largely controlled by their main effect and thus they could be considered as primary targets in future breeding programs. In addition, QE effects were observed for 18 QTL with no significant main effect. The contribution of phenotype variations explained by QE interactions for these QTL ranged from 0.2 to 9.1%. The highest number of QE interactions (8 loci) was detected for the SNI, suggesting that this trait was highly sensitive to environment. The QSNI-3-2 with the largest effect, accounting for 9.1% of the phenotypic variance, demonstrated changes in magnitude of effect; while QSNI-5-2 had differential expression between E1 and E2, demonstrating negative and positive QE effects, respectively.

Candidate genes within the 3H QTL cluster

A number of coincidental QTL affecting two seed quality traits and five malt characters were identified on the short arm of chromosome 3H. These loci explained a relatively large proportion of phenotypic variation, in particular QSNI-3-1, QFRI-3 and QBGL-3 accounted for 27, 34 and 37% of phenotypic variance, respectively. Defined QTL cluster occupied relatively small chromosomal region and comprised at once five loci controlling variation of important malting quality characteristics. That makes it valuable for using in future breeding programs through marker assisted selection. In order to reveal possible candidate genes underlying identified QTL, the marker BK_08 associated with majority of detected loci was located on the POPSEQ genetic map of Morex × Barke [40]. The morex_contig_1580005, bearing the sequence of BK_08, was mapped at position 45.4 cM. The list of all annotated genes in the region between 40.4–50.4 cM was downloaded from publicly available BARLEX database [41]. Even a rough annotation-based assay of the genes located in this region (Additional file 5: Table S3) revealed a number of well-known malting quality related genes. These are the genes involved in starch, protein, lipid and cell wall metabolism. For instance, two genes MLOC_68194.1 and MLOC_36529.1 encoding aspartyl and serine proteases were found to be present in selected region. These two proteinase types are known to be involved in storage protein solubilisation, as well as in the regulation of the starch degrading enzyme, ß-amylase [42, 43]. In addition, four genes (AK375791, MLOC_32229.1, AK355146, AK375970) coding for heat shock proteins (HSP) were present in the region of interest. HSP might play a crucial role in protection of enzymes from degradation during the malting process especially at the stage of kilning. Association of HSP genes with such malting quality traits as malt extract, ß-amylase activity, friability, free amino acids and Kolbach Index have been previously reported [44, 45]. Interestingly, the target region also contained the gene MLOC_6671.1 encoding serine carboxypeptidase (Ser-CP). Ser-CPs play an important role in hydrolysis of storage proteins in germinating barley seeds. At least 6 members in barley were previously identified, all of which were found to be active during germination [46]. Transcript profiling of the gene Cxp1 (MLOC_55542.1), mapped to long arm of chromosomes 3H, in the Steptoe/Morex population revealed that expression QTL for respective gene co-localized with QTL affecting diastatic power [47]. Moreover, the region contains 18 genes encoding different types of transcription factors (TFs). They also could be considered as probable candidate genes, as these TFs might control expression of structural or regulatory genes playing key role in malting process. The fact that the members of the MYB gene family were found to be present in the list of potential candidates is particularly interesting. It has been previously reported that a MYB protein (GAMYB) expressed in cereal aleurone cells triggered transcriptional activation of α-amylase gene through binding to the TAACAAA box in its promoter [48]. The presence of elements of gibberellin responsive complex (GARC) have been reported within promoter regions of some other hydrolytic enzymes such as (1–3, 1–4)-β-glucanase and cathepsin B-like protease [49, 50]. Moreover, Gubler et al. [51] in a series of transient expression experiments showed that GAMYB protein indeed was able to promote expression of reporter gene fused with promoters of genes encoding EII(1–3, 1–4)-β-glucanase and cathepsin B-like protease.

Conclusions

In this study we have investigated agronomic, malting and seed quality characteristics of two high yielding elite breeding lines that differ in response to drought stress. Our QTL mapping revealed several novel loci controlling malting quality which could be exploited in breeding programs. In this context, the complex QTL region on chromosome 3H seems of particular interest, as it harbors several large effect QTL, which need to be validated in different genetic backgrounds. Moreover, identification and validation of the candidate genes will help to select for intralocus recombinants, since recombination of favorable alleles is required for further optimization of the genetic make-up of this locus.

Methods

Plant materials and field experiments

Our study is based on a population of 100 unique DH lines developed from the cross between the elite spring barley lines Sofiara and Victoriana (KWS, Germany). The parental lines were selected because they largely differ with respect to YLD, TGW, SNI and, BGL, and their response to drought stress. The DH lines and their parents were evaluated in replicated field trials at three locations in 2 years (Table 6; Additional file 6: Figure S3). The experimental design followed a lattice (Wohlde, Walewice) or a randomized complete block design (Gatersleben). The plot sizes ranged from 0.4 m2 to 5.5 m2. In Gatersleben, terminal drought stress was applied 1 week after anthesis using a rain-out shelter in addition to a control variant with absence of drought stress. Soil moisture was monitored in the drought stress trials twice a week in 10, 20, 30 and 40 cm depths by using 82 evenly distributed soil moisture sensors (SM300, PR2/4; Delta T devices Ltd., England).

Table 6 Environments in which the Sofiara/Victoriana DH population was evaluated

Yield-related traits

Plant height was measured before harvest from the base to the top of the plant and recorded in centimeters. SL, SB, SA and TGW were assessed using a sample of 400–500 grains and a digital seed analyser (Marvin; GTA Sensorik GmbH, Germany).

Starch, carbon and nitrogen content

Soluble carbohydrates were extracted three times with 80% ethanol at 60 °C from 0.4 g of homogenized grains tissue and ethanol-insoluble pellet was used for the quantification of starch content. Pellet was solubilized in 1 N KOH for 1 h at 95 °C. After neutralization with 5 N HCl, starch hydrolysis was achieved by incubation with amyloglucosidase (6u/ml) and a-amylase (10 u/ml) at 37 °C for 90 min. Glucose was quantified using enzymatic assay kit (R-Biopharm AG, Darmstadt, Germany). Carbon and nitrogen content were determined in ground and oven dried samples of grains using an elemental analyzer (vario EL III; Elementar analysensysteme GmbH Hanau, Germany).

Malting quality analysis

Grain samples (80 g per genotype) were malted in the malting facilities of the Research Institute for Raw Materials at the VLB Berlin e.V.. Micromalting conditions were the following: 5 h wet steep, 19 h air-rest, 2 h wet steep, 22 h air-rest, spray steeping until reaching a final moisture content of 45%, then germination for 96 h at 14.5 °C and 95 to 98% relative air humidity. Germinated samples were kilned with the following regime: 16 h at 50 °C, 1 h at 60 °C, 1 h at 70 °C and 5 h at 80 °C.

To assess malt quality following twelve parameters such as malt moisture content (MMC), FGE, VIS, PRO, SNI, KOL, FLA, FRI, BGL, TML, REL and RTL were determined according to the methods of MEBAK (Mitteleuropäische Brautechnische Analysenkommission). A brief outline of some malting traits is given in Potokina et al. [44].

Genotyping and construction of the genetic map

DNA samples from leaf tissue of the DH lines were extracted using the DNeasy Plant DNA miniprep kit (Qiagen, Germany), and scored with barley 9 K SNP iSelect Illumina array. Genetic map was constructed using the software package JoinMap 4.0 [52]. Linkage groups were established using minimum LOD values of 4. Monte Carlo maximum likelihood (ML) mapping algorithm was applied to determine the orders of markers within each linkage group. Recombination frequencies were converted to centimorgans (cM) using Haldane’s mapping function. The order of markers was compared with a published map [40].

Phenotypic data analysis and QTL mapping

We performed a two-stage analysis of phenotypic data [53]. The BLUEs of all genotypes in each environment were combined and a linear mixed model across environments was fitted with random genotype, environment, and residual effects. All linear mixed models were implemented using ASReml-R [54]. For QTL mapping we applied a composite mapping (CM) approach [55, 56]. For more details see Additional file 7: Data S1.

Abbreviations

BGL:

β-glucan

FGE:

Fine grind extract

FLA:

Limit of attenuation

FRI:

Friability

FT:

Flowering time

KOL:

Kolbach index

LOD:

Logarithm of odds

MAS:

Marker-assisted selection

MMC:

Moisture content

PH:

Plant height

PRO:

Protein content

QTL:

Quantitative trait locus

REL:

Respiration losses

RTL:

Rootlet losses

SA:

Seed area

SB:

Seed breadth

SCC:

Seed carbon content

SL:

Seed length

SNC:

Seed nitrogen content

SNI:

Soluble nitrogen

SNP:

Single nucleotide polymorphism

SSC:

Seed starch content

TGW:

Thousand-grain weight

TML:

Total malting losses

VIS:

Viscosity

YLD:

Grain yield

References

  1. Newton AC, Flavell AJ, George TS, Leat P, Mullholland B, Ramsay L, Revoredo-Giha C, Russell J, Steffenson B, Swanston JS, Thomas WTB, Waugh R, White PJ Bingham IJ. Crops that feed the world 4 barley: a resilient crop? Strengths and weaknesses in the context of food security. Food Secur. 2011;3:141–78.

    Article  Google Scholar 

  2. Newman RK, Newman CW. Barley for food and health science, technology, and products. Hoboken: Willey; 2008.

    Book  Google Scholar 

  3. Fox GP, Panozzo JF, Li CD, Lance RCM, Inkerman PA, Henry RJ. Molecular basis of barley quality. Aust J Agric Res. 2003;54:1081–101.

    Article  CAS  Google Scholar 

  4. Molina-Cano JL, Francesch M, Pérez-Vendrell AM, Ramo T, Voltas J, Brufau J. Genetic and environmental variation in malting and feed barley quality. J Cereal Sci. 1997;25:37–47.

    Article  Google Scholar 

  5. Hayes PM, Jones BL. Malting quality from a QTL perspective. In 8th Int. Barley Genet. Symp., Adelaide Convention Centre, Adelaide, South Australia 22–27 Oct. 2000;8:99–105.

  6. Han F, Romagosa I, Ullrich SE, Jones BL, Hayes PM, Wesenberg DM. Molecular marker-assisted selection for malting quality traits in barley. Mol Breed. 1997;3:427–37.

    Article  CAS  Google Scholar 

  7. Cai SG, Yu G, Chen XH, Huang YC, Jiang XG, Zhang GP, Jin XL. Grain protein content variation and its association analysis in barley. BMC Plant Biol. 2013;13:35–45.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Gutiérrez L, Cuesta-Marcos A, Castro AJ, von Zitzewitz J, Schmitt M, Hayes PM. Association mapping of malting quality quantitative trait loci in winter barley: positive signals from small germplasm arrays. Plant Genom. 2011;4:256–72.

    Article  Google Scholar 

  9. Emebiri LC, Moody DB, Panozzo JF, Read BJ. Mapping of QTL for malting quality attributes in barley based on a cross of parents with low grain protein concentration. Field Crops Res. 2004;87:195–205.

    Article  Google Scholar 

  10. Elía M, Swanston JS, Moralejo M, Casas A, Pérez-Vendrell AM, Ciudad FJ, Thomas WTB, Smith PL, Ullrich SE, Molina-Cano JL. A model of the genetic differences in malting quality between European and north American barley cultivars based on a QTL study of the cross triumph × Morex. Plant Breed. 2010;129:280–90.

    Article  CAS  Google Scholar 

  11. Mohammadi M, Blake TK, Budde AD, Chao S, Hayes PM, Horsley RD, Obert DE, Ullrich SE, Smith KP. A genome-wide association study of malting quality across eight US barley breeding programs. Theor Appl Genet. 2015;128:705–21.

    Article  PubMed  CAS  Google Scholar 

  12. Swanston JS, Thomas WTB, Powell W, Young GR, Lawrence PE, Ramsey L, Waugh R. Using molecular markers to determine barleys most suitable for malt distilling. Mol Breed. 1999;5:103–9.

    Article  CAS  Google Scholar 

  13. Hayes PM, Lui BH, Knapp SJ, Chen F, Jones B, Blake T, Franckowiak J, Rasmusson D, Sorrells M, Ullrich SE, Wesenberg D, Kleinhofs A. Quantitative trait locus effects and environmental interaction in a sample of north American barley germplasm. Theor Appl Genet. 1993;87:392–401.

    Article  PubMed  CAS  Google Scholar 

  14. Schmalenbach I, Pillen K. Detection and verification of malting quality QTLs using wild barley introgression lines. Theor Appl Genet. 2009;118:1411–27.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Wang J, Yang J, Zhang Q, Zhu J, Jia Q, Hua W, Shang Y, Li C, Zhou M. Mapping a major QTL for malt extract of barley from a cross between TX9425 x Naso Nijo. Theor Appl Genet. 2015;128:943–52.

    Article  PubMed  Google Scholar 

  16. Oziel A, Hayes PM, Chen FQ, Jones B. Application of quantitative trait locus mapping to the development of winter-habit malting barley. Plant Breed. 1996;115:43–51.

    Article  CAS  Google Scholar 

  17. Mather DE, Tinker NA, Laberge DE, Edney M, Jones BL, Rossnagel BG, Legge WG, Briggs KK, Irvine RB, Falk DE. Regions of the genome that affect grain and malt quality in a north American two-row barley cross. Crop Sci. 1997;37:544–54.

    Article  CAS  Google Scholar 

  18. Islamovic E, Obert DE, Budde AD, Schmitt M, Brunickii R, Kilian A, Chao S, Lazo GR, Marshall JM, Jellen EN, Maughan PJ, Hu G, Klos KE, Brown RH, Jackson EW. Quantitative trait loci of barley malting quality trait components in the stellar/01Ab8219 mapping population. Mol Breed. 2014;34:59–73.

    Article  CAS  Google Scholar 

  19. Igartua E, Edney M, Rossnagel BG, Spaner D, Legge WG, Scoles GJ, Eckstein PE, Penner GA, Tinker NA, Briggs KG, Falk DE, Mather DE. Marker-based selection of QTL affecting grain and malt quality in two-row barley. Crop Sci. 2000;40:1426–33.

    Article  CAS  Google Scholar 

  20. Ayoub M, Armstrong E, Bridger G, Fortin MG, Mather DE. Marker-based selection in barley for a QTL region affecting alpha amylase activity of malt. Crop Sci. 2003;43:556–61.

    Article  CAS  Google Scholar 

  21. Laidò G, Barabaschi D, Tondelli A, Gianinetti A, Stanca AM, Li Destri Nicosia O, NDi F, Francia E, Pecchioni N. QTL alleles from a winter feed type can improve malting quality in barley. Plant Breed. 2009;128:598–605.

    Article  CAS  Google Scholar 

  22. Seiler C, Harshavardhan VT, Reddy PS, Hensel G, Kumlehn J, Eschen-Lippold L, Rajesh K, Korzun V, Wobus U, Lee J, Selvaraj G, Sreenivasulu N. Abscisic acid flux alterations result in differential ABA signalling responses and impact assimilation efficiency in barley under terminal drought stress. Plant Physiol. 2014;164:1677–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Li JZ, Huang XQ, Heinrichs F, Ganal MW, Röder MS. Analysis of QTLs for yield components, agronomic traits, and disease resistance in an advanced backcross population of spring barley. Genome. 2006;49:454–66.

    Article  PubMed  CAS  Google Scholar 

  24. Schmalenbach I, Leon J, Pillen K. Identification and verification of QTLs for agronomic traits using wild barley introgression lines. Theor Appl Genet. 2009;118:483–97.

    Article  PubMed  CAS  Google Scholar 

  25. Rode J, Ahlemeyer J, Friedt W, Ordon F. Identification of marker trait associations in the German winter barley breeding gene pool. (Hordeum vulgare L). Mol Breed. 2012;30:831–43.

    Article  Google Scholar 

  26. Schnaithmann F, Pillen K. Detection of exotic QTLs controlling nitrogen stress tolerance among wild barley introgression lines. Euphytica. 2013;189:67–88.

    Article  CAS  Google Scholar 

  27. Kalladan R, Worch S, Rolletschek H, Harshavardhan VT, Kuntze L, Seiler C, et al. Identification of quantitative trait loci contributing to yield and seed quality parameters under terminal drought in barley advanced backcross lines. Mol Breed. 2013;32(1):71–90.

    Article  Google Scholar 

  28. Worch S, Rajesh K, Harshavardhan VT, Pietsch C, Korzun V, Kuntze L, Börner A, Wobus U, Roder MS, Sreenivasulu N. Haplotyping, linkage mapping and expression analysis of barley genes regulated by terminal drought stress influencing seed quality. BMC Plant Biol. 2011;11(1)

  29. Muñoz-Amatriaín M, Cistué L, Xiong Y, Bilgic H, Budde AD, Schmitt MR, Smith KP, Hayes PM, Muehlbauer GJ. Structural and functional characterization of a winter malting barley. Theor Appl Genet. 2010;120:971–84.

    Article  PubMed  CAS  Google Scholar 

  30. Walker CK, Ford R, Munoz-Amatrianın M, Panozzo JF. The detection of QTLs in barley associated with endosperm hardiness, grain density, grain size and malting quality using rapid phenotyping tools. Theor Appl Genet. 2013;126:2533–51.

    Article  PubMed  Google Scholar 

  31. von Korff M, Wang H, Leon J, Pillen K. AB-QTL analysis in spring barley: III identification of exotic alleles for the improvement of malting quality in spring barley (H vulgare ssp spontaneum). Mol Breed. 2008;21:81–93.

    Article  CAS  Google Scholar 

  32. Fincher GB. Molecular and cellular biology associated with endosperm mobilization in germinating cereal-grains. Annu Rev Plant Physiol. 1989;40:305–46.

    Article  CAS  Google Scholar 

  33. Han F, Ullrich SE, Chirat S, Menteur S, Jestin L, Sarrafi A, Hayes PM, Jones BL, Blake TK, Wesenberg DM, Kleinhofs A, Kilian A. Mapping of β-glucan content and β-glucanase activity loci in barley grain and malt. Theor Appl Genet. 1995;91:921–7.

    PubMed  CAS  Google Scholar 

  34. Mohammadi M, Endelman JB, Nair S, Chao SM, Jones SS, Muehlbauer GJ, Ullrich SE, Baik BK, Wise ML, Smith KP. Association mapping of grain hardness, polyphenol oxidase, total phenolics, amylose content, and beta-glucan in US barley breeding germplasm. Mol Breed. 2014;34:1229–43.

    Article  CAS  Google Scholar 

  35. Islamovic E, Obert D, Oliver R, Harrison S, Ibrahim A, Marshall JM, Stuthman DD. Genetic dissection of grain beta-glucan and amylose content in barley (Hordeum vulgare L). Mol Breed. 2013;31:15–25.

    Article  CAS  Google Scholar 

  36. Edney MJ, Mather DE. Quantitative trait loci affecting germination traits and malt friability in a two-rowed by six-rowed barley cross. J Cereal Sci. 2004;39:283–90.

    Article  CAS  Google Scholar 

  37. Li JZ, Huang XQ, Heinrichs F, Ganal MW, Röder MS. Analysis of QTLs for yield, yield components, and malting quality in a BC3-DH population of spring barley. Theor Appl Genet. 2005;110:356–63.

    Article  PubMed  CAS  Google Scholar 

  38. Zhou T, Takashi I, Ryouichi K, Naohiko H, Makoto K, Takehiro H, Kazuhiro S. Malting quality quantitative trait loci on a high-density map of Mikamo golden x Harrington cross in barley (Hordeum vulgare L). Mol Breed. 2012;30:103–12.

    Article  CAS  Google Scholar 

  39. Matthies IE, Malosetti M, Röder MS, van Eeuwijk F. Genome-wide association mapping for kernel and malting quality traits using historical european barley records. PLoS One. 2014;9(11):e110046.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Mascher M, Muehlbauer GJ, Rokhsar DS, Chapman J, Schmutz J, Barry K, Muñoz-Amatriaín M, Close TJ, Wise RP, Schulman AH, Himmelbach A, Mayer KFX, Scholz U, Poland JA, Stein N, Waugh R. Anchoring and ordering NGS contig assemblies by population sequencing (POPSEQ). Plant J. 2013;76:718–27.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Colmsee C, Beier S, Himmelbach A, Schmutzer T, Stein N, Scholz U, Mascher M. BARLEX - the barley draft genome explorer. Mol Plant. 2015;8:964–6.

    Article  PubMed  CAS  Google Scholar 

  42. Jones BL, Budde AD. Various malt endoproteinase classes affect wort soluble protein levels. J Cereal Sci. 2005;41:95–106.

    Article  CAS  Google Scholar 

  43. Schmitt MR, Marinac L. Beta-amylase degradation by serine endoproteinases from green barley malt. J Cereal Sci. 2008;47:480–8.

    Article  CAS  Google Scholar 

  44. Potokina E, Caspers M, Prasad M, Kota R, Zhang H, Sreenivasulu N, Wang M, Graner A. Functional association between malting quality trait components and cDNA array based expression patterns in barley (Hordeum vulgare L). Mol Breed. 2004;14:153–70.

    Article  CAS  Google Scholar 

  45. Lapitan NLV, Hess A, Cooper B, Botha AM, Badillo D, Iyer H, Menert J, Close TJ, Wright L, Hanning G, Tahir M, Lawrence C. Differentially expressed genes during malting and correlation with malting quality phenotypes in barley (Hordeum vulgare L). Theor Appl Genet. 2009;118:937–52.

    Article  PubMed  CAS  Google Scholar 

  46. Dal Degan F, Rocher A, Cameron-Mills V, von Wettstein D. The expression of serine carboxypeptidases during maturation and germination of the barley grain. Proc Natl Acad Sci U S A. 1994;91:8209–13.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Potokina E, Prasad M, Malysheva L, Röder MS, Graner A. Expression genetics and haplotype analysis reveal cis regulation of serine carboxypeptidase I (Cxp1), a candidate gene for malting quality in barley (Hordeum vulgare L). Funct Integr Genomics. 2006;6:25–35.

    Article  PubMed  CAS  Google Scholar 

  48. Gubler F, Kalla R, Roberts JK, Jacobsen JV. Gibberellin-regulated expression of a myb gene in barley aleurone cells: evidence for MYB transactivation of a high-pI α-amylase gene promoter. Plant Cell. 1995;7:1879–91.

    PubMed  PubMed Central  CAS  Google Scholar 

  49. Saleski N, Fincher GB. Developmental regulation of (1–3, 1–4) -b-glucanase gene expression in barley. Plant Physiol. 1992;99:1226–31.

    Article  Google Scholar 

  50. Cejudo FJ, Ghose TK, Stabel P, Baulcombe DC. Analysis of the gibberellin-responsive promoter of a cathepsin B-like gene from wheat. Plant Mol Biol. 1992;20:849–56.

    Article  PubMed  CAS  Google Scholar 

  51. Gubler F, Raventos D, Keys M, Watts R, Mundy J, Jacobsen JV. Target genes and regulatory domains of the GAMYB transcriptional activator in cereal aleurone. Plant J. 1999;17:1–9.

    Article  PubMed  CAS  Google Scholar 

  52. Van Ooijen JW. JoinMap® 4, software for the calculation of genetic linkage maps in experimental populations. Wageningen: Kyazma BV; 2006.

    Google Scholar 

  53. Möhring J, Piepho HP. Comparison of weighting in two-stage analysis of plant breeding trials. Crop Sci. 2009;49:1977–88.

    Article  Google Scholar 

  54. Gilmour AR, Gogel B, Cullis BR, Thompson R. ASReml user guide release 3.0. Hemel Hempstead: VSN International Ltd; 2009.

    Google Scholar 

  55. Jansen RC. Controlling the type I and type II errors in mapping quantitative trait loci. Genetics. 1994;138:871–81.

    PubMed  PubMed Central  CAS  Google Scholar 

  56. Zeng ZB. Precision mapping of quantitative trait loci. Genetics. 1994;136:1457–68.

    PubMed  PubMed Central  CAS  Google Scholar 

  57. Jansen RC, Van Ooijen JW, Stam P, Liste C, Dean C. Genotype-by-environment interaction in genetic mapping of multiple quantitative trait loci. Theor Appl Genet. 1995;91:33–7.

    Article  PubMed  CAS  Google Scholar 

  58. Schwarz G. Estimating the dimension of a model. Ann Stat. 1978;6:461–4.

    Article  Google Scholar 

  59. Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138:963–71.

    PubMed  PubMed Central  CAS  Google Scholar 

  60. SAS Institute. SAS/STAT: user’s guide, version 93 SAS Institute Inc Cary. NC. 2009;

  61. R Core Team (2014) R: A language and environment for statistical computing R Foundation for Statistical Computing, Vienna, Austria. URL http://wwwR-projectorg.

Download references

Acknowledgments

We are grateful to Mandy Püffeld and Jana Lorenz for excellent technical assistance. We thank two anonymous reviewers for their valuable comments on this manuscript.

Funding

This research was funded by the Bundesministerium für Bildung and Forschung (BMBF). The funding bodies had no role in the experiment design, data analysis, preparation of the manuscript and decision to publish.

Availability of data and materials

All data generated during the current study are included in this published article.

Author information

Authors and Affiliations

Authors

Contributions

AK performed phenotypic measurements, genotypic data analysis, statistical analysis and drafted the manuscript; YJ performed statistical analysis and gave a significant support during the phenotypic and QTL analyses; CS, KS performed phenotypic measurements and revised manuscript; SK participated in design and carried out the collection of data; JCR critically revised study and contributed to the interpretation of the results; VK and AG designed, coordinated the study and revised manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Andreas Graner.

Ethics declarations

Ethics approval and consent to participate

The plant materials come from germplasm resources of the KWS LOCHOW GmbH. All experimental trials were conducted according to the national legislation.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Frequency distribution of the BLUEs of the 100 DH for the measured traits. (PDF 257 kb)

Additional file 2:

Table S1. Genetic positions of markers. (XLSX 48 kb)

Additional file 3:

Table S2. The collective percent of phenotypic variation explained by main-effect QTL and QE interactions for the trait. (XLSX 12 kb)

Additional file 4:

Figure S2. Locations of main-effect QTL on 3H (a) and 4H (b) chromosomes. (PDF 183 kb)

Additional file 5:

Table S3. The list of candidate genes associated with largest QTL cluster on chromosome 3H. (XLSX 34 kb)

Additional file 6:

Figure S3. Climate conditions in rainfed field located in Wohlde (a-b), Walewice (c-d) and Gatersleben (e). (PDF 409 kb)

Additional file 7:

Data S1. Experimental procedures related to phenotypic data analysis, and QTL mapping [57,58,59,60,61]. (DOCX 14 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kochevenko, A., Jiang, Y., Seiler, C. et al. Identification of QTL hot spots for malting quality in two elite breeding lines with distinct tolerance to abiotic stress. BMC Plant Biol 18, 106 (2018). https://doi.org/10.1186/s12870-018-1323-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-018-1323-4

Keywords