SEQanswers

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



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: Single Read and Paired End mRNA-Seq Illumina Libraries from 10 Nanograms Tot Newsbot! Literature Watch 0 11-09-2011 03:10 AM
How to count Paired-End reads for RNA-Seq read counts epistatic RNA Sequencing 0 08-31-2011 03:44 PM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-05-2011 12:51 AM
searching 454 (single read) + Illumina paired-end or matepair data boetsie General 0 10-04-2010 02:45 AM
Difference in paired-end and single-end read ? darshan Bioinformatics 1 10-01-2009 12:44 AM

Reply
 
Thread Tools
Old 08-03-2011, 04:18 PM   #1
fbarreto
Member
 
Location: San Diego

Join Date: Jan 2010
Posts: 15
Default RNA-seq: read counts in single- vs paired-end sequences

Hi,

I've read several threads with very related topics, but none have addressed the one simple question I have (which stems from my being new to this).

We already have an assembled transcriptome (from 454 data), and would now like to do an experiment to assess differential gene expression with Illumina.

My question is (assume one Illumina lane is being sequenced): will a 2 x 100 bp paired-end lane generate approximately HALF the number of reads of a 1 x 100 bp single end lane? (By read I mean a copy of a transcript, i.e. the two ends of a paired-end sequence count as only ONE read, for my purposes)

I would assume so, since the sequencing effort (in terms of base pairs) for such a paired-end is twice that of the single end.

Since I'm interested only in maximizing read counts, the extra information that comes in the paired-end sequence is not of interest to me at this point.

Any insight regarding the comparison of read counts per sequencing effort between the two methods would be greatly appreciated.

Thanks!
fbarreto is offline   Reply With Quote
Old 08-03-2011, 04:43 PM   #2
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

You produce the same number of reads. In a paired end run, you sequence BOTH ends of the same fragment, so both reads are coming from the same cluster on the flow cell. So with a paired-end run, you are getting exactly the same number, so this has no effect on your total read count (if you count the paired end reads as one read).

The advantage of a paired end run is that you can align the reads with greater confidence.
chadn737 is offline   Reply With Quote
Old 08-03-2011, 06:00 PM   #3
fbarreto
Member
 
Location: San Diego

Join Date: Jan 2010
Posts: 15
Default

Thanks for insight! But wouldn't a a 2 x 100bp paired-end run be sequencing a total of 200bp from each sequence (i.e. transcript copy), 100bp from each end?

So if a run sequences 1,000,000 bp (I know it's more, but just to keep it simple), I would get 5000 paired-end and 10,000 single end reads. Am I thinking about this wrong?

Thanks again!
fbarreto is offline   Reply With Quote
Old 08-03-2011, 06:09 PM   #4
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

I think you are confusing base pairs sequenced with reads count. The number of reads you get is a factor of the number of clusters on a flow cell. I think the recommended cluster density on an Illumina Hi-seq typically gives you ~65 millions per lane (somebody correct me if I am wrong). From these 65 million reads, doing a paired end sequencing run vs a single end run means that you will get 65 million pairs (or 130 million reads) vs 65 million single end reads. However, since you sequenced 200 bps versus 100 bps, you will get double the number of base pairs sequenced.

If all you are doing is differential gene expression, then you are not really concerned with base pairs sequenced, just number of reads. In this, paired end versus single end doesn't affect the number of reads.
chadn737 is offline   Reply With Quote
Old 08-03-2011, 06:19 PM   #5
fbarreto
Member
 
Location: San Diego

Join Date: Jan 2010
Posts: 15
Default

That's great to know! Thanks a lot!
My confusion was that, in my simple mind, the sequencing effort (in bp) versus number of reads had to balance out such that getting longer sequences would mean getting fewer. But if not, like you said, then paired-end seems like the way to go.

Much appreciated!
fbarreto is offline   Reply With Quote
Reply

Tags
paired-end, read count, rna-seq, single end

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 11:28 AM.


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