Single-Cell Transcriptomic Analysis Revealed a Critical Role of SPP1/CD44-Mediated Crosstalk Between Macrophages and Cancer Cells in Glioma

Single-Cell and Bulk RNA Sequencing Reveal Tumor Cell Characteristics and Communication Features of Primary and Lymphatic Metastatic Hypopharyngeal Squamous Cell Carcinoma

BACKGROUND

Lymph node metastasis (LNM) is strongly associated with poor prognosis in hypopharyngeal squamous carcinoma (HPSCC). Identifying key drivers of LNM and potential therapeutic targets in HPSCC is therefore essential for the early detection of high-risk patients and for informing personalized treatment strategies.

METHODS

Single-cell RNA sequencing data were used to characterize malignant epithelial cells (maECs) in HPSCC primary tumors (PT) and LNM, as well as differences in cell-to-cell communication. Concurrently, combined with bulk RNA sequencing data, a ligand receptor pairs (LRs) model was developed to predict the prognosis of HPSCC patients.

RESULTS

PT and LNM maECs have different gene expression characteristics, with genes involved in interferon signaling and TGF-β response pathways enriched in LNM maECs, suggesting potential immunosuppressive reprogramming. Cell communication analysis revealed distinct interactions and signaling features in PT and LNM microenvironments. Subsequently, a 4-LRs model was constructed to stratify HPSCC patients into low-or high-risk groups, with the high-risk group demonstrating significantly worse overall survival (OS) outcomes compared with the low-risk group in the training (p < 0.0001), testing (p = 0.0021), and entire (p < 0.0001) cohorts. Receiver operating characteristic curves showed that this risk model can effectively predict the 1-, 3-, and 5-year OS of HPSCC patients. Notably, the risk score effectively discriminated LNM status (area under the curve [AUC] = 0.927) among HPSCC patients, highlighting its potential as a HPSCC metastasis prediction tool.

CONCLUSIONS

These results provide biomarkers of LNM maECs as well as potential mechanisms of HPSCC metastasis, which may help with the precision treatment, diagnosis, and prognosis of HPSCC.

Investigative needle core biopsies support multimodal deep-data generation in glioblastoma

Glioblastoma (GBM) is an aggressive primary brain cancer with few effective therapies. Stereotactic needle biopsies are routinely used for diagnosis; however, the feasibility and utility of investigative biopsies to monitor treatment response remains ill-defined. Here, we demonstrate the depth of data generation possible from routine stereotactic needle core biopsies and perform highly resolved multi-omics analyses, including single-cell RNA sequencing, spatial transcriptomics, metabolomics, proteomics, phosphoproteomics, T-cell clonotype analysis, and MHC Class I immunopeptidomics on standard biopsy tissue obtained intra-operatively. We also examine biopsies taken from different locations and provide a framework for measuring spatial and genomic heterogeneity. Finally, we investigate the utility of stereotactic biopsies as a method for generating patient-derived xenograft (PDX) models. Multimodal dataset integration highlights spatially mapped immune cell-associated metabolic pathways and validates inferred cell-cell ligand-receptor interactions. In conclusion, investigative biopsies provide data-rich insight into disease processes and may be useful in evaluating treatment responses.

Microglia and brain macrophages are differentially associated with tumor necrosis in glioblastoma: A link to tumor progression

Background

Microglia and brain macrophages contribute significantly to the tumor microenvironment in highly malignant glioblastoma where they are considered important drivers of tumor progression. A better understanding of the role of the brain macrophages present in glioblastoma appears crucial for improving therapeutic outcomes, especially in the context of novel immunotherapeutic approaches.

Methods

We investigated the regulation of two well-established markers for microglia and brain macrophages, IBA1 and CD163, in relation to glioblastoma tumor necrosis using immunohistochemistry and modality fusion heatmaps of whole slide images obtained from adjacent tissue sections.

Results

IBA1 and CD163 showed remarkable differences in relation to glioblastoma tumor necrosis. Generally, IBA1 immunoreactive cells were far less common in necrotic tissue areas than CD163-expressing cells. We also found extensive and frequently diffuse extracellular CD163 deposition, especially in hypocellular necrobiotic tumor regions where IBA1 was typically absent.

Conclusions

Resident microglia seem more likely to be important for the diffuse infiltration of glioma cells in hypercellular tissue areas, whereas myeloid macrophages may be the main macrophage population in the wake of tumor necrosis. Since the necrotic niche with its interactions between microglia, brain macrophages, and glioblastoma/glioma stem cells is increasingly recognised as an important factor in tumor progression, further detailed studies of the macrophage populations in glioblastoma are warranted.

Single-cell RNA sequencing identifies molecular biomarkers predicting late progression to CDK4/6 inhibition in patients with HR+/HER2- metastatic breast cancer

BACKGROUND

Cyclin-dependent kinase 4/6 inhibitors (CDK4/6is) in combination with endocrine therapy are the standard treatment for patients with hormone receptor-positive, HER2-negative metastatic breast cancer (mBC). Despite the efficacy of CDK4/6is, intrinsic resistance occurs in approximately one-third of patients, highlighting the need for reliable predictive biomarkers.

METHODS

Single-cell RNA sequencing analyzed metastatic tumors from HR+/HER2- mBC patients pre-CDK4/6i treatment at baseline (BL) and/or at disease progression. BL samples were from CDK4/6i responders (median progression-free survival [mPFS] = 25.5 months), while progressors were categorized as early-progressors (EP, mPFS = 3 months) and late-progressors (LP, mPFS = 11 months). Metastatic sites included liver, pleural effusions, ascites, and bone. InferCNV distinguished tumor cells, and functional analysis utilized the Molecular Signatures Database.

RESULTS

LP tumors displayed enhanced Myc, EMT, TNF-α, and inflammatory pathways compared to those EP tumors. Samples from BL and LP responders showed increased tumor-infiltrating CD8+ T cells and natural killer (NK) cells compared to EP non-responders. Notably, despite a high frequency of CD8+ T cells in responding tumors, a functional analysis revealed significant upregulation of genes associated with stress and apoptosis in proliferative CD4+ and CD8+ T cells in BL tumors compared to in EP and LP tumors. These genes, including HSP90 and HSPA8, are linked to resistance to PD1/PD-L1 immune checkpoint inhibitors. A ligand-receptor analysis showed enhanced interactions associated with inhibitory T-cell proliferation (SPP1-CD44) and suppression of immune activity (MDK-NCL) in LP tumors. Longitudinal biopsies consistently revealed dynamic NK cell expansion and enhanced cytotoxic T cell activity, alongside upregulation of immune activity inhibition, in LP tumors compared to in BL tumors. Notably, the predictive biomarker panel from BL tumor cells was validated in 2 independent cohorts, where it consistently predicted a significant improvement in mPFS duration in signature-high versus -low groups.

CONCLUSION

This study underscores the significance of molecular biomarkers in predicting clinical outcomes to CDK4/6i. Tumor-infiltration CD8+ T and NK cells may also serve as baseline predictors. These insights pave the way for optimizing therapeutic strategies based on microenvironment-specific changes, providing a personalized and effective approach for managing HR+/HER2- mBC and improving patient outcomes.

STAT3-controlled CHI3L1/SPP1 positive feedback loop demonstrates the spatial heterogeneity and immune characteristics of glioblastoma

Proneural-mesenchymal transition (PMT) is a phenotypic alteration and contributes to the malignant progression of glioblastoma (GBM). Macrophages, as a main infiltrating component of the tumor immune microenvironment (TIM), control the biological processes of PMT; however, the mechanisms driving this process remain largely unknown. Here, the overall landscape of tumor and nontumor cells was described by scMulti-omics technology. Then, we demonstrated that chitinase-3-like protein 1 (CHI3L1) played a critical role in maintaining mesenchymal (MES) status and reprogramming macrophage phenotype using C57BL/6 and NSG mice models derived from PN20 cells. Mechanistically, osteopontin (OPN)/ITGB1 maintained the activation of nuclear factor κB (NF-κB) and signal transducer and activator of transcription 3 (STAT3) pathways by establishing a positive feedback loop with the CHI3L1-STAT3 axis, resulting in PMT. CHI3L1 enhanced the phosphorylation, nuclear localization, and transcriptional activity of STAT3 via directly binding its coiled-coil domain (CCD). Importantly, we screened and validated that hygromycin B (HB), an inhibitor of the STAT3-CCD domain, disrupted the CHI3L1-STAT3 interaction, thereby reducing the tumor burden in vitro and in vivo.

Gliomagenesis mimics an injury response orchestrated by neural crest-like cells

Glioblastoma is an incurable brain malignancy. By the time of clinical diagnosis, these tumours exhibit a degree of genetic and cellular heterogeneity that provides few clues to the mechanisms that initiate and drive gliomagenesis1,2. Here, to explore the early steps in gliomagenesis, we utilized conditional gene deletion and lineage tracing in tumour mouse models, coupled with serial magnetic resonance imaging, to initiate and then closely track tumour formation. We isolated labelled and unlabelled cells at multiple stages-before the first visible abnormality, at the time of the first visible lesion, and then through the stages of tumour growth-and subjected cells of each stage to single-cell profiling. We identify a malignant cell state with a neural crest-like gene expression signature that is highly abundant in the early stages, but relatively diminished in the late stage of tumour growth. Genomic analysis based on the presence of copy number alterations suggests that these neural crest-like states exist as part of a heterogeneous clonal hierarchy that evolves with tumour growth. By exploring the injury response in wounded normal mouse brains, we identify cells with a similar signature that emerge following injury and then disappear over time, suggesting that activation of an injury response program occurs during tumorigenesis. Indeed, our experiments reveal a non-malignant injury-like microenvironment that is initiated in the brain following oncogene activation in cerebral precursor cells. Collectively, our findings provide insight into the early stages of glioblastoma, identifying a unique cell state and an injury response program tied to early tumour formation. These findings have implications for glioblastoma therapies and raise new possibilities for early diagnosis and prevention of disease.

Gemistocytic tumor cells programmed for glial scarring characterize T cell confinement in IDH-mutant astrocytoma

Isocitrate dehydrogenase 1/2 mutant (IDHmt) astrocytoma is considered a T cell-deprived tumor, yet little is known regarding the phenotypes underlying T cell exclusion. Using bulk, single nucleus and spatial RNA and protein profiling, we demonstrate that a distinct spatial organization underlies T cell confinement to the perivascular space (T cell cuff) in IDHmt astrocytoma. T cell cuffs are uniquely characterized by a high abundance of gemistocytic tumor cells (GTC) in the surrounding stroma. Integrative analysis shows that GTC-high tumors are enriched for lymphocytes and tumor associated macrophages (TAM) and express immune cell migration and activation programs. Specifically, GTCs constitute a distinct sub-cluster of the astrocyte-like tumor cell state that co-localizes with immune reactive TAMs. Neighboring GTCs and TAMs express receptor-ligand pairs characteristic of reactive astrogliosis and glial scarring, such as SPP1/CD44 and IL-1β/IL1R1. Collectively, we reveal that T cell confinement in IDHmt astrocytomas associates with GTC-TAM networks that mimic glial scarring mechanisms.

Single-nucleus and spatial landscape of the sub-ventricular zone in human glioblastoma

The sub-ventricular zone (SVZ) is the most well-characterized neurogenic area in the mammalian brain. We previously showed that in 65% of patients with glioblastoma (GBM), the SVZ is a reservoir of cancer stem-like cells that contribute to treatment resistance and the emergence of recurrence. Here, we build a single-nucleus RNA-sequencing-based microenvironment landscape of the tumor mass and the SVZ of 15 patients and two histologically normal SVZ samples as controls. We identify a ZEB1-centered mesenchymal signature in the tumor cells of the SVZ. Moreover, the SVZ microenvironment is characterized by tumor-supportive microglia, which spatially coexist and establish crosstalks with tumor cells. Last, differential gene expression analyses, predictions of ligand-receptor and incoming/outgoing interactions, and functional assays reveal that the interleukin (IL)-1β/IL-1RAcP and Wnt-5a/Frizzled-3 pathways represent potential therapeutic targets in the SVZ. Our data provide insights into the biology of the SVZ in patients with GBM and identify potential targets of this microenvironment.

Transformation of brain myeloid cell populations by SIV in rhesus macaques revealed by multiomics

The primary immune constituents in the brain, microglia and macrophages, are the target for HIV in people and simian immunodeficiency virus (SIV) in nonhuman primates. This infection can lead to neurological dysfunction, known as HIV-associated neurocognitive disorder (HAND). Given the gaps in our knowledge on how these cells respond in vivo to CNS infection, we perform single-cell multiomic sequencing, including gene expression and ATAC-seq, on myeloid cells from the brains of rhesus macaques with SIV-induced encephalitis (SIVE) as well as uninfected controls. We find that SIVE significantly changes the myeloid cell populations. In SIVE, microglia-like cells express high levels of chemoattractants capable of recruiting highly activated CAM-like cells to the site of infection/inflammation. A unique population of microglia-like cells is found in which the chromatin accessibility of genes diverges from their RNA expression. Additionally, we observe a dramatic shift of upstream gene regulators and their targets in brain myeloid cells during SIVE. This study further uncovers the transcriptome, gene regulatory events, and potential roles of different brain myeloid phenotypes in SIVE. This might deepen the understanding of SIVE/HIVE and enlighten the therapeutic development.

Exploring the role of Disulfidptosis in glioma progression: insights into tumor heterogeneity and therapeutic potential through single-cell RNA sequencing

BACKGROUND

Gliomas, particularly glioblastoma (GBM), are the most common and aggressive primary brain tumors in adults, characterized by high malignancy and frequent recurrence. Despite standard treatments, including surgery, radiotherapy, and chemotherapy, the prognosis for GBM remains poor, with a median survival of less than 15 months and a five-year survival rate below 10%. Tumor heterogeneity and resistance to treatment create significant challenges in controlling glioma progression. Therefore, there is an urgent need for new therapeutic targets and strategies.

OBJECTIVE

This study investigates the role of Disulfidptosis, a recently discovered form of programmed cell death, in gliomas. Unlike apoptosis and necrosis, Disulfidptosis is driven by the abnormal accumulation of intracellular disulfide bonds, leading to protein misfolding and cytoskeletal collapse, particularly in cancer cells with metabolic dysregulation. We aim to explore how glioma cells respond to Disulfidptosis and identify potential therapeutic targets by analyzing the heterogeneity of gliomas at the single-cell level using single-cell RNA sequencing (scRNA-seq).

METHODS

scRNA-seq data from glioma patients were analyzed to uncover differences in ferroptosis-related pathways, including iron metabolism and lipid peroxidation. Cellular subpopulations within gliomas were profiled to assess their sensitivity to Disulfidptosis and the underlying mechanisms. Survival analysis was conducted to evaluate the clinical relevance of Disulfidptosis-related gene expression.

RESULTS

Multiple cell subpopulations within gliomas exhibit varying sensitivities to Disulfidptosis, influenced by their metabolic properties. Dysregulated iron metabolism and antioxidant mechanisms were identified as key factors impacting Disulfidptosis sensitivity. Glioma microenvironment signaling pathways also play a role in regulating Disulfidptosis. These findings suggest that activating Disulfidptosis pathways may provide novel therapeutic strategies to overcome treatment resistance in gliomas.

CONCLUSION

This study offers new insights into the role of Disulfidptosis in glioma progression and highlights its potential as a therapeutic target. By leveraging single-cell sequencing data, the research uncovers tumor heterogeneity and identifies specific cell populations resistant to Disulfidptosis. These findings may pave the way for personalized treatment strategies to improve survival outcomes in glioma patients.

Identification of SPP1 + macrophages as an immune suppressor in hepatocellular carcinoma using single-cell and bulk transcriptomics

Introduction

Macrophages and T cells play crucial roles in liver physiology, but their functional diversity in hepatocellular carcinoma (HCC) remains largely unknown.

Methods

Two bulk RNA-sequencing (RNA-seq) cohorts for HCC were analyzed using gene co-expression network analysis. Key gene modules and networks were mapped to single-cell RNA-sequencing (scRNA-seq) data of HCC. Cell type fraction of bulk RNA-seq data was estimated by deconvolution approach using single-cell RNA-sequencing data as a reference. Survival analysis was carried out to estimate the prognosis of different immune cell types in bulk RNA-seq cohorts. Cell-cell interaction analysis was performed to identify potential links between immune cell types in HCC.

Results

In this study, we analyzed RNA-seq data from two large-scale HCC cohorts, revealing a major and consensus gene co-expression cluster with significant implications for immunosuppression. Notably, these genes exhibited higher enrichment in liver macrophages than T cells, as confirmed by scRNA-seq data from HCC patients. Integrative analysis of bulk and single-cell RNA-seq data pinpointed SPP1 + macrophages as an unfavorable cell type, while VCAN + macrophages, C1QA + macrophages, and CD8 + T cells were associated with a more favorable prognosis for HCC patients. Subsequent scRNA-seq investigations and in vitro experiments elucidated that SPP1, predominantly secreted by SPP1 + macrophages, inhibits CD8 + T cell proliferation. Finally, targeting SPP1 in tumor-associated macrophages through inhibition led to a shift towards a favorable phenotype.

Discussion

This study underpins the potential of SPP1 as a translational target in immunotherapy for HCC.

Cellular MSI-H score: a robust predictive biomarker for immunotherapy response and survival in gastrointestinal cancer

Microsatellite instability-high (MSI-H) is a critical biomarker for immunotherapy, yet primary resistance remains a significant challenge. Current MSI-H detection methods evaluate the proportion of MSI-H loci, termed molecular MSI-H score, which can be affected by intratumoral heterogeneity (ITH). To address this limitation, we propose evaluating MSI-H at the cellular level to improve the prediction of immunotherapy outcomes. Using bulk tissue (TCGA-CRC) and cell line (CCLE-CRC) datasets, we identified genes highly expressed in MSI-H and MSS samples. These signatures were applied to a single-cell RNA sequencing (scCRC) dataset for enrichment analysis, enabling classification of tumor cells into MSI-H, MSS, and microsatellite dual (MSD) clusters using a Gaussian finite mixture model. Validation showed that MSI-H and MSS enrichment scores were higher in mismatch repair-deficient (MMRd) and mismatch repair-proficient (MMRp) patients, respectively. Functional enrichment analysis revealed that MSI-H cells were associated with pathways such as carboxylic acid catabolism, inflammatory responses, and IL-6/JAK2/STAT3 signaling. We developed a cellular MSI-H signature using genes specifically expressed in the MSI-H cell cluster and transformed the scCRC dataset into a cell-type-specific pseudobulk expression matrix. Using this matrix as a reference, we performed reference-based deconvolution on TCGA-CRC data. We defined the deconvolution score of MSI-H cell as cellular MSI-H score. This score strongly correlated with the molecular MSI-H score (R = 0.55, P < 0.001) and showed modest correlations with macrophage (MoMac, R = 0.14) and CD8+ T-cell (R = 0.11). To investigate its potential for clinical application, we applied the cellular MSI-H signature to the BJ-cohort, comprising 97 immunotherapy-treated gastrointestinal patients sequenced with a 395-gene panel. The cellular MSI-H score was significantly higher in responders (P = 0.002), positively correlated with tumor reduction percentage (R = 0.29, P = 0.006), and associated with improved progression-free survival (PFS) (HR: 0.00, 95% CI: 0.00-0.31, P = 0.021). In summary, the cellular MSI-H score reflects the MSI-H cell level within a tumor and demonstrates superior accuracy compared to molecular MSI-H status in predicting immunotherapy response and PFS. This underscores its potential as a more robust biomarker for guiding immunotherapy decisions.

Exploring the prognostic value of BRMS1 + microglia based on single-cell anoikis regulator patterns in the immunologic microenvironment of GBM

BACKGROUND

Anoikis is a specialized form of programmed cell death induced by the loss of cell adhesion to the extracellular matrix (ECM). Acquisition of anoikis resistance is a significant marker for cancer cell invasion, metastasis, therapy resistance, and recurrence. Although current research has identified multiple factors that regulate anoikis resistance, the pathological mechanisms of anoikis-mediated tumor microenvironment (TME) in glioblastoma (GBM) remain largely unexplored.

METHODS

Utilizing single-cell RNA sequencing (scRNA-seq) data and employing non-negative matrix factorization (NMF), we identified and characterized TME cell clusters with distinct anoikis-associated gene signatures. Prognostic and therapeutic response analyses were conducted using TCGA and CGGA datasets to assess the clinical significance of different TME cell clusters. The spatial relationship between BRMS1 + microglia and tumor cells was inferred from spatial transcriptome RNA sequencing (stRNA-seq) data. To simulate the tumor immune microenvironment, co-culture experiments were performed with microglia (HMC3) and GBM cells (U118/U251), and microglia were transfected with a BRMS1 overexpression lentivirus. Western blot or ELISA were used to detect BRMS1, M2 macrophage-specific markers, PI3K/AKT signaling proteins, and apoptosis-related proteins. The proliferation and apoptosis capabilities of tumor cells were evaluated using CCK-8, colony formation, and apoptosis assays, while the invasive and migratory abilities of tumor cells were assessed using Transwell assays.

RESULTS

NMF-based analysis successfully identified CD8 + T cell and microglia cell clusters with distinct gene signature characteristics. Trajectory analysis, cell communication, and gene regulatory network analyses collectively indicated that anoikis-mediated TME cell clusters can influence tumor cell development through various mechanisms. Notably, BRMS1 + AP-Mic exhibited an M2 macrophage phenotype and had significant cell communication with malignant cells. Moreover, high expression of BRMS1 + AP-Mic in TCGA and CGGA datasets was associated with poorer survival outcomes, indicating its detrimental impact on immunotherapy. Upregulation of BRMS1 in microglia may lead to M2 macrophage polarization, activate the PI3K/AKT signaling pathway through SPP1/CD44-mediated cell interactions, inhibit tumor cell apoptosis, and promote tumor proliferation and invasion.

CONCLUSION

This pioneering study used NMF-based analysis to reveal the important predictive value of anoikis-regulated TME in GBM for prognosis and immunotherapeutic response. BRMS1 + microglial cells provide a new perspective for a deeper understanding of the immunosuppressive microenvironment of GBM and could serve as a potential therapeutic target in the future.

SPP1+ macrophages in HR+ breast cancer are associated with tumor-infiltrating lymphocytes

Breast cancer categorized into hormone receptor-positive (HR+), HER2-positive (HER2+), and triple-negative (TNBC) subtypes, exhibits varied outcomes based on the number of tumor-infiltrating lymphocytes (TILs). To explore the divergent roles of TIL levels across different subtypes, we employed single-cell RNA sequencing on 31 patients with breast cancer. HR+ breast cancer with high TIL levels (TIL-high) revealed increased SPP1+ macrophages, increased SPP1 expression in other monocytes/macrophages (mono/macro) subgroups, and enriched pathways associated with extracellular matrix (ECM) remodeling in mono/macro. Moreover, cell-cell interaction analyses revealed enhanced SPP1, MIF, and FN1 signaling in the interaction between SPP1+ macrophages and T-cells in TIL-high HR+ breast cancer. Spatial transcriptomics data highlighted the close proximity of SPP1+ macrophages, CD8+ T-cells, and CD4+ T-cells in TIL-high HR+ breast cancer. Our findings unveil the novel influence of SPP1+ macrophages on T-cells in TIL-high HR+ breast cancer, potentially explaining the poor prognosis and offering insights for targeted interventions.

Elevated α-1,2-mannosidase MAN1C1 in glioma stem cells and its implications for immunological changes and prognosis in glioma patients

Glioblastoma multiforme (GBM) is the most aggressive type of primary brain tumor, and the presence of glioma stem cells (GSCs) has been linked to its resistance to treatments and recurrence. Additionally, aberrant glycosylation has been implicated in the aggressiveness of cancers. However, the influence and underlying mechanism of N-glycosylation on the GSC phenotype and GBM malignancy remain elusive. Here, we performed an in-silico analysis approach on publicly available datasets to examine the function of N-glycosylation-related genes in GSCs and gliomas, accompanied by a qRT-PCR validation experiment. We found that high α-1,2-mannosidase MAN1C1 is associated with immunological functions and worse survival of glioma patients. Differential gene expression analysis and qRT-PCR validation revealed that MAN1C1 is highly expressed in GSCs. Furthermore, higher MAN1C1 expression predicts worse outcomes in glioma patients. Also, MAN1C1 expression is increased in the perinecrotic region of GBM and is associated with immunological and inflammatory functions, a hallmark of the GBM mesenchymal subtype. Further analysis confirmed that MAN1C1 expression is closely associated with infiltrating immune cells and disrupted immune response in the GBM microenvironment. These suggest that MAN1C1 is a potential biomarker for gliomas and may be important as an immunotherapeutic target for GBM.

Transformation of brain myeloid cell populations by SIV in rhesus macaques revealed by multiomics

The primary immune constituents in the brain, microglia and macrophages, are the target for HIV in people and simian immunodeficiency virus (SIV) in nonhuman primates. This infection can lead to neurological dysfunction, known as HIV-associated neurocognitive disorder (HAND). Given the gaps in our knowledge on how these cells respond in vivo to CNS infection, we performed single-cell multiomic sequencing, including gene expression and ATAC-seq, on myeloid cells from the brains of rhesus macaques with SIV-induced encephalitis (SIVE) as well as uninfected controls. We found that the myeloid cell populations were significantly changed by SIVE. In SIVE microglia-like cells express high levels of chemoattractants capable of recruiting highly activated CAM-like cells to the site of infection/inflammation. A unique population of microglia-like cells was found in which the chromatin accessibility of genes diverged from their RNA expression. Additionally, we observed a dramatic shift of upstream gene regulators and their targets in brain myeloid cells during SIVE. In summary, this study further uncovers the transcriptome, gene regulatory events and potential roles of different brain myeloid phenotypes in SIVE.

A Single-Cell Transcriptome Profiling of Triptolide-Induced Nephrotoxicity in Mice

Triptolide (TP), an active component isolated from the traditional Chinese herb Tripterygium wilfordii Hook F (TWHF), shows great promise for treating inflammation-related diseases. However, its potential nephrotoxic effects remain concerning. The mechanism underlying TP-induced nephrotoxicity is inadequately elucidated, particularly at single-cell resolution. Hence, single-cell RNA sequencing (scRNA-seq) of kidney tissues from control and TP-treated mice is performed to generate a thorough description of the renal cell atlas upon TP treatment. Heterogeneous responses of nephron epithelial cells are observed after TP exposure, attributing differential susceptibility of cell subtypes to excessive reactive oxygen species and increased inflammatory responses. Moreover, TP disrupts vascular function by activating endothelial cell immunity and damaging fibroblasts. Severe immune cell damage and the activation of pro-inflammatory Macro_C1 cells are also observed with TP treatment. Additionally, ligand-receptor crosstalk analysis reveals that the SPP1 (osteopontin) signaling pathway targeting Macro_C1 cells is triggered by TP treatment, which may promote the infiltration of Macro_C1 cells to exacerbate renal toxicity. Overall, this study provides comprehensive information on the transcriptomic profiles and cellular composition of TP-associated nephrotoxicity at single-cell resolution, which can strengthen the understanding of the pathogenesis of TP-induced nephrotoxicity and provide valuable clues for the discovery of new therapeutic targets to ameliorate TP-associated nephrotoxicity.

Single-cell and spatial transcriptome assays reveal heterogeneity in gliomas through stress responses and pathway alterations

Background

Glioma is a highly heterogeneous malignancy of the central nervous system. This heterogeneity is driven by various molecular processes, including neoplastic transformation, cell cycle dysregulation, and angiogenesis. Among these biomolecular events, inflammation and stress pathways in the development and driving factors of glioma heterogeneity have been reported. However, the mechanisms of glioma heterogeneity under stress response remain unclear, especially from a spatial aspect.

Methods

This study employed single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) to explore the impact of oxidative stress response genes in oligodendrocyte precursor cells (OPCs). Our analysis identified distinct pathways activated by oxidative stress in two different types of gliomas: high- and low- grade (HG and LG) gliomas.

Results

In HG gliomas, oxidative stress induced a metabolic shift from oxidative phosphorylation to glycolysis, promoting cell survival by preventing apoptosis. This metabolic reprogramming was accompanied by epithelial-to-mesenchymal transition (EMT) and an upregulation of stress response genes. Furthermore, SCENIC (Single-Cell rEgulatory Network Inference and Clustering) analysis revealed that oxidative stress activated the AP1 transcription factor in HG gliomas, thereby enhancing tumor cell survival and proliferation.

Conclusion

Our findings provide a novel perspective on the mechanisms of oxidative stress responses across various grades of gliomas. This insight enhances our comprehension of the evolutionary processes and heterogeneity within gliomas, potentially guiding future research and therapeutic strategies.

Monitoring melanoma patients on treatment reveals a distinct macrophage population driving targeted therapy resistance

Resistance to targeted therapy remains a major clinical challenge in melanoma. To uncover resistance mechanisms, we perform single-cell RNA sequencing on fine-needle aspirates from resistant and responding tumors of patients undergoing BRAFi/MEKi treatment. Among the genes most prominently expressed in resistant tumors is POSTN, predicted to signal to a macrophage population associated with targeted therapy resistance (TTR). Accordingly, tumors from patients with fast disease progression after therapy exhibit high POSTN expression levels and high numbers of TTR macrophages. POSTN polarizes human macrophages toward a TTR phenotype and promotes resistance to targeted therapy in a melanoma mouse model, which is associated with a phenotype change in intratumoral macrophages. Finally, polarized TTR macrophages directly protect human melanoma cells from MEKi-induced killing via CD44 receptor expression on melanoma cells. Thus, interfering with the protective activity of TTR macrophages may offer a strategy to overcome resistance to targeted therapy in melanoma.

Medulloblastoma Spatial Transcriptomics Reveals Tumor Microenvironment Heterogeneity with High-Density Progenitor Cell Regions Correlating with High-Risk Disease

The tumor microenvironment (TME) of medulloblastoma (MB) influences progression and therapy response, presenting a promising target for therapeutic advances. Prior single-cell analyses have characterized the cellular components of the TME but lack spatial context. To address this, we performed spatial transcriptomic sequencing on sixteen pediatric MB samples obtained at diagnosis, including two matched diagnosis-relapse pairs. Our analyses revealed inter- and intra-tumoral heterogeneity within the TME, comprised of tumor-associated astrocytes (TAAs), macrophages (TAMs), stromal components, and distinct subpopulations of MB cells at different stages of neuronal differentiation and cell cycle progression. We identified dense regions of quiescent progenitor-like MB cells enriched in patients with high-risk (HR) features and an increase in TAAs, TAMs, and dysregulated vascular endothelium following relapse. Our study presents novel insights into the spatial architecture and cellular landscape of the medulloblastoma TME, highlighting spatial patterns linked to HR features and relapse, which may serve as potential therapeutic targets.

Multi-omics analysis of macrophage-associated receptor and ligand reveals a strong prognostic signature and subtypes in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a significant contributor to morbidity and mortality worldwide. The interaction between receptors and ligands is the primary mode of intercellular signaling and plays a vital role in the progression of HCC. This study aimed to identify the macrophage-related receptor ligand marker genes associated with HCC and further explored the molecular immune mechanisms attributed to altered biomarkers. Single-cell RNA sequencing data containing primary and recurrent samples were downloaded from the China National GeneBank. Cell types were first identified to explore differences between immune cells from different sample sources. CellChat analysis was used to infer and analyze intercellular communication networks quantitatively. Three molecular subtypes were constructed based on the screened twenty macrophage-associated receptor ligand genes. Bulk RNA-Seq data were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. After the screening, the minor absolute shrinkage and selection operator (LASSO) regression model was employed to identify key markers. After collecting peripheral blood and clinical information from patients, an enzyme-linked immunosorbent assay (ELISA) was used to detect the correlation between key markers and IL-10, one of the macrophage markers. After developing a new HCC risk adjustment model and conducting analysis, it was found that there were significant differences in immune status and gene mutations between the high-risk and low-risk groups of patients based on macrophage-associated receptor and ligand genes. This study identified SPP1, ANGPT2, and NCL as key biological targets for HCC. The drug-gene interaction network analysis identified wortmannin, ribavirin, and tarnafloxin as potential therapeutic drugs for the three key markers. In a clinical cohort study, patients with immune checkpoint inhibitor (ICI) resistance had significantly higher expression levels of OPN, ANGPT2, NCL, and IL-10 than patients with ICI-responsiveness. These three key markers were positively correlated with the expression level of IL-10. The signature based on macrophage-associated receptor and ligand genes can accurately predict the prognosis of patients with HCC and the sensitivity to immunotherapy. These results may help guide the development of targeted prevention and personalized treatment of HCC.

Single-nucleus and spatial landscape of the sub-ventricular zone in human glioblastoma

The sub-ventricular zone (SVZ) is the most well-characterized neurogenic area in the mammalian brain. We previously showed that in 65% of patients with glioblastoma (GBM), the SVZ is a reservoir of cancer stem-like cells that contribute to treatment resistance and emergence of recurrence. Here, we built a single-nucleus RNA-sequencing-based microenvironment landscape of the tumor mass (T_Mass) and the SVZ (T_SVZ) of 15 GBM patients and 2 histologically normal SVZ (N_SVZ) samples as controls. We identified a mesenchymal signature in the T_SVZ of GBM patients: tumor cells from the T_SVZ relied on the ZEB1 regulatory network, whereas tumor cells in the T_Mass relied on the TEAD1 regulatory network. Moreover, the T_SVZ microenvironment was predominantly characterized by tumor-supportive microglia, which spatially co-exist and establish heterotypic interactions with tumor cells. Lastly, differential gene expression analyses, predictions of ligand-receptor and incoming/outgoing interactions, and functional assays revealed that the IL-1β/IL-1RAcP and Wnt-5a/Frizzled-3 pathways are therapeutic targets in the T_SVZ microenvironment. Our data provide insights into the biology of the SVZ in GBM patients and identify specific targets of this microenvironment.

Immune cell infiltration and drug response in glioblastoma multiforme: insights from oxidative stress-related genes

BACKGROUND

GBM, also known as glioblastoma multiforme, is the most prevalent and lethal type of brain cancer. The cell proliferation, invasion, angiogenesis, and treatment of gliomas are significantly influenced by oxidative stress. Nevertheless, the connection between ORGs and GBM remains poorly comprehended. The objective of this research is to investigate the predictive significance of ORGs in GBM and their potential as targets for therapy.

METHODS

We identified differentially expressed genes in glioma and ORGs from public databases. A risk model was established using LASSO regression and Cox analysis, and its performance was evaluated with ROC curves. We then performed consistent cluster analysis on the model, examining its correlation with immunity and drug response. Additionally, PCR, WB and IHC were employed to validate key genes within the prognostic model.

RESULTS

9 ORGs (H6PD, BMP2, SPP1, HADHA, SLC25A20, TXNIP, ACTA1, CCND1, EEF1A1) were selected via differential expression analysis, LASSO and Cox analysis, and incorporated into the risk model with high predictive accuracy. Enrichment analyses using GSVA and GSEA focused predominantly on malignancy-associated pathways. Subtype C of GBM had the best prognosis with the lowest risk score. Furthermore, the model exhibited a strong correlation with the infiltration of immune cells and had the capability to pinpoint potential targeted therapeutic medications for GBM. Ultimately, we selected HADHA for in vitro validation. The findings indicated that GBM exhibits a significant upregulation of HADHA. Knockdown of HADHA inhibited glioma cell proliferation and diminished their migration and invasion capacities and influenced the tumor growth in vivo.

CONCLUSION

The risk model, built upon 9 ORGs and the identification of GBM subtypes, suggests that ORGs have a broad application prospect in the clinical immunotherapy and targeted drug treatment of GBM. HADHA significantly influences the development of gliomas, both in vivo and in vitro.

Pan-cancer and single-cell analyses identify CD44 as an immunotherapy response predictor and regulating macrophage polarization and tumor progression in colorectal cancer

Introduction

Cluster of differentiation (CD) 44 is a non-kinase cell surface transmembrane glycoprotein critical for tumor maintenance and progression.

Methods

We conducted a systematic analysis of the expression profile and genomic alteration profile of CD44 in 33 types of cancer. The immune characteristics of CD44 were comprehensively explored by TIMER2.0 and CIBERSORT. In addition, the CD44 transcriptional landscape was examined at the single-cell level. Then, Pseudotime trajectory analysis of CD44 gene expression was performed using Monocle 2, and CellChat was utilized to compare the crosstalk differences between CD44+monocytes and CD44- monocytes. Tumor immune dysfunction and exclusion (TIDE) was used to evaluate the predictive ability of CD44 for immune checkpoint blockade (ICB) responses. The effects of CD44 on colorectal cancer (CRC) and macrophage polarization were investigated by knocking down the expression of CD44 in HCT-116 cell and macrophages in vitro.

Results

The expression of CD44 elevated in most cancers, predicting unfavorable prognosis. In addditon, CD44 was correlation with immune cell infiltration and key immune regulators. CD44+ monocytes had a higher information flow intensity than CD44- monocytes. CD44 had good predictive ability for immune checkpoint blockade responses. Knockdown of CD44 inhibited the proliferation, migration, and invasion of HCT-116 cell in vitro. Knockdown of CD44 inhibited M2 macrophage polarization.

Discussion

These findings suggest that CD44 is involved in regulating tumor development, macrophage polarization, and has certain predictive value for patient clinical prognosis and response to immunotherapy.

Single-cell and spatial analyses revealed the co-location of cancer stem cells and SPP1+ macrophage in hypoxic region that determines the poor prognosis in hepatocellular carcinoma

In hepatocellular carcinoma (HCC), classical cancer stem cells (CSC) markers were shared by normal stem cells, targeting which may hinder hepatic regeneration and cause liver failure. Additionally, the spatial structure of CSC still remained elusive. To address these limitations, we undertook a comprehensive study combining single-cell data (56,022 cells from 20 samples) and spatial data (38,191 spots from eight samples) to obtain CSC signature and uncover its spatial structure. Utilizing the CytoTRACE algorithm, we discretely identified CSC, which displayed upregulated proliferation pathways regulated by HIF1A. A CSC signature of 107 genes was then developed using Weighted Gene Co-expression Network Analysis (WGCNA). Notably, HCC patients with high CSC levels exhibited an accumulation of SPP1+ macrophages (Macro_SPP1) expressing metalloproteinases (MMP9, MMP12, and MMP7) regulated by HIF1A, suggesting a hypoxic tumor region connecting Macro_SPP1 and CSC. Both CSC and Macro_SPP1 correlated with worse prognosis and undesirable immunotherapy response. Spatial analysis revealed the co-location of CSC and Macro_SPP1, with CD8 T cells excluded from the tumor region. The co-location area and non-tumor area of boundary exhibited a high level of hypoxia, with the HAVRC2 checkpoint highly expressed. Within the co-location area, the SPP1 signaling pathway was most active in cell-cell communication, with SPP1-CD44 and SPP1-ITGA/ITGB identified as the main ligand-receptor pairs. This study successfully constructed a CSC signature and demonstrated the co-location of CSC and Macro_SPP1 in a hypoxic region that exacerbates the tumor microenvironment in HCC.

Revealing the role of SPP1+ macrophages in glioma prognosis and therapeutic targeting by investigating tumor-associated macrophage landscape in grade 2 and 3 gliomas

BACKGROUND

Glioma is a highly heterogeneous brain tumor categorized into World Health Organization (WHO) grades 1-4 based on its malignancy. The suppressive immune microenvironment of glioma contributes significantly to unfavourable patient outcomes. However, the cellular composition and their complex interplays within the glioma environment remain poorly understood, and reliable prognostic markers remain elusive. Therefore, in-depth exploration of the tumor microenvironment (TME) and identification of predictive markers are crucial for improving the clinical management of glioma patients.

RESULTS

Our analysis of single-cell RNA-sequencing data from glioma samples unveiled the immunosuppressive role of tumor-associated macrophages (TAMs), mediated through intricate interactions with tumor cells and lymphocytes. We also discovered the heterogeneity within TAMs, among which a group of suppressive TAMs named TAM-SPP1 demonstrated a significant association with Epidermal Growth Factor Receptor (EGFR) amplification, impaired T cell response and unfavourable patient survival outcomes. Furthermore, by leveraging genomic and transcriptomic data from The Cancer Genome Atlas (TCGA) dataset, two distinct molecular subtypes with a different constitution of TAMs, EGFR status and clinical outcomes were identified. Exploiting the molecular differences between these two subtypes, we developed a four-gene-based prognostic model. This model displayed strong associations with an elevated level of suppressive TAMs and could be used to predict anti-tumor immune response and prognosis in glioma patients.

CONCLUSION

Our findings illuminated the molecular and cellular mechanisms that shape the immunosuppressive microenvironment in gliomas, providing novel insights into potential therapeutic targets. Furthermore, the developed prognostic model holds promise for predicting immunotherapy response and assisting in more precise risk stratification for glioma patients.

Single-cell RNA sequencing of cervical exfoliated cells reveals potential biomarkers and cellular pathogenesis in cervical carcinogenesis

Cervical cancer (CC) is a common gynecological malignancy. Despite the current screening methods have been proved effectively and significantly decreased CC morbidity and mortality, deficiencies still exist. Single-cell RNA sequencing (scRNA-seq) approach can identify the complex and rare cell populations at single-cell resolution. By scRNA-seq, the heterogeneity of tumor microenvironment across cervical carcinogenesis has been mapped and described. Whether these alterations could be detected and applied to CC screening is unclear. Herein, we performed scRNA-seq of 56,173 cervical exfoliated cells from 15 samples, including normal cervix, low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), and malignancy. The present study delineated the alteration of immune and epithelial cells derived during the cervical lesion progression. A subset of lipid-associated macrophage was identified as a tumor-promoting element and could serve as a biomarker for predicting the progression of LSIL into HSIL, which was then verified by immunofluorescence. Furthermore, cell-cell communication analysis indicated the SPP1-CD44 axis might exhibit a protumor interaction between epithelial cell and macrophage. In this study, we investigated the cervical multicellular ecosystem in cervical carcinogenesis and identified potential biomarkers for early detection.

Novel tumor-associated macrophage populations and subpopulations by single cell RNA sequencing

Tumor-associated macrophages (TAMs) are present in almost all solid tumor tissues. 16They play critical roles in immune regulation, tumor angiogenesis, tumor stem cell activation, tumor invasion and metastasis, and resistance to therapy. However, it is unclear how TAMs perform these functions. With the application of single-cell RNA sequencing (scRNA-seq), it has become possible to identify TAM subpopulations associated with distinct functions. In this review, we discuss four novel TAM subpopulations in distinct solid tumors based on core gene signatures by scRNA-seq, including FCN1 +, SPP1 +, C1Q + and CCL18 + TAMs. Functional enrichment and gene expression in scRNA-seq data from different solid tumor tissues found that FCN1 + TAMs may induce inflammation; SPP1 + TAMs are potentially involved in metastasis, angiogenesis, and cancer cell stem cell activation, whereas C1Q + TAMs participate in immune regulation and suppression; And CCL18 + cells are terminal immunosuppressive macrophages that not only have a stronger immunosuppressive function but also enhance tumor metastasis. SPP1 + and C1Q + TAM subpopulations can be further divided into distinct populations with different functions. Meanwhile, we will also present emerging evidence highlighting the separating macrophage subpopulations associated with distinct functions. However, there exist the potential disconnects between cell types and subpopulations identified by scRNA-seq and their actual function.

Examining the function of macrophage oxidative stress response and immune system in glioblastoma multiforme through analysis of single-cell transcriptomics

Background

Glioblastoma (GBM), a prevalent malignant neoplasm within the neuro-oncological domain, has been a subject of considerable scrutiny. Macrophages, serving as the principal immunological constituents, profoundly infiltrate the microenvironment of GBM. However, investigations elucidating the intricate immunological mechanisms governing macrophage involvement in GBM at the single-cell level remain notably limited.

Methods

We conducted a comprehensive investigation employing single-cell analysis, aiming to redefine the intricate cellular landscape within both the core and peripheral regions of GBM tumors. Our analytical focus extended to the profound study of macrophages, elucidating their roles within the context of oxidative stress, intercellular information exchange, and cellular trajectories concerning GBM and its assorted subpopulations. We pursued the identification of GBM prognostic genes intricately associated with macrophages. Utilizing experimental research to investigate the relevance of MANBA in the context of GBM.

Results

Our investigations have illuminated the central role of macrophages in the intricate interplay among various subpopulations within the GBM microenvironment. Notably, we observed a pronounced intensity of oxidative stress responses within macrophages when compared to their GBM counterparts in other subpopulations. Moreover, macrophages orchestrated intricate cellular communication networks, facilitated by the SPP1-CD44 axis, both internally and with neighboring subpopulations. These findings collectively suggest the potential for macrophage polarization from an M1 to an M2 phenotype, contributing to immune suppression within the tumor microenvironment. Furthermore, our exploration unearthed GBM prognostic genes closely associated with macrophages, most notably MANBA and TCF12. Remarkably, MANBA appears to participate in the modulation of neuroimmune functionality by exerting inhibitory effects on M1-polarized macrophages, thereby fostering tumor progression. To bolster these assertions, experimental validations unequivocally affirmed the promotional impact of MANBA on GBM, elucidated through its capacity to curb cell proliferation, invasiveness, and metastatic potential.

Conclusion

These revelations represent a pivotal step towards unraveling the intricate immunological mechanisms governing the interactions between macrophages and diverse subpopulations within the GBM milieu. Furthermore, they lay the foundation for the development of an innovative GBM prognostic model, with MANBA at its epicenter, and underscore the potential for novel immunotherapeutic targets in the ongoing pursuit of enhanced treatment modalities for this formidable malignancy.

CILP is a potential pan-cancer marker: combined silico study and in vitro analyses

CILP (Cartilage intermediate layer protein), an ECM (extracellular matrix) glycoprotein, is found to be associated with intervertebral disc degeneration, chronic heart failure, obese and cardiac fibrosis. However, there are few reports on the role of CILP in tumors. Thus, in this study, we mainly explored the function of CILP in the occurrence and development of tumors and whether it could be a potential pan-cancer marker. Pan-cancer data in this study were obtained from UCSC Xena. Single-cell data were obtained from GSE152938. ROC (Receiver operating characteristic) curves were used to evaluate the accuracy of CILP in predicting the occurrence of different tumor types. The Kaplan-Meier plots were used to assess the relationship between CILP expression and survival prognosis in different tumor types by COX regression analysis. Pseudotime analysis and cell communication analysis were used to further explore the function of CILP at Single cell level. The human RCC (renal cell carcinoma) cell lines ACHN and 786-O were used for further experimental verification. Bulk RNA-seq showed differences in CILP expression in several tumors. ROC curves showed that 14 tumors have AUC > 0.7. Kaplan-Meier plots indicated that CILP is a risk factor for patients in 3 kinds of tumors. ScRNA-seq (Single cell RNA sequencing) suggested that CILP might influence tumors through fibroblasts and cell-cell communication. Finally, we verified the function of CILP at the cellular level by using RCC cell lines ACHN and 786-O and found that knockdown of CILP could significantly inhibit migration and invasion. This finding supports that CILP could be a risk factor as well as a pan-cancer predictor for patients.

State transition and intercellular communication of synovial fibroblasts in response to chronic and acute shoulder injuries unveiled by single-cell transcriptomic analyses

PURPOSE

We aimed to investigate the heterogeneity of synovial fibroblasts and their potential to undergo cell state transitions at the resolution of single cells.

MATERIALS AND METHODS

We employed the single-cell RNA sequencing (scRNA-seq) approach to comprehensively map the cellular landscape of the shoulder synovium in individuals with chronic rotator cuff tears (RCTs) and acute proximal humerus fractures (PHFs). Utilizing unbiased clustering analysis, we successfully identified distinct subpopulations of fibroblasts within the synovial environment. We utilized Monocle 3 to delineate the trajectory of synovial fibroblast state transition. And we used CellPhone DB v2.1.0 to predict cell-cell communication patterns within the synovial microenvironment.

RESULTS

We identified eight main cell clusters in the shoulder synovium. Unbiased clustering analysis identified four synovial fibroblast subpopulations, with diverse biological functions associated with protein secretion, ECM remodeling, inflammation regulation and cell division. Lining, mesenchymal, pro-inflammatory and proliferative fibroblasts subsets were identified. Combining the results from StemID and characteristic gene features, mesenchymal fibroblasts exhibited characteristics of fibroblast progenitor cells. The trajectory of synovial fibroblast state transition showed a transition from mesenchymal to pro-inflammatory and lining phenotypes. In addition, the cross talk between fibroblast subclusters increased in degenerative shoulder diseases compared to acute trauma.

CONCLUSION

We successfully generated the scRNA-seq transcriptomic atlas of the shoulder synovium, which provides a comprehensive understanding of the heterogeneity of synovial fibroblasts, their potential to undergo state transitions, and their intercellular communication in the context of chronic degenerative and acute traumatic shoulder diseases.

Single-cell deconvolution reveals high lineage- and location-dependent heterogeneity in mesenchymal multivisceral stage 4 colorectal cancer

Metastasized colorectal cancer (CRC) is associated with a poor prognosis and rapid disease progression. Besides hepatic metastasis, peritoneal carcinomatosis is the major cause of death in Union for International Cancer Control (UICC) stage IV CRC patients. Insights into differential site-specific reconstitution of tumor cells and the corresponding tumor microenvironment are still missing. Here, we analyzed the transcriptome of single cells derived from murine multivisceral CRC and delineated the intermetastatic cellular heterogeneity regarding tumor epithelium, stroma, and immune cells. Interestingly, we found an intercellular site-specific network of cancer-associated fibroblasts and tumor epithelium during peritoneal metastasis as well as an autologous feed-forward loop in cancer stem cells. We furthermore deciphered a metastatic dysfunctional adaptive immunity by a loss of B cell-dependent antigen presentation and consecutive effector T cell exhaustion. Furthermore, we demonstrated major similarities of this murine metastatic CRC model with human disease and - based on the results of our analysis - provided an auspicious site-specific immunomodulatory treatment approach for stage IV CRC by intraperitoneal checkpoint inhibition.

Tumor macrophage functional heterogeneity can inform the development of novel cancer therapies

Macrophages represent a key component of the tumor microenvironment (TME) and are largely associated with poor prognosis. Therapeutic targeting of macrophages has historically focused on inhibiting their recruitment or reprogramming their phenotype from a protumor (M2-like) to an antitumor (M1-like) one. Unfortunately, this approach has not provided clinical breakthroughs that have changed practice. Emerging studies utilizing single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics have improved our understanding of the ontogeny, phenotype, and functional plasticity of macrophages. Overlaying the wealth of current information regarding macrophage molecular subtypes and functions has also identified novel therapeutic vulnerabilities that might drive better control of tumor-associated macrophages (TAMs). Here, we discuss the functional profiling of macrophages and provide an update of novel macrophage-targeted therapies in development.

Unexplained Causes of Glioma-Associated Epilepsies: A Review of Theories and an Area for Research

Approximately 30% of glioma patients are able to survive beyond one year postdiagnosis. And this short time is often overshadowed by glioma-associated epilepsy. This condition severely impairs the patient's quality of life and causes great suffering. The genetic, molecular and cellular mechanisms underlying tumour development and epileptogenesis remain incompletely understood, leading to numerous unanswered questions. The various types of gliomas, namely glioblastoma, astrocytoma and oligodendroglioma, demonstrate distinct seizure susceptibility and disease progression patterns. Patterns have been identified in the presence of IDH mutations and epilepsy, with tumour location in cortical regions, particularly the frontal lobe, showing a more frequent association with seizures. Altered expression of TP53, MGMT and VIM is frequently detected in tumour cells from individuals with epilepsy associated with glioma. However, understanding the pathogenesis of these modifications poses a challenge. Moreover, hypoxic effects induced by glioma and associated with the HIF-1a factor may have a significant impact on epileptogenesis, potentially resulting in epileptiform activity within neuronal networks. We additionally hypothesise about how the tumour may affect the functioning of neuronal ion channels and contribute to disruptions in the blood-brain barrier resulting in spontaneous depolarisations.

An integrated study to decipher immunosuppressive cellular communication in the PDAC environment

Pancreatic ductal adenocarcinoma (PDAC) is one the most aggressive cancers and characterized by a highly rigid and immunosuppressive tumor microenvironment (TME). The extensive cellular interactions are known to play key roles in the immune evasion, chemoresistance, and poor prognosis. Here, we used the spatial transcriptomics, scRNA-seq, and bulk RNA-seq datasets to enhance the insights obtained from each to decipher the cellular communication in the TME. The complex crosstalk in PDAC samples was revealed by the single-cell and spatial transcriptomics profiles of the samples. We show that tumor-associated macrophages (TAMs) are the central cell types in the regulation of microenvironment in PDAC. They colocalize with the cancer cells and tumor-suppressor immune cells and take roles to provide an immunosuppressive environment. LGALS9 gene which is upregulated in PDAC tumor samples in comparison to healthy samples was also found to be upregulated in TAMs compared to tumor-suppressor immune cells in cancer samples. Additionally, LGALS9 was found to be the primary component in the crosstalk between TAMs and the other cells. The widespread expression of P4HB gene and its interaction with LGALS9 was also notable. Our findings point to a profound role of TAMs via LGALS9 and its interaction with P4HB that should be considered for further elucidation as target in the combinatory immunotherapies for PDAC.

scRNA-seq and proteomics reveal the distinction of M2-like macrophages between primary and recurrent malignant glioma and its critical role in the recurrence

AIMS

Tumor-associated macrophages (TAMs) in the immune microenvironment play an important role in the increased drug resistance and recurrence of malignant glioma, but the mechanism remains incompletely inventoried. The focus of this study was to investigate the distinctions of M2-like TAMs in the immune microenvironment between primary and recurrent malignant glioma and its influence in the recurrence.

METHODS

We employed single-cell RNA sequencing to construct a single-cell atlas for a total of 23,010 individual cells from 6 patients with primary or recurrent malignant glioma and identified 5 cell types, including TAMs and malignant cells. Immunohistochemical techniques and proteomics analysis were performed to investigate the role of intercellular interaction between malignant cells and TAMs in the recurrence of malignant glioma.

RESULTS

Six subgroups of TAMs were annotated and M2-like TAMs were found to increase in recurrent malignant glioma significantly. A pseudotime trajectory and a dynamic gene expression profiling during the recurrence of malignant glioma were reconstructed. Up-regulation of several cancer pathways and intercellular interaction-related genes are associated with the recurrence of malignant glioma. Moreover, the M2-like TAMs can activate the PI3K/Akt/HIF-1α/CA9 pathway in the malignant glioma cells via SPP1-CD44-mediated intercellular interaction. Interestingly, high expression of CA9 can trigger the immunosuppressive response in the malignant glioma, thus promoting the degree of malignancy and drug resistance.

CONCLUSION

Our study uncovers the distinction of M2-like TAMs between primary and recurrent glioma, which offers unparalleled insights into the immune microenvironment of primary and recurrent malignant glioma.

Machine learning-based identification of lower grade glioma stemness subtypes discriminates patient prognosis and drug response

Glioma stem cells (GSCs) remodel their tumor microenvironment to sustain a supportive niche. Identification and stratification of stemness related characteristics in patients with glioma might aid in the diagnosis and treatment of the disease. In this study, we calculated the mRNA stemness index in bulk and single-cell RNA-sequencing datasets using machine learning methods and investigated the correlation between stemness and clinicopathological characteristics. A glioma stemness-associated score (GSScore) was constructed using multivariate Cox regression analysis. We also generated a GSC cell line derived from a patient diagnosed with glioma and used glioma cell lines to validate the performance of the GSScore in predicting chemotherapeutic responses. Differentially expressed genes (DEGs) between GSCs with high and low GSScores were used to cluster lower-grade glioma (LGG) samples into three stemness subtypes. Differences in clinicopathological characteristics, including survival, copy number variations, mutations, tumor microenvironment, and immune and chemotherapeutic responses, among the three LGG stemness-associated subtypes were identified. Using machine learning methods, we further identified genes as subtype predictors and validated their performance using the CGGA datasets. In the current study, we identified a GSScore that correlated with LGG chemotherapeutic response. Through the score, we also identified a novel classification of the LGG subtype and associated subtype predictors, which might facilitate the development of precision therapy.

Harnessing type I interferon-mediated immunity to target malignant brain tumors

Type I interferons have long been appreciated as a cytokine family that regulates antiviral immunity. Recently, their role in eliciting antitumor immune responses has gained increasing attention. Within the immunosuppressive tumor microenvironment (TME), interferons stimulate tumor-infiltrating lymphocytes to promote immune clearance and essentially reshape a "cold" TME into an immune-activating "hot" TME. In this review, we focus on gliomas, with an emphasis on malignant glioblastoma, as these brain tumors possess a highly invasive and heterogenous brain TME. We address how type I interferons regulate antitumor immune responses against malignant gliomas and reshape the overall immune landscape of the brain TME. Furthermore, we discuss how these findings can translate into future immunotherapies targeting brain tumors in general.

Machine learning-based glycolysis-associated molecular classification reveals differences in prognosis, TME, and immunotherapy for colorectal cancer patients

Background

Aerobic glycolysis is a process that metabolizes glucose under aerobic conditions, finally producing pyruvate, lactic acid, and ATP for tumor cells. Nevertheless, the overall significance of glycolysis-related genes in colorectal cancer and how they affect the immune microenvironment have not been investigated.

Methods

By combining the transcriptome and single-cell analysis, we summarize the various expression patterns of glycolysis-related genes in colorectal cancer. Three glycolysis-associated clusters (GAC) were identified with distinct clinical, genomic, and tumor microenvironment (TME). By mapping GAC to single-cell RNA sequencing analysis (scRNA-seq), we next discovered that the immune infiltration profile of GACs was similar to that of bulk RNA sequencing analysis (bulk RNA-seq). In order to determine the kind of GAC for each sample, we developed the GAC predictor using markers of single cells and GACs that were most pertinent to clinical prognostic indications. Additionally, potential drugs for each GAC were discovered using different algorithms.

Results

GAC1 was comparable to the immune-desert type, with a low mutation probability and a relatively general prognosis; GAC2 was more likely to be immune-inflamed/excluded, with more immunosuppressive cells and stromal components, which also carried the risk of the poorest prognosis; Similar to the immune-activated type, GAC3 had a high mutation rate, more active immune cells, and excellent therapeutic potential.

Conclusion

In conclusion, we combined transcriptome and single-cell data to identify new molecular subtypes using glycolysis-related genes in colorectal cancer based on machine-learning methods, which provided therapeutic direction for colorectal patients.

Cuproptosis combines immune landscape providing prognostic biomarker in head and neck squamous carcinoma

Head and neck squamous carcinomas (HNSC) are the seventh most common cancer around the world. Treatment options available today have considerable limitations in terms of efficacy. Identifying novel therapeutic targets for HNSC is, therefore, urgently needed. As a novel determined regulated cell death (RCD), Cuproptosis is correlated with the development, treatment response, and prognosis of various cancer. However, the potential role of Cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of HNSC remains unclear. To figure out whether TME cells and Cuproptosis could better predict prognosis, in this study, we analyzed the expression, mutation status, and other clinical information of 502 HNSC patients by dividing them into four clusters based on their CRGs and TME cell expression. Utilizing the LASSO-Cox method and bootstrap, we established Prognostic Cuproptosis and TME classifier, which were significantly associated with prognosis, pathways, clinical features, and immune cell infiltration in TME of HNSC. To go further, the subgroup Cup low/TMEhigh displayed a better prognosis than any others. Two GEO datasets demonstrated the proposed risk model's clinical applicability. Our GO enrichment analyses proved the conjoint effect of Cuproptosis and TME on tumor angiogenesis, proliferation, and so on. Single-cell analysis and Immunotherapy profile then provided a foundation for determining the molecular mechanisms. It revealed the prognostic risk score positively correlated with T cell activation and natural killer (NK) recruiting. As far as we know, this study is the first time to explore the involvement of CRGs regulation in the TME of HNSC. In a word, it is vital to use these findings to develop new therapeutic strategies.

The Significance of SPP1 in Lung Cancers and Its Impact as a Marker for Protumor Tumor-Associated Macrophages

Macrophages are a representative cell type in the tumor microenvironment. Macrophages that infiltrate the cancer microenvironment are referred to as tumor-associated macrophages (TAMs). TAMs exhibit protumor functions related to invasion, metastasis, and immunosuppression, and an increased density of TAMs is associated with a poor clinical course in many cancers. Phosphoprotein 1 (SPP1), also known as osteopontin, is a multifunctional secreted phosphorylated glycoprotein. Although SPP1 is produced in a variety of organs, at the cellular level, it is expressed on only a few cell types, such as osteoblasts, fibroblasts, macrophages, dendritic cells, lymphoid cells, and mononuclear cells. SPP1 is also expressed by cancer cells, and previous studies have demonstrated correlations between levels of circulating SPP1 and/or increased SPP1 expression on tumor cells and poor prognosis in many types of cancer. We recently revealed that SPP1 expression on TAMs is correlated with poor prognosis and chemoresistance in lung adenocarcinoma. In this review, we summarize the significance of TAMs in lung cancers and discuss the importance of SPP1 as a new marker for the protumor subpopulation of monocyte-derived TAMs in lung adenocarcinoma. Several studies have shown that the SPP1/CD44 axis contribute to cancer chemoresistance in solid cancers, so the SPP1/CD44 axis may represent one of the most critical mechanisms for cell-to-cell communication between cancer cells and TAMs.

Glioblastoma glycolytic signature predicts unfavorable prognosis, immunological heterogeneity, and ENO1 promotes microglia M2 polarization and cancer cell malignancy

Glioblastomas are the most malignant brain tumors, whose progress was promoted by aberrate aerobic glycolysis. The immune environment was highly engaged in glioblastoma formation, while its interaction with aerobic glycolysis remained unclear. Herein, we build a 7-gene Glycolytic Score (GS) by Elastic Net in the training set and two independent validating sets. The GS predicted malignant features and poor survival with good performances. Immune functional analyses and Cibersort calculation identified depressed T cells, B cells, natural killer cells immunity, and high immunosuppressive cell infiltration in the high-GS group. Also, high expressions of the immune-escape genes were discovered. Subsequently, the single-cell analyses validated the glycolysis-related immunosuppression. The functional results manifested the high-GS neoplastic cells' association with T cells, NK cells, and macrophage function regulation. The intercellular cross-talk showed strong associations between high-GS neoplastic cells and M2 macrophages/microglia in several immunological pathways. We finally confirmed that ENO1, the key gene of the GS, promoted M2 microglia polarization and glioblastoma cell malignant behaviors via immunofluorescence, clone formation, CCK8, and transwell rescue experiments. These results indicated the interactions between cancerous glycolysis and immunosuppression and glycolysis' role in promoting glioblastoma progression. Conclusively, we built a robust model and discovered strong interaction between GS and immune, shedding light on prognosis management improvement and therapeutic strategies development for glioblastoma patients.

RUNX1/CD44 axis regulates the proliferation, migration, and immunotherapy of gliomas: A single-cell sequencing analysis

BACKGROUND

Glioma is one of the most common, primary, and lethal adult brain tumors because of its extreme aggressiveness and poor prognosis. Several recent studies relevant to the immune function of CD44, a transmembrane glycoprotein as a significant hyaluronic acid receptor, have achieved great success, revealing the critical role of CD44 in immune infiltration in gliomas. The overexpression of CD44 has been verified to correlate with cancer aggressiveness and migration, while the clinical and immune features of CD44 expression have not yet been thoroughly characterized in gliomas.

METHODS

Molecular and clinical data of glioma collected from publicly available genomic databases were analyzed.

RESULTS

CD44 was up-expressed in malignant gliomas, notably in the 1p/19q non-codeletion cases, isocitrate dehydrogenase (IDH) wild-type, and mesenchymal subtypes in GBM samples. CD44 expression level strongly correlates with stromal and immune cells, mainly infiltrating the glioma microenvironment by single-cell sequencing analysis. Meanwhile, CD44 can be a promising biomarker in predicting immunotherapy responses and mediating the expression of PD-L1. Finally, RUNX1/CD44 axis could promote the proliferation and migration of gliomas.

CONCLUSIONS

Therefore, CD44 was responsible for glioma growth and progression. It could potentially lead to a novel target for glioma immunotherapy or a prognostic biomarker.

Multiplex immunofluorescence and single-cell transcriptomic profiling reveal the spatial cell interaction networks in the non-small cell lung cancer microenvironment

BACKGROUND

Conventional immunohistochemistry technologies were limited by the inability to simultaneously detect multiple markers and the lack of identifying spatial relationships among cells, hindering understanding of the biological processes in cancer immunology.

METHODS

Tissue slices of primary tumours from 553 IA∼IIIB non-small cell lung cancer (NSCLC) cases were stained by multiplex immunofluorescence (mIF) assay for 10 markers, including CD4, CD38, CD20, FOXP3, CD66b, CD8, CD68, PD-L1, CD133 and CD163, evaluating the amounts of 26 phenotypes of cells in tumour nest and tumour stroma. StarDist depth learning model was utilised to determine the spatial location of cells based on mIF graphs. Single-cell RNA sequencing (scRNA-seq) on four primary NSCLC cases was conducted to investigate the putative cell interaction networks.

RESULTS

Spatial proximity among CD20+ B cells, CD4+ T cells and CD38+ T cells (r2  = 0.41) was observed, whereas the distribution of regulatory T cells was associated with decreased infiltration levels of CD20+ B cells and CD38+ T cells (r2  = -0.45). Univariate Cox analyses identified closer proximity between CD8+ T cells predicted longer disease-free survival (DFS). In contrast, closer proximity between CD133+ cancer stem cells (CSCs), longer distances between CD4+ T cells and CD20+ B cells, CD4+ T cells and neutrophils, and CD20+ B cells and neutrophils were correlated with dismal DFS. Data from scRNA-seq further showed that spatially adjacent N1-like neutrophils could boost the proliferation and activation of T and B lymphocytes, whereas spatially neighbouring M2-like macrophages showed negative effects. An immune-related risk score (IRRS) system aggregating robust quantitative and spatial prognosticators showed that high-IRRS patients had significantly worse DFS than low-IRRS ones (HR 2.72, 95% CI 1.87-3.94, p < .001).

CONCLUSIONS

We developed a framework to analyse the cell interaction networks in tumour microenvironment, revealing the spatial architecture and intricate interplays between immune and tumour cells.

Molecular Crosstalk between Chromatin Remodeling and Tumor Microenvironment in Multiple Myeloma

Multiple myeloma (MM) is a complex disease driven by numerous genetic and epigenetic alterations that are acquired over time. Despite recent progress in the understanding of MM pathobiology and the availability of innovative drugs, which have pronounced clinical outcome, this malignancy eventually progresses to a drug-resistant lethal stage and, thus, novel therapeutic drugs/models always play an important role in effective management of MM. Modulation of tumor microenvironment is one of the hallmarks of cancer biology, including MM, which affects the myeloma genomic architecture and disease progression subtly through chromatin modifications. The bone marrow niche has a prime role in progression, survival, and drug resistance of multiple myeloma cells. Therefore, it is important to develop means for targeting the ecosystem between multiple myeloma bone marrow microenvironment and chromatin remodeling. Extensive gene expression profile analysis has indeed provided the framework for new risk stratification of MM patients and identifying novel molecular targets and therapeutics. However, key tumor microenvironment factors/immune cells and their interactions with chromatin remodeling complex proteins that drive MM cell growth and progression remain grossly undefined.

A Cuproptosis Activation Scoring model predicts neoplasm-immunity interactions and personalized treatments in glioma

Gliomas are malignant tumors in the central nervous system. Cuproptosis is a newly discovered cell death mechanism targeting lipoylated tricarboxylic acid cycle proteins. Previous studies have found that cuproptosis participates in tumor progression, but its role in gliomas is still elusive. Here, we systematically explored the bulk-tumor and single-cell transcriptome data to reveal its role in gliomas. The cuproptosis activity score (CuAS) was constructed based on cuproptosis-related genes, and machine learning techniques validated the score stability. High CuAS gliomas were more likely to have a poor prognosis and an aggressive mesenchymal (MES) subtype. Subsequently, the SCENIC algorithm predicted 20 CuAS-related transcription factors (TFs) in gliomas. Function enrichment and microenvironment analyses found that CuAS was associated with tumor immune infiltration. Accordingly, intercellular communications between neoplasm and immunity were explored by the R package "Cellchat". Five signaling pathways and 8 ligand-receptor pairs including ICAM1, ITGAX, ITGB2, ANXA1-FRR1, and the like, were identified to suggest how cuproptosis activity connected neoplastic and immune cells. Critically, 13 potential drugs targeting high CuAs gliomas were predicted according to the CTRP and PRISM databases, including oligomycin A, dihydroartemisinin, and others. Taken together, cuproptosis is involved in glioma aggressiveness, neoplasm-immune interactions, and may be used to assist in drug selection.

Indocarbocyanine nanoparticles extravasate and distribute better than liposomes in brain tumors

Glioblastoma (GBM) is the most devastating and aggressive brain tumor in adults. Hidden behind the blood-brain and blood-tumor barriers (BBTB), this invasive type of brain tumor is not readily accessible to nano-sized particles. Here we demonstrate that fluorescent indocarbocyanine lipids (ICLs: DiD, DiI) formulated in PEGylated lipid nanoparticle (PLN) exhibit highly efficient penetration and accumulation in GBM. PLN-formulated ICLs demonstrated more efficient penetration in GBM spheroids and organoids in vitro than liposomal ICLs. Over 82% of the tumor's extravascular area was positive for ICL fluorescence in the PLN group versus 13% in the liposomal group just one hour post-systemic injection in the intracranial GBM model. Forty-eight hours post-injection, PLN-formulated ICLs accumulated in 95% of tumor myeloid-derived suppressor cells and macrophages, 70% of tumor regulatory T cells, 50% of tumor-associated microglia, and 65% of non-immune cells. PLN-formulated ICLs extravasated better than PEGylated liposomal doxorubicin and fluorescent dextran and efficiently accumulated in invasive tumor margins and brain-invading cells. While liposomes were stable in serum in vitro and in vivo, PLNs disassembled before entering tumors, which could explain the differences in their extravasation efficiency. These findings offer an opportunity to improve therapeutic cargo delivery to invasive GBM.

Single-cell RNA sequencing reveals the cellular and molecular changes that contribute to the progression of lung adenocarcinoma

Pure ground glass nodules (GGNs) and solid nodules (SNs) represent early and relatively late stages of lung adenocarcinoma (LUAD) in radiology, respectively. The cellular and molecular characteristics of pure GGNs and SNs have not been comprehensively elucidated. Additionally, the mechanism driving the progression of lung adenocarcinoma from pure GGN to SN in radiology is also elusive. In this study, by analyzing the single-cell transcriptomic profiles of 76,762 cells from four pure GGNs, four SNs, and four normal tissues, we found that anti-tumor immunity mediated by NK and CD8+T cells gradually weakened with the progression of LUAD and humoral immunity mediated by plasma B cells was more active in SNs. Additionally, the proliferation ability of some special epithelial cell increased during the progression process from pure GGN to SN. Furthermore, stromal cells and M2 macrophages could assist the progression of LUAD. Through comprehensive analyses, we revealed dynamic changes in cellular components and intercellular interactions during the progression of LUAD. These findings could facilitate our understanding of LUAD and discovery of novel therapeutic targets.

A novel prognostic 7-methylguanosine signature reflects immune microenvironment and alternative splicing in glioma based on multi-omics analysis

Glioma is the most common type of central nervous system tumor with increasing incidence. 7-methylguanosine (m7G) is one of the diverse RNA modifications that is known to regulate RNA metabolism and its dysregulation was associated with various cancers. However, the expression pattern of m7G regulators and their roles in regulating tumor immune microenvironments (TIMEs) as well as alternative splicing events (ASEs) in glioma has not been reported. In this study, we showed that m7G regulators displayed a close correlation with each other and most of them were differentially expressed between normal and glioma tissues. Two m7G signatures were then constructed to predict the overall survival of both GBM and LGG patients with moderate predictive performance. The risk score calculated from the regression coefficient and expression level of signature genes was proved to be an independent prognostic factor for patients with LGG, thus, a nomogram was established on the risk score and other independent clinical parameters to predict the survival probability of LGG patients. We also investigated the correlation of m7G signatures with TIMEs in terms of immune scores, expression levels of HLA and immune checkpoint genes, immune cell composition, and immune-related functions. While exploring the correlation between signature genes and the ASEs in glioma, we found that EIF4E1B was a key regulator and might play dual roles depending on glioma grade. By incorporating spatial transcriptomic data, we found a cluster of cells featured by high expression of PTN exhibited the highest m7G score and may communicate with adjacent cancer cells via SPP1 and PTN signaling pathways. In conclusion, our work brought novel insights into the roles of m7G modification in TIMEs and ASEs in glioma, suggesting that evaluation of m7G in glioma could predict prognosis. Moreover, our data suggested that blocking SPP1 and PTN pathways might be a strategy for combating glioma.

The Role of Metabolic Plasticity of Tumor-Associated Macrophages in Shaping the Tumor Microenvironment Immunity

Cancer cells possess a high metabolic demand for their rapid proliferation, survival, and progression and thus create an acidic and hypoxic tumor microenvironment (TME) deprived of nutrients. Moreover, acidity within the TME is the central regulator of tumor immunity that influences the metabolism of the immune cells and orchestrates the local and systemic immunity, thus, the TME has a major impact on tumor progression and resistance to anti-cancer therapy. Specifically, myeloid cells, which include myeloid-derived suppressor cells (MDSC), dendritic cells, and tumor-associated macrophages (TAMs), often reprogram their energy metabolism, resulting in stimulating the angiogenesis and immunosuppression of tumors. This review summarizes the recent findings of glucose, amino acids, and fatty acid metabolism changes of the tumor-associated macrophages (TAMs), and how the altered metabolism shapes the TME and anti-tumor immunity. Multiple proton pumps/transporters are involved in maintaining the alkaline intracellular pH which is necessary for the glycolytic metabolism of the myeloid cells and acidic TME. We highlighted the roles of these proteins in modulating the cellular metabolism of TAMs and their potential as therapeutic targets for improving immune checkpoint therapy.