Abstract:Abstract Background High temperature seriously limits the annual production of fresh cut lilies, which is one of the four major cut flowers in the global cut flower market. There were few transcriptomes focused on the gene expression of lilies under heat stress. In order to reveal the potential heat response patterns in bulbous plants and provide important genes for further genetic engineering techniques to improve thermotolerance of lily, RNA sequencing of lilies under heat treatments were conducted. Results In this study, seedlings of Lilium longiflorum ‘White Heaven’ were heat-treated at 37 °C for different lengths of time (0 h, 0.5 h, 1 h, 3 h, 6 h, and 12 h with a 12 h-light/12 h-dark cycle). The leaves of these lily seedlings were immediately collected after heat treatments and quickly put into liquid nitrogen for RNA sequencing. 109,364,486–171,487,430 clean reads and 55,044 unigenes including 21,608 differentially expressed genes (DEGs) (fold change ≥2) were obtained after heat treatment. The number of DEGs increased sharply during the heat treatments of 0.5 h–1 h and 1 h–3 h compared to that of other periods. Genes of the heat stress transcription factor (HSF) family and the small heat shock proteins (small HSPs, also known as HSP20) family responded to heat stress early and quickly. Compared to that of the calcium signal and hormone pathways, DEGs of the HSF-HSP pathway and reactive oxygen species (ROS) pathway were significantly and highly induced. Moreover, they had the similar expression pattern in response to heat stress. Small HSPs family genes were the major components in the 50 most highly induced genes at each heat stress treatment and involved in ROS pathway in the rapid response to heat stress. Furthermore, the barley stripe mosaic virus induced gene silencing (BSMV-VIGS) of LlHsfA2 caused a significantly reduced thermotolerance phenotype in Lilium longiflorum ‘White Heaven’, meanwhile decreasing the expression of small HSPs family genes and increasing the ROS scavenging enzyme ascorbate peroxidase (APX) genes, indicating the potential interplay between these two pathways. Conclusions Based on our transcriptomic analysis, we provide a new finding that small HSPs play important roles in crosstalk between HSF-HSP and ROS pathways in heat stress response of lily, which also supply the groundwork for understanding the mechanism of heat stress in bulbous plants.

An AP2/ERF member LlERF012 confers thermotolerance via activation of HSF pathway in lily

Ethylene Response Factor LlERF110 Mediates Heat Stress Response via Regulation of LlHsfA3A Expression and Interaction with LlHsfA2 in Lilies (Lilium longiflorum)

LlbHLH87 interacts with LlSPT to modulate thermotolerance via activation of LlHSFA2 and LlEIN3 in lily

Lily LlHSFC2 coordinates with HSFAs to balance heat stress response and improve thermotolerance

Overexpression of a Novel Heat-Inducible Ethylene-Responsive Factor Gene LlERF110 from Lilium Longiflorum Decreases Thermotolerance

Lily HD-Zip I Transcription Factor LlHB16 Promotes Thermotolerance by Activating LlHSFA2 and LlMBF1c

The Heat Stress Transcription Factor LlHsfA4 Enhanced Basic Thermotolerance through Regulating ROS Metabolism in Lilies (Lilium Longiflorum)

Over-expression of LlHsfA2b , a lily heat shock transcription factor lacking trans-activation activity in yeast, can enhance tolerance to heat and oxidative stress in transgenic Arabidopsis seedlings

A LlMYB305-LlC3H18-LlWRKY33 Module Regulates Thermotolerance in Lily

LlHSFA1, a Novel Heat Stress Transcription Factor in Lily (lilium Longiflorum), Can Interact with LlHSFA2 and Enhance the Thermotolerance of Transgenic Arabidopsis Thaliana.

Lily WRKY factor LlWRKY22 promotes thermotolerance through autoactivation and activation of LlDREB2B

LlWRKY39 is involved in thermotolerance by activating LlMBF1c and interacting with LlCaM3 in lily (Lilium longiflorum)

HD-Zip I Protein LlHOX6 Antagonizes Homeobox Protein LlHB16 to Attenuate Basal Thermotolerance in Lily.

Alternative splicing provides a mechanism to regulate LlHSFA3 function in response to heat stress in lily.

A novel R2R3-MYB transcription factor LlMYB305 from <i>Lilium longiflorum</i> plays a positive role in thermotolerance via activating heat-protective genes

Cloning and Characterization of HsfA2 from Lily (lilium Longiflorum)

A Lily Heat-Inducible Multi-Protein Bridging Factor, LlMBF1c, Plays a Crucial Role in Plant Thermotolerance

A Lily Membrane-Associated NAC Transcription Factor, LlNAC014, is Involved in Thermotolerance Via Activation of the DREB2-HSFA3 Module

Involvement of Ca 2+ and CaM3 in Regulation of Thermotolerance in Lily ( Lilium longiflorum )

Small HSPs Play an Important Role in Crosstalk Between HSF-HSP and ROS Pathways in Heat Stress Response Through Transcriptomic Analysis in Lilies (lilium Longiflorum)

Comparative transcriptome analysis reveals heat stress-responsive genes and their signalling pathways in lilies (Lilium longiflorum vs. Lilium distichum)