BACKGROUND
barley  is one of the first cereals domesticated in the fertile crescent  <cit>  and today it is the fourth most important crop species concerning acreage next to wheat, maize, and rice  <cit> .
worldwide, barley is mainly used for animal feed and malting and only a very small amount is used for direct human consumption and bakery.
average yield of barley on the worldwide level is  <dig>  t/ha but in some european countries, e.g.
germany average yield is up to  <dig>  t/ha  <cit> .
barley yield in many parts of the world is reduced by biotic stress but also by abiotic stress e.g.
heat, salt, deficits in nitrogen nutrition and drought .
especially, in the juvenile stages from sowing to tillering, drought can severely influence barley development already reducing the potential yield  <cit> .
research on drought stress tolerance has become more important worldwide as due to climate change the number of drought periods will increase in the future  <cit> .
up to now, most studies conducted in barley focused on effects of terminal drought stress whereas drought in juvenile stages is less well documented  <cit> .
drought tolerance is a complex quantitative trait, that is controlled by various mechanisms  <cit> .
abscisic acid  is a key phytohormone involved in adaption to environmental stresses and regulation of plant development.
it promotes the closure of stomata under drought stress conditions initiated by a loss of turgor  <cit> .
furthermore, it increases the hydraulic conductivity of water, promotes chlorophyll breakdown and leads to leaf senescence  <cit> .
another relevant protein is ubiquitin which regulates the degradation of proteins  <cit> .
moreover, late embryo abundant  proteins and heat shock proteins that are involved in the protection of functional proteins are induced in response to various abiotic stresses .
for example protein kinases and protein phosphatases which activate or deactivate proteins by phosphorylation and dephosphorilation  <cit> , as well as the lea protein dehydrin which is described to have different functions in different stresses  <cit>  is often present along with drought stress  <cit> .
a lot of parameters indicative for drought stress influenced by these and additional genes were analysed in different crops  <cit> .
for example, biomass production  <cit> , yield  <cit> , photosynthesis rate  <cit> , as well as the content of free proline  <cit> , total content of soluble sugars  <cit> , or osmolality  <cit>  are parameters which are affected by drought stress in barley.
another factor relevant for yield improvement is leaf senescence  <cit> , which is a natural degradation process at the final stage of the development of organs and plants.
this process is divided into three steps and starts with reprogramming of gene expression to turn on senescence activating genes.
before programmed cell death in the terminal phase occurs, the second step in which nutrients and metabolites are transported from source  to sink  is important for yield and quality of the seeds harvested  <cit> .
during this reorganization phase, degradation of chlorophyll and a decrease in photosynthesis is observed  <cit> .
because of degradation of chlorophyll, yellowing of the leaves is a symptom of leaf senescence  <cit> , which in many studies is rated visually  <cit> , but can be more precisely determined by soil plant analysis development  readings which estimate leaf greenness  <cit> .
degradation of chlorophyll is regulated by chlorophyllase, pheophorbide α oxygenase and red chlorophyll catabolite reductase among others  <cit> , but so far for regulation of leaf senescence only a few genes are known  <cit> .
leaf senescence is a process which is influenced by a lot of external stress conditions e.g.
drought stress  <cit> .
stress often results in premature induction of leaf senescence and therefore leads to an inefficient recycling of resources and a massive yield loss .
in contrast, plants showing delayed leaf senescence under stress, represented by a “stay green effect”, minimize yield loss  <cit> .
genome wide association studies  are a powerful tool to subdivide such complex pathways as drought stress and leaf senescence by the detection of quantitative trait loci  out of the regression analysis of genotypic and phenotypic data .
up to now, some qtls involved in drought stress response were published in barley , whereas for leaf senescence only few qtls are known  <cit> .
besides this, qtls which are involved in the response to other abiotic stresses, e.g.
salt stress were identified in barley  <cit> .
molecular markers, such as single nucleotide polymorphisms  flanking qtls having a significant influence on the respective trait can be used for efficient marker assisted selection and smart breeding procedures  <cit> .
the aim of the present study is therefore the identification of qtls for drought stress and drought stress induced leaf senescence in early developmental stages of barley suited to be used in future barley breeding programs using gwas followed by the identification of the function of these qtls associated to respective traits.
material & methods
plant material and experimental setup
a set of  <dig> winter barley genotypes  consisting of  <dig> german winter barley cultivars [ <dig> two-rowed and  <dig> six-rowed,  <cit>  and  <dig> accessions of the spanish barley core collection   <cit>  were used to investigate drought stress induced leaf senescence in juvenile barley plants.
drought stress was applied in greenhouses of the julius kühn-institut in groß lüsewitz, germany according to honsdorf et al.
<cit> .
trials were conducted in a split plot design with three replications per genotype and variant .
ten seeds of each accession were sown per plastic pot  containing  <dig>  g of a mixed clay soil ed <dig> .
after germination, seedlings were reduced to seven plants per pot.
plants were grown under semi-controlled long day conditions in a temperature range from  <dig> to 22 °c at day  and  <dig> to 19 °c at night .
if natural radiation was below  <dig> klx, additional light was applied from 6 a.m. to 10 p.m. drought stress started at the primary leaf stage  seven days after sowing .
at this time watering of the stress variant was stopped till the soil reached 20 % of the maximal soil water capacity, and then this level was kept by weighting each pot and re-watering.
control plants were continuously watered to 70 % of the maximal soil water capacity.
water capacity was calculated of the saturated soil weight and drought weight according to din iso  <dig> 1996– <dig>  <cit> .
at the end of a four weeks stress period  physiological traits were determined and above ground biomass was harvested .
experimental setup a was repeated in three years.
a modified experimental setup was conducted to optimize lightening conditions of the primary leaves for the measurement of leaf colour  and electron transport rate at photosystem ii  .
general settings were the same as in setup a, but only four plants were grown in smaller pots  to allow wider spacing of pots and all leaves except the primary leaves were tied up to reduce shading.
experimental setup b was repeated in two years.
physiological parameters determined
six physiological traits, i.e.
leaf colour , electron transport rate at photosystem ii , content of free proline , content of soluble sugars , osmolality  and the aboveground biomass yield  were determined in the control and stress treatment.
measurement and sampling respectively were conducted on primary leaves.
chlorophyll content which was used as the main indicator for drought stress induced leaf senescence was measured 33– <dig> das by minolta spad readings , which gives a value for leaf colour.
three primary leaves of three plants for each pot were measured at five positions per leaf.
these spad readings turned out to be correlated to the chlorophyll content analysed photometrically  <cit> .
because of this relation, the chlorophyll content can be indirectly measured by spad  <cit> .
at 34– <dig> das chlorophyll fluorescence was measured in all genotypes using light adapted plants with the os1p-chlorophyll fluorometer  in the middle of three primary leaves per pot at one position per leaf.
the relative electron transport rate at photosystem ii   * par *  <dig>  84 *  <dig>  5) was calculated including the photosynthetically active radiation , as well as the quantum photosynthetic yield of psii ) and constants representing light which is absorbed by the leaf  and light which is equally absorbed by psi and psii   <cit> .
at  <dig> das five primary leaves per pot were harvested and cut in pieces of 1 cm length for the analysis of cfp, css and oa.
these samples were frozen in liquid nitrogen immediately and samples for cfp and css measurement were freeze dried.
for cfp measurement the ninhydrin method  <cit>  was applied, and for css measurement the anthron method  <cit>  was used.
both traits were measured photometrical using a spectrophotometer.
the concentration of these ingredients was determined with a standard curve calculated on a dry weight basis.
to assess oa, frozen leaf samples were grinded in a swing mill , filled up with 200 μl water and centrifuged at  <dig>  rpm for 15 min to get cell sap for the measurement of osmolality with a freezing-point osmometer .
osmolality was corrected for the water content of fresh and dry weight.
above ground biomass was harvested  <dig> das, too.
leaf material was dried in a compartment dryer at 105 °c and weighted.
for all traits an outlier test was calculated to exclude extreme deviations  <cit> .
to get information on the stability of all analysed traits in response to drought stress compared to the control, the drought susceptibility index  was calculated  <cit>  for each trait and across treatments according to the formula: dsi=1–lsmeantraitgenotype,stress/lsmeantraitgenotype,control/d,with  d=1–lsmeantraitassortment,stress/lsmeantraitassortment,control
dsi is a relative value estimated for each genotype and trait.
according to the formula genotypes revealing a dsi close to one are highly susceptible to drought and those close to zero or showing a negative value are tolerant.
statistical analysis of phenotypic data
statistical analyses were performed with sas  <dig>   <cit> .
least square means  were calculated with glm procedure for the replications of each genotype in the respective years and for both control and drought conditions.
descriptive statistics was calculated out of lsmeans by proc univariate.
analysis of variance  was calculated using proc mixed with genotypes , drought stress treatment  and gxt as fixed effects.
replication , year  and row type are chosen as random factors.
the heritability  was calculated with sas in two steps.
first, the variance components for the genotypes , variance associated with the genotype by year interaction  and ve which is the error variance were calculated with proc varcomp.
next, h <dig> was calculated with the following formula: vg/ for both well watered and drought stress conditions.
furthermore, the coefficient of correlation  by pearson was calculated with sas based on lsmeans.
genotyping and genome wide association study
for genotyping the whole set of genotypes was analysed with the barley illumina 9 k iselect snp-chip  <cit> .
population structure was calculated with structure  <dig> .
<dig>  <cit>  based on  <dig> simple sequence repeat  markers covering the whole genome.
the structure programme was run  <dig> times for pre-defined k  from  <dig> to  <dig> each.
to get the number of calculated subpopulations  with highest likelihood the procedure of evanno et al.
<cit>  was applied.
an independent run with  <dig>  iterations of a monte carlo markov chain with a length of the burn in period of  <dig>  was conducted for the k with the highest likelihood to obtain the q-matrix.
kinship was calculated with spagedi  <dig> d  <cit>  based on  <dig> ssrs and allele size correlation coefficient  <cit>  with  <dig>  permutations.
out of  <dig>  genetically mapped snp markers   <cit> ,  <dig>  polymorphic markers with minor allele frequencies higher than 5 % were taken into account.
based on these data and respective phenotypic data  gwas was conducted applying a mixed linear model  using tassel  <dig>   <cit> .
all results with p values < <dig>  were considered as significant marker trait associations.
linkage disequilibrium  was calculated on mapped polymorphic snps with r  <cit>  by an estimate of the average decay  <cit>  over all barley chromosomes.
sequences of significantly associated snp-markers  were downloaded from the james hutton institute  and respective sequences were compared against the plant proteome in the uniprotkb/swiss-prot protein database by blastx  to get information on the proteins coded by these sequences  <cit> .
using uniprot  <cit>  the involvement of respective proteins in drought stress and leaf senescence processes was analysed.
in a last step a genetic map with all significantly associated snps in genes coding for proteins known to be involved in drought stress tolerance was generated using mapchart  <dig>   <cit> .
RESULTS
phenotyping
the experiments revealed variability for genotype and treatment in all analysed traits as shown in table  <dig>  for the traits biomass yield , leaf colour  and the electron transport rate  the mean values for the stress treatment were lower than in the well watered variant.
an exception are some genotypes of the sbcc  showing no decrease in spad and etr or even an increase, represented by negative values across treatments .
in contrast to the above mentioned traits, osmolality , content of free proline  and total content of soluble sugars  increased under drought stress.table  <dig> descriptive statistics, heritability  and number of significant  quantitative trait loci
traita
aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars
bcontrol and drought stress treatment, as well as dsi: drought susceptibility index across treatments by fischer & maurer
cminimum, maximum, mean, standard deviation in % , coefficient of variation   and least significant difference
the coefficient of variation  was comparable for control and drought stress treatment  for all six traits.
heritabilities  estimated ranged between  <dig> for oa to  <dig>  for by in the control treatment and  <dig>  for oa and  <dig>  for spad in the stress treatment.
generally, h <dig> was higher for the stress treatment except for by and spad.
analysis of variance  revealed significant  genotype and treatment effects for all investigated traits and genotype x treatment interactions for by, cfp and css .table  <dig> analysis of variance  of analysed traits showing f and p values
traita
aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars
bgenotype, treatment and gxt: genotype x treatment interaction effect
to get information on the influence of the physiological parameters estimated on biomass yield as the indicator for drought stress and spad as the indicator for drought stress induced leaf senescence, correlations to these traits were calculated .
for control and stress treatment by is significantly correlated to spad with r =  <dig>  and r =  <dig> , respectively.
a significant correlation was also determined for css to spad with r =  <dig>  and for cfp with r =  <dig>  in the drought stress treatment whereas for the control treatment significantly negative correlations were found.
low but nevertheless significant correlations to spad were also detected for etr and oa in the control treatment.
similar correlations were detected for by.
high and significant correlations were found between by and css  and cfp  for the drought stress treatment.
under control conditions the sbcc being a sub-population of its own, influences the correlation by producing less by which results in reduced shading of the primary leaves and a negative correlation especially to cfp and css.
by correlating only the german cultivars, these effects are excluded and no correlations  were observed.
etr and oa were not significantly correlated to by.table  <dig> coefficient of correlation  for control and drought stress treatment
aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars
r is significant with *p < <dig> , **p < <dig>  and ***p < <dig>
genotyping
the set of genotypes was analysed with the 9 k iselect snp-chip available for barley.
in summary  <dig>  snps turned out to be polymorphic.
out of these,  <dig>  snps are mapped on the seven barley chromosomes  <cit> , showing a minor allele frequency  >5 %.
this set of snps was used for the calculation of the linkage disequilibrium decay , which turned out be on average  <dig>  cm for this set of genotypes.
the number of subpopulations was estimated at k =  <dig> .fig.
<dig> optimal k of the population structure.
the number of subpopulations within the set of barley genotypes was estimated at k =  <dig> by calculation described in evanno et al.
<cit>
genome wide association study
results of gwas are shown in detail in additional file  <dig> and summarized in tables  <dig> and  <dig>   <dig> snps significantly  associated to traits estimated in the control variant,  <dig> significantly associated snps in the stress treatment and  <dig> significantly associated snps across treatments , were detected using the mlm analysis in tassel.
significant associations were found on all barley chromosomes.
most of the significant marker trait associations were located on barley chromosome 2h and 5h.
a large number of snps on chromosome 5h around 45 cm turned out to be associated to spad and by in the stress treatment .
most significant associations for these traits were detected on chromosome 5h.
the highest number of marker trait associations was detected for by.
no associations were observed for cfp in control and stress treatment, whereas across treatments one significant association was found.
for etr two significant associations and for css one significant association in the stress treatment explaining  <dig>  %  and  <dig>  %  of the phenotypic variance respectively were detected.
the strongest association was observed on chromosome 1h for by across treatments with a  value of  <dig>  explaining  <dig>  % of the phenotypic variance.
for oa  <dig> significant associations were detected in the stress treatment located over all barley chromosomes, whereas in the control treatment no associations and across treatments only one association was found on chromosome 4h.table  <dig> significant markers traits associations detected under drought stress conditions at a significance of p < <dig>
traita
by
<dig>
spad
<dig>
etr
<dig>
oa
<dig>
css
<dig>
total qtl
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, css: total content of soluble sugars
bone genomic region up to  <dig>  cm ; the chromosomal position in cm was taken from the respective snp with the highest r2
cchromosome positions are based on comadran et al.
afunction of the proteins related to drought stress  or leaf senescence , given by uniprot
bmarkers and chromosome positions are based on comadran et al.
cby: biomass yield, css: total content of soluble sugars, etr: electron transport rate at psii, cfp: content of free proline, oa: osmolality, spad: leaf colour
dc: control treatment, s: stress treatment, dsi: drought susceptibility index across treatments by fischer & maurer
fig.
<dig> manhattan plots.
showing –log p of association between illumina snps and the analysed traits biomass yield  and leaf colour  against the position of the markers on all seven linkage groups  for stress treatment.
threshold indicates a significance level of p < <dig>
for spad as an indicator of leaf senescence, twelve significant  marker trait associations in the control treatment, six under stress conditions and two across treatments were detected.
for the drought stress indicator biomass yield  <dig> significant marker trait associations were found in the control treatment,  <dig> in the stress treatment and  <dig> across treatments.
significant marker trait associations for by were evenly distributed over all chromosomes.
out of these, eight were identical in the stress treatment, control treatment and across treatments and a high number of  <dig> marker trait associations were identical in control treatment and across treatments.
in total the  <dig> significant associations estimated in the control treatment account for  <dig> genomic regions , the  <dig> significant associations determined in the stress treatment represent  <dig> genomic regions and the  <dig> significant associations across treatments  account for  <dig> genomic regions.
by comparing the localisation of marker trait associations detected, it turned out that significant associations were found for different traits at same positions, e.g.
on chromosome 2h at 50 cm for spad and oa in the stress treatment and at 120 cm for by and spad in the control treatment.
furthermore, on chromosome 5h at 45 cm significant marker trait associations for by, spad and oa were detected in the stress treatment, and at 95 cm significant associations for by and spad in the control treatment.
summarizing, overlapping of qtl was found across treatments and for different traits, especially for by and spad, which are also significantly correlated .
one interesting qtl was observed on chromosome 5h at 45 cm where a significant association to by and spad in the drought stress treatment was found, which is also within the ld of a significant association to oa .
therefore, at this position a putative major qtl for drought stress and leaf senescence may be located.
from the  <dig> significantly associated snps detected in the control treatment, stress treatment and across treatments ,  <dig> proteins were identified by an ncbi blast of the marker sequences.
out of these  <dig> proteins turned out to be related to drought stress,  <dig> proteins related to leaf senescence and  <dig> proteins turned out to be not related to drought stress or leaf senescence.
out of the  <dig> proteins for drought stress and leaf senescence ,  <dig> revealed associations under drought stress conditions.
these were in a next step assigned to the barley chromosomes by the known genetic localization of respective snps .
most of these were located at barley chromosome 2h and 5h, none were mapped on chromosome 6h.fig.
<dig> genetic map.
shows qtl including the significant associated snp marker positions for significant blasted proteins  linked to drought stress or leaf senescence, related traits for drought stress treatment and percentage of phenotypic variance  of the snps for all linkage groups
discussion
using the experimental design described above a clear cut influence of drought stress on biomass yield and physiological parameters was observed .
this effect of reduced biomass under drought conditions in juvenile stages was also reported in honsdorf et al.
<cit>  with an even higher reduction due to drought in early developmental stages of barley by jamieson et al.
<cit> .
in experiments on terminal drought stress application in barley this effect was not so pronounced  <cit>  giving hint that barley is most susceptible to drought stress in early developmental stages.
furthermore, a reduction of the chlorophyll content under drought stress conditions in barley has been observed  <cit> , but there are also reports on adverse effects, e.g.
for rapeseed or potato  <cit> , which may be due to a reduced leaf growth under drought stress conditions resulting in a reduced cell expansion leading to a relatively higher chlorophyll density in the leaves.
in the present study measurement was done on primary leaves, which were fully expanded at the initiation of drought stress, so that this effect was excluded.
for biomass yield a significant correlation to the leaf senescence parameter chlorophyll content was observed .
this correlation may be based on a true genetic relationship between these parameters, as it is also reported in drought stress field studies on wheat  <cit>  and in glasshouse experiments  <cit> .
the electron transport rate at ps ii as a parameter for the chlorophyll fluorescence decreased under drought stress indicating the degradation of chlorophyll during drought stress induced leaf senescence , as already shown by fang et al.
<cit> , li et al.
<cit> , netto et al.
<cit>  and silva et al.
<cit> .
as expected, osmolality increased under drought stress thereby protecting cells against a turgor loss  <cit> .
this is also reported in other drought stress studies on barley  <cit> , but is more often detected in barley under salt stress conditions  <cit> .
for oa no correlation was found to by, as also shown in a drought stress study on spring barley  <cit> .
the amino acid proline has been described as an osmo-protectant  <cit>  and is accumulated along with several abiotic stresses, such as drought stress, as seen in the present study .
this effect was also found in other drought stress studies, for example on barley in pot experiments  <cit> , in greenhouse  <cit>  and on arabidopsis thaliana in climate chamber experiments  <cit> .
the role of proline accumulation is still controversially discussed as it is described to function as a radical scavenger, antioxidant and is involved in the regulation of apoptosis and in seed development  <cit> .
high correlations were detected for cfp to spad and by giving hint that this trait is involved in drought stress tolerance and leaf senescence .
the correlation of the proline content to spad was also found in winter survival studies of barley  <cit>  and in studies on salt stress tolerance in trigonella foenum-graecum  <cit> .
up to now no correlations of cfp to by under drought stress have been described for barley, but positive correlations to yield in wheat were observed under drought stress  <cit> .
soluble sugars are acting also as osmo-protectants and consequently like in our study an increase was detected in several drought stress experiments on barley  <cit> , wheat  <cit> , potato  <cit>  and also pea  <cit> .
furthermore, studies showed that an increase of soluble sugars occurs along with leaf senescence  <cit>  and that css was correlated to leaf senescence and biomass production  <cit> .
interactions between sugar and aba signalling may be responsible for the induction of senescence during drought stress  <cit> .
quite high values for the heritability of respective traits estimated under drought stress conditions  give hint that such an experimental design is suited together with a set of diverse genotypes and the respective number of snp-markers to detect qtls using a genome wide association approach .
like in other studies  <cit> , the highest number of associations was detected for the traits with the highest heritability.
in the present study these were spad and by.
most associations were found on barley chromosomes 2h and 5h on which qtls were located at 50 cm and at 45 cm, respectively.
also in other gwas studies of barley significant qtls for spad and by were located on these chromosomes.
close to the qtl for spad located on chromosome 2h  a qtl for spad under drought stress was also mapped by li et al.
<cit> .
moreover, on chromosome 2h at 115 cm a qtl for spad was identified in a pot experiment with post-flowering drought stress  <cit> .
qtls for spad were also located on chromosome 2h at  <dig>  cm and on chromosome 5h at  <dig>  cm in mediterranean dry land experiments , but no significant marker trait associations were detected at these positions in our experiments.
the same holds true for a qtl for spad on chromosome 5h at 139 cm  <cit> .
varshney et al.
<cit>  also detected a qtl for biomass yield on chromosome 5h at 95 cm and 156 cm, which is near to associations, which were found in our study on chromosome 5h at 152 cm and 167 cm.
in addition, in the present study a qtl for spad and by under drought stress treatment was detected on chromosome 5h at 45 cm, which has not been described before.
furthermore, a lot of significant marker trait associations were observed for osmolality under stress treatment distributed over all barley chromosomes.
this was also reported for barley based on growth chamber drought experiments  <cit> .
proteins involved in drought stress and leaf senescence were detected by a blast of snp marker sequences and it turned out that they are distributed over all barley chromosomes, except 6h with a focus on chromosomes 2h and 5h .
most interesting proteins detected in the drought stress treatment  are discussed in detail.
on chromosome 1h an atp-dependent zinc metalloprotease  which is a regulator of heat shock proteins turned out to be associated to by under drought stress.
this protein is involved in the thylakoid formation and in the removal of damaged d <dig> in the photosystem ii, preventing cell death under high-intensity light conditions.
in interaction with heat shock proteins it reduces chlorophyll a/b ratios in heat tolerance regulation in arabidopsis thaliana  <cit> .
heat stress often occurs simultaneously with drought and also leads to leaf senescence.
besides this, on chromosome 1h the pectin methylesterase  <dig>  was found by sequence alignment to be associated to css.
this protein acts in the modification of cell walls via demethylesterification of pectin and turned out to be up-regulated by drought stress in rice  <cit> .
it influences the mechanical stability of cell walls and thereby also of leaves.
on chromosome 2h a sucrose synthase  <dig>  turned out to be associated with spad.
this is a sucrose-cleaving enzyme that provides udp-glucose and fructose for various metabolic pathways and is involved in nucleic acid break down during leaf senescence as revealed by expression analysis e.g.
in cucumber and rapeseed  <cit> .
the snp marker with the homolog sequences to this protein was associated to the leaf senescence parameter spad at  <dig>  cm.
so there may be a direct relationship between the spad values and the activation of this enzyme, especially because this snp marker explains  <dig>  % of the phenotypic variance.
furthermore, a probable glutamate carboxypeptidase  revealed an association to oa which plays an important role in shoot apical meristem development and phytohormone homeostasis.
by microarray analysis it turned out that amp mediates aba production and is involved in abiotic stress response such as drought stress in arabidopsis thaliana  <cit> .
moreover, a marker with a sequence homologue to a phospholipid hydroperoxide glutathione peroxidase  associated to by was found which protects cells and enzymes from oxidative damage.
photometrical analyses of protein quantity and activity showed that the expression of gpx <dig> and gpx <dig> is reduced under drought stress and restored after recovery in winter wheat  <cit> .
on chromosome 3h an f-box protein was detected  associated to oa.
these proteins are ubiquitin related and negatively regulate aba mediated drought stress response in arabidopsis thaliana  <cit> .
furthermore, an association of oa to 1-aminocyclopropane-1-carboxylate oxidase  was detected which limits leaf growth by inhibiting the ethylene biosynthesis and so leads to drought tolerance.
this was figured out in barley by expression analyses of protein related genes  <cit> .
surprisingly, in the control treatment associated to by, a well known drought stress related protein, i.e.
dehydrin   <cit>  was found which was also in another study located on chromosome 3h  <cit> .
this protein belongs to the family of late embryogenesis abundant  proteins and is reported to be up-regulated in the protection mechanisms activated by plants in response to drought stress in wheat  <cit> .
on chromosome 4h an aba receptor  was located and associated to oa that activates aba signalling and aba-mediated responses such as stomatal closure and germination inhibition.
immuno-detection experiments of protein extracts revealed that aba signalling is involved in several stresses for example drought stress in arabidopsis thaliana  <cit> .
on chromosome 5h nucleotide pyrophosphatase/phosphodiesterase  was found associated to by which facilitates auxin transport by modulating apoplastic ph and regulating auxin-mediated developmental processes.
increased expression of protein related genes in transgenic barley confers tolerance to nacl and to drought by increasing ion retention  <cit> .
furthermore, three protein kinases associated to spad were located on this chromosome regulating protein activity by phosphorylation.
first a serine/threonine-protein kinase  which leads to stress induced programmed cell death, shown in arabidopsis thaliana by expression profiles of protein related genes  <cit> , second a serine/threonine-protein kinase  which is activated by hyperosmotic stress in rice  <cit>  and third the abscisic acid-inducible protein kinase  observed in wheat, which is also involved in dehydration stress response  <cit> .
moreover, an anthocyanin regulatory protein  was detected associated to by.
anthocyanin is often accumulated in abiotic stress response, among others in drought stress with a photoprotective function as shown in arabidopsis thaliana  <cit> .
all of these proteins showed homologies to sequences of snps at chromosome 5h around 45 cm and were associated significantly to by or spad, representing an interesting candidate qtl for drought stress and leaf senescence.
in addition, a dehydration-responsive element-binding protein  was found associated to oa, which delays water stress symptoms and promotes expression of drought tolerance genes in transgenic wheat  <cit> .
on chromosome 7h another ethylene responsive protein was found to be associated to by under drought stress.
the transcription factor erf <dig> is involved in the regulation of gene expression by stress factors  and progression of leaf senescence in arabidopsis thaliana  <cit> .
besides, on chromosome 7h the cation/h antiporter  <dig>  was detected associated to etr, which is important for ph gradients in the cell.
this protein plays a vital role in maintaining both cellular and intercellular ionic balances under stresses such as drought stress as observed in arabidopsis thaliana  <cit> .
in summary the blast of the associated snps to protein data bases revealed many proteins which are known to be involved in drought stress response or leaf senescence, respectively giving hint that the gwas approach is well suited for the genetic dissection of these traits in barley.
out of the qtl detected, the ones on chromosome 2h at 50 cm and chromosome 5h at 45 cm are of prime importance and may be involved in breeding barley for drought tolerance in the future due to the quite high amount of phenotypic variance explained.
CONCLUSIONS
by gwas marker trait associations for above ground biomass and physiological traits involved in drought stress tolerance and leaf senescence in early developmental stages of barley were detected.
major qtl for by and spad under drought stress were located at chromosome 2h at 50 cm and chromosome 5h at 45 cm, giving hint that in these regions putative major qtls for drought stress and leaf senescence are located.
with respect to the qtl on chromosome 2h, qtls for drought stress and leaf senescence were located at comparable positions in other gwas studies while the one on chromosome 5h was detected for the first time.
by blastx of respective snp carrying sequences,  <dig> proteins were identified being involved in drought stress or leaf senescence, respectively.
respective qtls may be the starting point for marker based selection in barley for drought stress tolerance in the juvenile stage.
additional files
additional file 1:
overview of the  <dig> analysed genotypes.
asbcc: spanish barley core collection.
additional file 2:
significant associations  of the genome wide association study and blasted proteins out of the marker sequences.
aby: biomass yield, css: total content of soluble sugars, etr: electron transport rate at psii, cfp: content of free proline, oa: osmolality, spad: leaf colour.
bc: control treatment, s: stress treatment, dsi: drought susceptibility index across treatments by fischer & maurer .
cmarkers and chromosome positions are based on comadran et al.
.
abbreviations
e.g.for example
i.e.id est
gwasgenome wide association study
qtlquantitative trait locus
snpsingle nucleotide polymorphism
abaabscisic acid
lealate embryogenesis abundant protein
spadsoil plant analysis development; measurement of chlorophyll content by colour
sbccspanish barley core collection
dasdays after sowing
etrelectron transport rate
cfpcontent of free proline
csstotal content of soluble sugars
oaosmolality
bybiomass yield
psiiphotosystem two
parphotosynthetic active radiation
dsidrought susceptibility index; value across treatments
lsmeanslast square means
h2heritability
cvcoefficient of variation
ssrsimple sequence repeat
mafminor allele frequency
knumber of subpopulations
mlmmixed linear model
ldlinkage disequilibrium
blastbasic local alignment search tool
competing interests
the authors declare that they have no competing interests.
authors' contributions
gw conducted all experiments, including statistical and bioinformatics analyses, as well as the sequence alignment and mainly wrote the manuscript.
me participated in the gwas and ld calculation.
cb, kh and fo designed the research, supervised the experimental design, data analysis and participated in writing the manuscript.
all authors approved the final manuscript.

 barley is one of the first cereals domesticated in the fertile crescent and today it is the fourth most important crop species concerning acreage next to wheat, maize, and rice. worldwide, barley is mainly used for animal feed and malting and only a very small amount is used for direct human consumption and bakery. average yield of barley on the worldwide level is  t/ha but in some european countries, e.g. germany average yield is up to  t/ha. barley yield in many parts of the world is reduced by biotic stress but also by abiotic stress e.g. heat, salt, deficits in nitrogen nutrition and drought. especially, in the juvenile stages from sowing to tillering, drought can severely influence barley development already reducing the potential yield. research on drought stress tolerance has become more important worldwide as due to climate change the number of drought periods will increase in the future. up to now, most studies conducted in barley focused on effects of terminal drought stress whereas drought in juvenile stages is less well documented. drought tolerance is a complex quantitative trait, that is controlled by various mechanisms. abscisic acid is a key phytohormone involved in adaption to environmental stresses and regulation of plant development. it promotes the closure of stomata under drought stress conditions initiated by a loss of turgor. furthermore, it increases the hydraulic conductivity of water, promotes chlorophyll breakdown and leads to leaf senescence. another relevant protein is ubiquitin which regulates the degradation of proteins. moreover, late embryo abundant proteins and heat shock proteins that are involved in the protection of functional proteins are induced in response to various abiotic stresses. for example protein kinases and protein phosphatases which activate or deactivate proteins by phosphorylation and dephosphorilation, as well as the lea protein dehydrin which is described to have different functions in different stresses is often present along with drought stress. a lot of parameters indicative for drought stress influenced by these and additional genes were analysed in different crops. for example, biomass production, yield, photosynthesis rate, as well as the content of free proline, total content of soluble sugars, or osmolality are parameters which are affected by drought stress in barley. another factor relevant for yield improvement is leaf senescence, which is a natural degradation process at the final stage of the development of organs and plants. this process is divided into three steps and starts with reprogramming of gene expression to turn on senescence activating genes. before programmed cell death in the terminal phase occurs, the second step in which nutrients and metabolites are transported from source to sink is important for yield and quality of the seeds harvested. during this reorganization phase, degradation of chlorophyll and a decrease in photosynthesis is observed. because of degradation of chlorophyll, yellowing of the leaves is a symptom of leaf senescence, which in many studies is rated visually, but can be more precisely determined by soil plant analysis development readings which estimate leaf greenness. degradation of chlorophyll is regulated by chlorophyllase, pheophorbide α oxygenase and red chlorophyll catabolite reductase among others, but so far for regulation of leaf senescence only a few genes are known. leaf senescence is a process which is influenced by a lot of external stress conditions e.g. drought stress. stress often results in premature induction of leaf senescence and therefore leads to an inefficient recycling of resources and a massive yield loss. in contrast, plants showing delayed leaf senescence under stress, represented by a “stay green effect”, minimize yield loss. genome wide association studies are a powerful tool to subdivide such complex pathways as drought stress and leaf senescence by the detection of quantitative trait loci out of the regression analysis of genotypic and phenotypic data. up to now, some qtls involved in drought stress response were published in barley, whereas for leaf senescence only few qtls are known. besides this, qtls which are involved in the response to other abiotic stresses, e.g. salt stress were identified in barley. molecular markers, such as single nucleotide polymorphisms flanking qtls having a significant influence on the respective trait can be used for efficient marker assisted selection and smart breeding procedures. the aim of the present study is therefore the identification of qtls for drought stress and drought stress induced leaf senescence in early developmental stages of barley suited to be used in future barley breeding programs using gwas followed by the identification of the function of these qtls associated to respective traits. material & methods plant material and experimental setup a set of winter barley genotypes consisting of german winter barley cultivars [ two-rowed and six-rowed, and accessions of the spanish barley core collection were used to investigate drought stress induced leaf senescence in juvenile barley plants. drought stress was applied in greenhouses of the julius kühn-institut in groß lüsewitz, germany according to honsdorf et al.. trials were conducted in a split plot design with three replications per genotype and variant. ten seeds of each accession were sown per plastic pot containing  g of a mixed clay soil ed. after germination, seedlings were reduced to seven plants per pot. plants were grown under semi-controlled long day conditions in a temperature range from to 22 °c at day and to 19 °c at night. if natural radiation was below klx, additional light was applied from 6 a.m. to 10 p.m. drought stress started at the primary leaf stage seven days after sowing. at this time watering of the stress variant was stopped till the soil reached 20 % of the maximal soil water capacity, and then this level was kept by weighting each pot and re-watering. control plants were continuously watered to 70 % of the maximal soil water capacity. water capacity was calculated of the saturated soil weight and drought weight according to din iso 1996–. at the end of a four weeks stress period physiological traits were determined and above ground biomass was harvested. experimental setup a was repeated in three years. a modified experimental setup was conducted to optimize lightening conditions of the primary leaves for the measurement of leaf colour and electron transport rate at photosystem ii. general settings were the same as in setup a, but only four plants were grown in smaller pots to allow wider spacing of pots and all leaves except the primary leaves were tied up to reduce shading. experimental setup b was repeated in two years. physiological parameters determined six physiological traits, i.e. leaf colour, electron transport rate at photosystem ii, content of free proline, content of soluble sugars, osmolality and the aboveground biomass yield were determined in the control and stress treatment. measurement and sampling respectively were conducted on primary leaves. chlorophyll content which was used as the main indicator for drought stress induced leaf senescence was measured 33– das by minolta spad readings, which gives a value for leaf colour. three primary leaves of three plants for each pot were measured at five positions per leaf. these spad readings turned out to be correlated to the chlorophyll content analysed photometrically. because of this relation, the chlorophyll content can be indirectly measured by spad. at 34– das chlorophyll fluorescence was measured in all genotypes using light adapted plants with the os1p-chlorophyll fluorometer in the middle of three primary leaves per pot at one position per leaf. the relative electron transport rate at photosystem ii  * par *   84 *   5) was calculated including the photosynthetically active radiation, as well as the quantum photosynthetic yield of psii ) and constants representing light which is absorbed by the leaf and light which is equally absorbed by psi and psii. at das five primary leaves per pot were harvested and cut in pieces of 1 cm length for the analysis of cfp, css and oa. these samples were frozen in liquid nitrogen immediately and samples for cfp and css measurement were freeze dried. for cfp measurement the ninhydrin method was applied, and for css measurement the anthron method was used. both traits were measured photometrical using a spectrophotometer. the concentration of these ingredients was determined with a standard curve calculated on a dry weight basis. to assess oa, frozen leaf samples were grinded in a swing mill, filled up with 200 μl water and centrifuged at  rpm for 15 min to get cell sap for the measurement of osmolality with a freezing-point osmometer. osmolality was corrected for the water content of fresh and dry weight. above ground biomass was harvested das, too. leaf material was dried in a compartment dryer at 105 °c and weighted. for all traits an outlier test was calculated to exclude extreme deviations. to get information on the stability of all analysed traits in response to drought stress compared to the control, the drought susceptibility index was calculated for each trait and across treatments according to the formula: dsi=1–lsmeantraitgenotype,stress/lsmeantraitgenotype,control/d,with d=1–lsmeantraitassortment,stress/lsmeantraitassortment,control dsi is a relative value estimated for each genotype and trait. according to the formula genotypes revealing a dsi close to one are highly susceptible to drought and those close to zero or showing a negative value are tolerant. statistical analysis of phenotypic data statistical analyses were performed with sas. least square means were calculated with glm procedure for the replications of each genotype in the respective years and for both control and drought conditions. descriptive statistics was calculated out of lsmeans by proc univariate. analysis of variance was calculated using proc mixed with genotypes, drought stress treatment and gxt as fixed effects. replication, year and row type are chosen as random factors. the heritability was calculated with sas in two steps. first, the variance components for the genotypes, variance associated with the genotype by year interaction and ve which is the error variance were calculated with proc varcomp. next, h was calculated with the following formula: vg/ for both well watered and drought stress conditions. furthermore, the coefficient of correlation by pearson was calculated with sas based on lsmeans. genotyping and genome wide association study for genotyping the whole set of genotypes was analysed with the barley illumina 9 k iselect snp-chip. population structure was calculated with structure. based on simple sequence repeat markers covering the whole genome. the structure programme was run times for pre-defined k from to each. to get the number of calculated subpopulations with highest likelihood the procedure of evanno et al. was applied. an independent run with iterations of a monte carlo markov chain with a length of the burn in period of was conducted for the k with the highest likelihood to obtain the q-matrix. kinship was calculated with spagedi d based on ssrs and allele size correlation coefficient with permutations. out of genetically mapped snp markers, polymorphic markers with minor allele frequencies higher than 5 % were taken into account. based on these data and respective phenotypic data gwas was conducted applying a mixed linear model using tassel. all results with p values < were considered as significant marker trait associations. linkage disequilibrium was calculated on mapped polymorphic snps with r by an estimate of the average decay over all barley chromosomes. sequences of significantly associated snp-markers were downloaded from the james hutton institute and respective sequences were compared against the plant proteome in the uniprotkb/swiss-prot protein database by blastx to get information on the proteins coded by these sequences. using uniprot the involvement of respective proteins in drought stress and leaf senescence processes was analysed. in a last step a genetic map with all significantly associated snps in genes coding for proteins known to be involved in drought stress tolerance was generated using mapchart.  phenotyping the experiments revealed variability for genotype and treatment in all analysed traits as shown in table  for the traits biomass yield, leaf colour and the electron transport rate the mean values for the stress treatment were lower than in the well watered variant. an exception are some genotypes of the sbcc showing no decrease in spad and etr or even an increase, represented by negative values across treatments. in contrast to the above mentioned traits, osmolality, content of free proline and total content of soluble sugars increased under drought stress.table descriptive statistics, heritability and number of significant quantitative trait loci traita aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars bcontrol and drought stress treatment, as well as dsi: drought susceptibility index across treatments by fischer & maurer cminimum, maximum, mean, standard deviation in %, coefficient of variation and least significant difference the coefficient of variation was comparable for control and drought stress treatment for all six traits. heritabilities estimated ranged between for oa to for by in the control treatment and for oa and for spad in the stress treatment. generally, h was higher for the stress treatment except for by and spad. analysis of variance revealed significant genotype and treatment effects for all investigated traits and genotype x treatment interactions for by, cfp and css.table analysis of variance of analysed traits showing f and p values traita aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars bgenotype, treatment and gxt: genotype x treatment interaction effect to get information on the influence of the physiological parameters estimated on biomass yield as the indicator for drought stress and spad as the indicator for drought stress induced leaf senescence, correlations to these traits were calculated. for control and stress treatment by is significantly correlated to spad with r =  and r = , respectively. a significant correlation was also determined for css to spad with r =  and for cfp with r =  in the drought stress treatment whereas for the control treatment significantly negative correlations were found. low but nevertheless significant correlations to spad were also detected for etr and oa in the control treatment. similar correlations were detected for by. high and significant correlations were found between by and css and cfp for the drought stress treatment. under control conditions the sbcc being a sub-population of its own, influences the correlation by producing less by which results in reduced shading of the primary leaves and a negative correlation especially to cfp and css. by correlating only the german cultivars, these effects are excluded and no correlations were observed. etr and oa were not significantly correlated to by.table coefficient of correlation for control and drought stress treatment aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, cfp: content of free proline, css: total content of soluble sugars r is significant with *p <, **p < and ***p < genotyping the set of genotypes was analysed with the 9 k iselect snp-chip available for barley. in summary snps turned out to be polymorphic. out of these, snps are mapped on the seven barley chromosomes, showing a minor allele frequency >5 %. this set of snps was used for the calculation of the linkage disequilibrium decay, which turned out be on average  cm for this set of genotypes. the number of subpopulations was estimated at k = .fig. optimal k of the population structure. the number of subpopulations within the set of barley genotypes was estimated at k =  by calculation described in evanno et al. genome wide association study results of gwas are shown in detail in additional file and summarized in tables  and snps significantly associated to traits estimated in the control variant, significantly associated snps in the stress treatment and significantly associated snps across treatments, were detected using the mlm analysis in tassel. significant associations were found on all barley chromosomes. most of the significant marker trait associations were located on barley chromosome 2h and 5h. a large number of snps on chromosome 5h around 45 cm turned out to be associated to spad and by in the stress treatment. most significant associations for these traits were detected on chromosome 5h. the highest number of marker trait associations was detected for by. no associations were observed for cfp in control and stress treatment, whereas across treatments one significant association was found. for etr two significant associations and for css one significant association in the stress treatment explaining  % and  % of the phenotypic variance respectively were detected. the strongest association was observed on chromosome 1h for by across treatments with a value of explaining  % of the phenotypic variance. for oa significant associations were detected in the stress treatment located over all barley chromosomes, whereas in the control treatment no associations and across treatments only one association was found on chromosome 4h.table significant markers traits associations detected under drought stress conditions at a significance of p < traita by spad etr oa css total qtl aby: biomass yield, spad: leaf colour, etr: electron transport rate at psii, oa: osmolality, css: total content of soluble sugars bone genomic region up to  cm; the chromosomal position in cm was taken from the respective snp with the highest r2 cchromosome positions are based on comadran et al. afunction of the proteins related to drought stress or leaf senescence, given by uniprot bmarkers and chromosome positions are based on comadran et al. cby: biomass yield, css: total content of soluble sugars, etr: electron transport rate at psii, cfp: content of free proline, oa: osmolality, spad: leaf colour dc: control treatment, s: stress treatment, dsi: drought susceptibility index across treatments by fischer & maurer fig. manhattan plots. showing –log p of association between illumina snps and the analysed traits biomass yield and leaf colour against the position of the markers on all seven linkage groups for stress treatment. threshold indicates a significance level of p < for spad as an indicator of leaf senescence, twelve significant marker trait associations in the control treatment, six under stress conditions and two across treatments were detected. for the drought stress indicator biomass yield significant marker trait associations were found in the control treatment, in the stress treatment and across treatments. significant marker trait associations for by were evenly distributed over all chromosomes. out of these, eight were identical in the stress treatment, control treatment and across treatments and a high number of marker trait associations were identical in control treatment and across treatments. in total the significant associations estimated in the control treatment account for genomic regions, the significant associations determined in the stress treatment represent genomic regions and the significant associations across treatments account for genomic regions. by comparing the localisation of marker trait associations detected, it turned out that significant associations were found for different traits at same positions, e.g. on chromosome 2h at 50 cm for spad and oa in the stress treatment and at 120 cm for by and spad in the control treatment. furthermore, on chromosome 5h at 45 cm significant marker trait associations for by, spad and oa were detected in the stress treatment, and at 95 cm significant associations for by and spad in the control treatment. summarizing, overlapping of qtl was found across treatments and for different traits, especially for by and spad, which are also significantly correlated. one interesting qtl was observed on chromosome 5h at 45 cm where a significant association to by and spad in the drought stress treatment was found, which is also within the ld of a significant association to oa. therefore, at this position a putative major qtl for drought stress and leaf senescence may be located. from the significantly associated snps detected in the control treatment, stress treatment and across treatments, proteins were identified by an ncbi blast of the marker sequences. out of these proteins turned out to be related to drought stress, proteins related to leaf senescence and proteins turned out to be not related to drought stress or leaf senescence. out of the proteins for drought stress and leaf senescence, revealed associations under drought stress conditions. these were in a next step assigned to the barley chromosomes by the known genetic localization of respective snps. most of these were located at barley chromosome 2h and 5h, none were mapped on chromosome 6h.fig. genetic map. shows qtl including the significant associated snp marker positions for significant blasted proteins linked to drought stress or leaf senescence, related traits for drought stress treatment and percentage of phenotypic variance of the snps for all linkage groups discussion using the experimental design described above a clear cut influence of drought stress on biomass yield and physiological parameters was observed. this effect of reduced biomass under drought conditions in juvenile stages was also reported in honsdorf et al. with an even higher reduction due to drought in early developmental stages of barley by jamieson et al.. in experiments on terminal drought stress application in barley this effect was not so pronounced giving hint that barley is most susceptible to drought stress in early developmental stages. furthermore, a reduction of the chlorophyll content under drought stress conditions in barley has been observed, but there are also reports on adverse effects, e.g. for rapeseed or potato, which may be due to a reduced leaf growth under drought stress conditions resulting in a reduced cell expansion leading to a relatively higher chlorophyll density in the leaves. in the present study measurement was done on primary leaves, which were fully expanded at the initiation of drought stress, so that this effect was excluded. for biomass yield a significant correlation to the leaf senescence parameter chlorophyll content was observed. this correlation may be based on a true genetic relationship between these parameters, as it is also reported in drought stress field studies on wheat and in glasshouse experiments. the electron transport rate at ps ii as a parameter for the chlorophyll fluorescence decreased under drought stress indicating the degradation of chlorophyll during drought stress induced leaf senescence, as already shown by fang et al., li et al., netto et al. and silva et al.. as expected, osmolality increased under drought stress thereby protecting cells against a turgor loss. this is also reported in other drought stress studies on barley, but is more often detected in barley under salt stress conditions. for oa no correlation was found to by, as also shown in a drought stress study on spring barley. the amino acid proline has been described as an osmo-protectant and is accumulated along with several abiotic stresses, such as drought stress, as seen in the present study. this effect was also found in other drought stress studies, for example on barley in pot experiments, in greenhouse and on arabidopsis thaliana in climate chamber experiments. the role of proline accumulation is still controversially discussed as it is described to function as a radical scavenger, antioxidant and is involved in the regulation of apoptosis and in seed development. high correlations were detected for cfp to spad and by giving hint that this trait is involved in drought stress tolerance and leaf senescence. the correlation of the proline content to spad was also found in winter survival studies of barley and in studies on salt stress tolerance in trigonella foenum-graecum. up to now no correlations of cfp to by under drought stress have been described for barley, but positive correlations to yield in wheat were observed under drought stress. soluble sugars are acting also as osmo-protectants and consequently like in our study an increase was detected in several drought stress experiments on barley, wheat, potato and also pea. furthermore, studies showed that an increase of soluble sugars occurs along with leaf senescence and that css was correlated to leaf senescence and biomass production. interactions between sugar and aba signalling may be responsible for the induction of senescence during drought stress. quite high values for the heritability of respective traits estimated under drought stress conditions give hint that such an experimental design is suited together with a set of diverse genotypes and the respective number of snp-markers to detect qtls using a genome wide association approach. like in other studies, the highest number of associations was detected for the traits with the highest heritability. in the present study these were spad and by. most associations were found on barley chromosomes 2h and 5h on which qtls were located at 50 cm and at 45 cm, respectively. also in other gwas studies of barley significant qtls for spad and by were located on these chromosomes. close to the qtl for spad located on chromosome 2h a qtl for spad under drought stress was also mapped by li et al.. moreover, on chromosome 2h at 115 cm a qtl for spad was identified in a pot experiment with post-flowering drought stress. qtls for spad were also located on chromosome 2h at  cm and on chromosome 5h at  cm in mediterranean dry land experiments, but no significant marker trait associations were detected at these positions in our experiments. the same holds true for a qtl for spad on chromosome 5h at 139 cm. varshney et al. also detected a qtl for biomass yield on chromosome 5h at 95 cm and 156 cm, which is near to associations, which were found in our study on chromosome 5h at 152 cm and 167 cm. in addition, in the present study a qtl for spad and by under drought stress treatment was detected on chromosome 5h at 45 cm, which has not been described before. furthermore, a lot of significant marker trait associations were observed for osmolality under stress treatment distributed over all barley chromosomes. this was also reported for barley based on growth chamber drought experiments. proteins involved in drought stress and leaf senescence were detected by a blast of snp marker sequences and it turned out that they are distributed over all barley chromosomes, except 6h with a focus on chromosomes 2h and 5h. most interesting proteins detected in the drought stress treatment are discussed in detail. on chromosome 1h an atp-dependent zinc metalloprotease which is a regulator of heat shock proteins turned out to be associated to by under drought stress. this protein is involved in the thylakoid formation and in the removal of damaged d in the photosystem ii, preventing cell death under high-intensity light conditions. in interaction with heat shock proteins it reduces chlorophyll a/b ratios in heat tolerance regulation in arabidopsis thaliana. heat stress often occurs simultaneously with drought and also leads to leaf senescence. besides this, on chromosome 1h the pectin methylesterase was found by sequence alignment to be associated to css. this protein acts in the modification of cell walls via demethylesterification of pectin and turned out to be up-regulated by drought stress in rice. it influences the mechanical stability of cell walls and thereby also of leaves. on chromosome 2h a sucrose synthase turned out to be associated with spad. this is a sucrose-cleaving enzyme that provides udp-glucose and fructose for various metabolic pathways and is involved in nucleic acid break down during leaf senescence as revealed by expression analysis e.g. in cucumber and rapeseed. the snp marker with the homolog sequences to this protein was associated to the leaf senescence parameter spad at  cm. so there may be a direct relationship between the spad values and the activation of this enzyme, especially because this snp marker explains  % of the phenotypic variance. furthermore, a probable glutamate carboxypeptidase revealed an association to oa which plays an important role in shoot apical meristem development and phytohormone homeostasis. by microarray analysis it turned out that amp mediates aba production and is involved in abiotic stress response such as drought stress in arabidopsis thaliana. moreover, a marker with a sequence homologue to a phospholipid hydroperoxide glutathione peroxidase associated to by was found which protects cells and enzymes from oxidative damage. photometrical analyses of protein quantity and activity showed that the expression of gpx and gpx is reduced under drought stress and restored after recovery in winter wheat. on chromosome 3h an f-box protein was detected associated to oa. these proteins are ubiquitin related and negatively regulate aba mediated drought stress response in arabidopsis thaliana. furthermore, an association of oa to 1-aminocyclopropane-1-carboxylate oxidase was detected which limits leaf growth by inhibiting the ethylene biosynthesis and so leads to drought tolerance. this was figured out in barley by expression analyses of protein related genes. surprisingly, in the control treatment associated to by, a well known drought stress related protein, i.e. dehydrin was found which was also in another study located on chromosome 3h. this protein belongs to the family of late embryogenesis abundant proteins and is reported to be up-regulated in the protection mechanisms activated by plants in response to drought stress in wheat. on chromosome 4h an aba receptor was located and associated to oa that activates aba signalling and aba-mediated responses such as stomatal closure and germination inhibition. immuno-detection experiments of protein extracts revealed that aba signalling is involved in several stresses for example drought stress in arabidopsis thaliana. on chromosome 5h nucleotide pyrophosphatase/phosphodiesterase was found associated to by which facilitates auxin transport by modulating apoplastic ph and regulating auxin-mediated developmental processes. increased expression of protein related genes in transgenic barley confers tolerance to nacl and to drought by increasing ion retention. furthermore, three protein kinases associated to spad were located on this chromosome regulating protein activity by phosphorylation. first a serine/threonine-protein kinase which leads to stress induced programmed cell death, shown in arabidopsis thaliana by expression profiles of protein related genes, second a serine/threonine-protein kinase which is activated by hyperosmotic stress in rice and third the abscisic acid-inducible protein kinase observed in wheat, which is also involved in dehydration stress response. moreover, an anthocyanin regulatory protein was detected associated to by. anthocyanin is often accumulated in abiotic stress response, among others in drought stress with a photoprotective function as shown in arabidopsis thaliana. all of these proteins showed homologies to sequences of snps at chromosome 5h around 45 cm and were associated significantly to by or spad, representing an interesting candidate qtl for drought stress and leaf senescence. in addition, a dehydration-responsive element-binding protein was found associated to oa, which delays water stress symptoms and promotes expression of drought tolerance genes in transgenic wheat. on chromosome 7h another ethylene responsive protein was found to be associated to by under drought stress. the transcription factor erf is involved in the regulation of gene expression by stress factors and progression of leaf senescence in arabidopsis thaliana. besides, on chromosome 7h the cation/h antiporter was detected associated to etr, which is important for ph gradients in the cell. this protein plays a vital role in maintaining both cellular and intercellular ionic balances under stresses such as drought stress as observed in arabidopsis thaliana. in summary the blast of the associated snps to protein data bases revealed many proteins which are known to be involved in drought stress response or leaf senescence, respectively giving hint that the gwas approach is well suited for the genetic dissection of these traits in barley. out of the qtl detected, the ones on chromosome 2h at 50 cm and chromosome 5h at 45 cm are of prime importance and may be involved in breeding barley for drought tolerance in the future due to the quite high amount of phenotypic variance explained.  by gwas marker trait associations for above ground biomass and physiological traits involved in drought stress tolerance and leaf senescence in early developmental stages of barley were detected. major qtl for by and spad under drought stress were located at chromosome 2h at 50 cm and chromosome 5h at 45 cm, giving hint that in these regions putative major qtls for drought stress and leaf senescence are located. with respect to the qtl on chromosome 2h, qtls for drought stress and leaf senescence were located at comparable positions in other gwas studies while the one on chromosome 5h was detected for the first time. by blastx of respective snp carrying sequences, proteins were identified being involved in drought stress or leaf senescence, respectively. respective qtls may be the starting point for marker based selection in barley for drought stress tolerance in the juvenile stage. additional files additional file 1: overview of the analysed genotypes. asbcc: spanish barley core collection. additional file 2: significant associations of the genome wide association study and blasted proteins out of the marker sequences. aby: biomass yield, css: total content of soluble sugars, etr: electron transport rate at psii, cfp: content of free proline, oa: osmolality, spad: leaf colour. bc: control treatment, s: stress treatment, dsi: drought susceptibility index across treatments by fischer & maurer. cmarkers and chromosome positions are based on comadran et al.. abbreviations e.g.for example i.e.id est gwasgenome wide association study qtlquantitative trait locus snpsingle nucleotide polymorphism abaabscisic acid lealate embryogenesis abundant protein spadsoil plant analysis development; measurement of chlorophyll content by colour sbccspanish barley core collection dasdays after sowing etrelectron transport rate cfpcontent of free proline csstotal content of soluble sugars oaosmolality bybiomass yield psiiphotosystem two parphotosynthetic active radiation dsidrought susceptibility index; value across treatments lsmeanslast square means h2heritability cvcoefficient of variation ssrsimple sequence repeat mafminor allele frequency knumber of subpopulations mlmmixed linear model ldlinkage disequilibrium blastbasic local alignment search tool competing interests the authors declare that they have no competing interests. authors' contributions gw conducted all experiments, including statistical and bioinformatics analyses, as well as the sequence alignment and mainly wrote the manuscript. me participated in the gwas and ld calculation. cb, kh and fo designed the research, supervised the experimental design, data analysis and participated in writing the manuscript. all authors approved the final manuscript. 
major qtls for drought stress and leaf senescence were located on chromosome 5h and 2h. gwas resulted in the identification of many qtls involved in drought stress and leaf senescence of which two major qtls for drought stress and leaf senescence were located on chromosome 5h and 2h. results may be the basis to incorporate breeding for tolerance to drought stress or leaf senescence in barley breeding via marker based selection procedures. drought stress, is a main factor for yield losses in barley. research in drought stress tolerance has become more important as due to climate change the number of drought periods will increase and tolerance to drought stress has become a goal of high interest in barley breeding. therefore, the aim is to identify quantitative trait loci involved in drought stress induced leaf senescence and drought stress tolerance in early developmental stages of barley by applying genome wide association studies on a set of winter barley genotypes. after a four weeks stress period leaf colour as an indicator of leaf senescence, electron transport rate at photosystem ii, content of free proline, content of soluble sugars, osmolality and the aboveground biomass indicative for drought stress response were determined in the control and stress variant in greenhouse pot experiments. blastx search for associated marker sequences revealed that respective snps are in some cases located in proteins related to drought stress or leaf senescence, e.g. barleyleaf senescencedrought stressgwasqtlissue-copyright-statement© the author 2015 based on these phenotypic data and polymorphic single nucleotide polymorphisms with a minor allele frequency >5 % derived from the illumina 9 k iselect snp chip, qtl were detected for all traits analysed. heritabilities ranged between for osmolality and for leaf colour in stress treatment and significant effects of genotype, treatment and genotype x treatment were estimated for most traits analysed. premature leaf senescence induced by external stress conditions, e.g.
 premature leaf senescence induced by external stress conditions, e.g. drought stress, is a main factor for yield losses in barley. research in drought stress tolerance has become more important as due to climate change the number of drought periods will increase and tolerance to drought stress has become a goal of high interest in barley breeding. therefore, the aim is to identify quantitative trait loci involved in drought stress induced leaf senescence and drought stress tolerance in early developmental stages of barley by applying genome wide association studies on a set of winter barley genotypes.  after a four weeks stress period leaf colour as an indicator of leaf senescence, electron transport rate at photosystem ii, content of free proline, content of soluble sugars, osmolality and the aboveground biomass indicative for drought stress response were determined in the control and stress variant in greenhouse pot experiments. significant phenotypic variation was observed for all traits analysed. heritabilities ranged between for osmolality and for leaf colour in stress treatment and significant effects of genotype, treatment and genotype x treatment were estimated for most traits analysed. based on these phenotypic data and polymorphic single nucleotide polymorphisms with a minor allele frequency >5 % derived from the illumina 9 k iselect snp chip, qtl were detected for all traits analysed. major qtls for drought stress and leaf senescence were located on chromosome 5h and 2h. blastx search for associated marker sequences revealed that respective snps are in some cases located in proteins related to drought stress or leaf senescence, e.g. nucleotide pyrophosphatase or serine/ threonin protein kinase.  gwas resulted in the identification of many qtls involved in drought stress and leaf senescence of which two major qtls for drought stress and leaf senescence were located on chromosome 5h and 2h. results may be the basis to incorporate breeding for tolerance to drought stress or leaf senescence in barley breeding via marker based selection procedures. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords barleyleaf senescencedrought stressgwasqtlissue-copyright-statement© the author 2015 