Genetic Diversity, Population Structure, and Historical Gene Flow Patterns of Nine Indigenous Greek Sheep Breeds

Ιn this study, we evaluated the genetic resources of nine Greek sheep breeds. The genotyping data of 292 animals were acquired from Illumina's OvineSNP50 Genotyping BeadChip. The genetic diversity and inbreeding levels were evaluated using the observed and expected heterozygosity indices, the F_IS inbreeding coefficient, and runs of homozygosity (ROH). The genetic differentiation of breeds was assessed using the F_ST index, whereas their population structure was analyzed using admixture and principal components analysis (PCA). Historical recombination patterns and genetic drift were evaluated based on linkage disequilibrium, effective population sizes, and gene flow analysis to reveal migration patterns. PCA revealed distinct clusters mostly separating mountainous, insular, and lowland breeds. The F_ST value was the lowest between Serres and Karagouniko breeds (0.050). Admixture analysis revealed a genetic substructure for Serres and Kalarritiko breeds, while Chios, followed by Katsika, demonstrated the highest within-breed genetic uniformity. ROH analysis revealed low levels of inbreeding for all breeds. Genetic introgression from both Anatolia and Eastern Europe has been evidenced for Greek sheep breeds. The results also revealed that Greek sheep breeds maintain adequate levels of genetic diversity, without signs of excessive inbreeding, and can serve as valuable resources for the conservation of local biodiversity.