Inferring very recent population growth rate from population-scale sequencing data: Using a large-sample coalescent estimator

Molecular Biology and Evolution

Volumn 32, Issue 11, 2015, Pages 2996-3011

  Chen, Hua ^a   Hey, Jody ^b   Chen, Kun ^c  

^a Chinese Academy of Sciences   (China)

^b TEMPLE UNIVERSITY   (United States)

^c HARVARD MEDICAL SCHOOL   (United States)

Author keywords

coalescent; genetics diversity; large sample sequencing; population growth rate

Indexed keywords

AFRICAN; ARTICLE; COMPUTER SIMULATION; EUROPEAN; EXOME; GENE SEQUENCE; GENEALOGY; GROWTH CURVE; GROWTH RATE; HAPLOTYPE; HISTORY; HUMAN; MUTATION RATE; POPULATION GROWTH; POPULATION MODEL; POPULATION SIZE; SAMPLE; SAMPLE SIZE; TIME; ALGORITHM; BIOLOGICAL MODEL; COMPARATIVE STUDY; DNA SEQUENCE; GENETIC VARIATION; POPULATION DENSITY; POPULATION GENETICS; PROCEDURES; STATISTICAL MODEL;

ALGORITHMS; GENETIC VARIATION; GENETICS, POPULATION; HUMANS; MODELS, GENETIC; MODELS, STATISTICAL; POPULATION DENSITY; POPULATION GROWTH; SEQUENCE ANALYSIS, DNA;

EID: 84946145659     PISSN: 07374038     EISSN: 15371719     Source Type: Journal    
DOI: 10.1093/molbev/msv158     Document Type: Article

References (51)

1. Adams A, Hudson R. 2004. Maximum-likelihood estimation of demographic parameters using the frequency spectrum of unlinked single-nucleotide polymorphisms. Genetics 168(3):1699-1712.
2. Akaike H. 1974. A new look at the statistical model identification. IEEE Trans Automat Contr. 19:716-723.
3. Altshuler D, Lander E, Ambrogio L, Bloom T, Cibulskis K, Fennell T, Gabriel S, Jaffe D, Shefler E, Sougnez C, et al. 2010. A map of human genome variation from population scale sequencing. Nature 467(7319):1061-1073.
4. Bhaskar A, Clark AG, Song YS. 2014. Distortion of genealogical properties when the sample is very large. Proc Natl Acad Sci U S A. 111(6):2385-2390.
5. Chen H. 2012. The joint allele frequency spectrum of multiple populations: a coalescent theory approach. Theor Popul Biol. 81(2):179-195.
6. Chen H. 2013. Intercoalescence time distribution of incomplete genealogies in temporally varying populations, and applications in population genetic inference. Ann Hum Genet. 77(2):158-173.
7. Chen H, Chen K. 2013. Asymptotic distributions of coalescence times and ancestral lineage numbers for populations with temporally varying size. Genetics 194(3):721-736.
8. Chen H, Green RE, Peaeabo S, Slatkin M. 2007. The joint allele-frequency spectrum in closely related species. Genetics 177(1):387-398.
9. Coventry A, Bull-Otterson LM, Liu X, Clark AG, Maxwell TJ, Crosby J, Hixson JE, Rea TJ, Muzny DM, Lewis LR, et al. 2010. Deep resequencing reveals excess rare recent variants consistent with explosive population growth. Nat Commun. 1:131.
10. Drummond AJ, Rambaut A, Shapiro B, Pybus OG. 2005. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 22(5):1185-1192.
11. Efron B. 1982. The jackknife, the bootstrap, and other resampling plans. Vol. 38. Philadelphia: SIAM.
12. Evans S, Shvets Y, Slatkin M. 2007. Non-equilibrium theory of the allele frequency spectrum. Theor Popul Biol. 71(1):109-119.
13. Fu W, O?Connor TD, Jun G, Kang HM, Abecasis G, Leal SM, Gabriel S, Altshuler D, Shendure J, Nickerson DA, et al. 2013. Analysis of 6, 515 exomes reveals the recent origin of most human protein-coding variants. Nature 493:216-220.
14. Fu YX. 1995. Statistical properties of segregating sites. Theor Popul Biol. 48(2):172-197.
15. Fu YX. 2006. Exact coalescent for the Wright-Fisher model. Theor Popul Biol. 69(4):385-394.
16. Galassi M, Davies J, Theiler J, Gough B, Jungman G, Booth M, Rossi F. 2009. GNU Scientific Library Reference Manual. Bristol (United Kingdom): Network Theory.
17. Gazave E, Ma L, Chang D, Coventry A, Gao F, Muzny D, Boerwinkle E, Gibbs RA, Sing CF, Clark AG, et al. 2014. Neutral genomic regions refine models of recent rapid human population growth. Proc Natl Acad Sci U S A. 111(2):757-762.
18. Genomes Project Consortium, et al. 2010. A map of human genome variation from population-scale sequencing. Nature 467(7319): 1061-1073.
19. Gravel S, Henn B, Gutenkunst R, Indap A, Marth G, Clark A, Yu F, Gibbs R, The 1000 Genomes Project, Bustamante C. 2011. Demographic history and rare allele sharing among human populations. Proc Natl Acad Sci U S A. 108(29):11983-11988.
20. Griffiths RC, Tavar-e S. 1994. Sampling theory for neutral alleles in a varying environment. Philos Trans R Soc Lond B Biol Sci. 344:403-410.
21. Griffiths RC, Tavar-e S. 1998. The age of a mutation in a general coalescent tree. Stoch Models. 14(1-2):273-295.
22. Gutenkunst RN, Hernandez RD, Williamson SH, Bustamante CD. 2009. Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data. PLoS Genet. 5(10):e1000695.
23. Harris K, Nielsen R. 2013. Inferring demographic history from a spectrum of shared haplotype lengths. PLoS Genet. 9(6):e1003521.
24. Hudson RR. 1990. Gene genealogies and the coalescent process. In: Futuyma DJ, Antonovics JD, editors. Oxford surveys in evolutionary biology. Vol. 7. New York: Oxford University Press. p. 1-44.
25. Keinan A, Clark AG. 2012. Recent explosive human population growth has resulted in an excess of rare genetic variants. Science 336(6082):740-743.
26. Keinan A, Mullikin JC, Patterson N, Reich D. 2007. Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in Europeans. Nat Genet. 39(10):1251-1255.
27. Li H. 2011. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27(21):2987-2993.
28. Li H, Durbin R. 2011. Inference of human population history from individual whole-genome sequences. Nature 475(7357):493-496.
29. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. 2009. The sequence alignment/map format and samtools. Bioinformatics 25(16):2078-2079.
30. Lukic S, Hey J. 2012. Demographic inference using spectral methods on SNP data, with an analysis of the human out-of-Africa expansion. Genetics 192(2):619-639.
31. Lukic S, Hey J, Chen K. 2011. Non-equilibrium allele frequency spectra via spectral methods. Theor Popul Biol. 79(4):203-219.
32. Marth GT, Czabarka E, Murvai J, Sherry ST. 2004. 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.
33. Maruvka Y, Shnerb N, Bar-Yam Y, Wakeley J. 2011. Recovering population parameters from a single gene genealogy: an unbiased estimator of the growth rate. Mol Biol Evol. 28(5):1617-1631.
34. McCullagh P, Nelder JA. 1989. Generalized linear models. London: CRC Press.
35. Nei M, Chakraborty R, Fuerst PA. 1976. Infinite allele model with varying mutation rate. Proc Natl Acad Sci U S A. 73(11):4164-4168.
36. Nelson M, Wegmann D, Ehm M, Kessner D, Jean P, Verzilli C, Shen J, Tang Z, Bacanu S, Fraser D, et al. 2012. An abundance of rare functional variants in 202 drug target genes sequenced in 14, 002 people. Science 337(6090):100-104.
37. Nielsen R, Korneliussen T, Albrechtsen A, Li Y, Wang J. 2012. SNP calling, genotype calling, and sample allele frequency estimation from newgeneration sequencing data. PLoS One 7(7):e37558.
38. Pan W. 2001. Akaike?s information criterion in generalized estimating equations. Biometrics 57(1):120-125.
39. Polanski A, Bobrowski A, Kimmel M. 2003. A note on distributions of times to coalescence, under time-dependent population size. Theor Popul Biol. 63(1):33-40.
40. Polanski A, Kimmel M. 2003. New explicit expressions for relative frequencies of single-nucleotide polymorphisms with application to statistical inference on population growth. Genetics 165(1):427-436.
41. Press W, Flannery B, Teukolsky S, Vetterling W. 1992. Numerical recipes in C: the art of scientific programming. Section 10:408-412.
42. Pybus OG, Rambaut A, Harvey PH. 2000. An integrated framework for the inference of viral population history from reconstructed genealogies. Genetics 155(3):1429-1437.
43. Schaffner SF, Foo C, Gabriel S, Reich D, Daly MJ, Altshuler D. 2005. Calibrating a coalescent simulation of human genome sequence variation. Genome Res. 15:1576-1583.
44. Song YS, Steinreucken M. 2012. A simple method for finding explicit analytic transition densities of diffusion processes with general diploid selection. Genetics 190(3):1117-1129.
45. Tennessen JA, Bigham AW, O?Connor TD, Fu W, Kenny EE, Gravel S, McGee S, Do R, Liu X, Jun G, et al. 2012. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science 337(6090):64-69.
46. Tyekucheva S, Makova KD, Karro JE, Hardison RC, Miller W, Chiaromonte F, et al. 2008. Human-macaque comparisons illuminate variation in neutral substitution rates. Genome Biol. 9(4):R76.
47. Wakeley J. 1993. Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA. J Mol Evol. 37(6):613-623.
48. Wakeley J, Takahashi T. 2003. Gene genealogies when the sample size exceeds the effective size of the population. Mol Biol Evol. 20(2):208-213.
49. Watterson G. 1975. On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 7(2):256-276.
50. Yang Z. 1994. Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J Mol Evol. 39(3):306-314.
51. Zivkovic D, Stephan W. 2011. Analytical results on the neutral nonequilibrium allele frequency spectrum based on diffusion theory. Theor Popul Biol. 79(4):184-191.