Integrative Multi-Omics Analysis of Stem Growth Habit Divergence in Wild Soybean (Glycine soja)

Stem architecture is a major determinant of lodging resistance, biomass accumulation, and harvest efficiency in soybean. However, the molecular features associated with contrasting stem growth habits in wild soybean remain incompletely characterised. Here, we performed an integrated transcriptomic, metabolomic, and epigenomic analysis of stem growth-habit divergence in wild soybean, comparing the wild-type accession ZYD7068 with contrasting vining and erect mutant lines derived from carbon-ion beam mutagenesis. Pairwise transcriptomic comparisons identified between 20 311 and 28 705 differentially expressed genes per contrast, with a core set of 2672 genes consistently altered across the comparisons. Functional enrichment, gene set variation analysis, and gene set enrichment analysis converged on xylem and phloem pattern formation as a prominent molecular pathway associated with growth-habit divergence. Random forest analysis identified BBR-BPC and ARF transcription factor families as major molecular discriminators, while metabolomic profiling revealed distinct metabolic profiles involving amino-acid-derived and lipid-associated metabolites. Whole-genome bisulfite sequencing revealed context-specific DNA methylation differences, including substantial variation in CHG methylation among erect mutant lines. Integrated network and in silico perturbation analyses prioritised four candidate genes associated with vascular development for future functional validation. Together, these results provide a multi-layer molecular resource for investigating stem growth-habit divergence in G. soja and establish testable candidate pathways and genes for subsequent functional studies and soybean improvement.