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.     

Genetic diversity measures of local European beef cattle breeds for conservation purposes

This study was undertaken to determine the genetic structure, evolutionary relationships, and the genetic diversity among 18 local cattle breeds from Spain, Portugal, and France using 16 microsatellites. Heterozygosities, estimates of Fst, genetic distances, multivariate and diversity analyses, and assignment tests were performed. Heterozygosities ranged from 0.54 in the Pirenaica breed to 0.72 in the Barrosã breed. Seven percent of the total genetic variability can be attributed to differences among breeds (mean F_st = 0.07; P < 0.01). Five different genetic distances were computed and compared with no correlation found to be significantly different from 0 between distances based on the effective size of the population and those which use the size of the alleles. The Weitzman recursive approach and a multivariate analysis were used to measure the contribution of the breeds diversity. The Weitzman approach suggests that the most important breeds to be preserved are those grouped into two clusters: the cluster formed by the Mirandesa and Alistana breeds and that of the Sayaguesa and Tudanca breeds. The hypothetical extinction of one of those clusters represents a 17% loss of diversity. A correspondence analysis not only distinguished four breed groups but also confirmed results of previous studies classifying the important breeds contributing to diversity. In addition, the variation between breeds was sufficiently high so as to allow individuals to be assigned to their breed of origin with a probability of 99% for simulated samples.     

Genetic characterization of Uruguayan Pampa Rocha pigs with microsatellite markers

In this study, we genetically characterized the Uruguayan pig breed Pampa Rocha. Genetic variability was assessed by analyzing a panel of 25 microsatellite markers from a sample of 39 individuals. Pampa Rocha pigs showed high genetic variability with observed and expected heterozygosities of 0.583 and 0.603, respectively. The mean number of alleles was 5.72. Twenty-four markers were polymorphic, with 95.8% of them in Hardy Weinberg equilibrium. The level of endogamy was low (F_IS = 0.0475). A factorial analysis of correspondence was used to assess the genetic differences between Pampa Rocha and other pig breeds; genetic distances were calculated, and a tree was designed to reflect the distance matrix. Individuals were also allocated into clusters. This analysis showed that the Pampa Rocha breed was separated from the other breeds along the first and second axes. The neighbour-joining tree generated by the genetic distances D_A showed clustering of Pampa Rocha with the Meishan breed. The allocation of individuals to clusters showed a clear separation of Pampa Rocha pigs. These results provide insights into the genetic variability of Pampa Rocha pigs and indicate that this breed is a well-defined genetic entity.     

Genetic diversity of eleven European pig breeds

A set of eleven pig breeds originating from six European countries, and including a small sample of wild pigs, was chosen for this study of genetic diversity. Diversity was evaluated on the basis of 18 microsatellite markers typed over a total of 483 DNA samples collected. Average breed heterozygosity varied from 0.35 to 0.60. Genotypic frequencies generally agreed with Hardy-Weinberg expectations, apart from the German Landrace and Schwäbisch-Hällisches breeds, which showed significantly reduced heterozygosity. Breed differentiation was significant as shown by the high among-breed fixation index (overall F_ST = 0.27), and confirmed by the clustering based on the genetic distances between individuals, which grouped essentially all individuals in 11 clusters corresponding to the 11 breeds. The genetic distances between breeds were first used to construct phylogenetic trees. The trees indicated that a genetic drift model might explain the divergence of the two German breeds, but no reliable phylogeny could be inferred among the remaining breeds. The same distances were also used to measure the global diversity of the set of breeds considered, and to evaluate the marginal loss of diversity attached to each breed. In that respect, the French Basque breed appeared to be the most "unique" in the set considered. This study, which remains to be extended to a larger set of European breeds, indicates that using genetic distances between breeds of farm animals in a classical taxonomic approach may not give clear resolution, but points to their usefulness in a prospective evaluation of diversity.     

Genetic diversity and differentiation in Portuguese cattle breeds using microsatellites

Genotype data from 30 microsatellites were used to assess genetic diversity and relationships among 10 native Portuguese cattle breeds, American Charolais and the Brazilian Caracú. Hardy–Weinberg equilibrium was observed for all loci/population combinations except for five loci in Brava de Lide and one locus in Alentejana that exhibited heterozygote deficiency. Estimates of average observed and expected heterozygosities, total number of alleles (TNA) per breed and mean number of alleles (MNA) per locus/population were obtained. A total of 390 alleles were detected. TNA among Iberian cattle ranged from 170 to 237 and MNA ranged from 5.67 to 8.07. The highest observed heterozygosities were found in the Caracú, Maronesa, Garvonesa and Arouquesa and the lowest in Brava de Lide and Mirandesa. Estimation of population subdivision using Wright's F_ST index showed that the average proportion of genetic variation explained by breed differences was 9%. Neighbour‐joining phylogenetic trees based on D_A distances showed that the genetic relationships of present‐day Portuguese native breeds are consistent with historical origins in the Brown Concave (Arouquesa, Mirandesa, Marinhoa) and Red Convex (Mertolenga, Alentejana, Garvonesa, Minhota) evolutionary groups. The Iberian Black Orthoide group, represented by Brava de Lide and Maronesa, and the Barrosã breed appeared to be more closely related to the Brown Concave group but may represent a separate lineage. The Caracú breed was not found to be closely associated with any of the native Portuguese breeds.     

Genetic diversity and population structure of 10 Chinese indigenous egg-type duck breeds assessed by microsatellite polymorphism

The genetic structure and diversity of 10 Chinese indigenous egg-type duck breeds were investigated using 29 microsatellite markers. The total number of animals examined were 569, on average 57 animals per breed were selected. The microsatellite marker set analysed provided 177 alleles (mean 6.1 alleles per locus, ranging from 3 to 10). All populations showed high levels of heterozygosity with the lowest estimate of 0.539 for the Jinding ducks, and the highest 0.609 observed for Jingjiang partridge ducks. The global heterozygote deficit across all populations (F _IT) amounted to −0.363. About 10% of the total genetic variability originated from differences among breeds, with all loci contributing significantly. An unrooted consensus tree was constructed using the NeighborNet tree based on the Reynold’s genetic distance. The structure software was used to assess genetic clustering of these egg-type duck breeds. Clustering analysis provided an accurate representation of the current genetic relations among the breeds. An integrated analysis was undertaken to obtain information on the population dynamics in Chinese indigenous egg-type duck breeds, and to better determine the conservation priorities.     

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.     

Multilocus genotypic data reveal high genetic diversity and low population genetic structure of Iranian indigenous sheep

Iranian livestock diversity is still largely unexplored, in spite of the interest in the populations historically reared in this country located near the Fertile Crescent, a major livestock domestication centre. In this investigation, the genetic diversity and differentiation of 10 Iranian indigenous fat‐tailed sheep breeds were investigated using 18 microsatellite markers. Iranian breeds were found to host a high level of diversity. This conclusion is substantiated by the large number of alleles observed across loci (average 13.83, range 7–22) and by the high within‐breed expected heterozygosity (average 0.75, range 0.72–0.76). Iranian sheep have a low level of genetic differentiation, as indicated by the analysis of molecular variance, which allocated a very small proportion (1.67%) of total variation to the between‐population component, and by the small fixation index (F_ST = 0.02). Both Bayesian clustering and principal coordinates analysis revealed the absence of a detectable genetic structure. Also, no isolation by distance was observed through comparison of genetic and geographical distances. In spite of high within‐breed variation, signatures of inbreeding were detected by the F_IS indices, which were positive in all and statistically significant in three breeds. Possible factors explaining the patterns observed, such as considerable gene flow and inbreeding probably due to anthropogenic activities in the light of population management and conservation programmes, are discussed.     

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F _st values against genome location. The highest F _st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.     

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.     

Genetic variation and relationships of eighteen Chinese indigenous pig breeds

Chinese indigenous pig breeds are recognized as an invaluable component of the world''s pig genetic resources and are divided traditionally into six types. Twenty-six microsatellite markers recommended by the FAO (Food and Agriculture Organization) and ISAG (International Society of Animal Genetics) were employed to analyze the genetic diversity of 18 Chinese indigenous pig breeds with 1001 individuals representing five types, and three commercial breeds with 184 individuals. The observed heterozygosity, unbiased expected heterozygosity and the observed and effective number of alleles were used to estimate the genetic variation of each indigenous breed. The unbiased expected heterozygosity ranged between 0.700 (Mashen) and 0.876 (Guanling), which implies that there is an abundant genetic variation stored in Chinese indigenous pig breeds. Breed differentiation was shown by fixation indices (F_IT, F_IS, and F_ST). The F_ST per locus varied from 0.019 (S0090) to 0.170 (SW951), and the average F_ST of all loci was 0.077, which means that most of the genetic variation was kept within breeds and only a little of the genetic variation exists between populations. The Neighbor-Joining tree was constructed based on the Nei D_A (1978) distances and one large cluster with all local breeds but the Mashen breed, was obtained. Four smaller sub-clusters were also found, which included two to four breeds each. These results, however, did not completely agree with the traditional type of classification. A Neighbor-Joining dendrogram of individuals was established from the distance of – ln(proportions of shared alleles); 92.14% of the individuals were clustered with their own breeds, which implies that this method is useful for breed demarcation. This extensive research on pig genetic diversity in China indicates that these 18 Chinese indigenous breeds may have one common ancestor, helps us to better understand the relative distinctiveness of pig genetic resources, and will assist in developing a national plan for the conservation and utilization of Chinese indigenous pig breeds.     

Genomic diversity and population structure of Mexican and Spanish bovine Lidia breed

The Lidia bovine breed is distinguished for its low genetic exchangeability given its selection on aggressive behavior, its management uniqueness and its subdivided structure. In this study, we present a comprehensive genome‐wide analysis of genetic diversity, population structure and admixture of 468 animals from Mexican and Spanish Lidia breed populations and 64 samples belonging to 10 Spanish native and American‐creole breeds using 37 148 single nucleotide polymorphisms. We found similar average inbreeding values in the Lidia breed, with different distributions within groups; variability of inbreeding values among Spanish lineages was significant and no differences were found among the Mexican sub‐populations. Together, the high F_IS of the lineages and the behavior of the runs of homozygosity are consequences of the lineage's small effective population sizes, contributing to their inbreeding increase. Population admixture analysis discarded any influence on the genetic structure of the Lidia populations from the Spanish native and American‐creole breeds. In addition, both Lidia populations depicted different genetic origins, with the exception of some Mexican individuals whose origins traced back to recent Spanish importations.     

Population genetic structure of the gynodioecious Thymus vulgaris L. (Labiatae) in southern France

In a self-compatible gynodioecious species, the abundance of female plants (which are obligate outcrossers) relative to hermaphrodites (which may self and outcross) may be a critical factor influencing genetic diversity and population structure. In the gynodioecious Thymus vulgaris L., female frequency varies from 5 to 95%, providing a suitable model to examine this issue. In this study, we use allozyme markers to (1) evaluate the relationship between female frequency, genetic diversity and population structure, (2) determine whether females and hermaphrodites vary in heterozygote deficiency and (3) examine whether other factors such as plant density are related to heterozygote deficiency. Twenty three natural populations, with female frequencies ranging from 11 to 92%, were sampled in and around the St-Martin-de-Londres basin in southern France. Based on four polymorphic allozyme loci, we found no significant correlation between female frequency and heterozygote deficiency. A significant (P < 0.05) F_IS value over loci and over populations of 0.11 was detected. The F_IS value per population showed a significant heterozygote deficiency in 11 of the 23 populations. However, no significant difference between female and hermaphrodite F_IS values was found. A significant heterozygote deficiency only occurred in populations of intermediate density. There was little differentiation among populations (F_ST = 0.038) nor among subpopulations within each population. The significant F_IS values are thus mostly due to inbreeding effects. The lack of a correlation between F_IS values and female frequency may be due to outcrossing in hermaphrodites and/or restoration of male fertility which may occur to a greater extent at low female frequency. The similarity of female and hermaphrodite F_IS values indicates that females may occasion high levels of biparental inbreeding.     

Population genetics of <Emphasis Type="Italic">Collisella subrugosa</Emphasis> (Patellogastropoda: Acmaeidae): evidence of two scales of population structure

Marine invertebrate populations usually show high levels of genetic variability that has frequently been associated with spatial and temporal environmental heterogeneity. One of the most heterogeneous marine environments is the intertidal zone, the habitat of Collisella subrugosa, the most widespread and abundant Brazilian limpet. C. subrugosa has planktonic larvae that can disperse over long distances, what can promote gene flow among shores, working against interpopulational differentiation. In this study we investigated the genetic variability and populational substructure of C. subrugosa through analysis of 24 allozyme loci in 14 samples (590 individuals) collected along 2,700 km of the Brazilian coast. The genetic variability was high ( and ), as expected for intertidal species. Genetic differentiation among samples was low (F _ST = 0.03) what may reflect intensive gene flow associated with larval dispersal. However, we detected an isolation-by-distance pattern of population substructure in one sampled region. High levels of heterozygote deficiency were also observed for many loci in each sample. Alternative hypothesis are discussed, and the “breeding groups” is suggested to explain these pattern, indicating the main cause as environmental heterogeneity.     

Genetic profile of the locally developed Meatmaster sheep breed in South Africa based on microsatellite analysis

This report is on the development and genetic distinctiveness of the Meatmaster, a locally developed South African composite sheep breed. Genotypes of 175 Meatmaster sheep were compared to genotypes of 242 individuals from six putative parent breeds and an outgroup. All comparisons were based on 10 microsatellite loci. No significant breed-specific alleles were found for the Meatmaster or most of the parental breeds. Nevertheless, a Bayesian-based assignment test placed 69.4% of Meatmaster individuals in clusters distinct from the well-defined breed-specific clusters of the parent breeds. Population differentiation from R_ST values showed low differentiation among the four Meatmaster populations (0.004–0.028), higher values for pair-wise combinations of Meatmaster populations and the parent breeds (0.051–0.194), and the highest values among Meatmaster populations and the outgroup: Namakwa Afrikaner (0.237–0.301). Clustering in a dendrogram placed all Meatmaster populations in distinct clusters with strong bootstrap support (97%). Trends from a hierarchical analysis of total variation confirmed some distinctiveness in the Meatmaster, with 2.315% of total variation found among Meatmaster populations compared to 11.455% among breeds. High levels of heterozygosity and numbers of alleles in Meatmaster populations indicated good retention of genetic diversity during the potential bottleneck created by selection in its early history. Results are discussed with reference to genetic distinctiveness in this developing breed, the genetic contribution of the Damara as a common ancestor, conservation of genetic diversity and an application for future genetic management for the Meatmaster breed.     

Genome‐wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.     

Sequence characterization of the melanocortin 1 receptor (MC1R) gene in sheep with different coat colours and identification of the putative e allele at the ovine Extension locus

Sequence of the melanocortin 1 receptor (MC1R) gene (the Extension locus) was obtained from a panel of 73 animals belonging to 9 Italian sheep breeds or populations (Appenninica, Bergamasca, Comisana, Cornigliese-like, Delle Langhe, Massese, Merinizzata Italiana, Sarda and Valle del Belice) with different coat colours. Evaluation of the identified polymorphisms on this phenotype was reported with in silico predictions and comparative approaches within and across breeds and across species. Five novel single nucleotide polymorphisms (SNPs), organized in three haplotypes, were detected. Another haplotype, including the two missense mutations already described for the E^D allele, was identified in few Massese sheep. One SNP (c.199C > T) caused a predicted amino acid substitution (p.R67C) in a highly conserved position of the first intracellular loop of the MC1R protein. The same substitution causes recessive pheomelanism in other species. We propose that the p.67C allele represents the recessive e allele at the ovine Extension series that was, so far, not completely recognized in sheep by classical genetic studies. This polymorphism was analysed in a total of 388 sheep of the 9 investigated breeds. The p.67C allele was identified only in the Valle del Belice breed (allele frequency of 21.3% in 176 analysed animals of this breed) in which the presence of epistatic white-determining loci might mask, at least in part, its effects. Confirming the effect of this novel allele on coat colour will lead to new perspectives on the composition of specialized coloured sheep lines.     

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 lineages reveal intense gene flow between Iberian wild boars and South Iberian pig breeds

The phylogeography of wild boars (WB) and domestic pigs (Sus scrofa) has contributed important insights into where and when domestication occurred. The geographic distribution of two core haplotypes (E1a and E1c) of the main European phylogenetic clade suggests that Central Europe was an early domestication centre, although the complexity of the pattern does not exclude the possibility that multiple domestication events occurred in different regions. To investigate the relationships among WB and domestic pig breeds in Iberia, a fragment of the mitochondrial DNA control region from a large sample (n = 409) of WB and local pig breeds was co‐analysed with published sequences from other European populations. The Iberian sample revealed a high frequency of a sub‐cluster (E1c) of the European haplogroup E1 in 77% of total Iberian samples, 96% of WB, 90% of Alentejano (Portugal) and 87% of Iberian breed pigs (Spain; Black Hairy, Black Hairless and Red varieties). Low genetic distance (F’_ST = 0.105) was observed between Alentejano (Portugal) and Iberian breed pigs (Spain). Alentejano and Iberian breed pigs showed low genetic distances to both Iberian and Central European WB (average F’_ST = 0.345 and 0.215, respectively). This pattern suggests that early pig husbandry in the Iberian Peninsula did not solely rely on imported Central European stock, but also included the recruitment of local WB.