距离法构建进化树讲课文档.pptVIP

距离法构建进化树第1页，共31页。

变形距离法如果假定速率恒定在各谱系间不成立则UPGMA法可能会给出错误的拓扑图第2页，共31页。

第3页，共31页。

不过该拓扑图错误也许能用被称为变形距离法(transformeddistancemethod)的校正方法来订正。简单地说此法是用一个组外单位outgroup作参照物来对被研究谱系间进化速率不等的情况作一些校正，然后对新得到的距离矩阵应用UPGMA法从而推测出该树的拓扑图。第4页，共31页。

组外单位是一个我们对其有外在认识的(OTU)。这些外在认识，如分类学的或古生物学的知识能清楚地表明该(OTU)已先于所有其他被研究的(OTU)而从其共同祖先分岐出来。在现在的例子中，我们假定分类单位对所有别的分类单位来说是一个组外单位，那么即可当作参照物，并用下式来变换距离。第5页，共31页。

这里dij是变形距离i=A,B或C，d=(dAD+dBD+dCD)/3d项的引入是为了保证所有dij的值都是正的在我们的例子里d=37/3而新距离矩阵中关于分类单位A、B和C的值为所以A和B最先被聚类在一起,然后,C再加入到该树中。第6页，共31页。

Least-squaresmethod最小二乘法（LS）

包括一般LS法和加权LS法容许各个分支核苷酸替代速率有所不同算法：实践中，可能的拓扑结构非常多，逐个计算不可能。一般采用计算一个拓扑结构，然后利用分支交换法搜索。第7页，共31页。

在标准的LS法中，所有可能的拓朴结构的相比后，最小的拓朴结构即为最终的树第8页，共31页。

第9页，共31页。

分支长度的估计第10页，共31页。

第11页，共31页。

有时某一被估出的枝长可能会是负的由于实际枝长绝不可能是负的所以最好用来代替这样的估值第12页，共31页。

第13页，共31页。

第14页，共31页。

第15页，共31页。

第16页，共31页。

最小进化（EM）法原理在本方法中，所有分支长度估计的和为：所有可能的拓扑结构都要计算S值，具有最小S值的拓扑结构被挑选作为最优树我们已提及分支长度的LS估计是由距离估计dij的函数给出的，即第17页，共31页。

第18页，共31页。

第19页，共31页。

Neighborjoiningmethod在一个无根两分叉树中,如果两个OTU通过一个内部节点联结,则它们就被说成是近邻(neighbors)比如4个OTU，i、j、m、n。、假定第一个和的值最小，那么我们把i和j对以及m和n对都记1分,其余对记0分。在对所有可能的四单位组都评过分以后,得总分最高的对子即被选为第一个近邻对第20页，共31页。

算法用NJ方法构建一棵星状树，这就是假设不存在聚合群如果估计一棵星状树的分支长度并计算出所有分支的总和（S0），那么这个数目应该比真实树或最终的树的Sf要大第21页，共31页。

第22页，共31页。

第23页，共31页。

第24页，共31页。

第25页，共31页。

第26页，共31页。

第27页，共31页。

Neighborjoiningmethod邻接法的一般步骤如下1．采用了Jukes-Cautor单参数模型，计算出第i条和第j条序列距离dij计算公式：其中q为两个序列中相应位置上相同碱基的概率第28页，共31页。

Neighborjoiningmethod2．计算第i个叶子结点（即第i个序列）的净分歧度ri其中N是叶子结点的个数，dik为叶子结点i和叶子结点k之间的距离第29页，共31页。

Neighborjoiningmethod3．计算任意两两结点i和j之间的速率校正距离（rate-correcteddistance）Mij：4．挑选出最小的速率校正距离Mij5．定义一个新的结点t，t的左、右孩子分别是第i个和第j个结点。结点t到i，j的距离为：第30页，共31页。

谢谢大家！ThankYou！