蛋白互作网络分析同样也是 SCI 文章中出现频率最高的分析方法,下面就给大家梳理一下流程!

另外大家都是私信我,现在有需要学习生信等问题可以扫码进群交流,如下:

前言

蛋白质互作网络(PPI, Protein-Protein Interaction Networks)是由蛋白通过彼此之间的相互作用构成,来参与生物信号传递、基因表达调节、能量和物质代谢及细胞周期调控等生命过程的各个环节。系统分析大量蛋白在生物系统中的相互作用关系,对了解生物系统中蛋白质的工作原理,了解疾病等特殊生理状态下生物信号和能量物质代谢的反应机制,以及了解蛋白之间的功能联系都有重要意义。

蛋白互作分析的工具包括 string 数据库也是一个在线工具,另外就是经典的网络作图软件 Cytoscape。

实例解析

我们就根据 RNA 12. SCI 文章中肿瘤免疫浸润计算方法之 CIBERSORT 结果来学习怎么绘制蛋白互作网络图。

1 数据读取

读取表达数据,由于我们挑选的都是免疫细胞,所以更加方便作图,选取其中的10个基因,如下:

exp<-read.table("exp.txt",header=T,check.names = F)
exp[1:5,1:3]
#SYMBOL TCGA-3L-AA1B-01A-11R-A37K-07 TCGA-4N-A93T-01A-11R-A37K-07
#1 DBNDD1                  145.7848372                     932.6654
#2  HSPB6                  463.5507914                     125.3295
#3   PDK4                  199.3329095                     223.2314
#4  ABCB5                    0.7485274                       0.0000
#5 PRSS22                  441.2341966                    1428.1244exp$SYMBOL[1:10]
#[1] "DBNDD1"    "HSPB6"     "PDK4"      "ABCB5"     "PRSS22"    "ASB4"      "TMEM132A"
#[8] "TNFRSF12A" "DLX6"      "CCL26"

2 STRING 数据库

string (vision 11.5) 数据库是一个搜索已知蛋白质之间和预测蛋白质之间相互作用的数据库,该数据库可应用于2031个物种,包含960万种蛋白和1380万种蛋白质之间的相互作用。

STRING https://cn.string-db.org


在官网的 help 里面有详细的说明,如下:

help https://cn.string-db.org/cgi/help?sessionId=bfsfWQPO6oix

Starting Point

输入文件格式必须在 STRING 起始页使用输入表单,左侧包括输入的格式,可以是蛋白名称,或者是序列,当然可以是单个也可以是多个蛋白,我们输入1个 m6A基因 FTO, 之后选择 Select, 之后在点击 SEARCH,如下:

这样我们就获得了与 FTO 相关的蛋白,如下:

如果以上图作为差异蛋白互作网络不够美观,也并不是所有差异蛋白均可形成互作网络,因此需将分析结果导出使用Cytoscape绘图。

点击Exports,选择…as simple tabular text output,点击download:

Excel表格打开保存好的结果,其中包括节点信息、node1_string_internal_id与combined_score等诸多信息:

点击某个蛋白小球会得到关于这个蛋白的信息以及蛋白结构,非常好用,自己研究吧!

接下来我们就利用 Cytoscape 做出更美观的图形。

3 Cytoscape

Cytoscape是一款免费图形化显示网络并进行分析和编辑的软件,用起来非常方便,最重要的是能够

其最新版本为 Cytoscape3.9.1,下载地址为:

Cytoscape https://cytoscape.org/

使用过程中如遇到不会,可以直接到 cytoscape 官网上查看说明,网页上有详细的使用说明:

manual https://manual.cytoscape.org/en/stable/

导入互作网络

打开软,我们看到该软件提供了一些例子,如下:

点击下图Import Networks from File System,导入Exports处下载好的数据,点击OK,导入数据完成,如下:

导入完成后得到可形成互作关系的差异蛋白网络,不能形成互作关系的单个差异蛋白不再出现:

设置节点间连接线的粗细

combined_score表示数据的支持度,combined_score越粗两种蛋白间互作关系越强,默认筛选阈值为0.4。

节点Edges的粗细由combined_score渐变调整,选择Tool里的Network Analysis工具中的Generate Style from Statistics:

这时我们可以对蛋白之间的关系执行一系列的分析,如下:

设置参数,根据需求自行设置,怎么好看怎么设置,如下:

蛋白互作这块我们就简单说下,具体一些细节性的设置软件包本身非常好用,并且带有说明,自行设置即可!

关注公众号,每日更新,扫码进群交流不停歇,马上就出视频版,关注我,您最佳的选择!

桓峰基因

生物信息分析,SCI文章撰写及生物信息基础知识学习:R语言学习,perl基础编程,linux系统命令,Python遇见更好的你

52篇原创内容

公众号

RNA 14. SCI 文章中差异表达基因之 蛋白互作网络 (PPI)相关推荐

  1. RNA 27 SCI文章中转录因子结合motif富集到调控网络 (RcisTarget)

    点击关注,桓峰基因 桓峰基因公众号推出转录组分析和临床预测模型教程,有需要生信的老师可以联系我们!首选看下转录分析教程整理如下: RNA 1. 基因表达那些事–基于 GEO RNA 2. SCI文章中 ...

  2. RNA 13. SCI 文章中加权基因共表达网络分析之 WGCNA

    WGCNA 分析流程 2008 年发表在 BMC 之后的影响力还是很高的,先后在各大期刊都能看到,但是就其分析的过程来看,还是需要有一定 R 语言的基础才能完整的复现出来文章中的结果,这期就搞出来供大 ...

  3. RNA 25. SCI文章中只有生信没有实验该怎么办?

    今天来介绍一个使用非常方便的在线免疫组化分析工具--PHA (The Human atlas),免疫组化时单基因干湿结合生信中最常见的补充实验的方法之一,性价比较高.但是如果没有条件进行自己样本的免疫 ...

  4. RNA 24. SCI文章中基于TCGA的免疫浸润细胞分析的在线小工具——TIMER

    点击关注,桓峰基因 桓峰基因 生物信息分析,SCI文章撰写及生物信息基础知识学习:R语言学习,perl基础编程,linux系统命令,Python遇见更好的你 135篇原创内容 公众号 今天来介绍一个使 ...

  5. RNA 25. SCI文章中估计组织浸润免疫细胞和基质细胞群的群体丰度(MCP-counter)

    点击关注,桓峰基因 今天来介绍一个利用基因表达估计组织浸润免疫细胞和基质细胞群的群体丰度的软件包--MCP-counter,亲试,非常好用. 桓峰基因的教程不但教您怎么使用,还会定期分析一些相关的文章 ...

  6. RNA 30. SCI文章中基于TCGA和GTEx数据挖掘神器(GEPIA2)

    这期介绍一个基于TCGA和GTEx数据挖掘神器(GEPIA2),个人觉得如果没有编程基础的可以直接利用这个在线小工具分析自己的研究的单个基因或者多个基因,效果还是蛮好的! 桓峰基因公众号推出转录组分析 ...

  7. RNA 29. SCI文章中基于TCGA的免疫浸润细胞分析 (TIMER2.0)

    桓峰基因公众号推出转录组分析教程,有需要生信的老师可以联系我们!转录分析教程整理如下: RNA 1. 基因表达那些事--基于 GEO RNA 2. SCI文章中基于GEO的差异表达基因之 limma ...

  8. 如何计算环境遗传相关 | 育种中的基因与环境互作

    基因与环境互作 1. 环境 参考: http://www.isbreeding.net/common/UploadFiles/file/teaching/%E6%95%B0%E9%87%8F%E9%8 ...

  9. mirna富集分析_GEO芯片数据下载,矩阵提取,差异基因分析,差异miRNA分析,miRNA靶基因预测,GO、KEGG功能,蛋白互作网络构建 - 生物信息学讨论版 -丁香园论坛...

    一.芯片差异基因分析 1. 芯片数据收集 在 NCBI GEO数据库下载 . 2.做差异分析 使用limma R包计算正常组织和病组织的差异表达情况 3.绘制火山图 4.绘制热图 使用pheatmap ...

  10. RNA 4. SCI 文章中基于TCGA 差异表达之 edgeR

    01. 前言 在不久的将来,预计新兴的数字基因表达(DGE) 技术在许多功能基因组学应用方面将超越微阵列技术.其中的一个基础数据分析任务,尤其是基因表达研究,涉及确定是否有证据表明对于转录本或外显子在 ...

最新文章

  1. 收缩 tempdb 数据库
  2. vue-router的基本使用
  3. 设置elf文件链接库的路径
  4. 腾达W303R v3 无线路由器使用移动光纤无法打开视频的设置方法
  5. IP网络设计系列之-局域网设计
  6. 用python快速合并代码(方便软著申请)
  7. Webpack 常用命令总结以及常用打包压缩方法
  8. android ImageView 之 android:scaleTye=
  9. GPS NMEA-0183协议常用数据格式及解析攻略
  10. 解密昇腾AI处理器--DaVinci架构(总览)
  11. GIS概念的总结(一)什么是GIS
  12. OpenCart框架运行流程介绍opencart资料链接
  13. 思科交换机命令大全 一
  14. 2022手机号段大全、归属运营商整理—2022.01.04更新(包含三大运营商)
  15. u盘里的东西删除怎么还原,分享2个指南
  16. 蛋糕究竟是怎样做大的
  17. 二叉树的镜像和对称二叉树
  18. 面试利器,精心整理了份Python数据分析,知识点高清速查表!
  19. 【商城报错】-商品上架功能报错
  20. 【CF917E】Upside Down(哈希二分)(后缀数组)(AC自动机)

热门文章

  1. 三维扫描仪在工业生产中的作用
  2. 关于谷歌浏览器不能打开页面问题的解决
  3. 律师视角下网络爬虫技术的罪与罚
  4. 五线谱音名和组别对照表_钢琴音区名称及琴键名称对照表
  5. canvas绘制网易云logo
  6. 计算机教案制作电子表格,制作电子表格教案
  7. Kettle使用 js 文件生成节假日表文件 附带2019-2020节假日文件
  8. zabbix邮箱告警配置
  9. “校园知网”5月5日冲刺计划书
  10. 谈谈公司内部的技术分享