- Research article
- Open Access
Identification of QTL underlying physiological and morphological traits of flag leaf in barley
BMC Genetics volume 16, Article number: 29 (2015)
Physiological and morphological traits of flag leaf play important roles in determining crop grain yield and biomass. In order to understand genetic basis controlling physiological and morphological traits of flag leaf, a double haploid (DH) population derived from the cross of Huaai 11 × Huadamai 6 was used to detect quantitative trait locus (QTL) underlying 7 physiological and 3 morphological traits at the pre-filling stage in year 2012 and 2013.
Total of 38 QTLs distributed on chromosome 1H, 2H, 3H, 4H, 6H and 7H were detected, and explained 6.53% - 31.29% phenotypic variation. The QTLs flanked by marker Bmag829 and GBM1218 on chromosome 2H were associated with net photosynthetic rate (Pn), stomatal conductance (Gs), flag leaf area (LA), flag leaf length (FLL), flag leaf width (FLW), relative chlorophyll content (SPD) and leaf nitrogen concentration (LNC).
Two QTL cluster regions associated with physiological and morphological traits, one each on the chromosome 2H and 7H, were observed. The two markers (Bmag829 and GBM1218) may be useful for marker assisted selection (MAS) in barley breeding.
Barley (Hordeum vulgare L.) is the fourth cereal crop in world production . High yield is always one of the important barley breeding aims . However, grain yield was controlled by complex biochemical and physiological processes, and closely related to physiological and morphological traits [3-7]. The top three leaves on a stem, especially the flag leaf, absorb most irradiation light, and were the primary source of carbohydrate production . In barley, importance of flag leaf on increasing grain yield has widely been studied [6,7,9]. However, previous studies have mainly focused on either morphological traits [10-12] or physiological traits of flag leaf [13-18] determining grain yield. Few QTLs associated with these traits have been applied to barley breeding due to complicated measurement procedure, inconsistency and dynamic process of physiological and morphological traits in barley developmental stage. Thus, comprehensive understanding the role of physiological and morphological traits of flag leaf on yield will provide a new insight in crop growth and development. Meanwhile, application of molecular marker and genetic map made it possible to map the region controlling quantitative traits [11,19,20].
Increasing photosynthetic capacity of leaf is one of the most important approaches to increase crop biomass . It was estimated that leaf photosynthesis contributing 30% biomass . Photosynthesis is an essential process to maintain crop growth and development. Photosynthetic capacity during reproductive stage is positively correlated with crop yield . Four main physiological parameters: net photosynthetic rate, stomatal conductance, intercellular CO2 concentration and transpiration rate, have been used to evaluate photosynthetic capacity. Teng et al.  reported that net photosynthetic rate in rice was controlled by multiple genes. In barley, QTL underlying net photosynthetic rate has been analyzed in two DH populations . According to Jiang et al. , stomatal conductance significantly affected net photosynthetic rate, and is a key parameter to assess limitation of photosynthesis in barley. Rybiński et al.  found significant linear relationship between transpiration rate and net photosynthetic rate in different irradiated times under laser light. However, the QTLs underlying stomatal conductance, intercellular CO2 concentration and transpiration rate have not been reported in barley.
Chlorophyll absorbs light energy and converts it into chemical energy. Maintaining higher level of chlorophyll content in leaf is one of the strategies for increasing photosynthesis and crop production . The structure and function of chloroplasts determine photosynthetic activity . Von Kroff et al.  reported a positive correlation between relative chlorophyll fluorescence in leaf and grain yield. The chlorophyll content was suggested as a reliable indicator for evaluating metabolic balance between photosynthesis and yield performance . Recently, chlorophyll content in barley leaf has widely been studied [11,14,26,28].
Nitrogen uptake and metabolism of flag leaf at the pre-filling stage provide main energy source to grain yield . The photosynthetically active leaf cells of chloroplasts contain most nitrogen . The most of assimilated nitrogen mainly come from photosynthesis. Leaf CO2 assimilation rate and nitrogen content per unit area was highly correlated . Depending on physiological status, nitrogen can be stored and assimilated in both leaves and roots . In fully developed leaves, about 75% nitrogen is allocated to chloroplasts, and mostly used for synthesizing components of photosynthetic apparatus . A positive correlation was found between photosynthetic capacity of leaves and their nitrogen content . In past few years, some studies have reported that nitrogen content in leaves was quantitative trait and controlled by multiple genes in barley. Stable QTLs were detected, but phenotypic contribution from each QTL was small [12,15,29].
Plant water status plays an important role in plant growth, development, and keeping yield stability . The physiological and morphological traits such as photosynthesis, transpiration of flag leaves and grain yield are closely correlated with plant water status [35,36]. In water deficit environment, crop must increases water use efficiency to resist drought, and sustains normal growth . Relative water content (RWC) was widely used to measure water status in barley . RWC is an important determinant of leaf metabolic activity, and reflects water balance in tissues . Maintenance of certain level of RWC can increase yield and its stability in cereals . As RWC is related to plant water-status, it can be used to evaluate water level in plant at a specific growth stage. It has been reported that RWC has a positive relationship with yield in cereals . QTLs associated with RWC were detected on chromosome 6H in different water conditions and developmental stages [16,40,41].
In present study, a DH population derived from the cross of Huaai 11 × Huadamai 6 was used to identify QTLs underlying physiological and morphological traits of flag leaf at the pre-filling stage. The identified QTLs can be used for molecular assisted selection (MAS) in barley breeding.
Phenotype analysis of the double population and parents
The statistics of 7 physiological and 3 morphological traits of flag leaf at the pre-filling stage were shown in Table 1. The values of Pn, Gs, Ci, Tr, RWC, SPD and LNC in Huaai 11 were higher than those in Huadamai 6. The values of LA, FLL and FLW were higher in Huadamai 6 than those in Huaai 11. The t-test showed that two parents were significant difference on all traits (p < 0.05). All traits displayed a normal distribution with the skewness and kurtosis among −1 and 1 (Table 1). Analysis of variance showed that genotype effects were significant (P < 0.01) for all traits studied. Effects between years were not significant (P > 0.05) except Pn, Gs and Tr traits. Genotype × year interactions were significant (P < 0.05) for all traits except LA, FLL and FLW (Table 2). In addition, all 7 physiological and 3 morphological traits at the pre-filling stage showed highly phenotypic variation in the DH population. The variable coefficients ranged from 5.22% to 30.91% in 2012, and 11.50% to 28.50% in 2013. Transgressive segregation in both directions was observed for all traits (Table 1). Heritability (Table 1) ranged from 44.13% to 80.67% and 52.66% to 85.57% in 2012 and 2013, respectively.
Correlations among Pn, Gs, Ci, and Tr were significant positive (P < 0.01, Table 3). Three morphological traits, LA, FLL and FLW, were also significantly positive correlated with each other (P < 0.01, Table 3). Significant positive correlation between Pn and SPD was detected with correlation coefficient of 0.335 in 2012 and 0.265 in 2013 (P < 0.01). LNC was significantly correlated with SPD (r = 0.283 in 2012 and 0.381 in 2013, P < 0.01). A negative correlation between Pn and LA was observed with r = −0.515 (year 2012) and −0.225 (year 2013) (P < 0.05). RWC was not significantly (P > 0.05) correlated with other traits except LA in 2013.
A total of 38 QTLs for 7 physiological and 3 morphological traits were detected and mapped (Figure 1; Table 4). 18 and 15 QTLs were detected in 2012 and 2013, respectively. Five QTLs based on mean value of each trait were detected for LA, FLL and FLW. The detected QTLs accounted for 7.14% - 24.58% and 6.53% - 25.36% phenotypic variation in 2012 and 2013, respectively. The QTLs based on mean values of LA, FLL and FLW explained 14.23% - 31.29% phenotypic variation.
Net photosynthetic rate
Three QTL underlying Pn trait were detected. Two QTLs, qPn2-10 and qPn4-17, were detected on chromosome 2H and 4H in 2012. They accounted for 8.66% and 12.63% total phenotypic variation, respectively. The QTL, qPn7-8 on chromosome 7H was detected in 2013, and accounted for 13.56% total phenotypic variation. Both qPn2-10 and qPn7-8 QTLs have alleles from Huaai 11 to increase net photosynthetic rate, the QTL qPn4-17 has allele from Huadamai 6 to increase net photosynthetic rate (Figure 1; Table 4).
Four QTLs associated with Gs trait were detected. Of them, three QTLs, qGs2-10, qGs3-13 and qGs7-6, were detected in 2012 and mapped on chromosome 2H, 3H and 7H, and accounted for 7.78%, 12.58% and 13.92% total phenotypic variation, respectively. In 2013, one QTL qGs2-13 was detected on chromosome 2H, and accounted for 7.47% total phenotypic variation. All these QTLs have alleles from Huaai 11 to increase stomatal conductance, their values ranged from 0.04 to 0.07 (Figure 1; Table 4).
Intercellular CO2 concentration
Three QTLs for Ci trait were detected. Of them, two QTLs, qCi2-16 and qCi7-3, were mapped on chromosome 2H and 7H in 2012, and accounted for 13.75% and 13.98% total phenotypic variation, respectively. One QTL qCi2-14 was identified in 2013, and accounted for 10.69% total phenotypic variation. These QTLs have alleles from Huaai 11 to increase intercellular CO2 concentration (Figure 1; Table 4).
Two QTLs underlying Tr trait were identified in 2012. The QTL qTr3-13 and qTr7-6 accounted for 14.00% and 14.02% total phenotypic variation, respectively. The additive effects of the two QTLs were 0.69 and 0.71, respectively, indicating that the alleles from Huaai 11 increased transpiration rate (Figure 1; Table 4).
Flag leaf area
Four QTLs underlying LA trait were detected on chromosome 2H and 3H. The QTL, qLA2-12 close to the marker GBM1218, was detected in both years and mean value, and accounted for 18.80% (year 2012), 12.48% (year 2013) and 29.83% (mean value from two years) phenotypic variation. The alleles from Huadamai 6 increased flag leaf area. Another QTL qLA3-9 detected in 2013 accounted for 8.72% phenotypic variation. The allele of QTL qLA3-9 from Huaai 11 increased flag leaf area (Figure 1; Table 4).
Flag leaf length
Seven QTLs associated with FLL trait were detected. The QTL, qFLL2-12 close to the marker GBM1218 on chromosome 2H, was detected in both years and mean value, and accounted for 24.58% (year 2012), 25.36% (year 2013) and 31.29% (mean value from two years) phenotypic variation. The alleles of the QTL, which increased flag leaf length, came from Huadamai 6. Other four QTLs, qFLL7-10, qFLL3-11, qFLL7-6 and qFLL7-8, accounted for 13.04%, 9.76%, 7.07% and 16.66% total phenotypic variation, respectively. The positive alleles of QTL qFLL7-10, qFLL3-11, qFLL7-6 and qFLL7-8 from Huadamai 6 contributed to the increase in flag leaf length by 1.06, 0.98, 0.79 and 1.14, respectively (Figure 1; Table 4).
Flag leaf width
For FLW trait, five putative QTLs were identified. The QTL, qFLW2-12 close to the marker GBM1218 on chromosome 2H, was detected in both years and mean value, and accounted for 13.63% (year 2012), 20.93% (year 2013) and 14.23% (mean value from two years) total phenotypic variation. The positive alleles of QTL qFLW2-12 from Huadamai 6 increased flag leaf width. Another QTL qFLW4-18 detected in 2013 and mean value was located on chromosome 4H, and accounted for 7.11% and 22.06% total phenotypic variation, respectively. The alleles of qFLW4-18 from Huaai 11 contributed to the increase in flag leaf width (Figure 1; Table 4).
Relative water content
Three QTLs underlying RWC were found. The QTL qRWC6-6 nearby the marker GMS6 on chromosome 6H was detected in both years, and accounted for 21.43% (year 2012) and 11.76% (year 2013) phenotypic variation. Their alleles from Huadamai 6 increased relative water content. Another QTL, qRWC7-9 was detected in year 2012 and mapped on chromosome 7H, which accounted for 15.31% phenotypic variation. The allele from Huaai 11 increased relative water content (Figure 1; Table 4).
Relative chlorophyll content
Four QTLs underlying SPD trait were found. The QTL qSPD2-10 was detected in both years and close to the marker Bmag829 on chromosome 2H, and accounted for 17.28% (year 2012) and 15.44% (year 2013) total phenotypic variation. Two QTLs, qSPD7-7 and qSPD7-9, were mapped on chromosome 7H and close to the marker Bmac167 (year 2012) and Bmag746 (year 2013). They accounted for 10.78% and 10.64% total phenotypic variation in year 2012 and 2013, respectively. All these QTLs have alleles from Huaai 11 contributed to the increase in relative chlorophyll content (Figure 1; Table 4).
Total nitrogen content
Three QTLs associated with LNC trait were detected. Of them, one QTL, qLNC1-10 on chromosome 1H, was detected in 2012 and accounted for 7.14% phenotypic variation. Two QTLs qLNC1-8 and qLNC2-10 were mapped on chromosome 1H and 2H in 2013, and accounted for 8.46% and 6.53% phenotypic variation, respectively. All these QTLs have alleles from Huaai 11 contributed to the increase in total nitrogen content (Figure 1; Table 4).
QTL analysis is a useful approach to discover and identify favorable alleles in barley . Ren et al.  have studied the correlation and QTL of agronomic and quality traits associated with grain yield in a barley DH population. However, QTL associated with physiological and morphological traits of flag leaf at the pre-filling stage have not been systematically analyzed.
Leaf net photosynthetic rate was easily affected by environment factors. It was reported the net photosynthetic rate was different in different environments including illumination intensity, temperature, content of CO2 and moisture in the air . In our experiment, we selected 9:00–11:00 am and 2:00–4:00 pm to measure photosynthesis based on the daily change rule of photosynthesis and our operational experience that photosynthesis was stable at these two time periods. In plant developmental stage, the four traits Pn, Gs, Ci and Tr index reflect plant photosynthetic capacity. The all four traits were closely related to grain yield. QTLs underlying Pn, Gs and Tr have been analyzed in rice . Wójcik-Jagła et al.  analyzed QTL underlying net photosynthetic rate in barley, and found one QTL nearby the marker bPb-8013 on chromosome 4H in the Suweren × MOB12055 population, one QTL on chromosome 5H in the STH754 × STH836 population. In our study, we detected one QTL nearby the marker EBmac788 on chromosome 4H. The consensus map of Wenzl et al.  showed that the marker bPb-8013 is far from EBmac788, indicating that the qPn4-17 was a new QTL identified here. In rice, QTL analysis of several physiological traits related to photosynthesis had been performed . In our study, 9 QTLs controlling Gs, Ci and Tr traits in barley flag leaf were detected. The identified QTLs may be useful for MAS in barley breeding.
To sustain crop growth and development, crop must produce abundant nutrition. The amount of nutrition produced mainly depends on flag leaf associated with Pn, SPD, LNC and LA, which were closely related to grain yield and biomass [3,7,9]. Four QTLs associated with relative chlorophyll content were detected. QTL qSPD2-10 was detected at 75.9 cM in 2012 and 2013, indicating this QTL was stable and less affected by environments. In barley, This et al.  detected 12 QTLs underlying chlorophyll content on chromosome 2H, 4H, 5H, 6H and 7H. Xue et al.  detected two QTLs underlying chlorophyll content on chromosome 2H. One QTL related to SPD trait has mapped on chromosome 2H . The high density consensus map  indicated the qSPD2-10 was close to the QTL (qFC2.2) , between marker Bmag0518 and Bmac0093. The QTL qSPD7-7 and qSPD7-9 were close to the centromere of chromosome 7H, and different from the QTL on chromosome 7H reported previously [17,28]. Five QTLs controlling nitrogen content of flag leaf were detected on chromosome 2H, 3H, 5H and 7H . Mickelson et al.  detected 19 QTLs on chromosome 3H, 4H, 5H, 6H and 7H associated with nitrogen concentration in flag leaf. Three QTLs underlying LNC trait were detected on chromosome 1H and 2H in our study, indicating that the two QTLs on chromosome 1H may be new QTL underlying nitrogen concentration in flag leaf. The QTL qLNC2-10 on centromere region of chromosome 2H is different from the QTL on chromosome 2H reported previously . Four QTLs associated with flag leaf area were identified. The QTL qLA2-12 on chromosome 2H located at 77.2 cM was detected in both years and mean value. Previous studies reported QTL underlying leaf area on chromosome 1H, 2H, 3H, 4H, 5H and 7H [12,45]. The qLA2-12 on 2HL is different from the QTL reported on 2HS . In our study, one region on chromosome 2H flanked by Bmag829 and GBM1218 contained the qPn2-10, qLA2-12, qSPD2-10 and qLNC2-10 (Figure 1), suggesting that there might be QTL cluster for controlling grain yield on chromosome 2H, and these molecular makers can be used for MAS to improve breeding efficiency.
Since year effects and genotype × year interactions were not significant (p > 0.05) for three morphological traits (LA, FLL, FLW), QTL analysis was performed for data from each year and mean value of two years. In our study, 16 QTLs associated with the 3 morphological traits (LA, FLL and FLW) were identified in two years and mean values, which located on chromosome 2H, 3H, 4H and 7H, respectively. Elberse et al.  detected 6 QTLs underlying leaf length on chromosome 1H, 2H, 4H and 5H, 3 QTLs controlling leaf width on chromosome 2H, 4H and 6H. Li et al.  reported a chromosome region on 3HS underlying leaf length and leaf area. Gyenis et al.  reported 3 QTLs controlling flag leaf length on chromosome 3H, 5H and 7H, and 3 QTLs underlying flag leaf width on 2H, 4H and 5H. Xue et al.  detected 2 QTLs controlling flag leaf length on chromosome 5H and 7H, and 2 QTLs controlling flag leaf width on chromosome 5H. The QTL qFLL2-12 located on chromosome 2HL, and is different from the QTL reported on 2HS . The QTL, qFLW2-12 located on chromosome 2HL, and is different from those QTLs reported on 2HS [10,46]. The 3 morphological traits were significantly correlated with each other (Table 3), a common QTL close to the marker GBM1218 on chromosome 2H controlled these traits (Figure 1; Table 4). Phenotypic correlations among traits and identification of QTL were generally in good agreement. QTLs controlling LA, FLL and FLW were detected on the same region of chromosome 2H in both years and mean values. This region was close to the marker GBM1218, and contained the qLA2-12, qFLL2-12 and qFLW2-12 (Figure 1), indicating that this region is important for controlling morphological trait in barley. Moreover, all QTL positive alleles except qLA3-9 and qFLW4-18 were contributed by Huadamai 6.
Photosynthesis process assimilates H2O and CO2 to produce carbohydrates, and can be influenced by plant water status. Relative water content of flag leaf is one important assessment criterion about plant water status . In our study, one common QTL on the chromosome 6H is close to marker GMS6. Teulat et al.  detected one QTL on the chromosome 6H under two different water treatments. Another study also detected two QTLs on the long arm of chromosome 6H . Previous studies on QTL underlying RWC trait of barley flag leaf found 2 genome regions on the chromosome 6H associated with RWC, which were close to BCD348B and BCD1, respectively [13,16,40,41]. These suggested that there might be a QTL cluster in this region. Chromosome 7H have 3 genome regions associated with RWC, which are nearby RZ123, Acl3 and Bass1B, respectively [13,16,40,41]. The QTL qRWC6-6 detected in present study was close to the marker BCD348B, and the QTL qRWC7-9 was close to the marker RZ123.
In our study, two QTL cluster regions associated with physiological and morphological traits, one each on the chromosome 2H and 7H, were observed (Figure 1). The head type trait was mapped on chromosome 2H between marker GBM1218 and Bmac93, which is close to the QTL cluster region (Figure 1). The heading date trait was also mapped on chromosome 2H close to marker GBM1218 in the QTL cluster region . The dwarfing gene was mapped on chromosome 7H in the QTL cluster region . The head type, heading date and plant height traits were considered to be significantly associated with grain yield [43,49,50]. The vrs1 locus controlling head type was mapped on chromosome 2H [51,52]. From http://wheat.pw.usda.gov/GG2/index.shtml, we found that the marker GBM1218 was close to vrs1 locus. Considering all information here, we suggested that the head type, heading date and plant height traits might be highly associated with these physiological and morphological traits, and could be considered as important factors to control grain yield. Pleiotropy and linkage were present in some important traits associated with yield parameters . In present study, there exist widely co-localized QTL between physiological and morphological traits, such as Pn, Gs, SPD, LNC traits on chromosome 2H nearby the marker Bmag829, and LA, FLL, FLW traits on chromosome 2H nearby the marker GBM1218, where the vrs1 locus was mapped to. There is always a concentration of QTL effects in the vrs1 locus. The co-localization of these QTL is most likely due to pleiotropic effect or gene linkage. Distinguishing linkage from pleiotropy is important for breeding purposes, especially if both desirable and undesirable traits are associated with the same locus or QTL region . Thus, in order to distinguish linkage and pleiotropy, further study is needed.
In this study, physiological and morphological traits showed significant difference in two parents Huaai 11 and Huadamai 6. We found that chromosome 2H and 7H each contained a QTL cluster region controlling grain yield. The molecular makers (Bmag829 and GBM1218) identified here can be used for marker assisted selection to improve breeding efficiency.
Plant materials and field experiments
A barley DH population consisting of 122 DH lines was derived from a cross between dwarfing barley cultivar Huaai 11 (six-rowed and dwarfing) and common feed barley cultivar Huadamai 6 (two-rowed and tall plant) using anther culture. The two parents Huaai 11 and Huadamai 6 are significant difference in plant height , physiological and morphological traits of flag leaf. Experiment was conducted in a rain shelter of the Huazhong Agricultural University, Wuhan, China. Side window of the rain shelter was open to make inside temperature and radiation similar to outside condition. The experiments were performed in year 2012 and 2013. The DH lines and parents were grown in a plot of 1.5 m long with interval of 0.6 m and 3 replications using a randomized complete block design. Twenty seeds from each DH line and parent were sown in two rows per plot. Prior to seeding, compound fertilizer (60 g/m2) was applied, and 20 g/m2 of urea were applied at the elongation stage. At the pre-filling stage, fully expanded flag leaves from main spike were sampled and used to measure 7 physiological and 3 morphological traits.
Quantification of physiological traits of flag leaf at the pre-filling stage
Four physiological traits, net photosynthetic rate (Pn, umol CO2 m−2 s−1), stomatal conductance (Gs, mol H2O m−2 s−1), intercellular CO2 concentration (Ci, μmol CO2 mol−1) and transpiration rate (Tr, mmol H2O m−2 s−1), were measured using LI6400 XT Portable photosynthesis system according to the methods described in . Measuring time was selected during 9:00–11:00 am and 2:00–4:00 pm. Three fully expanded and sun-exposed topmost flag leaves on main stem from each replication were measured. The parameters were set as follow: LeafFan at Fast, Flow at 500 μmols−1, Mixer at 400 ppm, Temp at off and Lamp according to the light intensity. The data was recorded after these parameters reading became relatively stable (usually about 1 min).
Weighing method was applied to measure relative water content (RWC) in flag leaves . A flag leaf was sampled from each replication and measured 3 times. After fresh leaves weighted (fw), leaves were immersed in a sealed bag containing distilled water, and kept for 24 hours to achieve completely rehydration. Then the turgid leaves were weighted (tw), and dried to constant weight (dw). RWC was calculated as: RWC = (fw-dw)/(tw-dw) × 100%.
SPAD-502 chlorophyll photometer was used to measure relative chlorophyll content (SPD) of flag leaves at the pre-filling stage. Four flag leaves from each replication were measured. SPD values in the top, medium and bottom part of flag leaf were averaged from three replications.
Leaf nitrogen concentration (LNC) was measured using the Kjeldahl Nitrogen determination method. Ten flag leaves from each replication were collected at the pre-filling stage, immediately dried at 105°C in an oven for at least 4 h and then ground into powder using Whirlwind grinding JFS-13A, and stored at 80°C until use. Hanon SH220 was used to digest 0.2 g flag leaf powder. The digestive juice was put in distillation Hanon K9840 Kjeldahl Auto Analyzer to measure consumed volume of standard HCL. Total nitrogen in flag leaf (%) was calculated using the formula:
Where: C is concentration of standard HCL in the titration (mol/L); V is consumed volume of standard HCL in the titration sample (ml); V 0 is consumed volume of standard HCL in the titration blank group (ml); 14 is the atomic mass of nitrogen (g); 100 is total volume of digestive juice (ml); 10 is extract volume of digestive juice (ml); M is powder weight of sample (g).
Quantification of morphological traits
Flag leaf area (LA, area of total leaf, in cm2), flag leaf length (FLL, from base of ligula to tip of leaf, in cm) and flag leaf width (FLW, widest part of leaf, in cm) were measured using LI-3000C Portable Area Meter. Four flag leaves of main spike from each replication were measured.
Statistics, correlation and QTL analyses were performed for the data from each year. Mean value from two years was also used for QTL analysis if genotype × year interaction did not reach significant level for that trait. Homogeneity of variance and normality of distribution were tested before analysis of variance (ANOVA). Heritability was calculated for each trait using ANOVA analysis. The General Linear Model was used for analysis of variance. All analyses were performed using IBM SPSS Statistics 19 software. P value less than 0.05 was considered as significance.
Linkage map was constructed using the software MAPMAKER version 3.0 . Genetic distance (centiMorgans, cM) was derived from Kosambi function. The software MapChart 2.2 was used to draw QTL location on the map.
Total of 153 SSR markers evenly distributed on 7 barley chromosomes were used to construct a barley linkage map as previous described [43,48]. The most likely location of QTL and their genetic effects were detected by composite interval mapping (CIM) using QTL Cartographer version 2.5 . After performing 1000 permutation test, a LOD threshold of 3.0 was used to declare presence of a putative QTL in a given genomic region . Composite interval mapping (CIM) was employed to identify QTL using Model 6 of the Standard module. Cofactors were chosen using the forward-backward method of stepwise regression. The genome was scanned at 2 cM intervals and the window size set at 10 cM. Percentage of phenotypic variation explained and additive effect of each QTL were also calculated by QTL Cartographer 2.5. QTL name was composed of q, the abbreviation of trait, the location of chromosome and the marker position on chromosome.
Quantitative trait locus
Marker assisted selection
Net photosynthetic rate
Intercellular CO2 concentration
Flag leaf area
Flag leaf length
Flag leaf width
Relative water content
Relative chlorophyll content
Leaf nitrogen concentration
Horsley RD, Franckowiak JD, Schwarz PB. Barley. In: Carena MJ, editor. Cereals. US: Springer; 2009. p. 227–50.
Teng S, Qian Q, Zeng D, Kunihiro Y, Fujimoto K, Huang D, et al. QTL analysis of leaf photosynthetic rate and related physiological traits in rice (Oryza sativa L). Euphytica. 2004;135:1–7.
Berdahl JD, Rasmusson DC, Moss DN. Effects of leaf area on photosynthetic rate, light penetration, and grain yield in barley. Crop Sci. 1972;12:177–80.
Flood PJ, Harbinson J, Aarts MG. Natural genetic variation in plant photosynthesis. Trends Plant Sci. 2011;16:327–35.
Sarrafi A, Planchon C, Ecochard R, Sioud Y. Inheritance of some physiological factors of productivity in barley. Genome. 1987;29:846–9.
Tungland L, Chapko LB, Wiersma JV, Rasmusson DC. Effect of erect leaf angle on grain yield in barley. Crop Sci. 1987;27:37–40.
Yap TC, Harvey BL. Relations between grain yield and photosynthetic parts above the flag leaf node in barley. Can J Plant Sci. 1972;52:241–6.
Sicher RC. Assimilate partitioning within leaves of small grain cereals. In: Abrol YP, Mohanty P, Govindjee, editors. Photosynthesis: Photoreactions to Plant Productivity. Netherlands: Springer; 1993. p. 351–60.
Thorne GN. Photosynthesis of ears and flag leaves of wheat and barley. Ann Bot. 1965;29:317–29.
Gyenis L, Yun SJ, Smith KP, Steffenson BJ, Bossolini E, Sanguineti MC, et al. Genetic architecture of quantitative trait loci associated with morphological and agronomic trait differences in a wild by cultivated barley cross. Genome. 2007;50:714–23.
Xue D, Chen M, Zhou M, Chen S, Mao Y, Zhang G. QTL analysis of flag leaf in barley (Hordeum vulgare L.) for morphological traits and chlorophyll content. J Zhejiang Uni Sci B. 2008;9:938–43.
Yin X, Kropff MJ, Stam P. The role of ecophysiological models in QTL analysis: the example of specific leaf area in barley. Heredity. 1999;82:415–21.
Diab AA, Teulat-Merah B, This D, Ozturk NZ, Benscher D, Sorrells ME. Identification of drought-inducible genes and differentially expressed sequence tags in barley. Theo Appl Genet. 2004;109:1417–25.
Guo P, Baum M, Varshney RK, Graner A, Grando S, Ceccarelli S. QTLs for chlorophyll and chlorophyll fluorescence parameters in barley under post-flowering drought. Euphytica. 2008;163:203–14.
Mickelson S, See D, Meyer FD, Garner JP, Foster CR, Blake TK, et al. Mapping of QTL associated with nitrogen storage and remobilization in barley (Hordeum vulgare L.) leaves. J Exp Bot. 2003;54:801–12.
Teulat B, Zoumarou-Wallis N, Rotter B, Salem MB, Bahri H, This D. QTL for relative water content in field-grown barley and their stability across Mediterranean environments. Theor Appl Genet. 2003;108:181–8.
This D, Borries C, Souyris I, Teulat B. QTL study of chlorophyll content as a genetic parameter of drought tolerance in barley. Barley Genet Newsl. 2000;30:20–3.
Wójcik-Jagła M, Rapacz M, Tyrka M, Kościelniak J, Crissy K, Żmuda K. Comparative QTL analysis of early short-time drought tolerance in Polish fodder and malting spring barleys. Theor Appl Genet. 2013;126:3021–34.
Qi X, Stam P, Lindhout P. Comparison and integration of four barley genetic maps. Genome. 1996;39:379–94.
Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, et al. A high-density consensus map of barley linking DArT markers to SSR. RFLP and STS loci and agricultural traits. BMC Genomics. 2006;7:206.
Horton P. Prospects for crop improvement through the genetic manipulation of photosynthesis: morphological and biochemical aspects of light capture. J Exp Bot. 2000;51:475–85.
Rawson HM, Constable GA. Carbon production of sunflower cultivars in field and controlled environments. I. Photosynthesis and transpiration of leaves, stems and heads. Funct Plant Biol. 1980;7:555–73.
Jiang Q, Roche D, Monaco TA, Hole D. Stomatal conductance is a key parameter to assess limitations to photosynthesis and growth potential in barley genotypes. Plant Biol. 2006;8:515–21.
Rybiński W, Garczyński S. Influence of laser light on leaf area and parameters of photosynthetic activity in DH lines of spring barley (Hordeum vulgare L.). Int Agrophys. 2004;18:261–8.
Rhodes MJC, Yemm EW. The development of chloroplasts and photosynthetic activities in young barley leaves. New Phytol. 1966;65:331–42.
Von Korff M, Grando S, Del Greco A, This D, Baum M, Ceccarelli S. Quantitative trait loci associated with adaptation to Mediterranean dryland conditions in barley. Theor Appl Genet. 2008;117:653–69.
Araus JL, Amaro T, Voltas J, Nakkoul H, Nachit MM. Chlorophyll fluorescence as a selection criterion for grain yield in durum wheat under Mediterranean conditions. Field Crops Res. 1998;55:209–23.
Siahsar BA, Aminfar Z. Mapping QTLs of physiological traits associated with salt tolerance in ‘Steptoe’ × ‘Morex’ doubled haploid lines of barley at seedling stage. J Food Agric Environ. 2010;8:751–9.
Yang L, Mickelson S, See D, Blake TK, Fischer AM. Genetic analysis of the function of major leaf proteases in barley (Hordeum vulgare L.) nitrogen remobilization. J Exp Bot. 2004;55:2607–16.
Sinclair TR, Horie T. Leaf nitrogen, photosynthesis, and crop radiation use efficiency: a review. Crop Sci. 1989;29:90–8.
Lewis OAM, James DM, Hewitt EJ. Nitrogen assimilation in barley (Hordeum vulgare L. cv. Mazurka) in response to nitrate and ammonium nutrition. Ann Bot. 1982;49:39–49.
Shangguan Z, Shao M, Dyckmans J. Effects of nitrogen nutrition and water deficit on net photosynthetic rate and chlorophyll fluorescence in winter wheat. J Plant Physiol. 2000;156:46–51.
Sage RF, Pearcy RW. The nitrogen use efficiency of C3 and C4 plants II. Leaf nitrogen effects on the gas exchange characteristics of Chenopodium album (L.) and Amaranthus retroflexus (L.). Plant Physiol. 1987;84:959–63.
Teulat B, Monneveux P, Wery J, Borries C, Souyris I, Charrier A, et al. Relationships between relative water content and growth parameters under water stress in barley: a QTL study. New Phytol. 1997;137:99–107.
Johnson RR, Frey NM, Moss DN. Effect of water stress on photosynthesis and transpiration of flag leaves and spikes of barley and wheat. Crop Sci. 1974;14:728–31.
González A, Martín I, Ayerbe L. Yield and osmotic adjustment capacity of barley under terminal water‐stress conditions. J Agron Crop Sci. 2008;194:81–91.
Zhao J, Sun H, Dai H, Zhang G, Wu F. Difference in response to drought stress among Tibet wild barley genotypes. Euphytica. 2010;172:395–403.
Matin MA, Brown JH, Ferguson H. Leaf water potential, relative water content, and diffusive resistance as screening techniques for drought resistance in barley. Agron J. 1989;81:100–5.
Sinclair TR, Ludlow MM. Who taught plants thermodynamics? The unfulfilled potential of plant water potential. Aust J Plant Physiol. 1985;12:213–7.
Teulat B, This D, Khairallah M, Borries C, Ragot C, Sourdille P, et al. Several QTLs involved in osmotic-adjustment trait variation in barley (Hordeum vulgare L.). Theor Appl Genet. 1998;96:688–98.
Teulat B, Borries C, This D. New QTLs identified for plant water status, water-soluble carbohydrate and osmotic adjustment in a barley population grown in a growth-chamber under two water regimes. Theor Appl Genet. 2001;103:161–70.
Varshney RK, Marcel TC, Ramsay L, Russell J, Röder MS, Stein N, et al. A high density barley microsatellite consensus map with 775 SSR loci. Theor Appl Genet. 2007;114:1091–103.
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.
Murchie EH, Pinto M, Horton P. Agriculture and the new challenges for photosynthesis research. New Phytol. 2009;181:532–52.
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.
Elberse IAM, Vanhala TK, Turin JHB, Stam P, van Damme JMM, van Tienderen PH. Quantitative trait loci affecting growth-related traits in wild barley (Hordeum spontaneum) grown under different levels of nutrient supply. Heredity. 2004;93:22–33.
Forster BP, Ellis RP, Moir J, Talame V, Sanguineti MC, Tuberosa R, et al. Genotype and phenotype associations with drought tolerance in barley tested in North Africa. Ann Appl Biol. 2004;144:157–68.
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.
Cuesta-Marcos A, Casas AM, Hayes PM, Gracia MP, Lasa JM, Ciudad F, et al. Yield QTL affected by heading date in Mediterranean grown barley. Plant Breed. 2009;128:46–53.
Garcı́a del Moral LF, Garcı́a del Moral MB, Molina-Cano JL, Slafer GA. Yield stability and development in two-and six-rowed winter barleys under Mediterranean conditions. Field Crops Res. 2003;81:109–19.
Pourkheirandish M, Wicker T, Stein N, Fujimura T, Komatsuda T. Analysis of the barley chromosome 2 region containing the six-rowed spike gene vrs1 reveals a breakdown of rice–barley micro collinearity by a transposition. Theor Appl Genet. 2007;114:1357–65.
Komatsuda T, Li W, Takaiwa F, Oka S. High resolution map around the vrs1 locus controlling two- and six-rowed spike in barley, Hordeum vulgare. Genome. 1999;42:248–53.
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.
Rapacz M, Kościelniak J, Jurczyk B, Adamska A, Wójcik M. Different patterns of physiological and molecular response to drought in seedlings of malt- and feed-type barleys (Hordeum vulgare). J Agron Crop Sci. 2010;196:9–19.
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, et al. MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics. 1987;1:174–81.
Wang S, Basten CJ, Zeng ZB. Windows QTL Cartographer 2.5. Raleigh, NC: Department of Statistics, North Carolina State Univ; 2007.
Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138:963–71.
This project was supported in part by the National Natural Science Foundation of China (31301310 and 31228017) and the earmarked fund for China Agriculture Research System (CARS-5).
The authors declare that they have no competing interests.
LL performed this study, statistical analysis and manuscript writing. XR assisted in phenotyping and software analysis. DS and GS conceived this study, coordinated the experiments, and wrote the manuscript. CL produced the Huaai 11 and Huadamai 6 DH population. All authors have read and approved the final version of this manuscript.
About this article
Cite this article
Liu, L., Sun, G., Ren, X. et al. Identification of QTL underlying physiological and morphological traits of flag leaf in barley. BMC Genet 16, 29 (2015). https://doi.org/10.1186/s12863-015-0187-y
- Net photosynthetic rate
- Stomatal conductance
- Flag leaf area
- Flag leaf length
- Flag leaf width
- Relative chlorophyll content
- Leaf nitrogen concentration