Abstract:Rapid advances in single-cell assays have outpaced methods for analysis of those data types. Different single-cell assays show extensive variation in sensitivity and signal to noise levels. In particular, scATAC-seq generates extremely sparse and noisy datasets. Existing methods developed to analyze this data require cells amenable to pseudo-time analysis or require datasets with drastically different cell-types. We describe a novel approach using self-organizing maps (SOM) to link scATAC-seq regions with scRNA-seq genes that overcomes these challenges and can generate draft regulatory networks. Our SOMatic package generates chromatin and gene expression SOMs separately and combines them using a linking function. We applied SOMatic on a mouse pre-B cell differentiation time-course using controlled Ikaros over-expression to recover gene ontology enrichments, identify motifs in genomic regions showing similar single-cell profiles, and generate a gene regulatory network that both recovers known interactions and predicts new Ikaros targets during the differentiation process. The ability of linked SOMs to detect emergent properties from multiple types of highly-dimensional genomic data with very different signal properties opens new avenues for integrative analysis of heterogeneous data.Gene expression is a tightly controlled process occurring in all cells during all stages of organismal life. How much and when genes are expressed is determined by gene regulatory networks (GRNs), which encode the biological programs that cells can perform. Each cell in an organism is constantly running through these networks to carry out its particular function. New techniques allow us to measure gene expression and chromatin accessibility using single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq). However, these techniques have relatively poor and different signal-to-noise ratios. In this work, we use a form of unsupervised learning called Self-Organizing Maps (SOMs) to analyze one step of B cell differentiation by linking separately trained scRNA-seq and scATAC-seq SOMs. We mine the linked SOMs to reconstruct the underlying GRN using a top-down approach. The resulting GRN not only recapitulated known regulatory linkages but also identified a large number of potential regulatory connections to the system. These methods should be generally applicable to linking heterogeneous high-throughput data with different signal-to-noise profiles.

GAGAM: a genomic annotation-based enrichment of scATAC-seq data for Gene Activity Matrix

Scart: Recognizing Cell Clusters and Constructing Trajectory from Single-Cell Epigenomic Data

Cross-Omic Transcription Factor Analysis: An Insight on Transcription Factor Accessibility and Expression Correlation

Cross-omic Transcription Factors meta-analysis: an insight on TFs accessibility and expression correlation

Deciphering cell types by integrating scATAC-seq data with genome sequences

Enhancement and Imputation of Peak Signal Enables Accurate Cell-Type Classification in scATAC-seq

Highly sensitive single-cell chromatin accessibility assay and transcriptome coassay with METATAC

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

GenoGAM: genome-wide generalized additive models for ChIP-Seq analysis

Benchmarking algorithms for gene set scoring of single-cell ATAC-seq data

Atlas-scale single-cell chromatin accessibility using nanowell-based combinatorial indexing

A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility

GFETM: Genome Foundation-based Embedded Topic Model for scATAC-seq Modeling

A Comprehensive Benchmarking Study on Computational Tools for Cross-omics Label Transfer from Single-cell RNA to ATAC Data

Scgat: A Cell-Type Annotation Framework for Single-Cell Transcriptomics Using Graph Attention Network and Meta Learning

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

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

Going beyond cell clustering and feature aggregation: Is there single cell level information in single-cell ATAC-seq data?

A Biophysical Model for ATAC-seq Data Analysis

TEA-seq: a trimodal assay for integrated single cell measurement of transcription, epitopes, and chromatin accessibility

Building gene regulatory networks from scATAC-seq and scRNA-seq using Linked Self Organizing Maps