MACS2: Missing Peaks in Bed file but Present in Bedgraph/Bigwig?
Entering edit mode
3.0 years ago
Ankit ▴ 350

Hi everyone,

I am performing allele specific ChIP-Seq analysis. First I split allele specific reads using SNP positions. This I did both for Input and ChIP sample. Then I used MACS2 for peak calling as follows:

macs2 callpeak -t Sample_rep1_allele1.bed -c Input_rep1_allele1.bed -B -f BED -g mm -n Sample_rep1_allele1
macs2 callpeak -t Sample_rep1_allele2.bed -c Input_rep1_allele2.bed -B -f BED -g mm -n Sample_rep1_allele2
macs2 callpeak -t Sample_rep2_allele1.bed -c Input_rep2_allele1.bed -B -f BED -g mm -n Sample_rep2_allele1
macs2 callpeak -t Sample_rep2_allele2.bed -c Input_rep2_allele2.bed -B -f BED -g mm -n Sample_rep2_allele2

#Merge both the replicates 

Later, I tried to look for monoallelic peaks (peaks unique to one of the allele). I used bedtools to get the peaks which are not shared between the two alleles:

bedtools intersect -wa -wb -a Sample_rep1-rep2-allele1_peaks.bed -b Sample_rep1-rep2-allele2_peaks.bed -v  > Sample_allele1_specific_peaks.bed
bedtools intersect -wa -wb -b Sample_rep1-rep2-allele1_peaks.bed -a Sample_rep1-rep2-allele2_peaks.bed -v  > Sample_allele2_specific_peaks.bed

When I visualised some of the random peaks on UCSC browser by uploading bigwig file, I noticed some of the peaks which were unique to allele1 was also present in allele2 with almost equal height and not present in either of the bed files (which made me consider them as monoallelic).

At some of the region's it was completely unique and doesn't overlap means they were completely monoallelic.

I am not sure why they all were not unique (monoallelic) when I tried to extract peaks only specific to one of the allele using bedtools (-v).

I thought it might be possible in either of peak calling step (allele1 or allele2) that peak was not called as the peak, may be because of some internal MACS2 calculations or during normalisation with respective controls. However I am not sure for the exact issue and the solution for this.

I hope I am clear with my question.

It will be very helpful if someone can suggest me what I am doing wrong or missing something.


macs2 ChIP-Seq bed bigwig allele-specific • 1.3k views
Entering edit mode

Can you add a representative screenshot?


Login before adding your answer.

Traffic: 1230 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6