Cytochromes p450

Overlooked and misunderstood: can glutathione conjugates be clues to understanding plant glutathione transferases?

Plant glutathione transferases (GSTs) constitute a large and diverse family of enzymes that are involved in plant stress response, metabolism and defence, yet their physiological functions remain largely elusive. Consistent with the traditional view on GSTs across organisms as detoxification enzymes, in vitro most plant GSTs catalyse glutathionylation, conjugation of the tripeptide glutathione (GSH; γ-Glu-Cys-Gly) onto reactive molecules. However, when it comes to elucidating GST functions, it remains a key challenge that the endogenous plant glutathione conjugates (GS-conjugates) that would result from such glutathionylation reactions are rarely reported. Furthermore, GSTs often display high substrate promiscuity, and their proposed substrates are prone to spontaneous chemical reactions with GSH; hence, single-gene knockouts rarely provide clear chemotypes or phenotypes. In a few cases, GS-conjugates are demonstrated to be biosynthetic intermediates that are rapidly further metabolized towards a pathway end product, explaining their low abundance and rare detection. In this review, we summarize the current knowledge of plant GST functions and how and possibly why evolution has resulted in a broad and extensive expansion of the plant GST family. Finally, we demonstrate that endogenous GS-conjugates are more prevalent in plants than assumed and suggest they are overlooked as clues towards the identification of plant GST functions. This article is part of the theme issue 'The evolution of plant metabolism'.

Contrasting defence mechanisms against spider mite infestation in cyanogenic and non-cyanogenic legumes

Understanding the complex interactions between plants and herbivores is essential for improving crop resistance. Aiming to expand the role of cyanogenesis in plant defence, we investigated the response of the cyanogenic Phaseolus lunatus (lima bean) and the non-cyanogenic Phaseolus vulgaris (common bean) to Tetranychus urticae (spider mite) infestation. Despite mite infesting both legumes, leaf damage infringed by this feeder was reduced in lima bean. Comparative transcriptome analyses revealed that both species exhibited substantial metabolic and transcriptional changes upon infestation, although alterations in P. lunatus were significantly more pronounced. Specific differences in amino acid homeostasis and key genes associated with the cyanogenic pathway were observed in these species, as well as the upregulation of the mandelonitrile lyase gene (PlMNL1) following T. urticae feeding. Concomitantly, the PIMNL1 activity increased. Lima bean plants also displayed an induction of β-cyanoalanine synthase (PlCYSC1), a key enzyme for cyanide detoxification, suggesting an internal regulatory mechanism to manage the toxicity of their defence responses. These findings contribute to our understanding of the legume-herbivore interactions and underscore the potential role of cyanogenesis in the elaboration of specific defensive responses, even within the same genus, which may reflect distinctive evolutionary adaptations or varying metabolic capabilities between species.

Dual Regulation of Cytochrome P450 Gene Expression by Two Distinct Small RNAs, a Novel tasiRNA and miRNA, in Marchantia polymorpha

The miR390-derived TAS3 trans-acting short-interfering RNAs (tasiRNAs) module represents a conserved RNA silencing pathway in the plant kingdom; however, its characterization in the bryophyte Marchantia polymorpha is limited. This study elucidated that MpDCL4 processes MpTAS3 double-stranded RNA (dsRNA) to generate tasiRNAs, primarily from the 5'- and 3'-ends of dsRNA. Notably, we discovered a novel tasiRNA, tasi78A, which can negatively regulate a cytochrome P450 gene, MpCYP78A101. Additionally, tasi78A was abundant in MpAGO1, and transient expression assays underscored the role of tasi78A in repressing MpCYP78A101. A microRNA, miR11700, also regulates MpCYP78A101 expression. This coordinate regulation suggests a role in modulating auxin signaling at apical notches of gemma, influencing the growth and sexual organ development of M. polymorpha and emphasizing the significance of RNA silencing in MpCYP78A101 regulation. However, phylogenetic analysis identified another paralog of the CYP78 family, Mp1g14150, which may have a redundant role with MpCYP78A101, explaining the absence of noticeable morphological changes in loss-of-function plants. Taken together, our findings provide new insights into the combined regulatory roles of miR390/MpTAS3/miR11700 in controlling MpCYP78A101 and expand our knowledge about the biogenesis and regulation of tasiRNAs in M. polymorpha.

Identification and Expression Analysis of Cytochrome P450 Genes Probably Involved in Triterpenoid Saponins Biosynthesis in Astragalus mongholicus

Cytochromes P450 (P450s) are one of the largest enzymatic protein families and play critical roles in the synthesis and metabolism of plant secondary metabolites. Astragaloside IV (AS-IV) is one of the primary active components in Astragalus herbs, exhibiting diverse biological activities and pharmacological effects. However, P450s involved in the astragaloside biosynthesis have not been systematically analyzed in Astragalus mongholicus (A. mongholicus). In this study, we identified 209 P450 genes from the genome of A. mongholicus (AmP450s), which were classified into nine clans and 47 families and performed a systematic overview of their physical and chemical properties, phylogeny, gene structures and conserved motifs. Weighted gene co-expression network analysis (WGCNA) revealed that AmP450s are critical in the astragaloside biosynthesis pathway. The expression levels of these AmP450s were verified by quantitative real-time PCR (qRT-PCR) analysis in the root, stem and leaf, showing that most AmP450s are abundant in the root. Additionally, the correlation analysis between gene expressions and AS-IV content showed that twelve AmP450s, especially CYP71A28, CYP71D16 and CYP72A69, may have significant potential in the biosynthesis of astragaloside. This study systematically investigates the P450s of A. mongholicus and offers valuable insights into further exploring the functions of CYP450s in the astragaloside biosynthesis pathway.

Unraveling the role of cytochrome P450 enzymes in oleanane triterpenoid biosynthesis in arjuna tree

Triterpenoids (C30-isoprenoids) represent a major group of natural products with various physiological functions in plants. Triterpenoids and their derivatives have medicinal uses owing to diverse bioactivities. Arjuna (Terminalia arjuna) tree bark accumulates highly oxygenated β-amyrin-derived oleanane triterpenoids (e.g., arjunic acid, arjungenin, and arjunolic acid) with cardioprotective roles. However, biosynthetic routes and enzymes remain poorly understood. We mined the arjuna transcriptome and conducted cytochrome P450 monooxygenase (P450) assays using Saccharomyces cerevisiae and Nicotiana benthamiana to identify six P450s and two P450 reductases for oxidative modifications of oleanane triterpenoids. P450 assays using oleananes revealed a greater substrate promiscuity of C-2α and C-23 hydroxylases/oxidases than C-28 oxidases. CYP716A233 and CYP716A432 catalyzed β-amyrin/erythrodiol C-28 oxidation to produce oleanolic acid. C-2α hydroxylases (CYP716C88 and CYP716C89) converted oleanolic acid and hederagenin to maslinic acid and arjunolic acid. CYP716C89 also hydroxylated erythrodiol and oleanolic aldehyde. However, CYP714E107a and CYP714E107b catalyzed oleanolic acid/maslinic acid/arjunic acid, C-23 hydroxylation to form hederagenin, arjunolic acid and arjungenin, and hederagenin C-23 oxidation to produce gypsogenic acid, but at a lower rate than oleanolic acid C-23 hydroxylation. Overall, P450 substrate selectivity suggested that C-28 oxidation is the first P450-catalyzed oxidative modification in the arjuna triterpenoid pathway. However, the pathway might branch thereafter through C-2α/C-23 hydroxylation of oleanolic acid. Taken together, these results provided new insights into substrate range of P450s and unraveled biosynthetic routes of triterpenoids in arjuna. Moreover, complete elucidation and reconstruction of arjunolic acid pathway in S. cerevisiae and N. benthamiana suggested the utility of arjuna P450s in heterologous production of cardioprotective compounds.

Enhancing drought tolerance in Malva parviflora plants through metabolic and genetic modulation using Beauveria bassiana inoculation

BACKGROUND

Enhancing crops' drought resilience is necessary to maintain productivity levels. Plants interact synergistically with microorganisms like Beauveria bassiana to improve drought tolerance. Therefore, the current study investigates the effects of biopriming with B. bassiana on drought tolerance in Malva parviflora plants grown under regular irrigation (90% water holding capacity (WHC)), mild (60% WHC), and severe drought stress (30% WHC).

RESULTS

The results showed that drought stress reduced the growth and physiological attributes of M. parviflora. However, those bioprimed with B. bassiana showed higher drought tolerance and enhanced growth, physiological, and biochemical parameters: drought stress enriched malondialdehyde and H2O2 contents. Conversely, exposure to B. bassiana reduced stress markers and significantly increased proline and ascorbic acid content under severe drought stress; it enhanced gibberellic acid and reduced ethylene. Bioprimed M. parviflora, under drought conditions, improved antioxidant enzymatic activity and the plant's nutritional status. Besides, ten Inter-Simple Sequence Repeat primers detected a 25% genetic variation between treatments. Genomic DNA template stability (GTS) decreased slightly and was more noticeable in response to drought stress; however, for drought-stressed plants, biopriming with B. bassiana retained the GTS.

CONCLUSION

Under drought conditions, biopriming with B. bassiana enhanced Malva's growth and nutritional value. This could attenuate photosynthetic alterations, up-regulate secondary metabolites, activate the antioxidant system, and maintain genome integrity.

Cytochrome b5 diversity in green lineages preceded the evolution of syringyl lignin biosynthesis

Lignin production marked a milestone in vascular plant evolution, and the emergence of syringyl (S) lignin is lineage specific. S-lignin biosynthesis in angiosperms, mediated by ferulate 5-hydroxylase (F5H, CYP84A1), has been considered a recent evolutionary event. F5H uniquely requires the cytochrome b5 protein CB5D as an obligatory redox partner for catalysis. However, it remains unclear how CB5D functionality originated and whether it coevolved with F5H. We reveal here the ancient evolution of CB5D-type function supporting F5H-catalyzed S-lignin biosynthesis. CB5D emerged in charophyte algae, the closest relatives of land plants, and is conserved and proliferated in embryophytes, especially in angiosperms, suggesting functional diversification of the CB5 family before terrestrialization. A sequence motif containing acidic amino residues in Helix 5 of the CB5 heme-binding domain contributes to the retention of CB5D function in land plants but not in algae. Notably, CB5s in the S-lignin-producing lycophyte Selaginella lack these residues, resulting in no CB5D-type function. An independently evolved S-lignin biosynthetic F5H (CYP788A1) in Selaginella relies on NADPH-dependent cytochrome P450 reductase as sole redox partner, distinct from angiosperms. These results suggest that angiosperm F5Hs coopted the ancient CB5D, forming a modern cytochrome P450 monooxygenase system for aromatic ring meta-hydroxylation, enabling the reemergence of S-lignin biosynthesis in angiosperms.

Decoding the Role of CYP450 Enzymes in Metabolism and Disease: A Comprehensive Review

Cytochrome P450 (CYP450) is a group of enzymes that play an essential role in Phase I metabolism, with 57 functional genes classified into 18 families in the human genome, of which the CYP1, CYP2, and CYP3 families are prominent. Beyond drug metabolism, CYP enzymes metabolize endogenous compounds such as lipids, proteins, and hormones to maintain physiological homeostasis. Thus, dysregulation of CYP450 enzymes can lead to different endocrine disorders. Moreover, CYP450 enzymes significantly contribute to fatty acid metabolism, cholesterol synthesis, and bile acid biosynthesis, impacting cellular physiology and disease pathogenesis. Their diverse functions emphasize their therapeutic potential in managing hypercholesterolemia and neurodegenerative diseases. Additionally, CYP450 enzymes are implicated in the onset and development of illnesses such as cancer, influencing chemotherapy outcomes. Assessment of CYP450 enzyme expression and activity aids in evaluating liver health state and differentiating between liver diseases, guiding therapeutic decisions, and optimizing drug efficacy. Understanding the roles of CYP450 enzymes and the clinical effect of their genetic polymorphisms is crucial for developing personalized therapeutic strategies and enhancing drug responses in diverse patient populations.

High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation

High-quality genome of rosemary (Salvia rosmarinus) represents a valuable resource and tool for understanding genome evolution and environmental adaptation as well as its genetic improvement. However, the existing rosemary genome did not provide insights into the relationship between antioxidant components and environmental adaptability. In this study, by employing Nanopore sequencing and Hi-C technologies, a total of 1.17 Gb (97.96%) genome sequences were mapped to 12 chromosomes with 46 121 protein-coding genes and 1265 non-coding RNA genes. Comparative genome analysis reveals that rosemary had a closely genetic relationship with Salvia splendens and Salvia miltiorrhiza, and it diverged from them approximately 33.7 million years ago (MYA), and one whole-genome duplication occurred around 28.3 MYA in rosemary genome. Among all identified rosemary genes, 1918 gene families were expanded, 35 of which are involved in the biosynthesis of antioxidant components. These expanded gene families enhance the ability of rosemary adaptation to adverse environments. Multi-omics (integrated transcriptome and metabolome) analysis showed the tissue-specific distribution of antioxidant components related to environmental adaptation. During the drought, heat and salt stress treatments, 36 genes in the biosynthesis pathways of carnosic acid, rosmarinic acid and flavonoids were up-regulated, illustrating the important role of these antioxidant components in responding to abiotic stresses by adjusting ROS homeostasis. Moreover, cooperating with the photosynthesis, substance and energy metabolism, protein and ion balance, the collaborative system maintained cell stability and improved the ability of rosemary against harsh environment. This study provides a genomic data platform for gene discovery and precision breeding in rosemary. Our results also provide new insights into the adaptive evolution of rosemary and the contribution of antioxidant components in resistance to harsh environments.

Two Cytochrome P450s, CYP709B1 and CYP704C1, Play Essential Roles in Metabolism-Based Multiple Herbicide Resistance in American Sloughgrass (Beckmannia syzigachne (Steud.) Fernald)

This study confirmed a field population of American sloughgrass (Beckmannia syzigachne (Steud.) Fernald) that developed simultaneously high levels of resistance (resistance index >10) to three divergent modes of action herbicides: fenoxaprop-P-ethyl, mesosulfuron-methyl, and isoproturon. The resistance phenotype observed in this population was not attributed to target-site alterations; rather, the resistant plants exhibited a significant increase in the activity of cytochrome P450s (P450s) and enhanced metabolism rates for all three herbicides. RNA sequencing revealed significant upregulation of two P450s, CYP709B1 and CYP704C1, in the resistant plants both before and after herbicide treatments. Molecular docking predicted that the homology models of these P450s should exhibit a binding affinity for a range of herbicides. The heterologous expression of the identified P450s in yeast cells indicated improved growth in the presence of all three of the aforementioned herbicides. Collectively, the increased expression of CYP709B1 and CYP704C1 likely contributed to the P450s-mediated enhanced metabolism, thereby conferring multiple herbicide resistance in B. syzigachne.

Regulation by distinct MYB transcription factors defines the roles of OsCYP86A9 in anther development and root suberin deposition

Cytochrome P450 proteins (CYPs) play critical roles in plant development and adaptation to fluctuating environments. Previous reports have shown that CYP86A proteins are involved in the biosynthesis of suberin and cutin in Arabidopsis. However, the functions of these proteins in rice remain obscure. In this study, a rice mutant with incomplete male sterility was identified. Cytological analyses revealed that this mutant was defective in anther development. Cloning of the mutant gene indicated that the responsible mutation was on OsCYP86A9. OsMYB80 is a core transcription factor in the regulation of rice anther development. The expression of OsCYP86A9 was abolished in the anther of osmyb80 mutant. In vivo and in vitro experiments showed that OsMYB80 binds to the MYB-binding motifs in OsCYP86A9 promoter region and regulates its expression. Furthermore, the oscyp86a9 mutant exhibited an impaired suberin deposition in the root, and was more susceptible to drought stress. Interestingly, genetic and biochemical analyses revealed that OsCYP86A9 expression was regulated in the root by certain MYB transcription factors other than OsMYB80. Moreover, mutations in the MYB genes that regulate OsCYP86A9 expression in the root did not impair the male fertility of the plant. Taken together, these findings revealed the critical roles of OsCYP86A9 in plant development and proposed that OsCYP86A9 functions in anther development and root suberin formation via two distinct tissue-specific regulatory pathways.

The Dark Side of the pollen: BSA-seq identified genomic regions linked to male sterility in globe artichoke

Globe artichoke (Cynara cardunculus var. scolymus; 2n = 2x = 34) is a food crop consumed for its immature flower heads. Traditionally, globe artichoke varietal types are vegetatively propagated. However, seed propagation makes it possible to treat the crop as annual, increasing field uniformity and reducing farmers costs, as well as pathogens diffusion. Despite globe artichoke's significant agricultural value and the critical role of heterosis in the development of superior varieties, the production of hybrids remains challenging without a reliable system for large-scale industrial seed production. Male sterility (MS) presents a promising avenue for overcoming these challenges by simplifying the hybridization process and enabling cost-effective seed production. However, within the Cynara genus, genic male sterility has been linked to three recessive loci in globe artichoke, with no definitive genetic mechanism elucidated to date. A 250 offsprings F2 population, derived from a cross between a MS globe artichoke and a male fertile (MF) cultivated cardoon (C. cardunculus var. altilis) and fitting a monogenic segregation model (3:1), was analyzed through BSA-seq, aiming at the identification of genomic regions/genes affecting male sterility. Four QTL regions were identified on chromosomes 4, 12, and 14. By analyzing the sequence around the highest pick on chromosome 14, a cytochrome P450 (CYP703A2) was identified, carrying a deleterious substitution (R/Q) fixed in the male sterile parent. A single dCAPS marker was developed around this SNP, allowing the discrimination between MS and MF genotypes within the population, suitable for applications in plant breeding programs. A 3D model of the protein was generated by homology modeling, revealing that the mutated amino acid is part of a highly conserved motif crucial for protein folding.

Discovery of active mouse, plant and fungal cytochrome P450s in endogenous proteomes and upon expression in planta

Eukaryotes produce a large number of cytochrome P450s that mediate the synthesis and degradation of diverse endogenous and exogenous metabolites. Yet, most of these P450s are uncharacterized and global tools to study these challenging, membrane-resident enzymes remain to be exploited. Here, we applied activity profiling of plant, mouse and fungal P450s with chemical probes that become reactive when oxidized by P450 enzymes. Identification by mass spectrometry revealed labeling of a wide range of active P450s, including six plant P450s, 40 mouse P450s and 13 P450s of the fungal wheat pathogen Zymoseptoria tritici. We next used transient expression of GFP-tagged P450s by agroinfiltration to show ER-targeting and NADPH-dependent, activity-based labeling of plant, mouse and fungal P450s. Both global profiling and transient expression can be used to detect a broad range of active P450s to study e.g. their regulation and discover selective inhibitors.

Overexpressing CYP81D11 enhances 2,4,6-trinitrotoluene tolerance and removal efficiency in Arabidopsis

Phytoremediation is a promising technology for removing the high-toxic explosive 2,4,6-trinitrotoluene (TNT) pollutant from the environment. Mining dominant genes is the key research direction of this technology. Most previous studies have focused on the detoxification of TNT rather than plants' TNT tolerance. Here, we conducted a transcriptomic analysis of wild type Arabidopsis plants under TNT stress and found that the Arabidopsis cytochrome P450 gene CYP81D11 was significantly induced in TNT-treated plants. Under TNT stress, the root length was approximately 1.4 times longer in CYP81D11-overexpressing transgenic plants than in wild type plants. The half-removal time for TNT was much shorter in CYP81D11-overexpressing transgenic plants (1.1 days) than in wild type plants (t1/2 = 2.2 day). In addition, metabolic analysis showed no difference in metabolites in transgenic plants compared to wild type plants. These results suggest that the high TNT uptake rates of CYP81D11-overexpressing transgenic plants were most likely due to increased tolerance and biomass rather than TNT degradation. However, CYP81D11-overexpressing plants were not more tolerant to osmotic stresses, such as salt or drought. Taken together, our results indicate that CYP81D11 is a promising target for producing bioengineered plants with high TNT removing capability.

Hordeum vulgare CYP76M57 catalyzes C2 shortening of tryptophan side chain by C-N bond rearrangement in gramine biosynthesis

The indole alkaloid gramine, 3-(dimethylaminomethyl)indole, is a defensive specialized metabolite found in some barley cultivars. In its biosynthetic process, the tryptophan (Trp) side chain is shortened by two carbon atoms to produce 3-(aminomethyl)indole (AMI), which is then methylated by N-methyltransferase (HvNMT) to produce gramine. Although side chain shortening is one of the crucial scaffold formation steps of alkaloids originating from aromatic amino acids, the gene and enzyme involved in the Trp-AMI conversion reactions are unknown. In this study, through RNA-seq analysis, 35 transcripts were shown to correlate with gramine production; among them, an uncharacterized cytochrome P450 (CYP) gene, CYP76M57, and HvNMT were identified as candidate genes for gramine production. Transgenic Arabidopsis thaliana and rice overexpressing CYP and HvNMT accumulate AMI, N-methyl-AMI, and gramine. CYP76M57, heterologously expressed in Pichia pastoris, was able to act on Trp to produce AMI. Furthermore, the amino group nitrogen of Trp was retained during the CYP76M57-catalyzed reaction, indicating that the C2 shortening of Trp proceeds with an unprecedented biosynthetic process, the removal of the carboxyl group and Cα and the rearrangement of the nitrogen atom to Cβ. In some gramine-non-accumulating barley cultivars, arginine 104 in CYP76M57 is replaced by threonine, which abolished the catalytic activity of CYP76M57 to convert Trp into AMI. These results uncovered the missing committed enzyme of gramine biosynthesis in barley and contribute to the elucidation of the potential functions of CYPs in plants and undiscovered specialized pathways.

Exploration and utilization of novel aldoxime, nitrile, and nitro compounds metabolizing enzymes from plants and arthropods

Aldoxime (R1R2C=NOH) and nitrile (R-C≡N) are nitrogen-containing compounds that are found in species representing all kingdoms of life. The enzymes discovered from the microbial "aldoxime-nitrile" pathway (aldoxime dehydratase, nitrile hydratase, amidase, and nitrilase) have been thoroughly studied because of their industrial importance. Although plants utilize cytochrome P450 monooxygenases to produce aldoxime and nitrile, many biosynthetic pathways are yet to be studied. Cyanogenic millipedes accumulate various nitrile compounds, such as mandelonitrile. However, no such aldoxime- and nitrile-metabolizing enzymes have been identified in millipedes. Here, I review the exploration of novel enzymes from plants and millipedes with characteristics distinct from those of microbial enzymes, the catalysis of industrially useful reactions, and applications of these enzymes for nitrile compound production.

Genome-wide identification and expression analyses of CYP450 genes in sweet potato (Ipomoea batatas L.)

BACKGROUND

Cytochrome P450 monooxygenases (CYP450s) play a crucial role in various biochemical reactions involved in the synthesis of antioxidants, pigments, structural polymers, and defense-related compounds in plants. As sweet potato (Ipomoea batatas L.) holds significant economic importance, a comprehensive analysis of CYP450 genes in this plant species can offer valuable insights into the evolutionary relationships and functional characteristics of these genes.

RESULTS

In this study, we successfully identified and categorized 95 CYP450 genes from the sweet potato genome into 5 families and 31 subfamilies. The predicted subcellular localization results indicate that CYP450s are distributed in the cell membrane system. The promoter region of the IbCYP450 genes contains various cis-acting elements related to plant hormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) have been identified in the IbCYP450 family proteins, with 5 genes lacking introns and only one exon. We observed extensive duplication events within the CYP450 gene family, which may account for its expansion. The gene duplication analysis results showed the presence of 15 pairs of genes with tandem repeats. Interaction network analysis reveals that IbCYP450 families can interact with multiple target genes and there are protein-protein interactions within the family. Transcription factor interaction analysis suggests that IbCYP450 families interact with multiple transcription factors. Furthermore, gene expression analysis revealed tissue-specific expression patterns of CYP450 genes in sweet potatoes, as well as their response to abiotic stress and plant hormones. Notably, quantitative real-time polymerase chain reaction (qRT‒PCR) analysis indicated the involvement of CYP450 genes in the defense response against nonbiological stresses in sweet potatoes.

CONCLUSIONS

These findings provide a foundation for further investigations aiming to elucidate the biological functions of CYP450 genes in sweet potatoes.

Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis

As an economically important tree, Gleditsia sinensis Lam. is widely planted. A lack of background genetic information on G. sinensis hinders molecular breeding. Based on PacBio single-molecule real-time (SMRT) sequencing and analysis of G. sinensis, a total of 95,183 non-redundant transcript sequences were obtained, of which 93,668 contained complete open reading frames (ORFs), 2,858 were long non-coding RNAs (LncRNAs) and 18,855 alternative splicing (AS) events were identified. Genes orthologous to different Gleditsia species pairs were identified, stress-related genes had been positively selected during the evolution. AGA, AGG, and CCA were identified as the universal optimal codon in the genus of Gleditsia. EIF5A was selected as a suitable fluorescent quantitative reference gene. 315 Cytochrome P450 monooxygenases (CYP450s) and 147 uridine diphosphate (UDP)-glycosyltransferases (UGTs) were recognized through the PacBio SMRT transcriptome. Randomized selection of GsIAA14 for cloning verified the reliability of the PacBio SMRT transcriptome assembly sequence. In conclusion, the research data lay the foundation for further analysis of the evolutionary mechanism and molecular breeding of Gleditsia.

Distinct heat response molecular mechanisms emerge in cassava vasculature compared to leaf mesophyll tissue under high temperature stress

With growing concerns over global warming, cultivating heat-tolerant crops has become paramount to prepare for the anticipated warmer climate. Cassava (Manihot esculenta Crantz), a vital tropical crop, demonstrates exceptional growth and productivity under high-temperature (HT) conditions. Yet, studies elucidating HT resistance mechanisms in cassava, particularly within vascular tissues, are rare. We dissected the leaf mid-vein from leaf, and did the comparative transcriptome profiling between mid-vein and leaf to figure out the cassava vasculature HT resistance molecular mechanism. Anatomical microscopy revealed that cassava leaf veins predominantly consisted of vasculature. A thermal imaging analysis indicated that cassava experienced elevated temperatures, coinciding with a reduction in photosynthesis. Transcriptome sequencing produced clean reads in total of 89.17G. Using Venn enrichment, there were 65 differentially expressed genes (DEGs) and 93 DEGs had been found highly specifically expressed in leaf and mid-vein. Further investigation disclosed that leaves enhanced pyruvate synthesis as a strategy to withstand high temperatures, while mid-veins fortified themselves by bolstering lignin synthesis by comprehensive GO and KEGG analysis of DEGs. The identified genes in these metabolic pathways were corroborated through quantity PCR (QPCR), with results aligning with the transcriptomic data. To verify the expression localization of DEGs, we used in situ hybridization experiments to identify the expression of MeCCoAMT(caffeoyl-coenzyme A-3-O-methyltransferase) in the lignin synthesis pathway in cassava leaf veins xylem. These findings unravel the disparate thermotolerance mechanisms exhibited by cassava leaves and mid-veins, offering insights that could potentially inform strategies for enhancing thermotolerance in other crops.

Comparative Transcriptome Combined with Morphophysiological Analyses Revealed Carotenoid Biosynthesis for Differential Chilling Tolerance in Two Contrasting Rice (Oryza sativa L.) Genotypes

Early spring cold spells can lead to leaf chlorosis during the rice seedling greening process. However, the physiological and molecular mechanisms underlying the rice greening process under low-temperature conditions remain unknown. In this study, comparative transcriptome and morphophysiological analyses were performed to investigate the mechanisms mediating the responses of the Koshihikari (Kos) and Kasalath (Kas) rice cultivars to chilling stress. According to their growth-related traits, electrolyte leakage, and chlorophyll fluorescence parameters, Kos was more tolerant to low-temperature stress than Kas. Moreover, chloroplast morphology was more normal (e.g., oval) in Kos than in Kas at 17 °C. The comparative transcriptome analysis revealed 610 up-regulated differentially expressed genes that were common to all four comparisons. Furthermore, carotenoid biosynthesis was identified as a critical pathway for the Kos response to chilling stress. The genes in the carotenoid biosynthesis pathway were expressed at higher levels in Kos than in Kas at 17 °C, which was in accordance with the higher leaf carotenoid content in Kos than in Kas. The lycopene β-cyclase and lycopene ε-cyclase activities increased more in Kos than in Kas. Additionally, the increases in the violaxanthin de-epoxidase and carotenoid hydroxylase activities in Kos seedlings resulted in the accumulation of zeaxanthin and lutein and mitigated the effects of chilling stress on chloroplasts. These findings have clarified the molecular mechanisms underlying the chilling tolerance of rice seedlings during the greening process.

Oxidation of Naringenin, Apigenin, and Genistein by Human Family 1 Cytochrome P450 Enzymes and Comparison of Interaction of Apigenin with Human P450 1B1.1 and Scutellaria P450 82D.1

Naringenin, an initial synthesized flavanone in various plant species, is further utilized for production of many biologically active flavonoids, e.g., apigenin, eriodictyol, and genistein, by various plant enzymes including cytochrome P450s (P450s or CYPs). We examined how these flavonoids are oxidized by human P450 family 1 and 2A enzymes. Naringenin was principally oxidized at the 3'-position to form eriodictyol by CYP1 enzymes more efficiently than by CYP2A enzymes, and the resulting eriodictyol was further oxidized to two penta-hydroxylated products. In contrast to plant P450 enzymes, these human P450s did not mediate the desaturation of naringenin and eriodictyol to give apigenin and luteolin, respectively. Apigenin was oxidized at the C3' and C6 positions to form luteolin and scutellarein by these P450s. CYP1B1.1 and 1B1.3 had high activities in apigenin 6-hydroxylation with a homotropic cooperative manner, as has been observed previously in chrysin 6-hydroxylation (Nagayoshi et al., Chem. Res. Toxicol. 2019, 32, 1268-1280). Molecular docking analysis suggested that CYP1B1 had two apigenin binding sites and showed similarities in substrate recognition sites to plant CYP82D.1, one of the enzymes in catalyzing apigenin and chrysin 6-hydroxylations in Scutellaria baicalensis. The present results suggest that human CYP1 enzymes and CYP2A13 in some reactions have important roles in the oxidation of naringenin, eriodictyol, apigenin, and genistein and that human CYP1B1 and Scutellaria CYP82D.1 have similarities in their SRS regions, catalyzing 6-hydroxylation of both apigenin and chrysin.

CYP72D19 from Tripterygium wilfordii catalyzes C-2 hydroxylation of abietane-type diterpenoids

KEY MESSAGE

CYP72D19, the first functional gene of the CYP72D subfamily, catalyzes the C-2 hydroxylation of abietane-type diterpenoids. The abietane-type diterpenoids, e.g., triptolide, tripdiolide, and 2-epitripdiolide, are the main natural products for the anti-tumor, anti-inflammatory, and immunosuppressive activities of Tripterygium wilfordii, while their biosynthetic pathways are not resolved. Here, we cloned and characterized the CYP72D19-catalyzed C-2 hydroxylation of dehydroabietic acid, a compound that has been proven to be a biosynthetic intermediate in triptolide biosynthesis. Through molecular docking and site-directed mutagenesis, L386, L387, and I493 near the active pocket were found to have an important effect on the enzyme activity, which also indicates that steric hindrance of residues plays an important role in function. In addition, CYP72D19 also catalyzed a variety of abietane-type diterpenoids with benzene ring, presumably because the benzene ring of the substrate molecule stabilized the C-ring, allowing the protein and the substrate to form a relatively stable spatial structure. This is the first demonstration of CYP72D subfamily gene function. Our research provides important genetic elements for the structural modification of active ingredients and the heterologous production of other 2-hydroxyl abietane-type natural products.

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Enhanced production of difficult-to-express proteins through knocking down rnpA gene expression

Escherichia coli has been employed as a workhorse for the efficient production of recombinant proteins. However, some proteins were found to be difficult to produce in E. coli. The stability of mRNA has been considered as one of the important factors affecting recombinant protein production. Here we report a generally applicable and simple strategy for enhancing mRNA stability, and consequently improving recombinant protein production in E. coli. RNase P, a ribozyme comprising an RNA subunit (RnpB) and a protein subunit (RnpA), is involved in tRNA maturation. Based on the finding that purified RnpA can digest rRNA and mRNA in vitro, it was reasoned that knocking down the level of RnpA might enhance recombinant protein production. For this, the synthetic small regulatory RNA-based knockdown system was applied to reduce the expression level of RnpA. The developed RnpA knockdown system allowed successful overexpression of 23 different recombinant proteins of various origins and sizes, including Cas9 protein, antibody fragment, and spider silk protein. Notably, a 284.9-kDa ultra-high molecular weight, highly repetitive glycine-rich spider silk protein, which is one of the most difficult proteins to produce, could be produced to 1.38 g L-1 , about two-fold higher than the highest value previously achieved, by a fed-batch culture of recombinant E. coli strain employing the RnpA knockdown system. The RnpA-knockdown strategy reported here will be generally useful for the production of recombinant proteins including those that have been difficult to produce.

Combinatorial engineering of betalain biosynthesis pathway in yeast Saccharomyces cerevisiae

BACKGROUND

Betalains, comprising red-violet betacyanins and yellow-orange betaxanthins, are the hydrophilic vacuolar pigments that provide bright coloration to roots, fruits, and flowers of plants of the Caryophyllales order. Betanin extracted from red beets is permitted quantum satis as a natural red food colorant (E162). Due to antioxidant activity, betanin has potential health benefits.

RESULTS

We applied combinatorial engineering to find the optimal combination of a dozen tyrosine hydroxylase (TyH) and 4,5-dopa-estradiol-dioxygenase (DOD) variants. The best-engineered Saccharomyces cerevisiae strains produced over six-fold higher betaxanthins than previously reported. By genome-resequencing of these strains, we found out that two copies of DOD enzyme from Bougainvillea glabra together with TyH enzymes from Abronia nealleyi, Acleisanthes obtusa, and Cleretum bellidiforme were present in the three high-betaxanthin-producing isolates. Next, we expressed four variants of glucosyltransferases from Beta vulgaris for betanin biosynthesis. The highest titer of betanin (30.8 ± 0.14 mg/L after 48 h from 20 g/L glucose) was obtained when completing the biosynthesis pathway with UGT73A36 glucosyltransferase from Beta vulgaris. Finally, we investigated betalain transport in CEN.PK and S288C strains of Saccharomyces cerevisiae and identified a possible role of transporter genes QDR2 and APL1 in betanin transport.

CONCLUSIONS

This study shows the potential of combinatorial engineering of yeast cell factories for the biotechnological production of betanin.

Beauveria bassiana rewires molecular mechanisms related to growth and defense in tomato

Plant roots can exploit beneficial associations with soil-inhabiting microbes, promoting growth and expanding the immune capacity of the host plant. In this work, we aimed to provide new information on changes occurring in tomato interacting with the beneficial fungus Beauveria bassiana. The tomato leaf proteome revealed perturbed molecular pathways during the establishment of the plant-fungus relationship. In the early stages of colonization (5-7 d), proteins related to defense responses to the fungus were down-regulated and proteins related to calcium transport were up-regulated. At later time points (12-19 d after colonization), up-regulation of molecular pathways linked to protein/amino acid turnover and to biosynthesis of energy compounds suggests beneficial interaction enhancing plant growth and development. At the later stage, the profile of leaf hormones and related compounds was also investigated, highlighting up-regulation of those related to plant growth and defense. Finally, B. bassiana colonization was found to improve plant resistance to Botrytis cinerea, impacting plant oxidative damage. Overall, our findings further expand current knowledge on the possible mechanisms underlying the beneficial role of B. bassiana in tomato plants.

Decreased metabolic diversity in common beans associated with domestication revealed by untargeted metabolomics, information theory, and molecular networking

The process of crop domestication leads to a dramatic reduction in the gene expression associated with metabolic diversity. Genes involved in specialized metabolism appear to be particularly affected. Although there is ample evidence of these effects at the genetic level, a reduction in diversity at the metabolite level has been taken for granted despite having never been adequately accessed and quantified. Here we leveraged the high coverage of ultra high performance liquid chromatography-high-resolution mass spectrometry based metabolomics to investigate the metabolic diversity in the common bean (Phaseolus vulgaris). Information theory highlights a shift towards lower metabolic diversity and specialization when comparing wild and domesticated bean accessions. Moreover, molecular networking approaches facilitated a broader metabolite annotation than achieved to date, and its integration with gene expression data uncovers a metabolic shift from specialized metabolism towards central metabolism upon domestication of this crop.

Genome-Wide Identification and Expression Profiling of Cytochrome P450 Monooxygenase Superfamily in Foxtail Millet

The cytochrome P450 monooxygenases (CYP450) are the largest enzyme family in plant metabolism and widely involved in the biosynthesis of primary and secondary metabolites. Foxtail millet (Setaria italica (L.) P. Beauv) can respond to abiotic stress through a highly complex polygene regulatory network, in which the SiCYP450 family is also involved. Although the CYP450 superfamily has been systematically studied in a few species, the research on the CYP450 superfamily in foxtail millet has not been completed. In this study, three hundred and thirty-one SiCYP450 genes were identified in the foxtail millet genome by bioinformatics methods, which were divided into four groups, including forty-six subgroups. One hundred and sixteen genes were distributed in thirty-three tandem duplicated gene clusters. Chromosome mapping showed that SiCYP450 was distributed on seven chromosomes. In the SiCYP450 family of foxtail millet, 20 conserved motifs were identified. Cis-acting elements in the promoter region of SiCYP450 genes showed that hormone response elements were found in all SiCYP450 genes. Of the three hundred and thirty-one SiCYP450 genes, nine genes were colinear with the Arabidopsis thaliana genes. Two hundred SiCYP450 genes were colinear with the Setaria viridis genes, including two hundred and forty-five gene duplication events. The expression profiles of SiCYP450 genes in different organs and developmental stages showed that SiCYP450 was preferentially expressed in specific tissues, and many tissue-specific genes were identified, such as SiCYP75B6, SiCYP96A7, SiCYP71A55, SiCYP71A61, and SiCYP71A62 in the root, SiCYP78A1 and SiCYP94D9 in leaves, and SiCYP78A6 in the ear. The RT-PCR data showed that SiCYP450 could respond to abiotic stresses, ABA, and herbicides in foxtail millet. Among them, the expression levels of SiCYP709B4, SiCYP71A11, SiCYP71A14, SiCYP78A1, SiCYP94C3, and SiCYP94C4 were significantly increased under the treatment of mesotrione, florasulam, nicosulfuron, fluroxypyr, and sethoxydim, indicating that the same gene might respond to multiple herbicides. The results of this study will help reveal the biological functions of the SiCYP450 family in development regulation and stress response and provide a basis for molecular breeding of foxtail millet.

Differential RNA-Seq Analysis Predicts Genes Related to Terpene Tailoring in Caryopteris × clandonensis

Enzymatic terpene functionalization is an essential part of plant secondary metabolite diversity. Within this, multiple terpene-modifying enzymes are required to enable the chemical diversity of volatile compounds essential in plant communication and defense. This work sheds light on the differentially transcribed genes within Caryopteris × clandonensis that are capable of functionalizing cyclic terpene scaffolds, which are the product of terpene cyclase action. The available genomic reference was subjected to further improvements to provide a comprehensive basis, where the number of contigs was minimized. RNA-Seq data of six cultivars, Dark Knight, Grand Bleu, Good as Gold, Hint of Gold, Pink Perfection, and Sunny Blue, were mapped on the reference, and their distinct transcription profile investigated. Within this data resource, we detected interesting variations and additionally genes with high and low transcript abundancies in leaves of Caryopteris × clandonensis related to terpene functionalization. As previously described, different cultivars vary in their modification of monoterpenes, especially limonene, resulting in different limonene-derived molecules. This study focuses on predicting the cytochrome p450 enzymes underlying this varied transcription pattern between investigated samples. Thus, making them a reasonable explanation for terpenoid differences between these plants. Furthermore, these data provide the basis for functional assays and the verification of putative enzyme activities.

Identification and characterization of cytochrome P450 CYP77A59 of loquat (Rhaphiolepis bibas) responsible for biosynthesis of phenylacetonitrile, a floral nitrile compound

MAIN CONCLUSION

Cytochrome P450 CYP77A59 is responsible for the biosynthesis of phenylacetonitrile in loquat flowers. Flowers of some plants emit volatile nitrile compounds, but the biosynthesis of these compounds is unclear. Loquat (Rhaphiolepis bibas) flowers emit characteristic N-containing volatiles, such as phenylacetonitrile (PAN), (E/Z)-phenylacetaldoxime (PAOx), and (2-nitroethyl)benzene (NEB). These volatiles likely play a defense role against pathogens and insects. PAN and NEB are commonly biosynthesized from L-phenylalanine via (E/Z)-PAOx. Two cytochrome P450s-CYP79D80 and "promiscuous fatty acid ω-hydroxylase" CYP94A90, which catalyze the formation of (E/Z)-PAOx from L-phenylalanine and NEB from (E/Z)-PAOx, respectively-are involved in NEB biosynthesis. However, the enzymes catalyzing the formation of PAN from (E/Z)-PAOx in loquat have not been identified. In this study, we aimed to identify candidate cytochrome P450s catalyzing PAN formation in loquat flowers. Yeast whole-cell biocatalyst assays showed that among nine candidate cytochrome P450s, CYP77A58 and CYP77A59 produced PAN from (E/Z)-PAOx. CYP77As catalyzed the dehydration of aldoximes, which is atypical of cytochrome P450; the reaction was NADPH-dependent, with an optimum temperature and pH of 40 °C and 8.0, respectively. CYP77As acted on (E/Z)-PAOx, (E/Z)-4-hydroxyphenylacetaldoxime, and (E/Z)-indole-3-acetaldoxime. Previously characterized CYP77As are known to hydroxylate fatty acids; loquat CYP77As did not act on tested fatty acids. We observed higher expression of CYP77A59 in flowers than in buds; expression of CYP77A58 was remarkably reduced in the flowers. Because the flowers, but not buds, emit PAN, CYP77A59 is likely responsible for the biosynthesis of PAN in loquat flowers. This study will help us understand the biosynthesis of floral nitrile compounds.

Chromosome-scale genome assembly of Prunus pusilliflora provides novel insights into genome evolution, disease resistance, and dormancy release in Cerasus L

Prunus pusilliflora is a wild cherry germplasm resource distributed mainly in Southwest China. Despite its ornamental and economic value, a high-quality assembled P. pusilliflora genome is unavailable, hindering our understanding of its genetic background, population diversity, and evolutionary processes. Here, we de novo assembled a chromosome-scale P. pusilliflora genome using Oxford Nanopore, Illumina, and chromosome conformation capture sequencing. The assembled genome size was 309.62 Mb, with 76 scaffolds anchored to eight pseudochromosomes. We predicted 33 035 protein-coding genes, functionally annotated 98.27% of them, and identified repetitive sequences covering 49.08% of the genome. We found that P. pusilliflora is closely related to Prunus serrulata and Prunus yedoensis, having diverged from them ~41.8 million years ago. A comparative genomic analysis revealed that P. pusilliflora has 643 expanded and 1128 contracted gene families. Furthermore, we found that P. pusilliflora is more resistant to Colletotrichum viniferum, Phytophthora capsici, and Pseudomonas syringae pv. tomato (Pst) DC3000 infections than cultivated Prunus avium. P. pusilliflora also has considerably more nucleotide-binding site-type resistance gene analogs than P. avium, which explains its stronger disease resistance. The cytochrome P450 and WRKY families of 263 and 61 proteins were divided into 42 and 8 subfamilies respectively in P. pusilliflora. Furthermore, 81 MADS-box genes were identified in P. pusilliflora, accompanying expansions of the SVP and AGL15 subfamilies and loss of the TM3 subfamily. Our assembly of a high-quality P. pusilliflora genome will be valuable for further research on cherries and molecular breeding.

Lactucin Synthase Inactivation Boosts the Accumulation of Anti-inflammatory 8-Deoxylactucin and Its Derivatives in Chicory (Cichorium intybus L.)

For several sesquiterpene lactones (STLs) found in Asteraceae plants, very interesting biomedical activities have been demonstrated. Chicory roots accumulate the guaianolide STLs 8-deoxylactucin, lactucin, and lactucopicrin predominantly in oxalated forms in the latex. In this work, a supercritical fluid extract fraction of chicory STLs containing 8-deoxylactucin and 11β,13-dihydro-8-deoxylactucin was shown to have anti-inflammatory activity in an inflamed intestinal mucosa model. To increase the accumulation of these two compounds in chicory taproots, the lactucin synthase that takes 8-deoxylactucin as the substrate for the regiospecific hydroxylation to generate lactucin needs to be inactivated. Three candidate cytochrome P450 enzymes of the CYP71 clan were identified in chicory. Their targeted inactivation using the CRISPR/Cas9 approach identified CYP71DD33 to have lactucin synthase activity. The analysis of the terpene profile of the taproots of plants with edits in CYP71DD33 revealed a nearly complete elimination of the endogenous chicory STLs lactucin and lactucopicrin and their corresponding oxalates. Indeed, in the same lines, the interruption of biosynthesis resulted in a strong increase of 8-deoxylactucin and its derivatives. The enzyme activity of CYP71DD33 to convert 8-deoxylactucin to lactucin was additionally demonstrated in vitro using yeast microsome assays. The identified chicory lactucin synthase gene is predominantly expressed in the chicory latex, indicating that the late steps in the STL biosynthesis take place in the latex. This study contributes to further elucidation of the STL pathway in chicory and shows that root chicory can be positioned as a crop from which different health products can be extracted.

Down-regulation of tomato STEROL GLYCOSYLTRANSFERASE 1 perturbs plant development and facilitates viroid infection

Potato spindle tuber viroid (PSTVd) is a plant pathogen naturally infecting economically important crops such as tomato (Solanum lycopersicum). Here, we aimed to engineer tomato plants highly resistant to PSTVd and developed several S. lycopersicum lines expressing an artificial microRNA (amiRNA) against PSTVd (amiR-PSTVd). Infectivity assays revealed that amiR-PSTVd-expressing lines were not resistant but instead hypersusceptible to the viroid. A combination of phenotypic, molecular, and metabolic analyses of amiRNA-expressing lines non-inoculated with the viroid revealed that amiR-PSTVd was accidentally silencing the tomato STEROL GLYCOSYLTRANSFERASE 1 (SlSGT1) gene, which caused late developmental and reproductive defects such as leaf epinasty, dwarfism, or reduced fruit size. Importantly, two independent transgenic tomato lines each expressing a different amiRNA specifically designed to target SlSGT1 were also hypersusceptible to PSTVd, thus demonstrating that down-regulation of SlSGT1 was responsible for the viroid-hypersusceptibility phenotype. Our results highlight the role of sterol glycosyltransferases in proper plant development and indicate that the imbalance of sterol glycosylation levels favors viroid infection, most likely by facilitating viroid movement.

Genes encoding cytochrome P450 monooxygenases and glutathione S-transferases associated with herbicide resistance evolved before the origin of land plants

Cytochrome P450 (CYP) monooxygenases and glutathione S-transferases (GST) are enzymes that catalyse chemical modifications of a range of organic compounds. Herbicide resistance has been associated with higher levels of CYP and GST gene expression in some herbicide-resistant weed populations compared to sensitive populations of the same species. By comparing the protein sequences of 9 representative species of the Archaeplastida-the lineage which includes red algae, glaucophyte algae, chlorophyte algae, and streptophytes-and generating phylogenetic trees, we identified the CYP and GST proteins that existed in the common ancestor of the Archaeplastida. All CYP clans and all but one land plant GST classes present in land plants evolved before the divergence of streptophyte algae and land plants from their last common ancestor. We also demonstrate that there are more genes encoding CYP and GST proteins in land plants than in algae. The larger numbers of genes among land plants largely results from gene duplications in CYP clans 71, 72, and 85 and in the GST phi and tau classes [1,2]. Enzymes that either metabolise herbicides or confer herbicide resistance belong to CYP clans 71 and 72 and the GST phi and tau classes. Most CYP proteins that have been shown to confer herbicide resistance are members of the CYP81 family from clan 71. These results demonstrate that the clan and class diversity in extant plant CYP and GST proteins had evolved before the divergence of land plants and streptophyte algae from a last common ancestor estimated to be between 515 and 474 million years ago. Then, early in embryophyte evolution during the Palaeozoic, gene duplication in four of the twelve CYP clans, and in two of the fourteen GST classes, led to the large numbers of CYP and GST proteins found in extant land plants. It is among the genes of CYP clans 71 and 72 and GST classes phi and tau that alleles conferring herbicide resistance evolved in the last fifty years.

Structure-function analysis of CYP719As involved in methylenedioxy bridge-formation in the biosynthesis of benzylisoquinoline alkaloids and its de novo production

Benzylisoquinoline alkaloids (BIAs) are a type of secondary metabolite with clinical application value. (S)-stylopine is a special BIA which contains methylenedioxy bridge structures. CYP719As could catalyze the methylenedioxy bridge-formation on the A or D rings of protoberberine alkaloids, while displaying significant substrate regiospecificity. To explore the substrate preference of CYP719As, we cloned and identified five CyCYP719A candidates from Corydalis yanhusuo. Two CyCYP719As (CyCYP719A39 and CyCYP719A42) with high catalytic efficiency for the methylenedioxy bridge-formation on the D or A rings were characterized, respectively. The residues (Leu 294 for CyCYP719A42 and Asp 289 for CyCYP719A39) were identified as the key to controlling the regioselectivity of CYP719As affecting the methylenedioxy bridge-formation on the A or D rings by homology modeling and mutation analysis. Furthermore, for de novo production of BIAs, CyCYP719A39, CyCYP719A42, and their mutants were introduced into the (S)-scoulerine-producing yeast to produce 32 mg/L (S)-stylopine. These results lay a foundation for understanding the structure-function relationship of CYP719A-mediated methylenedioxy bridge-formation and provide yeast strains for the BIAs production by synthetic biology.

Genome-wide analysis of cytochrome P450 genes in Citrus clementina and characterization of a CYP gene encoding flavonoid 3'-hydroxylase

Cytochrome P450s (CYPs) are the largest family of enzymes in plant and play multifarious roles in development and defense but the available information about the CYP superfamily in citrus is very limited. Here we provide a comprehensive genome-wide analysis of the CYP superfamily in Citrus clementina genome, identifying 301 CYP genes grouped into ten clans and 49 families. The characteristics of both gene structures and motif compositions strongly supported the reliability of the phylogenetic relationship. Duplication analysis indicated that tandem duplication was the major driving force of expansion for this superfamily. Promoter analysis revealed numerous cis-acting elements related to various responsiveness. RNA-seq data elucidated their expression patterns in citrus fruit peel both during development and in response to UV-B. Furthermore, we characterize a UV-B-induced CYP gene (Ciclev10019637m, designated CitF3'H) as a flavonoid 3'-hydroxylase for the first time. CitF3'H catalyzed numerous flavonoids and favored naringenin in yeast assays. Virus-induced silencing of CitF3'H in citrus seedlings significantly reduced the levels of 3'-hydroxylated flavonoids and their derivatives. These results together with the endoplasmic reticulum-localization of CitF3'H in plant suggest that this enzyme is responsible for the biosynthesis of 3'-hydroxylated flavonoids in citrus. Taken together, our findings provide extensive information about the CYP superfamily in citrus and contribute to further functional verification.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

Downregulation of barley ferulate 5-hydroxylase dramatically alters straw lignin structure without impact on mechanical properties

Barley is a major cereal crop for temperate climates, and a diploid genetic model for polyploid wheat. Cereal straw biomass is an attractive source of feedstock for green technologies but lignin, a key determinant of feedstock recalcitrance, complicates bio-conversion processes. However, manipulating lignin content to improve the conversion process could negatively affect agronomic traits. An alternative approach is to manipulate lignin composition which influences the physical and chemical properties of straw. This study validates the function of a barley ferulate 5-hydroxylase gene and demonstrates that its downregulation using the RNA-interference approach substantially impacts lignin composition. We identified five barley genes having putative ferulate 5-hydroxylase activity. Downregulation of HvF5H1 substantially reduced the lignin syringyl/guaiacyl (S/G) ratio in straw while the lignin content, straw mechanical properties, plant growth habit, and grain characteristics all remained unaffected. Metabolic profiling revealed significant changes in the abundance of 173 features in the HvF5H1-RNAi lines. The drastic changes in the lignin polymer of transgenic lines highlight the plasticity of barley lignification processes and the associated potential for manipulating and improving lignocellulosic biomass as a feedstock for green technologies. On the other hand, our results highlight some differences between the lignin biosynthetic pathway in barley, a temperate climate grass, and the warm climate grass, rice, and underscore potential diversity in the lignin biosynthetic pathways in grasses.

Tissue-preferential recruitment of electron transfer chains for cytochrome P450-catalyzed phenolic biosynthesis

Cytochrome P450 system consists of P450 monooxygenase and redox pattern(s). While the importance of monooxygenases in plant metabolism is well documented, the metabolic roles of the related redox components have been largely overlooked. Here, we show that distinct electron transfer chains are recruited in phenylpropanoid-monolignol P450 systems to support the synthesis and distribution of different classes of phenolics in different plant tissues. While Arabidopsis cinnamate 4-hydroxylase adopts conventional NADPH-cytochrome P450 oxidoreductase (CPR) electron transfer chain for its para-hydroxylation reaction, ferulate 5-hydroxylase uses both NADPH-CPR-cytochrome b₅ (CB5) and NADH-cytochrome b₅ reductase-CB5 chains to support benzene ring 5-hydroxylation, in which the former route is primarily recruited in the stem for syringyl lignin synthesis, while the latter dominates in the syntheses of 5-hydroxylated phenolics in seeds and seed coat suberin. Our study unveils an additional layer of complexity and versatility of P450 system that the plants evolved for diversifying phenolic repertoires.

Genetic architecture of inducible and constitutive metabolic profile related to drought resistance in qingke (Tibetan hulless barley)

Qingke (Tibetan hulless barley, Hordeum vulgare L. var. nudum) is the primary food crop on the Tibet Plateau, the long-term drought and other harsh environments makes qingke an important resource for the study of abiotic resistance. Here, we evaluated the drought sensitivity of 246 qingke varieties. Genome-wide association studies (GWAS) found that root-specific expressed gene CYP84 may be involved in the regulation of drought resistance. Based on widely targeted metabolic profiling, we identified 2,769 metabolites in qingke leaves, of which 302 were significantly changed in response to drought stress, including 4-aminobutyric acid (GABA), proline, sucrose and raffinose. Unexpectedly, these drought-induced metabolites changed more violently in drought-sensitive qingkes, while the constitutive metabolites that had little response to drought stress, such as C-glycosylflavonoids and some amino acids, accumulated excessively in drought-resistant qingkes. Combined with metabolite-based genome-wide association study (mGWAS), a total of 1,006 metabolites under optimal condition and 1,031 metabolites under mild drought stress had significant associated loci. As a marker metabolite induced by drought stress, raffinose was significantly associated with two conservatively adjacent α-galactosidase genes, qRT-PCR suggests that these two genes may jointly regulate the raffinose content in qingke. Besides, as constituent metabolites with stable differences between drought-sensitive and drought-resistant qingkes, a class of C-glycosylflavonoids are simultaneously regulated by a UDP-glucosyltransferase gene. Overall, we performed GWAS for sensitivity and widely targeted metabolites during drought stress in qingke for the first time, which provides new insights into the response mechanism of plant drought stress and drought resistance breeding.

Phylogenomic analysis of cytochrome P450 multigene family and its differential expression analysis in pepper (Capsicum annuum L.)

Plant cytochrome P450 is a multifamily enzyme widely involved in biochemical reactions for the synthesis of antioxidants, pigments, structural polymers, and defense-related compounds. Pepper (Capsicum annuum L.) is an economically important plant. A comprehensive identification and characterization of P450 genes would provide valuable information on the evolutionary relationships of genes and their functional characteristics. In this study, we identified P450 genes in pepper with the aid of bioinformatics methods to investigate the phylogenetic relation, gene structure, chromosomal localization, duplicated events, and collinearity among Solanaceae species. We identified and classified 478 genes of P450 from the pepper genome into two major clades and nine subfamilies through phylogenetic analysis. Massive duplication events were found in the P450 gene family, which may explain the expansion of the P450 gene family. In addition, we also found that these duplication genes may have undergone strict purification selection during evolution. Gene expression analysis showed that some P450 genes that belong to clan 71 in pepper may play an important role in placenta and pericarp development. Through quantitative real-time polymerase chain reaction and transcriptome analysis, we also found that many P450 genes were related to defensive and phytohormone response in pepper. These findings provide insight for further studies to identify the biological functions of the P450 genes in pepper.

Differential Water Deficit in Leaves Is a Principal Factor Modifying Barley Response to Drought Stress

In response to environmental stress, plants activate complex signalling, including being dependent on reactive oxygen-nitrogen-sulphur species. One of the key abiotic stresses is drought. As a result of drought, changes in the level of hydration of the plant occur, which obviously entails various metabolic alternations. The primary aim of this study was to determine the relationship between the response of barley to drought and the intensity of stress, therefore investigations were performed under various levels of water saturation deficit (WSD) in leaves at 15%, 30%, and 50%. In barley subjected to drought, most significant changes occurred under a slight dehydration level at 15%. It was observed that the gene expression of 9-cis-epoxycarotenoid dioxygenases, enzymes involved in ABA biosynthesis, increased significantly, and led to a higher concentration of ABA. This was most likely the result of an increase in the gene expression and enzyme activity of L-cysteine desulfhydrase, which is responsible for H₂S synthesis. Our results suggest that the differential water deficit in leaves underlies the activation of an appropriate defence, with ABA metabolism at the centre of these processes. Furthermore, at 15% WSD, a dominant contribution of H₂O₂-dependent signalling was noted, but at 30% and 50% WSD, significant NO-dependent signalling occurred.

Protocorm-like-body extract of Phalaenopsis aphrodite combats watermelon fruit blotch disease

Bacterial fruit blotch, caused by the seedborne gram-negative bacterium Acidovorax citrulli, is one of the most destructive bacterial diseases of cucurbits (gourds) worldwide. Despite its prevalence, effective and reliable means to control bacterial fruit blotch remain limited. Transcriptomic analyses of tissue culture-based regeneration processes have revealed that organogenesis-associated cellular reprogramming is often associated with upregulation of stress- and defense-responsive genes. Yet, there is limited evidence supporting the notion that the reprogrammed cellular metabolism of the regenerated tissued confers bona fide antimicrobial activity. Here, we explored the anti-bacterial activity of protocorm-like-bodies (PLBs) of Phalaenopsis aphrodite. Encouragingly, we found that the PLB extract was potent in slowing growth of A. citrulli, reducing the number of bacteria attached to watermelon seeds, and alleviating disease symptoms of watermelon seedlings caused by A. citrulli. Because the anti-bacterial activity can be fractionated chemically, we predict that reprogrammed cellular activity during the PLB regeneration process produces metabolites with antibacterial activity. In conclusion, our data demonstrated the antibacterial activity in developing PLBs and revealed the potential of using orchid PLBs to discover chemicals to control bacterial fruit blotch disease.

Molecular Characterization of an Isoflavone 2'-Hydroxylase Gene Revealed Positive Insights into Flavonoid Accumulation and Abiotic Stress Tolerance in Safflower

Flavonoids with significant therapeutic properties play an essential role in plant growth, development, and adaptation to various environments. The biosynthetic pathway of flavonoids has long been studied in plants; however, its regulatory mechanism in safflower largely remains unclear. Here, we carried out comprehensive genome-wide identification and functional characterization of a putative cytochrome P45081E8 gene encoding an isoflavone 2'-hydroxylase from safflower. A total of 15 CtCYP81E genes were identified from the safflower genome. Phylogenetic classification and conserved topology of CtCYP81E gene structures, protein motifs, and cis-elements elucidated crucial insights into plant growth, development, and stress responses. The diverse expression pattern of CtCYP81E genes in four different flowering stages suggested important clues into the regulation of secondary metabolites. Similarly, the variable expression of CtCYP81E8 during multiple flowering stages further highlighted a strong relationship with metabolite accumulation. Furthermore, the orchestrated link between transcriptional regulation of CtCYP81E8 and flavonoid accumulation was further validated in the yellow- and red-type safflower. The spatiotemporal expression of CtCYP81E8 under methyl jasmonate, polyethylene glycol, light, and dark conditions further highlighted its likely significance in abiotic stress adaption. Moreover, the over-expressed transgenic Arabidopsis lines showed enhanced transcript abundance in OE-13 line with approximately eight-fold increased expression. The upregulation of AtCHS, AtF3'H, and AtDFR genes and the detection of several types of flavonoids in the OE-13 transgenic line also provides crucial insights into the potential role of CtCYP81E8 during flavonoid accumulation. Together, our findings shed light on the fundamental role of CtCYP81E8 encoding a putative isoflavone 2'-hydroxylase via constitutive expression during flavonoid biosynthesis.

Functional genomics-enabled characterization of CYP81B140 and CYP81B141 from Plumbago zeylanica L. substantiates their involvement in plumbagin biosynthesis

Novel cytochrome P450s, CYP81B140 and CYP81B141 from Plumbago zeylanica were functionally characterized to understand their involvement in polyketide plumbagin biosynthesis. Further, we propose 3-methyl-1-8-naphthalenediol and isoshinanolone as intermediates for plumbagin biosynthesis. Plumbago zeylanica L. (P. zeylanica) is a medicinally important plant belonging to the family Plumbaginaceae. It comprises the most abundant naphthoquinone plumbagin having anti-cancer activity. Only the polyketide synthase (PKS) enzyme has been identified from the biosynthetic pathway which catalyzes iterative condensation of acetyl-CoA and malonyl-CoA molecules. The plumbagin biosynthesis involves hydroxylation, oxidation, hydration and dehydration of intermediate compounds which are expected to be catalyzed by cytochrome P450s (CYPs). To identify the CYPs, co-expression analysis was carried out using PKS as a candidate gene. Out of the eight identified CYPs, CYP81B140 and CYP81B141 have similar expression with PKS and belong to the CYP81 family. Phylogenetic analysis suggested that CYP81B140 and CYP81B141 cluster with CYPs from CYP81B, CYP81D, CYP81E and CYP81AA subfamilies which are known to be involved in the hydroxylation and oxidation reactions. Moreover, artificial microRNA-mediated transient individual silencing and co-silencing of CYP81B140 and CYP81B141 significantly reduced plumbagin and increased the 3-methyl-1-8-naphthalenediol and isoshinanolone content. Based on metabolite analysis, we proposed that 3-methyl-1-8-naphthalenediol and isoshinanolone function as intermediates for plumbagin biosynthesis. Transient silencing, over-expression and docking analysis revealed that CYP81B140 is involved in C-1 oxidation, C-4 hydroxylation and [C2-C3] hydration of 3-methyl-1-8-naphthalenediol to form isoshinanolone, whereas CYP81B141 is catalyzing [C2-C3] dehydration and C-4 oxidation of isoshinanolone to form plumbagin. Our results indicated that both CYP81B140 and CYP81B141 are promiscuous and necessary for plumbagin biosynthesis. This is the first report of identification and functional characterization of P. zeylanica-specific CYPs involved in plumbagin biosynthetic pathway and in general hexaketide synthesis in plants.

Cytochrome b 5: A versatile electron carrier and regulator for plant metabolism

Cytochrome b ₅ (CB5) is a small heme-binding protein, known as an electron donor delivering reducing power to the terminal enzymes involved in oxidative reactions. In plants, the CB5 protein family is substantially expanded both in its isoform numbers and cellular functions, compared to its yeast and mammalian counterparts. As an electron carrier, plant CB5 proteins function not only in fatty acid desaturation, hydroxylation and elongation, but also in the formation of specialized metabolites such as flavonoids, phenolic esters, and heteropolymer lignin. Furthermore, plant CB5s are found to interact with different non-catalytic proteins such as ethylene signaling regulator, cell death inhibitor, and sugar transporters, implicating their versatile regulatory roles in coordinating different metabolic and cellular processes, presumably in respect to the cellular redox status and/or carbon availability. Compared to the plentiful studies on biochemistry and cellular functions of mammalian CB5 proteins, the cellular and metabolic roles of plant CB5 proteins have received far less attention. This article summarizes the fragmentary information pertaining to the discovery of plant CB5 proteins, and discusses the conventional and peculiar functions that plant CB5s might play in different metabolic and cellular processes. Gaining comprehensive insight into the biological functions of CB5 proteins could offer effective biotechnological solutions to tailor plant chemodiversity and cellular responses to environment stimuli.

Cytochrome P450 monooxygenase-mediated tailoring of triterpenoids and steroids in plants

The cytochrome P450 monooxygenase (CYP) superfamily comprises hemethiolate enzymes that perform remarkable regio- and stereospecific oxidative chemistry. As such, CYPs are key agents for the structural and functional tailoring of triterpenoids, one of the largest classes of plant natural products with widespread applications in pharmaceuticals, food, cosmetics, and agricultural industries. In this review, we provide a full overview of 149 functionally characterised CYPs involved in the biosynthesis of triterpenoids and steroids in primary as well as in specialised metabolism. We describe the phylogenetic distribution of triterpenoid- and steroid-modifying CYPs across the plant CYPome, present a structure-based summary of their reactions, and highlight recent examples of particular interest to the field. Our review therefore provides a comprehensive up-to-date picture of CYPs involved in the biosynthesis of triterpenoids and steroids in plants as a starting point for future research.

Comparison of buckwheat genomes reveals the genetic basis of metabolomic divergence and ecotype differentiation

Golden buckwheat (Fagopyrum dibotrys or Fagopyrum cymosum) and Tartary buckwheat (Fagopyrum tataricum) belong to the Polygonaceae and the Fagopyrum genus is rich in flavonoids. Golden buckwheat is a wild relative of Tartary buckwheat, yet golden buckwheat is a traditional Chinese herbal medicine and Tartary buckwheat is a food crop. The genetic basis of adaptive divergence between these two buckwheats is poorly understood. Here, we assembled a high-quality chromosome-level genome of golden buckwheat and found a one-to-one syntenic relationship with the chromosomes of Tartary buckwheat. Two large inversions were identified that differentiate golden buckwheat and Tartary buckwheat. Metabolomic and genetic comparisons of golden buckwheat and Tartary buckwheat indicate an amplified copy number of FdCHI, FdF3H, FdDFR, and FdLAR gene families in golden buckwheat, and a parallel increase in medicinal flavonoid content. Resequencing of 34 wild golden buckwheat accessions across the two morphologically distinct ecotypes identified candidate genes, including FdMYB44 and FdCRF4, putatively involved in flavonoid accumulation and differentiation of plant architecture, respectively. Our comparative genomic study provides abundant genomic resources of genomic divergent variation to improve buckwheat with excellent nutritional and medicinal value.

Identification of key amino acid residues toward improving the catalytic activity and substrate specificity of plant-derived cytochrome P450 monooxygenases CYP716A subfamily enzyme for triterpenoid production in Saccharomyces cerevisiae

Triterpenoids constitute a group of specialized plant metabolites with wide structural diversity and high therapeutic value for human health. Cytochrome P450 monooxygenases (CYP) are a family of enzymes important for generating the structural diversity of triterpenoids by catalyzing the site-specific oxidization of the triterpene backbone. The CYP716 enzyme family has been isolated from various plant families as triterpenoid oxidases; however, their experimental crystal structures are not yet available and the detailed catalytic mechanism remains elusive. Here, we address this challenge by integrating bioinformatics approaches with data from other CYP families. Medicago truncatula CYP716A12, the first functionally characterized CYP716A subfamily enzyme, was chosen as the model for this study. We performed homology modeling, structural alignment, in silico site-directed mutagenesis, and molecular docking analysis to search and screen key amino acid residues relevant to the catalytic activity and substrate specificity of the CYP716A subfamily enzyme in triterpenoid biosynthesis. An in vivo functional analysis using engineered yeast that endogenously produced plant-derived triterpenes was performed to elucidate the results. When the amino acids in the signature region and substrate recognition sites (SRSs) were substituted, the product profile of CYP716A12 was modified. We identified amino acid residues that control the substrate contraction of the enzyme (D292) and engineered the enzyme to improve its catalytic activity and substrate specificity (D122, I212, and Q358) for triterpenoid biosynthesis. In addition, we demonstrated the versatility of this strategy by changing the properties of key residues in SRSs to improve the catalytic activity of Arabidopsis thaliana CYP716A1 (S356) and CYP716A2 (M206, F210) at C-28 on the triterpene backbone. This research has the potential to help in the production of desired triterpenoids in engineered yeast by increasing the catalytic activity and substrate specificity of plant CYP716A subfamily enzymes.

Genome-wide Association Study for Yield and Yield-Related Traits in Diverse Blackgram Panel (Vigna mungo L. Hepper) Reveals Novel Putative Alleles for Future Breeding Programs

Blackgram (Vigna mungo L. Hepper) is an important tropical and sub-tropical short-duration legume that is rich in dietary protein and micronutrients. Producing high-yielding blackgram varieties is hampered by insufficient genetic variability, absence of suitable ideotypes, low harvest index and susceptibility to biotic-abiotic stresses. Seed yield, a complex trait resulting from the expression and interaction of multiple genes, necessitates the evaluation of diverse germplasm for the identification of novel yield contributing traits. Henceforth, a panel of 100 blackgram genotypes was evaluated at two locations (Ludhiana and Gurdaspur) across two seasons (Spring 2019 and Spring 2020) for 14 different yield related traits. A wide range of variability, high broad-sense heritability and a high correlation of grain yield were observed for 12 out of 14 traits studied among all environments. Investigation of population structure in the panel using a set of 4,623 filtered SNPs led to identification of four sub-populations based on ad-hoc delta K and Cross entropy value. Using Farm CPU model and Mixed Linear Model algorithms, a total of 49 significant SNP associations representing 42 QTLs were identified. Allelic effects were found to be statistically significant at 37 out of 42 QTLs and 50 known candidate genes were identified in 24 of QTLs.