基本配置做好，接下来就是将各种各样的数据整合进circos环图里了，本章主要介绍links、highlights、text数据作图。
Links
Links主要试图利用连线关系，表明染色体上不同位置之间的联系。这在实际生物学问题中很常用，因为基因组往往在功能调控、空间结构上存在相互作用关系，而网络图就是最佳的展示方式之一。
我们可以利用<<include>>单独引入links配置文件，也可以在主配置文件circos.conf中添加<links>模块来实现这一功能：
show_links= yes<links># 调用show_links的值show= conf(show_links)# Circos提供两种link方式，一种是bezier适合点对点连线，一种是ribbon适合片对片链接ribbon= yes# ribbon的画法是start1->end1->end2->start2，故可能出现twist，设置flat = yes取消twist# 同样地，可以设置twist = yes强行制造twistflat= yes# links的起始终止位置radius= 0.45r-55p# Circos提供bezier_radius、crest、bezier_radius_purity等参数来控制link的弯曲形式，具体见下图# 若希望设置为undef，建议在<link>子模块中设置bezier_radius= 0rcolor= black_a5# 对于ribbon，还可以通过stroke_color和stroke_thickness来设置边缘<link># 指定数据文件位置file = data/links.txt# 这里不额外讲rules的使用方法，举些例子供参考吧<rules><rule> # 只有link数据才存在chr、size等前缀+1、2的情况# 表示当size1（link起始结点）的大小小于1Mb时condition = var(size1) < 1e6# 符合上述条件的link线不显示show = no</rule># 注意这里，Circos中是从前至后依次进行rule判断的，若满足第一个rule，则不继续判断后续rule了# 可通过flow参数改变这一规则<rule># 这里相当于总为truecondition = 1# eval表示赋值# 这里即把chr2（link末端结点）的染色体名称去掉前两个字母后，组合进"chr%s_a3"，再赋给colorcolor = eval(sprintf("chr%s_a3",substr(var(chr2),2)))# remap表示取值范围映射，这里即把size2（link末端结点）的值由[0,10e6]映射至[0,-100]，然后赋给zz = eval(remap_int(var(size2),0,10e6,0,-100))</rule></rules></link></links>
这里看一下官网提供的关于曲线弯曲设置的示意图：
Circos很有意思地提供了perturb等参数，它们在指定范围内随机取值，从而可以让不同<link>子模块随机呈现不同曲线效果：
# 设置是否使用perturbperturb= yes# 设为0表示不使用perturbperturb_crest= 0# new_value = value * [ pmin + (pmax-pmin)*urd ]，其中urd在[0，1)间平均分布perturb_bezier_radius= 0.5,1.2perturb_bezier_radius_purity= 0.5,1
接下来，我们看一下links模块的展示效果：
Highlights
高亮显示是一个很常用的功能，在Circos中，你既可以用它来标注感兴趣的不同区域，也可以仅仅用来装饰整个环图。Highlight包括两种使用方式，一种是利用<highlights>模块，再一种是设置type = highlight的<plot>子模块。从图形展示的角度来说区别不大，唯一就是<highlights>模块所呈现的数据在其它数据层的下面，<plot>模块则反之。不过<highlights>模块还有个方便之处，是可以更好地设置在染色体上显示（利用ideogram = yes，而不需设置r0和r1，从而显示得更完美）。
Highlight的基本数据格式只需要染色体名称和起始终止位置即可，也可以通过补充设置fill_color、r0、r1等覆盖默认值，达到特殊显示的需要。下面，我们设置两个type = highlight的<plot>子模块作为例子：
show_highlights= yes<plots># highlight1<plot>show= conf(show_highlights)# 表示这是一个highlight类型的plottype= highlight# 指定数据文件位置file= data/highlights.txt# highlight区块的内外半径，实际不严格区分哪个在内哪个在外r0= 1r+180pr1= 1r+181p# 边缘的颜色、粗细stroke_color= blackstroke_thickness= 1# 填充色fill_color= red</plot># highlight2<plot>show= conf(show_highlights)type= highlightfile= data/highlights.txtr0= 0.45r-50pr1= 0.45r-45p# 通过rules，设置颜色和karyotype保持一致<rules><rule># Circos使用Perl语法规则，这里即正则表达式# 设置所有人类染色体的颜色condition= var(chr) =~ /hs/fill_color= eval(lc sprintf("chr%s",substr(var(chr),2)))</rule># 设置小鼠染色体的颜色<rule>condition= var(chr) =~ /mm/fill_color= blues-5-seq-4</rule># 设置大鼠染色体的颜色<rule>condition= var(chr) =~ /rn/fill_color= reds-5-seq-4</rule></rules></plot></plots>
下面是展示效果，可以看到，这里就仅仅利用highlight模块作为装饰：
Text
基因组数据的分析和作图少不了对一些区域进行注释，所以text的显示是很重要的内容。Text是一个type = text的<plot>子模块，其基本数据格式除了要求染色体名称和起始终止位置外，还应有value一栏，即text的内容。下面例子使用了人类基因组中部分omim和cancer相关基因名称的数据。值得注意的是，通常在大量数据情况下，text不能完整显示，特别是邻近位置，所以对于text排版的设置尤为重要：
show_text= yes<plots><plot>show= conf(show_text)# 表示这是一个text类型的plottype= textfile= data/text.txtcolor= black# 通过r0和r1来设定显示半径，这里r0为内，r1为外r1= 0.73rr0= 0.415rlabel_size= 9plabel_font= condensed# label之间径方向间隔rpadding= 0r# label之间角方向间隔padding= 0r# 为邻近或者同一位置的多个label设置snuggle显示方式label_snuggle= yes# 这里的距离为角方向max_snuggle_distance= 1r# 与snuggle排列显示相关的参数，值越大，计算速度越快，但准确性下降snuggle_sampling= 1# 值越大，准确性越低snuggle_tolerance= 0.25r# 这里的link是指引label的连线，具体见下图show_links= yeslink_dims= 4p,2p,5p,2p,2plink_thickness= 2plink_color= grey<rules><rule># 为omim相关基因设置颜色condition= var(value) =~ /omim/color= blue# 继续判断后续ruleflow= continue</rule><rule># 为cancer相关基因设置颜色condition= var(value) =~ /cancer/color= redlabel_font= boldflow= continue</rule><rule>condition= 1# 匹配正则表达式value= eval(replace(var(value),qr/_.*/,""))</rule></rules></plot></plots>
我们首先看一下关于label的link参数设置：
最后看一下我们text模块的效果图：
这里有一个小技巧，text不止可以是文字，也可以是符号，如设置label_font = glyph，再将相应数据的value进行更改即可，比如可在上述配置中更改规则：
<rule>condition= var(value) =~ /cancer/label_font= glyphvalue= N</rule>
原文链接https://wenlongshen.github.io/2017/03/26/Circos-3/
https://blog.sciencenet.cn/blog-543513-1042282.html
上一篇：
Circos基础（2）下一篇：
Circos基础（4）
