Before coming to this forum I searched for quite lots of tutorials in the net, including code examples, videos, manuals etc, but I did not find what I was looking for and this is my last resource, so I would really appreciate if someone could give me some hints for my analysis.
My problem: I performed ChIP-seq with H3K27me3 and H3 (control) antibodies, and I have 2 conditions (wild type and mutant). I also have RNA-seq data, in which I observed downregulation of my genes of interest in the mutant.
What I would like to see if those downregulated genes in the mutant (in respect to WT) are also more enriched in H3K27me3, in comparison to WT.
I observed by using deeptools that in fact, some of my genes of interest are enriched in H3K27me3 in the mutant, but now I would like to get numbers.
So I want to 1)normalize each condition to H3, 2) Compare then mutant and WT (here I´m not looking for presence or absence of peak, but rather how much enrichment there is).
What some people I know are doing is: after mapping and QC, sorting, extending reads, they use bedtools coverage to get a table with counts. From this table they calculate RPKM value for each gene, compare conditions by setting an arbitrary log fold change threshold and with this they get calculate whether there is enrichment or not, and then correlate those genes with RNA seq data. But they don´t include normalization against H3 control.
What is your opinion about this pipeline? Do you have any advises on how to proceed with my analysis?
Thank you very much and kind regards,