【三維基因組】多組學(xué)可視化之washU

說到三維基因組可視化,就不得不說一下washU(http://epigenomegateway.wustl.edu/)了。washU相當(dāng)于IGV的升級(jí)版,一向被三維基因組CNS級(jí)別的文章所鐘情,能夠產(chǎn)生極其fancy的效果,是三維基因的可視化神器。

washU 目前支持包括人,小鼠,黑猩猩,斑馬魚等一系列物種。如下圖所示。

你可以選擇你需要的基因和版本進(jìn)行可視化分析。

在選取了基因組版本號(hào)之后,你可以在搜索??框中,選擇你所關(guān)注的基因或者SNP(如下圖):

輸入關(guān)注的基因


輸入關(guān)注的snp

在選擇完了關(guān)注的位點(diǎn)之后,你就可以通過view Local Tracks 導(dǎo)入你本地的文件。

washU 支持的文件包括bigwig,bed,bam等等(如下所示)。

像bam,bigwig,bed文件可以導(dǎo)入washU,以峰圖的形式呈現(xiàn)。

那么如何導(dǎo)入呢?有兩種方式,一方面是剛剛所說的 本地上傳 ,另一方面是通過URL提供數(shù)據(jù),如下圖:

知道了數(shù)據(jù)類型和數(shù)據(jù)的上傳方式,那么我們?nèi)绾螠?zhǔn)備數(shù)據(jù)呢?

首先看bed文件(如下圖):

包括4列: chr ,start,end ,value(如下圖)

處理命令如下:


sort -k1,1 -k2,2n track.bedgraph > track.bedgraph.sorted? #對(duì)bedgraph進(jìn)行排序bgzip? track.bedgraph.sorted? # 對(duì)bedgraph進(jìn)行壓縮

tabix -p bed track.bedgraph.sorted.gz? #對(duì)壓縮后的bedgraph 建立索

bgzip安裝如下:

wget https://github.com/samtools/bcftools/releases/download/1.2/bcftools-1.2.tar.bz2?

tar xvf bcftools-1.2.tar.bz2

cd bcftools-1.2

makemake install

samtools sort test.bam > test.sorted.bam ? ? ? ? ? ? ? ? ? ?

? ?samtools ?index test.sorted.bam

bam文件有兩種顯示方式,一方面是顯示每條reads 的信息,另一方面是顯示峰圖。

bam文件必須先sort再導(dǎo)入washU。

甲基化數(shù)據(jù)也可以用washU來展示(如下圖所示)


甲基化數(shù)據(jù)格式如下:

第一到三列代表甲基化位點(diǎn)的position

第四列代表甲基化的類型

第五列代表甲基化水平

第六列代表 正負(fù)鏈

第七列代表位點(diǎn)的覆蓋深度

Hi-C數(shù)據(jù)可視化主要分為兩種,一種是loop結(jié)構(gòu),一種是三角熱圖。其中三角熱圖數(shù)據(jù)采用.hic 文件格式。

.hic 文件主要是Juicer 產(chǎn)生的文件格式。



longrange

The?longrange?track is a?bed?format-like file type. Each row contains columns from left to right:?chromosome,?start?position?(0-based), and?end?position?(not included), interaction target in this format?chr2:333-444,55. As an example, interval “chr1:111-222” interacts with interval “chr2:333-444” on a score of 55, we will use following two lines to represent this interaction:

第一列: loop 左端的染色體編號(hào)

第二列:loop左端的起始位置

第三列:loop左端的終止位置

第四列:loop右端的位置以及交互數(shù)值。

此外,washU還支持calling card track 格式如下:

那么針對(duì)以上文件總結(jié)而言:

bgzip track.bedgraph.sorted?

tabix -p bed track.bedgraph.sorted.gz


針對(duì)bedgraph文件都要用bgzip壓縮并建立索引。

那么最后舉個(gè)例子,在導(dǎo)入longrange,以及bam或者bigwig,我們基本就有了下面這張類似與文獻(xiàn)里的圖。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容