SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-seq: Strategy for filtering low count genes? LeonDK RNA Sequencing 1 11-04-2014 12:22 AM
Different normalization methods with count data greener Bioinformatics 2 09-01-2013 06:27 AM
ChIP-Seq: A normalization strategy for comparing tag count data. Newsbot! Literature Watch 0 04-06-2012 06:00 AM
Quantile normalization for RNA seq data? Boel Bioinformatics 3 03-26-2010 04:07 PM

Reply
 
Thread Tools
Old 01-09-2015, 04:11 PM   #1
qliu
Junior Member
 
Location: Tucson,AZ

Join Date: Jan 2015
Posts: 1
Default RNA seq data filter low count before or after normalization

Hi there,

I just created login, my first time post here : )

I am trying to figure out which step should I filter the gene with low count.
before normalization, after normalization? and how to determine the optimal cutoff?

If you know any literature and document that address my question, please let me know.

thanks a lot.

Q
qliu is offline   Reply With Quote
Old 01-10-2015, 01:54 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,476
Default

Normally you do it after normalization, though there's typically not much difference to doing it before vs. after. Regarding the cut off, please see the genefilter package and the accompanying paper in PNAS.
dpryan is offline   Reply With Quote
Old 01-10-2015, 01:56 AM   #3
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 310
Default

Quote:
Originally Posted by qliu View Post
Hi there,

I just created login, my first time post here : )
Hello!

Quote:
If you know any literature and document that address my question, please let me know.
(Assuming you want to do differential gene expression expression analysis and you have RNAseq or similar type of data)

The edegR vignettes is a document where filtering is applied. I don't know of systemic studies addressing where to set the cutoff.

Quote:
I am trying to figure out which step should I filter the gene with low count.
before normalization, after normalization? and how to determine the optimal cutoff?
The vignettes quoted above apply the filter after transforming raw counts to counts per million.
In practice I think most of us apply a common sense threshold. If say a gene has 1 cpm (counts per million) in all the libraries, that gene can't be biologically interesting. The reason to reject genes early, i.e. before testing for de, is to make the adjustment for multiple testing less "aggressive".
dariober is offline   Reply With Quote
Old 01-10-2015, 02:11 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,476
Default

Quote:
Originally Posted by dariober View Post
I don't know of systemic studies addressing where to set the cutoff.
See the Bourgon et al. 2010 paper in PNAS. That described microarrays, but the same applies to RNAseq. Detection power is optimized by performing the filtering on the p-values after testing. I do occasionally prefilter, but only with limma/voom when I get weird fits (sometimes that happens and you have to remove low expressors to get appropriate results).

BTW, this is why I use DESeq2, it does all of this for me.
dpryan is offline   Reply With Quote
Old 01-10-2015, 02:21 AM   #5
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 310
Default

Quote:
Originally Posted by dpryan View Post
See the Bourgon et al. 2010 paper in PNAS. That described microarrays, but the same applies to RNAseq. Detection power is optimized by performing the filtering on the p-values after testing. I do occasionally prefilter, but only with limma/voom when I get weird fits (sometimes that happens and you have to remove low expressors to get appropriate results).

BTW, this is why I use DESeq2, it does all of this for me.
Thanks - I'll have a look!
dariober 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:45 PM.


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