Montenegro has a great diversity of indigenous and locally developed sheep breeds, adapted to different regions and climates. However, the current trend of the steady decline in sheep populations means that some of them are threatened with extinction. The aim of this study was the investigation of the genetic diversity and population structure of five Montenegrin sheep populations, using the OvineSNP50K BeadChip. Data from the studied sheep were supplemented with publicly available genotypes of worldwide breeds to determine admixture contributions and genetic relationships. Higher genetic diversity and low inbreeding were observed in Pivska Pramenka and Sjenicka sheep, with Ho = 0.411, He = 0.423, and FROH>2Mb = 0.056 and Ho = 0.421, He = 0.421, and FROH>2Mb = 0.041, respectively. Zetska Zuja had much lower FIS (0.001) and high FROH inbreeding (0.090), indicating a small population size and genetic drift. An analysis of the genetic population structure showed that Montenegrin sheep populations are genetically distinct populations. This analysis revealed the presence of admixture in most sheep populations (high heterogeneity), while genetic purity was observed in Zetska Zuja. The obtained results will help in the conservation management of endangered Montenegrin sheep populations.
