A systematic regulatory network related to bulbil formation in Lilium lancifolium based on metabolome and transcriptome analyses.

Background: Lilium lancifolium is a special wild triploid species native to China and can produce abundant bulbils on its stem under natural conditions, which is very valuable to study bulbil organogenesis in plants. Although similar to the lateral and tillering principles, the molecular mechanism underlying bulbil formation has remained incompletely understood.

Results: The metabolome and transcriptome of L. lancifolium bulbils across four development stages were analyzed. The pairwise comparison of metabolomes across the four stages identified 17 differential hormones, predominantly auxin (IAA), cytokinin (CK), and jasmonic acid (JA). Short Time-series Expression Miner (STEM) trend analysis of differential genes revealed four significant trends across these stages. The KEGG enrichment analysis of the four clusters highlighted pathways, such as plant hormone signal transduction, which were speculated to play a crucial role in development stages. these pathways were speculated to play a crucial role in development stages. To explore the key differential expressed genes and transcription factors associated with bulbil occurrence, two periods were focused on: Ll_UN and Ll_DN, which represented the stages with and without bulbils, respectively. Through correlation analysis and qRT-PCR analysis, 11 candidate differentially expressed genes and 27 candidate transcription factors were selected. By spraying exogenous hormones to validate these candidates, LlbHLH128, LlTIFY10A, LlbHLH93, and LlMYB108, were identified as the key genes for L. lancifolium bulbils.

Conclusion: A regulatory network of L. lancifolium bulbil development was predicted. LlTIFY10A and LlbHLH93 might be involved in the JA and auxin signal transduction pathways, which jointly formed a regulatory network to affect the occurrence of L. lancifolium bulbil. This study not only provided more information about the differentially expressed genes and metabolites through transcriptome and metabolomics analyses, but also provided a clearer understanding of the effect of hormones on bulbil formation in lily.