Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • DESeq2 multi-factor design

    Hi all,

    I have just started using DESeq2 using a multi-factor approach. As I'm new to this (and to statistics), I was hoping to get feedback / corrections on what I think the analysis should look like:

    I have the following colData:
    Code:
    sample	genotype	treatment	library_batch
    S1	WT	vehicle	B1
    S2	WT	vehicle	B2
    S3	WT	vehicle	B1
    S4	WT	vehicle	B2
    S5	WT	Drug	B2
    S6	WT	Drug	B1
    S7	WT	Drug	B2
    S8	KO	vehicle	B1
    S9	KO	vehicle	B2
    S10	KO	vehicle	B1
    S11	KO	vehicle	B2
    S12	KO	Drug	B1
    S13	KO	Drug	B2
    S14	KO	Drug	B1
    S15	KO	Drug	B2
    Blocking for library batch effects, I am interested in finding the following DE genes:
    Q1) DE genes between vehicle-treated "WT" and "KO" (i.e. the 'baseline' difference);
    Q2) Genes that are DE following treatment within each genotype ("WT" and "KO");
    Q3) Genes that show a genotype-dependent response to treatment.

    Looking through the DESeq2 tutorials and form posts, Q1 and Q2 would probably be easiest addressed by created a 'genotype_treatment' factor and working with contrast:

    Code:
    colData$condition <- paste(colData$genotype, colData$treatment, sep="_")
    and using a design
    Code:
    design = ~ library_batch + condition
    to extract results as:
    Q1:
    Code:
    res1 <- results(dds, contrast=c("condition", "KO_vehicle", "WT_vehicle"))
    Q2.1:
    Code:
    res2wt <- results(dds, contrast=c("condition", "WT_Drug", "WT_vehicle"))
    Q2.2:
    Code:
    res2ko <- results(dds, contrast=c("condition", "KO_Drug", "KO_vehicle"))
    As for Q3, I would then repeat the analysis including an interaction term:

    Code:
    design = ~ library_batch + genotype*treatment
    to extract results as:

    Code:
    res3 <- results(dds, name="genotypeKO.treatmentvehicle")
    Is this correct? I was wondering if there is a way to do this without running 2 separeate analysis designs...?

    On a more general note, can anyone recommend a resource to get more familiar with these design formulas and what they 'mean'?
    Last edited by PMRoberts; 11-16-2015, 03:48 PM.

  • #2
    Q1/2 looks very similar to what we've done with our RNASeq experiments. Mike Love seemed to think the approach was correct:

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X