Diversity analysis of sheep breeds from Southern peninsular and Eastern regions of India

Information is presented on the genetic diversity and relationship among six Indian sheep breeds/populations belonging to the Southern peninsular and Eastern agroecological zones, based on microsatellite markers. Parameters of genetic variation, viz., allele diversity, observed heterozygosity, gene diversity and population inbreeding estimates, were calculated for the six breeds. The allele diversity ranged from 6.40 to 7.92, whereas the gene diversity varied from 0.617 to 0.727. The highest allele and gene diversity was observed for Nellore sheep, while the lowest was exhibited by Garole breed. Within population inbreeding estimate (F _IS) revealed a significant deficit of heterozygotes in Deccani, Madgyal, Nellore and Garole, whereas Ganjam and Chhotanagpuri sheep showed an excess of heterozygotes. The contribution of each breed to the total diversity of the breeds was quantified by the Weitzman approach. The marginal loss of diversity incurred with removal of Nellore and Garole breeds was higher (>27%), whereas removal of Deccani breed resulted in lowest loss of diversity (3.84%) from the set. Estimation of the genetic differentiation (F _ST) and genetic distance (D _A) between the pairs of breeds revealed a close relationship between Deccani and Madgyal sheep (F _ST = 0.017; D _A = 0.080) and greatest demarcation between Madgyal and Garole breeds (F _ST = 0.110; D _A = 0.622). The information generated would help in shaping genetic management and conservation programs for the sheep breeds under consideration.     

Genetic diversity of Indian native cattle breeds as analysed with 20 microsatellite loci

The aim of the present study was to assess the genetic variation and establish the relationship amongst the three Indian zebu cattle breeds using 20 bovine‐specific microsatellite markers. A total of 136 unrelated DNA samples from Sahiwal (SC), Hariana (HC) and Deoni (DC) breeds of cattle were genotyped to estimate within and between breed genetic diversity indices. The estimated mean allelic diversity was 5.2, 6.5 and 5.9 in SC, HC and DC, respectively, with a total of 167 alleles. The average observed and expected heterozygosity for the population varied from 0.42 (SC) to 0.59 (DC), and from 0.61 (SC) to 0.70 (DC), respectively. Low values of genetic variability estimates were observed in SC when compared with DC and HC, indicating some loss of variability because of its relatively small population size. From global F‐statistics a significant deficit of heterozygotes of 24.2% (p < 0.05) was observed for each one of the analysed breeds whereas the total population had a 32.8% (p < 0.05) deficit of heterozygotes. The F_ST estimates demonstrated that approximately 88.7% of the total genetic variation was because of the genetic differentiation within each breed. Pair‐wise breed differentiation, Nei's standard and D_A genetic distance estimates revealed relatively close genetic similarity between HC and DC in comparison with SC. In the UPGMA‐based phylogenetic tree constructed from the genetic distances, HC and DC were grouped together in one cluster and SC in the other. The estimated time of divergence suggested a separation time of approximately 776 years between DC and HC, and a comparatively longer period (1296 years) between DC and SC.     

Molecular characterization and differentiation of five horse breeds raised in Algeria using polymorphic microsatellite markers

In this study, genetic analyses of diversity and differentiation were performed on five horse breeds raised in Algeria (Barb, Arab‐Barb, Arabian, Thoroughbred and French Trotter). All microsatellite markers were highly polymorphic in all the breeds. A total of 123 alleles from 14 microsatellite loci were detected in 201 horses. The average number of alleles per locus was the highest in the Arab‐Barb horses (7.86) and lowest in the thoroughbred breed (5.71), whereas the observed and expected heterozygosities per breed ranged from 0.71 (Thoroughbred) to 0.752 (Barb) and 0.71 (Thoroughbred) to 0.77 (Arab‐Barb), respectively. The genetic differentiation between the breeds was significant (p < 0.01) based on the infinitesimal model (F_ST). Three different approaches for evaluating the genetic relationships were applied. Genetic distances, the factorial correspondence analysis and structure analysis showed that a significant amount of genetic variation is maintained in the native horse populations and the other breeds. The Barb and Arab‐Barb breeds seem to be the most genetically related and support the decision to consider the breeds as same population.     

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein‐coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.     

Genetic diversity and differentiation of five Cuban cattle breeds using 30 microsatellite loci

Conservation and improvement strategies in farm animals should be based on a combination of genetic and phenotypic characteristics. Genotype data from 30 microsatellites were used to assess the genetic diversity and relationships among five Cuban cattle breeds (Siboney de Cuba, Criollo Cubano, Cebú Cubano, Mambí de Cuba and Taíno de Cuba). All microsatellite markers were highly polymorphic in all the breeds. The expected heterozygosity ranged from 0.67 ± 0.02 in the Taíno de Cuba breed to 0.75 ± 0.02 in the Mambí de Cuba breed, and the observed heterozygosity ranged from 0.66 ± 0.03 in the Cebú Cubano breed to 0.73 ± 0.02 in the Siboney de Cuba breed. The genetic differentiation between the breeds was significant (p < 0.01) based on the infinitesimal model (F_ST). The exact test for Hardy‐Weinberg equilibrium within breeds showed a significant deviation in each breed (p < 0.0003) for one or more loci. The genetic distance and structure analysis showed that a significant amount of genetic variation is maintained in the local cattle population and that all breeds studied could be considered genetically distinct. The Siboney de Cuba and Mambí de Cuba breeds seem to be the most genetically related among the studied five breeds.     

Evaluation of genetic diversity in Japanese Brown cattle population by pedigree analysis

The Japanese Brown is the second most common domestic beef breed in Japan. However, nowadays this breed is facing reduction in numbers because of pressure from a profitable domestic breed. This breed is uniformly characterized by its brown coat colour, but is comprised of two isolated sub‐breeds, Kumamoto and Kouchi, each possessing a different gene pool. Pedigree analyses were carried out for the two sub‐breeds using the pedigree records of animals born from 1970 to 2000. The effective population size has been found to be consistently reducing during the last three decades in both sub‐breeds. The current effective sizes were estimated to be 25.5 and 6.0 for the Kumamoto and Kouchi sub‐breeds, respectively. The estimate of the effective number of founders (N_ef) in the Kumamoto sub‐breed decreased from 152.1 to 74.4; that of non‐founders (N_enf), from 41.7 to 5.3; and that of founder genome equivalents (N_ge), from 32.7 to 4.9. The corresponding changes in the Kouchi sub‐breed were from 108.2 to 79.4, 16.2 to 4.1, and 14.1 to 3.9. Increasing differences between the two genetic diversity indices in the sub‐breeds indicate that the greater part of the reduction of genetic diversity can be attributed to genetic drift that accumulated in the non‐founder generations. A comparison with published estimates for several cattle breeds suggests the extremely limited genetic diversity of Japanese Brown. In addition to the avoidance of further reduction of genetic diversity, it will be important to counteract the process of breed decline by establishing a production system to efficiently utilize the unique characteristics of this breed and by developing links between the breed and products with market value.     

Genetic diversity,population structure,and correlations between locally adapted zebu and taurine breeds in Brazil using SNP markers

Genetic diversity is one of the most important issues in studies on conservation of cattle breeds and endangered species. The objective of this study was to estimate the levels of genetic differentiation between locally adapted taurine (Bos taurus taurus) and zebu (Bos taurus indicus) breeds in Brazil, which were genotyped for more than 777,000 SNPs. The fixation index (F _ST), principal component analysis (PCA), and Bayesian clustering were estimated. The F _ST highlighted genetic differentiation between taurine and zebu breeds. The taurine lines, Caracu and Caracu Caldeano, had significant genetic differentiation (F _ST close to 5%) despite their recent selection for different uses (meat and milk). This genetic variability can be used for conservation of locally adapted animals, as well as for breeding programs on zebu breeds. Introgression of zebu in locally adapted breeds was identified, especially in Curraleiro Pé-Duro breed. The Gyr breed, however, had low breed purity at genomic level due to its very heterogeneous mixing pattern.     

Genetic diversity measures of the bovine Alberes breed using microsatellites: variability among herds and types of coat colour*

The Alberes population is a native bovine breed of Catalonia with an unclear origin, which historically some authors have assumed as being composed of two different colour varieties (black and fawn). Sixteen microsatellite loci were analysed, all of them included in the AIRE2066 European Concerted Action list. Overall expected and observed heterozygosities reached values of 0.649 and 0.662, respectively. Genetic differences among black and fawn varieties were not significant (F_ST = 0.007), and therefore the population is a single variety with a great colour gradation. On the contrary, we detected significant genetic differences among herds (F_ST = 0.026; p < 0.001), showing a genetic heterogeneity over short geographical distances. The number of migrants per generation among pairs of herds oscillates between 1.46 (Roig and Freixe herds) and 5.62 (Castanyers and Roig herds). Moreover, inbreeding and bottleneck situations can be rejected. The Alberes breed has been grouped within the Cantabrian trunk, closely related to the Asturiana de la Montaña and Alistana breeds, although some other breeds may also have influenced the population along its history.     

Combined approaches to identify genomic regions involved in phenotypic differentiation between low divergent breeds: Application in Sardinian sheep populations

Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case–control genome‐wide association study (GWAS), the Bayesian population differentiation analysis (F_ST) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and F_ST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (<2 Mb) to highly significant regions identified by GWAS and F_ST analyses. Moreover, one candidate marker identified by GWAS and F_ST approaches was located in a run of homozygosity island that was shared by both breeds. We identified several genes involved in many phenotypic differences (such as stature and growth, reproduction, ear size, coat colour, behaviour) between the two analysed breeds. This study shows that combining several genome‐wide approaches could improve discovery of regions involved in the variability of breeding traits and responsible for the phenotypic diversity even between closely related breeds. Overall, the combination of such genome‐wide methods can be extended to other livestock breeds that share between them a similar genetic background, to understand the process that shapes the patterns of genetic variability between closely related populations.     

Adult body weight of indigenous and crossbred ewes under semi-arid conditions in India

Data on adult body weight of ewes belonging to Nali, Chokla breeds and their crosses (F1) with Rambouillet and Soviet Merino raised under semi-arid conditions in India were analysed. The ewes belonged to 5 distinct age groups ranging from 16.5 to 64 months. Nali and Chokla were significant lighter than all the 4 crossbred groups. There were no significant differences in body weight between crossbred groups except in the case of Soviet Merino x Chokla group which was the lightest and differed significantly from all other groups. Rambouillet crosses were heavier than Merino crosses but the difference was statistically significant only on Rambouillet x Chokla and Merino x Chokla. The improvement in body weight in crossbreds over native sheep ranged between 9.3 and 21.88%. The increase in live weight was significant only up to 49 months.     

Sperm Motion Characteristics of Garole Rams Raised for a Prolonged Period in a Semi-arid Tropical Environment

The Garole is a prolific but less well known and rare breed of small sheep found in the hot and humid Sunderban region of West Bengal. An ability to breed throughout the year and to graze in knee-deep water, resistance to foot rot and a strong mothering instinct are some of the special features of this breed. Garole rams could provide germplasm to incorporate prolificacy traits by artificial insemination of the nonprolific sheep breeds found in abundance in the semi-arid and arid tropical climates of India. The aim of the present study was to evaluate the semen production by Garole rams maintained in a semi-arid climate for three years and to objectively assess their semen quality by a computer-assisted sperm analysis technique. The donor rams were randomly selected each year from the original flock procured from their natural habitat or from the offspring born at the Institute farm. Semen was collected weekly for three weeks each autumn for three consecutive years from 8 rams each year. The overall means (SD) of the traits that did not differ significantly with age or year were volume, concentration, curvilinear velocity, average path velocity, amplitude of lateral head displacement, beat frequency, motility and the percentages of rapid motile sperms and of slow motile sperms. The age of the rams had a significant effect (p<0.05) on the straight-line velocity but this was not significantly affected by the length of exposure to the semi-arid climate. However, the age and year had significant effects (p<0.05) on linearity, straightness and the percentage of medium motile sperms. It was concluded that Garole rams are capable of producing good-quality semen even after a prolonged period of exposure to a semi-arid tropical climate.     

Genetic structure and differentiation of 12 African Bos indicus and Bos taurus cattle breeds, inferred from protein and microsatellite polymorphisms

Level of genetic differentiation, gene flow and genetic structuring of nine Bos indicus and three Bos taurus cattle breeds in Cameroon and Nigeria were estimated using the genetic information from 16 microsatellite, five blood protein and seven milk protein markers. The global heterozygote deficit across all populations (F_it) amounted to 11.7% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds (F_is) amounted to 6.1%. The breeds were moderately differentiated (F_st = 6%, p < 0.001) with all loci except CSN1S2 contributing significantly to the F_st value. The 12 populations belong to two genetic clusters, a zebu and a taurine cluster. While inferred sub‐clusters within the taurine group corresponded extremely well to predefined breed categorizations, no real sub‐clusters, corresponding to predefined breeds, existed within the zebu cluster. With the application of prior population information, cluster analysis achieved posterior probabilities from 0.962 to 0.994 of correctly assigning individuals to their rightful populations. High gene flow was evident between the zebu populations. Positive and negative implications of the observed genetic structure of the breeds on their development, improvement and conservation are discussed. The study shows that the breeds are threatened by uncontrolled breeding and therefore are at risk to become genetically uniform in the future. This situation can be avoided by putting in place effective breeding and management measures aimed at limiting uncontrolled mating between the breeds and to preserve special characteristics, genetic as well as breed biodiversity. The first step towards realizing these goals might be to geographically demarcate the breeds.     

Population structure and genetic differentiation of livestock guard dog breeds from the Western Balkans

Livestock guard dog (LGD) breeds from the Western Balkans are a good example of how complex genetic diversity pattern observed in dog breeds has been shaped by transition in dog breeding practices. Despite their common geographical origin and relatively recent formal recognition as separate breeds, the Karst Shepherd, Sarplaninac and Tornjak show distinct population dynamics, assessed by pedigree, microsatellite and mtDNA data. We genotyped 493 dogs belonging to five dog breeds using a set of 18 microsatellite markers and sequenced mtDNA from 94 dogs from these breeds. Different demographic histories of the Karst Shepherd and Tornjak breeds are reflected in the pedigree data with the former breed having more unbalanced contributions of major ancestors and a realized effective population size of less than 20 animals. The highest allelic richness was found in Sarplaninac (5.94), followed by Tornjak (5.72), whereas Karst Shepherd dogs exhibited the lowest allelic richness (3.33). Similarly, the highest mtDNA haplotype diversity was found in Sarplaninac, followed by Tornjak and Karst Shepherd, where only one haplotype was found. Based on F_ST differentiation values and high percentages of animals correctly assigned, all breeds can be considered genetically distinct. However, using microsatellite data, common ancestry between the Karst Shepherd and Sarplaninac could not be reconstructed, despite pedigree and mtDNA evidence of their historical admixture. Using neighbour‐joining, STRUCTURE or DAPC methods, Sarplaninac and Caucasian Shepherd breeds could not be separated and additionally showed close proximity in the NeighborNet tree. STRUCTURE analysis of the Tornjak breed demonstrated substructuring, which needs further investigation. Altogether, results of this study show that the official separation of these dog breeds strongly affected the resolution of genetic differentiation and thus suggest that the relationships between breeds are not only determined by breed relatedness, but in small populations even more importantly by stochastic effects.     

Aggregate diversity: New approach combining within- and between-breed genetic diversity

Between-breed genetic diversity is classically considered as a major criterion to be taken into account when setting priorities for conservation of domestic animal breeds. However, it has been argued that methods based on the between-breed component of genetic diversity may not be optimal because they ignore the within-breed component of variation. The paper considers the most common methods used to evaluate those two components when genetic diversity is evaluated on the basis of genetic markers, and proposes to define an aggregate diversity combining linearly the two components. This implies defining for each breed (or population) its contributions to the between-breed and to the within-breed diversity. When defining an aggregate diversity, one can weight these contributions by F_ST and 1−F_ST, respectively, since the fixation index F_ST of Wright represents the proportion of the total genetic variation which is due to differences in allelic frequencies between populations. Such an approach is valid when the objective is genetic improvement by selection within a so-called “meta-population”. However, in a more general context of animal breeding, when heterosis and complementarities between breeds have to be considered, as well as adaptation to specific environments, more weight should be given to the between-breed variation. The proper weight to apply may require solutions adapted to each particular situation. In a long-term conservation perspective, priorities should also take into account the degree of endangerment of each breed. By combining diversity contributions and probability of extinction, a cryopreservation potential (or priority) may be estimated for each breed. The problem is illustrated on a sample of 11 European pig breeds typed for 18 microsatellite loci.     

An insight into population structure and gene flow within pure‐bred cats

Investigation of genetic structure on the basis of pedigree information requires indicators adapted to the specific context of the populations studied. On the basis of pedigree‐based estimates of diversity, we analysed genetic diversity, mating practices and gene flow among eight cat populations raised in France, five of them being single breeds and three consisting of breed groups with varieties that may interbreed. When computed on the basis of coancestry rate, effective population sizes ranged from 127 to 1406, while the contribution of founders from other breeds ranged from 0.7 to 16.4%. In the five breeds, F_IS ranged between 0.96 and 1.83%, with this result being related to mating practices such as close inbreeding (on average 5% of individuals being inbred within two generations). Within the three groups of varieties studied, F_IT ranged from 1.59 to 3%, while values were estimated between 0.04 and 0.91%, which was linked to various amounts of gene exchanges between subpopulations at the parental level. The results indicate that cat breeds constitute populations submitted to low selection intensity, contrasting with relatively high individual inbreeding level caused by close inbreeding practices.     

The assessment of genetic diversity within and among the eight subpopulations of Japanese Black cattle using 52 microsatellite markers

Japanese Black cattle are at risk for genetic homogeneity due to intensive use of a few sires. Therefore, assessment of the actual genetic diversity of this breed is important for future breeding plans. In the present study, we investigated the genetic diversity within and among eight subpopulations of Japanese Black cattle using 52 microsatellite markers. The parameters for genetic diversity of Japanese Black cattle were comparable to those of other cattle breeds, suggesting that the relatively high genetic diversity of the breed. However, upon comparison among the eight subpopulations, the Hyogo subpopulation showed markedly low genetic diversity. The results of the pairwise F_ST values, phylogenetic network and structure analysis indicated that the Hyogo population has remarkably high level of genetic differentiation from other populations, while Yamagata, Niigata, Hiroshima and Kagawa populations have low levels of genetic differentiation. Furthermore, multidimensional scaling plots indicated that individuals in some subpopulations were separated from individuals in the other subpopulations. We conclude that while the overall genetic diversity of Japanese Black cattle is still maintained at a relatively high level, that of a particular subpopulation is significantly reduced, and therefore the effective population size of the breed needs to be controlled by correct mating strategies.     

Pedigree information reveals moderate to high levels of inbreeding and a weak population structure in the endangered Catalonian donkey breed

The Catalonian donkey is one of the most endangered donkey breeds in the world. At present, five main subpopulations exist: AFRAC, which consists of many genetically connected Catalonian localities; Berga, which consists of a single herd located also in Catalunya but under private management; and three minor non‐Catalonian subpopulations (Huesca, Sevilla and Toledo). In this study, we analysed the pedigree information of the Catalonian donkey herdbook to assess the genetic diversity and population structure of the breed. We found that the Catalonian donkey has suffered an important loss of genetic diversity and moderate to high increases of inbreeding because of the abuse of a few individuals in matings. This scenario is mainly characterized by the fact that both the effective number of founders and ancestors for the whole population was 70.6 and 27, respectively, while the equivalent number of founders was 146.5 and the number of ancestors explaining overall genetic variability was 93. In addition, only 14% of animals born between the 1960s and 1970s were significantly represented in the pedigree. Our results also show that subpopulations where breeders exchanged reproductive individuals had low levels of inbreeding and average relatedness. One subpopulation, Berga, was reproductively isolated and showed high levels of inbreeding (F = 7.22%), with average relatedness (AR = 6.61%) playing an important role in increasing the values of these coefficients in the whole pedigree. Using genealogical F‐statistics we have found little evidence of population structuring (F_ST = 0.0083) with major genetic differences among non‐Catalonian subpopulations.     

Genetic population structure of red snapper (Lutjanus malabaricus) and orange‐spotted grouper (Epinephelus coiodes) in Brunei and Sabah

The genetic population structure of red snapper Lutjanus malabaricus and orange‐spotted grouper Epinephelus coiodes in Brunei and Sabah was investigated using allozyme electrophoresis. Samples were collected from three sites in Brunei for E. coiodes and from three sites in Brunei and Sabah for L. malabaricus. A total of 22 loci and 16 loci were scored, respectively. The index of fixation (F_ST) for the E. coiodes population was 0.176 but, in general, it lacked within‐population structuring. The F_ST was particularly high between Brunei Bay/Jerudong and Brunei Bay/ Kuala Belait, suggesting genetic subdivision on a small spatial scale. Isolation of Brunei Bay from the South China Sea may constrain the movement of adult fishes and larval dispersal, thereby reducing homogeneity among subpopulations. All variable loci for E. coiodes were in Hardy‐Weinberg equilibrium except for MDH* and GTDHP* (P < 0.01), in which two subpopulations showed an excess of heterozygotes (P < 0.01). The study on the L. malabaricus population showed a heterozygote deficit of approximately 60% in variable loci (F_ST genetic variation within population = 0.45; P < 0.05); however, the mean observed heterozygosity for the population far exceeded L. malabaricus populations in Australia and Indonesia. A F_ST value of 0.076 revealed moderate genetic differentiation among subpopulations of L. malabaricus. The genotypes were likely to be drawn from the same distribution in Jerudong and Kuala Belait. This study infers that sustainable management of snapper and grouper resources in Brunei waters must take into account the presence of a single stock and two stocks, respectively.     

Population structure of pigs determined by single nucleotide polymorphisms observed in assembled expressed sequence tags

We have collected more than 190 000 porcine expressed sequence tags (ESTs) from full‐length complementary DNA (cDNA) libraries and identified more than 2800 single nucleotide polymorphisms (SNPs). In this study, we tentatively chose 222 SNPs observed in assembled ESTs to study pigs of different breeds; 104 were selected by comparing the cDNA sequences of a Meishan pig and samples of three‐way cross pigs (Landrace, Large White, and Duroc: LWD), and 118 were selected from LWD samples. To evaluate the genetic variation between the chosen SNPs from pig breeds, we determined the genotypes for 192 pig samples (11 pig groups) from our DNA reference panel with matrix‐assisted laser desorption ionization time‐of‐flight mass spectrometry. Of the 222 reference SNPs, 186 were successfully genotyped. A neighbor‐joining tree showed that the pig groups were classified into two large clusters, namely, Euro‐American and East Asian pig populations. F‐statistics and the analysis of molecular variance of Euro‐American pig groups revealed that approximately 25% of the genetic variations occurred because of intergroup differences. As the F_IS values were less than the F_ST values_, the clustering, based on the Bayesian inference, implied that there was strong genetic differentiation among pig groups and less divergence within the groups in our samples.     

The phylogeographic system survey of native sheep breeds in the eastern and southern Central Asia

The genetic diversity and phylogenetic survey of native sheep breeds in the eastern and southern Central Asia were assessed in the present study. The clustering, principal components, structure and F statistics all demonstrate that the native sheep breeds in these regions be classified into two genetic groups: Mongolia‐Tibetan sheep group and South‐Southeast Asia sheep group. The Mongolia sheep group and the Tibetan sheep group had a certain degree of gene communication from the ancient times. In the present study we demonstrated that the Chinese native sheep populations belonged to Mongolia‐Tibetan sheep group. However, the relationships among the sheep populations in Mongolia sheep group in China were not closely related to the geographical distance among sheep populations.