SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Did I make directional or non-directional BS-seq libraries? shawpa Sample Prep / Library Generation 8 11-22-2017 03:44 AM
NEBNext Ultra-Directional RNA Library Prep Kit anna_m Sample Prep / Library Generation 2 10-18-2016 12:01 PM
Directional/non-directional libraries vs paired-end/single sequencing Snake75 RNA Sequencing 2 03-07-2015 11:46 AM
Directional information from PE reads? tboothby Illumina/Solexa 11 12-01-2011 06:39 PM
Tophat library-type about the Directional Ligation sequencing (DeLi-seq) strategy evergreen Bioinformatics 0 09-08-2011 07:12 PM

Reply
 
Thread Tools
Old 05-23-2018, 06:36 PM   #1
arctan
Junior Member
 
Location: United States

Join Date: Dec 2017
Posts: 8
Default Seeing reads in both directions from a directional library?

I made bacterial RNA-seq libraries with ScriptSeq v2 and did a small sequencing test (MiSeq Nano, Single Read, 50 bp).

ScriptSeq v2 prepares directional libraries, as per the manual: "The sequence produced by the Read 1 Sequencing Primer is that of the sense strand of the original fragmented RNA molecule."

After trimming and QC, I aligned the reads against the reference genome using BowTie.

When I visualize the resulting SAM files on a genome viewer (IGV), I see reads overlapping annotated CDS features going in _both_ directions. This is definitely _not_ what I expected, since the library should be directional.

Since I am new to this, I am hoping to understand might be going on. Is there something wrong with my library?

Is it possible that there was some genomic DNA contamination in my RNA preps (although I did treat with DNase and I didn't see a large peak on Bioanalyzer) which then became fragmented and incorporated as a sequence-able molecule in the library?

What are the possible reasons for seeing bidirectional reads from what was supposed to be a directional library?
arctan is offline   Reply With Quote
Old 05-24-2018, 04:40 AM   #2
kmcarr
Senior Member
 
Location: USA, Midwest

Join Date: May 2008
Posts: 1,148
Default

Not surprising for bacterial RNA-Seq. Bacteria will transcribe large tracks of their genome from both strands, in polycistronic primary transcripts which may overlap. The implication is that reads aligning to the reverse strand of an annotated CDS may represent the primary transcript of a nearby gene on the opposite strand downstream of the annotated gene you are looking at.
kmcarr is offline   Reply With Quote
Old 05-24-2018, 08:30 AM   #3
arctan
Junior Member
 
Location: United States

Join Date: Dec 2017
Posts: 8
Default

Thanks for the reply! I will keep that in mind and a quick literature search does show significant antisense transcription in bacteria.

How would I distinguish real antisense transcription from genomic DNA contamination in the RNA-seq reads?
arctan is offline   Reply With Quote
Old 05-29-2018, 02:15 PM   #4
asanair
Junior Member
 
Location: Australia

Join Date: May 2018
Posts: 1
Default

Hi
I am Asha, We are planning to do Miseq WGS on food samples and completely new to this field. Our aim is to identify the food content , species specific to look for meat adulteration etc. which will be the a suitable pipeline for this, it will be great if you can put some inputs into this.
thanks
Asha
asanair is offline   Reply With Quote
Old 06-01-2018, 07:31 AM   #5
arctan
Junior Member
 
Location: United States

Join Date: Dec 2017
Posts: 8
Default

Well, I ended up getting the results of my run and the directionality (strand bias) is so bad and I reported about it here: http://seqanswers.com/forums/showpos...10&postcount=5

Note to Asha (asanair): I'm pretty new to this stuff. It's probably a good idea to post your new question on a new thread with a suitable subject line, instead of posting it as a reply in this thread. It's probably also better to be as specific as possible about your questions. Good luck! I know I certainly need it
arctan is offline   Reply With Quote
Reply

Tags
alignments, directional rna-seq, miseq, scriptseq

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 01:09 AM.


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