number of grouping affect log2 fold change in DESeq2 analysis?
1
1
Entering edit mode
8.2 years ago
star ▴ 350

I have RNA-seq data for samples in 10 different groups and I would like to find deferentially expressed genes. I have used Deseq2 R package. The deferentially expressed genes have FDR between 0.03 and 0.05. I checked log2FC and found it is low value in most deferentially expressed genes. How I can explain this? Why genes with low log2FC are found to be deferentially expressed? Also I have checked the normalized count and I couldn't find any reason!

R DESeq2 edgeR log2FC • 2.2k views
ADD COMMENT
0
Entering edit mode

How come no differentially expressed gene has fdr below 0.03? I guess you have some kind of filtering going on ..

ADD REPLY
3
Entering edit mode
8.2 years ago

The number of groups you're using suggests that there's a lot of variance shrinkage going on, which is why you're picking up low fold-change differences. My guess is that you have relatively few replicates per group, which is causing aggressive fold-change shrinkage.

ADD COMMENT

Login before adding your answer.

Traffic: 1682 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6