Trouble using getBM with ecaballus_gene_ensembl
1
0
Entering edit mode
@stephenwilliams-15198
Last seen 6.3 years ago

I'm having trouble querying biomart using the biomaRt package for a non human species. for example, this works

library('biomaRt')
mart <- useMart("ensembl",dataset="hsapiens_gene_ensembl")
getBM(attributes = c('entrezgene','hgnc_symbol', 'ensembl_gene_id'), 
       filters = 'hgnc_symbol', 
       values = 'RAI1', 
       mart = mart)
entrezgene hgnc_symbol ensembl_gene_id
1      10743        RAI1 ENSG00000108557

However when using the ecaballus_gene_ensembl mart the same thing does not work

library('biomaRt')
horse = useMart("ensembl", dataset = "ecaballus_gene_ensembl")
getBM(attributes = c('entrezgene','hgnc_symbol', 'ensembl_gene_id''), 
       filters = 'hgnc_symbol', 
       values = 'Rai1', 
       mart = horse)
[1] entrezgene      hgnc_symbol     ensembl_gene_id
<0 rows> (or 0-length row.names)

I have confirmed that this gene does exist in horse and that all the fields and filters do as well (confirmed using listFilters(horse)) so I'm very confused. Any help would be greatly appreciated. 

biomart biomartr • 1.2k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 2 minutes ago
United States

Does this help?

> getBM(attributes = c('entrezgene','hgnc_symbol', 'hgnc_id','ensembl_gene_id','external_gene_name'), filters = 'ensembl_gene_id', values = 'ENSECAG00000020530', mart = horse)
  entrezgene hgnc_symbol hgnc_id    ensembl_gene_id external_gene_name
1  100050870          NA      NA ENSECAG00000020530               RAI1
ADD COMMENT
0
Entering edit mode

Yes! This helps a lot. I didn't know about "external_gene_name".

ADD REPLY

Login before adding your answer.

Traffic: 1018 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