Search

Substantial growth losses in Douglas-fir can be the result of Swiss needle cast, a foliar disease caused by the fungus Nothophaeocryptopus gaeumannii. Although the disease is native to western North America, it is expected to become increasingly problematic in regions where climates become warmer in winter and wetter in spring, such as in coastal British Columbia (BC), Canada. Previous research suggests that tolerance to this disease is under partial genetic control. We therefore aim to screen for tolerance and resistance to the disease in the tree breeding population of coastal Douglas-fir (Pseudotsuga menziesii var. menziesii) in BC. We evaluated if early screening for resistance or tolerance to this disease is possible. We worked with 32 families grown for 18-years on two full-sibling genetic field trial sites representing different climates. We assessed >900 trees for disease signs and symptoms from 2017 to 2019. Needle retention was assessed in the field. In the laboratory, the proportion of stomata occluded with pseudothecia was measured, deoxyribonucleic acid (DNA) was extracted, and the proportion of fungal DNA in the needles was assessed using quantitative polymerase chain reaction. Strong differences were observed among families for needle retention, stomatal occlusion and fungal load. The milder and wetter site closer to the coast was generally more severely affected. Families showed rank changes between the two sites for all response variables. Higher needle retention was correlated with increased tree volume. No correlation was found between the proportion of stomata occluded with pseudothecia, fungal DNA load and needle retention. These results indicate that a more complex pathology is involved in causing needle loss. We conclude that screening for Swiss needle cast tolerance in the coastal BC Douglas-fir breeding population is possible if needle retention can be assessed and that area-specific deployment strategies may be needed given family rank changes in different environments.

The salicinoids are anti-herbivore phenolic glycosides unique to the Salicaceae (Populus and Salix). They consist of a salicyl alcohol glucoside core, which is usually further acylated with benzoic, cinnamic or phenolic acids. While salicinoid structures are well known, their biosynthesis remains enigmatic. Recently, two enzymes from poplar, salicyl alcohol benzoyl transferase and benzyl alcohol benzoyl transferase, were shown to catalyze the production of salicyl benzoate, a predicted potential intermediate in salicinoid biosynthesis. Here, we used transcriptomics and co-expression analysis with these two genes to identify two UDP-glucose-dependent glycosyltransferases (UGT71L1 and UGT78M1) as candidate enzymes in this pathway. Both recombinant enzymes accepted only salicyl benzoate, salicylaldehyde and 2-hydroxycinnamic acid as glucose acceptors. Knocking out the UGT71L1 gene by CRISPR/Cas9 in poplar hairy root cultures led to the complete loss of salicortin, tremulacin and tremuloidin, and a partial reduction of salicin content. This demonstrated that UGT71L1 is required for synthesis of the major salicinoids, and suggested that an additional route can lead to salicin. CRISPR/Cas9 knockouts for UGT78M1 were not successful, and its in vivo role thus remains to be determined. Although it has a similar substrate preference and predicted structure as UGT71L1, it appears not to contribute to the synthesis of salicortin, tremulacin and tremuloidin, at least in roots. The demonstration of UGT71L1 as an enzyme of salicinoid biosynthesis will open up new avenues for the elucidation of this pathway.

Douglas-fir transcriptomics. Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) is economically important with extensive breeding programs and seed trade. However, the molecular genetics of its seed development are largely unknown. We developed a transcriptome resource covering key developmental stages of megagametophytes over time: prefertilization, fertilization, embryogenesis, and early, unfertilized abortion. RNA sequencing reads were assembled de novo into 105,505 predicted high-confidence transcripts derived from 34,521 predicted genes. Expression levels were estimated based on alignment of the original reads to the reference. Megagametophytes express a distinct set of genes compared to those of vegetative tissues. Transcripts related to signaling, protein turnover, and RNA biogenesis have lower expression values in vegetative tissues, whereas cell wall remodeling, solute transport, and seed storage protein transcripts have higher expression values in megagametophytes. Seed storage protein transcripts become very abundant in both pollinated and unpollinated megagametophytes over time, even in aborting ovules. However, the absence of protein storage bodies in unfertilized megagametophytes suggests extensive posttranscriptional mechanisms that either inhibit storage protein translation or their aggregation into protein bodies. This novel transcriptome resource provides a foundation for further important insights into conifer seed development.

This work was funded by the Strategic Project Partnership Grants Program of NSERC and the Agence Nationale de Recherche - Programme Blanc 456 International.

Occurrence of stomata on both leaf surfaces (amphistomaty) promotes higher stomatal conductance and photosynthesis while simultaneously increasing exposure to potential disease agents in black cottonwood (Populus trichocarpa). A genome-wide association study (GWAS) with 2.2M single nucleotide polymorphisms generated through whole-genome sequencing found 280 loci associated with variation in adaxial stomatal traits, implicating genes regulating stomatal development and behavior. Strikingly, numerous loci regulating plant growth and response to biotic and abiotic stresses were also identified. The most significant locus was a poplar homologue of SPEECHLESS (PtSPCH1). Individuals possessing PtSPCH1 alleles associated with greater adaxial stomatal density originated primarily from environments with shorter growing seasons (e.g. northern latitudes, high elevations) or with less precipitation. PtSPCH1 was expressed in developing leaves but not developing stem xylem. In developing leaves, RNA sequencing showed patterns of coordinated expression between PtSPCH1 and other GWAS-identified genes. The breadth of our GWAS results suggests that the evolution of amphistomaty is part of a larger, complex response in plants. Suites of genes underpin this response, retrieved through genetic association to adaxial stomata, and show coordinated expression during development. We propose that the occurrence of amphistomaty in P. trichocarpa involves PtSPCH1 and reflects selection for supporting rapid growth over investment in immunity.

The shikimate pathway synthesizes aromatic amino acids essential for protein biosynthesis. Shikimate dehydrogenase (SDH) is a central enzyme of this primary metabolic pathway, producing shikimate. The structurally similar quinate is a secondary metabolite synthesized by quinate dehydrogenase (QDH). SDH and QDH belong to the same gene family, which diverged into two phylogenetic clades after a defining gene duplication just prior to the angiosperm/gymnosperm split. Non-seed plants that diverged before this duplication harbour only a single gene of this family. Extant representatives from the chlorophytes (Chlamydomonas reinhardtii), bryophytes (Physcomitrella patens) and lycophytes (Selaginella moellendorfii) encoded almost exclusively SDH activity in vitro. A reconstructed ancestral sequence representing the node just prior to the gene duplication also encoded SDH activity. Quinate dehydrogenase activity was gained only in seed plants following gene duplication. Quinate dehydrogenases of gymnosperms, represented here by Pinus taeda, may be reminiscent of an evolutionary intermediate since they encode equal SDH and QDH activities. The second copy in P. taeda maintained specificity for shikimate similar to the activity found in the angiosperm SDH sister clade. The codon for a tyrosine residue within the active site displayed a signature of positive selection at the node defining the QDH clade, where it changed to a glycine. Replacing the tyrosine with a glycine in a highly shikimate-specific angiosperm SDH was sufficient to gain some QDH function. Thus, very few mutations were necessary to facilitate the evolution of QDH genes.

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P

In plants, genes may sustain extensive pleiotropic functional properties by individually affecting multiple, distinct traits. We discuss results from three genome-wide association studies of approximately 400 natural poplar (Populus trichocarpa) accessions phenotyped for 60 ecological/biomass, wood quality, and rust fungus resistance traits. Single-nucleotide polymorphisms (SNPs) in the poplar ortholog of the class III homeodomain-leucine zipper transcription factor gene REVOLUTA (PtREV) were significantly associated with three specific traits. Based on SNP associations with fungal resistance, leaf drop, and cellulose content, the PtREV gene contains three potential regulatory sites within noncoding regions at the gene’s 3′ end, where alternative splicing and messenger RNA processing actively occur. The polymorphisms in this region associated with leaf abscission and cellulose content are suggested to represent more recent variants, whereas the SNP associated with leaf rust resistance may be more ancient, consistent with REV’s primary role in auxin signaling and its functional evolution in supporting fundamental processes of vascular plant development.

Insights into the physiology of nitrogen (N) uptake help us to understand the adaption of boreal coniferous forests to their environment. We compared fluxes of nitrate and ammonium in white spruce [Picea glauca (Moench) Voss] roots, measured using a non-invasive microelectrode ion flux measurement system (MIFE), and transcript abundance of ammonium and nitrate transporter genes in roots, determined by real time PCR. Seedlings were pretreated with water, or 50 μM or 1,500 μM NH4NO3 + 200 μM CaSO4 + 25 μM KH2PO4. Measurements were made on seedling roots 0–5, 5–10, 10–20 and 20–30 mm from the root tip. As ammonium and nitrate transporter family members in spruce are still uncharacterized, primers for real time PCR were designed to cover one family with each set of primers (AMT1, AMT2, NRT1, NRT2). The expression patterns obtained by real time PCR differed significantly among transporter family, treatments and root segments. Expression of AMT1 did not show a relationship with distance from the root tip, but the expression of AMT2 was generally greater 0–5 mm from the root tip than in segments farther from the tip. Expression of NRT1 was greatest 10–30 mm from the root tip, while expression of NRT2 was greatest 5–10 mm from the tip in all treatments, except the 1,500 μM NH4NO3 treatment. MIFE measurements showed the highest N uptake and proton efflux near the root tip and declining fluxes with increased distance from the root tip in the 50 μM N treatment. Significant net ammonium efflux was observed from some root segments in the 1,500 μM N treatment. Transporter gene expression and ion fluxes were not correlated. Though the measured net fluxes of ammonium were greater than the measured net fluxes of nitrate, the nitrate transporters were, in general, more highly expressed than the ammonium transporters.

The primary inflorescence stem of Arabidopsis thaliana is rich in lignified cell walls, in both vascular bundles and interfascicular fibres. Previous gene expression studies demonstrated a correlation between expression of phenylpropanoid biosynthetic genes and a subset of genes encoding ATP-binding cassette (ABC) transporters, especially in the ABCB/multi-drug resistance/P-glycoprotein (ABCB/MDR/PGP) and ABCG/pleiotropic drug resistance (ABCG/PDR) subfamilies. The objective of this study was to characterize these ABC transporters in terms of their gene expression and their function in development of lignified cells. Based on in silico analyses, four ABC transporters were selected for detailed investigation: ABCB11/MDR8, ABCB14/MDR12, ABCB15/MDR13, and ABCG33/PDR5. Promoter::glucuronidase reporter assays for each gene indicated that promoters of ABCB11, ABCB14, ABCB15, and ABCG33 transporters are active in the vascular tissues of primary stem, and in some cases in interfascicular tissues as well. Homozygous T-DNA insertion mutant lines showed no apparent irregular xylem phenotype or alterations in interfascicular fibre lignification or morphology in comparison with wild type. However, in abcb14-1 mutants, stem vascular morphology was slightly disorganized, with decreased phloem area in the vascular bundle and decreased xylem vessel lumen diameter. In addition, abcb14-1 mutants showed both decreased polar auxin transport through whole stems and altered auxin distribution in the procambium. It is proposed that both ABCB14 and ABCB15 promote auxin transport since inflorescence stems in both mutants showed a reduction in polar auxin transport, which was not observed for any of the ABCG subfamily mutants tested. In the case of ABCB14, the reduction in auxin transport is correlated with a mild disruption of vascular development in the inflorescence stem.

An approach based on an in silico analysis predicted that CYP77A4, a cytochrome P450 that so far has no identified function, might be a fatty acid-metabolizing enzyme. CYP77A4 was heterologously expressed in a Saccharomyces cerevisiae strain (WAT11) engineered for cytochrome P450 expression. Lauric acid (C12:0) was converted into a mixture of hydroxylauric acids when incubated with microsomes from yeast expressing CYP77A4. A variety of physiological C18 fatty acids were tested as potential substrates. Oleic acid (cis-Δ9C18:1) was converted into a mixture of ω-4- to ω-7-hydroxyoleic acids (75%) and 9,10-epoxystearic acid (25%). Linoleic acid (cis,cis-Δ9,Δ12C18:2) was exclusively converted into 12,13-epoxyoctadeca-9-enoic acid, which was then converted into diepoxide after epoxidation of the Δ9 unsaturation. Chiral analysis showed that 9,10-epoxystearic acid was a mixture of 9S/10R and 9R/10S in the ratio 33 : 77, whereas 12,13-epoxyoctadeca-9-enoic acid presented a strong enantiomeric excess in favor of 12S/13R, which represented 90% of the epoxide. Neither stearic acid (C18:0) nor linolelaidic acid (trans,trans-Δ9,Δ12C18:2) was metabolized, showing that CYP77A4 requires a double bond, in the cis configuration, to metabolize C18 fatty acids. CYP77A4 was also able to catalyze the in vitro formation of the three mono-epoxides of α-linolenic acid (cis,cis,cis-Δ9,Δ12,Δ15C18:3), previously described as antifungal compounds. Epoxides generated by CYP77A4 are further metabolized to the corresponding diols by epoxide hydrolases located in microsomal and cytosolic subcellular fractions from Arabidopsis thaliana. The concerted action of CYP77A4 with epoxide hydrolases and hydroxylases allows the production of compounds involved in plant–pathogen interactions, suggesting a possible role for CYP77A4 in plant defense.

Forest trees generally show high levels of local adaptation and efforts focusing on understanding adaptation to climate will be crucial for species survival and management. Here, we address fundamental questions regarding the molecular basis of adaptation in undomesticated forest tree populations to past climatic environments by employing an integrative quantitative genetics and landscape genomics approach. Using this comprehensive approach, we studied the molecular basis of climate adaptation in 433 Populus trichocarpa (black cottonwood) genotypes originating across western North America. Variation in 74 field-assessed traits (growth, ecophysiology, phenology, leaf stomata, wood, and disease resistance) was investigated for signatures of selection (comparing Q ST -F ST) using clustering of individuals by climate of origin (temperature and precipitation). 29,354 SNPs were investigated employing three different outlier detection methods and marker-inferred relatedness was estimated to obtain the narrow-sense estimate of population differentiation in wild populations. In addition, we compared our results with previously assessed selection of candidate SNPs using the 25 topographical units (drainages) across the P. trichocarpa sampling range as population groupings. Narrow-sense Q ST for 53% of distinct field traits was significantly divergent from expectations of neutrality (indicating adaptive trait variation); 2,855 SNPs showed signals of diversifying selection and of these, 118 SNPs (within 81 genes) were associated with adaptive traits (based on significant Q ST). Many SNPs were putatively pleiotropic for functionally uncorrelated adaptive traits, such as autumn phenology, height, and disease resistance. Evolutionary quantitative genomics in P. trichocarpa provides an enhanced understanding regarding the molecular basis of climate-driven selection in forest trees and we highlight that important loci underlying adaptive trait variation also show relationship to climate of origin. We consider our approach the most comprehensive, as it uncovers the molecular mechanisms of adaptation using multiple methods and tests. We also provide a detailed outline of the required analyses for studying adaptation to the environment in a population genomics context to better understand the species’ potential adaptive capacity to future climatic scenarios.

Three independent cytochrome P450 enzyme families catalyze the three rate-limiting hydroxylation steps in the phenylpropanoid pathway leading to the biosynthesis of lignin and numerous other phenolic compounds in plants. Their characterization at the molecular and enzymatic level has revealed an unexpected complexity of phenolic metabolism as the major route involves shikimate/quinate esters and alcohol/aldehyde intermediates. Engineering expression of CYP73s (encoding cinnamate 4-hydroxylase), CYP98s (encoding 4-coumaroylshikimate 3′-hydroxylase) or CYP84s (encoding coniferaldehyde 5-hydroxylase) leads to modified lignin and seed phenolic composition. In particular CYP73s and CYP98s also play essential roles in plant growth and development, while CYP84 constitutes a check-point for the synthesis of syringyl lignin and sinapate esters. Although recent data shed new light on the main path for lignin synthesis, they also raised new questions. Mutants and engineered plants revealed the existence of (an) alternative pathway(s), which most likely involve(s) different precursors and oxygenases. On the other hand, phylogenetic analysis of plant genomes show the existence of P450 gene duplications in each family, which may have led to the acquisition of novel or additional physiological functions in planta. In addition to the main lignin pathway, P450s contribute to the biosynthesis of many bioactive phenolic derivatives, with potential applications in medicine and plant defense, including lignans, phenylethanoids, benzoic acids, xanthones or quinoid compounds. A very small proportion of these P450s have been characterized so far, and rarely at a molecular level. The possible involvement of P450s in salicylic acid is discussed.

Cytochrome P450 mono-oxygenases play prominent roles in a diverse set of metabolic pathways, but the function of most of these enzymes remains obscure. A bottleneck in the functional genomics of this superfamily constitutes hypothesis generation to identify potential substrates (or substrate classes) individual P450s may act on. We used publicly available large-scale expression data to perform co-expression analysis comparing the expression matrix of each P450 with those from more than 4000 selected genes across thousands of microarrays. Based on functional annotations of co-expressed genes from a diverse set of databases, co-expressed pathways were thus identified for each P450. Using this approach, most P450s with known functions were placed into their respective pathways, thereby proofing the concept. As examples, pathway mapping results identifying novel P450s potentially acting on flower-specific monoterpenes and root-specific triterpenes are described. Co-expression results for all Arabidopsis P450s will be presented as a web resource on the ‘CYPedia’ web pages (http://ibmp.u-strasbg.fr/~CYPedia/).

Sequenced plant genomes provide a large reservoir of known genes with potential for use in crop and tree improvement, but assignment of specific functions to annotated genes in sequenced plant genomes remains a challenge. Furthermore, most plant genes belong to families encoding proteins with related but distinct functions. In this commentary, we discuss our development of Arabidopsis spotted whole genome longmer oligonucleotide microarrays, and their use in global transcription profiling. We show that longmer array based transcriptome analysis in Arabidopsis can be used as an efficient and effective gene discovery and functional genomics tool, particularly for functional analyses of members of large gene families. We discuss experiments that focus on gene families involved in phenylpropanoid natural product biosynthesis and fiber differentiation. These analyses have helped to elucidate functions of individual gene family members, and have identified new candidate genes involved in fiber development and differentiation. Results obtained by these studies in Arabidopsis can be used as the basis for gene discovery in commercially important plants, and we have focused our attention on Populus trichocarpa (poplar), a species important in forestry and agroforestry for which complete genome sequence information is available.

Different stages of vascular and interfascicular fiber differentiation can be identified along the axis of bolting stems in Arabidopsis. To gain insights into the metabolic, developmental, and regulatory events that control this pattern, we applied global transcript profiling employing an Arabidopsis full-genome longmer microarray. More than 5000 genes were differentially expressed, among which more than 3000 changed more than twofold, and were placed into eight expression clusters based on polynomial regression models. Within these, 182 upregulated transcription factors represent candidate regulators of fiber development. A subset of these candidates has been associated with fiber development and/or secondary wall formation and lignification in the literature, making them targets for functional studies and comparative genomic analyses with woody plants. Analysis of differentially expressed phenylpropanoid genes identified a set known to be involved in lignin biosynthesis. These were used to anchor co-expression analyses that allowed us to identify candidate genes encoding proteins involved in monolignol transport and monolignol dehydrogenation and polymerization. Similar analyses revealed candidate genes encoding enzymes that catalyze missing links in the shikimate pathway, namely arogenate dehydrogenase and prephenate aminotransferase.

NADPH:cytochrome P450 reductase (CPR) provides reducing equivalents to diverse cytochrome P450 monooxygenases. We isolated cDNAs for three CPR genes (CPR1,CPR2, and CPR3) from hybrid poplar (Populus trichocarpa × Populus deltoides). Deduced CPR2 and CPR3 amino acid sequences were 91% identical, but encoded isoforms divergent from CPR1 (72% identity). CPR1 and CPR2 were co-expressed together with the P450 enzyme cinnamate-4-hydroxylase (C4H) in yeast (Saccharomyces cerevisiae). Microsomes isolated from strains expressing CPR1/C4H or CPR2/C4H enhanced C4H activities approximately 10-fold relative to the C4H-only control strain, and catalyzed NADPH-dependent cytochrome c reduction. The divergent CPR isoforms (CPR1 and CPR2/3) contained entirely different N-terminal sequences, which are conserved in other plant CPRs and are diagnostic for two distinct classes of CPRs within the angiosperms. C-terminal green fluorescent protein fusions to CPR1 and CPR2 were constructed and expressed in both yeast and Arabidopsis. The fusion proteins expressed in yeast retained the ability to support C4H activity and, thus, were catalytically active. Both CPR::green fluorescent protein fusion proteins were strictly localized to the endoplasmic reticulum in transgenic Arabidopsis. The lack of localization of either isoform to chloroplasts, where P450s are known to be present, suggests that an alternative P450 reduction system may be operative in this organelle. Transcripts for the three poplar CPR genes were present ubiquitously in all tissues examined, but CPR2 showed highest expression in young leaf tissue.

Contrasting with the wealth of information available on the multiple roles of jasmonates in plant development and defense, knowledge about the functions and the biosynthesis of hydroxylated oxylipins remains scarce. By expressing the caleosin RESPONSIVE TO DESSICATION20 (RD20) in Saccharomyces cerevisiae, we show that the recombinant protein possesses an unusual peroxygenase activity with restricted specificity toward hydroperoxides of unsaturated fatty acid. Accordingly, Arabidopsis (Arabidopsis thaliana) plants overexpressing RD20 accumulate the product 13-hydroxy-9,11,15-octadecatrienoic acid, a linolenate-derived hydroxide. These plants exhibit elevated levels of reactive oxygen species (ROS) associated with early gibberellin-dependent flowering and abscisic acid hypersensitivity at seed germination. These phenotypes are dependent on the presence of active RD20, since they are abolished in the rd20 null mutant and in lines overexpressing RD20, in which peroxygenase was inactivated by a point mutation of a catalytic histidine residue. RD20 also confers tolerance against stress induced by Paraquat, Rose Bengal, heavy metal, and the synthetic auxins 1-naphthaleneacetic acid and 2,4-dichlorophenoxyacetic acid. Under oxidative stress, 13-hydroxy-9,11,15-octadecatrienoic acid still accumulates in RD20-overexpressing lines, but this lipid oxidation is associated with reduced ROS levels, minor cell death, and delayed floral transition. A model is discussed where the interplay between fatty acid hydroxides generated by RD20 and ROS is counteracted by ethylene during development in unstressed environments.

The shikimate pathway leads to the biosynthesis of aromatic amino acids essential for protein biosynthesis and the production of a wide array of plant secondary metabolites. Among them, quinate is an astringent feeding deterrent that can be formed in a single step reaction from 3-dehydroquinate catalyzed by quinate dehydrogenase (QDH). 3-Dehydroquinate is also the substrate for shikimate biosynthesis through the sequential actions of dehydroquinate dehydratase (DQD) and shikimate dehydrogenase (SDH) contained in a single protein in plants. The reaction mechanism of QDH resembles that of SDH. The poplar genome encodes five DQD/SDH-like genes (Poptr1 to Poptr5), which have diverged into two distinct groups based on sequence analysis and protein structure prediction. In vitro biochemical assays proved that Poptr1 and -5 are true DQD/SDHs, whereas Poptr2 and -3 instead have QDH activity with only residual DQD/SDH activity. Poplar DQD/SDHs have distinct expression profiles suggesting separate roles in protein and lignin biosynthesis. Also, the QDH genes are differentially expressed. In summary, quinate (secondary metabolism) and shikimate (primary metabolism) metabolic activities are encoded by distinct members of the same gene family, each having different physiological functions.

4-coumarate:CoA ligase (4CL), the last enzyme of the general phenylpropanoid pathway, provides precursors for the biosynthesis of a large variety of plant natural products. 4 CL catalyzes the formation of CoA thiol esters of 4-coumarate and other hydroxycinnamates in a two step reaction involving the formation of an adenylate intermediate. 4 CL shares conserved peptide motifs with diverse adenylate-forming enzymes such as firefly luciferases, non-ribosomal peptide synthetases, and acyl:CoA synthetases. Amino acid residues involved in 4 CL catalytic activities have been identified, but domains involved in determining substrate specificity remain unknown. To address this question, we took advantage of the difference in substrate usage between the Arabidopsis thaliana 4 CL isoforms At4CL1 and At4CL2. While both enzymes convert 4-coumarate, only At4CL1 is also capable of converting ferulate. Employing a domain swapping approach, we identified two adjacent domains involved in substrate recognition. Both substrate binding domain I (sbd I) and sbd II of At4CL1 alone were sufficient to confer ferulate utilization ability upon chimeric proteins otherwise consisting of At4CL2 sequences. In contrast, sbd I and sbd II of At4CL2 together were required to abolish ferulate utilization in the context of At4CL1. Sbd I corresponds to a region previously identified as the substrate binding domain of the adenylation subunit of bacterial peptide synthetases, while sbd II centers on a conserved domain of so far unknown function in adenylate-forming enzymes (GEI/LxIxG). At4CL1 and At4CL2 differ in nine amino acids within sbd I and four within sbd II, suggesting that these play roles in substrate recognition.

The cytochrome P450 family encompasses the largest family of enzymes in plant metabolism, and the functions of many of its members in Arabidopsis thaliana are still unknown. Gene coexpression analysis pointed to two P450s that were coexpressed with two monoterpene synthases in flowers and were thus predicted to be involved in monoterpenoid metabolism. We show that all four selected genes, the two terpene synthases (TPS10 and TPS14) and the two cytochrome P450s (CYP71B31 and CYP76C3), are simultaneously expressed at anthesis, mainly in upper anther filaments and in petals. Upon transient expression in Nicotiana benthamiana, the TPS enzymes colocalize in vesicular structures associated with the plastid surface, whereas the P450 proteins were detected in the endoplasmic reticulum. Whether they were expressed in Saccharomyces cerevisiae or in N. benthamiana, the TPS enzymes formed two different enantiomers of linalool: (−)-(R)-linalool for TPS10 and (+)-(S)-linalool for TPS14. Both P450 enzymes metabolize the two linalool enantiomers to form different but overlapping sets of hydroxylated or epoxidized products. These oxygenated products are not emitted into the floral headspace, but accumulate in floral tissues as further converted or conjugated metabolites. This work reveals complex linalool metabolism in Arabidopsis flowers, the ecological role of which remains to be determined.

Summary Actin is the best-studied member of a family of ATPases that originated, presumably via gene duplications, in the common ancestor of all living things.The discovery of actin-related proteins (ARPs) and the resolution of the three-dimensional structure of actin has refueled interest in the phylogeny of this coding region.In comparison to vertebrates and land plants, relatively little is known, however, about actin and actin gene families in protists.In this paper, the origin of actin is reviewed and sequence analyses of this coding region are presented.Phylogenetic analyses of actins are used to probe the evolutionary relationships between actin isoforms in three lineages, the prymnesiophytes (Le.Emiliana huxleyi), the heterokonts and the green~,algae/land plants, and the relationships between these sequences with actins from diverse eukaryotes.

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost-effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre-ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

The phenylpropanoid metabolism channels carbon from phenylalanine to the three monolignols and numerous other phenolic compounds. Our understanding of the pathway has changed tremendously over the last decade, which was driven largely by the biochemical and genetic characterization of the cytochrome P450s catalysing the hydroxylation of the aromatic ring. The first, cinnamate 4-hydroxylase, is the rate-limiting step into the phenylpropanoid pathway and is highly specific for cinnamate. Blocking this step impairs the ability of plants to produce lignin. The 3-hydroxylation occurs primarily on the level of the 4-coumaroyl-shikimate level, rather than on the free acid or CoA-ester as previously expected, thereby linking the far upstream shikimate pathway with the committed step towards S and G lignin. Finally, the last hydroxylation step occurs on the level of the aldehyde or alcohol and defines flow into S lignin as indicated by mutant and over-expressing lines. Here, we summarize our current understanding of the phenylpropanoid pathway, describe the phenotypic consequences of miss-regulating the rate limiting cytochrome P450s involved in the pathway on lignin structure and quantity, and discuss the flexibility of plants in channelling carbon through the phenylpropanoid grid.

An approach based on an in silico analysis predicted that CYP77A4, a cytochrome P450 that so far has no identified function, might be a fatty acid-metabolizing enzyme. CYP77A4 was heterologously expressed in a Saccharomyces cerevisiae strain (WAT11) engineered for cytochrome P450 expression. Lauric acid (C(12:0)) was converted into a mixture of hydroxylauric acids when incubated with microsomes from yeast expressing CYP77A4. A variety of physiological C(18) fatty acids were tested as potential substrates. Oleic acid (cis-Delta(9)C(18:1)) was converted into a mixture of omega-4- to omega-7-hydroxyoleic acids (75%) and 9,10-epoxystearic acid (25%). Linoleic acid (cis,cis-Delta(9),Delta(12)C(18:2)) was exclusively converted into 12,13-epoxyoctadeca-9-enoic acid, which was then converted into diepoxide after epoxidation of the Delta(9) unsaturation. Chiral analysis showed that 9,10-epoxystearic acid was a mixture of 9S/10R and 9R/10S in the ratio 33 : 77, whereas 12,13-epoxyoctadeca-9-enoic acid presented a strong enantiomeric excess in favor of 12S/13R, which represented 90% of the epoxide. Neither stearic acid (C(18:0)) nor linolelaidic acid (trans,trans-Delta(9),Delta(12)C(18:2)) was metabolized, showing that CYP77A4 requires a double bond, in the cis configuration, to metabolize C(18) fatty acids. CYP77A4 was also able to catalyze the in vitro formation of the three mono-epoxides of alpha-linolenic acid (cis,cis,cis-Delta(9),Delta(12),Delta(15)C(18:3)), previously described as antifungal compounds. Epoxides generated by CYP77A4 are further metabolized to the corresponding diols by epoxide hydrolases located in microsomal and cytosolic subcellular fractions from Arabidopsis thaliana. The concerted action of CYP77A4 with epoxide hydrolases and hydroxylases allows the production of compounds involved in plant-pathogen interactions, suggesting a possible role for CYP77A4 in plant defense.

From Retrogene to Phenolic Metabolism Metabolic plasticity, which involves the creation of new genes, is an essential feature of plant adaptation and speciation. Studying plants from the mustard family, Matsuno et al. (p. 1688 ) show that variants of the cytochrome P450 enzyme family were derived through retroposition, duplication, and subsequent mutaton. Evolutionary changes increased the volume of the substrate pocket altering with what sorts of substrates the enzymes could interact. The enzymes formed the basis for a new metabolic pathway, the products of which include constituents of pollen and of phenylpropanoid metabolism.

Within the last decade, the strategy for pathway discovery in plant secondary metabolism has reversed from a metabolite to protein to gene approach to a gene to pathway strategy. Different genome-wide “omics” strategies apply to gene discovery in model plants with sequenced genomes and plants of pharmaceutical or industrial interest. In this chapter, we provide a brief description as well as a few examples of the main approaches that have so far been applied to plant metabolism. Combination of such global approaches leads to the new field of ‘integrative biology’, which highlights metabolic networks connecting the different branches of primary and “secondary” metabolism.

Background Plants 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. Results Using 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. Conclusion Arabidopsis 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.

Background Sequencing of the first plant genomes has revealed that cytochromes P450 have evolved to become the largest family of enzymes in secondary metabolism. The proportion of P450 enzymes with characterized biochemical function(s) is however very small. If P450 diversification mirrors evolution of chemical diversity, this points to an unexpectedly poor understanding of plant metabolism. We assumed that extensive analysis of gene expression might guide towards the function of P450 enzymes, and highlight overlooked aspects of plant metabolism. Results We have created a comprehensive database, 'CYPedia', describing P450 gene expression in four data sets: organs and tissues, stress response, hormone response, and mutants of Arabidopsis thaliana, based on public Affymetrix ATH1 microarray expression data. P450 expression was then combined with the expression of 4,130 re-annotated genes, predicted to act in plant metabolism, for co-expression analyses. Based on the annotation of co-expressed genes from diverse pathway annotation databases, co-expressed pathways were identified. Predictions were validated for most P450s with known functions. As examples, co-expression results for P450s related to plastidial functions/photosynthesis, and to phenylpropanoid, triterpenoid and jasmonate metabolism are highlighted here. Conclusion The large scale hypothesis generation tools presented here provide leads to new pathways, unexpected functions, and regulatory networks for many P450s in plant metabolism. These can now be exploited by the community to validate the proposed functions experimentally using reverse genetics, biochemistry, and metabolic profiling.

4-coumarate::CoA ligase (4CL) gene family members are involved in channeling carbon flow into branch pathways of phenylpropanoid metabolism. Transgenic Arabidopsis plants containing the At4CL1 or At4CL2 promoter fused to the beta-glucuronidase (GUS) reporter gene show developmentally regulated GUS expression in the xylem tissues of the root and shoot. To identify regulatory genes involved in the developmental regulation of At4CL and other phenylpropanoid-specific genes, we generated ethyl methyl sulfate mutagenized populations of At4CL1::GUS and At4CL2::GUS transgenic lines and screened approximately 16,000 progeny for reduced or altered GUS expression. Several lines with reproducible patterns of reduced GUS expression were identified. However, the GUS-expression phenotype segregated in a non-Mendelian manner in all of the identified lines. Also, GUS expression was restored by 5-azacytidine (aza) treatment, suggesting inhibitory DNA methylation of the transgene. Southern analysis confirmed DNA methylation of the proximal promoter sequences of the transgene only in the mutant lines. In addition, retransformation of At4CL::GUS lines with further At4CL promoter constructs enhanced the GUS-silencing phenotype. Taken together, these results suggest that the isolated mutants are epimutants. Apparently, two different modes of silencing were engaged in the At4CL1::GUS and At4CL2::GUS silenced lines. While silencing in the seedlings of the At4CL1::GUS lines was root specific in seedlings, it affected all organs in the At4CL2::GUS lines. Also, At4CL1::GUS transgene silencing was confined to the transgene but At4CL2::GUS silencing extended to the endogenous At4CL2 gene. Organ-specific silencing of the At4CL1::GUS transgene cannot be explained by current models in the literature.

This chapter reconstructs the whole shikimate biosynthetic pathway, the central pathway of plant metabolism controlling metabolic flux from carbohydrate metabolism into the phenylpropanoid pathway, on the genomic level. Since the chapter focuses on fully sequenced genomes, all families of the three species rice, Arabidopsis , and poplar should be represented in their full extent. The chapter also includes expression data for Arabidopsis and poplar that support functional classifications, providing an unprecedented view of the structure and expression of the gene families that encode shikimate pathway enzymes. This comparative data also allows to build on previous observations to highlight promising candidate genes that encode yet unknown enzymes for missing links in shikimate pathway required for phenylalanine biosynthesis.

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 Arabidopsis thaliana genome contains at least 32 terpenoid synthase ( AtTPS ) genes [Aubourg et al., Mol. Genet. Genom. 267 (2002) 730] a few of which have recently been characterized. Based on hierarchical cluster analysis of AtTPS gene expression, measured by microarray profiling and validated with published expression data, we identified two groups of predominantly root expressed AtTPS genes containing five members with previously unknown biochemical functions ( At4g13280 , At4g13300 , At5g48110 , At1g33750 , and At3g29410 ). Among the root expressed AtTPS genes, a pair of tandem-organized genes, At4g13280 (AtTPS12) and At4g13300 (AtTPS13), shares 91% predicted amino acid identity indicating recent gene duplication. Bacterial expression of cDNAs and enzyme assays showed that both At4g13280 and At4g13300 encode sesquiterpene synthases catalyzing the conversion of farnesyl diphosphate to ( Z )-γ-bisabolene and the additional minor products E -nerolidol and α-bisabolol. Expression of β-glucuronidase ( GUS ) reporter gene fused to upstream genomic regions of At4g13280 or At4g13300 showed constitutive promoter activities in the cortex and sub-epidermal layers of Arabidopsis roots. In addition, highly localized promoter activities were found in leaf hydathodes and flower stigmata. Mechanical wounding of Arabidopsis leaves induced local expression of At4g13280 and At4g13300 . The functional characterization of At4g13280 gene product AtTPS12 and At4g13230 gene product AtTPS13 as ( Z )-γ-bisabolene synthases, together with the recent characterization of two flower-specific AtTPS [ At5g23960 and At5g44630 ; Tholl et al., Plant J. 42 (2005) 757], concludes the biochemical functional annotation of all four predicted Arabidopsis sesquiterpene synthase genes. Our data suggest biological functions for At4g13280 and At4g13300 in the rhizosphere with additional roles in aerial plant tissues.

Different stages of vascular and interfascicular fiber differentiation can be identified along the axis of bolting stems in Arabidopsis. To gain insights into the metabolic, developmental, and regulatory events that control this pattern, we applied global transcript profiling employing an Arabidopsis full-genome longmer microarray. More than 5000 genes were differentially expressed, among which more than 3000 changed more than twofold, and were placed into eight expression clusters based on polynomial regression models. Within these, 182 upregulated transcription factors represent candidate regulators of fiber development. A subset of these candidates has been associated with fiber development and/or secondary wall formation and lignification in the literature, making them targets for functional studies and comparative genomic analyses with woody plants. Analysis of differentially expressed phenylpropanoid genes identified a set known to be involved in lignin biosynthesis. These were used to anchor co-expression analyses that allowed us to identify candidate genes encoding proteins involved in monolignol transport and monolignol dehydrogenation and polymerization. Similar analyses revealed candidate genes encoding enzymes that catalyze missing links in the shikimate pathway, namely arogenate dehydrogenase and prephenate aminotransferase.

The phenylpropanoid enzyme 4-coumarate:coenzyme A ligase (4CL) plays a key role in general phenylpropanoid metabolism. 4CL is related to a larger class of prokaryotic and eukaryotic adenylate-forming enzymes and shares several conserved peptide motifs with these enzymes. In order to better characterize the nature of 4CL gene families in poplar, parsley, and tobacco, we used degenerate primers to amplify 4CL sequences from these species. In each species additional, divergent 4CL genes were found. Complete cDNA clones for the two new poplar 4CL genes were obtained, allowing examination of their expression patterns and determination of the substrate utilization profile of a xylem-specific isoform. Phylogenetic analysis of these genes and gene fragments confirmed previous results showing that 4CL proteins fall into two evolutionarily ancient subgroups . A comparative phylogenetic analysis of enzymes in the adenylate-forming superfamily showed that 4CLs, luciferases, and acetate CoA ligases each form distinct clades within the superfamily. According to this analysis, four Arabidopsis 4CL-like genes identified from the Arabidopsis Genome Project are only distantly related to bona fide 4CLs or are more closely related to fatty acid CoA ligases, suggesting that the three Arabidopsis 4CL genes previously characterized represent the extent of the 4CL gene family in this species.

4-Coumarate:coenzyme A ligase (4CL) is a key enzyme of general phenylpropanoid metabolism which provides the precursors for a large variety of important plant secondary products, such as lignin, flavonoids, or phytoalexins. To identify amino acids important for 4CL activity, eight mutations were introduced into Arabidopsis thaliana At4CL2. Determination of specific activities and K m values for ATP and caffeate of the heterologously expressed and purified proteins identified four distinct classes of mutants: enzymes with little or no catalytic activity; enzymes with greatly reduced activity but wild-type K m values; enzymes with drastically altered K m values; and enzymes with almost wild-type properties. The latter class includes replacement of a cysteine residue which is strictly conserved in 4CLs and had previously been assumed to be directly involved in catalysis. These results substantiate the close relationship between 4CL and other adenylate-forming enzymes such as luciferases, peptide synthetases, and fatty acyl-CoA synthetases.

Summary The enzyme 4-coumarate:CoA ligase (4CL) plays a key role in channelling carbon flow into diverse branch pathways of phenylpropanoid metabolism which serve important functions in plant growth and adaptation to environmental perturbations. Here we report on the cloning of the 4CL gene family from Arabidopsis thaliana and demonstrate that its three members, At4CL1, At4CL2 and At4CL3, encode isozymes with distinct substrate preference and specificities. Expression studies revealed a differential behaviour of the three genes in various plant organs and upon external stimuli such as wounding and UV irradiation or upon challenge with the fungus, Peronospora parasitica. Phylogenetic comparisons indicate that, in angiosperms, 4CL can be classified into two major clusters, class I and class II, with the At4CL1 and At4CL2 isoforms belonging to class I and At4CL3 to class II. Based on their enzymatic properties, expression characteristics and evolutionary relationships, At4CL3 is likely to participate in the biosynthetic pathway leading to flavonoids whereas At4CL1 and At4CL2 are probably involved in lignin formation and in the production of additional phenolic compounds other than flavonoids.

Lignin is an important biopolymer that is deposited in secondary cell walls of plant cells (e.g., tracheary elements) and in response to stresses such as wounding. Biosynthesis of lignin monomers occurs via the phenylpropanoid pathway, in which the enzyme 4-coumarate:CoA ligase (4CL) plays a key role by catalyzing the formation of hydroxycinnamoyl-CoA esters, subsequently reduced to the corresponding monolignols (hydroxycinnamoyl alcohols). 4CL is encoded by a family of four genes in Arabidopsis thaliana (At4CL1-At4CL4), which are developmentally regulated and co-expressed with other phenylpropanoid genes. We investigated in detail the wound-induced expression of At4CL1-At4CL4, and found that At4CL1 and At4CL2 mRNA accumulation follows biphasic kinetics over a period of 72 h, while At4CL4 expression is rapidly activated for a period of at least 12 h before declining. In order to localize cis-regulatory elements involved in the developmental and wound-induced regulation of the At4CL gene family members, At4CL promoter-beta-glucuronidase (GUS) reporter gene fusions were constructed and transferred into Arabidopsis plants. Analysis of these plants revealed that the promoter fragments direct discrete and distinct patterns of expression, some of which did not recapitulate expected patterns of wound-induced expression. The locations of regulatory elements associated with the At4CL2 gene were investigated in detail using a series of transgenic Arabidopsis plants containing promoter fragments and parts of the transcribed region of the gene fused to GUS. Positive and negative regulatory elements effective in modulating developmental expression or wound responsiveness of the gene were located both in the promoter and transcribed regions of the At4CL2 gene.

Background Nitrogen uptake, reallocation within the plant, and between subcellular compartments involves ammonium, nitrate and peptide transporters. Ammonium transporters are separated into two distinct families (AMT1 and AMT2), each comprised of five members on average in angiosperms. Nitrate transporters also form two discrete families (NRT1 and NRT2), with angiosperms having four NRT2s, on average. NRT1s share an evolutionary history with peptide transporters (PTRs). The NRT1/PTR family in land plants usually has more than 50 members and contains also members with distinct activities, such as glucosinolate and abscisic acid transport. Results Phylogenetic reconstructions of each family across 20 land plant species with available genome sequences were supplemented with subcellular localization and transmembrane topology predictions. This revealed that both AMT families diverged prior to the separation of bryophytes and vascular plants forming two distinct clans, designated as supergroups, each. Ten supergroups were identified for the NRT1/PTR family. It is apparent that nitrate and peptide transport within the NRT1/PTR family is polyphyletic, that is, nitrate and/or peptide transport likely evolved multiple times within land plants. The NRT2 family separated into two distinct clans early in vascular plant evolution. Subsequent duplications occurring prior to the eudicot/monocot separation led to the existence of two AMT1, six AMT2, 31 NRT1/PTR, and two NRT2 clans, designated as groups. Conclusion Phylogenetic separation of groups suggests functional divergence within the angiosperms for each family. Distinct groups within the NRT1/PTR family appear to separate peptide and nitrate transport activities as well as other activities contained within the family, for example nitrite transport. Conversely, distinct activities, such as abscisic acid and glucosinolate transport, appear to have recently evolved from nitrate transporters.

Plants produce a rich variety of natural products to face environmental constraints. Enzymes of the cytochrome P450 CYP98 family are key actors in the production of phenolic bioactive compounds. They hydroxylate phenolic esters for lignin biosynthesis in angiosperms, but also produce various other bioactive phenolics. We characterized CYP98s from a moss, a lycopod, a fern, a conifer, a basal angiosperm, a monocot and from two eudicots. We found that substrate preference of the enzymes has changed during evolution of land plants with typical lignin-related activities only appearing in angiosperms, suggesting that ferns, similar to lycopods, produce lignin through an alternative route. A moss CYP98 knock-out mutant revealed coumaroyl-threonate as CYP98 substrate in vivo and showed a severe phenotype. Multiple CYP98s per species exist only in the angiosperms, where we generally found one isoform presumably involved in the biosynthesis of monolignols, and additional isoforms, resulting from independent duplications, with a broad range of functions in vitro.; Les plantes produisent une grande variété de produits naturels pour faire face aux conditions environnementales. Les enzymes de la famille CYP98 des cytochromes P450 sont des enzymes clés dans la production des composés dérivés de la voie des phénylpropanoïdes. Ces enzymes sont impliquées dans l'hydroxylation des esters phénoliques pour la biosynthèse des monolignols chez les angiospermes, mais elles sont également impliquées dans la production de divers autres composés phénoliques solubles. Nous avons caractérisé des CYP98 représentatifs des mousses, Lycopodes, fougères, Gymnospermes, Angiospermes basales, Monocotylédones et Eudicotylédones et démontré que leur préférence de substrat a changé au cours de l'évolution. Un mutant knock-out de CYP98 de mousse a révélé un phénotype sévère et que le p-coumaroyl-thréonate est substrat de l’enzyme in vivo. Une duplication des CYP98s ne peut être observée que dans le génome des Angiospermes, qui présentent généralement une isoforme potentiellement impliquée dans la biosynthèse de la lignine et autres isoformes, résultant de duplications indépendantes, dont le spectre de substrats est plus large in vitro.