INRAE - BioSP, Centre INRAE PACA, 228 route de l'aérodrome, Domaine St-Paul - Site Agroparc, 84914, Avignon Cedex 9, France.
Department of Archaeogenetics, Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, 04103, Leipzig, Germany.
Theor Popul Biol. 2024 Oct;159:108-124. doi: 10.1016/j.tpb.2024.03.002. Epub 2024 Mar 16.
We introduce a modified spatial Λ-Fleming-Viot process to model the ancestry of individuals in a population occupying a continuous spatial habitat divided into two areas by a sharp discontinuity of the dispersal rate and effective population density. We derive an analytical formula for the expected number of shared haplotype segments between two individuals depending on their sampling locations. This formula involves the transition density of a skew diffusion which appears as a scaling limit of the ancestral lineages of individuals in this model. We then show that this formula can be used to infer the dispersal parameters and the effective population density of both regions, using a composite likelihood approach, and we demonstrate the efficiency of this method on a range of simulated data sets.
我们引入了一个修正的空间 Λ-Fleming-Viot 过程来模拟个体的祖先在一个连续的空间栖息地中的分布情况，该栖息地被扩散率和有效种群密度的急剧不连续分为两个区域。我们推导出了一个关于两个个体之间共享单倍型片段数量的预期的解析公式，这取决于它们的采样位置。这个公式涉及到一个斜扩散的转移密度，它出现在这个模型中个体的祖先谱系的一个标度极限中。然后，我们通过复合似然方法表明，这个公式可以用于推断两个区域的扩散参数和有效种群密度，并在一系列模拟数据集上展示了这种方法的效率。
