Simultaneous measurements of transcriptomic and epigenomic profiles in the same individual cells provide an unprecedented opportunity to understand cell fates. However, effective approaches for the integrative analysis of such data are lacking. Here, we present a single-cell aggregation and integration (scAI) method to deconvolute cellular heterogeneity from parallel transcriptomic and epigenomic profiles. Through iterative learning, scAI aggregates sparse epigenomic signals in similar cells learned in an unsupervised manner, allowing coherent fusion with transcriptomic measurements. Simulation studies and applications to three real datasets demonstrate its capability of dissecting cellular heterogeneity within both transcriptomic and epigenomic layers and understanding transcriptional regulatory mechanisms.

Simultaneous measurements of transcriptomic and epigenomic profiles in the same individual cells provide an unprecedented opportunity to understand cell fates. However, effective approaches for the integrative analysis of such data are lacking. Here, we present a single-cell aggregation and integration (scAI) method to deconvolute cellular heterogeneity from parallel transcriptomic and epigenomic profiles. Through iterative learning, scAI aggregates sparse epigenomic signals in similar cells learned in an unsupervised manner, allowing coherent fusion with transcriptomic measurements. Simulation studies and applications to three real datasets demonstrate its capability of dissecting cellular heterogeneity within both transcriptomic and epigenomic layers and understanding transcriptional regulatory mechanisms.

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the "fate decisions" of transition cells within the transition. The present study provides new insights into transition events during differentiation and development.

Cell-cell communication is a fundamental process that shapes biological tissue. Historically, studies of cell-cell communication have been feasible for one or two cell types and a few genes. With the emergence of single-cell transcriptomics, we are now able to examine the genetic profiles of individual cells at unprecedented scale and depth. The availability of such data presents an exciting opportunity to construct a more comprehensive description of cell-cell communication. This review discusses the recent explosion of methods that have been developed to infer cell-cell communication from non-spatial and spatial single-cell transcriptomics, two promising technologies which have complementary strengths and limitations. We propose several avenues to propel this rapidly expanding field forward in meaningful ways.

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the "fate decisions" of transition cells within the transition. The present study provides new insights into transition events during differentiation and development.

In vivo calcium imaging enables simultaneous recording of large neuronal ensembles engaged in complex operations. Many experiments require monitoring and identification of cell populations across multiple sessions. Population cell tracking across multiple sessions is complicated by non-rigid transformations induced by cell movement and imaging field shifts. We introduce SCOUT (Single-Cell spatiOtemporal longitUdinal Tracking), a fast, robust cell-tracking method utilizing multiple cell-cell similarity metrics, probabilistic inference, and an adaptive clustering methodology, to perform cell identification across multiple sessions. By comparing SCOUT with earlier cell-tracking algorithms on simulated, 1-photon, and 2-photon recordings, we show that our approach significantly improves cell-tracking quality, particularly when recordings exhibit spatial footprint movement between sessions or sub-optimal neural extraction quality.

Motivation

Single-cell RNA-sequencing (scRNA-seq) offers unprecedented resolution for studying cellular decision-making processes. Robust inference of cell state transition paths and probabilities is an important yet challenging step in the analysis of these data.

Results

Here we present scEpath, an algorithm that calculates energy landscapes and probabilistic directed graphs in order to reconstruct developmental trajectories. We quantify the energy landscape using 'single-cell energy' and distance-based measures, and find that the combination of these enables robust inference of the transition probabilities and lineage relationships between cell states. We also identify marker genes and gene expression patterns associated with cell state transitions. Our approach produces pseudotemporal orderings that are-in combination-more robust and accurate than current methods, and offers higher resolution dynamics of the cell state transitions, leading to new insight into key transition events during differentiation and development. Moreover, scEpath is robust to variation in the size of the input gene set, and is broadly unsupervised, requiring few parameters to be set by the user. Applications of scEpath led to the identification of a cell-cell communication network implicated in early human embryo development, and novel transcription factors important for myoblast differentiation. scEpath allows us to identify common and specific temporal dynamics and transcriptional factor programs along branched lineages, as well as the transition probabilities that control cell fates.

Availability and implementation

A MATLAB package of scEpath is available at https://github.com/sqjin/scEpath.

Supplementary information

Supplementary data are available at Bioinformatics online.

The interfollicular epidermis (IFE) forms a water-tight barrier that is often disrupted in inflammatory skin diseases. During homeostasis, the IFE is replenished by stem cells in the basal layer that differentiate as they migrate toward the skin surface. Conventionally, IFE differentiation is thought to be stepwise as reflected in sharp boundaries between its basal, spinous, granular and cornified layers. The transcription factor GRHL3 regulates IFE differentiation by transcriptionally activating terminal differentiation genes. Here we use single cell RNA-seq to show that murine IFE differentiation is best described as a single step gradualistic process with a large number of transition cells between the basal and spinous layer. RNA-velocity analysis identifies a commitment point that separates the plastic basal and transition cell state from unidirectionally differentiating cells. We also show that in addition to promoting IFE terminal differentiation, GRHL3 is essential for suppressing epidermal stem cell expansion and the emergence of an abnormal stem cell state by suppressing Wnt signaling in stem cells.

Delayed and often impaired wound healing in the elderly presents major medical and socioeconomic challenges. A comprehensive understanding of the cellular/molecular changes that shape complex cell-cell communications in aged skin wounds is lacking. Here, we use single-cell RNA sequencing to define the epithelial, fibroblast, immune cell types, and encompassing heterogeneities in young and aged skin during homeostasis and identify major changes in cell compositions, kinetics, and molecular profiles during wound healing. Our comparative study uncovers a more pronounced inflammatory phenotype in aged skin wounds, featuring neutrophil persistence and higher abundance of an inflammatory/glycolytic Arg1^Hi macrophage subset that is more likely to signal to fibroblasts via interleukin (IL)-1 than in young counterparts. We predict systems-level differences in the number, strength, route, and signaling mediators of putative cell-cell communications in young and aged skin wounds. Our study exposes numerous cellular/molecular targets for functional interrogation and provides a hypothesis-generating resource for future wound healing studies.

Spinal motor neurons (MNs) integrate sensory stimuli and brain commands to generate movements. In vertebrates, the molecular identities of the cardinal MN types such as those innervating limb versus trunk muscles are well elucidated. Yet the identities of finer subtypes within these cell populations that innervate individual muscle groups remain enigmatic. Here we investigate heterogeneity in mouse MNs using single-cell transcriptomics. Among limb-innervating MNs, we reveal a diverse neuropeptide code for delineating putative motor pool identities. Additionally, we uncover that axial MNs are subdivided into three molecularly distinct subtypes, defined by mediolaterally-biased Satb2, Nr2f2 or Bcl11b expression patterns with different axon guidance signatures. These three subtypes are present in chicken and human embryos, suggesting a conserved axial MN expression pattern across higher vertebrates. Overall, our study provides a molecular resource of spinal MN types and paves the way towards deciphering how neuronal subtypes evolved to accommodate vertebrate motor behaviors.