Perl版Circos-基础篇-核型和刻度、标签的绘制

在绘制之前,先向大家介绍一下circos绘图的基本组成,circos绘图组成示意图如下:

circos.conf主配置文件主要包含两方面的信息:(1) 配置文件block;(2)绘图文件block。circos.conf 是(1) 和(2)两部分组成的,这个除必要block之外,其余均可根据自身需要进行个性化配置。

 

今天就展示一下核型、刻度和标签的绘制。

(一)conf 文件的配置

CIRCOS.CONF

karyotype = data/karyotype/karyotype.human.txt #核型文件的导入

chromosomes_units = 1000000 #展示刻度的最小长度单元

<<include ideogram.conf>> #调用核型文件。

<<include ticks.conf>> #调用ticks.conf文件。

<image>

<<include etc/image.conf>> #调用circos安装目录下etc/image.conf。

</image>

<<include etc/colors_fonts_patterns.conf>> #调用circos安装目录下etc/colors_fonts_patterns.conf。

<<include etc/housekeeping.conf>> #调用circos安装目录下etc/housekeeping.conf。

circos.conf 文件中调用两个本地配置文件,其余未更改的配置在安装目录中搜索。

IDEOGRAM.CONF

<ideogram>

<spacing>

default = 0.005r # 设定染色体间的距离,设置圈图中染色体之间的空隙大小,以下设置为每个空隙大小为周长的 0.5%。

</spacing>

# Ideogram position, fill and outline

radius           = 0.90r # 设定ideograms在图离圆心的 90% 处。

thickness        = 20p #设定ideograms 的厚度为20p(像素)

fill             = yes #设定 ideograms 填充颜色,填充的颜色为核型文件最后一列。

stroke_color     = dgrey #设定轮廓颜色为dgrey。

stroke_thickness = 2p #设定轮廓厚度为2p。

# Minimum definition for ideogram labels.

show_label       = yes #设定是否显示 label(chr/scaffold/gene) 。 label对应着核型文件的第 4 列,即你展示出来的名称。如果其值为 yes,则必须要有 label_radius 参数来设定 label 的位置,否则不能生成结果。

label_font       = default #设定label字体,这里为默认字体。

label_radius     = 1r+60p #设定label位置,也可采用dims函数设定。

label_size       = 30 #设定label的大小。

label_parallel   = yes #设置label的浏览方向。

</ideogram>

TICKS.CONF

#全局设置

show_ticks          = yes #设定显示刻度线。

show_tick_labels    = yes #设定显示刻度线上的标签。

<ticks>

radius           = 1r #设定刻度线显示的位置。

color            = black #设定刻度线为黑色。

thickness        = 2p #设定刻度线粗细为2p。

# the tick label is derived by multiplying the tick position# by ‘multiplier’ and casting it in ‘format’:## sprintf(format,position*multiplier)#

 

multiplier       = 1e-6 #设定显示的比例,label=位置实际大小*multiplier。

format           = %d #设定格式化label的显示,用法与perl sprintf 相同。

<tick> #定义小刻度

spacing        = 5u #设定刻度线间隔5个最小单位。

size           = 10p #设定刻度线粗细10p。

</tick>

<tick> #定义大刻度

spacing        = 25u #设定刻度线间隔25个最小单位。

size           = 15p #设定刻度线粗细15p

show_label     = yes #设定显示刻度线对应的label。

label_size     = 20p #设定刻度线对应的label的大小。

label_offset   = 10p #设定刻度线和label之间的间隔为10p。

format         = %d #设定格式化label的显示,用法与perl sprintf 相同。

</tick>

</ticks>

(二)核型文件

文件位置:/opt/circos/circos-0.69-6/example/data/karyotype/karyotype.human.txt

类别 染色体 类别名称 显示名称 起始位置 终止位置 颜色

 

(三)运行命令

/opt/circos/circos-0.69-6/bin/circos -conf  circos.conf

(四)运行结果

 

 

初次运行的小伙伴,可以将文中配置文件的参数进行调试来加深印象,那么今天的初级核型作图就结束了,如有疑问欢迎留言。

统计与绘图

Perl版Circos-基础篇-colors、fonts和patterns的配置

2020-8-28 0:35:29

统计与绘图

Perl版Circos-基础篇-图形文件结构和文件格式

2020-8-28 0:37:10

加入Q群
0 条回复 A文章作者 M管理员
    暂无讨论,说说你的看法吧
个人中心
今日签到
有新私信 私信列表
搜索