The research team, led by Xinkang Feng, Daqian Sun, and senior authors Lijuan Qiu, Xiaobo Wang, and Yinghui Li, began with an unusually deep phenotyping effort. They grew 311 cultivated soybean accessions in Beijing over three consecutive growing seasons, from 2009 to 2011, and measured six distinct isoflavone components in the harvested seeds. The data revealed a clear hierarchy: glycoside-type isoflavones such as genistin, glycitin, and daidzin accumulated at far higher levels than their aglycone counterparts, with genistin alone accounting for 32.75 percent of total isoflavones while genistein made up just 4.23 percent. Statistical analysis using best linear unbiased prediction showed that glycoside-type components correlated strongly with total isoflavone content, with Pearson coefficients between 0.71 and 0.89, and that total content and genistin were substantially heritable, with broad-sense heritability estimates of 0.52 and 0.57 respectively. Aglycones, by contrast, showed heritability of only 0.12 to 0.30, suggesting they are far more sensitive to environmental conditions.

With a reliable phenotypic foundation in place, the researchers turned to a dual-pronged genetic strategy. First, a genome-wide association study scanned 4,765,428 single nucleotide polymorphisms across the panel using the FarmCPU model, with population structure controlled through principal components and significance assessed by Bonferroni correction. This analysis mapped 22 quantitative trait loci for isoflavone traits, seven of which overlapped with previously reported regions, lending confidence to the results. A particularly strong signal converged at the distal end of chromosome 5, where association signals for four different isoflavones clustered within a 74-kilobase window in a region the team designated qISO.Gm05. Second, a transcriptome-wide association study examined exon-level expression data from developing seeds at the R5.5 stage, testing 154,119 exons across 204 accessions. Among 71 significant genes, one stood out: Glyma.05G244100, repeatedly detected for four isoflavones, whose exon 2 usage correlated positively with seed isoflavone content and which sat squarely within the major qISO.Gm05 locus. The team named this convergence candidate GmISO5.

What makes GmISO5 remarkable is its identity. The gene encodes a member of the phosphatidylethanolamine-binding protein family, related to the well-known MOTHER OF FT AND TFL1, or MFT, proteins that govern flowering time, seed development, and germination in diverse plants. In soybean, this locus had previously been characterized as GmMFT, or GmST05, a major determinant of seed oil and protein content. No PEBP-family protein had ever been directly linked to isoflavone biosynthesis, making GmISO5 an unexpected addition to a regulatory repertoire previously dominated by MYB transcription factors and zinc-finger proteins. Expression analysis confirmed that GmISO5 is specifically active in developing seeds, peaking around 30 days after flowering, with the strongest signal in the cotyledon and embryonic axis rather than the seed coat.

To test whether GmISO5 truly drives isoflavone accumulation, the researchers generated two independent CRISPR/Cas9 knockout lines in the Wm82 background, which naturally carries the low-isoflavone allele. Both lines carried small frameshift deletions in exon 1 that introduced premature stop codons. The metabolic consequences were dramatic. Total seed isoflavones plummeted from 1,375.52 micrograms per gram in wild-type seeds to 273.67 micrograms per gram in the mutants, a reduction of 80.1 percent, with steep declines in glycitin, daidzin, daidzein, glycitein, and genistein. RNA sequencing of developing seeds revealed thousands of differentially expressed genes, and strikingly, 30 of the 65 known isoflavone biosynthetic genes were among those regulated by GmISO5, indicating that this single protein exerts broad transcriptional control over the entire pathway network rather than tweaking a single step.

Molecular dissection pinpointed two direct targets. GmCHS1 encodes chalcone synthase, the enzyme that catalyzes an early, committed step in flavonoid and isoflavonoid production, while GmHIDH encodes 2-hydroxyisoflavanone dehydratase, which operates in the isoflavone-specific branch to produce the aglycones daidzein and genistein. Yeast one-hybrid assays showed that GmISO5 associates with both promoters, and dual-luciferase reporter assays in tobacco leaves demonstrated that the protein transactivates them, producing markedly stronger luminescence signals than controls. Electrophoretic mobility shift assays then confirmed direct, sequence-specific binding: recombinant GmISO5 bound two fragments of the GmHIDH promoter and one fragment of the GmCHS1 promoter, with binding competed away by excess unlabeled probe and abolished by mutations. The three bound fragments shared a distinctive A/T-rich arrangement, TTAAAA separated by four to six nucleotides from ATAAT, which the authors propose may contribute to promoter recognition, though further work is needed to confirm it as a general binding motif.

Population genetics added an evolutionary twist. Analyzing sequence variation in 1,566 re-sequenced soybean accessions, the team found that GmISO5 variants fall into two major haplotypes: GmISO5-High, associated with elevated isoflavone accumulation, and GmISO5-Low, associated with reduced levels. The frequency of the low-isoflavone haplotype rose sharply from wild soybean through landraces to improved cultivars, reaching 91 percent in modern varieties. Nucleotide diversity in the surrounding 400-kilobase region declined progressively from wild soybeans to cultivars, Tajima’s D dropped in improved material, and fixation index values rose between populations, all signatures of a selective sweep. Intriguingly, accessions carrying the low haplotype also had significantly higher 100-seed weight and seed oil content, suggesting that breeders selecting for larger, oilier seeds inadvertently dragged the low-isoflavone allele along, depleting the high-isoflavone haplotype as a correlated consequence rather than by design.

Co-expression network analysis using independent component analysis reinforced the picture. Among 58 modules, one, designated IC37, showed the strongest association with isoflavone content and contained GmISO5 alongside its validated target GmHIDH and core pathway genes including CHS1, CHS11, CHR, CHI, 4CL, UGT72B8, and TT2. Subdividing IC37 revealed that Submodule 2, specifically enriched for flavonoid and isoflavonoid biosynthesis, correlated most strongly with seed isoflavone levels, marking it as the functional core of the network. The authors caution that co-expression evidence is correlative and that the module should be treated as a prioritization tool for future candidate-gene validation rather than proof of regulatory relationships.

The study’s most provocative finding concerns viral disease. Because isoflavones have been implicated in plant-pathogen interactions, and because multiplex CRISPR editing of flavonoid pathway genes has previously been shown to increase both isoflavone content and resistance to soybean mosaic virus, the team inoculated wild-type and GmISO5 knockout plants with SMV. Transcript levels of GmISO5 and GmHIDH rose in wild-type leaves two days after inoculation, but the knockout lines developed more severe symptoms and accumulated significantly more virus, as quantified by DAS-ELISA, with leaf isoflavone levels correspondingly reduced. The result implicates GmISO5 in the soybean response to SMV infection, yet the authors are careful to note that the experiments do not establish causation: reduced isoflavones may not be the direct cause of heightened susceptibility, and other GmISO5-regulated defense pathways could contribute. Genetic complementation, overexpression lines, and metabolite-rescue experiments will be needed to resolve the mechanism.

Taken together, the findings expand the known regulatory architecture of soybean isoflavone biosynthesis into an entirely new protein family and demonstrate the power of integrating genome-wide and transcriptome-wide association studies with network analysis and genome editing. For breeders, the rarity of the GmISO5-High haplotype in major production regions of the United States and northern China highlights an underutilized genetic resource for boosting the nutritional quality of soybean seeds, though any deployment must reckon with the allele’s links to seed size and oil traits that breeders have long favored. For plant virologists, the connection between a metabolic regulator and SMV susceptibility opens a fresh line of inquiry into how specialized metabolism shapes antiviral defense, one that the authors themselves flag as the study’s most important open question.

Subject of Research: Genetic regulation of seed isoflavone biosynthesis and soybean mosaic virus response in soybean

Article Title: Multi-omics integration identifies GmISO5 as a key regulator of seed isoflavone biosynthesis and implicates it in the response to soybean mosaic virus

Article References: Feng, X., Sun, D., Liu, H., Chu, L., Li, J., Li, D., Zhang, H., Wang, X., Ge, T., Zheng, H., Guo, S., Zhao, X., Li, J., Zhao, T., Li, Y., Wang, X., & Qiu, L. (2026). Multi-omics integration identifies GmISO5 as a key regulator of seed isoflavone biosynthesis and implicates it in the response to soybean mosaic virus. Journal of Advanced Research. https://doi.org/10.1016/j.jare.2026.09.011

Image Credits: AI Generated

DOI: 10.1016/j.jare.2026.09.011

Keywords: soybean, isoflavones, GmISO5, GWAS, TWAS, CRISPR, soybean mosaic virus, PEBP family, selective sweep, seed quality, chalcone synthase, plant virology