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.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

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

Genetic Diversity and Structure of Two Prominent Zebu Cattle Breeds Adapted to the Arid Region of India Inferred from Microsatellite Polymorphism

This study aims to assess the genetic diversity and population structure of two major zebu dairy breeds (Tharparkar and Rathi) adapted to the arid region of Rajasthan state of India. Various variability estimates indicate the existence of sufficient within-breed genetic diversity. Mean estimates of F-statistics are significantly different from zero: F _IS = 0.112 ± 0.029, F _IT = 0.169 ± 0.033, F _ST = 0.065 ± 0.017. The overall positive value of F _IS (0.112) and an F _IT value (0.169) that is more than the F _ST (0.065) indicate departure from random mating. The drift-based estimates reflect a moderate yet significant level of breed differentiation between the Tharparkar and Rathi breeds. The evaluation of an exact test, showing that allele frequencies across all the loci differed significantly, supports the population differentiation. This is paralleled by the outcome of neighbor-joining clustering based on allele-sharing distance measures. The allocation of a high percentage of individuals (95.7%) to their population of origin and correspondence analysis further substantiates the existence of a cohesive genetic structure in both the breeds.     

Variability of canine microsatellites within and between different dog breeds

Polymorphic animal microsatellites have proved valuable genetic markers. For this project, the variability of 19 canine microsatellite loci was examined within and between three pure breeds of dog: Greyhounds, Labradors, and German Shepherds. The number of alleles, absolute and relative frequencies, and the statistics that express polymorphism within a breed were determined. The evolutionary relationships among these closely related dog breeds were estimated by genetic distance measures developed for use with microsatellite loci. According to the pairwise genetic distances, Greyhounds and German Shepherds had longer diverse evolutionary histories than Greyhounds and Labradors or Labradors and German Shepherds. Although a few breed-specific alleles were observed, the significant differences between breeds are in their relative frequencies and distribution of the alleles across a locus. None of the three pure dog breeds corresponds to Hardy-Weinberg equilibrium. A considerable reduction in intrapopulation variation was observed within three pure breeds, compared with the population of individuals belonging to 15 dog breeds. This reduction was especially pronounced in the Greyhound breed, which expressed the lowest degree of variation. Intrapopulation variations of Labradors and German Shepherds did not differ significantly, that of Labradors being only slightly higher. The intra-species variation of dogs is lower than in humans, mouse, or rat, but similar to that in domestic animals, probably reflecting similarly high inbreeding coefficients. However, some highly informative loci were common to all dog breeds tested so far. Such population data are necessary for mapping studies and linkage analysis in dogs. Received: 31 July 1996 / Accepted: 21 October 1996     

Genetic Structure of an Endangered Portuguese Semiferal Pony Breed, the Garrano

The present study intends to survey the genetic variability of an endangered semiferal Portuguese native pony breed, the Garrano. Thirteen microsatellite markers were examined in 277 animals born in 1998, belonging to eight subpopulations corresponding to eight northern Portuguese geographic regions. Mean heterozygosity (H_o) in the Garrano breed was 0.732, ranging from 0.531 to 0.857 across subpopulations. Allelic frequencies and diversity differed significantly between regions, suggesting the existence of genetic differentiation within the breed confirmed by the population differentiation estimator F_ST. Allele sharing genetic distance (Dps) was used to determine the relationship between the analyzed subpopulations, some of which are diverging significantly from the others. Relationships among six Iberian horse breeds (including the Garrano) were assessed through the restricted maximum likelihood method, which clusters the Garrano with another Celtic pony, the Spanish Asturcon, traditionally bred in geographical continuity with the Garrano. Results reveal a statistically significant deficit of heterozygotes within the Garrano breed (F_IT = 0.031, p < 0.05). Regular monitoring should therefore be implemented so that the effect of genetic drift within subpopulations, enhanced by inbreeding, may be successfully minimized.     

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 structure and admixture in sheep from terminal breeds in the United States

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s F_ST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s F_ST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture , and based on cross-validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.     

Tracking the emergence of a new breed using 49,034 SNP in sheep

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F _ST. The divergence separating the Florida and Louisiana Native (F _ST = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F _ST = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.     

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.     

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

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

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.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

Genome-Wide Genetic Diversity and Differentially Selected Regions among Suffolk,Rambouillet, Columbia,Polypay, and Targhee Sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F _ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F _ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F _ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.     

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.     

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.