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.

In soybeans (Glycine max (L.) Merr.), their growth periods, DSF (days of sowing-to-flowering), and DFM (days of flowering-to-maturity) are determined by their required accumulative day-length (ADL) and active temperature (AAT). A sample of 354 soybean varieties from five world eco-regions was tested in four seasons in Nanjing, China. The ADL and AAT of DSF and DFM were calculated from daily day-lengths and temperatures provided by the Nanjing Meteorological Bureau. The improved restricted two-stage multi-locus genome-wide association study using gene–allele sequences as markers (coded GASM-RTM-GWAS) was performed. (i) For DSF and its related ADL_DSF and AAT_DSF, 130–141 genes with 384–406 alleles were explored, and for DFM and its related ADL_DFM and AAT_DFM, 124–135 genes with 362–384 alleles were explored, in a total of six gene–allele systems. DSF shared more ADL and AAT contributions than DFM. (ii) Comparisons between the eco-region gene–allele submatrices indicated that the genetic adaptation from the origin to the geographic sub-regions was characterized by allele emergence (mutation), while genetic expansion from primary maturity group (MG)-sets to early/late MG-sets featured allele exclusion (selection) without allele emergence in addition to inheritance (migration). (iii) Optimal crosses with transgressive segregations in both directions were predicted and recommended for breeding purposes, indicating that allele recombination in soybean is an important evolutionary drive. (iv) Genes of the six traits were mostly trait-specific involved in four categories of 10 groups of biological functions. GASM-RTM-GWAS showed potential in detecting directly causal genes with their alleles, identifying differential trait evolutionary drives, predicting recombination breeding potentials, and revealing population gene networks.