circos 可視化手冊-ideogram 篇

circos 主要用于展示染色體上的相關(guān)數(shù)據(jù),根據(jù)在染色上的位置進(jìn)行不同方式的可視化。

首先我們需要一個染色體的位置文件藐鹤。在circos中致燥,染色體的位置保存在一個文件當(dāng)中登疗,通過karyotype參數(shù)進(jìn)行設(shè)置,比如

染色體的位置保存在文件karyotype.human.txt中嫌蚤,部分內(nèi)容如下:

采用\t分隔辐益,共7列內(nèi)容。前兩列內(nèi)容是規(guī)定的脱吱,都是chr\t-,表示這部分內(nèi)容是定義染色體相關(guān)信息的智政,第3列是染色體的ID, ID必須是唯一的,用于區(qū)分不同的染色體箱蝠,第四列代表染色體的名字续捂,這個名字會顯示在最終生成的圖片中;第五列和第六列分別代表起始位置和終止位置宦搬,這里的長度都是染色體的總長度牙瓢,最后一列代表的是染色體的顏色,只不過采用了chr前綴來表示顏色间校,注意不要和染色體的ID和name 搞混淆了矾克。

circos中,內(nèi)置了許多的顏色撇簿,相關(guān)的配置保存在軟件安裝目錄的etc/colors.conf文件中聂渊,其中etc/colors.ucsc.conf文件中采用RGB標(biāo)準(zhǔn),定義了chr的各種顏色值四瘫。

從這里也能夠看出汉嗽,上面的`karyotype.human.txt``文件中抬驴,顏色是如何通過最后一列的內(nèi)容定義的等舔。

染色體的基本參數(shù)在ideogram中進(jìn)行定義。所有的參數(shù)可以分成以下3大類別

  1. 染色體的位置

  2. 染色體的顏色

  3. 染色體的標(biāo)簽

染色體的位置

circos提供了一個圓形的整體布局狡忙,用于展示染色體的數(shù)據(jù)洗做,每條染色體可以看做圓上的一段弧弓叛。通過染色體到圓心的距離來定義染色體的位置。具體的就是通過radius參數(shù)進(jìn)行定義诚纸。

circos中撰筷,關(guān)于這個參數(shù)的值,提供了兩種定義的方式畦徘;

  1. 絕對值定義
    絕對值的是通過像素定義毕籽,對應(yīng)的后綴為p,代表pixels抬闯。比如radius ?= 1000p

  2. 相對值定義
    相對值對應(yīng)的后綴為r,代表relative的意思。在etc/image.generic.conf文件中关筒,定義了參照的radius

radius of inscribed circle in image
radius ? ? ? ? = 1500p

所以我們在配置文件中定義的radius ?= 0.80r溶握, 實際等于0.8 * 1500 = 1200 像素。

染色體之間的間距通過spacing 這個block 進(jìn)行定義蒸播,default參數(shù)設(shè)置所有染色體之間的默認(rèn)距離

如果你希望改變某兩條染色體之間的距離睡榆,可以通過pairwise這個block, 用法如下

通過染色體的ID 指定具體的兩條染色體,之間用;分隔袍榆,然后通過spacing參數(shù)進(jìn)行定義它們之間的距離胀屿。

需要注意的是,這里采用的相對值的定義方法蜡塌,上面例子中的 20r代表的是相對default是20倍的距離碉纳,所以使用相對值時,一定要理解相對的參照是哪一個馏艾。

染色體的顏色

通過thickness定義染色體的寬度,單位和radius類似奴愉,有rp兩種琅摩,通常情況下我們使用p,更加方便。比如thickness = 20p锭硼。

對于染色體的顯示房资,有以下兩種方式

  1. 線段模式
    這個模式是默認(rèn)模式。在該模式下檀头,染色體就當(dāng)做一條線段轰异,顏色采用karyotype中的最后一列。

  2. 矩形模式
    在該模式下暑始,每條染色體被當(dāng)做一個可以填充的矩形搭独,所以會有邊框和填充兩部分構(gòu)成。

stroke_thicknessstroke_color 定義邊框?qū)挾群皖伾染怠H绻M刑畛渖栏危枰O(shè)置fill=yes, 此時會用karyotype中的最后一列定義的顏色填充這個矩形。

染色體的標(biāo)簽

karyotype中嗤朴,我們定義了染色體的名字配椭,也叫做label, 關(guān)于label 有對應(yīng)的一系列參數(shù),調(diào)整它的顯示雹姊。

默認(rèn)情況下股缸,這個名字并不會顯示在圖上;可以通過show_label進(jìn)行顯示吱雏,比如show_label = yes敦姻。

label 通過文字的形式展現(xiàn)在圖上瘾境,label_font設(shè)置字體,label_radius設(shè)置位置替劈,label_size設(shè)置大小寄雀, label_parallel設(shè)置label的方向,yes代表和染色體方向平行陨献,no代表和染色體方向垂直盒犹。

一個完整的ideogram例子如下:


?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市眨业,隨后出現(xiàn)的幾起案子急膀,更是在濱河造成了極大的恐慌,老刑警劉巖龄捡,帶你破解...
    沈念sama閱讀 216,372評論 6 498
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件卓嫂,死亡現(xiàn)場離奇詭異,居然都是意外死亡聘殖,警方通過查閱死者的電腦和手機晨雳,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,368評論 3 392
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來奸腺,“玉大人餐禁,你說我怎么就攤上這事⊥徽眨” “怎么了帮非?”我有些...
    開封第一講書人閱讀 162,415評論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長讹蘑。 經(jīng)常有香客問我末盔,道長,這世上最難降的妖魔是什么座慰? 我笑而不...
    開封第一講書人閱讀 58,157評論 1 292
  • 正文 為了忘掉前任陨舱,我火速辦了婚禮,結(jié)果婚禮上角骤,老公的妹妹穿的比我還像新娘隅忿。我一直安慰自己,他們只是感情好邦尊,可當(dāng)我...
    茶點故事閱讀 67,171評論 6 388
  • 文/花漫 我一把揭開白布背桐。 她就那樣靜靜地躺著,像睡著了一般蝉揍。 火紅的嫁衣襯著肌膚如雪链峭。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,125評論 1 297
  • 那天又沾,我揣著相機與錄音弊仪,去河邊找鬼熙卡。 笑死,一個胖子當(dāng)著我的面吹牛励饵,可吹牛的內(nèi)容都是我干的驳癌。 我是一名探鬼主播,決...
    沈念sama閱讀 40,028評論 3 417
  • 文/蒼蘭香墨 我猛地睜開眼役听,長吁一口氣:“原來是場噩夢啊……” “哼颓鲜!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起典予,我...
    開封第一講書人閱讀 38,887評論 0 274
  • 序言:老撾萬榮一對情侶失蹤甜滨,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后瘤袖,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體衣摩,經(jīng)...
    沈念sama閱讀 45,310評論 1 310
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,533評論 2 332
  • 正文 我和宋清朗相戀三年捂敌,在試婚紗的時候發(fā)現(xiàn)自己被綠了艾扮。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,690評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡占婉,死狀恐怖栏渺,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情锐涯,我是刑警寧澤,帶...
    沈念sama閱讀 35,411評論 5 343
  • 正文 年R本政府宣布填物,位于F島的核電站纹腌,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏滞磺。R本人自食惡果不足惜升薯,卻給世界環(huán)境...
    茶點故事閱讀 41,004評論 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望击困。 院中可真熱鬧涎劈,春花似錦、人聲如沸阅茶。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,659評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽脸哀。三九已至蹦浦,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間撞蜂,已是汗流浹背盲镶。 一陣腳步聲響...
    開封第一講書人閱讀 32,812評論 1 268
  • 我被黑心中介騙來泰國打工侥袜, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人溉贿。 一個月前我還...
    沈念sama閱讀 47,693評論 2 368
  • 正文 我出身青樓枫吧,卻偏偏與公主長得像,于是被迫代替她去往敵國和親宇色。 傳聞我的和親對象是個殘疾皇子九杂,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,577評論 2 353

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