Hi all,
I’m analyzing ChIP-seq data and currently performing a de novo motif analysis using MEME. The output shows a high number of repeat motifs (e.g. GAGAGAGAGAGA) and only as the 10th motif I find the motif of the TF I chipped. The data is very repeat rich, but repeat masking also results in loss of my TF motif as it is within these repeats. I’m afraid that this will not enable me to identify any other co-occurring motifs of interest, so I was wondering if there is a better approach.
Is there a way to tell MEME not to recognize these simple repeats as motifs? I’ve not found out how to do this just yet.
Do you recommend another tool which is more suitable for repeat-rich sequences?
Thank you very much!
Rob
3 answers
You probably have the peak coordinates and then used something like bedtools getfasta. What you can do is to first identify the genomic coordinates of these repeats, e.g. using any of the solutions from A: Code golf: detecting homopolymers of length N in the (human) genome (modified to match these tandem patterns you encountered) and then use these coordinates to filter out any peaks that intersect with these blacklisted coordinates e.g. using bedtools intersect. Then get sequences from the remaining peaks and re-run the motif search.
Seems like dust might be a helpful tool for this. You could also browse the excellent MEME suite Q&A page or post your own question there.
Try GimmeMotifs. It combines different motif prediction tools (including MEME) and compares the identified motifs to a background set of sequences. It usually work very well for ChIP-seq data (disclaimer: I wrote the software).
Log in to answer this question.
Can't this be biologically meaningful? Which TF is this?
It sure could be biologically relevant, but I expected the number one motif to be that of the Glucocorticoid Receptor (the chipped TF).