生信分析32:deeptools軟件介紹

Fig 1

對(duì)于各種表觀組學(xué)你辣,如ChIP-seq、ATAC-seq拦盹、CUT&TAG鹃祖、MNase-seq,deeptools是下游分析中非常常用的軟件之一普舆。該軟件發(fā)表于2014年(Fig 1)恬口,目前已經(jīng)引用超1000次。

Fig 2

Fig 3

今天的推送主要介紹deeptools軟件的兩個(gè)主要使用場(chǎng)景沼侣,數(shù)據(jù)質(zhì)控-相關(guān)性分析(Fig 2)和聚類分析(Fig 3)祖能。

01 輸入文件的準(zhǔn)備

Fig 4

對(duì)于比對(duì)的sam文件,我們轉(zhuǎn)為bam文件并完成基本過濾和排序华临,然后轉(zhuǎn)成bw文件芯杀。

Deeptools內(nèi)置了bamCoverage功能,可以完成bam2bw的操作(Fig 4)雅潭。

--binSize Size of the bins, in bases, for the output of the bigwig/bedgraph file. (Default: 50)

--normalizeUsing {RPKM,CPM,BPM,RPGC,None}

02??樣本相關(guān)性分析

Fig 5

Fig 6

multiBigwigSummary選項(xiàng)可以對(duì)多個(gè)bw文件進(jìn)行整合(Fig 5)揭厚,然后結(jié)合plotCorrelation和plotPCA可以實(shí)現(xiàn)樣本相關(guān)性分析(Fig 6-10)。

Fig 7

Fig 8

Fig 9

03 聚類分析

Fig 10

Fig 11

computeMatrix可以展示指定位置附近的seq信號(hào)強(qiáng)度扶供,它有“scale-regions”和“reference-point”兩種模式筛圆,前者將指定的區(qū)間縮放到指定大小,后者將指定的區(qū)間縮放成一個(gè)點(diǎn)椿浓,實(shí)例中使用scale-regions模式太援,輸入了所有基因的位置信息,將每個(gè)基因縮放成2k扳碍,展示上下游2k的信息(Fig 10)提岔。Bed文件為bed6格式,記得加入正負(fù)鏈信息(Fig 11)

Fig 12

利用plotHeatmap繪圖(Fig 12-15)

Fig 13 test1

Fig 14 test2

Fig 15 test3

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末笋敞,一起剝皮案震驚了整個(gè)濱河市碱蒙,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖赛惩,帶你破解...
    沈念sama閱讀 206,723評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件哀墓,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡喷兼,警方通過查閱死者的電腦和手機(jī)篮绰,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,485評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來季惯,“玉大人吠各,你說我怎么就攤上這事∶阕ィ” “怎么了走孽?”我有些...
    開封第一講書人閱讀 152,998評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)琳状。 經(jīng)常有香客問我,道長(zhǎng)盒齿,這世上最難降的妖魔是什么念逞? 我笑而不...
    開封第一講書人閱讀 55,323評(píng)論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮边翁,結(jié)果婚禮上翎承,老公的妹妹穿的比我還像新娘。我一直安慰自己符匾,他們只是感情好叨咖,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,355評(píng)論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著啊胶,像睡著了一般甸各。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上焰坪,一...
    開封第一講書人閱讀 49,079評(píng)論 1 285
  • 那天趣倾,我揣著相機(jī)與錄音,去河邊找鬼某饰。 笑死儒恋,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的黔漂。 我是一名探鬼主播诫尽,決...
    沈念sama閱讀 38,389評(píng)論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼炬守!你這毒婦竟也來了牧嫉?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 37,019評(píng)論 0 259
  • 序言:老撾萬榮一對(duì)情侶失蹤劳较,失蹤者是張志新(化名)和其女友劉穎驹止,沒想到半個(gè)月后浩聋,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,519評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡臊恋,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,971評(píng)論 2 325
  • 正文 我和宋清朗相戀三年衣洁,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片抖仅。...
    茶點(diǎn)故事閱讀 38,100評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡坊夫,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出撤卢,到底是詐尸還是另有隱情环凿,我是刑警寧澤,帶...
    沈念sama閱讀 33,738評(píng)論 4 324
  • 正文 年R本政府宣布放吩,位于F島的核電站智听,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏渡紫。R本人自食惡果不足惜到推,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,293評(píng)論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望惕澎。 院中可真熱鬧莉测,春花似錦、人聲如沸唧喉。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,289評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽八孝。三九已至谬运,卻和暖如春阅悍,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,517評(píng)論 1 262
  • 我被黑心中介騙來泰國(guó)打工栏饮, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留逢唤,地道東北人迈着。 一個(gè)月前我還...
    沈念sama閱讀 45,547評(píng)論 2 354
  • 正文 我出身青樓袍祖,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親浅辙。 傳聞我的和親對(duì)象是個(gè)殘疾皇子扭弧,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,834評(píng)論 2 345

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