Question: Quantification of acetylation enriched regions at transcription factor binding sites (WT vs KO)
gravatar for massimofaggiani
6 weeks ago by
massimofaggiani0 wrote:

I am still at the beginning of ChIP-seq data analysis and am stuck in the following problem: I want to quantify the enriched peaks of H3K27ac at the proximity (5kb) of each binding site of a transcription factor under study, and make a comparison with the knock-out of the same TF. I have already calculated the total number of H3K27ac peaks that fall in the 5kb range between the 2 conditions but not for each binding site. How could I accomplish this, given also the total number of TF binding sites? And then, how could I represent it graphically? Thank you.

chip-seq • 167 views
ADD COMMENTlink modified 6 weeks ago by Prakash2.0k • written 6 weeks ago by massimofaggiani0
gravatar for Prakash
6 weeks ago by
Prakash2.0k wrote:

One suggestion would be to quantify the H3K27ac enrichment ± 5kb (or ± 2kb) relative to peak summit for the given set of identified TF binding sites of your interest. This can be performed in HOMER using option e.g -size 5000 -hist 10 . More detail can be found here . Further graphically it can be represented in R using basic plot function or ggplot2

ADD COMMENTlink modified 5 weeks ago • written 6 weeks ago by Prakash2.0k

Thank you for your answer, but I managed to do that simply by adding -c (count) option in bedtools window command. Now I would like to filter only the regions for which there is an increase in acetylation peaks in the TF knock-out with respect to wild-type, having two separate bed files. Any suggestions?

ADD REPLYlink written 5 weeks ago by massimofaggiani0
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: 1179 users visited in the last hour