
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39256456
71847
10.1038/s41598-024-71847-9
Article
Comparative physiological and transcriptomic analyses reveal genotype specific response to drought stress in Siberian wildrye (Elymus sibiricus)
An Yongping 1
Wang Qian 2
Cui Yannong 3
Liu Xin 4
Wang Ping 1
Zhou Yue 1
Kang Peng 5
Chen Youjun 1
Wang Zhiwei 2
Zhou Qingping 1
Wang Pei wangpei@swun.edu.cn

1
1 https://ror.org/04gaexw88 grid.412723.1 0000 0004 0604 889X Sichuan Zoige Alpine Wetland Ecosystem National Observation and Research Station, College of Grassland Resources, Institute of Qinghai-Tibetan Plateau, Southwest Minzu University, Chengdu, 610225 China
2 https://ror.org/00ev3nz67 grid.464326.1 0000 0004 1798 9927 Guizhou Institute of Prataculture, Guizhou Academy of Agricultural Sciences, Guiyang, 550006 China
3 https://ror.org/0051rme32 grid.144022.1 0000 0004 1760 4150 College of Grassland Agriculture, Northwest A&F University, Yangling, 712100 China
4 grid.35155.37 0000 0004 1790 4137 National Key Laboratory of Crop Genetic Improvement, Ministry of Agriculture Key Laboratory of Crop Ecophysiology and Farming System in the Middle Reaches of the Yangtze River, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, 430070 China
5 https://ror.org/05xjevr11 grid.464238.f 0000 0000 9488 1187 College of Biological Science and Engineering, North Minzu University, Yinchuan, 750021 China
10 9 2024
10 9 2024
2024
14 2106011 4 2024
31 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Siberian wildrye (Elymus sibiricus) is a xero-mesophytic forage grass with high nutritional quality and stress tolerance. Among its numerous germplasm resources, some possess superior drought resistance. In this study, we firstly investigated the physiological differences between the leaves of drought-tolerant (DT) and drought-sensitive (DS) genotypes under different field water contents (FWC) in soil culture. The results showed that, under drought stress, DT maintained a lower leaf water potential for water absorption, sustained higher photosynthetic efficiency, and reduced oxidative damage in leaves by efficiently maintaining the ascorbic acid-glutathione (ASA-GSH) cycle to scavenge reactive oxygen species (ROS) compared to DS. Secondly, using RNA sequencing (RNA-seq), we analyzed the gene expression profiles of DT and DS leaves under osmotic stress of hydroponics induced by PEG-6000. Through differential analysis, we identified 1226 candidate unigenes, from which we subsequently screened out 115/212 differentially expressed genes (DEGs) that were more quickly induced/reduced in DT than in DS under osmotic stress. Among them, Unigene0005863 (EsSnRK2), Unigene0053902 (EsLRK10) and Unigene0031985 (EsCIPK5) may be involved in stomatal closure induced by abscisic acid (ABA) signaling pathway. Unigene0047636 (EsCER1) may positively regulates the synthesis of very-long-chain (VLC) alkanes in cuticular wax biosynthesis, influencing plant responses to abiotic stresses. Finally, the contents of wax and cutin were measured by GC–MS under osmotic stress of hydroponics induced by PEG-6000. Corresponding to RNA-seq, contents of wax monomers, especially alkanes and alcohols, showed significant induction by osmotic stress in DT but not in DS. It is suggested that limiting stomatal and cuticle transpiration under drought stress to maintain higher photosynthetic efficiency and water use efficiency (WUE) is one of the critical mechanisms that confer stronger drought resistance to DT. This study provides some insights into the molecular mechanisms underlying drought tolerance in E. sibiricus. The identified genes may provide a foundation for the selection and breeding of drought-tolerant crops.

Keywords

Siberian wildrye
Osmotic resistance
Differentially expressed genes (DEGs)
Stress tolerance
ROS
ABA
Cuticular wax
Subject terms

Plant sciences
Plant molecular biology
Plant physiology
Plant stress responses
National Natural Science Foundation of China31802122&32360424 Wang Pei Sichuan Science and Technology Program2024NSFSC0310 Wang Pei Southwest Minzu University Double World-Class ProjectCX2023013 Wang Pei issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Drought is the most common abiotic stress that inhibits plant growth and development and poses a serious threat to crop production sustainability1,2. Every year, numerous cases of reduced crop productivity due to droughts result in substantial economic losses, presenting a significant challenge for agricultural production3,4. Furthermore, with global climate change and the continued rise in greenhouse gas concentrations, droughts have become more frequent, severe, and prolonged than in recent decades5–7. Enhancing plant resilience is imperative for agricultural production and environmental sustainability8. Therefore, gaining a comprehensive understanding of the molecular mechanisms governing plant responses and adaptations to drought stress and breeding superior varieties with high drought resistance are promising approaches for ensuring global food security and promoting ecological restoration worldwide9–11.

In general, virtually every aspect of plant physiology, biochemistry, and cellular metabolism is affected by drought stress12–16. Drought results in hyperosmotic stress in plant cells, inhibiting root system water absorption, leading to abnormal growth, an imbalance in water relations, and decreased water-use efficiency (WUE)4,8,17–19. In addition, drought disrupts ionic homeostasis, osmotic balance, the cell membrane system, and photosynthesis13,14,20,21. Drought also leads to a reduction in stomatal conductance, leaf CO2 assimilation rate, and enzyme activity13,15,20. Moreover, the dry environment induced the production of reactive oxygen species (ROS) in biochemical reactions, including superoxide anion (O2·-), hydrogen peroxide (H2O2), hydroxyl radical (HO·), and singlet oxygen (1O2)13,22. ROS production and consumption are strictly controlled by the antioxidant systems15,18,20. Under drought stress, the activities of many oxidative protective enzymes, such as superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX) and glutathione reductase (GR), are altered3. Ascorbate/Ascorbic acid (ASA) and glutathione (GSH) are crucial for maintaining ROS levels within plant cells, as they prevent oxidative damage by scavenging these free radicals13,23. Dehydroascorbate (DHA) is the oxidized form of ASA and can be converted back into ASA through reduction processes in plant cells, thereby maintaining the levels of ASA in the plant24. During the process of plant adaptation to drought stress, the ASA/DHA system protects cellular structures and maintains normal physiological functions through its antioxidant capabilities25. However, as the drought progressively worsens, ROS production overwhelms the scavenging function of the antioxidant system, leading to severe cell damage and death26.

H2O2 has been identified as a key signaling molecule that regulates stomatal closure in plants27,28. Stomatal closure is the first defense mechanism by which plants respond to drought by preventing transpiration loss, but at the cost of reducing photosynthetic rates13. This reduction in photosynthesis can have significant consequences on plant growth under drought conditions. In addition to the ROS pathway, plants regulate stomatal closure through various pathways, including the abscisic acid (ABA) and Ca2+ signaling pathways. ABA-mediated stomatal closure is a crucial drought resistance mechanism in plants29. Moreover, the cuticle, which is composed of insoluble cutin and waxes, is another factor that limits transpiration under drought stress and serves as an important barrier for plants against various biotic and abiotic stresses30,31. Cutin plays a vital role in the prevention of moisture loss32. Leaf epidermal wax increases as soil moisture decreases, and the wax layer reduces stomatal conductance, transpiration, and photosynthesis, thereby facilitating plant adaptation to drought conditions33. Nonetheless, a recent study in Populus revealed no alteration in the overall amount of leaf waxes, but rather an elevation in their transpiration rate due to modifications in their composition34.

In addition to ROS scavenging and stomatal closure, plants undergo significant alterations in physiological and biochemical metabolism, and regulate root growth and morphology as adaptive responses to drought environments. These adaptations are intricately regulated by multiple signal networks, including signal perception, signal transduction, transcription, translation, and post-translational modifications2,35–37. It has been identified that several key components within these signaling networks play pivotal roles in regulating drought responses in various plant species, such as Arabidopsis, rice, and maize38–40. These components include receptors, phosphatases, kinases, and a wide array of transcription factors (TFs). These TFs encompass those involved in phytohormone signaling, such as PYR/PYL/RCAR (pyrabactin resistance 1/PYR1-like/regulatory component of ABA receptors), PP2C (clade A type 2C protein phosphatases), SnRK2s/SAPKs (SNF1/AMPK-related protein kinases), calcium-dependent protein kinases (CDPKs), as well as abscisic acid responsive elements binding factor (ABF), abscisic acid responsive elements (ABRE), myeloblastosis (MYB), WRKY, DREB, NAC, and BHLH TFs8,36,41–45. However, our knowledge of the genome-wide molecular mechanisms underlying drought tolerance in plants, especially in gramineous crops and forage grasses, remains limited.

Many native grass species grown in arid and semi-arid areas possess special stress tolerance mechanisms and abundant gene resources. Research on these grasses can potentially significantly contribute to the genetic improvement of crop stress resistance, thereby benefiting agriculture and animal husbandry. Siberian wildrye (Elymus sibiricus), a typical perennial meso-xerophytic forage grass widely distributed at high altitudes worldwide, offers good forage quality and environmental adaptability. It is extensively used in natural grasslands, cultivated pastures, and the ecological restoration of degraded grasslands46–49. Furthermore, E. sibiricus exhibits remarkable resistance to barren conditions, drought, cold, salinity, and diseases48,50. However, despite its impressive resilience, the mechanisms underlying its drought tolerance remain poorly understood.

Drought stress leads to water deficit in plants, while osmotic stress leads to water absorption difficulties, Therefore, many studies use osmotic stress to simulate drought stress. Polyethylene Glycol 6000 (PEG-6000), a nonionic osmotic agent, effectively simulates drought conditions by reducing water potential and is a valuable tool for studying plant responses to drought51. Drought, resulting from insufficient precipitation or soil moisture under natural conditions, leads to reduced water availability for plants. Similarly, PEG-6000-induced osmotic stress decreases soil water potential, limiting plant water absorption52. While these two phenomena differ in mechanism, physiological impact, adaptation strategies, and research applications, both focus on how plants adapt to water scarcity53. In our previous study, we screened drought-tolerant (DT, I-1-1-46) and drought-sensitive (DS, I-1-1-2) genotypes of E. sibiricus, and revealed that the Casparian strip and suberin lamellae in the root endodermis played a pivotal role in the difference in drought resistance between the DT and DS genotypes54. Here, we performed a comparative analysis of the physiological, biochemical and transcriptomic changes in the leaves of DT and DS genotypes in response to drought stress to clarify two major issues: (1) What are the differences in the physiological responses to drought stress in the leaves of different genotypes? (2) What is the molecular basis underlying the distinct drought tolerances of the two genotypes? The results of this study are expected to shed new light on the understanding of plant drought resistance, offering potential pathways for the development of more drought-tolerant forage grass and crop cultivars through molecular breeding methods.

Results

Effects of drought on leaf water status in DT and DS

Under the control condition (75% field water content [FWC]), there was no significant differences in leaf water potential (Ψw), osmotic potential (Ψs) and turgor pressure (Ψt) between DT and DS (Table 1). With increasing drought intensity, both DT and DS leaves exhibited significant decreases in Ψw and Ψs. Under mild (60% FWC) and moderate (45% FWC) drought stress, although there was no significant difference in Ψw between the two genotypes, Ψs of DT was significantly lower than that of DS (Table 1). Under severe drought stress (30% FWC), Ψw and Ψs of DT leaves were significantly lower than those of DS leaves (Table 1). With the exacerbation of drought stress, Ψt in both genotypes initially increased and then decreased, and it was significantly lower in DT than in DS under severe drought stress (Table 1). These results indicate that E. sibiricus can absorb water from drought soils by reducing the osmotic potential and, consequently, water potential, with DT exhibiting a more advantageous capacity in this regard.Table 1 Effects of drought on leaf water status of two E. sibiricus genotypes.

Genotype	Moisture gradient
Field water content (FWC) %	Water potential Ψw (Mpa)	Osmotic potential Ψs (Mpa)	Turgor pressure Ψt (Mpa)	
DT	75	− 0.82 ± 0.04 e	− 1.29 ± 0.01 f	0.46 ± 0.03 d	
60	− 1.19 ± 0.04 d	− 1.99 ± 0.02 d	0.80 ± 0.03 b	
45	− 2.06 ± 0.04 c	− 2.98 ± 0.01 b	0.91 ± 0.03 a	
30	− 2.83 ± 0.06 a	− 3.18 ± 0.06 a	0.35 ± 0.13 e	
DS	75	− 0.82 ± 0.04 e	− 1.27 ± 0.01 f	0.45 ± 0.03 de	
60	− 1.19 ± 0.04 d	− 1.85 ± 0.02 e	0.66 ± 0.06 c	
45	− 1.98 ± 0.02 c	− 2.89 ± 0.03 c	0.92 ± 0.04 a	
30	− 2.53 ± 0.08 b	− 2.99 ± 0.03 b	0.46 ± 0.05 d	
Values are presented as the means ± standard deviation (SD) (n = 6) and bars indicate SD; columns with different letters indicate significant differences at P < 0.05 (Tukey’s test). DT drought-tolerant genotype, DS drought-sensitive genotype, FWC field water content.

Effects of drought stress on photosynthesis of two E. sibiricus genotypes

Under the control conditions, there were no significant differences in net photosynthetic rate (Pn), stomatal conductance (Gs), and water use efficiency (WUE) between DT and DS, whereas the transpiration rate (Tr), intercellular carbon dioxide concentration (Ci), and limiting value of stomata (Ls) in DT were significantly higher than those in DS (Fig. 1). With decreasing soil water content, the Ci values of both DT and DS gradually increased, while the other five photosynthetic parameters showed a decreasing trend. The increase in Ci was as pronounced as the reduction in Tr and Ls, with Pn exhibiting the largest magnitude of change (Fig. 1). Overall, under drought stress, DT consistently exhibited higher values for all six photosynthetic indices compared to DS. The most pronounced differences were observed under severe drought stress, where DT had a 46, 30, 13, 17, and 19% higher Pn, Gs, Ci, WUE, and Ls than DS. Notably, at this point, the Tr of the DT leaves was was not significantly higher than that of the DS leaves (Fig. 1).Fig. 1 Responses of leaf photosynthetic characteristics in two E. sibiricus genotypes to drought stress. (A) Net photosynthetic rate (Pn); (B) transpiration rate (Tr); (C) stomatal conductance (Gs); (D) intercellular carbon dioxide concentration (Ci); (E) water use efficiency (WUE); (F) limiting value of stomata (Ls). The data are presented as the means of three replicates (± SD) based on Tukey's post hoc test. Means denoted with the different letters indicate significant differences at P < 0.05. The horizontal coordinates indicate different levels of field water content (FWC). DT drought-tolerant genotype, DS drought-sensitive genotype.

Differences in leaf ROS accumulation and ascorbate–glutathione cycling in drought responses

The superoxide anion radical (O2·-) content did not significantly differ between DT and DS in either the control or mild to moderate drought stress (≥ 45% FWC). However, under severe drought stress, the O2·- content in DS was notably higher (38%) than that in DT (Fig. 2A). Hydrogen peroxide (H2O2) levels in both DT and DS leaves increased with decreasing soil water content. Under mild drought stress (60% FWC), there was no significant difference in H2O2 levels between DT and DS; however, under other drought conditions, H2O2 levels in DT were significantly lower than those in DS (Fig. 2B). GSH and ASA are important antioxidants involved in ROS scavenging. There was little difference in GSH content between DT and DS across the 30–75% drought gradient (Fig. 2C). However, ASA content increased with increasing drought intensity and was significantly higher in DT than in DS under mild to severe drought stress (Fig. 2D).Fig. 2 Responses of leaf ascorbic acid-glutathione cycle in two E. sibiricus genotypes to drought stress. (A) Superoxide anion radical (O2·-); (B) hydrogen peroxide (H2O2); (C) glutathione (GSH); (D) ascorbic acid (ASA); (E) mono-dehydroascorbic acid reductase (MDHAR); (F) dehydroascorbate reductase (DHAR); (G) ascorbate peroxidase (APX); (H) ascorbic acid/dehydroascorbic acid (ASA/DHA). The data are presented as the means of three replicates (± SD) based on Tukey's post hoc test. Means denoted with the different letters indicate significant differences at P < 0.05. The horizontal coordinates indicate different levels of field water content (FWC). DT drought-tolerant genotype, DS drought-sensitive genotype.

Therefore, we analyzed ASA-related indices, including mono-dehydroascorbic acid reductase (MDHAR), dehydroascorbate reductase (DHAR), APX, and ASA/DHA. The results showed that the MDHAR, DHAR contents, and ASA/DHA continuously increased with increasing drought intensity (Fig. 2E,F,H), whereas the APX content initially decreased drastically and then slightly increased (Fig. 2G). In the control treatment (75% FWC), the contents of MDHAR, DHAR, and ASA/DHA were not significantly different between DT and DS, whereas the APX content was significantly higher in DT than in DS (Fig. 2E–H). The MDHAR content was significantly higher in DT than in DS under severe drought stress, but there was no significant difference between them in the remaining treatments (Fig. 2E). Under mild drought stress (60% FWC), the DHAR content was not significantly different between DT and DS, while under moderate (45% FWC) and severe (30% FWC) drought stress, the DHAR content in DT was significantly higher than in DS (Fig. 2F). The APX content and ASA/DHA of DT were significantly higher than those of DS under all drought treatments (Fig. 2G,H). These findings highlight the importance of ASA metabolism in the response of DT to drought stress.

Transcriptome analysis of DT and DS in response to osmotic stress

Libraries were constructed and sequenced on an Illumina HiSeq™ 4000 deep sequencing platform, and detailed RNA sequencing results are shown in Supplementary Table S1. Clean reads were assembled to obtain 136,699 unigenes; detailed information on the evaluation indicators can be found in Supplementary Table S2. These unigenes were annotated in four protein databases: 80,610 in Nr, 67,301 in KEGG, 31,763 in COG, and 39,932 in SwissProt (Supplementary Table S3 and Supplementary Fig. S1). The RNA-seq data were credible and reproducible using PCA and Pearson’s correlation coefficient analyses (Supplementary Figs. S2, and S3).

Volcano plots identified significantly upregulated or downregulated genes during osmotic stress at each time point (Fig. 3A–D). In DT, there were 3002 differentially expressed genes (DEGs) (1664 upregulated and 1338 downregulated) and 2464 DEGs (986 upregulated and 1478 downregulated) were observed at 3 and 24 h, respectively, compared with the control (Fig. 3E). In DS, there were 1555 DEGs (481 upregulated and 1074 down-regulated) at 3 h and 1691 DEGs (873 upregulated and 818 down-regulated) at 24 h compared with the control (Fig. 3E). Venn diagrams showed in Supplementary Fig. S4 illustrated the number of DEGs that were common and unique to two genotypes at the same treatment time (Supplementary Fig. S4). At 3 h, 2572 and 1125 DEGs were uniquely expressed in DT and DS, respectively, when compared to the control, 430 DEGs were common expressed in both DT and DS (Supplementary Fig. S4A). At 24 h, 1735 and 962 DEGs were uniquely expressed in DT and DS, respectively, when compared to the control, 729 DEGs were common expressed in both DT and DS (Supplementary Fig. S4B). Of which 197 DEGs were co-expressed at two treatment times (Supplementary Fig. S4C). Another common form of analysis prioritizes the number of common and uniquely expressed DEGs in the same genotype at different treatment times (Fig. 3F–G), 2234 and 1696 DEGs were unique in DT at 3 and 24 h, respectively, compared to the control (Fig. 3F), whereas 900 and 1036 DEGs were uniquely expressed in DS at 3 and 24 h, respectively, compared with the control (Fig. 3G). In addition, 768 and 655 DEGs were common at both time points for both DT and DS, respectively (Fig. 3F,G), indicating that a higher number of unigenes were regulated in DT than in DS. Further analysis revealed 197 DEGs were co-expressed in DS and DT (Fig. 3H).Fig. 3 Analysis of differentially expressed genes in E. sibiricus genotypes in response to osmotic stress. (A–D) Volcano plots of all expressed genes in drought-tolerant (DT) and drought-sensitive (DS) genotypes at different time points compared to the control (CK), with log2FC values drawn against -log10 (FDR) adjusted P-values. Orange and blue dots represent upregulated and downregulated differentially expressed genes (DEGs) based on |log2(foldchange)|≥ 2 and FDR ≤ 0.001; green dots represent genes that are not differentially expressed; (E) number of upregulated and downregulated DEGs identified under different osmotic-stress treatments compared to the control DT and DS. Orange indicates upregulated DEGs, while blue indicates downregulated DEGs; (F) common and unique DEGs under 3 and 24 h osmotic-stress treatments in DT; (G) common and unique DEGs under 3 and 24 h osmotic-stress treatments in DS; (H) common and unique DEGs under both 3 and 24 h osmotic-stress treatments in DT and DS.

GO enrichment analysis of the upregulated and downregulated genes revealed key GO term associated with osmotic stress. GO enrichment for DEGs was higher in DS than in DT, the number of genes annotated to GO was higher in DS (fg_num = 1380) than in DT (fg_num = 978) (Supplementary Fig. S5 and Supplementary Table S4). DS and DT are mainly enriched in the two GO terms associated with hydrolase activity in molecular function, hydrolase activity, acting on glycosyl bonds (GO:0016798) and hydrolase activity, hydrolyzing O-glycosyl compounds (GO:0004553), respectively (Supplementary Fig. S5 and Supplementary Table S4). Using the KEGG database to investigate the key genes involved in the KEGG pathways, these DEGs were found to be significantly enriched in plant hormone signal transduction (ko04075), biosynthesis of secondary metabolites (ko01110), and the MAPK signaling pathway (ko04016) (Supplementary Fig. S6 and Supplementary Table S5).

Screening of differential candidate genes in DT and DS under osmotic stress

We performed hierarchical clustering for 1226 DEGs (log2FC value) in response to osmotic stress in the two E. sibiricus genotypes (Fig. 3H). The heat map identified two clustering modules for the DEGs. Cluster1 consists of 833 DEGs that exhibited high expression levels in DT and DS under normal conditions, but low expression levels under osmotic stress. In contrast, cluster2 includes 393 DEGs that showed the opposite expression trend compared to cluster1 (Fig. 4A). Further analysis of cluster1 identified five notable trend blocks linked to profiles 3, 10, 2, 6, and 0, with 212, 167, 153, 134, and 95 unigenes, respectively. Additionally, cluster2 exhibited four significant trend blocks attributed to profiles 9, 16, 17, and 13, comprising 118, 115, 55, and 54 unigenes, respectively (Fig. 4B,C).Fig. 4 Clustering and trend analysis of differentially expressed genes (DEGs) in DT and DS genotypes under osmotic stress. (A) Hierarchical clustering analysis indicates DEG expression patterns in DT and DS after 0, 3, and 24 h of osmotic stress treatment; (B,C) trend analysis shows that the lines reflect expression patterns for each DEG at osmotic stress time points. Figures (B) and (C) respectively demonstrate the expression trends of DEGs in Cluster 1 and Cluster 2. The horizontal axes of Figures (B) and (C) are the same as that of Figure (A). Black lines in each trend blocks indicate mean changes in DEG expression. The different profiles represent the numbers of the different trend blocks, and the numbers within the trend blocks indicate the number of enriched genes. Colored trend blocks indicate significant enrichment (p < 0.05), while uncolored trend blocks do not show significant enrichment. Each inflection point represents a set of sample data.

To efficiently screen candidate genes associated with DT osmotic tolerance, we focused on profile3 within cluster1 and profile 16 within cluster2. In profile3 of cluster1, 212 DEGs exhibit high expression levels in both DT and DS genotypes in the absence of stress treatment. Expression levels decreased significantly in both DT and DS by osmotic stress, with a faster reduction in DT (from 3 to 24 h) compared to DS (24 h) under treatment. At 24 h, both exhibited similar low expression levels (Figs. 4B and 5A,B). Heatmap analysis showed that the expression of 212 DEGs in profile 3 was lower in the DT genotype compared to the DS genotype under osmotic stress (Fig. 5B). In profile16 of cluster1, 115 DEGs exhibit lower expression levels in both DT and DS genotypes without treatment. However, DT was highly induced compared to DS in the early stages of the 3-h osmotic stress treatment, and both were highly expressed in the 24-h stress treatment (Figs. 4C and 5C). Heatmap analysis showed that the expression of 115 DEGs in profile 16 was higher in the DT genotype than in the DS genotype under osmotic stress (Fig. 5D). Immediately after that, we adopted the following strategy to further screen the candidate unigenes: by analyzing the expression differences of DEGs between two genotypes under osmotic stress treatments and combining the information from related literature, we screened the unigenes that were specifically low-expressed in DT under osmotic stress in profile3 of Cluster1, and also screened the unigenes that were specifically high-expressed in DT under osmotic stress in profile16 of Cluster2. Following this strategy, eight candidate unigenes that were highly expressed specifically in DT under osmotic stress were screened from 327 DEGs: Unigene0005863 (EsSnRK2), Unigene0047636 (EsCER1), Unigene0005889 (EsGOLS1), Unigene0066379 (EsCAD8D), Unigene0005868 (EsSAHH), Unigene0031985 (EsCIPK5), Unigene0041137 (EsAZF2), and Unigene0053902 (EsLRK10) (Fig. 5E). Interestingly, all 8 candidate unigenes are from profile 16 of cluster 2. The expression of these eight candidate genes in DT was significantly upregulated under osmotic stress, whereas they remained unchanged in DS in response to osmotic stress. These candidate unigenes were enriched in cutin, suberin, and wax biosynthesis, galactose metabolism, cysteine and methionine metabolism, MAPK signaling pathway, plant hormone signal transduction, phenylpropanoid biosynthesis, biosynthesis of secondary metabolites and metabolic pathways (Fig. 5F). Among these, Unigene0005863 (EsSnRK2) encodes a key protein kinase in the MAPK and plant hormone (ABA) signal transduction pathways (Supplementary Figs. S8 and S9). In addition, our study found that unigenes associated with ROS scavenging in the ASA-GSH system were differentially expressed in DT and DS (Supplementary Fig. S10).Fig. 5 The trend analysis plots, expression heatmap, and enrichment bubble plots of DEGs in DT and DS genotypes at different time points. (A) Trend analysis graph for DT and DS in profile 3 trend block from cluster1; (B) heatmap of DT and DS expression in the profile 3 trend block; (C) Trend analysis graph for DT and DS in profile 16 trend block from cluster2; (D) heatmap of DT and DS expression in the profile 16 trend block; (E) heatmap expression of candidate genes; (F) enrichment bubble map of candidate genes.

Effects of osmotic stress on leaf cuticular wax and cutin monomer contents

Among the DEGs presented in Fig. 5E, we noted that Unigene0047636, a homologous gene of CER1 in Arabidopsis through sequence alignment which may encode an enzyme involved in forming alkane monomers during wax biosynthesis, was significantly induced by osmotic stress in DT, but not in DS (Supplementary Fig. S7). Therefore, we analyzed the cuticular wax and cutin monomer contents. Results showed that there was little difference in the wax and cutin monomer contents between the two genotypes without osmotic stress (Fig. 6A,B). However, osmotic stress (20% PEG-6000) significantly increased the contents of alkanes and alcohols aldehydes in the wax monomers of DT, without increasing any monomers in DS, resulting in the two major wax monomers being higher pronouncedly in DT than in DS under osmotic stress (Fig. 6A). However, only one DEG, Unigene0047636 (CER1), was detected in the wax biosynthesis pathway (Supplementary Fig. S7), suggesting that CER1 may be a key factor in wax biosynthesis in the DT response to osmotic stress. In addition, although osmotic stress slightly increased the contents of some cutin monomers in DT, such as alcohols, ω-OH fatty acids, and polyhydroxy fatty acids (PFAs), which were not induced by osmotic stress in DS, transcriptome sequencing did not detect any DEGs between the two genotypes in the cutin biosynthesis pathway (Fig. 6B, Supplementary Fig. S7).Fig. 6 Difference in leaf cuticle between DT and DS genotypes under osmotic stress. (A) Wax monomers content; (B) Cutin monomers content. Wax and cutin data are presented as the means ± SD based on Tukey's post hoc test. Means denoted with the different letters indicate significant differences at P < 0.05. DT drought-tolerant genotype, DS drought-sensitive genotype.

Discussion

The drought tolerance of plants is strongly related to their ability to use’ water17. The difference in the water potential gradient between the inside and outside of plant cells is the key to water transport across the membrane55. Under drought stress, plants can reduce cellular water potential through osmoregulation to maintain the water potential difference between the root system and the soil environment, and to absorb water from dry soils with low water potential56. Under drought stress, DT had a lower leaf Ψw and Ψs than DS (Table 1), which indicated the higher drought resistance of DT. Furthermore, our results showed that the Ψt of DT was significantly lower than that of DS under severe drought stress (Table 1). Cell turgor pressure plays a pivotal role in governing the stomatal aperture, with reduced turgor pressure predominantly leading to stomatal closure57. However, the Gs of DT was significantly higher than that of DS, corresponding to a higher Pn of DT than that of DS (Fig. 1A,C), indicating a higher photosynthetic capacity of DT. In addition, Ci and WUE of DT were higher than those of DS under different drought stress treatments (Fig. 1D,E). The factors causing reduced photosynthesis under drought stress can be divided into two categories: stomatal limitation and non-stomatal limitation; the former is due to a limitation of Gs, causing Ci to fail to meet photosynthetic demand, and the latter is due to reduced chloroplast and photosynthetic enzyme activities58,59. In the present study, the Gs of both genotypes gradually decreased with increasing drought intensity, whereas the Ci gradually increased and Ls decreased (Fig. 1C,D). This result indicated that the stomatal limitation of E. sibiricus gradually decreased with increasing drought intensity, and the factors leading to the inhibition of photosynthesis gradually changed to non-stomatal limitation.

In response to osmotic stress, plants control water balance through stomatal movement60. Plants employ multiple pathways to close stomata under drought stress, with the ABA signaling pathway being one of the most crucial mechanisms29. In addition to the ABA signaling pathway, other signaling pathways, such as ROS, can also regulate stomatal opening and closure during drought stress61. These mechanisms work synergistically to reduce water transpiration and enhance drought tolerance62. Our transcriptome study screened eight DEGs involved in the response of E. sibiricus to osmotic stress, which were significantly upregulated in DT, but not in DS, under osmotic stress (Fig. 5E). Among these eight genes, Unigene0005863 (EsSnRK2), Unigene0053902 (EsLRK10) and Unigene0031985 (EsCIPK5) may be involved in stomatal closure induced by the ABA signaling pathway. SnRK2 is a key protein kinase involved in the ABA signal transduction pathway, which regulates stomatal closure to prevent leaf water loss in Arabidopsis63. The wheat SnRK2s family member TaSnRK2.1 is Involved in stomatal closure, osmotic regulator biosynthesis, cellular ROS scavenging, homeostasis, and other stress-related biological processes during drought stress and ABA induction64. Unigene0053902 (EsLRK10) is homologous to LPK10L1.2 in Arabidopsis. Arabidopsis mutants lacking LPK10L1.2 exhibit ABA-insensitive and drought-sensitive phenotypes65. CdtCIPK5 from Cynodon dactylon (L.), the homologous gene of Unigene0031985 (EsCIPK5), is induced by multiple abiotic stresses such as ABA, dehydration, low temperature, and salt stress66. The high expression of these eight candidate genes in DT under osmotic stress may be the key to its drought resistance, but how they regulate drought resistance at the molecular level in E. sibiricus DT genotypes requires further study.

Stomata and cuticle are the two pathways that regulate transpiration water loss67,68. Plant cuticular waxes are the key components that reduce cuticle water loss through the cuticle69, acting as a crucial barrier to protect plants against environmental stresses30,70. These waxes, located in the leaf epidermis, are believed to contribute to the reduction of stomatal conductance, transpiration rates, and photosynthesis, thereby enhancing drought resistance of plants32,33. However, in the current research, we found that the wax content of DT was significantly higher than that of DS. Its stomatal conductance, transpiration rates, and photosynthetic activity were also significantly higher than those of DS. This may be because, under drought stress, the lower wax content of DS necessitates a drastic reduction in transpiration water loss by decreasing Gs, while DT limits cuticle transpiration due to its higher wax content, which in turn helps to maintain a higher Gs. This promotes CO2 assimilation and balances plant growth under drought stress. In addition, in our previous studies, DT has stronger apoplastic barrieres in the root endodermis, which can more effectively limit the return of water to the medium under drought stress. This is also one of the important strategies for maintaining photosynthesis54.

Waxes are primarily composed of very-long-chain (VLC) aliphatic lipids, such as alkanes, aldehydes, primary and secondary alcohols, ketones, and esters30,70. This protective function is particularly important in xerophytes. Shellakkutti et al.71 discovered that osmotic stress in both cultivated and wild barley led to an increase in wax crystals accumulation, total wax amounts, and wax gene expression, while cuticular conductance remained unaffected. Similarly, a recent study in Populus revealed that alterations in the compositions of wax, despite without altering the total amounts, led to elevated transpiration rates34. Alkanes, the predominant constituents of plant cuticle waxes72, are thought to play a crucial role in conferring drought tolerance73. Grünhofer et al.34 pinpointed a rise in the relative proportion of esters in the Populus × canescens cer6_1,2 mutant as the primary factor to the heightened transpiration rate, in fact, alkanes in the mutant decreased significantly. Notably, our results showed that the content of alkane, alcohol, and aldehyde monomer in DT was significantly higher than that in DS (Fig. 6A), suggesting that DT may enhance drought tolerance by accumulating more wax alkanes.

The wax biosynthetic pathway involves multiple steps and pathways, among which the modification of aldehydes, alkanes, secondary alcohols, and ketones by very-long-chain fatty acids (VLCFAs) through the alkane formation pathway is an extremely critical part74. This process is regulated by multiple genes that encode enzymes and transporters necessary for the deposition of cuticle constituents75. Among these genes, CER1 plays a pivotal role in the drought stress-induced biosynthesis of cuticular wax76,77. In a study on alkane synthesis in wheat epidermal waxes, researchers found that TaCER1-1A, a TaCER1 homolog, increased the alkane contents in both transgenic rice and Arabidopsis, leading to an improved epidermal barrier and enhanced drought tolerance in the leaf epidermis of transgenic lines78. Additionally, TaCER1-6A is a crucial gene in wheat, significantly involed in alkane biosynthesis. Knockout mutants of TaCER1-6A exhibited reduced alkane contents in their leaves, whereas overexpression of TaCER1-6A increased alkane contents, thereby improving drought tolerance in wheat73. Furthermore, in tomato and Poa pratensis, the SlCER1-1 and PpCER1-2, respectively, positively regulate the synthesis of wax alkanes, reducing cuticle permeability and enhancing drought tolerance79,80. Wen et al.81 found that the downregulation of two TaCER1 genes led to decreased long-chain alkanes synthesis in low-wax wheat mutants. Drought stress caused a significant increase in very-long-chain alkanes in the wax components of Brachypodium distachyon. BdCER1-8, homologous gene of CER1, is highly expressed in leaves82. In Arabidopsis, CER1 encodes an alkane-forming enzyme that catalyses the decarbonylation of aliphatic aldehydes to alkanes83. Our results showed that the expression of Unigene0047636, homologous to CER1, was upregulated in DT but remained unchanged in DS under drought stress (Fig. 5E). Harb et al.84 found that in barley, HvCER1 was specifically induced under drought stress in the drought-tolerant genotype, but not in the drought-sensitive genotype, which is consistent with our results. EsCER1 positively regulates the conversion of long-chain aldehyde monomers to long-chain alkane monomers during wax biosynthesis, which consequently leads to a significant increase in the contents of alkanes monomers in wax within DT (Fig. 6A). Furthermore, in the wax biosynthesis pathway, long-chain secondary alcohols, as downstream products of long-chain alkanes (Supplementary Fig. S7), increased correspondingly with the increased contents of long-chain alkanes (Fig. 6A), although no DEGs regulating the synthesis of alcohols from alkanes were detected (Supplementary Fig. S7). Therefore, we speculate that EsCER1 is involved in the alkanes synthesis of the wax biosynthesis pathway and plays a certain regulatory role in the synthesis of alcohols, which may be an important reason for the difference in drought resistance among the different genotypes of E. sibiricus.

The reactive oxide species (ROS) are essential signaling molecules in various physiological processes, playing a crucial role in regulating plant growth and development85. However, excessive ROS accumulation can be toxic to plants and can lead to oxidative damage86,87. To counteract ROS toxicity, plants possess a complex detoxification system, that includes a non-enzymatic ASA-GSH pathway88. The ASA-GSH pathway serves as a vital component in reducing oxidative damage by scavenging ROS through the detoxification mechanisms of ASA, GSH, and their associated enzymes89,90. The ASA-GSH cycle is comprised of four key enzymes: ascorbate peroxidase (APX), mono-dehydroascorbate reductase (MDHAR), dehydroascorbate reductase (DHAR), and glutathione reductase (GR)90. In the ASA-GSH cycle, ASA provides electrons for the APX-catalyzed reaction to complete the effective scavenging of H2O291. In this study, we found that under drought stress, the contents of ASA and APX in DT were significantly higher than those in DS (Fig. 2D,G), but its H2O2 content was significantly lower than that in DS(Fig. 2B), indicating that DT has a stronger H2O2 scavenging capacity. In the ASA-GSH cycle, GSH provides electrons for DHAR and promotes the reduction of DHA to regenerate ASA92. In this study, we found that the DHAR content in DT was significantly higher than that in DS under drought stress (Fig. 2F), suggesting stronger ASA regeneration in DT. Meanwhile, as GSH is oxidized to Glutathione Disulfide (GSSG) during ASA regeneration, the degree of GSH oxidation in DT might be greater than in DS, which could also lead to no significant difference in GSH content between DT and DS (Fig. 2C). The above studies suggest that the greater drought tolerance of DT compared to DS may be due to its ability to maintain lower ROS levels through higher ASA, DHAR, and APX concentrations in the ASA-GSH cycle.

APX1 and APX2 are cytoplasmic enzymes that play crucial roles in ROS scavenging93,94. Previous studies have demonstrated their significance in mitigating the ROS levels in plants. For instance, in a study on ROS accumulation in Arabidopsis, it was found that the overexpression of APX1 and APX2 resulted in a significant decrease in ROS levels95. Research on Stevia rebaudiana has demonstrated that the plant's response to abiotic stress involves the upregulation of key genes encoding ROS-scavenging enzymes such as glutathione S-transferase (GST) and APX, contributing to the modulation of its stress response96. Under osmotic stress, RNA-seq results revealed that genes relevant to the ASA-GSH cycle (APX, GR, DHAR, and GST) were differentially expressed in DT and DS (Supplementary Fig. S10). In addition, the overall expression of Unigene0106502 (EsAPX1), Unigene0005132 (EsAPX2), Unigene0013814 (EsGR1) and Unigene0023285 (EsGST4) were higher in DT than in DS (Supplementary Fig. S10). These findings may suggest a potential role for these unigenes in the regulation of ROS accumulation in DT and DS.

Materials and methods

Plant materials and treatments

Two genotypes (DT and DS) of E. sibiricus were used in this study54. Plant materials used to determine physiological phenotypes were grown in soil, materials for cutin, wax and transcriptome sequencing were cultured by hydroponics.

For soil culture: Seeds were vernalized for 3 days in the dark at 4 °C, then germinated on wet filter paper. After 2 days, the seedlings were transferred to pots (diameter, 10 cm; depth, 8.5 cm; 9 seedlings/pot) containing peat soil, and incubated at 23 and 19 °C during the day and night, respectively, with 16 h of light, and 40–60% relative humidity. A compound fertilizer with a 1:1:1 N:P:K ratio was applied every 2 days. After 10 days, when the plants reached the third-leaf stage, they were divided into four experimental groups: (1) control: watered with distilled water until water content was maintained at 75% of the field water content (FWC) by weight; (2) mild drought treatment: stopped watering until water content of the medium decreased to 60% FWC, and distilled water was used to maintain this content; (3) moderate drought treatment: maintained at 45% FWC; (4) severe drought treatment: maintained at 30% FWC. The FWC was calculated by assessing the moisture level of the peat soil using a cutting ring technique. Correlation indices were measured after 12 days of treatment, and six biological replicates were used for each treatment (three plants/replicate).

In hydroponics, seed germination was the same as that observed in the soil culture experiment. Two days after germination, the seedlings were transferred to hydroponic boxes filled with half-strength Hoagland solution (H353, Phyto TechLabs, Lenexa, KS, USA), which was aerated using an air pump. The environmental conditions in the climatic chamber were consistent with those of the soil culture. Subsequently, 28-day-old seedlings were treated as follows: (1) Seedlings were transferred to half-strength Hoagland solution containing 0 (control) or 20% PEG-6000 (BioFroxx, Einhausen, Germany) for 3 days. Cutin and wax of the leaf tissue were then measured after harvest. Each treatment consisted of six biological replicates. (2) The seedlings were transferred to half-strength Hoagland solution containing 20% PEG-6000 for treatment. Fresh leaves were collected at 0, 3, and 24 h after treatment and immediately frozen in liquid nitrogen, then stored at – 80 °C for transcriptome sequencing. Each treatment had three biological replicates.

Measurement of leaf water potential, osmotic potential and turgor pressure

The leaf water potential (Ψw) was measured using a water potential meter (PSYPRO; DIANJIANG TECH, USA). For the osmotic potential (Ψs), leaf samples (0.5 g) were washed with deionized water, dried with filter paper, and immediately frozen in liquid nitrogen. Subsequently, the samples were placed in a syringe the sap was squeezed out. The sap was then centrifuged, and 50 μL of the supernatant was used to measure the ic value with a freezing point osmometer (OSMOMAT 3000, Gonotec GmbH, Germany). Ψs was calculated according to the formula:Ψs=-icRT.

where R is the gas constant 0.008314 and T is the measured thermodynamic temperature (298.8 K). The turgor pressure (Ψt) is calculated according to the formula:Ψt=Ψw-Ψs.

Measurement of photosynthetic parameters

Photosynthetic parameters were measured using a portable photosynthesis system instrument (LI-6800, LI-COR Biosciences, USA). To ensure that the plants received sufficient light before measurement. The measurements were conducted between 8:00 and 11:30 a.m. The photosynthetic system parameters were set as described below: leaf chamber area of 2 cm2, photosynthetic effective radiation of 1000 lx, and temperature of 25 °C. The following photosynthetic indicators were measured: the net photosynthetic rate (Pn), stomatal conductance (Gs), transpiration rate (Tr), and intercellular carbon dioxide concentration (Ci). Based on the above measured indices, the water use efficiency was calculated according to the formula: WUE = Pn/Gs97 and the limiting value of stomata (Ls) was calculated using the formula: Ls = (Ca − Ci)/Ca98, where Ca is the atmospheric carbon dioxide concentration.

Determination of ROS and antioxidant contents

Plant leaves (0.1 g) were used to measure the content of ROS (O2·- and H2O2) in the control and drought treatment groups using O2·- and H2O2 kits, respectively. Ascorbic acid (ASA), glutathione (GSH) and dehydroascorbic acid (DHA) levels were determined using ASA, GSH and DHA kits, respectively. The activity levels of circulatory metabolism-related enzymes, including dehydroascorbate reductase (DHAR), monodehydroascorbic acid reductase (MDHAR) and ascorbate peroxidase (APX), were measured using DHAR, MDHAR and APX kits, respectively. All the kits were obtained from Suzhou Comin Biotechnology Co., Ltd. (Suzhou, China). The measurement of these indices were performed according to the manufacturer’s instructions.

Determination of cutin and wax contents

Leaves were used to determine the cutin and wax content. Cutin monomers were extracted via tissue delipidation, ester depolymerization, and silylation, as described by Jenkin and Molina99. Waxes were determined as follows: leaf tissues were extracted by placing leaves in 10 mL of chromatographically pure chloroform for 30 s; 25 µL C24 alkane (1 μg/µL) was added as an internal standard, then the solutions were blow-dried with nitrogen, silylated with 100 µL of N, O Bis-trimethylsilyl-trifluoroacetamid (BSTFA) and 100 µL of pyridine for 30 min at 100 °C, and blow-dried with nitrogen again, and finally dissolved in 100 µL of hexane.

The cutin and wax monomer contents were determined using a gas chromatograph-mass spectrometer (GC–MS; 8890-7000D, Agilent Technologies, USA) fitted with an HP-5MS capillary column (length, 30 m; i.d., 0.25 mm; film thickness, 0.25 mm). For cutin, the injector was set at 250 °C, and the injected split ratio was 1:10 and helium was the carrier gas at a constant flow rate of 1.5 mL min−1. The oven was initially set at 50 °C for 2 min, increased by 40 °C min−1 increments to 200 °C, maintained for 5 min, then increased by 5 °C min−1 increments to 310 °C, and maintained for 20 min. The temperature of the MS detector was 325 °C, and the MS was set to scan mode > 40–600 amu (electron impact ionization). For wax, the injector was set at 250 °C, the oven was initially set at 80 °C for 1 min, increased by 15 °C min−1 increments to 260 °C, maintained for 10 min, then increased by 5 °C min−1 increments to 320 °C, and maintained for 24 min. The temperature of the MS detector was 320 °C, and the MS was set to scan mode > 40–600 amu. Six biological replicates were assessed in each experiment.

RNA extraction, library construction and sequencing

Total RNA was extracted using a TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and verified using RNase-free agarose gel electrophoresis. Conventional kits were used to remove rRNA and enrich mRNA. The enriched mRNA was then cut into short fragments using a fragmentation buffer and reverse transcribed into cDNA using random primers. The cDNA fragments were purified using the QiaQuick PCR extraction kit (Qiagen, Venlo, Netherlands), end-repaired, poly(A) added, and ligated to Illumina sequencing adapters. The ligation products were size-selected by agarose gel electrophoresis, PCR-amplified, and sequenced using an Illumina HiSeqTM 4000 by Gene Denovo Biotechnology Co. (Guangzhou, China). Each treatment had three biological replicates.

RNA-seq data filtering of clean reads, denovo assembly and basic annotation of unigenes

Raw reads obtained from sequencing machines, including raw reads containing adapters or low-quality bases, affect subsequent assembly and analysis. Thus, to obtain high-quality clean reads, the reads were further filtered using Fastp software (version 0.18.0)100. The parameters were as follows: (1) reads containing adapters wer removed; (2) reads containing more than 10% of unknown nucleotides (N) were reomved; (3) low-quality reads containing more than 50% of low quality (Q-value ≤ 20) bases were removed. Transcriptome denovo assembly was conducted using the short reads assembly program, Trinity101. To annotate the unigenes, we used the BLASTx program (http://www.ncbi.nlm.nih.gov/BLAST/) with an E-value threshold of 1e-5 in the NCBI non-redundant protein (Nr) (http://www.ncbi.nlm.nih.gov), Swiss-Prot protein (http://www.expasy.ch/sprot), Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg), and COG/KOG databases (http://www.ncbi.nlm.nih.gov/COG). The protein with the highest sequence similarity to the given unigene was identified and the protein function annotation information of the unigene was obtained.

Differential expression and target genes analysis

The Omicsmart online data analysis platform (https://www.omicsmart.com/home.html#/) was used to calculate the differentially expressed genes (DEGs) between DT and DS genotypes under the control condition and 20% PEG-6000 treatment. Significant DEGs between the osmotic stress and control conditions were identified at |log2(foldchange) |> 2 and false discovery rate (FDR) adjusted P < 0.001. Subsequently, GO and KEGG pathway enrichment analysis was performed on the Omicsmart online data analysis platform (https://www.omicsmart.com/home.html#/) for differential basic genes in DT and DS under PEG-induced osmotic stress for 0, 3, and 24 h. Furthermore, a cluster heat map and trend analysis of the co-expressed genes were performed. Candidate genes with similar expression patterns were grouped together using the Omicsmart heatmap analysis tool (https://www.omicsmart.com/home.html#/). The clustering results were then divided into different trend blocks. Analysis of candidate genes trends using the Omicsmart trend analysis tool (https://www.omicsmart.com/home.html#/). Significant trend blocks were identified when |log2(foldchange)|> 2 and P < 0.05. The different profiles represent the numbers of the different trend blocks. Finally, pathway enrichment analysis was performed for candidate genes under osmotic stress.

Statistical analysis

Statistical analyses were performed using Microsoft Excel 2016. Analysis of physiological and phenotypic characterizations data was performed using IBM SPSS Statistics (version 26.0). Physiological and phenotypic characterizations values are shown as the mean ± standard deviation (SD) of biological replicates. Tukey's post hoc test was performed with a P-value < 0.05, considered significantly different. Origin 2021 was used for data visualization. The Omicsmart online data analysis platform (https://www.omicsmart.com/home.html#/) was used to visualize transcriptome sequencing data.

Conclusions

This study combined physiological and transcriptomic analyses to reveal differences in drought tolerance between two genotypes of E. sibiricus, DT and DS. The physiological results confirmed that DT had greater drought tolerance than DS. Transcriptome analysis identified multiple drought tolerance genes. By combining the physiological and biochemical results, we found that these genes may participate in several drought tolerance pathways in E. sibiricus, including ABA signaling pathway-induced stomatal closure, leaf cuticular wax deposition, ASA metabolism, and ROS scavenging. Unigene0005863 (EsSnRK2), Unigene0053902 (EsLRK10) and Unigene0031985 (EsCIPK5) may be involved in stomatal closure induced by abscisic acid (ABA) signaling pathway. Unigene0047636 (EsCER1) may positively regulates the synthesis of very-long-chain (VLC) alkanes in cuticular wax biosynthesis, influencing plant responses to abiotic stresses. This study provides some ideas for research on the mechanisms of drought tolerance in meso-xerophytic forage grasses and gives some clues for genetic improvement of drought tolerance in forages and crops.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71847-9.

Acknowledgements

We are very grateful to Sichuan Zoige Alpine Wetland Ecosystem National Observation and Research Station and College of Grassland Resources, Institute of Qinghai-Tibetan Plateau, Southwest Minzu University for providing us with experimental equipment and venue. Thanks to all authors for their hard work on this manuscript. We are grateful thank Guangzhou Gene Denovo Biotechnology Co., Ltd for assisting in sequencing and bioinformatics analysis. Thank you to Editage for polishing our manuscript! This study was supported by grants from the National Natural Science Foundation of China (31802122&32360424), Sichuan Science and Technology Program (2024NSFSC0310), and Southwest Minzu University Double World-Class Project (CX2023013).

Author contributions

Conceptualization, Y.A. and Pe.W.; methodology, X.L., Pi.W. and Pe.W.; software, Y.A. and Y.Z.; investigation, X.L., Pi.W., Y.Z. and Y.A.; resources, Q.Z., Y.Ch. and Pe.W.; data curation, Y.A., Y.Z. and Pe.W.; writing—original draft preparation, Y.A. and Pe.W.; writing—review and editing, Y.A., Q.W., Y.C., P.K., Z.W. and Pe.W. All authors have read and agreed to the published version of the manuscript.

Data availability

Sequence data that support the findings of this study have been deposited in the National Center for Biotechnology Information with the primary accession code PRJNA1101572.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Moore JP Vicre-Gibouin M Farrant JM Driouich A Adaptations of higher plant cell walls to water loss: Drought vs desiccation Physiol. Plant. 2008 134 237 245 10.1111/j.1399-3054.2008.01134.x 18494857
Moore, J. P., Vicre-Gibouin, M., Farrant, J. M. & Driouich, A. Adaptations of higher plant cell walls to water loss: Drought vs desiccation. Physiol. Plant. 134, 237–245 (2008).18494857 10.1111/j.1399-3054.2008.01134.x
2. Zhang H Sun X Dai M Improving crop drought resistance with plant growth regulators and rhizobacteria: Mechanisms, applications, and perspectives Plant Commun. 2022 3 100228 10.1016/j.xplc.2021.100228 35059626
Zhang, H., Sun, X. & Dai, M. Improving crop drought resistance with plant growth regulators and rhizobacteria: Mechanisms, applications, and perspectives. Plant Commun. 3, 100228 (2022).35059626 10.1016/j.xplc.2021.100228
3. Fang Y Xiong L General mechanisms of drought response and their application in drought resistance improvement in plants Cell Mol. Life Sci. 2015 72 673 689 10.1007/s00018-014-1767-0 25336153
Fang, Y. & Xiong, L. General mechanisms of drought response and their application in drought resistance improvement in plants. Cell Mol. Life Sci. 72, 673–689 (2015).25336153 10.1007/s00018-014-1767-0
4. Razi K Muneer S Drought stress-induced physiological mechanisms, signaling pathways and molecular response of chloroplasts in common vegetable crops Crit. Rev. Biotechnol. 2021 41 669 691 10.1080/07388551.2021.1874280 33525946
Razi, K. & Muneer, S. Drought stress-induced physiological mechanisms, signaling pathways and molecular response of chloroplasts in common vegetable crops. Crit. Rev. Biotechnol. 41, 669–691 (2021).33525946 10.1080/07388551.2021.1874280
5. Intergovernmental Panel on Climate, C. Climate Change 2021—The Physical Science Basis. (2023).
6. Ault TR On the essentials of drought in a changing climate Science. 2020 368 256 260 10.1126/science.aaz5492 32299944
Ault, T. R. On the essentials of drought in a changing climate. Science. 368, 256–260 (2020).32299944 10.1126/science.aaz5492
7. Muller LM Bahn M Drought legacies and ecosystem responses to subsequent drought Glob. Change Biol. 2022 28 5086 5103 10.1111/gcb.16270
Muller, L. M. & Bahn, M. Drought legacies and ecosystem responses to subsequent drought. Glob. Change Biol. 28, 5086–5103 (2022).10.1111/gcb.16270
8. Zhu JK Abiotic stress signaling and responses in plants Cell. 2016 167 313 324 10.1016/j.cell.2016.08.029 27716505
Zhu, J. K. Abiotic stress signaling and responses in plants. Cell. 167, 313–324 (2016).27716505 10.1016/j.cell.2016.08.029
9. Chaves MM Maroco JP Pereira JS Understanding plant responses to drought—From genes to the whole plant Funct. Plant Biol. 2003 30 239 264 10.1071/FP02076 32689007
Chaves, M. M., Maroco, J. P. & Pereira, J. S. Understanding plant responses to drought—From genes to the whole plant. Funct. Plant Biol. 30, 239–264 (2003).32689007 10.1071/FP02076
10. Qi J Reactive oxygen species signaling and stomatal movement in plant responses to drought stress and pathogen attack J. Integr. Plant Biol. 2018 60 805 826 10.1111/jipb.12654 29660240
Qi, J. et al. Reactive oxygen species signaling and stomatal movement in plant responses to drought stress and pathogen attack. J. Integr. Plant Biol. 60, 805–826 (2018).29660240 10.1111/jipb.12654
11. Geng MF Genome-wide investigation on transcriptional responses to drought stress in wild and cultivated rice Environ. Exp. Bot. 2021 189 104555 10.1016/j.envexpbot.2021.104555
Geng, M. F. et al. Genome-wide investigation on transcriptional responses to drought stress in wild and cultivated rice. Environ. Exp. Bot. 189, 104555 (2021).10.1016/j.envexpbot.2021.104555
12. Zhu JK Salt and drought stress signal transduction in plants Annu. Rev. Plant Biol. 2002 53 247 273 10.1146/annurev.arplant.53.091401.143329 12221975
Zhu, J. K. Salt and drought stress signal transduction in plants. Annu. Rev. Plant Biol. 53, 247–273 (2002).12221975 10.1146/annurev.arplant.53.091401.143329
13. Mahajan S Tuteja N Cold, salinity and drought stresses: An overview Arch. Biochem. Biophys. 2005 444 139 158 10.1016/j.abb.2005.10.018 16309626
Mahajan, S. & Tuteja, N. Cold, salinity and drought stresses: An overview. Arch. Biochem. Biophys. 444, 139–158 (2005).16309626 10.1016/j.abb.2005.10.018
14. Shinozaki K Yamaguchi-Shinozaki K Gene networks involved in drought stress response and tolerance J. Exp. Bot. 2007 58 221 227 10.1093/jxb/erl164 17075077
Shinozaki, K. & Yamaguchi-Shinozaki, K. Gene networks involved in drought stress response and tolerance. J. Exp. Bot. 58, 221–227 (2007).17075077 10.1093/jxb/erl164
15. Shao HB Chu LY Jaleel CA Zhao CX Water-deficit stress-induced anatomical changes in higher plants Cr Biol. 2008 331 215 225 10.1016/j.crvi.2008.01.002
Shao, H. B., Chu, L. Y., Jaleel, C. A. & Zhao, C. X. Water-deficit stress-induced anatomical changes in higher plants. Cr Biol. 331, 215–225 (2008).10.1016/j.crvi.2008.01.002
16. Pandey V Shukla A Acclimation and tolerance strategies of rice under drought stress Rice Sci. 2015 22 147 161 10.1016/j.rsci.2015.04.001
Pandey, V. & Shukla, A. Acclimation and tolerance strategies of rice under drought stress. Rice Sci. 22, 147–161 (2015).10.1016/j.rsci.2015.04.001
17. Farooq M Wahid A Kobayashi N Fujita D Basra SMA Plant drought stress: Effects, mechanisms and management Agron. Sustain. Dev. 2009 29 185 212 10.1051/agro:2008021
Farooq, M., Wahid, A., Kobayashi, N., Fujita, D. & Basra, S. M. A. Plant drought stress: Effects, mechanisms and management. Agron. Sustain. Dev. 29, 185–212 (2009).10.1051/agro:2008021
18. Zhang H Zhu J Gong Z Zhu JK Abiotic stress responses in plants Nat. Rev. Genet. 2022 23 104 119 10.1038/s41576-021-00413-0 34561623
Zhang, H., Zhu, J., Gong, Z. & Zhu, J. K. Abiotic stress responses in plants. Nat. Rev. Genet. 23, 104–119 (2022).34561623 10.1038/s41576-021-00413-0
19. Verma S Negi NP Pareek S Mudgal G Kumar D Auxin response factors in plant adaptation to drought and salinity stress Physiol. Plant. 2022 174 e13714 10.1111/ppl.13714 35560231
Verma, S., Negi, N. P., Pareek, S., Mudgal, G. & Kumar, D. Auxin response factors in plant adaptation to drought and salinity stress. Physiol. Plant. 174, e13714 (2022).35560231 10.1111/ppl.13714
20. Reddy AR Chaitanya KV Vivekanandan M Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants J. Plant Physiol. 2004 161 1189 1202 10.1016/j.jplph.2004.01.013 15602811
Reddy, A. R., Chaitanya, K. V. & Vivekanandan, M. Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants. J. Plant Physiol. 161, 1189–1202 (2004).15602811 10.1016/j.jplph.2004.01.013
21. Nezhadahmadi A Prodhan ZH Faruq G Drought tolerance in wheat Sci. World J. 2013 2013 610721 10.1155/2013/610721
Nezhadahmadi, A., Prodhan, Z. H. & Faruq, G. Drought tolerance in wheat. Sci. World J. 2013, 610721 (2013).10.1155/2013/610721
22. Waszczak C Carmody M Kangasjärvi J Reactive oxygen species in plant signaling Annu. Rev. Plant Biol. 2018 69 209 236 10.1146/annurev-arplant-042817-040322 29489394
Waszczak, C., Carmody, M. & Kangasjärvi, J. Reactive oxygen species in plant signaling. Annu. Rev. Plant Biol. 69, 209–236 (2018).29489394 10.1146/annurev-arplant-042817-040322
23. Gill SS Tuteja N Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants Plant Physiol. Biochem. 2010 48 909 930 10.1016/j.plaphy.2010.08.016 20870416
Gill, S. S. & Tuteja, N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930 (2010).20870416 10.1016/j.plaphy.2010.08.016
24. Saitoh Y Umezaki T Yonekura N Nakawa A Resveratrol potentiates intracellular ascorbic acid enrichment through dehydroascorbic acid transport and/or its intracellular reduction in HaCaT cells Mol. Cell Biochem. 2020 467 57 64 10.1007/s11010-020-03700-2 32080778
Saitoh, Y., Umezaki, T., Yonekura, N. & Nakawa, A. Resveratrol potentiates intracellular ascorbic acid enrichment through dehydroascorbic acid transport and/or its intracellular reduction in HaCaT cells. Mol. Cell Biochem. 467, 57–64 (2020).32080778 10.1007/s11010-020-03700-2
25. Anjum NA Metal/metalloid stress tolerance in plants: Role of ascorbate, its redox couple, and associated enzymes Protoplasma. 2014 251 1265 1283 10.1007/s00709-014-0636-x 24682425
Anjum, N. A. et al. Metal/metalloid stress tolerance in plants: Role of ascorbate, its redox couple, and associated enzymes. Protoplasma. 251, 1265–1283 (2014).24682425 10.1007/s00709-014-0636-x
26. Cruz de Carvalho MH Drought stress and reactive oxygen species: Production, scavenging and signaling Plant Signal Behav. 2008 3 156 165 10.4161/psb.3.3.5536 19513210
Cruz de Carvalho, M. H. Drought stress and reactive oxygen species: Production, scavenging and signaling. Plant Signal Behav. 3, 156–165 (2008).19513210 10.4161/psb.3.3.5536
27. Desikan R ABA, hydrogen peroxide and nitric oxide signalling in stomatal guard cells J. Exp. Bot. 2004 55 205 212 10.1093/jxb/erh033 14673026
Desikan, R. et al. ABA, hydrogen peroxide and nitric oxide signalling in stomatal guard cells. J. Exp. Bot. 55, 205–212 (2004).14673026 10.1093/jxb/erh033
28. Karpinska B Wingsle G Karpinski S Antagonistic effects of hydrogen peroxide and glutathione on acclimation to excess excitation energy in Arabidopsis Iubmb Life. 2000 50 21 26 10.1080/15216540050176548 11087116
Karpinska, B., Wingsle, G. & Karpinski, S. Antagonistic effects of hydrogen peroxide and glutathione on acclimation to excess excitation energy in Arabidopsis. Iubmb Life. 50, 21–26 (2000).11087116 10.1080/15216540050176548
29. Hsu PK Dubeaux G Takahashi Y Schroeder JI Signaling mechanisms in abscisic acid-mediated stomatal closure Plant J. 2021 105 307 321 10.1111/tpj.15067 33145840
Hsu, P. K., Dubeaux, G., Takahashi, Y. & Schroeder, J. I. Signaling mechanisms in abscisic acid-mediated stomatal closure. Plant J. 105, 307–321 (2021).33145840 10.1111/tpj.15067
30. Bi H The impact of drought on wheat leaf cuticle properties BMC Plant Biol. 2017 17 1 13 10.1186/s12870-017-1033-3 28049439
Bi, H. et al. The impact of drought on wheat leaf cuticle properties. BMC Plant Biol. 17, 1–13 (2017).28049439 10.1186/s12870-017-1033-3
31. Tapia G Morales-Quintana L Parra C Berbel A Alcorta M Study of nsLTPs in Lotus japonicus genome reveal a specific epidermal cell member (LjLTP10) regulated by drought stress in aerial organs with a putative role in cutin formation Plant Mol. Biol. 2013 82 485 501 10.1007/s11103-013-0080-x 23733601
Tapia, G., Morales-Quintana, L., Parra, C., Berbel, A. & Alcorta, M. Study of nsLTPs in Lotus japonicus genome reveal a specific epidermal cell member (LjLTP10) regulated by drought stress in aerial organs with a putative role in cutin formation. Plant Mol. Biol. 82, 485–501 (2013).23733601 10.1007/s11103-013-0080-x
32. Chen G A functional cutin matrix is required for plant protection against water loss Plant Signal Behav. 2011 6 1297 1299 10.4161/psb.6.9.17507 22019635
Chen, G. et al. A functional cutin matrix is required for plant protection against water loss. Plant Signal Behav. 6, 1297–1299 (2011).22019635 10.4161/psb.6.9.17507
33. Mohammadian MA Watling JR Hill RS The impact of epicuticular wax on gas-exchange and photoinhibition in Leucadendron lanigerum (Proteaceae) Acta Oecol. 2007 31 93 101 10.1016/j.actao.2006.10.005
Mohammadian, M. A., Watling, J. R. & Hill, R. S. The impact of epicuticular wax on gas-exchange and photoinhibition in Leucadendron lanigerum (Proteaceae). Acta Oecol. 31, 93–101 (2007).10.1016/j.actao.2006.10.005
34. Grünhofer P Changes in wax composition but not amount enhance cuticular transpiration Plant Cell Environ. 2023 47 91 105 10.1111/pce.14719 37718770
Grünhofer, P. et al. Changes in wax composition but not amount enhance cuticular transpiration. Plant Cell Environ. 47, 91–105 (2023).37718770 10.1111/pce.14719
35. Joshi R Transcription factors and plants response to drought stress: Current understanding and future directions Front Plant Sci. 2016 7 1029 10.3389/fpls.2016.01029 27471513
Joshi, R. et al. Transcription factors and plants response to drought stress: Current understanding and future directions. Front Plant Sci. 7, 1029 (2016).27471513 10.3389/fpls.2016.01029
36. Zenda T Key Maize drought-responsive genes and pathways revealed by comparative transcriptome and physiological analyses of contrasting inbred lines Int. J. Mol. Sci. 2019 20 1268 10.3390/ijms20061268 30871211
Zenda, T. et al. Key Maize drought-responsive genes and pathways revealed by comparative transcriptome and physiological analyses of contrasting inbred lines. Int. J. Mol. Sci. 20, 1268 (2019).30871211 10.3390/ijms20061268
37. Wang Y Transcriptomic and physiological responses of contrasting maize genotypes to drought stress Front. Plant Sci. 2022 13 928897 10.3389/fpls.2022.928897 35991451
Wang, Y. et al. Transcriptomic and physiological responses of contrasting maize genotypes to drought stress. Front. Plant Sci. 13, 928897 (2022).35991451 10.3389/fpls.2022.928897
38. Wilkins O Bräutigam K Campbell MM Time of day shapes Arabidopsis drought transcriptomes Plant J. 2010 63 715 727 10.1111/j.1365-313X.2010.04274.x 20553421
Wilkins, O., Bräutigam, K. & Campbell, M. M. Time of day shapes Arabidopsis drought transcriptomes. Plant J. 63, 715–727 (2010).20553421 10.1111/j.1365-313X.2010.04274.x
39. Quan R Overexpression of an ERF transcription factor TSRF1 improves rice drought tolerance Plant Biotechnol. J. 2010 8 476 488 10.1111/j.1467-7652.2009.00492.x 20233336
Quan, R. et al. Overexpression of an ERF transcription factor TSRF1 improves rice drought tolerance. Plant Biotechnol. J. 8, 476–488 (2010).20233336 10.1111/j.1467-7652.2009.00492.x
40. Cao L Transcriptional regulatory networks in response to drought stress and rewatering in maize (Zea mays L.) Mol. Genet. Genom. 2021 296 1203 1219 10.1007/s00438-021-01820-y
Cao, L. et al. Transcriptional regulatory networks in response to drought stress and rewatering in maize (Zea mays L.). Mol. Genet. Genom. 296, 1203–1219 (2021).10.1007/s00438-021-01820-y
41. Sakuma Y Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression Plant Cell. 2006 18 1292 1309 10.1105/tpc.105.035881 16617101
Sakuma, Y. et al. Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell. 18, 1292–1309 (2006).16617101 10.1105/tpc.105.035881
42. Fujita Y AREB1 is a transcription activator of novel ABRE-dependent ABA signaling that enhances drought stress tolerance in Arabidopsis Plant Cell. 2005 17 3470 3488 10.1105/tpc.105.035659 16284313
Fujita, Y. et al. AREB1 is a transcription activator of novel ABRE-dependent ABA signaling that enhances drought stress tolerance in Arabidopsis. Plant Cell. 17, 3470–3488 (2005).16284313 10.1105/tpc.105.035659
43. Liu S Global transcriptome and weighted gene co-expression network analyses of growth-stage-specific drought stress responses in maize Front. Genet. 2021 12 645443 10.3389/fgene.2021.645443 33574835
Liu, S. et al. Global transcriptome and weighted gene co-expression network analyses of growth-stage-specific drought stress responses in maize. Front. Genet. 12, 645443 (2021).33574835 10.3389/fgene.2021.645443
44. Seki M Umezawa T Urano K Shinozaki K Regulatory metabolic networks in drought stress responses Curr. Opin. Plant Biol. 2007 10 296 302 10.1016/j.pbi.2007.04.014 17468040
Seki, M., Umezawa, T., Urano, K. & Shinozaki, K. Regulatory metabolic networks in drought stress responses. Curr. Opin. Plant Biol. 10, 296–302 (2007).17468040 10.1016/j.pbi.2007.04.014
45. Umezawa T Yoshida R Maruyama K Yamaguchi-Shinozaki K Shinozaki K SRK2C, a SNF1-related protein kinase 2, improves drought tolerance by controlling stress-responsive gene expression in Arabidopsis thaliana Proc. Natl. Acad. Sci. U. S. A. 2004 101 17306 17311 10.1073/pnas.0407758101 15561775
Umezawa, T., Yoshida, R., Maruyama, K., Yamaguchi-Shinozaki, K. & Shinozaki, K. SRK2C, a SNF1-related protein kinase 2, improves drought tolerance by controlling stress-responsive gene expression in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U. S. A. 101, 17306–17311 (2004).15561775 10.1073/pnas.0407758101
46. Ma X Chen S Zhang X Bai S Zhang C Assessment of worldwide genetic diversity of Siberian wild rye (Elymus sibiricus L.) germplasm based on gliadin analysis Molecules. 2012 17 4424 4434 10.3390/molecules17044424 22499189
Ma, X., Chen, S., Zhang, X., Bai, S. & Zhang, C. Assessment of worldwide genetic diversity of Siberian wild rye (Elymus sibiricus L.) germplasm based on gliadin analysis. Molecules. 17, 4424–4434 (2012).22499189 10.3390/molecules17044424
47. Xie W Zhang J Zhao X Zhang Z Wang Y Transcriptome profiling of Elymus sibiricus, an important forage grass in Qinghai-Tibet plateau, reveals novel insights into candidate genes that potentially connected to seed shattering BMC Plant Biol. 2017 17 1 15 10.1186/s12870-017-1026-2 28049439
Xie, W., Zhang, J., Zhao, X., Zhang, Z. & Wang, Y. Transcriptome profiling of Elymus sibiricus, an important forage grass in Qinghai-Tibet plateau, reveals novel insights into candidate genes that potentially connected to seed shattering. BMC Plant Biol. 17, 1–15 (2017).28049439 10.1186/s12870-017-1026-2
48. De Y Shi F Gao F Mu H Yan W Siberian Wildrye (Elymus sibiricus L.) abscisic acid-insensitive 5 gene is involved in abscisic acid-dependent salt response Plants. 2021 10 1351 10.3390/plants10071351 34371554
De, Y., Shi, F., Gao, F., Mu, H. & Yan, W. Siberian Wildrye (Elymus sibiricus L.) abscisic acid-insensitive 5 gene is involved in abscisic acid-dependent salt response. Plants. 10, 1351 (2021).34371554 10.3390/plants10071351
49. Zhang Z High-altitude genetic selection and genome-wide association analysis of yield-related traits in Elymus sibiricus L. using SLAF sequencing Front. Plant Sci. 2022 13 874409 10.3389/fpls.2022.874409 35800604
Zhang, Z. et al. High-altitude genetic selection and genome-wide association analysis of yield-related traits in Elymus sibiricus L. using SLAF sequencing. Front. Plant Sci. 13, 874409 (2022).35800604 10.3389/fpls.2022.874409
50. Yu Q Comparative metabolomic studies of siberian wildrye (Elymus sibiricus L.): A new look at the mechanism of plant drought resistance Int. J. Mol. Sci. 2022 24 452 10.3390/ijms24010452 36613896
Yu, Q. et al. Comparative metabolomic studies of siberian wildrye (Elymus sibiricus L.): A new look at the mechanism of plant drought resistance. Int. J. Mol. Sci. 24, 452 (2022).36613896 10.3390/ijms24010452
51. Zhou Q SNP application improves drought tolerance in soybean Sci. Rep. 2023 13 10911 10.1038/s41598-023-38088-8 37407630
Zhou, Q. et al. SNP application improves drought tolerance in soybean. Sci. Rep. 13, 10911 (2023).37407630 10.1038/s41598-023-38088-8
52. Li C Wan Y Shang X Fang S Integration of transcriptomic and metabolomic analysis unveils the response mechanism of sugar metabolism in Cyclocarya paliurus seedlings subjected to PEG-induced drought stress Plant Physiol. Biochem. 2023 201 107856 10.1016/j.plaphy.2023.107856 37354727
Li, C., Wan, Y., Shang, X. & Fang, S. Integration of transcriptomic and metabolomic analysis unveils the response mechanism of sugar metabolism in Cyclocarya paliurus seedlings subjected to PEG-induced drought stress. Plant Physiol. Biochem. 201, 107856 (2023).37354727 10.1016/j.plaphy.2023.107856
53. Rajeswar S Narasimhan SJ Peg-induced drought stress in plants: A review Res. J. Pharm. Technol. 2021 14 6173 6178 10.52711/0974-360X.2021.01069
Rajeswar, S. & Narasimhan, S. J. Peg-induced drought stress in plants: A review. Res. J. Pharm. Technol. 14, 6173–6178 (2021).10.52711/0974-360X.2021.01069
54. Liu X Endodermal apoplastic barriers are linked to osmotic tolerance in meso-xerophytic grass Elymus sibiricus Front. Plant Sci. 2022 13 1007494 10.3389/fpls.2022.1007494 36212320
Liu, X. et al. Endodermal apoplastic barriers are linked to osmotic tolerance in meso-xerophytic grass Elymus sibiricus. Front. Plant Sci. 13, 1007494 (2022).36212320 10.3389/fpls.2022.1007494
55. Fricke W Water transport and energy Plant Cell Environ. 2016 40 977 994 10.1111/pce.12848 27756100
Fricke, W. Water transport and energy. Plant Cell Environ. 40, 977–994 (2016).27756100 10.1111/pce.12848
56. Martinez-Vilalta J Garcia-Forner N Water potential regulation, stomatal behaviour and hydraulic transport under drought: Deconstructing the iso/anisohydric concept Plant Cell Environ. 2017 40 962 976 10.1111/pce.12846 27739594
Martinez-Vilalta, J. & Garcia-Forner, N. Water potential regulation, stomatal behaviour and hydraulic transport under drought: Deconstructing the iso/anisohydric concept. Plant Cell Environ. 40, 962–976 (2017).27739594 10.1111/pce.12846
57. Franks PJ Higher rates of leaf gas exchange are associated with higher leaf hydrodynamic pressure gradients Plant Cell Environ. 2006 29 584 592 10.1111/j.1365-3040.2005.01434.x 17080609
Franks, P. J. Higher rates of leaf gas exchange are associated with higher leaf hydrodynamic pressure gradients. Plant Cell Environ. 29, 584–592 (2006).17080609 10.1111/j.1365-3040.2005.01434.x
58. Ampomah-Dwamena C A kiwifruit (Actinidia deliciosa) R2R3-MYB transcription factor modulates chlorophyll and carotenoid accumulation New Phytol. 2019 221 309 325 10.1111/nph.15362 30067292
Ampomah-Dwamena, C. et al. A kiwifruit (Actinidia deliciosa) R2R3-MYB transcription factor modulates chlorophyll and carotenoid accumulation. New Phytol. 221, 309–325 (2019).30067292 10.1111/nph.15362
59. Rouhi V Samson R Lemeur R Damme PV Photosynthetic gas exchange characteristics in three different almond species during drought stress and subsequent recovery Environ. Exp. Bot. 2007 59 117 129 10.1016/j.envexpbot.2005.10.001
Rouhi, V., Samson, R., Lemeur, R. & Damme, P. V. Photosynthetic gas exchange characteristics in three different almond species during drought stress and subsequent recovery. Environ. Exp. Bot. 59, 117–129 (2007).10.1016/j.envexpbot.2005.10.001
60. Li J Defense-related transcription factors WRKY70 and WRKY54 modulate osmotic stress tolerance by regulating stomatal aperture in Arabidopsis New Phytol. 2013 200 457 472 10.1111/nph.12378 23815736
Li, J. et al. Defense-related transcription factors WRKY70 and WRKY54 modulate osmotic stress tolerance by regulating stomatal aperture in Arabidopsis. New Phytol. 200, 457–472 (2013).23815736 10.1111/nph.12378
61. Liu H Signaling transduction of ABA, ROS, and Ca2+ in plant stomatal closure in response to drought Int. J. Mol. Sci. 2022 23 14824 10.3390/ijms232314824 36499153
Liu, H. et al. Signaling transduction of ABA, ROS, and Ca2+ in plant stomatal closure in response to drought. Int. J. Mol. Sci. 23, 14824 (2022).36499153 10.3390/ijms232314824
62. Lim C Baek W Jung J Kim J-H Lee S Function of ABA in stomatal defense against biotic and drought stresses Int. J. Mol. Sci. 2015 16 15251 15270 10.3390/ijms160715251 26154766
Lim, C., Baek, W., Jung, J., Kim, J.-H. & Lee, S. Function of ABA in stomatal defense against biotic and drought stresses. Int. J. Mol. Sci. 16, 15251–15270 (2015).26154766 10.3390/ijms160715251
63. Takahashi, F., Kuromori, T., Sato, H. & Shinozaki, K. Regulatory gene networks in drought stress responses and resistance in plants. Surv. Strateg. Extreme Cold Desiccation. 189–214 (2018).
64. Zhang Y Wheat ABA receptor TaPYL5 constitutes a signaling module with its downstream partners TaPP2C53/TaSnRK2.1/TaABI1 to modulate plant drought response Int. J. Mol. Sci. 2023 24 7969 10.3390/ijms24097969 37175676
Zhang, Y. et al. Wheat ABA receptor TaPYL5 constitutes a signaling module with its downstream partners TaPP2C53/TaSnRK2.1/TaABI1 to modulate plant drought response. Int. J. Mol. Sci. 24, 7969 (2023).37175676 10.3390/ijms24097969
65. Lim CW Yang SH Shin KH Lee SC Kim SH The AtLRK10L1.2, Arabidopsis ortholog of wheat LRK10, is involved in ABA-mediated signaling and drought resistance Plant Cell Rep. 2015 34 447 455 10.1007/s00299-014-1724-2 25533478
Lim, C. W., Yang, S. H., Shin, K. H., Lee, S. C. & Kim, S. H. The AtLRK10L1.2, Arabidopsis ortholog of wheat LRK10, is involved in ABA-mediated signaling and drought resistance. Plant Cell Rep. 34, 447–455 (2015).25533478 10.1007/s00299-014-1724-2
66. Huang S CBL4-CIPK5 pathway confers salt but not drought and chilling tolerance by regulating ion homeostasis Environ. Exp. Bot. 2020 179 104230 10.1016/j.envexpbot.2020.104230
Huang, S. et al. CBL4-CIPK5 pathway confers salt but not drought and chilling tolerance by regulating ion homeostasis. Environ. Exp. Bot. 179, 104230 (2020).10.1016/j.envexpbot.2020.104230
67. Hasanuzzaman M Zhou M Shabala S How does stomatal density and residual transpiration contribute to osmotic stress tolerance? Plants. 2023 12 494 10.3390/plants12030494 36771579
Hasanuzzaman, M., Zhou, M. & Shabala, S. How does stomatal density and residual transpiration contribute to osmotic stress tolerance?. Plants. 12, 494 (2023).36771579 10.3390/plants12030494
68. Yeats TH Rose JK The formation and function of plant cuticles Plant Physiol. 2013 163 5 20 10.1104/pp.113.222737 23893170
Yeats, T. H. & Rose, J. K. The formation and function of plant cuticles. Plant Physiol. 163, 5–20 (2013).23893170 10.1104/pp.113.222737
69. Grünhofer P Schreiber L Cutinized and suberized barriers in leaves and roots: Similarities and differences J. Plant Physiol. 2023 282 153921 10.1016/j.jplph.2023.153921 36780757
Grünhofer, P. & Schreiber, L. Cutinized and suberized barriers in leaves and roots: Similarities and differences. J. Plant Physiol. 282, 153921 (2023).36780757 10.1016/j.jplph.2023.153921
70. Kunst L Biosynthesis and secretion of plant cuticular wax Prog. Lipid Res. 2003 42 51 80 10.1016/S0163-7827(02)00045-0 12467640
Kunst, L. Biosynthesis and secretion of plant cuticular wax. Prog. Lipid Res. 42, 51–80 (2003).12467640 10.1016/S0163-7827(02)00045-0
71. Shellakkutti N Cuticular transpiration is not affected by enhanced wax and cutin amounts in response to osmotic stress in barley Physiol. Plant. 2022 174 e13735 10.1111/ppl.13735 35716005
Shellakkutti, N. et al. Cuticular transpiration is not affected by enhanced wax and cutin amounts in response to osmotic stress in barley. Physiol. Plant. 174, e13735 (2022).35716005 10.1111/ppl.13735
72. Samuels L Kunst L Jetter R Sealing plant surfaces: Cuticular wax formation by epidermal cells Annu. Rev. Plant Biol. 2008 59 683 707 10.1146/annurev.arplant.59.103006.093219 18251711
Samuels, L., Kunst, L. & Jetter, R. Sealing plant surfaces: Cuticular wax formation by epidermal cells. Annu. Rev. Plant Biol. 59, 683–707 (2008).18251711 10.1146/annurev.arplant.59.103006.093219
73. He J ECERIFERUM1-6A is required for the synthesis of cuticular wax alkanes and promotes drought tolerance in wheat Plant Physiol. 2022 190 1640 1657 10.1093/plphys/kiac394 36000923
He, J. et al. ECERIFERUM1-6A is required for the synthesis of cuticular wax alkanes and promotes drought tolerance in wheat. Plant Physiol. 190, 1640–1657 (2022).36000923 10.1093/plphys/kiac394
74. Lee SB Suh MC Recent advances in cuticular wax biosynthesis and its regulation in Arabidopsis Mol. Plant. 2013 6 246 249 10.1093/mp/sss159 23253604
Lee, S. B. & Suh, M. C. Recent advances in cuticular wax biosynthesis and its regulation in Arabidopsis. Mol. Plant. 6, 246–249 (2013).23253604 10.1093/mp/sss159
75. Lewandowska M Keyl A Feussner I Wax biosynthesis in response to danger: Its regulation upon abiotic and biotic stress New Phytol. 2020 227 698 713 10.1111/nph.16571 32242934
Lewandowska, M., Keyl, A. & Feussner, I. Wax biosynthesis in response to danger: Its regulation upon abiotic and biotic stress. New Phytol. 227, 698–713 (2020).32242934 10.1111/nph.16571
76. Sanjari S Molecular, chemical, and physiological analyses of sorghum leaf wax under post-flowering drought stress Plant Physiol. Biochem. 2021 159 383 391 10.1016/j.plaphy.2021.01.001 33450508
Sanjari, S. et al. Molecular, chemical, and physiological analyses of sorghum leaf wax under post-flowering drought stress. Plant Physiol. Biochem. 159, 383–391 (2021).33450508 10.1016/j.plaphy.2021.01.001
77. Dimopoulos N Drought stress modulates cuticular wax composition of the grape berry J. Exp. Bot. 2020 71 3126 3141 10.1093/jxb/eraa046 31985780
Dimopoulos, N. et al. Drought stress modulates cuticular wax composition of the grape berry. J. Exp. Bot. 71, 3126–3141 (2020).31985780 10.1093/jxb/eraa046
78. Li T TaCER1-1A is involved in cuticular wax alkane biosynthesis in hexaploid wheat and responds to plant abiotic stresses Plant Cell Environ. 2019 42 3077 3091 10.1111/pce.13614 31306498
Li, T. et al. TaCER1-1A is involved in cuticular wax alkane biosynthesis in hexaploid wheat and responds to plant abiotic stresses. Plant Cell Environ. 42, 3077–3091 (2019).31306498 10.1111/pce.13614
79. Wang D Ni Y Liao L Xiao Y Guo Y Poa pratensis ECERIFERUM1 (PpCER1) is involved in wax alkane biosynthesis and plant drought tolerance Plant Physiol. Biochem. 2021 159 312 321 10.1016/j.plaphy.2020.12.032 33421907
Wang, D., Ni, Y., Liao, L., Xiao, Y. & Guo, Y. Poa pratensis ECERIFERUM1 (PpCER1) is involved in wax alkane biosynthesis and plant drought tolerance. Plant Physiol. Biochem. 159, 312–321 (2021).33421907 10.1016/j.plaphy.2020.12.032
80. Wu H Tomato SlCER1–1 catalyzes the synthesis of wax alkanes, increasing drought tolerance and fruit storability Hortic. Res.-Engl. 2022 9 uhac004 10.1093/hr/uhac004
Wu, H. et al. Tomato SlCER1–1 catalyzes the synthesis of wax alkanes, increasing drought tolerance and fruit storability. Hortic. Res.-Engl. 9, uhac004 (2022).10.1093/hr/uhac004
81. Wen H Analysis of wheat wax regulation mechanism by liposome and transcriptome Front. Genet. 2021 12 757920 10.3389/fgene.2021.757920 34938312
Wen, H. et al. Analysis of wheat wax regulation mechanism by liposome and transcriptome. Front. Genet. 12, 757920 (2021).34938312 10.3389/fgene.2021.757920
82. Wu H Expression analysis and functional characterization of CER1 family genes involved in very-long-chain alkanes biosynthesis in Brachypodium distachyon Front Plant Sci. 2019 10 1389 10.3389/fpls.2019.01389 31737015
Wu, H. et al. Expression analysis and functional characterization of CER1 family genes involved in very-long-chain alkanes biosynthesis in Brachypodium distachyon. Front Plant Sci. 10, 1389 (2019).31737015 10.3389/fpls.2019.01389
83. Pascal S Arabidopsis CER1-LIKE1 functions in a cuticular very-long-chain alkane-forming complex Plant Physiol. 2019 179 415 432 10.1104/pp.18.01075 30514726
Pascal, S. et al. Arabidopsis CER1-LIKE1 functions in a cuticular very-long-chain alkane-forming complex. Plant Physiol. 179, 415–432 (2019).30514726 10.1104/pp.18.01075
84. Harb A The effect of drought on transcriptome and hormonal profiles in barley genotypes with contrasting drought tolerance Front. Plant Sci. 2020 11 618491 10.3389/fpls.2020.618491 33424910
Harb, A. et al. The effect of drought on transcriptome and hormonal profiles in barley genotypes with contrasting drought tolerance. Front. Plant Sci. 11, 618491 (2020).33424910 10.3389/fpls.2020.618491
85. Mittler R Zandalinas SI Fichman Y Van Breusegem F Reactive oxygen species signalling in plant stress responses Nat. Rev. Mol. Cell Biol. 2022 23 663 679 10.1038/s41580-022-00499-2 35760900
Mittler, R., Zandalinas, S. I., Fichman, Y. & Van Breusegem, F. Reactive oxygen species signalling in plant stress responses. Nat. Rev. Mol. Cell Biol. 23, 663–679 (2022).35760900 10.1038/s41580-022-00499-2
86. Mittler R ROS are good Trends Plant Sci. 2017 22 11 19 10.1016/j.tplants.2016.08.002 27666517
Mittler, R. ROS are good. Trends Plant Sci. 22, 11–19 (2017).27666517 10.1016/j.tplants.2016.08.002
87. Apel K Hirt H Reactive oxygen species: Metabolism, oxidative stress, and signal transduction Annu. Rev. Plant Biol. 2004 55 373 399 10.1146/annurev.arplant.55.031903.141701 15377225
Apel, K. & Hirt, H. Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu. Rev. Plant Biol. 55, 373–399 (2004).15377225 10.1146/annurev.arplant.55.031903.141701
88. Moustafa-Farag M Role of melatonin in plant tolerance to soil stressors: Salinity, pH and heavy metals Molecules. 2020 25 5359 10.3390/molecules25225359 33212772
Moustafa-Farag, M. et al. Role of melatonin in plant tolerance to soil stressors: Salinity, pH and heavy metals. Molecules. 25, 5359 (2020).33212772 10.3390/molecules25225359
89. Hong X Ascorbate peroxidase 1 allows monitoring of cytosolic accumulation of effector-triggered reactive oxygen species using a luminol-based assay Plant Physiol. 2023 191 1416 1434 10.1093/plphys/kiac551 36461917
Hong, X. et al. Ascorbate peroxidase 1 allows monitoring of cytosolic accumulation of effector-triggered reactive oxygen species using a luminol-based assay. Plant Physiol. 191, 1416–1434 (2023).36461917 10.1093/plphys/kiac551
90. Amancio S Photosynthetic and ascorbate-glutathione metabolism in the flag leaves as compared to spikes under drought stress of winter wheat (Triticum aestivum L.) PLoS One. 2018 13 e0194625 10.1371/journal.pone.0194625 29566049
Amancio, S. et al. Photosynthetic and ascorbate-glutathione metabolism in the flag leaves as compared to spikes under drought stress of winter wheat (Triticum aestivum L.). PLoS One. 13, e0194625 (2018).29566049 10.1371/journal.pone.0194625
91. Cheng H Ma X Jia S Li M Mao P Transcriptomic analysis reveals the changes of energy production and AsA-GSH cycle in oat embryos during seed ageing Plant Physiol. Biochem. 2020 153 40 52 10.1016/j.plaphy.2020.03.054 32474385
Cheng, H., Ma, X., Jia, S., Li, M. & Mao, P. Transcriptomic analysis reveals the changes of energy production and AsA-GSH cycle in oat embryos during seed ageing. Plant Physiol. Biochem. 153, 40–52 (2020).32474385 10.1016/j.plaphy.2020.03.054
92. Chang HY Enzymatic characterization and crystal structure analysis of Chlamydomonas reinhardtii dehydroascorbate reductase and their implications for oxidative stress Plant Physiol. Biochem. 2017 120 144 155 10.1016/j.plaphy.2017.09.026 29028546
Chang, H. Y. et al. Enzymatic characterization and crystal structure analysis of Chlamydomonas reinhardtii dehydroascorbate reductase and their implications for oxidative stress. Plant Physiol. Biochem. 120, 144–155 (2017).29028546 10.1016/j.plaphy.2017.09.026
93. Han SH Park YJ Park CM Light primes the thermally induced detoxification of reactive oxygen species during development of thermotolerance in Arabidopsis Plant Cell Physiol. 2019 60 230 241 10.1093/pcp/pcy206 30329122
Han, S. H., Park, Y. J. & Park, C. M. Light primes the thermally induced detoxification of reactive oxygen species during development of thermotolerance in Arabidopsis. Plant Cell Physiol. 60, 230–241 (2019).30329122 10.1093/pcp/pcy206
94. Li HB A cotton ascorbate peroxidase is involved in hydrogen peroxide homeostasis during fibre cell development New Phytol. 2007 175 462 471 10.1111/j.1469-8137.2007.02120.x 17635221
Li, H. B. et al. A cotton ascorbate peroxidase is involved in hydrogen peroxide homeostasis during fibre cell development. New Phytol. 175, 462–471 (2007).17635221 10.1111/j.1469-8137.2007.02120.x
95. Sun Y WRKY33-PIF4 loop is required for the regulation of H2O2 homeostasis Biochem. Biophys. Res. Commun. 2020 527 922 928 10.1016/j.bbrc.2020.05.041 32423827
Sun, Y. et al. WRKY33-PIF4 loop is required for the regulation of H2O2 homeostasis. Biochem. Biophys. Res. Commun. 527, 922–928. 10.1016/j.bbrc.2020.05.041 (2020).32423827 10.1016/j.bbrc.2020.05.041
96. Pal P Genome-wide transcriptional profiling and physiological investigation elucidating the molecular mechanism of multiple abiotic stress response in Stevia rebaudiana Bertoni Sci. Rep.-Uk. 2023 13 19853 10.1038/s41598-023-46000-7
Pal, P. et al. Genome-wide transcriptional profiling and physiological investigation elucidating the molecular mechanism of multiple abiotic stress response in Stevia rebaudiana Bertoni. Sci. Rep.-Uk. 13, 19853 (2023).10.1038/s41598-023-46000-7
97. Zhao, F., Wu, Y., Ma, S., Lei, X. & Liao, W. Increased water use efficiency in China and its drivers during 2000–2016. Ecosystems. 1–17 (2022).
98. Song X Zhou G He Q Zhou H Stomatal limitations to photosynthesis and their critical Water conditions in different growth stages of maize under water stress Agric. Water Manag. 2020 241 106330 10.1016/j.agwat.2020.106330
Song, X., Zhou, G., He, Q. & Zhou, H. Stomatal limitations to photosynthesis and their critical Water conditions in different growth stages of maize under water stress. Agric. Water Manag. 241, 106330 (2020).10.1016/j.agwat.2020.106330
99. Jenkin, S. & Molina, I. Isolation and compositional analysis of plant cuticle lipid polyester monomers. J. Vis. Exp. e53386 (2015).
100. Chen S Zhou Y Chen Y Gu J fastp: An ultra-fast all-in-one FASTQ preprocessor Bioinformatics. 2018 34 i884 i890 10.1093/bioinformatics/bty560 30423086
Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 34, i884–i890 (2018).30423086 10.1093/bioinformatics/bty560
101. Grabherr MG Full-length transcriptome assembly from RNA-Seq data without a reference genome Nat. Biotechnol. 2011 29 644 652 10.1038/nbt.1883 21572440
Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652 (2011).21572440 10.1038/nbt.1883
