Single-Cell RNA-Seq Analysis of Infiltrating Neoplastic Cells at the Migrating Front of Human Glioblastoma. Cell Rep. 2017;21:1399-1410. [PMID: 29091775 DOI: 10.1016/j.celrep.2017.10.030]

Gene choice in cancer cells is exclusive in ion transport but concurrent in DNA replication

Cancers share common cellular and physiological features. Little is known about whether distinctive gene expression patterns can be displayed at the single-cell level by gene families in cancer cells. The expression of gene homologs within a family can exhibit concurrence and exclusivity. Concurrence can promote all-or-none expression patterns of related genes and underlie alternative physiological states. Conversely, exclusive gene families express the same or similar number of homologs in each cell, allowing a broad repertoire of cell identities to be generated. We show that gene families involved in the cell-cycle and antigen presentation are expressed concurrently. Concurrence in the DNA replication complex MCM reflects the replicative status of cells, including cell lines and cancer-derived organoids. Exclusive expression requires precise regulatory mechanism, but cancer cells retain this form of control for ion homeostasis and extend it to gene families involved in cell migration. Thus, the cell adhesion-based identity of healthy cells is transformed to an identity based on migration in the population of cancer cells, reminiscent of epithelial-mesenchymal transition.

Co-expression of immune checkpoints in glioblastoma revealed by single-nucleus RNA sequencing and spatial transcriptomics

Glioblastoma (GBM) is one of the most malignant tumors of the central nervous system. The pattern of immune checkpoint expression in GBM remains largely unknown. We performed snRNA-Seq and spatial transcriptomic (ST) analyses on untreated GBM samples. 8 major cell types were found in both tumor and adjacent normal tissues, with variations in infiltration grade. Neoplastic cells_6 was identified in malignant cells with high expression of invasion and proliferator-related genes, and analyzed its interactions with microglia, MDM cells and T cells. Significant alterations in ligand-receptor interactions were observed, particularly between Neoplastic cells_6 and microglia, and found prominent expression of VISTA/VSIG3, suggesting a potential mechanism for evading immune system attacks. High expression of TIM-3, VISTA, PSGL-1 and VSIG-3 with similar expression patterns in GBM, may have potential as therapeutic targets. The prognostic value of VISTA expression was cross-validated in 180 glioma patients, and it was observed that patients with high VISTA expression had a poorer prognosis. In addition, multimodal cross analysis integrated SnRNA-seq and ST, revealing complex intracellular communication and mapping the GBM tumor microenvironment. This study reveals novel molecular characteristics of GBM, co-expression of immune checkpoints, and potential therapeutic targets, contributing to improving the understanding and treatment of GBM.

Probing the glioma micro-environment: Analysis using biopsy in combination with ultra-fast cyclic immunolabeling

The interaction between gliomas and the immune system is poorly understood and thus hindering development of effective immunotherapies for glioma patients. The immune response is highly variable during tumor development, and affected by therapies such as surgery, radiation, and chemotherapy. Currently, analysis of these local changes is difficult due to poor accessibility of the tumor and high-morbidity of sampling. In this study, we developed a model for repeat-biopsy in mice to study these local immunological changes over time. Using fine needle biopsy we were able to safely and repeatedly collect cells from intracranial tumors in mice. Ultra-fast cycling technology (FAST) was used for multi-cycle immunofluorescence of retrieved cells, and provided insights in the changing immune response over time. The combination of these techniques can be utilized to study changes in the immune response in glioma or other intracranial diseases over time, and in response to treatment within the same animal.

Spatiotemporal modeling reveals high-resolution invasion states in glioblastoma

BACKGROUND

Diffuse invasion of glioblastoma cells through normal brain tissue is a key contributor to tumor aggressiveness, resistance to conventional therapies, and dismal prognosis in patients. A deeper understanding of how components of the tumor microenvironment (TME) contribute to overall tumor organization and to programs of invasion may reveal opportunities for improved therapeutic strategies.

RESULTS

Towards this goal, we apply a novel computational workflow to a spatiotemporally profiled GBM xenograft cohort, leveraging the ability to distinguish human tumor from mouse TME to overcome previous limitations in the analysis of diffuse invasion. Our analytic approach, based on unsupervised deconvolution, performs reference-free discovery of cell types and cell activities within the complete GBM ecosystem. We present a comprehensive catalogue of 15 tumor cell programs set within the spatiotemporal context of 90 mouse brain and TME cell types, cell activities, and anatomic structures. Distinct tumor programs related to invasion align with routes of perivascular, white matter, and parenchymal invasion. Furthermore, sub-modules of genes serving as program network hubs are highly prognostic in GBM patients.

CONCLUSION

The compendium of programs presented here provides a basis for rational targeting of tumor and/or TME components. We anticipate that our approach will facilitate an ecosystem-level understanding of the immediate and long-term consequences of such perturbations, including the identification of compensatory programs that will inform improved combinatorial therapies.

Exploring the significance and potential mechanisms of hippo pathway-associated genes in prognosis of glioma patients

PURPOSE

Despite the efforts of countless researchers to develop glioma treatment strategies, the current therapeutic effect of glioma is still not ideal, and it is necessary to further explore the mechanism to guide treatment. Thus, this study aims to introduce a novel approach for predicting patient prognosis and guiding further treatment interventions.

METHODS

Initially, we conducted a differential gene expression analysis to identify Hippo pathway-associated genes overexpressed in tumors and determined genes correlated with prognosis. Subsequently, employing cluster analysis, we categorized samples into two groups and performed further analyses including prediction, immune cell infiltration abundance, and drug response rates. We utilized weighted gene co-expression analysis to reveal gene sets with high co-variation, delineate inter-sample gene correlation patterns, and conduct enrichment analysis. Prognostic models were built using ten machine learning algorithms combined in 101 different combinations, followed by evaluation and validation. Immune infiltration analysis, differential expression analysis of depleted T cell-related markers, drug sensitivity analysis, and exploration of pathway dysregulation were performed for different risk groups. Quality control and batch integration were performed, and single-cell data were analyzed using dimensionality reduction clustering algorithms and annotation tools to evaluate the activity of the prognostic model in malignant cells.

RESULTS

We conducted data filtering to identify genes overexpressed in tumors, intersecting these genes with Hippo pathway-related genes, identifying 62 genes correlated with prognosis, and performing cluster analysis to divide tumor tissues into two groups. Cluster 2 exhibited a poorer prognosis and demonstrated differences in immune cell infiltration. Utilizing weighted gene co-expression analysis on Cluster 2, we identified gene modules, conducted functional enrichment analysis, and delineated pathways. Employing a combined model based on ten machine learning algorithm combinations, we selected the optimal prognostic model system and validated the model's predictive ability within the dataset. Through immune-related analysis and drug sensitivity analysis, we uncovered differences in immune infiltration and varying sensitivities to chemotherapy drugs. Additionally, the enrichment analysis of gene set revealed discrepancies in upregulation within relevant pathways between the high and low-risk groups. Finally, annotation and evaluation of malignant cells via single-cell analysis showed increased activity of the prognostic model and variations in distribution across different prognostic levels in malignant cells.

CONCLUSION

This study introduces a novel approach utilizing the Hippo pathway and associated genes for glioma prognosis research, demonstrating the potential and significance of this method in evaluating the outcome for patients with glioma. These findings hold substantial clinical significance in guiding therapy and predicting outcomes for individuals diagnosed with glioma, offering significant clinical utility.

Effect of the mRNA decapping enzyme scavenger (DCPS) inhibitor RG3039 on glioblastoma

BACKGROUND

Patients with glioblastoma (GBM) have a poor prognosis and limited treatment options. The mRNA decapping enzyme scavenger (DCPS) is a cap-hydrolyzing enzyme. The DCPS inhibitor RG3039 exhibited excellent central nervous system bioavailability in vivo and was safe and well tolerated in healthy volunteers in a phase 1 clinical trial. In this study, we investigated the expression of DCPS in GBM and the anti-tumor activity of RG3039 in various preclinical models of GBM.

METHODS

DCPS expression was examined in human GBM and paired peritumoral tissues. Its prognostic role was evaluated together with clinicopathological characteristics of patients. The anti-GBM effect of RG3039 was determined using GBM cell lines, patient-derived organoids, and orthotopic mouse models. The therapeutic mechanisms of DCPS inhibition were explored.

RESULTS

DCPS is overexpressed in GBM and is associated with poor survival of patients with GBM. The DCPS inhibitor RG3039 exhibited robust anti-GBM activities in GBM cell lines, patient-derived organoids and orthotopic mouse models, with drug exposure achievable in humans. Mechanistically, RG3039 downregulated STAT5B expression, thereby suppressing proliferation, survival and colony formation of GBM cells.

CONCLUSIONS

DCPS is a promising target for GBM. Inhibition of DCPS with RG3039 at doses achievable in humans downregulates STAT5B expression and reduces proliferation, survival and colony formation of GBM cells. Given the excellent anti-cancer activity and central nervous system bioavailability in vivo and good tolerance in humans, RG3039 warrants further study as a potential GBM therapy.

Patient-based multilevel transcriptome exploration highlights relevant chemokines and chemokine receptor axes in glioblastoma

Chemokines and their receptors form a complex interaction network, crucial for precise leukocyte positioning and trafficking. In cancer, they promote malignant cell proliferation and survival but are also critical for immune cell infiltration in the tumor microenvironment. Glioblastoma (GBM) is the most common and lethal brain tumor, characterized by an immunosuppressive TME, with restricted immune cell infiltration. A better understanding of chemokine-receptor interactions is therefore essential for improving tumor immunogenicity. In this study, we assessed the expression of all human chemokines in adult-type diffuse gliomas, with particular focus on GBM, based on patient-derived samples. Publicly available bulk RNA sequencing datasets allowed us to identify the chemokines most abundantly expressed in GBM, with regard to disease severity and across different tumor subregions. To gain insight into the chemokines-receptor network at the single cell resolution, we explored GBmap, a curated resource integrating multiple scRNAseq datasets from different published studies. Our study constitutes the first patient-based handbook highlighting the relevant chemokine-receptor crosstalks, which are of significant interest in the perspective of a therapeutic modulation of the TME in GBM.

MALT1 protease inhibition restrains glioblastoma progression by reversing tumor-associated macrophage-dependent immunosuppression

MALT1 protease is an intracellular signaling molecule that promotes tumor progression via cancer cell-intrinsic and cancer cell-extrinsic mechanisms. MALT1 has been mostly studied in lymphocytes, and little is known about its role in tumor-associated macrophages. Here, we show that MALT1 plays a key role in glioblastoma (GBM)-associated macrophages. Mechanistically, GBM tumor cells induce a MALT1-NF-κB signaling axis within macrophages, leading to macrophage migration and polarization toward an immunosuppressive phenotype. Inactivation of MALT1 protease promotes transcriptional reprogramming that reduces migration and restores a macrophage "M1-like" phenotype. Preclinical in vivo analysis shows that MALT1 inhibitor treatment results in increased immuno-reactivity of GBM-associated macrophages and reduced GBM tumor growth. Further, the addition of MALT1 inhibitor to temozolomide reduces immunosuppression in the tumor microenvironment, which may enhance the efficacy of this standard-of-care chemotherapeutic. Together, our findings suggest that MALT1 protease inhibition represents a promising macrophage-targeted immunotherapeutic strategy for the treatment of GBM.

Graphical abstract

Integrated proteomic and glycoproteomic analysis reveals heterogeneity and molecular signatures of brain metastases from lung adenocarcinomas

Brain metastasis is a major cause of poor prognosis and death in lung adenocarcinoma (LUAD); however, the understanding of therapeutic strategies and mechanisms for brain metastases from LUAD (BM-LUAD) remains notably limited, especially at the proteomics levels. To address this issue, we conducted integrated proteomic and glycoproteomic analyses on 49 BM-LUAD tumors, revealing two distinct subtypes of the disease: BM-S1 and BM-S2. Whole exome sequencing analysis revealed that somatic mutations in STK11 and KEAP1, as well as copy number deletions on chr19p13.3, such as STK11, UQCR11, and SLC25A23, were more frequently detected in BM-S2. In BM-S1 tumors, we observed significant infiltration of GFAP + astrocytes, as evidenced by elevated levels of GFAP, GABRA2, GABRG1 and GAP43 proteins and an enrichment of astrocytic signatures in both our proteomic data and external spatial transcriptomic data. Conversely, BM-S2 tumors demonstrated higher levels of PD-1 immune cell infiltration, supported by the upregulation of PD-1 and LAG-3 genes. These findings suggest distinct microenvironmental adaptations required by the different BM-LUAD subtypes. Additionally, we observed unique glycosylation patterns between the subtypes, with increased fucosylation in BM-S1 and enhanced sialylation in BM-S2, primarily affected by glycosylation enzymes such as FUT9, B4GALT1, and ST6GAL1. Specifically, in BM-S2, these sialylation modifications are predominantly localized to the lysosomes, underscoring the critical role of N-glycosylation in the tumor progression of BM-LUAD. Overall, our study not only provides a comprehensive multi-omic data resource but also offers valuable biological insights into BM-LUAD, highlighting potential mechanisms and therapeutic targets for further investigation.

THOC1 complexes with SIN3A to regulate R-loops and promote glioblastoma progression

Glioblastoma (GBM), the most common and aggressive malignant brain tumor in adults, has a median survival of 21 months. To identify drivers of GBM proliferation, we conducted a CRISPR-knockout screen, which revealed THO Complex 1 (THOC1) as a key driver. Knocking down THOC1 significantly reduced GBM cell viability across patient-derived xenograft (PDX) lines, enhancing survival (p<0.01) in primary PDX models. Conversely, overexpressing THOC1 in non-cancerous cells bolstered viability, decreasing survival and causing tumor engraftment in vivo (p<0.01). Further investigation revealed THOC1's robust interaction with SIN3A, a histone deacetylase complex. Histone deacetylation has been previously shown to prevent the buildup of R-loops, structures that form normally during transcription but can be lethal in excess. We found that THOC1-knockdown leads to elevated R-loop levels and reduced histone deacetylation levels. Next, to understand the networks specifically regulated by THOC1-mediated R-loop prevention, we conducted unbiased RNA-sequencing on control and THOC1-knockdown GBM cells. We found that THOC1's role in R-loop prevention primarily affects telomeres, critical regions for cell replication. We further show that THOC1-knockdown results in significantly increased telomeric R-loop levels and shortened telomeres. Ultimately, this study suggests that targeting THOC1 shows promise as a therapeutic strategy to disrupt the delicate R-loop landscape and undermine GBM's replicative potential.

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.

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.

Glioblastoma cells increase expression of notch signaling and synaptic genes within infiltrated brain tissue

Glioblastoma remains one of the deadliest brain malignancies. First-line therapy consists of maximal surgical tumor resection, accompanied by chemotherapy and radiotherapy. Malignant cells escape surgical resection by migrating into the surrounding healthy brain tissue, where they give rise to the recurrent tumor. Based on gene expression, tumor cores can be subtyped into mesenchymal, proneural, and classical tumors, each being associated with differences in genetic alterations and cellular composition. In contrast, the adjacent brain parenchyma where infiltrating malignant cells escape surgical resection is less characterized in patients. Using spatial transcriptomics (n = 11), we show that malignant cells within proneural or mesenchymal tumor cores display spatially organized differences in gene expression, although such differences decrease within the infiltrated brain tissue. Malignant cells residing in infiltrated brain tissue have increased expression of genes related to neurodevelopmental pathways and glial cell differentiation. Our findings provide an updated view of the spatial landscape of glioblastomas and further our understanding of the malignant cells that infiltrate the healthy brain, providing new avenues for the targeted therapy of these cells after surgical resection.

StableMate: a statistical method to select stable predictors in omics data

Identifying statistical associations between biological variables is crucial to understanding molecular mechanisms. Most association studies are based on correlation or linear regression analyses, but the identified associations often lack reproducibility and interpretability due to the complexity and variability of omics datasets, making it difficult to translate associations into meaningful biological hypotheses. We developed StableMate, a regression framework, to address these challenges through a process of variable selection across heterogeneous datasets. Given datasets from different environments, such as experimental batches, StableMate selects environment-agnostic (stable) and environment-specific predictors in predicting the response of interest. Stable predictors represent robust functional dependencies with the response, and can be used to build regression models that make generalizable predictions in unseen environments. We applied StableMate to (i) RNA sequencing data of breast cancer to discover genes that consistently predict estrogen receptor expression across disease status; (ii) metagenomics data to identify microbial signatures that show persistent association with colon cancer across study cohorts; and (iii) single-cell RNA sequencing data of glioblastoma to discern signature genes associated with the development of pro-tumour microglia regardless of cell location. Our case studies demonstrate that StableMate is adaptable to regression and classification analyses and achieves comprehensive characterization of biological systems for different omics data types.

Brain macrophage senescence in glioma

Gliomas are a diverse group of primary central nervous system neoplasms with no curative therapies available. Brain macrophages comprise microglia in the brain parenchyma, border-associated macrophages in the meningeal-choroid plexus-perivascular space and monocyte-derived macrophages infiltrating the brain. With the great improvement of our recognition of brain macrophages, diverse macrophage populations have been found in the context of glioma, which exhibit functional and phenotypic heterogeneity. We have long thought that brain macrophage senescence is detrimental, manifested by specialized forms of persistent cell cycle arrest and chronic low-grade inflammation. Persistent senescence of macrophages may result in immune dysfunction, potentially contributing to glioma initiation and development. Given the crucial roles played by brain macrophages in glioma, we unravel how brain macrophages undergo reprogramming and their contribution to glioma. We outline general molecular alterations and specific biomarkers in senescent brain macrophages, as well as functional changes (such as metabolism, autophagy, phagocytosis, antigen presentation, and infiltration and recruitment). In addition, recent advances in genetic regulation and mechanisms linked to senescent brain macrophages are discussed. In particular, this review emphasizes the contribution of senescent brain macrophages to glioma, which may drive translational efforts to utilize brain macrophages as a prognostic marker or/and treatment target in glioma. An in-depth comprehending of how brain macrophage senescence functionally influences the tumor microenvironment will be key to our development of innovative therapeutics for glioma.

Exploring public cancer gene expression signatures across bulk, single-cell and spatial transcriptomics data with signifinder Bioconductor package

Understanding cancer mechanisms, defining subtypes, predicting prognosis and assessing therapy efficacy are crucial aspects of cancer research. Gene-expression signatures derived from bulk gene expression data have played a significant role in these endeavors over the past decade. However, recent advancements in high-resolution transcriptomic technologies, such as single-cell RNA sequencing and spatial transcriptomics, have revealed the complex cellular heterogeneity within tumors, necessitating the development of computational tools to characterize tumor mass heterogeneity accurately. Thus we implemented signifinder, a novel R Bioconductor package designed to streamline the collection and use of cancer transcriptional signatures across bulk, single-cell, and spatial transcriptomics data. Leveraging publicly available signatures curated by signifinder, users can assess a wide range of tumor characteristics, including hallmark processes, therapy responses, and tumor microenvironment peculiarities. Through three case studies, we demonstrate the utility of transcriptional signatures in bulk, single-cell, and spatial transcriptomic data analyses, providing insights into cell-resolution transcriptional signatures in oncology. Signifinder represents a significant advancement in cancer transcriptomic data analysis, offering a comprehensive framework for interpreting high-resolution data and addressing tumor complexity.

NKG2C/KLRC2 tumor cell expression enhances immunotherapeutic efficacy against glioblastoma

BACKGROUND

Activating and inhibitory receptors of natural killer (NK) cells such as NKp, NKG2, or CLEC are highly relevant to cold tumors including glioblastoma (GBM). Here, we aimed to characterize the expression of these receptors in GBM to gain insight into their potential role as modulators of the intratumoral microenvironment.

METHODS

We performed a transcriptomic analysis of several NK receptors with a focus on the activating receptor encoded by KLRC2, NKG2C, among bulk and single-cell RNA sequencing GBM data sets. We also evaluated the effects of KLRC2-overexpressing GL261 cells in mice treated with or without programmed cell death protein-1 (PD-1) monoclonal antibody (mAb). Finally, we analyzed samples from two clinical trials evaluating PD-1 mAb effects in patients with GBM to determine the potential of NKG2C to serve as a biomarker of response.

RESULTS

We observed significant expression of several inhibitory NK receptors on GBM-infiltrating NK and T cells, which contrasts with the strong expression of KLRC2 on tumor cells, mainly at the infiltrative margin. Neoplastic KLRC2 expression was associated with a reduction in the number of myeloid-derived suppressor cells and with a higher level of tumor-resident lymphocytes. A stronger antitumor activity after PD-1 mAb treatment was observed in NKG2Chigh-expressing tumors both in mouse models and patients with GBM whereas the expression of inhibitory NK receptors showed an inverse association.

CONCLUSIONS

This study explored the role of neoplastic NKG2C/KLRC2 expression in shaping the immune profile of GBM and suggests that it is a predictive biomarker for positive responses to immune checkpoint inhibitor treatment in patients with GBM. Future studies could further validate this finding in prospective trials.

The microenvironmental factors induced invasive tumor cells in glioblastoma

Glioblastoma (GBM) cells have the potential to switch from being "proliferative cells" to peritumoral "invasive cells". Peritumoral GBM cells have highly invasive properties that allow them to survive surgery, leading to recurrence. The mechanisms underlying the manner in which the tumor microenvironment (TME) regulates the invasiveness of GBM remain unclear. Single-cell RNA sequencing analysis revealed heterogeneity in GBM cells, microglia and macrophages. In this study, the Oncostatin M receptor (OSMR) and leukemia inhibitory factor receptor (LIFR) expression indicated higher invasiveness in core GBM cells. Under environmental stress, the expression of OSMR and LIFR were up-regulated with the effect of hypoxic, acidic, and low-glucose conditions in vitro. Functional experiments revealed that TME stress significantly influences the proliferation, migration and invasion of GBM cells. The differences in core/peripheral TMEs in GBM affected the invasive properties, indicating the significant role of OSMR expression within the TME in tumor progression and postoperative therapy.

Sirtuin 5 (SIRT5) Suppresses Tumor Growth by Regulating Mitochondrial Metabolism and Synaptic Remodeling in Gliomas

Sirtuin 5 (SIRT5) is increasingly recognized as a key regulator of cellular metabolism, which is commonly dysregulated in cancer cells, resulting in enhanced proliferation and tumor progression. To investigate the clinicopathologic implications of SIRT5 dysregulation in glioblastoma, we performed comprehensive analyses of transcriptomic data and functional verifications using in vitro and in vivo glioblastoma models. We found that higher SIRT5 expression levels were associated with a favorable prognosis in glioma patients. Knockdown of SIRT5 significantly enhanced glioblastoma cell growth. Our data suggest its potential role in regulating mitochondrial metabolism in gliomas. Furthermore, SIRT5 is also significantly correlated with synaptic remodeling pathways. Our findings indicate a tumor-suppressive role for SIRT5 that extends beyond regulating cancer metabolism, by which it may function through modulating neuroplasticity. Understanding these cellular interactions provides nuanced insights into the multifaceted role of SIRT5 and the broader therapeutic implications of this for the development of novel treatment strategies.

The role of progranulin in macrophages of a glioblastoma model

PURPOSE

Glioblastoma (GBM), characterized by astrocytic tumorigenesis, remains one of the most prognostically challenging tumor types. Targeting entire GBM microenvironment using novel therapeutic factors is currently desired investigation approach. In this study, we focused on progranulin (PGRN), a regulator of diverse cellular functions. Recent studies implicated PGRN in the poor prognostics of GBM patients. However, the specific role of PGRN in the GBM microenvironment remains elusive.

METHODS

We utilized public databases of GBM patient and previous single-cell RNA sequence to examine association between PGRN expression and patient survival/grade, and expression levels of PGRN in each cell constituting the tumor microenvironment. To clarify the role of PGRN in Tumor-associated macrophage (TAM), we examined cell proliferation and expression of some proteins in murine GBM cells when cell supernatants derived from TAM of PGRN knockout (Grn-/-) or wild type mice were treated with murine GBM cells.

RESULTS

Our results reveal significant PGRN expression in macrophages within the GBM environment, suggesting an association between increased PGRN expression in macrophages and tumor malignancy. TAM induction led to PGRN expression enhancement. Treatment with Grn-/- mouse -derived bone marrow-derived macrophage (BMDM) supernatant resulted in diminished GBM cell proliferation and cell cycle- and mesenchymal GBM subtype-associated reduced protein expression. Furthermore, the Grn-/- mouse-derived BMDM supernatant treatment reduced the phosphorylated STAT3 expression in GBM cells, while the expression of IL-6 and IL-10, known STAT3 pathway activators, diminished in Grn-/- mouse-derived BMDMs.

CONCLUSION

Our results suggest that macrophage-derived PGRN is pivotal for fostering malignant transformations within the tumor microenvironment.

ARF4-mediated retrograde trafficking as a driver of chemoresistance in glioblastoma

BACKGROUND

Cellular functions hinge on the meticulous orchestration of protein transport, both spatially and temporally. Central to this process is retrograde trafficking, responsible for targeting proteins to the nucleus. Despite its link to many diseases, the implications of retrograde trafficking in glioblastoma (GBM) are still unclear.

METHODS

To identify genetic drivers of TMZ resistance, we conducted comprehensive CRISPR-knockout screening, revealing ADP-ribosylation factor 4 (ARF4), a regulator of retrograde trafficking, as a major contributor.

RESULTS

Suppressing ARF4 significantly enhanced TMZ sensitivity in GBM patient-derived xenograft (PDX) models, leading to improved survival rates (P < .01) in both primary and recurrent lines. We also observed that TMZ exposure stimulates ARF4-mediated retrograde trafficking. Proteomics analysis of GBM cells with varying levels of ARF4 unveiled the influence of this pathway on EGFR signaling, with increased nuclear trafficking of EGFR observed in cells with ARF4 overexpression and TMZ treatment. Additionally, spatially resolved RNA-sequencing of GBM patient tissues revealed substantial correlations between ARF4 and crucial nuclear EGFR (nEGFR) downstream targets, such as MYC, STAT1, and DNA-PK. Decreased activity of DNA-PK, a DNA repair protein downstream of nEGFR signaling that contributes to TMZ resistance, was observed in cells with suppressed ARF4 levels. Notably, treatment with DNA-PK inhibitor, KU-57788, in mice with a recurrent PDX line resulted in prolonged survival (P < .01), highlighting the promising therapeutic implications of targeting proteins reliant on ARF4-mediated retrograde trafficking.

CONCLUSIONS

Our findings demonstrate that ARF4-mediated retrograde trafficking contributes to the development of TMZ resistance, cementing this pathway as a viable strategy to overcome chemoresistance in GBM.

Microtubule-associated NAV3 regulates invasive phenotypes in glioblastoma cells

Glioblastomas are aggressive brain tumors for which effective therapy is still lacking, resulting in dismal survival rates. These tumors display significant phenotypic plasticity, harboring diverse cell populations ranging from tumor core cells to dispersed, highly invasive cells. Neuron navigator 3 (NAV3), a microtubule-associated protein affecting microtubule growth and dynamics, is downregulated in various cancers, including glioblastoma, and has thus been considered a tumor suppressor. In this study, we challenge this designation and unveil distinct expression patterns of NAV3 across different invasion phenotypes. Using glioblastoma cell lines and patient-derived glioma stem-like cell cultures, we disclose an upregulation of NAV3 in invading glioblastoma cells, contrasting with its lower expression in cells residing in tumor spheroid cores. Furthermore, we establish an association between low and high NAV3 expression and the amoeboid and mesenchymal invasive phenotype, respectively, and demonstrate that overexpression of NAV3 directly stimulates glioblastoma invasive behavior in both 2D and 3D environments. Consistently, we observed increased NAV3 expression in cells migrating along blood vessels in mouse xenografts. Overall, our results shed light on the role of NAV3 in glioblastoma invasion, providing insights into this lethal aspect of glioblastoma behavior.

Intratumoral Cell Heterogeneity in Patient-Derived Glioblastoma Cell Lines Revealed by Single-Cell RNA-Sequencing

Glioblastoma cell lines derived from different patients are widely used in tumor biology research and drug screening. A key feature of glioblastoma is the high level of inter- and intratumor heterogeneity that accounts for treatment resistance. Our aim was to investigate whether intratumor heterogeneity is maintained in cell models. Single-cell RNA sequencing was used to investigate the cellular composition of a tumor sample and six patient-derived glioblastoma cell lines. Three cell lines preserved the mutational profile of the original tumor, whereas three others differed from their precursors. Copy-number variation analysis showed significantly rearranged genomes in all the cell lines and in the tumor sample. The tumor had the most complex cell composition, including cancer cells and microenvironmental cells. Cell lines with a conserved genome had less diverse cellularity, and during cultivation, a relative increase in the stem-cell-derived progenitors was noticed. Cell lines with genomes different from those of the primary tumors mainly contained neural progenitor cells and microenvironmental cells. The establishment of cell lines without the driver mutations that are intrinsic to the original tumors may be related to the selection of clones or cell populations during cultivation. Thus, patient-derived glioblastoma cell lines differ substantially in their cellular profile, which should be taken into account in translational studies.

A syngeneic spontaneous zebrafish model of tp53-deficient, EGFRvIII, and PI3KCAH1047R-driven glioblastoma reveals inhibitory roles for inflammation during tumor initiation and relapse in vivo

High-throughput vertebrate animal model systems for the study of patient-specific biology and new therapeutic approaches for aggressive brain tumors are currently lacking, and new approaches are urgently needed. Therefore, to build a patient-relevant in vivo model of human glioblastoma, we expressed common oncogenic variants including activated human EGFRvIII and PI3KCAH1047R under the control of the radial glial-specific promoter her4.1 in syngeneic tp53 loss-of-function mutant zebrafish. Robust tumor formation was observed prior to 45 days of life, and tumors had a gene expression signature similar to human glioblastoma of the mesenchymal subtype, with a strong inflammatory component. Within early stage tumor lesions, and in an in vivo and endogenous tumor microenvironment, we visualized infiltration of phagocytic cells, as well as internalization of tumor cells by mpeg1.1:EGFP+ microglia/macrophages, suggesting negative regulatory pressure by pro-inflammatory cell types on tumor growth at early stages of glioblastoma initiation. Furthermore, CRISPR/Cas9-mediated gene targeting of master inflammatory transcription factors irf7 or irf8 led to increased tumor formation in the primary context, while suppression of phagocyte activity led to enhanced tumor cell engraftment following transplantation into otherwise immune-competent zebrafish hosts. Altogether, we developed a genetically relevant model of aggressive human glioblastoma and harnessed the unique advantages of zebrafish including live imaging, high-throughput genetic and chemical manipulations to highlight important tumor-suppressive roles for the innate immune system on glioblastoma initiation, with important future opportunities for therapeutic discovery and optimizations.

Glioblastoma Neurovascular Progenitor Orchestrates Tumor Cell Type Diversity

Glioblastoma (GBM) is the deadliest form of primary brain tumor with limited treatment options. Recent studies have profiled GBM tumor heterogeneity, revealing numerous axes of variation that explain the molecular and spatial features of the tumor. Here, we seek to bridge descriptive characterization of GBM cell type heterogeneity with the functional role of individual populations within the tumor. Our lens leverages a gene program-centric meta-atlas of published transcriptomic studies to identify commonalities between diverse tumors and cell types in order to decipher the mechanisms that drive them. This approach led to the discovery of a tumor-derived stem cell population with mixed vascular and neural stem cell features, termed a neurovascular progenitor (NVP). Following in situ validation and molecular characterization of NVP cells in GBM patient samples, we characterized their function in vivo. Genetic depletion of NVP cells resulted in altered tumor cell composition, fewer cycling cells, and extended survival, underscoring their critical functional role. Clonal analysis of primary patient tumors in a human organoid tumor transplantation system demonstrated that the NVP has dual potency, generating both neuronal and vascular tumor cells. Although NVP cells comprise a small fraction of the tumor, these clonal analyses demonstrated that they strongly contribute to the total number of cycling cells in the tumor and generate a defined subset of the whole tumor. This study represents a paradigm by which cell type-specific interrogation of tumor populations can be used to study functional heterogeneity and therapeutically targetable vulnerabilities of GBM.

ID2-ETS2 axis regulates the transcriptional acquisition of pro-tumoral microglia phenotype in glioma

Glioblastoma is a highly aggressive brain tumour that creates an immunosuppressive microenvironment. Microglia, the brain's resident immune cells, play a crucial role in this environment. Glioblastoma cells can reprogramme microglia to create a supportive niche that promotes tumour growth. However, the mechanisms controlling the acquisition of a transcriptome associated with a tumour-supportive microglial reactive state are not fully understood. In this study, we investigated changes in the transcriptional profile of BV2 microglia exposed to C6 glioma cells. RNA-sequencing analysis revealed a significant upregulation of microglial inhibitor of DNA binding 1 (Id1) and Id2, helix-loop-helix negative transcription regulatory factors. The concomitant regulation of microglial ETS proto-oncogene 2, transcription factor (ETS2)-target genes, i.e., Dusp6, Fli1, Jun, Hmox1, and Stab1, led us to hypothesize that ETS2 could be regulated by ID proteins. In fact, ID2-ETS2 protein interactions increased in microglia exposed to glioma cells. In addition, perturbation of the ID2-ETS2 transcriptional axis influenced the acquisition of a microglial tumour-supportive phenotype. ID2 and ETS2 genes were found to be expressed by the tumour-associated microglia isolated from human glioblastoma tumour biopsies. Furthermore, ID2 and ETS2 gene expressions exhibited inverse prognostic values in patients with glioma in cohorts from The Cancer Genome Atlas. Collectively, our findings indicate that the regulation of ETS2 by ID2 plays a role in the transcriptional regulation of microglia in response to stimuli originating from glioblastoma cells, information that could lead to developing therapeutic strategies to manipulate microglial tumour-trophic functions.

Endothelial cell ferroptosis influences IDH wild-type glioblastoma growth in recurrent glioblastoma multiforme patients

Glioblastomas are known for their poor clinical prognosis, with recurrent tumors often exhibiting greater invasiveness and faster growth rates compared to primary tumors. To understand the intratumoral changes driving this phenomenon, we employed single-cell sequencing to analyze the differences between two pairs of primary and recurrent glioblastomas. Our findings revealed an upregulation of ferroptosis in endothelial cells within recurrent tumors, identified by the significant overexpression of the NOX4 gene. Further analysis indicated that knocking down NOX4 in endothelial cells reduced the activity of the ferroptosis pathway. Utilizing conditioned media from endothelial cells with lower ferroptosis activity, we observed a decrease in the growth rate of glioblastoma cells. These results highlighted the complex role of ferroptosis within tumors and suggested that targeting ferroptosis in the treatment of glioblastomas requires careful consideration of its effects on endothelial cells, as it may otherwise produce counterproductive outcomes.

DeepPIG: deep neural network architecture with pairwise connected layers and stochastic gates using knockoff frameworks for feature selection

Selecting relevant feature subsets is essential for machine learning applications. Among the feature selection techniques, the knockoff filter procedure proposes a unique framework that minimizes false discovery rates (FDR). However, employing a deep neural network architecture for a knockoff filter framework requires higher detection power. Using the knockoff filter framework, we present a Deep neural network with PaIrwise connected layers integrated with stochastic Gates (DeepPIG) for the feature selection model. DeepPIG exhibited better detection power in synthetic data than the baseline and recent models such as Deep feature selection using Paired-Input Nonlinear Knockoffs (DeepPINK), Stochastic Gates (STG), and SHapley Additive exPlanations (SHAP) while not violating the preselected FDR level, especially when the signal of the features were weak. The selected features determined by DeepPIG demonstrated superior classification performance compared with the baseline model in real-world data analyses, including the prediction of certain cancer prognosis and classification tasks using microbiome and single-cell datasets. In conclusion, DeepPIG is a robust feature selection approach even when the signals of features are weak. Source code is available at https://github.com/DMCB-GIST/DeepPIG .

GBM Immunotherapy: Macrophage Impacts

BACKGROUND

Glioblastoma (GBM) is an extremely aggressive form of brain tumor with low survival rates. Current treatments such as chemotherapy, radiation, and surgery are problematic due to tumor growth, invasion, and tumor microenvironment. GBM cells are resistant to these standard treatments, and the heterogeneity of the tumor makes it difficult to find a universal approach. Progression of GBM and acquisition of resistance to therapy are due to the complex interplay between tumor cells and the TME. A significant portion of the TME consists of an inflammatory infiltrate, with microglia and macrophages being the predominant cells.

METHODS

Analysis of the literature data over a course of 5 years suggest that the tumor-associated macrophages (TAMs) are capable of releasing cytokines and growth factors that promote tumor proliferation, survival, and metastasis while inhibiting immune cell function at the same time.

RESULTS

Thus, immunosuppressive state, provided with this intensively studied kind of TME cells, is supposed to promote GBM development through TAMs modulation of tumor treatment-resistance and aggressiveness. Therefore, TAMs are an attractive therapeutic target in the treatment of glioblastoma.

CONCLUSION

This review provides a comprehensive overview of the latest research on the nature of TAMs and the development of therapeutic strategies targeting TAMs, focusing on the variety of macrophage properties, being modulated, as well as molecular targets.

Glioblastoma microenvironment-from biology to therapy

Glioblastoma (GBM) is the most aggressive primary brain cancer. These tumors exhibit high intertumoral and intratumoral heterogeneity in neoplastic and nonneoplastic compartments, low lymphocyte infiltration, and high abundance of myeloid subsets that together create a highly protumorigenic immunosuppressive microenvironment. Moreover, heterogeneous GBM cells infiltrate adjacent brain tissue, remodeling the neural microenvironment to foster tumor electrochemical coupling with neurons and metabolic coupling with nonneoplastic astrocytes, thereby driving growth. Here, we review heterogeneity in the GBM microenvironment and its role in low-to-high-grade glioma transition, concluding with a discussion of the challenges of therapeutically targeting the tumor microenvironment and outlining future research opportunities.

Recurrent Glioblastoma-Molecular Underpinnings and Evolving Treatment Paradigms

Glioblastoma is the most common and lethal central nervous system malignancy with a median survival after progression of only 6-9 months. Major biochemical mechanisms implicated in glioblastoma recurrence include aberrant molecular pathways, a recurrence-inducing tumor microenvironment, and epigenetic modifications. Contemporary standard-of-care (surgery, radiation, chemotherapy, and tumor treating fields) helps to control the primary tumor but rarely prevents relapse. Cytoreductive treatment such as surgery has shown benefits in recurrent glioblastoma; however, its use remains controversial. Several innovative treatments are emerging for recurrent glioblastoma, including checkpoint inhibitors, chimeric antigen receptor T cell therapy, oncolytic virotherapy, nanoparticle delivery, laser interstitial thermal therapy, and photodynamic therapy. This review seeks to provide readers with an overview of (1) recent discoveries in the molecular basis of recurrence; (2) the role of surgery in treating recurrence; and (3) novel treatment paradigms emerging for recurrent glioblastoma.

Immunotherapy for pediatric low-grade gliomas

Pediatric low-grade gliomas (pLGGs) are the most common brain tumor types affecting children. Although gross-total resection remains the treatment of choice, many tumors are not amenable to complete removal, because they either involve midline structures, such as the optic chiasm or hypothalamus, and are not conducive to aggressive resection, or have diffuse biological features and blend with the surrounding brain. Historically, radiation therapy was used as the second-line option for disease control, but with the recognition that this often led to adverse long-term sequelae, particularly in young children, conventional chemotherapy assumed a greater role in initial therapy for unresectable tumors. A variety of agents demonstrated activity, but long-term disease control was suboptimal, with more than 50% of tumors exhibiting disease progression within 5 years. More recently, it has been recognized that a high percentage of these tumors in children exhibit constitutive activation of the mitogen-activated protein kinase (MAPK) pathway because of BRAF translocations or mutations, NFI mutations, or a host of other anomalies that converged on MAPK. This led to phase 1, 2, and 3 trials that explored the activity of blocking this signaling pathway, and the efficacy of this approach compared to conventional chemotherapy. Despite initial promise of these strategies, not all children tolerate this therapy, and many tumors resume growth once MAPK inhibition is stopped, raising concern that long-term and potentially life-long treatment will be required to maintain tumor control, even among responders. This observation has led to interest in other treatments, such as immunotherapy, that may delay or avoid the need for additional treatments. This chapter will summarize the place of immunotherapy in the current armamentarium for these tumors and discuss prior results and future options to improve disease control, with a focus on our prior efforts and experience in this field.

Probing the glioma micro-environment: analysis using biopsy in combination with ultra-fast cyclic immunolabeling

The interaction between gliomas and the immune system is poorly understood and thus hindering development of effective immunotherapies for glioma patients. The immune response is highly variable during tumor development, and affected by therapies such as surgery, radiation, and chemotherapy. Currently, analysis of these local changes is difficult due to poor accessibility of the tumor and high-morbidity of sampling. In this study, we developed a model for repeat-biopsy in mice to study these local immunological changes over time. Using fine needle biopsy we were able to safely and repeatedly collect cells from intracranial tumors in mice. Ultra-fast cycling technology (FAST) was used for multi-cycle immunofluorescence of retrieved cells, and provided insights in the changing immune response over time. The combination of these techniques can be utilized to study changes in the immune response in glioma or other intracranial diseases over time, and in response to treatment within the same animal.

Identification of JUN gene and cellular microenvironment in response to PD-1 blockade treatment in lung cancer patients via single-cell RNA sequencing

Exploring the molecular mechanisms of PD-1/PDL-1 blockade for non-small cell lung cancer (NSCLC) would facilitate understanding for tumor microenvironment (TME) and development of individualized medicine. To date, biomarkers of response to PD-1 blockade therapy were still limited. In this study, we hypothesize that cell type in the tumor microenvironment can influence the effect of PD-1 blockade immunotherapy through specific genes. Therefore, we re-analyze the single-cell RNA sequencing data and validation in tissue from lung adenocarcinoma patients. Dynamic changes of cellular subpopulation were observed after anti-PD-1 immunotherapy among TMEs between primary/metastasis or good/poor response patients. Non-exhausted CD8 T cells and dysregulated genes were observed in responsing patients from PD-1 blockade therapy. Among all changed genes, JUN, involved in PD-1 blockade immunotherapy pathway, and could be considered as a PD-1 responsing biomarker.

Revolutionizing Glioblastoma Treatment: A Comprehensive Overview of Modern Therapeutic Approaches

Glioblastoma is the most common malignant primary brain tumor in the adult population, with an average survival of 12.1 to 14.6 months. The standard treatment, combining surgery, radiotherapy, and chemotherapy, is not as efficient as we would like. However, the current possibilities are no longer limited to the standard therapies due to rapid advancements in biotechnology. New methods enable a more precise approach by targeting individual cells and antigens to overcome cancer. For the treatment of glioblastoma, these are gamma knife therapy, proton beam therapy, tumor-treating fields, EGFR and VEGF inhibitors, multiple RTKs inhibitors, and PI3K pathway inhibitors. In addition, the increasing understanding of the role of the immune system in tumorigenesis and the ability to identify tumor-specific antigens helped to develop immunotherapies targeting GBM and immune cells, including CAR-T, CAR-NK cells, dendritic cells, and immune checkpoint inhibitors. Each of the described methods has its advantages and disadvantages and faces problems, such as the inefficient crossing of the blood-brain barrier, various neurological and systemic side effects, and the escape mechanism of the tumor. This work aims to present the current modern treatments of glioblastoma.

Targeting inflammation in glioblastoma: An updated review from pathophysiology to novel therapeutic approaches

Glioblastoma (GBM) is a highly aggressive primary malignant brain tumor with a dismal prognosis despite current treatment strategies. Inflammation plays an essential role in GBM pathophysiology, contributing to tumor growth, invasion, immunosuppression, and angiogenesis. As a result, pharmacological intervention with anti-inflammatory drugs has been used as a potential approach for the management of GBM. To provide an overview of the current understanding of GBM pathophysiology, potential therapeutic applications of anti-inflammatory drugs in GBM, conventional treatments of glioblastoma and emerging therapeutic approaches currently under investigation. A narrative review was carried out, scanning publications from 2000 to 2023 on PubMed and Google Scholar. The search was not guided by a set research question or a specific search method but rather focused on the area of interest. Conventional treatments such as surgery, radiotherapy, and chemotherapy have shown some benefits, but their effectiveness is limited by various factors such as tumor heterogeneity and resistance.

WITHDRAWN: LonP1 Drives Proneural Mesenchymal Transition in IDH1-R132H Diffuse Glioma

The authors have withdrawn their manuscript owing to massive revision and data validation. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author.

Microglia in Glioblastomas: Molecular Insight and Immunotherapeutic Potential

Glioblastoma (GBM) is one of the most aggressive and devastating primary brain tumors, with a median survival of 15 months following diagnosis. Despite the intense treatment regimen which routinely includes maximal safe neurosurgical resection followed by adjuvant radio- and chemotherapy, the disease remains uniformly fatal. The poor prognosis associated with GBM is multifactorial owing to factors such as increased proliferation, angiogenesis, and metabolic switching to glycolytic pathways. Critically, GBM-mediated local and systemic immunosuppression result in inadequate immune surveillance and ultimately, tumor-immune escape. Microglia-the resident macrophages of the central nervous system (CNS)-play crucial roles in mediating the local immune response in the brain. Depending on the specific pathological cues, microglia are activated into either a pro-inflammatory, neurotoxic phenotype, known as M1, or an anti-inflammatory, regenerative phenotype, known as M2. In either case, microglia secrete corresponding pro- or anti-inflammatory cytokines and chemokines that either promote or hinder tumor growth. Herein, we review the interplay between GBM cells and resident microglia with a focus on contemporary studies highlighting the effect of GBM on the subtypes of microglia expressed, the associated cytokines/chemokines secreted, and ultimately, their impact on tumor pathogenesis. Finally, we explore how understanding the intricacies of the tumor-immune landscape can inform novel immunotherapeutic strategies against this devastating disease.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte- like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Hypoxia coordinates the spatial landscape of myeloid cells within glioblastoma to affect survival

Myeloid cells are highly prevalent in glioblastoma (GBM), existing in a spectrum of phenotypic and activation states. We now have limited knowledge of the tumor microenvironment (TME) determinants that influence the localization and the functions of the diverse myeloid cell populations in GBM. Here, we have utilized orthogonal imaging mass cytometry with single-cell and spatial transcriptomic approaches to identify and map the various myeloid populations in the human GBM tumor microenvironment (TME). Our results show that different myeloid populations have distinct and reproducible compartmentalization patterns in the GBM TME that is driven by tissue hypoxia, regional chemokine signaling, and varied homotypic and heterotypic cellular interactions. We subsequently identified specific tumor subregions in GBM, based on composition of identified myeloid cell populations, that were linked to patient survival. Our results provide insight into the spatial organization of myeloid cell subpopulations in GBM, and how this is predictive of clinical outcome.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte-like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Neutrophils in glioma microenvironment: from immune function to immunotherapy

Glioma is a malignant tumor of the central nervous system (CNS). Currently, effective treatment options for gliomas are still lacking. Neutrophils, as an important member of the tumor microenvironment (TME), are widely distributed in circulation. Recently, the discovery of cranial-meningeal channels and intracranial lymphatic vessels has provided new insights into the origins of neutrophils in the CNS. Neutrophils in the brain may originate more from the skull and adjacent vertebral bone marrow. They cross the blood-brain barrier (BBB) under the action of chemokines and enter the brain parenchyma, subsequently migrating to the glioma TME and undergoing phenotypic changes upon contact with tumor cells. Under glycolytic metabolism model, neutrophils show complex and dual functions in different stages of cancer progression, including participation in the malignant progression, immune suppression, and anti-tumor effects of gliomas. Additionally, neutrophils in the TME interact with other immune cells, playing a crucial role in cancer immunotherapy. Targeting neutrophils may be a novel generation of immunotherapy and improve the efficacy of cancer treatments. This article reviews the molecular mechanisms of neutrophils infiltrating the central nervous system from the external environment, detailing the origin, functions, classifications, and targeted therapies of neutrophils in the context of glioma.

Microglia and macrophages in glioblastoma: landscapes and treatment directions

Glioblastoma is the most common primary malignant tumour of the central nervous system and remains uniformly and rapidly fatal. The tumour-associated macrophage (TAM) compartment comprises brain-resident microglia and bone marrow-derived macrophages (BMDMs) recruited from the periphery. Immune-suppressive and tumour-supportive TAM cell states predominate in glioblastoma, and immunotherapies, which have achieved striking success in other solid tumours have consistently failed to improve survival in this 'immune-cold' niche context. Hypoxic and necrotic regions in the tumour core are found to enrich, especially in anti-inflammatory and immune-suppressive TAM cell states. Microglia predominate at the invasive tumour margin and express pro-inflammatory and interferon TAM cell signatures. Depletion of TAMs, or repolarisation towards a pro-inflammatory state, are appealing therapeutic strategies and will depend on effective understanding and classification of TAM cell ontogeny and state based on new single-cell and spatial multi-omic in situ profiling. Here, we explore the application of these datasets to expand and refine TAM characterisation, to inform improved modelling approaches, and ultimately underpin the effective manipulation of function.

Understanding glioblastoma at the single-cell level: Recent advances and future challenges

Glioblastoma, the most aggressive and prevalent form of primary brain tumor, is characterized by rapid growth, diffuse infiltration, and resistance to therapies. Intrinsic heterogeneity and cellular plasticity contribute to its rapid progression under therapy; therefore, there is a need to fully understand these tumors at a single-cell level. Over the past decade, single-cell transcriptomics has enabled the molecular characterization of individual cells within glioblastomas, providing previously unattainable insights into the genetic and molecular features that drive tumorigenesis, disease progression, and therapy resistance. However, despite advances in single-cell technologies, challenges such as high costs, complex data analysis and interpretation, and difficulties in translating findings into clinical practice persist. As single-cell technologies are developed further, more insights into the cellular and molecular heterogeneity of glioblastomas are expected, which will help guide the development of personalized and effective therapies, thereby improving prognosis and quality of life for patients.

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

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.

ID1high/activin Ahigh glioblastoma cells contribute to resistance to anti-angiogenesis therapy through malformed vasculature

Although bevacizumab (BVZ), a representative drug for anti-angiogenesis therapy (AAT), is used as a first-line treatment for patients with glioblastoma (GBM), its efficacy is notably limited. Whereas several mechanisms have been proposed to explain the acquisition of AAT resistance, the specific underlying mechanisms have yet to be sufficiently ascertained. Here, we established that inhibitor of differentiation 1 (ID1)high/activin Ahigh glioblastoma cell confers resistance to BVZ. The bipotent effect of activin A during its active phase was demonstrated to reduce vasculature dependence in tumorigenesis. In response to a temporary exposure to activin A, this cytokine was found to induce endothelial-to-mesenchymal transition via the Smad3/Slug axis, whereas prolonged exposure led to endothelial apoptosis. ID1 tumors showing resistance to BVZ were established to be characterized by a hypovascular structure, hyperpermeability, and scattered hypoxic regions. Using a GBM mouse model, we demonstrated that AAT resistance can be overcome by administering therapy based on a combination of BVZ and SB431542, a Smad2/3 inhibitor, which contributed to enhancing survival. These findings offer valuable insights that could contribute to the development of new strategies for treating AAT-resistant GBM.

MiRNAs from the Dlk1-Dio3 locus and miR-224/452 cluster contribute to glioblastoma tumor heterogeneity

Glioblastoma is one of the most common and aggressive brain tumors and has seen few improvements in patient outcomes. Inter-tumor heterogeneity between tumors of different patients as well as intra-tumor heterogeneity of cells within the same tumor challenge the development of effective drugs. MiRNAs play an essential role throughout the developing brain and regulate many key genes involved in oncogenesis, yet their role in driving many of the processes underlying tumor heterogeneity remains unclear. In this study, we highlight miRNAs from the Dlk1-Dio3 and miR-224/452 clusters which may be expressed cell autonomously and have expression that is associated with cell state genes in glioblastoma, most prominently in neural progenitor-like and mesenchymal-like states respectively. These findings implicate these miRNA clusters as potential regulators of glioblastoma intra-tumoral heterogeneity and may serve as valuable biomarkers for cell state identification.

The ins and outs of microglial cells in brain health and disease

Microglia are the brain's resident macrophages that play pivotal roles in immune surveillance and maintaining homeostasis of the Central Nervous System (CNS). Microglia are functionally implicated in various cerebrovascular diseases, including stroke, aneurysm, and tumorigenesis as they regulate neuroinflammatory responses and tissue repair processes. Here, we review the manifold functions of microglia in the brain under physiological and pathological conditions, primarily focusing on the implication of microglia in glioma propagation and progression. We further review the current status of therapies targeting microglial cells, including their re-education, depletion, and re-population approaches as therapeutic options to improve patient outcomes for various neurological and neuroinflammatory disorders, including cancer.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts

BACKGROUND

A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials.

METHODS

Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry, and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors.

RESULTS

We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components.

CONCLUSIONS

Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.