CyTOF workflow: differential discovery in high-throughput high-dimensional cytometry datasets

Plasmodium falciparum infection induces T cell tolerance that is associated with decreased disease severity upon re-infection

Immunity to severe malaria is acquired quickly, operates independently of pathogen load, and represents a highly effective form of disease tolerance. The mechanism that underpins tolerance remains unknown. We used a human rechallenge model of falciparum malaria in which healthy adult volunteers were infected three times over a 12 mo period to track the development of disease tolerance in real-time. We found that parasitemia triggered a hardwired innate immune response that led to systemic inflammation, pyrexia, and hallmark symptoms of clinical malaria across the first three infections of life. In contrast, a single infection was sufficient to reprogram T cell activation and reduce the number and diversity of effector cells upon rechallenge. Crucially, this did not silence stem-like memory cells but instead prevented the generation of cytotoxic effectors associated with autoinflammatory disease. Tolerized hosts were thus able to prevent collateral tissue damage in the absence of antiparasite immunity.

Comprehensive analysis of B cell repopulation in ocrelizumab-treated patients with multiple sclerosis by mass cytometry and proteomics

Ocrelizumab, an anti-CD20 antibody, depletes CD20+ B cells, which subsequently repopulate over months. Little is known about changes in other immune cell populations and molecular markers associated with B cell repopulation. Here, we performed a comprehensive characterization of immune cells from ocrelizumab-treated patients with multiple sclerosis (MS) using mass cytometry. About 50% of patients showed naive B cell repopulation after 6 months mainly with a transitional phenotype, whereas CD27+ memory B cells only rarely repopulated. This repopulation was associated with a reduction of memory T cells and activated myeloid cells, as well as reduced expression of activation/migration markers in both cell types. A plasma proteomics analysis identified proteins including TNFRSF13C, associated with B cell depletion and repopulation. Plasma levels of neurofilament light-chain protein declined after ocrelizumab treatment was not linked with B cell repopulation. These findings identify potential soluble markers for monitoring of ocrelizumab treatment in MS.

Understanding rheumatic disease through continuous cell state analysis

Autoimmune rheumatic diseases are a heterogeneous group of conditions, including rheumatoid arthritis (RA) and systemic lupus erythematosus. With the increasing availability of large single-cell datasets, novel disease-associated cell types continue to be identified and characterized at multiple omics layers, for example, 'T peripheral helper' (TPH) (CXCR5- PD-1hi) cells in RA and systemic lupus erythematosus and MerTK+ myeloid cells in RA. Despite efforts to define disease-relevant cell atlases, the very definition of a 'cell type' or 'lineage' has proven controversial as higher resolution assays emerge. This Review explores the cell types and states involved in disease pathogenesis, with a focus on the shifting perspectives on immune and stromal cell taxonomy. These understandings of cell identity are closely related to the computational methods adopted for analysis, with implications for the interpretation of single-cell data. Understanding the underlying cellular architecture of disease is also crucial for therapeutic research as ambiguity hinders translation to the clinical setting. We discuss the implications of different frameworks for cell identity for disease treatment and the discovery of predictive biomarkers for stratified medicine - an unmet clinical need for autoimmune rheumatic diseases.

CD147-high classical monocytes: a cellular biomarker for COVID-19 disease severity and treatment response

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can lead to severe coronavirus disease 2019 (COVID-19), which is characterized by cytokine storm and organ dysfunction. The spike S1 subunit induces inflammatory cytokine production, but the immune cell subsets that respond to S1 stimulation and contribute to disease severity remain unclear.

METHODS

We analyzed serum samples and peripheral blood mononuclear cells (PBMCs) from patients with COVID-19 (moderate: n = 7; severe: n = 25) and healthy controls (n = 38). Using mass cytometry (cytometry by time-of-flight; CyTOF), we analyzed immune cell responses to S1 subunit stimulation in PBMCs from healthy donors and patients with COVID-19. We examined correlations among identified cell populations, serum cytokine levels, and clinical parameters.

RESULTS

Serum S1 subunit levels correlated with disease severity and inflammatory cytokine concentrations. S1 subunit stimulation induced dose-dependent cytokine production from PBMCs, predominantly from myeloid cells. CyTOF analysis identified classical monocytes with high CD147 expression (CD147hi cMono) as the primary source of S1-induced cytokines. The proportion of CD147hi cMono increased significantly in severe COVID-19 and decreased with clinical improvement. The frequency of CD147hi cMono showed a stronger positive correlation with clinical severity markers in younger patients compared to older patients.

CONCLUSIONS

CD147hi cMono are the primary cellular source of S1-induced inflammatory cytokines and may serve as potential biomarkers for monitoring COVID-19 severity and treatment response.

Learning and teaching biological data science in the Bioconductor community

Modern biological research is increasingly data-intensive, leading to a growing demand for effective training in biological data science. In this article, we provide an overview of key resources and best practices available within the Bioconductor project-an open-source software community focused on omics data analysis. This guide serves as a valuable reference for both learners and educators in the field.

Immunodominant antiviral T cell responses outcompete immuno-subdominant antitumor responses to reduce the efficacy of oncolytic viroimmunotherapy

The paradigm in the field of oncolytic virotherapy proposes that tumor cell killing by an oncolytic virus (OV) culminates in the priming of antitumor CD8 T cells. However, this ignores the impact a highly immunodominant antiviral response against the OV has on the antitumor response which has been weakened by mechanisms of central tolerance. Here, we show that inflammatory Vesicular Stomatitis Virus (VSV) failed to prime an adoptively transferred, or pre-existing, population of tumor-reactive T cells. Combination with αPD1 immune checkpoint blockade therapy improved survival only when VSV expressed tumor associated antigens (TAA). These data show that, in this model, the highly inflammatory OV VSV alone actively outcompetes antitumor immunity. However, we also show that viral expression of a mutant near-self TAA can break central tolerance expanding heteroclitic self-reactive and near-self-reactive T cells, thus overcoming viral immunodominance by promoting tumor-specific T cell proliferation in parallel with expanding antiviral T cells.

Learning and teaching biological data science in the Bioconductor community

Modern biological research is increasingly data-intensive, leading to a growing demand for effective training in biological data science. In this article, we provide an overview of key resources and best practices available within the Bioconductor project - an open-source software community focused on omics data analysis. This guide serves as a valuable reference for both learners and educators in the field.

Tumor-specific CD8+ T cells from the bone marrow resist exhaustion and exhibit increased persistence in tumor-bearing hosts as compared with tumor-infiltrating lymphocytes

BACKGROUND

Immunotherapy is now an integral aspect of cancer therapy. Strategies employing adoptive cell therapy (ACT) have seen the establishment of chimeric antigen receptor (CAR)-T cells using peripheral blood lymphocytes as well as tumor-infiltrating lymphocytes (TILs) with significant clinical results. The bone marrow (BM) is an immunological niche housing T cells with specificity for previously encountered antigens, including tumor-associated antigens from certain solid cancers. This study sought to improve our understanding of tumor-specific BM T cells in the context of solid tumors by comparing them with TILs, and to assess whether there is a rationale for using the BM as a source of T cells for ACT against solid malignancies.

METHODS

We used the murine B16 melanoma model examining both the endogenous OVA-specific T cell response using an OVA-specific tetramer or examining the OVA-specific response with OVA-specific transgenic CD8+ (OT-1) T cells. Specifically, we compared baseline intrinsic properties of TILs or BM T cells from tumor-bearing mice and their changes following adoptive transfer in the tumor and bone marrow (as well as other compartments when indicated).

RESULTS

In tumor-bearing mice, endogenous tumor-specific T cells could be detected in the BM early in the course of tumor progression and possessed a more stem-cell-like and memory phenotype in an unsupervised cluster analysis compared with TILs which appeared more exhausted. The BM and tumor microenvironments significantly impact the fate of T cells. Naïve OT-1 transferred T cells acquired an exhausted phenotype in the tumor but maintained a more memory-like phenotype in the BM with tumor progression. Importantly, in a competitive transfer experiment, BM T cells infiltrated the tumor more efficiently than TILs, displayed a higher polyfunctionality with interleukin-2, interferon-γ, tumor necrosis factor-α production and showed greater persistence compared with TILs.

CONCLUSIONS

T cells from the BM appear superior to TILs as a source of cells for cellular therapy. They possess a memory-enriched phenotype and exhibit improved effector function, greater persistence within a tumor-bearing host, and the capacity for increased tumor infiltration. These data provide a foundation for further exploring the BM as a source of tumor-specific T cells for ACT in solid malignancies.

A new model measuring bacterial phagocytosis and phagolysosomal oxidation in humans using the intradermal injection of methylene blue-labeled Escherichia coli

Phagocytosis is an important leukocyte function; however, using existing models it cannot be measured in human tissues in vivo. To address this, we characterized a new phagocytosis model using intradermal methylene blue-labeled Escherichia coli injection (MBEC). Methylene blue (MB) is a licensed human medicine and bacterial stain potentially useful for labeling E. coli that is safe for human injection. Ex vivo coculture of leukocytes with MBEC caused MB to transfer into neutrophils and macrophages by phagocytosis. During this, a "red shift" in MB fluorescence was shown to be caused by phagolysosomal oxidation. Hence, MBEC coculture could be used to measure phagocytosis and phagolysosomal oxidation in humans, ex vivo. In healthy volunteers, inflammatory exudate sampling using suction blisters 2 to 24 h after intradermal MBEC injection showed that tissue-acquired neutrophils and monocytes contained more MB than their circulating counterparts, whereas blood and inflamed tissue T, B, and natural killer cells were MBlo. This was validated with spectral flow cytometry by visualizing the MB emission spectrum in tissue-acquired neutrophils. Neutrophil MB emission spectra demonstrated more red shift at 24 h compared with earlier time points, in keeping with progressive phagolysosomal MB oxidation in neutrophils over time in vivo. This new MBEC model can therefore measure bacterial phagocytosis and phagolysosomal oxidation in human skin, in vivo. This has a number of important research applications, e.g. in studying human phagocyte biology, testing novel antimicrobials, and understanding why certain groups such as males, the elderly or those with diabetes, recent surgery, or malnutrition are at increased risk of bacterial infection.

Integrative deep immune profiling of the elderly reveals systems-level signatures of aging, sex, smoking, and clinical traits

BACKGROUND

Aging increases disease susceptibility and reduces vaccine responsiveness, highlighting the need to better understand the aging immune system and its clinical associations. Studying the human immune system, however, remains challenging due to its complexity and significant inter-individual variability.

METHODS

We conducted an immune profiling study of 550 elderly participants (≥60 years) and 100 young controls (20-40 years) from the RESIST Senior Individuals (SI) cohort. Extensive demographic, clinical, and laboratory data were collected. Multi-color spectral flow cytometry and 48-plex plasma cytokine assays were used for deep immune phenotyping. Data were analyzed using unsupervised clustering and multi-dataset integration approaches.

FINDINGS

We studied 97 innate and adaptive immune cell populations, revealing intricate age- and sex-related changes in the elderly immune system. Our large sample size allowed detection of even subtle changes in cytokines and immune cell clusters. Integrative analysis combining clinical, laboratory, and immunological data revealed systems-level aging signatures, including shifts in specific immune cell subpopulations and cytokine concentrations (e.g., HGF and CCL27). Additionally, we identified unique immune signatures associated with smoking, obesity, and diseases such as osteoporosis, heart failure, and gout.

INTERPRETATION

This study provides one of the most comprehensive immune profiles of elderly individuals, uncovering high-resolution immune changes associated with aging. Our findings highlight clinically relevant immune signatures that enhance our understanding of aging-related diseases and could guide future research into new treatments, offering translational insights into human health and aging.

FUNDING

Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy-EXC 2155-project number 390874280.

Anti-Inflammatory Resveratrol Protects Mice From Early Mortality After Haematopoietic Stem Cell Transplantation

The occurrence of inflammation subsequent to haematopoietic stem cell transplantation is associated with an elevated risk of transplant-related mortality (TRM). However, the duration of inflammation and the potential efficacy of anti-inflammatory agents in reducing TRM remain uncertain. We performed a comprehensive investigation to examine the post-transplantation alterations of inflammatory mediators and to ascertain the correlation between inflammation level and TRM through the neutrophil-lymphocyte ratio, ELISAs and cytometric bead array. The findings revealed that the 30-day interval following transplantation is characterised by the most pronounced inflammatory response in both human and murine subjects, thereby elevating the risk of TRM. The inflammation is primarily caused by myeloid bias during haematopoietic reconstitution, which is a commonly overlooked aspect in clinical transplantation, additionally, a lesser extent of irradiation-induced injury. The administration of the anti-inflammatory agent resveratrol has the potential to reduce systemic inflammation and TRM by suppressing the NOD-like receptor signalling pathway and slowing down granulocyte implantation in HSCT mice. This approach did not impair the differentiation potential of haematopoietic stem cells. These findings demonstrate that the 30-day post-transplant period represents an opportunity to facilitate HSCT colonisation, mitigate transplant-related adverse effects, and potentially reap the benefits of anti-inflammatory treatments.

CXCL12 ameliorates neutrophilia and disease severity in SARS-CoV-2 infection

Neutrophils, particularly low-density neutrophils (LDNs), are believed to contribute to acute COVID-19 severity. Here, we showed that neutrophilia can be detected acutely and even months after SARS-CoV-2 infection in patients and mice, while neutrophil depletion reduced disease severity in mice. A key factor in neutrophilia and severe disease in infected mice was traced to the chemokine CXCL12 secreted by bone marrow cells and unexpectedly, endothelial cells. CXCL12 levels were negatively correlated with LDN numbers in longitudinal analyses of patient blood samples. CXCL12 blockade in SARS-CoV-2-infected mice increased blood/lung neutrophil numbers, thereby accelerating disease progression without changing lung virus titers. The exaggerated mortality caused by CXCL12 blockade could be reversed by neutrophil depletion. In addition, blocking interactions between SARS-CoV-2 and angiotensin-converting enzyme 2 (ACE2) reduced CXCL12 levels, suggesting a signal transduction from virus-mediated ACE2 ligation to increased CXCL12 secretion. Collectively, these results demonstrate a previously unappreciated role of CXCL12 in diminishing neutrophilia, including low-density neutrophilia, and its deleterious effects in SARS-CoV-2 infections. The results also support the involvement of SARS-CoV-2-endothelial cell interactions in viral pathogenesis.

Viral and immune dynamics of genital human papillomavirus infections in young women with high temporal resolution

Human papillomavirus (HPV) infections drive one in 20 new cancer cases, exerting a particularly high burden on women. Most anogenital HPV infections are cleared in less than two years, but the underlying mechanisms that favour persistence in around 10% of women remain largely unknown. Notwithstanding, it is precisely this information that is crucial for improving treatment, screening, and vaccination strategies. To understand viral and immune dynamics in non-persisting HPV infections, we set up an observational longitudinal cohort study with frequent on-site visits for biological sample collection. We enrolled 189 women aged from 18 to 25 and living in the area of Montpellier (France) between 2016 and 2020. We performed 974 on-site visits for a total of 1,619 months of follow-up. We collected data on virus load, local immune cell populations, local concentrations of cytokines, and circulating antibody titres. Using hierarchical Bayesian statistical modelling to simultaneously analyse the data from 164 HPV infections from 76 participants, we show that in two months after infection, HPV viral load in non-persisting infections reaches a plateau that lasts on average for 13 to 20 months (95% credibility interval) and is then followed by a rapid clearance phase. This first description of the dynamics of HPV infections comes with the identification of immune correlates associated with infection clearance, especially gamma-delta T cells and CXCL10 concentration. A limitation of this study on HPV kinetics is that many infection follow-ups are censored. Furthermore, some immune cell populations are difficult to label because cervical immunity is less well characterised than systemic immunity. These results open new perspectives for understanding the frontier between acute and chronic infections, and for controlling HPV-associated diseases, as well as for research on human cancers of infectious origin. Trial Registration: This trial was registered is registered at ClinicalTrials.gov under the ID NCT02946346. This study has been approved by the Comité de Protection des Personnes (CPP) Sud Méditerranée I (reference number 2016-A00712-49); by the Comité Consultatif sur le Traitement de l'Information en matière de Recherche dans le domaine de la Santé (reference number 16.504); by the Commission Nationale Informatique et Libertés (reference number MMS/ABD/ AR1612278, decision number DR-2016-488), by the Agence Nationale de Sécurité du Médicament et des Produits de Santé (reference 20160072000007).

Unveiling the power of high-dimensional cytometry data with cyCONDOR

High-dimensional cytometry (HDC) is a powerful technology for studying single-cell phenotypes in complex biological systems. Although technological developments and affordability have made HDC broadly available in recent years, technological advances were not coupled with an adequate development of analytical methods that can take full advantage of the complex data generated. While several analytical platforms and bioinformatics tools have become available for the analysis of HDC data, these are either web-hosted with limited scalability or designed for expert computational biologists, making their use unapproachable for wet lab scientists. Additionally, end-to-end HDC data analysis is further hampered due to missing unified analytical ecosystems, requiring researchers to navigate multiple platforms and software packages to complete the analysis. To bridge this data analysis gap in HDC we develop cyCONDOR, an easy-to-use computational framework covering not only all essential steps of cytometry data analysis but also including an array of downstream functions and tools to expand the biological interpretation of the data. The comprehensive suite of features of cyCONDOR, including guided pre-processing, clustering, dimensionality reduction, and machine learning algorithms, facilitates the seamless integration of cyCONDOR into clinically relevant settings, where scalability and disease classification are paramount for the widespread adoption of HDC in clinical practice. Additionally, the advanced analytical features of cyCONDOR, such as pseudotime analysis and batch integration, provide researchers with the tools to extract deeper insights from their data. We use cyCONDOR on a variety of data from different tissues and technologies demonstrating its versatility to assist the analysis of high-dimensional data from preprocessing to biological interpretation.

Multi-modal analysis reveals tumor and immune features distinguishing EBV-positive and EBV-negative post-transplant lymphoproliferative disorders

The oncogenic Epstein-Barr virus (EBV) can drive tumorigenesis with disrupted host immunity, causing malignancies including post-transplant lymphoproliferative disorders (PTLDs). PTLD can also arise in the absence of EBV, but the biological differences underlying EBV(+) and EBV(-) B cell PTLD and the associated host-EBV-tumor interactions remain poorly understood. Here, we reveal the core differences between EBV(+) and EBV(-) PTLD, characterized by increased expression of genes related to immune processes or DNA interactions, respectively, and the augmented ability of EBV(+) PTLD B cells to modulate the tumor microenvironment through elaboration of monocyte-attracting cytokines/chemokines. We create a reference resource of proteins distinguishing EBV(+) B lymphoma cells from EBV(-) B lymphoma including the immunomodulatory molecules CD300a and CD24, respectively. Moreover, we show that CD300a is essential for maximal survival of EBV(+) PTLD B lymphoma cells. Our comprehensive multi-modal analyses uncover the biological underpinnings of PTLD and offer opportunities for precision therapies.

Acute Myeloid Leukemia Skews Therapeutic WT1-specific CD8 TCR-T Cells Towards an NK-like Phenotype that Compromises Function and Persistence

Acute myeloid leukemia (AML) that is relapsed and/or refractory post-allogeneic hematopoietic cell transplantation (HCT) is usually fatal. In a prior study, we demonstrated that AML relapse in high-risk patients was prevented by post-HCT immunotherapy with Epstein-Barr virus (EBV)-specific donor CD8+ T cells engineered to express a high-affinity Wilms Tumor Antigen 1 (WT1)-specific T-cell receptor (TTCR-C4). However, in the present study, infusion of EBV- or Cytomegalovirus (CMV)-specific TTCR-C4 did not clearly improve outcomes in fifteen patients with active disease post-HCT. TCRC4-transduced EBV-specific T cells persisted longer post-transfer than CMV-specific T cells. Persisting TTCR-C4 skewed towards dysfunctional natural killer-like terminal differentiation, distinct from the dominant exhaustion programs reported for T-cell therapies targeting solid tumors. In one patient with active AML post-HCT, a sustained TTCR-C4 effector-memory profile correlated with long-term TTCR-C4 persistence and disease control. These findings reveal complex mechanisms underlying AML-induced T-cell dysfunction, informing future therapeutic strategies for addressing post-HCT relapse.

Single-Cell Analysis of Bone Marrow CD8+ T Cells in Myeloid Neoplasms Reveals Pathways Associated with Disease Progression and Response to Treatment with Azacitidine

Immunophenotypic analysis identified a BM CD57+CXCR3+ subset of CD8+ T cells associated with response to AZA in patients with MDS and AML. Single-cell RNA sequencing analysis revealed that IFN signaling is linked to the response to treatment, whereas TGF-β signaling is associated with treatment failure, providing insights into new therapeutic approaches.

Systemic longitudinal immune profiling identifies proliferating Treg cells as predictors of immunotherapy benefit: biomarker analysis from the phase 3 CONTINUUM and DIPPER trials

The identification of predictors for immunotherapy is often hampered by the absence of control groups in many studies, making it difficult to distinguish between prognostic and predictive biomarkers. This study presents biomarker analyses from the phase 3 CONTINUUM trial (NCT03700476), the first to show that adding anti-PD-1 (aPD1) to chemoradiotherapy (CRT) improves event-free survival (EFS) in patients with locoregionally advanced nasopharyngeal carcinoma. A dynamic single-cell atlas was profiled using mass cytometry on peripheral blood mononuclear cell samples from 12 pairs of matched relapsing and non-relapsing patients in the aPD1-CRT arm. Using a supervised representation learning algorithm, we identified a Ki67+ proliferating regulatory T cells (Tregs) population expressing high levels of activated and immunosuppressive molecules including FOXP3, CD38, HLA-DR, CD39, and PD-1, whose abundance correlated with treatment outcome. The frequency of these Ki67+ Tregs was significantly higher at baseline and increased during treatment in patients who relapsed compared to non-relapsers. Further validation through flow cytometry (n = 120) confirmed the predictive value of this Treg subset. Multiplex immunohistochemistry (n = 249) demonstrated that Ki67+ Tregs in tumors could predict immunotherapy benefit, with aPD1 improving EFS only in patients with low baseline levels of Ki67+ Tregs. These findings were further validated in the multicenter phase 3 DIPPER trial (n = 262, NCT03427827) and the phase 3 OAK trial of anti-PD-L1 immunotherapy in NSCLC, underscoring the predictive value of Ki67+ Treg frequency in identifying the beneficiaries of immunotherapy and potentially guiding personalized treatment strategies.

Parameter optimization for stable clustering using FlowSOM: a case study from CyTOF

High-dimensional cell phenotyping is a powerful tool to study molecular and cellular changes in health and diseases. CyTOF enables high-dimensional cell phenotyping using tens of surface and intra-cellular markers. To utilize the full potential of CyTOF, we need advanced clustering and machine learning methodologies to enable automated gating of the complex data. Here we show that critical modifications to a machine learning based FlowSOM package and precise parameter optimization can enable us to reliably analyze the complex CyTOF data. We show the impact of key parameters on clustering outcomes while addressing bugs within the publicly available package. We modified the FlowSOM pipeline to fix the bugs, enable scalability to handle large datasets and perform parameter optimization. We further validated this modified pipeline on a substantial external immunological dataset demonstrating the need of data-specific tailored parameter optimization to ensure reliable definition and interrogation of immune cell populations associated with immune disorders.

Mass cytometry identifies imbalance of multiple immune-cell subsets associated with biologics treatment in ankylosing spondylitis

OBJECTIVE

This study aims to comprehensively investigate immune-cell landscapes in ankylosing spondylitis (AS) patients and explore longitudinal immunophenotyping changes induced by biological agents.

METHODS

We employed mass cytometry with 35 cellular markers to analyze blood samples from 34 AS patients and 13 healthy controls (HC). Eleven AS patients were re-evaluated 1 month (4 patients) and 3 months (7 patients) after treatment with biological agents. Flow Self-Organizing Maps (FlowSOM) clustering was performed to identify specific cellular metaclusters. We compared cellular abundances across distinct subgroups and validated subset differences using gating strategies in flow cytometry scatter plots, visualized with FlowJo software. The proportions of differential subsets were then used for intercellular and clinical correlation analysis, as well as for constructing diagnostic models based on the random forest algorithm.

RESULTS

In AS patients, we identified and validated nine different immune-cell subsets compared to HC. Three subsets increased: helper T-cell 17 (Th17), mucosa-associated invariant T-cell (MAIT), and classical monocytes (CM). Six subsets decreased: effector memory T-cell (TEM), naïve B cells, transitional B cells, IL10+ memory B cells, non-classical monocytes (NCM), and neutrophils. Treatments with biological agents could rectify cellular abnormalities, particularly the imbalance of CM/NCM. Furthermore, these subsets may serve as biomarkers for assessing disease activity and constructing effective diagnostic models for AS.

CONCLUSION

These findings provide novel insights into the specific patterns of immune cell in AS, facilitating the further development of novel biomarkers and potential therapeutic targets for AS patients.

CyCadas: accelerating interactive annotation and analysis of clustered cytometry data

MOTIVATION

Single cell profiling by cytometry has emerged as a key technology in biology, immunology and clinical-translational medicine. The correct annotation, which refers to the identification of clusters as specific cell populations based on their marker expression, of clustered high-dimensional cytometry data, is a critical step of the analysis. Its accuracy determines the correct interpretation of the biological data. Despite the progress in various clustering algorithms, the annotation of clustered data still remains a manual, time consuming and error-prone task. We developed a user-friendly cluster annotation and differential abundance detection tool that can be applied on data generated with Self Organizing Map clustering algorithms, thus simplifying the annotation process of datasets that consist of hundreds or thousands of clusters.

RESULTS

We present Cytometry Cluster Annotation and Differential Abundance Suite (CyCadas), a semi-automated software tool that facilitates cluster annotation in cytometry data by offering both visual and computational guidance. CyCadas addresses the critical need for efficient and accurate annotation of high-resolution clustered cytometry data, significantly reducing the time needed to perform the analysis compared to both manual gating approaches and manual annotation of clustered data. The tool features a user-friendly interface, visual tools enabling data exploration and automated threshold estimation to separate negative and positive marker expression. It facilitates the definition and annotation of cell phenotypes among multiple clusters in a tree-based data structure. Finally, it calculates the abundance of various cell populations across the conditions with statistical interpretation. It is an ideal resource for researchers aiming to streamline their cytometry workflow.

AVAILABILITY AND IMPLEMENTATION

CyCadas is available as open source at: https://github.com/DII-LIH-Luxembourg/cycadas.

Low circulating tumor cell levels correlate with favorable outcomes and distinct biological features in multiple myeloma

There is growing interest in multiple myeloma (MM) circulating tumor cells (CTCs), but their rareness in peripheral blood (PB) and inconsistency in cutoffs question their clinical utility. Herein, we applied next-generation flow cytometry in 550 bone marrow (BM) and matched PB samples to define an optimal CTC cutoff for both transplant-eligible and transplant-ineligible newly diagnosed MM (NDMM) patients. Deep phenotyping was performed to investigate unique microenvironmental features associated with CTC dissemination. CTCs were detected in 90% of patients (median 0.01%; range: 0.0002%-12.6%) and increased levels associated with adverse features. Correlations were observed between high CTC percentages and a diffused MRI pattern, a distinct BM composition characterized by altered B-cell differentiation together with an expansion of effector cells and tumor-associated macrophages, as well as a greater phenotypic dissimilarity between BM and PB clonal cells. Progression-free survival (PFS) and overall survival (OS) gradually worsened with each logarithmic increment of CTCs. Conversely, NDMM patients without CTCs showed unprecedented outcomes, with 5-year PFS and OS rates of 83% and 97%, respectively. A cutoff of 0.02% CTCs was independent of the ISS, LDH, and cytogenetics in a multivariate analysis of risk factors for PFS. The 0.02% CTC cutoff synergized with the MGUS-like phenotype and the R-ISS for improving the risk stratification systems. MRD negativity was less frequent if CTCs were ≥0.02% at diagnosis, but whenever achieved, the poor prognosis of these patients was abrogated. This study shows the clinical utility of CTC assessment in MM and provides evidence toward a consensus cutoff for risk stratification.

ImmCellTyper facilitates systematic mass cytometry data analysis for deep immune profiling

Mass cytometry is a cutting-edge high-dimensional technology for profiling marker expression at the single-cell level, advancing clinical research in immune monitoring. Nevertheless, the vast data generated by cytometry by time-of-flight (CyTOF) poses a significant analytical challenge. To address this, we describe ImmCellTyper (https://github.com/JingAnyaSun/ImmCellTyper), a novel toolkit for CyTOF data analysis. This framework incorporates BinaryClust, an in-house developed semi-supervised clustering tool that automatically identifies main cell types. BinaryClust outperforms existing clustering tools in accuracy and speed, as shown in benchmarks with two datasets of approximately 4 million cells, matching the precision of manual gating by human experts. Furthermore, ImmCellTyper offers various visualisation and analytical tools, spanning from quality control to differential analysis, tailored to users' specific needs for a comprehensive CyTOF data analysis solution. The workflow includes five key steps: (1) batch effect evaluation and correction, (2) data quality control and pre-processing, (3) main cell lineage characterisation and quantification, (4) in-depth investigation of specific cell types; and (5) differential analysis of cell abundance and functional marker expression across study groups. Overall, ImmCellTyper combines expert biological knowledge in a semi-supervised approach to accurately deconvolute well-defined main cell lineages, while maintaining the potential of unsupervised methods to discover novel cell subsets, thus facilitating high-dimensional immune profiling.

Probing Dermal Immunity to Mycobacteria through a Controlled Human Infection Model

Cutaneous mycobacterial infections cause substantial morbidity and are challenging to diagnose and treat. An improved understanding of the dermal immune response to mycobacteria may inspire new therapeutic approaches. We conducted a controlled human infection study with 10 participants who received 2 × 106 CFUs of Mycobacterium bovis bacillus Calmette-Guérin (Tice strain) intradermally and were randomized to receive isoniazid or no treatment. Peripheral blood was collected at multiple time points for flow cytometry, bulk RNA sequencing (RNA-seq), and serum Ab assessments. Systemic immune responses were detected as early as 8 d postchallenge in this M. bovis bacillus Calmette-Guérin-naive population. Injection-site skin biopsies were performed at days 3 and 15 postchallenge and underwent immune profiling using mass cytometry and single-cell RNA-seq, as well as quantitative assessments of bacterial viability and burden. Molecular viability testing and standard culture results correlated well, although no differences were observed between treatment arms. Single-cell RNA-seq revealed various immune and nonimmune cell types in the skin, and communication between them was inferred by ligand-receptor gene expression. Day 3 communication was predominantly directed toward monocytes from keratinocyte, muscle, epithelial, and endothelial cells, largely via the migration inhibitory factor pathway and HLA-E-KLRK1 interaction. At day 15, communication was more balanced between cell types. These data reveal the potential role of nonimmune cells in the dermal immune response to mycobacteria and the utility of human challenge studies to augment our understanding of mycobacterial infections.

Myeloid-T cell interplay and cell state transitions associated with checkpoint inhibitor response in melanoma

BACKGROUND

The treatment of melanoma, the deadliest form of skin cancer, has greatly benefited from immunotherapy. However, many patients do not show a durable response, which is only partially explained by known resistance mechanisms.

METHODS

We performed single-cell RNA sequencing of tumor immune infiltrates and matched peripheral blood mononuclear cells of 22 checkpoint inhibitor (CPI)-naive stage III-IV metastatic melanoma patients. After sample collection, the same patients received CPI treatment, and their response was assessed.

FINDINGS

CPI responders showed high levels of classical monocytes in peripheral blood, which preferentially transitioned toward CXCL9-expressing macrophages in tumors. Trajectories of tumor-infiltrating CD8+ T cells diverged at the level of effector memory/stem-like T cells, with non-responder cells progressing into a state characterized by cellular stress and apoptosis-related gene expression. Consistently, predicted non-responder-enriched myeloid-T/natural killer cell interactions were primarily immunosuppressive, while responder-enriched interactions were supportive of T cell priming and effector function.

CONCLUSIONS

Our study illustrates that the tumor immune microenvironment prior to CPI treatment can be indicative of response. In perspective, modulating the myeloid and/or effector cell compartment by altering the described cell interactions and transitions could improve immunotherapy response.

FUNDING

This research was funded by Roche Pharma Research and Early Development.

PTP1B phosphatase dampens iPSC-derived neutrophil motility and antimicrobial function

Neutrophils are rapidly recruited to sites of infection and are critical for pathogen clearance. Therapeutic use of primary neutrophils has been limited, as they have a short lifespan and are not amenable to genetic manipulation. Human induced pluripotent stem cells (iPSCs) can provide a robust source of neutrophils for infusion and are genetically tractable. However, current work has indicated that dampened intracellular signaling limits iPSC-derived neutrophil (iNeutrophil) cellular activation and antimicrobial response. Here, we show that protein tyrosine phosphatase 1B (PTP1B) inhibits intracellular signaling and dampens iNeutrophil effector function. Deletion of the PTP1B phosphatase increased PI3K and ERK signaling and was associated with increased F-actin polymerization, cell migration, and phagocytosis. In contrast, other effector functions like NETosis and reactive oxygen species production were reduced. PTP1B-deficient neutrophils were more responsive to Aspergillus fumigatus and displayed rapid recruitment and control of hyphal growth. Accordingly, depletion of PTP1B increased production of inflammatory factors including the neutrophil chemokine interleukin-8. Taken together, these findings suggest that PTP1B limits iNeutrophil motility and antimicrobial function.

Peripheral immunophenotyping reveals lymphocyte stimulation in healthy women living with hereditary breast and ovarian cancer syndrome

Germline pathogenic variants in BRCA1 and BRCA2 (gpath(BRCA1/2)) represent genetic susceptibility for hereditary breast and ovarian cancer syndrome. Tumor-immune interactions are key contributors to breast cancer pathogenesis. Although earlier studies confirmed pro-tumorigenic immunological alterations in breast cancer patients, data are lacking in healthy carriers of gpath(BRCA1/2). Peripheral blood mononuclear cells of 66 women with or without germline predisposition or breast cancer were studied with a mass cytometry panel that identified 4 immune subpopulations of altered frequencies between healthy controls and healthy gpath(BRCA1) carriers, while no difference was observed in healthy gpath(BRCA2) carriers compared to controls. Moreover, 3 (one IgD-CD27+CD95+ B cell subpopulation and two CD45RA-CCR7+CD38+ CD4+ T cell subpopulations) out of these 4 subpopulations were also elevated in triple-negative breast cancer patients compared to controls. Our results reveal an activated peripheral immune phenotype in healthy carriers of gpath(BRCA1) that needs to be further elucidated to be leveraged in risk-reducing strategies.

TLR agonists polarize interferon responses in conjunction with dendritic cell vaccination in malignant glioma: a randomized phase II Trial

In this randomized phase II clinical trial, we evaluated the effectiveness of adding the TLR agonists, poly-ICLC or resiquimod, to autologous tumor lysate-pulsed dendritic cell (ATL-DC) vaccination in patients with newly-diagnosed or recurrent WHO Grade III-IV malignant gliomas. The primary endpoints were to assess the most effective combination of vaccine and adjuvant in order to enhance the immune potency, along with safety. The combination of ATL-DC vaccination and TLR agonist was safe and found to enhance systemic immune responses, as indicated by increased interferon gene expression and changes in immune cell activation. Specifically, PD-1 expression increases on CD4+ T-cells, while CD38 and CD39 expression are reduced on CD8+ T cells, alongside an increase in monocytes. Poly-ICLC treatment amplifies the induction of interferon-induced genes in monocytes and T lymphocytes. Patients that exhibit higher interferon response gene expression demonstrate prolonged survival and delayed disease progression. These findings suggest that combining ATL-DC with poly-ICLC can induce a polarized interferon response in circulating monocytes and CD8+ T cells, which may represent an important blood biomarker for immunotherapy in this patient population.Trial Registration: ClinicalTrials.gov Identifier: NCT01204684.

LectoScape: A Highly Multiplexed Imaging Platform for Glycome Analysis and Biomedical Diagnosis

Glycosylation, a fundamental biological process, involves the attachment of glycans to proteins, lipids, and RNA, and it plays a crucial role in various biological pathways. It is of great significance to obtain the precise spatial distribution of glycosylation modifications at the cellular and tissue levels. Here, we introduce LectoScape, an innovative method enabling detailed imaging of tissue glycomes with up to 1 μm resolution through image mass cytometry (IMC). This method utilizes 12 distinct, nonoverlapping lectins selected via microarray technology, enabling the multiplexed detection of a wide array of glycans. Furthermore, we developed an efficient labeling strategy for these lectins. Crucially, our approach facilitates the concurrent imaging of diverse glycan motifs, including N-glycan and O-glycan, surpassing the capabilities of existing technologies. Using LectoScape, we have successfully delineated unique glycan structures in various cell types, enhancing our understanding of the glycan distribution across human tissues. Our method has identified specific glycan markers, such as α2,3-sialylated Galβ1, 3GalNAc in O-glycan, and terminal GalNAc, as diagnostic indicators for cervical intraepithelial neoplasia. This highlights the potential of LectoScape in cancer diagnostics through the detection of abnormal glycosylation patterns.

Comparative single-cell multiplex immunophenotyping of therapy-naive patients with rheumatoid arthritis, systemic sclerosis, and systemic lupus erythematosus shed light on disease-specific composition of the peripheral immune system

Introduction

Systemic autoimmune diseases (SADs) are a significant burden on the healthcare system. Understanding the complexity of the peripheral immunophenotype in SADs may facilitate the differential diagnosis and identification of potential therapeutic targets.

Methods

Single-cell mass cytometric immunophenotyping was performed on peripheral blood mononuclear cells (PBMCs) from healthy controls (HCs) and therapy-naive patients with rheumatoid arthritis (RA), progressive systemic sclerosis (SSc), and systemic lupus erythematosus (SLE). Immunophenotyping was performed on 15,387,165 CD45+ live single cells from 52 participants (13 cases/group), using an antibody panel to detect 34 markers.

Results

Using the t-SNE (t-distributed stochastic neighbor embedding) algorithm, the following 17 main immune cell types were determined: CD4+/CD57- T cells, CD4+/CD57+ T cells, CD8+/CD161- T cells, CD8+/CD161+/CD28+ T cells, CD8dim T cells, CD3+/CD4-/CD8- T cells, TCRγ/δ T cells, CD4+ NKT cells, CD8+ NKT cells, classic NK cells, CD56dim/CD98dim cells, B cells, plasmablasts, monocytes, CD11cdim/CD172dim cells, myeloid dendritic cells (mDCs), and plasmacytoid dendritic cells (pDCs). Seven of the 17 main cell types exhibited statistically significant frequencies in the investigated groups. The expression levels of the 34 markers in the main populations were compared between HCs and SADs. In summary, 59 scatter plots showed significant differences in the expression intensities between at least two groups. Next, each immune cell population was divided into subpopulations (metaclusters) using the FlowSOM (self-organizing map) algorithm. Finally, 121 metaclusters (MCs) of the 10 main immune cell populations were found to have significant differences to classify diseases. The single-cell T-cell heterogeneity represented 64MCs based on the expression of 34 markers, and the frequency of 23 MCs differed significantly between at least twoconditions. The CD3- non-T-cell compartment contained 57 MCs with 17 MCs differentiating at least two investigated groups. In summary, we are the first to demonstrate the complexity of the immunophenotype of 34 markers over 15 million single cells in HCs vs. therapy-naive patients with RA, SSc, and SLE. Disease specific population frequencies or expression patterns of peripheral immune cells provide a single-cell data resource to the scientific community.

Clinical and cellular phenotypes resulting from a founder mutation in IL10RB

Inborn errors of immunity are a group of rare genetically determined diseases that impair immune system development or function. Many of these diseases include immune dysregulation, autoimmunity, or autoinflammation as prominent clinical features. In some children diagnosed with very early onset inflammatory bowel disease (VEOIBD), monogenic inborn errors of immune dysregulation underlie disease. We report a case of VEOIBD caused by a novel homozygous loss of function mutation in IL10RB. We use cytometry by time-of-flight with a broad panel of antibodies to interrogate the immunophenotype of this patient and detect reduced frequencies of CD4 and CD8 T cells with additional defects in some populations of T helper cells, innate-like T cells, and memory B cells. Finally, we identify the patient's mutation as a founder allele in an isolated indigenous population and estimate the age of this variant by studying the shared ancestral haplotype.

PPARγ attenuates cellular senescence of alveolar macrophages in asthma-COPD overlap

BACKGROUND

Asthma-chronic obstructive pulmonary disease (COPD) overlap (ACO) represents a complex condition characterized by shared clinical and pathophysiological features of asthma and COPD in older individuals. However, the pathophysiology of ACO remains unexplored. We aimed to identify the major inflammatory cells in ACO, examine senescence within these cells, and elucidate the genes responsible for regulating senescence.

METHODS

Bioinformatic analyses were performed to investigate major cell types and cellular senescence signatures in a public single-cell RNA sequencing (scRNA-Seq) dataset derived from the lung tissues of patients with ACO. Similar analyses were carried out in an independent cohort study Immune Mechanisms Severe Asthma (IMSA), which included bulk RNA-Seq and CyTOF data from bronchoalveolar lavage fluid (BALF) samples.

RESULTS

The analysis of the scRNA-Seq data revealed that monocytes/ macrophages were the predominant cell type in the lung tissues of ACO patients, constituting more than 50% of the cells analyzed. Lung monocytes/macrophages from patients with ACO exhibited a lower prevalence of senescence as defined by lower enrichment scores of SenMayo and expression levels of cellular senescence markers. Intriguingly, analysis of the IMSA dataset showed similar results in patients with severe asthma. They also exhibited a lower prevalence of senescence, particularly in airway CD206 + macrophages, along with increased cytokine expression (e.g., IL-4, IL-13, and IL-22). Further exploration identified alveolar macrophages as a major subtype of monocytes/macrophages driving cellular senescence in ACO. Differentially expressed genes related to oxidation-reduction, cytokines, and growth factors were implicated in regulating senescence in alveolar macrophages. PPARγ (Peroxisome Proliferator-Activated Receptor Gamma) emerged as one of the predominant regulators modulating the senescent signature of alveolar macrophages in ACO.

CONCLUSION

The findings suggest that senescence in macrophages, particularly alveolar macrophages, plays a crucial role in the pathophysiology of ACO. Furthermore, PPARγ may represent a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO.Key words ACO, Asthma, COPD, Macrophages, Senescence, PPARγ.

SUMOylation inhibitors activate anti-tumor immunity by reshaping the immune microenvironment in a preclinical model of hepatocellular carcinoma

PURPOSE

High levels of heterogeneity and immunosuppression characterize the HCC immune microenvironment (TME). Unfortunately, the majority of hepatocellular carcinoma (HCC) patients do not benefit from immune checkpoint inhibitors (ICIs) therapy. New small molecule therapies for the treatment of HCC are the goal of our research.

METHODS

SUMOylation inhibitors (TAK-981 and ML-792) were evaluated for the treatment of preclinical mouse HCC models (including subcutaneous and orthotopic HCC models). We profile immune cell subsets from tumor samples after SUMOylation inhibitors treatment using single-cell RNA sequencing (scRNA-seq), mass cytometry (CyTOF), flow cytometry, and multiple immunofluorescences (mIF).

RESULTS

We discover that SUMOylation is higher in HCC patient samples compared to normal liver tissue. TAK-981 and ML-792 decrease SUMOylation at nanomolar levels in HCC cells and also successfully reduced the tumor burden. Analysis combining scRNA-seq and CyTOF demonstrate that treatment with SUMOylation inhibitors reduces the exhausted CD8+T (Tex) cells while enhancing the cytotoxic NK cells, M1 macrophages and cytotoxic T lymphocytes (CTL) in preclinical mouse HCC model. Furthermore, SUMOylation inhibitors have the potential to activate innate immune signals from CD8+T, NK and macrophages while promoting TNFα and IL-17 secretion. Most notably, SUMOylation inhibitors can directly alter the TME by adjusting the abundance of intestinal microbiota, thereby restoring anti-tumor immunity in HCC models.

CONCLUSIONS

This preclinical study suggests that SUMO signaling inhibitors may be beneficial for the treatment of HCC.

Phase II Study of Eribulin plus Pembrolizumab in Metastatic Soft-tissue Sarcomas: Clinical Outcomes and Biological Correlates

PURPOSE

Eribulin modulates the tumor-immune microenvironment via cGAS-STING signaling in preclinical models. This non-randomized phase II trial evaluated the combination of eribulin and pembrolizumab in patients with soft-tissue sarcomas (STS).

PATIENTS AND METHODS

Patients enrolled in one of three cohorts: leiomyosarcoma (LMS), liposarcomas (LPS), or other STS that may benefit from PD-1 inhibitors, including undifferentiated pleomorphic sarcoma (UPS). Eribulin was administered at 1.4 mg/m2 i.v. (days 1 and 8) with fixed-dose pembrolizumab 200 mg i.v. (day 1) of each 21-day cycle, until progression, unacceptable toxicity, or completion of 2 years of treatment. The primary endpoint was the 12-week progression-free survival rate (PFS-12) in each cohort. Secondary endpoints included the objective response rate, median PFS, safety profile, and overall survival (OS). Pretreatment and on-treatment blood specimens were evaluated in patients who achieved durable disease control (DDC) or progression within 12 weeks [early progression (EP)]. Multiplexed immunofluorescence was performed on archival LPS samples from patients with DDC or EP.

RESULTS

Fifty-seven patients enrolled (LMS, n = 19; LPS, n = 20; UPS/Other, n = 18). The PFS-12 was 36.8% (90% confidence interval: 22.5-60.4) for LMS, 69.6% (54.5-89.0) for LPS, and 52.6% (36.8-75.3) for UPS/Other cohorts. All 3 patients in the UPS/Other cohort with angiosarcoma achieved RECIST responses. Toxicity was manageable. Higher IFNα and IL4 serum levels were associated with clinical benefit. Immune aggregates expressing PD-1 and PD-L1 were observed in a patient that completed 2 years of treatment.

CONCLUSIONS

The combination of eribulin and pembrolizumab demonstrated promising activity in LPS and angiosarcoma.

Fat-to-blood recirculation of partially dysfunctional PD-1+CD4 Tconv cells is associated with dysglycemia in human obesity

Obesity is characterized by the accumulation of T cells in insulin-sensitive tissues, including the visceral adipose tissue (VAT), that can interfere with the insulin signaling pathway eventually leading to insulin resistance (IR) and type 2 diabetes. Here, we found that PD-1+CD4 conventional T (Tconv) cells, endowed with a transcriptomic and functional profile of partially dysfunctional cells, are diminished in VAT of obese patients with dysglycemia (OB-Dys), without a concomitant increase in apoptosis. These cells showed enhanced capacity to recirculate into the bloodstream and had a non-restricted TCRβ repertoire divergent from that of normoglycemic obese and lean individuals. PD-1+CD4 Tconv were reduced in the circulation of OB-Dys, exhibited an altered migration potential, and were detected in the liver of patients with non-alcoholic steatohepatitis. The findings suggest a potential role for partially dysfunctional PD-1+CD4 Tconv cells as inter-organ mediators of IR in obese patients with dysglycemic.

PPARγ Attenuates Cellular Senescence of Alveolar Macrophages in Asthma- COPD Overlap

Asthma-chronic obstructive pulmonary disease (COPD) overlap (ACO) represents a complex condition characterized by shared clinical and pathophysiological features of asthma and COPD in older individuals. However, the pathophysiology of ACO remains unexplored. We aimed to identify the major inflammatory cells in ACO, examine senescence within these cells, and elucidate the genes responsible for regulating senescence. Bioinformatic analyses were performed to investigate major cell types and cellular senescence signatures in a public single-cell RNA sequencing (scRNA-Seq) dataset derived from the lung tissues of patients with ACO. Similar analyses were carried out in an independent cohort study Immune Mechanisms Severe Asthma (IMSA), which included bulk RNA-Seq and CyTOF data from bronchoalveolar lavage fluid (BALF) samples. The analysis of the scRNA-Seq data revealed that monocytes/ macrophages were the predominant cell type in the lung tissues of ACO patients, constituting more than 50% of the cells analyzed. Lung monocytes/macrophages from patients with ACO exhibited a lower prevalence of senescence as defined by lower enrichment scores of SenMayo and expression levels of cellular senescence markers. Intriguingly, analysis of the IMSA dataset showed similar results in patients with severe asthma. They also exhibited a lower prevalence of senescence, particularly in airway CD206 + macrophages, along with increased cytokine expression (e.g., IL-4, IL-13, and IL-22). Further exploration identified alveolar macrophages as a major subtype of monocytes/macrophages driving cellular senescence in ACO. Differentially expressed genes related to oxidation-reduction, cytokines, and growth factors were implicated in regulating senescence in alveolar macrophages. PPARγ (Peroxisome Proliferator-Activated Receptor Gamma) emerged as one of the predominant regulators modulating the senescent signature of alveolar macrophages in ACO. Collectively, the findings suggest that senescence in macrophages, particularly alveolar macrophages, plays a crucial role in the pathophysiology of ACO. Furthermore, PPARγ may represent a potential therapeutic target for interventions aimed at modulating senescence-associated processes in ACO.

Altered immune cell phenotypes within chronically ischemic human kidneys distal to occlusive renal artery disease

Renal artery stenosis (RAS) is a major cause of ischemic kidney disease, which is largely mediated by inflammation. Mapping the immune cell composition in ischemic kidneys might provide useful insight into the disease pathogenesis and uncover therapeutic targets. We used mass cytometry (CyTOF) to explore the single-cell composition in a unique data set of human kidneys nephrectomized due to chronic occlusive vascular disease (RAS, n = 3), relatively healthy donor kidneys (n = 6), and unaffected sections of kidneys with renal cell carcinoma (RCC, n = 3). Renal fibrosis and certain macrophage populations were also evaluated in renal sections. Cytobank analysis showed in RAS kidneys decreased cell populations expressing epithelial markers (CD45-/CD13+) and increased CD45+ inflammatory cells, whereas scattered tubular-progenitor-like cells (CD45-/CD133+/CD24+) increased compared with kidney donors. Macrophages switched to proinflammatory phenotypes in RAS, and the numbers of IL-10-producing dendritic cells (DC) were also lower. Compared with kidney donors, RAS kidneys had decreased overall DC populations but increased plasmacytoid DC. Furthermore, senescent active T cells (CD45+/CD28+/CD57+), aged neutrophils (CD45+/CD15+/CD24+/CD11c+), and regulatory B cells (CD45+/CD14-/CD24+/CD44+) were increased in RAS. RCC kidneys showed a distribution of cell phenotypes comparable with RAS but less pronounced, accompanied by an increase in CD34+, CD370+, CD103+, and CD11c+/CD103+ cells. Histologically, RAS kidneys showed significantly increased fibrosis and decreased CD163+/CD141+ cells. The single-cell platform CyTOF enables the detection of significant changes in renal cells, especially in subsets of immune cells in ischemic human kidneys. Endogenous pro-repair cell types in RAS warrant future study for potential immune therapy.NEW & NOTEWORTHY The single-cell platform mass cytometry (CyTOF) enables detection of significant changes in one million of renal cells, especially in subsets of immune cells in ischemic human kidneys distal to renal artery stenosis (RAS). We found that pro-repair cell types such as scattered tubular-progenitor-like cells, aged neutrophils, and regulatory B cells show a compensatory increase in RAS. Immune cell phenotype changes may reflect ongoing inflammation and impaired immune defense capability in the kidneys.

Ursa: A Comprehensive Multiomics Toolbox for High-Throughput Single-Cell Analysis

The burgeoning amount of single-cell data has been accompanied by revolutionary changes to computational methods to map, quantify, and analyze the outputs of these cutting-edge technologies. Many are still unable to reap the benefits of these advancements due to the lack of bioinformatics expertise. To address this issue, we present Ursa, an automated single-cell multiomics R package containing 6 automated single-cell omics and spatial transcriptomics workflows. Ursa allows scientists to carry out post-quantification single or multiomics analyses in genomics, transcriptomics, epigenetics, proteomics, and immunomics at the single-cell level. It serves as a 1-stop analytic solution by providing users with outcomes to quality control assessments, multidimensional analyses such as dimension reduction and clustering, and extended analyses such as pseudotime trajectory and gene-set enrichment analyses. Ursa aims bridge the gap between those with bioinformatics expertise and those without by providing an easy-to-use bioinformatics package for scientists in hoping to accelerate their research potential. Ursa is freely available at https://github.com/singlecellomics/ursa.

Associations of myeloid cells with cellular and humoral responses following vaccinations in patients with neuroimmunological diseases

Disease-modifying therapies (DMTs) are widely used in neuroimmunological diseases such as multiple sclerosis (MS) and neuromyelitis optica spectrum disorder (NMOSD). Although these treatments are known to predispose patients to infections and affect their responses to vaccination, little is known about the impact of DMTs on the myeloid cell compartment. In this study, we use mass cytometry to examine DMT-associated changes in the innate immune system in untreated and treated patients with MS (n = 39) or NMOSD (n = 23). We also investigated the association between changes in myeloid cell phenotypes and longitudinal responsiveness to homologous primary, secondary, and tertiary SARS-CoV-2 mRNA vaccinations. Multiple DMT-associated myeloid cell clusters, in particular CD64+HLADRlow granulocytes, showed significant correlations with B and T cell responses induced by vaccination. Our findings suggest the potential role of myeloid cells in cellular and humoral responses following vaccination in DMT-treated patients with neuroimmunological diseases.

High-dimensional mass cytometry reveals systemic and local immune signatures in necrotizing enterocolitis

Objective

Patients with necrotizing enterocolitis display severe gastrointestinal complications of prematurity, but the mechanism driving this clinical profile remains unknown. We used mass cytometry time-of-flight to characterize and compare immune cell populations in the blood and intestine tissue from patients with and without (controls) necrotizing enterocolitis at single-cell resolution.

Methods

We completed a deep mapping of the immune system of the peripheral blood mononuclear cells and intestinal mucosa tissue using mass cytometry to evaluate immune cell types, which revealed global immune dysregulation characteristics underlying necrotizing enterocolitis.

Results

Compared with controls, natural killer cells display signs of heightened activation and increased cytotoxic potential in the peripheral blood and mucosa of patients with necrotizing enterocolitis. Furthermore, CD4+ T effector memory cells, non-classical monocytes, active dendritic cells, and neutrophils were specifically enriched in the mucosa, suggesting trafficking from the periphery to areas of inflammation. Moreover, we mapped the systemic and local distinct immune signatures suggesting patterns of cell localization in necrotizing enterocolitis.

Conclusion

We used mass cytometry time-of-flight technology to identify immune cell populations specific to the peripheral blood and intestinal mucosa tissue from patients with necrotizing enterocolitis and controls. This information might be used to develop precise diagnosis and therapies that target specific cell populations in patients with necrotizing enterocolitis.

Interleukin-10-producing monocytes contribute to sex differences in pain resolution in mice and humans

Pain is closely associated with the immune system, which exhibits sexual dimorphism. For these reasons, neuro-immune interactions are suggested to drive sex differences in pain pathophysiology. However, our understanding of peripheral neuro-immune interactions on sex differences in pain resolution remains limited. Here, we have shown, in both a mouse model of inflammatory pain and in humans following traumatic pain, that males had higher levels of interleukin (IL)-10 than females, which were correlated with faster pain resolution. Following injury, we identified monocytes (CD11b+ Ly6C+ Ly6G-F4/80 mid ) as the primary source of IL-10, with IL-10-producing monocytes being more abundant in males than females. In a mouse model, neutralizing IL-10 signaling through antibodies, genetically ablating IL-10R1 in sensory neurons, or depleting monocytes with clodronate all impaired the resolution of pain hypersensitivity in both sexes. Furthermore, manipulating androgen levels in mice reversed the sexual dimorphism of pain resolution and the levels of IL-10-producing monocytes. These results highlight a novel role for androgen-driven peripheral IL-10-producing monocytes in the sexual dimorphism of pain resolution. These findings add to the growing concept that immune cells play a critical role in resolving pain and preventing the transition into chronic pain.

Graphical abstract

Expansion of extrafollicular B and T cell subsets in childhood-onset systemic lupus erythematosus

Introduction

Most childhood-onset SLE patients (cSLE) develop lupus nephritis (cLN), but only a small proportion achieve complete response to current therapies. The prognosis of children with LN and end-stage renal disease is particularly dire. Mortality rates within the first five years of renal replacement therapy may reach 22%. Thus, there is urgent need to decipher and target immune mechanisms that drive cLN. Despite the clear role of autoantibody production in SLE, targeted B cell therapies such as rituximab (anti-CD20) and belimumab (anti-BAFF) have shown only modest efficacy in cLN. While many studies have linked dysregulation of germinal center formation to SLE pathogenesis, other work supports a role for extrafollicular B cell activation in generation of pathogenic antibody secreting cells. However, whether extrafollicular B cell subsets and their T cell collaborators play a role in specific organ involvement in cLN and/or track with disease activity remains unknown.

Methods

We analyzed high-dimensional mass cytometry and gene expression data from 24 treatment naïve cSLE patients at the time of diagnosis and longitudinally, applying novel computational tools to identify abnormalities associated with clinical manifestations (cLN) and disease activity (SLEDAI).

Results

cSLE patients have an extrafollicular B cell expansion signature, with increased frequency of i) DN2, ii) Bnd2, iii) plasmablasts, and iv) peripheral T helper cells. Most importantly, we discovered that this extrafollicular signature correlates with disease activity in cLN, supporting extrafollicular T/B interactions as a mechanism underlying pediatric renal pathogenesis.

Discussion

This study integrates established and emerging themes of extrafollicular B cell involvement in SLE by providing evidence for extrafollicular B and peripheral T helper cell expansion, along with elevated type 1 IFN activation, in a homogeneous cohort of treatment-naïve cSLE patients, a point at which they should display the most extreme state of their immune dysregulation.

A systematic analysis of the human immune response to Plasmodium vivax

BACKGROUNDThe biology of Plasmodium vivax is markedly different from that of P. falciparum; how this shapes the immune response to infection remains unclear. To address this shortfall, we inoculated human volunteers with a clonal field isolate of P. vivax and tracked their response through infection and convalescence.METHODSParticipants were injected intravenously with blood-stage parasites and infection dynamics were tracked in real time by quantitative PCR. Whole blood samples were used for high dimensional protein analysis, RNA sequencing, and cytometry by time of flight, and temporal changes in the host response to P. vivax were quantified by linear regression. Comparative analyses with P. falciparum were then undertaken using analogous data sets derived from prior controlled human malaria infection studies.RESULTSP. vivax rapidly induced a type I inflammatory response that coincided with hallmark features of clinical malaria. This acute-phase response shared remarkable overlap with that induced by P. falciparum but was significantly elevated (at RNA and protein levels), leading to an increased incidence of pyrexia. In contrast, T cell activation and terminal differentiation were significantly increased in volunteers infected with P. falciparum. Heterogeneous CD4+ T cells were found to dominate this adaptive response and phenotypic analysis revealed unexpected features normally associated with cytotoxicity and autoinflammatory disease.CONCLUSIONP. vivax triggers increased systemic interferon signaling (cf P. falciparum), which likely explains its reduced pyrogenic threshold. In contrast, P. falciparum drives T cell activation far in excess of P. vivax, which may partially explain why falciparum malaria more frequently causes severe disease.TRIAL REGISTRATIONClinicalTrials.gov NCT03797989.FUNDINGThe European Union's Horizon 2020 Research and Innovation programme, the Wellcome Trust, and the Royal Society.

Inherited human ezrin deficiency impairs adaptive immunity

BACKGROUND

Inborn errors of immunity (IEI) are a group of monogenic diseases that confer susceptibility to infection, autoimmunity, and cancer. Despite the life-threatening consequences of some IEI, their genetic cause remains unknown in many patients.

OBJECTIVE

We investigated a patient with an IEI of unknown genetic etiology.

METHODS

Whole-exome sequencing identified a homozygous missense mutation of the gene encoding ezrin (EZR), substituting a threonine for an alanine at position 129.

RESULTS

Ezrin is one of the subunits of the ezrin, radixin, and moesin (ERM) complex. The ERM complex links the plasma membrane to the cytoskeleton and is crucial for the assembly of an efficient immune response. The A129T mutation abolishes basal phosphorylation and decreases calcium signaling, leading to complete loss of function. Consistent with the pleiotropic function of ezrin in myriad immune cells, multidimensional immunophenotyping by mass and flow cytometry revealed that in addition to hypogammaglobulinemia, the patient had low frequencies of switched memory B cells, CD4+ and CD8+ T cells, MAIT, γδ T cells, and centralnaive CD4+ cells.

CONCLUSIONS

Autosomal-recessive human ezrin deficiency is a newly recognized genetic cause of B-cell deficiency affecting cellular and humoral immunity.

Akkermansia muciniphila exacerbates food allergy in fibre-deprived mice

Alterations in the gut microbiome, including diet-driven changes, are linked to the rising prevalence of food allergy. However, little is known about how specific gut bacteria trigger the breakdown of oral tolerance. Here we show that depriving specific-pathogen-free mice of dietary fibre leads to a gut microbiota signature with increases in the mucin-degrading bacterium Akkermansia muciniphila. This signature is associated with intestinal barrier dysfunction, increased expression of type 1 and 2 cytokines and IgE-coated commensals in the colon, which result in an exacerbated allergic reaction to food allergens, ovalbumin and peanut. To demonstrate the causal role of A. muciniphila, we employed a tractable synthetic human gut microbiota in gnotobiotic mice. The presence of A. muciniphila within the microbiota, combined with fibre deprivation, resulted in stronger anti-commensal IgE coating and innate type-2 immune responses, which worsened symptoms of food allergy. Our study provides important insights into how gut microbes can regulate immune pathways of food allergy in a diet-dependent manner.

Immunomodulation by subcutaneously injected methacrylic acid-based hydrogels and tolerogenic dendritic cells in a mouse model of autoimmune diabetes

Type 1 diabetes is an autoimmune disease associated with the destruction of insulin-producing β cells. Immunotherapies are being developed to mitigate autoimmune diabetes. One promising option is the delivery of tolerogenic dendritic cells (DCs) primed with specific β-cell-associated autoantigens. These DCs can combat autoreactive cells and promote expansion of β-cell-specific regulatory immune cells, including Tregs. Tolerogenic DCs are typically injected systemically (or near target lymph nodes) in suspension, precluding control over the microenvironment surrounding tolerogenic DC interactions with the host. In this study we show that degradable, synthetic methacrylic acid (MAA)-based hydrogels are an inherently immunomodulating delivery vehicle that enhances tolerogenic DC therapy in the context of autoimmune diabetes. MAA hydrogels were found to affect the local recruitment and activation state of macrophages, DCs, T cells and other cells. Delivering tolerogenic DCs in the MAA hydrogel improved the local host response (e.g., fewer cytotoxic T cells) and enhanced peripheral Treg expansion. Non obese diabetic (NOD) mice treated with tolerogenic DCs subcutaneously injected in MAA hydrogels showed a delay in onset of autoimmune diabetes compared to control vehicles. Our findings further demonstrate the usefulness of MAA-based hydrogels as platforms for regenerative medicine in the context of type 1 diabetes.

High-dimensional analysis of T-cell profiling variations following belimumab treatment in systemic lupus erythematosus

OBJECTIVE

This study sought to elucidate the molecular impacts of belimumab (BEL) treatment on T-cell immune profiling in SLE.

METHODS

We used mass cytometry with 25 marker panels for T-cell immune profiling in peripheral blood T cells (CD3+) from 22 patients with BEL-treated SLE and 20 controls with non-BEL-treated SLE. An unsupervised machine-learning clustering, FlowSOM, was used to identify 39 T-cell clusters (TCLs; TCL01-TCL39). TCLs (% of CD3+) showing significant (p<0.05) associations with BEL treatment (BEL-TCL) were selected by a linear mixed-effects model for comparing groups of time-series data. Furthermore, we analysed the association between BEL treatment and variations in regulatory T-cell (Treg) phenotypes, and the ratio of other T-cell subsets to Treg as secondary analysis.

RESULTS

Clinical outcomes: BEL treatment was associated with a decrease in daily prednisolone use (coef=-0.1769, p=0.00074), and an increase in serum CH50 (coef=0.4653, p=0.003), C3 (coef=1.1047, p=0.00001) and C4 (coef=0.2990, p=0.00157) levels. Molecular effects: five distinct BEL-TCLs (TCL 04, 07, 11, 12 and 27) were identified. Among these, BEL-treated patients exhibited increased proportions in the Treg-like cluster TCL11 (coef=0.404, p=0.037) and two naïve TCLs (TCL04 and TCL07). TCL27 showed increased levels (coef=0.222, p=0.037) inversely correlating with baseline C3 levels. Secondary analyses revealed associations between BEL treatment and an increase in Tregs (coef=1.749, p=0.0044), elevated proportions of the fraction of Tregs with inhibitory function (fTregs, coef=0.7294, p=0.0178) and changes in peripheral helper T cells/fTreg (coef=-4.475, p=0.0319) and T helper 17/fTreg ratios (coef=-6.7868, p=0.0327). Additionally, BEL was linked to variations in T-cell immunoglobulin and mucin domain-containing protein-3 expression (coef=0.2422, p=0.039).

CONCLUSIONS

The study suggests an association between BEL treatment and variations in T cells, particularly Tregs, in SLE pathologies involving various immune cells.

GMM-Based Expanded Feature Space as a Way to Extract Useful Information for Rare Cell Subtypes Identification in Single-Cell Mass Cytometry

Cell subtype identification from mass cytometry data presents a persisting challenge, particularly when dealing with millions of cells. Current solutions are consistently under development, however, their accuracy and sensitivity remain limited, particularly in rare cell-type detection due to frequent downsampling. Additionally, they often lack the capability to analyze large data sets. To overcome these limitations, a new method was suggested to define an extended feature space. When combined with the robust clustering algorithm for big data, it results in more efficient cell clustering. Each marker's intensity distribution is presented as a mixture of normal distributions (Gaussian Mixture Model, GMM), and the expanded space is created by spanning over all obtained GMM components. The projection of the initial flow cytometry marker domain into the expanded space employs GMM-based membership functions. An evaluation conducted on three established cellular identification algorithms (FlowSOM, ClusterX, and PARC) utilizing the most substantial publicly available annotated dataset by Samusik et al. demonstrated the superior performance of the suggested approach in comparison to the standard. Although our approach identified 20 cell clusters instead of the expected 24, their intra-cluster homogeneity and inter-cluster differences were superior to the 24-cluster FlowSOM-based solution.

Dendritic Cell Vaccination in Conjunction with a TLR Agonist Polarizes Interferon Immune Responses in Malignant Glioma Patients

Autologous tumor lysate-pulsed dendritic cell (ATL-DC) vaccination is a promising immunotherapy for patients with high grade gliomas, but responses have not been demonstrated in all patients. To determine the most effective combination of autologous tumor lysate-pulsed DC vaccination, with or without the adjuvant toll-like receptor (TLR) agonists poly-ICLC or resiquimod, we conducted a Phase 2 clinical trial in 23 patients with newly diagnosed or recurrent WHO Grade III-IV malignant gliomas. We then performed deep, high-dimensional immune profiling of these patients to better understand how TLR agonists may influence the systemic immune responses induced by ATL-DC vaccination. Bulk RNAseq data demonstrated highly significant upregulation of type 1 and type 2 interferon gene expression selectively in patients who received adjuvant a TLR agonist together with ATL-DC. CyTOF analysis of patient peripheral blood mononuclear cells (PBMCs) showed increased expression of PD-1 on CD4+ T-cells, decreases in CD38 and CD39 on CD8+ T cells and elevated proportion of monocytes after ATL-DC + TLR agonist administration. In addition, scRNA-seq demonstrated a higher expression fold change of IFN-induced genes with poly-ICLC treatment in both peripheral blood monocytes and T lymphocytes. Patients who had higher expression of interferon response genes lived significantly longer and had longer time to progression compared to those with lower expression. The results suggest that ATL-DC in conjunction with adjuvant poly-ICLC induces a polarized interferon response in circulating monocytes and specific activation of a CD8+ T cell population, which may represent an important blood biomarker for immunotherapy in this patient population.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT01204684.

Early Infiltration of Innate Immune Cells to the Liver Depletes HNF4α and Promotes Extrahepatic Carcinogenesis

Multiple studies have identified metabolic changes within the tumor and its microenvironment during carcinogenesis. Yet, the mechanisms by which tumors affect the host metabolism are unclear. We find that systemic inflammation induced by cancer leads to liver infiltration of myeloid cells at early extrahepatic carcinogenesis. The infiltrating immune cells via IL6-pSTAT3 immune-hepatocyte cross-talk cause the depletion of a master metabolic regulator, HNF4α, consequently leading to systemic metabolic changes that promote breast and pancreatic cancer proliferation and a worse outcome. Preserving HNF4α levels maintains liver metabolism and restricts carcinogenesis. Standard liver biochemical tests can identify early metabolic changes and predict patients' outcomes and weight loss. Thus, the tumor induces early metabolic changes in its macroenvironment with diagnostic and potentially therapeutic implications for the host.

SIGNIFICANCE

Cancer growth requires a permanent nutrient supply starting from early disease stages. We find that the tumor extends its effect to the host's liver to obtain nutrients and rewires the systemic and tissue-specific metabolism early during carcinogenesis. Preserving liver metabolism restricts tumor growth and improves cancer outcomes. This article is highlighted in the In This Issue feature, p. 1501.