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

A high-throughput colorimetric screening assay for terpene synthase activity based on substrate consumption

Terpene synthases catalyze the formation of a variety of terpene chemical structures. Systematic mutagenesis studies have been effective in providing insights into the characteristic and complex mechanisms of C-C bond formations and in exploring the enzymatic potential for inventing new chemical structures. In addition, there is growing demand to increase terpene synthase activity in heterologous hosts, given the maturation of metabolic engineering and host breeding for terpenoid synthesis. We have developed a simple screening method for the cellular activities of terpene synthases by scoring their substrate consumption based on the color loss of the cell harboring carotenoid pathways. We demonstrate that this method can be used to detect activities of various terpene synthase or prenyltransferase genes in a high-throughput manner, irrespective of the product type, enabling the mutation analysis and directed evolution of terpene synthases. We also report the possibility for substrate-specific screening system of terpene synthases by taking advantage of the substrate-size specificity of C₃₀ and C₄₀ carotenoid pathways.

The sex-linked region in Populus tremuloides Turesson 141 corresponds to a pericentromeric region of about two million base pairs on P. trichocarpa chromosome 19

In the dioecious genus Populus, sex determination has been located to chromosome 19. However, despite a high degree of genome collinearity, various Populus species seem to differ with regard to the location of the sex-determining region on the respective chromosome and the apparent heterogametic sex. In this study, the boundaries of the recombination-suppressed, sex-linked region of the male P. tremuloides clone Turesson 141 were localised by genetic mapping using new SNP and InDel markers. The respective region seems to be located in a pericentromeric position. The corresponding P. trichocarpa genome region spans about two million bp and comprises 65 gene loci, which were bioinformatically evaluated for their potential as candidate genes for sex determination. Three putative transcription factor genes and four genes that are potentially involved in flower development processes, e.g. meristem transition from the vegetative to the reproductive phase, were identified. Populus tremuloides sequence data of the sex-linked region is required for a final search for candidate genes.

The evolutionary paths towards complexity: a metabolic perspective

As sessile organisms, land plants have exploited their metabolic systems to produce a panoply of structurally and functionally diverse natural chemicals and polymers to adapt to challenging ecosystems. Many of these core and specialized metabolites confer chemical shields against a multitude of abiotic stresses, while others play important roles in plants' interactions with their biotic environments. Plant specialized metabolites can be viewed as complex traits in the sense that the biosynthesis of these molecules typically requires multistep metabolic pathways comprising numerous specific enzymes belonging to diverse protein fold families. Resolving the evolutionary trajectories underlying the emergence of these specialized metabolic pathways will impact a fundamental question in biology – how do complex traits evolve in a Darwinian fashion? Here, I discuss several general patterns observed in rapidly evolving specialized metabolic systems in plants, and surmise mechanistic features at enzyme, pathway and organismal levels that rationalize the remarkable malleability of these systems through stepwise evolution. Future studies, focused on fine sampling of metabolic enzymes and pathways in phylogenetically related plant species, or employing directed evolution strategies in synthetic systems, will significantly broaden our perspective on how biological complexity arises at the metabolic level.

Overexpression of an isoprenyl diphosphate synthase in spruce leads to unexpected terpene diversion products that function in plant defense

Spruce (Picea spp.) and other conifers employ terpenoid-based oleoresin as part of their defense against herbivores and pathogens. The short-chain isoprenyl diphosphate synthases (IDS) are situated at critical branch points in terpene biosynthesis, producing the precursors of the different terpenoid classes. To determine the role of IDS and to create altered terpene phenotypes for assessing the defensive role of terpenoids, we overexpressed a bifunctional spruce IDS, a geranyl diphosphate and geranylgeranyl diphosphate synthase in white spruce (Picea glauca) saplings. While transcript level (350-fold), enzyme activity level (7-fold), and in planta geranyl diphosphate and geranylgeranyl diphosphate levels (4- to 8-fold) were significantly increased in the needles of transgenic plants, there was no increase in the major monoterpenes and diterpene acids of the resin and no change in primary isoprenoids, such as sterols, chlorophylls, and carotenoids. Instead, large amounts of geranylgeranyl fatty acid esters, known from various gymnosperm and angiosperm plant species, accumulated in needles and were shown to act defensively in reducing the performance of larvae of the nun moth (Lymantria monacha), a conifer pest in Eurasia. These results show the impact of overexpression of an IDS and the defensive role of an unexpected accumulation product of terpenoid biosynthesis with the potential for a broader function in plant protection.

Jasmonate induction of the monoterpene linalool confers resistance to rice bacterial blight and its biosynthesis is regulated by JAZ protein in rice

Jasmonic acid (JA) is involved in the regulation of host immunity in plants. Recently, we demonstrated that JA signalling has an important role in resistance to rice bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo) in rice. Here, we report that many volatile compounds accumulate in response to exogenous application of JA, including the monoterpene linalool. Expression of linalool synthase was up-regulated by JA. Vapour treatment with linalool induced resistance to Xoo, and transgenic rice plants overexpressing linalool synthase were more resistance to Xoo, presumably due to the up-regulation of defence-related genes in the absence of any treatment. JA-induced accumulation of linalool was regulated by OsJAZ8, a rice jasmonate ZIM-domain protein involving the JA signalling pathway at the transcriptional level, suggesting that linalool plays an important role in JA-induced resistance to Xoo in rice.

Valencene synthase from the heartwood of Nootka cypress (Callitropsis nootkatensis) for biotechnological production of valencene

Nootkatone is one of the major terpenes in the heartwood of the Nootka cypress Callitropsis nootkatensis. It is an oxidized sesquiterpene, which has been postulated to be derived from valencene. Both valencene and nootkatone are used for flavouring citrus beverages and are considered among the most valuable terpenes used at commercial scale. Functional evaluation of putative terpene synthase genes sourced by large-scale EST sequencing from Nootka cypress wood revealed a valencene synthase gene (CnVS). CnVS expression in different tissues from the tree correlates well with nootkatone content, suggesting that CnVS represents the first dedicated gene in the nootkatone biosynthetic pathway in C. nootkatensis The gene belongs to the gymnosperm-specific TPS-d subfamily of terpenes synthases and its protein sequence has low similarity to known citrus valencene synthases. In vitro, CnVS displays high robustness under different pH and temperature regimes, potentially beneficial properties for application in different host and physiological conditions. Biotechnological production of sesquiterpenes has been shown to be feasible, but productivity of microbial strains expressing valencene synthase from Citrus is low, indicating that optimization of valencene synthase activity is needed. Indeed, expression of CnVS in Saccharomyces cerevisiae indicated potential for higher yields. In an optimized Rhodobacter sphaeroides strain, expression of CnVS increased valencene yields 14-fold to 352 mg/L, bringing production to levels with industrial potential.

Explaining intraspecific diversity in plant secondary metabolites in an ecological context

Plant secondary metabolites (PSMs) are ubiquitous in plants and play many ecological roles. Each compound can vary in presence and/or quantity, and the composition of the mixture of chemicals can vary, such that chemodiversity can be partitioned within and among individuals. Plant ontogeny and environmental and genetic variation are recognized as sources of chemical variation, but recent advances in understanding the molecular basis of variation may allow the future deployment of isogenic mutants to test the specific adaptive function of variation in PSMs. An important consequence of high intraspecific variation is the capacity to evolve rapidly. It is becoming increasingly clear that trait variance linked to both macro- and micro-environmental variation can also evolve and may respond more strongly to selection than mean trait values. This research, which is in its infancy in plants, highlights what could be a missing piece of the picture of PSM evolution. PSM polymorphisms are probably maintained by multiple selective forces acting across many spatial and temporal scales, but convincing examples that recognize the diversity of plant population structures are rare. We describe how diversity can be inherently beneficial for plants and suggest fruitful avenues for future research to untangle the causes and consequences of intraspecific variation.

Multi-platform metabolomic analyses of ergosterol-induced dynamic changes in Nicotiana tabacum cells

Metabolomics is providing new dimensions into understanding the intracellular adaptive responses in plants to external stimuli. In this study, a multi-technology-metabolomic approach was used to investigate the effect of the fungal sterol, ergosterol, on the metabolome of cultured tobacco cells. Cell suspensions were treated with different concentrations (0-1000 nM) of ergosterol and incubated for different time periods (0-24 h). Intracellular metabolites were extracted with two methods: a selective dispersive liquid-liquid micro-extraction and a general methanol extraction. Chromatographic techniques (GC-FID, GC-MS, GC × GC-TOF-MS, UHPLC-MS) and (1)H NMR spectroscopy were used for quantitative and qualitative analyses. Multivariate data analyses (PCA and OPLS-DA models) were used to extract interpretable information from the multidimensional data generated from the analytical techniques. The results showed that ergosterol triggered differential changes in the metabolome of the cells, leading to variation in the biosynthesis of secondary metabolites. PCA scores plots revealed dose- and time-dependent metabolic variations, with optimal treatment conditions being found to be 300 nM ergosterol and an 18 h incubation period. The observed ergosterol-induced metabolic changes were correlated with changes in defence-related metabolites. The 'defensome' involved increases in terpenoid metabolites with five antimicrobial compounds (the bicyclic sesquiterpenoid phytoalexins: phytuberin, solavetivone, capsidiol, lubimin and rishitin) and other metabolites (abscisic acid and phytosterols) putatively identified. In addition, various phenylpropanoid precursors, cinnamic acid derivatives and - conjugates, coumarins and lignin monomers were annotated. These annotated metabolites revealed a dynamic reprogramming of metabolic networks that are functionally correlated, with a high complexity in their regulation.

The Amborella genome and the evolution of flowering plants

Amborella trichopoda is strongly supported as the single living species of the sister lineage to all other extant flowering plants, providing a unique reference for inferring the genome content and structure of the most recent common ancestor (MRCA) of living angiosperms. Sequencing the Amborella genome, we identified an ancient genome duplication predating angiosperm diversification, without evidence of subsequent, lineage-specific genome duplications. Comparisons between Amborella and other angiosperms facilitated reconstruction of the ancestral angiosperm gene content and gene order in the MRCA of core eudicots. We identify new gene families, gene duplications, and floral protein-protein interactions that first appeared in the ancestral angiosperm. Transposable elements in Amborella are ancient and highly divergent, with no recent transposon radiations. Population genomic analysis across Amborella's native range in New Caledonia reveals a recent genetic bottleneck and geographic structure with conservation implications.

To gibberellins and beyond! Surveying the evolution of (di)terpenoid metabolism

The diterpenoids are classically defined by their composition--four isoprenyl units (20 carbons)--and are generally derived from [E,E,E]-geranylgeranyl diphosphate (GGPP). Such metabolism seems to be ancient and has been extensively diversified, with ∼12,000 diterpenoid natural products known. Particularly notable are the gibberellin phytohormones, whose requisite biosynthesis has provided a genetic reservoir that gave rise to not only a large superfamily of ∼7,000 diterpenoids but also, to some degree, all plant terpenoid natural products. This review focuses on the diterpenoids, particularly the defining biosynthetic characteristics of the major superfamilies defined by the cyclization and/or rearrangement of GGPP catalyzed by diterpene synthases/cyclases, although it also includes some discussion of the important subsequent elaboration in the few cases where sufficient molecular genetic information is available. It additionally addresses the array of biological activity providing the selective pressures that drive the observed gene family expansion and diversification, along with biosynthetic gene clustering.

Biosynthesis of terpenoid natural products in fungi

Tens of thousands of terpenoid natural products have been isolated from plants and microbial sources. Higher fungi (Ascomycota and Basidiomycota) are known to produce an array of well-known terpenoid natural products, including mycotoxins, antibiotics, antitumor compounds, and phytohormones. Except for a few well-studied fungal biosynthetic pathways, the majority of genes and biosynthetic pathways responsible for the biosynthesis of a small number of these secondary metabolites have only been discovered and characterized in the past 5-10 years. This chapter provides a comprehensive overview of the current knowledge on fungal terpenoid biosynthesis from biochemical, genetic, and genomic viewpoints. Enzymes involved in synthesizing, transferring, and cyclizing the prenyl chains that form the hydrocarbon scaffolds of fungal terpenoid natural products are systematically discussed. Genomic information and functional evidence suggest differences between the terpenome of the two major fungal phyla--the Ascomycota and Basidiomycota--which will be illustrated for each group of terpenoid natural products.

High-throughput testing of terpenoid biosynthesis candidate genes using transient expression in Nicotiana benthamiana

To respond to the rapidly growing number of genes putatively involved in terpenoid metabolism, a robust high-throughput platform for functional testing is needed. An in planta expression system offers several advantages such as the capacity to produce correctly folded and active enzymes localized to the native compartments, unlike microbial or prokaryotic expression systems. Two inherent drawbacks of plant-based expression systems, time-consuming generation of transgenic plant lines and challenging gene-stacking, can be circumvented by transient expression in Nicotiana benthamiana. In this chapter we describe an expression platform for rapid testing of candidate terpenoid biosynthetic genes based on Agrobacterium mediated gene expression in N. benthamiana leaves. Simultaneous expression of multiple genes is facilitated by co-infiltration of leaves with several engineered Agrobacterium strains, possibly making this the fastest and most convenient system for the assembly of plant terpenoid biosynthetic routes. Tools for cloning of expression plasmids, N. benthamiana culturing, Agrobacterium preparation, leaf infiltration, metabolite extraction, and automated GC-MS data mining are provided. With all steps optimized for high throughput, this in planta expression platform is particularly suited for testing large panels of candidate genes in all possible permutations.

Heterologous stable expression of terpenoid biosynthetic genes using the moss Physcomitrella patens

Heterologous and stable expression of genes encoding terpenoid biosynthetic enzymes in planta is an important tool for functional characterization and is an attractive alternative to expression in microbial hosts for biotechnological production. Despite improvements to the procedure, such as streamlining of large scale Agrobacterium infiltration and upregulation of the upstream pathways, transient in planta heterologous expression quickly reaches limitations when used for production of terpenoids. Stable integration of transgenes into the nuclear genome of the moss Physcomitrella patens has already been widely recognized as a viable alternative for industrial-scale production of biopharmaceuticals. For expression of terpenoid biosynthetic genes, and reconstruction of heterologous pathways, Physcomitrella has unique attributes that makes it a very promising biotechnological host. These features include a high native tolerance to terpenoids, a simple endogenous terpenoid profile, convenient genome editing using homologous recombination, and cultivation techniques that allow up-scaling from single cells in microtiter plates to industrial photo-bioreactors. Beyond its use for functional characterization of terpenoid biosynthetic genes, engineered Physcomitrella can be a green biotechnological platform for production of terpenoids. Here, we describe two complementary and simple procedures for stable nuclear transformation of Physcomitrella with terpenoid biosynthetic genes, selection and cultivation of transgenic lines, and metabolite analysis of terpenoids produced in transgenic moss lines. We also provide tools for metabolic engineering through genome editing using homologous recombination.

Determination of residues responsible for substrate and product specificity of Solanum habrochaites short-chain cis-prenyltransferases

Isoprenoids are diverse compounds that have their biosynthetic origin in the initial condensation of isopentenyl diphosphate and dimethylallyl diphosphate to form C10 prenyl diphosphates that can be elongated by the addition of subsequent isopentenyl diphosphate units. These reactions are catalyzed by either cis-prenyltransferases (CPTs) or trans-prenyltransferases. The synthesis of volatile terpenes in plants typically proceeds through either geranyl diphosphate (C10) or trans-farnesyl diphosphate (C15), to yield monoterpenes and sesquiterpenes, respectively. However, terpene biosynthesis in glandular trichomes of tomato (Solanum lycopersicum) and related wild relatives also occurs via the cis-substrates neryl diphosphate (NPP) and 2Z,6Z-farnesyl diphosphate (Z,Z-FPP). NPP and Z,Z-FPP are synthesized by neryl diphosphate synthase1 (NDPS1) and Z,Z-farnesyl diphosphate synthase (zFPS), which are encoded by the orthologous CPT1 locus in tomato and Solanum habrochaites, respectively. In this study, comparative sequence analysis of NDPS1 and zFPS enzymes from S. habrochaites accessions that synthesize either monoterpenes or sesquiterpenes was performed to identify amino acid residues that correlate with the ability to synthesize NPP or Z,Z-FPP. Subsequent structural modeling, coupled with site-directed mutagenesis, highlighted the importance of four amino acids located within conserved domain II of CPT enzymes that form part of the second α-helix, for determining substrate and product specificity of these enzymes. In particular, the relative positioning of aromatic amino acid residues at positions 100 and 107 determines the ability of these enzymes to synthesize NPP or Z,Z-FPP. This study provides insight into the biochemical evolution of terpene biosynthesis in the glandular trichomes of Solanum species.

Mushroom hunting by using bioinformatics: application of a predictive framework facilitates the selective identification of sesquiterpene synthases in basidiomycota

The Basidiomycota fungi represent a diverse source of natural products, particularly the sesquiterpenoids. Recently, genome sequencing, genome mining, and the subsequent discovery of a suite of sesquiterpene synthases in Omphalotus olearius was described. A predictive framework was developed to facilitate the discovery of sesquiterpene synthases in Basidiomycota. Phylogenetic analyses indicated a conservation of both sequence and initial cyclization mechanisms used. Here, the first robust application of this predictive framework is reported. It was used to selectively identify sesquiterpene synthases that follow 1,6-, 1,10-, and 1,11-cyclization mechanisms in the crust fungus Stereum hirsutum. The successful identification and characterization of a 1,6- and a 1,10-cyclizing sesquiterpene synthase, as well as three 1,11-cyclizing Δ(6) -protoilludene synthases, is described. This study verifies the accuracy and utility of the predictive framework as a roadmap for the discovery of specific sesquiterpene synthases from Basidiomycota, and thus represents an important step forward in natural product discovery.

Profiling microRNAs and their targets in an important fleshy fruit: tomato (Solanum lycopersicum)

Tomato (Solanum lycopersicum) is an important and the most useful plant based diet. It is widely used for its antioxidant property. Presently, only two digits, tomato microRNAs (miRNAs) are reported in miRBase: a miRNA database. This study is aimed to profile and characterize more miRNAs and their targets in tomato. A comprehensive comparative genomic approach is applied and a total of 109 new miRNAs belonging to 106 families are identified and characterized from the tomato expressed sequence tags (ESTs). All these potential miRNAs are profiled for the first time in tomato. The profiled miRNAs are also observed with stable stem-loop structures (Precursor-miRNAs), whose length ranges from 45 to 329 nucleotides (nt) with an average of 125 nt. The mature miRNAs are found in the stem of pre-miRNAs and their length ranges from 19 to 24 nt with an average of 21 nt. Furthermore, twelve miRNAs are randomly selected and experimentally validated through RT-PCR. A total of 406 putative targets are also predicted for the newly 109 tomato miRNAs. These targets are involved in structural protein, metabolism, transcription factor, growth & development, stress related, signaling pathways, storage proteins and other vital processes. Some important proteins like; 9-cisepoxycarotenoid dioxygenase (NCED), transcription factor MYB, ATP-binding cassette transporters, terpen synthase, 14-3-3 and TIR-NBS proteins are also predicted as putative targets for tomato miRNAs. These findings improve a baseline data of miRNAs and their targets in tomato. This baseline data can be utilized to fine tune this important fleshy fruit for nutritional & antioxidant properties and also under biotic & abiotic stresses.

Isolation and characterization of terpene synthases in cotton (Gossypium hirsutum)

Cotton plants accumulate gossypol and related sesquiterpene aldehydes, which function as phytoalexins against pathogens and feeding deterrents to herbivorous insects. However, to date little is known about the biosynthesis of volatile terpenes in this crop. Herein is reported that 5 monoterpenes and 11 sesquiterpenes from extracts of a glanded cotton cultivar, Gossypium hirsutum cv. CCRI12, were detected by gas chromatography-mass spectrometry (GC-MS). By EST data mining combined with Rapid Amplification of cDNA Ends (RACE), full-length cDNAs of three terpene synthases (TPSs), GhTPS1, GhTPS2 and GhTPS3 were isolated. By in vitro assays of the recombinant proteins, it was found that GhTPS1 and GhTPS2 are sesquiterpene synthases: the former converted farnesyl pyrophosphate (FPP) into β-caryophyllene and α-humulene in a ratio of 2:1, whereas the latter produced several sesquiterpenes with guaia-1(10),11-diene as the major product. By contrast, GhTPS3 is a monoterpene synthase, which produced α-pinene, β-pinene, β-phellandrene and trace amounts of other monoterpenes from geranyl pyrophosphate (GPP). The TPS activities were also supported by Virus Induced Gene Silencing (VIGS) in the cotton plant. GhTPS1 and GhTPS3 were highly expressed in the cotton plant overall, whereas GhTPS2 was expressed only in leaves. When stimulated by mechanical wounding, Verticillium dahliae (Vde) elicitor or methyl jasmonate (MeJA), production of terpenes and expression of the corresponding synthase genes were induced. These data demonstrate that the three genes account for the biosynthesis of volatile terpenes of cotton, at least of this Upland cotton.

Cloning of a sesquiterpene synthase from Lavandula x intermedia glandular trichomes

The essential oil (EO) of Lavandula is dominated by monoterpenes, but can also contain small amounts of sesquiterpenes, depending on species and environmental conditions. For example, the sesquiterpene 9-epi-caryophyllene can make up to 8 % of the EO in a few species, including those commercially propagated for EO production. Here, we report the cloning and functional characterization of 9-epi-caryophyllene synthase (LiCPS) from the glandular trichomes of Lavandula x intermedia, cv. Grosso. The 1,617 bp open reading frame of LiCPS, which did not encode a transit peptide, was expressed in Escherichia coli and the recombinant protein purified by Ni-NTA agarose affinity chromatography. The ca. 60 kDa recombinant protein specifically converted farnesyl diphosphate to 9-epi-caryophyllene. LiCPS also produced a few monoterpenes when assayed with the monoterpene precursor geranyl diphosphate (GPP), but--unlike most monoterpene synthases--was not able to derive detectable amounts of any products from the cis isomer of GPP, neryl diphosphate. The LiCPS transcripts accumulated in developing L. x intermedia flowers and were highly enriched in glandular trichomes, but were not detected in leaves suggesting that the transcriptional expression of this gene is spatially and developmentally regulated.

Transcriptome sequencing and expression analysis of terpenoid biosynthesis genes in Litsea cubeba

BACKGROUND

Aromatic essential oils extracted from fresh fruits of Litsea cubeba (Lour.) Pers., have diverse medical and economic values. The dominant components in these essential oils are monoterpenes and sesquiterpenes. Understanding the molecular mechanisms of terpenoid biosynthesis is essential for improving the yield and quality of terpenes. However, the 40 available L. cubeba nucleotide sequences in the public databases are insufficient for studying the molecular mechanisms. Thus, high-throughput transcriptome sequencing of L. cubeba is necessary to generate large quantities of transcript sequences for the purpose of gene discovery, especially terpenoid biosynthesis related genes.

RESULTS

Using Illumina paired-end sequencing, approximately 23.5 million high-quality reads were generated. De novo assembly yielded 68,648 unigenes with an average length of 834 bp. A total of 38,439 (56%) unigenes were annotated for their functions, and 35,732 and 25,806 unigenes could be aligned to the GO and COG database, respectively. By searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG), 16,130 unigenes were assigned to 297 KEGG pathways, and 61 unigenes, which contained the mevalonate and 2-C-methyl-D-erythritol 4-phosphate pathways, could be related to terpenoid backbone biosynthesis. Of the 12,963 unigenes, 285 were annotated to the terpenoid pathways using the PlantCyc database. Additionally, 14 terpene synthase genes were identified from the transcriptome. The expression patterns of the 16 genes related to terpenoid biosynthesis were analyzed by RT-qPCR to explore their putative functions.

CONCLUSION

RNA sequencing was effective in identifying a large quantity of sequence information. To our knowledge, this study is the first exploration of the L. cubeba transcriptome, and the substantial amount of transcripts obtained will accelerate the understanding of the molecular mechanisms of essential oils biosynthesis. The results may help improve future genetic and genomics studies on the molecular mechanisms behind the chemical composition of essential oils in L. cubeba fruits.

Transcriptome analyses of a Chinese hazelnut species Corylus mandshurica

BACKGROUND

Corylus was renowned for its production of hazelnut and taxol. To understand the local adaptation of Chinese species and speed up breeding efforts in China, we analyzed the leaf transcriptome of Corylus mandshurica, which had a high tolerance to fungal infections and cold.

RESULTS

A total of 12,255,030 clean pair-end reads were generated and then assembled into 37,846 Expressed Sequence Tag (EST) sequences. During functional annotation, 26,565 ESTs were annotated with Gene Ontology (GO) terms using Blast2go and 11,056 ESTs were grouped into the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using KEGG Automatic Annotation Server (KAAS). We identified 45 ESTs that were homologous to enzymes and transcription factors responsible for taxol synthesis. The most differentiated orthologs between C. mandshurica and a European congener, C. avellana, were enriched in stress tolerance to fungal resistance and cold.

CONCLUSIONS

In this study, we detected a set of genes related to taxol synthesis in a taxol-producing angiosperm species for the first time and found a close relationship between most differentiated genes and different adaptations to fungal infection and cold in C. mandshurica and C. avellana. These findings provided tools to improve our understanding of local adaptation, genetic breeding and taxol production in hazelnut.

Identification and characterization of a novel monoterpene synthase from soybean restricted to neryl diphosphate precursor

Terpenes are important defensive compounds against herbivores and pathogens. Here, we report the identification of a new monoterpene synthase gene, GmNES, from soybean. The transcription of GmNES was up-regulated in soybean plants that were infested with cotton leafworm (Prodenia litura), mechanically wounded or treated with salicylic acid (SA). Gas chromatography-mass spectrometry (GC-MS) analysis revealed that recombinant GmNES enzyme exclusively produced nerol, generated from a newly identified substrate for monoterpene synthase: neryl diphosphate (NPP). This finding indicates that GmNES is a nerol synthase gene in soybean. Subcellular localization using GFP fusions showed that GmNES localized to the chloroplasts. Transgenic tobacco overexpressing GmNES was generated. In dual-choice assays, the GmNES-expressing tobacco lines significantly repelled cotton leafworm. In feeding tests with transgenic plants, the growth and development of cotton leafworm were significantly retarded. This study confirms the ecological role of terpenoids and provides new insights into their metabolic engineering in transgenic plants.

Bioengineering of plant (tri)terpenoids: from metabolic engineering of plants to synthetic biology in vivo and in vitro

Terpenoids constitute a large and diverse class of natural products that serve many functions in nature. Most of the tens of thousands of the discovered terpenoids are synthesized by plants, where they function as primary metabolites involved in growth and development, or as secondary metabolites that optimize the interaction between the plant and its environment. Several plant terpenoids are economically important molecules that serve many applications as pharmaceuticals, pesticides, etc. Major challenges for the commercialization of plant-derived terpenoids include their low production levels in planta and the continuous demand of industry for novel molecules with new or superior biological activities. Here, we highlight several synthetic biology methods to enhance and diversify the production of plant terpenoids, with a foresight towards triterpenoid engineering, the least engineered class of bioactive terpenoids. Increased or cheaper production of valuable triterpenoids may be obtained by 'classic' metabolic engineering of plants or by heterologous production of the compounds in other plants or microbes. Novel triterpenoid structures can be generated through combinatorial biosynthesis or directed enzyme evolution approaches. In its ultimate form, synthetic biology may lead to the production of large amounts of plant triterpenoids in in vitro systems or custom-designed artificial biological systems.

Transcriptome sequencing and differential gene expression analysis of delayed gland morphogenesis in Gossypium australe during seed germination

The genus Gossypium is a globally important crop that is used to produce textiles, oil and protein. However, gossypol, which is found in cultivated cottonseed, is toxic to humans and non-ruminant animals. Efforts have been made to breed improved cultivated cotton with lower gossypol content. The delayed gland morphogenesis trait possessed by some Australian wild cotton species may enable the widespread, direct usage of cottonseed. However, the mechanisms about the delayed gland morphogenesis are still unknown. Here, we sequenced the first Australian wild cotton species (Gossypiumaustrale) and a diploid cotton species (Gossypiumarboreum) using the Illumina Hiseq 2000 RNA-seq platform to help elucidate the mechanisms underlying gossypol synthesis and gland development. Paired-end Illumina short reads were de novo assembled into 226,184, 213,257 and 275,434 transcripts, clustering into 61,048, 47,908 and 72,985 individual clusters with N50 lengths of 1,710 bp, 1544 BP and 1,743 bp, respectively. The clustered Unigenes were searched against three public protein databases (TrEMBL, SwissProt and RefSeq) and the nucleotide and protein sequences of Gossypiumraimondii using BLASTx and BLASTn. A total of 21,987, 17,209 and 25,325 Unigenes were annotated. Of these, 18,766 (85.4%), 14,552 (84.6%) and 21,374 (84.4%) Unigenes could be assigned to GO-term classifications. We identified and analyzed 13,884 differentially expressed Unigenes by clustering and functional enrichment. Terpenoid-related biosynthesis pathways showed differentially regulated expression patterns between the two cotton species. Phylogenetic analysis of the terpene synthases family was also carried out to clarify the classifications of TPSs. RNA-seq data from two distinct cotton species provide comprehensive transcriptome annotation resources and global gene expression profiles during seed germination and gland and gossypol formation. These data may be used to further elucidate various mechanisms and help promote the usage of cottonseed.

Functional characterization of two clusters of Brachypodium distachyon UDP-glycosyltransferases encoding putative deoxynivalenol detoxification genes

Plant small-molecule UDP-glycosyltransferases (UGT) glycosylate a vast number of endogenous substances but also act in detoxification of metabolites produced by plant-pathogenic microorganisms. The ability to inactivate the Fusarium graminearum mycotoxin deoxynivalenol (DON) into DON-3-O-glucoside is crucial for resistance of cereals. We analyzed the UGT gene family of the monocot model species Brachypodium distachyon and functionally characterized two gene clusters containing putative orthologs of previously identified DON-detoxification genes from Arabidopsis thaliana and barley. Analysis of transcription showed that UGT encoded in both clusters are highly inducible by DON and expressed at much higher levels upon infection with a wild-type DON-producing F. graminearum strain compared with infection with a mutant deficient in DON production. Expression of these genes in a toxin-sensitive strain of Saccharomyces cerevisiae revealed that only two B. distachyon UGT encoded by members of a cluster of six genes homologous to the DON-inactivating barley HvUGT13248 were able to convert DON into DON-3-O-glucoside. Also, a single copy gene from Sorghum bicolor orthologous to this cluster and one of three putative orthologs of rice exhibit this ability. Seemingly, the UGT genes undergo rapid evolution and changes in copy number, making it difficult to identify orthologs with conserved substrate specificity.

Transorganellar complementation redefines the biochemical continuity of endoplasmic reticulum and chloroplasts

Tocopherols are nonpolar compounds synthesized and localized in plastids but whose genetic elimination specifically impacts fatty acid desaturation in the endoplasmic reticulum (ER), suggesting a direct interaction with ER-resident enzymes. To functionally probe for such interactions, we developed transorganellar complementation, where mutated pathway activities in one organelle are experimentally tested for substrate accessibility and complementation by active enzymes retargeted to a companion organelle. Mutations disrupting three plastid-resident activities in tocopherol and carotenoid synthesis were complemented from the ER in this fashion, demonstrating transorganellar access to at least seven nonpolar, plastid envelope-localized substrates from the lumen of the ER, likely through plastid:ER membrane interaction domains. The ability of enzymes in either organelle to access shared, nonpolar plastid metabolite pools redefines our understanding of the biochemical continuity of the ER and chloroplast with profound implications for the integration and regulation of organelle-spanning pathways that synthesize nonpolar metabolites in plants.

Genomic characterization, molecular cloning and expression analysis of two terpene synthases from Thymus caespititius (Lamiaceae)

The identification, isolation and functional characterization of two genes encoding two monoterpene synthases-γ-terpinene synthase (Tctps2) and α-terpineol synthase (Tctps5)-from three chemically distinct Thymus caespititius (Lamiaceae) genotypes were performed. Genomic exon-intron structure was also determined for both terpene synthase genes, revealing an organization with seven exons and six introns. The cDNA of Tctps2 was 2,308 bp long and had an open reading frame of 1,794 bp encoding for a protein with 598 amino acids. Tctps5 was longer, mainly due to intron sequences, and presented high intraspecific variability on the plants analyzed. It encoded for a protein of 602 amino acids from an open reading frame of 1,806 bp comprising a total of 2,507 bp genomic sequence. The amino acid sequence of these two active Tctps genes shared 74 % pairwise identity, ranging between 42 and 94 % similarity with about 50 known terpene synthases of other Lamiaceae species. Gene expression revealed a multi-product Tctps2 and Tctps5 enzymes, producing γ-terpinene and α-terpineol as major components, respectively. These enzymatic results were consistent with the monoterpene profile present in T. caespititius field plants, suggesting a transcriptional regulation in leaves. Herewith reported for the first time for this species, these two newly characterized Tctps genes improve the understanding of the molecular mechanisms of reaction responsible for terpene biosynthesis and chemical diversity found in T. caespititius.

Analysis of the expression of terpene synthase genes in relation to aroma content in two aromatic Vitis vinifera varieties

Grape (Vitis vinifera L.) flavour management in the vineyard requires knowledge of the derivation of individual flavour and aroma characteristics. Some of the most prevalent wine grape aroma constituents are terpenoids and this study represents a wide report about grape terpene synthase (TPS) gene transcript profiling in different tissues of two aromatic grapevine varieties, particularly flowers and developing berries, correlated with the accumulation patterns of free aroma compounds. All investigated genes belonging to the TPS-a and TPS-b subfamilies reached the highest expression in accordance with the peak of accumulation of the respective compounds. In the TPS-g subfamily, only one of the genes characterised for linalool synthases showed major transcript abundance in ripening berries, whereas the only geraniol synthase had a peak of expression in green berries and at the beginning of ripening, when geraniol concentration started to increase and overcome the linalool concentration. The genes identified in this study as being mainly responsible for linalool and geraniol synthesis during berry development, and the phenological phases in which they are mostly expressed, should be of interest to viticulturists and wine makers to improve decision making along the chain of production.

Evolution of a complex locus for terpene biosynthesis in solanum

Functional gene clusters, containing two or more genes encoding different enzymes for the same pathway, are sometimes observed in plant genomes, most often when the genes specify the synthesis of specialized defensive metabolites. Here, we show that a cluster of genes in tomato (Solanum lycopersicum; Solanaceae) contains genes for terpene synthases (TPSs) that specify the synthesis of monoterpenes and diterpenes from cis-prenyl diphosphates, substrates that are synthesized by enzymes encoded by cis-prenyl transferase (CPT) genes also located within the same cluster. The monoterpene synthase genes in the cluster likely evolved from a diterpene synthase gene in the cluster by duplication and divergence. In the orthologous cluster in Solanum habrochaites, a new sesquiterpene synthase gene was created by a duplication event of a monoterpene synthase followed by a localized gene conversion event directed by a diterpene synthase gene. The TPS genes in the Solanum cluster encoding cis-prenyl diphosphate-utilizing enzymes are closely related to a tobacco (Nicotiana tabacum; Solanaceae) diterpene synthase encoding Z-abienol synthase (Nt-ABS). Nt-ABS uses the substrate copal-8-ol diphosphate, which is made from the all-trans geranylgeranyl diphosphate by copal-8-ol diphosphate synthase (Nt-CPS2). The Solanum gene cluster also contains an ortholog of Nt-CPS2, but it appears to encode a nonfunctional protein. Thus, the Solanum functional gene cluster evolved by duplication and divergence of TPS genes, together with alterations in substrate specificity to utilize cis-prenyl diphosphates and through the acquisition of CPT genes.

Gene discovery of modular diterpene metabolism in nonmodel systems

Plants produce over 10,000 different diterpenes of specialized (secondary) metabolism, and fewer diterpenes of general (primary) metabolism. Specialized diterpenes may have functions in ecological interactions of plants with other organisms and also benefit humanity as pharmaceuticals, fragrances, resins, and other industrial bioproducts. Examples of high-value diterpenes are taxol and forskolin pharmaceuticals or ambroxide fragrances. Yields and purity of diterpenes obtained from natural sources or by chemical synthesis are often insufficient for large-volume or high-end applications. Improvement of agricultural or biotechnological diterpene production requires knowledge of biosynthetic genes and enzymes. However, specialized diterpene pathways are extremely diverse across the plant kingdom, and most specialized diterpenes are taxonomically restricted to a few plant species, genera, or families. Consequently, there is no single reference system to guide gene discovery and rapid annotation of specialized diterpene pathways. Functional diversification of genes and plasticity of enzyme functions of these pathways further complicate correct annotation. To address this challenge, we used a set of 10 different plant species to develop a general strategy for diterpene gene discovery in nonmodel systems. The approach combines metabolite-guided transcriptome resources, custom diterpene synthase (diTPS) and cytochrome P450 reference gene databases, phylogenies, and, as shown for select diTPSs, single and coupled enzyme assays using microbial and plant expression systems. In the 10 species, we identified 46 new diTPS candidates and over 400 putatively terpenoid-related P450s in a resource of nearly 1 million predicted transcripts of diterpene-accumulating tissues. Phylogenetic patterns of lineage-specific blooms of genes guided functional characterization.

Transcriptome resources and functional characterization of monoterpene synthases for two host species of the mountain pine beetle, lodgepole pine (Pinus contorta) and jack pine (Pinus banksiana)

BACKGROUND

The mountain pine beetle (MPB, Dendroctonus ponderosae) epidemic has affected lodgepole pine (Pinus contorta) across an area of more than 18 million hectares of pine forests in western Canada, and is a threat to the boreal jack pine (Pinus banksiana) forest. Defence of pines against MPB and associated fungal pathogens, as well as other pests, involves oleoresin monoterpenes, which are biosynthesized by families of terpene synthases (TPSs). Volatile monoterpenes also serve as host recognition cues for MPB and as precursors for MPB pheromones. The genes responsible for terpene biosynthesis in jack pine and lodgepole pine were previously unknown.

RESULTS

We report the generation and quality assessment of assembled transcriptome resources for lodgepole pine and jack pine using Sanger, Roche 454, and Illumina sequencing technologies. Assemblies revealed transcripts for approximately 20,000 - 30,000 genes from each species and assembly analyses led to the identification of candidate full-length prenyl transferase, TPS, and P450 genes of oleoresin biosynthesis. We cloned and functionally characterized, via expression of recombinant proteins in E. coli, nine different jack pine and eight different lodgepole pine mono-TPSs. The newly identified lodgepole pine and jack pine mono-TPSs include (+)-α-pinene synthases, (-)-α-pinene synthases, (-)-β-pinene synthases, (+)-3-carene synthases, and (-)-β-phellandrene synthases from each of the two species.

CONCLUSION

In the absence of genome sequences, transcriptome assemblies are important for defence gene discovery in lodgepole pine and jack pine, as demonstrated here for the terpenoid pathway genes. The product profiles of the functionally annotated mono-TPSs described here can account for the major monoterpene metabolites identified in lodgepole pine and jack pine.

Diversification of the monoterpene synthase gene family (TPSb) in Protium, a highly diverse genus of tropical trees

Plant monoterpenes are a diverse class of secondary metabolites mediating biotic and abiotic interactions with direct effects on plant fitness. To evaluate the hypothesis that monoterpene diversity is related to functional diversification after gene duplication, we reconstructed the evolutionary history of monoterpene synthases (TPSb)--the genes underlying monoterpene synthesis--in Protium, a taxonomically and chemically diverse genus of tropical trees. We isolated multiple copies of TPSb genes from chemically divergent Protium species, reconstructed the phylogeny of this gene family, used maximum-likelihood estimation of selection coefficients, and inferred residues evolving under positive selection. We found evidence for one ancient and multiple more recent duplication events giving rise to three, and potentially five, copies of TPSb genes currently present in Protium. There was evidence for adaptive evolution in one copy with a positively selected residue likely involved in protein folding and product specificity. All other copies were inferred to be evolving under a combination of stabilizing and/or relaxed selection. Although gene copy number is consistent with the extensive phenotypic diversity in monoterpenes shown in Protium, selection analyses suggest that not all copies are undergoing divergent selection consistent with a coevolutionary arms race with enemies, but instead may be under stabilizing and relaxed selection consistent with signaling or physiological stress functionality.

Functional characterization of three Coffea arabica L. monoterpene synthases: insights into the enzymatic machinery of coffee aroma

The chemical composition of the coffee beverage is extremely complex, being made up of hundreds of volatile and non-volatile compounds, many of which are generated in the thermal reactions that occur during the roasting process. However, in the raw coffee bean there are also compounds that survive roasting and are therefore extracted into the beverage. Monoterpenes are an example of this category, as their presence has been reported in the coffee flower, fruit, seed, roasted bean and in the beverage aroma. The present work describes the isolation, heterologous expression and functional characterization of three Coffea arabica cDNAs coding for monoterpene synthases. RNA was purified from C. arabica (cv. Catuai Red) flowers, seeds and fruits at 4 successive ripening stages. Degenerate primers were designed on the most conserved regions of the monoterpene synthase gene family, and then used to isolate monoterpene synthase-like sequences from the cDNA libraries. After 5'- and 3'-RACE, the complete transcripts of 4 putative C. arabica monoterpene synthases (CofarTPS) were obtained. Gene expression in different tissues and developmental stages was analysed. After heterologous expression in Escherichia coli, enzyme activity and substrate specificity were evaluated in vitro by incubation of the recombinant proteins with geranyl pyrophosphate (GPP), geranylgeranyl pyrophosphate (GGPP) and farnesyl pyrophosphate (FPP), precursors respectively of mono-, di- and sesquiterpenes. The reaction products were characterized by HS-SPME GC-MS. CofarTPS1 was classified as a limonene synthase gene, while CofarTPS2 and 3 showed lower activity with the production of linalool and β-myrcene.

Integrating the protein and metabolic engineering toolkits for next-generation chemical biosynthesis

Through microbial engineering, biosynthesis has the potential to produce thousands of chemicals used in everyday life. Metabolic engineering and synthetic biology are fields driven by the manipulation of genes, genetic regulatory systems, and enzymatic pathways for developing highly productive microbial strains. Fundamentally, it is the biochemical characteristics of the enzymes themselves that dictate flux through a biosynthetic pathway toward the product of interest. As metabolic engineers target sophisticated secondary metabolites, there has been little recognition of the reduced catalytic activity and increased substrate/product promiscuity of the corresponding enzymes compared to those of central metabolism. Thus, fine-tuning these enzymatic characteristics through protein engineering is paramount for developing high-productivity microbial strains for secondary metabolites. Here, we describe the importance of protein engineering for advancing metabolic engineering of secondary metabolism pathways. This pathway integrated enzyme optimization can enhance the collective toolkit of microbial engineering to shape the future of chemical manufacturing.

Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest

Plants produce a vast array of specialized metabolites, many of which are used as pharmaceuticals, flavors, fragrances, and other high-value fine chemicals. However, most of these compounds occur in non-model plants for which genomic sequence information is not yet available. The production of a large amount of nucleotide sequence data using next-generation technologies is now relatively fast and cost-effective, especially when using the latest Roche-454 and Illumina sequencers with enhanced base-calling accuracy. To investigate specialized metabolite biosynthesis in non-model plants we have established a data-mining framework, employing next-generation sequencing and computational algorithms, to construct and analyze the transcriptomes of 75 non-model plants that produce compounds of interest for biotechnological applications. After sequence assembly an extensive annotation approach was applied to assign functional information to over 800,000 putative transcripts. The annotation is based on direct searches against public databases, including RefSeq and InterPro. Gene Ontology (GO), Enzyme Commission (EC) annotations and associated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps are also collected. As a proof-of-concept, the selection of biosynthetic gene candidates associated with six specialized metabolic pathways is described. A web-based BLAST server has been established to allow public access to assembled transcriptome databases for all 75 plant species of the PhytoMetaSyn Project (www.phytometasyn.ca).

Metabolomic analysis of methyl jasmonate-induced triterpenoid production in the medicinal herb Centella asiatica (L.) urban

Centella asiatica is an important source of biologically active pentacyclic triterpenoids. The enhancement of the biosynthesis of the centellosides by manipulation of associated metabolic pathways is receiving much attention. Jasmonates play critical roles in plant metabolism by up-regulating the expression of genes related to secondary metabolites. Here, we investigated the effect of methyl jasmonate (MeJa) in C. asiatica through targeted metabolomic profiling of asiaticoside and madecassoside as well as their aglycones, asiatic acid and madecassic acid. Cell suspensions were treated with 0.2 mM MeJa for 2, 4 and 6 days. Liquid chromatography coupled to mass spectrometry (LC-MS) was used to explore induced changes in metabolite profiles, both qualitatively and quantitatively. Principal component analysis (PCA)-derived scores plots revealed clusters of sample replicates for control and treated samples at 2, 4 and 6 days while loading plots aided in identifying signatory biomarkers (asiatic acid and madecassic acid, as well as asiaticoside and madecassoside) that clearly demonstrate the variability between samples. In addition to increased biosynthesis of the targeted centelloids, other differential changes in the intracellular metabolite profiles reflected the response of the C. asiatica cells to the MeJa-treatment as a reprogramming of the metabolome.

Identification of genes in Thuja plicata foliar terpenoid defenses

Thuja plicata (western redcedar) is a long-lived conifer species whose foliage is rarely affected by disease or insect pests, but can be severely damaged by ungulate browsing. Deterrence to browsing correlates with high foliar levels of terpenoids, in particular the monoterpenoid α-thujone. Here, we set out to identify genes whose products may be involved in the production of α-thujone and other terpenoids in this species. First, we generated a foliar transcriptome database from which to draw candidate genes. Second, we mapped the storage of thujones and other terpenoids to foliar glands. Third, we used global expression profiling to identify more than 600 genes that are expressed at high levels in foliage with glands, but can either not be detected or are expressed at low levels in a natural variant lacking foliar glands. Fourth, we used in situ RNA hybridization to map the expression of a putative monoterpene synthase to the epithelium of glands and used enzyme assays with recombinant protein of the same gene to show that it produces sabinene, the monoterpene precursor of α-thujone. Finally, we identified candidate genes with predicted enzymatic functions for the conversion of sabinene to α-thujone. Taken together, this approach generated both general resources and detailed functional characterization in the identification of genes of foliar terpenoid biosynthesis in T. plicata.

Volatile profiles of fungi--chemotyping of species and ecological functions

Fungi emit a large spectrum of volatile organic compounds (VOCs). In the present study, we characterized and compared the odor profiles of ectomycorrhizal (EM), pathogenic and saprophytic fungal species with the aim to use these patterns as a chemotyping tool. Volatiles were collected from the headspace of eight fungal species including nine strains (four EM, three pathogens and two saprophytes) using the stir bar sorptive extraction method and analyzed by gas chromatography-mass spectrometry (GC-MS). After removal of VOCs released from the growth system, 54 VOCs were detected including 15 novel compounds not reported in fungi before. Principle component and cluster analyses revealed that fungal species differ in their odor profiles, particularly in the pattern of sesquiterpenes. The functional groups and species could be chemotyped by using their specific emission patterns. The different ecological groups could be predicted with probabilities of 90-99%, whereas for the individual species the probabilities varied between 55% and 83%. This study strongly supports the concept that the profiling of volatile compounds can be used for non-invasive identification of different functional fungal groups.

Why only some plants emit isoprene

Isoprene (2-methyl-1,3-butadiene) is emitted from many plants and it appears to have an adaptive role in protecting leaves from abiotic stress. However, only some species emit isoprene. Isoprene emission has appeared and been lost many times independently during the evolution of plants. As an example, our phylogenetic analysis shows that isoprene emission is likely ancestral within the family Fabaceae (= Leguminosae), but that it has been lost at least 16 times and secondarily gained at least 10 times through independent evolutionary events. Within the division Pteridophyta (ferns), we conservatively estimate that isoprene emissions have been gained five times and lost two times through independent evolutionary events. Within the genus Quercus (oaks), isoprene emissions have been lost from one clade, but replaced by a novel type of light-dependent monoterpene emissions that uses the same metabolic pathways and substrates as isoprene emissions. This novel type of monoterpene emissions has appeared at least twice independently within Quercus, and has been lost from 9% of the individuals within a single population of Quercus suber. Gain and loss of gene function for isoprene synthase is possible through relatively few mutations. Thus, this trait appears frequently in lineages; but, once it appears, the time available for evolutionary radiation into environments that select for the trait is short relative to the time required for mutations capable of producing a non-functional isoprene synthase gene. The high frequency of gains and losses of the trait and its heterogeneous taxonomic distribution in plants may be explained by the relatively few mutations necessary to produce or lose the isoprene synthase gene combined with the assumption that isoprene emission is advantageous in a narrow range of environments and phenotypes.

Formation of the unusual semivolatile diterpene rhizathalene by the Arabidopsis class I terpene synthase TPS08 in the root stele is involved in defense against belowground herbivory

Secondary metabolites are major constituents of plant defense against herbivore attack. Relatively little is known about the cell type-specific formation and antiherbivore activities of secondary compounds in roots despite the substantial impact of root herbivory on plant performance and fitness. Here, we describe the constitutive formation of semivolatile diterpenes called rhizathalenes by the class I terpene synthase (TPS) 08 in roots of Arabidopsis thaliana. The primary enzymatic product of TPS08, rhizathalene A, which is produced from the substrate all-trans geranylgeranyl diphosphate, represents a so far unidentified class of tricyclic diterpene carbon skeletons with an unusual tricyclic spiro-hydrindane structure. Protein targeting and administration of stable isotope precursors indicate that rhizathalenes are biosynthesized in root leucoplasts. TPS08 expression is largely localized to the root stele, suggesting a centric and gradual release of its diterpene products into the peripheral root cell layers. We demonstrate that roots of Arabidopsis tps08 mutant plants, grown aeroponically and in potting substrate, are more susceptible to herbivory by the opportunistic root herbivore fungus gnat (Bradysia spp) and suffer substantial removal of peripheral tissue at larval feeding sites. Our work provides evidence for the in vivo role of semivolatile diterpene metabolites as local antifeedants in belowground direct defense against root-feeding insects.

Differences in gene expression within a striking phenotypic mosaic Eucalyptus tree that varies in susceptibility to herbivory

BACKGROUND

Long-lived trees can accumulate mutations throughout their lifetimes that may influence biotic and abiotic interactions. For example, some Eucalyptus trees display marked variation in herbivore defence within a single canopy. These "mosaic" trees support foliage with distinct chemotypes which are differentially favoured by insect and vertebrate herbivores, resulting in susceptible and resistant branches within a single canopy. These mosaic trees provide a unique opportunity to explore the biosynthesis and genetic regulation of chemical defences in the foliage. The biosynthesis of the principal defence compounds, terpenoid-dominated essential oils, is well understood. However, the regulation of the genes involved and thus the control of phenotypic variation within a single tree canopy remains a mystery.

RESULTS

We sequenced the transcriptomes of the leaves of the two different chemotypes of a chemically mosaic Eucalyptus melliodora tree using 454 pyrosequencing technology. We used gene set enrichment analysis to identify differentially expressed transcripts and found the proportion of differentially expressed genes in the resistant and susceptible foliage similar to the transcript difference between functionally distinct tissues of the same organism, for example roots and leaves. We also investigated sequence differences in the form of single nucleotide polymorphisms and found 10 nucleotides that were different between the two branches. These are likely true SNPs and several occur in regulatory genes.

CONCLUSION

We found three lines of evidence that suggest changes to a 'master switch' can result in large scale phenotypic changes: 1. We found differential expression of terpene biosynthetic genes between the two chemotypes that could contribute to chemical variation within this plant. 2. We identified many genes that are differentially expressed between the two chemotypes, including some unique genes in each branch. These genes are involved in a variety of processes within the plant and many could contribute to the regulation of secondary metabolism, thus contributing to the chemical variation. 3. We identified 10 SNPs, some of which occur in regulatory genes that could influence secondary metabolism and thus contribute to chemical variation. Whilst this research is inherently limited by sample size, the patterns we describe could be indicative of other plant genetic mosaics.

Evolution of conifer diterpene synthases: diterpene resin acid biosynthesis in lodgepole pine and jack pine involves monofunctional and bifunctional diterpene synthases

Diterpene resin acids (DRAs) are major components of pine (Pinus spp.) oleoresin. They play critical roles in conifer defense against insects and pathogens and as a renewable resource for industrial bioproducts. The core structures of DRAs are formed in secondary (i.e. specialized) metabolism via cycloisomerization of geranylgeranyl diphosphate (GGPP) by diterpene synthases (diTPSs). Previously described gymnosperm diTPSs of DRA biosynthesis are bifunctional enzymes that catalyze the initial bicyclization of GGPP followed by rearrangement of a (+)-copalyl diphosphate intermediate at two discrete class II and class I active sites. In contrast, similar diterpenes of gibberellin primary (i.e. general) metabolism are produced by the consecutive activity of two monofunctional class II and class I diTPSs. Using high-throughput transcriptome sequencing, we discovered 11 diTPS from jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta). Three of these were orthologous to known conifer bifunctional levopimaradiene/abietadiene synthases. Surprisingly, two sets of orthologous PbdiTPSs and PcdiTPSs were monofunctional class I enzymes that lacked functional class II active sites and converted (+)-copalyl diphosphate, but not GGPP, into isopimaradiene and pimaradiene as major products. Diterpene profiles and transcriptome sequences of lodgepole pine and jack pine are consistent with roles for these diTPSs in DRA biosynthesis. The monofunctional class I diTPSs of DRA biosynthesis form a new clade within the gymnosperm-specific TPS-d3 subfamily that evolved from bifunctional diTPS rather than monofunctional enzymes (TPS-c and TPS-e) of gibberellin metabolism. Homology modeling suggested alterations in the class I active site that may have contributed to their functional specialization relative to other conifer diTPSs.

Metabolic engineering of plant monoterpenes, sesquiterpenes and diterpenes--current status and future opportunities

Terpenoids (a.k.a. isoprenoids) represent the most diverse class of natural products found in plants, with tens of thousands of reported structures. Plant-derived terpenoids have a multitude of pharmaceutical and industrial applications, but the natural resources for their extraction are often limited and, in many cases, synthetic routes are not commercially viable. Some of the most valuable terpenoids are not accumulated in model plants or crops, and genetic resources for breeding of terpenoid natural product traits are thus poorly developed. At present, metabolic engineering, either in the native producer or a heterologous host, is the only realistic alternative to improve yield and accessibility. In this review article, we will evaluate the state of the art of modulating the biosynthetic pathways for the production of mono-, sesqui- and diterpenes in plants.

Functional genomics reveals that a compact terpene synthase gene family can account for terpene volatile production in apple

Terpenes are specialized plant metabolites that act as attractants to pollinators and as defensive compounds against pathogens and herbivores, but they also play an important role in determining the quality of horticultural food products. We show that the genome of cultivated apple (Malus domestica) contains 55 putative terpene synthase (TPS) genes, of which only 10 are predicted to be functional. This low number of predicted functional TPS genes compared with other plant species was supported by the identification of only eight potentially functional TPS enzymes in apple 'Royal Gala' expressed sequence tag databases, including the previously characterized apple (E,E)-α-farnesene synthase. In planta functional characterization of these TPS enzymes showed that they could account for the majority of terpene volatiles produced in cv Royal Gala, including the sesquiterpenes germacrene-D and (E)-β-caryophyllene, the monoterpenes linalool and α-pinene, and the homoterpene (E)-4,8-dimethyl-1,3,7-nonatriene. Relative expression analysis of the TPS genes indicated that floral and vegetative tissues were the primary sites of terpene production in cv Royal Gala. However, production of cv Royal Gala floral-specific terpenes and TPS genes was observed in the fruit of some heritage apple cultivars. Our results suggest that the apple TPS gene family has been shaped by a combination of ancestral and more recent genome-wide duplication events. The relatively small number of functional enzymes suggests that the remaining terpenes produced in floral and vegetative and fruit tissues are maintained under a positive selective pressure, while the small number of terpenes found in the fruit of modern cultivars may be related to commercial breeding strategies.

The plant ontology as a tool for comparative plant anatomy and genomic analyses

The Plant Ontology (PO; http://www.plantontology.org/) is a publicly available, collaborative effort to develop and maintain a controlled, structured vocabulary ('ontology') of terms to describe plant anatomy, morphology and the stages of plant development. The goals of the PO are to link (annotate) gene expression and phenotype data to plant structures and stages of plant development, using the data model adopted by the Gene Ontology. From its original design covering only rice, maize and Arabidopsis, the scope of the PO has been expanded to include all green plants. The PO was the first multispecies anatomy ontology developed for the annotation of genes and phenotypes. Also, to our knowledge, it was one of the first biological ontologies that provides translations (via synonyms) in non-English languages such as Japanese and Spanish. As of Release #18 (July 2012), there are about 2.2 million annotations linking PO terms to >110,000 unique data objects representing genes or gene models, proteins, RNAs, germplasm and quantitative trait loci (QTLs) from 22 plant species. In this paper, we focus on the plant anatomical entity branch of the PO, describing the organizing principles, resources available to users and examples of how the PO is integrated into other plant genomics databases and web portals. We also provide two examples of comparative analyses, demonstrating how the ontology structure and PO-annotated data can be used to discover the patterns of expression of the LEAFY (LFY) and terpene synthase (TPS) gene homologs.

Non-pathogenic rhizobacteria interfere with the attraction of parasitoids to aphid-induced plant volatiles via jasmonic acid signalling

Beneficial soil-borne microbes, such as mycorrhizal fungi or rhizobacteria, can affect the interactions of plants with aboveground insects at several trophic levels. While the mechanisms of interactions with herbivorous insects, that is, the second trophic level, are starting to be understood, it remains unknown how plants mediate the interactions between soil microbes and carnivorous insects, that is, the third trophic level. Using Arabidopsis thaliana Col-0 and the aphid Myzus persicae, we evaluate here the underlying mechanisms involved in the plant-mediated interaction between the non-pathogenic rhizobacterium Pseudomonas fluorescens and the parasitoid Diaeretiella rapae, by combining ecological, chemical and molecular approaches. Rhizobacterial colonization modifies the composition of the blend of herbivore-induced plant volatiles. The volatile blend from rhizobacteria-treated aphid-infested plants is less attractive to an aphid parasitoid, in terms of both olfactory preference behaviour and oviposition, than the volatile blend from aphid-infested plants without rhizobacteria. Importantly, the effect of rhizobacteria on both the emission of herbivore-induced volatiles and parasitoid response to aphid-infested plants is lost in an Arabidopsis mutant (aos/dde2-2) that is impaired in jasmonic acid production. By modifying the blend of herbivore-induced plant volatiles that depend on the jasmonic acid-signalling pathway, root-colonizing microbes interfere with the attraction of parasitoids of leaf herbivores.

The biosynthetic origin of irregular monoterpenes in Lavandula: isolation and biochemical characterization of a novel cis-prenyl diphosphate synthase gene, lavandulyl diphosphate synthase

Lavender essential oils are constituted predominantly of regular monoterpenes, for example linalool, 1,8-cineole, and camphor. However, they also contain irregular monoterpenes including lavandulol and lavandulyl acetate. Although the majority of genes responsible for the production of regular monoterpenes in lavenders are now known, enzymes (including lavandulyl diphosphate synthase (LPPS)) catalyzing the biosynthesis of irregular monoterpenes in these plants have not been described. Here, we report the isolation and functional characterization of a novel cis-prenyl diphosphate synthase cDNA, termed Lavandula x intermedia lavandulyl diphosphate synthase (LiLPPS), through a homology-based cloning strategy. The LiLPPS ORF, encoding for a 305-amino acid long protein, was expressed in Escherichia coli, and the recombinant protein was purified by nickel-nitrilotriacetic acid affinity chromatography. The approximately 34.5-kDa bacterially produced protein specifically catalyzed the head-to-middle condensation of two dimethylallyl diphosphate units to LPP in vitro with apparent Km and kcat values of 208 ± 12 μm and 0.1 s(-1), respectively. LiLPPS is a homodimeric enzyme with a sigmoidal saturation curve and Hill coefficient of 2.7, suggesting a positive co-operative interaction among its catalytic sites. LiLPPS could be used to modulate the production of lavandulol and its derivatives in plants through metabolic engineering.

Plant P450s as versatile drivers for evolution of species-specific chemical diversity

The irreversible nature of reactions catalysed by P450s makes these enzymes landmarks in the evolution of plant metabolic pathways. Founding members of P450 families are often associated with general (i.e. primary) metabolic pathways, restricted to single copy or very few representatives, indicative of purifying selection. Recruitment of those and subsequent blooms into multi-member gene families generates genetic raw material for functional diversification, which is an inherent characteristic of specialized (i.e. secondary) metabolism. However, a growing number of highly specialized P450s from not only the CYP71 clan indicate substantial contribution of convergent and divergent evolution to the observed general and specialized metabolite diversity. We will discuss examples of how the genetic and functional diversification of plant P450s drives chemical diversity in light of plant evolution. Even though it is difficult to predict the function or substrate of a P450 based on sequence similarity, grouping with a family or subfamily in phylogenetic trees can indicate association with metabolism of particular classes of compounds. Examples will be given that focus on multi-member gene families of P450s involved in the metabolic routes of four classes of specialized metabolites: cyanogenic glucosides, glucosinolates, mono- to triterpenoids and phenylpropanoids.