SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
sharp peak in library - TruSeq Stranded Total RNA (rRNA-depleted E. coli) Seqasaurus Sample Prep / Library Generation 1 11-05-2015 02:02 AM
Majority of reads counted in stranded=reverse using HTSeq? sdlmark RNA Sequencing 2 03-17-2015 09:55 AM
Remove reverse complement redundancy in stranded transcriptome evt8 RNA Sequencing 4 11-05-2014 04:28 PM
End repair of TruSeq RNA and TruSeq mRNA stranded ykgenome Sample Prep / Library Generation 2 05-15-2014 12:01 AM
DNA - Assembly: Algorithmic recognition of reverse complement reads BerGer23 Bioinformatics 0 11-12-2012 02:47 AM

Reply
 
Thread Tools
Old 02-24-2016, 03:07 AM   #1
lordjonwald
Junior Member
 
Location: Belfast

Join Date: Oct 2014
Posts: 5
Question Unmerged reads from truseq stranded RNA library, reverse complement or not?

Hi all,

I'm trying to process metatranscriptomic reads (2x75) from a nextseq500 run, the libraries were prepared using a Truseq stranded RNA kit.

It appears that many of my paired ends do not merge, which is fine, but I'd like to concatenate the unmerged and merged reads to align against some protein databases and for de novo assembly and annotation if possible,

so my question is: is it necessary to reverse complement the second paired ends from the unmerged reads before concatenating? I imagine this wouldn’t be necessary for assemblers which like paired end reads as input but what if I wanted to input the reads as one single end read file?

Additionally (if I'm reading the truseq stranded protocol correctly) it seems that it is the antisense strand which is sequenced, if that's the case, is it actually the first paired end file which would need to be reverse complemented?

I've never worked with transcriptomic data before and I'm not entirely sure what the best practices are yet, so apologies if these are silly questions,

Anyway, any help at all would be greatly appreciated

Many thanks
lordjonwald is offline   Reply With Quote
Old 02-24-2016, 11:01 PM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Aligners and assemblers don't care which strand a read comes from; you don't need to reverse-complement anything in this case. If you were mapping to a genome to discover which strand transcripts came from then it could be useful, but not otherwise.

However, you'll get more accurate mapping results by running the merged reads as single-ended and the unmerged reads as pairs in a different process, then combining the outputs.
Brian Bushnell is offline   Reply With Quote
Old 02-25-2016, 05:40 AM   #3
lordjonwald
Junior Member
 
Location: Belfast

Join Date: Oct 2014
Posts: 5
Default

Thanks for your reply Brian, that really helps clear things up, I'll bear that in mind if I ever get around to mapping to a genome!
lordjonwald is offline   Reply With Quote
Reply

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 06:32 PM.


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