日韩黑丝制服一区视频播放|日韩欧美人妻丝袜视频在线观看|九九影院一级蜜桃|亚洲中文在线导航|青草草视频在线观看|婷婷五月色伊人网站|日本一区二区在线|国产AV一二三四区毛片|正在播放久草视频|亚洲色图精品一区

分享

UCSC瀏覽器的可視化

 yjt2004us 2018-06-22


引言

在生信的研究過(guò)程中,數(shù)據(jù)的展示往往是非常重要的一步,而選擇以哪種方式展示更是令人頭疼。在這里,我們介紹UCSC瀏覽器中實(shí)現(xiàn)可視化的過(guò)程,可以在參考基因組上展示出數(shù)據(jù)的信息,如感興趣的基因附近的甲基化值、表達(dá)值或差異的區(qū)域等,當(dāng)然更重要的是方便自己對(duì)數(shù)據(jù)的了解。


1 創(chuàng)建一個(gè)項(xiàng)目


地址:http://genome./

首先點(diǎn)擊My Data中的My Sessions進(jìn)入頁(yè)面。


新建一個(gè)項(xiàng)目。

點(diǎn)擊新建的項(xiàng)目進(jìn)入可視化的頁(yè)面,在圖片的下方點(diǎn)擊add custom tracks就是上傳自己的數(shù)據(jù)了。

也可以設(shè)置數(shù)據(jù)類型的物種,基因組和參考基因組的信息,參考基因組信息如果需要轉(zhuǎn)換,可以通過(guò)LiftOver來(lái)實(shí)現(xiàn)(hg19與hg38的相互轉(zhuǎn)換)。


1 數(shù)據(jù)的上傳


瀏覽器支持bigBed, bigChain, bigGenePred, bigMaf, bigPsl, bigWig, barChart,bigBarChart, BAM, VCF, BED, BED detail, bedGraph, broadPeak, CRAM, GFF, GTF,MAF, narrowPeak, Personal Genome SNP, PSL, 或 WIG等文件格式的上傳,主要分為通過(guò)網(wǎng)頁(yè)直接提交和服務(wù)器的URL鏈接上傳。

由于數(shù)據(jù)類型較多,通過(guò)可以直接上傳的BED文件為例介紹:

BED文件由3列必需的信息和9個(gè)額外可選擇的信息組成,三列必須信息就是:染色體,起始位置,終止位置。其他9個(gè)額外的信息如下:

BED文件的第一行是對(duì)這個(gè)條目的設(shè)置,包括名字、描述和顏色等,Color的選擇按照RGB顏色的參數(shù)可以選擇不同的顏色。上傳BED文件的數(shù)據(jù)是不能體現(xiàn)出數(shù)據(jù)的變化、高低等現(xiàn)象,只能顯示數(shù)據(jù)的區(qū)域信息,因此需要更高級(jí)的設(shè)置,即設(shè)置區(qū)域的顏色來(lái)區(qū)別數(shù)據(jù)的大小等。這個(gè)功能就是BED文件的第9列itemRgb選項(xiàng),需要注意的是,如果要使用第9列的信息則前面的幾列信息也要存在不能空缺。

保存后點(diǎn)擊go進(jìn)行在基因組上的可視化,可以尋找感興趣區(qū)域的數(shù)據(jù),根據(jù)顏色的深淺判斷甲基化值的大小。

而最簡(jiǎn)單的格式只需要前三列就能上傳查看,并且其他格式的數(shù)據(jù)也能很方面的轉(zhuǎn)換成BED格式進(jìn)行上傳。


當(dāng)文件大小超過(guò)50M時(shí),則不推薦通過(guò)網(wǎng)頁(yè)上傳數(shù)據(jù),這時(shí)候需要一個(gè)URL鏈接,通常是bigBed, bigWig, bigGenePred, BAM 和VCF格式的數(shù)據(jù),我們以bigWig格式為例。

bigWig是通過(guò)wig格式的文件轉(zhuǎn)換的二進(jìn)制壓縮文件,而wig文件也是和BED文件類似的包含區(qū)域信息的文件,一般使用MACS峰值探測(cè)后可以產(chǎn)生wig格式的文件。wig格式是可以直接上傳的但超過(guò)50m就推薦轉(zhuǎn)換成bigWig文件上傳。轉(zhuǎn)換命令如下:

wigToBigWig input.wig chrom.sizes myBigWig.bw(若出現(xiàn)超出染色體長(zhǎng)度的錯(cuò)誤,需編輯染色體信息使其長(zhǎng)度增加)

wigToBigWig程序下載地址:http://hgdownload.soe./admin/exe/linux.x86_64/wigToBigWig

chrom.sizes染色體信息下載地址:http://hgdownload.soe./admin/exe/

上傳的例子如下,需要在bigDataUrl后增加數(shù)據(jù)的在服務(wù)器中的鏈接。

wig和bigWig文件的優(yōu)勢(shì)在于可以體現(xiàn)出數(shù)據(jù)大小的變化和高低,例如組蛋白修飾的峰值等。


1 圖像的設(shè)置


首先是最上方的設(shè)置,對(duì)圖中的區(qū)域進(jìn)行放大和縮小(基于堿基數(shù)目),下方的輸入框可以輸入具體基因組的位置(tab分隔)或一個(gè)條目的名字跳轉(zhuǎn)到該處,也可以點(diǎn)擊圖像左右拖拽進(jìn)行微調(diào)。

圖中的每個(gè)track可以點(diǎn)擊左側(cè)灰色條框拖拽進(jìn)行上下移動(dòng)。

對(duì)于每一個(gè)上傳后的track在圖中都可以右鍵灰色條框設(shè)置其顯示效果。如果上傳的文件對(duì)每段區(qū)域有定義的名字,則會(huì)多出pack和squish的選項(xiàng)。

例如以瀏覽器中自帶的CpG島為例,選擇了pack后則會(huì)在每個(gè)CpG島前顯示其名字,這是比較合適顯示效果。用戶可以根據(jù)自己上傳的數(shù)據(jù)選擇不同的設(shè)置進(jìn)行調(diào)整。


對(duì)于有峰值顯示的數(shù)據(jù),可以左鍵灰色條框設(shè)置其顯示的高度。


圖像下方的一行設(shè)置中tracksearch可以搜索感興趣的條目,manage custom tracks即是上傳文件,configure是設(shè)置圖像的寬度,resize是調(diào)整圖像的寬度和瀏覽器寬度一致。


關(guān)于最下方的設(shè)置,是瀏覽器自帶的基因組中的參考信息,可以對(duì)沒(méi)有幫助的進(jìn)行隱藏,如果沒(méi)有顯示出感興趣的信息,可以使用track search進(jìn)行搜索添加。


圖像的保存可以點(diǎn)擊上方的PDF/PS進(jìn)行保存。


最后則是這個(gè)項(xiàng)目的保存,網(wǎng)頁(yè)中沒(méi)有保存的選項(xiàng),因此需要重復(fù)項(xiàng)目創(chuàng)建的過(guò)程,在界面上方點(diǎn)擊My Data中的My Sessions,使用當(dāng)前項(xiàng)目的名字進(jìn)行保存覆蓋即可。


小結(jié):本文僅提供了最簡(jiǎn)單直接的可視化實(shí)現(xiàn)方法,更詳細(xì)的介紹在UCSC官網(wǎng):http://genome./FAQ/FAQformat.html#format1。Genome Browser的可視化更多的是方面自己對(duì)數(shù)據(jù)的理解,更直觀的查看識(shí)別出的差異區(qū)域是否顯著,或是感興趣的例如CpG島、啟動(dòng)子等區(qū)域的信息,有助于更進(jìn)一步的分析。



    本站是提供個(gè)人知識(shí)管理的網(wǎng)絡(luò)存儲(chǔ)空間,所有內(nèi)容均由用戶發(fā)布,不代表本站觀點(diǎn)。請(qǐng)注意甄別內(nèi)容中的聯(lián)系方式、誘導(dǎo)購(gòu)買等信息,謹(jǐn)防詐騙。如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,請(qǐng)點(diǎn)擊一鍵舉報(bào)。
    轉(zhuǎn)藏 分享 獻(xiàn)花(0

    0條評(píng)論

    發(fā)表

    請(qǐng)遵守用戶 評(píng)論公約

    類似文章 更多