Transcriptome and metabolome reveal the accumulation of secondary metabolites in different varieties of Cinnamomum longepaniculatum

A Transcriptomic and Metabolomic Study on the Biosynthesis of Iridoids in Phlomoides rotata from the Qinghai-Tibet Plateau

Phlomoides rotata is a traditional Chinese herbal medicine that grows in the Qinghai-Tibet Plateau region at a 3100-5000 m altitude. Iridoid compounds are the main active compounds of the P. rotata used as medical ingredients and display anti-inflammatory, analgesic, and hepatoprotective properties. To better understand the biological mechanisms of iridoid compounds in this species, we performed a comprehensive analysis of the transcriptome and metabolome of P. rotata leaves from four different regions (3540-4270 m). Global metabolome profiling detected 575 metabolites, and 455 differentially accumulated metabolites (DAMs) were detected in P. rotata leaves from the four regions. Eight major DAMs related to iridoid metabolism in P. rotata leaves were investigated: shanzhiside methyl ester, 8-epideoxyloganic acid, barlerin, shanzhiside, geniposide, agnuside, feretoside, and catalpin. In addition, five soil physical and chemical indicators in P. rotata rhizosphere soils were analyzed. Four significant positive correlations were observed between alkaline nitrogen and geniposide, exchangeable calcium and geniposide, available potassium and shanzhiside, and available phosphorus and shanzhiside methyl ester. The transcriptome data showed 12 P. rotata cDNA libraries with 74.46 Gb of clean data, which formed 29,833 unigenes. Moreover, 78.91% of the unigenes were annotated using the eight public databases. Forty-one candidate genes representing 23 enzymes involved in the biosynthesis of iridoid compounds were identified in P. rotata leaves. Moreover, the DXS1, IDI1, 8-HGO1, and G10H2 genes associated with iridoid biosynthesis were specifically expressed in P. rotata. The integration of transcriptome and metabolome analyses highlights the crucial role of soil physical and chemical indicators and major gene expression related to iridoid metabolism pathways in P. rotata from different areas. Our findings provide a theoretical foundation for exploring the molecular mechanisms underlying iridoid compound accumulation in P. rotata.

Transcriptomic and Metabolomic Analyses Reveal the Response Mechanism of Ophiopogon japonicus to Waterlogging Stress

Ophiopogon japonicus, a plant that thrives in river alluvial dams, often faces waterlogging stress due to sustained rainfall and flood seasons, which significantly impacts its growth and development. Currently, the mechanisms of waterlogging tolerance in Ophiopogon japonicus are still unclear. This study analyzed the transcriptome and metabolome data for Ophiopogon japonicus in the Sichuan region (referred to as CMD) under varying degrees of waterlogging stress: mild, moderate, and severe. The results indicate that the group exposed to flooding stress exhibited a higher number of differentially expressed genes (DEGs) compared to the control group. Notably, most DEGs were downregulated and primarily enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant hormone signal transduction pathways. A total of 5151 differentially accumulated metabolites (DAMs) were identified, with significantly upregulated DAMs annotated to two clusters, namely flavonoids such as apiin, pelargonin, and others. Furthermore, our study revealed significant upregulation in the expression of C2H2 (C2H2 zinc finger proteins) and AP2/ERF-ERF (the subfamily ERF proteins of APETALA2/ethylene-responsive element binding factors) transcription factors in CMD under flooding stress, suggesting their critical roles in enabling CMD to adapt to these conditions. In conclusion, this research provides insights into the intricate molecular mechanisms underlying CMD's response to flooding stress and reports valuable genetic data for the development of transgenic plants with improved waterlogging tolerance.

Integrative transcriptome and metabolome analysis reveals the discrepancy in the accumulation of active ingredients between Lycium barbarum cultivars

MAIN CONCLUSION

The combined analysis of transcriptome and metabolome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum. Lycium barbarum L. has a high concentration of active ingredients and is well known in traditional Chinese herbal medicine for its therapeutic properties. However, there are many Lycium barbarum cultivars, and the content of active components varies, resulting in inconsistent quality between Lycium barbarum cultivars. At present, few research has been conducted to reveal the difference in active ingredient content among different cultivars of Lycium barbarum at the molecular level. Therefore, the transcriptome of 'Ningqi No.1' and 'Qixin No.1' during the three development stages (G, T, and M) was constructed in this study. A total of 797,570,278 clean reads were obtained. Between the two types of wolfberries, a total of 469, 2394, and 1531 differentially expressed genes (DEGs) were obtained in the 'G1 vs. G10,' 'T1 vs. T10,' and 'M1 vs. M10,' respectively, and were annotated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology identifiers. Using these transcriptome data, most DEGs related to the metabolism of the active ingredients in 'Ningqi No.1' and 'Qixin No.1' were identified. Moreover, a widely targeted metabolome analysis of the metabolites of 'Ningqi 1' and 'Qixin 1' fruits at the maturity stage revealed 1,135 differentially expressed metabolites (DEMs) in 'M1 vs. M10,' and many DEMs were associated with active ingredients such as flavonoids, alkaloids, terpenoids, and so on. We further quantified the flavonoid, lignin, and carotenoid contents of the two Lycium barbarum cultivars during the three developmental stages. The present outcome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum, which would provide the basic data for the formation of Lycium barbarum fruit quality and the breeding of outstanding strains.

The phenotypic variation mechanisms of Atractylodes lancea post-cultivation revealed by conjoint analysis of rhizomic transcriptome and metabolome

The wild Atractylodes lancea rhizomes have been traditionally used as herbal medicine. As the increasingly exhaustion of wild A. lancea, the artificial cultivation mainly contributed to the medicinal material production. However, besides the phenotypic variation of rhizome phenotypic trait alteration, the qualities of cultivated A. lancea decrease compared with the wild counterpart. To unveil the physiological and molecular mechanism beneath the phenotypic variation, GC-MS-based volatile organic compounds (VOCs) profiling and RNAseq-based transcriptome analysis were conducted. The volatile metabolomics profiling revealed 65 differentially accumulated metabolites (DAMs) while the transcriptomic profiling identified 12 009 differentially expressed unigenes (DEGs) post-cultivation. The volatile active compounds including atractylone, and eudesmol accumulated more in wild rhizome than in the cultivated counterpart, and several unigenes in terpene synthesis were downregulated under cultivated condition. Compared with the wild A. lancea rhizome, the contents of bioactive Jasmonic Acid (JAs) in cultivated A. lancea rhizome were higher, and evidences that JAs negatively regulate the terpenes biosynthesis in the cultivated A. lancea rhizome were also provided. The combinational omics analysis further indicated the high correlation between the ten cultivation-suppressed VOCs and the cultivation-altered genes for sesquiterpenoids biosynthesis in A. lancea. The network of the cultivation-altered transcription factors (TFs) and the ten VOCs suggested TFs (e.g. Arabidopsis ERF13 homologs and WRKY50) are involved in the regulation of terpenes biosynthesis. These results laid a theoretical basis for developing geo-herbalism medicinal plants with "high quality and optimal shape".

Comparative Analysis and Identification of Terpene Synthase Genes in Convallaria keiskei Leaf, Flower and Root Using RNA-Sequencing Profiling

The 'Lilly of the Valley' species, Convallaria, is renowned for its fragrant white flowers and distinctive fresh and green floral scent, attributed to a rich composition of volatile organic compounds (VOCs). However, the molecular mechanisms underlying the biosynthesis of this floral scent remain poorly understood due to a lack of transcriptomic data. In this study, we conducted the first comparative transcriptome analysis of C. keiskei, encompassing the leaf, flower, and root tissues. Our aim was to investigate the terpene synthase (TPS) genes and differential gene expression (DEG) patterns associated with essential oil biosynthesis. Through de novo assembly, we generated a substantial number of unigenes, with the highest count in the root (146,550), followed by the flower (116,434) and the leaf (72,044). Among the identified unigenes, we focused on fifteen putative ckTPS genes, which are involved in the synthesis of mono- and sesquiterpenes, the key aromatic compounds responsible for the essential oil biosynthesis in C. keiskei. The expression of these genes was validated using quantitative PCR analysis. Both DEG and qPCR analyses revealed the presence of ckTPS genes in the flower transcriptome, responsible for the synthesis of various compounds such as geraniol, germacrene, kaurene, linalool, nerolidol, trans-ocimene and valencene. The leaf transcriptome exhibited genes related to the biosynthesis of kaurene and trans-ocimene. In the root, the identified unigenes were associated with synthesizing kaurene, trans-ocimene and valencene. Both analyses indicated that the genes involved in mono- and sesquiterpene biosynthesis are more highly expressed in the flower compared to the leaf and root. This comprehensive study provides valuable resources for future investigations aiming to unravel the essential oil-biosynthesis-related genes in the Convallaria genus.

Transcriptomic analysis revealed the candidate metabolic pathways and genes associated with cold tolerance in a mutant without anthocyanin accumulation in common vetch (Vicia sativa L.)

Common vetch (Vicia sativa L.) is a leguminous crop used to feed livestock with vegetative organs or fertilize soils by returning to the field. Survival of fall-seeded plants is often affected by freezing damage during overwintering. This study aims to investigate the transcriptomic profiling in response to cold in a mutant with reduced accumulation of anthocyanins under normal growth and low-temperature conditions for understanding the underlying mechanisms. The mutant had increased cold a tolerance with higher survival rate and biomass during overwintering compared to the wild type, which led to increased forage production. Transcriptomic analysis in combination with qRT-PCR and physiological measurements revealed that reduced anthocyanins accumulation in the mutant resulted from reduced expression of serial genes involving in anthocyanin biosynthesis, which led to the altered metabolism, with an increased accumulation of free amino acids and polyamines. The higher levels of free amino acids and proline in the mutant under low temperature were associated with improved cold tolerance. The altered expression of some genes involved in ABA and GA signaling was also associated with increased cold tolerance in the mutant.

Cytogenetic Characterization and Metabolomic Differences of Full-Sib Progenies of Saccharum spp

Sugarcane smut is a worldwide fungal disease. Disease resistance breeding is the most economical and effective measure to prevent and control sugarcane smut. The cytogenetic characteristics and metabolomic differences of sugarcane F₁s are closely related to disease resistance. Zhongzhe 1 and G160 sugarcane from the same parents (ROC25 and Yunzhe89-7) were used; the plants were grown in accordance with the barrel method. When the seedlings had 4-5 leaves, genomic in situ hybridization (GISH) was performed; digoxigenin (DIG)-labeled female parental (ROC25)DNA and biotin-labeled male parental (Yunzhe89-7) DNA were used as probes, and the karyotypes of two hybrids were analyzed. The new sugarcane smut-resistant variety (Zhongzhe 1) and the susceptible variety (G160) derived from the same parent were analyzed via gas chromatography-mass spectrometry technology (GC-MS) to compare the metabolomic differences between them. GISH analysis revealed that the chromosome ploidy number of Zhongzhe 1 sugarcane and G160 sugarcane were 114 and 110, respectively. However, the two contain different numbers of chromosomes from the female (ROC25) and male (Yunzhe89-7) parents. Moreover, 258 significantly changed metabolites were identified in smut-resistant Zhongzhe 1, as compared with the smut-susceptible G160 sugarcane: 56 flavonoids, 52 phenolic acids, 30 lipids, 26 organic acids, 26 amino acids and derivatives, 19 nucleotides and derivatives, 5 alkaloids, 9 terpenoids, and 35 others. Multivariate statistical analysis revealed a distinct difference in metabolic pathways between Zhongzhe 1 sugarcane and G160, and both of these varieties had unique functional metabolites. Differences in chromosome composition may constitute the genetic basis for the difference in resistance to smut disease between Zhongzhe 1 sugarcane and G160 sugarcane, and a high accumulation of flavonoids, lipids, terpenoids and tannins may constitute the basis of resistance to smut disease for the Zhongzhe 1 variety.

Full-length transcriptome sequencing reveals the molecular mechanism of monoterpene and sesquiterpene biosynthesis in Cinnamomum burmannii

Essential oil of Cinnamomum burmannii is rich in monoterpenes and sesquiterpenes and is widely used in cosmetics and medicines. Knowledge about the enzymes that catalyze the formation of monoterpenes and sesquiterpenes in C. burmannii is insufficient. Therefore, anatomy observation of C. burmannii at the four developmental stages (7 days, CBS1; 14 days, CBS2; 21 days, CBS3, and 28 days, CBS4) were conducted to elucidate the origins of essential oil production. Twelve full-length transcriptomes of C. burmannii leaves at the four stages were generated using Oxford Nanopore Technologies. GC-MS analysis revealed 15 monoterpene and sesquiterpenes dramatically increased from CBS1 to CBS4. A weighted correlation network analysis (WGCNA) in association and differentially expressed genes across four developmental stages were performed. A total of 44 differentially expressed genes (DEGs) were involved in terpenoid syntheses during leaf development. Among them, the DEGs of the mevalonate acid (MVA) pathway were predominantly expressed at CBS1, while those of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway showed increased expression from CBS2 to CBS4. Besides, fourteen genes were associated with monoterpene synthesis and nine with sesquiterpene synthesis. Functions of these DEGs were further predicted with regard to gene expression profile and phylogenetic relationship with those characterized in previous studies. In addition, 922 long noncoding RNAs (lncRNAs) were detected, of which twelve were predicted to regulate monoterpene and sesquiterpene biosynthesis. The present study provided new insights the molecular mechanisms of monoterpenoid and sesquiterpenoid syntheses of C. burmannii.

Functional genomics-enabled characterization of CYP81B140 and CYP81B141 from Plumbago zeylanica L. substantiates their involvement in plumbagin biosynthesis

Novel cytochrome P450s, CYP81B140 and CYP81B141 from Plumbago zeylanica were functionally characterized to understand their involvement in polyketide plumbagin biosynthesis. Further, we propose 3-methyl-1-8-naphthalenediol and isoshinanolone as intermediates for plumbagin biosynthesis. Plumbago zeylanica L. (P. zeylanica) is a medicinally important plant belonging to the family Plumbaginaceae. It comprises the most abundant naphthoquinone plumbagin having anti-cancer activity. Only the polyketide synthase (PKS) enzyme has been identified from the biosynthetic pathway which catalyzes iterative condensation of acetyl-CoA and malonyl-CoA molecules. The plumbagin biosynthesis involves hydroxylation, oxidation, hydration and dehydration of intermediate compounds which are expected to be catalyzed by cytochrome P450s (CYPs). To identify the CYPs, co-expression analysis was carried out using PKS as a candidate gene. Out of the eight identified CYPs, CYP81B140 and CYP81B141 have similar expression with PKS and belong to the CYP81 family. Phylogenetic analysis suggested that CYP81B140 and CYP81B141 cluster with CYPs from CYP81B, CYP81D, CYP81E and CYP81AA subfamilies which are known to be involved in the hydroxylation and oxidation reactions. Moreover, artificial microRNA-mediated transient individual silencing and co-silencing of CYP81B140 and CYP81B141 significantly reduced plumbagin and increased the 3-methyl-1-8-naphthalenediol and isoshinanolone content. Based on metabolite analysis, we proposed that 3-methyl-1-8-naphthalenediol and isoshinanolone function as intermediates for plumbagin biosynthesis. Transient silencing, over-expression and docking analysis revealed that CYP81B140 is involved in C-1 oxidation, C-4 hydroxylation and [C2-C3] hydration of 3-methyl-1-8-naphthalenediol to form isoshinanolone, whereas CYP81B141 is catalyzing [C2-C3] dehydration and C-4 oxidation of isoshinanolone to form plumbagin. Our results indicated that both CYP81B140 and CYP81B141 are promiscuous and necessary for plumbagin biosynthesis. This is the first report of identification and functional characterization of P. zeylanica-specific CYPs involved in plumbagin biosynthetic pathway and in general hexaketide synthesis in plants.

Preparation, Optimization, and Characterization of Inclusion Complexes of Cinnamomum longepaniculatum Essential Oil in β-Cyclodextrin

Cinnamomum longepaniculatum essential oil (CLEO) possesses antibacterial, anti-inflammatory, and antioxidant activities. However, CLEO shows volatilization and poor solubility, which limits its application field. In this research, inclusion complexes of β-cyclodextrin (β-CD) with CLEO were produced, and its physicochemical properties were characterized. Response surface methodology was used to obtain optimum preparation conditions. A statistical model was generated to define the interactions among the selected variables. Results show that the optimal conditions were an H2O/β-CD ratio of 9.6:1 and a β-CD/CLEO ratio of 8:1, with the stirring temperature of 20 °C for the maximal encapsulation efficiency values. The physicochemical properties of CLEO/β-CD inclusion complexes (CLEO/β-CD-IC) were investigated. Fourier transform infrared spectroscopy showed that correlative characteristic bands of CLEO disappeared in the inclusion complex. X-ray diffraction presented different sharp peaks at the diffraction angle of CLEO/β-CD-IC. The thermogravimetric analysis demonstrated the thermal stability of CLEO was enhanced after encapsulation. Tiny aggregates with a smaller size of CLEO/β-CD-IC particles were observed by scanning electron microscopy. The comparison of β-CD, CLEO, and physical mixtures with CLEO/β-CD-IC confirmed the formation of inclusion complexes.