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. As applications, 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 also combine our demographic model with a previously estimated distribution of selective effects among newly arising amino acid mutations to accurately predict the frequency spectrum of nonsynonymous variants across three continental populations (YRI, CHB, CEU).

17 pages, 4 figures, supporting information included with source

Cited by in corpus (42)

• Inference of population splits and mixtures from genome-wide allele frequency data
• The genetic prehistory of southern Africa
• Population genetics models of local ancestry
• Sequence Capture Versus Restriction Site Associated DNA Sequencing for Phylogeography
• WFABC: a Wright-Fisher ABC-based approach for inferring effective population sizes and selection coefficients from time-sampled data
• Spatial Population Genetics: It's About Time
• Fundamental limits on the accuracy of demographic inference based on the sample frequency spectrum
• Neutral genomic regions refine models of recent rapid human population growth
• Inference of Natural Selection from Interspersed Genomic Elements Based on Polymorphism and Divergence
• Descartes' rule of signs and the identifiability of population demographic models from genomic variation data
• A novel spectral method for inferring general diploid selection from time series genetic data
• Distortion of genealogical properties when the sample is very large
• Efficient computation of the joint sample frequency spectra for multiple populations
• Population Genetics of Rare Variants and Complex Diseases
• A sequentially Markov conditional sampling distribution for structured populations with migration and recombination
• Exact simulation of the Wright-Fisher diffusion
• The site frequency spectrum for general coalescents
• An explicit transition density expansion for a multi-allelic Wright-Fisher diffusion with general diploid selection
• General triallelic frequency spectrum under demographic models with variable population size
• Separation of the largest eigenvalues in eigenanalysis of genotype data from discrete subpopulations
• A path integral formulation of the Wright-Fisher process with genic selection
• Computing the joint distribution of the total tree length across loci in populations with variable size
• Secure numerical simulations using fully homomorphic encryption
• Progress and open problems in evolutionary dynamics
• Demography and the age of rare variants
• The expected neutral frequency spectrum of linked sites
• The impact of population demography and selection on the genetic architecture of complex traits
• Can one hear the shape of a population history?
• Prediction and replication from case-control sequencing studies using custom genotyping and additional sequencing
• Genomic identification of founding haplotypes reveals the history of the selfing species Capsella rubella
• Recent selective sweeps in North American Drosophila melanogaster show signatures of soft sweeps
• Reproductive isolation between phylogeographic lineages scales with divergence
• Reconstructing Native American Migrations from Whole-genome and Whole-exome Data
• Non-equilibrium allele frequency spectra via spectral methods
• A note on the distribution of admixture segment lengths and ancestry proportions under pulse and two-wave admixture models
• The Rescaled Polya Urn and the Wright-Fisher process with mutation
• Minimal-assumption inference from population-genomic data
• Predicting discovery rates of genomic features
• High burden of private mutations due to explosive human population growth and purifying selection
• De novo genomic analyses for non-model organisms: an evaluation of methods across a multi-species data set
• Error-prone polymerase activity causes multinucleotide mutations in humans
• Estimating heterozygosity from a low-coverage genome sequence, leveraging data from other individuals sequenced at the same sites