Question: MACS2: Missing Peaks in Bed file but Present in Bedgraph/Bigwig?
0
gravatar for Ankit
6 months ago by
Ankit120
Ankit120 wrote:

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:

#Rep1
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
#Rep2
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 
Sample_rep1-rep2-allele1_peaks.bed 
Sample_rep1-rep2-allele2_peaks.bed

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.

Thanks

ADD COMMENTlink written 6 months ago by Ankit120

Can you add a representative screenshot?

ADD REPLYlink written 6 months ago by ATpoint26k
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.3.0
Traffic: 1268 users visited in the last hour