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.

Functional Inference of Gene Regulation Using Single-Cell Multi-Omics

Single-cell Multi-Scale Footprinting Reveals the Modular Organization of DNA Regulatory Elements

Chromatin Potential Identified by Shared Single Cell Profiling of RNA and Chromatin

Single-cell gene regulation network inference by large-scale data integration

Gene regulatory network inference in the era of single-cell multi-omics

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

A single-cell multimodal view on gene regulatory network inference from transcriptomics and chromatin accessibility data

Epiregulon: Inference of single-cell transcription factor activity to predict drug response and drivers of cell state

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets

Functional Cis -Regulatory Genomics for Systems Biology

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

Dissecting and improving gene regulatory network inference using single-cell transcriptome data

Inferring gene regulation from stochastic transcriptional variation across single cells at steady state

Understanding Transcriptional Regulation Using De-novo Sequence Motif Discovery, Network Inference and Interactome Data

Single-cell and spatial multiomic inference of gene regulatory networks using SCRIPro

Using single cell atlas data to reconstruct regulatory networks

SimiC enables the inference of complex gene regulatory dynamics across cell phenotypes

Multiplex profiling of developmental cis- regulatory elements with quantitative single-cell expression reporters

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

Robustness and applicability of transcription factor and pathway analysis tools on single-cell RNA-seq data

Single-cell transcriptomics unveils gene regulatory network plasticity