This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to decide H3K27me3 peaks in promoter region

Hi all!

I have a H3K27me3 bed file that is predicted by SICER and H3K4me3 narrowPeak file predicted by MACS2.

I want to know what are those monovalent promoter(that is with only one mark in the promoter region)

This is what I did:

1.I use bedtools to find the H3K4me3 summits that has no overlap with H3K27me3 broad peaks(bedtools intersect -v)

2.Then I find the overlap between step 1 result and the slopped annotation bed file that have the coordinate from -500 bp and +500 bp of TSS of genes

3.Similar to step 1 & 2, I try to find the promoter that only marked by H3K27me3 broad peaks

Here is what I got: The result of monovalent promoter of H3K4me3 seems good and reasonable, while the result of H3K27me3 is not quit good. I got many results that H3K27me3 peaks spread across gene body but not in the promoter region. (I check various transcript and it seems that there are no H3K27me3 peaks that fall in promoter region).

ABCA1 ABO ADD2

Any hints about this?

Thanks a lot !

chip-seq
• 3,643 views • link
updated by Sinji Sinji • written by izzy.yichao.cai izzy.yichao.cai

1 answer

I think I can help with this, but you'll have to walk me through what you want to do since i'm having trouble fully understanding the question.

  1. You're looking for genes that contain H3K4me3, but no H3K27me3, in a -500/+500 bp region around the TSS?

  2. You're then looking for genes that contain H3K27me3, but no H3K4me3, in a -500/+500 bp region around the TSS?

You should take into account that 'promoter' region is arbitrary and depending on the literature can range from -250/+250 to -3000/+3000 from the TSS of a gene.

• 0 views • link
written by Sinji Sinji

Log in to answer this question.

More posts like this