Likelihood-free inference of population structure and local adaptation in a Bayesian hierarchical model.
We address the problem of finding evidence of natural selection from genetic data, accounting for the confounding effects of demographic history. In the absence of natural selection, gene genealogies should all be sampled from the same underlying distribution, often approximated by a coalescent model. Selection at a particular locus will lead to a modified genealogy, and this motivates a number of recent approaches for detecting the effects of natural selection in the genome as "outliers" under some models. The demographic history of a population affects the sampling distribution of genealogies, and therefore the observed genotypes and the classification of outliers. Since we cannot see genealogies directly, we have to infer them from the observed data under some model of mutation and demography. Thus the accuracy of an outlier-based approach depends to a greater or a lesser extent on the uncertainty about the demographic and mutational model. A natural modeling framework for this type of problem is provided by Bayesian hierarchical models, in which parameters, such as mutation rates and selection coefficients, are allowed to vary across loci. It has proved quite difficult computationally to implement fully probabilistic genealogical models with complex demographies, and this has motivated the development of approximations such as approximate Bayesian computation (ABC). In ABC the data are compressed into summary statistics, and computation of the likelihood function is replaced by simulation of data under the model. In a hierarchical setting one may be interested both in hyperparameters and parameters, and there may be very many of the latter--for example, in a genetic model, these may be parameters describing each of many loci or populations. This poses a problem for ABC in that one then requires summary statistics for each locus, which, if used naively, leads to a consequent difficulty in conditional density estimation. We develop a general method for applying ABC to Bayesian hierarchical models, and we apply it to detect microsatellite loci influenced by local selection. We demonstrate using receiver operating characteristic (ROC) analysis that this approach has comparable performance to a full-likelihood method and outperforms it when mutation rates are variable across loci.
Bazin
Eric
E
School of Biological Sciences, University of Reading, Whiteknights, Reading RG6 6BX, United Kingdom.
Dawson
Kevin J
KJ
Beaumont
Mark A
MA
D001699:Biometry; D003710:Demography; D016678:Genome; D006801:Humans; D016013:Likelihood Functions; D018895:Microsatellite Repeats; D008954:Models, Biological; D015233:Models, Statistical; D044382:Population Groups; D011336:Probability; D012372:ROC Curve; D012641:Selection, Genetic
Heredity (Edinb). 1986 Dec;57 ( Pt 3):357-76
3804765
J R Soc Interface. 2009 Feb 6;6(31):187-202
19205079
Mol Ecol. 2004 Apr;13(4):969-80
15012769
Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15324-8
14663152
Genetics. 2008 Mar;178(3):1817-29
18245358
Stat Appl Genet Mol Biol. 2008;7(1):Article26
18764775
Genetics. 2009 Aug;182(4):1207-18
19506307
Proc Natl Acad Sci U S A. 2007 Nov 6;104(45):17614-9
17978179
J Am Stat Assoc. 2009 Mar 1;104(485):142-154
19623271
Genetics. 1973 May;74(1):175-95
4711903
PLoS Genet. 2007 Apr 20;3(4):e66
17447846
Genetics. 1931 Mar;16(2):97-159
17246615
Forensic Sci Int. 1994 Feb;64(2-3):125-40
8175083
Bioinformatics. 2004 Dec 12;20(18):3673-5
15271777
Mol Biol Evol. 1999 Dec;16(12):1791-8
10605120
Curr Biol. 2005 Mar 8;15(5):R159-60
15753023
Proc Natl Acad Sci U S A. 2009 Jun 30;106(26):10576-81
19525398
Bioinformatics. 2004 Oct 12;20(15):2485-7
15117750
Nature. 2007 Oct 18;449(7164):913-8
17943131
Genet Res. 1996 Apr;67(2):147-58
8801187
Genetics. 2001 Aug;158(4):1811-23
11514464
Curr Opin Genet Dev. 2006 Dec;16(6):624-9
17008091
Genetics. 2009 Apr;181(4):1507-19
19189952
Proc R Soc Lond B Biol Sci. 1966 Mar 22;164(995):362-79
4379525
Genetics. 2002 Dec;162(4):2025-35
12524368
Bioinformatics. 2005 Oct 15;21(20):3940-1
16096348
Genome Res. 2007 Oct;17(10):1505-19
17712021
Mol Ecol. 2008 Aug;17(15):3425-7
19160473
BMC Evol Biol. 2008 Nov 27;8:322
19038027
Evolution. 2002 Jan;56(1):154-66
11913661
Mol Ecol. 2006 Jan;15(1):209-23
16367841
Genet Res. 1997 Oct;70(2):155-74
9449192
Annu Rev Genet. 2002;36:721-50
12359738
Theor Popul Biol. 2003 May;63(3):221-30
12689793
Bioinformatics. 2002 Feb;18(2):337-8
11847089
Genetics. 1995 Jan;139(1):457-62
7705646
Genetics. 1989 Dec;123(4):887-99
2612899
Genome Res. 2006 Jun;16(6):702-12
16687733
Genome Res. 2005 Nov;15(11):1468-76
16251456
Genetics. 1996 Apr;142(4):1357-62
8846911
Comput Appl Biosci. 1997 Jun;13(3):235-8
9183526
Evolution. 1984 Nov;38(6):1358-1370
28563791
Theor Popul Biol. 1998 Apr;53(2):166-74
9615475
Proc Natl Acad Sci U S A. 2007 Feb 6;104(6):1760-5
17264216
Genetics. 2005 May;170(1):409-17
15716508
Genetics. 2008 Oct;180(2):977-93
18780740
Heredity (Edinb). 2009 Oct;103(4):285-98
19623208
Theor Popul Biol. 1983 Jun;23(3):300-13
6623407