Characterization of the Cytochrome P450 CYP716C52 in Celastrol Biosynthesis and Its Applications in Engineered Saccharomyces cerevisiae.

Celastrol is a bioactive pentacyclic triterpenoid with promising therapeutic effects that is mainly distributed in Celastraceae plants. Although some enzymes involved in the celastrol biosynthesis pathway have been reported, many biosynthetic steps remain unknown. Herein, transcriptomics and metabolic profiles of multiple species in Celastraceae were integrated to screen for cytochrome P450s (CYPs) that are closely related to celastrol biosynthesis. The CYP716 enzyme, TwCYP716C52, was found to be able to oxidize the C-2 position of polpunonic acid, a precursor of celastrol, to form the wilforic acid C. RNAi-mediated repression of TwCYP716C52 in Tripterygium wilfordii suspension cells further confirmed its involvement in celastrol biosynthesis. The C-2 catalytic mechanisms of TwCYP716C52 were further explored by using molecular docking and site-directed mutagenesis experiments. Moreover, a modular optimization strategy was used to construct an engineered yeast to produce wilforic acid C at a titer of 5.8 mg·L-1. This study elucidates the celastrol biosynthetic pathway and provides important functional genes and sufficient precursors for further enzyme discovery.

Probing the functions of friedelane-type triterpene cyclases from four celastrol-producing plants.

Triterpenes are among the most diverse plant natural products, and their diversity is closely related to various triterpene skeletons catalyzed by different 2,3-oxidosqualene cyclases (OSCs). Celastrol, a friedelane-type triterpene with significant bioactivities, is specifically distributed in higher plants, such as Celastraceae species. Friedelin is an important precursor for the biosynthesis of celastrol, and it is synthesized through the cyclization of 2,3-oxidosqualene, with the highest number of rearrangements being catalyzed by friedelane-type triterpene cyclases. However, the molecular mechanisms underlying the catalysis of friedelin production by friedelane-type triterpene cyclases have not yet been fully elucidated. In this study, transcriptome data of four celastrol-producing plants from Celastraceae were used to identify a total of 21 putative OSCs. Through functional characterization, the friedelane-type triterpene cyclases were separately verified in the four plants. Analysis of the selection pressure showed that purifying selection acted on these OSCs, and the friedelane-type triterpene cyclases may undergo weaker selective restriction during evolution. Molecular docking and site-directed mutagenesis revealed that changes in some amino acids that are unique to friedelane-type triterpene cyclases may lead to variations in catalytic specificity or efficiency, thereby affecting the synthesis of friedelin. Our research explored the functional diversity of triterpene synthases from a multispecies perspective. It also provides some references for further research on the relative mechanisms of friedelin biosynthesis.

CYP72D19 from Tripterygium wilfordii catalyzes C-2 hydroxylation of abietane-type diterpenoids.

KEY MESSAGE: CYP72D19, the first functional gene of the CYP72D subfamily, catalyzes the C-2 hydroxylation of abietane-type diterpenoids. The abietane-type diterpenoids, e.g., triptolide, tripdiolide, and 2-epitripdiolide, are the main natural products for the anti-tumor, anti-inflammatory, and immunosuppressive activities of Tripterygium wilfordii, while their biosynthetic pathways are not resolved. Here, we cloned and characterized the CYP72D19-catalyzed C-2 hydroxylation of dehydroabietic acid, a compound that has been proven to be a biosynthetic intermediate in triptolide biosynthesis. Through molecular docking and site-directed mutagenesis, L386, L387, and I493 near the active pocket were found to have an important effect on the enzyme activity, which also indicates that steric hindrance of residues plays an important role in function. In addition, CYP72D19 also catalyzed a variety of abietane-type diterpenoids with benzene ring, presumably because the benzene ring of the substrate molecule stabilized the C-ring, allowing the protein and the substrate to form a relatively stable spatial structure. This is the first demonstration of CYP72D subfamily gene function. Our research provides important genetic elements for the structural modification of active ingredients and the heterologous production of other 2-hydroxyl abietane-type natural products.

Structural and Catalytic Insight into the Unique Pentacyclic Triterpene Synthase TwOSC.

The oxidosqualene cyclase (OSC) catalyzed cyclization of the linear substrate (3S)-2,3-oxidosqualene to form diverse pentacyclic triterpenoid (PT) skeletons is one of the most complex reactions in nature. Friedelin has a unique PT skeleton involving a fascinating nine-step cation shuttle run (CSR) cascade rearrangement reaction, in which the carbocation formed at C2 moves to the other side of the skeleton, runs back to C3 to yield a friedelin cation, which is finally deprotonated. However, as crystal structure data of plant OSCs are lacking, it remains unknown why the CSR cascade reactions occur in friedelin biosynthesis, as does the exact catalytic mechanism of the CSR. In this study, we determined the first cryogenic electron microscopy structure of a plant OSC, friedelin synthase, from Tripterygium wilfordii Hook. f (TwOSC). We also performed quantum mechanics/molecular mechanics simulations to reveal the energy profile for the CSR cascade reaction and identify key residues crucial for PT skeleton formation. Furthermore, we semirationally designed two TwOSC mutants, which significantly improved the yields of friedelin and ß-amyrin, respectively.

[Seed morphology and germination characteristics of Tripterygium wilfordii].

The seeds of Tripterygium wilfordii are characterized by dormancy and a long germination cycle under natural sowing conditions. In this study, we developed a method for rapid germination of T. wilfordii seeds by analyzing the size, morphology, thousand-grain weight, viability, moisture content, physicochemical properties, and seed germination rates under different germination conditions. The seeds of T. wilfordii were fine columnar with a thick and hard outer seed coat. They had the length of 6.69 mm, the width of 2.14 mm, the thickness of 1.68 mm, the thousand-grain weight of 8.99 g, the moisture content of 8.86%, the soluble sugar content of 21.3 mg·g~(-1), the starch content of 28.9 mg·g~(-1), the soluble protein content of 44.2 mg·g~(-1), and the seed viability of only 54.0%. The seeds were respectively treated with distilled water, ultrasonication, low-temperature storage, 50 ℃ water, 100 mg·L~(-1) 6-BA, 0.6% KMnO_4, 1% KNO_3, 50 mg·L~(-1) NAA, and 100 mg·L~(-1) GA_3 solution. The results showed that soaking the seeds in 100 mg·L~(-1) GA_3 solution significantly promoted the germination. Further, the seeds were soaked in 50, 100, 250, 500, and 1 000 mg·L~(-1) GA_3 solutions, which demonstrated that high concentration(500 mg·L~(-1), 1 000 mg·L~(-1)) of GA_3 solutions increased the germination rate and speed and shortened the germination cycle from more than 3 months to less than 15 days. The findings of this study are of great significance to the breeding of T. wilfordii and lay a foundation for the large-scale propagation of T. wilfordii seeds and the excavation of T. wilfordii germplasm resources.

A cytochrome P450 CYP81AM1 from Tripterygium wilfordii catalyses the C-15 hydroxylation of dehydroabietic acid.

MAIN CONCLUSION: A novel cytochrome P450 from Tripterygium wilfordii, CYP81AM1, specifically catalyses the C-15 hydroxylation of dehydroabietic acid. This is the first CYP450 to be found in plants with this function. Cytochrome P450 oxygenases (CYPs) play an important role in the post-modification in biosynthesis of plant bioactive terpenoids. Here, we found that CYP81AM1 can catalyze the formation of 15-hydroxydehydroabietic acid by in vitro enzymatic reactions and in vivo yeast feeding assays. This is the first study to show that CYP81 family enzymes are involved in the hydroxylation of abietane diterpenoids. At the same time, we found that CYP81AM1 could not catalyse abietatriene and dehydroabietinol, suggesting that the occurrence of the reaction may be related to the carboxyl group. Through molecular docking and site mutations, it was found that some amino acid sites (F104, K107) near the carboxyl group had an important effect on enzyme activity, also suggesting that the carboxyl group played an important role in the occurrence of the reaction.

Engineering chimeric diterpene synthases and isoprenoid biosynthetic pathways enables high-level production of miltiradiene in yeast.

Miltiradiene is a key intermediate in the biosynthesis of many important natural diterpene compounds with significant pharmacological activity, including triptolide, tanshinones, carnosic acid and carnosol. Sufficient accumulation of miltiradiene is vital for the production of these medicinal compounds. In this study, comprehensive engineering strategies were applied to construct a high-yielding miltiradiene producing yeast strain. First, a chassis strain that can accumulate 2.1 g L-1 geranylgeraniol was constructed. Then, diterpene synthases from various species were evaluated for their ability to produce miltiradiene, and a chimeric miltiradiene synthase, consisting of class II diterpene synthase (di-TPS) CfTPS1 from Coleus forskohlii (Plectranthus barbatus) and class I di-TPS SmKSL1 from Salvia miltiorrhiza showed the highest efficiency in the conversion of GGPP to miltiradiene in yeast. Moreover, the miltiradiene yield was further improved by protein modification, which resulted in a final yield of 550.7 mg L-1 in shake flasks and 3.5 g L-1 in a 5-L bioreactor. This work offers an efficient and green process for the production of the important intermediate miltiradiene, and lays a foundation for further pathway reconstruction and the biotechnological production of valuable natural diterpenes.

Identification and functional characterization of squalene epoxidases and oxidosqualene cyclases from Tripterygium wilfordii.

KEY MESSAGE: We cloned two squalene epoxidases and five oxidosqualene cyclases, and identified their function using CRISPR/Cas9 tool and yeast heterologous expression. Triterpenes are the main active ingredients of Tripterygium wilfordii Hook.f., a traditional Chinese medicinal plant with many encouraging preclinical applications. However, the biosynthetic pathways of triterpenes in this plant are poorly understood. Here, we report on the isolation and identification of two squalene epoxidases (SQE6 and SQE7) and five oxidosqualene cyclases (OSC4-8) from T. wilfordii. Yeast complementation assays showed that TwSQE6 and TwSQE7 can functionally complement an erg1 yeast mutant that was constructed using the CRISPR/Cas9 system. The putative OSC genes were functionally characterized by heterologous expression in yeast. GC/MS analysis of the fermentation products of the transgenic yeast showed that both TwOSC4 and TwOSC6 are cycloartenol synthases, while TwOSC8 is a ß-amyrin synthase. The discovery of these genes expands our knowledge of key enzymes in triterpenoid biosynthesis, and provides additional target genes for increasing the production of triterpenes in T. wilfordii tissue cultures by disrupting competing pathways, or in chassis cells by reconstituting the triterpenoid biosynthetic pathway.

Isolation and characterization of a glycosyltransferase with specific catalytic activity towards flavonoids from Tripterygium wilfordii.

Flavonoids are important secondary metabolites that exist in many medicinal plants. Flavonoid glycosyltransferases can transfer sugar moieties to their parent rings, producing various flavonoid glycosides with significant pharmacological activities. Here, we report the molecular cloning of the O-glycosyltransferase TwUGT2 from Tripterygium wilfordii and its catalytic activity was explored by heterologous expression in E. coli. The results showed that TwUGT2 has specific glycosyltransferase activity towards C-3 and 7 hydroxyl groups of flavonoids, thereby converting quercetin and pinocembrin into isoquercitrin and pinocembrin 7-O-beta-D-glucoside, respectively. The identification of TwUGT2 will provide a useful molecular tool for synthetic biology and contribute to drug discovery.[Formula: see text].

Friedelane-type triterpene cyclase in celastrol biosynthesis from Tripterygium wilfordii and its application for triterpenes biosynthesis in yeast.

Celastrol is a promising bioactive compound isolated from Tripterygium wilfordii and has been shown to possess many encouraging preclinical applications. However, the celastrol biosynthetic pathway is poorly understood, especially the key oxidosqualene cyclase (OSC) enzyme responsible for cyclisation of the main scaffold. Here, we report on the isolation and characterisation of three OSCs from T. wilfordii: TwOSC1, TwOSC2 and TwOSC3. Both TwOSC1 and TwOSC3 were multiproduct friedelin synthases, while TwOSC2 was a ß-amyrin synthase. We further found that TwOSC1 and TwOSC3 were involved in the biosynthesis of celastrol and that their common product, friedelin, was a precursor of celastrol. We then reconstituted the biosynthetic pathway of friedelin in engineered yeast constructed by the CRISPR/Cas9 system, with protein modification and medium optimisation, leading to heterologous production of friedelin at 37.07 mg l-1 in a shake flask culture. Our study was the first to identify the genes responsible for biosynthesis of the main scaffold of celastrol and other triterpenes in T. wilfordii. As friedelin has been found in many plants, the results and approaches described here have laid a solid foundation for further explaining the biosynthesis of celastrol and related triterpenoids. Moreover, our results provide insights for metabolic engineering of friedelane-type triterpenes.

Probing the function of protein farnesyltransferase in Tripterygium wilfordii.

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

[Effect of GR24 on accumulation of diterpenoids in Tripterygium wilfordii suspension cells].

Terpenoids are main bioactive components in Tripterygium wilfordii,but the contents of some terpenoids are relatively low. In order to provide scientific evidence for the regulation of terpenoids in T. wilfordii,this research explored the effect of GR24 on accumulations of four diterpenoids( triptolide,tripterifordin,triptophenolide,and triptinin B) in T. wilfordii suspension cells by biological technology and UPLC-QQQ-MS/MS. The results indicated that 100 µmol·L-1 GR24 inhibited the accumulations of triptolide,tripterifordin,triptophenolide,and triptinin B to different degrees. Compared with the control group,the contents of 4 diterpenoids( in the induced group) were down to 96.59%,63.80%,61.02% and 33.59% in 240 h,respectively. Among them,the accumulation of triptinin B iswas significantly inhibited. In addition,the key time point of inhibitory effect was 120 h after induction with GR24 in some diterpenoids. This is the first systematic study focusing on the effect of GR24 on the accumulations of diterpenoids in T. wilfordii suspension cells. The dynamic accumulation ruleregularity of four diterpenoids after induced by GR24 was summarized,which laid a foundation for further study on the chemical response mechanism of terpenoids to GR24.

[Molecular authentication of Jinyinhua formula granule by using allele-specific PCR].

Traditional authentication method is hard to identify herb's authenticity of traditional Chinese medicine(TCM) formula granules because they have lost all their morphological characteristics. In this study, a new allele-specific PCR method was established for identifying the authentication of Jinyinhua formula granule (made from Lonicerae Japonicae Flos) based on an SNP site in trnL-trnF fragment. Genomic DNA was successfully extracted from Lonicerae Japonicae Flos and its formula granules by using an improved spin column method and then PCR was performed with the designed primer. Approximately 110 bp specific bands was obtained only in the authentic Lonicerae Japonicae Flos and its formula granules, while no bands were found in fake mixed products. In addition, the PCR product sequence was proved from Lonicerae Japonicae Flos trnL-trnF sequence by using BLAST method. Therefore, DNA molecular authentication method could make up the limitations of character identification method and microscopic identification, and quickly identify herb's authenticity of TCM formula granules, with enormous potential for market supervision and quality control.

[Cloning and protein expression analysis of geranyl diphosphate synthase genes in Tripterygium wilfordii].

Based on the transcriptome data, the study cloned full-length cDNA of TwGPPS1 and TwGPPS2 genes from Tripterygium wilfordii suspension cells and then analyzed the bioinformation of the sequence and protein expression. The cloned TwGPPS1 has a 1 278 bp open reading frame (ORF) encoding a polypeptide of 425 amino acids. The deduced isoelectric point (pI) was 6.68, a calculated molecular weight was about 47.189 kDa. The full-length cDNA of the TwGPPS2 contains a 1 269 bp open reading frame (ORF) encoding a polypeptide of 422 amino acids. The deduced isoelectric point (pI) was 6.71, a calculated molecular weight was about 46.774 kDa.The entire reading frame of TwGPPS1,2 was cloned into the pET-32a(+) vector and expressed in E. coli BL21 (DE3) cells to obtain the TwGPPS protein, which laid a basis for further study on the regulation of terpenoid secondary metabolism and biological synthesis.

Tandemly duplicated CYP82Ds catalyze 14-hydroxylation in triptolide biosynthesis and precursor production in Saccharomyces cerevisiae.

Triptolide is a valuable multipotent antitumor diterpenoid in Tripterygium wilfordii, and its C-14 hydroxyl group is often selected for modification to enhance both the bioavailability and antitumor efficacy. However, the mechanism for 14-hydroxylation formation remains unknown. Here, we discover 133 kb of tandem duplicated CYP82Ds encoding 11 genes on chromosome 12 and characterize CYP82D274 and CYP82D263 as 14-hydroxylases that catalyze the metabolic grid in triptolide biosynthesis. The two CYP82Ds catalyze the aromatization of miltiradiene, which has been repeatedly reported to be a spontaneous process. In vivo assays and evaluations of the kinetic parameters of CYP82Ds indicate the most significant affinity to dehydroabietic acid among multiple intermediates. The precursor 14-hydroxy-dehydroabietic acid is successfully produced by engineered Saccharomyces cerevisiae. Our study provides genetic elements for further elucidation of the downstream biosynthetic pathways and heterologous production of triptolide and of the currently intractable biosynthesis of other 14-hydroxyl labdane-type secondary metabolites.

Genome-wide analysis of MYB family genes in Tripterygium wilfordii and their potential roles in terpenoid biosynthesis.

Terpenoids are a class of significant bioactive components in the woody vine of Tripterygium wilfordii. Previous studies have shown that MYB transcription factors play important roles in plant secondary metabolism, growth, and developmental processes. However, the MYB involved in terpenoid biosynthesis in Tripterygium wilfordii are unknown. To identify Tripterygium wilfordii MYB (TwMYB) genes that are involved in terpenoid biosynthesis, we conducted the genome-wide analysis of the TwMYB gene family. A total of 207 TwMYBs were identified including 84 1R-TwMYB, 117 R2R3-TwMYB, four 3R-TwMYB, and two 4R-TwMYB genes. The most abundant R2R3-TwMYBs together with their Arabidopsis homologs were categorized into 26 subgroups. Intraspecific collinearity analysis found that the 74.9% of the TwMYBs may be generated by segmental duplication events, and 36.7% of duplicated gene pairs were derived from the specific whole genome duplication (WGD) event in Tripterygium wilfordii. In addition, interspecies collinearity analysis found that 16 TwMYB genes formed homologous gene pairs with MYB genes in seven representative species, which indicated they may have a key role in evolution. Notably, we found that the TwMYB genes were differentially expressed in various tissues by expression pattern analysis. In order to further select the candidate genes related to terpenoid biosynthesis, the assay of Methyl jasmonate (MeJA) induction and analysis of phylogenetic tree was conducted. It was speculated that six candidate TwMYB genes (TwMYB33, TwMYB34, TwMYB45, TwMYB67, TwMYB102, and TwMYB103) are involved in regulating terpenoid biosynthesis. This study is the first systematic analysis of the TwMYB gene family and will lay a foundation for the functional characterization of TwMYB genes in the regulation of terpenoid biosynthesis.

Mechanistic analysis for the origin of diverse diterpenes in Tripterygium wilfordii.

Tripterygium wilfordii is a valuable medicinal plant rich in biologically active diterpenoids, but there are few studies on the origins of these diterpenoids in its secondary metabolism. Here, we identified three regions containing tandemly duplicated diterpene synthase genes on chromosomes (Chr) 17 and 21 of T. wilfordii and obtained 11 diterpene synthases with different functions. We further revealed that these diterpene synthases underwent duplication and rearrangement at approximately 2.3-23.7 million years ago (MYA) by whole-genome triplication (WGT), transposon mediation, and tandem duplication, followed by functional divergence. We first demonstrated that four key amino acids in the sequences of TwCPS3, TwCPS5, and TwCPS6 were altered during evolution, leading to their functional divergence and the formation of diterpene secondary metabolites. Then, we demonstrated that the functional divergence of three TwKSLs was driven by mutations in two key amino acids. Finally, we discovered the mechanisms of evolution and pseudogenization of miltiradiene synthases in T. wilfordii and elucidated that the new function in TwMS1/2 from the terpene synthase (TPS)-b subfamily was caused by progressive changes in multiple amino acids after the WGT event. Our results provide key evidence for the formation of diverse diterpenoids during the evolution of secondary metabolites in T. wilfordii.

Functional characterization and substrate promiscuity of sesquiterpene synthases from Tripterygium wilfordii.

Acyclic terpenes, commonly found in plants, are of high physiological importance and commercial value, and their diversity was controlled by different terpene synthases. During the screen of sesquiterpene synthases from Tripterygium wilfordii, we observed that Ses-TwTPS1-1 and Ses-TwTPS2 promiscuously accepted GPP, FPP, and GGPP to produce corresponding terpene alcohols (linalool/nerolidol/geranyllinalool). The Ses-TwTPS1-2, Ses-TwTPS3, and Ses-TwTPS4 also showed unusual substrate promiscuity by catalyzing GGPP or GPP in addition to FPP as substrate. Furthermore, key residues for the generation of diterpene product, (E, E)-geranyllinalool, were screened depending on mutagenesis studies. The functional analysis of Ses-TwTPS1-1:V199I and Ses-TwTPS1-2:I199V showed that Val in 199 site assisted the produce of diterpene product geranyllinalool by enzyme mutation studies, which indicated that subtle differences away from the active site could alter the product outcome. Moreover, an engineered sesquiterpene high-yielding yeast that produced 162 mg/L nerolidol in shake flask conditions was constructed to quickly identify the function of sesquiterpene synthases in vivo and develop potential applications in microbial fermentation. Our functional characterization of acyclic sesquiterpene synthases will give some insights into the substrate promiscuity of diverse acyclic terpene synthases and provide key residues for expanding the product portfolio.

Cytochrome P450 catalyses the 29-carboxyl group formation of celastrol.

Celastrol, a potent anticancer and anti-obesity drug, was first isolated from Tripterygium wilfordii Hook. f. and it is produced in small quantities in many members of the Celastraceae family. The heterologous reconstitution of celastrol biosynthesis could be a promising method for the efficient production of celastrol and natural and unnatural derivatives thereof, yet only part of the biosynthetic pathway is known. Here, we report a cytochrome P450 monooxygenase (TwCYP712K1) from T. wilfordii that performs the three-step oxidation of friedelin to polpunonic acid in the celastrol pathway. Heterologous expression of TwCYP712K1 showed that TwCYP712K1 catalyses not only the transformation of friedelin to polpunonic acid but also the oxidation of ß-amyrin or α-amyrin. The role of TwCYP712K1 in the biosynthesis of celastrol was further revealed via RNA interference. Some key residues of TwCYP712K1 were also screened by molecular docking and site-directed mutagenesis. Our results lay a solid foundation for further elucidating the biosynthesis of celastrol and related triterpenoids.

The chromosome-level reference genome assembly for Panax notoginseng and insights into ginsenoside biosynthesis.

Panax notoginseng, a perennial herb of the genus Panax in the family Araliaceae, has played an important role in clinical treatment in China for thousands of years because of its extensive pharmacological effects. Here, we report a high-quality reference genome of P. notoginseng, with a genome size up to 2.66 Gb and a contig N50 of 1.12 Mb, produced with third-generation PacBio sequencing technology. This is the first chromosome-level genome assembly for the genus Panax. Through genome evolution analysis, we explored phylogenetic and whole-genome duplication events and examined their impact on saponin biosynthesis. We performed a detailed transcriptional analysis of P. notoginseng and explored gene-level mechanisms that regulate the formation of characteristic tubercles. Next, we studied the biosynthesis and regulation of saponins at temporal and spatial levels. We combined multi-omics data to identify genes that encode key enzymes in the P. notoginseng terpenoid biosynthetic pathway. Finally, we identified five glycosyltransferase genes whose products catalyzed the formation of different ginsenosides in P. notoginseng. The genetic information obtained in this study provides a resource for further exploration of the growth characteristics, cultivation, breeding, and saponin biosynthesis of P. notoginseng.