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 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.     

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 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.     

Genetic variations in blood proteins within and between Spanish dairy sheep breeds

Eleven flocks of three Spanish sheep breeds have been studied: four Churra, four Lacha and three Manchega. Genetic variations between flocks of the same and different breeds have been analysed, using eight genetic blood systems and a wide range of statistical methods. The results show that the Churra, Lacha and Manchega breeds present significant differences between each other. Significant intraracial variations have been observed in Churras and Lachas, although they were almost non-existent in Manchegas. In some cases the genetic distances between flocks of the same breed were greater than the distances between flocks of different breeds. The data indicate that the closest breeds genetically are Lacha and Manchega.     

Factors influencing the success of vaginal and laparoscopic artificial insemination in churra ewes: a field assay

Pregnancy rate following artificial insemination (AI) in sheep is variable depending on several factors. The Churra breed (milk breed of the North-West of Spain) yields lower fertility results compared with other local and European breeds (Manchega, Latxa, Merino, Lacaune, Sarde, etc.). In this work we studied the influence of many factors on the fertility of the Churra breed (insemination technique, year, farm, age, male, number of inseminations per ewe, lambing-insemination interval and technician), analyzing lambing data obtained after 44448 inseminations (39.67% cervical AI via vagina, AIV, and 60.33% intrauterine AI using laparoscopy, AIL) in a categorical model. The most important factors influencing fertility after AI were farm, year, season, AI technique, and technician. AIL showed significantly higher fertility results than AIV (44.89% versus 31.25%). Season significantly affected fertility in both cases, but differences were more evident in AIV. Fertility dropped 1.74% (AIV) and 2.07% (AIL) per year as the ewes aged. Finally, AI fertility decreased when the lambing-insemination interval was lower than 10 weeks.     

Hemotypes in Spanish sheep breeds

The structure and the genetic diversity of the Churra, Lacha and Manchega sheep breeds have been analysed using hemotypes observed in eight loci. The three breeds are different in their hemotypes in terms both of quantity and quality. The proportions of unique hemotypes in Churra (60%), Lacha (61%) and Manchega (67%) revealed a high level of individual diversity within each breed. Racial genetic diversity follows the descending order of: Manchega-Lacha-Churra. The value of N:H (number of animals: total hemotypes) in the multi-racial population was 3.05.     

Red cell reduced glutathione concentrations in Spanish Churra sheep

Erythrocyte reduced glutathione (GSH) levels were investigated in Spanish Churra sheep. GSH deficiency appeared in a high frequency, a clear bimodal distribution being apparent. No significant concentrations of amino acids were detected in the samples and no significant differences were found in potassium concentrations between the low-GSH and the high-GSH type animals. Such results indicate that erythrocyte GSH deficiency in Churra sheep may be similar to the ‘Merino type’ GSH deficiency. Furthermore, limited inheritance data suggested that a second type of GSH deficiency might be present also in Churra sheep.     

Genetic variation at microsatellite loci in Spanish sheep

Genetic variation at 18 microsatellite loci was analysed in six indigenous Spanish sheep: Churra; Latxa; Manchega; Rasa-Aragonesa; Castellana and Merino. Merinos had frequently the highest number of alleles per locus, whereas Latxas showed the lowest one at many loci. Markers ordered decreasingly according to the number of variants differentiated in the whole population were: MAF70; TGLA13; CSSM66; BM143, BM6444; MAF36; MAF64; CSSM6; TGLA53; OarFCB11; MAF33; BM4621; MAF48; MAF65; BM1258; ILSTS002; ADCYC and OarCP34. Parameters of variability such as effective number of alleles and gene diversities corroborated the high level of variation frequently displayed by microsatellite markers. Comparison of allele distributions among populations and loci did not reveal consistent shapes. Distributions were centralised in some cases, whereas in others some kind of skewness was evident. Breed-specific alleles were detected at most loci, being frequent in Merinos and rare in Churras.     

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.     

Economic weights of fertility, prolificacy, milk yield and longevity in dairy sheep

Economic weights have been estimated in two breeds (Latxa and Manchega) using economic and technical data collected in 41 Latxa and 12 Manchega dairy sheep flocks. The traits considered were fertility (lambing per year), prolificacy (number of lambs), milk yield (litres) and longevity (as productive life, in years). A linear function was used, relating these traits to the different costs in the flock. The variable costs involved in the profit function were feed and labour. From this function, economic weights were obtained. Labour is considered in the Latxa breed to be a constraint. Moreover, farm profits are unusually high, which probably means that some costs were not included according to the economic theory. For that reason, a rescaling procedure was applied constraining total labour time at the farm. Genetic gains were estimated with the resulting economic weights to test if they give any practical difference. Milk yield only as selection criterion was also considered. The medians of the estimated economic weights for fertility, prolificacy, milk yield and longevity were 138.60 € per lambing, 40.00 € per lamb, 1.18 € per l, 1.66 € per year, and 137.66 € per lambing, 34.17 € per lamb, 0.73 € per l, 2.16 € per year under the linear approach in the Latxa and Manchega breeds respectively. Most differences between breeds can be related to differences in production systems. As for the genetic gains, they were very similar for all economic weights, except when only milk yield was considered, where a correlated decrease in fertility led to a strong decrease in profit. It is concluded that the estimates are robust for practical purposes and that breeding programmes should consider inclusion of fertility. More research is needed to include other traits such as somatic cell score, milk composition and udder traits.     

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.     

Reproductive response and LH secretion in ewes treated with melatonin implants and induced to ovulate with the ram effect

Two experiments were conducted to examine the effects of treating seasonally anoestrous ewes with melatonin before ram introduction on reproductive response, and on LH secretion in anoestrous ewes induced to ovulate by rams.In Experiment 1, a total of 667 ewes from three flocks involving Merino (Flock 1, N = 149), Merino entrefino (Flock 2, N = 325) and Rasa Aragonesa (Flock 3, N = 203) breeds were used. Within each flock, ewes isolated from rams since the previous lambing were assigned at random to receive melatonin implants of Regulin (75, 175 and 105 in Merino, Merino entrefino and Rasa Aragonesa flocks, respectively) or to serve as untreated controls (74 in Merino, 150 in Merino entrefino and 98 in Rasa Aragonesa flocks). Fertile rams were introduced into all flocks 5 weeks after implantation in March (Flocks 1 and 2) or April (Flock 3), and remained with the ewes for a 50 day mating period. Percentage of ewes with luteal activity at ram introduction did not differ between melatonin treated and control ewes in any flock. There were no significant differences in either the mean interval from ram introduction to lambing or the distribution of lambing. Implantation with melatonin resulted in an improvement of prolificacy in all three flocks, although this only reached statistical significance in the Merino flock (1.15 vs. 1.03 in treated and control ewes, respectively, P < 0.05). Fertility was increased significantly (P < 0.05) in the Merino entrefino flock (64.5% in treated vs. 51.3% in control ewes).In Experiment 2, two trials were undertaken utilizing a total of 63 ewes. Trial 1 involved 24 mature Manchega ewes and Trial 2 involved 39 Merino ewe lambs. Half of the animals in each trial received a Regulin implant on 28 February (Trial 1) or 12 March (Trial 2) and the remaining half acted as controls. Rams were introduced 5 weeks after implantation and remained with the ewes for a 25 day period. In both trials, anoestrous ewes at ram introduction were bled at 20 min intervals for 3 h before and 5 h after ram introduction and then at 3 h intervals over the next 24 h for assessment of plasma concentrations of LH. Secretion of LH before or following introduction of rams was not affected by melatonin. Both treated and control anoestrous ewes in each trial responded to introduction of rams with an increase in the frequency of the LH pulses (P < 0.05), but no significant changes were detected in pulse amplitude or mean levels of LH. A preovulatory surge of LH was detected between 8 and 26 h after ram introduction, but neither mean interval from ram introduction to the peak of LH surge, nor the magnitude of the LH peak, was influenced by melatonin treatment.Results from this study show that: (1) melatonin implants administered during early seasonal anoestrus have the potential to improve reproductive performance in Spanish breeds of sheep, but the response is conditioned by breed, management system and environmental factors; (2) melatonin did not modify the secretion of LH in anoestrous ewes induced to ovulate by the ram effect under our experimental conditions.     

Genetic diversity in nine native Turkish sheep breeds based on microsatellite analysis

The genetic variability of 204 sheep from nine breeds (Awassi, Çine Çapar?, Karaka?, Karya, Karayaka, Morkaraman, Norduz, Sak?z and Tuj) growing in four different regions (western, northern, eastern and south‐eastern Anatolia) was assessed using 18 microsatellite loci. High mean number of alleles, allelic richness and factorial correspondence analysis showed the degree of admixture between native sheep breeds of Turkey. The Karya and Çine Çapar? breeds were observed as the most distinct of the breeds, and possible introgressions were detected in other breeds. It was found that 8.9% of genetic variation resulted from the difference between the populations. The genetic variation in Turkish breeds was not much higher than that of European breeds, which might be a consequence of the recent sharp decrease in sheep numbers.     

Influence of breed and age on morphometry and depth of inseminating catheter penetration in the ewe cervix: a postmortem study

A detailed examination of the cervical canal in the ewe was carried out. This analysis could be used to design new catheters for artificial insemination (AI) to achieve deeper cervical penetration and therefore better fertility results. Three hundred and sixty-five cervices from four sheep breeds (Churra, Assaf, Merino, Castellana) obtained postmortem were used. Cervix morphometry and depth of cervical penetration using two types of catheters were determined. A conventional straight catheter for ovine artificial insemination (IMV), and a bent catheter, ending in a stainless steel needle, 9 cm in length and with an 8 mm tip bent 45 degrees , were used. The results showed that the morphometry of the cervix depends on breed and age of the ewe. The cervices of Churra breed were shorter and narrower, and had a higher number of folds than those of other breeds. Postmortem cervical penetration was deeper when the cervices were longer and wider, and with fewer folds (Merino and Castellana breeds). In ageing ewes, the cervix tended to become longer and wider, with loose folds. This decreased structural complexity and significantly improved cervical penetration. The bent catheter allowed significantly greater cervical penetration than the straight IMV one.     

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.     

Mitochondrial diversity and the origin of Iberian sheep

Mitochondrial DNA diversity was analysed in 19 Iberian and six foreign sheep breeds. Three mtDNA lineages (B, A and C) were found in the Iberian sheep, with type B clearly predominating over the others. The results were analysed for each of the morphologically determined breed groups in Iberian sheep: Merino, Entrefino, Churro and Iberian trunks. MtDNA lineage C was found only in the Iberian trunk composed of Montesina and Ojalada. These two populations had high mtDNA variability, and in the Iberian sheep only Merino Branco had more variation. The other three Merino types studied showed moderate variability, including the most authentic Merino, the Spanish Merino. These three Merinos clustered closely in a multidimensional scaling representation of distances, while the fourth breed (Merino Branco) showed a clear separation. As for the other two trunks, breeds from the Churro group showed greater maternal uniformity while results for populations included in the so-called Entrefino trunk seemed to have a more heterogeneous maternal origin. The results obtained are discussed with available data from nuclear markers and with morphological classifications, and all this information is analysed in relation to the origin of the different Iberian sheep breeds.     

Sheep blood polymorphism and genetic divergence between French Rambouillet and Spanish Merino: role of genetic drift

Rambouillet sheep originating from Spanish Merino have been maintained in France as a small and closed flock since their importation. After 190 years of independent evolution, the flock has markedly differentiated from its Spanish parental population. The observed differences between them were characterized by the fixation in Rambouillet of the Mb and F30 alleles, which occurred in Spanish Merino with frequencies of 0.90 and 0.80 respectively (at two distinct blood group loci M and F30) and by the absence in Rambouillet of other alleles or phenogroups (at the Tf and the A, B, C blood group loci) which were observed in Spanish Merino with frequencies ranging from 0.10 to 0.28. On the basis of their phenotypic distributions at 11 blood polymorphic loci, the two populations differed significantly from each other (total chi 2 values = 352.62, 23 df, P less than 0.001). By comparing the observed magnitude of gene frequency differences between Rambouillet and Spanish Merino with the estimate of inbreeding coefficient for Rambouillet obtained from pedigrees, it appeared that the observed genetic differences could be attributed to the evolutionary change due to random drift in the small and closed flock of Rambouillet.     

Single-strand conformational polymorphism and sequence polymorphism of Mhc-DRB in Latxa and Karrantzar sheep: implications for Caprinae phylogeny

Single-strand conformational polymorphism analysis and DNA sequencing were used to characterize Mhc-DRB second exon variability in the Latxa and Karrantzar breeds of sheep. The presence of more than two sequences in some animals indicates that alleles of two different loci have been amplified. Six new alleles were identified by sequencing. The allele frequency distribution of the DRB1 gene is striking, with two alleles accounting for half of the gene pool in both breeds under study. The most frequent allele in both breeds was the same (named DRB1*0702), with some specific amino acids: Tyr in position 31 and Thr in 51. A species variability analysis was also performed including the entire set of sheep DRB exon 2 sequences. Based on the patchwork patterns of different alleles, interallelic recombination appears to be playing a significant role in the generation of allelic diversity at this locus in sheep. The phylogenetic tree of all known Caprinae DRB sequences shows that certain alleles from one species are more closely related to those from other species than they are to each other. Allele DRB1*0702 merits special attention due to its high similarity to the Mufflon allele. As this is the most frequent in both breeds analyzed, one can hypothesize that in sheep, both Mufflon and Argali have had different influences depending on the sheep breed under study and that the relationship between domestic sheep and Mufflon is greater than previously thought. The data generated in this study can serve as a basis for developing a typing assay for the sheep DRB genes in the Latxa and Karrantzar populations.