Drawing ChIPseq graphics

Question

[NHEJ]

In a number of ChIPseq papers I have seen graphs that have on the y-axis plotted "peak frequency" and on the x-axis plotted "log 10 (distance to transcription start site)"

I have also seen common ChIPseq plots showing on the y-axis "reads per million mapped reads" and on the x-axis "Distance to transcription start site."

However, none of these papers mentioned what software was used to generate said plots, and I was wondering if anybody here might know this information?

It would also be beneficial if anyone had a recommendation as to what kinds of plots they frequently encounter in day to day ChIPseq papers that they consider to be almost necessary to include in any ChIPseq paper.

[seidel]

Can you provide some specific references? Your description mentions simply x,y plots, which can be drawn with virtually every graphing program ever written. Is your question more about graphing programs? Or more about software for pushing around ChIP Seq data?

Answer 1

[Madelaine Gogol]

R is pretty common, but if you don't know R, you could try ngsplot. It's an R script you call from the command line with various flags and config files. It can create a lot of the average gene type plots.

Answer 2

[Ryan Dale]

I ran into the same problem a while ago -- I saw a lot of fancy ChIP-seq figures with no information in the methods to describe how they were made.

So I wrote metaseq (paper, github, docs). At the same time, others wrote ngsplot and deepTools to address the same lack of tools.

ngsplot and deepTools are great for use from the command line; metaseq is more of a toolbox. Example 1 walks you through making a custom ChIP-seq figure and gives a taste for what you can do with the package.

Answer 3

[Ming Tommy Tang]

You need to know some basic programming skills.

Annotate your peaks first ( how far is the peak from the TSSs), then you count the peak numbers in each different categories(0-5kb, 5-10kb etc) and then plot a bar graph.

Easier way, try to use GREAT tools

Answer 4

[madk00k]

A lot of graphs of this type can be created using CEAS (http://liulab.dfci.harvard.edu/CEAS/)

CEAS is also a part of Cistrome pipeline, which can be accessed from Galaxy web interface (http://cistrome.org/ap/) or from desktop application Unipro UGENE (http://ugene.unipro.ru)

Answer 5

[Guangchuang Yu]

ChIPseeker can creates similar figures: http://ygc.name/2014/04/30/visualization-methods-in-chipseeker/

Answer 6

[Ming Tommy Tang]

I prefer Genomation.

High quality output and very flexible and easy to use! If you know some R.