期刊
GENETICS
卷 189, 期 3, 页码 977-U398出版社
GENETICS SOCIETY AMERICA
DOI: 10.1534/genetics.111.129569
关键词
-
资金
- European Research Council [250152]
- United Kingdom Natural Environment Research Council [NE/I020288/1]
- MRC [G0701786] Funding Source: UKRI
- NERC [NE/I020288/1] Funding Source: UKRI
- Medical Research Council [G0701786] Funding Source: researchfish
- Natural Environment Research Council [NE/I020288/1] Funding Source: researchfish
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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据