The genus Lilium contains more than 100 wild species and numerous hybrid varieties. Some species of them have been used as medicine and food since ancient times. However, the research on the active components and the medical properties of lilies has only focused on a few species. In this study, the total phenolic acid content (TPC), total flavonoid content (TFC), and antioxidant capacity of 22 representative lilies were systematically investigated. The results showed that the TPC, TFC and antioxidant activity were highly variable among different lilies, but they were significantly positively correlated. Hierarchical cluster analysis indicated that L. henryi and L. regale were arranged in one group characterized by the highest TPC, TFC and antioxidant capacity, followed by Oriental hybrids and Trumpet and Oriental hybrids. The traditional edible and medicinal lilies were clustered in low TPC, TFC and antioxidant capacity group. A total of 577 secondary metabolites, including 201 flavonoids, 153 phenolic acids, were identified in the five species with great differences in antioxidant capacity by extensive targeted metabonomics. Differentially accumulated metabolites (DAMs) analysis reviewed that the DAMs were mainly enriched in secondary metabolic pathways such as isoflavonoid, folate, flavonoid, flavone, flavonol, phenylpropanoid, isoquinoline alkaloid biosynthesis, nicotinate and nicotinamide metabolism and so on. Correlation analysis identified that 64 metabolites were significantly positively correlated with antioxidant capacity (r ≥ 0.9 and p < 0.0001). These results suggested that the genus Lilium has great biodiversity in bioactive components. The data obtained greatly expand our knowledge of the bioactive constituents of Lilium spp. Additionally, it also highlights the potential application of Lilium plants as antioxidants, functional ingredients, cosmetic products and nutraceuticals.

Aerial bulbils are an important propagative organ, playing an important role in population expansion. However, the detailed gene regulatory patterns and molecular mechanism underlying bulbil formation remain unclear. Triploid Lilium lancifolium, which develops many aerial bulbils on the leaf axils of middle-upper stem, is a useful species for investigating bulbil formation. To investigate the mechanism of bulbil formation in triploid L. lancifolium, we performed histological and transcriptomic analyses using samples of leaf axils located in the upper and lower stem of triploid L. lancifolium during bulbil formation. Histological results indicated that the bulbils of triploid L. lancifolium are derived from axillary meristems that initiate de novo from cells at the base of petiole. Transcriptomic analysis generated approximately 650 million high-quality reads and 11,871 differentially expressed genes (DEGs). Functional analysis showed that the DEGs were significantly enriched in starch and sucrose metabolism and plant hormone signal transduction. Starch synthesis and accumulation likely promoted the initiation of upper bulbils in triploid L. lancifolium. Hormone-associated pathways exhibited distinct patterns of change in each sample. Auxin likely promoted the initiation of bulbils and then inhibited further bulbil formation. High biosynthesis and low degradation of cytokinin might have led to bulbil formation in the upper leaf axil. The present study achieved a global transcriptomic analysis focused on gene expression changes and pathways’ enrichment during upper bulbil formation in triploid L. lancifolium, laying a solid foundation for future molecular studies on bulbil formation.

The Asiatic hybrid lily (Lilium spp.) is a horticultural crop with high commercial value and diverse anthocyanin pigmentation patterns. However, the regulatory mechanism underlying lily flower color has been largely unexplored. Here, we identified a WRKY transcription factor from lily tepals, LhWRKY44, whose expression was closely associated with anthocyanin accumulation. Functional verification indicated that LhWRKY44 positively regulated anthocyanin accumulation. LhWRKY44 physically interacted with LhMYBSPLATTER and directly bound to the LhMYBSPLATTER promoter, which enhanced the effect of the LhMYBSPLATTER-LhbHLH2 MBW complex activator on anthocyanin accumulation. Moreover, EMSA and dual-luciferase assays revealed that LhWRKY44 activated and bound to the promoters of gene LhF3H and the intracellular anthocyanin-related glutathione S-transferase gene LhGST. Interestingly, our further results showed that LhWRKY44 participated in light and drought-induced anthocyanin accumulation, and improved the drought tolerance in lily via activating stress-related genes. These results generated a multifaceted regulatory mechanism for the LhWRKY44-meditaed enhancement by the environmental signal pathway of anthocyanin accumulation and expanded our understanding of the WRKY-mediated transcriptional regulatory hierarchy modulating anthocyanin accumulation in Asiatic hybrid lilies.

The bicolor Asiatic hybrid lily cultivar ‘Tiny Padhye’ is an attractive variety because of its unique color pattern. During its bicolor tepal development, the upper tepals undergo a rapid color change from green to white, while the tepal bases change from green to purple. However, the molecular mechanisms underlying these changes remain largely uncharacterized. To systematically investigate the dynamics of the lily bicolor tepal transcriptome during development, we generated 15 RNA-seq libraries from the upper tepals (S2-U) and basal tepals (S1-D, S2-D, S3-D, and S4-D) of Lilium ‘Tiny Padhye’. Utilizing the Illumina platform, a total of 295,787 unigenes were obtained from 713.12 million high-quality paired-end reads. A total of 16,182 unigenes were identified as differentially expressed genes during tepal development. Using Kyoto Encyclopedia of Genes and Genomes pathway analysis, candidate genes involved in the anthocyanin biosynthetic pathway (61 unigenes) and chlorophyll metabolic pathway (106 unigenes) were identified. Further analyses showed that most anthocyanin biosynthesis genes were transcribed coordinately in the tepal bases, but not in the upper tepals, suggesting that the bicolor trait of ‘Tiny Padhye’ tepals is caused by the transcriptional regulation of anthocyanin biosynthetic genes. Meanwhile, the high expression level of chlorophyll degradation genes and low expression level of chlorophyll biosynthetic genes resulted in the absence of chlorophylls from ‘Tiny Padhye’ tepals after flowering. Transcription factors putatively involved in the anthocyanin biosynthetic pathway and chlorophyll metabolism in lilies were identified using a weighted gene co-expression network analysis and their possible roles in lily bicolor tepal development were discussed. In conclusion, these extensive transcriptome data provide a platform for elucidating the molecular mechanisms of bicolor tepals in lilies and provide a basis for similar research in other closely related species.