SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
paired-end RNA-Seq of short transcripts fbarreto Illumina/Solexa 3 11-23-2011 11:43 AM
RNA-Seq: Identification of novel transcripts in annotated genomes using RNA-Seq. Newsbot! Literature Watch 0 06-29-2011 01:40 PM
RNA-Seq: FusionHunter: identifying fusion transcripts in cancer using paired-end RNA- Newsbot! Literature Watch 1 05-19-2011 12:56 AM
RNA-Seq: Composite Transcriptome Assembly of RNA-seq data in a Sheep Model for Delaye Newsbot! Literature Watch 0 03-26-2011 02:02 AM
RNA SEQ Bias toward short transcripts (RPKM) vruotti Bioinformatics 4 01-06-2010 10:09 AM

Reply
 
Thread Tools
Old 02-09-2012, 07:07 AM   #1
StopCodon
Junior Member
 
Location: India

Join Date: Feb 2012
Posts: 7
Default Transcripts from RNA-seq assembly

I assembled the transcriptome of an unsequenced organism using Trinity. I would like to do some downstream analysis on the assembled transcripts. When I look at the transcripts and translate them, I find they are littered with stop codons, does that indicate there is some issue with sequencing? How can one go from the transcript to the amino acid sequence? Can we assume that the beginning of the transcript would be the beginning of the transcription start site for the gene? I would love any pointers on this.

Thanks,
StopCodon is offline   Reply With Quote
Old 02-09-2012, 08:28 AM   #2
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 817
Default

Quote:
Can we assume that the beginning of the transcript would be the beginning of the transcription start site for the gene?
This is unlikely. Trinity assembles the whole RNA, including untranslated regions. If you want to look at protein predictions de-novo, then you need something like frameDP to correct for frame shift errors and identify the most likely start/stop sites for translation.
gringer is offline   Reply With Quote
Old 02-09-2012, 11:35 AM   #3
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 991
Default

Are you sure you were in the right reading frame when checking for stop codons? Look for start codons to fix the reading frame.

Quote:
Can we assume that the beginning of the transcript would be the beginning of the transcription start site for the gene?
I would say, possibly yes, if you meant to say transcription start site. But gringer might have been right in assuming that you meant the translation start.
Simon Anders is offline   Reply With Quote
Old 02-09-2012, 02:59 PM   #4
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 817
Default

Quote:
Originally Posted by Simon Anders View Post
I would say, possibly yes, if you meant to say transcription start site. But gringer might have been right in assuming that you meant the translation start.
Whoops, yes, I read that wrong. Just to add to this, depending on the mapping quality the start point may not necessarily be exactly at the transcript start location.
gringer is offline   Reply With Quote
Old 02-10-2012, 02:48 AM   #5
StopCodon
Junior Member
 
Location: India

Join Date: Feb 2012
Posts: 7
Default

Thanks. That clarifies some concepts. Yes, I had meant the Translation start site instead of transcription. Now my data makes more sense as well since very few of the assembled transcripts were starting with a start codon. Is FrameDP a good choice for finding peptide sequences? I saw it had only 10 citations to date. Are there any other tools for doing this?

Thanks.
StopCodon is offline   Reply With Quote
Old 02-10-2012, 02:52 AM   #6
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 817
Default

Quote:
Is FrameDP a good choice for finding peptide sequences? I saw it had only 10 citations to date. Are there any other tools for doing this?
Have a look here for other options:
http://seqanswers.com/forums/showthread.php?t=4036

FrameDP generates huge amounts of files, something like 5-10 times the number of transcripts, so you need to be careful with that. My computer claimed to have run out of space (with ~300GB free) due to me not being careful enough running that program.
gringer is offline   Reply With Quote
Old 07-08-2015, 02:58 AM   #7
GillermoPonz
Junior Member
 
Location: Madrid

Join Date: Oct 2014
Posts: 6
Default

I'm having the same problem just now. When I translate the assembled transcripts (in the 6 pssible reading frames), they are full of stop codons and I don't know what to do. I've been thinking about not normalizing the libraries when assembling in Trinity but I don't know if it make sense...

Did you finally solve the problem?

Thank you
GillermoPonz is offline   Reply With Quote
Reply

Tags
denovo assembly, rna-seq, transcriptome, trinity

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:39 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO