Abstract:The increasing availability of large-scale single-cell datasets has enabled the detailed description of cell states across multiple biological conditions and perturbations. In parallel, recent advances in unsupervised machine learning, particularly in transfer learning, have enabled fast and scalable mapping of these new single-cell datasets onto reference atlases. The resulting large-scale machine learning models however often have millions of parameters, rendering interpretation of the newly mapped datasets challenging. Here, we propose expiMap, a deep learning model that enables interpretable reference mapping using biologically understandable entities, such as curated sets of genes and gene programs. The key concept is the substitution of the uninterpretable nodes in an autoencoder’s bottleneck by labeled nodes mapping to interpretable lists of genes, such as gene ontologies, biological pathways, or curated gene sets, for which activities are learned as constraints during reconstruction. This is enabled by the incorporation of predefined gene programs into the reference model, and at the same time allowing the model to learn de novo new programs and refine existing programs during reference mapping. We show that the model retains similar integration performance as existing methods while providing a biologically interpretable framework for understanding cellular behavior. We demonstrate the capabilities of expiMap by applying it to 15 datasets encompassing five different tissues and species. The interpretable nature of the mapping revealed unreported associations between interferon signaling via the RIG-I/MDA5 and GPCRs pathways, with differential behavior in CD8+ T cells and CD14+ monocytes in severe COVID-19, as well as the role of annexins in the cellular communications between lymphoid and myeloid compartments for explaining patient response to the applied drugs. Finally, expiMap enabled the direct comparison of a diverse set of pancreatic beta cells from multiple studies where we observed a strong, previously unreported correlation between the unfolded protein response and asparagine N-linked glycosylation. Altogether, expiMap enables the interpretable mapping of single cell transcriptome data sets across cohorts, disease states and other perturbations. ### Competing Interest Statement Fabian J. Theis consults for Immunai Inc., Singularity Bio B.V., CytoReason Ltd, and Omniscope Ltd, and has ownership interest in Dermagnostix GmbH and Cellarity.

A deep generative model integrating single-cell time-frequency characteristics transformed from electrophysiological data with transcriptomic features

Using deep learning to quantify neuronal activation from single-cell and spatial transcriptomic data

Generative modeling and latent space arithmetics predict single-cell perturbation response across cell types, studies and species

Biologically Informed Deep Learning to Infer Gene Program Activity in Single Cells

MEA-seqX: High-resolution Profiling of Large-scale Electrophysiological and Transcriptional Network Dynamics

Deep inverse modeling reveals dynamic-dependent invariances in neural circuit mechanisms

SynapsNet: Enhancing Neuronal Population Dynamics Modeling via Learning Functional Connectivity

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

Interpretable modeling of time-resolved single-cell gene-protein expression with CrossmodalNet

A deep generative model for gene expression profiles from single-cell RNA sequencing

Targeting EEG/LFP Synchrony with Neural Nets

Integrating physiological and transcriptomic analyses at the single-neuron level

scNODE : Generative Model for Temporal Single Cell Transcriptomic Data Prediction

Inferring predictive genetic models and regulatory elements by deep learning of cross-species single-cell gene expression landscapes

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

Reconstructing Circuit Connectivity from in vivo Spike Trains Using Deep Domain-Adaptive Matching

Towards a "universal translator" for neural dynamics at single-cell, single-spike resolution

Electrophysiology prediction of single neurons based on their morphology

Functional Connectivity: Continuous-Time Latent Factor Models for Neural Spike Trains

Modeling nascent transcription from chromatin landscape and structure

Multimodal Charting of Molecular and Functional Cell States Via in Situ Electro-Sequencing.