通过使用溯祖方法对n个亚群中的迁移矩阵和有效种群大小进行最大似然估计。
Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach.
作者信息
Beerli P, Felsenstein J
机构信息
Department of Genetics, University of Washington, Seattle, WA 98195-7360, USA.
出版信息
Proc Natl Acad Sci U S A. 2001 Apr 10;98(8):4563-8. doi: 10.1073/pnas.081068098. Epub 2001 Apr 3.
Abstract
A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source-sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.
摘要
开发了一种基于合并的最大似然估计器，用于估计不等迁移率和不同亚种群大小的情况。该方法采用马尔可夫链蒙特卡罗方法来研究具有分支长度和迁移事件的可能系谱。通过使用来自四群体n岛模型和源 - 汇种群模型的模拟数据，展示了新方法的特性。我们在migrate中编码的估计方法与genetree进行了测试比较；两个程序产生的似然表面非常相似。即使马尔可夫链从不利参数开始，该算法也能相当快地收敛到估计值。该方法被用于通过使用来自三个人类群体的线粒体DNA数据来估计尼罗河流域的基因流。
相似文献
1
Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach.通过使用溯祖方法对n个亚群中的迁移矩阵和有效种群大小进行最大似然估计。
Proc Natl Acad Sci U S A. 2001 Apr 10;98(8):4563-8. doi: 10.1073/pnas.081068098. Epub 2001 Apr 3.
2
Exact computation of coalescent likelihood for panmictic and subdivided populations under the infinite sites model.无限位点模型下，均匀群体和分裂群体的合并似然的精确计算。
IEEE/ACM Trans Comput Biol Bioinform. 2010 Oct-Dec;7(4):611-8. doi: 10.1109/TCBB.2010.2.
3
Effect of unsampled populations on the estimation of population sizes and migration rates between sampled populations.未抽样群体对抽样群体之间种群大小估计和迁移率的影响。
Mol Ecol. 2004 Apr;13(4):827-36. doi: 10.1111/j.1365-294x.2004.02101.x.
4
Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Coalescent-based estimation of population parameters when the number of demes changes over time.2005年SMBE三国青年研究者研讨会会议记录。当种群数量随时间变化时基于溯祖理论的种群参数估计。
Mol Biol Evol. 2006 May;23(5):988-96. doi: 10.1093/molbev/msj111. Epub 2006 Feb 22.
5
Comparison of Bayesian and maximum-likelihood inference of population genetic parameters.群体遗传参数的贝叶斯推断与最大似然推断比较
Bioinformatics. 2006 Feb 1;22(3):341-5. doi: 10.1093/bioinformatics/bti803. Epub 2005 Nov 29.
6
Estimating effective population size and mutation rate from sequence data using Metropolis-Hastings sampling.使用Metropolis-Hastings抽样法从序列数据估计有效种群大小和突变率。
Genetics. 1995 Aug;140(4):1421-30. doi: 10.1093/genetics/140.4.1421.
7
An efficient Monte Carlo method for estimating Ne from temporally spaced samples using a coalescent-based likelihood.一种基于溯祖似然性，利用时间间隔样本估计有效种群大小（Ne）的高效蒙特卡罗方法。
Genetics. 2005 Jun;170(2):955-67. doi: 10.1534/genetics.104.038349. Epub 2005 Apr 16.
8
A continuous method for gene flow.连续基因流方法。
Genetics. 2013 Jul;194(3):687-96. doi: 10.1534/genetics.113.150904. Epub 2013 May 11.
9
Markov chain Monte Carlo without likelihoods.无似然马尔可夫链蒙特卡罗方法。
Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15324-8. doi: 10.1073/pnas.0306899100. Epub 2003 Dec 8.
10
Maximum-likelihood and markov chain monte carlo approaches to estimate inbreeding and effective size from allele frequency changes.基于等位基因频率变化估计近亲繁殖和有效种群大小的最大似然法和马尔可夫链蒙特卡罗方法
Genetics. 2003 Jul;164(3):1189-204. doi: 10.1093/genetics/164.3.1189.
引用本文的文献
1
Climate change-driven northward expansion of the mediterranean orchid Ophrys apifera from genetic and ecological perspectives.从遗传和生态角度看气候变化驱动的地中海兰花蜂兰向北扩张
Sci Rep. 2025 Aug 13;15(1):29634. doi: 10.1038/s41598-025-15167-6.
2
Robust phylodynamic inference and model specification for HIV transmission dynamics.用于HIV传播动力学的稳健系统发育动力学推断和模型规范。
Epidemics. 2025 Jul 16;52:100846. doi: 10.1016/j.epidem.2025.100846.
3
Genetic diversity and structure of Oncomelania hupensis snails in an area where Schistosoma japonicum transmission has been interrupted for nearly 30 years.日本血吸虫传播已中断近30年地区的湖北钉螺遗传多样性与结构
Parasite. 2025;32:38. doi: 10.1051/parasite/2025031. Epub 2025 Jun 24.
4
Bayesian Phylodynamic Inference of Multitype Population Trajectories Using Genomic Data.使用基因组数据对多类型群体轨迹进行贝叶斯系统发育动力学推断
Mol Biol Evol. 2025 Jun 4;42(6). doi: 10.1093/molbev/msaf130.
5
Genetic insights for enhancing conservation strategies in captive and wild Asian elephants through improved non-invasive DNA-based individual identification.通过改进基于非侵入性DNA的个体识别来增强圈养和野生亚洲象保护策略的遗传学见解。
PLoS One. 2025 May 12;20(5):e0320480. doi: 10.1371/journal.pone.0320480. eCollection 2025.
6
Effects of Social Structure on Effective Population Size Change Estimates.社会结构对有效种群大小变化估计的影响。
Evol Appl. 2025 Jan 14;18(1):e70063. doi: 10.1111/eva.70063. eCollection 2025 Jan.
7
Anthropogenic activities along the Lake Nyasa catchments alter the habitat and genetic diversity of a Lake Salmon, .尼亚萨湖流域的人为活动改变了湖鲑的栖息地和遗传多样性。
Heliyon. 2024 Oct 9;10(20):e39051. doi: 10.1016/j.heliyon.2024.e39051. eCollection 2024 Oct 30.
8
Population connectivity and size reductions in the Anthropocene: the consequence of landscapes and historical bottlenecks in white forsythia fragmented habitats.人类世的种群连通性和规模缩小：景观和历史瓶颈对白荚蒾破碎栖息地的影响。
BMC Ecol Evol. 2024 Oct 10;24(1):123. doi: 10.1186/s12862-024-02308-0.
9
Inference and applications of ancestral recombination graphs.祖先重组图的推断与应用
Nat Rev Genet. 2025 Jan;26(1):47-58. doi: 10.1038/s41576-024-00772-4. Epub 2024 Sep 30.
10
Anatomy of a pest control failure: introgression of cytochrome P450 337B3 alleles from invasive old-world bollworm into native corn earworm (Lepidoptera: Noctuidae).害虫防治失败的剖析：来自入侵的旧世界棉铃虫的细胞色素 P450 337B3 等位基因向本地玉米穗虫（鳞翅目：夜蛾科）的渗入。
J Insect Sci. 2024 Jul 1;24(4). doi: 10.1093/jisesa/ieae094.
本文引用的文献
1
The Theoretical Variance within and among Subdivisions of a Population That Is in a Steady State.处于稳态的种群各亚群内部及之间的理论方差。
Genetics. 1952 May;37(3):312-21. doi: 10.1093/genetics/37.3.312.
2
Evolution in Mendelian Populations.孟德尔群体中的进化。
Genetics. 1931 Mar;16(2):97-159. doi: 10.1093/genetics/16.2.97.
3
Maximum likelihood estimation of recombination rates from population data.基于群体数据的重组率最大似然估计。
Genetics. 2000 Nov;156(3):1393-401. doi: 10.1093/genetics/156.3.1393.
4
Sampling among haplotype resolutions in a coalescent-based genealogy sampler.基于溯祖的谱系采样器中单体型分辨率下的采样
Genet Epidemiol. 2000;19 Suppl 1:S15-21. doi: 10.1002/1098-2272(2000)19:1+<::AID-GEPI3>3.0.CO;2-V.
5
Usefulness of single nucleotide polymorphism data for estimating population parameters.单核苷酸多态性数据在估计群体参数方面的实用性。
Genetics. 2000 Sep;156(1):439-47. doi: 10.1093/genetics/156.1.439.
6
Likelihood analysis of ongoing gene flow and historical association.正在进行的基因流动和历史关联的似然分析。
Evolution. 2000 Feb;54(1):44-50. doi: 10.1111/j.0014-3820.2000.tb00006.x.
7
Inference from gene trees in a subdivided population.在一个细分群体中从基因树进行推断。
Theor Popul Biol. 2000 Mar;57(2):79-95. doi: 10.1006/tpbi.1999.1447.
8
A simple method of removing the effect of a bottleneck and unequal population sizes on pairwise genetic distances.一种消除瓶颈效应和种群大小不均对成对遗传距离影响的简单方法。
Proc Biol Sci. 2000 Jan 7;267(1438):81-7. doi: 10.1098/rspb.2000.0970.
9
Substitution rate variation among sites in mitochondrial hypervariable region I of humans and chimpanzees.人类和黑猩猩线粒体高变区I中各位点间的替换率差异。
Mol Biol Evol. 1999 Oct;16(10):1357-68. doi: 10.1093/oxfordjournals.molbev.a026046.
10
Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach.使用溯祖方法对两个种群的迁移率和有效种群数量进行最大似然估计。
Genetics. 1999 Jun;152(2):763-73. doi: 10.1093/genetics/152.2.763.
Zotero 插件