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.

Genetic architecture of terpene chemistry and growth traits and the impact of inbreeding on these traits in western redcedar (Thuja plicata).

Western redcedar (WRC; Thuja plicata) is a conifer of the Pacific Northwest of North America prized for its durable and rot-resistant wood. WRC has naturally low outcrossing rates and readily self-fertilizes in nature. Challenges faced in WRC breeding and propagation involve selecting trees for accelerated growth while also ensuring enhanced heartwood rot resistance and resistance to ungulate browsing, as well as mitigating potential effects of inbreeding depression. Terpenes, a large and diverse class of specialized metabolites, confer both rot and browse resistance in the wood and foliage of WRC, respectively. Using a Bayesian modelling approach, we isolated single nucleotide polymorphism (SNP) markers estimated to be associated with three different foliar terpene traits and four different heartwood terpene traits, as well as two growth traits. We found that all traits were complex, being associated with between 1700 and 3600 SNPs linked with putatively causal loci, with significant polygenic components. Growth traits tended to have a larger polygenic component while terpene traits had larger major gene components; SNPs with small or polygenic effect were spread across the genome, while larger-effect SNPs tended to be localized to specific linkage groups. To determine whether there was inbreeding depression for terpene chemistry or growth traits, we used mixed linear models for a genomic selection training population to estimate the effect of the inbreeding coefficient F on foliar terpenes, heartwood terpenes and several growth and dendrochronological traits. We did not find significant inbreeding depression for any assessed trait. We further assessed inbreeding depression across four generations of complete selfing and found that not only was inbreeding depression not significant but that selection for height growth was the only significant predictor for growth during selfing, suggesting that inbreeding depression due to selfing during operational breeding can be mitigated by increased selection intensity.

Genomic selection reveals hidden relatedness and increased breeding efficiency in western redcedar polycross breeding.

Western redcedar (WRC) is an ecologically and economically important forest tree species characterized by low genetic diversity with high self-compatibility and high heartwood durability. Using sequence capture genotyping of target genic and non-genic regions, we genotyped 44 parent trees and 1520 offspring trees representing 26 polycross (PX) families collected from three progeny test sites using 45,378 SNPs. Trees were phenotyped for eight traits related to growth, heartwood and foliar chemistry associated with wood durability and deer browse resistance. We used the genomic realized relationship matrix for paternity assignment, maternal pedigree correction, and to estimate genetic parameters. We compared genomics-based (GBLUP) and two pedigree-based (ABLUP: polycross and reconstructed full-sib [FS] pedigrees) models. Models were extended to estimate dominance genetic effects. Pedigree reconstruction revealed significant unequal male contribution and separated the 26 PX families into 438 FS families. Traditional maternal PX pedigree analysis resulted in up to 51% overestimation in genetic gain and 44% in diversity. Genomic analysis resulted in up to 22% improvement in offspring breeding value (BV) theoretical accuracy, 35% increase in expected genetic gain for forward selection, and doubled selection intensity for backward selection. Overall, all traits showed low to moderate heritability (0.09-0.28), moderate genotype by environment interaction (type-B genetic correlation: 0.51-0.80), low to high expected genetic gain (6.01%-55%), and no significant negative genetic correlation reflecting no large trade-offs for multi-trait selection. Only three traits showed a significant dominance effect. GBLUP resulted in smaller but more accurate heritability estimates for five traits, but larger estimates for the wood traits. Comparison between all, genic-coding, genic-non-coding and intergenic SNPs showed little difference in genetic estimates. In summary, we show that GBLUP overcomes the PX limitations, successfully captures expected historical and hidden relatedness as well as linkage disequilibrium (LD), and results in increased breeding efficiency in WRC.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

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.

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.

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.

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.

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.

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.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains.

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below -40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

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.

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.

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.

Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest.

BACKGROUND: The mountain pine beetle, Dendroctonus ponderosae Hopkins, is the most serious insect pest of western North American pine forests. A recent outbreak destroyed more than 15 million hectares of pine forests, with major environmental effects on forest health, and economic effects on the forest industry. The outbreak has in part been driven by climate change, and will contribute to increased carbon emissions through decaying forests. RESULTS: We developed a genome sequence resource for the mountain pine beetle to better understand the unique aspects of this insect's biology. A draft de novo genome sequence was assembled from paired-end, short-read sequences from an individual field-collected male pupa, and scaffolded using mate-paired, short-read genomic sequences from pooled field-collected pupae, paired-end short-insert whole-transcriptome shotgun sequencing reads of mRNA from adult beetle tissues, and paired-end Sanger EST sequences from various life stages. We describe the cytochrome P450, glutathione S-transferase, and plant cell wall-degrading enzyme gene families important to the survival of the mountain pine beetle in its harsh and nutrient-poor host environment, and examine genome-wide single-nucleotide polymorphism variation. A horizontally transferred bacterial sucrose-6-phosphate hydrolase was evident in the genome, and its tissue-specific transcription suggests a functional role for this beetle. CONCLUSIONS: Despite Coleoptera being the largest insect order with over 400,000 described species, including many agricultural and forest pest species, this is only the second genome sequence reported in Coleoptera, and will provide an important resource for the Curculionoidea and other insects.