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.

Analysis of multi-condition single-cell data with latent embedding multivariate regression

Disentangling cell type and state transcriptional programs

A machine learning one-class logistic regression model to predict stemness for single cell transcriptomics and spatial omics

A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

scEMB: Learning context representation of genes based on large-scale single-cell transcriptomics

scLEGA: an attention-based deep clustering method with a tendency for low expression of genes on single-cell RNA-seq data

Learning Consistent Subcellular Landmarks to Quantify Changes in Multiplexed Protein Maps.

LTMG: a novel statistical modeling of transcriptional expression states in single-cell RNA-Seq data

An interpretable single-cell RNA sequencing data clustering method based on latent Dirichlet allocation

Towards the Latent Transcriptome

LMD: Cluster-Independent Multiscale Marker Identification in Single-cell RNA-seq Data

Continuous lifelong learning for modeling of gene regulation from single cell multiome data by leveraging atlas-scale external data

Characterizing tissue composition through combined analysis of single-cell morphologies and transcriptional states

scELMo: Embeddings from Language Models are Good Learners for Single-cell Data Analysis

Joint Modeling of Cellular Heterogeneity and Condition Effects with scPCA in Single-Cell RNA-Seq

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

A flexible model for correlated count data, with application to multi-condition differential expression analyses of single-cell RNA sequencing data

Method of moments framework for differential expression analysis of single-cell RNA sequencing data

An expandable percutaneous catheter pump for left ventricular support: proof of concept.

A multi-view latent variable model reveals cellular heterogeneity in complex tissues for paired multimodal single-cell data

Regression-Based Analysis of Multimodal Single-Cell Data Integration Strategies