Question: findmotifsgenome in homer
0
gravatar for sophialovechan
2.7 years ago by
sophialovechan40 wrote:

Hi, I have some ATAC seq data along with RNA seq data. I have generated bed files of common peaks and unregulated peaks between two cell types. I was using homer to find common and upregulated motifs which have mapped upregulated genes from one cell type to another. I read the tutorial of Homer. It says that find motifgenome.pl aims at finding differential motifs compared to the background. But there was just one peak file in the script so I am not sure what the background is. Is it just the whole genome (hg19) or what ? Is it ok to use findmotifgenome. pl to find the common motifs?

sequencing • 1.8k views
ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by sophialovechan40
1
gravatar for ejm32
2.7 years ago by
ejm32440
Boston, MA
ejm32440 wrote:

Homer generates a background set of sequences based on the characteristics of the peaks you supplied (size, GC, etc).

As to finding common motifs, I'm not entirely sure what you mean... Common motifs between what exactly? A little info about your expt design and data will help with the second part of your question.

ADD COMMENTlink written 2.7 years ago by ejm32440

I have mapped upregulated genes (from RNA seq) with common peaks between two cell types and then did findmotifgenome.pl to find the motifs from that file. I have some interested motifs which I want to know which genes that they might target, so I used annotatePeak.pl to find out. The peaks score is very high and Focus Ratio/Region Size are gene names.

ADD REPLYlink written 2.7 years ago by sophialovechan40

if you are looking for de novo peaks as well, you might want to run it a few times since results may be slightly different between subsequent runs with few peaks. Just something to keep in mind.

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by steve2.0k
0
gravatar for sophialovechan
2.7 years ago by
sophialovechan40 wrote:

And another problems is that even though I can generate the motif files that have common motifs with possibly targeted up-regulated genes but when I used annotatePeaks.pl to find whether specific motifs are in those common peaks with targeted genes which are unregulated, the generated files are weird. Firstly, the peak score are very high and the row which supposed to show Focus Ratio/Region Size is the gene names rather than NA like the normal format. How can I fix that or are there any other ways that I can reach what I want to do ?

ADD COMMENTlink written 2.7 years ago by sophialovechan40

What command are you running to annotate the peaks? The "peak scores" for motifs will be different than your atac peak scores because the scoring scheme is different.

ADD REPLYlink written 2.7 years ago by ejm32440
0
gravatar for sophialovechan
2.7 years ago by
sophialovechan40 wrote:

example of my question

ADD COMMENTlink written 2.7 years ago by sophialovechan40
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: 1574 users visited in the last hour