从零开始学Circos绘制圈图(一)

一般基因组文章都会有下面这种酷炫图,用来描述基因组的基因密度分布,转座子的密度分布,和其他物种或者多倍体的多套染色体间的共线性关系,以及其他各种你只要测序就能加上的信息,比如说你要是测了ATAC-seq,加上全基因组开放状态,要是测了多个组织,多个时期的RNA-seq,那就加上热图展现这种变化关系。

circos绘制基因组

当然除了基因组文章,其他类型的文章也可以考虑这种图。接下来我将会写一些列教程(可能有视频),通过教别人学Circos的方式来自学Circos。

环境配置

建议在Linux环境下配置Circos,之后只要用conda就能配置好分析环境

# 安装
## circos
conda create -c bioconda -n circos circos

测试软件安装结果

# 测试circos
conda activate circos
# 确认安装
circos -V
# 显示如下
# circos | v 0.69-8 | 15 Jun 2019 | Perl 5.026002

可以从http://circos.ca/software/download/下载官方的教程文件,分别是

处理过程

Circos依赖于一些列的配置文件,用来定义复杂图形的各个部分,最终加工成图形。

因此,用Circos画图是一个不断增添内容的过程,你要不断根据输出结果来调整输入参数。

并且整个分析中,你还要拥有过关的数据预处理的能力,这是因为Circos不是数据处理工具,它只是展示你已有的数据。

处理过程

快速开始

不管怎么样,先快速绘制出一个Circos图再说。

果子老师说过,我们不是先成为了老司机才开车,而是开车多了才成为了老司机。

第一步,先新建一个文件夹,用于存放本次分析的所有数据和配置文件

mkdir -p my_first_circos && cd my_first_circos

然后用vim karyotype.tair10.txt编辑文本,新增如下内容

chr - chr1 chr1 0 30427617 black
chr - chr2 chr2 0 19698289 black
chr - chr3 chr3 0 23459830 black
chr - chr4 chr4 0 18585056 black
chr - chr5 chr5 0 26975502 black

之后创建一个circos.conf文件,用于增加各类配置参数

touch circos.conf

vim circos.conf,增加我们的第一条记录,染色体信息

karyotype = karyotype.tair10.txt

然而要想真正的出图,还需要增加至少以下配置语句才行

<ideogram>
<spacing>
default = 0.005r
</spacing>
radius           = 0.90r
thickness        = 20p
fill             = yes
stroke_color     = dgrey
stroke_thickness = 2p
</ideogram>

<image>
<<include etc/image.conf>>
</image>

<<include etc/colors_fonts_patterns.conf>>
<<include etc/housekeeping.conf>>

在当前路径下运行circos -conf circos.conf, 最终效果图如下

第一张图

虽然图比较丑,但是至少我们成功运行了人生第一次的circos, 这就相当于买了一套毛坯房,后面要做的事情就是不断装修。

比如说,我们至少可以让不同染色体拥有不同的颜色,修改之前的karyotype.tair10.txt中的最后一列

chr - chr1 chr1 0 30427617 chr1
chr - chr2 chr2 0 19698289 chr2
chr - chr3 chr3 0 23459830 chr3
chr - chr4 chr4 0 18585056 chr4
chr - chr5 chr5 0 26975502 chr5

在当前路径下运行circos -conf circos.conf, 最终效果图如下

1564117790558

这就引出了第一个知识点,配色

为了实现配色,需要circos.conf文件了有一个和配色有关的语句

<<include etc/colors_fonts_patterns.conf>>

这里<<>>表示通过相对路径的方式加载另外一个配置文件,它的实际路径是和circos所在目录同级的etc,可用下面语句看到colors_fonts_patterns.conf的内容

circos_path=$(dirname `which circos`)
less ${circos_path%bin}/etc/colors_fonts_patterns.conf

你会发现,这个文件里还嵌套其他的配置文件。最终通过层层排查,你才知道etc/colors.ucsc.conf才是实际定义我们填写的颜色名的文件,而颜色的定义如下:

chr1  = 153,102,0
chr2  = 102,102,0
chr3  = 153,153,30
chr4  = 204,0,0
chr5  = 255,0,0

还有一个问题,为什么这里用的是两个尖括号<<,而不是一个尖括号<呢?这是因为<已经被用于分隔不同的语句块,如下语句就表示etc/image.conf里的配置信息是用来调整和image有关的配置,而不是去调整ideogram的配置。

<image>
<<include etc/image.conf>>
</image>

以上是快速开始部分,后续将会在此基础上,做出发表级别的图。

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

推荐阅读更多精彩内容

  • 一、Python简介和环境搭建以及pip的安装 4课时实验课主要内容 【Python简介】: Python 是一个...
    _小老虎_阅读 5,719评论 0 10
  • circos 是一款perl 语言开发的画图软件,提供了染色体相关数据的一种可视化方式。其制作的图表精美,被科研工...
    BeautifulSoulpy阅读 5,700评论 1 9
  • 刘小泽上篇写于18.10.3 ,中篇写于18.10.4,下篇写于18.10.6学习障碍的规律往往是:不是因为它多难...
    刘小泽阅读 2,702评论 3 13
  • circos是用来展示数据的绘图工具,其数据展示方式主要通过配置文件来决定。不同的展示方式可通过不同的配置文件来进...
    井底蛙蛙呱呱呱阅读 11,710评论 2 23
  • 欢迎关注天下博客:http://blog.genesino.com/2012/04/circos-usage/Ci...
    生信宝典阅读 879评论 0 5