This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Profile Plot looks strange

New to CHiP-Seq. Made a profile plot via diffbind, and it looks very different compared to the ones presented in the vignette. These are just all the samples and sites plotted. Even when I do this for a contrast it still looks odd like this. I was wondering if I did something wrong. Could someone please take a look and possibly give me pointers on how to make this look better and not have the very dark red bar at the top? If anyone could also provide guidance on how to interpret this that would be great too. Thanks. enter image description here

profile-plot chip-seq diffbind

1 answer

The dark red bar at the top is the result of a subset of sites that have a strong, continuous (non-localized) signal over the full 3000bp range. This signal is skewing the scale so that the remainder of sites are showing only faintly.

It might be worth trying to identify the sites showing this non-localized, saturated signal. I suspect an MA plot would show a cluster of sites with very high concentrations (X-axis). You can retrieve the binding matrix and look for the sites with the highest counts, then see what they look like in a browser. You may want to concentrate on the fourth sample which has the highest read concentrations. You may also want to check the read counts normalized and non-nomalized to make sure this isn't an artefact of normalization.

Log in to answer this question.