SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Strand specific RNA-seq data from dUTP protocol ZhengXia Bioinformatics 8 10-16-2014 04:36 AM
TopHat: Aligning SOLiD strand specific RNAseq Ender985 RNA Sequencing 3 07-09-2012 10:37 PM
Strand-specific library appears not strand-specific oligo Illumina/Solexa 7 12-08-2011 09:54 AM
Cufflinks library-type for SOLiD strand specific single data gengen Bioinformatics 1 10-19-2010 08:41 AM
What is this strand specific pair-end RNA-seq data? xhuister Bioinformatics 2 09-14-2010 06:34 AM

Reply
 
Thread Tools
Old 02-17-2012, 04:51 AM   #1
Comosus
Junior Member
 
Location: UK

Join Date: Mar 2011
Posts: 2
Default Mapping strand-specific data with TopHat

I have RNA-seq data that I want to map using TopHat. The library prep was done using ScriptSeq v2, and the sequencing using Illumina.

The R1 sequences the sense cDNA from start to finish. The R2 sequences the antisense cDNA from finish to start.

The --library-type option fr-unstranded says that "Reads from the left-most end of the fragment (in transcript coordinates) map to the transcript strand, and the right-most end maps to the opposite strand.", which would seem to be just what I need.

However how does it deal with sense cDNA that has been expressed from the negative strand of the DNA? If an R1 read is just that, then wouldn't the read have to be reversed to map to the genome?

Has anyone dealt with anything similar, and if so, do you know which settings I should use for TopHat?
Comosus is offline   Reply With Quote
Old 03-13-2012, 07:43 AM   #2
Oliviervg
Member
 
Location: Brussels

Join Date: Sep 2011
Posts: 13
Default

I have the same problem.
I prepared my library with Scriptseq V2 (directional).

The sequencing was done with Illumina PE 2*76bp.


Can anyone explain more precisely the --library-type option of TopHat ?
I don't understand the difference between "fr-unstranded" and "fr-secondstrand"


What's the best parameters in my case ?

Thank-you.

Last edited by Oliviervg; 03-13-2012 at 08:00 AM.
Oliviervg is offline   Reply With Quote
Old 03-15-2012, 12:46 AM   #3
Oliviervg
Member
 
Location: Brussels

Join Date: Sep 2011
Posts: 13
Default

Any help would be very appreciated !
Thank you !
Oliviervg is offline   Reply With Quote
Old 03-15-2012, 04:18 AM   #4
Thomas Doktor
Senior Member
 
Location: University of Southern Denmark (SDU), Denmark

Join Date: Apr 2009
Posts: 105
Default

TopHat doesn't have those options, are you sure you're not thinking of Cufflinks?
For the scriptseq v2 kit I would use the fr-secondstrand librarytype since the first read will be the left-most in transcript coordinates.
Thomas Doktor is offline   Reply With Quote
Old 03-15-2012, 04:41 AM   #5
Oliviervg
Member
 
Location: Brussels

Join Date: Sep 2011
Posts: 13
Default

Thank you for your answer, I m going to use fr-secondstrand.
And, my TopHat 1.4 has those options !
Oliviervg is offline   Reply With Quote
Old 03-15-2012, 07:20 AM   #6
Thomas Doktor
Senior Member
 
Location: University of Southern Denmark (SDU), Denmark

Join Date: Apr 2009
Posts: 105
Default

You're right, they were just well hidden (and I never had use for anything other than the default).
Thomas Doktor is offline   Reply With Quote
Old 08-13-2012, 01:25 AM   #7
harryzs
Member
 
Location: Germany

Join Date: Dec 2010
Posts: 29
Default

Just a reminder, if you did not see this:

http://tophat.cbcb.umd.edu/faq.html
harryzs is offline   Reply With Quote
Reply

Tags
illumina, rna-seq, strand-specific, tophat

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 03:18 PM.


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