Question: How to identify differential expressed genes if the normalization assumptions are violated
0
gravatar for afli
10 months ago by
afli190
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: https://support.bioconductor.org/p/94661/ & 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!

Aifu.

edger deseq2 • 332 views
ADD COMMENTlink modified 10 months ago by Devon Ryan92k • written 10 months ago by afli190
1

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 10 months ago by sangram_keshari220

Thank you for your contribution.

ADD REPLYlink written 10 months ago by afli190
1
gravatar for Devon Ryan
10 months ago by
Devon Ryan92k
Freiburg, Germany
Devon Ryan92k 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 10 months ago by Devon Ryan92k

Thank you Ryan for your comment.

ADD REPLYlink written 10 months ago by afli190
Please log in to add an answer.

Help
Access

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