RNA Velocity allows the inference of cellular differentiation trajectories from single-cell RNA sequencing (scRNA-seq) data. It would be highly interesting to study these differentiation dynamics in the spatial context of tissues. Estimating spatial RNA velocities is, however, limited by the inability to spatially capture spliced and unspliced mRNA molecules in high-resolution spatial transcriptomics. We present SIRV, a method to spatially infer RNA velocities at the single-cell resolution by enriching spatial transcriptomics data with the expression of spliced and unspliced mRNA from reference scRNA-seq data. We used SIRV to infer spatial differentiation trajectories in the developing mouse brain, including the differentiation of midbrain-hindbrain boundary cells and marking the forebrain origin of the cortical hem and diencephalon cells. Our results show that SIRV reveals spatial differentiation patterns not identifiable with scRNA-seq data alone. Additionally, we applied SIRV to mouse organogenesis data and obtained robust spatial differentiation trajectories. Finally, we verified the spatial RNA velocities obtained by SIRV using 10x Visium data of the developing chicken heart and MERFISH data from human osteosarcoma cells. Altogether, SIRV allows the inference of spatial RNA velocities at the single-cell resolution to facilitate studying tissue development.

Multi-omic analyses are necessary to understand the complex biological processes taking place at the tissue and cell level, but also to make reliable predictions about, for example, disease outcome. Several linear methods exist that create a joint embedding using paired information per sample, but recently there has been a rise in the popularity of neural architectures that embed paired -omics into the same non-linear manifold. This work describes a head-to-head comparison of linear and non-linear joint embedding methods using both bulk and single-cell multi-modal datasets. We found that non-linear methods have a clear advantage with respect to linear ones for missing modality imputation. Performance comparisons in the downstream tasks of survival analysis for bulk tumor data and cell type classification for single-cell data lead to the following insights: First, concatenating the principal components of each modality is a competitive baseline and hard to beat if all modalities are available at test time. However, if we only have one modality available at test time, training a predictive model on the joint space of that modality can lead to performance improvements with respect to just using the unimodal principal components. Second, -omic profiles imputed by neural joint embedding methods are realistic enough to be used by a classifier trained on real data with limited performance drops. Taken together, our comparisons give hints to which joint embedding to use for which downstream task. Overall, product-of-experts performed well in most tasks and was reasonably fast, while early integration (concatenation) of modalities did quite poorly.

Objective: Mass cytometry (MC) immunoprofiling allows high-parameter phenotyping of immune cells. We set to investigate the potential of MC immuno-monitoring of axial spondyloarthritis (axSpA) patients enrolled in the Tight Control SpondyloArthritis (TiCoSpA) trial. Methods: Fresh, longitudinal PBMCs samples (baseline, 24, and 48 weeks) from 9 early, untreated axSpA patients and 7 HLA-B27⁺ controls were analyzed using a 35-marker panel. Data were subjected to HSNE dimension reduction and Gaussian mean shift clustering (Cytosplore), followed by Cytofast analysis. Linear discriminant analyzer (LDA), based on initial HSNE clustering, was applied onto week 24 and 48 samples. Results: Unsupervised analysis yielded a clear separation of baseline patients and controls including a significant difference in 9 T cell, B cell, and monocyte clusters (cl), indicating disrupted immune homeostasis. Decrease in disease activity (ASDAS score; median 1.7, range 0.6–3.2) from baseline to week 48 matched significant changes over time in five clusters: cl10 CD4 T_nai cells median 4.7 to 0.02%, cl37 CD4 T_em cells median 0.13 to 8.28%, cl8 CD4 T_cm cells median 3.2 to 0.02%, cl39 B cells median 0.12 to 2.56%, and cl5 CD38⁺ B cells median 2.52 to 0.64% (all p<0.05). Conclusions: Our results showed that a decrease in disease activity in axSpA coincided with normalization of peripheral T- and B-cell frequency abnormalities. This proof of concept study shows the value of MC immuno-monitoring in clinical trials and longitudinal studies in axSpA. MC immunophenotyping on a larger, multi-center scale is likely to provide crucial new insights in the effect of anti-inflammatory treatment and thereby the pathogenesis of inflammatory rheumatic diseases. Key Points • Longitudinal immuno-monitoring of axSpA patients through mass cytometry indicates that normalization of immune cell compartments coincides with decrease in disease activity. • Our proof of concept study confirms the value of immune-monitoring utilizing mass cytometry.

The prognosis of high-grade serous ovarian carcinoma (HGSOC) is poor, and treatment selection is challenging. A heterogeneous tumor microenvironment (TME) characterizes HGSOC and influences tumor growth, progression, and therapy response. Better characterization with multidimensional approaches for simultaneous identification and categorization of the various cell populations is needed to map the TME complexity. While mass cytometry allows the simultaneous detection of around 40 proteins, the CyTOFmerge MATLAB algorithm integrates data sets and extends the phenotyping. This pilot study explored the potential of combining two datasets for improved TME phenotyping by profiling single-cell suspensions from ten chemo-naïve HGSOC tumors by mass cytometry. A 35-marker pan-tumor dataset and a 34-marker pan-immune dataset were analyzed separately and combined with the CyTOFmerge, merging 18 shared markers. While the merged analysis confirmed heterogeneity across patients, it also identified a main tumor cell subset, additionally to the nine identified by the pan-tumor panel. Furthermore, the expression of traditional immune cell markers on tumor and stromal cells was revealed, as were marker combinations that have rarely been examined on individual cells. This study demonstrates the potential of merging mass cytometry data to generate new hypotheses on tumor biology and predictive biomarker research in HGSOC that could improve treatment effectiveness.

Motivation: Single-cell technologies allow deep characterization of different molecular aspects of cells. Integrating these modalities provides a comprehensive view of cellular identity. Current integration methods rely on overlapping features or cells to link datasets measuring different modalities, limiting their application to experiments where different molecular layers are profiled in different subsets of cells. Results: We present scTopoGAN, a method for unsupervised manifold alignment of single-cell datasets with non-overlapping cells or features. We use topological autoencoders (topoAE) to obtain latent representations of each modality separately. A topology-guided Generative Adversarial Network then aligns these latent representations into a common space. We show that scTopoGAN outperforms state-of-the-art manifold alignment methods in complete unsupervised settings. Interestingly, the topoAE for individual modalities also showed better performance in preserving the original structure of the data in the low-dimensional representations when compared to other manifold projection methods. Taken together, we show that the concept of topology preservation might be a powerful tool to align multiple single modality datasets, unleashing the potential of multi-omic interpretations of cells.

In spatial transcriptomics (ST) data, biologically relevant features such as tissue compartments or cell-state transitions are reflected by gene expression gradients. Here, we present SpaceWalker, a visual analytics tool for exploring the local gradient structure of 2D and 3D ST data. The user can be guided by the local intrinsic dimensionality of the high-dimensional data to define seed locations, from which a flood-fill algorithm identifies transcriptomically similar cells on the fly, based on the high-dimensional data topology. In several use cases, we demonstrate that the spatial projection of these flooded cells highlights tissue architectural features and that interactive retrieval of gene expression gradients in the spatial and transcriptomic domains confirms known biology. We also show that SpaceWalker generalizes to several different ST protocols and scales well to large, multi-slice, 3D whole-brain ST data while maintaining real-time interaction performance.

Due to the increase in bacterial resistance, improving the anti-infectious immunity of the host is rapidly becoming a new strategy for the prevention and treatment of bacterial pneumonia. However, the specific lung immune responses and key immune cell subsets involved in bacterial infection are obscure. Actinobacillus pleuropneumoniae (APP) can cause porcine pleuropneumonia, a highly contagious respiratory disease that has caused severe economic losses in the swine industry. Here, using high-dimensional mass cytometry, the major immune cell repertoire in the lungs of mice with APP infection was profiled. Various phenotypically distinct neutrophil subsets and Ly-6C+ inflammatory monocytes/macrophages accumulated post-infection. Moreover, a linear differentiation trajectory from inactivated to activated to apoptotic neutrophils corresponded with the stages of uninfected, onset, and recovery of APP infection. CD14+ neutrophils, which mainly increased in number during the recovery stage of infection, were revealed to have a stronger ability to produce cytokines, especially IL-10 and IL-21, than their CD14- counterparts. Importantly, MHC-II+ neutrophils with antigen-presenting cell features were identified, and their numbers increased in the lung after APP infection. Similar results were further confirmed in the lungs of piglets infected with APP and Klebsiella pneumoniae infection by using a single-cell RNA-seq technique. Additionally, a correlation analysis between cluster composition and the infection process yielded a dynamic and temporally associated immune landscape where key immune clusters, including previously unrecognized ones, marked various stages of infection. Thus, these results reveal the characteristics of key neutrophil clusters and provide a detailed understanding of the immune response to bacterial pneumonia.

Motivation: Single-cell multi-omics assays simultaneously measure different molecular features from the same cell. A key question is how to benefit from the complementary data available and perform cross-modal clustering of cells. Results: We propose Single-Cell Multi-omics Clustering (scMoC), an approach to identify cell clusters from data with comeasurements of scRNA-seq and scATAC-seq from the same cell. We overcome the high sparsity of the scATAC-seq data by using an imputation strategy that exploits the less-sparse scRNA-seq data available from the same cell. Subsequently, scMoC identifies clusters of cells by merging clusterings derived from both data domains individually. We tested scMoC on datasets generated using different protocols with variable data sparsity levels. We show that scMoC (i) is able to generate informative scATAC-seq data due to its RNA-guided imputation strategy and (ii) results in integrated clusters based on both RNA and ATAC information that are biologically meaningful either from the RNA or from the ATAC perspective.

SummaryFactors that govern the complex formation of memory T cells are not completelyunderstood. A better understanding of thedevelopment of memory Tcell hetero-geneity is however required to enhance vaccination and immunotherapy ap-proaches. Here we examined the impact of pathogen- and tissue-specific cueson memory CD8+T cell heterogeneity using high-dimensional single-cell mass cy-tometry and a tailored bioinformatics pipeline. We identified distinct populationsof pathogen-specific CD8+T cells that uniquely connected to a specific pathogenor associated to multiple types of acute and persistent infections. In addition, thetissue environment shaped the memory CD8+T cell heterogeneity, albeit to alesser extent than infection. The programming of memory CD8+T cell differenti-ation during acute infection is eventually superseded by persistent infection.Thus, the plethora of distinct memory CD8+T cell subsets that arise upon infec-tion is dominantly sculpted by the pathogen-specific cues and further shaped by the tissue environment.

Single-cell technologies are emerging fast due to their ability to unravel the heterogeneity of biological systems. While scRNA-seq is a powerful tool that measures whole-transcriptome expression of single cells, it lacks their spatial localization. Novel spatial transcriptomics methods do retain cells spatial information but some methods can only measure tens to hundreds of transcripts. To resolve this discrepancy, we developed SpaGE, a method that integrates spatial and scRNA-seq datasets to predict whole-transcriptome expressions in their spatial configuration. Using five dataset-pairs, SpaGE outperformed previously published methods and showed scalability to large datasets. Moreover, SpaGE predicted new spatial gene patterns that are confirmed independently using in situ hybridization data from the Allen Mouse Brain Atlas.

Mass cytometry by time-of-flight (CyTOF) is a valuable technology for high-dimensional analysis at the single cell level. Identification of different cell populations is an important task during the data analysis. Many clustering tools can perform this task, which is essential to identify “new” cell populations in explorative experiments. However, relying on clustering is laborious since it often involves manual annotation, which significantly limits the reproducibility of identifying cell-populations across different samples. The latter is particularly important in studies comparing different conditions, for example in cohort studies. Learning cell populations from an annotated set of cells solves these problems. However, currently available methods for automatic cell population identification are either complex, dependent on prior biological knowledge about the populations during the learning process, or can only identify canonical cell populations. We propose to use a linear discriminant analysis (LDA) classifier to automatically identify cell populations in CyTOF data. LDA outperforms two state-of-the-art algorithms on four benchmark datasets. Compared to more complex classifiers, LDA has substantial advantages with respect to the interpretable performance, reproducibility, and scalability to larger datasets with deeper annotations. We apply LDA to a dataset of ~3.5 million cells representing 57 cell populations in the Human Mucosal Immune System. LDA has high performance on abundant cell populations as well as the majority of rare cell populations, and provides accurate estimates of cell population frequencies. Further incorporating a rejection option, based on the estimated posterior probabilities, allows LDA to identify previously unknown (new) cell populations that were not encountered during training. Altogether, reproducible prediction of cell population compositions using LDA opens up possibilities to analyze large cohort studies based on CyTOF data.