[BioC] why biomaRt cannot extract 3UTR sequences for 1941 ENSGxxxxx ?

michael watson (IAH-C) michael.watson at bbsrc.ac.uk
Fri May 28 22:54:52 CEST 2010


The following (small) code works for me:

library(biomaRt)
mart <- useMart("ensembl","hsapiens_gene_ensembl")
ids <- c("ENST00000262187","ENST00000296271")
seq <- getSequence(id=ids, type="ensembl_transcript_id", mart=mart, seqType="3utr")
seq
________________________________________
From: bioconductor-bounces at stat.math.ethz.ch [bioconductor-bounces at stat.math.ethz.ch] On Behalf Of mauede at alice.it [mauede at alice.it]
Sent: 28 May 2010 21:41
To: Bioconductor  List
Subject: [BioC] why biomaRt cannot extract 3UTR sequences for 1941 ENSGxxxxx    ?

I executed the following lines several times from a script as well as pasting them in an R shell.
Systematically biomaRt is failing.
The problem is to extract the 3UTR sequences corresponding to a vector containing 1941
Ensembl Transcript numbers (some are duplicated ... is this s problem ?)
Please, find the failing instructions in the following including the ENST vector

Any suggestion is welcome. Thank you,
Maura

> hmart <- useMart('ensembl', dataset='hsapiens_gene_ensembl')
Checking attributes ... ok
Checking filters ... ok

> genes_map[,"ensembl_transcript_id"]
   [1] "ENST00000262187" "ENST00000296271" "ENST00000346166" "ENST00000381570"

<snip>

[1937] "ENST00000400907" "ENST00000400908" "ENST00000440864" "ENST00000309042"
[1941] "ENST00000254325"

> genes_seq <- getSequence (id=genes_map[,"ensembl_transcript_id"],type="ensembl_transcript_id",
+ seqType="3utr",mart=hmart)
Error in value[[3L]](cond) :
  Request to BioMart web service failed. Verify if you are still connected to the internet.  Alternatively the BioMart web service is temporarily down.



tutti i telefonini TIM!


        [[alternative HTML version deleted]]

_______________________________________________
Bioconductor mailing list
Bioconductor at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list