Question: Differential accessibility analysis with Genrich peaks
gravatar for eyleen.corrales
4 months ago by
eyleen.corrales20 wrote:


I'm attempting to do differential accessibility analysis on ATAC-seq paired-end data. I used Generich for peak calling, using ATAC-seq mode, which calls the peaks on the Tn5 cut sites. Afterwards, I tried to count the reads from my bam files within the regions of the consensus peaks and proceed to do differential analysis with edgeR. I have the doubt of how accurate this could be, as the intervals used for peak calling would be centered at the ends of the reads, so I guess the called peaks would not fully match with the pileup of reads from my bam files. Any ideas on how to address this? Should I be shifting the reads from my bam files?

Here is a small peace of code I'm using for counting reads from my bam files (bamFiles) in the peak ranges called by Genrich (peakList), in case that it's of any help:

consensus <- soGGi:::runConsensusRegions(GRangesList(peakList), method = "none")

blacklist <- rtracklayer::import("/ATAC_blacklists/ENCFF001TDO.hg19.bed")
standard.chr <- paste0("chr", c(1:22, "X", "Y"))
param <- readParam(max.frag=1000, pe="both", discard=blacklist, restrict=standard.chr)

peak.counts <- regionCounts(bamFiles, consensus, param=param)
genrich atac-seq • 247 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by eyleen.corrales20

I have not looked at the details of your code at all, but have you checked out the DiffBind package that provides a wrapper for DESeq2 and edgeR based analyses for ATAC-seq/ChIP-seq samples?

ADD REPLYlink written 4 months ago by Friederike6.5k

Thank you for your reply. I haven't tried it yet, but I think my concern would be the same. Do you normally use DiffBind to analyse differential accessibility of peaks called by Generich on ATAC-seq mode, without the need to shifting your bam file reads to the Tn5 cut site? Maybe I'm just worried about something that is not a real problem

ADD REPLYlink written 4 months ago by eyleen.corrales20
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: 1060 users visited in the last hour