【三维基因组】多组学可视化之washU

说到三维基因组可视化,就不得不说一下washU(http://epigenomegateway.wustl.edu/)了。washU相当于IGV的升级版,一向被三维基因组CNS级别的文章所钟情,能够产生极其fancy的效果,是三维基因的可视化神器。

washU 目前支持包括人,小鼠,黑猩猩,斑马鱼等一系列物种。如下图所示。

你可以选择你需要的基因和版本进行可视化分析。

在选取了基因组版本号之后,你可以在搜索🔍框中,选择你所关注的基因或者SNP(如下图):

输入关注的基因


输入关注的snp

在选择完了关注的位点之后,你就可以通过view Local Tracks 导入你本地的文件。

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

像bam,bigwig,bed文件可以导入washU,以峰图的形式呈现。

那么如何导入呢?有两种方式,一方面是刚刚所说的 本地上传 ,另一方面是通过URL提供数据,如下图:

知道了数据类型和数据的上传方式,那么我们如何准备数据呢?

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

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

处理命令如下:


sort -k1,1 -k2,2n track.bedgraph > track.bedgraph.sorted  #对bedgraph进行排序bgzip  track.bedgraph.sorted  # 对bedgraph进行压缩

tabix -p bed track.bedgraph.sorted.gz  #对压缩后的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再导入washU。

甲基化数据也可以用washU来展示(如下图所示)


甲基化数据格式如下:

第一到三列代表甲基化位点的position

第四列代表甲基化的类型

第五列代表甲基化水平

第六列代表 正负链

第七列代表位点的覆盖深度

Hi-C数据可视化主要分为两种,一种是loop结构,一种是三角热图。其中三角热图数据采用.hic 文件格式。

.hic 文件主要是Juicer 产生的文件格式。



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 左端的染色体编号

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

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

第四列:loop右端的位置以及交互数值。

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

那么针对以上文件总结而言:

bgzip track.bedgraph.sorted 

tabix -p bed track.bedgraph.sorted.gz


针对bedgraph文件都要用bgzip压缩并建立索引。

那么最后举个例子,在导入longrange,以及bam或者bigwig,我们基本就有了下面这张类似与文献里的图。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,294评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,780评论 3 391
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,001评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,593评论 1 289
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,687评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,679评论 1 294
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,667评论 3 415
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,426评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,872评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,180评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,346评论 1 345
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,019评论 5 340
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,658评论 3 323
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,268评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,495评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,275评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,207评论 2 352

推荐阅读更多精彩内容