Plasticity and evolution of (+)-3-carene synthase and (-)-sabinene synthase functions of a sitka spruce monoterpene synthase gene family associated with weevil resistance.

The monoterpene (+)-3-carene is associated with resistance of Sitka spruce against white pine weevil, a major North American forest insect pest of pine and spruce. High and low levels of (+)-3-carene in, respectively, resistant and susceptible Sitka spruce genotypes are due to variation of (+)-3-carene synthase gene copy number, transcript and protein expression levels, enzyme product profiles, and enzyme catalytic efficiency. A family of multiproduct (+)-3-carene synthase-like genes of Sitka spruce include the three (+)-3-carene synthases, PsTPS-3car1, PsTPS-3car2, PsTPS-3car3, and the (-)-sabinene synthase PsTPS-sab. Of these, PsTPS-3car2 is responsible for the relatively higher levels of (+)-3-carene in weevil-resistant trees. Here, we identified features of the PsTPS-3car1, PsTPS-3car2, PsTPS-3car3, and PsTPS-sab proteins that determine different product profiles. A series of domain swap and site-directed mutations, supported by structural comparisons, identified the amino acid in position 596 as critical for product profiles dominated by (+)-3-carene in PsTPS-3car1, PsTPS-3car2, and PsTPS-3car3, or (-)-sabinene in PsTPS-sab. A leucine in this position promotes formation of (+)-3-carene, whereas phenylalanine promotes (-)-sabinene. Homology modeling predicts that position 596 directs product profiles through differential stabilization of the reaction intermediate. Kinetic analysis revealed position 596 also plays a role in catalytic efficiency. Mutations of position 596 with different side chain properties resulted in a series of enzymes with different product profiles, further highlighting the inherent plasticity and potential for evolution of alternative product profiles of these monoterpene synthases of conifer defense against insects.

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.

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.

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.

An integrated genomic, proteomic and biochemical analysis of (+)-3-carene biosynthesis in Sitka spruce (Picea sitchensis) genotypes that are resistant or susceptible to white pine weevil.

Conifers are extremely long-lived plants that have evolved complex chemical defenses in the form of oleoresin terpenoids to resist attack from pathogens and herbivores. In these species, terpenoid diversity is determined by the size and composition of the terpene synthase (TPS) gene family and the single- and multi-product profiles of these enzymes. The monoterpene (+)-3-carene is associated with resistance of Sitka spruce (Picea sitchensis) to white pine weevil (Pissodes strobi). We used a combined genomic, proteomic and biochemical approach to analyze the (+)-3-carene phenotype in two contrasting Sitka spruce genotypes. Resistant trees produced significantly higher levels of (+)-3-carene than susceptible trees, in which only trace amounts were detected. Biosynthesis of (+)-3-carene is controlled, at the genome level, by a small family of closely related (+)-3-carene synthase (PsTPS-3car) genes (82-95% amino acid sequence identity). Transcript profiling identified one PsTPS-3car gene (PsTPS-3car1) that is expressed in both genotypes, one gene (PsTPS-3car2) that is expressed only in resistant trees, and one gene (PsTPS-3car3) that is expressed only in susceptible trees. The PsTPS-3car2 gene was not detected in genomic DNA of susceptible trees. Target-specific selected reaction monitoring confirmed this pattern of differential expression of members of the PsTPS-3car family at the proteome level. Kinetic characterization of the recombinant PsTPS-3car enzymes identified differences in the activities of PsTPS-3car2 and PsTPS-3car3 as a factor contributing to the different (+)-3-carene profiles of resistant and susceptible trees. In conclusion, variation of the (+)-3-carene phenotype is controlled by copy number variation of PsTPS-3car genes, variation of gene and protein expression, and variation in catalytic efficiencies.