This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Different Results With Goose And Biomart

I'm using the biomaRt package for R in order to get all the GO terms associated with with a particular Gene ID (Entrez ID). The script is quite simple and similar to the tutorial example:

ensembl = useMart("ensembl",dataset="hsapiens_gene_ensembl")
goids = getBM(attributes=c('entrezgene','go_id') ,filters='entrezgene', values=geneList, mart=ensembl)

This gives me a list of GO terms for each gene in geneList. In order to check these results I queried the GO dataset using GOOSE. So I translated Entrez ID into Uniprot ID and I executed a simple basic query (SQL Query by Uniprot ID). As result, usually the list obtained with biomaRt is longer than the one obtained using GOOSE (which directly query to GO dataset). How is this possible? Which method should I trust more? I think that GOOSE is more reliable, but why biomaRt list more GO terms?

biomart go r
• 3,030 views • link
updated by secretjess secretjess • written by TitoPullo TitoPullo

Can you please give an example?

• 1 views • link
written by Bert Overduin Bert Overduin

1 answer

I'm wondering if it's something to do with the multiple products from a single gene or perhaps the mapping of GO terms to Entrez/Uniprot IDs? If you repeat your Ensembl query using "uniprot_swissprot" as your filter do you get the same result?

• 0 views • link
written by secretjess secretjess

Log in to answer this question.

More posts like this