Patrilineal Perspective on the Austronesian Diffusion in Mainland Southeast Asia

Abstract

The Cham people are the major Austronesian speakers of Mainland Southeast Asia (MSEA) and the reconstruction of the Cham population history can provide insights into their diffusion. In this study, we analyzed non-recombining region of the Y chromosome markers of 177 unrelated males from four populations in MSEA, including 59 Cham, 76 Kinh, 25 Lao, and 17 Thai individuals. Incorporating published data from mitochondrial DNA (mtDNA), our results indicated that, in general, the Chams are an indigenous Southeast Asian population. The origin of the Cham people involves the genetic admixture of the Austronesian immigrants from Island Southeast Asia (ISEA) with the local populations in MSEA. Discordance between the overall patterns of Y chromosome and mtDNA in the Chams is evidenced by the presence of some Y chromosome lineages that prevail in South Asians. Our results suggest that male-mediated dispersals via the spread of religions and business trade might play an important role in shaping the patrilineal gene pool of the Cham people.

About the Authors:

Jun-Dong He

Contributed equally to this work with: Jun-Dong He, Min-Sheng Peng

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China, School of Life Sciences, University of Science and Technology of China, Hefei, China, KIZ/CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming, China

Min-Sheng Peng

Contributed equally to this work with: Jun-Dong He, Min-Sheng Peng

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China, KIZ/CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming, China, Graduate University of the Chinese Academy of Sciences, Beijing, China

Huy Ho Quang

Affiliation: Department of Immunology, Hanoi Medical University, Hanoi, Vietnam

Khoa Pham Dang

Affiliation: Department of Immunology, Hanoi Medical University, Hanoi, Vietnam

An Vu Trieu

Affiliation: Department of Immunology, Hanoi Medical University, Hanoi, Vietnam

Shi-Fang Wu

Affiliation: State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

Jie-Qiong Jin

Affiliation: State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

Robert W. Murphy

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China, Centre for Biodiversity and Conservation Biology, Royal Ontario Museum, Toronto, Canada

Yong-Gang Yao

Affiliation: Key Laboratory of Animal Models and Human Disease Mechanisms, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

Ya-Ping Zhang

* E-mail: [email protected]

Affiliations State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China, KIZ/CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming, China, Laboratory for Conservation and Utilization of Bio-resource, Yunnan University, Kunming, China

Introduction

The Austronesian language family is one of the largest and most widespread language families. It is spoken by more than 350 million people on islands from Madagascar to Easter Island [1], [2]. Nevertheless, the languages in this family have a rather limited distribution on the mainland. Chamic, the representative language of the family, is spoken by the Cham people. In Mainland Southeast Asia (MSEA), Chamic exists as a “linguistic enclave”, because it is surrounded by non-Austronesian-speaking groups (e.g. Mon-Khmers) [3], [4], [5]. Many studies investigate the diffusion of Austronesian in MSEA by tracing the origin of the Cham people. The “Out-of-Taiwan” hypothesis regards the Cham ancestors as the Austronesian immigrants from Island Southeast Asia (ISEA) and immigration is dated to around 500 BC [4], [6], [7]. Before the arrival of the Austronesian immigrants, southern Vietnam appears to have been occupied by the local Austro-Asiatic speakers, especially Mon-Khmers [8]. There is a high chance of admixture between the Chams and Mon-Khmer groups. Previously linguistic analyses of the Chamic report that some loan-words from Mon-Khmer languages form indigenous cultural contributions [4], [6]. The “Nusantao Maritime Trading and Communication Networks” hypothesis states that cultural diffusion through trading and communication networks played an important or even dominant role in the ethnogenesis of the Cham [9]. Because the origin of the Cham people is open to debate, the demographic history of the Austronesians in Southeast Asia requires further investigation.

Analyses of mitochondrial DNA (mtDNA) variation of the Cham population resolve a closer relationship with populations in MSEA rather than with those from ISEA, and this occurs despite that recent gene flow from ISEA [10]. This result suggests that the origin of the Cham people likely involves the massive assimilation of local Mon-Khmer populations, and this is accompanied with language shift. Thus, the Austronesian diffusion in MSEA appears to mediated mainly by cultural diffusion [10]. Because mtDNA data only offer a maternal perspective, only half of the story is known. Does patrilineal history reveal the same story? We address this question by evaluating non-recombining region of the Y chromosome (NRY) markers, including 26 single-nucleotide polymorphisms (Y-SNPs) and eight short tandem repeats (Y-STRs), in 59 male Cham individuals whose matrilineal histories are known [10]. For comparison, the NRY markers of 76 Kinh, 25 Lao, and 17 Thai males were also surveyed (Figure 1; Table 1).

[Figure omitted. See PDF.]

Figure 1. Populations from southern China and Southeast Asia analyzed in this study.

Diamonds denote the location of populations newly sampled in this study; populations reported are indicated by bold circles. Population ID numbers are shown in Table 1. MSEA, Mainland Southeast Asia; WISEA, western Island Southeast Asia; EISEA, eastern Island Southeast Asia.

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

[Figure omitted. See PDF.]

Table 1. General information for 57 populations in southern China and Southeast Asia.

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

Results

Phylogeny of Y chromosomes

Based on 26 Y-SNPs, all 177 newly genotyped males from the four populations were assigned to specific (sub-)haplogroups (paragroups) defined in the phylogeny (Figure 2; Table S1). Nearly 60% of the Chams' Y chromosomes belonged to P191-derived haplogroups. Within this group, O-M95* predominated and accounted for around 30% of all samples. Haplogroup C-M216, consisting of C-M217 and C-M216*, comprised 10.2% of the patrilineal lineages. One Cham individual (∼1.7%) rooted near the base of the tree as haplogroup F-M213* and six individuals (∼10.2%) rooted at the base of the tree as haplogroup K-P131*. Notably, South Asian-prevailing haplogroups R-M17 (∼13.6%), R-M124 (∼3.4%), and H-M69 (∼1.7%) are identified with the Chams.

[Figure omitted. See PDF.]

Figure 2. Classification tree of 26 NRY haplogroups along with their frequencies (%) in four populations.

Haplogroup defining markers are given along the branches; corresponding markers genotyped in Karafet et al. [27] are noted in brackets. The names of haplogroups are shown to the right of the branches using the mutation-based nomenclature of the Karafet et al. [35].

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

[Figure omitted. See PDF.]

Figure 3. PCA plot based on NRY haplogroup frequencies of 45 populations in southern China and Southeast Asia.

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

[Figure omitted. See PDF.]

Figure 4. MDS plot of 53 populations with RST genetic distances based on eight common Y-STRs.

For population information, see Table 1. Because of severe genetic drift [27], populations Taiwan, Nias, and Mentawai that were resolved as the outliers in the initial analyses and were excluded.

https://doi.org/10.1371/journal.pone.0036437.g004

Population structure

Genetic relationships between the Cham and other Southeast Asian populations were discerned with the aid of additional published Y-chromosomal datasets (Figure 1; Table 1). We employed a principal component analysis (PCA) based on the NRY haplogroup distribution frequencies of 45 populations (Table S2) to show the overall clustering pattern of the populations. Populations from eastern ISEA (EISEA) and from Laos formed two clusters in the first PC (Figure 3) and this pattern was mainly owed to haplogroups C-M216, K-P131*, and O-M95* (Figure S1). The second PC resolved a close affinity between the Kinh and Vietnamese (most likely, the Kinh) populations with those from mainland southern China due to the high frequency of haplogroup O-M88 (Figure S1). The Cham population showed a close affinity to some but not all populations from western ISEA (WISEA; Figure 3). The clustering pattern revealed by PC1 and PC2 was statistically significant (P<0.05) in AMOVA based on the same profiles of haplogroup distribution frequencies (Table S2). Nevertheless, in terms of the linguistic affinities, the difference between Austronesian (i.e. Cham and WISEA populations) and non-Austronesian (i.e. other MSEA populations) was not statistically significant according to AMOVA (p = 0.08). We incorporated data for eight common Y-STRs (DYS19, DYS389-I, DYS389-II, DYS390, DYS391, DYS392, DYS393, and DYS439) from additional populations in MSEA [11], [12], Multidimensional scaling (MDS) based on RST genetic distances for these Y-STRs did not associate the Chams with populations from WISEA (Figure 4).

[Figure omitted. See PDF.]

Figure 5. Median-joining networks of eight Y-STRs within NRY haplogroups C-M216, F-M213*, K-P131*, O-M7, O-M88, O-M95*, O-M119*, O-M134, O-P191*, O-P200*, and P-P27.1.

Sizes of the circles are proportional to haplotypes frequencies. The lengths of the lines are proportional to the mutational steps.

https://doi.org/10.1371/journal.pone.0036437.g005

[Figure omitted. See PDF.]

Figure 6. Median-joining networks of eight Y-STRs within NRY haplogroups R-M17 and R-M124.

Sizes of the circles are proportional to haplotypes frequencies. The lengths of the lines are proportional to the number of mutational steps.

https://doi.org/10.1371/journal.pone.0036437.g006

[Figure omitted. See PDF.]

Figure 7. NRY and mtDNA haplogroup profiles for the Chams and the Kinhs.

For mtDNA haplogroups, M* includes M17, M20, M21d, M22, M33c, M50, M51, M71, M72, M73, and M77; N* includes N21 and N23; R* includes R22 and R23 [10].

https://doi.org/10.1371/journal.pone.0036437.g007

Admixture in the Cham population

The origin of the Chams could not be simply explained as a demic diffusion of Austronesian immigrants from WISEA. The genetic patterns between the Cham and other Southeast Asian populations, as detected in PCA and MDS, suggested a more complex history. The complex demographic process likely involved genetic admixture with local non-Austronesian speakers in MSEA. Therefore, we performed the admixture analysis [13], [14] to quantify the proportion of genetic contribution from WISEA and MSEA to the Chams (Table 2). The patrilineal contribution from WISEA to the Chams (0.37595) was less than that from MSEA (0.62405). Comparatively, the Vietnamese (most likely, the Kinh) population from southern Vietnam had a dominant proportion of the MSEA contribution (0.842972; Table 2), although the large standard deviation values made the results should be treated with caution.

[Figure omitted. See PDF.]

Table 2. Admixture analysis of the two populations from southern Vietnam.

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

Haplotype diversity analyses

To discern the relationship between the Y-STR haplotypes in the Chams and other Southeast Asians, median-joining networks [15] were constructed using eight common Y-STRs for each of the 11 haplogroups found in the Cham population (Figure 5). In the networks of haplogroups O-M95* and P-P27.1, some haplotypes were exclusively shared between the Cham and ISEA populations. In the networks of C-M216, F-M213*, and K-P131*, some haplotypes in the Chams were derived directly from those in ISEA populations. These lineages in the Chams were most likely introduced by recent gene flow from ISEA. In contrast, the networks for haplogroups O-M7, O-M88, O-M134, O-P191*, and O-P200* indicated closer associations between the Chams and MSEA populations. Most Cham lineages either had identical counterparts or were linked to those haplotypes in MSEA populations; the numbers of mutations between the Chams and MSEA were less than those between the Chams and ISEA (Table S3). These patterns would suggest that these Chams lineages had an in situ origin from MSEA. Among the 48 haplotypes identified in the Chams, 11 and 18 were shared with those in ISEA and MSEA, respectively (Table S3). Nevertheless, the counts for shared haplotype did not differ significantly (two-tailed Fisher's exact test, P = 0.303; Table S4). Moreover, six haplotypes belonging to haplogroup O-M95* were shared by both ISEA and MSEA groups. The exact origin of these lineages in the Chams remains elusive.

To trace the source of the exotic South Asian prevailing components, we incorporated published data (Table S5) from India [16], [17], Pakistan [16], and West Asia [18], [19] and reconstructed median-joining networks of haplogroups R-M17 and R-M124 (Figure 6). All haplotypes in the Chams were scattered in the networks, which implied that these lineages had an origin via recent gene flow rather than deeply rooted ancestry. Two Cham lineages of R-M124 were shared the same haplotype with those from North India. This observation suggested that North India might be the original source of the R-M124 lineages in the Chams. The relationships among lineages of R-M17 were complex in the network, which suggested multiple geographic/ethnic sources for the R-M17 lineages in the Chams.

Discussion

Integrating the information from two uniparentally inherited markers (NRY and mtDNA) is a powerful means of disentangling the human population histories [20], and especially for elucidating sex-biased migrations and social-cultural effects [21]. Compared with our previous study for mtDNA variation in the Chams [10], the current assessment for NRY variation facilitates a better understanding into the origin of the Cham people. Both NRY and mtDNA haplogroup profiles (Figure 7) suggest that, in general, the Chams are indigenous to Southeast Asia. Characteristic East and Southeast Asian lineages, viz., NRY haplogroups O-P191 and C-M217, together with mtDNA haplogroups B, F, M7, and R9, accounted for the majority of the patrilineal (∼67.8%) and matrilineal (∼60.1%) gene pools of the Chams, respectively. Some ancient Southeast Asian components (NRY haplogroups: C-M216*, F-M213*, and K-P131*; mtDNA haplogroups: M*, N*, and R*) were also identified in the Chams.

The origin of the Chams appears to be much more complex, at least based on the results of PCA, MDS, AMOVA, and haplotype (near-) matching analyses. Recent gene flow from ISEA is detected in the patrilineal pool of the Chams, most likely via the dispersal of Austronesian speakers. Further, the Cham population also contains a significant amount of local genetic contributions from non-Austronesian populations in MSEA. This pattern corresponds with our previous study based on mtDNA [10]. Taken together, the origin of the Chams is mainly a result of admixture between the Austronesian immigrants from ISEA with the indigenous populations (most likely, Mon-Khmers) in MSEA.

South Asian NRY haplogroups R-M17, R-M124, and H-M69 [16], [22] are common in the Chams (∼18.6%; Figures 2) yet no mtDNA haplotypes are known [10]. Male South Asians contribute to the genetic makeup of Chams, but not South Asian females. The existence of these South Asian patrilineal lineages was in good accordance with the archaeological and historical records. The dominant religion of the Cham people is known to have been Hinduism (overwhelmingly Shaivism) and their culture was deeply influenced by that of India [3], [7]. Both Indian and Cham people appear to have played important roles in Southeast Asian maritime trade [23], [24]. Contact between the two peoples makes gene flow between them inevitable. The discordance between NRY and mtDNA contributions in the Chams (Figure 7) is well explained by the male-mediated dispersals, most likely through the spread of religions and business trade. In particular, the admixture between alien males and local females is compatible with the matrilocal residence in the Cham people [25], [26].

Patrilineal genetic structuring differs between the Chams and Kinhs. For instance, in contrast to the Chams, frequently the Kinhs have lineages (8/76, ∼10.5%) from the characteristic Chinese haplogroup O-M7 [27] yet only one lineage from the South Asian haplogroup R-M17 (Figure 7). In addition to the Sinicized cultures, substantial Chinese assimilation into the Kinh people via immigration is suggested for northern Vietnam [3], [7]. Thus, the different ethnohistories of the Chams and Kinhs are reflected by their unique mtDNA and NRY patterns.

In summary, this study expands our knowledge on the complex history of the Austronesian diffusion in MSEA. Further improvements to the resolution of the NRY tree [20], [28] will help to unravel the story of the Cham people. This initiative will also benefit from the employment of genome-wide autosomal markers [29], [30], [31]. In the future, a comprehensive study involving extensive sampling will pinpoint more details about the demographic history, such as the source and route for migration, the timing for admixture and expansion.

Materials and Methods

Samples and data collection

Blood samples of 177 unrelated males were collected from four populations (Table 1; Figure 1). Among them, samples from 59 Cham individuals were collected from Binh Thuan province, southern Vietnam. Binh Thuan was part of the Cham principality of Panduranga, the last Cham territory that had been annexed by Nguyen Vietnam in 1832 AD [3], [7], and it was said to harbor a significant number of Chamic speakers [7]. The mtDNA data of the Cham, Kinh, and Thai populations were previously reported [10], [32]. This study was approved by the Institutional Review Board of Kunming Institute of Zoology. All subjects were interviewed to obtain informed written consent before sample collection.

Comparative NRY data from southern China and Southeast Asia (Figure 1; Table 1) were taken from previously published literature [11], [12], [27], [33]. We uniformed all Y-SNPs and Y-STRs data into the same resolution to include as more populations as possible. This truncation of some data caused the NRY haplogroups collapsed into 16 clusters (Table S2), and Y-STRs were reduced to eight loci (DYS19, DYS389-I, DYS389-II, DYS390, DYS391, DYS392, DYS393, and DYS439). Additional data of haplogroups R-M17 and R-124 were collected from published South and West Asian datasets [16], [17], [18], [19] (Table S5).

DNA extraction and genotyping

Genomic DNA was extracted by the standard phenol/chloroform methods. Seventeen Y-SNPs (Table S1) were genotyped by the GenomeLabTM SNPstream® (Beckman Coulter). We used three panels of multiplex PCR reactions following manufacturer's recommendation (Protocol S1). The primers for multiplex PCR and single base extension reactions were designed by Autoprimer software (Beckman Coulter) [34]. To improve the resolution of phylogeny, we further screened nine Y-SNPs by direct sequencing some individuals (Table S1). The PCR amplification and sequencing primers were previously reported [35]. Using described methods [36], [37], [38], we genotyped eight Y-STRs (DYS19, DYS389-I, DYS389-II, DYS390, DYS391, DYS392, DYS393, and DYS439) on an ABI 3730 DNA Analyzer (Applied Biosystems). For DYS389-I and DYS389-II, we used the genotyped data of DYS389-I, and DYS389-II minus DYS389-I in our analyses.

Data analysis

Arlequin 3.5 (http://cmpg.unibe.ch/software/arlequin35/) was used to calculate AMOVA and RST distances [39]. Principal component analysis (PCA) and multidimensional scaling (MDS) were performed using SPSS 13.0 software (SPSS). In PCA, the original haplogroup frequency data were transformed to standardize against the different effect of genetic drift on haplogroups of different frequencies [40]. Admix 2.0 (http://web.unife.it/progetti/genetica/Isabelle/admix2_0.html) was used to estimate the level of admixture of MSEA and WISEA groups in the Cham and Vietnamese populations [13], [14]. The average haplogroup frequencies of MSEA and WISEA were taken for the two parental populations, respectively. Median-joining networks [15] of Y-STRs within certain haplogroups were constructed with NETWORK 4.6 (http://www.fluxus-engineering.com/network_terms.htm).

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Plot of haplogroup contribution of the first and second PC. The contribution of each haplogroup is calculated as the factor scores for PC1 and PC2 with regression (REGR) method in SPSS 13.0 software (SPSS).

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

(TIF)

Table S1.

The variations of Y-SNPs and Y-STRs in 177 unrelated males analyzed in this study.

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

(XLS)

Table S2.

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

Table S3.

https://doi.org/10.1371/journal.pone.0036437.s004

Table S4.

Fisher's exact test for haplotype (near-) matching analysis.

https://doi.org/10.1371/journal.pone.0036437.s005

(DOC)

Table S5.

General information of additional populations analyzed in haplogroups R-M17 and R-M124.

https://doi.org/10.1371/journal.pone.0036437.s006

(XLS)

Protocol S1.

The primers and protocols for the SNPs genotyped by the GenomeLabTM SNPstream®.

https://doi.org/10.1371/journal.pone.0036437.s007

(DOC)

Acknowledgments

We thank Chun-Ling Zhu, Gui-Mei Li, Tai-Cheng Zhou, and Shuai Hua for technical assistance. We are grateful to all volunteers for providing DNA samples. We thank Drs. Jatupol Kampuansai and Tatiana Karafet for supplying their NRY data. We also thank the two anonymous referees for their valuable comments.

Author Contributions

Conceived and designed the experiments: Y-PZ J-DH M-SP. Performed the experiments: J-DH J-QJ. Analyzed the data: M-SP J-DH. Contributed reagents/materials/analysis tools: M-SP J-DH HHQ KPD AVT S-FW J-QJ. Wrote the paper: Y-PZ J-DH M-SP Y-GY RWM.

Citation: He J-D, Peng M-S, Quang HH, Dang KP, Trieu AV, Wu S-F, et al. (2012) Patrilineal Perspective on the Austronesian Diffusion in Mainland Southeast Asia. PLoS ONE7(5): e36437. https://doi.org/10.1371/journal.pone.0036437

References

1. Bellwood P, Fox JJ, Tryon D (2006) The Austronesians in history: common origins and diverse transformations. In: Bellwood P, Fox JJ, Tryon D, editors. The Austronesians: historical and comparative perspectives. Canberra: ANU E Press. pp. 1–14.P. BellwoodJJ FoxD. Tryon2006The Austronesians in history: common origins and diverse transformations.P. BellwoodJJ FoxD. TryonThe Austronesians: historical and comparative perspectivesCanberraANU E Press114

2. Lewis MP, editor. (2009) MP Lewis2009Ethnologue: languages of the world. Sixteenth edition. Dallas, Tex: SIL International. Online version: http://www.ethnologue.com/. Ethnologue: languages of the world. Sixteenth edition. Dallas, Tex: SIL International. Online version: http://www.ethnologue.com/.

3. He P (2006) The origin and evolution of nationalities in the Indochina Peninsula. Beijing: The Ethnic Publishing House. P. He2006The origin and evolution of nationalities in the Indochina Peninsula.BeijingThe Ethnic Publishing House

4. Southworth WA (2004) The coastal states of Champa. In: Glover I, Bellwood PS, editors. Southeast Asia: from prehistory to history. London: RoutledgeCurzon. pp. 209–233.WA Southworth2004The coastal states of Champa.I. GloverPS BellwoodSoutheast Asia: from prehistory to historyLondonRoutledgeCurzon209233

5. Thurgood G (2005) A preliminary sketch of Phan Rang Cham. In: Adelaar KA, Himmelmann N, editors. The Austronesian languages of Asia and Madagascar. London: Routledge. pp. 489–512.G. Thurgood2005A preliminary sketch of Phan Rang Cham.KA AdelaarN. HimmelmannThe Austronesian languages of Asia and MadagascarLondonRoutledge489512

6. Thurgood G (1999) From ancient Cham to modern dialects: two thousand years of language contact and change. Honolulu, Hawaii: University of Hawai'i Press. xvii: 407.G. Thurgood1999From ancient Cham to modern dialects: two thousand years of language contact and change. Honolulu, Hawaii: University of Hawai'i Press.xvii407

7. Higham C (2002) Early cultures of Mainland Southeast Asia. Bangkok: River Books. C. Higham2002Early cultures of Mainland Southeast Asia.BangkokRiver Books

8. Bellwood P (2006) Austronesian prehistory in Southeast Asia: homeland, expansion and transformation. In: Bellwood P, Fox JJ, Tryon D, editors. The Austronesians: historical and comparative perspectives. Canberra: ANU E Press. pp. 103–114.P. Bellwood2006Austronesian prehistory in Southeast Asia: homeland, expansion and transformation.P. BellwoodJJ FoxD. TryonThe Austronesians: historical and comparative perspectivesCanberraANU E Press103114

9. Solheim WG, Bulbeck D, Flavel A (2006) Archaeology and culture in Southeast Asia: unraveling the Nusantao. Quezon City: University of the Philippines Press. WG SolheimD. BulbeckA. Flavel2006Archaeology and culture in Southeast Asia: unraveling the Nusantao.Quezon CityUniversity of the Philippines Press

10. Peng MS, Quang HH, Dang KP, Trieu AV, Wang HW, et al. (2010) Tracing the Austronesian footprint in Mainland Southeast Asia: a perspective from mitochondrial DNA. Mol Biol Evol 27: 2417–2430.MS PengHH QuangKP DangAV TrieuHW Wang2010Tracing the Austronesian footprint in Mainland Southeast Asia: a perspective from mitochondrial DNA.Mol Biol Evol2724172430

11. Kim SH, Kim KC, Shin DJ, Jin HJ, Kwak KD, et al. (2011) High frequencies of Y-chromosome haplogroup O2b-SRY465 lineages in Korea: a genetic perspective on the peopling of Korea. Investig Genet 2: 10.SH KimKC KimDJ ShinHJ JinKD Kwak2011High frequencies of Y-chromosome haplogroup O2b-SRY465 lineages in Korea: a genetic perspective on the peopling of Korea.Investig Genet210

12. Kutanan W, Kampuansai J, Fuselli S, Nakbunlung S, Seielstad M, et al. (2011) Genetic structure of the Mon-Khmer speaking groups and their affinity to the neighbouring Tai populations in Northern Thailand. BMC Genet 12: 56.W. KutananJ. KampuansaiS. FuselliS. NakbunlungM. Seielstad2011Genetic structure of the Mon-Khmer speaking groups and their affinity to the neighbouring Tai populations in Northern Thailand.BMC Genet1256

13. Bertorelle G, Excoffier L (1998) Inferring admixture proportions from molecular data. Mol Biol Evol 15: 1298–1311.G. BertorelleL. Excoffier1998Inferring admixture proportions from molecular data.Mol Biol Evol1512981311

14. Dupanloup I, Bertorelle G (2001) Inferring admixture proportions from molecular data: extension to any number of parental populations. Mol Biol Evol 18: 672–675.I. DupanloupG. Bertorelle2001Inferring admixture proportions from molecular data: extension to any number of parental populations.Mol Biol Evol18672675

15. Bandelt HJ, Forster P, Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48.HJ BandeltP. ForsterA. Röhl1999Median-joining networks for inferring intraspecific phylogenies.Mol Biol Evol163748

16. Sengupta S, Zhivotovsky LA, King R, Mehdi SQ, Edmonds CA, et al. (2006) Polarity and temporality of high-resolution Y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of central Asian pastoralists. Am J Hum Genet 78: 202–221.S. SenguptaLA ZhivotovskyR. KingSQ MehdiCA Edmonds2006Polarity and temporality of high-resolution Y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of central Asian pastoralists.Am J Hum Genet78202221

17. Thangaraj K, Naidu BP, Crivellaro F, Tamang R, Upadhyay S, et al. (2010) The influence of natural barriers in shaping the genetic structure of Maharashtra populations. PLoS ONE 5: e15283.K. ThangarajBP NaiduF. CrivellaroR. TamangS. Upadhyay2010The influence of natural barriers in shaping the genetic structure of Maharashtra populations.PLoS ONE5e15283

18. Mohammad T, Xue Y, Evison M, Tyler-Smith C (2009) Genetic structure of nomadic Bedouin from Kuwait. Heredity 103: 425–433.T. MohammadY. XueM. EvisonC. Tyler-Smith2009Genetic structure of nomadic Bedouin from Kuwait.Heredity103425433

19. Cadenas AM, Zhivotovsky LA, Cavalli-Sforza LL, Underhill PA, Herrera RJ (2008) Y-chromosome diversity characterizes the Gulf of Oman. Eur J Hum Genet 16: 374–386.AM CadenasLA ZhivotovskyLL Cavalli-SforzaPA UnderhillRJ Herrera2008Y-chromosome diversity characterizes the Gulf of Oman.Eur J Hum Genet16374386

20. Underhill PA, Kivisild T (2007) Use of Y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet 41: 539–564.PA UnderhillT. Kivisild2007Use of Y chromosome and mitochondrial DNA population structure in tracing human migrations.Annu Rev Genet41539564

21. Wilkins JF, Marlowe FW (2006) Sex-biased migration in humans: what should we expect from genetic data? Bioessays 28: 290–300.JF WilkinsFW Marlowe2006Sex-biased migration in humans: what should we expect from genetic data?Bioessays28290300

22. Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, et al. (2003) The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations. Am J Hum Genet 72: 313–332.T. KivisildS. RootsiM. MetspaluS. MastanaK. Kaldma2003The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations.Am J Hum Genet72313332

23. Reid A (1988) Southeast Asia in the age of commerce, 1450–1680 Volume One: the lands below the winds. New Heaven: Yale University Press. A. Reid1988Southeast Asia in the age of commerce, 1450–1680 Volume One: the lands below the winds.New HeavenYale University Press

24. Reid A (1988) Southeast Asia in the age of commerce, 1450–1680 Volume Two: expansion and crisis. New Heaven: Yale University Press. A. Reid1988Southeast Asia in the age of commerce, 1450–1680 Volume Two: expansion and crisis.New HeavenYale University Press

25. Jordan FM, Gray RD, Greenhill SJ, Mace R (2009) Matrilocal residence is ancestral in Austronesian societies. Proc R Soc B 276: 1957–1964.FM JordanRD GraySJ GreenhillR. Mace2009Matrilocal residence is ancestral in Austronesian societies.Proc R Soc B27619571964

26. Nakamura R (2009) Awar-Ahier: Two Keys to Understanding the Cosmology and Ethnicity of the Cham People (Ninh Thuân Province, Vietnam). In: Hardy A, Cucarzi M, Zolese P, editors. Champa and the archaeology of Mỹ So’n (Vietnam). Singapore: NUS Press. pp. 78–106.R. Nakamura2009Awar-Ahier: Two Keys to Understanding the Cosmology and Ethnicity of the Cham People (Ninh Thuân Province, Vietnam).A. HardyM. CucarziP. ZoleseChampa and the archaeology of Mỹ So’n (Vietnam)SingaporeNUS Press78106

27. Karafet TM, Hallmark B, Cox MP, Sudoyo H, Downey S, et al. (2010) Major east-west division underlies Y chromosome stratification across Indonesia. Mol Biol Evol 27: 1833–1844.TM KarafetB. HallmarkMP CoxH. SudoyoS. Downey2010Major east-west division underlies Y chromosome stratification across Indonesia.Mol Biol Evol2718331844

28. Cruciani F, Trombetta B, Massaia A, Destro-Bisol G, Sellitto D, et al. (2011) A revised root for the human y chromosomal phylogenetic tree: the origin of patrilineal diversity in Africa. Am J Hum Genet 88: 814–818.F. CrucianiB. TrombettaA. MassaiaG. Destro-BisolD. Sellitto2011A revised root for the human y chromosomal phylogenetic tree: the origin of patrilineal diversity in Africa.Am J Hum Genet88814818

29. Friedlaender JS, Friedlaender FR, Reed FA, Kidd KK, Kidd JR, et al. (2008) The genetic structure of Pacific Islanders. PLoS Genet 4: e19.JS FriedlaenderFR FriedlaenderFA ReedKK KiddJR Kidd2008The genetic structure of Pacific Islanders.PLoS Genet4e19

30. Kayser M, Lao O, Saar K, Brauer S, Wang XY, et al. (2008) Genorne-wide analysis indicates more Asian than melanesian ancestry of polynesians. Am J Hum Genet 82: 194–198.M. KayserO. LaoK. SaarS. BrauerXY Wang2008Genorne-wide analysis indicates more Asian than melanesian ancestry of polynesians.Am J Hum Genet82194198

31. The HUGO Pan-Asian Consortium (2009) Mapping human genetic diversity in Asia. Science 326: 1541–1545.Consortium The HUGO Pan-Asian2009Mapping human genetic diversity in Asia.Science32615411545

32. Yao YG, Nie L, Harpending H, Fu YX, Yuan ZG, et al. (2002) Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity. Am J Phys Anthropol 118: 63–76.YG YaoL. NieH. HarpendingYX FuZG Yuan2002Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity.Am J Phys Anthropol1186376

33. Cai X, Qin Z, Wen B, Xu S, Wang Y, et al. (2011) Human migration through bottlenecks from Southeast Asia into East Asia during Last Glacial Maximum revealed by Y chromosomes. PLoS ONE 6: e24282.X. CaiZ. QinB. WenS. XuY. Wang2011Human migration through bottlenecks from Southeast Asia into East Asia during Last Glacial Maximum revealed by Y chromosomes.PLoS ONE6e24282

34. Debnath M, Palanichamy MG, Mitra B, Jin JQ, Chaudhuri TK, et al. (2011) Y-chromosome haplogroup diversity in the sub-Himalayan Terai and Duars populations of East India. J Hum Genet 56: 765–771.M. DebnathMG PalanichamyB. MitraJQ JinTK Chaudhuri2011Y-chromosome haplogroup diversity in the sub-Himalayan Terai and Duars populations of East India.J Hum Genet56765771

35. Karafet TM, Mendez FL, Meilerman MB, Underhill PA, Zegura SL, et al. (2008) New binary polymorphisms reshape and increase resolution of the human Y chromosomal haplogroup tree. Genome Res 18: 830–838.TM KarafetFL MendezMB MeilermanPA UnderhillSL Zegura2008New binary polymorphisms reshape and increase resolution of the human Y chromosomal haplogroup tree.Genome Res18830838

36. Thomas MG, Bradman N, Flinn HM (1999) High throughput analysis of 10 microsatellite and 11 diallelic polymorphisms on the human Y-chromosome. Hum Genet 105: 577–581.MG ThomasN. BradmanHM Flinn1999High throughput analysis of 10 microsatellite and 11 diallelic polymorphisms on the human Y-chromosome.Hum Genet105577581

37. Ayub Q, Mohyuddin A, Qamar R, Mazhar K, Zerjal T, et al. (2000) Identification and characterisation of novel human Y-chromosomal microsatellites from sequence database information. Nucleic Acids Res 28: e8.Q. AyubA. MohyuddinR. QamarK. MazharT. Zerjal2000Identification and characterisation of novel human Y-chromosomal microsatellites from sequence database information.Nucleic Acids Res28e8

38. Prinz M, Boll K, Baum H, Shaler B (1997) Multiplexing of Y chromosome specific STRs and performance for mixed samples. Forensic Sci Int 85: 209–218.M. PrinzK. BollH. BaumB. Shaler1997Multiplexing of Y chromosome specific STRs and performance for mixed samples.Forensic Sci Int85209218

39. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567.L. ExcoffierHEL Lischer2010Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows.Mol Ecol Resour10564567

40. Richards M, Macaulay V, Torroni A, Bandelt HJ (2002) In search of geographical patterns in European mitochondrial DNA. Am J Hum Genet 71: 1168–1174.M. RichardsV. MacaulayA. TorroniHJ Bandelt2002In search of geographical patterns in European mitochondrial DNA.Am J Hum Genet7111681174