Diversification of P450 genes during land plant evolution

Identification and analysis of flavonoid pathway genes in responsive to drought and salinity stress in Medicago truncatula

Flavonoid compounds are widely present in various organs and tissues of different plants, playing important roles when plants are exposed to abiotic stresses. Different types of flavonoids are biosynthesized by a series of enzymes that are encoded by a range of gene families. In this study, a total of 63 flavonoid pathway genes were identified from the genome of Medicago truncatula. Gene structure analysis revealed that they all have different gene structure, with most CHS genes containing only one intron. Additionally, analysis of promoter sequences revealed that many cis-acting elements responsive to abiotic stress are located in the promoter region of flavonoid pathway genes. Furthermore, analysis on M. truncatula gene chip data revealed significant changes in expression level of most flavonoid pathway genes under the induction of salt or drought treatment. qRT-PCR further confirmed significant increase in expression level of several flavonoid pathway genes under NaCl and mannitol treatments, with CHS1, CHS9, CHS10, F3'H4 and F3'H5 genes showing significant up-regulation, indicating they are key genes in response to abiotic stress in M. truncatula. In summary, our study identified key flavonoid pathway genes that were involved in salt and drought response, which provides important insights into possible modification of flavonoid pathway genes for molecular breeding of forage grass with improved abiotic resistance.

An ancient role for CYP73 monooxygenases in phenylpropanoid biosynthesis and embryophyte development

The phenylpropanoid pathway is one of the plant metabolic pathways most prominently linked to the transition to terrestrial life, but its evolution and early functions remain elusive. Here, we show that activity of the t-cinnamic acid 4-hydroxylase (C4H), the first plant-specific step in the pathway, emerged concomitantly with the CYP73 gene family in a common ancestor of embryophytes. Through structural studies, we identify conserved CYP73 residues, including a crucial arginine, that have supported C4H activity since the early stages of its evolution. We further demonstrate that impairing C4H function via CYP73 gene inactivation or inhibitor treatment in three bryophyte species-the moss Physcomitrium patens, the liverwort Marchantia polymorpha and the hornwort Anthoceros agrestis-consistently resulted in a shortage of phenylpropanoids and abnormal plant development. The latter could be rescued in the moss by exogenous supply of p-coumaric acid, the product of C4H. Our findings establish the emergence of the CYP73 gene family as a foundational event in the development of the plant phenylpropanoid pathway, and underscore the deep-rooted function of the C4H enzyme in embryophyte biology.

Cytochrome P450 CYP81A104 in Eleusine indica confers resistance to multiherbicide with different modes of action

BACKGROUND

Developing herbicide-resistant (HR) crop cultivars is an efficient way to control weeds and minimize crop yield losses. However, widespread and long-term herbicide application has led to the evolution of resistant weeds. Here, we established a resistant (R) E. indica population, collected from imidazolinone-resistant rice cultivar fields.

RESULTS

The R population evolved 4.5-fold resistance to imazamox. Acetolactate synthase (ALS) gene sequencing and ALS activity assays excluded the effect of target-site resistance in this population. P450 inhibitor malathion pretreatment significantly reversed resistance to imazamox. RNA sequencing showed that a P450 gene CYP81A104 was expressed higher in R versus susceptible (S) plants. Arabidopsis overexpressing CYP81A104 showed resistance to ALS inhibitors (imazamox, tribenuron-methyl, penoxsulam and flucarbazone-sodium), PSII inhibitor (bentazone), hydroxyphenyl pyruvate dioxygenase inhibitor (mesotrione) and auxin mimics (MCPA), which was generally consistent with the results presented in the R population.

CONCLUSION

This study confirmed that the CYP81A104 gene endowed resistance to multiherbicides with different modes-of-action. Our findings provide an insight into the molecular characteristics of resistance and contribute to formulating an appropriate strategy for weed management in HR crops. © 2024 Society of Chemical Industry.

Characterization of the ddt1 Mutant in Rice and Its Impact on Plant Height Reduction and Water Use Efficiency

Rice (Oryza sativa L.), a fundamental global staple, nourishes over half of the world's population. The identification of the ddt1 mutant in rice through EMS mutagenesis of the indica cultivar Shuhui527 revealed a dwarf phenotype, characterized by reduced plant height, smaller grain size, and decreased grain weight. Detailed phenotypic analysis and map-based cloning pinpointed the mutation to a single-base transversion in the LOC_Os03g04680 gene, encoding a cytochrome P450 enzyme, which results in a premature termination of the protein. Functional complementation tests confirmed LOC_Os03g04680 as the DDT1 gene responsible for the observed phenotype. We further demonstrated that the ddt1 mutation leads to significant alterations in gibberellic acid (GA) metabolism and signal transduction, evidenced by the differential expression of key GA-related genes such as OsGA20OX2, OsGA20OX3, and SLR1. The mutant also displayed enhanced drought tolerance, as indicated by higher survival rates, reduced water loss, and rapid stomatal closure under drought conditions. This increased drought resistance was linked to the mutant's improved antioxidant capacity, with elevated activities of antioxidant enzymes and higher expression levels of related genes. Our findings suggest that DDT1 plays a crucial role in regulating both plant height and drought stress responses. The potential for using gene editing of DDT1 to mitigate the dwarf phenotype while retaining improved drought resistance offers promising avenues for rice improvement.

SlCYP94B18 and SlCYP94B19 monooxygenases for the catabolic turnover of jasmonates in tomato leaves

(3R,7S)-Jasmonoyl-L-isoleucine (JA-Ile) is a plant hormone that regulates plant defense responses and other physiological functions. The mechanism of attenuation of JA-Ile signaling in the plant body is essential because prolonged JA-Ile signaling can be detrimental to plant survival. In Arabidopsis thaliana, the cytochrome P450 monooxygenases, CYP94B1/B3/C1, inactivate JA-Ile by converting it into 12-hydroxy-jasmonoyl-L-isoleucine (12-OH-JA-Ile), and CYP94C1 converts 12-OH-JA-Ile into 12-carboxy-jasmonoyl-L-isoleucine (12-COOH-JA-Ile). In the present study, we aimed to identify the cytochrome P450 monooxygenases involved in the catabolic pathway of JA-Ile in tomato leaves. Based on a gene expression screening of SlCYP94 subfamily monooxygenases using qPCR and the time-course of JA-Ile catabolism, we identified SlCYP94B18 and SlCYP94B19 expressed in tomato leaves as candidate monooxygenases catalyzing the two-step catabolism of JA-Ile. An in vitro enzymatic assay using a yeast expression system revealed that these enzymes efficiently converted JA-Ile to 12-OH-JA-Ile, and then to 12-COOH-JA-Ile. SlCYP94B18 and SlCYP94B19 also catalyzed the oxidative catabolism of several JA-amino acid conjugates (JA-AAs), JA-Leu and JA-Val, in tomatoes. These results suggest that SlCYP94B18 and SlCYP94B19 plays a role in the two-step oxidation of JA-AAs, suggesting their broad involvement in regulating jasmonate signaling in tomatoes. Our results contribute to a deeper understanding of jasmonate signaling in tomatoes and may help to improve tomato cultivation and quality.

Multisubstrate specificity shaped the complex evolution of the aminotransferase family across the tree of life

Aminotransferases (ATs) are an ancient enzyme family that play central roles in core nitrogen metabolism, essential to all organisms. However, many of the AT enzyme functions remain poorly defined, limiting our fundamental understanding of the nitrogen metabolic networks that exist in different organisms. Here, we traced the deep evolutionary history of the AT family by analyzing AT enzymes from 90 species spanning the tree of life (ToL). We found that each organism has maintained a relatively small and constant number of ATs. Mapping the distribution of ATs across the ToL uncovered that many essential AT reactions are carried out by taxon-specific AT enzymes due to wide-spread nonorthologous gene displacements. This complex evolutionary history explains the difficulty of homology-based AT functional prediction. Biochemical characterization of diverse aromatic ATs further revealed their broad substrate specificity, unlike other core metabolic enzymes that evolved to catalyze specific reactions today. Interestingly, however, we found that these AT enzymes that diverged over billion years share common signatures of multisubstrate specificity by employing different nonconserved active site residues. These findings illustrate that AT family enzymes had leveraged their inherent substrate promiscuity to maintain a small yet distinct set of multifunctional AT enzymes in different taxa. This evolutionary history of versatile ATs likely contributed to the establishment of robust and diverse nitrogen metabolic networks that exist throughout the ToL. The study provides a critical foundation to systematically determine diverse AT functions and underlying nitrogen metabolic networks across the ToL.

Functional Validation of the Cytochrome P450 Family PgCYP309 Gene in Panax ginseng

Ginseng (Panax ginseng C. A. Meyer) is an ancient and valuable Chinese herbal medicine, and ginsenoside, as the main active ingredient of ginseng, has received wide attention because of its various pharmacological active effects. Cytochrome P450 is the largest family of enzymes in plant metabolism and is involved in the biosynthesis of terpenoids, alkaloids, lipids, and other primary and secondary plant metabolites. It is significant to explore more PgCYP450 genes with unknown functions and reveal their roles in ginsenoside synthesis. In this study, based on the five PgCYP450 genes screened in the pre-laboratory, through the correlation analysis with the content of ginsenosides and the analysis of the interactions network of the key enzyme genes for ginsenoside synthesis, we screened out those highly correlated with ginsenosides, PgCYP309, as the target gene from among the five PgCYP450 genes. Methyl jasmonate-induced treatment of ginseng adventitious roots showed that the PgCYP309 gene responded to methyl jasmonate induction and was involved in the synthesis of ginsenosides. The PgCYP309 gene was cloned and the overexpression vector pBI121-PgCYP309 and the interference vector pART27-PgCYP309 were constructed. Transformation of ginseng adventitious roots by the Agrobacterium fermentum-mediated method and successful induction of transgenic ginseng hairy roots were achieved. The transformation rate of ginseng hairy roots with overexpression of the PgCYP309 gene was 22.7%, and the transformation rate of ginseng hairy roots with interference of the PgCYP309 gene was 40%. Analysis of ginseng saponin content and relative gene expression levels in positive ginseng hairy root asexual lines revealed a significant increase in PPD, PPT, and PPT-type monomeric saponins Re and Rg2. The relative expression levels of PgCYP309 and PgCYP716A53v2 genes were also significantly increased. PgCYP309 gene promotes the synthesis of ginsenosides, and it was preliminarily verified that PgCYP309 gene can promote the synthesis of dammarane-type ginsenosides.

Cytochrome P450 CitCYP97B modulates carotenoid accumulation diversity by hydroxylating β-cryptoxanthin in Citrus

Carotenoids in plant foods provide health benefits by functioning as provitamin A. One of the vital provitamin A carotenoids, β-cryptoxanthin, is typically plentiful in citrus fruit. However, little is known about the genetic basis of β-cryptoxanthin accumulation in citrus. Here, we performed a widely targeted metabolomic analysis of 65 major carotenoids and carotenoid derivatives to characterize carotenoid accumulation in Citrus and determine the taxonomic profile of β-cryptoxanthin. We used data from 81 newly sequenced representative accessions and 69 previously sequenced Citrus cultivars to reveal the genetic basis of β-cryptoxanthin accumulation through a genome-wide association study. We identified a causal gene, CitCYP97B, which encodes a cytochrome P450 protein whose substrate and metabolic pathways in land plants were undetermined. We subsequently demonstrated that CitCYP97B functions as a novel monooxygenase that specifically hydroxylates the β-ring of β-cryptoxanthin in a heterologous expression system. In planta experiments provided further evidence that CitCYP97B negatively regulates β-cryptoxanthin content. Using the sequenced Citrus accessions, we found that two critical structural cis-element variations contribute to increased expression of CitCYP97B, thereby altering β-cryptoxanthin accumulation in fruit. Hybridization/introgression appear to have contributed to the prevalence of two cis-element variations in different Citrus types during citrus evolution. Overall, these findings extend our understanding of the regulation and diversity of carotenoid metabolism in fruit crops and provide a genetic target for production of β-cryptoxanthin-biofortified products.

Integrated Metabolomics and Transcriptomics Analysis Reveals New Insights into Triterpene Biosynthesis in Rosa rugosa

Rosa rugosa is highly regarded for its aesthetic and therapeutic qualities. In particular, R. rugosa's flowers are known to produce essential oils containing a mixture of volatile terpenes, phenylpropanoids, and other compounds. Despite this, extensive research exists on volatile terpenes in flowers, while the knowledge of non-volatile terpenes in distinct tissues is still limited. Using UPLC-ESI-MS/MS, a comprehensive analysis of the terpene metabolites in five different tissues of R. rugosa was conducted. These metabolites accumulated in distinct tissues, and the majority of them were triterpenoids. Transcriptome data were collected from five tissues using RNA-seq. Transcriptomics and metabolomics were utilized to evaluate the triterpene biosynthesis pathway, resulting in new insights into its regulation and biosynthesis. The RrOSC10 was identified as a key enzyme in converting 2,3-oxidosqualene into α-amyrin, potentially contributing to the triterpene biosynthesis pathway. Furthermore, the expression of the RrOSC10 gene was upregulated by salinity for 0.5 h and 1 h, with subsequent downregulation at 2 h. This study lays a foundation for future research on the biosynthesis and accumulation of triterpenes in R. rugosa.

Radical ligand transfer: mechanism and reactivity governed by three-component thermodynamics

Here, we demonstrate that the relationship between reactivity and thermodynamics in radical ligand transfer chemistry can be understood if this chemistry is dissected as concerted ion-electron transfer (cIET). Namely, we investigate radical ligand transfer reactions from the perspective of thermodynamic contributions to the reaction barrier: the diagonal effect of the free energy of the reaction, and the off-diagonal effect resulting from asynchronicity and frustration, which we originally derived from the thermodynamic cycle for concerted proton-electron transfer (cPET). This study on the OH transfer reaction shows that the three-component thermodynamic model goes beyond cPET chemistry, successfully capturing the changes in radical ligand transfer reactivity in a series of model FeIII-OH⋯(diflouro)cyclohexadienyl systems. We also reveal the decisive role of the off-diagonal thermodynamics in determining the reaction mechanism. Two possible OH transfer mechanisms, in which electron transfer is coupled with either OH- and OH+ transfer, are associated with two competing thermodynamic cycles. Consequently, the operative mechanism is dictated by the cycle yielding a more favorable off-diagonal effect on the barrier. In line with this thermodynamic link to the mechanism, the transferred OH group in OH-/electron transfer retains its anionic character and slightly changes its volume in going from the reactant to the transition state. In contrast, OH+/electron transfer develops an electron deficiency on OH, which is evidenced by an increase in charge and a simultaneous decrease in volume. In addition, the observations in the study suggest that an OH+/electron transfer reaction can be classified as an adiabatic radical transfer, and the OH-/electron transfer reaction as a less adiabatic ion-coupled electron transfer.

Cytochrome P450 CYP736A12 is crucial for Trichoderma asperellum-induced alleviation of phoxim phytotoxicity and reduction of pesticide residue in tomato roots

Trichoderma can enhance the metabolism of organophosphate pesticides in plants, but the mechanism is unclear. Here, we performed high-throughput transcriptome sequencing of roots upon Trichoderma asperellum (TM) inoculation and phoxim (P) application in tomato (Solanum lycopersicum L.). A total of 4059 differentially expressed genes (DEGs) were obtained, including 2110 up-regulated and 1949 down-regulated DEGs in P vs TM+P. COG and KOG analysis indicated that DEGs were mainly enriched in signal transduction mechanisms. We then focused on the pesticide detoxification pathway and screened out cytochrome P450 CYP736A12 as a putative gene for functional analysis. We suppressed the expression of CYP736A12 in tomato plants by virus-induced gene silencing and analyzed tissue-specific phoxim residues, oxidative stress markers, glutathione pool, GST activity and related gene expression. Silencing CYP736A12 significantly increased phoxim residue and induced oxidative stress in tomato plants, by attenuating the TM-induced increased activity of antioxidant and detoxification enzymes, redox homeostasis and transcripts of detoxification genes including CYP724B2, GSH1, GSH2, GR, GPX, GST1, GST2, GST3, and ABC. The study revealed a critical mechanism by which TM promotes the metabolism of phoxim in tomato roots, which can be useful for further understanding the Trichoderma-induced xenobiotic detoxification and improving food safety.

Employing Genomic Tools to Explore the Molecular Mechanisms behind the Enhancement of Plant Growth and Stress Resilience Facilitated by a Burkholderia Rhizobacterial Strain

The rhizobacterial strain BJ3 showed 16S rDNA sequence similarity to species within the Burkholderia genus. Its complete genome sequence revealed a 97% match with Burkholderia contaminans and uncovered gene clusters essential for plant-growth-promoting traits (PGPTs). These clusters include genes responsible for producing indole acetic acid (IAA), osmolytes, non-ribosomal peptides (NRPS), volatile organic compounds (VOCs), siderophores, lipopolysaccharides, hydrolytic enzymes, and spermidine. Additionally, the genome contains genes for nitrogen fixation and phosphate solubilization, as well as a gene encoding 1-aminocyclopropane-1-carboxylate (ACC) deaminase. The treatment with BJ3 enhanced root architecture, boosted vegetative growth, and accelerated early flowering in Arabidopsis. Treated seedlings also showed increased lignin production and antioxidant capabilities, as well as notably increased tolerance to water deficit and high salinity. An RNA-seq transcriptome analysis indicated that BJ3 treatment significantly activated genes related to immunity induction, hormone signaling, and vegetative growth. It specifically activated genes involved in the production of auxin, ethylene, and salicylic acid (SA), as well as genes involved in the synthesis of defense compounds like glucosinolates, camalexin, and terpenoids. The expression of AP2/ERF transcription factors was markedly increased. These findings highlight BJ3's potential to produce various bioactive metabolites and its ability to activate auxin, ethylene, and SA signaling in Arabidopsis, positioning it as a new Burkholderia strain that could significantly improve plant growth, stress resilience, and immune function.

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

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

Genome-wide identification of the CYP82 gene family in cucumber and functional characterization of CsCYP82D102 in regulating resistance to powdery mildew

The cytochrome P450 (CYP450) gene family plays a vital role in basic metabolism, hormone signaling, and enhances plant resistance to stress. Among them, the CYP82 gene family is primarily found in dicots, and they are typically activated in response to various specific environmental stresses. Nevertheless, their roles remain considerably obscure, particularly within the context of cucumber. In the present study, 12 CYP82 subfamily genes were identified in the cucumber genome. Bioinformatics analysis included gene structure, conserved motif, cis-acting promoter element, and so on. Subcellular localization predicted that all CYP82 genes were located in the endoplasmic reticulum. The results of cis element analysis showed that CYP82s may significantly affect the response to stress, hormones, and light exposure. Expression patterns of the CYP82 genes were characterized by mining available RNA-seq data followed by qRT-PCR (quantitative real-time polymerase chain reaction) analysis. Members of CYP82 genes display specific expression profiles in different tissues, and in response to PM and abiotic stresses in this study, the role of CsCYP82D102, a member of the CYP82 gene family, was investigated. The upregulation of CsCYP82D102 expression in response to powdery mildew (PM) infection and treatment with methyl jasmonate (MeJA) or salicylic acid (SA) was demonstrated. Further research found that transgenic cucumber plants overexpressing CsCYP82D102 display heightened resistance against PM. Wild-type (WT) leaves exhibited average lesion areas of approximately 29.7% at 7 dpi upon powdery mildew inoculation. In contrast, the two independent CsCYP82D102 overexpression lines (OE#1 and OE#3) displayed significantly reduced necrotic areas, with average lesion areas of approximately 13.4% and 5.7%. Additionally, this enhanced resistance is associated with elevated expression of genes related to the SA/MeJA signaling pathway in transgenic cucumber plants. This study provides a theoretical basis for further research on the biological functions of the P450 gene in cucumber plants.

Genome-Wide Association Study and Identification of Candidate Genes Associated with Seed Number per Pod in Soybean

Soybean (Glycine max [L.] Merr.) is one of the primary sources of plant protein and oil for human foods, animal feed, and industrial processing. The seed number per pod generally varies from one to four and is an important component of seed number per unit area and seed yield. We used natural variation in 264 landraces and improved cultivars or lines to identify candidate genes involved in the regulation of seed number per pod in soybean. Genome-wide association tests revealed 65 loci that are associated with seed number per pod trait. Among them, 11 could be detected in multiple environments. Candidate genes were identified for seed number per pod phenotype from the most significantly associated loci, including a gene encoding protein argonaute 4, a gene encoding histone acetyltransferase of the MYST family 1, a gene encoding chromosome segregation protein SMC-1 and a gene encoding exocyst complex component EXO84A. In addition, plant hormones were found to be involved in ovule and seed development and the regulation of seed number per pod in soybean. This study facilitates the dissection of genetic networks underlying seed number per pod in soybean, which will be useful for the genetic improvement of seed yield in soybean.

Transcriptomic Insights and Cytochrome P450 Gene Analysis in Kadsura coccinea for Lignan Biosynthesis

Kadsura coccinea is a medicinal plant from the Schisandraceae family that is native to China and has great pharmacological potential due to its lignans. However, there are significant knowledge gaps regarding the genetic and molecular mechanisms of lignans. We used transcriptome sequencing technology to analyze root, stem, and leaf samples, focusing on the identification and phylogenetic analysis of Cytochrome P450 (CYP) genes. High-quality data containing 158,385 transcripts and 68,978 unigenes were obtained. In addition, 36,293 unigenes in at least one database, and 23,335 across five databases (Nr, KEGG, KOG, TrEMBL, and SwissProt) were successfully annotated. The KEGG pathway classification and annotation of these unigenes identified 10,825 categorized into major metabolic pathways, notably phenylpropanoid biosynthesis, which is essential for lignan synthesis. A key focus was the identification and phylogenetic analysis of 233 Cytochrome P450 (CYP) genes, revealing their distribution across 38 families in eight clans, with roots showing specific CYP gene expression patterns indicative of their role in lignan biosynthesis. Sequence alignment identified 22 homologous single genes of these CYPs, with 6 homologous genes of CYP719As and 1 of CYP81Qs highly expressed in roots. Our study significantly advances the understanding of the biosynthesis of dibenzocyclooctadiene lignans, offering valuable insights for future pharmacological research and development.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis

BACKGROUND

Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies.

RESULTS

A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins.

CONCLUSIONS

A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Mutation to a cytochrome P450 -like gene alters the leaf color by affecting the heme and chlorophyll biosynthesis pathways in Brassica napus

The regulated biosynthesis of chlorophyll is important because of its effects on plant photosynthesis and dry biomass production. In this study, a map-based cloning approach was used to isolate the cytochrome P450 -like gene BnaC08g34840D (BnCDE1) from a chlorophyll-deficient mutant (cde1) of Brassica napus obtained by ethyl methanesulfonate (EMS) mutagenization. Sequence analyses revealed that BnaC08g34840D in the cde1 mutant (BnCDE1I320T ) encodes a substitution at amino acid 320 (Ile320Thr) in the conserved region. The over-expression of BnCDE1I320T in ZS11 (i.e., gene-mapping parent with green leaves) recapitulated a yellow-green leaf phenotype. The CRISPR/Cas9 genome-editing system was used to design two single-guide RNAs (sgRNAs) targeting BnCDE1I320T in the cde1 mutant. The knockout of BnCDE1I320T in the cde1 mutant via a gene-editing method restored normal leaf coloration (i.e., green leaves). These results indicate that the substitution in BnaC08g34840D alters the leaf color. Physiological analyses showed that the over-expression of BnCDE1I320T leads to decreases in the number of chloroplasts per mesophyll cell and in the contents of the intermediates of the chlorophyll biosynthesis pathway in leaves, while it increases heme biosynthesis, thereby lowering the photosynthetic efficiency of the cde1 mutant. The Ile320Thr mutation in the highly conserved region of BnaC08g34840D inhibited chlorophyll biosynthesis and disrupted the balance between heme and chlorophyll biosynthesis. Our findings may further reveal how the proper balance between the chlorophyll and heme biosynthesis pathways is maintained.

Genome-based identification of the CYP75 gene family in Orchidaceae and its expression patterns in Cymbidium goeringii

With a great diversity of species, Orchidaceae stands out as an essential component of plant biodiversity, making it a primary resource for studying angiosperms evolution and genomics. This study focuses on 13 published orchid genomes to identify and analyze the CYP75 gene family belonging to the cytochrome P450 superfamily, which is closely related to flavonoid biosynthetic enzymes and pigment regulation. We found 72 CYP75s in the 13 orchid genomes and further classified them into two classes: CYP75A and CYP75B subfamily, the former synthesizes blue anthocyanins, while the latter is involved in the production of red anthocyanins. Furthermore, the amount of CYP75Bs (53/72) greatly exceeds the amount of CYP75As (19/72) in orchids. Our findings suggest that CYP75B genes have a more important evolutionary role, as red plants are more common in nature than blue plants. We also discovered unique conserved motifs in each subfamily that serve as specific recognition features (motif 19 belong to CYP75A; motif 17 belong to CYP75B). Two diverse-colored varieties of C. goeringii were selected for qRT-PCR experiments. The expression of CgCYP75B1 was significantly higher in the purple-red variant compared to the yellow-green variant, while CgCYP75A1 showed no significant difference. Based on transcriptomic expression analysis, CYP75Bs are more highly expressed than CYP75As in floral organs, especially in colorful petals and lips. These results provide valuable information for future studies on CYP75s in orchids and other angiosperms.

Genome-wide analysis of the P450 gene family in tea plant (Camellia sinensis) reveals functional diversity in abiotic stress

BACKGROUND

Cytochrome P450 (Cytochrome P450s) genes are involved in the catalysis of various reactions, including growth, development, and secondary metabolite biosynthetic pathways. However, little is known about the characteristics and functions of the P450 gene family in Camellia sinensis (C. sinensis).

RESULTS

To reveal the mechanisms of tea plant P450s coping with abiotic stresses, analyses of the tea plant P450 gene family were conducted using bioinformatics-based methods. In total, 273 putative P450 genes were identified from the genome database of C. sinensis. The results showed that P450s were well-balanced across the chromosomes I to XV of entire genome, with amino acid lengths of 268-612 aa, molecular weights of 30.95-68.5 kDa, and isoelectric points of 4.93-10.17. Phylogenetic analysis divided CsP450s into 34 subfamilies, of which CYP71 was the most abundant. The predicted subcellular localization results showed that P450 was distributed in a variety of organelles, with chloroplasts, plasma membrane,,and cytoplasm localized more frequently. The promoter region of CsP450s contained various cis-acting elements related to phytohormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) were identified in the CsP450 family proteins, with 27 genes lacking introns and only one exon. The results of genome large segment duplication showed that there were 37 pairs of genes with tandem duplication. Interaction network analysis showed that CsP450 could interact with multiple types of target genes, and there are protein interactions within the family. Tissue expression analysis showed that P450 was highly expressed in roots and stems. Moreover, qPCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results.

CONCLUSIONS

This study lays the foundation for resolving the classification and functional study of P450 family genes and provides a reference for the molecular breeding of C. sinensis.

A Comprehensive Genome-Wide Investigation of the Cytochrome 71 (OsCYP71) Gene Family: Revealing the Impact of Promoter and Gene Variants (Ser33Leu) of OsCYP71P6 on Yield-Related Traits in Indica Rice (Oryza sativa L.)

The cytochrome P450 (CYP450) gene family plays a critical role in plant growth and developmental processes, nutrition, and detoxification of xenobiotics in plants. In the present research, a comprehensive set of 105 OsCYP71 family genes was pinpointed within the genome of indica rice. These genes were categorized into twelve distinct subfamilies, where members within the same subgroup exhibited comparable gene structures and conserved motifs. In addition, 105 OsCYP71 genes were distributed across 11 chromosomes, and 36 pairs of OsCYP71 involved in gene duplication events. Within the promoter region of OsCYP71, there exists an extensive array of cis-elements that are associated with light responsiveness, hormonal regulation, and stress-related signaling. Further, transcriptome profiling revealed that a majority of the genes exhibited responsiveness to hormones and were activated across diverse tissues and developmental stages in rice. The OsCYP71P6 gene is involved in insect resistance, senescence, and yield-related traits in rice. Hence, understanding the association between OsCYP71P6 genetic variants and yield-related traits in rice varieties could provide novel insights for rice improvement. Through the utilization of linear regression models, a total of eight promoters were identified, and a specific gene variant (Ser33Leu) within OsCYP71P6 was found to be linked to spikelet fertility. Additionally, different alleles of the OsCYP71P6 gene identified through in/dels polymorphism in 131 rice varieties were validated for their allelic effects on yield-related traits. Furthermore, the single-plant yield, spikelet number, panicle length, panicle weight, and unfilled grain per panicle for the OsCYP71P6-1 promoter insertion variant were found to contribute 20.19%, 13.65%, 5.637%, 8.79%, and 36.86% more than the deletion variant, respectively. These findings establish a robust groundwork for delving deeper into the functions of OsCYP71-family genes across a range of biological processes. Moreover, these findings provide evidence that allelic variation in the promoter and amino acid substitution of Ser33Leu in the OsCYP71P6 gene could potentially impact traits related to rice yield. Therefore, the identified promoter variants in the OsCYP71P6 gene could be harnessed to amplify rice yields.

Integrated analysis of metabolome and transcriptome reveals key candidate genes involved in flavonoid biosynthesis in Pinellia ternata under heat stress

Pinellia ternata (Thunb.) Breit. is an important traditional Chinese medicinal herb and very sensitive to high temperatures. To gain a better understanding of flavonoid biosynthesis under heat stress in P. ternata, we performed integrated analyses of metabolome and transcriptome data. P. ternata plants were subjected to a temperature of 38 °C, and samples were collected after 10 d of treatment. A total of 502 differential accumulated metabolites and 5040 different expressed transcripts were identified, with flavonoid biosynthesis predominantly enriched. Integrated metabolomics and transcriptome analysis showed that high temperature treatment upregulated the expression of CYP73A and downregulated the expression of other genes (such as HCT, CCoAOMT, DFR1, DFR2), which might inhibit the biosynthesis of the downstream metabolome, including such metabolites as chlorogenic acid, pelargonidin, cyanidin, and (-)-epigallocatechin in the flavonoid biosynthesis pathway. The transcription expression levels of these genes were validated by real-time PCR. Our results provide valuable insights into flavonoid composition and accumulation patterns and the candidate genes participating in the flavonoid biosynthesis pathways under heat stress in P. ternata.

Genome-wide investigation of Cytochrome P450 superfamily of Aquilaria agallocha: Association with terpenoids and phenylpropanoids biosynthesis

Agarwood is a dark resinous wood, produced when Aquilaria tree responds to wounding and microbial infection resulting in the accumulation of fragrant metabolites. Sesquiterpenoids and 2-(2-phenylethyl) chromones are the major phytochemicals in agarwood and Cytochrome P450s (CYPs) are one of the important enzymes in the biosynthesis of these fragrant chemicals. Thus, understanding the repertoire of CYP superfamily in Aquilaria can not only give insights into the fundamentals of agarwood formation, but can also provide a tool for the overproduction of the aroma chemicals. Therefore, current study was designed to investigate CYPs of an agarwood producing plant, Aquilaria agallocha. We identified 136 CYP genes from A. agallocha genome (AaCYPs) and classified them into 8 clans and 38 families. The promoter regions had stress and hormone-related cis-regulatory elements which indicate their participation in the stress response. Duplication and synteny analysis revealed segmental and tandem duplicated and evolutionary related CYP members in other plants. Potential members involved in the biosynthesis of sesquiterpenoids and phenylpropanoids were identified and found to be upregulated in methyl jasmonate-induced callus and infected Aquilaria trees by real-time quantitative PCR analyses. This study highlights the possible involvement of AaCYPs in agarwood resin development and their complex regulation during stress exposure.

Physiological and transcriptome analysis of Dendrobium officinale under low nitrogen stress

Nitrogen (N) is the main nutrient of plants, and low nitrogen usually affects plant growth and crop yield. The traditional Chinese herbal medicine Dendrobium officinale Kimura et. Migo is a typical low nitrogen-tolerant plant, and its mechanism in response to low nitrogen stress has not previously been reported. In this study, physiological measurements and RNA-Seq analysis were used to analyse the physiological changes and molecular responses of D. officinale under different nitrogen concentrations. The results showed that under low nitrogen levels, the growth, photosynthesis and superoxide dismutase activity were found to be significantly inhibited, while the activities of peroxidase and catalase, the content of polysaccharides and flavonoids significantly increased. Differentially expressed genes (DEGs) analysis showed that nitrogen and carbon metabolisms, transcriptional regulation, antioxidative stress, secondary metabolite synthesis and signal transduction all made a big difference in low nitrogen stress. Therefore, copious polysaccharide accumulation, efficient assimilation and recycling of nitrogen, as well as rich antioxidant components play critical roles. This study is helpful for understanding the response mechanism of D. officinale to low nitrogen levels, which might provide good guidance for practical production of high quality D. officinale .

Transcriptomic profiling reveals candidate allelopathic genes in rice responsible for interactions with barnyardgrass

Echinochloa crus-galli (barnyardgrass) is one of the most damaging weeds in rice fields worldwide. Allelopathy has been considered a possible application for weed management. Thus understanding its molecular mechanisms is important for rice production. This study generated transcriptomes from rice under mono- and co-culture with barnyardgrass at two-time points to identify the candidate genes controlling allelopathic interactions between rice and barnyardgrass. A total of 5,684 differentially expressed genes (DEGs) were detected, amongst which 388 genes were transcription factors. These DEGs include genes associated with momilactone and phenolic acid biosynthesis, which play critical roles in allelopathy. Additionally, we found significantly more DEGs at 3 hours than at 3 days, suggesting a quick allelopathic response in rice. Up-regulated DEGs involve diverse biological processes, such as response to stimulus and pathways related to phenylpropanoid and secondary metabolites biosynthesis. Down-regulated DEGs were involved in developmental processes, indicating a balance between growth and stress response to allelopathy from barnyardgrass. Comparison of DEGs between rice and barnyardgrass shows few common genes, suggesting different mechanisms underlying allelopathic interaction in these two species. Our results offer an important basis for identifying of candidate genes responsible for rice and barnyardgrass interactions and contribute valuable resources for revealing its molecular mechanisms.

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.

Identification of a Major Locus for Lodging Resistance to Typhoons Using QTL Analysis in Rice

We detected a new target quantitative trait locus (QTL) for lodging resistance in rice by analyzing lodging resistance to typhoons (Maysak and Haishen) using a scale from 0 (no prostrating) to 1 (little prostrating or prostrating) to record the resistance score in a Cheongcheong/Nagdong double haploid rice population. Five quantitative trait loci for lodging resistance to typhoons were detected. Among them, qTyM6 and qTyH6 exhibited crucial effects of locus RM3343-RM20318 on chromosome 6, which overlaps with our previous rice lodging studies for the loci qPSLSA6-2, qPSLSB6-5, and qLTI6-2. Within the target locus RM3343-RM20318, 12 related genes belonging to the cytochrome P450 protein family were screened through annotation. Os06g0599200 (OsTyM/Hq6) was selected for further analysis. We observed that the culm and panicle lengths were positively correlated with lodging resistance to typhoons. However, the yield was negatively correlated with lodging resistance to typhoons. The findings of this study improve an understanding of rice breeding, particularly the culm length, early maturing, and heavy panicle varieties, and the mechanisms by which the plant's architecture can resist natural disasters such as typhoons to ensure food safety. These results also provide the insight that lodging resistance in rice may be associated with major traits such as panicle length, culm length, tiller number, and heading date, and thereby improvements in these traits can increase lodging resistance to typhoons. Moreover, rice breeding should focus on maintaining suitable varieties that can withstand the adverse effects of climate change in the future and provide better food security.

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.

Genome-Wide Identification of the Cytochrome P450 Superfamily Genes and Targeted Editing of BnCYP704B1 Confers Male Sterility in Rapeseed

The cytochrome P450 (CYP450) monooxygenase superfamily, which is involved in the biosynthesis pathways of many primary and secondary metabolites, plays prominent roles in plant growth and development. However, systemic information about CYP450s in Brassica napus (BnCYP450) was previously undiscovered and their biological significance are far from understood. Members of clan 86 CYP450s, such as CYP704Bs, are essential for the formation of pollen exine in plant male reproduction, and the targeted mutagenesis of CYP704B genes has been used to create new male sterile lines in many crops. In the present study, a total of 687 BnCYP450 genes were identified in Brassica napus cultivar "Zhongshuang 11" (ZS11), which has nearly 2.8-fold as many CYP450 members as in Arabidopsis thaliana. It is rationally estimated since Brassica napus is a tetraploid oil plant with a larger genome compared with Arabidopsis thaliana. The BnCYP450 genes were divided into 47 subfamilies and clustered into nine clans. Phylogenetic relationship analysis reveals that CYP86 clan consists of four subfamilies and 109 BnCYP450s. Members of CYP86 clan genes display specific expression profiles in different tissues and in response to ABA and abiotic stresses. Two BnCYP450s within the CYP704 subfamily from CYP86 clan, BnCYP704B1a and BnCYP704B1b, display high similarity to MS26 (Male Sterility 26, also known as CYP704B1). These two BnCYP704B1 genes were specifically expressed in young buds. We then simultaneously knocked-out these two BnCYP704B1 genes through a clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) genome engineering system. The edited plants displayed a pollenless, sterile phenotype in mature anthers, suggesting that we successfully reproduced genic male sterility (GMS, also known as nuclear male sterility) lines in Brassica napus. This study provides a systemic view of BnCYP450s and offers a strategy to facilitate the commercial utility of the CRISPR/Cas9 system for the rapid generation of GMS in rapeseed via knocking-out GMS controlling genes.

Single-cell RNA sequencing provides a high-resolution roadmap for understanding the multicellular compartmentation of specialized metabolism

Monoterpenoid indole alkaloids (MIAs) are among the most diverse specialized metabolites in plants and are of great pharmaceutical importance. We leveraged single-cell transcriptomics to explore the spatial organization of MIA metabolism in Catharanthus roseus leaves, and the transcripts of 20 MIA genes were first localized, updating the model of MIA biosynthesis. The MIA pathway was partitioned into three cell types, consistent with the results from RNA in situ hybridization experiments. Several candidate transporters were predicted to be essential players shuttling MIA intermediates between inter- and intracellular compartments, supplying potential targets to increase the overall yields of desirable MIAs in native plants or heterologous hosts through metabolic engineering and synthetic biology. This work provides not only a universal roadmap for elucidating the spatiotemporal distribution of biological processes at single-cell resolution, but also abundant cellular and genetic resources for further investigation of the higher-order organization of MIA biosynthesis, transport and storage.

The protein conformational basis of isoflavone biosynthesis

Isoflavonoids play important roles in plant defense and also exhibit a range of mammalian health-promoting activities. Their biosynthesis is initiated by two enzymes with unusual catalytic activities; 2-hydroxyisoflavanone synthase (2-HIS), a membrane-bound cytochrome P450 catalyzing a coupled aryl-ring migration and hydroxylation, and 2-hydroxyisoflavanone dehydratase (2-HID), a member of a large carboxylesterase family that paradoxically catalyzes dehydration of 2-hydroxyisoflavanones to isoflavone. Here we report the crystal structures of 2-HIS from Medicago truncatula and 2-HID from Pueraria lobata. The 2-HIS structure reveals a unique cytochrome P450 conformation and heme and substrate binding mode that facilitate the coupled aryl-ring migration and hydroxylation reactions. The 2-HID structure reveals the active site architecture and putative catalytic residues for the dual dehydratase and carboxylesterase activities. Mutagenesis studies revealed key residues involved in substrate binding and specificity. Understanding the structural basis of isoflavone biosynthesis will facilitate the engineering of new bioactive isoflavonoids.

The structure and function of two isoflavone biosynthetic enzymes are reported revealing a novel cytochrome P450 conformation and identification of key residues for dual dehydratase and carboxylesterase activities

Manipulating GA-Related Genes for Cereal Crop Improvement

The global population is projected to experience a rapid increase in the future, which poses a challenge to global food sustainability. The "Green Revolution" beginning in the 1960s allowed grain yield to reach two billion tons in 2000 due to the introduction of semi-dwarfing genes in cereal crops. Semi-dwarfing genes reduce the gibberellin (GA) signal, leading to short plant stature, which improves the lodging resistance and harvest index under modern fertilization practices. Here, we reviewed the literature on the function of GA in plant growth and development, and the role of GA-related genes in controlling key agronomic traits that contribute to grain yield in cereal crops. We showed that: (1) GA is a significant phytohormone in regulating plant development and reproduction; (2) GA metabolism and GA signalling pathways are two key components in GA-regulated plant growth; (3) GA interacts with other phytohormones manipulating plant development and reproduction; and (4) targeting GA signalling pathways is an effective genetic solution to improve agronomic traits in cereal crops. We suggest that the modification of GA-related genes and the identification of novel alleles without a negative impact on yield and adaptation are significant in cereal crop breeding for plant architecture improvement. We observed that an increasing number of GA-related genes and their mutants have been functionally validated, but only a limited number of GA-related genes have been genetically modified through conventional breeding tools and are widely used in crop breeding successfully. New genome editing technologies, such as the CRISPR/Cas9 system, hold the promise of validating the effectiveness of GA-related genes in crop development and opening a new venue for efficient and accelerated crop breeding.

Plant Specialised Glycosides (PSGs): their biosynthetic enzymatic machinery, physiological functions and commercial potential

Plants, the primary producers of our planet, have evolved from simple aquatic life to very complex terrestrial habitat. This habitat transition coincides with evolution of enormous chemical diversity, collectively termed as 'Plant Specialised Metabolisms (PSMs)', to cope the environmental challenges. Plant glycosylation is an important process of metabolic diversification of PSMs to govern their in planta stability, solubility and inter/intra-cellular transport. Although, individual category of PSMs (terpenoids, phenylpropanoids, flavonoids, saponins, alkaloids, phytohormones, glucosinolates and cyanogenic glycosides) have been well studied; nevertheless, deeper insights of physiological functioning and genomic aspects of plant glycosylation/deglycosylation processes including enzymatic machinery (CYPs, GTs, and GHs) and regulatory elements are still elusive. Therefore, this review discussed the paradigm shift on genomic background of enzymatic machinery, transporters and regulatory mechanism of 'Plant Specialised Glycosides (PSGs)'. Current efforts also update the fundamental understanding about physiological, evolutionary and adaptive role of glycosylation/deglycosylation processes during the metabolic diversification of PSGs. Additionally, futuristic considerations and recommendations for employing integrated next-generation multi-omics (genomics, transcriptomics, proteomics and metabolomics), including gene/genome editing (CRISPR-Cas) approaches are also proposed to explore commercial potential of PSGs.

Melatonin enhances stress tolerance in pigeon pea by promoting flavonoid enrichment, particularly luteolin in response to salt stress

Melatonin improves plant resistance to multiple stresses by participating in the biosynthesis of metabolites. Flavonoids are an important family of plant secondary metabolites and are widely recognized to be involved in resistance; however, the crosstalk between melatonin and flavonoid is largely unknown. We found that the resistance of pigeon pea (Cajanus cajan) to salt, drought, and heat stresses were significantly enhanced by pre-treatment with melatonin. Combined transcriptome and LC-ESI-MS/MS metabolomics analyses showed that melatonin significantly induced the enrichment of flavonoids and mediated the reprogramming of biosynthetic pathway genes. The highest fold-increase in expression in response to melatonin treatment was observed for the CcF3´H family, which encodes an enzyme that catalyses the biosynthesis of luteolin, and the transcription factor CcPCL1 directly bonded to the CcF3´H-5 promoter to enhance its expression. In addition, salt stress also induced the expression of CcPCL1 and CcF3´H-5, and their overexpression in transgenic plants greatly enhanced salt tolerance by promoting the biosynthesis of luteolin. Overall, our results indicated that pre-treatment of pigeon pea with melatonin promoted luteolin biosynthesis through the CcPCL1 and CcF3´H-5 pathways, resulting in salt tolerance. Our study shows that melatonin enhances plant tolerance to multiple stresses by mediating flavonoid biosynthesis, providing new avenues for studying the crosstalk between melatonin and flavonoids.

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.

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.

Chromosome-scale genome assembly of Rhododendron molle provides insights into its evolution and terpenoid biosynthesis

BACKGROUND

Rhododendron molle (Ericaceae) is a traditional Chinese medicine, which has been used to treat rheumatism and relieve pain since ancient times. The characteristic grayanoids of this plant have been demonstrated to be the chemical basis for the analgesic activity. Moreover, unlike morphine, these diterpenoids are non-addictive. Grayanoids mainly distribute in the leaves, flowers, roots, and fruits of R. molle, with low content. Currently the research on the biosynthesis of grayanoids is hindered, partially due to lack of the genomic information.

RESULTS

In the present study, a total of 744 Mb sequences were generated and assembled into 13 chromosomes. An ancient whole-genome duplication event (Ad-β) was discovered that occurred around 70 million years ago. Tandem and segmental gene duplications led to specific gene expansions in the terpene synthase and cytochrome P450 (CYP450) gene families. Two diterpene synthases were demonstrated to be responsible for the biosynthesis of 16α-hydroxy-ent-kaurane, the key precursor for grayanoids. Phylogenetic analysis revealed a species-specific bloom of the CYP71AU subfamily, which may involve the candidate CYP450s responsible for the biosynthesis of grayanoids. Additionally, three putative terpene biosynthetic gene clusters were found.

CONCLUSIONS

We reported the first genome assembly of R. molle and investigated the molecular basis underpinning terpenoids biosynthesis. Our work provides a foundation for elucidating the complete biosynthetic pathway of grayanoids and studying the terpenoids diversity in R. molle.

APIP5 functions as a transcription factor and an RNA-binding protein to modulate cell death and immunity in rice

Many transcription factors (TFs) in animals bind to both DNA and mRNA, regulating transcription and mRNA turnover. However, whether plant TFs function at both the transcriptional and post-transcriptional levels remains unknown. The rice (Oryza sativa) bZIP TF AVRPIZ-T-INTERACTING PROTEIN 5 (APIP5) negatively regulates programmed cell death and blast resistance and is targeted by the effector AvrPiz-t of the blast fungus Magnaporthe oryzae. We demonstrate that the nuclear localization signal of APIP5 is essential for APIP5-mediated suppression of cell death and blast resistance. APIP5 directly targets two genes that positively regulate blast resistance: the cell wall-associated kinase gene OsWAK5 and the cytochrome P450 gene CYP72A1. APIP5 inhibits OsWAK5 expression and thus limits lignin accumulation; moreover, APIP5 inhibits CYP72A1 expression and thus limits reactive oxygen species production and defense compounds accumulation. Remarkably, APIP5 acts as an RNA-binding protein to regulate mRNA turnover of the cell death- and defense-related genes OsLSD1 and OsRac1. Therefore, APIP5 plays dual roles, acting as TF to regulate gene expression in the nucleus and as an RNA-binding protein to regulate mRNA turnover in the cytoplasm, a previously unidentified regulatory mechanism of plant TFs at the transcriptional and post-transcriptional levels.

Genome-wide transcriptome analysis and characterization of the cytochrome P450 flavonoid biosynthesis genes in pigeon pea (Cajanus cajan)

226 CcCYP450 genes were identified at the genomic level and were classified into 45 clades based on phylogenetic analysis. CcCYP75B165 gene was found that might play important roles in the biosynthesis of flavonoids in pigeon pea, and was significantly induced by methyl jasmonate (MeJA). The cytochrome P450 mono-oxygenase (CYP450) superfamily plays a key role in the flavonoid biosynthesis pathway and resists different kinds of stresses. Several CYP450 genes have been identified to be involved in the biosynthesis of crop protection agents. However, the CcCYP450 genes from pigeon pea have not been identified. Here, 226 CcCYP450 genes were identified at the genomic level by analysing the gene structure, distribution on chromosomes, gene duplication, and conserved motifs and were classified into 45 clades based on phylogenetic analysis. RNA-seq analysis revealed clear details of CcCYP450 genes that varied with time of MeJA (methyl jasmonate) induction. Among them, six CcCYP450 subfamily genes were found that might play important roles in the biosynthesis of flavonoids in pigeon pea. The overexpression of CcCYP75B165 in pigeon pea significantly induced the accumulation of genistin and downregulated the contents of cajaninstilbene acid, apigenin, isovitexin, and genistein and the expression of flavonoid synthase genes. This study provides theoretical guidance and plant genetic resources for cultivating new pigeon pea varieties with high flavonoid contents and exploring the molecular mechanisms of the biosynthesis of flavonoids under MeJA treatment.

Enhanced reactive oxygen detoxification occurs in salt-stressed soybean roots expressing GmSALT3

Soybean (Glycine max) is an important crop globally for food and edible oil production. Soybean plants are sensitive to salinity (NaCl), with significant yield decreases reported under saline conditions. GmSALT3 is the dominant gene underlying a major QTL for salt tolerance in soybean. GmSALT3 encodes a transmembrane protein belonging to the plant cation/proton exchanger (CHX) family, and is predominately expressed in root phloem and xylem associated cells under both saline and non-saline conditions. It is currently unknown through which molecular mechanism(s) the ER-localised GmSALT3 contributes to salinity tolerance, as its localisation excludes direct involvement in ion exclusion. In order to gain insights into potential molecular mechanism(s), we used RNA-seq analysis of roots from two soybean NILs (near isogenic lines); NIL-S (salt-sensitive, Gmsalt3), and NIL-T (salt-tolerant, GmSALT3), grown under control and saline conditions (200 mM NaCl) at three time points (0 h, 6 h, and 3 days). Gene ontology (GO) analysis showed that NIL-T has greater responses aligned to oxidation reduction. ROS were less abundant and scavenging enzyme activity was greater in NIL-T, consistent with the RNA-seq data. Further analysis indicated that genes related to calcium signalling, vesicle trafficking and Casparian strip (CS) development were upregulated in NIL-T following salt treatment. We propose that GmSALT3 improves the ability of NIL-T to cope with saline stress through preventing ROS overaccumulation in roots, and potentially modulating Ca2+ signalling, vesicle trafficking and formation of diffusion barriers.

Biochemical Investigation of Membrane-Bound Cytochrome b5 and the Catalytic Domain of Cytochrome b5 Reductase from Arabidopsis thaliana

The endoplasmic reticulum (ER) membrane of plant cells contains several enzymes responsible for the biosynthesis of a diverse range of molecules essential for plant growth and holds potential for industrial applications. Many of these enzymes are dependent on electron transfer proteins to sustain their catalytic cycles. In plants, two crucial ER-bound electron transfer proteins are cytochrome b5 and cytochrome b5 reductase, which catalyze the stepwise transfer of electrons from NADH to redox enzymes such as fatty acid desaturases, cytochrome P450s, and plant aldehyde decarbonylase. Despite the high significance of plant cytochrome b5 and cytochrome b5 reductase, they have eluded detailed characterization to date. Here, we overexpressed the full-length membrane-bound cytochrome b5 isoform B from the model plant Arabidopsis thaliana in Escherichia coli, purified the protein employing detergents as well as styrene-maleic acid (SMA) copolymers, and biochemically characterized the protein. The SMA-encapsulated cytochrome b5 exhibits a discoidal shape and the characteristic features of the active heme-bound state. We also overexpressed and purified the soluble domain of cytochrome b5 reductase from A. thaliana, establishing its activity, stability, and kinetic parameters. Further, we demonstrated that the plant cytochrome b5, purified in detergents and styrene maleic acid lipid particles (SMALPs), readily accepts electrons from the cognate plant cytochrome b5 reductase and distant electron mediators such as plant NADPH-cytochrome P450 oxidoreductase and cyanobacterial NADPH-ferredoxin reductase. We also measured the kinetic parameters of cytochrome b5 reductase for cytochrome b5. Our studies are the first to report the purification and detailed biochemical characterization of the plant cytochrome b5 and cytochrome b5 reductase from the bacterial overexpression system.

TaKLU Plays as a Time Regulator of Leaf Growth via Auxin Signaling

The growth of leaves is subject to strict time regulation. Several genes influencing leaf growth have been identified, but little is known about how genes regulate the orderly initiation and growth of leaves. Here, we demonstrate that TaKLU/TaCYP78A5 contributes to a time regulation mechanism in leaves from initiation to expansion. TaKLU encodes the cytochrome P450 CYP78A5, and its homolog AtKLU has been described whose deletion is detrimental to organ growth. Our results show that TaKLU overexpression increases leaf size and biomass by altering the time of leaf initiation and expansion. TaKLU-overexpressing plants have larger leaves with more cells. Further dynamic observations indicate that enlarged wheat leaves have experienced a longer expansion time. Different from AtKLU inactivation increases leaf number and initiation rates, TaKLU overexpression only smooths the fluctuations of leaf initiation rates by adjusting the initiation time of local leaves, without affecting the overall leaf number and initiation rates. In addition, complementary analyses suggest TaKLU is functionally conserved with AtKLU in controlling the leaf initiation and size and may involve auxin accumulation. Our results provide a new insight into the time regulation mechanisms of leaf growth in wheat.

Cytochrome P450 CYP709C56 metabolizing mesosulfuron-methyl confers herbicide resistance in Alopecurus aequalis

Multiple herbicide resistance in diverse weed species endowed by enhanced herbicide detoxification or degradation is rapidly growing into a great threat to herbicide sustainability and global food safety. Although metabolic resistance is frequently documented in the economically damaging arable weed species shortawn foxtail (Alopecurus aequalis Sobol.), relevant molecular knowledge has been lacking. Previously, we identified a field population of A. aequalis (R) that had evolved metabolic resistance to the commonly used acetolactate synthase (ALS)-inhibiting herbicide mesosulfuron-methyl. RNA sequencing was used to discover potential herbicide metabolism-related genes, and four cytochrome P450s (CYP709C56, CYP71R18, CYP94C117, and CYP94E14) were identified with higher expressions in the R vs. susceptible (S) plants. Here the full-length P450 complementary DNA transcripts were each cloned with identical sequences between the S and R plants. Transgenic Arabidopsis overexpressing CYP709C56 became resistant to the sulfonylurea herbicide mesosulfuron-methyl and the triazolo-pyrimidine herbicide pyroxsulam. This resistance profile generally but does not completely in accordance with what is evident in the R A. aequalis. Transgenic lines exhibited enhanced capacity for detoxifying mesosulfuron-methyl into O-demethylated metabolite, which is in line with the detection of O-demethylated herbicide metabolite in vitro in transformed yeast. Structural modeling predicted that mesosulfuron-methyl binds to CYP709C56 involving amino acid residues Thr-328, Thr-500, Asn-129, Gln-392, Phe-238, and Phe-242 for achieving O-demethylation. Constitutive expression of CYP709C56 was highly correlated with the metabolic mesosulfuron-methyl resistance in A. aequalis. These results indicate that CYP709C56 degrades mesosulfuron-methyl and its up-regulated expression in A. aequalis confers resistance to mesosulfuron-methyl.

Genome-Wide Identification and Expression Profiles of 13 Key Structural Gene Families Involved in the Biosynthesis of Rice Flavonoid Scaffolds

Flavonoids are a class of key polyphenolic secondary metabolites with broad functions in plants, including stress defense, growth, development and reproduction. Oryza sativa L. (rice) is a well-known model plant for monocots, with a wide range of flavonoids, but the key flavonoid biosynthesis-related genes and their molecular features in rice have not been comprehensively and systematically characterized. Here, we identified 85 key structural gene candidates associated with flavonoid biosynthesis in the rice genome. They belong to 13 families potentially encoding chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol synthase (FLS), leucoanthocyanidin dioxygenase (LDOX), anthocyanidin synthase (ANS), flavone synthase II (FNSII), flavanone 2-hydroxylase (F2H), flavonoid 3′-hydroxylase (F3′H), flavonoid 3′,5′-hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR). Through structural features, motif analyses and phylogenetic relationships, these gene families were further grouped into five distinct lineages and were examined for conservation and divergence. Subsequently, 22 duplication events were identified out of a total of 85 genes, among which seven pairs were derived from segmental duplication events and 15 pairs were from tandem duplications, demonstrating that segmental and tandem duplication events play important roles in the expansion of key flavonoid biosynthesis-related genes in rice. Furthermore, these 85 genes showed spatial and temporal regulation in a tissue-specific manner and differentially responded to abiotic stress (including six hormones and cold and salt treatments). RNA-Seq, microarray analysis and qRT-PCR indicated that these genes might be involved in abiotic stress response, plant growth and development. Our results provide a valuable basis for further functional analysis of the genes involved in the flavonoid biosynthesis pathway in rice.

Functional and structural insight into the flexibility of cytochrome P450 reductases from Sorghum bicolor and its implications for lignin composition

Plant NADPH-dependent cytochrome P450 reductase (CPR) is a multidomain enzyme that donates electrons for hydroxylation reactions catalyzed by class II cytochrome P450 monooxygenases involved in the synthesis of many primary and secondary metabolites. These P450 enzymes include trans-cinnamate-4-hydroxylase, p-coumarate-3′-hydroxylase, and ferulate-5-hydroxylase involved in monolignol biosynthesis. Because of its role in monolignol biosynthesis, alterations in CPR activity could change the composition and overall output of lignin. Therefore, to understand the structure and function of three CPR subunits from sorghum, recombinant subunits SbCPR2a, SbCPR2b, and SbCPR2c were subjected to X-ray crystallography and kinetic assays. Steady-state kinetic analyses demonstrated that all three CPR subunits supported the oxidation reactions catalyzed by SbC4H1 (CYP73A33) and SbC3′H (CYP98A1). Furthermore, comparing the SbCPR2b structure with the well-investigated CPRs from mammals enabled us to identify critical residues of functional importance and suggested that the plant flavin mononucleotide–binding domain might be more flexible than mammalian homologs. In addition, the elucidated structure of SbCPR2b included the first observation of NADP⁺ in a native CPR. Overall, we conclude that the connecting domain of SbCPR2, especially its hinge region, could serve as a target to alter biomass composition in bioenergy and forage sorghums through protein engineering.