ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide Buenrostro et al. 2015 Methodology
?
Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data Cheng et al. 2024 Research Article
?
Analytical Approaches for ATAC-seq Data Analysis Smith et al. 2021 Research Article
"Here, we explain the fundamentals of ATAC-seq data processing, summarize common analysis approaches, and review computational tools to provide recommendations for different research questions."

Generic ATAC-seq Analysis Workflow

General Data Processing
Convert Format 1.
?
Raw Read Quality 1.5
?
Trim Adapters 2.
?
Align Reads 3.
?
Remove Mitochondrial Reads 4.
?
Deduplicate 5.
?
Signal Tracks 6.
?
Call Peaks 7.
?
QC Plots 7a.
?
Summary Statistics 7b.
?
Project Specific Analysis
Motif Enhancement 9.
?
Differentially Accessible Regions 9.
?
Region Enrichment 9.
?
Nucleosome Positioning 9.
?

Guides & Tools

Table 1: Step by Step Guides to Performing ATAC-seq Data Analysis
Title and authorNotes and linkLast Updated
ATAC-seq data analysis: from FASTQ to peaks
Yiwei Niu
Blog style walkthrough of generalized ATAC-seq data analysis.2019
BIOINF525 Lab 3.2
Steve Parker
Minimal standard ATAC-seq analysis walkthrough.2016
Analysis of ATAC-seq data in R and Bioconductor
Rockefeller Bioinformatics Resource
Bioconductor ATAC-seq analysis course.2018
ATAC-seq
John M. Gaspar
Generalized ATAC-seq analysis walkthrough with included custom scripts.2019
ATAC-seq data analysis
Delisle L; Doyle M; & Heyl F
Galaxy training walkthrough of generalized ATAC-seq analysis.2020
Analytical Approaches for ATAC-seq Data Analysis
Table 2: Raw ATAC-seq Data Processing Pipelines
PipelineLanguageNotesDocsCitation
AIAPBash; R; PythonOptimized analysis with novel QC metrics++Liu et al. (2019)
Last updated: 2019
ATAC2GRNBash; PythonParameter optimized ATAC-seq pipeline+Pranzatelli, Michael, & Chiorini (2018)
Last updated: 2018
ATAC-pipePython; RAnalysis pipeline for ATAC-seq data including TF footprinting; cell-type classification; and regulatory network creation+++Zuo et al. (2019)
Last updated: 2019
ATACProcBash; Python; RComplete pipeline with additional downstream analyses included++Unpublished
Last updated: 2019
BasepairNACommercial. Web-based GUI for complete analysis?Unpublished
CIPHERR; Perl; PythonA data processing platform for ChIP-seq; RNA-seq; MNase-seq; DNase-seq; ATAC-seq; and GRO-seq datasets+Guzman & D’Orso (2017)
Last updated: 2017
ENCODEPython; BashComplete pipeline following ENCODE standards for ATAC/DNase-seq analysis++Unpublished
Last updated: 2020
esATACRComplete pipeline including downstream analyses+++Wei, Zhang, Fang, Li, & Wang (2018)
Last updated: 2019
GUAVAJava; Python; RGUI based complete ATAC-seq pipeline+Divate & Cheung (2018)
Last updated: 2019
I-ATACJavaGUI based interactive ATAC-seq pipeline+Ahmed & Ucar (2017)
Last updated: 2017
nfcore/atacseqPython; RComplete pipeline build using Nextflow+++Ewels et al. (2019)
Last updated: 2019
PEPATACPython; R; PerlComplete pipeline with unique analytical approaches and QC metrics+++Unpublished
Last updated: 2019
pyflow-ATACseqBash; PythonATAC-seq snakemake pipeline with included nucleosome positioning and TF footprinting++Unpublished
Last updated: 2019
snakePipes ATAC-seqPythonWorkflow system including ATAC-seq analysis+++Bhardwaj et al. (2019)
Last updated: 2019
Tobias RauschBash; R; PythonComplete pipeline with emphasis on downstream analyses++Rausch et al. (2019)
Last updated: 2020
Analytical Approaches for ATAC-seq Data Analysis
Table 3: ATAC-seq Advanced Quality Control Metric Tools
ToolLanguagesNotesDocsCitation
ATAqCBash; PythonGenerate ATAC-seq specific quality control metrics.+Unpublished
Last updated: 2017
ATACseqQCRProvides ATAC-seq specific quality control metrics and transcription factor footprinting.+++Ou et al. (2018)
Last updated: 2018
ataqvC++; BashATAC-seq QC and visualization.+++Orchard, Kyono, Hensley, Kitzman, & Parker (2020)
Last updated: 2020
Analytical Approaches for ATAC-seq Data Analysis
Table 4: Peak Calling Tools
ToolLanguagesNotesDocsCitation
F-SeqJavaCan be used as general peak caller to identify regions of open chromatin.++Boyle et al. (2008)
Last updated: 2016
GenrichCPeak caller for genomic enrichment assays with specific ATAC-seq mode.+++unpublished
Last updated: 2020
HMMRATACJavaIdentify nucleosome positioning and leverage ATAC-seq specific read outs to call peaks.+++Tarbell & Liu (2019)
Last updated: 2020
Hotspot2C++Identify significantly enriched genomic regions.++Unpublished
Last updated: 2019
HOMERPerl; C++Suite of tools that include the ability to call peaks from DNA enrichment assays.+++Heinz et al. (2010)
MACS2PythonSpecifically designed for ChIP-seq but broadly applicable to any DNA enrichment assay to call peaks.+++Zhang et al. (2020)
Last updated: 2020
PeaKDEckPerlPeak calling program for DNase-seq data.+++McCarthy & O’Callaghan (2014)
Last updated: 2014
Analytical Approaches for ATAC-seq Data Analysis
Table 5: Tools to Investigate Differentially Accessible Regions
ToolLanguagesNotesDocsCitation
DAStkPythonIdentifies changes in transcription factor activity by looking at changes in chromatin accessibility+++Tripodi et al. (2018)
Last updated: 2020
diffTFPython; RIdentifies differential transcription factors. Can operate in basic mode with just chromatin accessibility or in classification mode where it integrates RNA-seq.+++Berest et al. (2019)
Last updated: 2020
Analytical Approaches for ATAC-seq Data Analysis
Table 6: Motif Enrichment and Transcription Factor Footprinting Tools
ToolLanguagesNotesDocsCitation
BiFETRIdentify overrepresented transcription factor footprints.++Youn et al. (2019)
Last updated: 2019
BinDNaseRTranscription factor binding prediction using DNase-seq.+Kähärä & Lähdesmäki (2015)
Last updated: 2015
CENTIPEDERTranscription factor footprinting and binding site prediction.++Pique-Regi et al. (2011)
Last updated: 2010
DeFCoMPythonDetecting transcription factor footprints and underlying motifs using supervised learning.+++Quach & Furey (2017)
Last updated: 2017
DNase2TFRIdentify footprint candidates from DNase-seq data on user-specified regions.+Sung et al. (2014)
Last updated: 2017
HINT-ATACPythonUse open chromatin data to identify transcription factor footprints with modifications specific to ATAC-seq data.+++Li et al. (2019)
Last updated: 2019
HOMERPerl; C++A suite of tools for motif discovery and enrichment.+++Heinz et al. (2010)
Last updated: 2019
MEME SuitePerl; PythonSuite of tools for motif discovery; enrichment; and GO term analyses.+++Bailey et al. (2009)
Last updated: 2020
PIQBash; RModels genome-wide DNase profiles to identify transcription factor binding sites.++Sherwood et al. (2014)
Last updated: 2016
TOBIASPythonIdentify transcription factor footprints.++Bentsen et al. (2019)
Last updated: 2020
TRACEPythonTranscription factor footprinting.++Ouyang & Boyle (2019)
Last updated: 2020
WellingtonPythonIdentify TF footprints using DNase-seq data.+++Piper et al. (2013)
Last updated: 2019
Analytical Approaches for ATAC-seq Data Analysis
Table 7: Tools to Investigate Nucleosome Positioning
ToolLanguagesNotesDocsCitation
HMMRATACJavaIdentify nucleosome positioning and leverage ATAC-seq specific read outs to call peaks.+++Tarbell & Liu (2019)
Last updated: 2020
NucleoATACPython; RCall nucleosomes using ATAC-seq data.+++Schep et al. (2015)
Last updated: 2019
NucToolsPerl; RCalculate nucleosome occupancy profiles on chromatin accessibility data.+++Vainshtein et al. (2017)
Last updated: 2019
Analytical Approaches for ATAC-seq Data Analysis
Table 8: Tools to Investigate Region Enrichment
ToolLanguagesNotesDocsCitation
AnnotatrRAnnotate summarize and visualize genomic regions.+++Cavalcante & Sartor (2017)
Last updated: 2019
BART/BARTwebPythonPredict factors that bind at cis-regulatory regions.+++Wang et al. (2018)
Last updated: 2020
chipenrichRPerform gene set enrichment testing using genomic regions.+++Welch et al. (2014)
Last updated: 2020
coloc-statsPythonPerform co-localization analysis of genomic regions.+++Simovski et al. (2018)
Last updated: 2019
COLOJSPIdentify genomic features in close proximity to user-submitted genomic regions.++Kim et al. (2015)
Last updated: 2015
FEATnotatorPerl; RAnnotate genomic regions.++Podicheti & Mockaitis (2015)
Last updated: 2018
GenomeRunner.NETPerform annotation and enrichment of genomic regions against default or custom regulatory regions.++Dozmorov et al. (2016)
Last updated: 2016
GenometriCorrRDetermine spatial correlation between region sets.++Favorov et al. (2012)
Last updated: 2020
Genomic Association TesterPythonCalculate the significance of overlaps between multiple genomic region sets.+++Heger et al. (2013)
Last updated: 2019
GIGGLECGenomics search engine to uncover significantly shared genomic loci (regions) between data.+++Layer et al. (2018)
Last updated: 2019
GLANETJava; PerlGenomic loci annotation and enrichment tool between sets of genomic regions.+++Otlu et al. (2017)
Last updated: 2019
GREATCAnnotate genomic regions.+++McLean et al. (2010)
Last updated: 2019
LOLA/LOLAwebRDetermine significant enrichment between region sets to inform on biological meaning.+++Sheffield & Bock (2016)
Last updated: 2019
regioneRREvaluate significant associations between region sets using permutation testing.+++Gel et al. (2016)
Last updated: 2020
StereoGeneC++; REstimate genome-wide correlation between pairs of genomic features.++Stavrovskaya et al. (2017)
Last updated: 2019
Analytical Approaches for ATAC-seq Data Analysis
Table 9: Available Tools for Single-Cell ATAC-seq Data Processing
ToolLanguagesNotesDocsCitation
BAPR; PythonBead-based scATAC-seq data processing.++Lareau et al. (2019)
Last updated: 2019
BROCKMANR; Bash; RubyConvert genomics data into K-mer words associated with chromatin marks used to compare and identify changes across samples.++de Boer & Regev (2018)
Last updated: 2018
Cell Ranger ATACNACommercial. Set of analysis pipelines for Chromium single cell ATAC-seq.+++Unpublished
chromVARRIdentify transcription factor accessibility in single-cell data. Enables clustering of single-cell ATAC-seq data.+++Schep et al. (2017)
Last updated: 2019
CiceroRPredict cis-regulatory DNA interactions using single-cell chromatin accessibility data.+++Pliner et al. (2018)
Last updated: 2019
cisTopicRIdentify cell states and cis-regulatory topics from single-cell data.+++Bravo González-Blas et al. (2019)
Last updated: 2019
scABCRClassify single-cell ATAC using unsupervised clustering and identify chromatin regions specific to cell identity.+Zamanighomi et al. (2018)
Last updated: 2019
SCALEPythonClustering and visualization of single-cell ATAC-seq data into interpretable cell populations.++Xiong et al. (2019)
Last updated: 2019
ScasatBash; Python; RComplete pipeline to process scATAC-seq data with simple steps.+++Baker et al. (2019)
Last updated: 2019
scATAC-proR; PythonComprehensive pipeline for single cell ATAC-seq analysis.+++Yu et al. (2019)
Last updated: 2020
scOpenPythonChromatin-accessibility estimation of single-cell ATAC data.+Li et al. (2019)
Last updated: 2020
SCRATRUseful for studying single cell heterogeneity. Can identify changes in gene sets or transcription factor binding sites. Includes GUI and web-based service.+++Ji et al. (2017)
Last updated: 2018
SnapATACR; PythonSingle Nucleus Analysis Pipeline for ATAC-seq.+++Fang et al. (2019)
Last updated: 2019
Analytical Approaches for ATAC-seq Data Analysis