Calculating expected DNA remnants from ancient founding events in human population genetics
© Stacey et al; licensee BioMed Central Ltd. 2008
Received: 10 July 2008
Accepted: 17 October 2008
Published: 17 October 2008
Recent advancements in sequencing and computational technologies have led to rapid generation and analysis of high quality genetic data. Such genetic data have achieved wide acceptance in studies of historic human population origins and admixture. However, in studies relating to small, recent admixture events, genetic factors such as historic population sizes, genetic drift, and mutation can have pronounced effects on data reliability and utility. To address these issues we conducted genetic simulations targeting influential genetic parameters in admixed populations.
We performed a series of simulations, adjusting variable values to assess the affect of these genetic parameters on current human population studies and what these studies infer about past population structure. Final mean allele frequencies varied from 0.0005 to over 0.50, depending on the parameters.
The results of the simulations illustrate that, while genetic data may be sensitive and powerful in large genetic studies, caution must be used when applying genetic information to small, recent admixture events. For some parameter sets, genetic data will not be adequate to detect historic admixture. In such cases, studies should consider anthropologic, archeological, and linguistic data where possible.
In the past 20 years, DNA sequence data and advanced computational techniques have provided an unparalleled resource in the study of human origins and migration. These tools have demonstrated a Pleistocene colonization of America by Asian populations[3, 4] and have even prompted calculations of the size of the original human founding populations. Similarly, DNA sequence data have helped demonstrate the dynamics of large human populations such as primitive human migration out of Africa, the American migration, the Lemba migration in Africa, the migratory history of the Baltic States, and many others. Researchers have even used the population genetics of human disease vectors to trace human migration events. It may be difficult to underestimate the value genetic data have played and will continue to play on our ability to reconstruct historic population events.
But while sequence data have been used to study many forms of human migration, their utility in the study of small-scale migration is still in question. Research into small migrations like the Norse settlements in Greenland, a possible Polynesian migration to the New World, the North African Slave migration to America, and the pre-Columbian European migration to America[13, 14], have traditionally been based primarily on evidence other than DNA sequence information. However, recently, researchers have begun to apply genetic data to these smaller historical migrations and make conclusions about small historic populations using current DNA. For example, DNA information has recently been used to study the small indigenous populations of Tierra del Fuego, and to analyze Caucasian admixture in specific African American populations[16, 17]. It should be noted that genetic data have been used to study the large Norse migration to Ireland, but are an afterthought when researching their short-lived occupation of Canada[19, 20].
This raises questions about the utility of genetic data in providing evidence for historic migrations and inferences of unknown past events. While genetic studies can provide considerable information, they are also accompanied by variation and stochasticity. Because of these limitations, even the most complete studies of human populations have been called "not unequivocal" or "sobering" by those conducting the research. Recent reports have also addressed the limited depth of current genetic studies, indicating that most studies make conclusions after sequencing less than 1% of subjects' genomes, and sampling only small numbers of a population. Such methods can be especially problematic when dealing with historic admixture events that are very small. The difficulty is a function of the current architecture of genetic studies: researchers sample loci from a group of individuals and categorize individuals into groups based on which alleles they have at the loci tested[24, 25]. These categorizations are determined based on the most prevalent or probable genetic markers in an individual's genome. The results of these studies, then, can overlook genetic markers that simply are not sampled, which is common in small admixture events. Additionally, stochastic events can lead to allele fixation and further complicate matters, particularly in small populations. It has been suggested that studies of even the largest migrations should couple genetic information with archeological, anthropological, and linguistic data.
As our ability to collect and analyze DNA sequence data increases, understanding the probabilities and variability associated with admixture becomes especially important. In this study, we explore the utility of DNA sequence data in small, recent human migration studies. We use forward-based genetic simulation to explore three questions: 1) what variables contribute to the presence (or absence) of historic markers in today's genomes, 2) how do these variables affect the probability of finding historically admixed DNA in today's populations, and 3) how can studies be designed to maximize information from genetic data? These questions are answered through genetic simulation and a sample size study aimed at suggesting the numbers of subjects and loci that should be sampled to successfully detect small-scale admixture. In our simulations, we assume that migrant allele frequencies are known a priori. The simulations test our ability to detect these known migrant alleles in admixed descended populations. We find that genetic parameters, the stochasticity of genetic drift, and experimental design all play an important role in the ability to find historic DNA in current admixed populations.
We used the simuPOP software package for forward-based genetic simulations. In each simulation, a "migrant" population with distinct, known alleles was admixed with a "native" population. We followed the combined population through time and recorded the frequency of migrant alleles at each generation. Because migrant genetic parameters were known a priori, these simulated allele frequencies allow us to assess how parameters affect the ability of detecting migrant alleles in an admixed descendant population. We used a generation time of 23 years as a compromise among differing estimates of human generation times [28–30]. The simuPOP module allows numerous genetic variables to be altered and studied independently. The variables of interest in these initial simulations are basic genetic variables: native population size, migrant population size, mutation rate, time since admixture event, and initial allele frequencies. These variables allow the assessment of the role that population sizes, mutation, genetic drift, and allele frequency have on the amount of migrant DNA present in the admixed population after a number of generations. Our simulations have been designed so that total population sizes are as analogous to effective population sizes (Ne) as possible. We assume that each individual has an equal expectation of obtaining progeny, that there are equal sex ratios, and that the population remains constant over time. These assumptions allow the population size used in our study to be interpreted as an effective population size, though under some definitions of Ne our numbers will have different values of Ne than those assigned. The statistics and results in this study are based on the allele frequencies retrieved from the simuPOP software. We imported these numbers into the R statistical package for numeric and graphical analysis.
In our genetic simulations, we make a number of assumptions about the populations: random mating, absence of selection, no gene flow, and constant population size from time of the migratory event to the present. Actual populations experience some gene flow with neighboring populations[32, 33]; however, in our simulations, we do not consider this in an attempt to create a best-case scenario for the migrant allele. If such gene flow did occur, it could only decrease the chances of detecting the migration event by lowering the frequency of the migratory allele in the admixed population. In addition, real populations often experience growth following admixture. However, assuming that the migrant allele is growing at the same rate as the other alleles (random mating), the allele frequency should not be changed directly by population size increase, although the effects of drift could become less pronounced as a result of a greater population size. Further studies and simulations using population growth rates may be helpful in addressing the effects of population growth.
Our simulations can be grouped into two separate categories. The first is a series of simulations designed to assess how the parameters mentioned above can influence the presence of historic migrant DNA in today's populations. More concretely, these simulations answer this question: how does each genetic parameter affect the frequency of migrant alleles in an admixed population? Our simulations tested the effect of 4 variables: size of migrant population, size of native population, time since admixture event, and mutation rate at the locus of interest. We assigned each variable a high value and a low value based on current literature and ran a total of 16 simulations using a full factorial experimental design, altering only one variable at a time. This allowed us to study variables independently and assess how they affect the frequency of the migrant allele over time. We compare the impact of each variable by holding other variables constant and comparing the frequencies of the migrant allele.
In the first simulations, we modeled only one locus per individual and assumed no recombination; one locus is adequate to assess the role of these parameters on allele frequencies. We also initialized the migrant population with the migrant allele fixed (all migrant individuals possessed the migrant allele). This is unrealistic, but provides a best-case scenario for detecting the migrant allele. We replicated each simulation 250 times.
The second simulation category was a single simulation designed to mimic the genetic landscape of a true admixed population. We assigned mid-range values for migrant population size (200), native population size (5,000), and generations (100). In order to more realistically model a current study, we followed 1,000 loci on 20 different chromosomes on each individual. This represents a sample much larger than the recommended number needed in order to detect large human admixture. A standard recombination rate of 1.26 cM/Mb was used, though the human recombination rate has been shown to be negligible over 100 generations. At the beginning of the simulation, an initial migrant allele frequency and a mutation rate were randomly generated for each of the 1,000 loci on each individual, in order to model the DNA seen in actual human genetics research. The methods of random generation are outlined below.
Sample size study
In order to understand what must be done to successfully study data from historic admixture, we constructed a sample size study using the data from simulation 2. Small human genetics studies test approximately 50 loci when studying populations. Given the calculated frequency of migrant alleles in our simulated population, we calculated the number of migrant alleles that would be seen, on average, in each human subject of a genetic study. This is accomplished using the cumulative density function (CDF) of a binomially distributed random variable where the size parameter is 50 and the probability parameter is the expected migrant allele frequency. In comparison, one of the larger human genetic studies to date sequenced 993 loci in each human subject. Accordingly, we followed the same protocol to investigate a study of this magnitude, using the binomial CDF with a size parameter of 993 and the same probability parameter.
The most recent studies have again raised the bar as far as loci per subject, sampling 650,000 loci in each individual. Although sampling more loci will find a larger number of migrant alleles, the proportion of such markers in the population does not change when more samples are taken. The study conducted by Li et al. (2008) samples about 20 individuals per population group, a number similar to previous studies. Accordingly, we investigated the sample size necessary to find at least one migrant allele at each of the loci sequenced in a large genetic study.
Simulation study 1
Simulation study 2
Sample size study
In the case of a large study with as many as 933 loci, based upon the expected migrant allele frequency of 1.017%, almost every subject would demonstrate at least one migrant allele (Figure 7b). In fact, most subjects would demonstrate more than 9 migrant alleles. However, while large studies would expect to succeed in finding more migrant alleles in today's population, this alone cannot link the admixed population to the migrant population. The migrant alleles will still only represent, on average, 1% of every allele sequenced in the entire study. Therefore, although 9 migrant alleles may, on average, be found in each subject, it is hard to know if the migrant alleles will be redundant among loci and subjects or spread evenly throughout all the loci in the study. Additionally, these numbers could be considerably lower depending on the allele frequency in the migrating population.
Our results provide some important insights in detecting historic admixture. The simulations we present illustrate the effect that initial parameters have on the outcome of human admixture. Simple adjustments in the parameters in our simulation series changed the expected allele frequency outcome from as low as 0.0005 to over 0.50, an increase of three orders of magnitude. The results of any admixture study using genetic data, then, are highly dependent on the variables presented in these simulations (e.g., mutation rate, population sizes, and time since admixture (number of generations)).
High mutation rates can decrease the expected migrant allele frequency and the variability by more than 50 percent, especially in populations that experienced earlier migrations. For example, an increased mutation rate can change the mean final allele frequency from .0243 to .0128, or from .5016 to .2699 (depending on other variables, as reported in Figure 3). Researchers should keep this in mind when selecting loci for analysis. Because some DNA mutation rates are highly variable, choice of locus can have a profound impact on the number of migrant alleles detected years later. Many studies advocate the use of mtDNA due to data collecting feasibility and other factors. However, because the mutation rate is generally higher in mtDNA, it could corrupt signal in studies addressing historic admixture, even when the time frame is relatively recent.
The sizes of the migrant and native populations are fundamental for an understanding of expected allele frequency. With time since admixture as low as those we consider in our simulations, the most important factors are the sizes of the migrating and native populations. In our simulations, if the native population is large, changing the migrating population size results in a change of mean final allele frequency from .0243 to .0010. If the native population is small, those numbers change to .5016 and .0407. These are the most significant differences illustrated by our simulations and they attest to the important role of population sizes. Researchers should not expect to find many alleles from a small migratory group of 50 individuals in a large population today, even if sampling methods are exhaustive.
Additionally, we see that time plays an important role. The standard deviations presented in Table 1 demonstrate that allelic frequencies vary widely, particularly as the number of generations increases. High mutation rates combined with large time spans can reduce migrant allele frequencies significantly. When the mutation rate is low, however, the time since admixture does not affect the final mean allele frequency much (or at all), but it still has a profound impact on the standard deviation. For example, a change in time since admixture in one parameter set almost doubles the standard deviation from .0525 to .1044. As time increases, genetic drift causes the spread of final allele frequencies to increase, particularly when the population sizes are small. Thus, as the time since the admixture event increases, sample size for both loci and subjects becomes increasingly important.
In our second simulation, most of the migrant alleles are present in less than 2% of the population. In a study of a population where few subjects from many human populations are studied, alleles from a small-scale admixture will usually not be recovered at all. And these rare alleles could easily be ignored in favor of haplotypes that better categorize the population into clusters.
Our results demonstrate a profound and general fact: the values of these genetic parameters can drastically alter the expected frequency of migrant alleles in today's populations. Even in our simulations, where steps have been taken to ensure a best-case scenario for the migrant allele, there is often a large spread of possible outcomes. DNA data have been touted as a panacea for recovering information about the past, but their use depends so extensively on factors that are beyond our control that their application is not always appropriate. It is imperative, therefore, that researchers understand the implications of the variables we have presented and not rely solely on DNA sequence data when researching small, recent human migrations. We can only hope to understand basic details of population history when quantifying genetic data and even valid results derived from genetic data may still be misleading if viewed unilaterally, as demonstrated by Harpending et al [46, 47].
Improving probability of detecting historic admixture
• Large Migrant Population
• Identify informative migrant alleles
• Small Native Population
• Test large number of loci
• Low mutation rate at loci of interest
• Large sample size for each population
• Fewer generations since admixture event
• Establish methods for detecting rare alleles
• Collaborative approach (Archeology, Anthropology, Linguistics)
Perhaps most importantly, it must be remembered that drift is stochastic and that historic genetic parameters are, for the most part, unknown. Thus, the absence of specific genetic data is not conclusive evidence against historic admixture. Our results illustrate several parameter sets that would cause admixture to be either completely or practically undetectable today. To address the inconsistent results found in DNA all but the largest genetic studies need to continue to consider anthropologic, archeological, and linguistic data in order to formulate conclusions. Finally, our study demonstrates the utility of simulation studies to put bounds on parameter values and sample sizes for studies of human migration events.
The ability to detect historic admixture and make correct inferences based on genetic data depends on the interplay between population sizes, mutation rates, time, and other parameters. We explore the parameter space of historic alleles in current populations and demonstrate the broad implications of each of these genetic parameters on modern allele frequencies. Our results provide guidelines with respect to the population genetic parameters and their values needed to detect migrant alleles in an admixed population. While studies that focus on large admixture events should be able to draw specific and valid conclusions, we suggest that genetic data be used with caution when studying small admixture events. The random nature of admixed genetic data seen in these simulations demonstrates that the utility of genetic data is dependent on the context of each individual study. Increasing the number of loci and the number of individuals sampled will increase the probability of detecting small traces of signal, but other sources of evidence should always be considered where possible.
We thank Ryan Parr for comments on an earlier draft of this manuscript. This work was supported by an Eliza R. Snow Fellowship from Brigham Young University.
- Ebersberger I, Galgoczy P, Taudien S, Taenzer S, Platzer M, von Haeseler A: Mapping human genetic ancestry. Molecular Biology and Evolution. 2007, 24 (10): 2266-2276. 10.1093/molbev/msm156.View ArticlePubMedGoogle Scholar
- Templeton A: Out of Africa again and again. Nature. 2002, 416: 45-51. 10.1038/416045a.View ArticlePubMedGoogle Scholar
- Greenberg JH, II CGT, Zegura SL, Campbell L, Fox JA, Laughlin WS, Szathmary EJE, Weiss KM, Wollford E: The Settlement of the Americas: A Comparison of the Linguistic, Dental, and Genetic Evidence. Current Anthropology. 1986, 27 (5): 477-497. 10.1086/203472.View ArticleGoogle Scholar
- Goebel T: Pleistocene Human Colonization of Siberia and Peopling of the Americas: An Ecological Approach. Evolutionary Anthropology. 1999, 8 (6): 208-227. 10.1002/(SICI)1520-6505(1999)8:6<208::AID-EVAN2>3.0.CO;2-M.View ArticleGoogle Scholar
- Hey J: On the Number of New World Founders: A Population Genetic Portrait of the Peopling of the Americas. PlOS Biology. 2005, 3 (6): e193-10.1371/journal.pbio.0030193.PubMed CentralView ArticlePubMedGoogle Scholar
- Fagundes LJR, Ray N, Beaumont M, Neuenschwaner S, Salzano FM, Bonatto SL, Excoffier L: Statistical Evaluation of Alternative Models of Human Evolution. PNAS. 2007, 104 (45): 17614-17619. 10.1073/pnas.0708280104.PubMed CentralView ArticlePubMedGoogle Scholar
- Thomas MG, Parfill T, Weiss DA, Skorecki K, Wilson JF, Roux Ml, Bradman N, Goldstein DB: Y Chromosomes Traveling Sough: The Cohen Modal Haplotype and the Origins of the Lemba – the "Black Jews of Southern Africa". Am J Hum Genet. 2000, 66: 674-686. 10.1086/302749.PubMed CentralView ArticlePubMedGoogle Scholar
- Beckman L, Sikstrom C, Midelsaar AV, Krumina A, Ambrasiene D, Kucinskas V, Beckman G: Tranferrin Varients as Markers of Migrations and Admixtures between Populations in the Baltic Sea Region. Human Heredity. 1998, 48: 185-191. 10.1159/000022800.View ArticlePubMedGoogle Scholar
- Falush D, Wirth T, Linz B, Pritchard JK, Stephens M, Kidd M, Blaser MJ, Graham DY, Vacher S, Perez-Perez GI, et al: Traces of human migrations in Helicobacter pylori populations. Science. 2003, 299: 1582-1585. 10.1126/science.1080857.View ArticlePubMedGoogle Scholar
- McGhee R: Contact between Native Americans and the Medieval Norse: A Review of the Evidence. American Antiquity. 1984, 49 (1): 4-26. 10.2307/280509.View ArticleGoogle Scholar
- Klar KA: Linguistic Evidence for a Prehistoric Polynesia-Southern California Contact Event. Anthropological LInguistics. 2005, 47 (4): 369-400.Google Scholar
- Allen SD: More on the Free Black Population of the Southern Appalachian Mountains: Speculations on the North African Connection. Journal of Black Studies. 1995, 25 (6): 651-671. 10.1177/002193479502500601.View ArticleGoogle Scholar
- Epstein JF, Buchanan DB, Buttrey TV, Carter GF, Cook; WL, Covey C, Jett SC, Thomas A, Lee J, Mundkur B, Paulsen AC, et al: Pre-Columbian Old World Coins in America: An Examination of the Evidence. Current Anthropology. 1980, 21 (1): 1-20. 10.1086/202398.View ArticleGoogle Scholar
- Quaife MM: The Myth of the Kensington Rune Stone: The Norse Discovery of Minnesota 1362. The New England Quarterly. 1934, 7 (4): 613-645. 10.2307/359189.View ArticleGoogle Scholar
- Lalueza C, Perez-Perez A, Prats E, Cornudella L, Turbon D: Lack of founding Amerindian midochondrial DNA lineages in extinct Aborigines from Tierra del Fuego-Patagonia. Human Molecular Genetics. 1997, 6 (1): 41-46. 10.1093/hmg/6.1.41.View ArticlePubMedGoogle Scholar
- Parra EJ, Kittles RA, Argyropoulos G, Pfaff CL, Hiester K, Bonilla C, Sylvester N, Parrish-Gause D, Garvey WT, Jin L, et al: Ancestral Proportions and Admixture Dynamics in Geographically Defined African Americans Living in South Carolina. American Journal of Physical Anthropology. 2001, 114: 18-29. 10.1002/1096-8644(200101)114:1<18::AID-AJPA1002>3.0.CO;2-2.View ArticlePubMedGoogle Scholar
- Parra EJ, Marcini A, Akey J, Martinson J, Batzer MA, Cooper R, Forrester T, Allison DB, Deka R, Ferrell RE, et al: Estimating African American Admixture Proportions by Use of Population Specific Alleles. Am J Hum Genet. 1998, 63: 1839-1851. 10.1086/302148.PubMed CentralView ArticlePubMedGoogle Scholar
- McEvoy B, Brady C, Moore LT, Bradley DG: The Scale and Nature of Viking Settlement in Ireland from Y-chromosome Admixture Analysis. European Journal of Human Genetics. 2006, 14: 1288-1294. 10.1038/sj.ejhg.5201709.View ArticlePubMedGoogle Scholar
- Nansen F: The Norsemen in America. The Geographical Journal. 1911, 38 (6): 557-575. 10.2307/1778837.View ArticleGoogle Scholar
- Wallace B: The Norse in Newfoundland: L'Anse aux Meadows and Vinland. Newfoundland and Labrador Studies. 2003, 19 (1):
- Eshleman JA, Malhi RS, Smith DG: Mitochondrial DNA Studies of Native Americans: Conceptions and MIsconceptions of the Population Prehistory of the Americas. Evolutionary Anthropology. 2003, 12: 7-18. 10.1002/evan.10048.View ArticleGoogle Scholar
- Cann RL: Genetic Clues to Dispersal in Human Populations: Retracing the Past from the Present. Science. 2001, 291: 1742-1748. 10.1126/science.1058948.View ArticlePubMedGoogle Scholar
- Bolnick DA, Fullwiley D, Duster T, Cooper RS, Fujimura JH, Kahn J, Kaufman JS, Marks J, Morning A, Nelson A, et al: The Science and Business of Genetic Ancestry Testing. Science. 2007, 318: 399-400. 10.1126/science.1150098.View ArticlePubMedGoogle Scholar
- Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldmen MW: Clines, Clusters, and the Effect of Study Design on the Inferene of Human Population Structure. PLOS Genetics. 2005, 1 (6): e70-10.1371/journal.pgen.0010070.PubMed CentralView ArticlePubMedGoogle Scholar
- Li JZ, Absher DM, Tang H, Southwick AM, Castro AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, et al: Worldwide Human Relationships Inferred from Genome-Wide Patterns of Variation. Science. 2008, 319: 1100-1103. 10.1126/science.1153717.View ArticlePubMedGoogle Scholar
- Schurr TG: The Peopling of The New World: Perspectives from Molecular Anthropology. Annu Rev Anthropol. 2004, 33: 551-583. 10.1146/annurev.anthro.33.070203.143932.View ArticleGoogle Scholar
- Peng B, Kimmel M: simuPOP: a forward-time population genetics simulation envirnment. Bioinformatics. 2005, 21 (18): 3686-3687. 10.1093/bioinformatics/bti584.View ArticlePubMedGoogle Scholar
- Sherry ST, Harpending HC, Batzer MA, Stoneking M: Alu Evolution in Human Populations: Using the Coalescent to Estimate Effective Population Size. Genetics. 1997, 147: 1977-1982.PubMed CentralPubMedGoogle Scholar
- Hobolth A, Christensen OF, Mailund T, Schierup MH: Genomic Relationships and Speciation Times of Human, CHimpanzeee, and Gorilla Inferred from a Coalescent Hidden Markov Model. PLOS Genetics. 2007, 3 (2): e7-10.1371/journal.pgen.0030007.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Y-X: Estimating Mutation Rate and Generation Time from Longitudinal Samples of DNA Sequences. Mol Biol Evol. 2001, 18 (4): 620-626.View ArticlePubMedGoogle Scholar
- Hedrick P: Genetics of Populations. 2000, Jones & Bartlett Publishers, 2Google Scholar
- Wright S: Isolation by Distance. Genetics. 1943, 28 (2): 114-138.PubMed CentralPubMedGoogle Scholar
- Templeton AR: Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history. Molecular Ecology. 1998, 7: 381-397. 10.1046/j.1365-294x.1998.00308.x.View ArticlePubMedGoogle Scholar
- Evans SN, Shvets Y, Slatkin M: Non-equilibrium theory of the allele frequency spectrum. Theoretical Population Biology. 2006, 71: 109-119. 10.1016/j.tpb.2006.06.005.View ArticlePubMedGoogle Scholar
- Haviland WA: Family Size, Prehistoric Population Estimates, and the Ancient Maya. American Antiquity. 1972, 37 (1): 135-139. 10.2307/278895.View ArticleGoogle Scholar
- Spurdle AB, Jenkins T: The Origins of the Lemba "Black Jews" of Southern Africa: Evidence from p12E2 and Other T-Chromosome Markers. Am J Hum Genet. 1996, 59: 1126-1133.PubMed CentralPubMedGoogle Scholar
- Bedoya G, Montoya P, Garcý J, Soto I, Bourgeois S, Carvajal L, Labuda D, Alvarez V, Ospina J, Hedrick PW, et al: Admixture dynamics in Hispanics: A shift in the nuclear genetic ancestry of a South American population isolate. PNAS. 2006, 103 (19): 7734-7239. 10.1073/pnas.0508716103.View ArticleGoogle Scholar
- Zegura SL, Karafet TM, Zhivotovsky LA, Hammer MF: High-Resolution SNPs and Microsatellite Haplotypes Point to a Single, Recent Entry of Native American Y Chromosomes into the Americas. Mol Biol Evol. 2004, 21 (1): 164-175. 10.1093/molbev/msh009.View ArticlePubMedGoogle Scholar
- Horai S, Kondo R, Nakagawa-Hattori Y, Hayashi S, Sonoda S, Tajima K: Peopling of the Americas, Founded by Four Major Lineages of Mitochondrial DNA. Mol Biol Evol. 1993, 10 (1): 23-47.PubMedGoogle Scholar
- Bonatto SL, Salzano FM: A single and early migration for the peopling of the Americas supported by mitochondrial DNA sequence data. Proc Natl Acad Sci. 1997, 94: 1866-1871. 10.1073/pnas.94.5.1866.PubMed CentralView ArticlePubMedGoogle Scholar
- Wall JD: Detecting Ancient Admixture in Humans Using Sequence Polymorphism Data. Genetics. 2000, 154 (3): 1271-1279.PubMed CentralPubMedGoogle Scholar
- Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, Chen C-F, Thomas MA, Haussler DH, Jacob HJ: Comparative Recombination Rates in the Rat, Mouse, and Human Genomes. Genome Research. 2004, 14: 528-538. 10.1101/gr.1970304.PubMed CentralView ArticlePubMedGoogle Scholar
- Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES: High-resolution haplotype structure in the human genome. Nature Genetics. 2001, 29: 229-232. 10.1038/ng1001-229.View ArticlePubMedGoogle Scholar
- Basu A, Mukherjee N, Roy S, Sengupta S, Banerjee S, Chakraborty M, Dey B, Roy M, Roy B, Bhattacharyya NP, et al: Ethnic India: A Genomic View, With Special Reference to Peopling and Structure. Genome Research. 2007, 13: 2277-2290. 10.1101/gr.1413403.View ArticleGoogle Scholar
- Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldman MM: Clines, clusters, and the effect of study design on the inference of human population structure. PLoS Genetics. 2005, 1 (6): 660-671. 10.1371/journal.pgen.0010070.View ArticleGoogle Scholar
- Harpending H, Rogers A: Genetic Perspectives on Human Origins and Differentiation. Annu Rev Genomics Hum Genet. 2000, 1: 361-385. 10.1146/annurev.genom.1.1.361.View ArticlePubMedGoogle Scholar
- Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, Sherry ST: Genetic Traces of Ancient Demography. Proc Natl Acad Sci. 1998, 95: 1961-1967. 10.1073/pnas.95.4.1961.PubMed CentralView ArticlePubMedGoogle Scholar
- Nachman MW, Crowell SL: Estimate of the Mutation Rate per Nucleotide in Humans. Genetics. 2000, 156: 297-304.PubMed CentralPubMedGoogle Scholar
- Heyer E, Puymirat J, Dieltjes P, Bakker E, de Knijff P: Estimating Y chromosome specific microsatellite mutation frequencies using deep rooting pedigrees. Human Molecular Genetics. 1997, 6 (5): 799-803. 10.1093/hmg/6.5.799.View ArticlePubMedGoogle Scholar
- Kayser M, Roewer L, Hedman M, Henke L, Henke J, Brauer S, Kruger C, Krawczak M, Nagy M, Dobosz T, et al: Characteristics and Frequency of Germline Mutations at Microsatellite Loci from the Human Y Chromosome, as Revealed by Direct Observation in Father/Son Pairs. Am J Hum Genet. 2000, 66: 1580-1588. 10.1086/302905.PubMed CentralView ArticlePubMedGoogle Scholar
- Anagnostopoulos T, Green PM, Rowley G, Lewis CM, Giannelli F: DNA Variation in a 5-Mb Region of the X Chromosome and Estimates of Sex-Specific/Type-Specific Mutation Rates. Am J Hum Genet. 1999, 64: 508-517. 10.1086/302250.PubMed CentralView ArticlePubMedGoogle Scholar
- Calafell F, Shuster A, Speed WC, Kidd JR, Kidd KK: Short tandem repeat polymorphism evolution in humans. European Journal of Human Genetics. 1998, 6: 38-49. 10.1038/sj.ejhg.5200151.View ArticlePubMedGoogle Scholar
- Brinkmann B, Klintschar M, Neuhuber F, Huhne J, Rolf B: Mutation Rate in Human Microsatellites: Influence of the Struture and Length of Tandem Repeat. Am J Hum Genet. 1998, 62: 1408-1415. 10.1086/301869.PubMed CentralView ArticlePubMedGoogle Scholar
- Whittaker JC, Harbord RM, Boxall N, Mackay I, Dawson G, Sibly RM: Likelihood-Based Estimation of Microsatellite Mutation Rates. Genetics. 2003, 164: 781-787.PubMed CentralPubMedGoogle Scholar
- Slatkin M: A Measure of Population Subdivision Based on Microsatellite Allele Frequencies. Genetics. 1995, 139: 457-462.PubMed CentralPubMedGoogle Scholar
- Sigurðardottir S, Helgason A, Gulcher JR, Stefansson K, Donnelly P: The Mutation Rate in the Human mtDNA Control Region. Am J Hum Genet. 2000, 66: 1599-1609. 10.1086/302902.View ArticleGoogle Scholar
- Khrapko K, Coller HA, Andre PC, Li X-C, Hanekamp JS, Thilly WG: Mitochondrial mutational spectra in human cells and tissues. PNAS. 1997, 94: 13798-13803. 10.1073/pnas.94.25.13798.PubMed CentralView ArticlePubMedGoogle Scholar