Methods for detecting the genomic signatures of natural selection have been heavily studied, and they have been successful in identifying many selective sweeps. For most of these sweeps, the favored allele remains unknown, making it difficult to distinguish carriers of the sweep from non-carriers. In an ongoing selective sweep, carriers of the favored allele are likely to contain a future most recent common ancestor. Therefore, identifying them may prove useful in predicting the evolutionary trajectory—for example, in contexts involving drug-resistant pathogen strains or cancer subclones. The main contribution of this paper is the development and analysis of a new statistic, the Haplotype Allele Frequency (HAF) score. The HAF score, assigned to individual haplotypes in a sample, naturally captures many of the properties shared by haplotypes carrying a favored allele. We provide a theoretical framework for computing expected HAF scores under different evolutionary scenarios, and we validate the theoretical predictions with simulations. As an application of HAF score computations, we develop an algorithm (PreCIOSS: Predicting Carriers of Ongoing Selective Sweeps) to identify carriers of the favored allele in selective sweeps, and we demonstrate its power on simulations of both hard and soft sweeps, as well as on data from well-known sweeps in human populations.
References
[1]
Fu W, Akey JM. Selection and adaptation in the human genome. Annu Rev Genomics Hum Genet. 2013;14:467–489. doi: 10.1146/annurev-genom-091212-153509. pmid:23834317
[2]
Lachance J, Tishkoff SA. Population Genomics of Human Adaptation. Annu Rev Ecol Evol Syst. 2013 Nov;44:123–143. doi: 10.1146/annurev-ecolsys-110512-135833. pmid:25383060
[3]
Vitti JJ, Grossman SR, Sabeti PC. Detecting natural selection in genomic data. Annu Rev Genet. 2013;47:97–120. doi: 10.1146/annurev-genet-111212-133526. pmid:24274750
[4]
Nielsen R, Williamson S, Kim Y, Hubisz MJ, Clark AG, Bustamante C. Genomic scans for selective sweeps using SNP data. Genome Res. 2005 Nov;15(11):1566–1575. doi: 10.1101/gr.4252305. pmid:16251466
[5]
Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009 May;19(5):826–837. doi: 10.1101/gr.087577.108. pmid:19307593
[6]
Chen H, Patterson N, Reich D. Population differentiation as a test for selective sweeps. Genome Res. 2010 Mar;20(3):393–402. doi: 10.1101/gr.100545.109. pmid:20086244
[7]
Berg JJ, Coop G. A population genetic signal of polygenic adaptation. PLoS Genet. 2014 Aug;10(8):e1004412. doi: 10.1371/journal.pgen.1004412. pmid:25102153
[8]
Jeong C, Di Rienzo A. Adaptations to local environments in modern human populations. Curr Opin Genet Dev. 2014 Dec;29C:1–8. doi: 10.1016/j.gde.2014.06.011.
[9]
Tekola-Ayele F, Adeyemo A, Chen G, Hailu E, Aseffa A, Davey G, et al. Novel genomic signals of recent selection in an Ethiopian population. Eur J Hum Genet. 2014 Nov; advance online publication. doi: 10.1038/ejhg.2014.233. pmid:25370040
[10]
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, et al. Sequencing of 50 Human Exomes Reveals Adaptation to High Altitude. Science. 2010;329(5987):75–78. Available from: . doi: 10.1126/science.1190371. pmid:20595611
[11]
Simonson TS, Yang Y, Huff CD, Yun H, Qin G, Witherspoon DJ, et al. Genetic evidence for high-altitude adaptation in Tibet. Science. 2010 Jul;329(5987):72–75. doi: 10.1126/science.1189406. pmid:20466884
[12]
Scheinfeldt LB, Soi S, Thompson S, Ranciaro A, Woldemeskel D, Beggs W, et al. Genetic adaptation to high altitude in the Ethiopian highlands. Genome Biol. 2012;13(1):R1. doi: 10.1186/gb-2012-13-1-r1. pmid:22264333
[13]
Alkorta-Aranburu G, Beall CM, Witonsky DB, Gebremedhin A, Pritchard JK, Di Rienzo A. The genetic architecture of adaptations to high altitude in Ethiopia. PLoS Genet. 2012;8(12):e1003110. doi: 10.1371/journal.pgen.1003110. pmid:23236293
[14]
Huerta-Sanchez E, Degiorgio M, Pagani L, Tarekegn A, Ekong R, Antao T, et al. Genetic signatures reveal high-altitude adaptation in a set of ethiopian populations. Mol Biol Evol. 2013 Aug;30(8):1877–1888. doi: 10.1093/molbev/mst089. pmid:23666210
[15]
Udpa N, Ronen R, Zhou D, Liang J, Stobdan T, Appenzeller O, et al. Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes. Genome Biol. 2014 Feb;15(2):R36. doi: 10.1186/gb-2014-15-2-r36. pmid:24555826
[16]
Zhou D, Udpa N, Ronen R, Stobdan T, Liang J, Appenzeller O, et al. Whole-genome sequencing uncovers the genetic basis of chronic mountain sickness in Andean highlanders. Am J Hum Genet. 2013 Sep;93(3):452–462. doi: 10.1016/j.ajhg.2013.07.011. pmid:23954164
Pavlidis P, Jensen JD, Stephan W. Searching for footprints of positive selection in whole-genome SNP data from nonequilibrium populations. Genetics. 2010 Jul;185(3):907–922. doi: 10.1534/genetics.110.116459. pmid:20407129
[22]
Lin K, Li H, Schlotterer C, Futschik A. Distinguishing positive selection from neutral evolution: boosting the performance of summary statistics. Genetics. 2011 Jan;187(1):229–244. doi: 10.1534/genetics.110.122614. pmid:21041556
[23]
Ronen R, Udpa N, Halperin E, Bafna V. Learning natural selection from the site frequency spectrum. Genetics. 2013 Sep;195(1):181–193. doi: 10.1534/genetics.113.152587. pmid:23770700
[24]
Simonsen KL, Churchill GA, Aquadro CF. Properties of statistical tests of neutrality for DNA polymorphism data. Genetics. 1995 Sep;141(1):413–429. pmid:8536987
[25]
Braverman JM, Hudson RR, Kaplan NL, Langley CH, Stephan W. The hitchhiking effect on the site frequency spectrum of DNA polymorphisms. Genetics. 1995 Jun;140(2):783–796. pmid:7498754
[26]
Hudson RR, Bailey K, Skarecky D, Kwiatowski J, Ayala FJ. Evidence for positive selection in the superoxide dismutase (Sod) region of Drosophila melanogaster. Genetics. 1994 Apr;136(4):1329–1340. pmid:8013910
[27]
Depaulis F, Mousset S, Veuille M. Haplotype tests using coalescent simulations conditional on the number of segregating sites. Mol Biol Evol. 2001 Jun;18(6):1136–1138. doi: 10.1093/oxfordjournals.molbev.a003885. pmid:11371602
[28]
Innan H, Zhang K, Marjoram P, Tavare S, Rosenberg NA. Statistical tests of the coalescent model based on the haplotype frequency distribution and the number of segregating sites. Genetics. 2005 Mar;169(3):1763–1777. doi: 10.1534/genetics.104.032219. pmid:15654103
[29]
Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002 Oct;419(6909):832–837. doi: 10.1038/nature01140. pmid:12397357
[30]
Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006 Mar;4(3):e72. doi: 10.1371/journal.pbio.0040072. pmid:16494531
[31]
Toomajian C, Hu TT, Aranzana MJ, Lister C, Tang C, Zheng H, et al. A nonparametric test reveals selection for rapid flowering in the Arabidopsis genome. PLoS Biol. 2006 May;4(5):e137. doi: 10.1371/journal.pbio.0040137. pmid:16623598
[32]
Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007 Oct;449(7164):913–918. doi: 10.1038/nature06250. pmid:17943131
[33]
Kim Y, Stephan W. Selective sweeps in the presence of interference among partially linked loci. Genetics. 2003 May;164(1):389–398. pmid:12750349
[34]
Messer PW, Petrov DA. Population genomics of rapid adaptation by soft selective sweeps. Trends Ecol Evol (Amst). 2013 Nov;28(11):659–669. doi: 10.1016/j.tree.2013.08.003.
[35]
Hermisson J, Pennings PS. Soft sweeps: molecular population genetics of adaptation from standing genetic variation. Genetics. 2005 Apr;169(4):2335–2352. doi: 10.1534/genetics.104.036947. pmid:15716498
[36]
Pennings PS, Hermisson J. Soft sweeps II–molecular population genetics of adaptation from recurrent mutation or migration. Mol Biol Evol. 2006 May;23(5):1076–1084. doi: 10.1093/molbev/msj117. pmid:16520336
[37]
Ferrer-Admetlla A, Liang M, Korneliussen T, Nielsen R. On detecting incomplete soft or hard selective sweeps using haplotype structure. Mol Biol Evol. 2014 May;31(5):1275–1291. doi: 10.1093/molbev/msu077. pmid:24554778
[38]
Garud NR, Messer PW, Buzbas EO, Petrov DA. Recent selective sweeps in North American Drosophila melanogaster show signatures of soft sweeps. PLoS Genet. 2015 Feb;11(2):e1005004. doi: 10.1371/journal.pgen.1005004. pmid:25706129
[39]
Peter BM, Huerta-Sanchez E, Nielsen R. Distinguishing between selective sweeps from standing variation and from a de novo mutation. PLoS Genet. 2012;8(10):e1003011. doi: 10.1371/journal.pgen.1003011. pmid:23071458
[40]
Schrider DR, Mendes FK, Hahn MW, Kern AD. Soft Shoulders Ahead: Spurious Signatures of Soft and Partial Selective Sweeps Result from Linked Hard Sweeps. Genetics. 2015 Feb; advance online publication.
[41]
Wilson BA, Petrov DA, Messer PW. Soft selective sweeps in complex demographic scenarios. Genetics. 2014 Oct;198(2):669–684. doi: 10.1534/genetics.114.165571. pmid:25060100
[42]
Fu YX. Statistical properties of segregating sites. Theor Popul Biol. 1995 Oct;48(2):172–197. doi: 10.1006/tpbi.1995.1025. pmid:7482370
[43]
Hudson RR. Gene genealogies and the coalescent process. In: Futuyma D, Antonovics J, editors. Oxford Surveys in Evolutionary Biology. Oxford: Oxford University Press; 1990. p. 1–44.
[44]
Slatkin M, Hudson RR. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991 Oct;129(2):555–562. pmid:1743491
[45]
Graham R, Knuth DE, Patashnik O. Concrete Mathematics: A Foundation for Computer Science. 2nd ed. Reading, Mass: Addison-Wesley; 1994.
[46]
Nordborg M. Coalescent Theory. In: Balding DJ, Bishop M, Cannings C, editors. Handbook of statistical genetics. 3rd ed. John Wiley & Sons, Ltd; 2008. p. 843–877.
[47]
Ewing G, Hermisson J. MSMS: a coalescent simulation program including recombination, demographic structure and selection at a single locus. Bioinformatics. 2010 Aug;26(16):2064–2065. doi: 10.1093/bioinformatics/btq322. pmid:20591904
[48]
Brodersen KH, Ong CS, Stephan KE, Buhmann JM. The Balanced Accuracy and Its Posterior Distribution. In: Pattern Recognition (ICPR), 2010 20th International Conference on; 2010. p. 3121–3124.
[49]
Grossman SR, Shlyakhter I, Shylakhter I, Karlsson EK, Byrne EH, Morales S, et al. A composite of multiple signals distinguishes causal variants in regions of positive selection. Science. 2010 Feb;327(5967):883–886. doi: 10.1126/science.1183863. pmid:20056855
[50]
Gravel S, Henn BM, Gutenkunst RN, Indap AR, Marth GT, Clark AG, et al. Demographic history and rare allele sharing among human populations. Proc Natl Acad Sci USA. 2011 Jul;108(29):11983–11988. doi: 10.1073/pnas.1019276108. pmid:21730125
[51]
Altshuler DM, et al. Integrating common and rare genetic variation in diverse human populations. Nature. 2010 Sep;467(7311):52–58. doi: 10.1038/nature09298. pmid:20811451
[52]
Sequencing TC, Consortium A. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005 Sep;437(7055):69–87. doi: 10.1038/nature04072.
[53]
Kuokkanen M, Enattah NS, Oksanen A, Savilahti E, Orpana A, Jarvela I. Transcriptional regulation of the lactase-phlorizin hydrolase gene by polymorphisms associated with adult-type hypolactasia. Gut. 2003 May;52(5):647–652. doi: 10.1136/gut.52.5.647. pmid:12692047
[54]
Olds LC, Sibley E. Lactase persistence DNA variant enhances lactase promoter activity in vitro: functional role as a cis regulatory element. Hum Mol Genet. 2003 Sep;12(18):2333–2340. doi: 10.1093/hmg/ddg244. pmid:12915462
[55]
Troelsen JT, Olsen J, M？ller J, Sj？str？m H. An upstream polymorphism associated with lactase persistence has increased enhancer activity. Gastroenterology. 2003 Dec;125(6):1686–1694. doi: 10.1053/j.gastro.2003.09.031. pmid:14724821
[56]
Akey JM, Eberle MA, Rieder MJ, Carlson CS, Shriver MD, Nickerson DA, et al. Population history and natural selection shape patterns of genetic variation in 132 genes. PLoS Biol. 2004 Oct;2(10):e286. doi: 10.1371/journal.pbio.0020286. pmid:15361935
[57]
Stajich JE, Hahn MW. Disentangling the effects of demography and selection in human history. Mol Biol Evol. 2005 Jan;22(1):63–73. doi: 10.1093/molbev/msh252. pmid:15356276
[58]
Akey JM, Swanson WJ, Madeoy J, Eberle M, Shriver MD. TRPV6 exhibits unusual patterns of polymorphism and divergence in worldwide populations. Hum Mol Genet. 2006 Jul;15(13):2106–2113. doi: 10.1093/hmg/ddl134. pmid:16717058
[59]
Bhatia G, Patterson N, Pasaniuc B, Zaitlen N, Genovese G, Pollack S, et al. Genome-wide comparison of African-ancestry populations from CARe and other cohorts reveals signals of natural selection. Am J Hum Genet. 2011 Sep;89(3):368–381. doi: 10.1016/j.ajhg.2011.07.025. pmid:21907010
[60]
Sakamoto H, Yoshimura K, Saeki N, Katai H, Shimoda T, Matsuno Y, et al. Genetic variation in PSCA is associated with susceptibility to diffuse-type gastric cancer. Nat Genet. 2008 Jun;40(6):730–740. doi: 10.1038/ng.152. pmid:18488030
[61]
Wu X, Ye Y, Kiemeney LA, Sulem P, Rafnar T, Matullo G, et al. Genetic variation in the prostate stem cell antigen gene PSCA confers susceptibility to urinary bladder cancer. Nat Genet. 2009 Sep;41(9):991–995. doi: 10.1038/ng.421. pmid:19648920
[62]
Whitfield JB. Alcohol dehydrogenase and alcohol dependence: variation in genotype-associated risk between populations. Am J Hum Genet. 2002 Nov;71(5):1247–1250. doi: 10.1086/344287. pmid:12452180
[63]
Peng Y, Shi H, Qi XB, Xiao CJ, Zhong H, Ma RL, et al. The ADH1B Arg47His polymorphism in east Asian populations and expansion of rice domestication in history. BMC Evol Biol. 2010;10:15. doi: 10.1186/1471-2148-10-15. pmid:20089146
[64]
Osier MV, Pakstis AJ, Soodyall H, Comas D, Goldman D, Odunsi A, et al. A global perspective on genetic variation at the ADH genes reveals unusual patterns of linkage disequilibrium and diversity. Am J Hum Genet. 2002 Jul;71(1):84–99. doi: 10.1086/341290. pmid:12050823
[65]
Eng MY, Luczak SE, Wall TL. ALDH2, ADH1B, and ADH1C genotypes in Asians: a literature review. Alcohol Res Health. 2007;30(1):22–27. pmid:17718397
[66]
Li H, Mukherjee N, Soundararajan U, Tarnok Z, Barta C, Khaliq S, et al. Geographically separate increases in the frequency of the derived ADH1B*47His allele in eastern and western Asia. Am J Hum Genet. 2007 Oct;81(4):842–846. doi: 10.1086/521201. pmid:17847010
[67]
McGovern PE, Zhang J, Tang J, Zhang Z, Hall GR, Moreau RA, et al. Fermented beverages of pre- and proto-historic China. Proc Natl Acad Sci USA. 2004 Dec;101(51):17593–17598. doi: 10.1073/pnas.0407921102. pmid:15590771
[68]
Fujimoto A, Ohashi J, Nishida N, Miyagawa T, Morishita Y, Tsunoda T, et al. A replication study confirmed the EDAR gene to be a major contributor to population differentiation regarding head hair thickness in Asia. Hum Genet. 2008 Sep;124(2):179–185. doi: 10.1007/s00439-008-0537-1. pmid:18704500
[69]
Kimura R, Yamaguchi T, Takeda M, Kondo O, Toma T, Haneji K, et al. A common variation in EDAR is a genetic determinant of shovel-shaped incisors. Am J Hum Genet. 2009 Oct;85(4):528–535. doi: 10.1016/j.ajhg.2009.09.006. pmid:19804850
[70]
Bryk J, Hardouin E, Pugach I, Hughes D, Strotmann R, Stoneking M, et al. Positive selection in East Asians for an EDAR allele that enhances NF-kappaB activation. PLoS ONE. 2008;3(5):e2209. doi: 10.1371/journal.pone.0002209. pmid:18493316
[71]
Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007 Oct;449(7164):913–918. doi: 10.1038/nature06250. pmid:17943131
[72]
Williamson SH, Hernandez R, Fledel-Alon A, Zhu L, Nielsen R, Bustamante CD. Simultaneous inference of selection and population growth from patterns of variation in the human genome. Proc Natl Acad Sci USA. 2005 May;102(22):7882–7887. doi: 10.1073/pnas.0502300102. pmid:15905331
[73]
Luksza M, Lassig M. A predictive fitness model for influenza. Nature. 2014 Mar;507(7490):57–61. doi: 10.1038/nature13087. pmid:24572367
[74]
Lee MC, Lopez-Diaz FJ, Khan SY, Tariq MA, Dayn Y, Vaske CJ, et al. Single-cell analyses of transcriptional heterogeneity during drug tolerance transition in cancer cells by RNA sequencing. Proc Natl Acad Sci USA. 2014 Nov;111(44):E4726–4735. doi: 10.1073/pnas.1404656111. pmid:25339441
[75]
Nachman MW, Crowell SL. Estimate of the mutation rate per nucleotide in humans. Genetics. 2000 Sep;156(1):297–304. pmid:10978293
[76]
Campbell CD, Chong JX, Malig M, Ko A, Dumont BL, Han L, et al. Estimating the human mutation rate using autozygosity in a founder population. Nat Genet. 2012 Nov;44(11):1277–1281. doi: 10.1038/ng.2418. pmid:23001126
[77]
Hey J, Wakeley J. A coalescent estimator of the population recombination rate. Genetics. 1997 Mar;145(3):833–846. pmid:9055092
[78]
Szpiech ZA, Hernandez RD. selscan: An Efficient Multithreaded Program to Perform EHH-Based Scans for Positive Selection. Mol Biol Evol. 2014 Oct;31(10):2824–2827. doi: 10.1093/molbev/msu211. pmid:25015648
[79]
Frazer KA, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007 Oct;449(7164):851–861. doi: 10.1038/nature06258. pmid:17943122
