LAUSR

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.