DESeq2: interaction term as in limma
1
1
Entering edit mode
sbcn ▴ 70
@sbcn-4752
Last seen 12 weeks ago
Spain

Hi,

In limma's documentation, the interaction term is described as:
"The first of these questions relates to the WT.S vs WT.U comparison and the second to Mu.S vs Mu.U . The third relates to the difference of differences, i.e., (Mu.S-Mu.U)-(WT.S-WT.U) , which is called the interaction term."

I want to perform the same analysis using DESeq2 (on RNA-seq data).
I have 2 cell types: A and B
for each cell type, I have a KO of a protein and a control (no KO)
 
I looked at the vignette and at ?results, but I have some doubt regarding what would be equivalent to limma's interaction term (A_ko-A_control)-(B_ko-B_control)

Could someone tell me if the following attempt is correct:

dds <- makeExampleDESeqDataSet(n=100,m=12)

dds$celltype <- factor(c(rep("A",6), rep("B", 6)))
dds$protein <- factor(rep(c(rep("ko",3), rep("control", 3)), 2))

design(dds) <- ~ celltype + protein + celltype:protein

dds <- DESeq(dds)
#resultsNames(dds)

res <- results(dds, contrast=list( c("protein_ko_vs_control","celltypeB.proteinko") ))

Thanks a lot!

 

deseq2 deseq limma interaction term • 1.3k views
ADD COMMENT
4
Entering edit mode
Simon Anders ★ 3.7k
@simon-anders-3855
Last seen 14 months ago
Zentrum für Molekularbiologie, Universi…

Nearly correct.

The last line should be:

res <- results( dds, name="celltypeB.proteinko" )

 

Have a look at resultsNames:

> resultsNames(dds)
[1] "Intercept"             "celltype_B_vs_A"       "protein_ko_vs_control"
[4] "celltypeB.proteinko"  

 

Celltype_B_vs_A is the difference between the cell lines A and B in control (B_ctrl - A_ctrl), and protein_ko_vs_control is the difference between the ko and ctrl in cell line A (A_ko - A_ctrl). The interaction term, celltypeB.proteinko, is the difference between the effect of the knockdown in cell line B over that in A, i.e., (B_ko - B_ctrl) - (A_ko - A_ctrl), which is what you want.

What you looked at ist the sum the last two terms, i.e., which turns out to be (B_ko - B_ctrl).

R's way of naming interaction term is a bit peculiar. Read "celltypeB.proteinko" as "the additional effect that the knockout has if done on cell line B rather than A", or, equivalently, "the additional part of the difference between A and B, when assessed in ko rather than ctrl".

 

 

ADD COMMENT
0
Entering edit mode

Thanks a lot Simon for your prompt answer!

It is a lot more clear now to me what each term means.

If a gene is differentially expressed in B_ko - B_ctrl but not in A_ko - A_ctrl (change in expression triggered by the KO in cell line B but not in cell line A), it will come out of the interaction analysis also as significant, correct?

 

ADD REPLY
3
Entering edit mode

Not necessarily, because "absence of evidence is not evidence of absence" of an effect.

So, if the p value for (B_ko-B_ctrl) is high, this could either mean that there is no effect or that the variance is too high to tell. Hence, seeing a significant p value for (A_ko-A_ctrl) but not for (B_ko-B_ctrl) does not imply that there is a difference between (A_ko-A_ctrl) and (B_ko-B_ctrl).

This is precisely the reason why you have to test for significance of interaction and why it would be incorrect to simply run the analysis independently on the A and the B samples and look at the genes appearing only in one of the two results lists.

ADD REPLY
0
Entering edit mode

Thank you Simon for the clarification. Then the interaction is indeed what I am looking for.

Best,

 

ADD REPLY

Login before adding your answer.

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