Genome-Wide Identification of Chalcone Reductase Gene Family in Soybean: Insight into Root-Specific GmCHRs and Phytophthora sojae Resistance

Integrative analysis of metabolome and transcriptome reveals regulatory mechanisms of flavonoid biosynthesis in soybean under salt stress

Introduction

Salt stress is a major environmental factor that constrains soybean growth, development, and productivity. Flavonoids are key secondary metabolites that play a crucial role in enhancing plant resistance to both biotic and abiotic stress. However, a comprehensive understanding of the regulatory mechanisms underlying flavonoid biosynthesis under salt stress in soybean is lacking.

Methods

In this study, an integrative analysis of soybean metabolome and transcriptome was conducted using two soybean lines, FQ03 (salt-sensitive, SS) and FQ07 (salt-tolerant, ST).

Results

A total of 650 significantly changed metabolites were identified in SS and ST after salt stress treatment. Among them, 151 flavonoids were categorized into nine classes, with flavones and flavonols being the predominant flavonoid types in soybean. Heatmap analysis showed higher contents of most flavonoid metabolites in ST than in SS under salt stress, and the total flavonoid content in ST was significantly higher than that in SS. In addition, transcriptome analysis revealed a higher number of differentially expressed genes (DEGs) in ST than in SS under salt stress. KEGG enrichment analysis revealed that DEGs were mainly enriched in pathways related to phenylpropanoid biosynthesis, isoflavonoid biosynthesis, flavonoid biosynthesis, as well as flavone and flavonol biosynthesis. Notably, 55 DEGs that were mapped to the flavonoid biosynthetic pathway were identified, with most showing higher expression levels in ST than in SS. Weighted gene correlation network analysis identified eight structural genes and six transcription factor genes as key regulators of flavonoid biosynthesis within the blue module. Furthermore, qRT-PCR results confirmed the accuracy of the transcriptomic data and reliability of the identified candidate genes.

Discussion

This study provides insights into the regulatory mechanisms underlying salt stress responses in soybean and highlights hub genes as potential targets for developing salt-tolerant soybean varieties.

Soybean AROGENATE DEHYDRATASES (GmADTs): involvement in the cytosolic isoflavonoid metabolon or trans-organelle continuity?

Soybean (Glycine max) produces a class of phenylalanine (Phe) derived specialized metabolites, isoflavonoids. Isoflavonoids are unique to legumes and are involved in defense responses in planta, and they are also necessary for nodule formation with nitrogen-fixing bacteria. Since Phe is a precursor of isoflavonoids, it stands to reason that the synthesis of Phe is coordinated with isoflavonoid production. Two putative AROGENATE DEHYDRATASE (ADT) isoforms were previously co-purified with the soybean isoflavonoid metabolon anchor ISOFLAVONE SYNTHASE2 (GmIFS2), however the GmADT family had not been characterized. Here, we present the identification of the nine member GmADT family. We determined that the GmADTs share sequences required for enzymatic activity and allosteric regulation with other characterized plant ADTs. Furthermore, the GmADTs are differentially expressed, and multiple members have dual substrate specificity, also acting as PREPHENATE DEHYDRATASES. All GmADT isoforms were detected in the stromules of chloroplasts, and they all interact with GmIFS2 in the cytosol. In addition, GmADT12A interacts with multiple other isoflavonoid metabolon members. These data substantiate the involvement of GmADT isoforms in the isoflavonoid metabolon.

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.

Cloning of a CHS gene of Poncirus trifoliata and its expression in response to soil water deficit and arbuscular mycorrhizal fungi

Flavonoids are secondary metabolites widely found in plants with antioxidants, of which chalcone synthase (CHS) is a key enzyme required in flavonoid synthesis pathways. The objective of this study was to clone a CHS gene from trifoliate orange (Poncirus trifoliata) and analyze its biological information and partial functions. A PtCHS gene (NCBI accession: MZ350874) was cloned from the genome-wide of trifoliate orange, which has 1156 bp in length, encoding 391 amino acids, with a predicted protein relative molecular mass of 42640.19, a theoretical isoelectric point of 6.28, and a lipid coefficient of 89.82. The protein is stable, hydrophilic, and high sequence conservation (92.49% sequence homology with CHS gene of other species). PtCHS was highly expressed in stems, leaves and flowers, but very low expression in roots and seeds. Soil water deficit could up-regulate expressions of PtCHS in leaves. An arbuscular mycorrhizal fungus, Funneliformis mosseae, significantly increased plant biomass production, CHS activity, expressions of PtCHS, and total flavonoid content in leaves and roots, independent of soil water status. Total flavonoids were significantly positively correlated with PtCHS expression in leaves only and also positively with root mycorrhizal colonization. Such results provide insight into the important functions of PtCHS in trifoliate orange.

Functions and biosynthesis of plant signaling metabolites mediating plant-microbe interactions

Covering: 2015-2022Plants and microbes have coevolved since their appearance, and their interactions, to some extent, define plant health. A reasonable fraction of small molecules plants produced are involved in mediating plant-microbe interactions, yet their functions and biosynthesis remain fragmented. The identification of these compounds and their biosynthetic genes will open up avenues for plant fitness improvement by manipulating metabolite-mediated plant-microbe interactions. Herein, we integrate the current knowledge on their chemical structures, bioactivities, and biosynthesis with the view of providing a high-level overview on their biosynthetic origins and evolutionary trajectory, and pinpointing the yet unknown and key enzymatic steps in diverse biosynthetic pathways. We further discuss the theoretical basis and prospects for directing plant signaling metabolite biosynthesis for microbe-aided plant health improvement in the future.

The SoyaGen Project: Putting Genomics to Work for Soybean Breeders

The SoyaGen project was a collaborative endeavor involving Canadian soybean researchers and breeders from academia and the private sector as well as international collaborators. Its aims were to develop genomics-derived solutions to real-world challenges faced by breeders. Based on the needs expressed by the stakeholders, the research efforts were focused on maximizing realized yield through optimization of maturity and improved disease resistance. The main deliverables related to molecular breeding in soybean will be reviewed here. These include: (1) SNP datasets capturing the genetic diversity within cultivated soybean (both within a worldwide collection of > 1,000 soybean accessions and a subset of 102 short-season accessions (MG0 and earlier) directly relevant to this group); (2) SNP markers for selecting favorable alleles at key maturity genes as well as loci associated with increased resistance to key pathogens and pests (Phytophthora sojae, Heterodera glycines, Sclerotinia sclerotiorum); (3) diagnostic tools to facilitate the identification and mapping of specific pathotypes of P. sojae; and (4) a genomic prediction approach to identify the most promising combinations of parents. As a result of this fruitful collaboration, breeders have gained new tools and approaches to implement molecular, genomics-informed breeding strategies. We believe these tools and approaches are broadly applicable to soybean breeding efforts around the world.

Changes in the Density and Composition of Rhizosphere Pathogenic Fusarium and Beneficial Trichoderma Contributing to Reduced Root Rot of Intercropped Soybean

The dynamic of soil-borne disease is closely related to the rhizosphere microbial communities. Maize–soybean relay strip intercropping has been shown to significantly control the type of soybean root rot that tends to occur in monoculture. However, it is still unknown whether the rhizosphere microbial community participates in the regulation of intercropped soybean root rot. In this study, rhizosphere Fusarium and Trichoderma communities were compared in either healthy or root-rotted rhizosphere soil from monocultured and intercropped soybean, and our results showed the abundance of rhizosphere Fusarium in intercropping was remarkably different from monoculture. Of four species identified, F. oxysporum was the most aggressive and more frequently isolated in diseased soil of monoculture. In contrast, Trichoderma was largely accumulated in healthy rhizosphere soil of intercropping rather than monoculture. T. harzianum dramatically increased in the rhizosphere of intercropping, while T. virens and T. afroharzianum also exhibited distinct isolation frequency. For the antagonism test in vitro, Trichoderma strains had antagonistic effects on F. oxysporum with the percentage of mycelial inhibition ranging from 50.59–92.94%, and they displayed good mycoparasitic abilities against F. oxysporum through coiling around and entering into the hyphae, expanding along the cell–cell lumen and even dissolving cell walls of the target fungus. These results indicate maize–soybean relay strip intercropping significantly increases the density and composition proportion of beneficial Trichoderma to antagonize the pathogenic Fusarium species in rhizosphere, thus potentially contributing to the suppression of soybean root rot under the intercropping.

Molecular Assessment of Pathotype Diversity of Phytophthora sojae in Canada Highlights Declining Sources of Resistance in Soybean

The large-scale deployment of resistance to Phytophthora sojae (Rps) genes in soybean has led to the rapid evolution of the virulence profile (pathotype) of P. sojae populations. Determining the pathotypes of P. sojae isolates is important in selecting soybean germplasm carrying the proper Rps, but this process is fastidious and requires specific expertise. In this work, we used a molecular assay to assess the pathotypes of P. sojae isolates obtained throughout the provinces of Québec, Ontario, and Manitoba. In preliminary assays, the molecular tool showed equivalent prediction of the pathotypes as a phenotyping assay and proved to be much faster to apply while eliminating intermediate values. Upon analysis of nearly 300 isolates, 24 different pathotypes were detected in Québec and Ontario, compared with only eight in Manitoba, where soybean culture is more recent. Pathotypes 1a, 1c, and 1d was predominant in Québec, while 1a, 1b, 1c, 1d, and 1k pathotypes were the most common in Manitoba. Overall, the results showed that 98 and 86% of the isolates carried pathotype 1a or 1c, respectively, suggesting that Rps1a and Rps1c were no longer effective in Canada. Based on the history of soybean varieties used in surveyed fields, it was found that 84% of them contained Rps genes that were no longer resistant against the pathotypes of the isolates found in the fields. While highlighting an easier and more precise option to assess pathotypes, this study presents the first pan-Canadian survey of P. sojae and stresses the importance of carefully managing the declining sources of resistance.

The Induction of the Isoflavone Biosynthesis Pathway Is Associated with Resistance to Common Bacterial Blight in Phaseolus vulgaris L

Xanthomonas axonopodis infects common bean (Phaseolus vulgaris L.) causing the disease common bacterial blight (CBB). The aim of this study was to investigate the molecular and metabolic mechanisms underlying CBB resistance in P. vulgaris. Trifoliate leaves of plants of a CBB-resistant P. vulgaris recombinant inbred line (RIL) and a CBB-susceptible RIL were inoculated with X. axonopodis or water (mock treatment). Leaves sampled at defined intervals over a 48-h post-inoculation (PI) period were monitored for alterations in global transcript profiles. A total of 800 genes were differentially expressed between pathogen and mock treatments across both RILs; approximately half were differentially expressed in the CBB-resistant RIL at 48 h PI. Notably, there was a 4- to 32-fold increased transcript abundance for isoflavone biosynthesis genes, including several isoflavone synthases, isoflavone 2'-hydroxylases and isoflavone reductases. Ultra-high performance liquid chromatography-tandem mass spectrometry assessed leaf metabolite levels as a function of the PI period. The concentrations of the isoflavones daidzein and genistein and related metabolites coumestrol and phaseollinisoflavan were increased in CBB-resistant RIL plant leaves after exposure to the pathogen. Isoflavone pathway transcripts and metabolite profiles were unaffected in the CBB-susceptible RIL. Thus, induction of the isoflavone pathway is associated with CBB-resistance in P. vulgaris.

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.

The Effects of Domestication on Secondary Metabolite Composition in Legumes

Legumes are rich in secondary metabolites, such as polyphenols, alkaloids, and saponins, which are important defense compounds to protect the plant against herbivores and pathogens, and act as signaling molecules between the plant and its biotic environment. Legume-sourced secondary metabolites are well known for their potential benefits to human health as pharmaceuticals and nutraceuticals. During domestication, the color, smell, and taste of crop plants have been the focus of artificial selection by breeders. Since these agronomic traits are regulated by secondary metabolites, the basis behind the genomic evolution was the selection of the secondary metabolite composition. In this review, we will discuss the classification, occurrence, and health benefits of secondary metabolites in legumes. The differences in their profiles between wild legumes and their cultivated counterparts will be investigated to trace the possible effects of domestication on secondary metabolite compositions, and the advantages and drawbacks of such modifications. The changes in secondary metabolite contents will also be discussed at the genetic level to examine the genes responsible for determining the secondary metabolite composition that might have been lost due to domestication. Understanding these genes would enable breeding programs and metabolic engineering to produce legume varieties with favorable secondary metabolite profiles for facilitating adaptations to a changing climate, promoting beneficial interactions with biotic factors, and enhancing health-beneficial secondary metabolite contents for human consumption.

Identification of Resistance Genes to Phytophthora sojae in Domestic Soybean Cultivars from China Using Particle Bombardment

Phytophthora root and stem rot caused by Phytophthora sojae is a destructive disease that afflicts soybean plants throughout the world. The use of resistant soybean cultivars is the primary means of managing this disease, as well as the most effective and economical approach. There are abundant soybean germplasm resources in China that could be deployed for breeding programs; however, the resistance genes (Rps genes) in most cultivars are unknown, leading to uncertainty concerning which are resistant cultivars for use. The resistance genes Rps1a, Rps1c, and Rps1k prevent root and stem rot caused by most P. sojae isolates within a Chinese field population. This study identified three Rps genes in Chinese domestic soybean cultivars using three related avirulence genes by particle bombardment. The complex genetic diversity of soybean cultivars and P. sojae strains has made it difficult to define single Rps genes without molecular involvement. Gene cobombardment is a method for identifying Rps genes quickly and specifically. We showed that cultivars Dongnong 60 and Henong 72 contained Rps1a, while Hedou 19, Henong 76, 75-3, Wandou 21020, Zheng 196, Wandou 28, Heinong 71, and Wandou 29 all contained Rps1c. The cultivars Jidou 12, Henong 72, Heinong 71, and Wandou 29 contained Rps1k. The cultivar Henong 72 contained both Rps1a and Rps1k, while Wandou 29 and Heinong 71 contained both Rps1c and Rps1k. We then evaluated the phenotype of 11 domestic soybean cultivars reacting to P. sojae using the isolates P6497 and Ps1. The 11 domestic cultivars were all resistant to P6497 and Ps1. This research provides source materials and parent plant strains containing Rps1a, Rps1c, and Rps1k for soybean breeding programs.

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.

Formation of Flavonoid Metabolons: Functional Significance of Protein-Protein Interactions and Impact on Flavonoid Chemodiversity

Flavonoids are a class of plant specialized metabolites with more than 6,900 known structures and play important roles in plant survival and reproduction. These metabolites are derived from p-coumaroyl-CoA via the sequential actions of a variety of flavonoid enzymes, which have been proposed to form weakly bound, ordered protein complexes termed flavonoid metabolons. This review discusses the impacts of the formation of flavonoid metabolons on the chemodiversity of flavonoids. Specific protein-protein interactions in the metabolons of Arabidopsis thaliana and other plant species have been studied for two decades. In many cases, metabolons are associated with the ER membrane, with ER-bound cytochromes P450 hypothesized to serve as nuclei for metabolon formation. Indeed, cytochromes P450 have been found to be components of flavonoid metabolons in rice, snapdragon, torenia, and soybean. Recent studies illustrate the importance of specific interactions for the efficient production and temporal/spatial distribution of flavonoids. For example, in diverse plant species, catalytically inactive type-IV chalcone isomerase-like protein serves as an enhancer of flavonoid production via its involvement in flavonoid metabolons. In soybean roots, a specific isozyme of chalcone reductase (CHR) interacts with 2-hydroxyisoflavanone synthase, to which chalcone synthase (CHS) can also bind, providing a mechanism to prevent the loss of the unstable CHR substrate during its transfer from CHS to CHR. Thus, diversification in chemical structures and temporal/spatial distribution patterns of flavonoids in plants is likely to be mediated by the formation of specific flavonoid metabolons via specific protein-protein interactions.

Genome-wide identification and localization of chalcone synthase family in soybean (Glycine max [L]Merr)

BACKGROUND

Soybean is a paleopolyploid that has undergone two whole genome duplication events. Gene duplication is a type of genomic change that can lead to novel functions of pre-existing genes. Chalcone synthase (CHS) is the plant-specific type III polyketide synthase that catalyzes the first committed step in (iso)flavonoid biosynthesis in plants.

RESULTS

Here we performed a genome-wide search of CHS genes in soybean, and identified 21 GmCHS loci containing 14 unique GmCHS (GmCHS1-GmCHS14) that included 5 newly identified GmCHSs (GmCHS10-GmCHS14). Furthermore, 3 copies of GmCHS3 and 2 copies of GmCHS4 were found in soybean. Analysis of gene structure of GmCHSs revealed the presence of a single intron in protein-coding regions except for GmCHS12 that contained 3 introns. Even though GmCHS genes are located on 8 different chromosomes, a large number of these genes are present on chromosome 8 where they form 3 distinct clusters. Expression analysis of GmCHS genes revealed tissue-specific expression pattern, and that some GmCHS isoforms localize in the cytoplasm and the nucleus while other isoforms are restricted to cytoplasm only.

CONCLUSION

Overall, we have identified 21 GmCHS loci with 14 unique GmCHS genes in the soybean genome. Their gene structures and genomic organization together with the spatio-temporal expression and protein localization suggest their importance in the production of downstream metabolites such as (iso)flavonoids and their derived phytoalexins.

Isoflavonoid-specific prenyltransferase gene family in soybean: GmPT01, a pterocarpan 2-dimethylallyltransferase involved in glyceollin biosynthesis

Phytoalexin glyceollins are soybean-specific antimicrobial compounds that are derived from the isoflavonoid pathway. They are synthesized by soybean in response to extrinsic stress such as pathogen attack or injury, thereby conferring partial resistance if synthesized rapidly at the site of infection and at the required concentration. Soybean produces multiple forms of glyceollins that result from the differential prenylation reaction catalyzed by prenyltransferases (PTs) on either the C-2 or C-4 carbon of a pterocarpan glycinol. The soybean genome contains 77 PT-encoding genes (GmPTs) where at least 11 are (iso)flavonoid-specific. Transcript accumulation of five candidates GmPTs was increased in response to Phytophthora sojae infection, suggesting their role in phytoalexin synthesis. The induced GmPTs localize to plastids and display tissue-specific expression. We have in this study identified two additional GmPTs: an isoflavone dimethylallyltransferase 3 (IDT3); and a glycinol 2-dimethylallyl transferase GmPT01. GmPT01 prenylates (-)-glycinol at the C-2 position, localizes in the plastid, and exhibits root-specific gene expression. Furthermore, its expression is induced rapidly in response to stress, and is associated with a quantitative trait loci linked with resistance to P. sojae. Based on these results, we conclude that GmPT01 are possibly one of the loci involved in conferring partial resistance against stem and root rot disease in soybean.

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

Soybean (Glycine max) 5-deoxyisoflavonoids (daidzein and its conjugates) are precursors of glyceollin phytoalexins. They are also converted to equol by microbes in the human intestine, resulting in health benefits. 5-Deoxyisoflavonoids accumulate in the roots (93% mol/mol of the total root isoflavonoids) and seeds of unstressed soybean plants. Chalcone reductase (CHR) is a key enzyme mediating 5-deoxyisoflavonoid biosynthesis because it catalyzes the production of 6'-deoxychalcone through its effects on the chalcone synthase (CHS)-catalyzed reaction. The soybean genome encodes at least 11 CHR-related homologs, but it is unclear which ones are functionally important for daidzein accumulation in unstressed plants. Among the CHR homologs, the temporal and spatial expression patterns of GmCHR5 were the most correlated with the distribution patterns of 5-deoxyisoflavonoids. The CHR activity of GmCHR5 was confirmed in vitro and in planta. In the in vitro assays, the ratio of CHR products (6'-deoxychalcone) to total CHS products (R value) was dependent on GmCHR5 and CHS concentrations, with higher concentrations resulting in higher R values (i.e. approaching 90%). Subcellular localization analyses revealed that GmCHR5 was present in the cytoplasm and nucleus. Protein-protein interaction assays indicated that GmCHR5, but not GmCHR1 and GmCHR6, interacted with 2-hydroxyisoflavanone synthase (IFS) isozymes. The CHS isozymes also interacted with IFS isozymes but not with GmCHR5. The proposed micro-compartmentalization of isoflavone biosynthesis through the formation of an IFS-mediated metabolon is probably involved in positioning GmCHR5 close to CHS, resulting in an R value that is high enough for the accumulation of abundant 5-deoxyisoflavonoids in soybean roots.