ATAC-seq: Principles, Methods, and Applications in Genomics

By Dr. Zubair Khalid, DVM, MS, PhD · 2026-08-11

Key Takeaways

• ATAC-seq leverages a hyperactive Tn5 transposase to simultaneously fragment and tag nucleosome-depleted, accessible chromatin regions with sequencing adapters, enabling genome-wide mapping of regulatory elements.
• The assay is distinguished by its low input requirement (500–5,000 cells), rapid protocol (<3 hours), and high signal-to-noise ratio compared to predecessors like DNase-seq and MNase-seq.
• Key experimental parameters include optimizing the transposase-to-nuclei ratio, controlling tagmentation time (30 minutes at 37°C is standard), and determining the appropriate number of PCR cycles to avoid over-amplification and duplicate reads.
• Computational analysis necessitates applying the Tn5 shift correction ( +4 bp on the positive strand, -5 bp on the negative strand) and filtering mitochondrial reads, which can comprise a significant portion of the sequencing output.
• Quality control metrics such as fragment size distribution (showing periodicity at ~100 bp, ~200 bp, and ~400 bp) and Transcription Start Site (TSS) enrichment score are critical for assessing library quality.
• Single-cell ATAC-seq (scATAC-seq) extends the methodology to individual cells, facilitating the analysis of chromatin accessibility heterogeneity and the identification of rare cell populations within complex tissues.

Introduction to ATAC-seq

What is ATAC-seq?

Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is a high-throughput method for mapping genome-wide chromatin accessibility. The technique leverages a hyperactive Tn5 transposase to simultaneously fragment and tag accessible regions of chromatin with sequencing adapters, thereby identifying genomic loci that are nucleosome-depleted and bound by regulatory proteins. Because accessible chromatin marks active promoters, enhancers, insulators, and other regulatory elements, ATAC-seq provides a genome-wide readout of the regulatory landscape of a cell.

The method was first described in 2013 and rapidly became the method of choice for chromatin accessibility profiling due to its simplicity, speed, and low input requirements. Unlike older approaches that require millions of cells, ATAC-seq can generate high-quality profiles from as few as 500–5,000 cells, and even single cells when appropriately scaled. The entire experimental protocol, from cells to sequencing-ready libraries, can be completed in under 3 hours, making it the fastest chromatin accessibility assay available.

Comparison with other chromatin accessibility assays

ATAC-seq is one of several methods for assessing chromatin structure genome-wide. The two principal predecessors are DNase-seq and MNase-seq, each with distinct molecular mechanisms and biases.

DNase-seq uses the endonuclease DNase I to preferentially cleave accessible chromatin. After digestion, small DNase I-hypersensitive fragments are size-selected, adapter-ligated, and sequenced. The method provides high-resolution maps of open chromatin but requires 1–50 million cells and multiple days of processing. The enzyme's sequence preference for certain dinucleotides also introduces bias that must be computationally corrected.

MNase-seq (micrococcal nuclease digestion followed by sequencing) maps nucleosome positions rather than accessibility per se. MNase digests linker DNA between nucleosomes, leaving mononucleosomal fragments that are sequenced to determine nucleosome occupancy. The method reports the positions of nucleosomes themselves, which is complementary to, but distinct from, accessibility mapping. MNase-seq typically requires 1–10 million cells.

ATAC-seq differs from both in that the Tn5 transposase performs fragmentation and adapter insertion in a single enzymatic step. This eliminates separate ligation reactions and reduces sample loss. The transposase also has a strong preference for open chromatin, integrating into nucleosome-free regions with high specificity. The result is a method with superior signal-to-noise ratio, lower input requirements, and shorter protocol time compared to DNase-seq and MNase-seq.

Mechanism of ATAC-seq

Tn5 transposase and the transposome complex

The molecular engine of ATAC-seq is Tn5 transposase, a bacterial enzyme originally isolated from Thermus aquaticus. In nature, Tn5 catalyzes the movement of transposable elements through a "cut-and-paste" mechanism. The hyperactive variant used in ATAC-seq contains two point mutations (E54K and L372P) that increase its catalytic activity and reduce its sequence specificity, making it suitable for genome-wide applications.

The active form of the enzyme is a transposome complex: a dimer of Tn5 transposase bound to two mosaic end (ME) DNA sequences. In the commercial preparation used for ATAC-seq, the ME sequences are pre-loaded with Illumina sequencing adapters. The transposome is assembled by incubating the hyperactive Tn5 enzyme with adapter-containing ME duplexes at room temperature for 30 minutes, after which the complex is stable and catalytically active.

The Tn5 transposase exhibits a strong structural preference for nucleosome-free DNA. This preference arises from the physical constraint imposed by nucleosomes: the DNA wrapped around a histone octamer is sterically inaccessible to the transposase's DNA-binding domain. In contrast, linker DNA and nucleosome-depleted regions present a relatively unobstructed substrate. The enzyme's affinity for accessible DNA is further enhanced by its ability to recognize DNA structural features, such as minor groove width, that correlate with nucleosome exclusion.

Insertion bias and tagmentation process

The tagmentation reaction—a portmanteau of "tagmentation" and "fragmentation"—proceeds through a well-characterized series of steps. The transposome binds to accessible DNA and introduces a double-strand break with a 9-base-pair staggered cut. Simultaneously, the adapter sequences carried by the transposase are covalently joined to the 5′ ends of the fragmented DNA. The result is a population of DNA fragments that are immediately compatible with PCR amplification, without the need for separate end-repair, A-tailing, or adapter ligation steps.

The transposition reaction is distributive: each transposome molecule catalyzes one insertion event and then remains bound to the product DNA. This property is exploited in the standard ATAC-seq protocol, where the reaction is stopped by adding a chaotropic salt (typically EDTA) that denatures the transposase and releases the fragmented DNA.

Tn5 exhibits a modest sequence insertion bias, favoring integration into regions with a weak dinucleotide periodicity. This bias is largely corrected by the Tn5 enzyme's preference for DNA with a narrow minor groove, which is a feature of nucleosome-free regions. However, the bias is not negligible, and computational methods such as Tn5 shifting (adjusting read positions to account for the 9-bp duplication) and peak-calling algorithms that model the insertion distribution are used to mitigate its effects.

The fragment size distribution produced by tagmentation is informative. Because Tn5 integrates into nucleosome-free regions, the resulting fragments correspond to the distance between two accessible sites. Short fragments (<100 bp) represent nucleosome-free regions, while longer fragments (180–247 bp, 315–400 bp, and larger) correspond to DNA wrapped around one, two, or more nucleosomes. This size distribution allows ATAC-seq data to be used not only for accessibility mapping but also for inferring nucleosome positioning.

Experimental Workflow

Cell lysis and nuclei preparation

The ATAC-seq protocol begins with the isolation of intact nuclei. The quality of this step is critical: over-lysis releases genomic DNA that is then accessible to Tn5, producing high background; under-lysis leaves cytoplasmic components that can inhibit the transposition reaction.

For cultured cells, the standard protocol uses a hypotonic lysis buffer containing 10 mM Tris-HCl (pH 7.4), 10 mM NaCl, 3 mM MgCl₂, and 0.1% (v/v) IGEPAL CA-630. Cells are resuspended in this buffer at a concentration of approximately 50,000 cells per 50 µL and incubated on ice for 10 minutes. The lysis is then quenched by adding 1 mL of wash buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl₂) and centrifuging at 500 × g for 5 minutes at 4°C. The supernatant is discarded, and the nuclei pellet is resuspended in the transposition reaction mix.

For tissue samples, a more involved dissociation step is required. Fresh tissue is minced and digested with collagenase or a similar enzyme cocktail, followed by filtration through a 40-µm cell strainer. The resulting single-cell suspension is then subjected to the same lysis protocol. Importantly, ATAC-seq requires fresh or cryopreserved cells; fixed cells are not compatible because cross-linking prevents Tn5 access to chromatin.

The optimal nuclei number depends on the application. For bulk ATAC-seq, 50,000 cells is the standard starting amount, though the protocol can be scaled down to 500 cells with appropriate adjustments to reaction volumes. The key is to maintain a constant ratio of transposase to nuclei; over-titration leads to excessive fragmentation and short fragments, while under-titration produces large fragments and low complexity.

Tagmentation reaction

The tagmentation reaction is assembled in a 50 µL volume containing the isolated nuclei, 25 µL of 2× TD buffer (20 mM Tris-acetate pH 7.6, 10 mM Mg-acetate, 20% dimethylformamide), and 2.5 µL of Tn5 transposase (100 nM final concentration). The reaction is incubated at 37°C for 30 minutes in a thermocycler.

The temperature and duration of the tagmentation reaction are the primary variables controlling fragment size. Shorter incubation times (10–15 minutes) produce larger fragments and are preferred when starting with low cell numbers, as they reduce the extent of over-fragmentation. Conversely, longer incubation times (up to 60 minutes) generate a higher proportion of short fragments, which can improve the signal-to-noise ratio for accessibility mapping but at the cost of reduced library complexity.

After tagmentation, the reaction is stopped by adding 5 µL of 0.5 M EDTA and incubating at 65°C for 10 minutes. This step denatures the transposase and releases the fragmented DNA. The DNA is then purified using a commercial spin column (e.g., Qiagen MinElute) and eluted in 20 µL of elution buffer (10 mM Tris-HCl pH 8.0).

Library amplification and purification

The purified tagmented DNA is amplified by PCR to add the full adapter sequences and barcode indices. The PCR reaction uses primers that anneal to the partial adapter sequences introduced by Tn5. A typical reaction contains 20 µL of purified DNA, 2.5 µL of a 25 µM stock of each primer (forward and reverse), and 25 µL of 2× NEBNext High-Fidelity PCR Master Mix, in a final volume of 50 µL.

The PCR cycling conditions are: 72°C for 5 minutes (to extend the adapter sequences), 98°C for 30 seconds (initial denaturation), followed by 5–15 cycles of 98°C for 10 seconds, 63°C for 30 seconds, and 72°C for 1 minute. A final extension at 72°C for 5 minutes completes the reaction.

The number of PCR cycles is a critical parameter. Because Tn5 inserts adapters at both ends of each fragment, the library contains a heterogeneous population of molecules with varying adapter configurations. Only fragments with adapters at both ends are amplifiable. The optimal cycle number is determined empirically by qPCR: a 5 µL aliquot of the PCR reaction is monitored in a real-time instrument, and the reaction is stopped during the exponential phase, typically when the fluorescence reaches approximately one-third of the plateau. This typically corresponds to 8–12 cycles for 50,000 cells and 12–15 cycles for 500 cells. Over-amplification produces duplicate reads and reduces library complexity.

Following PCR, the amplified library is purified using AMPure XP beads (Beckman Coulter) at a 1.2× bead-to-sample ratio to remove primers and short fragments. A second size selection can be performed to enrich for fragments below 400 bp, which represent nucleosome-free and mononucleosomal regions. This is achieved by adding an additional 0.4× volume of beads and retaining the supernatant. The final library is quantified by Qubit fluorometry and assessed for quality on a Bioanalyzer or TapeStation, where a successful library shows a characteristic fragment size distribution with a prominent peak below 100 bp.

Data Analysis Pipeline

Preprocessing and QC

The raw sequencing data from ATAC-seq are processed through a bioinformatics pipeline that shares many steps with other genomic assays. The first step is quality control using FastQC to assess per-base quality scores, GC content, adapter contamination, and duplication rates. Low-quality bases and adapter sequences are removed with Trimmomatic or cutadapt.

A unique feature of ATAC-seq data is the presence of mitochondrial reads, which can constitute 30–80% of all reads depending on cell type and lysis efficiency. These reads are typically removed by aligning to the mitochondrial genome and filtering them out, or by aligning to a combined reference that includes the mitochondrial chromosome and then excluding those alignments from downstream analysis.

Another QC metric specific to ATAC-seq is the fragment size distribution. The insertion site of Tn5 creates a characteristic periodicity in the fragment sizes, with peaks at approximately 100 bp (nucleosome-free), 200 bp (mononucleosome), 400 bp (dinucleosome), and 600 bp (trinucleosome). A successful ATAC-seq library should show this periodic pattern. The ratio of reads in nucleosome-free regions to reads in nucleosome-occupied regions (the "NFR score") is a useful quality metric; a high NFR score indicates good signal.

The transcription start site (TSS) enrichment score is another important QC metric. This score measures the ratio of reads at annotated TSSs (±50 bp) to reads in flanking regions (±1–2 kb). A high TSS enrichment score (typically >6 for high-quality data) indicates that the library is enriched for accessible chromatin at promoters, which is the expected pattern.

Alignment and peak calling

Reads are aligned to the reference genome using a splice-aware aligner such as Bowtie2 or BWA-MEM. For ATAC-seq, the alignment parameters should be adjusted to allow for the short fragment sizes typical of the assay. After alignment, reads mapping to the mitochondrial genome, reads with low mapping quality (MAPQ < 30), and PCR duplicates are removed.

A critical step in ATAC-seq data processing is the Tn5 shift. Because Tn5 integrates as a dimer and creates a 9-bp staggered cut, the actual binding site of the transposase is offset from the read start position. For reads aligned to the positive strand, the read start position is shifted +4 bp; for reads on the negative strand, the read start is shifted −5 bp. This correction centers the reads on the actual insertion site and is essential for high-resolution peak calling and footprinting analysis.

Peak calling is performed using MACS2, which was originally developed for ChIP-seq but is widely used for ATAC-seq. The key parameters are __MASK_1 (to disable the shifting model, since the Tn5 shift has already been applied), MASK_2 and MASK_3 (to extend reads to a uniform 200-bp fragment size for peak detection), and MASK_4 (or MASK_5 for more stringent calling). The MASK_6__ option can be used to call broad domains of accessibility, which are characteristic of enhancer regions.

For differential accessibility analysis between conditions, tools such as DESeq2 or edgeR are used on the count matrix of reads overlapping peaks. The counts are normalized using the "library size" method (total reads in peaks) rather than total reads, because the latter is dominated by mitochondrial contamination and non-specific background.

Downstream analysis: motif and footprinting

Once peaks are called, the next step is to annotate them with respect to genomic features. Peaks are classified as promoter-proximal (within ±2 kb of an annotated TSS), intronic, exonic, or intergenic. The distribution of peaks across these categories provides insight into the regulatory architecture of the cell type under study.

Motif analysis identifies transcription factor binding sites within peaks. The HOMER software suite is commonly used for de novo motif discovery, which finds sequence motifs enriched in peaks relative to genomic background. Known motif analysis, using databases such as JASPAR or TRANSFAC, then matches the discovered motifs to specific transcription factors. This analysis reveals which transcription factors are likely to be active in the cell type, based on the accessibility of their binding sites.

Footprinting analysis provides a higher-resolution view of transcription factor occupancy. Because Tn5 cannot integrate into DNA that is tightly bound by a protein, the insertion frequency at transcription factor binding sites is reduced, creating a "footprint" of protected DNA flanked by regions of high accessibility. Tools such as HINT-ATAC or TOBIAS use the Tn5 insertion profile to detect these footprints and infer transcription factor occupancy. The depth of the footprint correlates with the occupancy of the factor, allowing quantitative comparisons between conditions.

Applications of ATAC-seq

Mapping regulatory elements

The most direct application of ATAC-seq is the genome-wide identification of regulatory elements. In a typical human cell line, ATAC-seq identifies 50,000–150,000 peaks, of which roughly 20–30% are at promoters and the remainder at enhancers, insulators, and other distal regulatory elements. By comparing ATAC-seq profiles across cell types or developmental stages, researchers can identify cell-type-specific regulatory elements and the transcription factors that define cellular identity.

ATAC-seq has been particularly valuable for identifying enhancers, which are difficult to predict from sequence alone. The accessibility of an enhancer correlates with its activity, and ATAC-seq peaks at enhancers are often associated with the binding of lineage-defining transcription factors. For example, in hematopoietic differentiation, ATAC-seq has revealed that enhancer accessibility changes dynamically as cells commit to specific lineages, with master regulators such as GATA1 and PU.1 establishing new enhancer landscapes.

The method is also used to study the effects of genetic variants on chromatin accessibility. Quantitative trait locus (QTL) mapping using ATAC-seq data (caQTLs) can identify single nucleotide polymorphisms that alter chromatin state, providing a mechanistic link between non-coding genetic variation and gene expression differences.

Single-cell ATAC-seq

Single-cell ATAC-seq (scATAC-seq) extends the method to individual cells, enabling the characterization of chromatin accessibility heterogeneity within complex tissues. The first scATAC-seq protocols, such as the microfluidics-based Fluidigm C1 approach, were limited to hundreds of cells. The development of droplet-based methods (10x Genomics Chromium, Bio-Rad ddSEQ) and combinatorial indexing (sci-ATAC-seq) has increased throughput to tens of thousands of cells per experiment.

In scATAC-seq, individual nuclei are compartmentalized in droplets or wells, and each nucleus is barcoded with a unique identifier during the tagmentation or amplification step. The resulting data are extremely sparse—each cell yields only 1,000–10,000 unique fragments, representing a small fraction of the accessible genome. Analysis therefore requires specialized tools that aggregate information across cells, such as latent semantic indexing (LSI) for dimensionality reduction and clustering.

scATAC-seq has been used to map the regulatory landscape of the developing brain, identify rare cell types in tumors, and reconstruct differentiation trajectories. A key advantage over scRNA-seq is that scATAC-seq provides direct information about the regulatory state of a cell, which is often more stable and cell-type-specific than the transcriptome.

Integrating with other omics

ATAC-seq data are most powerful when integrated with other genomic measurements. The most common integration is with RNA-seq: genes with accessible promoters and enhancers tend to be expressed, and the combination of accessibility and expression data can identify regulatory elements that control cell-type-specific gene expression programs. Tools such as Cicero link distal enhancers to their target promoters by analyzing co-accessibility patterns across single cells.

Integration with __MASK_7__ data for histone modifications and transcription factors provides a more complete picture of the regulatory landscape. For example, ATAC-seq peaks that overlap H3K27ac (a mark of active enhancers) are more likely to be functional than peaks lacking this modification. Conversely, peaks that overlap H3K27me3 (a repressive mark) may represent poised regulatory elements.

ATAC-seq can also be combined with __MASK_8__ to relate chromatin accessibility to DNA methylation. Accessible regions are typically hypomethylated, and the joint analysis of these two data types can identify regions where methylation changes precede or follow chromatin remodeling.

Advantages and Limitations

Advantages over other methods

ATAC-seq offers several practical advantages that have made it the dominant method for chromatin accessibility profiling. The most significant is the low input requirement. The standard protocol uses 50,000 cells, but the method works with as few as 500 cells, and even single cells with appropriate scaling. This makes ATAC-seq applicable to rare cell populations, such as circulating tumor cells, early embryos, and flow-sorted subpopulations.

The speed of the protocol is another major advantage. The entire experiment, from cells to sequencing-ready libraries, takes less than 3 hours. This is in contrast to DNase-seq, which requires multiple days and multiple purification steps. The reduced handling also means fewer opportunities for sample loss and technical variation.

The signal-to-noise ratio of ATAC-seq is superior to DNase-seq. Because Tn5 integrates only into accessible DNA and the reaction is distributive (each transposase molecule catalyzes one insertion), the background from non-specific cleavage is minimal. The fragment size distribution provides an internal control for the quality of the experiment, and the periodic pattern of nucleosomal fragments can be used to assess the extent of over-digestion.

Common limitations and how to mitigate

The most significant limitation of ATAC-seq is mitochondrial contamination. Because mitochondria lack nucleosomes, their DNA is highly accessible to Tn5 and can constitute the majority of sequencing reads. This is particularly problematic for samples with low nuclear content, such as red blood cells or platelets. Mitigation strategies include using a lysis buffer that selectively disrupts the nuclear membrane while leaving mitochondria intact, and computationally filtering mitochondrial reads after alignment.

The requirement for fresh or cryopreserved cells is another limitation. Unlike ChIP-seq, which can be performed on cross-linked cells, ATAC-seq requires intact nuclei with native chromatin structure. Cells that have been fixed with formaldehyde or stored for extended periods produce poor results. For clinical samples, this means that ATAC-seq must be performed shortly after collection, which can be logistically challenging.

Tn5 insertion bias is a third limitation. Although the hyperactive Tn5 variant has reduced sequence specificity, it still exhibits a preference for certain DNA structural features. This bias can affect the quantitative comparison of accessibility across loci, particularly for regions with unusual DNA geometry. Computational correction methods, such as the insertion bias correction implemented in HINT-ATAC, can partially mitigate this issue.

Finally, ATAC-seq provides a snapshot of accessibility but does not directly measure transcription factor binding or nucleosome occupancy. The interpretation of ATAC-seq data relies on the assumption that accessible regions are regulatory elements, which is generally true but not always. Integration with other data types, such as ChIP-seq and RNA-seq, is often necessary to draw functional conclusions.

Common Pitfalls and Troubleshooting

Optimizing nuclei number

The ratio of Tn5 transposase to nuclei is the most critical parameter in the ATAC-seq protocol. Too much transposase relative to nuclei results in over-tagmentation, producing libraries with a high proportion of very short fragments and reduced complexity. Too little transposase produces libraries with large fragments and low signal.

The standard protocol uses 2.5 µL of Tn5 (100 nM) for 50,000 cells. When scaling down to 500 cells, the reaction volume should be reduced proportionally (to 5 µL total) rather than keeping the volume constant and reducing the cell number. This maintains the transposase-to-nuclei ratio. For single-cell ATAC-seq, the transposase concentration is typically increased 10-fold to ensure that each nucleus receives sufficient enzyme.

If the library shows an unusually high proportion of fragments >1 kb, this indicates under-tagmentation. The remedy is to increase the amount of Tn5 or extend the incubation time. Conversely, if the library shows a single sharp peak at ~50 bp with no nucleosomal periodicity, this indicates over-tagmentation, and the Tn5 amount or incubation time should be reduced.

Controlling tagmentation time

The tagmentation reaction is typically incubated at 37°C for 30 minutes. However, the optimal time depends on the cell type and the downstream application. Cells with highly condensed chromatin, such as sperm or terminally differentiated neurons, may require longer incubation (45–60 minutes) to achieve sufficient fragmentation. Cells with open chromatin, such as embryonic stem cells, may require shorter incubation (15–20 minutes) to avoid over-digestion.

The temperature of the reaction also affects the outcome. Incubation at 37°C is standard, but some protocols use 30°C to slow the reaction and increase reproducibility. The reaction can be stopped at any time by adding EDTA, so it is possible to perform a time-course to determine the optimal conditions for a new cell type.

Dealing with mitochondrial reads

Mitochondrial reads are the most common source of wasted sequencing in ATAC-seq. In some cell types, particularly those with high mitochondrial content such as hepatocytes or muscle cells, mitochondrial reads can exceed 80% of the total. This is not only wasteful but also reduces the effective coverage of nuclear chromatin.

Several strategies can reduce mitochondrial contamination. The lysis buffer can be modified to include digitonin, which selectively permeabilizes the plasma membrane while leaving the mitochondrial membrane intact. A concentration of 0.01% digitonin in the lysis buffer has been shown to reduce mitochondrial reads without affecting nuclear accessibility. Alternatively, the nuclei can be purified by density gradient centrifugation (e.g., using a 30% Percoll gradient) to remove mitochondria before tagmentation.

Computationally, mitochondrial reads can be filtered by aligning to a reference that includes the mitochondrial chromosome and excluding those alignments. This is straightforward but does not recover the lost sequencing depth. For samples with severe mitochondrial contamination, it may be necessary to increase the total sequencing depth to achieve sufficient nuclear coverage.

Summary and Best Practices

Experimental checklist

1. Use fresh or freshly thawed cells; avoid fixed cells.
2. Optimize the lysis buffer for the cell type; consider digitonin for mitochondrial-rich samples.
3. Titrate the Tn5-to-nuclei ratio for each new cell type.
4. Monitor the tagmentation reaction; adjust time based on the fragment size distribution.
5. Determine the optimal PCR cycle number by qPCR; stop during the exponential phase.
6. Perform size selection to enrich for fragments <400 bp.
7. Verify library quality by Bioanalyzer or TapeStation before sequencing.
8. Apply the Tn5 shift and filter mitochondrial reads during data analysis.
9. Use the TSS enrichment score and fragment size periodicity as QC metrics.
10. Integrate ATAC-seq data with RNA-seq or ChIP-seq for functional interpretation.

Frequently Asked Questions

What is ATAC-seq?

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) is a method for mapping genome-wide chromatin accessibility. It uses a hyperactive Tn5 transposase to fragment and tag accessible chromatin regions with sequencing adapters, allowing the identification of regulatory elements such as promoters, enhancers, and insulators.

How does ATAC-seq work?

ATAC-seq works by incubating isolated nuclei with a Tn5 transposase pre-loaded with sequencing adapters. The transposase integrates into nucleosome-free regions of chromatin, simultaneously fragmenting the DNA and adding adapters. The tagged fragments are then amplified by PCR and sequenced. The resulting reads identify regions of accessible chromatin.

What are the applications of ATAC-seq?

ATAC-seq is used to map regulatory elements, identify cell-type-specific chromatin states, study the effects of genetic variants on accessibility, and characterize the regulatory landscape of complex tissues. Single-cell ATAC-seq extends these applications to individual cells, enabling the identification of rare cell types and the reconstruction of differentiation trajectories.

What is the difference between ATAC-seq and DNase-seq?

ATAC-seq uses Tn5 transposase to fragment and tag accessible chromatin in a single step, requires only 500–5,000 cells, and can be completed in hours. DNase-seq uses DNase I to cleave accessible chromatin, requires millions of cells, and takes several days. ATAC-seq also has a higher signal-to-noise ratio and supports single-cell analysis.

How many cells are needed for ATAC-seq?

The standard ATAC-seq protocol uses 50,000 cells, but the method works with as few as 500 cells. For single-cell ATAC-seq, individual nuclei are barcoded and analyzed separately, allowing the profiling of thousands to tens of thousands of cells per experiment.

What is single-cell ATAC-seq?

Single-cell ATAC-seq (scATAC-seq) is a variant of ATAC-seq that profiles chromatin accessibility in individual cells. Each nucleus is compartmentalized and barcoded, allowing the accessibility landscape of each cell to be determined separately. This enables the characterization of cellular heterogeneity and the identification of rare cell types.

Why are mitochondrial reads a problem in ATAC-seq?

Mitochondrial DNA lacks nucleosomes and is therefore highly accessible to Tn5 transposase. As a result, mitochondrial reads can constitute 30–80% of all sequencing reads, wasting sequencing depth and reducing coverage of nuclear chromatin. Mitochondrial reads can be reduced by modifying the lysis buffer or filtered computationally after alignment.

What are common pitfalls in ATAC-seq?

Common pitfalls include over-tagmentation (producing libraries with low complexity), under-tagmentation (producing libraries with large fragments and low signal), excessive PCR amplification (producing duplicate reads), and mitochondrial contamination. These issues can be mitigated by optimizing the transposase-to-nuclei ratio, controlling the tagmentation time, determining the optimal PCR cycle number, and using appropriate lysis conditions.

Key Takeaways

• ATAC-seq maps chromatin accessibility genome-wide using Tn5 transposase, which simultaneously fragments and tags accessible DNA with sequencing adapters.
• The method requires only 500–5,000 cells for bulk analysis and can be completed in under 3 hours, making it faster and more sensitive than DNase-seq or MNase-seq.
• The tagmentation reaction is controlled by the transposase-to-nuclei ratio, incubation time, and temperature; these parameters must be optimized for each cell type.
• Data analysis requires the Tn5 shift (+4 bp on the positive strand, −5 bp on the negative strand), filtering of mitochondrial reads, and peak calling with MACS2.
• The fragment size distribution provides an internal quality control, with peaks at ~100 bp (nucleosome-free), ~200 bp (mononucleosome), and ~400 bp (dinucleosome).
• Single-cell ATAC-seq enables the characterization of chromatin accessibility heterogeneity and the identification of rare cell types in complex tissues.
• ATAC-seq data are most powerful when integrated with other omics data, such as RNA-seq for gene expression and ChIP-seq for histone modifications and transcription factor binding.

Further Reading

• Sun Y, Miao N, Sun T. Detect accessible chromatin using ATAC-sequencing, from principle to applications. Hereditas. 2019. PubMed 31427911
• Marshall AS, Jones NS. Discovering Cellular Mitochondrial Heteroplasmy Heterogeneity with Single Cell RNA and ATAC Sequencing. Biology. 2021. PubMed 34198745
• Baek S, Lee I. Single-cell ATAC sequencing analysis: From data preprocessing to hypothesis generation. Computational and structural biotechnology journal. 2020. PubMed 32637041
• Craig AJ et al. Genome-wide profiling of transcription factor activity in primary liver cancer using single-cell ATAC sequencing. Cell reports. 2023. PubMed 37980571
• Wang D et al. Integrated single-cell RNA and ATAC sequencing of B-cell lymphoma-3(Bcl3) and endothelin-2(Edn2) proteins as targets to prevent glaucoma progression. International journal of biological macromolecules. 2025. PubMed 40562139

Related Topics

• Sequencing Coverage
• Deep Sequencing
• Gene Sequencing

Related Clinical & Scientific Guides

• MAPK Pathway: Mechanism, Function, and Clinical Relevance
• Mammalian Cell Culture Bioreactors: A Practical Guide
• Nucleotide Formation: Biosynthesis and Assembly of DNA/RNA Building Blocks