Functional Characterization of ent-Copalyl Diphosphate Synthase and Kaurene Synthase Genes from Coffea arabica L

The biochemical profile of coffee beans translates directly into quality traits, nutraceutical and health promoting properties of the coffee beverage. Ent-kaurene is the ubiquitous precursor for gibberellin biosynthesis in plants, but it also serves as an intermediate in specialized (i.e., secondary) diterpenoid metabolism that leads to a diversity of more than 1,000 different metabolites. Nutraceutical effects on human health attributed to diterpenes include antioxidant, anticarcinogenic, and anti-inflammatory properties. Cafestol (CAF) and kahweol (KAH) are two diterpenes found exclusively in the Coffea genus. Our objective was to identify and functionally characterize genes involved in the central step of ent-kaurene production. We identified 17 putative terpene synthase genes in the transcriptome of Coffea arabica. Two ent-copalyl diphosphate synthase (CaCPS) and three kaurene synthase (CaKS) were selected and manually annotated. Transcript expression profiles of CaCPS1 and CaKS3 best matched the CAF and KAH metabolite profiles in different tissues. CaCPS1 and CaKS3 proteins were heterologously expressed and functionally characterized. CaCPS1 catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to ent-copalyl diphosphate (ent-CPP), which is converted to ent-kaurene by CaKS3. Knowledge about the central steps of diterpene formation in coffee provides a foundation for future characterization of the subsequent enzymes involved in CAF and KAH biosynthesis.

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.

The genome of the forest insect pest Pissodes strobi reveals genome expansion and evidence of a Wolbachia endosymbiont

The highly diverse insect family of true weevils, Curculionidae, includes many agricultural and forest pests. Pissodes strobi, commonly known as the spruce weevil or white pine weevil, is a major pest of spruce and pine forests in North America. Pissodes strobi larvae feed on the apical shoots of young trees, causing stunted growth and can destroy regenerating spruce or pine forests. Here, we describe the nuclear and mitochondrial Pissodes strobi genomes and their annotations, as well as the genome of an apparent Wolbachia endosymbiont. We report a substantial expansion of the weevil nuclear genome, relative to other Curculionidae species, possibly driven by an abundance of class II DNA transposons. The endosymbiont observed belongs to a group (supergroup A) of Wolbachia species that generally form parasitic relationships with their arthropod host.

Constitutive and insect-induced transcriptomes of weevil-resistant and susceptible Sitka spruce

Spruce weevil (Pissodes strobi) is a significant pest of regenerating spruce (Picea) and pine (Pinus) forests in North America. Weevil larvae feed in the bark, phloem, cambium, and outer xylem of apical shoots, causing stunted growth or mortality of young trees. We identified and characterized constitutive and weevil-induced patterns of Sitka spruce (Picea sitchensis) transcriptomes in weevil-resistant (R) and susceptible (S) trees using RNA sequencing (RNA-seq) and differential expression (DE) analyses. We developed a statistical model for the analysis of RNA-seq data from treatment experiments with a 2 × 3 factorial design to differentiate insect-induced responses from the effects of mechanical damage. Across the different comparisons, we identified two major transcriptome contrasts: A large set of genes that was constitutively DE between R and S trees, and another set of genes that was DE in weevil-induced S-trees. The constitutive transcriptome unique to R trees appeared to be attuned to defense, while the constitutive transcriptome unique to S trees was enriched for growth-related transcripts. Notably, a set of transcripts annotated as "fungal" was detected consistently in the transcriptomes. Fungal transcripts were identified as DE in the comparison of R and S trees and in the weevil-affected DE transcriptome of S trees, suggesting a potential microbiome role in this conifer-insect interaction.

Genetic elucidation of interconnected antibiotic pathways mediating maize innate immunity

Specialized metabolites constitute key layers of immunity that underlie disease resistance in crops; however, challenges in resolving pathways limit our understanding of the functions and applications of these metabolites. In maize (Zea mays), the inducible accumulation of acidic terpenoids is increasingly considered to be a defence mechanism that contributes to disease resistance. Here, to understand maize antibiotic biosynthesis, we integrated association mapping, pan-genome multi-omic correlations, enzyme structure-function studies and targeted mutagenesis. We define ten genes in three zealexin (Zx) gene clusters that encode four sesquiterpene synthases and six cytochrome P450 proteins that collectively drive the production of diverse antibiotic cocktails. Quadruple mutants in which the ability to produce zealexins (ZXs) is blocked exhibit a broad-spectrum loss of disease resistance. Genetic redundancies ensuring pathway resiliency to single null mutations are combined with enzyme substrate promiscuity, creating a biosynthetic hourglass pathway that uses diverse substrates and in vivo combinatorial chemistry to yield complex antibiotic blends. The elucidated genetic basis of biochemical phenotypes that underlie disease resistance demonstrates a predominant maize defence pathway and informs innovative strategies for transferring chemical immunity between crops.

Synthetic Biology of Cannabinoids and Cannabinoid Glucosides in Nicotiana benthamiana and Saccharomyces cerevisiae

Phytocannabinoids are a group of plant-derived metabolites that display a wide range of psychoactive as well as health-promoting effects. The production of pharmaceutically relevant cannabinoids relies on extraction and purification from cannabis (Cannabis sativa) plants yielding the major constituents, Δ9-tetrahydrocannabinol and cannabidiol. Heterologous biosynthesis of cannabinoids in Nicotiana benthamiana or Saccharomyces cerevisiae may provide cost-efficient and rapid future production platforms to acquire pure and high quantities of both the major and the rare cannabinoids as well as novel derivatives. Here, we used a meta-transcriptomic analysis of cannabis to identify genes for aromatic prenyltransferases of the UbiA superfamily and chalcone isomerase-like (CHIL) proteins. Among the aromatic prenyltransferases, CsaPT4 showed CBGAS activity in both N. benthamiana and S. cerevisiae. Coexpression of selected CsaPT pairs and of CHIL proteins encoding genes with CsaPT4 did not affect CBGAS catalytic efficiency. In a screen of different plant UDP-glycosyltransferases, Stevia rebaudiana SrUGT71E1 and Oryza sativa OsUGT5 were found to glucosylate olivetolic acid, cannabigerolic acid, and Δ9-tetrahydrocannabinolic acid. Metabolic engineering of N. benthamiana for production of cannabinoids revealed intrinsic glucosylation of olivetolic acid and cannabigerolic acid. S. cerevisiae was engineered to produce olivetolic acid glucoside and cannabigerolic acid glucoside.

Effects of forced taxonomic transitions on metabolic composition and function in microbial microcosms

Surveys of microbial systems indicate that in many situations taxonomy and function may constitute largely independent ('decoupled') axes of variation. However, this decoupling is rarely explicitly tested experimentally, partly because it is hard to directly induce taxonomic variation without affecting functional composition. Here we experimentally evaluate this paradigm using microcosms resembling lake sediments and subjected to two different levels of salinity (0 and 19) and otherwise similar environmental conditions. We used DNA sequencing for taxonomic and functional profiling of bacteria and archaea and physicochemical measurements to monitor metabolic function, over 13 months. We found that the taxonomic composition of the saline systems gradually but strongly diverged from the fresh systems. In contrast, the metabolic composition (in terms of proportions of various genes) remained nearly identical across treatments and over time. Oxygen consumption rates and methane concentrations were substantially lower in the saline treatment, however, their similarity either increased (for oxygen) or did not change significantly (for methane) between the first and last sampling time, indicating that the lower metabolic activity in the saline treatments was directly and immediately caused by salinity rather than the gradual taxonomic divergence. Our experiment demonstrates that strong taxonomic shifts need not directly affect metabolic rates.

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.

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.

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.

Multiple genes recruited from hormone pathways partition maize diterpenoid defences

Duplication and divergence of primary pathway genes underlie the evolution of plant specialized metabolism; however, mechanisms partitioning parallel hormone and defence pathways are often speculative. For example, the primary pathway intermediate ent-kaurene is essential for gibberellin biosynthesis and is also a proposed precursor for maize antibiotics. By integrating transcriptional coregulation patterns, genome-wide association studies, combinatorial enzyme assays, proteomics and targeted mutant analyses, we show that maize kauralexin biosynthesis proceeds via the positional isomer ent-isokaurene formed by a diterpene synthase pair recruited from gibberellin metabolism. The oxygenation and subsequent desaturation of ent-isokaurene by three promiscuous cytochrome P450s and a new steroid 5α reductase indirectly yields predominant ent-kaurene-associated antibiotics required for Fusarium stalk rot resistance. The divergence and differential expression of pathway branches derived from multiple duplicated hormone-metabolic genes minimizes dysregulation of primary metabolism via the circuitous biosynthesis of ent-kaurene-related antibiotics without the production of growth hormone precursors during defence.

Terpenes in Cannabis sativa - From plant genome to humans

Cannabis sativa (cannabis) produces a resin that is valued for its psychoactive and medicinal properties. Despite being the foundation of a multi-billion dollar global industry, scientific knowledge and research on cannabis is lagging behind compared to other high-value crops. This is largely due to legal restrictions that have prevented many researchers from studying cannabis, its products, and their effects in humans. Cannabis resin contains hundreds of different terpene and cannabinoid metabolites. Many of these metabolites have not been conclusively identified. Our understanding of the genomic and biosynthetic systems of these metabolites in cannabis, and the factors that affect their variability, is rudimentary. As a consequence, there is concern about lack of consistency with regard to the terpene and cannabinoid composition of different cannabis 'strains'. Likewise, claims of some of the medicinal properties attributed to cannabis metabolites would benefit from thorough scientific validation.

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.

Contribution of isopentenyl phosphate to plant terpenoid metabolism

Plant genomes encode isopentenyl phosphate kinases (IPKs) that reactivate isopentenyl phosphate (IP) via ATP-dependent phosphorylation, forming the primary metabolite isopentenyl diphosphate (IPP) used generally for isoprenoid/terpenoid biosynthesis. Therefore, the existence of IPKs in plants raises unanswered questions concerning the origin and regulatory roles of IP in plant terpenoid metabolism. Here, we provide genetic and biochemical evidence showing that IP forms during specific dephosphorylation of IPP catalysed by a subset of Nudix superfamily hydrolases. Increasing metabolically available IP by overexpression of a bacterial phosphomevalonate decarboxylase (MPD) in Nicotiana tabacum resulted in significant enhancement in both monoterpene and sesquiterpene production. These results indicate that perturbing IP metabolism results in measurable changes in terpene products derived from both the methylerythritol phosphate (MEP) and mevalonate (MVA) pathways. Moreover, the unpredicted peroxisomal localization of bacterial MPD led us to discover that the step catalysed by phosphomevalonate kinase (PMK) imposes a hidden constraint on flux through the classical MVA pathway. These complementary findings fundamentally alter conventional views of metabolic regulation of terpenoid metabolism in plants and provide new metabolic engineering targets for the production of high-value terpenes in plants.

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.

An extended model of heartwood secondary metabolism informed by functional genomics

The development of heartwood (HW) and the associated accumulation of secondary metabolites, which are also known as 'specialized metabolites' or 'extractives', is an important feature of tree biology. Heartwood development can affect tree health with broader implications for forest health. Heartwood development also defines a variety of wood quality traits that are important in the forest industry such as durability and colour of wood products. In the bioproducts industry, HW provides a source of high-value small molecules such as fragrances and antimicrobials. The HW properties of decay resistance in living trees, durability and colour of wood products, and small molecule bioproducts are largely defined by secondary metabolites, the biosynthesis of which appears to be activated during the onset of HW formation. Traditionally, it is thought that HW formation involves a spike in the activity of secondary metabolism in parenchyma cells in a transition zone between sapwood and HW, followed by programmed cell-death. The resulting HW tissue is thought to consist entirely of dead cells. Here, we discuss a variation of existing models of HW formation, based on the recent discovery of HW-specific transcriptome signatures of terpenoid biosynthesis in sandalwood (Santalum album L.) that invokes the activity of living cells in HW.

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.

Genomics-Based Discovery of Plant Genes for Synthetic Biology of Terpenoid Fragrances: A Case Study in Sandalwood oil Biosynthesis

Terpenoid fragrances are powerful mediators of ecological interactions in nature and have a long history of traditional and modern industrial applications. Plants produce a great diversity of fragrant terpenoid metabolites, which make them a superb source of biosynthetic genes and enzymes. Advances in fragrance gene discovery have enabled new approaches in synthetic biology of high-value speciality molecules toward applications in the fragrance and flavor, food and beverage, cosmetics, and other industries. Rapid developments in transcriptome and genome sequencing of nonmodel plant species have accelerated the discovery of fragrance biosynthetic pathways. In parallel, advances in metabolic engineering of microbial and plant systems have established platforms for synthetic biology applications of some of the thousands of plant genes that underlie fragrance diversity. While many fragrance molecules (eg, simple monoterpenes) are abundant in readily renewable plant materials, some highly valuable fragrant terpenoids (eg, santalols, ambroxides) are rare in nature and interesting targets for synthetic biology. As a representative example for genomics/transcriptomics enabled gene and enzyme discovery, we describe a strategy used successfully for elucidation of a complete fragrance biosynthetic pathway in sandalwood (Santalum album) and its reconstruction in yeast (Saccharomyces cerevisiae). We address questions related to the discovery of specific genes within large gene families and recovery of rare gene transcripts that are selectively expressed in recalcitrant tissues. To substantiate the validity of the approaches, we describe the combination of methods used in the gene and enzyme discovery of a cytochrome P450 in the fragrant heartwood of tropical sandalwood, responsible for the fragrance defining, final step in the biosynthesis of (Z)-santalols.

Terpene synthases from Cannabis sativa

Cannabis (Cannabis sativa) plants produce and accumulate a terpene-rich resin in glandular trichomes, which are abundant on the surface of the female inflorescence. Bouquets of different monoterpenes and sesquiterpenes are important components of cannabis resin as they define some of the unique organoleptic properties and may also influence medicinal qualities of different cannabis strains and varieties. Transcriptome analysis of trichomes of the cannabis hemp variety ‘Finola’ revealed sequences of all stages of terpene biosynthesis. Nine cannabis terpene synthases (CsTPS) were identified in subfamilies TPS-a and TPS-b. Functional characterization identified mono- and sesqui-TPS, whose products collectively comprise most of the terpenes of ‘Finola’ resin, including major compounds such as β-myrcene, (E)-β-ocimene, (-)-limonene, (+)-α-pinene, β-caryophyllene, and α-humulene. Transcripts associated with terpene biosynthesis are highly expressed in trichomes compared to non-resin producing tissues. Knowledge of the CsTPS gene family may offer opportunities for selection and improvement of terpene profiles of interest in different cannabis strains and varieties.

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.

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.

Polymer Coatings of Cochlear Implant Electrode Surface - An Option for Improving Electrode-Nerve-Interface by Blocking Fibroblast Overgrowth

Overgrowth of connective tissue and scar formation induced by the electrode array insertion increase the impedance and, thus, diminish the interactions between neural probes as like cochlear implants (CI) and the target tissue. Therefore, it is of great clinical interest to modify the carrier material of the electrodes to improve the electrode nerve interface for selective cell adhesion. On one side connective tissue growth needs to be reduced to avoid electrode array encapsulation, on the other side the carrier material should not compromise the interaction with neuronal cells. The present in vitro-study qualitatively and quantitatively characterises the interaction of fibroblasts, glial cells and spiral ganglion neurons (SGN) with ultrathin poly(N,N-dimethylacrylamide) (PDMAA), poly(2-ethyloxazoline) (PEtOx) and poly([2-methacryloyloxy)ethyl]trimethylammoniumchlorid) (PMTA) films immobilised onto glass surfaces using a photoreactive anchor layer. The layer thickness and hydrophilicity of the polymer films were characterised by ellipsometric and water contact angle measurement. Moreover the topography of the surfaces was investigated using atomic force microscopy (AFM). The neuronal and non-neuronal cells were dissociated from spiral ganglions of postnatal rats and cultivated for 48 h on top of the polymer coatings. Immunocytochemical staining of neuronal and intermediary filaments revealed that glial cells predominantly attached on PMTA films, but not on PDMAA and PEtOx monolayers. Hereby, strong survival rates and neurite outgrowth were only found on PMTA, whereas PDMAA and PEtOx coatings significantly reduced the SG neuron survival and neuritogenesis. As also shown by scanning electron microscopy (SEM) SGN strongly survived and retained their differentiated phenotype only on PMTA. In conclusion, survival and neuritogenesis of SGN may be associated with the extent of the glial cell growth. Since PMTA was the only of the polar polymers used in this study bearing a cationic charge, it can be assumed that this charge favours adhesion of both glial cells and SG neurons glial cells and SGN.

Gene expression analysis of overwintering mountain pine beetle larvae suggests multiple systems involved in overwintering stress, cold hardiness, and preparation for spring development

Cold-induced mortality has historically been a key aspect of mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae), population control, but little is known about the molecular basis for cold tolerance in this insect. We used RNA-seq analysis to monitor gene expression patterns of mountain pine beetle larvae at four time points during their overwintering period—early-autumn, late-autumn, early-spring, and late-spring. Changing transcript profiles over the winter indicates a multipronged physiological response from larvae that is broadly characterized by gene transcripts involved in insect immune responses and detoxification during the autumn. In the spring, although transcripts associated with developmental process are present, there was no particular biological process dominating the transcriptome.

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.

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.

Quantitative metabolome, proteome and transcriptome analysis of midgut and fat body tissues in the mountain pine beetle, Dendroctonus ponderosae Hopkins, and insights into pheromone biosynthesis

Bark beetles (Coleoptera: Scolytinae) are pests of many forests around the world. The mountain pine beetle (MPB), Dendroctonus ponderosae Hopkins, is a significant pest of western North American pine forests. The MPB is able to overcome the defences of pine trees through pheromone-assisted aggregation that results in a mass attack of host trees. These pheromones, both male and female produced, are believed to be biosynthesized in the midgut and/or fat bodies of these insects. We used metabolite analysis, quantitative proteomics (iTRAQ) and transcriptomics (RNA-seq) to identify proteins and transcripts differentially expressed between sexes and between tissues when treated with juvenile hormone III. Juvenile hormone III induced frontalin biosynthesis in males and trans-verbenol biosynthesis in females, as well as affected the expression of many proteins and transcripts in sex- and tissue-specific ways. Based on these analyses, we identified candidate genes involved in the biosynthesis of frontalin, exo-brevicomin, and trans-verbenol pheromones.

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.

Enzymes for synthetic biology of ambroxide-related diterpenoid fragrance compounds

Ambrox and related ambroxides are highly priced in the fragrance industry, and valued for their delicate odor and fixative properties. Historically, ambrox was obtained from ambergris, a waxy excretion produced by sperm whales, now an endangered species. Synthetic ambroxides have replaced ambergris in perfume manufacture. Plant labdane diterpenoids can serve as starting material for ambroxide synthesis. Among these, the diterpene alcohol sclareol is the major industrial precursor obtained from cultivated clary sage (Salvia sclarea). In plants, a large family of diterpene synthase (diTPS) enzymes controls key reactions in diterpenoid biosynthesis. Advanced metabolite profiling and high-throughput sequencing of fragrant and medicinal plants have accelerated discovery of novel diTPS functions, providing a resource for combinatorial synthetic biology and metabolic engineering approaches. This chapter highlights recent progress on the discovery, characterization, and engineering of plant diTPSs with potential uses in ambroxide production. It features biosynthesis of sclareol, cis-abienol, and diterpene resin acids, as sources of genes and enzymes for diterpenoid bioproducts.

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.

Insights into conifer giga-genomes

Insights from sequenced genomes of major land plant lineages have advanced research in almost every aspect of plant biology. Until recently, however, assembled genome sequences of gymnosperms have been missing from this picture. Conifers of the pine family (Pinaceae) are a group of gymnosperms that dominate large parts of the world's forests. Despite their ecological and economic importance, conifers seemed long out of reach for complete genome sequencing, due in part to their enormous genome size (20-30 Gb) and the highly repetitive nature of their genomes. Technological advances in genome sequencing and assembly enabled the recent publication of three conifer genomes: white spruce (Picea glauca), Norway spruce (Picea abies), and loblolly pine (Pinus taeda). These genome sequences revealed distinctive features compared with other plant genomes and may represent a window into the past of seed plant genomes. This Update highlights recent advances, remaining challenges, and opportunities in light of the publication of the first conifer and gymnosperm genomes.

Proteomics indicators of the rapidly shifting physiology from whole mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae), adults during early host colonization

We developed proteome profiles for host colonizing mountain pine beetle adults, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). Adult insects were fed in pairs on fresh host lodgepole pine, Pinus contorta Dougl. ex Loud, phloem tissue. The proteomes of fed individuals were monitored using iTRAQ and compared to those of starved beetles, revealing 757 and 739 expressed proteins in females and males, respectively, for which quantitative information was obtained. Overall functional category distributions were similar for males and females, with the majority of proteins falling under carbohydrate metabolism (glycolysis, gluconeogenesis, citric acid cycle), structure (cuticle, muscle, cytoskeleton), and protein and amino acid metabolism. Females had 23 proteins with levels that changed significantly with feeding (p<0.05, FDR<0.20), including chaperones and enzymes required for vitellogenesis. In males, levels of 29 proteins changed significantly with feeding (p<0.05, FDR<0.20), including chaperones as well as motor proteins. Only two proteins, both chaperones, exhibited a significant change in both females and males with feeding. Proteins with differential accumulation patterns in females exhibited higher fold changes with feeding than did those in males. This difference may be due to major and rapid physiological changes occurring in females upon finding a host tree during the physiological shift from dispersal to reproduction. The significant accumulation of chaperone proteins, a cytochrome P450, and a glutathione S-transferase, indicate secondary metabolite-induced stress physiology related to chemical detoxification during early host colonization. The females' activation of vitellogenin only after encountering a host indicates deliberate partitioning of resources and a balancing of the needs of dispersal and reproduction.

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.

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.

Cloning and characterization of chitinases from interior spruce and lodgepole pine

Chitinases have been implicated in the defence of conifers against insects and pathogens. cDNA for six chitinases were cloned from interior spruce (Picea glauca x engelmannii) and four from lodgepole pine (Pinus contorta). The cloned interior spruce chitinases were annotated class I PgeChia1-1 and PgeChia1-2, class II PgeChia2-1, class IV PgeChia4-1, and class VII PgeChia7-1 and PgeChia7-2; lodgepole pine chitinases were annotated class I PcChia1-1, class IV PcChia4-1, and class VII PcChia7-1 and PcChia7-2. Chitinases were expressed in Escherichia coli with maltose-binding-protein tags and soluble proteins purified. Functional characterization demonstrated chitinolytic activity for the three class I chitinases PgeChia1-1, PgeChia1-2 and PcChia1-1. Transcript analysis established strong induction of most of the tested chitinases, including all three class I chitinases, in interior spruce and lodgepole pine in response to inoculation with bark beetle associated fungi (Leptographium abietinum and Grosmannia clavigera) and in interior spruce in response to weevil (Pissodes strobi) feeding. Evidence of chitinolytic activity and inducibility by fungal and insect attack support the involvement of these chitinases in conifer defense.

Manoyl oxide (13R), the biosynthetic precursor of forskolin, is synthesized in specialized root cork cells in Coleus forskohlii

Forskolin, a complex labdane diterpenoid found in the root of Coleus forskohlii (Lamiaceae), has received attention for its broad range of pharmacological activities, yet the biosynthesis has not been elucidated. We detected forskolin in the root cork of C. forskohlii in a specialized cell type containing characteristic structures with histochemical properties consistent with oil bodies. Organelle purification and chemical analysis confirmed the localization of forskolin and of its simplest diterpene precursor backbone, (13R) manoyl oxide, to the oil bodies. The labdane diterpene backbone is typically synthesized by two successive reactions catalyzed by two distinct classes of diterpene synthases. We have recently described the identification of a small gene family of diterpene synthase candidates (CfTPSs) in C. forskohlii. Here, we report the functional characterization of four CfTPSs using in vitro and in planta assays. CfTPS2, which synthesizes the intermediate copal-8-ol diphosphate, in combination with CfTPS3 resulted in the stereospecific formation of (13R) manoyl oxide, while the combination of CfTPS1 and CfTPS3 or CfTPS4 led to formation of miltiradiene, precursor of abietane diterpenoids in C. forskohlii. Expression profiling and phylogenetic analysis of the CfTPS family further support the functional diversification and distinct roles of the individual diterpene synthases and the involvement of CfTPS1 to CfTPS4 in specialized metabolism and of CfTPS14 and CfTPS15 in general metabolism. Our findings pave the way toward the discovery of the remaining components of the pathway to forskolin, likely localized in this specialized cell type, and support a role of oil bodies as storage organelles for lipophilic bioactive metabolites.

CYP345E2, an antenna-specific cytochrome P450 from the mountain pine beetle, Dendroctonus ponderosae Hopkins, catalyses the oxidation of pine host monoterpene volatiles

The mountain pine beetle (MPB, Dendroctonus ponderosae Hopkins) is a significant pest of western North American pine forests. This beetle responds to pheromones and host volatiles in order to mass attack and thus overcome the terpenoid chemical defences of its host. The ability of MPB antennae to rapidly process odorants is necessary to avoid odorant receptor saturation and thus the enzymes responsible for odorant clearance are an important aspect of host colonization. An antenna-specific cytochrome P450, DponCYP345E2, is the most highly expressed transcript in adult MPB antenna. In in vitro assays with recombinant enzyme, DponCYP345E2 used several pine host monoterpenes as substrates, including (+)-(3)-carene, (+)-β-pinene, (-)-β-pinene, (+)-limonene, (-)-limonene, (-)-camphene, (+)-α-pinene, (-)-α-pinene, and terpinolene. The substrates were epoxidized or hydroxylated, depending upon the substrate. To complement DponCYP345E2, we also functionally characterized the NADPH-dependent cytochrome P450 reductase and the cytochrome b5 from MPB. DponCYP345E2 is the first cytochrome P450 to be functionally characterized in insect olfaction and in MPB.

Disentangling detoxification: gene expression analysis of feeding mountain pine beetle illuminates molecular-level host chemical defense detoxification mechanisms

The mountain pine beetle, Dendroctonus ponderosae, is a native species of bark beetle (Coleoptera: Curculionidae) that caused unprecedented damage to the pine forests of British Columbia and other parts of western North America and is currently expanding its range into the boreal forests of central and eastern Canada and the USA. We conducted a large-scale gene expression analysis (RNA-seq) of mountain pine beetle male and female adults either starved or fed in male-female pairs for 24 hours on lodgepole pine host tree tissues. Our aim was to uncover transcripts involved in coniferophagous mountain pine beetle detoxification systems during early host colonization. Transcripts of members from several gene families significantly increased in insects fed on host tissue including: cytochromes P450, glucosyl transferases and glutathione S-transferases, esterases, and one ABC transporter. Other significantly increasing transcripts with potential roles in detoxification of host defenses included alcohol dehydrogenases and a group of unexpected transcripts whose products may play an, as yet, undiscovered role in host colonization by mountain pine beetle.

Biosynthesis of Sandalwood Oil: Santalum album CYP76F cytochromes P450 produce santalols and bergamotol

Sandalwood oil is one of the world’s most highly prized essential oils, appearing in many high-end perfumes and fragrances. Extracted from the mature heartwood of several Santalum species, sandalwood oil is comprised mainly of sesquiterpene olefins and alcohols. Four sesquiterpenols, α-, β-, and epi-β-santalol and α-exo-bergamotol, make up approximately 90% of the oil of Santalum album. These compounds are the hydroxylated analogues of α-, β-, and epi-β-santalene and α-exo-bergamotene. By mining a transcriptome database of S. album for candidate cytochrome P450 genes, we cloned and characterized cDNAs encoding a small family of ten cytochrome P450-dependent monooxygenases annotated as SaCYP76F37v1, SaCYP76F37v2, SaCYP76F38v1, SaCYP76F38v2, SaCYP76F39v1, SaCYP76F39v2, SaCYP76F40, SaCYP76F41, SaCYP76F42, and SaCYP76F43. Nine of these genes were functionally characterized using in vitro assays and yeast in vivo assays to encode santalene/bergamotene oxidases and bergamotene oxidases. These results provide a foundation for production of sandalwood oil for the fragrance industry by means of metabolic engineering, as demonstrated with proof-of-concept formation of santalols and bergamotol in engineered yeast cells, simultaneously addressing conservation challenges by reducing pressure on supply of sandalwood from native forests.

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.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data

White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20,356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies.

AVAILABILITY

The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

The genome and transcriptome of the pine saprophyte Ophiostoma piceae, and a comparison with the bark beetle-associated pine pathogen Grosmannia clavigera

BACKGROUND

Ophiostoma piceae is a wood-staining fungus that grows in the sapwood of conifer logs and lumber. We sequenced its genome and analyzed its transcriptomes under a range of growth conditions. A comparison with the genome and transcriptomes of the mountain pine beetle-associated pathogen Grosmannia clavigera highlights differences between a pathogen that colonizes and kills living pine trees and a saprophyte that colonizes wood and the inner bark of dead trees.

RESULTS

We assembled a 33 Mbp genome in 45 scaffolds, and predicted approximately 8,884 genes. The genome size and gene content were similar to those of other ascomycetes. Despite having similar ecological niches, O. piceae and G. clavigera showed no large-scale synteny. We identified O. piceae genes involved in the biosynthesis of melanin, which causes wood discoloration and reduces the commercial value of wood products. We also identified genes and pathways involved in growth on simple carbon sources and in sapwood, O. piceae's natural substrate. Like the pathogen, the saprophyte is able to tolerate terpenes, which are a major class of pine tree defense compounds; unlike the pathogen, it cannot utilize monoterpenes as a carbon source.

CONCLUSIONS

This work makes available the second annotated genome of a softwood ophiostomatoid fungus, and suggests that O. piceae's tolerance to terpenes may be due in part to these chemicals being removed from the cells by an ABC transporter that is highly induced by terpenes. The data generated will provide the research community with resources for work on host-vector-fungus interactions for wood-inhabiting, beetle-associated saprophytes and pathogens.

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.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data

White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20 356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies.

Availability: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

Contact:ibirol@bcgsc.ca

Supplementary information:Supplementary data are available at Bioinformatics online.

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.

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.