Using multiple filters in BiomaRt's getSequence
1
0
Entering edit mode
7.2 years ago
spiral01 ▴ 110

Hi. I have a data frame of ensembl gene IDs and their corresponding canonical transcript IDs. Below is a small subset of this data frame:

             geneID    transcriptID
ENSG00000186092 ENST00000335137
ENSG00000187634 ENST00000420190
ENSG00000187642 ENST00000341290
ENSG00000188157 ENST00000379370
ENSG00000186891 ENST00000379265

I would like to retrieve the protein coding sequences for these genes using the R package biomaRt. Using the code below I can do this:

mart <- useMart("ensembl", dataset="hsapiens_gene_ensembl")
seq = getSequence(id =  "ENSG00000186092",
                  type = c("ensembl_gene_id"),
                  seqType = "coding",
                  mart = mart)

My issue is that I would like to add a second filter of transcript ID so that I am only getting the gene sequence for the specific transcript. However, I am not sure how to do this for the getSequence function. I have tried to add a value argument (as with getBM) but this gives the unused argument error message. Does anyone know how to do this?

r biomart ensembl • 3.3k views
ADD COMMENT
1
Entering edit mode
7.2 years ago
Emily 24k

Change your filter so that you're filtering by the Ensembl transcript ID:

seq = getSequence(id =  "ENST00000335137",
               type = c("ensembl_transcript_id"),
               seqType = "coding",
               mart = mart)
ADD COMMENT
0
Entering edit mode

Hi Emily, thanks that does help. However is it actually possible to use two different filters?

ADD REPLY
0
Entering edit mode

You can filter by multiple types of ID, but you can filter by a list of IDs, plus a biotype or something.

ADD REPLY
0
Entering edit mode

Hi Emily, did you mean "You can <<not>> filter by multiple types of ID, but ..." ?

ADD REPLY

Login before adding your answer.

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