如何查询一个基因和一个通路的相关性!
怎么查看一个基因和某一个通路之间的相关性?我以ACE2和Cell Cycle信号通路为例,两步:
一、肯定是要查询两者在之前的研究当中是否有关系。
对于之前研究结果的确定的话,我们可以通过genecards来查询。genecards当中总结了这个基因应该参与哪些经典的通路。在genecards里面汇总了KEGG等多个通路数据库当中的信息。所以在这里,我们基本上能清楚这个基因主要是参与哪些基因。我们可以先进入生物医学之家(swyxzj.com,这个是生物医学必备网站),然后选择genecards进入:
进入网站之后,输入基因,随后点击pathway:
经过查询,我们发现ACE2和Cell Cycle信号通路是没有联系的。经过上面的查询,我们发现这个基因并不参与细胞周期信号通路的调控。那这个基因不参与这个通路并不能代表这个基因和通路内的基因没有关系。因此,我们查看这个基因和通路内的基因是否存在关系。
那么首先第一步,就是需要知道通路内都有哪些基因。我们可以通过KEGG等通路数据库来查询某一个通路内的基因都有哪些。这里我推荐另外一个通路查询数据库: PathCards(https://pathcards.genecards.org/)。这个数据库和我们刚刚提到的genecards是一个机构的。这个数据库汇总了KEGG等多个数据库的一个综合性通路查询数据库。一样,生物医学之家进入这个网站:
在这个数据库,我们输入Cell Cycle即可获得符合检索条件的通路:
进一步的我们可以看到这个通路内都有哪些基因:
就这样,我们就获得了目标通路的所有基因了。在我们获得基因之后,最基本的相互作用分析,能想到的就是我们之前介绍的蛋白相互所有分析(PPI),我们可以把获得的基因内的所有基因以及目标基因(ACE2),统一放到STRING里面,就可以查看ACE2和其他基因有没有关系了。一样,从生物医学之家进入:
STRING的使用以及网络分析:
然后一路默认即可:
经过分析,我们发现,ACE2只是在文本挖掘的结果当中与CDK4基因可能存在相互作用关系。这个结果照这个情况来看,还不是特别确定。因此在目前的蛋白相互相互作用的层面,有可能这这个基因和细胞周期相关基因是没有关系的。
二、预测两者是否相关。
我们做的是在蛋白层面观察基因之前是否存在相互作用关系。但是在mRNA层面的话,还没有查看具体的相关性。这个时候,我们用到的就是类似芯片或者RNA-seq方面的东西了。之前我们在单基因开展研究的帖子当中,提到了过最好是能在GEO里面找到有人做过的相关敲除或者过表达这个基因的芯片。这样的话,我们查看一下差异基因就可以了。
经过查询,我们发现并没有ACE2相关的过表达/敲除芯片。
既然没有的话,那我们就只能使用疾病相关的数据来提取目标基因的表达,来进行相关分析,进而来查看这些基因的是否存在共表达关系了。如果是肿瘤研究的话,那直接就使用TCGA的数据就可以来进行研究。如果是其他疾病的话,那可能就需要去GEO里面搜索自己想有的数据集了。对于TCGA的数据,也是有很好的数据库来直接得到分析结果的。例如,我们这里想要看在肠癌当中ACE2和这些基因的关系,那么我们就可以在cbioportal数据库当中,寻找共表达基因。这样,我们就可以获得和ACE2存在共表达关系的基因了。
由于细胞周期相关的基因有100+,我们不能一个一个来检索相关结果。所以我们就把相关的结果全部下载下来之后,在excel当中进行一下vlookup一下。最后我们发现,ACE2和16个细胞周期相关基因存在共表达关系。
所以,这里涉及到了GEO数据库和TCGA数据库的挖掘,现在绝大部分纯生信发低分SCI的数据都来源于这里。
我们今天就将原价1999元的:
“手把手教你实操挖掘GEO数据库发SCI课程”
免费发送给大家。如下即可获取:
🔽长按下方二维码关注🔽
对话框输入关键词:GEO
资源截图如下:
1、从0做到8的200+纯生信文献套路----范文:生信分析,最关键的就是有没有思路了,所以范文很重要啊。
2、手把手教你按照范文发纯生信SCI论文视频教程:
赶快免费领取吧,先下载收藏起来,以后学生信的时候学起来就方便啦。
“手把手教你按照范文发纯生信SCI论文视频教程”
免费领取方式:🔽①长按下方二维码关注🔽
②对话框输入关键词:GEO
更多资源将继续更新,请扫码关注:
辛苦整理,全文未插入任何广告!
觉得有用的话,您就点个在看、点赞!