Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia

Background  

Previously reported evidence indicates that pigs were independently domesticated in multiple places throughout the world. However, a detailed picture of the origin and dispersal of domestic pigs in East Asia has not yet been reported.     

De novo assembly of mitochondrial genomes provides insights into genetic diversity and molecular evolution in wild boars and domestic pigs

Up to date, the scarcity of publicly available complete mitochondrial sequences for European wild pigs hampers deeper understanding about the genetic changes following domestication. Here, we have assembled 26 de novo mtDNA sequences of European wild boars from next generation sequencing (NGS) data and downloaded 174 complete mtDNA sequences to assess the genetic relationship, nucleotide diversity, and selection. The Bayesian consensus tree reveals the clear divergence between the European and Asian clade and a very small portion (10 out of 200 samples) of maternal introgression. The overall nucleotides diversities of the mtDNA sequences have been reduced following domestication. Interestingly, the selection efficiencies in both European and Asian domestic pigs are reduced, probably caused by changes in both selection constraints and maternal population size following domestication. This study suggests that de novo assembled mitogenomes can be a great boon to uncover the genetic turnover following domestication. Further investigation is warranted to include more samples from the ever-increasing amounts of NGS data to help us to better understand the process of domestication.     

Mitochondrial DNA sequence variation and phylogenetic relationships among Iberian pigs and other domestic and wild pig populations

Nucleotide sequences of mitochondrial DNA (mtDNA) cytochrome B gene (1140 bp) and control region (707 bp) were used to determine the phylogenetic relationships among 51 pig samples representing ancient and current varieties of Iberian pigs (26), Spanish wild boars (seven) and other domestic pigs (18) of cosmopolitan (Duroc, Large White, Landrace, Pietrain and Meishan) and local (Spotted Black Jabugo, Basque and Mangalitza) breeds. A neighbour-joining tree constructed from pairwise distances provide evidence of the European origin of both Iberian pigs and Spanish wild boars. The introgression of Asian mtDNA haplotypes in the genetic pool of the Iberian breed seems unlikely. Four estimates of sequence divergence between European and Asian clades were calculated from the two main domains of the D-loop region and the synonymous and nonsynonymous nucleotide substitutions in the cytochrome B gene. The time since the divergence of pig ancestors was estimated at about 600,000 years before present.     

Identification of mitochondrial markers for genetic traceability of European wild boars and Iberian and Duroc pigs

Iberian pigs and wild boars are the source of highly priced meat and dry-cured products. Iberian maternal origin is mandatory for labeled Iberian products, making necessary the authentication of their maternal breed origin. Discrimination between wild and domestic pig maternal origin may be useful to distinguish labeled wild boar meat obtained from hunting or farming. In order to detect useful polymorphisms to trace Iberian, Duroc and wild boar maternal lineages, we herein investigated the complete porcine mitochondrial DNA (mtDNA) using three complementary approaches. Near-complete mtDNA sequences (16989 bp), excluding the minisatellite present in the displacement loop region (D-loop), were successfully determined in six Iberian pigs, two Duroc and six European wild boars. To complete the mtDNA analysis, the D-loop minisatellite region was also analyzed in the same set of samples by amplification and capillary electrophoresis detection. Finally, the frequencies of Asian and European Cytochrome B (Cyt B) haplotypes were estimated in Iberian (n = 96) and Duroc (n = 125) breeds. Comparison of near-complete mtDNA sequences revealed a total of 57 substitutions and two Indels. Out of them, 32 polymorphisms were potential Iberian markers, 10 potential Duroc markers and 16 potential wild boar markers. Fourteen potential markers (five Iberian and nine Duroc), were selected to be genotyped in 96 Iberian and 91 Duroc samples. Five wild boar potential markers were selected and tested in samples of wild boars (73) and domestic pigs including: 96 Iberian, 16 Duroc, 16 Large White and 16 Landrace. Genotyping results showed three linked markers (m.7998C>T, m.9111T>C, m.14719A>G) absent in Duroc and present in Iberian pigs with a frequency 0.72. Six markers (m.8158C>T, m.8297T>C, m.9230G>A, m.11859A>G, m.13955T>C, m.16933T>C), three of them linked, were absent in Iberian pigs and present in Duroc with a joint frequency of almost 0.50. Finally three linked markers (m.7188G>A, m.9224T>C, m.15823A>G) were solely detected in wild boars with a frequency 0.22. The D-loop minisatellite results showed overlapping ranges of fragment sizes and suggested heteroplasmy, a result that nullify the use of this region for the development of breed diagnostic markers. The Cyt B haplotype results showed the presence of European haplotypes in Iberian while one of the Asian haplotypes was detected in Duroc with a frequency 0.22, linked to the Duroc marker m.9230G>A. Our results are valuable to resolve the problems of Iberian and wild boar maternal origin determination but additional markers are required to achieve totally useful genetic tests.     

Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs

We completed phylogenetic analysis of the major non-coding region of the mitochondrial DNA (mtDNA) from 159 animals of eight Euro-American and six East Asian domesticated pig breeds and 164 Japanese and five European wild boars. A total of 62 mtDNA haplotypes were detected. Alignment of these regions revealed nucleotide variations (including gaps) at 73 positions, including 58 sites with transition nucleotide substitutions, and two transversion substitutions. Phylogenetic analysis of the sequences could not organize domestic pig breeds into discrete clusters. In addition, many of the haplotypes found in members of diverged clustering groups were found primarily in Euro-American pig breeds, indicating extensive introgression of Asian domestic pigs into European breeds. Furthermore, phylogenetic analysis allocated the DNA sequences of non-coding regions into two different groups, and the deepest branchpoint of this porcine phylogeny corresponded to 86 000-136 000 years before present. This time of divergence would predate the historical period when the pig is thought to have been domesticated from the wild boar.     

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.     

Relationship between wild greylag and European domestic geese based on mitochondrial DNA

The origins of the European domestic goose are uncertain. The available information comes from archaeological findings and historical literature, but genetic evidence has hitherto been scarce. The domestic goose in Europe is derived from the greylag goose (Anser anser), but it is not known where the initial domestication took place and which of the two subspecies of greylag goose was ancestral. We aimed to determine the amount and geographical distribution of genetic diversity in modern populations of greylag geese as well as in different breeds of the domestic goose to make inferences about goose domestication. We studied DNA sequence variation in the mitochondrial control region of greylag geese from multiple populations across Europe and western Asia as well as specimens of domestic geese representing 18 modern breeds and individuals not belonging to any recognised breed. Our results show notable differences in genetic diversity between different greylag goose populations and the presence of six mitochondrial haplogroups which show a degree of geographical partitioning. The genetic diversity of the domestic goose is low, with 84% of sampled individuals having one of two major closely related haplotypes, suggesting that modern European domestic geese may derive from a narrow genetic base. The site of domestication remains unresolved, but domestic geese in Turkey were unusually diverse, indicating the importance of further sampling in the vicinity of the eastern Mediterranean and the Near East. There appears to be past or ongoing hybridisation between greylags and domestic geese in particular areas, consistent with field observations.     

Multi-locus analysis reveals a different pattern of genetic diversity for mitochondrial and nuclear DNA between wild and domestic pigs in East Asia

Background

A major reduction of genetic diversity in mtDNA occurred during the domestication of East Asian pigs. However, the extent to which genetic diversity has been lost in the nuclear genome is uncertain. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationships and demographic history of domestic pigs and their ancestors, wild boars, we investigated 14 nuclear markers (including 8 functional genes, 2 pseudogenes and 4 intergenic regions) from 11 different chromosomes in East Asia-wide samples and pooled them with previously obtained mtDNA data for a combined analysis.

Principal Findings

The results indicated that domestic pigs and wild boars possess comparable levels of nucleotide diversity across the nuclear genome, which is inconsistent with patterns that have been found in mitochondrial genome.

Conclusions

This incongruence between the mtDNA and nuclear genomes is suggestive of a large-scale backcross between male wild boars and female domestic pigs in East Asia. Our data reveal the impacts of founder effects and backcross on the pig genome and help us better understand the complex demographic histories of East Asian pigs, which will be useful for future work on artificial selection.     

Ancient DNA reveals the maternal genetic history of East Asian domestic pigs

Zoo-archaeological and genetic evidence suggest that pigs were domesticated independently in Central China and Eastern Anatolia along with the development of agricultural communities and civilizations. However, the genetic history of domestic pigs, especially in China, has not been fully explored. In this study, we generate 42 complete mitochondrial DNA sequences from ～7500- to 2750-year-old individuals from the Yellow River basin. Our results show that the maternal genetic continuity of East Asian domestic pigs dates back to at least the Early to Middle Neolithic. In contrast, the Near Eastern ancestry in European domestic pigs saw a near-complete genomic replacement by the European wild boar. The majority of East Asian domestic pigs share close haplotypes, and the most recent common ancestor of most branches dates back to less than 20,000 years before present, inferred using new substitution rates of whole mitogenomes or combined protein-coding regions. Two major population expansion events of East Asian domestic pigs coincided with changes in climate, widespread adoption of introduced crops, and the development of agrarian societies. These findings add to our understanding of the maternal genetic composition and help to complete the picture of domestic pig evolutionary history in East Asia.     

Genetic evidence for introgression between domestic pigs and wild boars (Sus scrofa) in Belgium and Luxembourg: a comparative approach with multiple marker systems

Hybridization between wild species and their domestic relatives can be an important conservation and management problem. Genetic purity of the wild species is desirable per se and the phenomenon can have unpredictable evolutionary consequences. Declining European wild boar populations were frequently restocked with farmed wild boars that sometimes had been crossed with domestic pigs. We used simple polymerase chain reaction‐based diagnostic tests to detect the presence of mitochondrial DNA and coat colour alleles of domestic origin in wild boars from Belgium, Luxembourg, and western Germany. Microsatellite genotypes were used to test for genetic admixture between the wild boars and domestic pigs. Although almost one‐third of all Luxembourg wild boars carried Asian mitochondrial DNA haplotypes originating from domestic pigs, microsatellite‐based clustering only identified four putatively admixed individuals in Luxembourg. By contrast, clustering identified wild boar × domestic hybrids in most sampling locations in Belgium. We interpret these results as evidence of releases of hybrid captive‐reared wild boars. Our results emphasize the need (if working with classical markers) to use different systems to obtain an understanding as to whether hybridization between wild and domestic relatives might have affected the genetic make‐up of a local population. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 104–115.     

Phylogeographical analyses of domestic and wild yaks based on mitochondrial DNA: new data and reappraisal

Aim We aimed to examine the phylogeographical structure and demographic history of domestic and wild yaks (Bos grunniens) based on a wide range of samples and complete mitochondrial genomic sequences. Location The Qinghai‐Tibetan Plateau (QTP) of western China. Methods All available D‐loop sequences for 405 domesticated yaks and 47 wild yaks were examined, including new sequences from 96 domestic and 34 wild yaks. We further sequenced the complete mitochondrial genomes of 48 domesticated and 21 wild yaks. Phylogeographical analyses were performed using the mitochondrial D‐loop and the total genome datasets. Results We recovered a total of 123 haplotypes based on the D‐loop sequences in wild and domestic yaks. Phylogenetic analyses of this dataset and the mitochondrial genome data suggested three well‐supported and divergent lineages. Two lineages with six D‐loop haplogroups were recovered for all morphological breeds of domestic yaks across their distributions in the QTP, while one more lineage and more endemic haplogroups or haplotypes were found for wild yaks. Based on the mitochondrial genome data, the divergences of the three lineages were estimated to have occurred around 420,000 and 580,000 years ago, consistent with the geological records of two large glaciation events experienced in the QTP. Main conclusions There are distinct phylogeographical differences between wild and domestic yaks. However, there is no apparent geographical correlation between identified haplogroups and distributions of domestic yaks. Three differentiated lineages of yaks probably evolved allopatrically in different regions during the Pleistocene glaciation events, then reunited into a single gene pool during post‐glacial population expansion and migrations before the start of the domestication of yaks in the Holocene.     

Photoperiod and luteinizing hormone secretion in domestic and wild pigs

During seasonal anoestrus (long-days), oestradiol can effectively inhibit the pulsatile secretion of luteinizing hormone (LH) in sheep. The aim of our trial was to determine whether the same regulatory mechanism exists in the pig. Altogether, 20 ovariectomized and oestradiol-implanted gilts (16 domestic pigs, 4 European wild boars) were randomly allocated to two treatment groups. The first group was kept under a short-day light-dark cycle of 8L:16D, and the second group under a long-day light regime of 16L:8D. After a 6-week treatment period, blood samples were taken at 20-min intervals for 12h. After sampling, the light regimens were switched. Sampling was then repeated following another 6 weeks of treatment. In both treatment groups, 2.3 LH pulses occurred every 12h. The basal LH level was 0.7+/-0.4 ng/ml for the short-day group and 1.0+/-0.5 ng/ml for the long-day group. The mean LH level was 0.9+/-0.4 and 1.3+/-0.6 ng/ml and the LH pulse amplitude 0.5+/-0.4 and 0.6+/-0.5 ng/ml, respectively. The basal and mean LH levels were therefore lower in short-day gilts (P<0.05), while LH pulse amplitude and frequency remained unaffected by treatment. In conclusion, the 6-week period under two different light regimes resulted in higher basal LH concentration in long-day gilts but was not able to produce changes in LH frequency in prepubertal gilts.     

Mitochondrial haplotypes of European wild boars with 2n = 36 are closely related to those of European domestic pigs with 2n = 38

Wild boars from Western Europe have a 2n = 36 karyotype, in contrast to a karyotype of 2n = 38 in wild boars from Central Europe and Asia and in all domestic pigs. The phylogenetic status of this wild boar population is unclear, and it is not known if it has contributed to pig domestication. We have now sequenced the mtDNA control region from 30 European wild boars (22 with a confirmed 2n = 36 karyotype) and six Asian wild boars (two Hainan and four Dongbei wild boars) to address this question. The results revealed a close genetic relationship between mtDNA haplotypes from wild boars with 2n = 36 to those from domestic pigs with 2n = 38. Thus, we cannot exclude the possibility that wild boars with 2n = 36 may have contributed to pig domestication despite the karyotype difference. One of the European wild boars carried an Asian mtDNA haplotype, and this most likely reflects gene flow from domestic pigs to European wild boars. However, this gene flow does not appear to be extensive because the frequency of Asian haplotypes detected among European wild boars (c. 3%) were 10-fold lower than among European domestic pigs (c. 30%). Previous studies of mtDNA haplotypes have indicated that pig populations in Europe and Asia have experienced a population expansion, but it is not clear if the expansion occurred before or after domestication. The results of the present study are consistent with an expansion that primarily occurred prior to domestication because the mtDNA haplotypes found in European and Asian wild boars did not form their own clusters but were intermingled with haplotypes found in domestic pigs, indicating that they originated from the same population expansion.     

Phylogenetic relationships of Asian and European pig breeds determined by mitochondrial DNA D-loop sequence polymorphism

Phylogenetic relationships among Asian and European pig breeds were assessed using 1036 bp of mitochondrial DNA (mtDNA) D-loop sequences. An unweighted pair-group method with arithmetic mean (UPGMA) tree was constructed on the basis of maximum likelihood distances using sequences determined for three Cheju (Korea), 11 Chinese, one Westran (Australian feral origin) and two European pigs (Berkshire and Welsh), and also published sequences for four Japanese (including two Wild Boars), one Yucatan miniature, five European (including Large White, Landrace, Duroc, Swedish and Wild Boar) and two Meishan pigs. The Colombian collared peccary (Tayassu tajacu) sequence was also determined and used as an outgroup. The maximum parsimony with heuristic search method was used to determine bootstrap support values. Asian-type pigs clustered together (bootstrap support 33%), but were separate from European-type pigs that also clustered together (93%). The Westran pig, derived from the feral descendants of pigs inhabiting Kangaroo Island of South Australia, clustered with Asian pigs, demonstrating Asian origin of their mitochondria. Berkshire and Large White clustered with Asian pigs, indicating that Asian pigs were involved in the development of these breeds. Our findings clearly demonstrate that pigs indigenous to China, Korea and Japan are only recently diverged from each other and distinctly different from European-type pigs. European pig breeds consist of pigs with mitochondria of Asian and non-Asian type, some of which were formed from closely related maternal ancestors, if not from a single ancestor.     

Geographic population structure and sequence divergence in the mitochondrial DNA control region of the Japanese wild boar (Sus scrofa leucomystax), with reference to those of domestic pigs

Mitochondrial DNA (mtDNA) control regions from 40 Japanese wild boars were examined by direct sequencing after amplification by PCR. From the DNA sequences obtained, we found eight haplotypes, whose differences arose via transitions. The geographical distribution of these different haplotypes indicated that wild boar populations inhabited limited areas and that there was some restricted gene flow between local populations. Eight mtDNA haplotypes from Eastern and Western domestic pigs and the Ryukyu wild boar were also analyzed as references to those from Japanese wild boars. The cluster analyses of the control-region sequences showed that those from Japanese wild boards belong to the Asian type as do those from Eastern domestic pigs and the Ryukyu wild boar, which differed from the European type (Western domestic pigs).     

Introgressive hybridisation between domestic pigs (Sus scrofa domesticus) and endemic Corsican wild boars (S. s. meridionalis): effects of human-mediated interventions

Owing to the intensified domestication process with artificial trait selection, introgressive hybridisation between domestic and wild species poses a management problem. Traditional free-range livestock husbandry, as practiced in Corsica and Sardinia, is known to facilitate hybridisation between wild boars and domestic pigs (Sus scrofa). Here, we assessed the genetic distinctness and genome-wide domestic pig ancestry levels of the Corsican wild boar subspecies S. s. meridionalis, with reference to its Sardinian conspecifics, employing a genome-wide single nucleotide polymorphism (SNP) assay and mitochondrial control region (mtCR) haplotypes. We also assessed the reliance of morphological criteria and the melanocortin-1 receptor (MC1R) coat colour gene to identify individuals with domestic introgression. While Corsican wild boars showed closest affinity to Sardinian and Italian wild boars compared to other European populations based on principal component analysis, the observation of previously undescribed mtCR haplotypes and high levels of nuclear divergence (Weir’s θ > 0.14) highlighted the genetic distinctness of Corsican S. s. meridionalis. Across three complementary analyses of mixed ancestry (i.e., STRUCTURE, PCADMIX, and ELAI), proportions of domestic pig ancestry were estimated at 9.5% in Corsican wild boars, which was significantly higher than in wild boars in Sardinia, where free-range pig keeping was banned in 2012. Comparison of morphologically pure- and hybrid-looking Corsican wild boars suggested a weak correlation between morphological criteria and genome-wide domestic pig ancestry. The study highlights the usefulness of molecular markers to assess the direct impacts of management practices on gene flow between domestic and wild species.Subject terms: Conservation genomics, Genetic hybridization, Structural variation     

European wild boars and domestic pigs display different polymorphic patterns in the Toll-like receptor (TLR) 1, TLR2, and TLR6 genes

During the last decade, the Toll-like receptors (TLRs) have been extensively studied, and their immense importance in innate immunity is now being unveiled. Here, we report pronounced differences—probably reflecting the domestication process and differences in selective pressure—between wild boars and domestic pigs regarding single nucleotide polymorphisms (SNPs) in TLR genes. The open reading frames of TLR1, TLR2, and TLR6 were sequenced in 25 wild boars, representing three populations, and in 15 unrelated domestic pigs of Hampshire, Landrace, and Large White origin. In total, 20, 27, and 26 SNPs were detected in TLR1, TLR2, and TLR6, respectively. In TLR1 and TLR2, the numbers of SNPs detected were significantly lower (P?≤?0.05, P?≤?0.01) in the wild boars than in the domestic pigs. In the wild boars, one major high frequency haplotype was found in all three genes, while the same pattern was exhibited only by TLR2 in the domestic pigs. The relative frequency of non-synonymous (dN) and synonymous (dS) SNPs was lower for the wild boars than for the domestic pigs in all three genes. In addition, differences in diversity between the genes were revealed: the mean heterozygosity at the polymorphic positions was markedly lower in TLR2 than in TLR1 and TLR6. Because of its localization—in proximity of the bound ligand—one of the non-synonymous SNPs detected in TLR6 may represent species-specific function on the protein level. Furthermore, the codon usage pattern in the genes studied deviated from the general codon usage pattern in Sus scrofa.