This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to extract PRDM9 binding motifs in bed format?

I am seeking advice on relevant solution/database for extracting all PRDM9 motifs across human genome. Ideally I'd like the coordinates of start position and end position for such motifs across all human chromosomes.

chip-seq genome
• 3,594 views • link
updated by Alex Reynolds Alex Reynolds • written by reza.jabal reza.jabal

3 answers

You might look for a PRDM9 motif model matrix from one of these TF databases: http://floresta.eead.csic.es/footprintdb/?databases

Once you have the model's matrix in MEME format, you could run that through FIMO to find putative binding sites (in BED format): http://meme-suite.org/doc/fimo.html

• 1 views • link
written by Alex Reynolds Alex Reynolds

Alex thanks a lot for your suggestion. It is super helpful!

• 1 views • link
written by reza.jabal reza.jabal

EMBOSS fuzznuc might solve your issue. The output isn't exactly in bed format, but it does return start/end coordinates in a table-formatted file and also allows for fuzzy searching, which is probably handy considering the degeneracy of the PRDM9 binding motif. (Caveat being that I haven't used it myself, however)

Edit: apologies, I realize your question is also asking about a database from which these motifs could be pulled. If I'm understanding that correctly, would the fasta sequence of the human reference genome work?

• 1 views • link
written by aays aays

Thanks a lot for your suggestion. Seems like a handy package that helps solving the problem!

• 1 views • link
written by reza.jabal reza.jabal

Not entirely sure if this will help but maybe you can try using the HOMER software's scanMotifGenomeWide.pl

• 1 views • link
written by AB AB

Thank you for your suggestion, I am looking for a database that such motifs could be pulled from in bed format?

• 1 views • link
written by reza.jabal reza.jabal

Log in to answer this question.

More posts like this