Question: DiffBind DESeq2 MA plot skewed
gravatar for rbronste
14 months ago by
rbronste60 wrote:


Wondering what can cause skewing in an MAplot as seen below?

ADD COMMENTlink modified 14 months ago by Rory Stark2.6k • written 14 months ago by rbronste60
gravatar for Michael Love
14 months ago by
Michael Love20k
United States
Michael Love20k wrote:

The MA plot for ChIP-seq is always going to look different than RNA-seq, for various reasons. Here on the horizontal you have comparable peaks, and the diagonal pointing up and right shows those peaks which are only present in whatever group is in the numerator of the LFC, with little to no signal in the group that's in the denominator. (For normalization for ChIP-seq, I could see making small alterations to the default RNA-seq approach to estimating size factors. But I don't have any code for this.)

ADD COMMENTlink written 14 months ago by Michael Love20k
gravatar for Rory Stark
14 months ago by
Rory Stark2.6k
CRUK, Cambridge, UK
Rory Stark2.6k wrote:

As Michael Love says, the diagonals represent sites that are at or near zero on one condition but have higher binding in the other. What is unusual is that the most dense part of the plot lies on one of these diagonals (sites with binding in group1 but no binding in group2). The usual dense area in the middle of the plot, while less prominent here, is actually skewed a bit towards greater binding in group2. To see these two observations together is unusual.

It would be interesting to compare this to a non-normalized plot (set bNormalized=FALSE in the call to dba.plotMA). This would be to make sure that the normalization is not altering the binding characteristics too much. In RNA-seq, it is rare to see a large population of genes that are completely unexpressed in one group but expressed in the other, so normalization methods derived for RNA-seq can sometimes distort a ChIP-seq signal in ways that lead to misleading conclusions.


ADD COMMENTlink written 14 months ago by Rory Stark2.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 354 users visited in the last hour