Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 May 31;18(5):e0286099.
doi: 10.1371/journal.pone.0286099. eCollection 2023.

Genetic diversity and population structure analysis in cultivated soybean (Glycine max [L.] Merr.) using SSR and EST-SSR markers

Reena Rani  1   2 , Ghulam Raza  1   2 , Muhammad Haseeb Tung  1   2 , Muhammad Rizwan  3 , Hamza Ashfaq  1   2 , Hussein Shimelis  4 , Muhammad Khuram Razzaq  5 , Muhammad Arif  1   2
Affiliations

Genetic diversity and population structure analysis in cultivated soybean (Glycine max [L.] Merr.) using SSR and EST-SSR markers

Reena Rani et al. PLoS One. .

Abstract

Soybean (Glycine max) is an important legume that is used to fulfill the need of protein and oil of large number of population across the world. There are large numbers of soybean germplasm present in the USDA germplasm resources. Finding and understanding genetically diverse germplasm is a top priority for crop improvement programs. The current study used 20 functional EST-SSR and 80 SSR markers to characterize 96 soybean accessions from diverse geographic backgrounds. Ninety-six of the 100 markers were polymorphic, with 262 alleles (average 2.79 per locus). The molecular markers had an average polymorphic information content (PIC) value of 0.44, with 28 markers ≥ 0.50. The average major allele frequency was 0.57. The observed heterozygosity of the population ranged from 0-0.184 (average 0.02), while the expected heterozygosity ranged from 0.20-0.73 (average 0.51). The lower value for observed heterozygosity than expected heterozygosity suggests the likelihood of a population structure among the germplasm. The phylogenetic analysis and principal coordinate analysis (PCoA) divided the total population into two major groups (G1 and G2), with G1 comprising most of the USA lines and the Australian and Brazilian lines. Furthermore, the phylogenetic analysis and PCoA divided the USA lines into three major clusters without any specific differentiation, supported by the model-based STRUCTURE analysis. Analysis of molecular variance (AMOVA) showed 94% variation among individuals in the total population, with 2% among the populations. For the USA lines, 93% of the variation occurred among individuals, with only 2% among lines from different US states. Pairwise population distance indicated more similarity between the lines from continental America and Australia (189.371) than Asia (199.518). Overall, the 96 soybean lines had a high degree of genetic diversity.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. World map representing geographic regions of the 96 accessions used in the study.
Created using map chart (mapchart.com).
Fig 2
Fig 2. Phylogenetic tree of 96 soybean accessions using the UPGMA method.
Accessions in group 1 (G1) are represented with red color while accessions in group 2 (G2) are represented with green color.
Fig 3
Fig 3. Phylogenetic tree of 59 USA accessions using the UPGMA method.
Fig 4
Fig 4. Principal coordinate analysis of 96 soybean accessions using 96 markers to identify variation among accessions based on their country of origin: USA (red), China (brown), Brazil (light green), Pakistan (pink), India (dark green), Iran (black), Afghanistan (blue), Australia (violet), and Japan (gray).
Fig 5
Fig 5. Principal coordinate analysis showing variation between soybean accessions from different states in the USA.
Fig 6
Fig 6
(A) Graph of estimated membership fraction for K = 2. (B) Graphical representation of 96 soybean accessions using 96 markers for K = 2.
Fig 7
Fig 7
(A) Graph of estimated K value for 59 USA soybean accessions with K = 9. (B) Population structure of 59 USA soybean accessions using SSR and EST-SSR markers.
See this image and copyright information in PMC

References

    1. Pembele Ibanda A, Karungi J, Malinga GM, Adjumati Tanzito G, Ocan D, Badji A, et al. Influence of environment on soybean [Glycine max (L.) Merr.] resistance to groundnut leaf miner, Aproaerema modicella (Deventer) in Uganda. Journal of Plant Breeding and Crop Science. 2018;10(12):336–46. 10.5897/JPBCS2018.0764. - DOI
    1. Kumar SJ, Kumar A, Ramesh K, Singh C, Agarwal DK, Pal G, et al. Wall bound phenolics and total antioxidants in stored seeds of soybean (Glycine max) genotypes. Indian Journal of Agricultural Sciences. 2020;90:118–222.
    1. Hartman GL, Hill CB. 13 Diseases of Soybean and Their Management. The soybean: botany, production and uses: CABI Publishing Cambridge, USA; 2010.
    1. Singh G. The soybean: botany, production and uses: CABI; 2010.
    1. Gupta S, Manjaya J. Genetic diversity and population structure of Indian soybean [Glycine max (L.) Merr.] revealed by simple sequence repeat markers. Journal of Crop Science and Biotechnology. 2017;20(3):221–31. 10.1007/s12892-017-0023-0. - DOI