使用ggtree对进化树进行分组美化

最近,小编的老师想要在一张进化树上标出不同来源的材料,平时这种“累活”小编都是使用iTOL、evolview等工具进行手动美化的,但是这次进化树里的样本足足有800多个,想要清楚、准确的标完所有的颜色,势必要耗费大量的时间,在这危急关头,小编发现了传说中Y叔写的“ggtree"!拯救小编于水火之中!

ggtree可以像ggplot2一样,用图层化的语法绘制进化树,通过不同的图层组合即可绘制出更为复杂的进化树。

首先,我们要安装好“ggtree"。

source("https://bioconductor.org/biocLite.R")
biocLite("ggtree")

接着,我们需要准备好树文件和分组文件。

树文件就是我们通过Treebest、PAML、RAxML等工具获得的结果。

分组文件的格式如下:

准备好以后,我们就可以开始美化啦!

# 加载R包
library("ggtree")
# 读取树文件(ggtree针对不同工具生成的树文件有不同的函数进行读取,如没有对应的函数,可以使用通用的办法读取)
tree <- read.tree("file.tree")
# 读取分组信息
group_file <- read.table("group_file.txt",header = T,row.names = 1)
# 按类分组
groupInfo <- split(row.names(group_file), group_file$Group)
# 将分组信息添加到树中
tree <- groupOTU(tree, groupInfo)
# 绘制进化树
ggtree(tree, layout="fan", ladderize = FALSE, branch.length = "none",aes(color=group)) + geom_tiplab2(size=3) + theme(legend.position = "right")


美美的!分好组的进化树就画好啦!


当然,除了圆形进化树,我们还可以画各种形状的进化树,具体参数大家可以去ggtree的文档里寻找!

ggtree里还有各种函数去完成各种任务,剩下的就需要大家自己去探索啦!

另外:作者开通了知乎专栏 “小王子的生信笔记” ,欢迎大家关注!共同学习!共同进步!


参考资料:

https://www.bioconductor.org/packages/release/bioc/vignettes/ggtree/inst/doc/ggtree.html

https://guangchuangyu.github.io/cn/2016/03/ggtree-for-tree-visualization-annotation/

http://blog.sciencenet.cn/home.php?mod=space&uid=255662&do=blog&id=969228

https://shengxin.ren/article/190

http://blog.sina.com.cn/s/blog_83f77c940102vpyt.html


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

推荐阅读更多精彩内容

  • Swift1> Swift和OC的区别1.1> Swift没有地址/指针的概念1.2> 泛型1.3> 类型严谨 对...
    cosWriter阅读 11,093评论 1 32
  • 参考文档:https://pan.baidu.com/s/1HaQRu8c8bNfKTSbxF5__Sw相同内容网...
    liboxiang阅读 457评论 0 1
  • 踏雪寻芳未见梅 茶花开放埋雪堆 南天竹籽红映雪 藤攀菇蹲古树肥 2018-12-30.2019-01-05
    微风LG阅读 275评论 2 12
  • 秋日朗朗秋气爽,秋风瑟瑟秋叶黄。 秋花脉脉争秋艳,秋色萧萧映秋江。 秋鸟翙翙乱秋霞,秋草萋萋染秋霜。 秋果累累秋收...
    周万亮阅读 267评论 0 1
  • 你思或不思 记忆时不时浮现 你笑也好哭也罢 脸上都会有呈现 你想与不想 时间总是在流淌 皱纹在你的不经意间刻下 爱...
    思念随风而去阅读 510评论 1 17