Abstract:Ultrafast mapping of short reads to transcriptomic and metagenomic references via lightweight mapping techniques such as pseudoalignment has demonstrated success in substantially accelerating several types of analyses without much loss in accuracy compared to alignment-based approaches. The application of pseudoalignment to large reference sequences - like the genome - is, however, not trivial, due to the large size of the references or "targets" (i.e. chromosomes) and the presence of repetitive sequences within an individual reference sequence. This can lead to multiple matching locations for a k-mer within a single reference, which in turn can lead to false positive mappings and incorrect reference assignments for a read when the colors across the k-mer matches for a read are aggregated. Even when the read is determined to map to the appropriate reference, the increased occurrence of k-mer multi-matches within a reference can prevent the determination of the correct approximate position of the read, which is often critical in applications that map short reads to the genome. We propose a new and modified pseudoalignment scheme that partitions each reference into "virtual colors". These are essentially overlapping bins of fixed maximal extent on the reference sequences that are treated as distinct "colors" from the perspective of the pseudoalignment algorithm. A mapped k-mer is assigned a virtual color id that encodes the combination of the reference and within-reference bin in which the k-mer occurs. When the k-mers across a read are aggregated, the intersection is performed on virtual colors instead of the original colors (references), to determine the compatible set of targets (bins). The virtual colors can then be mapped back to the original references to provide the final mappings. The projection of the original reference sequences into virtual color space, and the corresponding modifications to the pseudoalignment procedure, can be applied dynamically at program invocation and without any modification of the underlying index itself. This makes the setting and modification of instance-appropriate parameters efficient and straightforward and the approach widely applicable. We apply this modified pseudoalignment procedure to process and map single-cell ATAC-seq data in our new tool alevin-fry-atac. We compare alevin-fry-atac to both Chromap and Cell Ranger ATAC. Alevin-fry-atac is highly scalable and, when using 32 threads, is approximately 1.78 times faster than Chromap (the second fastest approach) while using approximately 3 times less memory and mapping slightly more reads. The resulting peaks and clusters generated from alevin-fry-atac show high concordance with those obtained from both Chromap and the Cell Ranger ATAC pipeline, demonstrating that virtual color-enhanced pseudoalignment directly to the genome provides a fast, memory-frugal, and accurate alternative to existing approaches for single-cell ATAC-seq processing. The development of alevin-fry-atac brings single-cell ATAC-seq processing into a unified ecosystem with single-cell RNA-seq processing (via alevin-fry) to work toward providing a truly open alternative to many of the varied capabilities of CellRanger. Furthermore, our modified pseudoalignment approach should be easily applicable and extendable to other genome-centric mapping-based tasks and modalities such as standard DNA-seq, DNase-seq, Chip-seq and Hi-C.

ATAC-DEA: A Web-Based ATAC-Seq Data Differential Peak and Annotation Analysis Application

ATACgraph: Profiling Genome-Wide Chromatin Accessibility From ATAC-seq

PEPATAC: an optimized pipeline for ATAC-seq data analysis with serial alignments

Comparison of differential accessibility analysis strategies for ATAC-seq data

EAP: a versatile cloud-based platform for comprehensive and interactive analysis of large-scale ChIP/ATAC-seq data sets

ATAC-pipe: General Analysis of Genome-Wide Chromatin Accessibility

ATACdb: a Comprehensive Human Chromatin Accessibility Database.

W-ChIPeaks: a Comprehensive Web Application Tool for Processing ChIP-chip and ChIP-seq Data

CloudATAC: a cloud-based framework for ATAC-Seq data analysis

Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data

Cactus: a user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

An automated data analysis pipeline for GC-TOF-MS metabonomics studies.

DEApp: an interactive web interface for differential expression analysis of next generation sequence data

SnapATAC: A Comprehensive Analysis Package for Single Cell ATAC-seq

Incorporating network diffusion and peak location information for better single-cell ATAC-seq data analysis

DiffChIPL: a differential peak analysis method for high-throughput sequencing data with biological replicates based on limma

Fundamental and Practical Approaches for Single-Cell ATAC-seq Analysis

CoBRA: Containerized Bioinformatics Workflow for Reproducible ChIP/ATAC-seq Analysis

ATAC‐seq: A Method for Assaying Chromatin Accessibility Genome‐Wide

Alevin-fry-atac enables rapid and memory frugal mapping of single-cell ATAC-seq data using virtual colors for accurate genomic pseudoalignment

ShinyArchR.UiO: user-friendly,integrative and open-source tool for visualization of single-cell ATAC-seq data using ArchR