小男孩‘自慰网亚洲一区二区,亚洲一级在线播放毛片,亚洲中文字幕av每天更新,黄aⅴ永久免费无码,91成人午夜在线精品,色网站免费在线观看,亚洲欧洲wwwww在线观看

分享

UCSC瀏覽器的可視化

 yjt2004us 2018-06-22


引言

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


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


地址:http://genome./

首先點擊My Data中的My Sessions進入頁面。


新建一個項目。

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

也可以設(shè)置數(shù)據(jù)類型的物種,基因組和參考基因組的信息,參考基因組信息如果需要轉(zhuǎn)換,可以通過LiftOver來實現(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等文件格式的上傳,主要分為通過網(wǎng)頁直接提交和服務(wù)器的URL鏈接上傳。

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

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

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

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

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


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

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

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

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)勢在于可以體現(xiàn)出數(shù)據(jù)大小的變化和高低,例如組蛋白修飾的峰值等。


1 圖像的設(shè)置


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

圖中的每個track可以點擊左側(cè)灰色條框拖拽進行上下移動。

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

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


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


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


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


圖像的保存可以點擊上方的PDF/PS進行保存。


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


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



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

    0條評論

    發(fā)表

    請遵守用戶 評論公約

    類似文章 更多