Genetic diversity of three Chinese native sheep breeds

The sheep (Ovis aries L.) has been an important farm animal species since its domestication. A wide array of indigenous sheep breeds with abundant phenotypic diversity exists for domestication and selection. Therefore, assessing the genetic diversity of a local sheep resource using a multi-molecular system is helpful for maintaining and conserving those breeds. This study aimed to investigate the genetic diversity of three native Chinese sheep breeds (Tibetan sheep, Sishui Fur sheep, and Small-tailed Han sheep) using 15 microsatellite markers and the second exon of the DRA gene. In regards to the microsatellites, on average, 19 alleles per loci were observed among all individuals. Across loci, the HO within the population was 0.652 ± 0.022 in Tibetan sheep, 0.603 ± 0.023 in Small-tailed Han sheep and 0.635 ± 0.022 in SFS, and for most populations, the H _E and H _O were inconsistent. In addition, affluent private alleles within the breed indicated that the breeds have different domestication histories or sites. In regards to the 2 exon of the DRA gene, three haplotypes were constructed by seven single-nucleotide polymorphisms (SNPs), which were identified in the second DRA exon and inferred the potential for phenotypic variety in these Chinese native sheep. In summary, the current study reveals the importance of implementing effective conservation strategies for these three native Chinese sheep.     

Nucleotide variants in prion-related protein (testis-specific) gene (PRNT) and effects on Chinese and Mongolian sheep phenotypes

ABSTRACT

Studies of the ovine prion-related protein (testis-specific) gene (PRNT), including studies of genetic diversity, have highlighted its potential relationship to scrapie infection and economically important ovine traits. PRNT was previously reported to be highly polymorphic in Portuguese sheep. To characterize genetic polymorphisms in this gene in Asian sheep, a direct sequencing method was used to detect polymorphic loci in PRNT in 285 individual sheep from four Chinese and one Mongolian breeds. Seven SNP variants in PRNT were identified, including three novel variants (g.93G>A, g.162G>T, and g.190A>G) and four previously reported variants (g.17 C>T, g.112G>C, g.129C>T, and g.144A>G). In the five breeds that we analyzed, the mutation frequencies of g.190A>G in Lanzhou Fat-tail sheep (LFTS) and g.129C>T in the other four varieties were high (F>0.5). Moreover, thirteen different haplotypes that had a comparable distribution in the tested breeds were also identified; ‘C-G-G-C-A-G-A’ occurred at the highest frequency in the five sheep breeds. Additionally, we previously explored the significance of relationships between polymorphisms in PRNP or PRND and ovine growth performance. Here, we also performed correlation analysis in all tested loci. These loci polymorphisms were significantly associated with ten different growth traits (P<0.05), except for g.93G>A. Meanwhile, in contrast to a previous study, there was no significant association between the seven SNP loci analyzed and our previously reported sheep PRND or PRNP insertion/deletion mutations. Our findings may provide new insights into polymorphic variation in ovine PRNT, which may contribute to genetic improvements in economic traits that are important for sheep breeding.     

我国6个地方绵羊品种微卫星DNA多态性研究

利用聚丙烯酰胺凝胶电泳技术研究了我国蒙古羊、乌珠穆沁羊、哈萨克羊、阿勒泰羊、滩羊和藏绵羊 6个地方绵羊品种 17个微卫星标记的多态性 ,以探讨其遗传多样性、起源分化及群体间的遗传亲缘关系。结果表明微卫星标记不同位点间遗传多样性差异极显著 (P <0 0 1) ,群体间多态信息含量 (PIC)、近交程度 (Fis)和观察杂合度 (Obs .Het)差异不显著 ,但基因多样性 (genediversity)和期望杂合度 (Exp .Het)差异显著 (P <0 0 5 )。所研究的我国 6个地方绵羊品种与欧洲品种具有相似的遗传多样性 ,但具有较高的近交系数。个体和群体的聚类分析结果提示我国地方绵羊品种可能起源于两类祖先。群体间的聚类分析结果还表明 ,蒙古羊与乌珠穆沁羊分化不明显且具有较近的遗传亲缘关系 ,蒙古羊与藏绵羊间分化明显且具有较远的遗传亲缘关系。滩羊、阿勒泰羊以及藏绵羊间也具有较近的遗传亲缘关系。所研究的我国 6个地方绵羊品种的遗传分化 (Fst)与西班牙绵羊品种接近 ,但明显小于欧洲其他绵羊品种     

Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds

Background

Traditionally, Chinese indigenous sheep were classified geographically and morphologically into three groups: Mongolian, Kazakh and Tibetan. Herein, we aimed to evaluate the population structure and genome selection among 140 individuals from ten representative Chinese indigenous sheep breeds: Ujimqin, Hu, Tong, Large-Tailed Han and Lop breed (Mongolian group); Duolang and Kazakh (Kazakh group); and Diqing, Plateau-type Tibetan, and Valley-type Tibetan breed (Tibetan group).

Results

We analyzed the population using principal component analysis (PCA), STRUCTURE and a Neighbor-Joining (NJ)-tree. In PCA plot, the Tibetan and Mongolian groups were clustered as expected; however, Duolang and Kazakh (Kazakh group) were segregated. STRUCTURE analyses suggested two subpopulations: one from North China (Kazakh and Mongolian groups) and the other from the Southwest (Tibetan group). In the NJ-tree, the Tibetan group formed an independent branch and the Kazakh and Mongolian groups were mixed. We then used the d_i statistic approach to reveal selection in Chinese indigenous sheep breeds. Among the 599 genome sequence windows analyzed, sixteen (2.7%) exhibited signatures of selection in four or more breeds. We detected three strong selection windows involving three functional genes: RXFP2, PPP1CC and PDGFD. PDGFD, one of the four subfamilies of PDGF, which promotes proliferation and inhibits differentiation of preadipocytes, was significantly selected in fat type breeds by the Rsb (across pairs of populations) approach. Two consecutive selection regions in Duolang sheep were obviously different to other breeds. One region was in OAR2 including three genes (NPR2, SPAG8 and HINT2) the influence growth traits. The other region was in OAR 6 including four genes (PKD2, SPP1, MEPE, and IBSP) associated with a milk production quantitative trait locus. We also identified known candidate genes such as BMPR1B, MSRB3, and three genes (KIT, MC1R, and FRY) that influence lambing percentage, ear size and coat phenotypes, respectively.

Conclusions

Based on the results presented here, we propose that Chinese native sheep can be divided into two genetic groups: the thin type (Tibetan group), and the fat type (Mongolian and Kazakh group). We also identified important genes that drive valuable phenotypes in Chinese indigenous sheep, especially PDGFD, which may influence fat deposition in fat type sheep.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1384-9) contains supplementary material, which is available to authorized users.     

Conservation genetics in Chinese sheep: diversity of fourteen indigenous sheep (Ovis aries) using microsatellite markers

The domestic sheep (Ovis aries) has been an economically and culturally important farm animal species since its domestication around the world. A wide array of sheep breeds with abundant phenotypic diversity exists including domestication and selection as well as the indigenous breeds may harbor specific features as a result of adaptation to their environment. The objective of this study was to investigate the population structure of indigenous sheep in a large geographic location of the Chinese mainland. Six microsatellites were genotyped for 611 individuals from 14 populations. The mean number of alleles (±SD) ranged from 7.00 ± 3.69 in Gangba sheep to 10.50 ± 4.23 in Tibetan sheep. The observed heterozygote frequency (±SD) within a population ranged from 0.58 ± 0.03 in Gangba sheep to 0.71 ± 0.03 in Zazakh sheep and Minxian black fur sheep. In addition, there was a low pairwise difference among the Minxian black fur sheep, Mongolian sheep, Gansu alpine merino, and Lanzhou fat‐tailed sheep. Bayesian analysis with the program STRUCTURE showed support for 3 clusters, revealing a vague genetic clustering pattern with geographic location. The results of the current study inferred high genetic diversity within these native sheep in the Chinese mainland.     

Identification of differentially expressed genes in Mongolian sheep ovaries by suppression subtractive hybridization

Fecundity is an important trait in sheep. Because it is directly related to production costs and efficiency, it has great economic impact in sheep husbandry. Because Mongolian sheep are a longstanding, indigenous breed, they are genetically related to most other breeds of sheep in China. The study of genes related to reproductive traits is essential to improving the fecundity of Mongolian sheep. In the present study, suppression subtractive hybridization (SSH) was performed using forward and reverse nested primers on cDNA libraries from ovarian tissue of single-bearing (S) and biparous (B) Mongolian sheep (MS). This yielded 768 clones. The length of the inserted fragments ranged from 150 to 1000 bp. From these, dot blot hybridization followed by sequencing and homology blast search in GenBank resolved 373 differentially expressed clones, representing 185 gene sequences (homology >85% and length >200 bp), 10 expressed sequence tags (ESTs; homology >95% and length >100 bp), and 4 unknown ESTs. The analysis of the differentially expressed gene functions allowed these genes to be categorized into seven groups: cell/body or immune defense, metabolism, transportation, nucleic acid modification, cell development, signal transduction, and cell structure. Four differentially expressed genes, a disintegrin and metalloproteinase with thrombospondin motifs 1 (ADAMTS1), inhibitor of DNA binding 3 (ID3), bone morphogenetic protein 6 (BMP6), and integrin beta 1 (ITGB1), were randomly selected and verified using relative quantitative real-time polymerase chain reaction (RQ-PCR). The expression of these genes in BMS ovaries was 30.06, 11.55, 0.82, and 1.12-fold that of SMS ovaries, respectively.     

Genetic diversity of Chinese indigenous sheep breeds inferred from microsatellite markers

To determine the genetic diversity and phylogenetic relationships among Chinese sheep, 10 indigenous breeds and one introduced breed were genotyped for 19 microsatellite loci. The mean number of alleles per breed ranged from 5.44 (Guide Black Fur sheep) to 9.13 (Ujumqin sheep and Hulunbeier sheep), the expected heterozygosity varied from 0.623 (Guide Black Fur sheep) to 0.737 (Zhaotong sheep), and the allelic richness ranged from 5.169 (Guide Black Fur sheep) to 7.610 (Zhaotong sheep). The deviation from Hardy–Weinberg equilibrium (HWE) was statistically significant (P < 0.05) at three loci (SRCRSP5, OarAE129 and DYMS1) in most of the breeds. Chinese sheep breeds had maintained a high level of within-population genetic differentiation (95.23%), with the remainder explained by differentiation among populations (4.77%). The genetic differentiation pattern and genetic relationships among Chinese sheep breeds displayed a high consistency with the traditional classification. Both the Bayesian cluster and principal component analyses showed a reliable clustering pattern, which revealed three major clusters in Chinese indigenous sheep (Mongolian sheep, Kazakh sheep and Tibetan sheep), except Zhaotong and Guide Black Fur sheep. There were probably caused by different breeding history, geography isolation and different levels of inbreeding. This study will help to interpret the genetic characters of Chinese indigenous sheep and benefit to the future conservation programs.     

Genetic effects of PRNP gene insertion/deletion (indel) on phenotypic traits in sheep

Prion protein (PRNP) gene is well known for affecting mammal transmissible spongiform encephalopathies (TSE), and is also reported to regulate phenotypic traits (e.g. growth traits) in healthy ruminants. To identify the insertion/deletion (indel) variations of the PRNP gene and evaluate their effects on growth traits, 768 healthy individuals from five sheep breeds located in China and Mongolia were identified and analyzed. Herein, four novel indel polymorphisms, namely, Intron-1-insertion-7bp (I1-7bp), Intron-2-insertion-15bp (I2-15bp), Intron-2-insertion-19bp (I2-19bp), and 3′ UTR-insertion-7bp (3′ UTR-7bp), were found in the sheep PRNP gene. In five analyzed breeds, the minor allelic frequencies (MAF) of the above indels were in the range of 0.008 to 0.986 (I1-7bp), 0.113 to 0.336 (I2-15bp), 0.281 to 0.510 (I2-19bp), and 0.040 to 0.238 (3′ UTR-7bp). Additionally, there were 15 haplotypes and the haplotype ‘I_I2-15bp-D_3’UTR-7bp-D_I2-19bp-D_I1-7bp’ had the highest frequency, which varied from 0.464 to 0.629 in five breeds. Moreover, association analysis revealed that all novel indel polymorphisms were significantly associated with 13 different growth traits (P < 0.05). Particularly, the influences of I2-15bp on chest width (P = 0.001) in Small Tail Han sheep (ewe), 3′ UTR-7bp on chest circumference (P = 0.003) in Hu sheep, and I2-19bp on tail length (P = 0.001) in Tong sheep, were highly significant (P < 0.01). These findings may be a further step toward the detection of indel-based typing within and across sheep breeds, and of promising target loci for accelerating the progress of marker-assisted selection in sheep breeding.     

Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

Background

Artificial selection has caused rapid evolution in domesticated species. The identification of selection footprints across domesticated genomes can contribute to uncover the genetic basis of phenotypic diversity.

Methodology/Main Findings

Genome wide footprints of pig domestication and selection were identified using massive parallel sequencing of pooled reduced representation libraries (RRL) representing ∼2% of the genome from wild boar and four domestic pig breeds (Large White, Landrace, Duroc and Pietrain) which have been under strong selection for muscle development, growth, behavior and coat color. Using specifically developed statistical methods that account for DNA pooling, low mean sequencing depth, and sequencing errors, we provide genome-wide estimates of nucleotide diversity and genetic differentiation in pig. Widespread signals suggestive of positive and balancing selection were found and the strongest signals were observed in Pietrain, one of the breeds most intensively selected for muscle development. Most signals were population-specific but affected genomic regions which harbored genes for common biological categories including coat color, brain development, muscle development, growth, metabolism, olfaction and immunity. Genetic differentiation in regions harboring genes related to muscle development and growth was higher between breeds than between a given breed and the wild boar.

Conclusions/Significance

These results, suggest that although domesticated breeds have experienced similar selective pressures, selection has acted upon different genes. This might reflect the multiple domestication events of European breeds or could be the result of subsequent introgression of Asian alleles. Overall, it was estimated that approximately 7% of the porcine genome has been affected by selection events. This study illustrates that the massive parallel sequencing of genomic pools is a cost-effective approach to identify footprints of selection.     

Detection of nitrogen deficiency QTL in juvenile wild barley introgression linesgrowing in a hydroponic system

Background

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

Results

Two independent experiments using either Iranian or Ovine HapMap genotyping data contrasted thin and fat tail breeds. Population differentiation using F_ST in Iranian thin and fat tail breeds revealed seven genomic regions. Almost all of these regions overlapped with QTLs that had previously been identified as affecting fat and carcass yield traits in beef and dairy cattle. Study of selection sweep signatures using F_ST in thin and fat tail breeds sampled from the Ovine HapMap project confirmed three of these regions located on Chromosomes 5, 7 and X. We found increased homozygosity in these regions in favour of fat tail breeds on chromosome 5 and X and in favour of thin tail breeds on chromosome 7.

Conclusions

In this study, we were able to identify three novel regions associated with fat deposition in thin and fat tail sheep breeds. Two of these were associated with an increase of homozygosity in the fat tail breeds which would be consistent with selection for mutations affecting fat tail size several thousand years after domestication.     

Genetic variation of 5 SNPs of MC1R gene in Chinese indigenous sheep breeds

The purpose of this study was to assess genetic diversity, genetic differentiation relationship and population structure among 10 Chinese sheep populations using 5 single nucleotide polymorphisms (SNPs) in MC1R gene. The genetic diversity indices suggested that the intra-population variation levels of Chinese Merino and Large-tailed Han breeds were lowest than Kazakh Fat-Rumped. Chinese sheep breeds have maintained a high intra-population variation levels (95.23%). The genetic differentiation patterns and genetic relationships among Chinese sheep breeds displayed a high consistency with the traditional classification. The cluster trees were constructed by UPMGA method. The results showed that Chinese indigenous sheep populations have distinct genetic differentiation. The inter-population variation levels in Chinese sheep populations indicated three geographically independent domestication events have occurred. The Bayesian cluster analyses also showed a reliable clustering pattern, which revealed three major clusters in Chinese indigenous sheep populations (Mongolian group, Kazakh group and Tibetan group), except for Duolang and Minxian Black-fur. There were probably caused by different breeding history, geography isolation and different levels of inbreeding. The findings supported the related records in literature, ten sheep populations originated on different time stage from the primogenitor population and communicated genetically with each other in the process of natural and artificial selection, and in different ecological environment. It is concluded that Chinese indigenous sheep have higher genetic variation and diversity, genetic differentiation exist between Chinese sheep populations.     

Genome‐wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

Pigs from Asia and Europe were independently domesticated from c. 9000 years ago. During this period, strong artificial selection has led to dramatic phenotypic changes in domestic pigs. However, the genetic basis underlying these morphological and behavioural adaptations is relatively unknown, particularly for indigenous Chinese pigs. Here, we performed a genome‐wide analysis to screen 196 regions with selective sweep signals in Tongcheng pigs, which are a typical indigenous Chinese breed. Genes located in these regions have been found to be involved in lipid metabolism, melanocyte differentiation, neural development and other biological processes, which coincide with the evolutionary phenotypic changes in this breed. A synonymous substitution, c.669T>C, in ESR1, which colocalizes with a major quantitative trait locus for litter size, shows extreme differences in allele frequency between Tongcheng pigs and wild boars. Notably, the variant C allele in this locus exhibits high allele frequency in most Chinese populations, suggesting a consequence of positive selection. Five genes (PRM1, PRM2, TNP2, GPR149 and JMJD1C) related to reproductive traits were found to have high haplotype similarity in Chinese breeds. Two selected genes, MITF and EDNRB, are implied to shape the two‐end black colour trait in Tongcheng pig. Subsequent SNP microarray studies of five Chinese white‐spotted breeds displayed a concordant signature at both loci, suggesting that these two genes are responsible for colour variations in Chinese breeds. Utilizing massively parallel sequencing, we characterized the candidate sites that adapt to artificial and environmental selections during the Chinese pig domestication. This study provides fundamental proof for further research on the evolutionary adaptation of Chinese pigs.     

中国蒙系绵羊品种间的基因流动

利用结构基因座和微卫星标记分析了中国蒙系绵羊4个品种间的基因流动情况。结果显示, 由结构基因座得出的遗传分化系数在0.0164 ~ 0.0455之间, 由微卫星位点得出的遗传分化系数在0.0107 ~ 0.0239之间, 说明遗传变异绝大部分存在于品种内, 蒙系绵羊品种间的遗传分化水平很低。结构基因座(N_m = 7.971)和微卫星位点(N_m = 15.732)都反映出品种间的基因流通畅。品种间的遗传差异与地理距离之间无直接相关关系。初步推断我国蒙系绵羊品种间的遗传分化主要是自然选择(生境异质性)作用的结果。     

A genome-wide scan for signatures of directional selection in domesticated pigs

Background

Animal domestication involved drastic phenotypic changes driven by strong artificial selection and also resulted in new populations of breeds, established by humans. This study aims to identify genes that show evidence of recent artificial selection during pig domestication.

Results

Whole-genome resequencing of 30 individual pigs from domesticated breeds, Landrace and Yorkshire, and 10 Asian wild boars at ~16-fold coverage was performed resulting in over 4.3 million SNPs for 19,990 genes. We constructed a comprehensive genome map of directional selection by detecting selective sweeps using an F_ST-based approach that detects directional selection in lineages leading to the domesticated breeds and using a haplotype-based test that detects ongoing selective sweeps within the breeds. We show that candidate genes under selection are significantly enriched for loci implicated in quantitative traits important to pig reproduction and production. The candidate gene with the strongest signals of directional selection belongs to group III of the metabolomics glutamate receptors, known to affect brain functions associated with eating behavior, suggesting that loci under strong selection include loci involved in behaviorial traits in domesticated pigs including tameness.

Conclusions

We show that a significant proportion of selection signatures coincide with loci that were previously inferred to affect phenotypic variation in pigs. We further identify functional enrichment related to behavior, such as signal transduction and neuronal activities, for those targets of selection during domestication in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1330-x) contains supplementary material, which is available to authorized users.     

Signatures of Diversifying Selection in European Pig Breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.     

Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

绵羊线粒体DNA D-loop区串联重复序列变异研究

对来自69个我国地方绵羊品种和8个国外引入品种共计77个个体线粒体DNA控制区长度为75bp的串联重复序列进行了测序分析。在309个重复序列中检测到28个变异位点,其中7个为具有2个变异体的单现突变,1个为具有3个变异体的单现突变,20个为具有2个变异体的简约位点。由28个变异位点中归纳出63个单倍型,其中单倍型Ⅰ和单倍型Ⅲ具有较高的比例,分别为12．94％和30.42％。研究结果揭示我国地方绵羊可能起源于两个母系祖先。哈萨克羊和阿勒泰羊间以及蒙古羊和乌珠穆沁羊间分别具有较近的亲缘关系且没有明显的遗传分化。藏绵羊、蒙古羊和乌珠穆沁羊相对哈萨克羊和阿勒泰羊而言具有较低的遗传多样性。     

Study on Tandem Repeat Sequence Variation in Sheep mtDNA D-loop Region

The 75-nt-long tandem repeat sequence in the control region of mtDNA of 77 individuals, of which 69 were from different indigenous sheep breeds in China and 8 were from imported breeds, was sequenced and analyzed to investigate the origin and differentiation of Chinese indigenous sheep breeds and also the genetic diversities and relationships among them. A total of 28 variable sites were detected within 309 repeated sequences, among which 7 sites were singleton variable sites with two variants, 1 site was a singleton variable site with three variants, and 20 sites were parsimony informative sites with two variants. A total of 63 haplotypes were sorted from 28 polymorphic sites, among which two main and basic haplotypes, namely, Hap 1 and Hap 3 were present at a much higher proportion, at 12.94% and 30.42%, respectively. It could be inferred that Chinese indigenous sheep breeds originated from two maternal ancestors because of the maternal inheritance characteristics of the mtDNA. Altay sheep and Kazakstan sheep are closely related and do not differentiate significantly. Mongolian sheep and Ujumuqin sheep also share a close relationship. Tibetan sheep, Mongolian sheep, and Ujumuqin sheep have lower genetic diversity than Altay sheep and Kazakstan sheep.