RNA seq time course design using DESeq2
Entering edit mode
Last seen 5 months ago

Hello! I am having a difficult time with the design element of my Deseq2 timecourse following the vignette in (Love MI, Anders S, Kim V and Huber W. RNA-Seq workflow: gene-level exploratory analysis and differential expression [version 1; referees: 2 approved]. F1000Research 2015, 4:1070 (doi: 10.12688/f1000research.7035.1)) as a template.

My experiment is looking at Controls (wild type mice) Vs Experimental (knock out) with 4 biological replicates each, at 4 time points (p0, p7, p15, p30) These are not repeated measures on the same group of mice (each mouse was sac'ed at each timepoint). My hope is to look at condition specific differences over time. Here is what I have so far:

sampleFiles<- c("CB-1712-P0-M-KO.counts.processed.txt",

sampleNames <- c("CB-1712-P0-M-KO",

Genotype <- c("Exp","Control","Control","Control","Exp",
time<- c("0","0","0","0","0","0","0","0","7","7","7","7","7","7","7","7",

sampleTable<-data.frame(sampleName=sampleNames, fileName=sampleFiles, condition=Genotype, time=time)

ddsHTSeq<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design=~ condition + time + condition:time)
keep <- rowSums(counts(ddsHTSeq)) > 10
dds <- ddsHTSeq[keep,]

dds$condition <- relevel(dds$condition, ref = "Control")
ddsTC <- DESeq(dds, test="LRT", reduced = ~ condition + time)
resTC <- results(ddsTC)
resTC$symbol <- mcols(ddsTC)$symbol

"Intercept"     "condition_Exp_vs_Control" "time_15_vs_0"            
[4] "time_30_vs_0"             "time_7_vs_0"              "conditionExp.time15"     
[7] "conditionExp.time30"      "conditionExp.time7"

This seems to not be comparing what I want it to. I would like it to compare controls vs exp with time as an interaction term but it seems that it is not grouping the analysis the way I would like. Any help would be greatly appreciated!

deseq2 timecourse • 141 views
Entering edit mode
Last seen 1 minute ago
United States

"I would like it to compare controls vs exp with time as an interaction term but it seems that it is not grouping the analysis the way I would like"

It looks correct to me, so I think if the terms are unfamiliar, you should collaborate with a local statistician to help with the interpretation of results. I unfortunately don't have sufficient time to provide statistical design or interpretation support here, but I have to reserve my time only for software support.

Entering edit mode

Thanks for your quick response! Is there a part of the manual/vignette that explains the output of "resultsNames(ddsTC)" so I can better understand what comparison each header is referring to? Thanks

Entering edit mode

The vignette and workflow have quite a bit of information about how to interpret results, there is an interaction section in the vignette.

But I would not recommend to attempt to use interaction models without any background on what interactions do in linear models. At some point, we have to assume that the user is aware of what the model terms are doing. It's much better to consult with a statistical collaborator to make sure you are interpreting the results correctly.


Login before adding your answer.

Traffic: 243 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6