进化树的可视化软件非常多,其中R包 ggtree 功能非常强大,非常灵活,简单记录自己的学习笔记

第一步:使用 mafft 比对氨基酸序列

mafft --auto ggtree_pra_2.fasta > ggtree_pra_aligned.fasta

第二步:使用 FastTree 构建最大似然树

./FastTree ggtree_pra_aligned.fasta > ggtree_pra.tree

(以最简单的方式完成建树过程,因为主要目的是获得树文件来进行ggtree的学习,实际的建树过程应该比这个复杂的多,自己还有好多地方需要进一步的学习)

getwd()#查看当前工作路径

setwd("../Python_practice/a_kind_of/")#更改至树文件所在路径

library(ggtree)#加载ggtree

tree<-read.tree("ggtree_pra.tree")#读入树文件

ggtree(tree)#绘制进化树

(以上就是最基本的步骤)

接下来就是美化进化树需要用到的参数,比如用color,size,linetype调整线的颜色,粗细,线型等(和ggplot2很像或者说一样),可以加上branch.length="none"让枝的末端齐平,用layout=""参数更改树的形状,最常用的应该是layout="circular",更多地参数可用通过help(package="ggtree")命令来查看帮助文档,帮助文档非常详细;可以通过 geom_tiplab() 函数添加用于构树序列的名称(应该有专业名词,但自己想不起来叫啥了),geom_tiplab()函数可以用size参数来调节字体的大小,用hjust参数调节文字与枝末端的距离;绘制圆形树的时候需要加上aes(angle=angle)参数改变文字的角度。大家可能会发现还有一个函数是 geom_tiplab2(), 它和geom_tiplab()的主要区别就是绘制圆形的进化树时文字的方向; 有时用Rstudio查看绘图结果时可能会显示不完全,这时可以用 xlim() 函数来调节图片右边的空白区域大小。可以比较一下以下命令出图的差别:

ggtree(tree)+geom_tiplab()

ggtree(tree,branch.length="none")+geom_tiplab()

ggtree(tree,branch.length="none")+geom_tiplab(hjust=-0.05)

ggtree(tree,branch.length="none")+geom_tiplab()+xlim(NA,18)#xlim()函数里数值的大小可能得根据自己的树文件来确定

ggtree(tree,layout="circular")

ggtree(tree,layout="circular")+geom_tiplab()

ggtree(tree,layout="circular")+geom_tiplab(aes(angle=angle))

ggtree(tree,layout="circular",branch.length="none")+geom_tiplab(aes(angle=angle))

ggtree(tree,layout="circular",branch.length="none")+geom_tiplab2(aes(angle=angle))

接下来是给按不同的分组添加不同的颜色

首先是在树上显示出节点的位置

p1<-ggtree(tree,branch.length="none")+geom_tiplab(hjust=-0.1)+xlim(NA,18)

p1+geom_text2(aes(label=node))

p1+geom_text2(aes(subset=!isTip,label=node))#注意结果图的区别

使用geom_text2()函数添加 bootstrap values

p1+geom_text2(aes(subset=!isTip,label=support))#添加bootstrap values在最开始读入树文件的时候可能得该用read.newick()函数;tree<-read.newick("input.tree",node.label="support")

p1+geom_hilight(node=51,fill="red")+geom_hilight(node=43,fill="blue")#geom_hilight()函数添加色块

p1+geom_cladelabel(node=51,label="label one")#geom_cladelabel()函数在枝的末端添加线段和文字,node参数设置线段的位置,label参数设置的是添加文字的内容,color参数改变线和字的颜色,barsize参数改变线的粗细,angle参数改变文字的角度,offset参数改变线段和枝末端的距离,hjust参数改变添加的label和线段的距离,fontsize改变label的大小

p1+geom_strip(39,49,barsize=4,color="red")#geom_strip()函数在节点之间添加线段,参数和geom_cladelabel()很多一样的地方,alpha参数用来设置线段的透明度

p1+geom_point2(aes(subset=!isTip))#geom_point2()函数可以在节点上添加点,可以通过shape和color等参数改变点的形状和颜色;subset=!isTip这个参数设置的是只在内部节点打点,而枝的末端不添加点

接下来重复 ggtree 的作者写的一个教程 ​​bootstrap value 分段标记​

原教程贴出了完整的代码,但是没有给出树文件,贴出自己用到的树文件百度云链接​​密码:vs4p​

(本次用到的文件与原教程的bootstrap value范围不一样,本次用到的树文件的bootstrap value在0:1所以用cut函数分割的时候需要注意)

(遇到的问题,在一些论文中,进化树树上只添加bootstrap value大于75或者80的节点,这一点用ggtree如何实现暂时还没有找到解决方案)

ggtree笔记_bootstrap 重复的结果(原谅自己的审美观)



ggtree笔记_帮助文档_02 重复代码


(问题:图例上不知为什么会有NA存在)

原文的教程将点的颜色设置成了灰白黑,看起来可能更美观!

还有在枝顶点根据物种不同添加不同的形状,今天有点困了,明天补充上!