Hello everyone!
We do ChIP-seq across 2 conditions each with 2 replicates.
The number of called narrow peaks (MACS2 2.2.7.1, FDR < 0.1) are:
cond_1-1: 8,915
cond_1-2: 10,099
cond_2-1: 1,446
cond_2-2: 801
At first, I suspected that there are many small peaks detected in condition_1.
However, I observed the opposite as the size of called narrow peaks tend to be larger in condition_1.
(1) ~24% of peaks are longer than 600 bps in condition_1 while that number is ~13% for condition_2.
(2)the median peak size is ~400 bps in condition_1 but only ~330 bps in condition_2.
So, usually what could be the reason to call many more peaks in one sample than the other? the library sizes are similar across samples.
Thanks in advance!
Pei
•
written
by
Pei
1 answer
Either its biological, where your factor's binding is greatly decreased in condition two. Or it's technical where the IP did not work as well for samples in condition 2. If all samples were not prepared in parallel, I would be more wary of the latter.
I think it would be important to check other QC aspects of the libraries, like fragment sizes and replicate correlation. Also always good to visualize the signal on a genome browser.
Log in to answer this question.
-
Overestimate of ATAC-seq differential peaks
written by Darryl •
-
I do not understand how to get total peaks from MACS2 shown to me. What mistake(s) am I making here?
written by Ronin •
-
MACS2 peak caller omitting large chunks of peaks
written by Nevergivingup •
-
Comparing peaks between multiple samples
written by brisbio •
-
Macs2: predicted fragment length less than 2*tag length
written by langya •
-
Peak calling using macs2?
written by star
-
Diffbind - Very high proportion of significant peaks
written by kautilya
-
ChIP-seq differential peak calling (Any way to set the window size for calling differential enriche…
written by chiefcat
-
Difference in number of peaks called
written by anu014
-
How do I select top 100,000 non-overlapping peaks from MACS2 narrow peaks output?
written by BioinfGuru
