Search
Question: Reactome_gene IDs are pathway IDs in biomaRt
0
gravatar for b.nota
12 months ago by
b.nota300
Netherlands
b.nota300 wrote:

Hello,

I've been trying to use biomaRt to get reactome IDs (from pathways) of given reactome_gene IDs, however I found out that the reactome_gene IDs ARE the pathway IDs. It is now not possible to get the matching pathways by these gene IDs.

 

For example

library(biomaRt)

ensembl = useMart("ensembl",dataset="hsapiens_gene_ensembl")

reactomeID <- c("R-HSA-49155", "R-HSA-199420")

getBM(attributes=c('reactome_gene', 'reactome'),
      filters = 'reactome_gene',
      values = reactomeID,
      mart = ensembl)

[1] reactome_gene reactome
<0 rows> (or 0-length row.names)

getBM(attributes=c('reactome_gene', 'entrezgene'),
      filters = 'reactome_gene',
      values = 'R-HSA-49155',
      mart = ensembl)

[1] reactome_gene reactome
<0 rows> (or 0-length row.names)

However, when I use a pathway (reactome) ID instead of reactome_gene ID it does work, but this is not what I am asking to find.

getBM(attributes=c('reactome_gene', 'entrezgene'),
      filters = 'reactome_gene',
      values = 'R-HSA-110331',
      mart = ensembl)

 reactome_gene entrezgene
1  R-HSA-110331       4968
2  R-HSA-110331       4595
3  R-HSA-110331       4350


Is this a bug? Or am I missing something?

 

 

 

 

 

ADD COMMENTlink modified 12 months ago • written 12 months ago by b.nota300
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 121 users visited in the last hour