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.     

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

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

Genetic distances within and across cattle breeds as indicated by biallelic AFLP markers

We tested the use of biallelic Amplified fragment length polymorphism (AFLP) polymorphisms for the estimation of relative genetic distances of cattle individuals within or across breeds. An allele permutation procedure was developed to estimate the stochastic variation of the genetic distance that is inherent to a given dataset. In a panel of 47 Holstein-Friesian cattle analysed with 248 polymorphic markers, the average genetic distance of bulls selected for breeding was slightly lower than the distance of the cows. The observed standard deviation (SD) of the distance indicated genetic subdivision, which for the bulls was explained by variation in the additive relationship derived from herdbook data. Animals from three different breeds, the highly selected Holstein-Friesian, the Italian Brown and the historic Maremmana, were compared on the basis of 106 polymorphic markers. No breed-specific fragments were observed. The mean pair-wise genetic distance within breeds was 85% of the value across breeds, but principal coordinates analysis clustered the animals according to their breed of origin. Calculation of distances between the breeds indicated a relatively divergent position of the Maremmana, relative to the two other breeds. However, biallelic markers indicate that the process of breed formation had only a limited effect on the diversity at marker loci.     

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.     

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.     

Mitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds

This survey represents the first characterization of mitochondrial DNA diversity within three breeds of Indian sheep (two strains of the Deccani breed, as well as the Bannur and Garole breeds) from different geographic regions and with divergent phenotypic characteristics. A 1061-bp fragment of the mitochondrial genome spanning the control region, a portion of the 12S rRNA gene and the complete phenyl tRNA gene, was sequenced from 73 animals and compared with the corresponding published sequence from European and Asian breeds and the European Mouflon (Ovis musimon). Analysis of all 156 sequences revealed 73 haplotypes, 52 of which belonged to the Indian breeds. The three Indian breeds had no haplotypes in common, but one Indian haplotype was shared with European and other Asian breeds. The highest nucleotide and haplotype diversity was observed in the Bannur breed (0.00355 and 0.981 respectively), while the minimum was in the Sangamneri strain of the Deccani breed (0.00167 and 0.882 respectively). All 52 Indian haplotypes belonged to mitochondrial lineage A. Therefore, these Indian sheep are distinct from other Asian and European breeds studied so far. The relationships among the haplotypes showed strong breed structure and almost no introgression among these Indian breeds, consistent with Indian sheep husbandry, which discourages genetic exchange between breeds. These results have implications for the conservation of India's ovine biodiversity and suggest a common origin for the breeds investigated.     

Population structure in Indian sheep ascertained using microsatellite information

This study attempts to provide a comprehensive insight into the prevailing genetic status of Indian sheep breeds using microsatellite markers. Seventeen Indian sheep breeds from 3 agroecological zones were analysed using a panel of 25 microsatellite markers. All of the sheep breeds investigated were genetically diverse, as evident from the high allele (>6) and gene (>0.6) diversity values. The gene diversity values for all breeds ranged from 0.621 to 0.780. The within-population heterozygote deficit (F(IS)) varied from -0.098 to 0.234, reflecting significant levels for 12 of the 17 breeds investigated. The average genetic differentiation between all breeds (F(ST)) was 11.1%, revealing moderate discrimination between the indigenous sheep breeds. The genetic distance and principal component analysis revealed a separation of sheep breeds based on geographical propinquity. The Bayesian clustering approach suggested poor breed differentiation in the north-western arid and semi-arid region when compared to the breeds from the eastern and southern peninsular regions. The observed results mirror the divergent management strategies in the different agroecological regions, lack of specific selection policies, and intermixing of breeds in close proximity. Immediate steps to curb the intermixing and erosion of breed purity for some of these breeds need to be implemented, for example, by introducing measures like making proven rams available and ensuring their frequent exchange between flocks. The data generated here provides valuable information about the genetic structure of the 17 Indian sheep breeds and this can be used for designating priorities for their conservation.     

Runs of homozygosity reveal genome‐wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium‐density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH = 0.099), whereas Comisana showed the lowest (F_ROH = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.     

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.     

Microsatellite Analysis Characterizes Burkina Faso as a Genetic Contact Zone Between Sahelian and Djallonké Sheep

A total of 123 sheep belonging to the Djallonké, Mossi, and Burkina-Sahel breeds, along with 41 Spanish Xalda sheep were genotyped for 27 microsatellites. The pair Djallonké-Mossi had the highest between breeds molecular coancestry. Admixture analysis informed on the parental role of the Burkina-Sahel and Djallonké breeds. The Mossi breed was a hybrid population nearer to the Djallonké breed. Only half of the Mossi individuals were correctly assigned to their breed. The Burkina-Sahel and Djallonké breeds can be considered ancestrally different genetic entities. Differentiation between the Djallonké and Mossi breeds may be due to introgression of Sahelian sheep.     

Marker-derived phylogeny of European cattle supports demic expansion of agriculture

Using frequencies of 86 genes from 23 loci of blood group systems, blood and milk proteins, the genetic relationships among 14 cattle breeds including four native Balkan and four synthetic Balkan-Alp breeds were studied. The dendrogram and nonlinear map construction shows a consensus ‘Balkan breed cluster’, an ‘Alp breed cluster’, an unstable position of synthetic breeds and well-separated American breeds. Positive partial correlations between genetic distance and time elapsed since introduction of farming while keeping geographical distances constant, and regular patterns over thousands of kilometres indicate that large-scale cattle population movements together with human migration (in the Neolithic age) from the Near East into Europe across the Balkans are the most likely explanation for the genetic distances observed in our data. More recent breed differentiation and selection do not yet blur this initial pattern of European cattle populations.     

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.     

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 relationships among European cattle breeds

Genetic relationships among 37 European cattle breeds were investigated using blood group and serum protein polymorphisms. The 18 859 animals included in the study represented a random sample from pedigree populations in the UK. Within-breed variation was estimated by average heterozygosity and number of alleles observed, and breed relationships were evaluated by genetic distance. Standard errors of the heterozygosity, number of alleles and genetic distance were obtained by bootstrapping. The significance of breed differences was tested using an exact test of differentiation. French, Italian and Channel Island breeds were found to have generally higher heterozygosities and a greater number of alleles than breeds from mainland Britain and North Europe. Genetic distances ranged between 0·011 (±0·005) and 0·309 (±0·071). Two major breed groups were identified; a group of French, Italian and Channel Island breeds together with the Simmental and Gelbvieh, and a second group consisting of the mainland British and North European breeds. The exact test of breed differentiation showed all breeds to be significantly different from one another ( P < 0·0001). Overall relationships among breeds reflected their geographical origin and common ancestry rather than the agricultural use for which the breeds have been selected.     

Genetic diversity and admixture analysis of Sanfratellano and three other Italian horse breeds assessed by microsatellite markers

Sanfratellano is a native Sicilian horse breed, mainly reared in the north east of the Island, developed in the 19th century from local dams and sires with a restricted introgression of Oriental, African and, more recently, Maremmano stallions. In this study, the genetic relationships and admixture among Sanfratellano, the other two Sicilian autochthonous breeds and Maremmano breed were assessed using a set of microsatellites. The main goals were to infer the impact of Maremmano breed in the current Sanfratellano horse and to provide genetic information useful to improve the selection strategies of the Sanfratellano horse. The whole sample included 384 horses (238 Sanfratellano, 50 Sicilian Oriental Purebred, 30 Sicilian Indigenous and 66 Maremmano), chosen avoiding closely related animals. A total of 111 alleles from 11 microsatellite loci were detected, from four at HTG7 to 15 at ASB2 locus. The mean number of alleles was the lowest in Oriental Purebred (6.7), the highest in Sanfratellano (8.3). All the breeds showed a high level of gene diversity (He) ranging from 0.71 ± 0.04 in Sicilian Oriental Purebred to 0.81 ± 0.02 in Sicilian Indigenous. The genetic differentiation index was low; only about 6% of the diversity was found among breeds. Nei's standards (DS) and Reynolds' (DR) genetic distances reproduced the same population ranking. Individual genetic distances and admixture analysis revealed that: (a) nowadays Maremmano breed does not significantly influence the current Sanfratellano breed; (b) within Sanfratellano breed, it is possible to distinguish two well-defined groups with different proportions of Indigenous blood.     

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.     

How do introgression events shape the partitioning of diversity among breeds: a case study in sheep

Background

From domestication to the current pattern of differentiation, domestic species have been influenced by reticulate evolution with multiple events of migration, introgression, and isolation; this has resulted in a very large number of breeds. In order to manage these breeds and their genetic diversity, one must know the current genetic structure of the populations and the relationships among these. This paper presents the results of a genetic diversity analysis on an almost exhaustive sample of the sheep breeds reared in France. Molecular characterization was performed with a set of 21 microsatellite markers on a collection of 49 breeds that include five breed types: meat, hardy meat, dairy, high prolificacy and patrimonial breeds.

Results

Values of expected heterozygosity ranged from 0.48 to 0.76 depending on the breed, with specialized meat breeds exhibiting the lowest values. Neighbor-Net, multidimensional analysis or clustering approaches revealed a clear differentiation of the meat breeds compared to the other breed types. Moreover, the group that clustered meat breeds included all the breeds that originated from the United Kingdom (UK) and those that originated from crossbreeding between UK breeds and French local breeds. We also highlighted old genetic introgression events that were related to the diffusion of Merino rams to improve wool production. As a result of these introgression events, especially that regarding the UK breeds, the breeds that were clustered in the ‘meat type cluster’ exhibited the lowest contribution to total diversity. That means that similar allelic combinations could be observed in different breeds of this group.

Conclusions

The genetic differentiation pattern of the sheep breeds reared in France results from a combination of factors, i.e. geographical origin, historic gene flow, and breed use. The Merino influence is weaker than that of UK breeds, which is consistent with how sheep use changed radically at the end of 19^th century when wool-producing animals (Merino-like) were replaced by meat-producing breeds. These results are highly relevant to monitor and manage the genetic diversity of sheep and can be used to set priorities in conservation programs when needed.

Electronic supplementary material

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

Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers

Very little research into genetic diversity of Italian native dog breeds has been carried out so far. In this study we aimed to estimate and compare the genetic diversity of four native Italian shepherd dog breeds: the Maremma, Bergamasco, Lupino del Gigante and Oropa shepherds. Therefore, some cosmopolitan dog breeds, which have been widely raised in Italy for a long time past, have also been considered to check possible influence of these dog populations on the Italian autochthonous breeds considered here. A total of 212 individuals, belonging to 10 different dog breeds, were sampled and genotyped using 18 autosomal microsatellite loci. We analyzed the genetic diversity of these breeds, within breed diversity, breed relationship and population structure. The 10 breeds considered in this study were clearly genetically differentiated from each other, regardless of current population sizes and the onset of separate breeding history. The level of genetic diversity explained 20% of the total genetic variation. The level of H_E found here is in agreement with that found by other studies. The native Italian breeds showed generally higher genetic diversity compared with the long established, well-defined cosmopolitan dog breeds. As the Border Collie seems closer to the Italian breeds than the other cosmopolitan shepherd dogs considered here, a possible utilization of this breed to improve working performance in Italian traditional working shepherd dogs cannot be ignored. The data and information found here can be utilized in the organization of conservation programs planned to reduce inbreeding and to minimize loss of genetic variability.     

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.