This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding UTR regions in the RNA sequencing data from a non-model organism

Hi all friends,

I'm working on a RNA-seq project of a non-model plant, the library is generated from mRNA fraction (enriched with polyA) and sequenced as PE, stranded-specific. I have done de novo transcriptome assembly and annotation. Now, I would like to know if there is any way to determine the 3 and 5 UTR region of genes? Regarding 3 UTR, Since the library enrichment was done by oligo-dT primers, I'm not concerned about it, but I don't know what is the right procedure to determine these regions? Could you please advise me on this issue?

Thanks in advance

utr rna-seq non-model organism
• 5,004 views • link
updated by EVR EVR • written by seta seta

Hi!, I have to run a similar task...Did you find any solution?

Thanks!

• 0 views • link
written by etarisal • etarisal

2 answers

Hi Seta,

Use Transdecoder. It will output gff3 file which includes all the transcripts and their respective UTRs location and CDS location. Its very effective and reliable.

• 0 views • link
written by EVR EVR

You want to use an ORF finder to find the longest open reading frame in each transcript. EMBOSS has tools for this. One thing to be careful of is that if you have missed the 5' end of the transcript, you might not find a start codon for all your transcripts. The UTR is then most probably the sequence after the longest ORF.

• 0 views • link
written by i.sudbery i.sudbery

Log in to answer this question.

More posts like this