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.

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.

The bulbil is an important vegetative reproductive organ in triploid tiger lily (Lilium lancifolium). Based on our previously obtained transcriptome data, we screened two WUSCHEL-related homeobox (WOX) genes closely related to bulbil formation, LlWOX9 and LlWOX11. However, the biological functions and regulatory mechanisms of LlWOX9 and LlWOX11 are unclear. In this study, we cloned the full-length coding sequences of LlWOX9 and LlWOX11. Transgenic Arabidopsis (Arabidopsis thaliana) showed increased branch numbers, and the overexpression of LlWOX9 and LlWOX11 in stem segments promoted bulbil formation, while the silencing of LlWOX9 and LlWOX11 inhibited bulbil formation, indicating that LlWOX9 and LlWOX11 are positive regulators of bulbil formation. Cytokinin type-B response regulators could bind to the promoters of LlWOX9 and LlWOX11 and promote their transcription. LlWOX11 could enhance cytokinin pathway signaling by inhibiting the transcription of type-A LlRR9. Our study enriches the understanding of the regulation of plant development by the WOX gene family and lays a foundation for further research on the molecular mechanism of bulbil formation in lily.

The application of cosmetics is indispensable in our current society. In recent years, with an increasing awareness of the long-term health benefits of naturally sourced ingredients, plant-based cosmetic products have gained increasing attention. Lilium belongs to the Liliaceae family, which is one of the main plant families used in cosmetics for skin care treatment. A large number of studies have shown that Lilium plants are rich in components such as phenolic acids, flavonoids, and polysaccharides, with high potential for cosmetic applications. However, the application of lilies in cosmetics has not been systematically reported. This knowledge gap can easily lead to the neglect of its application in cosmetics because lilies are most familiar as ornamental plants. Integrating academic papers and patent publications, we analyzed the potential cosmetic application ingredients in lily, as well as their applications in cosmetics and related efficacy. Patent analysis showed that applications for lily-related cosmetic patents are mainly concentrated in East Asia, including China, Korea, and Japan. The application of lilies involves all aspects of cosmetics, such as sunscreens, facial cleansers, facial masks, conditioners, and so on. Its functions are also rich and diverse, including antiaging, radiation protective, whitening, moisturizing, freckle removal, acne treatment, and hair regeneration promotion. In addition, lilies are compatible with the application of other herbs. Moreover, with a change in people’s consumption concepts and the consideration of long-term health benefits, lily-based food and medicine innovation with health care and beautification effects may be a promising direction.

The bulbil is an important vegetative reproductive organ in triploid Lilium lancifolium whose development is promoted by cytokinins. Type-B response regulators (RRs) are critical regulators that mediate primary cytokinin responses and promote cytokinin-induced gene expression. However, the function of cytokinin type-B Arabidopsis RRs (ARRs) in regulating bulbil formation is unclear. In this study, we identified five type-B LlRRs, LlRR1, LlRR2, LlRR10, LlRR11 and LlRR12, in L. lancifolium for the first time. The five LlRRs encode proteins of 715, 675, 573, 582 and 647 amino acids. All of the regulators belong to the B-I subfamily, whose members typically contain a conserved CheY-homologous receiver (REC) domain and an Myb DNA-binding (MYB) domain at the N-terminus. As transcription factors, all five type-B LlRRs localize at the nucleus and are widely expressed in plant tissues, especially during axillary meristem (AM) formation. Functional analysis showed that type-B LlRRs are involved in bulbil formation in a functionally redundant manner and can activate LlRR9 expression. In summary, our study elucidates the process by which cytokinins regulate bulbil initiation in L. lancifolium through type-B LlRRs and lays a foundation for research on the molecular mechanism of bulbil formation in the lily.

Previous studies have suggested that multidrug and toxic compound extrusion (MATE) proteins might be involved in flavonoid transportation. However, whether MATE proteins are involved in anthocyanin accumulation in Lilium is unclear. Here, a flavonoid transport-related MATE candidate gene, LhDTX35, was cloned from the Asiatic hybrid lily cultivar ‘Tiny Padhye’ by rapid amplification of 5’ and 3’ cDNA ends (RACE) and found to encode 507 amino acids. BLASTx results indicated that LhDTX35 showed high homology to the DTX35 genes of other species. Bioinformatics analysis predicted that the protein encoded by LhDTX35 possessed 12 typical transmembrane segments and had functional domains typical of the MATE-like superfamily. Phylogenetic analysis grouped LhDTX35 in the same clade as the DTX35 of other species. Notably, the expression pattern of LhDTX35 was positively correlated with floral anthocyanin accumulation in ‘Tiny Padhye’. A subcellular localization assay showed that the protein encoded by LhDTX35 was plasmalemma localized but not nuclear, indicating that the LhDTX35 gene may function as a carrier protein to transport anthocyanins in Lilium. Functional complementation of the Arabidopsis DTX35 gene demonstrated that LhDTX35 could restore silique-infertility and the anthocyaninless phenotype of an Arabidopsis DTX35 mutant. These results indicated that LhDTX35 might be involved in anthocyanin accumulation in Lilium.