This is a test version of Biostars. For the public version, visit https://www.biostars.org.
can not find certain genes: Ensembl_Gene_Id To Entrezgene using Biomart

Hi, this is probably a silly question. I'm using getBM function to convert gene name to Entrezgene ID, but find that not all genes converted successfully. For those that didn't get their Entrezgene ID, I checked on NCBI website and they do have a Entrezgene ID.

Any advice/suggestions are appreciated. Here is my code: (I'm using biomaRt version 2.50.3 and R version 4.1.3)

test<-c("ABI3","CCL3L1","GNG5")

mapping <- getBM(
attributes =  c('entrezgene_id', 'hgnc_symbol'), 
filters = 'hgnc_symbol', values = test, mart = hsmart)

the first gene ABI3 (ID=51225) is good, the last two genes CCL3L1 (ID=6349) and GNG5 (ID=2787) cannot be assigned to their Entrezgene ID.

enter image description here

gene r biomart id
• 2,092 views • link
updated by GenoMax GenoMax • written by catherine catherine

Do you see any warnings?

• 0 views • link
written by Ram Ram

no any warnings

• 0 views • link
written by catherine catherine

Using EntrezDirect: so the ID's exist.

$ esearch -db gene -query "CCL3L1 [gene] and 9606 [taxID]" | esummary | xtract -pattern DocumentSummary -element Id
6349

$ esearch -db gene -query "GNG5 [gene] and 9606 [taxID]" | esummary | xtract -pattern DocumentSummary -element Id
2787

$ esearch -db gene -query "ABI3 [gene] and 9606 [taxID]" | esummary | xtract -pattern DocumentSummary -element Id
51225
• 0 views • link
written by GenoMax GenoMax

It could be that the ensembl version is outdated. OP, can you show the code that creates the hsmart object please?

• 0 views • link
written by Ram Ram

thank you, but it could be better if using R

• 0 views • link
written by catherine catherine 
hsmart <- useMart(dataset = "hsapiens_gene_ensembl", biomart = "ensembl")
• 0 views • link
written by catherine catherine

0 answers

No answers yet.

Log in to answer this question.

More posts like this