Combined metabolome and transcriptome profiling provides new insights into diterpene biosynthesis in S. pomifera glandular trichomes

Chemical and transcriptomic analyses of leaf trichomes from Cistus creticus subsp. creticus reveal the biosynthetic pathways of certain labdane-type diterpenoids and their acetylated forms

Labdane-related diterpenoids (LRDs), a subgroup of terpenoids, exhibit structural diversity and significant commercial and pharmacological potential. LRDs share the characteristic decalin-labdanic core structure that derives from the cycloisomerization of geranylgeranyl diphosphate (GGPP). Labdanes derive their name from the oleoresin known as 'Labdanum', 'Ladano', or 'Aladano', used since ancient Greek times. Acetylated labdanes, rarely identified in plants, are associated with enhanced biological activities. Chemical analysis of Cistus creticus subsp. creticus revealed labda-7,13(E)-dien-15-yl acetate and labda-7,13(E)-dien-15-ol as major constituents. In addition, novel labdanes such as cis-abienol, neoabienol, ent-copalol, and one as yet unidentified labdane-type diterpenoid were detected for the first time. These compounds exhibit developmental regulation, with higher accumulation observed in young leaves. Using RNA-sequencing (RNA-seq) analysis of young leaf trichomes, it was possible to identify, clone, and eventually functionally characterize labdane-type diterpenoid synthase (diTPS) genes, encoding proteins responsible for the production of labda-7,13(E)-dien-15-yl diphosphate (endo-7,13-CPP), labda-7,13(E)-dien-15-yl acetate, and labda-13(E)-ene-8α-ol-15-yl acetate. Moreover, the reconstitution of labda-7,13(E)-dien-15-yl acetate and labda-13(E)-ene-8α-ol-15-yl acetate production in yeast is presented. Finally, the accumulation of LRDs in different plant tissues showed a correlation with the expression profiles of the corresponding genes.

Multi-omics profiling reveal responses of three major Dendrobium species from different growth years to medicinal components

Dendrobium is a perennial herb found in Asia that is known for its medicinal and ornamental properties. Studies have shown that the stem is the primary medicinal component of Dendrobium spp. To investigate the effect of the species and age of Dendrobium (in years) on the content of its medicinal components, we collected the stems of 1-to-4-year-old D. officinale, D. moniliforme, and D. huoshanense, sequenced the transcriptome, metabolome, and microbiome, and analyzed the data in a comprehensive multi-omics study. We identified 10,426 differentially expressed genes (DEGs) with 644 differentially accumulated metabolites (DAMs) from 12 comparative groups and mapped the flavonoid pathway based on DEGs and DAMs. Transcriptomic and metabolomic data indicated a general trend of the accumulation of flavonoids exhibiting pharmacological effects in the three Dendrobium species. In addition, joint metabolome and microbiome analyses showed that actinobacteria was closely associated with flavonoid synthesis with increasing age. Our findings provide novel insights into the interactions of flavonoids of Dendrobium with the transcriptome and microbiome.

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

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

Phytochemical screening and biological evaluation of Greek sage (Salvia fruticosa Mill.) extracts

This study explores the influence of extraction solvents on the composition and bioactivity of Salvia fruticosa extracts. Ultrasound-assisted extraction with water, ethanol and their mixtures in variable proportions was used to produce four different extracts. An untargeted UPLC/MS‑based metabolomics was performed to discover metabolites profile variation between the extracts. In the analyzed samples, 2704 features had been detected, of which 95 were tentatively identified. The concentrations of the important metabolites, namely, caffeic acid, carnosic acid, carnosol, rosmarinic acid, salvianolic acid B and scutellarin, were determined, using UPLC-PDA methods. Rosmarinic acid was the dominant metabolite and antioxidant in all tested extracts, except the aqueous extract, in which scutellarin was the most abundant compound. The extracts and standards were examined for antioxidant activity and xanthine oxidase (XO) inhibitory activity. The most diverse in terms of chemical composition and rich in antioxidant compounds was 70% ethanolic extract and the strongest antioxidant was caffeic acid. All analyzed extracts showed the ability to inhibit XO activity, but the highest value was recorded for 30% ethanolic extract. Among tested standards, the most potent XO inhibitor was caffeic acid. The results suggest that the leaves of Greek sage are a source of natural XO inhibitors and may be an alternative to drugs produced by chemical synthesis.

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

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

Biosynthesis and the Transcriptional Regulation of Terpenoids in Tea Plants (Camellia sinensis)

Terpenes, especially volatile terpenes, are important components of tea aroma due to their unique scents. They are also widely used in the cosmetic and medical industries. In addition, terpene emission can be induced by herbivory, wounding, light, low temperature, and other stress conditions, leading to plant defense responses and plant-plant interactions. The transcriptional levels of important core genes (including HMGR, DXS, and TPS) involved in terpenoid biosynthesis are up- or downregulated by the MYB, MYC, NAC, ERF, WRKY, and bHLH transcription factors. These regulators can bind to corresponding cis-elements in the promoter regions of the corresponding genes, and some of them interact with other transcription factors to form a complex. Recently, several key terpene synthesis genes and important transcription factors involved in terpene biosynthesis have been isolated and functionally identified from tea plants. In this work, we focus on the research progress on the transcriptional regulation of terpenes in tea plants (Camellia sinensis) and thoroughly detail the biosynthesis of terpene compounds, the terpene biosynthesis-related genes, the transcription factors involved in terpene biosynthesis, and their importance. Furthermore, we review the potential strategies used in studying the specific transcriptional regulation functions of candidate transcription factors that have been discriminated to date.

GC-MS and LC-DAD-MS Phytochemical Profiling for Characterization of Three Native Salvia Taxa from Eastern Mediterranean with Antiglycation Properties

Salvia fruticosa and S. pomifera subsp. calycina are native to Eastern Mediterranean and S. pomifera subsp. pomifera is endemic to Greece. The primary aim of this study was to develop an analytical methodology for metabolomic profiling and to study their efficacy in combating glycation, the major biochemical complication of diabetes. After sequential ultrasound-assisted extraction of 2 g of leaves with petroleum ether and 70% methanol, the volatile metabolites in the petroleum ether extracts were studied with GC-MS (Gas Chromatography-Mass Spectrometry), whereas the polar metabolites in the hydroalcoholic extracts were determined and quantified by UHPLC-DAD-ESI-MS (Ultra-High Performance Liquid Chromatography-Diode Array Detector-Mass Spectrometry). This methodology was applied to five populations belonging to the three native taxa. 1,8-Cineole was the predominant volatile (34.8-39.0%) in S. fruticosa, while S. pomifera had a greater content of α-thujone (19.7-41.0%) and β-thujone (6.0-39.1%). Principal Component Analysis (PCA) analysis of the volatiles could discriminate the different taxa. UHPLC-DAD-ESI-MS demonstrated the presence of 50 compounds, twenty of which were quantified. PCA revealed that not only the taxa but also the populations of S. pomifera subsp. pomifera could be differentiated. All Salvia samples inhibited advanced glycation end-product formation in a bovine serum albumin/2-deoxyribose assay; rosmarinic and carnosic acid shared this activity. This study demonstrates the antiglycation activity of S. fruticosa and S. pomifera extracts for the first time and presents a miniaturized methodology for their metabolomic profiling, which could aid chemotaxonomic studies and serve as a tool for their authentication and quality control.

Engineering a powerful green cell factory for robust photoautotrophic diterpenoid production

Diterpenoids display a large and structurally diverse class of natural compounds mainly found as specialized plant metabolites. Due to their diverse biological functions they represent an essential source for various industrially relevant applications as biopharmaceuticals, nutraceuticals, and fragrances. However, commercial production utilizing their native hosts is inhibited by low abundances, limited cultivability, and challenging extraction, while the precise stereochemistry displays a particular challenge for chemical synthesis. Due to a high carbon flux through their native 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway towards photosynthetically active pigments, green microalgae hold great potential as efficient and sustainable heterologous chassis for sustainable biosynthesis of plant-derived diterpenoids. In this study, innovative synthetic biology and efficient metabolic engineering strategies were systematically combined to re-direct the metabolic flux through the MEP pathway for efficient heterologous diterpenoid synthesis in C. reinhardtii. Engineering of the 1-Deoxy-D-xylulose 5-phosphate synthase (DXS) as the main rate-limiting enzyme of the MEP pathway and overexpression of diterpene synthase fusion proteins increased the production of high-value diterpenoids. Applying fully photoautotrophic high cell density cultivations demonstrate potent and sustainable production of the high-value diterpenoid sclareol up to 656 mg L^-1 with a maximal productivity of 78 mg L^-1 day^-1 in a 2.5 L scale photobioreactor, which is comparable to sclareol titers reached by highly engineered yeast. Consequently, this work represents a breakthrough in establishing a powerful phototrophic green cell factory for the competetive use in industrial biotechnology.

Tanshinones: Leading the way into Lamiaceae labdane-related diterpenoid biosynthesis

Tanshinones are the bioactive diterpenoid constituents of the traditional Chinese medicinal herb Danshen (Salvia miltiorrhiza), and are examples of the phenolic abietanes widely found within the Lamiaceae plant family. Due to the significant interest in these labdane-related diterpenoid natural products, their biosynthesis has been intensively investigated. In addition to providing the basis for metabolic engineering efforts, this work further yielded pioneering insights into labdane-related diterpenoid biosynthesis in the Lamiaceae more broadly. This includes stereochemical foreshadowing of aromatization, with novel protein domain loss in the relevant diterpene synthase, as well as broader phylogenetic conservation of the relevant enzymes. Beyond such summary of more widespread metabolism, formation of the furan ring that characterizes the tanshinones also has been recently elucidated. Nevertheless, the biocatalysts for the pair of demethylations remain unknown, and the intriguing potential connection of these reactions to the further aromatization observed in the tanshinones are speculated upon here.

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.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

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.

Transcriptome analysis of near-isogenic lines for glume hairiness of wheat

Hairiness, which is a phenotypic trait common among land plants, primarily affects the stem, leaf, and floral organs. Plant hairiness is associated with complex functions. For example, glume hairiness in wheat is related to the resistance to biotic and abiotic stresses, and may also influence human health. In the present study, two pairs of near-isogenic lines (NILs) for glume hairiness, which were derived from a cross between a Tibetan semi-wild wheat accession (Triticum aestivum ssp. tibetanum Q1028) and a common wheat cultivar (T. aestivum 'Zhengmai 9023'), underwent a glume transcriptome analysis. We detected 27,935 novel genes, of which 18,027 were annotated. Additionally, 488 and 600 differentially expressed genes (DEGs) were detected in NIL1 and NIL2, respectively, with 37 DEGs detected in both NIL pairs. Moreover, 987 and 1584 single nucleotide polymorphisms (SNPs) were detected in NIL1 and NIL2, respectively, with 39 SNPs detected in both NIL pairs, of which most were located in the Hairy glume (Hg) gene region on chromosome arm 1AS. The annotation of the DEGs with gene ontology terms revealed that genes associated with hairiness in Arabidopsis and rice were similarly enriched. The possible functions of these genes related to glume hairiness were examined. The study results provide useful information for identifying candidate genes and the fine-mapping of Hg in the wheat genome.

Integrative Transcriptomic and Metabolic Analyses Provide Insights into the Role of Trichomes in Tea Plant (Camellia Sinensis)

Trichomes, which develop from epidermal cells, are regarded as one of the key features that are involved in the evaluation of tea quality and tea germplasm resources. The metabolites from trichomes have been well characterized in tea products. However, little is known regarding the metabolites in fresh tea trichomes and the molecular differences in trichomes and tea leaves per se. In this study, we developed a method to collect trichomes from tea plant tender shoots, and their main secondary metabolites, including catechins, caffeine, amino acids, and aroma compounds, were determined. We found that the majority of these compounds were significantly less abundant in trichomes than in tea leaves. RNA-Seq was used to investigate the differences in the molecular regulatory mechanism between trichomes and leaves to gain further insight into the differences in trichomes and tea leaves. In total, 52.96 Gb of clean data were generated, and 6560 differentially expressed genes (DEGs), including 4471 upregulated and 2089 downregulated genes, were identified in the trichomes vs. leaves comparison. Notably, the structural genes of the major metabolite biosynthesis pathways, transcription factors, and other key DEGs were identified and comparatively analyzed between trichomes and leaves, while trichome-specific genes were also identified. Our results provide new insights into the differences between tea trichomes and leaves at the metabolic and transcriptomic levels, and open up new doors to further recognize and re-evaluate the role of trichomes in tea quality formation and tea plant growth and development.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Lipid Transfer Proteins (AaLTP3 and AaLTP4) Are Involved in Sesquiterpene Lactone Secretion from Glandular Trichomes in Artemisia annua

In Artemisia annua plants, glandular trichomes (GTs) are responsible for the biosynthesis and secretion of sesquiterpene lactones including artemisinin/arteannuin B. Nonspecific lipid transfer proteins (LTPs) in plants bind and carry lipid molecules across the cell membrane and are also known as secretary proteins. Interestingly, the transcripts of LTP genes are exceptionally abundant in the GTs of A. annua. In the present study, we isolated two trichome-specific LTP genes (AaLTP3 and AaLTP4) from a Korean ecotype of A. annua. AaLTP3 was expressed abundantly in shoots, whereas AaLTP4 was expressed in flowers. The GUS signal driven by the AaLTP3 or AaLTP4 promoter in transgenic A. annua plants revealed that the AaLTP3 promoter was active on hair-like non-GTs and that the AaLTP4 promoter was active on GTs. Analysis of enhanced cyan fluorescent protein (ECFP) fluorescence fused with the AaLTP3 or AaLTP4 protein in transgenic tobacco revealed that ECFP florescence was very bright on secreted lipids of long GTs. Moreover, the florescence was also bright on the head cells of short trichomes and their secreted granules. Immunoblotting analysis of GT exudates in petioles of A. annua revealed a strong positive signal against the AaLTP4 antibody. Overexpression of AaLTP3 or AaLTP4 in transgenic A. annua plants resulted in enhanced production of sesquiterpene lactones (arteannuin B, artemisinin, dihydroartemisinic acid and artemisinic acid) compared with those of wild type. The present study shows that LTP genes (AaLTP3 or AaLTP4) play important roles in the sequestration and secretion of lipids in GTs of A. annua, which is useful for the enhanced production of sesquiterpene lactones by genetic engineering.

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

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

Methyl jasmonate-induced compositional changes of volatile organic compounds in Polygonum minus leaves

Polygonum minus Huds. is a medicinal aromatic plant rich in terpenes, aldehydes, and phenolic compounds. Methyl jasmonate (MeJA) is a plant signaling molecule commonly applied to elicit stress responses to produce plant secondary metabolites. In this study, the effects of exogenous MeJA treatment on the composition of volatile organic compounds (VOCs) in P. minus leaves were investigated by using a metabolomic approach. Time-course changes in the leaf composition of VOCs on days 1, 3, and 5 after MeJA treatment were analyzed through solid-phase microextraction (SPME) and gas chromatography-mass spectrometry (GC-MS). The VOCs found in MeJA-elicited leaves were similar to those found in mock-treated leaves but varied in quantity at different time points. We focused our analysis on the content and composition of monoterpenes, sesquiterpenes, and green leaf volatiles (GLVs) within the leaf samples. Our results suggest that MeJA enhances the activity of biosynthetic pathways for aldehydes and terpenes in P. minus. Hence, the production of aromatic compounds in this medicinal herb can be increased by MeJA elicitation. Furthermore, the relationship between MeJA elicitation and terpene biosynthesis in P. minus was shown through SPME-GC-MS analysis of VOCs combined with transcriptomic analysis of MeJA-elicited P. minus leaves from our previous study.

Orthogonal monoterpenoid biosynthesis in yeast constructed on an isomeric substrate

Synthetic biology efforts for the production of valuable chemicals are frequently hindered by the structure and regulation of the native metabolic pathways of the chassis. This is particularly evident in the case of monoterpenoid production in Saccharomyces cerevisiae, where the canonical terpene precursor geranyl diphosphate is tightly coupled to the biosynthesis of isoprenoid compounds essential for yeast viability. Here, we establish a synthetic orthogonal monoterpenoid pathway based on an alternative precursor, neryl diphosphate. We identify structural determinants of isomeric substrate selectivity in monoterpene synthases and engineer five different enzymes to accept the alternative substrate with improved efficiency and specificity. We combine the engineered enzymes with dynamic regulation of metabolic flux to harness the potential of the orthogonal substrate and improve the production of industrially-relevant monoterpenes by several-fold compared to the canonical pathway. This approach highlights the introduction of synthetic metabolism as an effective strategy for high-value compound production.

Gene Networks Underlying Cannabinoid and Terpenoid Accumulation in Cannabis

Glandular trichomes are specialized anatomical structures that accumulate secretions with important biological roles in plant-environment interactions. These secretions also have commercial uses in the flavor, fragrance, and pharmaceutical industries. The capitate-stalked glandular trichomes of Cannabis sativa (cannabis), situated on the surfaces of the bracts of the female flowers, are the primary site for the biosynthesis and storage of resins rich in cannabinoids and terpenoids. In this study, we profiled nine commercial cannabis strains with purportedly different attributes, such as taste, color, smell, and genetic origin. Glandular trichomes were isolated from each of these strains, and cell type-specific transcriptome data sets were acquired. Cannabinoids and terpenoids were quantified in flower buds. Statistical analyses indicated that these data sets enable the high-resolution differentiation of strains by providing complementary information. Integrative analyses revealed a coexpression network of genes involved in the biosynthesis of both cannabinoids and terpenoids from imported precursors. Terpene synthase genes involved in the biosynthesis of the major monoterpenes and sesquiterpenes routinely assayed by cannabis testing laboratories were identified and functionally evaluated. In addition to cloning variants of previously characterized genes, specifically CsTPS14CT [(-)-limonene synthase] and CsTPS15CT (β-myrcene synthase), we functionally evaluated genes that encode enzymes with activities not previously described in cannabis, namely CsTPS18VF and CsTPS19BL (nerolidol/linalool synthases), CsTPS16CC (germacrene B synthase), and CsTPS20CT (hedycaryol synthase). This study lays the groundwork for developing a better understanding of the complex chemistry and biochemistry underlying resin accumulation across commercial cannabis strains.

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.

Eukaryotic microalgae as hosts for light-driven heterologous isoprenoid production

Eukaryotic microalgae hold incredible metabolic potential for the sustainable production of heterologous isoprenoid products. Recent advances in algal engineering have enabled the demonstration of prominent examples of heterologous isoprenoid production. Isoprenoids, also known as terpenes or terpenoids, are the largest class of natural chemicals, with a vast diversity of structures and biological roles. Some have high-value in human-use applications, although may be found in their native contexts in low abundance or be difficult to extract and purify. Heterologous production of isoprenoid compounds in heterotrophic microbial hosts such as bacteria or yeasts has been an active area of research for some time and is now a mature technology. Eukaryotic microalgae represent sustainable alternatives to these hosts for biotechnological production processes as their cultivation can be driven by light and freely available CO₂ as a carbon source. Their photosynthetic lifestyles require metabolic architectures structured towards the generation of associated isoprenoids (carotenoids, phytol) which participate in photon capture, energy dissipation, and electron transfer. Eukaryotic microalgae should, therefore, contain inherently high capacities for the generation of heterologous isoprenoid products. Although engineering strategies in eukaryotic microalgae have lagged behind the more genetically tractable bacteria and yeasts, recent advances in algal engineering concepts have demonstrated prominent examples of light-driven heterologous isoprenoid production from these photosynthetic hosts. This work seeks to provide practical insights into the choice of eukaryotic microalgae as biotechnological chassis. Recent reports of advances in algal engineering for heterologous isoprenoid production are highlighted as encouraging examples that promote their expanded use as sustainable green-cell factories. Current state of the art, limitations, and future challenges are also discussed.

Non-volatile natural products in plant glandular trichomes: chemistry, biological activities and biosynthesis

The investigation methods, chemistry, bioactivities, and biosynthesis of non-volatile natural products involving 489 compounds in plant glandular trichomes are reviewed.

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.

Coactivation of MEP-biosynthetic genes and accumulation of abietane diterpenes in Salvia sclarea by heterologous expression of WRKY and MYC2 transcription factors

Plant abietane diterpenoids (e.g. aethiopinone, 1- oxoaethiopinone, salvipisone and ferruginol), synthesized in the roots of several Salvia spp, have antibacterial, antifungal, sedative and anti-proliferative properties. Recently we have reported that content of these compounds in S. sclarea hairy roots is strongly depending on transcriptional regulation of genes belonging to the plastidial MEP-dependent terpenoid pathway, from which they mostly derive. To boost the synthesis of this interesting class of compounds, heterologous AtWRKY18, AtWRKY40, and AtMYC2 TFs were overexpressed in S. sclarea hairy roots and proved to regulate in a coordinated manner the expression of several genes encoding enzymes of the MEP-dependent pathway, especially DXS, DXR, GGPPS and CPPS. The content of total abietane diterpenes was enhanced in all overexpressing lines, although in a variable manner due to a negative pleiotropic effect on HR growth. Interestingly, in the best performing HR lines overexpressing the AtWRKY40 TF induced a significant 4-fold increase in the final yield of aethiopinone, for which we have reported an interesting anti-proliferative activity against resistant melanoma cells. The present results are also informative and instrumental to enhance the synthesis of abietane diterpenes derived from the plastidial MEP-derived terpenoid pathway in other Salvia species.

Assessing Flux Distribution Associated with Metabolic Specialization of Glandular Trichomes

Many aromatic plants accumulate mixtures of secondary (or specialized) metabolites in anatomical structures called glandular trichomes (GTs). Different GT types may also synthesize different mixtures of secreted metabolites, and this contributes to the enormous chemical diversity reported to occur across species. Over the past two decades, significant progress has been made in characterizing the genes and enzymes that are responsible for the unique metabolic capabilities of GTs in different lineages of flowering plants. Less is known about the processes that regulate flux distribution through precursor pathways toward metabolic end-products. We discuss here the results from a meta-analysis of genome-scale models that were developed to capture the unique metabolic capabilities of different GT types.

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

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

Plant Glandular Trichomes: Natural Cell Factories of High Biotechnological Interest

Multicellular glandular trichomes are epidermal outgrowths characterized by the presence of a head made of cells that have the ability to secrete or store large quantities of specialized metabolites. Our understanding of the transcriptional control of glandular trichome initiation and development is still in its infancy. This review points to some central questions that need to be addressed to better understand how such specialized cell structures arise from the plant protodermis. A key and unique feature of glandular trichomes is their ability to synthesize and secrete large amounts, relative to their size, of a limited number of metabolites. As such, they qualify as true cell factories, making them interesting targets for metabolic engineering. In this review, recent advances regarding terpene metabolic engineering are highlighted, with a special focus on tobacco (Nicotiana tabacum). In particular, the choice of transcriptional promoters to drive transgene expression and the best ways to sink existing pools of terpene precursors are discussed. The bioavailability of existing pools of natural precursor molecules is a key parameter and is controlled by so-called cross talk between different biosynthetic pathways. As highlighted in this review, the exact nature and extent of such cross talk are only partially understood at present. In the future, awareness of, and detailed knowledge on, the biology of plant glandular trichome development and metabolism will generate new leads to tap the largely unexploited potential of glandular trichomes in plant resistance to pests and lead to the improved production of specialized metabolites with high industrial or pharmacological value.

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

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

Functional Diversification of Kaurene Synthase-Like Genes in Isodon rubescens

Ent-kaurene diterpenoids are the largest group of known Isodon diterpenoids. Among them, oridonin is accumulated in the leaves, and is the most frequently studied compound because of its antitumor and antibacterial activities. We have identified five copalyl diphosphate synthase (CPS) and six kaurene synthase-like (KSL) genes by transcriptome profiling of Isodon rubescens leaves. An in vitro assay assigns ten of them to five different diterpene biosynthesis pathways, except IrCPS3 that has a mutation in the catalytic motif. The Lamiaceae-specific clade genes (IrCPS1 and IrCPS2) synthesize the intermediate copalyl diphosphate (normal-CPP), while IrCPS4 and IrCPS5 synthesize the intermediate ent-copalyl diphosphate (ent-CPP). IrKSL2, IrKSL4, and IrKSL5 react with ent-CPP to produce an ent-isopimaradiene-like compound, ent-atiserene and ent-kaurene, respectively. Correspondingly, the Lamiaceae-specific clade genes IrKSL1 or IrKSL3 combined with normal-CPP led to the formation of miltiradiene. The compound then underwent aromatization and oxidization with a cytochrome P450 forming two related compounds, abietatriene and ferruginol, which were detected in the root bark. IrKSL6 reacts with normal-CPP to produce isopimaradiene. IrKSL3 and IrKSL6 have the γβα tridomain structure, as these proteins tend to possess the bidomain structure of IrKSL1, highlighting the evolutionary history of KSL gene domain loss and further elucidating chemical diversity evolution from a macroevolutionary stance in Lamiaceae.

Eigenvector metabolite analysis reveals dietary effects on the association among metabolite correlation patterns, gene expression, and phenotypes

INTRODUCTION

'Multi-omics' datasets obtained from an organism of interest reared under different environmental treatments are increasingly common. Identifying the links among metabolites and transcripts can help to elucidate our understanding of the impact of environment at different levels within the organism. However, many methods for characterizing physiological connections cannot address unidentified metabolites.

OBJECTIVES

Here, we use Eigenvector Metabolite Analysis (EvMA) to examine links between metabolomic, transcriptomic, and phenotypic variation data and to assess the impact of environmental factors on these associations. Unlike other methods, EvMA can be used to analyze datasets that include unidentified metabolites and unannotated transcripts.

METHODS

To demonstrate the utility of EvMA, we analyzed metabolomic, transcriptomic, and phenotypic datasets produced from 20 Drosophila melanogaster genotypes reared on four dietary treatments. We used a hierarchical distance-based method to cluster the metabolites. The links between metabolite clusters, gene expression, and overt phenotypes were characterized using the eigenmetabolite (first principal component) of each cluster.

RESULTS

EvMA recovered chemically related groups of metabolites within the clusters. Using the eigenmetabolite, we identified genes and phenotypes that significantly correlated with each cluster. EvMA identifies new connections between the phenotypes, metabolites, and gene transcripts. Conclusion EvMA provides a simple method to identify correlations between metabolites, gene expression, and phenotypes, which can allow us to partition multivariate datasets into meaningful biological modules and identify under-studied metabolites and unannotated gene transcripts that may be central to important biological processes. This can be used to inform our understanding of the effect of environmental mechanisms underlying physiological states of interest.

Carnosic acid biosynthesis elucidated by a synthetic biology platform

Synthetic biology approaches achieving the reconstruction of specific plant natural product biosynthetic pathways in dedicated microbial "chassis" have provided access to important industrial compounds (e.g., artemisinin, resveratrol, vanillin). However, the potential of such production systems to facilitate elucidation of plant biosynthetic pathways has been underexplored. Here we report on the application of a modular terpene production platform in the characterization of the biosynthetic pathway leading to the potent antioxidant carnosic acid and related diterpenes in Salvia pomifera and Rosmarinus officinalis.Four cytochrome P450 enzymes are identified (CYP76AH24, CYP71BE52, CYP76AK6, and CYP76AK8), the combined activities of which account for all of the oxidation events leading to the biosynthesis of the major diterpenes produced in these plants. This approach develops yeast as an efficient tool to harness the biotechnological potential of the numerous sequencing datasets that are increasingly becoming available through transcriptomic or genomic studies.

Production of the forskolin precursor 11β-hydroxy-manoyl oxide in yeast using surrogate enzymatic activities

Background

Several plant diterpenes have important biological properties. Among them, forskolin is a complex labdane-type diterpene whose biological activity stems from its ability to activate adenylyl cyclase and to elevate intracellular cAMP levels. As such, it is used in the control of blood pressure, in the protection from congestive heart failure, and in weight-loss supplements. Chemical synthesis of forskolin is challenging, and production of forskolin in engineered microbes could provide a sustainable source. To this end, we set out to establish a platform for the production of forskolin and related epoxy-labdanes in yeast.

Results

Since the forskolin biosynthetic pathway has only been partially elucidated, and enzymes involved in terpene biosynthesis frequently exhibit relaxed substrate specificity, we explored the possibility of reconstructing missing steps of this pathway employing surrogate enzymes. Using CYP76AH24, a Salvia pomifera cytochrome P450 responsible for the oxidation of C-12 and C-11 of the abietane skeleton en route to carnosic acid, we were able to produce the forskolin precursor 11β-hydroxy-manoyl oxide in yeast. To improve 11β-hydroxy-manoyl oxide production, we undertook a chassis engineering effort involving the combination of three heterozygous yeast gene deletions (mct1/MCT1, whi2/WHI2, gdh1/GDH1) and obtained a 9.5-fold increase in 11β-hydroxy-manoyl oxide titers, reaching 21.2 mg L⁻¹.

Conclusions

In this study, we identify a surrogate enzyme for the specific and efficient hydroxylation of manoyl oxide at position C-11β and establish a platform that will facilitate the synthesis of a broad range of tricyclic (8,13)-epoxy-labdanes in yeast. This platform forms a basis for the heterologous production of forskolin and will facilitate the elucidation of subsequent steps of forskolin biosynthesis. In addition, this study highlights the usefulness of using surrogate enzymes for the production of intermediates of complex biosynthetic pathways. The combination of heterozygous deletions and the improved yeast strain reported here will provide a useful tool for the production of numerous other isoprenoids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0440-8) contains supplementary material, which is available to authorized users.