Functional Diversification of Kaurene Synthase-Like Genes in Isodon rubescens

Gibberellin-biosynthetic ent-kaurene synthases in higher plants do not require their non-catalytic domains for the catalysis

ent-Kaurene is a biosynthetic intermediate diterpene of phytohormone gibberellins, and is biosynthesized from geranylgeranyl diphosphate via ent-copalyl diphosphate (ent-CDP). The successive cyclization is catalyzed by two distinct diterpene synthases, ent-CDP synthase (ent-CPS) and ent-kaurene synthase (KS). Homologs of these diterpene synthase genes have been reported to be involved in the biosynthesis of specialized-metabolic diterpenoids for defense in several plant species, including rice (Oryza sativa). These diterpene synthases consist of three domains, αβγ domains. Active sites of ent-CPS exist at the interface of β and γ domain, while those of KS are located within the α domain. We herein carried out domain-deletion experiments using several KSs and KS like enzymes (KSLs) to obtain insights into the roles of domains other than active-site domains. As previously reported in taxadiene synthase, deletion of γ or βγ domains drastically decreased activities of specialized-metabolic OsKSL5, OsKSL8, OsKSL7 and OsKSL10 in O. sativa. However, unexpectedly, only α domains of several gibberellin-biosynthetic KSs, including OsKS1 in O. sativa, AtKS in Arabidopsis thaliana, TaKS in wheat (Triticum aestivum) and BdKS1 in Brachypodium distachyon, retained their original functions. Additionally, the specialized-metabolic OsKSL4, which is closely related to OsKS1, also functioned without its βγ domains. Domain-swapping experiments showed that replacing βγ domains in OsKSL7 with those from other KS/KSLs retained the OsKSL7 activity. Moreover, deletion of βγ domains of bifunctional PpCPS/KS in moss (Physcomitrella patens) drastically impaired its KS-related activity. Thus, we demonstrate that monofunctional gibberellin-biosynthetic KSs are the unique diterpene synthases that retain their functions without βγ domains.

Heterologous Biosynthesis of Kauralexin A1 in Saccharomyces cerevisiae through Metabolic and Enzyme Engineering

Kauralexin A1 (KA1) is a key intermediate of the kauralexin A series metabolites of maize phytoalexins. However, their application is severely limited by their low abundance in maize. In this study, an efficient biosynthetic pathway was constructed to produce KA1 in Saccharomyces cerevisiae. Also, metabolic and enzyme engineering strategies were applied to construct the high-titer strains, such as chassis modification, screening synthases, the colocalization of enzymes, and multiple genomic integrations. First, the KA1 precursor ent-kaurene was synthesized using the efficient diterpene synthase GfCPS/KS from Fusarium fujikuroi, and optimized to reach 244.36 mg/L in shake flasks, which displayed a 200-fold increase compared to the initial strain. Then, the KA1 was produced under the catalysis of ZmCYP71Z18 from Zea mays and SmCPR1 from Salvia miltiorrhiza, and the titer was further improved by integrating the fusion protein into the genome. Finally, an ent-kaurene titer of 763.23 mg/L and a KA1 titer of 42.22 mg/L were achieved through a single-stage fed-batch fermentation in a 5 L bioreactor. This is the first report of the heterologous biosynthesis of maize diterpene phytoalexins in S. cerevisiae, which lays a foundation for further pathway reconstruction and biosynthesis of the kauralexin A series maize phytoalexins.

Multi-omics analysis of small RNA, transcriptome, and degradome to identify putative miRNAs linked to MeJA regulated and oridonin biosynthesis in Isodon rubescens

Isodon rubescens has garnered much attention due to its anti-tumor or anti-cancer properties. However, little is known about the molecular mechanism of oridonin biosynthesis leveraging the regulatory network between small RNAs and mRNAs. In this study, the regulatory networks of miRNAs and targets were examined by combining mRNA, miRNA, and degradome. A total of 348 miRNAs, including 287 known miRNAs and 61 novel miRNAs, were identified. Among them, 51 miRNAs were significantly expressed, and 36 miRNAs responded to MeJA. A total of 3066 target genes were associated with 228 miRNAs via degradome sequencing. Multi-omics analysis demonstrated that 27 miRNA-mRNA pairs were speculated to be involved in MeJA regulation, and 36 miRNA-mRNA pairs were hypothesized to be involved in the genotype-dependence of I. rubescens. Furthermore, 151 and 7 miRNA-mRNA modules were likely engaged in oridonin biosynthesis as identified by psRNATarget and degradome sequencing, respectively. Some miRNA-mRNA modules were confirmed via RT-qPCR. Moreover, miRNAs targeting plant hormone signal transduction pathway genes were identified, such as miR156, miR167, miR393, and PC-3p-19822_242. Collectively, our results demonstrate for the first time that miRNAs are identified in I. rubescens, and laid a solid foundation for further research on the molecular mechanism of oridonin biosynthesis mediated by miRNA.

Sequence-Structure Analysis Unlocking the Potential Functional Application of the Local 3D Motifs of Plant-Derived Diterpene Synthases

Plant-derived diterpene synthases (PdiTPSs) play a critical role in the formation of structurally and functionally diverse diterpenoids. However, the specificity or functional-related features of PdiTPSs are not well understood. For a more profound insight, we collected, constructed, and curated 199 functionally characterized PdiTPSs and their corresponding 3D structures. The complex correlations among their sequences, domains, structures, and corresponding products were comprehensively analyzed. Ultimately, our focus narrowed to the geometric arrangement of local structures. We found that local structural alignment can rapidly localize product-specific residues that have been validated by mutagenesis experiments. Based on the 3D motifs derived from the residues around the substrate, we successfully searched diterpene synthases (diTPSs) from the predicted terpene synthases and newly characterized PdiTPSs, suggesting that the identified 3D motifs can serve as distinctive signatures in diTPSs (I and II class). Local structural analysis revealed the PdiTPSs with more conserved amino acid residues show features unique to class I and class II, whereas those with fewer conserved amino acid residues typically exhibit product diversity and specificity. These results provide an attractive method for discovering novel or functionally equivalent enzymes and probing the product specificity in cases where enzyme characterization is limited.

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.

Tandem duplication and sub-functionalization of clerodane diterpene synthase originate the blooming of clerodane diterpenoids in Scutellaria barbata

Scutellaria barbata is a traditional Chinese herb medicine and a major source of bioactive clerodane diterpenoids. However, barely clerodanes have been isolated from the closely related S. baicalensis. Here we assembled a chromosome-level genome of S. barbata and identified three class II clerodane diterpene synthases (SbarKPS1, SbarKPS2 and SbaiKPS1) from these two organisms. Using in vitro and in vivo assays, SbarKPS1 was characterized as a monofunctional (-)-kolavenyl diphosphate synthases ((-)-KPS), while SbarKPS2 and SbaiKPS1 produced major neo-cleroda-4(18),13E-dienyl diphosphate with small amount of (-)-KPP. SbarKPS1 and SbarKPS2 shared a high protein sequence identity and formed a tandem gene pair, indicating tandem duplication and sub-functionalization probably led to the evolution of monofunctional (-)-KPS in S. barbata. Additionally, SbarKPS1 and SbarKPS2 were primarily expressed in the leaves and flowers of S. barbata, which was consistent with the distribution of major clerodane diterpenoids scutebarbatine A and B. In contrast, SbaiKPS1 was barely expressed in any tissue of S. baicalensis. We further explored the downstream class I diTPS by functional characterizing of SbarKSL3 and SbarKSL4. Unfortunately, no dephosphorylated product was detected in the coupled assays with SbarKSL3/KSL4 and four class II diTPSs (SbarKPS1, SbarKPS2, SbarCPS2 and SbarCPS4) when a phosphatase inhibitor cocktail was included. Co-expression of SbarKSL3/KSL4 with class II diTPSs in yeast cells did not increase the yield of the corresponding dephosphorylated products, either. Together, these findings elucidated the involvement of two class II diTPSs in clerodane biosynthesis in S. barbata, while the class I diTPS is likely not responsible for the subsequent dephosphorylation step.

Functional diversity of diterpene synthases in Aconitum plants

Members of the Aconitum genus within the Ranunculaceae family are known to accumulate a broad array of medicinal and toxic diterpenoid alkaloids (DAs). Historically, ent-copalyl diphosphate (ent-CPP) was considered the sole precursor in DAs biosynthesis. However, the recent discovery of ent-8,13-CPP synthase in A. gymnandrum Maxim., which participates in ent-atiserene biosynthesis, raises the question of whether this gene is conserved throughout the Aconitum genus. In this study, RNA sequencing and PacBio Iso-sequencing were employed to identify diterpene synthases (diTPSs) in four additional Aconitum species with distinct DA compositions. In vitro and in vivo analyses functionally characterized a diverse array of 10 class II and 9 class I diTPSs. In addition to the identification of seven class II diTPSs as ent-CPP synthases, three other synthases generating ent-8,13-CPP, 8,13-CPP, and 8α-hydroxy-CPP were also discovered. Four class I kaurene synthases-like (KSLs) were observed to react with ent-CPP to yield ent-kaurene. Three KSLs not only reacted with ent-CPP but also ent-8,13-CPP to produce ent-atiserene. AsiKSL2-1 was found to react with 8α-hydroxy-CPP to produce Z-abienol and AsiKSL2-2 exhibited no activity with any of the four intermediates. This research delineates the known diterpene biosynthesis pathways in six Aconitum species and explores the highly divergent diterpene synthases within the genus, which are consistent with their phylogeny and may be responsible for the differential distribution of diterpenoid alkaloids in root and aerial parts. These findings contribute valuable insights into the diversification of diterpene biosynthesis and establish a solid foundation for future investigation into DA biosynthetic pathways in Aconitum.

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.

Structural and mechanistic insights into the precise product synthesis by a bifunctional miltiradiene synthase

Selaginella moellendorffii miltiradiene synthase (SmMDS) is a unique bifunctional diterpene synthase (diTPS) that catalyses the successive cyclization of (E,E,E)-geranylgeranyl diphosphate (GGPP) via (+)-copalyl diphosphate (CPP) to miltiradiene, which is a crucial precursor of important medicinal compounds, such as triptolide, ecabet sodium and carnosol. Miltiradiene synthetic processes have been studied in monofunctional diTPSs, while the precise mechanism by which active site amino acids determine product simplicity and the experimental evidence for reaction intermediates remain elusive. In addition, how bifunctional diTPSs work compared to monofunctional enzymes is attractive for detailed research. Here, by mutagenesis studies of SmMDS, we confirmed that pimar-15-en-8-yl⁺ is an intermediate in miltiradiene synthesis. Moreover, we determined the apo-state and the GGPP-bound state crystal structures of SmMDS. By structure analysis and mutagenesis experiments, possible contributions of key residues both in class I and II active sites were suggested. Based on the structural and functional analyses, we confirmed the copal-15-yl⁺ intermediate and unveiled more details of the catalysis process in the SmMDS class I active site. Moreover, the structural and experimental results suggest an internal channel for (+)-CPP produced in the class II active site moving towards the class I active site. Our research is a good example for intermediate identification of diTPSs and provides new insights into the product specificity determinants and intermediate transport, which should greatly facilitate the precise controlled synthesis of various diterpenes.

Functional Characterization of Three Diterpene Synthases Responsible for Tetracyclic Diterpene Biosynthesis in Scoparia dulcis

Scoparia dulcis produces unique biologically active diterpenoids such as scopadulcic acid B (SDB). They are biosynthesized from geranylgeranyl diphosphate (GGPP) via syn-copalyl diphosphate (syn-CPP) and scopadulanol as an important key intermediate. In this paper, we functionally characterized three diterpene synthases, SdCPS2, SdKSL1 and SdKSL2, from S. dulcis. The SdCPS2 catalyzed a cyclization reaction from GGPP to syn-CPP, and SdKSL1 did from syn-CPP to scopadulan-13α-ol. On the other hand, SdKSL2 was found to incorporate a non-sense mutation at 682. Therefore, we mutated the nucleotide residue from A to G in SdKSL2 to produce SdKSL2mut, and it was able to recover the catalytic function from syn-CPP to syn-aphidicol-16-ene, the precursor to scopadulin. From our results, SdCPS2 and SdKSL1 might be important key players for SDB biosynthesis in S. dulcis.

Validation of suitable reference genes by various algorithms for gene expression analysis in Isodon rubescens under different abiotic stresses

Isodon rubescens (Hemsley) H. Hara (Lamiaceae) is a traditional Chinese medicine plant that has been used to treat various human diseases. Oridonin is one of the main active ingredients, and the route of its molecular biosynthesis remains to be determined. The study of gene expression patterns can provide clues toward the understanding of its biological functions. The selection of suitable reference genes for normalizing target gene expression is the first steps in any quantitative real-time PCR (RT-qPCR) gene expression study. Therefore, validation of suitable reference genes is necessary for obtaining reliable results in RT-qPCR analyses of I. rubescens. Here, 12 candidate reference genes were chosen, and their expression stability in different tissues of I. rubescens and in leaves under different abiotic stresses (NaCl, dehydration, SA, MeJA, and ABA) was evaluated using the ∆Ct, NormFinder, GeNorm, BestKeeper, and RankAggreg statistical tools. Analysis using the comprehensive tools of RankAggreg algorithm showed that GADPH, 18S and eIF were stably expressed in different tissues; UBQ, Apt, and HIS; Cycl, UBQ, and PP2A; GADPH, 18S, and eIF; eIF, UBQ, and PP2A; TUB, Cycl, and UBQ; were the best three candidate reference genes for the samples of Dehydration, NaCl, SA, MeJA, and ABA treatment, respectively. While for the concatenated sets of ND (NaCl and dehydration) and SMA (SA, MeJA, and ABA), UBQ, HIS, and TUA; UBQ, eIF and Apt were the three appropriate candidate reference genes, respectively. In addition, the expression patterns of HMGR in different tissues and under different treatments were used to confirm the reliability of the selected reference genes, indicating that the use of an inappropriate reference gene as the internal control will cause results with a large deviation. This work is the first study on the expression stability of reference genes in I. rubescens and will be particularly useful for gene functional research in this species.

Integrating Metabolomics and Transcriptomics to Unveil Atisine Biosynthesis in Aconitum gymnandrum Maxim

Diterpene alkaloids (DAs) are characteristic compounds in Aconitum, which are classified into four skeletal types: C18, C19, C20, and bisditerpenoid alkaloids. C20-DAs are thought to be the precursor of the other types. Their biosynthetic pathway, however, is largely unclear. Herein, we combine metabolomics and transcriptomics to unveil the methyl jasmonate (MJ) inducible biosynthesis of DAs in the sterile seedling of A. gymnandrum, the only species in the Subgenus Gymnaconitum (Stapf) Rapaics. Target metabolomics based on root and aerial portions identified 51 C19-DAs and 15 C20-DAs, with 40 inducible compounds. The highest content of C20-DA atisine was selected for further network analysis. PacBio Isoform sequencing integrated with RNA sequencing not only provided the full-length transcriptome but also their response to induction, revealing 1994 genes that exhibited up-regulated expression. Further, 38 genes involved in terpenoid biosynthesis were identified, including 7 diterpene synthases. In addition to the expected function of the four diterpene synthases, AgCPS5 was identified to be a new ent-8,13-CPP synthase in Aconitum and could also combine with AgKSL1 to form the C20-DAs precursor ent-atiserene. Combined with multiple network analyses, six CYP450 and seven 2-ODD genes predicted to be involved in the biosynthesis of atisine were also identified. This study not only sheds light on diterpene synthase evolution in Aconitum but also provides a rich dataset of full-length transcriptomes, systemic metabolomes, and gene expression profiles, setting the groundwork for further investigation of the C20-DAs biosynthesis pathway.

Design, Synthesis and Biological Evaluation of Lophanic Acid Derivatives as Antifungal and Antibacterial Agents

In order to discover more promising antifungal and antibacterial agents, a series of new derivatives were designed and synthesized by structure modification based on the naturally occurring antimicrobial compound lophanic acid. The structures of all the target compounds were well characterized by spectroscopic data. The stereochemistry of these compounds was further determined through the X-ray diffraction analysis of 6a. The synthetic compounds were evaluated for their antimicrobial activities against filamentous fungi (T. rubrum, T. mentagrophytes), yeasts (C. neoformans, C. albicans) and Gram-positive and Gram-negative bacteria (MRSA, S. mutans, S. sobrinus, and E. coli). Among them, 3d and 3i are found as the most promising leads that showed potent inhibitory effects against all the tested fungal and bacterial strains except for E. coli. The presence of the C-20 carboxylic ester groups and the free hydroxy group at C-13 was found to be essential for the antifungal and antibacterial activities of the lophanic acid derivatives.

Molecular Cloning and Functional Analysis of IrUGT86A1-like Gene in Medicinal Plant Isodon rubescens (Hemsl.) Hara

The synthesis of secondary metabolites in plants often includes glycosylation modifications. Often, the final step of constructing plant secondary metabolites is completed by glycosylation transferases, which are also involved in many cell processes. In this study, a UDP-glycosyltransferase gene (UGT) was amplified from Isodon rubescens (Hemsl.) Hara with RT-PCR and named IrUGT86A1-like (GenBank: MZ913258). Here, we found that IrUGT86A1-like gene is 1450 bp in length and encodes for 479 amino acids. Bioinformatics analysis revealed that IrUGT86A1-like is a stable and hydrophilic protein, located in the cytoplasm with a transmembrane domain. Phylogenetic analysis showed that IrUGT86A1-like protein has the closest genetic relationship with the UDP-glycosyltransferase 86A1-like protein (XP_042054241.1) of Salvia splendens. RT-qPCR analysis demonstrated that the expression of IrUGT86A1-like gene varied in different tissues; leaves had the highest expression followed by flowers, stems, and roots had the lowest expression. This expression trend is similar to the distribution of oridonin content in different tissues of I. rubescens. Additionally, IrUGT86A1-like gene was found to be positively enhanced by NaCl and MeJA treatment, and in contrast was down-regulated by ABA treatment. Finally, the prokaryotic expression vector pEASY®-Blunt E1-IrUGT86A1 was successfully used to express about 53 KD of IrUGT86A1-like protein. This research builds a foundation for further investigation on the function of this gene in the synthesis and modification of secondary metabolites.

Diterpene synthases from Leonurus japonicus elucidate epoxy-bridge formation of spiro-labdane diterpenoids

Spiro-9,13-epoxy-labdane diterpenoids are commonly found in Leonurus species, particularly in Leonurus japonicus Houtt., which is a medicinal herb of long-standing use in Asia and in which such spiro-heterocycles are present in at least 38 diterpenoids. Here, through generation of a transcriptome and functional characterization of six diterpene synthases (diTPSs) from L. japonicus, including three class II diTPSs (LjTPS1, LjTPS3, and LjTPS4) and three class I diTPSs (LjTPS5, LjTPS6, and LjTPS7), formation of the spiro-9,13-epoxy-labdane backbone was elucidated, along with identification of the relevant diTPSs for production of other labdane-related diterpenes. Similar to what has been found with diTPSs from other plant species, while LjTPS3 specifically produces the carbon-9 (C9) hydroxylated bicycle peregrinol diphosphate (PPP), the subsequently acting LjTPS6 yields a mixture of four products, largely labda-13(16),14-dien-9-ol, but with substantial amounts of viteagnusin D and the C13-S/R epimers of 9,13-epoxy-labda-14-ene. Notably, structure-function analysis identified a critical residue in LjTPS6 (I420) in which single site mutations enable specific production of the 13S epimer. Indeed, extensive mutagenesis demonstrated that LjTPS6:I420G reacts with PPP to both specifically and efficiently produce 9,13S-epoxy-labda-14-ene, providing a specialized synthase for further investigation of derived diterpenoid biosynthesis. The results reported here provide a strong foundation for future studies of the intriguing spiro-9,13-epoxy-labdane diterpenoid metabolism found in L. japonicus.

Mechanistic analysis for the origin of diverse diterpenes in Tripterygium wilfordii

Tripterygium wilfordii is a valuable medicinal plant rich in biologically active diterpenoids, but there are few studies on the origins of these diterpenoids in its secondary metabolism. Here, we identified three regions containing tandemly duplicated diterpene synthase genes on chromosomes (Chr) 17 and 21 of T. wilfordii and obtained 11 diterpene synthases with different functions. We further revealed that these diterpene synthases underwent duplication and rearrangement at approximately 2.3–23.7 million years ago (MYA) by whole-genome triplication (WGT), transposon mediation, and tandem duplication, followed by functional divergence. We first demonstrated that four key amino acids in the sequences of TwCPS3, TwCPS5, and TwCPS6 were altered during evolution, leading to their functional divergence and the formation of diterpene secondary metabolites. Then, we demonstrated that the functional divergence of three TwKSLs was driven by mutations in two key amino acids. Finally, we discovered the mechanisms of evolution and pseudogenization of miltiradiene synthases in T. wilfordii and elucidated that the new function in TwMS1/2 from the terpene synthase (TPS)-b subfamily was caused by progressive changes in multiple amino acids after the WGT event. Our results provide key evidence for the formation of diverse diterpenoids during the evolution of secondary metabolites in T. wilfordii.

Discovery and Functional Characterization of a Diverse Diterpene Synthase Family in the Medicinal Herb Isodon lophanthoides Var. gerardiana

Isodon lophanthoides var. gerardiana (Lamiaceae), also named xihuangcao, is a traditional Chinese medicinal herb that exhibits a broad range of pharmacological activities. Abietane-type diterpenoids are the characteristic constituents of I. lophanthoides, yet their biosynthesis has not been elucidated. Although the aerial parts are the most commonly used organs of I. lophanthoides, metabolite profiling by gas chromatography-mass spectrometry showed the underground parts also contain large amounts of labdane diterpenoids including abietatriene, miltiradiene and ferruginol, which are distinct from the 13-hydroxy-8(14)-abietene detected in the aerial parts. Comparative transcriptome analysis of root and leaf samples identified a diverse diterpene synthase family including 6 copalyl diphosphate synthase (IlCPS1-6) and 5 kaurene synthase-like (IlKSL1-5). Here we report the functional characterization of six of these enzymes using yeast heterologous expression system. Both IlCPS1 and IlCPS3 synthesized (+)-copalyl diphosphate (CPP), in combination with IlKSL1 resulted in miltiradiene, precursor of abietane-type diterpenoids, while coupling with IlKSL5 led to the formation of hydroxylated diterpene scaffold nezukol. Expression profiling and phylogenetic analysis further support the distinct evolutionary relationship and spatial distribution of IlCPS1 and IlCPS3. IlCPS2 converted GGPP into labda-7,13E-dien-15-ol diphosphate. IlCPS6 was identified as ent-CPS, indicating a role in gibberellin metabolism. We further identified a single residue that determined the water addition of nezukol synthase IlKSL5. Substitution of alanine 513 with isoleucine completely altered the product outcome from hydroxylated nezukol to isopimara-7,15-diene. Together, these findings elucidated the early steps of bioactive abietane-type diterpenoid biosynthesis in I. lophanthoides and the catalytic mechanism of nezukol synthase.

Functional identification of the terpene synthase family involved in diterpenoid alkaloids biosynthesis in Aconitum carmichaelii

Aconitum carmichaelii is a high-value medicinal herb widely used across China, Japan, and other Asian countries. Aconitine-type diterpene alkaloids (DAs) are the characteristic compounds in Aconitum. Although six transcriptomes, based on short-read next generation sequencing technology, have been reported from the Aconitum species, the terpene synthase (TPS) corresponding to DAs biosynthesis remains unidentified. We apply a combination of Pacbio isoform sequencing and RNA sequencing to provide a comprehensive view of the A. carmichaelii transcriptome. Nineteen TPSs and five alternative splicing isoforms belonging to TPS-b, TPS-c, and TPS-e/f subfamilies were identified. In vitro enzyme reaction analysis functional identified two sesqui-TPSs and twelve diTPSs. Seven of the TPS-c subfamily genes reacted with GGPP to produce the intermediate ent-copalyl diphosphate. Five AcKSLs separately reacted with ent-CPP to produce ent-kaurene, ent-atiserene, and ent-13-epi-sandaracopimaradie: a new diterpene found in Aconitum. AcTPSs gene expression in conjunction DAs content analysis in different tissues validated that ent-CPP is the sole precursor to all DAs biosynthesis, with AcKSL1, AcKSL2s and AcKSL3-1 responsible for C20 atisine and napelline type DAs biosynthesis, respectively. These data clarified the molecular basis for the C20-DAs biosynthetic pathway in A. carmichaelii and pave the way for further exploration of C19-DAs biosynthesis in the Aconitum species.

A review: biosynthesis of plant-derived labdane-related diterpenoids

Plant-derived labdane-related diterpenoids (LRDs) represent a large group of terpenoids. LRDs possess either a labdane-type bicyclic core structure or more complex ring systems derived from labdane-type skeletons, such as abietane, pimarane, kaurane, etc. Due to their various pharmaceutical activities and unique properties, many of LRDs have been widely used in pharmaceutical, food and perfume industries. Biosynthesis of various LRDs has been extensively studied, leading to characterization of a large number of new biosynthetic enzymes. The biosynthetic pathways of important LRDs and the relevant enzymes (especially diterpene synthases and cytochrome P450 enzymes) were summarized in this review.

Functional kaurene-synthase-like diterpene synthases lacking a gamma domain are widely present in Oryza and related species

Various diterpene synthases have been functionally identified in cultivated rice (Oryza sativa). These are the homologs of ent-copalyl diphosphate (ent-CDP) synthase and ent-kaurene synthase (KS) that are responsible for the biosynthesis of gibberellins, diterpenoid phytohormones. We isolated a cDNA encoding full-length OsKSL12, a previously uncharacterized KS like (KSL) enzyme that consists of a β-domain and an α-domain with an active center, but lacks an N-terminal γ-domain. Functional analysis using a bacterial expression system showed that recombinant OsKSL12 converted ent-CDP into ent-manool or ent-13-epi-manool. Comparative genomics revealed that functional OsKSL12 homologs exist in diverse wild species in the Oryzeae- Oryza nivara (Oryza rufipogon), Oryza coarctata, Oryza granulata, Leersia perrieri and Leersia tisseranti. KSL12 homologs in O. granulata, L. perrieri and L. tisseranti preferentially reacted with GGDP rather than ent-CDP, resulting in geranyllinalool rather than ent-manool or ent-13-epi-manool as the main product, meaning that KSL12 functionally diversified during evolution in the Oryzeae.

Recent progress and new perspectives for diterpenoid biosynthesis in medicinal plants

Diterpenoids, including more than 18,000 compounds, represent an important class of metabolites that encompass both phytohormones and some industrially relevant compounds. These molecules with complex, diverse structures and physiological activities, have high value in the pharmaceutical industry. Most medicinal diterpenoids are extracted from plants. Major advances in understanding the biosynthetic pathways of these active compounds are providing unprecedented opportunities for the industrial production of diterpenoids by metabolic engineering and synthetic biology. Here, we summarize recent developments in the field of diterpenoid biosynthesis from medicinal herbs. An overview of the pathways and known biosynthetic enzymes is presented. In particular, we look at the main findings from the past decade and review recent progress in the biosynthesis of different groups of ringed compounds. We also discuss diterpenoid production using synthetic biology and metabolic engineering strategies, and draw on new technologies and discoveries to bring together many components into a useful framework for diterpenoid production.

Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis in Salvia miltiorrhiza

Tanshinones are the bioactive nor-diterpenoid constituents of the Chinese medicinal herb Danshen (Salvia miltiorrhiza). These groups of chemicals have the characteristic furan D-ring, which differentiates them from the phenolic abietane-type diterpenoids frequently found in the Lamiaceae family. However, how the 14,16-epoxy is formed has not been elucidated. Here, we report an improved genome assembly of Danshen using a highly homozygous genotype. We identify a cytochrome P450 (CYP71D) tandem gene array through gene expansion analysis. We show that CYP71D373 and CYP71D375 catalyze hydroxylation at carbon-16 (C16) and 14,16-ether (hetero)cyclization to form the D-ring, whereas CYP71D411 catalyzes upstream hydroxylation at C20. In addition, we discover a large biosynthetic gene cluster associated with tanshinone production. Collinearity analysis indicates a more specific origin of tanshinones in Salvia genus. It illustrates the evolutionary origin of abietane-type diterpenoids and those with a furan D-ring in Lamiaceae.

ent-Atisane diterpenoids: isolation, structure and bioactivity

Covering: up to 2020 ent-Atisane diterpenoids are a class of over 150 members with diverse structures and valuable bioactivities. These compounds share a curious history in which the synthesis of the archetypal member preceded its isolation from natural sources. In this review, we provide a comprehensive summary of the isolation, structure, and bioactivity of ent-atisane diterpenoids from their discovery in 1965 to the present day.

Pathway elucidation and engineering of plant-derived diterpenoids

Plant-derived diterpenoids are indispensable to plant development, stress-resistance and interaction with environmental microorganisms. Besides significant roles in plant fitness and adaption, many bioactivities beneficial to human beings are also found in diterpenoids from terrestrial plants. However, these high-value compounds are always present in limited species with low-abundance. Complicated chemosynthesis hardly meets the needs of sufficient supplies. To overcome these obstacles, it is necessary to investigate how diterpenoids are biosynthesized in planta, and followed by engineering the biosynthetic pathway to achieve high yield production. This review will summarize the recent progress of plant diterpenoid biosynthetic pathway discovery and engineering, hoping to offer an inspiration for concerned researchers.

Divergent Evolution of the Diterpene Biosynthesis Pathway in Tea Plants (Camellia sinensis) Caused by Single Amino Acid Variation of ent-Kaurene Synthase

Most plant terpenoids are classified as secondary metabolites. A small portion of them are products of primary metabolism biosynthesized by relatively conserved pathways. Gibberellins (GAs), which are essential for plant growth and development, are diterpenoid phytohormones. (E,E,E)-Geranylgeranyl diphosphate (GGPP) is the precursor for both GAs and other diterpenoids of secondary metabolism. ent-Kaurene biosynthesis from GGPP is a key step of GA formation, which is catalyzed by two sequential and dedicated diterpene synthases (diTPSs): ent-copalyl diphosphate synthase (CPS) and ent-kaurene synthase (KS) of the terpene synthase gene family. Sharing a common evolutionary origin, CPS and KS belong to different TPS subfamilies. Tea plant (Camellia sinensis), the subject of this study, is a leaf-based economic crop. Budbreak mainly manipulated by GAs is a primary factor for targeted tea breeding. The key genes for gibberellin biosynthesis are known; however, they have not yet been characterized in tea plants. Here, we identified and functionally characterized three diterpene biosynthesis-related genes, including one CPS and two highly similar KSs in tea plants. These genes were initially identified through transcriptome sequencing. The functional characterization determined by enzymatic activity assay indicated that CsCPS could catalyze GGPP to form ent-copalyl diphosphate (ent-CPP), which was further used as the substrate by CsKS1 to produce ent-kaurene or by CsKS2 to produce 16α-hydroxy-ent-kaurane with ent-kaurene as a minor product, respectively. We demonstrated that the divergent evolution of diterpene biosynthesis in tea plants resulted from gene duplication of KSs, followed by functional divergence caused by single amino acid variation. This study would provide an insight into the diterpenoid metabolism and GA biosynthesis in tea plants to further understand leaf bud development or insect resistance and to provide a genetic basis for tea plant breeding.

Divergent synthesis of complex diterpenes through a hybrid oxidative approach

Polycyclic diterpenes exhibit many important biological activities, but de novo synthetic access to these molecules is highly challenging because of their structural complexity. Semisynthetic access has also been limited by the lack of chemical tools for scaffold modifications. We report a chemoenzymatic platform to access highly oxidized diterpenes by a hybrid oxidative approach that strategically combines chemical and enzymatic oxidation methods. This approach allows for selective oxidations of previously inaccessible sites on the parent carbocycles and enables abiotic skeletal rearrangements to additional underlying architectures. We synthesized a total of nine complex natural products with rich oxygenation patterns and skeletal diversity in 10 steps or less from ent-steviol.

Transcriptome analysis reveals novel enzymes for apo-carotenoid biosynthesis in saffron and allows construction of a pathway for crocetin synthesis in yeast

Crocus sativus is generally considered the source of saffron spice which is rich in apo-carotenoid compounds such as crocins, crocetin, picrocrocin, and safranal, which possess effective pharmacological activities. However, little is known about the exact genes involved in apo-carotenoid biosynthesis in saffron and the potential mechanism of specific accumulation in the stigma. In this study, we integrated stigmas at different developmental stages to perform in-depth transcriptome and dynamic metabolomic analyses to discover the potential key catalytic steps involved in apo-carotenoid biosynthesis in saffron. A total of 61 202 unigenes were obtained, and 28 regulators and 32 putative carotenogenic genes were captured after the co-expression network analysis. Moreover, 15 candidate genes were predicted to be closely related to safranal and crocin production, in which one aldehyde dehydrogenase (CsALDH3) was validated to oxidize crocetin dialdehyde into crocetin and a crocetin-producing yeast strain was created. In addition, a new branch pathway that catalyses the conversion of geranyl-geranyl pyrophosphate to copalol and ent-kaurene by the class II diterpene synthase CsCPS1 and three class I diterpene synthases CsEKL1/2/3 were investigated for the first time. Such gene to apo-carotenoid landscapes illuminate the synthetic charactersistics and regulators of apo-carotenoid biosynthesis, laying the foundation for a deep understanding of the biosynthesis mechanism and metabolic engineering of apo-carotenoids in plants or microbes.

An alternative splicing alters the product outcome of a class I terpene synthase in Isodon rubescens

The labdane-related diterpenoids are an important superfamily of natural products. Their structural diversity mainly depends on diterpene synthases, which generate the hydrocarbon skeletal structures. Isodon rubescens contains an expanded family of class I terpene synthases with different functions. Here we report a novel class I terpene synthase cDNA (IrKSL3a) with loss of 18 nucleotides compared with the reported cDNA sequence (IrKSL3). Inspection of IrKSL3 genomic sequence indicated that IrKSL3a and IrKSL3 transcripts may be generated by an alternative splicing event that utilizes different 3' splice site. In vitro assays showed that IrKSL3a produced isopimaradiene and miltiradiene, while IrKSL3 only produced miltiradiene. Protein sequence alignment found the six residues encoded by the alternative exon was unique to IrKSL3, which are 17 residues away from the conserved DDXXD motif. A deletion mutant of IrKSL3 showed that maintaining two residues within the six-amino acid is sufficient for miltiradiene production, while the other mutants lost nearly all enzymatic function. Our results illustrated how product outcomes can be changed by alternative splicing, and further gave an interesting example for studying the loop conformation in tuning product outcome in class I terpene synthase.

Diterpene synthases facilitating production of the kaurane skeleton of eriocalyxin B in the medicinal plant Isodon eriocalyx

The Isodon plants (Lamiaceae) have been used in traditional Chinese medicine to alleviate sufferings from inflammations and cancers. This feature has been attributed to the presence of pharmacologically active ent-kaurane diterpenoids such as eriocalyxin B and oridonin. The Isodon eriocalyx (Dunn) Kudô species native to southwest China can accumulate a particularly high content of ent-kaurane diterpenoids (∼1.5% w/w of dried leaves). We previously identified diterpene synthases IeCPS1 and IeCPS2 as ent-copalyl diphosphate synthases (ent-CPS) potentially involved in Isodon ent-kaurane diterpenoids biosynthesis. In this study, analysis of RNA-seq transcriptome of the I. eriocalyx plant revealed three other diterpene synthase genes (IeCPS3, IeKS1, and IeKSL1). Their functional characterization through coupled in vitro enzyme assays has confirmed that IeCPS3 is an ent-CPS specifically producing ent-copalyl diphosphate (ent-CPP). IeKS1 accepted ent-CPP to produce exclusively ent-kaurene and may thus be defined as an ent-kaurene synthase (ent-KS). When IeKSL1 was combined with IeCPS2 or IeCPS3, no product was detected. Based on tissue-specific expression and metabolic localization studies, the IeCPS3 and IeKS1 transcripts were significantly accumulated in leaves where the ent-kaurane diterpenoid eriocalyxin B dominates, whereas weak expression of both were observed in germinating seeds in which gibberellin biosynthetic pathway is normally active. Our findings suggest that both IeCPS3 and IeKS1 possess dual roles in general (gibberellins) and specialized diterpenoid metabolism, such as that of the Isodon ent-kaurane diterpenoids.

A database-driven approach identifies additional diterpene synthase activities in the mint family (Lamiaceae)

Members of the mint family (Lamiaceae) accumulate a wide variety of industrially and medicinally relevant diterpenes. We recently sequenced leaf transcriptomes from 48 phylogenetically diverse Lamiaceae species. Here, we summarize the available chemotaxonomic and enzyme activity data for diterpene synthases (diTPSs) in the Lamiaceae and leverage the new transcriptomes to explore the diTPS sequence and functional space. Candidate genes were selected with an intent to evenly sample the sequence homology space and to focus on species in which diTPS transcripts were found, yet from which no diterpene structures have been previously reported. We functionally characterized nine class II diTPSs and 10 class I diTPSs from 11 distinct plant species and found five class II activities, including two novel activities, as well as a spectrum of class I activities. Among the class II diTPSs, we identified a neo-cleroda-4(18),13E-dienyl diphosphate synthase from Ajuga reptans, catalyzing the likely first step in the biosynthesis of a variety of insect-antifeedant compounds. Among the class I diTPSs was a palustradiene synthase from Origanum majorana, leading to the discovery of specialized diterpenes in that species. Our results provide insights into the diversification of diterpene biosynthesis in the mint family and establish a comprehensive foundation for continued investigation of diterpene biosynthesis in the Lamiaceae.

Genome-wide detection of terpene synthase genes in holy basil (Ocimum sanctum L.)

Holy basil (Ocimum sanctum L.) and sweet basil (Ocimum basilicum L.) are the most commonly grown basil species in India for essential oil production and biosynthesis of potentially volatile and non-volatile phytomolecules with commercial significance. The aroma, flavor and pharmaceutical value of Ocimum species is a significance of its essential oil, which contains most of the monoterpenes and sesquiterpenes. A large number of plants have been studied for characterization and identification of terpene synthase genes, involved in terpenoids biosynthesis. The goal of this study is to discover and identify the putative functional terpene synthase genes in O. sanctum. HMMER search was performed by using a set of 13 well sequenced and annotated plant genomes including the newly sequenced genome of O. sanctum with Pfam-A database locally, using HMMER 3.0 hmmsearch for the two Pfam domains (PF01397 and PF03936). Using this search method 81 putative terpene synthases genes (OsaTPS) were identified in O. sanctum; the study further reveals 47 OsaTPS were putatively functional genes, 19 partial OsaTPS, and 15 OsaTPS as probably pseudogenes. All these identified OsaTPS genes were compared with other plant species, and phylogenetic analysis reveals the subfamily classification of OsaTPS in TPS-a, -b, -c, -e, -f and TPS-g subfamilies clusters. This genome-wide identification of OsaTPS genes, their phylogenetic analysis and secondary metabolite pathway mapping predictions together provide a comprehensive understanding of the TPS gene family in Ocimum sanctum and offer opportunities for the characterization and functional validation of numbers of terpene synthase genes.

Transcriptomic Insight into Terpenoid Biosynthesis and Functional Characterization of Three Diterpene Synthases in Scutellaria barbata

Scutellaria barbata (Lamiaceae) is an important medicinal herb widely used in China, Korea, India, and other Asian countries. Neo-clerodane diterpenoids are the largest known group of Scutellaria diterpenoids and show promising cytotoxic activity against several cancer cell lines. Here, Illumina-based deep transcriptome analysis of flowers, the aerial parts (leaf and stem), and roots of S. barbata was used to explore terpenoid-related genes. In total, 121,958,564 clean RNA-sequence reads were assembled into 88,980 transcripts, with an average length of 1370 nt and N50 length of 2144 nt, indicating high assembly quality. We identified nearly all known terpenoid-related genes (33 genes) involved in biosynthesis of the terpenoid backbone and 14 terpene synthase genes which generate skeletons for different terpenoids. Three full length diterpene synthase genes were functionally identified using an in vitro assay. SbTPS8 and SbTPS9 were identified as normal-CPP and ent-CPP synthase, respectively. SbTPS12 reacts with SbTPS8 to produce miltiradiene. Furthermore, SbTPS12 was proven to be a less promiscuous class I diterpene synthase. These results give a comprehensive understanding of the terpenoid biosynthesis in S. barbata and provide useful information for enhancing the production of bioactive neo-clerodane diterpenoids through genetic engineering.

Comparative transcriptome analysis of roots, stems and leaves of Isodon amethystoides reveals candidate genes involved in Wangzaozins biosynthesis

BACKGROUND

Isodon amethystoides (Ben-th) Cy Wu et Hsuan is an important traditional medicinal plant endowed with pharmacological properties effective in the treatment of various diseases, including pulmonary tuberculosis. The tetracyclic diterpenoids, Wangzaozins (Wangzaozin A, glaucocalyxin A, glaucocalyxin B), are the major bioactive compounds of I. amethystoides. However, the molecular information about the biosynthesis of these compounds still remains unclear.

RESULTS

An examination of the accumulated levels of Wangzaozins in I. amethystoides revealed considerable variations in the root, stem, and leaf tissues of this plant, indicating possible differences in metabolite biosynthesis and accumulation among various tissues. To better elucidate the tetracyclic diterpenoid biosynthesis pathway, we generated transcriptome sequences from the root, stem, and leaf tissues, and performed de novo sequence assembly, yielding 230,974 transcripts and 114,488 unigenes, with average N50 lengths of 1914 and 1241 bp, respectively. Putative functions could be assigned to 73,693 transcripts (31.9%) based on BLAST searches against annotation databases, including GO, KEGG, Swiss-Prot, NR, and Pfam. Moreover, the candidate genes involving in the diterpenoid biosynthesis, such as CPS, KSL, were also analyzed. The expression profiles of eight transcripts, involving the tetracyclic diterpenoid biosynthesis, were validated in different I. amethystoides tissues by qRT-PCR, unraveling the gene expression profile of the pathway. The differential expressions of ISPD, ISPF and ISPH (MEP pathway), and IaCPS and IaKSL (diterpenoid pathway) candidate genes in leaves and roots, may contribute to the high accumulation of Wangzaozins in I. amethystoides leaves.

CONCLUSION

The genomic dataset and analyses reported here lay the foundations for further research on this important medicinal plant.