ggtree实现系统发育树可视化

背景1:Newick格式

Newick是最常用的存储进化树的文件格式,如上面这个树,拓朴结构用newick格式可以表示为:

(B,(A,C,E),D);

括号最外层是根节点,它有三个子节点,B, (A,C,E)和D,而节点(A,C,E)也有三个子节点A,C和E。

加上分支长度,使用 : 来分隔:

(B:6.0,(A:5.0,C:3.0,E:4.0):5.0,D:11.0);

比如A:5.0代表的是A与其父节点的距离是5.0。

内部节点也可以有label,写在相应的括号外面,如下所示:

(B:6.0,(A:5.0,C:3.0,E:4.0)Ancestor1:5.0,D:11.0);

这是最为广泛支持的文件格式,很多进化树可视软件只支持newick格式。

背景2: 如何判断无根树有根树?

假设有下面三种形式的进化树:

[1] (a, b, c);      无根树

[2] (a, b, c):1;  有根树

[3] ((a,b),c);      有根树

如果最外层大括号内只有两个分枝,即为有根树,如[3];

如果最外层大括号内有三个或以上分枝,一般为无根树,如[1];

但是,如果大括号外存在枝长参数,如[2]中的:1,这种情况认为三个分枝以一定的枝长连接到根上,为有根树

ggtree使用

ggtree只映射节点到坐标系统中,而线条在geom_tree图层中计算并画出来。这是与其它软件最根本的不同,也是ggtree能够简单地用图层加注释信息的基础。

1.ggtree安装

if (!requireNamespace("BiocManager", quietly = TRUE))

    install.packages("BiocManager")

BiocManager::install("ggtree")

2.加载包

library("ggplot2")

library("ggtree")

3. 读取数据

tree=read.tree("D:/xxx.tre")

group_file = read.csv("D:/xx.csv",header = T, row.names =1) #读取物种group信息文件,注意要与载入的Tree文件的lable相同,否则无法加载分组颜色!

4. 添加分组信息数据至变量

groupInfo <- split(row.names(group_file), group_file$family) # 提取分组信息

tree1 <- groupOTU(tree, groupInfo)  #将分组信息赋予给树

ggtree(tree, layout = "circular")

5. 运行

ggtree(tree, layout = "circular")  #简单

ggtree(tree,layout = "circular") + geom_tiplab() + geom_point(color='firebrick') #

ggtree(tree1,layout = "circular") +geom_highlight(node = 1872,fill="blue ",alpha=0.5)+  geom_highlight(node=1081,fill=" orange", alpha=0.5)+ geom_highlight(node=1079,fill=" red ")+ geom_cladelabel(node=1872,label="Pinopsida",offset=1.2,barsize = 4,color=" blue ")+geom_cladelabel(node=1081,label="Magnoliopsida",offset=1.2,barsize = 4,color=" orange",)+geom_cladelabel(node=1079,label="Ginkgoopsida",offset=1.2,barsize = 4,color="red")  #自定义


参考:

G Yu, DK Smith, H Zhu, Y Guan, TTY Lam*.ggtree: an R package for visualization and annotation of phylogenetic treeswith their covariates and other associated data. Methods in Ecology andEvolution. 2017, 8(1):28-36.

使用ggtree实现进化树的可视化和注释(转载)。作者:余光创。开发过多个R/Bioconductor包,包括ChIPseeker, clusterProfiler, DOSE,ggtree,GOSemSim和ReactomePA. 链接:https://guangchuangyu.github.io/cn/2016/03/ggtree-for-tree-visualization-annotation/ 

怎样判断newick格式的进化树是有根树,还是无根树?http://blog.sciencenet.cn/blog-656335-926722.html

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

推荐阅读更多精彩内容