Functional characterization of nine Norway Spruce TPS genes and evolution of gymnosperm terpene synthases of the TPS-d subfamily

Z/E configuration controlled by a Taxus sesquiterpene synthase facilitating the biosynthesis of (3Z,6E)-α-farnesene

Plant enzymes often present advantages in the synthesis of natural products with specific configurations. Farnesene is a pharmacologically active sesquiterpene with three natural Z/E configurations, among which the enzyme selectively responsible for the biosynthesis of (3Z,6E)-α-farnesene remains elusive. Herein, a sesquiterpene synthase TwSTPS1 biosynthesizing (3Z,6E)-α-farnesene as the major product was identified from Taxus wallichiana through genome mining. Utilizing molecular dynamics simulations and mutation analysis, the catalytic mechanism of TwSTPS1, especially Z/E configuration control, was explored. Moreover, the crucial residues associated with the specific catalytic activity of TwSTPS1 was elucidated through mutagenesis experiments. The findings contribute to our understanding of the Z/E configuration control by plant terpene synthases and also provide an alternative tool for manipulating (3Z,6E)-α-farnesene production using synthetic biology.

Plant specialized metabolism: Diversity of terpene synthases and their products

Terpenoids are ubiquitous to all kingdoms of life and are one of the most diverse groups of compounds, both structurally and functionally. Despite being derived from common precursors, isopentenyl diphosphate and dimethylallyl diphosphate, their exceptional diversity is partly driven by the substrate and product promiscuity of terpene synthases that produce a wide array of terpene skeletons. Plant terpene synthases can be subdivided into different subfamilies based on sequence homology and function. However, in many cases, structural architecture of the enzyme is more essential to product specificity than primary sequence alone, and distantly related terpene synthases can often mediate similar reactions. As such, the focus of this brief review is on some of the recent progress in understanding terpene synthase function and diversity.

Integrating GC-MS and comparative transcriptome analysis reveals that TsERF66 promotes the biosynthesis of caryophyllene in Toona sinensis tender leaves

Introduction

The strong aromatic characteristics of the tender leaves of Toona sinensis determine their quality and economic value.

Methods and results

Here, GC-MS analysis revealed that caryophyllene is a key volatile compound in the tender leaves of two different T. sinensis varieties, however, the transcriptional mechanisms controlling its gene expression are unknown. Comparative transcriptome analysis revealed significant enrichment of terpenoid synthesis pathway genes, suggesting that the regulation of terpenoid synthesis-related gene expression is an important factor leading to differences in aroma between the two varieties. Further analysis of expression levels and genetic evolution revealed that TsTPS18 is a caryophyllene synthase, which was confirmed by transient overexpression in T. sinensis and Nicotiana benthamiana leaves. Furthermore, we screened an AP2/ERF transcriptional factor ERF-IX member, TsERF66, for the potential regulation of caryophyllene synthesis. The TsERF66 had a similar expression trend to that of TsTPS18 and was highly expressed in high-aroma varieties and tender leaves. Exogenous spraying of MeJA also induced the expression of TsERF66 and TsTPS18 and promoted the biosynthesis of caryophyllene. Transient overexpression of TsERF66 in T. sinensis significantly promoted TsTPS18 expression and caryophyllene biosynthesis.

Discussion

Our results showed that TsERF66 promoted the expression of TsTPS18 and the biosynthesis of caryophyllene in T. sinensis leaves, providing a strategy for improving the aroma of tender leaves.

Biochemical analysis of the TPS-a subfamily in Medicago truncatula

Terpenes are important mediators of plant chemical response to environmental cues. Here, we describe the genome-wide identification and biochemical characterization of TPS-a members in Medicago truncatula, a model legume crop. Genome mining identified thirty-nine full-length terpene synthases with a significant number predicted to produce monoterpenes and sesquiterpenes. Biochemical characterization of the TPS-a subfamily associated with sesquiterpene biosynthesis revealed such compounds, that exhibit substantial biological activity in other plants. Gene expression analysis using qPCR and the Medicago gene atlas illustrated distinct tissue and time-based variation in expression in leaves and roots. Together our work establishes the gene-to-metabolite relationships for sesquiterpene synthases in M. truncatula. Understanding the biosynthetic capacity is a foundational step to defining the ecological roles of this important family of compounds.

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.

MYB24 orchestrates terpene and flavonol metabolism as light responses to anthocyanin depletion in variegated grape berries

Variegation is a rare type of mosaicism not fully studied in plants, especially fruits. We examined red and white sections of grape (Vitis vinifera cv. 'Béquignol') variegated berries and found that accumulation of products from branches of the phenylpropanoid and isoprenoid pathways showed an opposite tendency. Light-responsive flavonol and monoterpene levels increased in anthocyanin-depleted areas in correlation with increasing MYB24 expression. Cistrome analysis suggested that MYB24 binds to the promoters of 22 terpene synthase (TPS) genes, as well as 32 photosynthesis/light-related genes, including carotenoid pathway members, the flavonol regulator HY5 HOMOLOGUE (HYH), and other radiation response genes. Indeed, TPS35, TPS09, the carotenoid isomerase gene CRTISO2, and HYH were activated in the presence of MYB24 and MYC2. We suggest that MYB24 modulates ultraviolet and high-intensity visible light stress responses that include terpene and flavonol synthesis and potentially affects carotenoids. The MYB24 regulatory network is developmentally triggered after the onset of berry ripening, while the absence of anthocyanin sunscreens accelerates its activation, likely in a dose-dependent manner due to increased radiation exposure. Anthocyanins and flavonols in variegated berry skins act as effective sunscreens but for different wavelength ranges. The expression patterns of stress marker genes in red and white sections of 'Béquignol' berries strongly suggest that MYB24 promotes light stress amelioration but only partly succeeds during late ripening.

Characterization of sandalwood (E,E)-α-farnesene synthase whose overexpression enhances cold tolerance through jasmonic acid biosynthesis and signaling in Arabidopsis

MAIN CONCLUSION

Santalum album (E,E)-α-farnesene synthase catalyzes FPP into (E,E)-α-farnesene. Overexpression of the SaAFS gene positively improved cold stress tolerance through JA biosynthesis and signaling pathways in Arabidopsis. Volatile terpenoids are released from plants that suffer negative effects following exposure to various biotic and abiotic stresses. Recent studies revealed that (E,E)-α-farnesene synthase (AFS) plays a significant role in a plant's defence against biotic attack. However, little is known about whether AFS contributes to plant resistance to cold stress. In this study, a SaAFS gene was isolated from Indian sandalwood (Santalum album L.) and functionally characterized. The SaAFS protein mainly converts farnesyl diphosphate to (E,E)-α-farnesene. SaAFS was clustered into the AFS clade from angiosperms, suggesting a highly conserved enzyme. SaAFS displayed a significant response to cold stress and methyl jasmonate. SaAFS overexpression (OE) in Arabidopsis enhanced cold tolerance by increasing proline content, reducing malondialdehyde content, electrolyte leakage, and accumulating reactive oxygen species. Transcriptomic analysis revealed that upregulated genes related to stress response and JA biosynthesis and signaling were detected in SaAFS-OE lines compared with wild type plants that were exposed to cold stress. Endogenous JA and jasmonoyl-isoleucine content increased significantly in SaAFS-OE lines exposed to cold stress. Collectively considered, these results suggest that the SaAFS gene is a positive regulator during cold stress tolerance via JA biosynthesis and signaling pathways.

Genome-wide identification of terpenoid synthase family genes in Gossypium hirsutum and functional dissection of its subfamily cadinene synthase A in gossypol synthesis

Plant terpenoid synthase (TPS) family genes participate in metabolite synthesis, hormones, gossypol, etc. Here, we genome-widely identified TPS family genes in 12 land plant species. Four hundred and thirty TPS-related genes were divided into seven subfamilies. The TPS-c in Bryophytes was suggested to be the earliest subfamily, followed by the TPS-e/f and TPS-h presence in ferns. TPS-a, the largest number of genes, was derived from monocotyledonous and dicotyledonous plants. Collinearity analysis showed that 38 out of the 76 TPS genes in G. hirsutum were collinear within G. arboreum and G. raimondii. Twenty-one GhTPS-a genes belong to the cadinene synthase (GhCDN) subfamily and were divided into five groups, A, B, C, D, and E. The special cis-elements in the promoters of 12 GhCDN-A genes suggested that the JA and ethylene signaling pathways may be involved in their expression regulation. When 12 GhCDN-A genes were simultaneously silenced through virus-induced gene silencing, the glandular color of GhCDN-A-silenced plants was lighter than that of the control, supported by a gossypol content decrease based on HPLC testing, suggesting that GhCDN-A subgroup genes participate in gossypol synthesis. According to RNA-seq analysis, gossypol synthesis-related genes and disease-resistant genes in the glandular variety exhibited upregulated expression compared to the glandless variety, whereas hormone signaling-related genes were downregulated. All in all, these results revealed plant TPS gene evolution rules and dissected the TPS subfamily, GhCDN-A, function in gossypol synthesis in cotton.

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

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

Dynamics of internal isoprenoid metabolites in young Picea abies (Norway spruce) shoots during drought stress conditions in springtime

Currently, large areas of Picea abies (Norway spruce) stands in Europe are increasingly affected by drought and heat waves. Moreover, early spring drought has occurred with much higher frequency. Our work focuses on physiological changes induced by drought in four-year-old spruce seedlings during shoot elongation. We investigated drought effect on photosynthetic rate, concentration of abscisic acid and its metabolites, amount and composition of monoterpenes in needles of seedlings from five different provenances (altitude range 550-1280 m above sea level) in Western Carpathians. Spruce seedlings subjected to one-month drought stress of moderate intensity (about 50% of soil water content at the end of experiment) showed significant reduction of CO₂ uptake and increased concentration of ABA related to untreated controls. Induced drought affected needle monoterpene content and composition. Observed changes in drought-induced physiological parameters were influenced by seedling provenance. The provenance from 920 m above sea level showed the greatest sensitivity to drought with significantly highest ABA content and, at the same time, a clear decline of CO₂ uptake and amounts of total monoterpenes. Our results indicating intra-specific provenance-related variability in physiological response of spruce seedlings to drought may provide a basis for improved reforestation strategies in drought risk areas.

Metabolome and transcriptome analysis of terpene synthase genes and their putative role in floral aroma production in Litchi chinensis

Volatile organic compounds (VOCs) are essential traits of flowers since they attract pollinators, aid in seed distribution, protect the plant from internal and external stimuli, and are involved in plant-plant and plant-environment interactions. Apart from their role in plants, VOCs are used in pharmaceuticals, fragrances, cosmetics, and flavorings. Litchi (Litchi chinensis Sonn.) is a popular fruit due to its enticing red appearance, exotic taste, and high nutritional qualities. Litchi flowers bloom as inflorescences primarily on the shoot terminals. There are three distinct flower types, two male and one female, all of which are produced on the same panicle and rely on insect pollination. Herein, we used a comprehensive metabolomic approach to examine the volatile profile of litchi fruit (green pericarp, yellow pericarp, and red pericarp) as well as male and female flowers (bud stage, half open and full bloom). From a quantitative examination of the volatiles in L. chinensis, a total of 19, 22, and 21 VOCs were discovered from female flowers, male flowers, and fruits, with the majority of them belonging to sesquiterpenes. Multivariate analysis revealed that the volatile profiles of fruits differ from those of male and female flowers. Three VOCs were unique to male flowers and ten to the fruit, while eight VOCs were shared by both male and female flowers and eleven by both male and female flowers and the fruit. Furthermore, for the first time, we identified and comprehensively studied the TERPENE SYNTHASE genes (TPS) using the litchi genome and transcriptome database, which revealed 38 TPS genes unevenly distributed across the 15 chromosomes. A phylogenetic study showed that LcTPS were grouped into TPS-b, TPS-c, TPS-e, TPS-f, and TPS-g subfamilies, with TPS-b having the most genes. The conserved motifs (RR_X8 W, NSE/DTE, and DD_XX D) were studied in LcTPSs, and significant variation between subfamilies was discovered. Furthermore, after integrating the metabolome and transcriptome datasets, several VOCs were shown to be development-specific and highly linked with distinct LcTPS genes, making them promising biomarkers. Interestingly, LcTPS17/20/23/24/31 were associated with monoterpene edges, while the rest were connected to sesquiterpene edges, indicating their probable participation in the aroma biosynthesis mechanism of certain compounds.

Full-Length Transcriptomic Sequencing and Temporal Transcriptome Expression Profiling Analyses Offer Insights into Terpenoid Biosynthesis in Artemisia argyi

Artemisiae argyi Folium is a traditional herbal medicine used for moxibustion heat therapy in China. The volatile oils in A.argyi leaves are closely related to its medicinal value. Records suggest that the levels of these terpenoids components within the leaves vary as a function of harvest time, with June being the optimal time for A. argyi harvesting, owing to the high levels of active ingredients during this month. However, the molecular mechanisms governing terpenoid biosynthesis and the time-dependent changes in this activity remain unclear. In this study, GC–MS analysis revealed that volatile oil levels varied across four different harvest months (April, May, June, and July) in A. argyi leaves, and the primarily terpenoids components (including both monoterpenes and sesquiterpenes) reached peak levels in early June. Through single-molecule real-time (SMRT) sequencing, corrected by Illumina RNA-sequencing (RNA-Seq), 44 full-length transcripts potentially involved in terpenoid biosynthesis were identified in this study. Differentially expressed genes (DEGs) exhibiting time-dependent expression patterns were divided into 12 coexpression clusters. Integrated chemical and transcriptomic analyses revealed distinct time-specific transcriptomic patterns associated with terpenoid biosynthesis. Subsequent hierarchical clustering and correlation analyses ultimately identified six transcripts that were closely linked to the production of these two types of terpenoid within A. argyi leaves, revealing that the structural diversity of terpenoid is related to the generation of the diverse terpene skeletons by prenyltransferase (TPS) family of enzymes. These findings can guide further studies of the molecular mechanisms underlying the quality of A. argyi leaves, aiding in the selection of optimal timing for harvests of A. argyi.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Genome of Lindera glauca provides insights into the evolution of biosynthesis genes for aromatic compounds

Lindera glauca is a crucial source of diverse industrial oil and medicines. The spicy aroma of tender leaves is caused by the presence of abundant aromatic compounds. Here, we present its chromosome-level genome assembly comprising 12 pseudochromosomes (2,092.2 Mb; scaffold N50: 186.5 Mb), which was predicted to have 65,145 protein-coding genes. Comparative genomic analyses indicated two whole-genome duplication (WGD) events in the Lauraceae family, contributing to the production of numerous terpene synthase (TPS) genes. We identified 138 TPS genes in L. glauca. Comparative transcriptomic analyses revealed high expression of genes Lg03G2346 and Lg08G140 in TPS-a and Lg07G2961 and Lg12G971 in TPS-b subfamilies, which regulated the biosynthesis of the monoterpenoid β-ocimene and sesquiterpenoid D-germacrene in L. glauca. The results suggested a molecular basis for species-specific terpenoid biosynthesis and provided a foundation for molecular breeding to produce desired characteristics and a valuable reference genome.

Cloning and functional characterization of three sesquiterpene synthase genes from Chamaecyparis formosensis Matsumura

Terpene synthase (TPS) analysis may contribute to a better understanding of terpenoids biosynthesis and the evolution of phylogenetic taxonomy. Chamaecyparis formosensis Matsumura is an endemic and valuable conifer of Taiwan. Its excellent wood quality, fragrance, and durability make it become the five precious conifers in Taiwan. In this study, three sesquiterpene synthase genes that belong to the TPS-d2 clade were isolated and characterized through in vitro reaction of recombinant protein and in vivo reaction of Escherichia coli heterologous expression system. The main product of Cf-GerA was germacrene A using GC/MS analysis, while the product of Cf-Aco and Cf-Gor were identified as acora-4(14),8-diene and (5R,6R,10S)-α-gorgonene by using NMR analysis. These are the first reported enzymes that biosynthesize acora-4(14),8-diene and (5 R,6 R,10 S)-α-gorgonene. Both sesquiterpene synthases may isomerize the farnesyl pyrophosphate substrate to nerolidyl pyrophosphate for further cyclization. Cf-Aco may catalyze 1,6-cyclization of nerolidyl cation while Cf-Gor may catalyze through an uncharged intermediate, isogermacrene A.

Identification and Defensive Characterization of PmCYP720B11v2 from Pinus massoniana

Pinus massoniana is a pioneer species for afforestation timber and oleoresin, while epidemics of pinewood nematode (PWN; Bursaphelenchus xylophilus) are causing a serious biotic disaster for P. massoniana in China. Importantly, resistant P. massoniana could leak copious oleoresin terpenoids to build particular defense fronts for survival when attacked by PWN. However, the defense mechanisms regulating this process remain unknown. Here, PmCYP720B11v2, a cytochrome P450 monooxygenase gene, was first identified and functionally characterized from resistant P. massoniana following PWN inoculation. The tissue-specific expression pattern and localization of PmCYP720B11v2 at the transcript and protein levels in resistant P. massoniana indicated that its upregulation in the stem supported its involvement in the metabolic processes of diterpene biosynthesis as a positive part of the defense against PWN attack. Furthermore, overexpression of PmCYP720B11v2 may enhance the growth and development of plants. In addition, PmCYP720B11v2 activated the metabolic flux of antioxidases and stress-responsive proteins under drought conditions and improved drought stress tolerance. Our results provide new insights into the favorable role of PmCYP720B11v2 in diterpene defense mechanisms in response to PWN attack in resistant P. massoniana and provide a novel metabolic engineering scenario to reform the stress tolerance potential of tobacco.

Bark Beetle Attack History Does Not Influence the Induction of Terpene and Phenolic Defenses in Mature Norway Spruce (Picea abies) Trees by the Bark Beetle-Associated Fungus Endoconidiophora polonica

Terpenes and phenolics are important constitutive and inducible conifer defenses against bark beetles and their associated fungi. In this study, the inducible defenses of mature Norway spruce (Picea abies) trees with different histories of attack by the spruce bark beetle, Ips typographus were tested by inoculation with the I. typographus-associated fungus Endoconidiophora polonica. We compared trees that had been under previous attack with those under current attack and those that had no record of attack. After fungal inoculation, the concentrations of mono-, sesqui-, and diterpenes in bark increased 3- to 9-fold. For the phenolics, the flavan-3-ols, catechin, and gallocatechin, increased significantly by 2- and 5-fold, respectively, while other flavonoids and stilbenes did not. The magnitudes of these inductions were not influenced by prior bark beetle attack history for all the major compounds and compound classes measured. Before fungal inoculation, the total amounts of monoterpenes, diterpenes, and phenolics (constitutive defenses) were greater in trees that had been previously attacked compared to those under current attack, possibly a result of previous induction. The transcript levels of many genes involved in terpene formation (isoprenyl diphosphate synthases and terpene synthases) and phenolic formation (chalcone synthases) were significantly enhanced by fungal inoculation suggesting de novo biosynthesis. Similar inductions were found for the enzymatic activity of isoprenyl diphosphate synthases and the concentration of their prenyl diphosphate products after fungal inoculation. Quantification of defense hormones revealed a significant induction of the jasmonate pathway, but not the salicylic acid pathway after fungal inoculation. Our data highlight the coordinated induction of terpenes and phenolics in spruce upon infection by E. polonica, a fungal associate of the bark beetle I. typographus, but provide no evidence for the priming of these defense responses by prior beetle attack.

A Chromosome-Level Genome of the Camphor Tree and the Underlying Genetic and Climatic Factors for Its Top-Geoherbalism

Camphor tree [Cinnamomum camphora (L.) J. Presl], a species in the magnoliid family Lauraceae, is known for its rich volatile oils and is used as a medical cardiotonic and as a scent in many perfumed hygiene products. Here, we present a high-quality chromosome-scale genome of C. camphora with a scaffold N50 of 64.34 Mb and an assembled genome size of 755.41 Mb. Phylogenetic inference revealed that the magnoliids are a sister group to the clade of eudicots and monocots. Comparative genomic analyses identified two rounds of ancient whole-genome duplication (WGD). Tandem duplicated genes exhibited a higher evolutionary rate, a more recent evolutionary history and a more clustered distribution on chromosomes, contributing to the production of secondary metabolites, especially monoterpenes and sesquiterpenes, which are the principal essential oil components. Three-dimensional analyses of the volatile metabolites, gene expression and climate data of samples with the same genotype grown in different locations showed that low temperature and low precipitation during the cold season modulate the expression of genes in the terpenoid biosynthesis pathways, especially TPS genes, which facilitates the accumulation of volatile compounds. Our study lays a theoretical foundation for policy-making regarding the agroforestry applications of camphor tree.

Overexpression of geranyl diphosphate synthase (PmGPPS1) boosts monoterpene and diterpene production involved in the response to pine wood nematode invasion

Outbreaks of pine wood nematode (PWN; Bursaphelenchus xylophilus) represent a severe biotic epidemic for the Pinus massoniana in China. When invaded by the PWN, the resistant P. massoniana might secret abundant oleoresin terpenoid to form certain defensive fronts for survival. However, the regulatory mechanisms of this process remain unclear. Here, the geranyl diphosphate synthase (PmGPPS1) gene was identified from resistant P. massoniana. Tissue-specific expression patterns of PmGPPS1 at transcript and protein level in resistant P. massoniana were determined by quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry. Functional characteristics analysis of PmGPPS1 was performed on transgenic Nicotiana benthamiana by overexpression, as genetic transformation of P. massoniana is, so far, not possible. In summary, we identified and functionally characterized PmGPPS1 from the resistant P. massoniana following PWN inoculation. Tissue-specific expression patterns and localization of PmGPPS1 indicated that it may play a positive role involved in the metabolic and defensive processes of oleoresin terpenes production in response to PWN attack. Furthermore, overexpression of PmGPPS1 may enhance the production of monoterpene, among which limonene reduced the survival of PWN in vitro. In addition, PmGPPS1 upregulated the expression level of key genes involved in mevalonic acid (MVA) pathway, the methylerythritol phosphate (MEP) pathway and gibberellins (GAs) biosynthesis to boost the growth and development of tobacco through a feedback regulation mechanism. Our results offered new insights into the pivotal role of the PmGPPS1 involved in terpene-based defense mechanisms responding to the PWN invasion in resistant P. massoniana and provided a new metabolic engineering scenario to improve monoterpene production in tobacco.

Monoterpene Synthase Genes and Monoterpene Profiles in Pinus nigra subsp. laricio

In the present study, we carried out a quantitative analysis of the monoterpenes composition in different tissues of the non-model conifer Pinus nigra J.F. Arnold subsp. laricio Palib. ex Maire (P. laricio, in short). All the P. laricio tissues examined showed the presence of the same fourteen monoterpenes, among which the most abundant were β-phellandrene, α-pinene, and β-pinene, whose distribution was markedly tissue-specific. In parallel, from the same plant tissues, we isolated seven full-length cDNA transcripts coding for as many monoterpene synthases, each of which was found to be attributable to one of the seven phylogenetic groups in which the d1-clade of the canonical classification of plants' terpene synthases can be subdivided. The amino acid sequences deduced from the above cDNA transcripts allowed to predict their putative involvement in the biosynthesis of five of the monoterpenes identified. Transcripts profiling revealed a differential gene expression across the different tissues examined, and was found to be consistent with the corresponding metabolites profiles. The genomic organization of the seven isolated monoterpene synthase genes was also determined.

Transcriptional up-regulation of host-specific terpene metabolism in aphid-induced galls of Pistacia palaestina

Galling insects gain food and shelter by inducing specialized anatomical structures in their plant hosts. Such galls often accumulate plant defensive metabolites protecting the inhabiting insects from predation. We previously found that, despite a marked natural chemopolymorphism in natural populations of Pistacia palaestina, the monoterpene content in Baizongia pistaciae-induced galls is substantially higher than in leaves of their hosts. Here we show a general up-regulation of key structural genes in both the plastidial and cytosolic terpene biosynthetic pathways in galls as compared with non-colonized leaves. Novel prenyltransferases and terpene synthases were functionally expressed in Escherichia coli to reveal their biochemical function. Individual Pistacia trees exhibiting chemopolymorphism in terpene compositions displayed differential up-regulation of selected terpene synthase genes, and the metabolites generated by their gene products in vitro corresponded to the monoterpenes accumulated by each tree. Our results delineate molecular mechanisms responsible for the formation of enhanced monoterpene in galls and the observed intraspecific monoterpene chemodiversity displayed in P. palaestina. We demonstrate that gall-inhabiting aphids transcriptionally reprogram their host terpene pathways by up-regulating tree-specific genes, boosting the accumulation of plant defensive compounds for the protection of colonizing insects.

Recent Advances in the Biosynthesis of Farnesene Using Metabolic Engineering

Farnesene, as an important sesquiterpene isoprenoid polymer of acetyl-CoA, is a renewable feedstock for diesel fuel, polymers, and cosmetics. It has been widely applied in agriculture, medicine, energy, and other fields. In recent years, farnesene biosynthesis is considered a green and economical approach because of its mild reaction conditions, low environmental pollution, and sustainability. Metabolic engineering has been widely applied to construct cell factories for farnesene biosynthesis. In this paper, the research progress, common problems, and strategies of farnesene biosynthesis are reviewed. They are mainly described from the perspectives of the current status of farnesene biosynthesis in different host cells, optimization of the metabolic pathway for farnesene biosynthesis, and key enzymes for farnesene biosynthesis. Furthermore, the challenges and prospects for future farnesene biosynthesis are discussed.

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

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

Unveiling Monoterpene Biosynthesis in Taiwania cryptomerioides via Functional Characterization

Taiwania cryptomerioides is a monotypic species, and its terpenoid-rich property has been reported in recent years. To uncover monoterpene biosynthesis in T. cryptomerioides, this study used transcriptome mining to identify candidates with tentative monoterpene synthase activity. Along with the phylogenetic analysis and in vitro assay, two geraniol synthases (TcTPS13 and TcTPS14), a linalool synthase (TcTPS15), and a β-pinene synthase (TcTPS16), were functionally characterized. Via the comparison of catalytic residues, the Cys/Ser at region 1 might be crucial in determining the formation of α-pinene or β-pinene. In addition, the Cupressaceae monoterpene synthases were phylogenetically clustered together; they are unique and different from those of published conifer species. In summary, this study aimed to uncover the ambiguous monoterpenoid network in T. cryptomerioide, which would expand the landscape of monoterpene biosynthesis in Cupressaceae species.

Diterpene Resin Acids and Olefins in Calabrian Pine (Pinus nigra subsp. laricio (Poiret) Maire) Oleoresin: GC-MS Profiling of Major Diterpenoids in Different Plant Organs, Molecular Identification and Expression Analysis of Diterpene Synthase Genes

A quali-quantitative analysis of diterpenoid composition in tissues obtained from different organs of Pinus nigra subsp. laricio (Poiret) Maire (Calabrian pine) was carried out. Diterpene resin acids were the most abundant diterpenoids across all the examined tissues. The same nine diterpene resin acids were always found, with the abietane type prevailing on the pimarane type, although their quantitative distribution was found to be remarkably tissue-specific. The scrutiny of the available literature revealed species specificity as well. A phylogeny-based approach allowed us to isolate four cDNAs coding for diterpene synthases in Calabrian pine, each of which belonging to one of the four groups into which the d3 clade of the plants' terpene synthases family can be divided. The deduced amino acid sequences allowed predicting that both monofunctional and bifunctional diterpene synthases are involved in the biosynthesis of diterpene resin acids in Calabrian pine. Transcript profiling revealed differential expression across the different tissues and was found to be consistent with the corresponding diterpenoid profiles. The isolation of the complete genomic sequences and the determination of their exon/intron structures allowed us to place the diterpene synthase genes from Calabrian pine on the background of current ideas on the functional evolution of diterpene synthases in Gymnosperms.

Mining of the Catharanthus roseus Genome Leads to Identification of a Biosynthetic Gene Cluster for Fungicidal Sesquiterpenes

Characterization of cryptic biosynthetic gene clusters (BGCs) from microbial genomes has been proven to be a powerful approach to the discovery of new natural products. However, such a genome mining approach to the discovery of bioactive plant metabolites has been muted. The plant BGCs characterized to date encode pathways for antibiotics important in plant defense against microbial pathogens, providing a means to discover such phytoalexins by mining plant genomes. Here is reported the discovery and characterization of a minimal BGC from the medicinal plant Catharanthus roseus, consisting of an adjacent pair of genes encoding a terpene synthase (CrTPS18) and cytochrome P450 (CYP71D349). These two enzymes act sequentially, with CrTPS18 acting as a sesquiterpene synthase, producing 5-epi-jinkoheremol (1), which CYP71D349 further hydroxylates to debneyol (2). Infection studies with maize revealed that 1 and 2 exhibit more potent fungicidal activity than validamycin. Accordingly, this study demonstrates that characterization of such cryptic plant BGCs is a promising strategy for the discovery of potential agrochemical leads. Moreover, despite the observed absence of 1 and 2 in C. roseus, the observed transcriptional regulation is consistent with their differential fungicidal activity, suggesting that such conditional coexpression may be sufficient to drive BGC assembly in plants.

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.

The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Magnolia officinalis, a representative tall aromatic tree of the Magnoliaceae family, is a medicinal plant that is widely used in diverse industries from medicine to cosmetics. We report a chromosome-scale draft genome of M. officinalis, in which ∼99.66% of the sequences were anchored onto 19 chromosomes with the scaffold N50 of 76.62 Mb. We found that a high proportion of repetitive sequences was a common feature of three Magnoliaceae with known genomic data. Magnoliids were a sister clade to eudicots-monocots, which provided more support for understanding the phylogenetic position among angiosperms. An ancient duplication event occurred in the genome of M. officinalis and was shared with Lauraceae. Based on RNA-seq analysis, we identified several key enzyme-coding gene families associated with the biosynthesis of lignans in the genome. The construction of the M. officinalis genome sequence will serve as a reference for further studies of Magnolia, as well as other Magnoliaceae.

Study of Terpenoid Synthesis and Prenyltransferase in Roots of Rehmannia glutinosa Based on iTRAQ Quantitative Proteomics

Rehmannia glutinosa has important medicinal value; terpenoid is one of the main active components in R. glutinosa. In this study, iTRAQ technique was used to analyze the relative abundance of proteins in roots of R. glutinosa, and 6,752 reliable proteins were quantified. GO enrichment results indicated that most proteins were involved in metabolic process or cellular process, 57.63% proteins had catalytic activity, and 65.80% proteins were enriched in membrane-bounded organelle. In roots of R. glutinosa, there were 38 KEGG enrichments with significance, more DEPs were found in some pathways, especially the proteasome pathway and TCA cycle with 15.0% DEPs between elongation stage and expansion stage of roots. Furthermore, five KEGG pathways of terpenoid synthesis were found. Most prenyltransferases belong to FPP/GGPP synthase family, involved in terpenoid backbone biosynthesis, and all interacted with biotin carboxylase CAC2. Compared with that at the elongation stage, many prenyltransferases exhibited higher expression at the expansion stage or maturation stage of roots. In addition, eight FPP/GGPP synthase encoding genes were cloned from R. glutinosa, namely FPPS, FPPS1, GGPS, GGPS3, GGPS4, GGPS5, GPPS and GPPS2, introns were also found in FPPS, FPPS1, GGPS5 and GGPS2, and FPP/GPP synthases were more conservative in organisms, especially in viridiplantae, in which the co-occurrence of GPPS or GPPS2 was significantly higher in plants. Further analysis found that FPP/GGPP synthases of R. glutinosa were divided into three kinds, GGPS, GPPS and FPPS, and their gene expression was significantly diverse in different varieties, growth periods, or tissues of R. glutinosa. Compared with that of GGPS, the expression of GPPS and FPPS was much higher in R. glutinosa, especially at the expansion stage and maturation stage. Thus, the synthesis of terpenoids in roots of R. glutinosa is intricately regulated and needs to be further studied.

Transcriptome sequencing and functional characterization of new sesquiterpene synthases from Curcuma wenyujin

Tubers of Curcuma wenyujin are rich in essential oils, mainly various sesquiterpenes, showing antibacterial, anti-viral and anti-tumor effects. However, the molecular mechanism of C. wenyujin is deficient and related sesquiterpene synthases are still unclear. In this study, the transcriptome data of tubers and leaves from C. wenyujin were obtained and assembled into 78 092 unigenes. Of them, 244 unigenes were predicted to be involved in terpenoid biosynthesis while 131 unigenes were categorized as the "Terpenoid backbone biosynthesis" (TBB) term. Twenty-two unigenes possessed terpene synthase domain; five were predicted to be sesquiterpene synthases. Of the 208 unigenes annotated as cytochromes P450, 8 unigenes with full-length coding sequences were part of the CYP71 clade that primarily may perform hydroxylations of specialized metabolites. Furthermore, Ten DEGs related to the C5 precursor supply and sesquiterpene synthesis were validated by Real-time PCR; that showed a close correspondence with transcriptome sequence. A novel germacrene B synthase (CwGBS) and α-santalene synthase (CwSS) were identified in metabolically engineering E. coli. This study provided the first de novo transcriptome comparative analysis of leaf and tuber tissues from C. wenyujin, aiming to understand genetic mechanisms. Key genes involved in the biosynthesis of sesquiterpene will help for revealing the underlying mechanisms of C. wenyujin.

Gut Bacteria Associated With Monochamus saltuarius (Coleoptera: Cerambycidae) and Their Possible Roles in Host Plant Adaptations

Monochamus saltuarius (Coleoptera: Cerambycidae) is an important native pest in the pine forests of northeast China and a dispersing vector of an invasive species Bursaphelenchus xylophilus. To investigate the bacterial gut diversity of M. saltuarius larvae in different host species, and infer the role of symbiotic bacteria in host adaptation, we used 16S rRNA gene Illumina sequencing and liquid chromatography-mass spectrometry metabolomics processing to obtain and compare the composition of the bacterial community and metabolites in the midguts of larvae feeding on three host tree species: Pinus koraiensis, Pinus sylvestris var. mongolica, and Pinus tabuliformis. Metabolomics in xylem samples from the three aforementioned hosts were also performed. Proteobacteria and Firmicutes were the predominant bacterial phyla in the larval gut. At the genus level, Klebsiella, unclassified_f__Enterobacteriaceae, Lactococcus, and Burkholderia–Caballeronia–Paraburkholderia were most dominant in P. koraiensis and P. sylvestris var. mongolica feeders, while Burkholderia–Caballeronia–Paraburkholderia, Dyella, Pseudoxanthomonas, and Mycobacterium were most dominant in P. tabuliformis feeders. Bacterial communities were similar in diversity in P. koraiensis and P. sylvestris var. mongolica feeders, while communities were highly diverse in P. tabuliformis feeders. Compared with the other two tree species, P. tabuliformis xylems had more diverse and abundant secondary metabolites, while larvae feeding on these trees had a stronger metabolic capacity for secondary metabolites than the other two host feeders. Correlation analysis of the association of microorganisms with metabolic features showed that dominant bacterial genera in P. tabuliformis feeders were more negatively correlated with plant secondary metabolites than those of other host tree feeders.

Molecular and Functional Evolution of the Spermatophyte Sesquiterpene Synthases

Sesquiterpenes are important defense and signal molecules for plants to adapt to the environment, cope with stress, and communicate with the outside world, and their evolutionary history is closely related to physiological functions. In this study, the information of plant sesquiterpene synthases (STSs) with identified functions were collected and sorted to form a dataset containing about 500 members. The phylogeny of spermatophyte functional STSs was constructed based on the structural comparative analysis to reveal the sequence–structure–function relationships. We propose the evolutionary history of plant sesquiterpene skeletons, from chain structure to small rings, followed by large rings for the first time and put forward a more detailed function-driven hypothesis. Then, the evolutionary origins and history of spermatophyte STSs are also discussed. In addition, three newly identified STSs CaSTS2, CaSTS3, and CaSTS4 were analyzed in this functional evolutionary system, and their germacrene D products were consistent with the functional prediction. This demonstrates an application of the structure-based phylogeny in predicting STS function. This work will help us to understand evolutionary patterns and dynamics of plant sesquiterpenes and STSs and screen or design STSs with specific product profiles as functional elements for synthetic biology application.

Dynamics of Isomeric and Enantiomeric Fractions of Pinene in Essential Oil of Picea abies Annual Needles during Growing Season

Norway spruce (Picea abies (L.) H. Karst.) is one of the most important commercial tree species distributed naturally in the Boreal and subalpine forest zone of Europe. All parts of spruce trees, including needles, accumulate essential oils with a variety of chemical properties and ecological functions, such as modulating plant–insect communication. Annual needle samples from 15 trees (five from each of three habitats) of 15–17 years old were assayed for essential oils and their major compounds, including α-pinene, β-pinene, (1S)-(−)-α-pinene, and (1R)-(+)-α-pinene across a growing season. Results showed strong positive correlation between percentages of α- and β-pinene isomers (r = 0.69, p < 0.05) and between pinene isomers and essential oils: α-pinene correlated with essential oil stronger (r = 0.62, p < 0.05) than β-pinene (r = 0.33, p < 0.05). Correlation analyses performed with some weather conditions, including average monthly temperature, growing sum of effective temperatures over 5 °C, duration of sunshine, accumulated precipitation, relative humidity, and pressure, showed that temperature is the most important weather condition related to pinene dynamics: negative correlations of moderate strength were established between percentages of α- and β- pinenes and average monthly temperatures (r = −0.36, p < 0.01, n = 75 and r = −0.33, p < 0.01, n = 75, respectively). Out of pinene enantiomers, only (1S)-(−)-α-pinene showed some negative correlation with monthly temperature (r = −0.26, p < 0.05, n = 75). Different patterns of essential oil and pinene dynamics during growing season within separate habitats suggested that some genetic variables of Picea abies might be involved.

Discovery and characterization of diterpene synthases in Chamaecyparis formosensis Matsum. which participated in an unprecedented diterpenoid biosynthesis route in conifer

Chamaecyparis formosensis Matsum. is an endemic and precious coniferous species of Taiwan, and is known for a high abundance of specialized metabolites, which contributes to the excellent timber durability. Several terpenoids were identified and isolated from C. formosensis wood and needles, and exhibit anti-fungal and anti-bacterial bioactivities, which may participate in plant defense against pathogens. In various identified compounds, not only cadinene and ferruginol, were identified in C. formosensis extracts but also unique diterpenoids, which include pisferal, totarol, and derivates of isoabienol. To understand the biosynthesis of these specific diterpenoids, we conducted a series of functional characterization of the C. formosensis diterpene synthases (CfdiTPSs), which participate in skeleton formation and differentiation of diterpenes. In this study, we identified eight diTPSs from C. formosensis transcriptome, and they all contain either class I or class II motif, which indicates they are all monofunctional enzymes. These candidates consist of three class II diTPSs and five class I diTPSs, and after conducting in vivo and in vitro assays, class II diTPS CfCPS1 was characterized as a (+)-copalyl diphosphate synthase ((+)-CPS), and class I diTPSs CfKSL1 could further convert (+)-copalyl diphosphate ((+)-CPP) to levopimaradiene. Meanwhile, CfKSL1 also accepted labda-13-en-8-ol diphosphate (LPP) as substrate and formed monoyl oxide. Another class I diTPS, CfKSL4, exhibits a strong enzymatic ability of isoabienol synthase, which is firstly reported in conifer. This finding provides potential participants in the biosynthesis of unique diterpenoids, and with this knowledge, we can further expand our understanding of diterpenoid metabolism in Cupressaceae and their potential role in plant defense.

Functional characterization of an Indian sandalwood (Santalum album L.) dual-localized bifunctional nerolidol/linalool synthase gene involved in stress response

Essential oils extracted from the heartwood of Indian sandalwood (Santalum album L.) contain linalool and nerolidol as minor components. However, nerolidol/linalool synthase (NES/LIS), which produce linalool and nerolidol, have yet to be characterized in sandalwood. Using a transcriptomic-based approach, a terpene synthase gene was screened from unigenes of transcriptome data derived from S. album seedlings exposed to low temperature (4 °C). The enzyme encoded by these complementary DNAs belongs to the TPS-b clade. Recombinant SaNES/LIS is a bifunctional enzyme that can catalyze the formation of (E)-nerolidol from farnesyl diphosphate and linalool from geranyl diphosphate, respectively. Whereas SaNES/LIS was primarily localized in chloroplastids, both as granular fluorescence and as diffuse fluorescence, it was also detected in the cytosol of a limited number of cells. Agrobacterium tumefaciens-mediated transient gene expression in planta produced the same terpene products as those obtained in vitro. Real-time PCR analysis showed the highest expression of SaNES/LIS in fruits, with about a three-fold higher level than in leaves, followed by flowers, heartwood and roots. SaNES/LIS transcripts were differentially activated in different tissues in response to methyl jasmonate, cold, high temperature, strong illumination, and drought stress. Our results provide novel insight into the role of sandalwood terpenoids in response to various environmental stresses.

Biochemical Characterization and Function of Eight Microbial Type Terpene Synthases from Lycophyte Selaginella moellendorffii

Selaginella moellendorffii is a lycophyte, a member of an ancient vascular plant lineage. Two distinct types of terpene synthase (TPS) genes were identified from this species, including S. moellendorffii TPS genes (SmTPSs) and S. moellendorffii microbial TPS-like genes (SmMTPSLs). The goal of this study was to investigate the biochemical functions of SmMTPSLs. Here, eight full-length SmMTPSL genes (SmMTPSL5, -15, -19, -23, -33, -37, -46, and -47) were functionally characterized from S. moellendorffii. Escherichia coli-expressed recombinant SmMTPSLs were tested for monoterpenes synthase and sesquiterpenes synthase activities. These enzymatic products were typical monoterpenes and sesquiterpenes that have been previous shown to be generated by typical plant TPSs when provided with geranyl diphosphate (GPP) and farnesyl diphosphate (FPP) as the substrates. Meanwhile, SmMTPSL23, -33, and -37 were up-regulated when induced by alamethicin (ALA) and methyl jasmonate (MeJA), suggesting a role for these genes in plants response to abiotic stresses. Furthermore, this study pointed out that the terpenoids products of SmMTPSL23, -33, and -37 have an antibacterial effect on Pseudomonas syringae pv. tomato DC3000 and Staphylococcus aureus. Taken together, these results provide more information about the catalytic and biochemical function of SmMTPSLs in S. moellendorffii plants.

Two terpene synthases in resistant Pinus massoniana contribute to defence against Bursaphelenchus xylophilus

Pine wood nematode (PWN; Bursaphelenchus xylophilus), a destructive pest of Pinus massoniana, is causing a severe epidemic of pine wilt disease in China. When invaded by PWN, resistant P. massoniana secretes an abundance of oleoresin terpenoids as a defensive strategy. However, regulatory mechanisms of this defence in resistant P. massoniana have yet to be elucidated. Here, we characterized two terpene synthase genes, α-pinene synthase (PmTPS4) and longifolene synthase (PmTPS21), identified in resistant P. massoniana and investigate the contribution of these genes to the oleoresin defence strategy in resistant masson pines. Up-regulation of these two genes in the stem supported their involvement in terpene biosynthesis as part of the defence against PWN. Recombinant protein expression revealed catalytic activity for the two PmTPSs, with PmTPS4 primarily producing α-pinene, while PmTPS21 produced α-pinene and longifolene simultaneously. The major enzymatic products of the two terpene synthases had inhibitory effects on PWN in vitro. We demonstrated that PmTPS4 and PmTPS21 played positive roles in terpene-defence mechanisms against PWN infestation. The major products of these terpene synthases could directly inhibit the survival rate of PWN in vitro. We revealed that PmTPS21 was a novel bifunctional enzyme capable of simultaneous production of both monoterpene and sesquiterpene.

Identification and Functional Characterization of Tissue-Specific Terpene Synthases in Stevia rebaudiana

In addition to the well-known diterpenoid steviol glycosides, Stevia rebaudiana (Stevia) produces many labdane-type diterpenoids and a wide range of mono- and sesquiterpenoids. However, biosynthesis of mono- and sesquiterpenoids in Stevia remains unknown. Here we analyzed the extracts of Stevia leaves, flowers, stems, and roots by Gas Chromatography–Mass Spectrometry and putatively identified a total of 69 volatile organic compounds, most of which were terpenoids with considerably varied quantities among the four tissues of Stevia. Using Stevia transcriptomes, we identified and functionally characterized five terpene synthases (TPSs) that produced major mono- and sesquiterpenoids in Stevia. Transcript levels of these Stevia TPSs and levels of corresponding terpenoids correlated well in Stevia tissues. Particularly, the root-specific SrTPS4 and SrTPS5 catalyzed the formation of γ-curcumene/zingiberene/β-sesquiphellandrene and α-longipinene/β-himachalene/himachalol as multifunctional sesqui-TPSs, respectively. Most of the SrTPSs were highly responsive to various environmental stresses in a tissue-specific manner. Taken together, our results provide new insights into how Stevia produces diverse terpenoids to confer differential responses to various environmental factors in each tissue.

Production of multiple terpenes of different chain lengths by subcellular targeting of multi-substrate terpene synthase in plants

Multi-substrate terpene synthases (TPSs) are distinct from typical TPSs that react with a single substrate. Although in vitro activity of few multi-substrate TPSs have been reported, in vivo characterization has not been well investigated for most of them. Here, a new TPS from Cananga odorata, CoTPS5, belonging to TPS-f subfamily was functionally characterized in vitro as well as in vivo. CoTPS5 reacted with multiple prenyl-pyrophosphate substrates of various chain lengths as a multi-substrate TPS. It catalyzed the formation of (E)-β-ocimene, (E,E)-α-farnesene and α-springene from geranyl pyrophosphate, (E,E)-farnesyl pyrophosphate and geranylgeranyl pyrophosphate, respectively. Upon transient expression in Nicotiana benthamiana, CoTPS5 localized to cytosol and produced only (E,E)-α-farnesene. However, expression of plastid-targeted CoTPS5 in N. benthamiana resulted in biosynthesis of all three compounds, (E)-β-ocimene, (E,E)-α-farnesene and α-springene. Similarly, transgenic Arabidopsis plants overexpressing plastid-targeted CoTPS5 showed stable and sustainable production of (E)-β-ocimene, (E,E)-α-farnesene and α-springene. Moreover, their production did not affect the growth and development of transgenic Arabidopsis plants. Our results demonstrate that redirecting multi-substrate TPS to a different intracellular compartment could be an effective way to prove in vivo activity of multi-substrate TPSs and thereby allowing for the production of multiple terpenoids simultaneously in plants.

Evolution of Labdane-Related Diterpene Synthases in Cereals

Gibberellins (GAs) are labdane-related diterpenoid phytohormones that regulate various aspects of higher plant growth. A biosynthetic intermediate of GAs is ent-kaurene, a tetra-cyclic diterpene that is produced through successive cyclization of geranylgeranyl diphosphate catalyzed by the two distinct monofunctional diterpene synthases—ent-copalyl diphosphate synthase (ent-CPS) and ent-kaurene synthase (KS). Various homologous genes of the two diterpene synthases have been identified in cereals, including rice (Oryza sativa), wheat (Triticum aestivum) and maize (Zea mays), and are believed to have been derived from GA biosynthetic ent-CPS and KS genes through duplication and neofunctionalization. They play roles in specialized metabolism, giving rise to diverse labdane-related diterpenoids for defense because a variety of diterpene synthases generate diverse carbon-skeleton structures. This review mainly describes the diterpene synthase homologs that have been identified and characterized in rice, wheat and maize and shows the evolutionary history of various homologs in rice inferred by comparative genomics studies using wild rice species, such as Oryza rufipogon and Oryza brachyantha. In addition, we introduce labdane-related diterpene synthases in bryophytes and gymnosperms to illuminate the macroscopic evolutionary history of diterpene synthases in the plant kingdom—bifunctional enzymes possessing both CPS and KS activities are present in bryophytes; gymnosperms possess monofunctional CPS and KS responsible for GA biosynthesis and also possess bifunctional diterpene synthases facilitating specialized metabolism for defense.

A transcriptomic variation map provides insights into the genetic basis of Pinus massoniana Lamb. evolution and the association with oleoresin yield

BACKGROUND

Masson pine (Pinus massoniana Lamb.), the dominant native coniferous species in southern China, is commercially important for supplying timber and oleoresin. However, knowledge of the genetic variability of masson pine germplasm is still limited. In this study, the genetic diversity and population structure of masson pine germplasm were assessed using 204 wild accessions from 10 main distribution regions using 94,194 core single-nucleotide polymorphisms (SNPs) obtained from transcriptome sequencing data.

RESULTS

The average expected heterozygosity was 0.2724, implying abundant genetic diversity within masson pine germplasm. Analysis of molecular variance (AMOVA) revealed that 3.29% of the variation was sourced from genetic differentiation. Structure analysis identified two geographically distinct groups. Discriminant analysis of principal components (DAPC) showed that one of those groups was further divided into two clusters. Sichuan and Chongqing provenance is the geographical origin, which diffused outward along two different lines. Oleoresin yield is reflected in the evolution of the two groups, and exhibits two different trends along the two lines of diffusion. The oleoresin yield may be associated with the genes of chitinase, CYP720B, cytochrome P450, ABC transporter, and AP2/ethylene-responsive transcription factor (ERF) based on SNPs and expression.

CONCLUSIONS

SNP markers from transcriptome sequencing are highly capable of evaluating genetic diversity within different species, as well as the genetic control of objective traits. The functions of these genes will be verified in future studies, and those genes strongly associated with oleoresin yield will be used to improve yields by means of early genotype selection and genetic engineering.

Identification of Three Monofunctional Diterpene Synthases with Specific Enzyme Activities Expressed during Heartwood Formation in Western Redcedar (Thuja plicata) Trees

Upon harvest, Western redcedar (WRC; Thuja plicata) trees have a high incidence and extent of heartwood rot. While monoterpenoids and lignans have been linked to rot resistance in this species, other specialized metabolites, such as diterpenes, are likely to contribute to rot resistance. Here we report the cloning and functional assessment of three putative diterpene synthase (TpdiTPS) genes expressed during heartwood formation in WRC. The predicted proteins of the three genes lack either of the two catalytically independent active sites typical of most diTPS, indicating monofunctional rather than bifunctional activity. To identify potential catalytic activities of these proteins, we expressed them in genetically engineered Escherichia coli strains that produce four potential substrates, geranylgeranyl diphosphate (GGDP), ent, syn, and normal stereoisomers of copalyl diphosphate (CDP). We found that TpdiTPS3 used GGDP to produce CDP. TpdiTPS2 used normal CDP to produce levopimaradiene. TpdiTPS1 showed stereoselectivity as it used normal CDP to produce sandaracopimaradiene and syn-CDP to produce syn-stemod-13(17)-ene. These genes and protein enzymatic activities have not been previously reported in WRC and provide an opportunity to assess their potential roles in heartwood rot resistance in this economically important species.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

Exploring novel bacterial terpene synthases

Terpenes are the largest class of natural products with extensive structural diversity and are widely used as pharmaceuticals, herbicides, flavourings, fragrances, and biofuels. While they have mostly been isolated from plants and fungi, the availability and analysis of bacterial genome sequence data indicates that bacteria also possess many putative terpene synthase genes. In this study, we further explore this potential for terpene synthase activity in bacteria. Twenty two potential class I terpene synthase genes (TSs) were selected to represent the full sequence diversity of bacterial synthase candidates and recombinantly expressed in E. coli. Terpene synthase activity was detected for 15 of these enzymes, and included mono-, sesqui- and diterpene synthase activities. A number of confirmed sesquiterpene synthases also exhibited promiscuous monoterpene synthase activity, suggesting that bacteria are potentially a richer source of monoterpene synthase activity then previously assumed. Several terpenoid products not previously detected in bacteria were identified, including aromandendrene, acora-3,7(14)-diene and longiborneol. Overall, we have identified promiscuous terpene synthases in bacteria and demonstrated that terpene synthases with substrate promiscuity are widely distributed in nature, forming a rich resource for engineering terpene biosynthetic pathways for biotechnology.

A Comprehensive Survey on the Terpene Synthase Gene Family Provides New Insight into Its Evolutionary Patterns

Terpenes are organic compounds and play important roles in plant growth and development as well as in mediating interactions of plants with the environment. Terpene synthases (TPSs) are the key enzymes responsible for the biosynthesis of terpenes. Although some species were employed for the genome-wide identification and characterization of the TPS family, limited information is available regarding the evolution, expansion, and retention mechanisms occurring in this gene family. We performed a genome-wide identification of the TPS family members in 50 sequenced genomes. Additionally, we also characterized the TPS family from aromatic spearmint and basil plants using RNA-Seq data. No TPSs were identified in algae genomes but the remaining plant species encoded various numbers of the family members ranging from 2 to 79 full-length TPSs. Some species showed lineage-specific expansion of certain subfamilies, which might have contributed toward species or ecotype divergence or environmental adaptation. A large-scale family expansion was observed mainly in dicot and monocot plants, which was accompanied by frequent domain loss. Both tandem and segmental duplication significantly contributed toward family expansion and expression divergence and played important roles in the survival of these expanded genes. Our data provide new insight into the TPS family expansion and evolution and suggest that TPSs might have originated from isoprenyl diphosphate synthase genes.

On the Evolution and Functional Diversity of Terpene Synthases in the Pinus Species: A Review

In the biosynthesis of terpenoids, the ample catalytic versatility of terpene synthases (TPS) allows the formation of thousands of different molecules. A steadily increasing number of sequenced plant genomes invariably show that the TPS gene family is medium to large in size, comprising from 30 to 100 functional members. In conifers, TPSs belonging to the gymnosperm-specific TPS-d subfamily produce a complex mixture of mono-, sesqui-, and diterpenoid specialized metabolites, which are found in volatile emissions and oleoresin secretions. Such substances are involved in the defence against pathogens and herbivores and can help to protect against abiotic stress. Oleoresin terpenoids can be also profitably used in a number of different fields, from traditional and modern medicine to fine chemicals, fragrances, and flavours, and, in the last years, in biorefinery too. In the present work, after summarizing the current views on the biosynthesis and biological functions of terpenoids, recent advances on the evolution and functional diversification of plant TPSs are reviewed, with a focus on gymnosperms. In such context, an extensive characterization and phylogeny of all the known TPSs from different Pinus species is reported, which, for such genus, can be seen as the first effort to explore the evolutionary history of the large family of TPS genes involved in specialized metabolism. Finally, an approach is described in which the phylogeny of TPSs in Pinus spp. has been exploited to isolate for the first time mono-TPS sequences from Pinus nigra subsp. laricio, an ecologically important endemic pine in the Mediterranean area.

Priming of inducible defenses protects Norway spruce against tree-killing bark beetles

Plants can form an immunological memory known as defense priming, whereby exposure to a priming stimulus enables quicker or stronger response to subsequent attack by pests and pathogens. Such priming of inducible defenses provides increased protection and reduces allocation costs of defense. Defense priming has been widely studied for short-lived model plants such as Arabidopsis, but little is known about this phenomenon in long-lived plants like spruce. We compared the effects of pretreatment with sublethal fungal inoculations or application of the phytohormone methyl jasmonate (MeJA) on the resistance of 48-year-old Norway spruce (Picea abies) trees to mass attack by a tree-killing bark beetle beginning 35 days later. Bark beetles heavily infested and killed untreated trees but largely avoided fungus-inoculated trees and MeJA-treated trees. Quantification of defensive terpenes at the time of bark beetle attack showed fungal inoculation induced 91-fold higher terpene concentrations compared with untreated trees, whereas application of MeJA did not significantly increase terpenes. These results indicate that resistance in fungus-inoculated trees is a result of direct induction of defenses, whereas resistance in MeJA-treated trees is due to defense priming. This work extends our knowledge of defense priming from model plants to an ecologically important tree species.

Full-Length Transcriptome Sequencing and Different Chemotype Expression Profile Analysis of Genes Related to Monoterpenoid Biosynthesis in Cinnamomum porrectum

Leaves of C. porrectum are rich in essential oils containing monoterpenes, sesquiterpenes and aromatic compounds, but the molecular mechanism of terpenoid biosynthesis in C. porrectum is still unclear. In this paper, the differences in the contents and compositions of terpenoids among three chemotypes were analyzed using gas chromatography mass spectrometry (GC/MS). Furthermore, the differential expression of gene transcripts in the leaf tissues of the three C. porrectum chemotypes were analyzed through a comparison of full-length transcriptomes and expression profiles. The essential oil of the three C. porrectum chemotypes leaves was mainly composed of monoterpenes. In the full-length transcriptome of C. porrectum, 104,062 transcripts with 306,337,921 total bp, an average length of 2944 bp, and an N50 length of 5449 bp, were obtained and 94025 transcripts were annotated. In the eucalyptol and linalool chemotype, the camphor and eucalyptol chemotype, and the camphor and linalool chemotype comparison groups, 21, 22 and 18 terpene synthase (TPS) unigenes were identified respectively. Three monoterpene synthase genes, CpTPS3, CpTPS5 and CpTPS9, were upregulated in the eucalyptol chemotype compared to the linalool chemotype and camphor chemotype. CpTPS1 was upregulated in the camphor chemotype compared to the linalool chemotype and the eucalyptol chemotype. CpTPS4 was upregulated in the linalool chemotype compared to the camphor chemotype and the eucalyptol chemotype. Different unigenes had different expression levels among the three chemotypes, but the unigene expression levels of the 2-C-methyl-D-erythritol 4phosphate (MEP) pathway were generally higher than those of the mevalonate acid (MVA) pathway. Quantitative reverse transcription PCR(qRT-PCR) further validated these expression levels. The present study provides new clues for the functional exploration of the terpenoid synthesis mechanism and key genes in different chemotypes of C. porrectum.

Phylogenetically distant group of terpene synthases participates in cadinene and cedrane-type sesquiterpenes accumulation in Taiwania cryptomerioides

Along with the species evolution, plants have evolved ways to produce a different collection of terpenoids to accommodate its biotic and abiotic environment, and terpene synthase (TPS) is one of the major contributors to various terpene compounds. The timber of a monotypic and relictual conifer species of Cupressace, Taiwania cryptomerioides, has excellent durability, and one of the essential factors for Taiwania to resist decay and insect pests is sesquiterpene. Compared to other conifers, Taiwania has much higher abundance of cadinene-type sesquiterpenes, and the presence of cedrene-type sesquiterpenes. To understand sesquiterpene biosynthesis in Taiwania, we functionally characterized 10 T. cryptomerioides TPSs (TcTPSs) in vivo or in planta, which could catalyze sesquiterpene formation and potentially are involved in biosynthesis of diverse sesquiterpenoids in Taiwania. The distant phylogenetic relationship and the intron loss event of TcTPSs correlate to the differentiation of chemical profile Taiwania compared to other conifers. Furthermore, we identified TcTPS3 and TcTPS12 as δ-cadinene synthase, and TcTPS6 as cedrol synthase, which demonstrates the important contributions of dynamic evolution in TPSs to the chemical diversity in plants. Combining with functional characterization and comparison of catalytic residues, we conclude at least three catalytic routes for sesquiterpene biosynthesis in this species, and the skeleton diversity has been expended in T. cryptomeriodes.