Trascription factor binding sites prediction using FIMO
0
1
Entering edit mode
5.3 years ago
Eisuan ▴ 20

Hello everyone!

In these days I was performing TF binding sites prediction using FIMO (MEME Suite). I am trying to check if there are TF binding sites for TF related to biological processes at the promoter level of the gene in my interest.

I am performing the analysis using;

  • Proximal promoter sequence (-800 / +200 from TSS)
  • promoter sequence region annotated in the ENSEMBL regulatory feature (-1600 / +1800 from TSS)
  • JASPAR 2018 PWM

In some cases I have obtained results with very high q values (like 0,50; but this was not always the case) and, in general, I know that this kind of predictions may result in many false positives. However, by contextualizing the position of the predicted binding sites (and others that are known to be near to each other) I got results that make sense.

I was wondering if there are other computational methods for increasing prediction reliability. By trying to visualize ChIP-seq experiment peaks I have seen that there were correspondences in predictions regions and peaks. The same thing for checking identity with the promoter of the orthologous gene in mouse. However, I have to admit that my approach in doing that was very naive and I am not sure it is correct.

Thank you for your support! Daniele

gene promoter • 2.4k views
ADD COMMENT
1
Entering edit mode

Is this human? If so you can intersect your promoter region with the ReMap database to see if there is experimental (ChIP-seq) evidence for binding at the promoter.

ADD REPLY
0
Entering edit mode

Thank you for this suggestion. I've just tried using that tool but unfortunately, it doesn't reveal any TF in the region of interest neither using the ENSEMBL gene ID (and modifying upstream/downstream analysis settings).

PS. Yes, it's a human gene and I am using hg38 reference genome.

ADD REPLY
1
Entering edit mode

That is unfortunate. You can always try to see if there is published ChIP-seq data towards the TFs whose motifs you found. ReMap is extensive but not comprehensive as far as I understand things. Beyond that it will probably come down to experments doing ChIP-(seq) on the factors you think that bind there.

ADD REPLY
0
Entering edit mode

It is funny because all the time I try to apply some tools that search matches in a database using this promoter region I find few / no results. I had the same experience using AME (Meme Suite) for this region.

I am pretty sure there are no problems with the promoter sequence I am using.

ADD REPLY
0
Entering edit mode

If you reproducibly get no results this may be because there are none?

ADD REPLY
1
Entering edit mode

Sounds like you've been pretty thorough to me.

ADD REPLY

Login before adding your answer.

Traffic: 1564 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6