Differentiation among Spanish sheep breeds using microsatellites

Genetic variability at 18 microsatellites was analysed on the basis of individual genotypes in five Spanish breeds of sheep – Churra, Latxa, Castellana, Rasa-Aragonesa and Merino -, with Awassi also being studied as a reference breed. The degree of population subdivision calculated between Spanish breeds from F_ST diversity indices was around 7% of total variability. A high degree of reliability was obtained for individual-breed assignment from the 18 loci by using different approaches among which the Bayesian method provided to be the most efficient, with an accuracy for nine microsatellites of over 99%. Analysis of the Bayesian assignment criterion illustrated the divergence between any one breed and the others, which was highest for Awassi sheep, while no great differences were evident among the Spanish breeds. Relationships between individuals were analysed from the proportion of shared alleles. The resulting dendrogram showed a remarkable breed structure, with the highest level of clustering among members of the Spanish breeds in Latxa and the lowest in Merino sheep, the latter breed exhibiting a peculiar pattern of clustering, with animals grouped into several closely set nodes. Analysis of individual genotypes provided valuable information for understanding intra- and inter-population genetic differences and allowed for a discussion with previously reported results using populations as taxonomic units.     

Genetic variation within the Merino sheep breed: analysis of closely related populations using microsatellites

Genetic relationships among six populations of Merino sheep were investigated using microsatellites. The history of the six populations is relatively well documented, with all being derived from the Spanish Merino breed within the last 400 years. Genetic variation was highest amongst the Spanish and Portuguese populations, although the preservation of genetic diversity within the other populations was high. By a variety of different statistical tests the French Mutton, German Mutton and New Zealand Merino populations could be differentiated from each other and the Iberian Merinos, indicating that microsatellites are able to track relatively recent changes in the population structure of sheep breeds. The dendrograms constructed on the basis of microsatellite allelic frequencies showed that populations that have shared selection criteria (meat vs. wool) tend to cluster together.     

Merino and Merino-derived sheep breeds: a genome-wide intercontinental study

Background

Merino and Merino-derived sheep breeds have been widely distributed across the world, both as purebred and admixed populations. They represent an economically and historically important genetic resource which over time has been used as the basis for the development of new breeds. In order to examine the genetic influence of Merino in the context of a global collection of domestic sheep breeds, we analyzed genotype data that were obtained with the OvineSNP50 BeadChip (Illumina) for 671 individuals from 37 populations, including a subset of breeds from the Sheep HapMap dataset.

Results

Based on a multi-dimensional scaling analysis, we highlighted four main clusters in this dataset, which corresponded to wild sheep, mouflon, primitive North European breeds and modern sheep (including Merino), respectively. The neighbor-network analysis further differentiated North-European and Mediterranean domestic breeds, with subclusters of Merino and Merino-derived breeds, other Spanish breeds and other Italian breeds. Model-based clustering, migration analysis and haplotype sharing indicated that genetic exchange occurred between archaic populations and also that a more recent Merino-mediated gene flow to several Merino-derived populations around the world took place. The close relationship between Spanish Merino and other Spanish breeds was consistent with an Iberian origin for the Merino breed, with possible earlier contributions from other Mediterranean stocks. The Merino populations from Australia, New Zealand and China were clearly separated from their European ancestors. We observed a genetic substructuring in the Spanish Merino population, which reflects recent herd management practices.

Conclusions

Our data suggest that intensive gene flow, founder effects and geographic isolation are the main factors that determined the genetic makeup of current Merino and Merino-derived breeds. To explain how the current Merino and Merino-derived breeds were obtained, we propose a scenario that includes several consecutive migrations of sheep populations that may serve as working hypotheses for subsequent studies.

Electronic supplementary material

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

Genetic relationships among Spanish sheep using microsatellites

Five indigenous Spanish breeds of sheep, Churra, Latxa, Manchega, Rasa–Aragonesa and Merino, with Awassi sheep as a reference breed were genotyped for 19 DNA microsatellites. Allele frequencies and mean heterozygosities revealed the greatest genetic variation in Merino sheep and the lowest in Awassis. Differences in variability were not great in the other breeds studied. The dendrograms obtained based on genetic distances showed a large differentiation between Awassi sheep and the Spanish breeds, as was to be expected from their distinct genetic origin. Merinos appeared separated from the other four breeds, of which, according to a classification based on the fleece characteristics, Churra and Latxa belong to the churro type and Manchega and Rasa–Aragonesa to the so called entrefino type, though no clear separation was evident between the two types. These results suggest that morphological data alone are insufficient for determining relationships between breeds and that studies involving genetic markers may be of great assistance.     

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.     

Mitochondrial D-loop sequence variation among Italian horse breeds

The genetic variability of the mitochondrial D-loop DNA sequence in seven horse breeds bred in Italy (Giara, Haflinger, Italian trotter, Lipizzan, Maremmano, Thoroughbred and Sarcidano) was analysed. Five unrelated horses were chosen in each breed and twenty-two haplotypes were identified. The sequences obtained were aligned and compared with a reference sequence and with 27 mtDNA D-loop sequences selected in the GenBank database, representing Spanish, Portuguese, North African, wild horses and an Equus asinus sequence as the outgroup. Kimura two-parameter distances were calculated and a cluster analysis using the Neighbour-joining method was performed to obtain phylogenetic trees among breeds bred in Italy and among Italian and foreign breeds. The cluster analysis indicates that all the breeds but Giara are divided in the two trees, and no clear relationships were revealed between Italian populations and the other breeds. These results could be interpreted as showing the mixed origin of breeds bred in Italy and probably indicate the presence of many ancient maternal lineages with high diversity in mtDNA sequences.     

Genetic diversity in the Churra tensina and Churra lebrijana endangered Spanish sheep breeds and relationship with other Churra group breeds and Spanish mouflon

The aim of the present study was to estimate the genetic intra-breed variability of Churra tensina and Churra lebrijana endangered breeds and to establish genetic relationships with Churra, Latxa and Merino breeds, as well as Spanish mouflon, by using 28 microsatellite markers, to provide useful information for their conservation. Allele frequencies and heterozygosity revealed high genetic variation in the two endangered breeds despite their small population size. Estimates of inbreeding coefficient (F_IS) were significant for all breeds studied, except for Churra lebrijana breed. The highest inbreeding coefficient (F_IS = 0.143) was found in the Spanish mouflon. Genetic differentiation tests (F_ST = 0.121) and assignment of individuals to populations indicated the existence of defined breed populations, and low genetic flow between these breeds. The highest pairwise Reynolds distance (D_R) values were observed between Mouflon and the domestic sheep breeds. Considering only domestic sheep breeds, the Churra lebrijana breed showed the highest pairwise D_R values_. The lowest values were found between Latxa and the other domestic sheep, except for Churra lebrijana. Results of pairwise D_R values, as well as phylogenetic tree and bottleneck analysis showed an important genetic isolation of the Churra lebrijana breed from the other Churra types, and genetic signatures of a demographic bottleneck. Finally, structure analysis of populations detected a population subdivision in the Latxa sheep breed. In conclusion, this study presents valuable insight into the existing genetic variability of two Spanish endangered breeds, as well as the first study in Spanish mouflon based on microsatellite analysis. The high degree of variability demonstrated in Churra tensina and Churra lebrijana implies that these populations are rich reservoirs of genetic diversity.     

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.     

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.     

Evidence of three maternal lineages in Near Eastern sheep supporting multiple domestication events

The variability of mtDNA was analysed in local sheep breeds reared throughout Turkey, for which a fragment of the D-loop region and the complete cytochrome b were sequenced. Phylogenetic analyses performed independently for the D-loop and the Cyt b gene revealed three clearly separated clusters indicating three major maternal lineages, two of which had been previously described as types B and A. The new type, C, was present in all the breeds analysed and showed considerable mtDNA variability. Divergence time was obtained on the basis of Cyt b gene and was estimated to be around 160,000-170,000 years ago for lineages B and A, whereas the divergence of lineage C proved to have occurred earlier (between 450,000 and 750,000 years ago). These times greatly predate domestication and suggest that the origin of modern sheep breeds was more complex than previously thought and that at least three independent sheep domestication events occurred. Our results, together with archaeological information and the current wild sheep populations in the Near East region support the high importance of this area in the sheep domestication process. Finally, the evidence of a third maternal lineage has important implications regarding the history of modern sheep.     

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.     

The origins of Iberian horses assessed via mitochondrial DNA

Despite a number of recent studies that have focused on the origin of domestic horses, genetic relationships between major geographical clusters still remain poorly understood. In this study we analyzed a 296 bp mtDNA fragment from the HVI region of 171 horses representing 11 native Iberian, Barb, and Exmoor breeds to assess the maternal phylogeography of Iberian horses. The mtDNA haplogroup with a CCG motif (nucleotide position 15,494 to 15,496) was the most frequent in Iberian and Barb breeds (0.42 and 0.57, respectively), regardless of geographic location or group of breeds. This finding supports the close genetic relationship between the ancestral mare populations of the Iberian Peninsula and Northern Africa. Phenotypic differences among the Northern and Southern Iberian groups of breeds are not explained by population subdivision based on maternal lineages. Our results also suggest that Northern Iberian ponies--which are phenotypically close to British ponies, especially Exmoor--are the result of an introgression rather than population replacement.     

Genetic and phylogenetic studies of Chinese native sheep breeds (Ovis aries) based on mtDNA D-loop sequences

To determine the genetic diversity and the origin of Chinese sheep, we analyzed 83 complete sequences of mtDNA D-loop from nine Chinese sheep breeds and a foreign breed, together with nine sheep and cattle available sequences from GenBank. The length of the sequences was considerably variable between 1103 and 1225 bp. The hapolotype diversity was 92.7%. The nucleotide diversity was 3.058%. And the mean nucleotide composition of the 83 sequences was 32.9% A, 29.8% T, 22.9% C and 14.4% G, respectively. The NJ phylogenetic tree (the number of replications of bootstrap test is 1000) revealed that there were three distinct major domestic sheep lineages (termed as lineages A–C) in the 10 breeds. The result indicated that Chinese native sheep breeds derive from three different maternal sources. The mismatch distribution analysis showed that the Fs values were −25.15, −12.28, −8.60 for the lineages A–C, respectively (P < 0.01), which suggested that atleast one population expansion events occur in the demographic history of Chinese sheep breeds.     

Genetic signatures of a Mediterranean influence in Iberian Peninsula sheep husbandry

Highly adaptable and versatile populations of domestic sheep, the result of millennia of intense husbandry, are found in almost every corner of the world. Here we describe a genetic survey of sheep from the western fringe of its European distribution. We studied the mitochondrial DNA control region sequences from 161 individuals belonging to 7 Portuguese sheep breeds. Our study revealed a high level of genetic diversity, with an average breed haplotype diversity of 0.983, substantially above that observed in central European breeds, as well as the presence of maternal lineages until now only found in the Middle East and Asia. A broad north-south pattern describes the most important trend in the Portuguese sheep population with a southern population clearly distinct from most other breeds. A recurrent influx of new genetic diversity, probably via the Mediterranean Sea, may explain these patterns and appears to corroborate the importance of this maritime route in the history of both mankind and livestock. Zooarchaeological studies of sheep bones from southern Portugal indicate a marked size increase during the Moslem period that may reflect an improvement of this animal--perhaps part of the well known "Arab agricultural revolution" in Andalusia. This could have been a time when the gene pool of Iberian sheep was substantially enriched and may help to explain the history of modern sheep breeds in this peninsula.     

Impact of strong selection for the PrP major gene on genetic variability of four French sheep breeds (Open Access publication)

Effective selection on the PrP gene has been implemented since October 2001 in all French sheep breeds. After four years, the ARR "resistant" allele frequency increased by about 35% in young males. The aim of this study was to evaluate the impact of this strong selection on genetic variability. It is focussed on four French sheep breeds and based on the comparison of two groups of 94 animals within each breed: the first group of animals was born before the selection began, and the second, 3–4 years later. Genetic variability was assessed using genealogical and molecular data (29 microsatellite markers). The expected loss of genetic variability on the PrP gene was confirmed. Moreover, among the five markers located in the PrP region, only the three closest ones were affected. The evolution of the number of alleles, heterozygote deficiency within population, expected heterozygosity and the Reynolds distances agreed with the criteria from pedigree and pointed out that neutral genetic variability was not much affected. This trend depended on breed, i.e. on their initial states (population size, PrP frequencies) and on the selection strategies for improving scrapie resistance while carrying out selection for production traits.     

On-station comparative analysis of reproductive and survival performance between Red Maasai,Dorper, and Merino sheep breeds

The reproductive performance of ewes and the survivability of lambs to weaning have a critical economic impact on sheep farming worldwide. Further, knowledge of major mortality causes allows an opportunity for improved flock management to evade financial losses. The maximum likelihood estimates for generalised linear mixed models and chi-square test methods were used to examine 971 mating records, 839 and 763 lambs born and weaned (singles or twins) from the Naivasha Sheep and Goats station in Kenya for the years 2011 to 2020 consisting of Dorper, Red Maasai (RedM), and Merino breeds. The RedM (P < 0.05) outperformed Dorper and Merino in weaning rate, whereas reproductive performance between the three breeds was not significantly different (P > 0.05) in litter size and multiple lambings per ewe lambing. On the one hand, Dorper significantly (P < 0.05) outperformed the other two breeds only in weaning weight per lamb born. In addition, among all the major causes of death, pneumonia appeared to be the one to which Dorper breeds were most susceptible (chi-square test, P < 0.05). According to the findings of this study, neither the Dorper nor the Merino sheep breeds were reproductively superior to the RedM in an extensive semi-arid production environment. In addition, Dorper's susceptibility to the leading causes of mortality, particularly pneumonia and sheep pox, were relatively high compared to other breeds and could be a precursor to massive economic losses for Dorper sheep producers. In contrast to the indigenous RedM breed, imported sheep breeds appeared to be more susceptible to major mortality-related under an extensive production system. Therefore, regardless of weaning weight, RedM breed production appears to be a more viable investment for small-scale farmers, particularly in semi-arid regions.     

Conservation priorities for Ethiopian sheep breeds combining threat status,breed merits and contributions to genetic diversity

Prioritizing livestock breeds for conservation needs to incorporate both genetic and non-genetic aspects important for the survival of the breeds. Here, we apply a maximum-utility-strategy to prioritize 14 traditional Ethiopian sheep breeds based on their threat status, contributions to farmer livelihoods (current breed merits) and contributions to genetic diversity. Contributions of the breeds to genetic diversity were quantified using Eding''s marker-estimated kinship approaches. Non-genetic aspects included threats (e.g. low population size, low preferences by farmers) and current merits (economic, ecological and cultural merits). Threat analysis identified eight of the 14 breeds as threatened. Analysis of current merits showed that sub-alpine and arid-lowland breeds contribute most to farmer livelihoods in comparison to other breeds. The highest contribution to the genetic diversity conserved was from the Simien breed. Simien showed high between-breed (low between-breed kinship = 0.04) as well as high within-breed diversity (low within-breed kinship = 0.09 and high H_E = 0.73 and allelic richness = 6.83). We combined the results on threat status, current breed merits and contributions to genetic diversity to produce a ranking of the 14 breeds for conservation purposes. Our results balance the trade-offs between conserving breeds as insurance against future uncertainties and current sustainable utilization. The ranking of breeds provides a basis for conservation strategies for Ethiopian sheep and contributes to a regional or global conservation plan.     

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.     

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

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