SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing
Similar Threads
Thread Thread Starter Forum Replies Last Post
vcfutils.pl -specifying minimum read depth? Lspoor Bioinformatics 3 05-27-2013 01:20 AM
Maximum read Depth in Samtools Anjali Bioinformatics 2 01-16-2012 02:15 AM
What is read depth -D100 means in samtools? ketan_bnf Bioinformatics 1 07-28-2011 06:54 PM
Read distribution at high sequence depth ForeignMan General 10 05-26-2011 03:50 AM
About the read depth of coverage El Mariachi Illumina/Solexa 2 12-30-2010 12:22 AM

Reply
 
Thread Tools
Old 09-27-2011, 12:44 AM   #1
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default Read depth recommendations

Hi all,

We'd like to perform some RNA-seq to look at gene expression level changes in mouse hippocampus due to a treatment of interest to us. We're not interested in finding new transcripts or looking for differences in splice junctions or anything of that sort. Consequently, I'm curious what people are recommending these days in terms of read depth.

On a related note, I've read a number of people here suggesting that paired end reads are probably not required for our sort of project. If that's the case, I'm curious what sort of read lengths (36bp, 50bp, etc.) people having been using that give them meaningful results.

Any suggestions you might have would be appreciated.
dpryan is offline   Reply With Quote
Old 09-30-2011, 10:30 AM   #2
mbblack
Senior Member
 
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245
Default

I was at a meeting a couple of weeks ago (the 2011 TIES meeting at UNC) and in a talk by Wendell Jones (a statistician with the company Expression Analysis) he talked briefly about this.

An Illumina white paper from a few years ago argued that 2-10 million mapped reads should be in the range of equal or better sensitivity than microarrays for differential expression estimation. Wendell, however, mentioned that his experience with experimental data over the years has seen that number climb, to where most of his clients are more often using 20-50 million mapped reads in order to be "comparable" or better to array data.

I think though, that most of these kinds of estimates are based on human data. We work mostly on rat and mouse models, and I honestly am not convinced of just what we need in terms of RNAseq coverage to get results equal to or better than our array results. For our first direct comparison, I have greater than 60 million mapped reads per sample (3 controls, 3 treatment animals, all mouse livers), but I get much less sensitivity for gene expression than with microarray data (same samples used too). We're trying another direct comparison soon (mouse liver samples already run with affy titan arrays) soon to be run on an ABI SoLid 5500xl, shooting for 10-20 million reads per sample.

Wendell also mentioned in his talk how differential expression significance has occasionally been seen to appear to be fine at low coverage, but suddenly drops out at high coverage, but he did not offer an explanation for that observation nor elaborate on the specificis.

Thus far in our research, we've been using 50bp single end reads, but I don't really think that 36bp reads would be a problem.

P.S. There is an FDA-led initiative called SEQC underway (a followup to the MAQC initiative - http://www.fda.gov/ScienceResearch/B...ls/default.htm ) - http://www.genomeweb.com/sequencing/...rna-sequencing which is intended to put some real numbers to issues like this, based on real comparison data.

<edit> actually, SEQC is also really MAQCIII, the third phase of the whole MAQC long term initiative. Some of the sequencing is done, some in the works right now and still some more to be done in the next few months. Data analysis is really just in the very initial stage.
__________________
Michael Black, Ph.D.
ScitoVation LLC. RTP, N.C.

Last edited by mbblack; 09-30-2011 at 10:36 AM.
mbblack is offline   Reply With Quote
Old 09-30-2011, 11:15 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Thanks mbblack, that's extremely helpful! I'll have to look more into SEQC and MAQC, they sound interesting.
dpryan 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 11:23 AM.


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