挑战7天光速完成一篇药物靶向MR,Day 4-5!
进度汇报:寻找靶基因+跑代码+出图片
药物靶向的代码非常的简单
也非常的好理解,但是要注意的是要搞清楚整个过程中的思路
第1步:明确蛋白与疾病表型存在明显的因果关系或者已知该蛋白能够有效作用于该疾病的治疗。
第2步:寻找调节该蛋白表达的一类基因来模拟药物作用
第3步:寻找到这个靶基因后获取其工具变量与疾病进行单变量MR分析。
这就是我们的药靶的文章设计思路
很重要的一部就是寻找靶基因
一般我们通过eQTL的网站无寻找eqtlgen
因为基因表达的过程会受到很多的调控
所以不仅要关注该基因,还要关注该基因的上下游1兆的基因信息
我的eGTL的数据来源于Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders,里面包括了来自人体血浆、CSF等的eQTL数据
我对这些下载来的数据进行了预处理
实际上本质就是本地数据的数据处理
然后再用这些数据进行药靶MR的分析
然后再做一个共定位分析
其实就是看暴露X与结局Y的GWAS summary中该区域的所有SNP的位置是不是基本一致的
今日图表弄完啦,就要开始写作啦!