JOINTLY: interpretable joint clustering of single-cell transcriptomes.

Single-cell and single-nucleus RNA-sequencing (sxRNA-seq) is increasingly being used to characterise the transcriptomic state of cell types at homeostasis, during development and in disease. However, this is a challenging task, as biological effects can be masked by technical variation. Here, we present JOINTLY, an algorithm enabling joint clustering of sxRNA-seq datasets across batches. JOINTLY performs on par or better than state-of-the-art batch integration methods in clustering tasks and outperforms other intrinsically interpretable methods. We demonstrate that JOINTLY is robust against over-correction while retaining subtle cell state differences between biological conditions and highlight how the interpretation of JOINTLY can be used to annotate cell types and identify active signalling programs across cell types and pseudo-time. Finally, we use JOINTLY to construct a reference atlas of white adipose tissue (WATLAS), an expandable and comprehensive community resource, in which we describe four adipocyte subpopulations and map compositional changes in obesity and between depots.

A single-cell, time-resolved profiling of Xenopus mucociliary epithelium reveals nonhierarchical model of development.

The specialized cell types of the mucociliary epithelium (MCE) lining the respiratory tract enable continuous airway clearing, with its defects leading to chronic respiratory diseases. The molecular mechanisms driving cell fate acquisition and temporal specialization during mucociliary epithelial development remain largely unknown. Here, we profile the developing Xenopus MCE from pluripotent to mature stages by single-cell transcriptomics, identifying multipotent early epithelial progenitors that execute multilineage cues before specializing into late-stage ionocytes and goblet and basal cells. Combining in silico lineage inference, in situ hybridization, and single-cell multiplexed RNA imaging, we capture the initial bifurcation into early epithelial and multiciliated progenitors and chart cell type emergence and fate progression into specialized cell types. Comparative analysis of nine airway atlases reveals an evolutionary conserved transcriptional module in ciliated cells, whereas secretory and basal types execute distinct function-specific programs across vertebrates. We uncover a continuous nonhierarchical model of MCE development alongside a data resource for understanding respiratory biology.

Predicting gene regulatory networks from cell atlases.

Recent single-cell RNA-sequencing atlases have surveyed and identified major cell types across different mouse tissues. Here, we computationally reconstruct gene regulatory networks from three major mouse cell atlases to capture functional regulators critical for cell identity, while accounting for a variety of technical differences, including sampled tissues, sequencing depth, and author assigned cell type labels. Extracting the regulatory crosstalk from mouse atlases, we identify and distinguish global regulons active in multiple cell types from specialised cell type-specific regulons. We demonstrate that regulon activities accurately distinguish individual cell types, despite differences between individual atlases. We generate an integrated network that further uncovers regulon modules with coordinated activities critical for cell types, and validate modules using available experimental data. Inferring regulatory networks during myeloid differentiation from wild-type and Irf8 KO cells, we uncover functional contribution of Irf8 regulon activity and composition towards monocyte lineage. Our analysis provides an avenue to further extract and integrate the regulatory crosstalk from single-cell expression data.