Question: Finding differences in motif enrichment between ATAC peak samples?
1
gravatar for a.rex
4 months ago by
a.rex170
a.rex170 wrote:

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 • 359 views
ADD COMMENTlink modified 4 months ago by James Ashmore2.6k • written 4 months ago by a.rex170

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

ADD REPLYlink written 4 months ago by a.rex170
0
gravatar for James Ashmore
4 months ago by
James Ashmore2.6k
UK/Edinburgh/MRC Centre for Regenerative Medicine
James Ashmore2.6k wrote:

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.

ADD COMMENTlink written 4 months ago by James Ashmore2.6k
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: 836 users visited in the last hour