Skip to main content
Download PDF
  • Research Article
  • Open access
  • Published:

Genome-wide SNP data unveils the globalization of domesticated pigs

Genetics Selection Evolution volume 49, Article number: 71 (2017) Cite this article

A Correction to this article was published on 04 June 2020

This article has been updated

Abstract

Background

Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection.

Results

A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection.

Conclusions

Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

Background

Domestication of pigs from wild boars occurred independently in Asia and Europe about 10,000 years ago [1]. Due to the biogeographic difference between wild ancestral populations, which results from 1.2 million years of separation, Asian and European pigs are genetically highly divergent [2,3,4,5]. Sus scrofa is native to Eurasia and North Africa, but was introduced into other parts of the world, i.e. into the Americas, primarily in its domesticated form, during the time of the European colonization in the sixteenth century, and later in Australia and New Zealand [6]. Both demographic processes, and natural as well as artificial selection, have led to the formation of a multitude of pig breeds around the world that vary in coat color, ear shape, body size, snout bluntness, behavior, growth rate, fatness, and prolificacy and other economically important traits.

In addition to domestication, crossbreeding between Asian and European indigenous pigs mediated by humans are significant landmarks in pig breeding history. Although anecdotal evidence exists even from the classical era, admixture between Western and Eastern pigs only started to become common in the mid- to late eighteenth century [7]. Introduction of Chinese pigs into Britain is documented from then and its aim was to improve the production characteristics of local pigs, which led to the creation of modern breeds such as Yorkshire (i.e., Large White), Berkshire and Hampshire [8]. In the late eighteenth century, Chinese pigs may also have been imported to America, and crossed with local pigs of European ancestry there [8], although most likely the Asian influence in American village pigs was through crosses with international breeds [9]. Reciprocally, at least since the 1840s, modern breeds such as Berkshire, Hampshire, Russian local pigs, Duroc, Large White and Landrace were introduced into China [10]. Such domestic animals were traded, loaded onto ships and released elsewhere. This is well documented, e.g., during the exploration of the Pacific by Captain Cook, who is credited for having released the first pigs on the New Zealand islands [11]. As in Europe, these imported pigs were used for crossbreeding with local breeds. However, in China, the introduction of pigs from outside and trading of pigs within China, appear to have been less widespread until recently, as is apparent from the high degree of geographic structure that remains in the Chinese traditional pig breeds [3]. Nevertheless, historical records and genetic evidence point to the contribution of European pigs to some East Asian breeds. For instance, the modern Korean Native pig is a cross between a local, traditional Korean pig and Berkshire. More recently, since the 1980s, Chinese pig breeders began programs to improve local breeds using Western stock [10] by creating synthetic breeds. In Africa, although advocated as an additional center of domestication, most of the evidence points to introgression from foreign breeds. Interestingly, Asian haplotypes predominate in East Africa, whereas European haplotypes predominate in West Africa [12].

Today, pig is a major livestock species, which in 2012 represented about 36.3% of the total meat production for human societies (www.fao.org), with major contributions from only a few international commercial breeds (i.e. Duroc, Large White, and Landrace). Nevertheless, hundreds of domesticated pig breeds worldwide [13] are still important for local meat production by small farmers. Many of these pig breeds have unique characteristics that differ from those of the international commercial breeds. Conservation of agrodiversity is one of the pillars to maintain food security, particularly in a rapidly changing world where consolidation of international plant and animal breeds is resulting in an increasingly narrow genetic basis for food production [14]. Thus, indigenous pig breeds, together with their wild relatives, are valued resources for the human society, not only for food, but also as genetic reservoirs. In addition, they constitute cultural and historical value since certain breeds are highly connected to local identity and specific agricultural practices. Finally, breed diversity can be leveraged for understanding the genetic basis of complex traits and adaptive evolution [15]. Large-scale genotyping technologies have enabled the analysis of the genetic ancestry and admixture of many domestic animals, including dogs [16], cattle [17], and pigs [18] and have also enabled the characterization of the genetic basis of phenotypic changes during domestication in chicken [19], dogs [20], rabbits [21] and pigs [22].

To date, population genetic studies using genomic data in pigs had a limited, usually regional, scope [9, 23,24,25]. Compared to previous generations of molecular markers, particularly microsatellites, single nucleotide polymorphism (SNP) markers allow for relatively straightforward data integration across studies since SNP genotypes can be compared unambiguously across studies. The aim of the current study was to perform a truly global integration of pig genotype data through the analysis of 1839 domestic pigs from 122 indigenous pig breeds that were collected in 29 different countries, together with 215 wild boars and 39 out-group individuals. As a result, our findings constitute a big leap in understanding the population structure, admixture, demographic history, and characterization of genetic loci involved in the domestication of pigs globally.

Methods

Samples and data

The raw Illumina 60K SNP data [26] of 3482 pigs, which include 3443 Sus scrofa and 39 non Sus scrofa suids (outgroups) (Table 1), were mainly obtained from three sources (see Additional file 1: Table S1): Wageningen University in The Netherlands (2464 individuals that encompass pig populations from Europe, Asia, Africa, Oceania, North America, international commercial pig populations, as well as outgroup suids), Jiangxi Agricultural University in China (821 individuals, which mainly consisted of pig populations from China, Russia and Ukraine), and the Autonomous University of Barcelona in Spain (197 individuals, which mainly represent pig populations from South America and Iberian pigs). Genomic SNP positions are based on the genome assembly Sscrofa 10.2 (EnsEMBL db version 83) [18].

Table 1 Number of populations and samples by continent
Full size table

We conducted a series of quality control procedures on the raw data using PLINK v1.9 [27]. First, we excluded the breeds with less than five individuals. For the breeds or populations with more than 20 individuals, we randomly removed one individual from a pair of highly related animals (identity by state score > 0.95), and then kept the top 20 samples ranked by the SNP call rate. Next, we removed SNPs with a minor allele frequency (MAF) lower than 0.01, a call rate lower than 90% and individuals with a call rate lower than 90%, which resulted in a dataset of 55,072 SNPs for 2093 individuals that was used to estimate ROH, haplotype diversity and effective population size. We further removed SNPs with a MAF lower than 0.05, and in high linkage disequilibrium (r2 > 0.2) by using—maf 0.05 and—indep-pairwise 50 10 0.2 in PLINK v1.9 [27], respectively, and generated a dataset of 15,427 SNPs for 2093 individuals for subsequent multi-dimensional scaling (MDS), neighbor joining (NJ) tree and admixture analyses. See Additional file 2: Table S2 for further details.

Statistical analysis

Population structure

MDS was carried out using—mds-plot and—cluster options in PLINK v1.9 [27] and visualized by R programming language [28]. The NJ tree was constructed using PHYLIP v3.69 [29] based on the identical by states matrix obtained by PLINK v1.9 [27], and visualized using FigTree v1.4 (http://beast.bio.ed.ac.uk/figtree). To facilitate visualization, we randomly selected six individuals from each population to build up the NJ tree. The geographical maps were plotted using R package MAPS [30] and MAPPLOTS (https://cran.r-project.org/web/packages/mapplots/). The coordinates of longitude and latitude of each population were set according to where the pigs were sampled (see Additional file 1: Table S1). The geographical distances between each pair of breeds were computed using distm function in R package GEOSPHERE (https://cran.r-project.org/web/packages/geosphere/). The proportion of mixed ancestry in the populations analyzed was evaluated by the ADMIXTURE 1.22 program [31]. We evaluated different K values with the mixed ancestry model (K = 2 to 17).

Runs of homozygosity, haplotype diversity and effective population size

Runs of homozygosity (ROH) of each breed were identified using PLINK v.1.07 by a 5-Mb sliding window process across the genome with at least 50 SNPs, allowing five missing calls and one heterozygous SNP. The minimum length for ROH was set to 500 kb. ROH statistics were then transformed to FROH. We inferred haplotypes of autosomes for all individuals using SHAPEIT v2 [32]. Haplotype diversity was calculated for populations with a minimum of 10 individuals. For each population, we randomly selected 10 individuals for the analysis. The haplotype diversity of a population was measured as the average number of haplotypes in windows of 5, 10, and 15 SNPs, respectively, which is similar to the method described in [16]. LD between adjacent SNPs was measured by the genotype correlation coefficient (r2) calculated by the—r2—ld-window 99999—ld-window-r2 0 command in PLINK v1.9. We used the same equation to fit the relationship between LD and genetic distance as previously described [33]. The SNPs used to calculate the LD within a population were filtered by applying the following criteria: a MAF higher than 0.05 and a P value for Hardy–Weinberg equilibrium higher than 1 × 10−6. The effective population size was estimated according to Sved [34], based on the equation r2 = 1/(4Ne c + 1), where r2 is the linkage disequilibrium between a pair of SNPs, Ne is the effective population size, c is the genetic distance in Morgan between a pair of SNPs, which was obtained by multiplying their physical distance and recombination rate [35]. The Ne at generation T, were obtained by the equation T = 1/2c, the same as described in [36].

Domestication loci

To detect loci that may have been selected during domestication, we calculated the fixation index (Fst) [37] between domestic and wild boars in Asia and Europe, separately. To avoid the influence of introgression, we used 42,808 SNPs (MAF > 0.05) in 782 individuals that have more than 90% of Asian or European ancestry in the analysis. We ranked the Fst values of genome-wide SNPs, and genes within a 100-kb region of high-Fst SNPs were identified as candidate genes that may have been involved in past selection. We selected candidate genes according to their functional relevance to phenotypes, such as e.g. behavior, development, energy metabolism, which may confer differences between domestic pigs and wild boars.

Results

Samples

After quality control, a dataset of 2093 samples representing 122 domestic pig breeds (1839 individuals), 19 wild boar populations (215 individuals) and five out-group populations (39 individuals) was available. The 122 domestic breeds included 104 local breeds or synthetic populations and 18 international commercial populations (Duroc, Large White, Landrace and Pietrain from different countries were considered as different breeds) (see Additional file 1: Table S1). Among the 104 local breeds, 39 originated from Europe, 40 from Asia, 22 from the Americas and three from other parts of the world. The wild boars were from widespread regions around the world and the five out-group populations are Babyrousa babyrussa, Sus barbatus, Sus celebensis and Sus verrucosus from the islands of Southeast Asia and Phacochoerus africanus from Africa. A more detailed description of many of these samples was reported in previous studies [9, 24, 25].

Global population structure

We performed multi-dimensional scaling (MDS) analysis on 15,427 SNPs for the 2093 individuals to investigate the genetic relationships between populations of domesticated pigs, wild boars and out-group populations. The first principal component separates Asian and European breeds (Fig. 1) and (see Additional file 3: Figure S1), in agreement with independent domestication and evolution of pigs in Asia and Europe. The American, African, and international commercial populations are much closer to European than to Asian pigs, which indicates a predominant contribution of European ancestry in the formation of these derived populations. Nevertheless, several populations are positioned in between the Asian and European main clusters. These include Oceania populations, populations from China (Lichahei and Sutai), the Korean native pig breed, Minisibs pigs from Russia and a Large White × Meishan F1 cross (Fig. 1) and (see Additional file 3: Figure S1). The intermediate position in the MDS reflects the fact that these populations are derived from both Asian and European pigs, and we validated the hybrid nature of these populations in subsequent admixture analysis. The Lichahei and Sutai pigs from China, the LW × Meishan F1 cross, and Minisib pigs from Russia consist of approximately 50% Asian and 50% European ancestries (see Additional file 4: Figure S2).Thus, the first MDS axis represents the Asian and European ancestries.

Fig. 1
figure 1

Global genetic structure of pig populations in this study. Multi-dimensional scaling analysis of pig populations from the five continents. Each point represents breed-average coordinates of eigenvalues of principle components 1 and 2. Points are mainly colored based on geographic origin of breeds; blue Asia, red Europe, green America, yellow Africa, purple Oceania, brown outgroup of Suids, black commercial Breeds; the pig populations in the middle of the graph representing admixed Asian and European ancestries are annotated

Full size image

The second axis separates domestic pigs from wild boars in both Asia and Europe (Fig. 1) and (see Additional file 3: Figure S1), which indicates that domestication and subsequent artificial selection also resulted in the differentiation between wild and domesticated populations. Regarding the international commercial pig breeds, the Duroc breed tends to cluster with American domestic breeds and to separate from other commercial breeds (Landrace, Large White and Pietrain) that tend to cluster together. These results agree with the fact that Duroc pigs were originally developed in the USA, whereas, the other international commercial breeds (e.g. Large White-England, Landrace-Denmark, and Pietrain-Belgium) originated in Europe.

The results from hierarchical clustering (neighbor joining based on identity by state distance metric) generally agree with those from the MDS analysis. An important observation is that, even on a global scale, all populations of the major commercial breeds, still cluster together, i.e. breed identity has been maintained for both commercial and non-commercial populations (see Additional file 5: Figure S3).

Regional population structures in Asia and Europe

To infer geographic region-specific details, we performed separate MDS analyses on pig populations within Asia and Europe, separately (Fig. 2). In Asia, most of the pig populations originate from China. In previous studies [3, 38], Chinese pigs were grouped into six categories that included pig types from Central China, the Yangtze River basin (East China), South China, Southwest China, North China, and Plateau (West and Northwest China), according to their external traits and geographical distributions [38]. The genetic clusters revealed in the MDS analyses are broadly concordant with this assignment into six categories (Fig. 2a, b). Pig breeds from South China and those from East and North China are located at each end of the first axis, while pig breeds from Central and West China are located in between (Fig. 2a, b). As expected under a model of gene flow between populations that is inversely proportional to physical distance, genetic distances were significantly correlated with geographical distances (Pearson correlation, P value = 1.5 × 10−62) (Fig. 2c). This is in sharp contrast to the results for pig breeds from the Americas where no concordance between genetic similarity and geographical distances was observed due to the complex colonization and breeding history of American pigs [9].

Fig. 2
figure 2

Regional genetic structure of indigenous pig populations in Asia and Europe. a, d show the results of the MDS analysis of pig populations from Asia (40 breeds) and Europe (26 breeds). Each point represents a breed, colors are assigned to each breed according to their geographical distributions, which are visualized in (b) and (e) for Asian and European pigs, respectively. c, f show the correlation between genetic and geographic distances among pig breeds in Asia and Europe, respectively. The legends of pig breeds are shown on the right. The upper legends in the blue box are for Asian breeds and the legends below in red box are for European breeds

Full size image

In Europe, the pig breeds from Southern Europe (Italy, Spain and Hungary) and those from North and middle Europe (Netherlands, Sweden, Poland, Germany and the Czech Republic) were genetically distinct. This distinction is represented by the first axis (Fig. 2d, e). Within the UK, British Lop, Leicoma, Middle White and Welsh breeds differ from Large Black, Gloucester Old Spot, Saddleback and Tamworth breeds. However, we observed no correlation with geographical distances among pig breeds in Europe (Fig. 2f), which is consistent with previous results based on microsatellites [3]. The absence of population structure in European pigs is explained, at least in part, by the Asian introgression and subsequent influence of highly productive “international” breeds on local pig diversity. In Europe, many ‘local’ or ‘traditional’ breeds have effectively become (partially) extinct due to such extensive crossbreeding [36, 39, 40].

Global genetic ancestries

We further examined the genetic ancestry of pig populations worldwide by varying the number of ancestries (K) in ADMIXTURE v 1.2 [31] (Fig. 3) and (see Additional file 4: Figure S2). The population structure generally agreed with the MDS results. At K = 2, the two ancestries clearly reflect Asian and European origins. The pig populations from America, Africa, Russia and neighboring countries (including Ukraine, Belorussia and Kazakhstan) are mainly of European ancestry and (see Additional file 4: Figure S2). Even the international commercial pig breeds are mostly of European origin although they have a large Asian ancestry component [23]. At K = 8, we found two distinct Asian ancestries that are represented by pig breeds from East (Meishan, CNMS) and South (Luchuan, CNLU) China. The other six ancestries are represented by European wild boars, Hampshire (UKHS) and Berkshire (UKBK), and four international commercial breeds including Duroc, Large White, Landrace, and Pietrain and (see Additional file 4: Figure S2), (K = 8). At a higher K value (K = 17), we found various breed-specific ancestries that reflect a recent isolated breeding history (Gansu Zang (CNGS), Meishan (CNMS), Luchuan (CNLU), Jinhua (CNJH) and Congjiangxiang (CNCJ) pigs from China, Tamworth (UKTA), Hampshire (UKHS) and Berkshire (UKBK) from Europe, and Mulefoot pig (USMU) from USA).

Fig. 3
figure 3

Landscape of worldwide admixture of pig populations. a Bar plot of admixture analysis (K = 17). Each vertical bar stands for an individual, the colors represent different ancestries. b Worldwide map of admixture for pig populations, the pie plots represent different breeds, the compositions of ancestries for a breed were calculated from the averages of ancestry composition of individuals within that breed. c, d Regional plot for Asian and European pigs, respectively

Full size image

Admixture between Asian and European ancestries

An additional interesting pattern is the widespread admixed ancestries that we observed in pig populations across the different continents. Note that many breeds showed evidence of both Asian and European ancestries. Since these two lineages evolved independently, this admixture is a result of human-mediated activities. In Europe, importation of Chinese pigs is well accredited since the eighteenth century.

In Asia, 23 (57.5%) of the 40 breeds analyzed have 5% or more inferred European ancestry (see Additional file 4: Figure S2 and Additional file 6: Figure S4), (K = 2). Most of the European introgression was inferred to be from Duroc, Landrace, Large White, Berkshire and Hampshire breeds (see Additional file 4: Figure S2 and Additional file 6: Figure S4), (K = 13). Eight Asian breeds have more than 20% of genomic introgression from European pigs. These breeds include a Korean local breed (KPKO), a Thailand local breed (THCD), and six breeds from China (Lichahei (CNLC) from Shandong Province, Sutai (CNST) from Jiangsu Province, Kele (CNKL) and Guanling (CNGU) pigs from Guizhou Province, Leanhua (CNLAH) from Jiangxi Province, Minzhu (CNMZ) from Northeast China. Among these breeds with a large degree of European introgression (>20%), the Korean pig, a known East–West synthetic breed, formed the largest European ancestry group (see Additional file 4: Figure S2 and Additional file 6: Figure S4), (K = 2 and K = 13), including ancestry from Berkshire, Hampshire, Landrace and Duroc, which reflects the complex breeding history of this breed (Fig. 3a) and (see Additional file 4: Figure S2) (K = 13 and K = 17). This is largely in line with the known origin of the Korean local pigs. Sutai and Lichahei have been mainly admixed with Duroc, while Min pigs have a considerable contribution from Berkshire (Fig. 3a) and (see Additional file 4: Figure S2 and Additional file 6: Figure S4) (K = 8 and K = 17). It is interesting that the admixture with European pigs occurred mainly in Western and Northern Chinese pig breeds (Gongbujiangda (CNXZ) and Milin Tibetan (CNML) pigs, Kele and Guanling pigs from Guizhou Province, Mingguangxiaoer (CNMG) pigs from Yunnan province, Bamei (CNBM) pigs from Gansu province, Laiwuhei (CNLH) pigs from Shandong Province, Hetaodaer (CNHT) from Inner Mongolia and Min (CNMZ) pigs from Heilongjiang Province); the European ancestries that are involved encompass Large White, Landrace, Berkshire or other European breeds (Fig. 3a–c) and (see Additional file 4: Figure S2 and Additional file 6: Figure S4) (K = 13 and K = 17). In comparison, the pig breeds from South and Central China, including Erhualian, Xiang, Dongshan, Shaziling, Congjiangxiang, Lantang, Jinhua, Litang Tibetan and Luchuan, show no or negligible introgression from European pigs (Fig. 3a, c).

Iberian pigs from Spain, Cinta Senese and Nera Siciliana pigs from Italy, and Mangalica pigs from Hungary showed little evidence of influence from Asian pigs (see Additional file 4: Figure S2 and Additional file 7: Figure S5). By contrast, there is evidence of introgression from Asian pigs for all other pig breeds from Europe, including those from Ukraine and Russia (see Additional file 4: Figure S2 and Additional file 7: Figure S5). These results confirm the widespread Asian influences in European breeds.

The North and South American samples consisted mainly of village and feral pigs from eight countries [9]. Consistent with previous studies on pigs from the Americas using 60K SNP [9], and mitochondrial DNA data [41], pig populations from rural areas have mosaic genetic compositions that consist of multiple ancestries from both Europe and Asia. The largest ancestry components were similar to Iberian pigs (ESIB), in agreement with a primigenious origin from the Iberian Peninsula. Other European components are related to Duroc, Landrace, Berkshire, Hampshire, and European Wild boars (Fig. 3a) and (see Additional file 4: Figure S2 and Additional file 8: Figure S6). Intriguingly, a considerable contribution from pigs from both east and south China was observed in most of the American Village pigs (Additional file 8: Figure S6). In general, the village pigs from Brazil, Mexico and Cuba have larger Asian components than the pigs from other American countries (Fig. 3a).

In Africa, the Tunisian wild boar shows a high degree of similarity to the European wild boar, and specifically to the wild boar from the Iberian Peninsula. A local breed from Kenya was inferred to contain both Asian and European ancestries (Fig. 3a) and (see Additional file 4: Figure S2 and Additional file 8: Figure S6). This is in agreement with mtDNA studies, which showed that Asian haplotypes were abundant in East Africa but completely absent in the Northern African pigs (i.e. Tunisian wild boars) [41].

In Oceania, the Australian feral pigs also show admixed ancestry from Asian and European pigs (Fig. 3a) and (see Additional file 4: Figure S2 and Additional file 8: Figure S6).

The four major international commercial pig breeds, i.e. Duroc, Landrace, Large White, and Pietrain, have a considerable percentage of Asian ancestry (see Additional file 4: Figure S2 and Additional file 9: Figure S7) (K = 2). At K = 8, these four breeds form four distinct ancestries. The Landrace and Large White breeds also showed a diversity of genetic ancestries related to various pig populations including Berkshire, Hampshire, South European local pigs or Asian pigs.

Genetic diversity

We analyzed runs of homozygosity (FROH), haplotype diversity, and effective population size for each pig population to assess their inbreeding history and effective population size. Previous studies showed that 60K SNP data provide reasonably accurate estimates of long FROH [15]. We calculated the total length of ROH with a minimum length of 500 kbp for each individual. Considerable variation in FROH occurs within and across populations, which reflects the complex breeding history of pigs (Fig. 4). The cumulative length of ROH ranged from 4.98 Mb for the Dutch LW × Meishan F1 population from the Netherlands to 591.57 Mb for the Mora Romagnola pigs from Italy, and represented between 0.2 and 20.8% of the genome. Since the Dutch LW × Meishan is an F1 cross, it was expected to have few if any ROH. The 10 populations with the highest FROH included the Mora Romagnola (ITMR) and Cinta Senese pigs (ITCS) from Italy, the Mangalica breed from Hungary (HUMA), the Korea local breed (KPKO), the Mulefoot (USMU) and Yucatan Mini pigs (USYU) from USA, the Creole pigs (CRCR) from Costa Rica, the Gloucester old spot pigs (UKGO) and Tamworth pigs (UKTA) from UK, and the Leanhua pigs (CNLA) from China, which indicates that these populations have recently experienced considerable inbreeding (Fig. 4) and (see Additional file 1: Table S1). In addition, the total length of ROH was negatively correlated with haplotype diversity (Pearson correlation coefficient = -0.71, P value = 1.6 × 10−20) (see Additional file 10: Figure S8). Therefore, populations with a high FROH normally have a low haplotype diversity (Fig. 4) and (see Additional file 11: Figure S9). The effective population size (Ne) for each population was estimated using linkage disequilibrium following the method described in [36] (see Additional file 1: Table S1). Considering only the populations with a minimum of 10 individuals, the estimated Ne of the past five generations is between 26 and 67. Even after accounting for a small systematic bias towards lower estimated Ne in populations that a have smaller sample size, it is clear that the Ne of indigenous pig populations is generally smaller than that of the international commercial breeds (see Additional file 12: Figure S10 and Additional file 1: Table S1).

Fig. 4
figure 4

Distribution of FROH for pig populations across the world. a Box plot showing the global distribution of total length of FROH for pig populations worldwide, each point represents the total length of FROH for one individual, each box represents one breed, the colors and the order of breeds are the same as those described in Fig. 3. Regional view of FROH distribution for pig populations in Asia (b), Europe (c), America, Oceania and Africa (d), and international commercial breeds (e)

Full size image

Loci involved in domestication

Domestication and artificial selection have resulted in a wide range of phenotypes across domestic pig breeds that differ from their wild relatives. These are related to behavior, body size, fertility, locomotion ability and adaptation to feed provided by humans. To detect genetic loci that could be involved in the transition from wild to domestic, we calculated the genome-wide fixation index (Fst) between domestic pigs and wild boars in Asia and Europe, separately (see “Methods” section) (Fig. 5). Empirically, we considered the 428 (1%) SNPs with the highest Fst values as potential loci under recent (domestication) selection. Only six outlier SNPs were shared between Asia and Europe, which only slightly exceeds the number expected based on re-sampling of SNPs (see Additional file 13: Figure S11). Thus, we found no evidence for specific loci being under selection during the independent domestication processes in Asia and Europe. We examined the genes that are located within 100 kb to the top outlier SNPs with extreme Fst values, and made the assumption that many of the genes around the top 30 outlier SNPs are involved in functions that are associated with phenotypic changes from wild to domestic (Fig. 5; Table 2). For Asian pigs, we identified genes related to muscle development (MSTN [42]), energy balance (NMU [43], LEP [44] and GSK3A [45]), social behavior (TBX19 [46] and PAFAH1B3 [47]), puberty and reproduction (GNRHR [48], ESR1 [49] and PATZ1 [50]) and perception of smell (Olfr466 [51]) (Table 2). For European pigs, we identified genes related to growth and body development (SOX2-OT [52]), cardiac system development (TBX20 [53]), metabolism of protein, glucose or fatty acid (TMEM67 [54], FOXA1 [55] and INSIG2 [56]), central nervous system (LRRC4 [57], VEPH1 [58] and CDH9 [59]), immune system (LAIR1 [60]), and reproduction (PLSCR4 [61]) (Table 2).

Fig. 5
figure 5

Genome-wide analysis of global Fst between domestic pigs and wild boars. Manhattan plot of genome-wide Fst values between domestic pigs and wild boars in Asia (a) and Europe (b). Fst values are shown on the y axis, and genomic positions on the x axis. The different chromosomes are represented by different colors

Full size image
Table 2 Candidate genes for domestication loci in Asia and Europe
Full size table

Discussion

Pig is one of the most important livestock species for humans as a valued, global, resource for meat production and as an excellent animal model to understand the genetic mechanisms that underlie complex traits [6, 18]. Its long domestication history, originating from a large diversity of wild ancestors throughout Eurasia, and selection for economic and cultural purposes have resulted in a large number of breeds globally, which show a wide phenotypic diversity. Our worldwide survey on SNP data from 122 breeds/populations and 215 wild boars worldwide, reveals genetic ancestries, introgression and inbreeding histories of pigs at a global scale and at an unprecedented detail. Although there are potential issues regarding ascertainment bias [40] associated with the SNP assay used in this study, admixture analyses using 60K SNP and 30 million SNPs called from whole-genome sequence data provided very similar results (see Additional file 14: Figure S12), which indicates that robust conclusions can be drawn from the 60K SNP assay data for a wide population study as presented here.

Population structure

The high degree of geographic structure observed here in the Asian domesticated pigs agrees with a previous report based on microsatellite markers [3], and differs substantially from that observed in pig populations from Europe and the Americas, in which almost no correlation between genetic and geographical distances exist [9]. The strong concordances between genetic and geographical distances for the pig populations in Asia may be attributed to the fact that pig populations within certain eco-geographical regions are more likely to have common ancestries, and that most of the breeds in Asia did not migrate over large distances. Furthermore, introgressions from European populations did not mask the identity of most Asian breeds, at least not to a large extent. Removing the breeds with more than 20% European introgression resulted in an increased correlation between genetic and geographical distances (see Additional file: 15 Figure S13). This indicates that admixture with geographically distant populations could be a major force in breaking regional genetic-geography concordance, as has been the case in Europe. Recent breed interchanges have largely masked an underlying geographic signal. However, it is interesting to note that some breeds have remained relatively unchanged for centuries. For instance, Ramirez et al. [62] showed that the modern Iberian breed is genetically very similar to a sixteenth century Spanish pig.

Contribution of Chinese populations to worldwide pigs

The Chinese ancestries in European pigs observed in this study confirmed, on a broader population scale, the findings of previous genetic studies [2, 23, 41, 63]. These results are consistent with the historical record that South Chinese pigs were brought to England from Guangzhou in South China, the only treaty port city in China at that time, and have contributed to local British breeds, such as Berkshire and Yorkshire around 200 years ago [7, 8, 63]. Interestingly, our analyses revealed that ancestry represented by Lantang pigs from Guangdong Province is likely the major source of introgression in American pigs (Fig. 3a). The Meishan pigs of Eastern China, a breed famous for its high prolificacy, were imported to Western countries including France, England and USA in the 1980s [64]. Meishan pigs were used in experimental crosses to study the genetic basis of complex traits [65]. Recent studies showed that many Asian alleles with favorable phenotypic effects reached a high frequency in European pigs. These included MC1R alleles that are associated with black coat color [66], an IGF2 allele for muscle growth [67], and AHR alleles for sow reproduction traits [23]. These studies underscored the importance of Asian pigs as vital genetic resources for international pig breeding and pork production.

Contribution of European ancestry to worldwide pig populations

Both MDS and admixture analysis showed that European pigs were the major contributors to pig populations in those regions of the world where Sus scrofa does not occur natively (America, Africa, and Oceania). These results are consistent with mitochondrial and Y chromosome polymorphisms [41]. This contribution is due, in part, to the waves of colonization by Europeans since the sixteenth century. In addition, recent increase in worldwide trading of commercial, improved, pigs throughout the globe, and the desire of local farmers to improve their pigs using these western breeds, have likely contributed to this process as well. Populations in the Americas, Africa, and Oceania tend to harbor multiple ancestries of Mediterranean countries and/or international commercial breeds such as Berkshire, Hampshire, and Duroc (Fig. 3), which indicates a very dynamic process of global mixing of populations during several centuries.

More recently, the global process of mixing has become ‘full circle’ by the introduction of European pigs, themselves heavily influenced by Asian pigs, in Asia. In fact, one of the main original findings of our study is the widespread European influence in many Asian populations, the extent of which was mostly unknown until now. In Asia, we observed widespread and complex gene flow from European pigs, which indicates that many Asian indigenous pig breeds are no longer strictly Asian, but also contain a genetic component of European origin. Occurrence of European introgression in Japan [68], Korea [69] and Vietnam [70] was reported before. In China, there are over 80 pig breeds and a high diversity in phenotypes [64]. Historical documents indicate at least three waves of introgression from European pigs since the 1840s. The first wave of introgression may have occurred around the 1840s, when European pig breeds including Berkshire, Large White, Duroc, pigs from Russia, and Tamworth were brought to China by Germans and Japanese [10]. Subsequently, starting from the early twentieth century, probably since the 1910s, large-scale importation of Western European pigs, such as Berkshire in the Hebei, Sichuan, and Jiangsu provinces, and of Russian pigs in Northwest China, took place to improve local breeds. This study demonstrates that many pig breeds from West and North China contain ancestries from European pigs, notably Berkshire, Hampshire, Large White and Russian pigs, which is in agreement with historical records. Since 1937, war and civic and economic upheavals hampered systematic breeding, which may explain why most of the pig populations in China maintained their geographic identity in spite of admixture with European pigs. Since the 1980s, due to changes in Chinese policies regarding the introduction of foreign agricultural germplasms, many international commercial breeds were introduced into China, which gave rise to several synthetic breeds, such as the Lichahei from Shandong Province, and Sutai pigs from Jiangsu Province. Both breeds currently display considerable ancestry from Duroc.

Indications for conservation of indigenous pigs

Inbreeding and decrease in effective population size may reduce the fitness of a population in response to challenges from changing environments or infectious diseases. This study provides an overview on the inbreeding, demography and admixture history of pig populations worldwide. First, our analysis revealed that 40 breeds or populations have substantial cumulative ROH (>200 Mb), and also exhibit low haplotype diversity, which indicate that these populations underwent recent inbreeding. This may reflect the fact that all domesticated and many wild populations are de facto under population management, deliberate or not. Second, we show that many of the indigenous pig breeds have smaller Ne than those of international commercial breeds. Since the commercial breeds are also the breeds that are the most admixed, this is not surprising. Finally, we found prevalent admixture of Asian and European ancestries in the indigenous pig populations, which suggests that many breeds have become less representative of the original local ancestries. The admixture of populations, particularly between East and West, has resulted in a re-shaping of the nucleotide diversity in the genomes of modern pig breeds. Because of that, only some of the current least admixed breeds may represent the original nucleotide and haplotype diversity in Europe. These results could help to make decisions on the conservation and management of pig populations. For example, Mangalica pigs from Hungary (HUMA) and Mora Romagnola pigs from Italy (ITMR) present the most extensive ROH in their genomes and the highest European ancestry, which indicate that these two breeds have undergone intensive inbreeding, and require special attention regarding conservation measures.

Genetic basis of domestication

Domestication of plants and animals has been one of the major transitions in human history. Farming practices have not only altered the human societies but the interactions with nature, especially for domesticated plants and animals. Domestication of pigs has led to dramatic phenotypic changes transforming the wild boar into pigs by altering their behavior, morphology, coat color, reproduction and physiology. The admixture and MDS analyses presented in this study confirm the close relationship between wild and domesticated Sus scrofa in the geographic areas where domestication took place. Therefore, the genomic regions that show a much higher than average differentiation between wild and domesticated pigs should be enriched for loci under selection during domestication. We identified a number of genes that are located near loci with extreme Fst values and that have functions that match the phenotypic changes from wild boars to domestic pigs. For instance, domestic pigs receive a stable feed supply from humans, while wild boars need to endure starvation if they cannot find food in the wild. We identified a number of genes with functions related to energy balance and metabolism (NUM, LEP FOXA1 and INSIG2), which could have contributed to the adaptation of pigs to food scarcity or abundance. Genes involved in growth (MSTN and SOX2-OT) and reproduction (GNRHR, PATZ1, ESR1, and PLSCR4) could be associated with improved meat production and reproduction traits in domestic pigs that have undergone strong artificial selection. Lastly, genes related to nervous system and behavior (TBX19, LRRC4, VEPH1 and CDH9) could be associated with changes in the behavior of domestic pigs compared to wild boars. The absence of signatures of selection found in previous studies [22,23,24] can be attributed to the higher density of SNPs, the specific selection-detection method, or the application of specific population contrasts in those studies. While further studies are needed to validate the role of the genes that we identified here in the domestication process, our findings confirm that this long-standing genetic experiment—i.e. domestication—is continuing to yield insights into biology and evolution.

Conclusions

We present the largest population study on pigs and their wild ancestors to date, which investigates the population structure and introgression of worldwide pig populations globally. We demonstrate regional and global mixing of pig diversity, which reflect that this species has essentially followed many of the globalization events over the past centuries. Population diversity statistics such as ROH provided insight on inbreeding history and effective population sizes that allow us to recommend guidelines for breeding and conservation programs. Similar to other domesticated species, pigs represent an excellent model to study adaptation. We have identified a number of candidate genes that could have been under positive selection during domestication.

Change history

  • 04 June 2020

    An amendment to this paper has been published and can be accessed via the original article.

References

  1. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, Robins J, et al. Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science. 2005;307:1618–21.

    Article  CAS  PubMed  Google Scholar 

  2. Giuffra E, Kijas JM, Amarger V, Carlborg O, Jeon JT, Andersson L. The origin of the domestic pig: independent domestication and subsequent introgression. Genetics. 2000;154:1785–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Megens HJ, Crooijmans RP, San Cristobal M, Hui X, Li N, Groenen MA. Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas of domestication. Genet Sel Evol. 2008;40:103–28.

    PubMed  PubMed Central  Google Scholar 

  4. Ai H, Fang X, Yang B, Huang Z, Chen H, Mao L, et al. Adaptation and possible ancient interspecies introgression in pigs identified by whole-genome sequencing. Nat Genet. 2015;47:217–25.

    Article  CAS  PubMed  Google Scholar 

  5. Frantz LA, Schraiber JG, Madsen O, Megens HJ, Bosse M, Paudel Y, et al. Genome sequencing reveals fine scale diversification and reticulation history during speciation in Sus. Genome Biol. 2013;14:R107.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Rothschild MF, Ruvinsky A. The genetics of the pig. 2nd ed. Wallingford: CAB International; 2001.

    Google Scholar 

  7. Darwin C. The variation of animals and plants under domestication. London: John Murray; 1875.

    Google Scholar 

  8. White S. From globalized pig breeds to capitalist pigs: a study in animal cultures and evolutionary history. Environ Hist Durh NC. 2011;16:94–120.

    Article  Google Scholar 

  9. Burgos-Paz W, Souza CA, Megens HJ, Ramayo-Caldas Y, Melo M, Lemus-Flores C, et al. Porcine colonization of the Americas: a 60k SNP story. Heredity (Edinb). 2013;110:321–30.

    Article  CAS  Google Scholar 

  10. Xu W. Introduction and domestication of European breeds of pig in modern China. Anc Mod Agric. 2004;1:54–62.

    Google Scholar 

  11. Gascoigne J. Captain cook: voyager between worlds. London: Hambledon Continuum; 2007.

    Google Scholar 

  12. Noce A, Amills M, Manunza A, Muwanika V, Muhangi D, Aliro T, et al. East African pigs have a complex Indian, Far Eastern and Western ancestry. Anim Genet. 2015;46:433–6.

    Article  CAS  PubMed  Google Scholar 

  13. Chen K, Baxter T, Muir WM, Groenen MA, Schook LB. Genetic resources, genome mapping and evolutionary genomics of the pig (Sus scrofa). Int J Biol Sci. 2007;3:153–65.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Godfray HC, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, et al. Food security: the challenge of feeding 9 billion people. Science. 2010;327:812–8.

    Article  CAS  PubMed  Google Scholar 

  15. Bosse M, Megens HJ, Madsen O, Paudel Y, Frantz LA, Schook LB, et al. Regions of homozygosity in the porcine genome: consequence of demography and the recombination landscape. PLoS Genet. 2012;8:e1003100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Vonholdt BM, Pollinger JP, Lohmueller KE, Han E, Parker HG, Quignon P, et al. Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication. Nature. 2010;464:898–902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Decker JE, McKay SD, Rolf MM, Kim J, Molina Alcala A, Sonstegard TS, et al. Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle. PLoS Genet. 2014;10:e1004254.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491:393–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Rubin CJ, Zody MC, Eriksson J, Meadows JR, Sherwood E, Webster MT, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464:587–91.

    Article  CAS  PubMed  Google Scholar 

  20. Axelsson E, Ratnakumar A, Arendt ML, Maqbool K, Webster MT, Perloski M, et al. The genomic signature of dog domestication reveals adaptation to a starch-rich diet. Nature. 2013;495:360–4.

    Article  CAS  PubMed  Google Scholar 

  21. Carneiro M, Rubin CJ, Di Palma F, Albert FW, Alfoldi J, Barrio AM, et al. Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication. Science. 2014;345:1074–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Rubin CJ, Megens HJ, Martinez Barrio A, Maqbool K, Sayyab S, Schwochow D, et al. Strong signatures of selection in the domestic pig genome. Proc Natl Acad Sci USA. 2012;109:19529–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bosse M, Megens HJ, Frantz LA, Madsen O, Larson G, Paudel Y, et al. Genomic analysis reveals selection for Asian genes in European pigs following human-mediated introgression. Nat Commun. 2014;5:4392.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Wilkinson S, Lu ZH, Megens HJ, Archibald AL, Haley C, Jackson IJ, et al. Signatures of diversifying selection in European pig breeds. PLoS Genet. 2013;9:e1003453.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ai H, Yang B, Li J, Xie X, Chen H, Ren J. Population history and genomic signatures for high-altitude adaptation in Tibetan pigs. BMC Genomics. 2014;15:834.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Ramos AM, Crooijmans RP, Affara NA, Amaral AJ, Archibald AL, et al. Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS One. 2009;4:e6524.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. R Core Team. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.

    Google Scholar 

  29. Felsenstein J. PHYLIP—phylogeny inference package. Cladistics. 1989;5:164–6.

    Google Scholar 

  30. Becker RA, Wilkes AR. Maps in S. AT\&T Bell Laboratories Statistics Research Report. 1993;96:93.2.

  31. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Delaneau O, Marchini J, Zagury JF. A linear complexity phasing method for thousands of genomes. Nat Methods. 2012;9:179–81.

    Article  CAS  Google Scholar 

  33. Amaral AJ, Megens HJ, Crooijmans RP, Heuven HC, Groenen MA. Linkage disequilibrium decay and haplotype block structure in the pig. Genetics. 2008;179:569–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Sved JA. Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theor Popul Biol. 1971;2:125–41.

    Article  CAS  PubMed  Google Scholar 

  35. Tortereau F, Servin B, Frantz L, Megens HJ, Milan D, Rohrer G, et al. A high density recombination map of the pig reveals a correlation between sex-specific recombination and GC content. BMC Genomics. 2012;13:586.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Herrero-Medrano JM, Megens HJ, Groenen MA, Ramis G, Bosse M, Perez-Enciso M, et al. Conservation genomic analysis of domestic and wild pig populations from the Iberian Peninsula. BMC Genet. 2013;14:106.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Cockerham BS, Wa CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:14.

    Google Scholar 

  38. Zhang Z. Pig breeds in China. Shanghai: Shanghai Scientific and Technical Publishers; 1986.

    Google Scholar 

  39. Herrero-Medrano JM, Megens HJ, Crooijmans RP, Abellaneda JM, Ramis G. Farm-by-farm analysis of microsatellite, mtDNA and SNP genotype data reveals inbreeding and crossbreeding as threats to the survival of a native Spanish pig breed. Anim Genet. 2013;44:259–66.

    Article  CAS  PubMed  Google Scholar 

  40. Herrero-Medrano JM, Megens HJ, Groenen MA, Bosse M, Perez-Enciso M, Crooijmans RP. Whole-genome sequence analysis reveals differences in population management and selection of European low-input pig breeds. BMC Genomics. 2014;15:601.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ramirez O, Ojeda A, Tomas A, Gallardo D, Huang LS, Folch JM, et al. Integrating Y-chromosome, mitochondrial, and autosomal data to analyze the origin of pig breeds. Mol Biol Evol. 2009;26:2061–72.

    Article  CAS  PubMed  Google Scholar 

  42. Wagner KR, McPherron AC, Winik N, Lee SJ. Loss of myostatin attenuates severity of muscular dystrophy in mdx mice. Ann Neurol. 2002;52:832–6.

    Article  CAS  PubMed  Google Scholar 

  43. Graham ES, Turnbull Y, Fotheringham P, Nilaweera K, Mercer JG, Morgan PJ, et al. Neuromedin U and Neuromedin U receptor-2 expression in the mouse and rat hypothalamus: effects of nutritional status. J Neurochem. 2003;87:1165–73.

    Article  CAS  PubMed  Google Scholar 

  44. Mankowska M, Szydlowski M, Salamon S, Bartz M, Switonski M. Novel polymorphisms in porcine 3′UTR of the leptin gene, including a rare variant within target sequence for MIR-9 gene in Duroc breed, not associated with production traits. Anim Biotechnol. 2015;26:156–63.

    Article  CAS  PubMed  Google Scholar 

  45. Waraich RS, Weigert C, Kalbacher H, Hennige AM, Lutz SZ, Häring HU, et al. Phosphorylation of Ser357 of rat insulin receptor substrate-1 mediates adverse effects of protein kinase C-delta on insulin action in skeletal muscle cells. J Biol Chem. 2008;283:11226–33.

    Article  CAS  PubMed  Google Scholar 

  46. Wasserman D, Geijer T, Sokolowski M, Rozanov V, Wasserman J. Genetic variation in the hypothalamic-pituitary-adrenocortical axis regulatory factor, T-box 19, and the angry/hostility personality trait. Genes Brain Behav. 2007;6:321–8.

    Article  CAS  PubMed  Google Scholar 

  47. Adachi H, Tsujimoto M, Hattori M, Arai H, Inoue K. cDNA cloning of human cytosolic platelet-activating factor acetylhydrolase gamma-subunit and its mRNA expression in human tissues. Biochem Biophys Res Commun. 1995;214:180–7.

    Article  CAS  PubMed  Google Scholar 

  48. Beneduzzi D, Trarbach EB, Min L, Jorge AA, Garmes HM, Renk AC, et al. Role of gonadotropin-releasing hormone receptor mutations in patients with a wide spectrum of pubertal delay. Fertil Steril. 2014;102:838–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Luo Y, Liu Q, Lei X, Wen Y, Yang YL, Zhang R, et al. Association of estrogen receptor gene polymorphisms with human precocious puberty: a systematic review and meta-analysis. Gynecol Endocrinol. 2015;31:516–21.

    Article  CAS  PubMed  Google Scholar 

  50. Yang WL, Ravatn R, Kudoh K, Alabanza L, Chin KV. Interaction of the regulatory subunit of the cAMP-dependent protein kinase with PATZ1 (ZNF278). Biochem Biophys Res Commun. 2010;391:1318–23.

    Article  CAS  PubMed  Google Scholar 

  51. Young JM, Shykind BM, Lane RP, Tonnes-Priddy L, Ross JA, Walker M, et al. Odorant receptor expressed sequence tags demonstrate olfactory expression of over 400 genes, extensive alternate splicing and unequal expression levels. Genome Biol. 2003;4:R71.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Amaral PP, Neyt C, Wilkins SJ, Askarian-Amiri ME, Sunkin SM, Perkins AC, et al. Complex architecture and regulated expression of the Sox2ot locus during vertebrate development. RNA. 2009;15:2013–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Shen T, Zhu Y, Patel J, Ruan Y, Chen B, Zhao G, et al. T-box20 suppresses oxidized low-density lipoprotein-induced human vascular endothelial cell injury by upregulation of PPAR-gamma. Cell Physiol Biochem. 2013;32:1137–50.

    Article  CAS  PubMed  Google Scholar 

  54. Wang M, Bridges JP, Na CL, Xu Y, Weaver TE. Meckel-Gruber syndrome protein MKS3 is required for endoplasmic reticulum-associated degradation of surfactant protein C. J Biol Chem. 2009;284:33377–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Kaestner KH, Katz J, Liu Y, Drucker DJ, Schutz G. Inactivation of the winged helix transcription factor HNF3alpha affects glucose homeostasis and islet glucagon gene expression in vivo. Genes Dev. 1999;13:495–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Yabe D, Brown MS, Goldstein JL. Insig-2, a second endoplasmic reticulum protein that binds SCAP and blocks export of sterol regulatory element-binding proteins. Proc Natl Acad Sci USA. 2002;99:12753–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Zhang Q, Wang J, Fan S, Wang L, Cao L, Tang K, et al. Expression and functional characterization of LRRC4, a novel brain-specific member of the LRR superfamily. FEBS Lett. 2005;579:3674–82.

    Article  CAS  PubMed  Google Scholar 

  58. Muto E, Tabata Y, Taneda T, Aoki Y, Muto A, Arai K, et al. Identification and characterization of Veph, a novel gene encoding a PH domain-containing protein expressed in the developing central nervous system of vertebrates. Biochimie. 2004;86:523–31.

    Article  CAS  PubMed  Google Scholar 

  59. Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, et al. Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009;459:528–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Tan J, Pieper K, Piccoli L, Abdi A, Foglierini M, Geiger R, et al. A LAIR1 insertion generates broadly reactive antibodies against malaria variant antigens. Nature. 2016;529:105–9.

    Article  CAS  PubMed  Google Scholar 

  61. Onteru SK, Fan B, Du ZQ, Garrick DJ, Stalder KJ, Rothschild MF. A whole-genome association study for pig reproductive traits. Anim Genet. 2012;43:18–26.

    Article  CAS  PubMed  Google Scholar 

  62. Ramirez O, Burgos-Paz W, Casas E, Ballester M, Bianco E, Olalde I, et al. Genome data from a sixteenth century pig illuminate modern breed relationships. Heredity (Edinb). 2015;114:175–84.

    Article  CAS  Google Scholar 

  63. Fang M, Andersson L. Mitochondrial diversity in European and Chinese pigs is consistent with population expansions that occurred prior to domestication. Proc Biol Sci. 2006;273:1803–10.

    Article  PubMed  PubMed Central  Google Scholar 

  64. China National Commission of Animal Genetic Resources. Animal genetic resources in China pigs. Beijing: China Agricultural Press; 2011.

    Google Scholar 

  65. Rothschild M, Jacobson C, Vaske D, Tuggle C, Wang L, et al. The estrogen receptor locus is associated with a major gene influencing litter size in pigs. Proc Natl Acad Sci USA. 1996;93:201–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kijas JM, Wales R, Tornsten A, Chardon P, Moller M, Andersson L. Melanocortin receptor 1 (MC1R) mutations and coat color in pigs. Genetics. 1998;150:1177–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Ojeda A, Huang LS, Ren J, Angiolillo A, Cho IC, Soto H, et al. Selection in the making: a worldwide survey of haplotypic diversity around a causative mutation in porcine IGF2. Genetics. 2008;178:1639–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Murakami K, Yoshikawa S, Konishi S, Ueno Y, Watanabe S, Mizoguchi Y. Evaluation of genetic introgression from domesticated pigs into the Ryukyu wild boar population on Iriomote Island in Japan. Anim Genet. 2014;45:517–23.

    Article  CAS  PubMed  Google Scholar 

  69. Edea Z, Kim SW, Lee KT, Kim TH, Kim KS. Genetic Structure of and evidence for admixture between Western and Korean native pig breeds revealed by single nucleotide polymorphisms. Asian-Aust J Anim Sci. 2014;27:1263–9.

    Article  Google Scholar 

  70. Pham LD, Do DN, Nam LQ, Van Ba N, Minh LT, Hoan TX, et al. Molecular genetic diversity and genetic structure of Vietnamese indigenous pig populations. J Anim Breed Genet. 2014;131:379–86.

    Article  CAS  PubMed  Google Scholar 

Download references

Author’s contributions

LH, MAMG, HJM and BY conceived and coordinated the study. HJM, LH, MAMG, MPE, ZN, RC, LI, MS, SG, PU, PD DD, AD, OH, GS, CK, GL, AA, LBS, ATri, PA and ATra provided the samples and generated the data. LC, BY and HJM analyzed the data. BY, HJM, LH, MPE and MAMG interpreted the results. BY and LC wrote the manuscript. LH, MAMG, MPE and HJM revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Marco Apollonio from Sassari University for contributing to the sampling of Italian wild boars, and Greger Larson from Oxford University for valuable comments on the manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Data available from the Dryad Digital Repository: DOI: doi:10.5061/dryad.30tk6; link: http://dx.doi.org/10.5061/dryad.30tk6.

Ethical statement

All authors declare that animal samples were obtained according to local/national laws at the time of sampling. Exchange of samples and data was in accordance with national and international regulations, and approved by data and sample owners.

Funding

This study is supported by National Production Technology System for the Pig Industry in China (nycytx-008) and Outstanding Talents and Innovation Team of Agricultural Science (2011-81) to LH, AGL2010-14822 and AGL2013-41834-R (Ministry of Economy and Science, Spain) to MPE. LI received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 656697. HJM received funding from the IMAGE project (Horizon 2020, No. 677353).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Bin Yang and Leilei Cui have contributed equally to this work.

Authors and Affiliations

  1. National Key Laboratory for Pig Genetic Improvement and Production Technology, Nanchang, China

    Bin Yang, Leilei Cui & Lusheng Huang

  2. Animal Breeding and Genomics, Wageningen University, Wageningen, The Netherlands

    Richard P. M. A. Crooijmans, Martien A. M. Groenen & Hendrik-Jan Megens

  3. Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB Consortium, Bellaterra, Barcelona, Spain

    Miguel Perez-Enciso

  4. Institut Catala de Recerca i Estudis Avancats (ICREA), Carrer de Lluís Companys, Barcelona, Spain

    Miguel Perez-Enciso

  5. All-Russian Research Institute of Animal Husbandry named after Academy Member L.K. Ernst, Dubrovitzy, Moscow Region, Russia

    Aleksei Traspov & Natalia Zinovieva

  6. Institute of Genomic Biology, University of Illinois, Urbana, Champaign, IL, USA

    Lawrence B. Schook

  7. Division of Genetics and Genomics, The Roslin Institute, R(D)SVS, University of Edinburgh, Edinburgh, UK

    Alan Archibald

  8. Animal and Aquatic Sciences, Chiang Mai University, Chiang Mai, Thailand

    Kesinee Gatphayak

  9. Division of Biotechnology and Reproduction of Livestock, Department of Animal Sciences, Georg-August-University, Göttingen, Germany

    Christophe Knorr

  10. Department of Genetics, Development and Molecular Biology, Aristotle University of Thessaloníki, Thessaloniki, Greece

    Alex Triantafyllidis & Panoraia Alexandri

  11. Research Centre for Biology- Zoology Division, LIPI, Bogor, Indonesia

    Gono Semiadi

  12. School of Biology, University of Nottingham, Notttingham, UK

    Olivier Hanotte

  13. Faculdade de Ciências and CESAM, Universidade de Lisboa, Lisbon, Portugal

    Deodália Dias

  14. Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia

    Peter Dovč

  15. Animal Breeding, Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland

    Pekka Uimari

  16. Department of Chemistry and Bioscience, Aalborg University, Aalborg East, Denmark

    Laura Iacolina

  17. Department of Science for Nature and Environmental Resources, University of Sassari, Sassari, Italy

    Laura Iacolina & Massimo Scandura

Authors
  1. Bin Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Leilei Cui
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Miguel Perez-Enciso
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Aleksei Traspov
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Richard P. M. A. Crooijmans
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Natalia Zinovieva
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lawrence B. Schook
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Alan Archibald
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Kesinee Gatphayak
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Christophe Knorr
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Alex Triantafyllidis
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Panoraia Alexandri
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Gono Semiadi
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Olivier Hanotte
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Deodália Dias
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Peter Dovč
    View author publications

    You can also search for this author in PubMed Google Scholar

  17. Pekka Uimari
    View author publications

    You can also search for this author in PubMed Google Scholar

  18. Laura Iacolina
    View author publications

    You can also search for this author in PubMed Google Scholar

  19. Massimo Scandura
    View author publications

    You can also search for this author in PubMed Google Scholar

  20. Martien A. M. Groenen
    View author publications

    You can also search for this author in PubMed Google Scholar

  21. Lusheng Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  22. Hendrik-Jan Megens
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Lusheng Huang or Hendrik-Jan Megens.

Additional files

12711_2017_345_MOESM1_ESM.xlsx

Additional file 1: Table S1. Detailed information and parameters of population genetics for pig populations in this study.

Additional file 2: Table S2. List of number of individuals and SNPs used in each step of analysis.

Additional file 3: Figure S1. MDS plot for all pig populations with detailed breed information.

Additional file 4: Figure S2. Neighbor-joining tree of pig populations under study.

12711_2017_345_MOESM5_ESM.tif

Additional file 5: Figure S3. Population structure of each population revealed by the ADMIXTURE software at K = 2, 3, 5, 7, 10, 16. We marked the grouping of pig populations on the top of the graph to improve its readability. The first layer of the legend marks the five groups: Asia, Europe, America, Commercial and AF&OC (Africa and Oceania). The second layer of the legend further denotes more specific regions in corresponding continents: ASWB (Asian wild boars), CNDM_W (Chinese western domestic pigs), CNDM_N (Chinese northern domestic pigs), CNDM_S (Chinese southern domestic pigs), ① (Chinese hybrid pigs), ② (Southeast Asian pigs); EUWB (European wild boars), EUDM_S (European southern domestic pigs), UKDM (English domestic pigs), EUDM_N (European northern domestic pigs), RUDM (domestic pigs in Russia and its neighbor countries); AMFR (American feral pigs), AMDM (American domestic pigs); LD (Landrace), LW (Large White), PI (Pietrain), DU (Duroc); AF&OC (Pigs from Africa and Oceania).

12711_2017_345_MOESM6_ESM.tif

Additional file 6: Figure S4. Expanded regional plot of Figure S2 showing the scenario of admixture for pig breeds and populations in Asia.

12711_2017_345_MOESM7_ESM.tif

Additional file 7: Figure S5. Expanded regional plot of Figure S2 showing scenario of admixture for pig breeds and populations in Europe and Russia.

12711_2017_345_MOESM8_ESM.tif

Additional file 8: Figure S6. Expanded regional plot of Figure S2 showing scenario of admixture for pig breeds and populations in North and South America.

12711_2017_345_MOESM9_ESM.tif

Additional file 9: Figure S7. Expanded regional plot of Figure S2 showing scenario of admixture for pig breeds and populations in African and Oceanian countries.

12711_2017_345_MOESM10_ESM.tif

Additional file 10: Figure S8. Scatterplot of the correlation between ROH length and haplotype diversity. Note that the total length of ROH and haplotype diversity are negatively correlated.

12711_2017_345_MOESM11_ESM.tif

Additional file 11: Figure S9. Distribution of haplotype diversity for pig populations across the world. Diamonds and vertical bars represent means and standard deviations of number of haplotypes respectively in 5-SNP s (A), 10-SNP (B), and 15-SNP (C) windows across the genome for each population with a minimum of 10 individuals.

12711_2017_345_MOESM12_ESM.tif

Additional file 12: Figure S10. Comparison of effective population size for domestic pigs and international commercial breeds.

12711_2017_345_MOESM13_ESM.tif

Additional file 13: Figure S11. Observed number of shared SNPs between the top 1% SNPs with the highest Fst values identified in Asia and Europe does not exceed the number of shared SNPs at random. The grey bars show the distribution of number of shared top SNPs at random, the red vertical line represents the observed number of shared top SNPs.

12711_2017_345_MOESM14_ESM.tif

Additional file 14: Figure S12. Results of admixture analysis using 30.4 million SNPs called from whole-genome sequence data (A) were similar to those results obtained using 60K SNP data (B). Whole-genome sequence raw data from 188 individuals mainly obtained from [4, 18, 70], were used in the analysis. The whole-genome SNPs were called using GATK best practice workflow (www.broadinstitute.org/gatk). A total of 30.4 million SNP with a MAF >0.02 and a call rate >70% were kept for admixture analysis (A). A total of 44,988 SNPs with genome positions that were concordant with those of the Illumina 60K SNPs were extracted from the 30.4 million SNP data to represent the results of 60K SNPs (B).

12711_2017_345_MOESM15_ESM.tif

Additional file 15: Figure S13. Scatter plot of geographical distances among pig breeds in China against their genetic distances after removing pigs breeds with more than 20% introgression from European ancestry as revealed by admixture analysis (K = 2).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, B., Cui, L., Perez-Enciso, M. et al. Genome-wide SNP data unveils the globalization of domesticated pigs. Genet Sel Evol 49, 71 (2017). https://doi.org/10.1186/s12711-017-0345-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-017-0345-y

Keywords

Genetics Selection Evolution

ISSN: 1297-9686