The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom

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.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Origin and early evolution of the plant terpene synthase family

As a midsized gene family conserved more by lineage than function, the typical plant terpene synthases (TPSs) could be a valuable tool to examine plant evolution. TPSs are pivotal in biosynthesis of gibberellins and related phytohormones as well as in formation of the extensive arsenal of specialized plant metabolites mediating ecological interactions whose production is often lineage specific. Yet the origin and early evolution of the TPS family is not well understood. Systematic analysis of an array of transcriptomes and sequenced genomes indicated that the TPS family originated after the divergence of land plants from charophytic algae. Phylogenetic and biochemical analyses support the hypothesis that the ancestral TPS gene encoded a bifunctional class I and II diterpene synthase producing the ent-kaurene required for phytohormone production in all extant lineages of land plants. Moreover, the ancestral TPS gene likely underwent duplication at least twice early in land plant evolution. Together these two gave rise to three TPS lineages leading to the extant TPS-c, TPS-e/f, and the remaining TPS (h/d/a/b/g) subfamilies, with the latter dedicated to secondary rather than primary metabolism while the former two contain those genes involved in ent-kaurene production. Nevertheless, parallel evolution from the ent-kaurene–producing class I and class II diterpene synthases has led to roles for TPS-e/f and -c subfamily members in secondary metabolism as well. These results clarify TPS evolutionary history and provide context for the role of these genes in producing the vast diversity of terpenoid natural products observed today in various land plant lineages.

Genome-wide association study identifies SNP markers and putative candidate genes for terpene traits important for Leptocybe invasa resistance in Eucalyptus grandis

Terpenes are an important group of plant specialized metabolites influencing, amongst other functions, defence mechanisms against pests. We used a genome-wide association study to identify single nucleotide polymorphism (SNP) markers and putative candidate genes for terpene traits. We tested 15,387 informative SNP markers derived from genotyping 416 Eucalyptus grandis individuals for association with 3 terpene traits, 1,8-cineole, γ-terpinene, and p-cymene. A multilocus mixed model analysis identified 21 SNP markers for 1,8-cineole on chromosomes 2, 4, 6, 7, 8, 9, 10, and 11, that individually explained 3.0%–8.4% and jointly 42.7% of the phenotypic variation. Association analysis of γ-terpinene found 32 significant SNP markers on chromosomes 1, 2, 4, 5, 6, 9, and 11, explaining 3.4–15.5% and jointly 54.5% of phenotypic variation. For p-cymene, 28 significant SNP markers were identified on chromosomes 1, 2, 3, 5, 6, 7, 10, and 11, explaining 3.4–16.1% of the phenotypic variation and jointly 46.9%. Our results show that variation underlying the 3 terpene traits is influenced by a few minor loci in combination with a few major effect loci, suggesting an oligogenic nature of the traits.

Deep roots and many branches: Origins of plant-specialized metabolic enzymes in general metabolism

Collectively, plants produce hundreds of thousands of specialized metabolites from simple building blocks such as amino acids, fatty acids, and isoprenoids. As additional specialized metabolic enzymes are described, there is increasing recognition of the importance of cooption of general metabolic enzymes to specialized metabolism by gene duplication, narrowing of expression, and alteration of enzymatic activities. Here, we examine how several classes of enzymes were each coopted multiple times. We demonstrate the simplicity of achieving the synthesis of analogous chemicals by coopting existing enzymes and summarize emerging insights that could inform rational metabolic engineering of both general and specialized metabolic enzymes.

Elucidation of the essential oil biosynthetic pathways in Cinnamomum burmannii through identification of six terpene synthases

Cinnamomum burmannii is a traditional plant that has long been used as a spice, food preservative, and food flavoring. Essential oils in C. burmannii, which mainly consist of mono- and sesquiterpenes such borneol, linalool, and caryophyllene, have impressive pharmaceutical properties. Although the transcriptome-based discovery of (+)-bornyl diphosphate synthase (CbTPS1) from C. burmannii was reported in our previous study, the remaining terpene synthases (TPSs) corresponding to various terpene biosynthesis pathways remain unidentified. In this study, we report the results of RNA-sequencing of a borneol type plant and functional characterization of six additional full-length candidate TPS genes (named CbTPS2-7). Phylogenetic analysis revealed that CbTPS2 and CbTPS3 together with the previously identified CbTPS1 protein belong to the TPS-b subfamily, and enzyme assays using geranyl diphosphate (GPP) and farnesyl diphosphate (FPP) as substrates revealed that CbTPS1, CbTPS2 and CbTPS3 catalyze the formation of monoterpenes. CbTPS4, CbTPS5, and CbTPS6, which belong to the TPS-a clade, generated monoterpenes and sesquiterpenes. CbTPS7, which belongs to the TPS-g clade, showed linalool/nerolidol synthase activity. These CbTPSs identified in C. burmannii produced a total of 10 monoterpenes and 14 sesquiterpenes in an in vitro assay. These findings clarify the biosynthesis pathways of 13 monoterpenoids and 12 sesquiterpenoids in the leaf essential oil of C. burmannii and shed light on terpene biosynthesis in Cinnamomum.

The Cycas genome and the early evolution of seed plants

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

Phytocompounds as an Alternative Antimicrobial Approach in Aquaculture

Despite culturing the fastest-growing animal in animal husbandry, fish farmers are often adversely economically affected by pathogenic disease outbreaks across the world. Although there are available solutions such as the application of antibiotics to mitigate this phenomenon, the excessive and injudicious use of antibiotics has brought with it major concerns to the community at large, mainly due to the rapid development of resistant bacteria. At present, the use of natural compounds such as phytocompounds that can be an alternative to antibiotics is being explored to address the issue of antimicrobial resistance (AMR). These phytocompounds are bioactive agents that can be found in many species of plants and hold much potential. In this review, we will discuss phytocompounds extracted from plants that have been evidenced to contain antimicrobial, antifungal, antiviral and antiparasitic activities. Further, it has also been found that compounds such as terpenes, phenolics, saponins and alkaloids can be beneficial to the aquaculture industry when applied. This review will focus mainly on compounds that have been identified between 2000 and 2021. It is hoped this review will shed light on promising phytocompounds that can potentially and effectively mitigate AMR.

Tissue-Specific Expression of the Terpene Synthase Family Genes in Rosa chinensis and Effect of Abiotic Stress Conditions

Rose (Rosa chinensis) is one of the most famous ornamental plants worldwide, with a variety of colors and fragrances. Terpene synthases (TPSs) play critical roles in the biosynthesis of terpenes. In this work, we report a comprehensive study on the genome-wide identification and characterization of the TPS family in R. chinensis. We identified 49 TPS genes in the R. chinensis genome, and they were grouped into five subfamilies (TPS-a, TPS-b, TPS-c, TPS-g and TPS-e/f). Phylogenetics, gene structure and conserved motif analyses indicated that the RcTPS genes possessed relatively conserved gene structures and the RcTPS proteins contained relatively conserved motifs. Multiple putative cis-acting elements involved in the stress response were identified in the promoter region of RcTPS genes, suggesting that some could be regulated by stress. The expression profile of RcTPS genes showed that they were predominantly expressed in the petals of open flowers, pistils, leaves and roots. Under osmotic and heat stresses, the expression of most RcTPS genes was upregulated. These data provide a useful foundation for deciphering the functional roles of RcTPS genes during plant growth as well as addressing the link between terpene biosynthesis and abiotic stress responses in roses.

Floral Scents and Fruit Aromas: Functions, Compositions, Biosynthesis, and Regulation

Floral scents and fruit aromas are crucial volatile organic compounds (VOCs) in plants. They are used in defense mechanisms, along with mechanisms to attract pollinators and seed dispersers. In addition, they are economically important for the quality of crops, as well as quality in the perfume, cosmetics, food, drink, and pharmaceutical industries. Floral scents and fruit aromas share many volatile organic compounds in flowers and fruits. Volatile compounds are classified as terpenoids, phenylpropanoids/benzenoids, fatty acid derivatives, and amino acid derivatives. Many genes and transcription factors regulating the synthesis of volatiles have been discovered. In this review, we summarize recent progress in volatile function, composition, biosynthetic pathway, and metabolism regulation. We also discuss unresolved issues and research perspectives, providing insight into improvements and applications of plant VOCs.

Diversity and Functional Evolution of Terpene Synthases in Rosaceae

Terpenes are organic compounds and play important roles in plant development and stress response. Terpene synthases (TPSs) are the key enzymes for the biosynthesis of terpenes. For Rosaceae species, terpene composition represents a critical quality attribute, but limited information is available regarding the evolution and expansion occurring in the terpene synthases gene family. Here, we selected eight Rosaceae species with sequenced and annotated genomes for the identification of TPSs, including three Prunoideae, three Maloideae, and two Rosoideae species. Our data showed that the TPS gene family in the Rosaceae species displayed a diversity of family numbers and functions among different subfamilies. Lineage and species-specific expansion of the TPSs accompanied by frequent domain loss was widely observed within different TPS clades, which might have contributed to speciation or environmental adaptation in Rosaceae. In contrast to Maloideae and Rosoideae species, Prunoideae species owned less TPSs, with the evolution of Prunoideae species, TPSs were expanded in modern peach. Both tandem and segmental duplication significantly contributed to TPSs expansion. Ka/Ks calculations revealed that TPSs genes mainly evolved under purifying selection except for several pairs, where the divergent time indicated TPS-e clade was diverged relatively anciently. Gene function classification of TPSs further demonstrated the function diversity among clades and species. Moreover, based on already published RNA-Seq data from NCBI, the expression of most TPSs in Malus domestica, Prunus persica, and Fragaria vesca displayed tissue specificity and distinct expression patterns either in tissues or expression abundance between species and TPS clades. Certain putative TPS-like proteins lacking both domains were detected to be highly expressed, indicating the underlying functional or regulatory potentials. The result provided insight into the TPS family evolution and genetic information that would help to improve Rosaceae species quality.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases

Monoterpene synthases are often promiscuous enzymes, yielding product mixtures rather than pure compounds due to the nature of the branched reaction mechanism involving reactive carbocations. Two previously identified bacterial monoterpene synthases, a linalool synthase (bLinS) and a cineole synthase (bCinS), produce nearly pure linalool and cineole from geranyl diphosphate, respectively. We used a combined experimental and computational approach to identify critical residues involved in bacterial monoterpenoid synthesis. Phe77 is essential for bCinS activity, guiding the linear carbocation intermediate towards the formation of the cyclic α-terpinyl intermediate; removal of the aromatic ring results in variants that produce acyclic products only. Computational chemistry confirmed the importance of Phe77 in carbocation stabilisation. Phe74, Phe78 and Phe179 are involved in maintaining the active site shape in bCinS without a specific role for the aromatic ring. Phe295 in bLinS, and the equivalent Ala301 in bCinS, are essential for linalool and cineole formation, respectively. Where Phe295 places steric constraints on the carbocation intermediates, Ala301 is essential for bCinS initial cyclisation and activity. Our multidisciplinary approach gives unique insights into how carefully placed amino acid residues in the active site can direct carbocations down specific paths, by placing steric constraints or offering stabilisation via cation-π interactions.

Oxylipin biosynthesis in spikemoss Selaginella moellendorffii: Identification of allene oxide synthase (CYP74L2) and hydroperoxide lyase (CYP74L1)

Nonclassical P450s of the CYP74 family catalyse the secondary conversions of fatty acid hydroperoxides to bioactive oxylipins in plants. The model organism, spikemoss Selaginella moellendorffii Hieron, possesses at least ten CYP74 genes of novel J, K, L, and M subfamilies. The cloning of three CYP74L genes and catalytic properties of recombinant proteins are described in the present work. The CYP74L1 possessed mainly hydroperoxide lyase (HPL) activity towards the 13(S)-hydroperoxide of α-linolenic acids (13-HPOT) and nearly equal HPL and allene oxide synthase (AOS) activities towards the 13(S)-hydroperoxide of linoleic acids (13-HPOD). The 9-hydroperoxides were poor substrates for CYP74L1 and led to the production of mainly the α-ketols (AOS products) and minorities of HPL and epoxyalcohol synthase (EAS) products. The CYP74L2 possessed the AOS activity towards all tested hydroperoxides. CYP74L3 possessed low HPL/EAS activity. Besides, the aerial parts of S. moellendorffii plants possessed complex oxylipins patterns including divinyl ethers, epoxyalcohols, and 12-oxo-phytodienoic acid. Characterization of the CYP74L enzymes and oxylipin pattern updates the knowledge on the complex oxylipin biosynthetic machinery in the surviving oldest taxa of vascular plants.

Full-length transcriptome analysis of multiple organs and identification of adaptive genes and pathways in Mikania micrantha

Mikania micrantha is a notorious invasive weed that has caused huge economic loss and negative ecological consequences in invaded areas. This species can adapt well to invasive environments with various stress factors. The identification of gene families and functional pathways related to environmental adaptability is lack in M. micrantha at the multi-organ full-length transcriptome level. In this study, we sequenced the transcriptomes of five M. micrantha organs using PacBio single-molecule real-time sequencing and Illumina RNA sequencing technologies. Based on the transcriptome data, full-length transcripts were captured and gene expression patterns among the five organs were analyzed. KEGG enrichment analysis of genes with higher expression indicated their special roles in environmental stress response and adversity adaptation in the various five organs. The gene families and pathways related to biotic and abiotic factors, including terpene synthases, glutathione S-transferases, antioxidant defense system, and terpenoid biosynthesis pathway, were characterized. The expression levels of most differentially expressed genes in the antioxidant defense system and terpenoid biosynthesis pathway were higher in root, stem, and leaf than in the other two organs, suggesting that root, stem, and leaf have strong ability to respond to adverse stresses and form the important organs of terpenoid synthesis and accumulation. Additionally, a large number of transcription factors and alternative splicing events were predicted. This study provides a comprehensive transcriptome resource for M. micrantha, and our findings facilitate further research on the adaptive evolution and functional genomics of this species.

A chromosome-scale genome and transcriptomic analysis of the endangered tropical tree Vatica mangachapoi (Dipterocarpaceae)

Vatica mangachapoi is a tropical tree species native to Southeast Asia. It has long been valued as a timber species because the wood resists decay, but it is now considered vulnerable to extinction due to habitat loss and overexploitation. Here, we present the first chromosome-level genome assembly of V. mangachapoi that we created by combining data from PacBio long read sequencing with Hi-C proximity ligation and Illumina short-read sequencing. The assembled genome was 456.21 Mb, containing 11 chromosome and a BUSCO score of 93.4%. From the newly assembled genome, 46,811 protein-coding genes were predicted. Repetitive DNA accounted for 53% of the genome. Phylogenomic and gene family analyses showed that V. mangachapoi diverged from a common ancestor of Gossypium raimondii 70 million years ago. Transcriptome analyses found 227 genes that were differentially expressed in the leaves of plants grown in normal soil relative to plants grown in dry, coastal, sandy soil. For these genes, we identified three significantly enriched with GO terms: responses to organonitrogen compounds, chitin-triggered immunity, and wound response. This genome provides an important comparative benchmark not only for future conservation work on V. mangachapoi but also for phylogenomics work on Dipterocarpaceae.

Prunus persica Terpene Synthase PpTPS1 Interacts with PpABI5 to Enhance Salt Resistance in Transgenic Tomatoes

Terpene synthase (TPS) is related to the production of aromatic substances, but there are few studies on the impact of abiotic stress on TPS and its molecular mechanism, especially in peaches. This study found that salt resistance and abscisic acid (ABA) sensitivity of transgenic tomatoes were enhanced by overexpression of PpTPS1. Moreover, it was found that PpTPS1 interacted with and antagonized the expression of the bZIP transcription factor ABA INSENSITIVE 5 (PpABI5), which is thought to play an important role in salt suitability. In addition, PpTCP1, PpTCP13, and PpTCP15 were found to activate the expression of PpTPS1 by yeast one-hybrid (Y1H) and dual-luciferase assays, and they could also be induced by ABA. In summary, PpTPS1 may be involved in the ABA signaling regulatory pathway and play an important role in salt acclimation, providing a new reference gene for the improvement of salt resistance in peaches.

Differential gene expression associated with a floral scent polymorphism in the evening primrose Oenothera harringtonii (Onagraceae)

Background

Plant volatiles play an important role in both plant-pollinator and plant-herbivore interactions. Intraspecific polymorphisms in volatile production are ubiquitous, but studies that explore underlying differential gene expression are rare. Oenothera harringtonii populations are polymorphic in floral emission of the monoterpene (R)-(−)-linalool; some plants emit (R)-(−)-linalool (linalool+ plants) while others do not (linalool- plants). However, the genes associated with differential production of this floral volatile in Oenothera are unknown. We used RNA-Seq to broadly characterize differential gene expression involved in (R)-(−)-linalool biosynthesis. To identify genes that may be associated with the polymorphism for this trait, we used RNA-Seq to compare gene expression in six different Oenothera harringtonii tissues from each of three linalool+ and linalool- plants.

Results

Three clusters of differentially expressed genes were enriched for terpene synthase activity: two were characterized by tissue-specific upregulation and one by upregulation only in plants with flowers that produce (R)-(−)-linalool. A molecular phylogeny of all terpene synthases identified two putative (R)-(−)-linalool synthase transcripts in Oenothera harringtonii, a single allele of which is found exclusively in linalool+ plants.

Conclusions

By using a naturally occurring polymorphism and comparing different tissues, we were able to identify candidate genes putatively involved in the biosynthesis of (R)-(−)-linalool. Expression of these genes in linalool- plants, while low, suggests a regulatory polymorphism, rather than a population-specific loss-of-function allele. Additional terpene biosynthesis-related genes that are up-regulated in plants that emit (R)-(−)-linalool may be associated with herbivore defense, suggesting a potential economy of scale between plant reproduction and defense.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08370-6.

The genome of the Paleogene relic tree Bretschneidera sinensis: insights into trade-offs in gene family evolution, demographic history, and adaptive SNPs

Among relic species, genomic information may provide the key to inferring their long-term survival. Therefore, in this study, we investigated the genome of the Paleogene relic tree species, Bretschneidera sinensis, which is a rare endemic species within southeastern Asia. Specifically, we assembled a high-quality genome for B. sinensis using PacBio high-fidelity and high-throughput chromosome conformation capture reads and annotated it with long and short RNA sequencing reads. Using the genome, we then detected a trade-off between active and passive disease defences among the gene families. Gene families involved in salicylic acid and MAPK signalling pathways expanded as active defence mechanisms against disease, but families involved in terpene synthase activity as passive defences contracted. When inferring the long evolutionary history of B. sinensis, we detected population declines corresponding to historical climate change around the Eocene–Oligocene transition and to climatic fluctuations in the Quaternary. Additionally, based on this genome, we identified 388 single nucleotide polymorphisms (SNPs) that were likely under selection, and showed diverse functions in growth and stress responses. Among them, we further found 41 climate-associated SNPs. The genome of B. sinensis and the SNP dataset will be important resources for understanding extinction/diversification processes using comparative genomics in different lineages.

Toward improved terpenoids biosynthesis: strategies to enhance the capabilities of cell factories

Terpenoids form the most diversified class of natural products, which have gained application in the pharmaceutical, food, transportation, and fine and bulk chemical industries. Extraction from naturally occurring sources does not meet industrial demands, whereas chemical synthesis is often associated with poor enantio-selectivity, harsh working conditions, and environmental pollutions. Microbial cell factories come as a suitable replacement. However, designing efficient microbial platforms for isoprenoid synthesis is often a challenging task. This has to do with the cytotoxic effects of pathway intermediates and some end products, instability of expressed pathways, as well as high enzyme promiscuity. Also, the low enzymatic activity of some terpene synthases and prenyltransferases, and the lack of an efficient throughput system to screen improved high-performing strains are bottlenecks in strain development. Metabolic engineering and synthetic biology seek to overcome these issues through the provision of effective synthetic tools. This review sought to provide an in-depth description of novel strategies for improving cell factory performance. We focused on improving transcriptional and translational efficiencies through static and dynamic regulatory elements, enzyme engineering and high-throughput screening strategies, cellular function enhancement through chromosomal integration, metabolite tolerance, and modularization of pathways.

Characterization of Terpene synthase variation in flowers of wild aquilegia species from Northeastern Asia

There are several causes for the great diversity in floral terpenes. The terpene products are determined by the catalytic fidelity, efficiency and plasticity of the active sites of terpene synthases (TPSs). However, the molecular mechanism of TPS in catalyzing terpene biosynthesis and its evolutionary fate in wild plant species remain largely unknown. In this study, the functionality of terpene synthases and their natural variants were assessed in two Northeastern Asia endemic columbine species and their natural hybrid. Synoptically, TPS7, TPS8, and TPS9 were highly expressed in these Aquilegia species from the Zuojia population. The in vitro and in vivo enzymatic assays revealed that TPS7 and TPS8 mainly produced (+)-limonene and β-sesquiphellandrene, respectively, whereas TPS9 produced pinene, similar to the major components released from Aquilegia flowers. Multiple sequence alignment of Aquilegia TPS7 and TPS8 in the Zuojia population revealed amino acid polymorphisms. Domain swapping and amino acid substitution assays demonstrated that 413A, 503I and 529D had impacts on TPS7 catalytic activity, whereas 420G, 538F and 545 L affected the ratio of β-sesquiphellandrene to β-bisabolene in TPS8. Moreover, these key polymorphic amino acid residues were found in Aquilegia species from the Changbai Mountain population. Interestingly, amino acid polymorphisms in TPSs were present in individuals with low expression levels, and nonsynonymous mutations could impact the catalytic activity or product specificity of these genes. The results of this study will shed new light on the function and evolution of TPS genes in wild plant species and are beneficial to the modification of plant fragrances.

The reference genome of camellia chekiangoleosa provides insights into camellia evolution and tea oil biosynthesis

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.

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.

MrTPS3 and MrTPS20 Are Responsible for β-Caryophyllene and α-Pinene Production, Respectively, in Red Bayberry (Morella rubra)

Red bayberry is a sweet, tart fruit native to China and grown widely in the south. The key organic compounds forming the distinctive aroma in red bayberry, are terpenoids, mainly β-caryophyllene and α-pinene. However, the key genes responsible for different terpenoids are still unknown. Here, transcriptome analysis on samples from four cultivars, during fruit development, with different terpenoid production, provided candidate genes for volatile organic compound (VOC) production. Terpene synthases (TPS) are key enzymes regulating terpenoid biosynthesis, and 34 TPS family members were identified in the red bayberry genome. MrTPS3 in chromosome 2 and MrTPS20 in chromosome 7 were identified as key genes regulating β-caryophyllene and α-pinene synthesis, respectively, by qRT-PCR. Subcellular localization and enzyme activity assay showed that MrTPS3 was responsible for β-caryophyllene (sesquiterpenes) production and MrTPS20 for α-pinene (monoterpenes). Notably, one amino acid substitution between dark color cultivars and light color cultivars resulted in the loss of function of MrTPS3, causing the different β-caryophyllene production. Our results lay the foundation to study volatile organic compounds (VOCs) in red bayberry and provide potential genes for molecular breeding.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

Identification of (-)-bornyl diphosphate synthase from Blumea balsamifera and its application for (-)-borneol biosynthesis in Saccharomyces cerevisiae

Borneol is a precious monoterpenoid with two chiral structures, (-)-borneol and (+)-borneol. Bornyl diphosphate synthase is the key enzyme in the borneol biosynthesis pathway. Many (+)-bornyl diphosphate synthases have been reported, but no (-)-bornyl diphosphate synthases have been identified. Blumea balsamifera leaves are rich in borneol, almost all of which is (-)-borneol. In this study, we identified a high-efficiency (-)-bornyl diphosphate synthase (BbTPS3) from B. balsamifera that converts geranyl diphosphate (GPP) to (-)-bornyl diphosphate, which is then converted to (-)-borneol after dephosphorylation in vitro. BbTPS3 exhibited a K_m value of 4.93 ± 1.38 μM for GPP, and the corresponding k_cat value was 1.49 s⁻¹. Multiple strategies were applied to obtain a high-yielding (-)-borneol producing yeast strain. A codon-optimized BbTPS3 protein was introduced into the GPP high-yield strain MD, and the resulting MD-B1 strain produced 1.24 mg·L^-1 (-)-borneol. After truncating the N-terminus of BbTPS3 and adding a Kozak sequence, the (-)-borneol yield was further improved by 4-fold to 4.87 mg·L^-1. Moreover, the (-)-borneol yield was improved by expressing the fusion protein module of ERG20^F96W-N127W-YRSQI-t14-BbTPS3K2, resulting in a final yield of 12.68 mg·L^-1 in shake flasks and 148.59 mg·L^-1 in a 5-L bioreactor. This work is the first reported attempt to produce (-)-borneol by microbial fermentation.

Cloning and Characterization of 1,8-Cineole Synthase (SgCINS) Gene From the Leaves of Salvia guaranitica Plant

Monoterpenes are one of the most common groups belonging to the terpenoid family, with a C10 structure comprising of two isoprene units. Most of monoterpenes are volatile plant compounds, and they act as signaling molecules between plants and the environment, particularly as defensive compounds against herbivores and pathogens. In this study, 1,8-cineole synthase (SgCINS) gene was identified and cloned from the leaves of Salvia guaranitica plant. To examine the role of SgCINS in insect resistance, we transformed and expressed this gene into tobacco leaves. The metabolic analysis revealed that the production of various types and amount of terpenoid was increased and decreased in SgCINS overexpression and control lines, respectively, suggesting that overexpressing SgCINS in transgenic tobacco plants lead to an increase in the production of various types of terpenoids and other phytochemical compounds. These results indicated why transgenic tobacco was highly resistant against cotton worm than the highly susceptible control plants. Our results demonstrate that the SgCINS gene can play an important role in plants against cotton worm insect attack, and pave the way for using terpenoids genes for improving resistance to insect attack in higher plants.

Terpene synthases in cucumber (Cucumis sativus) and their contribution to herbivore-induced volatile terpenoid emission

Terpenoids play important roles in flavour, pollinator attraction and defence of plants. In cucumber (Cucumis sativus) they are important components of the herbivore-induced plant volatile blend that attracts natural enemies of herbivores. We annotated the cucumber TERPENE SYNTHASE gene (CsTPS) family and characterized their involvement in the response towards herbivores with different feeding guilds using a combined molecular and biochemical approach. Transcripts of multiple CsTPS genes were upregulated in leaves upon herbivory and the products generated by the expressed proteins match the terpenoids recorded in the volatile blend released by herbivore-damaged leaves. Spatial and temporal analysis of the promoter activity of CsTPS genes showed that cell content-feeding spider mites (Tetranychus urticae) and thrips (Frankliniella occidentalis) induced promoter activity of CsTPS9 and CsTPS19 within hours after initiation of infestation, while phloem-feeding aphids (Myzus persicae) induced CsTPS2 promoter activity. Our findings offer detailed insights into the involvement of the TPS gene family in the dynamics and fine-tuning of the emission of herbivore-induced plant volatiles in cucumber, and open a new avenue to understand molecular mechanisms that affect plant-herbivore interactions.

Application of High-Throughput Sequencing on the Chinese Herbal Medicine for the Data-Mining of the Bioactive Compounds

The Chinese Herbal Medicine (CHM) has been used worldwide in clinic to treat the vast majority of human diseases, and the healing effect is remarkable. However, the functional components and the corresponding pharmacological mechanism of the herbs are unclear. As one of the main means, the high-throughput sequencing (HTS) technologies have been employed to discover and parse the active ingredients of CHM. Moreover, a tremendous amount of effort is made to uncover the pharmacodynamic genes associated with the synthesis of active substances. Here, based on the genome-assembly and the downstream bioinformatics analysis, we present a comprehensive summary of the application of HTS on CHM for the synthesis pathways of active ingredients from two aspects: active ingredient properties and disease classification, which are important for pharmacological, herb molecular breeding, and synthetic biology studies.

A Bioinformatics Tool for Efficient Retrieval of High-Confidence Terpene Synthases (TPS) and Application to the Identification of TPS in Coffea and Quillaja

Terpenoids are a class of compounds that are found in all living organisms. In plants, some terpenoids are part of primary metabolism, but most terpenes found in plants are classified as specialized metabolites, encoded by terpene synthases (TPS). It is not obvious how to assign the putative product of a given TPS using bioinformatics tools. Phylogenetic analyses easily assign TPS into families; however members of the same TPS family can synthetize more than one terpenoid-and, in many biotechnological applications, researchers are more interested in the product of a given TPS rather than its phylogenetic profile. Automated protein annotation can be used to classify TPS based on their products, despite the family they belong to. Here, we implement an automated bioinformatics method, search_TPS, to identify TPS proteins that synthesize mono, sesqui and diterpenes in Angiosperms. We verified the applicability of the method by classifying wet lab validated TPS and applying it to find TPS proteins in Coffea arabica, C. canephora, C. eugenioides, and Quillaja saponaria. Search_TPS is a computational tool based on PERL scripts that carries out a series of HMMER searches against a curated database of TPS profile hidden Markov models. The tool is freely available at https://github.com/liliane-sntn/TPS .

The biochemical and molecular investigation of flower color and scent sheds lights on further genetic modification of ornamental traits in Clivia miniata

Clivia miniata is renowned for its evergreen and strap-like leaves, whereas its floral color and scent are lacking diversity. Here, anthocyanin, volatile terpene, and carotenoid metabolisms were integrally investigated in C. miniata flowers. The results showed that pelargonidins and lutein might cooperate to confer orange or yellow color to C. miniata flowers, but only a trace amount of (+)-limonene was detected. The expression levels of CmF3'H and CmDFR appeared to be responsible for the ratio of cyanidin and pelargonidin derivatives in C. miniata, and the low expression of CmF3'H was responsible for the lack of cyanidins in flowers. Moreover, the CmF3'H promoter could not be activated by CmMYBAs, suggesting that it was controlled by novel regulators. Only two CmTPSs were functional, with CmTPS2 responsible for (+)-limonene synthesis, contributing to the monotonous flower volatile terpenes of C. miniata. CmCCD1a and CmCCD1b were able to cleave carotenoids at the 5,6 (5',6'), and 9,10 (9',10') positions to generate volatile apocarotenoids, whereas the substrates found in low-quantities or specific subcellular localizations of CmCCD1s might constrain volatile apocarotenoid release. Consequently, activating F3'H and introducing novel F3'5'H or versatile TPS may be effective ways to modify the floral color and scent, respectively. Alternatively, modifying the carotenoid flux or CCD1 localization might affect floral color and scent simultaneously. Taking these results together, the present study provides a preliminary deciphering of the genetic constraints underlying flower color and scent development, and proposes possible schemes for further genetic modification of ornamental traits in C. miniata and other plants.

Paeoniflorin in Paeoniaceae: Distribution, influencing factors, and biosynthesis

Paeoniflorin, a monoterpene glucoside, is increasingly used in the clinical treatment of many diseases because it has a variety of pharmacological activities. Besides, paeoniflorin has been considered the characteristic chemical constituent of Paeoniaceae plants since it was first reported in 1963. Therefore, how to better develop and utilize paeoniflorin in Paeoniaceae has always been a research hotspot. We reviewed the current knowledge on paeoniflorin in Paeoniaceae, with particular emphasis on its distribution and influencing factors. Moreover, the limited understanding of the biosynthesis pathway has restricted the production of paeoniflorin by synthetic biology. This review provides insights into the post-modification pathway of paeoniflorin biosynthesis and proposes directions for further analysis in the future.

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.

Short-term changes related to autotetraploidy in essential oil composition of Eucalyptus benthamii Maiden & Cambage and its applications in different bioassays

Some forest trees have been polyploidized to improve their traits and to supply new germplasms for breeding programs. As trees have a long juvenile stage, the early characterization of the chromosome set doubling effects is crucial for previous selection. Thus, we aimed to characterize the chemical variability of essential oils from diploid and autotetraploid germplasms (autotetraploid A and B) of Eucalyptus benthamii, as well as to evaluate their larvicidal and allelopathic effects. Autotetraploid A showed a higher essential oil yield than diploid and autotetraploid B, which did not differ quantitatively. Aromadendrene, viridiflorol and α-pinene were the major compounds in the diploid essential oil. In contrast, compounds were present in autotetraploids, such as 1,8-cineole, limonene, α-terpineol, and α-terpinyl-acetate. Essential oils from the diploid at 50-200 ppm were twice as larvicidal than those from autotetraploids against Aedes aegypti larvae. Considering the phytotoxicity bioassays using Lactuca sativa, essential oils from both ploidy levels affected root growth. Moreover, the essential oils inhibited shoot growth at all concentrations tested (187.5; 375; 750; 1500; and 3000 ppm). Autotetraploid A and B had the same effect on shoot growth as glyphosate. The essential oils had no cytogenotoxic effect on root meristematic cells of L. sativa, whereas phytotoxic potential was identified mainly in shoot growth. This work demonstrated a dramatic change in secondary metabolism (terpene composition) related to an increase in the ploidy level in Eucalyptus germplasms. In addition, we report the novelty of the chemical composition of essential oils among germplasms and their potential use as larvicidal and post-emergence weed control agents.

High-quality evergreen azalea genome reveals tandem duplication-facilitated low-altitude adaptability and floral scent evolution

Azalea belongs to Rhododendron, which is one of the largest genera of flowering plants and is well known for the diversity and beauty in its more than 1000 woody species. Rhododendron contains two distinct groups: the most high-altitude and a few low-altitude species; however, the former group is difficult to be domesticated for urban landscaping, and their evolution and adaptation are little known. Rhododendron ovatum has broad adaptation in low-altitude regions but possesses evergreen characteristics like high-altitude species, and it has floral fragrance that is deficient in most cultivars. Here we report the chromosome-level genome assembly of R. ovatum, which has a total length of 549 Mb with scaffold N50 of 41 Mb and contains 41 264 predicted genes. Genomic micro-evolutionary analysis of R. ovatum in comparison with two high-altitude Rhododendron species indicated that the expansion genes in R. ovatum were significantly enriched in defence responses, which may account for its adaptability in low altitudes. The R. ovatum genome contains much more terpene synthase genes (TPSs) compared with the species that lost floral fragrance. The subfamily b members of TPS are involved in the synthesis of sesquiterpenes as well as monoterpenes and play a major role in flora scent biosynthesis and defence responses. Tandem duplication is the primary force driving expansion of defence-responsive genes for extensive adaptability to the low-altitude environments. The R. ovatum genome provides insights into low-altitude adaptation and gain or loss of floral fragrance for Rhododendron species, which are valuable for alpine plant domestication and floral scent breeding.

The genome of Cymbidium sinense revealed the evolution of orchid traits

The Orchidaceae is of economic and ecological importance and constitutes ˜10% of all seed plant species. Here, we report a genome physical map for Cymbidium sinense, a well-known species belonging to genus Cymbidium that has thousands of natural variation varieties of flower organs, flower and leaf colours and also referred as the King of Fragrance, which make it arose into a unique cultural symbol in China. The high-quality chromosome-scale genome assembly was 3.52 Gb in size, 29 638 protein-coding genes were predicted, and evidence for whole-genome duplication shared with other orchids was provided. Marked amplification of cytochrome- and photosystem-related genes was observed, which was consistent with the shade tolerance and dark green leaves of C. sinense. Extensive duplication of MADS-box genes, and the resulting subfunctional and expressional differentiation, was associated with regulation of species-specific flower traits, including wild-type and mutant-type floral patterning, seasonal flowering and ecological adaption. CsSEP4 was originally found to positively regulate gynostemium development. The CsSVP genes and their interaction proteins CsAP1 and CsSOC1 were significantly expanded and involved in the regulation of low-temperature-dependent flowering. Important genetic clues to the colourful leaf traits, purple-black flowers and volatile trait in C. sinense were also found. The results provide new insights into the molecular mechanisms of important phenotypic traits of Cymbidium and its evolution and serve as a powerful platform for future evolutionary studies and molecular breeding of orchids.

Genome-Wide Identification and Expression Analysis of Terpene Synthase Genes in Cymbidium faberi

Terpene synthases (TPSs) are essential for forming terpenes, which play numerous functional roles in attracting pollinators, defending plants, and moderating the interaction between plants. TPSs have been reported in some orchids, but genome-wide identification of terpenes in Cymbidium faberi is still lacking. In this study, 32 putative TPS genes were classified in C. faberi and divided into three subfamilies (TPS-a, TPS-b, and TPS-e/f). Motif and gene structure analysis revealed that most CfTPS genes had the conserved aspartate-rich DDxxD motif. TPS genes in the TPS-a and TPS-b subfamilies had variations in the RRX₈W motif. Most cis-elements of CfTPS genes were found in the phytohormone responsiveness category, and MYC contained most of the numbers associated with MeJA responsiveness. The Ka/Ks ratios of 12/13 CfTPS gene pairs were less than one, indicated that most CfTPS genes have undergone negative selection. The tissue-specific expression patterns showed that 28 genes were expressed in at least one tissue in C. faberi, and TPS genes were most highly expressed in flowers, followed by leaves and pseudobulbs. In addition, four CfTPS genes were selected for the real-time reverse transcription quantitative PCR (RT-qPCR) experiment. The results revealed that CfTPS12, CfTPS18, CfTPS23, and CfTPS28 were mainly expressed in the full flowering stage. CfTPS18 could convert GPP to β-myrcene, geraniol, and α-pinene in vitro. These findings of CfTPS genes of C. faberi may provide valuable information for further studies on TPSs in orchids.

Terpene Synthase Gene OtLIS Confers Wheat Resistance to Sitobion avenae by Regulating Linalool Emission

Sitobion avenae (Fabricius) is a major insect pest of wheat worldwide that reduces crop yield and quality annually. Few germplasm resources with resistant genes to aphids have been identified and characterized. Here, octoploid Trititrigia, a species used in wheat distant hybridization breeding, was found to be repellent to S. avenae after 2 year field investigations and associated with physiological and behavioral assays. Linalool monoterpene was identified to accumulate dominantly in plants in response to S. avenae infestation. We cloned the resistance gene OtLIS by assembling the transcriptome of aphid-infested or healthy octoploid Trititrigia. Functional characterization analysis indicated that OtLIS encoded a terpene synthase and conferred resistance to S. avenae by linalool emission before and after aphid feeding. Our study suggests that the octoploid Trititrigia with the aphid-resistant gene OtLIS may have potential as a target resource for further breeding aphid-resistant wheat cultivars.

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.

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.

RNASeq analysis of drought-stressed guayule reveals the role of gene transcription for modulating rubber, resin, and carbohydrate synthesis

The drought-adapted shrub guayule (Parthenium argentatum) produces rubber, a natural product of major commercial importance, and two co-products with potential industrial use: terpene resin and the carbohydrate fructan. The rubber content of guayule plants subjected to water stress is higher compared to that of well-irrigated plants, a fact consistently reported in guayule field evaluations. To better understand how drought influences rubber biosynthesis at the molecular level, a comprehensive transcriptome database was built from drought-stressed guayule stem tissues using de novo RNA-seq and genome-guided assembly, followed by annotation and expression analysis. Despite having higher rubber content, most rubber biosynthesis related genes were down-regulated in drought-stressed guayule, compared to well-irrigated plants, suggesting post-transcriptional effects may regulate drought-induced rubber accumulation. On the other hand, terpene resin biosynthesis genes were unevenly affected by water stress, implying unique environmental influences over transcriptional control of different terpene compounds or classes. Finally, drought induced expression of fructan catabolism genes in guayule and significantly suppressed these fructan biosynthesis genes. It appears then, that in guayule cultivation, irrigation levels might be calibrated in such a regime to enable tunable accumulation of rubber, resin and fructan.

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

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

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.

Predictive Engineering of Class I Terpene Synthases Using Experimental and Computational Approaches

Terpenoids are a highly diverse group of natural products with considerable industrial interest. Increasingly, engineered microbes are used for the production of terpenoids to replace natural extracts and chemical synthesis. Terpene synthases (TSs) show a high level of functional plasticity and are responsible for the vast structural diversity observed in natural terpenoids. Their relatively inert active sites guide intrinsically reactive linear carbocation intermediates along one of many cyclisation paths via exertion of subtle steric and electrostatic control. Due to the absence of a strong protein interaction with these intermediates, there is a remarkable lack of sequence‐function relationship within the TS family, making product‐outcome predictions from sequences alone challenging. This, in combination with the fact that many TSs produce multiple products from a single substrate hampers the design and use of TSs in the biomanufacturing of terpenoids. This review highlights recent advances in genome mining, computational modelling, high‐throughput screening, and machine‐learning that will allow more predictive engineering of these fascinating enzymes in the near future.

Metabolic Perturbation and Synthetic Biology Strategies for Plant Terpenoid Production-An Updated Overview

Terpenoids represent one of the high-value groups of specialized metabolites with vast structural diversity. They exhibit versatile human benefits and have been successfully exploited in several sectors of day-to-day life applications, including cosmetics, foods, and pharmaceuticals. Historically, the potential use of terpenoids is challenging, and highly hampered by their bioavailability in their natural sources. Significant progress has been made in recent years to overcome such challenges by advancing the heterologous production platforms of hosts and metabolic engineering technologies. Herein, we summarize the latest developments associated with analytical platforms, metabolic engineering, and synthetic biology, with a focus on two terpenoid classes: monoterpenoids and sesquiterpenoids. Accumulated data showed that subcellular localization of both the precursor pool and the introduced enzymes were the crucial factors for increasing the production of targeted terpenoids in plants. We believe this timely review provides a glimpse of current state-of-the-art techniques/methodologies related to terpenoid engineering that would facilitate further improvements in terpenoids research.

Sesquiterpene biosynthesis in a leafy liverwort Radula lindenbergiana Gottsche ex C. Hartm

Liverworts (Marchantiophyta) are among the earliest diverging lineages of extant land plants. Among their unique features, most liverworts contain membrane-bound oil bodies, organelles that accumulate diverse secondary metabolites, especially terpenoids. In contrast to the rich information on liverwort terpenoid chemistry, little is known about their biosynthesis. Recently, terpenoid biosynthesis was studied in a model thalloid species Marchantiapolymorpha, in which sesquiterpenes and monoterpenes are biosynthesized by a new type of terpene synthases termed microbial terpene synthase-like (MTPSL) proteins. Here we study terpenoid biosynthesis in a leafy liverwort Radula lindenbergiana. Vegetative plants of R.lindenbergiana were found to contain a mixture of sesquiterpenes, with (E,E)-α-farnesene/β-curcumene and (Z)-β-bisabolene being the most abundant constituents. From the analysis of the R. lindenbergiana transcriptome, five full-length MTPSL genes were identified. They were designated RlMTPSL1-5, respectively. Recombinant RlMTPSL proteins were produced in Escherichia coli and tested for sesquiterpene synthase activities using farnesyl diphosphate (FPP) as substrate. All except RlMTPSL5 were demonstrated to catalyze the formation of different sesquiterpenes. RlMTPSL1 produced multiple sesquiterpenes with eremophilene and an unidentified sesquiterpene as major products. The major products of RlMTPSL2 and RlMTPSL3 were β-elemene and an unidentified sesquiterpene, respectively. RlMTPSL4 was also a multi-product sesquiterpene synthase with an unidentified sesquiterpene being the major product. Homology-based structural modeling was performed to understand the structural basis underlying different product profiles of the RlMTPSLs proteins. Most of the sesquiterpene products of the four active RlMTPSLs were also detected in R. lindenbergiana plants. Expression levels of the four RlMTPSL genes encoding active enzymes in vegetative plants were compared. In phylogenetic analysis, RlMTPSL genes were found to cluster together, indicating lineage-specific expansion of MTPSL genes in lineages leading to R.lindenbergiana and M. polymorpha. This study strengthens evidence for the contribution of MTPSL genes to terpenoid biosynthesis in liverworts.

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

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

Analgesic Potential of Terpenes Derived from Cannabis sativa

Pain prevalence among adults in the United States has increased 25% over the past two decades, resulting in high health-care costs and impacts to patient quality of life. In the last 30 years, our understanding of pain circuits and (intra)cellular mechanisms has grown exponentially, but this understanding has not yet resulted in improved therapies. Options for pain management are limited. Many analgesics have poor efficacy and are accompanied by severe side effects such as addiction, resulting in a devastating opioid abuse and overdose epidemic. These problems have encouraged scientists to identify novel molecular targets and develop alternative pain therapeutics. Increasing preclinical and clinical evidence suggests that cannabis has several beneficial pharmacological activities, including pain relief. Cannabis sativa contains more than 500 chemical compounds, with two principle phytocannabinoids, Δ9-tetrahydrocannabinol (Δ9-THC) and cannabidiol (CBD). Beyond phytocannabinoids, more than 150 terpenes have been identified in different cannabis chemovars. Although the predominant cannabinoids, Δ9-THC and CBD, are thought to be the primary medicinal compounds, terpenes including the monoterpenes β-myrcene, α-pinene, limonene, and linalool, as well as the sesquiterpenes β-caryophyllene and α-humulene may contribute to many pharmacological properties of cannabis, including anti-inflammatory and antinociceptive effects. The aim of this review is to summarize our current knowledge about terpene compounds in cannabis and to analyze the available scientific evidence for a role of cannabis-derived terpenes in modern pain management. SIGNIFICANCE STATEMENT: Decades of research have improved our knowledge of cannabis polypharmacy and contributing phytochemicals, including terpenes. Reform of the legal status for cannabis possession and increased availability (medicinal and recreational) have resulted in cannabis use to combat the increasing prevalence of pain and may help to address the opioid crisis. Better understanding of the pharmacological effects of cannabis and its active components, including terpenes, may assist in identifying new therapeutic approaches and optimizing the use of cannabis and/or terpenes as analgesic agents.

Elucidation of the Regular Emission Mechanism of Volatile β-Ocimene with Anti-insect Function from Tea Plants (Camellia sinensis) Exposed to Herbivore Attack

Herbivore-induced plant volatiles (HIPVs) play an important role in insect resistance. As a common HIPV in tea plants (Camellia sinensis), β-ocimene has shown anti-insect function in other plants. However, whether β-ocimene in tea plants also provides insect resistance, and its mechanism of synthesis and emission are unknown. In this study, β-ocimene was confirmed to interfere with tea geometrid growth via signaling. Light was identified as the key factor controlling regular emission of β-ocimene induced by the wounding from tea geometrids. β-Ocimene synthase (CsBOS1) was located in plastids and catalyzed β-ocimene formation in overexpressed tobacco. CsBOS1 expression in tea leaves attacked by tea geometrids showed a day-low and night-high variation pattern, while CsABCG expression involved in volatile emission showed the opposite pattern. These two genes might regulate the regular β-ocimene emission from tea plants induced by tea geometrid attack. This study advances the understanding on HIPV emission and signaling in tea plants.

Whole-genome assembly of Ganoderma leucocontextum (Ganodermataceae, Fungi) discovered from the Tibetan Plateau of China

Ganoderma leucocontextum, a newly discovered species of Ganodermataceae in China, has diverse pharmacological activities. G. leucocontextum was widely cultivated in southwest China, but the systematic genetic study has been impeded by the lack of a reference genome. Herein, we present the first whole-genome assembly of G. leucocontextum based on the Illumina and Nanopore platform from high-quality DNA extracted from a monokaryon strain (DH-8). The generated genome was 50.05 Mb in size with a N50 scaffold size of 3.06 Mb, 78,206 coding sequences and 13,390 putative genes. Genome completeness was assessed using the Benchmarking Universal Single-Copy Orthologs (BUSCO) tool, which identified 96.55% of the 280 Fungi BUSCO genes. Furthermore, differences in functional genes of secondary metabolites (terpenoids) were analyzed between G. leucocontextum and G. lucidum. G. leucocontextum has more genes related to terpenoids synthesis compared to G. lucidum, which may be one of the reasons why they exhibit different biological activities. This is the first genome assembly and annotation for G. leucocontextum, which would enrich the toolbox for biological and genetic studies in G. leucocontextum.