Deciphering past human population movements in Oceania: provably optimal trees of 127 mtDNA genomes

Aboriginal Australian mitochondrial genome variation - an increased understanding of population antiquity and diversity

Aboriginal Australians represent one of the oldest continuous cultures outside Africa, with evidence indicating that their ancestors arrived in the ancient landmass of Sahul (present-day New Guinea and Australia) ~55 thousand years ago. Genetic studies, though limited, have demonstrated both the uniqueness and antiquity of Aboriginal Australian genomes. We have further resolved known Aboriginal Australian mitochondrial haplogroups and discovered novel indigenous lineages by sequencing the mitogenomes of 127 contemporary Aboriginal Australians. In particular, the more common haplogroups observed in our dataset included M42a, M42c, S, P5 and P12, followed by rarer haplogroups M15, M16, N13, O, P3, P6 and P8. We propose some major phylogenetic rearrangements, such as in haplogroup P where we delinked P4a and P4b and redefined them as P4 (New Guinean) and P11 (Australian), respectively. Haplogroup P2b was identified as a novel clade potentially restricted to Torres Strait Islanders. Nearly all Aboriginal Australian mitochondrial haplogroups detected appear to be ancient, with no evidence of later introgression during the Holocene. Our findings greatly increase knowledge about the geographic distribution and phylogenetic structure of mitochondrial lineages that have survived in contemporary descendants of Australia’s first settlers.

Mitochondrial DNA diversity of present-day Aboriginal Australians and implications for human evolution in Oceania

Aboriginal Australians are one of the more poorly studied populations from the standpoint of human evolution and genetic diversity. Thus, to investigate their genetic diversity, the possible date of their ancestors' arrival and their relationships with neighboring populations, we analyzed mitochondrial DNA (mtDNA) diversity in a large sample of Aboriginal Australians. Selected mtDNA single-nucleotide polymorphisms and the hypervariable segment haplotypes were analyzed in 594 Aboriginal Australians drawn from locations across the continent, chiefly from regions not previously sampled. Most (~78%) samples could be assigned to mtDNA haplogroups indigenous to Australia. The indigenous haplogroups were all ancient (with estimated ages >40 000 years) and geographically widespread across the continent. The most common haplogroup was P (44%) followed by S (23%) and M42a (9%). There was some geographic structure at the haplotype level. The estimated ages of the indigenous haplogroups range from 39 000 to 55 000 years, dates that fit well with the estimated date of colonization of Australia based on archeological evidence (~47 000 years ago). The distribution of mtDNA haplogroups in Australia and New Guinea supports the hypothesis that the ancestors of Aboriginal Australians entered Sahul through at least two entry points. The mtDNA data give no support to the hypothesis of secondary gene flow into Australia during the Holocene, but instead suggest long-term isolation of the continent.

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.

From cheek swabs to consensus sequences: an A to Z protocol for high-throughput DNA sequencing of complete human mitochondrial genomes

BACKGROUND

Next-generation DNA sequencing (NGS) technologies have made huge impacts in many fields of biological research, but especially in evolutionary biology. One area where NGS has shown potential is for high-throughput sequencing of complete mtDNA genomes (of humans and other animals). Despite the increasing use of NGS technologies and a better appreciation of their importance in answering biological questions, there remain significant obstacles to the successful implementation of NGS-based projects, especially for new users.

RESULTS

Here we present an 'A to Z' protocol for obtaining complete human mitochondrial (mtDNA) genomes - from DNA extraction to consensus sequence. Although designed for use on humans, this protocol could also be used to sequence small, organellar genomes from other species, and also nuclear loci. This protocol includes DNA extraction, PCR amplification, fragmentation of PCR products, barcoding of fragments, sequencing using the 454 GS FLX platform, and a complete bioinformatics pipeline (primer removal, reference-based mapping, output of coverage plots and SNP calling).

CONCLUSIONS

All steps in this protocol are designed to be straightforward to implement, especially for researchers who are undertaking next-generation sequencing for the first time. The molecular steps are scalable to large numbers (hundreds) of individuals and all steps post-DNA extraction can be carried out in 96-well plate format. Also, the protocol has been assembled so that individual 'modules' can be swapped out to suit available resources.

The Q2 mitochondrial haplogroup in Oceania

Many details surrounding the origins of the peoples of Oceania remain to be resolved, and as a step towards this we report seven new complete mitochondrial genomes from the Q2a haplogroup, from Papua New Guinea, Fiji and Kiribati. This brings the total to eleven Q2 genomes now available. The Q haplogroup (that includes Q2) is an old and diverse lineage in Near Oceania, and is reasonably common; within our sample set of 430, 97 are of the Q haplogroup. However, only 8 are Q2, and we report 7 here. The tree with all complete Q genomes is proven to be minimal. The dating estimate for the origin of Q2 (around 35 Kya) reinforces the understanding that humans have been in Near Oceania for tens of thousands of years; nevertheless the Polynesian maternal haplogroups remain distinctive. A major focus now, with regard to Polynesian ancestry, is to address the differences and timing of the ‘Melanesian’ contribution to the maternal and paternal lineages as people moved further and further into Remote Oceania. Input from other fields such as anthropology, history and linguistics is required for a better understanding and interpretation of the genetic data.

Complete mitochondrial DNA genome sequences from the first New Zealanders

The dispersal of modern humans across the globe began ~65,000 y ago when people first left Africa and culminated with the settlement of East Polynesia, which occurred in the last 1,000 y. With the arrival of Polynesian canoes only 750 y ago, Aotearoa/New Zealand became the last major landmass to be permanently settled by humans. We present here complete mitochondrial genome sequences of the likely founding population of Aotearoa/New Zealand recovered from the archaeological site of Wairau Bar. These data represent complete mitochondrial genome sequences from ancient Polynesian voyagers and provide insights into the genetic diversity of human populations in the Pacific at the time of the settlement of East Polynesia.

Ancient migration routes of Austronesian-speaking populations in oceanic Southeast Asia and Melanesia might mimic the spread of nasopharyngeal carcinoma

Mitochondrial DNA (mtDNA) and non-recombining Y chromosome (NRY) are inherited uni-parentally from mother to daughter or from father to son respectively. Their polymorphism has initially been studied throughout populations of the world to demonstrate the "Out of Africa" hypothesis. Here, to correlate the distribution of nasopharyngeal carcinoma (NPC) in different populations of insular Asia, we analyze the mtDNA information (lineages) obtained from genotyping of the hyper variable region (HVS I & II) among 1400 individuals from island Southeast Asia (ISEA), Taiwan and Fujian and supplemented with the analysis of relevant coding region polymorphisms. Lineages that best represented a clade (a branch of the genetic tree) in the phylogeny were further analyzed using complete genomic mtDNA sequencing. Finally, these complete mtDNA sequences were used to construct a most parsimonious tree which now constitutes the most up-to-date mtDNA dataset available on ISEA and Taiwan. This analysis has exposed new insights of the evolutionary history of insular Asia and has strong implications in assessing possible correlations with linguistic, archaeology, demography and the NPC distribution in populations within these regions. To obtain a more objective and balanced genetic point of view, slowly evolving biallelic Y single nucleotide polymorphism (Y-SNP) was also analyzed. As in the first step above, the technique was first applied to determine affinities (macro analysis) between populations of insular Asia. Secondly, sixteen Y short tandem repeats (Y-STR) were used as they allow deeper insight (micro analysis) into the relationship between individuals of a same region. Together, mtDNA and NRY allowed a better definition of the relational, demographic, cultural and genetic components that constitute the make up of the present day peoples of ISEA. Outstanding findings were obtained on the routes of migration that occurred along with the spread of NPC during the settlement of insular Asia. The results of this analysis will be discussed using a conceptual approach.

Narrow genetic basis for the Australian dingo confirmed through analysis of paternal ancestry

The dingo (Canis lupus dingo) is an iconic animal in the native culture of Australia, but archaeological and molecular records indicate a relatively recent history on the continent. Studies of mitochondrial DNA (mtDNA) imply that the current dingo population was founded by a small population of already tamed dogs from Southeast Asia. However, the maternal genetic data might give a unilateral picture, and the gene pool has yet to be screened for paternal ancestry. We sequenced 14,437 bp of the Y-chromosome (Y-chr) from two dingoes and one New Guinea Singing Dog (NGSD). This positioned dingo and NGSD within the domestic dog Y-chr phylogeny, and produced one haplotype not detected before. With this data, we characterized 47 male dingoes in 30 Y-chr single-nucleotide polymorphism sites using protease-mediated allele-specific extension technology. Only two haplotypes, H3 and H60, were found among the dingoes, at frequencies of 68.1 and 31.9 %, respectively, compared to 27 haplotypes previously established in the domestic dog. While H3 is common among Southeast Asian dogs, H60 was specifically found in dingoes and the NGSD, but was related to Southeast Asian dog Y-chr haplotypes. H3 and H60 were observed exclusively in the western and eastern parts of Australia, respectively, but had a common range in Southeast. Thus, the Y-chr diversity was very low, similar to previous observations for d-loop mtDNA. Overall genetic evidence suggests a very restricted introduction of the first dingoes into Australia, possibly from New Guinea. This study further confirms the dingo as an isolated feral dog.

Complete mitochondrial genome sequencing reveals novel haplotypes in a Polynesian population

The high risk of metabolic disease traits in Polynesians may be partly explained by elevated prevalence of genetic variants involved in energy metabolism. The genetics of Polynesian populations has been shaped by island hoping migration events which have possibly favoured thrifty genes. The aim of this study was to sequence the mitochondrial genome in a group of Maoris in an effort to characterise genome variation in this Polynesian population for use in future disease association studies. We sequenced the complete mitochondrial genomes of 20 non-admixed Maori subjects using Affymetrix technology. DNA diversity analyses showed the Maori group exhibited reduced mitochondrial genome diversity compared to other worldwide populations, which is consistent with historical bottleneck and founder effects. Global phylogenetic analysis positioned these Maori subjects specifically within mitochondrial haplogroup - B4a1a1. Interestingly, we identified several novel variants that collectively form new and unique Maori motifs – B4a1a1c, B4a1a1a3 and B4a1a1a5. Compared to ancestral populations we observed an increased frequency of non-synonymous coding variants of several mitochondrial genes in the Maori group, which may be a result of positive selection and/or genetic drift effects. In conclusion, this study reports the first complete mitochondrial genome sequence data for a Maori population. Overall, these new data reveal novel mitochondrial genome signatures in this Polynesian population and enhance the phylogenetic picture of maternal ancestry in Oceania. The increased frequency of several mitochondrial coding variants makes them good candidates for future studies aimed at assessment of metabolic disease risk in Polynesian populations.

Complete mitochondrial DNA analysis of eastern Eurasian haplogroups rarely found in populations of northern Asia and eastern Europe

With the aim of uncovering all of the most basal variation in the northern Asian mitochondrial DNA (mtDNA) haplogroups, we have analyzed mtDNA control region and coding region sequence variation in 98 Altaian Kazakhs from southern Siberia and 149 Barghuts from Inner Mongolia, China. Both populations exhibit the prevalence of eastern Eurasian lineages accounting for 91.9% in Barghuts and 60.2% in Altaian Kazakhs. The strong affinity of Altaian Kazakhs and populations of northern and central Asia has been revealed, reflecting both influences of central Asian inhabitants and essential genetic interaction with the Altai region indigenous populations. Statistical analyses data demonstrate a close positioning of all Mongolic-speaking populations (Mongolians, Buryats, Khamnigans, Kalmyks as well as Barghuts studied here) and Turkic-speaking Sojots, thus suggesting their origin from a common maternal ancestral gene pool. In order to achieve a thorough coverage of DNA lineages revealed in the northern Asian matrilineal gene pool, we have completely sequenced the mtDNA of 55 samples representing haplogroups R11b, B4, B5, F2, M9, M10, M11, M13, N9a and R9c1, which were pinpointed from a massive collection (over 5000 individuals) of northern and eastern Asian, as well as European control region mtDNA sequences. Applying the newly updated mtDNA tree to the previously reported northern Asian and eastern Asian mtDNA data sets has resolved the status of the poorly classified mtDNA types and allowed us to obtain the coalescence age estimates of the nodes of interest using different calibrated rates. Our findings confirm our previous conclusion that northern Asian maternal gene pool consists of predominantly post-LGM components of eastern Asian ancestry, though some genetic lineages may have a pre-LGM/LGM origin.

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.

Ancient voyaging and Polynesian origins

The "Polynesian motif" defines a lineage of human mtDNA that is restricted to Austronesian-speaking populations and is almost fixed in Polynesians. It is widely thought to support a rapid dispersal of maternal lineages from Taiwan ~4000 years ago (4 ka), but the chronological resolution of existing control-region data is poor, and an East Indonesian origin has also been proposed. By analyzing 157 complete mtDNA genomes, we show that the motif itself most likely originated >6 ka in the vicinity of the Bismarck Archipelago, and its immediate ancestor is >8 ka old and virtually restricted to Near Oceania. This indicates that Polynesian maternal lineages from Island Southeast Asia gained a foothold in Near Oceania much earlier than dispersal from either Taiwan or Indonesia 3-4 ka would predict. However, we find evidence in minor lineages for more recent two-way maternal gene flow between Island Southeast Asia and Near Oceania, likely reflecting movements along a "voyaging corridor" between them, as previously proposed on archaeological grounds. Small-scale mid-Holocene movements from Island Southeast Asia likely transmitted Austronesian languages to the long-established Southeast Asian colonies in the Bismarcks carrying the Polynesian motif, perhaps also providing the impetus for the expansion into Polynesia.

High-throughput sequencing of complete human mtDNA genomes from the Philippines

Because of the time and cost associated with Sanger sequencing of complete human mtDNA genomes, practically all evolutionary studies have screened samples first to define haplogroups and then either selected a few samples from each haplogroup, or many samples from a particular haplogroup of interest, for complete mtDNA genome sequencing. Such biased sampling precludes many analyses of interest. Here, we used high-throughput sequencing platforms to generate, rapidly and inexpensively, 109 complete mtDNA genome sequences from random samples of individuals from three Filipino groups, including one Negrito group, the Mamanwa. We obtained on average ∼55-fold coverage per sequence, with <1% missing data per sequence. Various analyses attest to the accuracy of the sequences, including comparison to sequences of the first hypervariable segment of the control region generated by Sanger sequencing; patterns of nucleotide substitution and the distribution of polymorphic sites across the genome; and the observed haplogroups. Bayesian skyline plots of population size change through time indicate similar patterns for all three Filipino groups, but sharply contrast with such plots previously constructed from biased sampling of complete mtDNA genomes, as well as with an artificially constructed sample of sequences that mimics the biased sampling. Our results clearly demonstrate that the high-throughput sequencing platforms are the methodology of choice for generating complete mtDNA genome sequences.

European admixture on the Micronesian island of Kosrae: lessons from complete genetic information

The architecture of natural variation present in a contemporary population is a result of multiple population genetic forces, including population bottleneck and expansion, selection, drift, and admixture. We seek to untangle the contribution of admixture to genetic diversity on the Micronesian island of Kosrae. Toward this goal, we used a complete genetic approach by combining a dense genome-wide map of 100,000 single-nucleotide polymorphisms (SNPs) with data from uniparental markers from the mitochondrial genome and the nonrecombining portion of the Y chromosome. These markers were typed in approximately 3200 individuals from Kosrae, representing 80% of the adult population of the island. We developed novel software that uses SNP data to delineate ancestry for individual segments of the genome. Through this analysis, we determined that 39% of Kosraens have some European ancestry. However, the vast majority of admixed individuals (77%) have European alleles spanning less than 10% of their genomes. Data from uniparental markers show most of this admixture to be male, introduced in the late nineteenth century. Furthermore, pedigree analysis shows that the majority of European admixture on Kosrae is because of the contribution of one individual. This approach shows the benefit of combining information from autosomal and uniparental polymorphisms and provides new methodology for determining ancestry in a population.

Complete mitochondrial DNA sequences provide new insights into the Polynesian motif and the peopling of Madagascar

More than a decade of mitochondrial DNA (mtDNA) studies have given the 'Polynesian motif' renowned status as a marker for tracing the late-Holocene expansion of Austronesian speaking populations. Despite considerable research on the Polynesian motif in Oceania, there has been little equivalent work on the western edge of its expansion - leaving major issues unresolved regarding the motif's evolutionary history. This has also led to considerable uncertainty regarding the settlement of Madagascar. In this study, we assess mtDNA variation in 266 individuals from three Malagasy ethnic groups: the Mikea, Vezo, and Merina. Complete mtDNA genome sequencing reveals a new variant of the Polynesian motif in Madagascar; two coding region mutations define a Malagasy-specific sub-branch. This newly defined 'Malagasy motif' occurs at high frequency in all three ethnic groups (13-50%), and its phylogenetic position, geographic distribution, and estimated age all support a recent origin, but without conclusively identifying a specific source region. Nevertheless, the haplotype's limited diversity, similar to those of other mtDNA haplogroups found in our Malagasy groups, best supports a small number of initial settlers arriving to Madagascar through the same migratory process. Finally, the discovery of this lineage provides a set of new polymorphic positions to help localize the Austronesian ancestors of the Malagasy, as well as uncover the origin and evolution of the Polynesian motif itself.

Correcting for purifying selection: an improved human mitochondrial molecular clock

There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at approximately 15 kya that-unlike the uncorrected clock-matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55-70 kya, 5-20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses.

Migration of Chadic speaking pastoralists within Africa based on population structure of Chad Basin and phylogeography of mitochondrial L3f haplogroup

Background

Chad Basin, lying within the bidirectional corridor of African Sahel, is one of the most populated places in Sub-Saharan Africa today. The origin of its settlement appears connected with Holocene climatic ameliorations (aquatic resources) that started ~10,000 years before present (YBP). Although both Nilo-Saharan and Niger-Congo language families are encountered here, the most diversified group is the Chadic branch belonging to the Afro-Asiatic language phylum. In this article, we investigate the proposed ancient migration of Chadic pastoralists from Eastern Africa based on linguistic data and test for genetic traces of this migration in extant Chadic speaking populations.

Results

We performed whole mitochondrial genome sequencing of 16 L3f haplotypes, focused on clade L3f3 that occurs almost exclusively in Chadic speaking people living in the Chad Basin. These data supported the reconstruction of a L3f phylogenetic tree and calculation of times to the most recent common ancestor for all internal clades. A date ~8,000 YBP was estimated for the L3f3 sub-haplogroup, which is in good agreement with the supposed migration of Chadic speaking pastoralists and their linguistic differentiation from other Afro-Asiatic groups of East Africa. As a whole, the Afro-Asiatic language family presents low population structure, as 92.4% of mtDNA variation is found within populations and only 3.4% of variation can be attributed to diversity among language branches. The Chadic speaking populations form a relatively homogenous cluster, exhibiting lower diversification than the other Afro-Asiatic branches (Berber, Semitic and Cushitic).

Conclusion

The results of our study support an East African origin of mitochondrial L3f3 clade that is present almost exclusively within Chadic speaking people living in Chad Basin. Whole genome sequence-based dates show that the ancestral haplogroup L3f must have emerged soon after the Out-of-Africa migration (around 57,100 ± 9,400 YBP), but the "Chadic" L3f3 clade has much less internal variation, suggesting an expansion during the Holocene period about 8,000 ± 2,500 YBP. This time period in the Chad Basin is known to have been particularly favourable for the expansion of pastoralists coming from northeastern Africa, as suggested by archaeological, linguistic and climatic data.

Revealing the hidden complexities of mtDNA inheritance

Mitochondrial DNA (mtDNA) is a pivotal tool in molecular ecology, evolutionary and population genetics. The power of mtDNA analyses derives from a relatively high mutation rate and the apparent simplicity of mitochondrial inheritance (maternal, without recombination), which has simplified modelling population history compared to the analysis of nuclear DNA. However, in biology things are seldom simple, and advances in DNA sequencing and polymorphism detection technology have documented a growing list of exceptions to the central tenets of mitochondrial inheritance, with paternal leakage, heteroplasmy and recombination now all documented in multiple systems. The presence of paternal leakage, recombination and heteroplasmy can have substantial impact on analyses based on mtDNA, affecting phylogenetic and population genetic analyses, estimates of the coalescent and the myriad of other parameters that are dependent on such estimates. Here, we review our understanding of mtDNA inheritance, discuss how recent findings mean that established ideas may need to be re-evaluated, and we assess the implications of these new-found complications for molecular ecologists who have relied for decades on the assumption of a simpler mode of inheritance. We show how it is possible to account for recombination and heteroplasmy in evolutionary and population analyses, but that accurate estimates of the frequencies of biparental inheritance and recombination are needed. We also suggest how nonclonal inheritance of mtDNA could be exploited, to increase the ways in which mtDNA can be used in analyses.

mtDNAmanager: a Web-based tool for the management and quality analysis of mitochondrial DNA control-region sequences

Background

For the past few years, scientific controversy has surrounded the large number of errors in forensic and literature mitochondrial DNA (mtDNA) data. However, recent research has shown that using mtDNA phylogeny and referring to known mtDNA haplotypes can be useful for checking the quality of sequence data.

Results

We developed a Web-based bioinformatics resource "mtDNAmanager" that offers a convenient interface supporting the management and quality analysis of mtDNA sequence data. The mtDNAmanager performs computations on mtDNA control-region sequences to estimate the most-probable mtDNA haplogroups and retrieves similar sequences from a selected database. By the phased designation of the most-probable haplogroups (both expected and estimated haplogroups), mtDNAmanager enables users to systematically detect errors whilst allowing for confirmation of the presence of clear key diagnostic mutations and accompanying mutations. The query tools of mtDNAmanager also facilitate database screening with two options of "match" and "include the queried nucleotide polymorphism". In addition, mtDNAmanager provides Web interfaces for users to manage and analyse their own data in batch mode.

Conclusion

The mtDNAmanager will provide systematic routines for mtDNA sequence data management and analysis via easily accessible Web interfaces, and thus should be very useful for population, medical and forensic studies that employ mtDNA analysis. mtDNAmanager can be accessed at .

Characterizing the time dependency of human mitochondrial DNA mutation rate estimates

Previous research has established a discrepancy of nearly an order of magnitude between pedigree-based and phylogeny-based (human vs. chimpanzee) estimates of the mitochondrial DNA (mtDNA) control region mutation rate. We characterize the time dependency of the human mitochondrial hypervariable region one mutation rate by generating 14 new phylogeny-based mutation rate estimates using within-human comparisons and archaeological dates. Rate estimates based on population events between 15,000 and 50,000 years ago are at least 2-fold lower than pedigree-based estimates. These within-human estimates are also higher than estimates generated from phylogeny-based human-chimpanzee comparisons. Our new estimates establish a rapid decay in evolutionary mutation rate between approximately 2,500 and 50,000 years ago and a slow decay from 50,000 to 6 Ma. We then extend this analysis to the mtDNA-coding region. Our within-human coding region mutation rate estimates display a similar, though less rapid, time-dependent decay. We explore the possibility that multiple hits explain the discrepancy between pedigree-based and phylogeny-based mutation rates. We conclude that whereas nucleotide substitution models incorporating multiple hits do provide a possible explanation for the discrepancy between pedigree-based and human-chimpanzee mutation rate estimates, they do not explain the rapid decline of within-human rate estimates. We propose that demographic processes such as serial bottlenecks prior to the Holocene could explain the difference between rates estimated before and after 15,000 years ago. Our findings suggest that human mtDNA estimates of dates of population and phylogenetic events should be adjusted in light of this time dependency of the mutation rate estimates.

Mitochondrial DNA variation in Karkar Islanders

We analyzed 375 base pairs (bp) of the first hypervariable region (HVS-I) of the mitochondrial DNA (mtDNA) control region and intergenic COII/tRNALys 9-bp deletion from 47 Karkar Islanders (north coast of Papua New Guinea) belonging to the Waskia Papuan language group. To address questions concerning the origin and evolution of this population we compared the Karkar mtDNA haplotypes and haplogroups to those of neighbouring East Asians and Oceanic populations. The results of the phylogeographic analysis show grouping in three different clusters of the Karkar Islander mtDNA lineages: one group of lineages derives from the first Pleistocene settlers of New Guinea-Island Melanesia, a second set derives from more recent arrivals of Austronesian speaking populations, and the third contains lineages specific to the Karkar Islanders, but still rooted to Austronesian and New Guinea-Island Melanesia populations. Our results suggest (i) the absence of a strong association between language and mtDNA variation and, (ii) reveal that the mtDNA haplogroups F1a1, M7b1 and E1a, which probably originated in Island Southeast Asia and may be considered signatures of Austronesian population movements, are preserved in the Karkar Islanders but absent in other New Guinea-Island Melanesian populations. These findings indicate that the Karkar Papuan speakers retained a certain degree of their own genetic uniqueness and a high genetic diversity. We present a hypothesis based on archaeological, linguistic and environmental datasets to argue for a succession of (partial) depopulation and repopulation and expansion events, under conditions of structured interaction, which may explain the variability expressed in the Karkar mtDNA.

The genetic structure of Pacific Islanders

Human genetic diversity in the Pacific has not been adequately sampled, particularly in Melanesia. As a result, population relationships there have been open to debate. A genome scan of autosomal markers (687 microsatellites and 203 insertions/deletions) on 952 individuals from 41 Pacific populations now provides the basis for understanding the remarkable nature of Melanesian variation, and for a more accurate comparison of these Pacific populations with previously studied groups from other regions. It also shows how textured human population variation can be in particular circumstances. Genetic diversity within individual Pacific populations is shown to be very low, while differentiation among Melanesian groups is high. Melanesian differentiation varies not only between islands, but also by island size and topographical complexity. The greatest distinctions are among the isolated groups in large island interiors, which are also the most internally homogeneous. The pattern loosely tracks language distinctions. Papuan-speaking groups are the most differentiated, and Austronesian or Oceanic-speaking groups, which tend to live along the coastlines, are more intermixed. A small “Austronesian” genetic signature (always <20%) was detected in less than half the Melanesian groups that speak Austronesian languages, and is entirely lacking in Papuan-speaking groups. Although the Polynesians are also distinctive, they tend to cluster with Micronesians, Taiwan Aborigines, and East Asians, and not Melanesians. These findings contribute to a resolution to the debates over Polynesian origins and their past interactions with Melanesians. With regard to genetics, the earlier studies had heavily relied on the evidence from single locus mitochondrial DNA or Y chromosome variation. Neither of these provided an unequivocal signal of phylogenetic relations or population intermixture proportions in the Pacific. Our analysis indicates the ancestors of Polynesians moved through Melanesia relatively rapidly and only intermixed to a very modest degree with the indigenous populations there.

The origins and current genetic relationships of Pacific Islanders have been the subjects of interest and controversy for many decades. By analyzing the variation of a large number (687) of genetic markers in almost 1,000 individuals from 41 Pacific populations, and comparing these with East Asians and others, we contribute to the clarification and resolution of many of these issues. To judge by the populations in our survey, we find that Polynesians and Micronesians have almost no genetic relation to Melanesians, but instead are strongly related to East Asians, and particularly Taiwan Aborigines. A minority of Island Melanesian populations have indications of a small shared genetic ancestry with Polynesians and Micronesians (the ones that have this tie all speak related Austronesian languages). Inland groups who speak Papuan languages are particularly divergent and internally homogeneous. The genetic divergence among Island Melanesian populations, which is neatly organized by island, island size/topography, as well as their coastal or inland locations, is remarkable for such a small region, and enlarges our understanding of the texture of contemporary human variation.

Prehistoric inter-archipelago trading of Polynesian tree snails leaves a conservation legacy

Inter-archipelago exchange networks were an important aspect of prehistoric Polynesian societies. We report here a novel genetic characterization of a prehistoric exchange network involving an endemic Pacific island tree snail, Partula hyalina. It occurs in the Society (Tahiti only), Austral and Southern Cook Islands. Our genetic data, based on museum, captive and wild-caught samples, establish Tahiti as the source island. The source lineage is polymorphic in shell coloration and contains a second nominal species, the dark-shelled Partula clara, in addition to the white-shelled P. hyalina. Prehistoric inter-island introductions were non-random: they involved white-shelled snails only and were exclusively inter-archipelago in scope. Partulid shells were commonly used in regional Polynesian jewellery, and we propose that the white-shelled P. hyalina, originally restricted to Tahiti, had aesthetic value throughout these archipelagoes. Demand within the Society Islands could be best met by trading dead shells, but a low rate of inter-archipelago exchange may have prompted the establishment of multiple founder populations in the Australs and Southern Cooks. The alien carnivorous land snail Euglandina rosea has recently devastated populations of all 61 endemic species of Society Island partulid snails. Southern Cooks and Australs P. hyalina now represent the only unscathed wild populations remaining of this once spectacular land snail radiation.

Genome-wide analysis indicates more Asian than Melanesian ancestry of Polynesians

Analyses of mitochondrial DNA (mtDNA) and nonrecombining Y chromosome (NRY) variation in the same populations are sometimes concordant but sometimes discordant. Perhaps the most dramatic example known of the latter concerns Polynesians, in which about 94% of Polynesian mtDNAs are of East Asian origin, while about 66% of Polynesian Y chromosomes are of Melanesian origin. Here we analyze on a genome-wide scale, to our knowledge for the first time, the origins of the autosomal gene pool of Polynesians by screening 377 autosomal short tandem repeat (STR) loci in 47 Pacific Islanders and compare the results with those obtained from 44 Chinese and 24 individuals from Papua New Guinea. Our data indicate that on average about 79% of the Polynesian autosomal gene pool is of East Asian origin and 21% is of Melanesian origin. The genetic data thus suggest a dual origin of Polynesians with a high East Asian but also considerable Melanesian component, reflecting sex-biased admixture in Polynesian history in agreement with the Slow Boat model. More generally, these results also demonstrate that conclusions based solely on uniparental markers, which are frequently used in population history studies, may not accurately reflect the history of the autosomal gene pool of a population.

Revealing the prehistoric settlement of Australia by Y chromosome and mtDNA analysis

Published and new samples of Aboriginal Australians and Melanesians were analyzed for mtDNA (n=172) and Y variation (n=522), and the resulting profiles were compared with the branches known so far within the global mtDNA and the Y chromosome tree. (i) All Australian lineages are confirmed to fall within the mitochondrial founder branches M and N and the Y chromosomal founders C and F, which are associated with the exodus of modern humans from Africa approximately 50-70,000 years ago. The analysis reveals no evidence for any archaic maternal or paternal lineages in Australians, despite some suggestively robust features in the Australian fossil record, thus weakening the argument for continuity with any earlier Homo erectus populations in Southeast Asia. (ii) The tree of complete mtDNA sequences shows that Aboriginal Australians are most closely related to the autochthonous populations of New Guinea/Melanesia, indicating that prehistoric Australia and New Guinea were occupied initially by one and the same Palaeolithic colonization event approximately 50,000 years ago, in agreement with current archaeological evidence. (iii) The deep mtDNA and Y chromosomal branching patterns between Australia and most other populations around the Indian Ocean point to a considerable isolation after the initial arrival. (iv) We detect only minor secondary gene flow into Australia, and this could have taken place before the land bridge between Australia and New Guinea was submerged approximately 8,000 years ago, thus calling into question that certain significant developments in later Australian prehistory (the emergence of a backed-blade lithic industry, and the linguistic dichotomy) were externally motivated.

Mitochondrial diversity within modern human populations

With the recent increase in the available number of high-quality, full-length mitochondrial sequences, it is now possible to construct and analyze a comprehensive human mitochondrial consensus sequence. Using a data set of 827 carefully selected sequences, it is shown that modern humans contain extremely low levels of divergence from the mitochondrial consensus sequence, differing by a mere 21.6 nt sites on average. Fully 84.1% of the mitochondrial genome was found to be invariant and ‘private’ mutations accounted for 43.8% of the variable sites. Ninety eight percent of the variant sites had a primary nucleotide with an allele frequency of 0.90 or greater. Interestingly, the few truly ambiguous nucleotide sites could all be reliably assigned to either a purine or pyrimidine ancestral state. A comparison of this consensus sequence to several ancestral sequences derived from phylogenetic studies reveals a great deal of similarity, where, as expected, the most phylogenetically informative nucleotides in the ancestral studies tended to be the most variable nucleotides in the consensus. Allowing for this fact, the consensus approach provides variation data on the positions that do not contribute to phylogenetic reconstructions, and these data provide a baseline for measuring human mitochondrial variation in populations worldwide.

Melanesian mtDNA complexity

Melanesian populations are known for their diversity, but it has been hard to grasp the pattern of the variation or its underlying dynamic. Using 1,223 mitochondrial DNA (mtDNA) sequences from hypervariable regions 1 and 2 (HVR1 and HVR2) from 32 populations, we found the among-group variation is structured by island, island size, and also by language affiliation. The more isolated inland Papuan-speaking groups on the largest islands have the greatest distinctions, while shore dwelling populations are considerably less diverse (at the same time, within-group haplotype diversity is less in the most isolated groups). Persistent differences between shore and inland groups in effective population sizes and marital migration rates probably cause these differences. We also add 16 whole sequences to the Melanesian mtDNA phylogenies. We identify the likely origins of a number of the haplogroups and ancient branches in specific islands, point to some ancient mtDNA connections between Near Oceania and Australia, and show additional Holocene connections between Island Southeast Asia/Taiwan and Island Melanesia with branches of haplogroup E. Coalescence estimates based on synonymous transitions in the coding region suggest an initial settlement and expansion in the region at approximately 30-50,000 years before present (YBP), and a second important expansion from Island Southeast Asia/Taiwan during the interval approximately 3,500-8,000 YBP. However, there are some important variance components in molecular dating that have been overlooked, and the specific nature of ancestral (maternal) Austronesian influence in this region remains unresolved.

A mitochondrial stratigraphy for island southeast Asia

Island Southeast Asia (ISEA) was first colonized by modern humans at least 45,000 years ago, but the extent to which the modern inhabitants trace their ancestry to the first settlers is a matter of debate. It is widely held, in both archaeology and linguistics, that they are largely descended from a second wave of dispersal, proto-Austronesian-speaking agriculturalists who originated in China and spread to Taiwan approximately 5,500 years ago. From there, they are thought to have dispersed into ISEA approximately 4,000 years ago, assimilating the indigenous populations. Here, we demonstrate that mitochondrial DNA diversity in the region is extremely high and includes a large number of indigenous clades. Only a fraction of these date back to the time of first settlement, and the majority appear to mark dispersals in the late-Pleistocene or early-Holocene epoch most likely triggered by postglacial flooding. There are much closer genetic links to Taiwan than to the mainland, but most of these probably predated the mid-Holocene "Out of Taiwan" event as traditionally envisioned. Only approximately 20% at most of modern mitochondrial DNAs in ISEA could be linked to such an event, suggesting that, if an agriculturalist migration did take place, it was demographically minor, at least with regard to the involvement of women.