[BioC] intron sequences from biomaRt
Durinck, Steffen (NIH/NCI) [F]
durincks at mail.nih.gov
Wed Jun 20 15:29:04 CEST 2007
If the BioMart web server would implement this, it would become readily available in R, but the are currently no plans for them to do this and neither do we currently intend to implement this ourselves in biomaRt (unless you want to contribute this function;) ). Would my previous suggestion of subtracting the exonic sequences from the unspliced transcript sequences not work? If not you'll indeed have to use the Ensembl Core Perl API and create a script to retrieve the intronic sequences.
From: Dario Greco [mailto:dario.greco at helsinki.fi]
Sent: Mon 6/18/2007 7:14 AM
To: Durinck, Steffen (NIH/NCI) [F]
Cc: bioconductor at stat.math.ethz.ch
Subject: Re: [BioC] intron sequences from biomaRt
thank you very much for your email. i asked to the ensembl helpdesk and
this is what they replied:
"...The only way to retrieve intronic sequences in batch mode is by using
the Ensembl Core Perl API:
I hope this answers your question. Please let us know if you have any
other questions or problems..."
is there any idea of implementing/using this from R? shall i go by scratch with bioperl?
Durinck, Steffen (NIH/NCI) [F] wrote:
> Dear Dario,
> No there is currently no possibility to select intronic sequences directly. You could request this at helpdesk at ensembl.org and see if they want to add this feature in future versions of Ensembl. For now the only way to do it with biomaRt would be to retrieve the transcript or gene sequences and the exon sequences and then find the intronic sequences by splitting the transcript sequences on the different exons.
> Hope this helps,
> -----Original Message-----
> From: Dario Greco [mailto:dario.greco at helsinki.fi]
> Sent: Fri 6/15/2007 9:13 AM
> To: bioconductor at stat.math.ethz.ch
> Subject: [BioC] intron sequences from biomaRt
> dear list,
> i need to get the intron sequences for a group of entrez gene ids.
> is there any way to do it using biomaRt? apparently there is no option in the
> getSequence() function.
> R version 2.5.0 (2007-04-23)
> biomaRt RCurl XML
> "1.10.0" "0.8-1" "1.7-3"
> any suggestions?
> thanks for your help.
Institute of Biotechnology - University of Helsinki
Building Cultivator II
P.O.Box 56 Viikinkaari 4
Office: +358 9 191 58951
Fax: +358 9 191 58952
Mobile: +358 44 023 5780
More information about the Bioconductor