MBE Advance Access originally published online on December 18, 2007
Molecular Biology and Evolution 2008 25(2):468-474; doi:10.1093/molbev/msm277
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
mtDNA Variation Predicts Population Size in Humans and Reveals a Major Southern Asian Chapter in Human Prehistory
* Institute of Cognitive and Evolutionary Anthropology, University of Oxford, 58A Banbury Rd, Oxford, OX2 6QS, United Kingdom
Psychology Department, University of Auckland, Auckland, New Zealand
Bioinformatics Institute and Department of Computer Science, University of Auckland, Auckland, New Zealand
E-mail: quentin.atkinson{at}anthro.ox.ac.uk.
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
The relative timing and size of regional human population growth following our expansion from Africa remain unknown. Human mitochondrial DNA (mtDNA) diversity carries a legacy of our population history. Given a set of sequences, we can use coalescent theory to estimate past population size through time and draw inferences about human population history. However, recent work has challenged the validity of using mtDNA diversity to infer species population sizes. Here we use Bayesian coalescent inference methods, together with a global data set of 357 human mtDNA coding-region sequences, to infer human population sizes through time across 8 major geographic regions. Our estimates of relative population sizes show remarkable concordance with the contemporary regional distribution of humans across Africa, Eurasia, and the Americas, indicating that mtDNA diversity is a good predictor of population size in humans. Plots of population size through time show slow growth in sub-Saharan Africa beginning 143–193 kya, followed by a rapid expansion into Eurasia after the emergence of the first non-African mtDNA lineages 50–70 kya. Outside Africa, the earliest and fastest growth is inferred in Southern Asia
52 kya, followed by a succession of growth phases in Northern and Central Asia (49 kya), Australia (48 kya), Europe (42 kya), the Middle East and North Africa (40 kya), New Guinea (39 kya), the Americas (18 kya), and a second expansion in Europe (10–15 kya). Comparisons of relative regional population sizes through time suggest that between approximately 45 and 20 kya most of humanity lived in Southern Asia. These findings not only support the use of mtDNA data for estimating human population size but also provide a unique picture of human prehistory and demonstrate the importance of Southern Asia to our recent evolutionary past.
Key Words: mitochondrial DNA • human evolution • population genetics • Bayesian skyline plot, anthropology
|
Introduction |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
How humans colonized the globe remains "one of the greatest untold stories in the history of humankind" (Goebel 2007
). Mitochondrial DNA (mtDNA) has been a crucial line of evidence in developing the current understanding of our genetic prehistory. Phylogenetic studies of human mtDNA variation support a late Pleistocene expansion of modern humans from Africa (Cann et al. 1987; Vigilant et al. 1991; Watson et al. 1997; Ingman et al. 2000). More recent work suggests a "southern" migration route from sub-Saharan Africa along the Indian Ocean coast into Eurasia (Forster and Matsumura 2005; Macaulay et al. 2005; Thangaraj et al. 2005) and Sahul (Ingman and Gyllensten 2003; Friedlaender et al. 2005; van Holst Pellekaan et al. 2006) and a later migration from the Levant into North Africa and Europe (Olivieri et al. 2006). However, despite these advances, we lack a clear understanding of the timing and size of regional human population expansions.
Past population size can be estimated from genetic data using coalescent theory (Hudson 1990). For any locus, the information that can be extracted about population size depends on the accuracy with which the underlying genealogical tree can be resolved and the extent to which the historical demographic processes can be accommodated by the coalescent model. mtDNA's high copy number, absence of recombination (Elson et al. 2001), and rapid substitution rate (Howell et al. 1996), mean that its genealogical tree is well resolved, thus making mtDNA very well suited for coalescent inference of population size. Recently, however, Bazin et al. (2006) have argued that mtDNA variation is a poor indicator of population size in animals, probably due to the effects of selection acting to reduce diversity (known as genetic draft). Bazin et al. compared intraspecies mtDNA and allozyme diversity across major taxonomic divisions and found no association between mtDNA diversity and population size, whereas allozyme diversity did predict population size. Although the pattern Bazin et al. (2006) identify appears not to hold for human genetic diversity (human mtDNA and autosomal diversity both imply a similar population size [Eyre-Walker 2006]) human mtDNA does show evidence of selection at some loci (Mishmar et al. 2003; Kivisild et al. 2006). This raises the question of whether mtDNA is in fact a reliable predictor of human population size.
Previous work has used coalescent analyses of pairwise mtDNA sequence differences to argue for exponential growth in human population size beginning 30–130 kya (Harpending et al. 1993; Sherry et al. 1994; Rogers 1995; Ingman et al. 2000). Recent improvements in the available inference methods and complete mtDNA sequence data mean that population genetic parameters can now be estimated from sequence data with greater accuracy, without requiring a simple parametric growth model, and with an explicit framework for quantifying uncertainty in parameter estimates. The Bayesian skyline plot (BSP) (Drummond and Rambaut 2003; Drummond et al. 2005) uses Bayesian coalescent inference of phylogeny and a Markov Chain Monte Carlo (MCMC) (Metropolis et al. 1953) sampling algorithm to simultaneously estimate a posterior probability distribution for the ancestral genealogy, branch lengths, substitution model parameters, and population parameters through time, given a set of gene sequences. The resulting BSP represents a credibility interval for effective population size that incorporates uncertainty in nuisance parameters such as the underlying ancestral gene tree, branch lengths, and rate parameters. Unlike previous coalescent inference methods, by generating a sample of piecewise population size estimates, results are not contingent on a prespecified parametric growth model and the posterior distribution can be used to directly assess uncertainty in parameter estimates given the assumptions of the model. In addition, information loss is minimized by inferring population parameters directly from the sequence data using a coalescent framework, rather than converting sequence data to pairwise distances between sequences (Steel et al. 1988).
In this report, we use the BSP together with a global data set of 357 human mtDNA–coding region sequences to estimate effective population sizes through time for 8 major geographic regions. By comparing our coalescent-based estimates of regional effective population size with anthropological estimates of census population sizes, we show that mtDNA diversity is a remarkably good predictor of regional human population size. As a result, we are able to use our regional population size estimates to gain novel insight into the human colonization of the globe during the last 100,000 years. Most strikingly, we find evidence of a major Southern Asian chapter in our recent evolutionary past.
|
Materials and Methods |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Sequence Data
A data set of 357 human mtDNA–coding region sequences were assembled from the >3,000 sequences available on GenBank and divided into 8 geographic regions (supplementary table S1, Supplementary Material online). Only sequences that could be unambiguously assigned to regional groupings were considered for analysis. We sought a balanced sample of between approximately 25 and 50 sequences from each of the regions, with a larger sampling of the more diverse African population. Data were machine aligned in ClustalX (Thompson et al. 1997
) using reference sequence J01415.1. Regional boundaries (fig. 2a) were defined based on geographical barriers and independent genetic evidence for population subdivisions (Cavalli-Sforza et al. 1994). Results were found to be robust across a range of different regional definitions (supplementary figs. 1–3, Supplementary Material online). Alternative regional sequence groupings are shown in supplementary table S2 (Supplementary Material online).
|
Bayesian Skyline Plots
BSPs (Drummond et al. 2005
) of effective population size were produced from the data using MCMC (Metropolis et al. 1953) sampling in the program BEAST (version 1.3) (Drummond and Rambaut 2003). Effective population size is a compound population genetic parameter generally considered linearly proportional to census population size—here, the population of breeding females. It is influenced by many factors, including local extinction and recolonization and various forms of nonrandom mating (Wakeley 2000), and assumes that regional populations are isolated (i.e., there is no migration between populations). The underlying population size function of the BSP can be specified either as a stepwise or as a piecewise linear function of population size change. The main figures were obtained with a piecewise linear model. Estimates of effective population size are derived from the inferred lineage coalescent rate through time. Ancestral gene trees were inferred based on a general time-reversible (GTR) substitution model with site-specific rates for first, second, and third codon positions. A Bayes factor computed via importance sampling (Newton et al. 1994) indicated that this model fitted the data better than the commonly used GTR + 
+ I model (log Bayes factor = 19.1). Each MCMC sample was based on a run of 40,000,000 generations, sampled every 4,000, with the first 4,000,000 generations discarded as burn-in. Examination of autocorrelation times of the MCMC plots indicated runs had converged to the equilibrium distribution and provided adequate posterior samples by this time. All runs had an effective sample size of the coalescent prior of at least 1,000.
In order to plot population size with respect to time, it was necessary to calibrate rates of molecular evolution. Recent evidence of systematic variation in observed rates of molecular evolution through time (Ho et al. 2005) suggests that rates should be calibrated using divergence events of a similar time depth to the timescale of interest. For this reason, an internal rate calibration was used. Rates were fixed to 1.691 x 10–8 substitutions/site/year, based on a 45-kya age for haplotype Q in New Guinea. Haplogroup Q occurs at high frequencies only in New Guinea (parts of Melanesia have a closely related subtype [Q2], which also occurs at very low frequencies in Micronesia and Polynesia, probably as a result of more recent admixture [Friedlaender et al. 2005]). This is consistent with archaeological evidence for modern human arrival on the super continent Sahul (comprising what is now New Guinea and Australia) by approximately 42–45 kya (O'Connell and Allen 2004). In order to isolate uncertainty in the coalescent process and facilitate regional comparison, uncertainty in the calibration time itself was not incorporated into our analyses. Inferred times do, however, scale linearly with respect to the rate calibration, such that if, for example, a 50-kya age for the New Guinea clade were assumed, the timescale would be increased by approximately 11%.
|
Results and Discussion |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Figure 1 shows BSPs of effective population size through time for each of the 8 regions examined here. The growth curves are derived from changes in the inferred lineage coalescent rate through time. These changes are generally taken to reflect changes in census population size but could also be influenced by changes in population structure, selection pressures, or some combination of these. We can test whether our effective population size estimates provide a good indication of relative census population sizes by comparing relative size estimates for the "present" with contemporary regional population sizes. Figure 2b shows that relative regional population size estimates at time t = 0 are remarkably consistent with independent anthropological and historical estimates of relative precolonial population sizes across Africa, Eurasia, and the Americas (Biraben 1979
). Only estimates for Australia and New Guinea differ significantly from the anthropological estimates. This is likely to be due to differences in population level processes between Australia and New Guinea and the other 6 regions. Genetic distance between populations in Australia and New Guinea rises quickly with increasing geographic distance (Cavalli-Sforza et al. 1994). Studies of Australian Aboriginals show that tribes are genetically very isolated (Huoponen et al. 2001), with the majority of marriages within the tribe (Birdsell 1993). Highly structured populations with low-migration rates and little extinction/recolonization can be expected to preserve variation that would be lost due to drift in a larger metapopulation and hence produce relatively high effective population size estimates (Wakeley 2000). Excluding Australia and New Guinea, there is a striking relationship between effective and census population size estimates (r = 0.98). This relationship was found to be robust across different definitions of the regional group boundaries (supplementary figs. 1–3, Supplementary Material online). These results indicate that mtDNA variation can be used to infer relative regional population sizes in humans.
|
The finding of Bazin et al. (2006)
that population size does not influence species’ mtDNA diversity across the animal kingdom appears not to apply, at least at the population level, in humans. This is most likely because the historically small human population size and short timescales involved make it unlikely that a selective sweep has been able to "reset" regional mtDNA diversity since our expansion from Africa. We were therefore able to combine results from figure 1 to estimate the global population distribution of modern humans over the last 100,000 years (fig. 2c). An important caveat in interpreting these figures is that the population demographic history inferred from sequences sampled from a particular region does not necessarily represent the population history of the region but rather the history of the genetic lineages now inhabiting the region. For example, regional interference can occur if a major migration event brings the genetic signature from a region with a long history to a newly inhabited region. However, we note that even in cases where we know this has happened, such as the expansion from Africa itself and the colonization of the Americas, our demographic estimates appear to correctly infer a recent expansion event rather than being mislead by the history of the parent population.
Figures 1 and 2 together provide a unique picture of the human colonization of the globe. The plots all show evidence of significant population growth and reveal clear differences in regional population histories. The timescale in these figures was calibrated by assuming an mtDNA evolutionary rate of 1.691 x 10–8 substitutions/site/year. Using a slower rate would lead to older dates than presented. However, the relative times of the population expansions across regions would not be altered by a different calibration. We infer slow, roughly exponential population growth in sub-Saharan Africa from an most recent common ancestor (MRCA) 143–193 kya (95% highest posterior density [HPD]; fig. 1a), followed by a period of rapid growth in Eurasia shortly after the emergence of the first non-African mtDNA lineages (N-haplogroup MRCA 53–69 kya, 95% HPD; M haplogroup MRCA 50–64 kya, 95% HPD). This result is consistent with analyses of human microsatellite and SNPs data suggesting slow growth in Africa and more recent rapid growth following a bottleneck in Asia and Europe (Gonser et al. 2000; Zhivotovsky et al. 2000; Alonso and Armour 2001; Rogers 2001; Marth et al. 2004). The congruence between our findings and those from unlinked genetic loci provides further evidence that the observed trends in effective population size are not merely an artifact of selection, which would require that identical selection pressures had acted simultaneously on unrelated genes.
Southern Asia shows the earliest and most pronounced population expansion outside Africa, with a 5-fold increase in population size estimated by 52 kya (fig. 1b). An early Southern Asian growth phase is consistent with high mtDNA diversity in populations from the Indian Ocean Coast (Kivisild et al. 2003; Macaulay et al. 2005; Thangaraj et al. 2005; Sun et al. 2006), as well as recent archaeological evidence (Mellars 2006a), and strongly supports a rapid coastal migration along the "Southern Route" from Africa into Southern Asia. A southern expansion route is also supported by the analyses in supplementary figure S1 (Supplementary Material online), with growth inferred first in Southwestern Asia (India/Middle East) and slightly later in South East Asia. Perhaps most striking, however, is the magnitude of the Southern Asian growth phase, which implies that between approximately 45 and 20 kya over half of the global human population lived on the Indian subcontinent and what is now the Thai and Malay Peninsulas. Our population size estimates indicate that the proportion of people living in the region peaked at over 60% approximately 38 kya. This figure may even be an underestimate given the potential for regional interference along the lines mentioned above because early Southern Asian growth could cause us to infer a premature increase in the population size estimates for the other regions.
We infer less rapid, progressively later population growth in Northern and Central Asia (49 kya; fig. 1c), Europe (42 kya; fig. 1e), and the Middle East and North Africa (40 kya; fig. 1f). Slower growth rates in these regions may reflect less favorable conditions for human occupation at this time. The later timing of growth in Europe and the Middle East and North Africa is consistent with radio-carbon dating of the Upper Palaeolithic transition (Mellars 2006b) and mtDNA lineage analysis (Olivieri et al. 2006) indicating an expansion from the Levant into Europe and North Africa 40–45 kya. We also see later growth in Australia (48 kya; fig. 1d) and New Guinea (39 kya; fig. 1g). In the Americas, we find evidence of a rapid population expansion beginning 18 kya (fig. 1h). This result is consistent with a hypothesized colonization event after the last glacial maximum 20 kya, when lower sea levels made it possible to cross from Siberia to North America via the Bering land bridge (Forster 2004).
By 35 kya global growth outside Africa begins to slow, most markedly in Southern Asia and Australia, whereas in Europe estimated effective population size actually declines slightly. This leveling of growth may reflect the impact of Malthusian environmental constraints and, particularly in Europe, increasing glaciation during this time. Climate also seems the most likely explanation for a second major growth phase inferred in Europe beginning 10–15 kya (fig. 1e), shortly after the end of the last ice age 15 kya. The estimated timing of this growth phase fits neatly with archaeological evidence for an expansion of hunter-gatherers from glacial refugia in southern Europe (Housley et al. 1997). Such an interpretation is also supported by previous work into the relative Neolithic and Palaeolithic contributions to the European mtDNA gene pool (Richards et al. 2000) and the discovery among the Saami (in northern Scandinavia) of southwestern European mtDNA haplotypes dated to the late Pleistocene (Achilli et al. 2005).
On a more recent timescale, the advent of agriculture is thought to have driven human population growth from 5–10 kya. Uncertainty in the timing of the second European expansion means, we cannot rule out the possibility that this growth was linked, at least in part, to the spread of agriculture from Anatolia beginning 9 kya (Gkiasta et al. 2003; Gray and Atkinson 2003). However, the inferred expansion time is centered somewhat earlier than would be expected under this scenario and an equally strong growth signal is not inferred for Africa, Papua New Guinea, the Americas, and the rest of Eurasia, where we would expect agriculture to have had a similar effect. This does not mean that no significant growth occurred with the development of early agriculture (although evidence of poor health in early agricultural populations makes this a possibility [Larsen 2006]). Rather, a more likely explanation is that due to the limits of the temporal resolution of the method, very recent expansion in human population size associated with agriculture has simply not left a growth signal that we can detect here. Although the resolution of the method is ultimately limited by the rate of molecular evolution, in the future multilocus analyses will allow a finer-grained analysis.
The results reported here extend our understanding of the origin and evolution of modern humans. By verifying our regional population size estimates against independent anthropological and historical evidence, we have shown that mtDNA variation can be used to predict population size in humans. Our estimates of effective population size through time show that Southern Asia was not only a key waypoint in the human expansion from Africa but also a major chapter in human prehistory. In addition, we find complex patterns of regional population growth that can be linked with archaeological evidence and known climatic events. As the amount of available sequence data continues to grow and results from a range of loci are combined, we anticipate that coalescent-based methods will reveal an increasingly detailed picture of human prehistory.
|
Supplementary Material |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Supplementary materials are available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/).
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
We are grateful to Mark Beaumont, Athena Ferreira, Andrew Meade, Mark Pagel, Marcel van de Steeg, and Chris Venditti for helpful advice and comments. We also thank the many contributors to GenBank. We also acknowledge support from Leverhulme Trust grant F/00 239/T.
|
Footnotes |
|---|
Beth Shapiro, Associate Editor
|
References |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResults and DiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Achilli A, Rengo C, Battaglia V, et al, (13 co-authors). Saami and Berbers–an unexpected mitochondrial DNA link. Am J Hum Genet (2005) 76:883–886.[CrossRef][Web of Science][Medline]
Alonso S, Armour JA. A highly variable segment of human subterminal 16p reveals a history of population growth for modern humans outstide Africa. Proc Natl Acad Sci USA (2001) 98:864–869.
Bazin E, Glemin S, Galtier N. Population size does not influence mitochondrial genetic diversity in animals. Science (2006) 312:570–572.
Biraben JN. Essay on the evolution of numbers of mankind. Population (1979) 34:13–25.
Birdsell JB. Microevolutionary patterns in aboriginal Australia: a gradient analysis of clines (1993) New York: Oxford University Press.
Cann RL, Stoneking M, Wilson AC. Mitochondrial DNA and human evolution. Nature (1987) 325:31–36.[CrossRef]
Cavalli-Sforza LL, Menozzi P, Piazza A. The history and geography of human genes (1994) Princeton (NJ): Princeton University Press.
Drummond AJ, Rambaut A. (2003) BEAST. Version 1.3, Oxford University Press, Oxford [computer programme]. Available from: http://beast.bio.ed.ac.uk/. Accessed 20 January 2007.
Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol (2005) 22:1185–1192.
Elson JL, Andrews RM, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Analysis of European mtDNAs for recombination. Am J Hum Genet (2001) 68:145–153.[CrossRef][Web of Science][Medline]
Eyre-Walker A. Evolution—size does not matter for mitochondrial DNA. Science (2006) 312:537–538.
Forster P. Ice ages and the mitochondrial DNA chronology of human dispersals: a review. Philos Trans R Soc Lond B Biol Sci (2004) 359:255–264. discussion 264.
Forster P, Matsumura S. Evolution. Did early humans go north or south? Science (2005) 308:965–966.
Friedlaender J, Schurr T, Gentz F, et al, (13 co-authors). Expanding southwest pacific mitochondrial haplogroups P and Q. Mol Biol Evol (2005) 22:1506–1517.
Gkiasta M, Russell T, Shennan S, Steele J. Neolithic transition in Europe: the radiocarbon record revisited. Antiquity (2003) 77:45–62.[Web of Science]
Goebel T. Anthropology. The missing years for modern humans. Science (2007) 315:194–196.
Gonser R, Donnelly P, Nicholson G, Di Rienzo A. Microsatellite mutations and inferences about human demography. Genetics (2000) 154:1793–1807.
Gray RD, Atkinson QD. Language-tree divergence times support the Anatolian theory of Indo-European origin. Nature (2003) 426:435–439.[CrossRef][Medline]
Harpending HC, Sherry ST, Rogers AR, Stoneking M. The genetic-structure of ancient human-populations. Curr Anthropol (1993) 34:483–496.[CrossRef][Web of Science]
Ho SY, Phillips MJ, Cooper A, Drummond AJ. Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol (2005) 22:1561–1568.
Housley RA, Gamble CS, Street M, Pettitt P. Radiocarbon evidence for the Lateglacial human recolonisation of northern Europe. Proc Prehistoric Soc (1997) 63:25–26.
Howell N, Kubacka I, Mackey DA. How rapidly does the human mitochondrial genome evolve? Am J Hum Genet (1996) 59:501–509.[Web of Science][Medline]
Hudson RR. Gene genealogies and the coalescent process. Oxf Surv Evol Biol (1990) 7:1–14.
Huoponen K, Schurr TG, Chen Y, Wallace DC. Mitochondrial DNA variation in an aboriginal Australian population: evidence for genetic isolation and regional differentiation. Hum Immunol (2001) 62:954–969.[CrossRef][Web of Science][Medline]
Ingman M, Gyllensten U. Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines. Genome Res (2003) 13:1600–1606.
Ingman M, Kaessmann H, Paabo S, Gyllensten U. Mitochondrial genome variation and the origin of modern humans. Nature (2000) 408:708–713.[CrossRef][Medline]
Kivisild T, Rootsi S, Metspalu M, et al, (18 co-authors). The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations. Am J Hum Genet (2003) 72:313–332.[CrossRef][Web of Science][Medline]
Kivisild T, Shen P, Wall DP, et al, (17 co-authors). The role of selection in the evolution of human mitochondrial genomes. Genetics (2006) 172:373–387.
Larsen CS. The agricultural revolution as environmental catastrophe: implications for health and lifestyle in the Holocene. Quatern Int (2006) 150:12–20.[CrossRef]
Macaulay V, Hill C, Achilli A, et al, (21 co-authors). Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science (2005) 308:1034–1036.
Marth GT, Czabarka E, Murvai J, Sherry ST. The allele frequency spectrum in genome-wide human variation data reveals signals of differential demographic history in three large world populations. Genetics (2004) 166:351–372.
Mellars P. Going east: new genetic and archaeological perspectives on the modern human colonization of Eurasia. Science (2006a) 313:796–800.
Mellars P. A new radiocarbon revolution and the dispersal of modern humans in Eurasia. Nature (2006b) 439:931–935.[CrossRef][Medline]
Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E. Equations of state calculations by fast computing machines. J Chem Phy (1953) 21:1087–1091.[CrossRef]
Mishmar D, Ruiz-Pesini E, Golik P, et al, (13 co-authors). Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci USA (2003) 100:171–176.
Newton MA, Raftery AE, Davison AC, et al, (33 co-authors). Approximate bayesian-inference with the weighted likelihood bootstrap. J R Stat Soc Ser B Methodol (1994) 56:3–48.
O'Connell JF, Allen J. Dating the colonization of Sahul (Pleistocene Australia-New Guinea): a review of recent research. J Arch Sci (2004) 31:835–853.[CrossRef][Web of Science]
Olivieri A, Achilli A, Pala M, et al, (15 co-authors). The mtDNA legacy of the Levantine early Upper Palaeolithic in Africa. Science (2006) 314:1767–1770.
Richards M, Macaulay V, Hickey E, et al, (37 co-authors). Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet (2000) 67:1251–1276.[Web of Science][Medline]
Rogers AR. Genetic-evidence for a pleistocene population explosion. Evolution (1995) 49:608–615.[CrossRef][Web of Science]
Rogers AR. Order emerging from chaos in human evolutionary genetics. Proc Natl Acad Sci USA (2001) 98:779–780.
Sherry ST, Rogers AR, Harpending H, Soodyall H, Jenkins T, Stoneking M. Mismatch distributions of mtDNA reveal recent human population expansions. Hum Biol (1994) 66:761–775.[Web of Science][Medline]
Steel MA, Hendy MD, Penny D. Loss of information in genetic distances. Nature (1988) 336:118.[CrossRef][Medline]
Sun C, Kong QP, Palanichamy MG, Agrawal S, Bandelt HJ, Yao YG, Khan F, Zhu CL, Chaudhuri TK, Zhang YP. The dazzling array of basal branches in the mtDNA macrohaplogroup M from India as inferred from complete genomes. Mol Biol Evol (2006) 23:683–690.
Thangaraj K, Chaubey G, Kivisild T, Reddy AG, Singh VK, Rasalkar AA, Singh L. Reconstructing the origin of Andaman Islanders. Science (2005) 308:996.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res (1997) 25:4876–4882.
van Holst Pellekaan SM, Ingman M, Roberts-Thomson J, Harding RM. Mitochondrial genomics identifies major haplogroups in Aboriginal Australians. Am J Phys Anthropol (2006) 131:282–294.[CrossRef][Web of Science][Medline]
Vigilant L, Stoneking M, Harpending H, Hawkes K, Wilson AC. African populations and the evolution of human mitochondrial DNA. Science (1991) 253:1503–1507.
Wakeley J. The effects of subdivision on the genetic divergence of populations and species. Evol Int J Org Evol (2000) 54:1092–1101.
Watson E, Forster P, Richards M, Bandelt HJ. Mitochondrial footprints of human expansions in Africa. Am J Hum Genet (1997) 61:691–704.[Web of Science][Medline]
Zhivotovsky LA, Bennett L, Bowcock AM, Feldman MW. Human population expansion and microsatellite variation. Mol Biol Evol (2000) 17:757–767.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  W. Shi, Q. Ayub, M. Vermeulen, R.-g. Shao, S. Zuniga, K. van der Gaag, P. de Knijff, M. Kayser, Y. Xue, and C. Tyler-Smith A Worldwide Survey of Human Male Demographic History Based on Y-SNP and Y-STR Data from the HGDP-CEPH Populations Mol. Biol. Evol., February 1, 2010; 27(2): 385 - 393. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 R. Foley and C. Gamble The ecology of social transitions in human evolution Phil Trans R Soc B, November 12, 2009; 364(1533): 3267 - 3279. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. Petraglia, C. Clarkson, N. Boivin, M. Haslam, R. Korisettar, G. Chaubey, P. Ditchfield, D. Fuller, H. James, S. Jones, et al. Population increase and environmental deterioration correspond with microlithic innovations in South Asia ca. 35,000 years ago PNAS, July 28, 2009; 106(30): 12261 - 12266. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 A. W. Briggs, J. M. Good, R. E. Green, J. Krause, T. Maricic, U. Stenzel, C. Lalueza-Fox, P. Rudan, D. Brajkovic, Z. Kucan, et al. Targeted Retrieval and Analysis of Five Neandertal mtDNA Genomes Science, July 17, 2009; 325(5938): 318 - 321. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Powell, S. Shennan, and M. G. Thomas Late Pleistocene Demography and the Appearance of Modern Human Behavior Science, June 5, 2009; 324(5932): 1298 - 1301. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
K. B. Schroeder, M. Jakobsson, M. H. Crawford, T. G. Schurr, S. M. Boca, D. F. Conrad, R. Y. Tito, L. P. Osipova, L. A. Tarskaia, S. I. Zhadanov, et al. Haplotypic Background of a Private Allele at High Frequency in the Americas Mol. Biol. Evol., May 1, 2009; 26(5): 995 - 1016. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 Q. D Atkinson, R. D Gray, and A. J Drummond Bayesian coalescent inference of major human mitochondrial DNA haplogroup expansions in Africa Proc R Soc B, January 22, 2009; 276(1655): 367 - 373. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
B. M. Henn, C. R. Gignoux, M. W. Feldman, and J. L. Mountain Characterizing the Time Dependency of Human Mitochondrial DNA Mutation Rate Estimates Mol. Biol. Evol., January 1, 2009; 26(1): 217 - 230. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||