Genetic continuity across a deeply divergent linguistic contact zone in North Maluku, Indonesia
© Wilder et al; licensee BioMed Central Ltd. 2011
Received: 4 July 2011
Accepted: 18 November 2011
Published: 18 November 2011
The islands of North Maluku, Indonesia occupy a central position in the major prehistoric dispersal streams that shaped the peoples of Island Southeast Asia and the Pacific. Within this region a linguistic contact zone exists where speakers of Papuan and Austronesian languages reside in close proximity. Here we use population genetic data to assess the extent to which North Maluku populations experienced admixture of Asian genetic material, and whether linguistic boundaries reflect genetic differentiation today.
Autosomal and X-linked markers reveal overall Asian admixture of 67% in North Maluku, demonstrating a substantial contribution of genetic material into the region from Asia. We observe no evidence of population structure associated with ethnicity or language affiliation.
Our data support a model of widespread Asian admixture in North Maluku, likely mediated by the expansion of Austronesian-speaking peoples into the region during the mid Holocene. In North Maluku there is no genetic differentiation in terms of Austronesian- versus Papuan-speakers, suggesting extensive gene flow across linguistic boundaries. In a regional context, our results illuminate a major genetic divide at the Molucca Sea, between the islands of Sulawesi and North Maluku. West of this divide, populations exhibit predominantly Asian ancestry, with very little contribution of Papuan genetic material. East of the Molucca Sea, populations show diminished rates of Asian admixture and substantial persistence of Papuan genetic diversity.
The islands of North Maluku, Indonesia (also known as the Moluccas or Spice Islands) lie in a pivotal position on the migration routes of humans through the Indo-Pacific. Two major prehistoric movements of people are particularly important with respect to the linguistic, cultural, and genetic diversity of the region. First, the initial Pleistocene colonization of the western Pacific flowed from west to east through Indonesia to New Guinea and Australia, bringing the ancestors of modern Papuans and Australians into the region. Although the exact route(s) of this migration are not known, North Maluku lies along a probable pathway of inter-island travel between Sulawesi and New Guinea . Indeed, archaeological evidence confirms human habitation in North Maluku as far back as 32,500 years before present . The second major migration event with a significant impact on North Maluku was the spread of Asian Austronesian-speaking populations from the north via the Philippines into island Southeast Asia and Oceania. The oldest pottery found in North Maluku and associated with the Austronesian culture dates to 3,500 years before present . Several linguistic and cultural aspects of the Austronesian expansion appear to have originated in the islands of North Maluku, suggesting it played an important role in the spread of Austronesians throughout the region . As such, the prehistoric contributions to present-day communities in this region can largely be framed around these two main events.
The mechanism by which Austronesian languages became established in North Maluku remains unknown. The diversity of Papuan languages in the area suggests that they were well-established prior to the arrival of Austronesians [2, 7]. Two simple models of the establishment of Austronesian languages that make opposing genetic predictions are a Replacement Model and an Adoption Model. The first of these posits that Austronesian-speakers partially colonized North Maluku, displacing speakers of Papuan languages in the process. Even if there was some intermarriage between Papuan- and Austronesian-speakers, this model suggests that populations representing these two language groups will be genetically distinct from one another . Conversely, the Adoption Model posits that speakers of Papuan languages in North Maluku adopted Austronesian languages through cultural contact with Austronesian speakers, not intermarriage. This model predicts no major change in the genetic composition of the resident pre-existing population [8, 9], thereby proposing that there will be a minimal signature of Austronesian admixture in North Maluku and little genetic differentiation between Papuan and Austronesian-speaking populations. A genetic pattern consistent with the Adoption Model has been observed in many Austronesian-speaking populations in eastern Indonesia and Melanesia, and has been interpreted as evidence of widespread language borrowing in this region . Of course, the Replacement and Adoption Models reflect opposite extremes on a spectrum of more complex scenarios.
Given the remarkable linguistic diversity of North Maluku, and the obvious importance of its geographic position for both Pleistocene and Austronesian migration, we have surprisingly little knowledge of the genetic composition of populations living in this region today. Because languages in North Maluku have both Asian and Papuan origins, our goal here is to test whether there is a genetic signature that also reflects this division.
One approach to this question would be to genotype large numbers (i.e., hundreds of thousands) of single nucleotide polymorphisms (SNPs) across a range of samples (e.g., using a SNP chip). However, from the perspective of information content, this approach is expensive. The vast majority of SNPs in a population are unlikely to show frequency differences between Asian and Papuan groups, and therefore have no statistical power to infer Asian-Papuan ancestry. Instead, information about ancestry proportions is carried by a relatively small number of SNPs with high frequency differences between the parental groups. Such SNPs are called ancestry informative markers (AIMs). AIMs are considerably cheaper to screen than genotyping with standard SNP chips, but they retain extremely high power to infer ancestry proportions. One key limitation is that AIMs only have statistical power to answer the questions they were designed to address - in this case, proportions of Asian versus Papuan ancestry. Nevertheless, this is the main axis of population history in the Indo-Pacific region today.
The research presented here genotypes 27 previously identified AIMs in 340 individuals representing 11 different ethnicities in North Maluku (Figure 1). This sampling scheme includes multiple populations speaking both Papuan and Austronesian languages. The study has two main aims: 1) to estimate the amount of Asian admixture among populations in the region; and 2) to test whether linguistically differentiated groups are also genetically distinct from one another. Addressing these questions will enhance our understanding of how the remarkable contact zone in North Maluku originally formed, as well as shed light on the dynamics of the expansion of Austronesian speakers through the region.
Asian Admixture in North Maluku
To detect the genetic signature of Asian admixture into North Maluku, we genotyped 11 autosomal and 16 X-chromosomal AIMs in a sample of approximately 340 individuals. These unlinked and statistically independent markers were chosen because they have alleles that are highly differentiated between Asian and Melanesian source populations and are therefore useful to detect contributions from each of these potential parental sources . The source populations used to identify these AIMs were the Southern Han Chinese and highland Papua New Guineans. Other potential source populations have similar allele frequencies at these AIMs, and the choice of parental population has little effect on the power to detect admixture along this Asian-Papuan axis . Because we typed numerous markers from both the X chromosome and the autosomes, we are also able to explore potential sex biases in rates of admixture.
Asian admixture estimates
All Loci Combined
0.61 - 0.68
0.66 - 0.74
0.65 - 0.70
0.58 - 0.73
0.69 - 0.84
0.66 - 0.77
0.38 - 0.67
0.66 - 1.00
0.58 - 1.00
0.61 - 0.69
0.65 - 0.75
0.64 - 0.71
0.58 - 0.70
0.71 - 0.85
0.67 - 0.77
0.63 - 0.71
0.62 - 0.72
0.64 - 0.70
0.60 - 0.78
0.65 - 0.87
0.65 - 0.81
0.58 - 0.65
0.66 - 0.75
0.64 - 0.70
0.63 - 0.80
0.58 - 0.76
0.63 - 0.76
0.57 - 0.72
0.62 - 0.80
0.62 - 0.74
0.50 - 0.80
0.62 - 0.97
0.61 - 0.89
0.50 - 0.61
0.64 - 0.76
0.60 - 0.68
0.48 - 0.80
0.40 - 0.84
0.48 - 0.79
Genetic Differentiation of North Maluku Populations
Hierarchical analysis of genetic differentiation (AMOVA)
Among Language Groups
Among Populations, Within Groups
Generally speaking, there was very little variability among markers in terms of the AMOVA results (Additional File 2). With one exception, the among-group and among-population components of variation for each AIM did not significantly exceed values expected by chance under panmixia. The sole exception was for X-linked marker rs5987967, which had an among-populations component of 5.19% (p < 0.0001); however, the among-group component for this marker was not significantly higher than expected by chance (and was, in fact, statistically indistinguishable from zero). These results suggest that for the X-linked and autosomal AIMs examined here, there is no meaningful among-group or among-population component of sampled genetic variation (i.e., little genetic variation is apportioned according to linguistic groups or ethnicities).
Populations in North Maluku have been substantially impacted by gene flow from Asia. At the ancestry informative loci tested here, we estimate the overall Asian admixture fraction to be ~67% (using two different methodologies). This fraction appears slightly higher on the X-chromosome, indicating a possible Asian female bias during the admixture process. In a regional context, these results mirror previously surveyed populations in eastern Indonesia. Several islands have broadly similar admixture profiles (e.g., Sumba, Alor and Flores), which also show a trend toward female-biased Asian admixture using the same AIMs we screened here . It is interesting to note that Asian admixture estimates for Sulawesi, which is the major island immediately to the west of North Maluku, are much higher than observed in the Maluku islands (97% among three Sulawesi populations). Indeed, from Sulawesi west, all Island Southeast Asian populations exhibit near-complete Asian ancestry at these genetic markers, while populations surveyed to the east of North Maluku (including coastal New Guinea) show a substantially lower Asian contribution. As such, the Molucca Sea, between Sulawesi and North Maluku, represents a geographic breakpoint, to the east of which chromosomes of Papuan ancestry have been retained at appreciable frequency. Interestingly, this phylogeographic pattern is mirrored by the distribution of Austronesian languages. West of the Molucca Sea, indigenous languages have been completely replaced by Austronesian languages; to the east of the Molucca Sea, indigenous languages have tended more often to persist [5, 13]. Analyses of multiple genetic marker systems suggest a similar genetic breakpoint in the Nusa Tenggara region of Indonesia [11, 14], an area that also marks a linguistic boundary between Papuan and Austronesian languages . These patterns suggest that there is a major east-west divide in Asian influence, with respect to both genes and language, in Island Southeast Asia that runs from the Molucca Sea south to Nusa Tenggara. One possible explanation for this pattern is that speakers of Austronesian languages, carrying the Asian alleles at the AIMs we tested here, became less likely to replace existing populations (genetically or linguistically) as they moved east through Island Southeast Asia. Another alternative is that Papuan back-migration occurred subsequent to the Austronesian expansion, reaching its limits at the observed breakpoints. Both of these processes appear to be ongoing today .
It is remarkable that the genetic divide we observe at the Molucca Sea (and which has been observed previously in Nusa Tenggara) precisely follows a breakpoint in human morphological phenotypes first observed by the nineteenth century biologist Alfred R. Wallace . Wallace described this phenotypic boundary as the line dividing people with Malay versus Papuan physical characteristics. This breakpoint is displaced to the east of his more famous biogeographical divide (i.e., Wallace's Line), which he described based on the region's non-human fauna. Our data suggest that differing rates of Asian admixture, likely mediated by the geographical expansion of Austronesian-speakers, may be responsible for Wallace's phenotypic boundary observed for human populations in the region. Within North Maluku, Wallace also noted that individuals displayed what he deemed to be a composite of Papuan and Malay characteristics . Similar observations have been reiterated by modern anthropologists [17–19]. The results we report here suggest that these observations may reflect genuine genetic admixture between Asian and Papuan populations. This diverse ancestry of individuals in North Maluku suggests that they may harbour more deeply divergent genetic lineages than might be expected for small island populations, which likely has implications for health and medical care in this area. In particular, a number of alleles that confer protection against severe malaria, and also sometimes have negative pleiotropic effects, have been associated with either Austronesian or Melanesian origins [20–25]. Given the lack of population structure that we observe, these functional genetic variants of diverse evolutionary origins may now be geographically widespread across North Maluku.
Despite the relatively strong contribution of Asian genetic material to the peoples of North Maluku, we observe no correlation between genes and language among our surveyed populations. Whether Papuan- or Austronesian-speaking, all populations shared similar levels of Asian admixture and showed no evidence of genetic differentiation. This result is not easy to reconcile with either the Replacement or Adoption Models. Unless swamped by subsequent gene flow, the Replacement hypothesis predicts some degree of genetic differentiation between populations speaking Austronesian versus Papuan languages. Likewise, the Adoption hypothesis predicts little genetic evidence of Austronesian admixture. Neither of these patterns is consistent with our observations. Instead, our results suggest a process of extensive Asian admixture with heterogeneous linguistic replacement. Austronesian languages made inroads in North Maluku, but this process only sporadically resulted in replacement of indigenous languages (predominantly on the southern and eastern coasts of Halmahera). Why some Papuan languages persisted in the face of widespread gene flow, and others did not, represents an intriguing open question.
This outcome contrasts strongly with other nearby island systems, including many populations in coastal Melanesia, where replacement of indigenous languages with Austronesian ones is complete, despite the persistence of Papuan genetic variation [10, 26]. Indeed, within Island Southeast Asia, previous studies have found a heterogeneous relationship between genes and languages in different settings. For instance, Lansing et al.  found that on the island of Sumba there is a positive correlation between villages that have retained linguistic elements that likely trace to the single founding Austronesian language of the island and the frequency of Austronesian Y chromosome lineages. In contrast, Mona et al.  found no relationship between speakers of Austronesian versus Papuan languages with respect to Y chromosome or mtDNA variation in a multi-island comparison in the Indonesian province of East Nusa Tenggara. The island of New Guinea has produced especially mixed results suggestive of substantial local-scale heterogeneity in gene-language correlations; neither mtDNA or Y chromosome show associations with language in Indonesian West New Guinea, while studies of classical markers have sometimes shown a very strong association and sometimes not [29–32]. Regardless of the mechanism by which Austronesian languages replaced Papuan ones in North Maluku, our study suggests that deeply divergent languages do not necessarily represent a barrier to gene flow in this region, a pattern that is mirrored in many other parts of eastern Indonesia and the wider Melanesian area .
One important caveat of our AIM-based analysis is that it is designed expressly to tease out patterns of population genetic variation that have been shaped by Asian introgression into predominantly Papuan regions. Our finding of a lack of differentiation could be caused by either high rates of gene flow among populations or uniformly high rates of Asian admixture into individual populations, but not necessarily both. In other words, it would be premature to infer from our results that there exists no genetic structure that differentiates populations within North Maluku. However, if such population structure exists, it has not been caused by differential rates of Asian admixture. Other genetic systems may reveal important patterns of subdivision that result from alternative historical processes . For instance, North Maluku was the major commercial hub for the long-distance spice trade to Asia and Europe during the last two millennia . The genetic legacy of this global interaction may be better addressed using a different suite of genetic markers.
Populations in North Maluku inhabit a remarkable linguistic contact zone reflecting deeply divergent Papuan and Austronesian languages. Our data suggest that the arrival of Austronesian languages was accompanied by extensive genetic admixture across North Maluku together with sporadic replacement of indigenous Papuan languages by Austronesian ones. In the broader context of the Austronesian expansion into this region, we have identified the Molucca Sea as an important phylogeographic and linguistic divide, east of which indigenous languages and genes resisted replacement more effectively than to the west. The cause of this pattern remains unclear, but likely reflects a geographical or social transition point in the expansion dynamics of the Austronesian dispersal .
DNA samples were collected from unrelated adults of both sexes attending secondary school in the city of Ternate, North Maluku. Data on natal language and self-identified ethnicity were obtained from each individual. Because DNA collection occurred at a centralized location, the number of samples representing each population varied widely. For population-based analyses, we included in this study only those for which the sample size was at least five individuals. All samples were collected with informed consent by staff of the Eijkman Institute for Molecular Biology, using protocols approved by the Institute.
A total of 28 ancestry informative markers (AIMs) (see Additional File 1) were selected from previous work . Briefly, these AIMs were identified using a SNP screening process designed to find those that showed high FST between southern Han Chinese and PNG highlanders. Candidate markers were chosen from two sources: the HOMINID dataset, a collection of re-sequenced putatively neutral regions distributed across the human genome , and the Jakobsson et al. dataset, a collection of 500,000 SNPs typed in the HGDP-CEPH panel . To address the question of sex-specific admixture, markers were selected from both the autosomes and the X chromosome. To minimize the effects of natural selection, all AIMs were located away from genes (including introns, UTRs and immediate flanking regions) and are more than 1 cM distant from other markers in the panel (i.e., they are evolutionarily and statistically independent).
AIM genotype data were collected using a MassARRAY® iPLEX Gold SNP-typing platform (Sequenom) at the University of Arizona Genetics Core facility (Additional File 3). Genotype data from the 28 AIMs were initially checked for agreement with Hardy-Weinberg proportions and evidence of minimal sample dropout. One autosomal marker, rs12613102, showed significant deviations from expected Hard-Weinberg proportions and a high number of samples produced no genotype calls. The affected samples were not associated with missing data at other markers, and this marker was eliminated from all further analyses (data are not shown in Additional File 1).
Admixture estimates and confidence intervals were generated for each population using a variant of the admixture estimation method introduced by Chakraborty et al. , as modified by Cox et al. . Briefly, we inferred admixture rates using a weighted least-squares estimator  altered to account for the sampling error in each of the 'parental' (P1 and P2) and 'hybrid' (H) populations by inferring a frequency density for P1, P2 and H at each marker. Using a re-sampling approach, random variables were drawn from the three frequency distributions, admixture was calculated using the least-squares estimator and the process was repeated 105 times. The median admixture rate with 95% confidence intervals (i.e., 0.025 and 0.975 quantiles) was calculated from the distribution of re-sampled admixture rates. Note that this approach explicitly accounts for variation in sample sizes and effective population sizes, even when they are quite small and/or differ widely between populations.
Population differentiation was assessed using several methods. First, we employed a hierarchical analysis of variance (AMOVA), as described by Weir . For this analysis, we considered populations, which were determined by self-identified ethnicity (described above), as a subordinate classification to language group (Austronesian- versus Papuan-speakers). Language group affiliations were determined using data compiled in . Significance of individual components of variance (i.e., FST, FCT and FSC) was assessed using 1,000 random permutations of the data. We implemented the AMOVA using the program Arlequin version 188.8.131.52 . Our second examination of differentiation among populations used the Bayesian clustering algorithm STRUCTURE (v. 2.3.3) . We implemented STRUCTURE using a dataset containing AIM genotypes from our test samples from North Maluku together with the "parental" populations from which our ancestry informative markers were initially chosen (southern Han Chinese and highland Papua New Guinea). Data for the parental populations was taken from . We ran STRUCTURE with a burn-in period of at least 45,000 steps and 5 × 106 steps per run, varying K (the number of populations) from 1 to 5. We ran each condition multiple times to ensure similarity of results among runs. To improve the detection of structure in our data we allowed prior information about individuals to inform the clustering process, as described by Hubisz et al. . All results shown here use language affiliation (i.e., Papuan versus Austronesian) as prior information. Use of ethnicity produced very similar results (data not shown).
As a final examination of population structure, we performed an individual-based principal components analysis (PCA), implemented using the program SMARTPCA , to test for differentiation between Austronesian- and Papuan-speakers. In addition to providing a graphical depiction of the position of each individual along the principal eigenvectors, this method applies a formal statistical test based on Tracy-Widom theory to detect structure in the dataset.
Chase Ridenour and John Zadrozny assisted with data collection at Northern Arizona University. Ryan Sprissler and Elizabeth Zukowski assisted with marker genotyping at the University of Arizona. Comments from three anonymous reviewers helped to improve an earlier draft of this work. Funding that enabled this study was provided to JAW from the Department of Biological Sciences at Northern Arizona University and the National Science Foundation (BCS-1062258). MPC was supported by a Rutherford Fellowship from the Royal Society of New Zealand (RDF-10-MAU-001).
- Birdsell JB: The recalibration of a paradigm for the first peopling of Greater Australia. Sunda and Sahul : prehistoric studies in Southeast Asia, Melanesia and Australia. 1977, London, New York: Academic Press, 113-167.Google Scholar
- Bellwood P, Nitihaminoto G, Irwin G, Gunadi , Waluyo A, Tanudirjo D: 35,000 years of prehistory in the northern Moluccas. Bird's Head approaches : Irian Jaya studies, a programme for interdisciplinary research. Edited by: Bartstra G-J. Rotterdam: Balkema. 1998, 233-275.Google Scholar
- Spriggs M: Research questions in Maluku Archaeology. Cakelele. 1998, 9 (2): 51-64.Google Scholar
- Voorhoeve C: The languages of the North Halmaheran stock. Pacific Linguistics. 1988, A-76 (26): 181-209.Google Scholar
- Lewis MP, (ed.): Ethnologue: Languages of the World. 2009, Dallas: SIL International, 16Google Scholar
- Blust RA: Eastern Malayo-Polynesian: A subgrouping argument. Second International Conference of Austronesian Linguistics. Edited by: Wurm SA, Carrington L. 1798, Canberra: Australia National University, 181-234.Google Scholar
- Voorhoeve C: Contact-induced change in the non-Austronesian languages in the north Moluccas, Indonesia. Language Contact and Change in the Austronesian World. Edited by: Dutton T, Tryon DT. 1994, Berlin, New York: Mouton de Gruyter, 649-674.Google Scholar
- Diamond J, Bellwood P: Farmers and their languages: the first expansions. Science. 2003, 300 (5619): 597-603. 10.1126/science.1078208.View ArticlePubMedGoogle Scholar
- Cavalli-Sforza LL, Piazza A, Menozzi P, Mountain J: Reconstruction of human evolution: bringing together genetic, archaeological, and linguistic data. Proc Natl Acad Sci USA. 1988, 85 (16): 6002-6006. 10.1073/pnas.85.16.6002.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedlaender JS, Friedlaender FR, Reed FA, Kidd KK, Kidd JR, Chambers GK, Lea RA, Loo JH, Koki G, Hodgson JA: The genetic structure of Pacific Islanders. PLoS genetics. 2008, 4 (1): e19-10.1371/journal.pgen.0040019.PubMed CentralView ArticlePubMedGoogle Scholar
- Cox MP, Karafet TM, Lansing JS, Sudoyo H, Hammer MF: Autosomal and X-linked single nucleotide polymorphisms reveal a steep Asian-Melanesian ancestry cline in eastern Indonesia and a sex bias in admixture rates. Proc Biol Sci. 2010, 277 (1687): 1589-1596. 10.1098/rspb.2009.2041.PubMed CentralView ArticlePubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14 (8): 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Bellwood P: Comment on "Farming and Language in Island Southeast Asia: Reframing Austronesian History". Curr Anthropol. 2010, 51 (2): 240-241.Google Scholar
- Karafet TM, Hallmark B, Cox MP, Sudoyo H, Downey S, Lansing JS, Hammer MF: Major east-west division underlies Y chromosome stratification across Indonesia. Mol Biol Evol. 2010, 27 (8): 1833-1844. 10.1093/molbev/msq063.View ArticlePubMedGoogle Scholar
- Lansing JS, Cox MP, de Vet TA, Downey SS, Hallmark B, Sudoyo H: An ongoing Austronesian expansion in Island Southeast Asia. J Anthropol Archaeol. 2011, 30 (3): 262-272. 10.1016/j.jaa.2011.06.004.View ArticleGoogle Scholar
- Wallace AR: The Malay Archipelago. 2004, London; New York: RoutledgeGoogle Scholar
- Ishige N: The Galela of Halmahera : a preliminary survey. 1980, Suita, Osaka, Japan: National Museum of EthnologyGoogle Scholar
- Bellwood PS: Prehistory of the Indo-Malaysian Archipelago. 2007, Canberra, Australia: ANU E Press, 3Google Scholar
- LeBar FM, Appell GN: Ethnic groups of insular Southeast Asia. 1972, New Haven,: Human Relations Area Files PressGoogle Scholar
- Tsukahara T, Hombhanje FW, Lum JK, Hwaihwanje I, Masta A, Kaneko A, Kobayakawa T: Austronesian origin of the 27-bp deletion of the erythrocyte band 3 gene in East Sepik, Papua New Guinea inferred from mtDNA analysis. J Hum Genet. 2006, 51 (3): 244-248. 10.1007/s10038-005-0352-7.View ArticlePubMedGoogle Scholar
- Hill AV, Bowden DK, Trent RJ, Higgs DR, Oppenheimer SJ, Thein SL, Mickleson KN, Weatherall DJ, Clegg JB: Melanesians and Polynesians share a unique alpha-thalassemia mutation. Am J Hum Genet. 1985, 37 (3): 571-580.PubMed CentralPubMedGoogle Scholar
- Booth PB, Albrey JA, Whittaker J, Sanger R: Gerbich blood group system: a useful genetic marker in certain Melanesians of Papua and New Guinea. Nature. 1970, 228 (5270): 462-View ArticlePubMedGoogle Scholar
- Mgone CS, Koki G, Paniu MM, Kono J, Bhatia KK, Genton B, Alexander ND, Alpers MP: Occurrence of the erythrocyte band 3 (AE1) gene deletion in relation to malaria endemicity in Papua New Guinea. Trans R Soc Trop Med Hyg. 1996, 90 (3): 228-231. 10.1016/S0035-9203(96)90223-0.View ArticlePubMedGoogle Scholar
- Kimura M, Shimizu Y, Settheetham-Ishida W, Soemantri A, Tiwawech D, Romphruk A, Duangchan P, Ishida T: Twenty-seven base pair deletion in erythrocyte band 3 protein gene responsible for Southeast Asian ovalocytosis is not common among Southeast Asians. Hum Biol. 1998, 70 (6): 993-1000.PubMedGoogle Scholar
- Kimura M, Tamam M, Soemantri A, Nakazawa M, Ataka Y, Ohtsuka R, Ishida T: Distribution of a 27-bp deletion in the band 3 gene in South Pacific islanders. J Hum Genet. 2003, 48 (12): 642-645. 10.1007/s10038-003-0093-4.View ArticlePubMedGoogle Scholar
- Kayser M, Choi Y, van Oven M, Mona S, Brauer S, Trent RJ, Suarkia D, Schiefenhovel W, Stoneking M: The impact of the Austronesian expansion: evidence from mtDNA and Y chromosome diversity in the Admiralty Islands of Melanesia. Mol Biol Evol. 2008, 25 (7): 1362-1374. 10.1093/molbev/msn078.View ArticlePubMedGoogle Scholar
- Lansing JS, Cox MP, Downey SS, Gabler BM, Hallmark B, Karafet TM, Norquest P, Schoenfelder JW, Sudoyo H, Watkins JC: Coevolution of languages and genes on the island of Sumba, eastern Indonesia. Proc Natl Acad Sci USA. 2007, 104 (41): 16022-16026. 10.1073/pnas.0704451104.PubMed CentralView ArticlePubMedGoogle Scholar
- Mona S, Grunz KE, Brauer S, Pakendorf B, Castri L, Sudoyo H, Marzuki S, Barnes RH, Schmidtke J, Stoneking M: Genetic admixture history of Eastern Indonesia as revealed by Y-chromosome and mitochondrial DNA analysis. Mol Biol Evol. 2009, 26 (8): 1865-1877. 10.1093/molbev/msp097.View ArticlePubMedGoogle Scholar
- Serjeantson SW, Kirk RL, Booth PB: Linguistic and Genetic Differentiation in New-Guinea. J Hum Evol. 1983, 12 (1): 77-92. 10.1016/S0047-2484(83)80014-1.View ArticleGoogle Scholar
- Giles E, Ogan E, Steinberg AG: Gammaglobulin factors (Gm and Inv) in New Guinea: anthropological significance. Science. 1965, 150 (700): 1158-1160. 10.1126/science.150.3700.1158.View ArticlePubMedGoogle Scholar
- Kayser M, Brauer S, Weiss G, Schiefenhovel W, Underhill P, Shen P, Oefner P, Tommaseo-Ponzetta M, Stoneking M: Reduced y-chromosome, but not mitochondrial DNA, diversity in human populations from west new Guinea. Am J Hum Genet. 2003, 72 (2): 281-302. 10.1086/346065.PubMed CentralView ArticlePubMedGoogle Scholar
- Mona S, Tommaseo-Ponzetta M, Brauer S, Sudoyo H, Marzuki S, Kayser M: Patterns of Y-chromosome Diversity Intersect with the Trans-New Guinea Hypothesis. Mol Biol Evol. 2007Google Scholar
- Hunley K, Dunn M, Lindstrom E, Reesink G, Terrill A, Healy ME, Koki G, Friedlaender FR, Friedlaender JS: Genetic and linguistic coevolution in Northern Island Melanesia. PLoS genetics. 2008, 4 (10): e1000239-10.1371/journal.pgen.1000239.PubMed CentralView ArticlePubMedGoogle Scholar
- Donohue M, Denham T: Farming and Language in Island Southeast Asia Reframing Austronesian History. Curr Anthropol. 2010, 51 (2): 223-256. 10.1086/650991.View ArticleGoogle Scholar
- Tibbetts GR: A study of the Arabic texts containing material on South-East Asia. 1979, Leiden: BrillGoogle Scholar
- Cox MP: The Genetic Environment of Melanesia: Clines, Clusters and Contact. Population Genetics Research Progress. Edited by: Koven VT. 2008, Nova Science Publishers, 45-83.Google Scholar
- Wall JD, Cox MP, Mendez FL, Woerner A, Severson T, Hammer MF: A novel DNA sequence database for analyzing human demographic history. Genome Res. 2008, 18 (8): 1354-1361. 10.1101/gr.075630.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Jakobsson M, Scholz SW, Scheet P, Gibbs JR, VanLiere JM, Fung HC, Szpiech ZA, Degnan JH, Wang K, Guerreiro R: Genotype, haplotype and copy-number variation in worldwide human populations. Nature. 2008, 451 (7181): 998-1003. 10.1038/nature06742.View ArticlePubMedGoogle Scholar
- Chakraborty R, Kamboh MI, Nwankwo M, Ferrell RE: Caucasian genes in American blacks: new data. Am J Hum Genet. 1992, 50 (1): 145-155.PubMed CentralPubMedGoogle Scholar
- Weir BS: Genetic data analysis II : methods for discrete population genetic data. 1996, Sunderland, Mass.: Sinauer AssociatesGoogle Scholar
- Excoffier L, Lischer HE: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10 (3): 564-567. 10.1111/j.1755-0998.2010.02847.x.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMed CentralPubMedGoogle Scholar
- Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9 (5): 1322-1332. 10.1111/j.1755-0998.2009.02591.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Patterson N, Price AL, Reich D: Population structure and eigenanalysis. PLoS genetics. 2006, 2 (12): e190-10.1371/journal.pgen.0020190.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.