Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparative physiological and proteomic analysis indicates lower shock response to drought stress conditions in a self-pollinating perennial ryegrass

  • Fatemeh Raeisi Vanani,

    Roles Methodology, Software

    Affiliation Department of Plant Science, Faculty of Science, Shahrekord University, Shahrekord, Iran

  • Leila Shabani ,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    shabani-l@sku.ac.ir, lshabani@gmail.com

    Affiliations Department of Plant Science, Faculty of Science, Shahrekord University, Shahrekord, Iran, Research Institute of Biotechnology, Shahrekord University, Shahrekord, Iran

  • Mohammad R. Sabzalian,

    Roles Conceptualization, Writing – review & editing

    Affiliation Department of Agronomy and Plant Breeding, College of Agriculture, Isfahan University of Technology, Isfahan, Iran

  • Fariba Dehghanian,

    Roles Data curation, Formal analysis

    Affiliation Department of Cell and Molecular Biology & Microbiology, Faculty of Biological Science and Technology, University of Isfahan, Isfahan, Iran

  • Lisa Winner

    Roles Data curation

    Affiliation Core Facility Proteomics, Center for Biological Systems Analysis (ZBSA), University of Freiburg, Freiburg, Germany

Abstract

We investigated the physiological and proteomic changes in the leaves of three Lolium perenne genotypes, one Iranian putative self-pollinating genotype named S10 and two commercial genotypes of Vigor and Speedy, subjected to drought stress conditions. The results of this study indeed showed higher RWC (relative water content), SDW (shoot dry weight), proline, ABA (abscisic acid), nitrogen and amino acid contents, and antioxidant enzymes activities of S10 under drought stress in comparison with the two other genotypes. A total of 915 proteins were identified using liquid chromatography-mass spectrometry (LC/MS) analysis, and the number of differentially abundant proteins between normal and stress conditions was 467, 456, and 99 in Vigor, Speedy, and S10, respectively. Proteins involved in carbon and energy metabolism, photosynthesis, TCA cycle, redox, and transport categories were up-regulated in the two commercial genotypes. We also found that some protein inductions, including those involved in amino acid and ABA metabolisms, aquaporin, HSPs, photorespiration, and increases in the abundance of antioxidant enzymes, are essential responses of the two commercial genotypes to drought stress. In contrast, we observed only slight changes in the protein profile of the S10 genotype under drought stress. Higher homozygosity due to self-pollination in the genetic background of the S10 genotype may have led to a lower variation in response to drought stress conditions.

1 Introduction

Cool-season grasses, emerging from and adapted to cool climatic conditions are economically and ecologically among the most important species due to their critical role in CO2 assimilation as well as forage production for livestock farming utilization [1]. Perennial ryegrass (Lolium perenne L.) is the most frequently cultivated species considered properly for farming [2] as it combines high growth rates and supplies the nutritive value of fodder. However, Norris [3] and Turgeon [4] mentioned that perennial ryegrass is a drought-susceptible grass species, although there seems to be a noticeable genotypic variation for drought stress responses.

Drought stress is the main abiotic element limiting plant growth and crop output throughout the world [5], which is roughly limiting plant production in 25% of the world's land. Due to its geographical location, Iran has an arid (65%) to semi-arid (25%) climate condition and is generally considered a dry country. Therefore, drought stress is the most critical and common environmental stress that limits agricultural production in the country and reduces the efficiency of using semi-arid and rain-fed areas. In cool-season grass species, drought stress adversely affects leaf water potential, photosynthetic and carbohydrate accumulation, respiration, production of reactive oxygen species, lipid peroxidation, and denaturation of proteins (reviewed by Loka et al. [6]). However, many cool-season grasses can resist moderate drought and maintain their aerial growth to avoid and tolerate leaf dehydration. Also, they may be equipped with the strategies facilitating their survival of severe drought, mainly associated with dehydration avoidance and tolerance, which primarily occurs in meristematic tissues [7]. In some species and genotypes of grasses, summer dormancy is also another strategy that confers efficient survival of meristematic tissues through dehydration avoidance and tolerance [8].

The grass family (Poaceae) exhibits different reproductive strategies among more than 10,000 species, some promoting self-pollination, and others cross-pollination, due to mechanisms such as self-incompatibility and being dioecious [9]. Lolium perenne is commonly considered as an out-crossing species [10]. This breeding system maintains a high genetic variation within the progenies and keeps the plasticity of the variety against environmental instabilities [11, 12]. However, this will also prohibit the repeatability of results in scientific studies since keeping the background of genetic material is impossible. Recently, Torkian et al. [13] reported that some of the Iranian perennial ryegrasses have appropriate reproductive characters of seed production as well as turf characteristics, which are also self-pollinating. Therefore, they may be considered as a new source of germplasm for future turfgrass breeding programs. Before this, it needs to confirm the potential of this newly emerged material, especially their performance in drought stress environments.

Under drought stress conditions, numerous alterations are induced in plants' physiological and molecular characteristics. Physiological characters of plants, particularly those related to plant water status, have a great deal of importance in the growth and development of plants during water-stress conditions. In the present study, relative water content (RWC) has been used as an example of the characters which determine drought-stress tolerance of Lolium genotypes. Biomass is a crucial parameter that influences plant growth; therefore, the accumulation and distribution of dry mass are significant considerations when investigating the effect of drought stress on plant growth. The role of glutathione in the antioxidative defense system also provides a rationale for its use as a stress marker. In general, an initial stress response is related to changes in the glutathione redox state. By this, acclimation is marked by increased glutathione concentrations, increase in related enzyme activities, and a more reduced redox state of glutathione.

Among phytohormones, abscisic acid is the most important one that confers abiotic stress tolerance in crop plants [14, 15]. In drought conditions, ABA content in plants increases considerably, inspiring stress-tolerance effects that help plants adapt and survive under these stressful situations [16]. The 'omics' approaches are often applied to reveal molecular changes due to drought stress and to investigate samples collected from different populations or upon different physiological states, aiming to find molecules that would differentiate between these classes of samples. These molecular studies in L. perenne for understanding some underlying molecular mechanisms in response to drought stress have been somewhat limited. In other species, recently proteome analyses have opened the way to the identification of proteins involved in drought stress tolerance in rice, Oryza sativa [17], barley, Hordeum vulgare [18, 19], Stipa purpurea [20], Poa pratensis [21] and wheat, Triticum aestivum [22]. The effort described here related to the comparative physiological and proteomic analysis was applied to investigate the molecular events–due to drought stress tolerance- in three perennial ryegrass genotypes including S10 as a newly introduced breeding line from Iran with self-pollinating reproductive habit with superior seed and turf characteristics compared to Speedy (Speedy Green) and Vigor as cross-pollinating commercial genotypes.

2 Materials and methods

2.1 Plant material and experimental design

Lolium perenne plant material of this experiment consisted of one well-performed genotype of S10, selected from the previous experiment [13] and two commercial populations of Vigor and Speedy Green or in short Speedy from Barenbrug (Barenbrug Co., the Netherland). Thus, three genotypes in this research were examined in a pot culture in a greenhouse according to a completely randomized design with three replications. Greenhouse conditions consisted of 16h/8h of light/dark photoperiod at a temperature range from 18.3 to 25.5°C. Ten tillers (with approximately equal sizes) of each genotype previously grown in a field were selected and transplanted into every 60 pots (size: 20*15 cm) filled with a 1:3 (v/v) mixture of sand and soil. The soil was sandy loam in texture with pH = 7.78 and Electrical Conductivity of 1.34 Ds.m-1. All three genotypes were then propagated for two months before being used in the experiment. Then, plants were regularly irrigated and grown in controlled greenhouse conditions for two weeks before starting the drought stress treatment. For this treatment, pots were divided into two groups: well-watered plants were irrigated every two days (control), and water-stressed plants. Drought stress was induced on plants by limiting irrigation to ensure water depletion to 20% of the field capacity of the soil using theta probe (AT Delta-T Devices SM300, Cambridge, England) before re-irrigation for 30 days. After seven days of the first period of drought stress treatment, leaves were cut off and stored at -70°C for later protein extraction. Biochemical measurements for the three genotypes of Lolium under drought stress were carried out on 1 g of fresh leaves. Antioxidant enzymes activities, ABA hormone and glutathione contents were measured one week after treatment; and proline, total N and amino acid contents, and shoot dry weight were calculated three weeks after the initial drought stress. At the end of the experiment, the harvested plants in pots were washed with distilled water, followed by surface drying with filter paper. For each plant, the leaves and roots were separated. Then leaves and roots were dried for 15 min at 105°C and then at 65°C to a constant weight in an oven. The dry weights were then determined.

2.2 Relative water content (RWC)

Leaf relative water contents (RWC) were calculated with the method of Barrs and Weatherley [23]. The relative leaf water content (RWC) was calculated as 100*(FW-DW)/(TW-DW), where FW = fresh weight, DW = dry weight, and TW = turgid weight. TW was measured after saturating the water content of leaf discs for 24 hours. DW was measured after desiccating leaf discs at 60°C to air-dry state.

2.3 Proline content

A 0.1 g sample of the fresh leaf was ground in liquid nitrogen from all the Lolium genotypes which were subjected to drought treatment. Leaf proline content was analyzed at the 3rd week post-stress treatment. The content of free proline in leaves was determined as described by Troll and Lindsley [24].

2.4 Estimation of total nitrogen

For estimation of total N, 1 g of fine-ground leaf dry samples was digested with sulfuric acid, and assays were carried out according to the Kjeldahl method [25].

2.5 Extraction and measurement of ABA

ABA was extracted according to the procedure developed by Kelen et al. [26] with some modifications. Leaf samples (1g) were used for extraction. Leaf tissue was ground with a mortar and a pestle in liquid nitrogen and homogenized in 10 cm3 extraction solution containing butylated hydroxy-toluene (0.25 g) and ascorbic acid (0.24 g) dissolved in 90% methanol and stirred overnight at 4°C. The extract was filtered through a Whatman filter, and methanol was evaporated under vacuum. Then, pH of the aqueous phase was adjusted to 8.5 with 0.1 M phosphate buffer and then partitioned with ethyl acetate three times. After removal of the ethyl acetate phase, pH of the aqueous phase was adjusted to 2.5 with 0.2 M HCl. The solution was partitioned with ethyl acetate three times and then passed through anhydrous sodium sulfate. After, the ethyl acetate phase was evaporated under vacuum, and the dry residue containing hormones was dissolved in 0.5 cm3 of HPLC grade methanol and stored in vials at 4°C.

An aliquot (20 mm3) of the filtered samples was injected into a Waters Symmetry-C18 chromatographic column (250 × 4.6 mm) (USA) with isocratic elution at a flow rate of 0.7 cm3 min-1 at 25°C using a mobile phase containing acetic acid solution (0.2%) and methanol (100%) (50:50 v/v). Detection of ABA was carried out at 265 nm by a UV detector (Waters 2487) with known concentrations of (±)-abscisic acid (Sigma A1049).

2.6 Antioxidant enzymes and glutathione assays

Activities of catalase (CAT), superoxide dismutase (SOD), ascorbate peroxidase (APX), glutathione reductase (GR) and glutathione S-transferase (GST) enzymes were detected in leaf tissue extracts spectrophotometrically as described in our previous work [27]. Glutathione (GSH) content of the leaves was measured according to May and Leaver [28].

2.7 Amino acids assay

Amino acids were extracted from one gram of dry leaves according to the procedure described by Lisiewska et al. [29]. The chromatographic separation of 13 amino acids was achieved by the ion exchange chromatography in an amino acid analyzer (ARACUS, Germany). The identification of amino acids was carried out based on the retention time at the chromatographic separation. All determinations were carried out in two replications for each sample.

2.8 Proteomic analysis

Approximately 100 mg of plant material (fresh weight) was ground in liquid nitrogen using a pre-cooled mortar and pestle. The dry powder was transferred to 1.5 cm3 of precipitation solution (10% TCA w/v, 0.07% DTT w/v), according to Yang et al. [20]. The extract was left to precipitate at -20°C overnight. Precipitated proteins were then pelleted at 16,000 ×g for 15 min at 4°C. The pellet was washed at least three times with acetone containing 0.07% (w/v) DTT at -20°C for 30 min, then, the mixture was vortexed and centrifuged at 16,000 ×g for 20 min, at 4°C. The pellet was recovered and air-dried at room temperature and dissolved in 1 cm3 of 1x SDS sample buffer [30]. The solutions were vigorously vortex-mixed and kept at room temperature for one h and centrifuged at 14,000 ×g for 10 min at 4°C. The supernatant was transformed into a new tube and stored at 4°C overnight and run on a 5% stacking gel (in an SDS PAGE system with an equal volume of stacking (5%) and resolving (10%) gel) according to Raorane et al. [30] (S1 Fig). The single bands formed in stacking gel, containing condensed and purified total soluble proteins, were stained with coomassie brilliant blue (G-250) as described by Neuhoff et al. [31] and used in the following steps.

2.8.1 Mass spectrometry and data analysis.

After gel staining using coomassie, the gel was cut in 10 slices. Next, gel digestion was performed using trypsin, and peptides were extracted and desalted. The protocol reported by Rappsilber et al. was used for peptide extraction [32]. Then, liquid chromatography-mass spectrometry (LC/MS) analysis was performed on an LTQ Orbitrap XL (Thermo Fisher Scientific, Bremen, Germany) coupled to an Eksigent 2D nano flow-HPLC equipped with in-house packed C18 columns of approximately 20 cm length (Reprosil-Pur 120 C-18-AQ, 3 μm, Ammerbuch, Germany) without a pre-column. Because of an unavailable database for the proteins of L. perenne and limited protein database of Lolium genus, Brachypodium distachyon protein database was selected for the identification of proteins in this study. The analysis of MS raw data was performed using MaxQuant (v. 1.5.3.30 [33]) against the UniProt Brachypodium distachyon protein sequence database (released June 20, 2018; 44,786 entries in total). A false discovery rate of 1% was applied for both protein lists and peptide-spectrum matches (on modified peptides, separately). Analysis of protein groups was performed with the Perseus software [34]. Normalized ratios were log2-transformed, and a mean log2 was calculated across all three replicates. Finally, proteins with a p-value below 0.05 and fold change higher or lower than +1.5 and -0.5, respectively, were considered. The differentially abundant proteins were used for COG (Cluster of Orthologous Groups of proteins) category and annotation.

2.8.2 Bioinformatics analysis.

The proteomic data visualization in schematic metabolism pathways, as well as functional classification of proteins, was performed by using the MapMan (Version 3.6.0RC1, downloaded from http://mapman.gabipd.org/ Web Site) software program.

2.9 Statistical analyses

Statistical analyses were carried out with three biological and two technical replicates for proteomic analyses. The results of protein abundance were statistically evaluated by student's t-test by using the Multi Experiment Viewer (MEW) software considering P<0.05 as the critical significance level. Relative water content, shoot dry weight, total nitrogen, proline and glutathione contents, antioxidant enzymes activities, and ABA content were analyzed using ANOVA by the SAS software [35]. Significant differences between means were defined by LSD post-hoc test at P<0.05. PCA (Principal Component Analysis) was conducted utilizing STATGRAPHICS Centurion ver. 17 software (Statpoint Technologies, Inc., Warrenton, Virginia, USA).

3 Results

3.1 Physiological characteristics of ryegrass genotypes under drought stress

The results of statistical analyses for RWC, shoot dry weight (SDW), proline, ABA, total N, and glutathione contents of the three genotypes of Lolium perenne are presented in Table 1. RWC content of the three genotypes decreased significantly under drought stress (P<0.05). However, the RWC content of the S10 genotype was much higher compared to the other genotypes under both stress and control conditions. Also, the SDW of the genotypes decreased significantly under drought stress (P<0.05). However, like RWC, SDW was significantly higher in S10 than that of the other two genotypes under both stress and control conditions. At stress conditions, dry weights of Vigor, Speedy, and S10 genotypes decreased by 9.7%, 12.5%, and 5.5%, respectively, compared to the control condition.

thumbnail
Table 1. Relative water content (RWC), shoot dry weight (SDW), proline content, ABA concentration, total Nitrogen, and glutathione content of three genotypes of Lolium perenne under drought stress.

Values represent means of n = 3 replicates.

https://doi.org/10.1371/journal.pone.0234317.t001

Proline content of the three genotypes increased under drought stress (P<0.05). However, proline accumulation in S10 genotype was much higher than Vigor and Speedy genotypes under both drought and control conditions. The ABA concentration in the genotypes also increased under drought stress (P<0.05) and, like proline, the concentration of this hormone in S10 was significantly higher than that in the two commercial genotypes under both drought and control conditions. In all genotypes, the content of leaf total N decreased under drought stress (P<0.05). In contrast, drought stress significantly increased glutathione content in Vigor and S10 genotypes.

All of the tested antioxidant enzymes, including APX, SOD, CAT, GR, and GST, showed similar activity changes under drought stress in the three genotypes. However, the induction of APX activity in Vigor and Speedy was much higher than that of the S10 genotype under drought stress. In contrast, CAT and GR activities were higher in S10 genotype under drought stress conditions (Fig 1).

thumbnail
Fig 1. Antioxidant enzyme activities of three genotypes of Lolium perenne under drought stress.

Bars represent standard error. Mean values followed by the same letter are not significantly different (p < 0.05). Different letters indicate significant differences by LSD test at P < 0.05, n = 3. V = vigor genotype, S = speedy genotype.

https://doi.org/10.1371/journal.pone.0234317.g001

The content of amino acids, including aspartic acid, threonine, serine, glutamic acid, valine, methionine, isoleucine, leucine, tyrosine, arginine, phenylalanine, histidine and lysine of the three genotypes was measured three weeks after drought stress (Fig 2). The results showed that the majority of the amino acid content of all genotypes increased under drought conditions. However, the total amino acid content of S10 was much higher than that of Vigor and Speedy under both control and drought stress conditions. Under stress conditions, the total amino acid content of S10 was 262 (mg/100 g DW) as against 219 and 217 (mg/100 g DW) in Vigor and Speedy genotypes, respectively.

thumbnail
Fig 2. Amino acids content including Aspartic acid, Threonine, Serine, Glutamic acid, Valine, Methionine, Isoleucine, Leucine, Tyrosine, Arginine, Phenylalanine, Histidine and Lysine of three genotypes of Lolium perenne under drought stress.

https://doi.org/10.1371/journal.pone.0234317.g002

3.2 Proteomic characteristics of ryegrass genotypes under drought stress

The small grass species, Brachypodium distachyon, is potentially an ideal model plant system for grass research. Therefore, the Brachypodium protein database was used for protein identification to dissect tolerance to drought stress in Lolium perenne genotypes. Based on the Brachypodium protein database, a total of 915 proteins were identified in each genotype. The number of proteins whose relative abundance level was significantly altered under drought stress conditions was 467, 456, and 97 in Vigor, Speedy, and S10 genotypes, respectively.

Both “up-regulated” and “on (New)” proteins are considered together as “differentially increased proteins." Similarly, the term “differentially decreased proteins” included both “down-regulated” and “off (disappeared)” proteins. Among the proteins, 230, 282 and 40 differentially increased proteins (fold change >1.5, P < 0.05) and 237, 174 and 57 differentially decreased proteins (fold change <0.5, P < 0.05) were found under drought stress compared to control condition in Vigor, Speedy, and S10 genotypes, respectively, (Fig 3, S1 Table).

thumbnail
Fig 3. The number of up and down-regulated proteins in drought-stressed seedlings of Vigor, Speedy, and S10 compared to control plants of each genotype.

“on (New)”: proteins in experimental plants that were not found in controls. “off (disappeared)”: Some proteins found in controls were absent in the experimental samples.

https://doi.org/10.1371/journal.pone.0234317.g003

3.3 Functional annotation of identified and differentially abundant proteins

To simplify the analysis of protein expression, differentially abundant proteins were classified into 12 COG categories using Mapman ontology and the Brachypodium mapping file (Fig 4). Among 230 differentially increased proteins in the Vigor genotype, the most are involved in carbon and energy metabolism, protein, redox, and transport categories. Similarly, in Speedy genotype, among 282 differentially increased proteins, the most are involved in carbon and energy metabolism, protein, redox, stress, and transport categories. In contrast, in the S10 genotype, most of the proteins are involved in carbon and energy metabolism, and protein and stress categories decreased under drought stress (Fig 4).

thumbnail
Fig 4. COG annotation of differentially abundant proteins of three genotypes of Lolium perenne (fold increase>1.5, fold decrease<0.5, P < 0.05).

A: Vigor genotype, B: Speedy genotype, C: S10 genotype.

https://doi.org/10.1371/journal.pone.0234317.g004

Based on the COG annotation, drought stress induced changes in the abundance of proteins which are constituting the photosystem I (PSI) and PSII reaction centers and PSI and PSII polypeptide subunits. We found that proteins constituting reaction centers of the two photosystems (no. 6, 9, 17, 31, 54) increased in both Vigor and Speedy genotypes. In contrast, some proteins of reaction centers decreased in Vigor (no. 455) and S10 (no. 54) genotypes.

Several polypeptide subunits of photosystems highly decreased in the three genotypes of Lolium under drought stress. We observed a drought-induced increase in the levels of small and large subunits of RuBisCO in Vigor and Speedy genotypes, whereas, in S10, such a change was not detected. In the Vigor genotype, RuBisCO activase exhibited drought-induced accumulation. The Expression level of the protein sedoheptulose-bisphosphatase, which participates in the Calvin cycle, was up-regulated under drought in Speedy genotype (Fig 5A, S1 Table). The chloroplastic isoform of fructose-1,6-bisphosphate aldolase (FBA) was down-regulated in Vigor and Speedy genotypes; however, no changes were detected in S10 genotype. In contrast, the cytoplasmic isoform of FBA was up-regulated in Vigor and Speedy genotypes, whereas it was down-regulated in the S10 genotype.

thumbnail
Fig 5.

Drought-induced changes in the accumulation levels of selected proteins involved in carbon and energy metabolism category (A), cell wall, lipid metabolism, and N and amino acid metabolism categories (B), redox category (C) and hormone metabolism, stress, protein, and transport categories (D) identified in leaves of three genotypes of Lolium perenne. V = vigor genotype, S = speedy genotype.

https://doi.org/10.1371/journal.pone.0234317.g005

Alpha and beta subunits of pyruvate dehydrogenase complex, malate dehydrogenase, malic enzyme, citrate synthase, succinyl-CoA ligase, succinate dehydrogenase, isocitrate dehydrogenase, and aconitate hydratase, which are involved in the glycolysis pathway and Krebs cycle, increased in both Vigor and Speedy genotypes but did not change in the S10 genotype. However, the NADP dependent malic enzyme (NADP-ME) increased in all three genotypes.

Based on the COG annotation, the differentially abundant proteins involved in starch and sucrose metabolisms were up- or down-regulated in the three genotypes under drought stress. Regarding sucrose synthase-related proteins, we found that in S10 genotype, sucrose synthase (SS) increased, while fructose-1,6-bisphosphatase showed reduced levels (Fig 5A, S1 Table). We also found the up-regulation of starch-degradation enzymes, starch-related R1 protein, and starch phosphorylase in the three genotypes of Lolium in drought stress conditions (Table 1). On the contrary, the isoforms of ADP-glucose pyrophosphorylase involved in the pathway of starch synthesis were up-regulated only in Vigor and Speedy genotypes. In the S10 genotype, no changes were detected in some of the isoforms, and other isoforms showed down-regulated expression levels (Fig 5A, S1 Table).

Based on the COG annotation, the differentially abundant proteins in the cluster of cell walls, including cellulose, pectinase, β-xylosidase and β-1, 4-glucanase enzymes were down-regulated in Vigor and Speedy genotypes (Fig 5B, S1 Table). A significant increase in the protein abundance of arabinogalactan-rich proteins (AGPs) has also been found in the S10 genotype.

Our proteomic analysis revealed that drought-induced changes in concentration of proteins involved in lipid metabolism and two isoforms of the phospholipase enzyme increased in both commercial genotypes. Interestingly, opposite to these genotypes, proteins involved in lipid metabolism did not change or even were reduced in the S10 genotype under drought stress conditions (Fig 5B, S1 Table).

In the functional group of nitrogen and amino acid metabolism, some of the differentially abundant proteins increased in Vigor and Speedy genotypes and the proteins involved in the metabolism of methionine, cysteine, aspartate, glycine, and proline were up-regulated in these two commercial genotypes. Δ-1-pyrroline-5-carboxylate synthase (P5CS) enzyme, which catalyzes the first two steps in proline biosynthesis in plants, did not exhibit significant changes in the three genotypes of Lolium. However, drought-induced increase in the level of pyrroline-5-carboxylate reductase (P5CR), which is a part of the L-proline biosynthesis pathway, was significantly higher in Vigor and Speedy genotypes. Also, the abundance of 1-pyrroline-5-carboxylate dehydrogenase enzyme, which is a part of the pathway of L-proline degradation into L-glutamate, significantly increased in these two genotypes. We observed drought-induced reduction in the levels of nitrite reductase 1, glutamine synthetase, and glutamate dehydrogenase (involved in nitrogen flow) in Vigor and Speedy genotypes (Fig 5B, S1 Table).

Concerning the ROS-generation proteins, we found two differentially increased proteins, glycolate oxidase, and oxalate oxidase, in the two commercial genotypes. However, they decreased or did no change in drought-stressed S10 genotype. Under drought stress, plants respond to ROS-over abundance in cells by producing enzymatic and nonenzymatic ROS scavenging antioxidants. The increased level of dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR), catalase, ascorbate peroxidase (APX) and SOD was observed in leaves of both Vigor and Speedy genotypes. However, GR and glutathione peroxidase was up-regulated only in the Vigor genotype. Under drought stress, the up-regulation of glutathione S-transferase was found in the three genotypes of Lolium. The other hydrogen peroxide-decomposing enzymes, glutathione peroxidase, and peroxiredoxin and thioredoxin reductase increased only in Vigor and Speedy genotypes in response to drought stress (Fig 5C, S1 Table).

Given that ABA metabolism is a crucial aspect of the plant response to drought stress, proteins involved in ABA metabolism were searched in the regulated proteins. We found that two isoforms of 9-cis epoxy carotenoid dioxygenase (NCED1 and NCED4) increased during drought stress in Vigor and Speedy genotypes; however, the S10 genotype presented no alteration in this respect (Fig 5D, S1 Table).

Based on the COG annotation, the levels of four heat shock proteins (HSPs) (HSP60, chaperonin-60ALPHA, chaperonin 20 and DNAJ) increased under drought stress only in Vigor and Speedy genotypes; and no protein abundance changes of HSP-like proteins were detected in S10 genotype (Fig 5D, S1 Table).

Also, based on the COG annotation, some proteins of the transport system were up-regulated only in the two commercial genotypes. The PIP1;4 and PIP2B aquaporin proteins (AQPs) were up-regulated in Speedy and Vigor genotypes; however, the up-regulation level of these proteins was higher in Speedy genotype. Furthermore, the PIP2A protein and two isoforms of PIP were also up-regulated in the Speedy genotype (Fig 5D, S1 Table).

3.4 Multivariate data analysis

Principal Component Analysis (PCA) was conducted to determine the correlation between the three genotypes of ryegrass, drought stress, and all the examined characteristics. Noteworthy, the PCA analysis displayed a high correlation between biomass, total N, and RWC content (Fig 6); however, no positive correlation between biomass and total N with amino acids total content and antioxidant enzymes was found. The separation of S10 genotype from the two commercial genotypes (Vigor and Speedy) can be easily deduced from the data provided. Also, we found that the S10 genotype in drought conditions had a combination of higher dry weight (second component) and antioxidant enzymes (first component) than those in the other two genotypes. Therefore, under drought, better biomass and antioxidant potential along with higher drought tolerance are undoubtedly ascribed to this genotype in comparison to the two commercial genotypes. The PCA results also revealed the separation of control and drought conditions for all tested genotypes and indicated that drought stress was particularly distinctive by induction in the production of antioxidant enzymes. This PCA analysis may suggest that the defense strategy utilized by the two commercial genotypes to assist plants against drought is different from the mechanism used by the well-performed S10 genotype.

thumbnail
Fig 6. The two-dimensional biplot of PCA (principal component analysis) indicating the correlation between the three genotypes of ryegrass, drought stress, and all examined characteristics.

C: control condition, D: drought condition, V: Vigor genotype, Sp: Speedy genotype, S10: S10 genotype, CAT: Catalase activity, APX: Ascorbate peroxidase activity, SOD: Superoxide dismutase activity, GST: Glutathione-S-transferase activity, GR: Glutathione reductase activity, G: Glutathione content, ABA: ABA concentration, W: RWC content, Dry: shoot dry weight, P: proline content, N: Nitrogen total, aa: total amino acids content.

https://doi.org/10.1371/journal.pone.0234317.g006

4 Discussion

The decisive role of plant genotype for improving drought tolerance in Lolium perenne was reflected in higher RWC, dry weight and N content in S10 genotype under drought stress than those in Vigor and Speedy genotypes, indicating increased drought tolerance in this genotype of Lolium.

In this study, the accumulation of proline and some other amino acids increased in the three genotypes of Lolium under drought stress. However, the content of proline, aspartate, glutamate, phenylalanine, and tyrosine was higher in S10 genotypes than that in the other two genotypes (Fig 2). These results were consistent with the previous studies on tall fescue that showed up-regulation of proline and some amino acids under drought stress [36, 37]. Therefore, the results of this study indicate an increase in physiological response capacity of S10 genotype to drought stress, which may enable the plant to have excellent stomatal control and accordingly reduce water utilization.

Based on the Brachypodium protein database, we identified 467, 456, and 97 differentially abundant proteins in Vigor, Speedy, and S10 genotypes, respectively, under drought stress (S1 Table). These results demonstrate that drought stress coincides with changes in proteins. We also found that the number of differentially decreased proteins (57) was higher than that of differentially increased proteins (40) in the S10 genotype, showing that protein breakdown is the primary function during drought stress in this genotype. Some previous studies also showed that drought resistance was involved in an alteration in the gene expression responsive to a decrease in the transcript abundance [38, 39].

In the present study, we found that the increased proteins in the two commercial genotypes were mostly involved in carbon and energy metabolism, photosynthesis, TCA cycle, redox, and transport categories. However, in the S10 genotype, they were in N and amino acid metabolism category (Fig 5B, S1 Table). Also, in Vigor and Speedy genotypes, the decreased proteins were related to cell wall, N, and amino acid metabolism category and redox. However, in the S10 genotype, they were in carbon and energy metabolisms.

Preventing photosynthesis due to the closure of stomata is one of the destructive effects of drought stress which in turn results in a total decrease in the expression level of the proteins involved in photosynthesis [40, 41]. In contrast, in our study, the majority of photosynthesis-related proteins increased as an influence of drought stress in the two commercial genotypes. This increase in abundance mainly involved carbon fixation enzymes, including RuBisCO, RuBisCO activase, FBA, glyceraldehyde-3-phosphate dehydrogenase, sedoheptulose-1,7-bisphosphatase, and phosphoribulokinase. While light reaction-involved proteins were down-regulated under drought in all three genotypes of Lolium; therefore, the increase in carbon fixation enzymes may remain useless.

Damage to photosystems usually occurs when plants are subjected to drought stress [42]. This statement was confirmed in our study by the identification of reduction in light-harvesting complexes in Vigor and Speedy genotypes, and somehow in S10 genotype. The decrease in abundance of these polypeptides has been reported before in several species under drought stress [4345].

All of the enzymes involved in the Krebs cycle were enhanced under drought stress in the two commercial genotypes. The up-regulation of the proteins involved in the Krebs cycle under drought stress is in accordance with the results of previous researches [46, 47]. Enhancing TCA cycle activity may ensure energy for various processes towards tolerance to drought stress. The TCA pathway is one of the main pathways of energy production, and the up-regulation of the enzymes involved in this pathway can be due to the cell's demand for energy.

Sucrose synthase (SS), fructose-1,6-bisphosphatase, and ADP-glucose pyrophosphorylase are the key enzymes involved in the regulation of sucrose synthesis in the cytoplasm and starch synthesis in the chloroplast. Our data showed that under drought stress, the abundance of SS in the S10 genotype increased to 3-fold compared to the control condition. However, the abundance of fructose-1,6-bisphosphatase was reduced in this genotype. SS belongs to a large family of glycosyltransferase enzymes in plants that catalyzes the synthesis of sucrose from fructose and UDP-glucose, providing carbon for respiration and the biosynthesis of starch and cellulose [48].

During drought stress, biosynthesis of starch is generally emptied in several plant species [49, 50]. The ADP-glucose pyrophosphorylase enzyme plays a crucial role in the biosynthesis regulation of starch from glucose-1-phosphate. In the genotype S10, we found a reduction in the abundance of the ADP-glucose pyrophosphorylase. In contrast, the abundance of this protein and starch synthase increased in the two commercial genotypes of Lolium. There is evidence that at the beginning of drought stress, a short-term increase in starch synthesis may happen [51]. Then, the up-regulated expression levels of starch-degradation enzymes could lead to a reduced level of starch in the plants.

On the one hand, the plant requires energy to cope with the stress; on the other hand, the rate of photosynthesis in the plant decreases under stress conditions. Therefore, to provide the required energy for various activated pathways against stress, the plant starts metabolizing starch. Also, cell wall invertase, which is involved in sucrose degradation, decreased under drought stress in the two commercial genotypes of Lolium. Invertase catalyzes the hydrolysis of sucrose into glucose and fructose, plays a crucial role in primary metabolism and plant development. Down-regulation of cell wall invertase and sucrose synthase was reported from Piriformospora indica-inoculated barley exposed to drought stress compared to non-inoculated plants [52]. Our results imply a crucial role of SS and invertase in Lolium genotypes in the regulation of sugar biosynthesis when exposed to drought stress. Since sucrose is involved in the stability of proteins as a compatible osmolyte, thus in stressed plants, enhancement of sucrose may play an essential regulatory effect in the alleviation of damage [53].

The amounts of most of the proteins in the two commercial genotypes were down-regulated at the cell wall level in order to cope with drought stress. Plants may down-regulate some enzymes, such as cell wall hydrolysates, in order to reduce the hydrolysis of cell wall polysaccharides, which may be an approach to save their energy and preserve their carbohydrate reserves, and thus to survive under stress. An increase in cell wall synthesis, under drought stress conditions, is probably due to the increased mechanical strength as a strategy to reduce dehydration [54].

Arabinogalactan proteins (AGPs) belonging to the hydroxyproline-rich cell wall glycoprotein superfamily, are implicated in various aspects of plant growth and development, including cell differentiation, cell expansion and abiotic stress response through modulating cell wall expansion [55] and possible interaction of the plant with microbial endophytes [56]. We identified FLA11 as a fasciclin-like arabinogalactan protein that was up-regulated in the S10 genotype under drought stress. Similarly, the up-regulation of five FLA genes in Populus trichocarpa was reported under salt stress [57]. Accumulation of AGPs in the S10 genotype could be stimulating for extreme amounts of cell walls expansion under drought stress. The enlargement of cell wall, in turn, can cause a considerably higher water absorption expanding the cell volume of plants.

In plants, phospholipases hydrolyze membrane phospholipids to phosphatidic acid (a membranous second messenger molecule); the target of phosphatidic acid is ABI1 that is implicated in ABA signaling in stomatal responses [58]. The latest studies demonstrate that phospholipase plays an essential role in plant drought stress tolerance [59, 60,61]. In this study, we found that phospholipases were considerably accumulated in the two commercial genotypes of Lolium under drought stress. The activity and transcript levels of PLD (phospholipase D) in cowpea and peanut increased under drought stress in cultivars that were drought susceptible compared to drought-tolerant cultivars [62, 63]. Hence, we propose that phospholipases may be implicated in the membrane rearrangement and first happenings in ABA signal transduction in stomatal movement in Vigor and Speedy genotypes under drought stress.

Various studies have reported that drought stress alters the volume of free amino acids in plants' cells [64, 36]. These alterations can be due to increased protein degradation under stress [65, 66]. The participation of metabolic pathways of amino acids in growth and regulation of tolerance and adaption to stresses has been shown in plants. In this study, drought stress enhanced synthesis and degradative pathways of some amino acids, especially methionine, cysteine, aspartate, glycine, and proline. It strongly enhanced the synthesis of proline, as biochemical data also confirm that proline content was up-regulated in drought-stressed plants in all three genotypes (Table 1). Drought stress can induce proteolysis of protein and thereby leads to an accumulation of amino acids [67]. In our study, drought led to a rise in amino acid content in all three genotypes. The difference in amino acid contents resulted from drought stress is more likely to reflect the turnover of amino acids. The accumulation of P5CR, which participates in proline biosynthesis, was significantly higher under drought stress in the two commercial genotypes. Notwithstanding with no overall change at the protein level of P5CS in the examined genotypes, our biochemical analysis showed enhancement of proline level in leaves of all three genotypes. In the S10 genotype, we found a more significant increase than 2.1-fold in the accumulation of this amino acid, but only 1.2-fold alteration was recorded for Vigor and Speedy genotypes (Table 1).

Drought stress decreased photosynthesis, coupled with the reduction in nitrogen metabolism [68, 69]. Nitrogen assimilation occurs in plants via two key enzymes, glutamine synthetase (GS) and glutamate synthase (GOGAT). Singh et al. [70] showed that soil drying decreased the activities of these enzymes together with a decrease in N absorption and content of nitrogenous compounds in leaves of sunflower seedlings. Our results show that the abundance of GS and GOGAT decreased in the two commercial genotypes under stress, but no change was observed in the abundance of these two enzymes in the S10 genotype. Similarly, Ghaffari et al. [52] reported that drought stress declined NR, GS, and GOGAT enzymes in non-inoculated barley, but this was not the case for inoculated plants with Piriformospora indica. It was demonstrated in numerous studies that leaf nitrogen reduces progressively as the drought proceeds, and this can be correlated to photosynthetic machinery damage [71, 69]. Our biochemical data show that total nitrogen content decreased in Vigor and Speedy genotypes under drought stress, which confirms the lower activity of GS and GOGAT enzymes in this condition. Meanwhile, total nitrogen content did not change in the S10 genotype during drought stress, which is in accordance with the proteomic data. Numerous studies have reported crop productivity reduction due to drought stress as one of the most significant challenges [72, 73]. One of the reasons for the decreased growth rate of the plants exposed to drought stress can be due to the prevention of nitrogen fixation by drought stress. In this study, a down-regulation in the expression level of nitrite reductase (NR) protein was also detected in the two genotypes of Speedy and Vigor. The relative reduction of shoot dry weight under drought stress was similar for these two genotypes (about 11%), probably showing a lower NR activity in both genotypes. However, this reduction for the S10 genotype under drought was about 5.4 times smaller than the reduction in the two commercial genotypes.

Abscisic acid (ABA) is a phytohormone that plays a crucial role in plant tolerance to drought stress. ABA regulates many aspects of plant development under abiotic stresses, permitting the plant to resist drought conditions. Drought-induced increase in ABA content usually results in leaf stomata closure, stimulation of stress-related proteins, and provides different metabolites for the response to stress [74]. An increased ABA level has been formerly informed in several species under drought stress [20, 75,36]. We found an approximately 6-fold increase of ABA concentration in the two genotypes of Speedy and Vigor two days after drought stress. The ABA content is controlled by the relative values of biosynthesis, catabolism, conjugation, and its distribution through the plant. In our results, both the content of endogenous ABA and the expression of rate-limiting ABA biosynthetic protein (NCED) increased in the two commercial genotypes under drought stress. Similar proteomic analysis findings have been documented for several plant species, for instance, in Phaseolus [76], Arabidopsis [77], and Stipa [20] showing the up-regulation of NCED genes under drought stress. However, our findings showed that although no protein abundance changes of ABA biosynthesis-related enzymes were identified in the S10 genotype, the endogenous ABA increased in the leaves of this genotype, suggesting that the increase could be possibly due to the decreased catabolism of this hormone.

ROSs are generated as a routine by-product of aerobic metabolism in plant cells. Besides, drought stress is recognized to accelerate the accumulation of ROS in plant cells, lastly resulting in cell damage. Therefore, ROS levels need to be actively controlled for the protection of plants from drought-induced oxidative stress. The redox signaling also regulates a large number of transcripts that are necessary for stress acclimatization. Ascorbic acid (AsA), which is a constituent of the AsA-glutathione (GSH) cycle [78], maintains redox homeostasis in the plant cell. In order to maintain the antioxidative potential of AsA, rapid restitution of AsA is arranged by DHAR and MDHAR [7880]. Activated enhancement of DHAR and MDHAR in these stress responses has been found in various plants. In our study, up-regulation of DHAR and MDHAR in the two commercial genotypes of Lolium under stress suggests that these genotypes may increase AsA content by increasing DHAR and MDHAR enzymes. Also, glutathione (GSH; γ-glutamyl-cysteinyl-glycine) plays an essential function in sequestration and transport of reduced sulfur participating directly or indirectly in the detoxification of ROS [81].

Regarding glutathione-related proteins, we identified that glutathione reductase and glutathione peroxidase were up-regulated only in commercial genotypes; but not in the S10 genotype. However, the up-regulation of glutathione S-transferase was observed in all three genotypes. We observed that reduced glutathione content increased in the three genotypes of Lolium, probably due to the action of glutathione-related proteins. Thioredoxins (Trxs) are well-conserved disulfide reductases that regulate the redox state of target proteins. In plants, Trxs are proved to be essential for plant tolerance to oxidative stresses [82]. In Vigor and Speedy genotypes, the up-regulation of thioredoxin M-type was detected under drought stress. Peroxiredoxins (Prx) are a family of thiol dependent peroxidases found in plant cells, which use thioredoxin as their electron donor for the catalysis of H2O2. They can act as redox sensors, signal transducers, and molecular chaperones through the alterations in their oligomeric structures [83]. Also, we found up-regulation of the two isoforms of peroxiredoxin type 2 in Vigor and Speedy genotypes under drought stress. Generally, high numbers of detected leaf proteins related to ROS scavengings such as APX, SOD, GR and CAT, glutathione peroxidase, and peroxiredoxin and thioredoxin reductase were up-regulated. Up-regulation of most antioxidant enzymes has been formerly documented in different species under drought stress [20, 84,47]. However, in the S10 genotype under drought stress, the antioxidant system did not respond as vigorously as expected, except for SOD and GST enzymes, which were enhanced to play defensive roles against ROS. Ghaffari et al. [52] reported that four isoforms of GST and three isoforms of ascorbate peroxidase were up-regulated only in the non-inoculated barley; but not in the inoculated ones with mycorrhiza. Glutathione S-transferase was also up-regulated in all three genotypes of Lolium. Possibly because of the up-regulation of ROS-generating proteins in commercial genotypes, it seems that the amount of produced ROS in these genotypes was higher than that in the S10 genotype. Thus they increased the number of activated enzymes in order to maintain their homeostasis.

The up-regulated proteins included heat shock proteins (HSPs) and aquaporins (AQP), which are known to be critical for stress acclimatization [85] and water transport [86]. Heat shock proteins are implicated in plant acclimatization to different stresses [85, 87]. Members of various heat shock protein families play substantial roles in trafficking and folding of proteins and assembling of macromolecular in the cell [88]. Our results showed induction of heat shock protein 60, and increased abundance was found for the chaperonin-60ALPHA, chaperonin 20, and DNAJ, which is named HSP40, under drought treatment only in the two commercial genotypes. The expression of HSP has been indicated to increase in many plant species under drought stress [89, 20, 84]. Conversely, Ghaffari et al. [52] have reported that under drought stress, the expression of HSPs decreased in the non-mycorrhizal plants of barley. Our findings highlight that the HSPs may play an essential role in maintaining protein stability in some genotypes under drought stress.

Strikingly, in the two genotypes of Vigor and Speedy under drought stress, five plasma membrane intrinsic proteins were found, and all these exhibited increased levels. However, these proteins in the S10 genotype remained unchanged under stress. Several investigators have formerly found the up-regulation of aquaporins in response to water stress in different plant species [9092]. Also, the transcript level of aquaporin was enhanced by Trichoderma harzianum colonization in rice [93].

In our study, the evaluation of RWC content during drought stress indicated significantly higher values of RWC in the S10 genotype in comparison to the two commercial genotypes (Table 1). The better water status of the S10 genotype ensured that the expression of aquaporins gene is unnecessary in this genotype and permitted to maintain higher leaf net photosynthesis and water use efficiency than the two commercial genotypes when exposed to drought stress.

Comparably, low genetic gains characterize current population-based breeding methods for essential traits such as biomass yield in many plants. However, grass species have adopted a range of various breeding systems, some promoting self-pollination, some cross-pollination, and some asexual propagation. Our results demonstrated that the self-pollinating genotype of Lolium perenne (S10) might provide a better buffer in response to drought stress. Therefore, self-pollinating genotypes of the plant might be advantageous over cross-pollinating ones due to having more gain for desired traits, including drought stress resistance after fixation in genetic loci and selection for the well-performed characteristics [13]. In our study, it seems that the S10 genotype is a potential new genetic resource with a self-pollinating reproductive system that could be introduced to different genetic backgrounds for better resistance to biotic and abiotic stresses.

5 Conclusion

This study provides perceptions into the effects of drought stress on the abundance of proteins in the leaves of three genotypes of Lolium perenne. The results of this study indeed showed that RWC and SDW, proline, ABA, N, and amino acid contents, and antioxidant enzyme activities were significantly higher in the drought-tolerant genotype of S10 under drought stress in comparison with the two commercial genotypes. Results of the proteomics analysis revealed significant differences between the examined genotypes of Lolium under drought stress conditions (Fig 7), and this result seems to be related to the higher ability of the S10 genotype in response to drought stress. Intriguingly, we found only slight alterations in the protein profiles of S10 genotypes under drought stress, proposing that protein abundance was less influenced, possibly due to the tolerance in this genotype. That consequence may be related to the self-compatibility mechanism that can improve grass characteristics and alleviate some adverse effects of drought stress. Also, the results showed that acclimatization to drought stress in the two commercial genotypes involved specific responses that resembled between Vigor and Speedy. In this study, some underlying mechanisms were identified through which the plant responses may be improved to drought. These mechanisms included (1) significant increases in the abundance of antioxidant enzymes to scavenge ROS; and (2) changes in the expression of various proteins, including those involved in TCA cycle, amino acid and ABA metabolisms, aquaporins, and HSPs. These changes increased stress tolerance through ABA-mediated regulation of stomatal control to diminish water loss and assist in decreasing the adverse effects of drought by inhibiting damage to fundamental cell building blocks and improving osmoregulation. Additional analyses may be also required to clarify the root proteome that seems to participate in this response.

thumbnail
Fig 7. A schematic diagram of differences in drought response strategies between the S10 genotype and Vigor and Speedy genotypes of perennial ryegrass.

Abbreviations: PDC: pyruvate dehydrogenase complex, MDH: malate dehydrogenase, CS: citrate synthase, SCS: succinyl-CoA ligase, SDH: succinate dehydrogenase and IDH: isocitrate dehydrogenase.

https://doi.org/10.1371/journal.pone.0234317.g007

Supporting information

S1 Table. Proteins showing significant changes in responses to drought stress in Vigor, Speedy, and S10 genotypes.

B.C. is the abbreviation of the Bin Code (i.e., major functional categories). No.: identification number. ID: Δ: fold changes in drought-stressed plants (D) with respect to the control ones (C) (up: D/C, down:—C/D). new: not present in C; D.: disappeared, not present in drought stress.

https://doi.org/10.1371/journal.pone.0234317.s001

(DOCX)

S1 Fig. One-dimensional SDS-Polyacrylamide Gel Electrophoresis (1D SDS-PAGE) of total protein of three genotypes of Lolium perenne.

Left to right: Protein marker (lane 1), Control sample of Vigor genotype (lane 2,3), Drought sample of Vigor genotype (lane 4,5), Drought sample of Speedy genotype (lane 6), Control sample of Speedy genotype (lane 7), Control sample of S10 genotype (lane 8,9), Drought sample of S10 genotype (lane 10, 11).

https://doi.org/10.1371/journal.pone.0234317.s002

(TIF)

Acknowledgments

We acknowledge the University of Shahrekord, SKU, for supporting this research (grant number: 98GRD30M1032).

References

  1. 1. Falloon P, Betts R. Climate impacts on European agriculture and water management in the context of adaptation and mitigation—the importance of an integrated approach. Sci Total Environ. 2010; 408: 5667–5687. pmid:19501386
  2. 2. Furet PM, Berthier A, Decau ML, Morvan-Bertrand A, Prud'homme MP, Noiraud-Romy N, et al. Differential regulation of two sucrose transporters by defoliation and light conditions in perennial ryegrass. Plant Physiol Biochem. 2012; 61: 88–96. pmid:23085586
  3. 3. Norris I., 1985. Relationships between growth and measured weather factors among contrasting varieties of Lolium, Dactylis and Festuca species. Grass Forage Sci. 40, 151–159.
  4. 4. Turgeon A. Turfgrass Management. Copyright by Pearson Education. Inc. Uppere Saddle River, New; 2002.
  5. 5. Kramer P. Water relations of plants. Inc New York: Academic Press; 1983.
  6. 6. Loka D, Harper J, Humphreys M, Gasior D, Wootton-Beard P, Gwynn-Jones D, et al. Impacts of abiotic stresses on the physiology and metabolism of cool-season grasses: A review. Food Energy Secur. 2019; 8: e00152.
  7. 7. Norton MR, Malinowski DP, Volaire F. Plant drought survival under climate change and strategies to improve perennial grasses. Agron Sustain Dev. 2016; 36: 29.
  8. 8. Volaire F, Norton M. Summer dormancy in perennial temperate grasses. Ann. Bot. 2006; 98: 927–933. pmid:17028299
  9. 9. Connor HE. Breeding systems in the grasses: a survey. New Zeal J Bot. 1979; 17: 547–574.
  10. 10. Cornish MA, Hayward MD, Lawrence MJ. Self-incompatibility in ryegrass. Heredity. 1979; 43: 95–106.
  11. 11. Bellon MR, Adato M, Becerril J, Mindek D. Poor farmers’ perceived benefits from different types of maize germplasm: the case of creolization in lowland tropical Mexico. World Dev. 2006; 34: 113–129.
  12. 12. Ceccarelli S, Grando S, Maatougui M, Michael M, Slash M, Haghparast R, et al. Plant breeding and climate changes. J. Agric. Sci. 2010; 148: 627–637.
  13. 13. Torkian M, Sabzalian MR, Ehtemam MH. A simultaneous effect of selfing and Epichloë endophyte on forage, seed yield and turf characteristics of perennial ryegrass (Lolium perenne L.). Grass Forage Sci. 2019; 00: 1–12.
  14. 14. Shinozaki K, Yamaguchi-Shinozaki K. Molecular responses to dehydration and low temperature: differences and cross-talk between two stress signaling pathways. Curr. Opin. Plant Biol. 2000; 3: 217–223. pmid:10837265
  15. 15. Schroeder JI, Kwak JM, Allen GJ. Guard cell abscisic acid signaling and engineering drought hardiness in plants. Nature. 2001; 410: 327–330. pmid:11268200
  16. 16. Ng LM, Melcher K, Teh BT, Xu HE. Abscisic acid perception and signaling: structural mechanisms and applications. Acta Pharmacol. Sin. 2014; 35: 567–584. pmid:24786231
  17. 17. Mirzaei M, Soltani N, Sarhadi E, Pascovici D, Keighley T, Salekdeh GH, et al. Shotgun proteomic analysis of long-distance drought signaling in rice roots. J Proteome Res. 2011; 11: 348–358. pmid:22047206
  18. 18. Kausar R, Arshad M, Shahzad A, Komatsu S. Proteomics analysis of sensitive and tolerant barley genotypes under drought stress. Amino Acids. 2013; 44: 345–359. pmid:22707152
  19. 19. Wang N, Zhao J, He X, Sun H, Zhang G, Wu F. Comparative proteomic analysis of drought tolerance in the two contrasting Tibetan wild genotypes and cultivated genotype. BMC Genomics. 2015; 16: 432. pmid:26044796
  20. 20. Yang Y, Dong C, Yang S, Li X, Sun X, Yang Y. Physiological and proteomic adaptation of the alpine grass Stipa purpurea to a drought gradient. PLoS ONE. 2015; 10: e0117475. pmid:25646623
  21. 21. Xu C, Huang B. Differential proteomic responses to water stress induced by PEG in two creeping bentgrass cultivars differing in stress tolerance. J. Plant Physiol. 2010; 167: 1477–1485. pmid:20674080
  22. 22. Cheng L, Wang Y, He Q, Li H, Zhang X, Zhang F. Comparative proteomics illustrates the complexity of drought resistance mechanisms in two wheat (Triticum aestivum L.) cultivars under dehydration and rehydration. BMC Plant Biol. 2016; 16: 188. pmid:27576435
  23. 23. Barrs H, Weatherley P. A re-examination of the relative turgidity technique for estimating water deficits in leaves. Aust. J. Biol. Sci. 1962; 15: 413–428.
  24. 24. Troll W, Lindsley J. A photometric method for the determination of proline. J. Biol. Chem. 1955; 215: 655–660. pmid:13242563
  25. 25. Kjeldahl J. Neue Methode zur Bestimmung des Stickstoffs in organischen Körpern. Fresenius J Anal Chem. 1883; 22: 366–382.
  26. 26. Kelen M, Demiralay EC, Sen S, Özkan G. Separation of abscisic acid, indole-3-acetic acid, gibberellic acid in 99 R (Vitis berlandieri, Vitis rupestris) and rose oil (Rosa damascena) by reversed phase liquid chromatography. Turk J Chem. 2004; 28: 603–610.
  27. 27. Mirzahosseini Z, Shabani L, Sabzalian MR, Sharifi-Tehrani M. Neotyphodium endophytes may increase tolerance to Ni in tall fescue. Eur J Soil Biol. 2014; 63: 33–40.
  28. 28. May MJ, Leaver CJ Oxidative stimulation of glutathione synthesis in Arabidopsis thaliana suspension culture. Plant Physiol. 1993; 103: 157–166.
  29. 29. Lisiewska Z, Kmiecik W, Korus A. The amino acid composition of kale (Brassica oleracea L. var. acephala), fresh and after culinary and technological processing. Food Chem. 2008; 108: 642–648. pmid:26059143
  30. 30. Raorane ML, Narciso JO, Kohli A. Total soluble protein extraction for improved proteomic analysis of transgenic rice plant roots. In: MacDonald J, Kolotilin I, Menassa R, editors. Recombinant proteins from plants. Humana Press; 2016. p. 139–147.
  31. 31. Neuhoff V, Arold N, Taube DEW. Improved staining of proteins in polyacrylamide gels at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. Electrophoresis. 1988; 9: 255–62. pmid:2466658
  32. 32. Rappsilber J, Mann M, Ishihama Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat. Protoc. 2007; 2: 1896–906. pmid:17703201
  33. 33. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotech. 2008; 26: 1367–72.
  34. 34. Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods. 2016; 13: 731–40. pmid:27348712
  35. 35. SAS [computer program] (1999). Version 8. Cary, NC: SAS Institute Inc.
  36. 36. Man D, Bao YX, Han LB, Zhang XZ. Drought tolerance associated with proline and hormone metabolism in two tall fescue cultivars. HortSci. 2011; 46: 1027–1032.
  37. 37. Nagabhyru P, Dinkins RD, Wood CL, Bacon CW, Schardl CL. Tall fescue endophyte effects on tolerance to water-deficit stress. BMC Plant Biol. 2013; 13: 127. pmid:24015904
  38. 38. Ma Y, Shukla V, Merewitz EB. Transcriptome analysis of creeping bentgrass exposed to drought stress and polyamine treatment. PLoS ONE. 2017; 12: e0175848. pmid:28445484
  39. 39. Moon KB, Ahn DJ, Park JS, Jung WY, Cho HS, Kim HR, et al. Transcriptome profiling and characterization of drought-tolerant potato plant (Solanum tuberosum L.). Mol Cells. 2018; 41: 979–992. pmid:30396236
  40. 40. Patro L, Mohapatra PK, Biswal UC, Biswal B. Dehydration induced loss of photosynthesis in Arabidopsis leaves during senescence is accompanied by the reversible enhancement in the activity of cell wall beta-glucosidase. J Photochem Photobiol B. 2014; 137: 49–54. pmid:24793323
  41. 41. Ghotbi-Ravandi AA, Shahbazi M, Shariati M, Mulo P. Effects of mild and severe drought stress on photosynthetic efficiency in tolerant and susceptible barley (Hordeum vulgare L.) genotypes. J Agron Crop Sci. 2014; 200: 403–415.
  42. 42. Wang X, Cai X, Xu C, Wang Q, Dai S. Drought-responsive mechanisms in plant leaves revealed by proteomics. Int. J. Mol. Sci. 2016; 17: 1706.
  43. 43. Ye J, Wang S, Zhang F, Xie D, Yao Y. Proteomic analysis of leaves of different wheat genotypes subjected to PEG 6000 stress and rewatering. Plant Omics. 2013; 6: 286–94.
  44. 44. Ye T, Shi H, Wang Y, Yang F, Chan Z. Contrasting proteomic and metabolomic responses of Bermuda grass to drought and salt stresses. Front. Plant Sci. 2016; 7: 1694. pmid:27891145
  45. 45. Zhou S, Li M, Guan Q, Liu F, Zhang S, Chen W. Physiological and proteome analysis suggest critical roles for the photosynthetic system for high water-use efficiency under drought stress in Malus. Plant Sci. 2015; 236: 44–60. pmid:26025520
  46. 46. Deng X, Liu Y, Xu X, Liu D, Zhu G, Yan X, et al. Comparative proteome analysis of wheat flag leaves and developing grains under water deficit. Frontiers in Plant Sci. 2018; 9: 425.
  47. 47. Wang Y, Zhang X, Huang G, Feng F, Liu X, Guo R, et al. iTRAQ-based quantitative analysis of responsive proteins under PEG-induced drought stress in wheat leaves. Int. J. Mol. Sci. 2019; 20: 2621.
  48. 48. Coleman HD, Yan J, Mansfield SD. Sucrose synthase affects carbon partitioning to increase cellulose production and altered cell wall ultrastructure. Proc Natl Acad Sci USA. 2009; 106: 13118–23. pmid:19625620
  49. 49. Chaves MM. Effects of water deficits on carbon assimilation. J. Exp. Bot. 1991; 42: 1–16.
  50. 50. Lawlor DW, Cornic G. Photosynthetic carbon assimilation and associated metabolism in relation to water deficits in higher plants. Plant Cell Environ. 2002; 25: 275–94. pmid:11841670
  51. 51. Pinheiro C, Chaves MM. Photosynthesis and drought: can we make metabolic connections from available data? J. Exp. Bot. 2011; 63: 869–82.
  52. 52. Ghaffari MR, Mirzaei M, Ghabooli M, Khatabi B, Wu Y, Zabet-Moghaddam M, et al. Root endophytic fungus Piriformospora indica improves drought stress adaptation in barley by metabolic and proteomic reprogramming. Environ. Exp. Bot. 2019; 157: 197–210.
  53. 53. Lee ES, Kang CH, Park JH, Lee SY. Physiological significance of plant peroxiredoxins and the structure-related and multifunctional biochemistry of peroxiredoxin 1. Antioxid Redox Signal. 2018; 28: 625–639. pmid:29113450
  54. 54. Skirycz A, Inze D. More from less: Plant growth under limited water. Curr Opin Biotechnol. 2010; 21: 197–203. pmid:20363612
  55. 55. Olmos E, García De La Garma J, Gomez-Jimenez MC, Fernandez-Garcia N. Arabinogalactan proteins are involved in salt-adaptation and vesicle trafficking in Tobacco by-2 cell cultures. Front Plant Sci. 2017; 8: 1092. pmid:28676820
  56. 56. Nguema-Ona E, Vicre-Gibouin M, Cannesan MA, Driouich A. Arabinogalactan proteins in root-microbe interactions. Trends Plant Sci. 2013; 18: 440–449. pmid:23623239
  57. 57. Zang L, Zheng T, Chu Y, Ding C, Zhang W, Huang Q, et al. Genome-wide analysis of the fasciclin-like arabinogalactan protein gene family reveals differential expression patterns, localization, and salt stress response in Populus. Front. Plant Sci. 2015; 6: 1140. pmid:26779187
  58. 58. Zhang W, Qin C, Zhao J, Wang X. Phospholipase Da1-derived phosphatidic acid interacts with ABI1 phosphatase 2C and regulates abscisic acid signaling. Proc Natl Acad Sci USA. 2004; 101: 9508–9513. pmid:15197253
  59. 59. Hong Y, Zheng S, Wang X. Dual functions of phospholipase Dα1 in plant response to drought. Mol Plant. 2008; 1: 262–269. pmid:19825538
  60. 60. Bargmann BO, Laxalt AM, Riet B-ter, van Schooten B, Merquiol E, Testerink C, et al. Multiple PLDs required for high salinity and water deficit tolerance in plants. Plant Cell Physiol. 2009; 50: 78–89. pmid:19017627
  61. 61. Wang J, Ding B, Guo Y, Li M, Chen S, Huang G, et al. Overexpression of a wheat phospholipase D gene, TaPLDa, enhances tolerance to drought and osmotic stress in Arabidopsis thaliana. Planta. 2014; 240: 103–115. pmid:24705986
  62. 62. Maarouf HE, Zuily-Fodil Y, Gareil M, d'Arcy-Lameta A, Pham-Thi AT. Enzymatic activity and gene expression under water stress of phospholipase D in two cultivars of Vigna unguiculata L. Walp. differing in drought tolerance. Plant Mol Biol. 1999; 39: 1257–1265. pmid:10380811
  63. 63. Guo BZ, Xu G, Cao YG, Holbrook CC, Lynch RE. Identification and characterization of phospholipase D and its association with drought susceptibilities in peanut (Arachis hypogaea). Planta. 2006; 223: 512–520. pmid:16200410
  64. 64. Choudhary MK, Basu D, Datta A, Chakraborty N, Chakraborty S. Dehydration responsive nuclear proteome of rice (Oryza sativa L.) illustrates protein network, novel regulators of cellular adaptation, and evolutionary perspective. Mol Cell Proteomics. 2009; 8: 1579–98. pmid:19321431
  65. 65. Guo QF, Zhang J, Gao Q, Xing SC, Li F, Wang W. Drought tolerance through overexpression of monoubiquitin in transgenic tobacco. J. Plant Physiol. 2008; 165: 1745–1755. pmid:18280007
  66. 66. Miazek A, Zagdańska B. Involvement of exopeptidases in dehydration tolerance of spring wheat seedlings. Biol Plantarum. 2008; 52: 687–694.
  67. 67. Akinci S, Sösel DM. 2012. Plant water-stress response mechanisms. In: Rahman IMM, editor. Water stress. Intech, Rijeka, Croatia; 2012. p. 15–42.
  68. 68. Garg BK, Kathju S, Burman U. Influence of water stress on water relations, photosynthetic parameters and nitrogen metabolism of moth bean genotypes. Biol Plantarum. 2001; 44: 289–292.
  69. 69. Xu ZZ, Zhou GS. Nitrogen metabolism and photosynthesis in Leymus chinensis in response to long-term soil drought. J. Plant Growth Regul. 2006; 25: 252–266.
  70. 70. Singh NB, Singh D, Singh A. Biological seed priming mitigates the effects of water stress in sunflower seedlings. Physiol Mol Biol Plant. 2015; 21: 207–214.
  71. 71. Sinclair TR, Pinter PJ Jr, Kimball BA, Adamsen FJ, LaMorte RL, Wall GW, et al. Leaf nitrogen concentration of wheat subjected to elevated [CO2] and either water or N deficits. Agr Ecosyst Environ. 2000; 79: 53–60.
  72. 72. Fulda S, Mikkat S, Stegmann H, Horn R. Physiology and proteomics of drought stress acclimation in sunflower (Helianthus annuus L.). Plant Biol. 2010; 13: 632–642.
  73. 73. Shanker AK, Maheswari M, Yadav SK, Desai S, Bhanu D, Attal NB. Drought stress responses in crops. Funct Integr Genomics. 2014; 14: 11–22. pmid:24408129
  74. 74. Umezawa T, Nakashima K, Miyakawa T, Kuromori T, Tanokura M, Shinozaki K, et al. Molecular basis of the core regulatory network in ABA responses: sensing, signaling and transport. Plant Cell Physiol. 2010; 51: 1821–1839. pmid:20980270
  75. 75. Zhang X, Ervin EH, Evanylo GK, Haering KC. Impact of biosolids on hormone metabolism in drought-stressed tall fescue. Crop Sci. 2009; 49: 1893–1901.
  76. 76. Qin X, Zeevaart JAD. The 9-cis-epoxycarotenoid cleavage reaction is the key regulatory step of abscisic acid biosynthesis in water-stressed bean. Proc Natl Acad Sci USA. 1999; 96: 15354–15361. pmid:10611388
  77. 77. Iuchi S, Kobayashi M, Taji T, Naramoto M, Seki M, Kato T, et al. Regulation of drought tolerance by gene manipulation of 9-cis-epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis. Plant J. 2001; 27: 325–333. pmid:11532178
  78. 78. Asada K. The water-water cycle in chloroplasts: Scavenging of active oxygens and dissipation of excess photons. Annu Rev Plant Physiol Plant Mol Biol. 1999; 50: 601–639. pmid:15012221
  79. 79. Hossain MA, Nakano Y, Asada K. Monodehydroascorbate reductase in spinach-chloroplasts and its participation in regeneration of ascorbate for scavenging hydrogen-peroxide. Plant Cell Physiol. 1984; 25: 385–395.
  80. 80. Noctor G, Foyer CH. Ascorbate and glutathione: Keeping active oxygen under control. Annu Rev Plant Physiol Plant Mol Biol. 1998; 49: 249–279. pmid:15012235
  81. 81. Foyer CH, Noctor G. Ascorbate and glutathione: the heart of the redox hub. Plant Physiol. 2011; 155: 2–18. pmid:21205630
  82. 82. Vieira Dos Santos C, Rey P. Plant thioredoxins are key actors in the oxidative stress response. Trends Plant Sci. 2006; 11: 329e334.
  83. 83. Lee JC, Timashe SN. The stabilization of proteins by sucrose. J Biol Chem. 1981; 256: 7193–7201. pmid:7251592
  84. 84. Chmielewska K, Rodziewicz P, Swarcewicz B, Sawikowska A, Krajewski P, Marczak Ł, et al. Analysis of drought-induced proteomic and metabolomic changes in Barley (Hordeum vulgare L.) leaves and roots unravels some aspects of biochemical mechanisms involved in drought tolerance. Front. Plant Sci. 2016; 7: 1108. pmid:27512399
  85. 85. Swindell WR, Huebner M, Weber AP. Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics. 2007; 8: 125. pmid:17519032
  86. 86. Kaldenhoff R, Ribas-Carbo M, Sans JF, Lovisolo C, Heckwolf M, Uehlein N. Aquaporins and plant water balance. Plant Cell Environ. 2008; 31: 658e666.
  87. 87. Foyer CH, Noctor G. Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005; 7: 1866e1875.
  88. 88. Vierling E. The roles of heat-shock proteins in plants. Annu Rev Plant Biol. 1991; 42: 579–620.
  89. 89. Grigorova B, Vaseva I, Demirevska K, Feller U. Combined drought and heat stress in wheat: changes in some heat shock proteins. Biol Plantarum. 2011; 55: 105–111.
  90. 90. Alexandersson E, Fraysse L, Sjovall-Larsen S, Gustavsson S, Fellert M, Karlsson M, et al. Whole gene family expression and drought stress regulation of aquaporins. Plant Mol. Biol. 2005; 59: 469–484. pmid:16235111
  91. 91. Aroca R, Ferrante A, Vernieri P, Chrispeels MJ. Drought, abscisic acid and transpiration rate effects on the regulation of PIP aquaporin gene expression and abundance in Phaseolus vulgaris plants. Annal Bot. 2006; 98: 1301–1310.
  92. 92. Yue C, Cao H, Wang L, Zhou Y, Hao X, Zeng J, et al. Molecular cloning and expression analysis of tea plant aquaporin (AQP) gene family. Plant Physiol Biochem. 2014; 83: 65–76. pmid:25093260
  93. 93. Pandey V, Ansari MW, Tula S, Yadav S, Sahoo RK, Shukla N, et al. Dose-dependent response of Trichoderma harzianum in improving drought tolerance in rice genotypes. Planta. 2016; 243: 1251–1264. pmid:26898554