系统发育分析方法

系统发育分析方法第1页/共57页方法基本特征适用范围优点缺点NJ不需要分子钟假设，是基于最小进化原理，进行类的合并时，不仅要求待合并的类是相近的，而且要求待合并的类远离其他的类。远缘序列，进化距离不大，信息位点少的短序列假设少，树的构建相对准确，计算速度快，只得一颗树，可以分析较多的序列，运行速度优于最大简约法序列上的所有位点等同对待，且所分析的序列的进化距离不能太大MP基于进化过程中碱基替代数目最少这一假说，不需要替代模型，对所有可能的拓扑结构进行计算，并计算出所需替代数最小的那个拓扑结构，作为最优树近缘序列物种序列的数目≤12.善于分析某些特殊的分子数据如插入、缺失等序列有用。只适于序列数目N≤12。存在较多回复突变或平行突变时，结果较差。变异大的序列会出现长枝吸引而导致建树错误。ML依赖于某一个特定的替代模型来分析给定的一组序列数据，使得获得的每一个拓扑结构的似然率都为最大值，然后再挑出其中似然率最大的拓扑结构作为最优树。特定的替代的模，远缘序列很好的统计学基础，大样本时似然法可以获得参数统计的最小方差，在进化模型确定的情况下，ML法是与进化事实吻合最好的建树算法.所有可能的系统发育树都计算似然函数，计算量大，耗时时间长。依赖于合适的替代模型，BI基因进化模型的统计推论法，通过后验概率直观反映出各分支的可靠性而不需要自检法检验大而复杂的数据集具有坚实的数学和统计学基础，可以处理复杂和接近实际情况的进化模型对进化模型比较敏感，后验概率是建立在许多假说上，在现实中可能不成立第2页/共57页/phylip/software.html系统发育树构建的软件第3页/共57页常见软件软件名称用途DNAMAN序列分析的综合工具BioEdit

序列分析的综合工具DNASTAR序列分析的综合工具MAFFT多重序列比对工具Muscle多重序列比对工具ClustalX图形化的多序列比对工具;构建N-J系统树Gblocks冗余序列处理工具jModelTest,ModelTest,ModelGenerator进化模型选择工具PHYLIP集成的进化分析工具MEGA图形化、集成的进化分析工具PAUP集成的进化分析工具PHYML,PAML,RAxMLML建树工具MrBayes基于贝叶斯方法的建树工具TreeView进化树显示工具FigTree,AdobeIllustrator进化树显示和编辑工具第4页/共57页系统发育树构建的过程多序列比对(MAFFT)进化模型的选择(ModelTest)系统发育树的构建(RAxML,MrBayes,PAUP)系统发育树显示和编辑(FigTree,AdobeIllustrator）序列拼接(Mega)第5页/共57页序列拼接BioEditMegaSeqmanContigSequencer第6页/共57页多序列比对http://mafft.cbrc.jp/alignment/server/第7页/共57页速度：Muscle>MAFFT>Clustal比对准确性：MAFFT>Muscle>>Clustal第8页/共57页比对前第9页/共57页MAFFT7.0onlinealignment

http://mafft.cbrc.jp/alignment/server/index.html第10页/共57页第11页/共57页第12页/共57页第13页/共57页第14页/共57页PAUP软件使用流程（系统树构建）第15页/共57页1.将比对后的fasta格式文件转换成Nexus格式2.将paup命令粘贴到Nexus文件下方，在命令程序中指定外群，保存。beginpaup;

logfile=p_buffer.txt;

psetcollapse=minbrlen;[ctype1.5_1:all;]setmaxtrees=5000increase=no;

outgroup

****;setcriterion=parsimony;

hsearch

addseq=randomnreps=1000;roottrees

outroot=monophyl;

savetrees

brlens=yesfile=MP.tre;

pscoresALL/ci=yestl=yeshi=yesrc=yesri=yeskhtest=yes;

bootstrapnreps=1000Keepall=yes/AddSeq=randomnreps=10;

roottrees

outroot=monophyl;

savetreesfile=BT.trefrom=1to=1savebootp=bothmaxdec=0;

end;转换文件格式第16页/共57页.....第17页/共57页3.打开paup软件，打开Nexus文件然后运行即可。

第18页/共57页4.运行界面。MP树运行完后，点击“Stop”，继续运行BT树。第19页/共57页运行结果文件：第20页/共57页MP树第21页/共57页BT树第22页/共57页P-buffer文件第23页/共57页RAxML建树第24页/共57页Page

25程序自带的文件：raxmlHPC、raxmlHPC-PTHREADS、run三个

准备文件两个：phy格式的比对好的序列，txt格式的partition文件Run文件基因的名字和序列所在位置跑多少次比对好序列的文件名字输出结果的的名字Partition文件的名字Partition文件第25页/共57页MrBayes建树(Mrmodeltest2.3和Mrbayes3.2.2)第26页/共57页/projects/mrbayes/files/Version3.2.2fixesanumberofbugsinpreviousreleasesofversion3.2.

mrbayes_x86.exefor32bitsystem,mrbayes_x64.exefor64bitsystem第27页/共57页操作步骤1.Fasta文件转换成Nexus格式的文件2.把Mrmodelblock文件夹中对应的MrModelblock*loci文件中的内容粘贴到Nexus文件最下方（几个基因即对应几个loci的文件，比如4个基因，则需要复制MrModelblock4loci中的内容）

第28页/共57页3.修改刚刚粘贴的命令

参看文件红色框中是要修改的地方第29页/共57页4.Paup运行刚刚修改好的Nexus文件，得到SCORES文件5.将SCORES文件复制到MrModeltest2.3文件夹中，此时该文件夹包含以下文件：第30页/共57页6.运行cmd第31页/共57页得到txt文件第32页/共57页7.打开刚刚得到的txt文件，从每一个txt文件中找到如下Bayes的模型：第33页/共57页8.把得到的models复制到一个text文本中9.打开之前的NEXUS文件，删掉mrmodeltest命令,粘贴bayes命令:第34页/共57页beginmrbayes;[Thisblocksetsupseveraldifferentpartitionsthatcouldbeusedintheanalysisofthisdataset]

outgroupM_infuscans_CBS_869_96;[replacefungusXwithyouroutgrouptaxon][Whendefiningyourcharsetsbelow,thecharactersmustfolloweachotherdirectly,e.g.1-300,301-500,501-600andnot5-300,325-500,530-600.Youwillexcludedeverythingyoudonotwanttoincludeintheanalysis(e.g.1-4,301-324and501-529inthecharsetexcludedcharactersline.] charsetlocus1=1-286;[replacethexx'swithnumbersreflectingthecharacterspanningofyourgene1] charsetlocus2=287-605;[replacethexx'swithnumbersreflectingthecharacterspanningofyourgene2] charsetlocus3=606-1159;[replacethexx'swithnumbersreflectingthecharacterspanningofyourgene3]charsetlocus4=1160-1678;[replacethexx'swithnumbersreflectingthecharacterspanningofyourgene4] charsetexcludedchars=282-286601-6051155-11591674-1678;[listhereallofthecharactersthatyoudonotwanttoinclude.e.g.thebitsbetweentheloci] excludeexcludedchars; partitionAllLoci=4:locus1,locus2,locus3,locus4; logstartfilename=mydata.log;end;beginmrbayes;

outgroupM_infuscans_CBS_869_96; setpartition=AllLoci;

prsetapplyto=(1,2,4)statefreqpr=dirichlet(1,1,1,1);[Thismeanslocus1andlocus3aresame] prsetapplyto=(3)statefreqpr=fixed(equal);[Thisisthemodeloflocus2] lsetapplyto=(1,2)nst=2rates=gamma; lsetapplyto=(3,4)nst=2rates=propinv; unlinkshape=(all)pinvar=(all)statefreq=(all)revmat=(all); mcmcpngen=10000000relburnin=yesburninfrac=0.25printfreq=1000samplefreq=1000nchains=4savebrlens=yesstoprule=yesstopval=0.01; mcmc; sumt[conformat=simple];[usingMrbayes3.2.1shouldadd"conformat=simple"]end;红色字体文字是需要修改的地方；蓝色字体视情况修改，也可不修改。第35页/共57页10.把修改好的NEXUS文件放在Mrbayes文件夹中，打开Mrbayes，输入命定:exeNEXUS的文件名.nex11.Over~第36页/共57页CIPRES在线建树第37页/共57页CIPRES在线建ML树（RAxML）

/portal2/login!input.action第38页/共57页Createnewfolder第39页/共57页Uploaddata第40页/共57页第41页/共57页Createnewtask第42页/共57页第43页/共57页Selecttools第44页/共57页Setparam