Too many DEGs
1
0
Entering edit mode
Raymond ▴ 20
@raymond-14020
Last seen 4.8 years ago

Hi, I have 36 animals and got two tissues of each animal for RNAseq.

My design model is: sample+tissue

and I got DEG results like that:

out of 16445 with nonzero total read count adjusted p-value < 0.05 LFC > 0 (up) : 7276, 44% LFC < 0 (down) : 7508, 46% outliers [1] : 0, 0% low counts [2] : 0, 0% (mean count < 4)

90% of genes are DEGs at the level of 0.05. Is there any problem with my results?

Below I attached two figures:

1.Normalization figure. I plot the density function of the gene vst expressions in each sample. The figures showed that two tissues are different, but all normalized to the median.

Screen-Shot-2019-06-14-at-11-55-28-AM

  1. Gene expression data. I rank the average gene expression of genes in Tissue1, and the then plot all genes in Tissue2 based on the ranking order in Tissue1.

Screen-Shot-2019-06-14-at-12-46-03-PM

From this figure, it seems that a lot of genes between these two tissues are quite close to each other, there should not be so many DEGs.

Any suggestions?

Thanks for your time & regards,

Raymond

deseq2 • 522 views
ADD COMMENT
0
Entering edit mode

Can you edit this post and show the code you used for the analysis?

ADD REPLY
1
Entering edit mode
@mikelove
Last seen 18 hours ago
United States

You are testing a point null and have many samples. We discuss this situation in the DESeq2 paper. Why not specify an LFC threshold?

ADD COMMENT

Login before adding your answer.

Traffic: 682 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