Meningitis caused by Streptococcus suis serotype 2 (SS2) in humans and pigs is an acute nervous disorder associated with serious sequelae. Bacterial meningitis is tightly associated with immune cell responses and the local immune microenvironment. However, the dynamic changes of the immune system during the disease progression in the brain remains unclear. Here, single-cell mass cytometry analyses are used to comprehensively profile the composition and phenotypes of female mouse brain immune cells at different stages of SS2 meningitis. Ten major immune cell lineages are identified among which T cells and dendritic cells significantly increased during meningitis, with B cells increasing in the late stage. Specifically, SS2⁺PD-L1⁺ neutrophils with strong phagocytosis, bactericidal and apoptotic effects accumulate in the acute phase of SS2 infection. Microglia sequentially display the features of homeostasis, proliferation, and activation (enhanced MHCII and TLR2 signals and TNF-α secretion) during the process of meningitis. Both border-associated and monocyte-derived macrophages contribute to the process of SS2-induced meningitis, exhibiting upregulation of CD38 and MHCII. Interestingly, CD11c⁺CD8⁺T cells are the main contributor of IFN-γ and specifically appeared during SS2 infection. In addition, the appearance of other lymphocytes such as CCR6^+/lo B cells, CX3CR1⁺ NK and MHCII⁺ ILC3 are related to the progression of meningitis. Moreover, correlation analysis between the composition of immune cell clusters and the SS2 infection process yield a dynamic immune landscape in which key immune clusters, including some previously unidentified, mark different stages of infection. Together, these data reveal the unique infection-stage immune microenvironment during the progression of meningitis caused by SS2 and provide resources for the analysis of immunological pathogenesis, potential diagnostic markers and therapeutic targets for bacterial meningitis. (Figure presented.)

IL-21/IL-21R signaling is crucial in various immune diseases and cellular development, however, its role in bacterial pneumonia remains unclear. Here, IL-21R knockout (IL-21R^−/−) mice were more susceptible to Actinobacillus pleuropneumoniae (APP) than wild-type (WT) mice. High-dimensional mass cytometry analysis revealed that IL-21R deficiency inhibited neutrophil activation, decreased the numbers of monocytes and proinflammatory macrophages, and augmented the defective CD3^low T cells in the lungs. Intracellular cytokine staining showed decreased IFN-γ/TNF-α/IL-6 production in IL-21R^−/− mice, particularly in CD8⁺ T cells. Furthermore, a previously unrecognized Ly6C⁺Ly6G⁺CD4⁺ T cell subset emerged only in the lungs of WT mice post-APP infection, which was in an activated status with stronger secretion capacities of IL-10, IL-21, granzyme B, and perforin by flow cytometry. These cells polarized macrophages into M2- or M1- phenotype without/with infection, respectively, and enhanced proliferation, phagocytosis, and macrophage extracellular traps/ROS-mediated bactericidal activity of macrophages against-APP, Klebsiella pneumoniae, or Escherichia coli infection. Thus, our study demonstrated that IL-21 drives the differentiation of neutrophils, monocytes, and macrophages into pro-inflammatory subsets. IL-21-induced Ly6C⁺Ly6G⁺CD4⁺ T cells cooperate with macrophages to enhance bacterial clearance, providing a promising target for preventing bacterial pneumonia. (Figure presented.)

The integration of spatial omics technologies can provide important insights into the biology of tissues. Here we combined mass spectrometry imaging-based metabolomics and imaging mass cytometry-based immunophenotyping on a single tissue section to reveal metabolic heterogeneity at single-cell resolution within tissues and its association with specific cell populations such as cancer cells or immune cells. This approach has the potential to greatly increase our understanding of tissue-level interplay between metabolic processes and their cellular components.

Rationale: Chronic inflammation plays an important role in alveolar tissue damage in emphysema, but the underlying immune alterations and cellular interactions are incompletely understood. Objectives: To explore disease-specific pulmonary immune cell alterations and cellular interactions in emphysema. Methods: We used single-cell mass cytometry (CyTOF) to compare the immune compartment in alveolar tissue from 15 patients with severe emphysema and 5 control subjects. Imaging mass cytometry (IMC) was applied to identify altered cell-cell interactions in alveolar tissue from patients with emphysema (n = 12) compared with control subjects (n = 8). Measurements and Main Results: We observed higher percentages of central memory CD4 T cells in combination with lower proportions of effector memory CD4 T cells in emphysema. In addition, proportions of cytotoxic central memory CD8 T cells and CD127+CD27+CD69- T cells were higher in emphysema, the latter potentially reflecting an influx of circulating lymphocytes into the lungs. Central memory CD8 T cells, isolated from alveolar tissue from patients with emphysema, exhibited an IFN-γ response upon anti-CD3 and anti-CD28 activation. Proportions of CD1c+ dendritic cells, expressing migratory and costimulatory markers, were higher in emphysema. Importantly, IMC enabled us to visualize increased spatial colocalization of CD1c+ dendritic cells and CD8 T cells in emphysema in situ. Conclusions: Using CyTOF, we characterized the alterations of the immune cell signature in alveolar tissue from patients with chronic obstructive pulmonary disease stage III or IV emphysema versus control lung tissue. These data contribute to a better understanding of the pathogenesis of emphysema and highlight the feasibility of interrogating the immune cell signature using CyTOF and IMC in human lung tissue. Clinical trial registered with www.clinicaltrials.gov (NCT04918706).

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.

Duchenne muscular dystrophy is caused by mutations in the DMD gene, leading to lack of dystrophin. Chronic muscle damage eventually leads to histological alterations in skeletal muscles. The identification of genes and cell types driving tissue remodeling is a key step to developing effective therapies. Here we use spatial transcriptomics in two Duchenne muscular dystrophy mouse models differing in disease severity to identify gene expression signatures underlying skeletal muscle pathology and to directly link gene expression to muscle histology. We perform deconvolution analysis to identify cell types contributing to histological alterations. We show increased expression of specific genes in areas of muscle regeneration (Myl4, Sparc, Hspg2), fibrosis (Vim, Fn1, Thbs4) and calcification (Bgn, Ctsk, Spp1). These findings are confirmed by smFISH. Finally, we use differentiation dynamic analysis in the D2-mdx muscle to identify muscle fibers in the present state that are predicted to become affected in the future state.

With the advent of multiplex fluorescence in situ hybridization (FISH) and in situ RNA sequencing technologies, spatial transcriptomics analysis is advancing rapidly, providing spatial location and gene expression information about cells in tissue sections at single cell resolution. Cell type classification of these spatially-resolved cells can be inferred by matching the spatial transcriptomics data to reference atlases derived from single cell RNA-sequencing (scRNA-seq) in which cell types are defined by differences in their gene expression profiles. However, robust cell type matching of the spatially-resolved cells to reference scRNA-seq atlases is challenging due to the intrinsic differences in resolution between the spatial and scRNA-seq data. In this study, we systematically evaluated six computational algorithms for cell type matching across four image-based spatial transcriptomics experimental protocols (MERFISH, smFISH, BaristaSeq, and ExSeq) conducted on the same mouse primary visual cortex (VISp) brain region. We find that many cells are assigned as the same type by multiple cell type matching algorithms and are present in spatial patterns previously reported from scRNA-seq studies in VISp. Furthermore, by combining the results of individual matching strategies into consensus cell type assignments, we see even greater alignment with biological expectations. We present two ensemble meta-analysis strategies used in this study and share the consensus cell type matching results in the Cytosplore Viewer (https://viewer.cytosplore.org) for interactive visualization and data exploration. The consensus matching can also guide spatial data analysis using SSAM, allowing segmentation-free cell type assignment.

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.

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.

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.

Immune checkpoint therapy (ICT) has the power to eradicate cancer, but the mechanisms that determine effective therapy-induced immune responses are not fully understood. Here, using high-dimensional single-cell profiling, we interrogate whether the landscape of T cell states in the peripheral blood predict responses to combinatorial targeting of the OX40 costimulatory and PD-1 inhibitory pathways. Single-cell RNA sequencing and mass cytometry expose systemic and dynamic activation states of therapy-responsive CD4⁺ and CD8⁺ T cells in tumor-bearing mice with expression of distinct natural killer (NK) cell receptors, granzymes, and chemokines/chemokine receptors. Moreover, similar NK cell receptor-expressing CD8⁺ T cells are also detected in the blood of immunotherapy-responsive cancer patients. Targeting the NK cell and chemokine receptors in tumor-bearing mice shows the functional importance of these receptors for therapy-induced anti-tumor immunity. These findings provide a better understanding of ICT and highlight the use and targeting of dynamic biomarkers on T cells to improve cancer immunotherapy.

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.

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.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy in need of effective (immuno)therapeutic treatment strategies. For the optimal application and development of cancer immunotherapies, a comprehensive understanding of local and systemic immune profiles in patients with PDAC is required. Here, our goal was to decipher the interplay between local and systemic immune profiles in treatment-naïve patients with PDAC. METHODS: The immune composition of PDAC, matched non-malignant pancreatic tissue, regional lymph nodes, spleen, portal vein blood, and peripheral blood samples (collected before and after surgery) from 11 patients with PDAC was assessed by measuring 41 immune cell markers by single-cell mass cytometry. Furthermore, the activation potential of tumor-infiltrating lymphocytes as determined by their ability to produce cytokines was investigated by flow cytometry. In addition, the spatial localization of tumor-infiltrating innate lymphocytes in the tumor microenvironment was confirmed by multispectral immunofluorescence. RESULTS: We found that CD103+CD8+ T cells with cytotoxic potential are infrequent in the PDAC immune microenvironment and lack the expression of activation markers and checkpoint blockade molecule programmed cell death protein-1 (PD-1). In contrast, PDAC tissues showed a remarkable increased relative frequency of B cells and regulatory T cells as compared with non-malignant pancreatic tissues. Besides, a previously unappreciated innate lymphocyte cell (ILC) population (CD127-CD103+CD39+CD45RO+ ILC1-like) was discovered in PDAC tissues. Strikingly, the increased relative frequency of B cells and regulatory T cells in pancreatic cancer samples was reflected in matched portal vein blood samples but not in peripheral blood, suggesting a regional enrichment of immune cells that infiltrate the PDAC microenvironment. After surgery, decreased frequencies of myeloid dendritic cells were found in peripheral blood. CONCLUSIONS: Our work demonstrates an immunosuppressive landscape in PDAC tissues, generally deprived of cytotoxic T cells and enriched in regulatory T cells and B cells. The antitumor potential of ILC1-like cells in PDAC may be exploited in a therapeutic setting. Importantly, immune profiles detected in blood isolated from the portal vein reflected the immune cell composition of the PDAC microenvironment, suggesting that this anatomical location could be a source of tumor-associated immune cell subsets.

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.

Understanding the mechanisms and impact of booster vaccinations are essential in the design and delivery of vaccination programs. Here we show that a three dose regimen of a synthetic peptide vaccine elicits an accruing CD8⁺ T cell response against one SARS-CoV-2 Spike epitope. We see protection against lethal SARS-CoV-2 infection in the K18-hACE2 transgenic mouse model in the absence of neutralizing antibodies, but two dose approaches are insufficient to confer protection. The third vaccine dose of the single T cell epitope peptide results in superior generation of effector-memory T cells and tissue-resident memory T cells, and these tertiary vaccine-specific CD8⁺ T cells are characterized by enhanced polyfunctional cytokine production. Moreover, fate mapping shows that a substantial fraction of the tertiary CD8⁺ effector-memory T cells develop from re-migrated tissue-resident memory T cells. Thus, repeated booster vaccinations quantitatively and qualitatively improve the CD8⁺ T cell response leading to protection against otherwise lethal SARS-CoV-2 infection.