IQ-tree2构建系统发育树使用小结

IQ-TREE是一款用于系统发育推断的高效软件(Efficient software for phylogenomic inference),其主要具有以下特点:
1.高效的树重建
一种通过最大似然推断系统发育树的快速有效的随机算法。IQ-TREE在计算时间相似的可能性方面优于RAxMLPhyMLNguyen等人,2015)。

2.ModelFinder - 快速准确的模型选择
ModelFinderKalyaanamoorthy等人,2017)支持自由速率变化模型,并且比jModelTestProtTest快10到100倍。它还可以找到最合适的分区方案,如PartitionFinder

3.UFBoot - 超快引导近似
UFBoot提供近似无偏的分支支持值,运行速度比非参数引导程序快100倍,比RAxML快速引导程序快10到40倍(Minh等人,2013)。

4.支持大数据分析
通过检查点、安全数值和低内存模式支持具有数千个序列或数百万个比对位点的庞大数据集。多核CPU并行MPI系统用于加速分析。

准备工作

wget -c https://github.com/iqtree/iqtree2/releases/download/v2.1.3/iqtree-2.1.3-Linux.tar.gz
tar -zxvf iqtree-2.1.3-Linux.tar.gz
echo 'export PATH="~/iqtree-2.1.3-Linux/bin:$PATH"'>>~/.bashrc
source .bashrc
cd iqtree-2.1.3-Linux #该文件夹下有一个PHYLIP格式的示例比对文件example.phy

示例:

1. 根据ModelFinder自动选择的最适合模型,从序列比对(example.phy)中推断出最大似然树:
    iqtree2 -s example.phy #s指定 IQ-TREE 需要的比对文件的名称
    iqtree2 -s example.phy -m MFP 
    (# change -m MFP to -m TEST to resemble jModelTest/ProtTest)

-m是指定要在分析期间使用的模型名称的选项。特殊MFP关键字代表ModelFinder Plus,它告诉IQ-TREE使用所选模型执行ModelFinder和剩余分析。ModelFinder计算许多不同模型和Akaike信息准则(AIC)、校正Akaike信息准则(AICc)和贝叶斯信息准则(BIC)的初始简约树的对数似然。然后ModelFinder选择最小化BIC分数的模型(您也可以通过添加选项-AIC或-AICc分别更改为AIC或AICc)。当使用了ModelFinder模型时,IQ-TREE将写入一个附加文件:example.phy.model:所有测试模型的对数似然。它用作检查点文件以恢复中断的模型选择。如果您现在查看,example.phy.iqtree您将看到IQ-TREE被选TIM2+I+G4为此示例数据的最佳拟合模型。

正常情况下,在运行结束时IQ-TREE将写入几个输出文件,包括:
example.phy.iqtree:自读的主要报告文件。您可以查看此文件以查看计算结果。它还包含最终树的文本表示(见下文)。
example.phy.treefile:NEWICK格式的ML树,可以通过任何支持的树查看器程序(如FigTree或iTOL)进行可视化。
example.phy.log:整个运行的日志文件(也打印在屏幕上)。要报告错误,请将此日志文件和原始比对文件发送给作者。
对于此示例数据,生成的最大似然树可能如下所示(从.iqtree文件中提取):

NOTE: Tree is UNROOTED although outgroup taxon 'LngfishAu' is drawn at root.
注意:树是无根的,尽管外群分类单元'LngfishAu'在图中以根部出现。

+--------------LngfishAu
|
|        +--------------LngfishSA
+--------|
|        +--------------LngfishAf
|
|      +-------------------Frog
+------|
       |               +-----------------Turtle
       |         +-----|
       |         |     |      +-----------------------Sphenodon
       |         |     |   +--|
       |         |     |   |  +--------------------------Lizard
       |         |     +---|
       |         |         |      +---------------------Crocodile
       |         |         +------|
       |         |                +------------------Bird
       +---------|
                 |                  +----------------Human
                 |               +--|
                 |               |  |  +--------Seal
                 |               |  +--|
                 |               |     |   +-------Cow
                 |               |     +---|
                 |               |         +---------Whale
                 |          +----|
                 |          |    |         +------Mouse
                 |          |    +---------|
                 |          |              +--------Rat
                 +----------|
                            |   +----------------Platypus
                            +---|
                                +-------------Opossum

这很容易理解,哺乳动物(Human to Opossum)形成一个进化枝,而爬行动物(Turtle to Crocodile)和Bird形成一个单独的姐妹进化枝。这里的树画的是Lngfish的外群,在这个例子中它比其他物种更古老。然而,请注意IQ-TREE总是产生一棵无根的树,因为它对这个生物学背景一无所知;IQ-TREE简单地以这种方式绘制树,因为LngfishAu仅仅只是比对中出现的第一个序列而已。

其他示例:

2. 重新运行并且覆盖先前的输出
    iqtree -s example.phy -redo
    (Use `-redo` option if you really want to redo the analysis and overwrite all output files.)

3.更改输出前缀:
    iqtree -s example.phy --prefix myprefix
    输出文件的默认前缀是比对文件名,自定义前缀可以防止在同一文件夹内对同一路线执行多个分析时覆盖原先的输出文件。

4. 执行ModelFinder而不进行后续的树推理:
    iqtree2 -s example.phy -m MF
    (use '-m TEST' to resemble jModelTest/ProtTest)

5. 使用GTR+I+G模型推断最大似然树:
    iqtree2 -s example.phy -m GTR+I+G

6. 使用 4 个 CPU 内核来加速计算:
    iqtree2 -s example.phy -T 4

7. 确定在GTR+R4模型下使用的最佳内核数:
    iqtree2 -s example.phy -m GTR+R4 -T AUTO

8. 根据计算系统,可能需要设置可以自动分配的CPU内核的上限。使用该-ntmax选项执行此操作:
    iqtree2 -s example.phy -m TIM2+I+G -T AUTO -ntmax 8

9. 使用超快引导程序近似评估分支支持(UFBoot):
    iqtree2 -s example.phy -m TIM2+I+G -B 1000
    (-B指定引导复制的数量,其中1000是推荐的最小数量。)

10. 使用 UFBoot 减少模型严重违规的影响:
    iqtree2 -s example.phy -m TIM2+I+G -B 1000 -bnni
    -bnni降低由于严重模型违规而高估 UFBoot 分支支持的风险。
    使用此选项,UFBoot 将直接基于相应的引导程序比对使用爬山最近邻交换 (NNI) 搜索进一步优化每个引导程序树

11.使用单分支测试评估分支支持(SH-aLRT ):
    iqtree2 -s example.phy -m TIM2+I+G -alrt 1000
    -alrt指定SH-aLRT的引导复制次数,其中1000是推荐的最小数量。

12. 使用标准非参数引导程序评估分支支持:
    iqtree2 -s example.phy -m TIM2+I+G -b 100
    -b指定引导复制的数量,其中100是推荐的最小数量。输出文件类似于UFBoot程序生成的文件。

13. 在一次运行中结合ModelFinder,树搜索,超快bootstrap和SH-aLRT测试:
     iqtree2 -s example.phy --alrt 1000 -B 1000

#执行边连接的比例划分模型(example.nex):
    iqtree2 -s example.phy -p example.nex
    (replace '-p' by '-Q' for edge-unlinked model)
#通过合并分区找到最佳的分区方案:
    iqtree2 -s example.phy -p example.nex -m MF+MERGE
    (use '-m TESTMERGEONLY' to resemble PartitionFinder)
#找到最优分区方案,然后进行树推理和bootstrap:
    iqtree2 -s example.phy -p example.nex -m MFP+MERGE -B 1000

暂时先更到这里。。。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,684评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,143评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,214评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,788评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,796评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,665评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,027评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,679评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 41,346评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,664评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,766评论 1 331
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,412评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,015评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,974评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,203评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,073评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,501评论 2 343

推荐阅读更多精彩内容