- Research article
- Open Access
Transcriptome reprogramming during severe dehydration contributes to physiological and metabolic changes in the resurrection plant Haberlea rhodopensis
BMC Plant Biology volume 18, Article number: 351 (2018)
Water shortage is a major factor that harms agriculture and ecosystems worldwide. Plants display various levels of tolerance to water deficit, but only resurrection plants can survive full desiccation of their vegetative tissues. Haberlea rhodopensis, an endemic plant of the Balkans, is one of the few resurrection plants found in Europe. We performed transcriptomic analyses of this species under slight, severe and full dehydration and recovery to investigate the dynamics of gene expression and associate them with existing physiological and metabolomics data.
De novo assembly yielded a total of 142,479 unigenes with an average sequence length of 1034 nt. Among them, 18,110 unigenes were differentially expressed. Hierarchical clustering of all differentially expressed genes resulted in seven clusters of dynamic expression patterns. The most significant expression changes, involving more than 15,000 genes, started at severe dehydration (~ 20% relative water content) and were partially maintained at full desiccation (< 10% relative water content). More than a hundred pathways were enriched and functionally organized in a GO/pathway network at the severe dehydration stage. Transcriptomic changes in key pathways were analyzed and discussed in relation to metabolic processes, signal transduction, quality control of protein and DNA repair in this plant during dehydration and rehydration.
Reprograming of the transcriptome occurs during severe dehydration, resulting in a profound alteration of metabolism toward alternative energy supply, hormone signal transduction, and prevention of DNA/protein damage under very low cellular water content, underlying the observed physiological and metabolic responses and the resurrection behavior of H. rhodopensis.
Environmental stress, especially drought, place limitations on plant distribution and crop production worldwide. Intensive breeding efforts, based on crosses with wild drought-tolerant relatives as potential sources of useful traits, are underway. At the same time, interest in a small group called resurrection plants is increasing, as these plants can withstand long periods of almost full desiccation and recover their vegetative system activities quickly and fully upon rehydration. This makes them useful models to study desiccation tolerance and potential sources of genes, related to these important phenomena [1,2,3,4]. Further efforts are needed to elucidate the mechanisms of desiccation tolerance in resurrection plants as a sustained dynamic process.
The Balkan endemic Haberlea rhodopensis, one of the few species of resurrection plants, native to Europe, has been the subject of intensive genetic studies in recent years [5,6,7]. Transcriptomic analysis of H. rhodopensis gene expression under conditions of normal watering, dehydration (50% relative water content; RWC), full desiccation (5% RWC) and rehydration has been performed previously . We have shown that the dynamics of desiccation tolerance can be characterized by numerous changes in various processes at additional important time points during drying: the stages of slight (~ 65–75%) and severe (20–25% RWC) dehydration . For example, at the onset of drying (at 75% RWC), an approximately 25-fold increase in jasmonic acid (JA)  occurs, coinciding with the start of the decrease in photosynthetic performance . On the other hand, oxygen evolution in leaves of H. rhodopensis did not decrease until 20% RWC . Our studies revealed that the levels of numerous compounds changed dramatically during desiccation, particularly before or around the severe dehydration stage (25% RWC), including total phenols, sugars, total glutathione, malondialdehyde, sucrose/fructose ratio, abscisic acid (ABA), salicylic acid (SA) and starch levels [9, 11, 12]. Dramatic changes in photosynthetic performance and energy supply, such as ATP and glycolytic intermediates, were also observed at this stage of dehydration [9,10,11,12]. This time point appears to be of crucial importance, as non-resurrection plants related to H. rhodopensis tolerated the loss of almost two-thirds of their water content, whereas further drying to 20% RWC was irreversible .
In an attempt to dissect the dynamic transcriptional regulation events preceding the dramatic physiological and metabolic changes that occur during desiccation in H. rhodopensis, we performed de novo transcriptome sequencing at new time points during drying from the fully hydrated to desiccated state: slight and severe dehydration stages.
Transcriptome sequencing and assembly
Transcriptomic analysis of gene expression in plants that were regularly watered (fresh control, F), slightly dehydrated to 75% RWC (D75), severely dehydrated to 20% RWC (D20), fully dehydrated to < 10% RWC (DT) and fully rehydrated after DT (R) was performed using Illumina next-generation sequencing technology. De novo assembly yielded a total of 142,479 unigenes, with an average sequence length of 1034 nt and an N50 value of 1664 nt (Additional file 1: Table S1). Approximately 75% of all unigenes were novel to H. rhodopensis. The sequence length distribution of the unigenes is indicated in Additional file 2: Fig. S1a. The numbers of unigenes annotated using several databases, NCBI non-redundant protein sequences, NCBI nucleotide collection, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), the Clusters of Orthologous Group (COG) and Gene Ontology (GO), ranged from 91,753 to 51,046. Among the matching sequences, 61% showed the closest match to sequences of Sesamum indicum, a species belonging to the family Pedaliaceae, order Scrophulariales (Additional file 2: Fig. S1b) . There were 46,922 transcripts with no sequence similarity to known genes (orphan genes).
Analysis of differentially expressed genes (DEGs)
Differential gene expression analysis was performed to monitor transcriptomic variations in the plant under five dehydration treatments. From the 10 libraries, 41,279,554 to 55,613,670 clean sequence reads were obtained after trimming (Additional file 3: Table S2), which were sufficient for quantitative analysis of gene expression. Based on the fragments per kilobase of transcript per million mapped reads (FPKM) of the unigenes in each sample, 18,110 unigenes were identified as DEGs, of which approximately 40% (6877) were novel.
Higher numbers of DEGs were observed in the comparisons of D20 with F (7,675 up and 8894 down), D75 (3.712 up and 4.120 down) or R (8,496 up and 6958 down) and in the comparisons of DT with F (4,527 up and 5998 down) or R (5,517 up and 3477 down). The lowest number of DEGs was found in the comparison of F with R (514 up and 290 down) (Fig. 1a). Among these DEGs, the highest levels of transcripts encoding LEA (late embryogenesis abundant) proteins, early light-inducible proteins, galactinol synthases, beta-amylases, sucrose synthases, senescence-associated proteins and raffinose synthases were seen at D20, whereas catalases were induced after slight dehydration (Additional file 4: Table S3).
A Venn diagram of the DEGs illustrates their occurrence in each dehydration treatment (Fig. 1b). The largest number of unique genes was found in the D20 (7065), followed by DT (1210), R (129) and D75 (153) treatments, whereas only 119 genes were present in all dehydration treatments. More than 80% of novel DEGs were either unique to D20 or shared by D20 and DT (Additional file 2: Fig. S1c). Based on their FPKM values, a random selection of 20 genes encoding enzymes, transcription factors and regulators involved in metabolism and stress responses were subjected to quantitative real-time PCR (qPCR) analysis to validate the DEGs identified by RNA-seq. The expression levels of these genes were significantly correlated with those determined by qPCR (r = 0.99, p < 0.01) (Additional file 5: Fig. S2a, Additional file 6: Table S4).
Principal component analysis of transcript profiling of dehydrated and rehydrated plants
Principal component analysis was performed to capture the overall variance among samples. Samples in each treatment generated a cluster pattern according to water status (Fig. 2). The F, D75 and R samples were separated from the D20 and DT samples in principle component 1 (62.7%), whereas D75 was separated from F and R, and D20 from DT, in principle component 2 (11.7%).
Clustering of DEGs by expression profiles and their correlations with physiological changes during dehydration
The hierarchical clustering of all DEGs in the five treatments showed significant differences in gene expression profiles. Therefore, we defined seven clusters of DEGs, along with a representative curve describing each transcript accumulation pattern (Fig. 3a, b). Transcripts peaking at the DT, D20 and D75 stages were represented as Clusters 1, 2 and 7, respectively, whereas transcripts repressed during dehydration (D75, D20 and DT) were found in Clusters 3, 4, 5 and 6. Interestingly, most dehydration-repressed transcripts reached their lowest levels at the D20 stage, regardless of whether they returned to control levels after rehydration (Cluster 5), or were transiently up-regulated at certain stages, such as D75 (Cluster 4), DT (Cluster 3) and R (Cluster 6).
The GO term enrichment of each cluster indicates that both metabolic processes (catabolism and biosynthesis) and developmental growth were repressed upon slight dehydration (Cluster 3), whereas responses to stimuli and redox processes/homeostasis were actively maintained until the water content decreased to 20% (Clusters 2, 4 and 5) (Fig. 3c). On the other hand, the transcripts activated at D75 (Cluster 7, translation and gene expression related) differed markedly from those activated at D20 and DT (Clusters 2 and 1, catabolic process, response to stimulus, redox process/homeostasis, localization, and regulation of seed germination).
To understand the significance of these gene expression clusters and to test the robustness of our cluster analysis, the first principle component from each stage was subjected to correlation analysis with four previously published physiological parameters . We found that the decreases in net CO2 assimilation, stomatal conductance, transpiration and Fv/Fm during dehydration showed negative correlations with Clusters 1 and 2 and positive correlations with Clusters 5 and 6 (Fig. 3d).
GO and KEGG enrichment of DEGs in dehydrated plants
We performed KEGG analysis using the unigenes from each stage and observed the most significant enrichments in plant hormone signal transduction, spliceosome, starch and sucrose metabolism, and ubiquitin-mediated proteolysis pathways (Fig. 4a). Further GO and KEGG enrichment analyses were performed using the up- or down-regulated transcripts from each stage (Additional file 7: Table S5, Additional file 5: Fig. S2b). The GO terms electron transporter, chaperone binding, and iron-sulfur cluster binding were enriched among the transcripts induced by the D75, D20 and DT treatments, respectively, whereas developmental processes were enriched among the repressed transcripts during all dehydration stages (Additional file 7: Table S5). Metabolic pathways were also enriched among the transcripts that were repressed at all dehydration stages (Additional file 5: Fig. S2b). Interestingly, photosynthesis was enriched among not only the D20- and DT-repressed genes but also the D75-induced genes; plant hormone signal transduction was enriched among not only the transcripts that were decreased in all dehydrated samples but also in D20-induced transcripts. In addition, ubiquitin-mediated proteolysis, regulation of autophagy and protein processing in the endoplasmic reticulum (ER) were enriched among the transcripts induced by both D20 and DT. Notably, a large number of D20-specific DEGs associated with these pathways was not found in moderately dehydrated plants  (Additional file 2: Fig. S1d). This result prompted us to speculate that D20 (severe dehydration) is a crucial stage in transcriptional regulation of desiccation tolerance in H. rhodopensis.
Finally, we used these changes to construct a functionally organized GO/pathway term network based on the connectivity among GO terms (Fig. 4b, Additional file 8: Table S6). Numerous changes in this network were analyzed further, including primary and secondary metabolism, photosynthesis, hormone signal transduction, protein quality control and DNA repair.
DEGs related to energy and primary and secondary metabolism
The expression patterns of photosynthesis-related genes generally were in agreement with the inhibition of photosynthesis observed during D20 and DT, including genes encoding the subunits of photosystem I (PSI) and photosystem II (PSII), light-harvesting complexes II and I, cytochrome b6/f complex and electron transport. Approximately a quarter of these transcripts, including NADPH-quinone oxidoreductase subunits and ATP synthase, were up-regulated only during the early stage of dehydration (D75). However, a few transcripts encoding ferredoxin 3 (FDX3) and subunits of the PSI and PSII reaction centers (psaA and psbC) maintained elevated expression levels throughout the dehydration process (Additional file 9: Fig. S3).
The expression of genes related to metabolism was analyzed and closely compared with our previous metabolomics data [10,11,12]. A map was constructed of the D20-enriched KEGG pathways related to primary and secondary metabolism and the corresponding edges (significant co-expression relationship between two genes) (Additional file 10: Fig. S4). We found that the expression of phospholipase D (PLD), glycerol-3-phosphate dehydrogenase [NAD(+)], and manganese-dependent ADP-ribose/CDP-alcohol diphosphatase in the glycerophospholipid metabolic pathway was increased during desiccation. Two key enzymes involved in fatty acid degradation, 2-methylacyl-CoA dehydrogenase (MBCD) and acyl-CoA oxidase 2 (ACOX2), and enzymes involved in lipid degradation were also up-regulated. Branched-chain amino acid aminotransferase 3 (BCAL3) and hydroxymethylglutaryl-CoA lyase, involved in the valine, leucine and isoleucine degradation pathways, citrate synthase (CYSZ) involved in the citric acid (TCA) cycle, and enzymes involved in glycolysis/gluconeogenesis, such as triosephosphate isomerase (TPIC), phosphoglycerate kinase (PGKH) and hexokinase-1 (HXK1), were up-regulated, with the exception of the non-key enzyme pyrophosphate-fructose 6-phosphate 1-phosphotransferase subunit beta. Enzymes involved in sucrose synthesis, such as acid beta-fructofuranosidase and sucrose-phosphate synthase (SPSA1), were increased during dehydration, whereas sucrose synthase (SUS2) was elevated only in the D20 and DT treatments (Fig. 5a).
With respect to secondary metabolism, flavonoid 3′-monooxygenase (F3’H) and SRG1, which are involved in the accumulation of flavonoids, and flavanone 7-O-glucoside 2”-O-beta-L-rhamnosyltransferase (I7GT1), from the phenylpropanoid metabolism pathway, were up-regulated during dehydration. Two enzymes from the glutathione cycle, glutathione reductase (GSHRC) and phospholipid hydroperoxide glutathione peroxidase (GPX4), were also increased during desiccation (Fig. 5b).
Positive or negative correlations were found between the expression levels of most genes and the previously measured accumulation of their corresponding metabolites, depending on the reaction type (catabolic or anabolic) (Fig. 5c). For example, the levels of L-valine and leucine were negatively correlated with the expression of BCAL3, which plays a role in the catabolism of these amino acids, whereas SUS2 expression was correlated positively with sucrose content but negatively with UDP-glucose (Fig. 5c).
DEGs involved in plant hormone signaling pathways
Transcripts encoding PYR/PYL/RCAR-type ABA receptors (PYL1 and PYL4), which are key factors (ABF2, ABF3, DPBF3/AREB3, GBF4 and AIB) in the ABA signaling pathway, and NPR1, a key factor in the SA signaling pathway, were highly induced during the D20 and DT stages (Fig. 6). Brassinosteroid (BR) receptors (BRL1 and BRL3) and BZR1 (BEH4), key factors in the BR signaling pathway, a JA receptor (COI1), an ethylene receptor (EIN4) and key regulators (CTR1, EDR1, EIN3 and EIL1/EIN3-like) in the ethylene signaling pathway, a gibberellic acid (GA) receptor (GID1B), an auxin receptor (TIR1), IAA-acyl acid-amido synthetases (GH3.3 and GH3.6), and negative regulators (IAA27, IAA9 and IAA13) in the auxin signaling pathway were all induced during dehydration. In contrast, positive regulators of the cytokinin signaling pathway AHP1 and AHP4 were up-regulated only slightly during the early stage of dehydration.
DEGs related to protein quality control and DNA repair
Protein processing in the ER and ubiquitin-mediated proteolysis are both involved in protein-quality control. Among the related DEGs, most small heat shock proteins and HSP/HSC70s were most strongly up-regulated in the D20 treatment, whereas stromal HSP70c in chloroplasts was reduced in D20 and DT after transiently increasing in D75 (Fig. 7a). ER lumen protein-retaining receptor 2 (ERD2) and ER oxidoreductin-1 (ERO1) showed prominent up-regulation in the D20 and DT stages. The E3 ubiquitin-protein ligases, including the C terminus of HSC70-interacting protein, Rab interactor 2, RING membrane-anchor 1 homolog 1, and RING membrane-anchor 2 and 3, showed overall up-regulation from D75 onward, reaching the highest levels in D20 or DT.
Base excision repair (BER), nucleotide excision repair (NER) and mismatch repair (MMR) are the three main pathways involved in DNA repair. The majority of transcripts encoding proteins that participate in NER, including DNA damage binding proteins 1 and 2 (DDB1/DDB2), decrease in DNA methylation (DDM1), replication protein A 1 and 2 (RPA1/RPA2), and DNA repair helicase (XPB1) showed prominent patterns of up-regulation in D20 and DT (Fig. 7b). In contrast, most transcripts involved in BER and MMR were down-regulated during dehydration (Fig. 7b).
Most plants are not of the resurrecting type and therefore can tolerate water loss down to approximately 30% RWC only for relatively short periods [11, 15]. Below this water content, they are unable to recover even if water becomes available again; in other words, this level is the “point of no return” for most plants. Our previous physiological and metabolic studies described numerous changes that occurred in H. rhodopensis during severe dehydration [10, 12]. One of the aims of our present study was to understand whether and when (in terms of before or around this stage of dehydration) alterations occur at the transcriptome level in association with respective physiological changes.
It was previously shown that signaling and transcriptional changes induced by dehydration in H. rhodopensis regulate growth and photosynthesis (protein kinases and early light-inducible proteins), sugar metabolism (sucrose and raffinose synthase) and activation of diverse protectants (catalases, LEA and HSPs) . Our de novo analysis not only supported these results but also identified novel genes that extend our understanding of the roles that cyclic electron flow (CEF), carbon turnover, phenylpropanoids, stress hormone signal transduction, protein quality control and DNA repair play in desiccation tolerance. The additional dehydration time points used here allowed us to outline the dynamic transcriptomic regulation during the dehydration and rehydration processes. The largest numbers of both up- and down-regulated DEGs were observed at this stage (Fig. 1). Most dehydration-repressed transcripts reached their lowest level at D20, regardless of whether these genes were transiently up-regulated during certain stages of dehydration or upon recovery (Fig. 3). A significant portion of the D20 up-regulated DEGs were annotated in the same biological processes and KEGG pathways as those up-regulated in DT, but did not always match those up-regulated in D75, as a large proportion of the DEGs showed opposing up- and down-regulation when the expression patterns of D75 and D20 were compared (Fig. 3). These data clearly indicate that desiccation tolerance mechanisms undergo reprogramming at the transcriptome level at the D20 stage. This finding extends the present understanding of desiccation tolerance mechanisms in H. rhodopensis.
Energy and primary and secondary metabolic processes
Over-reduction of the electron transport chain and carbon starvation are among the main negative impacts of drought stress that result in reduced carbon fixation. Approximately 40% of the novel unigenes reported here are related to metabolic pathways and photosynthesis (Additional file 2: Fig. S1d). This large proportion allowed us to investigate the correlations of our transcriptomic data with physiological and metabolomic data (Fig. 3d; 5c). In this analysis, we identified changes in the photosynthetic apparatus and the mobilization of free and lipid-derived fatty acids and branched amino acids to ensure the availability of energy and alternative carbon sources to allow accumulation of osmoprotectants and antioxidants during desiccation. Many of these changes are similar to those that have been observed in plant seeds or during carbon starvation in non-resurrection plants [10, 16,17,18,19,20,21,22].
Resurrection plants can reduce linear electron transport flux through reorganization of the photosynthetic apparatus, thus preventing oxidative stress [10, 23]. A dehydration-triggered decrease in PSII performance was accompanied by reduced expression of associated genes, which belonged mostly to clusters 4, 5 and 6 (Fig. 3d), including chlorophyll a/b-binding protein and the light-harvesting complex Lhcb1, as reported previously . In addition, the decrease in transcription of the PSI-K subunit of PSI (Additional file 9: Fig. S3) is directly responsible for the altered distribution of excitation energy and uncoupling between the photosystems , possibly explaining the observed decrease in electron flow between PSI and PSII [10, 25]. On the other hand, our present transcriptomic evidences indicating maintenance of CEF during dehydration corroborated the results of our previous study . The up-regulation of ferredoxin and NADPH-quinone oxidoreductase subunits not only suggests the presence of the chloroplastic NAD(P)H dehydrogenase complex in H. rhodopensis but also highlights its role in desiccation tolerance as a component of CEF [26, 27].
The process of lipid turnover is very important in both tolerant and sensitive plants under drought stress. The up-regulated expression of PLD during dehydration (Fig. 5a) is likely to correspond with increased glycerophosphorylethanolamine (GPE) and glycerophosphorylcholine (GPC) and decreased phosphatidylethanolamine levels, reported in our previous study . PLD activity induction by dehydration has been reported in the resurrection plant Craterostigma plantagineum , as well as in the drought-tolerant Arabidopsis ecotype Columbia . A drastic decrease in lipid levels, especially membrane lipids such as phosphatidylethanolamine, has been observed in Ramonda serbica, a close relative of H. rhodopensis  and in the African resurrection plant Sporobolus stapfianus . Thus, the reported correlation between PLD and related metabolites  (Fig. 5c), combined with the increased levels of malondialdehyde [5, 11, 30], strongly confirm the occurrence of lipid degradation in H. rhodopensis under stress. In comparison with the massive degradation of membranes in desiccation-sensitive species , these processes are quite limited in resurrection plants. Here, the turnover of damaged lipids may have a protective role, ensuring the accumulation of osmoprotectants such as GPC , signaling molecules and carbon precursors needed for biosynthetic processes during drought stress [10, 19, 33]. The up-regulation of the main enzymes responsible for fatty-acid degradation, MBCD and ACOX2 (Fig. 5a), supports the decreased free fatty acid content found during desiccation (Fig. 5c) .
Degradation of branched-chain amino acids provides an alternative carbon source in plants under unfavorable conditions [17, 21]. The up-regulation of BCAL3 and HMGL, which are involved in the degradation of valine, leucine and isoleucine (Fig. 5a), may explain the decreasing levels of these amino acids during dehydration (Fig. 5c) [10, 12].
Increased expression of CYSZ, TPIC, PGKH, HXK1, SUS2 and SPSA1 (Fig. 5a), which are related to the citric acid cycle, glycolysis/gluconeogenesis and sucrose biosynthesis, is in agreement with previous reports of increased metabolite levels during these biological processes and massive accumulation of sucrose during dehydration (Fig. 5c) [5, 10, 12].
We also found significantly elevated expression of genes involved in the synthesis of phenylpropanoids and flavonoids (Fig. 5b). I7GT1 encodes a UDP-glycosyltransferase, a type of enzyme that is related to the accumulation of p-coumaryl alcohol 4-O-glucoside, coniferin, syringin, and coniferyl and sinapyl alcohol 4-O-glucoside via glucosylation of soluble intermediates in Arabidopsis . The up-regulation of F3’H and SRG1 during dehydration (Fig. 5b) is correlated directly with the accumulation of flavonoids, including quercetin, aromadendrin, taxifolin, ephedrine, leucocyanidin, cyanidin, cis-3,4-leucopelargonidin and pelargonidin [35, 36]. These phenylpropanoids and flavonoids are common plant secondary metabolites that function both as structural and signaling molecules .
Glutathione is involved in plant desiccation tolerance as an antioxidant [11, 37, 38]. DEGs related to glutathione metabolism have been reported previously in H. rhodopensis  and another resurrection plant species from the same family Gesneriaceae, Boea hygrometrica . Over-expression of two transcripts in this study (Fig. 5b), GPX4 and GSHRC, accorded with these data and strongly confirmed the importance of glutathione as a major factor in desiccation tolerance.
Plant hormone signal transduction
The plant hormone signaling network plays an integral role in the perception of and response to unfavorable environments. Previously, transient induction of JA and ABA has been reported, reaching their highest levels at approximately 55 and 25% RWC, respectively, in the dehydrated leaves of H. rhodopensis, whereas SA accumulates to its highest level at 25% RWC and is maintained during full dehydration . The increased expression of biosynthesis- and signaling-related genes for these stress hormones contributes to plant drought tolerance by inducing biosynthesis of defense proteins and protective secondary metabolites. In this study, we found that the JA receptor COI1 was induced during dehydration (Fig. 6), but that the transcriptional activator of the JA signaling pathway, JAR1, was significantly repressed at the D20 and DT stages, indicating conversion of the JA signaling pathway under continuous stress, in accordance with the reported pattern of JA dynamics . In contrast, at the transcriptional level, ABA and SA signaling transduction is activated at D20 and DT through the induction of key transcription factors such as ABFs, GBFs and AIB factors, and NPR1, thus confirming the roles of ABA and SA as stress signaling molecules and that of SA in the antioxidant response of the plant to dehydration .
Dynamic changes of ethylene in desiccation tolerance have not been reported; however, 1-aminocyclopropane-1-carboxylic acid oxidase transcripts have been shown to accumulate during dehydration . In this study, the increased transcription of ethylene receptors (ETR1, ETR2 and EIN4) and the key transcription factor EIN3/EIL1 during the early stage of dehydration suggests a possible role of ethylene in the early response to dehydration, which is followed by repression of signal transduction in H. rhodopensis, as indicated by up-regulation of the key negative regulator CTR1 at the D20 and DT stages. Ethylene negatively regulates cell proliferation at the root apical meristem . Thus, the transient increase of ethylene biosynthesis- and signaling-related genes in H. rhodopensis may contribute to the rapid inhibition of cell expansion during dehydration.
Crosstalk among SA, JA and ethylene-dependent signaling pathways regulates plant responses to abiotic and biotic stresses, which are frequently associated with hypersensitive response (HR)-like cell death [41,42,43]. The final response depends on the extent of crosstalk among these molecules. Inhibition of the ethylene and JA signaling pathways under severe and full dehydration and the crosstalk of these pathways with SA signal transduction likely contribute to prevention of HR-associated cell death, which is required for desiccation tolerance in resurrection plants.
Auxins, cytokinins, GAs and BRs have been recognized as crucial signaling molecules that control plant growth and development. The positive regulators of the cytokinin signaling pathway AHP1 and AHP4 were only slightly induced in the early stage of dehydration, whereas negative regulators of auxin pathways were induced throughout the dehydration process (Fig. 6), indicating inhibition of both signal transduction processes. This finding is in agreement with the observation of suspended growth of plants during dehydration. Notably, receptor genes for auxin, GA and BRs, including TIR1, GID1B, BRL1 and BRL3, were up-regulated throughout the dehydration and rehydration processes, in parallel with IAA-acyl acid amido conjugators (group II GH3s)  and the BZR1 homolog BEH2, indicating a possible mechanism for the rapid recovery of the resurrection plant upon rehydration.
Protein-quality control and DNA repair
The ability to maintain protein and DNA integrity is indispensable to cellular survival under constant exposure to adverse stresses that may result in mutagenesis or cell death . Despite this importance, data on protein quality control and DNA repair are very limited or lacking, although recent molecular studies show their potential involvement in desiccation tolerance in the budding yeast Saccharomyces cerevisiae, Arabidopsis seeds and the resurrection plant B. hygrometrica [46,47,48]. The mechanisms controlling these processes, such as the unfolded protein response, aid in protein folding or in degradation of misfolded secretory proteins by producing protein-folding and other factors [49, 50]. Activation of the unfolded protein response in H. rhodopensis was reflected by up-regulation of two key regulators during D20 and DT (Fig. 7a): ERD2, a transmembrane ER lumen protein-retaining receptor that controls traffic in the Golgi body and retrograde transport to reclaim ER proteins , and ERO1, a protein associated with protein disulfide isomerase oxidoreductases, which play a key role in ER redox regulation via protein stabilization by forming disulfide bonds between structures .
The ubiquitin–proteasome system is the major pathway for proteasomal degradation of damaged and misfolded proteins. Activation of this system in D20 and DT, indicated by up-regulation of the majority of genes encoding E3 ubiquitin ligase (Fig. 7a), hypothetically results in proteomic redistribution; however, the reasons for this redistribution are unclear and remain to be explored in the future. A regulatory role of protein ubiquitination has been revealed, based on the accumulation of important metabolites during acquisition of rapid desiccation tolerance in B. hygrometrica . E3 ubiquitin ligase may cooperatively interact with HSP/C70 for quality control [54, 55]. This interaction may occur intensively in H. rhodopensis during the D20 and DT stages.
Damaged DNA can be repaired via the processes of BER, MMR, and NER. In this study, we observed predominate activation of NER pathway genes, such as DDB1/DDB2, DDM1, RPA1/RPA2, XPB1 and PCNAs , indicating the importance of proofreading DNA repair to maintain genome integrity during severe dehydration. Studies of the DNA-damage response in seeds have suggested that particular proteins are required for certain repair mechanisms during the early stage of desiccation [47, 57]. Considering the similarity of desiccation tolerance in vegetative tissues of resurrection plants and seeds of non-recalcitrant species , a dynamic mechanism for maintaining genome integrity is therefore reasonably predicted in H. rhodopensis under desiccation. NER has been found to play an important role in reversing UV-induced DNA damage in plants . Building upon known protective responses such as LEA and HSPs , this is the first report to associate NER with desiccation tolerance in resurrection plants. Further studies are needed to clarify the detailed mechanism of the DNA-repair process in H. rhodopensis under desiccation conditions.
The underlying difference between resurrection and non-resurrection plants under extreme drought remains somewhat unclear. We believe that our transcriptomic data contribute to the elucidation of new desiccation tolerance mechanisms in the Balkan endemic plant H. rhodopensis, as we have identified a large number of novel dehydration-responsive genes (Fig. 8). During the early stages of stress, the increased transcript levels of stress hormone signaling pathways initiate crosstalk among these signaling molecules, which is followed by fine tuning to avoid an increase in hypersensitive reaction-induced cell death as dehydration becomes severe. Consequently, further protective mechanisms related to energy, primary metabolism and secondary metabolism, protein quality control and DNA repair were transcriptionally mobilized during this stage. Activation of sucrose synthesis, lipid and fatty-acid turnover, HSP levels, LEA, NER, and the flavonoid, phenylpropanoid and glutathione pathways are among the main transcriptomic changes observed during desiccation. Notably, CEF and the mobilization of reserve carbon sources are maintained throughout the late stages of desiccation, in contrast to the linear decreases in electron transport and carbon uptake.
Together, our results show that during severe dehydration, a reprograming of the transcriptome occurs that defines the desiccation tolerance and resurrection behavior of H. rhodopensis.
Materials and methods
Plant materials and stress treatments
In vitro-propagated plants of H. rhodopensis  were transferred to a standard soil–sand–gravel mixture in pots (6 cm diameter) and grown for approximately 1 year in a controlled environment at 22–24 °C, under a 16-h photoperiod, 60% relative humidity, and a photon flux density of 36 μmol m− 2 s− 1.
Dehydration was induced by withholding water from potted plants. Plants required at least 8 days to become fully desiccated. After 10 days in the fully dry state, the plants were re-watered. Rehydration was initiated by watering the plants, and full recovery took 1 week. A non-destructive method of monitoring water content based on in vivo fast fluorescence was applied to collect plant leaf samples, according to the established relationship between leaf water content and photosynthetic performance during stress and recovery . The Handy PEA (Hansatech Instrument Ltd., King’s Lynn, Norfolk, UK) device was used for measurement of the prompt chlorophyll a fluorescence of plants in the light-adapted state. Based on the SOM (self-organizing map) developed previously , we were able to collect samples at five time points corresponding to the following RWC stages: F (90% RWC), D75 (75% RWC), D20 (20% RWC), DT (6% RWC) and R (90% RWC). Mature, fully expanded leaves from the middle of rosettes of similar size and appearance were used in all measurements. At least four independent plants were used for each sampling point.
RNA extraction, cDNA library construction, and Illumina high-throughput sequencing
Plant leaves (0.4 g) were ground under liquid nitrogen to a fine powder using a cooled mortar and pestle. Samples were placed in 2-mL Eppendorf tubes, and total RNA was isolated following the protocol for the GeneMATRIX Universal RNA kit (EURx Ltd., Gdansk, Poland).
Extracted RNA was used to construct cDNA libraries for sequencing analysis. Ten cDNA libraries were prepared using the Truseq™ RNA sample preparation kit from Illumina (San Diego, CA, USA). Libraries were size-selected for cDNA target fragments of 200–300 bp on 2% low-range ultra-agarose gels followed by PCR amplification using Phusion DNA polymerase (NEB, Ipswich, MA, USA) for 15 cycles. After quantification using the TBS380 fluorometer, 150-bp paired-end RNA-seq libraries were constructed using the Illumina HiSeq 4000 platform (Illumina).
Transcriptome de novo assembly, sequence annotation, classification and alignment
The raw data were first processed using Seqprep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle) and the obtained high quality clean reads were subjected to de novo assembly using Trinity (https://github.com/trinityrnaseq/trinityrnaseq) . The unique assembled transcripts were then used for functional annotation using BLASTX and BLAST2GO (http://www.blast2go.com/b2ghome) . Metabolic pathway analysis was performed using KEGG (http://www.genome.jp/kegg) . Sequence alignment was performed using BLAT  with contig sequences from both this de novo assembly and published data  (doi: https://doi.org/10.1007/s00018-012-1155-6, 18_2012_1155_MOESM1_ESM.xlsx). Unigenes were considered to differ when the sequence homology was < 90% in identity and < 80% in coverage.
Differential expression analysis and functional enrichment
To identify DEGs between two different samples, the expression level of each transcript was calculated according to the FPKM method. RSEM (http://deweylab.biostat.wisc.edu/rsem/)  was used to quantify gene and isoform abundances. The R statistical software package Empirical analysis of Digital Gene Expression in R (http://www.bioconductor.org/packages/2.4/bioc/html/edgeR.html)  was utilized for differential gene expression analysis. Principle component analysis of DEGs was performed using SIMCA 13.0 (Umetrics, Satorius Stedim Biotech, Umeå, Sweden). In addition, functional enrichment analysis including GO and KEGG terms in comparison with the whole-transcriptome background was performed using a Bonferroni-corrected p-value < 0.05 threshold. GO functional enrichment and KEGG pathway enrichment were carried out using Goatools (https://github.com/tanghaibao/Goatools) and KOBAS (http://kobas.cbi.pku.edu.cn/download.php) . For visualization of the functionally grouped network of KEGG terms, the ClueGO plugin  was used with Cytoscape . To map all identified genes with altered expression in D20 and to visualize the pathway relationships among these genes, we used a general mapping method in which mapped genes represent a low percentage of all associated genes per pathway, and Bonferroni correction was omitted. Subsequently, the KEGG file for each pathway was downloaded, and all genes in the pathway were extracted from the ClueGO table. Then, network analysis was performed using Cytoscape software. The Pearson correlations between genes and metabolites were calculated using Cytoscape. The expression of each gene extracted from the ClueGO table for each pathway is expressed as the average log2 ratio for each treatment compared with the control.
All purified RNA samples described above were adjusted to the same concentration (100 ng/μL) using RNase-free ddH2O and then reverse-transcribed into cDNA using M-MLV reverse transcriptase (Promega, Madison, WI, USA), following the manufacturer’s instructions. Quantification of the transcripts was performed using SYBR® Green Realtime PCR Master Mix (TOYOBO, Osaka, Japan). Each reaction (10 μL) contained 0.4 μL of the cDNA template, 0.8 μL of the primer mixture (10 μM each, mixed to equivalent volume), 5.0 μL 2× PCR Master Mix and 3.8 μL ddH2O. Reactions were performed using a two-step method: 95 °C for 10 s, then 40 cycles of 95 °C for 10 s, 55 °C for 30 s, 72 °C for 30 s, and 95 °C for 15 s, followed by a final melt cycle from 55 °C to 95 °C. Gene-specific primers are shown in Additional file 6: Table S4. qPCR was performed for three technical replicates of each RNA sample (as well as two independent biological replicates). The normalized relative expression levels of the target genes were calculated by the 2-ΔΔCt (ΔCt = CtTarget – Ct18S, where Ct is the cycle threshold) using 18S rRNA as the internal standard.
The clustering of log-transformed expression data from each treatment was performed using Gene Cluster 3.0, with normalization of genes, the centered correlation similarity metric and complete linkage method. Clustered data were viewed, and the corresponding heatmaps were generated using JAVA TREEVIEW 1.6.4. Genes of interest were selected according to their relevance in particular biological processes or pathways, and then their expression data were pooled, log2 transformed and normalized to the untreated hydrated controls.
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Nucleotide excision repair
Principal component analysis
Costa MD, Artur MA, Maia J, Jonkheer E, Derks MF, Nijveen H, Williams B, Mundree SG, Jimenez-Gomez JM, Hesselink T, et al. A footprint of desiccation tolerance in the genome of Xerophyta viscosa. Nat Plants. 2017;3:17038.
Rodriguez MC, Edsgard D, Hussain SS, Alquezar D, Rasmussen M, Gilbert T, Nielsen BH, Bartels D, Mundy J. Transcriptomes of the desiccation-tolerant resurrection plant Craterostigma plantagineum. Plant J. 2010;63(2):212–28.
Xiao L, Yang G, Zhang L, Yang X, Zhao S, Ji Z, Zhou Q, Hu M, Wang Y, Chen M, et al. The resurrection genome of Boea hygrometrica: a blueprint for survival of dehydration. Proc Natl Acad Sci U S A. 2015;112(18):5833–7.
Yobi A, Schlauch KA, Tillett RL, Yim WC, Espinoza C, Wone BW, Cushman JC, Oliver MJ. Sporobolus stapfianus: insights into desiccation tolerance in the resurrection grasses from linking transcriptomics to metabolomics. BMC Plant Biol. 2017;17(1):67.
Gechev TS, Benina M, Obata T, Tohge T, Sujeeth N, Minkov I, Hille J, Temanni MR, Marriott AS, Bergstrom E, et al. Molecular mechanisms of desiccation tolerance in the resurrection glacial relic Haberlea rhodopensis. Cell Mol Life Sci. 2013;70(4):689–709.
Georgieva T, Christov NK, Djilianov D. Identification of desiccation-regulated genes by cDNA-AFLP in Haberlea rhodopensis: a resurrection plant. Acta Physiol Plant. 2012;34(3):1055–66.
Ivanova Z, Sablok G, Daskalova E, Zahmanova G, Apostolova E, Yahubyan G, Baev V. Chloroplast genome analysis of resurrection tertiary relict Haberlea rhodopensis highlights genes important for desiccation stress response. Front Plant Sci. 2017;8:204.
Georgieva K, Maslenkova L, Peeva V, Markovska Y, Stefanov D, Tuba Z. Comparative study on the changes in photosynthetic activity of the homoiochlorophyllous desiccation-tolerant Haberlea rhodopensis and desiccation-sensitive spinach leaves during desiccation and rehydration. Photosynth Res. 2005;85(2):191–203.
Djilianov DL, Dobrev PI, Moyankova DP, Vankova R, Georgieva DT, Gajdošová S, Motyka V. Dynamics of endogenous phytohormones during desiccation and recovery of the resurrection planspecies Haberlea rhodopensis. J Plant Growth Regul. 2013;32(3):564–74.
Mladenov P, Finazzi G, Bligny R, Moyankova D, Zasheva D, Boisson AM, Brugiere S, Krasteva V, Alipieva K, Simova S, et al. In vivo spectroscopy and NMR metabolite fingerprinting approaches to connect the dynamics of photosynthetic and metabolic phenotypes in resurrection plant Haberlea rhodopensis during desiccation and recovery. Front Plant Sci. 2015;6:564.
Djilianov D, Ivanov S, Moyankova D, Miteva L, Kirova E, Alexieva V, Joudi M, Peshev D, Van den Ende W. Sugar ratios, glutathione redox status and phenols in the resurrection species Haberlea rhodopensis and the closely related non-resurrection species Chirita eberhardtii. Plant Biol. 2011;13(5):767–76.
Moyankova D, Mladenov P, Berkov S, Peshev D, Georgieva D, Djilianov D. Metabolic profiling of the resurrection plant Haberlea rhodopensis during desiccation and recovery. Physiol Plantarum. 2014;152(4):675–87.
Wang L, Yu S, Tong C, Zhao Y, Liu Y, Song C, Zhang Y, Zhang X, Wang Y, Hua W. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. Genome Biol. 2014;15(2):R39.
Georgieva K, Sarvari E, Keresztes A. Protection of thylakoids against combined light and drought by a lumenal substance in the resurrection plant Haberlea rhodopensis. Ann Bot. 2010;105(1):117–26.
Oliver MJ, Guo L, Alexander DC, Ryals JA, Wone BW, Cushman JC. A sister group contrast using untargeted global metabolomic analysis delineates the biochemical regulation underlying desiccation tolerance in Sporobolus stapfianus. Plant Cell. 2011;23(4):1231–48.
Avelange-macherel MH, Ly-Vu B, Delaunay J, Richomme P, Leprince O. NMR metabolite profiling analysis reveals changes in phospholipid metabolism associated with the re-establishment of desiccation tolerance upon osmotic stress in germinated radicles of cucumber. Plant Cell Environ. 2006;29(4):471–82.
Diebold R, Schuster J, Däschner K, Binder S. The branched-chain amino acid transaminase gene family in Arabidopsis encodes plastid and mitochondrial proteins. Plant Physiol. 2002;129(2):540–50.
Dieuaide M, Couée I, Pradet A, Raymond P. Effects of glucose starvation on the oxidation of fatty acids by maize root tip mitochondria and peroxisomes: evidence for mitochondrial fatty acid β-oxidation and acyl-CoA dehydrogenase activity in a higher plant. Biochem J. 1993;296(1):199–207.
Fan J, Yu L, Xu C. A central role for triacylglycerol in membrane lipid breakdown, fatty acid β-oxidation, and plant survival under extended darkness. Plant Physiol. 2017;174(3):1517–30.
Leprince O, Hoekstra FA. The responses of cytochrome redox state and energy metabolism to dehydration support a role for cytoplasmic viscosity in desiccation tolerance. Plant Physiol. 1998;118(4):1253–64.
Schuster J, Binder S. The mitochondrial branched-chain aminotransferase (AtBCAT-1) is capable to initiate degradation of leucine, isoleucine and valine in almost all tissues in Arabidopsis thaliana. Plant Mol Biol. 2005;57(2):241–54.
Tarancón C, González-Grandío E, Oliveros JC, Nicolas M, Cubas P. A conserved carbon starvation response underlies bud dormancy in woody andherbaceous species. Front Plant Sci. 2017;8:788.
Challabathula D, Zhang Q, Bartels D. Protection of photosynthesis in desiccation-tolerant resurrection plants. J Plant Physiol. 2018;227:84–92.
Jensen PE, Gilpin M, Knoetzel J, Scheller HV. The PSI-K subunit of photosystem I is involved in the interaction between light-harvesting complex I and the photosystem I reaction center core. J Biol Chem. 2000;275(32):24701–8.
Goltsev V, Zaharieva I, Chernev P, Kouzmanova M, Kalaji HM, Yordanov I, Krasteva V, Alexandrov V, Stefanov D, Allakhverdiev SI. Drought-induced modifications of photosynthetic electron transport in intact leaves: analysis and use of neural networks as a tool for a rapid non-invasive estimation. Biochim Biophys Acta-Mol Cell Res. 2012;1817(8):1490–8.
Iwai M, Takizawa K, Tokutsu R, Okamuro A, Takahashi Y, Minagawa J. Isolation of the elusive supercomplex that drives cyclic electron flow in photosynthesis. Nature. 2010;464(7292):1210–3.
Suorsa M. Cyclic electron flow provides acclimatory plasticity for the photosynthetic machinery under various environmental conditions and developmental stages. Front Plant Sci. 2015;6:800.
Frank W, Munnik T, Kerkmann K, Salamini F, Bartels D. Water deficit triggers phospholipase D activity in the resurrection plant Craterostigma plantagineum. Plant Cell. 2000;12(1):111–23.
Gigon A, Matos A-R, Laffray D, Zuily-Fodil Y, Pham-Thi A-T. Effect of drought stress on lipid metabolism in the leaves of Arabidopsis thaliana (ecotype Columbia). Ann Bot. 2004;94(3):345–51.
Solti Á, Mihailova G, Sárvári É, Georgieva K. Antioxidative defence mechanisms contributes to desiccation tolerance in Haberlea rhodopensis population naturally exposed to high irradiation. Acta Biologica Szegediensis. 2014;58(1):11–4.
Li A, Wang D, Yu B, Yu X, Li W. Maintenance or collapse: responses of extraplastidic membrane lipid composition to desiccation in the resurrection plant Paraisometrum mileense. PLoS One. 2014;9(7):e103430.
Gallazzini M, Ferraris JD, Burg MB. GDPD5 is a glycerophosphocholine phosphodiesterase that osmotically regulates the osmoprotective organic osmolyte GPC. Proc Natl Acad Sci U S A. 2008;105(31):11026–31.
Rakic T, Lazarevic M, Jovanovic ZS, Radovic S, Siljak-Yakovlev S, Stevanovic B, Stevanovic V. Resurrection plants of the genus Ramonda: prospective survival strategies-unlock further capacity of adaptation, or embark on the path of evolution? Front Plant Sci. 2014;4:550.
Lanot A, Hodge D, Lim E-K, Vaistij FE, Bowles DJ. Redirection of flux through the phenylpropanoid pathway by increased glucosylation of soluble intermediates. Planta. 2008;228(4):609–16.
Saito Y. Bershad, Brian N, levy, Henry M. Manageability, availability and performance in Porcupine: a highly scalable, cluster-based mail service Acm Sigops Operating Systems Review. 1999;34(2):9–11.
Schoenbohm C, Martens S, Eder C, Forkmann G, Weisshaar B. Identification of the Arabidopsis thaliana flavonoid 3′-hydroxylase gene and functional expression of the encoded P450 enzyme. Biol Chem. 2000;381(8):749–53.
Georgieva K, Rapparini F, Bertazza G, Mihailova G, Sárvári É, Solti Á, Keresztes Á. Alterations in the sugar metabolism and in the vacuolar system of mesophyll cells contribute to the desiccation tolerance of Haberlea rhodopensis ecotypes. Protoplasma. 2017;254(1):193–201.
Toldi O, Tuba Z, Scott P. Vegetative desiccation tolerance: is it a goldmine for bioengineering crops? Plant Sci. 2009;176(2):187–99.
Wang L, Shang H, Liu Y, Zheng M, Wu R, Phillips J, Bartels D, Deng X. A role for a cell wall localized glycine-rich protein in dehydration and rehydration of the resurrection plant Boea hygrometrica. Plant Biol. 2009;11(6):837–48.
Street IH, Aman S, Zubo Y, Ramzan A, Wang X, Shakeel SN, Kieber JJ, Schaller GE. Ethylene inhibits cell proliferation of the Arabidopsis root meristem. Plant Physiol. 2015;169(1):338–50.
Clarke JD, Volko SM, Ledford H, Ausubel FM, Dong X. Roles of salicylic acid, jasmonic acid, and ethylene in cpr-induced resistance in Arabidopsis. Plant Cell. 2000;12(11):2175–90.
Dangl JL, Dietrich RA, Richberg MH. Death don't have no mercy: cell death programs in plant-microbe interactions. Plant Cell. 1996;8(10):1793–807.
Greenberg JT. Programmed cell death in plant-pathogen interactions. Annu Rev Plant Biol. 1997;48(1):525–45.
Okrent RA, Wildermuth MC. Evolutionary history of the GH3 family of acyl adenylases in rosids. Plant Mol Biol. 2011;76(6):489–505.
Waterworth WM, Bray CM, West CE. The importance of safeguarding genome integrity in germination and seed longevity. J Exp Bot. 2015;66(12):3549–58.
Calahan D, Dunham M, DeSevo C, Koshland DE. Genetic analysis of desiccation tolerance in Sachharomyces cerevisiae. Genetics. 2011;189(2):507–19.
Waterworth WM, Footitt S, Bray CM, Finch-Savage WE, West CE. DNA damage checkpoint kinase ATM regulates germination and maintains genome stability in seeds. Proc Natl Acad Sci U S A. 2016;113(34):9647–52.
Zhu Y, Wang B, Phillips J, Zhang ZN, Du H, Xu T, Huang LC, Zhang XF, Xu GH, Li WL, et al. Global transcriptome analysis reveals acclimation-primed processes involved in the acquisition of desiccation tolerance in Boea hygrometrica. Plant Cell Physiol. 2015;56(7):1429–41.
Ron D, Walter P. Signal integration in the endoplasmic reticulum unfolded protein response. Nat Rev Mol Cell Biol. 2007;8(7):519.
Vitale A, Boston RS. Endoplasmic reticulum quality control and the unfolded protein response: insights from plants. Traffic. 2008;9(10):1581–8.
Li J, Chu Z, Batoux M, Nekrasov V, Roux M, Chinchilla D, Zipfel C, Jones JD. Specific ER quality control components required for biogenesis of the plant innate immune receptor EFR. Proc Natl Acad Sci U S A. 2009;106(37):15973–8.
Oka OB, Bulleid NJ. Forming disulfides in the endoplasmic reticulum. Biochim Biophys Acta-Mol Cell Res. 2013;1833(11):2425–9.
Sun RZ, Lin CT, Zhang XF, Duan LX, Qi XQ, Gong YH, Deng X. Acclimation-induced metabolic reprogramming contributes to rapid desiccation tolerance acquisition in Boea hygrometrica. Environ Exp Bot. 2018;148:70–84.
Matsumura Y, Sakai J, Skach WR. Endoplasmic reticulum protein quality control is determined by cooperative interactions between Hsp/c70 protein and the CHIP E3 ligase. J Biol Chem. 2013;288(43):31069–79.
Zhou J, Zhang Y, Qi J, Chi Y, Fan B, Yu JQ, Chen Z. E3 ubiquitin ligase CHIP and NBR1-mediated selective autophagy protect additively against proteotoxicity in plant stress responses. PLoS Genet. 2014;10(1):e1004116.
Essers J, Theil AF, Baldeyron C, van Cappellen WA, Houtsmuller AB, Kanaar R, Vermeulen W. Nuclear dynamics of PCNA in DNA replication and repair. Mol Cell Biol. 2005;25(21):9350–9.
Osborne DJ, Boubriak II. DNA and desiccation tolerance. Seed Sci Res. 1994;4(02):175–85.
Farrant JM, Moore JP. Programming desiccation-tolerance: from plants to seeds to resurrection plants. Curr Opin Plant Biol. 2011;14(3):340–5.
Kunz BA, Anderson HJ, Osmond MJ, Vonarx EJ. Components of nucleotide excision repair and DNA damage tolerance in Arabidopsis thaliana. Environ Mol Mutagen. 2005;45(2–3):115–27.
Djilianov D, Genova G, Parvanova D, Zapryanova N, Konstantinova T, Atanassov A. In vitro culture of the resurrection plant Haberlea rhodopensis. Plant Cell Tissue Organ Cult. 2005;80(1):115–8.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Goto MKS. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kent WJ. BLAT--the BLAST-like alignment tool. Genome Res. 2002;12(4):656–64.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(Web Server issue):W316–22.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman W-H, Pagès F, Trajanoski Z, Galon J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
We thank Dr. Tao Xu (Institute of Botany, Chinese Academy of Sciences) for excellent assistance with data analysis.
This work was supported by China-Bulgaria’s Inter-Governmental S&T Cooperation project [grant number 14–10]; and Bulgarian Science Fund [grant number DNTS/China/01/7]. The funding bodies had no role in study design, data collection, analyses and interpretation, decision to publish or preparation of manuscript.
Availability of data and materials
The datasets reported in this paper are available at NCBI GEO (http://www.ncbi.nlm.nih.gov/geo/) with accession number GSE109958. All data generated or analyzed during this study are included in this article (and additional files) or available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of transcriptomic data of de novo assembly of fresh, dehydrated and rehydrated Haberlea rhodopensis. (XLSX 10 kb)
Figure S1. Unigene analysis of Haberlea rhodopensis transcriptomic data. (a) Length distribution of all unigenes. (b) Species similarity distribution. (c) Distribution of the novel DEGs found in this de novo assembly that are unique or commonly shared to different dehydration time points. (d) Distribution of D20 unique genes in KEGG pathways. (TIF 859 kb)
Table S2. Reads and base numbers of raw and clean sequences from ten libraries. (XLSX 10 kb)
Table S3. Expression of genes that were known as dehydration inducible genes (ELIP, GS, beta-amylase, catalase, LEA etc.). (XLSX 22 kb)
Figure S2. Confirmation of RNA-seq data using qRT-PCR and the KEGG analysis of up- and down-regulated transcripts. (a) Validation of RNA-seq data using qPCR. (b) KEGG analysis of up- and down-regulated unigenes of each stage. (TIF 721 kb)
Table S4. Designed primers for the qPCR experiment in this article. (XLSX 11 kb)
Table S5. GO enrichment of up- and down-regulated unigenes of each stage. (XLSX 13 kb)
Table S6. Pathways identified from the functionally organized GO/pathway term network using DEGs in D20 stage. (XLSX 13 kb)
Figure S3. Heatmap of expression patterns of transcripts encoding genes in photosynthesis. Data were log2 transformed and fold change is indicated by color bar. (TIF 1589 kb)
Figure S4. Metabolic network of primary and secondary metabolism constructed by the ClueGO analysis. Mapped genes with changes in D20 are represented by color bar corresponding to the log2 fold change ratio between control and D20. The genes of interest linked to previously published metabolomics data are shown as array on the right. (TIF 734 kb)
About this article
Cite this article
Liu, J., Moyankova, D., Lin, CT. et al. Transcriptome reprogramming during severe dehydration contributes to physiological and metabolic changes in the resurrection plant Haberlea rhodopensis. BMC Plant Biol 18, 351 (2018). https://doi.org/10.1186/s12870-018-1566-0
- Desiccation tolerance
- Haberlea rhodopensis
- Hormone signaling pathway
- Primary and secondary metabolism
- Resurrection plant