Identification of differentially expressed genes in flower, leaf and bulb scale of Lilium oriental hybrid 'Sorbonne' and putative control network for scent genes

The giant genome of lily provides insights into the hybridization of cultivated lilies

Lilies are economically important monocots known for their ornamental flowers, bulbs, and large genomes. The absence of their genomic information has impeded evolutionary studies and genome-based breeding efforts. Here, we present reference genomes for Lilium sargentiae (lily, 35.66 Gb) and Gloriosa superba (flame lily, 5.09 Gb). The giant lily genome is shaped by recent long terminal repeat retroelements. Phylogenetic analysis reveals diverse, independent origins of lily cultivars. Gene families involved in sucrose and starch metabolism are significantly expanded in the lily genome. Key homologs of XTH22, SOC1, and AP1/FUL-like genes regulate the development, bud growth transition, and floral bud growth transition of lily bulbs. Colchicine biosynthetic gene clusters are identified in G. superba but are absent in L. sargentiae, highlighting independent colchicine evolution in Colchicaceae. These genomic insights enhance understanding of Liliales evolution, providing a foundation for future breeding and molecular research.

Investigation of genes involved in scent and color production in Rosa damascena Mill

Rosa damascena Mill., commonly known as the King Flower, is a fragrant and important species of the Rosaceae family. It is widely used in the perfumery and pharmaceutical industries. The scent and color of the flowers are significant characteristics of this ornamental plant. This study aimed to investigate the relative expression of MYB1, CCD1, FLS, PAL, CER1, GT1, ANS and PAR genes under two growth stages (S1 and S2) in two morphs. The CCD1 gene pathway is highly correlated with the biosynthesis of volatile compounds. The results showed that the overexpression of MYB1, one of the important transcription factors in the production of fragrance and color, in the Hot pink morph of sample S2 increased the expression of PAR, PAL, FLS, RhGT1, CCD1, ANS, CER1, and GGPPS. The methyl jasmonate (MeJA) stimulant had a positive and cumulative effect on gene expression in most genes, such as FLS in ACC.26 of the S2 sample, RhGT1, MYB1, CCD1, PAR, ANS, CER1, and PAL in ACC.1. To further study, a comprehensive analysis was performed to evaluate the relationship between the principal volatile compounds and colors. Our data suggest that the rose with pink flowers had a higher accumulation content of flavonoids and anthocyanin. To separate essential oil compounds, GC/MS analysis identified 26 compounds in four samples. The highest amount of geraniol, one of the main components of damask rose, was found in the Hot pink flower, 23.54%, under the influence of the MeJA hormone.

Genome-Wide Characterization of Differentially Expressed Scent Genes in the MEP Control Network of the Flower of Lilium 'Sorbonne'

Fragrance is an important feature of ornamental lilies. Components of volatile substances and important genes for monoterpene synthesis in the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway were examined in this study. Twenty volatile compounds (2 in the budding stage, 3 in the initial flowering stage, 7 in the semi-flowering stage, 17 in the full-flowering stage, and 5 in withering stage) were detected in the Oriental lily 'Sorbonne' using gas chromatography-mass spectrometry. The semi- and full-flowering stages were key periods for volatile substance production and enzyme function. Sequence assembly from samples collected during all flowering stages resulted in the detection of 274,849 genes and 129,017 transcripts. RNA sequencing and heatmapping led to the detection of genes in the MEP monoterpene metabolism pathway. Through gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, we extracted key genes (LiDXS2, LiLIS, and LiMYS) and transcription factors (in the bHLH, MYB, HD-ZIP, and NAC families) associated with the MEP pathway. Tissue localization revealed that LiDXS2, LiLIS, and LiMYS were expressed in Lilium 'Sorbonne' petals in the full-flowering stage. Genes regulating the 1-deoxy-D-X-lignone-5-phosphate synthase family of rate-limiting enzymes, involved in the first step of monoterpene synthesis, showed high expression in the semi- and full-flowering stages. LiDXS2 was cloned and localized in chloroplast subcells. The relative expression of terpene-related genes in the MEP and mevalonic acid pathways of wild-type and LiLIS/LiMYS transgenic Arabidopsis thaliana, and changes in chemical composition, confirmed that LiLIS/LiMYS regulates the monoterpene synthesis pathway. The results of this study provide a theoretical basis for the synthesis of lily aromatic substances and the cultivation of new garden flower varieties.

A BAHD acyltransferase contributes to the biosynthesis of both ethyl benzoate and methyl benzoate in the flowers of Lilium oriental hybrid 'Siberia'

Lily is a popular flower worldwide due to its elegant appearance and pleasant fragrance. Floral volatiles of lily are predominated by monoterpenes and benzenoids. While a number of genes for monoterpene biosynthesis have been characterized, the molecular mechanism underlying floral benzenoid formation in lily remains unclear. Here, we report on the identification and characterization of a novel BAHD acyltransferase gene that contributes to the biosynthesis of two related floral scent benzoate esters, ethyl benzoate and methyl benzoate, in the scented Lilium oriental hybrid 'Siberia'. The emission of both methyl benzoate and ethyl benzoate in L. 'Siberia' was found to be tepal-specific, floral development-regulated and rhythmic. Through transcriptome profiling and bioinformatic analysis, a BAHD acyltransferase gene designated LoAAT1 was identified as the top candidate gene for the production of ethyl benzoate. In vitro enzyme assays and substrate feeding assays provide substantial evidence that LoAAT1 is responsible for the biosynthesis of ethyl benzoate. It was interesting to note that in in vitro enzyme assay, LoAAT1 can also catalyze the formation of methyl benzoate, which is typically formed by the action of benzoic acid methyltransferase (BAMT). The lack of an expressed putative BAMT gene in the flower transcriptome of L. 'Siberia', together with biochemical and expression evidence, led us to conclude that LoAAT1 is also responsible for, or at least contributes to, the biosynthesis of the floral scent compound methyl benzoate. This is the first report that a member of the plant BAHD acyltransferase family contributes to the production of both ethyl benzoate and methyl benzoate, presenting a new mechanism for the biosynthesis of benzoate esters.

Full-length transcriptome analysis revealed that 2,4-dichlorophenoxyacetic acid promoted in vitro bulblet initiation in lily by affecting carbohydrate metabolism and auxin signaling

Bulblet initiation, including adventitious bud initiation and bulblet formation, is a crucial process for lily and other bulbous flowers that are commercially propagated by vegetative means. Here, by a hybrid strategy combining Pacific Biosciences (PacBio) full-length sequencing and Illumina RNA sequencing (RNA-seq), high-quality transcripts of L. brownii (Lb) and its variety, L. brownii var. giganteum (Lbg), during in vitro bulblet initiation were obtained. A total of 53,576 and 65,050 high-quality non-redundant full-length transcripts of Lbg and Lb were generated, respectively. Morphological observation showed that Lbg possessed a stronger capacity to generate bulblets in vitro than Lb, and 1 mg L-1 2,4-dichlorophenoxyacetic acid (2,4-D) significantly increased bulblet regeneration rate in two lilies. Screening of differentially expressed transcripts (DETs) between different stages and Mfuzz analysis showed 0 DAT to 1 DAT was the crucial stage with the most complex transcriptional change, with carbohydrate metabolism pathway was significantly enriched. In addition, 6,218 and 8,965 DETs were screened between the 2,4-D-treated group and the control group in Lbg and Lb, respectively. 2,4-D application had evident effects on the expression of genes involved in auxin signaling pathway, such as TIRs, ARFs, Aux/IAAs, GH3s and SAURs. Then, we compared the expression profiles of crucial genes of carbohydrate metabolism between different stages and different treatments. SUSs, SUTs, TPSs, AGPLs, GBSSs and SSs showed significant responses during bulblet initiation. The expression of CWINs, SUTs and SWEETs were significantly upregulated by 2,4-D in two lilies. In addition, 2,4-D increased the expression of starch degradation genes (AMYs and BAMs) and inhibited starch synthesis genes (AGPLs, GBSSs and SSs). SBEs were significantly upregulated in Lbg but not in Lb. Significant co-expression was showed between genes involved in carbohydrate metabolism and auxin signaling, together with transcription factors such as bHLHs, MYBs, ERFs and C3Hs. This study indicates the coordinate regulation of bulblet initiation by carbohydrate metabolism and auxin signaling, serving as a basis for further studies on the molecular mechanism of bulblet initiation in lily and other bulbous flowers.

Changes of starch and sucrose content and related gene expression during the growth and development of Lanzhou lily bulb

As the main forms of carbohydrates, starch and sucrose play a vital role in the balance and coordination of various carbohydrates. Lanzhou lily is the most popular edible lily in China, mainly distributed in the central region of Gansu. To clarify the relationship between carbohydrate metabolism and bulb development of Lanzhou lily, so as to provide a basis for the promotion of the growth and development in Lanzhou lily and its important economic value, we studied lily bulbs in the squaring stage, flowering stage, half withering stage and withering stage. The plant height, fresh weight of mother and daughter bulbs continued to increase during the whole growth period and fresh weight of stem and leaf began to decrease in the half withering stage. The content of starch, sucrose and total soluble sugar in the lily mother bulb accumulated mostly in the flowering, withering and half withering stages, respectively. Starch, sucrose and total soluble sugar accumulated in the daughter bulb with the highest concentration during the withering stage. In the transcription level, sucrose synthase (SuSy1) and sucrose invertase (INV2) expressed the highest in squaring stage, and the expression was significantly higher in the mother bulb than in the daughter bulb. In flowering stage, the expression levels of soluble starch synthase (SSS1), starch-branching enzyme (SBE) and adenosine diphosphate-glucose pyrophosphorylase (AGP1) genes were higher in the mother bulb than in the daughter bulb. Altogether, our results indicate that starch and sucrose are important for the bulb growth and development of Lanzhou lily.

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.

Genome-Wide Transcriptomic Identification and Functional Insight of Lily WRKY Genes Responding to Botrytis Fungal Disease

Genome-wide transcriptome analysis using RNA-Seq of Lilium longiflorum revealed valuable genes responding to biotic stresses. WRKY transcription factors are regulatory proteins playing essential roles in defense processes under environmental stresses, causing considerable losses in flower quality and production. Thirty-eight WRKY genes were identified from the transcriptomic profile from lily genotypes, exhibiting leaf blight caused by Botrytis elliptica. Lily WRKYs have a highly conserved motif, WRKYGQK, with a common variant, WRKYGKK. Phylogeny of LlWRKYs with homologous genes from other representative plant species classified them into three groups- I, II, and III consisting of seven, 22, and nine genes, respectively. Base on functional annotation, 22 LlWRKY genes were associated with biotic stress, nine with abiotic stress, and seven with others. Sixteen unique LlWRKY were studied to investigate responses to stress conditions using gene expression under biotic and abiotic stress treatments. Five genes-LlWRKY3, LlWRKY4, LlWRKY5, LlWRKY10, and LlWRKY12-were substantially upregulated, proving to be biotic stress-responsive genes in vivo and in vitro conditions. Moreover, the expression patterns of LlWRKY genes varied in response to drought, heat, cold, and different developmental stages or tissues. Overall, our study provides structural and molecular insights into LlWRKY genes for use in the genetic engineering in Lilium against Botrytis disease.

Change in Sucrose Cleavage Pattern and Rapid Starch Accumulation Govern Lily Shoot-to-Bulblet Transition in vitro

In bulb crops, bulbing is a key progress in micropropagation and is the feature that most distinguishes bulbous crops from other plants. Generally, bulbing involves a shoot-to-bulblet transition; however, the underlying mechanism remains elusive. We explored this process by tracking the shoot-to-bulblet transition under different culture conditions. Rapid starch accumulation occurred at 15 days after transplanting (DAT) in the bulblet-inducing treatments as confirmed via histological observations and the significant elevation of starch synthesis related-gene transcription, including LohAGPS, LohAGPL, LohGBSS, LohSS, and LohSBE. However, for shoots that did not transition to bulblets and maintained the shoot status, much higher soluble sugars were detected. Interestingly, we observed a clear shift from invertase-catalyzed to sucrose synthase-catalyzed sucrose cleavage pattern based on the differential expression of LohCWIN and LohSuSy during the key transition stage (prior to and after bulbing at 0-15 DAT). Shoots that transitioned into bulblets showed significantly higher LohSuSy expression, especially LohSuSy4 expression, than shoots that did not transition. A symplastic phloem unloading pathway at the bulblet emergence stage (15 DAT) was verified via the 6(5)-carboxyfluorescein diacetate fluorescent tracer. We propose that starch is the fundamental compound in the shoot-to-bulblet transition and that starch synthesis is likely triggered by the switch from apoplastic to symplastic sucrose unloading, which may be related to sucrose depletion. Furthermore, this study is the first to provide a complete inventory of the genes involved in starch metabolism based on our transcriptome data. Two of these genes, LohAGPS1.2b and LohSSIIId, were verified by rapid amplification of cDNA ends cloning, and these data will provide additional support for Lilium research since whole genome is currently lacking.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

Composition and diurnal variation of floral scent emission in Rosa rugosa Thunb. and Tulipa gesneriana L

This study was aimed to explore the composition and diurnal variation analyses of floral scent emission from Rosa rugosa Thunb. and Tulipa gesneriana L. The floral scent from the fresh flower were collected at different time points (9:00, 12:00, 15:00, 18:00, and 21:00) using dynamic headspace collection and were analyzed using autothermal desorber-gas chromatography/mass spectrometry (ATD-GC/MS). The results showed that a total of 62 volatile flavor compounds were detected from Rosa rugosa Thunb and a total of 70 volatile flavor compounds were detected from Tulipa gesneriana L. They were identified with eight functional categories: alcohols, fatty hydrocarbons, terpenes, aldehydes, ketones, esters, and other substances. The total release amount first decreased, and then increased with time, and arrived at the lowest at 15:00. The release amounts of different categories present distinct change patterns. Among the components, phenylethyl alcohol, citronellol, methylene chloride, hexane, and acetone showed relatively higher release amounts and were thought as the main components in floral scent of Rosa rugosa Thunb. Alpha-Farnesene, ethanol, pentadecane, beta-ocimene, longifolene, caryophyllene, and acetone showed relatively higher release amounts and were thought as the main components in floral scent of Tulipa gesneriana L. Research of roses and tulips in aromatic in the garden provides a theoretical basis and research and improvement of the aroma components of aroma.

De novo transcriptome sequencing and gene expression profiling with/without B-chromosome plants of Lilium amabile

Supernumerary B chromosomes were found in Lilium amabile (2n = 2x = 24), an endemic Korean lily that grows in the wild throughout the Korean Peninsula. The extra B chromosomes do not affect the host-plant morphology; therefore, whole transcriptome analysis was performed in 0B and 1B plants to identify differentially expressed genes. A total of 154,810 transcripts were obtained from over 10 Gbp data by de novo assembly. By mapping the raw reads to the de novo transcripts, we identified 7,852 differentially expressed genes (log₂FC > |10|), in which 4,059 and 3,794 were up-and down-regulated, respectively, in 1B plants compared to 0B plants. Functional enrichment analysis revealed that various differentially expressed genes were involved in cellular processes including the cell cycle, chromosome breakage and repair, and microtubule formation; all of which may be related to the occurrence and maintenance of B chromosomes. Our data provide insight into transcriptomic changes and evolution of plant B chromosomes and deliver an informative database for future study of B chromosome transcriptomes in the Korean lily.

Probing the function of protein farnesyltransferase in Tripterygium wilfordii

We found two subunits FTase/GGTaseI-α and FTase-β formed a heterodimer to transfer a farnesyl group from FPP to protein N-dansyl-GCVLS, confirming they are responsible for protein farnesylation in planta. Tripterygium wilfordii is a medicinal plant with a broad spectrum of anti-inflammatory, immunosuppressive and anti-cancer activities. Recently, a number of studies have focused on investigating the biosynthetic pathways of its bioactive compounds, whereas little attention has been paid to the enzymes which play important roles in regulating diverse developmental processes of T. wilfordii. In this study, we report for the first time the identification and characterization of two subunits of farnesyltransferase (FTase), farnesyltransferase/geranylgeranyltransferase I-α (TwFTase/GGTase I-α) and farnesyltransferase-β (TwFTase-β), in this important medicinal plant. Cell-free in vivo assays, yeast two-hybrid (Y2H) and pull-down assays showed that the two subunits interact with each other to form a heterodimer to perform the role of specifically transferring a farnesyl group from FPP to the CAAX-box protein N-dansyl-GCVLS. Furthermore, we discovered that the two subunits had the same cytoplasmic localization pattern and displayed the same tissue expression pattern. These results indicated that we identified a functional TwFTase enzyme which contains two functionally complementary subunits TwFTase/GGTase I-α and TwFTase-β, which provides us promising genetic targets to construct transgenic plants or screen for more adaptable T. wilfordii mutants, which are able to survive in changing environments.

RNA viromes of the oriental hybrid lily cultivar "Sorbonne"

BACKGROUND

The lily is a perennial flowering plant belonging to the genus Lilium in the family Liliaceae. Most cultivated lily plants are propagated by bulbs. Therefore, numerous lily bulbs are frequently infected by diverse viruses causing viral diseases. To date, no study has examined the viromes of plants of one type with identical genetic backgrounds collected from different geographical regions.

RESULTS

Here, we examined different viromes of the lily cultivar "Sorbonne" using 172 gigabytes of transcriptome data composed of 23 libraries from four different projects for the cultivar "Sorbonne." We identified 396 virus-associated contigs from all but one library. We identified six different viruses, including Plantago asiatica mosaic virus (PlAMV), Cucumber mosaic virus (CMV), Lily symptomless virus (LSV), Tulip virus X (TVX), Lily mottle virus (LMoV), and Tobacco rattle virus (TRV). Of them, PlAMV was the most common virus infecting the lily. Scale and flower samples possessed a high number of virus-associated reads. We assembled 32 nearly complete genomes for the six identified viruses possessing the polyadenylate tails. Genomes of all six viruses were highly conserved in the lily cultivar "Sorbonne" based on mutation analysis. We identified defective RNAs from LSV, TVX, and PlAMV localized in the triple gene block region. Phylogenetic analyses showed that virus genomes are highly correlated with geographical regions and host plants.

CONCLUSIONS

We conducted comprehensive virome analyses of a single lily cultivar, "Sorbonne," using transcriptome data. Our results shed light on an array of lily virome-associated topics, including virus identification, the dominant virus, virus accumulation in different plant tissues, virus genome assembly, virus mutation, identification of defective RNAs, and phylogenetic relationships of identified viruses. Taken together, we provide very useful methods and valuable results that can be applied in other virome-associated studies.

Two-dimensional analysis provides molecular insight into flower scent of Lilium 'Siberia'

Lily is a popular flower around the world not only because of its elegant appearance, but also due to its appealing scent. Little is known about the regulation of the volatile compound biosynthesis in lily flower scent. Here, we conducted an approach combining two-dimensional analysis and weighted gene co-expression network analysis (WGCNA) to explore candidate genes regulating flower scent production. In the approach, changes of flower volatile emissions and corresponding gene expression profiles at four flower developmental stages and four circadian times were both captured by GC-MS and RNA-seq methods. By overlapping differentially-expressed genes (DEGs) that responded to flower scent changes in flower development and circadian rhythm, 3,426 DEGs were initially identified to be candidates for flower scent production, of which 1,270 were predicted as transcriptional factors (TFs). The DEGs were further correlated to individual flower volatiles by WGCNA. Finally, 37, 41 and 90 genes were identified as candidate TFs likely regulating terpenoids, phenylpropanoids and fatty acid derivatives productions, respectively. Moreover, by WGCNA several genes related to auxin, gibberellins and ABC transporter were revealed to be responsible for flower scent production. Thus, this strategy provides an important foundation for future studies on the molecular mechanisms involved in floral scent production.