UCSC-browser学习:创建自己的track hubs

(由于以前的几篇文章被平台锁了且申诉失败,只好重发一下)

在看paper的时候,经常会有这种图:


这是从UCSC browser上可视化数据然后截图下来的,我们可以将数据传到UCSC上,制作自己的track hubs,然后选择我们想看的区域截图就可以了。

1. Track Hubs的介绍

简单来说,track hubs相当于在UCSC browser上的一个可网络共享的一个基因组可视化工具。


在UCSC的页面选择track hubs,即可进入页面:


在国内一般会提示:


建议选择亚洲区,速度会快很多。

进去后在public hubs中可以看到别人做的hubs:


点connect,就可以看到别人做的hubs了:


在Hub development里可以检查我们创建的hub.txt是否有问题和直接可视化(此处xxx隐藏了部分网址),如果hub.txt中有错误的参数设置、找不到文件等这里会报错:


在My hubs中,贴上自己的hub.txt的url,就可以add进来了:


点击Assemblies下面的mm10即可进入,可以修改position看不同区域的情况:


2. 设置Track Hubs

核心的问题来了,上面的URL怎么得到呢?

首先我们需要创建类似于这样的文件夹:
Myhub/
   -- genomes.txt
   -- hub.txt
    mm10/
    -- trackDb.txt
    --Mydata.bw

在Myhub文件夹下,有genome.txt文件,hub.txt文件和mm10文件夹,mm10下有trackDb.txt文件和我们的基因组文件(bigwig, bed等)。

UCSC上的例子:


genome.txt必须的内容(基因组版本和trackDb.txt的路径):

genome mm10
trackDb mm10/trackDb.txt

hub.txt的内容(hub的名字;标签;genomes.txt;email;description):

hub ENCODE_mm10_kallisto
shortLabel kallisto
longLabel kallisto mm10 data
genomesFile genomes.txt
email xxx@xxx.com
descriptionUrl description.txt

(description可以不加,不过在Hub development中check的时候会出现一个warning)

mm10/trackHub.txt的内容(仅展示了一部分):

track neural-tube_15.5_1_forward
bigDataUrl neural-tube_15.5_1_forward.bw
shortLabel neural-tube 15.5 rep1 forward
longLabel neural-tube 15.5 rep1 forward
type bigWig
visibility full
autoScale on
maxHeightPixels 128:32:11
priority 1
color 255,0,0

其中:
bigDataUrl是我们要可视化的文件的名字;
type是格式,不仅仅支持bigwig还有bigBed、VCF、BAM等等格式。具体可以参考https://genome.ucsc.edu/goldenpath/help/hgTrackHubHelp.html中的介绍;
priority是优先级;
color不加默认是黑色,仅支持RGB格式。

然后我们将hub.txt的链接拷贝到UCSC的页面中就可以了。

另外,我们也可以通过链接访问到,比如如下链接:
http://genome.ucsc.edu/cgi-bin/hgTracks?db=hg19&hubUrl=http://ftp.ebi.ac.uk/pub/databases/ensembl/encode/integration_data_jan2011/hub.txt

db=hg19这里必须改成自己的基因组版本,后面的hubUrl也需要改成自己的。并且,最重要的是,如果我们使用亚洲版本的UCSC创建的hub的话,这个链接前面必须是http://genome-asia.ucsc.edu/,否则显示不出来。

3. 保存自己的hubs

在可视化界面,选择My sessions:


随后会跳转到如下界面,有账号的登录,没账号的需要注册一个:


登录了之后在Save Settings选择submit即可。


欢迎关注~

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