Skip to main content

Transcriptome profiling of Bergenia purpurascens under cold stress

Abstract

Bergenia purpurascens is an important medicinal, edible and ornamental plant. It generally grows in high-altitude areas with complex climates. There have been no reports about how B. purpurascens survives under cold stress. Here, the B. purpurascens under low temperature were subjected to transcriptomics analysis to explore the candidate genes and pathways that involved in the cold tolerance of B. purpurascens. Compared with the control treatment, we found 9,600 up-regulated differentially expressed genes (DEGs) and 7,055 down-regulated DEGs. A significant number of DEGs were involved in the Ca2+ signaling pathway, mitogen-activated protein kinase (MAPK) cascade, plant hormone signaling pathway, and lipid metabolism. A total of 400 transcription factors were found to respond to cold stress, most of which belonged to the MYB and AP2/ERF families. Five novel genes were found to be potential candidate genes involved in the cold tolerance of B. purpurascens. The study provide insights into further investigation of the molecular mechanism of how B. purpurascens survives under cold stress.

Peer Review reports

Introduction

Bergenia purpurascens (Hook.f.et Thoms.) Engl., also called Purple Bergenia, is an important medicinal, edible, and ornamental plant belonging to the genus Bergenia in the Saxifragaceae family. It grows extremely slowly and is mainly distributed in highland areas of southwest China and other Asian regions [1]. The whole grass and rhizome can be used as medicine. B. purpurascens is rich in active polyphenolic components like bergenin and arbutin. Bergenin has been successfully developed as a listing drug (Compound Bergenin Tablets) for treating chronic bronchitis [2]. It also treats other diseases such as pneumonia, tuberculosis, gastric ulcer, and acute [3]. At the same time, it has food and ornamental value. B. purpurascens not only has diverse functions but also has strong ecological adaptability. It grows at high altitudes of 2800–4800 m, where the temperature is characterized by a significant difference between day and night and a low average annual temperature [1]. Due to constant exposure to cold stress, B. purpurascens has evolved to tolerate cold stress and is likely to be more winter-hardy than other species of Bergenia [4]. It is a potential ideal material for the study of cold resistance in plants.

Several recent studies have explored the adaptation mechanism of plants grown in high-altitude areas to cold stress. For example, physiological and transcriptome analyses were performed for Ammopiptanthus mongolicus surviving in the Qinghai-Tibetan Plateau after freezing stress (-7.5 °C) [5]. He et al. [6] compared the cold resistance physiological and biochemical features of four Herba Rhodiola seedlings (grown at altitude over 3400 m) under low temperature. Metabolic analyses were conducted on Tibetan hulless barley (Hordeum distichon) exposed to six temperatures (24 °C, 12 °C, 5 °C, 0 °C, -5 °C, -8 °C) [7]. To date, research on B. purpurascens mainly focused on its pharmacological value and chemical components, as well as its morphology and growth characteristics [1, 8]. Few research on molecular biology or its genomic data, not to mention the cold resistance mechanism for B. purpurascens.

Low temperature (cold stress) affects plants’ geographical distribution, leading to reduced yields and threatening plant survival [9]. Cold stress can be classified as chilling stress (0 to 20 °C) and freezing stress (< 0 °C) [10]. Cell membrane structure and fluidity would change under low temperature condition [11]. The cell membrane changes can be sensed by calcium channels and proteins localized on the plasma membrane, leading to the changes in Ca2+ flow and the regulation of low-temperature responsive genes [12]. The CBF (C-repeat binding factors) pathway is one of the most widely reported pathways that responds to low temperature [13]. In addition, plant hormones are essential for plant cold acclimation, including abscisic acid (ABA), gibberellin (GA), ethylene, brassinosteroids (BR) and Jasmonic acid(JA) [14].

Currently, high-throughput sequencing technologies have been widely used in the study of the mechanism of plant response to biotic or abiotic stress [15]. The transcriptomic level of plants generally changes over a short period in response to stress [16]. To explore the mechanisms of B. purpurascens adaptation to cold stress, we performed transcriptome sequencing on B. purpurascens leaves treated at 4 °C for 6 h. Our study provides gene information resources for the analysis of the cold resistance mechanism of B. purpurascens, as well as a reference for the breeding of resistant species of plants.

Materials and methods

Plant growth and cold treatment

The B. purpurascens seedlings for this study were collected on October 13, 2022 in Jiulong County, Sichuan Province, China (~ 101.50°N, ~ 29.00°E), and planted into a mixture of peat soil and gravel at 22 °C under long-day conditions (16 h light of 70 µE m− 2 s− 1 and 8 h dark). The altitude of the sample collection site is 2900 m, where the annual average low temperature was 4 °C and the annual average high temperature was 18 °C in 2021 (meteorological data source: website: https://www.tianqi24.com/jiulong/history2021. html). The species information was confirmed by DNA barcoding technique with ITS2 and psbA-trnH sequences [17]. After being incubated at 4 and 22 °C for 6 h under light of 70 µE m− 2 s− 1, respectively, the first completely expanded leaves from the top were transferred into sample tubes followed by being immediately frozen in liquid nitrogen and finally stored at -80 °C. The 6 h time-point and 4 °C were selected based on the growing conditions of B. purpurascens and previous studies about cold stress on Sonneratia caseolaris, Populus tomentosa, Oxytropis ochrocephala and Arabidopsis [16, 18,19,20]. After these plants were transferred to 4 °C for 6 h, the transcriptome level changed dramatically.

cDNA library preparation and sequencing

The total RNA of the first completely spread leaves from the top was extracted using Trelief® RNAprep Pure Plant Kit (Tsingke, China). RNA concentration and purity were examined using Nanodrop 2000 (Thermo Fisher Scientific, USA), RNA integrity was examined through agarose gel electrophoresis, and RIN values were determined by Agilent 5300 Bioanalyzer (Agilent, USA). After the RNA samples reached the qualified level, the library assembly and RNA-seq analysis were performed on the Illumina Novaseq 6000 high-throughput sequencing platform with the assistance of Shanghai Meiji Biotechnology Company (Shanghai, China). Three biological replicates were performed.

De novo transcriptome assembly

To ensure the quality of the transcriptome sequence results and bioinformatics analysis, the raw sequences (raw reads) obtained from the sequencing were counted and quality controlled by fastp (v.0.19.5, https://github.com/OpenGene/fastp). Then Trinity (v.2.8.5, https://github.com/trinityrnaseq/trinityrnaseq) was used to splice the filtered clean reads to obtain the reference sequences for subsequent bioinformatics analysis. Gene function annotation was performed using six major databases: NR (NCBI non-redundant protein Sequences), Swiss-Prot (A manually annotated and reviewed protein sequence database), GO (Gene Ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes), eggNOG (Evolutionary Genealogy of Genes: Non-supervised Orthologous Groups) and Pfam (Protein family), with E-value ≤ 1e− 5.

Differential expression analysis and functional enrichment

Gene expression levels were calculated by RSEM (v1.3.1, http://deweylab.biostat.wisc.edu/rsem/) [21], and the expression level of the unigenes was converted to TPM (transcripts per million) value. Differential expression analysis was performed by DESeq2 (v.1.24.0, http://bioconductor.org/packages/stats/bioc/DESeq2/). Differentially expressed genes between samples were screened with log2|Fold Change| ≥ 1&Padjust < 0.05. The Python-based library GOATools was used for GO enrichment. The KEGG database was used for enrichment analysis of the signaling pathways of the differentially expressed genes to screen out the functional genes associated with cold stress [22].

Quantitative real-time PCR (qRT-PCR) validation

RNA was extracted from cold-treated and control-treated samples. First-strand cDNA was reverse transcribed from RNA using Evo M-MLV RT mix Kit with gDNA Clean for qRT-PCR Ver.2 (Accurate Biology, China). Ten genes (WRKY33, ZAT10, and eight randomly selected genes) were performed on qRT-PCR to validate the RNA-seq data. 18 S rRNA was used as the internal reference to normalize the expression data [23]. qRT-PCR was performed with SYBR GREEN I (Accurate Biology, China) on the Analytik Jena qTOWER3G fluorescent quantitative PCR instrument (Analytik Jena AG, Germany). The cycling conditions were 95 °C pre-denaturation for 30 s, 95 °C for 5 s, 60 °C for 30 s, 40 cycles. Primers used for qRT-PCR were designed by Primer-Blast and shown in Table S1. Three biological replicates were generated. The 2−ΔΔCt method was used to calculate the relative gene expression. RNA-seq value was based on the fold change of up-regulated or down-regulated DEGs.

Results

Transcriptome sequencing and de novo transcriptome assembly

RNA libraries were constructed and sequenced for a total of six samples from two groups of B. purpurascens samples, including control (CK) and cold (4 °C) treatment. Among them, there were 297,586,224 original fragments. After removing low-quality areas and adapters, 285,734,142 clean reads remained. Q30 (sequences with a sequencing error rate of less than 0.1%) were all above 93%, and the average GC content was 44.97% (Table 1). Clean data from 6 samples were mapped to the assembled unigenes using Trinity. The result showed that the number of unigenes obtained was 109,971.

Table 1 Assembly quality statistics

Functional annotation and DEGs screening

The obtained unigenes were subjected to NR, Swiss-prot, GO, KEGG, eggNOG, and Pfam public databases to obtain annotation information. 48,668 unigenes were annotated in at least one database within a predefined range of E-value. There were 47,573 (43.26%), 30,647 (27.87%), 39,987 (36.36%), 20,134 (18.31%), 33,838 (30.77%) and 24,055 (21.87%) unigenes annotated in NR, Swiss-prot, GO, KEGG, eggNOG, and Pfam, respectively (Table 2). Based on the read counts, DEGs were assessed with the DESeq2 program. According to Padjust < 0.05 and |log2FC| ≥ 1, a total of 16,655 differentially expressed genes (DEGs) were obtained after cold stress treatment, of which 9,600 were up-regulated and 7,055 were down-regulated (Fig. 1, Table S2).

Fig. 1
figure 1

Differentially expressed genes (DEGs) in B. purpurascens in response to cold stress compared to control (CK). (a) Column diagram. (b) Volcano plot. The genes meeting Padjust < 0.05 & |log2FC| ≥ 1 were considered as differentially expressed genes

Besides screening the homologues of the previously reported genes, transcriptomic data also have the potential to discover novel genes. Based on the fold changes of expression, we selected the genes with high expression level (average TPM > 50 in both control and cold-treated groups). Among the ten up-regulated and ten down-regulated genes with the most significant changes in expression, five genes (DN12077_c0_g1, DN12077_c1_g1, DN5468_c1_g1, DN7548_c0_g1, DN95039_c0_g1) were unknown and failed to be annotated (Table S3). These genes were proposed as potential novel candidate genes implicated in the cold tolerance of B. purpurascens.

Table 2 Statistics of functional annotation of non-redundant unigenes sequenced for B. purpurascens against public databases

GO and KEGG enrichment under cold stress

To uncover the function of DEGs, GO enrichment and KEGG enrichment were performed. For GO enrichment, GO terms were classified into three groups: cellular component, biological process, and molecular function. Among the top 20 significantly enriched GO terms (Fig. 2a), 14 were related to biological process, one to cellular component, and five to molecular function. In the biological process, hormone metabolic processes (GO: 0042445), hormone level regulation (GO: 0010817), and hormone-mediated signaling pathways (GO: 0009755) were the most significantly enriched GO terms. KEGG enrichment results showed that 2,093 DEGs under cold stress were enriched into 140 metabolic pathways. As shown in the top 20 items of KEGG enrichment analysis (Fig. 2b), plant-pathogen interaction (Ko04626), plant hormone signal transduction (Ko04075), phenylpropanoid biosynthesis (Ko00940), circadian rhythm-plant (Ko04712), and MAPK signaling pathway (Ko04016) were significantly enriched. In addition, five of the top 20 items were related to lipid metabolism, including alpha- Linolenic acid metabolism (Ko00592), cutin, suberine and wax biosynthesis (Ko00073), glycerolipid metabolism (Ko00561), glycerophospholipid metabolism (Ko00564), and linoleic acid metabolism (Ko00591).

Fig. 2
figure 2

GO and KEGG enrichment analysis of differentially expressed genes (DEGs) in response to cold stress compared to control. (a) GO enrichment of top 20 terms. (b) KEGG enrichment of top 20 pathways

Ca2+ signal transduction

The calcium ion (Ca2+) is an important secondary signaling molecule and plays a vital role in the cold stress response of plants [24]. Figure 3 and Table S4 showed the DEGs identified to be related to Ca2+ signal transduction under cold stress in B. purpurascens. Among them, there were 51 DEGs up-regulated and seven down-regulated. CNGCs are non-selective cation channels that can serve as a pathway for Ca2+ conduction into the cytoplasm [25]. In our study, among 12 CNGC-homologous DEGs, eight were up-regulated (Fig. 3b). Ca2+ sensors include calmodulins (CaMs), CaM-like proteins (CMLs), calcineurin B-like proteins (CBLs), Ca2+-dependent protein kinases (CPKs/CDPKs), and calcineurin B-like-interacting protein kinases (CIPKs) regulated by CBLs [26]. As shown in Fig. 3b, we found 46 DEGs related to Ca2+ sensors, including two CaMs, 17 CMLs, 14 CDPKs, and 13 CIPKs.

Fig. 3
figure 3

Ca2+ signal in response to cold stress. (a) Ca2+ signaling pathway schematic display. (b) Heatmap of DEGs involved in Ca2+ signaling pathway

MAPK signal transduction

In this study, the genes related to MAPK signal transduction led to significant change under cold stress (Table S5). The expression levels of the pathway genes MEKK1 (DN5274_c0_g1), MKK2 (DN10095_c0_g1, DN514_c0_g2), and MAPK4 (DN10834_c0_g2) were significantly increased when the samples were exposed to cold stress (Fig. 4a). In the MAPK signaling pathway activated by ABA, MAPKK17 / 18 (DN15645_c0_g1, DN5158_c0_g1) and MPK7 (DN47710_c0_g1) were up-regulated (Fig. 4b).

Fig. 4
figure 4

MAPK signal in response to cold stress. (a) MEKK1-MKK2-MPK4/6 pathway induced under cold stress. (b) MAPK signaling pathway induced by ABA in response to cold stress

Plant hormone signal transduction

Plant hormones play an important role in the regulation of plant growth and response to abiotic stress [27]. Under cold conditions, environmental stress factors induce the accumulation of plant hormones and trigger the plant hormone signaling pathway [14]. In our study, the results of KEGG enrichment showed a significant amount of DEGs enriched in plant hormone signal transduction (Table S6). Figure 5 showed DEGs in auxin, abscisic acid (ABA), and gibberellin (GA) signal transduction [28,29,30]. Additionally, Figure S3 showed part of the ABA synthesis pathway. In B. purpurascens, multiple homologues of ARF in auxin, ABF in ABA, and PIF3 in GA were down-regulated under cold stress. Moreover, multiple homologues of TIR1, IAA and GH3 in auxin, and PYR/PYL, PP2C and SnRK2 in ABA were up-regulated.

Fig. 5
figure 5

Plant hormone signal in response to cold stress. (a) Auxin, (b) Abscisic acid, (c) Gibberellin

Transcription factors responding to cold stress

In our study, we predicted 1,114 transcription factor genes from the B. purpurascens transcriptome. Among them, 400 transcription factors distributed across 26 gene families were identified as DEGs in response to cold stress, including 79 MYB transcription factors, 50 AP2/ERF transcription factors, 42 C2C2 transcription factors, 32 WRKY transcription factors, 31 bHLH transcription factors, and 27 NAC transcription factors (Table S7, Fig. 6). Based on the result of the hierarchical clustering (Fig. 6a), 400 transcription factors were divided into down-regulation group and up-regulation group with similar numbers. The top three up-regulated families were MYB, AP2/ERF and C2C2, and the top three down-regulated were MYB, C2C2 and bHLH (Fig. 6b). In addition, several well-known cold-related transcription factors, including CBF1 belonging to AP2 gene family and MYB15, were differentially expressed under cold stress as shown in Figure S1.

Fig. 6
figure 6

Differentially expressed transcription factors in B. purpurascens under cold stress. (a) Hierarchical clustering of 400 differentially expressed transcription factors. Expression values (TPM) are log2 - transformed. The colors indicate the gene expression levels from low (blue) to high (red). The detailed information with unigene IDs and annotations is shown in Table S6. (b) The distribution of the top 8 differentially expressed transcription factors in up-regulated and down-regulated categories

Validation of transcriptome data by qRT-PCR analyses

To confirm the DEGs obtained from Illumina Novaseq 6000 platform sequencing were credible, the qRT-PCR was performed on RNA extracted from cold and control-treated samples. Ten DEGs were selected for qRT-PCR, including two transcription factors (WRKY33 and ZAT10) and eight randomly selected genes (Figure S2). Comparison of qRT-PCR and RNA-seq showed that although there were small differences in the expression variability of these genes, the expression trends were consistent. The results suggested that the transcriptome data were reliable.

Discussion

Few studies have reported the genome and transcriptome of B. purpurascens, which hinder the understanding of the molecular mechanisms of its adaptation to severe highland environments and tolerance to abiotic stress. To our knowledge, this is the first transcriptomic analysis of B. purpurascens in response to cold stress. The discovery of cold-response genes and pathways in the study provided valuable insights into underlying the molecular mechanisms of cold tolerance of B. purpurascens. The study showed that a significant number of cold-response genes were involved in diverse biological processes, primarily encompassing lipid metabolism, Ca2+ signaling, MAPK, plant hormones, and transcription factors, which require further discussion.

DEGs involved in lipid metabolism

When the temperature decreases, the cellular membrane functions as the primary location for sensing temperature signals [9]. Lipids are one of the major components of biological membranes [31]. In the present study, KEGG enrichment suggested that numerous DEGs were significantly enriched in lipid metabolism pathways (Fig. 2b). These pathways might be related to the changes in membrane lipid composition in B. purpurascens under cold stress, which was similar to the membrane lipid change-related pathways in blueberry in response to cold stress [32]. Szymanski et al. [33] elucidated the association between the expression of lipid metabolism genes and membrane lipid composition. Studies conducted on potato [34] and pineapple [35] also showed that temperature alterations could induce changes in lipid components in the cell membrane to reduce cold damage. Meanwhile, the up-regulation of seven homologues of Fatty Acid Desaturase (FAD) genes (Table S8) suggested the enhanced production of unsaturated fatty acids and the consequent elevation of membrane fluidity in B. purpurascens. In previous studies, over-expression of Arabidopsis FAD7 in tobacco exhibited lesser chilling injuries compared to the wild-type [36]. Similarly, the peach fruit that was exposed to cold stress up-regulated FAD genes, delayed the degradation of phospholipids, and altered the degree of fatty acid unsaturation to mitigate the cold damage [37].

Ca2+ signal under cold stress

We noted that Ca2+ signaling related genes exhibit significant expression changes under cold stress, suggesting that B. purpurascens might stimulate Ca2+ signaling to transfer extracellular signals to the intracellular to improve its adaptive capacity to cope with the cold environment. To our data, 8 CNGC homologues in B. purpurascens were induced in response to cold stimulus, showing that low temperature may activate CNGC channels, allowing Ca2+ to enter the cell. Similarly, analysis of CNGCs expression patterns in Brassica oleracea revealed that half of the CNGCs were up-regulated after incubation at 4 °C [38]. CNGC9 was shown to regulate chilling tolerance by mediating cytoplasmic calcium elevation in rice [39]. In Arabidopsis, CNGC4 was also reported to be associated with low-temperature tolerance [40]. Ca2+ sensors (CaMs, CMLs, CDPKs, CBLs, and CIPKs) play a role in sensing the Ca2+ signaling and would be induced or inhibited by abiotic stress in various plants [41]. CaMs/CMLs and CBLs are considered to be non-catalytic relay sensors, while CDPKs are acknowledged as sensory responders with catalytic activity [24]. Our data indicated that a large number of genes encoding CAMs/CMLs, CDPKs, and CIPKs significantly responded to cold stress (Fig. 3, Table S4). Among these sensor genes, three CMLs homologues (DN42086_c0_g1, DN9345_c0_g2, and DN9345_c0_g1) were significantly increased and highly expressed in B. purpurascens under cold stress, which can be considered for further validation of their functions and exploration of downstream genes. The Arabidopsis CML gene family comprises 50 genes, several of which were found to respond to diverse biotic and abiotic stressors. [42]. In Medicago truncatula, CML5 positively regulated the expression of CBF42 and CBF1 to increase cold tolerance [43]. CML10 was identified to interact with cytosolic enzymes GSTU8 and FBA6 to regulate cold tolerance in Medicago sativa [44].

MAPK signal under cold stress

Through the analysis of transcriptome data, genes related to MAPK cascades (Fig. 4, Table S5) were induced by cold stress, which might enhance the cold tolerance of B. purpurascens. The MEKK1-MKK2-MPK4/MPK6 pathway is one of the well-studied MAPK cascades functioning under cold stress [45]. In Arabidopsis, it has been shown that cold temperature induces phosphorylation of MKK2 protein, MKK2 stimulates phosphorylation of MPK4/MPK6 protein, and MPK4/MPK6 regulates cold tolerance in plants by regulating ICE1 and CBF gene expression [46]. CRLK1, a calcium/calmodulin-regulated receptor-like kinase, has been shown to interact with MEKK1 to link calcium signaling and MAPK cascades [47]. In contrast to earlier findings, the CRLK1 homologue (DN14329_c0_g1) in B. purpurascens was down-regulated under cold stress (Table S4), which might be caused by negative feedback regulation in B. purpurascens or other unidentified factors. It is worth noting that in addition to the MEKK1-MKK2-MPK4/MPK6 pathway, the MAPK module consisting of MAP3K17/18-MKK3- MPK1/2/7/14 may also play a role in cold tolerance in B. purpurascens. This pathway could interact directly or indirectly with ABA signaling [48]. It has been shown that ABA induced the activation of a MAPK-like enzyme in epidermal tissue from peas [49]. Danquah et al. [50] revealed that MAP3K17 and MAP3K18 were regulated by the ABA core signaling module in Arabidopsis.

Plant hormones signal under cold stress

Hormones generally participate in regulating plant response to cold stress [51]. Under cold stress, plants increase ABA accumulation and activate the ABA signal pathway to tolerate stress [52]. 9-cis-epoxy carotenoid dioxygenase (NCED) is one of the rate-limiting enzymes in ABA biosynthesis [53]. It has been shown that the expression level of NCED was up-regulated after cold stress in some plants like blueberries [32], Pyrus ussuriensis [54], and Moso bamboo [55]. Similarly, we noticed that 3 NCED homologues (DN41962_c0_g1, DN9413_c0_g1, DN9413 c0_g2) in B. purpurascens were up-regulated under cold stress (Figure S3), suggesting that low temperature might increase endogenous ABA levels. PYL and PP2C are key genes in the ABA signaling pathway, and water dissipation will decrease when their expression levels increase [56]. In the present study, PYLs and PP2Cs were up-regulated in response to cold stress. Auxin also plays an important role under cold stress in plants such as Arabidopsis [57], rice [58], and Brassica napus [59]. The Aux/IAA proteins are a large family of auxin co-receptors involved in various physiological and developmental processes [60]. In Arabidopsis, the expression of IAA5 and IAA19 is directly promoted by CBF1 and DREB2A TFs in response to abiotic stress [61]. Similarly, five IAA homologues were up-regulated in B. purpurascens under cold stress, but whether they interact with CBF requires further study. The expression patterns of ABA and auxin signaling pathway genes (Fig. 5, Table S6) in B. purpurascens under cold stress were similar to those in Sonneratia caseolaris [18] and Malus sieversii [62] under cold stress. DELLA activity is essential in GA signaling [30]. Achard et al. [63] observed that the accumulation of DELLA proteins induced by CBF1 played a synergistic role in cold acclimation in conjunction with the activation of COR pathways induced by CBF1. In this study, one highly expressed DELLA homologue (DN628_c0_g1) was induced, and another DELLA (DN12432_c0_g1) with a low expression level was reduced. It indicated that the GA signaling pathway might also be involved in cold tolerance of B. purpurascens.

Transcription factors responding to cold stress

Numerous transcription factors from diverse families have been reported to be involved in plant cold resistance [64, 65]. To our data, most of the cold-induced transcription factors in B. purpurascens were consistent with maize [66] and asparagus bean [67]. Improving the regulatory capacity of a key transcription factor may alter the expression of multiple downstream functional genes and enhance plant cold resistance, like ERF105 in Arabidopsis [68] and WRKY63 in rice [69]. Among the transcription factors of B. purpurascens, ERF026 homologue (DN8708_c0_g1) with high expression was the most responsive to cold stress. However, whether it plays a critical role in cold resistance needs further studies to prove. The AP2/ERF family includes numerous important transcription factors with AP2 DNA-binding domains. Among them, the CBF factor is a key regulator in response to cold stress. In Arabidopsis, CBF1/DREB1B, CBF2/DREB1C, and CBF3/DREB1A are rapidly induced by low temperature to enhance freezing resistance [70]. In addition to model plants, it has been reported that CBF plays a positive role in cold tolerance in birch [71], soybean [72], and barley [73]. In Arabidopsis, ZAT10 and WRKY33 act as positive regulators, whereas MYB15 and EIN3 act as negative regulators of the CBF genes. To our data, we identified one CBF1 homologue, one ZAT10, one MYB15, and two WRKY33 in B. purpurascens (Figure S1), and their expression variations were consistent with previous studies. Additionally, the consistency of qRT-PCR and RNA-seq further indicated that WRKY33 and ZAT10 could be candidate genes for the following cold-tolerance functional studies (Figure S2).

Conclusion

The transcriptomic analysis of B. purpurascens under low temperature revealed that 46 Ca2+ sensors related genes (two CaMs, 17 CMLs, 14 CDPKs, and 13 CIPKs) and 400 transcription factor genes (79 MYBs, 50 AP2/ERFs, 42 C2C2s, 32 WRKYs, 31 bHLHs, 27 NACs, and others) were responding to cold stress. Two MAPK cascades, five lipid metabolism pathways, and three kinds of plant hormone signal pathways (ABA, auxin, and GA) may play vital roles in enhancing the cold resistance of B. purpurascens. Five novel genes were suggested to be potential candidate genes involved in the cold tolerance of B. purpurascens. Additionally, the enormous amount of transcriptome data generated here will fill the gap in B. purpurascens transcriptome data.

Data Availability

The sequence data have been deposited in the NCBI under BioProject accession number PRJNA981985 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA981985).

References

  1. Koul B, Kumar A, Yadav D, Jin JO. Bergenia Genus: traditional uses, Phytochemistry and Pharmacology. Molecules. 2020;25:5555.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Zhang SS, Liao ZX, Huang RZ, Gong CC, Ji LJ, Sun H. A new aromatic glycoside and its anti-proliferative activities from the leaves of Bergenia purpurascens. Nat Prod Res. 2018;32:668–75.

    Article  CAS  PubMed  Google Scholar 

  3. Chen XM, Yoshida T, Hatano T, Fukushima M, Okuda T. Galloylarbutin and other polyphenols from Bergenia purpurascens. Phytochemistry. 1987;26:515–7.

    Article  CAS  Google Scholar 

  4. Salimo ZM, Yakubu MN, da Silva EL, de Almeida ACG, Chaves YO, Costa EV, et al. Chemistry and Pharmacology of Bergenin or its derivatives: a promising molecule. Biomolecules. 2023;13:403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yang Z, Liu Y, Han H, Zhao X, Chen S, Li G, et al. Physiological and transcriptome analyses reveal the response of Ammopiptanthus mongolicus to extreme seasonal temperatures in a cold plateau desert ecosystem. Sci Rep. 2022;12:10630.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. He S, Zhao K, Ma L, Yang J, Chang Y, Ashraf MA. Comparison of cold resistance physiological and biochemical features of four Herba Rhodiola seedlings under low temperature. Saudi J Biol Sci. 2016;23:198–204.

    Article  CAS  PubMed  Google Scholar 

  7. Yang C, Yang H, Xu Q, Wang Y, Sang Z, Yuan H. Comparative metabolomics analysis of the response to cold stress of resistant and susceptible tibetan hulless barley (Hordeum distichon). Phytochemistry. 2020;174:112346.

    Article  CAS  PubMed  Google Scholar 

  8. Qu Y, Zhang C, Liu R, Wu H, Sun Y, Zhang N, et al. Rapid characterization the chemical constituents of Bergenia purpurascens and explore potential mechanism in treating osteoarthritis by Ultra high performance liquid chromatography coupled with quadrupole time-of‐flight mass spectrometry combined with network pharmacology. J Sep Sci. 2020;43:3333–48.

    Article  CAS  PubMed  Google Scholar 

  9. Ding Y, Shi Y, Yang S. Advances and challenges in uncovering cold tolerance regulatory mechanisms in plants. New Phytol. 2019;222:1690–704.

    Article  PubMed  Google Scholar 

  10. Guo X, Liu D, Chong K. Cold signaling in plants: insights into mechanisms and regulation. J Integr Plant Biol. 2018;60:745–56.

    Article  PubMed  Google Scholar 

  11. Rawat N, Singla-Pareek SL, Pareek A. Membrane dynamics during individual and combined abiotic stresses in plants and tools to study the same. Physiol Plant. 2021;171:653–76.

    Article  CAS  PubMed  Google Scholar 

  12. Ding Y, Shi Y, Yang S. Molecular regulation of plant responses to environmental temperatures. Mol Plant. 2020;13:544–64.

    Article  CAS  PubMed  Google Scholar 

  13. Park S, Lee CM, Doherty CJ, Gilmour SJ, Kim Y, Thomashow MF. Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J. 2015;82:193–207.

    Article  CAS  PubMed  Google Scholar 

  14. Manasa SL, Panigrahy M, Panigrahi KCS, Rout GR. Overview of cold stress regulation in plants. Bot Rev. 2022;88:359–87.

    Article  Google Scholar 

  15. Urano K, Kurihara Y, Seki M, Shinozaki K. Omics’ analyses of regulatory networks in plant abiotic stress responses. Curr Opin Plant Biol. 2010;13:132–8.

    Article  CAS  PubMed  Google Scholar 

  16. Usadel B, Bläsing O, Gibon Y, Poree F, Hohne M, Gunter M, et al. Multilevel genomic analysis of the response of transcripts, enzyme activities and metabolites in Arabidopsis rosettes to a Progressive decrease of temperature in the non-freezing range. Plant Cell Environ. 2008;31:518–47.

    Article  CAS  PubMed  Google Scholar 

  17. Wei Y, Zhang Y, Li Q, Shen K, Zhou J. Identification of DNA barcode of medicinal plants containing bergenin. West China J Pharm Sci. 2021;36:44–50.

    Google Scholar 

  18. Yang Y, Zheng C, Zhong C, Lu T, Gul J, Jin X, et al. Transcriptome analysis of Sonneratia Caseolaris seedlings under chilling stress. PeerJ. 2021;9:e11506.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Yang X, Zhao T, Rao P, Gao K, Yang X, Chen Z, et al. Transcriptome profiling of Populus tomentosa under cold stress. Ind Crops Prod. 2019;135:283–93.

    Article  CAS  Google Scholar 

  20. He W, Zhuang H, Fu Y, Guo L, Guo B, Guo L, et al. De novo Transcriptome Assembly of a Chinese locoweed (Oxytropis Ochrocephala) Species provides insights into genes Associated with Drought, Salinity, and Cold Tolerance. Front Plant Sci. 2015;6:1086.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ogata H, Goto S. KEGG: Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  Google Scholar 

  23. Yin D, Zhao Z, Guo F, Zhou G, Wang S, Li J. Screening of reference genes for real-time quantitative PCR in Bergenia purpurascens. Genom Appl Biol. 2017;36:4256–62.

    Google Scholar 

  24. Dodd AN, Kudla J, Sanders D. The language of calcium signaling. Annu Rev Plant Biol. 2010;61:593–620.

    Article  CAS  PubMed  Google Scholar 

  25. Duszyn M, Świeżawska B, Szmidt-Jaworska A, Jaworski K. Cyclic nucleotide gated channels (CNGCs) in plant signalling-current knowledge and perspectives. J Plant Physiol. 2019;241:153035.

    Article  CAS  PubMed  Google Scholar 

  26. DeFalco TA, Bender KW, Snedden WA. Breaking the code: Ca2+ sensors in plant signalling. Biochem J. 2009;425:27–40.

    Article  PubMed  Google Scholar 

  27. Waadt R, Seller CA, Hsu PK, Takahashi Y, Munemasa S, Schroeder JI. Plant hormone regulation of abiotic stress responses. Nat Rev Mol Cell Biol. 2022;23:680–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Leyser O. Auxin Signaling. Plant Physiol. 2018;176:465–79.

    Article  CAS  PubMed  Google Scholar 

  29. Chen K, Li GJ, Bressan RA, Song CP, Zhu JK, Zhao Y. Abscisic acid dynamics, signaling, and functions in plants. J Integr Plant Biol. 2020;62:25–54.

    Article  CAS  PubMed  Google Scholar 

  30. Gomi K, Matsuoka M. Gibberellin signalling pathway. Curr Opin Plant Biol. 2003;6:489–93.

    Article  CAS  PubMed  Google Scholar 

  31. Hou Q, Ufer G, Bartels D. Lipid signalling in plant responses to abiotic stress. Plant Cell Environ. 2016;39:1029–48.

    Article  CAS  PubMed  Google Scholar 

  32. Zhang F, Ji S, Wei B, Cheng S, Wang Y, Hao J, et al. Transcriptome analysis of postharvest blueberries (Vaccinium corymbosum ‘Duke’) in response to cold stress. BMC Plant Biol. 2020;20:80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Szymanski J, Brotman Y, Willmitzer L, Cuadros-Inostroza Á. Linking Gene expression and membrane lipid composition of Arabidopsis. Plant Cell. 2014;26:915–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. De Palma M, Grillo S, Massarelli I, Costa A, Balogh G, Vigh L, et al. Regulation of desaturase gene expression, changes in membrane lipid composition and freezing tolerance in potato plants. Mol Breed. 2007;21:15–26.

    Article  Google Scholar 

  35. Zhou Y, Pan X, Qu H, Underhill SJ. Low temperature alters plasma membrane lipid composition and ATPase activity of pineapple fruit during blackheart development. J Bioenerg Biomembr. 2014;46:59–69.

    Article  CAS  PubMed  Google Scholar 

  36. Kodama H, Horiguchi G, Nishiuchi T, Nishimura M, Iba K. Fatty acid desaturation during Chilling Acclimation is one of the factors involved in conferring low-temperature tolerance to young Tobacco leaves. Plant Physiol. 1995;107:1177–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Song C, Wang K, Xiao X, Liu Q, Yang M, Li X et al. Membrane lipid metabolism influences chilling injury during cold storage of peach fruit. Food Res Int. 2022;157.

  38. Kakar KU, Nawaz Z, Kakar K, Ali E, Almoneafy AA, Ullah R, et al. Comprehensive genomic analysis of the CNGC gene family in Brassica oleracea: novel insights into synteny, structures, and transcript profiles. BMC Genomics. 2017;18:869.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Wang J, Ren Y, Liu X, Luo S, Zhang X, Liu X, et al. Transcriptional activation and phosphorylation of OsCNGC9 confer enhanced chilling tolerance in rice. Mol Plant. 2021;14:315–29.

    Article  CAS  PubMed  Google Scholar 

  40. Thoen M, Davila Olivas N, Kloth K, Coolen S, Huang P, Aarts M, et al. Genetic architecture of plant stress resistance: multi-trait genome‐wide association mapping. New Phytol. 2016;213:1346–62.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Iqbal Z, Memon AG, Ahmad A, Iqbal MS. Calcium mediated cold acclimation in plants: Underlying Signaling and Molecular mechanisms. Front Plant Sci. 2022;13:855559.

    Article  PubMed  PubMed Central  Google Scholar 

  42. McCormack E, Tsai YC, Braam J. Handling calcium signaling: Arabidopsis CaMs and CMLs. Trends Plant Sci. 2005;10:383–9.

    Article  CAS  PubMed  Google Scholar 

  43. Sun Q, Huang R, Zhu H, Sun Y, Guo Z. A novel Medicago truncatula calmodulin-like protein (MtCML42) regulates cold tolerance and flowering time. Plant J. 2021;108:1069–82.

    Article  CAS  PubMed  Google Scholar 

  44. Yu S, Wu J, Sun Y, Zhu H, Sun Q, Zhao P, et al. A calmodulin-like protein (CML10) interacts with cytosolic enzymes GSTU8 and FBA6 to regulate cold tolerance. Plant Physiol. 2022;190:1321–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhao C, Wang P, Si T, Hsu C, Wang L, Zayed O et al. MAP kinase cascades regulate the Cold response by modulating ICE1 protein Stability. Dev Cell. 2017;43:618 – 29.e5.

  46. Teige M, Scheikl E, Eulgem T, Dóczi R, Ichimura K, Shinozaki K, et al. The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell. 2004;15:141–52.

    Article  CAS  PubMed  Google Scholar 

  47. Yang T, Shad Ali G, Yang L, Du L, Reddy AS, Poovaiah BW. Calcium/calmodulin-regulated receptor-like kinase CRLK1 interacts with MEKK1 in plants. Plant Signal Behav. 2010;5:991–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Danquah A, de Zelicourt A, Colcombet J, Hirt H. The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol Adv. 2014;32:40–52.

    Article  CAS  PubMed  Google Scholar 

  49. Burnett EC, Desikan R, Moser RC, Neill SJ. ABA activation of an MBP kinase in Pisum sativum epidermal peels correlates with stomatal responses to ABA. J Exp Bot. 2000;51:197–205.

    Article  CAS  PubMed  Google Scholar 

  50. Danquah A, de Zélicourt A, Boudsocq M, Neubauer J, Frei dit Frey N, Leonhardt N, et al. Identification and characterization of an ABA-activated MAP kinase cascade in Arabidopsis thaliana. Plant J. 2015;82:232–44.

    Article  CAS  PubMed  Google Scholar 

  51. Eremina M, Rozhon W, Poppenberger B. Hormonal control of cold stress responses in plants. Cell Mol Life Sci. 2016;73:797–810.

    Article  CAS  PubMed  Google Scholar 

  52. Verma V, Ravindran P, Kumar PP. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016;16:86.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Dong T, Park Y, Hwang I. Abscisic acid: biosynthesis, inactivation, homoeostasis and signalling. Essays Biochem. 2015;58:29–48.

    Article  PubMed  Google Scholar 

  54. Yang T, Huang XS. Deep sequencing-based characterization of transcriptome of Pyrus ussuriensis in response to cold stress. Gene. 2018;661:109–18.

    Article  CAS  PubMed  Google Scholar 

  55. Liu Y, Wu C, Hu X, Gao H, Wang Y, Luo H, et al. Transcriptome profiling reveals the crucial biological pathways involved in cold response in Moso bamboo (Phyllostachys edulis). Tree Physiol. 2020;40:538–56.

    Article  CAS  PubMed  Google Scholar 

  56. Hsu PK, Dubeaux G, Takahashi Y, Schroeder JI. Signaling mechanisms in abscisic acid-mediated stomatal closure. Plant J. 2020;105:307–21.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Rahman A. Auxin: a regulator of cold stress response. Physiol Plant. 2012;147:28–35.

    Article  PubMed  Google Scholar 

  58. Jain M, Khurana JP. Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. FEBS J. 2009;276:3148–62.

    Article  CAS  PubMed  Google Scholar 

  59. Gavelienė V, Novickienė L, Pakalniškytė L. Effect of auxin physiological analogues on rapeseed (Brassica napus) cold hardening, seed yield and quality. J Plant Res. 2013;126:283–92.

    Article  PubMed  Google Scholar 

  60. Luo J, Zhou JJ, Zhang JZ. Aux/IAA Gene Family in plants: molecular structure, regulation, and function. Int J Mol Sci. 2018;19:259.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Shani E, Salehin M, Zhang Y, Sanchez SE, Doherty C, Wang R, et al. Plant stress tolerance requires auxin-sensitive Aux/IAA Transcriptional Repressors. Curr Biol. 2017;27:437–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Zhou P, Li X, Liu X, Wen X, Zhang Y, Zhang D. Transcriptome profiling of Malus sieversii under freezing stress after being cold-acclimated. BMC Genomics. 2021;22.

  63. Achard P, Gong F, Cheminant S, Alioua M, Hedden P, Genschik P. The Cold-Inducible CBF1 factor–dependent signaling pathway modulates the Accumulation of the growth-repressing DELLA proteins via its Effect on Gibberellin Metabolism. Plant Cell. 2008;20:2117–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Baillo EH, Kimotho RN, Zhang Z, Xu P. Transcription Factors Associated with Abiotic and biotic stress tolerance and their potential for crops improvement. Genes (Basel). 2019;10:771.

    Article  CAS  PubMed  Google Scholar 

  65. Kidokoro S, Shinozaki K, Yamaguchi-Shinozaki K. Transcriptional regulatory network of plant cold-stress responses. Trends Plant Sci. 2022;27:922–35.

    Article  CAS  PubMed  Google Scholar 

  66. Waititu JK, Cai Q, Sun Y, Sun Y, Li C, Zhang C et al. Transcriptome profiling of Maize (Zea mays L.) leaves reveals Key Cold-responsive genes, transcription factors, and metabolic pathways regulating cold stress tolerance at the Seedling Stage. Genes. 2021;12.

  67. Miao M, Tan H, Liang L, Huang H, Chang W, Zhang J, et al. Comparative transcriptome analysis of cold-tolerant and -sensitive asparagus bean under chilling stress and recovery. PeerJ. 2022;10:e13167.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Bolt S, Zuther E, Zintl S, Hincha DK, Schmülling T. ERF105 is a transcription factor gene of Arabidopsis thaliana required for freezing tolerance and cold acclimation. Plant Cell Environ. 2016;40:108–20.

    Article  PubMed  Google Scholar 

  69. Zhang M, Zhao R, Huang K, Huang S, Wang H, Wei Z, et al. The OsWRKY63–OsWRKY76–OsDREB1B module regulates chilling tolerance in rice. Plant J. 2022;112:383–98.

    Article  CAS  PubMed  Google Scholar 

  70. Novillo F, Medina J, Salinas J. Arabidopsis CBF1 and CBF3 have a different function than CBF2 in cold acclimation and define different gene classes in the CBF regulon. Proc Natl Acad Sci U S A. 2007;104:21002–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Welling A, Palva ET. Involvement of CBF Transcription Factors in Winter Hardiness in Birch. Plant Physiol. 2008;147:1199–211.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Kidokoro S, Watanabe K, Ohori T, Moriwaki T, Maruyama K, Mizoi J, et al. Soybean DREB1/CBF-type transcription factors function in heat and drought as well as cold stress‐responsive gene expression. Plant J. 2015;81:505–18.

    Article  CAS  PubMed  Google Scholar 

  73. Marozsán-Tóth Z, Vashegyi I, Galiba G, Tóth B. The cold response of CBF genes in barley is regulated by distinct signaling mechanisms. J Plant Physiol. 2015;181:42–9.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was supported by grants from the National Natural Science Foundation of China (No.31900164), the Provincial Natural Science Foundation of Sichuan (No.2022NSFSC1612) and the Fundamental Research Funds for the Central Universities (No. 2682023ZTPY073).

Author information

Authors and Affiliations

Authors

Contributions

Zhu QK, Zhang XB, and Yu F designed the project and managed the project; Zhang XB and Yu F ran the experiments. Zhang XB, Yu F, Lyu X, Chen JY, Xu NM, Zeng HY and Wu YF analyzed the results. Zhang XB and Zhu QK wrote the manuscript.

Corresponding author

Correspondence to Qiankun Zhu.

Ethics declarations

Ethics approval and consent to participate

The Bergenia purpurascens used in this study were collected in their wild populations in non-protected areas. All plant materials used in this experiment were not endangered or protected species. The collection of plant material complies with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

12864_2023_9850_MOESM1_ESM.doc

Supplementary Material 1: Figure S1-S3. Figure S1. The expression profiles of differentially expressed transcription factors screened related to the CBF pathway. Figure S2. qRT-PCR performed on WRKY33, ZAT10 and 8 randomly selected genes. Figure S3. Expression levels of genes related to ABA synthesis in B. purpurascens after cold stress treatment.

12864_2023_9850_MOESM2_ESM.doc

Supplementary Material 2: Figure S4. Standard curves, amplification curves and melting curves of 18S rRNA and 10 selected DEGs for qRT-PCR.

12864_2023_9850_MOESM3_ESM.xlsx

Supplementary Material 3: Table S1-S8. Table S1. Sequences of primers validated by qRT-PCR in this study. Table S2. Identification of differentially expressed genes. Table S3. Ten up-regulated and ten down-regulated unigenes with the most significant changes in expression. Table S4. DEGs related to Ca2+ signal. Table S5. DEGs related to MAPK signaling. Table S6. DEGs related to plant hormones signal. Table S7. Differentially expressed transcription factors under cold stress. Table S8. FAD (Fatty acid desaturase) related genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, X., Yu, F., Lyu, X. et al. Transcriptome profiling of Bergenia purpurascens under cold stress. BMC Genomics 24, 754 (2023). https://doi.org/10.1186/s12864-023-09850-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09850-z

Keywords