Deep-sequencing transcriptome analysis of chilling tolerance mechanisms of a subnival alpine plant, Chorispora bungeana
- Zhiguang Zhao†1, 2,
- Lingling Tan†1,
- Chunyan Dang1,
- Hua Zhang1,
- Qingbai Wu2 and
- Lizhe An1Email author
© Zhao et al.; licensee BioMed Central Ltd. 2012
Received: 23 July 2012
Accepted: 19 November 2012
Published: 21 November 2012
The plant tolerance mechanisms to low temperature have been studied extensively in the model plant Arabidopsis at the transcriptional level. However, few studies were carried out in plants with strong inherited cold tolerance. Chorispora bungeana is a subnival alpine plant possessing strong cold tolerance mechanisms. To get a deeper insight into its cold tolerance mechanisms, the transcriptome profiles of chilling-treated C. bungeana seedlings were analyzed by Illumina deep-sequencing and compared with Arabidopsis.
Two cDNA libraries constructed from mRNAs of control and chilling-treated seedlings were sequenced by Illumina technology. A total of 54,870 unigenes were obtained by de novo assembly, and 3,484 chilling up-regulated and 4,571 down-regulated unigenes were identified. The expressions of 18 out of top 20 up-regulated unigenes were confirmed by qPCR analysis. Functional network analysis of the up-regulated genes revealed some common biological processes, including cold responses, and molecular functions in C. bungeana and Arabidopsis responding to chilling. Karrikins were found as new plant growth regulators involved in chilling responses of C. bungeana and Arabidopsis. However, genes involved in cold acclimation were enriched in chilling up-regulated genes in Arabidopsis but not in C. bungeana. In addition, although transcription activations were stimulated in both C. bungeana and Arabidopsis, no CBF putative ortholog was up-regulated in C. bungeana while CBF2 and CBF3 were chilling up-regulated in Arabidopsis. On the other hand, up-regulated genes related to protein phosphorylation and auto-ubiquitination processes were over-represented in C. bungeana but not in Arabidopsis.
We conducted the first deep-sequencing transcriptome profiling and chilling stress regulatory network analysis of C. bungeana, a subnival alpine plant with inherited cold tolerance. Comparative transcriptome analysis suggests that cold acclimation is not a major chilling tolerance mechanism of C. bungeana. Activation of protein phosphorylation and ubiquitination may confer chilling tolerance to C. bungeana in a more rapid and flexible way than cold acclimation. Such differences may have contributed to the differences in cold tolerance between C. bungeana and Arabidopsis. The results presented in this paper will be informative for gene discovery and the molecular mechanisms related to plant cold tolerance.
KeywordsAlpine plant Chorispora bungeana Chilling tolerance Cold acclimation Transcriptome
Chorispora bungeana Fisch. & C.A. Mey (C. bungeana) is a perennial subnival alpine plant that can survive freezing temperature . In the natural environments where C. bungeana is growing (origin of Urumqi River in Tianshan Mountains, Xinjiang Autonomous Region, China), snowing and hailing often occur during favorable growing seasons, and air temperature fluctuates frequently ranging from above 20°C to below −10°C. C. bungeana in local environment can survive, grow and flower even in snow. Our previous studies performed at physiological and molecular levels showed that this plant has strong cold (chilling and freezing) tolerance [1–6]. However, little is known about its tolerance mechanisms, if any, distinguishing C. bungeana from other tropical or temperate plants.
Not all plants are always ready to tolerate freezing temperatures. However, studies have shown many plants are tolerant of freezing temperature after exposure to non-freezing low temperature, a phenomenon called cold acclimation [7, 8]. In such a process, various physiological and biochemical changes occur in plant cells, which may confer subsequent acquired chilling and freezing tolerance to plants. For example, during cold acclimation, plants accumulate compatible solutes such as sucrose, raffinose and proline [9–12]; membrane compositions and behaviors are changed [13–16]; and the biosynthesis pathways of secondary metabolites such as flavonoids are activated [17, 18].
The physiological and biochemical changes during plant cold acclimation result mainly from expression changes of cold-responsive (COR) genes. A large number of studies demonstrate that gene expression changes occur in a wide range of plant species in cold responses, and it is believed that differences in COR gene expressions contribute to differences in plant cold tolerance. For example, considerable differences in the members of COR genes were found in Solanum commersonii and Solanum tuberosum, which are closely related species that differ in cold acclimation abilities .
The expressions of COR genes in plant cold responses are under the control of some key transcription factors (TFs). The best characterized TFs involved in plant cold responses are a class of AP2/EFR TFs known as DREB/CBF [20–23], which regulate COR gene expressions by binding to the DRE/CRT cis-elements in the promoter regions of COR genes. In Arabidopsis, there are three major CBFs - CBF1, CBF2 and CBF3 (also known as DREB1b, DREB1c, and DREB1a, respectively) . Constitutive expression of CBF1 and CBF3 can enhance freezing tolerance in non-acclimated Arabidopsis . Moreover, by studying the interactions with CBFs pathway, the roles of some cellular or environmental factors, such as calcium , light , and circadian rhythm , in plant cold tolerance are revealed. Nonetheless, CBFs may not represent all TFs that regulate the expressions of COR genes and confer cold tolerance to plants. Although CBF over-expression increases the freezing tolerance of Arabidopsis, potato  and poplar , it does not increase the freezing tolerance of tomato  and rice . Besides CBFs, some other TFs, such as ZAT12 and RAV1 [33, 34], are also discovered to regulate the expressions of COR genes.
Given the importance of COR genes in plant cold tolerance, studying the cold responses at transcription level may be a key step to identify specific tolerance mechanisms of plants. During the last two decades, numerous studies were carried out to reveal the transcriptional regulatory network of plants in response to cold stress. However, most of the studies were performed with Arabidopsis and others were conducted with crops such as Brassica napus, rice , barley  and potato . Some studies were performed with species adapted to arctic or alpine cold environments, such as Draba [38, 39] and Oxytropis , suggesting that plants may adapt to cold environments with different strategies and COR genes. However, due to lack of reference genome sequence, such studies are relatively few. Sequencing the genome of Coccomyxa subellipsoidea from the Antarctic suggested that gene losses and gains may contribute to low temperature adaptations , highlighting the importance of studying cold tolerance at whole genome or transcriptome level. Recently, the development of high-throughput deep-sequencing technologies makes it possible to study gene expressions at whole genome level without prior knowledge about reference genome sequence. In this study, we used Illumina deep-sequencing technology to study the transcriptome profiles of chilling-treated seedlings of C. bungeana.
C. bungeana is a Cruciferae species closely related to Arabidopsis. Our previous studies showed that the callus and suspension cells from C. bungeana were ready to endure freezing temperature (−4°C) without cold acclimation [3, 6]. The aim of this study is to examine what kinds of mechanisms contribute to the specific cold tolerance of C. bungeana. Our results showed a complicated regulatory network of C. bungeana responding to chilling stress. By comparative transcriptome analysis, a large number of common chilling responding processes, including a newly found karrikins responding process, were found in both C. bungeana and Arabidopsis. Furthermore, our results implied the differences between C. bungeana and Arabidopsis in cold acclimation and TF regulation networks. Importantly, our results suggested that protein phosphorylation and ubiquitination might serve as rapid and flexible mechanisms for cold tolerance regulations in C. bungeana.
Results and discussion
Sequencing and de novo assembly of C. bungeana transcriptome
Statistics of deep-sequencing
Total nucleotides (nt)
Statistics of the assembly (unigene number and percentage) with the Trinity method
> = 2000nt
Length of all Unigene (nt)
Statistics of the assembly (unigene number and percentage) with the SOAPdenovo software
> = 2000nt
Length of all Unigene (nt)
Functional annotation of all the unigenes of C. bungeana
Expression analysis, differential expression genes (DEGs) identification and qPCR verifications
The expressions of unigenes were analyzed with DEGseq R package. Firstly, we tried to identify DEGs by applying screening thresholds of 2 fold changes and Benjamini q value <0.001. We got 12,808 DEG candidates out of 52,753 expressed unigenes (Additional file 2). However, when we verified the expressions of the top 10 up-regulated and down-regulated unigenes by RT-PCR and qPCR, only 3 of them were amplified and none of them showed up or down-regulated trends in chilling-treated seedlings (data not shown). In addition, we found that 80% and 90% of the top 200 up and down-regulated unigenes presented only in one sample’s RNA-seq data, respectively. PCR amplification failures of the selected sequences suggested that such genes were most likely to be the artifacts of de novo assembly.
Top 50 up-regulated unigenes of C. bungeana by chilling stress. The homologs of Arabidopsis genes were presented for functional description of unigenes
log2 (Fold change)
cold-regulated 15a (COR15A)
ABA REPRESSOR1 (ABR1)
cytochrome P450, family 94, subfamily B, polypeptide 1 (CYP94B1)
Plant invertase/pectin methylesterase inhibitor superfamily protein
Integrase-type DNA-binding superfamily protein
mitogen-activated protein kinase kinase kinase 13 (MAPKKK13)
2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily protein
Haloacid dehalogenase-like hydrolase (HAD) superfamily protein
chlorophyllase 1 (CLH1)
cytochrome P450, family 94, subfamily B, polypeptide 1 (CYP94B1)
heat shock transcription factor A1E (HSFA1E)
PYR1-like 11 (PYL11)
UDP-glucosyl transferase 73B3 (UGT73B3)
Stigma-specific Stig1 family protein
gibberellin 2-oxidase 6 (GA2OX6)
myb domain protein 108 (MYB108)
S-adenosyl-L-methionine-dependent methyltransferases superfamily protein
O-Glycosyl hydrolases family 17 protein
EARLY LIGHT-INDUCIBLE PROTEIN 2 (ELIP2)
anthranilate synthase beta subunit 1 (ASB1)
LOW-TEMPERATURE-INDUCED 65 (LTI65)
Integrase-type DNA-binding superfamily protein
purine permease 18 (PUP18)
cytochrome P450, family 709, subfamily B, polypeptide 2 (CYP709B2)
Late embryogenesis abundant (LEA) hydroxyproline-rich glycoprotein family
receptor like protein 39 (RLP39)
disease resistance family protein / LRR family protein
UDP-glucosyl transferase 75B2 (UGT75B2)
UDP-glycosyltransferase 74 F1 (UGT74F1)
Integrase-type DNA-binding superfamily protein
FAD-binding Berberine family protein
Top 50 down-regulated unigenes of C. bungeana by chilling stress. The homologs of Arabidopsis genes were presented for functional description of unigenes
Protein of unknown function (DUF1442)
GATA transcription factor 15 (GATA15)
BARELY ANY MERISTEM 3 (BAM3)
ENHANCER OF ATNSI ACTIVITY (ENA)
Heavy metal transport/detoxification superfamily protein
S-adenosyl-L-methionine-dependent methyltransferases superfamily protein
Polynucleotidyl transferase, ribonuclease H-like superfamily protein
TORNADO 1 (TRN1)
histone H2A 11 (HTA11)
Protein of unknown function (DUF295)
lon protease 1 (LON1)
TOO MANY MOUTHS (TMM)
CBS domain-containing protein with a domain of unknown function (DUF21)
alpha/beta-Hydrolases superfamily protein
High throughput deep-sequencing is a powerful tool for DEGs screening, especially for species without available genomic information [45, 47, 48]. However, since Illumina sequencing is highly sensitive to templates presented in DNA samples, some traced transcripts or contaminants can be sequenced in one sample but not in other samples. This will have huge effects on the results of de novo assembly and increase false positive rate in DEGs identification. One strategy to reduce the false positive results is to set up biological repeats for sequencing and increase sequencing depth, but it will greatly increase the experimental costs. In this study, by simply applying an additional threshold (RPKM > =1) for DEGs screening without increasing costs, we got a high quality (confirmed by qPCR) list of chilling regulated DEGs.
GO network analysis of up-regulated DEGs of C. bungeana in response to chilling stress and comparison with Arabidopsis
Since both C. bungeana and Arabidopsis are Cruciferae species, it is more reliable to use the well-established GO and KEGG annotation systems of Arabidopsis to analyze the functions of C. bungeana DEGs. GO term and KEGG pathway enrichment analysis of DEGs were conducted with BiNGO , a Cytoscape plugin assessing overrepresentation of ontologies in biological networks, using the list of all unigenes with a minimal RPKM of 1 in both sequencing libraries as a reference set. To compare the chilling responding network of C. bungeana with Arabidopsis, the networks of chilling-regulated DEGs of Arabidopsis were constructed using previously published RNA-seq and microarray data (referred to ATH-SR and ATH-MA, respectively; see Methods for details).
Over-representative GO terms* in chilling-treated C. bungeana and Arabidopsis
GO functional description
response to karrikin
response to wounding
response to chitin
regulation of transcription, DNA-dependent
response to water deprivation
response to abscisic acid stimulus
response to cold
defense response to fungus
response to oxidative stress
response to salt stress
response to UV-B
flavonoid biosynthetic process
sequence-specific DNA binding transcription factor activity
Besides abscisic acid  and chitin responses , which were known to be involved in cold tolerance of plants, the biological process “response to karrikin” was found to be a common response to chilling stress in both C. bungeana and Arabidopsis. To our knowledge, no previous study reported the involvement of karrikins in cold tolerance of plants. Karrikins are a new group of plant growth regulators discovered in smoke that can stimulate seed germination . The biological and molecular functions of karrikins are largely unknown at present. Our results suggested that karrikins might play important roles in chilling tolerance of C. bungeana and Arabidopsis.
Nineteen biological processes were over-represented in chilling-treated C. bungeana but not in Arabidopsis. Nonetheless, it did not mean that such processes were specific to chilling responses of C. bungeana since most of them, such as salicylic acid [53, 54], jasmonic acid , and immune response , were reported to be involved in chilling response of Arabidopsis or other plants. However, two processes, “protein phosphorylation” and “protein autoubiquitination”, should be emphasized. Post-translational modifications of pre-existing proteins are believed to be a rapid pathway to get tolerance in plant responses to chilling stress and have important roles in plant cold acclimation . In alfafa, low temperature lead to rapid inhibition of PP2A activity, and in turn lead to phosphorylation of proteins involved in cold tolerance acquisitions [56, 57]. Transcriptional activation of genes of several kinase families were also found under low temperature stress, such as MAP kinase family genes MKK2, OsMEK1 and OsMAP1, CDPK family genes OsCDPK7[60, 61], OsCDPK13 and PaCDPK1, and CIPK family genes CIPK3 and CIPK7. Although many studies reported that certain protein kinases were activated and their transcriptional expression increased in response to cold stress, few studies reported that the expressions of protein kinases as a whole increased at transcriptome level. In our study, a large number of genes whose products were involved in protein phosphorylation were over-represented in chilling up-regulated DEGs in C. bungeana. Given the habitats of C. bungeana, in which the daytime temperatures fluctuate frequently and during almost the whole plant growing seasons, our results suggest that protein phosphorylation may be an important mechanism for rapid and flexible regulation of cold tolerance of C. bungeana.
Chilling up-regulated unigenes annotated with ubiquitination function
ubiquitin ligase, XB3 ortholog 2 in Arabidopsis thaliana (XBAT32)
ubiquitin ligase, XB3 ortholog 2 in Arabidopsis thaliana (XBAT32)
U-box domain E3 ubiquitin ligase protein, plant U-box 22 (PUB22)
U-box domain E3 ubiquitin ligase protein, plant U-box 23 (PUB23)
U-box domain E3 ubiquitin ligase protein, plant U-box 24 (PUB24)
A20/AN1-like zinc finger family protein
zinc finger (C3HC4-type RING finger) family protein
Chilling up-regulated TFs overlapped in C. bungeana and Arabidopsis
All gene symbols
Homologous to the flowering-time gene CONSTANS.
CZF1; SZF2; ZFAR1
Encodes a transcription factor that specifically binds to DRE/CRT cis elements (responsive to drought and low-temperature stress)
Member of Heat Stress Transcription Factor (Hsf) family
Related to AP2.7 (RAP2.7)
Encodes a zinc finger protein involved in high light and cold acclimation
Myb-like transcription factor that regulates hypocotyl growth by regulating free auxin levels in a time-of-day specific manner.
TEOSINTE BRANCHED 1; TCP2
TEOSINTE BRANCHED 1, cycloidea and PCF transcription factor 2 (TCP2)
Pseudo response regulator involved in the generation of circadian rhythms.
Encodes a member of WRKY Transcription Factor
Member of the plant WRKY transcription factor family
Pathogen-induced transcription factor
B-box type zinc finger family protein
Zinc finger (CCCH-type) family protein
Sequence-specific DNA binding transcription factors
B-box type zinc finger protein with CCT domain
Over-representative GO terms* in chilling stressed C. bungeana but not in Arabidopsis
Stimulus responses related:
response to salicylic acid stimulus
response to jasmonic acid stimulus
defense response by callose deposition in cell wall
defense response to bacterium
innate immune response
cellular response to abiotic stimulus
response to starvation
jasmonic acid biosynthetic process
indole glucosinolate metabolic process
aromatic amino acid family catabolic process
glutamine family amino acid catabolic process
camalexin biosynthetic process
L-phenylalanine metabolic process
recognition of pollen
organic substance transport
sequence-specific DNA binding
iron ion binding
electron carrier activity
structural constituent of cell wall
oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen
nutrient reservoir activity
protein serine/threonine kinase activity
carbon-nitrogen lyase activity
KEGG pathway analysis of up-regulated DEGs of C. bungeana in response to chilling stress and comparison with Arabidopsis
All over-represented pathways in C. bungeana, regardless of whether they were enriched in Arabidopsis, had proved to be important in plant cold tolerance. For instance, circadian rhythm regulates the expression of CBFs [28, 69], the core identified TFs that involved in plant cold tolerance. As another example, under chilling stress, plants preferentially accumulate polyunsaturated fatty acids such as linoleic and linolenic fatty acids [70–72], and genetically increasing of unsaturated fatty acids or lipids could enhance cold tolerance of transgenic plants, probably by maintaining membrane fluidity under cold stress [73, 74]. Our previous findings indicated that cold tolerance of C. bungeana was correlated with changes in membrane lipids and membrane-associated enzymes . Under chilling treatment, the proportion of unsaturated fatty acid in the plasma membrane increased significantly in callus of C. bungeana. Paralleling to these results, KEGG analysis in this study showed that unigenes involved in "alpha-Linolenic acid metabolism" were enriched significantly in chilling up-regulated DEGs, suggesting that lipid metabolism, especially linolenic acid metabolism, might play a role in chilling tolerance of C. bungeana.
GO network analysis of down-regulated DEGs of C. bungeana in response to chilling stress and comparison with Arabidopsis
C. bungeana is a perennial subnival alpine plant with high capacity of chilling and freezing resistance. In recent years, much effort has been taken in our research group to reveal the cold tolerance mechanisms of this plant at physiological and molecular levels. In this paper, we provide the first study on the transcriptome of chilling stressed seedlings of C. bungeana. We got 54,870 assembled unigenes using the Trinity de novo assembly method, and a number of chilling regulated genes were identified, providing useful resources for gene mining to improve cold tolerance of plants. Furthermore, the comparison of the functional networks of chilling regulated genes in C. bungeana and Arabidopsis provided informative results, which could help us tell the differences in cold tolerance mechanisms between C. bungeana and Arabidopsis. We found that karrikins might be new plant growth regulators involved in chilling tolerance of plants. Although gene expressions at the transcriptional level were stimulated by chilling in both C. bungeana and Arabidopsis, their activation networks were different as suggested by TFs analysis. Cold acclimation mechanism may be weak in or absent from C. bungeana because of lack of some CBFs orthologs. Alternatively, protein phosphorylation and ubiquitination may serve as more rapid and flexible cold tolerance mechanisms for C. bungeana to adapt to the harsh cold environments.
Plant material, growth conditions and treatments
Plant regeneration of C. bungeana via somatic embryogenesis was performed as described by Wang et al. . Callus was induced from matured seeds of C. bungeana on MS medium containing 4.0 mg l-1 GA3, 2.0 mg l-1 NAA, and 2.0 mg l-1 2,4-D. Seedlings were regenerated from callus on MS medium containing 3% sucrose in about 3 weeks. Regenerated plants were transferred to new MS medium containing 3% sucrose and grown at 22°C with a 14 h photoperiod under 80 μmol m-2 s-1 fluorescent light for further 7 days before treatments. For each treatment, ten plants (roots, shoots and leaves) were randomly pooled and treated in MS liquid medium containing 3% sucrose at 22°C or 2°C. Chilling stress was initiated 4 hours after dawn (zeitgeber time 4; ZT4). Upon the treatment time reaching 24 hours, both control and chilling stressed samples were collected at the same time point and frozen immediately with liquid nitrogen.
RNA extraction, cDNA library construction and RNA sequencing
For RNA sequencing, total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). The quality of total RNA was checked using the NanoDrop Spectrometer (ND-1000 Spectrophotometer, Peqlab) and the Agilent 2100 Bioanalyzer (RNA Nano Chip, Agilent). High quality RNA samples (20 μg each) were sent to Beijing Genomics Institute (BGI, Shenzhen) for cDNA libraries construction and sequencing using Illumina HiSeq™ 2000. The cDNA library construction method and Illumina deep-sequencing processes were the same as described by Xu et al. .
De novo assembly and sequences clustering
The Trinity method  was used for de novo assembly of the clean reads to generate Trinity unigenes, with optimized k-mer length of 25. Then, the Trinity unigenes of both libraries were clustered with TGICL software  to get sequences (final unigenes) that cannot be extended on either end. De novo assembly was also conducted with SOAPdenovo software  with optimized k-mer length of 41.
Files containing the raw read sequences and their quality scores are available from the National Center for Biotechnology Information (NCBI) Short Read Archive with the accession number: SRA054354. The Trinity unigenes have been deposited in the Transcriptome Shotgun Assembly Sequence Database (TSA) at NCBI [GenBank: JW988067-JW999999, KA000001-KA089547].
Expression analysis and identification of differentially expressed genes (DEGs)
Clean reads were mapped back to assembled unigenes with SOAPaligner (version 2.21)  allowing maximum 2 mismatches. The reads with unique best hits were counted for each unigene. The expression level of C. bungeana unigene was normalized by the number of RPKM (reads per kilobase exon region per million mapped reads) . Since Illumina sequencing method is highly sensitive, we only used a subset of unigenes which presented in both sequencing libraries with a minimal RPKM of 1 for DEGs analysis. Unigene expressions were analyzed using DEGseq R package  with MARS method. Chilling-regulated DEGs were identified with Benjamini q < 0.001  and normalized fold change > =2.
For comparisons, two public available data sets of Arabidopsis were used in our study. One data set (referred to ATH-SR, means Arabidopsis short reads) was RNA sequencing data downloaded from NCBI Sequence Read Archive (SRA) database (http://www.ncbi.nlm.nih.gov), including a chilling-treated sample (4°C; SRA accession: SRX006193) and a control (21°C; SRA accession: SRX006704) sample . After removing low quality reads (polyA/T/G/C sequences) and trimming off four NTs of both ends, all clean reads (28 NTs long) were mapped to Arabidopsis cDNAs (TAIR10) with SOAPaligner. DEGs identification was the same as described above. The DEGs and indentified gene with RPKM > =1 were listed in Additional file 6.
The other data set (referred to ATH-AR, means Arabidopsis array) was Affimetrix microarray data set (Expression Set: ME00325)  downloaded from TAIR (http://www.arabidopsis.org). Only cel files for 4 chilling-treated samples (2 for roots and 2 for shoots, 24-hour chilling-treated) and 4 control samples were used here. The cel files were imported into R and analyzed with Affy package . Root and shoot arrays were analyzed separately. Probes expressed in all root or shoot arrays were considered to be presented probes (by mas5 present calls). Differential expressed probes were identified using mas5 method of with FDR corrected p < 0.05 and fold change > =2 and mapped to Arabidopsis transcripts. The gene lists of roots and shoots were combined together to get chilling regulated DEGs and all expressed genes for further analysis (Additional file 7).
We used two methods for functional categorization of unigenes.
To get general gene ontology (GO) annotations for all unigenes, sequences longer than 200 bp were aligned to three public databases (NR, Swiss-Prot and KEGG) by BLASTX with E-value < =1e-5. The GO annotations for the top blast hits were retrieved with Blast2GO program  and used to annotate the C. bungeana transcripts. GO functional classification was performed by WEGO website tool .
For GO terms and KEGG pathways enrichment analysis, we used the Arabidopsis annotation systems. Briefly, the sequences of all unigenes were aligned against Arabidopsis peptide database (TAIR10) using BLASTX program with E-value < =1e-5. The top blast hits were considered to be putative orthologous genes (POGs). Then the C. bungeana unigenes were annotated with GO (downloaded from TAIR) and KEGG annotations (ath00001.keg, from http://www.kegg.jp/) for Arabidopsis POGs, respectively. The ontology (GO and KEGG) enrichment was analyzed with BiNGO plugins  for Cytoscape , using hypergeometric test for statistical analysis. For p value correction, we used rigorous Bonferroni correction method. The cutoff p value after correction was 0.05. For ATH-SR dataset, since the stressed sample was pooled from seedlings subjected to various periods of chilling-treated (1, 2, 5, 10, 24 hours of stressed) , the expressions of DEGs specific to a certain stage might have been “normalized”. Therefore, to get more information, we used FDR method instead of Bonferronic method for p value correction to find over-representative terms with BiNGO.
Quantitative real-time PCR (qPCR)
The gene-specific primers for real-time PCR analysis were designed using Primer Premier (version 5.0) software (PREMIER Biosoft). The specifities of primer pairs were confirmed by BLASTN with non-redundant unigene set of C. bungeana transcripts and the PCR products were checked by agrose electrophoresis to ensure single band amplifications. The primer sequences for all unigenes used in this study were listed in Additional file 8.
For qPCR analysis, total RNAs were extracted from control or chilling stressed C. bungeana seedlings (two biological repeats) with TRIZOL reagent and treated (20 μg RNA) with 1U DNase (TAKARA, Japan). cDNA was transcribed reversely from 1 μg of DNase-treated RNA with 200U M-MLV reverse transcriptase (Promega, USA) and analyzed with Platinum SYBR green qPCR supermix-UDG reagents (Invitrogen).
Before quantification of unigenes, the geNorm method was applied to select stable expressed unigenes in the four samples as reference genes . A total of 8 candidate reference unigenes were selected for reference genes screening, including an ACTIN2 ortholog, 3 unigenes showed stable expression levels in RNA-seq analysis and the other 4 unigenes were orthologs of recommended Arabidopsis reference genes . The information of reference gene candidates and the geNorm analysis results were shown Additional file 8. Three unigenes (CBT10872/AT3G60800, CBT28565/AT5G27630 and CBT12464/AT2G28390) expressed most stably in control and chilling-treated samples were selected and used for all qPCR analysis.
qPCR analysis was performed with three technical repeats for each sample. The relative expression levels of unigenes were normalized with the three selected reference genes with Pfaffl method [86, 88].
Availability of supporting data
The data sets supporting the results of this article are available in the NCBI GenBank repository, [http://www.ncbi.nlm.nih.gov/sites/nuccore?term=104929[BioProject], and in the NCBI SRA repository, [http://www.ncbi.nlm.nih.gov/sra?term=SRA054354].
CRT binding transcription factor
Differentially expressed gene
Kyoto Encyclopedia of Genes and Genomes
Putative orthologous gene
Quantitative real-time PCR
This work was funded by the Major Project of Cultivating New Varieties of Transgenic Organisms (2009ZX08009-029B), State Key Laboratory of Frozen Soil Engineering (SKLFSE201004) and the National Science Foundation of China (31070353).
- Fu XY, Chang JF, An LZ, Zhang MX, Xu SJ, Chen T, Liu YH, Xin H, Wang JH: Association of the cold-hardiness of Chorispora bungeana with the distribution and accumulation of calcium in the cells and tissues. Environ Exp Bot. 2006, 55: 282-293. 10.1016/j.envexpbot.2004.11.009.View ArticleGoogle Scholar
- Liu YJ, Zhao ZG, Si J, Di CX, Han J, An LZ: Brassinosteroids alleviate chilling-induced oxidative damage by enhancing antioxidant defense system in suspension cultured cells of Chorispora bungeana. Plant Growth Regul. 2009, 59: 207-214. 10.1007/s10725-009-9405-9.View ArticleGoogle Scholar
- Shi Y, An L, Zhang M, Huang C, Zhang H, Xu S: Regulation of the plasma membrane during exposure to low temperatures in suspension-cultured cells from a cryophyte (Chorispora bungeana). Protoplasma. 2008, 232: 173-181. 10.1007/s00709-008-0291-1.PubMedView ArticleGoogle Scholar
- Guo FX, Zhang MX, Chen Y, Zhang WH, Xu SJ, Wang JH, An LZ: Relation of several antioxidant enzymes to rapid freezing resistance in suspension cultured cells from alpine Chorispora bungeana. Cryobiology. 2006, 52: 241-250. 10.1016/j.cryobiol.2005.12.001.PubMedView ArticleGoogle Scholar
- Ding S, Huang CL, Sheng HM, Song CL, Li YB, An LZ: Effect of inoculation with the endophyte Clavibacter sp. strain Enf12 on chilling tolerance in Chorispora bungeana. Physiol Plant. 2011, 141: 141-151. 10.1111/j.1399-3054.2010.01428.x.PubMedView ArticleGoogle Scholar
- Wu JM, Zhao ZG, Xing H, Guo HP, Li WX, An LZ, Xu SJ, Chen T: Effects of freezing on plasma membrane H + −ATPase of the callus from Chorispora bungeana. Biol Plant. 2007, 51: 229-234. 10.1007/s10535-007-0046-4.View ArticleGoogle Scholar
- Guy CL: Cold accelimation and freezing stress tolerance: role of protein metabolism. Annu Rev Plant Biol. 1990, 41: 187-223. 10.1146/annurev.pp.41.060190.001155.View ArticleGoogle Scholar
- Thomashow MF: Plant cold acclimation: freezing tolerance genes and regulatory mechanisms. Ann Rev Plant Physiol Plant Mol Biol. 1999, 50: 571-599. 10.1146/annurev.arplant.50.1.571.View ArticleGoogle Scholar
- Steponkus PL: Cold acclimation of hedera helix: evidence for a two phase process. Plant Physiol. 1971, 47: 175-180. 10.1104/pp.47.2.175.PubMedPubMed CentralView ArticleGoogle Scholar
- Koster KL, Lynch DV: Solute accumulation and compartmentation during the cold acclimation of puma rye. Plant Physiol. 1992, 98: 108-113. 10.1104/pp.98.1.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamata T, Uemura M: Solute accumulation in heat seedlings during cold acclimation: contribution to increased freezing tolerance. Cryo Letters. 2004, 25: 311-322.PubMedGoogle Scholar
- Kushad MM, Yelenosky G: Evaluation of polyamine and proline levels during low temperature acclimation of citrus. Plant Physiol. 1987, 84: 692-695. 10.1104/pp.84.3.692.PubMedPubMed CentralView ArticleGoogle Scholar
- Evert DR, Weiser CJ: Relationship of electrical conductance at two frequencies to cold injury and acclimation in cornus stolonifera Michx. Plant Physiol. 1971, 47: 204-208. 10.1104/pp.47.2.204.PubMedPubMed CentralView ArticleGoogle Scholar
- Miller RW, de la Roche I, Pomeroy MK: Structural and functional responses of wheat mitochondrial membranes to growth at low temperatures. Plant Physiol. 1974, 53: 426-433. 10.1104/pp.53.3.426.PubMedPubMed CentralView ArticleGoogle Scholar
- Pomeroy MK, Raison JK: Maintenance of membrane fluidity during development of freezing tolerance of winter wheat seedlings. Plant Physiol. 1981, 68: 382-385. 10.1104/pp.68.2.382.PubMedPubMed CentralView ArticleGoogle Scholar
- Toivio-Kinnucan MA, Chen HH, Li PH, Stushnoff C: Plasma membrane alterations in callus tissues of tuber-bearing solanum species during cold acclimation. Plant Physiol. 1981, 67: 478-483. 10.1104/pp.67.3.478.PubMedPubMed CentralView ArticleGoogle Scholar
- Taulavuori E, Tahkokorpi M, Taulavuori K, Laine K: Anthocyanins and glutathione S-transferase activities in response to low temperature and frost hardening in Vaccinium myrtillus (L.). J Plant Physiol. 2004, 161: 903-911. 10.1016/j.jplph.2003.12.001.PubMedView ArticleGoogle Scholar
- Korn M, Peterek S, Mock HP, Heyer AG, Hincha DK: Heterosis in the freezing tolerance, and sugar and flavonoid contents of crosses between Arabidopsis thaliana accessions of widely varying freezing tolerance. Plant Cell Environ. 2008, 31: 813-827. 10.1111/j.1365-3040.2008.01800.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Carvallo MA, Pino MT, Jeknic Z, Zou C, Doherty CJ, Shiu SH, Chen TH, Thomashow MF: A comparison of the low temperature transcriptomes and CBF regulons of three plant species that differ in freezing tolerance: Solanum commersonii, Solanum tuberosum, and Arabidopsis thaliana. J Exp Bot. 2011, 62: 3807-3819. 10.1093/jxb/err066.PubMedPubMed CentralView ArticleGoogle Scholar
- Stockinger EJ, Gilmour SJ, Thomashow MF: Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci U S A. 1997, 94: 1035-1040. 10.1073/pnas.94.3.1035.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF: Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998, 280: 104-106. 10.1126/science.280.5360.104.PubMedView ArticleGoogle Scholar
- Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998, 10: 1391-1406.PubMedPubMed CentralView ArticleGoogle Scholar
- Medina J, Bargues M, Terol J, Perez-Alonso M, Salinas J: The Arabidopsis CBF gene family is composed of three genes encoding AP2 domain-containing proteins whose expression Is regulated by low temperature but not by abscisic acid or dehydration. Plant Physiol. 1999, 119: 463-470. 10.1104/pp.119.2.463.PubMedPubMed CentralView ArticleGoogle Scholar
- Shinwari ZK, Nakashima K, Miura S, Kasuga M, Seki M, Yamaguchi-Shinozaki K, Shinozaki K: An Arabidopsis gene family encoding DRE/CRT binding proteins involved in low-temperature-responsive gene expression. Biochem Biophys Res Commun. 1998, 250: 161-170. 10.1006/bbrc.1998.9267.PubMedView ArticleGoogle Scholar
- Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF: Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol. 2000, 124: 1854-1865. 10.1104/pp.124.4.1854.PubMedPubMed CentralView ArticleGoogle Scholar
- Doherty CJ, Van Buskirk HA, Myers SJ, Thomashow MF: Roles for Arabidopsis CAMTA transcription factors in cold-regulated gene expression and freezing tolerance. Plant Cell. 2009, 21: 972-984. 10.1105/tpc.108.063958.PubMedPubMed CentralView ArticleGoogle Scholar
- Franklin KA, Whitelam GC: Light-quality regulation of freezing tolerance in Arabidopsis thaliana. Nat Genet. 2007, 39: 1410-1413. 10.1038/ng.2007.3.PubMedView ArticleGoogle Scholar
- Dong MA, Farre EM, Thomashow MF: Circadian clock-associated 1 and late elongated hypocotyl regulate expression of the C-repeat binding factor (CBF) pathway in Arabidopsis. Proc Natl Acad Sci U S A. 2011, 108: 7241-7246. 10.1073/pnas.1103741108.PubMedPubMed CentralView ArticleGoogle Scholar
- Pino MT, Skinner JS, Park EJ, Jeknic Z, Hayes PM, Thomashow MF, Chen TH: Use of a stress inducible promoter to drive ectopic AtCBF expression improves potato freezing tolerance while minimizing negative effects on tuber yield. Plant Biotechnol J. 2007, 5: 591-604. 10.1111/j.1467-7652.2007.00269.x.PubMedView ArticleGoogle Scholar
- Benedict C, Skinner JS, Meng R, Chang Y, Bhalerao R, Huner NP, Finn CE, Chen TH, Hurry V: The CBF1-dependent low temperature signalling pathway, regulon and increase in freeze tolerance are conserved in Populus spp. Plant Cell Environ. 2006, 29: 1259-1272. 10.1111/j.1365-3040.2006.01505.x.PubMedView ArticleGoogle Scholar
- Hsieh TH, Lee JT, Yang PT, Chiu LH, Charng YY, Wang YC, Chan MT: Heterology expression of the Arabidopsis C-repeat/dehydration response element binding factor 1 gene confers elevated tolerance to chilling and oxidative stresses in transgenic tomato. Plant Physiol. 2002, 129: 1086-1094. 10.1104/pp.003442.PubMedPubMed CentralView ArticleGoogle Scholar
- Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 2003, 33: 751-763. 10.1046/j.1365-313X.2003.01661.x.PubMedView ArticleGoogle Scholar
- Vogel JT, Zarka DG, Van Buskirk HA, Fowler SG, Thomashow MF: Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J. 2005, 41: 195-211.PubMedView ArticleGoogle Scholar
- Fowler S, Thomashow MF: Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell. 2002, 14: 1675-1690. 10.1105/tpc.003483.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaglo KR, Kleff S, Amundsen KL, Zhang X, Haake V, Zhang JZ, Deits T, Thomashow MF: Components of the Arabidopsis C-repeat/dehydration-responsive element binding factor cold-response pathway are conserved in Brassica napus and other plant species. Plant Physiol. 2001, 127: 910-917. 10.1104/pp.010548.PubMedPubMed CentralView ArticleGoogle Scholar
- Yun KY, Park MR, Mohanty B, Herath V, Xu F, Mauleon R, Wijaya E, Bajic VB, Bruskiewich R, de Los Reyes BG: Transcriptional regulatory network triggered by oxidative signals configures the early response mechanisms of japonica rice to chilling stress. BMC Plant Biol. 2010, 10: 16-10.1186/1471-2229-10-16.PubMedPubMed CentralView ArticleGoogle Scholar
- Greenup AG, Sasani S, Oliver SN, Walford SA, Millar AA, Trevaskis B: Transcriptome analysis of the vernalization response in barley (Hordeum vulgare) seedlings. PLoS One. 2011, 6: e17900-10.1371/journal.pone.0017900.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou D, Zhou J, Meng L, Wang Q, Xie H, Guan Y, Ma Z, Zhong Y, Chen F, Liu J: Duplication and adaptive evolution of the COR15 genes within the highly cold-tolerant Draba lineage (Brassicaceae). Gene. 2009, 441: 36-44. 10.1016/j.gene.2008.06.024.PubMedView ArticleGoogle Scholar
- von Meijenfeldt N: Unraveling the cold response in Draba. PhD thesis. 2010, The Netherlands: Universiteit van Amsterdam, The Institute for Biodiversity and Ecosystem Dynamics (IBED) and the Swammerdam Institute for Life Sciences (SILS)Google Scholar
- Archambault A, Stromvik MV: PR-10, defensin and cold dehydrin genes are among those over expressed in Oxytropis (Fabaceae) species adapted to the arctic. Funct Integr Genomics. 2011, 11: 497-505. 10.1007/s10142-011-0223-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Blanc G, Agarkova I, Grimwood J, Kuo A, Brueggeman A, Dunigan DD, Gurnon J, Ladunga I, Lindquist E, Lucas S, et al: The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation. Genome Biol. 2012, 13: R39-10.1186/gb-2012-13-5-r39.PubMedPubMed CentralView ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.PubMedPubMed CentralView ArticleGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.PubMedView ArticleGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25: 1966-1967. 10.1093/bioinformatics/btp336.PubMedView ArticleGoogle Scholar
- Xu DL, Long H, Liang JJ, Zhang J, Chen X, Li JL, Pan ZF, Deng GB, Yu MQ: De novo assembly and characterization of the root transcriptome of Aegilops variabilis during an interaction with the cereal cyst nematode. BMC Genomics. 2012, 13: 133-10.1186/1471-2164-13-133.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakayama K, Okawa K, Kakizaki T, Honma T, Itoh H, Inaba T: Arabidopsis Cor15am is a chloroplast stromal protein that has cryoprotective activity and forms oligomers. Plant Physiol. 2007, 144: 513-523. 10.1104/pp.106.094581.PubMedPubMed CentralView ArticleGoogle Scholar
- Qiu Q, Ma T, Hu Q, Liu B, Wu Y, Zhou H, Wang Q, Wang J, Liu J: Genome-scale transcriptome analysis of the desert poplar, Populus euphratica. Tree Physiology. 2011, 31: 452-461. 10.1093/treephys/tpr015.PubMedView ArticleGoogle Scholar
- Wong MM, Cannon CH, Wickneswari R: Identification of lignin genes and regulatory sequences involved in secondary cell wall formation in Acacia auriculiformis and Acacia mangium via de novo transcriptome sequencing. BMC Genomics. 2011, 12: 342-10.1186/1471-2164-12-342.PubMedPubMed CentralView ArticleGoogle Scholar
- Maere S, Heymans K, Kuiper M: BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21: 3448-3449. 10.1093/bioinformatics/bti551.PubMedView ArticleGoogle Scholar
- Ishitani M, Xiong L, Stevenson B, Zhu JK: Genetic analysis of osmotic and cold stress signal transduction in Arabidopsis: interactions and convergence of abscisic acid-dependent and abscisic acid-independent pathways. Plant Cell. 1997, 9: 1935-1949.PubMedPubMed CentralView ArticleGoogle Scholar
- Yeh S, Moffatt BA, Griffith M, Xiong F, Yang DS, Wiseman SB, Sarhan F, Danyluk J, Xue YQ, Hew CL, et al: Chitinase genes responsive to cold encode antifreeze proteins in winter cereals. Plant Physiol. 2000, 124: 1251-1264. 10.1104/pp.124.3.1251.PubMedPubMed CentralView ArticleGoogle Scholar
- Nelson DC, Riseborough JA, Flematti GR, Stevens J, Ghisalberti EL, Dixon KW, Smith SM: Karrikins discovered in smoke trigger Arabidopsis seed germination by a mechanism requiring gibberellic acid synthesis and light. Plant Physiol. 2009, 149: 863-873.PubMedPubMed CentralView ArticleGoogle Scholar
- Scott IM, Clarke SM, Wood JE, Mur LA: Salicylate accumulation inhibits growth at chilling temperature in Arabidopsis. Plant Physiol. 2004, 135: 1040-1049. 10.1104/pp.104.041293.PubMedPubMed CentralView ArticleGoogle Scholar
- Ding CK, Wang CY, Gross KC, Smith DL: Jasmonate and salicylate induce the expression of pathogenesis-related-protein genes and increase resistance to chilling injury in tomato fruit. Planta. 2002, 214: 895-901. 10.1007/s00425-001-0698-9.PubMedView ArticleGoogle Scholar
- Huang X, Li J, Bao F, Zhang X, Yang S: A gain-of-function mutation in the Arabidopsis disease resistance gene RPP4 confers sensitivity to low temperature. Plant Physiol. 2010, 154: 796-809. 10.1104/pp.110.157610.PubMedPubMed CentralView ArticleGoogle Scholar
- Monroy AF, Sangwan V, Dhindsa RS: Low temperature signal transduction during cold acclimation: protein phosphatase 2A as an early target for cold-inactivation. Plant J. 1998, 13: 653-660. 10.1046/j.1365-313X.1998.00070.x.View ArticleGoogle Scholar
- Monroy AF, Sarhan F, Dhindsa RS: Cold-Induced Changes in Freezing Tolerance, Protein Phosphorylation, and Gene Expression (Evidence for a Role of Calcium). Plant Physiol. 1993, 102: 1227-1235.PubMedPubMed CentralView ArticleGoogle Scholar
- Teige M, Scheikl E, Eulgem T, Doczi R, Ichimura K, Shinozaki K, Dangl JL, Hirt H: The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell. 2004, 15: 141-152. 10.1016/j.molcel.2004.06.023.PubMedView ArticleGoogle Scholar
- Wen JQ, Oono K, Imai R: Two novel mitogen-activated protein signaling components, OsMEK1 and OsMAP1, are involved in a moderate low-temperature signaling pathway in rice. Plant Physiol. 2002, 129: 1880-1891. 10.1104/pp.006072.PubMedPubMed CentralView ArticleGoogle Scholar
- Saijo Y, Hata S, Kyozuka J, Shimamoto K, Izui K: Over-expression of a single Ca2 + −dependent protein kinase confers both cold and salt/drought tolerance on rice plants. Plant J. 2000, 23: 319-327. 10.1046/j.1365-313x.2000.00787.x.PubMedView ArticleGoogle Scholar
- Saijo Y, Kinoshita N, Ishiyama K, Hata S, Kyozuka J, Hayakawa T, Nakamura T, Shimamoto K, Yamaya T, Izui K: A Ca(2+)-dependent protein kinase that endows rice plants with cold- and salt-stress tolerance functions in vascular bundles. Plant Cell Physiol. 2001, 42: 1228-1233. 10.1093/pcp/pce158.PubMedView ArticleGoogle Scholar
- Abbasi F, Onodera H, Toki S, Tanaka H, Komatsu S: OsCDPK13, a calcium-dependent protein kinase gene from rice, is induced by cold and gibberellin in rice leaf sheath. Plant Mol Biol. 2004, 55: 541-552. 10.1007/s11103-004-1178-y.PubMedView ArticleGoogle Scholar
- Tsai TM, Chen YR, Kao TW, Tsay WS, Wu CP, Huang DD, Chen WH, Chang CC, Huang HJ: PaCDPK1, a gene encoding calcium-dependent protein kinase from orchid, Phalaenopsis amabilis, is induced by cold, wounding, and pathogen challenge. Plant Cell Rep. 2007, 26: 1899-1908. 10.1007/s00299-007-0389-5.PubMedView ArticleGoogle Scholar
- Kim KN, Cheong YH, Grant JJ, Pandey GK, Luan S: CIPK3, a calcium sensor-associated protein kinase that regulates abscisic acid and cold signal transduction in Arabidopsis. Plant Cell. 2003, 15: 411-423. 10.1105/tpc.006858.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang C, Ding S, Zhang H, Du H, An L: CIPK7 is involved in cold response by interacting with CBL1 in Arabidopsis thaliana. Plant Sci. 2011, 181: 57-64. 10.1016/j.plantsci.2011.03.011.PubMedView ArticleGoogle Scholar
- Dong CH, Agarwal M, Zhang Y, Xie Q, Zhu JK: The negative regulator of plant cold responses, HOS1, is a RING E3 ligase that mediates the ubiquitination and degradation of ICE1. Proc Natl Acad Sci U S A. 2006, 103: 8281-8286. 10.1073/pnas.0602874103.PubMedPubMed CentralView ArticleGoogle Scholar
- Michael TP, Salome PA, McClung CR: Two Arabidopsis circadian oscillators can be distinguished by differential temperature sensitivity. Proc Natl Acad Sci U S A. 2003, 100: 6878-6883. 10.1073/pnas.1131995100.PubMedPubMed CentralView ArticleGoogle Scholar
- Mikkelsen MD, Thomashow MF: A role for circadian evening elements in cold-regulated gene expression in Arabidopsis. Plant J. 2009, 60: 328-339. 10.1111/j.1365-313X.2009.03957.x.PubMedView ArticleGoogle Scholar
- Fowler SG, Cook D, Thomashow MF: Low temperature induction of Arabidopsis CBF1, 2, and 3 is gated by the circadian clock. Plant Physiol. 2005, 137: 961-968. 10.1104/pp.104.058354.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerloff ED, Richardson T, Stahmann MA: Changes in Fatty acids of alfalfa roots during cold hardening. Plant Physiol. 1966, 41: 1280-1284. 10.1104/pp.41.8.1280.PubMedPubMed CentralView ArticleGoogle Scholar
- Palta JP, Whitaker BD, Weiss LS: Plasma Membrane Lipids Associated with Genetic Variability in Freezing Tolerance and Cold Acclimation of Solanum Species. Plant Physiol. 1993, 103: 793-803.PubMedPubMed CentralGoogle Scholar
- Williams JP, Khan MU, Mitchell K, Johnson G: The effect of temperature on the level and biosynthesis of unsaturated fatty acids in diacylglycerols of Brassica napus leaves. Plant Physiol. 1988, 87: 904-910. 10.1104/pp.87.4.904.PubMedPubMed CentralView ArticleGoogle Scholar
- Kodama H, Hamada T, Horiguchi G, Nishimura M, Iba K: Genetic enhancement of cold tolerance by expression of a gene for chloroplast [omega]-3 fatty acid desaturase in transgenic tobacco. Plant Physiol. 1994, 105: 601-605.PubMedPubMed CentralGoogle Scholar
- Nishida I, Murata N: Chilling sensitivity in plants and cyanobacteria: the crucial contribution of membrane lipids. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 541-568. 10.1146/annurev.arplant.47.1.541.PubMedView ArticleGoogle Scholar
- Wu J, Zhao Z, An L, Liu Y, Xu S, Gao D, Zhang Y: Inhibition of glutathione synthesis decreases chilling tolerance in Chorispora bungeana callus. Cryobiology. 2008, 57: 9-17. 10.1016/j.cryobiol.2008.04.001.PubMedView ArticleGoogle Scholar
- Wang J, An L, Wang R, Yang D, Si J, Fu X, Chang J, Xu S: Plant regeneration of Chorispora bungeana via somatic embryogenesis. In Vitro Cell Dev Biol Plant. 2006, 42: 148-151. 10.1079/IVP2006748.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedView ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26: 136-138. 10.1093/bioinformatics/btp612.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.Google Scholar
- Filichkin SA, Priest HD, Givan SA, Shen R, Bryant DW, Fox SE, Wong WK, Mockler TC: Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Res. 2010, 20: 45-58. 10.1101/gr.093302.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D'Angelo C, Bornberg-Bauer E, Kudla J, Harter K: The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007, 50: 347-363. 10.1111/j.1365-313X.2007.03052.x.PubMedView ArticleGoogle Scholar
- Gautier L, Cope L, Bolstad BM, Irizarry RA: affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.PubMedView ArticleGoogle Scholar
- 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: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-297. 10.1093/nar/gkl031.PubMedPubMed CentralView ArticleGoogle Scholar
- 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: 2498-2504. 10.1101/gr.1239303.PubMedPubMed CentralView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034.PubMedPubMed CentralView ArticleGoogle Scholar
- Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR: Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005, 139: 5-17. 10.1104/pp.105.063743.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMedPubMed CentralView ArticleGoogle Scholar