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. 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 Ca2+ 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. 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.

Our previous studies have indicated that a partial NAC domain protein gene is strongly up-regulated by cold stress (4 °C) in tiger lily (Lilium lancifolium). In this study, we cloned the full-length of this NAC gene, LlNAC2, to further investigate the function of LlNAC2 in response to various abiotic stresses and the possible involvement in stress tolerance of the tiger lily plant. LlNAC2 was noticeably induced by cold, drought, salt stresses, and abscisic acid (ABA) treatment. Promoter analysis showed that various stress-related cis-acting regulatory elements were located in the promoter of LlNAC2; and the promoter was sufficient to enhance activity of GUS protein under cold, salt stresses and ABA treatment. DREB1 (dehydration-responsive binding protein1) from tiger lily (LlDREB1) was proved to be able to bind to the promoter of LlNAC2 by yeast one-hybrid (Y1H) assay. LlNAC2 was shown to physically interact with LlDREB1 and zinc finger-homeodomain ZFHD4 from the tiger lily (LlZFHD4) by bimolecular fluorescence complementation (BiFC) assay. Overexpressing LlNAC2 in Arabidopsis thaliana showed ABA hypersensitivity and enhanced tolerance to cold, drought, and salt stresses. These findings indicated LlNAC2 is involved in both DREB/CBF-COR and ABA signaling pathways to regulate stress tolerance of the tiger lily.

Most commercial cultivars of lily are sensitive to abiotic stresses. However, tiger lily (Lilium lancifolium L.), one of the most widely distributed wild lilies in Asia, has strong abiotic stresses resistance. Thus, it is indispensable to identify stress-responsive candidate genes in tiger lily for the stress resistance improvement of plants. In this study, a MYB related homolog (LlMYB3) from tiger lily was functionally characterized as a positive regulator in plant stress tolerance. LlMYB3 is a nuclear protein with transcriptional activation activity at C-terminus. The expression of LlMYB3 gene was induced by multiple stress treatments. Several stress-related cis-acting regulatory elements (MYBRS, MYCRS, LTRE and DRE/CRT) were located within the promoter of LlMYB3; however, the promoter activity was not induced sufficiently by various stresses treatments. Overexpressing LlMYB3 in Arabidopsis thaliana L. transgenic plants showed ABA hypersensitivity and enhanced tolerance to cold, drought, and salt stresses. Furthermore, we found LlMYB3 highly co-expressed with LlCHS2 gene under cold treatment; yeast one-hybrid (Y1H) assays demonstrated LlMYB3 was able to bind to the promoter of LlCHS2. These findings suggest that the stress-responsive LlMYB3 may be involved in anthocyanin biosynthesis pathway to regulate stress tolerance of tiger lily.

RcNAC72, a key transcription factor that may respond to drought stress in Rosa chinensis ‘Old Blush’, was selected in our previous study. In the present study, we found that RcNAC72 is localized in the nucleus and is a transcriptional activator. RcNAC72 expression could be significantly induced by drought, low temperature, salt as well as abscisic acid (ABA) treatment. Analysis of the promoter revealed that multiple abiotic stress and hormone response elements were located in the promoter region. The promoter could respond to drought, low temperature, salt and ABA treatments to activate GUS gene expression. Overexpressing RcNAC72 in Arabidopsis thaliana enhanced sensitivity to ABA and tolerance to drought stress. Silencing of RcNAC72 by virus-induced gene silencing (VIGS) in rose leaves significantly reduced leaf water loss tolerance and leaf extension capacity. Physical interaction of RcNAC72 with RcDREB2A was shown by means of the yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays. RcABF4 was demonstrated to be able to bind to the promoter of RcNAC72 by means of the yeast one-hybrid (Y1H) assay. These results provide new insights into the regulatory network of RcNAC72 response to drought stress in roses.

The quality of Lily cut flower was determined by the quality of bulbs. During the process of vernalization and flower bud differentiation, sugar massively accumulated in the bulb, which influenced the bulb development. However, the details of sugar genes’ regulation mechanism for these processes were not fully understood. Here, morphological physiology, transcriptomes and gene engineering technology were used to explore this physiological change. Seventy-two genes of 25 kinds of sugar metabolism-related genes were annotated after re-analyzing transcriptome data of Oriental hybrid lily ‘Sorbonne’ bulbs, which were generated on Hiseq Illumina 2000. The results showed that these genes were closely related to lily bulb vernalization and development. Combining gene expression pattern with gene co-expression network, five genes (Contig5669, Contig13319, Contig7715, Contig1420 and Contig87292) were considered to be the most potential signals, and the sucrose transporter gene (SUT) was the focus of this study. Carbohydrate transport pathway and genes’ regulation mechanism were inferred through a physiological and molecular test. SUT seemed to be the sugar sensor that could sense and regulate sugar concentration, which might have effects on other genes, such as FT, LFY and so on. LoSUT2 and LoSUT4 genes were cloned from Oriental hybrid lily ‘Sorbonne’ by RACE, which was the first time for these genes in Oriental hybrid lily ‘Sorbonne’. The physiological properties of these proteins were analyzed such as hydrophobicity and phosphorylation. In addition, secondary and tertiary structures of proteins were predicted, which indicated the two proteins were membrane proteins. Their cellular locations were verified through positioning the experiment of the fluorescent vector. They were highly expressed in cells around phloem, which illustrated the key role of these genes in sugar transport. Furthermore, transient expression assays showed that overexpressed LoSUT2 and LoSUT4 in Arabidopsis thaliana bloomed significantly earlier than the wild type and the expression of FT, SOC1 and LFY were also affected by LoSUT2 and LoSUT4, which indicated that LoSUT2 and LoSUT4 may regulate plants flowering time.

Sugar transport and distribution plays an important role in lily bulb development and resistance to abiotic stresses. In this study, a member of the Sugar Will Eventually be Exported Transporters (SWEET) gene family, LoSWEET14, from Oriental hybrid lily ‘Sorbonne’ was identified. LoSWEET14 encodes a protein of 278 amino acids and is capable of transporting sucrose and some types of hexoses. The transcript level of the LoSWEET14 gene was significantly increased under various stress conditions including drought, cold, salt stresses, and abscisic acid (ABA) treatment. Overexpression of LoSWEET14 in tobacco (Nicotiana tabacum) showed that the transgenic lines had larger leaves, accumulated more soluble sugars, and were more resistant to drought, cold, and salt stresses, while becoming more sensitive to ABA compared with wild-type lines. Promoter analysis revealed that multiple stress-related cis-acting elements were found in the promoter of LoSWEET14. According to the distribution of cis-acting elements, different lengths of 5′-deletion fragments were constructed and the LoSWEET14-pro3(-540 bp) was found to be able to drive GUS gene expression in response to abiotic stresses and ABA treatment. Furthermore, a yeast one hybrid (Y1H) assay proved that the AREB/ABF (ABRE-binding protein/ABRE-binding factor) from lilies (LoABF2) could bind to the promoter of LoSWEET14. These findings indicated that LoSWEET14 is induced by LoABF2 to participate in the ABA signaling pathway to promote soluble sugar accumulation in response to multiple abiotic stresses.