草庐IT

STRINGdb | 交互式可视化蛋白网络互作

研究生学生信 2023-10-09 原文

在STRING数据库中提交基因，获取PPI网络是比较常规的操作，包括下游使用cytoscape进行美化和hub网络的提取，但过程还是比较麻烦，本期介绍通过R语言来挖掘STRING数据库

1. 效果展示

通过利用差异分析结果，使用STRINGdb包挖掘蛋白交互关系，使用visNetwork构建互作网络，并利用RPubs将网络上传至服务器，实现蛋白互作网络的交互式可视化

image

2. 前期准备

2.1 R包

library(STRINGdb)
library(visNetwork)

2.2 差异分析结果


> head(diff_exp_example1)
     genes      logFC    logCPM       PValue          FDR
1   CYP1A1  -9.900908  6.373327 6.248616e-10 8.596222e-06
2     CSF3 -10.224338 10.540303 4.759447e-09 3.273786e-05
3     IL1B  -6.791297  7.881774 2.937997e-08 1.347267e-04
4 TNFSF13B   3.304202  5.305227 8.005570e-08 2.280509e-04
5     FOSB  -5.449454  7.196261 8.288541e-08 2.280509e-04
6    PTGDS   3.771284 11.508825 2.298475e-07 5.270021e-04

2.3 备注

本期内容主要借鉴了部分STRINGcb的教程：STRINGdb.pdf，里面也有很多其他功能，并未纳入进来，感兴趣的朋友可以看一下原文档~

3. 代码

将diffSig.xls放入file文件夹中，一定保证列名顺序一致

library(STRINGdb)
library(visNetwork)

string_db <- STRINGdb$new( version="11",# STRING版本
                           species=9606, # 物种
                           score_threshold=700, # 交互强度阈值，如果低于此阈值的关系将不被纳入
                           input_directory="./file") #可以设置，将数据库配置文件保存到本地
diff_exp_example1<-read.table("./file/diffSig.xls", #输入文件
                              sep = "\t",header = T)
example1_mapped <- string_db$map(diff_exp_example1,
                                  "genes",#gene列列名
                                 removeUnmappedRows = TRUE )

inter_raw<-string_db$get_interactions(example1_mapped$STRING_id)

plot_data<-merge(inter_raw,example1_mapped,by.x = "from",by.y = "STRING_id")

edges <- plot_data[,c(1:2)]

nodes <- example1_mapped
nodes$type <- ifelse(nodes$logFC>0,"up","down")
nodes$label <- nodes$gene
colnames(nodes)[c(1,4,2,6,7,8)]<-c("gene", #genesymbol列
                                   "pvalue", # P值列
                                   "logFC",
                                   "id", # ENSP列
                                   "group", # 上下调列
                                   "label") # genesymbol列
nodes <- nodes[nodes$id%in%c(edges$from,edges$to),]

visNetwork(nodes, edges) %>%
  visIgraphLayout()%>%
  # darkblue square with shadow for group "A"
  visGroups(groupname = "down", color = "#034889", shape = "square", 
            shadow = list(enabled = TRUE)) %>% 
  # red triangle for group "B"
  visGroups(groupname = "up", color = "#af002b", shape = "triangle")%>%
  visLegend()%>%
  visNodes(font=list(size = 50))%>%
  visEdges(
    shadow = FALSE,
    color = list(color = "#0085AF", highlight = "#C62F4B")
  ) %>%
  visOptions(highlightNearest = list(enabled = T, degree = 1, hover = T),
             selectedBy = "group") %>% 
  visLayout(randomSeed = 123)%>% 
  visEdges(arrows = 'from')%>% 
  visInteraction(navigationButtons = TRUE)

4. 将交互网络上传至公共服务器

首先点击这个像眼睛一样的图标


image

我这里使用的是RPubs

image

然后跟着流程注册账号，设置网页名称即可~

原始project后台联系获取


image.png

感谢观看，如果有用还请多多点赞，关注，在看，转发！

蛋白STRINGdbsectioncodestrong

有关STRINGdb | 交互式可视化蛋白网络互作的更多相关文章

  1. ruby-on-rails - 如何在 ruby​​ 交互式 shell 中有多行？ - 2

    这可能是个愚蠢的问题。但是，我是一个新手......你怎么能在交互式ruby​​shell中有多行代码？好像你只能有一条长线。按回车键运行代码。无论如何我可以在不运行代码的情况下跳到下一行吗？再次抱歉，如果这是一个愚蠢的问题。谢谢。 最佳答案 这是一个例子:2.1.2:053>a=1=>12.1.2:054>b=2=>22.1.2:055>a+b=>32.1.2:056>ifa>b#Thecode‘if..."startsthedefinitionoftheconditionalstatement.2.1.2:057?>puts"f

  2. ruby - Ruby 中的波形可视化 - 2

    我即将开始一个将录制和编辑音频文件的项目，我正在寻找一个好的库(最好是Ruby，但会考虑Java或.NET以外的任何库)以进行实时可视化波形。有人知道我应该从哪里开始搜索吗？ 最佳答案 要流入浏览器的数据量很大。Flash或Flex图表可能是唯一能提高内存效率的解决方案。Javascript图表往往会因大型数据集而崩溃。 关于ruby-Ruby中的波形可视化，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.c

  3. ruby - 用 Ruby 编写一个简单的网络服务器 - 2

    我想在Ruby中创建一个用于开发目的的极其简单的Web服务器(不，不想使用现成的解决方案)。代码如下:#!/usr/bin/rubyrequire'socket'server=TCPServer.new('127.0.0.1',8080)whileconnection=server.acceptheaders=[]length=0whileline=connection.getsheaders想法是从命令行运行这个脚本，提供另一个脚本，它将在其标准输入上获取请求，并在其标准输出上返回完整的响应。到目前为止一切顺利，但事实证明这真的很脆弱，因为它在第二个请求上中断并出现错误:/usr/b

  4. 网络编程套接字 - 2

    网络编程套接字网络编程基础知识理解源`IP`地址和目的`IP`地址理解源MAC地址和目的MAC地址认识端口号理解端口号和进程ID理解源端口号和目的端口号认识`TCP`协议认识`UDP`协议网络字节序socket编程接口`sockaddr``UDP`网络程序服务器端代码逻辑：需要用到的接口服务器端代码`udp`客户端代码逻辑`udp`客户端代码`TCP`网络程序服务器代码逻辑多个版本服务器单进程版本多进程版本多线程版本线程池版本服务器端代码客户端代码逻辑客户端代码TCP协议通讯流程TCP协议的客户端/服务器程序流程三次握手（建立连接）数据传输四次挥手（断开连接）TCP和UDP对比网络编程基础知识

  5. ruby - 检查网络文件是否存在，而不下载它？ - 2

    是否可以在不实际下载文件的情况下检查文件是否存在？我有这么大的(~40mb)文件，例如:http://mirrors.sohu.com/mysql/MySQL-6.0/MySQL-6.0.11-0.glibc23.src.rpm这与ruby​​不严格相关，但如果发件人可以设置内容长度就好了。RestClient.get"http://mirrors.sohu.com/mysql/MySQL-6.0/MySQL-6.0.11-0.glibc23.src.rpm",headers:{"Content-Length"=>100} 最佳答案

  6. ruby - 404 未找到，但可以从网络浏览器正常访问 - 2

    我在这方面尝试了很多URL，在我遇到这个特定的之前，它们似乎都很好:require'rubygems'require'nokogiri'require'open-uri'doc=Nokogiri::HTML(open("http://www.moxyst.com/fashion/men-clothing/underwear.html"))putsdoc这是结果:/Users/macbookair/.rvm/rubies/ruby-2.0.0-p481/lib/ruby/2.0.0/open-uri.rb:353:in`open_http':404NotFound(OpenURI::HT

  7. 深度学习12. CNN经典网络 VGG16 - 2

    深度学习12.CNN经典网络VGG16一、简介1.VGG来源2.VGG分类3.不同模型的参数数量4.3x3卷积核的好处5.关于学习率调度6.批归一化二、VGG16层分析1.层划分2.参数展开过程图解3.参数传递示例4.VGG16各层参数数量三、代码分析1.VGG16模型定义2.训练3.测试一、简介1.VGG来源VGG（VisualGeometryGroup）是一个视觉几何组在2014年提出的深度卷积神经网络架构。VGG在2014年ImageNet图像分类竞赛亚军，定位竞赛冠军；VGG网络采用连续的小卷积核（3x3）和池化层构建深度神经网络，网络深度可以达到16层或19层，其中VGG16和VGG

  8. 【网络】-- 网络基础 - 2

    （本文是网络的宏观的概念铺垫）目录计算机网络背景网络发展认识"协议"网络协议初识协议分层OSI七层模型TCP/IP五层(或四层)模型报头以太网碰撞路由器IP地址和MAC地址IP地址与MAC地址总结IP地址MAC地址计算机网络背景网络发展        是最开始先有的计算机，计算机后来因为多项技术的水平升高，逐渐的计算机变的小型化、高效化。后来因为计算机其本身的计算能力比较的快速：独立模式：计算机之间相互独立。    如：有三个人，每个人做的不同的事物，但是是需要协作的完成。    而这三个人所做的事是需要进行协作的，然而刚开始因为每一台计算机之间都是互相独立的。所以前面的人处理完了就需要将数据

  9. 常见网络安全产品汇总（私信发送思维导图） - 2

    安全产品安全网关类防火墙Firewall防火墙防火墙主要用于边界安全防护的权限控制和安全域的划分。防火墙•信息安全的防护系统，依照特定的规则，允许或是限制传输的数据通过。防火墙是一个由软件和硬件设备组合而成，在内外网之间、专网与公网之间的界面上构成的保护屏障。下一代防火墙•下一代防火墙，NextGenerationFirewall，简称NGFirewall，是一款可以全面应对应用层威胁的高性能防火墙，提供网络层应用层一体化安全防护。生产厂家•联想网御、CheckPoint、深信服、网康、天融信、华为、H3C等防火墙部署部署于内、外网编辑额，用于权限访问控制和安全域划分。UTM统一威胁管理（Un

  10. 【Linux操作系统】——网络配置与SSH远程 - 2

    Linux操作系统——网络配置与SSH远程安装完VMware与系统后，需要进行网络配置。第一个目标为进行SSH连接，可以从本机到VMware进行文件传送，首先需要进行网络配置。1.下载远程软件首先需要先下载安装一款远程软件：FinalShell或者xhell7FinalShellxhell7FinalShell下载：Windows下载http://www.hostbuf.com/downloads/finalshell_install.exemacOS下载http://www.hostbuf.com/downloads/finalshell_install.pkg2.配置CentOS网络安装好

随机推荐