Involvement of chalcone reductase in the soybean isoflavone metabolon: identification of GmCHR5, which interacts with 2-hydroxyisoflavanone synthase

Structures, activities, and putative biosynthetic pathways of characteristic polyphenolic compounds from Morus plants: A review

Morus plants played a pivotal role in ancient Chinese sericulture and silk production, which served as critical components of economy and culture. Besides, many parts of mulberry trees, including roots, leaves, stems, and fruits, hold various medicinal value, and have been utilized in traditional medicine for thousands of years. The chemical composition of mulberry has been reported in many literatures, while the characteristic compounds have not been systematically summarized. In this review, we focused on the polyphenolic compounds in mulberry, including flavonoids, 2-arylbenzofurans, and Diels-Alder (D-A) adducts, and summarized their structural features, structure-activity relationships, and potential biosynthetic pathways. The results revealed a characteristic class of 2'-hydroxylated flavonoids and stilbenes which played an important role in the biosynthesis of downstream 2-arylbenzofurans and D-A adducts in mulberry but had been overlooked by most studies. The prenylated modifications of different compounds were also discussed and their function as precursors of D-A adducts was emphasized. We also describe the effects of different modifications on biological activities. Besides, the chemical composition of Morus was most similar to that of Artocarpus in the Moraceae family in that they had almost identical characteristic compounds. Finally, a putative total biosynthetic pathway of D-A adducts in mulberry was proposed based on structure derivation and combination of verified reactions. This review contributes to the understanding of the biological activity and biosynthesis of the characteristic components of Morus plants.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

PagARGOS promotes low-lignin wood formation in poplar

Wood formation, which occurs mainly through secondary xylem development, is important not only for supplying raw material for the 'ligno-chemical' industry but also for driving the storage of carbon. However, the complex mechanisms underlying the promotion of xylem formation remain to be elucidated. Here, we found that overexpression of Auxin-Regulated Gene involved in Organ Size (ARGOS) in hybrid poplar 84 K (Populus alba × Populus tremula var. glandulosa) enlarged organ size. In particular, PagARGOS promoted secondary growth of stems with increased xylem formation. To gain further insight into how PagARGOS regulates xylem development, we further carried out yeast two-hybrid screening and identified that the auxin transporter WALLS ARE THIN1 (WAT1) interacts with PagARGOS. Overexpression of PagARGOS up-regulated WAT1, activating a downstream auxin response promoting cambial cell division and xylem differentiation for wood formation. Moreover, overexpressing PagARGOS caused not only higher wood yield but also lower lignin content compared with wild-type controls. PagARGOS is therefore a potential candidate gene for engineering fast-growing and low-lignin trees with improved biomass production.

Structural insights into catalytic promiscuity of chalcone synthase from Glycine max (L.) Merr.: Coenzyme A-induced alteration of product specificity

Catalytic promiscuity of enzymes plays a pivotal role in driving the evolution of plant specialized metabolism. Chalcone synthase (CHS) catalyzes the production of 2',4,4',6'-tetrahydroxychalcone (THC), a common precursor of plant flavonoids, from p-coumaroyl-coenzyme A (-CoA) and three malonyl-CoA molecules. CHS has promiscuous product specificity, producing a significant amount of p-coumaroyltriacetic lactone (CTAL) in vitro. However, mechanistic aspects of this CHS promiscuity remain to be clarified. Here, we show that the product specificity of soybean CHS (GmCHS1) is altered by CoA, a reaction product, which selectively inhibits THC production (IC50, 67 μM) and enhances CTAL production. We determined the structure of a ternary GmCHS1/CoA/naringenin complex, in which CoA is bound to the CoA-binding tunnel via interactions with Lys55, Arg58, and Lys268. Replacement of these residues by alanine resulted in an enhanced THC/CTAL production ratio, suggesting the role of these residues in the CoA-mediated alteration of product specificity. In the ternary complex, a mobile loop ("the K-loop"), which contains Lys268, was in a "closed conformation" placing over the CoA-binding tunnel, whereas in the apo and binary complex structures, the K-loop was in an "open conformation" and remote from the tunnel. We propose that the production of THC involves a transition of the K-loop conformation between the open and closed states, whereas synthesis of CTAL is independent of it. In the presence of CoA, an enzyme conformer with the closed K-loop conformation becomes increasingly dominant, hampering the transition of K-loop conformations to result in decreased THC production and increased CTAL production.

Structural and Interactional Analysis of the Flavonoid Pathway Proteins: Chalcone Synthase, Chalcone Isomerase and Chalcone Isomerase-like Protein

Chalcone synthase (CHS) and chalcone isomerase (CHI) catalyze the first two committed steps of the flavonoid pathway that plays a pivotal role in the growth and reproduction of land plants, including UV protection, pigmentation, symbiotic nitrogen fixation, and pathogen resistance. Based on the obtained X-ray crystal structures of CHS, CHI, and chalcone isomerase-like protein (CHIL) from the same monocotyledon, Panicum virgatum, along with the results of the steady-state kinetics, spectroscopic/thermodynamic analyses, intermolecular interactions, and their effect on each catalytic step are proposed. In addition, PvCHI's unique activity for both naringenin chalcone and isoliquiritigenin was analyzed, and the observed hierarchical activity for those type-I and -II substrates was explained with the intrinsic characteristics of the enzyme and two substrates. The structure of PvCHS complexed with naringenin supports uncompetitive inhibition. PvCHS displays intrinsic catalytic promiscuity, evident from the formation of p-coumaroyltriacetic acid lactone (CTAL) in addition to naringenin chalcone. In the presence of PvCHIL, conversion of p-coumaroyl-CoA to naringenin through PvCHS and PvCHI displayed ~400-fold increased Vmax with reduced formation of CTAL by 70%. Supporting this model, molecular docking, ITC (Isothermal Titration Calorimetry), and FRET (Fluorescence Resonance Energy Transfer) indicated that both PvCHI and PvCHIL interact with PvCHS in a non-competitive manner, indicating the plausible allosteric effect of naringenin on CHS. Significantly, the presence of naringenin increased the affinity between PvCHS and PvCHIL, whereas naringenin chalcone decreased the affinity, indicating a plausible feedback mechanism to minimize spontaneous incorrect stereoisomers. These are the first findings from a three-body system from the same species, indicating the importance of the macromolecular assembly of CHS-CHI-CHIL in determining the amount and type of flavonoids produced in plant cells.

Assessment and Partial Characterization of Candidate Genes in Dihydrochalcone and Arbutin Biosynthesis in an Apple-Pear Hybrid by De Novo Transcriptome Assembly

Apples (Malus × domestica Borkh.) and pears (Pyrus communis L.) are valuable crops closely related within the Rosaceae family with reported nutraceutical properties derived from secondary metabolites including phloridzin and arbutin, which are distinctive phenolic metabolites characterizing apples and pears, respectively. Here, we generated a de novo transcriptome assembly of an intergeneric hybrid between apple and pear, accumulating intermediate levels of phloridzin and arbutin. Combining RNA-seq, in silico functional annotation prediction, targeted gene expression analysis, and expression-metabolite correlations, we identified candidate genes for functional characterization, resulting in the identification of active arbutin synthases in the hybrid and parental genotypes. Despite exhibiting an active arbutin synthase in vitro, the natural lack of arbutin in apples is reasoned by the absence of the substrate and broad substrate specificity. Altogether, our study serves as the basis for future assessment of potential physiological roles of identified genes by genome editing of hybrids and pears.

Identification of two 6'-deoxychalcone 4'-glucosyltransferase genes in dahlia (Dahlia variabilis)

MAIN CONCLUSION

Two glycosyltransferase genes belonging to UGT88 family were identified to have 6'-deoxychalcone 4'-glucosyltransferase activity in dahlia. 6'-Deoxychalcones (isoliquiritigenin and butein) are important pigments for yellow and orange to red flower color. 6'-Deoxychalcones are glucosylated at the 4'-position in vivo, but the genes encoding 6'-deoxychalcone 4'-glucosyltransferase have not yet been identified. In our previous study, it was indicated that snapdragon (Antirrhinum majus) chalcone 4'-O-glucosyltransferase (Am4'CGT) has isoliquiritigenin 4'-glucosylation activity. Therefore, to identify genes encoding 6'-deoxychalcone 4'-glucosyltransferase in dahlia (Dahlia variabilis), genes expressed in ray florets that shared high homology with Am4'CGT were explored. As a result, c34671_g1_i1 and c35662_g1_i1 were selected as candidate genes for 6'-deoxychalcone 4'-glucosyltransferases in dahlia. We conducted transient co-overexpression of three genes (c34671_g1_i1 or c35662_g1_i1, dahlia aldo-keto reductase1 (DvAKR1) or soybean (Glycine max) chalcone reductase5 (GmCHR5), and chili pepper (Capsicum annuum) MYB transcription factor (CaMYBA)) in Nicotiana benthamiana by agroinfiltration. Transient overexpression of c34671_g1_i1, DvAKR1, and CaMYBA resulted in increase in the accumulation of isoliquiritigenin 4'-glucosides, isoliquiritigenin 4'-O-glucoside, and isoliquiritigenin 4'-O-[6-O-(malonyl)-glucoside]. However, transient overexpression of c35662_g1_i1, DvAKR1, and CaMYBA did not increase accumulation of isoliquiritigenin 4'-glucosides. Using GmCHR5 instead of DvAKR1 showed similar results suggesting that c34671_g1_i1 has isoliquiritigenin 4'-glucosyltransferase activity. In addition, we conducted co-overexpression of four genes (c34671_g1_i1, c35662_g1_i1 or Am4'CGT, DvAKR1 or GmCHR5, CaMYBA, and chalcone 3-hydroxylase from dahlia). Accumulation of butein 4'-O-glucoside and butein 4'-O-[6-O-(malonyl)-glucoside] was detected for c35662_g1_i1, suggesting that c35662_g1_i1 has butein 4'-glucosyltransferase activity. Recombinant enzyme analysis also supported butein 4'-glucosyltransferases activity of c35662_g1_i1. Therefore, our results suggested that both c34671_g1_i1 and c35662_g1_i1 are 6'-deoxychalcone 4'-glucosyltransferases but with different substrate preference.

Probing the Deoxyflavonoid Biosynthesis: Naringenin Chalcone Is a Substrate for the Reductase Synthesizing Isoliquiritigenin

Isoliquiritigenin formation is a key reaction during deoxyflavonoid biosynthesis, which is catalyzed by two enzymes, chalcone synthase (CHS) and reductase (CHR). The substrates for CHS are established. However, the substrate for CHR is unknown. In this study, an in vitro reaction was performed to confirm whether naringenin chalcone can be a substrate. Naringenin chalcone was used as a substrate during the CHR reaction. Analyzing the product revealed that isoliquiritigenin was produced from naringenin chalcone, indicating that naringenin chalcone is a substrate. This study is the first to identify a substrate for CHR, reveals that deoxyflavonoid biosynthesis diverges from naringenin chalcone, endorses the term "chalcone reductase," and answers the long-standing questions about doubly-labeled acetic acid uptake pattern in deoxyflavonoid biosynthesis.

Pour some sugar on me: The diverse functions of phenylpropanoid glycosylation

The phenylpropanoid metabolism is the source of a vast array of specialized metabolites that play diverse functions in plant growth and development and contribute to all aspects of plant interactions with their surrounding environment. These compounds protect plants from damaging ultraviolet radiation and reactive oxygen species, provide mechanical support for the plants to stand upright, and mediate plant-plant and plant-microorganism communications. The enormous metabolic diversity of phenylpropanoids is further expanded by chemical modifications known as "decorative reactions", including hydroxylation, methylation, glycosylation, and acylation. Among these modifications, glycosylation is the major driving force of phenylpropanoid structural diversification, also contributing to the expansion of their properties. Phenylpropanoid glycosylation is catalyzed by regioselective uridine diphosphate (UDP)-dependent glycosyltransferases (UGTs), whereas glycosyl hydrolases known as β-glucosidases are the major players in deglycosylation. In this article, we review how the glycosylation process affects key physicochemical properties of phenylpropanoids, such as molecular stability and solubility, as well as metabolite compartmentalization/storage and biological activity/toxicity. We also summarize the recent knowledge on the functional implications of glycosylation of different classes of phenylpropanoid compounds. A balance of glycosylation/deglycosylation might represent an essential molecular mechanism to regulate phenylpropanoid homeostasis, allowing plants to dynamically respond to diverse environmental signals.

Glyceollins from soybean: Their pharmacological effects and biosynthetic pathways

Flavonoids are a highly abundant class of secondary metabolites present in plants. Isoflavonoids, in particular, are primarily synthesized in leguminous plants within the subfamily Papilionoideae. Numerous reports have established the favorable role of isoflavonoids in preventing a range of human diseases. Among the isoflavonoid components, glyceollins are synthesized specifically in soybean plants and have displayed promising effects in mitigating the occurrence and progression of breast and ovarian cancers as well as other diseases. Consequently, glyceollins have become a sought-after natural component for promoting women's health. In recent years, extensive research has focused on investigating the molecular mechanism underlying the preventative properties of glyceollins against various diseases. Substantial progress has also been made toward elucidating the biosynthetic pathway of glyceollins and exploring potential regulatory factors. Herein, we provide a review of the research conducted on glyceollins since their discovery five decades ago (1972-2023). We summarize their pharmacological effects, biosynthetic pathways, and advancements in chemical synthesis to enhance our understanding of the molecular mechanisms of their function and the genes involved in their biosynthetic pathway. Such knowledge may facilitate improved glyceollin synthesis and the creation of health products based on glyceollins.

Discovering Dynamic Plant Enzyme Complexes in Yeast for Kratom Alkaloid Pathway Identification

Discovering natural product biosynthetic pathways of medicinal plants is challenging and laborious. Capturing the coregulation patterns of pathway enzymes, particularly transcriptomic regulation, has proven an effective method to accelerate pathway identification. In this study, we developed a yeast-based screening method to capture the protein-protein interactions (PPI) between plant enzymes, which is another useful pattern to complement the prevalent approach. Combining this method with plant multiomics analysis, we discovered four enzyme complexes and their organized pathways from kratom, an alkaloid-producing plant. The four pathway branches involved six enzymes, including a strictosidine synthase, a strictosidine β-D-glucosidase (MsSGD), and four medium-chain dehydrogenase/reductases (MsMDRs). PPI screening selected six MsMDRs interacting with MsSGD from 20 candidates predicted by multiomics analysis. Four of the six MsMDRs were then characterized as functional, indicating the high selectivity of the PPI screening method. This study highlights the opportunity of leveraging post-translational regulation features to discover novel plant natural product biosynthetic pathways.

Integrative analysis of the metabolome and transcriptome reveals the molecular regulatory mechanism of isoflavonoid biosynthesis in Ormosia henryi Prain

Flavonoids are important components of many phytopharmaceuticals, however, most studies on flavonoids and isoflavonoids have been conducted on herbaceous plants of the family Leguminosae, such as soybean, and less attention has been paid to woody plants. To fill this gap, we characterized the metabolome and transcriptome of five plant organs of Ormosia henryi Prain (OHP), a woody Leguminosae plant with great pharmaceutical value. Our results indicate that OHP possesses a relatively high content of isoflavonoids as well as significant diversity, with greater diversity of isoflavonoids in the roots. Combined with transcriptome data, the pattern of isoflavonoid accumulation was found to be highly correlated with differential expression genes. Furthermore, the use of trait-WGCNA network analysis identified OhpCHSs as a probable hub enzyme that directs the downstream isoflavonoid synthesis pathway. Transcription factors, such as MYB26, MYB108, WRKY53, RAV1 and ZFP3, were found to be involved in the regulation of isoflavonoid biosynthesis in OHP. Our findings will be beneficial for the biosynthesis and utilization of woody isoflavonoids.

Protein Modelling Highlighted Key Catalytic Sites Involved in Position-Specific Glycosylation of Isoflavonoids

Uridine diphosphate glycosyltransferases (UGTs) are known for promiscuity towards sugar acceptors, a valuable characteristic for host plants but not desirable for heterologous biosynthesis. UGTs characterized for the O-glycosylation of isoflavonoids have shown a variable efficiency, substrate preference, and OH site specificity. Thus, 22 UGTs with reported isoflavonoid O-glycosylation activity were analyzed and ranked for OH site specificity and catalysis efficiency. Multiple-sequence alignment (MSA) showed a 33.2% pairwise identity and 4.5% identical sites among selected UGTs. MSA and phylogenetic analysis highlighted a comparatively higher amino acid substitution rate in the N-terminal domain that likely led to a higher specificity for isoflavonoids. Based on the docking score, OH site specificity, and physical and chemical features of active sites, selected UGTs were divided into three groups. A significantly high pairwise identity (67.4%) and identical sites (31.7%) were seen for group 1 UGTs. The structural and chemical composition of active sites highlighted key amino acids that likely define substrate preference, OH site specificity, and glycosylation efficiency towards selected (iso)flavonoids. In conclusion, physical and chemical parameters of active sites likely control the position-specific glycosylation of isoflavonoids. The present study will help the heterologous biosynthesis of glycosylated isoflavonoids and protein engineering efforts to improve the substrate and site specificity of UGTs.

Machine learning uncovers accumulation mechanism of flavonoid compounds in Polygonatum cyrtonema Hua

The compositions and yield of flavonoid compounds of Polygonatum cyrtonema Hua (PC) are important indices of the quality of medicinal materials. However, the flavonoids compositions and accumulation mechanism are still unclear in PC. Here, we identified 22 flavonoids using widely-targeted metabolome analysis in 15 genotypes of PC. Then weighted gene co-expression network analysis based on 45 transcriptome samples was performed to construct 12 co-expressed modules, in which blue module highly correlated with flavonoids was identified. Furthermore, 4 feature genes including PcCHS1, PcCHI, PcCHS2 and PcCHR5 were identified from 94 hub genes in blue module via machine learning methods support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF), and their functions on metabolic flux of flavonoids pathway were confirmed by tobacco transient expression system. Our findings identified representative flavonoids and key enzymes in PC that provided new insight for elite breeding rich in flavonoids, and thus will be beneficial for rapid development of great potential economic and medicinal value of PC.

A Polyphenol Oxidase Catalyzes Aurone Synthesis in Marchantia polymorpha

Aurones constitute one of the major classes of flavonoids, with a characteristic furanone structure that acts as the C-ring of flavonoids. Members of various enzyme families are involved in aurone biosynthesis in different higher plants, suggesting that during evolution plants acquired the ability to biosynthesize aurones independently and convergently. Bryophytes also produce aurones, but the biosynthetic pathways and enzymes involved have not been determined. The present study describes the identification and characterization of a polyphenol oxidase (PPO) that acts as an aureusidin synthase (MpAS1) in the model liverwort, Marchantia polymorpha. Crude enzyme assays using an M. polymorpha line overexpressing MpMYB14 with high accumulation of aureusidin showed that aureusidin was biosynthesized from naringenin chalcone and converted to riccionidin A. This activity was inhibited by N-phenylthiourea, an inhibitor specific to enzymes of the PPO family. Of the six PPOs highly induced in the line overexpressing MpMyb14, one, MpAS1, was found to biosynthesize aureusidin from naringenin chalcone when expressed in Saccharomyces cerevisiae. MpAS1 also recognized eriodictyol chalcone, isoliquiritigenin and butein, showing the highest activity for eriodictyol chalcone. Members of the PPO family in M. polymorpha evolved independently from PPOs in higher plants, indicating that aureusidin synthases evolved in parallel in land plants.

Regulatory ligand binding in plant chalcone isomerase-like (CHIL) proteins

Chalcone isomerase-like (CHIL) is a noncatalytic protein that enhances flavonoid content in green plants by serving as a metabolite binder and a rectifier of chalcone synthase (CHS). Rectification of CHS catalysis occurs through direct protein-protein interactions between CHIL and CHS, which alter CHS kinetics and product profiles, favoring naringenin chalcone production. These discoveries raise questions about how CHIL proteins interact structurally with metabolites and how CHIL-ligand interactions affect interactions with CHS. Using differential scanning fluorimetry (DSF) on a CHIL protein from Vitis Vinifera (VvCHIL), we report positive thermostability effects are induced by the binding of naringenin chalcone and negative thermostability effects are induced by the binding of naringenin. Naringenin chalcone further causes positive changes to CHIL-CHS binding, while naringenin causes negative changes to CHIL-CHS binding. These results suggest that CHILs may act as sensors for ligand-mediated pathway feedback by influencing CHS function. The protein X-ray crystal structure of VvCHIL compared with the protein X-ray crystal structure of a CHIL from Physcomitrella patens, reveals key amino acid differences at a ligand binding site of VvCHIL that can be substituted to nullify the destabilizing effect caused by naringenin. Together these results support a role for CHIL proteins as metabolite sensors that modulate the committed step of the flavonoid pathway.

Daidzein Hydroxylation by CYP81E63 Is Involved in the Biosynthesis of Miroestrol in Pueraria mirifica

White Kwao Krua (Pueraria candollei var. mirifica), a Thai medicinal plant, is a rich source of phytoestrogens, especially isoflavonoids and chromenes. These phytoestrogens are well known; however, their biosynthetic genes remain largely uncharacterized. Cytochrome P450 (P450) is a large protein family that plays a crucial role in the biosynthesis of various compounds in plants, including phytoestrogens. Thus, we focused on P450s involved in the isoflavone hydroxylation that potentially participates in the biosynthesis of miroestrol. Three candidate P450s were isolated from the transcriptome libraries by considering the phylogenetic and expression data of each tissue of P. mirifica. The candidate P450s were functionally characterized both in vitro and in planta. Accordingly, the yeast microsome harboring PmCYP81E63 regiospecifically exhibited either 2' or 3' daidzein hydroxylation and genistein hydroxylation. Based on in silico calculation, PmCYP81E63 had higher binding energy with daidzein than with genistein, which supported the in vitro result of the isoflavone specificity. To confirm in planta function, the candidate P450s were then transiently co-expressed with isoflavone-related genes in Nicotiana benthamiana. Despite no daidzein in the infiltrated N. benthamiana leaves, genistein and hydroxygenistein biosynthesis were detectable by liquid Chromatography with tandem mass spectrometry (LC-MS/MS). Additionally, we demonstrated that PmCYP81E63 interacted with several enzymes related to isoflavone biosynthesis using bimolecular fluorescence complementation studies and a yeast two-hybrid analysis, suggesting a scheme of metabolon formation in the pathway. Our findings provide compelling evidence regarding the involvement of PmCYP81E63 in the early step of the proposed miroestrol biosynthesis in P. mirifica.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Discovering dynamic plant enzyme complexes in yeast for novel alkaloid pathway identification from a medicinal plant kratom

Discovering natural product biosynthetic pathways from medicinal plants is challenging and laborious, largely due to the complexity of the transcriptomics-driven pathway prediction process. Here we developed a novel approach that captures the protein-level connections between enzymes for pathway discovery with improved accuracy. We proved that heterologous protein-protein interaction screening in yeast enabled the efficient discovery of both dynamic plant enzyme complexes and the pathways they organize. This approach discovered complexes and pathways in the monoterpene indole alkaloid metabolism of a medicinal plant, kratom with high success rate. Screening using a strictosidine β-D-glucosidase (MsSGD1) against 19 medium-chain dehydrogenase/reductases (MsMDRs) identified five MsSGD1-MsMDR complexes. Three out of the five interacting MsMDRs were then proven functional, while the remaining 14 non-interacting candidates did not show obvious activities. The work discovered three branched pathways by combining transcriptomics, metabolomics, and heterologous PPI screening and demonstrated a new plant pathway discovery strategy.

Interaction of PKR with STCS1: an indispensable step in the biosynthesis of lunularic acid in Marchantia polymorpha

Unlike bibenzyls derived from the vascular plants, lunularic acid (LA), a key precursor for macrocyclic bisbibenzyl synthesis in nonvascular liverworts, exhibits the absence of one hydroxy group within the A ring. It was hypothesized that both polyketide reductase (PKR) and stilbenecarboxylate synthase 1 (STCS1) were involved in the LA biosynthesis, but the underlined mechanisms have not been clarified. This study used bioinformatics analysis with molecular, biochemical and physiological approaches to characterize STCS1s and PKRs involved in the biosynthesis of LA. The results indicated that MpSTCS1s from Marchantia polymorpha catalyzed both C2→C7 aldol-type and C6→C1 Claisen-type cyclization using dihydro-p-coumaroyl-coenzyme A (CoA) and malonyl-CoA as substrates to yield a C6-C2-C6 skeleton of dihydro-resveratrol following decarboxylation and the C6-C3-C6 type of phloretin in vitro. The protein-protein interaction of PKRs with STCS1 (PPI-PS) was revealed and proved essential for LA accumulation when transiently co-expressed in Nicotiana benthamiana. Moreover, replacement of the active domain of STCS1 with an 18-amino-acid fragment from the chalcone synthase led to the PPI-PS greatly decreasing and diminishing the formation of LA. The replacement also increased the chalcone formation in STCS1s. Our results highlight a previously unrecognized PPI in planta that is indispensable for the formation of LA.

Phylogenetic Analysis and Protein Modelling of Isoflavonoid Synthase Highlights Key Catalytic Sites towards Realising New Bioengineering Endeavours

Isoflavonoid synthase (IFS) is a critical enzyme for the biosynthesis of over 2400 isoflavonoids. Isoflavonoids are an important class of plant secondary metabolites that have a range of pharmaceutical and nutraceutical properties. With growing interest in isoflavonoids from both research and industrial perspectives, efforts are being forwarded to enhance isoflavonoid production in-planta and ex-planta; therefore, in-silico analysis and characterisation of available IFS protein sequences are needed. The present study is the first-ever attempt toward phylogenetic analysis and protein modelling of available IFS protein sequences. Phylogenetic analysis has shown that IFS amino acid sequences have 86.4% pairwise identity and 26.5% identical sites, and the sequences were grouped into six different clades. The presence of a β-hairpin and extra loop at catalytic sites of Trifolium pratense, Beta vulgaris and Medicago truncatula, respectively, compared with Glycyrrhiza echinata are critical structural differences that may affect catalytic function. Protein docking highlighted the preference of selected IFS for liquiritigenin compared with naringenin and has listed T. pratense as the most efficient candidate for heterologous biosynthesis of isoflavonoids. The in-silico characterisation of IFS represented in this study is vital in realising the new bioengineering endeavours and will help in the characterisation and selection of IFS candidate enzymes for heterologous biosynthesis of isoflavonoids.

Heterologous biosynthesis of isobavachalcone in tobacco based on in planta screening of prenyltransferases

Isobavachalcone (IBC) is a prenylated chalcone mainly distributed in some Fabaceae and Moraceae species. IBC exhibits a wide range of pharmacological properties, including anti-bacterial, anti-viral, anti-inflammatory, and anti-cancer activities. In this study, we attempted to construct the heterologous biosynthesis pathway of IBC in tobacco (Nicotiana tabacum). Four previously reported prenyltransferases, including GuILDT from Glycyrrhiza uralensis, HlPT1 from Humulus lupulus, and SfILDT and SfFPT from Sophora flavescens, were subjected to an in planta screening to verify their activities for the biosynthesis of IBC, by using tobacco transient expression with exogenous isoliquiritigenin as the substrate. Only SfFPT and HlPT1 could convert isoliquiritigenin to IBC, and the activity of SfFPT was higher than that of HlPT1. By co-expression of GmCHS8 and GmCHR5 from Glycine max, endogenous isoliquiritigenin was generated in tobacco leaves (21.0 μg/g dry weight). After transformation with a multigene vector carrying GmCHS8, GmCHR5, and SfFPT, de novo biosynthesis of IBC was achieved in transgenic tobacco T₀ lines, in which the highest amount of IBC was 0.56 μg/g dry weight. The yield of IBC in transgenic plants was nearly equal to that in SfFPT transient expression experiments, in which substrate supplement was sufficient, indicating that low IBC yield was not attributed to the substrate supplement. Our research provided a prospect to produce valuable prenylflavonoids using plant-based metabolic engineering.

Evaluation of pathways to the C-glycosyl isoflavone puerarin in roots of kudzu (Pueraria montana lobata)

Kudzu (Pueraria montana lobata) is used as a traditional medicine in China and Southeast Asia but is a noxious weed in the Southeastern United States. It produces both O- and C-glycosylated isoflavones, with puerarin (C-glucosyl daidzein) as an important bioactive compound. Currently, the stage of the isoflavone pathway at which the C-glycosyl unit is added remains unclear, with a recent report of direct C-glycosylation of daidzein contradicting earlier labeling studies supporting C-glycosylation at the level of chalcone. We have employed comparative mRNA sequencing of the roots from two Pueraria species, one of which produces puerarin (field collected P. montana lobata) and one of which does not (commercial Pueraria phaseoloides), to identify candidate uridine diphosphate glycosyltransferase (UGT) enzymes involved in puerarin biosynthesis. Expression of recombinant UGTs in Escherichia coli and candidate C-glycosyltransferases in Medicago truncatula were used to explore substrate specificities, and gene silencing of UGT and key isoflavone biosynthetic genes in kudzu hairy roots employed to test hypotheses concerning the substrate(s) for C-glycosylation. Our results confirm UGT71T5 as a C-glycosyltransferase of isoflavone biosynthesis in kudzu. Enzymatic, isotope labeling, and genetic analyses suggest that puerarin arises both from the direct action of UGT71T5 on daidzein and via a second route in which the C-glycosidic linkage is introduced to the chalcone isoliquiritigenin.

A novel aldo-keto reductase gene is involved in 6'-deoxychalcone biosynthesis in dahlia (Dahlia variabilis)

A novel gene belonging to the aldo-keto reductase 13 family is involved in isoliquiritigenin biosynthesis in dahlia. The yellow pigments of dahlia flowers are derived from 6'-deoxychalcones, which are synthesized via a two-step process, involving the conversion of 3-malonyl-CoA and 4-coumaloyl-CoA into isoliquiritigenin in the first step, and the subsequent generation of butein from isoliquiritigenin. The first step reaction is catalyzed by chalcone synthase (CHS) and aldo-keto reductase (AKR). AKR has been implicated in the isoflavone biosynthesis in legumes, however, isolation of butein biosynthesis related AKR members are yet to be reported. A comparative RNA-seq analysis between two dahlia cultivars, 'Shukuhai' and its butein-deficient lateral mutant 'Rinka', was used in this study to identify a novel AKR gene involved in 6'-deoxychalcone biosynthesis. DvAKR1 encoded a AKR 13 sub-family protein with significant differential expression levels, and was phylogenetically distinct from the chalcone reductases, which belongs to the AKR 4A sub-family in legumes. DNA sequence variation and expression profiles of DvAKR1 gene were correlated with 6'-deoxychalcone accumulation in the tested dahlia cultivars. A single over-expression analysis of DvAKR1 was not sufficient to initiate the accumulation of isoliquiritigenin in tobacco, in contrast, its co-overexpression with a chalcone 4'-O-glucosyltransferase (Am4'CGT) from Antirrhinum majus and a MYB transcription factor, CaMYBA from Capsicum annuum successfully induced isoliquiritigenin accumulation. In addition, DvAKR1 homologous gene expression was detected in Coreopsideae species accumulating 6'-deoxychalcone, but not in Asteraceae species lacking 6'-deoxychalcone production. These results not only demonstrate the involvement of DvAKR1 in the biosynthesis of 6'-deoxychalcone in dahlia, but also show that 6'-deoxychalcone occurrence in Coreopsideae species developed evolutionarily independent from legume species.

1-Octen-3-ol is formed from its primeveroside after mechanical wounding of soybean leaves

KEY MESSAGE

Hydrolysis of 1-octen-3-yl β-primeveroside implemented by a system with high structure-specificity is accountable for the rapid formation of 1-octen-3-ol from soybean leaves after mechanical wounding. 1-Octen-3-ol is a volatile compound ubiquitous in fungi; however, a subset of plant species also has the ability to form 1-octen-3-ol. Owing to its volatile nature, it has been anticipated that 1-octen-3-ol is associated with the effort of the emitter to control the behavior of the surrounding organisms; however, its ecological significance and the enzymes involved in its biosynthesis have not been fully elucidated, particularly in plants. We previously found that soybean (Glycine max) seeds contain 1-octen-3-yl β-primeveroside (pri). To elucidate the physiological significance and the biosynthesis of 1-octen-3-ol in plants, changes in the amount of 1-octen-3-yl pri during development of soybean plants was examined. A high 1-octen-3-yl pri level was found in young developing green organs, such as young leaves and sepals. Treatment of soybean leaves with methyl jasmonates resulted in a significant increase in the amount of 1-octen-3-yl pri; suggesting its involvement in defense responses. Although 1-octen-3-ol was below the detection limit in intact soybean leaves, mechanical damage to the leaves caused rapid hydrolysis of almost all 1-octen-3-yl pri to liberate volatile 1-octen-3-ol. Under the same conditions, the other glycosides, including isoflavone glycoside and linalool diglycoside, were hardly hydrolyzed. Therefore, the enzyme system to liberate aglycone from glycosides in soybean leaves should have strict substrate specificity. 1-Octen-3-yl pri might function as a storage form of volatile 1-octen-3-ol for immediate response against stresses accompanying tissue wounding.

Silencing of Dihydroflavonol 4-reductase in Chrysanthemum Ray Florets Enhances Flavonoid Biosynthesis and Antioxidant Capacity

Flavonoid biosynthesis requires the activities of several enzymes, which form weakly-bound, ordered protein complexes termed metabolons. To decipher flux regulation in the flavonoid biosynthetic pathway of chrysanthemum (Chrysanthemum morifolium Ramat), we suppressed the gene-encoding dihydroflavonol 4-reductase (DFR) through RNA interference (RNAi)-mediated post-transcriptional gene silencing under a floral-specific promoter. Transgenic CmDFR-RNAi chrysanthemum plants were obtained by Agrobacterium-mediated transformation. Genomic PCR analysis of CmDFR-RNAi chrysanthemums propagated by several rounds of stem cuttings verified stable transgene integration into the genome. CmDFR mRNA levels were reduced by 60–80% in CmDFR-RNAi lines compared to those in wild-type (WT) plants in ray florets, but not leaves. Additionally, transcript levels of flavonoid biosynthetic genes were highly upregulated in ray florets of CmDFR-RNAi chrysanthemum relative to those in WT plants, while transcript levels in leaves were similar to WT. Total flavonoid contents were high in ray florets of CmDFR-RNAi chrysanthemums, but flavonoid contents of leaves were similar to WT, consistent with transcript levels of flavonoid biosynthetic genes. Ray florets of CmDFR-RNAi chrysanthemums exhibited stronger antioxidant capacity than those of WT plants. We propose that post-transcriptional silencing of CmDFR in ray florets modifies metabolic flux, resulting in enhanced flavonoid content and antioxidant activity.

MULTI-SCULPT: Multiplex Integration via Selective, CRISPR-Mediated, Ultralong Pathway Transformation in Yeast for Plant Natural Product Synthesis

Yeast has been a versatile model host for complex and valuable natural product biosynthesis via the reconstruction of heterologous biosynthetic pathways. Recent advances in natural product pathway elucidation have uncovered many large and complicated plant pathways that contain 10-30 genes for the biosynthesis of structurally complex, valuable natural products. However, the ability to reconstruct ultralong pathways efficiently in yeast does not match the increasing demand for valuable plant natural product biomanufacturing. Here, we developed a one-pot, multigene pathway integration method in yeast, named MULTI-SCULPT for multiplex integration via selective, CRISPR-mediated, ultralong pathway transformation. Leveraging multilocus genomic disruption via CRISPR/Cas9, newly developed native and synthetic genetic parts, and fine-tuned gene integration and characterization methods, we managed to integrate 21 DNA inserts that contain a 12-gene plant isoflavone biosynthetic pathway into yeast with a 90-100% success rate in 12 days. This method enables fast and efficient ultralong biosynthetic pathway integration and can allow for the fast iterative integration of even longer pathways in the future. Ultimately, this method will accelerate combinatorial optimization of elucidated plant natural product pathways and accelerate putative pathway characterization heterologously.

Systems Metabolic Engineering of Escherichia coli Coculture for De Novo Production of Genistein

Genistein is a plant-derived isoflavone possessing various bioactivities to prevent aging, carcinogenesis, and neurodegenerative and inflammation diseases. As a typical complex flavonoid, its microbial production from sugar remains to be completed. Here, we use systems metabolic engineering stategies to design and develop a three-strain commensalistic Escherichia coli coculture that for the first time realized the de novo production of genistein. First, we reconstituted the naringenin module by screening and incorporating chalcone isomerase-like protein, an auxiliary component to rectify the chalcone synthase promiscuity. Furthermore, we devised and constructed the genistein module by N-terminal modifications of plant P450 enzyme 2-hydroxyisoflavanone synthase and cytochrome P450 enzyme reductase. When naringenin-producing strain was cocultivated with p-coumaric acid-overproducing strain (a phenylalanine-auxotroph), two-strain coculture worked as commensalism through a unidirectional nutrient flow, which favored the efficient production of naringenin with a titer of 206.5 mg/L from glucose. A three-strain commensalistic coculture was subsequently engineered, which produced the highest titer to date of 60.8 mg/L genistein from a glucose and glycerol mixture. The commensalistic coculture is a flexible and versatile platform for the production of flavonoids, indicating a promising future for production of complex natural products in engineered E. coli.

Biochemistry and regulation of aurone biosynthesis

Aurones are a group of flavonoids that confer a bright yellow coloration to certain ornamental flowers and are a promising structural target for the development of new therapeutic drugs. Since the first identification of the snapdragon aurone synthase as a polyphenol oxidase (PPO) in 2000, several important advances in the biochemistry and regulation of aurone biosynthesis have been achieved. For example, several other aurone synthases have been identified in distantly related plants, which not only include PPOs but also peroxidases. Elucidation of the subcellular localization of aurone biosynthesis in snapdragon led to the establishment of a method to genetically engineer novel yellow flowers. The crystal structure of an aurone-producing PPO was clarified and provided important insights into the structure-function relationship of aurone-producing PPOs. A locus (SULFUREA) that negatively regulates aurone biosynthesis in snapdragon was identified, illustrating the evolution of flower color pattern through selection on regulatory small RNAs.

Phylogenomic classification and synteny network analyses deciphered the evolutionary landscape of aldo-keto reductase (AKR) gene superfamily in the plant kingdom

Aldo-keto reductase-domain (PF00248) containing proteins (AKRs) are NAD(P)(H)-dependent oxidoreductases of a multigene superfamily that mediate versatile functions in plants ranging from detoxification, metal chelation, potassium ion efflux to specialized metabolism. To uncover the complete repertoire of AKR gene superfamily in plants, a systematic kingdom-wide identification, phylogeny reconstruction, classification and synteny network clustering analyses were performed in this study using 74 diverse plant genomes. Plant AKRs were omnipresent, legitimately classified into 4 groups (based on phylogeny) and 14 subgroups (based on the ≥ 60% of protein sequence identity). Species composition of AKR subgroups highlights their distinct emergence during plant evolution. Loss of AKR subgroups among plants was apparent and that various lineage-, order/family- and species-specific losses were observed. The subgroups IA, IVB and IVF were flourished and diversified well during plant evolution, likely related to the complexity of plant's specialized metabolism and environmental adaptation. About 65% of AKRs were in genomic synteny regions across the plant kingdom and the AKRs relevant to important functions (e.g. vitamin B₆ metabolism) were in profoundly conserved angiosperm-wide synteny communities. This study underscores the evolutionary landscape of plant AKRs and provides a comprehensive resource to facilitate the functional characterization of them.

CRISPR/Cas9-Mediated Targeted Mutagenesis of GmUGT Enhanced Soybean Resistance Against Leaf-Chewing Insects Through Flavonoids Biosynthesis

Leaf-chewing insects are important pests that cause yield loss and reduce seed quality in soybeans (Glycine max). Breeding soybean varieties that are resistant to leaf-chewing insects can minimize the need for insecticide use and reduce yield loss. The marker gene for QTL-M, Glyma.07g110300 (LOC100775351) that encodes a UDP-glycosyltransferase (UGT) is the major determinant of resistance against leaf-chewing insects in soybean; it exhibits a loss of function in insect-resistant soybean germplasms. In this study, Agrobacterium-mediated transformation introduced the CRISPR/Cas9 expression vector into the soybean cultivar Tianlong No. 1 to generate Glyma.07g110300-gene mutants. We obtained two novel types of mutations, a 33-bp deletion and a single-bp insertion in the GmUGT coding region, which resulted in an enhanced resistance to Helicoverpa armigera and Spodoptera litura. Additionally, overexpressing GmUGT produced soybean varieties that were more sensitive to H. armigera and S. litura. Both mutant and overexpressing lines exhibited no obvious phenotypic changes. The difference in metabolites and gene expression suggested that GmUGT is involved in imparting resistance to leaf-chewing insects by altering the flavonoid content and expression patterns of genes related to flavonoid biosynthesis and defense. Furthermore, ectopic expression of the GmUGT gene in the ugt72b1 mutant of Arabidopsis substantially rescued the phenotype of H. armigera resistance in the atugt72b1 mutant. Our study presents a strategy for increasing resistance against leaf-chewing insects in soybean through CRISPR/Cas9-mediated targeted mutagenesis of the UGT genes.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Identification and characterization of unique 5-hydroxyisoflavonoid biosynthetic key enzyme genes in Lupinus albus

5-Hydroxyisoflavonoids, no 5-deoxyisoflavonoids, in Lupinus species, are due to lack of CHRs and Type II CHIs, and the key enzymes of isoflavonoid biosynthetic pathway in white lupin were identified. White lupin (Lupinus albus) is used as food ingredients owing to rich protein, low starch, and rich bioactive compounds such as isoflavonoids. The isoflavonoids biosynthetic pathway in white lupin still remains unclear. In this study, only 5-hydroxyisoflavonoids, but no 5-deoxyisoflavonoids, were detected in white lupin and other Lupinus species. No 5-deoxyisoflavonoids in Lupinus species are due to lack of CHRs and Type II CHIs. We further found that the CHI gene cluster containing both Type I and Type II CHIs possibly arose after the divergence of Lupinus with other legume clade. LaCHI1 and LaCHI2 identified from white lupin metabolized naringenin chalcone to naringenin in yeast and tobacco (Nicotiana benthamiana), and were bona fide Type I CHIs. We further identified two isoflavone synthases (LaIFS1 and LaIFS2), catalyzing flavanone naringenin into isoflavone genistein and also catalyzing liquiritigenin into daidzein in yeast and tobacco. In addition, LaG6DT1 and LaG6DT2 prenylated genistein at the C-6 position into wighteone. Two glucosyltransferases LaUGT1 and LaUGT2 metabolized genistein and wighteone into its 7-O-glucosides. Taken together, our study not only revealed that exclusive 5-hydroxyisoflavonoids do exist in Lupinus species, but also identified key enzymes in the isoflavonoid biosynthetic pathway in white lupin.

Functional Characterization of a Flavone Synthase That Participates in a Kumquat Flavone Metabolon

Flavones predominantly accumulate as O- and C-glycosides in kumquat plants. Two catalytic mechanisms of flavone synthase II (FNSII) support the biosynthesis of glycosyl flavones, one involving flavanone 2-hydroxylase (which generates 2-hydroxyflavanones for C-glycosylation) and another involving the direct catalysis of flavanones to flavones for O-glycosylation. However, FNSII has not yet been characterized in kumquats. In this study, we identified two kumquat FNSII genes (FcFNSII-1 and FcFNSII-2), based on transcriptome and bioinformatics analysis. Data from in vivo and in vitro assays showed that FcFNSII-2 directly synthesized apigenin and acacetin from naringenin and isosakuranetin, respectively, whereas FcFNSII-1 showed no detectable catalytic activities with flavanones. In agreement, transient overexpression of FcFNSII-2 in kumquat peels significantly enhanced the transcription of structural genes of the flavonoid-biosynthesis pathway and the accumulation of several O-glycosyl flavones. Moreover, studying the subcellular localizations of FcFNSII-1 and FcFNSII-2 demonstrated that N-terminal membrane-spanning domains were necessary to ensure endoplasmic reticulum localization and anchoring. Protein-protein interaction analyses, using the split-ubiquitin yeast two-hybrid system and bimolecular fluorescence-complementation assays, revealed that FcFNSII-2 interacted with chalcone synthase 1, chalcone synthase 2, and chalcone isomerase-like proteins. The results provide strong evidence that FcFNSII-2 serves as a nucleation site for an O-glycosyl flavone metabolon that channels flavanones for O-glycosyl flavone biosynthesis in kumquat fruits. They have implications for guiding genetic engineering efforts aimed at enhancing the composition of bioactive flavonoids in kumquat fruits.

Identification of the Genes Coding for Carthamin Synthase, Peroxidase Homologs that Catalyze the Final Enzymatic Step of Red Pigmentation in Safflower (Carthamus tinctorius L.)

Carthamin, a dimeric quinochalcone that is sparingly soluble in water, is obtained from the yellow-orange corolla of fully blooming safflower (Carthamus tinctorius L.) florets. Carthamin is a natural red colorant, which has been used worldwide for more than 4500 years and is the major component of Japanese 'beni' used for dyeing textiles, in cosmetics and as a food colorant. The biosynthetic pathway of carthamin has long remained uncertain. Previously, carthamin was proposed to be derived from precarthamin (PC), a water-soluble quinochalcone, via a single enzymatic process. In this study, we identified the genes coding for the enzyme responsible for the formation of carthamin from PC, termed 'carthamin synthase' (CarS), using enzyme purification and transcriptome analysis. The CarS proteins were purified from the cream-colored corolla of safflower and identified as peroxidase homologs (CtPOD1, CtPOD2 and CtPOD3). The purified enzyme catalyzed the oxidative decarboxylation of PC to produce carthamin using O2, instead of H2O2, as an electron acceptor. In addition, CarS catalyzed the decomposition of carthamin. However, this enzymatic decomposition of carthamin could be circumvented by adsorption of the pigment to cellulose. These CtPOD isozymes were not only expressed in the corolla of the carthamin-producing orange safflower cultivars but were also abundantly expressed in tissues and organs that did not produce carthamin and PC. One CtPOD isozyme, CtPOD2, was localized in the extracellular space. Based on the results obtained, a model for the stable red pigmentation of safflower florets during flower senescence and the traditional 'beni' manufacturing process is proposed.

Diversion of metabolic flux towards 5-deoxy(iso)flavonoid production via enzyme self-assembly in Escherichia coli

5-Deoxy(iso)flavonoids are structural representatives of phenylpropanoid-derived compounds and play critical roles in plant ecophysiology. Recently, 5-deoxy(iso)flavonoids gained significant interest due to their potential applications as pharmaceuticals, nutraceuticals, and food additives. Given the difficulties in their isolation from native plant sources, engineered biosynthesis of 5-deoxy(iso)flavonoids in a microbial host is a highly promising alternative approach. However, the production of 5-deoxy(iso)flavonoids is hindered by metabolic flux imbalances that result in a product profile predominated by non-reduced analogues. In this study, GmCHS7 (chalcone synthase from Glycine max) and GuCHR (chalcone reductase from Glycyrrhizza uralensis) were preliminarily utilized to improve the CHR ratio (CHR product to total CHS product). The use of this enzyme combination improved the final CHR ratio from 39.7% to 50.3%. For further optimization, a protein-protein interaction strategy was employed, basing on the spatial adhesion of GmCHS7:PDZ and GuCHR:PDZlig. This strategy further increased the ratio towards the CHR-derived product (54.7%), suggesting partial success of redirecting metabolic flux towards the reduced branch. To further increase the total carbon metabolic flux, 15 protein scaffolds were programmed with stoichiometric arrangement of the three sequential catalysts GmCHS7, GuCHR and MsCHI (chalcone isomerase from Medicago sativa), resulting in a 1.4-fold increase in total flavanone production, from 69.4 mg/L to 97.0 mg/L in shake flasks. The protein self-assembly strategy also improved the production and direction of the lineage-specific compounds 7,4'-dihydroxyflavone and daidzein in Escherichia coli. This study presents a significant advancement of 5-deoxy(iso)flavonoid production and provides the foundation for production of value-added 5-deoxy(iso)flavonoids in microbial hosts.

Nutritional evaluation and transcriptome analyses of short-time germinated seeds in soybean (Glycine max L. Merri.)

Germination is a common practice for nutrition improvement in many crops. In soybean, the nutrient value and genome-wide gene expression pattern of whole seeds germinated for short-time has not been fully investigated. In this study, protein content (PC), water soluble protein content (WSPC), isoflavone compositions were evaluated at 0 and 36 h after germination (HAG), respectively. The results showed that at 36HAG, PC was slightly decreased (P > 0.05) in ZD41, J58 and JHD, WSPC and free isoflavone (aglycones: daidzein, genistein, and glycitein) were significantly increased (P < 0.05), while total isoflavone content was unchanged. Transcriptomic analysis identified 5240, 6840 and 15,766 DEGs in different time point comparisons, respectively. GO and KEGG analysis showed that photosynthesis process was significantly activated from 18HAG, and alternative splicing might play an important role during germination in a complex manner. Response to hydrogen peroxide (H2O2) was found to be down regulated significantly from 18 to 36HAG, suggesting that H2O2 might play an important role in germination. Expression pattern analysis showed the synthesis of storage proteins was slowing down, while the genes coding for protein degradation (peptidase and protease) were up regulated as time went by during germination. For genes involved in isoflavone metabolism pathway, UGT (7-O-glucosyltransferase) coding genes were significantly up regulated (40 up-DEGs vs 27 down-DEGs), while MAT (7-O-glucoside-6''-O-malonyltransferase) coding genes were down regulated, which might explain the increase of aglycones after germination. This study provided a universal transcriptomic atlas for whole soybean seeds germination in terms of nutrition and gene regulation mechanism.

Recent advances on the roles of flavonoids as plant protective molecules after UV and high light exposure

Flavonoids are plant specialized metabolites that consist of one oxygenated and two aromatic rings. Different flavonoids are grouped according to the oxidation degree of the carbon rings; they can later be modified by glycosylations, hydroxylations, acylations, methylations, or prenylations. These modifications generate a wide collection of different molecules which have various functions in plants. All flavonoids absorb in the UV wavelengths, they mostly accumulate in the epidermis of plant cells and their biosynthesis is generally activated after UV exposure. Therefore, they have been assumed to protect plants against exposure to radiation in this range. Some flavonoids also absorb in other wavelengths, for example anthocyanins, which absorb light in the visible part of the solar spectrum. Besides, some flavonoids show antioxidant properties, that is, they act as scavengers of reactive oxygen species that could be produced after high fluence UV exposure. However, to date most reports were based on in vitro studies, and there is very little in vivo evidence of how their roles are carried out. In this review we first summarize the biosynthetic pathway of flavonoids and their characteristics, and we describe recent advances on the investigation of the role of three of the most abundant flavonoids: flavonols, flavones, and anthocyanins, protecting plants against UV exposure and high light exposure. We also present examples of how using UV-B supplementation to increase flavonoid content, is possible to improve plant nutritional and pharmaceutical values.

Enzyme Complexes of Ptr4CL and PtrHCT Modulate Co-enzyme A Ligation of Hydroxycinnamic Acids for Monolignol Biosynthesis in Populus trichocarpa

Co-enzyme A (CoA) ligation of hydroxycinnamic acids by 4-coumaric acid:CoA ligase (4CL) is a critical step in the biosynthesis of monolignols. Perturbation of 4CL activity significantly impacts the lignin content of diverse plant species. In Populus trichocarpa, two well-studied xylem-specific Ptr4CLs (Ptr4CL3 and Ptr4CL5) catalyze the CoA ligation of 4-coumaric acid to 4-coumaroyl-CoA and caffeic acid to caffeoyl-CoA. Subsequently, two 4-hydroxycinnamoyl-CoA:shikimic acid hydroxycinnamoyl transferases (PtrHCT1 and PtrHCT6) mediate the conversion of 4-coumaroyl-CoA to caffeoyl-CoA. Here, we show that the CoA ligation of 4-coumaric and caffeic acids is modulated by Ptr4CL/PtrHCT protein complexes. Downregulation of PtrHCTs reduced Ptr4CL activities in the stem-differentiating xylem (SDX) of transgenic P. trichocarpa. The Ptr4CL/PtrHCT interactions were then validated in vivo using biomolecular fluorescence complementation (BiFC) and protein pull-down assays in P. trichocarpa SDX extracts. Enzyme activity assays using recombinant proteins of Ptr4CL and PtrHCT showed elevated CoA ligation activity for Ptr4CL when supplemented with PtrHCT. Numerical analyses based on an evolutionary computation of the CoA ligation activity estimated the stoichiometry of the protein complex to consist of one Ptr4CL and two PtrHCTs, which was experimentally confirmed by chemical cross-linking using SDX plant protein extracts and recombinant proteins. Based on these results, we propose that Ptr4CL/PtrHCT complexes modulate the metabolic flux of CoA ligation for monolignol biosynthesis during wood formation in P. trichocarpa.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

Metabolons, enzyme-enzyme assemblies that mediate substrate channeling, and their roles in plant metabolism

Metabolons are transient multi-protein complexes of sequential enzymes that mediate substrate channeling. They differ from multi-enzyme complexes in that they are dynamic, rather than permanent, and as such have considerably lower dissociation constants. Despite the fact that a huge number of metabolons have been suggested to exist in plants, most of these claims are erroneous as only a handful of these have been proven to channel metabolites. We believe that physical protein-protein interactions between consecutive enzymes of a pathway should rather be called enzyme-enzyme assemblies. In this review, we describe how metabolons are generally assembled by transient interactions and held together by both structural elements and non-covalent interactions. Experimental evidence for their existence comes from protein-protein interaction studies, which indicate that the enzymes physically interact, and direct substrate channeling measurements, which indicate that they functionally interact. Unfortunately, advances in cell biology and proteomics have far outstripped those in classical enzymology and flux measurements, rendering most reports reliant purely on interactome studies. Recent developments in co-fractionation mass spectrometry will likely further exacerbate this bias. Given this, only dynamic enzyme-enzyme assemblies in which both physical and functional interactions have been demonstrated should be termed metabolons. We discuss the level of evidence for the manifold plant pathways that have been postulated to contain metabolons and then list examples in both primary and secondary metabolism for which strong evidence has been provided to support these claims. In doing so, we pay particular attention to experimental and mathematical approaches to study metabolons as well as complexities that arise in attempting to follow them. Finally, we discuss perspectives for improving our understanding of these fascinating but enigmatic interactions.

Metabolic Engineering of Isoflavones: An Updated Overview

Isoflavones are ecophysiologically active secondary metabolites derived from the phenylpropanoid pathway. They were mostly found in leguminous plants, especially in the pea family. Isoflavones play a key role in plant-environment interactions and act as phytoalexins also having an array of health benefits to the humans. According to epidemiological studies, a high intake of isoflavones-rich diets linked to a lower risk of hormone-related cancers, osteoporosis, menopausal symptoms, and cardiovascular diseases. These characteristics lead to the significant advancement in the studies on genetic and metabolic engineering of isoflavones in plants. As a result, a number of structural and regulatory genes involved in isoflavone biosynthesis in plants have been identified and characterized. Subsequently, they were engineered in various crop plants for the increased production of isoflavones. Furthermore, with the advent of high-throughput technologies, the regulation of isoflavone biosynthesis gains attention to increase or decrease the level of isoflavones in the crop plants. In the review, we begin with the role of isoflavones in plants, environment, and its benefits in human health. Besides, the main theme is to discuss the updated research progress in metabolic engineering of isoflavones in other plants species and regulation of production of isoflavones in soybeans.

Identification of a Unique Type of Isoflavone O-Methyltransferase, GmIOMT1, Based on Multi-Omics Analysis of Soybean under Biotic Stress

Isoflavonoids are commonly found in leguminous plants. Glycitein is one of the isoflavones produced by soybean. The genes encoding the enzymes in the isoflavone biosynthetic pathway have mostly been identified and characterized. However, the gene(s) for isoflavone O-methyltransferase (IOMT), which catalyzes the last step of glycitein biosynthesis, has not yet been identified. In this study, we conducted multi-omics analyses of fungal-inoculated soybean and indicated that glycitein biosynthesis was induced in response to biotic stress. Moreover, we identified a unique type of IOMT, which participates in glycitein biosynthesis. Soybean seedlings were inoculated with Aspergillus oryzae or Rhizopus oligosporus and sampled daily for 8 d. Multi-omics analyses were conducted using liquid chromatography-tandem mass spectrometry and RNA sequencing. Metabolome analysis revealed that glycitein derivatives increased following fungal inoculation. Transcriptome co-expression analysis identified two candidate IOMTs that were co-expressed with the gene encoding flavonoid 6-hydroxylase (F6H), the key enzyme in glycitein biosynthesis. The enzymatic assay of the two IOMTs using respective recombinant proteins showed that one IOMT, named as GmIOMT1, produced glycitein. Unlike other IOMTs, GmIOMT1 belongs to the cation-dependent OMT family and exhibited the highest activity with Zn2+ among cations tested. Moreover, we demonstrated that GmIOMT1 overexpression increased the levels of glycitein derivatives in soybean hairy roots when F6H was co-expressed. These results strongly suggest that GmIOMT1 participates in inducing glycitein biosynthesis in response to biotic stress.

Managing enzyme promiscuity in plant specialized metabolism: A lesson from flavonoid biosynthesis: Mission of a "body double" protein clarified

Specificities of enzymes involved in plant specialized metabolism, including flavonoid biosynthesis, are generally promiscuous. This enzyme promiscuity has served as an evolutionary basis for new enzyme functions and metabolic pathways in land plants adapting to environmental challenges. This phenomenon may lead, however, to inefficiency in specialized metabolism and adversely affect metabolite-mediated plant survival. How plants manage enzyme promiscuity for efficient specialized metabolism is, thus, an open question. Recent studies of flavonoid biosynthesis addressing this issue have revealed a conserved strategy, namely, a homolog of chalcone isomerase with no catalytic activity binds to chalcone synthase, a key flavonoid pathway enzyme, to narrow (or rectify) the enzyme's highly promiscuous product specificity. Reducing promiscuity via specific protein-protein interactions among metabolic enzymes and proteins may be a solution adopted by land plants to achieve efficient operation of specialized metabolism, while the intrinsic promiscuity of enzymes has likely been retained incidentally.

Saccharomyces Cerevisiae-An Interesting Producer of Bioactive Plant Polyphenolic Metabolites

Secondary phenolic metabolites are defined as valuable natural products synthesized by different organisms that are not essential for growth and development. These compounds play an essential role in plant defense mechanisms and an important role in the pharmaceutical, cosmetics, food, and agricultural industries. Despite the vast chemical diversity of natural compounds, their content in plants is very low, and, as a consequence, this eliminates the possibility of the production of these interesting secondary metabolites from plants. Therefore, microorganisms are widely used as cell factories by industrial biotechnology, in the production of different non-native compounds. Among microorganisms commonly used in biotechnological applications, yeast are a prominent host for the diverse secondary metabolite biosynthetic pathways. Saccharomyces cerevisiae is often regarded as a better host organism for the heterologous production of phenolic compounds, particularly if the expression of different plant genes is necessary.

Crystal structure of chalcone synthase, a key enzyme for isoflavonoid biosynthesis in soybean

Isoflavonoid is one of the groups of flavonoids that play pivotal roles in the survival of land plants. Chalcone synthase (CHS), the first enzyme of the isoflavonoid biosynthetic pathway, catalyzes the formation of a common isoflavonoid precursor. We have previously reported that an isozyme of soybean CHS (termed GmCHS1) is a key component of the isoflavonoid metabolon, a protein complex to enhance efficiency of isoflavonoid production. Here, we determined the crystal structure of GmCHS1 as a first step of understanding the metabolon structure, as well as to better understand the catalytic mechanism of GmCHS1.

Single-base deletion in GmCHR5 increases the genistein-to-daidzein ratio in soybean seed

Novel mutant alleles related to isoflavone content are useful for breeding programs to improve the disease resistance and nutritional content of soybean. However, identification of mutant alleles from high-density mutant libraries is expensive and time-consuming because soybean has a large, complicated genome. Here, we identified the gene responsible for increased genistein-to-daidzein ratio in seed of the mutant line F333ES017D9. For this purpose, we used a time- and cost-effective approach based on selective genotyping of a small number of F₂ plants showing the mutant phenotype with nearest-neighboring-nucleotide substitution-high-resolution melting analysis markers, followed by alignment of short reads obtained by next-generation sequencing analysis with the identified locus. In the mutant line, GmCHR5 harbored a single-base deletion that caused a change in the substrate flow in the isoflavone biosynthetic pathway towards genistein. Mutated GmCHR5 was expressed at a lower level during seed development than wild-type GmCHR5. Ectopic overexpression of GmCHR5 increased the production of daidzein derivatives in both the wild-type and mutant plants. The present strategy will be useful for accelerating identification of mutant alleles responsible for traits of interest in agronomically important crops.

De novo biosynthesis of liquiritin in Saccharomyces cerevisiae

Liquiritigenin (LG), isoliquiritigenin (Iso-LG), together with their respective glycoside derivatives liquiritin (LN) and isoliquiritin (Iso-LN), are the main active flavonoids of Glycyrrhiza uralensis, which is arguably the most widely used medicinal plant with enormous demand on the market, including Chinese medicine prescriptions, preparations, health care products and even food. Pharmacological studies have shown that these ingredients have broad medicinal value, including anti-cancer and anti-inflammatory effects. Although the biosynthetic pathway of glycyrrhizin, a triterpenoid component from G. uralensis, has been fully analyzed, little attention has been paid to the biosynthesis of the flavonoids of this plant. To obtain the enzyme-coding genes responsible for the biosynthesis of LN, analysis and screening were carried out by combining genome and comparative transcriptome database searches of G. uralensis and homologous genes of known flavonoid biosynthesis pathways. The catalytic functions of candidate genes were determined by in vitro or in vivo characterization. This work characterized the complete biosynthetic pathway of LN and achieved the de novo biosynthesis of liquiritin in Saccharomyces cerevisiae using endogenous yeast metabolites as precursors and cofactors for the first time, which provides a possibility for the economical and sustainable production and application of G. uralensis flavonoids through synthetic biology.