Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait.

BACKGROUND: Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution. RESULTS: Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits. CONCLUSIONS: Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

Monoterpenyl esters in juvenile mountain pine beetle and sex-specific release of the aggregation pheromone trans-verbenol.

A recent outbreak of mountain pine beetle (MPB) has spread over more than 25 million hectares of pine forests in western North America, affecting pine species of sensitive boreal and mountain ecosystems. During initial host colonization, female MPB produce and release the aggregation pheromone trans-verbenol to coordinate a mass attack of individual trees. trans-Verbenol is formed by hydroxylation of α-pinene, a monoterpene of the pine oleoresin defense. It is thought that adult females produce and immediately release trans-verbenol when encountering α-pinene on a new host tree. Here, we show that both sexes of MPB accumulate the monoterpenyl esters verbenyl oleate and verbenyl palmitate during their development in the brood tree. Verbenyl oleate and verbenyl palmitate were retained in adult female MPB until the time of emergence from brood trees, but were depleted in males. Adult females released trans-verbenol in response to treatment with juvenile hormone III (JHIII). While both sexes produced verbenyl esters when exposed to α-pinene, only females responded to JHIII with release of trans-verbenol. Accumulation of verbenyl esters at earlier life stages may allow adult females to release the aggregation pheromone trans-verbenol upon landing on a new host tree, independent of access to α-pinene. Formation of verbenyl esters may be part of a general detoxification system to overcome host monoterpene defenses in both sexes, from which a specialized and female-specific system of pheromone biosynthesis and release may have evolved.

Genomic content of chemosensory genes correlates with host range in wood-boring beetles (Dendroctonus ponderosae, Agrilus planipennis, and Anoplophora glabripennis).

BACKGROUND: Olfaction and gustation underlie behaviors that are crucial for insect fitness, such as host and mate selection. The detection of semiochemicals is mediated via proteins from large and rapidly evolving chemosensory gene families; however, the links between a species' ecology and the diversification of these genes remain poorly understood. Hence, we annotated the chemosensory genes from genomes of select wood-boring coleopterans, and compared the gene repertoires from stenophagous species with those from polyphagous species. RESULTS: We annotated 86 odorant receptors (ORs), 60 gustatory receptors (GRs), 57 ionotropic receptors (IRs), 4 sensory neuron membrane proteins (SNMPs), 36 odorant binding proteins (OBPs), and 11 chemosensory proteins (CSPs) in the mountain pine beetle (Dendroctonus ponderosae), and 47 ORs, 30 GRs, 31 IRs, 4 SNMPs, 12 OBPs, and 14 CSPs in the emerald ash borer (Agrilus planipennis). Four SNMPs and 17 CSPs were annotated in the polyphagous wood-borer Anoplophora glabripennis. The gene repertoires in the stenophagous D. ponderosae and A. planipennis are reduced compared with those in the polyphagous A. glabripennis and T. castaneum, which is largely manifested through small gene lineage expansions and entire lineage losses. Alternative splicing of GR genes was limited in D. ponderosae and apparently absent in A. planipennis, which also seems to have lost one carbon dioxide receptor (GR1). A. planipennis has two SNMPs, which are related to SNMP3 in T. castaneum. D. ponderosae has two alternatively spliced OBP genes, a novel OBP "tetramer", and as many as eleven IR75 members. Simple orthology was generally rare in beetles; however, we found one clade with orthologues of putative bitter-taste GRs (named the "GR215 clade"), and conservation of IR60a from Drosophila melanogaster. CONCLUSIONS: Our genome annotations represent important quantitative and qualitative improvements of the original datasets derived from transcriptomes of D. ponderosae and A. planipennis, facilitating evolutionary analysis of chemosensory genes in the Coleoptera where only a few genomes were previously annotated. Our analysis suggests a correlation between chemosensory gene content and host specificity in beetles. Future studies should include additional species to consolidate this correlation, and functionally characterize identified proteins as an important step towards improved control of these pests.

Cytochromes P450 Preferentially Expressed in Antennae of the Mountain Pine Beetle.

The mountain pine beetle (MPB, Dendroctonus ponderosae) is a forest insect pest endemic to western North America. During dispersal and host colonization, MPB identify suitable host trees by olfaction of monoterpene volatiles, contend with host terpene defenses, and communicate with conspecifics using terpenoid and other pheromones. Cytochromes P450 (P450s) have been proposed to function in MPB olfaction, terpene detoxification, and pheromone biosynthesis. Here, we identified P450s that were abundant in the antennae transcriptome. Analysis of transcript levels across different life stages and tissues in females and males showed additional expression of most of these P450s in the midgut or fat body. These expression profiles suggest specific or overlapping functions in olfaction, detoxification, and pheromone biosynthesis.

Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism.

White spruce (Picea glauca), a gymnosperm tree, has been established as one of the models for conifer genomics. We describe the draft genome assemblies of two white spruce genotypes, PG29 and WS77111, innovative tools for the assembly of very large genomes, and the conifer genomics resources developed in this process. The two white spruce genotypes originate from distant geographic regions of western (PG29) and eastern (WS77111) North America, and represent elite trees in two Canadian tree-breeding programs. We present an update (V3 and V4) for a previously reported PG29 V2 draft genome assembly and introduce a second white spruce genome assembly for genotype WS77111. Assemblies of the PG29 and WS77111 genomes confirm the reconstructed white spruce genome size in the 20 Gbp range, and show broad synteny. Using the PG29 V3 assembly and additional white spruce genomics and transcriptomics resources, we performed MAKER-P annotation and meticulous expert annotation of very large gene families of conifer defense metabolism, the terpene synthases and cytochrome P450s. We also comprehensively annotated the white spruce mevalonate, methylerythritol phosphate and phenylpropanoid pathways. These analyses highlighted the large extent of gene and pseudogene duplications in a conifer genome, in particular for genes of secondary (i.e. specialized) metabolism, and the potential for gain and loss of function for defense and adaptation.

Frontalin pheromone biosynthesis in the mountain pine beetle, Dendroctonus ponderosae, and the role of isoprenyl diphosphate synthases.

The mountain pine beetle (Dendroctonus ponderosae Hopkins) is the most destructive pest of western North American pine forests. Adult males produce frontalin, an eight-carbon antiaggregation pheromone, via the mevalonate pathway, as part of several pheromones that initiate and modulate the mass attack of host trees. Frontalin acts as a pheromone, attractant, or kairomone in most Dendroctonus species, other insects, and even elephants. 6-Methylhept-6-en-2-one, a frontalin precursor, is hypothesized to originate from 10-carbon geranyl diphosphate (GPP), 15-carbon farnesyl diphosphate (FPP), or 20-carbon geranylgeranyl diphosphate (GGPP) via a dioxygenase- or cytochrome P450-mediated carbon-carbon bond cleavage. To investigate the role of isoprenyl diphosphate synthases in pheromone biosynthesis, we characterized a bifunctional GPP/FPP synthase and a GGPP synthase in the mountain pine beetle. The ratio of GPP to FPP produced by the GPP/FPP synthase was highly dependent on the ratio of the substrates isopentenyl diphosphate and dimethylallyl diphosphate used in the assay. Transcript levels in various tissues and life stages suggested that GGPP rather than GPP or FPP is used as a precursor to frontalin. Reduction of transcript levels by RNA interference of the isoprenyl diphosphate synthases identified GGPP synthase as having the largest effect on frontalin production, suggesting that frontalin is derived from a 20-carbon isoprenoid precursor rather than from the 10- or 15-carbon precursors.

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.

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.

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

UNLABELLED: 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.

Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae).

BACKGROUND: The European spruce bark beetle, Ips typographus, and the North American mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae), are severe pests of coniferous forests. Both bark beetle species utilize aggregation pheromones to coordinate mass-attacks on host trees, while odorants from host and non-host trees modulate the pheromone response. Thus, the bark beetle olfactory sense is of utmost importance for fitness. However, information on the genes underlying olfactory detection has been lacking in bark beetles and is limited in Coleoptera. We assembled antennal transcriptomes from next-generation sequencing of I. typographus and D. ponderosae to identify members of the major chemosensory multi-gene families. RESULTS: Gene ontology (GO) annotation indicated that the relative abundance of transcripts associated with specific GO terms was highly similar in the two species. Transcripts with terms related to olfactory function were found in both species. Focusing on the chemosensory gene families, we identified 15 putative odorant binding proteins (OBP), 6 chemosensory proteins (CSP), 3 sensory neuron membrane proteins (SNMP), 43 odorant receptors (OR), 6 gustatory receptors (GR), and 7 ionotropic receptors (IR) in I. typographus; and 31 putative OBPs, 11 CSPs, 3 SNMPs, 49 ORs, 2 GRs, and 15 IRs in D. ponderosae. Predicted protein sequences were compared with counterparts in the flour beetle, Tribolium castaneum, the cerambycid beetle, Megacyllene caryae, and the fruit fly, Drosophila melanogaster. The most notable result was found among the ORs, for which large bark beetle-specific expansions were found. However, some clades contained receptors from all four beetle species, indicating a degree of conservation among some coleopteran OR lineages. Putative GRs for carbon dioxide and orthologues for the conserved antennal IRs were included in the identified receptor sets. CONCLUSIONS: The protein families important for chemoreception have now been identified in three coleopteran species (four species for the ORs). Thus, this study allows for improved evolutionary analyses of coleopteran olfaction. Identification of these proteins in two of the most destructive forest pests, sharing many semiochemicals, is especially important as they might represent novel targets for population control.

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.

The primary diterpene synthase products of Picea abies levopimaradiene/abietadiene synthase (PaLAS) are epimers of a thermally unstable diterpenol.

The levopimaradiene/abietadiene synthase from Norway spruce (Picea abies; PaLAS) has previously been reported to produce a mixture of four diterpene hydrocarbons when incubated with geranylgeranyl diphosphate as the substrate: levopimaradiene, abietadiene, neoabietadiene, and palustradiene. However, variability in the assay products observed by GC-MS of this and orthologous conifer diterpene synthases over the past 15 years suggested that these diterpenes may not be the initial enzyme assay products but are rather the products of dehydration of an unstable alcohol. We have identified epimers of the thermally unstable allylic tertiary alcohol 13-hydroxy-8(14)-abietene as the products of PaLAS. The identification of these compounds, not previously described in conifers, as the initial products of PaLAS has considerable implications for our understanding of the complexity of the biosynthetic pathway of the structurally diverse diterpene resin acids of conifer defense.

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

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

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.

Chromosome-level genome assembly reveals genomic architecture of northern range expansion in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae).

Genome sequencing methods and assembly tools have improved dramatically since the 2013 publication of draft genome assemblies for the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). We conducted proximity ligation library sequencing and scaffolding to improve contiguity, and then used linkage mapping and recent bioinformatic tools for correction and further improvement. The new assemblies have dramatically improved contiguity and gaps compared to the originals: N50 values increased 26- to 36-fold, and the number of gaps were reduced by half. Ninety per cent of the content of the assemblies is now contained in 12 and 11 scaffolds for the female and male assemblies, respectively. Based on linkage mapping information, the 12 largest scaffolds in both assemblies represent all 11 autosomal chromosomes and the neo-X chromosome. These assemblies now have nearly chromosome-sized scaffolds and will be instrumental for studying genomic architecture, chromosome evolution, population genomics, functional genomics, and adaptation in this and other pest insects. We also identified regions in two chromosomes, including the ancestral-X portion of the neo-X chromosome, with elevated differentiation between northern and southern Canadian populations.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins.

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

Transcriptome mining, functional characterization, and phylogeny of a large terpene synthase gene family in spruce (Picea spp.).

BACKGROUND: In conifers, terpene synthases (TPSs) of the gymnosperm-specific TPS-d subfamily form a diverse array of mono-, sesqui-, and diterpenoid compounds, which are components of the oleoresin secretions and volatile emissions. These compounds contribute to defence against herbivores and pathogens and perhaps also protect against abiotic stress. RESULTS: The availability of extensive transcriptome resources in the form of expressed sequence tags (ESTs) and full-length cDNAs in several spruce (Picea) species allowed us to estimate that a conifer genome contains at least 69 unique and transcriptionally active TPS genes. This number is comparable to the number of TPSs found in any of the sequenced and well-annotated angiosperm genomes. We functionally characterized a total of 21 spruce TPSs: 12 from Sitka spruce (P. sitchensis), 5 from white spruce (P. glauca), and 4 from hybrid white spruce (P. glauca × P. engelmannii), which included 15 monoterpene synthases, 4 sesquiterpene synthases, and 2 diterpene synthases. CONCLUSIONS: The functional diversity of these characterized TPSs parallels the diversity of terpenoids found in the oleoresin and volatile emissions of Sitka spruce and provides a context for understanding this chemical diversity at the molecular and mechanistic levels. The comparative characterization of Sitka spruce and Norway spruce diterpene synthases revealed the natural occurrence of TPS sequence variants between closely related spruce species, confirming a previous prediction from site-directed mutagenesis and modelling.

Identification and functional characterization of monofunctional ent-copalyl diphosphate and ent-kaurene synthases in white spruce reveal different patterns for diterpene synthase evolution for primary and secondary metabolism in gymnosperms.

The biosynthesis of the tetracyclic diterpene ent-kaurene is a critical step in the general (primary) metabolism of gibberellin hormones. ent-Kaurene is formed by a two-step cyclization of geranylgeranyl diphosphate via the intermediate ent-copalyl diphosphate. In a lower land plant, the moss Physcomitrella patens, a single bifunctional diterpene synthase (diTPS) catalyzes both steps. In contrast, in angiosperms, the two consecutive cyclizations are catalyzed by two distinct monofunctional enzymes, ent-copalyl diphosphate synthase (CPS) and ent-kaurene synthase (KS). The enzyme, or enzymes, responsible for ent-kaurene biosynthesis in gymnosperms has been elusive. However, several bifunctional diTPS of specialized (secondary) metabolism have previously been characterized in gymnosperms, and all known diTPSs for resin acid biosynthesis in conifers are bifunctional. To further understand the evolution of ent-kaurene biosynthesis as well as the evolution of general and specialized diterpenoid metabolisms in gymnosperms, we set out to determine whether conifers use a single bifunctional diTPS or two monofunctional diTPSs in the ent-kaurene pathway. Using a combination of expressed sequence tag, full-length cDNA, genomic DNA, and targeted bacterial artificial chromosome sequencing, we identified two candidate CPS and KS genes from white spruce (Picea glauca) and their orthologs in Sitka spruce (Picea sitchensis). Functional characterization of the recombinant enzymes established that ent-kaurene biosynthesis in white spruce is catalyzed by two monofunctional diTPSs, PgCPS and PgKS. Comparative analysis of gene structures and enzyme functions highlights the molecular evolution of these diTPSs as conserved between gymnosperms and angiosperms. In contrast, diTPSs for specialized metabolism have evolved differently in angiosperms and gymnosperms.

Functional plasticity of paralogous diterpene synthases involved in conifer defense.

The diversity of terpenoid compounds produced by plants plays an important role in mediating various plant-herbivore, plant-pollinator, and plant-pathogen interactions. This diversity has resulted from gene duplication and neofunctionalization of the enzymes that synthesize and subsequently modify terpenes. Two diterpene synthases in Norway spruce (Picea abies), isopimaradiene synthase and levopimaradiene/abietadiene synthase, provide the hydrocarbon precursors for most of the diterpene resin acids found in the defensive oleoresin of conifers. Although these paralogous enzymes are 91% identical at the amino acid level, one is a single-product enzyme, whereas the other is a multiproduct enzyme that forms completely different products. We used a rational approach of homology modeling, protein sequence comparison, domain swapping, and a series of reciprocal site-directed mutagenesis to identify the specific residues that direct the different product outcomes. A one-amino acid mutation switched the levopimaradiene/abietadiene synthase into producing isopimaradiene and sandaracopimaradiene and none of its normal products. Four mutations were sufficient to reciprocally reverse the product profiles for both of these paralogous enzymes while maintaining catalytic efficiencies similar to the wild-type enzymes. This study illustrates how neofunctionalization can result from relatively minor changes in protein sequence, increasing the diversity of secondary metabolites important for conifer defense.

Cytochromes P450: terpene detoxification and pheromone production in bark beetles.

Bark beetles (family: Curculionidae; subfamily: Scolytinae) in the Dendroctonus and Ips genera are the most destructive forest pests in the Northern hemisphere. They use cytochromes P450 (P450s) to detoxify tree-produced terpenes to produce pheromones, in de novo pheromone production and to oxidize odorants on antennae. Many Dendroctonus spp. use trans-verbenol as an aggregation pheromone, and it is formed from host-tree produced α-pinene hydroxylated by CYP6DE1 during larval stages, stored as verbenyl ester of fatty acids, and then released when the female begins feeding on a new host tree. Ips spp. hydroxylate de novo produced myrcene to form ipsdienol. Subsequent steps form the appropriate enantiomeric composition of ipsdienol and convert ipsdienol to ipsenol. In this article we review recent progress in elucidating the functions of P450s in Ips and Dendroctonus species and in doing so provide insights into the role of these enzymes in host phytochemical detoxification and pheromone production.