Question: How to identify differential expressed genes if the normalization assumptions are violated
gravatar for afli
2.2 years ago by
China, Beijing, IGDB
afli190 wrote:

Hi, I do the DESeq2 analysis of RNA-seq samples from two different tissues, each with 3 replicates, and both PCA or sample correlation heatmap results is normal and as expected. The total diffrential gene number has reached ~60% (padj<0.1), this did not obey the DESeq2 assumption that most of genes are not diffrentially expressed. So using DESeq2 may not be appropriate.

I search for two related posts: & Differential expression for two very different samples. But I still cannot find the good solution for the case.

I know there must be a large proportion of genes differentially expressed between the tissues, the problem is how to identify them in a correct way? I think many of you may encounter this phenomenon when analysis RNA-seq samples between different tissues. How do you figure out this?

Thank you very much!


edger deseq2 • 1.1k views
ADD COMMENTlink modified 2.2 years ago by Devon Ryan98k • written 2.2 years ago by afli190

With padj<0.1 possibilities that, you are considering a high number of false positive genes. You may try with stringent cut-offs like padj< 0.05 and also Log2FC > 2 (mostly consider in this field, but that vary from research questions). Likewise, you can get some significant amount of genes.

ADD REPLYlink written 2.2 years ago by sangram_keshari240

Thank you for your contribution.

ADD REPLYlink written 2.2 years ago by afli190
gravatar for Devon Ryan
2.2 years ago by
Devon Ryan98k
Freiburg, Germany
Devon Ryan98k wrote:

You run into problems when there are global changes in expression that aren't (roughly) symmetric. It's OK if the majority of genes are DE, as long as they're not more or less changing in one direction. This is the case for all standard normalization methods and not particular to DESeq2. Have a look at the MA plot and next time consider adding some spike-ins to double check that the results are reasonable.

ADD COMMENTlink written 2.2 years ago by Devon Ryan98k

Thank you Ryan for your comment.

ADD REPLYlink written 2.2 years ago by afli190
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1654 users visited in the last hour