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

Hi,

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 • 142 views
ADD COMMENTlink modified 3 months ago by h.mon29k • written 3 months ago by dljone020
0
gravatar for ATpoint
3 months ago by
ATpoint35k
Germany
ATpoint35k 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 3 months ago • written 3 months ago by ATpoint35k

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 3 months ago by dljone020
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: 1295 users visited in the last hour