Sharon Jancsik

Genomics of hybrid poplar (Populus trichocarpa× deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences in poplar

As part of a genomics strategy to characterize inducible defences against insect herbivory in poplar, we developed a comprehensive suite of functional genomics resources including cDNA libraries, expressed sequence tags (ESTs) and a cDNA microarray platform. These resources are designed to complement the existing poplar genome sequence and poplar (Populus spp.) ESTs by focusing on herbivore‐ and elicitor‐treated tissues and incorporating normalization methods to capture rare transcripts. From a set of 15 standard, normalized or full‐length cDNA libraries, we generated 139 007 3′‐ or 5′‐end sequenced ESTs, representing more than one‐third of the c. 385 000 publicly available Populus ESTs. Clustering and assembly of 107 519 3′‐end ESTs resulted in 14 451 contigs and 20 560 singletons, altogether representing 35 011 putative unique transcripts, or potentially more than three‐quarters of the predicted c. 45 000 genes in the poplar genome. Using this EST resource, we developed a cDNA microarray containing 15 496 unique genes, which was utilized to monitor gene expression in poplar leaves in response to herbivory by forest tent caterpillars (Malacosoma disstria). After 24 h of feeding, 1191 genes were classified as up‐regulated, compared to only 537 down‐regulated. Functional classification of this induced gene set revealed genes with roles in plant defence (e.g. endochitinases, Kunitz protease inhibitors), octadecanoid and ethylene signalling (e.g. lipoxygenase, allene oxide synthase, 1‐aminocyclopropane‐1‐carboxylate oxidase), transport (e.g. ABC proteins, calreticulin), secondary metabolism [e.g. polyphenol oxidase, isoflavone reductase, (–)‐germacrene D synthase] and transcriptional regulation [e.g. leucine‐rich repeat transmembrane kinase, several transcription factor classes (zinc finger C3H type, AP2/EREBP, WRKY, bHLH)]. This study provides the first genome‐scale approach to characterize insect‐induced defences in a woody perennial providing a solid platform for functional investigation of plant–insect interactions in poplar.

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

BackgroundIn 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.ResultsThe 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.ConclusionsThe 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.

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.

Aminocyclopropane Carboxylic Acid Synthase Is a Regulated Step in Ethylene-Dependent Induced Conifer Defense. Full-Length cDNA Cloning of a Multigene Family, Differential Constitutive, and Wound- and Insect-Induced Expression, and Cellular and Subcellular Localization in Spruce and Douglas Fir

In conifer stems, formation of chemical defenses against insects or pathogens involves specialized anatomical structures of the phloem and xylem. Oleoresin terpenoids are formed in resin duct epithelial cells and phenolics accumulate in polyphenolic parenchyma cells. Ethylene signaling has been implicated in the induction of these chemical defenses. Recently, we reported the cloning of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) from spruce (Picea spp.) and Douglas fir (Pseudotsuga menziesii). ACO protein was constitutively expressed in Douglas fir and only weakly induced upon wounding. We now cloned seven full-length and one near full-length cDNA representing four distinct 1-aminocyclopropane-1-carboxylic acid synthases (ACS; ACS1, ACS2, ACS3, and ACS4) from spruce and Douglas fir. Cloning of ACS has not previously been reported for any gymnosperm. Using gene-specific, quantitative real-time polymerase chain reaction, we measured constitutive expression for the four ACS genes and the single-copy ACO gene in various tissues of Sitka spruce (Picea sitchensis) and in white spruce (Picea glauca) somatic embryos. ACO and ACS4 were ubiquitously expressed at high levels; ACS1 was predominantly expressed in developing embryos and ACS2 and ACS3 were expressed only at very low levels. Insect attack or mechanical wounding caused strong induction of ACS2 and ACS3 in Sitka spruce bark, a moderate increase in ACO transcripts, but had no effect on ACS1 and ACS4. ACS protein was also strongly induced following mechanical wounding in Douglas fir and was highly abundant in resin duct epithelial cells and polyphenolic parenchyma cells. These results suggest that ACS, but not ACO, is a regulated step in ethylene-induced conifer defense.

Evolution of Conifer Diterpene Synthases: Diterpene Resin Acid Biosynthesis in Lodgepole Pine and Jack Pine Involves Monofunctional and Bifunctional Diterpene Synthases

Summary: The biosynthesis of diterpenoid specialized (i.e. secondary) metabolites in conifers is emerging as a modular metabolic network that is built from multiple diterpene synthases and cytochromes P450 of variable functions and activities to produce a diverse and dynamic array of metabolites for plant defense. Diterpene resin acids (DRAs) are major components of pine (Pinus spp.) oleoresin. They play critical roles in conifer defense against insects and pathogens and as a renewable resource for industrial bioproducts. The core structures of DRAs are formed in secondary (i.e. specialized) metabolism via cycloisomerization of geranylgeranyl diphosphate (GGPP) by diterpene synthases (diTPSs). Previously described gymnosperm diTPSs of DRA biosynthesis are bifunctional enzymes that catalyze the initial bicyclization of GGPP followed by rearrangement of a (+)-copalyl diphosphate intermediate at two discrete class II and class I active sites. In contrast, similar diterpenes of gibberellin primary (i.e. general) metabolism are produced by the consecutive activity of two monofunctional class II and class I diTPSs. Using high-throughput transcriptome sequencing, we discovered 11 diTPS from jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta). Three of these were orthologous to known conifer bifunctional levopimaradiene/abietadiene synthases. Surprisingly, two sets of orthologous PbdiTPSs and PcdiTPSs were monofunctional class I enzymes that lacked functional class II active sites and converted (+)-copalyl diphosphate, but not GGPP, into isopimaradiene and pimaradiene as major products. Diterpene profiles and transcriptome sequences of lodgepole pine and jack pine are consistent with roles for these diTPSs in DRA biosynthesis. The monofunctional class I diTPSs of DRA biosynthesis form a new clade within the gymnosperm-specific TPS-d3 subfamily that evolved from bifunctional diTPS rather than monofunctional enzymes (TPS-c and TPS-e) of gibberellin metabolism. Homology modeling suggested alterations in the class I active site that may have contributed to their functional specialization relative to other conifer diTPSs.

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

Abstract 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.

Poplar defense against insects: genome analysis, full-length cDNA cloning, and transcriptome and protein analysis of the poplar Kunitz-type protease inhibitor family.

*Kunitz protease inhibitors (KPIs) feature prominently in poplar defense responses against insects. The increasing availability of genomics resources enabled a comprehensive analysis of the poplar (p)KPI family. *Using genome analysis, expressed sequence tag (EST) mining and full-length (FL)cDNA cloning we established an inventory and phylogeny of pKPIs. Microarray and real-time PCR analyses were used to profile pKPI gene expression following real or simulated insect attack. Proteomics of insect midgut content was used to monitor stability of pKPI protein. *We identified 31 pKPIs in the genome and validated gene models by EST mining and cloning of 41 unique FLcDNAs. Genome organization of the pKPI family, with six poplar-specific subfamilies, suggests that tandem duplications have played a major role in its expansion. pKPIs are expressed throughout the plant and many are strongly induced by insect attack, although insect-specific signals seem initially to suppress the tree pKPI response. We found substantial peptide coverage for a potentially intact pKPI protein in insect midgut after eating poplar leaves. *These results highlight the complexity of an important defense gene family in poplar with regard to gene family size, differential constitutive and insect-induced gene expression, and resilience of at least one pKPI protein to digestion by herbivores.

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)

BackgroundThe 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.ResultsWe 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.ConclusionIn 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.

A conifer UDP-sugar dependent glycosyltransferase contributes to acetophenone metabolism and defense against insects

A conifer glycosyltransferase is involved in the chemical resistance against spruce budworm, a major forest defoliator. Acetophenones are phenolic compounds involved in the resistance of white spruce (Picea glauca) against spruce budworm (Choristoneura fumiferiana), a major forest pest in North America. The acetophenones pungenol and piceol commonly accumulate in spruce foliage in the form of the corresponding glycosides, pungenin and picein. These glycosides appear to be inactive against the insect but can be cleaved by a spruce β-glucosidase, PgβGLU-1, which releases the active aglycons. The reverse glycosylation reaction was hypothesized to involve a family 1 UDP-sugar dependent glycosyltransferase (UGT) to facilitate acetophenone accumulation in the plant. Metabolite and transcriptome profiling over a developmental time course of white spruce bud burst and shoot growth revealed two UGTs, PgUGT5 and PgUGT5b, that glycosylate pungenol. Recombinant PgUGT5b enzyme produced mostly pungenin, while PgUGT5 produced mostly isopungenin. Both UGTs also were active in vitro on select flavonoids. However, the context of transcript and metabolite accumulation did not support a biological role in flavonoid metabolism but correlated with the formation of pungenin in growing shoots. Transcript levels of PgUGT5b were higher than those of PgUGT5 in needles across different genotypes of white spruce. These results support a role of PgUGT5b in the biosynthesis of the glycosylated acetophenone pungenin in white spruce.