P450s controlling metabolic bifurcations in plant terpene specialized metabolism

Plant Metabolic Network 16: expansion of underrepresented plant groups and experimentally supported enzyme data

The Plant Metabolic Network (PMN) is a free online database of plant metabolism available at https://plantcyc.org. The latest release, PMN 16, provides metabolic databases representing >1200 metabolic pathways, 1.3 million enzymes, >8000 metabolites, >10 000 reactions and >15 000 citations for 155 plant and green algal genomes, as well as a pan-plant reference database called PlantCyc. This release contains 29 additional genomes compared with PMN 15, including species listed by the African Orphan Crop Consortium and nonflowering plant species. Furthermore, 52 new enzymes with experimentally supported function information have been included in this release. The single-species databases contain a combination of experimental information from the literature and computationally predicted information obtained through PMN's database generation pipeline for a single species, while PlantCyc contains only experimental information but for any species within Viridiplantae. PMN is a comprehensive resource for querying, visualizing, analyzing and interpreting omics data with metabolic knowledge. It also serves as a useful and interactive tool for teaching plant metabolism.

Molecular responses of agroinfiltrated Nicotiana benthamiana leaves expressing suppressor of silencing P19 and influenza virus-like particles

The production of influenza vaccines in plants is achieved through transient expression of viral hemagglutinins (HAs), a process mediated by the bacterial vector Agrobacterium tumefaciens. HA proteins are then produced and matured through the secretory pathway of plant cells, before being trafficked to the plasma membrane where they induce formation of virus-like particles (VLPs). Production of VLPs unavoidably impacts plant cells, as do viral suppressors of RNA silencing (VSRs) that are co-expressed to increase recombinant protein yields. However, little information is available on host molecular responses to foreign protein expression. This work provides a comprehensive overview of molecular changes occurring in Nicotiana benthamiana leaf cells transiently expressing the VSR P19, or co-expressing P19 and an influenza HA. Our data identifies general responses to Agrobacterium-mediated expression of foreign proteins, including shutdown of chloroplast gene expression, activation of oxidative stress responses and reinforcement of the plant cell wall through lignification. Our results also indicate that P19 expression promotes salicylic acid (SA) signalling, a process dampened by co-expression of the HA protein. While reducing P19 level, HA expression also induces specific signatures, with effects on lipid metabolism, lipid distribution within membranes and oxylipin-related signalling. When producing VLPs, dampening of P19 responses thus likely results from lower expression of the VSR, crosstalk between SA and oxylipin pathways, or a combination of both outcomes. Consistent with the upregulation of oxidative stress responses, we finally show that reduction of oxidative stress damage through exogenous application of ascorbic acid improves plant biomass quality during production of VLPs.

Comparative transcriptome and metabolite profiling reveal diverse pattern of CYP-TS gene expression during corosolic acid biosynthesis in Lagerstroemia speciosa (L.) Pers

KEY MESSAGE

Extensive leaf transcriptome profiling and differential gene expression analysis of field grown and elicited shoot cultures of L. speciosa suggest that differential synthesis of CRA is mediated primarily by CYP and TS genes, showing functional diversity. Lagerstroemia speciosa L. is a tree species with medicinal and horticultural attributes. The pentacyclic triterpene, Corosolic acid (CRA) obtained from this species is widely used for the management of diabetes mellitus in traditional medicine. The high mercantile value of the compound and limited availability of innate resources entail exploration of alternative sources for CRA production. Metabolic pathway engineering for enhanced bioproduction of plant secondary metabolites is an attractive proposition for which, candidate genes in the pathway need to be identified and characterized. Therefore, in the present investigation, we focused on the identification of cytochrome P450 (CYP450) and oxidosqualene cyclases (OSC) genes and their differential expression during biosynthesis of CRA. The pattern of differential expression of these genes in the shoot cultures of L. speciosa, elicited with different epigenetic modifiers (azacytidine (AzaC), sodium butyrate (NaBu) and anacardic acid (AA)), was studied in comparison with field grown plant. Further, in vitro cultures with varying (low to high) concentrations of CRA were systematically assessed for the expression of CYP-TS and associated genes involved in CRA biosynthesis by transcriptome sequencing. The sequenced samples were de novo assembled into 180,290 transcripts of which, 92,983 transcripts were further annotated by UniProt. The results are collectively given in co-occurrence heat maps to identify the differentially expressed genes. The combined transcript and metabolite profiles along with RT-qPCR analysis resulted in the identification of CYP-TS genes with high sequence variation. Further, instances of concordant/discordant relation between CRA biosynthesis and CYP-TS gene expression were observed, indicating functional diversity in genes.

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

BACKGROUND

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

MATERIALS AND METHODS

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

RESULTS

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

Combinatorial biosynthesis in yeast leads to over 200 diterpenoids

Diterpenoids form a diverse group of natural products, many of which are or could become pharmaceuticals or industrial chemicals. The modular character of diterpene biosynthesis and the promiscuity of the enzymes involved make combinatorial biosynthesis a promising approach to generate libraries of diverse diterpenoids. Here, we report on the combinatorial assembly in yeast of ten diterpene synthases producing (+)-copalyl diphosphate-derived backbones and four cytochrome P450 oxygenases (CYPs) in diverse combinations. This resulted in the production of over 200 diterpenoids. Based on literature and chemical database searches, 162 of these compounds can be considered new-to-Nature. The CYPs accepted most substrates they were given but remained regioselective with few exceptions. Our results provide the basis for the systematic exploration of the diterpenoid chemical space in yeast using sequence databases.

Integrating Transcriptome and Chemical Analyses to Provide Insights into Biosynthesis of Terpenoids and Flavonoids in the Medicinal Industrial Crop Andrographis paniculate and Its Antiviral Medicinal Parts

Andrographis paniculata is a medicinal plant traditionally used to produce diterpene lactones and flavonoids, which possess various biological activities. Widely distributed in China, India, and other Southeast Asia countries, A. paniculata has become an important economic crop, significantly treating SARS-CoV-2, and is being cultivated on a large scale in southern China. The biosynthesis of active ingredients in A. paniculata are regulated and controlled by genes, but their specific roles are still not fully understood. To further explore the growth regulation factors and utilization of its medicinal parts of this industrial crop, chemical and transcriptome analyses were conducted on the roots, stems, and leaves of A. paniculata to identify the biosynthesis pathways and related candidate genes of the active ingredients. The chemical analysis revealed that the main components of A. paniculata were diterpene lactones and flavonoids, which displayed potential ability to treat SARS-CoV-2 through molecular docking. Moreover, the transcriptome sequencing annotated a total of 40,850 unigenes, including 7962 differentially expressed genes. Among these, 120 genes were involved in diterpene lactone biosynthesis and 60 genes were involved in flavonoid biosynthesis. The expression of diterpene lactone-related genes was the highest in leaves and the lowest in roots, consistent with our content determination results. It is speculated that these highly expressed genes in leaves may be involved in the biosynthesis pathway of diterpenes. Furthermore, two class Ⅰ terpene synthases in A. paniculata transcriptome were also annotated, providing reference for the downstream pathway of the diterpene lactone biosynthesis. With their excellent market value, our experiments will promote the study of the biosynthetic genes for active ingredients in A. paniculata and provide insights for subsequent in vitro biosynthesis.

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.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Identification of P450 Candidates Associated with the Biosynthesis of Physalin-Class Compounds in Physalis angulata

The Physalis genus has long been used as traditional medicine in the treatment of various diseases. Physalins, the characteristic class of compounds in this genus, are major bioactive constituents. To date, the biogenesis of physalins remains largely unknown, except for the recently established knowledge that 24-methyldesmosterol is a precursor of physalin. To identify the genes encoding P450s that are putatively involved in converting 24-methyldesmosterol to physalins, a total of 306 P450-encoding unigenes were retrieved from our recently constructed P. angulata transcriptome. Extensive phylogenetic analysis proposed 21 P450s that might participate in physalin biosynthesis. To validate the candidates, we developed a virus-induced gene silencing (VIGS) system for P. angulata, and four P450 candidates were selected for the VIGS experiments. The reduction in the transcripts of the four P450 candidates by VIGS all led to decreased levels of physalin-class compounds in the P. angulata leaves. Thus, this study provides a number of P450 candidates that are likely associated with the biosynthesis of physalin-class compounds, forming a strong basis to reveal the unknown physalin biosynthetic pathway in the future.

Peroxisome Proliferator Activated Receptor-γ Agonistic Compounds from the Jellyfish-Derived Fungus Cladosporium oxysporum

In our search for peroxisome proliferator-activated receptor (PPAR) agonists, five undescribed compounds, namely two acyclic diterpenes (1 and 2; cladopsol A and cladopsol B), two sesquiterpenes (3 and 4; cladopsol C and cladopsol D), and one C21-ecdysteroid (5; cladopsol E), and 15 known compounds were isolated from the jellyfish-derived fungus - Cladosporium oxysporum. The structures of the undescribed compounds were defined using UV, NMR, HR-ESI-MS, and electronic circular dichroism (ECD) spectroscopy and a modified Mosher's method. Luciferase reporter assay and docking analysis suggested that cladopsol B may function as a PPAR-γ partial agonist with a potential antidiabetic lead which may evade the side effects of full agonists. Moreover, cladopsol B stimulated glucose uptake in HepG2 cells with an efficacy comparable to that of rosiglitazone, but with less side effect induced by lipid accumulation in 3T3-L1 cells. Therefore, cladopsol B could serve as a molecular skeleton in a study of advanced antidiabetic lead with less side effect.

The known, unknown, and the intriguing about members of a critically endangered traditional medicinal plant genus Aconitum

Humanity will always be indebted to plants. In the ongoing scientific era, the 'Herbal Revolution' has helped discover several valuable medicinal plants and associated novel secondary metabolites from the diverse unexplored ecosystems, treating several diseases via phytotherapy. The Aconitum genus comprises several economically-important poisonous mountainous medicinal plant species whose unique biodiversity is on the verge of extinction due to illegal human intervention triggered habitat loss, over-harvesting, and unrestricted trading. Owing to its vast diversity of diterpene alkaloids, most species are extensively used to treat several ailments in rural parts of the world. Irrespective of this, many unexplored and intriguing prospects exist to understand and utilize this critical plant for human benefit. This systematic review tries to fill this gap by compiling information from the sporadically available literature known for ~300 Aconitum spp. regarding its nomenclature and classification, endangerment, plant morphology, ploidy, secondary metabolites, drug pharmacokinetics, conservation, and omics-based computational studies. We also depicted the disparity in the studied model organisms for this diverse genus. The absence of genomic/metagenomic data is becoming a limiting factor in understanding its plant physiology, metabolic pathways, and plant-microbes interactions, and therefore must be promoted. Additionally, government support and public participation are crucial in establishing conservation protocols to save this plant from endangerment.

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.

QTL Detection for Rice Grain Length and Fine Mapping of a Novel Locus qGL6.1

BACKGROUND

Grain length (GL) that is directly associated with appearance quality is a key target of selection in rice breeding. Although abundant quantitative trait locus (QTL) associated with GL have been identified, it was still relatively weak to identify QTL for GL from japonica genetic background, as the shortage of japonica germplasms with long grains. We performed QTLs analysis for GL using a recombinant inbred lines (RILs) population derived from the cross between japonica variety GY8 (short grains) and LX1 (long grains) in four environments.

RESULTS

A total of 197 RILs were genotyped with 285 polymorphic SNP markers. Three QTLs qGL5.3, qGL6.1 and qGL11 were detected to control GL by individual environmental analyses and multi-environment joint analysis. Of these, a major-effect and stable QTL qGL6.1 was identified to be a novel QTL, and its LX1 allele had a positive effect on GL. For fine-mapping qGL6.1, a BC₁F₂ population consisting of 2,487 individuals was developed from a backcross between GY8 and R176, one line with long grain. Eight key informative recombinants were identified by nine kompetitive allele specific PCR (KASP) markers. By analyzing key recombinants, the qGL6.1 locus was narrowed down to a 40.41 kb genomic interval on chromosome 6. One candidate gene LOC_Os06g43304.1 encoding cytochrome P450 (CYP71D55) was finally selected based on the difference in the transcriptional expression and variations in its upstream and downstream region.

CONCLUSIONS

Three QTLs qGL5.3, qGL6.1 and qGL11 were identified to control grain length in rice. One novel QTL qGL6.1 was fine mapped within 40.41 kb region, and LOC_Os06g43304.1 encoding cytochrome P450 (CYP71D55) may be its candidate gene. We propose that the further cloning of the qGL6.1 will facilitate improving appearance quality in japonica varieties.

Plant Specialised Glycosides (PSGs): their biosynthetic enzymatic machinery, physiological functions and commercial potential

Plants, the primary producers of our planet, have evolved from simple aquatic life to very complex terrestrial habitat. This habitat transition coincides with evolution of enormous chemical diversity, collectively termed as 'Plant Specialised Metabolisms (PSMs)', to cope the environmental challenges. Plant glycosylation is an important process of metabolic diversification of PSMs to govern their in planta stability, solubility and inter/intra-cellular transport. Although, individual category of PSMs (terpenoids, phenylpropanoids, flavonoids, saponins, alkaloids, phytohormones, glucosinolates and cyanogenic glycosides) have been well studied; nevertheless, deeper insights of physiological functioning and genomic aspects of plant glycosylation/deglycosylation processes including enzymatic machinery (CYPs, GTs, and GHs) and regulatory elements are still elusive. Therefore, this review discussed the paradigm shift on genomic background of enzymatic machinery, transporters and regulatory mechanism of 'Plant Specialised Glycosides (PSGs)'. Current efforts also update the fundamental understanding about physiological, evolutionary and adaptive role of glycosylation/deglycosylation processes during the metabolic diversification of PSGs. Additionally, futuristic considerations and recommendations for employing integrated next-generation multi-omics (genomics, transcriptomics, proteomics and metabolomics), including gene/genome editing (CRISPR-Cas) approaches are also proposed to explore commercial potential of PSGs.

Genome structure and evolutionary history of frankincense producing Boswellia sacra

Boswellia sacra Flueck (family Burseraceae) tree is wounded to produce frankincense. We report its de novo assembled genome (667.8 Mb) comprising 18,564 high-confidence protein-encoding genes. Comparing conserved single-copy genes across eudicots suggest >97% gene space assembly of B. sacra genome. Evolutionary history shows B. sacra gene-duplications derived from recent paralogous events and retained from ancient hexaploidy shared with other eudicots. The genome indicated a major expansion of Gypsy retroelements in last 2 million years. The B. sacra genetic diversity showed four clades intermixed with a primary genotype—dominating most resin-productive trees. Further, the stem transcriptome revealed that wounding concurrently activates phytohormones signaling, cell wall fortification, and resin terpenoid biosynthesis pathways leading to the synthesis of boswellic acid—a key chemotaxonomic marker of Boswellia. The sequence datasets reported here will serve as a foundation to investigate the genetic determinants of frankincense and other resin-producing species in Burseraceae.

•

Assembly and architecture of frankincense producing Boswellia sacra Flueck

•

Comparative genomics and evolutionary history of frankincense tree within orders

•

Transcriptome of stem part and gene expression patterns of wounding to the tree

•

Resin biosynthesis pathway and related CYP450 enzymes and gene families

Chemotaxonomic investigation of plant terpenoids with an established database (TeroMOL)

Terpenoids constitute the biggest class of plant-derived natural products with diverse chemical structures and extensive biological activities. Interpreting enzyme functions and mining new structures of terpenoids could be inspired by the cheminformatic and chemotaxonomic analysis, whereas it is hampered by the incompleteness of available data for terpenoids. Here a dedicated terpenoids database, TeroMOL, is developed to collect more than 170 000 terpenoids and their derivatives annotated with reported biological sources, along with a user-friendly and freely accessible webserver to visualise and analyse the terpenoids skeletons and organism sources. The quantitative distributions as well as the qualitative trends between terpenoid skeletons and organism sources in plant kingdom are revealed from a chemotaxonomic view, while no comparisons are attempted due to the inherent data biases. Nevertheless, the terpenoid chemomarkers in several organisms are discussed based on the available data with highly enriched and exclusive carbon skeletons. We believe that the TeroMOL database and its accessory computational tools will be very promising for exploring the chemical space and biological sources of terpenoids, and assisting the terpenoid research community in the future.

Bio-inspired chemical space exploration of terpenoids

Many computational methods are devoted to rapidly generating pseudo-natural products to expand the open-ended border of chemical spaces for natural products. However, the accessibility and chemical interpretation were often ignored or underestimated in conventional library/fragment-based or rule-based strategies, thus hampering experimental synthesis. Herein, a bio-inspired strategy (named TeroGen) is developed to mimic the two key biosynthetic stages (cyclization and decoration) of terpenoid natural products, by utilizing physically based simulations and deep learning models, respectively. The precision and efficiency are validated for different categories of terpenoids, and in practice, more than 30 000 sesterterpenoids (10 times as many as the known sesterterpenoids) are predicted to be linked in a reaction network, and their synthetic accessibility and chemical interpretation are estimated by thermodynamics and kinetics. Since it could not only greatly expand the chemical space of terpenoids but also numerate plausible biosynthetic routes, TeroGen is promising for accelerating heterologous biosynthesis, bio-mimic and chemical synthesis of complicated terpenoids and derivatives.

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.

Anti-viral triterpenes: a review

Triterpenes are naturally occurring derivatives biosynthesized following the isoprene rule of Ruzicka. The triterpenes have been reported to possess a wide range of therapeutic applications including anti-viral properties. In this review, the recent studies (2010-2020) concerning the anti-viral activities of triterpenes have been summarized. The structure activity relationship studies have been described as well as brief biosynthesis of these triterpenes is discussed.

Identification and Tissue-Specific Expression Analysis of CYP720B Subfamily Genes in Slash Pine and Loblolly Pine

Diterpene resin acids (DRAs) are major components of pine oleoresin that can effectively resist the invasion of insects and pathogenic microorganisms. The subfamily of cytochrome P450s, CYP720B, catalyzes diterpene products into DRAs. Identifying CYP720B subfamily members and revealing the characteristics of tissue-specific expression would help understand diterpene-rich structures and diverse types. Slash pine and loblolly pine are important pines that provide oleoresin products. In this study, we identified CYP720B candidate genes based on the Pinus taeda V2.0 genome and full-length transcriptome of slash pine by PacBio. A total of 17 genes in slash pine and 19 in loblolly pine were identified and classified into four main clades by phylogenetic analysis. An analysis of cis-acting elements showed that CYP720B genes were closely related to adversity resistance. The gene expression of these candidates in different tissues was quantified by real-time quantitative PCR (RT–qPCR) analysis. Most of the genes showed relatively higher expression levels in roots and stems than in the other tissues, corresponding with the results of DRA component detection by gas chromatography–mass spectrometry (GC–MS), which indicated that stems and roots might be important tissues in oleoresin biosynthesis. These results provide a valuable resource for a better understanding of the biological role of individual CYP720Bs in slash pine and loblolly pine.

Chloranthus genome provides insights into the early diversification of angiosperms

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

Cytochrome P450-catalyzed biosynthesis of furanoditerpenoids in the bioenergy crop switchgrass (Panicum virgatum L.)

Specialized diterpenoid metabolites are important mediators of plant-environment interactions in monocot crops. To understand metabolite functions in plant environmental adaptation that ultimately can enable crop improvement strategies, a deeper knowledge of the underlying species-specific biosynthetic pathways is required. Here, we report the genomics-enabled discovery of five cytochrome P450 monooxygenases (CYP71Z25-CYP71Z29) that form previously unknown furanoditerpenoids in the monocot bioenergy crop Panicum virgatum (switchgrass). Combinatorial pathway reconstruction showed that CYP71Z25-CYP71Z29 catalyze furan ring addition directly to primary diterpene alcohol intermediates derived from distinct class II diterpene synthase products. Transcriptional co-expression patterns and the presence of select diterpenoids in switchgrass roots support the occurrence of P450-derived furanoditerpenoids in planta. Integrating molecular dynamics, structural analysis and targeted mutagenesis identified active site determinants that contribute to the distinct catalytic specificities underlying the broad substrate promiscuity of CYP71Z25-CYP71Z29 for native and non-native diterpenoids.

Taxonomic Insights and Its Type Cyclization Correlation of Volatile Sesquiterpenes in Vitex Species and Potential Source Insecticidal Compounds: A Review

Sesquiterpenes (SS) are secondary metabolites formed by the bonding of 3 isoprene (C5) units. They play an important role in the defense and signaling of plants to adapt to the environment, face stress, and communicate with the outside world, and their evolutionary history is closely related to their physiological functions. This review considers their presence and extensively summarizes the 156 sesquiterpenes identified in Vitextaxa, emphasizing those with higher concentrations and frequency among species and correlating with the insecticidal activities and defensive responses reported in the literature. In addition, we classify the SS based on their chemical structures and addresses cyclization in biosynthetic origin. Most relevant sesquiterpenes of the Vitex genus are derived from the germacredienyl cation mainly via bicyclogermacrene and germacrene C, giving rise to aromadrendanes, a skeleton with the highest number of representative compounds in this genus, and 6,9-guaiadiene, respectively, indicating the production of 1.10-cyclizing sesquiterpene synthases. These enzymes can play an important role in the chemosystematics of the genus from their corresponding routes and cyclizations, constituting a new approach to chemotaxonomy. In conclusion, this review is a compilation of detailed information on the profile of sesquiterpene in the Vitex genus and, thus, points to new unexplored horizons for future research.

Differences in diterpenoid diversity reveal new evidence for separating the genus Coleus from Plectranthus

Covering: up to 2019The large and medicinally important tropical plant genus Plectranthus (Lamiaceae) was recently split into three separate genera on the basis of molecular and morphological evidence; Plectranthus sensu stricto, Coleus and Equilabium. We found striking differences between the diterpenoids which strongly support this taxonomic split. Coleus is characterised by abietanes oxygenated at C-14 such as royleanones, spirocoleons and acylhydroquinones, which could be useful chemotaxonomic markers to distinguish this genus from Plectranthus s.s. In contrast, the abietanes in Plectranthus s.s. lack C-14 oxygenation, but are often acylated with unusual acids. Equilabium species do not seem to produce diterpenoids. The structures of the nearly 240 abietanes so far reported from Coleus and Plectranthus and their distribution are presented. The aim of this Highlight is to provide an overview of the differences in diterpenoid diversity between these newly defined genera, which are relevant to predict which previously understudied species could hold untapped potential for their medicinal and other economic uses, and to underpin future research on how these plants have evolved to synthesise distinct abietane types.

Volatilome and Essential Oil of Ulomoides dermestoides: A Broad-Spectrum Medical Insect

Ulomoides dermestoides are used as a broad-spectrum medical insect in the alternative treatment of various diseases. Preliminary volatilome studies carried out to date have shown, as the main components, methyl-1,4-benzoquinone, ethyl-1,4-benzoquinone, 1-tridecene, 1-pentadecene, and limonene. This work focused on the production of metabolites and their metabolic variations in U. dermestoides under stress conditions to provide additional valuable information to help better understand the broad-spectrum medical uses. To this end, VOCs were characterized by HS-SPME with PEG and CAR/PDMS fibers, and the first reported insect essential oils were obtained. In HS-SMPE, we found 17 terpenes, six quinones, five alkenes, and four aromatic compounds; in the essential oils, 53 terpenes, 54 carboxylic acids and derivatives, three alkynes, 12 alkenes (1-Pentadecene, EOT1: 77.6% and EOT2: 57.9%), 28 alkanes, nine alkyl disulfides, three aromatic compounds, 19 alcohols, three quinones, and 12 aldehydes were identified. Between both study approaches, a total of 171 secondary metabolites were identified with no previous report for U. dermestoides. A considerable number of the identified metabolites showed previous studies of the activity of pharmacological interest. Therefore, considering the wide variety of activities reported for these metabolites, this work allows a broader vision of the therapeutic potential of U. dermestoides in traditional medicine.

The scarecrow-like transcription factor SlSCL3 regulates volatile terpene biosynthesis and glandular trichome size in tomato (Solanum lycopersicum)

Tomato (Solanum lycopersicum L.) type VI glandular trichomes that occur on the surface of leaves, stems, young fruits and flowers produce and store a blend of volatile monoterpenes and sesquiterpenes. These compounds play important roles in the interaction with pathogens and herbivorous insects. Although the function of terpene synthases in the biosynthesis of volatile terpenes in tomato has been comprehensively investigated, the deciphering of their transcriptional regulation is only just emerging. We selected transcription factors that are over-expressed in trichomes based on existing transcriptome data and silenced them individually by virus-induced gene silencing. Of these, SlSCL3, a scarecrow-like (SCL) subfamily transcription factor, led to a significant decrease in volatile terpene content and expression of the corresponding terpene synthase genes when its transcription level was downregulated. Overexpression of SlSCL3 dramatically increased both the volatile terpene content and glandular trichome size, whereas its homozygous mutants showed reduced terpene biosynthesis. However, its heterozygous mutants also showed a significantly elevated volatile terpene content and enlarged glandular trichomes, similar to the overexpression plants. SlSCL3 modulates the expression of terpene biosynthetic pathway genes by transcriptional activation, but neither direct protein-DNA binding nor interaction with known regulators was observed. Moreover, transcript levels of the endogenous copy of SlSCL3 were decreased in the overexpression plants but increased in the heterozygous and homozygous mutants, suggesting feedback repression of its own promoter. Taken together, our results provide new insights into the role of SlSCL3 in the complex regulation of volatile terpene biosynthesis and glandular trichome development in tomato.

Metabolomic and phenotypic implications of the application of fertilization products containing microcontaminants in lettuce (Lactuca sativa)

Cultivation practice using organic amendments is plausible to ensure global food security. However, plant abiotic stress due to the presence of metals and organic microcontaminants (OMCs) in fertilization products cannot be overlooked. In this study, we monitored lettuce metabolism and phenotypic response following the application of either sewage sludge (SS), the organic fraction of municipal solid waste, swine manure (SM), chemical fertilizers (CF), or no amendment (C) in a greenhouse facility. The experimental set-up consisted of five treatments with five replicates (25 experimental units randomly distributed). All fertilizers were supplied at the equivalent agronomic total nitrogen dose, but the occurrence of trace metals and/or OMCs was greater in the SS and SM than the rest. Non-target metabolomic analysis (high-resolution mass spectrometry coupled with partial least squares regression) identified more than 300 plant metabolites (amino acids, organic acids, sugar alcohols, and sugars), 55 of which showed significant changes in their relative abundances depending on the type of amendment. Functional analysis indicated that the use of CF or SS increased the levels of metabolites involved in carbohydrate and nitrogen metabolism. Therefore, although SS and SM fertilizers had a greater presence of heavy metals and/or OMCs, our results indicate that they did not induce measurable adverse effects in the lettuce phenotype or metabolism. Metabolic changes between fertilizers (CF and SS vs. C and SM) were mainly due to nitrogen availability.

Integrated proteomic and metabolomic analyses reveal the importance of aroma precursor accumulation and storage in methyl jasmonate-primed tea leaves

In response to preharvest priming with exogenous methyl jasmonate (MeJA), tea plants adjust their physiological behavior at the molecular level. The whole-organism reconfiguration of aroma formation from the precursor to storage is poorly understood. In this study, we performed iTRAQ proteomic analysis and identified 337, 246, and 413 differentially expressed proteins in tea leaves primed with MeJA for 12 h, 24 h, and 48 h, respectively. Furthermore, a total of 266 nonvolatile and 100 volatile differential metabolites were identified by utilizing MS-based metabolomics. A novel approach that incorporated the integration of extended self-organizing map-based dimensionality was applied. The vivid time-scale changes tracing physiological responses in MeJA-primed tea leaves are marked in these maps. Jasmonates responded quickly to the activation of the jasmonic acid pathway in tea leaves, while hydroxyl and glycosyl jasmonates were biosynthesized simultaneously on a massive scale to compensate for the exhausted defense. The levels of α-linolenic acid, geranyl diphosphate, farnesyl diphosphate, geranylgeranyl diphosphate, and phenylalanine, which are crucial aroma precursors, were found to be significantly changed in MeJA-primed tea leaves. Green leaf volatiles, volatile terpenoids, and volatile phenylpropanoids/benzenoids were spontaneously biosynthesized from responding precursors and subsequently converted to their corresponding glycosidic forms, which can be stably stored in tea leaves. This study elucidated the physiological response of tea leaves primed with exogenous methyl jasmonate and revealed the molecular basis of source and sink changes on tea aroma biosynthesis and catabolism in response to exogenous stimuli. The results significantly enhance our comprehensive understanding of tea plant responses to exogenous treatment and will lead to the development of promising biotechnologies to improve fresh tea leaf quality.

Full-Length Transcriptome Analyses of Genes Involved in Triterpenoid Saponin Biosynthesis of Psammosilene tunicoides Hairy Root Cultures With Exogenous Salicylic Acid

Triterpenoid saponins constitute a diverse class of bioactive compounds in medicinal plants. Salicylic acid (SA) is an efficient elicitor for secondary metabolite production, but a transcriptome-wide regulatory network of SA-promoted triterpenoid saponin biosynthesis remains little understood. In the current study, we described the establishment of the hairy root culture system for Psammosilene tunicoides, a triterpenoid saponin-producing medicinal herb in China, using genetic transformation by Agrobacterium rhizogenes. Compared to controls, we found that total saponin content was dramatically increased (up to 2.49-fold) by the addition of 5 mg/L SA in hairy roots for 1 day. A combination of single-molecule real-time (SMRT) and next-generation sequencing (Illumina RNA-seq) was generated to analyze the full-length transcriptome data for P. tunicoides, as well as the transcript profiles in treated (8 and 24 h) and non-treated (0 h) groups with 5 mg/L SA in hairy roots. A total of 430,117 circular consensus sequence (CCS) reads, 16,375 unigenes and 4,678 long non-coding RNAs (lncRNAs) were obtained. The average length of unigenes (2,776 bp) was much higher in full-length transcriptome than that derived from single RNA-seq (1,457 bp). The differentially expressed genes (DEGs) were mainly enriched in the metabolic process. SA up-regulated the unigenes encoding SA-binding proteins and antioxidant enzymes in comparison with controls. Additionally, we identified 89 full-length transcripts encoding enzymes putatively involved in saponin biosynthesis. The candidate transcription factors (WRKY, NAC) and structural genes (AACT, DXS, SE, CYP72A) might be the key regulators in SA-elicited saponin accumulation. Their expression was further validated by quantitative real-time PCR (qRT-PCR). These findings preliminarily elucidate the regulatory mechanisms of SA on triterpenoid saponin biosynthesis in the transcriptomic level, laying a foundation for SA-elicited saponin augmentation in P. tunicoides.

Stereocontrolled Radical Bicyclizations of Oxygenated Precursors Enable Short Syntheses of Oxidized Abietane Diterpenoids

The power of cation-initiated cyclizations of polyenes for the synthesis of polycyclic terpenoids cannot be overstated. However, a major limitation is the intolerance of many relevant reaction conditions toward the inclusion in the substrate of polar functionality, particularly in unprotected form. Radical polycyclizations are important alternatives to bioinspired cationic variants, in part owing to the range of possible initiation strategies, and in part for the functional group tolerance of radical reactions. In this article, we demonstrate that Co-catalyzed MHAT-initiated radical bicyclizations are not only tolerant of oxidation at virtually every position in the substrate, oftentimes in unprotected form, but these functional groups can also contribute to high levels of stereochemical control in these complexity-generating transformations. Specifically, we show the effects of protected or unprotected hydroxy groups at six different positions and their impact on stereoselectivity. Further, we show how multiply oxidized substrates perform in these reactions, and finally, we document the utility of these reactions in the synthesis of three aromatic abietane diterpenoids.

De novo RNA sequencing and analysis reveal the putative genes involved in diterpenoid biosynthesis in Aconitum vilmorinianum roots

In this study, the putative genes involved in diterpenoid alkaloids biosynthesis in A. vilmorinianum roots were revealed by transcriptome sequencing. 59.39 GB of clean bases and 119,660 unigenes were assembled, of which 69,978 unigenes (58.48%) were annotated. We identified 27 classes of genes (139 candidate genes) involved in the synthesis of diterpenoid alkaloids, including the mevalonate (MVA) pathway, the methylerythritol 4-phosphate (MEP) pathway, the farnesyl diphosphate regulatory pathway, and the diterpenoid scaffold synthetic pathway. 12 CYP450 genes were identified. We found that hydroxymethylglutaryl-CoA reductase was the key enzyme in MVA metabolism, which was regulated by miR6300. Transcription factors, such as bHLH, AP2/EREBP, and MYB, used to synthesize the diterpenes were analyzed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02646-6.

Genetic elucidation of interconnected antibiotic pathways mediating maize innate immunity

Specialized metabolites constitute key layers of immunity that underlie disease resistance in crops; however, challenges in resolving pathways limit our understanding of the functions and applications of these metabolites. In maize (Zea mays), the inducible accumulation of acidic terpenoids is increasingly considered to be a defence mechanism that contributes to disease resistance. Here, to understand maize antibiotic biosynthesis, we integrated association mapping, pan-genome multi-omic correlations, enzyme structure-function studies and targeted mutagenesis. We define ten genes in three zealexin (Zx) gene clusters that encode four sesquiterpene synthases and six cytochrome P450 proteins that collectively drive the production of diverse antibiotic cocktails. Quadruple mutants in which the ability to produce zealexins (ZXs) is blocked exhibit a broad-spectrum loss of disease resistance. Genetic redundancies ensuring pathway resiliency to single null mutations are combined with enzyme substrate promiscuity, creating a biosynthetic hourglass pathway that uses diverse substrates and in vivo combinatorial chemistry to yield complex antibiotic blends. The elucidated genetic basis of biochemical phenotypes that underlie disease resistance demonstrates a predominant maize defence pathway and informs innovative strategies for transferring chemical immunity between crops.

The foxtail millet (Setaria italica) terpene synthase gene family

Terpenoid metabolism plays vital roles in stress defense and the environmental adaptation of monocot crops. Here, we describe the identification of the terpene synthase (TPS) gene family of the panicoid food and bioenergy model crop foxtail millet (Setaria italica). The diploid S. italica genome contains 32 TPS genes, 17 of which were biochemically characterized in this study. Unlike other thus far investigated grasses, S. italica contains TPSs producing all three ent-, (+)- and syn-copalyl pyrophosphate stereoisomers that naturally occur as central building blocks in the biosynthesis of distinct monocot diterpenoids. Conversion of these intermediates by the promiscuous TPS SiTPS8 yielded different diterpenoid scaffolds. Additionally, a cytochrome P450 monooxygenase (CYP99A17), which genomically clustered with SiTPS8, catalyzes the C19 hydroxylation of SiTPS8 products to generate the corresponding diterpene alcohols. The presence of syntenic orthologs to about 19% of the S. italica TPSs in related grasses supports a common ancestry of selected pathway branches. Among the identified enzyme products, abietadien-19-ol, syn-pimara-7,15-dien-19-ol and germacrene-d-4-ol were detectable in planta, and gene expression analysis of the biosynthetic TPSs showed distinct and, albeit moderately, inducible expression patterns in response to biotic and abiotic stress. In vitro growth-inhibiting activity of abietadien-19-ol and syn-pimara-7,15-dien-19-ol against Fusarium verticillioides and Fusarium subglutinans may indicate pathogen defensive functions, whereas the low antifungal efficacy of tested sesquiterpenoids supports other bioactivities. Together, these findings expand the known chemical space of monocot terpenoid metabolism to enable further investigations of terpenoid-mediated stress resilience in these agriculturally important species.

Genome assembly of Chiococca alba uncovers key enzymes involved in the biosynthesis of unusual terpenoids

Chiococca alba (L.) Hitchc. (snowberry), a member of the Rubiaceae, has been used as a folk remedy for a range of health issues including inflammation and rheumatism and produces a wealth of specialized metabolites including terpenes, alkaloids, and flavonoids. We generated a 558 Mb draft genome assembly for snowberry which encodes 28,707 high-confidence genes. Comparative analyses with other angiosperm genomes revealed enrichment in snowberry of lineage-specific genes involved in specialized metabolism. Synteny between snowberry and Coffea canephora Pierre ex A. Froehner (coffee) was evident, including the chromosomal region encoding caffeine biosynthesis in coffee, albeit syntelogs of N-methyltransferase were absent in snowberry. A total of 27 putative terpene synthase genes were identified, including 10 that encode diterpene synthases. Functional validation of a subset of putative terpene synthases revealed that combinations of diterpene synthases yielded access to products of both general and specialized metabolism. Specifically, we identified plausible intermediates in the biosynthesis of merilactone and ribenone, structurally unique antimicrobial diterpene natural products. Access to the C. alba genome will enable additional characterization of biosynthetic pathways responsible for health-promoting compounds in this medicinal species.

Genomics-enabled analysis of specialized metabolism in bioenergy crops: current progress and challenges

Plants produce a staggering diversity of specialized small molecule metabolites that play vital roles in mediating environmental interactions and stress adaptation. This chemical diversity derives from dynamic biosynthetic pathway networks that are often species-specific and operate under tight spatiotemporal and environmental control. A growing divide between demand and environmental challenges in food and bioenergy crop production has intensified research on these complex metabolite networks and their contribution to crop fitness. High-throughput omics technologies provide access to ever-increasing data resources for investigating plant metabolism. However, the efficiency of using such system-wide data to decode the gene and enzyme functions controlling specialized metabolism has remained limited; due largely to the recalcitrance of many plants to genetic approaches and the lack of 'user-friendly' biochemical tools for studying the diverse enzyme classes involved in specialized metabolism. With emphasis on terpenoid metabolism in the bioenergy crop switchgrass as an example, this review aims to illustrate current advances and challenges in the application of DNA synthesis and synthetic biology tools for accelerating the functional discovery of genes, enzymes and pathways in plant specialized metabolism. These technologies have accelerated knowledge development on the biosynthesis and physiological roles of diverse metabolite networks across many ecologically and economically important plant species and can provide resources for application to precision breeding and natural product metabolic engineering.

Grape terpenoids: flavor importance, genetic regulation, and future potential

Terpenes significantly affect the flavor and quality of grapes and wine. This review summarizes recent research on terpenoids with regard to grape wine. Although, the grapevine terpene synthase gene family is the largest identified, genetic modifications involving terpenes to improve wine flavor have received little attention. Key enzyme modulation alters metabolite production. Over the last decade, the heterologous manipulation of grape glycosidase has been used to alter terpenoids, and cytochrome P450s may affect terpene synthesis. Metabolic and genetic engineering can further modify terpenoid metabolism, while using transgenic grapevines (trait transfer to the plant) could yield more flavorful wine. We also discuss traits involved in wine aroma quality, and the strategies that can be used to improve grapevine breeding technology.

Strategies for the production of biochemicals in bioenergy crops

Industrial crops are grown to produce goods for manufacturing. Rather than food and feed, they supply raw materials for making biofuels, pharmaceuticals, and specialty chemicals, as well as feedstocks for fabricating fiber, biopolymer, and construction materials. Therefore, such crops offer the potential to reduce our dependency on petrochemicals that currently serve as building blocks for manufacturing the majority of our industrial and consumer products. In this review, we are providing examples of metabolites synthesized in plants that can be used as bio-based platform chemicals for partial replacement of their petroleum-derived counterparts. Plant metabolic engineering approaches aiming at increasing the content of these metabolites in biomass are presented. In particular, we emphasize on recent advances in the manipulation of the shikimate and isoprenoid biosynthetic pathways, both of which being the source of multiple valuable compounds. Implementing and optimizing engineered metabolic pathways for accumulation of coproducts in bioenergy crops may represent a valuable option for enhancing the commercial value of biomass and attaining sustainable lignocellulosic biorefineries.

Specialized diterpenoid metabolism in monocot crops: Biosynthesis and chemical diversity

Among the myriad specialized metabolites that plants employ to mediate interactions with their environment, diterpenoids form a chemically diverse group with vital biological functions. A few broadly abundant diterpenoids serve as core pathway intermediates in plant general metabolism. The majority of plant diterpenoids, however, function in specialized metabolism as often species-specific chemical defenses against herbivores and microbial diseases, in below-ground allelopathic interactions, as well as abiotic stress responses. Dynamic networks of anti-microbial diterpenoids were first demonstrated in rice (Oryza sativa) over four decades ago, and more recently, unique diterpenoid blends with demonstrated antibiotic bioactivities were also discovered in maize (Zea mays). Enabled by advances in -omics and biochemical approaches, species-specific diterpenoid-diversifying enzymes have been identified in these and other Poaceous species, including wheat (Triticum aestivum) and switchgrass (Panicum virgatum), and are discussed in this article with an emphasis on the critical diterpene synthase and cytochrome P450 monooxygenase families and their products. The continued investigation of the biosynthesis, diversity, and function of terpenoid-mediated crop defenses provides foundational knowledge to enable the development of strategies for improving crop resistance traits in the face of impeding pest, pathogen, and climate pressures impacting global agricultural production.

Predicting Rearrangement-Competent Terpenoid Oxidation Levels

Results of density functional theory calculations on rearrangements of potential biosynthetic precursors to the sesquiterpenoid illisimonin A reveal that only some possible precursors, those with certain specific oxidation patterns, are rearrangement-competent.

Phytochemical and comparative transcriptome analyses reveal different regulatory mechanisms in the terpenoid biosynthesis pathways between Matricaria recutita L. and Chamaemelum nobile L

Background

Matricaria recutita (German chamomile) and Chamaemelum nobile (Roman chamomile) belong to the botanical family Asteraceae. These two herbs are not only morphologically distinguishable, but their secondary metabolites – especially the essential oils present in flowers are also different, especially the terpenoids. The aim of this project was to preliminarily identify regulatory mechanisms in the terpenoid biosynthetic pathways that differ between German and Roman chamomile by performing comparative transcriptomic and metabolomic analyses.

Results

We determined the content of essential oils in disk florets and ray florets in these two chamomile species, and found that the terpenoid content in flowers of German chamomile is greater than that of Roman chamomile. In addition, a comparative RNA-seq analysis of German and Roman chamomile showed that 54% of genes shared > 75% sequence identity between the two species. In particular, more highly expressed DEGs (differentially expressed genes) and TF (transcription factor) genes, different regulation of CYPs (cytochrome P450 enzymes), and rapid evolution of downstream genes in the terpenoid biosynthetic pathway of German chamomile could be the main reasons to explain the differences in the types and levels of terpenoid compounds in these two species. In addition, a phylogenetic tree constructed from single copy genes showed that German chamomile and Roman chamomile are closely related to Chrysanthemum nankingense.

Conclusion

This work provides the first insights into terpenoid biosynthesis in two species of chamomile. The candidate unigenes related to terpenoid biosynthesis will be important in molecular breeding approaches to modulate the essential oil composition of Matricaria recutita and Chamaemelum nobile.

A New Epi-neoverrucosane-type Diterpenoid from the Liverwort Pleurozia subinflata in Borneo

New bioactive 13-epi-neoverrucosane diterpenoid, 5β-acetoxy-13-epi-neoverrucosanic acid (1) along with three known secondary metabolites, 13-epi-neoverrucosan-5β-ol (2), chelodane (3) and (E)-β-farnesene (4) were isolated from the MeOH extract of east Malaysia's liverwort Pleurozia subinflata. The chemical structure of new compound was elucidated by the analyses of its spectroscopic data (FTIR, NMR and HR-ESI-MS). These epi-neoverrucosane-type compounds seem to be notable chemosystematic markers for P. subinflata in Borneo. Compound 3 was widespread in marine sponges however this is the first record for 3 to be found in liverwort. These metabolites were tested for their antifungal potentials against selected fungi from the marine environment. Compound 1 exhibited effective antifungal activity against Lagenidium thermophilum.

De novo transcriptome of Gymnema sylvestre identified putative lncRNA and genes regulating terpenoid biosynthesis pathway

Gymnema sylvestre is a highly valuable medicinal plant in traditional Indian system of medicine and used in many polyherbal formulations especially in treating diabetes. However, the lack of genomic resources has impeded its research at molecular level. The present study investigated functional gene profile of G. sylvestre via RNA sequencing technology. The de novo assembly of 88.9 million high quality reads yielded 23,126 unigenes, of which 18116 were annotated against databases such as NCBI nr database, gene ontology (GO), KEGG, Pfam, CDD, PlantTFcat, UniProt & GreeNC. Total 808 unigenes mapped to 78 different Transcription Factor families, whereas 39 unigenes assigned to CYP450 and 111 unigenes coding for enzymes involved in the biosynthesis of terpenoids including transcripts for synthesis of important compounds like Vitamin E, beta-amyrin and squalene. Among them, presence of six important enzyme coding transcripts were validated using qRT-PCR, which showed high expression of enzymes involved in methyl-erythritol phosphate (MEP) pathway. This study also revealed 1428 simple sequence repeats (SSRs), which may aid in molecular breeding studies. Besides this, 8 putative long non-coding RNAs (lncRNAs) were predicted from un-annotated sequences, which may hold key role in regulation of essential biological processes in G. sylvestre. The study provides an opportunity for future functional genomic studies and to uncover functions of the lncRNAs in G. sylvestre.

An integrated strategy to identify genes responsible for sesquiterpene biosynthesis in turmeric

Metabolic module, gene expression pattern and PLS modeling were integrated to precisely identify the terpene synthase responsible for sesquiterpene formation. Functional characterization confirmed the feasibility and sensitivity of this strategy. Plant secondary metabolite biosynthetic pathway elucidation is crucial for the production of these compounds with metabolic engineering. In this study, an integrated strategy was employed to predict the gene function of sesquiterpene synthase (STS) genes using turmeric as a model. Parallel analysis of gene expression patterns and metabolite modules narrowed the candidates into an STS group in which the STSs showed a similar expression pattern. The projections to latent structures by means of partial least squares model was further employed to establish a clear relationship between the candidate STS genes and metabolites and to predict three STSs (ClTPS16, ClTPS15 and ClTPS14) involved in the biosynthesis of several sesquiterpene skeletons. Functional characterization revealed that zingiberene and β-sesquiphellandrene were the major products of ClTPS16, and β-eudesmol was produced by ClTPS15, both of which indicated the accuracy of the prediction. Functional characterization of a control STS, ClTPS1, produced a small amount of β-sesquiphellandrene, as predicted, which confirmed the sensitivity of metabolite module analysis. This integrated strategy provides a methodology for gene function predictions, which represents a substantial improvement in the elucidation of biosynthetic pathways in nonmodel plants.

Cytochrome P450 enzymes: A driving force of plant diterpene diversity

In plant terpene biosynthesis, oxidation of the hydrocarbon backbone produced by terpene synthases is typically carried out by cytochrome P450 oxygenases (CYPs). The modifications introduced by CYPs include hydroxylations, sequential oxidations at one position and ring rearrangements and closures. These reactions significantly expand the structural diversity of terpenoids, but also provide anchoring points for further decorations by various transferases. In recent years, there has been a significant increase in reports of CYPs involved in plant terpene pathways. Plant diterpenes represent an important class of metabolites that includes hormones and a number of industrially relevant compounds such as pharmaceutical, aroma or food ingredients. In this review, we provide a comprehensive survey on CYPs reported to be involved in plant diterpene biosynthesis to date. A phylogenetic analysis showed that only few CYP clans are represented in diterpene biosynthesis, namely CYP71, CYP85 and CYP72. Remarkably few CYP families and subfamilies within those clans are involved, indicating specific expansion of these clades in plant diterpene biosynthesis. Nonetheless, the evolutionary trajectory of CYPs of specialized diterpene biosynthesis is diverse. Some are recently derived from gibberellin biosynthesis, while others have a more ancient history with recent expansions in specific plant families. Among diterpenoids, labdane-related diterpenoids represent a dominant class. The availability of CYPs from diverse plant species able to catalyze oxidations in specific regions of the labdane-related backbones provides opportunities for combinatorial biosynthesis to produce novel diterpene compounds that can be screened for biological activities of interest.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The genome of the medicinal plant Andrographis paniculata provides insight into the biosynthesis of the bioactive diterpenoid neoandrographolide

Andrographis paniculata is a herbaceous dicot plant widely used for its anti-inflammatory and anti-viral properties across its distribution in China, India and other Southeast Asian countries. A. paniculata was used as a crucial therapeutic treatment during the influenza epidemic of 1919 in India, and is still used for the treatment of infectious disease in China. A. paniculata produces large quantities of the anti-inflammatory diterpenoid lactones andrographolide and neoandrographolide, and their analogs, which are touted to be the next generation of natural anti-inflammatory medicines for lung diseases, hepatitis, neurodegenerative disorders, autoimmune disorders and inflammatory skin diseases. Here, we report a chromosome-scale A. paniculata genome sequence of 269 Mb that was assembled by Illumina short reads, PacBio long reads and high-confidence (Hi-C) data. Gene annotation predicted 25 428 protein-coding genes. In order to decipher the genetic underpinning of diterpenoid biosynthesis, transcriptome data from seedlings elicited with methyl jasmonate were also obtained, which enabled the identification of genes encoding diterpenoid synthases, cytochrome P450 monooxygenases, 2-oxoglutarate-dependent dioxygenases and UDP-dependent glycosyltransferases potentially involved in diterpenoid lactone biosynthesis. We further carried out functional characterization of pairs of class-I and -II diterpene synthases, revealing the ability to produce diversified labdane-related diterpene scaffolds. In addition, a glycosyltransferase able to catalyze O-linked glucosylation of andrograpanin, yielding the major active product neoandrographolide, was also identified. Thus, our results demonstrate the utility of the combined genomic and transcriptomic data set generated here for the investigation of the production of the bioactive diterpenoid lactone constituents of the important medicinal herb A. paniculata.

Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses

Loblolly pine (Pinus taeda L.) is one of the most important species for oleoresin (a mixture of terpenoids) in South China. The high oleoresin content of loblolly pine is associated with resistance to bark beetles and other economic benefits. In this study, we conducted transcriptome analyses of loblolly pine secondary xylem to gain insight into the genes involved in terpenoid biosynthesis. A total of 372 unigenes were identified as being critical for oleoresin production, including genes for ATP-binding cassette (ABC) transporters, the cytochrome P450 (CYP) protein family, and terpenoid backbone biosynthesis enzymes. Six key genes involved in terpenoid biosynthetic pathways were selected for multiple sequence alignment, conserved motif prediction, and phylogenetic and expression profile analyses. The protein sequences of all six genes exhibited a higher degree of sequence conservation, and upstream genes were relatively more conserved than downstream genes in terpenoid biosynthetic pathways. The N-terminal regions of these sequences were less conserved than the C-terminal ends, as the N-terminals were quite diverse in both length and composition. The phylogenetic analyses revealed that most genes originated from gene duplication after species divergence, and partial genes exhibited incomplete lineage sorting. In addition, the expression profile analyses showed that all six genes exhibited high expression levels during the high-oleoresin-yielding phase.

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.