DESeq2 analysis for low RNA input
2
0
Entering edit mode
@capricygcapricyg-17892
Last seen 21 months ago
United States

Hello, there,

Our cells were treated in two tough conditions (A and B), and as a result, numbers of the survival cells were pretty low, particularly for condition A. The total mapped reads for condition A was about 10% of that for condition B. I wonder if I can still run DESeq2 for differential analysis between A and B.

Thanks a lot!

C.

deseq2 • 599 views
ADD COMMENT
0
Entering edit mode

How about variance-mean dependence for low-reads sample? Does it still hold?

ADD REPLY
0
Entering edit mode

Yes VST stills works if that is what you are asking?

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 32 minutes ago
United States

What are the number of reads for the two groups (in millions)?

ADD COMMENT
0
Entering edit mode

Condition B, we had 0.8-2 millions and condition A, only 0.04-0.15 million.

The numbers of the detected genes for condition B was about 5000 (human genome) and for A about 1500...

ADD REPLY
0
Entering edit mode

I think this is a problem for any kind of DE analysis, because you can get a count of 20 in condition B, but a 0 in condition A under the null. We discuss this as a pathological case in the DESeq2 paper. One way to deal with this would be to filter out genes unless you see a positive count for all samples. This isn't ideal but otherwise you are subject to the problem that 0's which are compatible with the null will look like significant differences.

ADD REPLY

Login before adding your answer.

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