Cytochrome P450 mono-oxygenases in conifer genomes: discovery of members of the terpenoid oxygenase superfamily in spruce and pine

Identification of P450 Candidates Associated with the Biosynthesis of Physalin-Class Compounds in Physalis angulata

The Physalis genus has long been used as traditional medicine in the treatment of various diseases. Physalins, the characteristic class of compounds in this genus, are major bioactive constituents. To date, the biogenesis of physalins remains largely unknown, except for the recently established knowledge that 24-methyldesmosterol is a precursor of physalin. To identify the genes encoding P450s that are putatively involved in converting 24-methyldesmosterol to physalins, a total of 306 P450-encoding unigenes were retrieved from our recently constructed P. angulata transcriptome. Extensive phylogenetic analysis proposed 21 P450s that might participate in physalin biosynthesis. To validate the candidates, we developed a virus-induced gene silencing (VIGS) system for P. angulata, and four P450 candidates were selected for the VIGS experiments. The reduction in the transcripts of the four P450 candidates by VIGS all led to decreased levels of physalin-class compounds in the P. angulata leaves. Thus, this study provides a number of P450 candidates that are likely associated with the biosynthesis of physalin-class compounds, forming a strong basis to reveal the unknown physalin biosynthetic pathway in the future.

Combining QTL Mapping and Transcriptomics to Decipher the Genetic Architecture of Phenolic Compounds Metabolism in the Conifer White Spruce

Conifer forests worldwide are becoming increasingly vulnerable to the effects of climate change. Although the production of phenolic compounds (PCs) has been shown to be modulated by biotic and abiotic stresses, the genetic basis underlying the variation in their constitutive production level remains poorly documented in conifers. We used QTL mapping and RNA-Seq to explore the complex polygenic network underlying the constitutive production of PCs in a white spruce (Picea glauca) full-sib family for 2 years. QTL detection was performed for nine PCs and differentially expressed genes (DEGs) were identified between individuals with high and low PC contents for five PCs exhibiting stable QTLs across time. A total of 17 QTLs were detected for eight metabolites, including one major QTL explaining up to 91.3% of the neolignan-2 variance. The RNA-Seq analysis highlighted 50 DEGs associated with phenylpropanoid biosynthesis, several key transcription factors, and a subset of 137 genes showing opposite expression patterns in individuals with high levels of the flavonoids gallocatechin and taxifolin glucoside. A total of 19 DEGs co-localized with QTLs. Our findings represent a significant step toward resolving the genomic architecture of PC production in spruce and facilitate the functional characterization of genes and transcriptional networks responsible for differences in constitutive production of PCs in conifers.

High marker density GWAS provides novel insights into the genomic architecture of terpene oil yield in Eucalyptus

Terpenoid-based essential oils are economically important commodities, yet beyond their biosynthetic pathways, little is known about the genetic architecture of terpene oil yield from plants. Transport, storage, evaporative loss, transcriptional regulation and precursor competition may be important contributors to this complex trait. Here, we associate 2.39 million single nucleotide polymorphisms derived from shallow whole-genome sequencing of 468 Eucalyptus polybractea individuals with 12 traits related to the overall terpene yield, eight direct measures of terpene concentration and four biomass-related traits. Our results show that in addition to terpene biosynthesis, development of secretory cavities, where terpenes are both synthesized and stored, and transport of terpenes were important components of terpene yield. For sesquiterpene concentrations, the availability of precursors in the cytosol was important. Candidate terpene synthase genes for the production of 1,8-cineole and α-pinene, and β-pinene (which comprised > 80% of the total terpenes) were functionally characterized as a 1,8-cineole synthase and a β/α-pinene synthase. Our results provide novel insights into the genomic architecture of terpene yield and we provide candidate genes for breeding or engineering of crops for biofuels or the production of industrially valuable terpenes.

Four terpene synthases contribute to the generation of chemotypes in tea tree (Melaleuca alternifolia)

BACKGROUND

Terpene rich leaves are a characteristic of Myrtaceae. There is significant qualitative variation in the terpene profile of plants within a single species, which is observable as "chemotypes". Understanding the molecular basis of chemotypic variation will help explain how such variation is maintained in natural populations as well as allowing focussed breeding for those terpenes sought by industry. The leaves of the medicinal tea tree, Melaleuca alternifolia, are used to produce terpinen-4-ol rich tea tree oil, but there are six naturally occurring chemotypes; three cardinal chemotypes (dominated by terpinen-4-ol, terpinolene and 1,8-cineole, respectively) and three intermediates. It has been predicted that three distinct terpene synthases could be responsible for the maintenance of chemotypic variation in this species.

RESULTS

We isolated and characterised the most abundant terpene synthases (TPSs) from the three cardinal chemotypes of M. alternifolia. Functional characterisation of these enzymes shows that they produce the dominant compounds in the foliar terpene profile of all six chemotypes. Using RNA-Seq, we investigated the expression of these and 24 additional putative terpene synthases in young leaves of all six chemotypes of M. alternifolia.

CONCLUSIONS

Despite contributing to the variation patterns observed, variation in gene expression of the three TPS genes is not enough to explain all variation for the maintenance of chemotypes. Other candidate terpene synthases as well as other levels of regulation must also be involved. The results of this study provide novel insights into the complexity of terpene biosynthesis in natural populations of a non-model organism.

Transcriptome Assembly and Systematic Identification of Novel Cytochrome P450s in Taxus chinensis

Taxus spp. is a highly valuable medicinal plant with multiple pharmacological effects on various cancers. Cytochrome P450s (CYP450s) play important roles in the biosynthesis of active compounds in Taxus spp., such as the famous diterpenoid, Taxol. However, some specific CYP450 enzymes involved in the biosynthesis of Taxol remain unknown, and the systematic identification of CYP450s in Taxus has not been reported. In this study, 118 full-length and 175 partial CYP450 genes were identified in Taxus chinensis transcriptomes. The 118 full-length genes were divided into 8 clans and 29 families. The CYP71 clan included all A-type genes (52) belonging to 11 families. The other seven clans possessed 18 families containing 66 non-A-type genes. Two new gymnosperm-specific families were discovered, and were named CYP864 and CYP947 respectively. Protein sequence alignments revealed that all of the T. chinensis CYP450s hold distinct conserved domains. The expression patterns of all 118 CYP450 genes during the long-time subculture and MeJA elicitation were analyzed. Additionally, the expression levels of 15 novel CYP725 genes in different Taxus species were explored. Considering all the evidence, 6 CYP725s were identified to be candidates for Taxol biosynthesis. The cis-regulatory elements involved in the transcriptional regulation were also identified in the promoter regions of CYP725s. This study presents a comprehensive overview of the CYP450 gene family in T. chinensis and can provide important insights into the functional gene studies of Taxol 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.

Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families

The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.

Genus Cistus: a model for exploring labdane-type diterpenes' biosynthesis and a natural source of high value products with biological, aromatic, and pharmacological properties

The family Cistaceae (Angiosperm, Malvales) consists of 8 genera and 180 species, with 5 genera native to the Mediterranean area (Cistus, Fumara, Halimium, Helianthemum, and Tuberaria). Traditionally, a number of Cistus species have been used in Mediterranean folk medicine as herbal tea infusions for healing digestive problems and colds, as extracts for the treatment of diseases, and as fragrances. The resin, ladano, secreted by the glandular trichomes of certain Cistus species contains a number of phytochemicals with antioxidant, antibacterial, antifungal, and anticancer properties. Furthermore, total leaf aqueous extracts possess anti-influenza virus activity. All these properties have been attributed to phytochemicals such as terpenoids, including diterpenes, labdane-type diterpenes and clerodanes, phenylpropanoids, including flavonoids and ellagitannins, several groups of alkaloids and other types of secondary metabolites. In the past 20 years, research on Cistus involved chemical, biological and phylogenetic analyses but recent investigations have involved genomic and molecular approaches. Our lab is exploring the biosynthetic machinery that generates terpenoids and phenylpropanoids, with a goal to harness their numerous properties that have applications in the pharmaceutical, chemical and aromatic industries. This review focuses on the systematics, botanical characteristics, geographic distribution, chemical analyses, biological function and biosynthesis of major compounds, as well as genomic analyses and biotechnological approaches of the main Cistus species found in the Mediterranean basin, namely C. albidus, C. creticus, C. crispus, C. parviflorus, C. monspeliensis, C. populifolius, C. salviifolius, C. ladanifer, C. laurifolius, and C. clusii.

Mass spectrometric analysis of differentially expressed proteins in an endangered medicinal herb, Picrorhiza kurroa

Picrorhiza kurroa grown in the Northwestern Himalayan region is used in various herbal formulations but extensive harvesting of this plant has led it to near extinction. The active constituents responsible for the medicinal properties of P. kurroa have been identified as picroside-I and picroside-II which are present in a particular ratio (1:1.5) in herbal formulations like Picroliv. The biosynthetic pathway of picrosides has been partially deciphered till date and needs to be elucidated completely. Review of literature revealed that no information is available as of today on the proteome analysis of Picrorhiza kurroa w.r.t. picroside-II biosynthesis. Therefore, with the aim of identifying proteins associated with picroside biosynthesis in Picrorhiza kurroa, differential protein expression was studied under picroside accumulating versus nonaccumulating conditions using SDS-PAGE. A total of 19 differentially expressed proteins were identified using MALDI-TOF/TOF MS followed by MASCOT search. Proteins involved in diverse functions were identified amongst which the most important proteins were glyceraldehyde-3-phosphate dehydrogenase, 1-aminocyclopropane-1-carboxylate oxidase, photosystem I reaction centre subunit V, 2-oxoglutarate ferrous-dependent oxygenase and putative cytochrome P450 superfamily protein because of their role in picroside biosynthesis. These identified proteins provide an insight and a basic platform for thorough understanding of biosynthesis of secondary metabolites and various other physiological processes of P. kurroa.

Molecular cloning and characterization of a cytochrome P450 taxoid 9á-hydroxylase in Ginkgo biloba cells

Taxol is a well-known effective anticancer compound. Due to the inability to synthesize sufficient quantities of taxol to satisfy commercial demand, a biotechnological approach for a large-scale cell or cell-free system for its production is highly desirable. Several important genes in taxol biosynthesis are currently still unknown and have been shown to be difficult to isolate directly from Taxus, including the gene encoding taxoid 9α-hydroxylase. Ginkgo biloba suspension cells exhibit taxoid hydroxylation activity and provides an alternate means of identifying genes encoding enzymes with taxoid 9α-hydroxylation activity. Through analysis of high throughput RNA sequencing data from G. biloba, we identified two candidate genes with high similarity to Taxus CYP450s. Using in vitro cell-free protein synthesis assays and LC-MS analysis, we show that one candidate that belongs to the CYP716B, a subfamily whose biochemical functions have not been previously studied, possessed 9α-hydroxylation activity. This work will aid future identification of the taxoid 9α-hydroxylase gene from Taxus sp.

Identification of genes in Thuja plicata foliar terpenoid defenses

Thuja plicata (western redcedar) is a long-lived conifer species whose foliage is rarely affected by disease or insect pests, but can be severely damaged by ungulate browsing. Deterrence to browsing correlates with high foliar levels of terpenoids, in particular the monoterpenoid α-thujone. Here, we set out to identify genes whose products may be involved in the production of α-thujone and other terpenoids in this species. First, we generated a foliar transcriptome database from which to draw candidate genes. Second, we mapped the storage of thujones and other terpenoids to foliar glands. Third, we used global expression profiling to identify more than 600 genes that are expressed at high levels in foliage with glands, but can either not be detected or are expressed at low levels in a natural variant lacking foliar glands. Fourth, we used in situ RNA hybridization to map the expression of a putative monoterpene synthase to the epithelium of glands and used enzyme assays with recombinant protein of the same gene to show that it produces sabinene, the monoterpene precursor of α-thujone. Finally, we identified candidate genes with predicted enzymatic functions for the conversion of sabinene to α-thujone. Taken together, this approach generated both general resources and detailed functional characterization in the identification of genes of foliar terpenoid biosynthesis in T. plicata.

Plant P450s as versatile drivers for evolution of species-specific chemical diversity

The irreversible nature of reactions catalysed by P450s makes these enzymes landmarks in the evolution of plant metabolic pathways. Founding members of P450 families are often associated with general (i.e. primary) metabolic pathways, restricted to single copy or very few representatives, indicative of purifying selection. Recruitment of those and subsequent blooms into multi-member gene families generates genetic raw material for functional diversification, which is an inherent characteristic of specialized (i.e. secondary) metabolism. However, a growing number of highly specialized P450s from not only the CYP71 clan indicate substantial contribution of convergent and divergent evolution to the observed general and specialized metabolite diversity. We will discuss examples of how the genetic and functional diversification of plant P450s drives chemical diversity in light of plant evolution. Even though it is difficult to predict the function or substrate of a P450 based on sequence similarity, grouping with a family or subfamily in phylogenetic trees can indicate association with metabolism of particular classes of compounds. Examples will be given that focus on multi-member gene families of P450s involved in the metabolic routes of four classes of specialized metabolites: cyanogenic glucosides, glucosinolates, mono- to triterpenoids and phenylpropanoids.

Geranylgeraniol oxidase activity involved in oxidative formation of geranylgeranoic acid in human hepatoma cells

Geranylgeranoic acid (GGA), a 20-carbon acyclic polyprenoic acid (all-trans 3,7,11,15-tetramethyl- 2,4,6,10,14-hexadecatetraenoic acid) and its derivatives were developed as synthetic "acyclic retinoids" for cancer chemoprevention. Previously, we have shown the natural occurrence of GGA in various medicinal herbs and reported enzymatic formation of GGA from geranylgeraniol (GGOH) through geranylgeranial (GGal) by rat liver homogenates. Here, we present several lines of evidence that a putative GGOH oxidase is involved in GGA synthesis by human hepatoma cell lysates. First, conversion of GGOH to GGal did not require exogenous NAD(+), whereas the conversion from GGal to GGA absolutely required additional NAD(+). Second, GGal synthesis from GGOH was coupled with consumption of oxygen from the reaction mixture. Third, GGOH-dependent GGal synthesis activity was proteinase K-resistant and even enhanced by proteinase K treatment; GGOH oxidase activity was enriched in the mitochondrial fraction. Finally, recombinant human monoamine oxidase (MAO)-B, but not MAO-A catalyzed oxidation of GGOH to GGal. These data suggest that a putative mitochondrial GGOH oxidase is involved in the initial step of GGA synthesis from GGOH.

Genome-wide analysis of a TaLEA-introduced transgenic Populus simonii × Populus nigra dwarf mutant

A dwarf mutant (dwf1) was obtained among 15 transgenic lines, when TaLEA (Tamarix androssowii late embryogenesis abundant gene) was introduced into Populus simonii × Populus nigra by Agrobacterium tumefaciens-mediated transformation. Under the same growth conditions, dwf1 height was significantly reduced compared with the wild type and the other transgenic lines. Because only one transgenic line (dwf1) displayed the dwarf phenotype, we considered that T-DNA insertion sites may play a role in the mutant formation. The mechanisms underlying this effect were investigated using TAIL-PCR (thermal asymmetric interlaced PCR) and microarrays methods. According to the TAIL-PCR results, two flanking sequences located on chromosome IV and VIII respectively, were cloned. The results indicated the integration of two independent T-DNA copies. We searched for the potential genes near to the T-DNA insertions. The nearest gene was a putative poplar AP2 transcription factor (GI: 224073210). Expression analysis showed that AP2 was up-regulated in dwf1 compared with the wild type and the other transgenic lines. According to the microarrays results, a total of 537 genes involved in hydrolase, kinase and transcription factor activities, as well as protein and nucleotide binding, showed significant alterations in gene expression. These genes were expressed in more than 60 metabolic pathways, including starch, sucrose, galactose and glycerolipid metabolism and phenylpropanoids and flavonoid biosyntheses. Our transcriptome and T-DNA insertion sites analyses might provide some useful insights into the dwarf mutant formation.

Bifunctional cis-abienol synthase from Abies balsamea discovered by transcriptome sequencing and its implications for diterpenoid fragrance production

The labdanoid diterpene alcohol cis-abienol is a major component of the aromatic oleoresin of balsam fir (Abies balsamea) and serves as a valuable bioproduct material for the fragrance industry. Using high-throughput 454 transcriptome sequencing and metabolite profiling of balsam fir bark tissue, we identified candidate diterpene synthase sequences for full-length cDNA cloning and functional characterization. We discovered a bifunctional class I/II cis-abienol synthase (AbCAS), along with the paralogous levopimaradiene/abietadiene synthase and isopimaradiene synthase, all of which are members of the gymnosperm-specific TPS-d subfamily. The AbCAS-catalyzed formation of cis-abienol proceeds via cyclization and hydroxylation at carbon C-8 of a postulated carbocation intermediate in the class II active site, followed by cleavage of the diphosphate group and termination of the reaction sequence without further cyclization in the class I active site. This reaction mechanism is distinct from that of synthases of the isopimaradiene- or levopimaradiene/abietadiene synthase type, which employ deprotonation reactions in the class II active site and secondary cyclizations in the class I active site, leading to tricyclic diterpenes. Comparative homology modeling suggested the active site residues Asp-348, Leu-617, Phe-696, and Gly-723 as potentially important for the specificity of AbCAS. As a class I/II bifunctional enzyme, AbCAS is a promising target for metabolic engineering of cis-abienol production.

Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway

Background

In trees, a substantial amount of carbon is directed towards production of phenolics for development and defense. This metabolic pathway is also a major factor in resistance to insect pathogens in spruce. In such gene families, environmental stimuli may have an important effect on the evolutionary fate of duplicated genes, and different expression patterns may indicate functional diversification.

Results

Gene families in spruce (Picea) have expanded to superfamilies, including O-methyltransferases, cytochrome-P450, and dirigents/classIII-peroxidases. Neo-functionalization of superfamily members from different clades is reflected in expression diversification. Genetical genomics can provide new insights into the genetic basis and evolution of insect resistance in plants. Adopting this approach, we merged genotype data (252 SNPs in a segregating pedigree), gene expression levels (for 428 phenylpropanoid-related genes) and measures of susceptibility to Pissodes stobi, using a partial-diallel crossing-design with white spruce (Picea glauca). Thirty-eight expressed phenylpropanoid-related genes co-segregated with weevil susceptibility, indicating either causative or reactive effects of these genes to weevil resistance. We identified eight regulatory genomic regions with extensive overlap of quantitative trait loci from susceptibility and growth phenotypes (pQTLs) and expression QTL (eQTL) hotspots. In particular, SNPs within two different CCoAOMT loci regulate phenotypic variation from a common set of 24 genes and three resistance traits.

Conclusions

Pest resistance was associated with individual candidate genes as well as with trans-regulatory hotspots along the spruce genome. Our results showed that specific genes within the phenylpropanoid pathway have been duplicated and diversified in the conifer in a process fundamentally different from short-lived angiosperm species. These findings add to the information about the role of the phenylpropanoid pathway in the evolution of plant defense mechanisms against insect pests and provide substantial potential for the functional characterization of several not yet resolved alternative pathways in plant defenses.

Induction of isoprenyl diphosphate synthases, plant hormones and defense signalling genes correlates with traumatic resin duct formation in Norway spruce (Picea abies)

Norway spruce (Picea abies) defends itself against herbivores and pathogens by formation of traumatic resin ducts filled with terpenoid-based oleoresin. An important group of enzymes in terpenoid biosynthesis are the short-chain isoprenyl diphosphate synthases which produce geranyl diphosphate (C(10)), farnesyl diphosphate (C(15)), and geranylgeranyl diphosphate (C(20)) as precursors of monoterpenes, sesquiterpenes, and diterpene resin acids, respectively. After treatment with methyl jasmonate (MJ) we investigated the expression of all isoprenyl diphosphate synthase genes characterized to date from Norway spruce and correlated this with formation of traumatic resin ducts and terpene accumulation. Formation of traumatic resin ducts correlated with higher amounts of monoterpenes, sesquiterpenes and diterpene resin acids and an upregulation of isoprenyl diphosphate synthase genes producing geranyl diphosphate or geranylgeranyl diphosphate. Among defense hormones, jasmonate and jasmonate-isoleucine conjugate accumulated to higher levels in trees with extensive traumatic resin duct formation, whereas salicylate did not. Jasmonate and ethylene are likely to both be involved in formation of traumatic resin ducts based on elevated transcripts of genes encoding lipoxygenase and 1-aminocyclopropane-1-carboxylic acid oxidase associated with resin duct formation. Other genes involved in defense signalling in other systems, mitogen-activated protein kinase3 and nonexpressor of pathogenesis-related gene1, were also associated with traumatic resin duct formation. These responses were detected not only at the site of MJ treatment, but also systemically up to 60 cm above the site of treatment on the trunk.

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.

Cytochromes p450

There are 244 cytochrome P450 genes (and 28 pseudogenes) in the Arabidopsis genome. P450s thus form one of the largest gene families in plants. Contrary to what was initially thought, this family diversification results in very limited functional redundancy and seems to mirror the complexity of plant metabolism. P450s sometimes share less than 20% identity and catalyze extremely diverse reactions leading to the precursors of structural macromolecules such as lignin, cutin, suberin and sporopollenin, or are involved in biosynthesis or catabolism of all hormone and signaling molecules, of pigments, odorants, flavors, antioxidants, allelochemicals and defense compounds, and in the metabolism of xenobiotics. The mechanisms of gene duplication and diversification are getting better understood and together with co-expression data provide leads to functional characterization.

Global transcriptome analysis of constitutive resistance to the white pine weevil in spruce

Constitutive defense mechanisms are critical to the understanding of defense mechanisms in conifers because they constitute the first barrier to attacks by insect pests. In interior spruce, trees that are putatively resistant and susceptible to attacks by white pine weevil (Pissodes strobi) typically exhibit constitutive differences in traits such as resin duct size and number, bark thickness, and terpene content. To improve our knowledge of their genetic basis, we compared globally the constitutive expression levels of 17,825 genes between 20 putatively resistant and 20 putatively susceptible interior spruce trees from the British Columbia tree improvement program. We identified 54 upregulated and 137 downregulated genes in resistant phenotypes, relative to susceptible phenotypes, with a maximum fold change of 2.24 and 3.91, respectively. We found a puzzling increase of resistance by downregulated genes, as one would think that "procuring armaments" is the best defense. Also, although terpenes and phenolic compounds play an important role in conifer defense, we found few of these genes to be differentially expressed. We found 15 putative small heat-shock proteins (sHSP) and several other stress-related proteins to be downregulated in resistant trees. Downregulated putative sHSP belong to several sHSP classes and represented 58% of all tested putative sHSP. These proteins are well known to be involved in plant response to various kinds of abiotic stress; however, their role in constitutive resistance is not yet understood. The lack of correspondence between transcriptome profile clusters and phenotype classifications suggests that weevil resistance in spruce is a complex trait.

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.

Polished rice as natural sources of cancer-preventing geranylgeranoic acid

Geranylgeranoic acid, a 20-carbon polyprenoic acid (all-trans 3,7,11,15-tetramethyl-2,4,6,10,14-hexadecatetraenoic acid) and its derivatives were previously developed as synthetic “acyclic retinoids” for cancer chemoprevention. Recently, we demonstrated the natural occurrence of geranylgeranoic acid in various medicinal herbs (Shidoji and Ogawa, 2004). In this present study, we present several lines of evidence to demonstrate that geranylgeranyl diphosphate taken in foods could be metabolized to GGA through geranylgeraniol and geranylgeranyl aldehyde via the following steps: 1) The conversion from geranylgeranyl diphosphate to geranylgeraniol was demonstrated to occur by the action of bovine intestinal alkaline phosphatase, with a K_m of 46.1 µM. 2) Geranylgeraniol oxidase-mediated conversion of geranylgeraniol to geranylgeranyl aldehyde was revealed in rat liver homogenates, which activity was mainly localized in the mitochondrial fraction. The mitochondrial enzyme showed a K_m of 92.9 µM. 3) The conversion of geranylgeranyl aldehyde to geranylgeranoic acid by geranylgeranyl aldehyde dehydrogenase in rat liver homogenates was absolutely dependent on exogenously added NAD⁺ or NADP⁺. The K_m of the mitochondrial geranylgeranyl aldehyde dehydrogenase was 27.5 µM for geranylgeranyl aldehyde. Taken together, our data suggest that cancer preventive geranylgeranoic acid could be a physiological metabolite from commonly consumed foods.

Terpenoid biosynthesis and specialized vascular cells of conifer defense

Defense-related terpenoid biosynthesis in conifers is a dynamic process closely associated with specialized anatomical structures that allows conifers to cope with attack from many potential pests and pathogens. The constitutive and inducible terpenoid defense of conifers involves several hundred different monoterpenes, sesquiterpenes and diterpenes. Changing arrays of these many compounds are formed from the general isoprenoid pathway by activities of large gene families for two classes of enzymes, the terpene synthases and the cytochrome P450-dependent monooxygenases of the CYP720B group. Extensive studies have been conducted on the genomics, proteomics and molecular biochemical characterization of these enzymes. Many of the conifer terpene synthases are multi-product enzymes, and the P450 enzymes of the CYP720B group are promiscuous in catalyzing multiple oxidations, along homologous series of diterpenoids, from a broad spectrum of substrates. The terpene synthases and CYP720B genes respond to authentic or simulated insect attack with increased transcript levels, protein abundance and enzyme activity. The constitutive and induced oleoresin terpenoids for conifer defense accumulate in preformed cortical resin ducts and in xylem trauma-associated resin ducts. Formation of these resin ducts de novo in the cambium zone and developing xylem, following insect attack or treatment of trees with methyl jasmonate, is a unique feature of the induced defense of long-lived conifer trees.

Diversification of P450 genes during land plant evolution

Plant cytochromes P450 (P450s) catalyze a wide variety of monooxygenation/hydroxylation reactions in primary and secondary metabolism. The number of P450 genes in plant genomes is estimated to be up to 1% of total gene annotations of each plant species. This implies that diversification within P450 gene superfamilies has led to the emergence of new metabolic pathways throughout land plant evolution. The conserved P450 families contribute to chemical defense mechanisms under terrestrial conditions and several are involved in hormone biosynthesis and catabolism. Species-specific P450 families are essential for the biosynthetic pathways of species-specialized metabolites. Future genome-wide analyses of P450 gene clusters and coexpression networks should help both in identifying the functions of many orphan P450s and in understanding the evolution of this versatile group of enzymes.

Targeted isolation, sequence assembly and characterization of two white spruce (Picea glauca) BAC clones for terpenoid synthase and cytochrome P450 genes involved in conifer defence reveal insights into a conifer genome

BACKGROUND

Conifers are a large group of gymnosperm trees which are separated from the angiosperms by more than 300 million years of independent evolution. Conifer genomes are extremely large and contain considerable amounts of repetitive DNA. Currently, conifer sequence resources exist predominantly as expressed sequence tags (ESTs) and full-length (FL)cDNAs. There is no genome sequence available for a conifer or any other gymnosperm. Conifer defence-related genes often group into large families with closely related members. The goals of this study are to assess the feasibility of targeted isolation and sequence assembly of conifer BAC clones containing specific genes from two large gene families, and to characterize large segments of genomic DNA sequence for the first time from a conifer.

RESULTS

We used a PCR-based approach to identify BAC clones for two target genes, a terpene synthase (3-carene synthase; 3CAR) and a cytochrome P450 (CYP720B4) from a non-arrayed genomic BAC library of white spruce (Picea glauca). Shotgun genomic fragments isolated from the BAC clones were sequenced to a depth of 15.6- and 16.0-fold coverage, respectively. Assembly and manual curation yielded sequence scaffolds of 172 kbp (3CAR) and 94 kbp (CYP720B4) long. Inspection of the genomic sequences revealed the intron-exon structures, the putative promoter regions and putative cis-regulatory elements of these genes. Sequences related to transposable elements (TEs), high complexity repeats and simple repeats were prevalent and comprised approximately 40% of the sequenced genomic DNA. An in silico simulation of the effect of sequencing depth on the quality of the sequence assembly provides direction for future efforts of conifer genome sequencing.

CONCLUSION

We report the first targeted cloning, sequencing, assembly, and annotation of large segments of genomic DNA from a conifer. We demonstrate that genomic BAC clones for individual members of multi-member gene families can be isolated in a gene-specific fashion. The results of the present work provide important new information about the structure and content of conifer genomic DNA that will guide future efforts to sequence and assemble conifer genomes.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis)

Background

Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation.

Results

As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs.

Conclusion

This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

Terpenoid biomaterials

Terpenoids (isoprenoids) encompass more than 40 000 structures and form the largest class of all known plant metabolites. Some terpenoids have well-characterized physiological functions that are common to most plant species. In addition, many of the structurally diverse plant terpenoids may function in taxonomically more discrete, specialized interactions with other organisms. Historically, specialized terpenoids, together with alkaloids and many of the phenolics, have been referred to as secondary metabolites. More recently, these compounds have become widely recognized, conceptually and/or empirically, for their essential ecological functions in plant biology. Owing to their diverse biological activities and their diverse physical and chemical properties, terpenoid plant chemicals have been exploited by humans as traditional biomaterials in the form of complex mixtures or in the form of more or less pure compounds since ancient times. Plant terpenoids are widely used as industrially relevant chemicals, including many pharmaceuticals, flavours, fragrances, pesticides and disinfectants, and as large-volume feedstocks for chemical industries. Recently, there has been a renaissance of awareness of plant terpenoids as a valuable biological resource for societies that will have to become less reliant on petrochemicals. Harnessing the powers of plant and microbial systems for production of economically valuable plant terpenoids requires interdisciplinary and often expensive research into their chemistry, biosynthesis and genomics, as well as metabolic and biochemical engineering. This paper provides an overview of the formation of hemi-, mono-, sesqui- and diterpenoids in plants, and highlights some well-established examples for these classes of terpenoids in the context of biomaterials and biofuels.

Metabolic diversification--independent assembly of operon-like gene clusters in different plants

Operons are clusters of unrelated genes with related functions that are a feature of prokaryotic genomes. Here, we report on an operon-like gene cluster in the plant Arabidopsis thaliana that is required for triterpene synthesis (the thalianol pathway). The clustered genes are coexpressed, as in bacterial operons. However, despite the resemblance to a bacterial operon, this gene cluster has been assembled from plant genes by gene duplication, neofunctionalization, and genome reorganization, rather than by horizontal gene transfer from bacteria. Furthermore, recent assembly of operon-like gene clusters for triterpene synthesis has occurred independently in divergent plant lineages (Arabidopsis and oat). Thus, selection pressure may act during the formation of certain plant metabolic pathways to drive gene clustering.

Cloning and characterization of isoprenyl diphosphate synthases with farnesyl diphosphate and geranylgeranyl diphosphate synthase activity from Norway spruce (Picea abies) and their relation to induced oleoresin formation

The conifer Picea abies (Norway spruce) employs terpenoid-based oleoresins as part of its constitutive and induced defense responses to herbivores and pathogens. The isoprenyl diphosphate synthases are branch-point enzymes of terpenoid biosynthesis leading to the various terpene classes. We isolated three genes encoding isoprenyl diphosphate synthases from P. abies cDNA libraries prepared from the bark and wood of methyl jasmonate-treated saplings and screened via a homology-based PCR approach using degenerate primers. Enzyme assays of the purified recombinant proteins expressed in Escherichia coli demonstrated that one gene (PaIDS 4) encodes a farnesyl diphosphate synthase and the other two (PaIDS 5 and PaIDS 6) encode geranylgeranyl diphosphate synthases. The sequences have moderate similarity to those of farnesyl diphosphate and geranylgeranyl diphosphate synthases already known from plants, and the kinetic properties of the enzymes are not unlike those of other isoprenyl diphosphate synthases. Of the three genes, only PaIDS 5 displayed a significant increase in transcript level in response to methyl jasmonate spraying, suggesting its involvement in induced oleoresin biosynthesis.