SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Illumina Directional RNA-Seq flobpf Illumina/Solexa 58 11-19-2015 04:18 PM
mRNA-seq data analysis smsm Illumina/Solexa 6 10-24-2012 04:14 AM
directional mRNA-Seq Sample Preparation Claire1 Sample Prep / Library Generation 3 02-02-2012 01:47 AM
Sequence for directional mRNA seq adaptors? gaster Illumina/Solexa 7 10-04-2010 06:53 AM
mrna seq or directional mrna seq link1 Sample Prep / Library Generation 0 08-12-2010 06:58 PM

Reply
 
Thread Tools
Old 12-03-2010, 04:12 PM   #1
jmtepp
Junior Member
 
Location: Oakland, CA

Join Date: Jul 2010
Posts: 1
Post Analysis of Directional mRNA-seq data / Illumina

We are developing a library preparation for directional mRNA-seq. All sequences will read 5' -> 3' direction. Our protocol isolates mRNA, fragments the messages and then these are re-purified with oligodT - resulting in a library of short 3' fragments for sequencing. This eliminates the need to adjust for transcript length when doing quantitative expression analysis. We're using bowtie to align the raw reads to a (cDNA) reference transcriptome (minus introns and intergenic regions). BUT - approximately 12% of the loci in this organism (Arabidopsis) have overlapping 3' tails (Watson/Crick strands) leading to ambiguities when using BowTie ( which looks for alignments on both strands). Basically we're looking for a short read alignment algorithm that would take our directional data and align it ONLY with transcripts in the 5' -> 3' direction. Illumina is developing a directional mRNA-seq library prep method similar to the one we are using (without the purification of the 3' ends). I wonder if anyone has experience analyzing directional mRNA-seq data?
jmtepp is offline   Reply With Quote
Old 12-03-2010, 04:50 PM   #2
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 504
Default

Since your libraries are directional, you can use that information to resolve the ambiguity of overlapping 3' ends. The tag will be in the same 5'->3' direction as one of the two genes (conventionally indicated as +).

Harold
HESmith is offline   Reply With Quote
Old 12-04-2010, 05:10 AM   #3
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

Quote:
Originally Posted by jmtepp View Post
Basically we're looking for a short read alignment algorithm that would take our directional data and align it ONLY with transcripts in the 5' -> 3' direction.
As you are already using Bowtie, have a look at its manual again, especially at the description of the option '--norc'. It causes Bowtie to do precisely this.

Last edited by Simon Anders; 12-04-2010 at 05:10 AM. Reason: sp
Simon Anders is offline   Reply With Quote
Old 05-09-2012, 05:53 AM   #4
qqtwee
Member
 
Location: Beijing

Join Date: Feb 2011
Posts: 16
Default strand specific RNA-seq mapping

Does anyone know how do mapping for strand specific RNA-seq data? which tools can work well? what the parameter of tools for strand specific infomation of RNA-seq data?
Looking forward to your answer.
Thank you very much! Best wishes!
qqtwee is offline   Reply With Quote
Old 05-09-2012, 03:19 PM   #5
jujubix
Member
 
Location: Vancouver

Join Date: May 2011
Posts: 14
Default

Quote:
Originally Posted by qqtwee View Post
Does anyone know how do mapping for strand specific RNA-seq data? which tools can work well? what the parameter of tools for strand specific infomation of RNA-seq data?
Looking forward to your answer.
Thank you very much! Best wishes!
Maybe TopHat with the '--library-type' option which will treat the reads in a strand specific manner as described in the manual?
jujubix is offline   Reply With Quote
Old 05-09-2012, 07:08 PM   #6
qqtwee
Member
 
Location: Beijing

Join Date: Feb 2011
Posts: 16
Default

Quote:
Originally Posted by jujubix View Post
Maybe TopHat with the '--library-type' option which will treat the reads in a strand specific manner as described in the manual?
Thank you for your reply! Do you know BWA has option that can treat reads in a strand specific ?
qqtwee is offline   Reply With Quote
Old 05-09-2012, 07:14 PM   #7
jujubix
Member
 
Location: Vancouver

Join Date: May 2011
Posts: 14
Default

Quote:
Originally Posted by qqtwee View Post
Thank you for your reply! Do you know BWA has option that can treat reads in a strand specific ?
No, I do not know. Although the manual probably says no?
jujubix is offline   Reply With Quote
Old 02-27-2014, 12:13 AM   #8
stanikay
Junior Member
 
Location: Pretoria

Join Date: Oct 2013
Posts: 4
Default Mappind strand specific RNAseq bacteria pair-end reads

Just following up on this thread, which mapper will be most ideal for mapping strand specific RNAseq bacteria pair-end reads,

I know tophat has the library type option but its mostly ideal for eukaryotic genomes.

Any thoughts will be appreciated.

Thanx
stanikay is offline   Reply With Quote
Reply

Tags
alignment, bowtie, directional rna-seq

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 05:17 AM.


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