Cytochrome P450 promiscuity leads to a bifurcating biosynthetic pathway for tanshinones

Research advances in cytochrome P450-catalysed pharmaceutical terpenoid biosynthesis in plants

Terpenoids, the biggest class of plant secondary metabolites, have a wide range of significant physiological roles, while many of them are important natural drugs. Biosynthesis of pharmaceutical terpenoids in plants is a fairly complex process, most of which involves cytochrome P450 (CYP450) monooxygenases. CYP450 enzymes are versatile biocatalysts that play critical roles in terpenoid skeleton modification and structural diversity. Therefore, the discovery and identification of CYP450 genes is significant for elucidating the terpenoid biosynthetic pathway. This review summarizes the progress and cloning strategies relating to CYP450s in pharmaceutical terpenoid biosynthesis of the past decade.

Efficient production of glycyrrhetinic acid in metabolically engineered Saccharomyces cerevisiae via an integrated strategy

Background

Glycyrrhetinic acid (GA) is the most important ingredient in licorice due to its outstanding anti-inflammatory activity and wide application in the medicine and cosmetics industries. Contemporary industrial production of GA by acid hydrolysis of glycyrrhizin which was extracted from Glycyrrhiza plants, is not environment-friendly and devastates farmland since the Glycyrrhiza rhizomes grow up to 10 m underground.

Results

In this study, GA was produced through metabolically engineering Saccharomyces cerevisiae by introducing the entire heterogeneous biosynthetic pathway of GA. Codon optimized CYP88D6 and CYP72A154, combined with β-AS (β-amyrin synthase encoding gene) and the NADPH-cytochrome P450 reductase gene of Arabidopsis thaliana were introduced into S. cerevisiae. The resulting strain (Y1) produced 2.5 mg/L of β-amyrin and 14 μg/L of GA. The cytochrome b5 from G. uralensis (GuCYB5) was identified and the introduction of this novel GuCYB5 increased the efficiency of GA production by eightfold. The joint utilization of the GuCYB5 gene along with 10 known MVA pathway genes from S. cerevisiae were overexpressed in a stable chromosome integration to achieve higher GA production. Using the combined strategy, GA concentration improved by 40-fold during batch fermentation. The production was further improved to 8.78 mg/L in fed-batch fermentation, which was increased by a factor of nearly 630.

Conclusions

This study first investigated the influence of carbon flux in the upstream module and the introduction of a newly identified GuCYB5 on GA production. The newly identified GuCYB5 was highly effective in improving GA production. An integrated strategy including enzyme discovery, pathway optimization, and fusion protein construction was provided in improving GA production, achieving a 630 fold increase in GA production. The metabolically engineered yeast cell factories provide an alternative approach to glycyrrhetinic acid production, replacing the traditional method of plant extraction.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1138-5) contains supplementary material, which is available to authorized users.

Proteomic analysis reveals novel insights into tanshinones biosynthesis in Salvia miltiorrhiza hairy roots

Salvia miltiorrhiza is a medicinal plant highly appreciated by its content of tanshinones and salvianolic acids. Tanshinones are of particular relevance for their anti-oxidant, anti-tumoral and anti-inflammatory properties. Abiotic and biotic agents as silver nitrate and yeast extract have shown efficiently to stimulate tanshinone accumulation, but the underlying molecular mechanism remains essentially unknown. By using hairy roots as experimental material and the elicitors mentioned, were obtained up to 22 mg of tanshinones per gram of dry weight. Differential label-free quantitative proteomic analysis was applied to study the proteins involved in tanshinone biosynthesis. A total of 2650 proteins were identified in roots extracts, of which 893 showed statistically (p < 0.05) significant change in relative abundance compared to control roots, 251 proteins were upregulated and 642 downregulated. Among the upregulated proteins the predominant functional categories were metabolism (47%), stress defense (18%) and redox homeostasis (10%). Within the metabolism category, isoprenoid metabolism enzymes, cytochromes P450 and FAD-binding berberine proteins showed abundance profile linked to tanshinone concentration. The results presented here allowed to propose 5 new cytochromes P450 and 5 berberine enzymes as candidates to be involved into tanshinone biosynthesis, a novel finding that opens new avenues to improve tanshinone production through biotechnological approaches.

Functional characterization of the cytochrome P450 monooxygenase CYP71AU87 indicates a role in marrubiin biosynthesis in the medicinal plant Marrubium vulgare

BACKGROUND

Horehound (Marrubium vulgare) is a medicinal plant whose signature bioactive compounds, marrubiin and related furanoid diterpenoid lactones, have potential applications for the treatment of cardiovascular diseases and type II diabetes. Lack of scalable plant cultivation and the complex metabolite profile of M. vulgare limit access to marrubiin via extraction from plant biomass. Knowledge of the marrubiin-biosynthetic enzymes can enable the development of metabolic engineering platforms for marrubiin production. We previously identified two diterpene synthases, MvCPS1 and MvELS, that act sequentially to form 9,13-epoxy-labd-14-ene. Conversion of 9,13-epoxy-labd-14-ene by cytochrome P450 monooxygenase (P450) enzymes can be hypothesized to facilitate key functional modification reactions in the formation of marrubiin and related compounds.

RESULTS

Mining a M. vulgare leaf transcriptome database identified 95 full-length P450 candidates. Cloning and functional analysis of select P450 candidates showing high transcript abundance revealed a member of the CYP71 family, CYP71AU87, that catalyzed the hydroxylation of 9,13-epoxy-labd-14-ene to yield two isomeric products, 9,13-epoxy labd-14-ene-18-ol and 9,13-epoxy labd-14-ene-19-ol, as verified by GC-MS and NMR analysis. Additional transient Nicotiana benthamiana co-expression assays of CYP71AU87 with different diterpene synthase pairs suggested that CYP71AU87 is specific to the sequential MvCPS1 and MvELS product 9,13-epoxy-labd-14-ene. Although the P450 products were not detectable in planta, high levels of CYP71AU87 gene expression in marrubiin-accumulating tissues supported a role in the formation of marrubiin and related diterpenoids in M. vulgare.

CONCLUSIONS

In a sequential reaction with the diterpene synthase pair MvCPS1 and MvELS, CYP71AU87 forms the isomeric products 9,13-epoxy labd-14-ene-18/19-ol as probable intermediates in marrubiin biosynthesis. Although its metabolic relevance in planta will necessitate further genetic studies, identification of the CYP71AU87 catalytic activity expands our knowledge of the functional landscape of plant P450 enzymes involved in specialized diterpenoid metabolism and can provide a resource for the formulation of marrubiin and related bioactive natural products.

Tanshinones, Critical Pharmacological Components in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge, a member of the Lamiaceae family, is valued in traditional Chinese Medicine. Its dried root (named Danshen) has been used for hundreds of years, primarily for the treatment of cardiovascular and cerebrovascular diseases. Tanshinones are the main active ingredients in S. miltiorrhiza and exhibit significant pharmacological activities, such as antioxidant activity, anti-inflammatory activity, cardiovascular effects, and antitumor activity. Danshen dripping pill of Tianshili is an effective drug widely used in the clinical treatment of cardiovascular diseases. With the increasing demand for clinical drugs, the traditional method for extracting and separating tanshinones from medicinal plants is insufficient. Therefore, in combination with synthetic biological methods and strategies, it is necessary to analyze the biosynthetic pathway of tanshinones and construct high-yield functional bacteria to obtain tanshinones. Moreover, the biosynthesis of tanshinones has been studied for more than two decades but remains to be completely elucidated. This review will systematically present the composition, extraction and separation, pharmacological activities and biosynthesis of tanshinones from S. miltiorrhiza, with the intent to provide references for studies on other terpenoid bioactive components of traditional Chinese medicines and to provide new research strategies for the sustainable development of traditional Chinese medicine resources.

Biochemical engineering in China

Chinese biochemical engineering is committed to supporting the chemical and food industries, to advance science and technology frontiers, and to meet major demands of Chinese society and national economic development. This paper reviews the development of biochemical engineering, strategic deployment of these technologies by the government, industrial demand, research progress, and breakthroughs in key technologies in China. Furthermore, the outlook for future developments in biochemical engineering in China is also discussed.

CYP76B74 Catalyzes the 3''-Hydroxylation of Geranylhydroquinone in Shikonin Biosynthesis

Shikonin and its derivatives are the most abundant naphthoquinone pigments formed in species of the medicinally and economically valuable Boraginaceae. A key step in the shikonin biosynthetic pathway, namely the C-3'' hydroxylation of the prenylated phenolic intermediate geranylhydroquinone, is expected to be catalyzed by a cytochrome P450. To identify cytochrome P450 candidates with transcription profiles similar to those of genes known to be involved in shikonin biosynthesis, we carried out coexpression analysis of transcriptome data sets of shikonin-proficient and shikonin-deficient cell lines and examined the spatial expression of candidate genes in different organs of Arnebia euchroma In biochemical assays using geranylhydroquinone as the substrate, CYP76B74 exhibited geranylhydroquinone 3''-hydroxylase activity and produced 3''-hydroxy-geranylhydroquinone. In CYP76B74 RNA interference A. euchroma hairy roots, shikonin derivative accumulation decreased dramatically, which demonstrated that CYP76B74 is required for shikonin biosynthesis in the plant. Phylogenetic analysis confirmed that CYP76B74 belonged to the CYP76B subfamily and was most likely derived from an ancestral geraniol 10-hydroxylase. In a subcellular localization analysis, a GFP-CYP76B74 fusion localized to endoplasmic reticulum membranes. Our results demonstrate that CYP76B74 catalyzes the key hydroxylation step in shikonin biosynthesis with high efficiency. The characterization of the CYP76B74 described here paves the way for further exploration of the ring closure reactions generating the naphthoquinone skeleton as well as for the alternative metabolism of geranylhydroquinone 3''-hydroxylase to dihydroechinofuran.

Enhancement of tanshinone production in Salvia miltiorrhiza hairy root cultures by metabolic engineering

BACKGROUND

Tanshinones are diterpenoid compounds that are used to treat cardiovascular diseases. As current extraction methods for tanshinones are inefficient, there is a pressing need to improve the production of these bioactive compounds to meet increasing demand.

RESULTS

Overexpression of SmMDS (2-c-methyl-d-erythritol 2,4-cyclodiphosphate synthase, a tanshinone biosynthesis gene) in transgenic Salvia miltiorrhiza hairy roots significantly increased the tanshinone yield compared to the control, and total tanshinone content in SmMDS-overexpressing lines increased after elicitor treatment. Total tanshinones increased to 2.5, 2.3, and 3.2 mg/g DW (dry weight) following treatment with Ag⁺, YE (yeast extract), and MJ (methyl jasmonate), respectively, compared with the non-induced transgenic line (1.7 mg/g DW). Also, qRT-PCR analysis showed that the expression levels of two pathway genes was positively correlated with increased accumulation of tanshinone.

CONCLUSIONS

Our study provides an effective strategy for increasing the content of tanshinones and other natural compounds using a combination of genetic engineering and elicitor treatment.

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.

Engineering of CYP76AH15 can improve activity and specificity towards forskolin biosynthesis in yeast

Background

Forskolin is a high-value diterpenoid produced exclusively by the Lamiaceae plant Coleus forskohlii. Today forskolin is used pharmaceutically for its adenyl-cyclase activating properties. The limited availability of pure  forskolin is currently hindering its full utilization, thus a new environmentally friendly, scalable and sustainable strategy is needed for forskolin production. Recently, the entire biosynthetic pathway leading to forskolin was elucidated. The key steps of the pathway are catalyzed by cytochrome P450 enzymes (CYPs), which have been shown to be the limiting steps of the pathway. Here we study whether protein engineering of the substrate recognition sites (SRSs) of CYPs can improve their efficiency towards forskolin biosynthesis in yeast.

Results

As a proof of concept, we engineered the enzyme responsible for the first putative oxygenation step of the forskolin pathway: the conversion of 13R-manoyl oxide to 11-oxo-13R-manoyl oxide, catalyzed by the CYP76AH15. Four CYP76AH15 variants—engineered in the SRS regions—yielded at least a twofold increase of 11-oxo-13R-manoyl oxide when expressed in yeast cells grown in microtiter plates. The highest titers (5.6-fold increase) were observed with the variant A99I, mutated in the SRS1 region. Double or triple CYP76AH15 mutant variants resulted in additional enzymes with optimized performances. Moreover, in planta CYP76AH15 can synthesize ferruginol from miltiradiene. In this work, we showed that the mutants affecting 11-oxo-13R-manoyl oxide synthesis, do not affect ferruginol production, and vice versa.

The best performing variant, A99I, was utilized to reconstruct the forskolin biosynthetic pathway in yeast cells. Although these strains showed increased 11-oxo-manoyl oxide production and higher accumulation of other pathway intermediates compared to the native CYP76AH15, lower production of forskolin was observed.

Conclusions

As demonstrated for CYP76AH15, site-directed mutagenesis of SRS regions of plant CYPs may be an efficient and targeted approach to increase the performance of these enzymes. Although in this work we have managed to achieve higher efficiency and specificity of the first CYP of the pathway, further work is necessary in order to increase the overall production of forskolin in yeast cells.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1027-3) contains supplementary material, which is available to authorized users.

The AP2/ERF transcription factor SmERF1L1 regulates the biosynthesis of tanshinones and phenolic acids in Salvia miltiorrhiza

Tanshinones and phenolic acids are two important metabolites synthesized by the traditional Chinese medicinal plant Salvia miltiorrhiza. There is increasing market demand for these compounds. Here, we isolated and functionally characterized SmERF1L1, a novel JA (Jasmonic acid)-responsive gene encoding AP2/ERF transcription factor, from Salvia miltiorrhiza. SmERF1L1 was responsive to methyl jasmonate (MJ), yeast extraction (YE), salicylic acid (SA) and ethylene treatments. Subcellular localization assay indicated that SmERF1L1 located in the nucleus. Overexpression of SmERF1L1 significantly increased tanshinones production in transgenic S. miltiorrhiza hairy roots by comprehensively upregulating tanshinone biosynthetic pathway genes, especially SmDXR. Yeast one-hybrid (Y1H) and electrophoretic mobility shift assay (EMSA) showed that SmERF1L1 binds to the GCC-box of SmDXR promoter while dual luciferase (Dual-LUC) assay showed that SmERF1L1 positively regulated the expression of SmDXR. Our study suggested that the SmERF1L1 may be a good potential target for further metabolic engineering of bioactive component biosynthesis in S. miltiorrhiza.

Trends in herbgenomics

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

Engineering yeast metabolism for production of terpenoids for use as perfume ingredients, pharmaceuticals and biofuels

Terpenoids represent a large class of natural products with significant commercial applications. These chemicals are currently mainly obtained through extraction from plants and microbes or through chemical synthesis. However, these sources often face challenges of unsustainability and low productivity. In order to address these issues, Escherichia coli and yeast have been metabolic engineered to produce non-native terpenoids. With recent reports of engineering yeast metabolism to produce several terpenoids at high yields, it has become possible to establish commercial yeast production of terpenoids that find applications as perfume ingredients, pharmaceuticals and advanced biofuels. In this review, we describe the strategies to rewire the yeast pathway for terpenoid biosynthesis. Recent advances will be discussed together with challenges and perspectives of yeast as a cell factory to produce different terpenoids.

Transcriptional Profiles of SmWRKY Family Genes and Their Putative Roles in the Biosynthesis of Tanshinone and Phenolic Acids in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a Chinese traditional herb for treating cardiovascular and cerebrovascular diseases, and tanshinones and phenolic acids are the dominated medicinal and secondary metabolism constituents of this plant. WRKY transcription factors (TFs) can function as regulators of secondary metabolites biosynthesis in many plants. However, studies on the WRKY that regulate tanshinones and phenolics biosynthesis are limited. In this study, 69 SmWRKYs were identified in the transcriptome database of S. miltiorrhiza, and phylogenetic analysis indicated that some SmWRKYs had closer genetic relationships with other plant WRKYs, which were involved in secondary metabolism. Hairy roots of S. miltiorrhiza were treated by methyl jasmonate (MeJA) to detect the dynamic change trend of SmWRKY, biosynthetic genes, and medicinal ingredients accumulation. Base on those date, a correlation analysis using Pearson’s correlation coefficient was performed to construct gene-to-metabolite network and identify 9 SmWRKYs (SmWRKY1, 7, 19, 29, 45, 52, 56, 58, and 68), which were most likely to be involved in tanshinones and phenolic acids biosynthesis. Taken together, this study has provided a significant resource that could be used for further research on SmWRKY in S. miltiorrhiza and especially could be used as a cue for further investigating SmWRKY functions in secondary metabolite accumulation.

The Identification of Core Gene Expression Signature in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignancies, which causes serious financial burden worldwide. This study aims to investigate the potential mechanisms contributing to HCC and identify core biomarkers. The HCC gene expression profile GSE41804 was picked out to analyze the differentially expressed genes (DEGs). Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out using DAVID. We constructed a protein-protein interaction (PPI) network to visualize interactions of the DEGs. The survival analysis of these hub genes was conducted to evaluate their potential effects on HCC. In this analysis, 503 DEGs were captured (360 downregulated genes and 143 upregulated genes). Meanwhile, 15 hub genes were identified. GO analysis showed that the DEGs were mainly enriched in oxidative stress, cell cycle, and extracellular structure. KEGG analysis suggested the DEGs were enriched in the absorption, metabolism, and cell cycle pathway. PPI network disclosed that the top3 modules were mainly enriched in cell cycle, oxidative stress, and liver detoxification. In conclusion, our analysis uncovered that the alterations of oxidative stress and cell cycle are two major signatures of HCC. TOP2A, CCNB1, and KIF4A might promote the development of HCC, especially in proliferation and differentiation, which could be novel biomarkers and targets for diagnosis and treatment of HCC.

SmJAZ8 acts as a core repressor regulating JA-induced biosynthesis of salvianolic acids and tanshinones in Salvia miltiorrhiza hairy roots

Jasmonates (JAs) are important plant hormones that regulate a variety of plant development and defense processes, including biosynthesis of secondary metabolites. The JASMONATE ZIM DOMAIN (JAZ) proteins act as negative regulators in the JA signaling pathways of plants. We first verified that methyl jasmonate (MeJA) enhanced the accumulation of both salvianolic acids and tanshinones in Salvia miltiorrhiza (Danshen) hairy roots by inducing the expression of their biosynthetic pathway genes. Nine JAZ genes were cloned from Danshen and their expression levels in hairy roots were all increased by treatment with MeJA. When analyzed in detail, however, SmJAZ8 showed the strongest expression in the induced hairy roots. Overexpression or RNAi of SmJAZ8 deregulated or up-regulated the yields of salvianolic acids and tanshinones in the MeJA-induced transgenic hairy roots, respectively, and transcription factors and biosynthetic pathway genes showed an expression pattern that mirrored the production of the compounds. Genetic transformation of SmJAZ8 altered the expression of other SmJAZ genes, suggesting evidence of crosstalk occurring in JAZ-regulated secondary metabolism. Furthermore, the transcriptome analysis revealed a primary-secondary metabolism balance regulated by SmJAZ8. Altogether, we propose a novel role for SmJAZ8 as a negative feedback loop controller in the JA-induced biosynthesis of salvianolic acids and tanshinones.

Biosynthesis of bioactive diterpenoids in the medicinal plant Vitex agnus-castus

Vitex agnus-castus L. (Lamiaceae) is a medicinal plant historically used throughout the Mediterranean region to treat menstrual cycle disorders, and is still used today as a clinically effective treatment for premenstrual syndrome. The pharmaceutical activity of the plant extract is linked to its ability to lower prolactin levels. This feature has been attributed to the presence of dopaminergic diterpenoids that can bind to dopamine receptors in the pituitary gland. Phytochemical analyses of V. agnus-castus show that it contains an enormous array of structurally related diterpenoids and, as such, holds potential as a rich source of new dopaminergic drugs. The present work investigated the localisation and biosynthesis of diterpenoids in V. agnus-castus. With the assistance of matrix-assisted laser desorption ionisation-mass spectrometry imaging (MALDI-MSI), diterpenoids were localised to trichomes on the surface of fruit and leaves. Analysis of a trichome-specific transcriptome database, coupled with expression studies, identified seven candidate genes involved in diterpenoid biosynthesis: three class II diterpene synthases (diTPSs); three class I diTPSs; and a cytochrome P450 (CYP). Combinatorial assays of the diTPSs resulted in the formation of a range of different diterpenes that can account for several of the backbones of bioactive diterpenoids observed in V. agnus-castus. The identified CYP, VacCYP76BK1, was found to catalyse 16-hydroxylation of the diol-diterpene, peregrinol, to labd-13Z-ene-9,15,16-triol when expressed in Saccharomyces cerevisiae. Notably, this product is a potential intermediate in the biosynthetic pathway towards bioactive furan- and lactone-containing diterpenoids that are present in this species.

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.

Diverse responses of tanshinone biosynthesis to biotic and abiotic elicitors in hairy root cultures of Salvia miltiorrhiza and Salvia castanea Diels f. tomentosa

Salvia miltiorrhiza (S. miltiorrhiza) and Salvia castanea Diels f. tomentosa (S. castanea) are both used for treatment of cardiovascular diseases. They have the same bioactive compound tanshinones, but whose contents are hugely different. This study illustrated diverse responses of tanshinone biosynthesis to yeast extract (YE) and Ag⁺ in hairy roots of the two species. YE enhanced both the growth and tanshinone biosynthesis of two hairy roots, and contributed more to tanshinone accumulation in S. castanea than that in S. miltiorrhiza. Genes encoding 1-deoxy-d-xylulose 5-phosphate synthase (DXS2), geranylgeranyl diphosphatesynthase (GGPPS1), copalyl diphosphate synthase (CPS1), and two cytochromes P450 (CYP76AH1 and CYP76AH3) were also more responsive to YE in S. castanea than those in S. miltiorrhiza. Accumulations of dihydrotanshinone I and tanshinone I, and most biosynthetic genes in S. miltiorrhiza were more responsive to Ag⁺ than those in S. castanea. Accumulations of dihydrotanshinone I and cryptotanshinone were more responsive to YE, while tanshinone IIA accumulation was more responsive to Ag⁺ in S. miltiorrhiza. However, accumulations of other four tanshinones and related genes in S. castanea were more responsive to YE than Ag⁺. This study provides foundations for studying diverse specialized metabolism between the related species.

Comprehensive transcriptome profiling of Salvia miltiorrhiza for discovery of genes associated with the biosynthesis of tanshinones and phenolic acids

Tanshinones and phenolic acids are crucial bioactive compounds biosynthesized in Salvia miltiorrhiza. Methyl jasmonate (MeJA) is an effective elicitor to enhance the production of phenolic acids and tanshinones simultaneously, while yeast extract (YE) is used as a biotic elicitor that only induce tanshinones accumulation. However, little was known about the different molecular mechanism. To identify the downstream and regulatory genes involved in tanshinone and phenolic acid biosynthesis, we conducted comprehensive transcriptome profiling of S. miltiorrhiza hairy roots treated with either MeJA or YE. Total 55588 unigenes were assembled from about 1.72 billion clean reads, of which 42458 unigenes (76.4%) were successfully annotated. The expression patterns of 19 selected genes in the significantly upregulated unigenes were verified by quantitative real-time PCR. The candidate downstream genes and other cytochrome P450s involved in the late steps of tanshinone and phenolic acid biosynthesis pathways were screened from the RNA-seq dataset based on co-expression pattern analysis with specific biosynthetic genes. Additionally, 375 transcription factors were identified to exhibit a significant up-regulated expression pattern in response to induction. This study can provide us a valuable gene resource for elucidating the molecular mechanism of tanshinones and phenolic acids biosynthesis in hairy roots of S. miltiorrhiza.

Overcoming the plasticity of plant specialized metabolism for selective diterpene production in yeast

Plants synthesize numerous specialized metabolites (also termed natural products) to mediate dynamic interactions with their surroundings. The complexity of plant specialized metabolism is the result of an inherent biosynthetic plasticity rooted in the substrate and product promiscuity of the enzymes involved. The pathway of carnosic acid-related diterpenes in rosemary and sage involves promiscuous cytochrome P450s whose combined activity results in a multitude of structurally related compounds. Some of these minor products, such as pisiferic acid and salviol, have established bioactivity, but their limited availability prevents further evaluation. Reconstructing carnosic acid biosynthesis in yeast achieved significant titers of the main compound but could not specifically yield the minor products. Specific production of pisiferic acid and salviol was achieved by restricting the promiscuity of a key enzyme, CYP76AH24, through a single-residue substitution (F112L). Coupled with additional metabolic engineering interventions, overall improvements of 24 and 14-fold for pisiferic acid and salviol, respectively, were obtained. These results provide an example of how synthetic biology can help navigating the complex landscape of plant natural product biosynthesis to achieve heterologous production of useful minor metabolites. In the context of plant adaptation, these findings also suggest a molecular basis for the rapid evolution of terpene biosynthetic pathways.

Engineering microbial cell factories for the production of plant natural products: from design principles to industrial-scale production

Plant natural products (PNPs) are widely used as pharmaceuticals, nutraceuticals, seasonings, pigments, etc., with a huge commercial value on the global market. However, most of these PNPs are still being extracted from plants. A resource-conserving and environment-friendly synthesis route for PNPs that utilizes microbial cell factories has attracted increasing attention since the 1940s. However, at the present only a handful of PNPs are being produced by microbial cell factories at an industrial scale, and there are still many challenges in their large-scale application. One of the challenges is that most biosynthetic pathways of PNPs are still unknown, which largely limits the number of candidate PNPs for heterologous microbial production. Another challenge is that the metabolic fluxes toward the target products in microbial hosts are often hindered by poor precursor supply, low catalytic activity of enzymes and obstructed product transport. Consequently, despite intensive studies on the metabolic engineering of microbial hosts, the fermentation costs of most heterologously produced PNPs are still too high for industrial-scale production. In this paper, we review several aspects of PNP production in microbial cell factories, including important design principles and recent progress in pathway mining and metabolic engineering. In addition, implemented cases of industrial-scale production of PNPs in microbial cell factories are also highlighted.

SmMYB36, a Novel R2R3-MYB Transcription Factor, Enhances Tanshinone Accumulation and Decreases Phenolic Acid Content in Salvia miltiorrhiza Hairy Roots

Phenolic acids and tanshinones are two major bioactive components in Salvia miltiorrhiza Bunge. A novel endogenous R2R3-MYB transcription factor, SmMYB36, was identified in this research. This transcript factor can simultaneously influence the content of two types of components in SmMYB36 overexpression hairy roots. SmMYB36 was mainly localized in the nucleus of onion epidermis and it has transactivation activity. The overexpression of SmMYB36 promoted tanshinone accumulation but inhibited phenolic acid and flavonoid biosynthesis in Salvia miltiorrhiza hairy roots. The altered metabolite content was due to changed metabolic flow which was regulated by transcript expression of metabolic pathway genes. The gene transcription levels of the phenylpropanoid general pathway, tyrosine derived pathway, methylerythritol phosphate pathway and downstream tanshinone biosynthetic pathway changed significantly due to the overexpression of SmMYB36. The wide distribution of MYB binding elements (MBS, MRE, MBSI and MBSII) and electrophoretic mobility shift assay results indicated that SmMYB36 may be an effective tool to regulate metabolic flux shifts.

A Global Coexpression Network Approach for Connecting Genes to Specialized Metabolic Pathways in Plants

Plants produce diverse specialized metabolites (SMs), but the genes responsible for their production and regulation remain largely unknown, hindering efforts to tap plant pharmacopeia. Given that genes comprising SM pathways exhibit environmentally dependent coregulation, we hypothesized that genes within a SM pathway would form tight associations (modules) with each other in coexpression networks, facilitating their identification. To evaluate this hypothesis, we used 10 global coexpression data sets, each a meta-analysis of hundreds to thousands of experiments, across eight plant species to identify hundreds of coexpressed gene modules per data set. In support of our hypothesis, 15.3 to 52.6% of modules contained two or more known SM biosynthetic genes, and module genes were enriched in SM functions. Moreover, modules recovered many experimentally validated SM pathways, including all six known to form biosynthetic gene clusters (BGCs). In contrast, bioinformatically predicted BGCs (i.e., those lacking an associated metabolite) were no more coexpressed than the null distribution for neighboring genes. These results suggest that most predicted plant BGCs are not genuine SM pathways and argue that BGCs are not a hallmark of plant specialized metabolism. We submit that global gene coexpression is a rich, largely untapped resource for discovering the genetic basis and architecture of plant natural products.

Biosynthesis of the oxygenated diterpene nezukol in the medicinal plant Isodon rubescens is catalyzed by a pair of diterpene synthases

Plants produce an immense diversity of natural products (i.e. secondary or specialized metabolites) that offer a rich source of known and potentially new pharmaceuticals and other desirable bioproducts. The Traditional Chinese Medicinal plant Isodon rubescens (Lamiaceae) contains an array of bioactive labdane-related diterpenoid natural products. Of these, the ent-kauranoid oridonin is the most prominent specialized metabolite that has been extensively studied for its potent antimicrobial and anticancer efficacy. Mining of a previously established transcriptome of I. rubescens leaf tissue identified seven diterpene synthase (diTPSs) candidates. Here we report the functional characterization of four I. rubescens diTPSs. IrTPS5 and IrTPS3 were identified as an ent-copalyl diphosphate (CPP) synthase and a (+)-CPP synthase, respectively. Distinct transcript abundance of IrTPS5 and the predicted ent-CPP synthase IrTPS1 suggested a role of IrTPS5 in specialized ent-kaurene metabolism possibly en route to oridonin. Nicotiana benthamiana co-expression assays demonstrated that IrTPS4 functions sequentially with IrTPS3 to form miltiradiene. In addition, IrTPS2 converted the IrTPS3 product (+)-CPP into the hydroxylated tricyclic diterpene nezukol not previously identified in I. rubescens. Metabolite profiling verified the presence of nezukol in I. rubescens leaf tissue. The proposed IrTPS2-catalyzed reaction mechanism proceeds via the common ionization of the diphosphate group of (+)-CPP, followed by formation of an intermediary pimar-15-en-8-yl+ carbocation and neutralization of the carbocation by water capture at C-8 to yield nezukol, as confirmed by nuclear magnetic resonance (NMR) analysis. Oxygenation activity is rare for the family of class I diTPSs and offers new catalysts for developing metabolic engineering platforms to produce a broader spectrum of bioactive diterpenoid natural products.

Identification of a Novel (-)-5-Epieremophilene Synthase from Salvia miltiorrhiza via Transcriptome Mining

Salvia miltiorrhiza, a medicinal plant in China, has been used for thousands of years to treat coronary heart diseases. Although biosynthesis of tanshinones, a group of diterpenoids in S. miltiorrhiza, has been extensively investigated, to date we know little about the formation of monoterpenes and sesquiterpenes in this medicinal plant. Here, we report the characterization of three sesquiterpene synthases, named SmSTPS1, SmSTPS2, and SmSTPS3, which catalyzed the formation of a new compound, (-)-5-epieremophilene. Additionally, the (-)-5-epieremophilene biosynthesis activity of SmSTPS1 was confirmed by transient expression in Nicotiana benthamiana. Despite the similar enzyme activities of SmSTPS1, SmSTPS2, and SmSTPS3, the three (-)-5-epieremophilene synthase genes displayed different spatial expression patterns and responded differently to hormone treatments, implicating their specific roles in plant-environment interactions. Our results provide valuable data to understanding the biosynthesis and composition of terpenes in plant.

The 2-oxoglutarate-dependent dioxygenase superfamily participates in tanshinone production in Salvia miltiorrhiza

Highly oxidized tanshinones are pharmacological ingredients extracted from the medicinal model plant Salvia miltiorrhiza and are mainly used to treat cardiovascular diseases. Previous studies have confirmed that cytochrome P450 mono-oxygenases (CYP450s) have a key function in the biosynthesis of tanshinones; however, no solid evidence links oxidation to the 2-oxoglutarate-dependent dioxygenase (2OGD) superfamily. Here, we identified 132 members of the DOXB and DOXC subfamilies of 2OGD by scanning the 2OG-FeII Oxy domain using a genome-wide strategy in S. miltiorrhiza. The DOXC class was phylogenetically divided into twelve clades. Combining phylogenetic relationships, differential expression and co-expression from various organs and tissues revealed that two 2OGDs were directly related to flavonoid metabolism, and that 13 2OGDs from different clades were predicted to be involved in tanshinone biosynthesis. Based on this insight into tanshinone production, we experimentally detected significant decreases in miltirone, cryptotanshinone, and tanshinone IIA (0.16-, 0.56-, and 0.56-fold, respectively) in 2OGD5 RNAi transgenic lines relative to the control lines using a metabonomics analysis. 2OGD5 was found to play a crucial role in the downstream biosynthesis of tanshinones following the hydroxylation of CYPs. Our results highlight the evolution and diversification of 2OGD superfamily members and suggest that they contribute to the complexity of tanshinone metabolites.

Total biosynthesis of the cyclic AMP booster forskolin from Coleus forskohlii

Forskolin is a unique structurally complex labdane-type diterpenoid used in the treatment of glaucoma and heart failure based on its activity as a cyclic AMP booster. Commercial production of forskolin relies exclusively on extraction from its only known natural source, the plant Coleus forskohlii, in which forskolin accumulates in the root cork. Here, we report the discovery of five cytochrome P450s and two acetyltransferases which catalyze a cascade of reactions converting the forskolin precursor 13R-manoyl oxide into forskolin and a diverse array of additional labdane-type diterpenoids. A minimal set of three P450s in combination with a single acetyl transferase was identified that catalyzes the conversion of 13R-manoyl oxide into forskolin as demonstrated by transient expression in Nicotiana benthamiana. The entire pathway for forskolin production from glucose encompassing expression of nine genes was stably integrated into Saccharomyces cerevisiae and afforded forskolin titers of 40 mg/L.

Elucidation of terpenoid metabolism in Scoparia dulcis by RNA-seq analysis

Scoparia dulcis biosynthesize bioactive diterpenes, such as scopadulcic acid B (SDB), which are known for their unique molecular skeleton. Although the biosynthesis of bioactive diterpenes is catalyzed by a sequence of class II and class I diterpene synthases (diTPSs), the mechanisms underlying this process are yet to be fully identified. To elucidate these biosynthetic machinery, we performed a high-throughput RNA-seq analysis, and de novo assembly of clean reads revealed 46,332 unique transcripts and 40,503 two unigenes. We found diTPSs genes including a putative syn-copalyl diphosphate synthase (SdCPS2) and two kaurene synthase-like (SdKSLs) genes. Besides them, total 79 full-length of cytochrome P450 (CYP450) genes were also discovered. The expression analyses showed selected CYP450s associated with their expression pattern of SdCPS2 and SdKSL1, suggesting that CYP450 candidates involved diterpene modification. SdCPS2 represents the first predicted gene to produce syn-copalyl diphosphate in dicots. In addition, SdKSL1 potentially contributes to the SDB biosynthetic pathway. Therefore, these identified genes associated with diterpene biosynthesis lead to the development of genetic engineering focus on diterpene metabolism in S. dulcis.

Targeted mutagenesis in the medicinal plant Salvia miltiorrhiza

CRISPR/Cas9 is a powerful genome editing tool that has been extensively used in model plants and crops, such as Arabidopsis thaliana, rice, wheat, and soybean. Here, we report the use of CRISPR/Cas9 to precisely knock out the committed diterpene synthase gene (SmCPS1) involved in tanshinone biosynthesis in Salvia miltiorrhiza, a traditional Chinese medicinal herb with significant pharmacological activities, such as vasorelaxation, protection against ischemia-reperfusion injury, and antiarrhythmic effects. Three homozygous and eight chimeric mutants were obtained from 26 independent transgenic hairy root lines by Agrobacterium rhizogenes-mediated transformation. The metabolomic analysis based on LC-qTOF-MS and Q-TRAP-LC-MS/MS revealed that tanshinones, especially cryptotanshinone, tanshinone IIA and tanshinone I, are completely missing in homozygous mutants, without influencing other phenolic acid metabolites. By contrast, tanshinones are decreased but still detectable in chimeric mutants, which is similar to a previously-reported an RNAi study of SmCPS1. These results demonstrate that Agrobacterium rhizogenes- mediated transformation using CRISPR/Cas9 is a simple and efficient genome editing tool in S. miltiorrhiza, thus paving the way for large-scale genome editing in S. miltiorrhiza, which is important for pathway elucidation of secondary metabolites, quality improvement, and yield increases for this valuable traditional Chinese medicinal herb.

Elucidation of the biosynthesis of carnosic acid and its reconstitution in yeast

Rosemary extracts containing the phenolic diterpenes carnosic acid and its derivative carnosol are approved food additives used in an increasingly wide range of products to enhance shelf-life, thanks to their high anti-oxidant activity. We describe here the elucidation of the complete biosynthetic pathway of carnosic acid and its reconstitution in yeast cells. Cytochrome P450 oxygenases (CYP76AH22-24) from Rosmarinus officinalis and Salvia fruticosa already characterized as ferruginol synthases are also able to produce 11-hydroxyferruginol. Modelling-based mutagenesis of three amino acids in the related ferruginol synthase (CYP76AH1) from S. miltiorrhiza is sufficient to convert it to a 11-hydroxyferruginol synthase (HFS). The three sequential C20 oxidations for the conversion of 11-hydroxyferruginol to carnosic acid are catalysed by the related CYP76AK6-8. The availability of the genes for the biosynthesis of carnosic acid opens opportunities for the metabolic engineering of phenolic diterpenes, a class of compounds with potent anti-oxidant, anti-inflammatory and anti-tumour activities.

Diterpenes are plant products with high antioxidant properties and potential application as food additives and therapeutics. Here, the authors describe the complete biosynthetic pathway of carnosic acid and reconstruct it in yeast, opening the way to metabolic engineering of phenolic diterpenes.

A novel glucuronosyltransferase has an unprecedented ability to catalyse continuous two-step glucuronosylation of glycyrrhetinic acid to yield glycyrrhizin

Glycyrrhizin is an important bioactive compound that is used clinically to treat chronic hepatitis and is also used as a sweetener world-wide. However, the key UDP-dependent glucuronosyltransferases (UGATs) involved in the biosynthesis of glycyrrhizin remain unknown. To discover unknown UGATs, we fully annotated potential UGATs from Glycyrrhiza uralensis using deep transcriptome sequencing. The catalytic functions of candidate UGATs were determined by an in vitro enzyme assay. Systematically screening 434 potential UGATs, we unexpectedly found one unique GuUGAT that was able to catalyse the glucuronosylation of glycyrrhetinic acid to directly yield glycyrrhizin via continuous two-step glucuronosylation. Expression analysis further confirmed the key role of GuUGAT in the biosynthesis of glycyrrhizin. Site-directed mutagenesis revealed that Gln-352 may be important for the initial step of glucuronosylation, and His-22, Trp-370, Glu-375 and Gln-392 may be important residues for the second step of glucuronosylation. Notably, the ability of GuUGAT to catalyse a continuous two-step glucuronosylation reaction was determined to be unprecedented among known glycosyltransferases of bioactive plant natural products. Our findings increase the understanding of traditional glycosyltransferases and pave the way for the complete biosynthesis of glycyrrhizin.

Orthogonal Assays Clarify the Oxidative Biochemistry of Taxol P450 CYP725A4

Natural product metabolic engineering potentially offers sustainable and affordable access to numerous valuable molecules. However, challenges in characterizing and assembling complex biosynthetic pathways have prevented more rapid progress in this field. The anticancer agent Taxol represents an excellent case study. Assembly of a biosynthetic pathway for Taxol has long been stalled at its first functionalization, putatively an oxygenation performed by the cytochrome P450 CYP725A4, due to confounding characterizations. Here, through combined in vivo (Escherichia coli), in vitro (lipid nanodisc), and metabolite stability assays, we verify the presence and likely cause of this enzyme's inherent promiscuity. Thereby, we remove the possibility that promiscuity simply existed as an artifact of previous metabolic engineering approaches. Further, spontaneous rearrangement and the stabilizing effect of a hydrophobic overlay suggest a potential role for nonenzymatic chemistry in Taxol's biosynthesis. Taken together, this work confirms taxadiene-5α-ol as a primary enzymatic product of CYP725A4 and provides direction for future Taxol metabolic and protein engineering efforts.

Extreme promiscuity of a bacterial and a plant diterpene synthase enables combinatorial biosynthesis

Diterpenes are widely distributed across many biological kingdoms, where they serve a diverse range of physiological functions, and some have significant industrial utility. Their biosynthesis involves class I diterpene synthases (DTSs), whose activity can be preceded by that of class II diterpene cyclases (DTCs). Here, a modular metabolic engineering system was used to examine the promiscuity of DTSs. Strikingly, both a bacterial and plant DTS were found to exhibit extreme promiscuity, reacting with all available precursors with orthogonal activity, producing an olefin or hydroxyl group, respectively. Such DTS promiscuity enables combinatorial biosynthesis, with remarkably high yields for these unoptimized non-native enzymatic combinations (up to 15mg/L). Indeed, it was possible to readily characterize the 13 unknown products. Notably, 16 of the observed diterpenes were previously inaccessible, and these results provide biosynthetic routes that are further expected to enable assembly of more extended pathways to produce additionally elaborated 'non-natural' diterpenoids.