Identification and characterization of a major QTL underlying soybean isoflavone malonylglycitin content

Joint GWAS and WGCNA Identify Genes Regulating the Isoflavone Content in Soybean Seeds

Isoflavone is a secondary metabolite of the soybean phenylpropyl biosynthesis pathway with physiological activity and is beneficial to human health. In this study, the isoflavone content of 205 soybean germplasm resources from 3 locations in 2020 showed wide phenotypic variation. A joint genome-wide association study (GWAS) and weighted gene coexpression network analysis (WGCNA) identified 33 single-nucleotide polymorphisms and 11 key genes associated with soybean isoflavone content. Gene ontology enrichment analysis, gene coexpression, and haplotype analysis revealed natural variations in the Glyma.12G109800 (GmOMT7) gene and promoter region, with Hap1 being the elite haplotype. Transient overexpression and knockout of GmOMT7 increased and decreased the isoflavone content, respectively, in hairy roots. The combination of GWAS and WGCNA effectively revealed the genetic basis of soybean isoflavone and identified potential genes affecting isoflavone synthesis and accumulation in soybean, providing a valuable basis for the functional study of soybean isoflavone.

Cytochrome GmGLY1 is Involved in the Biosynthesis of Glycitein in Soybean

Isoflavones, the major secondary metabolites of interest due to their benefits to both human and plant health, are exclusively produced by legumes. In this study, we profiled the isoflavone content in dry seeds from 211 soybean [Glycine max (L.) Merr.] accessions grown across five environments. Broad and discernible phenotypic variations were observed among accessions, regions, and years of growth. Twenty-six single-nucleotide polymorphisms (SNPs) associated with the sum of glycitein (GLE), glycitin (GL), 6″-O-acetylglycitin (AGL), and 6″-O-malonylglycitin (MGL) contents were detected in multiple environments via a genome-wide association study (GWAS). These SNPs were located on chromosome 11 (8,148,438 bp to 8,296,956 bp, renamed qGly11-01). Glyma.11g108300 (GmGLY1), a gene that encodes a P450 family protein, was identified via sequence variation analysis, functional annotation, weighted gene coexpression network analysis (WGCNA), and expression profile analysis of candidate gene, and hairy roots transformation in soybean. Overexpression of GmGLY1 increased the glycitein content (GLC) in soybean hairy roots and transgenic seeds, while CRISPR/Cas9-generated mutants exhibited decreased GLC and increased daidzein content (DAC). Haplotype analysis revealed that GmGLY1 allelic variations significantly affect the GLC accumulation. These findings enhance our understanding of genes influencing GLC in soybean and may guide breeding for lines with high and stable GLC.

Isoflavonoid metabolism in leguminous plants: an update and perspectives

Isoflavonoids constitute a well-investigated category of phenylpropanoid-derived specialized metabolites primarily found in leguminous plants. They play a crucial role in legume development and interactions with the environment. Isoflavonoids usually function as phytoalexins, acting against pathogenic microbes in nature. Additionally, they serve as signaling molecules in rhizobial symbiosis. Notably, owing to their molecular structure resembling human estrogen, they are recognized as phytoestrogens, imparting positive effects on human health. This review comprehensively outlines recent advancements in research pertaining to isoflavonoid biosynthesis, transcriptional regulation, transport, and physiological functions, with a particular emphasis on soybean plants. Additionally, we pose several questions to encourage exploration into novel contributors to isoflavonoid metabolism and their potential roles in plant-microbe interactions.

Identification of Genes Responsible for the Synthesis of Glycitein Isoflavones in Soybean Seeds

Soybean (Glycine max (L.) Merrill) isoflavones are among the most important secondary metabolites, with functional benefits for human health. Soybeans accumulate three aglycone forms of isoflavones: genistein, daidzein, and glycitein. Soybean landrace Kumachi-1 does not accumulate malonylglycitin at all. Gene structure analysis indicated that Glyma.11G108300 (F6H4) of Kumachi-1 has a 3.8-kbp insertion, resulting in a truncated flavonoid 6-hydroxylase (F6H) sequence compared to the wild-type sequence in Fukuyutaka. Mapping experiments using a mutant line (MUT1246) with a phenotype similar to that of Kumachi-1, with a single-nucleotide polymorphism (SNP) in F6H4, revealed co-segregation of this mutation and the absence of glycitein isoflavones. We also identified a mutant line (K01) that exhibited a change in the HPLC retention time of glycitein isoflavones, accumulating glycoside and malonylglycoside forms of 6-hydroxydaidzein. K01 contains an SNP that produces a premature stop codon in Glyma.01G004200 (IOMT3), a novel soybean isoflavone O-methyltransferase (IOMT) gene. We further analyzed transgenic hairy roots of soybeans expressing Glyma.11G108300 (F6H4) and Glyma.01G004200 (IOMT3). Those overexpressing F6H4 accumulated malonylglycoside forms of 6-hydroxydaidzein (M_6HD), and co-expression of F6H4 and IOMT3 increased the level of malonylglycitin but not of M_6HD. These results indicate that F6H4 and IOMT3 are responsible for glycitein biosynthesis in soybean seed hypocotyl.

Modified High-Resolution Melting (HRM) Marker Systems Increasing Discriminability Between Homozygous Alleles

Targeted single-nucleotide polymorphism (SNP) genotyping, especially for functional nucleotide polymorphism, is widely used for current breeding programs in crops. One of the cost- and time-effective approaches for genotyping is high-resolution melting (HRM) analysis for polymerase chain reaction (PCR) amplicons, including target SNP. The reliability of a genotype obtained from an HRM marker depends on the difference in T_m values between two amplicons. Increasing the reliability of HRM marker genotypes could be archived with the selection of the best nearest neighboring nucleotide substitution (NNNs) in primer sequences surrounding SNPs. This chapter provides an easy-way protocol to design primer sequences for NNNs-HRM markers with table and web service, as well as several tips to develop HRM markers that distinguish between homozygous alleles (e.g., between A/A and C/C).

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Effect of Kefir on Soybean Isoflavone Aglycone Content in Soymilk Kefir

Kefir is a traditional fermented milk originating in the Caucasus area and parts of Eastern Europe. In this study, the kefir culture, which is modified upon the addition of lactic acid bacteria (LAB) cells, specifically for soymilk kefir fermentation with the highest capacity of isoflavone biotransformation, was successfully produced, and the metagenomics composition of soymilk or milk fermented using these kefir cultures was investigated. The metagenome analysis showed that the microbiota of kefir in M-K (milk inoculated with kefir), SM-K (equal volumes of soymilk and milk inoculated with kefir), and S-K (pure milk inoculated with kefir) were related to the addition of soymilk or not. Furthermore, the HPLC chromatogram revealed that Guixia 2 (Guangzhou, China) may be a good source of soymilk kefir fermentation due to its high isoflavone aglycone content (90.23 ± 1.26 μg/g in daidzein, 68.20 ± 0.74 μg/g in genistein). Importantly, the starter culture created by adding 1.5 g probiotics (Biostime®, Guangzhou, China) to Chinese kefir showed a significant increase in the levels of isoflavone aglycones (72.07 ± 0.53 μg/g in isoflavone aglycones). These results provided insight into understanding the suitable soybean cultivar and starter cultures, which exhibit promising results of isoflavone biotransformation and flavor promotion during soymilk kefir fermentation.

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.

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.