Transcriptome analysis of chrysanthemum (Dendranthema grandiflorum) in response to low temperature stress

Background Chrysanthemum is one kind of ornamental plant well-known and widely used in the world. However, its quality and production were severely affected by low temperature conditions in winter and early spring periods. Therefore, we used the RNA-Seq platform to perform a de novo transcriptome assembly to analyze chrysanthemum (Dendranthema grandiflorum) transcription response to low temperature. Results Using Illumina sequencing technology, a total of 86,444,237 high-quality clean reads and 93,837 unigenes were generated from four libraries: T01, controls; T02, 4 °C cold acclimation (CA) for 24 h; T03, − 4 °C freezing treatments for 4 h with prior CA; and T04, − 4 °C freezing treatments for 4 h without prior CA. In total, 7583 differentially expressed genes (DEGs) of 36,462 annotated unigenes were identified. We performed GO and KEGG pathway enrichment analyses, and excavated a group of important cold-responsive genes related to low temperature sensing and signal transduction, membrane lipid stability, reactive oxygen species (ROS) scavenging and osmoregulation. These genes encode many key proteins in plant biological processes, such as protein kinases, transcription factors, fatty acid desaturase, lipid-transfer proteins, antifreeze proteins, antioxidase and soluble sugars synthetases. We also verified expression levels of 10 DEGs using quantitative real-time polymerase chain reaction (qRT-PCR). In addition, we performed the determination of physiological indicators of chrysanthemum treated at low temperature, and the results were basically consistent with molecular sequencing results. Conclusion In summary, our study presents a genome-wide transcript profile of Dendranthema grandiflorum var. jinba and provides insights into the molecular mechanisms of D. grandiflorum in response to low temperature. These data contributes to our deeper relevant researches on cold tolerance and further exploring new candidate genes for chilling-tolerance and freezing-tolerance chrysanthemum molecular breeding. Electronic supplementary material The online version of this article (10.1186/s12864-018-4706-x) contains supplementary material, which is available to authorized users.


Background
Low temperature is a common environmental stress factor, which seriously affects plant growth and becomes an important factors limiting plant distribution. Plants will make a series of physiological and biochemical changes as defense systems in response to chilling (0-15°C) and freezing (< 0°C) injury [1]. These cold resistance mechanisms involve biofilm system, osmotic regulators, protective enzyme systems, and fatty acids. For many temperate plant species, exposed to chilling temperature for a period of time could increase their tolerance to subsequent freezing conditions, and this process is called 'cold acclimation (CA)' [2]. Cold acclimation can regulate many cold-responsive genes differentially expressed, induce the establishment of cold resistance mechanism, and finally improve the adaptability of plants to low temperature [3].
Dendranthema grandiflorum var. jinba is one of the best-selling cut chrysanthemum varieties in China. In order to realize the annual production and supply of cut chrysanthemum, in most areas of China, factories need to use the facilities for heating in winter, which greatly increases the cost of production. Thus, understanding the mechanism of chilling and freezing stress responses and improving the cold tolerance of chrysanthemum by gene transfer are of great importance.
Chrysanthemum has large genomes and lacks genomic information on the basis of cold tolerance. With the development of next-generation sequencing (NGS) technology, large scale transcriptome data has become availavle in various species [4]. In recent years, Illumina RNA-Seq technology has been successfully applied to many plant species, such as Beta vulgaris [5], Spartina pectinata [6], Camellia sinensis [7], Lilium lancifolium [8], and Camellia japonica [9], for its high accuracy and sensitivity of gene discovery.
In our study, based on Illumina NGS technology, a fully characterized chrysanthemum transcriptome was represented. Physiological experiments and molecular sequencing analysis combined to explore the cold mechanism of chrysanthemum under chilling and freezing stresses. Moreover, we compare and analyze the physiological and molecular aspects of chrysanthemum with and without cold acclimation to explore the effect of acclimation on chrysanthemum. Lots of cold-induced genes were identified, and some of them were of great importance in coldtolerance chrysanthemum molecular breeding.

Plant materials and low temperature treatments
Dendranthema grandiflorum var. jinba was used in this study. The buds raised from tissue-cultured seedlings were grown on MS medium (16 h photoperiod, 25°C/22°C day/ night temperature) for twenty days. Then twenty-day old chrysanthemum seedlings were transferred to pots filled with a 1:1 mixture of peat and perlite, and acclimated for three days at normal condition. Then the four groups of seedlings were respectively subjected to following treatments: (T01) normal condition as control, (T02) 4°C for 24 h, (T03) 4°C for 24 h, followed by − 4°C for 4 h, (T04) -4°C for 4 h. Each treatment collected a mixture of three biologically replicated leaves, then samples rapidly frozen with liquid nitrogen and stored at − 80°C. There were four samples in total used for RNA-Seq and differential expression analyses.

RNA preparation
Total RNA was extracted according to manufacturer's instructions. 1% agarose gels was used for the monitoring of RNA degradation and contamination; the Nano-Photometer® spectrophotometer (IMPLEN, CA, USA) was used for check of RNA purity; Qubit® RNA Assay Kit was used for measurement of RNA concentration; and the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA) was used for assessing of RNA integrity.

Library preparation for transcriptome sequencing
A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext®Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) according to manufacturer's recommendations. Briefly, mRNA was purified from total RNA by magnetic beads and cut randomly into short fragments by fragmentation buffer. First strand cDNA was synthesized by random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-), using mRNA as a template. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. The cDNA fragments of 150~200 bp in length were selected by AMPure XP beads, after end-repair and single nucleotide A (adenine) addition. Then, the cDNA library was obtained by PCR enrichment. At last, PCR products were purified and library quality was assessed on the Agilent Bioanalyzer 2100 system. The library preparations were sequenced on an Illumina Hiseq 4000 platform and paired-end reads were generated.

Transcriptome assembly and gene functional annotation
The clean data were obtained by removing reads containing adapter and low quality reads from raw data. Transcriptome assembly was performed using Trinity [10] based on clean data with high quality. The sequences of unigenes were compared with Nr (NCBI non-redundant), Pfam (Protein family), Swiss-Prot (a manually annotated and reviewed protein sequence database), GO (Gene Ontology), COG (Cluster of orthologous groups), KOG (euKaryotic Orthologous Groups) and eggNOG (evolutionary genealogy of genes: Nonsupervised Orthologous Groups) databases by BLAST, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthology results were obtained by comparing with KEGG using KOBAS2.0 [11]. After predicting the amino acid sequence of unigenes, we use HMMER [12] to compare them with Pfam database to get the annotation information of unigenes.

Differential expression analysis
EdgeR program package was used for adjusting read counts of each sequenced library, and DEGseq (2010) R package was used for differential expression analysis of two samples. P value was adjusted using q value [13]. In this study, q value < 0.01 and |log2 (fold change)| > 1 were set as the threshold for significantly differential expression.

Validation of RNA-Seq data by qRT-PCR
Expression of Ten DEGs in chrysanthemum was analyzed by qRT-PCR with three replicates to verify the RNA-Seq data. The qRT-PCR was performed using the SsoFast EvaGreen supermix (Bio-Rad, Hercules, CA, USA) and Bio-Rad CFX96™ detection system. A final 10 μL qRT-PCR reaction mixture contained 5 uL SsoFast EvaGreen supermix, 1 uL diluted cDNA sample, and 150 nM primers. Relative expression levels were calculated by the 2 −ΔΔCT method, and the Elongation Factor 1α (EF1α) gene (Forward primer: TTTTGGTATCTGGTCCTGGAG, Reverse primer: CCATTCAAGCGACAGACTCA) was used as a reference for quantitative expression analysis. The primers used in qRT-PCR are listed in Table 1.

Transcriptome sequencing and assembly
To comprehensively investigate the transcriptome and gene expression profiles of D. grandiflorum under normal and low temperature stress, four cDNA sample from leaves of chrysanthemum seedlings were prepared and sequenced using Illumina HiSeq 4000 platform. An overview of the RNA-Seq reads derived from the four libraries was presented in Table 2. After the low-quality reads were removed, 25.69Gb clean reads were obtained with an average of 6.42Gb (21.5 million) reads for each sample, and the percentage of Q30 base in each sample was not less than 96.33% (Additional file 1: Figure S1).
With the Trinity program [10], a total of 199,754 transcripts were obtained from the clean reads with an N50 length of 1321 bp and a mean length of 878.8 bp. Among them, 93,837 unigenes were generated with an average length of 719.19 bp. The length distributions of transcripts and unigenes are given in Table 3. The clean data of each sample was compared with the transcript or unigene library. The proportion of mapped reads per library ranged from 71.71% to 73.77%, and unique mapped reads ranged from 47.49% to 51.43% (Table 2).

Gene annotation and functional classification
To predict and analyze the function of the unigenes, we carried out functional annotation by using BLAST against multiple databases such as Nr, Swiss-Prot, GO (Fig. 1a).  GO assignments system was used to classify the possible functions of chrysanthemum genes. A total of 19,534 unigenes were successfully annotated and classified into 52 functional groups of three major GO categories (Cell Component (CC), Molecular Function (MF), and Biological Processes (BP)) ( Fig. 1b). The top three GO terms for classified genes were 'cell' (8319), 'cell part'(8319) and 'organelle' (6248) for cell component category; 'catalytic activity' (11,009), 'binding' (9622) and 'transporter activity' (1273) for molecular function; and 'metabolic process' (13,917), 'cellular process' (11,378) and 'single-organism process' (9623) for biological processes. In addition, genes involved in 'membrane' , 'response to stimulus' and 'biological regulation' were also highly represent.

Analysis of potential differentially expressed genes (DEGs)
The clean reads from each sample were compared to the unigene library, and gene expression levels were estimated using RSEM (RNA-Seq by Expectation Maximization) according to the results of comparison. FPKM (fragments per kilobase of transcript per million fragments mapped) values were used to indicate the expression abundance of unigenes.
From three comparison, including CP1 (T01 vs T02), CP2 (T01 vs T03) and CP3 (T01 vs T04), a large number of DEGs were identified. The number of DEGs detected was as follows: CP1 3028 (1846 up-and 1182 downregulated), CP2 3906 (2799 and 1107) and CP3 3245 (2270 and 975) (Fig. 2b). We found that more DEGs up regulated in treatment T03 which underwent CA than the treatment T04 without CA. Moreover, the venn diagram was used for analyses of unique and overlapping sets of DEGs in each treatment. As shown in Fig. 2a, there were 641 genes differentially expressed in all comparisons, indicating that these genes involved in the process of plants response to both chilling and freezing stresses, which were of great value for low temperature tolerance study. We also identified 1166, 1086 and 1324 DEGs expressed only in T02, T03 and T04, respectively. Thus it can be seen that there were 1166 genes only playing roles in the process of CA, not working in freezing period, and 1324 genes involved in response to rapid freezing. Of total 6557 DEGs, 4378 genes were successfully annotated with eight databases. The number of annotated DEGs in each treatment is shown in Table 4.

KEGG pathway analysis for DEGs
KEGG annotation was used for DEGs from the above comparisons. In the CP1, 493 DEGs were assigned to the KEGG database involved 108 pathways; for CP2, 519 DEGs were assigned to 113 pathways; and for CP3, 386 DEGs assigned to 105 pathways. The most reliable top ten significantly enriched pathways of DEGs under low temperature treatments were represented as Table 5.
Integrate these three comparisons, the major pathways involved in the cold response mechanism were 'Starch and sucrose metabolism' , 'Galactose metabolism' , 'Glycerophospholipid metabolism' , 'Biosynthesis of unsaturated fatty acids' , 'Plant hormone signal transduction' and 'Plantpathogen interaction' (Additional file 1: Figure S2). As shown in Table 6, for first four biological metabolism and synthesis pathways, the number of DEGs in CP1 was more than CP2 or CP3; but for last two signal transduction related pathways, DEGs number of CP1 was less than CP2 or CP3.
Protein kinase (PK) and transcription factors (TFs) responding to low temperature A total of 96, 131, and 129 DEGs encoding protein kinases belonging to 23 major families were found in CP1, CP2 and CP3 treatments, respectively. Among them, 51 (14 up-and 37 down-regulated), 71 (50 and 21) and 80 (71 and 9) receptor-like kinases (RLKs) were identified. And there were 4 (2 up-and 2 downregulated), 8 (7 and 1) and 10 (10 and 0) genes of CP1, CP2 and CP3 involved in the mitogen-activated protein kinase (MAPK) cascade. We found 3, 12 and 9 DEGs of calcium-dependent protein kinase (CDPK) family and 7, 5 and 4 DEGs of CBL-interacting protein kinase (CIPK) family were all up regulated in CP1, CP2 and CP3. Moreover, we identified 13, 12 and 11 protein phosphatase (PP) differentially expressed in three treatments. The details of these PK and PP genes are shown in Additional file 2: Table S1. These results indicated that protein phosphorylation may be induced by low temperature in chrysanthemum. In addition, it was interesting that there were more up-regulated PK in CP2 and CP3 than CP1 (Fig. 3a).
In this study, totally 123 (84 up-and 39 downregulated), 190 (157 and 33) and 193 (172 and 21) transcription factor (TF) genes were identified in CP1, CP2 and CP3. Obviously, the number of up-regulated expressed transcription TFs in CP2 and CP3 was greater than that of CP1 (Fig. 3b). A large number of TF families were involved in low temperature response. Generally, the major TF families were AP2/ERF, WRKY, bHLH, MYB, MYB-related, and NAC. In treatment CP1, the AP2/ERF family was the largest group (18%), followed by the bHLH family (8%) and the GARP-G2-like family (6%); for CP2, the top three family were AP2/ERF (17%),  Figure S3). Of AP2/ERF genes, 17, 29 and 33 were up-regulated in CP1, CP2 and CP3, respectively. And the members of WRKY and MYB families were more active in freezing condition rather than in chilling condition.

Gene related to the membrane system and osmoregulation
The increase of membrane lipids unsaturation level can increase the fluidity of membrane, accordingly improving the cold resistance of plants [19]. In our study, we identified 12 lipid-transfer proteins (LTPs) and 7 fatty acid desaturase (FADs) genes. Among them, 5 LTPs and 3 FADs were up-regulated in CP1, 7 LTPs and 5 FADs up-regulalted in CP2, and only 4 LTPs and 1 FADs were up-regulated in CP3. The details of genes related to the biofilm system are shown in Table 7.

Verification of DEGs using qRT-PCR
To further verify the expression levels of genes obtained from Illumina RNA-Seq platform, we have selected ten genes induced by low temperature to perform qRT-PCR, including genes encoding protein   Analyses of phenotypic and physiological changes of chrysanthemum at different low temperatures.
As shown in Fig. 5a, the phenotypes of chrysanthemum seedlings under chilling stress (4°C) was not significantly different from those under normal conditions, but only the youngest leaves of individual seedlings appeared dying. And under freezing stress (− 4°C), it was clear that the seedlings without prior CA (T04) were more wilted than those with CA (T03). Under low temperature, the membrane permeability of chrysanthemum leaves increased and electrolyte leakage occurred. The degree of electrolyte extravasation was higher in chrysanthemum under freezing conditions than that under the chilling condition (Fig. 5c). And the electrolyte conductivity of T03 was lower than T04, suggesting that CA could increase the tolerance of chrysanthemum to freezing stress. In order to further explore the oxidative damage of chrysantmeum seedlings, we carried out histochemical staining experiments on the leaves. As a whole, chrysanthemum leaves were subjected to more severe oxidative stress under the conditions of freezing than chilling; and under freezing conditions, seedling with prior CA (T03) accumulated less H 2 O 2 and O 2 − than seedlings without CA (T04) (Fig. 5b). Furthermore, we determined activities of antioxidant enzymes (POD and CAT) in   chrysanthemum. Under chilling stress, POD activity of seedlings was not significantly different than normal condition, but CAT activity increased obviously. Under freezing conditions, both POD and CAT activities significantly increased, and activities of T03 were remarkably higher than T04 (Fig. 5d, e). Moreover, we measured the changes in the osmoregulation substances of chrysanthemum at low temperatures. Under chilling stress, both proline and soluble sugar contents of chrysanthemum were increased. And under freezing condition, proline and soluble sugar contents of T03 were substantially increased, whereas T04 was low in proline and soluble sugar contents (Fig. 5f ). The result indicated that CA can induce chrysanthemum osmotic regulation systems to work better under freezing conditions.

Discussion
Low temperature is an important stress factor for plants. It can not only limit the widespread distribution of plants, but also seriously affect the growth and development of them, and cause a series of adverse metabolic reactions [9]. However, 'jinba' exactly is a kind of cutting-chrysanthemum which is difficult to massproduce in winter with normal form. And now, the data available on the molecular basis of chrysanthemum response to cold stress is not limited any longer with the development of high-throughput sequencing.
Based on the quantitative statistics of DEGs in each comparison, we found a lot of genes induced by low temperature and the majority of them involved in upregulation rather than down-regulation (Fig. 2b). From Fig. 2a, we found that some genes may only be induced by freezing and function on protecting chrysanthemum from freezing damage. As expected, more genes were up-regulated in T03 than T04 (Fig. 2b), supporting the statement that CA process may contribute to plants obtaining increased tolerance to freezing conditions [20].

KEGG analysis
As shown in Table 6, when exposed to chilling conditions, DEGs of chrysanthemum associated with biosynthesis and metabolism pathways were enriched more than freezing conditions, suggesting that chilling stress could induce more biosynthesis and metabolism genes than freezing. However, the pathways of 'Plant-pathogen interaction' and 'Plant hormone signal transduction' were greater enriched in freezing condition than chilling. Plant-pathogen interaction pathway involves lots of sensing and signaling genes related to Ca 2+ signal transducion pathway and MAPK cascade, it demonstrated that more genes participated in the signal sensing and transduction pathways were induced by freezing than chilling stress.
In addition, under freezing stress, whether biosynthesis and metabolism related pathways or signal related pathways in CP2 were greater enriched than CP3. It indicated that CA can relieve the inhibition of biosynthesis and metabolism by freezing injury, and strengthen the response of signal transduction pathway to freezing injury.

Low temperature signal transduction
For surviving in low temperature environment, plants developed a series of complex signaling pathways, and the activities of which are modulated by adding or removing phosphate groups [21]. Lots of protein kinases and protein phosphatases participate in the processes of phosphorylation and de-phosphorylation, such as RLKs, MAPKs, CDPKs, CIPKs, protein phosphatase P (PPP), protein phosphatase M (PPM), and protein tyrosine phosphatase (PTP) [22]. Larges of previous experiments implicated that RLKs work on hormone signaling events in response to environmental stimuli [23][24][25][26]. In our Fig. 4 The relative expression levels of ten DEGs identified in the comparison CP1 between RNA-Seq and qRT-PCR. The genes relative expression levels were determined by 2 −ΔΔCT as expressed, and were normalized to the expression level of EF1α experiment, lots of RLKs were found both in all treatments, and up-regulated RLKs in freezing stress were more than in chilling stress (Additional file 2: Table S1). The MAPK cascade activated by reactive oxygen species mediates cold stress signaling in Arabidopsis [27]. Plants over-expressing MKK2 increased freezing tolerance, whereas mkk2 null mutant plants were hypersensitive to cold stress [27]. Overexpression of OsMAPK5 conferred cold tolerance in rice seedlings [28]. We respectively identified 2, 8, 10 MAPKs genes up-regulated in CP1, CP2 and CP3, suggesting that MAPK cascade played a significant role in the low temperature stress, especially in freezing stress in chrysanthemum. In addition, protein phosphatase has been shown to be crucial components of MAPKs cascade [21]. Relevantly, several PPP, PPM and PTP differentially expressed in each treatment were identified. CDPKs and CIPKs are two key kinases of the Ca 2+ -signaling pathway, which is known to be involved in responses to cold stimuli [29,30]. And it has been reported that OsCIPK03, OsCIPK12, OsCDPK7 and OsCDPK13 played important roles in conferring cold tolerance in rice [31][32][33]. In this study, 3, 12, 9 CDPK and 7, 5, 4 CIPK genes were differentially expressed in three treatments, respectively, and all of them were upregulated. These genes may provide excellent genetic resources for breeding engineering of low-temperature tolerance chrysanthemum. The different expression of genes involving in the signal cascade mechanism can affect the expression of genes participating in the formation of plant hormones such as ABA, ET, SA and JA. And then, theses hormones may amplify the cascades or initiate some new signaling pathways [34]. In this experiment, ABA, ET and JA are top three responsive hormone signal transduction pathways in chrysanthemum, and there were numbers of DEGs involved in them. ABA is a plant hormone extensively involved in environmental stresses. ABA-responsive element (ABRE) and dehydrationresponsive element (DRE) were two cis-acting elements involved in ABA-mediated gene expression [21]. Besides DEGs encoding three key proteins (ABA receptor, PP2C and SnRK2) of ABA signaling pathway were identified, another one ABF and 11 DRE-binding factors (DREB) were obtained in our study (Additional file 3: Table S2). Both ABF and DREB can activate the down-stream genes expression in response to cold stress [35]. This indicated that low temperature-induced expression of chrysanthemum appeared to be ABA dependent.
Lipid molecules are also involved in signal transduction during low temperature stress. Phosphatidic acid (one kind of membranous secondary messenger molecule) can be produced by phospholipase D (PLD) [36]. Li et al. have shown that the cold-induced freezing tolerance of Atpldd T-DNA knock-out mutant Arabidopsis was impaired, whereas that of AtPLDd-overexpression plants was enhanced [37]. In our study, we totally found 6 PLD DEGs, 3, 4 and 2 of which represented upregulation in CP1, CP2 and CP3, respectively. The upregulation amplitude of PLD genes in CP2 was higher than CP3. This indicated that PLD was closely related to low-temperature tolerance of chrysanthemum, and CA might contribute to improve freezing tolerance of chrysanthemum.

TFs responding to low temperature
When subjected to low temperature stress, TFs could be activated through a series of signal transduction pathways. Then activated TFs bind specifically to the corresponding cis-acting element to activate the expressions of various downstream resistance-related genes, thereby enhancing the cold tolerance of plants [38]. AP2/ERF, WRKY, MYBrelated, MYB, NAC were five major TF classes identified as DEGs in response to low temperature stress. And some members of above TF families have been identified as regulators involved in low temperature stress responses [39][40][41]. Among various reported TFs, DREBs/CBFs played pivotal roles in improving low temperature tolerance of plants [34]. Overexpression of Arabidopsis DREBs enhanced chilling or freezing tolerance in many plant species, and overexpression DREB1 of several plant species in transgenic Arabidopsis also conferred tolerance to freezing stress [34]. In our study, a total of 10 DREBs were upregulated under chilling or freezing stress. Interestingly, there were only 3 WRKY and 3 MYB TFs up-regulated in the chilling treatment, whereas respectively about 17 were up-regulated in freezing conditions. It can be inferred that WRKY and MYB TFs may play vital roles in response to the freezing stress. Moreover, some other TF families were also found to be differentially expressed in our study, such as bHLH, bZIP, C2H2, GRAS and C3H. These results suggested that a large number of TFs were involved in the regulation mechanism of chrysanthemum in response to low temperature stress through different pathways.

Cold resistant genes associated with the biofilm system
When the temperature drops, phase changes firstly occurs in the biofilm membrane lipid, which leads to the increase of membrane permeability and the leakage of electrolyte. Physiologically, Fig. 5c showed the chrysanthemum biofilm system suffered more severe persecution under freezing conditions, especially T04 chrysanthemum. Therefore, we explored the molecular mechanisms of biological membranes. Increasing the proportion of unsaturated fatty acids in membrane lipids is propitious to membrane fluidity, thus decreasing the phase changes temperature of the membrane lipids [42]. Therefore, the content of membrane fatty acid and the fluidity of membrane are closely related to the cold resistance of plants. In molecular terms, we found a total of 12 genes encoding LTPs and 7 genes encoding FADs were differentially expressed from three libraries (CP1, CP2, and CP3). Among them, 7 FADs genes included omega-6 and omega-3 fatty acid desaturases genes, which are mainly involved in processes of unsaturated fatty acids (linoleic acid and α-linolenic acid) biosynthesis. And previous reports showed that increasing omega-3 desaturase genes expression could enhanced resistance of different plants to cold stress [43,44]. Meanwhile, PLD and secretory phospholipase A2 (SPLA2) were also two favorable factors of membrane fluidity [45,46]. Relevantly, we found numbers of DEGs encoding PLD and SPLA2 in three libraries. These results indicated that chrysanthemum may be able to tolerant low temperatures by enhancing the stabilization of cell membrane. In general, genes encoding LTPs, FADs, PLDs and SPLA2 were obviously much more active in chrysanthemum with CA (T03) than in chrysanthemum without CA (T04) ( Table 7) , suggesting that chrysanthemum acquired freezing tolerance from CA process.