Single-cell multimodal glioma analyses identify epigenetic regulators of cellular plasticity and environmental stress response

DNA methylation heterogeneity attributable to a complex tumor immune microenvironment prompts prognostic risk in glioma

Gliomas are malignant tumours of the human nervous system with different World Health Organization (WHO) classifications, glioblastoma (GBM) with higher grade and are more malignant than lower-grade glioma (LGG). To dissect how the DNA methylation heterogeneity in gliomas is influenced by the complex cellular composition of the tumour immune microenvironment, we first compared the DNA methylation profiles of purified human immune cells and bulk glioma tissue, stratifying three tumour immune microenvironmental subtypes for GBM and LGG samples from The Cancer Genome Atlas (TCGA). We found that more intermediate methylation sites were enriched in glioma tumour tissues, and used the Proportion of sites with Intermediate Methylation (PIM) to compare intertumoral DNA methylation heterogeneity. A larger PIM score reflected stronger DNA methylation heterogeneity. Enhanced DNA methylation heterogeneity was associated with stronger immune cell infiltration, better survival rates, and slower tumour progression in glioma patients. We then created a Cell-type-associated DNA Methylation Heterogeneity Contribution (CMHC) score to explore the impact of different immune cell types on heterogeneous CpG site (CpGct) in glioma tissues. We identified eight prognosis-related CpGct to construct a risk score: the Cell-type-associated DNA Methylation Heterogeneity Risk (CMHR) score. CMHR was positively correlated with cytotoxic T-lymphocyte infiltration (CTL), and showed better predictive performance for IDH status (AUC = 0.96) and glioma histological phenotype (AUC = 0.81). Furthermore, DNA methylation alterations of eight CpGct might be related to drug treatments of gliomas. In conclusion, we indicated that DNA methylation heterogeneity is associated with a complex tumour immune microenvironment, glioma phenotype, and patient's prognosis.

Immunotherapy drives mesenchymal tumor cell state shift and TME immune response in glioblastoma patients

BACKGROUND

Glioblastoma is a highly aggressive type of brain tumor for which there is no curative treatment available. Immunotherapies have shown limited responses in unselected patients, and there is an urgent need to identify mechanisms of treatment resistance to design novel therapy strategies.

METHODS

Here we investigated the phenotypic and transcriptional dynamics at single-cell resolution during nivolumab immune checkpoint treatment of glioblastoma patients.

RESULTS

We present the integrative paired single-cell RNA-seq analysis of 76 tumor samples from patients in a clinical trial of the PD-1 inhibitor nivolumab and untreated patients. We identify a distinct aggressive phenotypic signature in both tumor cells and the tumor microenvironment in response to nivolumab. Moreover, nivolumab-treatment was associated with an increased transition to mesenchymal stem-like tumor cells, and an increase in TAMs and exhausted and proliferative T cells. We verify and extend our findings in large external glioblastoma dataset (n = 298), develop a latent immune signature and find 18% of primary glioblastoma samples to be latent immune, associated with mesenchymal tumor cell state and TME immune response. Finally, we show that latent immune glioblastoma patients are associated with shorter overall survival following immune checkpoint treatment (P = .0041).

CONCLUSIONS

We find a resistance mechanism signature in one fifth of glioblastoma patients associated with a tumor-cell transition to a more aggressive mesenchymal-like state, increase in TAMs and proliferative and exhausted T cells in response to immunotherapy. These patients may instead benefit from neuro-oncology therapies targeting mesenchymal tumor cells.

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.

Genetic and epigenetic instability as an underlying driver of progression and aggressive behavior in IDH-mutant astrocytoma

In recent years, the classification of adult-type diffuse gliomas has undergone a revolution, wherein specific molecular features now represent defining diagnostic criteria of IDH-wild-type glioblastomas, IDH-mutant astrocytomas, and IDH-mutant 1p/19q-codeleted oligodendrogliomas. With the introduction of the 2021 WHO CNS classification, additional molecular alterations are now integrated into the grading of these tumors, given equal weight to traditional histologic features. However, there remains a great deal of heterogeneity in patient outcome even within these established tumor subclassifications that is unexplained by currently codified molecular alterations, particularly in the IDH-mutant astrocytoma category. There is also significant intercellular genetic and epigenetic heterogeneity and plasticity with resulting phenotypic heterogeneity, making these tumors remarkably adaptable and robust, and presenting a significant barrier to the design of effective therapeutics. Herein, we review the mechanisms and consequences of genetic and epigenetic instability, including chromosomal instability (CIN), microsatellite instability (MSI)/mismatch repair (MMR) deficits, and epigenetic instability, in the underlying biology, tumorigenesis, and progression of IDH-mutant astrocytomas. We also discuss the contribution of recent high-resolution transcriptomics studies toward defining tumor heterogeneity with single-cell resolution. While intratumoral heterogeneity is a well-known feature of diffuse gliomas, the contribution of these various processes has only recently been considered as a potential driver of tumor aggressiveness. CIN has an independent, adverse effect on patient survival, similar to the effect of histologic grade and homozygous CDKN2A deletion, while MMR mutation is only associated with poor overall survival in univariate analysis but is highly correlated with higher histologic/molecular grade and other aggressive features. These forms of genomic instability, which may significantly affect the natural progression of these tumors, response to therapy, and ultimately clinical outcome for patients, are potentially measurable features which could aid in diagnosis, grading, prognosis, and development of personalized therapeutics.

The CEBPB+ glioblastoma subcluster specifically drives the formation of M2 tumor-associated macrophages to promote malignancy growth

Rationale: The heterogeneity of tumor cells within the glioblastoma (GBM) microenvironment presents a complex challenge in curbing GBM progression. Understanding the specific mechanisms of interaction between different GBM cell subclusters and non-tumor cells is crucial. Methods: In this study, we utilized a comprehensive approach integrating glioma single-cell and spatial transcriptomics. This allowed us to examine the molecular interactions and spatial localization within GBM, focusing on a specific tumor cell subcluster, GBM subcluster 6, and M2-type tumor-associated macrophages (M2 TAMs). Results: Our analysis revealed a significant correlation between a specific tumor cell subcluster, GBM cluster 6, and M2-type TAMs. Further in vitro and in vivo experiments demonstrated the specific regulatory role of the CEBPB transcriptional network in GBM subcluster 6, which governs its tumorigenicity, recruitment of M2 TAMs, and polarization. This regulation involves molecules such as MCP1 for macrophage recruitment and the SPP1-Integrin αvβ1-Akt signaling pathway for M2 polarization. Conclusion: Our findings not only deepen our understanding of the formation of M2 TAMs, particularly highlighting the differential roles played by heterogeneous cells within GBM in this process, but also provided new insights for effectively controlling the malignant progression of GBM.

Constitutive type-1 interferons signaling activity in malignant gliomas

PURPOSE

Recent studies revealed a pro-tumor effect of constitutive Type-1 interferons (IFN-I) production and the downstream signaling activity in several malignancies. In contrast, heterogeneity and clinical significance of the signaling activity in gliomas remain unknown. Thus, we aimed to depict the heterogeneity and clinical significance of constitutive Type-1 interferon (IFN-I) production and the downstream signaling activity in gliomas.

METHODS

We utilized multiplex immunofluorescence (mIF) on a 364 gliomas tissue microarray from our cohort. Moreover, we conducted bioinformatic analyses on the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) databases to investigate the heterogeneity and clinical significance of constitutive IFN-I signaling activity in gliomas.

RESULTS

We observed high heterogeneity of the constitutive IFN-I signaling activity among glioma subtypes. Signaling increased with the WHO malignancy grade while decreasing in the gliomas with IDH mutations. Additionally, high IFN-I activity served as an independent predictor of unfavorable outcomes, and global DNA hypermethylation in IDH-mutant gliomas was associated with decreased IFN-I signaling activity. Positive correlations were observed between the IFN-I activity and glioma-associated inflammation, encompassing both anti-tumor and pro-tumor immune responses. Furthermore, the IFN-I activity varied significantly among tumor and immune cells in the glioma microenvironment (GME). Notably, a distinct pattern of IFN-I signaling activity distribution in GME cells was observed among glioma subtypes, and the pattern was independently associated with patient overall survival.

CONCLUSIONS

Constitutive IFN-I signaling activity varies significantly among glioma subtypes and represents a potential indicator for increased glioma inflammation and unfavorable clinical outcomes.

Cancer plasticity in therapy resistance: Mechanisms and novel strategies

Therapy resistance poses a significant obstacle to effective cancer treatment. Recent insights into cell plasticity as a new paradigm for understanding resistance to treatment: as cancer progresses, cancer cells experience phenotypic and molecular alterations, corporately known as cell plasticity. These alterations are caused by microenvironment factors, stochastic genetic and epigenetic changes, and/or selective pressure engendered by treatment, resulting in tumor heterogeneity and therapy resistance. Increasing evidence suggests that cancer cells display remarkable intrinsic plasticity and reversibly adapt to dynamic microenvironment conditions. Dynamic interactions between cell states and with the surrounding microenvironment form a flexible tumor ecosystem, which is able to quickly adapt to external pressure, especially treatment. Here, this review delineates the formation of cancer cell plasticity (CCP) as well as its manipulation of cancer escape from treatment. Furthermore, the intrinsic and extrinsic mechanisms driving CCP that promote the development of therapy resistance is summarized. Novel treatment strategies, e.g., inhibiting or reversing CCP is also proposed. Moreover, the review discusses the multiple lines of ongoing clinical trials globally aimed at ameliorating therapy resistance. Such advances provide directions for the development of new treatment modalities and combination therapies against CCP in the context of therapy resistance.

The Advancement and Application of the Single-Cell Transcriptome in Biological and Medical Research

Single-cell RNA sequencing technology (scRNA-seq) has been steadily developing since its inception in 2009. Unlike bulk RNA-seq, scRNA-seq identifies the heterogeneity of tissue cells and reveals gene expression changes in individual cells at the microscopic level. Here, we review the development of scRNA-seq, which has gone through iterations of reverse transcription, in vitro transcription, smart-seq, drop-seq, 10 × Genomics, and spatial single-cell transcriptome technologies. The technology of 10 × Genomics has been widely applied in medicine and biology, producing rich research results. Furthermore, this review presents a summary of the analytical process for single-cell transcriptome data and its integration with other omics analyses, including genomes, epigenomes, proteomes, and metabolomics. The single-cell transcriptome has a wide range of applications in biology and medicine. This review analyzes the applications of scRNA-seq in cancer, stem cell research, developmental biology, microbiology, and other fields. In essence, scRNA-seq provides a means of elucidating gene expression patterns in single cells, thereby offering a valuable tool for scientific research. Nevertheless, the current single-cell transcriptome technology is still imperfect, and this review identifies its shortcomings and anticipates future developments. The objective of this review is to facilitate a deeper comprehension of scRNA-seq technology and its applications in biological and medical research, as well as to identify avenues for its future development in alignment with practical needs.

Single nucleus transcriptomics, pharmacokinetics, and pharmacodynamics of combined CDK4/6 and mTOR inhibition in a phase 0/1 trial of recurrent high-grade glioma

Outcomes for adult patients with a high-grade glioma continue to be dismal and new treatment paradigms are urgently needed. To optimize the opportunity for discovery, we performed a phase 0/1 dose-escalation clinical trial that investigated tumor pharmacokinetics, pharmacodynamics, and single nucleus transcriptomics following combined ribociclib (CDK4/6 inhibitor) and everolimus (mTOR inhibitor) treatment in recurrent high-grade glioma. Patients with a recurrent high-grade glioma (n = 24) harboring 1) CDKN2A / B deletion or CDK4 / 6 amplification, 2) PTEN loss or PIK3CA mutations, and 3) wild-type retinoblastoma protein (Rb) were enrolled. Patients received neoadjuvant ribociclib and everolimus treatment and no dose-limiting toxicities were observed. The median unbound ribociclib concentrations in Gadolinium non-enhancing tumor regions were 170 nM (range, 65 - 1770 nM) and 634 nM (range, 68 - 2345 nM) in patients receiving 5 days treatment at the daily dose of 400 and 600 mg, respectively. Unbound everolimus concentrations were below the limit of detection (< 0.1 nM) in both enhancing and non-enhancing tumor regions at all dose levels. We identified a significant decrease in MIB1 positive cells suggesting ribociclib-associated cell cycle inhibition. Single nuclei RNAseq (snRNA) based comparisons of 17 IDH-wild-type on-trial recurrences to 31 IDH-wild-type standard of care treated recurrences data demonstrated a significantly lower fraction of cycling and neural progenitor-like (NPC-like) malignant cell populations. We validated the CDK4/6 inhibitor-directed malignant cell state shifts using three patient-derived cell lines. The presented clinical trial highlights the value of integrating pharmacokinetics, pharmacodynamics, and single nucleus transcriptomics to assess treatment effects in phase 0/1 surgical tissues, including malignant cell state shifts. ClinicalTrials.gov identifier: NCT03834740 .

Phenotypic noise and plasticity in cancer evolution

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

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.

A longer time to relapse is associated with a larger increase in differences between paired primary and recurrent IDH wild-type glioblastomas at both the transcriptomic and genomic levels

Glioblastoma (GBM) is the most common malignant brain tumor in adults, which remains incurable and often recurs rapidly after initial therapy. While large efforts have been dedicated to uncover genomic/transcriptomic alternations associated with the recurrence of GBMs, the evolutionary trajectories of matched pairs of primary and recurrent (P-R) GBMs remain largely elusive. It remains challenging to identify genes associated with time to relapse (TTR) and construct a stable and effective prognostic model for predicting TTR of primary GBM patients. By integrating RNA-sequencing and genomic data from multiple datasets of patient-matched longitudinal GBMs of isocitrate dehydrogenase wild-type (IDH-wt), here we examined the associations of TTR with heterogeneities between paired P-R GBMs in gene expression profiles, tumor mutation burden (TMB), and microenvironment. Our results revealed a positive correlation between TTR and transcriptomic/genomic differences between paired P-R GBMs, higher percentages of non-mesenchymal-to-mesenchymal transition and mesenchymal subtype for patients with a short TTR than for those with a long TTR, a high correlation between paired P-R GBMs in gene expression profiles and TMB, and a negative correlation between the fitting level of such a paired P-R GBM correlation and TTR. According to these observations, we identified 55 TTR-associated genes and thereby constructed a seven-gene (ZSCAN10, SIGLEC14, GHRHR, TBX15, TAS2R1, CDKL1, and CD101) prognostic model for predicting TTR of primary IDH-wt GBM patients using univariate/multivariate Cox regression analyses. The risk scores estimated by the model were significantly negatively correlated with TTR in the training set and two independent testing sets. The model also segregated IDH-wt GBM patients into two groups with significantly divergent progression-free survival outcomes and showed promising performance for predicting 1-, 2-, and 3-year progression-free survival rates in all training and testing sets. Our findings provide new insights into the molecular understanding of GBM progression at recurrence and potential targets for therapeutic treatments.

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.

Cellular diversity through space and time: adding new dimensions to GBM therapeutic development

The current median survival for glioblastoma (GBM) patients is only about 16 months, with many patients succumbing to the disease in just a matter of months, making it the most common and aggressive primary brain cancer in adults. This poor outcome is, in part, due to the lack of new treatment options with only one FDA-approved treatment in the last decade. Advances in sequencing techniques and transcriptomic analyses have revealed a vast degree of heterogeneity in GBM, from inter-patient diversity to intra-tumoral cellular variability. These cutting-edge approaches are providing new molecular insights highlighting a critical role for the tumor microenvironment (TME) as a driver of cellular plasticity and phenotypic heterogeneity. With this expanded molecular toolbox, the influence of TME factors, including endogenous (e.g., oxygen and nutrient availability and interactions with non-malignant cells) and iatrogenically induced (e.g., post-therapeutic intervention) stimuli, on tumor cell states can be explored to a greater depth. There exists a critical need for interrogating the temporal and spatial aspects of patient tumors at a high, cell-level resolution to identify therapeutically targetable states, interactions and mechanisms. In this review, we discuss advancements in our understanding of spatiotemporal diversity in GBM with an emphasis on the influence of hypoxia and immune cell interactions on tumor cell heterogeneity. Additionally, we describe specific high-resolution spatially resolved methodologies and their potential to expand the impact of pre-clinical GBM studies. Finally, we highlight clinical attempts at targeting hypoxia- and immune-related mechanisms of malignancy and the potential therapeutic opportunities afforded by single-cell and spatial exploration of GBM patient specimens.

Cancer stem cell hypothesis 2.0 in glioblastoma: Where are we now and where are we going?

Over the past 2 decades, the cancer stem cell (CSC) hypothesis has provided insight into many malignant tumors, including glioblastoma (GBM). Cancer stem cells have been identified in patient-derived tumors and in some mouse models, allowing for a deeper understanding of cellular and molecular mechanisms underlying GBM growth and therapeutic resistance. The CSC hypothesis has been the cornerstone of cellular heterogeneity, providing a conceptual and technical framework to explain this longstanding phenotype in GBM. This hypothesis has evolved to fit recent insights into how cellular plasticity drives tumor growth to suggest that CSCs do not represent a distinct population but rather a cellular state with substantial plasticity that can be achieved by non-CSCs under specific conditions. This has further been reinforced by advances in genomics, including single-cell approaches, that have used the CSC hypothesis to identify multiple putative CSC states with unique properties, including specific developmental and metabolic programs. In this review, we provide a historical perspective on the CSC hypothesis and its recent evolution, with a focus on key functional phenotypes, and provide an update on the definition for its use in future genomic studies.

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.

CACNA1B protects naked mole-rat hippocampal neuron from apoptosis via altering the subcellular localization of Nrf2 after 60Co irradiation

Although radiotherapy is the most effective treatment modality for brain tumors, it always injures the central nervous system, leading to potential sequelae such as cognitive dysfunction. Radiation induces molecular, cellular, and functional changes in neuronal and glial cells. The hippocampus plays a critical role in learning and memory; therefore, concerns about radiation-induced injury are widespread. Multiple studies have focused on this complex problem, but the results have not been fully elucidated. Naked mole rat brains were irradiated with 60Co at a dose of 10 Gy. On 7 days, 14 days, and 28 days after irradiation, hippocampi in the control groups were obtained for next-generation sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were subsequently performed. Venn diagrams revealed 580 differentially expressed genes (DEGs) that were common at different times after irradiation. GO and KEGG analyses revealed that the 580 common DEGs were enriched in molecular transducer activity. In particular, CACNA1B mediated regulatory effects after irradiation. CACNA1B expression increased significantly after irradiation. Downregulation of CACNA1B led to a reduction in apoptosis and reactive oxygen species levels in hippocampal neurons. This was due to the interaction between CACNA1B and Nrf2, which disturbed the normal nuclear localization of Nrf2. In addition, CACNA1B downregulation led to a decrease in the cognitive functions of naked mole rats. These findings reveal the pivotal role of CACNA1B in regulating radiation-induced brain injury and will lead to the development of a novel strategy to prevent brain injury after irradiation.

VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

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.

Hypoxia drives shared and distinct transcriptomic changes in two invasive glioma stem cell lines

Glioblastoma (GBM) is the most common primary malignant cancer of the central nervous system. Insufficient oxygenation (hypoxia) has been linked to GBM invasion and aggression, leading to poor patient outcomes. Hypoxia induces gene expression for cellular adaptations. However, GBM is characterized by high intertumoral (molecular subtypes) and intratumoral heterogeneity (cell states), and it is not well understood to what extent hypoxia triggers patient-specific gene responses and cellular diversity in GBM. Here, we surveyed eight patient-derived GBM stem cell lines for invasion phenotypes in 3D culture, which identified two GBM lines showing increased invasiveness in response to hypoxia. RNA-seq analysis of the two patient GBM lines revealed a set of shared hypoxia response genes concerning glucose metabolism, angiogenesis, and autophagy, but also a large set of patient-specific hypoxia-induced genes featuring cell migration and anti-inflammation, highlighting intertumoral diversity of hypoxia responses in GBM. We further applied the Shared GBM Hypoxia gene signature to single cell RNA-seq datasets of glioma patients, which showed that hypoxic cells displayed a shift towards mesenchymal-like (MES) and astrocyte-like (AC) states. Interestingly, in response to hypoxia, tumor cells in IDH-mutant gliomas displayed a strong shift to the AC state, whereas tumor cells in IDH-wildtype gliomas mainly shifted to the MES state. This distinct hypoxia response of IDH-mutant gliomas may contribute to its more favorable prognosis. Our transcriptomic studies provide a basis for future approaches to better understand the diversity of hypoxic niches in gliomas.

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.

Allele-specific transcriptional effects of subclonal copy number alterations enable genotype-phenotype mapping in cancer cells

Subclonal copy number alterations are a prevalent feature in tumors with high chromosomal instability and result in heterogeneous cancer cell populations with distinct phenotypes. However, the extent to which subclonal copy number alterations contribute to clone-specific phenotypes remains poorly understood. We develop TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population. TreeAlign accurately encodes dosage effects from subclonal copy number alterations, the impact of allelic imbalance on allele-specific transcription, and obviates the need to define genotypic clones from a phylogeny a priori, leading to highly granular definitions of clones with distinct expression programs. These improvements enable clone-clone gene expression comparisons with higher resolution and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.

Integrated profiling identifies ferredoxin 1 as an immune-related biomarker of malignant phenotype in glioma

Background

Glioma, a highly resistant and recurrent type of central nervous system tumor, poses a significant challenge in terms of effective drug treatments and its associated mortality rates. Despite the discovery of Ferredoxin 1 (FDX1) as a crucial participant in cuproptosis, an innovative mechanism of cellular demise, its precise implications for glioma prognosis and tumor immune infiltration remain inadequately elucidated.

Methods

To analyze pan-cancer data, we employed multiple public databases. Gene expression evaluation was performed using tissue microarray (TMA) and single-cell sequencing data. Furthermore, four different approaches were employed to assess the prognostic importance of FDX1 in glioma. We conducted the analysis of differential expression genes (DEGs) and Gene Set Enrichment Analysis (GSEA) to identify immune-related predictive signaling pathways. Somatic mutations were assessed using Tumor Mutation Burden (TMB) and waterfall plots. Immune cell infiltration was evaluated with five different algorithms. Furthermore, we performed in vitro investigations to evaluate the biological roles of FDX1 in glioma.

Results

Glioma samples exhibited upregulation of FDX1, which in turn predicted poor prognosis and was positively associated with unfavorable clinicopathological characteristics. Notably, the top four enriched signaling pathways were immune-related, and the discovery revealed a connection between the expression of FDX1 and the frequency of mutations or the TMB. The FDX1_high group exhibited heightened infiltration of immune cells, and there existed a direct association between the expression of FDX1 and the regulation of immune checkpoint. In vitro experiments demonstrated that FDX1 knockdown reduced proliferation, migration, invasion and transition from G2 to M phase in glioma cells.

Conclusion

In glioma, FDX1 demonstrated a positive association with the advancement of malignancy and changes in the infiltration of immune cells.

Glioblastoma: An Update in Pathology, Molecular Mechanisms and Biomarkers

Glioblastoma multiforme (GBM) is the most common and malignant type of primary brain tumor in adults. Despite important advances in understanding the molecular pathogenesis and biology of this tumor in the past decade, the prognosis for GBM patients remains poor. GBM is characterized by aggressive biological behavior and high degrees of inter-tumor and intra-tumor heterogeneity. Increased understanding of the molecular and cellular heterogeneity of GBM may not only help more accurately define specific subgroups for precise diagnosis but also lay the groundwork for the successful implementation of targeted therapy. Herein, we systematically review the key achievements in the understanding of GBM molecular pathogenesis, mechanisms, and biomarkers in the past decade. We discuss the advances in the molecular pathology of GBM, including genetics, epigenetics, transcriptomics, and signaling pathways. We also review the molecular biomarkers that have potential clinical roles. Finally, new strategies, current challenges, and future directions for discovering new biomarkers and therapeutic targets for GBM will be discussed.

Machine learning-based investigation of regulated cell death for predicting prognosis and immunotherapy response in glioma patients

Glioblastoma is a highly aggressive and malignant type of brain cancer that originates from glial cells in the brain, with a median survival time of 15 months and a 5-year survival rate of less than 5%. Regulated cell death (RCD) is the autonomous and orderly cell death under genetic control, controlled by precise signaling pathways and molecularly defined effector mechanisms, modulated by pharmacological or genetic interventions, and plays a key role in maintaining homeostasis of the internal environment. The comprehensive and systemic landscape of the RCD in glioma is not fully investigated and explored. After collecting 18 RCD-related signatures from the opening literature, we comprehensively explored the RCD landscape, integrating the multi-omics data, including large-scale bulk data, single-cell level data, glioma cell lines, and proteome level data. We also provided a machine learning framework for screening the potentially therapeutic candidates. Here, based on bulk and single-cell sequencing samples, we explored RCD-related phenotypes, investigated the profile of the RCD, and developed an RCD gene pair scoring system, named RCD.GP signature, showing a reliable and robust performance in predicting the prognosis of glioblastoma. Using the machine learning framework consisting of Lasso, RSF, XgBoost, Enet, CoxBoost and Boruta, we identified seven RCD genes as potential therapeutic targets in glioma and verified that the SLC43A3 highly expressed in glioma grades and glioma cell lines through qRT-PCR. Our study provided comprehensive insights into the RCD roles in glioma, developed a robust RCD gene pair signature for predicting the prognosis of glioma patients, constructed a machine learning framework for screening the core candidates and identified the SLC43A3 as an oncogenic role and a prediction biomarker in glioblastoma.

Identifying PLAUR as a Pivotal Gene of Tumor Microenvironment and Regulating Mesenchymal Phenotype of Glioblastoma

The mesenchymal (MES) phenotype of glioblastoma (GBM) is the most aggressive and therapy-resistant subtype of GBM. The MES phenotype transition during tumor progression results from both tumor-intrinsic genetic alterations and tumor-extrinsic microenvironmental factors. In this study, we sought to identify genes that can modulate the MES phenotype via both mechanisms. By integrating weighted gene co-expression network analysis (WGCNA) and the differential expression analysis of hypoxia-immunosuppression-related genes, we identified the plasminogen activator, urokinase receptor (PLAUR) as the hub gene. Functional enrichment analysis and GSVA analysis demonstrated that PLAUR was associated with the MES phenotype of glioma and the hypoxia-immunosuppression-related microenvironmental components. Single-cell sequencing analysis revealed that PLAUR mediated the ligand-receptor interaction between tumor-associated macrophages (TAMs) and glioma cells. Functional experiments in vitro with cell lines or primary glioma cells and xenograft models using BALB/c nude mice confirmed the role of PLAUR in promoting the MES phenotype of GBM. Our findings indicate that PLAUR regulates both glioma cells and tumor cell-extrinsic factors that favor the MES phenotype and suggest that PLAUR might be a potential target for GBM therapy.

Functional Contribution and Clinical Implication of Cancer-Associated Fibroblasts in Glioblastoma

PURPOSE

The abundance and biological contribution of cancer-associated fibroblasts (CAF) in glioblastoma (GBM) are poorly understood. Here, we aim to uncover its molecular signature, cellular roles, and potential tumorigenesis implications.

EXPERIMENTAL DESIGN

We first applied single-cell RNA sequencing (RNA-seq) and bioinformatics analysis to identify and characterize stromal cells with CAF transcriptomic features in human GBM tumors. Then, we performed functional enrichment analysis and in vitro assays to investigate their interactions with malignant GBM cells.

RESULTS

We found that CAF abundance was low but significantly correlated with tumor grade, poor clinical outcome, and activation of extracellular matrix remodeling using three large cohorts containing bulk RNA-seq data and clinical information. Proteomic analysis of a GBM-derived CAF line and its secretome revealed fibronectin (FN1) as a critical candidate factor mediating CAF functions. This was validated using in vitro cellular models, which demonstrated that CAF-conditioned media and recombinant FN1 could facilitate the migration and invasion of GBM cells. In addition, we showed that CAFs were more abundant in the mesenchymal-like state (or subtype) than in other states of GBMs. Interestingly, cell lines resembling the proneural state responded to the CAF signaling better for the migratory and invasive phenotypes.

CONCLUSIONS

Overall, this study characterized the molecular features and functional impacts of CAFs in GBM, alluding to novel cell interactions mediated by CAFs in the GBM microenvironment.

IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy

BACKGROUND

Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur.

RESULTS

Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation.

CONCLUSIONS

We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

A multidimensional atlas of human glioblastoma-like organoids reveals highly coordinated molecular networks and effective drugs

Recent advances in the genomics of glioblastoma (GBM) led to the introduction of molecular neuropathology but failed to translate into treatment improvement. This is largely attributed to the genetic and phenotypic heterogeneity of GBM, which are considered the major obstacle to GBM therapy. Here, we use advanced human GBM-like organoid (LEGO: Laboratory Engineered Glioblastoma-like Organoid) models and provide an unprecedented comprehensive characterization of LEGO models using single-cell transcriptome, DNA methylome, metabolome, lipidome, proteome, and phospho-proteome analysis. We discovered that genetic heterogeneity dictates functional heterogeneity across molecular layers and demonstrates that NF1 mutation drives mesenchymal signature. Most importantly, we found that glycerol lipid reprogramming is a hallmark of GBM, and several targets and drugs were discovered along this line. We also provide a genotype-based drug reference map using LEGO-based drug screen. This study provides new human GBM models and a research path toward effective GBM therapy.

Temporal change of DNA methylation subclasses between matched newly diagnosed and recurrent glioblastoma

The longitudinal transition of phenotypes is pivotal in glioblastoma treatment resistance and DNA methylation emerged as an important tool for classifying glioblastoma phenotypes. We aimed to characterize DNA methylation subclass heterogeneity during progression and assess its clinical impact. Matched tissues from 47 glioblastoma patients were subjected to DNA methylation profiling, including CpG-site alterations, tissue and serum deconvolution, mass spectrometry, and immunoassay. Effects of clinical characteristics on temporal changes and outcomes were studied. Among 47 patients, 8 (17.0%) had non-matching classifications at recurrence. In the remaining 39 cases, 28.2% showed dominant DNA methylation subclass transitions, with 72.7% being a mesenchymal subclass. In general, glioblastomas with a subclass transition showed upregulated metabolic processes. Newly diagnosed glioblastomas with mesenchymal transition displayed increased stem cell-like states and decreased immune components at diagnosis and exhibited elevated immune signatures and cytokine levels in serum. In contrast, tissue of recurrent glioblastomas with mesenchymal transition showed increased immune components but decreased stem cell-like states. Survival analyses revealed comparable outcomes for patients with and without subclass transitions. This study demonstrates a temporal heterogeneity of DNA methylation subclasses in 28.2% of glioblastomas, not impacting patient survival. Changes in cell state composition associated with subclass transition may be crucial for recurrent glioblastoma targeted therapies.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Pan-cancer ion transport signature reveals functional regulators of glioblastoma aggression

Ion channels, transporters, and other ion-flux controlling proteins, collectively comprising the "ion permeome", are common drug targets, however, their roles in cancer remain understudied. Our integrative pan-cancer transcriptome analysis shows that genes encoding the ion permeome are significantly more often highly expressed in specific subsets of cancer samples, compared to pan-transcriptome expectations. To enable target selection, we identified 410 survival-associated IP genes in 33 cancer types using a machine-learning approach. Notably, GJB2 and SCN9A show prominent expression in neoplastic cells and are associated with poor prognosis in glioblastoma, the most common and aggressive brain cancer. GJB2 or SCN9A knockdown in patient-derived glioblastoma cells induces transcriptome-wide changes involving neuron projection and proliferation pathways, impairs cell viability and tumor sphere formation in vitro, perturbs tunneling nanotube dynamics, and extends the survival of glioblastoma-bearing mice. Thus, aberrant activation of genes encoding ion transport proteins appears as a pan-cancer feature defining tumor heterogeneity, which can be exploited for mechanistic insights and therapy development.

Wnt inhibition alleviates resistance to immune checkpoint blockade in glioblastoma

Wnt signaling plays a critical role in the progression and treatment outcome of glioblastoma (GBM). Here, we identified WNT7b as a heretofore unknown mechanism of resistance to immune checkpoint inhibition (αPD1) in GBM patients and murine models. Acquired resistance to αPD1 was found to be associated with the upregulation of Wnt7b and β-catenin protein levels in GBM in patients and in a clinically relevant, stem-rich GBM model. Combining the porcupine inhibitor WNT974 with αPD1 prolonged the survival of GBM-bearing mice. However, this combination had a dichotomous response, with a subset of tumors showing refractoriness. WNT974 and αPD1 expanded a subset of DC3-like dendritic cells (DCs) and decreased the granulocytic myeloid-derived suppressor cells (gMDSCs) in the tumor microenvironment (TME). By contrast, monocytic MDSCs (mMDSCs) increased, while T-cell infiltration remained unchanged, suggesting potential TME-mediated resistance. Our preclinical findings warrant the testing of Wnt7b/β-catenin combined with αPD1 in GBM patients with elevated Wnt7b/β-catenin signaling.

Identification of efferocytosis-related subtypes in gliomas and elucidating their characteristics and clinical significance

Introduction: Gliomas, the most prevalent tumors of the central nervous system, are known for their aggressive nature and poor prognosis. The heterogeneity among gliomas leads to varying responses to the same treatments, even among similar glioma types. In our study, we efferocytosis-related subtypes and explored their characteristics in terms of immune landscape, intercellular communication, and metabolic processes, ultimately elucidating their potential clinical implications. Methods and Results: We first identified efferocytosis-related subtypes in Bulk RNA-seq using the NMF algorithm. We then preliminarily demonstrated the correlation of these subtypes with efferocytosis by examining enrichment scores of cell death pathways, macrophage infiltration, and the expression of immune ligands. Our analysis of single-cell RNA-seq data further supported the association of these subtypes with efferocytosis. Through enrichment analysis, we found that efferocytosis-related subtypes differ from other types of gliomas in terms of immune landscape, intercellular communication, and substance metabolism. Moreover, we found that the efferocytosis-related classification is a prognostic factor with robust predictive performance by calculating the AUC values. We also found that efferocytosis-related subtypes, when compared with other gliomas in drug sensitivity, survival, and TIDE scores, show a clear link to the effectiveness of chemotherapy, radiotherapy, and immunotherapy in glioma patients. Discussion: We identified efferocytosis-related subtypes in gliomas by analyzing the expression of 137 efferocytosis-associated genes, exploring their characteristics in immune landscape, intercellular communication, metabolic processes, and genomic variations. Moreover, we discovered that the classification of efferocytosis-related subtypes has a strong prognostic predictive power and holds potential significance in guiding clinical treatment.

Functional and Molecular Heterogeneity in Glioma Stem Cells Derived from Multiregional Sampling

Glioblastoma (GBM) is an aggressive and highly heterogeneous primary brain tumor. Glioma stem cells represent a subpopulation of tumor cells with stem cell traits that are presumed to be the cause of tumor relapse. There exists complex tumor heterogeneity in drug sensitivity patterns between glioma stem cell (GSC) cultures derived from different patients. Here, we describe that heterogeneity also exists between GSC cultures derived from multiple biopsies within a single tumor. From biopsies harvested within spatially distinct regions representing the entire tumor mass, we established seven GSC cultures and compared their stem cell properties, mutations, gene expression profiles, and drug sensitivity patterns against 115 different anticancer drugs. The results were compared to 14 GSC cultures derived from other patients. Between the multiregional-derived GSC cultures, we observed only minor differences in their phenotype, proliferative capacity, and global gene expression. Further, they displayed intratumoral heterogeneity in mutational profiles and sensitivity patterns to anticancer drugs. This heterogeneity, however, did not exceed the extensive heterogeneity found between GSC cultures derived from other GBM patients. Our results suggest that the use of GSC cultures from one single focal biopsy may underestimate the overall complexity of the GSC population and display the importance of including GSC cultures reflecting the entire tumor mass in drug screening strategies.

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.

The nanoprodrug of polytemozolomide combines with MGMT siRNA to enhance the effect of temozolomide in glioma

Temozolomide (TMZ) is a conventional chemotherapeutic drug for glioma, however, its clinical application and efficacy is severely restricted by its drug resistance properties. O6-methylguanine-DNA methyltransferase (MGMT) is a DNA repair enzyme, which can repair the DNA damage caused by TMZ. A large number of clinical data show that reducing the expression of MGMT can enhance the chemotherapeutic efficacy of TMZ. Therefore, in order to improve the resistance of glioma to TMZ, an angiopep-2 (A2) modified nanoprodrug of polytemozolomide (P(TMZ)n) that combines with MGMT siRNA (siMGMT) targeting MGMT was developed (A2/T/D/siMGMT). It not only increased the amount of TMZ within tumor lesion site, but also reduced MGMT expression in glioma. The in vitro experiments indicated that the A2/T/D/siMGMT effectively enhanced the cellular uptake of TMZ and siMGMT, and resulted in a significant cell apoptosis and cytotoxicity in the glioma cells. The in vivo experiments showed that glioma growth was inhibited and the survival time of animals were prolonged remarkably after A2/T/D/siMGMT was injected via tail vein. The results showed that the therapeutic effect of A2/T/D/siMGMT in the treatment of glioma was significantly improved.

The promoting effect and mechanism of MAD2L2 on stemness maintenance and malignant progression in glioma

BACKGROUND

Glioblastoma, the most common primary malignant tumor of the brain, is associated with poor prognosis. Glioblastoma cells exhibit high proliferative and invasive properties, and glioblastoma stem cells (GSCs) have been shown to play a crucial role in the malignant behavior of glioblastoma cells. This study aims to investigate the molecular mechanisms involved in GSCs maintenance and malignant progression.

METHODS

Bioinformatics analysis was performed based on data from public databases to explore the expression profile of Mitotic arrest deficient 2 like 2 (MAD2L2) and its potential function in glioma. The impact of MAD2L2 on glioblastoma cell behaviors was assessed through cell viability assays (CCK8), colony formation assays, 5-Ethynyl-2'-deoxyuridine (EDU) incorporation assays, scratch assays, and transwell migration/invasion assays. The findings from in vitro experiments were further validated in vivo using xenograft tumor model. GSCs were isolated from the U87 and LN229 cell lines through flow cytometry and the stemness characteristics were verified by immunofluorescence staining. The sphere-forming ability of GSCs was examined using the stem cell sphere formation assay. Bioinformatics methods were conducted to identified the potential downstream target genes of MAD2L2, followed by in vitro experimental validation. Furthermore, potential upstream transcription factors that regulate MAD2L2 expression were confirmed through chromatin immunoprecipitation (ChIP) and dual-luciferase reporter assays.

RESULTS

The MAD2L2 exhibited high expression in glioblastoma samples and showed significant correlation with patient prognosis. In vitro and in vivo experiments confirmed that silencing of MAD2L2 led to decreased proliferation, invasion, and migration capabilities of glioblastoma cells, while decreasing stemness characteristics of glioblastoma stem cells. Conversely, overexpression of MAD2L2 enhanced these malignant behaviors. Further investigation revealed that MYC proto-oncogene (c-MYC) mediated the functional role of MAD2L2 in glioblastoma, which was further validated through a rescue experiment. Moreover, using dual-luciferase reporter gene assays and ChIP assays determined that the upstream transcription factor E2F-1 regulated the expression of MAD2L2.

CONCLUSION

Our study elucidated the role of MAD2L2 in maintaining glioblastoma stemness and promoting malignant behaviors through the regulation of c-MYC, suggesting its potential as a therapeutic target.

Deciphering tumor-infiltrating dendritic cells in the single-cell era

Dendritic cells (DCs) serve as a pivotal link connecting innate and adaptive immunity by processing tumor-derived antigens and activating T cells. The advent of single-cell sequencing has revolutionized the categorization of DCs, enabling a high-resolution characterization of the previously unrecognized diversity of DC populations infiltrating the intricate tumor microenvironment (TME). The application of single-cell sequencing technologies has effectively elucidated the heterogeneity of DCs present in the tumor milieu, yielding invaluable insights into their subpopulation structures and functional diversity. This review provides a comprehensive summary of the current state of knowledge regarding DC subtypes in the TME, drawing from single-cell studies conducted across various human tumors. We focused on the categorization, functions, and interactions of distinct DC subsets, emphasizing their crucial roles in orchestrating tumor-related immune responses. Additionally, we delve into the potential implications of these findings for the identification of predictive biomarkers and therapeutic targets. Enhanced insight into the intricate interplay between DCs and the TME promises to advance our comprehension of tumor immunity and, in turn, pave the way for the development of more efficacious cancer immunotherapies.

Glioblastoma revisited: from neuronal-like invasion to pacemaking

In recent years, two developments have helped us to better understand the fundamental biology of glioblastoma: the description of a striking intratumoral heterogeneity including gene expression-based cell states, and the discovery that neuro-cancer interactions and cancer-intrinsic neurodevelopmental mechanisms are fundamental features of glioblastoma. In this opinion article, we aim to integrate both developments. We explain how two key disease features are characterized by different neural mechanisms related to distinct but plastic cancer cell states: first, the single cell-dominated invasive parts and second, the more solid parts which are dominated by communicating cell networks constantly activated by pacemaker-like glioblastoma cells. The resulting integrative roadmap of molecular and functional heterogeneity contributes to the Cancer Neuroscience of glioblastoma and suggests novel therapeutic strategies.

Overexpression of CD99 is associated with tumor adaptiveness and indicates the tumor recurrence and therapeutic responses in gliomas

Glioma undergoes adaptive changes, leading to poor prognosis and resistance to treatment. CD99 influences the migration and invasion of glioma cells and plays an oncogene role. However, whether CD99 can affect the adaptiveness of gliomas is still lacking in research, making its clinical value underestimated. Here, we enrolled our in-house and public multiomics datasets for bioinformatic analysis and conducted immunohistochemistry staining to investigate the role of CD99 in glioma adaptive response and its clinical implications. CD99 is expressed in more adaptative glioma subtypes and cell states. Under hypoxic conditions, CD99 is upregulated in glioma cells and is associated with angiogenesis and metabolic adaptations. Gliomas with over-expressed CD99 also increased the immunosuppressive tumor-associated macrophages. The relevance with tumor adaptiveness of CD99 presented clinical significance. We discovered that CD99 overexpression is associated with short-time recurrence and validated its prognostic value. Additionally, Glioma patients with high expression of CD99 were resistant to chemotherapy and radiotherapy. The CD99 expression was also related to anti-angiogenic and immune checkpoint inhibitor therapy response. Inhibitors of the PI3K-AKT pathway have therapeutic potential against CD99-overexpressing gliomas. Our study identified CD99 as a biomarker characterizing the adaptive response in glioma. Gliomas with high CD99 expression are highly tolerant to stress conditions such as hypoxia and antitumor immunity, making treatment responses dimmer and tumor progression. Therefore, for patients with CD99-overexpressing gliomas, tumor adaptiveness should be fully considered during treatment to avoid drug resistance, and closer clinical monitoring should be carried out to improve the prognosis.

Programs, Origins, and Niches of Immunomodulatory Myeloid Cells in Gliomas

Gliomas are incurable malignancies notable for an immunosuppressive microenvironment with abundant myeloid cells whose immunomodulatory properties remain poorly defined. Here, utilizing scRNA-seq data for 183,062 myeloid cells from 85 human tumors, we discover that nearly all glioma-associated myeloid cells express at least one of four immunomodulatory activity programs: Scavenger Immunosuppressive, C1Q Immunosuppressive, CXCR4 Inflammatory, and IL1B Inflammatory. All four programs are present in IDH1 mutant and wild-type gliomas and are expressed in macrophages, monocytes, and microglia whether of blood or resident myeloid cell origins. Integrating our scRNA-seq data with mitochondrial DNA-based lineage tracing, spatial transcriptomics, and organoid explant systems that model peripheral monocyte infiltration, we show that these programs are driven by microenvironmental cues and therapies rather than myeloid cell type, origin, or mutation status. The C1Q Immunosuppressive program is driven by routinely administered dexamethasone. The Scavenger Immunosuppressive program includes ligands with established roles in T-cell suppression, is induced in hypoxic regions, and is associated with immunotherapy resistance. Both immunosuppressive programs are less prevalent in lower-grade gliomas, which are instead enriched for the CXCR4 Inflammatory program. Our study provides a framework to understand immunomodulatory myeloid cells in glioma, and a foundation to develop more effective immunotherapies.

Single-Nucleus Transcriptional Profiling of GAD2-Positive Neurons From Mouse Lateral Habenula Reveals Distinct Expression of Neurotransmission- and Depression-Related Genes

Background

Glutamatergic projection neurons of the lateral habenula (LHb) drive behavioral state modulation by regulating the activity of midbrain monoaminergic neurons. Identifying circuit mechanisms that modulate LHb output is of interest for understanding control of motivated behaviors.

Methods

A small population of neurons within the medial subnucleus of the mouse LHb express the GABAergic (gamma-aminobutyric acidergic)-synthesizing enzyme GAD2, and they can inhibit nearby LHb projection neurons; however, these neurons lack markers of classic inhibitory interneurons, and they coexpress the vesicular glutamate transporter VGLUT2. To determine the molecular phenotype of these neurons, we genetically tagged the nuclei of GAD2-positive cells and used fluorescence-activated nuclear sorting to isolate and enrich these nuclei for single-nucleus RNA sequencing.

Results

Our data confirm that GAD2+/VGLUT2+ neurons intrinsic to the LHb coexpress markers of both glutamatergic and GABAergic transmission and that they are transcriptionally distinct from either GABAergic interneurons or habenular glutamatergic neurons. We identify gene expression programs within these cells that show sex-specific differences in expression and that are implicated in major depressive disorder, which has been linked to LHb hyperactivity. Finally, we identify the Ntng2 gene encoding the cell adhesion protein netrin-G2 as a marker of LHb GAD2+/VGLUT2+ neurons and a gene product that may contribute to their target projections.

Conclusions

These data show the value of using genetic enrichment of rare cell types for transcriptome studies, and they advance understanding of the molecular composition of a functionally important class of GAD2+ neurons in the LHb.

Glioma Metabolic Feedback In Situ: A First-In-Human Pharmacodynamic Trial of Difluoromethylornithine + AMXT-1501 Through High-Molecular Weight Microdialysis

BACKGROUND AND OBJECTIVES

No new drug has improved survival for glioblastoma since temozolomide in 2005, due in part to the relative inaccessibility of each patient's individualized tumor biology and its response to therapy. We have identified a conserved extracellular metabolic signature of enhancing high-grade gliomas enriched for guanidinoacetate (GAA). GAA is coproduced with ornithine, the precursor to protumorigenic polyamines through ornithine decarboxylase (ODC). AMXT-1501 is a polyamine transporter inhibitor that can overcome tumoral resistance to the ODC inhibitor, difluoromethylornithine (DFMO). We will use DFMO with or without AMXT-1501 to identify candidate pharmacodynamic biomarkers of polyamine depletion in patients with high-grade gliomas in situ . We aim to determine (1) how blocking polyamine production affects intratumoral extracellular guanidinoacetate abundance and (2) the impact of polyamine depletion on the global extracellular metabolome within live human gliomas in situ.

METHODS

DFMO, with or without AMXT-1501, will be administered postoperatively in 15 patients after clinically indicated subtotal resection for high-grade glioma. High-molecular weight microdialysis catheters implanted into residual tumor and adjacent brain will be used for postoperative monitoring of extracellular GAA and polyamines throughout therapeutic intervention from postoperative day (POD) 1 to POD5. Catheters will be removed on POD5 before discharge.

EXPECTED OUTCOMES

We anticipate that GAA will be elevated in tumor relative to adjacent brain although it will decrease within 24 hours of ODC inhibition with DFMO. If AMXT-1501 effectively increases the cytotoxic impact of ODC inhibition, we expect an increase in biomarkers of cytotoxicity including glutamate with DFMO + AMXT-1501 treatment when compared with DFMO alone.

DISCUSSION

Limited mechanistic feedback from individual patients' gliomas hampers clinical translation of novel therapies. This pilot Phase 0 study will provide in situ feedback during DFMO + AMXT-1501 treatment to determine how high-grade gliomas respond to polyamine depletion.

Epigenetic profiling reveals a strong association between lack of 5-ALA fluorescence and EGFR amplification in IDH-wildtype glioblastoma

Background

5-aminolevulinic acid (5-ALA) fluorescence-guided resection increases the percentage of complete CNS tumor resections and improves the progression-free survival of IDH-wildtype glioblastoma patients. A small subset of IDH-wildtype glioblastoma shows no 5-ALA fluorescence. An explanation for these cases is missing. In this study, we used DNA methylation profiling to further characterize non-fluorescent glioblastomas.

Methods

Patients with newly diagnosed and recurrent IDH-wildtype glioblastoma that underwent surgery were analyzed. The intensity of intraoperative 5-ALA fluorescence was categorized as non-visible or visible. DNA was extracted from tumors and genome-wide DNA methylation patterns were analyzed using Illumina EPIC (850k) arrays. Furthermore, 5-ALA intensity was measured by flow cytometry on human gliomasphere lines (BT112 and BT145).

Results

Of 74 included patients, 12 (16.2%) patients had a non-fluorescent glioblastoma, which were compared to 62 glioblastomas with 5-ALA fluorescence. Clinical characteristics were equally distributed between both groups. We did not find significant differences between DNA methylation subclasses and 5-ALA fluorescence (P = .24). The distribution of cells of the tumor microenvironment was not significantly different between the non-fluorescent and fluorescent tumors. Copy number variations in EGFR and simultaneous EGFRvIII expression were strongly associated with 5-ALA fluorescence since all non-fluorescent glioblastomas were EGFR-amplified (P < .01). This finding was also demonstrated in recurrent tumors. Similarly, EGFR-amplified glioblastoma cell lines showed no 5-ALA fluorescence after 24 h of incubation.

Conclusions

Our study demonstrates an association between non-fluorescent IDH-wildtype glioblastomas and EGFR gene amplification which should be taken into consideration for recurrent surgery and future studies investigating EGFR-amplified gliomas.

Decoding Cancer Evolution: Integrating Genetic and Non-Genetic Insights

The development of cancer begins with cells transitioning from their multicellular nature to a state akin to unicellular organisms. This shift leads to a breakdown in the crucial regulators inherent to multicellularity, resulting in the emergence of diverse cancer cell subpopulations that have enhanced adaptability. The presence of different cell subpopulations within a tumour, known as intratumoural heterogeneity (ITH), poses challenges for cancer treatment. In this review, we delve into the dynamics of the shift from multicellularity to unicellularity during cancer onset and progression. We highlight the role of genetic and non-genetic factors, as well as tumour microenvironment, in promoting ITH and cancer evolution. Additionally, we shed light on the latest advancements in omics technologies that allow for in-depth analysis of tumours at the single-cell level and their spatial organization within the tissue. Obtaining such detailed information is crucial for deepening our understanding of the diverse evolutionary paths of cancer, allowing for the development of effective therapies targeting the key drivers of cancer evolution.

The role of APOBEC3C in modulating the tumor microenvironment and stemness properties of glioma: evidence from pancancer analysis

Background

It is now understood that APOBEC3 family proteins (A3s) are essential in tumor progression, yet their involvement in tumor immunity and stemness across diverse cancer types remains poorly understood.

Methods

In the present study, comprehensive genome-wide statistical and bioinformatic analyses were conducted to elucidate A3 family expression patterns, establishing clinically relevant correlations with prognosis, the tumor microenvironment(TME), immune infiltration, checkpoint blockade, and stemness across cancers. Different experimental techniques were applied, including RT-qPCR, immunohistochemistry, sphere formation assays, Transwell migration assays, and wound-healing assays, to investigate the impact of A3C on low-grade glioma (LGG) and glioblastoma multiforme (GBM), as well as its function in glioma stem cells(GSCs).

Results

Dysregulated expression of A3s was observed in various human cancer tissues. The prognostic value of A3 expression differed across cancer types, with a link to particularly unfavorable outcomes in gliomas. A3s are associated with the the TME and stemness in multiple cancers. Additionally, we developed an independent prognostic model based on A3s expression, which may be an independent prognostic factor for OS in patients with glioma. Subsequent validation underscored a strong association between elevated A3C expression and adverse prognostic outcomes, higher tumor grades, and unfavorable histology in glioma. A potential connection between A3C and glioma progression was established. Notably, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses implicated A3C in immune system-related diseases, with heightened A3C levels contributing to an immunosuppressive tumor microenvironment (TME) in glioma. Furthermore, in vitro experiments substantiated the role of A3C in sustaining and renewing glioma stem cells, as A3C deletion led to diminished proliferation, invasion, and migration of glioma cells.

Conclusion

The A3 family exhibits heterogeneous expression across various cancer types, with its expression profile serving as a predictive marker for overall survival in glioma patients. A3C emerges as a regulator of glioma progression, exerting its influence through modulation of the tumor microenvironment and regulation of stemness.

Integrated clustering signature of genomic heterogeneity, stemness and tumor microenvironment predicts glioma prognosis and immunotherapy response

BACKGROUND

Glioma is the most frequent primary tumor of the central nervous system. The high heterogeneity of glioma tumors enables them to adapt to challenging environments, leading to resistance to treatment. Therefore, to detect the driving factors and improve the prognosis of glioma, it is essential to have a comprehensive understanding of the genomic heterogeneity, stemness, and immune microenvironment of glioma.

METHODS

We classified gliomas into various subtypes based on stemness, genomic heterogeneity, and immune microenvironment consensus clustering analysis. We identified risk hub genes linked to heterogeneous characteristics using WGCNA, LASSO, and multivariate Cox regression analysis and utilized them to create an effective risk model.

RESULTS

We thoroughly investigated the genomic heterogeneity, stemness, and immune microenvironment of glioma and identified the risk hub genes RAB42, SH2D4A, and GDF15 based on the TCGA dataset. We developed a risk model utilizing these genes that can reliably predict the prognosis of glioma patients. The risk signature showed a positive correlation with T cell exhaustion and increased infiltration of immunosuppressive cells, and a negative correlation with the response to immunotherapy. Moreover, we discovered that SH2D4A, one of the risk hub genes, could stimulate the migration and proliferation of glioma cells.

CONCLUSIONS

This study identified risk hub genes and established a risk model by analyzing the genomic heterogeneity, stemness, and immune microenvironment of glioma. Our findings will facilitate the diagnosis and prediction of glioma prognosis and may lead to potential treatment strategies for glioma.

Targeting mitochondrial energetics reverses panobinostat- and marizomib-induced resistance in pediatric and adult high-grade gliomas

In previous studies, we demonstrated that panobinostat, a histone deacetylase inhibitor, and bortezomib, a proteasomal inhibitor, displayed synergistic therapeutic activity against pediatric and adult high-grade gliomas. Despite the remarkable initial response to this combination, resistance emerged. Here, in this study, we aimed to investigate the molecular mechanisms underlying the anticancer effects of panobinostat and marizomib, a brain-penetrant proteasomal inhibitor, and the potential for exploitable vulnerabilities associated with acquired resistance. RNA sequencing followed by gene set enrichment analysis (GSEA) was employed to compare the molecular signatures enriched in resistant compared with drug-naïve cells. The levels of adenosine 5'-triphosphate (ATP), nicotinamide adenine dinucleotide (NAD)+ content, hexokinase activity, and tricarboxylic acid (TCA) cycle metabolites required for oxidative phosphorylation to meet their bioenergetic needs were analyzed. Here, we report that panobinostat and marizomib significantly depleted ATP and NAD+ content, increased mitochondrial permeability and reactive oxygen species generation, and promoted apoptosis in pediatric and adult glioma cell lines at initial treatment. However, resistant cells exhibited increased levels of TCA cycle metabolites, which required for oxidative phosphorylation to meet their bioenergetic needs. Therefore, we targeted glycolysis and the electron transport chain (ETC) with small molecule inhibitors, which displayed substantial efficacy, suggesting that resistant cell survival is dependent on glycolytic and ETC complexes. To verify these observations in vivo, lonidamine, an inhibitor of glycolysis and mitochondrial function, was chosen. We produced two diffuse intrinsic pontine glioma (DIPG) models, and lonidamine treatment significantly increased median survival in both models, with particularly dramatic effects in panobinostat- and marizomib-resistant cells. These data provide new insights into mechanisms of treatment resistance in gliomas.

Hypoxic niches attract and sequester tumor-associated macrophages and cytotoxic T cells and reprogram them for immunosuppression

Glioblastoma (GBM), a highly lethal brain cancer, is notorious for immunosuppression, but the mechanisms remain unclear. Here, we documented a temporospatial patterning of tumor-associated myeloid cells (TAMs) corresponding to vascular changes during GBM progression. As tumor vessels transitioned from the initial dense regular network to later scant and engorged vasculature, TAMs shifted away from perivascular regions and trafficked to vascular-poor areas. This process was heavily influenced by the immunocompetence state of the host. Utilizing a sensitive fluorescent UnaG reporter to track tumor hypoxia, coupled with single-cell transcriptomics, we revealed that hypoxic niches attracted and sequestered TAMs and cytotoxic T lymphocytes (CTLs), where they were reprogrammed toward an immunosuppressive state. Mechanistically, we identified chemokine CCL8 and cytokine IL-1β as two hypoxic-niche factors critical for TAM trafficking and co-evolution of hypoxic zones into pseudopalisading patterns. Therefore, perturbation of TAM patterning in hypoxic zones may improve tumor control.