More posts like this
-
how to define cis regulatory elements using histone modifications ChIPseq
written by BenHu •
-
What is the difference between tissues and cell lines in ENCODE>
written by burcakotlu •
-
Getting hints on protein function by correlating its ChIP-seq peaks with public datasets?
written by bsksln •
-
Low input ChIP-seq on frozen tissue
written by evakoole •
-
regulation of protein expression
written by Pappu
-
Find transcription factors that bind to inputted nucleotide sequence?
written by randalljellis
-
Finding histone modification enrichment specific loci using modENCODE data
written by dr_bantz
-
Plot of intergenomic distances between all bound TF sites?
written by bede.portz
-
Blog: Live Webinar on 26 Aug: Calling Narrow and Broad Peaks from ChIP-Seq data
written by Genomics
-
Chip-sequencing of histone modifications
written by Affan
This question doesn't make sense as written. Can you clarify? Histones can be modified in a number of ways (a distribution of ways perhaps). And then there are transcription factors, which can be covalently modified, but aren't typically thought of in this regard. Histones and Transcription Factors and the modifications each may undergo are separate entities. What are you trying to achieve? What are you trying to figure out? What would your plot consist of? What question would it address?
Dear seidel,
Thanks for your reply. I have a bed file of H3K4me1 and a bed file of MCM2, I would like to plot the H3K4me1 enrichment on MCM2 as shown in the image below. The abscissa represents the coordinates of the MCM2. I used a simple Python script to handle it, but the running time was extremely slow because of the huge data, so I wanted to ask if there was a better way to do this thing. Thank you very much. Sincerely.
deepTools plotProfile (or plotHeatmap which also adds detailed heatmaps) can be used to plot those enrichment profiles you are referring to. (https://deeptools.readthedocs.io/en/develop/content/tools/plotProfile.html)
However, your bed files won't provide enough information for that, as they only contain the genomic location and average enrichment of (probably) peaks of H3K4me1 and MCM2 and no detailed, binned enrichment values. You would have to run deepTools computeMatrix prior to plotProfile, which requires bigWigs with the signal of your histone modification and transcription factor besides bed files specifying the genomic regions.
Hi lily, Thank you so much.
Deeptool is a really powerful tool. I'm going to learn how to use it.
Thanks again.