Question: How to retrieve "UTR" sequence from "Exon" of "Canonical Transcript" from a list of gene symbols?
0
gravatar for Hyun Jin
4.9 years ago by
Hyun Jin20
United States
Hyun Jin20 wrote:

Hi, 

I'm trying to retrive "UTR" sequence (i.e. 5'UTR sequence) from "Exon" of "Canonical Transcript" using a list of gene symbols. 

 

First I tried Ensembl's martview to do the job (http://www.ensembl.org/biomart/martview/ ).

But output provides more 5UTR sequences (i.e 150 results) then than input gene numbers (i.e. 80)

I assume that one gene could correspond to multiple transcripts and this is why it gives more results than It was given. Note that my input was gene ID, not transcript ID. 

 

Therefore, I thought retriving the UTR exon sequence only from "Canonical transcript" using a list of gene symbol could avoid this problems.

Then I searched a bit, and tried the table browser from UCSC genome browser, selected KnownCanonical table, which I can select only one canonical transcripts, but it only provides coordinates of UTR, not the sequence itself. 

 

Please advice me or let me know any reference that I can look up, that would be very helpful.

I can do some basic R programming, but never used BiomaRt. but I'm willing to try if BiomaRt 'getSequence' is the way to go. 

Thank you very much!!!

 

genome browser biomart R • 1.8k views
ADD COMMENTlink modified 4.8 years ago by Devon Ryan92k • written 4.9 years ago by Hyun Jin20
1
gravatar for Devon Ryan
4.8 years ago by
Devon Ryan92k
Freiburg, Germany
Devon Ryan92k wrote:

If you have coordinates (ideally in BED format) then you can just use bedtools getfasta to retrieve the corresponding sequence.
 

ADD COMMENTlink written 4.8 years ago by Devon Ryan92k

Thank you Devon, I do have BED file with canonical transcripts 5UTR coordinates. I will try

ADD REPLYlink written 4.8 years ago by Hyun Jin20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1267 users visited in the last hour