Evolution of substrate recognition sites (SRSs) in cytochromes P450 from Apiaceae exemplified by the CYP71AJ subfamily

MicroRNA PC-5p-3991_515 mediates triflumezopyrim susceptibility in the small brown planthopper through regulating the post-transcriptional expression of P450 CYP417A2

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are recognized as a major contributor to metabolic resistance in insects to most insecticides, through gene overexpressions and protein mutations. MicroRNA (miRNA), an important post-transcriptional regulator, has been reported to promote insecticide resistance by mediating the expression of detoxification enzyme genes.

RESULTS

In the present study, we reported that a novel microRNA PC-5p-3991_515 was involved in the post-transcriptional regulation of CYP417A2 and mediated the triflumezopyrim susceptibility in the small brown planthopper (SBPH), Laodelphax striatellus (Fallén). The tissue expression profiles showed that CYP417A2 was highly expressed in fat body. CYP417A2 was significantly up-regulated at 12, 36, 60, 84 and 108 h after the triflumezopyrim treatment. RNA interference (RNAi) against CYP417A2 significantly increased triflumezopyrim susceptibility in SBPH. According to the prediction by miRanda and TargetScan software, three miRNAs were indicated to bind to CYP417A2. However, when oversupply of agomir, only two miRNAs, PC-3p-625_4405 and PC-5p-3991_515, significantly increased the susceptibility to triflumezopyrim and decreased CYP417A2 levels. Furthermore, PC-5p-3991_515 was confirmed to be involved in the post-transcriptional regulation of CYP417A2 by dual luciferase reporter assay. Meanwhile, PC-5p-3991_515 was co-localized with CYP417A2 in the midgut in situ hybridization.

CONCLUSION

Our findings revealed that the novel microRNA, PC-5p-3991_515, post-transcriptionally regulated CYP417A2 expression, which then mediated the triflumezopyrim susceptibility in SBPH. © 2023 Society of Chemical Industry.

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.

The genome of the toxic invasive species Heracleum sosnowskyi carries an increased number of genes despite absence of recent whole-genome duplications

Heracleum sosnowskyi, belonging to a group of giant hogweeds, is a plant with large effects on ecosystems and human health. It is an invasive species that contributes to the deterioration of grassland ecosystems. The ability of H. sosnowskyi to produce linear furanocoumarins (FCs), photosensitizing compounds, makes it very dangerous. At the same time, linear FCs are compounds with high pharmaceutical value used in skin disease therapies. Despite this high importance, it has not been the focus of genetic and genomic studies. Here, we report a chromosome-scale assembly of Sosnowsky's hogweed genome. Genomic analysis revealed an unusually high number of genes (55106) in the hogweed genome, in contrast to the 25-35 thousand found in most plants. However, we did not find any traces of recent whole-genome duplications not shared with its confamiliar, Daucus carota (carrot), which has approximately thirty thousand genes. The analysis of the genomic proximity of duplicated genes indicates on tandem duplications as a main reason for this increase. We performed a genome-wide search of the genes of the FC biosynthesis pathway and surveyed their expression in aboveground plant parts. Using a combination of expression data and phylogenetic analysis, we found candidate genes for psoralen synthase and experimentally showed the activity of one of them using a heterologous yeast expression system. These findings expand our knowledge on the evolution of gene space in plants and lay a foundation for further analysis of hogweed as an invasive plant and as a source of FCs.

Status of Fungicide Resistance and Physiological Characterization of Tebuconazole Resistance in Rhizocotonia solani in Sichuan Province, China

The resistance prevalence of chemical fungicides has caused increasingly serious agro-ecological environmental problems. However, there are few previous reports about resistance to succinate dehydrogenase (SDHI) or sterol demethylation inhibitor (DMI) in Rhizoctonia solani, one of the main agro-diseases. In this study, the fungicide resistance of 122 R. solani isolates in Sichuan Province was monitored by the mycelial growth rate method. Results showed that all isolates were susceptible to hexaconazole and most isolates were susceptible to thifluzamide, except for the field isolate MSRS-2-7 due to a moderate resistance to thifluzamide (16.43-fold resistance ratio, RR), compared to the sensitivity baseline of thifluzamide (0.042 μg/mL EC₅₀ values). On the contrary, many isolates showed moderate or high resistance to tebuconazole (10.59- to 60.78-fold RR), reaching EC₅₀ values of 0.54~3.10 μg/mL, especially for a highly resistant isolate LZHJ-1-8 displaying moderate resistance to epoxiconazole (35.40-fold RR due to a 3.54 μg/mL EC₅₀ value). The fitness determination found that the tebuconazole-resistant isolates showed higher fitness cost with these characteristics, including a lower growth rate, higher relative electric conductivity, an increased ability to tolerate tebuconazole, and high osmotic pressure. Four new mutations of cytochrome P450 sterol 14α-demethylase (CYP51), namely, S94A, N406S, H793R, and L750P, which is the target for DMI fungicides, was found in the tebuconazole-resistant isolates. Furthermore, the lowest binding energy with tebuconazole was also found in the LZHJ-1-8 isolate possessing all the mutations through analyses with Discovery Studio software. Therefore, these new mutation sites of CYP51 may be linked to the resistance against tebuconazole, and its application for controlling R. solani should be restricted in some areas.

Comparison of herbicide specificity of CYP81A cytochrome P450s from rice and a multiple-herbicide resistant weed, Echinochloa phyllopogon

BACKGROUND

CYP81A cytochrome P450s (CYP81As) play a key role in herbicide detoxification in Poaceae plants. Crop CYP81As confer natural tolerance to multiple herbicides, whereas CYP81As in weeds disrupt herbicide action. Identifying differences in CYP81A herbicide specificity between crops and weeds could provide valuable information for controlling weeds. In this study, we quantitatively compared herbicide specificity between CYP81A6 from rice (Oryza sativa) and CYP81A12 and CYP81A21 from a weed, Echinochloa phyllopogon, using transgenic Escherichia coli and Arabidopsis.

RESULTS

All three CYP81As metabolized the five tested herbicides and formed similar metabolites, with the highest relative activities of 400 to 580% toward bentazone compared to their activity on bensulfuron-methyl (defined as 100%). However, they showed differing activity toward propyrisulfuron. The relative activities of Echinochloa phyllopogon CYP81A12 (12.2%) and CYP81A21 (34.4%) toward propyrisulfuron were lower than that of rice CYP81A6 (98.5%). Additionally, rice CYP81A6 produced O-demethylated propyrisulfuron and hydroxylated products, whereas Echinochloa phyllopogon CYP81As produced only hydroxylated products. Arabidopsis expressing CYP81A12 and CYP81A21 exhibited lower levels of resistance against propyrisulfuron compared to that in Arabidopsis expressing CYP81A6. Homology modeling and in silico docking revealed that bensulfuron-methyl docked well into the active centers of all three CYP81As, whereas propyrisulfuron docked into rice CYP81A6 but not into Echinochloa phyllopogon CYP81As.

CONCLUSION

The differing herbicide specificity displayed by rice CYP81A6 and Echinochloa phyllopogon CYP81A12 and CYP81A21 will help design inhibitors (synergists) of weed CYP81As, as well as develop novel herbicide ingredients that are selectively metabolized by crop CYP81As, to overcome the problem of herbicide resistance. © 2022 Society of Chemical Industry.

Genome-Wide Evolutionary Analysis of Putative Non-Specific Herbicide Resistance Genes and Compilation of Core Promoters between Monocots and Dicots

Herbicides are key weed-control tools, but their repeated use across large areas has favored the evolution of herbicide resistance. Although target-site has been the most prevalent and studied type of resistance, non-target-site resistance (NTSR) is increasing. However, the genetic factors involved in NTSR are widely unknown. In this study, four gene groups encoding putative NTSR enzymes, namely, cytochrome-P450, glutathione-S-transferase (GST), uridine 5'-diphospho-glucuronosyltransferase (UDPGT), and nitronate monooxygenase (NMO) were analyzed. The monocot and dicot gene sequences were downloaded from publicly available databases. Phylogenetic trees revealed that most of the CYP450 resistance-related sequences belong to CYP81 (5), and in GST, most of the resistance sequences belonged to GSTU18 (9) and GSTF6 (8) groups. In addition, the study of upstream promoter sequences of these NTSR genes revealed stress-related cis-regulatory motifs, as well as eight transcription factor binding sites (TFBS) were identified. The discovered TFBS were commonly present in both monocots and dicots, and the identified motifs are known to play key roles in countering abiotic stress. Further, we predicted the 3D structure for the resistant CYP450 and GST protein and identified the substrate recognition site through the homology approach. Our description of putative NTSR enzymes may be used to develop innovative weed control techniques to delay the evolution of NTSR.

Characterization of Cytochrome P450s with Key Roles in Determining Herbicide Selectivity in Maize

Safeners such as metcamifen and benoxacor are widely used in maize to enhance the selectivity of herbicides through the induction of key detoxifying enzymes, notably cytochrome P450 monooxygenases (CYPs). Using a combination of transcriptomics, proteomics, and functional assays, the safener-inducible CYPs responsible for herbicide metabolism in this globally important crop have been identified. A total of 18 CYPs belonging to clans 71, 72, 74, and 86 were safener-induced, with the respective enzymes expressed in yeast and screened for activity toward thiadiazine (bentazon), sulfonylurea (nicosulfuron), and triketone (mesotrione and tembotrione) chemistries. Herbicide metabolism was largely restricted to family CYP81A members from clan 71, notably CYP81A9, CYP81A16, and CYP81A2. Quantitative transcriptomics and proteomics showed that CYP81A9/CYP81A16 were dominant enzymes in safener-treated field maize, whereas only CYP81A9 was determined in sweet corn. The relationship between CYP81A sequence and activities were investigated by splicing CYP81A2 and CP81A9 together as a series of recombinant chimeras. CYP81A9 showed wide ranging activities toward the three herbicide chemistries, while CYP81A2 uniquely hydroxylated bentazon in multiple positions. The plasticity in substrate specificity of CYP81A9 toward multiple herbicides resided in the second quartile of its N terminal half. Further phylogenetic analysis of CYP81A9 showed that the maize enzyme was related to other CYP81As linked to agrochemical metabolism in cereals and wild grasses, suggesting this clan 71 CYP has a unique function in determining herbicide selectivity in arable crops.

First purified recombinant CYP75B including transmembrane helix with unexpected high substrate specificity to (2R)-naringenin

Anthochlor pigments (chalcones and aurones) play an important role in yellow flower colourization, the formation of UV-honey guides and show numerous health benefits. The B-ring hydroxylation of chalcones is performed by membrane bound cytochrome P450 enzymes. It was assumed that usual flavonoid 3′-hydroxlases (F3′Hs) are responsible for the 3,4- dihydroxy pattern of chalcones, however, we previously showed that a specialized F3′H, namely chalcone 3-hydroxylase (CH3H), is necessary for the hydroxylation of chalcones. In this study, a sequence encoding membrane bound CH3H from Dahlia variabilis was recombinantly expressed in yeast and a purification procedure was developed. The optimized purification procedure led to an overall recovery of 30% recombinant DvCH3H with a purity of more than 84%. The enzyme was biochemically characterized with regard to its kinetic parameters on various substrates, including racemic naringenin, as well as its enantiomers (2S)-, and (2R)-naringenin, apigenin and kaempferol. We report for the first time the characterization of a purified Cytochrome P450 enzyme from the flavonoid biosynthesis pathway, including the transmembrane helix. Further, we show for the first time that recombinant DvCH3H displays a higher affinity for (2R)-naringenin than for (2S)-naringenin, although (2R)-flavanones are not naturally formed by chalcone isomerase.

Enhanced Accumulation of Biologically Active Coumarin and Furanocoumarins in Callus Culture and Field-grown Plants of Ruta chalepensis Through LED Light-treatment

Ruta chalepensis, a medicinal plant, produces biologically active coumarins (CRs) and furanocoumarins (FCRs). However, their yield is quite low in cultivated plants. In this work, the influence of light emitting diodes (LEDs) was investigated on the accumulation of CRs and FCRs in the callus cultures and field-grown plants of R. chalepensis. Among the various tested wavelengths of LED lights, maximum accumulation of CR and FCRs was recorded under blue LED treatment in both the callus cultures as well as field-grown plants as compared to respective controls treated with white LED. Metabolite analyses of LED-treated field-grown plants showed that highest concentrations of CR (umbelliferone, 2.8-fold), and FCRs (psoralen, 2.3-fold; xanthotoxin, 3.8-fold; bergapten, 1.16-fold) were accumulated upon blue LED-treatment for six days. CR and FCRs contents were also analyzed in the blue- and red-LED-treated in vitro callus tissue. Upon blue LED-treatment, callus accumulated significantly high levels of umbelliferone (48.6 ± 1.2 µg/g DW), psoralen (370.12 ± 10.6 µg/g DW) and xanthotoxin (10.16 ± 0.48 µg/g DW). These findings imply that blue LED-treatment is a viable option as a non-invasive and low-cost elicitation technology for the enhanced production of biologically active CR and FCRs in field-grown plants and callus cultures of R. chalepensis.

Cytochrome P450-catalyzed biosynthesis of furanoditerpenoids in the bioenergy crop switchgrass (Panicum virgatum L.)

Specialized diterpenoid metabolites are important mediators of plant-environment interactions in monocot crops. To understand metabolite functions in plant environmental adaptation that ultimately can enable crop improvement strategies, a deeper knowledge of the underlying species-specific biosynthetic pathways is required. Here, we report the genomics-enabled discovery of five cytochrome P450 monooxygenases (CYP71Z25-CYP71Z29) that form previously unknown furanoditerpenoids in the monocot bioenergy crop Panicum virgatum (switchgrass). Combinatorial pathway reconstruction showed that CYP71Z25-CYP71Z29 catalyze furan ring addition directly to primary diterpene alcohol intermediates derived from distinct class II diterpene synthase products. Transcriptional co-expression patterns and the presence of select diterpenoids in switchgrass roots support the occurrence of P450-derived furanoditerpenoids in planta. Integrating molecular dynamics, structural analysis and targeted mutagenesis identified active site determinants that contribute to the distinct catalytic specificities underlying the broad substrate promiscuity of CYP71Z25-CYP71Z29 for native and non-native diterpenoids.

A cytochrome P450 CYP81AM1 from Tripterygium wilfordii catalyses the C-15 hydroxylation of dehydroabietic acid

A novel cytochrome P450 from Tripterygium wilfordii, CYP81AM1, specifically catalyses the C-15 hydroxylation of dehydroabietic acid. This is the first CYP450 to be found in plants with this function. Cytochrome P450 oxygenases (CYPs) play an important role in the post-modification in biosynthesis of plant bioactive terpenoids. Here, we found that CYP81AM1 can catalyze the formation of 15-hydroxydehydroabietic acid by in vitro enzymatic reactions and in vivo yeast feeding assays. This is the first study to show that CYP81 family enzymes are involved in the hydroxylation of abietane diterpenoids. At the same time, we found that CYP81AM1 could not catalyse abietatriene and dehydroabietinol, suggesting that the occurrence of the reaction may be related to the carboxyl group. Through molecular docking and site mutations, it was found that some amino acid sites (F104, K107) near the carboxyl group had an important effect on enzyme activity, also suggesting that the carboxyl group played an important role in the occurrence of the reaction.

Impact of Genomic and Transcriptomic Resources on Apiaceae Crop Breeding Strategies

The Apiaceae taxon is one of the most important families of flowering plants and includes thousands of species used for food, flavoring, fragrance, medical and industrial purposes. This study had the specific intent of reviewing the main genomics and transcriptomic data available for this family and their use for the constitution of new varieties. This was achieved starting from the description of the main reproductive systems and barriers, with particular reference to cytoplasmic (CMS) and nuclear (NMS) male sterility. We found that CMS and NMS systems have been discovered and successfully exploited for the development of varieties only in Foeniculum vulgare, Daucus carota, Apium graveolens and Pastinaca sativa; whereas, strategies to limit self-pollination have been poorly considered. Since the constitution of new varieties benefits from the synergistic use of marker-assisted breeding in combination with conventional breeding schemes, we also analyzed and discussed the available SNP and SSR marker datasets (20 species) and genomes (8 species). Furthermore, the RNA-seq studies aimed at elucidating key pathways in stress tolerance or biosynthesis of the metabolites of interest were limited and proportional to the economic weight of each species. Finally, by aligning 53 plastid genomes from as many species as possible, we demonstrated the precision offered by the super barcoding approach to reconstruct the phylogenetic relationships of Apiaceae species. Overall, despite the impressive size of this family, we documented an evident lack of molecular data, especially because genomic and transcriptomic resources are circumscribed to a small number of species. We believe that our contribution can help future studies aimed at developing molecular tools for boosting breeding programs in crop plants of the Apiaceae family.

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.

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.

Roles of the variable P450 substrate recognition sites SRS1 and SRS6 in esfenvalerate metabolism by CYP6AE subfamily enzymes in Helicoverpa armigera

The cotton bollworm P450s of the clustered CYP6AE subfamily share high sequence identities but differ dramatically in their capacity to metabolize xenobiotics, especially esfenvalerate. Among them, CYP6AE17 has the highest sequence identity with CYP6AE18 but shows ~7-fold higher metabolic efficiency. CYP6AE11 is most active towards esfenvalerate but CYP6AE20 is inactive even though the enzymes share 54.8% sequence identity. Sequence analysis revealed the SRS1 (Substrate Recognition Site) and SRS6 between CYP6AE17 and CYP6AE18, and SRS1 between CYP6AE11 and CYP6AE20 are the most variable among all six SRSs. In order to identify the key factors that underlie the observed catalytic difference, we exchanged these SRS sequences between two pairs of P450s and studied the activity of the resulting hybrid mutants or chimeras. In vitro metabolism showed that the CYP6AE17/18 chimeras had 2- and 14-fold decreased activities and the CYP6AE18/17 chimeras had 6- and 10-fold increased activities to esfenvalerate. Meanwhile, after exchanging SRS1 with each other, the CYP6AE11/20 chimera folded incorrectly but the CYP6AE20/11 chimera gained moderate activity to esfenvalerate. Molecular modelling showed that amino acids variants within SRS1 or SRS6 change the shape and chemical environment of the active sites, which may affect the ligand-binding interactions. These results indicate that the protein structure variation resulting from the sequence diversity of SRSs promotes the evolution of insect chemical defense and contributes to the development of insect resistance to pesticides.

Convergent recruitment of 5'-hydroxylase activities by CYP75B flavonoid B-ring hydroxylases for tricin biosynthesis in Medicago legumes

Tricin (3',5'-dimethoxylated flavone) is a predominant flavonoid amongst monocots but occurs only in isolated and unrelated dicot lineages. Although tricin biosynthesis has been intensively studied in monocots, it has remained largely elusive in tricin-accumulating dicots. We investigated a subgroup of cytochrome P450 (CYP) 75B subfamily flavonoid B-ring hydroxylases (FBHs) from two tricin-accumulating legumes, Medicago truncatula and alfalfa (Medicago sativa), by phylogenetic, molecular, biochemical and mutant analyses. Five Medicago cytochrome P450 CYP75B FBHs are phylogenetically distant from other legume CYP75B members. Among them, MtFBH-4, MsFBH-4 and MsFBH-10 were expressed in tricin-accumulating vegetative tissues. In vitro and in planta analyses demonstrated that these proteins catalyze 3'- and 5'-hydroxylations critical to tricin biosynthesis. A key amino acid polymorphism, T492G, at their substrate recognition site 6 domain is required for the novel 5'-hydroxylation activities. Medicago truncatula mtfbh-4 mutants were tricin-deficient, indicating that MtFBH-4 is indispensable for tricin biosynthesis. Our results revealed that these Medicago legumes had acquired the tricin pathway through molecular evolution of CYP75B FBHs subsequent to speciation from other nontricin-accumulating legumes. Moreover, their evolution is independent of that of grass-specific CYP75B apigenin 3'-hydroxylases/chrysoeriol 5'-hydroxylases dedicated to tricin production and Asteraceae CYP75B flavonoid 3',5'-hydroxylases catalyzing the production of delphinidin-based pigments.

Two CYP71AJ enzymes function as psoralen synthase and angelicin synthase in the biosynthesis of furanocoumarins in Peucedanum praeruptorum Dunn

Psoralen synthase and angelicin synthase responsible for the formation of psoralen and angelicin in Peucedanum praeruptorum Dunn were identified and functionally characterized, respectively. Furanocoumarins were reported to possess several activities such as anticancer, anti-inflammatory and neuroprotective, and function as phytotoxin and allelochemical in plants. Furanocoumarins are the main bioactive ingredient in P. praeruptorum which is a commonly used traditional Chinese medicine. Phenylalanine ammonia lyase (PAL), 4-coumarate: CoA ligase (4CL), p-coumaroyl CoA 2'-hyfroxylase (C2'H) were cloned previously to elucidate the biosynthetic mechanism of coumarin lactone ring. However, the genes involved in complex coumarins in P. praeruptorum have not been explored. Herein, putative psoralen synthase CYP71AJ49 and angelicin synthase CYP71AJ51 were cloned from P. praeruptorum. In vivo and in vitro yeast assays were conducted to confirm their activities. Furthermore, the results of High Performance Liquid Chromatography-Electrospray Ionization Mass Spectrometry (HPLC-ESI-MS) verified that CYP71AJ49 catalyzed the conversion of marmesin to psoralen, and CYP71AJ51 catalyzed columbianetin to angelicin. Subsequently, the expression profile showed that CYP71AJ49 and CYP71AJ51 were easily affected by environmental conditions, especially UV and temperature. The genes tissue-specific expression and compounds tissue-specific distribution pattern indicated the existence of substance transport in P. praeruptorum. Phylogenetic analysis was conducted with 27 CYP71AJs, CYP71AJ49 and CYP71AJ51 were classified in I-4 and I-2, respectively. These results provide further insight to understand the biosynthetic mechanism of complex coumarins.

CYP71BL9, the missing link in costunolide synthesis of sunflower

Despite intensive research in recent years, the biosynthetic route to costunolide in sunflower so far remained obscured. Additional P450 sequences from public sunflower transcriptomic database were screened to search for candidate enzymes which are able to introduce the 6α-hydroxy-group required for the esterification with the carboxy group of germacarane A acid, the final step in costunolide formation. CYP71BL9, a new P450 enzyme from sunflower was shown to catalyze this hydroxylation, hence being identified as HaCOS. Phylogentically, HaCOS is closer related to HaG8H than to any other known costunolide synthase in Asteraceae.The enzyme was successfully employed to reconstruct the sunflower biosynthesis of costunolide in transformed tobacco. Contrary, in yeast, only minor amounts of sesquiterpene lactone was produced, while 5-hydroxyfarnesylic acid was formed instead. HaCOS in combination with HaG8H produced 8β-hydroxycostunolide (eupatolide) in transformed plants, thus indicating that sunflower possesses two independent modes of eupatolide synthesis via HaCOS and via HaES. The lack of HaCOS expression and of costunolide in trichomes suggests that the enzyme triggers the costunolied synthesis of the inner tissues of sunflower and might be linked to growth regulation processes.

Angelicin-A Furocoumarin Compound With Vast Biological Potential

Angelicin, a member of the furocoumarin group, is related to psoralen which is well known for its effectiveness in phototherapy. The furocoumarins as a group have been studied since the 1950s but only recently has angelicin begun to come into its own as the subject of several biological studies. Angelicin has demonstrated anti-cancer properties against multiple cell lines, exerting effects via both the intrinsic and extrinsic apoptotic pathways, and also demonstrated an ability to inhibit tubulin polymerization to a higher degree than psoralen. Besides that, angelicin too demonstrated anti-inflammatory activity in inflammatory-related respiratory and neurodegenerative ailments via the activation of NF-κB pathway. Angelicin also showed pro-osteogenesis and pro-chondrogenic effects on osteoblasts and pre-chondrocytes respectively. The elevated expression of pro-osteogenic and chondrogenic markers and activation of TGF-β/BMP, Wnt/β-catenin pathway confirms the positive effect of angelicin bone remodeling. Angelicin also increased the expression of estrogen receptor alpha (ERα) in osteogenesis. Other bioactivities, such as anti-viral and erythroid differentiating properties of angelicin, were also reported by several researchers with the latter even displaying an even greater aptitude as compared to the commonly prescribed drug, hydroxyurea, which is currently on the market. Apart from that, recently, a new application for angelicin against periodontitis had been studied, where reduction of bone loss was indirectly caused by its anti-microbial properties. All in all, angelicin appears to be a promising compound for further studies especially on its mechanism and application in therapies for a multitude of common and debilitating ailments such as sickle cell anaemia, osteoporosis, cancer, and neurodegeneration. Future research on the drug delivery of angelicin in cancer, inflammation and erythroid differentiation models would aid in improving the bioproperties of angelicin and efficacy of delivery to the targeted site. More in-depth studies of angelicin on bone remodeling, the pro-osteogenic effect of angelicin in various bone disease models and the anti-viral implications of angelicin in periodontitis should be researched. Finally, studies on the binding of angelicin toward regulatory genes, transcription factors, and receptors can be done through experimental research supplemented with molecular docking and molecular dynamics simulation.

Uncovering the evolutionary origin of blue anthocyanins in cereal grains

Functional divergence after gene duplication plays a central role in plant evolution. Among cereals, only Hordeum vulgare (barley), Triticum aestivum (wheat) and Secale cereale (rye) accumulate delphinidin-derived (blue) anthocyanins in the aleurone layer of grains, whereas Oryza sativa (rice), Zea mays (maize) and Sorghum bicolor (sorghum) do not. The underlying genetic basis for this natural occurrence remains elusive. Here, we mapped the barley Blx1 locus involved in blue aleurone to an approximately 1.13 Mb genetic interval on chromosome 4HL, thus identifying a trigenic cluster named MbHF35 (containing HvMYB4H, HvMYC4H and HvF35H). Sequence and expression data supported the role of these genes in conferring blue-coloured (blue aleurone) grains. Synteny analyses across monocot species showed that MbHF35 has only evolved within distinct Triticeae lineages, as a result of dispersed gene duplication. Phylogeny analyses revealed a shared evolution pattern for MbHF35 in Triticeae, suggesting that these genes have co-evolved together. We also identified a Pooideae-specific flavonoid 3',5'-hydroxylase (F3'5'H) lineage, termed here Mo_F35H2, which has a higher amino acid similarity with eudicot F3'5'Hs, demonstrating a scenario of convergent evolution. Indeed, selection tests identified 13 amino acid residues in Mo_F35H2 that underwent positive selection, possibly driven by protein thermostablility selection. Furthermore, through the interrogation of barley germplasm there is evidence that HvMYB4H and HvMYC4H have undergone human selection. Collectively, our study favours blue aleurone as a recently evolved trait resulting from environmental adaptation. Our findings provide an evolutionary explanation for the absence of blue anthocyanins in other cereals and highlight the importance of gene functional divergence for plant diversity and environmental adaptation.

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.

Catalysis by allene oxide synthases (CYP74A and CYP74C): Alterations by the Phe/Leu mutation at the SRS-1 region

The CYP74 family of cytochromes P450 includes four fatty acid hydroperoxide metabolizing enzymes: allene oxide synthase (AOS), hydroperoxide lyase (HPL), divinyl ether synthase, and epoxyalcohol synthase (EAS). All P450s have six substrate recognition sites (SRSs) in their structures. Some CYP74 mutations in SRSs leading to their interconversions including substitutions in "F/L toggle" (SRS-1 region) were reported before. For further elucidation of the role of this site in CYP74 catalysis, the effect of Phe/Leu mutation on the specificity of selected AOSs was studied in the present work. Mutant forms of ZmAOS1 (CYP74A19, Zea mays), LeAOS3 (CYP74C3, Lycopersicon esculentum), and PpAOS2 (CYP74A8, Physcomitrella patens) acquired partial EAS activity. Mutant forms of ZmAOS1 and PpAOS2 possessed additional HPL activities. The results validate the significance of the "F/L toggle" as a catalytic determinant of CYP74s, as well as the importance of the conserved Phe at this site for the AOS catalysis.

High-throughput metabolomic and transcriptomic analyses vet the potential route of cerpegin biosynthesis in two varieties of Ceropegia bulbosa Roxb

Exploration with high-throughput transcriptomics and metabolomics of two varieties of Ceropegia bulbosa identifies candidate genes, crucial metabolites and a potential cerpegin biosynthetic pathway. Ceropegia bulbosa is an important medicinal plant, used in the treatment of various ailments including diarrhea, dysentery, and syphilis. This is primarily attributed to the presence of pharmaceutically active secondary metabolites, especially cerpegin. As this plant belongs to an endemic threatened category, genomic resources are not available hampering exploration on the molecular basis of cerpegin accumulation till now. Therefore, we undertook high-throughput metabolomic and transcriptomic analyses using different tissues from two varieties namely, C. bulbosa var. bulbosa and C. bulbosa var. lushii. Metabolomic analysis revealed spatial and differential accumulation of various metabolites. We chemically synthesized and characterized the cerpegin and its derivatives by liquid chromatography tandem-mass spectrometry (LC-MS/MS). Importantly, these comparisons suggested the presence of cerpegin and 5-allyl cerpegin in all C. bulbosa tissues. Further, de novo transcriptome analysis indicated the presence of significant transcripts for secondary metabolic pathways through the Kyoto encyclopedia of genes and genomes database. Tissue-specific profiling of transcripts and metabolites showed a significant correlation, suggesting the intricate mechanism of cerpegin biosynthesis. The expression of potential candidate genes from the proposed cerpegin biosynthetic pathway was further validated by qRT-PCR and NanoString nCounter. Overall, our findings propose a potential route of cerpegin biosynthesis. Identified transcripts and metabolites have built a foundation as new molecular resources that could facilitate future research on biosynthesis, regulation, and engineering of cerpegin or other important metabolites in such non-model plants.

Insights into Heterologous Biosynthesis of Arteannuin B and Artemisinin in Physcomitrella patens

: Metabolic engineering is an integrated bioengineering approach, which has made considerable progress in producing terpenoids in plants and fermentable hosts. Here, the full biosynthetic pathway of artemisinin, originating from Artemisia annua, was integrated into the moss Physcomitrella patens. Different combinations of the five artemisinin biosynthesis genes were ectopically expressed in P. patens to study biosynthesis pathway activity, but also to ensure survival of successful transformants. Transformation of the first pathway gene, ADS, into P. patens resulted in the accumulation of the expected metabolite, amorpha-4,11-diene, and also accumulation of a second product, arteannuin B. This demonstrates the presence of endogenous promiscuous enzyme activity, possibly cytochrome P450s, in P. patens. Introduction of three pathway genes, ADS-CYP71AV1-ADH1 or ADS-DBR2-ALDH1 both led to the accumulation of artemisinin, hinting at the presence of one or more endogenous enzymes in P. patens that can complement the partial pathways to full pathway activity. Transgenic P. patens lines containing the different gene combinations produce artemisinin in varying amounts. The pathway gene expression in the transgenic moss lines correlates well with the chemical profile of pathway products. Moreover, expression of the pathway genes resulted in lipid body formation in all transgenic moss lines, suggesting that these may have a function in sequestration of heterologous metabolites. This work thus provides novel insights into the metabolic response of P. patens and its complementation potential for A. annua artemisinin pathway genes. Identification of the related endogenous P. patens genes could contribute to a further successful metabolic engineering of artemisinin biosynthesis, as well as bioengineering of other high-value terpenoids in P. patens.

Engineering of CYP76AH15 can improve activity and specificity towards forskolin biosynthesis in yeast

Background

Forskolin is a high-value diterpenoid produced exclusively by the Lamiaceae plant Coleus forskohlii. Today forskolin is used pharmaceutically for its adenyl-cyclase activating properties. The limited availability of pure  forskolin is currently hindering its full utilization, thus a new environmentally friendly, scalable and sustainable strategy is needed for forskolin production. Recently, the entire biosynthetic pathway leading to forskolin was elucidated. The key steps of the pathway are catalyzed by cytochrome P450 enzymes (CYPs), which have been shown to be the limiting steps of the pathway. Here we study whether protein engineering of the substrate recognition sites (SRSs) of CYPs can improve their efficiency towards forskolin biosynthesis in yeast.

Results

As a proof of concept, we engineered the enzyme responsible for the first putative oxygenation step of the forskolin pathway: the conversion of 13R-manoyl oxide to 11-oxo-13R-manoyl oxide, catalyzed by the CYP76AH15. Four CYP76AH15 variants—engineered in the SRS regions—yielded at least a twofold increase of 11-oxo-13R-manoyl oxide when expressed in yeast cells grown in microtiter plates. The highest titers (5.6-fold increase) were observed with the variant A99I, mutated in the SRS1 region. Double or triple CYP76AH15 mutant variants resulted in additional enzymes with optimized performances. Moreover, in planta CYP76AH15 can synthesize ferruginol from miltiradiene. In this work, we showed that the mutants affecting 11-oxo-13R-manoyl oxide synthesis, do not affect ferruginol production, and vice versa.

The best performing variant, A99I, was utilized to reconstruct the forskolin biosynthetic pathway in yeast cells. Although these strains showed increased 11-oxo-manoyl oxide production and higher accumulation of other pathway intermediates compared to the native CYP76AH15, lower production of forskolin was observed.

Conclusions

As demonstrated for CYP76AH15, site-directed mutagenesis of SRS regions of plant CYPs may be an efficient and targeted approach to increase the performance of these enzymes. Although in this work we have managed to achieve higher efficiency and specificity of the first CYP of the pathway, further work is necessary in order to increase the overall production of forskolin in yeast cells.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1027-3) contains supplementary material, which is available to authorized users.

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

The production of coumarins and furanocoumarins (FCs) in higher plants is widely considered a model illustration of the adaptation of plants to their environment. In this report, we show that the multiplication of cytochrome P450 variants within the CYP71AZ subfamily has contributed to the diversification of these molecules. Multiple copies of genes encoding this enzyme family are found in Apiaceae, and their phylogenetic analysis suggests that they have different functions within these plants. CYP71AZ1 from Ammi majus and CYP71AZ3, 4, and 6 from Pastinaca sativa were functionally characterized. While CYP71AZ3 merely hydroxylated esculetin, the other enzymes accepted both simple coumarins and FCs. Superimposing in silico models of these enzymes led to the identification of different conformations of three regions in the enzyme active site. These sequences were subsequently utilized to mutate CYP71AZ4 to resemble CYP71AZ3. The swapping of these regions lead to significantly modified substrate specificity. Simultaneous mutations of all three regions shifted the specificity of CYP71AZ4 to that of CYP71AZ3, exclusively accepting esculetin. This approach may explain the evolution of this cytochrome P450 family regarding the appearance of FCs in parsnip and possibly in the Apiaceae.

In vitro characterization of CYP102G4 from Streptomyces cattleya: A self-sufficient P450 naturally producing indigo

Self-sufficient CYP102As possess outstanding hydroxylating activity to fatty acids such as myristic acid. Other CYP102 subfamily members share substrate specificity of CYP102As, but, occasionally, unusual characteristics of its own subfamily have been found. In this study, only one self-sufficient cytochrome P450 from Streptomyces cattleya was renamed from CYP102A_scat to CYP102G4, purified and characterized. UV-Vis spectrometry pattern, FAD/FMN analysis, and protein sequence comparison among CYP102s have shown that CYP102 from Streptomyces cattleya belongs to CYP102G subfamily. It showed hydroxylation activity toward fatty acids generating ω-1, ω-2, and ω-3-hydroxyfatty acids, which is similar to the general substrate specificity of CYP102 family. Unexpectedly, however, expression of CYP102G4 showed indigo production in LB medium batch flask culture, and high catalytic activity (k_cat/K_m) for indole was measured as 6.14±0.10min^-1mM^-1. Besides indole, CYP102G4 was able to hydroxylate aromatic compounds such as flavone, benzophenone, and chloroindoles. Homology model has shown such ability to accept aromatic compounds is due to its bigger active site cavity. Unlike other CYP102s, CYP102G4 did not have biased cofactor dependency, which was possibly determined by difference in NAD(P)H binding residues (Ala984, Val990, and Tyr1064) compared to CYP102A1 (Arg966, Lys972 and Trp1046). Overall, a self-sufficient CYP within CYP102G subfamily was characterized using purified enzymes, which appears to possess unique properties such as an only prokaryotic CYP naturally producing indigo.

Localization and in-Vivo Characterization of Thapsia garganica CYP76AE2 Indicates a Role in Thapsigargin Biosynthesis

The Mediterranean plant Thapsia garganica (dicot, Apiaceae), also known as deadly carrot, produces the highly toxic compound thapsigargin. This compound is a potent inhibitor of the sarcoplasmic-endoplasmic reticulum Ca2+-ATPase calcium pump in mammals and is of industrial importance as the active moiety of the anticancer drug mipsagargin, currently in clinical trials. Knowledge of thapsigargin in planta storage and biosynthesis has been limited. Here, we present the putative second step in thapsigargin biosynthesis, by showing that the cytochrome P450 TgCYP76AE2, transiently expressed in Nicotiana benthamiana, converts epikunzeaol into epidihydrocostunolide. Furthermore, we show that thapsigargin is likely to be stored in secretory ducts in the roots. Transcripts from TgTPS2 (epikunzeaol synthase) and TgCYP76AE2 in roots were found only in the epithelial cells lining these secretory ducts. This emphasizes the involvement of these cells in the biosynthesis of thapsigargin. This study paves the way for further studies of thapsigargin biosynthesis.

A bacterial artificial chromosome (BAC) genomic approach reveals partial clustering of the furanocoumarin pathway genes in parsnip

Furanocoumarins are specialized metabolites that are involved in the defense of plants against phytophagous insects. The molecular and functional characterization of the genes involved in their biosynthetic pathway is only partially complete. Many recent reports have described gene clusters responsible for the biosynthesis of specialized metabolites in plants. To investigate possible co-localization of the genes involved in the furanocoumarin pathway, we sequenced parsnip BAC clones spanning two different gene loci. We found that two genes previously identified in this pathway, CYP71AJ3 and CYP71AJ4, were located on the same BAC, whereas a third gene, PsPT1, belonged to a different BAC clone. Chromosome mapping using fluorescence in situ hybridization (FISH) indicated that PsPT1 and the CYP71AJ3-CYP71AJ4 clusters are located on two different chromosomes. Sequencing the BAC clone harboring PsPT1 led to the identification of a gene encoding an Fe(II) α-ketoglutarate-dependent dioxygenase (PsDIOX) situated in the neighborhood of PsPT1 and confirmed the occurrence of a second gene cluster involved in the furanocoumarin pathway. This enzyme metabolizes p-coumaroyl CoA, leading exclusively to the synthesis of umbelliferone, an important intermediate compound in furanocoumarin synthesis. This work provides an insight into the genomic organization of genes from the furanocoumarin biosynthesis pathway organized in more than one gene cluster. It also confirms that the screening of a genomic library and the sequencing of BAC clones represent a valuable tool to identify genes involved in biosynthetic pathways dedicated to specialized metabolite synthesis.

Semirational Protein Engineering of CYP153AM.aq. -CPRBM3 for Efficient Terminal Hydroxylation of Short- to Long-Chain Fatty Acids

The regioselective terminal hydroxylation of alkanes and fatty acids is of great interest in a variety of industrial applications, such as in cosmetics, in fine chemicals, and in the fragrance industry. The chemically challenging activation and oxidation of non-activated C-H bonds can be achieved with cytochrome P450 enzymes. CYP153AM.aq. -CPRBM3 is an artificial fusion construct consisting of the heme domain from Marinobacter aquaeolei and the reductase domain of CYP102A1 from Bacillus megaterium. It has the ability to hydroxylate medium- and long-chain fatty acids selectively at their terminal positions. However, the activity of this interesting P450 construct needs to be improved for applications in industrial processes. For this purpose, the design of mutant libraries including two consecutive steps of mutagenesis is demonstrated. Targeted positions and residues chosen for substitution were based on semi-rational protein design after creation of a homology model of the heme domain of CYP153AM.aq. , sequence alignments, and docking studies. Site-directed mutagenesis was the preferred method employed to address positions within the binding pocket, whereas diversity was created with the aid of a degenerate codon for amino acids located at the substrate entrance channel. Combining the successful variants led to the identification of a double variant-G307A/S233G-that showed alterations of one position within the binding pocket and one position located in the substrate access channel. This double variant showed twofold increased activity relative to the wild type for the terminal hydroxylation of medium-chain-length fatty acids. This variant furthermore showed improved activity towards short- and long-chain fatty acids and enhanced stability in the presence of higher concentrations of fatty acids.

Molecular evolution of parsnip (Pastinaca sativa) membrane-bound prenyltransferases for linear and/or angular furanocoumarin biosynthesis

In Apiaceae, furanocoumarins (FCs) are plant defence compounds that are present as linear or angular isomers. Angular isomers appeared during plant evolution as a protective response to herbivores that are resistant to linear molecules. Isomeric biosynthesis occurs through prenylation at the C6 or C8 position of umbelliferone. Here, we report cloning and functional characterization of two different prenyltransferases, Pastinaca sativa prenyltransferase 1 and 2 (PsPT1 and PsPT2), that are involved in these crucial reactions. Both enzymes are targeted to plastids and synthesize osthenol and demethylsuberosin (DMS) using exclusively umbelliferone and dimethylallylpyrophosphate (DMAPP) as substrates. Enzymatic characterization using heterologously expressed proteins demonstrated that PsPT1 is specialized for the synthesis of the linear form, demethylsuberosin, whereas PsPT2 more efficiently catalyses the synthesis of its angular counterpart, osthenol. These results are the first example of a complementary prenyltransferase pair from a single plant species that is involved in synthesizing defensive compounds. This study also provides a better understanding of the molecular mechanisms governing the angular FC biosynthetic pathway in apiaceous plants, which involves two paralogous enzymes that share the same phylogenetic origin.