Evidence for a third allele at the beta-lactoglobulin (beta-Lg) locus of sheep milk and its occurrence in different breeds

Milk samples from 189 Merinoland Sheep, 145 Black Faced Mutton Sheep, 89 East Friesian Milk Sheep, 36 Rh?n, 36 Pleven, 23 Tsigaja, 25 Black Razka and 86 Hungarian Merino X Pleven (F1) sheep were analysed by polyacrylamide gel electrophoresis under acid conditions and isoelectric focusing in ultrathin layer polyacrylamide gels with carrier ampholytes. Six different beta-lactoglobulin (beta-Lg) phenotypes (A, AB, B, AC, BC and C) were observed by both methods. The occurrence of three codominant alleles (beta-LgA, beta-LgB, beta-LgC) at an autosomal locus (beta-Lg) was supported by family and population data on genetic equilibrium. Differences in gene frequencies between the breeds were observed.     

Evidence for a third allele at the β-lactoglobulin (β-Lg) locus of sheep milk and its occurrence in different breeds

Summary. Milk samples from 189 Merinoland Sheep, 145 Black Faced Mutton Sheep, 89 East Friesian Milk Sheep, 36 Rhön, 36 Pleven, 23 Tsigaja, 25 Black Razka and 86 Hungarian Merino x Pleven (F1) sheep were analysed by polyacrylamide gel electrophoresis under acid conditions and isoelectric focusing in ultrathin layer polyacrylamide gels with carrier ampholytes. Six different β-lactoglobulin (β-Lg) phenotypes (A, AB, B, AC, BC and C) were observed by both methods. The occurrence of three codominant alleles (β-Lg^A, β-Lg^B, β-Lg^c) at an autosomal locus (β-Lg) was supported by family and population data on genetic equilibrium. Differences in gene frequencies between the breeds were observed.     

Evaluation of the genetic variability of 23 bovine microsatellite markers in four Belgian cattle breeds

The polymorphism of 23 microsatellites in the four main cattle breeds in Belgium (Holstein Friesian, Belgian Blue, Belgian Red Pied and East Flemish) was analysed. Heterozygosity, polymorphism information content, the effective number of alleles, exclusion probability and the probability of genotypic identity for two random individuals were calculated for all microsatellites and all breeds. The Belgian Blue breed is generally a little less polymorphic in comparison with the other three breeds. Estimates of the genetic distances between these breeds confirmed the widely accepted proposition that the Belgian Blue is the most genetically distinct of these breeds. The three other breeds are likely to become one population, given current breeding strategies. Exclusion probabilities in parentage control cases are >0·9999 in all four breeds when all 23 microsatellites are used and >0·98 with only the two most polymorphic multiplexes.     

Predictive potential of microsatellite markers on heterosis of fecundity in crossbred sheep

Small Tail Han (STH) sheep is a famous Chinese local breed and has perfect prolificacy performance, but it is inferior to imported mutton sheep breeds on meat production. In this study, six imported male sheep populations (White Suffolk, Black Suffolk, Texel, Dorper, South African Mutton Merino and East Friesian) were crossbred with STH female sheep respectively. The heterosis values of litter size, average daily gain (ADG) and feed conversion ratio (FCR) of crossbred sheep were analyzed for seeking the optimal cross. Meanwhile 28 microsatellite markers were used to measure the genetic distance between imported populations and STH population. Regression between the genetic distance and heterosis was analyzed for evaluating potential of microsatellite on predicting heterosis. Results showed a significant positive linear correlation (r = 0.892, P < 0.05) between heterosis of litter size and genetic distance D _A of six crosses. This implied that these microsatellite markers had moderate potential to forecast heterosis of litter size in sheep. Results of this study also indicated that South African Mutton Merino and East Friesian sheep would be the optimal sire breeds for the litter size and might bring the greatest economic benefit in six imported populations; Suffolk sheep could be prior consideration as sire breeds when breeding objective focused on ADG. Finally these results provided valuable information for Chinese sheep industry.     

The pattern of runs of homozygosity and genomic inbreeding in world-wide sheep populations

Genome-wide pattern of runs of homozygosity (ROH) across ovine genome can provide a useful resource for studying diversity and demography history in sheep. We analyzed 50 k SNPs chip data of 2536 animals to identify pattern, distribution and level of ROHs in 68 global sheep populations. A total of 60,301 ROHs were detected in all breeds. The majority of the detected ROHs were <16 Mb and the average total number of ROHs per individual was 23.8 ± 13.8. The ROHs greater than 1 Mb covered on average 8.2% of the sheep autosomes, 1% of which was related to the ROHs with 1–4 Mb of length. The mean sum of ROH length in two-thirds of the populations was less than 250 Mb ranging from 21.7 to near 570 Mb. The level of genomic inbreeding was relatively low. The average of the inbreeding coefficients based on ROH (F_ROH) was 0.09 ± 0.05. It was rising in a stepwise manner with distance from Southwest Asia and maximum values were detected in North European breeds. A total of 465 ROH hotspots were detected in 25 different autosomes which partially surrounding 257 Refseq genes across the genome. Most of the detected genes were related to growth, body weight, meat production and quality, wool production and pigmentation. In conclusion, our analysis showed that the sheep genome, compared with other livestock species such as cattle and pig, displays low levels of homozygosity and appropriate genetic diversity for selection response and genetic merit gain.     

A review on the potential effects of environmental and economic factors on sheep genetic diversity: Consequences of climate change

Climate change has a significant effect on the productivity of livestock including milk, meat, and reproduction. This could be attributed to the internal diversion of energy resources towards adaptive mechanisms. Among the climate change variables, thermal stress seems to be the major limiting factor in animal agriculture. A better understanding of the effects of climate change-influenced ecological factors on the genetic diversity of livestock species is warranted. Sheep is an ideal livestock species to be used in investigating environmental adaptation due to its wide range of agroecological habitats, genetic and phenotypic variability. There is a heavy reliance on sheep genetic diversity for future animal protein security, but the implications of climate change on their genetic diversity receive less attention.Here, the potential environmental factors influencing natural selection in sheep populations are presented. We argue that prolonged exposure to these factors plays a major role in influencing the development of adaptation traits in indigenous sheep breeds, consequently leading to the alteration of genetic diversity at specific loci. The factors discussed include hot temperatures (heat stress), insufficient water, low quantity and quality of forage, and prevalence of parasites, pests, and diseases. In addition, genetic diversity, some signatures of selection for adaptation and economic angles of selection are also briefly discussed.A better understanding of environmental factors influencing the genetic diversity of sheep populations will inform breeding and management programs and may offer an opportunity for greater production efficiency with low input costs.     

Genetic characterization and breed assignment in five Italian sheep breeds using microsatellite markers

The knowledge of the genetic relationship and admixture among neighbouring populations is crucial for conservation efforts. The aim of this study was to analyse the genetic diversity of five Italian sheep breeds (Appenninica, Garfagnina Bianca, Massese, Pomarancina and Zerasca) using a panel of 24 microsatellite markers. Blood samples from 226 individuals belonging to the aforementioned populations were obtained and genotyped. All the investigated breeds showed a significant heterozygote deficiency caused by the high level of inbreeding indicated also by the high level of F_IS (0.146). Genetic differentiation between breeds was moderate (F_ST = 0.05) but significant and the individuals could be assigned to their breeds with an high success rate even if the inter-individual distances showed that few animals clustered separately from the other individuals of the same breed, especially for Pomarancina breed. The genetic distances reflect the historical knowledge of these breeds and some patterns of ancestral and recent gene flow between neighbour populations arise. The clustering analysis detects the presence of six clusters. Massese and Zerasca breeds were grouped together as well as Appenninica and Pomarancina with the latter forming two distinct clusters equally represented. The formation of this last breed is occurred with the absorption of individuals of the Appenninica breed and the gene flow probably continued in these recent years allowing the presence of a population substructure for Pomarancina breed. Such substructure supports the high level of heterozygote deficiency found for this breed despite the relatively high population size. The five populations analysed presented some genetic similarities but a clear uniqueness of the populations has been showed for almost all of them. Special attention to monitor genetic variability and to organize mating plans should be given especially for the three endangered breeds.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

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 (F_ROH) are generally small (F_ROH < 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.     

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale,Merino and Creole sheep

The aim of this study was to investigate the genetic diversity within and among three breeds of sheep: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip^®. Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (F_ST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (F_ST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (F_ST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.     

Genetic relationships between Spanish Assaf (Assaf.E) and Spanish native dairy sheep breeds

At present, the Assaf is the main dairy sheep in Spain. The Spanish Assaf (Assaf.E) was formed by male-mediated absorption of native Spanish sheep. Here we assess the genetic relationships among the Assaf.E and major native Spanish dairy breeds using microsatellites to contribute to the knowledge of the formation and within-population genetic variability of the breed. Blood samples from 44 unrelated Assaf.E individuals from 23 different Assaf.E flocks spread throughout 6 different Spanish provinces were obtained and genotyped using 14 microsatellites. Up to 312 additional samples belonging to the Awassi and Milchschaf sheep breeds and to six native Spanish dairy sheep breeds (Castellana, Churra, Latxa, Manchega, and Rubia de El Molar) as well as samples from Merino individuals to be used as the outgroup were also analysed observed (H_o) and expected (H_e) heterozygosity, rarefacted number of alleles per locus and distances based on molecular coancestry information were computed. Probabilities of assignment of the Assaf.E individuals to native Spanish dairy sheep breeds and cryptic genetic structure in the whole dataset were also assessed. It can be concluded that the Assaf.E breed has low genetic variability and high genetic distance with respect native Spanish dairy sheep breeds. From our results, the formation of the Assaf.E breed basically occurred via the absorption of individuals belonging to the Entrefino type, particularly to the Castellana and Manchega populations. Furthermore, Churra individuals may have participated in the formation of the Assaf.E breed at an early moment of the introduction of the breed into Spain.     

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of F_ST were higher in Wallachian sheep (F_ST = 0.14) than in Sumava sheep (F_ST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the N_e estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (F_HOM) among individuals, from ?0.04 to 0.16 in Sumava and from ?0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.     

Genome-wide scan reveals important additive and non-additive genetic effects associated with resistance to Haemonchus contortus in Florida Native sheep

Florida Native sheep is among the sheep breeds best adapted to humid and hot climatic conditions such as those of Florida, USA, and have shown a superior ability to regulate nematode burdens. This is one of the oldest sheep breeds in North America and is an endangered species. To ensure genetic diversity and long-term survival of the breed, protection of the current genetic stock is critical and conservation efforts are required to promote its breeding and production. The objective of the present study was to investigate the importance of additive and non-additive genetic effects on resistance to natural Haemonchus contortus infections in Florida Native sheep using a whole genome scan. A total of 200 sheep were evaluated in the present study. Phenotypic records included faecal egg count (FEC, eggs/gram), FAMACHA® score, packed cell volume (PCV, %), body condition score and average daily gain (ADG, kg). Sheep were genotyped using the GGP Ovine 50K SNP chip and 45.2 k single nucleotide polymorphism (SNP) markers spanning the entire genome were available for quality control procedures. Mixed models were used to analyse the response variables and included the identity by state matrix to control for population structure. Bonferroni correction was used to control for multiple testing and a second arbitrary threshold (0.1 × 10^?3) was used. Fifteen SNPs with additive and non-additive genetic effects and located in Ovis aries chromosome OAR1, 2, 3, 6, 8, 10, 11, 12, 13 and 21 were associated with FEC, FAMACHA® score, PCV and ADG. These SNPs could be potential genetic markers for resistance to natural H. contortus exposure in Florida Native sheep.     

Genome-wide DNA arrays profiling unravels the genetic structure of Iranian sheep and pattern of admixture with worldwide coarse-wool sheep breeds

Archaeological and genetic evidence show that sheep were originally domesticated in area around the North of Zagros mountains, North-west of Iran. The Persian plateau exhibits a variety of native sheep breeds with a common characteristic of coarse-wool production. Therefore, knowledge about the genetic structure and diversity of Iranian sheep and genetic connections with other sheep breeds is of great interest. To this end, we genotyped 154 samples from 11 sheep breeds distributed across Iran with the Ovine Infinium HD SNP 600 K BeadChip array, and analyzed this dataset combined with the retrieved data of 558 samples from 19 worldwide coarse-wool sheep breeds. The average genetic diversity ranged from 0.315 to 0.354, while the F_ST values ranged from 0.016 to 0.177 indicating a low differentiation of Iranian sheep. Analysis of molecular variance showed that 90.21 and 9.79% of the source of variation were related to differences within and between populations, respectively. Our results indicated that the coarse-wool sheep from Europe were clearly different from those of the Asia. Accordingly, the Asiatic mouflon was positioned between Asian and European countries. In addition, we found that the genetic background of Iranian sheep is present in sheep from China and Kyrgyzstan, as well as India. The revealed admixture patterns of the Iranian sheep and other coarse-wool sheep breeds probably resulted from the expansion of nomads and through the Silk Road trade network.     

Genome-wide analysis in endangered populations: a case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breed’s genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.     

Investigation of the Booroola (FecB) and Inverdale (FecX(I)) mutations in 21 prolific breeds and strains of sheep sampled in 13 countries

Twenty-one of the world's prolific sheep breeds and strains were tested for the presence of the FecB mutation of BMPR1B and the FecX(I) mutation of BMP15. The breeds studied were Romanov (2 strains), Finn (2 strains), East Friesian, Teeswater, Blueface Leicester, Hu, Han, D'Man, Chios, Mountain Sheep (three breeds), German Whiteheaded Mutton, Lleyn, Loa, Galician, Barbados Blackbelly (pure and crossbred) and St. Croix. The FecB mutation was found in two breeds, Hu and Han from China, but not in any of the other breeds. The 12 Hu sheep sampled were all homozygous carriers of FecB (FecB(B)/FecB(B)) whereas the sample of 12 Han sheep included all three genotypes (FecB(B)/FecB(B), FecB(B)/FecB+, FecB+/FecB+) at frequencies of 0.33, 0.58 and 0.08, respectively. There was no evidence of FecX(I) in any of the breeds sampled.     

Genetic diversity and variability in Alpine sheep breeds

The aim of this study was to analyze, through the use of microsatellite molecular markers, the genetic diversity and variability of nine Alpine sheep breeds reared in Italy (Bergamasca, Biellese, Schwarzbraunes Bergschaf, Tiroler Bergschaf, and Schnalserschaf), Germany (Brillenschaf and Weisses Bergschaf) and Slovenia (Bovška and Jezerzkǒ-Solčavska), and one Italian crossbreed population. Allelic richness was rather high in each breed highlighting a considerable genetic diversity. However, the study evidenced a significant departure from Hardy–Weinberg equilibrium in all analyzed breeds caused by a heterozygote deficiency. Such lack seems to be caused both to the presence of population substructure within breed and to a rather high level of inbreeding. The genetic differentiation among breed was rather low (F_ST = 0.057) but significant. The clustering analysis performed with STRUCTURE detected the presence of eight clusters; Schwarzbraunes Bergschaf and Tiroler Bergschaf were grouped together and Biellese and Bergamasca as well, evidencing a high similarity of their genetic make up. Reynolds’ genetic distance estimates confirmed the close relationship between these pairs of breeds. Moreover, Bovška, Jezerzkǒ-Solčavska and Brillenschaf resulted rather similar as expected according to their region of origin and to their common ancestors. Molecular coancestry confirmed as well such findings, suggesting that the studied breeds had origin from different ancestral populations. Concluding, microsatellite resulted a useful tool to investigate breed variability and to characterize Alpine sheep breeds. Obtained findings suggest the need to set up a conservation plan aiming to safeguard and increase the genetic variability of the studied breeds compromised by the high level of inbreeding. Microsatellites genotyping could help to monitor breed variability and to organize matings.     

小尾寒羊微卫星与RAPD标记的研究

小尾寒羊是世界上具有非季节性发情和多胎特性的高繁殖率绵羊品种之一。选择4个位于绵羊6号染色体上且与FecB基因紧密连锁的微卫星标记,对小尾寒羊基因组进行PCR扩增后,采用最小二乘法估计各等位基因片段对产羔数的影响。分析结果表明,等位基因OarJIA-5、OarJIA-10、BM143-12和OarHH55-11可以作为小尾寒羊多胎位点的分子标记。从100条随机引物中筛选出18条引物,对小尾寒羊、大尾寒羊、洼地绵羊、滩羊等4个绵羊品种和鲁北山羊的基因组进行扩增,共扩增出146条带,其中94条表现出多态性,占64．60%,同时扩增出每个品种的特异性条带。采用Nei氏标准距离和NJ聚类分析对不同品种的遗传关系进行分析,结果表明,4个绵羊品种亲缘关系很近,提示它们可能起源于共同的原始祖先。     

Genome‐wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds

Italy counts several sheep breeds, arisen over centuries as a consequence of ancient and recent genetic and demographic events. To finely reconstruct genetic structure and relationships between Italian sheep, 496 subjects from 19 breeds were typed at 50K single nucleotide polymorphism loci. A subset of foreign breeds from the Sheep HapMap dataset was also included in the analyses. Genetic distances (as visualized either in a network or in a multidimensional scaling analysis of identical by state distances) closely reflected geographic proximity between breeds, with a clear north–south gradient, likely because of high levels of past gene flow and admixture all along the peninsula. Sardinian breeds diverged more from other breeds, a probable consequence of the combined effect of ancient sporadic introgression of feral mouflon and long‐lasting genetic isolation from continental sheep populations. The study allowed the detection of previously undocumented episodes of recent introgression (Delle Langhe into the endangered Altamurana breed) as well as signatures of known, or claimed, historical introgression (Merino into Sopravissana and Gentile di Puglia; Bergamasca into Fabrianese, Appenninica and, to a lesser extent, Leccese). Arguments that would question, from a genomic point of view, the current breed classification of Bergamasca and Biellese into two separate breeds are presented. Finally, a role for traditional transhumance practices in shaping the genetic makeup of Alpine sheep breeds is proposed. The study represents the first exhaustive analysis of Italian sheep diversity in an European context, and it bridges the gap in the previous HapMap panel between Western Mediterranean and Swiss breeds.     

Population structure, genetic variation and morphological diversity in indigenous sheep of Ethiopia

We investigated genetic and morphological diversity and population structure of 14 traditional sheep populations originating from four ecological zones in Ethiopia (sub-alpine, wet highland, sub-humid lowland and arid lowland). All animals (n = 672) were genotyped for 17 microsatellite markers and scored for 12 morphological characters. The sheep were initially classified as fat-tailed (11 populations), thin-tailed (one population) and fat-rumped sheep (two populations). These classifications are thought to correspond to three consecutive introduction events of sheep from the Near-East into East Africa. For the 14 populations, allelic richness ranged from 5.87 to 7.51 and expected heterozygosity (H(E)) from 0.66 to 0.75. Genetic differentiations (F(ST) values) between all pairs of populations, except between sub-alpine populations, were significantly different from zero (P < 0.001). Cluster analysis of morphological characters and a dendrogram constructed from genetic distances were broadly consistent with the classification into fat-tailed, thin-tailed and fat-rumped sheep. Bayesian cluster analysis using microsatellite markers indicated that there has been further genetic differentiation after the initial introduction of sheep into Ethiopia. Investigation of factors associated with genetic variation showed that an isolation-by-distance model, independently of other factors, explained most of the observed genetic variation. We also obtained a strong indication of adaptive divergence in morphological characters, patterns of morphological variation being highly associated with ecology even when the effect of neutral genetic divergence (F(ST)) was parcelled out in partial Mantel tests. Using a combination of F(ST) values, Bayesian clustering analysis and morphological divergence, we propose a classification of Ethiopian sheep into six breed groups and nine breeds.