### INTRODUCTION

### MATERIALS AND METHODS

_{IS}), genetic distance (F

_{ST}) indices were calculated by Genalex 6.5 [18]. Data were crosschecked by FSTAT and Genepop.

^{5}, MCM steps: 5×10

^{5}, repetition: 5, model: admixture, allele frequencies correlated) and for calculation of membership probability of individuals. For estimation of K, Evanno’s method [19] was applied on Structure output. Bayesian stochastic partition-based approach implemented in BAPS 6.0 [20] was also applied to estimate K.

### RESULTS AND DISCUSSION

_{ST}values (p<0.000, Table 3.) revealed two herds (C and E) distinct from the majority of other Limousin herds (Figure 2). The pairwise F

_{ST}values of population C compared to the others (0.066 to 0.120) fell into the range of moderate genetic distance: 0.050 to 0.150 [25], while population E displayed also moderate genetic distance (F

_{ST}values in range 0.052 to 0.064) but only to six populations (G, H, J, L, N, and P). F

_{ST(C-E)}was 0.148, all other pairs—excluding C and E herds—displayed low genetic distance (F

_{ST}<0.049).

_{Evanno}= 2, Figure 3). Structure indicated only two major groups, where population A, C, N were separated from the remaining 13 populations.

_{ST}values. Population B and K remained together on the dendrogram. The remaining groups are the least divergent from each other. Plotting neighbour joining tree of genetic distance (data not shown) the A, C, N populations remained separated. E herd was more similar to the remaining herds and it shared a node with population F -as we see on the BAPS generated result (Figure 4) but was placed on the longest branch among populations (excluding A, C, N) in accordance with the PCA analysis.