A Comparison of the Flavonoid Biosynthesis Mechanisms of Dendrobium Species by Analyzing the Transcriptome and Metabolome

Full-length transcriptome sequencing provides new insights into the complexity of flavonoid biosynthesis in Glechoma longituba

Glechoma longituba has been frequently used in treating urolithiasis and cholelithiasis due to the presence of flavonoids, which are its major bioactive constituents. However, research on the molecular background of flavonoid biosynthesis in G. longituba is limited. In this study, we used single-molecule real-time combined with next-generation sequencing technologies to construct the complete transcriptome of G. longituba. We identified 404,648 non-redundant transcripts, including 249,697 coding sequences, 197,811 simple sequence repeats, 174,846 long noncoding RNA, and 176,554 coding RNA. Moreover, we functionally annotated 346,218 isoforms (85.56%) and identified 86,528 differentially expressed genes. We also identified 55 non-redundant full-length isoforms related to the flavonoid biosynthetic pathway. Pearson correlation analysis revealed that the expression levels of some key genes of the flavonoid biosynthesis pathway were significantly positively correlated with the flavonoid metabolites. Furthermore, we performed bioinformatics analysis (sequence and structural) of isoform_47029 (encoding flavanone 3-hydroxylase) and isoform_53692 (encoding flavonol synthase) to evaluate their potential biological functions. Finally, we validated gene expression levels of 12 flavonoid-related key enzyme genes using quantitative real-time PCR. Overall, this study provides full-length transcriptome information on G. longituba for the first time and valuable molecular resources for further research on the medicinal properties of this plant.

Multi-Omics Elucidates Difference in Accumulation of Bioactive Constituents in Licorice (Glycyrrhiza uralensis) under Drought Stress

Licorice is a frequently applied herb with potential edible and medicinal value based on various flavonoids and triterpenes. However, studies on detailed flavonoid and triterpene metabolism and the molecular basis of their biosynthesis in licorice are very limited, especially under drought conditions. In the present study, we carried out transcriptome, proteome, and metabolome experiments. To ultimately combine three omics for analysis, we performed a bioinformatics comparison, integrating transcriptome data and proteome data through a Cloud platform, along with a simplified biosynthesis of primary flavonoids and triterpenoids in the KEGG pathway based on metabolomic results. The biosynthesis pathways of triterpenes and flavonoids are enriched at both gene and protein levels. Key flavonoid-related genes (PAL, 4CL, CHS, CHI, CYP93C, HIDH, HI4OMT, and CYP81E1_7) and representative proteins (HIDH, CYP81E1_7, CYP93C, and VR) were obtained, which all showed high levels after drought treatment. Notably, one R2R3-MYB transcription factor (Glyur000237s00014382.1), a critical regulator of flavonoid biosynthesis, achieved a significant upregulated expression as well. In the biosynthesis of glycyrrhizin, both gene and protein levels of bAS and CYP88D6 have been found with upregulated expression under drought conditions. Most of the differentially expressed genes (DEGs) and proteins (DEPs) showed similar expression patterns and positively related to metabolic profiles of flavonoid and saponin. We believe that suitable drought stress may contribute to the accumulation of bioactive constituents in licorice, and our research provides an insight into the genetic study and quality breeding in this plant.

Uncovering nutritional metabolites and candidate genes involved in flavonoid metabolism in Houttuynia cordata through combined metabolomic and transcriptomic analyses

The perennial herb Houttuynia cordata has long been cultivated and used as medicinal and edible plant in Asia. Nowadays, increasing attention is attracted due to its numerous health benefits. Flavonoids are the main chemical constituents exerting pharmacological activities. In the present study, we investigated both metabolome and transcriptome of two H. cordata accessions (6# and 7#) with distinct flavonoids contents. In total 397 metabolites, i.e., 220 flavonoids, 92 amino acids and derivatives, 20 vitamins, and 65 saccharides were abundant in aboveground part. Cyanidin-3-O-rutinoside and quercetin-3-O-galactoside were the most abundant flavonoids, which can be categorized into seven classes, namely anthocyanidins, chalcones, flavanols, flavanones, flavanonols, flavones, and flavonols. Flavonols was the most abundant group. Contents of 112 flavonoids differed significantly between the two accessions, with catechin-(7,8-bc)-4α-(3,4-dihydroxyphenyl)-dihydro-2-(3H)-one, cinchonain Id, and cinchonain Ic being the dominant flavonoid metabolites among them. Pinocembrin-7-O-neohesperidoside, pinocembrin-7-O-rutinoside, and kaempferol-3-O-galactoside-4'-O-glucoside were uniquely abundant in accession 7. Transcriptome data revealed a total of 110 different expressed genes related to flavonoid metabolism, with more highly expressed genes observed in 7#. We annotated a total of 19 differential flavonoid metabolites and 34 differentially expressed genes that are associated with the flavonoid metabolic network. Based on the transcriptome and qPCR data a total of 8 key candidate genes involved in flavonoid metabolism were identified. The ANS gene were found to play an important role in the synthesis of cyanidin-3-O-glucoside, while the CHI, F3'H and FLS genes were mainly responsible for controlling the levels of flavanones, flavones, and flavonols, respectively. Collectively, the present study provides important insights into the molecular mechanism underlying flavonoid metabolism in H. cordata.

Selection of Reference Genes in Evodia rutaecarpa var. officinalis and Expression Patterns of Genes Involved in Its Limonin Biosynthesis

E. rutaecarpa var. officinalis is a traditional Chinese medicinal plant known for its therapeutic effects, which encompass the promotion of digestion, the dispelling of cold, the alleviation of pain, and the exhibition of anti-inflammatory and antibacterial properties. The principal active component of this plant, limonin, is a potent triterpene compound with notable pharmacological activities. Despite its significance, the complete biosynthesis pathway of limonin in E. rutaecarpa var. officinalis remains incompletely understood, and the underlying molecular mechanisms remain unexplored. The main purpose of this study was to screen the reference genes suitable for expression analysis in E. rutaecarpa var. officinalis, calculate the expression patterns of the genes in the limonin biosynthesis pathway, and identify the relevant enzyme genes related to limonin biosynthesis. The reference genes play a pivotal role in establishing reliable reference standards for normalizing the gene expression data, thereby ensuring precision and credibility in the biological research outcomes. In order to identify the optimal reference genes and gene expression patterns across the diverse tissues (e.g., roots, stems, leaves, and flower buds) and developmental stages (i.e., 17 July, 24 August, 1 September, and 24 October) of E. rutaecarpa var. officinalis, LC-MS was used to analyze the limonin contents in distinct tissue samples and developmental stages, and qRT-PCR technology was employed to investigate the expression patterns of the ten reference genes and eighteen genes involved in limonin biosynthesis. Utilizing a comprehensive analysis that integrated three software tools (GeNorm ver. 3.5, NormFinder ver. 0.953 and BestKeeper ver. 1.0) and Delta Ct method alongside the RefFinder website, the best reference genes were selected. Through the research, we determined that Act1 and UBQ served as the preferred reference genes for normalizing gene expression during various fruit developmental stages, while Act1 and His3 were optimal for different tissues. Using Act1 and UBQ as the reference genes, and based on the different fruit developmental stages, qRT-PCR analysis was performed on the pathway genes selected from the "full-length transcriptome + expression profile + metabolome" data in the limonin biosynthesis pathway of E. rutaecarpa var. officinalis. The findings indicated that there were consistent expression patterns of HMGCR, SQE, and CYP450 with fluctuations in the limonin contents, suggesting their potential involvement in the limonin biosynthesis of E. rutaecarpa var. officinalis. This study lays the foundation for further research on the metabolic pathway of limonin in E. rutaecarpa var. officinalis and provides reliable reference genes for other researchers to use for conducting expression analyses.

Flavonoid synthesis in Lamiophlomis rotata from Qinghai-Tibet Plateau is influenced by soil properties, microbial community, and gene expression

Lamiophlomis rotata is a medicinal plant in Qinghai-Tibet Plateau, in which flavonoid compounds are the major medicinal components. However, it remains unclear how flavonoid metabolism of L. rotata is influenced by soil properties and microbial community. In this study, we collected L. rotata seedlings and rhizosphere soils from five habitats ranging from 3750 to 4270 m of altitude and analyzed the effects of habitat conditions on flavonoid metabolism. The activities of peroxidase, cellulase, and urease were increased with altitude, while those of alkaline phosphatase, alkaline protease, and sucrase were decreased with altitude. Analysis of OTUs showed that the total number of bacterial genera was higher than that of fungal genera. The highest number of fungal genera was 132, and that of bacterial genera was 33 in Batang (BT) town in Yushu County at an altitude of 3880 m, suggesting that the fungal communities may play a critical role in L. rotata rhizosphere soils. Flavonoids in leaves and roots of L. rotata shared a similar pattern, with a trend of increasing levels with altitude. The highest flavonoid content measured, 12.94 mg/g in leaves and 11.43 mg/g in roots, was from Zaduo (ZD) County at an altitude of 4208 m. Soil peroxidases affected quercetin content in leaves of L. rotata, while the fungus Sebacina affected flavonoid content in leaves and roots of L. rotata. The expression of PAL, F3'H, FLS, and FNS genes showed a declining trend in leaves with altitude, while F3H showed an increasing trend in both leaves and roots. Overall, soil physicochemical properties and microbial community affect flavonoid metabolism in L. rotata in Qinghai-Tibet Plateau. The variations in flavonoid content and gene expression as well as their associations with soil factors revealed the complexity of the growth conditions and genetic makeup in L. rotata habitats of Qinghai-Tibet Plateau.

Transcriptomic and Metabolomic Analyses Reveal Differences in Flavonoid Pathway Gene Expression Profiles between Two Dendrobium Varieties during Vernalization

Dendrobium (Orchidaceae, Epidendoideae) plants have flowers with a wide variety of colors that persist for a long period throughout the year. The yellow coloration of Dendrobium flowers is mainly determined by the flavonol pathway and the flavone pathway, but the relevant biosynthesis mechanisms during vernalization remain unclear. To explore the similarities and differences in flavonoid biosynthesis in different tissues during vernalization, we selected two species of Dendrobium for a flower color study: Dendrobium capillipes Rchb (which has yellow flowers) and Dendrobium nobile Lindl (which has white flowers). We collected a total of 36 samples from six tissue types and both Dendrobium species during vernalization and subjected the samples to metabolic profiling and transcriptome sequencing. A total of 31,504 differentially expressed genes (DEGs) were identified between different tissues of the two Dendrobium species by transcriptomic analysis. However, many differentially accumulated metabolites (DAMs) and DEGs were enriched not only in the general pathway of "flavonoid biosynthesis" but also in multiple subpathways of "flavone and flavonol biosynthesis". According to a combined transcriptome and metabolome analysis, Putrescine hydroxycinnamoyl transferase 1 (LOC110093422) may be the main gene responsible for the differences in flavonoid accumulation during vernalization, which is closely associated with yellow flowers. Taken together, the results of our study preliminarily revealed the metabolites responsible for and the key genes regulating flavonoid biosynthesis during vernalization. These results provide a basis for the further study of the molecular mechanism of flavonoid synthesis during vernalization.

Integration of high-throughput omics technologies in medicinal plant research: The new era of natural drug discovery

Medicinal plants are natural sources to unravel novel bioactive compounds to satisfy human pharmacological potentials. The world's demand for herbal medicines is increasing year by year; however, large-scale production of medicinal plants and their derivatives is still limited. The rapid development of modern technology has stimulated multi-omics research in medicinal plants, leading to a series of breakthroughs on key genes, metabolites, enzymes involved in biosynthesis and regulation of active compounds. Here, we summarize the latest research progress on the molecular intricacy of medicinal plants, including the comparison of genomics to demonstrate variation and evolution among species, the application of transcriptomics, proteomics and metabolomics to explore dynamic changes of molecular compounds, and the utilization of potential resources for natural drug discovery. These multi-omics research provide the theoretical basis for environmental adaptation of medicinal plants and allow us to understand the chemical diversity and composition of bioactive compounds. Many medicinal herbs' phytochemical constituents and their potential health benefits are not fully explored. Given their large diversity and global distribution as well as the impacts of growth duration and environmental factors on bioactive phytochemicals in medicinal plants, it is crucial to emphasize the research needs of using multi-omics technologies to address basic and applied problems in medicinal plants to aid in developing new and improved medicinal plant resources and discovering novel medicinal ingredients.