This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding differences in motif enrichment between ATAC peak samples?

I have performed ATAC-seq on two different biological samples. I want to compare peaks between these files, as well as enrichment for particular motifs.

Firstly, I have mapped my PE reads, and then removed duplicates. I now have a bam file, that I use deeptools BamCoverage to normalise to RPKM, and the output is bedgraph. I do these for both my samples that I wish to compare.

Now, I want to compare the motifs that are enriched in the peaks between the two files. Is the following a sound methodology:

  • Run the resultant RPKM bedgraph files through Homer's findmotifsgenome.pl. Can the output knownResults.html and and homer de novo html, be compared between the two samples to check for differences in motif enrichment?
atac-seq homer
• 4,523 views • link
updated by danvoronov • danvoronov • written by a.rex a.rex

I would have to convert from bdg to bed if I wanted to use homer

• 1 views • link
written by a.rex a.rex

2 answers

Take a look at the recently released diffTF (https://difftf.readthedocs.io/en/latest/index.html) - it will allow you to infer differential transcription factor binding and motif enrichment.

• 1 views • link
written by James Ashmore James Ashmore

You can also use knownResults.txt file from the findMotifsGenome.pl so this could be used for comparisons.

Alternatively try RGT motif analysis tool for differential motif analysis, although I don't think it can do de-novo predictions. http://www.regulatory-genomics.org/motif-analysis/introduction/

Dan

• 0 views • link
written by danvoronov • danvoronov

Log in to answer this question.

More posts like this