Skip to main content
Download PDF

Transcriptome and metabolome analysis reveals key genes and secondary metabolites of Casuarina equisetifolia ssp. incana in response to drought stress

BMC Plant Biology volume 23, Article number: 200 (2023) Cite this article

Abstract

Casuarina equisetifolia is drought tolerant, salt tolerant, and able to grow in barren environments. It is often used to reduce wind damage, to prevent sand erosion, and to help establish plant communities in tropical and subtropical coastal zones. To determine the basis for its drought tolerance, we conducted transcriptomic and metabolic analyses of young branchlets under a non-drought treatment (D_0h) and 2-, 12-, and 24-h-long drought treatments (D_2h, D_12h, and D_24h). A total of 5033 and 8159 differentially expressed genes (DEGs) were identified in D_2h/D_0h and D_24h/D_0h. These DEGs were involved in plant hormone signal transduction, jasmonic acid (JA) biosynthesis, flavonoid biosynthesis, and phenylpropanoid biosynthesis. A total of 148 and 168 differentially accumulated metabolites (DAMs) were identified in D_12h/D_0h and D_24h/D_0h, which were mainly amino acids, phenolic acids, and flavonoids. In conclusion, C. equisetifolia responds to drought by regulating plant hormone signal transduction and the biosynthesis of JA, flavonoid, and phenylpropanoid. These results increase the understanding of drought tolerance in C. equisetifolia at both transcriptional and metabolic levels and provide new insights into coastal vegetation reconstruction and management.

Peer Review reports

Introduction

Drought is a common abiotic stress affecting plant growth and development [1]. Many physiological and biochemical processes of plants are affected by drought stress, resulting in reactive oxygen species accumulation, membrane structure damage, ion imbalance, and inhibition of enzyme activity, photosynthesis, and respiration [2, 3]. Plants have evolved various strategies to adapt to drought stress, such as morphological, cell physiological, metabolic, and molecular changes [4, 5]. As global climate change has increased the frequency of drought, the mechanisms of plant responses or tolerance to drought have attracted extensive attention of researchers [6, 7].

Plants have the ability to adapt to environmental stress by regulating physiological, molecular, and metabolic processes [3]. To reduce the osmotic and oxidative stress caused by drought, plants accumulate osmotic-regulating substances (e.g., soluble sugar and proline), promote the activity of protective enzymes (e.g., superoxide dismutase, catalase, and peroxidase) and increase the content of non-enzyme free radical scavengers (e.g., reduced glutathione and vitamin C) [8, 9]. Phytohormones help plants to improve their adaptability to drought. The abscisic acid (ABA)-mediated signaling pathway induces stomatal closure and regulates the drought responsive genes, which can improve plant drought tolerance [10, 11]. Jasmonate (JA) also mediates the response of plants to drought stress and participates in the biosynthesis of secondary metabolites involved in defense [12,13,14].

With the development of high-throughput sequencing technology, several recent studies have used RNA sequencing (RNA-seq) to reveal the expression patterns and signal regulation pathways of drought stress-related genes, and to identify the key genes and molecular mechanisms of plant responses to drought stress [1, 15]. In Tamarix psammophila, for example, researchers found 1618 and 2716 differentially expressed genes (DEGs) after 1 and 2 weeks of drought treatment, respectively; these DEGs were mainly involved in the mitogen-activated protein kinase (MAPK) signaling pathway, tryptophan and α-linolenic acid metabolism, and biosynthesis of flavonoid and phenylpropanoid [1]. Analysis of drought-stressed Seriphidium transiliense seedlings by the weighted gene co-expression network (WGCNA) method indicated that transcription factors mainly belonging to WRKY, bHLH, NAC, LEA, AP2/ERF, MYB, GRAS, C2H2, MADS, and bZIP families were important in the response to drought [11].

Complementary to transcriptomics, metabolomic analysis can reveal how the synthesis, decomposition, and transformation of metabolites change in response to external environmental stress [1, 16]. For example, the contents of lipids and organic acids in T. psammophila increased under drought [1]. Arabidopsis responded to drought stress by accumulating flavonoids, amino acids, and lipids [17, 18]. The important metabolites differing in drought-tolerant vs. drought-sensitive sesame plants under drought stress were ABA, amino acids, and organic acids [19]. Metabolomic analysis showed that phenylalanine, oxidized glutathione, and ABA were related to the drought tolerance of Sophora davidii [6]. Moreover, the biosynthesis and accumulation of unsaturated fatty acids were found to be involved in the response of maize and peanut to drought stress [20, 21].

Casuarina equisetifolia is a pioneer tree species commonly used for wind prevention and sand fixation in Australia, the Pacific Islands, tropical America and Southeast Asia [22, 23]. C. equisetifolia has strong drought and barren tolerance, therefore it can grow in impoverished and arid coastal zones [24]. Previous studies showed that the maximal photochemical efficiency and the activity of ribulose bisphosphate carboxylase oxygenase of C. equisetifolia were not affected by water stress [25, 26]. C. equisetifolia could adapt to nutrient deficiency and drought stress by increasing the ratios of total phenolics : nitrogen and extractable condensed tannins : nitrogen [27]. However, the molecular mechanisms of C. equisetifolia responses to drought stress are still unknown. In this study, we used transcriptomic and metabolomic analyses to study the drought response of C. equisetifolia branchlets. We identified several key genes and metabolic pathways involved in the drought response of C. equisetifolia. This study thereby increases our understanding of the mechanisms of C. equisetifolia response to drought and may also provide useful information for the management and restoration of zonal vegetation.

Materials and methods

Plant materials and experimental design

C. equisetifolia ssp. incana was selected for this study, because it has high drought, wind, and salt tolerance, and has been widely planted in southeast coastal areas of China [23]. The seedlings of C. equisetifolia were preserved and cultivated in South China Botanical Garden, Chinese Academy of Sciences (Tianhe District, Guangzhou, Guangdong Province, China). After the germination of C. equisetifolia seeds in a 1:1 (v/v) mixture of vermiculite and perlite, its seedlings were transferred to plastic pots (17.5 cm inner diameter and 17.5 cm height; 2 seedlings/pot) containing a 3:1 (v/v) mixture of sand and coconut coir. The seedlings with a height of 20 cm were used in this study. The used vermiculite, perlite, sand, and coconut coir were sterilized by gamma ray irradiation. As not all seedlings form mycorrhizal or nodule with fungi or Frankia, sterilization was carried in order to ensure that all seedlings grow under uniform biological medium environment. In addition, this study mainly focused on the drought resistance of C. equisetifolia excluding the contribution of root symbioses [28].

Drought treatment was implemented by controlling irrigation. The volume water content (VWC) of the medium was checked three times a day. The control and drought treatment corresponded to 45% and 10% respectively. The soil moisture sensor (Hydrosense II) equipped with a CS659 probe (Campbell Scientific, USA) was used to determine the VWC of the medium [29]. The branchlets of C. equisetifolia seedlings were collected at four time points: when the VWC reached 45% as the control (D_0h), and respectively 2 h (D_2h), 12 h (D_12h) and 24 h (D_24h) after the VWC reached 10% as the drought treatment. The time course of previous studies on tolerance of halophytes is within several hours (0–48 h), which is conducive to revealing the response of plants to abiotic stress at the gene level [30, 31]. C. equisetifolia has been proved to be excellent in drought tolerance and can make rapid response to environmental stress [32]. Moreover, our preliminary experiment (tested by quantitative real-time PCR, qRT-PCR) showed that that 2 and 24 h were sensitive time points for C. equisetifolia to respond to drought stress. Thus, in order to understand its molecular response to drought stress, the time course of this study refers to the previous studies of halophytes. There were three biological repetitions at each stage. After the branchlets of the C. equisetifolia seedlings were collected, they were immediately frozen in liquid nitrogen and stored at -80 ℃ for RNA extraction and metabolite profiling.

RNA extraction and sequencing

The total RNA of branchlets was extracted by TRIzol (Invitrogen, USA), and the genomic DNA was removed by DNase I (Takara, Japan). The quality of RNA samples was evaluated with a 2100 Bioanalyser (Agilent, USA) and a NanoDrop ND-2000 (Thermo Fisher Scientific, USA). The RNA samples meeting the quality criteria (OD260 / 280 = 1.8–2.2, OD 260 / 230 ≥ 2.0, RIN ≥ 6.5, 28 S: 18 S ≥ 1.0, > 2 µg) were used to construct libraries and sequences. Sequencing was completed at Shanghai Majorbio Bio-pharm Co., Ltd. (http://www.majorbio.com/) using an Illumina NovaSeq 6000 System (Illumina, USA). SepPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle) software were used for the quality control of the raw data. The data were analyzed on a Major Bio Cloud Platform (https://cloud.majorbio.com/), and the reference genome version was fafu_v1 (http://forestry.fafu.edu.cn/db/Casuarinaceae/index.php). Fragments per kilobase per million (FPKM) was used to normalize and calculate gene expression levels. Raw counts were analyzed for DEGs using DESeq2 based on a negative binomial distribution with P-adjust < 0.05 and |log2FC| ≥ 1. Gene ontology enrichment was analyzed by Goatools (https://github.com/tanghaibao/goatools) and the Fisher test with P-adjust < 0.05. The National Center for Biotechnology Information (NCBI), Non-redundant (Nr, http://ftp.ncbi.nih.gov/blast/db/FASTA/nr.gz), Kyoto Encyclopedia of Genes and Genomes (KEGG, https://www.kegg.jp/), and Gene Ontology (GO, http://geneontology.org/) databases were used for gene function annotation.

Metabolome extraction and analysis

The metabolites were extracted and analyzed by Wuhan Metware Co., Ltd. (https://www.metware.cn/). The freeze-dried branchlets were ground at 30 Hz in a mixing mill MM400 (Retsch, Germany) for 15 min. The powder (100 mg) was extracted overnight with 1 ml of 70% methanol (v/v) at 4 °C. After centrifugation at 10,000 g for 10 min, each supernatant was absorbed (Cnwbond Carbon-GCB SPE Cartridge, Anpel, China) and passed through a 0.22-µm pore size filter. The extracted samples were collected and subjected to metabolome analysis by an UPLC system (Shim-pack UFLC Shimadzu CBM30A, Shimadzu, Japan) and a tandem mass spectrometry (MS/MS) system (Applied Biosystems 6500 Q TRAP, Thermo Fisher Scientific, USA). The working parameters of the UPLC and MS were the method as described by [33]. Based on the Metware database, metabolites were identified according to secondary spectrum information, and metabolites were quantified by multi reaction monitoring of triple quadrupole mass spectrometry. Metabolite abundance was quantified according to peak area. The following thresholds were used to determine whether metabolites differed among treatments (D_0h, D_12h, and D_24h): variable importance in projection (VIP) ≥ 1, P-value < 0.05, fold-change ≥ 2, and fold-change ≤ 0.5]

Quantitative real-time PCR analysis

The remaining RNA samples (D_0h, D_12h, and D_24h) from RNA-seq experiment were used for qRT-PCR analysis. According to the manufacturer’s instructions, the cDNAs were synthesized by the GoScriptTM Reverse Transcription System (Promega, USA). The qRT-PCR was performed with the Unique AptamerTM qPCR SYBR® Green Master Mix (Novogene, China) in a LightCycler® 480 II real-time PCR system (Roche, Switzerland). The PCR cycling condition was 95 °C for 2 min, 40 cycles of 95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s, and the final extension was 60 °C for 15 s [34]. The elongation factor 1-alpha (EF1α) and ubiquitin (UBI) were used as reference genes [35, 36]. The expression level of the target gene in each sample was calculated by the 2−ΔΔCT method [37]. The expression level of each gene in all samples (D_0h, D_2h, and D_24h) was set to 1 in the control (D_0h). Three biological and technical replicates were performed for qRT-PCR. The qRT-PCR primers were designed by the Integrated DNA Technologies (http://www.idtdna.com/Primerquest/Home), and were listed in Table S1.

Results

Transcriptome of C. equisetifolia in response to drought stress

To determine the molecular response of C. equisetifolia to drought, we performed transcriptome sequencing of branchlets collected at three time points (0, 2, and 24 h after drought treatment) with three independent biological replicates for each time point. A total of 403,360,560 clean reads with GC content ranging from 47.15 to 47.88% were generated from RNA-seq, including 131,801,006; 130,745,702; and 140,813,852 clean reads from 0, 2, and 24 h, respectively. The Q30 base percentage of each library was > 93%, indicating that high-quality reads were obtained and could be used for further analysis. An average of 94.67% of clean reads could be mapped to the reference genome, and > 82% of the clean reads from each library were uniquely mapped reads (Table S2 ).

The correlation analysis of all transcriptome samples showed that D_2h samples were strongly and positively correlated with D_0h and D_24h samples (Fig. S1a). Principal component analysis (PCA) was used to reduce the dimension of the data and to visualize the relationship among samples. The first principal component (PC1) and the second principal component (PC2) explained 66.48% and 15.96% of the total variation, respectively (Fig. S1b). Compared with control (D_0h), a total of 5033 DEGs were identified at 2 h (D_2h) of simulated drought stress, with 2299 significantly up-regulated genes and 2734 significantly down-regulated genes. A total of 8159 DEGs were obtained from 0 h vs. 24 h comparisons, consisting of 3724 significantly up-regulated genes and 4435 significantly down-regulated genes (Fig. 1a). For overlap analysis, we identified 4206 DEGs in both the D_2h samples and the D_24h samples (Fig. 1b).

Fig. 1
figure 1

Numbers of DEGs and the main KEGG pathways enriched by DEGs in C. equisetifolia leaves that were exposed to drought stress for 2 h (D_2h) or 24 h (D_24h) relative to leaves that were not exposed to drought stress (D_0h). DEGs were selected based on a cut-off of P-adjust < 0.05 and |log2FC| ≥ 1. (a) Number of up- and down-regulated differentially expressed genes (DEGs) between D_0h and D_2h, and between D_0h and D_24h. (b) Venn analysis of DEGs between D_0h vs. D_2h and D_0h vs. D_24h. (c) Top 20 KEGG pathways enriched by DEGs at 0 h vs. 2 h. (d) Top 20 KEGG pathways enriched by DEGs at 0 h vs. 24 h

Full size image

To explore the function of DEGs, we performed KEGG enrichment analysis (Fig. 1c, d). Relative to D_0h samples, there were three highly representative KEGG pathways in both D_2h and D_24h samples, which were “plant hormone signal transduction” (P-adjust(2 h vs. 0 h) = 0.000, P-adjust(24 h vs. 0 h) = 0.000), “plant-pathogen interaction” (P-adjust(2 h vs. 0 h) = 0.007, P-adjust(24 h vs. 0 h) = 0.042), and “MAPK signaling pathway” (P-adjust(2 h vs. 0 h) = 0.000, P-adjust(24 h vs. 0 h) = 0.000). The numbers of DEGs in each pathway mentioned above were greater in D_24h samples than in D_2h samples. For example, only 75 DEGs were assigned to “plant hormone signal transduction” in D_2h, while 107 DEGs were assigned to “plant hormone signal transduction” in D_24h (Table S3 ). In total, 4206 genes were regarded as DEGs in both pairwise comparisons (0 h vs. 2 h and 0 h vs. 24 h) under drought stress treatment. The expression patterns of the common 4206 DEGs were visualized with a heatmap (Fig. 2a). The result showed that these DEGs could be divided into two main groups: one group had the highest transcript levels in the control (D_0h), while the other group had the highest transcript levels under drought treatment. To understand the general trend of these DEGs, we performed k-means clustering analysis and obtained 6 clusters (Fig. 2b). The top three clusters with the largest number of DEGs were cluster 1 (1789 genes), cluster 3 (1511 genes), and cluster 2 (840 genes) (Fig. 2b). In cluster 1, where transcript levels of 1789 genes were up-regulated after drought stress treatment (Fig. 2b), and the top two pathways were “MAPK signaling pathway” and “plant hormone signal transduction” (Table S4 ). In cluster 2, all of the DEGs exhibited a dramatic decrease during the treatment (Fig. 2b), and the top two pathways in cluster 2 were “plant hormone signal transduction” and “flavonoid biosynthesis” (Table S4 ).

Fig. 2
figure 2

Hierarchical cluster analysis of common DEGs in C. equisetifolia leaves that were exposed to drought stress for different times. DEGs were selected based on a cut-off of P-adjust < 0.05 and |log2FC| ≥ 1. (a) Hierarchical cluster analysis of common DEGs in D_0 h vs. D_2 h and D_0 h vs. D_24 h. Red represents down-regulation and blue represents up-regulation of gene expression. Expression values were Z-scaled (log10 FPKM). (b) The expression patterns of six clusters. The blue line indicates the average expression values of all genes in the gene set

Full size image

To verify the accuracy of transcriptome data during drought stress, six DEGs from pathways of plant hormone signal transduction, phenylpropanoid and flavonoid biosynthesis were selected for qRT-PCR. The results of RNA-seq showed that three DEGs were up-regulated while three DEGs were down-regulated (Fig. S2). The expression trends of six DEGs were consistent with the transcriptomic data, which verifies its reliability and accuracy.

Metabolomics of C. equisetifolia in response to drought stress

To understand the changes of metabolites under drought stress in C. equisetifolia, we conducted a wide analysis of metabolites on D_0h, D_12h, and D_24h samples, and identified the differentially accumulated metabolites (DAMs). In comparison to the control (D_0h), 148 DAMs were found in D_12h samples, with 123 significantly up-regulated and 25 significantly down-regulated metabolites (Fig. 3a). A total of 168 metabolites were identified as DAMs in D_24h samples, consisting of 145 significantly up-regulated and 23 significantly down-regulated metabolites (Fig. 3b). To reveal the pathways of DAMs, we performed a KEGG pathway enrichment analysis. In D_12h samples, the top five pathways were “biosynthesis of secondary metabolites”, “biosynthesis of amino acids”, “phenylpropanoid biosynthesis”, “aminoacyl-tRNA biosynthesis”, and “2-oxocarboxylic acid metabolism” (Fig. 3c). However, the top five pathways in D_24h samples were “biosynthesis of secondary metabolites”, “biosynthesis of amino acids”, “2-oxocarboxylic acid metabolism”, “ABC transporters”, “aminoacyl-tRNA biosynthesis”, and “phenylpropanoid biosynthesis” (Fig. 3d). The contents of many compounds increased significantly under drought stress, and these compounds were mainly involved in amino acid biosynthesis, phenylpropanoid biosynthesis, and flavonoid biosynthesis. In the metabolic pathway of amino acid biosynthesis, the contents of methionine, leucine, isoleucine, proline, and phenylalanine increased significantly.

Fig. 3
figure 3

Differentially accumulated metabolites (DAMs) at different times in C. equisetifolia under drought stress. DAMs were selected based on a cut-off of VIP ≥ 1, P-value < 0.05, fold-change ≥ 2, and fold-change ≤ 0.5. Volcano plot of DAMs in (a) D_0h vs. D_12h and in (b) D_0h vs. D_24h. (c) Top 20 enriched pathways for DAMs of D_0h vs. D_12h. (d) Top 20 enriched pathways for DAMs of D_0h vs. D_24h

Full size image

Plant hormone signal transduction pathways involved in the response to drought stress

The “plant hormone signal transduction” pathway was highly represented in the KEGG enrichment analysis (Fig. 1c, d). To analyze the relationship between plant hormones and the drought stress response in C. equisetifolia, we investigated the genes related to the plant hormone signal transduction under drought stress using RNA-seq data. Four plant hormone signal transduction pathways, including ABA, auxin, cytokinin, and JA, were involved in the drought stress response (Fig. 4). For the ABA signal transduction pathway, 4 DEGs encoding ABA receptors pyrabactin resistance/pyrabactin resistance-like (PRY/PYL) proteins were significantly down-regulated, while the 7 DEGs encoding type 2 C protein phosphatase (PP2C), 2 DEGs encoding ABA-activated SNF1-related protein kinases 2 (SnRK2), and 1 DEG encoding ABA-responsive element binding factor (ABF) were significantly up-regulated, suggesting the importance of ABA in the drought stress response. In the JA signal transduction pathway, 1 DEG encoding coronatine insensitive 1 (COI1, significantly down-regulated), 2 DEGs encoding jasmonate ZIM-domain (JAZ, significantly up-regulated), and 2 DEGs encoding MYC2 (1 significantly down-regulated and 1 significantly up-regulated) were expressed significantly under drought stress (Fig. 4). In addition, the 8 DEGs involved in JA biosynthesis were significantly up-regulated, including 3 DEGs encoding PLA2G/PLA2G16, 1 DEG encoding LOX2S, 1 DEG encoding 12-oxo-phytodienoic acid reductase (OPR), 1 DEG encoding OPC-8:0 CoA Ligase1 (OPCL1), 1 DEG encoding ACAA1, and 1 DEG encoding jasmonate O-methyltransferase (Fig. 5).

Fig. 4
figure 4

DEGs involved in plant hormone signal transduction at different times in C. equisetifolia under drought stress. The original figure of plant hormone signal transduction pathway was cited from the Kanehisa laboratories [69]. DEGs were selected based on a cut-off of P-adjust < 0.05 and |log2FC| ≥ 1. The two boxes represent log2FC (D_2h / D_0h) and log2FC (D_24h / D_0h), respectively. Red text indicates significantly up-regulated DEGs, green text indicates significantly down-regulated DEGs, and blue text indicates both significantly up-regulated and down-regulated DEGs.

AUX1: Auxin-Resistant1. AUX/IAA: Auxin/Indole-3-Acetic Acid. ARF: auxin response factor. GH3: Gretchen hagen 3. SAUR: Small auxin up RNA. CRE1: Cytokinin receptor (1) AHP: histidine phosphotransfer proteins. A-ARR: type-A response regulator. PYL: pyrabactin resistance-like. PP2C: phosphatase 2 C. SnRK2: ABA-activated SNF1-related protein kinases (2) ABF: ABA-responsive element binding factor. COI1: coronatine insensitive 1. JAZ: jasmonate ZIM-domain. MYC2: transcription factor

Full size image
Fig. 5
figure 5

DEGs and DAMs in the jasmonic acid (JA) biosynthetic pathway at different times in C. equisetifolia under drought stress. The original figure of JA biosynthetic pathway was cited from the Kanehisa laboratories [70, 71, 72]. DEGs were selected based on a cut-off of P-adjust < 0.05 and |log2FC| ≥ 1. DAMs were selected based on a cut-off of VIP ≥ 1, P-value < 0.05, fold-change ≥ 2, and fold-change ≤ 0.5. The two boxes represent log2FC (D_2h / D_0h) and log2FC (D_24h / D_0h) of DEGs, respectively. The two triangles represent log2FC (D_2h / D_0h) and log2FC (D_24h / D_0h) of DAMs, respectively. Red text represents significantly up-regulated DEGs/DAMs. AOS: allene oxide synthase. OPR: 12-oxo-phytodienoic acid reductase. OPCL1: OPC-8:0 CoA Ligase1. ACAA1: acetyl-CoA acyltransferase 1

Full size image

Secondary metabolism involved in the response to drought stress

Drought treatments (both 12 and 24 h) significantly increased the accumulation of methionine, leucine, isoleucine, and proline in C. equisetifolia. The significant accumulation of its compounds in phenylpropanoid and flavonoid biosynthesis was consistent with the DEGs whose expression levels were significantly up-regulated in these pathways. For example, the shikimate O-hydroxycinnamoyl transferase gene HCT was up-regulated more than 1.98-fold at 2 h and 1.02-fold at 24 h after drought stress in C. equisetifolia (Fig. 6). The product of this enzyme, namely p-coumaryl shikimic acid, was significantly up-regulated after drought stress (Fig. 6). The significant increase of dihydrokaempferol and sinapyl alcohol is likely the result of the significant up-regulation of these DEGs under drought stress (Fig. 6).

Fig. 6
figure 6

DEGs and DAMs in amino acid, phenylpropanoid, and flavonoid biosynthetic pathways at different times under drought stress. The original figures of amino acid, phenylpropanoid, and flavonoid biosynthetic pathways were cited from the Kanehisa laboratories [70, 71, 72]. DEGs were selected based on a cut-off of P-adjust < 0.05 and |log2FC| ≥ 1. DAMs were selected based on a cut-off of VIP ≥ 1, P-value < 0.05, fold-change ≥ 2, and fold-change ≤ 0.5. The two boxes represent log2FC (D_2h / D_0h) and log2FC (D_24h / D_0h) of DEGs, respectively. The two triangles represent log2FC (D_2h / D_0h) and log2FC (D_24h / D_0h) of DAMs, respectively. Red text represents significantly up-regulated DEGs/DAMs. 4CL: 4-coumarate: coenzyme A ligase. HCT: hydroxycinnamoyl transferase. CCR: cinnamoyl-coA reductase. CAD: cinnamoyl alcohol dehydrogenase. F3H: Flavanone 3-hydroxylase

Full size image

Discussion

Plant hormone signal transduction pathways in response to drought stress

Many studies of plants have shown that ABA- and JA-mediated signal transduction pathways are related to the response to drought stress [38,39,40]. It has been reported that the ABA and JA signaling pathways of C. equisetifolia and Casuarina glauca were activated under salt stress [41, 42]. ABA that accumulates under drought stress can activate specific serine / threonine kinase family SnRK2, which can phosphorylate downstream transcription factors and activate ABA response genes [43, 44]. In this study, there were 2 DEGs encoding SnRK2 significantly up-regulated in C. equisetifolia subjected to both 2 and 24 h of drought stress. PYL proteins are ABA receptors that are involved in the response of plants to ABA [45]. Overexpression of PYLs in Arabidopsis can promote its response to ABA, but some PYLs are down-regulated after exogenous ABA treatment [46, 47]. The expression pattern of PYLs in Arabidopsis can be explained by negative feedback regulation [46, 48]. Our results showed that 4 PYLs were significantly down-regulated after drought treatment, which was consistent with the previous studies. Researchers have also reported that drought stress activated downstream transcription factors involved in ABA response, and the activated transcription factors may in return regulate the expression of PYLs and PP2Cs [45]. In the current study, ABA response mediated downregulation of receptor PYLs, upregulation of negative regulator PP2Cs, and upregulation of the downstream transcription factor ABF. The negative feedback can reduce the adverse effects caused by ABA accumulation, which may help to increase the adaptability of C. equisetifolia to drought stress.

JAs not only regulate plant growth and development, but also improve plant tolerance to environmental stress via the JA signaling pathway [49]. JA accumulation induced by stress activates downstream transcription factors and then activates regulatory genes in response to stress [50]. The expression of COI1 in Poa pratensis leaves was recently found to be inhibited at the initial stage of drought stress [51]. We found that COI1 expression was significantly down-regulated in C. equisetifolia, suggesting that the JA signaling pathway was involved in the response to drought. Previous studies have shown that overexpression of JAZ6 in Oryza sativa and JAZs in Poa pratensis can increase the tolerance of those plants to mannitol and drought stress [51, 52]. There were 2 up-regulated DEGs encoding JAZs in treated C. equisetifolia branchlets, which may help improve the drought tolerance of the plant. It is reported that the MYC2 transcription level in Arabidopsis decreased significantly under drought stress [53]. Our results showed that the expression of 1 DEG encoding MYC2 was significantly down regulated, which was consistent with the results of previous studies. These results suggested that the plant hormones might be important in the drought stress response of C. equisetifolia.

The JA biosynthesis pathway in response to drought stress

JA can increase the tolerance of plants to water stress by regulating stomata, scavenging reactive oxygen species, and altering root development [50]. Previous studies have shown that continuous drought stress promoted the up-regulation of the JA biosynthesis genes OPR, ACX, and ACAA1 in P. pratensis [51] and AOS and OPR in Cicer arietinum [54]. The JA biosynthesis gene ACAA1 is related to peroxisome, and the up-regulation of its expression may help remove reactive oxygen species produced by drought stress [51]. Our results showed that the expressions of the JA biosynthesis genes OPR, AOS, and ACAA1 were significantly up-regulated, which is consistent with previous studies [51, 54]. One of the substrates for the synthesis of JA, 12-OPDA, is also a signal molecule with overlapping and different functions from JA [55, 56]. Previous studies have shown that the accumulation of 12-OPDA reduces stomatal apertures and increases plant drought tolerance [55]. We also found that 12-OPDA of C. equisetifolia accumulated significantly under drought stress, which is consistent with previous studies. In C. equisetifolia, the significantly up-regulated genes and significantly increased metabolites in JA biosynthesis pathway may improve drought tolerance in this species.

Secondary metabolism in response to drought stress

Amino acids are osmotic regulators that maintain and stabilize cell structure and that usually accumulate in plants in response to osmotic stress [57]. Previous studies have shown that plants (O. sativa, Zea mays, and Arabidopsis) accumulate a large number of amino acids (such as proline, leucine, tryptophan, isoleucine, and histidine) to increase drought tolerance [18, 58, 59]. In addition, C. glauca accumulates ornithine and proline under salt stress to regulate osmotic pressure [60]. In this study, the contents of proline, leucine, isoleucine, methionine, and phenylalanine were significantly increased in response to drought stress, indicating that amino acid metabolism was involved in the response of C. equisetifolia to drought stress.

Plants respond to oxidative damage induced by osmotic stress by accumulating active oxygen scavengers such as phenols and flavonoids [61]. Many studies have shown that the accumulation of phenols and flavonoids in plants can improve drought tolerance [62, 63]. It has been reported that the biosynthesis pathway of secondary metabolites of C. glauca was activated under salt stress, and the increase of flavonoid content alleviated oxidative damage [28, 64,65,66]. Our results showed that the pathways of phenylpropanoid and flavonoid biosynthesis in C. equisetifolia respond to drought stress at the transcriptomic and metabolomic levels. The accumulation of one flavonoid (dihydrokaempferol) and five phenolic acids (p-coumaroyl shikimic acid, caffeoyl shikimic acid, caffeic acid, ferulic acid, and sinapyl alcohol) in the current study was significantly increased under drought stress, which was consistent with the results of previous studies. Previous studies have shown that overexpression of F3H (an important gene in the flavonoid biosynthesis pathway) and 4CL (4-coumaric acid CoA ligase, a branch point in the phenylpropanoid biosynthesis pathway) can increase the drought tolerance of plants [67,68,69]. We therefore suspect that flavonoids, phenols, and related genes may be involved in the response of C. equisetifolia to drought stress.

Conclusion

We analyzed the transcript and metabolic profiles of C. equisetifolia branchlets to increase our understanding of the mechanisms underlying the plant’s responses to drought stress. Pathway analysis showed that DEGs under drought stress were related to ABA and JA signal transduction pathways, phenylpropanoid biosynthesis, and flavonoid biosynthesis. Drought stress not only induced the expression of genes related to JA biosynthesis but also increased the accumulation of flavonoids and phenols and the expression of related genes. Our results also showed that drought stress increased the accumulation of amino acids, such as leucine, isoleucine, proline, methionine, and phenylalanine, which could help C. equisetifolia regulate osmotic pressure and maintain cell structure under drought stress. These results increase our understanding of how C. equisetifolia tolerates or resists drought stress at both transcriptional and metabolic levels and provide reference data for the management and restoration of regional vegetation.

Data Availability

The sequenced raw reads generated in this study have been submitted to the National Center for Biotechnology Information (NCBI) with BioProject ID: PRJNA902157 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA902157?reviewer=fa84revdsbtmdef8f97ahav6q5).

References

  1. Sun T-T, Su Z-H, Wang R, Liu R, Yang T, Zuo W-T, Wen S-S, Wang L-Q. Transcriptome and metabolome analysis reveals the molecular mechanisms of Tamarix taklamakanensis under progressive drought and rehydration treatments. Environ Exp Bot. 2022;195:104766.

    Article  CAS  Google Scholar 

  2. Anjum SA, Ashraf U, Tanveer M, Khan I, Hussain S, Shahzad B, Zohaib A, Abbas F, Saleem MF, Ali I, Wang LC. Drought induced changes in growth, osmolyte accumulation and antioxidant metabolism of three maize hybrids. Front Plant Sci. 2017;8:69.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Shinozaki K, Yamaguchi-Shinozaki K. Gene networks involved in drought stress response and tolerance. J Exp Bot. 2007;58(2):221–7.

    Article  CAS  PubMed  Google Scholar 

  4. Li B, Fan R, Fan Y, Liu R, Zhang H, Chen T, Liu J, Li H, Zhao X, Song C-P. The flavonoid biosynthesis regulator PFG3 confers drought stress tolerance in plants by promoting flavonoid accumulation. Environ Exp Bot. 2022;196:104792.

    Article  Google Scholar 

  5. Chaves MM, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot. 2009;103(4):551–60.

    Article  CAS  PubMed  Google Scholar 

  6. Zhao X, Huang L-J, Sun X-F, Zhao L-L, Wang P-C. Transcriptomic and metabolomic analyses reveal key metabolites, pathways and candidate genes in Sophora davidii (Franch.) Skeels seedlings under drought stress. Front Plant Sci. 2022;13:785702.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Rahimi M, Kordrostami M, Mortezavi M. Evaluation of tea (Camellia sinensis L.) biochemical traits in normal and drought stress conditions to identify drought tolerant clones. Physiol Mol Biol Plants. 2019;25(1):59–69.

    Article  CAS  PubMed  Google Scholar 

  8. Miao Y, Lv D, Wang P, Wang X-C, Chen J, Miao C, Song C-P. An Arabidopsis glutathione peroxidase functions as both a redox transducer and a scavenger in abscisic acid and drought stress responses. Plant Cell. 2006;18(10):2749–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Lu X-P, Gao H-J, Zhang L, Wang Y-P, Shao K-Z, Zhao Q, Zhang J-L. Dynamic responses of Haloxylon ammodendron to various degrees of simulated drought stress. Plant Physiol Biochem. 2019;139:121–31.

    Article  CAS  PubMed  Google Scholar 

  10. Fujii H, Chinnusamy V, Rodrigues A, Rubio S, Antoni R, Park S-Y, Cutler SR, Sheen J, Rodriguez PL, Zhu J-K. In vitro reconstitution of an abscisic acid signalling pathway. Nature. 2009;462(7273):660–U138.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Liu X, Chen A, Wang Y, Jin G, Zhang Y, Gu L, Li C, Shao X, Wang K. Physiological and transcriptomic insights into adaptive responses of Seriphidium transiliense seedlings to drought stress. Environ Exp Bot. 2022;194:104736.

    Article  CAS  Google Scholar 

  12. Wasternack C, Feussner I. The oxylipin pathways: Biochemistry and function. Annu Rev Plant Biol. 2018;69:363–86.

    Article  CAS  PubMed  Google Scholar 

  13. Howe GA, Major IT, Koo AJ. Modularity in jasmonate signaling for multistress resilience. Annu Rev Plant Biol. 2018;69:387–415.

    Article  CAS  PubMed  Google Scholar 

  14. Chini A, Monte I, Fernandez-Barbero G, Boter M, Hicks G, Raikhel N, Solano R. A small molecule antagonizes jasmonic acid perception and auxin responses in vascular and nonvascular plants. Plant Physiol. 2021;187(3):1399–413.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Joshi R, Wani SH, Singh B, Bohra A, Dar ZA, Lone AA, Pareek A, Singla-Pareek SL. Transcription factors and plants response to drought stress: current understanding and future directions. Front Plant Sci. 2016;7:1029.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gargallo-Garriga A, Sardans J, Perez-Trujillo M, Oravec M, Urban O, Jentsch A, Kreyling J, Beierkuhnlein C, Parella T, Penuelas J. Warming differentially influences the effects of drought on stoichiometry and metabolomics in shoots and roots. New Phytol. 2015;207(3):591–603.

    Article  CAS  PubMed  Google Scholar 

  17. Tarazona P, Feussner K, Feussner I. An enhanced plant lipidomics method based on multiplexed liquid chromatography-mass spectrometry reveals additional insights into cold- and drought-induced membrane remodeling. Plant J. 2015;84(3):621–33.

    Article  CAS  PubMed  Google Scholar 

  18. Pires MV, Pereira Junior AA, Medeiros DB, Daloso DM, Pham PA, Barros KA, Engqvist MKM, Florian A, Krahnert I, Maurino VG, Araujo WL, Fernie AR. The influence of alternative pathways of respiration that utilize branched-chain amino acids following water shortage in Arabidopsis. Plant Cell Environ. 2016;39(6):1304–19.

    Article  CAS  PubMed  Google Scholar 

  19. You J, Zhang Y, Liu A, Li D, Wang X, Dossa K, Zhou R, Yu J, Zhang Y, Wang L, Zhang X. Transcriptomic and metabolomic profiling of drought-tolerant and susceptible sesame genotypes in response to drought stress. BMC Plant Biol. 2019;19:267.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Yang L, Fountain JC, Ji P, Ni X, Chen S, Lee RD, Kemerait RC, Guo B. Deciphering drought-induced metabolic responses and regulation in developing maize kernels. Plant Biotechnol J. 2018;16(9):1616–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Gundaraniya SA, Ambalam PS, Tomar RS. Metabolomic profiling of drought-tolerant and susceptible peanut (Arachis hypogaea L.) genotypes in response to drought stress. Acs Omega. 2020;5(48):31209–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wheeler GS, Taylor GS, Gaskin JF, Purcell MF. Ecology and Management of Sheoak (Casuarina spp.), an invader of Coastal Florida. USA J Coastal Res. 2011;27(3):485–92.

    Article  Google Scholar 

  23. Ye GF, Zhang HX, Chen BH, Nie S, Liu H, Gao W, Wang HY, Gao YB, Gu LF. De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth. Plant J. 2019;97(4):779–94.

    Article  CAS  PubMed  Google Scholar 

  24. Hu P, Zhong C, Zhang Y, Jiang Q, Chen Y, Chen Z, Pinyopusarerk K, Bush D. Geographic variation in seedling morphology of Casuarina equisetifolia subsp equisetifolia (Casuarinaceae). Aust J Bot. 2016;64(2):160–70.

    Article  Google Scholar 

  25. SanchezRodriguez J, MartinezCarrasco R, Perez P. Photosynthetic electron transport and carbon-reduction-cycle enzyme activities under long-term drought stress in Casuarina equisetifolia Forst. & Forst Photosynth Res. 1997;52(3):255–62.

    Article  CAS  Google Scholar 

  26. Martinez-Carrasco R, Sanchez-Rodriguez J, Perez P. Changes in chlorophyll fluorescence during the course of photoperiod and in response to drought in Casuarina equisetifolia Forst. and Forst. Photosynthetica 2002;40(3):363–368.

  27. Zhang LH, Shao HB, Ye GF, Lin YM. Effects of fertilization and drought stress on tannin biosynthesis of Casuarina equisetifolia seedlings branchlets. Acta Physiol Plant. 2012;34(5):1639–49.

    Article  CAS  Google Scholar 

  28. Ribeiro-Barros AI, Pawlowski K, Ramalho JC. Mechanisms of salt stress tolerance in Casuarina: a review of recent research. J For Res. 2022;27(2):113–6.

    Article  CAS  Google Scholar 

  29. Liu J, Shrestha P, Skabelund LR, Todd T, Decker A, Kirkham MB. Growth of prairie plants and sedums in different substrates on an experimental green roof in mid-continental USA. Sci Total Environ. 2019;697:134089.

    Article  CAS  PubMed  Google Scholar 

  30. Tiika RJ, Wei J, Cui GX, Ma YJ, Yang HS, Duan HR. Transcriptome-wide characterization and functional analysis of Xyloglucan endo-transglycosylase/hydrolase (XTH) gene family of Salicornia europaea L. under salinity and drought stress. BMC Plant Biol. 2021;21(1):491.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wang C, Gao C, Wang L, Zheng L, Yang C, Wang Y. Comprehensive transcriptional profiling of NaHCO3-stressed Tamarix hispida roots reveals networks of responsive genes. Plant Mol Biol. 2014;84(1–2):145–57.

    Article  CAS  Google Scholar 

  32. Xu H, Yu J, You L, Xiao S, Nie S, Li T, Ye G, Lin D. Drought resistance evaluation of Casuarina equisetifolia half-sib families at the seedling stage and the response of five NAC genes to drought stress. Forests. 2022;13(12):2037.

    Article  Google Scholar 

  33. Guo H, Guo H, Zhang L, Tang Z, Yu X, Wu J, Zeng F. Metabolome and transcriptome association analysis reveals dynamic regulation of purine metabolism and flavonoid synthesis in transdifferentiation during somatic embryogenesis in cotton. Int J Mol Sci. 2019;20(9):2070.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Si C, Zeng D, Yu Z, Teixeira da Silva JA, Duan J, He C, Zhang J. Transcriptomic and metabolomic analyses reveal the main metabolites in Dendrobium officinale leaves during the harvesting period. Plant physiol Biochem. 2022;190:24–34.

    Article  CAS  PubMed  Google Scholar 

  35. Fan C, Qiu Z, Zeng B, Liu Y, Li X, Guo G. Selection of reference genes for quantitative real-time PCR in Casuarina equisetifolia under salt stress. Biol Plant. 2017;61(3):463–72.

    Article  CAS  Google Scholar 

  36. Veluthakkal R, Dasgupta MG. Isolation and characterization of pathogen defence-related class I chitinase from the actinorhizal tree Casuarina equisetifolia. For Pathol. 2012;42(6):467–80.

    Article  Google Scholar 

  37. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  38. Liu H, Zhang Y-H, Yin H, Wang W-X, Zhao X-M, Du Y-G. Alginate oligosaccharides enhanced Triticum aestivum L. tolerance to drought stress. Plant Physiol Biochem. 2013;62:33–40.

    Article  CAS  PubMed  Google Scholar 

  39. Okamoto M, Peterson FC, Defries A, Park S-Y, Endo A, Nambara E, Volkman BF, Cutler SR. Activation of dimeric ABA receptors elicits guard cell closure, ABA-regulated gene expression, and drought tolerance. Proc. Natl. Acad. Sci. U. S. A. 2013;110(29):12132–12137.

  40. Lu G, Gao C, Zheng X, Han B. Identification of OsbZIP72 as a positive regulator of ABA response and drought tolerance in rice. Planta. 2009;229(3):605–15.

    Article  CAS  PubMed  Google Scholar 

  41. Wang Y, Zhang J, Qiu Z, Zeng B, Zhang Y, Wang X, Chen J, Zhong C, Deng R, Fan C. Transcriptome and structure analysis in root of Casuarina equisetifolia under NaCl treatment. PEERJ. 2021;9:e12133.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Fernandes I, Paulo OS, Marques I, Sarjkar I, Sen A, Graca I, Pawlowski K, Ramalho JC, Ribeiro-Barros AI. Salt stress tolerance in Casuarina glauca: insights from the branchlets transcriptome. Plants-Basel. 2022;11(21):2942.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Du X, Yu R, Shi C, Wang Y, Meng R, Shi W, Jin M, Wei X, Sun T. Comparative transcriptomics reveals osmotic and ionic stress key genes contributing to the difference in the salinity tolerance of two pak choi cultivars. Environ Exp Bot. 2021;191:104621.

    Article  CAS  Google Scholar 

  44. Shen X, Guo X, Zhao D, Zhang Q, Jiang Y, Wang Y, Peng X, Wei Y, Zhai Z, Zhao W, Li T. Cloning and expression profiling of the PacSnRK2 and PacPP2C gene families during fruit development, ABA treatment, and dehydration stress in sweet cherry. Plant Physiol Biochem. 2017;119:275–85.

    Article  CAS  PubMed  Google Scholar 

  45. Chen Y, Feng L, Wei N, Liu ZH, Hu S, Li XB. Overexpression of cotton PYL genes in Arabidopsis enhances the transgenic plant tolerance to drought stress. Plant Physiol Biochem. 2017;115:229–38.

    Article  CAS  PubMed  Google Scholar 

  46. Santiago J, Rodrigues A, Saez A, Rubio S, Antoni R, Dupeux F, Park SY, Marquez JA, Cutler SR, Rodriguez PL. Modulation of drought resistance by the abscisic acid receptor PYL5 through inhibition of clade A PP2Cs. Plant J. 2009;60(4):575–88.

    Article  CAS  PubMed  Google Scholar 

  47. Lee SC, Lim CW, Lan W, He K, Luan S. ABA signaling in guard cells entails a dynamic protein-protein interaction relay from the PYL-RCAR family receptors to ion channels. Mol Plant. 2013;6(2):528–38.

    Article  CAS  PubMed  Google Scholar 

  48. Merlot S, Gosti F, Guerrier D, Vavasseur A, Giraudat J. The ABI1 and ABI2 protein phosphatases 2 C act in a negative feedback regulatory loop of the abscisic acid signalling pathway. Plant J. 2001;25(3):295–303.

    Article  CAS  PubMed  Google Scholar 

  49. Wang Y, Mostafa S, Zeng W, Jin B. Function and mechanism of jasmonic acid in plant responses to abiotic and biotic stresses. Int J Mol Sci. 2021;22(16):8568.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Ababaf M, Omidi H, Bakhshandeh A. Changes in antioxidant enzymes activities and alkaloid amount of Catharanthus roseus in response to plant growth regulators under drought condition. Ind Crop Prod. 2021;167:113505.

    Article  CAS  Google Scholar 

  51. Chen Y, Chen Y, Shi Z, Jin Y, Sun H, Xie F, Zhang L. Biosynthesis and signal transduction of ABA, JA, and BRs in response to drought stress of Kentucky Bluegrass. Int J Mol Sci. 2019;20(6):1289.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Lima JM, Nath M, Dokku P, Raman KV, Kulkarni KP, Vishwakarma C, Sahoo SP, Mohapatra UB, Mithra SVA, Chinnusamy V, Robin S, Sarla N, Seshashayee M, Singh K, Singh AK, Singh NK, Sharma RP. Mohapatra. Physiological, anatomical and transcriptional alterations in a rice mutant leading to enhanced water stress tolerance. Aob Plants. 2015;7:plv023.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Liu N, Staswick PE, Avramova Z. Memory responses of jasmonic acid-associated Arabidopsis genes to a repeated dehydration stress. Plant Cell Environ. 2016;39(11):2515–29.

    Article  CAS  PubMed  Google Scholar 

  54. De Domenico S, Bonsegna S, Horres R, Pastor V, Taurino M, Poltronieri P, Imtiaz M, Kahl G, Flors V, Winter P, Santino A. Transcriptomic analysis of oxylipin biosynthesis genes and chemical profiling reveal an early induction of jasmonates in chickpea roots under drought stress. Plant Physiol Biochem. 2012;61:115–22.

    Article  PubMed  Google Scholar 

  55. Savchenko T, Kolla VA, Wang CQ, Nasafi Z, Hicks DR, Phadungchob B, Chehab WE, Brandizzi F, Froehlich J, Dehesh K. Functional convergence of oxylipin and abscisic acid pathways controls stomatal closure in response to drought. Plant Physiol. 2014;164(3):1151–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Kienow L, Schneider K, Bartsch M, Stuible H-P, Weng H, Miersch O, Wasternack C, Kombrink E. Jasmonates meet fatty acids: functional analysis of a new acyl-coenzyme A synthetase family from Arabidopsis thaliana. J Exp Bot. 2008;59(2):403–19.

    Article  CAS  PubMed  Google Scholar 

  57. Bartels D, Sunkar R. Drought and salt tolerance in plants. Crit Rev Plant Sci. 2005;24(1):23–58.

    Article  CAS  Google Scholar 

  58. Obata T, Witt S, Lisec J, Palacios-Rojas N, Florez-Sarasa I, Yousfi S, Luis Araus J, Cairns JE, Fernie AR. Metabolite profiles of maize leaves in drought, heat, and combined stress field trials reveal the relationship between metabolism and grain yield. Plant Physiol. 2015;169(4):2665–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Ma X, Xia H, Liu Y, Wei H, Zheng X, Song C, Chen L, Liu H, Luo L. Transcriptomic and metabolomic studies disclose key metabolism pathways contributing to well-maintained photosynthesis under the drought and the consequent drought-tolerance in rice. Front Plant Sci. 2016;7:1886.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Jorge TF, Duro N, da Costa M, Florian A, Ramalho JC, Ribeiro-Barros AI, Fernie AR, Antonio C. GC-TOF-MS analysis reveals salt stress-responsive primary metabolites in Casuarina glauca tissues. Metabolomics. 2017;13(8):95.

    Article  Google Scholar 

  61. Sarker U, Oba S. Drought stress effects on growth, ROS markers, compatible solutes, phenolics, flavonoids, and antioxidant activity in Amaranthus tricolor. Appl Biochem Biotechnol. 2018;186(4):999–1016.

    Article  CAS  PubMed  Google Scholar 

  62. Sanchez-Rodriguez E, Rubio-Wilhelmi M, Cervilla LM, Blasco B, Rios JJ, Rosales MA, Romero L, Ruiz JM. Genotypic differences in some physiological parameters symptomatic for oxidative stress under moderate drought in tomato plants. Plant Sci. 2010;178(1):30–40.

    Article  CAS  Google Scholar 

  63. Nakabayashi R, Mori T, Saito K. Alternation of flavonoid accumulation under drought stress in Arabidopsis thaliana. Plant Signal Behav. 2014;9(8):e29518.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Jorge TF, Tohge T, Wendenburg R, Ramalho JC, Lidon FC, Ribeiro-Barros AI, Fernie AR, Antonio C. Salt-stress secondary metabolite signatures involved in the ability of Casuarina glauca to mitigate oxidative stress. Environ Exp Bot. 2019;166:103808.

    Article  CAS  Google Scholar 

  65. Graca I, Mendes VM, Marques I, Duro N, da Costa M, Ramalho JC, Pawlowski K, Manadas B, Pinto Ricardo CP, Ribeiro-Barros AI. Comparative proteomic analysis of nodulated and non-nodulated Casuarina glauca Sieb. Ex Spreng. Grown under salinity conditions using sequential window acquisition of all theoretical mass spectra (SWATH-MS). Int J Mol Sci. 2020;21(1):78.

    Article  CAS  Google Scholar 

  66. Wang Y, Dong F, Tang M. Transcriptome analysis of arbuscular mycorrhizal Casuarina glauca in damage mitigation of roots on NaCl stress. Microorganisms. 2022;10(1):15.

    Article  Google Scholar 

  67. Sun S-C, Xiong X-P, Zhang X-L, Feng H-J, Zhu Q-H, Sun J, Li Y-J. Characterization of the Gh4CL gene family reveals a role of Gh4CL7 in drought tolerance. BMC Plant Biol. 2020;20(1):125.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Vogt T. Phenylpropanoid Biosynthesis. Mol Plant. 2010;3(1):2–20.

    Article  CAS  PubMed  Google Scholar 

  69. Andre CM, Schafleitner R, Legay S, Lefevre I, Alvarado Aliaga CA, Nomberto G, Hoffmann L, Hausman J-F, Larondelle Y, Evers D. Gene expression changes related to the production of phenolic compounds in potato tubers grown under drought stress. Phytochemistry. 2009;70(9):1107–16.

    Article  CAS  PubMed  Google Scholar 

  70. Kanehisa M, Goto SKEGG. Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Liping Wei and Ke Liu for technical assistance.

Funding

This study was supported the National Key R&D Program of China (2022YFC3103700, 2021YFC3100401), the Forestry Science and Technology Innovation Fund Project of Guangdong Province (2022KYXM09), and the Youth Innovation Promotion Association of the Chinese Academy of Sciences (Y202077).

Author information

Authors and Affiliations

  1. CAS Engineering Laboratory for Vegetation Ecosystem Restoration on Islands and Coastal Zones, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China

    Shike Zhang, Chunmei He, Shuguang Jian & Nan Liu

  2. College of Life Sciences, Gannan Normal University, Ganzhou, 341000, China

    Shike Zhang & Nan Liu

  3. University of Chinese Academy of Sciences, Beijing, 100049, China

    Shike Zhang

  4. Guangdong Provincial Key Laboratory of Silviculture, Protection and Utilization, Guangdong Academy of Forestry, Guangzhou, 510520, China

    Long Wei

Authors
  1. Shike Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chunmei He
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Long Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shuguang Jian
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nan Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Shuguang Jian, Nan Liu, Long Wei, and Chunmei He conceived and designed the project. Shike Zhang and Nan Liu performed the experiments. Shike Zhang and Chunmei He analyzed data. Shike Zhang wrote the manuscript. Nan Liu and Shuguang Jian supervised and revised the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Shuguang Jian or Nan Liu.

Ethics declarations

Ethics approval and consent to participate

C. equisetifolia seedlings were collected from South China Botanical Garden and were identified by professor Nan Liu. The methods involved in this study were carried out in compliance with local and national regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, S., He, C., Wei, L. et al. Transcriptome and metabolome analysis reveals key genes and secondary metabolites of Casuarina equisetifolia ssp. incana in response to drought stress. BMC Plant Biol 23, 200 (2023). https://doi.org/10.1186/s12870-023-04206-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04206-x

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us