如何查询基因和通路的相关性,确定研究思路

很多小伙伴私信喵学姐,说设计课题、设计实验的时候,找不到思路,有目标基因但是又无法确定相关的通路。

今天喵学姐就教大家怎么查找基因和通路的关系,减少试错成本。

查找现有的研究结果

对于之前研究结果的确定,首先我们打开genecards(www.genecards.org/)。genecards当中总结了这个基因参与哪些经典的通路。在genecards里面汇总了KEGG等多个通路数据库当中的信息。

这里我们就以TP53为例,输入想要检索的基因选择pathways进行快速定位可以查看与这个基因有关的通路

图片

确定基因和通路内基因的关系

01.查询通路内的基因

有时候我们想研究的目标通路不在检索列表,那我们还可以进一步探究候选基因和目标通路内其他基因之间的关系。这里就可以用到PathCards(https://pathcards.genecards.org/)对目标通路进行检索,获取参与目标通路的所有基因列表

图片

在这个数据库,我们以Cell Cycle为例,输入检索就可以看到该通路所有相关的基因

图片

图片

这样我们就获得了目标通路的所有基因了。

02.蛋白相互作用分析

在我们获得基因之后,就可以看一下他们之间的互作关系,我们可以把获得的所有基因以及目标基因,统一放到STRING(https://string-db.org/)里面,就可以查看TP53和其他基因有没有关系了。

图片

图片

分析完成后在下方下载结果文件

图片

在Excel中筛选与TP53互作的蛋白

图片

03.共表达分析

如果还想探究候选基因和目标通路中的基因在某个疾病中的共表达关系,我们可以借助cBioPortal数据库(www.cbioportal.org/)。选择想要分析的疾病队列(这里以以甲状腺癌的TCGA队列为例),并输入候选基因后,点击Co-expression模块获取基因的共表达关系。

图片

  

图片

△下载数据表格流程

我们把相关的结果全部下载下来之后,在excel中筛选一下,根据需求设定相关系数的阈值。在本示例中,我们以p<0.05和r>0.5为阈值,

图片

这里可以看出只有FGD1的相关系数是大于0.5的。那如果我们后续要做一些富集分析,可以把相关系数阈值灵活调整,如设置为0.4,筛出8000多个基因,我们可以选取前100个基因进行后续的富集分析

图片

通过这些方法,大家会发现要了解基因和通路的关系也并不复杂,除了看文献,花点时间在数据库查询,也能获得不少有价值的数据,更为后续的分析、实验设计做了基础

  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
基因通路评价的打分函数通常基于基因表达数据和通路数据库,可以用来评估基因通路的激活程度和相关性,以及其与疾病和药物的关联程度。以下是一个简单的基因通路评价打分函数的示例,其中使用了基因表达数据和KEGG通路数据库: ```python import numpy as np import pandas as pd from bioservices import KEGG # 定义基因通路评价打分函数 def pathway_score(pathway_id, gene_expr, kegg_db): # 从KEGG数据库中获取基因通路信息 genes, _, _, _, _, _ = kegg_db.get(pathway_id) genes = [g.split(':')[1] for g in genes] # 从基因表达数据中筛选出对应的基因 gene_expr = gene_expr.loc[genes] # 计算基因表达数据的平均值 avg_expr = gene_expr.mean() # 计算基因表达数据的标准差 std_expr = gene_expr.std() # 计算基因表达数据的标准分 z_score = (gene_expr - avg_expr) / std_expr # 计算基因通路的得分 pathway_score = np.mean(z_score) return pathway_score ``` 该函数的输入参数包括基因通路ID、基因表达数据和KEGG通路数据库。函数首先从KEGG数据库中获取指定基因通路基因列表,然后从基因表达数据中筛选出对应的基因,并计算基因表达数据的平均值和标准差。接下来,函数计算基因表达数据的标准分,并计算基因通路的得分,最后将得分作为函数的输出返回。 在实际应用中,可以使用该函数对多个基因通路进行评价,并比较它们之间的得分来确定与疾病或药物相关联的基因通路

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值