蛋白组科研信息学基因家族分析

蛋白ID转基因ID

2020-12-28  本文已影响0人  生信编程日常

将Ensembl 中的蛋白ID转化成基因ID,可以通过clusterProfiler这个包。
如以大鼠的基因与蛋白转化为例;
安装clusterProfiler与大鼠org.Rn.eg.db,如果是人的注释包为org.Hs.eg.db,小鼠的注释包为org.Mm.eg.db.

BiocManager::install("clusterProfiler")
BiocManager::install("org.Rn.eg.db")
library(clusterProfiler)
library(org.Rn.eg.db)

查看可以转化的ID:

keytypes(org.Rn.eg.db)

[1] "ACCNUM" "ALIAS" "ENSEMBL" "ENSEMBLPROT" "ENSEMBLTRANS"
[6] "ENTREZID" "ENZYME" "EVIDENCE" "EVIDENCEALL" "GENENAME"
[11] "GO" "GOALL" "IPI" "ONTOLOGY" "ONTOLOGYALL"
[16] "PATH" "PFAM" "PMID" "PROSITE" "REFSEQ"
[21] "SYMBOL" "UNIGENE" "UNIPROT"

将蛋白ID转为基因ID:

PROT2SYMBOL<-bitr(data$Protein.Group.Accessions,fromType = "ENSEMBLPROT",toType = "SYMBOL",OrgDb = org.Rn.eg.db)
head(PROT2SYMBOL)
上一篇 下一篇

猜你喜欢

热点阅读