The present study aimed at analyzing the genetic variability of indigenous goat breeds (Capra hircus) using the MHC-associated microsatellite markers BF1, BM1818, BM1258, DYMS1, and SMHCC1. The following breeds were included: Chinese Xuhuai, Indian Changthangi and Pashmina, Kenyan Small East African (SEA) and Galla, and Albanian Vendi. To examine genetic variability, the levels of heterozigosity, degrees of inbreeding, and genetic differences among the breeds were analyzed. The mean number of alleles ranged from nine in the Galla to 14.5 in the Vendi breed. The mean observed heterozygosity and mean expected heterozygosity varied from 0.483 in the Vendi to 0.577 in the Galla breed, and from 0.767 in the SEA to 0.879 in the Vendi breed, respectively. Significant loss of heterozygosity (p < 0.01) indicated that these loci were not in Hardy-Weinberg equilibrium. The mean F(IS) values ranged from 0.3299 in the SEA to 0.4605 in the Vendi breed with a mean value of 0.3623 in all breeds (p < 0.001). Analysis of molecular variance indicated that 7.14% and 4.74% genetic variation existed among the different breeds and geographic groups, whereas 92.86% and 95.26% existed in the breeds and the geographic groups, respectively (p < 0.001). The microsatellite marker analysis disclosed a high degree of genetic polymorphism. Loss of heterozygosity could be due to genetic drift and endogamy. The genetic variation among populations and geographic groups does not indicate a correlation of genetic differences with geographic distance.