The genomic diversity of Taiwanese Austronesian groups: Implications for the "Into- and Out-of-Taiwan" models

The origin and dispersal of the Austronesian language family, one of the largest and most widespread in the world, have long attracted the attention of linguists, archaeologists, and geneticists. Even though there is a growing consensus that Taiwan is the source of the spread of Austronesian languages, little is known about the migration patterns of the early Austronesians who settled in and left Taiwan, i.e. the "Into-Taiwan" and "out-of-Taiwan" events. In particular, the genetic diversity and structure within Taiwan and how this relates to the into-/out-of-Taiwan events are largely unexplored, primarily because most genomic studies have largely utilized data from just two of the 16 recognized Highland Austronesian groups in Taiwan. In this study, we generated the largest genome-wide data set of Taiwanese Austronesians to date, including six Highland groups and one Lowland group from across the island and two Taiwanese Han groups. We identified fine-scale genomic structure in Taiwan, inferred the ancestry profile of the ancestors of Austronesians, and found that the southern Taiwanese Austronesians show excess genetic affinities with the Austronesians outside of Taiwan. Our findings thus shed new light on the Into- and Out-of-Taiwan dispersals.

Revisiting informed consent in forensic genomics in light of current technologies and the times

Informed consent is based on basic ethical principles that should be considered when conducting biomedical and behavioral research involving human subjects. These principles-respect, beneficence, and justice-form the foundations of informed consent which in itself is grounded on three fundamental elements: information, comprehension, and voluntary participation. While informed consent has focused on human subjects and research, the practice has been adopted willingly in the forensic science arena primarily to acquire reference samples from family members to assist in identifying missing persons. With advances in molecular biology technologies, data mining, and access to metadata, it is important to assess whether the past informed consent process and in particular associated risks are concomitant with these increased capabilities. Given the state-of-the-art, areas in which informed consent may need to be modified and augmented are as follows: reference samples from family members in missing persons or unidentified human remains cases; targeted analysis of an individual(s) during forensic genetic genealogy cases to reduce an investigative burden; donors who provide their samples for validation studies (to include population studies and entry into databases that would be applied to forensic statistical calculations) to support implementation of procedures and operations of the forensic laboratory; family members that may contribute samples or obtain genetic information from a molecular autopsy; and use of medical and other acquired samples that could be informative for identification purposes. The informed consent process should cover (1) purpose for collection of samples; (2) process to analyze the samples (to include type of data); (3) benefits (to donor, target, family, community, etc. as applicable); (4) risks (to donor, target, family, community, etc. as applicable); (5) access to data/reports by the donor; (6) sample disposition; (7) removal of data process (i.e., expungement); (8) process to ask questions/assessment of comprehension; (9) follow-up processes; and (10) voluntary, signed, and dated consent. Issues surrounding these topics are discussed with an emphasis on addressing risk factors. Addressing informed consent will allow human subjects to make decisions voluntarily and with autonomy as well as secure the use of samples for intended use.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

Shared paternal ancestry of Han, Tai-Kadai-speaking, and Austronesian-speaking populations as revealed by the high resolution phylogeny of O1a-M119 and distribution of its sub-lineages within China

OBJECTIVES

The aim of this research was to explore the origin, diversification, and demographic history of O1a-M119 over the past 10,000 years, as well as its role during the formation of East Asian and Southeast Asian populations, particularly the Han, Tai-Kadai-speaking, and Austronesian-speaking populations.

MATERIALS AND METHODS

Y-chromosome sequences (n = 141) of the O1a-M119 lineage, including 17 newly generated in this study, were used to reconstruct a revised phylogenetic tree with age estimates, and identify sub-lineages. The geographic distribution of 12 O1a-M119 sub-lineages was summarized, based on 7325 O1a-M119 individuals identified among 60,009 Chinese males.

RESULTS

A revised phylogenetic tree, age estimation, and distribution maps indicated continuous expansion of haplogroup O1a-M119 over the past 10,000 years, and differences in demographic history across geographic regions. We propose several sub-lineages of O1a-M119 as founding paternal lineages of Han, Tai-Kadai-speaking, and Austronesian-speaking populations. The sharing of several young O1a-M119 sub-lineages with expansion times less than 6000 years between these three population groups supports a partial common ancestry for them in the Neolithic Age; however, the paternal genetic divergence pattern is much more complex than previous hypotheses based on ethnology, archeology, and linguistics.

DISCUSSION

Our analyses contribute to a better understanding of the demographic history of O1a-M119 sub-lineages over the past 10,000 years during the emergence of Han, Austronesians, Tai-Kadai-speaking populations. The data described in this study will assist in understanding of the history of Han, Tai-Kadai-speaking, and Austronesian-speaking populations from ethnology, archeology, and linguistic perspectives in the future.

Ancient DNA from Guam and the peopling of the Pacific

Humans reached the Mariana Islands in the western Pacific by ∼3,500 y ago, contemporaneous with or even earlier than the initial peopling of Polynesia. They crossed more than 2,000 km of open ocean to get there, whereas voyages of similar length did not occur anywhere else until more than 2,000 y later. Yet, the settlement of Polynesia has received far more attention than the settlement of the Marianas. There is uncertainty over both the origin of the first colonizers of the Marianas (with different lines of evidence suggesting variously the Philippines, Indonesia, New Guinea, or the Bismarck Archipelago) as well as what, if any, relationship they might have had with the first colonizers of Polynesia. To address these questions, we obtained ancient DNA data from two skeletons from the Ritidian Beach Cave Site in northern Guam, dating to ∼2,200 y ago. Analyses of complete mitochondrial DNA genome sequences and genome-wide SNP data strongly support ancestry from the Philippines, in agreement with some interpretations of the linguistic and archaeological evidence, but in contradiction to results based on computer simulations of sea voyaging. We also find a close link between the ancient Guam skeletons and early Lapita individuals from Vanuatu and Tonga, suggesting that the Marianas and Polynesia were colonized from the same source population, and raising the possibility that the Marianas played a role in the eventual settlement of Polynesia.

Paternal gene pool of Malays in Southeast Asia and its applications for the early expansion of Austronesians

OBJECTIVES

The origin and differentiation of Austronesian populations and their languages have long fascinated linguists, archeologists, and geneticists. However, the founding process of Austronesians and when they separated from their close relatives, such as the Daic and Austro-Asiatic populations in the mainland of Asia, remain unclear. In this study, we explored the paternal origin of Malays in Southeast Asia and the early differentiation of Austronesians.

MATERIALS AND METHODS

We generated whole Y-chromosome sequences of 50 Malays and co-analyzed 200 sequences from other Austronesians and related populations. We generated a revised phylogenetic tree with time estimation.

RESULTS

We identified six founding paternal lineages among the studied Malays samples. These founding lineages showed a surprisingly coincident expansion age at 5000 to 6000 years ago. We also found numerous mostly close related samples of the founding lineages of Malays among populations from Mainland of Asia.

CONCLUSION

Our analyses provided a refined phylogenetic resolution for the dominant paternal lineages of Austronesians found by previous studies. We suggested that the co-expansion of numerous founding paternal lineages corresponds to the initial differentiation of the most recent common ancestor of modern Austronesians. The splitting time and divergence pattern in perspective of paternal Y-chromosome evidence are highly consistent with the previous theories of ethnologists, linguists, and archeologists.

Admixture and natural selection shaped genomes of an Austronesian-speaking population in the Solomon Islands

People in the Solomon Islands today are considered to have derived from Asian- and Papuan-related ancestors. Papuan-related ancestors colonized Near Oceania about 47,000 years ago, and Asian-related ancestors were Austronesian (AN)-speaking population, called Lapita, who migrated from Southeast Asia about 3,500 years ago. These two ancestral populations admixed in Near Oceania before the expansion of Lapita people into Remote Oceania. To understand the impact of the admixture on the adaptation of AN-speaking Melanesians in Near Oceania, we performed the genome-wide single nucleotide polymorphism (SNP) analysis of 21 individuals from Munda, the main town of the New Georgia Islands in the western Solomon Islands. Population samples from Munda were genetically similar to other Solomon Island population samples. The analysis of genetic contribution from the two different ancestries to the Munda genome revealed significantly higher proportions of Asian- and Papuan-related ancestries in the region containing the annexin A1 (ANXA1) gene (Asian component > 82.6%) and in the human leukocyte antigen (HLA) class II region (Papuan component > 85.4%), respectively. These regions were suspected to have undergone natural selection since the time of admixture. Our results suggest that admixture had affected adaptation of AN-speaking Melanesians in the Solomon Islands.

Association study of CREBRF missense variant (rs373863828:G > A; p.Arg457Gln) with levels of serum lipid profile in the Pacific populations

BACKGROUND

A missense variant (rs373863828:G > A; p.Arg457Gln) of the CREBRF gene is strongly associated with a higher body mass index (BMI; kg/m²) in Polynesian populations. This variant has also been reported to be associated with lower total cholesterol in Samoans.

AIM

The aim of this study is to examine the association of rs373863828:G > A with levels of serum lipids in four Pacific populations.

METHODS

A total of 613 adult subjects were recruited from Tonga (Polynesians) and the Solomon Islands (Melanesians and Micronesians). Multiple regression analyses adjusted for age and sex were performed to examine the association of rs373863828 with levels of serum lipids in each population.

RESULTS

A significant association of rs373863828:G > A with lower level of HDL-cholesterol was detected in the Tonga population (β = -3.32 and p-value = 0.030). The expected change in HDL-cholesterol with respect to a single copy of the rs373863828-A allele was 3.32 mg/dL. However, the association between rs373863828-A and lower levels of HDL-cholesterol was not significant after further adjustment for BMI in the Tonga population (β = -2.32 and p-value = 0.13).

CONCLUSIONS

The rs373863828-A allele may not directly affect the level of serum HDL-cholesterol independent of BMI. To confirm the present findings, association studies with large sample sizes and functional analyses are required.

The Gateway from Near into Remote Oceania: New Insights from Genome-Wide Data

A widely accepted two-wave scenario of human settlement of Oceania involves the first out-of-Africa migration circa 50,000 years ago (ya), and the more recent Austronesian expansion, which reached the Bismarck Archipelago by 3,450 ya. Whereas earlier genetic studies provided evidence for extensive sex-biased admixture between the incoming and the indigenous populations, some archaeological, linguistic, and genetic evidence indicates a more complicated picture of settlement. To study regional variation in Oceania in more detail, we have compiled a genome-wide data set of 823 individuals from 72 populations (including 50 populations from Oceania) and over 620,000 autosomal single nucleotide polymorphisms (SNPs). We show that the initial dispersal of people from the Bismarck Archipelago into Remote Oceania occurred in a "leapfrog" fashion, completely by-passing the main chain of the Solomon Islands, and that the colonization of the Solomon Islands proceeded in a bidirectional manner. Our results also support a divergence between western and eastern Solomons, in agreement with the sharp linguistic divide known as the Tryon-Hackman line. We also report substantial post-Austronesian gene flow across the Solomons. In particular, Santa Cruz (in Remote Oceania) exhibits extraordinarily high levels of Papuan ancestry that cannot be explained by a simple bottleneck/founder event scenario. Finally, we use simulations to show that discrepancies between different methods for dating admixture likely reflect different sensitivities of the methods to multiple admixture events from the same (or similar) sources. Overall, this study points to the importance of fine-scale sampling to understand the complexities of human population history.

Investigating the origins of eastern Polynesians using genome-wide data from the Leeward Society Isles

The debate concerning the origin of the Polynesian speaking peoples has been recently reinvigorated by genetic evidence for secondary migrations to western Polynesia from the New Guinea region during the 2nd millennium BP. Using genome-wide autosomal data from the Leeward Society Islands, the ancient cultural hub of eastern Polynesia, we find that the inhabitants' genomes also demonstrate evidence of this episode of admixture, dating to 1,700-1,200 BP. This supports a late settlement chronology for eastern Polynesia, commencing ~1,000 BP, after the internal differentiation of Polynesian society. More than 70% of the autosomal ancestry of Leeward Society Islanders derives from Island Southeast Asia with the lowland populations of the Philippines as the single largest potential source. These long-distance migrants into Polynesia experienced additional admixture with northern Melanesians prior to the secondary migrations of the 2nd millennium BP. Moreover, the genetic diversity of mtDNA and Y chromosome lineages in the Leeward Society Islands is consistent with linguistic evidence for settlement of eastern Polynesia proceeding from the central northern Polynesian outliers in the Solomon Islands. These results stress the complex demographic history of the Leeward Society Islands and challenge phylogenetic models of cultural evolution predicated on eastern Polynesia being settled from Samoa.

Contrasting maternal and paternal genetic variation of hunter-gatherer groups in Thailand

The Maniq and Mlabri are the only recorded nomadic hunter-gatherer groups in Thailand. Here, we sequenced complete mitochondrial (mt) DNA genomes and ~2.364 Mbp of non-recombining Y chromosome (NRY) to learn more about the origins of these two enigmatic populations. Both groups exhibited low genetic diversity compared to other Thai populations, and contrasting patterns of mtDNA and NRY diversity: there was greater mtDNA diversity in the Maniq than in the Mlabri, while the converse was true for the NRY. We found basal uniparental lineages in the Maniq, namely mtDNA haplogroups M21a, R21 and M17a, and NRY haplogroup K. Overall, the Maniq are genetically similar to other negrito groups in Southeast Asia. By contrast, the Mlabri haplogroups (B5a1b1 for mtDNA and O1b1a1a1b and O1b1a1a1b1a1 for the NRY) are common lineages in Southeast Asian non-negrito groups, and overall the Mlabri are genetically similar to their linguistic relatives (Htin and Khmu) and other groups from northeastern Thailand. In agreement with previous studies of the Mlabri, our results indicate that the Malbri do not directly descend from the indigenous negritos. Instead, they likely have a recent origin (within the past 1,000 years) by an extreme founder event (involving just one maternal and two paternal lineages) from an agricultural group, most likely the Htin or a closely-related group.

Adaptation of human skin color in various populations

BACKGROUND

Skin color is a well-recognized adaptive trait and has been studied extensively in humans. Understanding the genetic basis of adaptation of skin color in various populations has many implications in human evolution and medicine.

DISCUSSION

Impressive progress has been made recently to identify genes associated with skin color variation in a wide range of geographical and temporal populations. In this review, we discuss what is currently known about the genetics of skin color variation. We enumerated several cases of skin color adaptation in global modern humans and archaic hominins, and illustrated why, when, and how skin color adaptation occurred in different populations. Finally, we provided a summary of the candidate loci associated with pigmentation, which could be a valuable reference for further evolutionary and medical studies.

CONCLUSION

Previous studies generally indicated a complex genetic mechanism underlying the skin color variation, expanding our understanding of the role of population demographic history and natural selection in shaping genetic and phenotypic diversity in humans. Future work is needed to dissect the genetic architecture of skin color adaptation in numerous ethnic minority groups around the world, which remains relatively obscure compared with that of major continental groups, and to unravel the exact genetic basis of skin color adaptation.

Phylogeography of Y-chromosome haplogroup O3a2b2-N6 reveals patrilineal traces of Austronesian populations on the eastern coastal regions of Asia

Austronesian diffusion is considered one of the greatest dispersals in human history; it led to the peopling of an extremely vast region, ranging from Madagascar in the Indian Ocean to Easter Island in Remote Oceania. The Y-chromosome haplogroup O3a2b*-P164(xM134), a predominant paternal lineage of Austronesian populations, is found at high frequencies in Polynesian populations. However, the internal phylogeny of this haplogroup remains poorly investigated. In this study, we analyzed -seventeen Y-chromosome sequences of haplogroup O3a2b*-P164(xM134) and generated a revised phylogenetic tree of this lineage based on 310 non-private Y-chromosome polymorphisms. We discovered that all available O3a2b*-P164(xM134) samples belong to the newly defined haplogroup O3a2b2-N6 and samples from Austronesian populations belong to the sublineage O3a2b2a2-F706. Additionally, we genotyped a series of Y-chromosome polymorphisms in a large collection of samples from China. We confirmed that the sublineage O3a2b2a2b-B451 is unique to Austronesian populations. We found that O3a2b2-N6 samples are widely distributed on the eastern coastal regions of Asia, from Korea to Vietnam. Furthermore, we propose- that the O3a2b2a2b-B451 lineage represents a genetic connection between ancestors of Austronesian populations and ancient populations in North China, where foxtail millet was domesticated about 11,000 years ago. The large number of newly defined Y-chromosome polymorphisms and the revised phylogenetic tree of O3a2b2-N6 will be helpful to explore the origin of proto-Austronesians and the early diffusion process of Austronesian populations.

Fine-Scale Human Population Structure in Southern Africa Reflects Ecogeographic Boundaries

Recent genetic studies have established that the KhoeSan populations of southern Africa are distinct from all other African populations and have remained largely isolated during human prehistory until ∼2000 years ago. Dozens of different KhoeSan groups exist, belonging to three different language families, but very little is known about their population history. We examine new genome-wide polymorphism data and whole mitochondrial genomes for >100 South Africans from the ≠Khomani San and Nama populations of the Northern Cape, analyzed in conjunction with 19 additional southern African populations. Our analyses reveal fine-scale population structure in and around the Kalahari Desert. Surprisingly, this structure does not always correspond to linguistic or subsistence categories as previously suggested, but rather reflects the role of geographic barriers and the ecology of the greater Kalahari Basin. Regardless of subsistence strategy, the indigenous Khoe-speaking Nama pastoralists and the N|u-speaking ≠Khomani (formerly hunter-gatherers) share ancestry with other Khoe-speaking forager populations that form a rim around the Kalahari Desert. We reconstruct earlier migration patterns and estimate that the southern Kalahari populations were among the last to experience gene flow from Bantu speakers, ∼14 generations ago. We conclude that local adoption of pastoralism, at least by the Nama, appears to have been primarily a cultural process with limited genetic impact from eastern Africa.

Resolving the ancestry of Austronesian-speaking populations

There are two very different interpretations of the prehistory of Island Southeast Asia (ISEA), with genetic evidence invoked in support of both. The "out-of-Taiwan" model proposes a major Late Holocene expansion of Neolithic Austronesian speakers from Taiwan. An alternative, proposing that Late Glacial/postglacial sea-level rises triggered largely autochthonous dispersals, accounts for some otherwise enigmatic genetic patterns, but fails to explain the Austronesian language dispersal. Combining mitochondrial DNA (mtDNA), Y-chromosome and genome-wide data, we performed the most comprehensive analysis of the region to date, obtaining highly consistent results across all three systems and allowing us to reconcile the models. We infer a primarily common ancestry for Taiwan/ISEA populations established before the Neolithic, but also detected clear signals of two minor Late Holocene migrations, probably representing Neolithic input from both Mainland Southeast Asia and South China, via Taiwan. This latter may therefore have mediated the Austronesian language dispersal, implying small-scale migration and language shift rather than large-scale expansion.

Early Lapita skeletons from Vanuatu show Polynesian craniofacial shape: Implications for Remote Oceanic settlement and Lapita origins

With a cultural and linguistic origin in Island Southeast Asia the Lapita expansion is thought to have led ultimately to the Polynesian settlement of the east Polynesian region after a time of mixing/integration in north Melanesia and a nearly 2,000-y pause in West Polynesia. One of the major achievements of recent Lapita research in Vanuatu has been the discovery of the oldest cemetery found so far in the Pacific at Teouma on the south coast of Efate Island, opening up new prospects for the biological definition of the early settlers of the archipelago and of Remote Oceania in general. Using craniometric evidence from the skeletons in conjunction with archaeological data, we discuss here four debated issues: the Lapita-Asian connection, the degree of admixture, the Lapita-Polynesian connection, and the question of secondary population movement into Remote Oceania.

MC1R diversity in Northern Island Melanesia has not been constrained by strong purifying selection and cannot explain pigmentation phenotype variation in the region

BACKGROUND

Variation in human skin pigmentation evolved in response to the selective pressure of ultra-violet radiation (UVR). Selection to maintain darker skin in high UVR environments is expected to constrain pigmentation phenotype and variation in pigmentation loci. Consistent with this hypothesis, the gene MC1R exhibits reduced diversity in African populations from high UVR regions compared to low-UVR non-African populations. However, MC1R diversity in non-African populations that have evolved under high-UVR conditions is not well characterized.

METHODS

In order to test the hypothesis that MC1R variation has been constrained in Melanesians the coding region of the MC1R gene was sequenced in 188 individuals from Northern Island Melanesia. The role of purifying selection was assessed using a modified McDonald Kreitman's test. Pairwise FST was calculated between Melanesian populations and populations from the 1000 Genomes Project. The SNP rs2228479 was genotyped in a larger sample (n = 635) of Melanesians and tested for associations with skin and hair pigmentation.

RESULTS

We observe three nonsynonymous and two synonymous mutations. A modified McDonald Kreitman's test failed to detect a significant signal of purifying selection. Pairwise FST values calculated between the four islands sampled here indicate little regional substructure in MC1R. When compared to African, European, East and South Asian populations, Melanesians do not exhibit reduced population divergence (measured as FST) or a high proportion of haplotype sharing with Africans, as one might expect if ancestral haplotypes were conserved across high UVR populations in and out of Africa. The only common nonsynonymous polymorphism observed, rs2228479, is not significantly associated with skin or hair pigmentation in a larger sample of Melanesians.

CONCLUSIONS

The pattern of sequence diversity here does not support a model of strong selective constraint on MC1R in Northern Island Melanesia This absence of strong constraint, as well as the recent population history of the region, may explain the observed frequencies of the derived rs2228479 allele. These results emphasize the complex genetic architecture of pigmentation phenotypes, which are controlled by multiple, possibly interacting loci. They also highlight the role that population history can play in influencing phenotypic diversity in the absence of strong natural selection.

Denisovan Ancestry in East Eurasian and Native American Populations

Although initial studies suggested that Denisovan ancestry was found only in modern human populations from island Southeast Asia and Oceania, more recent studies have suggested that Denisovan ancestry may be more widespread. However, the geographic extent of Denisovan ancestry has not been determined, and moreover the relationship between the Denisovan ancestry in Oceania and that elsewhere has not been studied. Here we analyze genome-wide single nucleotide polymorphism data from 2,493 individuals from 221 worldwide populations, and show that there is a widespread signal of a very low level of Denisovan ancestry across Eastern Eurasian and Native American (EE/NA) populations. We also verify a higher level of Denisovan ancestry in Oceania than that in EE/NA; the Denisovan ancestry in Oceania is correlated with the amount of New Guinea ancestry, but not the amount of Australian ancestry, indicating that recent gene flow from New Guinea likely accounts for signals of Denisovan ancestry across Oceania. However, Denisovan ancestry in EE/NA populations is equally correlated with their New Guinea or their Australian ancestry, suggesting a common source for the Denisovan ancestry in EE/NA and Oceanian populations. Our results suggest that Denisovan ancestry in EE/NA is derived either from common ancestry with, or gene flow from, the common ancestor of New Guineans and Australians, indicating a more complex history involving East Eurasians and Oceanians than previously suspected.

Genome-wide insights into the genetic history of human populations

Although mtDNA and the non-recombining Y chromosome (NRY) studies continue to provide valuable insights into the genetic history of human populations, recent technical, methodological and computational advances and the increasing availability of large-scale, genome-wide data from contemporary human populations around the world promise to reveal new aspects, resolve finer points, and provide a more detailed look at our past demographic history. Genome-wide data are particularly useful for inferring migrations, admixture, and fine structure, as well as for estimating population divergence and admixture times and fluctuations in effective population sizes. In this review, we highlight some of the stories that have emerged from the analyses of genome-wide SNP genotyping data concerning the human history of Southern Africa, India, Oceania, Island South East Asia, Europe and the Americas and comment on possible future study directions. We also discuss advantages and drawbacks of using SNP-arrays, with a particular focus on the ascertainment bias, and ways to circumvent it.

Human settlement history between Sunda and Sahul: a focus on East Timor (Timor-Leste) and the Pleistocenic mtDNA diversity

BACKGROUND

Distinct, partly competing, "waves" have been proposed to explain human migration in(to) today's Island Southeast Asia and Australia based on genetic (and other) evidence. The paucity of high quality and high resolution data has impeded insights so far. In this study, one of the first in a forensic environment, we used the Ion Torrent Personal Genome Machine (PGM) for generating complete mitogenome sequences via stand-alone massively parallel sequencing and describe a standard data validation practice.

RESULTS

In this first representative investigation on the mitochondrial DNA (mtDNA) variation of East Timor (Timor-Leste) population including >300 individuals, we put special emphasis on the reconstruction of the initial settlement, in particular on the previously poorly resolved haplogroup P1, an indigenous lineage of the Southwest Pacific region. Our results suggest a colonization of southern Sahul (Australia) >37 kya, limited subsequent exchange, and a parallel incubation of initial settlers in northern Sahul (New Guinea) followed by westward migrations <28 kya.

CONCLUSIONS

The temporal proximity and possible coincidence of these latter dispersals, which encompassed autochthonous haplogroups, with the postulated "later" events of (South) East Asian origin pinpoints a highly dynamic migratory phase.

Genetic structure among Fijian island populations

We examined nine Y chromosome short tandem repeats (Y-STRs) and the mitochondrial DNA (mtDNA) hypervariable segment 1 region in the Fijian island populations of Viti Levu, Vanua Levu, Kadavu, the Lau islands and Rotuma. We found significant genetic structure among these populations for the Y-STRs, both with and without the Rotumans, but not for the mtDNA. We also found that all five populations exhibited the sex-biased admixture associated with areas settled by Austronesian-speaking people, with paternal lineages more strongly associated with Melanesian populations and maternal lineages more strongly associated with Polynesian populations. We also found that the Rotumans in the north and the Lau Islanders in the east were genetically more similar to Polynesian populations than were the other Fijians, but only for the mtDNA. For the Y-STRs, the Rotumans and the Lau Islanders were genetically as similar to Melanesian populations as were the other three populations. Of the five populations, the Rotumans were the most different in almost every regard. Although past genetic studies treated the Fijians as being genetically homogenous despite known geographic, phenotypic, cultural and linguistic variation, our findings show significant genetic variation and a need for a closer examination of individual island populations within Fiji, particularly the Rotumans, in order to better understand the process of the peopling of Fiji and of the surrounding regions.

Human paternal and maternal demographic histories: insights from high-resolution Y chromosome and mtDNA sequences

Background

Comparisons of maternally-inherited mitochondrial DNA (mtDNA) and paternally-inherited non-recombining Y chromosome (NRY) variation have provided important insights into the impact of sex-biased processes (such as migration, residence pattern, and so on) on human genetic variation. However, such comparisons have been limited by the different molecular methods typically used to assay mtDNA and NRY variation (for example, sequencing hypervariable segments of the control region for mtDNA vs. genotyping SNPs and/or STR loci for the NRY). Here, we report a simple capture array method to enrich Illumina sequencing libraries for approximately 500 kb of NRY sequence, which we use to generate NRY sequences from 623 males from 51 populations in the CEPH Human Genome Diversity Panel (HGDP). We also obtained complete mtDNA genome sequences from the same individuals, allowing us to compare maternal and paternal histories free of any ascertainment bias.

Results

We identified 2,228 SNPs in the NRY sequences and 2,163 SNPs in the mtDNA sequences. Our results confirm the controversial assertion that genetic differences between human populations on a global scale are bigger for the NRY than for mtDNA, although the differences are not as large as previously suggested. More importantly, we find substantial regional variation in patterns of mtDNA versus NRY variation. Model-based simulations indicate very small ancestral effective population sizes (<100) for the out-of-Africa migration as well as for many human populations. We also find that the ratio of female effective population size to male effective population size (N_f/N_m) has been greater than one throughout the history of modern humans, and has recently increased due to faster growth in N_f than N_m.

Conclusions

The NRY and mtDNA sequences provide new insights into the paternal and maternal histories of human populations, and the methods we introduce here should be widely applicable for further such studies.

Coevolution of languages and genes

The evolution of languages shares certain characteristics with that of genes, such as the predominantly vertical line of transmission and the retention of traces of past events such as contact. Thus, studies of language phylogenies and their correlations with genetic phylogenies can enrich our understanding of human prehistory, while insights gained from genetic studies of past population contact can help shed light on the processes underlying language contact and change. As demonstrated by recent research, these evolutionary processes are more complex than simple models of gene-language coevolution predict, with linguistic boundaries only occasionally functioning as barriers to gene flow. More frequently, admixture takes place irrespective of linguistic differences, but with a detectable impact of contact-induced changes in the languages concerned.

Recent developments in the genetic history of East Asia and Oceania

Recent developments in our understanding of the genetic history of Asia and Oceania have been driven by technological advances. Specifically, our understanding of the past has been augmented by: genome sequences from ancient hominins and ancient modern humans; more comprehensive studies of existing populations (e.g., complete mtDNA genome sequences and genome-wide data) and the development of new statistics and analytical methods to interpret the abundance of new data. We review some of the new discoveries since we entered the age of archaic and modern genomics and how they have changed our understanding of the settlement and subsequent population dynamics in Asia and the Pacific.

Taiwan Y-chromosomal DNA variation and its relationship with Island Southeast Asia

Background

Much of the data resolution of the haploid non-recombining Y chromosome (NRY) haplogroup O in East Asia are still rudimentary and could be an explanatory factor for current debates on the settlement history of Island Southeast Asia (ISEA). Here, 81 slowly evolving markers (mostly SNPs) and 17 Y-chromosomal short tandem repeats were used to achieve higher level molecular resolution. Our aim is to investigate if the distribution of NRY DNA variation in Taiwan and ISEA is consistent with a single pre-Neolithic expansion scenario from Southeast China to all ISEA, or if it better fits an expansion model from Taiwan (the OOT model), or whether a more complex history of settlement and dispersals throughout ISEA should be envisioned.

Results

We examined DNA samples from 1658 individuals from Vietnam, Thailand, Fujian, Taiwan (Han, plain tribes and 14 indigenous groups), the Philippines and Indonesia. While haplogroups O1a*-M119, O1a1*-P203, O1a2-M50 and O3a2-P201 follow a decreasing cline from Taiwan towards Western Indonesia, O2a1-M95/M88, O3a*-M324, O3a1c-IMS-JST002611 and O3a2c1a-M133 decline northward from Western Indonesia towards Taiwan. Compared to the Taiwan plain tribe minority groups the Taiwanese Austronesian speaking groups show little genetic paternal contribution from Han. They are also characterized by low Y-chromosome diversity, thus testifying for fast drift in these populations. However, in contrast to data provided from other regions of the genome, Y-chromosome gene diversity in Taiwan mountain tribes significantly increases from North to South.

Conclusion

The geographic distribution and the diversity accumulated in the O1a*-M119, O1a1*-P203, O1a2-M50 and O3a2-P201 haplogroups on one hand, and in the O2a1-M95/M88, O3a*-M324, O3a1c-IMS-JST002611 and O3a2c1a-M133 haplogroups on the other, support a pincer model of dispersals and gene flow from the mainland to the islands which likely started during the late upper Paleolithic, 18,000 to 15,000 years ago. The branches of the pincer contributed separately to the paternal gene pool of the Philippines and conjointly to the gene pools of Madagascar and the Solomon Islands. The North to South increase in diversity found for Taiwanese Austronesian speaking groups contrasts with observations based on mitochondrial DNA, thus hinting to a differentiated demographic history of men and women in these populations.

Maternal history of Oceania from complete mtDNA genomes: contrasting ancient diversity with recent homogenization due to the Austronesian expansion

Archaeology, linguistics, and existing genetic studies indicate that Oceania was settled by two major waves of migration. The first migration took place approximately 40 thousand years ago and these migrants, Papuans, colonized much of Near Oceania. Approximately 3.5 thousand years ago, a second expansion of Austronesian-speakers arrived in Near Oceania and the descendants of these people spread to the far corners of the Pacific, colonizing Remote Oceania. To assess the female contribution of these two human expansions to modern populations and to investigate the potential impact of other migrations, we obtained 1,331 whole mitochondrial genome sequences from 34 populations spanning both Near and Remote Oceania. Our results quantify the magnitude of the Austronesian expansion and demonstrate the homogenizing effect of this expansion on almost all studied populations. With regards to Papuan influence, autochthonous haplogroups support the hypothesis of a long history in Near Oceania, with some lineages suggesting a time depth of 60 thousand years, and offer insight into historical interpopulation dynamics. Santa Cruz, a population located in Remote Oceania, is an anomaly with extreme frequencies of autochthonous haplogroups of Near Oceanian origin; simulations to investigate whether this might reflect a pre-Austronesian versus Austronesian settlement of the island failed to provide unequivocal support for either scenario.

Human genetics of the Kula Ring: Y-chromosome and mitochondrial DNA variation in the Massim of Papua New Guinea

The island region at the southeastern-most tip of New Guinea and its inhabitants known as Massim are well known for a unique traditional inter-island trading system, called Kula or Kula Ring. To characterize the Massim genetically, and to evaluate the influence of the Kula Ring on patterns of human genetic variation, we analyzed paternally inherited Y-chromosome (NRY) and maternally inherited mitochondrial (mt) DNA polymorphisms in >400 individuals from this region. We found that the nearly exclusively Austronesian-speaking Massim people harbor genetic ancestry components of both Asian (AS) and Near Oceanian (NO) origin, with a proportionally larger NO NRY component versus a larger AS mtDNA component. This is similar to previous observations in other Austronesian-speaking populations from Near and Remote Oceania and suggests sex-biased genetic admixture between Asians and Near Oceanians before the occupation of Remote Oceania, in line with the Slow Boat from Asia hypothesis on the expansion of Austronesians into the Pacific. Contrary to linguistic expectations, Rossel Islanders, the only Papuan speakers of the Massim, showed a lower amount of NO genetic ancestry than their Austronesian-speaking Massim neighbors. For the islands traditionally involved in the Kula Ring, a significant correlation between inter-island travelling distances and genetic distances was observed for mtDNA, but not for NRY, suggesting more male- than female-mediated gene flow. As traditionally only males take part in the Kula voyages, this finding may indicate a genetic signature of the Kula Ring, serving as another example of how cultural tradition has shaped human genetic diversity.

Taiwanese aborigines: genetic heterogeneity and paternal contribution to Oceania

In the present study, for the first time, 293 Taiwanese aboriginal males from all nine major tribes (Ami, Atayal, Bunun, Rukai, Paiwan, Saisat, Puyuma, Tsou, Yami) were genotyped with 17 YSTR loci in a attend to reveal migrational patterns connected with the Austronesian expansion. We investigate the paternal genetic relationships of these Taiwanese aborigines to 42 Asia-Pacific reference populations, geographically selected to reflect various locations within the Austronesian domain. The Tsou and Puyuma tribes exhibit the lowest (0.1851) and the highest (0.5453) average total genetic diversity, respectively. Further, the fraction of unique haplotypes is also relatively high in the Puyuma (86.7%) and low in Tsou (33.3%) suggesting different demographic histories. Multidimensional scaling (MDS) and analysis of molecular variance (AMOVA) revealed several notable findings: 1) the Taiwan indigenous populations are highly diverse. In fact, the level of inter-population heterogeneity displayed by the Taiwanese aboriginal populations is close to that exhibited among all 51 Asia-Pacific populations examined; 2) the asymmetrical contribution of the Taiwanese aborigines to the Oceanic groups. Ami, Bunun and Saisiyat tribes exhibit the strongest paternal links to the Solomon and Polynesian island communities, whereas most of the remaining Taiwanese aboriginal groups are more genetically distant to these Oceanic inhabitants; 3) the present YSTR analyses does not reveal a strong paternal affinity of the nine Taiwanese tribes to their continental Asian neighbors. Overall, our current findings suggest that, perhaps, only a few of the tribes were involved in the migration out of Taiwan.

Early Austronesians: into and out of Taiwan

A Taiwan origin for the expansion of the Austronesian languages and their speakers is well supported by linguistic and archaeological evidence. However, human genetic evidence is more controversial. Until now, there had been no ancient skeletal evidence of a potential Austronesian-speaking ancestor prior to the Taiwan Neolithic ~6,000 years ago, and genetic studies have largely ignored the role of genetic diversity within Taiwan as well as the origins of Formosans. We address these issues via analysis of a complete mitochondrial DNA genome sequence of an ~8,000-year-old skeleton from Liang Island (located between China and Taiwan) and 550 mtDNA genome sequences from 8 aboriginal (highland) Formosan and 4 other Taiwanese groups. We show that the Liangdao Man mtDNA sequence is closest to Formosans, provides a link to southern China, and has the most ancestral haplogroup E sequence found among extant Austronesian speakers. Bayesian phylogenetic analysis allows us to reconstruct a history of early Austronesians arriving in Taiwan in the north ~6,000 years ago, spreading rapidly to the south, and leaving Taiwan ~4,000 years ago to spread throughout Island Southeast Asia, Madagascar, and Oceania.

Distribution of an allele associated with blond hair color across Northern Island Melanesia

Pigmentation of the skin, hair, and eyes is a complex trait controlled by multiple genetic loci. Recently a non-synonymous mutation in the pigmentation candidate gene TYRP1 was shown to be significantly associated with a blond-hair phenotype in populations from the Solomon Islands. The distribution of this mutation in the islands of Northern Island Melanesia, where the blondism phenotype is also prevalent, was unknown. Here, we present data describing the distribution of this allele in 550 individuals sampled from across this region, and test for associations between genotype at this locus and quantitatively measured skin and hair pigmentation phenotype. We report that the frequency of the 93C allele is notably lower than observed in the Solomons (0.12 vs. 0.26). The allele exhibits significant geographic heterogeneity across the islands sampled (χ(2)  = 108.4, P < 0.0001). It is observed at its highest frequencies on the islands of New Ireland and New Hanover, while being almost completely absent from the large island of New Britain. Using linear regression with age, sex, and island as covariates we report that, as in the Solomons, the 93C allele is significantly associated with a decrease in hair pigmentation but not skin pigmentation. We discuss the distribution of the 93C allele across the Southwest Pacific in light of its possible place of origin and dispersal.

Investigating sex-specific dynamics using uniparental markers: West New Guinea as a case study

Mitochondrial DNA (mtDNA) and Y chromosome (NRY) genetic markers have been often contrasted to investigate sex-specific dynamics. Traditionally, isolation by distance, intrapopulation genetic diversity and population differentiation are estimated from both markers and compared. Two possible sources of bias are often neglected. First, kilometric distances are frequently used as predictor of the connectivity between groups, hiding the role played by environmental features at a microgeographic scale. Second, the comparison of intrapopulation diversity and population differentiation between mtDNA and NRY is hampered by their different mutational mechanisms and rates. Here, we show how to account for these biases by analyzing from a different perspective a published dataset of eight West New Guinea (WNG) populations for which mtDNA control region sequences and seven linked NRY microsatellites had been typed. First, we modeled the connectivity among sampled populations by computing the number of days required to travel between groups. Then, we investigated the differences between the two sexes accounting for the molecular characteristics of the markers examined to obtain estimates on the product of the effective population size and the migration rate among demes (Nm). We achieved this goal by studying the shape of the gene genealogy at several sampling levels and using spatial explicit simulations. Both the direction and the rate of migration differ between male and females, with an Nm estimated to be >6 times higher in the latter under many evolutionary scenarios. We finally highlight the importance of applying metapopulation models when analyzing the genetic diversity of a species. We have applied the prediction of the sampling theory in a meta-population and we have corroborated our finding using spatial explicit simulations. Both approaches are fundamentally meant to deal with structured populations: we strongly believe in the importance of tacking structure into account when inferring the demographic history of a species.

A highly unstable recent mutation in human mtDNA

An A-to-G transition at position 16247 in the human mtDNA genome denotes haplogroup B4a1a1a and its sublineages. Informally known as the "Polynesian motif," this haplogroup has been widely used as a marker in Oceania of genetic affiliation with the Austronesian expansion. The 16247G allele has arisen only once in the human mtDNA phylogeny, about 7,000 thousand years ago, and is nearly fixed in Remote Oceania. We analyzed 536 complete mtDNA genome sequences from the Solomon Islands from haplogroup B4a1a1 and associated subhaplogroups and found multiple independent back mutations from 16247G to 16247A. We also find elevated levels of heteroplasmy at this position in samples with the 16247G allele, suggesting the ongoing occurrence of somatic back-mutations and/or transmission of heteroplasmy. Moreover, the G allele is predicted to introduce a novel stem-loop structure in the DNA sequence that may be structurally unfavorable, thereby accounting for the remarkable number of back-mutations observed at the 16247G allele in this short evolutionary time span. More generally, haplogroup-calling scripts result in inaccurate haplogroup calls involving the back-mutation and need to be supplemented with other types of analyses; this may be true for other mtDNA lineages because no other lineage has been investigated to the same extent (over 500 complete mtDNA sequences).

Melanesian blond hair is caused by an amino acid change in TYRP1

Naturally blond hair is rare in humans and found almost exclusively in Europe and Oceania. Here, we identify an arginine-to-cysteine change at a highly conserved residue in tyrosinase-related protein 1 (TYRP1) as a major determinant of blond hair in Solomon Islanders. This missense mutation is predicted to affect catalytic activity of TYRP1 and causes blond hair through a recessive mode of inheritance. The mutation is at a frequency of 26% in the Solomon Islands, is absent outside of Oceania, represents a strong common genetic effect on a complex human phenotype, and highlights the importance of examining genetic associations worldwide.

MtDNA SNP multiplexes for efficient inference of matrilineal genetic ancestry within Oceania

Human mitochondrial DNA (mtDNA) is a convenient marker for tracing matrilineal bio-geographic ancestry and is widely applied in forensic, genealogical and anthropological studies. In forensic applications, DNA-based ancestry inference can be useful for finding unknown suspects by concentrating police investigations in cases where autosomal STR profiling was unable to provide a match, or can help provide clues in missing person identification. Although multiplexed mtDNA single nucleotide polymorphism (SNP) assays to infer matrilineal ancestry at a (near) continental level are already available, such tools are lacking for the Oceania region. Here, we have developed a hierarchical system of three SNaPshot multiplexes for genotyping 26 SNPs defining all major mtDNA haplogroups for Oceania (including Australia, Near Oceania and Remote Oceania). With this system, it was possible to conclusively assign 74% of Oceanian individuals to their Oceanian matrilineal ancestry in an established literature database (after correcting for obvious external admixture). Furthermore, in a set of 161 genotyped individuals collected in Australia, Papua New Guinea and Fiji, 87.6% were conclusively assigned an Oceanian matrilineal origin. For the remaining 12.4% of the genotyped samples either a Eurasian origin was detected indicating likely European admixture (1.9%), the identified haplogroups are shared between Oceania and S/SE-Asia (5%), or the SNPs applied did not allow a geographic inference to be assigned (5.6%). Sub-regional assignment within Oceania was possible for 32.9% of the individuals genotyped: 49.5% of Australians were assigned an Australian origin and 13.7% of the Papua New Guineans were assigned a Near Oceanian origin, although none of the Fijians could be assigned a specific Remote Oceanian origin. The low assignment rates of Near and Remote Oceania are explained by recent migrations from Asia via Near Oceania into Remote Oceania. Combining the mtDNA multiplexes for Oceania introduced here with those we developed earlier for all other continental regions, global matrilineal bio-geographic ancestry assignment from DNA is now achievable in a highly efficient way that is also suitable for applications with limited material such as forensic case work.