Question: What will lead DESeq2 shows no DEGs if workflow adopted is okay?
0
gravatar for gmumuxi
3 months ago by
gmumuxi0
gmumuxi0 wrote:

Hi,

ddsFullCountTable <- DESeqDataSetFromMatrix(countData = mouse_count, colData = mouse_pheno, design = ~ batch + age + gender + treatment)  

dds <- DESeq(ddsFullCountTable, test="LRT", reduced = ~ batch + age + gender)  

res_dds <- results(dds, contrast = c("treatment", "med", "ck"))

summary(res_dds) 
# out of 32102 with nonzero total read 
# adjusted p-value < 0.1
# LFC > 0 (up)       : 0, 0%
# LFC < 0 (down)     : 0, 0%
# outliers [1]       : 1269, 3.6%
# low counts [2]     : 0, 0%
# (mean count < 0)
# [1] see 'cooksCutoff' argument of ?results
# [2] see 'independentFiltering' argument of ?results

Here are 5 med mice and 7 check mice (so the ncols is 12 for countData) RNA-Seq data be used to do DE analysis, to find the DEGs affected by medicine, but seems no genes. This problem has confused me for a long time, I checked and tired many times, still haven't been solved. As far as I understand it, I guess:

  1. Too few samples (and too many covariates?)?
  2. The drug not statistically significant in this case?
  3. The code, I guess not?

Or other problems?

Thanks!

deseq2 • 136 views
ADD COMMENTlink modified 3 months ago by arfranco130 • written 3 months ago by gmumuxi0
Answer: What will lead DESeq2 shows no DEGs if workflow adopted is okay?
0
gravatar for Michael Love
3 months ago by
Michael Love24k
United States
Michael Love24k wrote:

Nothing wrong with the code.

There do seem to be a lot of outlier counts. Have you looked at plotPCA?

ADD COMMENTlink written 3 months ago by Michael Love24k

Thanks for your reply! Since outliers here is 1269, it mean the "outlier" are genes, not samples. However I considered the outlier problem, it seems no obvious steps(function or tool) that to quantitative identification and remove them in DESeq2 wokflow?

ADD REPLYlink written 3 months ago by gmumuxi0

Take a look at the vignette, these are removed by DESeq2. It explains in detail.

ADD REPLYlink written 3 months ago by Michael Love24k
Answer: What will lead DESeq2 shows no DEGs if workflow adopted is okay?
0
gravatar for arfranco
3 months ago by
arfranco130
European Union
arfranco130 wrote:

Run a plotPCA like Michael recommends you

I had an experiment with 2 conditions and 3 biological replicates. By running a PCA analysis, I noticed that two samples were exchanged (a control was actually a treated sample, and viceversa), and this was enough for not getting a single DE gene

ADD COMMENTlink modified 3 months ago • written 3 months ago by arfranco130
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: 211 users visited in the last hour