AB-QTL analysis reveals new alleles associated to proline accumulation and leaf wilting under drought stress conditions in barley (Hordeum vulgare L.)

  • Mohammed A Sayed1,

    Affiliated with

    • Henrik Schumann1,

      Affiliated with

      • Klaus Pillen2,

        Affiliated with

        • Ali A Naz1 and

          Affiliated with

          • Jens Léon1Email author

            Affiliated with

            BMC Genetics201213:61

            DOI: 10.1186/1471-2156-13-61

            Received: 20 December 2011

            Accepted: 26 June 2012

            Published: 20 July 2012



            Land plants have evolved several measures to maintain their life against abiotic stresses. The accumulation of proline is the most generalized response of plants under drought, heat or salt stress conditions. It is known as an osmoprotectant which also acts as an instant source of energy during drought recovery process. But, both its role and genetic inheritance are poorly understood in agriculture crops. In the present work, advanced backcross quantitative trait locus (AB-QTL) analysis was performed to elucidate genetic mechanisms controlling proline accumulation and leaf wilting in barley under drought stress conditions.


            The analysis revealed eight QTL associated to proline content (PC) and leaf wilting (WS). QTL for PC were localized on chromosome 3H, 4H, 5H and 6H. The strongest QTL effect QPC.S42.5H was detected on chromosome 5H where drought inducible exotic allele was associated to increase PC by 54%. QTL effects QPC.S42.3H, QPC.S42.4H and QPC.S42.6H were responsible to heighten PC due to the preeminence of elite alleles over the exotic alleles which ranged from 26% to 43%. For WS, QTL have been localized on chromosome 1H, 2H, 3H and 4H. Among these, QWS.S42.1H and QWS.S42.4H were associated to decrease in WS due to the introgression of exotic alleles. In addition, two digenic epistatic interaction effects were detected for WS where the additive effect of exotic alleles imparted a favorable increase in the trait value.


            The present data represents a first report on whole-genome mapping of proline accumulation and leaf wilting in barley. The detected QTL are linked to new alleles from both cultivated and wild accessions which bring out an initial insight on the genetic inheritance of PC and WS. These QTL alleles are fixed in the isogenic background of Scarlett, which will allow for positional cloning of underlying genes and to develop drought resilient barley cultivars.


            QTL analysis Drought stress tolerance Proline content Leaf wilting


            Water is fundamental to life and its shortage can cause unprecedented risks to the survival of flora and fauna. Land plants suffer more because of their sessile mode and hence, drought is by far the most devastating abiotic stress that limits agricultural production worldwide [13]. In this scenario, the development of drought-adaptive cultivars is essential to reduce crop losses in agriculture.

            Barley (Hordeum vulgare L.) is the fourth important cereal crop and considered as a genomic model for the tribe Triticeae. It is characterized by immense variation with respect to drought adaptation and has been cultivated from boreal to equatorial regions of the world [4]. In this regard, the wild progenitor of barley, H. vulgare ssp. spontaneum, has a real potential of maintaining a stable population in the desert conditions of Middle East. Genetic dissection of such novel adaptation is a prime question towards understanding drought stress tolerance in plants. The slow progress in this area is mainly due to quantitative inheritance of drought related traits. The advent of molecular tools has made it possible to dissect genetic inheritance of this trait-complex and several successful QTL analyses have been performed in barley and related species [516]. These studies revealed that crop plants evolved a number of drought adaptive traits to maintain their life in water deficit conditions. For instance, leaf rolling is an interesting adaptation to conserve internal water by reducing the transpirational losses under drought stress conditions [5]. Inability of this process may result in leaf wilting and death of leaves because of failure to cope with the transpiration demands of plants [17]. Price et al. [18] reported eleven QTL for leaf drying which contributed significantly to drought avoidance of upland rice. In sorghum, Sanchez et al. [7] identified four major QTL for stay-green, the most important agronomic trait for sorghum cultivation under drought conditions. These data showed that leaf wilting (drying) is a vital drought adaptive trait which offers a straightforward determination of drought tolerance in plants. Therefore, it has been used in large scale screenings of drought tolerance under field conditions [19]. No report of QTL mapping for leaf wilting was found in barley.

            The knowledge created in model species suggested that plants inherit cues of physiological responses which determine or regulate the development of different drought adaptive traits. Among these, accumulation of free proline is a major and the most studied cellular signal of plants under stress conditions. It is primarily synthesized from glutamate via two successive reductions catalyzed by Δl pyrroline-5-carboxylate synthetase (P5CS) and P5C reductase (P5CR), respectively [20]. It acts as a widely distributed osmosolute that protects plants against drought, salinity or low temperature etc. [21, 22]. Its higher reserves can be converted back to glutamate which in turn becomes a prompt source of energy during drought recovery period [23, 24]. It has also been reported that the accumulation of proline products and catabolism of glutamate can result in the expression of several drought inducible genes in rice [25]. Further, its significance has also been reported under control conditions in fulfilling diverse functions during plant development like rapid cell division, floral transition and embryo development [24, 2629]. These reports highlight a diverse role of proline metabolism in plant development under both control and stress conditions in model plants but little is known about its significance in the process of drought stress tolerance in agricultural crops. For instance, its levels were tested to screen drought stress tolerance in barley where higher proline level was found in drought susceptible genotypes [30, 31], which has created a significant confusion about its modulation in the mechanism of drought stress tolerance in crop species. One major reason behind this lacking is the poor understanding of its genetic inheritance in crops.

            The present study was focused on the dissection of the genetic inheritance of proline content (PC) and the associated trait leaf wilting (WS) in barley. We report the first whole-genome QTL map for PC and WS that reveals the identification of associated QTL alleles from both cultivated and exotic origins. Furthermore, the role of stable and drought inducible QTL alleles as well as the digenic epistatic interactions in the determination of QTL effects associated to PC and WS under drought stress conditions has been investigated.


            Variance analysis of PC and WS

            A variance analysis of proline content (PC) and wilting score (WS) among the population S42, across treatments and years, is presented in Table 1. It reveals significant variation for PC among BC2DH lines of the population S42, between treatments, S42 by treatment interaction, across years and block by (treatment x year). For WS, highly significant variation were found among population S42, between treatments, across years and block by (treatment x year) except for the S42 by treatment interaction effect. A lower heritability of PC was observed among the 301 BC2DH lines but these lines displayed heritable responses for WS at h2 =0.75.
            Table 1

            Mutli-factorial analysis of variance of PC and WS among population S42




















            S42 x Treatment












            Block (Treatment x Year)


















            S42 x Treatment












            Block (Treatment x Year)





            1)Traits; PC (proline content) and WS (leaf wilting score), 2)source of variance, 3)degree of freedom, 4)F-value, 5)P-value (<0.05, <0.01, <0.001), 6)heritability of traits among BC2DH population S42.

            Phenotypic characterization of PC and WS

            The parents, Scarlett and ISR42-8 showed significant variation for PC and WS under drought stress and control conditions (Figure 1). Scarlett revealed a remarkable increase of PC from 0.8 μmol/gDW (control) to 4.7 μmol/gDW under drought stress whereas ISR42-8 showed a modest increase in PC in drought block as compared to control. On average, Scarlett accumulated around 3.3 μmol/gDW more PC than ISR42-8 under drought stress conditions. Likewise, Scarlett showed a mean WS of 2.5 under control conditions that increase up to 5.0 under drought stress condition whereas ISR42-8 presented a slight increase in WS from 2.0 (control block) to 3.7 (drought stress block).
            Figure 1

            Comparison of PC (μmol/gDW) and WS (from 0 to 9) between ISR42-8 (ISR) and Scarlett (SCA) under control and drought stress conditions. Trait values have been averaged across years 2007, 2008 and 2009. PC was measured according to Bates et al. [32]. WS was assessed by a visual scoring from 0 to 9 according to de Datta et al. [33].

            The population S42 showed a range of PC values with a mean of 5.9 μmol/gDW under drought stress conditions (Figure 2A). A total of 15 BC2DH lines accumulated the least amount of PC (1.0 μmol/gDW) under drought stress. The largest group of 101 BC2DH lines produced 3.0 μmol/gDW of PC. The second largest group contained 54 BC2DH lines which synthesized around 5 μmol/gDW of PC. Altogether, 131 BC2DH lines showed values of PC that exceeded Scarlett and the mean value of population S42. In these genotypes, the PC ranged from 7.0 to 25.0 μmol/gDW. For WS, population S42 showed a significant variation and presented a mean WS of 5.1 under drought stress treatment (Figure 2B). Two BC2DH lines showed the lowest (resistant) WS at 3.0 while around 40 BC2DH lines displayed WS 4.0. A total of 125 and 96 BC2DH lines showed leaf wilting severity at 5.0 and 6.0 respectively. Around 38 BC2DH lines appeared to be highly susceptible under drought and accounted for WS 7.0 and 8.0.
            Figure 2

            Frequency distribution of population S42 for PC (A) and WS (B) under drought stress condition. The vertical lines represent mean trait value of Scarlett, ISR42-8 and population S42. PC was measured according to Bates et al. [32]. WS was assessed by a visual scoring from 0 to 9 according to de Datta et al. [33].

            To dissect the relationship among PC and WS, Pearson correlation coefficient (r) were calculated which indicated a significant but weak correlation between PC and WS traits (r = 0.2, P < 0.001).

            QTL analysis

            The phenotypic and genotypic data has been subjected to QTL analysis for the identification of QTL associated to PC and WS (Table 2).
            Table 2

            Summary of QTL for PC and WS








            R 2 8) (%)


            RP[Hsp]10) (%)

            Add 11)

            Proline content






            M, M x T


            < 0.01



            - 43

            - 0.77





            M, M x T


            < 0.05



            - 27

            - 0.47





            M x T


            < 0.01









            M, M x T


            < 0.01



            - 26

            - 0.44

            Wilting score





            94.9 - 123.9



            < 0.01



            - 17

            - 0.27




            38.9 - 44.7



            < 0.01








            118.7 - 141.9



            < 0.01








            140.2 - 146



            < 0.01



            - 17

            - 0.18

            1)Description of quantitative trait locus. 2)Chromosome. 3)Linked DNA marker revealing strongest F-value 4)CentiMorgan positions of associated DNA markers from the first to the last significant marker in a QTL region. 5)Effects showing by the QTL, marker main (M) and marker by treatment (M*T) interaction effects. 6)F-value of the given marker locus. 7)Probability of false discovery rate at P < 0.05. 8)Genetic variance explained by M or M x T. For calculation see von Korff et al. [34]. 9)Trait value of homozygous exotic genotypes Hsp in terms of Lsmeans across years.10)Relative performance of the homozygous exotic allele, RPHsp, a “-“indicates a reduction in trait value. 11)The additive effect is half the difference between the phenotypic means of the homozygous elite and exotic marker genotypes. Traits; PC (proline content) and WS (leaf wilting score).

            QTL for PC

            The QTL analysis revealed four QTL for PC located on chromosomes 3H, 4H, 5H and 6H. The strongest QTL effect, QPC.S42.5H was detected on chromosome 5H where an exotic allele accounted for a 54% increase in PC and asserted the highest positive additive effect (0.89). Linked marker to this QTL showed a marker by treatment (M x T) interaction effect and explained 4.1% of the explained genetic variance. The remaining three QTL, QPC.S42.3H, QPC.S42.4H and QPC.S42.6H showed a decreasing trend of PC due to the introgression of exotic alleles. These QTL effects revealed a preponderance of elite alleles over the exotic alleles for PC. At QTL, QPC.S42.3H the elite allele showed 43% increase in PC with respect to the exotic allele and explained 6.1% of the explained genetic variance. Similarly, the relative performance of exotic alleles at QTL, QPC.S42.4H and QPC.S42.6H was 26% and 27% less than the respective elite alleles. Here, each QTL represented around 4% of the genetic variance (Table 2, Figure 3).
            Figure 3

            Marker by treatment interaction effect for PC. Lsmeans of elite (Hv) and exotic (Hsp) alleles for PC under control and drought conditions were compared. Vertical lines indicate standard error of the mean.

            QTL for WS

            Reduced WS under drought is a desirable trait and therefore, QTL effects which accounted for lower WS, represent favorable leads for drought tolerance. QTL analysis revealed four QTL for WS on chromosomes 1H, 2H, 3H and 4H. Among these, two favorable exotic alleles at QTL, QWS.S42.1H and QWS.S42.4H were associated to almost 17% decrease in WS. These exotic alleles explained 12.0% and 9.4% of the explained genetic variance, respectively. In contrast, the exotic alleles at QTL, QWS.S42.2H and QWS.S42.3H were associated with an enhancement of WS as compared to elite alleles. Here, elite alleles appeared to contribute in decreasing WS of which the elite allele at QTL, QWS.S42.2H posed a 22% superior response as compared to the exotic allele and accounted for 5.6% of the explained genetic variance. Likewise, the relative performance of the elite allele at QWS.S42.3H was 35% higher than the respective exotic allele (Table 2).

            Epistatic interaction effects

            Digenic epistatic interactions have been tested for PC and WS among 301 BC2DH genotypes. No epistatic interaction effect was found for PC whereas two interaction effects were detected for WS (Figure 4). The first interaction effect was identified between marker locus bPb5339 (1H) and HvFT2 (3H). At these loci, the combination of Hv/Hv or Hv/Hsp resulted in mean WS 4.4 which reduced to 2.3 as the elite allele was substituted with exotic allele at bPb5339 (1H). In the second effect, a similar interaction of exotic alleles was detected at marker locus bPb-0353 (3H) and Bmac316 (6H) which dropped WS from 4.2 (Hv/Hv) to 1.9 (Hsp/Hsp).
            Figure 4

            Digenic epistatic interaction effects for WS. Lsmeans of four genotypes, Hv/Hv (elite allele at loci 1 and 2), Hv/Hsp (elite allele at locus 1 and exotic allele at locus 2), Hsp/Hv (exotic allele at locus 1 and elite allele at locus 2), Hsp/Hsp (exotic allele at loci 1 and 2). Genotype combinations with significantly different Lsmeans when compared with the control genotype Hv/Hv are marked with asterisks (***P < 0.001, **P < 0.01, *P < 0.05). Vertical lines indicate standard error of the mean.


            The present study reports on the genetic dissection of PC and WS by using 301 BC2DH lines of a cross between barley cultivar Scarlett and wild accession ISR42-8. An advanced backcross population was utilized for a straightforward detection and introgression of favorable exotic alleles in the Scarlett background according to Tanksley and Nelson [35]. Our data showed a significant variation of PC and WS between parents as well as among the BC2DH lines. This population showed skewed frequency distributions for both traits because normal distribution is not expected in a BC2DH population. Mather and Jinks [36] proposed a bias to one of the parent when two parents bear difference for a trait because of successive backcrossing of the recurrent parent. The present work reports on high resolution QTL mapping data by using the biggest double haploid population in barley. It is accepted that the strength of a QTL analysis primarily depends upon the size of mapping population and the density of markers on the genetic map [37].

            QTL analysis identified four QTL for PC where at three QTL the elite alleles from Scarlett were associated to heighten PC. These findings agree with higher PC in the susceptible parent Scarlett which suggests that higher proline levels may define its drought sensitivity. Proline level was used as a metabolic measure to screen drought stress tolerance in barley and higher PC were found in drought susceptible genotypes [30, 31] which agrees with the response of Scarlett in the present study. It has been reported that proline levels vary considerably in different plant organs. It was found higher in reproductive organs than vegetative as well as in the organs bearing endogenously controlled dehydration, e.g. in seeds or pollen [3840]. However, in barley leaves the preferential accumulation of proline in epidermis and vascular bundles was only observed under stress conditions [41], suggesting proline induction a reliable marker for the measurement of drought stress response in barley. Our data accounted for the preeminence of Scarlett's alleles at QTL effect QPC.S42.3H, QPC.S42.4H and QPC.S42.6H to enhance PC as compared to respective exotic alleles. Interestingly, these effects revealed both M and M x T interaction effects thus indicating their role under control and drought stress conditions. The knowledge gained in the model plant Arabidopsis suggests that proline biosynthesis occurs both in normal and stress conditions. Under normal conditions, it is synthesized to maintain the housekeeping function of the cell. Altogether, three enzyme coding genes, P5CS1 and P5CS2 and P5CR have been described to regulate proline biosynthesis in Arabidopsis. Under drought, the expression of P5CS1 and P5CR increased which result in higher proline synthesis in the chloroplast whereas P5CS2 is primarily linked to housekeeping proline synthesis in the cytosole [23, 24]. The map position of P5CS1 and P5CS2 and P5CR are not known in barley and hence the identification of QPC.S42.3H, QPC.S42.4H and QPC.S42.6H on chromosomes 3H, 4H and 6H may provide an initial knowledge of proline biosynthetic loci in barley.

            The strongest QTL effect was detected on chromosome 5H, where an exotic allele enhanced PC by 54%. This exotic allele showed an M x T interaction effect only which suggests that QPC.S42.5H may underlie a drought inducible gene of wild origin possibly similar to transcriptional activity of P5CS1 or P5CR which induced under drought stress conditions [23]. However, it is still an open question whether QPC.S42.5H underlies an exotic variant of P5CS1 or P5CR or a new component of proline biosynthesis in barley. To address this question candidate gene approach will be adapted to test DNA sequence polymorphism and expression differences of barley orthologs in Scarlett and ISR42-8. Alternatively, QTL bearing ILs are available for positional cloning of QPC.S42.5H. A single report was found on the genetic mapping of PC in barley where Siahsar and Narouei [42] reported two QTL for salt stress related to proline accumulation on chromosome 5H by using 72 double haploid lines of a cross between Steptoe and Morex. To see the effect of PC on plant performance we measured plant biomass which revealed major shoot weight QTL were associated to Scarlett's alleles suggesting a possible linkage drag of the exotic alleles on shoot biomass as expected from wild barley accession ISR42-8 being inferior in shoot development with respect to Scarlett (data not shown). Therefore, the effect of PC variation on biomass production may not be visible in such mixed background but candidate IL can be used to test its effect on plant performance. Although, the role of proline is still a debate in crop plants, its utility has been validated in the process of drought stress tolerance in Arabidopsis and rice via transgenic approaches [4345]. Therefore, more efforts are needed to generalize these effects in crops like barley.

            The extent of leaf WS is a straightforward criterion to measure a plant's ability to tolerate water limiting conditions which may be linked to drought inducible metabolites like proline. Our data revealed highly significant but weak correlation of PC and WS. Stewart [46] studied proline accumulation in wilted barley leaves and found that wilting caused 40 fold stimulation of proline biosynthesis in non-starved leaves when compared to starved leaves. He suggested the role of carbohydrates which appeared to supply precursors for proline bio-synthesis. These data indicate that although leaf wilting may stimulate proline biosynthesis, it was unable to modulate the process of its biosynthesis which agrees with the partial dependence of PC and WS. No co-localization of QTL was found which also suggest an independent genetic inheritance of both traits. Our data showed four QTL for WS where at QWS.S42.1H and QWS.S42.4H the introgression of exotic alleles was associated to reduce WS under drought stress conditions. Although these QTL are genetically independent to PC, they may still have significance in modulating proline metabolism which requires lower wilting environment. These findings support the idea of introgressing favorable exotic alleles from a tolerant wild accession ISR42-8 into Scarlett where higher PC may be ceased because of precocious leaf wilting and death under drought stress condition. Abscisic acid (ABA) is one major factor implicated for leaf drying (senescence) under severe drought [47]. Quarrie et al. [48] made a QTL analysis to dissect the genetics behind ABA accumulation and found ten QTL for low and high ABA accumulation in rice under drought stress conditions. In a previous study, they found a QTL for ABA accumulation on 5A of wheat in the region of vernalization gene Vrn1 [49]. Thus, the localization of QWS.S42.4H at marker locus VrnH2 may indicate a putative orthologous region of Vrn1 variant in barley. QWS.S42.1H was linked to HvABAIP where Tondelli et al. [10] described QTL for drought stress tolerance in barley. HvA1 (group 3 LEA protein), which is known as ABA induced barley effector gene, was mapped at the similar region on chromosome 1H whose overexpression confers dehydration tolerance in transgenic rice plants [50, 51]. QTL effects QWS.S42.2H and QWS.S42.3H were associated to superior performance of elite alleles for decrease in WS. These findings suggest that the susceptible parent Scarlett also harbors useful alleles for WS that may underlie valuable components of shoot development and their replacement with detrimental exotic alleles might be a reason behind such variation. Strikingly, all four QTL for WS showed no M x T interaction effect suggesting significant and stable QTL effects across control and drought stress blocks. A possible reason may lie in the genetic control of WS where the associated QTL alleles are active across both treatments. Alternatively, it is possible that QTL were unable to pass the criteria of M x T interaction effects due to error created by different environmental (E) conditions across years as phenotypic evaluations were made in each year in single replication. It is worthwhile to mention that QTL x E (years) was not the focus of the present study because the experiments were carried out inside a tunnel and not under real field conditions. Furthermore, two epistatic interactions were detected where exotic alleles presented an additive role in the development lower WS. Here, the drought sensitive elite allele seems epistatic to drought tolerant exotic alleles which overshadowed the performance of exotic alleles in the determination of drought tolerance. A tolerant WS is evident in the donor parent ISR42-8 of these exotic alleles. However, lower wilting due to the interaction of a flowering gene specific marker HvFT2 on chromosome 3H is an interesting outcome which suggests a putative role of flowering as a regulator in drought stress tolerance. Wang et al. [52] identified DNA polymorphism in the 3'UTR of HvFT2 between Scarlett and ISR42-8. Von Korff et al. [34] studied heading date variation in population S42 and found at most of the ten QTL identified that the introgression of the wild allele was associated with enhanced flowering. They indicated early and heterogeneous flowering habits in wild barley were presumably linked to adaptation in a water limiting environment. These data also highlight the necessity of digenic interactions effect for a high resolution QTL analyses. The role of epistatic interactions has been accepted crucial for the determination of a final phenotype of quantitative traits [53, 54].


            The present study brings out a QTL map that reveals a new insight into the genetic inheritance of PC and WS in barley (Figure 5). It suggests that higher proline accumulation merely defines the extent of drought sensitivity of Scarlett. However, the significance of higher proline content in drought tolerance may be conditioned by its utility during the drought recovery period. Unfortunately, most of drought susceptible genotypes like Scarlett fail to utilize proline reserves efficiently because of early leaf senescence (wilting) and leaf death which results in proline reduction during drought conditions. To complement this lacking natural adaptation ISR42-8 seems suitable since it bears low wilting QTL alleles which may facilitate proline metabolism during drought recovery period. Therefore, these exotic alleles seem important genetic resources to decrease drought sensitivity of cultivated varieties. However, further experiments are needed in a near isogenic background to dissect a precise influence of these QTL effects on plant performance. The above mentioned leads are fixed in an isogenic background and currently introgression lines are available for further phenotyping, marker assisted selection as well as for positional cloning of underlying genes.
            Figure 5

            QTL map for PC and WS in barley. QTL for PC and WS are marked with grey and black bars, respectively. Digenic epistatic interactions have been highlighted with dotted lines where arrow heads indicate associated markers on both sides. Linkage map was drawn using MapChart ver.2.2 where the markers and genetic positions are presented on right and left of the chromosome, respectively [55].


            Plant material

            A double haploid mapping population designated as S42 was used for QTL analysis. It consisted of 301 BC2DH lines which were achieved from 76 original BC2F1 lines of a cross between German spring barley cultivar Scarlett (H. vulgare L.) with the wild accession ISR42-8 (H. vulgare ssp. spontaneum) originating from Israel. The cultivar Scarlett was used as the recurrent parent whereas ISR42-8 was utilized as the donor of the drought related traits. The development of the BC2DH population was according to von Korff et al. [56]. The proportion of donor genome in this population is given in von Korff et al. [56] and Schmalenbach et al. [57].

            Genotyping of population S42

            The population S42 was genotyped with simple sequence repeats (SSRs), diversity array technology (DArT) and gene-specific marker systems. A linkage map of 371 genetic markers has been established that contains 106 SSRs, 255 DArT and 10 gene-specific DNA markers. The SSR markers and gene-specific markers were according to von Korff et al. [56] and Wang et al. [52], respectively. For DArT genotyping, DNA of population S42 and both parents was isolated by using DNeasy 96 Plant Kit (Qiagen, Hilden, Germany). The DNA was sent to Diversity Arrays Technology P/L-Triticarte P/L, Yarralumla ACT, Australia for genotyping [58]. At each locus, homozygous elite or homozygous exotic genotypes have been scored as (Hv) and (Hsp), respectively. The chromosomal positions of the DArT markers are according to Wenzl et al. [59]. A linkage map has been drawn by using MapChart ver.2.2 [55].

            Experimental setup

            Phenotypic characterization of the population S42 was repeated three times during the summer seasons 2007, 2008 and 2009 inside the plastic foliar tunnels at the Institute of Crop Science and Resource Conservation, Faculty of Agriculture, University of Bonn. A total of 12 seed of individual BC2DH line, Scarlett and ISR42-8 were sown in 22 x 22 x 26 cm plastic pots containing a mixture of top soil, silica sand, milled lava and peat dust (Terrasoil®, Cordel & Sohn, Salm, Germany). The pots were arranged in a split-plot design where individual BC2DH lines and parents were assigned randomly to the two irrigation blocks with the control and the drought treatments. The drought stress treatments were applied at growth stage ranging from 28 to 32 (BBCH) among BC2DH lines across 2007, 2008 and 2009.

            Watering was done with a drip irrigation system (Netafim ™, Adelaide, Australia) by watering each pot three times a day with a calculated amount of water. Echo2 sensors (Decagon Dev., Pullman, WA, USA) were used to determine the volumetric soil moisture content (VMC) digitally with the frequency domain technique. For this purpose, before sowing eight sensors were installed in the pots of the reference cultivar Scarlett within a depth of 10 cm, four in control pots and four in the pots of the drought treatment where data were collected manually and controlled every day. Comparing these moisture data with the mass and the developmental stage and the transpiration conditions for the next three days from the weather forecast, further irrigation was determined in such a way, that in the control block the VMC remains at field capacity and in the drought treatment the VMC follows a planned curve from the VMC of field capacity (40%) down to the VMC of near wilting point (15%). The plants were fertilized three times per season with 250 ml of NPK liquid fertilizer containing 7% N, 3% P2O5 and 6% K2O. The plants were sprayed against fungi and insects as recommended for barley cultivation if needed. Environmental conditions across 2007, 2008 and 2009, are presented in Table 3.
            Table 3

            The average temperature and relative humidity across year 2007, 2008 and 2009 at Poppelsdorf field station Bonn, Germany






            Temp (°C)

            RH (%)

            Temp (°C)

            RH (%)

            Temp (°C)

            RH (%)




































            * The values showed average temperature (°C) and relative humidity (RH) across each month in 2007, 2008 and 2009.

            Phenotypic evaluation

            For the phenotypic evaluation of the population S42 proline content (PC) and wilting score (WS) were measured. For PC, a total of eight fully expanded leaves of two main tillers were harvested (four leaves each plant) separately in each control and drought stress blocks at the end of the drought period (growth stage around BBCH 41–49). The leaves were cut and wrapped in plastic foil, then frozen in liquid nitrogen, freeze-dried and ground to a fine powder by the help of a milling device (Retsch, Hann, Germany). PC was measured by colorimetric procedure according to Bates et al. [32]. For this, a total of 30 mg of the powder was homogenized in 2 ml of 3% sulphosalicylic acid by vortexing 3 times for 15 sec and centrifuged at 14000 rpm for 10 min. A total of 500 μl of the supernatant was taken in a falcon tube and increased to a volume of 1000 μl by adding 3% sulphosalicylic acid. Later, 1000 μl of ninhydrin acid and 1000 μl of glacial acetic acid were added. Then falcon tubes were vortexed for 15 sec and the resulting mixture heated to 100°C for 1 hour in a water bath. The reaction was stopped by placing tubes in an ice bath. Afterwards, 2 ml of toluene was added to each tube and vortexed for 20 sec. Subsequently, tubes were kept at room temperature for 5 min for a phase separation and the absorbance of chromophore containing toluene (blank) was read at 520 nm using a spectrophotometer. The proline concentration was determined with a standard curve method, using known concentrations of L-proline (P-0380, Sigma-Aldrich®) and calculated on a dry weight basis by the formula:

            For WS, cumulative response of 12 plants of each genotype for leaf wilting (WS) was scored at the end of drought period by visual rating from 0 up to 9 where 0 represents all leaves green and 9 all leaves apparently dried [33].

            Statistical analysis

            Statistical analysis was performed by using SAS software where analysis of variance (ANOVA), heritability (h2) and least square means (Lsmeans) were calculated using PROC MIXED procedure [60].

            QTL analysis

            QTL analysis is to a great extent a model selection. Broman and Speed [61] and Bauer et al. [62] compared the forward selection strategy in REML analysis with Bayesian techniques and found that forward selection strategy is very effective to detect QTL associated to quantitative traits. We employed a multiple QTL model iteratively extended and reduced by forward selection and backward elimination, respectively using the PROC MIXED procedure in SAS [60]. This QTL model bears the ability to utilize individual observations of each trait value simultaneously across year, blocks and therefore, trait values were not averaged across years for marker trait analysis. For QTL analysis with fixed marker effects, it compares the marker variances and the error variance in the F-test. These effects showed a normal distribution there and hence we did not use any data transformation. According to this model, in each round of the forward selection (or backward elimination as the QTL can be excluded also) process, the most informative marker was added as a fixed factor (QTL) into the model and then all remaining markers were scanned with the respective model containing the previously found QTL. In this multiple QTL model, the following iterations were continued until no more additional QTL could be detected. Starting point was the following mixed hierarchical model:
            where the total phenotypic value was the sum of general mean μ, fixed effect M i of the i-th marker genotype, random effect L j (M i ) of the j-th DH line nested in the i-th marker genotype, fixed effect T k of the k-th treatment, random interaction effect L j *T k of the j-th DH line and the k-th treatment, fixed interaction effect M i *T k of the i-th marker genotype and the k-th treatment, fixed effect Y l of the l-th year, fixed interaction effect T k *Y l of the k-th treatment and l-th year, random effect B m (T k *Y l ) of m-th block nested in treatment and years, residue ε n(ijklm) of X ijklmn . P values from F-tests were adjusted genome-wide across all single marker tests using the probability of false discovery rate (PFDR), implemented in the SAS procedure PROC MULTTEST. The significant marker main effects as well as marker by treatment interaction with PFDR < 0.05 were accepted as putative QTL for the next iteration. The final model was:

            where ∑QTL represents the detected QTL from the forward/backward selection process.The contribution of a QTL to trait genotypic variance was estimated by the R 2 coefficient (percentage of the explained genotypic variance) according to von Korff et al. [34]. Relative performance of an exotic allele RPHsp to the corresponding elite allele was used as a measure of trait improvement that was calculated by the formula (RPHsp = Hsp - Hv/Hv * 100). Hsp and Hv represent the Lsmeans of trait value of the homozygous exotic and elite genotypes, respectively.

            Digenic epistatic effects

            The digenic epistatic interactions between marker pairs were tested in 301 BC2DH lines with the SAS procedure PROC MIXED using the following mixed hierarchical model:where M1 i and M2 j are the fixed effects of the i-th and j-th markers, respectively. M1 i *M2 j is the fixed interaction effect of the i-th M1 genotype with j-th M2 genotype. L k (M1 i *M2 j ) is the random interaction effect of i-th M1 and j-th M2 markers nested in k-th BC2DH line. Epistatic effects were accepted based on probability of false discovery rate (PFDR < 0.05) and has been calculated by PROC MULTTEST procedure in SAS [63].



            Advanced backcross quantitative trait locus


            Proline content


            Wilting score


            Marker main effect

            M x T: 

            Marker by treatment interaction effect

            M x E: 

            Marker by environment interaction effect


            Relative performance.



            We offer special thanks to Anne Reinders for her support in the successful completion of this research project. We are also grateful to Dr. Agim Ballvora for his valuable comments and Karin Woitol and Arifuzzaman Md. for correcting the manuscript. A sincere gratitude to the Ministry of Higher Education, Egypt for the award of PhD scholarship. Part of this work was funded by the German Plant Genome Research Initiative (GABI) of the Federal Ministry of Education and Research (BMBF, project 0312278A).

            Authors’ Affiliations

            Institute of Crop Science and Resource Conservation, Crop Genetics and Biotechnology Unit, University of Bonn
            Institute for Agricultural and Nutritional Sciences, Plant Breeding, Martin-Luther-University Halle-Wittenberg


            1. Boyer JS: Plant productivity and environment. Science 1982, 218:443–448.PubMedView Article
            2. Blum A: Plant Breeding for Water-limited Environments. New York: Springer; 2011.View Article
            3. Pennisi E: The blue revolution, drop by drop, gene by gene. Science 2008, 320:171–173.PubMedView Article
            4. Schulte D, Close TJ, Graner A, Langridge P, Matsumoto T, Muehlbauer G, Sato K, Schulman AH, Waugh R, Wise RP: The international barley sequencing consortium - at the threshold of efficient access to the barley genome. Plant Physiol 2009, 149:142–147.PubMedView Article
            5. Champoux MC, Wang G, Sarkarung S, Mackill DJ, O’Toole JC, Huang N, McCouch SR: Locating genes associated with root morphology and drought avoidance in rice via linkage to molecular markers. Theor Appl Genet 1995, 90:969–981.View Article
            6. 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–170.View Article
            7. Sanchez AC, Subudhi PK, Rosenow DT, Nguyen HT: Mapping QTLs associated with drought resistance in sorghum (Sorghum bicolor L. Moench). Plant Mol Biol 2002, 48:713–726.PubMedView Article
            8. Teulat B, Zoumarou-Wallia N, Rotter B, Ben Salem M, 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–188.PubMedView Article
            9. Talamé V, Sanguineti MC, Chiapparino E, Bahri H, Ben Salem M, Forster BP, Ellis RP, Rhouma S, Zoumarou W, Waugh R, Tuberosa R: Identification of Hordeum spontaneum QTL alleles improving field performance of barley grown under rainfed conditions. Ann Appl Biol 2004, 144:309–319.View Article
            10. Tondelli A, Francia E, Barabaschi D, Aprile A, Skinner JS, Stockinger EJ, Stanca AM, Pecchioni N: Mapping regulatory genes as candidates for cold and drought stress tolerance in barley. Theor Appl Genet 2006, 112:445–454.PubMedView Article
            11. Quarrie SA, Quarrie SP, Radosevic R: Dissecting a wheat QTL for yield present in a range of environments: From the QTL to candidate genes. J Exp Bot 2006, 11:2627–2637.View Article
            12. Maccaferri M, Sanguineti MC, Corneti S, Ortega JLA, Ben Salem M, Bort J, DeAmbrogio E, del Moral LFG, Demontis A, El-Ahmed A, Maalouf F, Machlab H, Martos V, Moragues M, Motawaj J, Nachit M, Nserallah N, Ouabbou H, Royo C, Slama A, Tuberosa R: Quantitative trait loci for grain yield and adaptation of durum wheat (Triticum durum Desf.) across a wide range of water availability. Genetics 2008, 178:489–511.PubMedView Article
            13. Mathews KL, Malosetti M, Chapman S, McIntyre L, Reynolds M, Shorter R, van Eeuwijk F: Multi-environment QTL mixed models for drought stress adaptation in wheat. Theor Appl Genet 2008, 117:1077–1091.PubMedView Article
            14. McKay JK, Richards JH, Nemali KS, Sen S, Mitchell-Olds T, Boles S, Stahl EA, Wayne T, Juenger TE: Genetics of drought adaptation in Arabidopsis thaliana II. QTL analysis of a new mapping population, Kas-1 x Tsu-1. Evolution 2008, 62:3014–3026.PubMedView Article
            15. Messmer R, Fracheboud Y, Bänziger M, Vargas M, Ribaut JM SP: Drought stress and tropical maize: QTL-by-environment interactions and stability of QTLs across environments for yield components and secondary traits. Theor Appl Genet 2009, 119:913–930.PubMedView Article
            16. El Soda M, Nadakuduti SS, Pillen K, Uptmoor R: Stability parameter and genotype mean estimates for drought stress effects on root and shoot growth of wild barley pre-introgression lines. Mol Breeding 2010, 26:583–593.View Article
            17. Blum A: Plant Breeding for Stress Environments. Florida. USA: CRC. Inc; 1988.
            18. Price AH, Townend J, Jones MP, Audebert A, Courtois B: Mapping QTLs associated with drought avoidance in upland rice grown in the Philippines and West Africa. Plant Mol Biol 2002, 48:683–695.PubMedView Article
            19. Clarke JM: Effect of leaf rolling on leaf water loss in Triticum ssp. Can J Plant Sci 1986, 66:885–891.View Article
            20. Bartels D, Sunker R: Drought and salt tolerance in plants. Criti Rev Plant Sci 2005, 24:23–58.View Article
            21. Bohnert HJ, Nelson DE, Jensen RG: Adaptations to environmental stresses. Plant Cell 1995, 7:1099–1111.PubMed
            22. Sleator RD, Hill C: Bacterial osmoadaptation: the role of osmolytes in bacterial stress and virulence. FEMS Microbiol Rev 2002, 26:49–71.PubMedView Article
            23. Szabados L, Savouré A: Proline: a multifunctional amino acid. Trend Plant Sci 2009, 15:89–97.View Article
            24. Lehmann S, Funck D, Szabados L, Rentsch D: Proline metabolism and transport in plant development. Amino Acids 2010, 39:949–962.PubMedView Article
            25. Iyer S, Caplan A: Products of proline catabolism can induce osmotically regulated genes. Plant Physiol 1998, 116:203–211.View Article
            26. Mattioli R, Marchese D, D’Angeli S, Altamura M, Costantino P, Trovato M: Modulation of intracellular proline levels affects flowering time and inflorescence architecture in Arabidopsis. Plant Mol Biol 2008, 66:277–288.PubMedView Article
            27. Mattioli R, Falasca G, Sabatini S, Altamura MM, Costantino P, Trovato M: The proline biosynthetic genes P5CS1 and P5CS2 play overlapping roles in Arabidopsis flower transition but not in embryo development. Physiol Plantarum 2009, 137:72–85.View Article
            28. Lansac AR, Sullivan CY, Johnson BE: Accumulation of free proline in sorghum ( Sorghum bicolor ) pollen. Can J Bot 1996, 74:40–45.View Article
            29. Schmidt R, Stransky H, Koch W: The amino acid permease AAP8 is important for early seed development in Arabidopsis thaliana. Planta 2007, 226:805–813.PubMedView Article
            30. Singh TN, Aspinall D, Paleg LG: Proline accumulation and varietal adaptability to drought in barley: a potential metabolic measure of drought resistance. Nature New Biol 1972, 236:188–190.PubMedView Article
            31. Hanson AD, Nelsen CE, Pedersen AR, Everson EH: Capacity for proline accumulation during water stress in barley and its implications for breeding for drought resistance. Crop Sci 1979, 19:489–493.View Article
            32. Bates L, Waldren R, Teare I: Rapid determination of free proline for water stress studies. Plant & Soil 1973, 39:205–207.View Article
            33. de Datta SK, Malabuyouc JA, Aragon EL: A field screening technique for evaluating rice germplasm for drought tolerance during the vegetative stage. Field Crop Res 1998, 19:123–134.View Article
            34. von Korff M, Wang H, Léon J, Pillen K: AB-QTL analysis in spring barley: II. Detection of favourable exotic alleles for agronomic traits introgressed from wild barley (H. vulgare ssp. spontaneum). Theor Appl Genet 2006, 112:1221–1231.PubMedView Article
            35. Tanksley SD, Nelson JC: Advanced backcross QTL analysis: a method for simultaneous discovery and transfer of valuable QTL from unadapted germplasm into elite breeding lines. Theor Appl Genet 1996, 92:191–203.View Article
            36. Mather K, Jinks JL: Biometrical Genetics. 3rd edition. London: Chapman and Hall; 1982.
            37. Collard BCY, Jahufer MZZ, Brouwer JB, Pang ECK: An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. Euphytica 2005, 142:169–196.View Article
            38. Chiang H, Dandekar AM: Regulation of proline accumulation in Arabidopsis thaliana (L.) Heynh during development and in response to desiccation. Plant Cell Environ 1995, 18:1280–1290.View Article
            39. Fujita T, Maggio A, Garcfa-Rfos M, Bressan RA, Csonka LN: Comparative analysis of the regulation of expression and structures of two evolutionarily divergent genes for Ll',-pyrroline-5-carboxylate synthetase from tomato. Plant Physiol 1998, 118:661–674.PubMedView Article
            40. Armengaud P, Thiery L, Buhot N, Grenier-de March G, Savoure A: Transcriptional regulation of proline biosynthesis in Medicc/go Iruncalu[a reveals developmental and environmental specific features. Physiol Plant 2004, 120:442–450.PubMedView Article
            41. Zúñiga G, Argandoiia VH, Corcuera LJ: Distribution of glycine-betaine and proline in water stressed and unstressed barley leaves. Phytochemistry 1989, 28:419–420.View Article
            42. 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, 8:751–759.
            43. KaviKishor PB, Hong Z, Miao GH, Hu CAA, Verma DPS: Over-expression of delta-pyrroline-5-carboxylate synthetase increases proline production and confers osmotolerance in transgenic plants. Plant Physiol 1995, 108:1387–1394.
            44. Zhu BC, Su J, Chan MC, Verma DPS, Fan YL, Wu R: Over-expression of a delta-pyrroline-5-carboxylate synthetase gene and analysis of tolerance to water-stress and salt-stress in transgenic rice. Plant Sci 1998, 139:41–48.View Article
            45. Hong Z, Lakkineni K, Zhang Z, Verma DPS: Removal of feedback inhibition of delta-pyrroline-5-carboxylate synthetase results in increased proline accumulation and protection of plants from osmotic stress. Plant Physiol 2000, 122:1129–1136.PubMedView Article
            46. Stewart CR: Role of carbohydrates in proline accumulation in wilted barley 1eaves. Plant Physiol 1978, 61:775–778.PubMedView Article
            47. Christmann A, Elmar WW, Ernst S, Erwin G: A hydraulic signal in root-to-shoot signalling of water shortage. Plant J 2007, 52:167–174.PubMedView Article
            48. Quarrie SA, Laurie DA, Zhu J, Lebreton C, Semikhodskii A, Steed A, Witsenboer H, Calestani C: QTL analysis to study the association between leaf size and abscisic acid accumulation in droughted rice leaves and comparisons across cereals. Plant Mol Biol 1997, 35:155–165.PubMedView Article
            49. Quarrie SA, Gulli M, Calestani C, Steed A, Marmiroli N: Location of a gene regulating drought-induced abscisic acid production on the long arm of chromosome 5A of wheat. Theor Appl Genet 1994, 89:794–800.View Article
            50. Casaretto J, Ho TD: The transcription factors HvABI5 and HvVP1 are required for the abscisic acid induction of gene expression in barley aleurone cells. Plant Cell 2003, 15:271–284.PubMedView Article
            51. Chandra Babu R, Zhang JA, Blum JA, Ho T-HD, Wu R, Nguyen HT: HVA1, a LEA gene from barley confers dehydration tolerance in transgenic rice (Oryza sativa L.) via cell membrane protection. Plant Sci 2004, 166:855–862.View Article
            52. Wang G, Schmalenbach I, von Korff M, Lèon J, Kilian B, Rode J, Pillen K: Association of barley photoperiod and vernalization genes with QTLs for flowering time and agronomic traits in a BC2DH population and a set of wild barley introgression lines. Theor Appl Genet 2010, 120:1559–1574.PubMedView Article
            53. Caicedo AL, Stinchcombe JR, Olsen KM, Schmitt J, Purugganan MD: Epistatic interaction between Arabidopsis FRI and FLC flowering time genes generates a latitudinal cline in a life history trait. Proc Natl Acad Sci USA 2004, 101:15670–15675.PubMedView Article
            54. Phillips PC: Epistasis-the essential role of gene interactions in the structure and evolution of genetic systems. Nature Rev Genet 2008, 9:855–867.PubMedView Article
            55. Voorrips RE: MapChart: software for the graphical representation of linkage maps and QTLs. J of Heredity 2002, 93:77–78.View Article
            56. von Korff M, Wang H, Léon J, Pillen K: Development of candidate introgression lines using an exotic barley accession (H. vulgare ssp. spontaneum) as donor. Theor Appl Genet 2004, 109:1736–1745.PubMedView Article
            57. Schmalenbach I, March TJ, Bringezu T, Waugh R, Pillen K: High-resolution genotyping of wild barley introgression lines and fine-mapping of the threshability locus thresh-1 using the Illumina GoldenGate assay. G3 2011, 1:187–196.PubMedView Article
            58. Jaccoud D, Peng KM, Feinstein D, Kilian A: Diversity arrays: a solid state technology for sequence information independent genotyping. Nucleic Acids Res 2001, 29:1–7.View Article
            59. Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, Hearnden P, Maier C, Xia L, Caig V, Ovesná J, Cakir M, Poulsen D, Wang J, Raman R, Smith KP, Muehlbauer GJ, Chalmers KJ, Kleinhofs A, Huttner E, Kilian A: A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and agricultural traits. BMC Genomics 2006, 7:206.PubMedView Article
            60. SAS Institute: The SAS System for Windows, release 9.2. Cary NC: SAS Institute; 2008.
            61. Broman KW, Speed TP: A model selection approach for identification of quantitative trait loci in experimental crosses. J R Stat Soc B 2002, 64:641–656.View Article
            62. Bauer AM, Hoti F, von Korff M, Pillen K, Léon J, Sillanpää MJ: Advanced backcross QTL analysis in spring barley (H. vulgare ssp. spontaneum) comparing a REML versus a Bayesian model in multi-environmental field trials. Theor Appl Genet 2009, 119:105–123.PubMedView Article
            63. Benjamini Y, Yekutieli D: Quantitative trait loci analysis using the false discovery rate. Genetics 2005, 171:783–790.PubMedView Article


            © Sayed et al.; licensee BioMed Central Ltd. 2012

            This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.