- Research article
- Open Access
Detection of QTLs for seedling characteristics in barley (Hordeum vulgare L.) grown under hydroponic culture condition
BMC Genetics volume 18, Article number: 94 (2017)
Seedling characteristics play significant roles in the growth and development of barley (Hordeum vulgare L.), including stable stand establishment, water and nutrients uptake, biotic resistance and abiotic stresses, and can influence yield and quality. However, the genetic mechanisms underlying seedling characteristics in barley are largely unknown and little research has been done. In the present work, 21 seedling-related characteristics are assessed in a barley double haploid (DH) population, grown under hydroponic conditions. Of them, leaf age (LAG), shoot height (SH), maximum root length (MRL), main root number (MRN) and seedling fresh weight (SFW) were investigated at the 13th, 20th, 27th, and 34th day after germination. The objectives were to identify quantitative trait loci (QTLs) underlying these seedling characteristics using a high-density linkage map and to reveal the QTL expression pattern by comparing the QTLs among four different seedling growth stages.
A total of 70 QTLs were distributed over all chromosomes except 4H, and, individually, accounted for 5.01%–77.78% of phenotypic variation. Out of the 70 detected QTLs, 23 showed a major effect on 14 seedling-related characteristics. Ten co-localized chromosomal regions on 2H (five regions), 3H (two regions) and 7H (three regions) involved 39 QTLs (55.71%), each simultaneously influenced more than one trait. Meanwhile, 9 co-localized genomic regions involving 22 QTLs for five seedling characteristics (LAG, SH, MRL, MRN and SFW) at the 13th, 20th, 27th and 34th day-old seedling were common for two or more growth stages of seedling. QTL in the vicinity of Vrs1 locus on chromosome 2H with the favorable alleles from Huadamai 6 was found to have the largest main effects on multiple seedling-related traits.
Six QTL cluster regions associated with 16 seedling-related characteristics were observed on chromosome 2H, 3H and 7H. The majority of the 29 regions identified for five seedling characteristics were selectively expressed at different developmental stages. The genetic effects of 9 consecutive expression regions displayed different developmental influences at different developmental stages. These findings enhanced our understanding of a genetic basis underlying seedling characteristics in barley. Some QTLs detected here could be used for marker-assisted selection (MAS) in barley breeding.
As the earliest domesticated crop and the fourth most important cereal crop in the world, barley (Hordeum vulgare L.) is not only widely used as human food and animal feed, but is also an ideal model of genetic research because it is early maturing, diploid, self- fertilizing and has a short growth period [1,2,3]. However, its growth and production are also greatly affected by seedling characteristics. The seedling stage (from seedling emergence to jointing) is considered as the most critical stage for barley production, and plays significant role in barley growth and development, including stable stand establishment, water and nutrients uptake, biotic resistance and abiotic stresses, and can influence yield and quality [4, 5]. In fact, strong seedling vigor and rapid seedling growth are important breeding targets in barley, as well as in other crops [6,7,8]. Seedling vigor has been identified to be associated with salt tolerance , drought tolerance [10, 11], nutrient uptake , weed competitiveness  and yield [13, 14]. It is difficult to assess seedling characteristics because a majority of important seedling characteristics are complex, quantitative traits determined by an array of developmental processes, genetic and environmental factors. The genetic mechanisms controlling seedling characteristics in barley are poorly known and little research has been done in this regard. Genetic variation for seedling characteristics in barley varieties does exist, and dissecting the genetic and molecular basis of seedling characteristics is necessary for genetic improvement of barley cultivars to enhance seedling vigor.
The information on the genetic basis underlying seedling characteristics, especially for root characteristics, is limited because it is difficult to obtain reliable phenotypic data in a large number of seedlings and in a complex external environment. It is also difficult to continuously measure these traits from the same seedling due to destruction of the seedlings in field experiments [15,16,17]. Therefore, an alternative approach is to examine seedling characteristics under controlled conditions using hydroponic culture, which has several advantages over field and other conventional techniques, including: (i) ease of investigation of seedling characteristics (as compared with field experiment), (ii) exclusion of soil and environment interference to increase measurement repeatability, (iii) investigation of large numbers of lines in short period of time, and (iv) precise control of the concentration of nutrient concentration . At present, this method has been widely used in rice [19,20,21], maize [15, 22, 23] and wheat [16, 24, 25], and many QTLs have been detected at various water and nutrient regimes. Some researchers have used this method to identify QTLs associated with salt tolerance [3, 26], waterlogging tolerance [27,28,29] and nitrogen stress tolerance  in barley.
With the development of molecular markers and continuous encryption of genetic map, QTL analysis has become a powerful tool to dissect complex traits and identify chromosomal regions harboring genes that control these quantitative traits [31, 32], and has been extensively used for genetic dissection of important traits in barley [33,34,35]. The QTLs underlying yield [34, 36, 37] and quality [35, 38, 39] traits at later growth stages have been characterized previously, and the QTLs for seedling characteristics reported were mainly identified under salt tolerance [29, 40, 41], waterlogging tolerance [29, 42], drought tolerance  and nitrogen stress tolerance , while seedling characteristics related to grain yield and quality at early developmental stages are not well investigated.
The objectives of our study were to identify QTLs underlying seedling characteristics in hydroponics at seedling growth stages, and compare its QTL expression patterns among four different seedling growth stages.
Evaluation of seedling characteristics
The DH population and two parents were grown hydroponically in a greenhouse as shown in Fig. 1a. Figure 1b showed phenotypic difference in the seedling of the Huadama 6 and Huaai 11. Huadama 6 and Huaai 11 have significant differences in seedling characteristics. All traits values, except FOLW, in Huadamai 6 were higher than those in Huaai 11. The T-test showed that two parents were significantly different on all traits (p < 0.05) but FILW and FOLW (Table 1).
The 21 seedling characteristics in DH population were also significantly different, with a highly phenotypic variation. The variation coefficient ranged from 6.88%–35.29%. Heritability ranged from 89.86% to 97.05% (Table 1). Frequency distribution histogram of the 21 seedling characteristics were shown in Additional file 1: Figure S1. Shapiro-Wilk test showed that some of the characteristics fitted the normal distribution, including LAG (S2), SH (S1-S4), MRL (S1), SFW (S1), LL (except FILL), LW and LA. In addition, transgressive segregation was observed for all 21 seedling characteristics.
Detection of QTLs
A total of 70 QTLs on all barley chromosomes (except 4H) were detected in 17 of the 21 seedling characteristics using ICIM mapping (Table 2; Fig. 2). No QTL was identified for SELL, FILW, THLW and FOLW. Of these, a single QTL contributed between 5.01% and 77.78% to the phenotypic variance, with LOD value ranging from 3.04 to 27.07 (Table 2; Additional file 2: Figure S2). Most QTLs primarily dispersed on chromosome 2H (35 QTLs), 3H (11 QTLs) and 7H (19 QTLs) (Figs. 2 and 3a). Ten co-localized chromosomal regions affected more than one seedling characteristic distributed on chromosome 2H, 3H and 7H (Table 3).
Dynamic QTL mapping on five seedling traits (LAG, SH, SFW, MRL and MRN) at four developmental stages identified 8 QTLs for LAG, 14 QTLs for SH, 8 QTLs for SFW, 8 QTLs for MRL and 4 QTLs for MRN (Table 2; Fig. 2). Phenotypic variation explained by single QTL ranged from 6.03% to 73.86% for LAG, from 5.49% to 24.90% for SH, from 7.85% to 32.39% for SFW, from 8.08% to 16.23% for MRL and from 7.71% to 31.28% for MRN. Interestingly, we found that only one QTL for SH was detected at all four stages, 8 QTLs for LAG, SH, SFW and MRL were detected for two or three stages, but also some QTLs were specific for either S1, S2, S3 and S4 (Table 4; Fig. 4). In addition, we also observed that some QTLs, such as the QTL at 123 cM on chromosome 2H, not only expressed at all stages for SH but also affected the SFW at the S2 stage. Similar QTLs were identified at 127/129 cM and 137/139 cM on chromosome 2H and at 30/32 cM and 80/81 cM on chromosome 3H. Two major QTLs for LAG (qLAG2–228) and SFW (qSFW2–202) explaining more than 20% of the phenotypic variation at more than one stage of seedling growth was identified on chromosome 2H. We analyzed the temporal-dynamics contribution of QTLs detected at any of the four stages to the proportion of explained phenotypic variation for five seedling characteristics (Fig. 5). The results showed that the proportion of explained phenotypic variance of some QTLs remained on a stable level at all stages of seedling growth; whereas, others substantially changed among the four stages of seedling growth.
We assessed four seedling characteristics (SLFW, RFW, SLDW and RDW) at S4 stage, and detected one QTL for SLFW, 2 QTLs for RFW, 3 QTLs for SLDW and 2 QTLs for RDW (Table 2; Fig. 2). A single QTL identified for these four seedling characteristics could individually explain 28.83% (SLFW), 9.81%–32.06% (RFW), 8.17%–23.62% (SLDW) and 10.58%–20.22% (RDW) of the phenotypic variation. These QTLs were almost invariably located on chromosome 2H except for QTL qSLDW7–45. The major QTL qSLFW2–202 was detected for SLFW on chromosome 2H around 127/129 cM and co-located with QTLs for SFW at S1, S3 and S4 (qSFW2–202), MRN at S2 (qMRN2–199) and SLDW at S4 (qSLDW2–202). Another major QTL was identified on chromosome 2H around 137/139 cM and had effect on MRN at S3 and S4 (qMRL2–214), RFW at S4 (qRFW2–209) and RDW at S4 (qRDW2–214).
For eight seedling characteristics (FILL, THLL, FOLL, SELW, FILA, SELA, THLA and FOLA), a total of 21 QTLs on chromosome 2H (nine QTLs), chromosome 3H (three QTLs), chromosome 6H (one QTL) and on 7H (eight QTLs) were identified. The LOD value ranged from 3.17 to 27.07, and individual QTL accounted for 5.01%–77.78% of the phenotypic variation (Table 2; Additional file 2: Figure S2). Three, two, two, two, one, four, five, one and two QTLs were detected for FILL, THLL, FOLL, SELW, FILA, SELA, THLA and FOLA, respectively. Interestingly, we found that six major QTLs, qFILL7–154, qTHLL7–145, qFOLL7–237, qFILA7–134, qSELA7–142 and qTHLA7–145, distributed on chromosome 7H, had positive alleles from Huadamai 6. In addition, we also observed a number of QTLs that affected multiple traits, such as 127/129 cM on chromosome 2H.
Seedling characteristics include stem-and-leaf and root systems. A good stem-and-leaf system in the seedling stage is the basis for seedlings growth potential of robust, strong plants that are resistant to stress. To sustain a strong seedling vigor and rapid seedling growth, seedling must produce abundant nutrition. The amount of nutrition production mainly depends on the first few leaves of the seedling . In addition, growth and development of healthy plants mainly rely on the strong root system (long MRL, high RDW and enough MRN), which were closely related to final grain yield and quality . Some studies indicated that MRL and MRN were one of the most desirable traits contributing to drought avoidance, and RDW was another important parameter for the roots associated with water use efficiency and long-term drought [46, 47]. Recent studies have also confirmed that most of the seedling traits, such as leaf area (LA), seedling height (SH), seedling fresh weight (FW), seedling dry weight (DW), shoot dry weight (SDW) and root dry weight (RDW), were closely and positively correlated with yield component traits . Previous genetic research on barley were focused mainly on the later stage; thus, the genetic mechanism underlying seedling characteristics is still lacking. In this study, we assessed QTLs underlying 21 seedling characteristics (including five seedling characteristics assessed at four stages of seedling growth) under hydroponic culture conditions and compared the QTL expression pattern among four different seedling growth stages.
Advantages of hydroponic culture
For larger populations, it is difficult to assess seedling characteristics, especially whole-root systems, and for multi-stage barley planted in soil. Hydroponics culture allowed us to assess the same seedling over multi-stages. Genes or QTLs expression under hydroponic conditions might be different from field trials. However, the expression of genes or QTLs is always affected by the environment, while complex environments and uneven fertilizing in the field trial will affect the genetic mechanism of seedling characteristics, and may not reflect the intrinsic genetic program of barley seedling characteristics . Thus, stable and consistent environmental control is needed to assess seedling characteristics. Recent studies have shown that QTLs identified under hydroponics conditions corresponded to QTLs detected in field trials, and suggested that a hydroponic system is a fast and cost-effective method for early QTL detection and marker-assisted allelic selection . Consequently, we established a hydroponic system with homogeneous-growth conditions to evaluate barley-seedling characteristics. Temperature, humidity and lighting conditions were controlled by the automatic control system, and the nutrient solution was renewed every seventh day to ensure adequate nutrition. QTLs identified in this environment could reflect the intrinsic genetic mechanisms underlying barley seedling characteristics.
QTL for seedling characteristics
With the same DH population used in this study, Ren et al. [49,50,51,52], Liu et al.  and Wang et al.  identified certain QTLs for physiological and morphological traits of flag leaf and some agronomic and quality traits (Additional file 3: Figure S3). Similar to their study, most of the QTLs identified in this study were also located on chromosome 2H and 7H. We found that some QTLs for seedling characteristics were co-located with QTLs for yield or yield-related traits identified in previous studies. For example, a significant QTL qSH2–191 for SH detected in all stages of seedling growth was close to the SNP marker 2_514353957, which is likely the same to the yield-related QTL qSms2–7 and qTgw2–1 reported by Wang et al.  (Table 2; Additional file 3: Figure S3). In addition, the co-localized QTLs for THLL, SELA and THLA were detected on chromosome 7H and are likely the same to qRWC7–9 and qSPD7–9 for physiological and morphological traits of flag leaf at the pre-filling stage reported by Liu et al.  (Table 2; Additional file 3: Figure S3). These results suggested that seedling characteristics and yield are related to a certain extent.
In previous studies, QTLs conferring seedling height (SH) [29, 30], root dry weight (RDW) [46, 53,54,55] and stem-and-leaf dry weight (SLDW) [30, 42, 46, 53, 55] were reported on seven linkage groups, while QTLs for seedling fresh weight (SFW), leaf age (LAG) and main root number (MRN) were rarely reported in barley. QTLs for root length were identified on chromosome 1H, 2H, 3H, 4H, 5H and 7H [46, 53,54,55], and QTLs conferring root fresh weight (RFW) were previously reported on chromosome 1H, 2H and 5H [29, 42]. Recently, Hoffmann et al.  analyzed QTL underlying number of leaves longer than 2 cm after 14 days and length of the youngest completely unfolded leaf after 14 days, and found six QTLs for leaf number on chromosome 2H, 4H, 6H and 7H, and nine QTLs for leaf length on all seven chromosomes in wild barley introgression lines (S42ILs). By using GrainGenes3.0 website (http://wheat.pw.usda.gov/GG3/) to compare the genetic markers with those results present in this study, we found that some QTLs identified in our study are likely same to the QTLs reported in the past research. Such as, qSH3–44 for SH is likely the same to the QTL QHei.S42IL-3H reported by Hoffmann et al. . Arifuzzaman et al.  reported that QTL QSdw.S42.2H.c for SLDW on chromosome 2H was linked with the marker bPb-8143, and since the bPb-8143 is near the morphological marker Vrs1, as inferred from GrainGenes3.0 (http://wheat.pw.usda.gov/GG3/), we suggested that a co-localized QTLs for SLDW, SFW, SLFW, MRN, FILA, SELW, SELA and FOLA identified in our study is likely the same to the QTL QSdw.S42.2H.c. This QTL may be an important locus for controlling stem-and-leaf related traits in barley. The SNP marker 2HL_34260490 may be useful for marker-assisted selection (MAS) in barley breeding.
In addition, some new QTLs were also detected in this study. QTL for SH (qSH7–184) on chromosome 7H and QTL for RFW (qRFW2–90 and qRFW2–209) on chromosome 2H, are likely different from the QTL reported by Xu et al. . The QTL for MRL (qMRL2–113 and qMRL3–117), QTL for RDW (qRDW2–130 and qRDW2–214) and QTL for MRN (qMRN3–120) are likely different from those QTLs reported previously [42, 46]. The QTL, qMRL7–27 (S1 and S2), located on chromosome 7H, is different from the QRl.S42IL-7H.a and QRl.S42IL-7H.b, as reported by Hoffmann et al. . The qMRL2–214 (S3 and S4) for MRL, mapped together with qRDW2–214 (S4) for RDW and qRFW2–209 (S4) for RFW on chromosome 2H, are different from those QTLs on chromosome 2H, as reported previously [55, 56]. This region may be important for controlling seedling root traits in barley.
Selective expression of QTLs
According to developmental genetics, different QTLs may have different expression dynamics during trait development [7, 57]. Previous studies focused mainly on late-growth stages, where analysis was limited to the performance of a trait at a fixed time or stage of ontogenesis [33, 34, 58]. In the present study, we assessed five seedling characteristics (LAG, SH, SFW, MRL and MRN) at four stages of seedling growth (13th, 20th, 27th and 34th day-after germination) to reveal the QTL expression pattern.
A total of 42 QTLs involved 29 regions were detected for these five traits at four stages of seedling growth (Table 2; Fig. 4). Although the number of QTL expressed in each stage was approximately the same, only 9 regions were continually expressed in two or more different stages of seedling growth (Table 4; Figs. 3b and 4). Of these, only one of the nine regions was persistently expressed at all stages (Table 4; Figs. 4 and 5). The majority of the 29 regions for these five seedling characteristics were selectively expressed at different developmental stages. This result is in agreement with other studies where different QTLs could be identified at different growth stages, where only a few co-localized QTLs were detected at all stages [7, 59, 60]. The regions that continually expressed at different growth stages might play an important role in the growth and development of these characteristics. However, these regions do not always produce the same effect at each stage but show different developmental influences in different seedling growth stages. For example, two co-localized regions for SH were identified on chromosome 2H and 3H; its contribution to the phenotypic variance undergoes temporal changes. On chromosome 2H, a co-localized region contribution to the phenotypic variance decreased from 18.61% at S1 to 17.74% at S2 and then gradually increased until S4 (22.05%). Differently, a co-localization region on chromosome 3H contribution to the phenotypic variance gradually decreased from 13.49% at S1 to 10.92% at S3 and disappeared at S4. Besides, the co-localized regions for MRL, SFW and LAG contribution to the phenotypic variance also showed similar temporal changes. Generally, our results revealed that the genetic architecture of the five seedling characteristics shows dynamic temporal changes during growth and development.
Pleiotropy or linkage of QTLs
We found here that many QTLs controlling multiple seedling characteristics were located at the same or overlapping marker interval on some chromosomes. This implied an existence of pleiotropic QTL or tightly linked QTL in our study. Xu et al.  suggested that if two QTL peaks are located very close to each other, and the 1-LOD support intervals completely or mostly overlapped, these two QTLs would be regarded as a single QTL having pleiotropic effect. In the present study, we delimited 1.5-LOD support intervals to confirm whether these two QTLs controlling multiple seedling characteristics can be identified as the pleiotropic QTL or tightly linked QTL. Based on this assumption, we found ten regions on chromosomes 2H, 3H and 7H that showed pleiotropic or tightly linkage effects on different seedling characteristics (Table 3).
Five regions on chromosome 2H were observed to simultaneously affect multiple seedling characteristics. Of these, a significant region with pleiotropic or tightly linkage effects on SFW, SLFW, SLDW, MRN, SELW, FILA, SELA and FOLA was close to both the SNP marker, 2HL_34260490 and the morphological marker, Vrs1 (Fig. 3C; Table 3). The Vrs1 locus (controlling row type) has a pleiotropic effect on multiple traits [34, 61,62,63]. For example, Marquez-Cedillo et al.  and Wang et al.  reported that the Vrs1 locus had a pleiotropic effect on thousand grain weight, spike length, plant height, spikelet number on main spike and grain number per plant. Furthermore, some malting quality traits, such as grain-protein-percentage, diastatic power and soluble/total protein ratio, were also coincident with the Vrs1 locus reported by Marquez-Cedillo et al. . Sato et al.  suggested that the fusarium head blight (FHB) resistance QTL in Vrs1 locus reported in two-row × six-row crosses might be pleiotropic effect of Vrs1 locus. Our study also showed a possible pleiotropic effect of the Vrs1 locus and provided evidence to further support previous suggestions.
On chromosome 3H and 7H, we also found that some regions showed pleiotropic effects on different seedling characteristics (Table 3). A pleiotropic QTL or tightly linked QTL, detected between SNP marker 3HL_29995337 and 3HL_32354397 near the SSR marker Bmag13, was related to SH and SELA. Examining the 3H linkage maps of Xue et al.  and GrainGenes3.0 (http://wheat.pw.usda.gov/GG3/) revealed that the marker bPb-6504 is near to the marker Bmag13, indicating that this QTL is likely the same to the qSPL3c associated with spikes per line. In the previous studies using the same DH population, QTL for heading date (Qhd3–13), culm length (Qcl3–13), two internode length(Qitw3–13) and three internode length (Qith3–13), as reported by Ren et al. [49, 51], and QTLs for stomatal conductance (qGs3–13) and transpiration rate (qTr3–13) reported by Liu et al.  were also near this region (Additional file 3: Figure S3). A QTL located at 93/94 cM on chromosome 7H showed that pleiotropic or tightly linkage effects on THLL, SELA and THLA is likely the same to the QTL (qRWC7–9) for relative water content and QTL (qSPD7–9) for relative chlorophyll content, as reported by Liu et al. . These results provided favorable support for pleiotropic QTL or tightly linked QTL identified in our study. Further, this pleiotropic QTL or tightly linked QTL may be used for marker-assisted selection (MAS) and genetic improvement of barley seedling characteristics in barley breeding.
The phenomenon of pleiotropic effect or linked for QTLs universally exists in various crop genomes. Pleiotropic effects of genetic loci are thought to play a critical role in evolution, reflecting functional and developmental relationships among phenotypes . The application of linked QTL and pleiotropic QTL for breeding purposes should be carried out with care, especially if the pleiotropy and linkage are in negative direction . For example, the Huadamai 6 alleles with the pleiotropic QTL or tightly linked QTL on chromosome 3H increased main roots number but decreased maximum root length (Table 3). However, a pleiotropic QTL with opposite effects on the different traits complicated its application in barley improvement.
QTL clusters in the genome
In our study, we found that most of QTLs were located on chromosomes 2H, 3H and 7H (Figs. 2 and 3a). It is noteworthy that many QTLs gathered in the same chromosome regions, even for the unrelated traits (Table 5). For instance, one QTL cluster region for RDW, SLDW, MRL, SFW and RFW was mapped on chromosome 2H close to the marker GBM1218. Liu et al.  reported QTL clusters for physiological and morphological traits of flag leaf (net photosynthetic rate, stomatal conductance, flag leaf area, flag leaf length, flag leaf width, relative chlorophyll content and leaf nitrogen concentration) in the similar region. Similarly, Ren et al.  detected QTLs for heading date traits in this region. Wang et al.  also detected QTLs for spikelet number on main spike in this region. Considering all information here, we suggested that this region may be a credible region for cluster of QTLs.
The most important QTL cluster was detected on chromosome 2H close to the morphological marker Vrs1, which was clustered with 22 QTLs (31.43%) for LAG, SH, MRL, MRN, SFW, SLFW, SLDW, RFW, RDW, FOLL, SELW, FILA, SELA and FOLA with contribution of 12.05%–32.39% (Table 5; Fig. 2). Wang et al.  reported QTL clusters for some agronomic traits (spikelet number on main spike, grain number per spike, spikelet number per plant, grain number per plant, grain weight per plant and thousand grain weight) in this chromosome region. Hori et al.  detected QTLs for plant height, spike exsertion length and thousand-kernel weight in this region. In addition, QTLs for spike number, floret number, grain number and hundred grain weight were also detected in this region by Chutimanitsakun et al. . QTLs for seedling traits in this region may be a result of the pleiotropic effect of the Vrs1 locus or expression of the closely linked genes for Vrs1.
Another significant QTL cluster underlying LAG, FILL, FILA and SELA was detected on chromosome 2H between SNP marker 2HL_18227895 and 2_598509820 (Table 5; Fig. 2). The effects from Huadamai 6 were positive for LAG, but negative for FILL, FILA and SELA. Castro et al.  detected QTLs for photoperiod response in this region. Wang et al.  also detected QTLs for length of main spike, spike density, thousand grain weight and grain weight per spike in this region, and suggested that the SNP marker 2_598509820 in this region could be used for marker-assisted selection. In addition, another three QTL clusters were obtained on chromosome 3H and 7H, all of which enhanced development of stem-and-leaf traits (Table 5; Fig. 2).
QTL clustering in barley was repeatedly reported in some studies [33, 34, 56, 69]. In our study, we demonstrated several significant QTL clusters of seedling characteristics in barley under hydroponic culture. The genetic mechanism for this widespread phenomenon may be a linkage of genes and of pleiotropic effect of a single QTL in the same genomic region [70, 71]. However, whether it is linkage or pleiotropy, further research using fine mapping and cloning of QTLs or genes are needed validation.
In this research, a number of genomic regions containing many seedling-related QTLs involved in different traits and diverse seedling growth stages were detected. The two parents used for mapping population construction are different in row types (six-rowed dwarfing barley cultivar Huaai 11 and two-rowed barley cultivar Huadamai 6). We found that chromosome 2H close to the morphological marker Vrs1 contained the most important QTL cluster regions controlling multiple seedling characteristics. Some QTLs detected in this research could be used as a potential target for marker-assisted selection (MAS) in barley breeding.
Plant materials and hydroponic culture experiments
The barley double haploid (DH) population used in the current study was derived by anther culture from a cross between Huaai 11 and Huadamai 6 in our research group. This population consisted of 122 DH lines, have been described in our previous studies [33, 34, 49,50,51,52]. In the present study, 122 lines and the two parental cultivars were evaluated at seedling stage using the hydroponic culture (Fig. 1a). The experiment was conducted in greenhouse during the year 2014–2015 at the College of Plant Science and Technology of Huazhong Agricultural University, Wuhan, China. The conditions of greenhouse were set for 25/15 °C during a 16/8 h light/dark cycle, a photon flux density of 300 μmol m−2 s−1 during the 16 h light period, and relative humidity of 65% for 24 h. Barley seeds of the two parents and 122 lines were surface sterilized for 15 min in a 5% solution of sodium hypochlorite, rinsed thoroughly with distilled water, soaked in deionised water for 2 h, and then germinated for 5 days on plastic board floating over the deionized water at room temperature. After 5 days, 6 uniform seedlings from each line were selected for transplanting into foam board (polystyrene) suspended in the plastic tanks (34 cm × 25 cm × 12 cm, length × width × height) with nutrient solution. Each seedling was anchored to holes 1.5 cm diameter in foam board with sterilized sponge. Each foam board had 5 rows of holes, and each row had 6 holes for a total of 30 seedlings per board. The distance between two seedlings was 3.6 cm (between rows) and 4.0 cm (within a row), respectively. Each tank contained eight liters of revised Hoagland’s nutrient solution , which was renewed every seventh day to prevent nutrient exhaustion. When renewed the nutrient solution, each tank moves randomly to ensure that same environmental conditions. The solution PH was adjusted to 6.5 using diluted NaOH and HCl before refreshing. The two parents and DH lines were grown in completely randomized design with six replicates. In each replication, one seedling was used for data collection.
Phenotyping of seedling characteristics
The phenotypes measured included leaf age (LAG), seedling height (SH, cm), maximum root length (MRL, cm), main roots number (MRN) and seedling fresh weight (SWF, g) at four growth stages (13th, 20th, 27th and 34th days after germination) on the same six seedlings for each line. Root fresh weight (RFW, g), stem-and-leaf fresh weight (SLFW, g), root dry weight (RDW, g) and stem-and-leaf dry weight (SLDW, g) were only measured at the 34th day after germination. For convenience, we used S1, S2, S3 and S4 to represent the 13th, 20th, 27th and 34th day after germination, respectively. In addition, leaf length (LL, cm), leaf width (LW, cm) and leaf area (LA, cm2) were measured in four leaves according to the leaf order of seedling. LL was measured for first leaf length (FILL, cm), second leaf length (SELL, cm), third leaf length (THLL, cm) and fourth leaf length (FOLL, cm). LW was measured for first leaf width (FILW, cm), second leaf width (SELW, cm), third leaf width (THLW, cm) and fourth leaf width (FOLW, cm). LA as measured for first leaf area (FILA, cm2), second leaf area (SELA, cm2), third leaf area (THLA, cm2) and fourth leaf area (FOLA, cm2). The methods of measurement are listed in Table 6.
The mean phenotypic values of the 21 seedling characteristics obtained from hydroponic culture were subjected to statistical analysis. Normality of distribution was tested using the method of Shapiro-Wilk. The descriptive statistics analyses were performed using SPSS programs (IBM SPSS Statistics, Chicago, IL, USA, http://www.ibm.com/analytics/us/en/technology/spss). The heritability was estimated . Frequency distribution and QTL-likelihood maps for the seedling characteristics were drawn using the Origin programs (OriginLab, Northampton, MA, USA, http://www.originlab.com). P-value higher than 0.05 was deemed as significance level.
The high-density genetic linkage map constructed by Ren et al.  for the ‘Huaai 11 × Huadamai 6’ population was used in the QTL analysis, which contains 1894 SNP markers and 68 SSR markers, covering all 7 chromosomes and spanning 1375.80 cM of the barley genome with an average marker distance of 0.7 cM. The estimates of QTL positions and effects were determined by inclusive composite interval mapping (ICIM)  using the software QTL IciMapping 4.1 . The mapping method was selected ICIM-ADD (additive effects) in ‘BIP (QTL mapping in biparental populations)’ function to analyze each trait from a single environment. Significant LOD (likelihood-of-odd) threshold for declaring a QTL of each trait was determined by 1000 permutations test with a Type 1 error of 0.05 . The scanning step size was set at 1.0 centimorgan (cM) and the probability in stepwise regression (PIN) was 0.001. The rule recommended by Liu et al.  was used to name the QTLs, which were mapped on linkage groups using the software MapChart 2.2 . The QTLs whose percentage of phenotypic variation explained (PVE) exceeded 20%, were considered as major QTL, and otherwise minor QTL. If the 1.5-LOD support interval of two QTLs overlapped, these two QTL would be defined as co-localized QTLs. However, we independently counted each QTL to clearly specify the number of detected QTLs. Graingenes website (http://wheat.pw.usda.gov/GG3/) was used to compare the maker information.
First leaf area
First leaf length
First leaf width
Fourth leaf area
Fourth leaf length
Fourth leaf width
Inclusive composite interval mapping
Likelihood of odd
Marker assisted selection
Quantitative trait locus
13th day after germination
20th day after germination
27th day after germination
34th day after germination
Second leaf area
Second leaf length
Second leaf width
Single nucleotide polymorphism
Simple sequence repeat
Third leaf area
Third leaf length
Third leaf width
Qin D, Dong J, Xu F, Guo G, Ge S, Xu Q, et al. Characterization and fine mapping of a novel barley stage green-revertible albino gene (HvSGRA) by bulked Segregant analysis based on SSR assay and specific length amplified fragment sequencing. BMC Genomics. 2015;16:838.
Heneen WK. Cytogenetics and molecular Cytogenetics of barley: a model cereal crop with a large genome. In: Ullrich SE, editor. Barley: production, improvement, and uses. UK: Wiley-Blackwell; 2011. p. 112–21.
Sbei H, Sato K, Shehzad T, Harrabi M, Okuno K. Detection of QTLs for salt tolerance in Asian barley (Hordeum vulgare L.) by association analysis with SNP markers. Breed Sci. 2014;64:378–88.
HTT V, Nguyen HTT, Tran KD, Khuat TH, Nakamura C. Genetic diversity of Vietnamese lowland rice germplasms as revealed by SSR markers in relation to seedling vigour under submergence. Biotechnol Biotechnol Equip. 2016;30:17–25.
Kan G, Ning L, Li Y, Hu Z, Zhang W, He X, et al. Identification of novel loci for salt stress at the seed germination stage in soybean. Breed Sci. 2016;66:530–41.
Nguyen TCT, Obermeier C, Friedt W, Abrams SR, Snowdon RJ. Disruption of germination and seedling development in Brassica napus by mutations causing severe seed hormonal imbalance. Front Plant Sci. 2016;7:322.
CG X, Li XQ, Xue Y, Huang YW, Gao J, Xing YZ. Comparison of quantitative trait loci controlling seedling characteristics at two seedling stages using rice recombinant inbred lines. Theor Appl Genet. 2004;109:640–7.
Bertholdsson NO. Early vigour and allelopathy–two useful traits for enhanced barley and wheat competitiveness against weeds. Weed Res. 2005;45:94–102.
Colmer TD, Munns R, Flowers TJ. Improving salt tolerance of wheat and barley: future prospects. Aust J Exp Agric. 2006;45:1425–43.
Abdulrahmani B, Ghassemi-Golezani K, Valizadeh M, Feizi-Asl V. Seed priming and seedling establishment of barley (Hordeum vulgare L.). J Food Agric Environ. 2007;5:179–84.
López-Castañeda C, Richards RA, Farquhar GD, Williamson RE. Seed and seedling characteristics contributing to variation in early vigor among temperate cereals. Crop Sci. 1996;36:1257–66.
Liao MT, Fillery IRP, Palta JA. Early vigorous growth is a major factor influencing nitrogen uptake in wheat. Funct Plant Biol. 2004;31:121–9.
Bertholdsson NO, Brantestam AK. A century of Nordic barley breeding–effects on early vigour root and shoot growth, straw length, harvest index and grain weight. Eur J Agron. 2009;30:266–74.
Ellis RH. Seed and seedling vigour in relation to crop growth and yield. Plant Growth Regul. 1992;11:249–55.
Song W, Wang B, Hauck AL, Dong X, Li J, Lai J. Genetic dissection of maize seedling root system architecture traits using an ultra-high density bin-map and a recombinant inbred line population. J Integr Plant Biol. 2016;58:266–79.
Atkinson JA, Wingen LU, Griffiths M, Pound MP, Gaju O, Foulkes MJ, et al. Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat. J Exp Bot. 2015;66:2283–92.
Tuberosa R, Sanguineti MC, Landi P, Giuliani MM, Salvi S, Conti S. Identification of QTLs for root characteristics in maize grown in hydroponics and analysis of their overlap with QTLs for grain yield in the field at two water regimes. Plant Mol Biol. 2002;48:697–712.
Kumar B, Abdel-Ghani AH, Pace J, Reyes-Matamoros J, Hochholdinger F, Lübberstedt T. Association analysis of single nucleotide polymorphisms in candidate genes with root traits in maize (Zea mays L.) seedlings. Plant Sci. 2014;224:9–19.
Wang H, Xu X, Zhan X, Zhai R, Wu W, Shen X, et al. Identification of qRL7, a major quantitative trait locus associated with rice root length in hydroponic conditions. Breed Sci. 2013;63:267–74.
Gamuyao R, Chin JH, Pariasca-Tanaka J, Pesaresi P, Catausan S, Dalid C, et al. The protein kinase Pstol1 from traditional rice confers tolerance of phosphorus deficiency. Nature. 2012;488:535–9.
Syed MA, Iftekharuddaula KM, Mian MAK, Rasul MG, Rahmam GKMM, Panaullah GM, et al. Main effect QTLs associated with arsenic phyto-toxicity tolerance at seedling stage in rice (Oryza sativa L.). Euphytica. 2016;209:805–14.
Burton AL, Johnson JM, Foerster JM, Hirsch CN, Buell CR, Hanlon MT, et al. QTL mapping and phenotypic variation for root architectural traits in maize (Zea mays L.). Theor Appl Genet. 2014;127:2293–311.
Abdel-Ghani AH, Kumar B, Reyes-Matamoros J, Gonzalez-Portilla PJ, Jansen C, Martin JPS, et al. Genotypic variation and relationships between seedling and adult plant traits in maize (Zea mays L.) inbred lines grown under contrasting nitrogen levels. Euphytica. 2013;189:123–33.
Liu X, Li R, Chang X, Jing R, Mapping QTL. For seedling root traits in a doubled haploid wheat population under different water regimes. Euphytica. 2013;189:51–66.
Kabir MR, Liu G, Guan P, Wang F, Khan AA, Ni Z, et al. Mapping QTLs associated with root traits using two different populations in wheat (Triticum aestivum L.). Euphytica. 2015;206:175–90.
Ellis RP, Forster BP, Gordon DC, Handley LL, Keith RP, Lawrence P, et al. Phenotype/genotype associations for yield and salt tolerance in a barley mapping population segregating for two dwarfing genes. J Exp Bot. 2002;53:1163–76.
Bertholdsson NO, Holefors A, Macaulay M, Crespo-Herrera LAQTL. For chlorophyll fluorescence of barley plants grown at low oxygen concentration in hydroponics to simulate waterlogging. Euphytica. 2015;201:357–65.
Bertholdsson NO. Screening for barley waterlogging tolerance in nordic barley cultivars (Hordeum vulgare L.) using chlorophyll fluorescence on hydroponically-grown plants. Agronomy. 2013;3:376–90.
Xu R, Wang J, Li C, Johnson P, Lu C, Zhou MA. Single locus is responsible for salinity tolerance in a chinese landrace barley (Hordeum vulgare L.). PLoS One. 2012;7:e43079.
Hoffmann A, Maurer A, Pillen K. Detection of nitrogen deficiency QTL in juvenile wild barley introgression lines growing in a hydroponic system. BMC Genet. 2012;13:88.
González AM, Yuste-Lisbona FJ, Rodiño AP, Ron AMD, Capel C, García-Alcázar M, et al. Uncovering the genetic architecture of Colletotrichum lindemuthianum resistance through QTL mapping and epistatic interaction analysis in common bean. Front Plant Sci. 2015;6:141.
Paterson AH, Lander ES, Hewitt JD, Peterson S, Lincoln SE, Tanksley SD. Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms. Nature. 1988;335:721–6.
Liu L, Sun G, Ren X, Li C, Sun D. Identification of QTL underlying physiological and morphological traits of flag leaf in barley. BMC Genet. 2015;16:29.
Wang J, Sun G, Ren X, Li C, Liu L, Wang Q, et al. QTL underlying some agronomic traits in barley detected by SNP markers. BMC Genet. 2016;17:103.
Zhou G, Panozzo J, Zhang XQ, Cakir M, Harasymow S, Li C. QTL mapping reveals genetic architectures of malting quality between Australian and Canadian malting barley (Hordeum vulgare L.). Mol Breed. 2016;36:70.
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.
Xue DW, Zhou MX, Zhang XQ, Chen S, Wei K, Zeng FR, et al. Identification of QTLs for yield and yield components of barley under different growth conditions. J Zheijang Univ Sci B. 2010;11:169–76.
Szűcs P, Blake VC, Bhat PR, Chao S, Close TJ, Cuesta-Marcos A, et al. An integrated resource for barley linkage map and malting quality QTL alignment. The Plant Genome. 2009;2:134–40.
Panozzo JF, Eckermann PJ, Mather DE, Moody DB, Black CK, Collins HM, et al. QTL analysis of malting quality traits in two barley populations. Aust J Agric Res. 2007;58:858–66.
Ahmadi-Ochtapeh H, Soltanloo H, Ramezanpour SS, Naghavi MR, Nikkhah HR, Rad SY. QTL mapping for salt tolerance in barley at seedling growth stage. Biol Plantarum. 2015;59:283–90.
Siahsar BA, Narouei M. Mapping QTLs of physiological traits associated with salt tolerance in ‘Steptoe’×‘Morex’doubled haploid lines of barley at seedling stage. J Food Agr Environ. 2010;88:751–9.
Broughton S, Zhou G, Teakle NL, Matsuda R, Zhou M, O’Leary RA, et al. Waterlogging tolerance is associated with root porosity in barley (Hordeum vulgare L.). Mol Breed. 2015;35:27.
Chen G, Krugman T, Fahima T, Chen K, Hu Y, Röder M, et al. Chromosomal regions controlling seedling drought resistance in Israeli wild barley, Hordeum spontaneum C. Koch. Genet Resour Crop Evol. 2010;57:85–99.
Asseng S, Turner NC, Botwright T, Condon AG. Evaluating the impact of a trait for increased specific leaf area on wheat yields using a crop simulation model. Agron J. 2003;95:10–9.
Nicola S. Understanding root systems to improve seedling quality. HortTechnology. 1998;8:544–9.
Arifuzzaman M, Sayed MA, Muzammil S, Pillen K, Schumann H, Naz AA, et al. Detection and validation of novel QTL for shoot and root traits in barley (Hordeum vulgare L.). Mol Breed. 2014;34:1373–87.
Songsri P, Jogloy S, Holbrook CC, Kesmala T, Vorasoot N, Akkasaeng C, et al. Association of root, specific leaf area and SPAD chlorophyll meter reading to water use efficiency of peanut under different available soil water. Agric Water Manag. 2009;96:790–8.
Sabiel SA, Huang S, Hu X, Ren X, Fu C, Peng J, et al. SNP-based association analysis for seedling traits in durum wheat (Triticum turgidum L. durum (Desf.)). Breed Sci. 2017;67:83–94.
Ren X, Sun D, Sun G, Li C, Dong W. Molecular detection of QTL for agronomic and quality traits in a doubled haploid barley population. Aust J Crop Sci. 2013;7:878–86.
Ren X, Wang J, Liu L, Sun G, Li C, Luo H, et al. SNP-based high density genetic map and mapping of btwd1 dwarfing gene in barley. Sci Rep. 2016;6:31741.
Ren X, Sun D, Dong W, Sun G, Li C. Molecular detection of QTL controlling plant height components in a doubled haploid barley population. Genet Mol Res. 2014;13:3089–99.
Ren X, Sun D, Guan W, Sun G, Li C. Inheritance and identification of molecular markers associated with a novel dwarfing gene in barley. BMC Genet. 2010;11:89.
Arifuzzaman M, Günal S, Bungartz A, Muzammil S, Afsharyan NP, Léon J, et al. Genetic mapping reveals broader role of Vrn-H3 gene in root and shoot development beyond heading in barley. PLoS One. 2016;11:e0158718.
Naz AA, Arifuzzaman M, Muzammil S, Pillen K, Léon J. Wild barley introgression lines revealed novel QTL alleles for root and related shoot traits in the cultivated barley (Hordeum vulgare L.). BMC Genet. 2014;15:107.
Reinert S, Kortz A, Léon J, Naz AA. Genome-wide association mapping in the global diversity set reveals new QTL controlling root system and related shoot variation in barley. Front Plant Sci. 2016;7:1061.
Schmalenbach I, Pillen K. Detection and verification of malting quality QTLs using wild barley introgression lines. Theor Appl Genet. 2009;118:1411–27.
WR W, Li WM, Tang DZ, HR L, Worland AJ. Time-related mapping of quantitative trait loci underlying tiller number in rice. Genetics. 1999;151:297–303.
Mansour E, Casas AM, Gracia MP, Molina-Cano JL, Moralejo M, Cattivelli L, et al. Quantitative trait loci for agronomic traits in an elite barley population for Mediterranean conditions. Mol Breed. 2014;33:249–65.
Qu Y, Mu P, Zhang H, Chen CY, Gao Y, Tian Y, et al. Mapping QTLs of root morphological traits at different growth stages in rice. Genetica. 2008;133:187–200.
Price AH, Steele KA, Moore BJ, Jones RGW. Upland rice grown in soil-filled chambers and exposed to contrasting water-deficit regimes: II. Mapping quantitative trait loci for root morphology and distribution. Field Crops Res. 2002;76:25–43.
Marquez-Cedillo LA, Hayes PM, Kleinhofs A, Legge WG, Rossnagel BG, Sato K, et al. 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:625–37.
Marquez-Cedillo LA, Hayes PM, Jones BL, Kleinhofs A, Legge WG, Rossnagel BG, 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.
Sato K, Hori K, Takeda K. Detection of Fusarium head blight resistance QTLs using five populations of top-cross progeny derived from two-row × two-row crosses in barley. Mol Breed. 2008;22:517–26.
Xue D, Huang Y, Zhang X, Kang W, Westcott S, Li C, et al. Identification of QTLs associated with salinity tolerance at late growth stage in barley. Euphytica. 2009;169:187–96.
Ehrich TH, Vaughn TT, Koreishi SF, Linsey RB, Pletscher LS, Cheverud JM. Pleiotropic effects on mandibular morphology I. Developmental morphological integration and differential dominance. J Exp Zool B Mol Dev Evol. 2003;296:58–79.
Hori K, Kobayashi T, Shimizu A, Sato K, Takeda K, Kawasaki S. Efficient construction of high-density linkage map and its application to QTL analysis in barley. Theor Appl Genet. 2003;107:806–13.
Chutimanitsakun Y, Nipper RW, Cuesta-Marcos A, Cistué L, Corey A, Filichkina T, et al. Construction and application for QTL analysis of a restriction site associated DNA (RAD) linkage map in barley. BMC Genomics. 2011;12:4.
Castro AJ, Hayes P, Viega L, Vales I. Transgressive segregation for phenological traits in barley explained by two major QTL alleles with additivity. Plant Breed. 2008;127:561–8.
Hori K, Sato K, Takeda K. Detection of seed dormancy QTL in multiple mapping populations derived from crosses involving novel barley germplasm. Theor Appl Genet. 2007;115:869–76.
Peng J, Ronin Y, Fahima T, Röder MS, Li Y, Nevo E, et al. Domestication quantitative trait loci in Triticum dicoccoides, the progenitor of wheat. Proc Natl Acad Sci U S A. 2003;100:2489–94.
Wang H, Huang C, Guo H, Li X, Zhao W, Dai B, et al. QTL mapping for fiber and yield traits in upland cotton under multiple environments. PLoS One. 2015;10:e0130742.
Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. Calif Agric Exp Stat Circ. 1950;347:1–32.
Hanson C, Robinson H, Comstock R. Biometrical studies of yield in segregating populations of Korean lespedeza. Agron J. 1956;48:268–72.
Li H, Ye G, Wang JA. Modified algorithm for the improvement of composite interval mapping. Genetics. 2007;175:361–74.
Meng L, Li H, Zhang L, Wang J. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. The Crop J. 2015;3:269–83.
Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138:963–71.
Voorrips RE. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93:77–8.
This project was supported by the earmarked fund for China Agriculture Research System (CARS-5).
Availability of data and materials
All relevant datasets supporting the conclusions of this article are available within the article and its additional files.
Ethics approval and consent to participate
Our research did not involve any human subjects, human material, or human data. All our studies have been performed in accordance with ethics committee polices. The hydroponic culture experiment in greenhouse at the College of Plant Science and Technology of Huazhong Agricultural had been allowed and not involved private interests.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Frequency distribution of 21 seedling characteristics in DH population. P value of Shapiro-Wilk test for each stage was shown, the hypothesis of normal distribution was accepted when P > 0.05 (significant at P = 0.05), and the trend lines of the accepted normal distribution were shown. Trait abbreviations refer to the Table 6. S1 to S4 represented the 13th, 20th, 27th and 34th day after germination, respectively. (JPEG 5396 kb)
QTL likelihood map for 17 seedling characteristics in the DH population grown under hydroponic culture conditions using inclusive composite interval mapping. Genetic maps (all chromosomes together) of barley linkage groups were shown in the abscissa and LOD scores of each trait in the ordinate. The significant LOD threshold was determined to be 3.0 by 1000 permutations test. Trait abbreviations refer to the Table 6. S1 to S4 represented the 13th, 20th, 27th and 34th day after germination, respectively. (JPEG 4454 kb)
QTL locations for previous studies and current studies detected in the Huaai 11 × Huadamai 6 DH population. QTLs were mapped to chromosomes using peak position and nearest marker. (JPEG 4137 kb)
About this article
Cite this article
Wang, Q., Sun, G., Ren, X. et al. Detection of QTLs for seedling characteristics in barley (Hordeum vulgare L.) grown under hydroponic culture condition. BMC Genet 18, 94 (2017) doi:10.1186/s12863-017-0562-y
- Barley (Hordeum vulgare L.)
- Seedling characteristics
- Hydroponic culture
- Quantitative trait loci (QTL)
- Dynamic QTL