关于用R语言提取KEGG通路基因可不可靠

于 2024-09-13 16:53:49 发布
阅读量170 收藏 2
点赞数 6
文章标签： r语言 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_68600081/article/details/142214935
版权

关于用R语言提取KEGG通路基因可不可靠

一直有个问题就是用R语言提取的KEGG通路的基因到底可不可靠
用的R包是’KEGGREST’

library(KEGGREST)

listDatabases()
keggList("pathway","hsa")
keggGet("hsa00020")

[[1]]$GENE
 [1] "1431"                                                                                            
 [2] "CS; citrate synthase [KO:K01647] [EC:2.3.3.1]"                                                   
 [3] "47"                                                                                              
 [4] "ACLY; ATP citrate lyase [KO:K01648] [EC:2.3.3.8]"                                                
 [5] "50"                                                                                              
 [6] "ACO2; aconitase 2 [KO:K01681] [EC:4.2.1.3]"                                                      
 [7] "48"                                                                                              
 [8] "ACO1; aconitase 1 [KO:K01681] [EC:4.2.1.3]"                                                      
 [9] "3417"                                                                                            
[10] "IDH1; isocitrate dehydrogenase (NADP(+)) 1 [KO:K00031] [EC:1.1.1.42]"                            
[11] "3418"                                                                                            
[12] "IDH2; isocitrate dehydrogenase (NADP(+)) 2 [KO:K00031] [EC:1.1.1.42]"                            
[13] "3420"                                                                                            
[14] "IDH3B; isocitrate dehydrogenase (NAD(+)) 3 non-catalytic subunit beta [KO:K00030] [EC:1.1.1.41]" 
[15] "3421"                                                                                            
[16] "IDH3G; isocitrate dehydrogenase (NAD(+)) 3 non-catalytic subunit gamma [KO:K00030] [EC:1.1.1.41]"
[17] "3419"                                                                                            
[18] "IDH3A; isocitrate dehydrogenase (NAD(+)) 3 catalytic subunit alpha [KO:K00030] [EC:1.1.1.41]"    
[19] "55753"                                                                                           
[20] "OGDHL; oxoglutarate dehydrogenase L [KO:K00164] [EC:1.2.4.2]"                                    
[21] "4967"                                                                                            
[22] "OGDH; oxoglutarate dehydrogenase [KO:K00164] [EC:1.2.4.2]"                                       
[23] "1743"                                                                                            
[24] "DLST; dihydrolipoamide S-succinyltransferase [KO:K00658] [EC:2.3.1.61]"                          
[25] "1738"                                                                                            
[26] "DLD; dihydrolipoamide dehydrogenase [KO:K00382] [EC:1.8.1.4]"                                    
[27] "8802"                                                                                            
[28] "SUCLG1; succinate-CoA ligase GDP/ADP-forming subunit alpha [KO:K01899] [EC:6.2.1.4 6.2.1.5]"     
[29] "8801"                                                                                            
[30] "SUCLG2; succinate-CoA ligase GDP-forming subunit beta [KO:K01900] [EC:6.2.1.4 6.2.1.5]"          
[31] "8803"                                                                                            
[32] "SUCLA2; succinate-CoA ligase ADP-forming subunit beta [KO:K01900] [EC:6.2.1.4 6.2.1.5]"          
[33] "6389"                                                                                            
[34] "SDHA; succinate dehydrogenase complex flavoprotein subunit A [KO:K00234] [EC:1.3.5.1]"           
[35] "6390"                                                                                            
[36] "SDHB; succinate dehydrogenase complex iron sulfur subunit B [KO:K00235] [EC:1.3.5.1]"            
[37] "6391"                                                                                            
[38] "SDHC; succinate dehydrogenase complex subunit C [KO:K00236]"                                     
[39] "6392"                                                                                            
[40] "SDHD; succinate dehydrogenase complex subunit D [KO:K00237]"                                     
[41] "2271"                                                                                            
[42] "FH; fumarate hydratase [KO:K01679] [EC:4.2.1.2]"                                                 
[43] "4190"                                                                                            
[44] "MDH1; malate dehydrogenase 1 [KO:K00025] [EC:1.1.1.37]"                                          
[45] "4191"                                                                                            
[46] "MDH2; malate dehydrogenase 2 [KO:K00026] [EC:1.1.1.37]"                                          
[47] "5091"                                                                                            
[48] "PC; pyruvate carboxylase [KO:K01958] [EC:6.4.1.1]"                                               
[49] "5105"                                                                                            
[50] "PCK1; phosphoenolpyruvate carboxykinase 1 [KO:K01596] [EC:4.1.1.32]"                             
[51] "5106"                                                                                            
[52] "PCK2; phosphoenolpyruvate carboxykinase 2, mitochondrial [KO:K01596] [EC:4.1.1.32]"              
[53] "5161"                                                                                            
[54] "PDHA2; pyruvate dehydrogenase E1 subunit alpha 2 [KO:K00161] [EC:1.2.4.1]"                       
[55] "5160"                                                                                            
[56] "PDHA1; pyruvate dehydrogenase E1 subunit alpha 1 [KO:K00161] [EC:1.2.4.1]"                       
[57] "5162"                                                                                            
[58] "PDHB; pyruvate dehydrogenase E1 subunit beta [KO:K00162] [EC:1.2.4.1]"                           
[59] "1737"                                                                                            
[60] "DLAT; dihydrolipoamide S-acetyltransferase [KO:K00627] [EC:2.3.1.12]"

最后得到这30个基因，那到底可不可靠呢，就让我们回到KEGG网站吧
在这里插入图片描述

我进行了逐个比较发现其实都是有的，比如1.1.1.41在分析结果有三个，但这里只有1个，点进去之后就能找到3个了
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
实在不行就用使用KEGG API进行对比，也可以看到这30个基因都在
在这里插入图片描述

所以用R包提取KEGG通路的基因是可靠的，其实这个包使用的就是KEGG的API接口

reference
434383719#:~:text=提取某条KEGG

Jorth_1
关注
R-下载某一条通路的所有基因名字（KEGG）
upyours00的博客
06-04 1万+
实例：如何拿到KEGG数据库中多巴胺通路相关的基因集 一、确定目标通路 打开KEGG选择pathway，在搜索框前输入物种，框内填入关键词。 筛选结果显示仅有hsa04728符合我们的研究目的 二、下载hsa04728通路中的全部基因 ##1.安装R包KEGGREST 首次安装时电脑可能会显示与当前R语言版本不配，可以从bioconductor 的官网下载安装 if(!requireNamespace("BiocManager", quietly = TRUE)) install.pac
R语者小case之——从KEGG原始网页批量生成通路的基因表格
weixin_49352797的博客
07-15 1039
R语者小case——从差异表达的结果读取GSEA的input 欢迎使用Markdown编辑器 你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。 新的改变 我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客： 全新的界面设计 ，将会带来全新的写作体验； 在创作中心设置你喜爱的代码高亮样式，Mark
差异基因富集分析（R语言——GO&KEGG&GSEA）
m0_73145613的博客
07-03 2494
通路富集分析气泡图
如何获取kegg通路的基因列表
TOP生物信息
02-07 2916
kegg基因列表
利用KEGGREST包提取所有代谢相关的基因
arrhythmia08的博客
05-19 2504
KEGG上的基因数量和通路数量都是不断更新的，所以KEGGREST包能提取的基因数量也在不断变化。我去年运行的时候总基因8192个，现在已经有8443个了，以后也会不断增加。KEGG官网上也有少量的代谢通路不以00开头，以01开头。数量太少，且不是重要的通路，所以我省略了。想要提取更全的基因基的朋友可以自己加上。写笔记主要是为了自己进步，也是方便自己查看。R语言新手，如有不对请指正。
获取KEGG通路的基因列表 做单细胞GSEA、GSVA分析
生信小博士的博客
03-22 2337
这通常涉及使用专门的R包，如KEGGREST或biomaRt，来查询KEGG数据库并检索特定通路的基因列表。这包括加载单细胞RNA-seq数据，通常使用Seurat或其他单细胞分析包进行预处理。：使用GSVA包对单细胞数据执行基因集变异分析（GSVA），根据KEGG通路的基因列表评估每个单细胞样本的通路活性。果：最后，基于GSVA分析结果，绘制热图或其他类型的图表来展示不同单细胞样本中通路活性的变化。
获得kegg通路中所有基因的名称，并将他们合并去重
BELLEsang的博客
12-12 963
获取kegg基因名称，Entry号是hsa04010。获得csv文件后使用python进行合并和去重。# 将合并后的数据保存到一个新的CSV文件。# 创建带有去重基因的DataFrame。# 保存去重后的数据到一个新的CSV文件。# 将四个数据框的第一列连接起来。# 将列表输出为CSV文件。#获取通路中gene信息。# 读取合并后的数据。
【R语言】——基因GO/KEGG富集分析!超级简单的保姆级教程！
热门推荐
weixin_54004950的博客
01-09 3万+
GO/KEGG功能富集分析中重要的是背景基因的选择，使用R clusterProfiler包对基因进行富集，需要导入目的基因（前景基因）相对应物种的参考基因组（背景基因），现阶段“bioconductor”已有十几种常见动物，如人类、小鼠等物种的OrgDb。”介绍如何使用DAVID在线分析工具对基因进行GO/KEGG功能富集分析和使用R ggplot包对获得的基因GO/KEGG功能富集结果进行可视化。：富集到该GO term/KEGG term中的基因数目/给定基因的总数目；
如何使用R语言优雅地解析KEGG KO 的htext文件？
xiaoafei1234的博客
07-24 961
于是小编上网搜各种分析流程和教程之类的，也算是找到几个python或者perl脚本之类的，但是使用过后，发现要么就是出错，要么分析出来的效果不尽人意，所以小编打算自己动手写代码进行分析。list格式也算是R里面最复杂的格式，不过也没那么复杂，只要仔细观察，你就可以找到这个list里面的一些规律，有兴趣的朋友可以自己点进去看看，这里我不多说了，直接上代码。​​​​​这段代码就是解析json文件的核心代码了，接下来的工作大家也可以使用自己的需求进行写，小编在这里演示一下自己的方法。这玩意要怎么用来数据分析？
R语言绘图.pdfR语言绘图.pdf
05-06
- 生物信息专业图的绘制，如KEGG通路图、GO富集图，可使用`ggplot2`或专门的软件如`BiocManager`。 5. **绘图软件结合**： - Photoshop和Illustrator可以用于后期图片的调整和整合，提升图像质量。 - ...
R语言|clusterprofile超几何分布富集分析 GO,KEGG富集分析，循环Fisher‘s test
mocchan925的博客
12-24 4241
超几何分布富集分析 GO,KEGG富集分析，循环Fisher‘s testID转换GOKEGG把KEGG里的geneid转回名称（readable table）超几何分布以及生成data frame生成用来循环fisher的table循环Fisher's test 这几天真的被这个超几何富集分析做的头大。原来单纯用找到的显著基因分析出来的pathway不一定是准确的，或者说不一定是具有显著性的。详细的原理推荐下面两个链接 https://www.jianshu.com/p/3d01a66e235b htt
GO/KEGG/自定义通路的富集方法
weixin_44070677的博客
09-09 6220
提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档 目录 文章目录 前言 一、clusterProfiler 进行GO、KEGG、Canonical Pathways 二、使用步骤 1.加载所需的R包 2.读入数据与处理表格 3.enrichGO函数进行GO富集 前言 clusterProfiler 是业界大神Y叔写的一个R包，可以用来做各种富集分析，如GO、KEGG、以及GSEA富集分析等，并且对富集分析结果进行可视化。 这里将使用clusterProfiler包对
R获取指定GO term和KEGG pathway的gene list基因集
weixin_30274627的博客
05-18 5188
clusterProfiler没有显性的接口，但是可以直接扣取clusterProfiler里的函数。 核心函数就是get_GO_data GO_DATA <- get_GO_data("org.Hs.eg.db", "BP", "SYMBOL") 　　 可以看到输入的是GO数据库，选定类别，基因名字类型，输出的就是整个数据库。 但是想调用这个函数没那么简单，得导入一系列...
【scMetabolism】一站式代谢通路基因集打分
DerekRusso的博客
04-26 2104
scMetabolism
R语言论文插图模板第9期—滑珠散点图
最新发布
qq_26447137的博客
09-11 363
R语言论文插图模板第9期—滑珠散点图
YOLOv5改进 | 模块缝合 | C3 融合RFAConv和CBAM注意力机制 【二次融合 小白必备】
kay_545
09-10 1340
YOLOv5，YOLOv5改进，YOLOv5模块缝合
Autosar(Davinci) --- 创建一个S/R类型的port（下）
wj博客
09-07 220
前面章节我们讲解了S/R类型的Port如何创建，这一章节，我们着重讲一下生成的代码，以及我们如何添加代码让这些门与灯之间的关系产生连接。
嵌入式鸿蒙系统开发语言与开发方法分析
嵌入式技术开发
09-10 1086
可复用的UI单元，可组合其它组件，如上述被@Component装饰的struct Hello。UI描述:声明式的方法来描述UI的结构，例如build()方法中的代码块。
r语言kegg通路富集分析
09-16
R语言中的KEGG通路富集分析是一种用于研究基因集合在KEGG通路中富集况的方法。首先，需要准备好R环境和所需的包，比如"clusterProfiler"和"pathview"。然后，将基因表达数据读入并进行相应的数据处理，如选择感兴趣的基因，并为每个基因赋予ID。接下来，使用enrichKEGG函数进行KEGG富集分析，设置所需的参数，比如基因集合、物种和显著性阈值。分析结果可以写入文件中以供后续使用。最后，使用pathview函数生成每个KEGG通路的可视化图形。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值