Search
Question: Differential binding of ATAC-seq data with no replicates
1
gravatar for kennethcondon2007
7 months ago by
kennethcondon200710 wrote:

Hi,

I have ATAC-seq data for 8 time points of 4 conditions --> 32 samples in total. The only replicates I had were technical which were pooled early in the pipeline. There are no biological replicates. I have completed MACS2 peak calling and now want the differential binding (DB). I would like to use HTseq to create the counts table and then DEseq2 to analyze for DB.

I have read a number of posts all replied with good advice from Micheal Love (DEseq2 extended to open chromatin anlaysis: normalization, dispersion fit, and "too many" differences, Replicate for DESeq2 and A: GFOLD file as input for DESeq2) regarding DEseq2, ATAC-seq and replicates....but none that address what to do when there are no biological replicates. Can HTseq and DEseq2 run with no replicates?

Also I'm a bit confused as to what needs to be counted. I have the BAM files and the output of MACS2 I have filtered to remove overlapping peaks (as per guidance). So do I simply count the number of reads in the original data that align to the remaining non-overlapping uniquely named peaks (separately for each sample) - is that what the count table (of each individual sample) should consist of? For example:

peak_name     read count

peak1                  452

peak2                  34

peak3                  458

 

I will be monitoring this post so let me know if I can clarify question(s) or provide additional information.

Thank you all for your help,

Kenneth

ADD COMMENTlink modified 6 months ago • written 7 months ago by kennethcondon200710
3
gravatar for Gordon Smyth
7 months ago by
Gordon Smyth31k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth31k wrote:

Actually the approach of calling peaks separately in each sample would lead to overestimation of significance in a DB analysis even if you did have replicates. See

https://doi.org/10.1093/nar/gku351

for an explanation of why this is so.

ADD COMMENTlink modified 7 months ago • written 7 months ago by Gordon Smyth31k

Thank you Gordon. That's my next read.

 

EDIT With 8 time points per condition, should I be pooling all 8 into a single file and peak call on that file?

ADD REPLYlink modified 7 months ago • written 7 months ago by kennethcondon200710

No, you would still be over-stating the differences between conditions. To follow the advice in our NAR article, you would pool all 32 libraries together for peak calling.

ADD REPLYlink written 7 months ago by Gordon Smyth31k

Hi Gordan,

I really need to flesh this out a bit. It is quite a headache. So...

1) For downstream processing I need a consensus peak set. So ALL 32 samples must be merged. Then call peaks on the merged data to identify a consensus set of peaks. Then for each sample, count the number of 5’ read ends that are located within each peak region of the consensus peak set. Differential analysis is then carried out by comparing the counts of each individual sample.

2) To merge sample before peak calling - I assume I can just pass multiple files to MACS2 (-t A B C) rather than running bedtools merge.

3) Considering time and the fact that I have already called the peaks for each sample individually. Could I not just merge the peak files to get a consensus peak set instead as suggested by Devon Ryan and Sukhdeep Singh.

Kenneth.

ADD REPLYlink modified 7 months ago • written 7 months ago by kennethcondon200710

May I ask if you came to an answer to this question?

ADD REPLYlink written 10 weeks ago by ATPoint0
1
gravatar for Michael Love
7 months ago by
Michael Love14k
United States
Michael Love14k wrote:

Honestly, I don't think it's much use to run DESeq2 without replicates. The software allows for this analysis with large caveats (see the relevant section in ?results) that such an analysis is only vaguely exploratory. But you could just as well plot the log ratios of peak counts across conditions and skip running DESeq2.

ADD COMMENTlink written 7 months ago by Michael Love14k

Thank you Micheal - I was quite alarmed when I was informed there were no replicates, but there is nothing I can do about it now except remind them to include a bioinformatician during the planning stage.

ADD REPLYlink written 7 months ago by kennethcondon200710
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.2.0
Traffic: 293 users visited in the last hour