SuperCellCyto: enabling efficient analysis of large scale cytometry datasets.

Advancements in cytometry technologies have enabled quantification of up to 50 proteins across millions of cells at single cell resolution. Analysis of cytometry data routinely involves tasks such as data integration, clustering, and dimensionality reduction. While numerous tools exist, many require extensive run times when processing large cytometry data containing millions of cells. Existing solutions, such as random subsampling, are inadequate as they risk excluding rare cell subsets. To address this, we propose SuperCellCyto, an R package that builds on the SuperCell tool which groups highly similar cells into supercells. SuperCellCyto is available on GitHub ( https://github.com/phipsonlab/SuperCellCyto ) and Zenodo ( https://doi.org/10.5281/zenodo.10521294 ).

CD8+ T Cells Mediate Lethal Lung Pathology in the Absence of PD-L1 and Type I Interferon Signalling following LCMV Infection.

CD8+ T cells are critical to the adaptive immune response against viral pathogens. However, overwhelming antigen exposure can result in their exhaustion, characterised by reduced effector function, failure to clear virus, and the upregulation of inhibitory receptors, including programmed cell death 1 (PD-1). However, exhausted T cell responses can be "re-invigorated" by inhibiting PD-1 or the primary ligand of PD-1: PD-L1. Further, the absence of the type I interferon receptor IFNAR1 also results in T cell exhaustion and virus persistence in lymphocytic choriomeningitis virus Armstrong (LCMV-Arm)-infected mice. In this study, utilizing single- and double-knockout mice, we aimed to determine whether ablation of PD-1 could restore T cell functionality in the absence of IFNAR1 signalling in LCMV-Arm-infected mice. Surprisingly, this did not re-invigorate the T cell response and instead, it converted chronic LCMV-Arm infection into a lethal disease characterized by severe lung inflammation with an infiltration of neutrophils and T cells. Depletion of CD8+ T cells, but not neutrophils, rescued mice from lethal disease, demonstrating that IFNAR1 is required to prevent T cell exhaustion and virus persistence in LCMV-Arm infection, and in the absence of IFNAR1, PD-L1 is required for survival. This reveals an important interplay between IFNAR1 and PD-L1 with implications for therapeutics targeting these pathways.

Clodronate is not protective in lethal viral encephalitis despite substantially reducing inflammatory monocyte infiltration in the CNS.

Bone marrow (BM)-derived monocytes induce inflammation and tissue damage in a range of pathologies. In particular, in a mouse model of West Nile virus (WNV) encephalitis (WNE), nitric oxide-producing, Ly6Chi inflammatory monocytes from the BM are recruited to the central nervous system (CNS) and contribute to lethal immune pathology. Reducing the migration of these cells into the CNS using monoclonal antibody blockade, immune-modifying particles or CSF-1R inhibitors reduces neuroinflammation, improving survival and/or clinical outcomes. Macrophages can also be targeted more broadly by administration of clodronate-encapsulated liposomes, which induce apoptosis in phagocytes. In this study, clodronate reduced the inflammatory infiltrate by 70% in WNE, however, surprisingly, this had no effect on disease outcome. More detailed analysis demonstrated a compensatory increase in neutrophils and enhanced activation status of microglia in the brain. In addition, we observed increased numbers of Ly6Chi BM monocytes with an increased proliferative capacity and expression of SCA-1 and CD16/32, potentially indicating output of immature cells from the BM. Once in the brain, these cells were more phagocytic and had a reduced expression of antigen-presenting molecules. Lastly, we show that clodronate also reduces non-myeloid cells in the spleen and BM, as well as ablating red blood cells and their proliferation. These factors likely impeded the therapeutic potential of clodronate in WNE. Thus, while clodronate provides an excellent system to deplete macrophages in the body, it has larger and broader effects on the phagocytic and non-phagocytic system, which must be considered in the interpretation of data.

A unique cytotoxic CD4+ T cell-signature defines critical COVID-19.

Objectives: SARS-CoV-2 infection causes a spectrum of clinical disease presentation, ranging from asymptomatic to fatal. While neutralising antibody (NAb) responses correlate with protection against symptomatic and severe infection, the contribution of the T-cell response to disease resolution or progression is still unclear. As newly emerging variants of concern have the capacity to partially escape NAb responses, defining the contribution of individual T-cell subsets to disease outcome is imperative to inform the development of next-generation COVID-19 vaccines. Methods: Immunophenotyping of T-cell responses in unvaccinated individuals was performed, representing the full spectrum of COVID-19 clinical presentation. Computational and manual analyses were used to identify T-cell populations associated with distinct disease states. Results: Critical SARS-CoV-2 infection was characterised by an increase in activated and cytotoxic CD4+ lymphocytes (CTL). These CD4+ CTLs were largely absent in asymptomatic to severe disease states. In contrast, non-critical COVID-19 was associated with high frequencies of naïve T cells and lack of activation marker expression. Conclusion: Highly activated and cytotoxic CD4+ T-cell responses may contribute to cell-mediated host tissue damage and progression of COVID-19. Induction of these potentially detrimental T-cell responses should be considered when developing and implementing effective COVID-19 control strategies.

Circulating effector γδ T cell populations are associated with acute coronavirus disease 19 in unvaccinated individuals.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection causes severe coronavirus disease 2019 (COVID-19) in a small proportion of infected individuals. The immune system plays an important role in the defense against SARS-CoV-2, but our understanding of the cellular immune parameters that contribute to severe COVID-19 disease is incomplete. Here, we show that populations of effector γδ T cells are associated with COVID-19 in unvaccinated patients with acute disease. We found that circulating CD27neg CD45RA+ CX3CR1+ Vδ1effector cells expressing Granzymes (Gzms) were enriched in COVID-19 patients with acute disease. Moreover, higher frequencies of GzmB+ Vδ2+ T cells were observed in acute COVID-19 patients. SARS-CoV-2 infection did not alter the γδ T cell receptor repertoire of either Vδ1+ or Vδ2+ subsets. Our work demonstrates an association between effector populations of γδ T cells and acute COVID-19 in unvaccinated individuals.

TrackSOM: Mapping immune response dynamics through clustering of time-course cytometry data.

Mapping the dynamics of immune cell populations over time or disease-course is key to understanding immunopathogenesis and devising putative interventions. We present TrackSOM, a novel method for delineating cellular populations and tracking their development over a time- or disease-course cytometry datasets. We demonstrate TrackSOM-enabled elucidation of the immune response to West Nile Virus infection in mice, uncovering heterogeneous subpopulations of immune cells and relating their functional evolution to disease severity. TrackSOM is easy to use, encompasses few parameters, is quick to execute, and enables an integrative and dynamic overview of the immune system kinetics that underlie disease progression and/or resolution.

A high-dimensional cytometry atlas of peripheral blood over the human life span.

Age can profoundly affect susceptibility to a broad range of human diseases. Children are more susceptible to some infectious diseases such as diphtheria and pertussis, while in others, such as coronavirus disease 2019 and hepatitis A, they are more protected compared with adults. One explanation is that the composition of the immune system is a major contributing factor to disease susceptibility and severity. While most studies of the human immune system have focused on adults, how the immune system changes after birth remains poorly understood. Here, using high-dimensional spectral flow cytometry and computational methods for data integration, we analyzed more than 50 populations of immune cells in the peripheral blood, generating an immune cell atlas that defines the healthy human immune system from birth up to 75 years of age. We focused our efforts on children under 18 years old, revealing major changes in immune cell populations after birth and in children of schooling age. Specifically, CD4+ T effector memory cells, Vδ2+ gamma delta (γδ)T cells, memory B cells, plasmablasts, CD11c+ B cells and CD16+ CD56bright natural killer (NK) cells peaked in children aged 5-9 years old, whereas frequencies of T helper 1, T helper 17, dendritic cells and CD16+ CD57+ CD56dim NK cells were highest in older children (10-18 years old). The frequency of mucosal-associated invariant T cells was low in the first several years of life and highest in adults between 19 and 30 years old. Late adulthood was associated with fewer mucosal-associated invariant T cells and Vδ2+ γδ T cells but with increased frequencies of memory subsets of B cells, CD4+ and CD8+ T cells and CD57+ NK cells. This human immune cell atlas provides a critical resource to understand changes to the immune system during life and provides a reference for investigating the immune system in the context of human disease. This work may also help guide future therapies that target specific populations of immune cells to protect at-risk populations.

SARS-CoV-2 infection results in immune responses in the respiratory tract and peripheral blood that suggest mechanisms of disease severity.

Respiratory tract infection with SARS-CoV-2 results in varying immunopathology underlying COVID-19. We examine cellular, humoral and cytokine responses covering 382 immune components in longitudinal blood and respiratory samples from hospitalized COVID-19 patients. SARS-CoV-2-specific IgM, IgG, IgA are detected in respiratory tract and blood, however, receptor-binding domain (RBD)-specific IgM and IgG seroconversion is enhanced in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples correlates with RBD-specific IgM and IgG levels. Cytokines/chemokines vary between respiratory samples and plasma, indicating that inflammation should be assessed in respiratory specimens to understand immunopathology. IFN-α2 and IL-12p70 in endotracheal aspirate and neutralization in sputum negatively correlate with duration of hospital stay. Diverse immune subsets are detected in respiratory samples, dominated by neutrophils. Importantly, dexamethasone treatment does not affect humoral responses in blood of COVID-19 patients. Our study unveils differential immune responses between respiratory samples and blood, and shows how drug therapy affects immune responses during COVID-19.

The cytokines interleukin-6 and interferon-α induce distinct microglia phenotypes.

BACKGROUND: Elevated production of the cytokines interleukin (IL)-6 or interferon (IFN)-α in the central nervous system (CNS) is implicated in the pathogenesis of neurological diseases such as neuromyelitis optica spectrum disorders or cerebral interferonopathies, respectively. Transgenic mice with CNS-targeted chronic production of IL-6 (GFAP-IL6) or IFN-α (GFAP-IFN) recapitulate important clinical and pathological features of these human diseases. The activation of microglia is a prominent manifestation found both in the human diseases and in the transgenic mice, yet little is known about how this contributes to disease pathology. METHODS: Here, we used a combination of ex vivo and in situ techniques to characterize the molecular, cellular and transcriptomic phenotypes of microglia in GFAP-IL6 versus GFAP-IFN mice. In addition, a transcriptomic meta-analysis was performed to compare the microglia response from GFAP-IL6 and GFAP-IFN mice to the response of microglia in a range of neurodegenerative and neuroinflammatory disorders. RESULTS: We demonstrated that microglia show stimulus-specific responses to IL-6 versus IFN-α in the brain resulting in unique and extensive molecular and cellular adaptations. In GFAP-IL6 mice, microglia proliferated, had shortened, less branched processes and elicited transcriptomic and molecular changes associated with phagocytosis and lipid processing. In comparison, microglia in the brain of GFAP-IFN mice exhibited increased proliferation and apoptosis, had larger, hyper-ramified processes and showed transcriptomic and surface marker changes associated with antigen presentation and antiviral response. Further, a transcriptomic meta-analysis revealed that IL-6 and IFN-α both contribute to the formation of a core microglia response in animal models of neurodegenerative and neuroinflammatory disorders, such as Alzheimer's disease, tauopathy, multiple sclerosis and lipopolysaccharide-induced endotoxemia. CONCLUSIONS: Our findings demonstrate that microglia responses to IL-6 and IFN-α are highly stimulus-specific, wide-ranging and give rise to divergent phenotypes that modulate microglia responses in neuroinflammatory and neurodegenerative diseases.

Peripheral B-cell dysregulation is associated with relapse after long-term quiescence in patients with multiple sclerosis.

B cells play a major role in multiple sclerosis (MS), with many successful therapeutics capable of removing them from circulation. One such therapy, alemtuzumab, is thought to reset the immune system without the need for ongoing therapy in a proportion of patients. The exact cells contributing to disease pathogenesis and quiescence remain to be identified. We utilized mass cytometry to analyze B cells from the blood of patients with relapse-remitting MS (RRMS) before and after alemtuzumab treatment, and during relapse. A complementary RRMS cohort was analyzed by single-cell RNA sequencing. The R package "Spectre" was used to analyze these data, incorporating FlowSOM clustering, sparse partial least squares-discriminant analysis and permutational multivariate analysis of variance. Immunoglobulin (Ig)A+ and IgG1 + B-cell numbers were altered, including higher IgG1 + B cells during relapse. B-cell linker protein (BLNK), CD40 and CD210 expression by B cells was lower in patients with RRMS compared with non-MS controls, with similar results at the transcriptomic level. Finally, alemtuzumab restored BLNK, CD40 and CD210 expression by IgA+ and IgG1 + B cells, which was altered again during relapse. These data suggest that impairment of IgA+ and IgG1 + B cells may contribute to MS pathogenesis, which can be restored by alemtuzumab.

High-parameter cytometry unmasks microglial cell spatio-temporal response kinetics in severe neuroinflammatory disease.

BACKGROUND: Differentiating infiltrating myeloid cells from resident microglia in neuroinflammatory disease is challenging, because bone marrow-derived inflammatory monocytes infiltrating the inflamed brain adopt a 'microglia-like' phenotype. This precludes the accurate identification of either cell type without genetic manipulation, which is important to understand their temporal contribution to disease and inform effective intervention in its pathogenesis. During West Nile virus (WNV) encephalitis, widespread neuronal infection drives substantial CNS infiltration of inflammatory monocytes, causing severe immunopathology and/or death, but the role of microglia in this remains unclear. METHODS: Using high-parameter cytometry and dimensionality-reduction, we devised a simple, novel gating strategy to identify microglia and infiltrating myeloid cells during WNV-infection. Validating our strategy, we (1) blocked the entry of infiltrating myeloid populations from peripheral blood using monoclonal blocking antibodies, (2) adoptively transferred BM-derived monocytes and tracked their phenotypic changes after infiltration and (3) labelled peripheral leukocytes that infiltrate into the brain with an intravenous dye. We demonstrated that myeloid immigrants populated only the identified macrophage gates, while PLX5622 depletion reduced all 4 subsets defined by the microglial gates. RESULTS: Using this gating approach, we identified four consistent microglia subsets in the homeostatic and WNV-infected brain. These were P2RY12hi CD86-, P2RY12hi CD86+ and P2RY12lo CD86- P2RY12lo CD86+. During infection, 2 further populations were identified as 'inflammatory' and 'microglia-like' macrophages, recruited from the bone marrow. Detailed kinetic analysis showed significant increases in the proportions of both P2RY12lo microglia subsets in all anatomical areas, largely at the expense of the P2RY12hi CD86- subset, with the latter undergoing compensatory proliferation, suggesting replenishment of, and differentiation from this subset in response to infection. Microglia altered their morphology early in infection, with all cells adopting temporal and regional disease-specific phenotypes. Late in disease, microglia produced IL-12, downregulated CX3CR1, F4/80 and TMEM119 and underwent apoptosis. Infiltrating macrophages expressed both TMEM119 and P2RY12 de novo, with the microglia-like subset notably exhibiting the highest proportional myeloid population death. CONCLUSIONS: Our approach enables detailed kinetic analysis of resident vs infiltrating myeloid cells in a wide range of neuroinflammatory models without non-physiological manipulation. This will more clearly inform potential therapeutic approaches that specifically modulate these cells.

Intrapulmonary vaccination with delta-inulin adjuvant stimulates non-polarised chemotactic signalling and diverse cellular interaction.

There is an urgent need for novel vaccination strategies to combat respiratory pathogens. Mucosal vaccine delivery is an attractive option as it directly targets the site of infection; however, preclinical development has been hindered by a lack of suitable mucosal adjuvants and a limited understanding of their immune effects in the lung environment. Herein, we define the early immune events following the intrapulmonary delivery of a vaccine incorporating the adjuvant delta-inulin. Analysis of the early inflammatory response showed vaccine-induced innate cell recruitment to lungs and local lymph nodes (LN) was transient and non-polarised, correlating with an increase in pulmonary chemotactic factors. Use of fluorescently labelled adjuvant revealed widespread tissue dissemination of adjuvant particles, coupled with broad cellular uptake and transit to the lung-draining LN by a range of innate immune cells. Mass cytometric analysis revealed extensive phenotypic changes in innate and adaptive cell subsets induced by vaccination; this included identification of unconventional lymphocytes such as γδ-T cells and MAIT cells that increased following vaccination and displayed an activated phenotype. This study details a comprehensive view of the immune response to intrapulmonary adjuvant administration and provides pre-clinical evidence to support delta-inulin as a suitable adjuvant for pulmonary vaccines.

Integrated immune dynamics define correlates of COVID-19 severity and antibody responses.

SARS-CoV-2 causes a spectrum of COVID-19 disease, the immunological basis of which remains ill defined. We analyzed 85 SARS-CoV-2-infected individuals at acute and/or convalescent time points, up to 102 days after symptom onset, quantifying 184 immunological parameters. Acute COVID-19 presented with high levels of IL-6, IL-18, and IL-10 and broad activation marked by the upregulation of CD38 on innate and adaptive lymphocytes and myeloid cells. Importantly, activated CXCR3+cTFH1 cells in acute COVID-19 significantly correlate with and predict antibody levels and their avidity at convalescence as well as acute neutralization activity. Strikingly, intensive care unit (ICU) patients with severe COVID-19 display higher levels of soluble IL-6, IL-6R, and IL-18, and hyperactivation of innate, adaptive, and myeloid compartments than patients with moderate disease. Our analyses provide a comprehensive map of longitudinal immunological responses in COVID-19 patients and integrate key cellular pathways of complex immune networks underpinning severe COVID-19, providing important insights into potential biomarkers and immunotherapies.

Using single-cell cytometry to illustrate integrated multi-perspective evaluation of clustering algorithms using Pareto fronts.

MOTIVATION: Many 'automated gating' algorithms now exist to cluster cytometry and single cell sequencing data into discrete populations. Comparative algorithm evaluations on benchmark datasets rely either on a single performance metric, or a few metrics considered independently of one another. However, single metrics emphasise different aspects of clustering performance and do not rank clustering solutions in the same order. This underlies the lack of consensus between comparative studies regarding optimal clustering algorithms and undermines the translatability of results onto other non-benchmark datasets. RESULTS: We propose the Pareto fronts framework as an integrative evaluation protocol, wherein individual metrics are instead leveraged as complementary perspectives. Judged superior are algorithms that provide the best trade-off between the multiple metrics considered simultaneously. This yields a more comprehensive and complete view of clustering performance. Moreover, by broadly and systematically sampling algorithm parameter values using the Latin Hypercube sampling method, our evaluation protocol minimises (un)fortunate parameter value selections as confounding factors. Furthermore, it reveals how meticulously each algorithm must be tuned in order to obtain good results, vital knowledge for users with novel data. We exemplify the protocol by conducting a comparative study between three clustering algorithms (ChronoClust, FlowSOM and Phenograph) using four common performance metrics applied across four cytometry benchmark datasets. To our knowledge, this is the first time Pareto fronts have been used to evaluate the performance of clustering algorithms in any application domain. AVAILABILITY: Implementation of our Pareto front methodology and all scripts to reproduce this article are available at https://github.com/ghar1821/ParetoBench.

Contribution of STAT1 to innate and adaptive immunity during type I interferon-mediated lethal virus infection.

Signal transducers and activators of transcription (STAT) 1 is critical for cellular responses to type I interferons (IFN-Is), with the capacity to determine the outcome of viral infection. We previously showed that while wildtype (WT) mice develop mild disease and survive infection with lymphocytic choriomeningitis virus (LCMV), LCMV infection of STAT1-deficient mice results in a lethal wasting disease that is dependent on IFN-I and CD4+ cells. IFN-Is are considered to act as a bridge between innate and adaptive immunity. Here, we determined the relative contribution of STAT1 on innate and adaptive immunity during LCMV infection. We show that STAT1 deficiency results in a biphasic disease following LCMV infection. The initial, innate immunity-driven phase of disease was characterized by rapid weight loss, thrombocytopenia, systemic cytokine and chemokine responses and leukocyte infiltration of infected organs. In the absence of an adaptive immune response, this first phase of disease largely resolved resulting in survival of the infected host. However, in the presence of adaptive immunity, the disease progressed into a second phase with continued cytokine and chemokine production, persistent leukocyte extravasation into infected tissues and ultimately, host death. Overall, our findings demonstrate the key contribution of STAT1 in modulating innate and adaptive immunity during type I interferon-mediated lethal virus infection.

Zika virus encephalitis in immunocompetent mice is dominated by innate immune cells and does not require T or B cells.

BACKGROUND: Until the end of the twentieth century, Zika virus (ZIKV) was thought to cause a mostly mild, self-limiting disease in humans. However, as the geographic distribution of ZIKV has shifted, so too has its pathogenicity. Modern-day ZIKV infection is now known to cause encephalitis, acute disseminated encephalomyelitis, and Guillain-Barré syndrome in otherwise healthy adults. Nevertheless, the underlying pathogenetic mechanisms responsible for this shift in virulence remain unclear. METHODS: Here, we investigated the contribution of the innate versus the adaptive immune response using a new mouse model involving intracranial infection of adult immunocompetent mice with a moderately low dose of ZIKV MR766. To determine the contribution of type I interferons (IFN-Is) and adaptive immune cells, we also studied mice deficient for the IFN-I receptor 1 (Ifnar1-/-) and recombination-activating gene 1 (Rag1-/-). RESULTS: We show that intracranial infection with ZIKV resulted in lethal encephalitis. In wild-type mice, ZIKV remained restricted predominantly to the central nervous system (CNS) and infected neurons, whereas astrocytes and microglia were spared. Histological and molecular analysis revealed prominent activation of resident microglia and infiltrating monocytes that were accompanied by an expression of pro-inflammatory cytokines. The disease was independent of T and B cells. Importantly, unlike peripheral infection, IFN-Is modulated but did not protect from infection and lethal disease. Lack of IFN-I signaling resulted in spread of the virus, generalized inflammatory changes, and accelerated disease onset. CONCLUSIONS: Using intracranial infection of immunocompetent wild-type mice with ZIKV, we demonstrate that in contrast to the peripheral immune system, the CNS is susceptible to infection and responds to ZIKV by initiating an antiviral immune response. This response is dominated by resident microglia and infiltrating monocytes and macrophages but does not require T or B cells. Unlike in the periphery, IFN-Is in the CNS cannot prevent the establishment of infection. Our findings show that ZIKV encephalitis in mice is dependent on the innate immune response, and adaptive immune cells play at most a minor role in disease pathogenesis.

Staining of Phosphorylated Signalling Markers Protocol for Mass Cytometry.

Mass cytometry is a multi-parametric technique that offers insight into functional and biological systems at a single cell level (Tanner et al., Cancer Immunol Immunother 62:955-965, 2013). One of the major advantages of mass cytometry is the ability to measure multiple intracellular markers, including phosphorylated proteins that are part of major signaling pathways, such as NF-κB, JAK/STAT, and ERK/MAPK. Here we describe an optimized mass cytometry protocol for staining human clinical blood samples with panels that include phosphorylated antibodies.

Analysis of the Murine Bone Marrow Hematopoietic System Using Mass and Flow Cytometry.

The hematopoietic system produces erythrocytes (red blood cells), leukocytes (white blood cells), and thrombocytes (platelets) throughout the life of an organism. Long-lived hematopoietic stem cells give rise to early progenitors with multi-lineage potential that progressively differentiate into lineage-specific progenitors. Following lineage commitment, these progenitors proliferate and expand, before eventually differentiating into their mature forms. This process drives the up- and downregulation of a wide variety of surface and intracellular markers throughout differentiation, making cytometric analysis of this interconnected system challenging. Moreover, during inflammation, the hematopoietic system can be mobilized to re-prioritize the production of various lineages, in order to match increased demand, often at the expense of other lineages. As such, the response of the hematopoietic system in the bone marrow (BM) is a critical component of both immunity and disease. Because of the complexity of the hematopoietic system in steady state and disease, high-dimensional cytometry technologies are well suited to the exploration of these complex systems. Here we describe a protocol for the extraction of murine bone marrow, and preparation for examination using high-dimensional flow or mass cytometry. Additionally, we describe methods for performing cell cycle assays using bromodeoxyuridine (BrdU) or iododeoxyuridine (IdU). Finally, we describe an analytical method that allows for a system-level analysis of the hematopoietic system in steady state or inflammatory scenarios.