Differential Expression with Factorial Design and adjust for Lane+Batch Effect
0
0
Entering edit mode
5.4 years ago
Sam ★ 3.4k

Hi, so recently we have got a data back with the following design

 Condition Treatment Lane Batch Case Treated 1 1 Case Treated 2 1 Case Treated 1 4 Control Treated 2 4 Control Treated 1 1 Control Treated 2 1 Case Untreated 1 4 Case Untreated 2 3 Case Untreated 2 2 Control Untreated 2 1 Control Untreated 1 3 Control Untreated 1 2
data.frame(Condition=rep(rep(c("Case","Control"), each=3),2), Treatment=rep(c("Treated", "Untreated"),each=6), Lane=c(1,2,1,2,1,2,1,2,2,2,1,1), Batch=c(1,1,4,4,1,1,4,3,2,1,3,2))

And we would like to test the following

1. Compare effect of treatment in case

2. Compare effect of treatment in control

3. Compare and contrast the effect of treatment in case when compared to control.

With DESeq, we have used the design ~Batch+Lane+Condition+Treatment+Condition:Treatment where

dds$Condition<- relevel(dds$Condition, ref="Control")
dds$Treatment<- relevel(dds$Treatment, ref="Untreated")

and from reading the documentation of DESeq2 and limma user guide, my understanding is that the following usage of results might provide the desire statistic outcome

1. Compare effect of treatment in case:

results(dds, list(c("Condition_Case_vs_Control","ConditionCase.TreatmentTreated")))

2. Compare effect of treatment in control:

results(dds, contrast=c("Treatment", "Treated", "Untreated"))

3. Compare and contrast the effect of treatment in case when compared to control.

results(dds, name="ConditionCase.TreatmentTreated")

Or for 1, will it be better if we

dds$Condition<- relevel(dds$Condition, ref="Case")

and do

results(dds, contrast=c("Treatment", "Treated", "Untreated"))

Is my understanding correct?

RNA-Seq deseq2 • 2.2k views
0
Entering edit mode

I'm pretty sure you want "batch" to be a factor (likewise with "lane", though that won't have a big effect in this case)...

0
Entering edit mode

Oh, right, forgot to set them to factor in this snappy. So other than that, is my use of those code correct??

So for example, in

results(dds, list(c("Condition_Case_vs_Control","ConditionCase.TreatmentTreated")))


I am not sure if DESeq2 will return the p-value of effect of treatment in case adjusting for batch and lane or will DESeq2 return the p-value of effect of treatment in case under the base case of batch and base case of lane. What we would like to get will be the former results but are worry that the later is what we get

1
Entering edit mode

If you include "batch" and "lane" in the design then you will only ever receive results corrected for them. I'd need to double check whether the contrasts you're using do what you want, though, since you're using a factorial design but then aren't asking factorial questions (yes, you can do this, but at least I personally find it more confusing to do that).

0
Entering edit mode

Thank you Devon, will it be better if the design is:

~Batch+Lane+Condition+Condition:Treatment