Question: Defining ChIP-seq target genes: is there a robust cut-off?
gravatar for YanO
4.5 years ago by
United Kingdom
YanO130 wrote:

I am trying to identify the target genes of my transcription factor after ChIP-seq. I have called the peaks and I have also looked at the number of reads falling around the Transcription Start Site of each gene. I am wondering how to define if a gene is bound by my TF.

Does it make sense to do the following: Pick a cutoff value for number reads (per million) falling into (for example) +-2kb of the TSS. I've made heatmaps so I could estimate by eye the number of genes bound, and take the cutoff value from this, but this seems a bad idea. I also don't want to simply assign a peak to it's nearest gene and, provided it's close to the TSS, call the gene a target, that seems much too liberal. Similarly, looking at my heatmaps and picking e.g. 'the top 200 genes' also seems flimsy.

Is there a robust definition for whether a gene is a target or not, using ChIP-seq data?

Any advice would be most apreciated.

chip-seq • 2.0k views
ADD COMMENTlink modified 4.5 years ago by Ido Tamir5.1k • written 4.5 years ago by YanO130
gravatar for Ido Tamir
4.5 years ago by
Ido Tamir5.1k
Ido Tamir5.1k wrote:

Have a look at the GREAT tool and the paper where they outline different strategies to be more comprehensive in assigning peaks to genes. Whether you are comfortable with a more conservative approach is up to you. I guess it depends on the transcription factor and how often you find it, how robust the peaks are etc... You could even include TADs or CTCF regulatory boundaries if you are more of the adventurous type and actually have this data.

ADD COMMENTlink written 4.5 years ago by Ido Tamir5.1k
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: 1711 users visited in the last hour