Vascular morphology and expression of angiogenic factors in sheep ovaries and their relation with the lambing rate.

The study aims to present the vascular anatomy of ovarian arterial system in sheep and to explore its relation with the angiogenic factors vascular endothelial growth factor (VEGF) and angiopoietin-2 (Ang-2), environmental factors and the lambing rate. Small-tailed Han sheep, Hu sheep, Aohan finewool sheep and Tibetan sheep who are from different habitats in China and with different lambing rates were investigated. The vascular corrosion casting technique and scanning electron microscopy were used to understand arterial architecture, and real-time quantitative PCR (qRT-PCR) and Western blotting were used for assessing VEGF and Ang-2 expression. Results showed that high-prolificity breeds had significantly larger vessel diameters in ovaries than low-prolificity ones. VEGF gene transcript and protein expression were related with the lambing rate, while Ang-2 gene transcript was negatively related with the lambing rate. Results indicated that high-prolificity breeds have more abundant arterial blood supply, which is associated with VEGF expression.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success.

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA-miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm-egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

Low genome-wide homozygosity in 11 Spanish ovine breeds.

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds.

A dataset consisting of 787 animals with high-density SNP chip genotypes (346 774 SNPs) and 939 animals with medium-density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between FPED and FROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean FROH estimates from medium-density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make-up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.

Inferring the population structure of the Maghreb sheep breeds using a medium-density SNP chip.

North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome-wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat-tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.

Detecting selection signatures in three Iranian sheep breeds.

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype-phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP-EHH (cross-population extend haplotype homozygosity) and FST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori-Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations.

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short-tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D5 ) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.

Identification of two TYRP1 loss-of-function alleles in Valais Red sheep.

The Valais Red sheep breed is a local breed of the Swiss canton Valais. Although the breed is characterised by its brown colour, black animals occasionally occur and the objective of this study was to identify the causative genetic variants responsible for the obvious difference. A GWAS using high-density SNP data to compare 51 brown and 38 black sheep showed a strong signal on chromosome 2 at the TYRP1 locus. Haplotype analyses revealed three different brown-associated alleles. The WGS of three sheep revealed four protein-changing variants within the TYRP1 gene. Three of these variants were associated with the recessively inherited brown coat colour. This includes the known missense variant TYRP1:c.869G>T designated as bS oay and two novel loss-of-function variants. We propose to designate the frame-shift variant TYRP1:c.86_87delGA as bVS 1 and the nonsense variant TYRP1:c.1066C>T as bVS 2 . Interestingly, the bVS 1 allele occurs only in local breeds of Switzerland whereas the bVS 2 allele seems to be more widespread across Europe.

Genome-wide association studies for milk production traits in Valle del Belice sheep using repeated measures.

Genome-wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test-day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.

Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep.

T-box transcription factor T (TBXT), encoding the brachyury protein, is an embryonic nuclear transcription factor involved in mesoderm formation and differentiation. Previous studies indicate that TBXT mutations are responsible for the tailless or short-tailed phenotype of many vertebrates. To verify whether the tailless phenotype in fat-rumped sheep is associated with TBXT mutations, exon 2 of the TBXT gene for 301 individuals belonging to 13 Chinese and Iranian sheep breeds was directly sequenced. Meanwhile, 380 samples were used to detect the genotypes of the candidate variations by mapping to their reads databases in the Sequence Read Archive repository of GenBank. The results showed that one missense mutation, c.334G>T (GGG>TGG) with a completely linked synonymous variant c.333G>C (CCG>CCC) was found to be associated with the 'tailless' characteristic in typical fat-rumped sheep breeds. The c.334G>T transversion led to the conversion of glycine to tryptophan at the 112th amino acid in the T-box domain of the brachyury protein. In addition, crossbreeding experiments for long-tailed and tailless sheep showed that CT/CT allele of nucleotides (nt) 333 and 334, a recessive mutation, would cause sheep tails to be shorter, suggesting that these two linked variants at nucleotides 333 and 334 in TBXT are probably causative mutations responsible for the tailless phenotype in sheep.

Genetic differentiation and phylogeny of 27 sheep populations based on structural gene loci.

To explore the genetic divergence and phylogeny of Chinese indigenous sheep breeds, in the current study, we analyzed the polymorphisms of 5 structural loci in ten sheep populations, including Sishui Fur, Sunite, Wurank, Bayinbuluke, Altay, Small-Tailed Han, Wadi, Tan, Tong and Hu sheep. The data were then compared with those from an additional 13 Asian and 4 European sheep populations acquired by the same experimental method. Based on the genetic distance and the results of a cluster analysis, we constructed the phylogenetic relationship of 27 populations. The results showed that the sheep populations in this study could be classified into four genetic groups: "Mongolian", "Tibetan", "South-Southeast Asian" and "European" sheep groups. All 10 Chinese sheep breeds belonged to the "Mongolian sheep" lineage; however, Finnish Landrace sheep and Yunnan sheep could not be classified into any of the four groups. These results could provide a good reference for the protection and utilization of primary breed resources in China and phylogenic research on Asian sheep populations.

Genetic selection for increased mean and reduced variance of twinning rate in Belclare ewes.

It is sometimes possible to breed for more uniform individuals by selecting animals with a greater tendency to be less variable, that is, those with a smaller environmental variance. This approach has been applied to reproduction traits in various animal species. We have evaluated fecundity in the Irish Belclare sheep breed by analyses of flocks with differing average litter size (number of lambs per ewe per year, NLB) and have estimated the genetic variance in environmental variance of lambing traits using double hierarchical generalized linear models (DHGLM). The data set comprised of 9470 litter size records from 4407 ewes collected in 56 flocks. The percentage of pedigreed lambing ewes with singles, twins and triplets was 30, 54 and 14%, respectively, in 2013 and has been relatively constant for the last 15 years. The variance of NLB increases with the mean in this data; the correlation of mean and standard deviation across sires is 0.50. The breeding goal is to increase the mean NLB without unduly increasing the incidence of triplets and higher litter sizes. The heritability estimates for lambing traits were NLB, 0.09; triplet occurrence (TRI) 0.07; and twin occurrence (TWN), 0.02. The highest and lowest twinning flocks differed by 23% (75% versus 52%) in the proportion of ewes lambing twins. Fitting bivariate sire models to NLB and the residual from the NLB model using a double hierarchical generalized linear model (DHGLM) model found a strong genetic correlation (0.88 ± 0.07) between the sire effect for the magnitude of the residual (VE ) and sire effects for NLB, confirming the general observation that increased average litter size is associated with increased variability in litter size. We propose a threshold model that may help breeders with low litter size increase the percentage of twin bearers without unduly increasing the percentage of ewes bearing triplets in Belclare sheep.

The genetic diversity and structure of 18 sheep breeds exposed to isolation and selection.

The phylogenetic layout of the genotyped (30 microsatellite) 18 sheep breeds in this study demands and provides the opportunity to evaluate both neutral and adaptive components of genetic diversity in a naturally and artificially selected and subdivided sheep population. Seven Pramenka strains from Bosnia and Herzegovina and Croatia characterized by a very low intensity of artificial selection, preserved the highest neutral genetic variability. Eight central and north-western European breeds under considerable artificial isolation and selection preserved the lowest genetic variability. Only combinations of various phylogenetic parameters offer a reasonable explanation for underlying evolutionary forces working in the investigated island and mainland sheep breeds under variable natural and artificial selection. More than 60% of total genetic, diversity was allocated to virtually unselected Pramenka strains, and an additional 25% to native moderately selected Graue Gehoernte Heidschnucke and intensively selected Ostfriesische Milchschafe. Some economically very important breeds and strains did not contribute to a pool with maximal genetic diversity, while they play an important role in the cultural heritage of respective countries.

Design of a low-density SNP chip for the main Australian sheep breeds and its effect on imputation and genomic prediction accuracy.

Genotyping sheep for genome-wide SNPs at lower density and imputing to a higher density would enable cost-effective implementation of genomic selection, provided imputation was accurate enough. Here, we describe the design of a low-density (12k) SNP chip and evaluate the accuracy of imputation from the 12k SNP genotypes to 50k SNP genotypes in the major Australian sheep breeds. In addition, the impact of imperfect imputation on genomic predictions was evaluated by comparing the accuracy of genomic predictions for 15 novel meat traits including carcass and meat quality and omega fatty acid traits in sheep, from 12k SNP genotypes, imputed 50k SNP genotypes and real 50k SNP genotypes. The 12k chip design included 12 223 SNPs with a high minor allele frequency that were selected with intermarker spacing of 50-475 kb. SNPs for parentage and horned or polled tests also were represented. Chromosome ends were enriched with SNPs to reduce edge effects on imputation. The imputation performance of the 12k SNP chip was evaluated using 50k SNP genotypes of 4642 animals from six breeds in three different scenarios: (1) within breed, (2) single breed from multibreed reference and (3) multibreed from a single-breed reference. The highest imputation accuracies were found with scenario 2, whereas scenario 3 was the worst, as expected. Using scenario 2, the average imputation accuracy in Border Leicester, Polled Dorset, Merino, White Suffolk and crosses was 0.95, 0.95, 0.92, 0.91 and 0.93 respectively. Imputation scenario 2 was used to impute 50k genotypes for 10 396 animals with novel meat trait phenotypes to compare genomic prediction accuracy using genomic best linear unbiased prediction (GBLUP) with real and imputed 50k genotypes. The weighted mean imputation accuracy achieved was 0.92. The average accuracy of genomic estimated breeding values (GEBVs) based on only 12k data was 0.08 across traits and breeds, but accuracies varied widely. The mean GBLUP accuracies with imputed 50k data more than doubled to 0.21. Accuracies of genomic prediction were very similar for imputed and real 50k genotypes. There was no apparent impact on accuracy of GEBVs as a result of using imputed rather than real 50k genotypes, provided imputation accuracy was >90%.

Y chromosome haplotype diversity of domestic sheep (Ovis aries) in northern Eurasia.

Variation in two SNPs and one microsatellite on the Y chromosome was analyzed in a total of 663 rams representing 59 breeds from a large geographic range in northern Eurasia. SNPA-oY1 showed the highest allele frequency (91.55%) across the breeds, whereas SNPG-oY1 was present in only 56 samples. Combined genotypes established seven haplotypes (H4, H5, H6, H7, H8, H12 and H19). H6 dominated in northern Eurasia, and H8 showed the second-highest frequency. H4, which had been earlier reported to be absent in European breeds, was detected in one European breed (Swiniarka), whereas H7, which had been previously identified to be unique to European breeds, was present in two Chinese breeds (Ninglang Black and Large-tailed Han), one Buryatian (Transbaikal Finewool) and two Russian breeds (North Caucasus Mutton-Wool and Kuibyshev). H12, which had been detected only in Turkish breeds, was also found in Chinese breeds in this work. An overall low level of haplotype diversity (median h = 0.1288) was observed across the breeds with relatively higher median values in breeds from the regions neighboring the Near Eastern domestication center of sheep. H6 is the dominant haplotype in northwestern and eastern China, in which the haplotype distribution could be explained by the historical translocations of the H4 and H8 Y chromosomes to China via the Mongol invasions followed by expansions to northwestern and eastern China. Our findings extend previous results of sheep Y chromosomal genetic variability and indicate probably recent paternal gene flows between sheep breeds from distinct major geographic regions.

A revised phylogeny of Antilopini (Bovidae, Artiodactyla) using combined mitochondrial and nuclear genes.

Antilopini (gazelles and their allies) are one of the most diverse but phylogenetically controversial groups of bovids. Here we provide a molecular phylogeny of this poorly understood taxon using combined analyses of mitochondrial (CYTB, COIII, 12S, 16S) and nuclear (KCAS, SPTBN1, PRKCI, MC1R, THYR) genes. We explore the influence of data partitioning and different analytical methods, including Bayesian inference, maximum likelihood and maximum parsimony, on the inferred relationships within Antilopini. We achieve increased resolution and support compared to previous analyses especially in the two most problematic parts of their tree. First, taxa commonly referred to as "gazelles" are recovered as paraphyletic, as the genus Gazella appears more closely related to the Indian blackbuck (Antilope cervicapra) than to the other two gazelle genera (Nanger and Eudorcas). Second, we recovered a strongly supported sister relationship between one of the dwarf antelopes (Ourebia) and the Antilopini subgroup Antilopina (Saiga, Gerenuk, Springbok, Blackbuck and gazelles). The assessment of the influence of taxon sampling, outgroup rooting, and data partitioning in Bayesian analyses helps explain the contradictory results of previous studies.

Mitochondrial analysis sheds light on the origin of hair sheep.

A total of 180 mtDNA sequences from hair Caribbean (93), West African (73) and Canarian-wooled (14) sheep were analysed to shed light on the origin of hair sheep. A comparison of 360 Iberian sheep sequences retrieved from GenBank was performed to assess a possible European origin of the Caribbean hair sheep. These 180 sequences gave 48 different haplotypes (16 in Caribbean sheep). All Caribbean and Canarian-wooled sequences and 91.8% of the West African samples belonged to haplogroup B. The sheep analysed showed wide haplotypic identity. Caribbean sheep shared roughly two-thirds of their samples with Canarian-wooled and West African samples, respectively. Principal component analysis showed that the Caribbean and the Canarian-wooled sheep clustered together. Additional analyses showed that hair and Iberian sheep had wide genetic identity. It was not possible to ascertain a single Canarian, African or European origin of the Caribbean hair sheep using mtDNA markers only. European, African and Caribbean hair sheep maternal genetic backgrounds likely result from related domestication events.

The ACACA gene is a potential candidate gene for fat content in sheep milk.

No major gene has yet been reported in sheep that explains the variation of milk fat content. The coding region of the acetyl-CoA carboxylase alpha (ACACA) gene, which plays an important role in de novo fatty acid synthesis, had been investigated, but no non-synonymous mutations have been reported. In this study, the genomic regions encoding the three promoters of the ACACA gene were directly sequenced in 264 sheep of three different breeds, and 10 SNPs were identified. Allele frequencies of most SNPs significantly differed (P = 0.05-0.0001) between breeds. The SNPs that potentially altered either gene regulatory elements or putative binding sites of transcription factors were made evident through in silico analysis. The association analysis with milk traits, performed for one SNP of PIII (GenBank AJ292286, g.1330G>T), showed a significant allelic substitution effect (+0.33%, P < 0.0001 and +0.35%, P < 0.01) in the Altamurana and Gentile breeds respectively. Because this SNP was located in the binding site of the paired box protein transcription factors, which was shown to function as an efficient promoter element, and because PIII transcripts are expressed in the mammary gland, the SNP in PIII of the ACACA gene might affect the variation of fat content in sheep milk.

Novel polymorphisms in ovine prion protein gene.

The aim of this study was to identify the PRNP polymorphisms outside the standard codons 136, 154 and 171 in 1110 sheep with no clinical sign of scrapie from all 18 Turkish native sheep breeds and compare our results with published data on ovine PRNP polymorphism from other regions of the world. Among the 22 amino acid polymorphisms and three silent mutations, 10 were novel for ovine PRNP: p.Gly94Gly, p.Leu128Ile, p.Met132Leu, p.Ser135Arg, p.Met137Val, p.Asn146Lys, p.Arg159Arg, p.Tyr160Asn, p.Gln163His and p.Thr193Ser. These data reveal that sheep breeds close to the historic center of small ruminant domestication have remained highly diverse in the prion gene locus, with distinctive genetic similarities to both Asian and European sheep breeds.

SNP variation in ADRB3 gene reflects the breed difference of sheep populations.

The ß3-adrenergic receptor (ADRB3), a G-protein coupled receptor, plays a major role in energy metabolism and regulation of lipolysis and homeostasis. We detect the single nucleotide polymorphism (SNP) variation in full-length sequence of ovine ADRB3 gene in 12 domestic sheep populations within four types by polymerase chain reaction-single strand conformation polymorphism and sequencing to reveal the breed difference. Twenty-two SNPs, 12 of which in the exon 1 and ten in the intron, were detected, and 12 new exonic and four new intronic SNPs were found. Most SNPs presented in Shanxi Dam Line and least ones in Dorset. The average SNP number in both meat and dual purpose for meat and wool breeds was significantly higher than general and dual purpose breeds for wool and meat. Frequency of each SNP in studied breeds or types was different. The 18C Del and 1617T Ins majorly existed in dual purpose breeds for wool and meat. The 25A Del, 119C>G and 130C>T were mostly found in the meat and dual purpose for meat and wool breeds. The 1764C>A more frequently presented in meat than in other types. The majority of variations came from within the populations as suggested by analysis of molecular variance. Close relationship presented among the Chinese and western breeds, respectively. In conclusion, SNPs of ovine ADRB3 gene can reflect the breed difference and within- and between-population variations, and to a great extent, the breed relationship.