【基因组学】使用IGV可视化查看基因组重测序比对结果

写在前面
IGV能够将基因组的变异情况可视化,因此广泛应用于基因组学的研究中。
下面我以windows下的IGV为例,介绍一下IGV在查看基因组变异方面的应用。

一.下载

下载链接如下:https://software.broadinstitute.org/software/igv/download
和一般的软件下载情况类似,不同的是,需要在电脑上安装JAVA,这个应该难不倒各位看官老爷,因此不再赘述。

二.用户使用手册

官网的使用手册包含以下18个方面,我主要关注的是 Viewing DataViewing Alignments。

1. Viewing Data

1.1 Default Display

IGV的一条track代表一个样本或试验,每条track,IGV都会展示track identifier,一个或多个属性,和track数据,如图所示。


image.png

当导入文件时,IGV会根据文件的后缀名自动识别文件。如下图所示。


image.png

1.2 Changing the Display

1.3 Segmented Data

1.4 GWAS Data

IGV可以像曼哈顿图那样展示GWAS数据。

1.5 RNA Secondary Structure

2. Viewing Alignments

看的最多的比对文件就是下机数据比对到基因组的sam/bam文件。
导入IGV的bam要附带着它的索引文件,这个可以由samtools生成。
Tracks
当你导入bam文件时,IGV会自动搜索该路径下的索引文件。
导入bam文件后,会生成3个相关的tracks。

  • Alignment Track(看单条比对reads)
  • Coverage Track(看测序覆盖的深度)
  • Splice Junction Track(提供跨越剪接连接的读取的另一种视角)
    默认情况下展示Alignment TrackCoverage Track。

Coverage Track
如果一个核苷酸与参考序列的差异大于20%的加权质量reads, IGV会根据每个碱基的reads数的比例赋予染色。

image.png

  • 通过右键单击跟踪并选择设置等位基因频率阈值来覆盖单个覆盖track的默认阈值。
  • 要更改所有覆盖轨迹的默认值,请在视图> Preferences下的Alignments选项卡中设置覆盖等位基因-分数阈值。
    将鼠标悬停在覆盖栏上查看计数细节。从右键菜单中复制计数细节到计算机的剪贴板。


    image.png

Alignment Track
默认情况下,匹配参考基因组的碱基展示为灰色。碱基的插入用紫色表示,碱基的删除用黑色横线表示。

  • Transparency for Mapped Reads
    注意,如屏幕截图所示,以浅灰色边框和透明或白色填充显示的对齐,其映射质量为零。
    在这种情况下,reads也映射到另一个同样好的位置。
    也有可能reads不能被唯一地放置,但其他的放置不一定能带来同样高质量的得分。


    image.png
  • Insertions
    IGV用紫色或红色表示碱基的插入,鼠标悬浮能够看到该碱基的详细信息。


    image.png
  • Deletions
    在有gap的reads中,IGV用黑色短线表示相对于参考基因组的删除。

    image.png

    Paired-End Alignments

  • View As Pairs
    默认情况下,IGV显示单独的reads,因为它们紧密地打包在一起。
    从右键菜单中选择View as pairs,以如下图所示的方式显示连接两端的成对。鼠标悬停(2)也会显示在普通视图(左)的单个reads或成对读取视图(右)的成对reads。


    image.png

2.1 Interpreting Color by Insert Size

IGV使用不同的颜色来标志异常的插入。当你在弹出菜单中通过插入大小选择颜色对齐>时,默认的着色方案是:

image.png

  • Deletions
    单端的缺失。


    image.png

    双端的缺失。


    image.png

    在IGV中展示的结果。
    image.png
  • Insertions
    单端示意图。


    image.png

双端示意图。


image.png

检测到异常插入大小受片段大小的限制。读长必须足够长,以跨越插入位置,并包含两端映射到参考基因组的序列,能检测到的最大插入大小为:fragment length - (2x read length)。
IGV示意图。


image.png
  • Inter-chromosomal Rearrangement
    下面这个例子中,reads一端在染色体1上,另一端在染色体6上。


    image.png

2.2 Interpreting Color by Pair Orientation

成对reads的方向可用于检测结构变异事件,包括染色体倒置,复制和易位。
绿色、蓝绿色和深蓝色的阴影表示倒置、复制和易位的结构事件;不同的颜色取决于不同的平台。


image.png

Inversions
染色体上的倒置。

image.png

双末端测序揭示的倒置。
image.png

在IGV上展示如下。
image.png

Inverted Duplication
当一大段DNA被复制并以与原始序列相反的构型插入基因组时,这被称为反向复制。
image.png

这将会有重叠的左边和右边的reads,并且有可能由覆盖深度/拷贝数所决定。
image.png

在IGV中展示如下。
image.png

Tandem Duplication

当一大段DNA被复制并插入到基因组中原始序列的旁边时,这被称为串联复制。

image.png

IGV展示如下。
image.png

Translocation on the Same Chromosome
当一大段DNA从一个位置移除并插入到其他位置时,这就是易位。
image.png

同一染色体上的易位可通过对取向的颜色编码检测,而两个染色体之间的易位可通过插入大小的染色检测。
image.png

2.3 Interpreting Color by Bisulfite Mode

2.4 Splice Junctions

2.5 Sashimi Plot

未完待续~

参考链接:
http://software.broadinstitute.org/software/igv/UserGuide

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

推荐阅读更多精彩内容

  • IGV(Integrative Genomics Viewer)是一款本地使用的基因组可视化神器,不管你用何种系统...
    EwanH阅读 4,036评论 0 14
  • 刘小泽写于19.6.8之前最开始学习IGV的时候写过两篇推文,当时简单做了了解,但有朋友说意犹未尽,嗯我也这么觉得...
    刘小泽阅读 26,027评论 6 67
  • 一. 简介 在基因组数据分析过程中,我们可以将fastq,fasta,ONT等不同格式的数据mapping到参考基...
    小飞虎阅读 6,316评论 0 19
  • 我是黑夜里大雨纷飞的人啊 1 “又到一年六月,有人笑有人哭,有人欢乐有人忧愁,有人惊喜有人失落,有的觉得收获满满有...
    陌忘宇阅读 8,532评论 28 53
  • 信任包括信任自己和信任他人 很多时候,很多事情,失败、遗憾、错过,源于不自信,不信任他人 觉得自己做不成,别人做不...
    吴氵晃阅读 6,187评论 4 8