Differential Water Deficit in Leaves Is a Principal Factor Modifying Barley Response to Drought Stress

In response to environmental stress, plants activate complex signalling, including being dependent on reactive oxygen–nitrogen–sulphur species. One of the key abiotic stresses is drought. As a result of drought, changes in the level of hydration of the plant occur, which obviously entails various metabolic alternations. The primary aim of this study was to determine the relationship between the response of barley to drought and the intensity of stress, therefore investigations were performed under various levels of water saturation deficit (WSD) in leaves at 15%, 30%, and 50%. In barley subjected to drought, most significant changes occurred under a slight dehydration level at 15%. It was observed that the gene expression of 9-cis-epoxycarotenoid dioxygenases, enzymes involved in ABA biosynthesis, increased significantly, and led to a higher concentration of ABA. This was most likely the result of an increase in the gene expression and enzyme activity of L-cysteine desulfhydrase, which is responsible for H2S synthesis. Our results suggest that the differential water deficit in leaves underlies the activation of an appropriate defence, with ABA metabolism at the centre of these processes. Furthermore, at 15% WSD, a dominant contribution of H2O2-dependent signalling was noted, but at 30% and 50% WSD, significant NO-dependent signalling occurred.


Introduction
In response to abiotic stress, plants trigger complex signalling cascades, including synthesis of reactive oxygen species (ROS), reactive nitrogen species (RNS) and reactive sulphur species (RSS) [1,2]. These signalling molecules regulate numerous physiological processes and enable plants to adapt to adverse conditions [3]. ROS, RNS and RSS (headed by hydrogen sulphide, H 2 S) engage in various molecular processes, including the modification of proteins containing thiol groups (-SH), which are particularly vulnerable to oxidation. One of the key elements of antioxidative defence is glutathione. The level of reduced glutathione (GSH), which is involved in the protection of -SH groups against oxidation, increases in plants under abiotic stress [4]. The level of GSH and the activities of the ascorbate-GSH cycle enzymes on the one hand are associated with ROS detoxification, but on the other hand, they are involved in redox signalling under drought stress [5]. Until now, the role of ROS and RNS during plant response against abiotic stress is better discerned than RSS, and thus there is a greater need to research the latter molecules.
One of the roles of H 2 S as a signalling molecule is the post-translational modification of reactive cysteine residues in proteins by their S-sulfhydration. H 2 S converts -SH to persulfide bridges (-SSH), leading to functional changes in proteins, which is an analogous mechanism to the S-nitrosylation induced by nitric oxide (NO) or other oxidative

Response of Barley to Water Deficiency
Water saturation deficit was marked for all the harvested plants. The control plants showed WSD at the level of 8%. Groups of plants treated with drought were harvested at WSD, which reached 15%, 30% and 50%, respectively. Figure 1 shows a canonical variate analysis (CVA) of the variability of the nine traits of four WSDs in terms of the first two canonical variates. In the graph, the coordinates of the point for particular WSDs are the values for the first and second canonical variates, respectively. The first two canonical variates accounted for 99.96% of the total multivariate variability between the individual WSDs ( Figure 1). variates accounted for 99.96% of the total multivariate variability between the individua WSDs ( Figure 1). The greatest variation in terms of all the nine traits jointly measured with Mahalanob distances was found for the control and 15% WSD (distance between them amounted to 364. 1 The greatest similarity was found between 15% WSD and 50% WSD (38.2) ( Table 1). Lipid peroxidation during a water deficit was estimated as the 2-thiobarbituric aci reactive substances (TBARS) content in the well-watered control and plants dehydrate to 15%, 30% and 50% WSD (Figure 2A). The TBARS content increased linearly durin dehydration, and it was the highest at the greatest dehydration at 50% WSD of the leave reaching 33.7 nmol of TBARS per g of fresh weight (FW), which was 6.6 times higher tha the lipid peroxidation in control plants.
Similarly, the concentration of NO2 − increased under dehydration, however reached the highest value of 26.5 µM NO2 − per g of FW at 30% of WSD and remained at comparable level at 50% of WSD ( Figure 2B). The content of NO3 − ions at 15% and 30% o WSD was comparable to the control, although the concentration of these ions was signi icantly lower at 30% WSD than at 15%. A significant increase in concentration was ob served only with deep dehydration at 50% WSD, which led to an increase of 52% of th NO3 − ions level in comparison to plants grown under control conditions ( Figure 2C).
The concentration of H2O2 rapidly increased under drought stress, and its highes values were recorded at the lowest dehydration level at 15% WSD in the initial phase o dehydration and was over 5.65 times higher than in the control plants. The H2O2 conten The greatest variation in terms of all the nine traits jointly measured with Mahalanobis distances was found for the control and 15% WSD (distance between them amounted to 364.1). The greatest similarity was found between 15% WSD and 50% WSD (38.2) ( Table 1). Lipid peroxidation during a water deficit was estimated as the 2-thiobarbituric acid reactive substances (TBARS) content in the well-watered control and plants dehydrated to 15%, 30% and 50% WSD ( Figure 2A). The TBARS content increased linearly during dehydration, and it was the highest at the greatest dehydration at 50% WSD of the leaves reaching 33.7 nmol of TBARS per g of fresh weight (FW), which was 6.6 times higher than the lipid peroxidation in control plants.
Similarly, the concentration of NO 2 − increased under dehydration, however it reached the highest value of 26.5 µM NO 2 − per g of FW at 30% of WSD and remained at a comparable level at 50% of WSD ( Figure 2B). The content of NO 3 − ions at 15% and 30% of WSD was comparable to the control, although the concentration of these ions was significantly lower at 30% WSD than at 15%. A significant increase in concentration was observed only with deep dehydration at 50% WSD, which led to an increase of 52% of the NO 3 − ions level in comparison to plants grown under control conditions ( Figure 2C).
The concentration of H 2 O 2 rapidly increased under drought stress, and its highest values were recorded at the lowest dehydration level at 15% WSD in the initial phase of dehydration and was over 5.65 times higher than in the control plants. The H 2 O 2 content remained high during the whole stress period, although the significant drop was noticed at 30% WSD (by 30% in comparison with 15% WSD) ( Figure 2D). remained high during the whole stress period, although the significant drop was noticed at 30% WSD (by 30% in comparison with 15% WSD) ( Figure 2D).

Variability of ABA and GSH/GSSG Contents on Dehydration
The level of ABA increased linearly with the intensity of dehydration. At the greatest dehydration at 50% WSD, its level increased more than six times compared to the control plants ( Figure 3A). The GSH content in the plant's response to dehydration was dependent on the WSD level. Reduced glutathione concentration increased by more than double in the first phase of drought stress at 15% WSD, as compared to the control. At 30%, a significant drop in the GSH content was observed. It decreased 12.8-fold in comparison to 15% WSD. At 50% WSD, a GSH content reached a concentration that was comparable with to control plants ( Figure 3B). Unlike the reduced form of glutathione, its oxidized form (GSSG) remained constant at all levels of tissue dehydration ( Figure 3B).  and dehydrated to 15%, 30%, and 50% WSD plants. Vertical bars indicate means (n = 9) ± standard deviation. Letters indicate Fisher's least significant differences, which were estimated at a significance level of α = 0.05.

Variability of ABA and GSH/GSSG Contents on Dehydration
The level of ABA increased linearly with the intensity of dehydration. At the greatest dehydration at 50% WSD, its level increased more than six times compared to the control plants ( Figure 3A). The GSH content in the plant's response to dehydration was dependent on the WSD level. Reduced glutathione concentration increased by more than double in the first phase of drought stress at 15% WSD, as compared to the control. At 30%, a significant drop in the GSH content was observed. It decreased 12.8-fold in comparison to 15% WSD. At 50% WSD, a GSH content reached a concentration that was comparable with to control plants ( Figure 3B). Unlike the reduced form of glutathione, its oxidized form (GSSG) remained constant at all levels of tissue dehydration ( Figure 3B).

Variability of ABA and GSH/GSSG Contents on Dehydration
The level of ABA increased linearly with the intensity of dehydration. At the greatest dehydration at 50% WSD, its level increased more than six times compared to the control plants ( Figure 3A). The GSH content in the plant's response to dehydration was dependent on the WSD level. Reduced glutathione concentration increased by more than double in the first phase of drought stress at 15% WSD, as compared to the control. At 30%, a significant drop in the GSH content was observed. It decreased 12.8-fold in comparison to 15% WSD. At 50% WSD, a GSH content reached a concentration that was comparable with to control plants ( Figure 3B). Unlike the reduced form of glutathione, its oxidized form (GSSG) remained constant at all levels of tissue dehydration ( Figure 3B).  . Effect of drought on abscisic acid (ABA) (A) and reduced (GSH) and oxidized (GSSG) glutathione (B) contents in the control and dehydrated to 15%, 30%, and 50% WSD plants. Vertical bars indicate means (n = 9) ± standard deviation. Letters indicate Fisher's least significant differences, which were estimated at a significance level of α = 0.05.

DES and NR Activities
The DES activity significantly increased at 15% of WSD by over 90% in comparison to plants grown in optimal conditions. At higher dehydration, the DES activity slightly decreased, however it remained higher by approximately 50% than in the control at both levels of dehydration: 30% and 50% of WSD ( Figure 4A). The NR activity was the same as the DES activity that reached the highest value at WSD 15%, and although it was higher by 41% than in the control, the observed difference was not statistically significant. A significant decrease in NR activity was observed at 30% WSD and it was 83% lower than at the initial stage of stress. At 50% WSD, NR activity returned to a level comparable to that of the control plants ( Figure 4B). bars indicate means (n = 9) ± standard deviation. Letters indicate Fisher's least significant differences, which were estimated at a significance level of α = 0.05.

DES and NR Activities
The DES activity significantly increased at 15% of WSD by over 90% in comparison to plants grown in optimal conditions. At higher dehydration, the DES activity slightly decreased, however it remained higher by approximately 50% than in the control at both levels of dehydration: 30% and 50% of WSD ( Figure 4A). The NR activity was the same as the DES activity that reached the highest value at WSD 15%, and although it was higher by 41% than in the control, the observed difference was not statistically significant. A significant decrease in NR activity was observed at 30% WSD and it was 83% lower than at the initial stage of stress. At 50% WSD, NR activity returned to a level comparable to that of the control plants ( Figure 4B). . Changes in the L-cysteine desulfhydrase (A) and nitrate reductase (B) activities in control and dehydrated to 15%, 30%, and 50% WSD plants. Vertical bars indicate means (n = 9) ± standard deviation. Letters indicate Fisher's least significant differences, which were estimated at a significance level of α = 0.05.

Analysis of Gene Expression of NCED1, NCED2, ABA8'OH1, ABA8'OH2, UGT1 and DES
Expression of genes encoding 9-cis-epoxycarotenoid dioxygenase was observed at the highest level in barley leaves at 15% WSD ( Figure 5A,B). The increase of transcript levels in comparison to the control was five-fold and two-fold higher for NCED1 and NCED2, respectively. With increasing WSD, a decrease in the expression of these genes was observed. Expression profiles of genes encoding ABA 8'-hydroxylase, responsible for ABA catabolism, were different. The highest level of the ABA8'OH1 transcript was observed in the control, and a decrease in the expression of this gene was noted in the barley leaves under drought ( Figure 5C). In turn, ABA8'OH2 expression increased in leaves with increasing WSD ( Figure 5D). The level of the UGT1 transcript, encoding the enzyme involved in the formation of ABA-GEs, was highest in leaves at 15% and 30% WSD ( Figure  5E). The expression of the gene encoding DES was the highest in barley leaves at 15% WSD and an over two-fold increase in the transcript level of this gene was observed, as compared to the control ( Figure 5F). Expression of genes encoding 9-cis-epoxycarotenoid dioxygenase was observed at the highest level in barley leaves at 15% WSD ( Figure 5A,B). The increase of transcript levels in comparison to the control was five-fold and two-fold higher for NCED1 and NCED2, respectively. With increasing WSD, a decrease in the expression of these genes was observed. Expression profiles of genes encoding ABA 8'-hydroxylase, responsible for ABA catabolism, were different. The highest level of the ABA8'OH1 transcript was observed in the control, and a decrease in the expression of this gene was noted in the barley leaves under drought ( Figure 5C). In turn, ABA8'OH2 expression increased in leaves with increasing WSD ( Figure 5D). The level of the UGT1 transcript, encoding the enzyme involved in the formation of ABA-GEs, was highest in leaves at 15% and 30% WSD ( Figure 5E). The expression of the gene encoding DES was the highest in barley leaves at 15% WSD and an over two-fold increase in the transcript level of this gene was observed, as compared to the control ( Figure 5F).

Correlation Analysis
The highest significant positive correlation was observed for ABA content and lipid peroxidation (TBARS) as well as for lipid peroxidation and nitrite ions content (correlation strength 0.96 at a significance level of α = 0.001), following correlation between DES activity and hydrogen peroxide concertation (strength 0.92 at a significance level of α = 0.001) and ABA content and nitrite ions content (strength 0.88 at a significance level of α = 0.001). A slightly lower positive correlation was noticed between reduced glutathione content and nitrate reductase activity (strength 0.73 at a significance level of α = 0.01) and between the ABA content and nitrate ions level (strength 0.58 at a significance level of α = 0.05). Apart from positive correlations, significant negative correlations were also observed. Such correlations were observed for the concentration of nitrite ions and the content of reduced glutathione as well as for the concentration of nitrite ions and the activity of nitrate reductase (strength −0.62 and −0.61, respectively, at a significance level of α = 0.05) ( Figure 6). Transcript levels of NCED1 (A), NCED2 (B), ABA8 OH1 (C), ABA8 OH2 (D), UGT1 (E) and DES (F) in the control and dehydrated plants to 15%, 30%, and 50% WSD. Vertical bars indicate means ± standard deviation of at least three technical replicates of semi-quantitative reverse transcription-PCR performed for each of the two biological replicates.

Correlation Analysis
The highest significant positive correlation was observed for ABA content and lipid peroxidation (TBARS) as well as for lipid peroxidation and nitrite ions content (correlation strength 0.96 at a significance level of α = 0.001), following correlation between DES activity and hydrogen peroxide concertation (strength 0.92 at a significance level of α = 0.001) and ABA content and nitrite ions content (strength 0.88 at a significance level of α = 0.001). A slightly lower positive correlation was noticed between reduced glutathione content and nitrate reductase activity (strength 0.73 at a significance level of α = 0.01) and between the ABA content and nitrate ions level (strength 0.58 at a significance level of α = 0.05). Apart from positive correlations, significant negative correlations were also observed. Such correlations were observed for the concentration of nitrite ions and the content of reduced glutathione as well as for the concentration of nitrite ions and the activity of nitrate reductase (strength −0.62 and −0.61, respectively, at a significance level of α = 0.05) ( Figure 6).

Discussion
Adequate water content in plant tissues is necessary for the proper plant growth and development, and in the case of arable crops, to produce high and good-quality crops. Maintaining the optimal water content in tissues requires coordinating plant growth with the water availability for the roots and the rate of water loss in the transpiration process. When this balance is disturbed because of soil drought, acclimatization responses are activated in the plant to adapt to unfavourable conditions, including e.g., closing the stomata, synthesis of compounds that protect cell structures and ensure the maintenance of the integrity of biological membranes [16]. The main research hypothesis of this study is that the varied level of hydration in spring barley leaves translates into contrasting responses that are triggered by water deficit in the plant in response to drought stress. We decided to verify this hypothesis, which we forged based on our previous research on the response of cereals to drought, using molecular and biochemical methods. For our research, we selected parameters that are known to be strongly associated with plant defence responses to diverse types of environmental stresses (abiotic and biotic). We have paid particular attention to the regulation of ABA content in leaves, which is known to be "the hormonal heart" of the plant response to drought stress [17].

Discussion
Adequate water content in plant tissues is necessary for the proper plant growth and development, and in the case of arable crops, to produce high and good-quality crops. Maintaining the optimal water content in tissues requires coordinating plant growth with the water availability for the roots and the rate of water loss in the transpiration process. When this balance is disturbed because of soil drought, acclimatization responses are activated in the plant to adapt to unfavourable conditions, including e.g., closing the stomata, synthesis of compounds that protect cell structures and ensure the maintenance of the integrity of biological membranes [16]. The main research hypothesis of this study is that the varied level of hydration in spring barley leaves translates into contrasting responses that are triggered by water deficit in the plant in response to drought stress. We decided to verify this hypothesis, which we forged based on our previous research on the response of cereals to drought, using molecular and biochemical methods. For our research, we selected parameters that are known to be strongly associated with plant defence responses to diverse types of environmental stresses (abiotic and biotic). We have paid particular attention to the regulation of ABA content in leaves, which is known to be "the hormonal heart" of the plant response to drought stress [17].
A rapid production of ROS, including H 2 O 2 , is one of the earliest activated plant responses against water deficit. Under drought condition, this molecule is among the others synthesized to perform the function of a secondary messenger in ABA-regulated stomata closure [18,19]. While analysing our results, we found several significant correlations between the analysed biochemical parameters. These dependencies indicate a specific phytohormone and signal molecule crosstalk during spring barley's response to water deficit. This is especially noticeable at 15% WSD. This would be coherent with the regularity that usually, in the first phase (an alarm phase) of the plant stress response (in this particular case 15% WSD), there is an intense liberation of signalling molecules to trigger defence responses against the stress detected by the plant.
It was noticed that, because of the increasing water deficit in barley leaves, and thus the incremental intensity of oxidative stress manifested by an elevated level of TBARS (markers of lipid oxidation), an enhanced content of H 2 O 2 and ABA and the enzyme activity of DES catabolizing L-cysteine to H 2 S were observed. This specific biochemical triad (H 2 S, H 2 O 2 and ABA) activated in barley plants struggling with progressive dehydration of leaf tissues may significantly promote the induction and development of drought acclimatization responses dependent on these three signalling molecules. It is highly likely that the end physiological effect is the closure of the stomata by barley plants, which reduces water loss through transpiration, thus maintaining a relatively optimal water level in the leaves is possible and water stressed barley plants can more effectively deal with the progressive secondary stress, namely oxidative stress. This claim is confirmed by numerous studies showing that H 2 S induced stomata closure by mediating in the ABA-dependent signalling, and what is especially important is that H 2 S stimulated the synthesis of H 2 O 2 in stomata, which led to stomata closure as well [7]. It has also been shown that the level of glutathione (which is part of the non-enzymatic antioxidant machinery) and its dependent responses to abiotic stresses are regulated by H 2 S [20,21]. We noted that GSH concentration increased more than double in the first phase of drought stress at 15% WSD compared to the control plants, thus suggesting its beneficial effect on the effectiveness of maintaining the redox balance in barley leaf cells at this stage of the response to drought.
As already mentioned above, ABA is a regulatory hub of phytohormone and other signalling molecule levels during plant water deficit stress responses, thus we took a closer look at the expression level of genes that participate in ABA metabolism. Many of the mechanisms of plant response to water scarcity depend on the ABA content in plant organs. In the tissues of plants undergoing drought, a several, or even several dozen, increase in the content of this phytohormone is observed [22]. Indeed, in barley leaves, the ABA content increased linearly with the exacerbation of dehydration, which allowed for the continued functioning of plants that were exposed to the lack of water.
Until recently, it was postulated that soil drought stimulates the biosynthesis of ABA in the roots, which is then transported through the xylem to shoots and leaves, where it reaches the guard cells and induces the closure of the stomata [23]. However, as studies have shown, when Arabidopsis roots and shoots were subjected to water deficit separately, an increase in ABA content was only observed in the shoots, indicating that they are the main sites for ABA synthesis in response to water deficit [24]. We observed that, in the leaves of barley plants subjected to drought, increased ABA biosynthesis takes place, which was indicated by an increased level of transcripts of genes encoding 9-cis-epoxycarotenoid dioxygenase-NCED1 and NCED2. It can be assumed that NCED1, whose increase in gene expression was higher in stressed plants at 15% WSD compared to the control, than that of NCED2, played a more significant role in the accumulation of ABA in the leaves of barley plants growing in the conditions of water deficit. Next, with increasing WSD, a decrease in the expression of these genes occurred. Similarly, in wheat leaves, dehydration caused a significant increase in the transcript level of one of the genes encoding NCED-TaNCED1 [25]. In turn, among the five genes encoding NCED in the Arabidopsis, expression of only one, AtNCED3, was strongly induced by drought, and overexpression of this gene resulted in increased ABA accumulation and increased drought tolerance [26]. Thus, in both cereals and model plants, the water deficit induced the expression of genes encoding NCED, hence this enzyme plays a key role in the regulation of ABA biosynthesis in response to drought.
The physiological effects mediated by ABA require the maintenance of an appropriate content of this phytohormone in plant tissues. Hence, in plants subjected to various stress factors, often intensive ABA biosynthesis, it is accompanied by an increase in the transcript level of genes encoding enzymes involved in ABA catabolism, wherein conversion of the ABA to an inactive form takes place due to ABA8'OH activity [27]. The genes encoding this enzyme belong to the cytochrome P450 family [28]. Four genes encoding ABA8'OH were found in the Arabidopsis. These genes showed tissue-specific patterns of expression that were altered by various environmental stimuli [27]. It has been shown that the transcript level of genes encoding ABA8'-OH increased in leaves of dehydrated Arabidopsis plants [29]. In turn, the ABA8'OH2 expression increased in leaves of barley plants with progressive water shortage, suggesting that ABA8'OH2 plays an essential role in the breakdown of ABA molecules during barley response to drought.
In leaves of barley plants growing under drought at 15% and 30% WSD, intensive ABA biosynthesis seems to coexist with the formation of inactive ABA-GEs, as evidenced by the increased level of the UGT1 transcript compared to control plants. ABA-GEs can be a transport form of this phytohormone and can also be stored in plant tissues [22,30]. ABA-GE are hydrolysed by β-glucosidase, resulting in the liberation of biologically active ABA molecules [31]. The hydrolysis of ABA-GE to ABA occurs in one step and compared to de novo ABA biosynthesis, it can lead to an increase in the content of this phytohormone in plant tissues much faster [32]. ABA-GEs formation and deconjugation may therefore play a vital role in increasing local ABA concentration in response to stress factors. The elevated level of the UGT1 transcript in the leaves of barley plants subjected to drought may indicate the accumulation of ABA in a conjugated form, which is a reserve of this phytohormone in the event of a deepening water deficit, which can be quickly activated in the process of deconjugation. Similarly, an increase in the transcript level of the UGT-encoding gene was observed in drought-treated Vigna angularis [33].
Plants need nitrogen (N) for developmental processes and defence against abiotic stress. A series of reactions catalysing by oxidoreductases (NO 3 − →NO 2 − →NH 4 + ), enabling N assimilation and biosynthesis of amino acids by transamination, is one of the main points of regulation of defence responses of plants against stresses [34]. The first of these reactions is catalysed by NR. We would like to bring to attention here that NO 3 − conversion to NO 2 − by NR is a restrictive step in N assimilation [35]. Thus, increased NR activity in barley leaves at 15% WSD, although the observed difference was not statistically significant in comparison with control plants, may promote the de novo synthesis of nutrients. Considering the participation of NR in response to abiotic stress, it is impossible not to emphasis its other physiological role, namely the ability to synthesise NO. We only mention this and hypothesize that the synthesis of NO by NR may take place, but we have no direct experimental evidence for it, so this issue requires further investigation soon.
Another interesting aspect related to NO metabolism is the observed coincidence revealed by statistical analysis and showing a significant positive correlation (strength 0.96 at a significance level of α = 0.001) for lipid peroxidation and nitrites, major end products of NO turnover. In addition, in barley leaves at 30% and 50% WSD, this trend was accompanied by a reduced level of H 2 O 2 compared to samples collected at 15% WSD. Thinking about the interpretation of these results, we conducted experiments to describe equivalent results in the field of vertebrate biochemistry. Van der Vliet et al. [36] revealed that NO 2 − can be oxidized by the haem peroxidases, such as horseradish peroxidase, myeloperoxidase (MPO), and lactoperoxidase (LPO), in the presence of H 2 O 2 , to nitrogen dioxide radical, which can promote tyrosine nitration. These authors concluded that NO 2 − , both at physiological and pathological contents, can be a substrate for the mammalian MPO and LPO and that production of nitrogen dioxide radicals by peroxidase-catalysed oxidation of nitrites can be a supplementary pathway inducing cytotoxicity or defence response associated with NO signalling. However, Byun et al. [37] showed that nitrogen dioxide radical produced by the MPO-H 2 O 2 − NO 2 − system was responsible for lipid peroxidation of low-density lipoprotein. We cannot rule out that this phenomenon, known from research on animal models, has been detected by us in our plant model. This is supported by the fact that the horseradish peroxidase mentioned above is also capable of oxidizing NO 2 − to nitrogen dioxide radical in the presence of H 2 O 2 . We have already undertaken research efforts to verify whether this mechanism actually works effectively in plants.

Plant Material
All experiments were conducted on spring barley (Hordeum vulgare L.) 'Airway'. It is a variety of fodder barley with high prolificacy and protein content in kernels. It is recommended for medium-intensive and intensive cultivation technology on moderatequality soils. Spring barley 'Airway' is highly resistant to Rhynchosporium secalis, Drechslera sorokiniana and Blumeria graminis f. sp. hordei and good to Pyrenophora teres and Puccinia hordei. Plants are of medium height with particularly good resistance to lodging. The advantage of this variety is high resistance to peduncle breakage.
Initially, grains were soaked in water for 2 h and then surface decontaminated in 5% sodium hypochlorite solution. Then, in order to eliminate microbial contamination, the grains were incubated for one hour on a magnetic stirrer in 0.2% Plant Preservative Mix-ture™ (PPM) (Plant Cell Technologies, Inc., Washington, DC, USA). The decon-taminated grains were placed in a petri dish on filter paper soaked in 0.2% PPM. After a two-day incubation of the grains at 23 • C in the darkness, they were transferred to pots (four grains per pot) filled with commercial potting mix. Plants were grown in a growth chamber at a temperature of 23 ± 2 • C and 18 • C with a 16 h day (100 ± 25 µmol m −2 s −1 light)/8 h night photoperiod, respectively, and under 70-80% humidity.
After fourteen days of growth under optimal conditions, drought was induced by the cessation of water supply to the pots. The research material consisted of plants growing under control conditions and subjected to drought stress at various stress levels at 15, 30% and 50% WSD. To put it in detail, after the cessation of irrigation, the plants were incubated in the growth chamber for a further six days, during which time the WSD in the leaves reached 15%, after another three days 30% and after another two days 50%. Appropriate control plants were sampled. All plants (control and stressed) were in the same stage of development, namely at the beginning of tillering.

Measurement of WSD
The dehydration of spring barley leaves was determined as the percentage of WSD, calculated according to Turner [38], as follows: Water saturated mass − actual fresh mass Water saturated mass − dry mass × 100% Water saturated mass was measured after submersion of leaves overnight in water in the dark, and the dry mass was determined after drying the tissue at 80 • C for 16 h. Actual fresh mass weight was measured right after sample collection.

L-Cysteine Desulfhydrase (DES) Activity
The activity of DES was determined using a modified method of Alvarez et al. [39].

Lipid Peroxidation
Lipid peroxidation was measured by the Heath and Packer [41] method. 100 mg of leaves were homogenized in 0.5 mL of 0.1% (w/v) trichloroacetic acid (TCA). The homogenate was centrifuged for 10 min at 13,200× g at 4 • C and the supernatant was collected. 1.5 mL of 0.5% (w/v) 2-thiobarbituric acid dissolved in 20% TCA was added to 0.5 mL of the supernatant. The mixture was incubated at 95 • C for 25 min and the reaction was stopped by transferring samples on ice. After cooling, 200 µL of each mixture was taken and pipetted into Nunc U-bottom 96-well plate (Thermo Scientific, Waltham, MA, USA). The absorbance of the supernatant was measured at 532 nm, with a reading at 600 nm subtracted (nonspecific turbidity) using a Varioskan LUX Multimode Microplate Reader (Thermo Scientific, Waltham, MA, USA). The TBARS content was determined using an extinction coefficient of 155 mM −1 cm −1 (to quantify lipid peroxides). It was expressed as nmol of TBARs per gram of FW.

Nitrate Reductase (NR) Activity
The activity of nitrate reductase was determined by the modified Jaworski method [42]. 100 mg of the leaves were homogenized in liquid nitrogen and 1 mL of 0.1 M Tris-HCl buffer (pH 7.5) was added. The homogenate was centrifuged for 15 min at 13,200× g at 4 • C and the supernatant was collected. To 45 µL of the reaction mixture which consisted of 5% (v/v) propanol and 30 mM KNO 3 dissolved in 0.1 M phosphate buffer (pH 7.5) and 5 µL of supernatant was added. The mixture was incubated for 30 min at 30 • C and 50 µL of a solution consisting of 1% (w/v) sulphanilamide and 0.02% (v/v) N-(1-naphthyl)ethylenediamine hydrochloride in 3 M HCl was added. The absorbance at 540 nm was measured in Nunc U-bottom 96-well plate (Thermo Scientific, Waltham, MA, USA) on a Varioskan LUX Multimode Microplate Reader (Thermo Scientific, Waltham, MA, USA). Protein content was determined by the Bradford [40] method. The NR activity was calculated using a standard curve and expressed in nM NO 2 − per min per mg of protein.

Glutathione Content
Reduced glutathione was assayed following the method described by Gronwald et al. [43]. 200 mg of leaves was ground to a powder in liquid nitrogen and then homogenized in 7.57 mM sodium ascorbate and incubated for 30 min at 4 • C. Next, the homogenate was centrifuged for 15 min at 16,000× g at 4 • C. The supernatant was used for the total glutathione (GSH + GSSG) assay and for the oxidized glutathione (GSSG) assay. N-ethylmaleimide (NEM), which masks thiol groups in GSH, was used to determine GSSG. 1 mL of the reaction mixture consisted of 0.4 mL of extract and 10 mM NEM dissolved in 125 mM sodium-phosphate buffer (pH 7.5) containing 5 mM EDTA. The mixture was incubated at room temperature for 70 min and then NEM was removed by extraction with ether, which was repeated 5 times. Glutathione was determined as described by Ellman [44] with the 5,5 -dithiobis (2-nitrobenzoic acid) (DTNB). The reaction mixture for the assay of GSSG or total glutathione consisted of 10 mM phosphate buffer (pH 7.0), 0.5 mM NADPH, 0.12 units glutathione reductase and supernatant. An increase in the absorbance at 412 nm was measured after the addition of DTNB. The total glutathione (GSH + GSSG) and GSSG contents were calculated using a standard curve and expressed in the µM of glutathione per gram of FW.

Determination of Endogenous ABA Content
The endogenous ABA concentration was measured using the Phytodetek ABA enzyme immunoassay test kit (Agdia, Inc., Elkhart, IN, USA). A sample of 250 mg of frozen barley leaves was homogenized in liquid nitrogen, then 1.5 mL of extraction buffer (80% methanol, 2% acetic acid and 20 mg/l butylated hydroxytoluene) was added to the powdered tissue. The ABA was extracted and described by Fidler et al. (2018) [45]. The ABA content was measured in triplicate.

Total RNA Extraction
Total RNA was isolated from the barley leaves using the GeneMATRIX Universal RNA Purification Kit (EURx, Gdańsk, Poland), including the RNase-free DNase I (EURx, Gdańsk, Poland) treatment according to the manufacturer protocol.

Semi-Quantitative Reverse Transcription-PCR (sqRT-PCR)
The transcript levels of the genes encoding enzymes engaged in ABA metabolism and DES in barley leaves were determined using Titanium One-Step RT-PCR Kit (Clontech Laboratories, Mountain View, CA, USA). sqRT-PCRs were performed with gene-specific primers and aliquots of 100 ng of total RNA were used as templates in each reaction. As an internal control amplification of the 18S rRNA gene was performed. The reactions were performer under following conditions: 60 min at 50  Figure S1 and Table S1 in Supplementary Materials.

Statistical Analysis
Three independent experiments in three replicates (n = 9) were performed. The normality of the distribution of the nine traits, i.e., H 2 O 2 , DES, NO 3 − , NO 2 − , TBARS, NR, GSH, GSSG and ABA, was verified with Shapiro-Wilk's normality test to check whether the analysis of variance (ANOVA) met the assumption that the ANOVA model residuals followed normal distribution. The homogeneity of variance was evaluated using Bartlett's test. Box's M test was used to check multivariate normality and homogeneity of variancecovariance matrices. All the traits had normal distribution. A one-way (WSD) multivariate analysis of variance (MANOVA) was performed. Following this, one-way analyses of variance (ANOVA) were performed to verify the null hypotheses of a lack of WSD effect on the nine observed traits, independently for each one. The arithmetic means, minimal and maximal values as well as standard deviations were calculated. Moreover, Fisher's least significant differences (LSDs) were estimated at a significance level of α = 0.05. The relationships between observed traits were estimated using Pearson's linear correlation coefficients. The results were also analysed using multivariate methods. A canonical variate analysis was applied to present a multi-trait assessment of the similarity of the tested WSDs in a lower number of dimensions with the least possible loss of information. The Mahalanobis distance was suggested as a measure of "polytrait" WSDs similarity [46], the significance of which was verified by means of critical value Dα called "the least significant distance" [47]. The Mahalanobis distances were calculated for all pairs of WSDs.
The GenStat v. 18 statistical software package (VSN International) was used for the analyses. All statistical analyses are presented in Table S2 in Supplementary Materials.

Conclusions
In conclusion, our study indicates that the specific responses of spring barley plants against drought are related to the intensity of water stress in planta measured as WSD in leaves. In barley plants exposed to drought, the most noticeable molecular and metabolic changes took place under slight dehydration level at 15% WSD. It was noted that the gene expression of 9-cis-epoxycarotenoid dioxygenases, enzymes involved in ABA biosynthesis, was significantly enhanced, which led to an elevated level of ABA. This was accompanied by an increase in the gene expression and enzyme activity of DES, which is responsible for H 2 S production. These discoveries are consistent with the outlook that during the first phase (an alarm phase) of the plant ecophysiological reaction (in our experimental case at 15% WSD) there is an intense biosynthesis of signalling molecules to trigger defence against stressful environmental conditions. Furthermore, results indicate a dominant contribution of H 2 O 2 signalling at 15% WSD, while at 30% and 50% WSD a significant role is played by NO signalling as increased level of nitrites (major end products of NO metabolism) was noted.