Phenotypic and genetic diversity of domestic yak (Bos grunniens) in high-altitude rangelands of Gilgit-Baltistan, Pakistan

Due to climate change, irrigated agriculture may become restricted in the mountain valleys in northern Pakistan in the future. Hence, the importance of yak (Bos grunniens) keeping in the mountain ranges as risk-mitigating strategy for mountain dwellers will potentially increase. However, little is known about the current status of the domestic yak in this region. We therefore used phenotypic characteristics and 13 microsatellite loci to determine the phenotypic differences and the level of genetic differentiation between populations of six valleys. Larger body measures and partially different physical appearance were observed in Shimshal and Khaplu yaks, especially when compared with yaks in the Chapurson valley. Overall, the mean observed heterozygosity was similar to the mean expected heterozygosity. Average genetic diversity was highest in the Hopar population and lowest in the Haramosh population. A low FIS value indicated that individuals were less related than expected under a model of random mating. Three distinct genetic clusters were found for the six yak populations under study. Genetic distances were largest between Shimshal and Khaplu populations, and lowest between populations of Phandar and Hopar. It is concluded that yaks of Shimshal, Khaplu and Haramosh valleys were genetically distinct from yak populations in Chapurson, Hopar and Phandar valleys, indicating that the free-range conditions and pastoral yak rearing system in the region have preserved the underlying genetic diversity of the yak populations.