title: Calibrating a coalescent simulation of human genome sequence variation.
abstract: Population genetic models play an important role in human genetic research, connecting empirical observations about sequence variation with hypotheses about underlying historical and biological causes. More specifically, models are used to compare empirical measures of sequence variation, linkage disequilibrium (LD), and selection to expectations under a "null" distribution. In the absence of detailed information about human demographic history, and about variation in mutation and recombination rates, simulations have of necessity used arbitrary models, usually simple ones. With the advent of large empirical data sets, it is now possible to calibrate population genetic models with genome-wide data, permitting for the first time the generation of data that are consistent with empirical data across a wide range of characteristics. We present here the first such calibrated model and show that, while still arbitrary, it successfully generates simulated data (for three populations) that closely resemble empirical data in allele frequency, linkage disequilibrium, and population differentiation. No assertion is made about the accuracy of the proposed historical and recombination model, but its ability to generate realistic data meets a long-standing need among geneticists. We anticipate that this model, for which software is publicly available, and others like it will have numerous applications in empirical studies of human genetics.
authors.affiliation: Program in Medical and Population Genetics, The Broad Institute, Cambridge, Massachusetts 02139, USA. sfs@broad.mit.edu
authors.forename: Stephen F
authors.initials: SF
authors.lastname: Schaffner
authors.forename: Catherine
authors.initials: C
authors.lastname: Foo
authors.forename: Stacey
authors.initials: S
authors.lastname: Gabriel
authors.forename: David
authors.initials: D
authors.lastname: Reich
authors.forename: Mark J
authors.initials: MJ
authors.lastname: Daly
authors.forename: David
authors.initials: D
authors.lastname: Altshuler
mesh_terms: D003198:Computer Simulation; D019143:Evolution, Molecular; D005787:Gene Frequency; D014644:Genetic Variation; D005828:Genetics, Population; D015894:Genome, Human / Q000235:genetics*; D023281:Genomics / Q000379:methods*; D006801:Humans; D015810:Linkage Disequilibrium; D008957:Models, Genetic