Question: ChIP-seq enrichment over specific genes
gravatar for science_lizard
13 months ago by
science_lizard0 wrote:

I'm wondering if anyone knows a good program, perhaps on Bioconductor, that will help determine ChIP-seq enrichment over particular genes (would be wonderful if they would also be able to bin that enrichment over gene compartments like promotors, introns, etc.). For example, if I knock down a histone regulator, I'd like to show over specific genes that the mark is reduced in treatment vs control samples. Of course, I'd be happy to do this manually (i.e., making a box plot in R after calling ChIP-seq peaks, though I'm not sure what value would be best to plot), but Bioconductor packages do tend to be a bit more user-friendly for programming newbies like myself. Thanks in advance!

chip-seq • 376 views
ADD COMMENTlink modified 13 months ago by Charles Warden7.6k • written 13 months ago by science_lizard0

Why don't you perform a standard differential analysis, e.g. with DiffBind or csaw, and then check if the histone mark is significantly depleted at these genes (or its promoters). To illustrate this, I would probably plot the log2 read counts or the enrichments over input over these regions as a boxplot, one box per condition.

ADD REPLYlink written 13 months ago by ATpoint29k
gravatar for Charles Warden
13 months ago by
Charles Warden7.6k
Duarte, CA
Charles Warden7.6k wrote:

If you have gene mappings, there may still be use for relatively standard gene enrichment tools (Enrichr, GATHER, DAVID, etc.)

Strictly speaking, I believe these programs are designed for ChIP-Seq enrichment:


ChIP-Enrich: web-interface or Bioconductor package

However, if downstream targets are included among "sets of genes," I think the motif enrichment programs could also be useful:



MEME Suite - includes GOMo for GO Motif analysis

Cistrome (SeqPos)

ADD COMMENTlink modified 13 months ago • written 13 months ago by Charles Warden7.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.3.0
Traffic: 1125 users visited in the last hour