首页 > 数据库 >你敢相信,我5分钟做了模式植物的GO和KEGG富集分析,并创建了orgDb数据库

你敢相信,我5分钟做了模式植物的GO和KEGG富集分析,并创建了orgDb数据库

时间:2024-08-12 19:59:40浏览次数:13  
标签:分析 富集 KEGG WGCNA ## 代码 教程 GO orgDb

原文教程:我5分钟做了模式植物的GO和KEGG富集分析,并创建了orgDb数据库

一边学习,一边总结,一边分享!

本期教程

GO富集

获得本教程Data and Code,请在后台回复:20240811

2022年教程总汇

https://mp.weixin.qq.com/s/Lnl258WhbK2a8pRZFuIyVg

2023年教程总汇

https://mp.weixin.qq.com/s/wCTswNP8iHMNvu5GQauHdg

写在前面

我们在前面的教程分享了模式植物构建orgDb数据库 | 以org.Slycompersicum.eg.db为例,这两天也做对应的分析,但是发现前面的代码速度太慢了,以及会出现报错的情况。

那么我们也做了对应的优化,可以大大减少花费时间,以及减少我们的报错。

现在我们的代码,可以直接读取eggnog分析获得数据,非常的方便。

此外,我们在这里依旧推荐大家使用本地的Egg-mapper进行注释,我自己使用下来,非常的快捷。虽然网页版的Egg-mapper工具,现在可以支持最多10000条的序列,但是速度依旧是比较慢的。

此外,本地使用Egg-mapper使用时,数据库的下载可能对新手小白具有挑战性,若是本推文阅读量超过5000+,我们也会视频讲解+文本,手把手教大家。


Code

本次,我们依旧使用我们的一步法进行数据库的构建。

source("../Set_OrgDb_Database.R")

# 使用函数
# 使用函数
Set_OrgDb_Database(
  emapper_file = "tomato.emapper.annotations",  ## 输入的eggnog结果文件
  json_file = "../ko00001.json",               ## 下载ko00001.json,下载网址:https://www.genome.jp/kegg-bin/get_htext?ko00001
  tax_id = "4081",                          ## 物种信息
  genus = "Solanum",
  species = "lycompersicum",
  versions = "4.0",                         ## 版本号
  maintainer = "du<16***@qq.com>",     ## 修改为你的名字和邮箱
  author = "du<16****@qq.com>",         ## 修改为你的名字和邮箱
  outputDir = "."## 保存路径
)

我们自己这边注注释番茄黄瓜甜瓜辣椒的数据库,基本每个数据库会在2分钟内完成。

加载数据包并进行分析

install.packages("org.Slycompersicum.eg.db/", repos = NULL,
                 type = 'sources')

library(org.Slycompersicum.eg.db)

GO 富集分析

Gene_list <- read.table("01.tomato_target_gene.id.txt",header = F)$V1
head(Gene_list)

go_result <- enrichGO(gene = Gene_list,
                      OrgDb = org.Slycompersicum.eg.db,
                      keyType = "GID",
                      ont = "ALL",
                      qvalueCutoff = 1,
                      pvalueCutoff = 1)

go_data <- as.data.frame(go_result)

dim(go_data)
head(go_data)

write.csv(go_data, "02.tomato.miRNA.GO.result.csv")

## 绘图
pdf("04.tomato.go.柱状图.pdf",width = 8, height = 10)
barplot(go_result,
        drop = TRUE,
        showCategory = 20)
dev.off()

KEGG 富集分析

##KEGG 富集分析
pathway2gene <- AnnotationDbi::select(org.Slycompersicum.eg.db,
                                      keys = keys(org.Slycompersicum.eg.db),
                                      columns = c("Pathway","Ko")) %>%
  na.omit() %>%
  dplyr::select(Pathway, GID)

必须加载的数据,在数据包制作过程中会生成。

load("kegg_info.RData")
kegg_result <- enricher(Gene_list,
                        TERM2GENE = pathway2gene, 
                        TERM2NAME = pathway2name,
                        pvalueCutoff = 1,
                        qvalueCutoff = 1)

kegg_data <- as.data.frame(kegg_result)
dim(kegg_data)

write.csv(kegg_data, "./03.tomato.kegg.result.csv")

pdf("05.kegg.tomato.柱状图.pdf",height = 4, width = 6)
barplot(kegg_result, showCategory= 10,
        drop=F,
        color="pvalue", 
        font.size=10) 
dev.off()

富集网络图

MF_ego <- enrichGO(
  gene  = gene_list,
  keyType = "GID",
  OrgDb   = org.Slycompersicum.eg.db,
  ont     = "MF", ##'@GO的种类,BP,CC,MF
  pAdjustMethod = "BH",
  pvalueCutoff  = 0.5##'@结合自己的需求进行调整
)
pdf("./功能富集/05.MF_网络图.pdf", width = 8, height = 8)
plotGOgraph(MF_ego)
dev.off()

MF_ego <- pairwise_termsim(MF_ego)
pdf("./功能富集/06.网络图.pdf",width = 8, height = 8)
emapplot(MF_ego, cex.params = list(category_label = 0.8, line = 0.5)) + 
  scale_fill_continuous(low = "#e06663", high = "#327eba", name = "p.adjust",
                        guide = guide_colorbar(reverse = TRUE, order=2), trans='log10')
dev.off()


写在最后

分析就是不停的折腾,不停的捣鼓。小杜在这里记录是自己分析中所学的、用的,以及存在的问题。我们的代码可能在运行时候会出现不同的问题,也会出现,在我这里正常运行,而在你那边程序中却报错的情况。

对于我自己而言,也会出现这样的情况。自己的代码或教程,在前面分析过程中可以运行,但是现在再去使用时也会报错。

好比,这个教程,我昨天的运行时,也各种报错,内心是发狂的,不知道为什么。最后,仅仅只是更改了某个地方的参数,导致在制作过程中缺失变量。

获得本教程Data and Code,请在后台回复:20240811

若我们的教程对你有所帮助,请点赞+收藏+转发,这是对我们最大的支持。

往期部分文章

1. 最全WGCNA教程(替换数据即可出全部结果与图形)


2. 精美图形绘制教程

3. 转录组分析教程

4. 转录组下游分析

小杜的生信筆記 ,主要发表或收录生物信息学教程,以及基于R分析和可视化(包括数据分析,图形绘制等);分享感兴趣的文献和学习资料!!

标签:分析,富集,KEGG,WGCNA,##,代码,教程,GO,orgDb
From: https://blog.csdn.net/kanghua_du/article/details/141140342

相关文章

  • 从AlphaGo到AlphaPong:AI乒乓球机器人将如何改变乒乓球竞技?
    在人工智能领域,谷歌DeepMind的AlphaGo曾以其超凡的围棋能力震惊了世界。而如今,DeepMind又将其技术延伸到了乒乓球领域,推出了乒乓球机器人,被网友戏称为“AlphaPong”。这一机器人使用ABB的IRB1100机械手臂,能够与普通人进行激烈的乒乓对抗。那么,中国在乒乓球领域的传统优势是......
  • How to change the google translate voice to male All In One
    HowtochangethegoogletranslatevoicetomaleAllInOneTL;DR:Thereisnowaytodothatfornow./Thereiscurrentlynowaytodothis.https://support.google.com/translate/thread/2017149/is-there-a-male-voice-for-google-translator-where-can-i-s......
  • 计算机毕业设计django+vue代驾服务【开题+论文+程序】
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着城市化进程的加速和汽车保有量的持续增长,代驾服务作为一种便捷、安全的出行方式,逐渐受到广大消费者的青睐。然而,传统的代驾服务模式往......
  • 计算机毕业设计django+vue民宿预定管理系统625l0【开题+论文+程序】
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着互联网技术的迅猛发展和旅游业的不断壮大,民宿作为旅游住宿的重要组成部分,其市场需求日益增加。然而,传统的民宿管理方式已难以满足日益......
  • Django 解除跨域限制
    Django解除跨域限制在Web开发中,跨域资源共享(CORS)是一个重要的安全特性,它限制了网页只能与其同源的服务器进行交互。然而,在开发过程中,我们经常需要前端(如Vue.js、React等)与后端(如Django)之间进行跨域请求。为了方便开发和测试,我们需要在Django中解除跨域限制。本文将详细介绍如何......
  • 基于django+vue基于web的园区车辆出入管理系统【开题报告+程序+论文】计算机毕设
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着城市化进程的加快,各类产业园区、住宅小区及商业综合体等园区规模不断扩大,车辆管理成为园区管理中的重要环节。传统的车辆出入管理方式......
  • 三十分钟入门基础Go(Java小子版)
    前言Go语言定义Go(又称Golang)是Google的RobertGriesemer,RobPike及KenThompson开发的一种静态、强类型、编译型语言。Go语言语法与C相近,但功能上有:内存安全,GC,结构形态及CSP-style并发计算。适用范围本篇文章适用于学习过其他面向对象语言(Java、Php),但没有学过......
  • Django项目与Vue的集成
    Django项目与Vue的集成在现代Web开发领域,前后端分离已成为一种主流趋势。Django,作为一个强大的PythonWeb框架,以其丰富的功能和高度的可扩展性而受到开发者的青睐。而Vue.js,作为一个轻量级的渐进式JavaScript框架,以其简洁的API和灵活的组件系统成为前端开发的热门选择。将Django......
  • 159.302 The 8-Puzzle: Search Algorithms
    159.302ArtificialIntelligenceAssignment#1The8-Puzzle:SearchAlgorithmsMaximumnumberofmemberspergroup:3studentsDeadlineforsubmission:9thofSeptemberInstructionsYourtaskistowriteaC++programthatwillsolvethe8-puzzleprob......
  • 记录JSch连接SFTP Exception:Algorithm negotiation fail问题解决
    问题描述:关于正式环境访问外网连接不成功 1、首先检查是否开放防火墙(已确认开放),策略开放后,通过命令连接是否畅通: 通过telnet命令,可以得出,访问畅通。telnet192.168.1.122 2、查看生产环境日志,观察生产环境访问外网服务器异常:抛出异常,提示:算法协商失败com.jcraft.j......