植物分子进化中系统发生树重建方法

1、植物分子进化中系统发生树重建方法摘要：在现代分子进化研究中，根据现有生物基因或物种多样性重建生物的进化史是一个非常重要的问题，一个可靠的系统发生的推断，将揭示出有关生物进化过程的顺序，有助于我们了解生物进化的历史和进化机制，本文就植物分子进化中系统发生树的重建方法进行了阐述。关键词:植物分子进化系统发生树重建方法运用dna序列测定已经普遍运用于系统与进化植物学的 研究中，从dna水平研究生物多样性与生物进化的分子生物 学技术作为一种研究手段，通过揭示dna分子中核昔酸的变 异研究植物的系统发生、种内分化及遗传多样性等，正在改 变植物系统这门古老而经典的学科，为解决过去依靠植物形 态和稀缺的化石

2、证据在分类学、系统发育、物种形成与进化 等方面的难题提供极为有效的途径。本文就植物分子进化中 系统发生树的常用重建方法做介绍。1.分子系统发生树的重建目前，利用分子生物学数据重建系统发生树的方法很 多。在重建时，对不同类型的数据应采取不同的重建方法。1. 1非加权组平均法(upgma)upgma法在算法上较简单。聚类时，首先将距离最小的 两个otu聚在一起，形成一个新的otu,其分支点位于两个 otu间距离的1/2处，然后计算新的otu与其他otu间的平 均距离，再找出其中最小的两个otu进行聚类。如此反复， 直到所有的otu都聚到一起，最终得到一个完整的系统发生 树。1.2邻结法（nj）与u

3、pgma法相比，nj法在算法上相对较复杂，它跟踪的 是树上的节点而不是otuso在聚类过程中，根据原始距离矩 阵，基于其他所有节点间的平均趋异程度而对每对节点间的 距离做调整，即将每个otu的趋异程度标准化，从而形成一 个新的距离矩阵。重建时将距离最小的两个终节点连接起 来，在树中增加一个共祖节点，同时去除原初的两个终节点 及其分支，即对整个树进行修剪。随后，新增加的共祖节点 被视成终节点，重复上一次循环。在每一次循环过程中，都 有两个终节点被一个新的共祖节点取代。整个循环直到只有 两个终节点时为止。从所得到的系统发生树来看，两个聚在 一起的out其所在的终节点到共祖节点的距离并不一定相 同。

4、1.3最大简约法（mp）mp方法利用的只是对简约分析能提供信息的特征。如在 dna序列数据中，利用的只是存在核昔酸序列差异（至少有 2种不同类型的核昔酸序列）的位点，这些位点称为简约信 息位点。利用mp方法重建系统发生树，实际上是一个对给 定otus,其所有可能的树进行比较的过程。对某一个可能的 树，首先对每个位点祖先序列的核昔酸组成做出推断，然后 统计每个位点用来阐明差异的核昔酸最小替换数目。在整个 树中，所有信息简约位点最小核昔酸替换数的总和称为树的 长度。通过比较所有可能树，选择其中长度最小的树作为最 终的系统发生树，即最大简约树。1.4最大似然法（ml）利用最大似然法推断一组序列的系统

5、发生树，需首先确 定序列进化的模型，目前使用较多的是一些相对较简单的模 型，如jukes cantor模型、kimura二参数模型及一般二参 数模型。这些模型都建立在一定假设基础上，然后基于一定 的模型考虑两个otu序列间的关系，找到支的长度。这个过 程需要寻找在某一进化距离上由第一种序列真正转换成第 二种序列的可能性，并确定在最大可能下的进化距离。接着 将多个otu构成的所有可能树作为最佳树，对重建每个树的 统计量进行似然估计，最后通过对树长度的优化，从而获得 最佳树各参数的最大似然估计。1. 5贝叶斯推断法最常用的是mcmc法，其基本思想是构造出一条马尔柯 夫链，该链的状态空间为统计模型参

6、数和不变后验分布参 数。链的构造由多步完成，每步状态空间的状态都被推荐为 链的下一个连接点。首先在状态空间中随机挑选一个状态作 为链的当前态，随机扰动当前态各参数，从状态空间中推荐 一个新态，计算推荐态的相对后验概率密度，若推荐态的后 验概率密度高于当前态，则链的移动被接受，推荐态则作为 下一循环的当前态。若推荐态的后验概率密度低于当前态， 则计算由metropolis等和hasting提出的推荐态与当前态 后验概率的比率，该值接近1时接受推荐态，接近0时则拒 绝推荐态，此时，当前态作为自身的下一个连接点。对上述 过程重复若干次，最终马尔柯夫链将停留在后验概率高的状 态，某态的后验概率就是马尔柯夫链停留该态的时间分值。1. 6不同构树方法的分析在上述方法中，farris和penny认为距离法得到的结果 要比离散特征法的差，felsenstein和nei则认为farris和 penny的论据是对距离法的错误理解。事实上，在有些情况 下距离法能得到比离散特征法更正确的系统发生树。在距离 法中，upgma比较简单且实用，当使用的距离数据是来源于 对含核昔酸数量较多的多个基因的分析结果时，利用upgma 法能得到可靠的系统发生树。在离散特征法中，在不同世系 间进化速率相差较大，且进化