Visualising the number of overlapping peaks in ChIP-seq studies
0
2
Entering edit mode
7.0 years ago
James Ashmore ★ 3.4k

I am usually asked if I can show the overlap between peaks in ChIP-seq studies using a venn diagram. However a venn diagram supposes the relationship is 1:1, where as a single peak in one sample could overlap multiple peaks in another sample. My question is how do you normally calculate the amount of peak overlap and how do you prefer to visualise/summarise this overlap?

ChIP-Seq • 6.5k views
ADD COMMENT
1
Entering edit mode

This is a common misconception that even the most seasoned analysts make.
A peak isn't a single entity. Its a region of the genome. You shouldn't use Venn Diagrams (or any other categorical visualization) for something with a beginning, end, and ideally intensity. Check out the deeptools package to see how to overlay positional information for 1 or more categories of sample/feature.

ADD REPLY
0
Entering edit mode

I haven't tried myself but venneuler for R seems like what you need: http://ieeexplore.ieee.org/xpls/icp.jsp?arnumber=5728808&tag=1

ADD REPLY
0
Entering edit mode

I have the same question as you have.

ADD REPLY
0
Entering edit mode

I have a post here about how to do this with mergePeaks in HOMER and the VennDiagram package in R, and also a post here about how to do the same with UpSet plots.

ADD REPLY

Login before adding your answer.

Traffic: 678 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6