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

Dissecting the immune infiltrate of primary luminal B-like breast carcinomas in relation to age

The impact of aging on the immune landscape of luminal breast cancer (Lum-BC) is poorly characterized. Understanding the age-related dynamics of immune editing in Lum-BC is anticipated to improve the therapeutic benefit of immunotherapy in older patients. To this end, here we applied the 'multiple iterative labeling by antibody neo-deposition' (MILAN) technique, a spatially resolved single-cell multiplex immunohistochemistry method. We created tissue microarrays by sampling both the tumor center and invasive front of luminal breast tumors collected from a cohort of treatment-naïve patients enrolled in the prospective monocentric IMAGE (IMmune system and AGEing) study. Patients were subdivided into three nonoverlapping age categories (35-45 = 'young', n = 12; 55-65 = 'middle', n = 15; ≥70 = 'old', n = 26). Additionally, depending on localization and amount of cytotoxic T lymphocytes, the tumor immune types 'desert' (n = 22), 'excluded' (n = 19), and 'inflamed' (n = 12) were identified. For the MILAN technique we used 58 markers comprising phenotypic and functional markers allowing in-depth characterization of T and B lymphocytes (T&B-lym). These were compared between age groups and tumor immune types using Wilcoxon's test and Pearson's correlation. Cytometric analysis revealed a decline of the immune cell compartment with aging. T&B-lym were numerically less abundant in tumors from middle-aged and old compared to young patients, regardless of the geographical tumor zone. Likewise, desert-type tumors showed the smallest immune-cell compartment and were not represented in the group of young patients. Analysis of immune checkpoint molecules revealed a heterogeneous geographical pattern of expression, indicating higher numbers of PD-L1 and OX40-positive T&B-lym in young compared to old patients. Despite the numerical decline of immune infiltration, old patients retained higher expression levels of OX40 in T helper cells located near cancer cells, compared to middle-aged and young patients. Aging is associated with important numerical and functional changes of the immune landscape in Lum-BC. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Identification of single-cell blasts in pediatric acute myeloid leukemia using an autoencoder

Pediatric acute myeloid leukemia (AML) is an aggressive blood cancer with a poor prognosis and high relapse rate. Current challenges in the identification of immunotherapy targets arise from patient-specific blast immunophenotypes and their change during disease progression. To overcome this, we present a new computational research tool to rapidly identify malignant cells. We generated single-cell flow cytometry profiles of 21 pediatric AML patients with matched samples at diagnosis, remission, and relapse. We coupled a classifier to an autoencoder for anomaly detection and classified malignant blasts with 90% accuracy. Moreover, our method assigns a developmental stage to blasts at the single-cell level, improving current classification approaches based on differentiation of the dominant phenotype. We observed major immunophenotype and developmental stage alterations between diagnosis and relapse. Patients with KMT2A rearrangement had more profound changes in their blast immunophenotypes at relapse compared to patients with other molecular features. Our method provides new insights into the immunophenotypic composition of AML blasts in an unbiased fashion and can help to define immunotherapy targets that might improve personalized AML treatment.

Artificial antigen-presenting cell system reveals CD28's role in modulating T cell functions during human immunodeficiency virus infection

T cell immune dysfunction is a prominent feature of chronic HIV infection. To evaluate non-specific dysfunction, a method involving both generic activation and T cell receptor (TCR) stimulation is necessary. We created a tunable artificial antigen-presenting cell (aAPC) system. This system consists of lipid bilayers on cytometry-compatible silica microbeads (5 μm). When only anti-CD3 is incorporated, T cell activation is limited. Introducing anti-CD28 agonists significantly elevates the cytokine expression and upregulation of activation-induced markers. CD28 co-stimulation modulates the response profile, preferentially promoting IL-2 expression relative to other cytokines. aAPCs-stimulated CD4+ and CD8+ T cells from untreated HIV-infected individuals exhibit altered effector functions and diminished CD28 dependence. These functions are skewed toward TNFα, IFNγ and CD107a, with reduced IL-2. Antiretroviral therapy partially normalizes this distorted profile in CD4+ T cells, but not in CD8+ T cells. Our findings show T cell intrinsic biases that may contribute to persistent systemic T cell dysfunction associated with HIV pathogenesis.

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

Lineage plasticity is a hallmark of cancer progression that impacts therapy outcomes, yet the mechanisms mediating this process remain unclear. Here, we introduce a versatile in vivo platform to interrogate neuroendocrine lineage transformation throughout prostate cancer progression. Transplanted mouse prostate organoids with human-relevant driver mutations (Rb1-/-; Trp53-/-; cMyc+ or Pten-/-; Trp53-/-; cMyc+) develop adenocarcinomas, but only those with Rb1 deletion advance to aggressive, ASCL1+ neuroendocrine prostate cancer (NEPC) resistant to androgen receptor signaling inhibitors. Notably, this transition requires an in vivo microenvironment not replicated by conventional organoid culture. Using multiplexed immunofluorescence and spatial transcriptomics, we reveal that ASCL1+ cells arise from KRT8+ luminal cells, progressing into transcriptionally heterogeneous ASCL1+;KRT8- NEPC. Ascl1 loss in established NEPC causes transient regression followed by recurrence, but its deletion before transplantation abrogates lineage plasticity, resulting in castration-sensitive adenocarcinomas. This dynamic model highlights the importance of therapy timing and offers a platform to identify additional lineage plasticity drivers.

A Bayesian framework for inferring dynamic intercellular interactions from time-series single-cell data

Characterizing cell-cell communication and tracking its variability over time are crucial for understanding the coordination of biological processes mediating normal development, disease progression, and responses to perturbations such as therapies. Existing tools fail to capture time-dependent intercellular interactions and primarily rely on databases compiled from limited contexts. We introduce DIISCO, a Bayesian framework designed to characterize the temporal dynamics of cellular interactions using single-cell RNA-sequencing data from multiple time points. Our method utilizes structured Gaussian process regression to unveil time-resolved interactions among diverse cell types according to their coevolution and incorporates prior knowledge of receptor-ligand complexes. We show the interpretability of DIISCO in simulated data and new data collected from T cells cocultured with lymphoma cells, demonstrating its potential to uncover dynamic cell-cell cross talk.

Envafolimab plus lenvatinib and transcatheter arterial chemoembolization for unresectable hepatocellular carcinoma: a prospective, single-arm, phase II study

Evidences regarding the feasibility of transcatheter arterial chemoembolization (TACE)-based therapy for unresectable hepatocellular carcinoma (uHCC) remains limited. This study aimed to investigate the efficacy and safety of TACE combined with envafolimab and lenvatinib for uHCC. Eligible patients with uHCC received envafolimab and lenvatinib after TACE until disease progression, conversion to surgery, intolerable toxicities, or death. The primary endpoint was the objective response rate (ORR) assessed according to Response Evaluation Criteria in Solid Tumors (RECIST) 1.1 criteria. Between March 2022 and July 2022, 38 patients were included for safety analysis, and 36 patients were included for efficacy analysis. As of the data cutoff (13 December 2023), the median follow-up was 16.9 months. The ORR was 50%, and disease control rate (DCR) was 83.3% per RECIST 1.1 (ORR and DCR of both 83.3% per modified RECIST (mRECIST)). The median progression-free survival (PFS) was 7.58 months. Of 36 patients, 17 patients were converted to resectable HCC with a surgical conversion rate of 47.2%, and 16 patients underwent surgery with R0 resection rate of 100%, pathologic complete response (pCR) rate of 31.3%. Overall incidences of treatment-related adverse events (TRAEs) of any grade was 97.4%. Grade ≥ 3 TRAEs were observed in 52.6% patients. No treatment-related deaths occurred. Image mass cytometry (IMC) analysis revealed that combined treatment improved the immune status of the tumor microenvironment, and resident macrophages had the potential to predict efficacy of this treatment. Envafolimab plus lenvatinib and TACE yielded promising survival outcomes and conversion efficiency with a tolerable safety profile. Trial registration Clinical trials: NCT05213221.

Phenotypic, transcriptomic, and spatial characterization of CD45RB+ naïve mature B cells: Implications in Sjögren's disease

The conventional classification of mature B cells overlooks the diversity within IgD+ CD27- naïve B cells. Here, to identify distinct mature naïve B cells, we categorized CD45RBMEM55- B cells (NA RB-) and CD45RBMEM55+ B cells (NA RB+) and explore their function and localization in circulation and tissues under physiological and pathological conditions. NA RB+ B cells, found in secondary lymphoid organs, differentiate into plasmablasts and secrete IgM. In Sjögren's disease, their numbers decrease, and they show over-activation and abnormal migration, suggesting an adaptive disease response. NA RB+ B cells also appear in inflamed salivary glands, indicating involvement in local immune responses. These findings highlight the distinct roles of NA RB+ B cells in health and Sjögren's disease.

Autoantigen-specific CD4+ T cells acquire an exhausted phenotype and persist in human antigen-specific autoimmune diseases

Pro-inflammatory autoantigen-specific CD4+ T helper (auto-Th) cells are central orchestrators of autoimmune diseases (AIDs). We aimed to characterize these cells in human AIDs with defined autoantigens by combining human leukocyte antigen (HLA)-tetramer-based and activation-based multidimensional ex vivo analyses. In aquaporin4-antibody-positive neuromyelitis optica spectrum disorder (AQP4-NMOSD) patients, auto-Th cells expressed CD154, but proliferative capacity and pro-inflammatory cytokines were strongly reduced. Instead, exhaustion-associated co-inhibitory receptors were expressed together with FOXP3, the canonical regulatory T cell (Treg) transcription factor. Auto-Th cells responded in vitro to checkpoint inhibition and provided potent B cell help. Cells with the same exhaustion-like (ThEx) phenotype were identified in soluble liver antigen (SLA)-antibody-autoimmune hepatitis and BP180-antibody-positive bullous pemphigoid, AIDs of the liver and skin, respectively. While originally described in cancer and chronic infection, our data point to T cell exhaustion as a common mechanism of adaptation to chronic (self-)stimulation across AID types and link exhausted CD4+ T cells to humoral autoimmune responses, with implications for therapeutic targeting.

Precursor central memory versus effector cell fate and naïve CD4+ T cell heterogeneity

Upon antigenic stimulation, naïve CD4+ T cells can give rise to phenotypically distinct effector T helper cells and long-lived memory T cells. We computationally reconstructed the in vivo trajectory of CD4+ T cell differentiation during a type I inflammatory immune response and identified two distinct differentiation paths for effector and precursor central memory T cells arising directly from naïve CD4+ T cells. Unexpectedly, our studies revealed heterogeneity among naïve CD4+ T cells, which are typically considered homogeneous save for their diverse T cell receptor usage. Specifically, a previously unappreciated population of naïve CD4+ T cells sensing environmental type I IFN exhibited distinct activation thresholds, suggesting that naïve CD4+ T cell differentiation potential may be influenced by environmental cues. This population was expanded in human viral infection and type I IFN response-lined autoimmunity. Understanding the relevance of naïve T cell heterogeneity to beneficial and maladaptive T cell responses may have therapeutic implications for adoptive T cell therapies in cancer immunotherapy and vaccination.

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.

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

The tumor microenvironment is widely recognized for its central role in driving cancer progression and influencing prognostic outcomes. There have been increasing efforts dedicated to characterizing this complex and heterogeneous environment, including developing potential prognostic tools by leveraging modern deep learning methods. However, the identification of generalizable data-driven biomarkers has been limited, in part due to the inability to interpret the complex, black-box predictions made by these models. In this study, we introduce a data-driven yet interpretable 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 and show that the identified tumor microenvironment patterns result in a risk stratification system that provides new complementary information with respect to standard stratification systems. Our results, which are validated in two independent cohorts, allow for new insights into the prognostic implications of the breast tumor microenvironment. This methodology could be applied to other cancer types more generally, providing insights into the cellular patterns of organization associated with different outcomes.

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.

The CD8+ T cell content of transbronchial biopsies from patients with a first episode of clinically stable grade A1 cellular rejection is associated with future chronic lung allograft dysfunction

BACKGROUND

T cells drive acute cellular rejection (ACR) and its progression to chronic lung allograft dysfunction (CLAD) following lung transplantation. International Society for Heart and Lung Transplantation grade A1 ACR without associated allograft dysfunction is often untreated, yet some patients develop progressive graft dysfunction. T-cell composition of A1 ACR lesions may have prognostic value; therefore, protein-level and epigenetic techniques were applied to transbronchial biopsy tissue to determine whether differential T-cell infiltration in recipients experiencing a first episode of stable grade A1 ACR (StA1R) is associated with early CLAD.

METHODS

Sixty-two patients experiencing a first episode of StA1R were divided into those experiencing CLAD within 2 years (n = 13) and those remaining CLAD-free for 5 or more years (n = 49). Imaging mass cytometry (IMC) was used to profile the spectrum and distribution of intragraft T cell phenotypes on a subcohort (n = 16; 8 early-CLAD and 8 no early-CLAD). Immunofluorescence was used to quantify CD4+, CD8+, and FOXP3+ cells. Separately, CD3+ cells were fluorescently labeled, micro-dissected, and the degree of Treg-specific demethylated region methylation was determined.

RESULTS

PhenoGraph unsupervised clustering on IMC revealed 50 unique immune cell subpopulations. Methylation and immunofluorescence analyses demonstrated no significant differences in Tregs between early-CLAD and no early-CLAD groups. Immunofluorescence revealed that patients who developed CLAD within 2 years of lung transplantation showed greater CD8+ T cell infiltration compared to those who remained CLAD-free for 5 or more years.

CONCLUSIONS

In asymptomatic patients with a first episode of A1 rejection, greater CD8+ T cell content may be indicative of worse long-term outlook.

Global and single-cell proteomics view of the co-evolution between neural progenitors and breast cancer cells in a co-culture model

BACKGROUND

Presence of nerves in tumours, by axonogenesis and neurogenesis, is gaining increased attention for its impact on cancer initiation and development, and the new field of cancer neuroscience is emerging. A recent study in prostate cancer suggested that the tumour microenvironment may influence cancer progression by recruitment of Doublecortin (DCX)-expressing neural progenitor cells (NPCs). However, the presence of such cells in human breast tumours has not been comprehensively explored.

METHODS

Here, we investigate the presence of DCX-expressing cells in breast cancer stromal tissue from patients using Imaging Mass Cytometry. Single-cell analysis of 372,468 cells across histopathological images of 107 breast cancers enabled spatial resolution of neural elements in the stromal compartment in correlation with clinicopathological features of these tumours. In parallel, we established a 3D in vitro model mimicking breast cancer neural progenitor-innervation and examined the two cell types as they co-evolved in co-culture by using mass spectrometry-based global proteomics.

FINDINGS

Stromal presence of DCX + cells is associated with tumours of higher histological grade, a basal-like phenotype, and shorter patient survival in tumour tissue from patients with breast cancer. Global proteomics analysis revealed significant changes in the proteomic landscape of both breast cancer cells and neural progenitors in co-culture.

INTERPRETATION

These results support that neural involvement plays an active role in breast cancer and warrants further studies on the relevance of nerve elements for tumour progression.

FUNDING

This work was supported by the Research Council of Norway through its Centre of Excellence funding scheme, project number 223250 (to L.A.A), the Norwegian Cancer Society (to L.A.A. and H.V.), the Regional Health Trust Western Norway (Helse Vest) (to L.A.A.), the Meltzer Research Fund (to H.V.) and the National Institutes of Health (NIH)/NIGMS grant R01 GM132129 (to J.A.P.).

Neonatal inflammation impairs developmentally-associated microglia and promotes a highly reactive microglial subset

Microglia and border-associated macrophages play critical roles in both immunity and neurodevelopment. The disruption of microglial development trajectories by neonatal inflammation is an important issue in research on neurodevelopmental disorders (NDDs), as models have suggested a strong association between inflammation and cognitive deficits. Here, we explored by single-cell RNA sequencing and flow cytometry the impact of neonatal inflammation in a mouse NDD model on brain myeloid cell subsets. A specific subset of microglia expressing the complement receptor C5ar1 has been identified, in which inflammatory pathways are most strongly activated. Based on transcriptional similarity, this subset appears to originate from the most mature and "homeostatic" microglia at this stage of development and demonstrated hypersensitivity to inflammation. Besides that, Spp1-microglia supporting oligodendrocyte differentiation, primitive and proliferative microglia were reduced by inflammation. These findings suggest major changes in microglial subsets developmental trajectories and reactivity contributing to NDDs induced by neonatal inflammation.

Protocol for high-plex, whole-slide imaging of human formalin-fixed paraffin-embedded tissue using PhenoCycler-Fusion

Single-cell spatial analysis of proteins is rapidly becoming increasingly important in revealing biological insights. Here, we present a protocol for automated high-plex multi-slide immunofluorescence staining and imaging of human head and neck cancer formalin-fixed paraffin-embedded (FFPE) sections using PhenoCycler-Fusion 2.0 technology. We describe steps for preparing human head and neck cancer FFPE tissues, staining with a panel of immunophenotyping markers, and Flow Cell assembly. We then detail procedures for setting up for a PhenoCycler-Fusion run, post-run Flow Cell removal, and downstream analyses. For complete details on the use and execution of this protocol, please refer to Jhaveri et al.1.

CelltypeR: A flow cytometry pipeline to characterize single cells from brain organoids

Motivated by the cellular heterogeneity in complex tissues, particularly in brain and induced pluripotent stem cell (iPSC)-derived brain models, we developed a complete workflow to reproducibly characterize cell types in complex tissues. Our approach combines a flow cytometry (FC) antibody panel with our computational pipeline CelltypeR, enabling dataset aligning, unsupervised clustering optimization, cell type annotating, and statistical comparisons. Applied to human iPSC derived midbrain organoids, it successfully identified the major brain cell types. We performed fluorescence-activated cell sorting of CelltypeR-defined astrocytes, radial glia, and neurons, exploring transcriptional states by single-cell RNA sequencing. Among the sorted neurons, we identified subgroups of dopamine neurons: one reminiscent of substantia nigra cells most vulnerable in Parkinson's disease. Finally, we used our workflow to track cell types across a time course of organoid differentiation. Overall, our adaptable analysis framework provides a generalizable method for reproducibly identifying cell types across FC datasets in complex tissues.

Spatiotemporal metabolomic approaches to the cancer-immunity panorama: a methodological perspective

Metabolic reprogramming drives the development of an immunosuppressive tumor microenvironment (TME) through various pathways, contributing to cancer progression and reducing the effectiveness of anticancer immunotherapy. However, our understanding of the metabolic landscape within the tumor-immune context has been limited by conventional metabolic measurements, which have not provided comprehensive insights into the spatiotemporal heterogeneity of metabolism within TME. The emergence of single-cell, spatial, and in vivo metabolomic technologies has now enabled detailed and unbiased analysis, revealing unprecedented spatiotemporal heterogeneity that is particularly valuable in the field of cancer immunology. This review summarizes the methodologies of metabolomics and metabolic regulomics that can be applied to the study of cancer-immunity across single-cell, spatial, and in vivo dimensions, and systematically assesses their benefits and limitations.

Flexible and robust cell type annotation for highly multiplexed tissue images

Identifying cell types in highly multiplexed images is essential for understanding tissue spatial organization. Current cell type annotation methods often rely on extensive reference images and manual adjustments. In this work, we present a tool, Robust Image-Based Cell Annotator (RIBCA), that enables accurate, automated, unbiased, and fine-grained cell type annotation for images with a wide range of antibody panels, without requiring additional model training or human intervention. Our tool has successfully annotated over 1 million cells, revealing the spatial organization of various cell types across more than 40 different human tissues. It is open-source and features a modular design, allowing for easy extension to additional cell types.

ESCHR: a hyperparameter-randomized ensemble approach for robust clustering across diverse datasets

Clustering is widely used for single-cell analysis, but current methods are limited in accuracy, robustness, ease of use, and interpretability. To address these limitations, we developed an ensemble clustering method that outperforms other methods at hard clustering without the need for hyperparameter tuning. It also performs soft clustering to characterize continuum-like regions and quantify clustering uncertainty, demonstrated here by mapping the connectivity and intermediate transitions between MNIST handwritten digits and between hypothalamic tanycyte subpopulations. This hyperparameter-randomized ensemble approach improves the accuracy, robustness, ease of use, and interpretability of single-cell clustering, and may prove useful in other fields as well.

SMAC mimetic drives microglia phenotype and glioblastoma immune microenvironment

Tumor-associated macrophages/microglia (TAMs) are highly plastic and heterogeneous immune cells that can be immune-supportive or tumor-supportive depending of the microenvironment. TAMs are the most abundant immune cells in glioblastoma (GB), and play a key role in immunosuppression. Therefore, TAMs reprogramming toward immune-supportive cells is a promising strategy to overcome immunosuppression. By leveraging scRNAseq human GB databases, we identified that Inhibitor of Apoptosis Proteins (IAP) were expressed by TAMs. To investigate their role in TAMs-related immunosuppression, we antagonized IAP using the central nervous system permeant SMAC mimetic GDC-0152 (SMg). On explants and cultured immune cells isolated from human GB samples, SMg modified TAMs activity. We showed that SMg treatment promoted microglia pro-apoptotic and anti-tumoral function via caspase-3 pro-inflammatory cleavage and the inhibition of tumoroids growth. Then we designed a relevant immunogenic mouse GB model to decipher the spatio-temporal densities, distribution, phenotypes and function of TAMs with or without SMg treatment. We used 3D imaging techniques, a transgenic mouse with fluorescent TAM subsets and mass cytometry. We confirmed that SMg promoted microglia activation, antigen-presenting function and tumor infiltration. In addition, we observed a remodeling of blood vessels, a decrease in anti-inflammatory macrophages and an increased level of monocytes and their mo-DC progeny. This remodeling of the TAM landscape is associated with an increase in CD8 T cell density and activation. Altogether, these results demonstrated that SMg drives the immunosuppressive basal microglia toward an active phenotype with pro-apoptotic and anti-tumoral function and modifies the GB immune landscape. This identifies IAP as targets of choice for a potential mechanism-based therapeutic strategy and SMg as a promising molecule for this application.

scBubbletree: computational approach for visualization of single cell RNA-seq data

BACKGROUND

Visualization approaches transform high-dimensional data from single cell RNA sequencing (scRNA-seq) experiments into two-dimensional plots that are used for analysis of cell relationships, and as a means of reporting biological insights. Yet, many standard approaches generate visuals that suffer from overplotting, lack of quantitative information, and distort global and local properties of biological patterns relative to the original high-dimensional space.

RESULTS

We present scBubbletree, a new, scalable method for visualization of scRNA-seq data. The method identifies clusters of cells of similar transcriptomes and visualizes such clusters as "bubbles" at the tips of dendrograms (bubble trees), corresponding to quantitative summaries of cluster properties and relationships. scBubbletree stacks bubble trees with further cluster-associated information in a visually easily accessible way, thus facilitating quantitative assessment and biological interpretation of scRNA-seq data. We demonstrate this with large scRNA-seq data sets, including one with over 1.2 million cells.

CONCLUSIONS

To facilitate coherent quantification and visualization of scRNA-seq data we developed the R-package scBubbletree, which is freely available as part of the Bioconductor repository at: https://bioconductor.org/packages/scBubbletree/.

Role of circulating T follicular helper subsets following Ty21a immunization and oral challenge with wild type S. Typhi in humans

Despite decades of intense research, our understanding of the correlates of protection against Salmonella Typhi (S. Typhi) infection and disease remains incomplete. T follicular helper cells (TFH), an important link between cellular and humoral immunity, play an important role in the development and production of high affinity antibodies. While traditional TFH cells reside in germinal centers, circulating TFH (cTFH) (a memory subset of TFH) are present in blood. We used specimens from a typhoid controlled human infection model whereby participants were immunized with Ty21a live attenuated S. Typhi vaccine and then challenged with virulent S. Typhi. Some participants developed typhoid disease (TD) and some did not (NoTD), which allowed us to assess the association of cTFH subsets in the development and prevention of typhoid disease. Of note, the frequencies of cTFH were higher in NoTD than in TD participants, particularly 7 days after challenge. Furthermore, the frequencies of cTFH2 and cTFH17, but not cTFH1 subsets were higher in NoTD than TD participants. However, we observed that ex-vivo expression of activation and homing markers were higher in TD than in NoTD participants, particularly after challenge. Moreover, cTFH subsets produced higher levels of S. Typhi-specific responses (cytokines/chemokines) in both the immunization and challenge phases. Interestingly, unsupervised analysis revealed unique clusters with distinct signatures for each cTFH subset that may play a role in either the development or prevention of typhoid disease. Importantly, we observed associations between frequencies of defined cTFH subsets and anti-S. Typhi antibodies. Taken together, our results suggest that circulating TFH2 and TFH17 subsets might play an important role in the development or prevention of typhoid disease. The contribution of these clusters was found to be distinct in the immunization and/or challenge phases. These results have important implications for vaccines aimed at inducing long-lived protective T cell and antibody responses.

Distinct tumor architectures and microenvironments for the initiation of breast cancer metastasis in the brain

Brain metastasis, a serious complication of cancer, hinges on the initial survival, microenvironment adaptation, and outgrowth of disseminated cancer cells. To understand the early stages of brain colonization, we investigated two prevalent sources of cerebral relapse, triple-negative (TNBC) and HER2+ (HER2BC) breast cancers. Using mouse models and human tissue samples, we found that these tumor types colonize the brain, with a preference for distinctive tumor architectures, stromal interfaces, and autocrine programs. TNBC models tend to form perivascular sheaths with diffusive contact with astrocytes and microglia. In contrast, HER2BC models tend to form compact spheroids driven by autonomous tenascin C production, segregating stromal cells to the periphery. Single-cell transcriptomics of the tumor microenvironment revealed that these architectures evoke differential Alzheimer's disease-associated microglia (DAM) responses and engagement of the GAS6 receptor AXL. The spatial features of the two modes of brain colonization have relevance for leveraging the stroma to treat brain metastasis.

A spatially resolved single-cell lung atlas integrated with clinical and blood signatures distinguishes COVID-19 disease trajectories

COVID-19 is characterized by a broad range of symptoms and disease trajectories. Understanding the correlation between clinical biomarkers and lung pathology during acute COVID-19 is necessary to understand its diverse pathogenesis and inform more effective treatments. Here, we present an integrated analysis of longitudinal clinical parameters, peripheral blood markers, and lung pathology in 142 Brazilian patients hospitalized with COVID-19. We identified core clinical and peripheral blood signatures differentiating disease progression between patients who recovered from severe disease compared with those who succumbed to the disease. Signatures were heterogeneous among fatal cases yet clustered into two patient groups: "early death" (<15 days until death) and "late death" (>15 days). Progression to early death was characterized systemically and in lung histopathological samples by rapid endothelial and myeloid activation and the presence of thrombi associated with SARS-CoV-2+ macrophages. In contrast, progression to late death was associated with fibrosis, apoptosis, and SARS-CoV-2+ epithelial cells in postmortem lung tissue. In late death cases, cytotoxicity, interferon, and T helper 17 (TH17) signatures were only detectable in the peripheral blood after 2 weeks of hospitalization. Progression to recovery was associated with higher lymphocyte counts, TH2 responses, and anti-inflammatory-mediated responses. By integrating antemortem longitudinal blood signatures and spatial single-cell lung signatures from postmortem lung samples, we defined clinical parameters that could be used to help predict COVID-19 outcomes.

Interleukin-2 immunotherapy reveals human regulatory T cell subsets with distinct functional and tissue-homing characteristics

Due to its stimulatory potential for immunomodulatory CD4+ regulatory T (Treg) cells, low-dose interleukin-2 (IL-2) immunotherapy has gained considerable attention for the treatment of autoimmune diseases. In this investigator-initiated single-arm non-placebo-controlled phase-2 clinical trial of low-dose IL-2 immunotherapy in systemic lupus erythematosus (SLE) patients, we generated a comprehensive atlas of in vivo human immune responses to low-dose IL-2. We performed an in-depth study of circulating and cutaneous immune cells by imaging mass cytometry, high-parameter flow cytometry, transcriptomics, and targeted serum proteomics. Low-dose IL-2 stimulated various circulating immune cells, including Treg cells with a skin-homing phenotype that appeared in the skin of SLE patients in close interaction with endothelial cells. Analysis of surface proteins and transcriptomes revealed different IL-2-driven Treg cell activation programs, including gut-homing CD38+, skin-homing HLA-DR+, and highly proliferative inflammation-homing CD38+ HLA-DR+ Treg cells. Collectively, these data define the distinct human Treg cell subsets that are responsive to IL-2 immunotherapy.

GPR37 processing in neurodegeneration: a potential marker for Parkinson's Disease progression rate

The orphan G protein-coupled receptor 37 (GPR37), widely associated with Parkinson's disease (PD), undergoes proteolytic processing under physiological conditions. The N-terminus domain is proteolyzed by a disintegrin and metalloproteinase 10 (ADAM-10), which generates various membrane receptor forms and ectodomain shedding (ecto-GPR37) in the extracellular environment. We investigated the processing and density of GPR37 in several neurodegenerative conditions, including Lewy body disease (LBD), multiple system atrophy (MSA), corticobasal degeneration (CBD), progressive supranuclear palsy (PSP), and Alzheimer's disease (AD). The presence of ecto-GPR37 peptides in the cerebrospinal fluid (CSF) of PD, MSA, CBD and PSP patients was assessed through an in-house nanoluciferase-based immunoassay. This study identified increased receptor processing in early-stage LBD within the PFC and striatum, key brain areas in neurodegeneration. In MSA only the 52 kDa form of GPR37 appeared in the striatum. This form was also significantly elevated in the striatum of AD necropsies. On the contrary, GPR37 processing remained unchanged in the brains of CBD and PSP patients. Furthermore, while CSF ecto-GPR37 increased in PD patients, its levels remained unchanged in MSA, CBD, and PSP subjects. Importantly, patients with PD with rapid progression of the disease did not have elevated ecto-GPR37 in the CSF, while those with slow progression showed a significant increase, suggesting a possible prognostic use of ecto-GPR37 in PD. This research underscores the distinctive processing and density patterns of GPR37 in neurodegenerative diseases, providing crucial insights into its potential role as an indicator of PD progression rates.

A novel GFP-based strategy to quantitate cellular spatial associations in HSV-1 viral pathogenesis

Periodic reactivation of herpes simplex virus type 1 (HSV-1) triggers immune responses that result in corneal scarring (CS), known as herpes stromal keratitis (HSK). Despite considerable research, fully understanding HSK and eliminating it remains challenging due to a lack of comprehensive analysis of HSV-1-infected immune cells in both corneas and trigeminal ganglia (TG). We engineered a recombinant HSV-1 expressing green fluorescent protein (GFP) in the virulent McKrae virus strain that does not require corneal scarification for efficient virus replication (GFP-McKrae). Next-generation sequencing (NGS) analysis, along with in vitro and in vivo assays, showed that GFP-McKrae virus was similar to WT-McKrae virus. Furthermore, corneal cells infected with GFP-McKrae were quantitatively analyzed using image mass cytometry (IMC). The single-cell reconstruction data generated cellular maps of corneas based on the expression of 25 immune cell markers in GFP-McKrae-infected mice. Corneas from mock control mice showed the presence of T cells and macrophages, whereas corneas from GFP-McKrae-infected mice on days 3 and 5 post-infection (PI) exhibited increased immune cells. Notably, on day 3 PI, increased GFP expression was observed in closely situated clusters of DCs, macrophages, and epithelial cells. By day 5 PI, macrophages and T cells became prominent. Finally, immunostaining methods detected HSV-1 or GFP and gD proteins in latently infected TG. This study presents a valuable strategy for identifying cellular spatial associations in viral pathogenesis and holds promise for future therapeutic applications.IMPORTANCEThe goal of this study was to establish quantitative approaches to analyze immune cell markers in HSV-1-infected intact corneas and trigeminal ganglia from primary and latently infected mice. This allowed us to define spatial and temporal interactions between specific immune cells and their potential roles in virus replication and latency. To accomplish this important goal, we took advantage of the utility of GFP-McKrae virus as a valuable research tool while also highlighting its potential to uncover previously unrecognized cell types that play pivotal roles in HSV-1 replication and latency. Such insights will pave the way for developing targeted therapeutic approaches to tackle HSV-1 infections more effectively.

ImmCellTyper facilitates systematic mass cytometry data analysis for deep immune profiling

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

Recipient tissue microenvironment determines developmental path of intestinal innate lymphoid progenitors

Innate lymphoid cells (ILCs) are critical in maintaining tissue homeostasis, and during infection and inflammation. Here we identify, by using combinatorial reporter mice, a rare ILC progenitor (ILCP) population, resident to the small intestinal lamina propria (siLP) in adult mice. Transfer of siLP-ILCP into recipients generates group 1 ILCs (including ILC1 and NK cells), ILC2s and ILC3s within the intestinal microenvironment, but almost exclusively group 1 ILCs in the liver, lung and spleen. Single cell gene expression analysis and high dimensional spectral cytometry analysis of the siLP-ILCPs and ILC progeny indicate that the phenotype of the group 1 ILC progeny is also influenced by the tissue microenvironment. Thus, a local pool of siLP-ILCP can contribute to pan-ILC generation in the intestinal microenvironment but has more restricted potential in other tissues, with a greater propensity than bone marrow-derived ILCPs to favour ILC1 and ILC3 production. Therefore, ILCP potential is influenced by both tissue of origin and the microenvironment during development. This may provide additional flexibility during the tuning of immune reactions.

Single-cell and spatial sequencing identifies senescent and germinal tumor cells in adamantinomatous craniopharyngiomas

Adamantinomatous craniopharyngioma (ACP) is a clinically aggressive tumor without effective treatment method. Previous studies proposed a paracrine tumorigenesis model, in which oncogenic β-catenin induces senescence in pituitary stem cells and the senescent cells lead the formation of paracrine tumors through secretion of pro-tumorigenic factors. However, there lacks characterization on senescent cells in ACPs. Here, we profiled 12 ACPs with single-cell RNA and TCR-sequencing to elucidate the cellular atlas in ACPs and 3 of them were also subject to spatial sequencing to localize different subpopulations of the tumor cells. In total, we obtained the transcriptome profiles of 70,682 cells. Tumor cells, which were unambiguously identified through the cellular mutation status of the driver CTNNB1 mutations, were clustered into 6 subsets. The whorl-like cluster (WC) cells show distinct molecular features from the other tumor cells and the palisading epithelium (PE) cells consists of a proliferating subset. Other than typical PE and WC, we identified two novel subpopulations of the tumor cells. In one subpopulation, the cells express a high level of cytokines, e.g., FDCSP and S100A8/A9, and are enriched with the senescence-associated secretory phenotype (SASP) factors. Hematoxylin and eosin staining reveals that these SASP cells lack an ordered structures and their nuclei are elongated. In the other subpopulation, the cell sizes are small and they are tightly packed together with an unusual high density expressing a high level of mitochondrial genes (median 10.9%). These cells are the origin of the tumor developmental trajectories revealed by RNA velocity and pseudo-time analysis. Single-cell RNA and TCR analysis reveals that some ACPs are infiltrated with clonally expanded cytotoxic T cells. We propose a hypothesis that WC and PE are formed via different negative regulation mechanisms of the overactivated WNT/β-catenin signaling which provides a new understanding on the tumorigenesis of ACPs. The study lays a foundation for future studies on targeting senescent cells in ACPs with senolytic compounds or other therapeutic agents.

A unique human cord blood CD8+CD45RA+CD27+CD161+ T-cell subset identified by flow cytometric data analysis using Seurat

Advances in single-cell level analytical techniques, especially cytometric approaches, have led to profound innovation in biomedical research, particularly in the field of clinical immunology. This has resulted in an expansion of high-dimensional data, posing great challenges for comprehensive and unbiased analysis. Conventional manual analysis is thus becoming untenable to handle these challenges. Furthermore, most newly developed computational methods lack flexibility and interoperability, hampering their accessibility and usability. Here, we adapted Seurat, an R package originally developed for single-cell RNA sequencing (scRNA-seq) analysis, for high-dimensional flow cytometric data analysis. Based on a 20-marker antibody panel and analyses of T-cell profiles in both adult blood and cord blood (CB), we showcased the robust capacity of Seurat in flow cytometric data analysis, which was further validated by Spectre, another high-dimensional cytometric data analysis package, and conventional manual analysis. Importantly, we identified a unique CD8+ T-cell population defined as CD8+CD45RA+CD27+CD161+ T cell that was predominantly present in CB. We characterised its IFN-γ-producing and potential cytotoxic properties using flow cytometry experiments and scRNA-seq analysis from a published dataset. Collectively, we identified a unique human CB CD8+CD45RA+CD27+CD161+ T-cell subset and demonstrated that Seurat, a widely used package for scRNA-seq analysis, possesses great potential to be repurposed for cytometric data analysis. This facilitates an unbiased and thorough interpretation of complicated high-dimensional data using a single analytical pipeline and opens a novel avenue for data-driven investigation in clinical immunology.

In vivo CRISPR screens identify a dual function of MEN1 in regulating tumor-microenvironment interactions

Functional genomic screens in two-dimensional cell culture models are limited in identifying therapeutic targets that influence the tumor microenvironment. By comparing targeted CRISPR-Cas9 screens in a two-dimensional culture with xenografts derived from the same cell line, we identified MEN1 as the top hit that confers differential dropout effects in vitro and in vivo. MEN1 knockout in multiple solid cancer types does not impact cell proliferation in vitro but significantly promotes or inhibits tumor growth in immunodeficient or immunocompetent mice, respectively. Mechanistically, MEN1 knockout redistributes MLL1 chromatin occupancy, increasing H3K4me3 at repetitive genomic regions, activating double-stranded RNA expression and increasing neutrophil and CD8+ T cell infiltration in immunodeficient and immunocompetent mice, respectively. Pharmacological inhibition of the menin-MLL interaction reduces tumor growth in a CD8+ T cell-dependent manner. These findings reveal tumor microenvironment-dependent oncogenic and tumor-suppressive functions of MEN1 and provide a rationale for targeting MEN1 in solid cancers.

Comprehensive Characterization of a Subfamily of Ca2+-Binding Proteins in Mouse and Human Retinal Neurons at Single-Cell Resolution

Ca2+-binding proteins (CaBPs; CaBP1-5) are a subfamily of neuronal Ca2+ sensors with high homology to calmodulin. Notably, CaBP4, which is exclusively expressed in rod and cone photoreceptors, is crucial for maintaining normal retinal functions. However, the functional roles of CaBP1, CaBP2, and CaBP5 in the retina remain elusive, primarily due to limited understanding of their expression patterns within inner retinal neurons. In this study, we conducted a comprehensive transcript analysis using single-cell RNA sequencing datasets to investigate the gene expression profiles of CaBPs in mouse and human retinal neurons. Our findings revealed notable similarities in the overall expression patterns of CaBPs across both species. Specifically, nearly all amacrine cell, ganglion cell, and horizontal cell types exclusively expressed CaBP1. In contrast, the majority of bipolar cell types, including rod bipolar (RB) cells, expressed distinct combinations of CaBP1, CaBP2, and CaBP5, rather than a single CaBP as previously hypothesized. Remarkably, mouse rods and human cones exclusively expressed CaBP4, whereas mouse cones and human rods coexpressed both CaBP4 and CaBP5. Our single-cell reverse transcription polymerase chain reaction analysis confirmed the coexpression CaBP1 and CaBP5 in individual RBs from mice of either sex. Additionally, all three splice variants of CaBP1, primarily L-CaBP1, were detected in mouse RBs. Taken together, our study offers a comprehensive overview of the distribution of CaBPs in mouse and human retinal neurons, providing valuable insights into their roles in visual functions.

Spatial host-microbiome sequencing reveals niches in the mouse gut

Mucosal and barrier tissues, such as the gut, lung or skin, are composed of a complex network of cells and microbes forming a tight niche that prevents pathogen colonization and supports host-microbiome symbiosis. Characterizing these networks at high molecular and cellular resolution is crucial for understanding homeostasis and disease. Here we present spatial host-microbiome sequencing (SHM-seq), an all-sequencing-based approach that captures tissue histology, polyadenylated RNAs and bacterial 16S sequences directly from a tissue by modifying spatially barcoded glass surfaces to enable simultaneous capture of host transcripts and hypervariable regions of the 16S bacterial ribosomal RNA. We applied our approach to the mouse gut as a model system, used a deep learning approach for data mapping and detected spatial niches defined by cellular composition and microbial geography. We show that subpopulations of gut cells express specific gene programs in different microenvironments characteristic of regional commensal bacteria and impact host-bacteria interactions. SHM-seq should enhance the study of native host-microbe interactions in health and disease.

scATAcat: cell-type annotation for scATAC-seq data

Cells whose accessibility landscape has been profiled with scATAC-seq cannot readily be annotated to a particular cell type. In fact, annotating cell-types in scATAC-seq data is a challenging task since, unlike in scRNA-seq data, we lack knowledge of 'marker regions' which could be used for cell-type annotation. Current annotation methods typically translate accessibility to expression space and rely on gene expression patterns. We propose a novel approach, scATAcat, that leverages characterized bulk ATAC-seq data as prototypes to annotate scATAC-seq data. To mitigate the inherent sparsity of single-cell data, we aggregate cells that belong to the same cluster and create pseudobulk. To demonstrate the feasibility of our approach we collected a number of datasets with respective annotations to quantify the results and evaluate performance for scATAcat. scATAcat is available as a python package at https://github.com/aybugealtay/scATAcat.

Macrophage Migration Inhibitory Factor Promotes Thromboinflammation and Predicts Fast Progression of Aortic Stenosis

BACKGROUND

Aortic stenosis (AS) is driven by progressive inflammatory and fibrocalcific processes regulated by circulating inflammatory and valve resident endothelial and interstitial cells. The impact of platelets, platelet-derived mediators, and platelet-monocyte interactions on the acceleration of local valvular inflammation and mineralization is presently unknown.

METHODS

We prospectively enrolled 475 consecutive patients with severe symptomatic AS undergoing aortic valve replacement. Clinical workup included repetitive echocardiography, analysis of platelets, monocytes, chemokine profiling, aortic valve tissue samples for immunohistochemistry, and gene expression analysis.

RESULTS

The patients were classified as fast-progressive AS by the median ∆Vmax of 0.45 m/s per year determined by echocardiography. Immunohistological aortic valve analysis revealed enhanced cellularity in fast-progressive AS (slow- versus fast-progressive AS; median [interquartile range], 247 [142.3-504] versus 717.5 [360.5-1234]; P<0.001) with less calcification (calcification area, mm2: 33.74 [27.82-41.86] versus 20.54 [13.52-33.41]; P<0.001). MIF (macrophage migration inhibitory factor)-associated gene expression was significantly enhanced in fast-progressive AS accompanied by significantly elevated MIF plasma levels (mean±SEM; 6877±379.1 versus 9959±749.1; P<0.001), increased platelet activation, and decreased intracellular MIF expression indicating enhanced MIF release upon platelet activation (CD62P, %: median [interquartile range], 16.8 [11.58-23.8] versus 20.55 [12.48-32.28], P=0.005; MIF, %: 4.85 [1.48-9.75] versus 2.3 [0.78-5.9], P<0.001). Regression analysis confirmed that MIF-associated biomarkers are strongly associated with an accelerated course of AS.

CONCLUSIONS

Our findings suggest a key role for platelet-derived MIF and its interplay with circulating and valve resident monocytes/macrophages in local and systemic thromboinflammation during accelerated AS. MIF-based biomarkers predict an accelerated course of AS and represent a novel pharmacological target to attenuate progression of AS.

A cell state-specific metabolic vulnerability to GPX4-dependent ferroptosis in glioblastoma

Glioma cells hijack developmental programs to control cell state. Here, we uncover a glioma cell state-specific metabolic liability that can be therapeutically targeted. To model cell conditions at brain tumor inception, we generated genetically engineered murine gliomas, with deletion of p53 alone (p53) or with constitutively active Notch signaling (N1IC), a pathway critical in controlling astrocyte differentiation during brain development. N1IC tumors harbored quiescent astrocyte-like transformed cell populations while p53 tumors were predominantly comprised of proliferating progenitor-like cell states. Further, N1IC transformed cells exhibited increased mitochondrial lipid peroxidation, high ROS production and depletion of reduced glutathione. This altered mitochondrial phenotype rendered the astrocyte-like, quiescent populations more sensitive to pharmacologic or genetic inhibition of the lipid hydroperoxidase GPX4 and induction of ferroptosis. Treatment of patient-derived early-passage cell lines and glioma slice cultures generated from surgical samples with a GPX4 inhibitor induced selective depletion of quiescent astrocyte-like glioma cell populations with similar metabolic profiles. Collectively, these findings reveal a specific therapeutic vulnerability to ferroptosis linked to mitochondrial redox imbalance in a subpopulation of quiescent astrocyte-like glioma cells resistant to standard forms of treatment.

A single-cell atlas of spatial and temporal gene expression in the mouse cranial neural plate

The formation of the mammalian brain requires regionalization and morphogenesis of the cranial neural plate, which transforms from an epithelial sheet into a closed tube that provides the structural foundation for neural patterning and circuit formation. Sonic hedgehog (SHH) signaling is important for cranial neural plate patterning and closure, but the transcriptional changes that give rise to the spatially regulated cell fates and behaviors that build the cranial neural tube have not been systematically analyzed. Here we used single-cell RNA sequencing to generate an atlas of gene expression at six consecutive stages of cranial neural tube closure in the mouse embryo. Ordering transcriptional profiles relative to the major axes of gene expression predicted spatially regulated expression of 870 genes along the anterior-posterior and mediolateral axes of the cranial neural plate and reproduced known expression patterns with over 85% accuracy. Single-cell RNA sequencing of embryos with activated SHH signaling revealed distinct SHH-regulated transcriptional programs in the developing forebrain, midbrain, and hindbrain, suggesting a complex interplay between anterior-posterior and mediolateral patterning systems. These results define a spatiotemporally resolved map of gene expression during cranial neural tube closure and provide a resource for investigating the transcriptional events that drive early mammalian brain development.

The interferon/STAT1 signaling axis is a common feature of tumor-initiating cells in breast cancer

A tumor cell subpopulation of tumor-initiating cells (TIC), or "cancer stem cells", are associated with therapeutic resistance, as well as both local and distant recurrence. Enriched populations of TIC are identified by markers including aldehyde dehydrogenase (ALDH1) activity, the cell surface marker combination CD44 + /CD24 - , or fluorescent reporters for signaling pathways that regulate TIC function. We showed previously that S ignal T ransducer and A ctivator of T ranscription (STAT)-mediated transcription allows enrichment for TIC in claudin-low models of human triple-negative breast cancer using a STAT-responsive reporter. However, the molecular phenotypes of STAT TIC are not well understood, and there is no existing method to lineage-trace TIC as they undergo cell state changes. Using a new STAT-responsive lineage-tracing (LT) system in conjunction with our original reporter, we enriched for cells with enhanced mammosphere-forming potential in some, but not all, basal-like triple-negative breast cancer (TNBC) xenograft models (TNBC) indicating TIC-related and TIC-independent functions for STAT signaling. Single-cell RNA sequencing (scRNAseq) of reporter-tagged xenografts and clinical samples identified a common interferon (IFN)/STAT1-associated transcriptional state, previously linked to inflammation and macrophage differentiation, in TIC. Surprisingly, most of the genes we identified are not present in previously published TIC signatures derived using bulk RNA sequencing. Finally, we demonstrated that bone marrow stromal cell antigen 2 (BST2), is a cell surface marker of this state, and that it functionally regulates TIC frequency. These results suggest TIC may exploit the IFN/STAT1 signaling axis to promote their activity, and that targeting this pathway may help eliminate TIC.

Significance

TIC differentially express interferon response genes, which were not previously reported in bulk RNA sequencing-derived TIC signatures, highlighting the importance of coupling single-cell transcriptomics with enrichment to derive TIC signatures.

Single-cell analysis of innate spinal cord regeneration identifies intersecting modes of neuronal repair

Adult zebrafish have an innate ability to recover from severe spinal cord injury. Here, we report a comprehensive single nuclear RNA sequencing atlas that spans 6 weeks of regeneration. We identify cooperative roles for adult neurogenesis and neuronal plasticity during spinal cord repair. Neurogenesis of glutamatergic and GABAergic neurons restores the excitatory/inhibitory balance after injury. In addition, a transient population of injury-responsive neurons (iNeurons) show elevated plasticity 1 week post-injury. We found iNeurons are injury-surviving neurons that acquire a neuroblast-like gene expression signature after injury. CRISPR/Cas9 mutagenesis showed iNeurons are required for functional recovery and employ vesicular trafficking as an essential mechanism that underlies neuronal plasticity. This study provides a comprehensive resource of the cells and mechanisms that direct spinal cord regeneration and establishes zebrafish as a model of plasticity-driven neural repair.

Spatial Single Cell Profiling Using Imaging Mass Cytometry: Inflammatory Versus Penetrating Crohn's Disease

BACKGROUND AND AIMS

Fistula formation is a major complication in Crohn's disease [CD] and the role of the immune cell compartment remains to be elucidated. Thus, we compared the immune cell compartment of CD fistula to inflammatory CD colitis using imaging mass cytometry [IMC] and immunofluorescence.

METHODS

A 36-marker panel including structural, functional, and lineage markers for use in IMC was established. This panel was applied to analyse paraffin-embedded CD fistula tract [n = 11], CD colitis [n = 10], and colon samples from non-inflamed controls [n = 12]. Computational methods for cell segmentation, dimensionality reduction, and cell type clustering were used to define cell populations for cell frequency, marker distribution, and spatial neighbourhood analysis. Multiplex immunofluorescence was used for higher resolution spatial analysis.

RESULTS

Analysis of cell frequencies in CD fistulas compared to CD colitis and control colonic samples revealed a significant increase in neutrophils, effector cytotoxic T cells, and inflammatory macrophages in CD fistula samples, whereas regulatory T cells were decreased. Neutrophils in CD fistula expressed significantly more matrix metalloproteinase 9 [MMP9], correlating with extracellular matrix remodelling. Neighbourhood analysis revealed a strong association between MMP9+ neutrophils and effector cytotoxic T cells in both CD fistulas and colitis.

CONCLUSIONS

This study presents the first highly multiplexed single cell analysis of the immune cell compartment of CD fistulas and their spatial context. It links immune cell dynamics, particularly MMP9+ neutrophils, to extracellular matrix remodelling in CD fistulas, offering insights into the complex network of cellular interactions and potential therapeutic targets for CD complications.

Fasting reshapes tissue-specific niches to improve NK cell-mediated anti-tumor immunity

Fasting is associated with improved outcomes in cancer. Here, we investigated the impact of fasting on natural killer (NK) cell anti-tumor immunity. Cyclic fasting improved immunity against solid and metastatic tumors in an NK cell-dependent manner. During fasting, NK cells underwent redistribution from peripheral tissues to the bone marrow (BM). In humans, fasting also reduced circulating NK cell numbers. NK cells in the spleen of fasted mice were metabolically rewired by elevated concentrations of fatty acids and glucocorticoids, augmenting fatty acid metabolism via increased expression of the enzyme CPT1A, and Cpt1a deletion impaired NK cell survival and function in this setting. In parallel, redistribution of NK cells to the BM during fasting required the trafficking mediators S1PR5 and CXCR4. These cells were primed by an increased pool of interleukin (IL)-12-expressing BM myeloid cells, which improved IFN-γ production. Our findings identify a link between dietary restriction and optimized innate immune responses, with the potential to enhance immunotherapy strategies.

Epithelium/imcDC2 axis facilitates the resistance of neoadjuvant anti-PD-1 in human NSCLC

BACKGROUND

Therapeutic resistance is a main obstacle to achieve long-term benefits from immune checkpoint inhibitors. The underlying mechanism of neoadjuvant anti-PD-1 resistance remains unclear.

METHODS

Multi-omics analysis, including mass cytometry, single-cell RNA-seq, bulk RNA-seq, and polychromatic flow cytometry, was conducted using the resected tumor samples in a cohort of non-small cell lung cancer (NSCLC) patients received neoadjuvant anti-PD-1 therapy. Tumor and paired lung samples acquired from treatment-naïve patients were used as a control. In vitro experiments were conducted using primary cells isolated from fresh tissues and lung cancer cell lines. A Lewis-bearing mouse model was used in the in vivo experiment.

RESULTS

The quantity, differentiation status, and clonal expansion of tissue-resident memory CD8+ T cells (CD8+ TRMs) are positively correlated with therapeutic efficacy of neoadjuvant anti-PD-1 therapy in human NSCLC. In contrast, the quantity of immature CD1c+ classical type 2 dendritic cells (imcDC2) and galectin-9+ cancer cells is negatively correlated with therapeutic efficacy. An epithelium/imDC2 suppressive axis that restrains the antitumor response of CD8+ TRMs via galectin-9/TIM-3 was uncovered. The expression level of CD8+ TRMs and galectin-9+ cancer cell-related genes predict the clinical outcome of anti-PD-1 neoadjuvant therapy in human NSCLC patients. Finally, blockade of TIM-3 and PD-1 could improve the survival of tumor-bearing mouse by promoting the antigen presentation of imcDC2 and CD8+ TRMs-mediated tumor-killing.

CONCLUSION

Galectin-9 expressing tumor cells sustained the primary resistance of neoadjuvant anti-PD-1 therapy in NSCLC through galectin-9/TIM-3-mediated suppression of imcDC2 and CD8+ TRMs. Supplement of anti-TIM-3 could break the epithelium/imcDC2/CD8+ TRMs suppressive loop to overcome anti-PD-1 resistance.

TRIAL REGISTRATION NUMBER

NCT03732664.

Anti-HLA serologic response to CD38-targeting desensitization therapy is challenged by peripheral memory B cells in highly sensitized kidney transplant candidates

High human leukocyte antigen (HLA) sensitization limits access to compatible transplantation. New CD38-targeting agents have been shown to reduce anti-HLA antibodies, although with important interpatient variability. Thus, pretreatment identification of responder and nonresponder (NR) patients is needed for treatment decision-making. We analyzed 26 highly sensitized (HS) patients from 2 desensitization trials using anti-CD38 monoclonal antibodies. Hierarchical clustering identified 3 serologic responder groups: high responders, low responders, and NR. Spectral flow cytometry and functional HLA-specific memory B cell (mBC) assessment were first conducted on peripheral blood mononuclear cells and bone marrow samples from 16 patients treated with isatuximab (NCT04294459). Isatuximab effectively depleted bone marrow plasma cells, peripheral CD38-expressing plasmablasts, plasma cells, transitional B cells, and class-switch mBCs, ultimately reducing frequencies of HLA-specific immunoglobulin G (IgG)-producing mBCs. Multidimensional spectral flow cytometry with partial least squares discriminant analysis revealed that pretreatment abundance of specific circulating mBC phenotypes, especially CD38neg class-switch mBCs, accurately distinguished between high serologic responders and low responders or NR (AUC 0.958, 0.860-1.000, P = .009), who also displayed significantly lower frequencies of HLA-specific IgG-producing mBCs (P < .0001). This phenotypical mBC signature predicting response to therapy was validated in an external HS patient cohort (n = 10) receiving daratumumab (NCT04204980). This study identifies critical circulating mBC subset phenotypes that distinguish HS patients with successful serologic responses to CD38-targeting desensitization therapies, potentially guiding treatment decision-making.

Single cell dual-omic atlas of the human developing retina

The development of the retina is under tight temporal and spatial control. To gain insights into the molecular basis of this process, we generate a single-nuclei dual-omic atlas of the human developing retina with approximately 220,000 nuclei from 14 human embryos and fetuses aged between 8 and 23-weeks post-conception with matched macular and peripheral tissues. This atlas captures all major cell classes in the retina, along with a large proportion of progenitors and cell-type-specific precursors. Cell trajectory analysis reveals a transition from continuous progression in early progenitors to a hierarchical development during the later stages of cell type specification. Both known and unrecorded candidate transcription factors, along with gene regulatory networks that drive the transitions of various cell fates, are identified. Comparisons between the macular and peripheral retinae indicate a largely consistent yet distinct developmental pattern. This atlas offers unparalleled resolution into the transcriptional and chromatin accessibility landscapes during development, providing an invaluable resource for deeper insights into retinal development and associated diseases.

scMaui: a widely applicable deep learning framework for single-cell multiomics integration in the presence of batch effects and missing data

The recent advances in high-throughput single-cell sequencing have created an urgent demand for computational models which can address the high complexity of single-cell multiomics data. Meticulous single-cell multiomics integration models are required to avoid biases towards a specific modality and overcome sparsity. Batch effects obfuscating biological signals must also be taken into account. Here, we introduce a new single-cell multiomics integration model, Single-cell Multiomics Autoencoder Integration (scMaui) based on variational product-of-experts autoencoders and adversarial learning. scMaui calculates a joint representation of multiple marginal distributions based on a product-of-experts approach which is especially effective for missing values in the modalities. Furthermore, it overcomes limitations seen in previous VAE-based integration methods with regard to batch effect correction and restricted applicable assays. It handles multiple batch effects independently accepting both discrete and continuous values, as well as provides varied reconstruction loss functions to cover all possible assays and preprocessing pipelines. We demonstrate that scMaui achieves superior performance in many tasks compared to other methods. Further downstream analyses also demonstrate its potential in identifying relations between assays and discovering hidden subpopulations.

Definition of a Multi-Omics Signature for Esophageal Adenocarcinoma Prognosis Prediction

Esophageal cancer is a highly lethal malignancy, representing 5% of all cancer-related deaths. The two main subtypes are esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While most research has focused on ESCC, few studies have analyzed EAC for transcriptional signatures linked to diagnosis or prognosis. In this study, we utilized single-cell RNA sequencing and bulk RNA sequencing to identify specific immune cell types that contribute to anti-tumor responses, as well as differentially expressed genes (DEGs). We have characterized transcriptional signatures, validated against a wide cohort of TCGA patients, that are capable of predicting clinical outcomes and the prognosis of EAC post-surgery with efficacy comparable to the currently accepted prognostic factors. In conclusion, our findings provide insights into the immune landscape and therapeutic targets of EAC, proposing novel immunological biomarkers for predicting prognosis, aiding in patient stratification for post-surgical outcomes, follow-up, and personalized adjuvant therapy decisions.

aKNNO: single-cell and spatial transcriptomics clustering with an optimized adaptive k-nearest neighbor graph

Typical clustering methods for single-cell and spatial transcriptomics struggle to identify rare cell types, while approaches tailored to detect rare cell types gain this ability at the cost of poorer performance for grouping abundant ones. Here, we develop aKNNO to simultaneously identify abundant and rare cell types based on an adaptive k-nearest neighbor graph with optimization. Benchmarking on 38 simulated and 20 single-cell and spatial transcriptomics datasets demonstrates that aKNNO identifies both abundant and rare cell types more accurately than general and specialized methods. Using only gene expression aKNNO maps abundant and rare cells more precisely compared to integrative approaches.

A spatial expression atlas of the adult human proximal small intestine

The mouse small intestine shows profound variability in gene expression along the crypt-villus axis1,2. Whether similar spatial heterogeneity exists in the adult human gut remains unclear. Here we use spatial transcriptomics, spatial proteomics and single-molecule fluorescence in situ hybridization to reconstruct a comprehensive spatial expression atlas of the adult human proximal small intestine. We describe zonated expression and cell type representation for epithelial, mesenchymal and immune cell types. We find that migrating enterocytes switch from lipid droplet assembly and iron uptake at the villus bottom to chylomicron biosynthesis and iron release at the tip. Villus tip cells are pro-immunogenic, recruiting γδ T cells and macrophages to the tip, in contrast to their immunosuppressive roles in mouse. We also show that the human small intestine contains abundant serrated and branched villi that are enriched at the tops of circular folds. Our study presents a detailed resource for understanding the biology of the adult human small intestine.