Genetic diversity of wild boar populations and domestic pig breeds (Sus scrofa) in South-western Europe

Genetic characterization of domestic pigs in the core zone of swine production of Argentina

Argentina is a small player in the global pork market, contributing only 0.7% of the total production. With increasing global demand for meat, there is an opportunity for countries with an agricultural profile to grow their pork production. However, there is a need to understand the current state of the pork production sector in all aspects to inform decision-making. The aim of this study was to genetically characterize pig herds from different production strata in the primary region for pork production in the country. For this purpose, phylogenetic and genetic variability analyses were performed using the mitochondrial control region marker (n=95 pig samples). Moreover, genotyping of ryr1 and PRKAG3 genes (n=108 pig samples) were performed to evaluate the frequency of deleterious alleles for meat quality traits in the region. The results showed high levels of genetic variability in the pig herds (Hd= 0.840 ± 0.031 and π= 0.010 ± 0.001), with a creole sow and Iberian lineage standing out in the phylogeny. The genotyping of the ryr1 marker revealed the presence of the deleterious t allele in all analyzed strata. However, the RN-allele of the PRKAG3 gene was detected only in the two lower strata. This study represents the first analysis of the phylogenetic relationships among domestic pigs from Argentina and provides an initial assessment of genetic variability in the region. Additionally, the results present, for the first time, the frequency of deleterious alleles for pig production in the productive core area, demonstrating their prevalence.

Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history

European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.

Introduced, Mixed, and Peripheral: Conservation of Mitochondrial-DNA Lineages in the Wild Boar (Sus scrofa L.) Population in the Urals

Translocations and introductions are important events that allow organisms to overcome natural barriers. The genetic background of colonization success and genetic consequences of the establishment of populations in new environments are of great interest for predicting species’ colonization success. The wild boar has been introduced into many parts of the world. We analyzed sequences of the mitochondrial-DNA control region in the wild boars introduced into the Ural region and compared them with sequences from founder populations (from Europe, the Caucasus, Central Asia, and the Far East). We found that the introduced population has high genetic diversity. Haplotypes from all the major phylogenetic clades were detected in the analyzed group of the animals from the Urals. In this group, no haplotypes identical to Far Eastern sequences were detectable despite a large number of founders from that region. The contribution of lineages originating from Eastern Europe was greater than expected from the proportions (%) of European and Asian animals in the founder populations. This is the first study on the genetic diversity and structure of a wild boar population of mixed origin at the northern periphery of this species’ geographical range.

Patterns of genetic variation on wild pig (Sus scrofa) populations over a complete range of the species in Argentina

In Argentina, wild pigs (Sus scrofa) are represented by wild boars and feral pigs. These populations inhabit a wide territory due to natural dispersion and human translocation. Previous studies in other countries have detected crossbreeding between the different morphotypes, even with domestic pigs. This crossing can lead to introgression of improved traits in the wild population, which would increase the invasive potential and, therefore, the damage to native ecosystems. The aim of this work was to obtain the patterns of genetic variability throughout its current distribution in Argentina, in order to elucidate genetic relationships between wild boar and feral pig populations through the molecular marker control region. For this purpose, studies of genetic variability and population structure were carried out using 377 sequences from Argentinian wild pigs. The high values of haplotype and nucleotide diversity (Hd = 0.866 and π = 0.00959) obtained, and the cluster analyses (SAMOVA and BAPS) could indicate mixing between wild pigs and/or with domestic pigs. The star-like shapes observed in the haplotype network and neutral tests (Fu’s Fs and Tajima’s D) are consistent with a recent population expansion, supporting previous reports that indicate crossbreeding increases invasive potential.

Population Genetic Structure of the Wild Boar (Sus scrofa) in the Carpathian Basin

In the Carpathian Basin the wild boar (Sus scrofa) belongs among the most important game species both ecologically and economically, therefore knowing more about the basics of the genetics of the species is a key factor for accurate and sustainable management of its population. The aim of this study was to estimate the genetic diversity and to elucidate the genetic structure and location of wild boar populations in the Carpathian Basin. A total of 486 samples were collected and genotyped using 13 STR markers. The number of alleles varied between 4 and 14, at 9 of the 13 loci the observed heterozygosity was significantly different (p < 0.05) from the expected value, showing remarkable introgression in the population. The population was separated into two groups, with an F_st value of 0.03, suggesting the presence of two subpopulations. The first group included 147 individuals from the north-eastern part of Hungary, whereas the second group included 339 samples collected west and south of the first group. The two subpopulations’ genetic indices are roughly similar. The lack of physical barriers between the two groups indicates that the genetic difference is most likely caused by the high reproduction rate and large home range of the wild boars, or by some genetic traces’ having been preserved from both the last ice age and the period before the Hungarian water regulation.

Spatial genetic structure and historical demography of East Asian wild boar

Pleistocene climatic fluctuations may have had a profound impact on the evolutionary history of many species. The geographical pattern of European wild boar (Sus scrofa) is clearly studied, and it was greatly influenced by ancient climatic events, especially the Last Glacial Maximum. Previous research on genetic variation has mainly focused on the origin and distribution histories of domestic pigs. However, some questions have not been answered, including those concerning the genetic diversity, geographical pattern and possible historic influence of climate on East Asian wild boar (EAWB). Employing the control region of mtDNA (511 bp), we investigated the contributions of historic climate, which possibly shaped the genetic pattern of wild boar. Given that the level of genetic diversity of wild boars is higher in East Asia than in Europe, 172 haplotypes were detected from 680 individuals. Phylogenetic analysis demonstrated the complex phylogeographic structure of EAWB. Mismatch analysis, neutrality tests and the Bayesian Skyline Plot results all retrieved signals of a rapid population expansion, which might have played an important role in driving the formation of complex spatial genetic structure. Genetic data and species distribution modelling showed that the Last Glacial Maximum had weak effect on the distribution of the EAWB. We suggest that, in shaping spatial genetic structure in East Asian, long-term gene flow and population history played more important roles than Pleistocene climate fluctuations.

Phylogeography and epidemiology of Brucella suis biovar 2 in wildlife and domestic swine

Swine brucellosis due to Brucella suis biovar 2 (bv2) is enzootic in wild boar and hare in continental Europe and may cause major economic losses to the pig industry, mainly in free-ranged pig farms. The high nucleotide identity found among the B. suis biovar 2 isolates has long hindered the full understanding of the epidemiology and the phylogeography of the disease. Here, we used multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) and whole-genome analysis to identify single-nucleotide polymorphisms (SNPs) in order to gain insights from the largest B. suis bv2 dataset analyzed so far composed of domestic pigs and wildlife isolates collected throughout Europe since the 1970s. We found four major clades with a specific phylogeographic pattern. The Iberian clade contains isolates exclusively from the Iberian Peninsula. The Central European clade includes most isolates from France, Northern Italy, Switzerland and an important proportion of those of Northern Spain. The Eastern European clade clustered isolates from Croatia and Hungary mainly but also from areas of France, Germany, Italy and Poland. Finally, a separated Sardinian clade grouped three isolates from this island. At fine scale, MLVA demonstrated an endemic status of the infection in Europe and it allowed tracking a large outbreak formed by different farms from Spain linked to the same infection source. The whole genome SNP analysis showed that the strains form genetically distinct clades, shared between wild boar and pigs, in agreement with the MLVA clades. Interestingly, all hare isolates clustered together within two groups composed exclusively of wildlife isolates. Our results support the hypothesis that maintenance and spread of B. suis bv2 in Europe is a dynamic process linked to the natural expansion of wild boar as the main wild reservoir of the infection, while spread over long distances is found largely dependent on anthropogenic activities.

Bronze Age meat industry: ancient mitochondrial DNA analyses of pig bones from the prehistoric salt mines of Hallstatt (Austria)

Objective

In the Bronze Age Hallstatt metropolis (‘Salzkammergut’ region, Upper Austria), salt richness enabled the preservation of pork meat to sustain people’s livelihood suggesting an organized meat production industry on a yearly basis of hundreds of pigs. To pattern the geographic and temporal framework of the early management of pig populations in the surrounding areas of Hallstatt, we want to gain insights into the phylogeographic network based on DNA sequence variation among modern pigs, wild boars and prehistoric (likely) domestic pigs.

Results

In this pilot study, we successfully adapted ancient DNA extraction and sequencing approaches for the analysis of mitochondrial DNA sequence variation in ten prehistoric porcine teeth specimens. Minimum-spanning network analyses revealed unique mitochondrial control region DNA haplotypes ranging within the variation of modern domestic pig and wild boar lineages and even shared haplotypes between prehistoric and modern domestic pigs and wild boars were observed.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3340-7) contains supplementary material, which is available to authorized users.

Demographic history, current expansion and future management challenges of wild boar populations in the Balkans and Europe

Wild boar (Sus scrofa), one of the most widespread wildlife species, has entered a stage of continuous growth in Europe, and could even be considered a pest species. We analysed microsatellite variability in 723 wild boars from across Europe, including the northern Dinaric Balkans. Our aims were: (1) to define the population structure of wild boars in the Balkans and its relation with other European populations; (2) to estimate effective populations sizes, levels of intra- and inter-population diversity, inbreeding migration and gene flow patterns; (3) to test subpopulations for bottlenecks; (4) to interpret these results in light of current knowledge about the demographic history of wild boars in Europe; and (5) to discuss the relevance of these findings for management and conservation. Strong population structuring was observed and 14 subpopulations were revealed. High genetic diversity was found, and besides the well-known identity of the Italian populations of Sardinia and Castelporziano, we bring new insights into other potential relevant, refugial populations such as Littoral Slovenia, South Portugal, North-western Iberia and an entire cluster in the Balkans. There was evidence of gene flow going from these refugial subpopulations towards less peripheral and more admixed subpopulations. Recent population bottlenecks and expansions were detected, mostly in the peninsular refuge subpopulations. The results are consistent with the fluctuations of wild boar numbers in Europe since the beginning of the twentieth century. These results should be taken into account in future conservation and management plans for wild boar populations in Europe.

Contact Zone of Asian and European Wild Boar at North West of Iran

Wild boar (Sus scrofa) are widely distributed throughout the Old World. Most studies have focused on Europe and East Asia with the genetic diversity of West Asia being less well studied. In particular, the genetic variability and genetic structure of the Iranian populations are not yet known; gaps which prevent scientists from resolving the genetic relationships of the Eurasian wild boar. This paper is the first attempt to provide information about genetic relationships among modern Iranian populations of the Eurasian wild boar (S. scrofa) by sequencing 572 bp of the mitochondrial (mt) DNA control region. As a result of this investigation, it was discovered that Iran contains not only Middle Eastern haplotypes, but also shares haplotypes with Europe and East Asia. The Italian clade, which is endemic in Italy, is not identified in Iran, while all other clades, including Asiatic, European, Near East 1, and Near East 2 are found based on the phylogenetic tree and median-joining network. The results of this study illustrate that north west of Iran (specifically Southwest Caspian Sea) is the contact zone between the Asian (Near Eastern and Far Eastern), and the European clades. In light of the fact that the domestication of pigs occurs in Anatolia, this finding is important.

Hybrid origin of European commercial pigs examined by an in-depth haplotype analysis on chromosome 1

Although all farm animals have an original source of domestication, a large variety of modern breeds exist that are phenotypically highly distinct from the ancestral wild population. This phenomenon can be the result of artificial selection or gene flow from other sources into the domesticated population. The Eurasian wild boar (Sus scrofa) has been domesticated at least twice in two geographically distinct regions during the Neolithic revolution when hunting shifted to farming. Prior to the establishment of the commercial European pig breeds we know today, some 200 years ago Chinese pigs were imported into Europe to improve local European pigs. Commercial European domesticated pigs are genetically more diverse than European wild boars, although historically the latter represents the source population for domestication. In this study we examine the cause of the higher diversity within the genomes of European commercial pigs compared to their wild ancestors by testing two different hypotheses. In the first hypothesis we consider that European commercial pigs are a mix of different European wild populations as a result of movement throughout Europe, hereby acquiring haplotypes from all over the European continent. As an alternative hypothesis, we examine whether the introgression of Asian haplotypes into European breeds during the Industrial Revolution caused the observed increase in diversity. By using re-sequence data for chromosome 1 of 136 pigs and wild boars, we show that an Asian introgression of about 20% into the genome of European commercial pigs explains the majority of the increase in genetic diversity. These findings confirm that the Asian hybridization, that was used to improve production traits of local breeds, left its signature in the genome of the commercial pigs we know today.

Untangling the hybrid nature of modern pig genomes: a mosaic derived from biogeographically distinct and highly divergent Sus scrofa populations

The merging of populations after an extended period of isolation and divergence is a common phenomenon, in natural settings as well as due to human interference. Individuals with such hybrid origins contain genomes that essentially form a mosaic of different histories and demographies. Pigs are an excellent model species to study hybridization because European and Asian wild boars diverged ~1.2 Mya, and pigs were domesticated independently in Europe and Asia. During the Industrial Revolution in England, pigs were imported from China to improve the local pigs. This study utilizes the latest genomics tools to identify the origin of haplotypes in European domesticated pigs that are descendant from Asian and European populations. Our results reveal fine-scale haplotype structure representing different ancient demographic events, as well as a mosaic composition of those distinct histories due to recently introgressed haplotypes in the pig genome. As a consequence, nucleotide diversity in the genome of European domesticated pigs is higher when at least one haplotype of Asian origin is present, and haplotype length correlates negatively with recombination frequency and nucleotide diversity. Another consequence is that the inference of past effective population size is influenced by the background of the haplotypes in an individual, but we demonstrate that by careful sorting based on the origin of haplotypes, both distinct demographic histories can be reconstructed. Future detailed mapping of the genomic distribution of variation will enable a targeted approach to increase genetic diversity of captive and wild populations, thus facilitating conservation efforts in the near future.

Contemporary genetic structure, phylogeography and past demographic processes of wild boar Sus scrofa population in Central and Eastern Europe

The wild boar (Sus scrofa) is one of the most widely distributed mammals in Europe. Its demography was affected by various events in the past and today populations are increasing throughout Europe. We examined genetic diversity, structure and population dynamics of wild boar in Central and Eastern Europe. MtDNA control region (664 bp) was sequenced in 254 wild boar from six countries (Poland, Hungary, Belarus, Ukraine, Moldova and the European part of Russia). We detected 16 haplotypes, all known from previous studies in Europe; 14 of them belonged to European 1 (E1) clade, including 13 haplotypes from E1-C and one from E1-A lineages. Two haplotypes belonged respectively to the East Asian and the Near Eastern clade. Both haplotypes were found in Russia and most probably originated from the documented translocations of wild boar. The studied populations showed moderate haplotype (0.714±0.023) and low nucleotide diversity (0.003±0.002). SAMOVA grouped the genetic structuring of Central and Eastern European wild boar into three subpopulations, comprising of: (1) north-eastern Belarus and the European part of Russia, (2) Poland, Ukraine, Moldova and most of Belarus, and (3) Hungary. The multimodal mismatch distribution, Fu's Fs index, Bayesian skyline plot and the high occurrence of shared haplotypes among populations did not suggest strong demographic fluctuations in wild boar numbers in the Holocene and pre-Holocene times. This study showed relatively weak genetic diversity and structure in Central and Eastern European wild boar populations and underlined gaps in our knowledge on the role of southern refugia and demographic processes shaping genetic diversity of wild boar in this part of Europe.

Conservation genomic analysis of domestic and wild pig populations from the Iberian Peninsula

BACKGROUND

Inbreeding is among the major concerns in management of local livestock populations. The effective population size of these populations tends to be small, which enhances the risk of fitness reduction and extinction. High-density SNP data make it possible to undertake novel approaches in conservation genetics of endangered breeds and wild populations.A total of 97 representative samples of domestic and wild pig populations from the Iberian Peninsula, subjected to different levels of threat with extinction, were genotyped with a 60 K SNP panel. Data analyses based on: (i) allele frequency differences; (ii) linkage disequilibrium and (iii) runs of homozygosity were integrated to study population relationships, inbreeding and demographic history.

RESULTS

The domestic pigs analyzed belonged to local Spanish and Portuguese breeds: Iberian ─ including the variants Retinto Iberian, Negro Iberian and Manchado de Jabugo ─, Bisaro and Chato Murciano. The population structure and persistence of phase analysis suggested high genetic relations between Iberian variants, with recent crossbreeding of Manchado de Jabugo with other pig populations. Chato Murciano showed a high frequency of long runs of homozygosity indicating recent inbreeding and reflecting the recent bottleneck reported by historical records. The Chato Murciano and the Manchado de Jabugo breeds presented the lowest effective population sizes in accordance with their status of highly inbred breeds. The Iberian wild boar presented a high frequency of short runs of homozygosity indicating past small population size but no signs of recent inbreeding. The Iberian breed showed higher genetic similarities with Iberian wild boar than the other domestic breeds.

CONCLUSIONS

High-density SNP data provided a consistent overview of population structure, demographic history and inbreeding of minority breeds and wild pig populations from the Iberian Peninsula. Despite the very different background of the populations used, we found a good agreement between the different analyses. Our results are also in agreement with historical reports and provide insight in the events that shaped the current genetic variation of pig populations from the Iberian Peninsula. The results exposed will aid to design and implement strategies for the future management of endangered minority pig breeds and wild populations.

How immunogenetically different are domestic pigs from wild boars: a perspective from single-nucleotide polymorphisms of 19 immunity-related candidate genes

The coexistence of wild boars and domestic pigs across Eurasia makes it feasible to conduct comparative genetic or genomic analyses for addressing how genetically different a domestic species is from its wild ancestor. To test whether there are differences in patterns of genetic variability between wild and domestic pigs at immunity-related genes and to detect outlier loci putatively under selection that may underlie differences in immune responses, here we analyzed 54 single-nucleotide polymorphisms (SNPs) of 19 immunity-related candidate genes on 11 autosomes in three pairs of wild boar and domestic pig populations from China, Iberian Peninsula, and Hungary. Our results showed no statistically significant differences in allele frequency and heterozygosity across SNPs between three pairs of wild and domestic populations. This observation was more likely due to the widespread and long-lasting gene flow between wild boars and domestic pigs across Eurasia. In addition, we detected eight coding SNPs from six genes as outliers being under selection consistently by three outlier tests (BayeScan2.1, FDIST2, and Arlequin3.5). Among four non-synonymous outlier SNPs, one from TLR4 gene was identified as being subject to positive (diversifying) selection and three each from CD36, IFNW1, and IL1B genes were suggested as under balancing selection. All of these four non-synonymous variants were predicted as being benign by PolyPhen-2. Our results were supported by other independent lines of evidence for positive selection or balancing selection acting on these four immune genes (CD36, IFNW1, IL1B, and TLR4). Our study showed an example applying a candidate gene approach to identify functionally important mutations (i.e., outlier loci) in wild and domestic pigs for subsequent functional experiments.

Introgression and isolation contributed to the development of Hungarian Mangalica pigs from a particular European ancient bloodline

Background

Mangalica breeds are indigenous to Hungary and their breeding history dates back to about 200–250 years ago. They are fat-type pigs and have a rare curly hair phenotype. The aim of our study was to establish the relationships between these unique breeds and other European breeds.

Results

Based on a core sequence of 382 bp present in 2713 mitochondrial D-loop sequences from pigs belonging to 38 local breeds from nine countries, five cosmopolitan breeds and wild boars from 14 countries, we identified 164 haplotypes. More than half of the 2713 sequences belonged to either four haplotypes characteristic of continental European breeds or two haplotypes characteristic of British/cosmopolitan breeds; each haplotype is present in more than 100 individuals. Most Mangalica individuals belonged either to one of these common continental European haplotypes or to two Mangalica-specific haplotypes that were absent in all other breeds. In addition, we identified the ancestral mitochondrial D-loop signature present in these 2713 sequences and found that ~ 80% carried the European ancient signatures, ANC-Aside and ANC-Cside or their closely related signatures, while most of the remaining sequences carried a modern Asian signature, ANC-Easia. Mangalica individuals carried the ANC-Aside signature, but not the ANC-Cside or ANC-Easia signatures.

Conclusions

In all the Mangalica individuals, a unique ancient European signature was found in the mitochondrial DNA D-loop region, but they belonged almost exclusively to either certain very abundant European or two Mangalica-specific D-loop haplotypes. This indicates that the present-day Mangalica population in Hungary evolved either by introgression of other European breeds and wild boars or via total isolation after the divergence of European ancient porcine bloodlines.

Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands

Background

Native pig breeds in the Iberian Peninsula are broadly classified as belonging to either the Celtic or the Mediterranean breed groups, but there are other local populations that do not fit into any of these groups. Most of the native pig breeds in Iberia are in danger of extinction, and the assessment of their genetic diversity and population structure, relationships and possible admixture between breeds, and the appraisal of conservation alternatives are crucial to adopt appropriate management strategies.

Methods

A panel of 24 microsatellite markers was used to genotype 844 animals representing the 17 most important native swine breeds and wild populations existing in Portugal and Spain and various statistical tools were applied to analyze the results.

Results

Genetic diversity was high in the breeds studied, with an overall mean of 13.6 alleles per locus and an average expected heterozygosity of 0.80. Signs of genetic bottlenecks were observed in breeds with a small census size, and population substructure was present in some of the breeds with larger census sizes. Variability among breeds accounted for about 20% of the total genetic diversity, and was explained mostly by differences among the Celtic, Mediterranean and Basque breed groups, rather than by differences between domestic and wild pigs. Breeds clustered closely according to group, and proximity was detected between wild pigs and the Mediterranean cluster of breeds. Most breeds had their own structure and identity, with very little evidence of admixture, except for the Retinto and Entrepelado varieties of the Mediterranean group, which are very similar. Genetic influence of the identified breed clusters extends beyond the specific geographical areas across borders throughout the Iberian Peninsula, with a very sharp transition from one breed group to another. Analysis of conservation priorities confirms that the ranking of a breed for conservation depends on the emphasis placed on its contribution to the between- and within-breed components of genetic diversity.

Conclusions

Native pig breeds in Iberia reveal high levels of genetic diversity, a solid breed structure and a clear organization in well-defined clusters.

Domestication does not narrow MHC diversity in Sus scrofa

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.

Genome-wide single nucleotide polymorphism analysis reveals recent genetic introgression from domestic pigs into Northwest European wild boar populations

Present-day genetic introgression from domestic pigs into European wild boar has been suggested in various studies. However, no hybrids have been identified beyond doubt mainly because available methods were unable to quantify the extent of introgression and rule out natural processes. Genetic introgression from domestic pigs may have far-reaching ecological consequences by altering traits like the reproduction rate or immunology of wild boar. In this study, we demonstrate a novel approach to investigate genetic introgression in a Northwest (NW) European wild boar data set using a genome-wide single nucleotide polymorphism (SNP) assay developed for domestic pigs. We quantified the extent of introgression using allele frequency spectrum analysis, in silico hybridization simulations and genome distribution patterns of introgressed SNPs. Levels of recent introgression in the study area were expected to be low, as pig farming practices are prevailingly intensive and indoors. However, evidence was found for geographically widespread presence of domestic pig SNPs in 10% of analysed wild boar. This was supported by the identification of two different pig mitochondrial DNA haplotypes in three of the identified hybrid wild boar, suggesting that introgression had occurred from multiple sources (pig breeds). In silico hybridization simulations showed that the level of introgression in the identified hybrid wild boar is equivalent to first-generation hybrids until fifth-generation backcrosses with wild boar. The distribution pattern of introgressed SNPs supported these assignments in four of nine hybrids. The other five hybrids are considered advanced-generation hybrids, resulting from interbreeding among hybrid individuals. Three of nine hybrids were genetically associated with a different wild boar population than the one in which they were sampled. This discrepancy suggests that genetic introgression has occurred through the escape or release of an already hybridized farmed wild boar stock. We conclude that genetic introgression from domestic pigs into NW European wild boar populations is more recent and more common than expected and that genome-wide SNP analysis is a promising tool to quantify recent hybridization in free-living populations.