Cytochrome P450 promiscuity leads to a bifurcating biosynthetic pathway for tanshinones

Integrated multi-omics analysis reveals genes involved in flavonoid biosynthesis and trichome development of Artemisia argyi

Artemisia argyi is an herbaceous plant of the genus Artemisia. Its young and mature leaves are used as food and medicine, respectively. Glandular trichomes (GTs) are distributed on the leaf surface in A. argyi and are generally considered the location of flavonoid biosynthesis and accumulation. However, the mechanism of flavonoid biosynthesis and accumulation in A. argyi remains unclear. In this study, the coregulatory genes involved in flavonoid biosynthesis and trichome development in this species were screened and evaluated, and the biosynthetic pathways for key flavonoids in A. argyi were uncovered. AaMYB1 and AaYABBY1 were screened using weighted gene co-expression network analysis, and both genes were then genetically transformed into Nicotiana tabacum L. cv. K326 (tobacco). Simultaneously, AaYABBY1 was also genetically transformed into Arabidopsis thaliana. The total flavonoid and rutin contents were increased in tobacco plants overexpressing AaMYB1 and AaYABBY1, and the expression levels of genes participating in the flavonoid synthesis pathway, such as PAL, FLS, and F3H, were significantly up-regulated in plants overexpressing these genes. These results indicated that AaMYB1 and AaYABBY1 promote flavonoid biosynthesis in tobacco. Furthermore, compared to that in the wild-type, the trichome density was significantly increased in tobacco and A. thaliana plants overexpressing AaYABBY1. These results confirm that AaYABBY1 might be involved in regulating trichome formation in A. argyi. This indicates the potential genes involved in and provides new insights into the development of trichome cellular factories based on the "development-metabolism" interaction network and the cultivation of high-quality A. argyi.

Genome-Wide Identification and Functional Analysis of Salvia miltiorrhiza MicroRNAs Reveal the Negative Regulatory Role of Smi-miR159a in Phenolic Acid Biosynthesis

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs in plants. They play critical functions in various biological processes during plant growth and development. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant with significant medicinal, economic, and academic values. In order to elucidate the role of miRNAs in S. miltiorrhiza, six small RNA libraries from mature roots, young roots, stems, mature leaves, young leaves and flowers of S. miltiorrhiza and one degradome library from mixed tissues were constructed. A total of 184 miRNA precursors, generating 137 known and 49 novel miRNAs, were genome-widely identified. The identified miRNAs were predicted to play diversified regulatory roles in plants through regulating 891 genes. qRT-PCR and 5' RLM-RACE assays validated the negative regulatory role of smi-miR159a in SmMYB62, SmMYB78, and SmMYB80. To elucidate the function of smi-miR159a in bioactive compound biosynthesis, smi-miR159a transgenic hairy roots were generated and analyzed. The results showed that overexpression of smi-miR159a caused a significant decrease in rosmarinic acid and salvianolic acid B contents. qRT-PCR analysis showed that the targets of smi-miR159a, including SmMYB62, SmMYB78, and SmMYB80, were significantly down-regulated, accompanied by the down-regulation of SmPAL1, SmC4H1, Sm4CL1, SmTAT1, SmTAT3, SmHPPR1, SmRAS, and SmCYP98A14 genes involved in phenolic acid biosynthesis. It suggests that smi-miR159a is a significant negative regulator of phenolic acid biosynthesis in S. miltiorrhiza.

Molecular characterization and structure basis of a malonyltransferase with both substrate promiscuity and catalytic regiospecificity from Cistanche tubulosa

Enzymatic malonylation of natural glycosides provides a promising alternative method for drug-like malonylated glycosides supply. However, the catalytic potential and structural basis of plant malonyltransferase are far from being fully elucidated. This work identified a new malonyltransferase CtMaT1 from Cistanche tubulosa. It displayed unprecedented mono- and/or di-malonylation activity toward diverse glucosides with different aglycons. A "one-pot" system by CtMaT1 and a malonyl-CoA synthetase was established to biosynthesize nine new malonylated glucosides. Structural investigations revealed that CtMaT1 possesses an adequately spacious acyl-acceptor pocket capable of accommodating diverse glucosides. Additionally, it recognizes malonyl-CoA through strong electrotactic and hydrogen interactions. QM/MM calculation revealed the H167-mediated SN2 reaction mechanism of CtMaT1, while dynamic simulations detected the formation of stable hydrogen bonds between the glucose-6-OH group and H167, resulting in its high malonylation regiospecificity. Calculated energy profiles of two isomeric glycosides highlighted lower reaction energy barriers towards glucoside substrates, emphasizing CtMaT1's preference for glucosides. Furthermore, a mutant CtMaT1H36A with notably increased di-malonylation activity was obtained. The underlying molecular mechanism was illuminated through MM/GBSA binding free energy calculation. This study significantly advances the understanding of plant acyltransferases from both functional and protein structural perspectives, while also providing a versatile tool for enzymatic malonylation applications in pharmacology.

Analysis of ginseng rusty root symptoms transcriptome and its pathogenesis directed by reactive oxygen species theory

Ginseng rusty root symptoms (GRS) is a primary disease of ginseng, which seriously decreases the yield and quality of ginseng and causes enormous losses to ginseng production. GRS prevention and control is still challenging due to its unclear etiology. In this study, the phloem tissue of healthy Panax ginseng (AG), the nonred tissue of the phloem epidermis around the lesion (BG), and the red lesion site tissue of GRS (CG) were extracted for mRNA transcriptomic analysis; 35,958 differentially expressed genes (DEGs) were identified and were associated with multiple stress resistance pathways, reactive oxygen species (ROS), and iron ion binding. Further study showed that the contents of O2 •-, H2O2, and malondialdehyde (MDA) were significantly increased in BG and CG tissues. Under anaerobic conditions caused by excessive soil moisture, the overproduction of ROS destroys cell membranes, simultaneously converting Fe2+ to Fe3+ and depositing it in the cell wall, which results in GRS, as evidenced by the success of the GRS induction test.

The Smi-miR858a-SmMYB module regulates tanshinone and phenolic acid biosynthesis in Salvia miltiorrhiza

Tanshinones and phenolic acids are two major classes of bioactive compounds in Salvia miltiorrhiza. Revealing the regulatory mechanism of their biosynthesis is crucial for quality improvement of S. miltiorrhiza medicinal materials. Here we demonstrated that Smi-miR858a-Smi-miR858c, a miRNA family previously known to regulate flavonoid biosynthesis, also played critical regulatory roles in tanshinone and phenolic acid biosynthesis in S. miltiorrhiza. Overexpression of Smi-miR858a in S. miltiorrhiza plants caused significant growth retardation and tanshinone and phenolic acid reduction. Computational prediction and degradome and RNA-seq analyses revealed that Smi-miR858a could directly cleave the transcripts of SmMYB6, SmMYB97, SmMYB111, and SmMYB112. Yeast one-hybrid and transient transcriptional activity assays showed that Smi-miR858a-regulated SmMYBs, such as SmMYB6 and SmMYB112, could activate the expression of SmPAL1 and SmTAT1 involved in phenolic acid biosynthesis and SmCPS1 and SmKSL1 associated with tanshinone biosynthesis. In addition to directly activating the genes involved in bioactive compound biosynthesis pathways, SmMYB6, SmMYB97, and SmMYB112 could also activate SmAOC2, SmAOS4, and SmJMT2 involved in the biosynthesis of methyl jasmonate, a significant elicitor of plant secondary metabolism. The results suggest the existence of dual signaling pathways for the regulation of Smi-miR858a in bioactive compound biosynthesis in S. miltiorrhiza.

DNA methylation regulates biosynthesis of tanshinones and phenolic acids during growth of Salvia miltiorrhiza

DNA methylation plays a crucial role in the regulation of plant growth and the biosynthesis of secondary metabolites. Danshen (Salvia miltiorrhiza) is a valuable Chinese herbal medicine commonly used to treat cardiovascular diseases; its active ingredients are tanshinones and phenolic acids, which primarily accumulate in roots. Here, we conducted a targeted metabolic analysis of S. miltiorrhiza roots at 3 distinct growth stages: 40 d old (r40), 60 d old (r60), and 90 d old (r90). The contents of tanshinones (cryptotanshinone, tanshinone I, tanshinone IIA, and rosmariquinone) and phenolic acids (rosmarinic acid and salvianolic acid B) gradually increased during plant development. Whole-genome bisulfite sequencing and transcriptome sequencing of roots at the 3 growth stages revealed an increased level of DNA methylation in the CHH context (H represents A, T, or C) context at r90 compared with r40 and r60. Increased DNA methylation levels were associated with elevated expression of various genes linked to epigenetic regulations, including CHROMOMETHYLASE2 (SmCMT2), Decrease in DNA Methylation 1 (SmDDM1), Argonaute 4 (SmAGO4), and DOMAINS REARRANGED METHYLTRANSFERASE 1 (SmDRM1). Moreover, expression levels of many genes involved in tanshinone and salvianolic acid biosynthesis, such as copalyldiphosphate synthase 5 (SmCPS5), cytochrome P450-related enzyme (SmCYP71D464), geranylgeranyl diphosphate synthase (SmGGPPS1), geranyl diphosphate synthase (SmGPPS), hydroxyphenylpyruvate reductase (SmHPPR), and hydroxyphenylpyruvate dioxygenase (SmHPPD), were altered owing to hyper-methylation, indicating that DNA methylation plays an important role in regulating tanshinone and phenolic acid accumulation. Our data shed light on the epigenetic regulation of root growth and the biosynthesis of active ingredients in S. miltiorrhiza, providing crucial clues for further improvement of active compound production via molecular breeding in S. miltiorrhiza.

Endophytic fungus Penicillium steckii DF33 promoted tanshinones biosynthesis in Salvia miltiorrhiza by regulating the expression of CYP450 genes

Salvia miltiorrhiza, a prominent traditional Chinese medicinal resource, has been extensively employed in the management of cardiovascular and cerebrovascular ailments. Ensuring the consistency of S. miltiorrhiza raw materials revolves around the imperative task of maintaining stable tanshinones content and composition. An effective approach in this regard involves the utilization of endophytic fungi as inducers. Within this context, our study spotlights an endophytic fungus, Penicillium steckii DF33, isolated from the roots of S. miltiorrhiza. Remarkably, this fungus has demonstrated a significant capacity to boost the biosynthesis and accumulation of tanshinones. The primary objective of this investigation is to elucidate the underlying regulatory mechanism by which DF33 enhances and regulates the biosynthesis and accumulation of tanshinones. This is achieved through its influence on the differential expression of crucial CYP450 genes within the S. miltiorrhiza hairy roots system. The results revealed that the DF33 elicitor not only promotes the growth of hairy roots but also enhances the accumulation of tanshinones. Notably, the content of cryptotanshinone was reached 1.6452 ± 0.0925 mg g-1, a fourfold increase compared to the control group. Our qRT-PCR results further demonstrate that the DF33 elicitor significantly up-regulates the expression of most key enzyme genes (GGPPS, CPS1, KSL1, CYP76AH1, CYP76AH3, CYP76AK1, CYP71D411) involved in the tanshinone biosynthesis pathway. This effect is particularly pronounced in certain critical CYP450 genes and Tanshinone ⅡA synthase (SmTⅡAS), with their expression levels peaking at 7 days or 14 days, respectively. In summary, endophytic P. steckii DF33 primarily enhances tanshinone biosynthesis by elevating the expression levels of pivotal enzyme genes associated with the modification and transformation stages within the tanshinone biosynthesis pathway. These findings underscore the potential of employing plant probiotics, specifically endophytic and root-associated microbes, to facilitate the biosynthesis and transformation of vital constituents in medicinal plants, and this approach holds promise for enhancing the quality of traditional Chinese medicinal materials.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Mining and functional characterization of NADPH-cytochrome P450 reductases of the DNJ biosynthetic pathway in mulberry leaves

BACKGROUND

1-Deoxynojirimycin (DNJ), the main active ingredient in mulberry leaves, with wide applications in the medicine and food industries due to its significant functions in lowering blood sugar, and lipids, and combating viral infections. Cytochrome P450 is a key enzyme for DNJ biosynthesis, its activity depends on the electron supply of NADPH-cytochrome P450 reductases (CPRs). However, the gene for MaCPRs in mulberry leaves remains unknown.

RESULTS

In this study, we successfully cloned and functionally characterized two key genes, MaCPR1 and MaCPR2, based on the transcriptional profile of mulberry leaves. The MaCPR1 gene comprised 2064 bp, with its open reading frame (ORF) encoding 687 amino acids. The MaCPR2 gene comprised 2148 bp, and its ORF encoding 715 amino acids. The phylogenetic tree indicates that MaCPR1 and MaCPR2 belong to Class I and Class II, respectively. In vitro, we found that the recombinant enzymes MaCPR2 protein could reduce cytochrome c and ferricyanide using NADPH as an electron donor, while MaCPR1 did not. In yeast, heterologous co-expression indicates that MaCPR2 delivers electrons to MaC3'H hydroxylase, a key enzyme catalyzing the production of chlorogenic acid from 3-O-p-coumaroylquinic acid.

CONCLUSIONS

These findings highlight the orchestration of hydroxylation process mediated by MaCPR2 during the biosynthesis of secondary metabolite biosynthesis in mulberry leaves. These results provided a foundational understanding for fully elucidating the DNJ biosynthetic pathway within mulberry leaves.

Comparative transcriptomics of two Salvia subg. Perovskia species contribute towards molecular background of abietane-type diterpenoid biosynthesis

Tanshinones, are a group of diterpenoid red pigments present in Danshen - an important herbal drug of Traditional Chinese Medicine which is a dried root of Salvia miltiorrhiza Bunge. Some of the tanshinones are sought after as pharmacologically active natural products. To date, the biosynthetic pathway of tanshinones has been only partially elucidated. These compounds are also present in some of the other Salvia species, i.a. from subgenus Perovskia, such as S. abrotanoides (Kar.) Sytsma and S. yangii B.T. Drew. Despite of the close genetic relationship between these species, significant qualitative differences in their diterpenoid profile have been discovered. In this work, we have used the Liquid Chromatography-Mass Spectrometry analysis to follow the content of diterpenoids during the vegetation season, which confirmed our previous observations of a diverse diterpenoid profile. As metabolic differences are reflected in different transcript profile of a species or tissues, we used metabolomics-guided transcriptomic approach to select candidate genes, which expression possibly led to observed chemical differences. Using an RNA-sequencing technology we have sequenced and de novo assembled transcriptomes of leaves and roots of S. abrotanoides and S. yangii. As a result, 134,443 transcripts were annotated by UniProt and 56,693 of them were assigned as Viridiplantae. In order to seek for differences, the differential expression analysis was performed, which revealed that 463, 362, 922 and 835 genes indicated changes in expression in four comparisons. GO enrichment analysis and KEGG functional analysis of selected DEGs were performed. The homology and expression of two gene families, associated with downstream steps of tanshinone and carnosic acid biosynthesis were studied, namely: cytochromes P-450 and 2-oxoglutarate-dependend dioxygenases. Additionally, BLAST analysis revealed existence of 39 different transcripts related to abietane diterpenoid biosynthesis in transcriptomes of S. abrotanoides and S. yangii. We have used quantitative real-time RT-PCR analysis of selected candidate genes, to follow their expression levels over the vegetative season. A hypothesis of an existence of a multifunctional CYP76AH89 in transcriptomes of S. abrotanoides and S. yangii is discussed and potential roles of other CYP450 homologs are speculated. By using the comparative transcriptomic approach, we have generated a dataset of candidate genes which provides a valuable resource for further elucidation of tanshinone biosynthesis. In a long run, our investigation may lead to optimization of diterpenoid profile in S. abrotanoides and S. yangii, which may become an alternative source of tanshinones for further research on their bioactivity and pharmacological therapy.

Oxetane Ring Formation in Taxol Biosynthesis Is Catalyzed by a Bifunctional Cytochrome P450 Enzyme

Taxol is a potent drug used in various cancer treatments. Its complex structure has prompted extensive research into its biosynthesis. However, certain critical steps, such as the formation of the oxetane ring, which is essential for its activity, have remained unclear. Previous proposals suggested that oxetane formation follows the acetylation of taxadien-5α-ol. Here, we proposed that the oxetane ring is formed by cytochrome P450-mediated oxidation events that occur prior to C5 acetylation. To test this hypothesis, we analyzed the genomic and transcriptomic information for Taxus species to identify cytochrome P450 candidates and employed two independent systems, yeast (Saccharomyces cerevisiae) and plant (Nicotiana benthamiana), for their characterization. We revealed that a single enzyme, CYP725A4, catalyzes two successive epoxidation events, leading to the formation of the oxetane ring. We further showed that both taxa-4(5)-11(12)-diene (endotaxadiene) and taxa-4(20)-11(12)-diene (exotaxadiene) are precursors to the key intermediate, taxologenic oxetane, indicating the potential existence of multiple routes in the Taxol pathway. Thus, we unveiled a long-elusive step in Taxol biosynthesis.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Engineering yeast for the production of plant terpenoids using synthetic biology approaches

Covering: 2011-2022The low amounts of terpenoids produced in plants and the difficulty in synthesizing these complex structures have stimulated the production of terpenoid compounds in microbial hosts by metabolic engineering and synthetic biology approaches. Advances in engineering yeast for terpenoid production will be covered in this review focusing on four directions: (1) manipulation of host metabolism, (2) rewiring and reconstructing metabolic pathways, (3) engineering the catalytic activity, substrate selectivity and product specificity of biosynthetic enzymes, and (4) localizing terpenoid production via enzymatic fusions and scaffolds, or subcellular compartmentalization.

Miltiradiene Production by Cytoplasmic Metabolic Engineering in Nicotiana benthamiana

Plant natural products are important sources of innovative drugs, but the extraction and isolation of medicinal natural products from plants is challenging as these compounds have complex structures that are difficult to synthesize chemically. Therefore, utilizing heterologous expression systems to produce medicinal natural products in plants is a novel, environmentally friendly, and sustainable method. In this study, Nicotiana benthamiana was used as the plant platform to successfully produce miltiradiene, the key intermediate of tanshinones, which are the bioactive constituents of the Chinese medicinal plant Salvia miltiorrhiza. The yield of miltiradiene was increased through cytoplasmic engineering strategies combined with the enhancement of isoprenoid precursors. Additionally, we discovered that overexpressing SmHMGR alone accelerated apoptosis in tobacco leaves. Due to the richer membrane systems and cofactors in tobacco compared to yeast, tobacco is more conducive to the expression of plant enzymes. Therefore, this study lays the foundation for dissecting the tanshinone biosynthetic pathway in tobacco, which is essential for subsequent research. Additionally, it highlights the potential of N. benthamiana as an alternative platform for the production of natural products in plants.

Functional characterization of CYP81C16 involved in the tanshinone biosynthetic pathway in Salvia miltiorrhiza

Danshen, the dried roots and rhizomes of Salvia miltiorrhiza Bunge (S. miltiorrhiza), is widely used in the treatment of cardiovascular and cerebrovascular diseases. Tanshinones, the bioactive compounds from Danshen, exhibit a wide spectrum of pharmacological properties, suggesting their potential for future therapeutic applications. Tanshinone biosynthesis is a complex process involving at least six P450 enzymes that have been identified and characterized, most of which belong to the CYP76 and CYP71 families. In this study, CYP81C16, a member of the CYP71 clan, was identified in S. miltiorrhiza. An in vitro assay revealed that it could catalyze the hydroxylation of four para-quinone-type tanshinones, namely neocryptotanshinone, deoxyneocryptotanshinone, and danshenxinkuns A and B. SmCYP81C16 emerged as a potential broad-spectrum oxidase targeting the C-18 position of para-quinone-type tanshinones with an impressive relative conversion rate exceeding 90%. Kinetic evaluations andin vivo assays underscored its highest affinity towards neocryptotanshinone among the tested substrates. The overexpression of SmCYP81C16 promoted the accumulation of (iso)tanshinone in hairy root lines. The characterization of SmCYP81C16 in this study accentuates its potential as a pivotal tool in the biotechnological production of tanshinones, either through microbial or plant metabolic engineering.

(±)-Salvicatone A: A Pair of C27-Meroterpenoid Enantiomers with Skeletons from the Roots and Rhizomes of Salvia castanea Diels f. tomentosa Stib

(±)-Salvicatone A (1), a C27-meroterpenoid featuring a unique 6/6/6/6/6-pentacyclic carbon skeleton with a 7,8,8a,9,10,10a-hexahydropyren-1 (6H)-one motif, was isolated from the roots and rhizomes of Salvia castanea Diels f. tomentosa Stib. Its structure was characterized by comprehensive spectroscopic analyses along with computer-assisted structure elucidation, including ACD/structure elucidator and quantum chemical calculations with 1H/13C NMR and electronic circular dichroism. Biogenetically, compound 1 was constructed from decarboxylation following [4 + 2] Diels-Alder cycloaddition reaction between caffeic acid and miltirone analogue. Bioassays showed that (-)-1 and (+)-1 inhibited nitric oxide production in lipopolysaccharide-induced RAW264.7 macrophage cells with an IC50 value of 6.48 ± 1.25 and 15.76 ± 5.55 μM, respectively. The structure-based virtual screening based on the pharmacophores in ePharmaLib, as well as the molecular docking and molecular dynamics simulations study, implied that (-)-1 and (+)-1 may potentially bind to retinoic acid receptor-related orphan receptor C to exert anti-inflammatory activities.

CYP72D19 from Tripterygium wilfordii catalyzes C-2 hydroxylation of abietane-type diterpenoids

KEY MESSAGE

CYP72D19, the first functional gene of the CYP72D subfamily, catalyzes the C-2 hydroxylation of abietane-type diterpenoids. The abietane-type diterpenoids, e.g., triptolide, tripdiolide, and 2-epitripdiolide, are the main natural products for the anti-tumor, anti-inflammatory, and immunosuppressive activities of Tripterygium wilfordii, while their biosynthetic pathways are not resolved. Here, we cloned and characterized the CYP72D19-catalyzed C-2 hydroxylation of dehydroabietic acid, a compound that has been proven to be a biosynthetic intermediate in triptolide biosynthesis. Through molecular docking and site-directed mutagenesis, L386, L387, and I493 near the active pocket were found to have an important effect on the enzyme activity, which also indicates that steric hindrance of residues plays an important role in function. In addition, CYP72D19 also catalyzed a variety of abietane-type diterpenoids with benzene ring, presumably because the benzene ring of the substrate molecule stabilized the C-ring, allowing the protein and the substrate to form a relatively stable spatial structure. This is the first demonstration of CYP72D subfamily gene function. Our research provides important genetic elements for the structural modification of active ingredients and the heterologous production of other 2-hydroxyl abietane-type natural products.

Identification of novel flavin-dependent monooxygenase from Strobilanthes Cusia reveals molecular basis of indoles' biosynthetic logic

BACKGROUND

Strobilanthes cusia (Nees) Kuntze is a traditional medical plant distributed widely in south China. The indole compounds that originated from the plant are responsible for its pharmacological activities. However, the reason why indole ingredients are accumulated in this herb and how it is biosynthesized has remained largely unknown.

RESULTS

In this study, metabolic and transcriptional profiling measurement experiments of different S. cusia organs were carried out to understand the underlying molecular basis of indoles' biosynthetic logic. A metabolic investigation demonstrated that the indoles are primarily accumulated mainly in aerial parts, particularly in leaves. RNA-seq was employed to reveal the organ specific accumulation of indoles in different S. cusia organs. Meanwhile, a flavin-dependent monooxygenase gene (ScFMO1) was found in S. cusia, and it has capacity to produce indoxyl from indole by the fermentation assay. Finally, we assessed the outcomes of transient expression experiment in tobacco and confirmed that ScFMO1 localizes in cytoplasm.

CONCLUSIONS

Our results suggest that ScFMO1 plays a key role in biosynthesis of indoles (Indigo, indirubin, indican, etc.), it will be useful for illuminating the molecular basis of the medicinal indoles' biosynthesis and developing strategies for improving their yields.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

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.

Unveiling the spatial distribution and molecular mechanisms of terpenoid biosynthesis in Salvia miltiorrhiza and S. grandifolia using multi-omics and DESI-MSI

Salvia miltiorrhiza and S. grandifolia are rich in diterpenoids and have therapeutic effects on cardiovascular diseases. In this study, the spatial distribution of diterpenoids in both species was analyzed by a combination of metabolomics and mass spectrometry imaging techniques. The results indicated that diterpenoids in S. miltiorrhiza were mainly abietane-type norditerpenoid quinones with a furan or dihydrofuran D-ring and were mainly distributed in the periderm of the roots, e.g. cryptotanshinone and tanshinone IIA. The compounds in S. grandifolia were mainly phenolic abietane-type tricyclic diterpenoids with six- or seven-membered C-rings, and were widely distributed in the periderm, phloem, and xylem of the roots, e.g. 11-hydroxy-sugiol, 11,20-dihydroxy-sugiol, and 11,20-dihydroxy-ferruginol. In addition, the leaves of S. grandifolia were rich in tanshinone biosynthesis precursors, such as 11-hydroxy-sugiol, while those of S. miltiorrhiza were rich in phenolic acids. Genes in the upstream pathway of tanshinone biosynthesis were highly expressed in the root of S. grandifolia, and genes in the downstream pathway were highly expressed in the root of S. miltiorrhiza. Here, we describe the specific tissue distributions and mechanisms of diterpenoids in two Salvia species, which will facilitate further investigations of the biosynthesis of diterpenoids in plant synthetic biology.

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

Analysis of the HD-Zip I transcription factor family in Salvia miltiorrhiza and functional research of SmHD-Zip12 in tanshinone synthesis

Background

The homeodomain-leucine zipper I (HD-Zip I) transcription factor is a plant-specific protein that plays an essential role in the abiotic stress response of plants. Research on the HD-Zip I family in Salvia miltiorrhiza is still lacking.

Methods and Results

In this study, a total of 25 SmHD-Zip I proteins were identified. Their characterizations, phylogenetic relationships, conserved motifs, gene structures, and cis-elements were analyzed comprehensively using bioinformatics methods. Expression profiling revealed that SmHD-Zip I genes exhibited distinctive tissue-specific patterns and divergent responses to ABA, PEG, and NaCl stresses. SmHD-Zip12 responded the most strongly to ABA, PEG, and NaCl, so it was used for transgenic experiments. The overexpression of SmHD-Zip12 significantly increased the content of cryptotanshinone, dihydrotanshinone I, tanshinone I, and tanshinone IIA by 2.89-fold, 1.85-fold, 2.14-fold, and 8.91-fold compared to the wild type, respectively. Moreover, in the tanshinone biosynthetic pathways, the overexpression of SmHD-Zip12 up-regulated the expression levels of SmAACT, SmDXS, SmIDS, SmGGPPS, SmCPS1, SmCPS2, SmCYP76AH1, SmCYP76AH3, and SmCYP76AK1 compared with the wild type.

Conclusions

This study provides information the possible functions of the HD-Zip I family and lays a theoretical foundation for clarifying the functional mechanism of the SmHD-Zip12 gene in regulating the synthesis of tanshinone in S. miltiorrhiza.

Caterpillars Detoxify Diterpenoid from Nepeta stewartiana by the Molting Hormone Gene CYP306A1

Herbivorous insects are well known for detoxifying a broad range of the defense compounds produced by the plants that they feed on, but knowledge of the mechanisms of detoxification is still very limited. Here, we describe a system in which two species of lepidopteran caterpillars metabolize an abietane diterpene from the plants of Nepeta stewartiana Diels to an oxygenated derivative that is less active biologically. We found that this transformation could be catalyzed by a cytochrome P450 enzyme in caterpillars, which are associated with molting. Most interestingly, abietane diterpene targets the molting-associated gene CYP306A1 to alter the content of molting hormones in the insect at specific developmental stages and competitively inhibit molting hormone metabolism. These findings identify the mechanism by which caterpillars are able to detoxify abietane diterpenoid through hydroxylation at the C-19 position, which may be opening up exciting research questions into the mechanisms of interaction between plants and insects.

Functional Study and Efficient Catalytic Element Mining of CYP76AHs in Salvia Plants

Salvia is a large genus with hundreds of species used in traditional Chinese medicine. Tanshinones are a highly representative class of exclusive compounds found in the Salvia genus that exhibit significant biological activity. Tanshinone components have been identified in 16 Salvia species. The CYP76AH subfamily (P450) is crucial for the synthesis of tanshinone due to its catalytic generation of polyhydroxy structures. In this study, a total of 420 CYP76AH genes were obtained, and phylogenetic analysis showed their clear clustering relationships. Fifteen CYP76AH genes from 10 Salvia species were cloned and studied from the perspectives of evolution and catalytic efficiency. Three CYP76AHs with significantly improved catalytic efficiency compared to SmCYP76AH3 were identified, providing efficient catalytic elements for the synthetic biological production of tanshinones. A structure-function relationship study revealed several conserved residues that might be related to the function of CYP76AHs and provided a new mutation direction for the study of the directed evolution of plant P450.

Chromosome-level genome assembly of Salvia miltiorrhiza with orange roots uncovers the role of Sm2OGD3 in catalyzing 15,16-dehydrogenation of tanshinones

Salvia miltiorrhiza is well known for its clinical practice in treating heart and cardiovascular diseases. Its roots, used for traditional Chinese medicine materials, are usually brick-red due to accumulation of red pigments, such as tanshinone IIA and tanshinone I. Here we report a S. miltiorrhiza line (shh) with orange roots. Compared with the red roots of normal S. miltiorrhiza plants, the contents of tanshinones with a single bond at C-15,16 were increased, whereas those with a double bond at C-15,16 were significantly decreased in shh. We assembled a high-quality chromosome-level genome of shh. Phylogenomic analysis showed that the relationship between two S. miltiorrhiza lines with red roots was closer than the relationship with shh. It indicates that shh could not be the mutant of an extant S. miltiorrhiza line with red roots. Comparative genomic and transcriptomic analyses showed that a 1.0 kb DNA fragment was deleted in shh Sm2OGD3m. Complementation assay showed that overexpression of intact Sm2OGD3 in shh hairy roots recovered furan D-ring tanshinone accumulation. Consistently, in vitro protein assay showed that Sm2OGD3 catalyzed the conversion of cyptotanshinone, 15,16-dihydrotanshinone I and 1,2,15,16-tetrahydrotanshinone I into tanshinone IIA, tanshinone I and 1,2-dihydrotanshinone I, respectively. Thus, Sm2OGD3 functions as tanshinone 15,16-dehydrogenase and is a key enzyme in tanshinone biosynthesis. The results provide novel insights into the metabolic network of medicinally important tanshinone compounds.

Plant terpene specialized metabolism: complex networks or simple linear pathways?

From the perspectives of pathway evolution, discovery and engineering of plant specialized metabolism, the nature of the biosynthetic routes represents a critical aspect. Classical models depict biosynthesis typically from an end-point angle and as linear, for example, connecting central and specialized metabolism. As the number of functionally elucidated routes increased, the enzymatic foundation of complex plant chemistries became increasingly well understood. The perception of linear pathway models has been severely challenged. With a focus on plant terpenoid specialized metabolism, we review here illustrative examples supporting that plants have evolved complex networks driving chemical diversification. The completion of several diterpene, sesquiterpene and monoterpene routes shows complex formation of scaffolds and their subsequent functionalization. These networks show that branch points, including multiple sub-routes, mean that metabolic grids are the rule rather than the exception. This concept presents significant implications for biotechnological production.

Identification and functional characterization of three cytochrome P450 genes for the abietane diterpenoid biosynthesis in Isodon lophanthoides

We identify two ferruginol synthases and a 11-hydroxyferruginol synthase from a traditional Chinese medicinal herb Isodon lophanthoides and propose their involvement in two independent abietane diterpenoids biosynthetic pathways. Isodon lophanthoides is a traditional Chinese medicinal herb rich in highly oxidized abietane-type diterpenoids. These compounds exhibit a wide range of pharmaceutical activities, yet the biosynthesis is barely known. Here, we describe the screening and functional characterization of P450s that oxidize the abietane skeleton abietatriene. We mainly focused on CYP76 family and identified 12 CYP76AHs by mining the RNA-seq data of I. lophanthoides. Among the 12 CYP76AHs, 6 exhibited similar transcriptional expression features as upstream diterpene synthases, including root or leaf-preferential expression pattern and highly MeJA inducibility. These six P450s were considered as first-tier candidates and functionally characterized in yeast and plant cells. In yeast assays showed that both CYP76AH42 and CYP76AH43 were ferruginol synthases hydroxylating the C12 position of abietatriene, whereas CYP76AH46 was characterized as a 11-hydroxyferruginol synthase which catalyzes two successive oxidations at C12 and C11 of abietatriene. Heterologous expression of three CYP76AHs in Nicotiana benthamiana resulted in the formation of ferruginol. qPCR analysis showed CYP76AH42 and CYP76AH43 were mainly expressed in the root, which was consistent with the distribution of ferruginol in the root periderms. CYP76AH46 was primarily expressed in the leaves where barely ferruginol or 11-hydroxyferruginol was detected. In addition to distinct organ-specific expression pattern, three CYP76AHs exhibited different genomic structures (w or w/o introns), low protein sequence identities (51-63%) and were placed in separate subclades in the phylogenetic tree. These results suggest that the identified CYP76AHs may be involved in at least two independent abietane biosynthetic pathways in the aerial and underground parts of I. lophanthoides.

A chromosome-level genome assembly reveals that tandem-duplicated CYP706V oxidase genes control oridonin biosynthesis in the shoot apex of Isodon rubescens

The ent-kaurenoids (e.g., oridonin and enmein) from the Isodon genus (Lamiaceae) are one class of diterpenoids with rich structural diversity and intriguing pharmaceutical activity. In contrast to the well-established gibberellin pathway, oxidative modifications diversifying the ent-kaurene skeleton in Isodon have remained undetermined for half a century. Here we report a chromosome-level genome assembly of I. rubescens, a well-recognized oridonin producer long favored by Asian people as a traditional herb with antitumor effects. The shoot apex was confirmed to be the actual region actively producing ent-kaurene diterpenoids. Through comparative genomics and phylogenetic analyses, we discovered a cluster of tandem-duplicated CYP706V oxygenase-encoding genes located on an ancient genomic block widely distributed in eudicots, whereas almost exclusively emerged in Isodon plants. In the shoot apex, IrCYP706V2 and IrCYP706V7 oxidized the ent-kaurene core in the initial stage of oridonin biosynthesis. Loss of CYP706Vs in other Lamiaceae plants offered an explanation for the specific kaurenoid production in Isodon plants. Moreover, we found that the Isodon genomes encode multiple diterpenoid synthases that are potentially involved in generating diterpenoid diversity. These findings provided new insights into the evolution of the lineage-specific diterpenoid pathway and laid a foundation for improving production of bioactive ent-kaurene-type diterpenoids by molecular breeding and synthetic biology approaches.

Plant (di)terpenoid evolution: from pigments to hormones and beyond

Covering: up to 2014-2022.Diterpenoid biosynthesis in plants builds on the necessary production of (E,E,E)-geranylgeranyl diphosphate (GGPP) for photosynthetic pigment production, with diterpenoid biosynthesis arising very early in land plant evolution, enabling stockpiling of the extensive arsenal of (di)terpenoid natural products currently observed in this kingdom. This review will build upon that previously published in the Annual Review of Plant Biology, with a stronger focus on enzyme structure-function relationships, as well as additional insights into the evolution of (di)terpenoid metabolism since generated.

Promiscuity, a Driver of Plant Cytochrome P450 Evolution?

Plant cytochrome P450 monooxygenases were long considered to be highly substrate-specific, regioselective and stereoselective enzymes, in this respect differing from their animal counterparts. The functional data that have recently accumulated clearly counter this initial dogma. Highly promiscuous P450 enzymes have now been reported, mainly in terpenoid pathways with functions in plant adaptation, but also some very versatile xenobiotic/herbicide metabolizers. An overlap and predictable interference between endogenous and herbicide metabolism are starting to emerge. Both substrate preference and permissiveness vary between plant P450 families, with high promiscuity seemingly favoring retention of gene duplicates and evolutionary blooms. Yet significant promiscuity can also be observed in the families under high negative selection and with essential functions, usually enhanced after gene duplication. The strategies so far implemented, to systematically explore P450 catalytic capacity, are described and discussed.

Cytochrome P450s in plant terpenoid biosynthesis: discovery, characterization and metabolic engineering

As the largest family of natural products, terpenoids play valuable roles in medicine, agriculture, cosmetics and food. However, the traditional methods that rely on direct extraction from the original plants not only produce low yields, but also result in waste of resources, and are not applicable at all to endangered species. Modern heterologous biosynthesis is considered a promising, efficient, and sustainable production method, but it relies on the premise of a complete analysis of the biosynthetic pathway of terpenoids, especially the functionalization processes involving downstream cytochrome P450s. In this review, we systematically introduce the biotech approaches used to discover and characterize plant terpenoid-related P450s in recent years. In addition, we propose corresponding metabolic engineering approaches to increase the effective expression of P450 and improve the yield of terpenoids, and also elaborate on metabolic engineering strategies and examples of heterologous biosynthesis of terpenoids in Saccharomyces cerevisiae and plant hosts. Finally, we provide perspectives for the biotech approaches to be developed for future research on terpenoid-related P450.

Leveraging yeast to characterize plant biosynthetic gene clusters

Plant biosynthetic gene clusters (BGCs) contain multiple physically clustered non-homologous genes that encode enzymes catalyzing diverse reactions in one plant natural product biosynthetic pathway. A growing number of plant BGCs have emerged as an underlying resource for understanding plant specialized metabolism and evolution, but the characterization remains challenging. Recent studies have demonstrated that baker's yeast can serve as a versatile platform for the characterization of plant BGCs, from single-gene characterization to multiple genes and hitherto unknown putative BGC validation and elucidation. In this review, we will summarize the strategies and examples of the applications of yeast in plant BGC characterization and share our perspective on the development of a systematic pipeline to fully leverage yeast to advance the understanding of plant BGCs and plant natural product biomanufacturing.

A tracking work on how Sm4CL2 re-directed the biosynthesis of salvianolic acids and tanshinones in Salvia miltiorrhiza hairy roots

Overexpression and antisense expression of Sm4CL2 re-directed the biosynthesis of salvianolic acids and tanshinones in Salvia miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a widely used traditional Chinese medicine and its main active ingredients are water-soluble phenolic acids and lipophilic diterpenoids which are produced through the phenylpropanoid pathway and terpenoid pathway, respectively. 4-Coumaric acid: Coenzyme A ligase (4CL) is a key enzyme in the phenylpropanoid metabolism. We had obtained Sm4CL2-overexpressing (Sm4CL2-OE) and antisense Sm4CL2-expressing (anti-Sm4CL2) danshen hairy roots over ten years ago. In the follow-up study, we found that total salvianolic acids in Sm4CL2-OE-4 hairy roots increased to 1.35 times of the control-3, and that in anti-Sm4CL2-1 hairy roots decreased to 37.32% of the control-3, but tanshinones in anti-Sm4CL2-1 was accumulated to 1.77 ± 0.16 mg/g of dry weight, compared to undetectable in Sm4CL2-OE-4 and the control-3 hairy roots. Interestingly, Sm4CL2-OE-4 hairy roots contained more lignin, 1.36 times of the control-3, and enhanced cell wall and xylem lignification. Transcriptomic analysis revealed that overexpression of Sm4CL2 caused the upregulation of other phenylpropanoid pathway genes and antisense Sm4CL2 expression resulted in the downregulation of other phenylpropanoid pathway genes but activated the expression of terpenoid pathway genes like SmCYP76AK5, SmGPPS.SSUII.1 and SmDXS2. Protein-protein interaction analysis suggested that Sm4CL2 might interact with PAL, PAL4, CSE, CCoAOMT and SmCYP84A60, and appeared to play a key role in the interaction network. The tracking work in this study proved that Sm4CL2 could redirect both salvianolic acids and tanshinones biosynthesis possibly through synergistically regulating other pathway genes. It also indicated that genetic modification of plant secondary metabolism with biosynthetic gene might cause other responses through protein-protein interactions.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Transcriptional Responses for Biosynthesis of Triterpenoids in Exogenous Inducers Treated Inonotus Hispidus Using RNA-Seq

Inonotus hispidus is a traditional medicinal that grows in Northeast China and produces various economically important compounds, including polysaccharide compounds and terpenoids; triterpenoid saponins is the main bioactive component. Our research group has found that the accumulation of triterpenoid was affected by exogenous inducers. Experimental results showed that treatment with methyl jasmonate (MeJA) and oleic acid significantly increased the triterpenoid content of I. hispidus. However, how exogenous inducers enhance production of secondary metabolites in I. hispidus is not well understood. In this study, metabolite changes were further investigated with UPLC-TOF/MS following exogenous inducer treatment. As a result, a total of eight types of triterpenoids in I. hispidus were identified. The RNA-seq analysis was used to evaluate the effects of exogenous inducers on the expression of triterpenoid-synthesis-related genes in I. hispidus in liquid fermentation. This study is the first exploration to profile the transcriptome of I. hispidus after adding exogenous inducers; the generated data and gene will facilitate further molecular studies on the physiology and metabolism in this fungi. By comparative transcriptomic analysis, a series of candidate genes involved in the biosynthetic pathway of triterpenoids are identified, providing new insights into their biosynthesis at the transcriptome level.

Spatiotemporal and Transcriptional Characterization on Tanshinone Initial Synthesis in Salvia miltiorrhiza Roots

Tanshinones are the bioactive constituents of Danshen (Salvia miltiorrhiza Bunge), which is used in Traditional Chinese Medicine to treat cardiovascular and other diseases, and they synthesize and accumulate in the root periderm of S. miltiorrhiza. However, there is no relevant report on the initial stage of tanshinone synthesis, as well as the root structure and gene expression characteristics. The present study aims to provide new insights into how these bioactive principles begin to synthesize by characterizing possible differences in their biosynthesis and accumulation during early root development from both spatial and temporal aspects. The morphological characteristics and the content of tanshinones in roots of S. miltiorrhiza were investigated in detail by monitoring the seedlings within 65 days after germination (DAGs). The ONT transcriptome sequencing was applied to investigate gene expression patterns. The periderm of the S. miltiorrhiza storage taproot initially synthesized tanshinone on about 30 DAGs. Three critical stages of tanshinone synthesis were preliminarily determined: preparation, the initial synthesis, and the continuous rapid synthesis. The difference of taproots in the first two stages was the smallest, and the differentially expressed genes (DEGs) were mainly enriched in terpene synthesis. Most genes involved in tanshinone synthesis were up regulated during the gradual formation of the red taproot. Plant hormone signal transduction and ABC transport pathways were widely involved in S. miltiorrhiza taproot development. Five candidate genes that may participate in or regulate tanshinone synthesis were screened according to the co-expression pattern. Moreover, photosynthetic ferredoxin (FD), cytochrome P450 reductase (CPR), and CCAAT binding transcription factor (CBF) were predicted to interact with the known downstream essential enzyme genes directly. The above results provide a necessary basis for analyzing the initial synthesis and regulation mechanism of Tanshinones.

The SmMYB36-SmERF6/SmERF115 module regulates the biosynthesis of tanshinones and phenolic acids in salvia miltiorrhiza hairy roots

Tanshinone and phenolic acids are the most important active substances of Salvia miltiorrhiza, and the insight into their transcriptional regulatory mechanisms is an essential process to increase their content in vivo. SmMYB36 has been found to have important regulatory functions in the synthesis of tanshinone and phenolic acid; paradoxically, its mechanism of action in S. miltiorrhiza is not clear. Here, we demonstrated that SmMYB36 functions as a promoter of tanshinones accumulation and a suppressor of phenolic acids through the generation of SmMYB36 overexpressed and chimeric SmMYB36-SRDX (EAR repressive domain) repressor hairy roots in combination with transcriptomic-metabolomic analysis. SmMYB36 directly down-regulate the key enzyme gene of primary metabolism, SmGAPC, up-regulate the tanshinones biosynthesis branch genes SmDXS2, SmGGPPS1, SmCPS1 and down-regulate the phenolic acids biosynthesis branch enzyme gene, SmRAS. Meanwhile, SmERF6, a positive regulator of tanshinone synthesis activating SmCPS1, was up-regulated and SmERF115, a positive regulator of phenolic acid biosynthesis activating SmRAS, was down-regulated. Furthermore, the seven acidic amino acids at the C-terminus of SmMYB36 are required for both self-activating domain and activation of target gene expression. As a consequence, this study contributes to reveal the potential relevance of transcription factors synergistically regulating the biosynthesis of tanshinone and phenolic acid.

Applications of protein engineering in the microbial synthesis of plant triterpenoids

Triterpenoids are a class of natural products widely used in fields related to medicine and health due to their biological activities such as hepatoprotection, anti-inflammation, anti-viral, and anti-tumor. With the advancement in biotechnology, microorganisms have been used as cell factories to produce diverse natural products. Despite the significant progress that has been made in the construction of microbial cell factories for the heterogeneous biosynthesis of triterpenoids, the industrial production of triterpenoids employing microorganisms has been stymied due to the shortage of efficient enzymes as well as the low expression and low catalytic activity of heterologous proteins in microbes. Protein engineering has been demonstrated as an effective way for improving the specificity, catalytic activity, and stability of the enzyme, which can be employed to overcome these challenges. This review summarizes the current progress in the studies of Oxidosqualene cyclases (OSCs), cytochrome P450s (P450s), and UDP-glycosyltransferases (UGTs), the key enzymes in the triterpenoids synthetic pathway. The main obstacles restricting the efficient catalysis of these key enzymes are analyzed, the applications of protein engineering for the three key enzymes in the microbial synthesis of triterpenoids are systematically reviewed, and the challenges and prospects of protein engineering are also discussed.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

Molecular Mechanism of Tanshinone against Prostate Cancer

Prostate cancer (PCa) is the most common malignant tumor of the male urinary system in Europe and America. According to the data in the World Cancer Report 2020, the incidence rate of PCa ranks second in the prevalence of male malignant tumors and varies worldwide between regions and population groups. Although early PCa can achieve good therapeutic results after surgical treatment, due to advanced PCa, it can adapt and tolerate androgen castration-related drugs through a variety of mechanisms. For this reason, it is often difficult to achieve effective therapeutic results in the treatment of advanced PCa. Tanshinone is a new fat-soluble phenanthraquinone compound derived from Salvia miltiorrhiza that can play a therapeutic role in different cancers, including PCa. Several studies have shown that Tanshinone can target various molecular pathways of PCa, including the signal transducer and activator of transcription 3 (STAT3) pathway, androgen receptor (AR) pathway, phosphatidylinositol-3-kinase (PI3K)/protein kinase B (Akt)/mammalian target of rapamycin (mTOR) pathway, and mitogen-activated protein kinase (MAPK) pathway, which will affect the release of pro-inflammatory cytokines and affect cell proliferation, apoptosis, tumor metabolism, genomic stability, and tumor drug resistance. Thus, the occurrence and development of PCa cells are inhibited. In this review, we summarized the in vivo and in vitro evidence of Tanshinone against prostate cancer and discussed the effect of Tanshinone on nuclear factor kappa-B (NF-κB), AR, and mTOR. At the same time, we conducted a network pharmacology analysis on the four main components of Tanshinone to further screen the possible targets of Tanshinone against prostate cancer and provide ideas for future research.

The sage genome provides insight into the evolutionary dynamics of diterpene biosynthesis gene cluster in plants

The widely cultivated medicinal and ornamental plant sage (Salvia officinalis L.) is an evergreen shrub of the Lamiaceae family, native to the Mediterranean. We assembled a high-quality sage genome of 480 Mb on seven chromosomes, and identified a biosynthetic gene cluster (BGC) encoding two pairs of diterpene synthases (diTPSs) that, together with the cytochromes P450 (CYPs) genes located inside and outside the cluster, form two expression cascades responsible for the shoot and root diterpenoids, respectively, thus extending BGC functionality from co-regulation to orchestrating metabolite production in different organs. Phylogenomic analysis indicates that the Salvia clades diverged in the early Miocene. In East Asia, most Salvia species are herbaceous and accumulate diterpenoids in storage roots. Notably, in Chinese sage S. miltiorrhiza, the diterpene BGC has contracted and the shoot cascade has been lost. Our data provide genomic insights of micro-evolution of growth type-associated patterning of specialized metabolite production in plants.

Transcriptome of the coralline alga Calliarthron tuberculosum (Corallinales, Rhodophyta) reveals convergent evolution of a partial lignin biosynthesis pathway

The discovery of lignins in the coralline red alga Calliarthron tuberculosum raised new questions about the deep evolution of lignin biosynthesis. Here we present the transcriptome of C. tuberculosum supported with newly generated genomic data to identify gene candidates from the monolignol biosynthetic pathway using a combination of sequence similarity-based methods. We identified candidates in the monolignol biosynthesis pathway for the genes 4CL, CCR, CAD, CCoAOMT, and CSE but did not identify candidates for PAL, CYP450 (F5H, C3H, C4H), HCT, and COMT. In gene tree analysis, we present evidence that these gene candidates evolved independently from their land plant counterparts, suggesting convergent evolution of a complex multistep lignin biosynthetic pathway in this red algal lineage. Additionally, we provide tools to extract metabolic pathways and genes from the newly generated transcriptomic and genomic datasets. Using these methods, we extracted genes related to sucrose metabolism and calcification. Ultimately, this transcriptome will provide a foundation for further genetic and experimental studies of calcifying red algae.

Functional Characterization of a 2OGD Involved in Abietane-Type Diterpenoids Biosynthetic Pathway in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most commonly used Chinese medicinal herbs. Tanshinones, the most abundant lipid-soluble bioactive constituents of S. miltiorrhiza, are a class of structural highly oxidized abietane-type diterpenoids with multiple pharmacological activities. Although several enzymes, including diterpene synthase, cytochrome P450, and Fe(II)/2-oxoglutarate-dependent dioxygenase (2OGD), have been functionally characterized in biosynthesis of abietane-type diterpenoids, the highly oxidized structure and complex secondary metabolic network of tanshinones imply that more oxidases should be characterized. Here, we identified a new 2OGD (Sm2OGD25) from S. miltiorrhiza. Molecular cloning and functional studies in vitro showed that Sm2OGD25 could catalyze the hydroxylation of sugiol at C-15 and C-16 positions to produce hypargenin B and crossogumerin C, respectively. The phylogenetic analysis of the DOXC family demonstrated that Sm2OGD25 belongs to the DOXC54 clade. Furthermore, structural modeling and site-directed mutagenesis characterization revealed the importance of the hydrogen-bonding residue Y339 and the hydrophobic residues (V122, F129, A144, A208, F303, and L344) in substrate binding and enzyme activity. This study will promote further studies on the catalytic characterization of plant 2OGDs and the secondary metabolic biosynthesis network of diterpenoids.

Identification of Abietane-Type Diterpenoids and Phenolic Acids Biosynthesis Genes in Salvia apiana Jepson Through Full-Length Transcriptomic and Metabolomic Profiling

Salvia apiana (S. apiana) Jepson is a medicinal plant that is frequently used by the Chumash Indians in southern California as a diaphoretic, calmative, diuretic, or antimicrobial agent. Abietane-type diterpenoids (ATDs) and phenolic acids (PAs) are the main bioactive ingredients in S. apiana. However, few studies have looked into the biosynthesis of ATDs and PAs in S. apiana. In this study, using metabolic profiling focused on the ATDs and PAs in the roots and leaves of S. apiana, we found a distinctive metabolic feature with all-around accumulation of ATDs, but absence of salvianolic acid B. To identify the candidate genes involved in these biosynthesis pathways, full-length transcriptome was performed by PacBio single-molecule real-time (SMRT) sequencing. A total of 50 and 40 unigenes were predicted to be involved in ATDs and PAs biosynthesis, respectively. Further transcriptional profile using Illumina HiSeq sequencing showed that the transcriptional variations of these pathways were consistent with the accumulation patterns of corresponding metabolites. A plant kingdom-wide phylogenetic analysis of cytochromes (CYPs) identified two CYP76AK and two CYP76AH subfamily genes that might contribute for the specific ATDs biosynthesis in S. apiana. We also noticed that the clade VII laccase gene family was significantly expanded in Salvia miltiorrhiza compared with that of S. apiana, indicating their involvements in the formation of salvianolic acid B. In conclusion, our results will enable the further understanding of ATDs and PAs biosynthesis in S. apiana and Salvia genus.

One-Step Regeneration of Hairy Roots to Induce High Tanshinone Plants in Salvia miltiorrhiza

Salvia miltiorrhiza is a traditional Chinese medicinal plant of Labiatae, which has been widely utilized to treat a variety of cardiovascular and cerebrovascular diseases. However, due to the long growth cycle, low content of active ingredients, and serious quality deterioration of S. miltiorrhiza, the use of biotechnology to improve S. miltiorrhiza to meet the growing demand for clinical applications has become a research hotspot. In this study, a novel one-step hairy root regeneration method was developed, which could rapidly obtain hairy roots and regenerated plants with high tanshinone content. By optimizing the parameters of Agrobacterium rhizogenes transformation in S. miltiorrhiza, it was finally established that the explants were infected in Ar.qual (OD₆₀₀ = 0.6) for 10 min, co-cultured for 3 days, and then screened on the screening medium containing 7.5 mg/l hygromycin, the maximum transformation frequency can reach 73.85%. GFP and PCR detection yielded a total of 9 positive transgenic hairy root lines and 11 positive transgenic regenerated plants. SmGGPPS1 was successfully overexpressed in positive transgenic regenerated plants, according to the results of qRT-PCR. The content of tanshinone IIA and cryptotanshinone were dramatically enhanced in transgenic regenerated plants and hairy roots by Ultra Performance Liquid Chromatography analysis. Based on the Agrobacterium-mediated transformation of S. miltiorrhiza, this study developed a new method for regenerating plants with transgenic hairy roots. This method provides a foundation for the breeding of S. miltiorrhiza and the sustainable development of medicinal plant resources, as well as provides a useful reference for the application of other species.