Cytochrome P450 Monooxygenase CYP716A141 is a Unique β-Amyrin C-16β Oxidase Involved in Triterpenoid Saponin Biosynthesis in Platycodon grandiflorus

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.

Identification and functional characterization of two trans-isopentenyl diphosphate synthases and one squalene synthase involved in triterpenoid biosynthesis in Platycodon grandiflorus

MAIN CONCLUSION

Two trans-isopentenyl diphosphate synthase and one squalene synthase genes were identified and proved to be involved in the triterpenoid biosynthesis in Platycodon grandiflorus. Platycodon grandiflorus is a commonly used traditional Chinese medicine. The main bioactive compounds of P. grandiflorus are triterpenoid saponins. The biosynthetic pathway of triterpenoid saponins in P. grandiflorus has been preliminarily explored. However, limited functional information on related genes has been reported. A total of three trans-isopentenyl diphosphate synthases (trans-IDSs) genes (PgFPPS, PgGGPPS1 and PgGGPPS2) and one squalene synthase (SQS) gene (PgSQS) in P. grandiflorus were screened and identified from transcriptome dataset. Subcellular localization of the proteins was defined based on the analysis of GFP-tagged. The activity of genes was verified in Escherichia coli, demonstrating that recombinant PgFPPS catalysed the production of farnesyl diphosphate. PgGGPPS1 produced geranylgeranyl diphosphate, whereas PgGGPPS2 did not exhibit catalytic activity. By structural identification of encoding genes, a transmembrane region was found at the C-terminus of the PgSQS gene, which produced an insoluble protein when expressed in E. coli but showed no apparent effect on the enzyme function. Furthermore, some triterpenoid saponin synthesis-related genes were discovered by combining the component content and the gene expression assays at the five growth stages of P. grandiflorus seedlings. The accumulation of active components in P. grandiflorus was closely associated with the expression level of genes related to the synthesis pathway.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of β-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

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.

Heterologous Expression of Platycodon grandiflorus PgF3'5'H Modifies Flower Color Pigmentation in Tobacco

Flavonoid-3',5'-hydroxylase (F3'5'H) is the key enzyme for the biosynthesis of delphinidin-based anthocyanins, which are generally required for purple or blue flowers. Previously, we isolated a full-length cDNA of PgF3'5'H from Platycodon grandiflorus, which shared the highest homology with Campanula medium F3'5'H. In this study, PgF3'5'H was subcloned into a plant over-expression vector and transformed into tobacco via Agrobacterium tumefaciens to investigate its catalytic function. Positive transgenic tobacco T0 plants were obtained by hygromycin resistance screening and PCR detection. PgF3'5'H showed a higher expression level in all PgF3'5'H transgenic tobacco plants than in control plants. Under the drive of the cauliflower mosaic virus (CaMV) 35S promoter, the over-expressed PgF3'5'H produced dihydromyricetin (DHM) and some new anthocyanin pigments (including delphinidin, petunidin, peonidin, and malvidin derivatives), and increased dihydrokaempferol (DHK), taxifolin, tridactyl, cyanidin derivatives, and pelargonidin derivatives in PgF3'5'H transgenic tobacco plants by ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) analysis, resulting in a dramatic color alteration from light pink to magenta. These results indicate that PgF3'5'H products have F3'5'H enzyme activity. In addition, PgF3'5'H transfer alters flavonoid pigment synthesis and accumulation in tobacco. Thus, PgF3'5'H may be considered a candidate gene for gene engineering to enhance anthocyanin accumulation and the molecular breeding project for blue flowers.

Plant Specialised Glycosides (PSGs): their biosynthetic enzymatic machinery, physiological functions and commercial potential

Plants, the primary producers of our planet, have evolved from simple aquatic life to very complex terrestrial habitat. This habitat transition coincides with evolution of enormous chemical diversity, collectively termed as 'Plant Specialised Metabolisms (PSMs)', to cope the environmental challenges. Plant glycosylation is an important process of metabolic diversification of PSMs to govern their in planta stability, solubility and inter/intra-cellular transport. Although, individual category of PSMs (terpenoids, phenylpropanoids, flavonoids, saponins, alkaloids, phytohormones, glucosinolates and cyanogenic glycosides) have been well studied; nevertheless, deeper insights of physiological functioning and genomic aspects of plant glycosylation/deglycosylation processes including enzymatic machinery (CYPs, GTs, and GHs) and regulatory elements are still elusive. Therefore, this review discussed the paradigm shift on genomic background of enzymatic machinery, transporters and regulatory mechanism of 'Plant Specialised Glycosides (PSGs)'. Current efforts also update the fundamental understanding about physiological, evolutionary and adaptive role of glycosylation/deglycosylation processes during the metabolic diversification of PSGs. Additionally, futuristic considerations and recommendations for employing integrated next-generation multi-omics (genomics, transcriptomics, proteomics and metabolomics), including gene/genome editing (CRISPR-Cas) approaches are also proposed to explore commercial potential of PSGs.

Cytochrome P450 monooxygenase-mediated tailoring of triterpenoids and steroids in plants

The cytochrome P450 monooxygenase (CYP) superfamily comprises hemethiolate enzymes that perform remarkable regio- and stereospecific oxidative chemistry. As such, CYPs are key agents for the structural and functional tailoring of triterpenoids, one of the largest classes of plant natural products with widespread applications in pharmaceuticals, food, cosmetics, and agricultural industries. In this review, we provide a full overview of 149 functionally characterised CYPs involved in the biosynthesis of triterpenoids and steroids in primary as well as in specialised metabolism. We describe the phylogenetic distribution of triterpenoid- and steroid-modifying CYPs across the plant CYPome, present a structure-based summary of their reactions, and highlight recent examples of particular interest to the field. Our review therefore provides a comprehensive up-to-date picture of CYPs involved in the biosynthesis of triterpenoids and steroids in plants as a starting point for future research.

Identification of key amino acid residues toward improving the catalytic activity and substrate specificity of plant-derived cytochrome P450 monooxygenases CYP716A subfamily enzyme for triterpenoid production in Saccharomyces cerevisiae

Triterpenoids constitute a group of specialized plant metabolites with wide structural diversity and high therapeutic value for human health. Cytochrome P450 monooxygenases (CYP) are a family of enzymes important for generating the structural diversity of triterpenoids by catalyzing the site-specific oxidization of the triterpene backbone. The CYP716 enzyme family has been isolated from various plant families as triterpenoid oxidases; however, their experimental crystal structures are not yet available and the detailed catalytic mechanism remains elusive. Here, we address this challenge by integrating bioinformatics approaches with data from other CYP families. Medicago truncatula CYP716A12, the first functionally characterized CYP716A subfamily enzyme, was chosen as the model for this study. We performed homology modeling, structural alignment, in silico site-directed mutagenesis, and molecular docking analysis to search and screen key amino acid residues relevant to the catalytic activity and substrate specificity of the CYP716A subfamily enzyme in triterpenoid biosynthesis. An in vivo functional analysis using engineered yeast that endogenously produced plant-derived triterpenes was performed to elucidate the results. When the amino acids in the signature region and substrate recognition sites (SRSs) were substituted, the product profile of CYP716A12 was modified. We identified amino acid residues that control the substrate contraction of the enzyme (D292) and engineered the enzyme to improve its catalytic activity and substrate specificity (D122, I212, and Q358) for triterpenoid biosynthesis. In addition, we demonstrated the versatility of this strategy by changing the properties of key residues in SRSs to improve the catalytic activity of Arabidopsis thaliana CYP716A1 (S356) and CYP716A2 (M206, F210) at C-28 on the triterpene backbone. This research has the potential to help in the production of desired triterpenoids in engineered yeast by increasing the catalytic activity and substrate specificity of plant CYP716A subfamily enzymes.

Whole-Genome Identification and Analysis of Multiple Gene Families Reveal Candidate Genes for Theasaponin Biosynthesis in Camellia oleifera

Camellia oleifera is an economically important oilseed tree. Seed meals of C. oleifera have a long history of use as biocontrol agents in shrimp farming and as cleaning agents in peoples’ daily lives due to the presence of theasaponins, the triterpene saponins from the genus Camellia. To characterize the biosynthetic pathway of theasaponins in C. oleifera, members of gene families involved in triterpenoid biosynthetic pathways were identified and subjected to phylogenetic analysis with corresponding members in Arabidopsis thaliana, Camellia sinensis, Actinidia chinensis, Panax ginseng, and Medicago truncatula. In total, 143 triterpenoid backbone biosynthetic genes, 1169 CYP450s, and 1019 UGTs were identified in C. oleifera. The expression profiles of triterpenoid backbone biosynthetic genes were analyzed in different tissue and seed developmental stages of C. oleifera. The results suggested that MVA is the main pathway for triterpenoid backbone biosynthesis. Moreover, the candidate genes for theasaponin biosynthesis were identified by WGCNA and qRT-PCR analysis; these included 11 CYP450s, 14 UGTs, and eight transcription factors. Our results provide valuable information for further research investigating the biosynthetic and regulatory network of theasaponins.

Rice (Oryza sativa L.) cytochrome P450 protein 716A subfamily CYP716A16 regulates disease resistance

Background

The sustainable development of rice production is facing severe threats by a variety of pathogens, such as necrotrophic Rhizoctonia solani and hemibiotrophic Xanthomonas oryzae pv. oryzae (Xoo). Mining and applying resistance genes to increase the durable resistance of rice is an effective method that can be used to control these diseases.

Results

In this research, we isolated and characterized CYP716A16, which is a positive regulator of rice to R. solani AG1-IA and Xoo, and belongs to the cytochrome P450 (CYP450) protein 716A subfamily. Overexpression (OE) of CYP716A16 resulted in enhanced resistance to R. solani AG1-IA and Xoo, while RNA interference (RNAi) of CYP716A16 resulted in increased susceptibility compared with wild-type (WT) plants. Additionally, jasmonic acid (JA)-dependent defense responses and reactive oxygen species (ROS) were activated in the CYP716A16-OE lines after R. solani AG1-IA inoculation. The comparative transcriptomic and metabolomics analysis of CYP716A16-OE and the WT lines showed that OE of CYP716A16 activated the biosynthesis of flavonoids and increased the amounts of narcissoside, methylophiopogonanone A, oroxin A, and amentoflavone in plants.

Conclusion

Based on these results, we suggest that JA-dependent response, ROS level, multiple resistance-related proteins, and flavonoid contents play an important role in CYP716A16-regulated R. solani AG1-IA and Xoo resistance. Our results broaden our knowledge regarding the function of a P450 protein 716A subfamily in disease resistance and provide new insight into the molecular mechanism of rice immune response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08568-8.

Cytochrome P450 Monooxygenase/Cytochrome P450 Reductase Bi-Enzymatic System Isolated From Ilex asprella for Regio-Specific Oxidation of Pentacyclic Triterpenoids

Ilex asprella is a plant from Aquifoliaceae. Its root is commonly used as folk medicinal materials in southern China. The chemical compositions of I. asprella are rich in pentacyclic triterpenoids, which show various biological activities and demonstrate a good prospect for drug development. The elucidation of biosynthesis mechanism of triterpenoids in I. asprella could lay important foundations for the production of these precious plant secondary metabolites by metabolic engineering. Our previous studies have revealed IaAO1 (a CYP716A210 homolog) responsible for the C-28 oxidation of α- and β-amyrin. Herein, we reported the identification of three more cytochrome P450 monooxygenase genes IaAO2 (a CYP716A212 homolog), IaAO4 (CYP714E88), IaAO5 (CYP93A220), and a cytochrome P450 reductase gene IaCPR by using Saccharomyces cerevisiae eukaryotic expression system and gas chromatography-mass spectrometry. Among them, the protein encoded by IaAO2 can catalyze the C-28 oxidation of α-amyrin and β-amyrin, IaAO4 can catalyze the C-23 oxidation of ursolic acid and oleanolic acid, while IaAO5 is responsible for the C-24 oxidation of β-amyrin. By introducing three genes IaAO1, IaAO4 and IaCPR into S. cerevisiae. We constructed an engineered yeast strain that can produce C-23 hydroxyl ursane-type triterpenoid derivatives. This study contributes to a thorough understanding of triterpenoid biosynthesis of medicinal plants and provides important tools for further metabolic engineering.

Identification and Functional Characterization of FLOWERING LOCUS T in Platycodon grandiflorus

Platycodon grandiflorus roots have been used as a foodstuff and traditional medicine for thousands of years in East Asia. In order to increase the root development of P. grandiflorus, cultivators removed the inflorescences, suggesting the possible negative effect of flowering on root development. This indicates that the genetic improvement of P. grandiflorus by late flowering is a potential approach to increase productivity. However, nothing is known about key genes integrating multiple flowering pathways in P. grandiflorus. In order to fill this gap, we identified potential homologs of the FLOWERING LOCUS T (FT) gene in P. grandiflorus. The alignment with other FT members and phylogenetic analysis revealed that the P. grandiflorus FT (PlgFT) protein contains highly conserved functional domains and belongs to the FT-like clade. The expression analysis revealed spatial variations in the transcription of PlgFT in different organs. In addition, the expression level of PlgFT was increased by high temperature but not by photoperiodic light input signals, presumably due to lacking the CONSTANS binding motif in its promoter region. Furthermore, PlgFT induced early flowering upon its overexpression in P. grandiflorus, suggesting the functional role of PlgFT in flowering. Taken together, we functionally characterized PlgFT as a master regulator of P. grandiflorus flowering under inductive high temperature, which will serve as an important target gene for improving the root productivity.

Transcriptome analysis identifies putative genes involved in triterpenoid biosynthesis in Platycodon grandiflorus

Comprehensive transcriptome analysis of different Platycodon grandiflorus tissues discovered genes related to triterpenoid saponin biosynthesis. Platycodon grandiflorus (Jacq.) A. DC. (P. grandiflorus), a traditional Chinese medicine, contains considerable triterpenoid saponins with broad pharmacological activities. Triterpenoid saponins are the major components of P. grandiflorus. Here, single-molecule real-time and next-generation sequencing technologies were combined to comprehensively analyse the transcriptome and identify genes involved in triterpenoid saponin biosynthesis in P. grandiflorus. We quantified four saponins in P. grandiflorus and found that their total content was highest in the roots and lowest in the stems and leaves. A total of 173,354 non-redundant transcripts were generated from the PacBio platform, and three full-length transcripts of β-amyrin synthase, the key synthase of β-amyrin, were identified. A total of 132,610 clean reads obtained from the DNBSEQ platform were utilised to explore key genes related to the triterpenoid saponin biosynthetic pathway in P. grandiflorus, and 96 differentially expressed genes were selected as candidates. The expression levels of these genes were verified by quantitative real-time PCR. Our reliable transcriptome data provide valuable information on the related biosynthesis pathway and may provide insights into the molecular mechanisms of triterpenoid saponin biosynthesis in P. grandiflorus.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Fritillaria hupehensis

One of the most commonly utilized medicinal plants in China is Fritillaria hupehensis (Hsiao et K.C. Hsia). However, due to a lack of genomic resources, little is known about the biosynthesis of relevant compounds, particularly the flavonoid biosynthesis pathway. A PacBio RS II sequencing generated a total of 342,044 reads from the bulb, leaf, root, and stem, of which 316,438 were full-length (FL) non-redundant reads with an average length of 1365 bp and a N50 of 1888 bp. There were also 38,607 long non-coding RNAs and 7914 simple sequence repeats detected. To improve our understanding of processes implicated in regulating secondary metabolite biosynthesis in F. hupehensis tissues, we evaluated potential metabolic pathways. Overall, this study provides a repertoire of FL transcripts in F. hupehensis for the first time, and it will be a valuable resource for marker-assisted breeding and research into bioactive compounds for medicinal and pharmacological applications.

Bayogenin 3-O-cellobioside confers non-cultivar-specific defence against the rice blast fungus Pyricularia oryzae

Rice cultivars from japonica and indica lineage possess differential resistance against blast fungus as a result of genetic divergence. Whether different rice cultivars also show distinct metabolomic changes in response to P. oryzae, and their role in host resistance, are poorly understood. Here, we examine the responses of six different rice cultivars from japonica and indica lineage challenged with P. oryzae. Both susceptible and resistant rice cultivars expressed several metabolites exclusively during P. oryzae infection, including the saponin Bayogenin 3-O-cellobioside. Bayogenin 3-O-cellobioside level in infected rice directly correlated with their resistant attributes. These findings reveal, for the first time to our knowledge that besides oat, other grass plants including rice produces protective saponins. Our study provides insight into the role of pathogen-mediated metabolomics reprogramming in host immunity. The correlation between Bayogenin 3-O-Cellobioside levels and blast resistance suggests that engineering saponin expression in cereal crops represents attractive and sustainable disease management.

A cellulose synthase-derived enzyme catalyses 3-O-glucuronosylation in saponin biosynthesis

Triterpenoid saponins are specialised metabolites distributed widely in the plant kingdom that consist of one or more sugar moieties attached to triterpenoid aglycones. Despite the widely accepted view that glycosylation is catalysed by UDP-dependent glycosyltransferase (UGT), the UGT which catalyses the transfer of the conserved glucuronic acid moiety at the C-3 position of glycyrrhizin and various soyasaponins has not been determined. Here, we report that a cellulose synthase superfamily-derived glycosyltransferase (CSyGT) catalyses 3-O-glucuronosylation of triterpenoid aglycones. Gene co-expression analyses of three legume species (Glycyrrhiza uralensis, Glycine max, and Lotus japonicus) reveal the involvement of CSyGTs in saponin biosynthesis, and we characterise CSyGTs in vivo using Saccharomyces cerevisiae. CSyGT mutants of L. japonicus do not accumulate soyasaponin, but the ectopic expression of endoplasmic reticulum membrane-localised CSyGTs in a L. japonicus mutant background successfully complement soyasaponin biosynthesis. Finally, we produced glycyrrhizin de novo in yeast, paving the way for sustainable production of high-value saponins.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of β-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

Representational Difference Analysis of Transcripts Involved in Jervine Biosynthesis

Veratrum-type steroidal alkaloids (VSA) are the major bioactive ingredients that strongly determine the pharmacological activities of Veratrum nigrum. Biosynthesis of VSA at the molecular and genetic levels is not well understood. Next-generation sequencing of representational difference analysis (RDA) products after elicitation and precursor feeding was applied to identify candidate genes involved in VSA biosynthesis. A total of 12,048 contigs with a median length of 280 bases were received in three RDA libraries obtained after application of methyl jasmonate, squalene and cholesterol. The comparative analysis of annotated sequences was effective in identifying candidate genes. GABAT2 transaminase and hydroxylases active at C-22, C-26, C-11, and C-16 positions in late stages of jervine biosynthesis were selected. Moreover, genes coding pyrroline-5-carboxylate reductase and enzymes from the short-chain dehydrogenases/reductases family (SDR) associated with the reduction reactions of the VSA biosynthesis process were proposed. The data collected contribute to better understanding of jervine biosynthesis and may accelerate implementation of biotechnological methods of VSA biosynthesis.

Lotus japonicus Triterpenoid Profile and Characterization of the CYP716A51 and LjCYP93E1 Genes Involved in Their Biosynthesis In Planta

Lotus japonicus is an important model legume plant in several fields of research, such as secondary (specialized) metabolism and symbiotic nodulation. This plant accumulates triterpenoids; however, less information regarding its composition, content and biosynthesis is available compared with Medicago truncatula and Glycine max. In this study, we analyzed the triterpenoid content and composition of L. japonicus. Lotus japonicus accumulated C-28-oxidized triterpenoids (ursolic, betulinic and oleanolic acids) and soyasapogenols (soyasapogenol B, A and E) in a tissue-dependent manner. We identified an oxidosqualene cyclase (OSC) and two cytochrome P450 enzymes (P450s) involved in triterpenoid biosynthesis using a yeast heterologous expression system. OSC9 was the first enzyme derived from L. japonicus that showed α-amyrin (a precursor of ursolic acid)-producing activity. CYP716A51 showed triterpenoid C-28 oxidation activity. LjCYP93E1 converted β-amyrin into 24-hydroxy-β-amyrin, a metabolic intermediate of soyasapogenols. The involvement of the identified genes in triterpenoid biosynthesis in L. japonicus plants was evaluated by quantitative real-time PCR analysis. Furthermore, gene loss-of-function analysis of CYP716A51 and LjCYP93E1 was conducted. The cyp716a51-mutant L. japonicus hairy roots generated by the genome-editing technique produced no C-28 oxidized triterpenoids. Likewise, the complete abolition of soyasapogenols and soyasaponin I was observed in mutant plants harboring Lotus retrotransposon 1 (LORE1) in LjCYP93E1. These results indicate that the activities of these P450 enzymes are essential for triterpenoid biosynthesis in L. japonicus. This study increases our understanding of triterpenoid biosynthesis in leguminous plants and provides information that will facilitate further studies of the physiological functions of triterpenoids using L. japonicus.

Multi-omics data-driven investigations of metabolic diversity of plant triterpenoids

The vast majority of structurally diverse metabolites play essential roles in mediating the interactions between plant and environment, and constitute a valuable resource for industrial applications. Recent breakthroughs in sequencing technology have greatly accelerated metabolic studies of natural plant products, providing opportunities to investigate the molecular basis underlying the diversity of specialized plant metabolites through large-scale analysis. Here, we focus on the biosynthesis of plant triterpenoids, especially the three diversifying reactions (cyclization, oxidation and glycosylation) that largely contribute to the structural diversity of triterpenoids. Gene mining through large-scale omics data and functional characterization of metabolic genes including enzymes, transcription factors and transporters could provide important insights into the evolution of specialized plant metabolism and pave the way for the production of high-value metabolites or derivatives using synthetic biology approaches.

Transcriptome sequencing and identification of cytochrome P450 monooxygenases involved in the biosynthesis of maslinic acid and corosolic acid in Avicennia marina

Grey mangrove (Avicennia marina) is a traditional medicine used for the treatment of various diseases, including rheumatism and ulcers; however, the compounds responsible for its curative effects remain largely unknown. Triterpenoids are a diverse group of plant-specialized metabolites derived from a common precursor, 2,3-oxidosqualene. Triterpenoids are potentially responsible for the beneficial effects of A. marina; however, the chemical profiles of triterpenoids in A. marina and their biosynthetic genes have not been identified. Cytochrome P450 monooxygenases (P450s) have key roles in the structural diversification of plant triterpenoids by catalyzing site-specific oxidation of triterpene scaffolds. Recent studies have revealed that the CYP716 family represents the most common clade of P450s involved in triterpenoid biosynthesis. In this study, we performed triterpenoid profiling and RNA sequencing of A. marina leaves. Mining of CYP716 family genes and enzymatic activity assays of encoded proteins revealed that CYP716A259 catalyzed oxidation at the C-28 position of the pentacyclic triterpene skeletons of β-amyrin, α-amyrin, and lupeol to produce oleanolic acid, ursolic acid, and betulinic acid, respectively. The other functionally defined P450, CYP716C53, catalyzed the C-2α hydroxylation of oleanolic acid and ursolic acid to produce maslinic acid and corosolic acid, respectively. The possible involvement of CYP716A259 and CYP716C53 in the biosynthesis of these health-benefiting compounds in A. marina leaves, and the possible contribution of the resulting compounds to the reported bioactivities of A. marina leaf extract, are discussed.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

Comparative analysis of CYP716A subfamily enzymes for the heterologous production of C-28 oxidized triterpenoids in transgenic yeast

Several enzymes of the CYP716A subfamily have been reported to be involved in triterpenoid biosynthesis. Members of this subfamily oxidize various positions along the triterpenoid backbone and the majority of them catalyze a three-step oxidation at the C-28 position. Interestingly, C-28 oxidation is a common feature in oleanolic acid, ursolic acid, and betulinic acid, which are widely distributed in plants and exhibit important biological activities. In this work, three additional CYP716A enzymes isolated from olive, sugar beet, and coffee, were characterized as multifunctional C-28 oxidases. Semi-quantitative comparisons of in vivo catalytic activity were made against the previously characterized enzymes CYP716A12, CYP716A15, and CYP716A52v2. When heterologously expressed in yeast, the isolated enzymes differed in both catalytic activity and substrate specificity. This study indicates that the screening of enzymes from different plants could be a useful means of identifying enzymes with enhanced catalytic activity and desired substrate specificity. Furthermore, we show that "naturally-evolved" enzymes can be useful in the heterologous production of pharmacologically and industrially important triterpenoids.

The Basic Helix-Loop-Helix Transcription Factor GubHLH3 Positively Regulates Soyasaponin Biosynthetic Genes in Glycyrrhiza uralensis

Glycyrrhiza uralensis (licorice) is a widely used medicinal plant belonging to the Fabaceae. Its main active component, glycyrrhizin, is an oleanane-type triterpenoid saponin widely used as a medicine and as a natural sweetener. Licorice also produces other triterpenoids, including soyasaponins. Recent studies have revealed various oxidosqualene cyclases and cytochrome P450 monooxygenases (P450s) required for the biosynthesis of triterpenoids in licorice. Of these enzymes, β-amyrin synthase (bAS) and β-amyrin C-24 hydroxylase (CYP93E3) are involved in the biosynthesis of soyasapogenol B (an aglycone of soyasaponins) from 2,3-oxidosqualene. Although these biosynthetic enzyme genes are known to be temporally and spatially expressed in licorice, the regulatory mechanisms underlying their expression remain unknown. Here, we identified a basic helix-loop-helix (bHLH) transcription factor, GubHLH3, that positively regulates the expression of soyasaponin biosynthetic genes. GubHLH3 preferentially activates transcription from promoters of CYP93E3 and CYP72A566, the second P450 gene newly identified and shown to be responsible for C-22β hydroxylation in soyasapogenol B biosynthesis, in transient co-transfection assays of promoter-reporter constructs and transcription factors. Overexpression of GubHLH3 in transgenic hairy roots of G. uralensis enhanced the expression levels of bAS, CYP93E3 and CYP72A566. Moreover, soyasapogenol B and sophoradiol (22β-hydroxy-β-amyrin), an intermediate between β-amyrin and soyasapogenol B, were increased in transgenic hairy root lines overexpressing GubHLH3. We found that soyasaponin biosynthetic genes and GubHLH3 were co-ordinately up-regulated by methyl jasmonate (MeJA). These results suggest that GubHLH3 regulates MeJA-responsive expression of soyasaponin biosynthetic genes in G. uralensis. The regulatory mechanisms of triterpenoid biosynthesis in legumes are compared and discussed.

A multifunctional oxidosqualene cyclase from Tripterygium regelii that produces both α- and β-amyrin

Tripterygium regelii is a rich source of triterpenoids, containing many types of triterpenes with high chemical diversity and interesting pharmacological properties. The cDNA of the multifunctional oxidosqualene cyclase (TrOSC, GenBank accession number: MH161182), consisting of a 2289 bp open reading frame and coding for 762 amino acids, was cloned from the stems and roots of Tripterygium regelii. Phylogenetic analysis using OSC genes from other plants suggested that TrOSC might be a mixed-amyrin synthase. The coding sequence was cloned into the expression vector pYES2 and transformed into the yeast Saccharomyces cerevisiae. The resulting products were analysed by GC-MS. Surprisingly, although it showed 76% sequence identity to lupeol synthase from Ricinus communis, TrOSC was found to be a multifunctional triterpene synthase producing both α- and β-amyrin, the precursors of ursane and oleanane type triterpenes, respectively. qRT-PCR analysis revealed that the transcript of TrOSC accumulated mainly in roots and stems. Taken together, our findings contribute to the knowledge of key genes in the pentacyclic triterpene biosynthesis pathway.

A multifunctional oxidosqualene cyclase was cloned from Tripterygium regelii and identified as a mixed-amyrin synthase, which can produce both α- and β-amyrin.

Researches on Transcriptome Sequencing in the Study of Traditional Chinese Medicine

Due to its incomparable advantages, the application of transcriptome sequencing in the study of traditional Chinese medicine attracts more and more attention of researchers, which greatly promote the development of traditional Chinese medicine. In this paper, the applications of transcriptome sequencing in traditional Chinese medicine were summarized by reviewing recent related papers.

Triterpene Structural Diversification by Plant Cytochrome P450 Enzymes

Cytochrome P450 monooxygenases (P450s) represent the largest enzyme family of the plant metabolism. Plants typically devote about 1% of the protein-coding genes for the P450s to execute primary metabolism and also to perform species-specific specialized functions including metabolism of the triterpenes, isoprene-derived 30-carbon compounds. Triterpenes constitute a large and structurally diverse class of natural products with various industrial and pharmaceutical applications. P450-catalyzed structural modification is crucial for the diversification and functionalization of the triterpene scaffolds. In recent times, a remarkable progress has been made in understanding the function of the P450s in plant triterpene metabolism. So far, ∼80 P450s are assigned biochemical functions related to the plant triterpene metabolism. The members of the subfamilies CYP51G, CYP85A, CYP90B-D, CYP710A, CYP724B, and CYP734A are generally conserved across the plant kingdom to take part in plant primary metabolism related to the biosynthesis of essential sterols and steroid hormones. However, the members of the subfamilies CYP51H, CYP71A,D, CYP72A, CYP81Q, CYP87D, CYP88D,L, CYP93E, CYP705A, CYP708A, and CYP716A,C,E,S,U,Y are required for the metabolism of the specialized triterpenes that might perform species-specific functions including chemical defense toward specialized pathogens. Moreover, a recent advancement in high-throughput sequencing of the transcriptomes and genomes has resulted in identification of a large number of candidate P450s from diverse plant species. Assigning biochemical functions to these P450s will be of interest to extend our knowledge on triterpene metabolism in diverse plant species and also for the sustainable production of valuable phytochemicals.