- Research article
- Open Access
Identification of QTLs controlling grain protein concentration using a high-density SNP and SSR linkage map in barley (Hordeum vulgare L.)
BMC Plant Biology volume 17, Article number: 122 (2017)
Grain protein concentration (GPC) is a major determinant of quality in barley (Hordeum vulgare L.). Breeding barley cultivars with high GPC has practical value for feed and food properties. The aim of the present study was to identify quantitative trait loci (QTLs) for GPC that could be detected under multiple environments.
A population of 190 recombinant inbred lines (RILs) deriving from a cross between Chinese landrace ZGMLEL with high GPC (> 20%) and Australian cultivar Schooner was used for linkage and QTL analyses. The genetic linkage map spanned 2353.48 cM in length with an average locus interval of 2.33 cM. GPC was evaluated under six environments for the RIL population and the two parental lines. In total, six environmentally stable QTLs for GPC were detected on chromosomes 2H (1), 4H (1), 6H (1), and 7H (3) and the increasing alleles were derived from ZGMLEL. Notably, the three QTLs on chromosome 7H (QGpc.ZiSc-7H.1, QGpc.ZiSc-7H.2, and QGpc.ZiSc-7H.3) that linked in coupling phase were firstly identified. Moreover, the genetic effects of stable QTLs on chromosomes 2H, 6H and 7H were validated using near isogenic lines (NILs).
Collectively, the identified QTLs expanded our knowledge about the genetic basis of GPC in barley and could be selected to develop cultivars with high grain protein concentration.
Protein is an essential nutrient for the survival of humans and animals [1, 2]. Protein in mature cereal grains, in particular, provides a substantial portion of the world’s plant protein, and its concentration determines the nutritional quality and end use properties of the grain [3, 4]. Barley (Hordeum vulgare L.) is one of the earliest domesticated crops in the world. Approximately 25% of its production has relatively lower GPC and is suitable for malting and brewing, while the remaining 75% with relatively higher GPC is used for feed and food (http://faostat.fao.org/). Hence, there is increasing need for breeding barley cultivars with high GPC, but this has been hindered by the relatively low heritability of GPC due to the significant interaction between environmental and genetic factors [5, 6]. Based on a statistical methodology, the genetic factors (quantitative trait loci, QTLs) that involved in determination of GPC can be elucidated . Thus, identification and utilization of environmentally stable QTLs associated with GPC will provide an alternative but promising strategy for high GPC barley breeding.
To date, numerous studies have been conducted on dissecting the genetic basis of GPC, and QTLs have been mapped on all seven barley chromosomes. In particular, several consensus QTLs mapped on chromosomes 2H, 4H, 5H, 6H, and 7H have been repeatedly detected by multiple studies [8,9,10,11,12,13,14,15]. For example, two QTLs on chromosomes 5HS and 6HS located in the Bmac0096-Bmag0323 and ABG458-HVM74 intervals, respectively, have been repeatedly detected [10,11,12, 16]. Moreover, these two loci have also been identified by genome-wide association studies (GWAS) [17,18,19]. In addition, two genes (HvNAM1 and HvNAM2) on chromosomes 6H and 2H in barley, which were suggested to be orthologous to TtNAM-B1, contributed a substantial effect on GPC [17, 20]. Notably, a recent study revealed that a single nucleotide polymorphism (SNP) within the second intron of HvNAM2 was associated with GPC, which is useful in developing high quality barley cultivars . Although these identified QTLs/genes for GPC that could be expressed under multiple environments might be valuable for GPC improvement in barley, most of the genetics studies focused on breeding and selection for low-protein barley [21, 22].
A saturated genetic linkage map will improve the precision of QTL localization and estimation of phenotypic variance, especially for some small and medium-sized QTLs . Due to the abundance of SNPs in plant genome, SNP markers have been widely used in genetic linkage map construction [24,25,26]. High-density SNP linkage maps have been largely used in QTL detection for yield and quality in barley [27,28,29]. However, QTL mapping for GPC based on a high-density SNP map has rarely been reported. Here, to identify QTLs for GPC, a RIL population including 190 lines derived from a cross between the Chinese landrace ZGMLEL with high GPC (> 20%) and the Australian cultivar Schooner was used for linkage and QTL analyses. Furthermore, near-isogenic line (NIL) populations were developed to validate the environmentally stable QTLs.
A RIL population (generations F9 to F11) containing 190 RILs derived from two spring barley varieties, ZGMLEL and Schooner, was employed to identify QTLs controlling for GPC. ZGMLEL is a hull-less landrace with high GPC, while Schooner is a hulled cultivar with low GPC. All the RILs and their parental lines were kindly provided by Dr. Yawen Zeng (Yunnan Academy of Agricultural Sciences, China).
For NIL development, one RIL line (RIL7) was crossed with the recurrent parent (Schooner). Because the QTLs on chromosomes 6H (i.e. QGpc.ZiSc-6H.2, QGpc.ZiSc-6H.3 and QGpc.ZiSc-6H.4) and 7H (i.e. QGpc.ZiSc-7H.1, QGpc.ZiSc-7H.2 and QGpc.ZiSc-7H.3) were linked in coupling phase, the QTL clusters were introgressed into the Schooner background as a unity, respectively. After three backcross generations (BC3), individuals that solely exhibited heterozygosity at one QTL region were self-pollinated to produce its corresponding BC3F2 populations. Finally, three NIL populations, that is BC3F2-I (region 2H), BC3F2-II (region 6H) and BC3F2-III (region 7H), were developed for the validation their corresponding QTLs. The number of progenies in BC3F2-I, BC3F2-II and BC3F2-III populations were 249, 205 and 213, respectively.
Field experiments were carried out in three locations, including Shangzhuang Experiment Station of CAU (China Agricultural University) in Beijing, Wangtaibao Experiment Station of NAAFS (Ningxia Academy of Agriculture and Forestry Sciences) in Ningxia Hui Autonomous Region, and Dishang Experiment Station of HAAFS (Hebei Academy of Agriculture and Forestry Science) in Hebei Province. The RIL population and the two parents were grown during three growing seasons from 2013 to 2015, providing data for six environments. Location-year information and corresponding weather data are presented in Additional file 1: Table S1. In field trials, each plot consisted of 2 rows that were 2 m long with approximately 20 plants per row. The middle ten plants in each line were bulk-harvested at maturity and measured for grain protein concentration (GPC).
Three BC3F2 populations for QTL validation were planted in Beijing (2016). Individuals were grown in 2-m-long rows with a 0.25-m row spacing. Within each row, 15 plants were evenly sown. At maturity, all the panicles were harvested from single-plant and sun-dried. Grain protein concentration (GPC), grain yield (GY) and thousand grain weight (TKW) were scored on a single-plant basis.
The field experiments were in accordance with local practice. All the trails were conducted under optimum irrigation. Nitrogen (N) was supplied at a rate of 220 kg/ha, including 70 kg/ha of N as diammonium phosphate and 80 kg/ha of N as urea applied before sowing. In addition, 70 kg/ha of N as urea was applied at booting stage.
Phenotypic evaluation and statistical analysis
Mature grains of RIL population and BC3F2 populations were ground to a powder using a Cyclotec 1093 sample mill (Hoganas City, Sweden). Then, the ground powder was dried to a constant mass in an 80 °C oven. The total nitrogen content was determined using the Kjeldahl method with a FOSS Kjeltec ™ 2300 and then the GPC was calculated using a factor of 5.83 . GY and TKW of the NIL populations were measured on a single-plant basis. TKW was determined using a camera-assisted phenotyping system, which was provided by Hangzhou Wanshen Detection Technology Co., Ltd. (Hangzhou, China).
The basic statistical analysis was performed using SPSS version 20.0 (SPSS, Chicago, IL, USA). The Shapiro-Wilk test was conducted using R software (V. 3.2.2) for the normality test. The best linear unbiased prediction (BLUP) for GPC across the six environments was calculated using SAS® V.8 (SAS Institute Inc. 2000) with the PROC MIXED procedure. Under the random-effect model, environments were treated as fixed, and genotype and genotype-environments interactions were considered as random factors. The broad sense heritability (h B 2) on a family basis was calculated using SAS® V.8 (SAS Institute Inc. 2000) with the PROC GLM procedure, which was calculated according to the following formula: h B 2 = ó g 2 / (ó g 2 + ó ge 2 /n + ó 2 /nr) where ó g 2 = genotypic variance, ó ge 2 = genotype by environmental variance, ó 2 = the residual error variance, n = the number of environments, and r = number of replicates.
Linkage and QTL analyses
The RIL population was genotyped using the barley 9 K SNP chip developed from the RNA-seq data of barley varieties . Additionally, a total of 21 polymorphic SSR markers were employed to genotype the RIL population, and most of the SSR sequences were obtained from http://wheat.pw.usda.gov/GG3/. Only markers with less than 5% missing data were selected for map construction. The genetic linkage map was constructed using RECORD 2.0  and JoinMap 4.0 . Markers with identical segregation were first removed using RECORD 2.0. After removing the redundancy, the unique markers were grouped using JoinMap 4.0 with a LOD value of 10. Finally, the marker order was established using the maximum likelihood mapping algorithm and the map distance was calculated using the Kosambi mapping function. The probe sequences of the SNP assigned to barley chromosomes were queried using the BLAST algorithm against barley reference genome sequence to locate chromosomal positions with a cutoff criterion of E-value ≤1e-10. The quality of the genetic map was validated using the alignments between SNP map and barley reference genome. Only the best hit of the SNP against the reference genome was selected for the collinearity analysis when the SNP was located to multiple paralogous positions in the genome.
The average GPC data in each environment and BLUP values across six environments were collected for QTL analysis. WinQTLCart2.5 software with the composite interval mapping (CIM) method was used to identify QTLs for GPC. The walking speed was set to 1 cM. Model 6 was chosen for QTL analysis, with 5 control markers and 10 cM window size defaults. The LOD threshold was set via 1000 permutations at P ≤ 0.05, and these QTLs were considered “identified QTLs”. A 2-LOD support with a 99% confidence level was chosen for each identified QTL. The identified QTLs detected in different environments with overlapping confidence intervals were regarded as the same in this study. The QTLs were named following the rules of Blake and Blake .
The flanking markers of the QTLs were employed to define the target region, which could be used to compare to the barley Genome Zipper developed by Mayer et al. . Gene sequences of three grasses (rice, sorghum, and Brachypodium) were used as queries to blast against the database “assembly_WGSMorex” at IPK Barley BLAST Server (http://webblast.ipk-gatersleben.de/barley_ibsc/). The Morex contigs with best hit were employed to search for simple sequence repeat using the SSR Hunter software. Finally, the selected sequences were used to design SSR markers using the Primer3 software (http://frodo.wi.mit.edu/primer3/).
Analysis of GPC
The basic statistics of minimum, maximum, mean standard deviation and coefficient of variation for GPC in six environments are listed in Table 1. The GPC of parent ZGMLEL ranged from 20.52% to 22.88% in the six environments evaluated, which was significantly higher (P < 0.01) than that of Schooner (16.35–17.20%) (Table 1; Additional file 2: Table S2). Moreover, the 190 RILs exhibited a wide range of variation in GPC, with coefficients of variation (CVs) ranging from 6.94% to 7.80% in the six environments. The Shapiro-Wilk for testing normality was performed for GPC based on the mean value collected from six environments (Fig. 1). In all of the six environments, GPC showed normal distribution, suggesting a quantitative nature of GPC in barley. Remarkably, the broad sense heritability (h B 2) for the GPC of the RILs was 80.67%, indicating that the GPC variance was mostly determined by genetic factors.
Construction of a high-density genetic linkage map
Of the 7864 SNP markers on the chip, 1526 (19.40%) were polymorphic between ZGMLEL and Schooner. After removing 53 SNPs with over 5% missing data, we used 1473 SNP markers and 21 polymorphic SSR markers to construct the linkage map. The resultant linkage map comprised nine linkage groups that contained 1011 unique loci and spanned 2353.48 cM. These linkage maps had an average locus interval of 2.33 cM (Table 2, Additional file 3: Table S3, Additional file 4: Figure S1). The identity and polarity of linkage groups were determined by BLAST against the barley reference sequence databases .
To further validate the quality of the map, SNP flanking sequences were employed to align with the barley reference sequence. Of 1473 SNP markers, 1411 (95.79%) were successfully assigned to the barley genome (Table 2; Additional file 5: Table S4). A good collinearity of the genetic map with the barley reference genome sequence was observed (Fig. 2), indicating a high quality of the genetic linkage map. However, several chromosome intervals were inconsistent with the reference genome sequence, i.e., chromosomes 2H at 76.43–204.46 cM, 4H at 125.91–127.33 cM, and 5H-1 at 100.55–146.27 cM.
Of the 5.1 Gb size of barley genome, 3.1 Gb has been successfully anchored to the physical map through population sequencing (POPSEQ) . We calculated the coverage ratio for each barley chromosome. Chromosomes 1H, 3H, 4H and 7H displayed similar coverage ratios at 99.93, 95.40, 99.36, and 99.64%, respectively, and chromosomes 2H (65.51%), 5H (63.13%) and 6H (63.52%) exhibited lower ratios (Table 2).
QTL mapping of GPC
A total of 17 QTLs were detected for GPC, which are randomly distributed among chromosomes 2H (3 QTLs), 4H (3 QTLs), 5H (3 QTLs), 6H (4 QTLs), 7H (4 QTLs) (Fig. 3; Table 3). With an exception of QGpc.ZgSc-6H.1, ZGMLEL contributed effects for increased GPC at other 16 QTLs. These 16 significant QTL had LOD values ranging from 2.51 to 15.51 and explained the GPC variation from 2.4% to 19.86%. Schooner contributed effects for increased GPC at the locus of QGpc.ZgSc-6H.1, which had a LOD value of 3.37 and accounted for 5.70% of GPC variation. This indicated that the favorable alleles for increased GPC were mainly inherited from the feed barley ZGMLEL. In the present study, significant QTLs that could be detected in no less than three environments as well as in the combined analysis were defined as environmentally stable QTLs. According to this criterion, 6 of 17 significant QTLs were environmentally stable QTLs, which were identified on chromosomes 2HL (1), 4HS (1), 6HL (1), and 7HS (3).
Three QTLs associated with GPC were detected on chromosome 2HL. Only the one flanked by SCRI_RS_171032 and BOPA2_12_30901, QGpc.ZgSc-2H.1, was considered being environmentally stable, which could be detected under three environments. QGpc.ZgSc-2H.1 had a LOD value of 3.78 and explained 4.47% of the GPC variation for the combined analysis. The other two putative QTLs, QGpc.ZgSc-2H.2 and QGpc.ZgSc-2H.3, were observed in one environment and explained 5.69–8.58% of GPC variation.
Among three significant QTLs for GPC detected on chromosome 4H, one stable QTL was identified on chromosome 4HS and designated QGpc.ZgSc-4H.1. This QTL had a minor effect, with a LOD value of 3.55, and explained 4.20% of the GPC variation for the combined analysis. QGpc.ZgSc-4H.2 and QGpc.ZgSc-4H.3 could be detected only in one environment, which had LOD and R 2 values ranging from 4.58 to 5.59 and 5.68% to 6.10%, respectively.
Chromosome 5HL carried three putative QTLs significantly associated with GPC, which were designated QGpc.ZgSc-5H.3, QGpc.ZgSc-5H.1, and QGpc.ZgSc-5H.2. QGpc.ZgSc-5H.3 had a LOD value of 3.31 and explained 3.37% of the GPC variation for the combined analysis. The other two QTLs had the LOD values ranging from 4.13 to 4.71, and explained the GPC variation from 6.56% to 7.08%.
Four significant QTLs were identified on chromosome 6HL and were designated QGpc.ZgSc-6H.1, QGpc.ZgSc-6H.2, QGpc.ZgSc-6H.3, and QGpc.ZgSc-6H.4. Schooner contributed the effects for increasing GPC at QGpc.ZgSc-6H.1, and ZGMLEL contributed increased GPC at the other three loci. QGpc.ZgSc-6H.3 was the most stable QTL for GPC, which could be detected in all six environments and explained as much as 17.90% of the GPC variation for the combined analysis. Two putative QTLs, QGpc.ZgSc-6H.2 and QGpc.ZgSc-6H.4, were identified in two environments and explained 8.45–18.78% of GPC variation.
Four QTLs were found to be significantly associated with GPC on chromosome 7HS. Among these significant QTLs, three were environmentally stable QTLs, and they were designated QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3. These three stable QTLs were detected in three to five environments, explaining 12.37–13.29% of GPC variation for the combined analysis. The last putative QTL, QGpc.ZgSc-7H.4, was detected at E6 and accounted for 7.48% of GPC variation.
To study the three genomic regions on chromosomes 2H, 6H and 7H that possessing environmentally stable QTLs for GPC in more depth, three BC3F2 populations were developed and named BC3F2-I, BC3F2-II and BC3F2-III, respectively. Accordingly, three sets of SSR markers were used for foreground selection, i.e., 2L10, 2L11 and 2L12 for BC3F2-I, 6L89, 6L155 and 6L147 for BC3F2-II, and 7S40, 7S69, 7S87 and 7S89 for BC3F2-III (Additional file 6: Figure S2, Additional file 7: Table S5). A total of 76 SSRs were used for background selection of the BC3F1 individuals. Finally, three BC3F1 individuals that exhibit heterozygosity solely at genomic regions 2H, 6H or 7H were selfed to produce their corresponding BC3F2 populations. These three BC3F1 individuals shared 93.42, 92.10 and 94.74% similarities in genetic background with the recurrent parent, respectively.
To determine whether the stable QTLs affect the GPC in NIL populations, we compared the GPC between two homozygous groups. Based on the genotype of flanking markers (2L10 and 2L11 on 2H, 6L89 and 6L147 on 6H, 7S87 and 7S40 on 7H), two homozygous groups, namely, ZGMLEL homozygous (ZZ) and Schooner homozygous (SS) were classified in each NIL population. The evaluation results for GPC showed that plants with ZZ genotype in BC3F2-I, BC3F2-II and BC3F2-III had an average GPC of 13.82, 14.18 and 14.20%, respectively. In contrast, plants with SS genotype in BC3F2-I, BC3F2-II and BC3F2-III had an average GPC of 13.15, 13.19 and 13.48%, respectively, which is similar to the recurrent parent, Schooner (13.32%) (Additional file 8: Table S6). Based on the GPC value, highly significant difference was found between two homozygous genotypes in each NIL population (P < 0.01) (Table 4, Additional file 8: Table S6). The allelic effects of the three populations were in the same direction as the original allele, with alleles from ZGMLEL increasing GPC. These results suggested that the stable QTLs on chromosomes 2H, 6H and 7H had significant effect on GPC, which was in agreement with the detection in RIL population.
Previous studies reported that there was negative relationship between GPC and grain yield . Thus, we measured thousand kernel weight (TKW) and grain yield per plant (GY) for the three BC3F2 populations but found no significant difference for TKW and GY (Table 4; Additional files 8: Table S6).
The advantages and disadvantages of the present genetic linkage map
QTL mapping is a reliable way to resolve the genetic basis of GPC, and a high-density map will increase the accuracy of QTL detection . In the present study, a high-density map comprised of 1473 SNP and 21 SSR, and spanned 2354.48 cM in length. Notably, our genetic map has a good collinearity with the barley reference genome, which is suitable for the identification of QTLs . However, several chromosome intervals were inconsistent with the reference genome sequence. This could be partially explained by the following reasons: 1) a suppressed recombination frequency at the centromere region, 2) the presence of partially homologous sequences or duplication, and 3) the deficiency of polymorphic markers. In addition, four chromosomes (1H, 3H, 4H and 7H) had high genome coverage (95.40–99.93%), while three (2H, 5H and 6H) showed low genome coverage (63.13–65.51%), which might be caused by the lack of polymorphic markers within the centromeric region. Due to the low recombination frequency in the centromeric region, we speculated that it would not influence the identification of the QTLs.
Compared with two SNP maps reported by Close et al.  and Muñoz-Amatriaín et al. , the whole genome of our map expanded in genetic distance by 41.75 and 111.51%, respectively, with individual chromosome extended by 26.44% to 83.89% and 24.30% to 160.97%, respectively (Additional file 9: Table S7). Missing genotype data of each line and large number of heterozygotes could lead to expanded genetic distance [39, 40]. Consistent with this, similar phenomenon was observed in our SNP genotype data, which could partially contribute to the large whole genetic distance. The casual reason will be an interestingly area to further investigation.
Extensive variation for GPC in barley
Determining the phenotypic variation of GPC in a segregating population is a prerequisite for elucidating its genetic foundation and for breeding barley cultivars with desirable GPC. Extensive variation in GPC in different barley genotypes has been reported previously. For example, analysis of 59 cultivated and 99 Tibetan wild barley accessions showed that the GPC ranged from 8.02% to 13.50% and Tibetan wild barley had much higher GPC than cultivated barley . QTL analysis provides an efficient way to look for associations between the phenotypic variance and the markers segregating in a bi-parental population [41, 42] and has been widely used in dissecting GPC variation in barley populations. However, the lack of parental lines with high GPC in most previous studies may have hindered the detection of possible major QTLs for GPC [10, 43]. In the current study, the rare accession ZGMLEL, with consistently high GPC (20.52–22.88%), and an Australian cultivar, Schooner, with relatively low GPC (16.35–17.20%), were used to construct the mapping population. A relatively high broad sense heritability (80.67%) was found, suggesting that QTLs/genes controlling GPC are less environmentally influenced in the ZGMLEL × Schooner population. Thus, the ZGMLEL × Schooner population is a perfect material for identifying QTLs for GPC. Environmentally stable QTLs detected in this way might be suitable for marker-assisted selection (MAS) in barley breeding, which is anticipated to increase efficiency of the genetic improvement for GPC.
Novel QTLs controlling GPC on chromosomes 4H and 7H
To enhance the GPC of barley, novel genes or QTLs with increased effects are of interest for breeding purposes. In our study, two genomic regions harboring four stable QTLs (QGpc.ZgSc-4H.1, QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3) appeared to be novel (Fig. 3, Additional file 10: Table S8).
A significant QTL QGpc.ZgSc-4H.1 for GPC was identified in the telomeric region of chromosome 4HS and was steadily expressed in three environments. QTLs affecting GPC have been identified on 4HS [8, 43, 44] and 4HL [15, 45, 46]. For example, Marquez-Cedillo et al. identified a QTL for GPC at the region of the intermedium-c (int-c) locus, which is obviously different from QGpc.ZgSc-4H.1 (Fig. 3) . Therefore, QGpc.ZgSc-4H.1 likely represents a new locus for GPC, although its contribution to the variation of GPC was relatively small.
Remarkably, three neighboring QTLs (QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3) were detected on chromosome 7HS. QTLs for GPC on chromosome 7H have been extensively reported [8, 12,13,14]. For example, Emebiri et al.  reported mapping of two QTLs for GPC in the telomeric and centromeric regions that are probably within the physical intervals of 15.8–40.0 and 261.8–277.6 Mb, respectively (Fig. 3) . Marquez-Cedillo et al.  and Abdel-Haleem et al.  identified a consensus QTL near the nud locus on chromosome 7HL [8, 14]. However, the location of these QTLs was different from that detected in the present study. Therefore, these three QTLs in adjacent intervals are likely to be new QTLs, which might be due to the utilization of specific genetic materials in the present study.
Consensus QTL regions for GPC on chromosomes 2H and 6H
An efficient method to introgress favorable alleles into elite germplasm is to select consensus QTLs that steadily affect GPC in different genetic backgrounds and environments . For example, Emebiri reported that pyramiding two consensus QTLs on chromosomes 6HS and 5HS could significantly decrease GPC levels by 4% compared to the commercial check . In this study, two genomic regions on chromosomes 2HL and 6HL for GPC were coincident with QTLs reported in previous studies (Fig. 3; Additional file 10: Table S8). For example, a stable QTL QGpc.ZgSc-2H.1 on chromosome 2HL was coincident with the locus reported by Marquez-Cedillo et al. . Another major QTL, QGpc.ZgSc-6H.3, explaining the highest GPC variance was mapped at a similar locus to Qpro6a detected in the Morex/Steptoe DH population . However, the additive effect of our loci (0.43–0.59%), however, is higher than Qpro6a (0.14%), which might be caused by the special materials used in our study.
To date, two homologous genes, HvNAM1 and HvNAM2 associating with GPC on the short arm of chromosomes 6H and 2H, respectively, have been widely studied [17, 48, 49]. For example, Cai et al.  performed a multi-platform candidate gene-based association analysis using 59 cultivated and 99 Tibetan wild barley genotypes and found that the haplotypes of HvNAM1 and HvNAM2 markers were associated with GPC in barley. In the present study, two identified QTLs, QGpc.ZgSc-6H.3 and QGpc.ZgSc-2H.1 associated with GPC were also detected on chromosomes 6H and 2H, respectively, while they were both located on the long arms, demonstrating that HvNAM1 and HvNAM2 were obviously different from the QTLs detected in this study.
QTLs for GPC linked in coupling phase on chromosomes 6H and 7H
Neighboring QTLs associated with many important traits, such as yield and quality, that are linked in coupling phase are commonly observed in primary QTL analysis [50, 51]. Previous studies have tried to dissect QTLs in coupling phase using nearly isogenic lines (NILs) or residual heterozygous lines (RHLs), and found that coupling QTLs were partially attributed to tightly linked independent QTLs [52,53,54,55]. For example, Han et al. identified two QTLs each for malt extract and for α-amylase and two to three for diastatic power in a complex QTL region using advanced segregation populations . In this study, we detected two genomic regions on chromosomes 6HL and 7HS, each of which harbored linked QTLs for GPC. Region 6H contains three neighboring QTLs, i.e., one environmentally stable QTL (QGpc.ZgSc-6H.3) and two putative QTLs (QGpc.ZgSc-6H.2 and QGpc.ZgSc-6H.4). These three QTLs with favorable alleles from one parent (ZGMLEL) were in coupling phase. A shadow QTL, significant but false, is caused by a real QTL in an adjacent marker interval . Since QGpc.ZgSc-6H.2 and QGpc.ZgSc-6H.4 were located close to the stable QTL QGpc.ZgSc-6H.3, it is difficult to determine whether these two loci were shadow or genuine QTLs. Similarly, region 7H also contains three linked QTLs that were in coupling phase. Unlike the region 6H, region 7H harbored three environmentally stable QTLs and showed similar effects on GPC. Further studies are needed to dissect these two complex regions using advanced population.
The contribution of stable QTLs to GPC
QTL effect was generally not precisely estimated in primary QTL analysis due to the genetic noise in mapping populations [56,57,58,59,60,61,62]. In view of this point, NILs were proposed and developed as an ideal population for QTL validation, especially for the QTL with a minor effect [63, 64]. In the present study, three genetic intervals harboring five stable QTLs, QGpc.ZgSc-2H.1, QGpc.ZgSc-6H.3, QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3 were identified and had additive effects of 0.23, 0.46, 0.38, 0.41, and 0.39% in the RIL population, respectively. Of these loci, three located on chromosome 7HS were linked together and resided in one genomic region. Three NIL populations, BC3F2-I, BC3F2-II and BC3F2-III, were developed according to a standard process of consecutive backcross, which targeted the genomic regions of chromosomes 2H, 6H and 7H, respectively. These stable QTLs were validated in the corresponding NIL population and their contribution to GPC could be directly compared between two homozygous groups in a similar genetic background. In populations I, II and III, the average GPC of plants carrying homozygous ZGMLEL were 0.66, 0.99 and 0.71% higher than that of plants carrying homozygous Schooner, respectively, providing further evidence for the reliability of these stable QTLs. Interestingly, the locus on chromosome 2HL (QGpc.ZgSc-2H.1) exhibited a strong potential increased in GPC in the Schooner background, which illustrated the conclusion that near-isogenic lines could be used to identify a quantitative locus even though it showed a relatively small effect on the phenotype . Thus, the minor QTL, QGpc.ZgSc-2H.1, is feasible for cloning using NILs with least genetic noise. The estimation of the combination of three loci (QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3) in BC3F2-III, however, was somewhat lower than expected by the sum of the individual effects of three loci, which might be caused by QTL × QTL interactions, QTL × environment interactions or QTL × new genetic background interactions.
A relatively lower GPC of BC3F2 populations (11.66–16.82%) grown in Beijing during year 2016 was observed as compared with the RIL population grown in Beijing during year 2013 (16.08–23.36%) or 2014 (14.08–23.75%). This provided further evidence that GPC is largely modified by environmental conditions, which may be due to the alteration of weather condition in 2016. Previous studies indicated high temperature during grain filling stage could result in enhanced GPC [66, 67]. Interestingly, the number of days after flowering with a maximum temperature above 30 °C in Beijing during 2013 and 2014 were 21 and 23 days, respectively, which is obviously more than that of 2016 (17 days) (Additional file 1: Table S1). Collectively, we speculated that the lower of GPC in 2016 may could be partially attributed to the alteration of temperature as compared to the other years.
Potential application of stable QTL for MAS in barley breeding
Since grain protein concentration is greatly influenced by environmental factors, breeding high-GPC cultivars only through phenotypic evaluation has been proved to be less effective [68, 69]. Hence, selection of genomic regions containing QTLs that could express steadily under multiple environments is an efficient way to cultivate barley varieties . Here, we report mapping of six environmentally stable QTLs for GPC that might be useful during barley breeding. Furthermore, we verified the effect of five stable QTLs located on chromosomes 2HL (QGpc.ZgSc-2H.1), 6HL (QGpc.ZgSc-6H.3) and 7HS (QGpc.ZgSc-7H.1, QGpc.ZgSc-7H.2, and QGpc.ZgSc-7H.3) using three NIL populations. In many cases, improvement of GPC is always accompanied by a significant reduction in grain yield [20, 71]. Notably, our preliminary data revealed that no significant difference was found for TKW and GY between two different homozygous groups for each QTL region. However, to further investigate the effect of QTL for GPC on yield, it is necessary to carry out experiment using lines, instead of single-plant strategy. Taken together, the identification of SSR marker intervals flanking these stable QTLs on chromosomes 2H (2L10-2L12), 6H (6L89-6L147) and 7H (7S87-7S40) may provide favorable regions for marker-assisted introgression into the elite barley germplasm.
Based on genotyping 190 RILs in a genome-wide scale and measuring of GPC collected from six environments, six environmentally stable QTLs were significant associated with GPC, among which four QTLs on chromosomes 4H and 7H were firstly identified in the present study. Furthermore, three genomic regions harboring five stable QTLs on chromosomes 2H, 6H and 7H were validated using NIL populations, suggesting the reliability of QTLs detected in primary population. The markers linked to the stable QTLs would be valuable for MAS in barley breeding.
Grain protein concentration
Quantitative trait locus
Recombinant inbred line
Single nucleotide polymorphism
Simple sequence repeat
Thousand grain weight
Friedman M. Nutritional value of proteins from different food sources. A review J Agr Food Chem. 1996;44(1):6–29.
Young VR, Pellett PL. Plant proteins in relation to human protein and amino acid nutrition. Am J Clin Nutr. 1994;59(5):1203S–12S.
Shewry PR, Halford NG. Cereal seed storage proteins: structures, properties and role in grain utilization. J Exp Bot. 2002;53(370):947–58.
Baik BK, Ullrich SE. Barley for food: characteristics, improvement, and renewed interest. J Cereal Sci. 2008;48(2):233–42.
Piper TE, Rasmusson DC. Selection for low protein in barley. Crop Sci. 1984;24(5):853–4.
Zeng ZB. Precision mapping of quantitative trait loci. Genetics. 1994;136(4):1457–68.
Goblirsch CA, Horsley RD, Schwarz PB. A strategy to breed low-protein barley with acceptable kernel color and diastatic power. Crop Sci. 1996;36(1):41–4.
Marquez-Cedillo L, Hayes P, Jones B, Kleinhofs A, Legge W, Rossnagel B, Sato K, et al. QTL analysis of malting quality in barley based on the doubled-haploid progeny of two elite north American varieties representing different germplasm groups. Theor Appl Genet. 2000;101:173–84.
Marquez-Cedillo LA, Hayes PM, Kleinhofs A, Legge WG, Rossnagel BG, Sato K, Ullrich SE, Wesenberg DM. QTL analysis of agronomic traits in barley based on the doubled haploid progeny of two elite north American varieties representing different germplasm groups. Theor Appl Genet. 2001;103(4):625–37.
See D, Kanazin V, Kephart K, Blake T. Mapping genes controlling variation in barley grain protein concentration. Crop Sci. 2002;42(3):680–5.
Mickelson S, See D, Meyer FD, Garner JP, Foster CR, Blake TK, Fischer AM. Mapping of QTL associated with nitrogen storage and remobilization in barley (Hordeum vulgare L.) leaves. J Exp Bot. 2003;54(383):801–12.
Emebiri LC, Moody DB, Panozzo JF, Chalmers KJ, Kretschmer JM, Ablett GA. Identification of QTLs associated with variations in grain protein concentration in two-row barley. Crop Pasture Sci. 2003;54:1211–21.
Emebiri LC, Moody DB, Horsley R, Panozzo J, Read BJ. The genetic control of grain protein content variation in a doubled haploid population derived from a cross between Australian and north American two-rowed barley lines. J Cereal Sci. 2005;41(1):107–14.
Abdel-Haleem H, Bowman J, Giroux M, Kanazin V, Talbert H, Surber L, Blake T. Quantitative trait loci of acid detergent fiber and grain chemical composition in hulled × hull-less barley population. Euphytica. 2010;172(3):405–18.
Abdel-Haleem H, Bowman J, Surber L, Blake T. Variation in feed quality traits for beef cattle in Steptoe × Morex barley population. Mol Breeding. 2012;29(2):503–14.
Moralejo M, Swanston JS, Muñoz P, Prada D, Elía M, Russell JR, Ramsay L, Cistue L, Codesal P, Casas AM, Romagosa I, Powell W, Molina-Cano JL. Use of new EST markers to elucidate the genetic differences in grain protein content between European and north American two-rowed malting barleys. Theor Appl Genet. 2004;110(1):116–25.
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.
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.
Pauli D, Brown-Guedira G, Blake TK. Identification of malting quality QTLs in advanced generation breeding germplasm. J Am Soc Brew Chem. 2015;73(1):29–40.
Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J. A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. 2006;314(5803):1298–301.
Moody DB, Emebiri LC, Spackman M, Mekuria GT. Breeding for qualitative and quantitative traits using markers. In: Proceedings of a joint meeting of the 11th Australian barley technical symposium and the 53rd Australian cereal chem conference, Glenelg, SA; 2003.
Emebiri LC. Breeding malting barley for consistently low grain protein to sustain production against predicted changes from global warming. Mol Breeding. 2015;35(1):1–10.
Stange M, Utz HF, Schrag TA, Melchinger AE, Würschum T. High-density genotyping: an overkill for QTL mapping? Lessons learned from a case study in maize and simulations. Theor Appl Genet. 2013;126(10):2563–74.
Close TJ, Bhat PR, Lonardi S, Wu Y, Rostoks N, Ramsay L, Druka A, Stein N, Svensson JT, Wanamaker S, Bozdag S, Roose ML, Moscou MJ, Chao S, Varshney R, Szucs P, Sato K, Hayes PM, Matthews DE, Kleinhofs A, Muehlbauer GJ, DeYoung J, Marshall DF, Madishetty K, Fenton RD, Condamine P, Graner A, Waugh R. Development and implementation of high-throughput SNP genotyping in barley. BMC Genomics. 2009;10:582.
Muñoz-Amatriaín M, Cuesta-Marcos A, Endelman JB, Comadran J, Bonman JM, Bockelman HE, Chao S, Russell J, Waugh R, Hayes PM, Muehlbauer GJ. The USDA barley core collection: genetic diversity, population structure, and potential for genome-wide association studies. PLoS One. 2014;9(4):e94688.
Muñoz-Amatriaín M, Moscou MJ, Bhat PR, Svensson JT, Bartoš J, Suchánková P, Šimková H, Endo TR, Fenton RD, Lonardi S, Castillo AM, Chao S, Cistué L, Cuesta-Marcos A, Forrest KL, Hayden MJ, Hayes PM, Horsley RD, Makoto K, Moody D, Sato K, Vallés MP, Wulff BBH, Muehlbauer GJ, Doležel J, Close TJ. An improved consensus linkage map of barley based on flow-sorted chromosomes and single nucleotide polymorphism markers. Plant Genome. 2011;4(3):238–49.
Comadran J, Kilian B, Russell J, Ramsay L, Stein N, Ganal M, Shaw P, Bayer M, Thomas W, Marshall D, Hedley P, Tondelli A, Pecchioni N, Francia E, Korzun V, Walther A, Waugh R. Natural variation in a homolog of antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat Gent. 2012;4(12):1388–92.
Walker CK, Ford R, Muñoz-Amatriaín M, Panozzo JF. The detection of QTLs in barley associated with endosperm hardness, grain density, grain size and malting quality using rapid phenotyping tools. Theor Appl Genet. 2013;126(10):2533–51.
Gong X, Wheeler R, Bovill WD, McDonald GK. QTL mapping of grain yield and phosphorus efficiency in barley in a Mediterranean-like environment. Theor Appl Genet. 2016;129(9):1657–72.
Guo BJ, Luan H, Lin S, Lv C, Zhang X, Xu R. Comparative proteomic analysis of two barley cultivars (Hordeum vulgare L.) with contrasting grain protein content. Front. Plant Sci. 2016;7:542.
Comadran J, Kilian B, Russell J, Ramsay L, Stein N, Ganal M, Shaw P, Bayer M, Thomas W, Marshall D, Hedley P, Tondelli A, Pecchioni N, Francia E, Korzun V, Walther A, Waugh R. Natural variation in a homolog of antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat Gent. 2012;44(12):1388–92.
Van Os H, Stam P, Visser RG, Van Eck HJ. RECORD: a novel method for ordering loci on a genetic linkage map. Theor Appl Genet. 2005;112(1):30–40.
Van Ooijen JW. JoinMap 4®, software for the calculation of genetic linkage maps in experimental populations. Kyazma BV, Wageningen, Netherlands. 2006.
Blake VC, Blake TK (2008) Barley QTL: a compilation of known loci and the discovery of new QTL for digestibility. Proc. 10th Intl. Barley Genetics Symp. April 5–10, 2008. Alexandria, Egypt.
Mayer KF, Martis M, Hedley PE, Šimková H, Liu H, Morris JA, et al. Unlocking the barley genome by chromosomal and comparative genomics. Plant Cell. 2011;23(4):1249–63.
The International Barley Genome Sequencing Consortium. A physical, genetic and functional sequence assembly of the barley genome. Nature. 2012;491:711–6.
Groos C, Robert N, Bervas E, et al. Genetic analysis of grain protein-content, grain yield and thousand-kernel weight in bread wheat. Theor Appl Genet. 2003;106(6):1032–40.
Wang X, Yu K, Li H, et al. High-density SNP map construction and QTL identification for the apetalous character in Brassica napus L. Front Plant Sci. 2015;6:1164.
Shifman S, Bell JT, Copley RR, Taylor MS. A high-resolution single nucleotide polymorphism genetic map of the mouse genome. PLoS Biol. 2006;4(12):e395.
Kong A, Gudbjartsson DF, Sainz J, et al. A high-resolution recombination map of the human genome. Nat Genet. 2002;31(3):241–7.
Kearsey MJ, Farquhar AG. QTL analysis in plants; where are we now? Heredity. 1998;80(2):137–42.
Andersson L. Genetic dissection of phenotypic diversity in farm animals. Nat Rev Genet. 2001;2(2):130–8.
Canci PC, Nduulu LM, Dill-Macky R, Muehlbauer GJ, Rasmusson DC, Smith KP. Genetic relationship between kernel discoloration and grain protein concentration in barley. Crop Sci. 2003;43(5):1671–9.
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(1):43–51.
Mather DE, Tinker NA, LaBerge DE, Edney M, Jones BL, Rossnagel BG, Legge WG, Briggs KG, Irvine RG, Falk DE, Kasha KJ. Regions of the genome that affect grain and malt quality in a north American two-row barley cross. Crop Sci. 1997;37(2):544–54.
Hayes PM, Liu BH, Knapp SJ, Chen F, Jones B, Blake T, Franckowiak J, Rasmusson D, Sorrells M, Ullrich SE, Wesenberg D, Wesenberg D. Quantitative trait locus effects and environmental interaction in a sample of north American barley germplasm. Theor Appl Genet. 1993;87(3):392–401.
Hao ZF, Li XH, Xie CX, Li MS, Zhang DG, Bai L, Zhang SH. Two consensus quantitative trait loci clusters controlling anthesis-silking interval, ear setting and grain yield might be related with drought tolerance in maize. Ann Appl Biol. 2008;153(1):73–83.
Jukanti AK, Fischer AM. A high-grain protein content locus on barley (Hordeum vulgare) chromosome 6 is associated with increased flag leaf proteolysis and nitrogen remobilization. Physiol Plantarum. 2008;132(4):426–39.
Distelfeld A, Korol A, Dubcovsky J, et al. Colinearity between the barley grain protein content (GPC) QTL on chromosome arm 6HS and the wheat Gpc-B1 region. Mol Breeding. 2008;22(1):25–38.
Zhao JY, Becker HC, Zhang DQ, Zhang YF, Ecke W. Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield. Theor Appl Genet. 2006;113(1):33–8.
Wu B, Mao DH, Liu TM, Li ZX, Xing YZ. Two quantitative trait loci for grain yield and plant height on chromosome 3 are tightly linked in coupling phase in rice. Mol Breeding. 2015;35(8):1–13.
Han F, Ullrich SE, Kleinhofs A, Jones BL, Hayes PM, Wesenberg DM. Fine structure mapping of the barley chromosome-1 centromere region containing malting-quality QTLs. Theor Appl Genet. 1997;95(5):903–10.
Cai H, Morishima H. QTL clusters reflect character associations in wild and cultivated rice. Theor Appl Genet. 2002;104(8):1217–28.
Du JH, Fan YY, Wu JR, Zhuang JY. Dissection of QTLs for yield traits on the short arm of rice chromosome 6. Agr Sci China. 2008;7(5):513–20.
Guo L, Wang K, Chen JY, Huang DR, Fan YY, Zhuang JY. Dissection of two quantitative trait loci for grain weight linked in repulsion on the long arm of chromosome 1 of rice (Oryza sativa L.). The Crop Journal. 2013;1(1):70–6.
Tuyen DD, Zhang HM, Xu DH. Validation and high-resolution mapping of a major quantitative trait locus for alkaline salt tolerance in soybean using residual heterozygous line. Mol Breeding. 2013;31(1):79–86.
Hamwieh A, Tuyen DD, Cong H, Benitez ER, Takahashi R, Xu DH. Identification and validation of a major QTL for salt tolerance in soybean. Euphytica. 2011;179(3):451–9.
Landi P, Sanguineti MC, Salvi S, Giuliani S, Bellotti M, Maccaferri M, Conti S, Tuberosa R. Validation and characterization of a major QTL affecting leaf ABA concentration in maize. Mol Breeding. 2005;15(3):291–303.
Paterson AH, DeVerna JW, Lanini B, Tanksley SD. Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics. 1990;124(3):735–42.
Alpert KB, Tanksley SD. High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: a major fruit weight quantitative trait locus in tomato. PNAS. 1996;93(26):15503–7.
Wang GW, He YQ, Xu CG, Zhang Q. Identification and confirmation of three neutral alleles conferring wide compatibility in inter-subspecific hybrids of rice (Oryza sativa L.) using near-isogenic lines. Theor Appl Genet. 2005;111(4):702–10.
Chu ZH, Fu BY, Yang H, Xu CG, ZK L, Sanchez A, Park YJ, Bennetzen JL, Zhang QF, Wang SP. Targeting xa13, a recessive gene for bacterial blight resistance in rice. Theor Appl Genet. 2005;112(3):455–61.
Zhang Y, Luo L, Liu T, Xu C, Xing Y. Four rice QTL controlling number of spikelets per panicle expressed the characteristics of single Mendelian gene in near isogenic backgrounds. Theor Appl Genet. 2009;118(6):1035–44.
Pumphrey MO, Bernardo R, Anderson JA. Validating the QTL for Fusarium head blight resistance in near-isogenic wheat lines developed from breeding populations. Crop Sci. 2007;47(1):200–6.
Takahashi Y, Shomura A, Sasaki T, Yano M. Hd6, a rice quantitative trait locus involved in photoperiod sensitivity, encodes the α subunit of protein kinase CK2. PNAS. 2001;98(14):7922–7.
Lin CJ, Li CY, Lin SK, et al. Influence of high temperature during grain filling on the accumulation of storage proteins and grain quality in rice (Oryza sativa L.). J Agr Food Chem. 2010;58(19):10545–52.
Ingvordsen CH, Gislum R, Jørgensen JR, et al. Grain protein concentration and harvestable protein under future climate conditions. A study of 108 spring barley accessions. J Exp Bot. 2016;67:2151–58.
Shewry PR. Improving the protein content and composition of cereal grain. J Cereal Sci. 2007;46(3):239–50.
Uauy C, Brevis JC, Chen X, Khan I, Jackson L, Chicaiza O, et al. High-temperature adult-plant (HTAP) stripe rust resistance gene Yr36 from Triticum turgidum ssp. dicoccoides is closely linked to the grain protein content locus Gpc-B1. Theor Appl Genet. 2005;112(1):97–105.
Zhou H, Liu S, Liu Y, Liu Y, You J, Deng M, Ma J, Chen G, Wei Y, Liu C, Zheng Y. Mapping and validation of major quantitative trait loci for kernel length in wild barley (Hordeum vulgare ssp. spontaneum). BMC Genet. 2016;17(1):130.
Löffler CM, Busch RH. Selection for grain protein, grain yield, and nitrogen partitioning efficiency in hard red spring wheat. Crop Sci. 1982;22(3):591–5.
Ramsay L, Comadran J, Druka A, et al. INTERMEDIUM-C, a modifier of lateral spikelet fertility in barley, is an ortholog of the maize domestication gene TEOSINTE BRANCHED 1. Nat Gent. 2011;43:169–72.
This work was financially supported by the National Natural Science Foundation of China (Grant No. 31271710 and No. 31260326).
Availability of data and materials
All data generated or analysed during this study are included in the supplementary files as well as in the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Location-year information and climate data for field trails. (DOC 39 kb)
Phenotypic data of two parental lines and 190 RILs under six individual environments and the combined analysis. (XLSX 31 kb)
Map locations of all markers and genotypic data of 190 RILs. (XLSX 1050 kb)
Illustration of nine linkage groups constructed using 190 RILs. Notes: A centiMorgan (cM) scale is shown on the left. Black ellipses represent the approximate position for chromosome centromeres. The marker names are not shown. Detail information of the linkage group is provided in Additional file 3. (DOC 159 kb)
Alignments of SNP markers to barley reference sequence. (XLSX 76 kb)
New genetic linkage maps in the target region of the stable QTLs detected in ZGMLEL × Schooner RIL population. (A) chromosome 2HL, (B) 6HL, and (C) 7HS. New SSR markers are showed in bold and underlined. (DOC 299 kb)
Markers developed to select the target regions of stable QTLs on chromosomes 2H, 6H and 7H. (DOC 38 kb)
Phenotypic data of two homozygous groups in each NIL population. (XLSX 36 kb)
Genomic regions harboring environmentally stable QTLs for grain protein concentration (GPC) in the ZGMLEL × Schooner population. Notes: QTLs in bold represent the environmentally stable QTLs. (DOC 35 kb)
About this article
Cite this article
Fan, C., Zhai, H., Wang, H. et al. Identification of QTLs controlling grain protein concentration using a high-density SNP and SSR linkage map in barley (Hordeum vulgare L.). BMC Plant Biol 17, 122 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-017-1067-6
- Barley (Hordeum vulgare L.)
- Grain protein concentration (GPC)
- SNP marker
- Near isogenic line