Question: DESeq2 comparing differential expression results
0
gravatar for rrcutler
2.1 years ago by
rrcutler50
rrcutler50 wrote:

Hi all,

I am attempting to meaningfully compare two resulting list of significantly (padj < 0.05) differentially expressed genes in an experiment with a shared control, although I am getting unexpected results. This experiment used 5 replicates for each treatment and 5 shared controls.

For example, when comparing two different treatments against the control, this yields 1025 and 579 DE genes respectively. When comparing these two different treatments against each other, this only yields 18 DE genes. In addition, when looking at the shared genes of the resulting list when comparing the aforementioned two different treatments against the control, this yields 329 DE genes.

This is confusing, as it seems the two treatments are not that different from each other when directly comparing them against each other. But when comparing to the shared control independently, this tells a much different story. My question is - what may be the reason for this unexpected difference?

In addition, is there any available software or a better method to deal with these types of comparisons?

Thanks,

-R

edger deseq2 • 422 views
ADD COMMENTlink modified 2.1 years ago by Michael Love25k • written 2.1 years ago by rrcutler50
Answer: DESeq2 comparing differential expression results
2
gravatar for Michael Love
2.1 years ago by
Michael Love25k
United States
Michael Love25k wrote:

hi R,

A and B can both be found to have significant differences when comparing to C, but might not be found to have significantly differences from each other. This isn't an inconsistency.

Imagine below I draw the range of sample gene expression within a group (just one gene):

A: [====]
B:     [====]
C:                   [====]

So here A and B might both be found significantly different than C, but their range of expression values overlap, so perhaps you don't find them significantly different from each other (you'd need more samples perhaps).

The numbers you have above don't need to add up to the same number, because there are other genes which may be different but not significantly so according to a test.

ADD COMMENTlink written 2.1 years ago by Michael Love25k

Thanks Micheal, this makes more sense now. 

ADD REPLYlink written 2.1 years ago by rrcutler50
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 16.09
Traffic: 231 users visited in the last hour