This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Peak length for TF motif discovery

I noticed several possibilities of using the whole peak fasta sequence, 100 or 250bp up and downstream of peak midpoint for motif discovery in MEME or HOMER. I am wondering what is the standard protocol for this. The same goes for selecting background sequences - should it come from shuffling the peaks or from random positions in the genome.

chip-seq
• 1,156 views • link
updated by khorms khorms • written by Pappu Pappu

1 answer

I don't think there is a standard protocol for this, try different lengths. Regarding the background: if you do shuffling, make sure to keep the dinucleotide content (it's a very important predictor for many biological features). If you choose random positions from the genome, pick the ones where (1) there is no peak, (2) dinucleotide content is close to dinucleotide content of the peaks you have (for example, for each peak find one non-peak with similar dinucleotide content)

• 0 views • link
written by khorms khorms

...or simply run Homer with default settings to do exactly that :) Default length is 200bp in this.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this