Metabolomics integrated with transcriptomics provides insights into the phenylpropanoids biosynthesis pathway in Lilium davidii var. unicolor and L. lancifolium Thunb.

Lilium spp. is a world-famous bulbous flower with outstanding ornamental, edible, and medicinal value. Evaluating the taste of edible lilies and identifying important active substances and genes are necessary for germplasm improvement, new variety breeding, and industrial application. To better understand the phenylpropanoids and regalosides biosynthesis, L. davidii var. unicolor and L. lancifolium Thunb. bulbs were used for transcriptome and metabolite analysis. Results showed that the phenols and flavonoid contents in JT were lower than in LT, while the saponins and alkaloid contents in JT were higher than in LT. A total of 20,520 differentially expressed genes (DEGs) and 383 differential metabolites were searched. Integrated transcriptomics and metabolomics analysis showed that phenylpropanoid biosynthesis and flavonoid biosynthesis were differentially altered. Ninety-nine unigenes encoding ten phenolic acids and two flavonoids were identified as candidate genes involved in phenylpropanoid and flavonoid biosynthesis. WGCNA analysis showed 76 phenylpropanoid and flavonoid biosynthesis-related unigenes were verified as likely to be involved in phenylpropanoid metabolism and regalosides accumulation. Among them, 15 genes were used for qRT-PCR, and four genes were utilized for tissue-specific expression pattern analysis. Down-regulation of LdPAL2 and LdC4H1 in bulbs of L. davidii var. unicolor via virus induced gene silence (VIGS) reduced the contents of p-coumaric acid and cinnamic acid. These results contribute to understanding phenylpropanoid metabolism and identifying potential functional genes for improving the regalosides and flavonoids content in Lilium bulbs.

Extraction, purification, characteristics, bioactivities, prospects, and toxicity of Lilium spp. polysaccharides.

The genus Lilium has been widely used worldwide as a food and medicinal ingredient in East Asia for over 2000 years due to its higher nutritional and medicinal value. Polysaccharide is the most important bioactive ingredient in Lilium spp. and has various health benefits. Recently, Lilium spp. polysaccharides (LSPs) have attracted significant attention from industries and researchers due to their various biological properties, such as antioxidant, immunomodulatory, antitumor, antibacterial, hypoglycaemic, and anti-radiation. However, the development and utilization of LSP-based functional biomaterials and medicines are limited by a lack of comprehensive understanding regarding the structure-activity relationships (SARs), industrial applications, and safety of LSPs. This review provides an inclusive overview of the extraction, purification, structural features, bioactivities, and mechanisms of LSPs. SARs, applications, toxicities, and influences of structural modifications on bioactivities are also highlighted, and the potential development and future study direction are scrutinized. This article aims to offer a complete understanding of LSPs and provide a foundation for further research and application of LSPs as therapeutic agents and multifunctional biomaterials.

Transcriptomic analysis reveals candidate genes associated with anther development in Lilium Oriental Hybrid 'Siberia'.

Lily (Lilium spp. and hybrids) is an important cut flower crop worldwide. Lily flowers have large anthers, which release a large amount of pollen that stains the tepals or clothing and thus can affect the commercial value of cut flowers. In this study, lily Oriental 'Siberia' was used to investigate the regulatory mechanism of lily anther development, which may provide information to prevent pollen pollution in the future. Based on the flower bud length, anther length and color, and anatomical observations, lily anther development was categorized into five stages: green (G), green-to-yellow 1 (GY1), green-to-yellow 2 (GY2), yellow (Y), and purple (P). Total RNA was extracted from the anthers at each stage for transcriptomic analysis. A total of 268.92-Gb clean reads were generated, and 81,287 unigenes were assembled and annotated. The number of differentially expressed genes (DEGs) and unique genes were largest for the pairwise comparison between the G and GY1 stages. The G and P samples were clustered separately, whereas the GY1, GY2, and Y samples were clustered together in scatter plots from a principal component analysis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses of DEGs detected in the GY1, GY2, and Y stages revealed that the pectin catabolic process, hormone levels, and phenylpropanoid biosynthesis were enriched. The DEGs associated with jasmonic acid biosynthesis and signaling were highly expressed at the early stages (G and GY1), whereas the DEGs associated with phenylpropanoid biosynthesis were mainly expressed in the intermediate stages (GY1, GY2, and Y). The DEGs involved in the pectin catabolic process were expressed at advanced stages (Y and P). Cucumber mosaic virus-induced gene silencing of LoMYB21 and LoAMS caused a strongly inhibited anther dehiscence phenotype, but without affecting the development of other floral organs. These results provide novel insights for understanding the regulatory mechanism of anther development in lily and other plants.

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.

SMRT and Illumina RNA Sequencing and Characterization of a Key NAC Gene LoNAC29 during the Flower Senescence in Lilium oriental 'Siberia'.

Lily (Lilium spp.) is an important cut flower around the world. Flower senescence in lilies is characterized by the wilting and abscission of tepals, which results in a decrease in flower quality and huge economic loss. However, the mechanism underlying flower senescence in lilies is largely unknown. In this study, single-molecule, real-time (SMRT) and Illumina sequencing were carried out in L. oriental 'Siberia'. Sequencing yielded 73,218 non-redundant transcripts, with an N50 of 3792 bp. These data were further integrated with three published transcriptomes through cogent analysis, which yielded 62,960 transcripts, with an increase in N50 of 3935 bp. Analysis of differentially expressed genes showed that 319 transcription factors were highly upregulated during flower senescence. The expression of twelve NAC genes and eleven senescence-associated genes (SAGs) showed that LoNAC29 and LoSAG39 were highly expressed in senescent flowers. Transient overexpression of LoNAC29 and LoSAG39 in tepals of lily notably accelerated flower senescence, and the promoter activity of LoSAG39 was strongly induced by LoNAC29. This work supported new evidence for the molecular mechanism of flower senescence and provided better sequence data for further study in lilies.

LhGST is an anthocyanin-related glutathione S-transferase gene in Asiatic hybrid lilies (Lilium spp.).

KEY MESSAGE: LhGST, an anthocyanin-related GST gene, was identified from Asiatic hybrid lilies. Expression and functional analyses demonstrated that LhGST might be involved in anthocyanin sequestration in lily tepals. Anthocyanins are responsible for the pink, red and purple pigmentation of flowers in Asiatic hybrid lilies, synthesized at the cytoplasmic surface of the endoplasmic reticulum (ER) and then transported to the vacuole. To date, the mechanism involved in the intracellular transport of anthocyanins in lilies has not been well elucidated. Here, full-length glutathione S-transferase gene (LhGST) was identified from lilies. Expression analysis revealed that LhGST was positively correlated with anthocyanin accumulation. Phylogenetic tree analysis showed that LhGST clustered with other anthocyanin-related GSTs in the same phi clade. Moreover, functional complementation of an Arabidopsis tt19 mutant demonstrated that LhGST might be involved in anthocyanin accumulation in lily tepals. Additionally, according to phenotype analysis, LhGST was found to be correlated with the transport of anthocyanin in lilies by virus-induced gene silencing (VIGS) assay. In addition, cis-element analysis of the LhGST promoter showed the presence of ABA-, auxin-, MeJA-, gibberellin-, light-, and stress-responsive elements and an MYB recognition site (MRS, CCGTTG). Yeast one-hybrid and dual-luciferase report assays revealed that the promoter of LhGST was activated by LhMYB12-lat, which is a key R2R3-MYB transcription factor that regulates anthocyanin biosynthesis in lilies. In conclusion, our results revealed that LhGST plays a key role in anthocyanin transport and accumulation in the tepals of lilies.

Evaluation of resistance to Botrytis elliptica in Lilium hybrid cultivars.

Gray mold disease, caused by the fungus Botrytis elliptica, is one of the major diseases affecting Lilium species, and it has become a limiting factor in the production of ornamental lilies. To support selecting and breeding Botrytis-resistant cultivars, a total of 50 Lilium cultivars belonging to seven hybrid types were evaluated using a detached leaf technique for resistance to B. elliptica. Through resistance evaluations, Oriental × Trumpet and Oriental hybrid cultivars were classified as resistant lines, while Asiatic and Trumpet hybrids were classified as susceptible lines. A highly resistant (HR) Oriental hybrid, 'Sorbonne', and a highly susceptible (HS) Asiatic hybrid, 'Tresor', were selected for further study of early host-parasite interactions. After infection, B. elliptica grew faster and more easily on the leaves of 'Tresor' than on those of 'Sorbonne' during initial infection; when 'Tresor' leaves were completely infected, only a few lesions were observed on 'Sorbonne' leaves. Biochemical differences between these two cultivars were found following inoculation with B. elliptica, as shown by studies of reactive oxygen species (ROS) and the enzymatic antioxidant system. Rapid accumulation of H2O2 and ·O2- to trigger a defense response was detected in HR 'Sorbonne'. Meanwhile, higher levels of antioxidant activity, including SOD, POD and CAT, were found in HR 'Sorbonne' than in HS 'Tresor' before 48 h post-inoculation, but antioxidant activity was reduced with subsequent infection progress. These large and timely increases in ROS and antioxidant activities could be the main contributors to the high resistance of the 'Sorbonne' cultivar.

Characterization of LhSorTGA2, a novel TGA2-like protein that interacts with LhSorNPR1 in oriental hybrid lily Sorbonne.

BACKGROUND: Non-expressor of pathogenesis-related genes 1 (NPR1) regulates expression of pathogenesis-related (PR) genes by interacting with TGA family proteins during systemic acquired resistance (SAR). However, no TGA-like proteins or their interacting partners have been characterized in lily. RESULTS: In the present study, LhSorTGA2, a novel TGA-like protein, was identified as an interacting partner of LhSorNPR1 (an NPR-like protein) by bimolecular fluorescence complementation (BIFC) and yeast two-hybrid assay (Y2H). Subcellular localization of GFP-tagged proteins targeted LhSorTGA2 to the nucleus, whereas GFP-labeled LhSorNPR1 was observed both in the nucleus and at the cytomembrane. Sequence alignment revealed that LhSorTGA2 was featured with a basic leucine zipper (bZIP) domain and two glutamine rich acid domains (QI and QII). Further phylogenetic analysis showed that TGA family proteins can be grouped into three subclades, within which LhSorTGA2 was clustered into subclade I, together with AtTGA2/5/6. Expression of LhSorTGA2 was investigated in different tissues by qPCR, and the highest expression level was observed in stem. Besides, when treated with phytohormones (SA, MeJA, ETH and ABA) or fungal pathogen Botrytis elliptica, LhSorTGA2 expression was also induced at different time points post treatments. CONCLUSIONS: Collectively, these results suggested that LhSorTGA2 was an interacting partner of LhSorNPR1, which might function in regulating expression of PR genes in lily during SAR.

Validation of Reference Genes for Quantitative Real-Time PCR during Bicolor Tepal Development in Asiatic Hybrid Lilies (Lilium spp.).

Quantitative real-time PCR (qRT-PCR) is a reliable and high-throughput technique for gene expression studies, but its accuracy depends on the expression stability of reference genes. To date, several reliable reference gene identifications have been reported in Lilium spp., but none has been obtained for lily tepals at different developmental stages. In this study, ten candidate reference genes were selected and evaluated for their expression stability in Lilium 'Tiny Padhye' during the process of bicolor tepal development. The expression stability of these candidates was evaluated by three software programs (geNorm, NormFinder, and BestKeeper) and the comparative ΔCt method, and comprehensive stability rankings were generated by RefFinder. As a result, TIP41-like family gene (TIP41) and actin (ACT) were the best combination of reference genes for tepals at different developmental stages; TIP41 and F-box family gene (F-box) for tepals under shading treatment; ACT, actin11 (ACT11), and elongation factor 1-α (EF1-α) for different tissues; and ACT, TIP41, and ACT11 for all samples. The selected optimal reference genes were further verified by analyzing the expression levels of flavonoid 3'-hydroxylase (LhF3'H) and anthocyanidin 3-O-glucosyltransfersae (LhUFGT) in tepals at different developmental stages. This study provides useful information for gene expression characterization in lilies under different experimental conditions, and can serve as a basis for similar research in other closely related species.

Characterization of LhSorP5CS, a gene catalyzing proline synthesis in Oriental hybrid lily Sorbonne: molecular modelling and expression analysis.

BACKGROUND: Abiotic stresses negatively affect plant growth and flower production. In plants, P5CS proteins are key enzymes that catalyzed the rate-limiting steps of proline synthesis, and proline is a well-known osmoprotectant that is closely related to abiotic stress tolerance. However, information about the P5CS genes, their effects on proline accumulation, and their role in abiotic stress tolerance in Lilium is still lacking. RESULTS: We isolated and characterized a novel gene (LhSorP5CS) from Oriental hybrid lily cultivar Sorbonne. Phylogenetic analysis indicated that LhSorP5CS is a member of the P5CS family. The three-dimensional structure of LhSorP5CS predicted by homology modeling showed high similarity to its correspondant human P5CS template. Further gene expression analysis revealed that LhSorP5CS expression was up-regulated by NaCl, mannitol, and ABA, and that stress-exposed plants accumulated proline at a significantly higher level than in the control. CONCLUSIONS: LhSorP5CS characterized in this study is involved in proline synthesis in lily, and that it might play an important role in abiotic stress tolerance. However, there should be other P5CS homologues in the lily genome, and some of them could be highly stress-induced and more important for proline accumulation. Future studies on P5CS family genes would be of great importance to proline-related stress tolerance in lily.

Spatiotemporal Transcriptome Analysis Provides Insights into Bicolor Tepal Development in Lilium "Tiny Padhye".

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.

Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies.

Lily tepals have a short lifespan. Once the tepals senesce, the ornamental value of the flower is lost. Some cultivars have attractive purple ovaries and fruits which greatly enhance the ornamental value of Asiatic hybrid lilies. However, little is known about the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries. To investigate the transcriptional network that governs purple ovary coloration in Asiatic hybrid lilies, we obtained transcriptome data from green ovaries (S1) and purple ovaries (S2) of Asiatic "Tiny Padhye". Comparative transcriptome analysis revealed 4228 differentially expressed genes. Differential expression analysis revealed that ten unigenes including four CHS genes, one CHI gene, one F3H gene, one F3'H gene, one DFR gene, one UFGT gene, and one 3RT gene were significantly up-regulated in purple ovaries. One MYB gene, LhMYB12-Lat, was identified as a key transcription factor determining the distribution of anthocyanins in Asiatic hybrid lily ovaries. Further qPCR results showed unigenes related to anthocyanin biosynthesis were highly expressed in purple ovaries of three purple-ovaried Asiatic hybrid lilies at stages 2 and 3, while they showed an extremely low level of expression in ovaries of three green-ovaried Asiatic hybrid lilies during all developmental stages. In addition, shading treatment significantly decreased pigment accumulation by suppressing the expression of several unigenes related to anthocyanin biosynthesis in ovaries of Asiatic "Tiny Padhye". Lastly, a total of 15,048 Simple Sequence Repeats (SSRs) were identified in 13,710 sequences, and primer pairs for SSRs were designed. The results could further our understanding of the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries.

Rapid and sensitive detection of Lily symptomless virus by reverse transcription loop-mediated isothermal amplification.

Lily symptomless virus (LSV) is one of the most prevalent viruses that infect lily plants worldwide. A rapid and sensitive reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay was developed for detection of LSV, using two primer pairs that specifically amplified the conserved sequence of LSV coat protein. The optimum reaction conditions were as follows: 4mM MgCl2 and 0.8M betaine with incubation at 64°C for 30min. The limit of detection of LSV from infected lily leaves was 10-fold higher for RT-LAMP than for conventional RT-PCR. Moreover, RT-LAMP detected LSV in not only symptomatic, but also in symptomless tissues of infected plants. These findings indicate that our RT-LAMP method for LSV can serve as a low-cost, simple, and rapid alternative to conventional detection assays.