Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat and cufflinks of non-strand-specific reads gfmgfm Bioinformatics 3 11-09-2012 04:35 AM
Mapping strand-specific data with TopHat Comosus RNA Sequencing 6 08-13-2012 01:25 AM
strand specific RNA-seq mapping qqtwee RNA Sequencing 1 07-10-2012 09:23 PM

Thread Tools
Old 01-25-2013, 01:09 AM   #1
Location: Canada

Join Date: Jun 2011
Posts: 47
Default mapping strand-specific RNA-seq reads using TopHat

I am using ScriptSeq Complete Gold Kit to generate pair-end RNAseq data, and wished to map the reads in a strand-specific manner. I referred TopHat FAQ, which suggested to try 'library-type' - 'fr-firststrand' and 'fr-secondstrand' separately, and the one resulting in many more splicing junctions in junctions.bed should be the right library type.

I checked the forum, and people said that the kit should go with 'fr-secondstrand', but I still wished to double-check it. So I ran 'unstranded','firststrand' and 'secondstrand' separately on my read pairs, but found the results are quite unexpectedly similar. Specifically the junctions.bed and acceptedhits.bam file are pretty much of the same size for all of them. As I am using pair end, presumably if I force strand-specific mapping, at least one option above should give really poor alignment. But the alignment was quite similar, as you can see from the attached figure (from 1st to 3rd panels are first-strand, second-strand and non-stranded, respectively).

Can you help see what's wrong here? Thanks!

tophat --library-type fr-secondstrand -o $3 -r $4 -p 4 -g 1 -G $GTF $Ref $read_A $read_B
Attached Images
File Type: png Screen Shot 2013-01-25 at 2.04.13 AM.png (40.4 KB, 31 views)
caswater is offline   Reply With Quote
Old 03-17-2015, 03:41 AM   #2
Location: USA

Join Date: Mar 2012
Posts: 38

Unstranded library of RNASeq can have strand bias, If yes then how can it be ruled out?
do we really look into it in case of fusion identification using tophat fusion?
bharati is offline   Reply With Quote
Old 03-17-2015, 08:27 AM   #3
Senior Member
Location: California

Join Date: Jul 2014
Posts: 198

Tophat by default aligns to the genome and retains discordant mappings as well. So that would explain why your three alignments are fairly similar. A better way to check the strand specificity and library-type is using RSeqC.

fanli is offline   Reply With Quote

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 02:26 PM.

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