The ancient CYP716 family is a major contributor to the diversification of eudicot triterpenoid biosynthesis

Advances in production and structural derivatization of the promising molecule ursolic acid

Ursolic acid (UA) is a ursane-type pentacyclic triterpenoid compound, naturally produced in plants via specialized metabolism and exhibits vast range of remarkable physiological activities and pharmacological manifestations. Owing to significant safety and efficacy in different medical conditions, UA may serve as a backbone to produce its derivatives with novel therapeutic functions. This review aims to provide ideas for exploring more diverse structures to improve UA pharmacological activity and increasing its biological yield to meet the industrial requirements by systematically reviewing the current research progress of UA. We first provides an overview of the pharmacological activities, acquisition methods and structural modifications of UA. Among them, we focused on the synthetic modifications of UA to yield valuable derivatives with enhanced therapeutic potential. Furthermore, harnessing the essential advances for green synthesis of UA and its derivatives by advent of metabolic engineering and synthetic biology are of great concern. In this regard, all pivotal advances for enhancing the production of UA have been discussed. In combination with the advantages of UA biosynthesis and transformation strategy, large-scale microbial production of UA is a promising platform for further exploration.

Immunostimulatory 6/6/6/6 Tetracyclic Triterpenoid Saponins with the Methyl-30 Incorporated Cyclization from the Root of Colquhounia elegans

Two novel triterpenoid saponins, colqueleganoids A (1) and B (2), with the first methyl-30 incorporated 6/6/6/6-cyclized carbon skeleton (named colquelegane), were isolated from the root of Colquhounia elegans. Their structures including absolute configuration were determined by spectroscopic methods and X-ray crystallographic analyses. Interestingly, both compounds significantly enhanced TNF-α production and 1 also increased the IL-6 production in RAW264.7 macrophages stimulated with lipopolysaccharide (LPS), suggesting their potential application as immunostimulants in immunotherapy and vaccination.

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.

BAHD acetyltransferase contributes to wound-induced biosynthesis of oleo-gum resin triterpenes in Boswellia

Triterpenes (30-carbon isoprene compounds) represent a large and highly diverse class of natural products that play various physiological functions in plants. The triterpene biosynthetic enzymes, particularly those catalyzing the late-stage regio-selective modifications are not well characterized. The bark of select Boswellia trees, e.g., B. serrata exudes specialized oleo-gum resin in response to wounding, which is enriched with boswellic acids (BAs), a unique class of C3α-epimeric pentacyclic triterpenes with medicinal properties. The bark possesses a network of resin secretory structures comprised of vertical and horizontal resin canals, and amount of BAs in bark increases considerably in response to wounding. To investigate BA biosynthetic enzymes, we conducted tissue-specific transcriptome profiling and identified a wound-responsive BAHD acetyltransferase (BsAT1) of B. serrata catalyzing the late-stage C3α-O-acetylation reactions in the BA biosynthetic pathway. BsAT1 catalyzed C3α-O-acetylation of αBA, βBA, and 11-keto-βBA in vitro and in planta assays to produce all the major C3α-O-acetyl-BAs (3-acetyl-αBA, 3-acetyl-βBA, and 3-acetyl-11-keto-βBA) found in B. serrata bark and oleo-gum resin. BsAT1 showed strict specificity for BA scaffold, whereas it did not acetylate the more common C3β-epimeric pentacyclic triterpenes. The analysis of steady-state kinetics using various BAs revealed distinct substrate affinity and catalytic efficiency. BsAT1 transcript expression coincides with increased levels of C3α-O-acetyl-BAs in bark in response to wounding, suggesting a role of BsAT1 in wound-induced biosynthesis of C3α-O-acetyl-BAs. Overall, the results provide new insights into the biosynthesis of principal chemical constituents of Boswellia oleo-gum resin.

De Novo Biosynthesis of the Oleanane-Type Triterpenoids of Tunicosaponins in Yeast

Tunicosaponins are natural products extracted from Psammosilene tunicoides, which is an important ingredient of Yunnan Baiyao Powder, an ancient and famous Asian herbal medicine. The representative aglycones of tunicosaponins are the oleanane-type triterpenoids of gypsogenin and quillaic acid, which were found to manipulate a broad range of virus-host fusion via wrapping the heptad repeat-2 (HR2) domain prevalent in viral envelopes. However, the unknown biosynthetic pathway and difficulty in chemical synthesis hinder the therapeutic use of tunicosaponins. Here, two novel cytochrome P450-dependent monooxygenases that take part in the biosynthesis of tunicosaponins, CYP716A262 (CYP091) and CYP72A567 (CYP099), were identified from P. tunicoides. In addition, the whole biosynthesis pathway of the tunicosaponin aglycones was reconstituted in yeast by transforming the platform strain BY-bAS with the CYP716A262 and CYP716A567 genes, the resulting strain could produce 146.84 and 314.01 mg/L of gypsogenin and quillaic acid, respectively. This synthetic biology platform for complicated metabolic pathways elucidation and microbial cell factories construction can provide alternative sources of important natural products, helping conserve natural plant resources.

Plant cytochrome P450 plasticity and evolution

The superfamily of cytochrome P450 (CYP) enzymes plays key roles in plant evolution and metabolic diversification. This review provides a status on the CYP landscape within green algae and land plants. The 11 conserved CYP clans known from vascular plants are all present in green algae and several green algae-specific clans are recognized. Clan 71, 72, and 85 remain the largest CYP clans and include many taxa-specific CYP (sub)families reflecting emergence of linage-specific pathways. Molecular features and dynamics of CYP plasticity and evolution are discussed and exemplified by selected biosynthetic pathways. High substrate promiscuity is commonly observed for CYPs from large families, favoring retention of gene duplicates and neofunctionalization, thus seeding acquisition of new functions. Elucidation of biosynthetic pathways producing metabolites with sporadic distribution across plant phylogeny reveals multiple examples of convergent evolution where CYPs have been independently recruited from the same or different CYP families, to adapt to similar environmental challenges or ecological niches. Sometimes only a single or a few mutations are required for functional interconversion. A compilation of functionally characterized plant CYPs is provided online through the Plant P450 Database (erda.dk/public/vgrid/PlantP450/).

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.

Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe

Whole-genome duplication (WGD) plays important roles in plant evolution and function, yet little is known about how WGD underlies metabolic diversification of natural products that bear significant medicinal properties, especially in nonmodel trees. Here, we reveal how WGD laid the foundation for co-option and differentiation of medicinally important ursane triterpene pathway duplicates, generating distinct chemotypes between species and between developmental stages in the apple tribe. After generating chromosome-level assemblies of a widely cultivated loquat variety and Gillenia trifoliata, we define differentially evolved, duplicated gene pathways and date the WGD in the apple tribe at 13.5 to 27.1 Mya, much more recent than previously thought. We then functionally characterize contrasting metabolic pathways responsible for major triterpene biosynthesis in G. trifoliata and loquat, which pre- and postdate the Maleae WGD, respectively. Our work mechanistically details the metabolic diversity that arose post-WGD and provides insights into the genomic basis of medicinal properties of loquat, which has been used in both traditional and modern medicines.

Biosynthetic diversity in triterpene cyclization within the Boswellia genus

This review is not intended to describe the triterpenes isolated from the Boswellia genus, since this information has been covered elsewhere. Instead, the aim is to provide insights into the biosynthesis of triterpenes in Boswellia. This genus, which has 24 species, displays fascinating structural diversity and produces a number of medicinally important triterpenes, particularly boswellic acids. Over 300 volatile components have been reported in the essential oil of Boswellia, and more than 100 diterpenes and triterpenes have been isolated from this genus. Given that no triterpene biosynthetic enzymes have yet been isolated from any members of the Boswellia genus, this review will cover the likely biosynthetic pathways as inferred from structures in nature and the probable types of biosynthetic enzymes based on knowledge of triterpene biosynthesis in other plant species. It highlights the importance of frankincense and the factors and threats affecting its production. It covers triterpene biosynthesis in the genus Boswellia, including dammaranes, tirucallic acids, lupanes, oleananes, ursanes and boswellic acids. Strategies for elucidating triterpene biosynthetic pathways in Boswellia are considered. Furthermore, the possible mechanisms behind wound-induced resin synthesis by the tree and related gene expression profiling are covered. In addition, the influence of the environment and the genotype on the biosynthesis of resin and on variations in the compositions and types of resins will also be reviewed.

Modulation of Arabidopsis root growth by specialized triterpenes

Plant roots are specialized belowground organs that spatiotemporally shape their development in function of varying soil conditions. This root plasticity relies on intricate molecular networks driven by phytohormones, such as auxin and jasmonate (JA). Loss-of-function of the NOVEL INTERACTOR OF JAZ (NINJA), a core component of the JA signaling pathway, leads to enhanced triterpene biosynthesis, in particular of the thalianol gene cluster, in Arabidopsis thaliana roots. We have investigated the biological role of thalianol and its derivatives by focusing on Thalianol Synthase (THAS) and Thalianol Acyltransferase 2 (THAA2), two thalianol cluster genes that are upregulated in the roots of ninja mutant plants. THAS and THAA2 activity was investigated in yeast, and metabolite and phenotype profiling of thas and thaa2 loss-of-function plants was carried out. THAA2 was shown to be responsible for the acetylation of thalianol and its derivatives, both in yeast and in planta. In addition, THAS and THAA2 activity was shown to modulate root development. Our results indicate that the thalianol pathway is not only controlled by phytohormonal cues, but also may modulate phytohormonal action itself, thereby affecting root development and interaction with the environment.

The genome of Tripterygium wilfordii and characterization of the celastrol biosynthesis pathway

Tripterygium wilfordii is a vine from the Celastraceae family that is used in traditional Chinese medicine (TCM). The active ingredient, celastrol, is a friedelane-type pentacyclic triterpenoid with putative roles as an antitumor, immunosuppressive, and anti-obesity agent. Here, we report a reference genome assembly of T. wilfordii with high-quality annotation using a hybrid sequencing strategy. The total genome size obtained is 340.12 Mb, with a contig N50 value of 3.09 Mb. We successfully anchored 91.02% of sequences into 23 pseudochromosomes using high-throughput chromosome conformation capture (Hi–C) technology. The super-scaffold N50 value was 13.03 Mb. We also annotated 31,593 structural genes, with a repeat percentage of 44.31%. These data demonstrate that T. wilfordii diverged from Malpighiales species approximately 102.4 million years ago. By integrating genome, transcriptome and metabolite analyses, as well as in vivo and in vitro enzyme assays of two cytochrome P450 (CYP450) genes, TwCYP712K1 and TwCYP712K2, it is possible to investigate the second biosynthesis step of celastrol and demonstrate that this was derived from a common ancestor. These data provide insights and resources for further investigation of pathways related to celastrol, and valuable information to aid the conservation of resources, as well as understand the evolution of Celastrales.

Cloning, Yeast Expression, and Characterization of a β-Amyrin C-28 Oxidase (CYP716A249) Involved in Triterpenoid Biosynthesis in Polygala tenuifolia

Polygala tenuifolia Willd. is a traditional Chinese herbal medicine that is widely used in treating nervous system disorders. Triterpene saponins in P. tenuifolia (polygala saponins) have excellent biological activity. As a precursor for the synthesis of presenegin, oleanolic acid (OA) plays an important role in the biosynthesis of polygala saponins. However, the mechanism behind the biosynthesis of polygala saponins remains to be elucidated. In this study, we found that CYP716A249 (GenBank: ASB17946) oxidized the C-28 position of β-amyrin to produce OA. Using quantitative real-time PCR, we observed that CYP716A249 had the highest expression in the roots of 2-year-old P. tenuifolia, which provided a basis for the selection of samples for gene cloning. To identify the function of CYP716A249, the strain R-BE-20 was constructed by expressing β-amyrin synthase in yeast. Then, CYP716A249 was co-expressed with β-amyrin synthase to construct the strain R-BPE-20 by using the lithium acetate method. Finally, we detected β-amyrin and OA by ultra-HPLC-Q Exactive hybrid quadrupole-Orbitrap high-resolution accurate mass spectrometry and GC-MS. The results of this study provide insights into the biosynthesis pathway of polygala saponins.

Cancer cells employ an evolutionarily conserved polyploidization program to resist therapy

Unusually large cancer cells with abnormal nuclei have been documented in the cancer literature since 1858. For more than 100 years, they have been generally disregarded as irreversibly senescent or dying cells, too morphologically misshapen and chromatin too disorganized to be functional. Cell enlargement, accompanied by whole genome doubling or more, is observed across organisms, often associated with mitigation strategies against environmental change, severe stress, or the lack of nutrients. Our comparison of the mechanisms for polyploidization in other organisms and non-transformed tissues suggest that cancer cells draw from a conserved program for their survival, utilizing whole genome doubling and pausing proliferation to survive stress. These polyaneuploid cancer cells (PACCs) are the source of therapeutic resistance, responsible for cancer recurrence and, ultimately, cancer lethality.

Dissection of the general two-step di-C-glycosylation pathway for the biosynthesis of (iso)schaftosides in higher plants

Schaftoside and isoschaftoside are bioactive natural products widely distributed in higher plants including cereal crops and medicinal herbs. Their biosynthesis may be related with plant defense. However, little is known on the glycosylation biosynthetic pathway of these flavonoid di-C-glycosides with different sugar residues. Herein, we report that the biosynthesis of (iso)schaftosides is sequentially catalyzed by two C-glycosyltransferases (CGTs), i.e., CGTa for C-glucosylation of the 2-hydroxyflavanone aglycone and CGTb for C-arabinosylation of the mono-C-glucoside. The two enzymes of the same plant exhibit high homology but remarkably different sugar acceptor and donor selectivities. A total of 14 CGTa and CGTb enzymes were cloned and characterized from seven dicot and monocot plants, including Scutellaria baicalensis, Glycyrrhiza uralensis, Oryza sativa ssp. japonica, and Zea mays, and the in vivo functions for three enzymes were verified by RNA interference and overexpression. Through transcriptome analysis, we found homologous genes in 119 other plants, indicating this pathway is general for the biosynthesis of (iso)schaftosides. Furthermore, we resolved the crystal structures of five CGTs and realized the functional switch of SbCGTb to SbCGTa by structural analysis and mutagenesis of key amino acids. The CGT enzymes discovered in this paper allow efficient synthesis of (iso)schaftosides, and the general glycosylation pathway presents a platform to study the chemical defense mechanisms of higher plants.

Differential Expression of Genes Involved in Saikosaponin Biosynthesis Between Bupleurum chinense DC. and Bupleurum scorzonerifolium Willd

Radix Bupleuri (roots of Bupleurum spp.) is an important medicinal herb. Triterpenoid saponins of saikosaponins generally constitute the main class of secondary metabolites of plants in the Bupleurum genus. However, the molecular regulatory mechanism underlying their biosynthesis remains elusive. In this study, we observed significantly different saikosaponin biosynthesis between Bupleurum chinense and Bupleurum scorzonerifolium at the seedling stage. The sequential and expression characterization of 232 genes in the triterpenoid saponin biosynthetic pathway, which includes the mevalonate (MVA) pathway and methylerythritol phosphate (MEP) pathway, between B. chinense and B. scorzonerifolium was also investigated. Sixty of these genes may be involved in saikosaponin biosynthesis. Manipulation of these genes, especially those of the β-AS, P450, and UGT families, may improve saikosaponin production.

Engineering Critical Enzymes and Pathways for Improved Triterpenoid Biosynthesis in Yeast

Triterpenoids represent a diverse group of phytochemicals that are widely distributed in the plant kingdom and have many biological activities. The heterologous production of triterpenoids in Saccharomyces cerevisiae has been successfully implemented by introducing various triterpenoid biosynthetic pathways. By engineering related enzymes as well as through yeast metabolism, the yield of various triterpenoids is significantly improved from the milligram per liter scale to the gram per liter scale. This achievement demonstrates that engineering critical enzymes is considered a potential strategy to overcome the main hurdles of the industrial application of these potent natural products. Here, we review strategies for designing enzymes to improve the yield of triterpenoids in S. cerevisiae in terms of three main aspects: 1, elevating the supply of the precursor 2,3-oxidosqualene; 2, optimizing triterpenoid-involved reactions; and 3, lowering the competition of the native sterol pathway. Then, we provide challenges and prospects for further enhancing triterpenoid production in S. cerevisiae.

Triterpenoid saponins in tea (Camellia sinensis) plants: biosynthetic gene expression, content variations, chemical identification and cytotoxicity

Very little is known about saponins in tea and their biosynthesis in tea plants despite of the importance. Here, we studied tea saponins and their biosynthesis genes. Saponins were promptly recovered in tea infusions. Cytotoxicity of tea saponin extracts on human tongue squamous and hepatocellular carcinoma lines showed respective IC50 values of 29.2 and 17.5 μg/mL, which may be attributable to over 40 saponins identified in green tea. Saponin contents varied in shoot tips of 42 tea plant varieties but did not change drastically during tea processing. Saponin biosynthetic gene expression was consistent with its contents in plant tissues. Thus, plant tips produce significant amounts of saponins, which are stable during tea processing, and ready to be recovered to tea infusions to provide potent health benefits to consumers. This study paves a road towards clarifying the biosynthesis and genetic improvement of saponins in tea plants.

Drivers of metabolic diversification: how dynamic genomic neighbourhoods generate new biosynthetic pathways in the Brassicaceae

Plants produce an array of specialized metabolites with important ecological functions. The mechanisms underpinning the evolution of new biosynthetic pathways are not well-understood. Here, we exploit available genome sequence resources to investigate triterpene biosynthesis across the Brassicaceae. Oxidosqualene cyclases (OSCs) catalyze the first committed step in triterpene biosynthesis. Systematic analysis of 13 sequenced Brassicaceae genomes was performed to identify all OSC genes. The genome neighbourhoods (GNs) around a total of 163 OSC genes were investigated to identify Pfam domains significantly enriched in these regions. All-vs-all comparisons of OSC neighbourhoods and phylogenomic analysis were used to investigate the sequence similarity and evolutionary relationships of the numerous candidate triterpene biosynthetic gene clusters (BGCs) observed. Functional analysis of three representative BGCs was carried out and their triterpene pathway products were elucidated. Our results indicate that plant genomes are remarkably plastic, and that dynamic GNs generate new biosynthetic pathways in different Brassicaceae lineages by shuffling the genes encoding a core palette of triterpene-diversifying enzymes, presumably in response to strong environmental selection pressure. These results illuminate a genomic basis for diversification of plant-specialized metabolism through natural combinatorics of enzyme families, which can be mimicked using synthetic biology to engineer diverse bioactive molecules.

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.

A Seed-Specific Regulator of Triterpene Saponin Biosynthesis in Medicago truncatula

Plants produce a vast array of defense compounds to protect themselves from pathogen attack or herbivore predation. Saponins are a specific class of defense compounds comprising bioactive glycosides with a steroidal or triterpenoid aglycone backbone. The model legume Medicago truncatula synthesizes two types of saponins, hemolytic saponins and nonhemolytic soyasaponins, which accumulate as specific blends in different plant organs. Here, we report the identification of the seed-specific transcription factor TRITERPENE SAPONIN ACTIVATION REGULATOR3 (TSAR3), which controls hemolytic saponin biosynthesis in developing M. truncatula seeds. Analysis of genes that are coexpressed with TSAR3 in transcriptome data sets from developing M. truncatula seeds led to the identification of CYP88A13, a cytochrome P450 that catalyzes the C-16α hydroxylation of medicagenic acid toward zanhic acid, the final oxidation step of the hemolytic saponin biosynthesis branch in M. truncatula In addition, two uridine diphosphate glycosyltransferases, UGT73F18 and UGT73F19, which glucosylate hemolytic sapogenins at the C-3 position, were identified. The genes encoding the identified biosynthetic enzymes are present in clusters of duplicated genes in the M. truncatula genome. This appears to be a common theme among saponin biosynthesis genes, especially glycosyltransferases, and may be the driving force of the metabolic evolution of saponins.

Identification and analysis of CYP450 and UGT supergene family members from the transcriptome of Aralia elata (Miq.) seem reveal candidate genes for triterpenoid saponin biosynthesis

BACKGROUND

Members of the cytochrome P450 (CYP450) and UDP-glycosyltransferase (UGT) gene superfamily have been shown to play essential roles in regulating secondary metabolite biosynthesis. However, the systematic identification of CYP450s and UGTs has not been reported in Aralia elata (Miq.) Seem, a highly valued medicinal plant.

RESULTS

In the present study, we conducted the RNA-sequencing (RNA-seq) analysis of the leaves, stems, and roots of A. elata, yielding 66,713 total unigenes. Following annotation and KEGG pathway analysis, we were able to identify 64 unigenes related to triterpenoid skeleton biosynthesis, 254 CYP450s and 122 UGTs, respectively. A total of 150 CYP450s and 92 UGTs encoding > 300 amino acid proteins were utilized for phylogenetic and tissue-specific expression analyses. This allowed us to cluster 150 CYP450s into 9 clans and 40 families, and then these CYP450 proteins were further grouped into two primary branches: A-type (53%) and non-A-type (47%). A phylogenetic analysis of 92 UGTs and other plant UGTs led to clustering into 16 groups (A-P). We further assessed the expression patterns of these CYP450 and UGT genes across A. elata tissues, with 23 CYP450 and 16 UGT members being selected for qRT-PCR validation, respectively. From these data, we identified CYP716A295 and CYP716A296 as the candidate genes most likely to be associated with oleanolic acid synthesis, while CYP72A763 and CYP72A776 were identified as being the most likely to play roles in hederagenin biosynthesis. We also selected five unigenes as the best candidates for oleanolic acid 3-O-glucosyltransferase. Finally, we assessed the subcellular localization of three CYP450 proteins within Arabidopsis protoplasts, highlighting the fact that they localize to the endoplasmic reticulum.

CONCLUSIONS

This study presents a systematic analysis of the CYP450 and UGT gene family in A. elata and provides a foundation for further functional characterization of these two multigene families.

A systematic comparison of triterpenoid biosynthetic enzymes for the production of oleanolic acid in Saccharomyces cerevisiae

Triterpenoids are high-value plant metabolites with numerous applications in medicine, agriculture, food, and home and personal care products. However, plants produce triterpenoids in low abundance, and their complex structures make their chemical synthesis prohibitively expensive and often impossible. As such, the yeast Saccharomyces cerevisiae has been explored as an alternative means of production. An important triterpenoid is oleanolic acid because it is the precursor to many bioactive triterpenoids of commercial interest, such as QS-21 which is being evaluated as a vaccine adjuvant in clinical trials against HIV and malaria. Oleanolic acid is derived from 2,3-oxidosqualene (natively produced by yeast) via a cyclisation and a multi-step oxidation reaction, catalysed by a β-amyrin synthase and a cytochrome P450 of the CYP716A subfamily, respectively. Although many homologues have been characterised, previous studies have used arbitrarily chosen β-amyrin synthases and CYP716As to produce oleanolic acid and its derivatives in yeast. This study presents the first comprehensive comparison of β-amyrin synthase and CYP716A enzyme activities in yeast. Strains expressing different homologues are compared for production, revealing 6.3- and 4.5-fold differences in β-amyrin and oleanolic acid productivities and varying CYP716A product profiles, which are important to consider when engineering strains for the production of bioactive oleanolic acid derivatives.

Advances in biosynthesis, regulation, and metabolic engineering of plant specialized terpenoids

Plant specialized terpenoids are natural products that have no obvious role in growth and development, but play many important functional roles to improve the plant's overall fitness. Besides, plant specialized terpenoids have immense value to humans due to their applications in fragrance, flavor, cosmetic, and biofuel industries. Understanding the fundamental aspects involved in the biosynthesis and regulation of these high-value molecules in plants not only paves the path to enhance plant traits, but also facilitates homologous or heterologous engineering for overproduction of target molecules of importance. Recent developments in functional genomics and high-throughput analytical techniques have led to unraveling of several novel aspects involved in the biosynthesis and regulation of plant specialized terpenoids. The knowledge thus derived has been successfully utilized to produce target specialized terpenoids of plant origin in homologous or heterologous host systems by metabolic engineering and synthetic biology approaches. Here, we provide an overview and highlights on advances related to the biosynthetic steps, regulation, and metabolic engineering of plant specialized terpenoids.

Production of the bioactive plant-derived triterpenoid morolic acid in engineered Saccharomyces cerevisiae

Morolic acid is a plant-derived triterpenoid that possesses pharmacological properties such as cytotoxicity, as well as anti-HIV, anti-HSV, anti-inflammatory, and antidiabetic effects. The significant therapeutic properties of morolic acid are desirable in the context of pharmacological and drug development research, but the low accessibility of morolic acid from natural resources limits its applications. In the present study, we developed a microbial system for the production of morolic acid. Using a combinatorial biosynthesis approach, a novel production pathway was constructed in Saccharomycescerevisiae by coexpressing BfOSC2 (germanicol synthase) from Bauhinia forficata and CYP716A49 (triterpene C-28 oxidase) from Beta vulgaris. Moreover, we reconstructed the cellular galactose regulatory network by introducing a chimeric transcriptional activator (fusion of Gal4dbd.ER.VP16) to overdrive the genes under the control of the galactose promoter. We also overexpressed truncated HMG1, encoding feedback-inhibition-resistant form of 3-hydroxy-3-methylglutaryl-coenzyme A reductase 1 and sterol-regulating transcription factor upc2-1, to increase the isoprenoid precursors in the mevalonate pathway. Using this yeast system, we achieved morolic acid production up to 20.7 ± 1.8 mg/L in batch culture. To our knowledge, this is the highest morolic acid titer reported from a heterologous host, indicating a promising approach for obtaining rare natural triterpenoids.

Transcriptome analysis and functional characterization of oxidosqualene cyclases of the arjuna triterpene saponin pathway

Arjuna (Terminalia arjuna) tree has been popular in Indian traditional medicine to treat cardiovascular ailments. The tree accumulates bioactive triterpene glycosides (saponins) and aglycones (sapogenins), in a tissue-preferential manner. Oleanane triterpenes/saponins (derived from β-amyrin) with potential cardioprotective function predominantly accumulate in the bark. However, arjuna triterpene saponin pathway enzymes remain to be identified and biochemically characterized. Here, we employed a combined transcriptomics, metabolomics and biochemical approach to functionally define a suite of oxidosqualene cyclases (OSCs) that catalyzed key reactions towards triterpene scaffold diversification. De novo assembly of 131 millions Illumina NextSeq500 sequencing reads obtained from leaf and stem bark samples led to a total of 156,650 reference transcripts. Four distinct OSCs (TaOSC1-4) with 54-71 % sequence identities were identified and functionally characterized. TaOSC1, TaOSC3 and TaOSC4 were biochemically characterized as β-amyrin synthase, cycloartenol synthase and lupeol synthase, respectively. However, TaOSC2 was found to be a multifunctional OSC producing both α-amyrin and β-amyrin, but showed a preference for α-amyrin product. Both TaOSC1 and TaOSC2 produced β-amyrin, the direct precursor for oleanane triterpene/saponin biosynthesis; but, TaOSC1 transcript expressed preferentially in bark, suggesting a major role of TaOSC1 in the biosynthesis of oleanane triterpenes/saponins in bark.

NAA and 6-BA promote accumulation of oleanolic acid by JA regulation in Achyranthes bidentata Bl

Application of plant growth regulators has become one of the most important means of improving yield and quality of medicinal plants. To understand the molecular basis of phytohormone-regulated oleanolic acid metabolism, RNA-seq was used to analyze global gene expression in Achyranthes bidentata treated with 2.0 mg/L 1-naphthaleneacetic acid (NAA) and 1.0 mg/L 6-benzyladenine (6-BA). Compared with untreated controls, the expression levels of 20,896 genes were significantly altered with phytohormone treatment. We found that 13071 (62.5%) unigenes were up-regulated, and a lot of differentially expressed genes involved in hormone or terpenoid biosynthesis, or transcription factors were significantly up-regulated. These results suggest that oleanolic acid biosynthesis induced by NAA and 6-BA occurs due to the expression of key genes involved in jasmonic acid signal transduction. This study is the first to analyze the production and hormonal regulation of medicinal A. bidentata metabolites at the molecular level. The results herein contribute to a better understanding of the regulation of oleanane-type triterpenoid saponins accumulation and define strategies to improve the yield of these useful metabolites.

Comparative proteome and metabolome analyses of latex-exuding and non-exuding Taraxacum koksaghyz roots provide insights into laticifer biology

Taraxacum koksaghyz has been identified as one of the most promising alternative rubber crops. Its high-quality rubber is produced in the latex of laticifers, a specialized cell type that is organized in a network of elongated tubules throughout the entire plant body. In order to gain insights into the physiological role(s) of latex and hence laticifer biology, we examine the effects of barnase-induced latex RNA degradation on the metabolite and protein compositions in the roots. We established high-quality datasets that enabled precise discrimination between cellular and physiological processes in laticifers and non-laticifer cell types of roots at different vegetative stages. We identified numerous latex-specific proteins, including a perilipin-like protein that has not been studied in plants yet. The barnase-expressing plants revealed a phenotype that did not exude latex, which may provide a valuable genetic basis for future studies of plant-environment interactions concerning latex and also help to clarify the evolution and arbitrary distribution of latex throughout the plant kingdom. The overview of temporal changes in composition and protein abundance provided by our data opens the way for a deeper understanding of the molecular interactions, reactions, and network relationships that underlie the different metabolic pathways in the roots of this potential rubber crop.

Identification of α-Amyrin 28-Carboxylase and Glycosyltransferase From Ilex asprella and Production of Ursolic Acid 28-O-β-D-Glucopyranoside in Engineered Yeast

Ilex asprella is a medicinal plant that is used extensively in southern China. The plant contains ursane-type triterpenoids and triterpenoid saponins which are known to be responsible for its pharmacological activities. Previously, a transcriptomic analysis of I. asprella was carried out and the gene IaAS1, which is important in the formation of the core structure α-amyrin, was identified. However, the genes related to the subsequent derivatization of the core structures of the triterpenoid remain largely unknown. Herein, we describe the cloning and functional characterization of an amyrin 28-carboxylase IaAO1 (designated as IaCYP716A210) and a glycosyltransferase IaAU1 (designated as UGT74AG5), based on transcriptomic data. The expression of IaAO1 in an α-amyrin producing yeast strain led to the accumulation of ursolic acid. An enzyme assay using recombinant protein IaAU1 purified from E. coli revealed that IaAU1 can catalyze the conversion of ursolic acid to ursolic acid 28-O-β-D-glucopyranoside. IaAU1 has regiospecificity for catalyzing the 28-O-glucosylation of ursane-/oleanane-type triterpene acids, as it can also catalyze the conversion of oleanolic acid, hederagenin, and ilexgenin A to their corresponding glycosyl compounds. Moreover, co-expression of IaAO1 and IaAU1 in the α-amyrin-producing yeast strain led to the production of ursolic acid 28-O-β-D-glucopyranoside, although in relatively low amounts. Our study reveals that IaAO1 and IaAU1 might play a role in the biosynthesis of pentacyclic triterpenoid saponins in I. asprella and provides insights into the potential application of metabolic engineering to produce ursane-type triterpene glycosides.

Recent advances of metabolic engineering strategies in natural isoprenoid production using cell factories

This review covers the strategies mostly developed in the last three years for microbial production of isoprenoid, classified according to the engineering targets.

Quantitative trait loci influencing pentacyclic triterpene composition in apple fruit peel

The chemical composition of pentacyclic triterpenes was analysed using a ‘Royal Gala’ x ‘Granny Smith’ segregating population in 2013 and 2015, using apple peels extracted from mature fruit at harvest and after 12 weeks of cold storage. In 2013, 20 compound isoforms from nine unique compound classes were measured for both treatments. In 2015, 20 and 17 compound isoforms from eight unique compound classes were measured at harvest and after cold storage, respectively. In total, 68 quantitative trait loci (QTLs) were detected on 13 linkage groups (LG). Thirty two and 36 QTLs were detected for compounds measured at harvest and after cold storage, respectively. The apple chromosomes with the most QTLs were LG3, LG5, LG9 and LG17. The largest effect QTL was for trihydroxy-urs-12-ene-28-oic acid, located on LG5; this was measured in 2015 after storage, and was inherited from the ‘Royal Gala’ parent (24.9% of the phenotypic variation explained).

Molecular Basis of C-30 Product Regioselectivity of Legume Oxidases Involved in High-Value Triterpenoid Biosynthesis

The triterpenes are structurally diverse group of specialized metabolites with important roles in plant defense and human health. Glycyrrhizin, with a carboxyl group at C-30 of its aglycone moiety, is a valuable triterpene glycoside, the production of which is restricted to legume medicinal plants belonging to the Glycyrrhiza species. Cytochrome P450 monooxygenases (P450s) are important for generating triterpene chemodiversity by catalyzing site-specific oxidation of the triterpene scaffold. CYP72A154 was previously identified from the glycyrrhizin-producing plant Glycyrrhiza uralensis as a C-30 oxidase in glycyrrhizin biosynthesis, but its regioselectivity is rather low. In contrast, CYP72A63 from Medicago truncatula showed superior regioselectivity in C-30 oxidation, improving the production of glycyrrhizin aglycone in engineered yeast. The underlying molecular basis of C-30 product regioselectivity is not well understood. Here, we identified two amino acid residues that control C-30 product regioselectivity and contribute to the chemodiversity of triterpenes accumulated in legumes. Amino acid sequence comparison combined with structural analysis of the protein model identified Leu149 and Leu398 as important amino acid residues for C-30 product regioselectivity. These results were further confirmed by mutagenesis of CYP72A154 homologs from glycyrrhizin-producing species, functional phylogenomics analyses, and comparison of corresponding residues of C-30 oxidase homologs in other legumes. These findings could be combined with metabolic engineering to further enhance the production of high-value triterpene compounds.

CYP712K4 Catalyzes the C-29 Oxidation of Friedelin in the Maytenus ilicifolia Quinone Methide Triterpenoid Biosynthesis Pathway

The native Brazilian plant Maytenus ilicifolia accumulates a set of quinone methide triterpenoids with important pharmacological properties, of which maytenin, pristimerin and celastrol accumulate exclusively in the root bark of this medicinal plant. The first committed step in the quinone methide triterpenoid biosynthesis is the cyclization of 2,3-oxidosqualene to friedelin, catalyzed by the oxidosqualene cyclase friedelin synthase (FRS). In this study, we produced heterologous friedelin by the expression of M. ilicifolia FRS in Nicotiana benthamiana leaves and in a Saccharomyces cerevisiae strain engineered using CRISPR/Cas9. Furthermore, friedelin-producing N. benthamiana leaves and S. cerevisiae cells were used for the characterization of CYP712K4, a cytochrome P450 from M. ilicifolia that catalyzes the oxidation of friedelin at the C-29 position, leading to maytenoic acid, an intermediate of the quinone methide triterpenoid biosynthesis pathway. Maytenoic acid produced in N. benthamiana leaves was purified and its structure was confirmed using high-resolution mass spectrometry and nuclear magnetic resonance analysis. The three-step oxidation of friedelin to maytenoic acid by CYP712K4 can be considered as the second step of the quinone methide triterpenoid biosynthesis pathway, and may form the basis for further discovery of the pathway and heterologous production of friedelanes and ultimately quinone methide triterpenoids.

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.

The TriForC database: a comprehensive up-to-date resource of plant triterpene biosynthesis

Triterpenes constitute a large and important class of plant natural products with diverse structures and functions. Their biological roles range from membrane structural components over plant hormones to specialized plant defence compounds. Furthermore, triterpenes have great potential for a variety of commercial applications such as vaccine adjuvants, anti-cancer drugs, food supplements and agronomic agents. Their biosynthesis is carried out through complicated, branched pathways by multiple enzyme types that include oxidosqualene cyclases, cytochrome P450s, and UDP-glycosyltransferases. Given that the number of characterized triterpene biosynthesis enzymes has been growing fast recently, the need for a database specifically focusing on triterpene enzymology became eminent. Here, we present the TriForC database (http://bioinformatics.psb.ugent.be/triforc/), encompassing a comprehensive catalogue of triterpene biosynthesis enzymes. This highly interlinked database serves as a user-friendly access point to versatile data sets of enzyme and compound features, enabling the scanning of a complete catalogue of experimentally validated triterpene enzymes, their substrates and products, as well as the pathways they constitute in various plant species. The database can be accessed by direct browsing or through convenient search tools including keyword, BLAST, plant species and substructure options. This database will facilitate gene mining and creating genetic toolboxes for triterpene synthetic biology.

Hybrid sequencing of the Gynostemma pentaphyllum transcriptome provides new insights into gypenoside biosynthesis

Background

Gypenosides are a group of triterpene saponins from Gynostemma pentaphyllum that are the same as or very similar to ginsenosides from the Panax species. Several enzymes involved in ginsenoside biosynthesis have been characterized, which provide important clues for elucidating the gypenoside biosynthetic pathway. We suppose that gypenosides and ginsenosides may have a similar biosynthetic mechanism and that the corresponding enzymes in the two pathways may have considerable similarity in their sequences. To further understand gypenoside biosynthesis, we sequenced the G. pentaphyllum transcriptome with a hybrid sequencing-based strategy and then determined the candidate genes involved in this pathway using phylogenetic tree construction and gene expression analysis.

Results

Following the PacBio standard analysis pipeline, 66,046 polished consensus sequences were obtained, while Illumina data were assembled into 140,601 unigenes with Trinity software. Then, these output sequences from the two analytical routes were merged. After removing redundant data with CD-HIT software, a total of 140,157 final unigenes were obtained. After functional annotation, five 2,3-oxidosqualene cyclase genes, 145 cytochrome P450 genes and 254 UDP-glycosyltransferase genes were selected for the screening of genes involved in gypenoside biosynthesis. Using phylogenetic analysis, several genes were divided into the same subfamilies or closely related evolutionary branches with characterized enzymes involved in ginsenoside biosynthesis. Using real-time PCR technology, their expression patterns were investigated in different tissues and at different times after methyl jasmonate induction. Since the genes in the same biosynthetic pathway are generally coexpressed, we speculated that GpOSC1, GpCYP89, and GpUGT35 were the leading candidates for gypenoside biosynthesis. In addition, six GpWRKYs and one GpbHLH might play a possible role in regulating gypenoside biosynthesis.

Conclusions

We developed a hybrid sequencing strategy to obtain longer length transcriptomes with increased accuracy, which will greatly contribute to downstream gene screening and characterization, thus improving our ability to elucidate secondary metabolite biosynthetic pathways. With this strategy, we found several candidate genes that may be involved in gypenoside biosynthesis, which laid an important foundation for the elucidation of this biosynthetic pathway, thus greatly contributing to further research in metabolic regulation, synthetic biology and molecular breeding in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6000-y) contains supplementary material, which is available to authorized users.

Integrated metabolomics identifies CYP72A67 and CYP72A68 oxidases in the biosynthesis of Medicago truncatula oleanate sapogenins

INTRODUCTION

Triterpene saponins are important bioactive plant natural products found in many plant families including the Leguminosae.

OBJECTIVES

We characterize two Medicago truncatula cytochrome P450 enzymes, MtCYP72A67 and MtCYP72A68, involved in saponin biosynthesis including both in vitro and in planta evidence.

METHODS

UHPLC-(-)ESI-QToF-MS was used to profile saponin accumulation across a collection of 106 M. truncatula ecotypes. The profiling results identified numerous ecotypes with high and low saponin accumulation in root and aerial tissues. Four ecotypes with significant differential saponin content in the root and/or aerial tissues were selected, and correlated gene expression profiling was performed.

RESULTS

Correlation analyses between gene expression and saponin accumulation revealed high correlations between saponin content with gene expression of β-amyrin synthase, MtCYP716A12, and two cytochromes P450 genes, MtCYP72A67 and MtCYP72A68. In vivo and in vitro biochemical assays using yeast microsomes containing MtCYP72A67 revealed hydroxylase activity for carbon 2 of oleanolic acid and hederagenin. This finding was supported by functional characterization of MtCYP72A67 using RNAi-mediated gene silencing in M. truncatula hairy roots, which revealed a significant reduction of 2β-hydroxylated sapogenins. In vivo and in vitro assays with MtCYP72A68 produced in yeast showed multifunctional oxidase activity for carbon 23 of oleanolic acid and hederagenin. These findings were supported by overexpression of MtCYP72A68 in M. truncatula hairy roots, which revealed significant increases of oleanolic acid, 2β-hydroxyoleanolic acid, hederagenin and total saponin levels.

CONCLUSIONS

The cumulative data support that MtCYP72A68 is a multisubstrate, multifunctional oxidase and MtCYP72A67 is a 2β-hydroxylase, both of which function during the early steps of triterpene-oleanate sapogenin biosynthesis.

Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering

BACKGROUND

Betulinic acid is a pentacyclic lupane-type triterpenoid and a potential antiviral and antitumor drug, but the amount of betulinic acid in plants is low and cannot meet the demand for this compound. Yarrowia lipolytica, as an oleaginous yeast, is a promising microbial cell factory for the production of highly hydrophobic compounds due to the ability of this organism to accumulate large amounts of lipids that can store hydrophobic products and supply sufficient precursors for terpene synthesis. However, engineering for the heterologous production of betulinic acid and related triterpenoids has not developed as systematically as that for the production of other terpenoids, thus the production of betulinic acid in microbes remains unsatisfactory.

RESULTS

In this study, we applied a multimodular strategy to systematically improve the biosynthesis of betulinic acid and related triterpenoids in Y. lipolytica by engineering four functional modules, namely, the heterogenous CYP/CPR, MVA, acetyl-CoA generation, and redox cofactor supply modules. First, by screening 25 combinations of cytochrome P450 monooxygenases (CYPs) and NADPH-cytochrome P450 reductases (CPRs), each of which originated from 5 different sources, we selected two optimal betulinic acid-producing strains. Then, ERG1, ERG9, and HMG1 in the MVA module were overexpressed in the two strains, which dramatically increased betulinic acid production and resulted in a strain (YLJCC56) that exhibited the highest betulinic acid yield of 51.87 ± 2.77 mg/L. Then, we engineered the redox cofactor supply module by introducing NADPH- or NADH-generating enzymes and the acetyl-CoA generation module by directly overexpressing acetyl-CoA synthases or reinforcing the β-oxidation pathway, which further increased the total triterpenoid yield (the sum of the betulin, betulinic acid, betulinic aldehyde yields). Finally, we engineered these modules in combination, and the total triterpenoid yield reached 204.89 ± 11.56 mg/L (composed of 65.44% betulin, 23.71% betulinic acid and 10.85% betulinic aldehyde) in shake flask cultures.

CONCLUSIONS

Here, we systematically engineered Y. lipolytica and achieved, to the best of our knowledge, the highest betulinic acid and total triterpenoid yields reported in microbes. Our study provides a suitable reference for studies on heterologous exploitation of P450 enzymes and manipulation of triterpenoid production in Y. lipolytica.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of β-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

The enzymes OSC1 and CYP716A263 produce a high variety of triterpenoids in the latex of Taraxacum koksaghyz

Only very little is known about the resin composition of natural rubber from the dandelion species Taraxacum koksaghyz, thus its full characterization could provide new insights into how the isoprenoid end-products influence the physical properties of natural rubber, and this resin might be a good source of highly diverse triterpenoids. Here, we present a comprehensive analysis of the triterpenoid composition in an acetone extract and identified 13 triterpenes and triterpenoids also including the so far unknown pentacyclic compounds lup-19(21)-en-3-ol (1) and its ketone lup-19(21)-en-3-one (2). We purified single triterpenes from the acetone extract by developing a two-step HPLC system that is adapted to the structural differences of the described triterpenoids. Furthermore, we isolated six different oxidosqualene cyclases (OSCs) and two P450 enzymes, and we functionally characterized TkOSC1 and CYP716A263 in Nicotiana benthamiana and Saccharomyces cerevisiae in detail. TkOSC1 is a multifunctional OSC that was capable of synthesizing at least four of the latex-predominant pentacyclic triterpenes (taraxasterol, α-, β-amyrin and lup-19(21)-en-3-ol) while CYP716A263 oxidized pentacyclic triterpenes at the C-3 position. The identified enzymes responsible for biosynthesis and modification of pentacyclic triterpenes in T. koksaghyz latex may represent excellent tools for bioengineering approaches to produce pentacyclic triterpenes heterologously.

Oxidosqualene cyclase and CYP716 enzymes contribute to triterpene structural diversity in the medicinal tree banaba

Pentacyclic triterpenes (PCTs) represent a major class of bioactive metabolites in banaba (Lagerstroemia speciosa) leaves; however, biosynthetic enzymes and their involvement in the temporal accumulation of PCTs remain to be studied. We use an integrated approach involving transcriptomics, metabolomics and gene function analysis to identify oxidosqualene cyclases (OSCs) and cytochrome P450 monooxygenases (P450s) that catalyzed sequential cyclization and oxidative reactions towards PCT scaffold diversification. Four monofunctional OSCs (LsOSC1,3-5) converted the triterpene precursor 2,3-oxidosqualene to either lupeol, β-amyrin or cycloartenol, and a multifunctional LsOSC2 formed α-amyrin as a major product along with β-amyrin. Two CYP716 family P450s (CYP716A265, CYP716A266) catalyzed C-28 oxidation of α-amyrin, β-amyrin and lupeol to form ursolic acid, oleanolic acid and betulinic acid, respectively. However, CYP716C55 catalyzed C-2α hydroxylation of ursolic acid and oleanolic acid to produce corosolic acid and maslinic acid, respectively. Besides, combined transcript and metabolite analysis suggested major roles for the LsOSC2, CYP716A265 and CYP716C55 in determining leaf ursane and oleanane profiles. Combinatorial expression of OSCs and CYP716s in Saccharomyces cerevisiae and Nicotiana benthamiana led to PCT pathway reconstruction, signifying the utility of banaba enzymes for bioactive PCT production in alternate plant/microbial hosts that are more easily tractable than the tree species.

Multi-tissue integration of transcriptomic and specialized metabolite profiling provides tools for assessing the common bean (Phaseolus vulgaris) metabolome

Common bean (Phaseolus vulgaris L.) is an important legume species with a rich natural diversity of landraces that originated from the wild forms following multiple independent domestication events. After the publication of its genome, several resources for this relevant crop have been made available. A comprehensive characterization of specialized metabolism in P. vulgaris, however, is still lacking. In this study, we used a metabolomics approach based on liquid chromatography-mass spectrometry to dissect the chemical composition at a tissue-specific level in several accessions of common bean belonging to different gene pools. Using a combination of literature search, mass spectral interpretation, 13 C-labeling, and correlation analyses, we were able to assign chemical classes and/or putative structures for approximately 39% of all measured metabolites. Additionally, we integrated this information with transcriptomics data and phylogenetic inference from multiple legume species to reconstruct the possible metabolic pathways and identify sets of candidate genes involved in the biosynthesis of specialized metabolites. A particular focus was given to flavonoids, triterpenoid saponins and hydroxycinnamates, as they represent metabolites involved in important ecological interactions and they are also associated with several health-promoting benefits when integrated into the human diet. The data are presented here in the form of an accessible resource that we hope will set grounds for further studies on specialized metabolism in legumes.

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.

Biosynthesis of Plant Triterpenoid Saponins in Microbial Cell Factories

Triterpenoid saponins are triterpenoid glycoside compounds which have been widely used in pharmaceutical, agricultural, and food industries. Traditionally, they are extracted from plants, which is time-consuming and environmentally unfriendly. Recently, de novo synthesis of triterpenoid saponins in microbial cell factories was realized, which provides a promising and green approach to alter the traditional supply way. However, the complex biosynthetic pathway and the poor suitability between the endogenous and heterogeneous pathways tremendously limit the yield of triterpenoid saponins. We introduce the biosynthetic pathways of triterpenoid saponins first, and we then summarize the microbial cell factories developed to produce these compounds. Further, we discuss the strategies applied to enhance the production. This paper systematically illustrates the biosynthesis of plant triterpenoid saponins in microbial cell factories.

Transcriptomics-based identification and characterization of 11 CYP450 genes of Panax ginseng responsive to MeJA

Cytochromes P450 (CYP450s), a superfamily of mono-oxygenases, are essential to generate highly functionalized secondary metabolites in plants and contribute to the diversification of specialized triterpenoid biosynthesis in eudicots. However, screening and identifying the exact CYP450 genes in ginsenoside biosynthesis is extremely challenging due to existence of large quantities of members in CYP450 superfamily. Therefore, to screen the CYP450 genes involved in ginsenoside biosynthesis, transcriptome dataset of Panax ginseng was created in our previous work using the technique of the next-generation sequencing. On the basis of bioinformatics analysis, 16 putative CYP450 genes with significant differential expression were screened from the dataset and submitted to GenBank, in which 11 of them have been cloned. Methyl jasmonate (MeJA) was used as an elicitor to analyze the expression profiles of candidate CYP450 genes by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The results of qRT-PCR analysis revealed that the expression of some CYP450 genes were strongly induced by MeJA and showed different transcription levels at different treatment time points. Homology analysis indicated that each putative CYP450 protein of P. ginseng has a conserved domain consisting of E-E-R-F-P-R-G. The CYP450 genes were screened and cloned here to enrich the resources of CYP450 genes, and the results of bioinformatics analysis provided a foundation to further identify the function of CYP450s involved in ginsenoside biosynthesis. Furthermore, this study facilitated the construction of microbial cell factories for increasing the production of ginsenosides by means of metabolic engineering.