Omics - Hunter

利用clusterProfiler实现蛋白名称、基因名称等ID快速转换

clusterProfiler 是一个非常好用的R包,可以用来做各种富集分析,如GO、KEGG、DO(Disease Ontology analysis)、Reactome pathway analysis以及GSEA富集分析等。而除了富集分析,它还可以非常方便的做蛋白、基因等ID之间的转换。

# 安装R包
if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")

BiocManager::install("clusterProfiler")
BiocManager::install("org.Hs.eg.db")

以人的蛋白、基因名之间的转换为例:

# 加载R包
library(clusterProfiler)
library(org.Hs.eg.db)

查看org.Hs.eg.db数据库中支持多少类ID的转换:

以Uniport ID转换为例

gene <- bitr(
  c("P04114", "P01024"),
  fromType = "UNIPROT",
  toType = c("UNIPROT", "ENSEMBL", "SYMBOL", "ENTREZID"),
  OrgDb = org.Hs.eg.db
)

当然,我们也可以尝试其他ID的转换,本文点到为止。

参考文章:

1.http://www.bioconductor.org/packages/release/bioc/vignettes/clusterProfiler/inst/doc/clusterProfiler.html


作者:陈浩


版权:本文版权归作者所有


免责声明:本文中使用的部分图片来自于网络或者参考资料,如有侵权,请联系博主:chenhao__@__evvail.com(发件请删除下划线)进行删除


转载注意:除非特别声明,本站点内容均为作者原创文章,转载须以链接形式标明本文链接


本文链接:https://evvail.com/2020/04/15/411.html

3 评论

  1. 比较好奇 我该怎么把转换好的ID作为一列加入原来的表达矩阵呢?

  2. 陈浩

    2020/5/8 在 19:55

    bitr用法,直接可以访问:https://yulab-smu.github.io/clusterProfiler-book/chapter14.html#bitr

发表回复

如果你有什么好的建议或者疑问请给我留言,谢谢!

Captcha Code