Transcriptome Analysis Reveals Coexpression Networks and Hub Genes Involved in Papillae Development in Lilium auratum

Lilium is a genus of important ornamental plants with many colouring pattern variations. Lilium auratum is the parent of Oriental hybrid lilies. A typical feature of L. auratum is the presence of red-orange special raised spots named papillae on the interior tepals. Unlike the usual raised spots, the papillae are slightly rounded or connected into sheets and usually have hairy tips. To elucidate the potential genes regulating papillae development in L. auratum, we performed high-throughput sequencing of its tepals at different stages. Genes involved in the flavonoid biosynthesis pathway were significantly enriched during the colouration of the papillae, and CHS, F3H, F3'H, FLS, DFR, ANS, and UFGT were significantly upregulated. To identify the key genes involved in the papillae development of L. auratum, we performed weighted gene coexpression network analysis (WGCNA) and further analysed four modules. In total, 51, 24, 1, and 6 hub genes were identified in four WGCNA modules, MEbrown, MEyellow, MEpurple, and MEred, respectively. Then, the coexpression networks were constructed, and important genes involved in trichome development and coexpressed with anthocyanin biosynthesis genes, such as TT8, TTG1, and GEM, were identified. These results indicated that the papillae are essentially trichomes that accumulate anthocyanins. Finally, we randomly selected 12 hub genes for qRT-PCR analysis to verify the accuracy of our RNA-Seq analysis. Our results provide new insights into the papillae development in L. auratum flowers.