Question: Deseq2: Time series experiment with 2 conditions and 5 time points
0
gravatar for stan
23 months ago by
stan50
Pretoria, South Africa
stan50 wrote:

Hi all,

I have done a time series analysis using the steps outlined below and I got differentially expressed genes showing a condition-specific effect (according to the example at http://www.bioconductor.org/help/workflows/rnaseqGene/#time)

> sampleTable<-data.frame(sampleName=sampleFiles, fileName=sampleFiles, condition=sampleCondition, time=time)
> sampleTable$condition <- relevel(sampleTable$condition, "Valor")
> sampleTable$time <- factor(sampleTable$time, levels=c("0","6","12","24","72"))
> ddsHTSeq<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design=~condition + time + condition:time)

> ddsTS2 <- DESeq(ddsHTSeq, test="LRT", reduced = ~ condition:time)

> resultsNames(ddsTS1)
 [1] "Intercept"              "condition_BP1_vs_Valor" "time_6_vs_0"           
 [4] "time_12_vs_0"           "time_24_vs_0"           "time_72_vs_0"          
 [7] "conditionBP1.time6"     "conditionBP1.time12"    "conditionBP1.time24"   
[10] "conditionBP1.time72"

> summary(resTS1)

out of 1840 with nonzero total read count
adjusted p-value < 0.1
LFC > 0 (up)     : 138, 7.5%
LFC < 0 (down)   : 113, 6.1%
outliers [1]     : 4, 0.22%
low counts [2]   : 0, 0%
(mean count < 0.2)

But now I am interested in picking the differential genes between the 2 conditions at each time-point, how do I get these? which contrast should i use specifically and which elements from my resultsNames(ddsTS1) output. i checked ?results but still not sure whether to use the condition effect or the interaction effect?

many thanks,

Stan

rna-seq deseq2 R • 4.4k views
ADD COMMENTlink modified 5 weeks ago by deni.ribicic0 • written 23 months ago by stan50
2
gravatar for stan
23 months ago by
stan50
Pretoria, South Africa
stan50 wrote:

So, I found the answer for my question in one of the Deseq2 threads at the bioconductor support forum, a nice explanation is given and the thread is found here: https://support.bioconductor.org/p/65676/#66860 incase it might be of use to someone else.

ADD COMMENTlink written 23 months ago by stan50
0
gravatar for deni.ribicic
5 weeks ago by
deni.ribicic0 wrote:

Hi Stan,

This is a bit late reply, but I am facing actually the same problem. So the answer is basically to separate your time-points and test each and every on certain condition? This seems to be tedious especially if you have a lot of time points!

Cheers, Deni

ADD COMMENTlink written 5 weeks ago by deni.ribicic0

Hi Deni, i doubt doing so will give you pairwise comparisons between the two conditions? perhaps i am missing something in your statement. The answers provided in the link above worked fine for me. Cheers Stan

ADD REPLYlink written 26 days ago by stan50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1495 users visited in the last hour