circos 可视化手册-ideogram 篇

circos 主要用于展示染色体上的相关数据,根据在染色上的位置进行不同方式的可视化。

首先我们需要一个染色体的位置文件。在circos中,染色体的位置保存在一个文件当中,通过karyotype参数进行设置,比如

染色体的位置保存在文件karyotype.human.txt中,部分内容如下:

采用\t分隔,共7列内容。前两列内容是规定的,都是chr\t-,表示这部分内容是定义染色体相关信息的,第3列是染色体的ID, ID必须是唯一的,用于区分不同的染色体,第四列代表染色体的名字,这个名字会显示在最终生成的图片中;第五列和第六列分别代表起始位置和终止位置,这里的长度都是染色体的总长度,最后一列代表的是染色体的颜色,只不过采用了chr前缀来表示颜色,注意不要和染色体的ID和name 搞混淆了。

circos中,内置了许多的颜色,相关的配置保存在软件安装目录的etc/colors.conf文件中,其中etc/colors.ucsc.conf文件中采用RGB标准,定义了chr的各种颜色值。

从这里也能够看出,上面的`karyotype.human.txt``文件中,颜色是如何通过最后一列的内容定义的。

染色体的基本参数在ideogram中进行定义。所有的参数可以分成以下3大类别

  1. 染色体的位置

  2. 染色体的颜色

  3. 染色体的标签

染色体的位置

circos提供了一个圆形的整体布局,用于展示染色体的数据,每条染色体可以看做圆上的一段弧。通过染色体到圆心的距离来定义染色体的位置。具体的就是通过radius参数进行定义。

circos中,关于这个参数的值,提供了两种定义的方式;

  1. 绝对值定义
    绝对值的是通过像素定义,对应的后缀为p,代表pixels。比如radius  = 1000p

  2. 相对值定义
    相对值对应的后缀为r,代表relative的意思。在etc/image.generic.conf文件中,定义了参照的radius

radius of inscribed circle in image
radius         = 1500p

所以我们在配置文件中定义的radius  = 0.80r, 实际等于0.8 * 1500 = 1200 像素。

染色体之间的间距通过spacing 这个block 进行定义,default参数设置所有染色体之间的默认距离

如果你希望改变某两条染色体之间的距离,可以通过pairwise这个block, 用法如下

通过染色体的ID 指定具体的两条染色体,之间用;分隔,然后通过spacing参数进行定义它们之间的距离。

需要注意的是,这里采用的相对值的定义方法,上面例子中的 20r代表的是相对default是20倍的距离,所以使用相对值时,一定要理解相对的参照是哪一个。

染色体的颜色

通过thickness定义染色体的宽度,单位和radius类似,有rp两种,通常情况下我们使用p,更加方便。比如thickness = 20p

对于染色体的显示,有以下两种方式

  1. 线段模式
    这个模式是默认模式。在该模式下,染色体就当做一条线段,颜色采用karyotype中的最后一列。

  2. 矩形模式
    在该模式下,每条染色体被当做一个可以填充的矩形,所以会有边框和填充两部分构成。

stroke_thicknessstroke_color 定义边框宽度和颜色。如果希望有填充色,需要设置fill=yes, 此时会用karyotype中的最后一列定义的颜色填充这个矩形。

染色体的标签

karyotype中,我们定义了染色体的名字,也叫做label, 关于label 有对应的一系列参数,调整它的显示。

默认情况下,这个名字并不会显示在图上;可以通过show_label进行显示,比如show_label = yes

label 通过文字的形式展现在图上,label_font设置字体,label_radius设置位置,label_size设置大小, label_parallel设置label的方向,yes代表和染色体方向平行,no代表和染色体方向垂直。

一个完整的ideogram例子如下:


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

推荐阅读更多精彩内容