Converting gene id to entrez id with clusterprofiler
0
0
Entering edit mode
9 months ago
sumitra.20 • 0

Hi,

Im new to R and i've been trying to use clusterprofiler to convert my gene_ids to entrez_id for further analysis. I work on the msu7 rice genome and my gene ids looks like this (LOC_Os08g09610, LOC_Os06g05660). Since, it was a non-model organism i tried creating new annotation and ended up with an error. I might be doing a very silly mistake here, but i get error that says there is no package in the path. Also, im not really sure if the kind of gene_id i have will be applicable to these functions, my final goal is to eventually perform KEGG enrichment. Any help will be appreciated. (On a side note, i tried posting this on bioconducter but for some reason im not able to post any questions there, not sure what's wrong)

```

library(AnnotationHub)
ah = AnnotationHub::AnnotationHub()  
SL<-AnnotationHub::query(ah, c("Oryza", "Japonica"))
SL_OrgDb <- SL[["AH91623"]]


library(clusterProfiler)
genes2<-read.table("./deg_gene_list.txt")
geneList<-genes2
geneList1<-(geneList$V1)  

gene.df <- bitr(geneList1, fromType = c("GENENAME"),
                toType = "ENTREZID",
                OrgDb = SL_OrgDb)

Error in loadNamespace(name) : 
  there is no package called ‘/home/combio7/.cache/R/AnnotationHub/92e64320d6b3d_98369’


> SL_OrgDb
                                                   AH91623 
"/home/combio7/.cache/R/AnnotationHub/92e64320d6b3d_98369" 
msu7 clusterProfiler annotationhub oryza_sativa • 547 views
ADD COMMENT
0
Entering edit mode

Based on locus search here: http://rice.uga.edu/analyses_search_locus.shtml this looks like Oryza sativa nipponbare genome. NCBI's official genome page is here: https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_001433935.1/

ADD REPLY

Login before adding your answer.

Traffic: 2776 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6