Gene-allele system of shade tolerance in southern China soybean germplasm revealed by genome-wide association study using gene-allele sequence as markers.

KEY MESSAGE: Fifty-three shade tolerance genes with 281 alleles in the SCSGP were identified directly using gene-allele sequence as markers in RTM GWAS, from which optimized crosses, evolutionary motivators, and gene-allele networks were explored. Shade tolerance is a key for optimal cultivation of soybean inter/relay-cropped with corn. To explore the shade tolerance gene-allele system in the southern China soybean germplasm, we proposed using gene-allele sequence markers (GASMs) in a restricted two-stage multi-locus model genome-wide association study (GASM-RTM-GWAS). A representative sample with 394 accessions was tested for their shade tolerance index (STI), in Nanning, China. Through whole-genome re-sequencing, 47,586 GASMs were assembled. From GASM-RTM-GWAS, 53 main-effect STI genes with 281 alleles (2-13 alleles/gene) (totally 63 genes with 308 alleles, including 38 G × E genes with 191 alleles) were identified and then organized into a gene-allele matrix composed of eight submatrices corresponding to geo-seasonal subpopulations. The population featured mild STI changes (1.69 → 1.56-1.82) and mild gene-allele changes (92.5% alleles inherited, 0% alleles excluded, 7.5% alleles emerged) from the primitive (SAIII) to the derived seven subpopulations, but large transgressive recombination potentials and optimal crosses were predicted. The 63 STI genes were annotated into six biological categories (metabolic process, catalytic activity, response to stresses, transcription and translation, signal transduction and transport and unknown functions), interacted as gene networks. From the STI gene-allele system, 38 important alleles of 22 genes were nominated for further in-depth study. GASM-RTM-GWAS performed powerful and efficient in germplasm population genetic study comparing to other procedures through facilitating direct and thorough identification of its gene-allele system, from which genome-wide breeding by design could be achieved, and evolutionary motivators and gene-allele networks could be explored.

An Improved Genome-Wide Association Procedure Explores Gene-Allele Constitutions and Evolutionary Drives of Growth Period Traits in the Global Soybean Germplasm Population.

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 ADLDSF and AATDSF, 130-141 genes with 384-406 alleles were explored, and for DFM and its related ADLDFM and AATDFM, 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.

Differential SW16.1 allelic effects and genetic backgrounds contributed to increased seed weight after soybean domestication.

Although seed weight has increased following domestication from wild soybean (Glycine soja) to cultivated soybean (Glycine max), the genetic basis underlying this change is unclear. Using mapping populations derived from chromosome segment substitution lines of wild soybean, we identified SW16.1 as the causative gene underlying a major quantitative trait locus controlling seed weight. SW16.1 encodes a nucleus-localized LIM domain-containing protein. Importantly, the GsSW16.1 allele from wild soybean accession N24852 had a negative effect on seed weight, whereas the GmSW16.1 allele from cultivar NN1138-2 had a positive effect. Gene expression network analysis, reverse-transcription quantitative polymerase chain reaction, and promoter-luciferase reporter transient expression assays suggested that SW16.1 regulates the transcription of MT4, a positive regulator of seed weight. The natural variations in SW16.1 and other known seed weight genes were analyzed in soybean germplasm. The SW16.1 polymorphism was associated with seed weight in 247 soybean accessions, showing much higher frequency of positive-effect alleles in cultivated soybean than in wild soybean. Interestingly, gene allele matrix analysis of the known seed weight genes revealed that G. max has lost 38.5% of the G. soja alleles and that most of the lost alleles had negative effects on seed weight. Our results suggest that eliminating negative alleles from G. soja led to a higher frequency of positive alleles and changed genetic backgrounds in G. max, which contributed to larger seeds in cultivated soybean after domestication from wild soybean. Our findings provide new insights regarding soybean domestication and should assist current soybean breeding programs.