Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data

Abstract

Demographic models built from genetic data play important roles in illuminating prehistorical events and serving as null models in genome scans for selection. We introduce an inference method based on the joint frequency spectrum of genetic variants within and between populations. For candidate models we numerically compute the expected spectrum using a diffusion approximation to the one-locus, two-allele Wright-Fisher process, involving up to three simultaneous populations. Our approach is a composite likelihood scheme, since linkage between neutral loci alters the variance but not the expectation of the frequency spectrum. We thus use bootstraps incorporating linkage to estimate uncertainties for parameters and significance values for hypothesis tests. Our method can also incorporate selection on single sites, predicting the joint distribution of selected alleles among populations experiencing a bevy of evolutionary forces, including expansions, contractions, migrations, and admixture. We model human expansion out of Africa and the settlement of the New World, using 5 Mb of noncoding DNA resequenced in 68 individuals from 4 populations (YRI, CHB, CEU, and MXL) by the Environmental Genome Project. We infer divergence between West African and Eurasian populations 140 thousand years ago (95% confidence interval: 40-270 kya). This is earlier than other genetic studies, in part because we incorporate migration. We estimate the European (CEU) and East Asian (CHB) divergence time to be 23 kya (95% c.i.: 17-43 kya), long after archeological evidence places modern humans in Europe. Finally, we estimate divergence between East Asians (CHB) and Mexican-Americans (MXL) of 22 kya (95% c.i.: 16.3-26.9 kya), and our analysis yields no evidence for subsequent migration. Furthermore, combining our demographic model with a previously estimated distribution of selective effects among newly arising amino acid mutations accurately predicts the frequency spectrum of nonsynonymous variants across three continental populations (YRI, CHB, CEU).

Full text

Its licence does not allow us to show the full text. Read it at PubMed Central:

References 61

1. 1.

   Mellars P. Going east: new genetic and archaeological perspectives on the modern human colonization of Eurasia. Science. 2006;313:796-800.

   Not in our library10.1126/science.1128402

2. 2.

   Goebel T, Waters MR, O'Rourke DH. The late Pleistocene dispersal of modern humans in the Americas. Science. 2008;319:1497-1502.

   Not in our library10.1126/science.1153569

3. 3.

   Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG. Recent and ongoing selection in the human genome. Nat Rev Genet. 2007;8:857-868.

   Not in our library10.1038/nrg2187

4. 4.

   Adams AM, Hudson RR. Maximum-likelihood estimation of demographic parameters using the frequency spectrum of unlinked single-nucleotide polymorphisms. Genetics. 2004;168:1699-1712.

   Not in our library10.1534/genetics.104.030171

5. 5.

   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.

   Not in our library10.1534/genetics.166.1.351

6. 6.

   Voight BF, Adams AM, Frisse LA, Qian Y, Hudson RR, et al. Interrogating multiple aspects of variation in a full resequencing data set to infer human population size changes. Proc Natl Acad Sci USA. 2005;102:18508-18513.

   Not in our library10.1073/pnas.0507325102

7. 7.
   On the Number of New World Founders: A Population Genetic Portrait of the Peopling of the Americas · 2005

   No notice
8. 8.

   Schaffner SF, Foo C, Gabriel S, Reich D, Daly MJ, et al. Calibrating a coalescent simulation of human genome sequence variation. Genome Res. 2005;15:1576-1583.

   Not in our library10.1101/gr.3709305

9. 9.

   Becquet C, Przeworski M. A new approach to estimate parameters of speciation models with application to apes. Genome Res. 2007;17:1505-1519.

   Not in our library10.1101/gr.6409707

10. 10.
    Genome-Wide Patterns of Nucleotide Polymorphism in Domesticated Rice · 2007

    No notice
11. 11.

    Keinan A, Mullikin JC, Patterson N, Reich D. Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in Europeans. Nat Genet. 2007;39:1251-1255.

    Not in our library10.1038/ng2116

12. 12.

    Garrigan D, Kingan SB, Pilkington MM, Wilder JA, Cox MP, et al. Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data. Genetics. 2007;177:2195-2207.

    Not in our library10.1534/genetics.107.077495

13. 13.

    Mulligan CJ, Kitchen A, Miyamoto MM. Updated three-stage model for the peopling of the Americas. PLoS ONE. 2008;3:e3199.

    Not in our library10.1371/journal.pone.0003199

14. 14.
    A three-stage colonization model for the peopling of the Americas. · 2008

    Not checked
15. 15.

    Cox M, Woerner A, Wall J, Hammer M. Intergenic DNA sequences from the human X chromosome reveal high rates of global gene flow. BMC Genetics. 2008;9:76.

    Not in our library10.1186/1471-2156-9-76

16. 16.

    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.

    Not in our library10.1093/molbev/msi103

17. 17.

    Hernandez RD, Williamson SH, Bustamante CD. Context dependence, ancestral misidentification, and spurious signatures of natural selection. Mol Biol Evol. 2007;24:1792-1800.

    Not in our library10.1093/molbev/msm108

18. 18.

    Nielsen R, Hubisz M, Hellmann I, Torgerson D, Andrés A, et al. Darwinian and demographic forces affecting human protein coding genes. 2009.

    Not in our library10.1101/gr.088336.108

19. 19.

    Hey J, Nielsen R. Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics. 2004;167:747-760.

    Not in our library10.1534/genetics.103.024182

20. 20.

    Sawyer SA, Hartl DL. Population genetics of polymorphism and divergence. Genetics. 1992;132:1161-1176.

    Not in our library10.1093/genetics/132.4.1161

21. 21.

    Bustamante CD, Wakeley J, Sawyer S, Hartl DL. Directional selection and the site-frequency spectrum. Genetics. 2001;159:1779-1788.

    Not in our library10.1093/genetics/159.4.1779

22. 22.

    Wakeley J. Coalescent Theory: an Introduction. 2008.

    Not in our library

23. 23.

    Williamson SH, Hernandez R, Fledel-Alon A, Zhu L, Nielsen R, et al. Simultaneous inference of selection and population growth from patterns of variation in the human genome. Proc Natl Acad Sci USA. 2005;102:7882-7887.

    Not in our library10.1073/pnas.0502300102

24. 24.

    Hernandez RD, Hubisz MJ, Wheeler DA, Smith DG, Ferguson B, et al. Demographic histories and patterns of linkage disequilibrium in Chinese and Indian rhesus macaques. Science. 2007;316:240-243.

    Not in our library10.1126/science.1140462

25. 25.

    Wiuf C. Consistency of estimators of population scaled parameters using composite likelihood. J Math Biol. 2006;53:821-841.

    Not in our library10.1007/s00285-006-0031-0

26. 26.

    Zhu L, Bustamante CD. A composite-likelihood approach for detecting directional selection from DNA sequence data. Genetics. 2005;170:1411-1421.

    Not in our library10.1534/genetics.104.035097

27. 27.

    Livingston RJ, von Niederhausern A, Jegga AG, Crawford DC, Carlson CS, et al. Pattern of sequence variation across 213 environmental response genes. Genome Res. 2004;14:1821-1831.

    Not in our library10.1101/gr.2730004

28. 28.

    Fischer RA. On the dominance ratio. Proc Roy Soc Edin. 1922;55:399-433.

    Not in our library

29. 29.

    Kimura M. Diffusion models in population genetics. J Appl Probab. 1964;1:177-232.

    Not in our library

30. 30.

    Ewens WJ. Mathematical Population Genetics: I. Theoretical Introduction. 2000.

    Not in our library

31. 31.

    Watterson GA. On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975;7:256-276.

    Not in our library10.1016/0040-5809(75)90020-9

32. 32.

    Nagylaki T. The strong-migration limit in geographically structured populations. J Math Biol. 1980;9:101-114.

    Not in our library10.1007/bf00275916

33. 33.

    Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R. Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005;15:1496-1502.

    Not in our library10.1101/gr.4107905

34. 34.

    Nielsen R, Hubisz MJ, Clark AG. Reconstituting the frequency spectrum of ascertained single-nucleotide polymorphism data. Genetics. 2004;168:2373-2382.

    Not in our library10.1534/genetics.104.031039

35. 35.

    Hudson RR. Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics. 2002;18:337-338.

    Not in our library10.1093/bioinformatics/18.2.337

36. 36.

    Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes: The Art of Scientific Computing. 2007.

    Not in our library

37. 37.

    Chang JS, Cooper G. A practical difference scheme for Fokker-Planck equations. J Comput Phys. 1970;6:1-16.

    Not in our library

38. 38.

    Oliphant TE. Guide to NumPy. 2006.

    Not in our library

39. 39.

    Oliphant TE. Python for scientific computing. Comput Sci Eng. 2007;9:10-20.

    Not in our library

40. 40.

    Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9:90-95.

    Not in our library

41. 41.

    NIEHS Environmental Genome Project.

    Not in our library

42. 42.
    Population History and Natural Selection Shape Patterns of Genetic Variation in 132 Genes · 2004

    No notice
43. 43.

    Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005;437:69-87.

    Not in our library10.1038/nature04072

44. 44.

    Hwang DG, Green P. Bayesian Markov chain Monte Carlo sequence analysis reveals varying neutral substitution patterns in mammalian evolution. Proc Natl Acad Sci USA. 2004;101:13994-14001.

    Not in our library10.1073/pnas.0404142101

45. 45.

    Kumar S, Filipski A, Swarna V, Walker A, Hedges SB. Placing confidence limits on the molecular age of the human-chimpanzee divergence. Proc Natl Acad Sci USA. 2005;102:18842-18847.

    Not in our library10.1073/pnas.0509585102

46. 46.

    Kondrashov AS. Direct estimates of human per nucleotide mutation rates at 20 loci causing Mendelian diseases. Hum Mutat. 2002;21:12-27.

    Not in our library10.1002/humu.10147

47. 47.

    Fenner JN. Cross-cultural estimation of the human generation interval for use in genetics-based population divergence studies. Am J Phys Anthropol. 2005;128:415-423.

    Not in our library10.1002/ajpa.20188

48. 48.

    Tremblay M, Vézina H. New estimates of intergenerational time intervals for the calculation of age and origin of mutations. Am J Hum Genet. 2000;66:651-658.

    Not in our library10.1086/302770

49. 49.
    Assessing the Evolutionary Impact of Amino Acid Mutations in the Human Genome · 2008

    No notice
50. 50.

    Heinrich JG. Can the likelihood-function value be used to measure goodness of fit? 2001.

    Not in our library

51. 51.

    Price AL, Patterson N, Yu F, Cox DR, Waliszewska A, et al. A genomewide admixture map for Latino populations. Am J Hum Genet. 2007;80:1024-1036.

    Not in our library10.1086/518313

52. 52.

    Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945-959.

    Not in our library10.1093/genetics/155.2.945

53. 53.

    Patterson N, Hattangadi N, Lane B, Lohmueller KE, Hafler DA, et al. Methods for high-density admixture mapping of disease genes. Am J Hum Genet. 2004;74:979-1000.

    Not in our library10.1086/420871

54. 54.

    Kryukov GV, Shpunt A, Stamatoyannopoulos JA, Sunyaev SR. Power of deep, all-exon resequencing for discovery of human trait genes. Proc Natl Acad Sci USA. 2009;106:3871-3876.

    Not in our library10.1073/pnas.0812824106

55. 55.

    Weale ME, Weiss DA, Jager RF, Bradman N, Thomas MG. Y chromosome evidence for Anglo-Saxon mass migration. Mol Biol Evol. 2002;19:1008-1021.

    Not in our library10.1093/oxfordjournals.molbev.a004160

56. 56.

    Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, et al. Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008;319:1100-1104.

    Not in our library10.1126/science.1153717

57. 57.

    Jakobsson M, Scholz SW, Scheet P, Gibbs JR, VanLiere JM, et al. Genotype, haplotype and copy-number variation in worldwide human populations. Nature. 2008;451:998-1003.

    Not in our library10.1038/nature06742

58. 58.

    Wall JD, Cox MP, Mendez FL, Woerner A, Severson T, et al. A novel DNA sequence database for analyzing human demographic history. Genome Res. 2008;18:1354-1361.

    Not in our library10.1101/gr.075630.107

59. 59.

    Braverman JM, Hudson RR, Kaplan NL, Langley CH, Stephan W. The hitchhiking effect on the site frequency spectrum of DNA polymorphisms. Genetics. 1995;140:783-796.

    Not in our library10.1093/genetics/140.2.783

60. 60.

    Myers S, Fefferman C, Patterson N. Can one learn history from the allelic spectrum? Theor Popul Biol. 2008;73:342-348.

    Not in our library10.1016/j.tpb.2008.01.001

61. 61.

    Pierce DA, Schafer DW. Residuals in generalized linear models. J Am Stat Assoc. 1986;81:977-986.

    Not in our library

Cited by papers we hold 55

• Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef · 2023
  Not checked
• On the limits of fitting complex models of population history to f-statistics · 2023
  No notice
• Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes · 2022
  Not checked
• Genotype–environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence · 2022
  Not checked
• Recommendations for improving statistical inference in population genomics · 2022
  No notice
• Demographic inference from multiple whole genomes using a particle filter for continuous Markov jump processes. · 2021
  Not checked
• Evidence for opposing selective forces operating on human-specific duplicated TCAF genes in Neanderthals and humans · 2021
  Not checked
• Detection of hard and soft selective sweeps from Drosophila melanogaster population genomic data. · 2021
  Not checked
• The Impact of Purifying and Background Selection on the Inference of Population History: Problems and Prospects · 2021
  No notice
• Inclusion of variants discovered from diverse populations improves polygenic risk score transferability · 2021
  No notice
• Evaluating insect-host interactions as a driver of species divergence in palm flower weevils · 2020
  Not checked
• Dating genomic variants and shared ancestry in population-scale sequencing data · 2020
  No notice
• A community-maintained standard library of population genetic models · 2020
  No notice
• Contrasting signatures of genomic divergence during sympatric speciation · 2020
  No notice
• Stairway Plot 2: demographic history inference with folded SNP frequency spectra · 2020
  Corrected
• Origin and adaptation to high altitude of Tibetan semi-wild wheat · 2020
  No notice
• An estimator of first coalescent time reveals selection on young variants and large heterogeneity in rare allele ages among human populations. · 2019
  Not checked
• A practical introduction to sequentially Markovian coalescent methods for estimating demographic history from genomic data · 2019
  No notice
• Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history · 2019
  No notice
• Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement · 2019
  No notice

and 35 more

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.