This is a test version of Biostars. For the public version, visit https://www.biostars.org.
annotation - biomaRt - getBM - multiple entrez ID

Hi

I am trying to follow what is said in here , however, without any success

https://support.bioconductor.org/p/52407/

when I do it as described here, I am fine 

ex <- c("ENSG00000215417", "ENSG00000224078", "ENSG00000198366",
"ENSG00000196176", "ENSG00000166012", "ENSG00000158406",
"ENSG00000196787")
mart <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene2genomeEx <- getBM(values = ex, filters = "ensembl_gene_id", mart
= mart, attributes = c("ensembl_gene_id", "entrezgene","hgnc_symbol",
"external_gene_id", "external_gene_db", "description",
"chromosome_name", "strand"))

however, when I try to use gene name, it crashes and gets me error 

ex <- c("ACTN4","TUBA1B","ACTN1","TP53")
mart <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
    gene2genomeEx <- getBM(values = ex, filters = "external_gene_name", mart
    = mart, attributes = c("external_gene_name", "entrezgene","hgnc_symbol",
    "external_gene_id", "external_gene_db", "description",
    "chromosome_name", "strand"))

any thought ?

genome biomart entrez
• 12,018 views • link
updated by Emily Emily • written by Learner Learner

1 answer

The problem is that there are no attributes called external_gene_id or external_gene_db.

Take a look:

require(biomaRt)

Look-up ENSEMBL gene IDs:

ex <- c("ENSG00000215417", "ENSG00000224078", "ENSG00000198366",
"ENSG00000196176", "ENSG00000166012", "ENSG00000158406",
"ENSG00000196787")
mart <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene2genomeEx <- getBM(values = ex,
  filters = "ensembl_gene_id",
  mart = mart,
  attributes = c("ensembl_gene_id", "entrezgene_id",
    "hgnc_symbol", "external_gene_name",
    "description", "chromosome_name",
    "strand"))
gene2genomeEx

  ensembl_gene_id entrezgene hgnc_symbol external_gene_name
1 ENSG00000158406       8365    HIST1H4H           HIST1H4H
2 ENSG00000166012      79101       TAF1D              TAF1D
3 ENSG00000196787       8969   HIST1H2AG          HIST1H2AG
4 ENSG00000215417     407975     MIR17HG            MIR17HG
5 ENSG00000224078       3653      SNHG14             SNHG14

Look-up 'external' gene names:

ex <- c("ACTN4","TUBA1B","ACTN1","TP53")
mart <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene2genomeEx <- getBM(values = ex,
  filters = "external_gene_name",
  mart = mart,
  attributes = c("external_gene_name", "entrezgene_id",
    "hgnc_symbol", "description",
    "chromosome_name", "strand"))
gene2genomeEx

  external_gene_name entrezgene hgnc_symbol
1               TP53       7157        TP53
2              ACTN4         81       ACTN4
3             TUBA1B      10376      TUBA1B
4              ACTN1         87       ACTN1
5              ACTN4         81       ACTN4
                                            description chromosome_name strand
1 tumor protein p53 [Source:HGNC Symbol;Acc:HGNC:11998]              17     -1
2     actinin alpha 4 [Source:HGNC Symbol;Acc:HGNC:166]              19      1
3  tubulin alpha 1b [Source:HGNC Symbol;Acc:HGNC:18809]              12     -1
4     actinin alpha 1 [Source:HGNC Symbol;Acc:HGNC:163]              14     -1
5     actinin alpha 4 [Source:HGNC Symbol;Acc:HGNC:166]  CHR_HG26_PATCH      1

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

@Kevin Blighe Kevin, how did you check for it? I checked the attributes but with no success. is there a way to unique them based on GO id?

• 0 views • link
written by Learner Learner

This is the first time that you mention GO id. What do you mean? Please try to explain in detail the issue that faces you.

• 0 views • link
written by Kevin Blighe Kevin Blighe

@Kevin Blighe instead the gene_name, please now try this one, you will see that you get error

ex <- c("GO:0000002","GO:0042254","GO:0042254","GO:0000022","GO:0000028","GO:0000028","GO:0000045")
• 0 views • link
written by Learner Learner

If you're trying to run exactly the same query, but providing a list of GO ids then it expected you will get an error (or at least no results). The filters argument defines the column in the database you want to search. If you're asking to search a column containing gene IDs or symbols, but looking for GO ids (which are completely different from gene names) then you wont find any results.

You'll need to include the full code that produced the error, otherwise it's very hard to see what exactly you're trying to do.

• 0 views • link
written by Mike Smith Mike Smith

@Mike Smith I get this error but I don't know the argument of 'listFilters'

Error in getBM(values = ex, filters = "go_id", mart = mart, attributes = c("go_id",  : 
  Invalid filters(s): go_id 
Please use the function 'listFilters' to get valid filter names
• 0 views • link
updated by Ram Ram • written by Learner Learner

To filter by GO IDs you have to use the filter for GO IDs. If you're just looking to get genes associated with those exact terms, use go, if you want to get the genes associated with all their child terms too, use go_parent_term.

• 0 views • link
written by Emily Emily

Hi Kevin, I've tried to copy paste exactly what you did 

ex <- c("ENSG00000215417", "ENSG00000224078", "ENSG00000198366",
"ENSG00000196176", "ENSG00000166012", "ENSG00000158406",
"ENSG00000196787")
mart <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene2genomeEx <- getBM(values = ex,
  filters = "ensembl_gene_id",
  mart = mart,
  attributes = c("ensembl_gene_id", "entrezgene",
    "hgnc_symbol", "external_gene_name",
    "description", "chromosome_name",
    "strand"))
gene2genomeEx

But it does't like the attrubute entrezgene:

Error in getBM(values = ex, filters = "ensembl_gene_id", mart = mart,  : 
  Invalid attribute(s): entrezgene 
Please use the function 'listAttributes' to get valid attribute names

Do you have a suggestion of why?

thank you

• 0 views • link
written by Morris_Chair Morris_Chair

Thanks Mike. I noticed this change myself and will update old code where needed, including my answer here (above).

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this
  • UniProt ID, genes
    written by Rob • Rob
    Hi friends, How can I download the UniProt ID for protein coding genes? I have about 20000 genes and it is hard to just search …
  • Missing genes in converting "ensembl_gene_id_version" to "ensembl_gene_id" in biomaRt package
    written by Maryam • Maryam
    Hi, I am using biomaRt package to convert "ensembl_gene_id_version" to "ensembl_gene_id" but unfortunately I have some missing genes. The empty values for the rest are …
  • gene ID RNAseq
    written by Rob • Rob
    Hi friends How can I get gene numeric ID and hugo ID by R script? what script should I use? I have this but does …
  • gene ensemble ID- no symbol
    written by Rob • Rob
    Hello friends, I used this code to convert ensemble gene ID to gene symbol: df <- read.csv("gene_Ens_ID.csv") library('biomaRt') mart <- useDataset("hsapiens_gene_ensembl", useMart("ensembl")) attributeNames <- c('ensembl_gene_id','external_gene_name','description','gene_biotype', …
  • annotation - biomaRt - Error in UseMethod
    written by ferreira.cristina • ferreira.cristina
    Hi, I try run the getBM from biomaRt but a error was returned. ensembl=useMart("ENSEMBL_MART_ENSEMBL") ensembl = useDataset("hsapiens_gene_ensembl",mart=ensembl) list <- c("DUOX2", "DIRAS3") x <- getBM(attributes = …
  • How to extract all attribute from Biomart?
    written by Learner Learner
    I want to extract all attributes , normally I do the following to extract 2 attributes, in this example: "ensembl_gene_id" and "hgnc_symbol": library("biomaRt") ensembl <- …
  • GetBM gives controversial results
    written by lizaveta • lizaveta
    I want to get chromosome location of SMARCA1 gene. library("biomaRt") human = useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset = "hsapiens_gene_ensembl",host="grch37.ensembl.org", path="/biomart/martservice") getBM(attributes = c("hgnc_symbol","chromosome_name", "start_position", "end_position","strand"), filters = "hgnc_symbol", …
  • Biomart error while convertint transcript ID.
    written by always_learning always_learning
    library("biomaRt") ensembl<- useMart("ensembl", dataset="hsapiens_gene_ensembl") value =read.csv("Refseq.txt", strip.white=TRUE) colnames(value) = c("name") values <- as.vector(value$name) getBM(attributes=c("refseq_mrna", "ensembl_gene_id", "hgnc_symbol"), filters = "refseq_mrna", values = values, mart= ensembl) and …
  • Ensembl Biomart Plant database
    written by #### ####
    I want to fetch Arabidospsis annotation via ensembl biomart R package, here is the example how I am trying to fetch it: library("biomaRt") species = …
  • trying to retrieve data on go term from ensembl through biomaRt
    written by chrisclarkson100 chrisclarkson100
    <p>Trying to retrieve data on genes pertaining to the GO term mentioned.... it won&#39;t work...</p> <pre> ensembl = useMart(&quot;ensembl&quot;,dataset=&quot;hsapiens_gene_ensembl&quot;) go=c(&quot;GO:0006306&quot;) getBM(attributes=&quot;hgnc_symbol&quot;, &quot;hgnc_id&quot; filters=c(&quot;go&quot;, &quot;chromosome_name&quot;), values=list(go, …