《进化树选择》由会员分享，可在线阅读，更多相关《进化树选择（2页珍藏版）》请在金锄头文库上搜索。
1、系统进化树构建及数据分析的简介Postedon08六月2009by柳城，阅读1,278简洁版繁體一、引言开始动笔写这篇短文之前，我问自己，为什么要写这样的文章？写这样的文章有实际的意义吗？我希望能够解决什么样的问题？带着这样的疑惑，我随手在丁香园（DXY）上以关键字“进化分析求助”进行了搜索，居然有289篇相关的帖子（2006年9月12日）。而以关键字“进化分析”和“进化”为关键字搜索，分别找到2,733和7,724篇相关的帖子。考虑到有些帖子的内容与分子进化无关，这里我保守的估计，大约有3,0004,000篇帖子的内容，是关于分子进化的。粗略地归纳一下，我大致将提出的问题分为下述的几类：1涉
2、及基本概念例如，“分子进化与生物进化是不是一个概念”，“关于微卫星进化模型有没有什么新的进展”以及“关于Kruglyak的模型有没有改进的出现”，等等。2关于构建进化树的方法的选择例如，“用boostrapNJ得到XX图，请问该怎样理解？能否应用于文章？用boostraptest中的ME法得到的是XXX树，请问与上个树比，哪个更好”，等等。3关于软件的选择例如，“想做一个进化树，不知道什么软件能更好的使用且可以说明问题，并且有没有说明如何做”，“拿到了16srRNA数据，打算做一个系统进化树分析，可是原来没有做过这方面的工作啊，都要什么软件”，“请问各位高手用ClustalX做出来的进化树与p
3、hylip做的有什么区别”，“请问有做过进化树分析的朋友，能不能提供一下，做树的时候参数的设置，以及代表的意思。还有各个分支等数值的意思，说明的问题等”，等等。4蛋白家族的分类问题例如，“搜集所有的关于一个特定domain的序列，共141条，做的进化树不知具体怎么分析”5新基因功能的推断例如，“根据一个新基因A氨基酸序列构建的系统发生树，这个进化树能否说明这个新基因A和B同源，属于同一基因家族”，等等。6计算基因分化的年代例如，“想在基因组水平比较两个或三个比较接近物种之间的进化年代的远近，具体推算出他们之间的分歧时间”，“如何估计病毒进化中变异所需时间”，等等。7进化树的编辑例如生成的进化树
4、图片，如何进行后续的编辑，比如希望在图片上标注某些特定的内容，等由于相关的帖子太多，作者在这里对无法阅读全部的相关内容而致以歉意。同时，作者归纳的这七个问题也并不完全代表所有的提问。对于问题1所涉及到的基本的概念，作者推荐读者可参考由MasatoshiNei与SudhirKumar所撰写的分子进化与系统发育（MolecularEvolutionandPhylogenetics）一书，以及相关的分子进化方面的最新文献。对于问题7，作者之一lylover一般使用Powerpoint进行编辑，而Photoshop、Illustrator及Windows自带的画图工具等都可以使用。这里，作者在这里对问
5、题2-6进行简要地解释和讨论，并希望能够初步地解答初学者的一些疑问。二、方法的选择首先是方法的选择。基于距离的方法有UPGMA、ME（MinimumEvolution，最小进化法）和NJ（Neighbor-Joining，邻接法）等。其他的几种方法包括MP（Maximumparsimony，最大简约法）、ML（Maximumlikelihood，最大似然法）以及贝叶斯（Bayesian）推断等方法。其中UPGMA法已经较少使用。一般来讲，如果模型合适，ML的效果较好。对近缘序列，有人喜欢MP，因为用的假设最少。MP般不用在远缘序列上，这时一般用NJ或ML。对相似度很低的序列，NJ往往出现Lon
6、g-branchattraction（LBA，长枝吸引现象），有时严重干扰进化树的构建。贝叶斯的方法则太慢。对于各种方法构建分子进化树的准确性，一篇综述（HallBG.MolBiolEvol2005,22（3）:792-802）认为贝叶斯的方法最好，其次是ML,然后是MP。其实如果序列的相似性较高，各种方法都会得到不错的结果，模型间的差别也不大。对于NJ和ML，是需要选择模型的。对于各种模型之间的理论上的区别，这里不作深入的探讨，可以参看Nei的书。对于蛋白质序列以及DNA序列，两者模型的选择是不同的。以作者的经验来说，对于蛋白质的序列，一般选择PoissonCorrection（泊松修正）这
7、一模型。而对于核酸序列，一般选择Kimura2-parameter（Kimura-2参数）模型。如果对各种模型的理解并不深入，作者并不推荐初学者使用其他复杂的模型。Bootstrap几乎是一个必须的选项。一般Bootstrap的值70，则认为构建的进化树较为可靠。如果Bootstrap的值太低，则有可能进化树的拓扑结构有错误，进化树是不可靠的。对于进化树的构建，如果对理论的了解并不深入，作者推荐使用缺省的参数。需要选择模型的时候（例如用NJ或者ML建树），对于蛋白序列使用PoissonCorrection模型，对于核酸序列使用Kimura-2参数模型。另外需要做Bootstrap检验，当Bootstrap值过低时，所构建的进化树其拓扑结构可能存在问题。并且，一般推荐用两种不同的方法构建进化树，如果所得到的进化树类似，则结果较为可靠。三、软件的选择表1中列出了一些与构建分子进化树相关的软件。构建NJ树，可以用PHYLIP（写得有点问题，例如比较慢，并且Bootstrap检验不方便）或者MEGAoMEGA是Nei开发的方法并设计的图形化的软件，使用非常方便。作者推荐MEGA软件为初学者的首选。虽然多雪列比对工具ClustalW/X自带了一个NJ的建树程序，但是该程序只有p-distance模型，而且构建的树不够准确，一般不用来构建进化树。
