这篇文章主要介绍如何使用iTOL在进化树中添加热图,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们一定要看完!

使用的工具依然是iTOL这款在线软件,可点这里 绘一棵超酷炫的系统发育树 了解iTOL的使用方法。基础的使用方法就不再讲述了,只对配置文件进行说明。

配置文件准备

首先来看一下所需要的文件的格式:

DATASET_HEATMAP

#SEPARATOR TAB
SEPARATOR SPACE
DATASET_LABEL example_heatmap
COLOR #ff0000

#define labels for each individual field column
FIELD_LABELS A1 A2 A3 B1 B2 B3

MARGIN 0
STRIP_WIDTH 25

COLOR_MIN #ff0000
COLOR_MAX #0000ff
USE_MID_COLOR 1
COLOR_MID #ffff00

DATA

155864 1.33118766397716 1.17857919899514 -0.409551070843932 -0.630306286554169 -1.04604589954572 -0.423863606028482
83334 -0.0506076910870376 1.3902493708319 1.06059120726857 -0.886946488704782 -0.820282432690725 -0.693003965617928
217992 0.215479234585023 1.4989464081584 0.688218558910058 -1.21085219033218 -0.764324705926616 -0.427467305394685
562 1.06564649867261 0.961867445609222 0.597829776500456 -0.461509960030431 -1.19966321785905 -0.964170542892806
198215 0.96900544380798 1.21939609883181 0.426791627672821 -0.838548688929931 -0.653211840822719 -1.12343264055996
623 0.727282455309188 1.17101033530621 0.806443403628944 -0.820220618604714 -0.920091135497558 -0.964424440142074
209261 0.34017826857092 1.66855815466279 0.357643216125776 -0.545542360482732 -1.00245964906156 -0.818377629815193
601 0.693183645180168 1.46683923827496 0.279576953582569 -1.04280232644501 -1.03511634861385 -0.36168116197883
602 0.70273801067879 1.37079662629187 0.366781728164418 -0.773552511347515 -1.33593177597361 -0.330832077813954
229193 0.877365674482512 1.19322948550412 0.554384448895058 -0.958849032442236 -1.14266481591487 -0.523465760524581

格式说明:

配置文件中以"#"开头的行是注释行,只起说明作用,iTOL不会读取该行的内容。

DATASET_HEATMAP:配置文件数据类型声明,这里声明数据集是对进化树添加热图信息。

SEPARATOR SPACE:指定数据分隔符,即数据列与列之间使用空格分隔。

COLOR #ff00ff:指定数据集初始颜色,可在后面进行更改。

FIELD_LABELS A1 A2 A3 B1 B2 B3:为DATA中的每一列数据定义标签(指定样品名称)。

MARGIN 0:指定热图与进化树之间的距离,可为负数,会导致与进化树重叠。

STRIP_WIDTH 25:指定热图中每一格的宽度。

COLOR_MIN #ff0000:指定热图中最小值的颜色。

COLOR_MAX #0000ff:指定热图中最大值的颜色。

USE_MID_COLOR 1:指定是否设置中间值颜色,0代表不设置,1代表设置。

COLOR_MID #ffff00:指定热图中中间值的颜色。

如果没有特殊设置,以上部分通常不需要改,直接拷贝使用即可。

DATA 之后的行就是所有的基因在各样品中的表达量数据,这里我用的数据是基因在样品中的表达量水平log2(FPKM)进行归一化转换(scale number)后的值。数据分为七列,第一列为进化树枝节点ID(基因ID),第二到第七列为基因在各样品中的表达量数据,注意数据要与 FIELD_LABELS 所指定样品名称的数量及顺序保持一致。

注释文件美化

按照以上格式准备好注释文件后,直接将其拖到iTOL进化树展示所在的位置,就能自动在进化树中添加热图了。如果对图片的效果不满意,还也可以在Datasets中对热图的颜色、宽度及位置进行调整。