Skip to main content

Transcriptional and physiological data reveal the dehydration memory behavior in switchgrass (Panicum virgatum L.)

Abstract

Background

Switchgrass (Panicum virgatum L.) is a model biofuel plant because of its high biomass, cellulose-richness, easy degradation to ethanol, and the availability of extensive genomic information. However, a little is currently known about the molecular responses of switchgrass plants to dehydration stress, especially multiple dehydration stresses.

Results

Studies on the transcriptional profiles of 35-day-old tissue culture plants revealed 741 dehydration memory genes. Gene Ontology and pathway analysis showed that these genes were enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant hormone signal transduction. Further analysis of specific pathways combined with physiological data suggested that switchgrass improved its dehydration resistance by changing various aspects of its responses to secondary dehydration stress (D2), including the regulation of abscisic acid (ABA) and jasmonic acid (JA) biosynthesis and signal transduction, the biosynthesis of osmolytes (l-proline, stachyose and trehalose), energy metabolism (i.e., metabolic process relating to photosynthetic systems, glycolysis, and the TCA cycle), and lignin biosynthesis. The transcriptional data and chemical substance assays showed that ABA was significantly accumulated during both primary (D1) and secondary (D2) dehydration stresses, whereas JA accumulated during D1 but became significantly less abundant during D2. This suggests the existence of a complicated signaling network of plant hormones in response to repeated dehydration stresses. A homology analysis focusing on switchgrass, maize, and Arabidopsis revealed the conservation and species-specific distribution of dehydration memory genes.

Conclusions

The molecular responses of switchgrass plants to successive dehydration stresses have been systematically characterized, revealing a previously unknown transcriptional memory behavior. These results provide new insights into the mechanisms of dehydration stress responses in plants. The genes and pathways identified in this study will be useful for the genetic improvement of switchgrass and other crops.

Background

Switchgrass (Panicum virgatum L.) is a warm season, cross-pollinated, perennial C4 grass that is one of the most important lignocellulosic biofuel feedstock and pastures [1,2,3]. It belongs to the family Gramineae, and has two ecotypes: a lowland ecotype (tetraploid) and an upland ecotype (hexaploids and octoploids) [4]. Because of its rapid biomass production, cellulose-richness, easy degradation to ethanol, and the availability of extensive genomic information, switchgrass is a model bioenergy plant for research and is widely used for bioethanol production [5,6,7,8,9].

The perennial switchgrass has robust roots, exuberant tillers, and excellent resistance to drought, cold, heat, and biotic stresses [3, 10]. However, it is propagated by seed in large-scale plantations, and the seedlings are sensitive to water-deficit stress because of their slow development and weak roots, especially in arid, salinized or desertified regions. To improve the stress resistance of seedlings and create new varieties, it is important to understand the physiological responses and molecular mechanisms that are activated in switchgrass seedlings exposed to multiple drought stresses.

Previous studies on the drought stress responses of switchgrass seedlings have focused on morphological and physiological indices; a little is known about the response mechanisms at the molecular level. The water potential, gas exchange rates and photochemical traits of switchgrass leaves declined significantly during drought stress, and their drought tolerance correlated with an increase in the levels of reactive oxygen species induced by abscisic acid (ABA) [11, 12]. The drought tolerances of 49 switchgrass genotypes have been assessed using physiological and morphological parameters, revealing that drought-tolerant genotypes tend to have higher levels of ABA, spermine, trehalose, and fructose [13]. In addition, microRNAs have been shown to play important roles in switchgrass responses to drought stress, which involved the biosynthesis of carbon compounds, glucose, starch, fatty acid and lignin [14]. In particular, miR156 miR162 exhibited pronounced changes in expression under severe drought stress, suggesting that they may contribute to the acclimatization process [15].

The dissection of dehydration stress response mechanisms at the molecular level is required to support genetic improvement and germplasm creation by subtly regulating pathways involved in the responses of interest. Numerous studies have shown that ABA plays key roles in plant drought stress responses. Two classes of drought-related signaling pathways have been identified: ABA-dependent and ABA-independent [16, 17]. ABA is synthesized from mevalonic acid in plants; the key enzyme in the process is 9-cis-epoxycarotenoid dioxygenase (NCED) [18]. ABA is degraded by abscisic acid 8′-hydroxylase (CYP707A3) [19] and ABA glucosyltransferase (ABA GTase) [20]. Proteins involved in ABA-dependent signaling pathways include protein phosphatase 2C (PP2C), ABA-responsive element binding factor (AREB/ABF), NAC transcriptional factor (RD26), and the downstream elements RD22 (responsive to dehydration 22), RD29B, RD20A [17, 21,22,23]. ABA-independent signaling proteins include members of the HD-ZIP and DREB2 families, as well as the downstream components early responsive to dehydration 1 (ERD1) and RD29A, among others [17, 24,25,26]. Most modern efforts to identify genetic modifications that improve plant drought resistance involve these two pathways.

However, there is another set of genes that respond to multiple dehydration stress but which have been neglected for a long time. These genes were named dehydration stress memory genes (DMGs) because they exhibit a “memory” of previous dehydration stresses [27, 28]. The expression of DMGs can be strongly and rapidly induced or suppressed when plants experience secondary dehydration [28]. Stress memory has been observed in many species, including Arabidopsis thaliana [28,29,30], Zea mays [31], Brassica rapa [32], Solanum tuberosum [33] and Aptenia cordifolia [34, 35]. These transcriptional memory genes are closely intertwined with histone methylation, plant hormone signaling, and osmotic adjustment [28, 34, 36,37,38,39,40]. In Arabidopsis and maize, four types of DMGs were identified using high-throughput RNA sequencing (RNA-Seq) [30, 31]. DMG expression allows plants to strongly and rapidly respond to repeated or long-term intermittent dehydration stresses and improves their survival rate and biomass [41].

To identify DMGs, explore their formation networks, and establish a firmer foundation for genetic improvement, we analyzed the response patterns of switchgrass plants subjected to repeated dehydration stresses. Switchgrass is a cross-pollinated plant, so tissue-cultured plants were used in this work to eliminate the effects of seedling genotypes. Transcriptomic data were acquired using RNA-Seq, and DMGs and regulatory networks were identified and analyzed carefully based on annotations and published papers. Physiological indices, chemical substance assays and transcriptional information suggested that switchgrass improved its response to secondary dehydration stress by changing the pathways of ABA and jasmonic acid (JA) biosynthesis and signal transduction, increasing the biosynthesis of osmolytes (l-proline, raffinose, stachyose and trehalose), and modifying its energy metabolism systems (specifically, those involved in photosynthesis, glycolysis, and the tricarboxylic acid cycle). Lignin biosynthesis, the persistence of the transcriptional memory, and the conservation and species-specificity of DMGs in switchgrass, maize and Arabidopsis were also analyzed.

Methods

Plant material and experimental design

Switchgrass is allogamous and has strong genetic self-incompatibility, so it is difficult to create homogenous materials. In this study, we acquired a homogeneous genotype of switchgrass by in vitro tissue culture from axillary buds of a single plant belonging to the lowland Alamo cultivar (introduced from the USA and domesticated at Northwest A&F University, China). The regenerated shoots were proliferated in MS medium supplemented with 13.3 μM L−1 6-benzylaminopurine (BAP) [42]. Cluster shoots were separated and transplanted to substrates of sand and soil (1:1) in a greenhouse illuminated with natural sunlight (day/night ~ 14/10 h; daily temperature 28–32 °C). Multiple dehydration stresses were imposed based on designs previously used with Arabidopsis [28, 30] and maize [31] (Fig. 1a). Five-week-old homogenous plants with three fully expanded leaves were removed from their substrates, washed to remove residual soil, and then acclimated overnight in a plant growth chamber with their roots in water. The next morning, the plants were blotted onto filter paper to remove water, and air-dried for 2 h (this corresponded to the first dehydration stress, D1). At the same time, non-stressed control plants designated C1 were sampled as controls. The drying was followed by a 22 h recovery period of full re-hydration at 25 °C (first recovery, R1). For the subsequent stress treatment, R1 plants were air-dried for 2 h (the second dehydration stress, D2) followed by another recovery period (second recovery, R2). The same procedures were repeated for D3, R3, and so on (Fig. 1a). To determine the duration of the transcriptional memory, stressed plants were replanted in soil for 3, 5, or 7 days before being re-exposed to dehydration stress [28] (Fig. 1b).

Fig. 1
figure 1

Experimental design and physiological assays. a Experimental design for repeated dehydration stresses and sample collection. b Experimental design for determining the persistence of dehydration memory and sample collection. c Survival rates of trained and non-trained switchgrass. d Relative water content of leaves during repeated dehydration stresses. e Water loss from leaves during the first, second and third dehydration stresses. fi, Determination of photosynthetic activity during repeated dehydration stresses. C control, D13 the first, second, and third dehydration stresses, R12 the first and second recovery periods, RNA-Seq RNA sequencing, qRT-PCR quantitative real-time PCR, RWC relative water content, PS I and PS II photosynthesis system I and II, ABA abscisic acid, JA jasmonic acid, Fv/Fm maximum photosynthetic potential, qN non-photochemical quenching, Yield (I) and Yield (II) real-time quantum yield of PSI and PSII. Values for RWC, water loss and photosynthesis parameters are means of n = 5 measurements ± SD

RNA isolation and RNA-Seq library construction

Second fully expanded leaves were collected for RNA isolation and for physiological and biochemical assays. The collected leaves were immediately frozen in liquid nitrogen and stored at − 80 °C. Total RNA was extracted with the Trizol reagent (Invitrogen Life Technologies, USA) according to a published protocol [43]. The isolated RNA was subsequently treated with RNase-Free DNase I (Takara, http://www.takarabiomed.com.cn/). The quantity, purity and integrity of the resulting RNA was confirmed using NanoDrop One (Thermo Scientific, USA) and Agilent 2100 (Agilent Technologies, USA) instruments; only RNAs classified as A-type (optimal RNA samples for RNA sequencing, BGI-Shenzhen, http://www.genomics.cn/) were used for cDNA library construction.

Two biological replicates of four treatments were sequenced (C1, D1, R1 and D2, Fig. 1a), giving a total of eight RNA sequencing (RNA-Seq) samples. Each treatment contained more than 12 plants. Total RNA was performed according to the standard transcriptome resequencing experimental process, and was sequenced using the Illumina GAII platform (HiSeq2500) at BGI-Shenzhen [43]. Expressed tags were calculated and normalized to FPKM (Fragments Per Kilobase of exon per Million fragments mapped) values for gene expression analysis using the Cufflinks program [44]. An overview of the RNA-Seq process and data on the repeatability of different biological replicates are provided in Additional file 1: Figure S1. The raw sequencing files of transcriptomic data were uploaded to sequence read archive (SRA) at NCBI (Accession Number PRJNA394165).

Bioinformatics analysis

Genes in the transcriptomic profiles were predicted using data from the switchgrass genomic database (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Pvirgatum) [8]. Differentially expressed genes (DEGs) were identified using the Cuffdiff [44] tool from the Cufflinks package (http://cole-trapnell-lab.github.io/cufflinks/) based on their FPKM values and the associated q values (false discovery rate adjusted p values) between different treatments. Specifically, a gene was identified as a DEG if it had a q value ≤ 0.05 and a | log (base2) fold change | ≥ 1. In addition, the gene’s FPKM value under the dehydration stress (for up-regulated genes) or control (for down-regulated genes) treatments was required to be at least 20% of the mean FPKM value for all expressed genes [45]. Genes exhibiting common or differential expression in different treatments were compared and identified using Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/index.html).

Dehydration stress memory genes (DMGs) were identified by comparing expressed fold changes between D1/C1 and D2/D1 within the DEG libraries. Genes that were up-regulated, down-regulated, and exhibited no significant change in expression are indicated by the symbols “+”, “−” and “=”, respectively. Changes in an individual gene’s expression are denoted in the format [a/b], where a is the symbol indicating the change between C1 and D1, and b is the symbol indicating the change between D1 and D2. Thus, for example, a gene of the [+/+] type was up-regulated between C1 and D1, and also between D1 and D2, whereas a gene of the [=/+] type exhibited no change in expression between C1 and D1, but was up-regulated between D1 and D2 [30, 31]. Eight types of DEGs were defined [+/+], [+/−], [−/+], [−/−], [+/=], [−/=], [=/+] and [=/−] (Table 1). The first four types were considered to be DMGs because their responses changed between D1 and D2. The [+/=] and [−/=] genes were considered non-memory genes, and [=/+] and [=/−] genes were defined as late-response genes. A list of all identified DMGs is provided in Additional file 2: Table S1.

Table 1 Dehydration response and transcriptional memory genes in switchgrass, maize and Arabidopsis

All the identified switchgrass genes were annotated using the NCBI Blast tools (version 2.2.8) in conjunction with the NCBI non-redundant protein sequences (nr) database, the Arabidopsis Information Resource Proteins database (release TAIR10), the maize B73 RefGen_v3 database, and the Rice Genome Annotation Project (RGAP) release 7. Gene Ontology (GO) analysis was performed using AgBase GORetriever and GOSlimViewer [46], and pathways were enriched using the NCBI Flink (https://www.ncbi.nlm.nih.gov/Structure/flink/flink.cgi) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. Specific pathways were verified using Blastp in conjunction with the switchgrass genomic database and the TAIR10 database. Orthologous genes belonging to larger families were analyzed with reference to previously published works.

Reverse transcription and quantitative real-time PCR

Total RNA was reverse-transcribed into cDNA using the PrimeScript™ RT reagent Kit with the gDNA Eraser (RR047A, Takara) according to the standard protocol. cDNA was detected using two housekeeping genes, PvACT2 (PvActin 2) [47] and PveEF- (Eukaryotic elongation factor 1α) [48]. Quantitative real-time PCR (qRT-PCR) was performed with QuantStudio™ 7 Flex Real-Time PCR System (Thermo Scientific, USA) using SYBR Premix Ex Taq™ II Kit (RR820A, Takara) with three technical replicates. Relative gene expressed levels were calculated using the ΔΔCT method with PveEF- as the internal control, and three biological replicates were performed for each experiment. The primers used in the qRT-PCR experiments were designed using Primer Premier 6 (http://www.premierbiosoft.com/primerdesign/index.html), and are listed in Additional file 3: Table S2.

Physiological indices and chemical substance assays

The plants used to determine the survival rates of stress trained or non-stress-trained plants were grown in substrates of sand and soil (1:1) in greenhouses as described in the plant material section. Relative water content (RWC) was calculated using the equation \({\text{RWC }}\left( \% \right) = [({\text{fresh}}\,{\text{weight}} - {\text{dry}}\,{\text{weight}})/({\text{turgid}}\,{\text{weight}}\, -\, {\text{dry}}\,{\text{weight}})] \times 100\%\) [28]. The water loss was determined by measuring the leaves’ weight at a fixed time interval after their detachment from the plants [28]. The Photosynthesis System I and II (PSI and PSII) parameters were measured with a DUAL-PAM-100 fluorometer (WALZ, Germany) by synchronously recording changes in P700 (PSI) and chlorophyll fluorescence (PSII) [49]. The contents of abscisic acid (ABA) and jasmonic acid (JA) in the leaves of plants subjected to repeated dehydration stress were determined by high performance liquid chromatography–mass spectrometry (HPLC–MS) with an Agilent 1290 Infinity II liquid chromatograph (Agilent Technologies, USA) and a QTRAP 6500 MS/MS System (AB SCIEX, USA). Reference standards (ABA and JA, LC grade) were bought from Sigma-Aldrich (USA). Proline was extracted from leaves with 3% sulfosalicylic acid, and then reacted with acid ninhydrin reagent. After extracting with methylbenzene, absorbances at 520 nm were measured to calculate proline content.

Five biological replicates were used in the determination of RWC, water loss and photosynthesis parameters, and three biological replicates were used for hormone determination. Each biological replicate involved at least 12 plants, and only the second fully expanded leaves were collected for the experiments. All data were analyzed with the Duncan test [50] using the SAS 9.2 software (http://www.sas.com/en_us/software/sas9.html).

Results

Physiological indices indicate that switchgrass exhibits dehydration memory

To determine whether switchgrass plants exhibit stress memory after a previous dehydration stress, we designed multiple dehydration stress experiments and measured the plants’ survival rates, leaf relative water content (RWC) and water loss as the experiments progressed (Fig. 1a). The survival rates of trained plants were significantly higher than those for non-trained plants (Fig. 1c), and the plants that had experienced dehydration stress also exhibited slower water loss than non-trained plants (Fig. 1e). The RWC of the leaves after the first recovery (R1) and second recovery (R2) periods was similar to that in control plants (Fig. 1d), indicating that the plants recovered fully from the previous dehydration stress during the recovery periods, and started each new stress period with a similar RWC. In addition, there were significant changes in the plants’ survival rates and water losses during successive dehydration periods, indicating that the training had improved the plants’ dehydration stress resistance and that the trained plants exhibited stress memory.

To validate this conclusion, selected photosynthetic parameters were measured for PSI and PSII, namely Fv/Fm, non-photochemical quenching (qN), Yield (I) and Yield (II). These parameters reflect the maximum photosynthetic potential of the leaves, the photoprotective capacity of the leaves, and the real-time quantum yields of PSI and PSII, respectively [51]. There were no appreciable differences in Fv/Fm or qN between plants during the first dehydration stress (D1) and the corresponding controls (C1), or between plants during the first recovery period (R1) and the second dehydration stress (D2) (Fig. 1f, g). However, Yield (I) and Yield (II) during R1 and D2 were significantly lower than in C1 and D1. These results indicated that switchgrass underwent appreciable changes after the first dehydration stress even though the mild stresses imposed in this work were not lethal.

Identification of dehydration memory genes (DMGs) in switchgrass

To characterize the mechanisms responsible for dehydration memory in switchgrass, the transcriptomes of C1, D1, R1, and D2 plants were studied by high-throughput RNA sequencing (Fig. 1a). 98,007 genes have been identified according to switchgrass genome, and 47,207 genes (48.17%) were expressed in the transcriptomes. The average proportion of clean reads in all samples was 97.14%, and the repetitive rates of gene expression numbers across two biological replicates were 89.73% (C1), 87.77% (D1), 92.40% (R1) and 91.08% (D2) (Additional file 1: Figure S1). Differentially expressed genes (DEGs) were identified as genes for which q ≤ 0.05, | log (base2) fold change | ≥ 1, and the FPKM value under dehydration stress (for up-regulated genes) or in the control treatment (for down-regulated genes) was at least 20% of the mean FPKM value for all expressed genes [45]. The coefficient values of DEGs in two biological replicates were 98.66% (C1), 99.23% (D1), 99.78% (R1) and 98.67% (D2) (Additional file 1: Figure S1), demonstrating the excellent repeatability of the results across biological replicates.

Analysis of the expressed fold changes for DEGs between D1 and C1, and between D2 and D1, revealed four types of DMGs, two types of non-memory genes, and two types of late-response genes in switchgrass. Switchgrass has more detected genes than Arabidopsis or maize, but its number of dehydration response genes (1566) is lower than that determined for either species (2062 in maize, and 6759 in Arabidopsis) (Table 1). The number of transcriptional memory genes in switchgrass (741) was similar to that in maize (816), but significantly lower than that in Arabidopsis (1963). The non-memory genes exhibited similar trends in all three species. Interestingly, switchgrass has twice as many late-response genes as maize, and three times as many as Arabidopsis (Table 1).

The memory types of switchgrass dehydration stress response genes were revealed by constructing Venn diagrams (Fig. 2). There were 811 and 3679 genes up-regulated in D1 and D2, respectively, of which 356 were common genes (i.e., genes up-regulated in both dehydration periods). Ninety-three of the memory genes were [+/+] memory genes, corresponding to 26.12% of the common up-regulated genes and 2.2% of all up-regulated genes (Fig. 2a). In addition, there were 755 and 3534 genes that were down-regulated in D1 and D2, respectively, and 517 common genes that were down-regulated in both dehydration periods. One hundred and thirteen of these genes were [−/−] memory genes, corresponding to 21.86% of the common down-regulated genes and 3.0% of all down-regulated genes (Fig. 2b). Moreover, 381 [+/−] memory genes and 154 [−/+] memory genes were identified, corresponding to 5 and 2% of all DEGs (7205), respectively (Fig. 2c). Finally, 337 (4.4%) [+/=] type and 488 (6.4%) [−/=] type non-memory genes were identified, along with 3158 (41.6%) [=/+] type and 2860 (37.7%) [=/−] type late-response genes (Fig. 2c). The numbers of [=/+] and [=/−] late-response genes identified in switchgrass were both significantly higher than the corresponding numbers for maize and Arabidopsis (Table 1). Among the dehydration memory genes, [+/−] memory genes were significantly more abundant than the other three types, accounting for 51% of all DMGs in switchgrass, 65% in maize, and 44% in Arabidopsis (Table 1).

Fig. 2
figure 2

a Venn diagrams of up-regulated genes in D1 and D2. b Venn diagrams of down-regulated genes in D1 and D2. c Venn diagrams of different types of dehydration memory genes. C control, D12 first and second dehydration stresses, R1 first recovery period, “+”, “−” and “=” denote up-regulated genes, down-regulated genes, and genes with no significant changes in expression, respectively

Verification of DMGs by quantitative real-time PCR (qRT-PCR)

The RNA-Seq results were validated using qRT-PCR to quantify the expression of two randomly selected genes from each DMG types, making a total of 8 genes or 32 treatments: Pavir.J32219 and Pavir.Aa02860 representing the [+/+] type, Pavir.J31344 and Pavir.Ea03187 representing the [−/−] type, Pavir.Ea00263 and Pavir.J23463 representing the [+/−] type, and Pavir.Ea00667 and Pavir.J10667 representing the [−/+] type. The measured expression of all eight genes was entirely consistent with the transcriptomic data (Fig. 3), indicating that the transcriptomic approach based on RNA-Seq is a reliable method for analyzing the response of switchgrass to multiple dehydration stresses.

Fig. 3
figure 3

Verification of four types of dehydration memory genes by quantitative real-time PCR. C control, D12 first and second dehydration stresses, R1 first recovery period. Relative gene expressed levels were calculated using the ΔΔCT method with PveEF- as the internal control, and three biological replicates were performed for each experiment

Functional analysis of switchgrass DMGs

To explore the underlying mechanism of dehydration memory, all the DMGs were annotated using Blastp in conjunction with the NCBI (nr) protein database and species-specific databases for Arabidopsis (TAIR10), maize (B73 RefGen_v3), and rice (RGAP release 7) (Additional file 2: Table S1). The Gene Ontology (GO) annotations of DEGs in D1 and D2, and DMGs, were collected and used to construct graphs (Fig. 4). The chosen genes were enriched in the following Cellular Component GO terms: membrane (GO: 0016020), intracellular (GO: 0005634) and nucleus (GO: 0005622). They were also enriched in the binding (GO: 0005488), catalytic activity (GO: 0003824) and transferase activity (GO: 0016740) terms from the molecular function category (Fig. 4). The biological process annotations revealed enrichment in the terms biosynthetic process (GO: 0009058), cellular protein modification process (GO: 0006464) and nucleobase-containing compound metabolic process (GO: 0006139) in D1, while GO: 0009058, GO: 0006139 and transport (GO: 0006810) were enriched in D2. In addition to these four terms, the carbohydrate metabolic process (GO: 0005975), response to stress (GO: 0006950) and signal transduction (GO: 0007165) terms were also enriched among the DMGs (Fig. 4). Finally, the pathway enrichment analysis showed that phenylpropanoid biosynthesis (ko00940), starch and sucrose metabolism (ko00500), plant hormone signal transduction (ko04075), photosynthesis—antenna proteins (ko00196), galactose metabolism (ko00052) were enriched among the DMGs (Additional file 4: Table S3). These results indicated that dehydration stress memory in switchgrass involves plant hormone biosynthesis and signal transduction, osmotic adjustment (ko00500, ko00052, etc.), photosynthesis, and lignin biosynthesis (ko00940). Therefore, differentially expressed genes associated with these pathways were investigated.

Fig. 4
figure 4

Gene Ontology analysis of dehydration memory genes with differentially expressed genes in D1 and D2. a Cellular component. b Molecular function. c Biological process. D1–2, first and second dehydration stresses

Plant hormone biosynthesis and signal transduction

Abscisic acid (ABA) biosynthesis and signal transduction

Plant hormone biosynthesis and signal transduction (GO: 0007165; ko04075) were enriched in both the GO and pathway analyses (Fig. 4; Additional file 4: Table S3), and more than 80 DEGs were identified as being involved in ABA and jasmonic acid (JA) biosynthesis and signal transduction (Additional file 5: Table S4). ABA plays important roles in the drought stress response, contributing to signal perception and transduction, stomatal closure, and osmotic adjustment [20, 38, 52]. In switchgrass subjected to multiple dehydration stresses, four orthologs of AtNCEDs, the rate-limiting enzyme of ABA biosynthesis, were significantly up-regulated in D1 and D2. One of them (Pavir.Ba03791) was a [+/+]-type memory gene. In addition, two orthologs of CYP707A3 (Pavir.Ab02775 and Pavir.J08122), a key enzyme of ABA catabolism, were down-regulated in D1 and slightly down-regulated in D2 (Fig. 5a; Additional file 5: Table S4). The qRT-PCR relative expression analysis also showed that the orthologs of NCED (Pavir.J24772 and Pavir.Ba03791) were strongly up-regulated in D1 and D2 (Fig. 5c). The transcriptional data suggest that ABA accumulated during the dehydration stresses. Therefore, we determined the ABA contents of leaves in different treatments using HPLC–MS, confirming that this hormone accumulated significantly in D1, R1 and D2, with very large accumulations occurring in the latter two periods (Fig. 5d). These results strongly suggest that ABA plays a key role in the response of switchgrass to repeated dehydration stresses. Moreover, the qRT-PCR and HPLC–MS results verified the reliability of the initial transcriptional data.

Fig. 5
figure 5

ABA and JA metabolism in multiple dehydration stresses. a ABA biosynthesis and catabolism. b JA biosynthesis. c Relative expression of orthologs of NCED and ABF. d ABA and JA concentrations in leaves determined by HPLC-MS. ABA abscisic acid, JA jasmonic acid, NSY neoxanthin synthase [EC 5.3.99.9], NCED 9-cis-epoxycarotenoid dioxygenase [EC 1.13.11.51], CYP707A [+]-abscisic acid 8′-hydroxylase [EC 1.14.13.93], ABF ABA responsive element binding factor, 13-LOX linoleate 13S-lipoxygenase [EC 1.13.11.12], AOS allene oxide synthase [EC 4.2.1.92], AOC allene oxide cyclase [EC 5.3.99.6], OPR 12-oxophytodienoate reductase [EC 1.3.1.42], ACX acyl-CoA oxidase [EC 1.3.3.6], MFP multifunctional—protein enoyl-CoA hydratase [EC 4.2.1.17], KAT 3-ketoacyl-CoA thiolase [EC 1.1.1.35], ACAT acetyl-CoA C-acyltransferase [EC 2.3.1.16], ACH acyl-CoA hydrolase [EC 3.1.2.20], FW fresh weight. Relative expression values and plant hormone concentrations are means of n = 3 measurements ± SD (Duncan’s test: **p < 0.01)

Clade A protein phosphatase 2C (PP2CA) enzymes function as co-receptors in ABA signaling and respond positively to increases in ABA levels or stress-induced ABA biosynthesis [53, 54]. Among the downstream enzymes of the PP2Cs are subfamily 2 of the SNF1-related kinases (SnRK2s), which are positive regulators of ABA signal transduction [55, 56]. Our results identified nine PP2CA orthologs that were strongly up-regulated in D1 and D2, and four orthologs of SnRK2s that were up-regulated in D1 and D2, especially in D2 (Fig. 6; Additional file 5: Table S4). Among these genes, three PP2CA orthologs and one SnRK2 ortholog were [+/+] type memory genes. In addition, two orthologs (Pavir.Aa00597, Pavir.J00256) of ABA responsive element binding factors (ABFs) belonging to the ABA-dependent pathway were significantly up-regulated in D1 and D2 (Fig. 6; Additional file 5: Table S4). Both the transcriptional and qRT-PCR data indicated that these two genes functioned as non-memory genes ([+/=]) (Fig. 5c). In addition, two orthologs (Pavir.Ea02724, Pavir.Eb03036) of responsive to desiccation 22 (RD22) and one ortholog (Pavir.J02116) of ras-related protein Rab-18 were significantly up-regulated in D2 (Additional file 5: Table S4).

Fig. 6
figure 6

ABA and JA signal transduction in multiple dehydration stresses. ABA abscisic acid, JA jasmonic acid, NCED 9-cis-epoxycarotenoid dioxygenase [EC 1.13.11.51], NAC transcription factor NAC, RD responsive to desiccation, ABF/AREB ABA responsive element binding factor, SnRK2 serine/threonine-protein kinase SRK2 [EC:2.7.11.1], AOS allene oxide synthase [EC 4.2.1.92], AOC allene oxide cyclase [EC 5.3.99.6], JAR1 jasmonic acid-amino synthetase, COI1 coronatine-insensitive protein 1, JAZ jasmonate ZIM domain-containing protein, MYC2 transcription factor MYC2

JA biosynthesis and signal transduction

JA is synthesized in the chloroplasts and peroxisomes, and a previous study demonstrated that it engages in crosstalk with ABA biosynthesis and signal transduction, and is also involved in biotic or abiotic stress responses [40, 57,58,59,60,61]. More than 50 of the DEGs identified in this work were annotated with JA biosynthesis and signal transduction terms (Additional file 6: Table S5). Twelve orthologs of upstream genes that function in the chloroplast (LOX, AOS and AOC) were down-regulated in D2, and 18 orthologs of downstream genes that function in the peroxisome (OPR, ACX, MFP, KAT, ACAT and ACH) were up-regulated in D2 (Fig. 5b; Additional file 6: Table S5). The HPLC–MS assay showed that JA accumulated significantly in D1 but returned to its initial level in R1 and D2 (Fig. 5d).

JAR1, COI1, and the JAZ transcription factors play important roles in JA signaling [58, 61]. In the present study, six orthologs of JAR1 and five orthologs of the JAZ transcription factor were down-regulated during the multiple dehydration stresses, especially in D2, while the COI1 orthologs Pavir.Ea03420 and Pavir.J07109 were strongly up-regulated in D2 (Fig. 6; Additional file 6: Table S5). Furthermore, two orthologs of MYC2, a basic helix-loop-helix (bHLH) transcription factor involved in AOC expression and the crosstalk between ABA and JA biosynthesis and signal transduction [40], were slightly up-regulated in D1 and down-regulated in D2. Another four orthologs of MYC2 were down-regulated in both D1 and D2 (Fig. 6; Additional file 6: Table S5).

Other plant hormones in multiple dehydration stresses

Except ABA and JA, 34 orthologs involved in auxins, brassinosteroids, cytokinins, ethylene, gibberellins and salicylic acid biosynthesis and signal transduction exhibited dehydration stress memory (Additional file 7: Table S6). The identified orthologs of auxins biosynthesis and signal transduction were up-regulated, especially in D2. Three orthologs involved in ethylene biosynthesis (Pavir.J11398, Pavir.Ca01179 and Pavir.Ca01180) were up-regulated in D1 and down-regulated in D2. Brassinosteroids, cytokinins, gibberellins and salicylic acid were also participated in the multiple dehydration stresses (Additional file 7: Table S6).

Transcription factors (TFs)

In addition to the above TFs (ABFs, JAZs) involved in ABA and JA signal transduction, 46 TFs showed dehydration stress memory in switchgrass, including eight AP2-EREBP genes, six bHLH genes, six MYB/MYB-related genes, six WRKY genes, five NAC genes, four C2H2 genes, three C3H genes and eight other TF genes (Table 2). As with DMGs in general, [+/−] type TFs were the most abundant dehydration memory type, accounting for 47.8% of all dehydration memory TFs. Interestingly, all the dehydration memory WRKY genes belonged to the [+/−] memory type (Table 2).

Table 2 Transcription factor involved in transcriptional memory

Osmotic adjustment

Osmotic adjustment is essential for water uptake and maintenance, membrane protection, and ROS scavenging in plants subjected to dehydration stress [62]. Switchgrass plants exposed to repeated dehydration stress exhibited differential expression of genes related to proline, raffinose family oligosaccharides (RFO) and trehalose metabolism (Fig. 7; Additional file 8: Table S7). Eight orthologs of P5CS and OAT, the key enzymes involved in the biosynthesis of l-proline from l-glutamate and l-ornithine, respectively, were up-regulated in D2, and an ortholog of proline dehydrogenase (ProDH), which is involved in proline degradation, was down-regulated in D1 and D2 (Fig. 7a; Additional file 8: Table S7). The proline determination assay showed that the proline content of switchgrass leaves had no significant change in D1, but was accumulated in R1 and D2, especially in D2 (Fig. 7b). Several genes related to RFO biosynthesis were up-regulated in D1 and D2, and expressed significantly more strongly in D2, including the AtGolS3 (AT1G09350) orthologs Pavir.Ba00293 and Pavir.Bb03717, AtRS2 (AT3G57520) ortholog Pavir.J05624, and AtRS5 (AT5G40390) ortholog Pavir.Eb00453. All these genes were [+/+] memory type genes. Another four orthologs of AtRS2 and two orthologs of STS were strongly up-regulated in D2 (Fig. 7c; Additional file 8: Table S7). Among genes related to trehalose biosynthesis, the AtTPS1 (AT1G78580) ortholog Pavir.Ca01680, Pavir.J00828, and AtTPS8 (AT1G70290) ortholog Pavir.Ea02872 were down-regulated in D2, while the AtTPPC (AT1G22210) ortholog Pavir.J12499, AtTPPF (AT4G12430) ortholog Pavir.Ab02625, and Pavir.J38619 were down-regulated in D1 and D2. Additionally, five orthologs of AtTPS5-7 and four orthologs of AtTPPJ were up-regulated in D1 and D2, especially during the second stress. Of these genes, the AtTPS5 (AT4G17770) ortholog Pavir.Ia04010 and AtTPPJ (AT5G65140) ortholog Pavir.Fb01190 were [+/+]-type DMGs (Fig. 7d; Additional file 8: Table S7). These results indicate that l-proline, raffinose, stachyose and trehalose play key roles in the resistance of switchgrass to dehydration stress. Moreover, the fact that several [+/+] type memory genes were related to RFO and trehalose biosynthesis demonstrates that switchgrass plants “remember” their first dehydration stress and improve their stress resistance by making a stronger osmotic adjustment during subsequent dehydration stresses.

Fig. 7
figure 7

Schematic overview of l-proline, raffinose family oligosaccharide, and trehalose biosynthesis. a The biosynthesis of l-proline. b The contents of proline in leaves during repeated dehydration stresses. c The biosynthesis of raffinose family oligosaccharides. d The biosynthesis of trehalose. P5CS glutamate 5-kinase [EC 2.7.2.11], OAT ornithine aminotransferase [EC 2.6.1.13], ProDH proline dehydrogenase [EC 1.5.5.2], FW fresh weight, GOLS inositol 3-α-galactosyltransferase/galactinol synthase [EC 2.4.1.123], RS galactinol–sucrose galactosyltransferase/raffinose synthase [EC 2.4.1.82], STS galactinol–raffinose galactosyltransferase/stachyose synthetase [EC 2.4.1.67], TPS trehalose 6-phosphate synthetase [EC 2.4.1.15], TPP trehalose 6-phosphate phosphatase [EC 3.1.3.12]

Photosynthesis, glycolysis and tricarboxylic acid cycle (TCA cycle)

Physiological assays and pathway enrichment analyses have shown that multiple dehydration stresses cause significant changes in photosynthesis-related energy metabolism and the basic physiological processes that produce ATP and related intermediates that are vital to plant metabolism (Figs. 1f–i; 4; Additional file 4: Table S3). We therefore investigated differentially expressed genes relating to photosynthesis, glycolysis and the TCA cycle. All 28 genes related to the light reactions of photosynthesis (photosystem I, photosystem II, plastoquinol–plastocyanin reductase and ferredoxin-NADP+ reductase) were significantly down-regulated in D2, acting as [=/−] type late-response genes (Additional file 9: Table S8). Among DEGs associated with the C4 photosynthetic carbon assimilation cycle, 15 (78.9%) were down-regulated in D2, and 38 (79.2%) genes related to Calvin cycle were repressed in D2 (Additional file 9: Table S8). These transcriptional data were consistent with the results of physiological PSI and PSII assays using a Dual-PAM 100 instrument (Fig. 1h, i), which indicated that the photosynthetic systems of the switchgrass plants were repressed in D2 and that their behavior during D2 differed significantly from that in D1.

Glycolysis and TCA cycle provide much of the ATP and related intermediates used in plant metabolism. In switchgrass subjected to dehydration stress, nine (90%) orthologs of three key enzymes in glycolysis—hexokinase (EC 2.7.1.1), 6-phosphofructokinase (EC 2.7.1.11) and pyruvate kinase (EC 2.7.1.40)—were up-regulated in D2, indicating that the rate of glycolysis probably increased during the secondary stress even though several other genes related to glycolysis were down-regulated in D2 (Additional file 10: Table S9). In addition, three orthologs of enzymes involved in the first two steps in TCA cycle, citrate (Si)-synthase (EC 2.3.3.1) and aconitate hydratase (EC 4.2.1.3), were up-regulated in D2, while three orthologs of malate dehydrogenase (EC 1.1.1.37) were down-regulated in D2 (Additional file 10: Table S9). These changes in the expression of genes related to glycolysis and the TCA cycle demonstrated that switchgrass plants improved their energy metabolism upon repeated dehydration stress.

Changes in lignin biosynthesis induced by multiple dehydration stresses

Switchgrass is a target lignocellulosic biofuel species, and water-deficit stress strongly affects its biomass accumulation [13, 14]. Pathway enrichment analysis of DMGs in switchgrass plants subjected to repeated dehydration stress revealed enrichment in phenylpropanoid biosynthesis (ko00940) (Additional file 4: Table S3). Based on the published pathway of switchgrass lignin biosynthesis [63,64,65,66,67] and the KEGG database (map00940), we investigated the expression of related genes during the multiple stress treatment. Of the 28 such genes that were detected, ten were differentially expressed: PvPAL1 and PvCCoAOMT2b were down-regulated in D1 and up-regulated in D2, belonging to the [−/+] memory gene type; eight genes were identified as [=/+] late-response genes, including PvC4H2, PvC3′H1a, PvCCoAOMT1a, PvHCT1a, PvHCT2a, PvF5H1, PvCCR-like1a and PvCAD1 (Fig. 8; Additional file 11: Table S10). These results demonstrate that lignin biosynthesis in switchgrass was repressed in D1 but increased in D2.

Fig. 8
figure 8

Lignin biosynthesis in multiple dehydration stresses. This model was established based on data from the KEGG database (map00940) and published papers (Shen et al. [63]; Xu et al. [64]; Escamilla-Treviño et al. [65]; Fu et al. [66, 67]). Thicker lines indicate the routes that are most strongly supported by both biochemical and genetic evidence. PAL phenylalanine ammonia-lyase, 4CL 4-coumarate CoA ligase, C4H cinnamate 4-hydroxylase, C3′H 4-coumaroyl shikimate 3′-hydroxylase, CCoAOMT caffeoyl-CoA 3-O-methyltransferase, HCT hydroxycinnamoyl CoA:shikimate hydroxycinnamoyl transferase, F5H ferulate 5-hydroxylase, CCR cinnamoyl CoA reductase, COMT caffeic acid 3-O-methyltransferase, CAD cinnamyl alcohol dehydrogenase

Persistence of transcriptional memory

To investigate the persistence of the transcriptional memory in the absence of dehydration stress, trained and untrained plants were replanted in soil for 3, 5, and 7 days (Fig. 1b), and the expression of four representative memory genes and two non-memory genes was monitored by qRT-PCR. After 3 and 5 days, the expressions of the memory genes in trained plants subjected to a second dehydration stress event (D2) differed significantly from those of non-trained plants subjected to their first dehydration stress (D1). However, the expression of these genes in trained plants started to return towards the baseline level after 7 days. The relative expression of the non-memory genes in D2 was always similar to that in D1 (Fig. 9). These results showed that dehydration stress memory persists for around 5 days in switchgrass plants, but gradually fades after 7 days in water.

Fig. 9
figure 9

Length of the transcriptional memory detected by quantitative real-time PCR. C control, D1–2 the first and second dehydration stresses. Relative gene expressed levels were calculated using the ΔΔCT method with PveEF- as the internal control, and three biological replicates were performed for each experiment

Dehydration memory genes of switchgrass and maize

Dehydration memory has been observed in maize and Arabidopsis, and is considered to be a common phenomenon in plants [30, 31]. To better characterize the conservation and species-specificity of dehydration memory genes, we investigated the similarities and differences between dehydration memory in switchgrass, maize and Arabidopsis. In total, 47,207 genes were detected in the switchgrass transcriptomes, which was more than in maize (39,635) or Arabidopsis (33,555). There were 1566 (3.3%) and 5639 (11.9%) switchgrass genes that were differentially expressed in D1 and D2; the corresponding values for maize were 2062 (5.2%) and 2924 (7.4%), while those for Arabidopsis were 6579 (19.6%) and 1371 (4.1%) (Table 1). The number of response genes in switchgrass was similar to that in maize, but significantly less than in Arabidopsis. In switchgrass, 47.3% (741) of all dehydration response genes were identified as DMGs, while the ratio was 39.6% (816) in maize and 29.8% (1963) in Arabidopsis. Remarkably, the 51% of DMGs were [+/−] type memory genes, which is similar to the proportions in maize (65%) and Arabidopsis (44%) (Table 1).

To clarify the functions of DMGs and identify key genes involved in dehydration memory, we compared switchgrass DMGs to maize and Arabidopsis databases using the Blast tools. Few DMGs in switchgrass were homologous to proteins in the TAIR10 Arabidopsis database, so the following section focuses on similarities and differences between switchgrass and maize. 741 memory genes had 537 orthologs in maize, including 76 (14.2%) DMGs, 99 (18.4%) non-memory genes, 77 (14.3%) late-responding genes and 285 non-responding genes. Among these 76 DMGs, 64 (84.2%) were [+/−] type genes (Fig. 10a).

Fig. 10
figure 10

Orthologous compare of response patterns to multiple dehydration stress in switchgrass and maize. a Comparison of orthologous dehydration memory genes in switchgrass and maize. b Distribution of switchgrass orthologs in maize according to their memory response types

We then investigated whether the orthologs of memory genes in switchgrass displayed similar types of memory responses in maize. The maize orthologs of the switchgrass memory genes included both conserved and non-conserved memory types, but the orthologs belonged to other types (non-memory genes, late-responding genes and non-responding genes) were the largest categories (Fig. 10b). Every type of switchgrass memory genes had orthologs belonging to the same memory type in maize, even though the numbers were significantly different: 3 [+/+], 31 [+/−], 1 [−/+] and 2 [−/−]. Remarkably, almost all the [+/−] type orthologs (93.9%) belonged to the same memory type (Fig. 10b; Additional file 12: Table S11). These conserved memory genes probably played important roles in the multiple stress response and are involved in plant growth, signal transduction, and other unknown pathways (30.2%). In addition, many memory orthologs belonged to different memory types (48.7%) in the two species, demonstrating the evolution of species-specific transcriptional memory responses (Additional file 13: Table S12).

Discussion

Biological significance of transcription memory in switchgrass

In natural environments, plants may suffer water-deficit stress on several occasions over a few days [68]. In this study, we designed experiments to alter the water balance with successive cycles of water-deficit stresses. We hypothesized that dehydration stress memory in switchgrass helped the trained plants to resist subsequent water-deficit stresses with stronger and faster dehydration responses.

Several gene families were found to include coexisting dehydration stress memory genes, non-memory genes, and late-response genes. This has also been observed in both Arabidopsis and maize [30, 31]. For instance, nine orthologs of NCED were detected by RNA-Seq, four of which were differentially expressed in D1 and D2. However, only one (Pavir.Ba03791) was a [+/+] memory gene, i.e., a gene expressed more strongly in D2 than in D1 (Additional file 5: Table S4). Trehalose 6-phosphate synthetase (TPS) is the key enzyme in trehalose biosynthesis [69]. Of its eight differentially expressed orthologs in switchgrass, three were down-regulated ([=/−]) between D1 and D2, four were up-regulated (=/+) between D1 and D2, and one was a [+/+] memory gene (Additional file 8: Table S7). These results suggest that the orthologs have developed different functions during evolutionary history, with several members of the family providing stress memory to reduce energy consumption and improve resistance to multiple drought stresses.

Functions of ABA in dehydration stress memory

ABA is strongly induced by drought stress and plays key roles in the drought stress response [17, 70]. In Arabidopsis, endogenous ABA accumulated in D1, became less abundant in R1, and accumulated again in D2. However, the extent of its accumulation in D1 was identical to that in D2 [28]. After 7 days of watering, ABA levels could be restored to their prior peak despite the loss of dehydration stress memory [28]. In the areb1/areb2/abf3 Arabidopsis triple mutant, the expressed levels of trained marker genes were strongly reduced but remained much higher than the transcript level in D1 [38]. These results indicate that ABA is not the key factor in the induction of transcriptional memory in Arabidopsis. In this work, ABA accumulated in D1, R1 and D2, and especially strongly in R1 and D2 (Fig. 5d). In contrast to the observations for Arabidopsis, there was significant ABA accumulation during R1 and D2 in plants subjected to repeated dehydration stress, indicating that ABA may play different roles in switchgrass and Arabidopsis, and revealing a species-specific differences between switchgrass (a C4 monocot plant) and Arabidopsis (a C3 eudicot plant).

JA accumulation and its crosstalk with ABA

JA is an important hormone that plays key roles in plant responses to biotic and abiotic stresses [61, 71]. It is synthesized rapidly when a plant suffers drought, cold or other stresses [59, 72]. In switchgrass, JA-related genes (especially downstream genes expressed in the peroxisome) were significantly up-regulated in D2 (Fig. 5b; Additional file 6: Table S5). This suggests that JA may play an important role in the response to repeated dehydration stresses. However, HPLC–MS data suggested that JA was accumulated in D1 but recovered to its initial levels in R1 and D2 (Fig. 5d). It is not clear why JA accumulated in D1 but not in D2, because the downstream genes of JA biosynthesis were up-regulated in D2.

JA is synthesized from α-linolenic acid in chloroplasts, and the step that commits the substrate to JA formation is the formation of 13-hydroperoxy-9,11,15-octadecatrienoic acid (13-HPOT) catalyzed by allene oxide synthase (AOS) [61]. The expression of AOS is controlled by a bHLH transcriptional factor, MYC2 [40]. In switchgrass, two orthologs of MYC2 were slightly up-regulated in D1 and down-regulated in D2, and another four orthologs of MYC2 were down-regulated in both D1 and D2 (Fig. 6; Additional file 6: Table S5). Previous studies revealed that transient JA accumulation is needed for the increase in ABA levels in citrus roots under drought stress [60]. Therefore, we inferred that when the switchgrass plants experienced the first dehydration stress, MYC2 was rapidly up-regulated and induced the expression of AOS, leading to the quick accumulation of JA. Conversely, during secondary dehydration stress, all the orthologs of MYC2 were down-regulated, suppressing the transcription of AOS and thus reducing JA accumulation even though the downstream genes of JA biosynthesis were up-regulated.

There is sophisticated crosstalk between JA and ABA, with both synergistic and antagonistic effects in plants. For example, dehydration stress and ABA can induce JA synthesis [73,74,75], and ABA plays a key role in inducing MYC2 expression in response to JA in Arabidopsis [40]. Conversely, in citrus roots, increased ABA levels were observed following transient JA accumulation [60]. In this work, JA accumulated during D1 and ABA accumulated in D1, R1 and D2, especially in the latter two periods (Fig. 5d). These results suggest that the transient accumulation of JA may induce strong ABA accumulation, but ABA accumulation is not sufficient to induce JA accumulation. In addition, MYC2 is critical for JA biosynthesis but it was first identified as a component of an ABA signal transduction system [17, 40], suggesting that it may play an important roles in the crosstalk between JA and ABA during repeated dehydration stress (Fig. 6). Further studies on this signaling network are warranted.

Osmotic adjustment, energy metabolism and lignin biosynthesis

Osmotic adjustment is a common plant response to water-deficit stress, but plant species differ in their patterns of osmolyte biosynthesis. For instance, proline, raffinose and galactinol accumulated in Arabidopsis and provided protection against drought, salinity, and cold stress [76, 77]. Conversely, raffinose and trehalose (but not proline) were synthesized in the leaves of the physic nut under severe drought stress [45]. In switchgrass subjected to dehydration stress, genes related to the biosynthesis of proline, raffinose family oligosaccharides, and trehalose were strongly up-regulated in leaves (Fig. 7; Additional file 8: Table S7), indicating that these substances play important roles in osmotic adjustment, especially in secondary dehydration stress. The results of proline content determination assay demonstrated that proline was significantly accumulated in D2, and these data are consistent with the related gene expression (Fig. 7a, b).

Plants may adjust their energy metabolism to increase the supply of ATP and related compounds [45, 78]. The physiological assays conducted in this work and the transcriptional data for PSI and PSII showed that photosynthesis and the Calvin cycle proceeded normally in D1, but were significantly suppressed in D2 (Fig. 1; Additional file 9: Table S8). In addition, orthologs of several key enzymes in glycolysis and the TCA cycle were differentially expressed during secondary dehydration stress, causing changes in the ATP supply and intermediate metabolism. The affected enzymes include hexokinase, 6-phosphofructokinase, pyruvate kinase, citrate (Si)-synthase, aconitate hydratase and malate dehydrogenase (Additional file 10: Table S9). The different responses during D1 and D2 demonstrate the transcriptional memory of switchgrass.

Switchgrass is a lignocellulose biofuel plant, and genes important in its biosynthesis of lignin have been identified in vivo and in vitro [63]. During drought stress, microRNAs regulate its lignin and biomass synthesis [14, 79], and the total plant lignin concentrations increased during drought [80]. In this work, genes related to lignin biosynthesis were slightly down-regulated in D1 but significantly up-regulated in D2, and thus exhibited [=/+] or [−/+] memory behavior (Fig. 8; Additional file 11: Table S10). The data on the transcription of lignin biosynthesis genes acquired in this work will provide an excellent foundation for the genetic modification of switchgrass.

Dehydration stress memory in Arabidopsis, maize and switchgrass

Transcriptional memory persisted for ~ 5 days in Arabidopsis [28] and maize [31], and was associated with trimethylated histone H3 Lys4 (H3K4me3) and a stalled RNA polymerase II [28]. In this work, the expressions of representative dehydration stress memory genes recovered to their initial levels after 7 days of rewatering (Fig. 9). It this appears that transcriptional memory in switchgrass persists for ~ 5 days, which is enough to accommodate the diurnal cycle of dehydration stress, and probably plays important roles in the response to long-term successive dehydration stresses.

Dehydration stress memory has been observed in several species, but DMGs were identified only in Arabidopsis and maize [30, 31]. In this study, we identified DMGs in switchgrass, and identified interesting similarities and differences between these three species: (i) the number of response genes in switchgrass was similar to that in maize, but much lower than that for Arabidopsis, whereas the number of late-response genes in switchgrass was significantly greater than in the other two species. (ii) The ratio of DMGs to all response genes in switchgrass (47.3%) was higher than the ratios in maize (39.6%) and Arabidopsis (29.8%). (iii) Orthologs of switchgrass DMGs in maize mainly belonged to other gene types (non-response, non-memory, or late-response genes); even those that contributed to dehydration stress memory in two species generally functioned as different types of memory gene in each. (iv) 14.2% of all orthologs of switchgrass DMGs in maize exhibited dehydration stress memory, and each type had orthologs belonging to the same memory type in maize. Nearly all the orthologs (93.9%) of [+/−]-type genes belonged to the same memory type in maize. These memory-conserved genes demonstrated the conservation of dehydration stress memory during evolution, and are likely to be excellent targets for genetic modification in switchgrass or related crops in the family Gramineae (Additional file 12: Table S11). With the methods of overexpression, drought-induced expression or gene editing, the [+/+] and [−/−] memory genes are considered as excellent targets to enhance drought tolerance of switchgrass and other crops.

Conclusions

In this study, we identified 741 dehydration stress memory genes in switchgrass plants subjected to multiple dehydration stresses. Analyses of specific pathways combined with physiological data allowed us to explore the molecular mechanisms of repeated dehydration stress responses in switchgrass at the transcriptional and physiological levels. Switchgrass plants acclimatized to secondary dehydration stresses by regulating the crosstalk of ABA and JA biosynthesis and signal transduction, adjusting the biosynthesis of osmotic substances (l-proline, stachyose and trehalose), changing their energy metabolism (by modifying the expression of genes relating to photosynthesis systems, glycolysis, and the TCA cycle), and through lignin biosynthesis. The transcriptional data showed that dehydration memory persisted for ~ 5 days in switchgrass. Analyses of homologous genes in switchgrass, maize and Arabidopsis revealed that some memory genes are conserved across species, but there are also important species-specific aspects of stress memory. These results revealed a previously unknown transcriptional memory behavior in switchgrass, and the genes and pathways identified in this study will be useful for the future studies of switchgrass and other crops.

Abbreviations

ABA:

abscisic acid

BLAST:

basic local alignment search tool

cDNA:

complementary DNA

D1–3:

the first to third dehydration stress

DEG:

differentially expressed genes

DMG:

dehydration stress memory gene

FPKM:

fragments per kilobase of exon per million fragments mapped

GO:

Gene Ontology

HPLC–MS:

high performance liquid chromatography–mass spectrometry

JA:

jasmonic acid

KEGG:

Kyoto encyclopedia of genes and genomes database

nr:

NCBI non-redundant protein sequences database

PSI and PSII:

photosynthesis system I and II

qN:

non-photochemical quenching

qRT-PCR:

quantitative real-time PCR

R1–2:

the first and second dehydration stress

RGAP:

rice genome annotation project

RNA-Seq:

RNA sequencing

RWC:

Relative water content

TAIR10:

Arabidopsis information resource proteins database

TCA cycle:

tricarboxylic acid cycle

References

  1. Mclaughlin SB, Kszos LA. Development of switchgrass (Panicum virgatum) as a bioenergy feedstock in the United States. Biomass Bioenerg. 2005;28:515–35.

    Article  Google Scholar 

  2. Wu Z, Cao Y, Yang R, Qi T, Hang Y, Lin H, Zhou G, Wang ZY, Fu C. Switchgrass SBP-box transcription factors PvSPL1 and 2 function redundantly to initiate side tillers and affect biomass yield of energy crop. Biotechnol Biofuels. 2016;9:101.

    Article  Google Scholar 

  3. David K, Ragauskas AJ. Switchgrass as an energy crop for biofuel production: a review of its ligno-cellulosic chemical properties. Energy Environ Sci. 2010;3:1182–90.

    Article  CAS  Google Scholar 

  4. Hultquist SJ, Vogel KP, Lee DJ, Arumuganathan K, Kaeppler S. Chloroplast DNA and nuclear DNA content variations among cultivars of switchgrass, Panicum virgatum L. Crop Sci. 1996;36:1049–52.

    Article  Google Scholar 

  5. Wright L, Turhollow A. Switchgrass selection as a “model” bioenergy crop: a history of the process. Biomass Bioenerg. 2010;34:851–68.

    Article  Google Scholar 

  6. Sanderson MA, Reed RL, McLaughlin SB, Wullschleger SD, Conger BV, Parrish DJ, Wolf DD, Taliaferro C, Hopkins AA, Ocumpaugh WR, et al. Switchgrass as a sustainable bioenergy crop. Bioresour Technol. 1996;56:83–93.

    Article  CAS  Google Scholar 

  7. Lewandowski I, Scurlock JMO, Lindvall E, Christou M. The development and current status of perennial rhizomatous grasses as energy crops in the US and Europe. Biomass Bioenergy Biomass Bioenergy. 2003;25:335–61.

    Article  Google Scholar 

  8. Casler MD, Tobias CM, Kaeppler SM, Buell CR, Wang ZY, Cao PJ, Schmutz J, Ronald P. The switchgrass genome: tools and strategies. Plant Genome. 2011;4:273–82.

    Article  CAS  Google Scholar 

  9. Zhang JY, Lee YC, Torres-Jerez I, Wang M, Yin Y, Chou WC, He J, Shen H, Srivastava AC, Pennacchio C, et al. Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L.). Plant J. 2013;74:160–73.

    Article  CAS  Google Scholar 

  10. Somleva MN. Switchgrass (Panicum virgatum L.). Methods Mol Biol. 2006;344:65–73.

    CAS  Google Scholar 

  11. Aimar D, Calafat M, Andrade AM, Carassay L, Bouteau F, Abdala G, Molas ML. Drought effects on the early development stages of Panicum virgatum L.: cultivar differences. Biomass Bioenergy. 2014;66:49–59.

    Article  CAS  Google Scholar 

  12. Meyer E, Aspinwall MJ, Lowry DB, Palacio-Mejia JD, Logan TL, Fay PA, Juenger TE. Integrating transcriptional, metabolomic, and physiological responses to drought stress and recovery in switchgrass (Panicum virgatum L.). BMC Genomics. 2014;15:527.

    Article  Google Scholar 

  13. Liu Y, Zhang X, Tran H, Shan L, Kim J, Childs K, Ervin EH, Frazier T, Zhao B. Assessment of drought tolerance of 49 switchgrass (Panicum virgatum) genotypes using physiological and morphological parameters. Biotechnol Biofuels. 2015;8:152.

    Article  Google Scholar 

  14. Xie FL, Stewart CN, Taki FA, He QL, Liu HW, Zhang BH. High-throughput deep sequencing shows that microRNAs play important roles in switchgrass responses to drought and salinity stress. Plant Biotechnol J. 2014;12:354–66.

    Article  CAS  Google Scholar 

  15. Sun GL, Stewart CN, Xiao P, Zhang BH. MicroRNA expression analysis in the cellulosic biofuel crop switchgrass (Panicum virgatum) under abiotic stress. PLoS ONE. 2012;7:e32017.

    Article  CAS  Google Scholar 

  16. Riera M, Valon C, Fenzi F, Giraudat J, Leung J. The genetics of adaptive responses to drought stress: abscisic acid-dependent and abscisic acid-independent signalling components. Physiol Plant. 2005;123:111–9.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  18. Behnam B, Iuchi S, Fujita M, Fujita Y, Takasaki H, Osakabe Y, Yamaguchi-Shinozaki K, Kobayashi M, Shinozaki K. Characterization of the promoter region of an Arabidopsis gene for 9-cis-epoxycarotenoid dioxygenase involved in dehydration-inducible transcription. DNA Res. 2013;20:315–24.

    Article  CAS  Google Scholar 

  19. Nambara E, Marion-Poll A. Abscisic acid biosynthesis and catabolism. Annu Rev Plant Biol. 2005;56:165–85.

    Article  CAS  Google Scholar 

  20. Seki M, Umezawa T, Urano K, Shinozaki K. Regulatory metabolic networks in drought stress responses. Curr Opin Plant Biol. 2007;10:296–302.

    Article  CAS  Google Scholar 

  21. Lu PL, Chen NZ, An R, Su Z, Qi BS, Ren F, Chen J, Wang XC. A novel drought-inducible gene, ATAF1, encodes a NAC family protein that negatively regulates the expression of stress-responsive genes in Arabidopsis. Plant Mol Biol. 2007;63:289–305.

    Article  CAS  Google Scholar 

  22. Pizzio GA, Rodriguez L, Antoni R, Gonzalez-Guzman M, Yunta C, Merilo E, Kollist H, Albert A, Rodriguez PL. The PYL4 A194T mutant uncovers a key role of PYR1-like4/protein phosphatase 2CA interaction for abscisic acid signaling and plant drought resistance. Plant Physiol. 2013;163:441–55.

    Article  CAS  Google Scholar 

  23. Choi H, Hong J, Ha J, Kang J, Kim SY. ABFs, a family of ABA-responsive element binding factors. J Biol Chem. 2000;275:1723–30.

    Article  CAS  Google Scholar 

  24. Simpson SD, Nakashima K, Narusaka Y, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Two different novel cis-acting elements of erd1, a clpA homologous Arabidopsis gene function in induction by dehydration stress and dark-induced senescence. Plant J. 2003;33:259–70.

    Article  CAS  Google Scholar 

  25. Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998;10:1391–406.

    Article  CAS  Google Scholar 

  26. Tran LS, Nakashima K, Sakuma Y, Osakabe Y, Qin F, Simpson SD, Maruyama K, Fujita Y, Shinozaki K, Yamaguchi-Shinozaki K. Co-expression of the stress-inducible zinc finger homeodomain ZFHD1 and NAC transcription factors enhances expression of the ERD1 gene in Arabidopsis. Plant J. 2007;49:46–63.

    Article  CAS  Google Scholar 

  27. Bruce TJA, Matthes MC, Napier JA, Pickett JA. Stressful “memories” of plants: evidence and possible mechanisms. Plant Sci. 2007;173:603–8.

    Article  CAS  Google Scholar 

  28. Ding Y, Fromm M, Avramova Z. Multiple exposures to drought ‘train’ transcriptional responses in Arabidopsis. Nat Commun. 2012;3:740.

    Article  Google Scholar 

  29. Jaskiewicz M, Conrath U, Peterhänsel C. Chromatin modification acts as a memory for systemic acquired resistance in the plant stress response. EMBO Rep. 2011;12:50–5.

    Article  CAS  Google Scholar 

  30. Ding Y, Liu N, Virlouvet L, Riethoven JJ, Fromm M, Avramova Z. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 2013;13:229.

    Article  Google Scholar 

  31. Ding Y, Virlouvet L, Liu N, Riethoven JJ, Fromm M, Avramova Z. Dehydration stress memory genes of Zea mays; comparison with Arabidopsis thaliana. BMC Plant Biol. 2014;14:141.

    Article  Google Scholar 

  32. Bilichak A, Ilnytskyy Y, Woycicki R, Kepeshchuk N, Fogen D, Kovalchuk I. The elucidation of stress memory inheritance in Brassica rapa plants. Front Plant Sci. 2015;6:5.

    Google Scholar 

  33. Ramirez DA, Rolando JL, Yactayo W, Monneveux P, Mares V, Quiroz R. Improving potato drought tolerance through the induction of long-term water stress memory. Plant Sci. 2015;238:26–32.

    Article  CAS  Google Scholar 

  34. Fleta-Soriano E, Pinto-Marijuan M, Munne-Bosch S. Evidence of drought stress memory in the facultative CAM, Aptenia cordifolia: possible role of phytohormones. PLoS ONE. 2015;10:e0135391.

    Article  Google Scholar 

  35. Pinto-Marijuan M, Cotado A, Fleta-Soriano E, Munne-Bosch S. Drought stress memory in the photosynthetic mechanisms of an invasive CAM species, Aptenia cordifolia. Photosynth Res. 2017;131:241–53.

    Article  CAS  Google Scholar 

  36. Lamke J, Brzezinka K, Altmann S, Baurle I. A hit-and-run heat shock factor governs sustained histone methylation and transcriptional stress memory. EMBO J. 2016;35:162–75.

    Article  Google Scholar 

  37. Liu N, Fromm M, Avramova Z. H3K27me3 and H3K4me3 chromatin environment at super-induced dehydration stress memory genes of Arabidopsis thaliana. Mol Plant. 2014;7:502–13.

    Article  CAS  Google Scholar 

  38. Virlouvet L, Ding Y, Fujii H, Avramova Z, Fromm M. ABA signaling is necessary but not sufficient for RD29B transcriptional memory during successive dehydration stresses in Arabidopsis thaliana. Plant J. 2014;79:150–61.

    Article  CAS  Google Scholar 

  39. Feng XJ, Li JR, Qi SL, Lin QF, Jin JB, Hua XJ. Light affects salt stress-induced transcriptional memory of P5CS1 in Arabidopsis. Proc Natl Acad Sci USA. 2016;113:E8335–43.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  41. Walter J, Nagy L, Hein R, Rascher U, Beierkuhnlein C, Willner E, Jentsch A. Do plants remember drought? Hints towards a drought-memory in grasses. Environ Exp Bot. 2011;71:34–40.

    Article  Google Scholar 

  42. Xi YJ, Fu CX, Ge YX, Nandakumar R, Hisano H, Bouton J, Wang ZY. Agrobacterium-mediated transformation of switchgrass and inheritance of the transgenes. Bioenergy Res. 2009;2:275–83.

    Article  Google Scholar 

  43. Xu K, Sun F, Wang YW, Chai G, Shi L, Liu S, Xi Y. De novo assembly and transcriptome analysis of two contrary tillering mutants to learn the mechanisms of tillers outgrowth in switchgrass (Panicum virgatum L.). Front Plant Sci. 2015;6:749.

    Google Scholar 

  44. Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31:46.

    Article  CAS  Google Scholar 

  45. Zhang C, Zhang L, Zhang S, Zhu S, Wu P, Chen Y, Li M, Jiang H, Wu G. Global analysis of gene expression profiles in physic nut (Jatropha curcas L.) seedlings exposed to drought stress. BMC Plant Biol. 2015;15:17.

    Article  CAS  Google Scholar 

  46. McCarthy FM, Wang N, Magee GB, Nanduri B, Lawrence ML, Camon EB, Barrell DG, Hill DP, Dolan ME, Williams WP. AgBase: a functional genomics resource for agriculture. BMC Genomics. 2006;7:229.

    Article  Google Scholar 

  47. Huang LK, Yan HD, Jiang XM, Zhang XQ, Zhang YW, Huang X, Zhang Y, Miao JM, Xu B, Frazier T, Zhao BY. Evaluation of candidate reference genes for normalization of quantitative RT-PCR in switchgrass under various abiotic stress conditions. Bioenergy Res. 2014;7:1201–11.

    Article  CAS  Google Scholar 

  48. Gimeno J, Eattock N, Deynze AV, Blumwald E. Selection and validation of reference genes for gene expression analysis in switchgrass (Panicum virgatum) using quantitative real-time RT-PCR. PLoS ONE. 2014;9:e91474.

    Article  Google Scholar 

  49. Gao S, Zheng Z, Gu W, Xie X, Huan L, Pan G, Wang G. Photosystem I shows a higher tolerance to sorbitol-induced osmotic stress than photosystem II in the intertidal macro-algae Ulva prolifera (Chlorophyta). Phys Plant. 2014;152:380–8.

    Article  CAS  Google Scholar 

  50. Duncan DB. Multiple range and multiple F tests. Biometrics. 1955;11:1–42.

    Article  Google Scholar 

  51. Baker NR. Chlorophyll fluorescence: a probe of photosynthesis in vivo. Annu Rev Plant Biol. 2008;59:89–113.

    Article  CAS  Google Scholar 

  52. Krasensky J, Jonak C. Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J Exp Bot. 2012;63:1593.

    Article  CAS  Google Scholar 

  53. Fuchs S, Grill E, Meskiene I, Schweighofer A. Type 2C protein phosphatases in plants. FEBS J. 2013;280:681.

    Article  CAS  Google Scholar 

  54. Fujita Y, Fujita M, Shinozaki K, Yamaguchishinozaki K. ABA-mediated transcriptional regulation in response to osmotic stress in plants. J Plant Res. 2011;124:509.

    Article  CAS  Google Scholar 

  55. Kulik A, Wawer I, Krzywińska E, Bucholc M, Dobrowolska G. SnRK2 protein kinases—key regulators of plant response to abiotic stresses. OMICS J Integr Biol. 2011;15:859.

    Article  CAS  Google Scholar 

  56. Xie T, Ren R, Zhang YY, Pang Y, Yan C, Gong X, He Y, Li W, Miao D, Hao Q. Molecular mechanism for inhibition of a critical component in the Arabidopsis thaliana abscisic acid signal transduction pathways, SnRK2.6, by protein phosphatase ABI1. J Biol Chem. 2012;287:794–802.

    Article  CAS  Google Scholar 

  57. Mueller MJ. Enzymes involved in jasmonic acid biosynthesis. Physiol Plant. 1997;100:653–63.

    Article  CAS  Google Scholar 

  58. Wasternack C. Jasmonates in stress, growth, and development. In: Hirt H, editor. Plant stress biology. Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA; 2010. p. 91–118.

    Google Scholar 

  59. Fang L, Su L, Sun X, Li X, Sun M, Karungo SK, Fang S, Chu J, Li S, Xin H. Expression of Vitis amurensis NAC26 in Arabidopsis enhances drought tolerance by modulating jasmonic acid synthesis. J Exp Bot. 2016;67:2829–45.

    Article  CAS  Google Scholar 

  60. de Ollas C, Hernando B, Arbona V, Gomez-Cadenas A. Jasmonic acid transient accumulation is needed for abscisic acid increase in citrus roots under drought stress conditions. Physiol Plant. 2013;147:296–306.

    Article  Google Scholar 

  61. Chehab EW, Braam J. Jasmonates in plant defense responses. In: Witzany G, Baluška F, editor. Biocommunication of plants. Berlin, Heidelberg: Springer; 2012. p. 67–88.

    Chapter  Google Scholar 

  62. Ashraf M, Akram NA, Al-Qurainy F, Foolad MR. Drought tolerance: roles of organic osmolytes, growth regulators, and mineral nutrients. Adv Agron. 2011;111:249–96.

    Article  CAS  Google Scholar 

  63. Shen H, Mazarei M, Hisano H, Escamilla-Trevino L, Fu CX, Pu YQ, Rudis MR, Tang YH, Xiao XR, Jackson L, et al. A genomics approach to deciphering lignin biosynthesis in switchgrass. Plant Cell. 2013;25:4342–61.

    Article  CAS  Google Scholar 

  64. Xu B, Escamilla-Trevino LL, Sathitsuksanoh N, Shen ZX, Shen H, Zhang YHP, Dixon RA, Zhao BY. Silencing of 4-coumarate:coenzyme A ligase in switchgrass leads to reduced lignin content and improved fermentable sugar yields for biofuel production. New Phytol. 2011;192:611–25.

    Article  CAS  Google Scholar 

  65. Escamilla-Treviño LL, Shen H, Uppalapati SR, Ray T, Tang Y, Hernandez T, Yin Y, Xu Y, Dixon RA. Switchgrass (Panicum virgatum) possesses a divergent family of cinnamoyl CoA reductases with distinct biochemical properties. New Phytol. 2010;185:143–55.

    Article  Google Scholar 

  66. Fu C, Mielenz JR, Xiao X, Ge Y, Hamilton CY, Rodriguez M Jr, Chen F, Foston M, Ragauskas A, Bouton J, et al. Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc Natl Acad Sci USA. 2011;108:3803–8.

    Article  CAS  Google Scholar 

  67. Fu C, Xiao X, Xi Y, Ge Y, Chen F, Bouton J, Dixon RA, Wang Z-Y. Downregulation of cinnamyl alcohol dehydrogenase (CAD) leads to improved saccharification efficiency in switchgrass. Bioenergy Res. 2011;4:153–64.

    Article  Google Scholar 

  68. Ruggiero C, Angelino G, Maggio A. Developmental regulation of water uptake in wheat. J Plant Physiol. 2007;164:1170–8.

    Article  CAS  Google Scholar 

  69. Henry C, Bledsoe SW, Siekman A, Kollman A, Waters BM, Feil R, Stitt M, Lagrimini LM. The trehalose pathway in maize: conservation and gene regulation in response to the diurnal cycle and extended darkness. J Exp Bot. 2014;65:5959–73.

    Article  CAS  Google Scholar 

  70. Micco VD, Aronne G. Plant responses to drought stress: from morphological to molecular features. In: Aroca R, editor. Chapter 2 Morpho-anatomical traits for plant adaptation to drought, vol. II. Berlin: Springer; 2012. p. 466.

    Google Scholar 

  71. Sasaki-Sekimoto Y, Taki N, Obayashi T, Aono M, Matsumoto F, Sakurai N, Suzuki H, Hirai MY, Noji M, Saito K. Coordinated activation of metabolic pathways for antioxidants and defence compounds by jasmonates and their roles in stress tolerance in Arabidopsis. Plant J. 2005;44:653–68.

    Article  CAS  Google Scholar 

  72. Li Q, Lei S, Du K, Li L, Pang X, Wang Z, Wei M, Fu S, Hu L, Xu L. RNA-seq based transcriptomic analysis uncovers alpha-linolenic acid and jasmonic acid biosynthesis pathways respond to cold acclimation in Camellia japonica. Sci Rep. 2016;6:36463.

    Article  CAS  Google Scholar 

  73. Creelman RA, Mullet JE. Jasmonic acid distribution and action in plants: regulation during development and response to biotic and abiotic stress. Proc Natl Acad Sci USA. 1995;92:4114–9.

    Article  CAS  Google Scholar 

  74. Adie B, Perez-Perez J, Perez-Perez M, Godoy M, Sanchez-Serrano J, Schmelz EA, Solano R. ABA is an essential signal for plant resistance to pathogens affecting JA biosynthesis and the activation of defenses in Arabidopsis. Plant Cell. 2007;19:1665–81.

    Article  CAS  Google Scholar 

  75. Munoz-Espinoza VA, Lopez-Climent MF, Casaretto JA, Gomez-Cadenas A. Water stress responses of tomato mutants impaired in hormone biosynthesis reveal abscisic acid, jasmonic acid and salicylic acid interactions. Front Plant Sci. 2015;6:997.

    Article  Google Scholar 

  76. Oba Y, Sanada Y, Wada K. Proline accumulation protects against drought stress tolerance in Arabidopsis thaliana. Plant Physiol. 1997;114:471–471.

    Google Scholar 

  77. Taji T, Ohsumi C, Iuchi S, Seki M, Kasuga M, Kobayashi M, Yamaguchi-Shinozaki K, Shinozaki K. Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 2002;29:417–26.

    Article  CAS  Google Scholar 

  78. Cramer GR, Ergul A, Grimplet J, Tillett RL, Tattersall EA, Bohlman MC, Vincent D, Sonderegger J, Evans J, Osborne C, et al. Water and salinity stress in grapevines: early and late changes in transcript and metabolite profiles. Funct Integr Genomics. 2007;7:111–34.

    Article  CAS  Google Scholar 

  79. Trumbo JL, Zhang BH, Stewart CN. Manipulating microRNAs for improved biomass and biofuels from plant feedstocks. Plant Biotechnol J. 2015;13:337–54.

    Article  CAS  Google Scholar 

  80. Jiang YW, Yao Y, Wang Y. Physiological response, cell wall components, and gene expression of switchgrass under short-term drought stress and recovery. Crop Sci. 2012;52:2718–27.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

YX, SL, FS and CZ conceived and designed the experiments. XP, XG and CZ carried out the experiments. GT and CZ analyzed the data. CZ, SL and YX wrote and revised the paper. All authors read and approved the final manuscript.

Acknowledgements

This research was supported by grants from the National Natural Science Foundation of China (Nos. 31601370, 31301384, 31371690), China Postdoctoral Science Foundation (No. 171305), PhD Start-up Foundation (Z109021616) and Basic Scientific Research Foundation (Z109021705) of Northwest A&F University.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The referenced genomic sequences are available at plant comparative genomics portal of the department of energy’s joint genome institute (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Pvirgatum), and the transcriptomic data were uploaded to sequence read archive (SRA) at NCBI (Accession Number PRJNA394165).

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yajun Xi.

Additional files

Additional file 1: Figure S1.

Overview of the RNA-Seq results and repeatability of different biological replicates. The first eight figures show the proportion of clean reads in the sequenced samples; the next three figures show the repeated expression of genes in two biological replicates; the last four figures showed the correlations between differentially expressed genes in two biological replicates. C, control; D1-2, the first and second dehydration stresses; R1, the first recovery period.

Additional file 2: Table S1.

List of dehydration memory genes (DMGs) identified in this work.

Additional file 3: Table: S2.

Primers used in qRT-PCR.

Additional file 4: Table S3.

Top 20 pathways involved in dehydration memory.

Additional file 5: Table S4.

Differentially expressed genes involved in ABA biosynthesis, degradation and signal transduction.

Additional file 6: Table S5.

Differentially expressed genes involved in JA biosynthesis and signal transduction.

Additional file 7: Table S6.

Differentially expressed genes involved in other plant hormones biosynthesis and signal transduction.

Additional file 8: Table S7.

Differentially expressed genes involved in the biosynthesis of L-proline, stachyose and trehalose.

Additional file 9: Table S8.

Differentially expressed genes involved in photosynthesis.

Additional file 10: Table S9.

Differentially expressed genes involved in glycolysis and TCA cycle.

Additional file 11: Table S10.

Differentially expressed genes involved in lignin biosynthesis.

Additional file 12: Table S11.

Common DMGs in switchgrass and maize.

Additional file 13: Table S12.

List of DMGs that are different between switchgrass and maize.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, C., Peng, X., Guo, X. et al. Transcriptional and physiological data reveal the dehydration memory behavior in switchgrass (Panicum virgatum L.). Biotechnol Biofuels 11, 91 (2018). https://doi.org/10.1186/s13068-018-1088-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-018-1088-x

Keywords