Evolution of diterpene metabolism: Sitka spruce CYP720B4 catalyzes multiple oxidations in resin acid biosynthesis of conifer defense against insects

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.

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.

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.

High-level diterpene production by transient expression in Nicotiana benthamiana

BACKGROUND

Characterization of plant terpene synthases is typically done by production of recombinant enzymes in Escherichia coli. This is often difficult due to solubility and codon usage issues. Furthermore, plant terpene synthases which are targeted to the plastids, such as diterpene synthases, have to be shortened in a more or less empirical approach to improve expression. We report here an optimized Agrobacterium-mediated transient expression assay in Nicotiana benthamiana for plant diterpene synthase expression and product analysis.

RESULTS

Agrobacterium-mediated transient expression of plant diterpene synthases in N. benthamiana led to the accumulation of diterpenes within 3 days of infiltration and with a maximum at 5 days. Over 50% of the products were exported onto the leaf surface, thus considerably facilitating the analysis by reducing the complexity of the extracts. The robustness of the method was tested by expressing three different plant enzymes, cembratrien-ol synthase from Nicotiana sylvestris, casbene synthase from Ricinus communis and levopimaradiene synthase from Gingko biloba. Furthermore, co-expression of a 1-deoxy-D-xylulose-5-phosphate synthase from tomato and a geranylgeranyl diphosphate synthase from tobacco led to a 3.5-fold increase in the amount of cembratrien-ol produced, with maximum yields reaching 2500 ng/cm2.

CONCLUSION

With this optimized method for diterpene synthase expression and product analysis, a single infiltrated leaf of N. benthamiana would be sufficient to produce quantities required for the structure elucidation of unknown diterpenes. The method will also be of general use for gene function discovery, pathway reconstitution and metabolic engineering of diterpenoid biosynthesis in plants.

Picking sides: distinct roles for CYP76M6 and CYP76M8 in rice oryzalexin biosynthesis

Natural products biosynthesis often requires the action of multiple CYPs (cytochromes P450), whose ability to introduce oxygen, increasing solubility, is critical for imparting biological activity. In previous investigations of rice diterpenoid biosynthesis, we characterized CYPs that catalyse alternative hydroxylation of ent-sandaracopimaradiene, the precursor to the rice oryzalexin antibiotic phytoalexins. In particular, CYP76M5, CYP76M6 and CYP76M8 were all shown to carry out C-7β hydroxylation, whereas CYP701A8 catalyses C-3α hydroxylation, with oxy groups found at both positions in oryzalexins A-D, suggesting that these may act consecutively in oryzalexin biosynthesis. In the present paper, we report that, although CYP701A8 only poorly reacts with 7β-hydroxy-ent-sandaracopimaradiene, CYP76M6 and CYP76M8 readily react with 3α-hydroxy-ent-sandaracopimaradiene. Notably, their activity yields distinct products, resulting from hydroxylation at C-9β by CYP76M6 or C-7β by CYP76M8, on different sides of the core tricyclic ring structure. Thus CYP76M6 and CYP76M8 have distinct non-redundant roles in orzyalexin biosynthesis. Moreover, the resulting 3α,7β- and 3α,9β-diols correspond to oryzalexins D and E respectively. Accordingly, the results of the present study complete the functional identification of the biosynthetic pathway underlying the production of these bioactive phytoalexins. In addition, the altered regiochemistry catalysed by CYP76M6 following C-3α hydroxylation has some implications for its active-site configuration, offering further molecular insight.

Association genetics of oleoresin flow in loblolly pine: discovering genes and predicting phenotype for improved resistance to bark beetles and bioenergy potential

Rapidly enhancing oleoresin production in conifer stems through genomic selection and genetic engineering may increase resistance to bark beetles and terpenoid yield for liquid biofuels. We integrated association genetic and genomic prediction analyses of oleoresin flow (g 24 h(-1)) using 4854 single nucleotide polymorphisms (SNPs) in expressed genes within a pedigreed population of loblolly pine (Pinus taeda) that was clonally replicated at three sites in the southeastern United States. Additive genetic variation in oleoresin flow (h(2) ≈ 0.12-0.30) was strongly correlated between years in which precipitation varied (r(a) ≈ 0.95), while the genetic correlation between sites declined from 0.8 to 0.37 with increasing differences in soil and climate among sites. A total of 231 SNPs were significantly associated with oleoresin flow, of which 81% were specific to individual sites. SNPs in sequences similar to ethylene signaling proteins, ABC transporters, and diterpenoid hydroxylases were associated with oleoresin flow across sites. Despite this complex genetic architecture, we developed a genomic prediction model to accelerate breeding for enhanced oleoresin flow that is robust to environmental variation. Results imply that breeding could increase oleoresin flow 1.5- to 2.4-fold in one generation.

CYP76AH1 catalyzes turnover of miltiradiene in tanshinones biosynthesis and enables heterologous production of ferruginol in yeasts

Cytochrome P450 enzymes (CYPs) play major roles in generating highly functionalized terpenoids, but identifying the exact biotransformation step(s) catalyzed by plant CYP in terpenoid biosynthesis is extremely challenging. Tanshinones are abietane-type norditerpenoid naphthoquinones that are the main lipophilic bioactive components of the Chinese medicinal herb danshen (Salvia miltiorrhiza). Whereas the diterpene synthases responsible for the conversion of (E,E,E)-geranylgeranyl diphosphate into the abietane miltiradiene, a potential precursor to tanshinones, have been recently described, molecular characterization of further transformation of miltiradiene remains unavailable. Here we report stable-isotope labeling results that demonstrate the intermediacy of miltiradiene in tanshinone biosynthesis. We further use a next-generation sequencing approach to identify six candidate CYP genes being coregulated with the diterpene synthase genes in both the rhizome and danshen hairy roots, and demonstrate that one of these, CYP76AH1, catalyzes a unique four-electron oxidation cascade on miltiradiene to produce ferruginol both in vitro and in vivo. We then build upon the previous establishment of miltiradiene production in Saccharomyces cerevisiae, with incorporation of CYP76AH1 and phyto-CYP reductase genes leading to heterologous production of ferruginol at 10.5 mg/L. As ferruginol has been found in many plants including danshen, the results and the approaches that were described here provide a solid foundation to further elucidate the biosynthesis of tanshinones and related diterpenoids. Moreover, these results should facilitate the construction of microbial cell factories for the production of phytoterpenoids.

Bacteria associated with a tree-killing insect reduce concentrations of plant defense compounds

Bark beetles encounter a diverse array of constitutive and rapidly induced terpenes when attempting to colonize living conifers. Concentrations of these compounds at entry sites can rapidly reach levels toxic to beetles, their brood, and fungal symbionts. Large numbers of beetles can overwhelm tree defenses via pheromone-mediated mass attacks, but the mechanisms are poorly understood. We show that bacteria associated with mountain pine beetles can metabolize monoterpenes and diterpene acids. The abilities of different symbionts to reduce concentrations of different terpenes appear complementary. Serratia reduced concentrations of all monoterpenes applied to media by 55-75 %, except for α-pinene. Beetle-associated Rahnella reduced (-)- and (+)-α-pinene by 40 % and 45 %, respectively. Serratia and Brevundimonas reduced diterpene abietic acid levels by 100 % at low concentrations. However, high concentrations exhausted this ability, suggesting that opposing rates of bacterial metabolism and plant induction of terpenes are critical. The two major fungal symbionts of mountain pine beetle, Grosmannia clavigera and Ophiostoma montium were highly susceptible to abietic acid. Grosmannia clavigera did not reduce total monoterpene concentrations in lodgepole pine turpentine. We propose the ability of bark beetles to exert landscape-scale impacts may arise partly from micro-scale processes driven by bacterial symbionts.

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.

Plant chemical defense: at what cost?

Plants are sessile organisms and dependent on deployment of secondary metabolites for their response to biotic and abiotic challenges. A trade-off is envisioned between resources allocated to growth, development, and reproduction and to the biosynthesis, storage, and maintenance of secondary metabolites. However, increasing evidence suggests that secondary metabolites serve auxiliary roles, including functions associated with primary metabolism. In this opinion article, we examine how the costs of plant chemical defense can be offset by multifunctional biosynthesis and the optimization of primary metabolism. These additional benefits may negate the trade-off between primary and secondary metabolism, and provide plants with an innate plasticity required for growth, development, and interactions with their environment.

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.

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.

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.

Discovery and functional characterization of two diterpene synthases for sclareol biosynthesis in Salvia sclarea (L.) and their relevance for perfume manufacture

BACKGROUND

Sclareol is a diterpene natural product of high value for the fragrance industry. Its labdane carbon skeleton and its two hydroxyl groups also make it a valued starting material for semisynthesis of numerous commercial substances, including production of Ambrox® and related ambergris substitutes used in the formulation of high end perfumes. Most of the commercially-produced sclareol is derived from cultivated clary sage (Salvia sclarea) and extraction of the plant material. In clary sage, sclareol mainly accumulates in essential oil-producing trichomes that densely cover flower calices. Manool also is a minor diterpene of this species and the main diterpene of related Salvia species.

RESULTS

Based on previous general knowledge of diterpene biosynthesis in angiosperms, and based on mining of our recently published transcriptome database obtained by deep 454-sequencing of cDNA from clary sage calices, we cloned and functionally characterized two new diterpene synthase (diTPS) enzymes for the complete biosynthesis of sclareol in clary sage. A class II diTPS (SsLPPS) produced labda-13-en-8-ol diphosphate as major product from geranylgeranyl diphosphate (GGPP) with some minor quantities of its non-hydroxylated analogue, (9 S, 10 S)-copalyl diphosphate. A class I diTPS (SsSS) then transformed these intermediates into sclareol and manool, respectively. The production of sclareol was reconstructed in vitro by combining the two recombinant diTPS enzymes with the GGPP starting substrate and in vivo by co-expression of the two proteins in yeast (Saccharomyces cerevisiae). Tobacco-based transient expression assays of green fluorescent protein-fusion constructs revealed that both enzymes possess an N-terminal signal sequence that actively targets SsLPPS and SsSS to the chloroplast, a major site of GGPP and diterpene production in plants.

CONCLUSIONS

SsLPPS and SsSS are two monofunctional diTPSs which, together, produce the diterpenoid specialized metabolite sclareol in a two-step process. They represent two of the first characterized hydroxylating diTPSs in angiosperms and generate the dihydroxylated labdane sclareol without requirement for additional enzymatic oxidation by activities such as cytochrome P450 monoxygenases. Yeast-based production of sclareol by co-expresssion of SsLPPS and SsSS was efficient enough to warrant the development and use of such technology for the biotechnological production of scareol and other oxygenated diterpenes.

Bifunctional cis-abienol synthase from Abies balsamea discovered by transcriptome sequencing and its implications for diterpenoid fragrance production

The labdanoid diterpene alcohol cis-abienol is a major component of the aromatic oleoresin of balsam fir (Abies balsamea) and serves as a valuable bioproduct material for the fragrance industry. Using high-throughput 454 transcriptome sequencing and metabolite profiling of balsam fir bark tissue, we identified candidate diterpene synthase sequences for full-length cDNA cloning and functional characterization. We discovered a bifunctional class I/II cis-abienol synthase (AbCAS), along with the paralogous levopimaradiene/abietadiene synthase and isopimaradiene synthase, all of which are members of the gymnosperm-specific TPS-d subfamily. The AbCAS-catalyzed formation of cis-abienol proceeds via cyclization and hydroxylation at carbon C-8 of a postulated carbocation intermediate in the class II active site, followed by cleavage of the diphosphate group and termination of the reaction sequence without further cyclization in the class I active site. This reaction mechanism is distinct from that of synthases of the isopimaradiene- or levopimaradiene/abietadiene synthase type, which employ deprotonation reactions in the class II active site and secondary cyclizations in the class I active site, leading to tricyclic diterpenes. Comparative homology modeling suggested the active site residues Asp-348, Leu-617, Phe-696, and Gly-723 as potentially important for the specificity of AbCAS. As a class I/II bifunctional enzyme, AbCAS is a promising target for metabolic engineering of cis-abienol production.