- Research article
AB-QTL analysis reveals new alleles associated to proline accumulation and leaf wilting under drought stress conditions in barley (Hordeum vulgareL.)
BMC Geneticsvolume 13, Article number: 61 (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.
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 [1–3]. 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 . 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 [5–16]. 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 . Inability of this process may result in leaf wilting and death of leaves because of failure to cope with the transpiration demands of plants . Price et al.  reported eleven QTL for leaf drying which contributed significantly to drought avoidance of upland rice. In sorghum, Sanchez et al.  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 . 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 . 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 . 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, 26–29]. 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.
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).
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.
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).
The phenotypic and genotypic data has been subjected to QTL analysis for the identification of QTL associated to PC and WS (Table 2).
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).
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).
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 . 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  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 .
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 [38–40]. However, in barley leaves the preferential accumulation of proline in epidermis and vascular bundles was only observed under stress conditions , 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 . 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  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 [43–45]. 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  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 . Quarrie et al.  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 . 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.  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.  identified DNA polymorphism in the 3'UTR of HvFT2 between Scarlett and ISR42-8. Von Korff et al.  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.
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. . The proportion of donor genome in this population is given in von Korff et al.  and Schmalenbach et al. .
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.  and Wang et al. , 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 . 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. . A linkage map has been drawn by using MapChart ver.2.2 .
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.
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. . 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 .
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 .
QTL analysis is to a great extent a model selection. Broman and Speed  and Bauer et al.  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 . 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 R2 coefficient (percentage of the explained genotypic variance) according to von Korff et al. . 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 .
Advanced backcross quantitative trait locus
Marker main effect
- M x T:
Marker by treatment interaction effect
- M x E:
Marker by environment interaction effect
Boyer JS: Plant productivity and environment. Science. 1982, 218: 443-448. 10.1126/science.218.4571.443.
Blum A: Plant Breeding for Water-limited Environments. 2011, New York: Springer
Pennisi E: The blue revolution, drop by drop, gene by gene. Science. 2008, 320: 171-173. 10.1126/science.320.5873.171.
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. 10.1104/pp.108.128967.
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.
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. 10.1007/s001220000503.
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. 10.1023/A:1014894130270.
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. 10.1007/s00122-003-1417-7.
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. 10.1111/j.1744-7348.2004.tb00346.x.
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. 10.1007/s00122-005-0144-7.
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.
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. 10.1534/genetics.107.077297.
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. 10.1007/s00122-008-0846-8.
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. 10.1111/j.1558-5646.2008.00474.x.
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. 10.1007/s00122-009-1099-x.
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. 10.1007/s11032-010-9393-5.
Blum A: Plant Breeding for Stress Environments. 1988, Florida. USA: CRC. Inc
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. 10.1023/A:1014805625790.
Clarke JM: Effect of leaf rolling on leaf water loss in Triticum ssp. Can J Plant Sci. 1986, 66: 885-891. 10.4141/cjps86-111.
Bartels D, Sunker R: Drought and salt tolerance in plants. Criti Rev Plant Sci. 2005, 24: 23-58. 10.1080/07352680590910410.
Bohnert HJ, Nelson DE, Jensen RG: Adaptations to environmental stresses. Plant Cell. 1995, 7: 1099-1111.
Sleator RD, Hill C: Bacterial osmoadaptation: the role of osmolytes in bacterial stress and virulence. FEMS Microbiol Rev. 2002, 26: 49-71. 10.1111/j.1574-6976.2002.tb00598.x.
Szabados L, Savouré A: Proline: a multifunctional amino acid. Trend Plant Sci. 2009, 15: 89-97.
Lehmann S, Funck D, Szabados L, Rentsch D: Proline metabolism and transport in plant development. Amino Acids. 2010, 39: 949-962. 10.1007/s00726-010-0525-3.
Iyer S, Caplan A: Products of proline catabolism can induce osmotically regulated genes. Plant Physiol. 1998, 116: 203-211. 10.1104/pp.116.1.203.
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. 10.1007/s11103-007-9269-1.
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. 10.1111/j.1399-3054.2009.01261.x.
Lansac AR, Sullivan CY, Johnson BE: Accumulation of free proline in sorghum (Sorghum bicolor) pollen. Can J Bot. 1996, 74: 40-45. 10.1139/b96-006.
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. 10.1007/s00425-007-0527-x.
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.
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. 10.2135/cropsci1979.0011183X001900040015x.
Bates L, Waldren R, Teare I: Rapid determination of free proline for water stress studies. Plant & Soil. 1973, 39: 205-207. 10.1007/BF00018060.
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.
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. 10.1007/s00122-006-0223-4.
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. 10.1007/BF00223376.
Mather K, Jinks JL: Biometrical Genetics. 1982, London: Chapman and Hall, 3
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. 10.1007/s10681-005-1681-5.
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. 10.1111/j.1365-3040.1995.tb00187.x.
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. 10.1104/pp.118.2.661.
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. 10.1111/j.0031-9317.2004.00251.x.
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. 10.1016/0031-9422(89)80024-X.
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.
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.
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. 10.1016/S0168-9452(98)00175-7.
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. 10.1104/pp.122.4.1129.
Stewart CR: Role of carbohydrates in proline accumulation in wilted barley 1eaves. Plant Physiol. 1978, 61: 775-778. 10.1104/pp.61.5.775.
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. 10.1111/j.1365-313X.2007.03234.x.
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. 10.1023/A:1005864202924.
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.
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. 10.1105/tpc.007096.
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. 10.1016/j.plantsci.2003.11.023.
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. 10.1007/s00122-010-1276-y.
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. 10.1073/pnas.0406232101.
Phillips PC: Epistasis-the essential role of gene interactions in the structure and evolution of genetic systems. Nature Rev Genet. 2008, 9: 855-867. 10.1038/nrg2452.
Voorrips RE: MapChart: software for the graphical representation of linkage maps and QTLs. J of Heredity. 2002, 93: 77-78. 10.1093/jhered/93.1.77.
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. 10.1007/s00122-004-1818-2.
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. 2011.
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. 10.1093/nar/29.1.1.
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-10.1186/1471-2164-7-206.
SAS Institute: The SAS System for Windows, release 9.2. 2008, Cary NC: SAS Institute
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. 10.1111/1467-9868.00354.
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. 10.1007/s00122-009-1021-6.
Benjamini Y, Yekutieli D: Quantitative trait loci analysis using the false discovery rate. Genetics. 2005, 171: 783-790. 10.1534/genetics.104.036699.
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).
The authors declare that they have no competing interests.
JL has conceptualized the research and MAS, HS carried out phenotypic experimentation. JL, KP and MAS contributed in achieving genotyping data. JL has written the program for QTL analysis and JL, MAS and AAN analyzed that data. JL, AAN and MAS have written the manuscript. All authors approved the final version of the manuscript.