Search PubMedSearch

PubMed · 42706246

Analysis of deep-resequencing data of 984 soybean accessions reveals structural variations underlying agronomic traits.

Abstract

Genomic structural variants (SVs) are major sources of genetic variation and have profound impacts on phenotypic traits. However, their functional effects remain largely unexplored in soybean. Here, we resequence 940 soybean accessions. Together with 44 publicly available datasets, we identify 602,281 SVs. Using a graph-based genome, we detect an additional 58,760 presence/absence variations (PAVs) that broadly affect gene expression. Population genomic analyses reveal that SVs serve as a core driving force for soybean domestication and improvement. Integrating SVs with QTLs for oil and protein content, and performing GWAS on 27 traits, we identify key functional SVs. These include transposable element insertions altering seed coat color, multiple insertions within a cytochrome P450 gene modifying flower and hypocotyl color, and a GmMATE1 deletion enhancing seed size. Together, our study establishes a comprehensive SV map of soybean, offering a valuable resource for dissecting the genetic basis of complex traits to accelerate molecular breeding.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Jie Li, Wenhao Yue, Guangqi He, Yubin Zou, Weidong Wang, James Whelan, Junmei Wang, Xianzhong Feng, Huixia Shou. 2026-07-29. Analysis of deep-resequencing data of 984 soybean accessions reveals structural variations underlying agronomic traits.. https://doi.org/10.1038/s41467-026-76055-9

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

GmMYB29 activates Gm4CL3 to enhance soybean resistance to Heterodera glycines.

Soybean cyst nematode is a devastating soil-borne pathogen that severely limits soybean yield worldwide. To uncover downstream target genes of the resistance-associated transcription factor GmMYB29, we combined ChIP-seq and RNA-seq data from T3-generation GmMYB29-overexpressing soybean plants, alongside Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, to screen candidate genes carrying transcription factor binding peaks within the 2000 bp region upstream of transcription start sites (TSS). Four orthogonal molecular assays-yeast one-hybrid (Y1H), electrophoretic mobility shift assay (EMSA), dual-luciferase reporter (LUC) system, and GUS histochemical staining-collectively confirmed the specific physical interaction between GmMYB29 and the promoter of Gm4CL3. We generated transgenic soybean hairy roots overexpressing Gm4CL3 (OX-Gm4CL3) and CRISPR-Cas9-mediated Gm4CL3 knockout lines (KO-Gm4CL3), with wild-type (WT) plants serving as controls. Inoculation assays using SCN 3 demonstrated that OX-Gm4CL3 roots displayed substantially improved SCN resistance, while KO-Gm4CL3 roots were hypersusceptible to nematode infection. Mechanistic investigations revealed that Gm4CL3 promotes lignin deposition in root tissues to block SCN penetration. Furthermore, GmMYB29 and Gm4CL3 act synergistically to activate lignin biosynthetic pathways and strengthen soybean resistance against SCN 3 (SCN Race 3, the dominant physiological race in Northeast China). In summary, this study functionally characterizes Gm4CL3 and defines a previously unreported GmMYB29-Gm4CL3 regulatory cascade that mediates plant defense against SCN. This module functions independent of classic SCN resistance loci rhg1/Rhg4, providing new genetic resources for SCN-resistant soybean molecular breeding.

Glycine max

Genome-wide identification, characterization, evolutionary analysis, and expression profiling of the FCS-like zinc finger (FLZ) gene family in soybean (Glycine max L.) under abiotic stresses.

Drought and salinity limit soybean yield. Despite their role in the SnRK1 energy-sensing complex, a systematic study of FCS-Like Zinc Finger (FLZ) proteins in soybean has not been reported. We performed a genome-wide identification of the GmFLZ gene family, identifying 40 members distributed across 18 of the 20 soybean chromosomes. Phylogenetic analysis of 87 FLZ proteins from Glycine max, Arabidopsis thaliana, and Oryza sativa revealed four major evolutionary clades, suggesting that diversification predates the separation of monocots and dicots. Structural analysis identified ten conserved motifs, with Motifs 1 and 2 present in all family members. Gene duplication analysis identified 304 paralogous pairs, most arising from segmental duplication. Ka/Ks analysis indicated localized positive selection in six gene pairs and purifying selection in 97.9% of pairs. Tissue-specific expression profiling across nine tissues showed that GmFLZ5, GmFLZ15, GmFLZ25, and GmFLZ34 had the highest expression levels detected across the GmFLZ family, with GmFLZ5 the most highly expressed member in leaves, nodules, and stem and showing moderate expression in pod, root, and root hairs, whereas GmFLZ18, GmFLZ23, and GmFLZ37 showed root-preferential expression. RT-qPCR validation under drought (20% PEG-6000) and salt (200 mM NaCl) treatments in the Giza 5 cultivar showed that 36 and 34 of the 40 GmFLZ genes, respectively, exhibited at least a two-fold change in expression, with GmFLZ21 and GmFLZ35 among the most strongly induced under salt stress. These findings provide an evolutionary and functional framework for the GmFLZ family and identify candidate genes for future functional studies in soybean stress tolerance.

Glycine max

Transgenic overexpression of GmAPC7-CT improves seed yield and reduces susceptibility to soybean mosaic virus and Meloidogyne incognita in soybean.

Stable transgenic soybean lines overexpressing the GmAPC7-CT gene have demonstrated increased seed yield and reduced susceptibility to the soybean mosaic virus and Meloidogyne incognita. The Anaphase-Promoting Complex subunit 7 (APC7) is a core structural component of the anaphase-promoting complex or cyclosome (APC/C). The terminal region of this AtAPC7 gene has been shown in Arabidopsis thaliana to accumulate more transcripts than the full-length gene. The AtAPC7-CT gene (terminal region of the AtAPC7) encodes a protein with significant homology to a tobacco viral replication inhibitor (IVR). Its stable overexpression in transgenic A. thaliana lines resulted in notable improvements in biomass, seed yield, earliness of vegetative-reproductive transitions, and reduced susceptibility to viruses. In this study, we generated stable transgenic soybean lines overexpressing the GmAPC7-CT gene (terminal region or 3' portion of Glyma.15G096000, corresponding to the AtAPC7-CT) and evaluated seed yield and susceptibility of these lines to soybean mosaic virus and Meloidogyne incognita. The GmAPC7-CT gene is 624 nucleotides long and encodes a 207-amino acid protein with two tetratricopeptide repeat (TPR) domains. GmAPC7-CT showed 100% amino acid identity with full-length GmAPC7, 81.16% identity with AtAPC7-CT, and 87.94% identity with tobacco IVR. Stable transgenic lines demonstrated significant advancements in plant development and seed yield, with the top three lines producing up to 43% more pods, 44% more seeds, and a 16% increase in seed weight. Furthermore, these soybean lines showed up to a 70% reduction in susceptibility to soybean mosaic virus and M. incognita, reflected by decreased viral RNA load and nematode reproduction factor. Collectively, these results support a conserved role of GmAPC7-CT in soybean and AtAPC7-CT in A. thaliana, acting similarly to the tobacco IVR. Thus, our findings underscore the strong biotechnological potential of the GmAPC7-CT gene to improve key agronomic traits in soybean through genetic engineering approaches, including conventional breeding, transgenesis, and genome editing.

Glycine max