Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata

Lonicera caerulea genome reveals molecular mechanisms of freezing tolerance and anthocyanin biosynthesis

INTRODUCTION

Lonicera caerulea L. (blue honeysuckle) is a noteworthy fleshy-fruited tree and a prominent medicinal plant, which possesses notable characteristics such as exceptional resilience to winter conditions and early maturation, and the richest source of functional anthocyanins, particularly cyanidin-3-glucoside. The molecular mechanisms responsible for its freezing tolerance and anthocyanin biosynthesis remain largely unknown.

OBJECTIVES

Here, a chromosome-scale genome of L. caerulea was presented, aiming to examine the genetic foundations that underlie these characteristics of blue honeysuckle.

METHODS

The PacBio HiFi reads and Hi-C data were used to construct high-quality genome of blue honeysuckle. Comparative genomic and transcriptomic analyses were conducted to elucidate the molecular mechanisms of freezing tolerance and anthocyanin biosynthesis.

RESULTS

Comparative genomics analysis between L. caerulea and L. japonica revealed that the dynamic changes of duplicated genes contributed to their phytochemical reconstruction and environmental adaptation. Moreover, the ABA and ICE-CBF-COR signaling pathways were closely correlated to the freezing tolerance of L. caerulea. Genome-wide identification and biochemical function indicated that three anthocyanin 3',5'-O-methyltransferases (LcOMT2, LcOMT14, and LcOMT20) and two 3'-O-glycosyltransferases (LcUGT78X1 and LcUGT95P1) were responsible for anthocyanin biosynthesis. In addition, LcUGT78X1 was regarded as the potent glycosyltransferase for the accumulation of cyanidin-3-glucoside in L. caerulea.

CONCLUSION

This research elucidates the crucial roles of the ABA and ICE-CBF-COR signaling pathways in enhancing freezing tolerance, while also identifying highly efficient anthocyanin biosynthetic enzymes in L. caerulea. These findings advance the understanding of environmental adaptation and phytochemical production in Lonicera species.

Genome-Wide Identification and Characterization of OSC Gene Family in Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum is a traditional Chinese medicinal plant of considerable application value and commercial potential, primarily due to its production of various bioactive compounds, particularly dammarane-type triterpenoid saponins that are structurally analogous to ginsenosides. Oxidosqualene cyclase (OSC), a pivotal enzyme in the biosynthesis of triterpenoid metabolites in plants, catalyzes the conversion of oxidosqualene into triterpenoid precursors, which are essential components of the secondary metabolites found in G. pentaphyllum. To elucidate the role of OSC gene family members in the synthesis of gypenosides within G. pentaphyllum, this study undertook a comprehensive genome-wide identification and characterization of OSC genes within G. pentaphyllum and compared their expression levels across populations distributed over different geographical regions by both transcriptome sequencing and qRT-PCR experimental validation. The results identified a total of 11 members of the OSC gene family within the genome of G. pentaphyllum. These genes encode proteins ranging from 356 to 767 amino acids, exhibiting minor variations in their physicochemical properties, and are localized in peroxisomes, cytoplasm, plasma membranes, and lysosomes. All GpOSCs contain highly conserved DCTAE and QW sequences that are characteristic of the OSC gene family. A phylogenetic analysis categorized the GpOSCs into four distinct subfamilies. A cis-element analysis of the GpOSC promoters revealed a substantial number of abiotic stress-related elements, indicating that these genes may respond to drought conditions, low temperatures, and anaerobic environments, thus potentially contributing to the stress resistance observed in G. pentaphyllum. Expression analyses across different G. pentaphyllum populations demonstrated significant variability in OSC gene expression among geographically diverse samples of G. pentaphyllum, likely attributable to genetic variation or external factors such as environmental conditions and soil composition. These differences may lead to the synthesis of various types of gypenosides within geographically distinct G. pentaphyllum populations. The findings from this study enhance our understanding of both the evolutionary history of the OSC gene family in G. pentaphyllum and the biosynthetic mechanisms underlying triterpenoid compounds. This knowledge is essential for investigating molecular mechanisms involved in forming dammarane-type triterpenoid saponins as well as comprehending geographical variations within G. pentaphyllum populations. Furthermore, this research lays a foundation for employing plant genetic engineering techniques aimed at increasing gypenoside content.

Integrated genomic, transcriptomic, and metabolomic analyses of Ilex hylonoma provide insights into the triterpenoid saponin biosynthesis

Ilex is known for its rich content of secondary metabolites, particularly triterpenoid saponins. These compounds hold significant value in natural remedies and herbal medicine. However, the molecular mechanisms responsible for triterpenoid biosynthesis in plants of this genus remain largely unexplored. In this study, we successfully generated the first chromosome-scale genome of Ilex hylonoma. The assembly, comprising 20 anchored chromosomes, has an N50 contig size of 2.13 Mb and a scaffold size of 33.68 Mb. Comparative genome analyses with two other congeners with available chromosome-level genomes suggested that an end-to-end chromosome fusion event likely contributed to the reduction in chromosome number from n = 20 to n = 19 within this genus. By integrating transcriptomic and metabolomic data, we identified the gene expression patterns and metabolite profiles of I. hylonoma across three commonly utilized medicinal tissues. We subsequently pinpointed candidate genes involved in the regulation of triterpenoid saponin biosynthesis, including CYP450 genes, UGT genes, and associated transcription factors. Furthermore, yeast heterologous expression analysis revealed that ihyl08363 catalyzed the conversion of β-amyrin into oleanolic acid, while ihyl04303 catalyzed the C-2α hydroxylation of oleanolic acid to produce maslinic acid. This integrated analysis provides valuable insights into the biosynthesis of important triterpenoid saponins in medicinal Ilex plants.

A Novel Biosynthetic Strategy for Ginsenoside Ro: Construction of a Metabolically Engineered Saccharomyces cerevisiae Strain Using a Newly Identified UGAT Gene from Panax ginseng as the Key Enzyme Gene and Optimization of Fermentation Conditions

Ginsenoside Ro, as one of the few oleanane-type ginsenosides, is well known for its unique molecular structure and biological activities. Currently, research on the biosynthesis of ginsenoside Ro is still in its early stages. Therefore, the establishment of a new ginsenoside Ro cell factory is of great significance for the in-depth development and utilization of genes related to ginsenoside Ro synthesis, as well as for the exploration of pathways to obtain ginsenoside Ro. In this study, we cloned endogenous constitutive promoters, terminators, and other genetic elements from S. cerevisiae BY4741. These elements were then sequentially assembled with the uridine diphosphate glucuronic acid transferase gene identified in our previously study (PgUGAT252645) and several other reported key enzyme genes, to construct DNA fragments used for integration into the genome of S. cerevisiae BY4741. By sequentially transferring these DNA fragments into chemically competent cells of engineering strains and conducting screening and target product detection, we successfully constructed an engineered S. cerevisiae strain (BY-Ro) for ginsenoside Ro biosynthesis using S. cerevisiae BY4741 as the host cell. Strain BY-Ro produced 253.32 μg/L of ginsenoside Ro under optimal fermentation conditions. According to subsequent measurements and calculations, this equates to 0.033 mg/g DCW, corresponding to approximately 31% of the ginsenoside Ro content found in plant samples. This study not only included a deeper investigation into the function of PgUGAT252645 but also provides a novel engineering platform for ginsenoside Ro biosynthesis.

Metabolomics combined with network pharmacology reveals the potential development value of Campanumoea javanica Bl. and its metabolite differences with Codonopsis Radix

Campanumoea javanica Bl. (CJ) traditionally used in Southwestern China, is now widely consumed as a health food across the nation. Due to its similar efficacy to Codonopsis Radix (CR) and their shared botanical family, CJ is often used as a substitute for CR. According to the Chinese Pharmacopoeia, Codonopsis pilosula var. modesta (Nannf.) L.T. Shen (CPM), Codonopsis pilosula (Franch.) Nannf. (CP), and Codonopsis tangshen Oliv. (CT) are the primary sources of CR. However, details on the differences in composition, effectiveness, and compositional between CJ and CR are still limited. Besides, there is little evidence to support the application of CJ as a drug. In this study, we employed widely targeted metabolomics, network pharmacology analysis, and molecular docking to explore the disparities in metabolite profiles between CJ and CR and to predict the pharmacological mechanisms of the dominant differential metabolites of CJ and their potential medicinal applications. The widely targeted metabolomics results indicated that 1,076, 1,102, 1,102, and 1,093 compounds, most phenolic acids, lipids, amino acids, and flavonoids, were characterized in CJ, CPM, CP, and CT, respectively. There were an average of 1061 shared compounds in CJ and CRs, with 95.07% similarity in metabolic profiles. Most of the metabolites in CJ were previously unreported. Twelve of the seventeen dominant metabolites found in CJ were directly associated with treating cancer and lactation, similar to the traditional medicinal efficacy. The molecular docking results showed that the dominant metabolites of CJ had good docking activity with the core targets PIK3R1, PIK3CA, ESR1, HSP90AA1, EGFR, and AKT1. This study provides a scientific basis for understanding the similarities and differences between CJ and CR at the metabolome level, offering a theoretical foundation for developing innovative medications from CJ. Additionally, it significantly enhances the metabolite databases for both CJ and CR.

Genome-Wide Identification of the CYP716 Gene Family in Platycodon grandiflorus (Jacq.) A. DC. and Its Role in the Regulation of Triterpenoid Saponin Biosynthesis

The main type of saponins occurring in the root of Platycodon grandiflorus (Jacq.) A. DC. are oleanolic acid glycosides. The CYP716 gene family plays a major role in catalyzing the conversion of β-amyrin into oleanolic acid. However, studies on the CYP716 genes in P. grandiflorus are limited, and its evolutionary history remains poorly understood. In this study, 22 PgCYP716 genes were identified, distributed among seven subfamilies. Cis-acting elements of the PgCYP716 promoters were mainly involved in plant hormone regulation and responses to abiotic stresses. PgCYP716A264, PgCYP716A391, PgCYP716A291, and PgCYP716BWv3 genes were upregulated in the root and during saponin accumulation, as shown by RNA-seq analysis, suggesting that these four genes play an important role in saponin synthesis. The results of subcellular localization indicated that these four genes encoded membrane proteins. Furthermore, the catalytic activity of these four genes was proved in the yeast, which catalyzed the conversion of β-amyrin into oleanolic acid. We found that the content of β-amyrin, platycodin D, platycoside E, platycodin D3, and total saponins increased significantly when either of the four genes was over expressed in the transgenic hair root. In addition, the expression of PgSS, PgGPPS2, PgHMGS, and PgSE was also upregulated while these four genes were overexpressed. These data support that these four PgCYP716 enzymes oxidize β-amyrin to produce oleanolic acid, ultimately promoting saponin accumulation by activating the expression of upstream pathway genes. Our results enhanced the understanding of the functional variation among the PgCYP716 gene family involved in triterpenoid biosynthesis and provided a theoretical foundation for improving saponin content and enriching the saponin biosynthetic pathway in P. grandiflorus.

Transcriptome and metabolome analysis revealed the dynamic change of bioactive compounds of Fructus Ligustri Lucidi

BACKGROUND

The Fructus Ligustri Lucidi, the fruit of Ligustrum lucidum, contains a variety of bioactive compounds, such as flavonoids, triterpenoids, and secoiridoids. The proportions of these compounds vary greatly during the different fruit development periods of Fructus Ligustri Lucidi. However, a clear understanding of how the proportions of the compounds and their regulatory biosynthetic mechanisms change across the different fruit development periods of Fructus Ligustri Lucidi is still lacking.

RESULTS

In this study, metabolite profiling and transcriptome analysis of six fruit development periods (45 DAF, 75 DAF, 112 DAF, 135 DAF, 170 DAF, and 195 DAF) were performed. Seventy compounds were tentatively identified, of which secoiridoids were the most abundant. Eleven identified compounds were quantified by high performance liquid chromatography. A total of 103,058 unigenes were obtained from six periods of Fructus Ligustri Lucidi. Furthermore, candidate genes involved in triterpenoids, phenylethanols, and oleoside-type secoiridoid biosynthesis were identified and analyzed. The in vitro enzyme activities of nine glycosyltransferases involved in salidroside biosynthesis revealed that they can catalyze trysol and hydroxytyrosol to salidroside and hydroxylsalidroside.

CONCLUSIONS

These results provide valuable information to clarify the profile and molecular regulatory mechanisms of metabolite biosynthesis, and also in optimizing the harvest time of this fruit.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

The high-quality genome of Grona styracifolia uncovers the genomic mechanism of high levels of schaftoside, a promising drug candidate for treatment of COVID-19

Recent study has evidenced that traditional Chinese medicinal (TCM) plant-derived schaftoside shows promise as a potential drug candidate for COVID-19 treatment. However, the biosynthetic pathway of schaftoside in TCM plants remains unknown. In this study, the genome of the TCM herb Grona styracifolia (Osbeck) H.Ohashi & K.Ohashi (GSO), which is rich in schaftoside, was sequenced, and a high-quality assembly of GSO genome was obtained. Our findings revealed that GSO did not undergo recent whole genome duplication (WGD) but shared an ancestral papilionoid polyploidy event, leading to the gene expansion of chalcone synthase (CHS) and isoflavone 2'-hydroxylase (HIDH). Furthermore, GSO-specific tandem gene duplication resulted in the gene expansion of C-glucosyltransferase (CGT). Integrative analysis of the metabolome and transcriptome identified 13 CGTs and eight HIDHs involved in the biosynthetic pathway of schaftoside. Functional studies indicated that CGTs and HIDHs identified here are bona fide responsible for the biosynthesis of schaftoside in GSO, as confirmed through hairy root transgenic system and in vitro enzyme activity assay. Taken together, the ancestral papilionoid polyploidy event expanding CHSs and HIDHs, along with the GSO-specific tandem duplication of CGT, contributes, partially if not completely, to the robust biosynthesis of schaftoside in GSO. These findings provide insights into the genomic mechanisms underlying the abundant biosynthesis of schaftoside in GSO, highlighting the potential of GSO as a source of bioactive compounds for pharmaceutical development.

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

Chromosome-level genome assembly of Prunella vulgaris L. provides insights into pentacyclic triterpenoid biosynthesis

Prunella vulgaris is one of the bestselling and widely used medicinal herbs. It is recorded as an ace medicine for cleansing and protecting the liver in Chinese Pharmacopoeia and has been used as the main constitutions of many herbal tea formulas in China for centuries. It is also a traditional folk medicine in Europe and other countries of Asia. Pentacyclic triterpenoids are a major class of bioactive compounds produced in P. vulgaris. However, their biosynthetic mechanism remains to be elucidated. Here, we report a chromosome-level reference genome of P. vulgaris using an approach combining Illumina, ONT, and Hi-C technologies. It is 671.95 Mb in size with a scaffold N50 of 49.10 Mb and a complete BUSCO of 98.45%. About 98.31% of the sequence was anchored into 14 pseudochromosomes. Comparative genome analysis revealed a recent WGD in P. vulgaris. Genome-wide analysis identified 35 932 protein-coding genes (PCGs), of which 59 encode enzymes involved in 2,3-oxidosqualene biosynthesis. In addition, 10 PvOSC, 358 PvCYP, and 177 PvUGT genes were identified, of which five PvOSCs, 25 PvCYPs, and 9 PvUGTs were predicted to be involved in the biosynthesis of pentacyclic triterpenoids. Biochemical activity assay of PvOSC2, PvOSC4, and PvOSC6 recombinant proteins showed that they were mixed amyrin synthase (MAS), lupeol synthase (LUS), and β-amyrin synthase (BAS), respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of pentacyclic triterpenoids in P. vulgaris.

A comprehensive analysis integrating phenotypic assessment uncovering thornless cultivar lineages in Aralia elata

Aralia elata is an Araliaceae woody plant species found in Northeastern Asia. To understand how genetic pools are distributed for A.elata clones, we were to analyze the population structure of A.elata cultivars and identify how these are correlated with thorn-related phenotype which determines the utility of A.elata. We found that the de novo assembled genome of 'Yeongchun' shared major genomic compartments with the public A.elata genome assembled from the wild-type from China. To identify the population structure of the 32 Korean and Japanese cultivars, we identified 44 SSR markers and revealed three main sub-clusters using ΔK analysis with one isolated cultivar. Machine-learning based clustering with thorn-related phenotype correlated moderately with population structure based on SSR analysis suggested multi-layered genetic regulation of thorn-related phenotypes. Thus, we revealed genetic lineage of A.elata and uncovered isolated cultivar which can provide new genetic material for further breeding.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Functional Identification of HhUGT74AG11-A Key Glycosyltransferase Involved in Biosynthesis of Oleanane-Type Saponins in Hedera helix

Hedera helix is a traditional medicinal plant. Its primary active ingredients are oleanane-type saponins, which have extensive pharmacological effects such as gastric mucosal protection, autophagy regulation actions, and antiviral properties. However, the glycosylation-modifying enzymes responsible for catalyzing oleanane-type saponin biosynthesis remain unidentified. Through transcriptome, cluster analysis, and PSPG structural domain, this study preliminarily screened four candidate UDP-glycosyltransferases (UGTs), including Unigene26859, Unigene31717, CL11391.Contig2, and CL144.Contig9. In in vitro enzymatic reactions, it has been observed that Unigene26859 (HhUGT74AG11) has the ability to facilitate the conversion of oleanolic acid, resulting in the production of oleanolic acid 28-O-glucopyranosyl ester. Moreover, HhUGT74AG11 exhibits extensive substrate hybridity and specific stereoselectivity and can transfer glycosyl donors to the C-28 site of various oleanane-type triterpenoids (hederagenin and calenduloside E) and the C-7 site of flavonoids (tectorigenin). Cluster analysis found that HhUGT74AG11 is clustered together with functionally identified genes AeUGT74AG6, CaUGT74AG2, and PgUGT74AE2, further verifying the possible reason for HhUGT74AG11 catalyzing substrate generalization. In this study, a novel glycosyltransferase, HhUGT74AG11, was characterized that plays a role in oleanane-type saponins biosynthesis in H. helix, providing a theoretical basis for the production of rare and valuable triterpenoid saponins.

A haplotype-resolved gap-free genome assembly provides novel insight into monoterpenoid diversification in Mentha suaveolens 'Variegata'

Mentha is a commonly used spice worldwide, which possesses medicinal properties and fragrance. These characteristics are conferred, at least partially, by essential oils such as menthol. In this study, a gap-free assembly with a genome size of 414.3 Mb and 31,251 coding genes was obtained for Mentha suaveolens 'Variegata'. Based on its high heterozygosity (1.5%), two complete haplotypic assemblies were resolved, with genome sizes of 401.9 and 405.7 Mb, respectively. The telomeres and centromeres of each haplotype were almost fully annotated. In addition, we detected a total of 41,135 structural variations. Enrichment analysis demonstrated that genes involved in terpenoid biosynthesis were affected by these structural variations. Analysis of volatile metabolites showed that M. suaveolens mainly produces piperitenone oxide rather than menthol. We identified three genes in the M. suaveolens genome which encode isopiperitenone reductase (ISPR), a key rate-limiting enzyme in menthol biosynthesis. However, the transcription levels of ISPR were low. Given that other terpenoid biosynthesis genes were expressed, M. suaveolens ISPRs may account for the accumulation of piperitenone oxide in this species. The findings of this study may provide a valuable resource for improving the detection rate and accuracy of genetic variants, thereby enhancing our understanding of their impact on gene function and expression. Moreover, our haplotype-resolved gap-free genome assembly offers novel insights into molecular marker-assisted breeding of Mentha.

Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

PgDDS Changes the Plant Growth of Transgenic Aralia elata and Improves the Production of Re and Rg3 in Its Leaves

Aralia elata (Miq.) Seem is a medicinal plant that shares a common pathway for the biosynthesis of triterpenoid saponins with Panax ginseng. Here, we transferred the dammarenediol-II synthase gene from P. ginseng (PgDDS; GenBank: AB122080.1) to A. elata. The growth of 2-year-old transgenic plants (L27; 9.63 cm) was significantly decreased compared with wild-type plants (WT; 74.97 cm), and the leaflet shapes and sizes of the transgenic plants differed from those of the WT plants. Based on a terpene metabolome analysis of leaf extracts from WT, L13, and L27 plants, a new structural skeleton for ursane-type triterpenoid saponins was identified. Six upregulated differentially accumulated metabolites (DAMs) were detected, and the average levels of Rg3 and Re in the leaves of the L27 plants were 42.64 and 386.81 μg/g, respectively, increased significantly compared with the WT plants (15.48 and 316.96 μg/g, respectively). Thus, the expression of PgDDS in A. elata improved its medicinal value.

The genome of the toxic invasive species Heracleum sosnowskyi carries an increased number of genes despite absence of recent whole-genome duplications

Heracleum sosnowskyi, belonging to a group of giant hogweeds, is a plant with large effects on ecosystems and human health. It is an invasive species that contributes to the deterioration of grassland ecosystems. The ability of H. sosnowskyi to produce linear furanocoumarins (FCs), photosensitizing compounds, makes it very dangerous. At the same time, linear FCs are compounds with high pharmaceutical value used in skin disease therapies. Despite this high importance, it has not been the focus of genetic and genomic studies. Here, we report a chromosome-scale assembly of Sosnowsky's hogweed genome. Genomic analysis revealed an unusually high number of genes (55106) in the hogweed genome, in contrast to the 25-35 thousand found in most plants. However, we did not find any traces of recent whole-genome duplications not shared with its confamiliar, Daucus carota (carrot), which has approximately thirty thousand genes. The analysis of the genomic proximity of duplicated genes indicates on tandem duplications as a main reason for this increase. We performed a genome-wide search of the genes of the FC biosynthesis pathway and surveyed their expression in aboveground plant parts. Using a combination of expression data and phylogenetic analysis, we found candidate genes for psoralen synthase and experimentally showed the activity of one of them using a heterologous yeast expression system. These findings expand our knowledge on the evolution of gene space in plants and lay a foundation for further analysis of hogweed as an invasive plant and as a source of FCs.

Identification of Functional Brassinosteroid Receptor Genes in Oaks and Functional Analysis of QmBRI1

Brassinosteroids (BRs) play important regulatory roles in plant growth and development, with functional BR receptors being crucial for BR recognition or signaling. Although functional BR receptors have been extensively studied in herbaceous plants, they remain largely under-studied in forest tree species. In this study, nine BR receptors were identified in three representative oak species, of which BRI1s and BRL1s were functional BR receptors. Dispersed duplications were a driving force for oak BR receptor expansion, among which the Brassinosteroid-Insensitive-1 (BRI1)-type genes diverged evolutionarily from most rosids. In oak BRI1s, we identified that methionine in the conserved Asn-Gly-Ser-Met (NGSM) motif was replaced by isoleucine and that the amino acid mutation occurred after the divergence of Quercus and Fagus. Compared with QmBRL1, QmBRI1 was relatively highly expressed during BR-induced xylem differentiation and in young leaves, shoots, and the phloem and xylem of young stems of Quercus mongolica. Based on Arabidopsis complementation experiments, we proved the important role of QmBRI1 in oak growth and development, especially in vascular patterning and xylem differentiation. These findings serve as an important supplement to the findings of the structural, functional and evolutionary studies on functional BR receptors in woody plants and provide the first example of natural mutation occurring in the conserved BR-binding region (NGSM motif) of angiosperm BRI1s.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Genomic and Transcriptomic Insights into the Evolution and Divergence of MIKC-Type MADS-Box Genes in Carica papaya

MIKC-type MADS-box genes, also known as type II genes, play a crucial role in regulating the formation of floral organs and reproductive development in plants. However, the genome-wide identification and characterization of type II genes as well as a transcriptomic survey of their potential roles in Carica papaya remain unresolved. Here, we identified and characterized 24 type II genes in the C. papaya genome, and investigated their evolutional scenario and potential roles with a widespread expression profile. The type II genes were divided into thirteen subclades, and gene loss events likely occurred in papaya, as evidenced by the contracted member size of most subclades. Gene duplication mainly contributed to MIKC-type gene formation in papaya, and the duplicated gene pairs displayed prevalent expression divergence, implying the evolutionary significance of gene duplication in shaping the diversity of type II genes in papaya. A large-scale transcriptome analysis of 152 samples indicated that different subclasses of these genes showed distinct expression patterns in various tissues, biotic stress response, and abiotic stress response, reflecting their divergent functions. The hub-network of male and female flowers and qRT-PCR suggested that TT16-3 and AGL8 participated in male flower development and seed germination. Overall, this study provides valuable insights into the evolution and functions of MIKC-type genes in C. papaya.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62 Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

Identification and Quantitation of the Bioactive Components in Wasted Aralia elata Leaves Extract with Endothelial Protective Activity

Aralia elata, a renowned medicinal plant with a rich history in traditional medicine, has gained attention for its potential therapeutic applications. However, the leaves of this plant have been largely overlooked and discarded due to limited knowledge of their biological activity and chemical composition. To bridge this gap, a comprehensive study was conducted to explore the therapeutic potential of the 70% ethanol extract derived from Aralia elata leaves (LAE) for the treatment of cardiovascular disease (CVD). Initially, the cytotoxic effects of LAE on human umbilical vein endothelial cells (HUVECs) were assessed, revealing no toxicity within concentrations up to 5 μg/mL. This suggests that LAE could serve as a safe raw material for the development of health supplements and drugs aimed at promoting cardiovascular well-being. Furthermore, the study found that LAE extract demonstrated anti-inflammatory properties in HUVECs by modulating the PI3K/Akt and MAPK signaling pathways. These findings are particularly significant as inflammation plays a crucial role in the progression of CVD. Moreover, LAE extract exhibited the ability to suppress the expression of adhesion molecules VCAM-1 and ICAM-1, which are pivotal in leukocyte migration to inflamed blood vessels observed in various pathological conditions. In conjunction with the investigation on therapeutic potential, the study also established an optimal HPLC-PDA-ESI-MS/MS method to identify and confirm the chemical constituents present in 24 samples collected from distinct regions in South Korea. Tentative identification revealed the presence of 14 saponins and nine phenolic compounds, while further analysis using PCA and PLS-DA allowed for the differentiation of samples based on their geographical origins. Notably, specific compounds such as chlorogenic acid, isochlorogenic acid A, and quercitrin emerged as marker compounds responsible for distinguishing samples from different regions. Overall, by unraveling its endothelial protective activity and identifying key chemical constituents, this research not only offers valuable insights for the development of novel treatments but also underscores the importance of utilizing and preserving natural resources efficiently.

Expression of PnSS Promotes Squalene and Oleanolic Acid (OA) Accumulation in Aralia elata via Methyl Jasmonate (MeJA) Induction

Aralia elata is an important herb due to the abundance of pentacyclic triterpenoid saponins whose important precursors are squalene and OA. Here, we found that MeJA treatment promoted both precursors accumulation, especially the latter, in transgenic A. elata, overexpressing a squalene synthase gene from Panax notoginseng(PnSS). In this study, Rhizobium-mediated transformation was used to express the PnSS gene. Gene expression analysis and high-performance liquid chromatography (HPLC) were used to identify the effect of MeJA on squalene and OA accumulation. The PnSS gene was isolated and expressed in A. elata. Transgenic lines showed a very high expression of the PnSS gene and farnesyl diphosphate synthase gene (AeFPS) and a slightly higher squalene content than the wild-type, but endogenous squalene synthase (AeSS), squalene epoxidase (AeSE), and β-amyrin synthase (Aeβ-AS) gene were decreased as well as OA content. Following one day of MeJA treatment, the expression levels of PeSS, AeSS, and AeSE genes increased significantly. On day 3, the maximum content of both products reached 17.34 and 0.70 mg·g^-1, which increased 1.39- and 4.90-fold than in the same lines without treatment. Transgenic lines expressing PnSS gene had a limited capability to promote squalene and OA accumulation. MeJA strongly activated their biosynthesis pathways, leading to enhance yield.

Novel insights into maize (Zea mays) development and organogenesis for agricultural optimization

In maize, intrinsic hormone activities and sap fluxes facilitate organogenesis patterning and plant holistic development; these hormone movements should be a primary focus of developmental biology and agricultural optimization strategies. Maize (Zea mays) is an important crop plant with distinctive life history characteristics and structural features. Genetic studies have extended our knowledge of maize developmental processes, genetics, and molecular ecophysiology. In this review, the classical life cycle and life history strategies of maize are analyzed to identify spatiotemporal organogenesis properties and develop a definitive understanding of maize development. The actions of genes and hormones involved in maize organogenesis and sex determination, along with potential molecular mechanisms, are investigated, with findings suggesting central roles of auxin and cytokinins in regulating maize holistic development. Furthermore, investigation of morphological and structural characteristics of maize, particularly node ubiquity and the alternate attachment pattern of lateral organs, yields a novel regulatory model suggesting that maize organ initiation and subsequent development are derived from the stimulation and interaction of auxin and cytokinin fluxes. Propositions that hormone activities and sap flow pathways control organogenesis are thoroughly explored, and initiation and development processes of distinctive maize organs are discussed. Analysis of physiological factors driving hormone and sap movement implicates cues of whole-plant activity for hormone and sap fluxes to stimulate maize inflorescence initiation and organ identity determination. The physiological origins and biogenetic mechanisms underlying maize floral sex determination occurring at the tassel and ear spikelet are thoroughly investigated. The comprehensive outline of maize development and morphogenetic physiology developed in this review will enable farmers to optimize field management and will provide a reference for de novo crop domestication and germplasm improvement using genome editing biotechnologies, promoting agricultural optimization.

Site-directed mutagenesis identified the key active site residues of 2,3-oxidosqualene cyclase HcOSC6 responsible for cucurbitacins biosynthesis in Hemsleya chinensis

Hemsleya chinensis is a Chinese traditional medicinal plant, containing cucurbitacin IIa (CuIIa) and cucurbitacin IIb (CuIIb), both of which have a wide range of pharmacological effects, including antiallergic, anti-inflammatory, and anticancer properties. However, few studies have been explored on the key enzymes that are involved in cucurbitacins biosynthesis in H. chinensis. Oxidosqualene cyclase (OSC) is a vital enzyme for cyclizing 2,3-oxidosqualene and its analogues. Here, a gene encoding the oxidosqualene cyclase of H. chinensis (HcOSC6), catalyzing to produce cucurbitadienol, was used as a template of mutagenesis. With the assistance of AlphaFold2 and molecular docking, we have proposed for the first time to our knowledge the 3D structure of HcOSC6 and its binding features to 2,3-oxidosqualene. Mutagenesis experiments on HcOSC6 generated seventeen different single-point mutants, showing that single-residue changes could affect its activity. Three key amino acid residues of HcOSC6, E246, M261 and D490, were identified as a prominent role in controlling cyclization ability. Our findings not only comprehensively characterize three key residues that are potentially useful for producing cucurbitacins, but also provide insights into the significant role they could play in metabolic engineering.

The Identification of SQS/SQE/OSC Gene Families in Regulating the Biosynthesis of Triterpenes in Potentilla anserina

The tuberous roots of Potentilla anserina (Pan) are an edible and medicinal resource in Qinghai-Tibetan Plateau, China. The triterpenoids from tuberous roots have shown promising anti-cancer, hepatoprotective, and anti-inflammatory properties. In this study, we carried out phylogenetic analysis of squalene synthases (SQSs), squalene epoxidases (SQEs), and oxidosqualene cyclases (OSCs) in the pathway of triterpenes. In total, 6, 26, and 20 genes of SQSs, SQEs, and OSCs were retrieved from the genome of Pan, respectively. Moreover, 6 SQSs and 25 SQEs genes expressed in two sub-genomes (A and B) of Pan. SQSs were not expanded after whole-genome duplication (WGD), and the duplicated genes were detected in SQEs. Twenty OSCs were divided into two clades of cycloartenol synthases (CASs) and β-amyrin synthases (β-ASs) by a phylogenetic tree, characterized with gene duplication and evolutionary divergence. We speculated that β-ASs and CASs may participate in triterpenes synthesis. The data presented act as valuable references for future studies on the triterpene synthetic pathway of Pan.

The triterpenoid saponin content difference is associated with the two type oxidosqualene cyclase gene copy numbers of Pulsatilla chinensis and Pulsatilla cernua

Pulsatilla chinensis is an important medicinal herb, its dried radix is used to treat the inflammation since ancient China. Triterpenoid saponins are proved to be the main active compounds of Pulsatilla genus. The triterpenoid saponin contents vary widely in different Pulsatilla species. But no enzyme involved in the triterpenoid saponin biosynthetic pathway was identified in Pulsitilla genus. This seriously limits the explanation of the triterpene content difference of Pulsatilla species. In this article, we obtained two oxidosqualene cyclase (OSC) genes from P. chinensis and P. cernua by touchdown PCR and anchored PCR. These two OSCs converted 2,3-oxidosqualene into different triterpenoids. The OSC from P. cernua is a monofunctional enzyme for β-amyrin synthesis, while the OSC from P. chinensis is a multifunctional enzyme for lupeol and β-amyrin synthesis, and the lupeol is the main product. Then we identified the 260th amino acid residue was the key site for the product difference by gene fusion and site-directed mutant technology. When the 260th amino acid residue was tryptophan (W260) and phenylalanine (F260), the main catalysate was β-amyrin and lupeol, respectively. Then we found that the expression of these two genes was strongly correlated with the lupeol-type and β-amyrin-type triterpenoid contents in P. cernua and P. chinensis. Finally, we found the gene copy number difference of these two genotypes leaded to the triterpenoid diversity in P. cernua and P. chinensis. This study provides useful information for the molecular breeding and quality improvement of P. chinensis and a molecular marker to identify the P. chinensis decoction pieces.

Genome-wide identification of MBD gene family members in Eleutherococcus senticosus with their expression motifs under drought stress and DNA demethylation

BACKGROUND

Methyl-binding domain (MBD) is a class of methyl-CpG-binding domain proteins that affects the regulation of gene expression through epigenetic modifications. MBD genes are not only inseparable from DNA methylation but have also been identified and validated in various plants. Although MBD is involved in a group of physiological processes and stress regulation in these plants, MBD genes in Eleutherococcus senticosus remain largely unknown.

RESULTS

Twenty EsMBD genes were identified in E. senticosus. Among the 24 chromosomes of E. senticosus, EsMBD genes were unevenly distributed on 12 chromosomes, and only one tandem repeat gene existed. Collinearity analysis showed that the fragment duplication was the main motif for EsMBD gene expansion. As the species of Araliaceae evolved, MBD genes also evolved and gradually exhibited different functional differentiation. Furthermore, cis-acting element analysis showed that there were numerous cis-acting elements in the EsMBD promoter region, among which light response elements and anaerobic induction elements were dominant. The expression motif analysis revealed that 60% of the EsMBDs were up-regulated in the 30% water content group.

CONCLUSIONS

By comparing the transcriptome data of different saponin contents of E. senticosus and integrating them with the outcomes of molecular docking analysis, we hypothesized that EsMBD2 and EsMBD5 jointly affect the secondary metabolic processes of E. senticosus saponins by binding to methylated CpG under conditions of drought stress. The results of this study laid the foundation for subsequent research on the E. senticosus and MBD genes.

Genome-wide identification and functional analysis of the TCP gene family in rye (Secale cereale L.)

TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP) proteins are plant-specific transcription factors that play significant roles in plant growth, development, and stress response. Rye is a high-value crop with strong resistance to adverse environments. However, the functions of TCP proteins in rye are rarely reported. Based on a genome-wide analysis, the present study identified 26 TCP genes (ScTCPs) in rye. Mapping showed an uneven distribution of the ScTCP genes on the seven rye chromosomes and detected three pairs of tandem duplication genes. Phylogenetic analysis divided these genes into PCF (Proliferrating Cell Factors), CIN (CINCINNATA), and CYC (CYCLOIDEA)/TB1 (Teosinte Branched1) classes, which showed the highest homology between rye and wheat genes. Analysis of miRNA targeting sites indicated that five ScTCP genes were identified as potential targets of miRNA319. Promoter cis-acting elements analysis indicated that ScTCPs were regulated by light signals. Further analysis of the gene expression patterns and functional annotations suggested the role of a few ScTCPs in grain development and stress response. In addition, two TB1 homologous genes (ScTCP9 and ScTCP10) were identified in the ScTCP family. Synteny analysis showed that TB1 orthologous gene pairs existed before the ancestral divergence. Finally, the yeast two-hybrid assay and luciferase complementation imaging assay proved that ScTCP9, localized in the nucleus, interacts with ScFT (Flowering locus T), indicating their role in regulating flowering time. Taken together, this comprehensive study of ScTCPs provides important information for further research on gene function and crop improvement.

Full-length transcriptome, proteomics and metabolite analysis reveal candidate genes involved triterpenoid saponin biosynthesis in Dipsacus asperoides

Dipsacus asperoides is a traditional medicinal herb widely used in inflammation and fracture in Asia. Triterpenoid saponins from D. asperoides are the main composition with pharmacological activity. However, the biosynthesis pathway of triterpenoid saponins has not been completely resolved in D. asperoides. Here, the types and contents of triterpenoid saponins were discovered with different distributions in five tissues (root, leaf, flower, stem, and fibrous root tissue) from D. asperoides by UPLC-Q-TOF-MS analysis. The discrepancy between five tissues in D. asperoides at the transcriptional level was studied by combining single-molecule real-time sequencing and next- generation sequencing. Meanwhile, key genes involved in the biosynthesis of saponin were further verified by proteomics. In MEP and MVA pathways, 48 differentially expressed genes were identified through co-expression analysis of transcriptome and saponin contents, including two isopentenyl pyrophosphate isomerase and two 2,3-oxidosqualene β-amyrin cyclase, etc. In the analysis of WGCNA, 6 cytochrome P450s and 24 UDP- glycosyltransferases related to the biosynthesis of triterpenoid saponins were discovered with high transcriptome expression. This study will provide profound insights to demonstrate essential genes in the biosynthesis pathway of saponins in D. asperoides and support for the biosynthetic of natural active ingredients in the future.

Microorganisms for Ginsenosides Biosynthesis: Recent Progress, Challenges, and Perspectives

Ginsenosides are major bioactive compounds present in the Panax species. Ginsenosides exhibit various pharmaceutical properties, including anticancer, anti-inflammatory, antimetastatic, hypertension, and neurodegenerative disorder activities. Although several commercial products have been presented on the market, most of the current chemical processes have an unfriendly environment and a high cost of downstream processing. Compared to plant extraction, microbial production exhibits high efficiency, high selectivity, and saves time for the manufacturing of industrial products. To reach the full potential of the pharmaceutical resource of ginsenoside, a suitable microorganism has been developed as a novel approach. In this review, cell biological mechanisms in anticancer activities and the present state of research on the production of ginsenosides are summarized. Microbial hosts, including native endophytes and engineered microbes, have been used as novel and promising approaches. Furthermore, the present challenges and perspectives of using microbial hosts to produce ginsenosides have been discussed.

The chromosome-level genome of female ginseng (Angelica sinensis) provides insights into molecular mechanisms and evolution of coumarin biosynthesis

Coumarins are natural products with important medicinal values, and include simple coumarins, furanocoumarins and pyranocoumarins. Female ginseng (Angelica sinensis) is a renowned herb with abundant coumarins, originated in China and known for the treatment of female ailments for thousands of years. The molecular basis of simple coumarin biosynthesis in A. sinensis and the evolutionary history of the genes involved in furanocoumarin biosynthesis are largely unknown. Here, we generated the first chromosome-scale genome of A. sinensis. It has a genome size of 2.37 Gb, which was generated by combining PacBio and Hi-C sequencing technologies. The genome was predicted to contain 43 202 protein-coding genes dispersed mainly on 11 pseudochromosomes. We not only provided evidence for whole-genome duplication (WGD) specifically occurring in the Apioideae subfamily, but also demonstrated the vital role of tandem duplication for phenylpropanoid biosynthesis in A. sinensis. Combined analyses of transcriptomic and metabolomic data revealed key genes and candidate transcription factors regulating simple coumarin biosynthesis. Furthermore, phylogenomic synteny network analyses suggested prenyltransferase genes involved in furanocoumarin biosynthesis evolved independently in the Moraceae, Fabaceae, Rutaceae and Apiaceae after ζ and ε WGD. Our work sheds light on coumarin biosynthesis, and provides a benchmark for accelerating genetic research and molecular breeding in A. sinensis.

Applications of protein engineering in the microbial synthesis of plant triterpenoids

Triterpenoids are a class of natural products widely used in fields related to medicine and health due to their biological activities such as hepatoprotection, anti-inflammation, anti-viral, and anti-tumor. With the advancement in biotechnology, microorganisms have been used as cell factories to produce diverse natural products. Despite the significant progress that has been made in the construction of microbial cell factories for the heterogeneous biosynthesis of triterpenoids, the industrial production of triterpenoids employing microorganisms has been stymied due to the shortage of efficient enzymes as well as the low expression and low catalytic activity of heterologous proteins in microbes. Protein engineering has been demonstrated as an effective way for improving the specificity, catalytic activity, and stability of the enzyme, which can be employed to overcome these challenges. This review summarizes the current progress in the studies of Oxidosqualene cyclases (OSCs), cytochrome P450s (P450s), and UDP-glycosyltransferases (UGTs), the key enzymes in the triterpenoids synthetic pathway. The main obstacles restricting the efficient catalysis of these key enzymes are analyzed, the applications of protein engineering for the three key enzymes in the microbial synthesis of triterpenoids are systematically reviewed, and the challenges and prospects of protein engineering are also discussed.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Comprehensive identification of bHLH transcription factors in Litsea cubeba reveals candidate gene involved in the monoterpene biosynthesis pathway

Litsea cubeba (Lour.) Person, an economically important aromatic plant producing essential oils, has lemon-like fragrance and 96.44-98.44% monoterpene contents. bHLH transcription factor plays an important role in plant secondary metabolism and terpene biosynthesis. In this study, we used bioinformatics to identify bHLH transcription factors in L. cubeba, 173 bHLH genes were identified from L. cubeba and divided these into 26 subfamilies based on phylogenetic analysis. The majority of bHLHs in each subfamily shared comparable structures and motifs. While LcbHLHs were unevenly distributed across 12 chromosomes, 10 tandem repeats were discovered. Expression profiles of bHLH genes in different tissues demonstrated that LcbHLH78 is a potential candidate gene for regulating monoterpene biosynthesis. LcbHLH78 and the terpene synthase LcTPS42 showed comparable expression patterns in various tissues and fruit development stages of L. cubeba. Subcellular localization analysis revealed that LcbHLH78 protein localizes to the nucleus, consistent with a transcription factor function. Importantly, transient overexpression of LcbHLH78 increased geraniol and linalol contents. Our research demonstrates that LcbHLH78 enhances terpenoid biosynthesis. This finding will be beneficial for improving the quality of L. cubeba and provides helpful insights for further research into the control mechanism of LcbHLH genes over terpenoid biosynthesis.