进化树分析及相关软件使用说明

《进化树分析及相关软件使用说明》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、进化树分析及相关软件使用说明我在此介绍几个进化树分析及其相关软件的使用和应用范围。这几个软件分别是PHYLIP、PUZZLE、PAUP、TREEVIEW、CLUSTALX和PHYLO-WIN（LINUX）。在介绍软件之前，我先简要地叙述一下有关进化树分析的一些方法学问题。进化树也称种系树，英文名叫“Phyligenetictree”。对于一个完整的进化树分析需要以下几个步骤：⑴要对所分析的多序列目标进行排列（Toalignsequences）。做ALIGNMENT的软件很多，最经常使用的有CLUSTALX和CLUSTALW，前者是在WINDOW下的而后者是在DOS下的。⑵要构建一个

2、进化树（Toreconstrutphyligenetictree）。构建进化树的算法主要分为两类：独立元素法（discretecharactermethods）和距离依靠法（distancemethods）。所谓独立元素法是指进化树的拓扑形状是由序列上的每个碱基/氨基酸的状态决定的（例如：一个序列上可能包含很多的酶切位点，而每个酶切位点的存在与否是由几个碱基的状态决定的，也就是说一个序列碱基的状态决定着它的酶切位点状态，当多个序列进行进化树分析时，进化树的拓扑形状也就由这些碱基的状态决定了）。而距离依靠法是指进化树的拓扑形状由两两序列的进化距离决定的。进化树枝条的长度代表着进化距离

3、。独立元素法包括最大简约性法（MaximumParsimonymethods）和最大可能性法（MaximumLikelihoodmethods）；距离依靠法包括除权配对法（UPGMAM）和邻位相连法（Neighbor-joining）。⑶对进化树进行评估。主要采用Bootstraping法。进化树的构建是一个统计学问题。我们所构建出来的进化树只是对真实的进化关系的评估或者模拟。如果我们采用了一个适当的方法，那么所构建的进化树就会接近真实的“进化树”。模拟的进化树需要一种数学方法来对其进行评估。不同的算法有不同的适用目标。一般来说，最大简约性法适用于符合以下条件的多序列：i所要比较的

4、序列的碱基差别小，ii对于序列上的每一个碱基有近似相等的变异率，iii没有过多的颠换/转换的倾向，iv所检验的序列的碱基数目较多（大于几千个碱基）；用最大可能性法分析序列则不需以上的诸多条件，但是此种方法计算极其耗时。如果分析的序列较多，有可能要花上几天的时间才能计算完毕。UPGMAM（Unweightedpairgroupmethodwitharithmeticmean）假设在进化过程中所有核苷酸/氨基酸都有相同的变异率，也就是存在着一个分子钟。这种算法得到的进化树相对来说不是很准确，现在已经很少使用。邻位相连法是一个经常被使用的算法，它构建的进化树相对准确，而且计算快捷。其缺点

5、是序列上的所有位点都被同等对待，而且，所分析的序列的进化距离不能太大。另外，需要特别指出的是对于一些特定多序列对象来说可能没有任何一个现存算法非常适合它。最好是我们来发展一个更好的算法来解决它。但无疑这是非常难的。我想如果有人能建立这样一个算法的话，那他（她）完全可以在Proc.Natl.Acad.Sci.USA.上发一篇高质量的文章。下面介绍几个软件的使用。首先是PHYLIP。其是多个软件的压缩包，下载后双击则自动解压。当你解压后就挥发现PHYLIP的功能极其强大，主要包括五个方面的功能软件：i，DNA和蛋白质序列数据的分析软件。ii，序列数据转变成距离数据后，对距离数据分析的软

6、件。iii，对基因频率和连续的元素分析的软件。iv，把序列的每个碱基/氨基酸独立看待（碱基/氨基酸只有0和1的状态）时，对序列进行分析的软件。v，按照DOLLO简约性算法对序列进行分析的软件。vi，绘制和修改进化树的软件。在此，我主要对前两种功能软件进行说明。我们现在有几个序列如下：Mo3ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGCACGGTACCATMo5ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGTACGGTACCATMo6ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGTACG

7、GTACCATMo7ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGTACAGTACCATMo8ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGTACAGTACCATMo9ATGTATCTCGTACATTACTGCCAGCCACCATGAATATTGTACGGTACCATMo12ATGTATTTCGTACATTACTGCCAGCCACCATGAATATTGTACGGTACCATMo13ATGTATCTCGTA