Genome-wide association analysis identified consistent QTL for seed yield in a soybean diversity panel tested across multiple environments.

Improving seed yield is one of the main targets of soybean [Glycine max (L.) Merr.] breeding. Identification of loci that influence productivity and understanding their genetic mechanism will help marker-assisted trait introgression. The present study evaluated a diverse panel of 541 soybean genotypes consisting of three maturity groups (MGs III-V) in four environments in Kansas, U.S. Data on seed yield, seed weight, shattering resistance, days to maturity, and plant height showed significant genotype, environmental, and genotype × environment interaction variations. Seed yield and shattering had moderate broad-sense heritability (<85%), while the rest of the traits showed high broad-sense heritability (>90%). The SoySNP50K iSelect BeadChip dataset was used to identify significantly associated loci via genome-wide association studies (GWAS). A total of 19 single-nucleotide polymorphisms (SNPs) were significantly associated with seed yield. Particularly, two stable seed yield quantitative trait loci (QTL) on chromosomes 9 and 17 were consistently detected in at least three out of four environments. Candidate gene analysis surrounding seed yield QTL on chromosome 9 showed that Glyma.09G048900, an oxygen binding protein, was the closest to the QTL peak. Similarly, Glyma.17G090200 and Glyma.17G090400 were within 20-kb region of the seed yield QTL on chromosome 17. The candidate genes warrant further analysis to determine their functional mechanisms and develop markers for seed yield improvement.