2024-01-16 SNP位点与转录因子结合分析记录1

本人既往研究的关键SNP位点信息:rs2243890 - SNP - NCBI (nih.gov)

Decreased activity of RCAN1.4 is a potential risk factor for congenital heart disease in a Han Chinese population (fudan.edu.cn)
前期完成该位点部分研究,实验证实启动子区的SNP位点影响转录因子结合,推测结合了转录抑制因子,因此需要进一步研究具体结合的转录因子。

全基因组关联分析(GWAS)发现的很多变异位点基本为非编码,这些变异位点:
1)要么调控基因表达(eQTL);
2)要么影响增强子活性;
3)要么影响转录因子(TF)结合特异性;
4)要么啥也不是。

针对以上四种情况:

1)是否调控基因表达(eQTL)可通过GTEx Portal查询。

无法打开链接

2)是否影响增强子活性查找相关资料找到之前的推文:感兴趣的SNP/区域上是否有增强子/转录因子?增强子/转录因子调控哪个靶基因?(EnhancerDB) (qq.com)

哇塞,好人,点赞,跳转链接,激动...

哦豁,挂了

3)是否影响转录因子(TF)结合特异性则可通过如下数据库进行查询:GVATdb (sdsc.edu)

Natures有一篇文献Systematic analysis of binding of transcription factors to noncoding variants。

该文献针对95,886个常见变异位点(SNPs,欧洲和亚洲人群的MAF> 1%)与270个转录因子的结合特异性进行了大量的SNP-SELEX实验,并以此构建了GVATdb数据库。

下图是针对转录因子或者SNP位点进行检索的实验结果图:


网站示例SNP

输入目标SNP:rs2243890


唉,也不知道是没有收录还是真的没有结合,忧伤

还是学习一下吧,每列的释义如下:

  • oligo_auc:转录因子与40bp DNA 序列的结合得分, 用AUC(Area under Curve)值表示;
  • oligo_pval:对OBS进行25,000次蒙特卡洛随机化后得到的p值。p<0.05表示TF与基因组片段的特异性“结合”;
  • Alt:SNP(hg19)的替代等位基因;
  • Ref:SNP(hg19)的参考等位基因;
  • ref_auc:ref 与 TF 的结合得分;
  • alt_auc:alt 与 TF 的结合得分;
  • pbs:结合倾向性得分,公式为:Ref 等位基因得分减去 Alt 等位基因得分,负值表示转录因子更倾向于结合 Alt 等位基因;
  • p-value:对 PBS 进行25,000次蒙特卡洛随机化后得到的 p 值。 p<0.01 表示 TF 与 Ref 等位基因和 Alt 等位基因结合“存在差异”。
    此外,对于没有纳入 SNP-SELEX 实验的 SNP 位点,作者还建立了 deltaSVM 模型,用于预测未纳入的 SNP 位点与 TF 的结合特异性,如下图所示:
提交了位置信息还是没有收到结果

又搜索了一下,找到另一个网站进行预测:两步预测SNP与转录因子结合的神器分享 (biodog.cn)
按照提示输入目标SNP

没有结果-_-

先记录到这里,后面再继续学习吧。

参考文献:Yan J, Qiu Y, Dos Santos A M R, et al. Systematic analysis of binding of transcription factors to noncoding variants[J]. Nature, 2021: 1-5.

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容