Exploring genetic diversity and population structure of Punjab goat breeds using Illumina 50 K SNP bead chip

Pakistan has 35 goat breeds. Moreover, the province of Punjab has highest goat population constituting 37% of country's total population with seven goat breeds including Beetal, Daira Deen Panah, Nachi, Barbari, Teddi, Pahari, and Pothwari. The diversity study of breeds warrants the documentation of breeds particularly using genome wide panel of markers, i.e., SNP chip. The objective of the current study was to fill this gap of information. Therefore, in current study we collected total of 879 unrelated goat blood samples along with data on body weight measurements; genomic DNA was extracted, and genotyping was carried out using 50 K SNP bead chip. Quality control measures were performed in Plink 1.07. Genetic diversity was observed among studied populations using heterozygosity and pairwise F_ST estimates, principal component analysis, admixture analysis in Plink software with visualization in Clumpak, and constructing phylogenetic tree in Mega 7 software. Moderate to high level of heterozygosity was observed among the studied populations. Coefficient of inbreeding varied from 0.0186 ± 0.0327 in Pahari to 0.183 ± 0.0715 in Barbari. Barbari and Daira Deen Panah had quite higher level of inbreeding coefficient as compared to all other breeds with value of 0.183 ± 0.0715 and 0.1378 ± 0.0741, respectively. PCA identified three steps of subdividing the seven goat breeds at various levels of K. All the seven breeds made independent clusters at various levels of PCA. Admixture analysis revealed the distinctness of Teddi and Barbari breeds. Genetic sub-structuring was observed in the admixture patterns of Beetal breed. Moreover, high level of genetic admixture was observed in Nachi, Pahari, Pothwari, and Daira Deen Panah breeds. Admixture results were further interpreted by calculating pairwise F_ST values. Our results provided first insights about genetic diversity of Pakistani goat breeds based on genomic data. To conclude, the enriched goat breed diversity in Pakistan could provide valuable genetic reservoir for national breeding schemes.