A general method for calculating likelihoods under the coalescent process

K. Lohse, R. Harrison, N.H. Barton, Genetics 189 (2011) 977–987.


Journal Article | Published | English
Author
; ;
Department
Abstract
Analysis of genomic data requires an efficient way to calculate likelihoods across very large numbers of loci. We describe a general method for finding the distribution of genealogies: we allow migration between demes, splitting of demes [as in the isolation-with-migration (IM) model], and recombination between linked loci. These processes are described by a set of linear recursions for the generating function of branch lengths. Under the infinite-sites model, the probability of any configuration of mutations can be found by differentiating this generating function. Such calculations are feasible for small numbers of sampled genomes: as an example, we show how the generating function can be derived explicitly for three genes under the two-deme IM model. This derivation is done automatically, using Mathematica. Given data from a large number of unlinked and nonrecombining blocks of sequence, these results can be used to find maximum-likelihood estimates of model parameters by tabulating the probabilities of all relevant mutational configurations and then multiplying across loci. The feasibility of the method is demonstrated by applying it to simulated data and to a data set previously analyzed by Wang and Hey (2010) consisting of 26,141 loci sampled from Drosophila simulans and D. melanogaster. Our results suggest that such likelihood calculations are scalable to genomic data as long as the numbers of sampled individuals and mutations per sequence block are small.
Publishing Year
Date Published
2011-11-01
Journal Title
Genetics
Volume
189
Issue
3
Page
977 - 987
IST-REx-ID

Cite this

Lohse K, Harrison R, Barton NH. A general method for calculating likelihoods under the coalescent process. Genetics. 2011;189(3):977-987. doi:10.1534/genetics.111.129569
Lohse, K., Harrison, R., & Barton, N. H. (2011). A general method for calculating likelihoods under the coalescent process. Genetics, 189(3), 977–987. https://doi.org/10.1534/genetics.111.129569
Lohse, Konrad, Richard Harrison, and Nicholas H Barton. “A General Method for Calculating Likelihoods under the Coalescent Process.” Genetics 189, no. 3 (2011): 977–87. https://doi.org/10.1534/genetics.111.129569.
K. Lohse, R. Harrison, and N. H. Barton, “A general method for calculating likelihoods under the coalescent process,” Genetics, vol. 189, no. 3, pp. 977–987, 2011.
Lohse K, Harrison R, Barton NH. 2011. A general method for calculating likelihoods under the coalescent process. Genetics. 189(3), 977–987.
Lohse, Konrad, et al. “A General Method for Calculating Likelihoods under the Coalescent Process.” Genetics, vol. 189, no. 3, Genetics Society of America, 2011, pp. 977–87, doi:10.1534/genetics.111.129569.

Link(s) to Main File(s)
Access Level
OA Open Access

Export

Marked Publications

Open Data IST Research Explorer

Search this title in

Google Scholar