Differential analysis of multiple samples with different time points
1
0
Entering edit mode
Firas • 0
@firas-17791
Last seen 12 months ago

Hello, 

I have the following samples, with three replicates in each of them: 

Time point Population Genotype1 Genotype2
00h Precursor WT KO
06h Subset A WT KO
18h Subset A WT KO
72h Subset A WT KO
06h Subset B WT KO
18h Subset B WT KO
72h Subset B WT KO

My biological setup:

  • I started with 2 samples of Precursor (considered also as 00h): KO or WT. 
  • Then differentiated them under two different conditions (Subset A and Subset B).
  • The samples at different time points are collected from separated wells.

What I would like to analyse:

  • the effect of WT vs KO across all samples
  • taking the WT alone, the effects of Subset A vs Subset B
  • taking the WT alone, the effects of Precursor vs Subset A or Subset B  

I would like to get some help on how to set up contrasts for this experiment. 

Many thanks

diffbind R atac-seq • 572 views
ADD COMMENT
3
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 19 hours ago
CRUK, Cambridge, UK

If you want to track the differences across time between two groups in a single comparison, DiffBind is not really set up to do that. You'd need to use one of the underlying analysis methods, like DESeq2, directly. There are a number of threads discussing how to set up time series analyses in DESeq2

There are several things you can do in DiffBind, provided you set up the metadata in a useful way.For example, if you set:

  • DBA_TISSUE to be Precursor, SubsetA, or SubsetB
  • DBA_CONDITION to be WT or KO (you can rename it Genotype for display purposes)
  • DBA_TREATMENT to be T00, T06, T18, or T72 (you can rename it Timepoint for display purposes)
  • DBA_REPLICATE to be 1, 2, or 3

then you can compare:

  • Overall differences between WT and KO, across all time points but independent of genotype:
dba.contrast(myDBA, group1=myDBA$masks$WT, group2=myDBA$masks$KO, 
                    block=DBA_TISSUE)
  • Overall differences between WT and KO, across either genotype but independent of time:
dba.contrast(myDBA, group1=myDBA$masks$WT, group2=myDBA$masks$KO, 
                    block=DBA_TREATMENT)
  • Differences between Genotypes A and B in WT, independent of time point:
dba.contrast(myDBA, group1=myDBA$masks$WT&myDBA$masks$SubsetA, 
                    group2=myDBA$masks$WT&myDBA$masks$SubsetB,
                    block=DBA_TREATMENT)
  • Differences between Genotypes A and B in WT at a specific timepoint:
dba.contrast(myDBA, group1=myDBA$masks$WT&myDBA$masks$SubsetA&myDBA$masks$T18,
                    group2=myDBA$masks$WT&myDBA$masks$SubsetB&myDBA$masks$T18)
  • Differences between Precursor and each of the Genotypes in WT, independent of time point:
dba.contrast(myDBA, group1=myDBA$masks$WT&myDBA$masks$Precursor, 
                    group2=myDBA$masks$WT&myDBA$masks$SubsetA,
                    block=DBA_TREATMENT)
dba.contrast(myDBA, group1=myDBA$masks$WT&myDBA$masks$Precursor, 
                    group2=myDBA$masks$WT&myDBA$masks$SubsetB,
                    block=DBA_TREATMENT)
  • And at a specific timepoint:
dba.contrast(myDBA, group1=myDBA$masks$WT&myDBA$masks$Precursor&
                           myDBA$masks$T18,
                    group2=myDBA$masks$WT&myDBA$masks$SubsetA&
                           myDBA$masks$T18)
ADD COMMENT
0
Entering edit mode

Thank you very much, Rory, for your swift reply. 

This is exactly what I want. Best wishes, Firas

ADD REPLY

Login before adding your answer.

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