距离法构建进化树.ppt

《距离法构建进化树.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、OUT:operationaltaxonomicunits分类学单位变形距离法如果假定速率恒定在各谱系间不成立则UPGMA法可能会给出错误的拓扑图不过该拓扑图错误也许能用被称为变形距离法(transformeddistancemethod)的校正方法来订正。简单地说此法是用一个组外单位outgroup作参照物来对被研究谱系间进化速率不等的情况作一些校正，然后对新得到的距离矩阵应用UPGMA法从而推测出该树的拓扑图。组外单位是一个我们对其有外在认识的(OTU)。这些外在认识，如分类学的或古生物学的知识能清楚地表明该(OTU)已先于

2、所有其他被研究的(OTU)而从其共同祖先分岐出来。在现在的例子中，我们假定分类单位对所有别的分类单位来说是一个组外单位，那么即可当作参照物，并用下式来变换距离。这里dij是变形距离i=A,B或C，d=(dAD+dBD+dCD)/3d项的引入是为了保证所有dij的值都是正的在我们的例子里d=37/3而新距离矩阵中关于分类单位A、B和C的值为所以A和B最先被聚类在一起,然后,C再加入到该树中。Least-squaresmethod最小二乘法（LS）包括一般LS法和加权LS法容许各个分支核苷酸替代速率有所不同算法：实践中，可能的拓扑结

3、构非常多，逐个计算不可能。一般采用计算一个拓扑结构，然后利用分支交换法搜索。在标准的LS法中，所有可能的拓朴结构的相比后，　　　最小的拓朴结构即为最终的树分支长度的估计有时某一被估出的枝长可能会是负的由于实际枝长绝不可能是负的所以最好用来代替这样的估值最小进化（EM）法原理在本方法中，所有分支长度估计的和为：所有可能的拓扑结构都要计算S值，具有最小S值的拓扑结构被挑选作为最优树我们已提及分支长度的LS估计是由距离估计dij的函数给出的，即Neighborjoiningmethod在一个无根两分叉树中,如果两个OTU通过一个内部节

4、点联结,则它们就被说成是近邻(neighbors)比如4个OTU，i、j、m、n。、假定第一个和的值最小，那么我们把i和j对以及m和n对都记1分,其余对记0分。在对所有可能的四单位组都评过分以后,得总分最高的对子即被选为第一个近邻对算法用NJ方法构建一棵星状树，这就是假设不存在聚合群如果估计一棵星状树的分支长度并计算出所有分支的总和（S0），那么这个数目应该比真实树或最终的树的Sf要大Neighborjoiningmethod邻接法的一般步骤如下1．采用了Jukes-Cautor单参数模型，计算出第i条和第j条序列距离dij计算

5、公式：其中q为两个序列中相应位置上相同碱基的概率Neighborjoiningmethod2．计算第i个叶子结点（即第i个序列）的净分歧度ri其中N是叶子结点的个数，dik为叶子结点i和叶子结点k之间的距离Neighborjoiningmethod3．计算任意两两结点i和j之间的速率校正距离（rate-correcteddistance）Mij：4．挑选出最小的速率校正距离Mij5．定义一个新的结点t，t的左、右孩子分别是第i个和第j个结点。结点t到i，j的距离为：Neighborjoiningmethod结点t与进化树其他结点

6、k的距离为dtk6．从叶子结点集合L中删除结点i和j，并在集合中添加t结点，N的值减去1。7．如果叶子结点集合L中剩余的叶子结点个数大于2，则重复步骤2继续计算直到叶子结点个数为2，即进化树完全建成