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.

Reverse engineering of regulatory networks in human B cells

Conditional Network Analysis Identifies Candidate Regulator Genes in Human B Cells

Genome-Wide Discovery of Modulators of Transcriptional Interactions in Human B Lymphocytes

Uncovering the regulatory landscape of early human B-cell lymphopoiesis and its implications in the pathogenesis of B-cell acute lymphoblastic leukemia

Stochastic modeling of a gene regulatory network driving B cell development in germinal centers

Identification of Gene Regulatory Networks in B-Cell Progenitor Differentiation and Leukemia

Characterization of Transcriptional Networks in Blood Stem and Progenitor Cells Using High-Throughput Single-Cell Gene Expression Analysis

Systematic discovery and perturbation of regulatory genes in human T cells reveals the architecture of immune networks

Identification of Hub Genes Associated with the Pathogenesis of Diffuse Large B-cell Lymphoma Subtype One Characterized by Host Response Via Integrated Bioinformatic Analyses.

Reverse Engineering Validation Using A Benchmark Synthetic Gene Circuit in Human Cells

Quantitative modelling of signaling in aggressive B cell lymphoma unveils conserved core network

Quantitative modeling of signaling in aggressive B cell lymphoma unveils conserved core network

Network-Constrained Eigen-Single-Cell Profile Estimation for Uncovering Crucial Immunogene Regulatory Systems in Human Bone Marrow

The Master Regulator Protein BAZ2B Can Reprogram Human Hematopoietic Lineage-Committed Progenitors into a Multipotent State

Temporal Genetic Program Following B-Cell Receptor Cross-Linking: Altered Balance Between Proliferation and Death in Healthy and Malignant B Cells

A Minimal Regulatory Network of Extrinsic and Intrinsic Factors Recovers Observed Patterns of CD4+ T Cell Differentiation and Plasticity

Quantitative modeling of the terminal differentiation of B cells and mechanisms of lymphomagenesis

Enforced MYC expression selectively redirects transcriptional programs during human plasma cell differentiation

A spatiotemporal map of co-receptor signaling networks underlying B cell activation

Leveraging chromatin accessibility for transcriptional regulatory network inference in T Helper 17 Cells

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