Question: Motif location for enriched motif in HOMER
0
gravatar for Prakash
8 months ago by
Prakash1.0k
India
Prakash1.0k wrote:

Hi,

A quick question, why homer gives different number of regions for a given motif when extracting motif location using findMotifsGenome.pl and annotatePeaks.pl

e.g from HOMER tutorial

findMotifsGenome.pl ERalpha.peaks hg18 MotifOutputDirectory/ -find motif1.motif > outputfile.txt

annotatePeaks.pl ERalpha.peaks hg18 -m motif1.motif > outputfile.txt

If I analyze for my peaks for a given motif , I get different no. of total regions having this motif.

chip-seq next-gen • 508 views
ADD COMMENTlink modified 8 months ago by genomax67k • written 8 months ago by Prakash1.0k
2
gravatar for Kevin Blighe
8 months ago by
Kevin Blighe42k
Republic of Ireland
Kevin Blighe42k wrote:

Just on face value, the programs have different functionality, so, I would not expect them to produce the exact same results.

annotatePeaks.pl is about annotating already-identified peaks with current knowledge.

findMotifsGenome.pl is specifically about "discovering motifs in ... genomic regions". One clear distinction that I see is that findMotifsGenome.pl specifically states that it will search for novel motifs.

There are undoubtedly other subtle differences in the command line parameter settings.

Kevin

ADD COMMENTlink written 8 months ago by Kevin Blighe42k

Thanks Kevin !

Actually I performed de novo motif analysis of peaks using findMotifsGenome.pl and it showed enrichment for a motif, say 20 % of total target. I was interested to know the genomic regions for this enriched motif. So I used findMotifsGenome.pl to extract the regions and it extracted genomic regions more than 20% whereas annotatePeaks.pl extracted ~20 % of total peaks.

ADD REPLYlink written 8 months ago by Prakash1.0k
1

Perhaps you could look at some individual examples of where they differ. Then, you may better understand in which way(s) the 2 functions behave differently

ADD REPLYlink written 8 months ago by Kevin Blighe42k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1116 users visited in the last hour