circos軟件學習筆記:在圓圈里添加lines(interaction)

前幾篇筆記記錄了如何用circos畫出帶有刻度以及刻度標簽的染色體骨架笨触。這篇筆記來學習如何在這個圓圈里添加其他的Data信息演侯。首先我們跟著官網(wǎng)教程來練習添加lines顺又。官網(wǎng):here

Links代表了兩個基因組位置之間的聯(lián)系(比如相似性、不同捧灰、或者fusion等等)钻注,這個聯(lián)系是通過一條曲線或者直線來表示的蚂且。

Links可視化的格式可以使用rules來設置,rules組成一個表達式的決策樹幅恋,對每個link進行評估杏死。根據(jù)數(shù)據(jù)值測試links,根據(jù)位置捆交、大小或其他參數(shù)來動態(tài)的改變其展示方式淑翼。

Links的數(shù)據(jù)在<link> blocks里設置。Links從一個定義的徑向位置開始品追,并在bezier_radius定義的徑向位置上有它們的控制點(調(diào)整曲率)玄括。比如在我們最開始下載的示例數(shù)據(jù)文件夾data/5/segdup.txt這個文件:

$ head segdup.txt
hs1 465 30596 hs2 114046768 114076456
hs1 486 76975 hs15 100263879 100338121
hs1 486 30596 hs9 844 30515
hs1 486 9707 hsY 57762276 57771573
hs1 486 9707 hsX 154903076 154912373
hs1 486 9707 hs16 427 9533
hs1 8256 76975 hs19 11001 79672
hs1 23908 30596 hs12 17641 24322
hs1 59871 76975 hs6 5001 22075
hs1 71096 76975 hs6 170824311 170830148

舉個例子:在你設置文件的時候就應該這樣寫:

<links>

<link>
file          = data/5/segdup.txt
radius        = 0.8r 
bezier_radius = 0r
color         = black_a4
thickness     = 2
</link>

</links>

在Rule塊可以被添加到任何links block或者plot blocks里,并且形成一個決策鏈肉瓦,用來更改如何展示data遭京,是線條胃惜,柱形圖,還是散點圖等等哪雕,一般格式是:

<rules> 

<rule>
...
</rule>

<rule>
...
</rule>

...

</rules>

每一個rule有一個條件船殉、格式化語句和一個可選的flow statement。如果條件為true斯嚎,則將這個rule應用于數(shù)據(jù)點利虫,不再檢查其他規(guī)則(除非flow=continue)。如果條件為false堡僻,則檢查下一個rule糠惫。
var(X)指的是數(shù)據(jù)點的變量X的值。這里intrachr是指染色體內(nèi)的:

<rule>
condition     = var(intrachr)
# Any links that are intra-chromosomal will not be shown. Further rules are not tested.
show          = no
</rule>

rules是按照順序運行的钉疫。你可以通過向上移動block或添加importance參數(shù)硼讽,在決策樹中向上移動rules。首先是帶有importance參數(shù)的rules被運行陌选,然后依次運行不帶參數(shù)的rules:

<rules>

<rule>
# 2nd
</rule>

<rule>
# 3rd
</rule>

<rule>
# 1st 
importance = 10
</rule>

</rules>

示例腳本

circos.conf主要配置文件:

karyotype = data/karyotype/karyotype.human.txt
chromosomes_units = 1000000

chromosomes_display_default = no
chromosomes                 = /hs[1-4]$/
chromosomes_reverse         = /hs[234]/
chromosomes_scale           = hs1=0.5r,/hs[234]/=0.5rn
chromosomes_radius          = hs4:0.9r

# 在上一篇筆記里理郑,我們是在ideogram.conf更改染色體顏色的蹄溉,現(xiàn)在我們用第二種方法來更改染色體顏色:
# 星號的后綴代表覆蓋顏色咨油,可以覆蓋colors_fonts_patterns.conf文件里的設置
<colors>
chr1* = red
chr2* = orange
chr3* = green
chr4* = blue
</colors>

<links> #links的開頭

<link>
file          = data/5/segdup.txt
radius        = 0.8r #從徑向位置的0.8r處開始線條
bezier_radius = 0r
color         = black_a4
thickness     = 2

<rules># 在links塊里添加rules
<rule> #第一個rule開始
condition     = var(intrachr) #條件是:在同一染色體內(nèi)
# 任何在同一染色體內(nèi)的links都不會被展示
show          = no
</rule> #第一個rule結束

<rule> #第二個rule
#這個rule是針對剩下的不在同一染色體內(nèi)的links,condition=1代表條件始終是true
condition     = 1
#設置links的線條顏色柒爵,所有結尾落在2號染色體上的links的顏色役电,和2號染色體本身的顏色一致
color         = eval(var(chr2))
#在執(zhí)行完這一個rule后,繼續(xù)檢查下一個rule棉胀,所以設置continue:
flow          = continue
</rule>#第二個rule結束

<rule>#第三個rule
# 如果links從1號染色體開始:
condition     = from(hs1)
# 那么links的起始點要非撤ㄉ靠近染色體骨架,0.99r
radius1       = 0.99r
</rule> #第三個rule結束

<rule>#第四個rule
# 如果links的結束點在1號染色體:
condition     = to(hs1)
# 那么links的結束點要非逞渖荩靠近染色體骨架霎挟,對同一個染色體進行設置,要用radius2:
radius2       = 0.99r
</rule>#第四個rule結束

</rules>#所有的rule結束

</link> #link結束
</links> #links的結尾

<<include ideogram.conf>>
<<include ticks.conf>>

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

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

ideogram.conf腳本(與上一篇筆記里的差不多):

<ideogram>

<spacing>
default = 0.005r
</spacing>

# Ideogram position, fill and outline
radius           = 0.90r
thickness        = 20p
fill             = yes
stroke_color     = dgrey
stroke_thickness = 2p

# Minimum definition for ideogram labels.

show_label       = yes
# see etc/fonts.conf for list of font names
label_font       = default
label_radius     = 1.075r  # if ideogram radius is constant, and you'd like labels close to image edge,
                           # use the dims() function to access the size of the image
                           # label_radius  = dims(image,radius) - 60p
label_size       = 30
label_parallel   = yes

</ideogram>

ticks.conf腳本(與上一篇筆記里的差不多):

show_ticks          = yes
show_tick_labels    = yes

<ticks>
radius           = 1r
color            = black
thickness        = 2p
multiplier       = 1e-6
format           = %d

<tick>
spacing        = 5u
size           = 10p
</tick>

<tick>
spacing        = 25u
size           = 15p
show_label     = yes
label_size     = 20p
label_offset   = 10p
format         = %d
</tick>

</ticks>

運行circos.conf:

$ circos -conf ./circos.conf -outputdir ./ -outputfile tutorial_1_4_practice.png
$ xdg-open tutorial_1_4_practice.png
這里可以看到麻掸,所有連接到1號染色體的links都和染色體離的很近酥夭,而連接到2,3,4號染色體的links都和染色體有一段距離,這就是上面我們進行設置rule的時候設置的
?著作權歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
禁止轉(zhuǎn)載脊奋,如需轉(zhuǎn)載請通過簡信或評論聯(lián)系作者熬北。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市诚隙,隨后出現(xiàn)的幾起案子讶隐,更是在濱河造成了極大的恐慌,老刑警劉巖久又,帶你破解...
    沈念sama閱讀 219,490評論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件巫延,死亡現(xiàn)場離奇詭異效五,居然都是意外死亡,警方通過查閱死者的電腦和手機炉峰,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,581評論 3 395
  • 文/潘曉璐 我一進店門火俄,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人讲冠,你說我怎么就攤上這事瓜客。” “怎么了竿开?”我有些...
    開封第一講書人閱讀 165,830評論 0 356
  • 文/不壞的土叔 我叫張陵谱仪,是天一觀的道長。 經(jīng)常有香客問我否彩,道長疯攒,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,957評論 1 295
  • 正文 為了忘掉前任列荔,我火速辦了婚禮敬尺,結果婚禮上,老公的妹妹穿的比我還像新娘贴浙。我一直安慰自己砂吞,他們只是感情好,可當我...
    茶點故事閱讀 67,974評論 6 393
  • 文/花漫 我一把揭開白布崎溃。 她就那樣靜靜地躺著蜻直,像睡著了一般。 火紅的嫁衣襯著肌膚如雪袁串。 梳的紋絲不亂的頭發(fā)上概而,一...
    開封第一講書人閱讀 51,754評論 1 307
  • 那天,我揣著相機與錄音囱修,去河邊找鬼赎瑰。 笑死,一個胖子當著我的面吹牛破镰,可吹牛的內(nèi)容都是我干的餐曼。 我是一名探鬼主播,決...
    沈念sama閱讀 40,464評論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼啤咽,長吁一口氣:“原來是場噩夢啊……” “哼晋辆!你這毒婦竟也來了?” 一聲冷哼從身側響起宇整,我...
    開封第一講書人閱讀 39,357評論 0 276
  • 序言:老撾萬榮一對情侶失蹤瓶佳,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后鳞青,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體霸饲,經(jīng)...
    沈念sama閱讀 45,847評論 1 317
  • 正文 獨居荒郊野嶺守林人離奇死亡为朋,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,995評論 3 338
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了厚脉。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片习寸。...
    茶點故事閱讀 40,137評論 1 351
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖傻工,靈堂內(nèi)的尸體忽然破棺而出霞溪,到底是詐尸還是另有隱情,我是刑警寧澤中捆,帶...
    沈念sama閱讀 35,819評論 5 346
  • 正文 年R本政府宣布鸯匹,位于F島的核電站,受9級特大地震影響泄伪,放射性物質(zhì)發(fā)生泄漏殴蓬。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,482評論 3 331
  • 文/蒙蒙 一蟋滴、第九天 我趴在偏房一處隱蔽的房頂上張望染厅。 院中可真熱鬧,春花似錦津函、人聲如沸肖粮。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,023評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽尿赚。三九已至,卻和暖如春蕉堰,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背悲龟。 一陣腳步聲響...
    開封第一講書人閱讀 33,149評論 1 272
  • 我被黑心中介騙來泰國打工屋讶, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人须教。 一個月前我還...
    沈念sama閱讀 48,409評論 3 373
  • 正文 我出身青樓皿渗,卻偏偏與公主長得像,于是被迫代替她去往敵國和親轻腺。 傳聞我的和親對象是個殘疾皇子乐疆,可洞房花燭夜當晚...
    茶點故事閱讀 45,086評論 2 355

推薦閱讀更多精彩內(nèi)容