eQTL hotspot : eQTL热点区域

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

众所周知，eQTL分为了cis-eQTL和trans-eQTL两种作用方式，cis模式下只能调控临近的基因，而trans模式突破了距离限制，在该模式下一个eQTL位点潜在的靶标基因数量大大增加。在eQTL-gene构成的调控网络中，eQTL节点的degree并不是均匀分布，往往少数几个位点的degree很高，表示这些位点调控了大多数的基因，我们将这些调控了多个基因的eQTL位点称之为eQTL hotspot。

关于eQTL hotspot的分析和可视化，可以参考以下文章

https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-4871-y

本文通过eQTL分析来识别与肌内脂肪含量， 即IMF这一性状相关的调控区域。首先通过matrixEQTL软件进行eQTL分析，采用默认的距离阈值，SNP与gene的距离小于1Mb以内的为cis-eQTL, 大于1Mb的为trans-eQTL，识别到了1268个cis-eQTL, 10334个trans-eQTL, 通过之前介绍的cis-trans plot展示所有eQTL在染色体上的分布

上图中横坐标和纵坐标分别表示SNP和gene在染色体上的位置，图中的散点表示eQTL-gene对，其中位于对角线附近的点为cis-eQTL, 其他区域的点为trans-eQTL。

以每个eQTL位点调控的基因数目为指标，来识别eQTL hotspot，如下图所示

水平的蓝色虚线表示基因个数的阈值，文章中采用的是119，横坐标为eQTL位点的染色体位置，纵坐标为每个eQTL位点对应的靶基因个数，箭头标记的地方就是预测出来的eQTL hotspot。

初始的eQTL hotspot就是几个独立的SNP位点，考虑到SNP位点之间的连锁，分析这些位点之间的LD程度，将连锁程度高的位点合并，归于同一个区域，这些区域就是最终识别到的eQTL hotspot region。

对于hotspot region, 需要进行后续的功能注释，进一步分析其功能。可以对这些区域覆盖的基因进行GO/KEGG等功能富集分析，也可以注释这些区域中包含的转录因子，通过转录因子来联通eQTL和其靶基因。本文就是分析了eQTL hotspot对应的转录因子，找到了3个转录因子，可视化如下

采用circos图来绘制eQTL和其靶基因的调控关系。对于eQTL分析而言，识别到eQTL位点只是第一步，通过eQTL hotspot区域的功能注释来研究eQTL的功能作用机制，才是真正具有生物学意义的数据分析。

·end·

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

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

  更多精彩

• KEGG数据库，除了pathway你还知道哪些

• DNA甲基化数据分析专题

• 突变检测数据分析专题

• mRNA数据分析专题

• lncRNA数据分析专题

• circRNA数据分析专题

• miRNA数据分析专题

• 单细胞转录组数据分析专题

• chip_seq数据分析专题

• Hi-C数据分析专题

• HLA数据分析专题

• TCGA肿瘤数据分析专题

• 基因组组装数据分析专题

• CNV数据分析专题

• GWAS数据分析专题

  写在最后

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

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

生物信息入门

只差这一个

公众号

点击下方“阅读原文”，我们一起进步！