A high throughput genotyping approach reveals distinctive autosomal genetic signatures for European and Near Eastern wild boar

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.

Resilience to Historical Human Manipulations in the Genomic Variation of Italian Wild Boar Populations

Human activities can globally modify natural ecosystems determining ecological, demographic and range perturbations for several animal species. These changes can jeopardize native gene pools in different ways, leading either to genetic homogenization, or conversely, to the split into genetically divergent demes. In the past decades, most European wild boar (Sus scrofa) populations were heavily managed by humans. Anthropic manipulations have strongly affected also Italian populations through heavy hunting, translocations and reintroductions that might have deeply modified their original gene pools. In this study, exploiting the availability of the well-mapped porcine genome, we applied genomic tools to explore genome-wide variability in Italian wild boar populations, investigate their genetic structure and detect signatures of possible introgression from domestic pigs and non-native wild boar. Genomic data from 134 wild boar sampled in six areas of peninsular Italy and in Sardinia were gathered using the Illumina Porcine SNP60 BeadChip (60k Single Nucleotide Polymorphisms – SNPs) and compared with reference genotypes from European specimens and from domestic pigs (both commercial and Italian local breeds), using multivariate and maximum-likelihood approaches. Pairwise FST values, multivariate analysis and assignment procedures indicated that Italian populations were highly differentiated from all the other analyzed European wild boar populations. Overall, a lower heterozygosity was found in the Italian population than in the other European regions. The most diverging populations in Castelporziano Presidential Estate and Maremma Regional Park can be the result of long-lasting isolation, reduced population size and genetic drift. Conversely, an unexpected similarity was found among Apennine populations, even at high distances. Signatures of introgression from both non-Italian wild boar and domestic breeds were very limited. To summarize, we successfully applied genome-wide procedures to explore, for the first time, the genomic diversity of Italian wild boar, demonstrating that they represent a strongly heterogeneous assemblage of demes with different demographic and manipulation histories. Nonetheless, our results suggest that a native component of genomic variation is predominant over exogenous ones in most populations.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East-Caucasian and Central Italian populations

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East-Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.

Ancient pigs reveal a near-complete genomic turnover following their introduction to Europe

Archaeological evidence indicates that domestic pigs arrived in Europe, alongside farmers from the Near East ∼8,500 y ago, yet mitochondrial genomes of modern European pigs are derived from European wild boars. To address this conundrum, we obtained mitochondrial and nuclear data from modern and ancient Near Eastern and European pigs. Our analyses indicate that, aside from a coat color gene, most Near Eastern ancestry in the genomes of European domestic pigs disappeared over 3,000 y as a result of interbreeding with local wild boars. This implies that pigs were not domesticated independently in Europe, yet the first 2,500 y of human-mediated selection applied by Near Eastern Neolithic farmers played little role in the development of modern European pigs.

Archaeological evidence indicates that pig domestication had begun by ∼10,500 y before the present (BP) in the Near East, and mitochondrial DNA (mtDNA) suggests that pigs arrived in Europe alongside farmers ∼8,500 y BP. A few thousand years after the introduction of Near Eastern pigs into Europe, however, their characteristic mtDNA signature disappeared and was replaced by haplotypes associated with European wild boars. This turnover could be accounted for by substantial gene flow from local European wild boars, although it is also possible that European wild boars were domesticated independently without any genetic contribution from the Near East. To test these hypotheses, we obtained mtDNA sequences from 2,099 modern and ancient pig samples and 63 nuclear ancient genomes from Near Eastern and European pigs. Our analyses revealed that European domestic pigs dating from 7,100 to 6,000 y BP possessed both Near Eastern and European nuclear ancestry, while later pigs possessed no more than 4% Near Eastern ancestry, indicating that gene flow from European wild boars resulted in a near-complete disappearance of Near East ancestry. In addition, we demonstrate that a variant at a locus encoding black coat color likely originated in the Near East and persisted in European pigs. Altogether, our results indicate that while pigs were not independently domesticated in Europe, the vast majority of human-mediated selection over the past 5,000 y focused on the genomic fraction derived from the European wild boars, and not on the fraction that was selected by early Neolithic farmers over the first 2,500 y of the domestication process.

Assessing Sus scrofa diversity among continental United States, and Pacific islands populations using molecular markers from a gene banks collection

Human migration and trade facilitated domesticated livestock movement, gene flow and development of diverse populations upon which agriculture is based. In addition, varying USA ecological conditions has led to a diverse set of livestock populations to utilize. Quantifying genetic diversity of these populations is incomplete. This paper quantifies genetic diversity captured by the National Animal Germplasm Program and explores genetic structure and differences among 19 pig populations (feral populations from Pacific islands, continental US, and Chinese breeds) using 70,231 SNP from 500 animal samples. Among continental US breeds F_is was consistently low suggesting genetic variability is sufficiently available for breeders to use. A unique population structure using principal component analysis illustrated clear distinctions between Duroc, Yorkshire, Hampshire, breeds of Chinese origin, and feral Pacific Island populations were identified. Five Y chromosome haplotypes were evaluated and demonstrated migration patterns from European, central Asia, and potentially Polynesian waves of gene flow. Quantifying diversity and potential origin of Pacific populations provides insight for future uses, and the need for preservation. Viewing gene bank holdings in context of diversity measures we found a lack of inbreeding within breeds, suggesting the collection represents a wide sampling of individual breeds.

Introgression of Chinese haplotypes contributed to the improvement of Danish Duroc pigs

The distribution of Asian ancestry in the genome of Danish Duroc pigs was investigated using whole-genome sequencing data from European wild boars, Danish Duroc, Chinese Meishan and Bamaxiang pigs. Asian haplotypes deriving from Meishan and Bamaxiang occur widely across the genome. Signatures of selection on Asian haplotypes are common in the genome, but few of these haplotypes have been fixed. By defining 50-kb windows with more than 50% Chinese ancestry, which did not exhibit extreme genetic differentiation between Meishan and Bamaxiang as candidate regions, the enrichment of quantitative trait loci in candidate regions supports that Asian haplotypes under selection play an important role in contributing genetic variation underlying production, reproduction, meat and carcass, and exterior traits. Gene annotation of regions with the highest proportion of Chinese ancestry revealed genes of biological interest, such as NR6A1. Further haplotype clustering analysis suggested that a haplotype of Chinese origin around the NR6A1 gene was introduced to Europe and then underwent a selective sweep in European pigs. Besides, functional genes in candidate regions, such as AHR and PGRMC2, associated with fertility, and SAL1, associated with meat quality, were identified. Our results demonstrate the contribution of Asian haplotypes to the genomes of European pigs. Findings herein facilitate further genomic studies such as genomewide association study and genomic prediction by providing ancestry information of variants.

Population admixture in Chinese and European Sus scrofa

Relationships between different populations were investigated using Porcine 60 K data from 1,135 domestic pigs and wild boars across Europe and China. The results indicate that most European breeds have been introgressed with Chinese ancestry, but the extent of introgression varies considerably among breeds. Moreover, the main source of this introgression is pigs from South China, closely related to Bamaxiang and Dongshan pigs. Contributions from East and Central Chinese pig breeds are also detectable. Phylogeny reconstruction places European wild boars among European domestic breeds. Coalescent simulations indicate that this may be the result of gene flow from European wild boars to European domestic pigs. These results will facilitate further genomic studies such as genome-wide association studies, selection signature detection and genomic prediction.

Origin and dispersal of early domestic pigs in northern China

It is widely accepted that modern pigs were domesticated independently at least twice, and Chinese native pigs are deemed as direct descendants of the first domesticated pigs in the corresponding domestication centers. By analyzing mitochondrial DNA sequences of an extensive sample set spanning 10,000 years, we find that the earliest pigs from the middle Yellow River region already carried the maternal lineages that are dominant in both younger archaeological populations and modern Chinese pigs. Our data set also supports early Neolithic pig utilization and a long-term in situ origin for northeastern Chinese pigs during 8,000–3,500 BP, suggesting a possibly independent domestication in northeast China. Additionally, we observe a genetic replacement in ancient northeast Chinese pigs since 3,500 BP. The results not only provide increasing evidence for pig origin in the middle Yellow River region but also depict an outline for the process of early pig domestication in northeast China.

Romanian wild boars and Mangalitza pigs have a European ancestry and harbour genetic signatures compatible with past population bottlenecks

We aimed to analyse the genetic diversity of Romanian wild boars and to compare it with that from other wild boar and pig populations from Europe and Asia. Partial sequencing of the mitochondrial encoded cytochrome b (MT-CYB) gene from 36 Romanian wild boars and 36 domestic pigs (Mangalitza, Bazna and Vietnamese breeds) showed that the diversity of Romanian wild boars and Mangalitza pigs is fairly reduced, and that most of the members of these two populations share a common MT-CYB haplotype. Besides, in strong contrast with the Bazna animals, Romanian wild boars and Mangalitza swine did not carry Asian variants at the MT-CYB locus. The autosomal genotyping of 18 Romanian wild boars with the Illumina Porcine SNP60 BeadChip revealed that their genetic background is fundamentally European, even though signs of a potential Near Eastern ancestry (~25%) were detectable at K = 4 (the most significant number of clusters), but not at higher K-values. Admixture analysis also showed that two wild boars are of a hybrid origin, which could be explained by the mating of feral animals with domestic pigs. Finally, a number of Romanian wild boars displayed long runs of homozygosity, an observation that is consistent with the occurrence of past population bottlenecks and the raise of inbreeding possibly due to overhunting or to the outbreak of infectious diseases.

Conservation priorities of Iberoamerican pig breeds and their ancestors based on microsatellite information

Criollo pig breeds are descendants from pigs brought to the American continent starting with Columbus second trip in 1493. Pigs currently play a key role in social economy and community cultural identity in Latin America. The aim of this study was to establish conservation priorities among a comprehensive group of Criollo pig breeds based on a set of 24 microsatellite markers and using different criteria. Spain and Portugal pig breeds, wild boar populations of different European geographic origins and commercial pig breeds were included in the analysis as potential genetic influences in the development of Criollo pig breeds. Different methods, differing in the weight given to within- and between-breed genetic variability, were used in order to estimate the contribution of each breed to global genetic diversity. As expected, the partial contribution to total heterozygosity gave high priority to Criollo pig breeds, whereas Weitzman procedures prioritized Iberian Peninsula breeds. With the combined within- and between-breed approaches, different conservation priorities were achieved. The Core Set methodologies highly prioritized Criollo pig breeds (Cr. Boliviano, Cr. Pacifico, Cr. Cubano and Cr. Guadalupe). However, weighing the between- and within-breed components with FST and 1-FST, respectively, resulted in higher contributions of Iberian breeds. In spite of the different conservation priorities according to the methodology used, other factors in addition to genetic information also need to be considered in conservation programmes, such as the economic, cultural or historical value of the breeds involved.

Microsatellite and Mitochondrial DNA Study of Native Eastern European Cattle Populations: The Case of the Romanian Grey

The Eastern European Grey cattle are regarded as the direct descendants of the aurochs (Bos taurus primigenius). Nowadays in Romania, less than 100 Grey animals are being reared and included in the national gene reserve. We examined the genetic diversity among Romanian Grey, Brown, Spotted and Black and White cattle breeds, with a particular focus on Romanian Grey through the use of (i) 11 bovine specific microsatellite markers on 83 animals and (ii) 638 bp length of mitochondrial DNA (mtDNA) D-loop region sequence data from a total of 81 animals. Both microsatellite and mtDNA analysis revealed a high level of genetic variation in the studied breeds. In Romanian Grey a total of 100 alleles were found, the mean number of observed alleles per locus was 9.091; the average observed heterozygosity was 0.940; the Wright's fixation index (FIS) was negative (-0.189) and indicates that there is no inbreeding and no selection pressure. MtDNA analysis revealed 52 haplotypes with 67 variable sites among the Romanian cattle breeds without any insertion or deletion. Haplotype diversity was 0.980 ± 0.007 and ranged from 0.883 ± 0.056 (Brown) to 0.990 ± 0.028 (Spotted and Black and White). The highest genetic variability of the mtDNA was recorded in the Grey breed, where 18 haplotypes were identified. The most frequent mtDNA D-loop region belonged to T3 haplogroup (80.247%), which was found across all studied breeds, while T2 haplotypes (16.049%) was only found in Grey, Spotted and Black and White genotypes. The T1 haplotypes (3.704%) were found in the Grey and Spotted. The current results contribute to the general knowledge on genetic diversity found in Eastern European cattle breeds and could prove a valuable tool for the conservation efforts of animal genetic resources (FAnGR).

Ancient DNA sheds light on the ancestry of pre-hispanic Canarian pigs

Background

Canarian Black (CB) pigs belong to an autochthonous and endangered breed, which is spread throughout the Canarian archipelago. It is commonly accepted that they represent a relic of the pig populations that were bred by the Berbers in North Africa over millennia. It is important to note that the geographic isolation of the Canary Islands has preserved this genetic legacy intact from foreign introgressions until the Spanish conquest of the archipelago in the 15^th century. Ten years ago, it was demonstrated that, in CB pigs, the frequency of the Asian A2 cytochrome-b haplogroup reached 73%. The current work aimed at investigating whether this observation is explained by either a recent or an ancient introgression of CB pigs with Far Eastern pigs.

Results

Genetic analyses of 23 ancient samples from pre-hispanic Canarian pigs (420 to 2500 years before present) showed that Near Eastern and Far Eastern genetic signatures were totally absent in the primitive Canarian pre-hispanic pigs. Indeed, the haplotypes detected in these pigs were closely related to those of North African and European wild boars.

Conclusions

Our results demonstrate that the high frequency of the Far Eastern mitochondrial cytochrome B A2 haplotype in modern Canarian Black pigs probably corresponds to a relatively recent introgression with British breeds.

Electronic supplementary material

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

Genome data from a sixteenth century pig illuminate modern breed relationships

Ancient DNA (aDNA) provides direct evidence of historical events that have modeled the genome of modern individuals. In livestock, resolving the differences between the effects of initial domestication and of subsequent modern breeding is not straight forward without aDNA data. Here, we have obtained shotgun genome sequence data from a sixteenth century pig from Northeastern Spain (Montsoriu castle), the ancient pig was obtained from an extremely well-preserved and diverse assemblage. In addition, we provide the sequence of three new modern genomes from an Iberian pig, Spanish wild boar and a Guatemalan Creole pig. Comparison with both mitochondrial and autosomal genome data shows that the ancient pig is closely related to extant Iberian pigs and to European wild boar. Although the ancient sample was clearly domestic, admixture with wild boar also occurred, according to the D-statistics. The close relationship between Iberian, European wild boar and the ancient pig confirms that Asian introgression in modern Iberian pigs has not existed or has been negligible. In contrast, the Guatemalan Creole pig clusters apart from the Iberian pig genome, likely due to introgression from international breeds.

Mining the pig genome to investigate the domestication process

Pig domestication began around 9000 YBP in the Fertile Crescent and Far East, involving marked morphological and genetic changes that occurred in a relatively short window of time. Identifying the alleles that drove the behavioural and physiological transformation of wild boars into pigs through artificial selection constitutes a formidable challenge that can only be faced from an interdisciplinary perspective. Indeed, although basic facts regarding the demography of pig domestication and dispersal have been uncovered, the biological substrate of these processes remains enigmatic. Considerable hope has been placed on new approaches, based on next-generation sequencing, which allow whole-genome variation to be analyzed at the population level. In this review, we provide an outline of the current knowledge on pig domestication by considering both archaeological and genetic data. Moreover, we discuss several potential scenarios of genome evolution under the complex mixture of demography and selection forces at play during domestication. Finally, we highlight several technical and methodological approaches that may represent significant advances in resolving the conundrum of livestock domestication.

Genetic structure of wild boar (Sus scrofa) populations from East Asia based on microsatellite loci analyses

Background

Wild boar, Sus scrofa, is an extant wild ancestor of the domestic pig as an agro-economically important mammal. Wild boar has a worldwide distribution with its geographic origin in Southeast Asia, but genetic diversity and genetic structure of wild boar in East Asia are poorly understood. To characterize the pattern and amount of genetic variation and population structure of wild boar in East Asia, we genotyped and analyzed microsatellite loci for a total of 238 wild boar specimens from ten locations across six countries in East and Southeast Asia.

Results

Our data indicated that wild boar populations in East Asia are genetically diverse and structured, showing a significant correlation of genetic distance with geographic distance and implying a low level of gene flow at a regional scale. Bayesian-based clustering analysis was indicative of seven inferred genetic clusters in which wild boars in East Asia are geographically structured. The level of genetic diversity was relatively high in wild boars from Southeast Asia, compared with those from Northeast Asia. This gradient pattern of genetic diversity is consistent with an assumed ancestral population of wild boar in Southeast Asia. Genetic evidences from a relationship tree and structure analysis suggest that wild boar in Jeju Island, South Korea have a distinct genetic background from those in mainland Korea.

Conclusions

Our results reveal a diverse pattern of genetic diversity and the existence of genetic differentiation among wild boar populations inhabiting East Asia. This study highlights the potential contribution of genetic variation of wild boar to the high genetic diversity of local domestic pigs during domestication in East Asia.

Whole-genome sequence analysis reveals differences in population management and selection of European low-input pig breeds

BACKGROUND

A major concern in conservation genetics is to maintain the genetic diversity of populations. Genetic variation in livestock species is threatened by the progressive marginalisation of local breeds in benefit of high-output pigs worldwide. We used high-density SNP and re-sequencing data to assess genetic diversity of local pig breeds from Europe. In addition, we re-sequenced pigs from commercial breeds to identify potential candidate mutations responsible for phenotypic divergence among these groups of breeds.

RESULTS

Our results point out some local breeds with low genetic diversity, whose genome shows a high proportion of regions of homozygosis (>50%) and that harbour a large number of potentially damaging mutations. We also observed a high correlation between genetic diversity estimates using high-density SNP data and Next Generation Sequencing data (r = 0.96 at individual level). The study of non-synonymous SNPs that were fixed in commercial breeds and also in any local breed, but with different allele, revealed 99 non-synonymous SNPs affecting 65 genes. Candidate mutations that may underlie differences in the adaptation to the environment were exemplified by the genes AZGP1 and TAS2R40. We also observed that highly productive breeds may have lost advantageous genotypes within genes involve in immune response--e.g. IL12RB2 and STAB1-, probably as a result of strong artificial in the intensive production systems in pig.

CONCLUSIONS

The high correlation between genetic diversity computed with the 60K SNP and whole genome re-sequence data indicates that the Porcine 60K SNP Beadchip provides reliable estimates of genomic diversity in European pig populations despite the expected bias. Moreover, this analysis gave insights for strategies to the genetic characterization of local breeds. The comparison between re-sequenced local pigs and re-sequenced commercial pigs made it possible to report candidate mutations to be responsible for phenotypic divergence among those groups of breeds. This study highlights the importance of low input breeds as a valuable genetic reservoir for the pig production industry. However, the high levels of ROHs, inbreeding and potentially damaging mutations emphasize the importance of the genetic characterization of local breeds to preserve their genomic variability.

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.