title: Inference of site frequency spectra from high-throughput sequence data: quantification of selection on nonsynonymous and synonymous sites in humans.
abstract: Sequencing errors and random sampling of nucleotide types among sequencing reads at heterozygous sites present challenges for accurate, unbiased inference of single-nucleotide polymorphism genotypes from high-throughput sequence data. Here, we develop a maximum-likelihood approach to estimate the frequency distribution of the number of alleles in a sample of individuals (the site frequency spectrum), using high-throughput sequence data. Our method assumes binomial sampling of nucleotide types in heterozygotes and random sequencing error. By simulations, we show that close to unbiased estimates of the site frequency spectrum can be obtained if the error rate per base read does not exceed the population nucleotide diversity. We also show that these estimates are reasonably robust if errors are nonrandom. We then apply the method to infer site frequency spectra for zerofold degenerate, fourfold degenerate, and intronic sites of protein-coding genes using the low coverage human sequence data produced by the 1000 Genomes Project phase-one pilot. By fitting a model to the inferred site frequency spectra that estimates parameters of the distribution of fitness effects of new mutations, we find evidence for significant natural selection operating on fourfold sites. We also find that a model with variable effects of mutations at synonymous sites fits the data significantly better than a model with equal mutational effects. Under the variable effects model, we infer that 11% of synonymous mutations are subject to strong purifying selection.
authors: Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JT, UK. keightley.genetics2011@gmail.com
authors: Peter D
authors: PD
authors: Keightley
authors: Daniel L
authors: DL
authors: Halligan
references: Theor Popul Biol. 1975 Apr;7(2):256-76
references: 1145509
references: Nature. 2010 Oct 28;467(7319):1061-73
references: 20981092
references: Genome Res. 2009 Jul;19(7):1195-201
references: 19439516
references: Mol Biol Evol. 2010 Jan;27(1):177-92
references: 19759235
references: PLoS Genet. 2008 May 30;4(5):e1000083
references: 18516229
references: Genetics. 2009 May;182(1):295-301
references: 19293142
references: Genome Res. 2005 Oct;15(10):1373-8
references: 16204190
references: Mol Biol Evol. 1992 Jul;9(4):744-52
references: 1630310
references: Philos Trans R Soc Lond B Biol Sci. 2010 Apr 27;365(1544):1187-93
references: 20308093
references: Genetics. 2006 Jun;173(2):891-900
references: 16547091
references: Genetics. 1991 Oct;129(2):513-23
references: 1743489
references: Mol Ecol. 2010 Mar;19 Suppl 1:277-84
references: 20331786
references: PLoS Genet. 2009 Aug;5(8):e1000592
references: 19662163
references: Mol Biol Evol. 2008 Jan;25(1):199-206
references: 17981928
references: PLoS Comput Biol. 2006 Jan;2(1):e5
references: 16410828
references: Mol Biol Evol. 2008 Nov;25(11):2409-19
references: 18725384
references: Nat Rev Genet. 2006 Feb;7(2):98-108
references: 16418745
references: Genetics. 2007 Dec;177(4):2251-61
references: 18073430
references: Genome Res. 2011 Jun;21(6):952-60
references: 20980557
references: Genome Res. 2008 Jul;18(7):1020-9
references: 18411405
references: Genetics. 1983 Oct;105(2):437-60
references: 6628982
references: Science. 2010 Jan 1;327(5961):92-4
references: 20044577
references: Proc Natl Acad Sci U S A. 2006 May 2;103(18):6940-5
references: 16632609
references: Genome Res. 2010 Jan;20(1):101-9
references: 19952140