The genetic variation and relationships among six Turkish water buffalo populations, typical of different regions, were assessed using a set of 26 heterologous (bovine) microsatellite markers. Between seven and 17 different alleles were identified per microsatellite in a total of 254 alleles. The average number of alleles across all loci in all the analysed populations was found to be 12.57. The expected mean heterozygosity (H(e)) per population ranged between 0.5 and 0.58. Significant departures from Hardy-Weinberg equilibrium were observed for 44 locus-population combinations. Population differentiation was analysed by estimation of the F(st) index (values ranging from 0.053 to 0.123) among populations. A principal component analysis of variation revealed the Merzifon population to show the highest differentiation compared with the others. In addition, some individuals of the Danamandira population appeared clearly separated, while the Afyon, Coskun, Pazar and Thural populations represented a single cluster. The assignment of individuals to their source populations, performed using the Bayesian clustering approach implemented in the structure 2.2 software, supports a high differentiation of Merzifon and Danamandira populations. The results of this study are useful for the development of conservation strategies for the Turkish buffalo.