Plant height (PH) is an important agronomic trait and is closely related to yield in soybean [Glycine max (L.) Merr.]. Previous studies have identified many QTLs for PH. Due to the complex genetic background of PH in soybean, there are few reports on its fine mapping. In this study, we used a mapping population derived from a cross between a chromosome segment substitution line CSSL3228 (donor N24852 (G. Soja), a receptor NN1138–2 (G. max)) and NN1138–2 to fine map a wild soybean allele of greater PH by QTL-seq and linkage mapping. We identified a QTL for PH in a 1.73 Mb region on soybean chromosome 13 through QTL-seq, which was confirmed by SSR marker-based classical QTL mapping in the mapping population. The linkage analysis showed that the QTLs of PH were located between the SSR markers BARCSOYSSR_13_1417 and BARCSOYSSR_13_1421 on chromosome 13, and the physical distance was 69.3 kb. RT-PCR and sequence analysis of possible candidate genes showed that Glyma.13 g249400 revealed significantly higher expression in higher PH genotypes, and the gene existed 6 differences in the amino acids encoding between the two parents. Data presented here provide support for Glyma.13 g249400 as a possible candidate genes for higher PH in wild soybean line N24852.

Drought is one of the most important factors affecting plant growth and productivity. The previous results on drought tolerance (DT) genetic system in soybean indicated a complex of genes not only few ones were involved in the trait. This study is featured with a relatively thorough identification of QTL-allele/candidate-gene system using an efficient restricted two-stage multi-locus multi-allele genome-wide association study, on two comprehensive DT indicators, membership index values of relative plant weight (MPW) and height (MPH), instead of a single biological characteristic, in a large sample (564 accessions) of the Chinese cultivated soybean population (CCSP). Based on 24,694 multi-allele markers, 75 and 64 QTL with 261 and 207 alleles (2–12/locus) were detected for MPW and MPH, explaining 54.7% and 47.1% of phenotypic variance, respectively. The detected QTL-alleles were organized into a QTL-allele matrix for each indicator, indicating DT is a super-trait conferred by two (even more) QTL-allele systems of sub-traits. Each CCSP matrix was separated into landrace (LR) and released cultivar (RC) sub-matrices, which showed significant differentiation in QTL-allele constitutions, with 58 LR alleles excluded and 16 new ones emerged in RC. Using the matrices, optimal crosses with great DT transgressive recombinants were predicted. From the detected QTL, 177 candidate genes were annotated and validated with quantitative Real-time PCR, and grouped into nine categories, with ABA and stress responders as the major parts. The key point of the above results is the establishment of relatively full QTL-allele matrices composed of numerous gene functions jointly conferring DT, therefore, demonstrates the complexity of DT genetic system and potential of CCSP in DT breeding.

Annual wild soybean (G. soja) is the ancestor of the cultivated soybean (G. max). To reveal the genetic changes from soja to max, an improved wild soybean chromosome segment substitution line (CSSL) population, SojaCSSLP5, composed of 177 CSSLs with 182 SSR markers (SSR-map), was developed based on SojaCSSLP1 generated from NN1138-2(max)×N24852(soja). The SojaCSSLP5 was genotyped further through whole-genome resequencing, resulting in a physical map with 1366 SNPLDBs (SNP linkage-disequilibrium blocks), which are composed of more markers/segments, shorter marker length and more recombination breakpoints than the SSR-map and caused 721 new wild substituted segments. Using the SNPLDB-map, two loci co-segregating with seed-coat color (SCC) and six loci for days to flowering (DTF) with 88.02% phenotypic contribution were identified. Integrated with parental RNA-seq and DNA-resequencing, two SCC and six DTF candidate genes, including three previously cloned (G, E2 and GmPRR3B) and five newly detected ones, were predicted and verified at nucleotide mutant level, and then demonstrated with the consistency between gene-alleles and their phenotypes in SojaCSSLP5. In total, six of the eight genes were identified with the parental allele-pairs coincided to those in 303 germplasm accessions, then were further demonstrated by the consistency between gene-alleles and germplasm phenotypes. Accordingly, the CSSL population integrated with parental DNA and RNA sequencing data was demonstrated to be an efficient platform in identifying candidate wild vs. cultivated gene-alleles.

Salinization of cultivated soils may result in either high salt levels or alkaline conditions, both of which stress crops and reduce performance. We sampled genotypes included in the Northeast China soybean germplasm population (NECSGP) to identify possible genes that affect tolerance to alkaline soil conditions. In this study, 361 soybean accessions collected in Northeast China were tested under 220 mM NaHCO₃:Na₂CO₃ = 9:1 (pH = 9.8) to evaluate the alkali-tolerance (ATI) at the seedling stage in Mudanjiang, Heilongjiang, China. The restricted two-stage multi-locus model genome-wide association study (RTM-GWAS) with gene–allele sequences as markers (6503 GASMs) based on simplified genome resequencing (RAD-sequencing) was accomplished. From this analysis, 132 main effect candidate genes with 359 alleles and 35 Gene × Environment genes with 103 alleles were identified, explaining 90.93% and 2.80% of the seedling alkali-tolerance phenotypic variation, respectively. Genetic variability of ATI in NECSGP was observed primarily within subpopulations, especially in ecoregion B, from which 80% of ATI-tolerant accessions were screened out. The biological functions of 132 candidate genes were classified into eight functional categories (defense response, substance transport, regulation, metabolism-related, substance synthesis, biological process, plant development, and unknown function). From the ATI gene–allele system, six key genes–alleles were identified as starting points for further study on understanding the ATI gene network.

Shade tolerance is a key trait for cultivars in inter/relay-cropped soybeans in maize fields. Our previous genome-wide association study (GWAS) results on southern China soybean germplasm revealed that the shade tolerance was conferred by a complex of genes with multiple alleles. To complete our understanding of the shade tolerance gene system, GWAS with gene–allele sequences as markers (designated GASM-RTM-GWAS) was conducted in a recombinant inbred line (RIL) population between two extreme parents using the shade tolerance index (STI) and relative pith cell length (RCL) as indicators. Altogether, 211 genes, comprising 99 and 119 genes (seven shared) for STI and RCL, respectively, were identified and then annotated into a similar set of five biological categories. Furthermore, transcriptome analysis detected 7837 differentially expressed genes (DEGs), indicating plentiful DEGs involved in the expression of regulatory/causal GWAS genes. Protein–protein interaction (PPI) analysis and gene functional analysis for both GWAS genes and DEGs showed a group of interrelated causal genes and a group of interrelated DEGs; the former were included in the latter and their functions were interconnected as a gene network. For further understanding of the response of soybean to shade stress in a sequential connection, six chronological gene modules were grouped as signal activation and transport, signal-transduction, signal amplification, gene expression, regulated metabolites, and material transport. From the modules, 12 key genes were selected as entry points for further analysis. Our study provides an overview of the shade tolerance gene network as a new insight into a complex-trait genetic system, rather than the usual way of starting from a hand-picked single gene.

Annual wild soybean, the ancestor of cultivated soybean, underwent a significant increase in seed oil content during domestication. To elucidate the genetic basis of this change, a chromosome segment substitution line population (177 lines) constructed with cultivated soybean NN1138-2 as recipient and wild soybean N24852 as donor was used in this study. Phenotypic evaluation across three distinct environments led to the identification of two major QTL/segments, qOC14 on chromosome 14 and qOC20 on chromosome 20, which collectively explained 39.46% of the phenotypic variation, with individual contributions of 17.87% and 21.59%, respectively. Both wild alleles exhibited negative additive effects, with values of −0.35% and −0.42%, respectively, consistent with the inherently low oil content of wild soybeans. Leveraging transcriptome and genome data from the two parents, two candidate genes were predicted. Notably, Glyma.14G179800 is a novel candidate gene encoding a PHD-type zinc finger domain-containing protein, and the hap-A haplotype exhibits a positive effect on oil content. In contrast, Glyma.20G085100 is a reported POWR1 gene, known to regulate protein and oil content. Our findings not only validate the role of known gene but, more importantly, unveil a new candidate gene, offering valuable genetic resources and theoretical targets for molecular breeding of high-oil soybean.