玩转基因组浏览器之展示RNA_seq中的基因表达量

欢迎关注”生信修炼手册”!

和CNV类似,IGV也可以以热图的形式展示基因表达量的数据,要求表达量文件的格式为gct, 示意如下

玩转基因组浏览器之展示RNA_seq中的基因表达量

该文件是\t分隔的纯文本文件,第一行内容总是为#1.2, 表示版本,第二行表示表达量矩阵的维度,第一个值对应探针探针个数,第二个数值代表样本个数,第三行是表达量矩阵的表头,前两列固定是NAMEDescription, NAME是探针ID,必须保证唯一,Description表示描述信息,如果没有,可以用na填充,后面每列对应一个样本。

IGV会根据Description的内容将探针映射到基因组区域,有以下两种写法

na |@RABGGTA|
na |@chr6:1950428-1950681|

第一种写法映射到基因,第二种写法映射到染色体区域,如果有多个区域的话,可以用逗号链接。GCT文件导入IGV之后,示意如下

玩转基因组浏览器之展示RNA_seq中的基因表达量

每个样本对应一个track, 根据表达量数值大小,颜色从蓝色过滤到红色,当然也可以通过菜单栏的View->Color Legends修改对应的图例

玩转基因组浏览器之展示RNA_seq中的基因表达量

·end·

—如果喜欢,快分享给你的朋友们吧—

原创不易,欢迎收藏,点赞,转发!生信知识浩瀚如海,在生信学习的道路上,让我们一起并肩作战!

本公众号深耕耘生信领域多年,具有丰富的数据分析经验,致力于提供真正有价值的数据分析服务,擅长个性化分析,欢迎有需要的老师和同学前来咨询。

  更多精彩

  写在最后

转发本文至朋友圈,后台私信截图即可加入生信交流群,和小伙伴一起学习交流。

扫描下方二维码,关注我们,解锁更多精彩内容!

玩转基因组浏览器之展示RNA_seq中的基因表达量

一个只分享干货的

生信公众号