Abstract
BACKGROUND: Lilium lancifolium, a very important cold-resistant wild flower for lily cold resistance breeding, is widely distributed in southwestern and northeastern China. To gain a better understanding of the cold signaling pathway and the molecular metabolic reactions involved in the cold response, we performed a genome-wide transcriptional analysis using RNA-Seq. RESULTS: Approximately 104,703 million clean 90- bp paired-end reads were obtained from three libraries (CK 0 h, Cold-treated 2 h and 16 h at 4 °C); 18,736 unigenes showed similarity to known proteins in the Swiss-Prot protein database, and 15,898, 13,705 and 1849 unigenes aligned to existing sequences in the KEGG and COG databases (comprising 25 COG categories) and formed 12 SOM clusters, respectively. Based on qRT-PCR results, we studied three signal regulation pathways--the Ca(2+) and ABA independent/dependent pathways--that conduct cold signals to signal transduction genes such as LlICE and LlCDPK and transcription factor genes such as LlDREB1/CBF, LlAP2/EREBP, LlNAC1, LlR2R3-MYB and LlBZIP, which were expressed highly in bulb. LlFAD3, Llβ-amylase, LlP5CS and LlCLS responded to cold and enhanced adaptation processes that involve changes in the expression of transcripts related to cellular osmoprotectants and carbohydrate metabolism during cold stress. CONCLUSIONS: Our study of differentially expressed genes involved in cold-related metabolic pathways and transcription factors facilitated the discovery of cold-resistance genes and the cold signal transcriptional networks, and identified potential key components in the regulation of the cold response in L lancifolium, which will be most beneficial for further research and in-depth exploration of cold-resistance breeding candidate genes in lily.