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!
•
written
by
science_lizard •
1 answer
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)
Log in to answer this question.
-
Suggestion for enrichment analysis for small microarray set.
written by fu_entomology •
-
Graphing Average Expression of Group of Genes
written by cthangav •
-
Best normalization method to use with ChIP-seq counts
written by science_lizard •
-
Correlate ChIP-seq to gene expression
written by science_lizard •
-
NGS Plot Heatmap Colors
written by cemil.a.kerimoglu •
-
Identify genes with chip signal > input prior to peak calling
written by SciDuck •
-
Tissue-specific Histone Mark bam files?
written by epaminonda •
-
Upsampling BAMs, or downsampling by A LOT?
written by science_lizard •
-
"Good" epigenetic histone modification ChIP-seq datasets
written by Bohdan Khomtchouk
-
Can I use AME - Analysis of Motif Enrichment - for RNA-Seq data?
written by Chris Gene
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.