The TASSLE v. 3.0 software was used to calculate the linkage disequilibrium (r 2) for all pairwise loci of the SNP markers . The General Linear Model (GLM) and the Q matrix from STRUCTURE software were used to identify the association of 172 SSR and 5,107 SNP markers with flower time . The Bonferroni-corrected thresholds for the p-value were used to determine the significance of association and were 2.90 × 10−4 (0.05/172), and 9.79 × 10−6 (0.05/5107) for SSR and SNP markers, respectively. Functional annotations of SNPs and SSRs were performed using the Phytozome database (https://phytozome.jgi.doe.gov) and SoyBase database (SoyBase (http://soybase.org).