Question: DESeq2 Design Formula for Time Series Experiment
1
gravatar for gkuffel
2.1 years ago by
gkuffel10
Loyola University
gkuffel10 wrote:

Hi everyone,

I have miRNA-Seq data for an experiment with the following parameters:

SampleID     Treatment     Timepoint     Patient

   P1_0                  TB                 0             1

   P1_1                  TB                 1             1

   P1_2                  TB                 3             1

   HC1_0                None            0              Control1

   HC1_2                None            3              Control1

   P2_0                  TB                 0              2

   P2_1                  TB                 1              2

  P2_2                   TB                 2              2

   P2_3                  TB                 3              2

   HC2_0                None            0              Control2

  HC2_3                 None            3              Control2

I am really struggling figuring out how to enter the design formula. There are 2 patients each with a baseline sample before treatment and 2 samples at time points after treatment. Each patient has an age-matched control with samples taken at two time points. One patient is a full responder to treatment and one is not but I am blinded to that info. The investigator wants to identify DEmiRNAs in the following:

1. Between patients 

2. As a result of treatment over time

Do I need to include interaction terms? So far I have identified 128 DEmiRNAs using this formula:

design=~timepoint + individual

Which I believe is testing for differences between individuals across all time points, can anyone provide some much needed guidance?

deseq2 • 457 views
ADD COMMENTlink modified 2.1 years ago by David Requena20 • written 2.1 years ago by gkuffel10
Answer: DESeq2 Design Formula for Time Series Experiment
2
gravatar for David Requena
2.1 years ago by
The Rockefeller University. New York, USA.
David Requena20 wrote:

To compare time points, you should consider "timepoint" in your formula, as you did.

To compare treated vs untreated, you should include "treatment" in your formula.

If you think that your treated patients will experience changes over time, which are different than the changes experienced by the untreated, you should add the interaction term "treatment:timepoint".

I guess that "individual" is your column "patient" or "treatment:age", which are equivalent.

 

I think you'll might find the explanation in the DEseq2 vignette useful (section "Interactions"):

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#interactions

And also these posts:

A: DESeq2 paired multifactor test

With DESeq2 "Not full rank" Error with design ~ line + time + condition

 

David

ADD COMMENTlink modified 2.1 years ago • written 2.1 years ago by David Requena20

 

David,

Thank you so much for your reply. I have a few more questions if you have time:

1. So with my current design (design =~ Timepoint + Patient) am I testing for differences between patients while controlling for variability caused by time?

2.  If I use the design: design = ~ Timepoint + Treatment, is that testing for differences between treatments while controlling for variability caused by time?

3. When you ask if I think the treated patients will experience changes over time different than changes experienced by untreated do you mean changes caused by the treatment or just naturally ocurring changes?

4. Last and maybe most important question. The most important comparsison the group wants to make is does the treatment have an effect on the patients over time. I still don't understand if I am using the correct design for this.

 

ADD REPLYlink written 2.1 years ago by gkuffel10
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 16.09
Traffic: 100 users visited in the last hour