欢迎光临散文网 会员登陆 & 注册

iqtree基于ML法建立系统发育树的步骤

2023-03-07 17:19 作者:鱼了个渔  | 我要投稿

IQtree是利用最大似然法构建系统发生树的,具有高效建树、模型选择、超快自展等特点,适用于大数据的系统发育分析

1.   首先在NCBI下载用于建树的序列,把这些序列放到同一个fasta文件里,每条序列的标签名要以“>”开始并且序列标签名中不要有空格和标点等,不然后续无法识别. IQtree支持多种序列数据格式,如:Phylip/ fasta/fas/nexus/clustlw等格式,程序也会自动转为Phylip格式,省去格式转换的步骤。如图所示:

2.将上述的fasta格式序列文件进行序列比对,我一般采用如下两种方法:

    采用Geneious进行Multiple Align(多序列比对),比对后保存为Nexus格式;再到Mesquite中进行序列剪切,保存文件格式为Simplified NEXUS;最后在Geneious中将文件转化为fasta格式。

     采用MEGA进行alignment序列处理,然后再保存成fasta格式。

如图所示(采用方法1):


3. 在应用ML法重建系统发育树前,需要先进行核苷酸替换模型的选择,可以直接用IQ-TREE的ModelFinder这个模块,能快速获得最佳进化模型,也可以使用jModelTest、ModelTest-NG等选择最佳模型。

4.iqtree是要使用命令运行的。在建树前,首先需要调用出Window的命令运行界面Win键+R,输入cmd打开命令提示符,点击确定。如图所示:

5.因为我的 iqtree2.exe 安装在了D盘,所以要先切换到D盘, 然后回车。

6.之后输入cd+iqtree2.exe的安装路径,再回车;把iqtree2.exe拖入,回车

看上边9点,需要什么就复制哪一条指令,找最适模型+建树,可选择1,例如:

iqtree2 -s example.phy —— iqtree2 -s +自己序列的名字,然后,回车就可以运算了

 除此之外,也可以直接在命名提示符里写指令,具体操作如下:前5步都一样,从第6步开始稍微有所不同

  • 转盘之后,输入cd+iqtree2.exe 的安装路径,再回车

  • 接着测式最佳模型,我自己在记事本中提前写好了指令,将测式最佳模型的指令复制粘贴

  • 然后将比对好的序列fasta文件拖入,回车

  • 这样就找到了最佳模型,有fasta序列的那个文件夹将会出现几个文件,先不用管,后期建树的时候这些文件都会被覆盖

  • 接着就开始建树了,复制粘贴建树指令,回车

  •  接着继续输入指令,然后开始跑树(序列多的需要的时间很长), 再做一次之前的操作,需要加-redo,防止结果被误覆盖.

 结果分析 

建树完成后,会在序列文件夹中生成多个文件,覆盖之前产生的文件,主要包括程序运行日志、一致树文件(example.fasta.contree)、ML树文件(example.fasta.treefile)即TREEFILE文件(含有Bootstrap/UFBoot/SH-aLRT值的BP评估分支置信度),可以用Figtree进行查看和赋根,使用Adode Photoshop 对进化树进行编辑和美化.

注:此方法的步骤是我看各种学习视频和互联网资料之后自己总结的,若有错误的地方或有更好的学习步骤,欢迎指出和分享。





iqtree基于ML法建立系统发育树的步骤的评论 (共 条)

分享到微博请遵守国家法律