IGV基因组浏览器可视化高通量测序数据

2018-03-30 13:24:56 浏览数 (1)

IGV是本地浏览测序数据功能最为强大的基因组浏览器,支持多种不同类型的输入格式和不同的显示方式,如峰图、线图、柱状图、Sashimi-plot。同时还可以配合bedtools使用。

新版的Windows IGV分发包中包含了Java运行环境,点击链接http://data.broadinstitute.org/igv/projects/downloads/IGV_Win_2.3.93.zip即可下载,解压后,点击igv.bat文件即可启动。

若启动失败,使用记事本打开并编辑igv.bat文件,在文件的最后新起一行输入pause,保存后,再尝试打开,就可以在Windows下的命令行界面(黑色背景的小框中)看到错误信息,根据信息提示去解决问题。

IGV具体使用如下:

如果基因组未存在于当前列表中,则选择更多从Broad IGV获取对应的基因组信息。若所研究的物种没有被Broad IGV收录,则需要自己构建基因组信息。具体参见前面发过的文章 (在公众号后台输入 IGV基因组索引 即可获取)。

IGV展示所用的bw文件可以直接从bam文件得来。下面列出几个工具和相应的代码辅助转换。考虑到不同样品的测序深度不同,需要根据样品的测序量做一个标准化处理,建议使用deepTools的bamCoverage命令。

IGV是一个界面可操作软件,学习的方式就是每个菜单点开看看,点选下看下效果,比如修改下Track的颜色,Track的高度,纵轴的坐标范围Auto scale, Log scale 或 Data range等方便比较。

如果想比较多个基因的表达峰图,可以先在IGV中添加相应的基因列表,然后针对特定的列表选择view,就可以在一个屏幕显示所有加入的基因。

如果我们要查询的区域很多,可以配合bedtools igv -i input.bed,生成一个脚本,批量捕获多个区域的比对信息或coverage信息。使用Tools - Run Batch Script - 选择bedtools igv的输出结果,即可以看到动画般的截屏界面。完成后,就可以一张张的浏览图片或者传给没有安装IGV的老师或朋友查看。

构建自己的基因组索引

  1. 在IGV的菜单栏点击 Genomes>Create .genome File,会弹出一个窗口填写需要的信息。
  2. 输入一个ID和基因组的描述信息,这个ID等同于常见的mm10,会作为构建的基因组索引的名字。
  3. 使用 samtools faidx genome.fasta索引基因组文件,然后输入基因组FASTA格式文件的路径。
  4. 提供基因注释的GTF文件,这个文件可以是拼装基因组预测的基因,也可以是cufflinks 等拼装的转录组的GTF文件。
  5. 其它的选项可以忽略。
  6. 点击 Save,选择基因组索引文件*.genome的存储路径。

0 人点赞