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.

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

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

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

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

Mapping enhancer-gene regulatory interactions from single-cell data

Linking enhancers to genes

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

Enhlink infers distal and context-specific enhancer–promoter linkages

Cell-type-specific mapping of enhancers and target genes from single-cell multimodal data

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

Single-cell transcriptomics unveils gene regulatory network plasticity

Global prediction of chromatin accessibility using small-cell-number and single-cell RNA-seq

Profiling the genetic determinants of chromatin accessibility with scalable single-cell CRISPR screens

MOCHA: advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human disease cohorts

Linking Genotypes with Multiple Phenotypes in Single-Cell CRISPR Screens

High-throughput single-cell chromatin accessibility CRISPR screens enable unbiased identification of regulatory networks in cancer

scooby: Modeling multi-modal genomic profiles from DNA sequence at single-cell resolution

Sensitive dissection of a genomic regulatory landscape using bulk and targeted single-cell activation

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

MOCHA's advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human cohorts

Nested epistasis enhancer networks for robust genome regulation