Transcriptome and Metabolome Analyses Provide Insights into the Stomium Degeneration Mechanism in Lily

Lily (Lilium spp.) is a widely cultivated horticultural crop that has high ornamental and commercial value but also the serious problem of pollen pollution. However, mechanisms of anther dehiscence in lily remain largely unknown. In this study, the morphological characteristics of the stomium zone (SZ) from different developmental stages of ‘Siberia’ lily anthers were investigated. In addition, transcriptomic and metabolomic data were analyzed to identify the differentially expressed genes (DEGs) and secondary metabolites involved in stomium degeneration. According to morphological observations, SZ lysis occurred when flower buds were 6–8 cm in length and was completed in 9 cm. Transcriptomic analysis identified the genes involved in SZ degeneration, including those associated with hormone signal transduction, cell structure, reactive oxygen species (ROS), and transcription factors. A weighted co-expression network showed strong correlations between transcription factors. In addition, TUNEL (TdT-mediated dUTP nick-end labeling) assays showed that programmed cell death was important during anther SZ degeneration. Jasmonates might also have key roles in anther dehiscence by affecting the expression of the genes involved in pectin lysis, water transport, and cysteine protease. Collectively, the results of this study improve our understanding of anther dehiscence in lily and provide a data platform from which the molecular mechanisms of SZ degeneration can be revealed.