Transcriptome features of stone cell development in weevil-resistant and susceptible Sitka spruce.

Stone cells are a specialized, highly lignified cell type found in both angiosperms and gymnosperms. In conifers, abundance of stone cells in the cortex provides a robust constitutive physical defense against stem feeding insects. Stone cells are a major insect-resistance trait in Sitka spruce (Picea sitchensis), occurring in dense clusters in apical shoots of trees resistant (R) to spruce weevil (Pissodes strobi) but being rare in susceptible (S) trees. To learn more about molecular mechanisms of stone cell formation in conifers, we used laser microdissection and RNA sequencing to develop cell-type-specific transcriptomes of developing stone cells from R and S trees. Using light, immunohistochemical, and fluorescence microscopy, we also visualized the deposition of cellulose, xylan, and lignin associated with stone cell development. A total of 1293 genes were differentially expressed at higher levels in developing stone cells relative to cortical parenchyma. Genes with potential roles in stone cell secondary cell wall formation (SCW) were identified and their expression evaluated over a time course of stone cell formation in R and S trees. The expression of several transcriptional regulators was associated with stone cell formation, including a NAC family transcription factor and several genes annotated as MYB transcription factors with known roles in SCW formation.

Terpene Synthases and Terpene Variation in Cannabis sativa.

Cannabis (Cannabis sativa) resin is the foundation of a multibillion dollar medicinal and recreational plant bioproducts industry. Major components of the cannabis resin are the cannabinoids and terpenes. Variations of cannabis terpene profiles contribute much to the different flavor and fragrance phenotypes that affect consumer preferences. A major problem in the cannabis industry is the lack of proper metabolic characterization of many of the existing cultivars, combined with sometimes incorrect cultivar labeling. We characterized foliar terpene profiles of plants grown from 32 seed sources and found large variation both within and between sets of plants labeled as the same cultivar. We selected five plants representing different cultivars with contrasting terpene profiles for clonal propagation, floral metabolite profiling, and trichome-specific transcriptome sequencing. Sequence analysis of these five cultivars and the reference genome of cv Purple Kush revealed a total of 33 different cannabis terpene synthase (CsTPS) genes, as well as variations of the CsTPS gene family and differential expression of terpenoid and cannabinoid pathway genes between cultivars. Our annotation of the cv Purple Kush reference genome identified 19 complete CsTPS gene models, and tandem arrays of isoprenoid and cannabinoid biosynthetic genes. An updated phylogeny of the CsTPS gene family showed three cannabis-specific clades, including a clade of sesquiterpene synthases within the TPS-b subfamily that typically contains mostly monoterpene synthases. The CsTPSs described and functionally characterized here include 13 that had not been previously characterized and that collectively explain a diverse range of cannabis terpenes.

Biosynthesis of the anti-diabetic metabolite montbretin A: glucosylation of the central intermediate mini-MbA.

Type 2 diabetes (T2D) affects over 320 million people worldwide. Healthy lifestyles, improved drugs and effective nutraceuticals are different components of a response against the growing T2D epidemic. The specialized metabolite montbretin A (MbA) is being developed for treatment of T2D and obesity due to its unique pharmacological activity as a highly effective and selective inhibitor of the human pancreatic α-amylase. MbA is an acylated flavonol glycoside found in small amounts in montbretia (Crocosmia × crocosmiiflora) corms. MbA cannot be obtained in sufficient quantities for drug development from its natural source or by chemical synthesis. To overcome these limitations through metabolic engineering, we are investigating the genes and enzymes of MbA biosynthesis. We previously reported the first three steps of MbA biosynthesis from myricetin to myricetin 3-O-(6'-O-caffeoyl)-glucosyl rhamnoside (mini-MbA). Here, we describe the sequence of reactions from mini-MbA to MbA, and the discovery and characterization of the gene and enzyme responsible for the glucosylation of mini-MbA. The UDP-dependent glucosyltransferase CcUGT3 (UGT703E1) catalyzes the 1,2-glucosylation of mini-MbA to produce myricetin 3-O-(glucosyl-6'-O-caffeoyl)-glucosyl rhamnoside. Co-expression of CcUGT3 with genes for myricetin and mini-MbA biosynthesis in Nicotiana benthamiana validated its biological function and expanded the set of genes available for metabolic engineering of MbA.

Biosynthesis of the microtubule-destabilizing diterpene pseudolaric acid B from golden larch involves an unusual diterpene synthase.

The diversity of small molecules formed via plant diterpene metabolism offers a rich source of known and potentially new biopharmaceuticals. Among these, the microtubule-destabilizing activity of pseudolaric acid B (PAB) holds promise for new anticancer agents. PAB is found, perhaps uniquely, in the coniferous tree golden larch (Pseudolarix amabilis, Pxa). Here we describe the discovery and mechanistic analysis of golden larch terpene synthase 8 (PxaTPS8), an unusual diterpene synthase (diTPS) that catalyzes the first committed step in PAB biosynthesis. Mining of the golden larch root transcriptome revealed a large TPS family, including the monofunctional class I diTPS PxaTPS8, which converts geranylgeranyl diphosphate into a previously unknown 5,7-fused bicyclic diterpene, coined "pseudolaratriene." Combined NMR and quantum chemical analysis verified the structure of pseudolaratriene, and co-occurrence with PxaTPS8 and PAB in P amabilis tissues supports the intermediacy of pseudolaratriene in PAB metabolism. Although PxaTPS8 adopts the typical three-domain structure of diTPSs, sequence phylogeny places the enzyme with two-domain TPSs of mono- and sesqui-terpene biosynthesis. Site-directed mutagenesis of PxaTPS8 revealed several catalytic residues that, together with quantum chemical calculations, suggested a substantial divergence of PxaTPS8 from other TPSs leading to a distinct carbocation-driven reaction mechanism en route to the 5,7-trans-fused bicyclic pseudolaratriene scaffold. PxaTPS8 expression in microbial and plant hosts provided proof of concept for metabolic engineering of pseudolaratriene.

Oleoresin defenses in conifers: chemical diversity, terpene synthases and limitations of oleoresin defense under climate change.

Conifers have evolved complex oleoresin terpene defenses against herbivores and pathogens. In co-evolved bark beetles, conifer terpenes also serve chemo-ecological functions as pheromone precursors, chemical barcodes for host identification, or nutrients for insect-associated microbiomes. We highlight the genomic, molecular and biochemical underpinnings of the large chemical space of conifer oleoresin terpenes and volatiles. Conifer terpenes are predominantly the products of the conifer terpene synthase (TPS) gene family. Terpene diversity is increased by cytochromes P450 of the CYP720B class. Many conifer TPS are multiproduct enzymes. Multisubstrate CYP720B enzymes catalyse multistep oxidations. We summarise known terpenoid gene functions in various different conifer species with reference to the annotated terpenoid gene space in a spruce genome. Overall, biosynthesis of terpene diversity in conifers is achieved through a system of biochemical radiation and metabolic grids. Expression of TPS and CYP720B genes can be specific to individual cell types of constitutive or traumatic resin duct systems. Induced terpenoid transcriptomes in resin duct cells lead to dynamic changes of terpene composition and quantity to fend off herbivores and pathogens. While terpenoid defenses have contributed much to the evolutionary success of conifers, under new conditions of climate change, these defences may become inconsequential against range-expanding forest pests.

Functions of stone cells and oleoresin terpenes in the conifer defense syndrome.

Conifers depend on complex defense systems against herbivores. Stone cells (SC) and oleoresin are physical and chemical defenses of Sitka spruce that have been separately studied in previous work. Weevil oviposit at the tip of the previous year's apical shoot (PYAS). We investigated interactions between weevil larvae and trees in controlled oviposition experiments with resistant (R) and susceptible (S) Sitka spruce. R trees have an abundance of SC in the PYAS cortex. SC are mostly absent in S trees. R trees and S trees also differ in the composition of oleoresin terpenes. Transcriptomes of R and S trees revealed differences in long-term weevil-induced responses. Performance of larvae was significantly reduced on R trees compared with S trees under experimental conditions that mimicked natural oviposition behavior at apical shoot tips and may be attributed to the effects of SC. In oviposition experiments designed for larvae to feed below the area of highest SC abundance, larvae showed an unusual feeding behavior and oleoresin appeared to function as the major defense. The results support a role for both SC and oleoresin terpenes and possible synergies between these traits in the defense syndrome of weevil-resistant Sitka spruce.

Discovery, Biosynthesis and Stress-Related Accumulation of Dolabradiene-Derived Defenses in Maize.

Terpenoids are a major component of maize (Zea mays) chemical defenses that mediate responses to herbivores, pathogens, and other environmental challenges. Here, we describe the biosynthesis and elicited production of a class of maize diterpenoids, named dolabralexins. Dolabralexin biosynthesis involves the sequential activity of two diterpene synthases, ENT-COPALYL DIPHOSPHATE SYNTHASE (ZmAN2) and KAURENE SYNTHASE-LIKE4 (ZmKSL4). Together, ZmAN2 and ZmKSL4 form the diterpene hydrocarbon dolabradiene. In addition, we biochemically characterized a cytochrome P450 monooxygenase, ZmCYP71Z16, which catalyzes the oxygenation of dolabradiene to yield the epoxides 15,16-epoxydolabrene (epoxydolabrene) and 3ß-hydroxy-15,16-epoxydolabrene (epoxydolabranol). The absence of dolabradiene and epoxydolabranol in Zman2 mutants under elicited conditions confirmed the in vivo biosynthetic requirement of ZmAN2. Combined mass spectrometry and NMR experiments demonstrated that much of the epoxydolabranol is further converted into 3ß,15,16-trihydroxydolabrene (trihydroxydolabrene). Metabolite profiling of field-grown maize root tissues indicated that dolabralexin biosynthesis is widespread across common maize cultivars, with trihydroxydolabrene as the predominant diterpenoid. Oxidative stress induced dolabralexin accumulation and transcript expression of ZmAN2 and ZmKSL4 in root tissues, and metabolite and transcript accumulation were up-regulated in response to elicitation with the fungal pathogens Fusarium verticillioides and Fusarium graminearum Consistently, epoxydolabranol significantly inhibited the growth of both pathogens in vitro at 10 µg mL-1, while trihydroxydolabrene-mediated inhibition was specific to Fverticillioides These findings suggest that dolabralexins have defense-related roles in maize stress interactions and expand the known chemical space of diterpenoid defenses as genetic targets for understanding and ultimately improving maize resilience.

A molecular and genomic reference system for conifer defence against insects.

Insect pests are part of natural forest ecosystems contributing to forest rejuvenation but can also cause ecological disturbance and economic losses that are expected to increase with climate change. The white pine or spruce weevil (Pissodes strobi) is a pest of conifer forests in North America. Weevil-host interactions with various spruce (Picea) species have been explored as a genomic and molecular reference system for conifer defence against insects. Interactions occur in two major phases of the insect life cycle. In the exophase, adult weevils are free-moving and display behaviour of host selection for oviposition that is affected by host traits. In the endophase, insects live within the host where mobility and development from eggs to young adults are affected by a complex system of host defences. Genetic resistance exists in several spruce species and involves synergism of constitutive and induced chemical and physical defences that comprise the conifer defence syndrome. Here, we review conifer defences that disrupt the weevil life cycle and mechanisms by which trees resist weevil attack. We highlight molecular and genomic aspects and a possible role for the weevil microbiome. Knowledge of this conifer defence system is supporting forest health strategies and tree breeding for insect resistance.

Cell-type- and tissue-specific transcriptomes of the white spruce (Picea glauca) bark unmask fine-scale spatial patterns of constitutive and induced conifer defense.

Plant defenses often involve specialized cells and tissues. In conifers, specialized cells of the bark are important for defense against insects and pathogens. Using laser microdissection, we characterized the transcriptomes of cortical resin duct cells, phenolic cells and phloem of white spruce (Picea glauca) bark under constitutive and methyl jasmonate (MeJa)-induced conditions, and we compared these transcriptomes with the transcriptome of the bark tissue complex. Overall, ~3700 bark transcripts were differentially expressed in response to MeJa. Approximately 25% of transcripts were expressed in only one cell type, revealing cell specialization at the transcriptome level. MeJa caused cell-type-specific transcriptome responses and changed the overall patterns of cell-type-specific transcript accumulation. Comparison of transcriptomes of the conifer bark tissue complex and specialized cells resolved a masking effect inherent to transcriptome analysis of complex tissues, and showed the actual cell-type-specific transcriptome signatures. Characterization of cell-type-specific transcriptomes is critical to reveal the dynamic patterns of spatial and temporal display of constitutive and induced defense systems in a complex plant tissue or organ. This was demonstrated with the improved resolution of spatially restricted expression of sets of genes of secondary metabolism in the specialized cell types.

Biosynthesis of the psychotropic plant diterpene salvinorin A: Discovery and characterization of the Salvia divinorum clerodienyl diphosphate synthase.

Salvia divinorum commonly known as diviner's sage, is an ethnomedicinal plant of the mint family (Lamiaceae). Salvia divinorum is rich in clerodane-type diterpenoids, which accumulate predominantly in leaf glandular trichomes. The main bioactive metabolite, salvinorin A, is the first non-nitrogenous natural compound known to function as an opioid-receptor agonist, and is undergoing clinical trials for potential use in treating neuropsychiatric diseases and drug addictions. We report here the discovery and functional characterization of two S. divinorum diterpene synthases (diTPSs), the ent-copalyl diphosphate (ent-CPP) synthase SdCPS1, and the clerodienyl diphosphate (CLPP) synthase SdCPS2. Mining of leaf- and trichome-specific transcriptomes revealed five diTPSs, two of which are class II diTPSs (SdCPS1-2) and three are class I enzymes (SdKSL1-3). Of the class II diTPSs, transient expression in Nicotiana benthamiana identified SdCPS1 as an ent-CPP synthase, which is prevalent in roots and, together with SdKSL1, exhibits a possible dual role in general and specialized metabolism. In vivo co-expression and in vitro assays combined with nuclear magnetic resonance (NMR) analysis identified SdCPS2 as a CLPP synthase. A role of SdCPS2 in catalyzing the committed step in salvinorin A biosynthesis is supported by its biochemical function, trichome-specific expression and absence of additional class II diTPSs in S. divinorum. Structure-guided mutagenesis revealed four catalytic residues that enabled the re-programming of SdCPS2 activity to afford four distinct products, thus advancing our understanding of how neo-functionalization events have shaped the array of different class II diTPS functions in plants, and may promote synthetic biology platforms for a broader spectrum of diterpenoid bioproducts.

Heartwood-specific transcriptome and metabolite signatures of tropical sandalwood (Santalum album) reveal the final step of (Z)-santalol fragrance biosynthesis.

Tropical sandalwood (Santalum album) produces one of the world's most highly prized fragrances, which is extracted from mature heartwood. However, in some places such as southern India, natural populations of this slow-growing tree are threatened by over-exploitation. Sandalwood oil contains four major and fragrance-defining sesquiterpenols: (Z)-α-santalol, (Z)-ß-santalol, (Z)-epi-ß-santalol and (Z)-α-exo-bergamotol. The first committed step in their biosynthesis is catalyzed by a multi-product santalene/bergamotene synthase. Sandalwood cytochromes P450 of the CYP76F sub-family were recently shown to hydroxylate santalenes and bergamotene; however, these enzymes produced mostly (E)-santalols and (E)-α-exo-bergamotol. We hypothesized that different santalene/bergamotene hydroxylases evolved in S. album to stereo-selectively produce (E)- or (Z)-sesquiterpenols, and that genes encoding (Z)-specific P450s contribute to sandalwood oil formation if co-expressed in the heartwood with upstream genes of sesquiterpene biosynthesis. This hypothesis was validated by the discovery of a heartwood-specific transcriptome signature for sesquiterpenoid biosynthesis, including highly expressed SaCYP736A167 transcripts. We characterized SaCYP736A167 as a multi-substrate P450, which stereo-selectively produces (Z)-α-santalol, (Z)-ß-santalol, (Z)-epi-ß-santalol and (Z)-α-exo-bergamotol, matching authentic sandalwood oil. This work completes the discovery of the biosynthetic enzymes of key components of sandalwood fragrance, and highlights the evolutionary diversification of stereo-selective P450s in sesquiterpenoid biosynthesis. Bioengineering of microbial systems using SaCYP736A167, combined with santalene/bergamotene synthase, has potential for development of alternative industrial production systems for sandalwood oil fragrances.

Modularity of Conifer Diterpene Resin Acid Biosynthesis: P450 Enzymes of Different CYP720B Clades Use Alternative Substrates and Converge on the Same Products.

Cytochrome P450 enzymes of the CYP720B subfamily play a central role in the biosynthesis of diterpene resin acids (DRAs), which are a major component of the conifer oleoresin defense system. CYP720Bs exist in families of up to a dozen different members in conifer genomes and fall into four different clades (I-IV). Only two CYP720B members, loblolly pine (Pinus taeda) PtCYP720B1 and Sitka spruce (Picea sitchensis) PsCYP720B4, have been characterized previously. Both are multisubstrate and multifunctional clade III enzymes, which catalyze consecutive three-step oxidations in the conversion of diterpene olefins to DRAs. These reactions resemble the sequential diterpene oxidations affording ent-kaurenoic acid from ent-kaurene in gibberellin biosynthesis. Here, we functionally characterized the CYP720B clade I enzymes CYP720B2 and CYP720B12 in three different conifer species, Sitka spruce, lodgepole pine (Pinus contorta), and jack pine (Pinus banksiana), and compared their activities with those of the clade III enzymes CYP720B1 and CYP720B4 of the same species. Unlike the clade III enzymes, clade I enzymes were ultimately found not to be active with diterpene olefins but converted the recently discovered, unstable diterpene synthase product 13-hydroxy-8(14)-abietene. Through alternative routes, CYP720B enzymes of both clades produce some of the same profiles of conifer oleoresin DRAs (abietic acid, neoabietic acid, levopimaric acid, and palustric acid), while clade III enzymes also function in the formation of pimaric acid, isopimaric acid, and sandaracopimaric acid. These results highlight the modularity of the specialized (i.e. secondary) diterpene metabolism, which produces conifer defense metabolites through variable combinations of different diterpene synthase and CYP720B enzymes.

Exploring diterpene metabolism in non-model species: transcriptome-enabled discovery and functional characterization of labda-7,13E-dienyl diphosphate synthase from Grindelia robusta.

Grindelia robusta or gumweed, is a medicinal herb of the sunflower family that forms a diverse suite of diterpenoid natural products. Its major constituents, grindelic acid and related grindelane diterpenoids accumulate in a resinous exudate covering the plants' surfaces, most prominently the unopened composite flower. Recent studies demonstrated potential pharmaceutical applications for grindelic acid and its synthetic derivatives. Mining of the previously published transcriptome of G. robusta flower tissue identified two additional diterpene synthases (diTPSs). We report the in vitro and in vivo functional characterization of an ent-kaurene synthase of general metabolism (GrTPS4) and a class II diTPS (GrTPS2) of specialized metabolism that converts geranylgeranyl diphosphate (GGPP) into labda-7,13E-dienyl diphosphate as verified by nuclear magnetic resonance (NMR) analysis. Tissue-specific transcript abundance of GrTPS2 in leaves and flowers accompanied by the presence of an endocyclic 7,13 double bond in labda-7,13E-dienyl diphosphate suggest that GrTPS2 catalyzes the first committed reaction in the biosynthesis of grindelic acid and related grindelane metabolites. With the formation of labda-7,13E-dienyl diphosphate, GrTPS2 adds an additional function to the portfolio of monofunctional class II diTPSs, which catalytically most closely resembles the bifunctional labda-7,13E-dien-15-ol synthase of the lycopod Selaginella moellendorffii. Together with a recently identified functional diTPS pair of G. robusta producing manoyl oxide, GrTPS2 lays the biosynthetic foundation of the diverse array of labdane-related diterpenoids in the genus Grindelia. Knowledge of these natural diterpenoid metabolic pathways paves the way for developing biotechnology approaches toward producing grindelic acid and related bioproducts.

The gymnosperm cytochrome P450 CYP750B1 catalyzes stereospecific monoterpene hydroxylation of (+)-sabinene in thujone biosynthesis in western redcedar.

Western redcedar (WRC; Thuja plicata) produces high amounts of oxygenated thujone monoterpenoids associated with resistance against herbivore feeding, particularly ungulate browsing. Thujones and other monoterpenoids accumulate in glandular structures in the foliage of WRC. Thujones are produced from (+)-sabinene by sabinol and sabinone. Using metabolite analysis, enzyme assays with WRC tissue extracts, cloning, and functional characterization of cytochrome P450 monooxygenases, we established that trans-sabin-3-ol but not cis-sabin-3-ol is the intermediate in thujone biosynthesis in WRC. Based on transcriptome analysis, full-length complementary DNA cloning, and characterization of expressed P450 proteins, we identified CYP750B1 and CYP76AA25 as the enzymes that catalyze the hydroxylation of (+)-sabinene to trans-sabin-3-ol. Gene-specific transcript analysis in contrasting WRC genotypes producing high and low amounts of monoterpenoids, including a glandless low-terpenoid clone, as well as assays for substrate specificity supported a biological role of CYP750B1 in α- and ß-thujone biosynthesis. This P450 belongs to the apparently gymnosperm-specific CYP750 family and is, to our knowledge, the first member of this family to be functionally characterized. In contrast, CYP76AA25 has a broader substrate spectrum, also converting the sesquiterpene farnesene and the herbicide isoproturon, and its transcript profiles are not well correlated with thujone accumulation.

Function of Sitka spruce stone cells as a physical defence against white pine weevil.

Stone cells are a physical defence of conifers against stem feeding insects such as weevils and bark beetles. In Sitka spruce, abundance of stone cells in the cortex of apical shoot tips is associated with resistance to white pine weevil. However, the mode of action by which stone cells interfere with growth and development of weevil larvae is unknown. We developed a bioassay system for testing potential effects of stone cells, which were isolated from resistant trees, on weevil larvae. Bioassays using artificial diet and controlled amounts of stone cells focused on physical defence. We evaluated the effects of stone cells on establishment of neonate larvae, mandible wear and changes in relative growth rates of third instar larvae. Establishment of neonates and relative growth rates of third instars were significantly reduced by stone cells. Stone cells appeared to be indigestible by weevil larvae. Our results suggest that stone cells affect weevil establishment and development by forming a physical feeding barrier against neonate larvae at the site of oviposition, and by reducing access to nutrients in the cortex of resistant trees, which contain an abundance of stone cells in place of a more nutrient rich tissue in susceptible trees.

Histology and cell wall biochemistry of stone cells in the physical defence of conifers against insects.

Conifers possess an array of physical and chemical defences against stem-boring insects. Stone cells provide a physical defence associated with resistance against bark beetles and weevils. In Sitka spruce (Picea sitchensis), abundance of stone cells in the cortex of apical shoots is positively correlated with resistance to white pine weevil (Pissodes strobi). We identified histological, biochemical and molecular differences in the stone cell phenotype of weevil resistant (R) or susceptible (S) Sitka spruce genotypes. R trees displayed significantly higher quantities of cortical stone cells near the apical shoot node, the primary site for weevil feeding. Lignin, cellulose, xylan and mannan were the most abundant components of stone cell secondary walls, respectively. Lignin composition of stone cells isolated from R trees contained a higher percentage of G-lignin compared with S trees. Transcript profiling revealed higher transcript abundance in the R genotype of coumarate 3-hydroxylase, a key monolignol biosynthetic gene. Developing stone cells in current year apical shoots incorporated fluorescent-tagged monolignol into the secondary cell wall, while mature stone cells of previous year apical shoots did not. Stone cell development is an ephemeral process, and fortification of shoot tips in R trees is an effective strategy against insect feeding.

Expanding the Landscape of Diterpene Structural Diversity through Stereochemically Controlled Combinatorial Biosynthesis.

Plant-derived diterpenoids serve as important pharmaceuticals, food additives, and fragrances, yet their low natural abundance and high structural complexity limits their broader industrial utilization. By mimicking the modularity of diterpene biosynthesis in plants, we constructed 51 functional combinations of class I and II diterpene synthases, 41 of which are "new-to-nature". Stereoselective biosynthesis of over 50 diterpene skeletons was demonstrated, including natural variants and novel enantiomeric or diastereomeric counterparts. Scalable biotechnological production for four industrially relevant targets was accomplished in engineered strains of Saccharomyces cerevisiae.

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.

Diterpene synthases of the biosynthetic system of medicinally active diterpenoids in Marrubium vulgare.

Marrubium vulgare (Lamiaceae) is a medicinal plant whose major bioactive compounds, marrubiin and other labdane-related furanoid diterpenoids, have potential applications as anti-diabetics, analgesics or vasorelaxants. Metabolite and transcriptome profiling of M. vulgare leaves identified five different candidate diterpene synthases (diTPSs) of the TPS-c and TPS-e/f clades. We describe the in vitro and in vivo functional characterization of the M. vulgare diTPS family. In addition to MvEKS ent-kaurene synthase of general metabolism, we identified three diTPSs of specialized metabolism: MvCPS3 (+)-copalyl diphosphate synthase, and the functional diTPS pair MvCPS1 and MvELS. In a sequential reaction, MvCPS1 and MvELS produce a unique oxygenated diterpene scaffold 9,13-epoxy-labd-14-ene en route to marrubiin and an array of related compounds. In contrast with previously known diTPSs that introduce a hydroxyl group at carbon C-8 of the labdane backbone, the MvCPS1-catalyzed reaction proceeds via oxygenation of an intermediate carbocation at C-9, yielding the bicyclic peregrinol diphosphate. MvELS belongs to a subgroup of the diTPS TPS-e/f clade with unusual ßα-domain architecture. MvELS is active in vitro and in vivo with three different prenyl diphosphate substrates forming the marrubiin precursor 9,13-epoxy-labd-14-ene, as identified by nuclear magnetic resonance (NMR) analysis, manoyl oxide and miltiradiene. MvELS fills a central position in the biosynthetic system that forms the foundation for the diverse repertoire of Marrubium diterpenoids. Co-expression of MvCPS1 and MvELS in engineered E. coli and Nicotiana benthamiana offers opportunities for producing precursors for an array of biologically active diterpenoids.

Comparative genomics of the pine pathogens and beetle symbionts in the genus Grosmannia.

Studies on beetle/tree fungal symbionts typically characterize the ecological and geographic distributions of the fungal populations. There is limited understanding of the genome-wide evolutionary processes that act within and between species as such fungi adapt to different environments, leading to physiological differences and reproductive isolation. Here, we assess genomic evidence for such evolutionary processes by extending our recent work on Grosmannia clavigera, which is vectored by the mountain pine beetle and jeffrey pine beetle. We report the genome sequences of an additional 11 G. clavigera (Gc) sensu lato strains from the two known sibling species, Grosmannia sp. (Gs) and Gc. The 12 fungal genomes are structurally similar, showing large-scale synteny within and between species. We identified 103,430 single-nucleotide variations that separated the Grosmannia strains into divergent Gs and Gc clades, and further divided each of these clades into two subclades, one of which may represent an additional species. Comparing variable genes between these lineages, we identified truncated genes and potential pseudogenes, as well as seven genes that show evidence of positive selection. As these variable genes are involved in secondary metabolism and in detoxifying or utilizing host-tree defense chemicals (e.g., polyketide synthases, oxidoreductases, and mono-oxygenases), their variants may reflect adaptation to the specific chemistries of the host trees Pinus contorta, P. ponderosa, and P. jeffreyi. This work provides a comprehensive resource for developing informative markers for landscape population genomics of these ecologically and economically important fungi, and an approach that could be extended to other beetle-tree-associated fungi.