Question: Using ATACseq to do motif analysis in the promoters of a gene set
gravatar for dljone02
10 months ago by
dljone020 wrote:


I have ATACseq data analyzed for my cell of interest. I also have RNAseq from my cells with/without stimulus. After doing differential gene expression from the RNAseq, I want to do motif enrichment in the promoters of this gene set (~4,000 genes). Specifically, I would like to use bedtools to intersect the DNA accessibility regions of the promoters of this gene set, to identify the regions with the highest likelihood of transcription factor binding.

Since RNA is stranded (can be on + or - strand) and ATAC is not, will there be any issue with using a standard bedtools command as follows?

bedtools intersect -a RNAseq_DEG_promoters.bed -b ATAC.bed > result.bed

Any help appreciated

rna-seq bedtools atacseq • 233 views
ADD COMMENTlink modified 10 months ago by h.mon32k • written 10 months ago by dljone020
gravatar for ATpoint
10 months ago by
ATpoint44k wrote:

Motifs are not stranded, neither is transcription factor binding a strand-specific event. They bind double stranded DNA recognizing a motif. If a sequence is overrepresented in your selection it will be overrepresented no matter if providing the top strand or its reverse complement. So no, you can use intersect as usual, only the coordinates matter. The only thing I would do though is to limit the peak to like 200bp from the peak center. THis is what motif tools such as Homer do by default to avoid motif matches by change and since the actual TF binding (and therefore the motif) is a narrow event. Depending on your peak calling you might get peaks that are notably larger than 200bp and therefore scan way more nucleotides than necessary, increasing the change of false-positives.

ADD COMMENTlink modified 10 months ago • written 10 months ago by ATpoint44k

Ok thanks. Was more or less just concerned about using "bedtools intersect" with a dataset that is stranded vs not stranded....this is ok?

ADD REPLYlink written 10 months ago by dljone020
Please log in to add an answer.


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