Transcript 系统进化树的构建方法与软件应用

系统进化树的构建方法与软件应用
姓名：张镜悬
E-mail：[email protected]
•什么是系统进化树
系统进化树又称为演化树，是表明被认为具有共同祖先的各物
种间演化关系的树.在树中每个节点代表其各个分支的最近共
同祖先，而节点的线段长度对应了其演化的距离。
(http://en.wikipedia.org/wiki/Phylogenetic_tree)
•直系同源和旁系同源
•直系同源：同源的基因是由共同的祖先基因进化而产生的。
•旁系同源：同源的基因是由于基因复制产生的。
•这也就告诉我们用于分子进化分析中的序列必须是直系同
源的才可以真实的反映其进化的过程。
•系统进化树的分类
根据树是否有根，进化树可以分为有根树和无根树两类。
•有根树和无根树的进化层面上的意义
有根树反应了树上物种或者基因进化的时间顺序,通过分析有
根树的长度，可以了解不同的物种或者基因以什么方式和速率
进化。
无根树只反映分类单元之间的距离,而不涉及谁是谁的祖先问
题
做有根树需要指定outgroup。所谓out group ， 就是你所分
析的东西之外的一个group。比如你分析人类的不同人种，就
选个chimpanzee，你要分析哺乳动物，就选个鳄鱼乌龟之类，
总之保证它在 你要分析的group之外，但又不太远就行了。将
你选定的东西指定为outgroup，做出来的树就是有根树。out
group可以不只一个，它是一个group。
•系统进化树的结构
branch
node
进化树的结构主要分为三部分：
• 树叶
• 树枝
• 节点
其中我们把从同一个节点上分出的两个分支叫做sister
group.
Sister group 从结构上可以理解为从进化史上看两者非常
接近，其次两者拥有唯一的共同的祖先。
•系统进化树的结构
a
b
c
d
a
b
d
•从结构上看，我们认为这三个树是等价的
c
a
c
d
b
•构建系统进化树的理论方法
•最大简约法(maximum parsimony，MP)最早源于形态性状研究，现在已经推
广到分子序列的进化分析中。最大简约法的理论基础是奥卡姆哲学原则，这
个原则认为：解释一个过程的最好理论是所需假设数目最少的那一个。对所
有可能的拓扑结构进行计算，并计算出所需替代数最小的那个拓扑结构，作
为最优树。
•优点：最大简约法对于分析某些特殊的分子数据如插入、缺失等序列有用。
在分析的序列位点上没有回复突变或平行突变，且被检验的序列位点数很大
的时候，最大简约法能够推导获得一个很好的进化树。
•缺点：在分析序列上存在较多 的回复突变或平行突变，而被检验的序列位
点数又比较少的时候，最大简约法可能会给出一个不合理的或者错误的进化
树推导结果。
•构建系统进化树的理论方法
最大似然法(maximum likelihood,ML)
•最早应用于系统发育分析是在对基因频率数据的分析上，后来基于分子序列
的分析中也已经引入了最大似然法的分析方法。当样本量很大的时候，似然
法可以获得参数统计的最小方差。
•最大似然法分析中，选取一个特定的替代模型来分析给定的一组序列数据，
使得获得的每一个拓扑结构的似然率都为最大值，然后再挑出其中似然率最
大的拓扑结构 作为最优树。
•最大似然进化模型
简单假设所有核苷酸（或者氨基酸）之间相互转变的概率是一样的
程序会把所有可能的核苷酸轮流置于进化树的内部节点上，并且计算每个这
样的序列产生实际数据的可能性。
所有可能再现的几率被加总，产生一个特定点的似然值，然后这个数据集的
所有比对位点的似然值的加和就是整个进化树的似然值。
•构建系统进化树的理论方法
• 邻近法（Neighbor-Joining Method，NJ ）
该方法通过确定距离最近(或相邻)的成对分类单位来使系统树的总距离达到
最小。相邻是指两个分类单位在某一无根分叉树中仅通过一个节点(node)相
连。通过循序地将相邻点合并成新的点，就可以建立一个相应的拓扑树。
•构建系统进化树的所涉及的工具
•PHYLIP
•MEGA
•R
•Matlab
•BioEdit
•TreeView
•PHYML
•ClustalX
•构建系统进化树的所涉及的工具
PHYLIP
由美国华盛顿大学Felsenstein开发，可
以免费下载，适用于绝大多数操作系统
PAUP
由美国simthsonion institute开发，仅
适用于Apple-Macintosh和UNIX操作系统
MEGA
美国宾夕法尼亚州立大学MasatoshiNei
开发的分子进化遗传学软件，图形化，
集成的进行分析工具，不包括ML
MOLPHY
日本国立统计数理研究所开发，最大似
然法构树
PAML
英国University college London开发，
最大似然法构树和分子进化模型
•构建系统进化树的所涉及的工具
PUZZLE
应用quarter puzzling方法（一种最大简
约法）构建系统树
TreeView
英国University of Glasgow开发，进化
树显示工具
Phylogeny
欧洲生物信息研究所（EBI）的系统发育
分析软件
PHYML
快速的ＭＬ建树工具
MrBayes
基于贝叶斯方法的建树工具
MAC5
基于贝叶斯方法的建树工具

构建ＮＪ树，可以用PHYLIP或者MEGA

构建MP树，可以使用PHYLIP或者MEGA

构建ML树可以使用PHYML，速度快，同时构建ML树还可以用PHYLIP，或者可
以使用BioEdit

贝叶斯的算法以MrBayes为代表，不过速度比较慢
关于系统发育分析的更多知识请参阅：
http://www.bioon.com/biology/bioinfo2/78842.shtml
•构建系统进化树的主要步骤
•大体来说构建系统进化树的步骤有三步：
1.序列比对 (ClustalX2)
•构建系统进化树的主要步骤
2. 掐头去尾 选取所需序列 转换格式（BioEdit or
ClustalX2）
Example：
3. 利用相关软件绘制系统进化树（BioEdit，MEGA）
•实例讲解
下面的内容将教大家如何来构建自己的系统进化树。
首先我们需要弄清楚一个很重要的问题，什么是Fasta 格式?
在生物信息学中，FASTA格式（又称为Pearson格式），是一种基于
文本用于表示核苷酸序列或氨基酸序列的格式。在这种格式中碱基对或
氨基酸用单个字母来编码，且允许在序列前添加序列名及注释。序列文
件的第一行是由大于号“>”或分号“;”打头的任意文字说明（习惯常
用“>”作为起始），用于序列标记。从第二行开始为序列本身，只允
许使用既定的核苷酸或氨基酸编码符号。
• 构建我们自己的Fasta 文件
很多情况下，Fasta文件
是直接可以从数据库中下
载得到的，但是根据实际
要求的不同，有时候我们
需要自己构建Fasta文件，
如果您已近有了想用来构
建进化树的序列，您可以
如右图所示构建自己的文
件，文件的保存格式是：
文件名.txt
•实例讲解
下面我们以禽流感病毒为例，构建系统进化树。
首先我们要下载我们所需的序列。
http://www.ncbi.nlm.nih.gov/genomes/FLU/Database/n
ph-select.cgi?go=genomeset
•实例讲解
请在Define search set: 中选择我们想要的禽流感病毒
的Type, Host, Country/Region, Subtype.
这里我们选在了A型
禽流感病毒，当然在
这次练习中您喜欢的
任意类型。
•实例讲解
请在Define search set: 中选择我们想要的禽流感病毒
的Type, Host, Country/Region, Subtype.
这里我们选在了A型禽
流感病毒，当然在这次
练习中您喜欢的任意类
型。当您确定之后请点
击 Show results
•实例讲解
当您点击完 Show results 之后你要做的就是选在我
们所需的序列了
•实例讲解
因为禽流
感病毒不像别
的很多别的病
毒只有核苷酸
序列，它拥有
八个或者七个
Negative sense RNA。
•实例讲解
这里我们只要选中其
中一种就可以了，比
如说ＨＡ，我们可以
选择２０个样本来构
建系统进化树。样本
选择完之后请点击
Download， 文件类型
选 择 Nucleotide
(Fasta)，并把文件保
存在计算机您熟悉的
地方。（当然根据需
求的不同您也可以选
在蛋白序列）
•实例讲解
文件下载完之后，下载的Fasta 文件直接用 ClustalX
2.0.12打开
•实例讲解
在进行多序列比对之前我们需要对软件进行一些设置
1.选择Alignment标签
2.选择Output format options
请将Clustalw sequences
numbers选项设置为 On
之后点击 Ok ，在返回主界
面之后请点击Alignment 标
签选择 Do Complete
Alignment选项
选择保存路径之后点击ok，剩下的时
间可以去喝点咖啡休息一下。
•实例讲解
点击主界面
中的 File标
签选择 Save
as 选 项 ， 并
按照例子设
置参数
从图中我们可以发现起
始序列最短的是从位置
22开始的，而尾端序列
最短的是在位置1738，
通过设置我们可以保存
这样一批已经经过掐头
去尾后的序列,保存格
式为：文件名.aln。当
然我们也可以直接保存
为Fasta format, 如果
选择前者我们需要用
BioEdit 转 换 格 式 ， 如
果是后者我们可以直接
进入建树阶段。
•实例讲解
经过ClustalX2掐头去尾后的序列可以用BioEdit软件打开，
选择FileSave as保存类型为:文件名.fasta. 当我们查
询结果的时候可以发现这和用ClustalX2保存的fasta文件是
一致的。
•实例讲解
下一步我们将介绍如何用MEGA构建我们的进化树，首先请
大家用MEGA软件将我们之前保留的Fasta文件打开。
•实例讲解
下一步我们将介绍如何用MEGA构建我们的进化树，首先请
大家用MEGA软件将我们之前保留的Fasta文件打开这时候会有
两个窗口，选择File标签-->Convert file format to Mega.
•实例讲解
选择File标签-->Convert file format to Mega.
当给出相应的文件路径之后点击ok ,然后制
定输出文件格式：文件名.meg
•实例讲解
双击刚才保存的meg文件.
选择数据类型，在本次测试中我们
用的是核苷酸序列，对于右边的参
数信息请点击help按钮。
更具实际的情况我
们这里选择No选项
•实例讲解
下一步进入建树的最后阶段
在Plylogeny中选择建树
方法，这里我们选择NJ法。
蛋白质序列一般选择
Poisson Correction（泊
松校正）,对于核苷酸序列
一般采用Kimura-2模型
参数设置好之后点
击compute.
•实例讲解
根据Mega的计算最终我们得到了序列中的进化关系。
Mega软件还可以自动提供一份简
要的分析报告，你只需要点击
Caption按钮报告便可以自动生
成。
如果Bootstrap Value >70我们
认为这个分支是可靠的

进化树评估优化方法简介：常用的两种方法就是Bootstrap和Jackknife。
所谓Bootstraping法 就是从整个序列的碱基（氨基酸）中任意
选取一半，剩下的一半序列随机补齐组成一个新的序列。这样，一个序
列就可以变成了许多序列，一个多序列组也就可以变 成许多个多序列组。
根据某种算法（最大简约性法、最大可能性法、邻位相连法）每个多序
列组都可以生成一个进化树。将生成的许多进化树进行比 较，按照多数
规则（majority-rule）我们就会得到一个最“逼真”的进化树。其数值
反应了该树枝的可信的百分比。
所谓Jackknife则是另外一种随机选取序列的方法。它与
Bootstrap法的区别是不将剩下的一半序列补齐，只生成一个缩短了一半
的新序列。

Double Check
通常情况下当我们用建树的一种方法获得了树图之后，我们建议大家
可以通过另外的方法建树，如果先后的两个树图大体一致，我们认为之
前构建的树是相对可靠的。
Thanks