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.

Using high-throughput multi-omics data to investigate structural balance in elementary gene regulatory network motifs

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

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

Network Reconstruction for Trans Acting Genetic Loci Using Multi-Omics Data and Prior Information

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

Deciphering regulatory architectures from synthetic single-cell expression patterns

A Modular Framework for Gene Set Analysis Integrating Multilevel Omics Data

Topological benchmarking of algorithms to infer Gene Regulatory Networks from Single-Cell RNA-seq Data

Functionality Encoded In Topology? Discovering Macroscopic Regulatory Modules from Large-Scale Protein-DNA Interaction Networks

Understanding Distal Transcriptional Regulation from Sequence Motif, Network Inference and Interactome Perspectives

A Biophysical Model for ATAC-seq Data Analysis

Construction and Analysis of an Integrated Regulatory Network Derived from High-Throughput Sequencing Data

Quantifying the relationship between co-expression, co-regulation and gene function

A multifactorial evaluation framework for gene regulatory network reconstruction

Detecting phenotype-driven transitions in regulatory network structure

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

Identifying targets of multiple co-regulating transcription factors from expression time-series by Bayesian model comparison

A Bayesian data fusion based approach for learning genome-wide transcriptional regulatory networks

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

Dissecting endogeneous genetic circuits from first principles

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