Novel triterpene oxidizing activity ofArabidopsis thalianaCYP716A subfamily enzymes

Genomic investigation of plant secondary metabolism: insights from synteny network analysis of oxidosqualene cyclase flanking genes

The clustered distribution of genes involved in metabolic pathways within the plant genome has garnered significant attention from researchers. By comparing and analyzing changes in the flanking regions of metabolic genes across a diverse array of species, we can enhance our understanding of the formation and distribution of biosynthetic gene clusters (BGCs). In this study, we have designed a workflow that uncovers and assesses conserved positional relationships between genes in various species by using synteny neighborhood networks (SNN). This workflow is then applied to the analysis of flanking genes associated with oxidosqualene cyclases (OSCs). The method allows for the recognition and comparison of homologous blocks with unique flanking genes accompanying different subfamilies of OSCs. The examination of the flanking genes of OSCs in 122 plant species revealed multiple genes with conserved positional relationships with OSCs in angiosperms. Specifically, the earliest adjacency of OSC genes and CYP716 genes first appeared in basal eudicots, and the nonrandom occurrence of CYP716 genes in the flanking region of OSC persists across different lineages of eudicots. Our study showed the substitution of genes in the flanking region of the OSC varies across different plant lineages, and our approach facilitates the investigation of flanking gene rearrangements in the formation of OSC-related BGCs.

Synergistic interaction between wheat streak mosaic virus and Triticum mosaic virus modulates wheat transcriptome to favor disease severity

Wheat streak mosaic virus (WSMV; Tritimovirus tritici) and Triticum mosaic virus (TriMV; Poacevirus tritici), the type members of the genera Tritimovirus and Poacevirus, respectively, in the family Potyviridae, are economically important wheat viruses in the Great Plains region of the USA. Co-infection of wheat by WSMV and TriMV results in disease synergism. Wheat transcriptome from singly (WSMV or TriMV) and doubly (WSMV+TriMV) infected upper uninoculated leaves were analyzed by RNA-Seq at 9, 12, and 21 days postinoculation. A total of 31,754 differentially expressed wheat genes were identified among all comparisons. Weighted gene co-expression network analysis resulted in 11 co-expression modules that broadly indicated gene expression profiles attributable to control, single, and double infections. Gene ontology, protein domain and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis revealed that genes specifically related to photosynthesis, growth, stress, senescence, and defense were differentially enriched. Analyses of transcription factor families indicated that genes encoding MADS-Box and ARFs were strongly enriched in control plants, moderately repressed in TriMV-infected plants, and more strongly repressed in WSMV- and doubly-infected plants, whereas genes encoding WRKYs and NACs were more enriched in WSMV or doubly infected plants. Synergistic interactions between WSMV and TriMV drastically enhanced disease phenotype compared to individual virus infections. The progression of disease phenotype was correlated to transcriptomic changes, indicating the strong disruption to plant metabolism and likely channeling of energy and metabolites for viral replication. There also appeared to be a connection between viral replication and plastid health, with stronger downregulation of genes needed for chloroplast functions and integrity and increased synergism between TriMV and WSMV. This study provides an overview of transcriptomic changes distinctly influenced by TriMV and WSMV either singly or in combination and provides a good correlation between specific transcription factors and genes associated with metabolism to observed phenotypic changes in plant growth and disease synergism.

Integrated metabolomic and transcriptomic analysis revealed the role of PbrCYP94B in wax accumulation in pear fruit after bagging treatment

Preharvest fruit bagging is a safe and environmentally friendly production measure. Cuticular wax, as the first protective layer on the fruit surface, has important functions. However, the effects of preharvest bagging on cuticular wax synthesis in pears and the related molecular mechanisms are still unclear. Here, the impact of fruit bagging with different materials on cuticular wax synthesis in pear fruit, and the underlying molecular mechanism, were revealed from metabolomic, transcriptomic, morphological, and molecular biological perspectives. Our results revealed that, compared with that in the not bagged (NB) treatment group (0.59 mg/cm2), the total wax concentration was 1.32- and 1.37-fold greater in the single-layered white paper bag (WPB, 1.37 mg/cm2) and double-layered yellow-white paper bag, (YWPB, 1.40 mg/cm2) treatment groups, while it was slightly lower in the double-layered yellow-black paper bag (YBPB, 0.45 mg/cm2) group, which was consistent with the scanning electron microscopy (SEM) results. Integrated metabolomic and transcriptomic analysis revealed 29 genes associated with cuticular wax synthesis. Overexpression of PbrCYP94B, which is a key gene in the wax synthesis pathway in pear fruit, increased the total wax and alkane contents. This study provides valuable insights for the creation of new pear germplasms with high wax contents.

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.

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.

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.

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

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

Gene identification and semisynthesis of the anti-inflammatory oleanane-type triterpenoid wilforlide A

Wilforlide A is one of the main active constituents produced in trace amounts in Tripterygium wilfordii Hook F, which has excellent anti-inflammatory and immune suppressive effects. Despite the seeming structural simplicity of the compound, the biosynthetic pathway of wilforlide A remains unknown. Gene-specific expression analysis and genome mining were used to identify the gene candidates, and their functions were studied in vitro and in vivo. A modularized two-step (M2S) technique and CRISPR-Cas9 methods were used to construct engineering yeast. Here, we identified a cytochrome P450, TwCYP82AS1, that catalyses C-22 hydroxylation during wilforlide A biosynthesis. We also found that TwCYP712K1 to K3 can further oxidize the C-29 carboxylation of oleanane-type triterpenes in addition to friedelane-type triterpenes. Reconstitution of the biosynthetic pathway in engineered yeast increased the precursor supply, and combining TwCYP82AS1 and TwCYP712Ks produced abrusgenic acid, which was briefly acidified to achieve the semisynthesis of wilforlide A. Our work presents an alternative metabolic engineering approach for obtaining wilforlide A without relying on extraction from plants.

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.

Hydroxylases involved in terpenoid biosynthesis: a review

Terpenoids are pervasive in nature and display an immense structural diversity. As the largest category of plant secondary metabolites, terpenoids have important socioeconomic value in the fields of pharmaceuticals, spices, and food manufacturing. The biosynthesis of terpenoid skeletons has made great progress, but the subsequent modifications of the terpenoid framework are poorly understood, especially for the functionalization of inert carbon skeleton usually catalyzed by hydroxylases. Hydroxylase is a class of enzymes that plays an important role in the modification of terpenoid backbone. This review article outlines the research progress in the identification, molecular modification, and functional expression of this class of enzymes in the past decade, which are profitable for the discovery, engineering, and application of more hydroxylases involved in the plant secondary metabolism.

Whole-genome Duplication Reshaped Adaptive Evolution in A Relict Plant Species, Cyclocarya paliurus

Cyclocarya paliurus is a relict plant species that survived the last glacial period and shows a population expansion recently. Its leaves have been traditionally used to treat obesity and diabetes with the well-known active ingredient cyclocaric acid B. Here, we presented three C. paliurus genomes from two diploids with different flower morphs and one haplotype-resolved tetraploid assembly. Comparative genomic analysis revealed two rounds of recent whole-genome duplication events and identified 691 genes with dosage effects that likely contribute to adaptive evolution through enhanced photosynthesis and increased accumulation of triterpenoids. Resequencing analysis of 45 C. paliurus individuals uncovered two bottlenecks, consistent with the known events of environmental changes, and many selectively swept genes involved in critical biological functions, including plant defense and secondary metabolite biosynthesis. We also proposed the biosynthesis pathway of cyclocaric acid B based on multi-omics data and identified key genes, in particular gibberellin-related genes, associated with the heterodichogamy in C. paliurus species. Our study sheds light on evolutionary history of C. paliurus and provides genomic resources to study the medicinal herbs.

Genomic, transcriptomic, and metabolomic analysis of Oldenlandia corymbosa reveals the biosynthesis and mode of action of anti-cancer metabolites

Plants accumulate a vast array of secondary metabolites, which constitute a natural resource for pharmaceuticals. Oldenlandia corymbosa belongs to the Rubiaceae family, and has been used in traditional medicine to treat different diseases, including cancer. However, the active metabolites of the plant, their biosynthetic pathway and mode of action in cancer are unknown. To fill these gaps, we exposed this plant to eight different stress conditions and combined different omics data capturing gene expression, metabolic profiles, and anti-cancer activity. Our results show that O. corymbosa extracts are active against breast cancer cell lines and that ursolic acid is responsible for this activity. Moreover, we assembled a high-quality genome and uncovered two genes involved in the biosynthesis of ursolic acid. Finally, we also revealed that ursolic acid causes mitotic catastrophe in cancer cells and identified three high-confidence protein binding targets by Cellular Thermal Shift Assay (CETSA) and reverse docking. Altogether, these results constitute a valuable resource to further characterize the biosynthesis of active metabolites in the Oldenlandia group, while the mode of action of ursolic acid will allow us to further develop this valuable compound.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

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

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

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

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

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

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

The lncRNA MARS modulates the epigenetic reprogramming of the marneral cluster in response to ABA

Clustered organization of biosynthetic non-homologous genes is emerging as a characteristic feature of plant genomes. The co-regulation of clustered genes seems to largely depend on epigenetic reprogramming and three-dimensional chromatin conformation. In this study, we identified the long non-coding RNA (lncRNA) MARneral Silencing (MARS), localized inside the Arabidopsis marneral cluster, which controls the local epigenetic activation of its surrounding region in response to abscisic acid (ABA). MARS modulates the POLYCOMB REPRESSIVE COMPLEX 1 (PRC1) component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) binding throughout the cluster in a dose-dependent manner, determining H3K27me3 deposition and chromatin condensation. In response to ABA, MARS decoys LHP1 away from the cluster and promotes the formation of a chromatin loop bringing together the MARNERAL SYNTHASE 1 (MRN1) locus and a distal ABA-responsive enhancer. The enrichment of co-regulated lncRNAs in clustered metabolic genes in Arabidopsis suggests that the acquisition of novel non-coding transcriptional units may constitute an additional regulatory layer driving the evolution of biosynthetic pathways.

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.

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 Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

Comparative Analysis of NADPH-Cytochrome P450 Reductases From Legumes for Heterologous Production of Triterpenoids in Transgenic Saccharomyces cerevisiae

Triterpenoids are plant specialized metabolites with various pharmacological activities. They are widely distributed in higher plants, such as legumes. Because of their low accumulation in plants, there is a need for improving triterpenoid production. Cytochrome P450 monooxygenases (CYPs) play critical roles in the structural diversification of triterpenoids. To perform site-specific oxidations, CYPs require the electrons that are transferred by NADPH-cytochrome P450 reductase (CPR). Plants possess two main CPR classes, class I and class II. CPR classes I and II have been reported to be responsible for primary and specialized (secondary) metabolism, respectively. In this study, we first analyzed the CPR expression level of three legumes species, Medicago truncatula, Lotus japonicus, and Glycyrrhiza uralensis, showing that the expression level of CPR class I was lower and more stable, while that of CPR class II was higher in almost all the samples. We then co-expressed different combinations of CYP716As and CYP72As with different CPR classes from these three legumes in transgenic yeast. We found that CYP716As worked better with CPR-I from the same species, while CYP72As worked better with any CPR-IIs. Using engineered yeast strains, CYP88D6 paired with class II GuCPR produced the highest level of 11-oxo-β-amyrin, the important precursor of high-value metabolites glycyrrhizin. This study provides insight into co-expressing genes from legumes for heterologous production of triterpenoids in yeast.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Integrating pathway elucidation with yeast engineering to produce polpunonic acid the precursor of the anti-obesity agent celastrol

BACKGROUND

Celastrol is a promising anti-obesity agent that acts as a sensitizer of the protein hormone leptin. Despite its potent activity, a sustainable source of celastrol and celastrol derivatives for further pharmacological studies is lacking.

RESULTS

To elucidate the celastrol biosynthetic pathway and reconstruct it in Saccharomyces cerevisiae, we mined a root-transcriptome of Tripterygium wilfordii and identified four oxidosqualene cyclases and 49 cytochrome P450s as candidates to be involved in the early steps of celastrol biosynthesis. Using functional screening of the candidate genes in Nicotiana benthamiana, TwOSC4 was characterized as a novel oxidosqualene cyclase that produces friedelin, the presumed triterpenoid backbone of celastrol. In addition, three P450s (CYP712K1, CYP712K2, and CYP712K3) that act downstream of TwOSC4 were found to effectively oxidize friedelin and form the likely celastrol biosynthesis intermediates 29-hydroxy-friedelin and polpunonic acid. To facilitate production of friedelin, the yeast strain AM254 was constructed by deleting UBC7, which afforded a fivefold increase in friedelin titer. This platform was further expanded with CYP712K1 to produce polpunonic acid and a method for the facile extraction of products from the yeast culture medium, resulting in polpunonic acid titers of 1.4 mg/L.

CONCLUSION

Our study elucidates the early steps of celastrol biosynthesis and paves the way for future biotechnological production of this pharmacologically promising compound in engineered yeast strains.

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

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

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

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

Discovery and Engineering of Cytochrome P450s for Terpenoid Biosynthesis

Terpenoids represent 60% of known natural products, including many drugs and drug candidates, and their biosynthesis is attracting great interest. However, the unknown cytochrome P450s (CYPs) in terpenoid biosynthetic pathways make the heterologous production of related terpenoids impossible, while the slow kinetics of some known CYPs greatly limit the efficiency of terpenoid biosynthesis. Thus, there is a compelling need to discover and engineer CYPs for terpenoid biosynthesis to fully realize their great potential for industrial application. This review article summarizes the current state of CYP discovery and engineering in terpenoid biosynthesis, focusing on recent synthetic biology approaches toward prototyping CYPs in heterologous hosts. We also propose several strategies for further accelerating CYP discovery and engineering.

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.

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

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

A Novel Multifunctional C-23 Oxidase, CYP714E19, is Involved in Asiaticoside Biosynthesis

Centella asiatica is widely used as a medicinal plant due to accumulation of the ursane-type triterpene saponins asiaticoside and madecassoside. The molecular structure of both compounds suggests that they are biosynthesized from α-amyrin via three hydroxylations, and the respective Cyt P450-dependent monooxygenases (P450 enzymes) oxidizing the C-28 and C-2α positions have been reported. However, a third enzyme hydroxylating C-23 remained elusive. We previously identified 40,064 unique sequences in the transcriptome of C. asiatica elicited by methyl jasmonate, and among them we have now found 149 unigenes encoding putative P450 enzymes. In this set, 23 full-length cDNAs were recognized, 13 of which belonged to P450 subfamilies previously implicated in secondary metabolism. Four of these genes were highly expressed in response to jasmonate treatment, especially in leaves, in accordance with the accumulation patterns of asiaticoside. The functions of these candidate genes were tested using heterologous expression in yeast cells. Gas chromatography-mass spectrometry (GC-MS) analysis revealed that yeast expressing only the oxidosqualene synthase CaDDS produced the asiaticoside precursor α-amyrin (along with its isomer β-amyrin), while yeast co-expressing CaDDS and CYP716A83 also contained ursolic acid along with oleanolic acid. This P450 enzyme thus acts as a multifunctional triterpenoid C-28 oxidase converting amyrins into corresponding triterpenoid acids. Finally, yeast strains co-expressing CaDDS, CYP716A83 and CYP714E19 produced hederagenin and 23-hydroxyursolic acid, showing that CYP714E19 is a multifunctional triterpenoid oxidase catalyzing the C-23 hydroxylation of oleanolic acid and ursolic acid. Overall, our results demonstrate that CaDDS, CYP716A83 and CYP714E19 are C. asiatica enzymes catalyzing consecutive steps in asiaticoside biosynthesis.

Identification of candidate genes mediating apple fruit-cracking resistance following the application of gibberellic acids 4 + 7 and the cytokinin 6-benzyladenine

Calyx-end cracking in 'Pink Lady' apple is treated by a solution of gibberellic acids 4 and 7 (GA₄₊₇) and the cytokinin 6-benzyladenine (BA). Although the GA₄₊₇ and BA mixture is applied early in apple fruit development, it mitigates cracking that becomes evident in the mature fruit, implying a long-term treatment effect. The reduced incidence of peel cracking is associated with increased epidermal cell density, which is maintained until fruit maturation. Presently, the expression of genes that have been previously reported to be associated with epidermal cell patterning and cuticle formation, or cracking resistance, was monitored in the peel during fruit development and following GA₄₊₇ and BA treatment. For most of the genes whose expression is naturally upregulated during fruit development, the early GA₄₊₇ and BA treatment maintained or further increased the high expression level in the mature peel. Where the expression of a gene was downregulated during development, no change was detected in the treated mature peel. Gene-networking analysis supported the interaction between gene clusters of cell-wall synthesis, cuticle formation and GA signaling. Overall, the data suggested that the GA₄₊₇ and BA treatment did not modify developmental cues, but promoted or enhanced the innate developmental program.

Transcriptomic Analysis of Kalopanax septemlobus and Characterization of KsBAS, CYP716A94 and CYP72A397 Genes Involved in Hederagenin Saponin Biosynthesis

Kalopanax septemlobus, commonly named the castor aralia tree, is a highly valued woody medicinal tree belonging to the family Araliaceae. Kalopanax septemlobus contains approximately 15 triterpenoid saponins primarily constituted of hederagenin aglycones. Hederagenin is a representative precursor for hemolytic saponin in plants. In the present study, transcriptome analysis was performed to discover genes involved in hederagenin saponin biosynthesis in K. septemlobus. De novo assembly generated 82,698 unique sequences, including 17,747 contigs and 64,951 singletons, following 454 pyrosequencing. Oxidosqualene cyclases (OSCs) are enzymes that catalyze the formation of diverse triterpene skeletons from 2,3-oxidosqualene. Heterologous expression of an OSC sequence in yeast revealed that KsBAS is a β-amyrin synthase gene. Cytochrome P450 genes (CYPs) make up a supergene family in the plant genome and play a key role in the biosynthesis of sapogenin aglycones. In total, 95 contigs and 110 singletons annotated as CYPs were obtained by sequencing the K. septemlobus transcriptome. By heterologous expression in yeast, we found that CYP716A94 was β-amyrin 28-oxidase involved in oleanolic acid production from β-amyrin, and CYP72A397 was oleanolic acid 23-hydroxylase involved in hederagenin production from oleanolic acid. Engineered yeast co-expressing KsBAS, CYP716A94 and CYP72A397 produced hederagenin. Kalopanax septemlobus CYP72A397 is a novel CYP enzyme that synthesizes hederagenin aglycone from oleanolic acid as a single product. In conclusion, we characterized three genes participating in sequential steps for hederagenin biosynthesis from β-amyrin, which are likely to play a major role in hederagenin saponin biosynthesis in K. septemlobus.

P450s controlling metabolic bifurcations in plant terpene specialized metabolism

ABSTRACT

Catalyzing stereo- and regio-specific oxidation of inert hydrocarbon backbones, and a range of more exotic reactions inherently difficult in formal chemical synthesis, cytochromes P450 (P450s) offer outstanding potential for biotechnological engineering. Plants and their dazzling diversity of specialized metabolites have emerged as rich repository for functional P450s with the advances of deep transcriptomics and genome wide discovery. P450s are of outstanding interest for understanding chemical diversification throughout evolution, for gaining mechanistic insights through the study of their structure-function relationship, and for exploitation in Synthetic Biology. In this review, we highlight recent developments and examples in the discovery of plant P450s involved in the biosynthesis of industrially relevant monoterpenoids, sesquiterpenoids, diterpenoids and triterpenoids, throughout 2016 and early 2017. Examples were selected to illustrate the spectrum of value from commodity chemicals, flavor and fragrance compounds to pharmacologically active terpenoids. We focus on a recently emerging theme, where P450s control metabolic bifurcations and chemical diversity of the final product profile, either within a pathway, or through neo-functionalization in related species. The implications may inform approaches for rational assembly of recombinant pathways, biotechnological production of high value terpenoids and generation of novel chemical entities.

A translational synthetic biology platform for rapid access to gram-scale quantities of novel drug-like molecules

Plants are an excellent source of drug leads. However availability is limited by access to source species, low abundance and recalcitrance to chemical synthesis. Although plant genomics is yielding a wealth of genes for natural product biosynthesis, the translation of this genetic information into small molecules for evaluation as drug leads represents a major bottleneck. For example, the yeast platform for artemisinic acid production is estimated to have taken >150 person years to develop. Here we demonstrate the power of plant transient transfection technology for rapid, scalable biosynthesis and isolation of triterpenes, one of the largest and most structurally diverse families of plant natural products. Using pathway engineering and improved agro-infiltration methodology we are able to generate gram-scale quantities of purified triterpene in just a few weeks. In contrast to heterologous expression in microbes, this system does not depend on re-engineering of the host. We next exploit agro-infection for quick and easy combinatorial biosynthesis without the need for generation of multi-gene constructs, so affording an easy entrée to suites of molecules, some new-to-nature, that are recalcitrant to chemical synthesis. We use this platform to purify a suite of bespoke triterpene analogs and demonstrate differences in anti-proliferative and anti-inflammatory activity in bioassays, providing proof of concept of this system for accessing and evaluating medicinally important bioactives. Together with new genome mining algorithms for plant pathway discovery and advances in plant synthetic biology, this advance provides new routes to synthesize and access previously inaccessible natural products and analogs and has the potential to reinvigorate drug discovery pipelines.

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

The roots of Platycodon grandiflorus are widely used as a crude drug. The active components include a variety of triterpenoid saponins. Recent studies have revealed that Cyt P450 monooxygenases (P450s) function as triterpene oxidases in triterpenoid saponin biosynthesis in many plant species. However, there have been no reports regarding triterpene oxidases in P. grandiflorus. In this study, we performed transcriptome analysis of three different P. grandiflorus tissues (roots, leaves and petals) using RNA sequencing (RNA-Seq) technology. We cloned six P450 genes that were highly expressed in roots, and classified them as belonging to the CYP716A, CYP716D and CYP72A subfamilies. We heterologously expressed these P450s in an engineered yeast strain that produces β-amyrin, one of the most common triterpenes in plants. Two of the CYP716A subfamily P450s catalyzed oxidation reactions of the β-amyrin skeleton. One of these P450s, CYP716A140v2, catalyzed a three-step oxidation reaction at C-28 on β-amyrin to produce oleanolic acid, a reaction performed by CYP716A subfamily P450s in a variety of plant species. The other P450, CYP716A141, catalyzed the hydroxylation of β-amyrin at C-16β. This reaction is unique among triterpene oxidases isolated to date. These results enhance our knowledge of functional variation among CYP716A subfamily enzymes involved in triterpenoid biosynthesis, and provide novel molecular tools for use in synthetic biology to produce triterpenoid saponins with pre-defined structures.

Transcriptomic Comparison Reveals Candidate Genes for Triterpenoid Biosynthesis in Two Closely Related Ilex Species

Native to Southern China, Ilex pubescens and Ilex asprella are frequently used in traditional Chinese medicine. Both of them produce a large variety of ursane-type triterpenoid saponins, which have been demonstrated to have different pharmacological effects. However, little is known about their biosynthesis. In this study, transcriptomic analysis of I. pubescens and comparison with its closely related specie I. asprella were carried out to identify potential genes involved in triterpenoid saponin biosynthesis. Through RNA sequencing (RNA-seq) and de novo transcriptome assembly of I. pubescens, a total of 68,688 UniGene clusters are obtained, of which 32,184 (46.86%) are successfully annotated by comparison with the sequences in major public databases (NCBI, Swiss-Prot, and KEGG). It includes 128 UniGenes related to triterpenoid backbone biosynthesis, 11 OSCs (oxidosqualene cyclases), 233 CYPs (cytochrome P450), and 269 UGTs (UDP-glycosyltransferases). By homology-based blast and phylogenetic analysis with well-characterized genes involved in triterpenoid saponin biosynthesis, 5 OSCs, 14 CYPs, and 1 UGT are further proposed as the most promising candidate genes. Transcriptomic comparison between two Ilex species using blastp and OrthoMCL method reveals high sequence similarity. All OSCs and UGTs as well as most CYPs are classified as orthologous genes, while only 5 CYPs in I. pubescens and 3 CYPs in I. asprella are species-specific. One of OSC candidates, named as IpAS1, was successfully cloned and expressed in Saccharomyces cerevisiae INVSc1. Analysis of the yeast extract by gas chromatography (GC) and gas chromatography-mass spectrometry (GC-MS) shows IpAS1 is a mixed amyrin synthase, producing α-amyrin and β-amyrin at ratio of 5:1, which is similar to its ortholog IaAS1 from I. asprella. This study is the first exploration to profile the transcriptome of I. pubescens, the generated data and gene models will facilitate further molecular studies on the physiology and metabolism in this plant. By comparative transcriptomic analysis, a series of candidate genes involved in the biosynthetic pathway of triterpenoid saponins are identified, providing new insight into their biosynthesis at transcriptome level.

Two CYP716A subfamily cytochrome P450 monooxygenases of sweet basil play similar but nonredundant roles in ursane- and oleanane-type pentacyclic triterpene biosynthesis

The medicinal plant sweet basil (Ocimum basilicum) accumulates bioactive ursane- and oleanane-type pentacyclic triterpenes (PCTs), ursolic acid and oleanolic acid, respectively, in a spatio-temporal manner; however, the biosynthetic enzymes and their contributions towards PCT biosynthesis remain to be elucidated. Two CYP716A subfamily cytochrome P450 monooxygenases (CYP716A252 and CYP716A253) are identified from a methyl jasmonate-responsive expression sequence tag collection and functionally characterized, employing yeast (Saccharomyces cerevisiae) expression platform and adapting virus-induced gene silencing (VIGS) in sweet basil. CYP716A252 and CYP716A253 catalyzed sequential three-step oxidation at the C-28 position of α-amyrin and β-amyrin to produce ursolic acid and oleanolic acid, respectively. Although CYP716A253 was more efficient than CYP716A252 for amyrin C-28 oxidation in yeast, VIGS revealed essential roles for both of these CYP716As in constitutive biosynthesis of ursolic acid and oleanolic acid in sweet basil leaves. However, CYP716A253 played a major role in elicitor-induced biosynthesis of ursolic acid and oleanolic acid. Overall, the results suggest similar as well as distinct roles of CYP716A252 and CYP716A253 for the spatio-temporal biosynthesis of PCTs. CYP716A252 and CYP716A253 might be useful for the alternative and sustainable production of PCTs in microbial host, besides increasing plant metabolite content through genetic modification.

CYP716A179 functions as a triterpene C-28 oxidase in tissue-cultured stolons of Glycyrrhiza uralensis

CYP716A179, a cytochrome P450 monooxygenase expressed predominantly in tissue-cultured stolons of licorice ( Glycyrrhiza uralensis ), functions as a triterpene C-28 oxidase in the biosynthesis of oleanolic acid and betulinic acid. Cytochrome P450 monooxygenases (P450s) play key roles in the structural diversification of plant triterpenoids. Among these, the CYP716A subfamily, which functions mainly as a triterpene C-28 oxidase, is common in plants. Licorice (Glycyrrhiza uralensis) produces bioactive triterpenoids, such as glycyrrhizin and soyasaponins, and relevant P450s (CYP88D6, CYP72A154, and CYP93E3) have been identified; however, no CYP716A subfamily P450 has been isolated. Here, we identify CYP716A179, which functions as a triterpene C-28 oxidase, by RNA sequencing analysis of tissue-cultured stolons of G. uralensis. Heterologous expression of CYP716A179 in engineered yeast strains confirmed the production of oleanolic acid, ursolic acid, and betulinic acid from β-amyrin, α-amyrin, and lupeol, respectively. The transcript level of CYP716A179 was about 500 times higher in tissue-cultured stolons than in intact roots. Oleanolic acid and betulinic acid were consistently detected only in tissue-cultured stolons. The discovery of CYP716A179 helps increase our understanding of the mechanisms of tissue-type-dependent triterpenoid metabolism in licorice and provides an additional target gene for pathway engineering to increase the production of glycyrrhizin in licorice tissue cultures by disrupting competing pathways.

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

Triterpenoids are widespread bioactive plant defence compounds with potential use as pharmaceuticals, pesticides and other high-value products. Enzymes belonging to the cytochrome P450 family have an essential role in creating the immense structural diversity of triterpenoids across the plant kingdom. However, for many triterpenoid oxidation reactions, the corresponding enzyme remains unknown. Here we characterize CYP716 enzymes from different medicinal plant species by heterologous expression in engineered yeasts and report ten hitherto unreported triterpenoid oxidation activities, including a cyclization reaction, leading to a triterpenoid lactone. Kingdom-wide phylogenetic analysis of over 400 CYP716s from over 200 plant species reveals details of their evolution and suggests that in eudicots the CYP716s evolved specifically towards triterpenoid biosynthesis. Our findings underscore the great potential of CYP716s as a source for generating triterpenoid structural diversity and expand the toolbox available for synthetic biology programmes for sustainable production of bioactive plant triterpenoids.

Cytochrome P450 family enzymes have an essential role in the creation of triterpenoid diversity in plants. Here, the authors describe triterpenoid synthesis as mediated by CYP716 enzymes in medicinal plant species, and perform phylogenetic analysis to describe CYP716 molecular evolution in plants.

Functional Characterization of CYP716 Family P450 Enzymes in Triterpenoid Biosynthesis in Tomato

Triterpenoids are a group of structurally diverse specialized metabolites that frequently show useful bioactivities. These chemicals are biosynthesized from the common precursor 2,3-oxidosqualene in plants. The carbon skeletons produced by oxidosqualene cyclase (OSC) are usually modified by cytochrome P450 monooxygenases (P450s) and UDP-dependent glycosyltransferases. These biosynthetic enzymes contribute to the structural diversification of plant triterpenoids. Until now, many P450 enzymes have been characterized as triterpenoid oxidases. Among them, the CYP716 family P450 enzymes, which have been isolated from a wide range of plant families, seem to contribute to the triterpenoid structural diversification. Many CYP716 family P450 enzymes have been characterized as the multifunctional triterpene C-28 oxidases, which oxidize α-amyrin and β-amyrin to the widely distributed triterpenoids ursolic and oleanolic acids, respectively. Tomato (Solanum lycopersicum) is one of the most important solanaceous crops in the world. However, little information is known regarding its triterpenoid biosynthesis. To understand the mechanism of triterpenoid biosynthesis in tomato, we focused on the function of CYP716 family enzymes as triterpenoid oxidases. We isolated all six CYP716 family genes from the Micro-Tom cultivar of tomato, and functionally characterized them in the heterologous yeast expression system. The in vivo enzymatic assays showed that CYP716A44 and CYP716A46 exhibited the ordinary C-28 oxidation activity against α-amyrin and β-amyrin to produce ursolic and oleanolic acids, respectively. Interestingly, one CYP716E subfamily enzyme, CYP716E26, exhibited the previously unreported C-6β hydroxylation activity against β-amyrin to produce a rare bioactive triterpenoid, daturadiol (olean-12-ene-3β,6β-diol). To determine the roles of the CYP716 family genes in tomato triterpenoid biosynthesis, we analyzed the gene expression and triterpenoid accumulation patterns in different plant tissues by performing the quantitative real-time polymerase chain reaction (qPCR) and gas chromatography-mass spectrometry (GC-MS) analyses, respectively. High levels of the CYP716A44 gene expression and the accumulation of C-28-oxidized triterpenoids, ursolic acid, and oleanolic acid were observed in the roots, indicating a significant contribution of the CYP716A44 gene in the triterpenoid biosynthesis in tomato. Thus, our study partially elucidated the mechanism of triterpenoid biosynthesis in tomato, and identified CYP716E26 as a novel C-6β hydroxylase for its subsequent use in the combinatorial biosynthesis of bioactive triterpenoids.

Triterpene Structural Diversification by Plant Cytochrome P450 Enzymes

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

Saponin determination, expression analysis and functional characterization of saponin biosynthetic genes in Chenopodium quinoa leaves

Quinoa (Chenopodium quinoa Willd.) is a highly nutritious pseudocereal with an outstanding protein, vitamin, mineral and nutraceutical content. The leaves, flowers and seed coat of quinoa contain triterpenoid saponins, which impart bitterness to the grain and make them unpalatable without postharvest removal of the saponins. In this study, we quantified saponin content in quinoa leaves from Ecuadorian sweet and bitter genotypes and assessed the expression of saponin biosynthetic genes in leaf samples elicited with methyl jasmonate. We found saponin accumulation in leaves after MeJA treatment in both ecotypes tested. As no reference genes were available to perform qPCR in quinoa, we mined publicly available RNA-Seq data for orthologs of 22 genes known to be stably expressed in Arabidopsis thaliana using geNorm, NormFinder and BestKeeper algorithms. The quinoa ortholog of At2g28390 (Monensin Sensitivity 1, MON1) was stably expressed and chosen as a suitable reference gene for qPCR analysis. Candidate saponin biosynthesis genes were screened in the quinoa RNA-Seq data and subsequent functional characterization in yeast led to the identification of CqbAS1, CqCYP716A78 and CqCYP716A79. These genes were found to be induced by MeJA, suggesting this phytohormone might also modulate saponin biosynthesis in quinoa leaves. Knowledge of the saponin biosynthesis and its regulation in quinoa may aid the further development of sweet cultivars that do not require postharvest processing.