Joerg Bohlmann

The Genome of Black Cottonwood, Populus trichocarpa (Torr. & Gray)

We report the draft genome of the black cottonwood tree, Populus trichocarpa. Integration of shotgun sequence assembly with genetic mapping enabled chromosome-scale reconstruction of the genome. More than 45,000 putative protein-coding genes were identified. Analysis of the assembled genome revealed a whole-genome duplication event; about 8000 pairs of duplicated genes from that event survived in the Populus genome. A second, older duplication event is indistinguishably coincident with the divergence of the Populus and Arabidopsis lineages. Nucleotide substitution, tandem gene duplication, and gross chromosomal rearrangement appear to proceed substantially more slowly in Populus than in Arabidopsis. Populus has more protein-coding genes than Arabidopsis, ranging on average from 1.4 to 1.6 putative Populus homologs for each Arabidopsis gene. However, the relative frequency of protein domains in the two genomes is similar. Overrepresented exceptions in Populus include genes associated with lignocellulosic wall biosynthesis, meristem development, disease resistance, and metabolite transport.

The Norway spruce genome sequence and conifer genome evolution

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

Use of Ecotilling as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa

Ecotilling was used as a simple nucleotide polymorphism (SNP) discovery tool to examine DNA variation in natural populations of the western black cottonwood, Populus trichocarpa, and was found to be more efficient than sequencing for large‐scale studies of genetic variation in this tree. A publicly available, live reference collection of P. trichocarpa from the University of British Columbia Botanical Garden was used in this study to survey variation in nine different genes among individuals from 41 different populations. A large amount of genetic variation was detected, but the level of variation appears to be less than in the related species, Populus tremula, based on reported statistics for that tree. Genes examined varied considerably in their level of variation, from PoptrTB1 which had a single SNP, to PoptrLFY which had more than 23 in the 1000‐bp region examined. Overall nucleotide diversity, measured as Total, was relatively low at 0.00184. Linkage disequilibrium, on the other hand, was higher than reported for some woody plant species, with mean r2 equal to 0.34. This study reveals the potential of Ecotilling as a rapid genotype discovery method to explore and utilize the large pool of genetic variation in tree species.

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

BackgroundThe 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.ResultsWe developed a genome sequence resource for the mountain pine beetle to better understand the unique aspects of this insects 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.ConclusionsDespite 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.

Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest.

Plants produce a vast array of specialized metabolites, many of which are used as pharmaceuticals, flavors, fragrances, and other high-value fine chemicals. However, most of these compounds occur in non-model plants for which genomic sequence information is not yet available. The production of a large amount of nucleotide sequence data using next-generation technologies is now relatively fast and cost-effective, especially when using the latest Roche-454 and Illumina sequencers with enhanced base-calling accuracy. To investigate specialized metabolite biosynthesis in non-model plants we have established a data-mining framework, employing next-generation sequencing and computational algorithms, to construct and analyze the transcriptomes of 75 non-model plants that produce compounds of interest for biotechnological applications. After sequence assembly an extensive annotation approach was applied to assign functional information to over 800,000 putative transcripts. The annotation is based on direct searches against public databases, including RefSeq and InterPro. Gene Ontology (GO), Enzyme Commission (EC) annotations and associated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps are also collected. As a proof-of-concept, the selection of biosynthetic gene candidates associated with six specialized metabolic pathways is described. A web-based BLAST server has been established to allow public access to assembled transcriptome databases for all 75 plant species of the PhytoMetaSyn Project (www.phytometasyn.ca).

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

BackgroundThe 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.ResultsGene 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.ConclusionsThe 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.

Synthetic biosystems for the production of high-value plant metabolites

Plants display an immense diversity of specialized metabolites, many of which have been important to humanity as medicines, flavors, fragrances, pigments, insecticides and other fine chemicals. Apparently, much of the variation in plant specialized metabolism evolved through events of gene duplications followed by neo- or sub-functionalization. Most of the catalytic diversity of plant enzymes is unexplored since previous biochemical and genomics efforts have focused on a relatively small number of species. Interdisciplinary research in plant genomics, microbial engineering and synthetic biology provides an opportunity to accelerate the discovery of new enzymes. The massive identification, characterization and cataloguing of plant enzymes coupled with their deployment in metabolically optimized microbes provide a high-throughput functional genomics tool and a novel strain engineering pipeline.

Evolution of Diterpene Metabolism: Sitka Spruce CYP720B4 Catalyzes Multiple Oxidations in Resin Acid Biosynthesis of Conifer Defense against Insects

Diterpene resin acids (DRAs) are specialized (secondary) metabolites of the oleoresin defense of conifers produced by diterpene synthases and cytochrome P450s of the CYP720B family. The evolution of DRA metabolism shares common origins with the biosynthesis of ent-kaurenoic acid, which is highly conserved in general (primary) metabolism of gibberellin biosynthesis. Transcriptome mining in species of spruce (Picea) and pine (Pinus) revealed CYP720Bs of four distinct clades. We cloned a comprehensive set of 12 different Sitka spruce (Picea sitchensis) CYP720Bs as full-length cDNAs. Spatial expression profiles, methyl jasmonate induction, and transcript enrichment in terpenoid-producing resin ducts suggested a role of CYP720B4 in DRA biosynthesis. CYP720B4 was characterized as a multisubstrate, multifunctional enzyme by the formation of oxygenated diterpenoids in metabolically engineered yeast, yeast in vivo transformation of diterpene substrates, in vitro assays with CYP720B4 protein produced in Escherichia coli, and alteration of DRA profiles in RNA interference-suppressed spruce seedlings. CYP720B4 was active with 24 different diterpenoid substrates, catalyzing consecutive C-18 oxidations in the biosynthesis of an array of diterpene alcohols, aldehydes, and acids. CYP720B4 was most active in the formation of dehydroabietic acid, a compound associated with insect resistance of Sitka spruce. We identified patterns of convergent evolution of CYP720B4 in DRA metabolism and ent-kaurene oxidase CYP701 in gibberellin metabolism and revealed differences in the evolution of specialized and general diterpene metabolism in a gymnosperm. The genomic and functional characterization of the gymnosperm CYP720B family highlights that the evolution of specialized metabolism involves substantial diversification relative to conserved, general metabolism.

Wound-Induced Terpene Synthase Gene Expression in Sitka Spruce That Exhibit Resistance or Susceptibility to Attack by the White Pine Weevil

We analyzed the expression pattern of various terpene synthase (TPS) genes in response to a wounding injury applied to the apical leader of Sitka spruce (Picea sitchensis Bong. Carr.) genotypes known to be resistant (R) or susceptible (S) to white pine weevil (Pissodes strobi Peck.) attack. The purpose was to test if differences in constitutive or wound-induced TPS expression can be associated with established weevil resistance. All wounding treatments were conducted on 9-year-old R and S trees growing under natural field conditions within the range of variation for weevil R and S genotypes. Representative cDNAs of the monoterpene synthase (mono-TPS), sesquiterpene synthase (sesqui-TPS), and diterpene synthase (di-TPS) classes were isolated from Sitka spruce to assess TPS transcript levels. Based on amino acid sequence similarity, the cDNAs resemble Norway spruce (Picea abies) (−)-linalool synthase (mono-TPS; PsTPS-Linl) and levopimaradiene/abietadiene synthase (di-TPS; PsTPS-LASl), and grand fir (Abies grandis) δ-selinene synthase (sesqui-TPS; PsTPS-Sell). One other mono-TPS was functionally identified as (−)-limonene synthase (PsTPS-Lim). No significant difference in constitutive expression levels for these TPSs was detected between R and S trees. However, over a postwounding period of 16 d, only R trees exhibited significant transcript accumulation for the mono- and sesqui-TPS tested. Both R and S trees exhibited a significant accumulation of PsTPS-LASl transcripts. An assessment of traumatic resin duct formation in wounded leaders showed that both R and S trees responded by forming traumatic resin ducts; however, the magnitude of this response was significantly greater in R trees. Collectively, our data imply that the induced resinosis response is an important aspect of defense in weevil R Sitka spruce trees growing under natural conditions.

Comparative transcriptome analysis of Arabidopsis thaliana infested by diamond back moth ( Plutella xylostella ) larvae reveals signatures of stress response, secondary metabolism, and signalling

BackgroundPlants are exposed to attack from a large variety of herbivores. Feeding insects can induce substantial changes of the host plant transcriptome. Arabidopsis thaliana has been established as a relevant system for the discovery of genes associated with response to herbivory, including genes for specialized (i.e. secondary) metabolism as well as genes involved in plant-insect defence signalling.ResultsUsing a 70-mer oligonulceotide microarray covering 26,090 gene-specific elements, we monitored changes of the Arabidopsis leaf transcriptome in response to feeding by diamond back moth (DBM; Plutella xylostella) larvae. Analysis of samples from a time course of one hour to 24 hours following onset of DBM feeding revealed almost three thousand (2,881) array elements (including 2,671 genes with AGI annotations) that were differentially expressed (>2-fold; p[t-test] < 0.05) of which 1,686 also changed more than twofold in expression between at least two time points of the time course with p(ANOVA) < 0.05. While the majority of these transcripts were up-regulated within 8 h upon onset of insect feeding relative to untreated controls, cluster analysis identified several distinct temporal patterns of transcriptome changes. Many of the DBM-induced genes fall into ontology groups annotated as stress response, secondary metabolism and signalling. Among DBM-induced genes associated with plant signal molecules or phytohormones, genes associated with octadecanoid signalling were clearly overrepresented. We identified a substantial number of differentially expressed genes associated with signal transduction in response to DBM feeding, and we compared there expression profiles with those of previously reported transcriptome responses induced by other insect herbivores, specifically Pieris rapae, Frankliniella occidentalis, Bemisia tabaci,Myzus persicae, and Brevicoryne brassicae.ConclusionArabidopsis responds to feeding DBM larvae with a drastic reprogramming of the transcriptome, which has considerable overlap with the response induced by other insect herbivores. Based on a meta-analysis of microarray data we identified groups of transcription factors that are either affected by multiple forms of biotic or abiotic stress including DBM feeding or, alternatively, were responsive to DBM herbivory but not to most other forms of stress.