SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
paired-end reads mapped to genome.. gene with only one direction of paired-end reads? danwiththeplan Bioinformatics 2 09-22-2011 02:06 AM
How to count Paired-End reads for RNA-Seq read counts epistatic RNA Sequencing 0 08-31-2011 02:44 PM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-04-2011 11:51 PM
RNA-Seq: The Drosophila melanogaster transcriptome by paired-end RNA sequencing. Newsbot! Literature Watch 0 12-24-2010 02:13 AM
RNA-Seq: A Probabilistic Framework for Aligning Paired-end RNA-seq Data. Newsbot! Literature Watch 11 10-16-2010 09:27 AM

Reply
 
Thread Tools
Old 11-19-2008, 08:32 PM   #1
xuying
Member
 
Location: Shanghai, China

Join Date: Mar 2008
Posts: 16
Default RNA-seq with paired-end reads

anyone's doing RNA-seq with paired-end reads now? how to use pairing information is a challenge. any good ideas?
xuying is offline   Reply With Quote
Old 11-19-2008, 09:55 PM   #2
dvh
Member
 
Location: london, uk

Join Date: Jul 2008
Posts: 35
Default

Yes. We've done paired end RNAseq - partly because of mixed sample type (some DNAseq etc) flow cells.

The paired end info is useful for novel gene expression and splice junction discovery.

For simple RPKM etc we just treated each read as single end and aligned to genome+splice junctions seperately - this is likely to be suboptimal although easy approach. I'd be very interested in people's comments...

david
dvh is offline   Reply With Quote
Old 12-08-2008, 10:03 PM   #3
sci_guy
Member
 
Location: Sydney, Australia

Join Date: Jan 2008
Posts: 83
Default Rna-pet

Paired-end ditags (PET) is an interesting strategy. It is very economical way of finding splice junction variants. Of course it doesn't offer the same resolution as full paired end RNA-Seq. I believe the team at GIS have migrated to Illumina and SOLiD in particular for RNA-PET now.
sci_guy is offline   Reply With Quote
Old 01-26-2009, 10:31 AM   #4
monad
Member
 
Location: California

Join Date: May 2008
Posts: 40
Default

PE and long read would be better to detect splicing junction and SNPs.
monad is offline   Reply With Quote
Old 06-04-2009, 01:41 AM   #5
ieuanclay
Member
 
Location: Basel, Switzerland

Join Date: Feb 2009
Posts: 27
Default

Hi all, We have had success with PE RNAseq (illumina).

Some of the main pitfalls that we found:
-when you make your cDNA, if you use random hexamers (best results for us), your first strand synthesis reaction will tend to produce sequences that start/end in GC rich sequences. We have found that these tend to be non-accurate annealling products, and you can improve your mapping efficiency by stripping the first/last few bases.
- TopHAT (see bowtie/ben langmeads posts) has also been very useful (and extremely fast for mapping!!
- Creating 'spliced' reference genomes is very useful for retaining the benefits of PE mapping, without having to introduce huge 'insert' sizes...
- Avoid PET tags - the extra ligation steps, if you are not careful, can introduce a high percentage of 'false' splice events... With the illumina PE strategy there is no need for this...

Has anyone had any luck with barcoding RNAseq samples?

Ieuan
ieuanclay is offline   Reply With Quote
Old 07-07-2009, 12:12 AM   #6
seq_GA
Senior Member
 
Location: Asiana

Join Date: Feb 2009
Posts: 124
Default

Quote:
Originally Posted by ieuanclay View Post
Hi all, We have had success with PE RNAseq (illumina).

Some of the main pitfalls that we found:
-when you make your cDNA, if you use random hexamers (best results for us), your first strand synthesis reaction will tend to produce sequences that start/end in GC rich sequences. We have found that these tend to be non-accurate annealling products, and you can improve your mapping efficiency by stripping the first/last few bases.
- TopHAT (see bowtie/ben langmeads posts) has also been very useful (and extremely fast for mapping!!
- Creating 'spliced' reference genomes is very useful for retaining the benefits of PE mapping, without having to introduce huge 'insert' sizes...
- Avoid PET tags - the extra ligation steps, if you are not careful, can introduce a high percentage of 'false' splice events... With the illumina PE strategy there is no need for this...

Has anyone had any luck with barcoding RNAseq samples?

Ieuan

Hi Ieuan,

It is interesting to see PE RNAseq.
I still ahve some doubts regarding mapping PE reads. Because ELAND can support eland_rna analysis for SR(single read) only.
Any tool to map tags with splice juction or can you please explain about the mapping part of PE RNA seq tags?

Thanks.
seq_GA is offline   Reply With Quote
Old 07-07-2009, 02:22 AM   #7
ieuanclay
Member
 
Location: Basel, Switzerland

Join Date: Feb 2009
Posts: 27
Default

Check the forums here for TopHat if you want to know about mapping RNAseq!

The actual way that you map is very dependent on the source of your RNA, so i would need to know more about your experiments! My main experience is with nascent RNA, so we PE map againstt a genomic reference. Spliced RNA can be mapped against spliced genomes (again see TopHat). Or a combination of both can be used, i.e. mapping against the genomic reference using a sensible PE insert limit, to find nascent RNA, then mapping the reads that don't map to the reference genome against a spliced genome. This way you can flag potentially spliced and unspliced transcripts.

There really are too many ways of proceeding, you just have to pick a processing pipeline that makes sense for your material/experiment and that you can justify when it comes to publishing.

Hope this helps!?
ieuanclay is offline   Reply With Quote
Old 08-22-2009, 08:05 PM   #8
baohua100
Senior Member
 
Location: Canada

Join Date: Jun 2008
Posts: 103
Default

no one show how to use paired-end info to infer junction and alternative splicing?
baohua100 is offline   Reply With Quote
Old 09-05-2009, 12:31 PM   #9
karthikprabu
Junior Member
 
Location: italy

Join Date: Sep 2009
Posts: 6
Default

hai ....I would like to know is it possible to use Bowtie to map the reads which does not have its genome sequenced?
karthikprabu is offline   Reply With Quote
Old 09-11-2009, 03:00 PM   #10
Pepe
Member
 
Location: Germany

Join Date: Mar 2009
Posts: 28
Default

You can't use Bowtie without a reference genome. But there are others that can (Velvet, etc.)... just search for them in this forum...
Pepe 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 03:14 PM.


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