The Local Origin of the Tibetan Pig and Additional Insights into the Origin of Asian Pigs

Abstract

Background

The domestic pig currently indigenous to the Tibetan highlands is supposed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. However, there is no direct genetic evidence of either the local origin or exotic migration of the Tibetan pig.

Methods and Findings

We analyzed mtDNA hypervariable segment I (HVI) variation of 218 individuals from seven Tibetan pig populations and 1,737 reported mtDNA sequences from domestic pigs and wild boars across Asia. The Bayesian consensus tree revealed a main haplogroup M and twelve minor haplogroups, which suggested a large number of small scale in situ domestication episodes. In particular, haplogroups D1 and D6 represented two highly divergent lineages in the Tibetan highlands and Island Southeastern Asia, respectively. Network analysis of haplogroup M further revealed one main subhaplogroup M1 and two minor subhaplogroups M2 and M3. Intriguingly, M2 was mainly distributed in Southeastern Asia, suggesting for a local origin. Similar with haplogroup D6, M3 was mainly restricted in Island Southeastern Asia. This pattern suggested that Island Southeastern Asia, but not Southeastern Asia, might be the center of domestication of the so-called Pacific clade (M3 and D6 here) described in previous studies. Diversity gradient analysis of major subhaplogroup M1 suggested three local origins in Southeastern Asia, the middle and downstream regions of the Yangtze River, and the Tibetan highlands, respectively.

Conclusions

We identified two new origin centers for domestic pigs in the Tibetan highlands and in the Island Southeastern Asian region.

About the Authors:

Shuli Yang

Contributed equally to this work with: Shuli Yang, Hao Zhang, Huaming Mao, Dawei Yan, Shaoxiong Lu

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Hao Zhang

Contributed equally to this work with: Shuli Yang, Hao Zhang, Huaming Mao, Dawei Yan, Shaoxiong Lu

Affiliation: College of Animal Science and Technology, China Agricultural University, Beijing, China

Huaming Mao

Contributed equally to this work with: Shuli Yang, Hao Zhang, Huaming Mao, Dawei Yan, Shaoxiong Lu

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Dawei Yan

Contributed equally to this work with: Shuli Yang, Hao Zhang, Huaming Mao, Dawei Yan, Shaoxiong Lu

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Shaoxiong Lu

Contributed equally to this work with: Shuli Yang, Hao Zhang, Huaming Mao, Dawei Yan, Shaoxiong Lu

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Linsheng Lian

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Guiying Zhao

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Yulin Yan

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Weidong Deng

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Xianwei Shi

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Shuxin Han

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Shuai Li

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Xiujuan Wang

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Xiao Gou

* E-mail: [email protected]

Affiliation: College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China

Introduction

Tibetan pigs are mainly distributed in the Tibetan highlands, the largest continuous high elevation ecosystem in the world with an average elevation of more than 4000 m. These animals show striking phenotypic and physiological differences from lowland pigs that have allowed them to adapt to the extreme conditions such as hypoxia. In addition, domestic pig bones discovered in the Tibet archaeological site, Kanuo, that are estimated to be around 5,000 years old indicate that Tibetan people could domesticate the Tibetan pig [1]. This discovery also supports that the Tibetan pig may have a local origin. However, there is no direct genetic evidence of either the local origin or exotic migration of these animals. Generally, the domestic animals currently indigenous to the Tibetan highlands are, with the exception of the yak, believed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. The earliest evidence of human activity and yak husbandry (Qiang Culture) has been found in this region,suggesting that wild yaks (Bos grunniens) were first domesticated in Tibet [2]. Previous studies using mitochondrial DNA (mtDNA) also suggest that Tibetan yak have a local origin [3]. Therefore, to use a comparative mtDNA analysis of the Tibetan pig with relevant populations across Asia we may be able to elucidate the origin of this breed and the degree and pattern of introgression of genetic material from exotic breeds into Tibetan pigs.

The rough genetic relationship between wild boars and domestic pigs has been studied by comparing the mtDNA control region between a relatively small sample of animals [4]–[10]. These earlier reports demonstrate that Asian-type pigs were only recently diverged from each other, whereas they indeed have significant differences from European-type breeds. Furthermore, studies using nuclear DNA markers have also documented a distinct genetic distance between the European and Asian pig populations [11]–[13]. Recently, comprehensive studies have been performed and showed a worldwide phylogeography of wild boars and domestic pigs [14]–[19]. The first systematic study [14] on mtDNA polymorphisms in wild boar and domestic pig populations covering almost the entire globe revealed that pigs were domesticated in at least six different centers. Further studies using appropriately sampled populations in the Island Southeastern Asian [15], Pacific Islands [17], Southern Asian [19], Southeastern Asian [18], [19] and Eastern Asian [18] regions have proposed that the Pacific clade of pigs was recently domesticated within Southeast Asia. Furthermore, Asian pig domestication appears to have occurred mainly in Northeastern India, the Mekong region, and the middle and downstream regions of the Yangtze River. Interestingly, the three centers of domestication are all near the Tibetan highlands.

The specific genetic effects of the Asian centers of origin (listed above) on the Tibetan pig are unclear. In this study, the commonly used hypervariable fragment (HVI) in the mtDNA control region was used to characterize domestic pig and wild boar diversity across a wide range of pig breeds. In addition to 1,714 sequences retrieved from GenBank, 241 novel individuals were sampled from different regions where the samples have not been adequately collected before. These regions contained 12 pig populations indigenous to Yunnan province and the Tibetan highlands of China. Therefore, we here provide the most comprehensive screening of mtDNA variations among domestic pigs and wild boars in Asia to obtain genetic information on the evolutionary domestication history of the Tibetan pig.

Materials and Methods

Sampling and sequencing

Prior approval by the Institutional Animal Care and Use Committee of Yunnan Agricultural University (Contract 2007-0081) was given for all experimental procedures involving animals in the present study. Blood or tissue samples of 241 individuals were collected from pigs indigenous to Yunnan province (n = 102) and the Tibetan highlands (n = 139) of China (Fig. 1). The animals sampled were the Yunnan native pigs (Diannanxiaoer n = 26, Baoshan n = 19, Saba n = 41, Gaoligongshan n = 10, Wujin n = 6), and the Tibetan pigs (Linzhi n = 23, Shannan n = 11, Changdu n = 17, Diqing n = 39, Aba n = 12, Ganzi n = 10, and Hezuo n = 27). DNA was isolated using standard SDS/proteinase K extraction methodology or the Qiagen blood and tissue extraction kit (Qiagen, Valencia, USA). Based on the complete pig mtDNA sequence [20], a primer set (5′-15336CCCAAAGCTGAAATTCTAAA-3′, 5′-15822GGTGAGATGGCCCTGAAGTAAG-3′) was designed to amplify a 431 bp segment of the commonly used mtDNA control region (HVI). PCR products were sequenced directly using the ABI PRISM Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems) and electrophoresed on an ABI 3700 PRISM DNA sequencer (Applied Biosystems). All 241 sequences were submitted to GenBank, with reference numbers HQ148311–HQ148551.

[Figure omitted. See PDF.]

Figure 1. Sampling locations used in the study.

Codes 1–7 represent the Tibetan pig populations of Shannan, Linzhi, Changdu, Diqing, Ganzi, Aba and Hezuo, respectively. Codes 8–12 represent the Yunnan pig populations of Baoshan, Diannanxiaoer, Gaoligongshan, Saba and Wujin, respectively.

https://doi.org/10.1371/journal.pone.0028215.g001

Data analysis

To obtain good coverage of pig breeds, the dataset of 241 samples was complemented with data on 1663 domestic pigs and 292 wild boars. A list of all samples, with GenBank accession number, breed and geographical origin, is provided in an editable format in Supplementary Table S1. The sequences were edited and aligned by the DNAstar package (DNASTAR, Madison, WI), and the polymorphisms in the analyzed segments were exported using MEGA3 [21]. The African warthog (Phacochoerus africanus) is known to be distinct from the Eurasian wild boars and was consequently used as an outgroup for the phylogenetic analyses. Gaps in the aligned sequences were excluded from subsequent analysis. Bayesian consensus (MB) trees were constructed using the Monte Carlo Markov chain (MCMC) approach by Bayes 3.1 [22]. Four independent runs were performed, each with 5 million generations. The posterior probabilities listed on the trees represent the lowest recorded values between all of the runs. The median joining network plots of mtDNA control region sequences [23] were constructed using the Network4.1 program (http://www.fluxus-engineering.com/sharenet.htm).

To test the relationship between the Tibetan pig and other lowland breeds, the distribution of pigs in Asia was partitioned into the 15 geographical regions referred to in previous studies [14], [18]: the Tibetan highlands, South Asia (SA), Yunnan, South East Asia (SEA), Island South East Asia (ISEA), Pacific Islands, the upstream region of the Yangtze River (URYZ), the middle and downstream region of the Yangtze River (MDYZ), South China (SC), Taiwan, the middle and upstream region of the Yellow River (MUYR), the downstream region of the Yellow River (DRYR), North East China (NEC), Japan and Korea.

Results

Sequence variations and haplotypes

We sequenced the mtDNA control region of 241 novel pig samples in the present study, and also downloaded 1,714 Asian pig sequences from GenBank. In total, 229 haplotypes were found using data (Table S1) from this sample of 1,663 domestic pigs and 292 wild boars (Table S2). Of these haplotypes, 162 existed only in domestic pigs, 85 existed only in wild boars, and 18 were shared between both (Table 1). Forty-three haplotypes were identified in the 218 individuals from the seven Tibetan pig populations. The number of haplotypes present in animals from the Tibetan highlands was higher than in animals from any other region (Table S2).

[Figure omitted. See PDF.]

Table 1. Haplotype distribution between domestic pigs and wild boars.

https://doi.org/10.1371/journal.pone.0028215.t001

Phylogenetic analysis of domestic pigs and wild boars in Asia

The tree topology for the 229 haplotypes obtained using the Bayesian (MCMC) phylogenetic method is presented in Figure S1. This Bayesian consensus tree revealed 20 divergent haplogroups, of which eight (labeled D1 through to D8) were found only in domestic specimens, seven (labeled from W1 to W7) were found only in wild boars, four (labeled M, DW1, DW2 and DW3) constituted a mix of wild and domestic samples, and the remaining one haplogroup (termed E) belonged to a European cluster (Figure S1 and Table 1). Most of haplogroups have been defined by high posterior probabilities (above 70%), implying reliable structure among these haplogroups.

The main haplogroup M is the most frequently represented group in Asian breeds with respect to both the number of samples and haplotypes (Tables 1 and S2). Among the 229 haplotypes found in the total sample of animals used in the present study, haplogroup M was represented by 141 haplotypes found in 1,497 domestic pigs and 191 wild boars. By contrast, the 19 minor haplogroups harbored only 88 haplotypes which were found in 166 domestic pigs and 101 wild boars (Tables 1 and S2). However, the D1 individuals were restricted to 15 haplotypes unique to the Hezuo pig of the Tibetan highlands. In addition, the DW1 is an independent group shared between domestic pigs and wild boars from SA, including India, Bhutan and Sri Lanka. Yet, D4 is another independent SA group that exists in Bhutanese pigs. The groups D5 and D7 were mainly found to be distributed in SEA. Interestingly, group D5 is restricted only to Cambodia and Laos, whereas group D7 is found mainly in Myanmar and Thailand. Likewise, groups DW2 and DW3 are also shared between domestic pigs and wild boars in Taiwan and Thailand respectively, whereas group D6 is mainly located in the Pacific Islands and ISEA region. By contrast, groups D2, D3 and D8 are dispersed throughout China. These 12 groups are regionally distributed despite the existence of only a small number of individuals and haplotypes, which suggests that a large number of small scale in situ domestication episodes occurred. Group E animals likely originated from European breeds as this group harbors 9 European haplotypes (representing 32 samples), and they are mainly in the regions of ISEA, SA, MDYZ, Korea and Japan. In addition, the seven wild boar groups also exhibit regional distributions, including one SA group (W1), two Korean groups (W2 and W7), two Japanese groups (W3 and W4), one SEA group (W5), and one Chinese group (W6).

Network analysis of Haplogroup M

Haplogroup M harbors three subgroups M1, M2 and M3. Only the regions containing subgroup M2 or M3 are shown in Figure 2 (the Pacific Islands in A, ISEA in B, Yunnan in C, and SEA in D). Nevertheless, the subgroup M2 in Figure 2D and M3 in Figure 2B exhibit a star-like distribution profile that is typical of exponential population growth. Each of them harbors core haplotypes that have a series of one, two, or greater than two mutation distance derivatives (derived haplotypes) detected in domestic pigs and wild boars. Therefore, taken together, haplotypes within subgroups M2 and M3 might have originated from their core haplotypes as a result of domestication events followed by subsequent expansion.

[Figure omitted. See PDF.]

Figure 2. Median joining networks of haplogroup M in the Pacific Islands (A), ISEA (B), Yunnan (C) and SEA (D) regions.

The haplotypes are symbolized by circles and separated by one substitutional step. The length of each branch is proportional to the number of associated mutations. The size of the circle is proportional to the frequency of the haplotype and the white and black areas of the circle represent the proportion of this frequency made up by domestic pigs and wild boars respectively. Subgroup M2 is identified in B, C and D, and subgroup M3 in A, B, C and D.

https://doi.org/10.1371/journal.pone.0028215.g002

The subgroup M2 is mainly distributed in the SEA region (Figure 2D). Of the 13 haplotypes representing 99 domestic pigs and eight wild boars in this subgroup, 10 haplotypes that accounted for 62 domestic pigs and one haplotype for four wild boars were found in this region. A small number of individuals were also found in such adjacent regions as Yunnan and ISEA (Figures 2C and 2B). This result would support for a SEA origin of the subgroup M2.

The subgroup M3 is mainly found in the ISEA region (Figure 2B). Of the 23 haplotypes identified in 69 domestic pigs and 24 wild boars in this subgroup, 16 haplotypes representing 49 domestic pigs and seven wild boars were found in the ISEA area. In addition, two domestic haplotypes harboring 20 individuals were found in the Pacific Islands (Figure 2A), and 11 wild boar haplotypes were found in the areas of SEA and Yunnan (Figures 2D and 2C). Thus the pattern is consistent with an ISEA origin of the subgroup M3.

The subgroup M1 is distributed across Asia and the Pacific Islands, and harbors 78 domestic haplotypes (representing 1,329 samples) and 34 wild haplotypes (representing 159 samples), between which 11 haplotypes are shared (Table S2). Obviously, subgroup M1 is the major subgroup of the haplogroup M. However, it is difficult to identify any potential domestical center(s) of this subgroup without the analysis of haplotype distributions among the 15 regions.

The divisity gradient analysis of the subgroup M1 between different geographical regions

The subgroup M1 possesses four widely distributed core haplotypes (Hap5, Hap7, Hap8 and Hap9; Table S2), which cover all of Asia and are represented by 820 domestic pigs and 76 wild boars in the current study. The haplotype network (Figure 3) highlights that the four core haplotypes form a parallelogram of one-mutation length, and possess considerable one, two, or greater than two mutation distance derivatives detectable in both domestic pigs and wild boars. The full extent of the diversity in the four core haplotypes within the subgroup M1 is best shown in regions such as MUYR, DRYR, URYZ, MDYZ, SEA, Yunnan, ISEA, Japan and the Tibetan highlands. However, the number of derived or unique haplotypes is the greatest in the SEA, Yunnan, MDYZ and Tibetan highland areas among these nine regions (Figure 3 and Table 2). If we presumably regard the Tibetan highlands, MDYZ and SEA regions as diversity centers, and Yunnan as their crossroads, the diversity of the M1 subgroup would then show a gradual decrease when moving outwards from these centers, whereas it could present a high value in Yunnan. Wild boars and domestic pigs across Asia share this subgroup and possess a roughly similar diversity gradient. Hereby, the result demonstrates that the domestication of the M1 subgroup could mainly occur in the three regions.

[Figure omitted. See PDF.]

Figure 3. Median joining networks of the M1 subgroup.

The haplotypes are symbolized by circles and are separated by one substitutional step. The length of each branch is proportional to the number of associated mutations. The four core haplotypes are indicated with bold lines. The size of the circle is proportional to the frequency of the haplotype and the white and black areas of the circle represent the proportion of this frequency made up by domestic pigs and wild boars. Three large red dots represent the three diversity centers of the Tibetan highlands, Southeastern Asia (SEA) and the middle and down stream region of the Yangtze River (MDYZ). The other 12 regions exhibited diversities that decreased along a gradient (see Table 2 for regional diversity). The abbreviations for the regions are explained in the main text.

https://doi.org/10.1371/journal.pone.0028215.g003

[Figure omitted. See PDF.]

Table 2. Regional distributions of core haplotypes, derived haplotypes and unique haplotypes of domestic pigs and wild boars within subgroup M1.

https://doi.org/10.1371/journal.pone.0028215.t002

Discussion

Small scale pig domestication episodes in Asia

At least six independent lineages are highlighted by the global pattern of mtDNA phylogeography of wild boars and domestic pigs [14]. The present study comprehensively analysing previously reported and novel mtDNA data provides valuable insights into the origins of Asian pigs and the associated small scale domestication episodes.

In Eastern Asia, haplogroups D2, D3 and D8 are dispersed across the Yangtze and Yellow River regions, suggesting that sporadic domestication episodes occurred in these two valley civilizations. Whereas the haplogroup D1 is an independent clade of Tibetan pigs from the highlands in the upstream region of the Yellow River and harbors 15 unique haplotypes. Evidence from ancient DNA also proves that a domestication episode appeared in this area [24]. Additionally, haplogroup DW2 was shared between one Taiwanese domestic pig and one Taiwanese wild boar, suggesting that wild introgression may have occurred.

Phylogenetic analysis of mtDNA haplotypes using domestic pigs and wild boars in the SEA region (Table S1) revealed three haplogroups DW3, D5 and D7, whilst the network analysis identified the subgroup M2. Notably, in this study we described the haplogroups D7 and DW3 for the first time. The haplogroup D5 and the subgroup M2 are identical to haplogoup MTSEA and are present in relatively high frequencies in domestic pigs in the mountainous areas of Cambodia and Laos [19]. Unexpectedly, the haplogroup D7 is mainly found in Myanmar and Thailand, whereas DW3 is shared by domestic pigs and wild boars within the borders of Thailand itself. In addition, the three haplogroups and subgroup M2 are restricted to SEA, indicating small scale domestication episodes in this region.

ISEA sequences recently submitted to GenBank (Table S1) were analyzed in this study. The subgroup M3 revealed by network analysis (Figure 2) and the divergent haplogroup D6 identified by phylogenetic analysis (Figure S1) are mainly restricted to this region. Importantly, the core haplotype Hap1 of the subgroup M3 was shared by four Indonesian wild boars and 26 native domestic pigs, notably, three further Indonesian wild boars possessed sequences that were only one- or two-mutation distance derivatives from Hap1. The subgroup M3 and haplogroup D6 were very similar to the Pacific clade described by Larson et al [14]. This result suggests that ISEA, but not SEA, was the center of domestication of the so-called Pacific clade (M3 and D6 here), although a total of 17 wild boars from the SEA region and Yunnan Province also fell into this subgroup. Thus, the present mtDNA evidence does not support previous results unfortunately that the Pacific clade was recently domesticated within SEA and subsequently dispersed throughout ISEA and Pacific Islands [15], [17], [24].

However, DW1 is an independent haplogroup located in SA and shared between domestic pigs and wild boars, which is also consistent with clade D3 as reported by Larson et al. [14]. Moreover, DW1 was found mainly in sympatric domestic and wild pigs in Northern India, suggesting that an independent domestication event occurred in the Indus Valley Civilization [14], [25]. In addition, haplogroup D4 was found only in domestic pigs in Bhutan, suggesting that a small scale domestication event occurred in the Himalayan foothills.

With the exception of haplogroup E, these 19 minor mtDNA lineages may correspond to distinct small scale domestication events, although the possibility of wild introgression cannot be excluded for those clades represented by only a few individuals.

Domestication center of major subgroup M1 in Asia

The major subgroup M1 harbors nearly half of the haplotypes identified in this study (domestic pigs, 72/162; wild boars, 39/85) and over half of the individuals (domestic pigs, 1329/1663; wild boars, 159/292) (Table S1). This subgroup was represented by samples taken from a range of breeds in all the 15 examined regions. There was considerable overlap in the haplotype status of domestic pig samples from SEA, SC, URYZ, MDYZ, MUYR, Yunnan and Japan, in addition to associations with local wild boars. It is therefore difficult to specify any potential domestication site for subgroup M1 based on limited haplotype distribution or the shared status of domestic pigs and wild boars. The diversity gradient [26] judged by the number of derived haplotypes, the number of unique haplotypes and the proportionate allocation to the four core haplotypes, may have value in furthering our understanding of the origin of this subgroup. The present mtDNA data show a gradient in the diversity of subgroup M1 from high values within the diversity centers (MDYZ, SEA and the Tibetan highlands) to low values in the other regions (Figure 3, Table 2). These three areas may therefore be the site of origin. Subgroup M1 actually contained clades D2 (Eastern Asia) and D5 (Southeastern Asia) and haplogroups D1a2 (Yangtze River), D1b and D1a1a (Mekong region) found previously by the fine grained mtDNA phylogenomic analysis of wild boars and domestic pigs [14], [18]. However, the Tibetan highlands were revealed to be a potential diversity center and origin for the Tibetan pig for the first time in the present study. Great diversity is found in Yunnan which lies at the crossroads between the populations from the three different domestication centers of MDYZ, SEA and Tibet. The Yellow River, URYZ, NEC, SC, ISEA, and Pacific Islands regions, together with Japan, Korea and Taiwan, represent dispersal areas for the initial diversification. The gradual decrease in the diversity from SEA through ISEA and into the Pacific Islands suggests that subgroup M1 pigs in the Pacific Islands were probably domesticated in SEA and subsequently migrated through ISEA and into the Pacific Islands due to human colonization of remote Oceania associated with the Lapita cultural complex. This conclusion was reached without reference to the Pacific clade described by Lum et al [17] and Larson et al [15], [24].

The origin of the Tibetan pig

Fifteen haplotypes (representing 17 individuals) of Tibetan pigs were found in the highly divergent D1 lineage, whilst the other Tibetan pigs belonged to subgroup M1 with the exception of one individual in group D3 and one in D7. The high divergence of the D1 haplogroup and association with one of the diversity centers of subgroup M1 provide two lines of evidence strongly suggesting that Tibetan pigs originated in the Tibetan highlands. Modern Tibetans are thought to be descendants of people who have occupied the Tibetan Plateau since the mid-Holocene (between 7,000 and 5,000 years ago) [27], or potentially even since the late Pleistocene (approximately 21,000 years ago) [28], [29]. 5000 years old domestic pig bones from the Kanuo archaeological site also suggest that Tibetan people domesticated the Tibetan pig [3]. All the evidence strongly support our hypothesis that the Tibetan pig has a local rigin.

Understandably, in most cases, the distribution of domestic animals may be expected to follow the migration of their owners. Previous genetic studies using maternal, paternal and autosomal markers suggest three separate origins of Tibetan people [29]–[34], namely Northern Mongolia, the Yellow River region and Central Asia. The Tibetan plateau is the largest and most geographically isolated area with a high altitude on earth and gradually descends to a low altitude only at its northeastern edge. The Tangbo Ancient Road was a highly important trade route snaking westward from Chang'an, and then across the Tanggula Mountains before reaching Lhasa. Unfortunately, no molecular data are now available to prove the three regions listed above as centers of origin of domestic pigs. Although ancient DNA samples from the genus Sus found from the Yellow River drainage basin also belonged to the most common haplotypes found throughout Asia [24], this could not be used to support the occurrence of a major domestication event in this region. It is suggested that only limited human-mediated dispersal of domestic pigs occurred from the three regions listed above into the Tibetan highlands via the Tangbo Ancient Road. Obviously, the Himalayan mountain range extends from Pakistan in the west to Burma in the east and forms a natural barrier between the Tibetan plateau and the Indian subcontinent. In addition, the diversity of SA region is very low. The genetic contribution from SA to the Tibetan pig may have been highly limited by this natural barrier. Meanwhile, the trade between Tibet and southwestern China was carried out through so called the Tea-horse Ancient Road including the Yunnan-Tibet and the Sichuan-Tibet routes, which together stretch across more than 4,000 kilometers of the Hengduan mountains. A genetic contribution to the Tibetan pig may still have occurred from SEA and the Yangtze River region but, in view of the extreme geographical isolation of the Tibetan highlands, this gene flow from the lowlands cannot readily explain the unexpectedly high diversity in the Tibetan highlands.

In conclusions, our results show that both the substantial genetic diversity and the highly divergent lineage present in the geographically isolated Tibetan highlands strongly support the local domestication of the Tibetan pig. Additionally, it may also be interesting to study the origins of other domestic animals indigenous to the Tibetan highlands in order to further confirm the conclusion.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Bayesian consensus tree of 229 pig mtDNA control region haplotypes. The haplotypes from Hap1 to Hap229 were found in 241 individuals sequenced in this study and 1714 reported sequences (Table S1). Haplotypes existing only in domestic populations, only in wild populations or existing in both were identified. The methods used to generate the tree are discussed in the main text. The digits at the nodes are posterior probabilities. The regional distributions of haplotypes and haplotype sharing between domestic pigs and wild boars within the same haplogroups are presented in Tables S2 and 1.

https://doi.org/10.1371/journal.pone.0028215.s001

(TIF)

Table S1.

Information on sampling of Asian domestic pigs and wild boar.

https://doi.org/10.1371/journal.pone.0028215.s002

(XLS)

Table S2.

Haplotype and geographical distributions across 15 regions.

https://doi.org/10.1371/journal.pone.0028215.s003

(XLS)

Acknowledgments

The authors thank staff of the Animal Science and Veterinary Bureaus of Shangri-la, Linzhi, Aba, Hezuo and Shannan Counties for collections of the Tibetan pig samples for this study.

Author Contributions

Conceived and designed the experiments: XG HZ SY HM DY S. Lu. Performed the experiments: SY HZ SH S. Li XW DY YY GZ. Analyzed the data: XG S. Lu XS LL. Contributed reagents/materials/analysis tools: WD. Wrote the paper: XG HZ SY HM DY S. Lu.

Citation: Yang S, Zhang H, Mao H, Yan D, Lu S, Lian L, et al. (2011) The Local Origin of the Tibetan Pig and Additional Insights into the Origin of Asian Pigs. PLoS ONE6(12): e28215. https://doi.org/10.1371/journal.pone.0028215

References

1. Xu W (2000) The relationship between Han and Tibetan revealed by agricultural origins. Nongyekaogu 1: 248–254.W. Xu2000The relationship between Han and Tibetan revealed by agricultural origins.Nongyekaogu1248254(in Chinese). (in Chinese).

2. Qian YW (1979) The origin of domesticated animals. Beijing: Biohistory Science Press. YW Qian1979The origin of domesticated animalsBeijingBiohistory Science Press

3. Guo S, Savolainen P, Su J, Zhang Q, Qi D, et al. (2006) Origin of mitochondrial DNA diversity of domestic yaks. BMC Evol Biol 6: 73.S. GuoP. SavolainenJ. SuQ. ZhangD. Qi2006Origin of mitochondrial DNA diversity of domestic yaks.BMC Evol Biol673

4. Watanobe T, Ishiguro N, Okumura N, Nakano M, Matsui A, et al. (2001) Ancient mitochondrial DNA reveals the origin of Sus scrofa from Rebun Island, Japan. J Mol Evol 52: 281–289.T. WatanobeN. IshiguroN. OkumuraM. NakanoA. Matsui2001Ancient mitochondrial DNA reveals the origin of Sus scrofa from Rebun Island, Japan.J Mol Evol52281289

5. Watanobe T, Ishiguro N, Nakano M, Takamiya H, Matsui A, et al. (2002) Prehistoric introduction of domestic pigs onto the Okinawa Islands: ancient mitochondrial DNA evidence. J Mol Evol 55: 222–231.T. WatanobeN. IshiguroM. NakanoH. TakamiyaA. Matsui2002Prehistoric introduction of domestic pigs onto the Okinawa Islands: ancient mitochondrial DNA evidence.J Mol Evol55222231

6. Watanobe T, Ishiguro N, Nakano M (2003) Phylogeography and population structure of the Japanese wild boar Sus scrofa leucomystax: mitochondrial DNA variation. Zoolog Sci 20: 1477–1489.T. WatanobeN. IshiguroM. Nakano2003Phylogeography and population structure of the Japanese wild boar Sus scrofa leucomystax: mitochondrial DNA variation.Zoolog Sci2014771489

7. Giuffra E, Kijas JM, Amarger V, Carlborg O, Jeon JT, et al. (2000) The origin of the domestic pig: independent domestication and subsequent introgression. Genetics 154: 1785–1791.E. GiuffraJM KijasV. AmargerO. CarlborgJT Jeon2000The origin of the domestic pig: independent domestication and subsequent introgression.Genetics15417851791

8. Okumura N, Kurosawa Y, Kobayashi E, Watanobe T, Ishiguro N, et al. (2001) Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs. Anim Genet 32: 139–147.N. OkumuraY. KurosawaE. KobayashiT. WatanobeN. Ishiguro2001Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs.Anim Genet32139147

9. Kijas JM, Andersson L (2001) A phylogenetic study of the origin of the domestic pig estimated from the near-complete mtDNA genome. J Mol Evol 52: 302–308.JM KijasL. Andersson2001A phylogenetic study of the origin of the domestic pig estimated from the near-complete mtDNA genome.J Mol Evol52302308

10. Kim KI, Lee JH, Li K, Zhang YP, Lee SS, et al. (2002) Phylogenetic relationships of Asian and European pig breeds determined by mitochondrial DNA D-loop sequence polymorphism. Anim Genet 33: 19–25.KI KimJH LeeK. LiYP ZhangSS Lee2002Phylogenetic relationships of Asian and European pig breeds determined by mitochondrial DNA D-loop sequence polymorphism.Anim Genet331925

11. Fan B, Wang ZG, Li YJ, Zhao XL, Liu B, et al. (2002) Genetic variation analysis within and among Chinese indigenous swine populations using microsatellite markers. Anim Genet 33: 422–427.B. FanZG WangYJ LiXL ZhaoB. Liu2002Genetic variation analysis within and among Chinese indigenous swine populations using microsatellite markers.Anim Genet33422427

12. Yang SL, Wang ZG, Liu B, Zhang GX, Zhao SH, et al. (2003) Genetic variation and relationships of eighteen Chinese indigenous pig breeds. Genet Sel Evol 35: 657–671.SL YangZG WangB. LiuGX ZhangSH Zhao2003Genetic variation and relationships of eighteen Chinese indigenous pig breeds.Genet Sel Evol35657671

13. Fang M, Hu X, Jiang T, Braunschweig M, Hu L, et al. (2005) The phylogeny of Chinese indigenous pig breeds inferred from microsatellite markers. Anim Genet 36: 7–13.M. FangX. HuT. JiangM. BraunschweigL. Hu2005The phylogeny of Chinese indigenous pig breeds inferred from microsatellite markers.Anim Genet36713

14. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, et al. (2005) Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science 307: 1618–1621.G. LarsonK. DobneyU. AlbarellaM. FangE. Matisoo-Smith2005Worldwide phylogeography of wild boar reveals multiple centers of pig domestication.Science30716181621

15. Larson G, Cucchi T, Fujita M, Matisoo-Smith E, Robins J, et al. (2007) Phylogeny and ancient DNA of Sus provides insights into neolithic expansion in Island Southeast Asia and Oceania. Proc Natl Acad Sci U S A 104: 4834–4839.G. LarsonT. CucchiM. FujitaE. Matisoo-SmithJ. Robins2007Phylogeny and ancient DNA of Sus provides insights into neolithic expansion in Island Southeast Asia and Oceania.Proc Natl Acad Sci U S A10448344839

16. Larson G, Albarella U, Dobney K, Rowley Conwy P, Schibler J, et al. Ancient DNA, pig domestication, and the spread of the Neolithic into Europe. Proc Natl Acad Sci U S A 104: 15276–15281.G. LarsonU. AlbarellaK. DobneyP. Rowley ConwyJ. SchiblerAncient DNA, pig domestication, and the spread of the Neolithic into Europe.Proc Natl Acad Sci U S A1041527615281

17. Lum JK, McIntyre JK, Greger DL, Huffman KW, Vilar MG (2006) Recent Southeast Asian domestication and Lapita dispersal of sacred male pseudohermaphroditic ‘tuskers’ and hairless pigs of Vanuatu. Proc Natl Acad Sci U S A 103: 17190–17195.JK LumJK McIntyreDL GregerKW HuffmanMG Vilar2006Recent Southeast Asian domestication and Lapita dispersal of sacred male pseudohermaphroditic ‘tuskers’ and hairless pigs of Vanuatu.Proc Natl Acad Sci U S A1031719017195

18. Wu GS, Yao YG, Qu KX, Ding ZL, Li H, et al. (2007) Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia. Genome Biol 8: R245.GS WuYG YaoKX QuZL DingH. Li2007Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia.Genome Biol8R245

19. Tanaka K, Iwaki Y, Takizawa T, Dorji T, Tshering G, et al. (2008) Mitochondrial diversity of native pigs in the mainland South and South-east Asian countries and its relationships between local wild boars. Anim Sci J 79: 417–434.K. TanakaY. IwakiT. TakizawaT. DorjiG. Tshering2008Mitochondrial diversity of native pigs in the mainland South and South-east Asian countries and its relationships between local wild boars.Anim Sci J79417434

20. Ursing BM, Arnason U (1998) The complete mitochondrial DNA sequence of the pig (Sus scrofa). J Mol Evol 47: 302–306.BM UrsingU. Arnason1998The complete mitochondrial DNA sequence of the pig (Sus scrofa).J Mol Evol47302306

21. Kumar S, Tamura K, Nei M (2004) MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 5: 150–163.S. KumarK. TamuraM. Nei2004MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment.Brief Bioinform5150163

22. Ronquist F, Huelsenbeck JP (2003) MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.F. RonquistJP Huelsenbeck2003MRBAYES 3: Bayesian phylogenetic inference under mixed models.Bioinformatics1915721574

23. Bandelt HJ, Forster P, Sykes BC, Richards MB (1995) Mitochondrial portraits of human populations using median networks. Genetics 141: 743–753.HJ BandeltP. ForsterBC SykesMB Richards1995Mitochondrial portraits of human populations using median networks.Genetics141743753

24. Larson G, Liu RR, Zhao X, Yuan J, Fuller D, et al. (2010) Patterns of East Asian pig domestication, migration, and turnover revealed by modern and ancient DNA. Proc Natl Acad Sci U S A 107: 7686–7691.G. LarsonRR LiuX. ZhaoJ. YuanD. Fuller2010Patterns of East Asian pig domestication, migration, and turnover revealed by modern and ancient DNA.Proc Natl Acad Sci U S A10776867691

25. Groves CP, Grubb P (1993) The Eurasian suids (Sus and Babyrousa): Taxonomy and Description. In: Oliver WLR, editor. Pigs, Peccaries, and Hippos: Status Survey and Conservation Action Plan. Switzerland: IUCN. pp. 107–111.CP GrovesP. Grubb1993The Eurasian suids (Sus and Babyrousa): Taxonomy and Description.WLR OliverPigs, Peccaries, and Hippos: Status Survey and Conservation Action PlanSwitzerlandIUCN107111

26. Pang J, Kluetsch C, Zou X, Zhang A, Luo L, et al. (2009) mtDNA data indicate a single origin for dogs south of Yangtze river, less than 16,300 years ago, from numerous wolves. Mol Biol Evol 26: 2849–2864.J. PangC. KluetschX. ZouA. ZhangL. Luo2009mtDNA data indicate a single origin for dogs south of Yangtze river, less than 16,300 years ago, from numerous wolves.Mol Biol Evol2628492864

27. Su B, Xiao C, Deka R, Seielstad T, Kangwanpong D, et al. (2000) Y chromosome haplotypes reveal prehistoric migrations to the Himalayas. Hum Genet 107: 582–590.B. SuC. XiaoR. DekaT. SeielstadD. Kangwanpong2000Y chromosome haplotypes reveal prehistoric migrations to the Himalayas.Hum Genet107582590

28. Moore LG (2001) Human genetic adaptation to high altitude. High Alt Med Biol 2: 257–279.LG Moore2001Human genetic adaptation to high altitude.High Alt Med Biol2257279

29. Zhao M, Kong Q, Wang H, Peng M, Xie X, et al. (2009) Mitochondrial genome evidence reveals successful Late Paleolithic settlement on the Tibetan Plateau. Proc Natl Acad Sci U S A 106: 21230–21235.M. ZhaoQ. KongH. WangM. PengX. Xie2009Mitochondrial genome evidence reveals successful Late Paleolithic settlement on the Tibetan Plateau.Proc Natl Acad Sci U S A1062123021235

30. Chu JY, Huang W, Kuang SQ, Wang JM, Xu JJ, et al. (1998) Genetic relationship of populations in China. Proc Natl Acad Sci U S A 95: 11763–11768.JY ChuW. HuangSQ KuangJM WangJJ Xu1998Genetic relationship of populations in China.Proc Natl Acad Sci U S A951176311768

31. Nie M, Roychoudhury AK (1993) Evolutionary relationships of human populations on a global scale. Mol Biol Evol 10: 927–943.M. NieAK Roychoudhury1993Evolutionary relationships of human populations on a global scale.Mol Biol Evol10927943

32. Torroni A, Miller JA, Moore LG, Zamudio S, Zhuang J, et al. (1994) Mitochondrial DNA analysis in Tibet: implications for the origin of the Tibetan population and its adaptation to high altitude. Am J Phys Anthropol 93: 189–199.A. TorroniJA MillerLG MooreS. ZamudioJ. Zhuang1994Mitochondrial DNA analysis in Tibet: implications for the origin of the Tibetan population and its adaptation to high altitude.Am J Phys Anthropol93189199

33. Qian Y, Qian B, Su B, Yu J, Ke Y, et al. (2000) Multiple origins of Tibetan Y chromosomes. Hum Genet 106: 453–454.Y. QianB. QianB. SuJ. YuY. Ke2000Multiple origins of Tibetan Y chromosomes.Hum Genet106453454

34. Van Driem G (2005) Tibeto-Burman vs Indo-Chinese: implications for population geneticists, archeologists and prehistorians. In: Sagart L, Blench R, Sanches-Mazas A, editors. The peopling of East Asia: putting together archaeology, linguistics and genetics. London: Routledge Curzon. pp. 81–106.G. Van Driem2005Tibeto-Burman vs Indo-Chinese: implications for population geneticists, archeologists and prehistorians.L. SagartR. BlenchA. Sanches-MazasThe peopling of East Asia: putting together archaeology, linguistics and geneticsLondonRoutledge Curzon81106