Single cell atac seq dimensionality reduction

single-cell-atac-seq-dimensionality-reduction

Summary

Reduce the dimensionality of scATAC-seq data using iterative latent semantic indexing (LSI) in ArchR to enable downstream clustering, visualization, and trajectory analysis. This skill is essential for handling the high-dimensional, sparse peak-by-cell matrices typical of single-cell chromatin accessibility data.

When to use

Apply this skill after loading and preprocessing raw scATAC-seq data into an ArchR project object when you need to compute low-dimensional embeddings for clustering, UMAP/tSNE visualization, or integrated multi-omic analysis. Use it as a prerequisite before trajectory analysis (Monocle3 or Slingshot), gene expression matrix integration, or any downstream analysis that requires dimensionality-reduced cell representations.

When NOT to use

• Input data is already a low-dimensional embedding or has been previously dimensionality-reduced by another method—apply this skill only on raw or minimally processed peak-by-cell matrices.
• scATAC-seq data has not been quality-filtered or peak-called; LSI performance depends on upstream preprocessing.
• Analysis goal is restricted to single-cell gene expression only without chromatin accessibility component; use gene expression-specific dimensionality reduction methods instead.

Inputs

• ArchR project object (with imported scATAC-seq peak matrix)
• Peak-by-cell count matrix (sparse, typically filtered for quality control)

Outputs

• ArchR project object with LSI embedding (addIterativeLSI result)
• Low-dimensional cell representation (LSI dimensions)
• Optional: combined LSI embedding (if using addCombinedDims for multiome data)

How to apply

Invoke the addIterativeLSI function on your ArchR project object to compute iterative latent semantic indexing across the peak-by-cell matrix. This function performs dimensionality reduction while accounting for sparsity inherent in scATAC-seq data. For paired scATAC-seq and scRNA-seq analysis, compute LSI embeddings for both modalities separately, then use addCombinedDims to integrate the reduced dimensions from both datasets into a unified embedding space suitable for joint clustering and downstream analysis.

Related tools

• ArchR (Primary package providing addIterativeLSI and addCombinedDims functions for scATAC-seq dimensionality reduction and multiome integration) — https://github.com/GreenleafLab/ArchR
• monocle3 (Downstream trajectory analysis tool that depends on dimensionality-reduced embeddings from LSI) — https://github.com/cole-trapnell-lab/monocle3
• Slingshot (Alternative downstream trajectory analysis tool that accepts LSI embeddings for pseudotime inference)
• R (Programming language and environment for executing ArchR functions)

Examples

library(ArchR); proj <- addIterativeLSI(ArchRProj = proj, useMatrix = 'TileMatrix', name = 'Iterative LSI')

Evaluation signals

• ArchR project object successfully contains LSI embedding(s) accessible via getReducedDims() or equivalent accessor.
• Dimensionality is reduced from tens of thousands of peaks to tens of LSI dimensions (typically 30–50 dimensions retained).
• Subsequent UMAP or tSNE visualization using the LSI embedding produces distinct cell clusters with biological meaning.
• If multiome: combined LSI embedding reflects both scATAC-seq and scRNA-seq signal without one modality dominating; validate via markers specific to each modality.
• Trajectory analysis (Monocle3 or Slingshot) downstream of LSI embedding converges and produces biologically coherent pseudotime ordering.

Limitations

• LSI assumes approximate low-rank structure in the peak-by-cell matrix; highly noisy or unfiltered datasets may produce poor embeddings.
• Iterative LSI is computationally intensive for >1 million cells, though ArchR is designed to handle this on standard laptops.
• Combined dimensionality reduction (addCombinedDims) may mask dataset-specific signals if the two modalities (scATAC and scRNA) have strong batch effects; batch correction before combination is recommended.
• ArchR is currently in beta with active development, so interface and behavior may change.

Evidence

• [intro] ArchR supports trajectory analysis using both monocle3 and Slingshot functions including addIterativeLSI as a prerequisite step: "ArchR now supports paired scATAC-seq and scRNA-seq Analysis! See updates with importFeatureMatrix, addGeneExpressionMatrix, addIterativeLSI, addCombinedDims"
• [readme] addIterativeLSI is a core dimensionality reduction workflow step in ArchR: "See updates with importFeatureMatrix, addGeneExpressionMatrix, addIterativeLSI, addCombinedDims"
• [readme] ArchR excels in speed and resource usage for large scATAC-seq datasets: "ArchR excels in both speed and resource usage, making it possible to analyze 1 million cells in 8 hours on a MacBook Pro laptop."
• [readme] ArchR is a comprehensive R package for scATAC-seq analysis: "ArchR is a full-featured R package for processing and analyzing single-cell ATAC-seq data."