Design for introducing batch effect in DESeq2
1
0
Entering edit mode
ag1805x ▴ 50
@ag1805x-15215
Last seen 8 days ago
University of Allahabad

I have RNAseq data taken from two different ENA project IDs (PrjA & PrjB). Both PrjA & PrjB contains test and control samples. PrjA was PE 50bp sequencing while PrjB was SE 100bp. Is it possible we can combine the data to identify DEG between test and control? Would batch effect removal by Combat be a better option or use the design = project+condition in DESeq2?

deseq2 combat batch effect rna-seq edger • 338 views
ADD COMMENT
0
Entering edit mode

If you DO know the batches, then what swbarnes2 recommended is the best option. If you don't know then you could use sva.

In your case, however, I don't think that it's wise to just combine the data because it seems that you use different platforms for the sequencing. Could you show us the MDS plot?

ADD REPLY
1
Entering edit mode
swbarnes2 ▴ 840
@swbarnes2-14086
Last seen 42 minutes ago
San Diego

Do what the DESeq software recommends. Give DESeq raw counts, not some kind of batch corrected, and include batch as a factor in the design.

ADD COMMENT
1
Entering edit mode

You can do that even when the data came from two different projects? Seems like the batch issue in this case is not simply multiple runs of the same instrument but other things are different as well.

ADD REPLY

Login before adding your answer.

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