The CYP71AZ P450 Subfamily: A Driving Factor for the Diversification of Coumarin Biosynthesis in Apiaceous Plants

Three types of enzymes complete the furanocoumarins core skeleton biosynthesis in Angelica sinensis

Furanocoumarins (FCs) are widely distributed secondary metabolites found in higher plants, including Apiaceae, Rutaceae, Moraceae, and Fabaceae. They play a crucial role in the physiological functions of plants and are well-known for their diverse pharmacological activities. As a representative plant of the Apiaceae family, Angelica sinensis is highly valued for its medicinal properties and FCs are one of the main ingredients of A. sinensis. However, the biosynthetic mechanism of FCs in A. sinensis remains poorly understood. In this study, we successfully cloned and verified three types of enzymes using genome analysis and in vitro functional verification, which complete the biosynthesis of the FCs core skeleton in A. sinensis. It includes a p-coumaroyl CoA 2'-hydroxylase (AsC2'H) responsible for umbelliferone formation, two UbiA prenyltransferases (AsPT1 and AsPT2) that convert umbelliferone to demethylsuberosin (DMS) and osthenol, respectively, and two CYP736 subfamily cyclases (AsDC and AsOD) that catalyze the formation of FCs core skeleton. Interestingly, AsOD was demonstrated to be a bifunctional cyclase and could catalyze both DMS and osthenol, but had a higher affinity to osthenol. The characterization of these enzymes elucidates the molecular mechanism of FCs biosynthesis, providing new insights and technologies for understanding the diverse origins of FCs biosynthesis.

Two types of coumarins-specific enzymes complete the last missing steps in pyran- and furanocoumarins biosynthesis

Pyran- and furanocoumarins are key representatives of tetrahydropyrans and tetrahydrofurans, respectively, exhibiting diverse physiological and medical bioactivities. However, the biosynthetic mechanisms for their core structures remain poorly understood. Here we combined multiomics analyses of biosynthetic enzymes in Peucedanum praeruptorum and in vitro functional verification and identified two types of key enzymes critical for pyran and furan ring biosynthesis in plants. These included three distinct P. praeruptorum prenyltransferases (PpPT1-3) responsible for the prenylation of the simple coumarin skeleton 7 into linear or angular precursors, and two novel CYP450 cyclases (PpDC and PpOC) crucial for the cyclization of the linear/angular precursors into either tetrahydropyran or tetrahydrofuran scaffolds. Biochemical analyses of cyclases indicated that acid/base-assisted epoxide ring opening contributed to the enzyme-catalyzed tetrahydropyran and tetrahydrofuran ring refactoring. The possible acid/base-assisted catalytic mechanisms of the identified cyclases were theoretically investigated and assessed using site-specific mutagenesis. We identified two possible acidic amino acids Glu303 in PpDC and Asp301 in PpOC as vital in the catalytic process. This study provides new enzymatic tools in the epoxide formation/epoxide-opening mediated cascade reaction and exemplifies how plants become chemically diverse in terms of enzyme function and catalytic process.

Transcriptome analysis provides insights into coumarin biosynthesis in the medicinal plant Angelica dahurica cv. Yubaizhi

Angelica dahurica roots have a long history of use in traditional Chinese medicine due to their high coumarin content. To address the increasing demand for these roots, a synthetic biology approach has been proposed. Nevertheless, our comprehension of coumarin biosynthesis and its regulation remains limited. In this study, we utilized Hiseq2500 sequencing to analyze the transcriptomes of A. dahurica at different growth stages while concurrently quantifying coumarin content. Differentially expressed gene (DEG) analysis was employed to identify key genes involved in coumarin and terpenoid backbone biosynthesis. Weighted gene co-expression network analysis (WGCNA) was applied to identify gene modules strongly associated with coumarin content, elucidating the regulatory relationships between transcription factors (TFs) and pathway genes. Furthermore, KEGG enrichment analysis was used to explore essential pathways governing coumarin biosynthesis, with the identification of hub genes. Our results indicated that total coumarin content was highest in the roots, followed by leaves and stems, across all three developmental stages. Transcriptome analysis identified a total of 92,478 genes, among which 215 and 30 genes were implicated in coumarin and terpenoid backbone biosynthesis, respectively. Within the 73 identified gene modules by WGCNA, three modules-namely aquamarine1 (comprising two OMTs, one CSE, one AACT, one HDS, two PSs, one 2OGO, four UGTs, and seven CYP450s), darkmagenta (containing one UGT and one HDR), and navajowhite2 (consisting of one HCT, three UGTs, one CYP71A25, one OMT, one CSE, one HDS, and one PT)-were strongly associated with imperatorin, oxypeucedanin, and isoimperatorin content, respectively. KEGG enrichment analysis highlighted significant enrichment of cytochrome P450, transporter, and ubiquitin system pathways. Moreover, TF-gene regulatory analysis unveiled the complexity of coumarin biosynthesis, with 17 TF families regulating 17 genes in the aquamarine1 module, 8 TF families regulating 2 genes in the darkmagenta module, and 8 TF families regulating 7 genes in the navajowhite2 module. These comprehensive findings provide valuable insights into coumarin biosynthesis in A. dahurica, facilitating future research and potential applications in traditional Chinese medicine and synthetic biology strategies.

Reconstitution and Optimization of the Marmesin Biosynthetic Pathway in Yeast

Marmesin is essential in plant defense systems and exhibits various biological activities. In this study, we reconstituted the marmesin biosynthetic pathway in the Saccharomyces cerevisiae BY4741 chassis. We engineered the aromatic amino acid (AAA) biosynthetic pathways by introducing Escherichia coli-derived ppsA to improve the availability of the AAA precursor phosphoenolpyruvate, overexpressing the feedback inhibition resistance genes ARO4K229L and ARO7G141S to direct the metabolic flux toward the tyrosine branch, and deleting ARO10, PDC5, and PDC6 to reduce the byproducts from the Ehrlich pathway. The umbelliferone 6-dimethylallyltransferase (U6DT) and marmesin synthase (MS) involved in marmesin synthesis were optimized to increase marmesin production. Marmesin production was improved by truncating the transmembrane domains of PcU6DT, FcMS, and AtCPR1 and increasing the copy numbers of the genes encoding the truncated enzymes. Finally, a marmesin titer of 27.7 mg/L was obtained in shake flasks using the engineered yeast strain MU5. The constructed marmesin-producing strain provides the foundation for the green and large-scale production of pharmaceutically important furanocoumarins.

Two O-Methyltransferases Mediate Multiple Methylation Steps in the Biosynthesis of Coumarins in Cnidium monnieri

Coumarins with methoxy groups such as osthole (1), xanthotoxin (2), bergapten (3), and isopimpinellin (4) are typical bioactive ingredients of many medicinal plants. The methylation steps remain widely unknown. Herein, we report the discovery of two methyltransferases in the biosynthesis of O-methyl coumarins in Cnidium monnieri by transcriptome mining, heterologous expression, and in vitro enzymatic assays. The results reveal that (i) CmOMT1 catalyzes the methylation of osthenol (8) as the final step in the biosynthesis of 1, (ii) CmOMT2 shows the highest efficiency and preference for methylating xanthotoxol (11) to form 2, and (iii) CmOMT1 and CmOMT2 also efficiently transform bergaptol (10) and 8-hydroxybergapten (13) into 3 or 4, suggesting the CmOMTs mediate multistep methylations in the biosynthesis of linear furanocoumarins in C. monnieri.

Exceptional evolutionary lability of flower-like inflorescences (pseudanthia) in Apiaceae subfamily Apioideae

PREMISE

Pseudanthia are widespread and have long been postulated to be a key innovation responsible for some of the angiosperm radiations. The aim of our study was to analyze macroevolutionary patterns of these flower-like inflorescences and their potential correlation with diversification rates in Apiaceae subfamily Apioideae. In particular, we were interested to investigate evolvability of pseudanthia and evaluate their potential association with changes in the size of floral display.

METHODS

The framework for our analyses consisted of a time-calibrated phylogeny of 1734 representatives of Apioideae and a morphological matrix of inflorescence traits encoded for 847 species. Macroevolutionary patterns in pseudanthia were inferred using Markov models of discrete character evolution and stochastic character mapping, and a principal component analysis was used to visualize correlations in inflorescence architecture. The interdependence between net diversification rates and the occurrence of pseudocorollas was analyzed with trait-independent and trait-dependent approaches.

RESULTS

Pseudanthia evolved in 10 major clades of Apioideae with at least 36 independent origins and 46 reversals. The morphospace analysis recovered differences in color and compactness between floral and hyperfloral pseudanthia. A correlation between pseudocorollas and size of inflorescence was also strongly supported. Contrary to our predictions, pseudanthia are not responsible for variation in diversification rates identified in this subfamily.

CONCLUSIONS

Our results suggest that pseudocorollas evolve as an answer to the trade-off between enlargement of floral display and costs associated with production of additional flowers. The high evolvability and architectural differences in apioid pseudanthia may be explained on the basis of adaptive wandering and evolutionary developmental biology.

A new P450 involved in the furanocoumarin pathway underlies a recent case of convergent evolution

Furanocoumarins are phytoalexins often cited as an example to illustrate the arms race between plants and herbivorous insects. They are distributed in a limited number of phylogenetically distant plant lineages, but synthesized through a similar pathway, which raised the question of a unique or multiple emergence in higher plants. The furanocoumarin pathway was investigated in the fig tree (Ficus carica, Moraceae). Transcriptomic and metabolomic approaches led to the identification of CYP76F112, a cytochrome P450 catalyzing an original reaction. CYP76F112 emergence was inquired using phylogenetics combined with in silico modeling and site-directed mutagenesis. CYP76F112 was found to convert demethylsuberosin into marmesin with a very high affinity. This atypical cyclization reaction represents a key step within the polyphenol biosynthesis pathway. CYP76F112 evolutionary patterns suggests that the marmesin synthase activity appeared recently in the Moraceae family, through a lineage-specific expansion and diversification. The characterization of CYP76F112 as the first known marmesin synthase opens new prospects for the use of the furanocoumarin pathway. It also supports the multiple acquisition of furanocoumarin in angiosperms by convergent evolution, and opens new perspectives regarding the ability of cytochromes P450 to evolve new functions related to plant adaptation to their environment.

Plant cytochrome P450 plasticity and evolution

The superfamily of cytochrome P450 (CYP) enzymes plays key roles in plant evolution and metabolic diversification. This review provides a status on the CYP landscape within green algae and land plants. The 11 conserved CYP clans known from vascular plants are all present in green algae and several green algae-specific clans are recognized. Clan 71, 72, and 85 remain the largest CYP clans and include many taxa-specific CYP (sub)families reflecting emergence of linage-specific pathways. Molecular features and dynamics of CYP plasticity and evolution are discussed and exemplified by selected biosynthetic pathways. High substrate promiscuity is commonly observed for CYPs from large families, favoring retention of gene duplicates and neofunctionalization, thus seeding acquisition of new functions. Elucidation of biosynthetic pathways producing metabolites with sporadic distribution across plant phylogeny reveals multiple examples of convergent evolution where CYPs have been independently recruited from the same or different CYP families, to adapt to similar environmental challenges or ecological niches. Sometimes only a single or a few mutations are required for functional interconversion. A compilation of functionally characterized plant CYPs is provided online through the Plant P450 Database (erda.dk/public/vgrid/PlantP450/).

The evolution of coevolution in the study of species interactions

The study of reciprocal adaptation in interacting species has been an active and inspiring area of evolutionary research for nearly 60 years. Perhaps owing to its great natural history and potential consequences spanning population divergence to species diversification, coevolution continues to capture the imagination of biologists. Here we trace developments following Ehrlich and Raven's classic paper, with a particular focus on the modern influence of two studies by Dr. May Berenbaum in the 1980s. This series of classic work presented a compelling example exhibiting the macroevolutionary patterns predicted by Ehrlich and Raven and also formalized a microevolutionary approach to measuring selection, functional traits, and understanding reciprocal adaptation between plants and their herbivores. Following this breakthrough was a wave of research focusing on diversifying macroevolutionary patterns, mechanistic chemical ecology, and natural selection on populations within and across community types. Accordingly, we breakdown coevolutionary theory into specific hypotheses at different scales: reciprocal adaptation between populations within a community, differential coevolution among communities, lineage divergence, and phylogenetic patterns. We highlight progress as well as persistent gaps, especially the link between reciprocal adaptation and diversification.

Signatures of selection and evolutionary relevance of cytochrome P450s in plant-insect interactions

Enzymes in the cytochrome P450 (P450) superfamily have important functions ranging from those that are essential for the physiology and development of the individual to those that mediate interactions between individuals and their biotic environment. Until recently the study of P450s had focused on single functions, substrates, or pathways. Recent advances in sequencing, genome assembly, and phylogenetic methods have returned emphasis to the adaptive value of these enzymes in the context of herbivory. Comparisons of whole repertoires of P450s across related species reveal that P450s capable of metabolizing xenobiotics have an increased rate of gains compared to losses after gene duplications. In plants, studies have focused on enzymes and end-functions that have converged to provide increased resistance to herbivory. This review summarizes the latest findings related to the ecological value of P450s in the interactions between phytophagous insects and their host plants.

Biosynthesis and heterologous production of furanocoumarins: perspectives and current challenges

Covering: up to October 2020 Furanocoumarins are plant secondary metabolites used to treat several skin disorders, such as psoriasis and vitiligo, and also with other potential therapeutic activities. Furanocoumarins are extracted from plants where they accumulate in low amounts over long growth periods. In addition, their extraction and purification are difficult in an environmentally unfriendly and expensive process. Hence, new sustainable and greener production schemes able to overcome such limitations ought to be developed. While the heterologous production of simple coumarins has been demonstrated, the biosynthesis of more complex furanocoumarins remains greatly unexplored. Although several important steps of the pathway have been elucidated in the last decade, the complete pathway has not been completely unravelled. In this paper, we review the natural conversion of amino acids into furanocoumarins, as well as the heterologous expression of each enzyme of the pathway. We also explore the challenges that need to be addressed so that their heterologous production can become a viable alternative.

A GDSL lipase-like from Ipomoea batatas catalyzes efficient production of 3,5-diCQA when expressed in Pichia pastoris

The synthesis of 3,5-dicaffeoylquinic acid (3,5-DiCQA) has attracted the interest of many researchers for more than 30 years. Recently, enzymes belonging to the BAHD acyltransferase family were shown to mediate its synthesis, albeit with notably low efficiency. In this study, a new enzyme belonging to the GDSL lipase-like family was identified and proven to be able to transform chlorogenic acid (5-O-caffeoylquinic acid, 5-CQA, CGA) in 3,5-DiCQA with a conversion rate of more than 60%. The enzyme has been produced in different expression systems but has only been shown to be active when transiently synthesized in Nicotiana benthamiana or stably expressed in Pichia pastoris. The synthesis of the molecule could be performed in vitro but also by a bioconversion approach beginning from pure 5-CQA or from green coffee bean extract, thereby paving the road for producing it on an industrial scale.

Miguel et al. identify a new enzyme belonging to the GDSL lipase-like family that is involved in the final stage of transformation of 5-CQA into 3,5-diCQA. This enzyme is able to realize an efficient transformation by over 60%, making the transformation process a valuable technological tool that can be easily transferred on an industrial scale.

Root-Secreted Coumarins and the Microbiota Interact to Improve Iron Nutrition in Arabidopsis

Plants benefit from associations with a diverse community of root-colonizing microbes. Deciphering the mechanisms underpinning these beneficial services are of interest for improving plant productivity. We report a plant-beneficial interaction between Arabidopsis thaliana and the root microbiota under iron deprivation that is dependent on the secretion of plant-derived coumarins. Disrupting this pathway alters the microbiota and impairs plant growth in iron-limiting soil. Furthermore, the microbiota improves iron-limiting plant performance via a mechanism dependent on plant iron import and secretion of the coumarin fraxetin. This beneficial trait is strain specific yet functionally redundant across phylogenetic lineages of the microbiota. Transcriptomic and elemental analyses revealed that this interaction between commensals and coumarins promotes growth by relieving iron starvation. These results show that coumarins improve plant performance by eliciting microbe-assisted iron nutrition. We propose that the bacterial root microbiota, stimulated by secreted coumarins, is an integral mediator of plant adaptation to iron-limiting soils.

Convergent evolution of the UbiA prenyltransferase family underlies the independent acquisition of furanocoumarins in plants

Furanocoumarins (FCs) are plant-specialized metabolites with potent allelochemical properties. The distribution of FCs is scattered with a chemotaxonomical tendency towards four distant families with highly similar FC pathways. The mechanism by which this pathway emerged and spread in plants has not been elucidated. Furanocoumarin biosynthesis was investigated in Ficus carica (fig, Moraceae), focusing on the first committed reaction catalysed by an umbelliferone dimethylallyltransferase (UDT). Comparative RNA-seq analysis among latexes of different fig organs led to the identification of a UDT. The phylogenetic relationship of this UDT to previously reported Apiaceae UDTs was evaluated. The expression pattern of F. carica prenyltransferase 1 (FcPT1) was related to the FC contents in different latexes. Enzymatic characterization demonstrated that one of the main functions of FcPT1 is UDT activity. Phylogenetic analysis suggested that FcPT1 and Apiaceae UDTs are derived from distinct ancestors, although they both belong to the UbiA superfamily. These findings are supported by significant differences in the related gene structures. This report describes the identification of FcPT1 involved in FC biosynthesis in fig and provides new insights into multiple origins of the FC pathway and, more broadly, into the adaptation of plants to their environments.

Convergent evolution leading to the appearance of furanocoumarins in citrus plants

Furanocoumarins are defense molecules mainly described in four plant families that are phylogenetically distant. Molecular characterization of the biosynthetic pathway has been started for many years in Apiaceae and Rutaceae. The results obtained thus far in Apiaceae indicated a major role of cytochromes P450 (P450s) in the CYP71 family. In the present work, we describe the importance of another subfamily of P450s, CYP82D, identified by using a deep analysis of the citrus (Rutaceae) genome and microarray database. CYP82D64 is able to hydroxylate xanthotoxin to generate 5-OH-xanthotoxin. Minor and limited amino acid changes in the CYP82D64 coding sequence between Citrus paradisi and Citrus hystrix provide the enzyme in the latter with the ability to hydroxylate herniarin, but with low efficiency. The kinetic constants of the enzyme are consistent with those of other enzymes of this type in plants and indicate that it may be the physiological substrate. The activity of the enzyme is identical to that of CYP71AZ6 identified in parsnip, showing possible evolutionary convergence between these two families of plants. It is highly possible that these molecules are derived from the synthesis of ubiquitous coumarins throughout the plant kingdom.

Substrate-specificity of cytochrome P450-mediated detoxification as an evolutionary strategy for specialization on furanocoumarin-containing hostplants: CYP6AE89 in parsnip webworms

The parsnip webworm, Depressaria pastinacella, is restricted to two hostplant genera containing six structurally diverse furanocoumarins. Of these, imperatorin is detoxified by a specialized cytochrome P450, CYP6AB3. A previous whole-larva transcriptome analysis confirmed the presence of nine transcripts that belong to the CYP6AE subfamily. Here, by examining midgut-specific gene expression patterns we determined that CYP6AE89 transcripts were highly expressed and furanocoumarin-inducible. Computer docking and energy-minimization of a CYP6AE89 model with all six furanocoumarins showed that 5-methoxylated bergapten and 8-methoxylated xanthotoxin had the smallest distances from the heme to the proton-donor residue in the catalytic I-helix, and that the 5,8-dimethoxylated isopimpinellin and bergapten had the smallest energy-minimized distance from the heme oxygen to the furan ring double bond. To evaluate this prediction, we expressed the CYP6AE89 protein in an Escherichia coli system, and used it to detect high catalytic activity against the two mono-methoxylated linear furanocoumarins - bergapten and xanthotoxin - and weak activity against isopimpinellin. Thus, CYP6AE89, like CYP6AB3, is probably specialized for detoxifying only a subset of hostplant furanocoumarins. A maximum-likelihood tree built with six representative lepidopterans with manually annotated cytochrome P450s shows that CYP6AE89 may have evolved much faster than the other CYP6AE proteins, possibly indicative of host selection pressure.

Untargeted Metabolomics Approach Reveals Diverse Responses of Pastinaca Sativa to Ozone and Wounding Stresses

Stresses such as wounding or atmospheric pollutant exposure have a significant impact on plant fitness. Since it has been widely described that the metabolome directly reflects plant physiological status, a way to assess this impact is to perform a global metabolomic analysis. In this study, we investigated the effect of two abiotic stresses (mechanical wounding and ozone exposure) on parsnip metabolic balance using a liquid chromatography-mass spectrometry-based untargeted metabolomic approach. For this purpose, parsnip leaves were submitted to an acute ozone exposure or were mechanically wounded and sampled 24, 48, and 72 h post-treatment. Multivariate and univariate statistical analyses highlighted numerous differentially-accumulated metabolic features as a function of time and treatment. Mechanical wounding led to a more differentiated response than ozone exposure. We found that the levels of coumarins and fatty acyls increased in wounded leaves, while flavonoid concentration decreased in the same conditions. These results provide an overview of metabolic destabilization through differentially-accumulated compounds and provide a better understanding of global plant metabolic changes in defense mechanisms.

Defence mechanisms of Ficus: pyramiding strategies to cope with pests and pathogens

Ficus species have adapted to diverse environments and pests by developing physical or chemical protection strategies. Physical defences are based on the accumulation of minerals such as calcium oxalate crystals, amorphous calcium carbonates and silica that lead to tougher plants. Additional cellular structures such as non-glandular trichomes or laticifer cells make the leaves rougher or sticky upon injury. Ficus have also established structures that are able to produce specialized metabolites (alkaloids, terpenoids, and phenolics) or proteins (proteases, protease inhibitors, oxidases, and chitinases) that are toxic to predators. All these defence mechanisms are distributed throughout the plant and can differ depending on the genotype, the stage of development or the environment. In this review, we present an overview of these strategies and discuss how these complementary mechanisms enable effective and flexible adaptation to numerous hostile environments.