Medicago truncatula CYP716A12 Is a Multifunctional Oxidase Involved in the Biosynthesis of Hemolytic Saponins

Roles of three cytochrome P450 monooxygenases in triterpene biosynthesis and their potential impact on growth and development

Pentacyclic triterpenoids, recognized for their natural bioactivity, display complex spatiotemporal accumulation patterns within the ecological model plant Nicotiana attenuata. Despite their ecological importance, the underlying biosynthetic enzymes and functional attributes of triterpenoid synthesis in N. attenuata remain unexplored. Here, we show that 3 cytochrome P450 monooxygenases (NaCYP716A419, NaCYP716C87, and NaCYP716E107) from N. attenuata oxidize the pentacyclic triterpene skeleton, as evidenced by heterologous expression in Nicotiana benthamiana. NaCYP716A419 catalyzed a consecutive 3-step oxidation reaction at the C28 position of β-amyrin/lupeol/lupanediol, yielding the corresponding alcohol, aldehyde, and carboxylic acid. NaCYP716C87 hydroxylated the C2α position of β-amyrin/lupeol/lupanediol/erythrodiol/oleanolic acid/betulinic acid, while NaCYP716E107 hydroxylated the C6β position of β-amyrin/oleanolic acid. The genes encoding these 3 CYP716 enzymes are highly expressed in flowers and respond to induction by ABA, MeJA, SA, GA3, and abiotic stress treatments. Using VIGS technology, we revealed that silencing of NaCYP716A419 affects the growth and reproduction of N. attenuata, suggesting the ecological significance of these specialized metabolite biosynthetic steps.

Transcriptome sequencing and identification of full-length genes involved in the biosynthesis of anticancer compounds Oleanolic acid and Ursolic acid in Achyranthes aspera L

Achyranthes aspera is renowned for its rich medicinal properties since the Ayurvedic era. This plant is known for the presence of experimentally validated anticancer compounds like oleanolic acid (OA) and ursolic acid (UA). Our study involved sequencing the RNA from the root tissue of A. aspera to elucidate the genes responsible for synthesizing these two critical secondary metabolites. Through RNA-Seq analysis, we assembled approximately 167,698 transcripts, averaging 847 base pairs in length, with an N50 value of 1509 bp. From this data, we mapped 604 sequences involved in the metabolism of terpenoids and polyketide pathways. Among them, 241 transcripts were mapped to the triterpenoid biosynthesis pathway, which included 127 transcripts involved in OA and UA biosynthesis. From these transcripts, we identified 22 full-length genes coding for all the 21 enzymes required for OA and UA biosynthesis. Identifying these full-length genes will lead to a better understanding of the pathway and adopting genetic engineering approaches.

Histochemical and ultrastructural localization of triterpene saponins in Medicago truncatula

In the Medicago genus, saponins are complex mixtures of triterpene pentacyclic glycosides extensively studied for their different and economically relevant biological and pharmaceutical properties. This research is aimed at determining for the first time the tissue and cellular localization of triterpene saponins in vegetative organs of Medicago truncatula, a model plant species for legumes, by histochemistry and transmission electron microscopy. The results showed that saponins are present mainly in the palisade mesophyll layer of leaves, whereas in stems they are mostly located in the primary phloem and the subepidermal cells of cortical parenchyma. In root tissue, saponins occur in the secondary phloem region. Transmission electron microscopy revealed prominent saponin accumulation within the leaf and stem chloroplasts, while in the roots the saponins are found in the vesicular structures. Our results demonstrate the feasibility of using histochemistry and transmission electron microscopy to localize M. truncatula saponins at tissue and cellular levels and provide important information for further studies on biosynthesis and regulation of valuable bioactive saponins on agronomic relevant Medicago spp., such as alfalfa (Medicago sativa L.). RESEARCH HIGHLIGHTS: The Medicago genus represents a valuable rich source of saponins, one of the most interesting groups of secondary plant metabolites, which possess relevant biological and pharmacological properties. Plant tissue and cellular localization of saponins is of great importance to better understand their biological functions, biosynthetic pathway, and regulatory mechanisms. We elucidate the localization of saponins in Medicago truncatula with histochemical and transmission electron microscopy studies.

Unraveling the role of cytochrome P450 enzymes in oleanane triterpenoid biosynthesis in arjuna tree

Triterpenoids (C30-isoprenoids) represent a major group of natural products with various physiological functions in plants. Triterpenoids and their derivatives have medicinal uses owing to diverse bioactivities. Arjuna (Terminalia arjuna) tree bark accumulates highly oxygenated β-amyrin-derived oleanane triterpenoids (e.g., arjunic acid, arjungenin, and arjunolic acid) with cardioprotective roles. However, biosynthetic routes and enzymes remain poorly understood. We mined the arjuna transcriptome and conducted cytochrome P450 monooxygenase (P450) assays using Saccharomyces cerevisiae and Nicotiana benthamiana to identify six P450s and two P450 reductases for oxidative modifications of oleanane triterpenoids. P450 assays using oleananes revealed a greater substrate promiscuity of C-2α and C-23 hydroxylases/oxidases than C-28 oxidases. CYP716A233 and CYP716A432 catalyzed β-amyrin/erythrodiol C-28 oxidation to produce oleanolic acid. C-2α hydroxylases (CYP716C88 and CYP716C89) converted oleanolic acid and hederagenin to maslinic acid and arjunolic acid. CYP716C89 also hydroxylated erythrodiol and oleanolic aldehyde. However, CYP714E107a and CYP714E107b catalyzed oleanolic acid/maslinic acid/arjunic acid, C-23 hydroxylation to form hederagenin, arjunolic acid and arjungenin, and hederagenin C-23 oxidation to produce gypsogenic acid, but at a lower rate than oleanolic acid C-23 hydroxylation. Overall, P450 substrate selectivity suggested that C-28 oxidation is the first P450-catalyzed oxidative modification in the arjuna triterpenoid pathway. However, the pathway might branch thereafter through C-2α/C-23 hydroxylation of oleanolic acid. Taken together, these results provided new insights into substrate range of P450s and unraveled biosynthetic routes of triterpenoids in arjuna. Moreover, complete elucidation and reconstruction of arjunolic acid pathway in S. cerevisiae and N. benthamiana suggested the utility of arjuna P450s in heterologous production of cardioprotective compounds.

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

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

Plant triterpenoid saponins function as susceptibility factors to promote the pathogenicity of Botrytis cinerea

The gray mold fungus Botrytis cinerea is a necrotrophic pathogen that causes diseases in hundreds of plant species, including high-value crops. Its polyxenous nature and pathogenic success are due to its ability to perceive host signals in its favor. In this study, we found that laticifer cells of Euphorbia lathyris are a source of susceptibility factors required by B. cinerea to cause disease. Consequently, poor-in-latex (pil) mutants, which lack laticifer cells, show full resistance to this pathogen, whereas lot-of-latex mutants, which produce more laticifer cells, are hypersusceptible. These S factors are triterpenoid saponins, which are widely distributed natural products of vast structural diversity. The downregulation of laticifer-specific oxydosqualene cyclase genes, which encode the first committed step enzymes for triterpene and, therefore, saponin biosynthesis, conferred disease resistance to B. cinerea. Likewise, the Medicago truncatula lha-1 mutant, compromised in triterpenoid saponin biosynthesis, showed enhanced resistance. Interestingly, the application of different purified triterpenoid saponins pharmacologically complemented the disease-resistant phenotype of pil and hla-1 mutants and enhanced disease susceptibility in different plant species. We found that triterpenoid saponins function as plant cues that signal transcriptional reprogramming in B. cinerea, leading to a change in its growth habit and infection strategy, culminating in the abundant formation of infection cushions, the multicellular appressoria apparatus dedicated to plant penetration and biomass destruction in B. cinerea. Taken together, these results provide an explanation for how plant triterpenoid saponins function as disease susceptibility factors to promote B. cinerea pathogenicity.

Cardiac glycosides protect wormseed wallflower (Erysimum cheiranthoides) against some, but not all, glucosinolate-adapted herbivores

The chemical arms race between plants and insects is foundational to the generation and maintenance of biological diversity. We asked how the evolution of a novel defensive compound in an already well-defended plant lineage impacts interactions with diverse herbivores. Erysimum cheiranthoides (Brassicaceae), which produces both ancestral glucosinolates and novel cardiac glycosides, served as a model. We analyzed gene expression to identify cardiac glycoside biosynthetic enzymes in E. cheiranthoides and characterized these enzymes via heterologous expression and CRISPR/Cas9 knockout. Using E. cheiranthoides cardiac glycoside-deficient lines, we conducted insect experiments in both the laboratory and field. EcCYP87A126 initiates cardiac glycoside biosynthesis via sterol side-chain cleavage, and EcCYP716A418 has a role in cardiac glycoside hydroxylation. In EcCYP87A126 knockout lines, cardiac glycoside production was eliminated. Laboratory experiments with these lines revealed that cardiac glycosides were highly effective defenses against two species of glucosinolate-tolerant specialist herbivores, but did not protect against all crucifer-feeding specialist herbivores in the field. Cardiac glycosides had lesser to no effect on two broad generalist herbivores. These results begin elucidation of the E. cheiranthoides cardiac glycoside biosynthetic pathway and demonstrate in vivo that cardiac glycoside production allows Erysimum to escape from some, but not all, specialist herbivores.

RNA-seq analysis reveals narrow differential gene expression in MEP and MVA pathways responsible for phytochemical divergence in extreme genotypes of Thymus daenensis Celak

BACKGROUND

Here, we investigated the underlying transcriptional-level evidence behind phytochemical differences between two metabolically extreme genotypes of Thymus daenensis. The genotypes 'Zagheh-11' (thymol/carvacrol type, poor in essential oil [EO] [2.9%] but rich in triterpenic acids) and 'Malayer-21' (thymol type and rich in EO [3.8%]) were selected from an ongoing breeding program and then clonally propagated for further experimental use.

MATERIALS AND METHODS

GC-MS, GC-FID, and HPLC-PDA were utilized to monitor the fluctuation of secondary metabolites at four phenological stages (vegetative, bud burst, early, and full-flowering stages). The highest phytochemical divergence was observed at early flowering stage. Both genotypes were subjected to mRNA sequencing (approximately 100 million paired reads) at the aforementioned stage. The expression patterns of four key genes involved in the biosynthesis of terpenoids were also validated using qRT-PCR.

RESULTS

Carvacrol content in 'Zagheh-11' (26.13%) was approximately 23 times higher than 'Malayer-21' (1.12%). Reciprocally, about 10% higher thymol was found in 'Malayer-21' (62.15%). Moreover, the concentrations of three major triterpenic acids in 'Zagheh-11' were approximately as twice as those found in 'Malayer-21'. Transcriptome analysis revealed a total of 1840 unigenes that were differentially expressed, including terpene synthases, cytochrome P450, and terpenoid backbone genes. Several differentially expressed transcription factors (such as MYB, bZIP, HB-HD-ZIP, and WRKY families) were also identified. These results suggest that an active cytosolic mevalonate (MVA) pathway may be linked to higher levels of sesquiterpenes, triterpenic acids, and carvacrol in 'Zagheh-11'. The chloroplastic pathway of methyl erythritol phosphate (MEP) may have also contributed to a higher accumulation of thymol in Malayer-21. Indeed, 'Zagheh-11' showed higher expression of certain genes (HMGR, CYP71D180, β-amyrin 28-monooxygenase, and sesquiterpene synthases) in the MVA pathway, while some genes in the MEP pathway (including DXR, ispG, and γ-terpinene synthase) were distinctly expressed in Malayer-21. Future efforts in metabolic engineering of MVA/MEP pathways may benefit from these findings to produce increased levels of desired secondary metabolites at commercial scale.

Cardiac glycosides protect wormseed wallflower (Erysimum cheiranthoides) against some, but not all, glucosinolate-adapted herbivores

The chemical arms race between plants and insects is foundational to the generation and maintenance of biological diversity. We asked how the evolution of a novel defensive compound in an already well-defended plant lineage impacts interactions with diverse herbivores. Erysimum cheiranthoides (Brassicaceae), which produces both ancestral glucosinolates and novel cardiac glycosides, served as a model.We analyzed gene expression to identify cardiac glycoside biosynthetic enzymes in E. cheiranthoides and characterized these enzymes via heterologous expression and CRISPR/Cas9 knockout. Using E. cheiranthoides cardiac glycoside-deficient lines, we conducted insect experiments in both the laboratory and field.EcCYP87A126 initiates cardiac glycoside biosynthesis via sterol side chain cleavage, and EcCYP716A418 has a role in cardiac glycoside hydroxylation. In EcCYP87A126 knockout lines, cardiac glycoside production was eliminated. Laboratory experiments with these lines revealed that cardiac glycosides were highly effective defenses against two species of glucosinolate-tolerant specialist herbivores but did not protect against all crucifer-feeding specialist herbivores in the field. Cardiac glycosides had lesser to no effect on two broad generalist herbivores.These results begin elucidation of the E. cheiranthoides cardiac glycoside biosynthetic pathway and demonstrate in vivo that cardiac glycoside production allows Erysimum to escape from some, but not all, specialist herbivores.

Combinatorial optimization and spatial remodeling of CYPs to control product profile

Activating inert substrates is a challenge in nature and synthetic chemistry, but essential for creating functionally active molecules. In this work, we used a combinatorial optimization approach to assemble cytochrome P450 monooxygenases (CYPs) and reductases (CPRs) to achieve a target product profile. By creating 110 CYP-CPR pairs and iteratively screening different pairing libraries, we demonstrated a framework for establishing a CYP network that catalyzes six oxidation reactions at three different positions of a chemical scaffold. Target product titer was improved by remodeling endoplasmic reticulum (ER) size and spatially controlling the CYPs' configuration on the ER. Out of 47 potential products that could be synthesized, 86% of the products synthesized by the optimized network was our target compound quillaic acid (QA), the aglycone backbone of many pharmaceutically important saponins, and fermentation achieved QA titer 2.23 g/L.

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.

Discovery of unique CYP716C oxidase involved in pentacyclic triterpene biosynthesis from Camptotheca acuminata

Dozens of triterpenes have been isolated from Camptotheca acuminata, however, triterpene metabolism in this plant remains poorly understood. The common C28 carboxy located in the oleanane-type and ursane-type triterpenes indicates the existence of a functionally active triterpene, C28 oxidase, in this plant. Thorough mining and screening of the CYP716 genes were initiated using the multi-omics database for C. acuminata. Two CYP716A (CYP716A394 and CYP716A395) and three CYP716C (CYP716C80-CYP716C82) were identified based on conserved domain analyses and hierarchical cluster analyses. CYP716 microsomal proteins were prepared and their enzymatic activities were evaluated in vitro. The CYP716 classified into the CYP716C subfamily displays β-amyrin oxidation activity, and CYP716A displays α-amyrin and lupeol oxidation activity, based on gas chromatography-mass spectrometry analyses. The oxidation products were determined based on their mass and nuclear magnetic resonance spectrums. The optimum reaction conditions and kinetic parameters for CYP716C were determined, and functions were verified in Nicotiana benthaminana. Relative quantitative analyses revealed that these CYP716C genes were enriched in the leaves of C. acuminata plantlets after 60 d. These results indicate that CYP716C plays a dominant role in oleanane-type triterpene metabolism in the leaves of C. acuminata via a substrate-specific manner, and CYP716A is responsible for ursane- and lupane-type triterpene metabolism in fruit. This study provides valuable insights into the unique CYP716C-mediated oxidation step of pentacyclic triterpene biosynthesis in C. acuminata.

Natural products of pentacyclic triterpenoids: from discovery to heterologous biosynthesis

Covering: up to 2022Pentacyclic triterpenoids are important natural bioactive substances that are widely present in plants and fungi. They have significant medicinal efficacy, play an important role in reducing blood glucose and protecting the liver, and have anti-inflammatory, anti-oxidation, anti-fatigue, anti-viral, and anti-cancer activities. Pentacyclic triterpenoids are derived from the isoprenoid biosynthetic pathway, which generates common precursors of triterpenes and steroids, followed by cyclization with oxidosqualene cyclases (OSCs) and decoration via cytochrome P450 monooxygenases (CYP450s) and glycosyltransferases (GTs). Many biosynthetic pathways of triterpenoid saponins have been elucidated by studying their metabolic regulation network through the use of multiomics and identifying their functional genes. Unfortunately, natural resources of pentacyclic triterpenoids are limited due to their low content in plant tissues and the long growth cycle of plants. Based on the understanding of their biosynthetic pathway and transcriptional regulation, plant bioreactors and microbial cell factories are emerging as alternative means for the synthesis of desired triterpenoid saponins. The rapid development of synthetic biology, metabolic engineering, and fermentation technology has broadened channels for the accumulation of pentacyclic triterpenoid saponins. In this review, we summarize the classification, distribution, structural characteristics, and bioactivity of pentacyclic triterpenoids. We further discuss the biosynthetic pathways of pentacyclic triterpenoids and involved transcriptional regulation. Moreover, the recent progress and characteristics of heterologous biosynthesis in plants and microbial cell factories are discussed comparatively. Finally, we propose potential strategies to improve the accumulation of triterpenoid saponins, thereby providing a guide for their future biomanufacturing.

Class I and II NADPH-cytochrome P450 reductases exhibit different roles in triterpenoid biosynthesis in Lotus japonicus

Cytochrome P450 monooxygenases (CYPs) are enzymes that play critical roles in the structural diversification of triterpenoids. To perform site-specific oxidations of the triterpene scaffold, CYPs require electrons transferred by NADPH-cytochrome P450 reductase (CPR), which is classified into two main classes, class I and class II, based on their structural difference. Lotus japonicus is a triterpenoids-producing model legume with one CPR class I gene (LjCPR1) and a minimum of two CPR class II genes (LjCPR2-1 and LjCPR2-2). CPR classes I and II from different plants have been reported to be involved in different metabolic pathways. By performing gene expression analyses of L. japonicus hairy root culture treated with methyl jasmonate (MeJA), this study revealed that LjCPR1, CYP716A51, and LUS were down-regulated which resulted in no change in betulinic acid and lupeol content. In contrast, LjCPR2s, bAS, CYP93E1, and CYP72A61 were significantly upregulated by MeJA treatment, followed by a significant increase of the precursors for soyasaponins, i.e. β-amyrin, 24-OH β-amyrin, and sophoradiol content. Triterpenoids profile analysis of LORE1 insertion and hairy root mutants showed that the loss of the Ljcpr2-1 gene significantly reduced soyasaponins precursors but not in Ljcpr1 mutants. However, Ljcpr1 and Ljcpr2-1 mutants showed a significant reduction in lupeol and oleanolic, ursolic, and betulinic acid contents. Furthermore, LjCPR1, but not LjCPR2, was crucial for seed development, supporting the previous notion that CPR class I might support plant basal metabolism. This study suggests that CPR classes I and II play different roles in L. japonicus triterpenoid biosynthesis.

Comparative transcriptome analysis to identify putative genes involved in carvacrol biosynthesis pathway in two species of Satureja, endemic medicinal herbs of Iran

Satureja is rich in phenolic monoterpenoids, mainly carvacrol, that is of interest due to diverse biological activities including antifungal and antibacterial. However, limited information is available regarding the molecular mechanisms underlying carvacrol biosynthesis and its regulation for this wonderful medicinal herb. To identify the putative genes involved in carvacrol and other monoterpene biosynthesis pathway, we generated a reference transcriptome in two endemic Satureja species of Iran, containing different yields (Satureja khuzistanica and Satureja rechingeri). Cross-species differential expression analysis was conducted between two species of Satureja. 210 and 186 transcripts related to terpenoid backbone biosynthesis were identified for S. khuzistanica and S. rechingeri, respectively. 29 differentially expressed genes (DEGs) involved in terpenoid biosynthesis were identified, and these DEGs were significantly enriched in monoterpenoid biosynthesis, diterpenoid biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, carotenoid biosynthesis and ubiquinone and other terpenoid-quinone biosynthesis pathways. Expression patterns of S. khuzistanica and S. rechingeri transcripts involved in the terpenoid biosynthetic pathway were evaluated. In addition, we identified 19 differentially expressed transcription factors (such as MYC4, bHLH, and ARF18) that may control terpenoid biosynthesis. We confirmed the altered expression levels of DEGs that encode carvacrol biosynthetic enzymes using quantitative real-time PCR (qRT-PCR). This study is the first report on de novo assembly and transcriptome data analysis in Satureja which could be useful for an understanding of the main constituents of Satureja essential oil and future research in this genus.

The saponin bomb: a nucleolar-localized β-glucosidase hydrolyzes triterpene saponins in Medicago truncatula

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a β-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned β-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Elucidation of the pathway for biosynthesis of saponin adjuvants from the soapbark tree

The Chilean soapbark tree (Quillaja saponaria) produces soap-like molecules called QS saponins that are important vaccine adjuvants. These highly valuable compounds are sourced by extraction from the bark, and their biosynthetic pathway is unknown. Here, we sequenced the Q. saponaria genome. Through genome mining and combinatorial expression in tobacco, we identified 16 pathway enzymes that together enable the production of advanced QS pathway intermediates that represent a bridgehead for adjuvant bioengineering. We further identified the enzymes needed to make QS-7, a saponin with excellent therapeutic properties and low toxicity that is present in low abundance in Q. saponaria bark extract. Our results enable the production of Q. saponaria vaccine adjuvants in tobacco and open the way for new routes to access and engineer natural and new-to-nature immunostimulants.

Betulin, Synthesized by PpCYP716A1, Is a Key Endogenous Defensive Metabolite of Peach against Aphids

Wild pest-resistant germplasms employ secondary metabolites to withstand insect attacks. A close wild relative of the cultivated peach, Prunus davidiana, displays strong resistance to green peach aphids by utilizing metabolites to cope with aphid infestation; however, the underlying mechanism of aphid resistance remains mostly unknown. Here, metabolomic analysis was performed to explore the changes in metabolite levels in P. davidiana after aphid infestation. The data revealed that betulin is a key defensive metabolite in peaches that protects against aphids and possesses potent aphidicidal activity. Further toxicity tests demonstrated that betulin was toxic to pests but not to beneficial insects. Additionally, transcriptomic and phylogenetic analyses revealed that the cytochrome P450 gene PpCYP716A1 was responsible for betulin synthesis─this finding was confirmed by the heterologous expression of this gene. This study revealed a strategy whereby plants harness defense metabolites to develop resistance to pests. These findings may facilitate controlling such pests.

Interference between ER stress-related bZIP-type and jasmonate-inducible bHLH-type transcription factors in the regulation of triterpene saponin biosynthesis in Medicago truncatula

Triterpene saponins (TS) are a structurally diverse group of metabolites that are widely distributed in plants. They primarily serve as defense compounds and their production is often triggered by biotic stresses through signaling cascades that are modulated by phytohormones such as the jasmonates (JA). Two JA-modulated basic helix-loop-helix (bHLH) transcription factors (TFs), triterpene saponin biosynthesis activating regulator 1 (TSAR1) and TSAR2, have previously been identified as direct activators of TS biosynthesis in the model legume Medicago truncatula. Here, we report on the involvement of the core endoplasmic reticulum (ER) stress-related basic leucine zipper (bZIP) TFs bZIP17 and bZIP60 in the regulation of TS biosynthesis. Expression and processing of M. truncatula bZIP17 and bZIP60 proteins were altered in roots with perturbed TS biosynthesis or treated with JA. Accordingly, such roots displayed an altered ER network structure. M. truncatula bZIP17 and bZIP60 proteins were shown to localize in the nucleus and appeared to be capable of interfering with the TSAR-mediated transactivation of TS biosynthesis genes. Furthermore, interference between ER stress-related bZIP and JA-modulated bHLH TFs in the regulation of JA-dependent terpene biosynthetic pathways may be widespread in the plant kingdom, as we demonstrate that it also occurs in the regulation of monoterpene indole alkaloid biosynthesis in the medicinal plant Catharanthus roseus.

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.

Engineering Critical Amino Acid Residues of Lanosterol Synthase to Improve the Production of Triterpenoids in Saccharomyces cerevisiae

Triterpenoids are a subgroup of terpenoids and have wide applications in the food, cosmetics, and pharmaceutical industries. The heterologous production of various triterpenoids in Saccharomyces cerevisiae, as well as other microbes, has been successfully implemented as these production hosts not only produce the precursor of triterpenoids 2,3-oxidosqualene by the mevalonate pathway but also allow simple expression of plant membrane-anchored enzymes. Nevertheless, 2,3-oxidosqualene is natively converted to lanosterol catalyzed by the endogenous lanosterol synthase (Erg7p), causing low production of recombinant triterpenoids. While simple deletion of ERG7 was not effective, in this study, the critical amino acid residues of Erg7p were engineered to lower this critical enzyme activity. The engineered S. cerevisiae indeed accumulated 2,3-oxidosqualene up to 180 mg/L. Engineering triterpenoid synthesis into the ERG7-modified strain resulted in 7.3- and 3-fold increases in the titers of dammarane-type and lupane-type triterpenoids, respectively. This study presents an efficient inducer-free strategy for lowering Erg7p activity, thereby providing 2,3-oxidosqualene for the enhanced production of various triterpenoids.

Engineering of Yarrowia lipolytica for the production of plant triterpenoids: Asiatic, madecassic, and arjunolic acids

Several plant triterpenoids have valuable pharmaceutical properties, but their production and usage is limited since extraction from plants can burden natural resources, and result in low yields and purity. Here, we engineered oleaginous yeast Yarrowia lipolytica to produce three valuable plant triterpenoids (asiatic, madecassic, and arjunolic acids) by fermentation. First, we established the recombinant production of precursors, ursolic and oleanolic acids, by expressing plant enzymes in free or fused versions in a Y. lipolytica strain previously optimized for squalene production. Engineered strains produced up to 11.6 mg/g DCW ursolic acid or 10.2 mg/g DCW oleanolic acid. The biosynthetic pathway from ursolic acid was extended by expressing the Centella asiatica cytochrome P450 monoxygenases CaCYP716C11p, CaCYP714E19p, and CaCYP716E41p, resulting in the production of trace amounts of asiatic acid and 0.12 mg/g DCW madecassic acid. Expressing the same C. asiatica cytochromes P450 in oleanolic acid-producing strain resulted in the production of oleanane triterpenoids. Expression of CaCYP716C11p in the oleanolic acid-producing strain yielded 8.9 mg/g DCW maslinic acid. Further expression of a codon-optimized CaCYP714E19p resulted in 4.4 mg/g DCW arjunolic acid. Lastly, arjunolic acid production was increased to 9.1 mg/g DCW by swapping the N-terminal domain of CaCYP714E19p with the N-terminal domain from a Kalopanax septemlobus cytochrome P450. In summary, we have demonstrated the production of asiatic, madecassic, and arjunolic acids in a microbial cell factory. The strains and fermentation processes need to be further improved before the production of these molecules by fermentation can be industrialized.

Genome-wide transcriptome analysis and characterization of the cytochrome P450 flavonoid biosynthesis genes in pigeon pea (Cajanus cajan)

226 CcCYP450 genes were identified at the genomic level and were classified into 45 clades based on phylogenetic analysis. CcCYP75B165 gene was found that might play important roles in the biosynthesis of flavonoids in pigeon pea, and was significantly induced by methyl jasmonate (MeJA). The cytochrome P450 mono-oxygenase (CYP450) superfamily plays a key role in the flavonoid biosynthesis pathway and resists different kinds of stresses. Several CYP450 genes have been identified to be involved in the biosynthesis of crop protection agents. However, the CcCYP450 genes from pigeon pea have not been identified. Here, 226 CcCYP450 genes were identified at the genomic level by analysing the gene structure, distribution on chromosomes, gene duplication, and conserved motifs and were classified into 45 clades based on phylogenetic analysis. RNA-seq analysis revealed clear details of CcCYP450 genes that varied with time of MeJA (methyl jasmonate) induction. Among them, six CcCYP450 subfamily genes were found that might play important roles in the biosynthesis of flavonoids in pigeon pea. The overexpression of CcCYP75B165 in pigeon pea significantly induced the accumulation of genistin and downregulated the contents of cajaninstilbene acid, apigenin, isovitexin, and genistein and the expression of flavonoid synthase genes. This study provides theoretical guidance and plant genetic resources for cultivating new pigeon pea varieties with high flavonoid contents and exploring the molecular mechanisms of the biosynthesis of flavonoids under MeJA treatment.

A High-Quality Reference Genome Sequence and Genetic Transformation System of Aralia elata

Aralia elata is a perennial woody plant of the genus Aralia in the family Araliaceae. It is rich in saponins and therefore has a wide range of pharmacological effects. Here, we report a high-quality reference genome of A. elata, with a genome size of 1.21 Gb and a contig N50 of 51.34 Mb, produced by PacBio HiFi sequencing technology. This is the first genome assembly for the genus Aralia. Through genome evolutionary analysis, we explored the phylogeny and whole genome duplication (WGD) events in the A. elata genome. The results indicated that a recent WGD event occurred in the A. elata genome. Estimation of the divergence times indicated that the WGD may be shared by Araliaceae. By analyzing the genome sequence of A. elata and combining the transcriptome data from three tissues, we discovered important genes related to triterpene saponins biosynthesis. Furthermore, based on the embryonic callus induction system of A. elata established in our laboratory, we set up the genetic transformation system of this plant. The genomic resources and genetic transformation system obtained in this study provide insights into A. elata and lays the foundation for further exploration of the A. elata regulatory mechanism.

High-yield bioactive triterpenoid production by heterologous expression in Nicotiana benthamiana using the Tsukuba system

Oleanolic acid is a pentacyclic triterpenoid found in numerous plant species and is a precursor to several bioactive triterpenoids with commercial potential. However, oleanolic acid accumulates at low levels in plants, and its chemical synthesis is challenging. Here, we established a method for producing oleanolic acid in substantial quantities via heterologous expression of pathway enzymes in Nicotiana benthamiana. The "Tsukuba system" is one of the most efficient agroinfiltration-based transient protein expression systems using the vector pBYR2HS, which contains geminiviral replication machinery and a double terminator for boosting expression. Additionally, the pBYR2HS vector contains an expression cassette for the gene-silencing suppressor p19 protein from tomato bushy stunt virus, which can also contribute to enhancing the expression of target proteins. In this study, we evaluated the applicability of this system to heterologous triterpenoid production in N. benthamiana. Medicago truncatula cytochrome P450 monooxygenase (CYP) 716A12 is the first enzyme to be functionally characterized as β-amyrin C-28 oxidase producing oleanolic acid. A mutant CYP716A12 (D122Q) with improved catalytic activity engineered in our previous study was co-expressed with other enzymes in N. benthamiana leaves. Using pBYR2HS, oleanolic acid yield was increased 13.1-fold compared with that using the conventional binary vector, indicating the advantage of the Tsukuba system. We also demonstrated the efficacy of co-expressing a mutant Arabidopsis thaliana HMGR1 catalytic domain, additional NADPH-cytochrome P450 reductase (CPR) transferring electrons to heterologous CYPs, and application of ascorbic acid for preventing leaf necrosis after agroinfiltration, to improve product yield. As a result, the product yields of both simple (β-amyrin) and oxidized (oleanolic acid and maslinic acid) triterpenoids were significantly improved compared with the previously reported yield in heterologous triterpenoid production in N. benthamiana leaves.

Genome-Wide Analysis of the Cytochrome P450 Gene Family Involved in Salt Tolerance in Gossypium hirsutum

Plant cytochrome P450 (P450) participates in a wide range of biosynthetic reactions and targets a variety of biological molecules. These reactions lead to various fatty acid conjugates, plant hormones, secondary metabolites, lignin, and various defensive compounds. In our previous research, transcriptome analysis was performed on the salt-tolerant upland cotton "Tongyan No. 1." Many differentially expressed genes (DEGs) belong to the P450 family, and their domains occur widely in plants. In this current research, P450 genes were identified in Gossypium hirsutum with the aid of bioinformatics methods for investigating phylogenetic relations, gene structure, cis-elements, chromosomal localization, and collinearity within a genome. qRT-PCR was conducted to analyze P450 gene expression patterns under salt stress. The molecular weights of the 156 P450 genes were in the range of 5,949.6-245,576.3 Da, and the length of the encoded amino acids for all the identified P450 genes ranged from 51 to 2,144. P450 proteins are divided into four different subfamilies based on phylogenetic relationship, gene structure, and chromosomal localization of gene replication. The length of P450 genes in upland cotton differs greatly, ranging from 1,500 to 13,000 bp. The number of exons in the P450 family genes ranged from 1 to 9, while the number of introns ranged from 0 to 8, and there were similar trends within clusters. A total of 31 cis-acting elements were identified by analyzing 1,500 bp promoter sequences. Differences were found in cis-acting elements among genes. The consistency between qRT-PCR and previous transcriptome analysis of salt tolerance DEGs indicated that they were likely to be involved in the salt tolerance of cotton seedlings. Our results provide valuable information on the evolutionary relationships of genes and functional characteristics of the gene family, which is beneficial for further study of the cotton P450 gene family.

Toward the Heterologous Biosynthesis of Plant Natural Products: Gene Discovery and Characterization

Plant natural products (PNPs) represent a vast and diverse group of natural products, which have wide applications such as emulsifiers in cosmetics, sweeteners in foods, and active ingredients in medicines. Large-scale production of certain PNPs (e.g., artemisinin, taxol) has been implemented by reconstruction of biosynthetic pathways in heterologous hosts. However, unknown biosynthetic pathways greatly restrict wide applications of heterologous production of PNPs of interest. With the rapid development of sequencing and multiomics analysis technologies, huge amounts of omics data, i.e., genomics, transcriptomics, and proteomics, have been deposited in public databases, which is a precious resource for identification of the unknown biosynthetic pathway of PNPs. Herein, we have enumerated the approaches which have been widely used to screen candidate genes involved in the biosynthesis of PNPs of interest. We also discuss recent developments in the characterization of putative genes and elucidation of the complete biosynthetic pathway in heterologous hosts.

Genome and transcriptome analysis to understand the role diversification of cytochrome P450 gene under excess nitrogen treatment

BACKGROUND

Panax notoginseng (Burk.) F. H. Chen (P. notoginseng) is a medicinal plant. Cytochrome P450 (CYP450) monooxygenase superfamily is involved in the synthesis of a variety of plant hormones. Studies have shown that CYP450 is involved in the synthesis of saponins, which are the main medicinal component of P. notoginseng. To date, the P. notoginseng CYP450 family has not been systematically studied, and its gene functions remain unclear.

RESULTS

In this study, a total of 188 PnCYP genes were identified, these genes were divided into 41 subfamilies and clustered into 9 clans. Moreover, we identified 40 paralogous pairs, of which only two had Ka/Ks ratio greater than 1, demonstrating that most PnCYPs underwent purification selection during evolution. In chromosome mapping and gene replication analysis, 8 tandem duplication and 11 segmental duplication events demonstrated that PnCYP genes were continuously replicating during their evolution. Gene ontology (GO) analysis annotated the functions of 188 PnCYPs into 21 functional subclasses, suggesting the functional diversity of these gene families. Functional divergence analyzed the members of the three primitive branches of CYP51, CYP74 and CYP97 at the amino acid level, and found some critical amino acid sites. The expression pattern of PnCYP450 related to nitrogen treatment was studied using transcriptome sequencing data, 10 genes were significantly up-regulated and 37 genes were significantly down-regulated. Combined with transcriptome sequencing analysis, five potential functional genes were screened. Quantitative real-time PCR (qRT-PCR) indicated that these five genes were responded to methyl jasmonate (MEJA) and abscisic acid (ABA) treatment.

CONCLUSIONS

These results provide a valuable basis for comprehending the classification and biological functions of PnCYPs, and offer clues to study their biological functions in response to nitrogen treatment.

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.

β-Amyrin Synthase1 Controls the Accumulation of the Major Saponins Present in Pea (Pisum sativum)

The use of pulses as ingredients for the production of food products rich in plant proteins is increasing. However, protein fractions prepared from pea or other pulses contain significant amounts of saponins, glycosylated triterpenes that can impart an undesirable bitter taste when used as an ingredient in foodstuffs. In this article, we describe the identification and characterization of a gene involved in saponin biosynthesis during pea seed development, by screening mutants obtained from two Pisum sativum TILLING (Targeting Induced Local Lesions IN Genomes) populations in two different genetic backgrounds. The mutations studied are located in a gene designated PsBAS1 (β-amyrin synthase1), which is highly expressed in maturing pea seeds and which encodes a protein previously shown to correspond to an active β-amyrin synthase. The first allele is a nonsense mutation, while the second mutation is located in a splice site and gives rise to a mis-spliced transcript encoding a truncated, nonfunctional protein. The homozygous mutant seeds accumulated virtually no saponin without affecting the seed nutritional or physiological quality. Interestingly, BAS1 appears to control saponin accumulation in all other tissues of the plant examined. These lines represent a first step in the development of pea varieties lacking bitterness off-flavors in their seeds. Our work also shows that TILLING populations in different genetic backgrounds represent valuable genetic resources for both crop improvement and functional genomics.

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.

Bacterial Skin Infections in Livestock and Plant-Based Alternatives to Their Antibiotic Treatment

Due to its large surface area, the skin is susceptible to various injuries, possibly accompanied by the entrance of infective agents into the body. Commensal organisms that constitute the skin microbiota play important roles in the orchestration of cutaneous homeostasis and immune competence. The opportunistic pathogen Staphylococcus aureus is present as part of the normal biota of the skin and mucous membranes in both humans and animals, but can cause disease when it invades the body either due to trauma or because of the impaired immune response of the host. Colonization of livestock skin by S. aureus is a precursor for majority of bacterial skin infections, which range from boils to sepsis, with the best-characterized being bovine mastitis. Antibiotic treatment of these infections can contribute to the promotion of resistant bacterial strains and even to multidrug resistance. The development of antibiotic resistance to currently available antibiotics is a worldwide problem. Considering the increasing ability of bacteria to effectively resist antibacterial agents, it is important to reduce the livestock consumption of antibiotics to preserve antibiotic effectiveness in the future. Plants are recognized as sources of various bioactive substances, including antibacterial activity towards clinically important microorganisms. This review provides an overview of the current knowledge on the major groups of phytochemicals with antibacterial activity and their modes of action. It also provides a list of currently known and used plant species aimed at treating or preventing bacterial skin infections in livestock.

CRISPR/Cas9-Mediated Targeted Mutagenesis of CYP93E2 Modulates the Triterpene Saponin Biosynthesis in Medicago truncatula

In the Medicago genus, triterpene saponins are a group of bioactive compounds extensively studied for their different biological and pharmaceutical properties. In this work, the CRISPR/Cas9-based approach with two single-site guide RNAs was used in Medicago truncatula (barrel medic) to knock-out the CYP93E2 and CYP72A61 genes, which are responsible for the biosynthesis of soyasapogenol B, the most abundant soyasapogenol in Medicago spp. No transgenic plants carrying mutations in the target CYP72A61 gene were recovered while fifty-two putative CYP93E2 mutant plant lines were obtained following Agrobacterium tumefaciens-mediated transformation. Among these, the fifty-one sequenced plant lines give an editing efficiency of 84%. Sequencing revealed that these lines had various mutation patterns at the target sites. Four T0 mutant plant lines were further selected and examined for their sapogenin content and plant growth performance under greenhouse conditions. The results showed that all tested CYP93E2 knock-out mutants did not produce soyasapogenols in the leaves, stems and roots, and diverted the metabolic flux toward the production of valuable hemolytic sapogenins. No adverse influence was observed on the plant morphological features of CYP93E2 mutants under greenhouse conditions. In addition, differential expression of saponin pathway genes was observed in CYP93E2 mutants in comparison to the control. Our results provide new and interesting insights into the application of CRISPR/Cas9 for metabolic engineering of high-value compounds of plant origin and will be useful to investigate the physiological functions of saponins in planta.

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.

Ginsenosides in Panax genus and their biosynthesis

Ginsenosides are a series of glycosylated triterpenoids which belong to protopanaxadiol (PPD)-, protopanaxatriol (PPT)-, ocotillol (OCT)- and oleanane (OA)-type saponins known as active compounds of Panax genus. They are accumulated in plant roots, stems, leaves, and flowers. The content and composition of ginsenosides are varied in different ginseng species, and in different parts of a certain plant. In this review, we summarized the representative saponins structures, their distributions and the contents in nearly 20 Panax species, and updated the biosynthetic pathways of ginsenosides focusing on enzymes responsible for structural diversified ginsenoside biosynthesis. We also emphasized the transcription factors in ginsenoside biosynthesis and non-coding RNAs in the growth of Panax genus plants, and highlighted the current three major biotechnological applications for ginsenosides production. This review covered advances in the past four decades, providing more clues for chemical discrimination and assessment on certain ginseng plants, new perspectives for rational evaluation and utilization of ginseng resource, and potential strategies for production of specific ginsenosides.

Plant cytochrome P450s: Role in stress tolerance and potential applications for human welfare

Cytochrome P450s (CYPs) are a versatile group of enzymes and one of the largest families of proteins, controlling various physiological processes via biosynthetic and detoxification pathways. CYPs perform multiple roles through a critical irreversible enzymatic reaction in which an oxygen atom is inserted within hydrophobic molecules, converting them into the reactive and hydro soluble components. During evolution, plants have acquired significantly more number of CYPs and represent about 1% of the encoded genes . CYPs are highly conserved proteins involved in growth, development and tolerance against biotic and abiotic stresses. Furthermore, CYPs reinforce plants' molecular and chemical defense mechanisms by regulating the biosynthesis of secondary metabolites, enhancing reactive oxygen species (ROS) scavenging and controlling biosynthesis and homeostasis of phytohormones, including abscisic acid (ABA) and jasmonates. Thus, they are the critical targets of metabolic engineering for enhancing plant defense against environmental stresses. Additionally, CYPs are also used as biocatalysts in the fields of pharmacology and phytoremediation. Herein, we highlight the role of CYPs in plant stress tolerance and their applications for human welfare.

Cleaning the Medicago Microarray Database to Improve Gene Function Analysis

Transcriptomics studies have been facilitated by the development of microarray and RNA-Seq technologies, with thousands of expression datasets available for many species. However, the quality of data can be highly variable, making the combined analysis of different datasets difficult and unreliable. Most of the microarray data for Medicago truncatula, the barrel medic, have been stored and made publicly accessible on the web database Medicago truncatula Gene Expression atlas (MtGEA). The aim of this work is to ameliorate the quality of the MtGEA database through a general method based on logical and statistical relationships among parameters and conditions. The initial 716 columns available in the dataset were reduced to 607 by evaluating the quality of data through the sum of the expression levels over the entire transcriptome probes and Pearson correlation among hybridizations. The reduced dataset shows great improvements in the consistency of the data, with a reduction in both false positives and false negatives resulting from Pearson correlation and GO enrichment analysis among genes. The approach we used is of general validity and our intent is to extend the analysis to other plant microarray databases.

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.

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

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

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.

Transcriptome analysis of Paris polyphylla var. yunnanensis illuminates the biosynthesis and accumulation of steroidal saponins in rhizomes and leaves

Paris polyphylla var. yunnanensis can synthesize Paris saponins with multiple effective therapies, and its rhizome has become an indispensable ingredient in many patented drugs. However, how Paris saponin content changes in tissues at different stages and the molecular mechanisms underlying the production and accumulation of the bioactive compounds are unclear. This study aimed to uncover the mechanisms underlying the biosynthesis and accumulation by integrating transcriptome sequencing and phytochemical investigation of the leaves and rhizomes at different growth stages. Paris saponin content in leaves was lower during the fruiting stage than the vegetative stage, whereas the content in rhizomes increased during the fruiting stage. The candidate genes related to Paris saponin biosynthesis were determined by transcriptome analyses. Most biosynthetic genes were found to be abundantly expressed in the leaves during the vegetative stage in the light of expression profiles and functional enrichment results. The expression patterns of the differentially expressed genes related to the biosynthesis were positively correlated with the accumulation of saponins in tissues. These findings suggest that both leaves and rhizomes are capable of biosynthesizing Paris saponins, and that aerial plant parts can be used to extract them. The different patterns of biosynthesis and accumulation in the leaves and rhizomes were also determined here. This study will help improve our understanding of the mechanisms underlying the biosynthesis and accumulation of Paris saponins, and aid in the comprehensive development and utilization of this medicinal plant.

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.

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.

Triterpenic saponins from Medicago marina L

The saponin composition of leaves and roots from Medicago marina L., sea medic, was investigated by a combination of chromatographic, spectroscopic and spectrometric (GC, LC, ESI-MS/MS, NMR) methods. Several compounds were detected and quantified by HPLC using the external standard method. Saponins from this plant species consist of a mixture of high molecular weight bidesmosidic derivatives of medicagenic and zanhic acid, containing up to six sugars in the molecules. Six of the detected saponins were previously isolated and reported as constituents of other Medicago spp.; one saponin was previously described in other plant species; four saponins are undescribed compounds in Medicago and never reported before in other plant species. These are: 3-O-β-D-glucopyranosyl-(1 → 2)-β-D-glucopyranosylzanhic acid 28-O-β-D-xylopyranosyl-(1 → 4)-[β-D-apiofuranosyl-(1 → 3)]-α-L-rhamnopyranosyl-(1 → 2)-α-L-arabinopyranosyl ester; 3-O-β-D-glucopyranosyl-(1 → 2)-β-D-glucopyranosylzanhic acid 28-O-β-D-xylopyranosyl-(1 → 4)-α-L-rhamnopyranosyl-(1 → 2)-α-L-arabinopyranosyl ester; 3-O-β-D-glucopyranosyl-(1 → 2)-β-D-glucopyranosylmedicagenic acid 28-O-β-D-xylopyranosyl-(1 → 4)-[α-L-arabinopyranosyl-(1 → 3)]-α-L-rhamnopyranosyl-(1 → 2)-α-L-arabinopyranosyl ester and 3-O-β-D-glucopyranosylmedicagenic acid 28-O-β-D-xylopyranosyl-(1 → 4)-[α-L-arabinopyranosyl-(1 → 3)]-α-L-rhamnopyranosyl-(1 → 2)-α-L-arabinopyranosyl ester. The specific saponins synthesized by M. marina may have a role in its tolerance to environment, representing a reservoir of osmolytic sugars.

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.