I'm using DiffBind to analyse some H3K4me3 ChIP-seq data. I've used DiffBind for H3K27ac ChIP-seq data as well, and when I drew MA plots using these data they were nicely symmetrical and centered around 0. However, this is not the case with my H3K4me3 data: It looks like the central line of the MA plot is not horizontal, but diagonal, with a lot of positive log fold changes at low "concentrations" (peak height) and negative log fold changes at high concentrations (see attached figure). When I use these data for differential binding analysis with DiffBind's default settings, both these groups (low conc/positive LFC and high conc/negative LFC) contain lots of differentially bound sites. However, when I use edgeR instead, this is not the case and only about a quarter of the sites are differentially bound (all low conc/positive LFC).
I have matching RNA-seq data, and these do not suggest that there is a global shift in transcription or anything. I also know that there is a difference in experimental efficiency between the two groups I am comparing (disease v. healthy). When I use
dba.analyze() with DESeq2 to compensate for this, I get results very similar to the edgeR result. Basically, I don't have other evidence to suggest there is a global change in signal between my two conditions (but I have not systematically tried to rule this out in the lab).
But even with stringent normalisation, the central line in the MA plot is diagonal. Does anyone know if this means anything, and what the consequences are for the normalisation method I should apply?
Any help would be much appreciated.