Data-Driven Phenotypic Dissection of AML Reveals Progenitor-like Cells that Correlate with Prognosis

Identification of clinically relevant subsets CD39+PD-1+CD8+ T cells and CD39+ regulatory T cells in intrahepatic cholangiocarcinoma using single-cell CyTOF

Intrahepatic cholangiocarcinoma (iCCA) is an aggressive liver malignancy with limited treatment options and a dismal prognosis. The tumor immune microenvironment (TIME) is crucial for iCCA progression, yet its comprehensive characterization remains incomplete. This study utilized mass cytometry by time of flight (CyTOF) to comprehensively analyze immune cell populations in fresh iCCA tumor samples and adjacent peritumor liver tissues. Notably, NK cell percentages significantly decreased in iCCA lesions compared to peritumor liver tissues. Conversely, an enrichment of immunosuppressive CD39+Foxp3+CD4+ regulatory T cells (CD39+T-regs) and exhausted-like CD8+T cells (with pronounced CD39 and PD-1 expression) within TIME was identified and confirmed by multiplex immunofluorescence staining in an independent patient cohort (n = 140). Crucially, tumor-infiltrating CD39+T-regs and CD39+PD-1+CD8+T cells emerged as independent prognostic indicators associated with an unfavorable prognosis in iCCA. These findings unveil the intricate immune landscape within iCCA, offering valuable insights for disease management and novel cancer immunotherapies.

Intratumoral antigen signaling traps CD8+ T cells to confine exhaustion to the tumor site

Immunotherapy advances have been hindered by difficulties in tracking the behaviors of lymphocytes after antigen signaling. Here, we assessed the behavior of T cells active within tumors through the development of the antigen receptor signaling reporter (AgRSR) mouse, fate-mapping lymphocytes responding to antigens at specific times and locations. Contrary to reports describing the ready egress of T cells out of the tumor, we find that intratumoral antigen signaling traps CD8+ T cells in the tumor. These clonal populations expand and become increasingly exhausted over time. By contrast, antigen-signaled regulatory T cell (Treg) clonal populations readily recirculate out of the tumor. Consequently, intratumoral antigen signaling acts as a gatekeeper to compartmentalize CD8+ T cell responses, even within the same clonotype, thus enabling exhausted T cells to remain confined to a specific tumor tissue site.

Mapping cell diversity in human sporadic cerebral cavernous malformations

BACKGROUND

Cerebral cavernous malformation (CCM) is a low-flow, bleeding-prone vascular disease that can cause cerebral hemorrhage, seizure and neurological deficits. Its inheritance mode includes sporadic or autosomal dominant inheritance with incomplete penetrance, namely sporadic CCM (SCCM) and familial CCM. SCCM is featured by single lesion and single affection in a family. Among CCM patients especially SCCM, the pathogenesis of the corresponding phenotypes and pathological features or candidate genes have not been fully elucidated yet.

METHODS

Here, we performed in-depth single-cell RNA sequencing (scRNA-Seq) and bulk assay for transposase-accessible chromatin sequencing (ATAC-Seq) in SCCM and control patients. Further validation was conducted for the gene of interest using qPCR and RNA in situ hybridization (RNA FISH) techniques to provide further atlas and evidence for SCCM generative process.

RESULTS

We identified six cell types in the SCCM and control vessels and found that the expression of NEK1, RNPC3, FBRSL1, IQGAP2, MCUB, AP3B1, ESCO1, MYO9B and PVT1 were up-regulated in SCCM tissues. Among the six cell types, we found that compared with control conditions, PVT1 showed a rising peak which followed the pseudo-time axis in endothelial cell clusters of SCCM samples, while showed an increasing trend in smooth muscle cell clusters of SCCM samples. Further experiments indicated that, compared with the control vessels, PVT1 exhibited significantly elevated expression in SCCM samples.

CONCLUSION

In SCCM conditions, We found that in the process of development from control to lesion conditions, PVT1 showed a rising peak in endothelial cells and showed an increasing trend in smooth muscle cells at the same time. Overall, there was a significantly elevated expression of NEK1, RNPC3, FBRSL1, IQGAP2, MCUB, AP3B1, ESCO1, MYO9B and PVT1 in SCCM specimens compared to control samples.

SifiNet: a robust and accurate method to identify feature gene sets and annotate cells

SifiNet is a robust and accurate computational pipeline for identifying distinct gene sets, extracting and annotating cellular subpopulations, and elucidating intrinsic relationships among these subpopulations. Uniquely, SifiNet bypasses the cell clustering stage, commonly integrated into other cellular annotation pipelines, thereby circumventing potential inaccuracies in clustering that may compromise subsequent analyses. Consequently, SifiNet has demonstrated superior performance in multiple experimental datasets compared with other state-of-the-art methods. SifiNet can analyze both single-cell RNA and ATAC sequencing data, thereby rendering comprehensive multi-omic cellular profiles. It is conveniently available as an open-source R package.

Vδ1 Effector and Vδ2 γδ T-Cell Subsets Shift in Frequency and Are Linked to Plasma Inflammatory Markers During Antiretroviral Therapy-Suppressed HIV Infection

BACKGROUND

Chronic inflammation is prevalent with antiretroviral therapy (ART)-suppressed human immunodeficiency virus (HIV) infection and one immune cell subset putatively driving this phenomenon is TIGIT+ γδ T cells.

METHODS

To elucidate γδ T-cell phenotypic diversity, spectral flow cytometry was performed on blood lymphocytes from individuals of a HIV and aging cohort and data were analyzed using bioinformatic platforms. Plasma inflammatory markers were measured and correlated with γδ T-cell subset frequencies.

RESULTS

Thirty-nine distinct γδ T-cell subsets were identified (22 Vδ1+, 14 Vδ2+, and 3 Vδ1-Vδ2-Vγ9+) and TIGIT was nearly exclusively found on the Vδ1+CD45RA+CD27- effector populations. People with ART-suppressed HIV infection (PWH) exhibited high frequencies of distinct clusters of Vδ1+ effectors distinguished via CD8, CD16, and CD38 expression. Among Vδ2+ cells, most Vγ9+ (innate-like) clusters were lower in PWH; however, CD27+ subsets were similar in frequency between participants with and without HIV. Comparisons by age revealed lower 'naive' Vδ1+CD45RA+CD27+ cells in older individuals, regardless of HIV status. Plasma inflammatory markers were selectively linked to subsets of Vδ1+ and Vδ2+ cells.

CONCLUSIONS

These results further elucidate γδ T-cell subset complexity and reveal distinct alterations and connections with inflammatory pathways of Vδ1+ effector and Vδ2+ innate-like subsets during ART-suppressed HIV infection.

Exploring the prognostic significance of arm-level copy number alterations in triple-negative breast cancer

Somatic copy number alterations (SCNAs) are prevalent in cancer and play a significant role in both tumorigenesis and therapeutic resistance. While focal SCNAs have been extensively studied, the impact of larger arm-level SCNAs remains poorly understood. Here, we investigated the association between arm-level SCNAs and overall survival in triple-negative breast cancer (TNBC), an aggressive subtype of breast cancer lacking targeted therapies. We identified frequent arm-level SCNAs, including 21q gain and 7p gain, which correlated with poor overall survival in TNBC patients. Further, we identified the expression of specific genes within these SCNAs associated with survival. Notably, we found that the expression of RIPK4, a gene located on 21q, exhibited a strong correlation with poor overall survival. In functional assays, we demonstrated that targeting Ripk4 in a murine lung metastatic TNBC model significantly reduced tumor burden, improved survival, and increased CD4+ and CD8+ T cell infiltration. RIPK4 enhanced the survival of triple-negative breast cancer cells at secondary sites, thereby facilitating the formation of metastatic lesions. Our findings highlight the significance of arm-level SCNAs in breast cancer progression and identify RIPK4 as a putative driver of TNBC metastasis and immunosuppression.

PD-1 Receptor (+) T cells are associated with the efficacy of the combined treatment with regulatory t cells and rituximab in type 1 diabetes children via regulatory t cells suppressive activity amelioration

An imbalance between exaggerated autoaggressive T cell responses, primarily CD8 + T cells, and impaired tolerogenic mechanisms underlie the development of type 1 diabetes mellitus. Disease-modifying strategies, particularly immunotherapy focusing on FoxP3 + T regulatory cells (Treg), and B cells facilitating antigen presentation for T cells, show promise. Selective depletion of B cells may be achieved with an anti-CD20 monoclonal antibody (mAb). In a 2-year-long flow cytometry follow-up, involving 32 peripheral blood T and B cell markers across three trial arms (Treg + rituximab N = 12, Treg + placebo N = 13, control N = 11), we observed significant changes. PD-1 receptor (+) CD4 + Treg, CD4 + effector T cells (Teffs), and CD8 + T cell percentages increased in the combined regimen group by the end of follow-up. Conversely, the control group exhibited a notable reduction in PD-1 receptor (+) CD4 + Teff percentages. Considering clinical endpoints, higher PD-1 receptor (+) expression on T cells correlated with positive responses, including a higher mixed meal tolerance test AUC, and reduced daily insulin dosage. PD-1 receptor (+) T cells emerged as a potential therapy outcome biomarker. In vitro validation confirmed that successful Teff suppression was associated with elevated PD-1 receptor (+) Treg levels. These findings support PD-1 receptor (+) T cells as a reliable indicator of treatment with combined immunotherapy consisting of Tregs and anti-CD20 mAb efficacy in type 1 diabetes mellitus.

Discrete latent embedding of single-cell chromatin accessibility sequencing data for uncovering cell heterogeneity

Single-cell epigenomic data has been growing continuously at an unprecedented pace, but their characteristics such as high dimensionality and sparsity pose substantial challenges to downstream analysis. Although deep learning models-especially variational autoencoders-have been widely used to capture low-dimensional feature embeddings, the prevalent Gaussian assumption somewhat disagrees with real data, and these models tend to struggle to incorporate reference information from abundant cell atlases. Here we propose CASTLE, a deep generative model based on the vector-quantized variational autoencoder framework to extract discrete latent embeddings that interpretably characterize single-cell chromatin accessibility sequencing data. We validate the performance and robustness of CASTLE for accurate cell-type identification and reasonable visualization compared with state-of-the-art methods. We demonstrate the advantages of CASTLE for effective incorporation of existing massive reference datasets in a weakly supervised or supervised manner. We further demonstrate CASTLE's capacity for intuitively distilling cell-type-specific feature spectra that unveil cell heterogeneity and biological implications quantitatively.

Integrative spatial analysis reveals a multi-layered organization of glioblastoma

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Endothelial deletion of p53 generates transitional endothelial cells and improves lung development during neonatal hyperoxia

Bronchopulmonary dysplasia (BPD), a prevalent and chronic lung disease affecting premature newborns, results in vascular rarefaction and alveolar simplification. Although the vasculature has been recognized as a main player in this disease, the recently found capillary heterogeneity and cellular dynamics of endothelial subpopulations in BPD remain unclear. Here, we show Cap2 cells are damaged during neonatal hyperoxic injury, leading to their replacement by Cap1 cells which, in turn, significantly decline. Single-cell RNA-seq identifies the activation of numerous p53 target genes in endothelial cells, including Cdkn1a (p21) . While global deletion of p53 results in worsened vasculature, endothelial-specific deletion of p53 reverses the vascular phenotype and improves alveolar simplification during hyperoxia. This recovery is associated with the emergence of a transitional EC state, enriched for oxidative stress response genes and growth factors. These findings implicate the p53 pathway in EC type transition during injury-repair and highlights the endothelial contributions to BPD.

AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs for atlas-scale datasets like Human Pancreas Analysis Program (HPAP), we develop AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX shows the higher performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulates known islet pathobiology and shows differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

Single-Cell Analysis Identifies NOTCH3-Mediated Interactions between Stromal Cells That Promote Microenvironment Remodeling and Invasion in Lung Adenocarcinoma

Cancer immunotherapy has revolutionized the treatment of lung adenocarcinoma (LUAD); however, a significant proportion of patients do not respond. Recent transcriptomic studies to understand determinants of immunotherapy response have pinpointed stromal-mediated resistance mechanisms. To gain a better understanding of stromal biology at the cellular and molecular level in LUAD, we performed single-cell RNA sequencing of 256,379 cells, including 13,857 mesenchymal cells, from 9 treatment-naïve patients. Among the mesenchymal cell subsets, FAP+PDPN+ cancer-associated fibroblasts (CAF) and ACTA2+MCAM+ pericytes were enriched in tumors and differentiated from lung-resident fibroblasts. Imaging mass cytometry revealed that both subsets were topographically adjacent to the perivascular niche and had close spatial interactions with endothelial cells (EC). Modeling of ligand and receptor interactomes between mesenchymal and ECs identified that NOTCH signaling drives these cell-to-cell interactions in tumors, with pericytes and CAFs as the signal receivers and arterial and PLVAPhigh immature neovascular ECs as the signal senders. Either pharmacologically blocking NOTCH signaling or genetically depleting NOTCH3 levels in mesenchymal cells significantly reduced collagen production and suppressed cell invasion. Bulk RNA sequencing data demonstrated that NOTCH3 expression correlated with poor survival in stroma-rich patients and that a T cell-inflamed gene signature only predicted survival in patients with low NOTCH3. Collectively, this study provides valuable insights into the role of NOTCH3 in regulating tumor stroma biology, warranting further studies to elucidate the clinical implications of targeting NOTCH3 signaling.

SIGNIFICANCE

NOTCH3 signaling activates tumor-associated mesenchymal cells, increases collagen production, and augments cell invasion in lung adenocarcinoma, suggesting its critical role in remodeling tumor stroma.

Spatial Biology of Breast Cancer

Spatial findings have shaped on our understanding of breast cancer. In this review, we discuss how spatial methods, including spatial transcriptomics and proteomics and the resultant understanding of spatial relationships, have contributed to concepts regarding cancer progression and treatment. In addition to discussing traditional approaches, we examine how emerging multiplex imaging technologies have contributed to the field and how they might influence future research.

PI3K/mTOR inhibition induces tumour microenvironment remodelling and sensitises pS6high uterine leiomyosarcoma to PD-1 blockade

BACKGROUND

Uterine leiomyosarcomas (uLMS) are aggressive tumours with poor prognosis and limited treatment options. Although immune checkpoint blockade (ICB) has proven effective in some 'challenging-to-treat' cancers, clinical trials showed that uLMS do not respond to ICB. Emerging evidence suggests that aberrant PI3K/mTOR signalling can drive resistance to ICB. We therefore explored the relevance of the PI3K/mTOR pathway for ICB treatment in uLMS and explored pharmacological inhibition of this pathway to sensitise these tumours to ICB.

METHODS

We performed an integrated multiomics analysis based on TCGA data to explore the correlation between PI3K/mTOR dysregulation and immune infiltration in 101 LMS. We assessed response to PI3K/mTOR inhibitors in immunodeficient and humanized uLMS patient-derived xenografts (PDXs) by evaluating tumour microenvironment modulation using multiplex immunofluorescence. We explored response to single-agent and a combination of PI3K/mTOR inhibitors with PD-1 blockade in humanized uLMS PDXs. We mapped intratumoural dynamics using single-cell RNA/TCR sequencing of serially collected biopsies.

RESULTS

PI3K/mTOR over-activation (pS6high) associated with lymphocyte depletion and wound healing immune landscapes in (u)LMS, suggesting it contributes to immune evasion. In contrast, PI3K/mTOR inhibition induced profound tumour microenvironment remodelling in an ICB-resistant humanized uLMS PDX model, fostering adaptive anti-tumour immune responses. Indeed, PI3K/mTOR inhibition induced macrophage repolarisation towards an anti-tumourigenic phenotype and increased antigen presentation on dendritic and tumour cells, but also promoted infiltration of PD-1+ T cells displaying an exhausted phenotype. When combined with anti-PD-1, PI3K/mTOR inhibition led to partial or complete tumour responses, whereas no response to single-agent anti-PD-1 was observed. Combination therapy reinvigorated exhausted T cells and induced clonal hyper-expansion of a cytotoxic CD8+ T-cell population supported by a CD4+ Th1 niche.

CONCLUSIONS

Our findings indicate that aberrant PI3K/mTOR pathway activation contributes to immune escape in uLMS and provides a rationale for combining PI3K/mTOR inhibition with ICB for the treatment of this patient population.

Mast cell-derived interleukin-4 mediates activation of dendritic cell during toll-like receptor 2-mediated inflammation

Introduction

During an innate inflammation, immune cells form distinct pro- and anti-inflammatory regions around pathogen-containing core-regions. Mast cells are localized in an anti-inflammatory microenvironment during the resolution of an innate inflammation, suggesting antiinflammatory roles of these cells.

Methods

High-content imaging was used to investigated mast cell-dependent changes in the regional distribution of immune cells during an inflammation, induced by the toll-like receptor (TLR)-2 agonist zymosan.

Results

The distance between the zymosan-containing core-region and the anti-inflammatory region, described by M2-like macrophages, increased in mast cell-deficient mice. Absence of mast cells abolished dendritic cell (DC) activation, as determined by CD86-expression and localized the DCs in greater distance to zymosan particles. The CD86- DCs had a higher expression of the pro-inflammatory interleukins (IL)-1β and IL-12/23p40 as compared to activated CD86+ DCs. IL-4 administration restored CD86 expression, cytokine expression profile and localization of the DCs in mast cell-deficient mice. The IL-4 effects were mast cell-specific, since IL-4 reduction by eosinophil depletion did not affect activation of DCs.

Discussion

We found that mast cells induce DC activation selectively at the site of inflammation and thereby determine their localization within the inflammation. Overall, mast cells have antiinflammatory functions in this inflammation model and limit the size of the pro-inflammatory region surrounding the zymosan-containing core region.

In situ single-cell therapeutic response imaging facilitated by the TRIPODD fluorescence imaging platform

Purpose: Small molecule drugs such as tyrosine kinase inhibitors (TKIs) targeting tumoral molecular dependencies have become standard of care for numerous cancer types. Notably, epidermal growth factor receptor (EGFR) TKIs (e.g., erlotinib, afatinib, osimertinib) are the current first-line treatment for non-small cell lung cancer (NSCLC) due to their improved therapeutic outcomes for EGFR mutated and overexpressing disease over traditional platinum-based chemotherapy. However, many NSCLC tumors develop resistance to EGFR TKI therapy causing disease progression. Currently, the relationship between in situ drug target availability (DTA), local protein expression and therapeutic response cannot be accurately assessed using existing analytical tools despite being crucial to understanding the mechanism of therapeutic efficacy. Procedure: We have previously reported development of our fluorescence imaging platform termed TRIPODD (Therapeutic Response Imaging through Proteomic and Optical Drug Distribution) that is capable of simultaneous quantification of single-cell DTA and protein expression with preserved spatial context within a tumor. TRIPODD combines two complementary fluorescence imaging techniques: intracellular paired agent imaging (iPAI) to measure DTA and cyclic immunofluorescence (cyCIF), which utilizes oligonucleotide conjugated antibodies (Ab-oligos) for spatial proteomic expression profiling on tissue samples. Herein, TRIPODD was modified and optimized to provide a downstream analysis of therapeutic response through single-cell DTA and proteomic response imaging. Results: We successfully performed sequential imaging of iPAI and cyCIF resulting in high dimensional imaging and biomarker assessment to quantify single-cell DTA and local protein expression on erlotinib treated NSCLC models. Pharmacodynamic and pharmacokinetic studies of the erlotinib iPAI probes revealed that administration of 2.5 mg/kg each of the targeted and untargeted probe 4 h prior to tumor collection enabled calculation of DTA values with high Pearson correlation to EGFR, the erlotinib molecular target, expression in the tumors. Analysis of single-cell biomarker expression revealed that a single erlotinib dose was insufficient to enact a measurable decrease in the EGFR signaling cascade protein expression, where only the DTA metric detected the presence of bound erlotinib. Conclusion: We demonstrated the capability of TRIPODD to evaluate therapeutic response imaging to erlotinib treatment as it relates to signaling inhibition, DTA, proliferation, and apoptosis with preserved spatial context.

Bi-level Graph Learning Unveils Prognosis-Relevant Tumor Microenvironment Patterns from Breast Multiplexed Digital Pathology

The tumor microenvironment is widely recognized for its central role in driving cancer progression and influencing prognostic outcomes. Despite extensive research efforts dedicated to characterizing this complex and heterogeneous environment, considerable challenges persist. In this study, we introduce a data-driven approach for identifying patterns of cell organizations in the tumor microenvironment that are associated with patient prognoses. Our methodology relies on the construction of a bi-level graph model: (i) a cellular graph, which models the intricate tumor microenvironment, and (ii) a population graph that captures inter-patient similarities, given their respective cellular graphs, by means of a soft Weisfeiler-Lehman subtree kernel. This systematic integration of information across different scales enables us to identify patient subgroups exhibiting unique prognoses while unveiling tumor microenvironment patterns that characterize them. We demonstrate our approach in a cohort of breast cancer patients, where the identified tumor microenvironment patterns result in a risk stratification system that provides complementary, new information with respect to alternative standards. Our results, which are validated in a completely independent cohort, allow for new insights into the prognostic implications of the breast tumor microenvironment, and this methodology could be applied to other cancer types more generally.

Tissue and cellular spatiotemporal dynamics in colon aging

Tissue structure and molecular circuitry in the colon can be profoundly impacted by systemic age-related effects, but many of the underlying molecular cues remain unclear. Here, we built a cellular and spatial atlas of the colon across three anatomical regions and 11 age groups, encompassing ~1,500 mouse gut tissues profiled by spatial transcriptomics and ~400,000 single nucleus RNA-seq profiles. We developed a new computational framework, cSplotch, which learns a hierarchical Bayesian model of spatially resolved cellular expression associated with age, tissue region, and sex, by leveraging histological features to share information across tissue samples and data modalities. Using this model, we identified cellular and molecular gradients along the adult colonic tract and across the main crypt axis, and multicellular programs associated with aging in the large intestine. Our multi-modal framework for the investigation of cell and tissue organization can aid in the understanding of cellular roles in tissue-level pathology.

An autophagy program that promotes T cell egress from the lymph node controls responses to immune checkpoint blockade

Lymphatic endothelial cells (LECs) of the lymph node (LN) parenchyma orchestrate leukocyte trafficking and peripheral T cell dynamics. T cell responses to immunotherapy largely rely on peripheral T cell recruitment in tumors. Yet, a systematic and molecular understanding of how LECs within the LNs control T cell dynamics under steady-state and tumor-bearing conditions is lacking. Intravital imaging combined with immune phenotyping shows that LEC-specific deletion of the essential autophagy gene Atg5 alters intranodal positioning of lymphocytes and accrues their persistence in the LNs by increasing the availability of the main egress signal sphingosine-1-phosphate. Single-cell RNA sequencing of tumor-draining LNs shows that loss of ATG5 remodels niche-specific LEC phenotypes involved in molecular pathways regulating lymphocyte trafficking and LEC-T cell interactions. Functionally, loss of LEC autophagy prevents recruitment of tumor-infiltrating T and natural killer cells and abrogates response to immunotherapy. Thus, an LEC-autophagy program boosts immune-checkpoint responses by guiding systemic T cell dynamics.

B cells in the pneumococcus-infected lung are heterogeneous and require CD4+ T cell help including CD40L to become resident memory B cells

Recovery from respiratory pneumococcal infections generates lung-localized protection against heterotypic bacteria, mediated by resident memory lymphocytes. Optimal protection in mice requires re-exposure to pneumococcus within days of initial infection. Serial surface marker phenotyping of B cell populations in a model of pneumococcal heterotypic immunity revealed that bacterial re-exposure stimulates the immediate accumulation of dynamic and heterogeneous populations of B cells in the lung, and is essential for the establishment of lung resident memory B (BRM) cells. The B cells in the early wave were activated, proliferating locally, and associated with both CD4+ T cells and CXCL13. Antagonist- and antibody-mediated interventions were implemented during this early timeframe to demonstrate that lymphocyte recirculation, CD4+ cells, and CD40 ligand (CD40L) signaling were all needed for lung BRM cell establishment, whereas CXCL13 signaling was not. While most prominent as aggregates in the loose connective tissue of bronchovascular bundles, morphometry and live lung imaging analyses showed that lung BRM cells were equally numerous as single cells dispersed throughout the alveolar septae. We propose that CD40L signaling from antigen-stimulated CD4+ T cells in the infected lung is critical to establishment of local BRM cells, which subsequently protect the airways and parenchyma against future potential infections.

Multiplexed single-cell lineage tracing of mitotic kinesin inhibitor resistance in glioblastoma

Glioblastoma (GBM) is a deadly brain tumor, and the kinesin motor KIF11 is an attractive therapeutic target with roles in proliferation and invasion. Resistance to KIF11 inhibitors, which has mainly been studied in animal models, presents significant challenges. We use lineage-tracing barcodes and single-cell RNA sequencing to analyze resistance in patient-derived GBM neurospheres treated with ispinesib, a potent KIF11 inhibitor. Similar to GBM progression in patients, untreated cells lose their neural lineage identity and become mesenchymal, which is associated with poor prognosis. Conversely, cells subjected to long-term ispinesib treatment exhibit a proneural phenotype. We generate patient-derived xenografts and show that ispinesib-resistant cells form less aggressive tumors in vivo, even in the absence of drug. Moreover, treatment of human ex vivo GBM slices with ispinesib demonstrates phenotypic alignment with in vitro responses, underscoring the clinical relevance of our findings. Finally, using retrospective lineage tracing, we identify drugs that are synergistic with ispinesib.

Spatiotemporal architecture of immune cells and cancer-associated fibroblasts in high-grade serous ovarian carcinoma

High-grade serous ovarian carcinoma (HGSOC), the deadliest form of ovarian cancer, is typically diagnosed after it has metastasized and often relapses after standard-of-care platinum-based chemotherapy, likely due to advanced tumor stage, heterogeneity, and immune evasion and tumor-promoting signaling from the tumor microenvironment. To understand how spatial heterogeneity contributes to HGSOC progression and early relapse, we profiled an HGSOC tissue microarray of patient-matched longitudinal samples from 42 patients. We found spatial patterns associated with early relapse, including changes in T cell localization, malformed tertiary lymphoid structure (TLS)-like aggregates, and increased podoplanin-positive cancer-associated fibroblasts (CAFs). Using spatial features to compartmentalize the tissue, we found that plasma cells distribute in two different compartments associated with TLS-like aggregates and CAFs, and these distinct microenvironments may account for the conflicting reports about the role of plasma cells in HGSOC prognosis.

IL-1 receptor 1 signaling shapes the development of viral antigen-specific CD4+ T cell responses following COVID-19 mRNA vaccination

BACKGROUND

The innate immune cytokine interleukin (IL)-1 can affect T cell immunity, a critical factor in host defense. In a previous study, we identified a subset of human CD4+ T cells which express IL-1 receptor 1 (IL-1R1). However, the expression of such receptor by viral antigen-specific CD4+ T cells and its biological implication remain largely unexplored. This led us to investigate the implication of IL-1R1 in the development of viral antigen-specific CD4+ T cell responses in humans, including healthy individuals and patients with primary antibody deficiency (PAD), and animals.

METHODS

We characterized CD4+ T cells specific for SARS-CoV-2 spike (S) protein, influenza virus, and cytomegalovirus utilizing multiplexed single cell RNA-seq, mass cytometry and flow cytometry followed by an animal study.

FINDINGS

In healthy individuals, CD4+ T cells specific for viral antigens, including S protein, highly expressed IL-1R1. IL-1β promoted interferon (IFN)-γ expression by S protein-stimulated CD4+ T cells, supporting the functional implication of IL-1R1. Following the 2nd dose of COVID-19 mRNA vaccines, S protein-specific CD4+ T cells with high levels of IL-1R1 increased, likely reflecting repetitive antigenic stimulation. The expression levels of IL-1R1 by such cells correlated with the development of serum anti-S protein IgG antibody. A similar finding of increased expression of IL-1R1 by S protein-specific CD4+ T cells was also observed in patients with PAD following COVID-19 mRNA vaccination although the expression levels of IL-1R1 by such cells did not correlate with the levels of serum anti-S protein IgG antibody. In mice immunized with COVID-19 mRNA vaccine, neutralizing IL-1R1 decreased IFN-γ expression by S protein-specific CD4+ T cells and the development of anti-S protein IgG antibody.

INTERPRETATION

Our results demonstrate the significance of IL-1R1 expression in CD4+ T cells for the development of viral antigen-specific CD4+ T cell responses, contributing to humoral immunity. This provides an insight into the regulation of adaptive immune responses to viruses via the IL-1 and IL-1R1 interface.

FUNDING

Moderna to HJP, National Institutes of Health (NIH) 1R01AG056728 and R01AG055362 to IK and KL2 TR001862 to JJS, Quest Diagnostics to IK and RB, and the Mathers Foundation to RB.

An esophagus cell atlas reveals dynamic rewiring during active eosinophilic esophagitis and remission

Coordinated cell interactions within the esophagus maintain homeostasis, and disruption can lead to eosinophilic esophagitis (EoE), a chronic inflammatory disease with poorly understood pathogenesis. We profile 421,312 individual cells from the esophageal mucosa of 7 healthy and 15 EoE participants, revealing 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unclear features of EoE. Active disease displays enrichment of ALOX15+ macrophages, PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Ligand-receptor expression uncovers eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+ TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures includes mast and CD4+ TRM cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors. These cellular alterations in EoE and remission advance our understanding of eosinophilic inflammation and opportunities for therapeutic intervention.

Reconstruction of macroglia and adult neurogenesis evolution through cross-species single-cell transcriptomic analyses

Macroglia fulfill essential functions in the adult vertebrate brain, producing and maintaining neurons and regulating neuronal communication. However, we still know little about their emergence and diversification. We used the zebrafish D. rerio as a distant vertebrate model with moderate glial diversity as anchor to reanalyze datasets covering over 600 million years of evolution. We identify core features of adult neurogenesis and innovations in the mammalian lineage with a potential link to the rarity of radial glia-like cells in adult humans. Our results also suggest that functions associated with astrocytes originated in a multifunctional cell type fulfilling both neural stem cell and astrocytic functions before these diverged. Finally, we identify conserved elements of macroglial cell identity and function and their time of emergence during evolution.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and γH2AX foci links to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-μm micropores. To search for pancreatic cancer cell of origin, analysis of single-cell data sets revealed that the extracellular matrix shapes an alternate route of acinar-ductal transdifferentiation of acinar cells into a central hub of elegantly restrained topoisomerase II α (TOP2A)-overexpressing cancer cells that spread out as unique cancer clusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in patients with The Cancer Genome Atlas-PAAD. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased VAV1 expression, and prolonged overall survival in the KPC mice. Decreases of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodel the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induces the onset of cells of origin of cancer by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment against cells-of-origin cancer.

The neuroendocrine transition in prostate cancer is dynamic and dependent on ASCL1

Lineage plasticity is a recognized hallmark of cancer progression that can shape therapy outcomes. The underlying cellular and molecular mechanisms mediating lineage plasticity remain poorly understood. Here, we describe a versatile in vivo platform to identify and interrogate the molecular determinants of neuroendocrine lineage transformation at different stages of prostate cancer progression. Adenocarcinomas reliably develop following orthotopic transplantation of primary mouse prostate organoids acutely engineered with human-relevant driver alterations (e.g., Rb1-/-; Trp53-/-; cMyc+ or Pten-/-; Trp53-/-; cMyc+), but only those with Rb1 deletion progress to ASCL1+ neuroendocrine prostate cancer (NEPC), a highly aggressive, androgen receptor signaling inhibitor (ARSI)-resistant tumor. Importantly, we show this lineage transition requires a native in vivo microenvironment not replicated by conventional organoid culture. By integrating multiplexed immunofluorescence, spatial transcriptomics and PrismSpot to identify cell type-specific spatial gene modules, we reveal that ASCL1+ cells arise from KRT8+ luminal epithelial cells that progressively acquire transcriptional heterogeneity, producing large ASCL1+;KRT8- NEPC clusters. Ascl1 loss in established NEPC results in transient tumor regression followed by recurrence; however, Ascl1 deletion prior to transplantation completely abrogates lineage plasticity, yielding adenocarcinomas with elevated AR expression and marked sensitivity to castration. The dynamic feature of this model reveals the importance of timing of therapies focused on lineage plasticity and offers a platform for identification of additional lineage plasticity drivers.

CD8 effector T cells enhance teclistamab response in BCMA-exposed and -naïve multiple myeloma

ABSTRACT

Teclistamab, a B-cell maturation antigen (BCMA)- and CD3-targeting bispecific antibody, is an effective novel treatment for relapsed/refractory multiple myeloma (R/RMM), but efficacy in patients exposed to BCMA-directed therapies and mechanisms of resistance have yet to be fully delineated. We conducted a real-world retrospective study of commercial teclistamab, capturing both clinical outcomes and immune correlates of treatment response in a cohort of patients (n = 52) with advanced R/RMM. Teclistamab was highly effective with an overall response rate (ORR) of 64%, including an ORR of 50% for patients with prior anti-BCMA therapy. Pretreatment plasma cell BCMA expression levels had no bearing on response. However, comprehensive pretreatment immune profiling identified that effector CD8+ T-cell populations were associated with response to therapy and a regulatory T-cell population associated with nonresponse, indicating a contribution of immune status in outcomes with potential utility as a biomarker signature to guide patient management.

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 ).

SifiNet: A robust and accurate method to identify feature gene sets and annotate cells

SifiNet is a robust and accurate computational pipeline for identifying distinct gene sets, extracting and annotating cellular subpopulations, and elucidating intrinsic relationships among these subpopulations. Uniquely, SifiNet bypasses the cell clustering stage, commonly integrated into other cellular annotation pipelines, thereby circumventing potential inaccuracies in clustering that may compromise subsequent analyses. Consequently, SifiNet has demonstrated superior performance in multiple experimental datasets compared with other state-of-the-art methods. SifiNet can analyze both single-cell RNA and ATAC sequencing data, thereby rendering comprehensive multiomic cellular profiles. It is conveniently available as an open-source R package.

Faecal microbial transfer and complex carbohydrates mediate protection against COPD

OBJECTIVE

Chronic obstructive pulmonary disease (COPD) is a major cause of global illness and death, most commonly caused by cigarette smoke. The mechanisms of pathogenesis remain poorly understood, limiting the development of effective therapies. The gastrointestinal microbiome has been implicated in chronic lung diseases via the gut-lung axis, but its role is unclear.

DESIGN

Using an in vivo mouse model of cigarette smoke (CS)-induced COPD and faecal microbial transfer (FMT), we characterised the faecal microbiota using metagenomics, proteomics and metabolomics. Findings were correlated with airway and systemic inflammation, lung and gut histopathology and lung function. Complex carbohydrates were assessed in mice using a high resistant starch diet, and in 16 patients with COPD using a randomised, double-blind, placebo-controlled pilot study of inulin supplementation.

RESULTS

FMT alleviated hallmark features of COPD (inflammation, alveolar destruction, impaired lung function), gastrointestinal pathology and systemic immune changes. Protective effects were additive to smoking cessation, and transfer of CS-associated microbiota after antibiotic-induced microbiome depletion was sufficient to increase lung inflammation while suppressing colonic immunity in the absence of CS exposure. Disease features correlated with the relative abundance of Muribaculaceae, Desulfovibrionaceae and Lachnospiraceae family members. Proteomics and metabolomics identified downregulation of glucose and starch metabolism in CS-associated microbiota, and supplementation of mice or human patients with complex carbohydrates improved disease outcomes.

CONCLUSION

The gut microbiome contributes to COPD pathogenesis and can be targeted therapeutically.

Immune Modulation with RANKL Blockade through Denosumab Treatment in Patients with Cancer

Denosumab is a fully human mAb that binds receptor activator of NFκB ligand (RANKL). It is routinely administered to patients with cancer to reduce the incidence of new bone metastasis. RANK-RANKL interactions regulate bone turnover by controlling osteoclast recruitment, development, and activity. However, these interactions also can regulate immune cells including dendritic cells and medullary thymic epithelial cells. Inhibition of the latter results in reduced thymic negative selection of T cells and could enhance the generation of tumor-specific T cells. We examined whether administering denosumab could modify modulate circulating immune cells in patients with cancer. Blood was collected from 23 patients with prostate cancer and 3 patients with renal cell carcinoma, all of whom had advanced disease and were receiving denosumab, prior to and during denosumab treatment. Using high-dimensional mass cytometry, we found that denosumab treatment by itself induced modest effects on circulating immune cell frequency and activation. We also found minimal changes in the circulating T-cell repertoire and the frequency of new thymic emigrants with denosumab treatment. However, when we stratified patients by whether they were receiving chemotherapy and/or steroids, patients receiving these concomitant treatments showed significantly greater immune modulation, including an increase in the frequency of natural killer cells early and classical monocytes later. We also saw broad induction of CTLA-4 and TIM3 expression in circulating lymphocytes and some monocyte populations. These findings suggest that denosumab treatment by itself has modest immunomodulatory effects, but when combined with conventional cancer treatments, can lead to the induction of immunologic checkpoints. See related Spotlight by Nasrollahi and Davar, p. 383.

The covariance environment defines cellular niches for spatial inference

A key challenge of analyzing data from high-resolution spatial profiling technologies is to suitably represent the features of cellular neighborhoods or niches. Here we introduce the covariance environment (COVET), a representation that leverages the gene-gene covariate structure across cells in the niche to capture the multivariate nature of cellular interactions within it. We define a principled optimal transport-based distance metric between COVET niches that scales to millions of cells. Using COVET to encode spatial context, we developed environmental variational inference (ENVI), a conditional variational autoencoder that jointly embeds spatial and single-cell RNA sequencing data into a latent space. ENVI includes two decoders: one to impute gene expression across the spatial modality and a second to project spatial information onto single-cell data. ENVI can confer spatial context to genomics data from single dissociated cells and outperforms alternatives for imputing gene expression on diverse spatial datasets.

Meloxicam treatment disrupts the regional structure of innate inflammation sites by targeting the pro-inflammatory effects of prostanoids

BACKGROUND AND PURPOSE

Non-steroidal anti-inflammatory drugs (NSAIDs) are the most widely prescribed drugs in the world due to their analgesic, antipyretic and anti-inflammatory effects. However, NSAIDs inhibit prostanoid synthesis, interfering with their pro-inflammatory and anti-inflammatory functions and potentially prolonging acute inflammation.

EXPERIMENTAL APPROACH

We used high-content immunohistochemistry to define the impact of meloxicam treatment on spatially separated pro-inflammatory and anti-inflammatory processes during innate inflammation in mice induced by zymosan. This allowed us to determine the effect of meloxicam treatment on the organization of pro-inflammatory and anti-inflammatory microenvironments, thereby identifying relevant changes in immune cell localization, recruitment and activation.

KEY RESULTS

Meloxicam treatment reduced zymosan-induced thermal hypersensitivity at early time points but delayed its resolution. High-content immunohistochemistry revealed that the pro-inflammatory area was smaller after treatment, diminishing neutrophil recruitment, M1-like macrophage polarization, and especially phagocytosis by neutrophils and macrophages. The polarization of macrophages towards the M2-like anti-inflammatory phenotype was unaffected, and the number of anti-inflammatory eosinophils actually increased.

CONCLUSION AND IMPLICATIONS

High-content immunohistochemistry was able to identify relevant meloxicam-mediated effects on inflammatory processes based on alterations in the regional structure of inflammation sites. Meloxicam delayed the clearance of pathogens by inhibiting pro-inflammatory processes, causing prolonged inflammation. Our data suggest that the prescription of NSAIDs as a treatment during an acute pathogen-driven inflammation should be reconsidered in patients with compromised immune systems.

Depleting myeloid-biased haematopoietic stem cells rejuvenates aged immunity

Ageing of the immune system is characterized by decreased lymphopoiesis and adaptive immunity, and increased inflammation and myeloid pathologies1,2. Age-related changes in populations of self-renewing haematopoietic stem cells (HSCs) are thought to underlie these phenomena3. During youth, HSCs with balanced output of lymphoid and myeloid cells (bal-HSCs) predominate over HSCs with myeloid-biased output (my-HSCs), thereby promoting the lymphopoiesis required for initiating adaptive immune responses, while limiting the production of myeloid cells, which can be pro-inflammatory4. Ageing is associated with increased proportions of my-HSCs, resulting in decreased lymphopoiesis and increased myelopoiesis3,5,6. Transfer of bal-HSCs results in abundant lymphoid and myeloid cells, a stable phenotype that is retained after secondary transfer; my-HSCs also retain their patterns of production after secondary transfer5. The origin and potential interconversion of these two subsets is still unclear. If they are separate subsets postnatally, it might be possible to reverse the ageing phenotype by eliminating my-HSCs in aged mice. Here we demonstrate that antibody-mediated depletion of my-HSCs in aged mice restores characteristic features of a more youthful immune system, including increasing common lymphocyte progenitors, naive T cells and B cells, while decreasing age-related markers of immune decline. Depletion of my-HSCs in aged mice improves primary and secondary adaptive immune responses to viral infection. These findings may have relevance to the understanding and intervention of diseases exacerbated or caused by dominance of the haematopoietic system by my-HSCs.

Ultrasound-responsive spherical nucleic acid against c-Myc/PD-L1 to enhance anti-tumoral macrophages in triple-negative breast cancer progression

Triple-negative breast cancer (TNBC) is the most challenging breast cancer subtype because of its aggressive behavior and limited therapeutic targets. c-Myc is hyperactivated in the majority of TNBC tissues, however, it has been considered an "undruggable" target due to its disordered structure. Herein, we developed an ultrasound-responsive spherical nucleic acid (SNA) against c-Myc and PD-L1 in TNBC. It is a self-assembled and carrier-free system composed of a hydrophilic small-interfering RNA (siRNA) shell and a hydrophobic core made of a peptide nucleic acid (PNA)-based antisense oligonucleotide (ASO) and a sonosensitizer. We accomplished significant enrichment in the tumor by enhanced permeability and retention (EPR) effect, the controllable release of effective elements by ultrasound activation, and the combination of targeted therapy, immunotherapy and physiotherapy. Our study demonstrated significant anti-tumoral effects in vitro and in vivo. Mass cytometry showed an invigorated tumor microenvironment (TME) characterized by a significant alteration in the composition of tumor-associated macrophages (TAM) and decreased proportion of PD-1-positive (PD-1+) T effector cells after appropriate treatment of the ultrasound-responsive SNA (USNA). Further experiments verified that tumor-conditioned macrophages residing in the TME were transformed into the anti-tumoral population. Our finding offers a novel therapeutic strategy against the "undruggable" c-Myc, develops a new targeted therapy for c-Myc/PD-L1 and provides a treatment option for the TNBC.

In vivo single-cell high-dimensional mass cytometry analysis to track the interactions between Klebsiella pneumoniae and myeloid cells

In vivo single-cell approaches have transformed our understanding of the immune populations in tissues. Mass cytometry (CyTOF), that combines the resolution of mass spectrometry with the ability to conduct multiplexed measurements of cell molecules at the single cell resolution, has enabled to resolve the diversity of immune cell subsets, and their heterogeneous functionality. Here we assess the feasibility of taking CyTOF one step further to immuno profile cells while tracking their interactions with bacteria, a method we term Bac-CyTOF. We focus on the pathogen Klebsiella pneumoniae interrogating the pneumonia mouse model. Using Bac-CyTOF, we unveil the atlas of immune cells of mice infected with a K. pneumoniae hypervirulent strain. The atlas is characterized by a decrease in the populations of alveolar and monocyte-derived macrophages. Conversely, neutrophils, and inflammatory monocytes are characterized by an increase in the subpopulations expressing markers of less active cells such as the immune checkpoint PD-L1. These are the cells infected. We show that the type VI secretion system (T6SS) contributes to shape the lung immune landscape. The T6SS governs the interaction with monocytes/macrophages by shifting Klebsiella from alveolar macrophages to interstitial macrophages and limiting the infection of inflammatory monocytes. The lack of T6SS results in an increase of cells expressing markers of active cells, and a decrease in the subpopulations expressing PD-L1. By probing Klebsiella, and Acinetobacter baumannii strains with limited ability to survive in vivo, we uncover that a heightened recruitment of neutrophils, and relative high levels of alveolar macrophages and eosinophils and the recruitment of a characteristic subpopulation of neutrophils are features of mice clearing infections. We leverage Bac-CyTOF-generated knowledge platform to investigate the role of the DNA sensor STING in Klebsiella infections. sting-/- infected mice present features consistent with clearing the infection including the reduced levels of PD-L1. STING absence facilitates Klebsiella clearance.

Single-cell analysis of 5-aminolevulinic acid intraoperative labeling specificity for glioblastoma

OBJECTIVE

Glioblastoma (GBM) is the most common and aggressive malignant primary brain tumor, and resection is a key part of the standard of care. In fluorescence-guided surgery (FGS), fluorophores differentiate tumor tissue from surrounding normal brain. The heme synthesis pathway converts 5-aminolevulinic acid (5-ALA), a fluorogenic substrate used for FGS, to fluorescent protoporphyrin IX (PpIX). The resulting fluorescence is believed to be specific to neoplastic glioma cells, but this specificity has not been examined at a single-cell level. The objective of this study was to determine the specificity with which 5-ALA labels the diversity of cell types in GBM.

METHODS

The authors performed single-cell optical phenotyping and expression sequencing-version 2 (SCOPE-seq2), a paired single-cell imaging and RNA sequencing method, of individual cells on human GBM surgical specimens with macroscopically visible PpIX fluorescence from patients who received 5-ALA prior to surgery. SCOPE-seq2 allowed the authors to simultaneously image PpIX fluorescence and unambiguously identify neoplastic cells from single-cell RNA sequencing. Experiments were also conducted in cell culture and co-culture models of glioma and in acute slice cultures from a mouse glioma model to investigate cell- and tissue-specific uptake and secretion of 5-ALA and PpIX.

RESULTS

SCOPE-seq2 analysis of human GBM surgical specimens revealed that 5-ALA treatment resulted in labeling that was not specific to neoplastic glioma cells. The cell culture further demonstrated that nonneoplastic cells could be labeled by 5-ALA directly or by PpIX secreted from surrounding neoplastic cells. Acute slice cultures from mouse glioma models showed that 5-ALA preferentially labeled GBM tumor tissue over nonneoplastic brain tissue with significant labeling in the tumor margins, and that this contrast was not due to blood-brain barrier disruption.

CONCLUSIONS

Together, these findings support the use of 5-ALA as an indicator of GBM tissue but question the main advantage of 5-ALA for specific intracellular labeling of neoplastic glioma cells in FGS. Further studies are needed to systematically compare the performance of 5-ALA to that of potential alternatives for FGS.

DELVE: feature selection for preserving biological trajectories in single-cell data

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .

SpiDe-Sr: blind super-resolution network for precise cell segmentation and clustering in spatial proteomics imaging

Spatial proteomics elucidates cellular biochemical changes with unprecedented topological level. Imaging mass cytometry (IMC) is a high-dimensional single-cell resolution platform for targeted spatial proteomics. However, the precision of subsequent clinical analysis is constrained by imaging noise and resolution. Here, we propose SpiDe-Sr, a super-resolution network embedded with a denoising module for IMC spatial resolution enhancement. SpiDe-Sr effectively resists noise and improves resolution by 4 times. We demonstrate SpiDe-Sr respectively with cells, mouse and human tissues, resulting 18.95%/27.27%/21.16% increase in peak signal-to-noise ratio and 15.95%/31.63%/15.52% increase in cell extraction accuracy. We further apply SpiDe-Sr to study the tumor microenvironment of a 20-patient clinical breast cancer cohort with 269,556 single cells, and discover the invasion of Gram-negative bacteria is positively correlated with carcinogenesis markers and negatively correlated with immunological markers. Additionally, SpiDe-Sr is also compatible with fluorescence microscopy imaging, suggesting SpiDe-Sr an alternative tool for microscopy image super-resolution.

Differences in the expression of cortex-wide neural dynamics are related to behavioral phenotype

Behavior differs across individuals, ranging from typical to atypical phenotypes.1 Understanding how differences in behavior relate to differences in neural activity is critical for developing treatments of neuropsychiatric and neurodevelopmental disorders. One hypothesis is that differences in behavior reflect individual differences in the dynamics of how information flows through the brain. In support of this, the correlation of neural activity between brain areas, termed "functional connectivity," varies across individuals2 and is disrupted in autism,3 schizophrenia,4 and depression.5 However, the changes in neural activity that underlie altered behavior and functional connectivity remain unclear. Here, we show that individual differences in the expression of different patterns of cortical neural dynamics explain variability in both functional connectivity and behavior. Using mesoscale imaging, we recorded neural activity across the dorsal cortex of behaviorally "typical" and "atypical" mice. All mice shared the same recurring cortex-wide spatiotemporal motifs of neural activity, and these motifs explained the large majority of variance in cortical activity (>75%). However, individuals differed in how frequently different motifs were expressed. These differences in motif expression explained differences in functional connectivity and behavior across both typical and atypical mice. Our results suggest that differences in behavior and functional connectivity are due to changes in the processes that select which pattern of neural activity is expressed at each moment in time.

Spatial Mass Cytometry-Based Single-Cell Imaging Reveals a Disrupted Epithelial-Immune Axis in Prurigo Nodularis

Prurigo nodularis (PN) is a chronic, inflammatory skin condition that disproportionately affects African Americans and features intensely pruritic, hyperkeratotic nodules on the extremities and trunk. PN is understudied compared with other inflammatory skin diseases, with the spatial organization of the cutaneous infiltrate in PN yet to be characterized. In this work, we employ spatial imaging mass cytometry to visualize PN lesional skin inflammation and architecture with single-cell resolution through an unbiased machine learning approach. PN lesional skin has increased expression of caspase 3, NF-kB, and phosphorylated signal transducer and activator of transcription 3 compared with healthy skin. Keratinocytes in lesional skin are subdivided into CD14+CD33+, CD11c+, CD63+, and caspase 3-positive innate subpopulations. CD14+ macrophage populations expressing phosphorylated extracellular signal-regulated kinase 1/2 correlate positively with patient-reported itch (P = .006). Hierarchical clustering reveals a cluster of patients with PN with greater atopy, increased NF-kB+ signal transducer and activator of transcription 3-positive phosphorylated extracellular signal-regulated kinase 1/2-positive monocyte-derived myeloid dendritic cells, and increased vimentin expression (P < .05). Neighborhood analysis finds interactions between CD14+ macrophages, CD3+ T cells, monocyte-derived myeloid dendritic cells, and keratinocytes expressing innate immune markers. These findings highlight phosphorylated extracellular signal-regulated kinase-positive CD14+ macrophages as contributors to itch and suggest an epithelial-immune axis in PN pathogenesis.

Exploring NK cell receptor dynamics in paediatric leukaemias: implications for immunotherapy and prognosis

Objectives

Immunotherapies targeting natural killer (NK) cell receptors have shown promise against leukaemia. Unfortunately, cancer immunosuppressive mechanisms that alter NK cell phenotype prevent such approaches from being successful. The study utilises advanced cytometry to examine how cancer immunosuppressive pathways affect NK cell phenotypic changes in clinical samples.

Methods

In this study, we conducted a high-dimensional examination of the cell surface expression of 16 NK cell receptors in paediatric patients with acute myeloid leukaemia and acute lymphoblastic leukaemia, as well as in samples of non-age matched adult peripheral blood (APB) and umbilical cord blood (UCB). An unsupervised analysis was carried out in order to identify NK cell populations present in paediatric leukaemias.

Results

We observed that leukaemia NK cells clustered together with UCB NK cells and expressed relatively higher levels of the NKG2A receptor compared to APB NK cells. In addition, CD56dimCD16+CD57- NK cells lacking NKG2A expression were mainly absent in paediatric leukaemia patients. However, CD56br NK cell populations expressing high levels of NKG2A were highly represented in paediatric leukaemia patients. NKG2A expression on leukaemia NK cells was found to be positively correlated with the expression of its ligand, suggesting that the NKG2A-HLA-E interaction may play a role in modifying NK cell responses to leukaemia cells.

Conclusion

We provide an in-depth analysis of NK cell populations in paediatric leukaemia patients. These results support the development of immunotherapies targeting immunosuppressive receptors, such as NKG2A, to enhance innate immunity against paediatric leukaemia.

Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Microbial cancer immunotherapy reprograms hematopoietic stem cells to enhance anti-tumor immunity

Mycobacterium bovis BCG is the vaccine against tuberculosis and an immunotherapy for bladder cancer. When administered intravenously, BCG reprograms bone marrow hematopoietic stem and progenitor cells (HSPCs), leading to heterologous protection against infections. Whether HSPC-reprogramming contributes to the anti-tumor effects of BCG administered into the bladder is unknown. We demonstrate that BCG administered in the bladder in both mice and humans reprograms HSPCs to amplify myelopoiesis and functionally enhance myeloid cell antigen presentation pathways. Reconstitution of naive mice with HSPCs from bladder BCG-treated mice enhances anti-tumor immunity and tumor control, increases intratumoral dendritic cell infiltration, and synergizes with checkpoint blockade. We conclude that bladder BCG acts systemically, reprogramming HSPC-encoded innate immunity, highlighting the broad potential of modulating HSPC phenotypes to improve tumor immunity.

One Sentence Summary

BCG administered in the bladder reprograms bone marrow HSPCs and contributes to tumor control via enhanced myeloid cells.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

Circulation immune cell landscape in canonical pathogenesis of colorectal adenocarcinoma by CyTOF analysis

Current studies on the immune microenvironment of colorectal cancer (CRC) were mostly limited to the tissue level, lacking relevant studies in the peripheral blood, and failed to describe its alterations in the whole process of adenocarcinoma formation, especially of adenoma carcinogenesis. Here, we constructed a large-scale population cohort and used the CyTOF to explore the changes of various immune cell subsets in peripheral blood of CRC. We found monocytes and basophils cells were significantly higher in adenocarcinoma patients. Compared with early-stage CRC, effector CD4+T cells and naive B cells were higher in patients with lymph node metastasis, whereas the basophils were lower. We also performed random forest algorithm and found monocytes play the key role in carcinogenesis. Our study draws a peripheral blood immune cell landscape of the occurrence and development of CRC at the single-cell level and provides a reference for other researchers.

Single cell atlas of human gastric muscle immune cells and macrophage-driven changes in idiopathic gastroparesis

Gastrointestinal immune cells, particularly muscularis macrophages (MM) interact with the enteric nervous system and influence gastrointestinal motility. Here we determine the human gastric muscle immunome and its changes in patients with idiopathic gastroparesis (IG). Single cell sequencing was performed on 26,000 CD45+ cells obtained from the gastric tissue of 20 subjects. We demonstrate 11 immune cell clusters with T cells being most abundant followed by myeloid cells. The proportions of cells belonging to the 11 clusters were similar between IG and controls. However, 9/11 clusters showed 578-11,429 differentially expressed genes. In IG, MM had decreased expression of tissue-protective and microglial genes and increased the expression of monocyte trafficking and stromal activating genes. Furthermore, in IG, IL12 mediated JAK-STAT signaling involved in the activation of tissue-resident macrophages and Eph-ephrin signaling involved in monocyte chemotaxis were upregulated. Patients with IG had a greater abundance of monocyte-like cells. These data further link immune dysregulation to the pathophysiology of gastroparesis.

Differentiation marker-negative CD4+ T cells persist after yellow fever virus vaccination and contribute to durable memory

Factors that contribute to durable immunological memory remain incompletely understood. In our longitudinal analyses of CD4+ T cell responses to the yellow fever virus (YFV) vaccine by peptide-MHC tetramers, we unexpectedly found naïve phenotype virus-specific CD4+ T cells that persisted months to years after immunization. These Marker negative T cells (TMN) lacked CD95, CXCR3, CD11a, and CD49d surface protein expression, distinguishing them from previously discovered stem-cell memory T cells. Functionally, they resembled genuine naïve T cells upon in vitro stimulation. Single-cell TCR sequencing detected expanded clonotypes within the TMN subset and identified a shared repertoire with memory and effector T cells. T cells expressing TMN-associated TCRs were rare before vaccination, suggesting their expansion following vaccination. Longitudinal tracking of YFV-specific responses over the subsequent years revealed superior stability of the TMN subset and their association with the longevity of the overall population. The identification of these long-lived, antigen-experienced T cells may inform the design of durable T cell-based vaccines and engineered T cell therapies.

Application of Mass Cytometry Platforms to Solid Organ Transplantation

Transplantation serves as the cornerstone of treatment for patients with end-stage organ disease. The prevalence of complications, such as allograft rejection, infection, and malignancies, underscores the need to dissect the complex interactions of the immune system at the single-cell level. In this review, we discuss studies using mass cytometry or cytometry by time-of-flight, a cutting-edge technology enabling the characterization of immune populations and cell-to-cell interactions in granular detail. We review the application of mass cytometry in human and experimental animal studies in the context of transplantation, uncovering invaluable contributions of the tool to understanding rejection and other transplant-related complications. We discuss recent innovations that have the potential to streamline and standardize mass cytometry workflows for application to multisite clinical trials. Additionally, we introduce imaging mass cytometry, a technique that couples the power of mass cytometry with spatial context, thereby mapping cellular interactions within tissue microenvironments. The synergistic integration of mass cytometry and imaging mass cytometry data with other omics data sets and high-dimensional data platforms to further define immune dynamics is discussed. In conclusion, mass cytometry technologies, when integrated with other tools and data, shed light on the intricate landscape of the immune response in transplantation. This approach holds significant potential for enhancing patient outcomes by advancing our understanding and facilitating the development of new diagnostics and therapeutics.