Glioblastoma epigenome profiling identifies SOX10 as a master regulator of molecular tumour subtype

Profiling triple-negative breast cancer-specific super-enhancers identifies high-risk mesenchymal development subtype and BETi-Targetable vulnerabilities

BACKGROUND

Super-enhancers (SEs) are critical regulators of tumorigenesis and represent promising targets for bromodomain and extra-terminal domain inhibitors (BETi). However, clinical studies across various solid tumors, including triple-negative breast cancer (TNBC), have demonstrated limited BETi efficacy. This study aims to investigate SE heterogeneity in TNBC and its influence on BETi effectiveness, with the goal of advancing BETi precision treatment strategies and enhancing therapeutic efficacy.

METHODS

We conducted a comprehensive analysis of H3K27ac ChIP-Seq data from TNBC cell lines and clinical samples, integrating multiple bulk RNA-Seq, scRNA-Seq, and stRNA-Seq datasets to characterize the SE landscape and heterogeneity in TNBC. Utilizing various bioinformatics algorithms, CERES scoring, and clinical prognostic data on transcription factors (TFs), we identified core transcriptional regulatory circuits (CRCs) composed of TNBC-specific SEs and master regulators, characterizing different TNBC subtypes. The biological significance of CRCs in these different TNBC subtypes and their influence on BETi sensitivity were evaluated using in vitro and in vivo models.

RESULTS

Our findings revealed a distinct SE landscape in TNBC compared to non-TNBC and normal breast epithelium, allowing TNBC to be classified into distinct subtypes based on TNBC-specific SEs. Importantly, we identified a high-risk mesenchymal development subtype, validated across cell lines and transcriptomic analyses, primarily driven by a CRC consisting of the master regulator VAX2 and a TNBC-specific SE. This SE-VAX2 CRC is essential for sustaining the malignant traits of this subtype and increasing its sensitivity to BETi.

CONCLUSIONS

Our research clarifies the heterogeneity of SEs in TNBC and identifies a high-risk mesenchymal development subtype driven by the SE-VAX2 CRC. The subtype shows more sensitivity to BETi, supporting its precision application in TNBC.

Individualized patient tumor organoids faithfully preserve human brain tumor ecosystems and predict patient response to therapy

Tumor organoids are important tools for cancer research, but current models have drawbacks that limit their applications for predicting response to therapy. Here, we developed a fast, efficient, and complex culture system (IPTO, individualized patient tumor organoid) that accurately recapitulates the cellular and molecular pathology of human brain tumors. Patient-derived tumor explants were cultured in induced pluripotent stem cell (iPSC)-derived cerebral organoids, thus enabling culture of a wide range of human tumors in the central nervous system (CNS), including adult, pediatric, and metastatic brain cancers. Histopathological, genomic, epigenomic, and single-cell RNA sequencing (scRNA-seq) analyses demonstrated that the IPTO model recapitulates cellular heterogeneity and molecular features of original tumors. Crucially, we showed that the IPTO model predicts patient-specific drug responses, including resistance mechanisms, in a prospective patient cohort. Collectively, the IPTO model represents a major breakthrough in preclinical modeling of human cancers, which provides a path toward personalized cancer therapy.

Single-cell RNA sequencing of neonatal cortical astrocytes reveals versatile cell clusters during astrocyte-neuron conversion

BACKGROUND

Astrocytes are extensively utilized as starting cells for neuronal conversion. Our previous study discovered that a portion of primary cultured mouse neonatal cortical astrocytes can be directly converted into neurons after exposure to a neurogenic induction condition. Recent in vivo studies have demonstrated astrocyte heterogeneity in terms of their developmental origin, molecular profile, physiology, and functional outputs. We hypothesized that the heterogeneity of primary astrocytes in our study could influence their conversion potential.

METHODS AND RESULTS

We performed single-cell RNA sequencing on cells harvested at key time points during in vitro astrocyte-to-neuron conversion, specifically on Day 1 and Day 9. Through single-cell RNA sequencing analysis, we identified several subpopulations of astrocytes, labeled as Astrocyte 1 to Astrocyte 3, based on distinct gene expression patterns. Pseudotime trajectory analysis predicted the existence of three distinct cell states throughout the conversion process. Astrocyte 3 exhibited a higher propensity for neuronal conversion, with proliferation genes like Mki67 being highly expressed. Additionally, several candidate genes were identified as potentially crucial in the conversion process. Astrocyte 3 is considered a unique subtype population of astrocytes.

CONCLUSIONS

Our investigation underscores the diversity of primary neonatal cortical astrocytes and provides critical insights into the potential for astrocyte-to-neuron conversion, which may be harnessed to enhance the efficiency of this astrocyte-neuron conversion process.

Longitudinal multimodal profiling of IDH-wildtype glioblastoma reveals the molecular evolution and cellular phenotypes underlying prognostically different treatment responses

BACKGROUND

Despite recent advances in the biology of IDH-wildtype glioblastoma, it remains a devastating disease with median survival of less than 2 years. However, the molecular underpinnings of the heterogeneous response to the current standard-of-care treatment regimen consisting of maximal safe resection, adjuvant radiation, and chemotherapy with temozolomide remain unknown.

METHODS

Comprehensive histopathologic, genomic, and epigenomic evaluation of paired initial and recurrent glioblastoma specimens from 106 patients was performed to investigate the molecular evolution and cellular phenotypes underlying differential treatment responses.

RESULTS

While TERT promoter mutation and CDKN2A homozygous deletion were early events during gliomagenesis shared by initial and recurrent tumors, most other recurrent genetic alterations (eg, EGFR, PTEN, and NF1) were commonly private to initial or recurrent tumors indicating acquisition later during clonal evolution. Furthermore, glioblastomas exhibited heterogeneous epigenomic evolution with subsets becoming more globally hypermethylated, hypomethylated, or remaining stable. Glioblastoma that underwent sarcomatous transformation had shorter interval to recurrence and were significantly enriched in NF1, TP53, and RB1 alterations and the mesenchymal epigenetic class. Patients who developed somatic hypermutation following temozolomide treatment had significantly longer interval to disease recurrence and prolonged overall survival, and increased methylation at 4 specific CpG sites in the promoter region of MGMT was significantly associated with this development of hypermutation. Finally, an epigenomic evolution signature incorporating change in DNA methylation levels across 347 critical CpG sites was developed that significantly correlated with clinical outcomes.

CONCLUSIONS

Glioblastoma undergoes heterogeneous genetic, epigenetic, and cellular evolution that underlies prognostically different treatment responses.

The Human Pathology Atlas for deciphering the prognostic features of human cancers

BACKGROUND

Cancer is one of the leading causes of mortality worldwide, highlighting the urgent need for a deeper molecular understanding and the development of personalized treatments. The present study aims to establish a solid association between gene expression and patient survival outcomes to enhance the utility of the Human Pathology Atlas for cancer research.

METHODS

In this updated analysis, we examined the expression profiles of 6918 patients across 21 cancer types. We integrated data from 10 independent cancer cohorts, creating a cross-validated, reliable collection of prognostic genes. We applied systems biology approach to identify the association between gene expression profiles and patient survival outcomes. We further constructed prognostic regulatory networks for kidney renal clear cell carcinoma (KIRC) and liver hepatocellular carcinoma (LIHC), which elucidate the molecular underpinnings associated with patient survival in these cancers.

FINDINGS

We observed that gene expression during the transition from normal to tumorous tissue exhibited diverse shifting patterns in their original tissue locations. Significant correlations between gene expression and patient survival outcomes were identified in KIRC and LIHC among the major cancer types. Additionally, the prognostic regulatory network established for these two cancers showed the indicative capabilities of the Human Pathology Atlas and provides actionable insights for cancer research.

INTERPRETATION

The updated Human Pathology Atlas provides a significant foundation for precision oncology and the formulation of personalized treatment strategies. These findings deepen our understanding of cancer biology and have the potential to advance targeted therapeutic approaches in clinical practice.

FUNDING

The Knut and Alice Wallenberg Foundation (72110), the China Scholarship Council (Grant No. 202006940003).

Unravelling the mosaic: Epigenetic diversity in glioblastoma

Glioblastoma is the most common primary malignant brain tumour. Despite decades of intensive research in the disease, its prognosis remains poor, with an average survival of only 14 months after diagnosis. The remarkable level of intra- and interpatient heterogeneity is certainly contributing to the lack of progress in tackling this tumour. Epigenetic dysregulation plays an important role in glioblastoma biology and significantly contributes to intratumour heterogeneity. However, it is becoming increasingly clear that it also contributes to intertumour heterogeneity, which historically had mainly been linked to diverse genetic events occurring in different patients. In this review, we explore how DNA methylation, chromatin remodelling, microRNA (miRNA) dysregulation, and long noncoding RNA (lncRNA) alterations contribute to intertumour heterogeneity in glioblastoma, including its implications for advanced tumour stratification, which is the essential first step for developing more effective patient-specific therapeutic approaches.

Analysis of transcription profiles for the identification of master regulators as the key players in glioblastoma

Glioblastoma (GBM) is the most common malignant brain tumor with poor overall survival. Current treatment management for GBM has low efficacy, mainly due to high inter-patient heterogeneity. The transcription profiles in GBM define cell properties essential for tumor progression. We have developed an approach for the identification of master regulators (MRs) that are responsible for the gene expression changes in GBM. The approach is based on transcription factor enrichment analysis with subsequent "upstream" analysis in the signaling network. The main feature of the approach is that all calculations are performed for transcription profiles from individual samples, which allows taking into account GBM transcription heterogeneity. We identified 451 MRs that were up-regulated or down-regulated and, thus, were important parts of positive feedback loops. The number of MRs in the samples correlated with the degree of tumor immune infiltration, while the differences in MR profiles were generally consistent with the known GBM subtypes: mesenchymal, classical, and proneural. MRs densely interact with each other in the signaling network that may be associated with the robustness to pharmacological intervention. We identified 102 receptors among MRs, which is coherent with the importance of cell-cell interactions for GBM progression. The role of some of them in GBM is not currently investigated: lysophosphatidic acid receptors 5 and 6, sphingosine-1-phosphate receptor 4, lysophosphatidylserine receptors GPR34 and GPR174, and G protein-coupled receptors 84 and 132 for fatty acids. Information on the revealed MRs can be used to search for novel therapeutic strategies to treat GBM.

Towards the genomic sequence code of DNA fragility for machine learning

Genomic DNA breakages and the subsequent insertion and deletion mutations are important contributors to genome instability and linked diseases. Unlike the research in point mutations, the relationship between DNA sequence context and the propensity for strand breaks remains elusive. Here, by analyzing the differences and commonalities across myriads of genomic breakage datasets, we extract the sequence-linked rules and patterns behind DNA fragility. We show the overall deconvolution of the sequence influence into short-, mid- and long-range effects, and the stressor-dependent differences in defining the range and compositional effects on DNA fragility. We summarize and release our feature compendium as a library that can be seamlessly incorporated into genomic machine learning procedures, where DNA fragility is of concern, and train a generalized DNA fragility model on cancer-associated breakages. Structural variants (SVs) tend to stabilize regions in which they emerge, with the effect most pronounced for pathogenic SVs. In contrast, the effects of chromothripsis are seen across regions less prone to breakages. We find that viral integration may bring genome fragility, particularly for cancer-associated viruses. Overall, this work offers novel insights into the genomic sequence basis of DNA fragility and presents a powerful machine learning resource to further enhance our understanding of genome (in)stability and evolution.

SOX10 mediates glioblastoma cell-state plasticity

Phenotypic plasticity is a cause of glioblastoma therapy failure. We previously showed that suppressing the oligodendrocyte-lineage regulator SOX10 promotes glioblastoma progression. Here, we analyze SOX10-mediated phenotypic plasticity and exploit it for glioblastoma therapy design. We show that low SOX10 expression is linked to neural stem-cell (NSC)-like glioblastoma cell states and is a consequence of temozolomide treatment in animal and cell line models. Single-cell transcriptome profiling of Sox10-KD tumors indicates that Sox10 suppression is sufficient to induce tumor progression to an aggressive NSC/developmental-like phenotype, including a quiescent NSC-like cell population. The quiescent NSC state is induced by temozolomide and Sox10-KD and reduced by Notch pathway inhibition in cell line models. Combination treatment using Notch and HDAC/PI3K inhibitors extends the survival of mice carrying Sox10-KD tumors, validating our experimental therapy approach. In summary, SOX10 suppression mediates glioblastoma progression through NSC/developmental cell-state transition, including the induction of a targetable quiescent NSC state. This work provides a rationale for the design of tumor therapies based on single-cell phenotypic plasticity analysis.

The peritumoral edema index and related mechanisms influence the prognosis of GBM patients

Background

Peritumoral brain edema (PTBE) represents a characteristic phenotype of intracranial gliomas. However, there is a lack of consensus regarding the prognosis and mechanism of PTBE. In this study, clinical imaging data, along with publicly available imaging data, were utilized to assess the prognosis of PTBE in glioblastoma (GBM) patients, and the associated mechanisms were preliminarily analyzed.

Methods

We investigated relevant imaging features, including edema, in GBM patients using ITK-SNAP imaging segmentation software. Risk factors affecting progression-free survival (PFS) and overall survival (OS) were assessed using a Cox proportional hazard regression model. In addition, the impact of PTBE on PFS and OS was analyzed in clinical GBM patients using the Kaplan-Meier survival analysis method, and the results further validated by combining data from The Cancer Imaging Archive (TCIA) and The Cancer Genome Atlas (TCGA). Finally, functional enrichment analysis based on TCIA and TCGA datasets identified several pathways potentially involved in the mechanism of edema formation.

Results

This study included a total of 32 clinical GBM patients and 132 GBM patients from public databases. Univariate and multivariate analyses indicated that age and edema index (EI) are independent risk factors for PFS, but not for OS. Kaplan-Meier curves revealed consistent survival analysis results between IE groups among both clinical patients and TCIA and TCGA patients, suggesting a significant effect of PTBE on PFS but not on OS. Furthermore, functional enrichment analysis predicted the involvement of several pathways related mainly to cellular bioenergetics and vasculogenic processes in the mechanism of PTBE formation. While these novel results warrant confirmation in a larger patient cohort, they support good prognostic value for PTBE assessment in GBM.

Conclusions

Our results indicate that a low EI positively impacts disease control in GBM patients, but this does not entirely translate into an improvement in OS. Multiple genes, signaling pathways, and biological processes may contribute to the formation of peritumoral edema in GBM through cytotoxic and vascular mechanisms.

HDAC7 drives glioblastoma to a mesenchymal-like state via LGALS3-mediated crosstalk between cancer cells and macrophages

Background: Glioblastoma multiforme (GBM) is an aggressive brain tumour for which current traditional treatment approaches have been unsuccessful, owing to the high genetic heterogeneity and immunosuppressive tumour microenvironment. Methods: Single-cell and spatial transcriptomic data revealed the niche-specific enrichment of mesenchymal-like (MES-like) GBM cells and monocyte-derived macrophages (MDMs); Gain- and loss-of-function assays of HDAC7 was confirmed both in vitro and in vivo assays. Mechanistically, mass spectrum, RNA immunoprecipitation (RIP), and co-immunoprecipitation assays were conducted. Results: We found that HDAC7, which upregulated by TRIM28-mediated sumoylation at the protein levels, inhibited SOX8 expression by mediating H3K27 deacetylation. And the down-regulated SOX8 facilitated the transcriptional activity of JUN, to induce LGALS3 secretion, which then bind to the membrane protein ITGB1 on GSC and MDMs in the autocrine and paracrine manners to facilitate the transformation of the mesenchymal phenotype of GBM and the M2 polarization of MDMs, respectively. In turn, LGALS3 could also secreted by M2 MDMs to promote MES transition of GBM in a paracrine manner, creating a positive feedback loop. In translational medicine, we found that blocking LGALS3 improved the therapeutic sensitivity of HDAC inhibitors. Conclusions: Our findings revealed the role of the novel HDAC7-H3K27ac-SOX8/JUN-LGALS3-ITGB1 axis in maintaining the crosstalk between MES GBM and M2 MDM, highlighting that HDAC7 and LGALS3 may serve as potential prognostic biomarkers and therapeutic targets in GBM.

Integration of multi-omics data revealed the orphan CpG islands and enhancer-dominated c is-regulatory network in glioma

The complex transcriptional regulatory network leads to the poor prognosis of glioma. The role of orphan CpG islands (oCGIs) in the transcriptional regulatory network has been overlooked. We conducted a comprehensive exploration of the cis-regulatory roles of oCGIs and enhancers by integrating multi-omics data. Direct regulation of target genes by oCGIs or enhancers is of great importance in the cis-regulatory network. Furthermore, based on single-cell multi-omics data, we found that the highly activated cis-regulatory network in cluster 2 (C2) sustains the high proliferative potential of glioma cells. The upregulation of oCGIs and enhancers related genes in C2 results in glioma patients exhibiting resistance to radiotherapy and chemotherapy. These findings were further validated through glioma cell line related experiments. Our study offers insight into the pathogenesis of glioma and provides a strategy to treat this challenging disease.

Unveiling the Molecular Mechanisms of Glioblastoma through an Integrated Network-Based Approach

Background/Objectives: Despite current treatments extending the lifespan of Glioblastoma (GBM) patients, the average survival time is around 15-18 months, underscoring the fatality of GBM. This study aims to investigate the impact of sample heterogeneity on gene expression in GBM, identify key metabolic pathways and gene modules, and explore potential therapeutic targets. Methods: In this study, we analysed GBM transcriptome data derived from The Cancer Genome Atlas (TCGA) using genome-scale metabolic models (GEMs) and co-expression networks. We examine transcriptome data incorporating tumour purity scores (TPSs), allowing us to assess the impact of sample heterogeneity on gene expression profiles. We analysed the metabolic profile of GBM by generating condition-specific GEMs based on the TPS group. Results: Our findings revealed that over 90% of genes showing brain and glioma specificity in RNA expression demonstrate a high positive correlation, underscoring their expression is dominated by glioma cells. Conversely, negatively correlated genes are strongly associated with immune responses, indicating a complex interaction between glioma and immune pathways and non-tumorigenic cell dominance on gene expression. TPS-based metabolic profile analysis was supported by reporter metabolite analysis, highlighting several metabolic pathways, including arachidonic acid, kynurenine and NAD pathway. Through co-expression network analysis, we identified modules that significantly overlap with TPS-correlated genes. Notably, SOX11 and GSX1 are upregulated in High TPS, show a high correlation with TPS, and emerged as promising therapeutic targets. Additionally, NCAM1 exhibits a high centrality score within the co-expression module, which shows a positive correlation with TPS. Moreover, LILRB4, an immune-related gene expressed in the brain, showed a negative correlation and upregulated in Low TPS, highlighting the importance of modulating immune responses in the GBM mechanism. Conclusions: Our study uncovers sample heterogeneity's impact on gene expression and the molecular mechanisms driving GBM, and it identifies potential therapeutic targets for developing effective treatments for GBM patients.

Microenvironment T-Type calcium channels regulate neuronal and glial processes to promote glioblastoma growth

Background

Glioblastoma (GBM) is the most common primary malignant brain tumor. The aim of this study was to elucidate the role of microenvironment and intrinsic T-type calcium channels (Cav3) in regulating tumor growth and progression.

Methods

We grafted syngeneic GBM cells into Cav3.2 knockout mice to assess the role of microenvironment T-Type calcium channels on GBM tumor growth. We performed single-cell RNA-seq (scRNA-seq) of tumors from WT and Cav3.2 KO mice to elucidate the regulation of tumors by the microenvironment. We used neurons from WT and Cav3.2 KO mice in co-culture with GBM stem cells (GSC) to assess the effects of Cav3.2 on neuron/GSC synaptic connections and tumor cell growth.

Results

Cav3.2 KO in the microenvironment led to significant reduction of GBM growth and prolongation of animal survival. scRNA-seq showed that microenvironment Cav3.2 regulates neuronal and glial biological processes. Microenvironment Cav3.2 downregulated numerous genes associated with regulating the OPC cell state in GBM tumors such as SOX10 and Olig2. Neuronal Cav3.2 promoted neuron/GSC synaptic connections and GSC growth. Treatment of GSCs with the Cav3 blocker mibefradil downregulated genes associated with neuronal processes. The Cav3 blocker drug mibefradil synergized with temozolomide (TMZ) and radiation to reduce in vivo tumor growth and prolong animal survival.

Conclusions

Together these data reveal a role for microenvironment Cav3 in promoting GBM tumor progression through regulating neuronal and glial processes particularly associated with the OPC-cell state. Targeting both intrinsic and microenvironment Cav3 with the inhibitor mibefradil significantly enhanced the anti-GBM effects of TMZ and radiation.

Navigating glioblastoma complexity: the interplay of neurotransmitters and chromatin

Glioblastoma is the most aggressive brain cancer with an unfavorable prognosis for patient survival. Glioma stem cells, a subpopulation of cancer cells, drive tumor initiation, self-renewal, and resistance to therapy and, together with the microenvironment, play a crucial role in glioblastoma maintenance and progression. Neurotransmitters such as noradrenaline, dopamine, and serotonin have contrasting effects on glioblastoma development, stimulating or inhibiting its progression depending on the cellular context and through their action on glioma stem cells, perhaps changing the epigenetic landscape. Recent studies have revealed that serotonin and dopamine induce chromatin modifications related to transcriptional plasticity in the mammalian brain and possibly in glioblastoma; however, this topic still needs to be explored because of its potential implications for glioblastoma treatment. Also, it is essential to consider that neurotransmitters' effects depend on the tumor's microenvironment since it can significantly influence the response and behavior of cancer cells. This review examines the possible role of neurotransmitters as regulators of glioblastoma development, focusing on their impact on the chromatin of glioma stem cells.

p300 KAT Regulates SOX10 Stability and Function in Human Melanoma

SOX10 is a lineage-specific transcription factor critical for melanoma tumor growth; on the other hand, SOX10 loss-of-function drives the emergence of therapy-resistant, invasive melanoma phenotypes. A major challenge has been developing therapeutic strategies targeting SOX10's role in melanoma proliferation while preventing a concomitant increase in tumor cell invasion. In this study, we report that the lysine acetyltransferase (KAT) EP300 and SOX10 gene loci on chromosome 22 are frequently co-amplified in melanomas, including UV-associated and acral tumors. We further show that p300 KAT activity mediates SOX10 protein stability and that the p300 inhibitor A-485 downregulates SOX10 protein levels in melanoma cells via proteasome-mediated degradation. Additionally, A-485 potently inhibits proliferation of SOX10+ melanoma cells while decreasing invasion in AXLhigh/MITFlow melanoma cells through downregulation of metastasis-related genes. We conclude that the SOX10/p300 axis is critical to melanoma growth and invasion and that inhibition of p300 KAT activity through A-485 may be a worthwhile therapeutic approach for SOX10-reliant tumors.

SIGNIFICANCE

The p300 KAT inhibitor A-485 blocks SOX10-dependent proliferation and SOX10-independent invasion in hard-to-treat melanoma cells.

Discovery of a Therapeutic Agent for Glioblastoma Using a Systems Biology-Based Drug Repositioning Approach

Glioblastoma (GBM), a highly malignant tumour of the central nervous system, presents with a dire prognosis and low survival rates. The heterogeneous and recurrent nature of GBM renders current treatments relatively ineffective. In our study, we utilized an integrative systems biology approach to uncover the molecular mechanisms driving GBM progression and identify viable therapeutic drug targets for developing more effective GBM treatment strategies. Our integrative analysis revealed an elevated expression of CHST2 in GBM tumours, designating it as an unfavourable prognostic gene in GBM, as supported by data from two independent GBM cohorts. Further, we pinpointed WZ-4002 as a potential drug candidate to modulate CHST2 through computational drug repositioning. WZ-4002 directly targeted EGFR (ERBB1) and ERBB2, affecting their dimerization and influencing the activity of adjacent genes, including CHST2. We validated our findings by treating U-138 MG cells with WZ-4002, observing a decrease in CHST2 protein levels and a reduction in cell viability. In summary, our research suggests that the WZ-4002 drug candidate may effectively modulate CHST2 and adjacent genes, offering a promising avenue for developing efficient treatment strategies for GBM patients.

Multi-scale signaling and tumor evolution in high-grade gliomas

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

CRISPRi screen of long non-coding RNAs identifies LINC03045 regulating glioblastoma invasion

INTRODUCTION

Glioblastoma (GBM) invasion studies have focused on coding genes, while few studies evaluate long non-coding RNAs (lncRNAs), transcripts without protein-coding potential, for role in GBM invasion. We leveraged CRISPR-interference (CRISPRi) to evaluate invasive function of GBM-associated lncRNAs in an unbiased functional screen, characterizing and exploring the mechanism of identified candidates.

METHODS

We implemented a CRISPRi lncRNA loss-of-function screen evaluating association of lncRNA knockdown (KD) with invasion capacity in Matrigel. Top screen candidates were validated using CRISPRi and oligonucleotide(ASO)-mediated knockdown in three tumor lines. Clinical relevance of candidates was assessed via The Cancer Genome Atlas(TCGA) and Genotype-Tissue Expression(GTEx) survival analysis. Mediators of lncRNA effect were identified via differential expression analysis following lncRNA KD and assessed for tumor invasion using knockdown and rescue experiments.

RESULTS

Forty-eight lncRNAs were significantly associated with 33-83% decrease in invasion (p<0.01) upon knockdown. The top candidate, LINC03045, identified from effect size and p-value, demonstrated 82.7% decrease in tumor cell invasion upon knockdown, while LINC03045 expression was significantly associated with patient survival and tumor grade(p<0.0001). RNAseq analysis of LINC03045 knockdown revealed that WASF3, previously implicated in tumor invasion studies, was highly correlated with lncRNA expression, while WASF3 KD was associated with significant decrease in invasion. Finally, WASF3 overexpression demonstrated rescue of invasive function lost with LINC03045 KD.

CONCLUSION

CRISPRi screening identified LINC03045, a previously unannotated lncRNA, as critical to GBM invasion. Gene expression is significantly associated with tumor grade and survival. RNA-seq and mechanistic studies suggest that this novel lncRNA may regulate invasion via WASF3.

A prognostic neural epigenetic signature in high-grade glioma

Neural-tumor interactions drive glioma growth as evidenced in preclinical models, but clinical validation is limited. We present an epigenetically defined neural signature of glioblastoma that independently predicts patients' survival. We use reference signatures of neural cells to deconvolve tumor DNA and classify samples into low- or high-neural tumors. High-neural glioblastomas exhibit hypomethylated CpG sites and upregulation of genes associated with synaptic integration. Single-cell transcriptomic analysis reveals a high abundance of malignant stemcell-like cells in high-neural glioblastoma, primarily of the neural lineage. These cells are further classified as neural-progenitor-cell-like, astrocyte-like and oligodendrocyte-progenitor-like, alongside oligodendrocytes and excitatory neurons. In line with these findings, high-neural glioblastoma cells engender neuron-to-glioma synapse formation in vitro and in vivo and show an unfavorable survival after xenografting. In patients, a high-neural signature is associated with decreased overall and progression-free survival. High-neural tumors also exhibit increased functional connectivity in magnetencephalography and resting-state magnet resonance imaging and can be detected via DNA analytes and brain-derived neurotrophic factor in patients' plasma. The prognostic importance of the neural signature was further validated in patients diagnosed with diffuse midline glioma. Our study presents an epigenetically defined malignant neural signature in high-grade gliomas that is prognostically relevant. High-neural gliomas likely require a maximized surgical resection approach for improved outcomes.

T cell-mediated tumor killing based signature to predict the prognosis and immunotherapy for glioblastoma

Despite the significant advancements in cancer treatment brought by immune checkpoint inhibitors (ICIs), their effectiveness in treating glioblastoma (GBM) remains highly dissatisfactory. Immunotherapy relies on the fundamental concept of T cell-mediated tumor killing (TTK). Nevertheless, additional investigation is required to explore its potential in prognostic prediction and regulation of tumor microenvironment (TME) in GBM. TTK sensitivity related genes (referred to as GSTTKs) were obtained from the TISIDB. The training cohort was available from the TCGA-GBM, while the independent validation group was gathered from GEO database. Firstly, we examined differentially expressed GSTTKs (DEGs) with limma package. Afterwards, the prognostic DEGs were identified and the TTK signature was established with univariate and LASSO Cox analyses. Next, we examined the correlation between the TTK signature and outcome of GBM as well as immune phenotypes of TME. Furthermore, the evaluation of TTK signature in predicting the effectiveness of immunotherapy has also been conducted. We successfully developed a TTK signature with an independent predictive value. Patients who had a high score experienced a worse prognosis compared to patients with low scores. The TTK signature showed a strong positive association with the infiltration degree of immunocyte and the presence of various immune checkpoints. Moreover, individuals with a lower score exhibited increased responsiveness to ICIs and experienced improved prognosis. In conclusions, we successfully developed and verified a TTK signature that has the ability to predict the outcome and immune characteristics of GBM. Furthermore, the TTK signature has the potential to direct the personalized immunotherapy for GBM.

Genomic landscape of glioblastoma without IDH somatic mutation in 42 cases: a comprehensive analysis using RNA sequencing data

PURPOSE

Glioblastoma is a malignant brain tumor with a poor prognosis. Genetic mutations associated with this disease are complex are not fully understood and require further elucidation for the development of new treatments. The purpose of this study was to comprehensively analyze genetic mutations in glioblastomas and evaluate the usefulness of RNA sequencing.

PATIENTS AND METHODS

We analyzed 42 glioblastoma specimens that were resected in routine clinical practice and found wild-type variants of the IDH1 and IDH2 genes. RNA was extracted from frozen specimens and sequenced, and genetic analyses were performed using the CLC Genomics Workbench.

RESULTS

The most common genetic alterations in the 42 glioblastoma specimens were TP53 mutation (28.6%), EGFR splicing variant (16.7%), EGFR mutation (9.5%), and FGFR3 fusion (9.5%). Novel genetic mutations were detected in 8 patients (19%). In 12 cases (28.6%), driver gene mutations were not detected, suggesting an association with PPP1R14A overexpression. Our findings suggest the transcription factors SOX10 and NKX6-2 are potential markers in glioblastoma.

CONCLUSION

RNA sequencing is a promising approach for genotyping glioblastomas because it provides comprehensive information on gene expression and is relatively cost-effective.

Epigenetic regulation of tumor-immune symbiosis in glioma

Glioma is a type of aggressive and incurable brain tumor. Patients with glioma are highly resistant to all types of therapies, including immunotherapies. Epigenetic reprogramming is a key molecular hallmark in tumors across cancer types, including glioma. Mounting evidence highlights a pivotal role of epigenetic regulation in shaping tumor biology and therapeutic responses through mechanisms involving both glioma cells and immune cells, as well as their symbiotic interactions in the tumor microenvironment (TME). In this review, we discuss the molecular mechanisms of epigenetic regulation that impacts glioma cell biology and tumor immunity in both a cell-autonomous and non-cell-autonomous manner. Moreover, we provide an overview of potential therapeutic approaches that can disrupt epigenetic-regulated tumor-immune symbiosis in the glioma TME.

Utility of combining OLIG2 and SOX10 IHC expression in CNS tumours: promising biomarkers for subtyping paediatric- and adult-type gliomas

AIMS

The SOX10 transcription factor is important for the maturation of oligodendrocytes involved in central nervous system (CNS) myelination. Currently, very little information exists about its expression and potential use in CNS tumour diagnoses. The aim of our study was to characterize the expression of SOX10 in a large cohort of CNS tumours and to evaluate its potential use as a biomarker.

METHODS

We performed immunohistochemistry (IHC) for SOX10 and OLIG2 in a series of 683 cases of adult- and paediatric-type CNS tumours from different subtypes. The nuclear immunostaining results for SOX10 and OLIG2 were scored as positive (≥10% positive tumour cells) or negative.

RESULTS

OLIG2 and SOX10 were positive in diffuse midline gliomas (DMG), H3-mutant, and EZHIP-overexpressed. However, in all DMG, EGFR-mutant, SOX10 was constantly negative. In diffuse paediatric-type high-grade gliomas (HGG), all RTK1 cases were positive for both OLIG2 and SOX10. RTK2 cases were all negative for both OLIG2 and SOX10. MYCN cases variably expressed OLIG2 and were all immunonegative for SOX10. In glioblastoma, IDH-wildtype, OLIG2 was mostly positive, but SOX10 was variably expressed, depending on the epigenetic subtype. All circumscribed astrocytic gliomas were positive for both OLIG2 and SOX10 except pleomorphic xanthoastrocytomas, astroblastomas, MN1-altered, and subependymal giant cell astrocytomas. SOX10 was negative in ependymomas, meningiomas, pinealoblastomas, choroid plexus tumours, intracranial Ewing sarcomas, and embryonal tumours except neuroblastoma, FOXR2-activated.

CONCLUSION

To conclude, SOX10 can be incorporated into the IHC panel routinely used by neuropathologists in the diagnostic algorithm of embryonal tumours and for the subtyping of paediatric and adult-type HGG.

Mean global DNA methylation serves as independent prognostic marker in IDH-wildtype glioblastoma

BACKGROUND

The IDH-wildtype glioblastoma (GBM) patients have a devastating prognosis. Here, we analyzed the potential prognostic value of global DNA methylation of the tumors.

METHODS

DNA methylation of 492 primary samples and 31 relapsed samples, each treated with combination therapy, and of 148 primary samples treated with radiation alone were compared with patient survival. We determined the mean methylation values and estimated the immune cell infiltration from the methylation data. Moreover, the mean global DNA methylation of 23 GBM cell lines was profiled and correlated to their cellular radiosensitivity as measured by colony formation assay.

RESULTS

High mean DNA methylation levels correlated with improved survival, which was independent from known risk factors (MGMT promoter methylation, age, extent of resection; P = 0.009) and methylation subgroups. Notably, this correlation was also independent of immune cell infiltration, as higher number of immune cells indeed was associated with significantly better OS but lower mean methylation. Radiosensitive GBM cell lines had a significantly higher mean methylation than resistant lines (P = 0.007), and improved OS of patients treated with radiotherapy alone was also associated with higher DNA methylation (P = 0.002). Furthermore, specimens of relapsed GBM revealed a significantly lower mean DNA methylation compared to the matching primary tumor samples (P = 0.041).

CONCLUSIONS

Our results indicate that mean global DNA methylation is independently associated with outcome in glioblastoma. The data also suggest that a higher DNA methylation is associated with better radiotherapy response and less aggressive phenotype, both of which presumably contribute to the observed correlation with OS.

p300 KAT regulates SOX10 stability and function in human melanoma

SOX10 is a lineage-specific transcription factor critical for melanoma tumor growth, while SOX10 loss-of-function drives the emergence of therapy-resistant, invasive melanoma phenotypes. A major challenge has been developing therapeutic strategies targeting SOX10's role in melanoma proliferation, while preventing a concomitant increase in tumor cell invasion. Here, we report that the lysine acetyltransferase (KAT) EP300 and SOX10 gene loci on Chromosome 22 are frequently co-amplified in melanomas, including UV-associated and acral tumors. We further show that p300 KAT activity mediates SOX10 protein stability and that the p300 inhibitor, A-485, downregulates SOX10 protein levels in melanoma cells via proteasome-mediated degradation. Additionally, A-485 potently inhibits proliferation of SOX10+ melanoma cells while decreasing invasion in AXLhigh/MITFlow melanoma cells through downregulation of metastasis-related genes. We conclude that the SOX10/p300 axis is critical to melanoma growth and invasion, and that inhibition of p300 KAT activity through A-485 may be a worthwhile therapeutic approach for SOX10-reliant tumors.

Diffuse Pediatric-type High-grade Glioma Arising in an Ovarian Mature Cystic Teratoma

Immature neuroectodermal tissue can be found in the ovary as part of an immature teratoma or as part of a teratoma with malignant neuroectodermal transformation. Such lesions may closely resemble central nervous system tumors, but their biologic similarity is unclear. We describe an 18-yr-old female who presented with abdominal pain caused by an ovarian mass with widespread metastases. Histology showed a primitive, high-grade tumor arising in the background of a mature teratoma. The tumor was SOX10 positive, with focal expression of GFAP, S100, NSE, and synaptophysin. Molecular analysis demonstrated co-amplification of PDGFRA and KIT , alterations common in high-grade gliomas. By whole-genome methylation profiling, it clustered into the "diffuse pediatric-type high-grade glioma, RTK1 subtype, subclass c" group. Despite progressing through 2 lines of chemotherapy with widespread metastatic disease, she achieved an excellent response to chemotherapy directed toward aggressive germ cell tumors. This case emphasizes the importance of immunohistochemical, genomic, and epigenetic analyses to accurately classify these exceedingly rare tumors and determine the optimal therapy.

Methylation-directed regulatory networks determine enhancing and silencing of mutation disease driver genes and explain inter-patient expression variation

BACKGROUND

Common diseases manifest differentially between patients, but the genetic origin of this variation remains unclear. To explore possible involvement of gene transcriptional-variation, we produce a DNA methylation-oriented, driver-gene-wide dataset of regulatory elements in human glioblastomas and study their effect on inter-patient gene expression variation.

RESULTS

In 175 of 177 analyzed gene regulatory domains, transcriptional enhancers and silencers are intermixed. Under experimental conditions, DNA methylation induces enhancers to alter their enhancing effects or convert into silencers, while silencers are affected inversely. High-resolution mapping of the association between DNA methylation and gene expression in intact genomes reveals methylation-related regulatory units (average size = 915.1 base-pairs). Upon increased methylation of these units, their target-genes either increased or decreased in expression. Gene-enhancing and silencing units constitute cis-regulatory networks of genes. Mathematical modeling of the networks highlights indicative methylation sites, which signified the effect of key regulatory units, and add up to make the overall transcriptional effect of the network. Methylation variation in these sites effectively describe inter-patient expression variation and, compared with DNA sequence-alterations, appears as a major contributor of gene-expression variation among glioblastoma patients.

CONCLUSIONS

We describe complex cis-regulatory networks, which determine gene expression by summing the effects of positive and negative transcriptional inputs. In these networks, DNA methylation induces both enhancing and silencing effects, depending on the context. The revealed mechanism sheds light on the regulatory role of DNA methylation, explains inter-individual gene-expression variation, and opens the way for monitoring the driving forces behind deferential courses of cancer and other diseases.

Molecular signature to predict quality of life and survival with glioblastoma using Multiview omics model

Glioblastoma multiforme (GBM) patients show a variety of signs and symptoms that affect their quality of life (QOL) and self-dependence. Since most existing studies have examined prognostic factors based only on clinical factors, there is a need to consider the value of integrating multi-omics data including gene expression and proteomics with clinical data in identifying significant biomarkers for GBM prognosis. Our research aimed to isolate significant features that differentiate between short-term (≤ 6 months) and long-term (≥ 2 years) GBM survival, and between high Karnofsky performance scores (KPS ≥ 80) and low (KPS ≤ 60), using the iterative random forest (iRF) algorithm. Using the Cancer Genomic Atlas (TCGA) database, we identified 35 molecular features composed of 19 genes and 16 proteins. Our findings propose molecular signatures for predicting GBM prognosis and will improve clinical decisions, GBM management, and drug development.

Drug resistance in glioblastoma: from chemo- to immunotherapy

As the most common and aggressive type of primary brain tumor in adults, glioblastoma is estimated to end over 10,000 lives each year in the United States alone. Stand treatment for glioblastoma, including surgery followed by radiotherapy and chemotherapy (i.e., Temozolomide), has been largely unchanged since early 2000. Cancer immunotherapy has significantly shifted the paradigm of cancer management in the past decade with various degrees of success in treating many hematopoietic cancers and some solid tumors, such as melanoma and non-small cell lung cancer (NSCLC). However, little progress has been made in the field of neuro-oncology, especially in the application of immunotherapy to glioblastoma treatment. In this review, we attempted to summarize the common drug resistance mechanisms in glioblastoma from Temozolomide to immunotherapy. Our intent is not to repeat the well-known difficulty in the area of neuro-oncology, such as the blood-brain barrier, but to provide some fresh insights into the molecular mechanisms responsible for resistance by summarizing some of the most recent literature. Through this review, we also hope to share some new ideas for improving the immunotherapy outcome of glioblastoma treatment.

A systematic approach identifies p53-DREAM pathway target genes associated with blood or brain abnormalities

p53 (encoded by Trp53) is a tumor suppressor, but mouse models have revealed that increased p53 activity may cause bone marrow failure, likely through dimerization partner, RB-like, E2F4/E2F5 and MuvB (DREAM) complex-mediated gene repression. Here, we designed a systematic approach to identify p53-DREAM pathway targets, the repression of which might contribute to abnormal hematopoiesis. We used Gene Ontology analysis to study transcriptomic changes associated with bone marrow cell differentiation, then chromatin immunoprecipitation-sequencing (ChIP-seq) data to identify DREAM-bound promoters. We next created positional frequency matrices to identify evolutionary conserved sequence elements potentially bound by DREAM. The same approach was developed to find p53-DREAM targets associated with brain abnormalities, also observed in mice with increased p53 activity. Putative DREAM-binding sites were found for 151 candidate target genes, of which 106 are mutated in a blood or brain genetic disorder. Twenty-one DREAM-binding sites were tested and found to impact gene expression in luciferase assays, to notably regulate genes mutated in dyskeratosis congenita (Rtel1), Fanconi anemia (Fanca), Diamond-Blackfan anemia (Tsr2), primary microcephaly [Casc5 (or Knl1), Ncaph and Wdr62] and pontocerebellar hypoplasia (Toe1). These results provide clues on the role of the p53-DREAM pathway in regulating hematopoiesis and brain development, with implications for tumorigenesis.

NLRP6 potentiates PI3K/AKT signalling by promoting autophagic degradation of p85α to drive tumorigenesis

The PI3K/AKT pathway plays an essential role in tumour development. NOD-like receptors (NLRs) regulate innate immunity and are implicated in cancer, but whether they are involved in PI3K/AKT pathway regulation is poorly understood. Here, we report that NLRP6 potentiates the PI3K/AKT pathway by binding and destabilizing p85α, the regulatory subunit of PI3K. Mechanistically, NLRP6 recruits the E3 ligase RBX1 to p85α and ubiquitinates lysine 256 on p85α, which is recognized by the autophagy cargo receptor OPTN, causing selective autophagic degradation of p85α and subsequent activation of the PI3K/AKT pathway by reducing PTEN stability. We further show that loss of NLRP6 suppresses cell proliferation, colony formation, cell migration, and tumour growth in glioblastoma cells in vitro and in vivo. Disruption of the NLRP6/p85α interaction using the Pep9 peptide inhibits the PI3K/AKT pathway and generates potent antitumour effects. Collectively, our results suggest that NLRP6 promotes p85α degradation via selective autophagy to drive tumorigenesis, and the interaction between NLRP6 and p85α can be a promising therapeutic target for tumour treatment.

Epigenetic neural glioblastoma enhances synaptic integration and predicts therapeutic vulnerability

Neural-tumor interactions drive glioma growth as evidenced in preclinical models, but clinical validation is nascent. We present an epigenetically defined neural signature of glioblastoma that independently affects patients' survival. We use reference signatures of neural cells to deconvolve tumor DNA and classify samples into low- or high-neural tumors. High-neural glioblastomas exhibit hypomethylated CpG sites and upregulation of genes associated with synaptic integration. Single-cell transcriptomic analysis reveals high abundance of stem cell-like malignant cells classified as oligodendrocyte precursor and neural precursor cell-like in high-neural glioblastoma. High-neural glioblastoma cells engender neuron-to-glioma synapse formation in vitro and in vivo and show an unfavorable survival after xenografting. In patients, a high-neural signature associates with decreased survival as well as increased functional connectivity and can be detected via DNA analytes and brain-derived neurotrophic factor in plasma. Our study presents an epigenetically defined malignant neural signature in high-grade gliomas that is prognostically relevant.

TGF-β based risk model to predict the prognosis and immune features in glioblastoma

Background

Transforming growth factor-β (TGF-β) is a multifunctional cytokine with an important role in tissue development and tumorigenesis. TGF-β can inhibit the function of many immune cells, prevent T cells from penetrating into the tumor center, so that the tumor cells escape from immune surveillance and lead to low sensitivity to immunotherapy. However, its potential roles in predicting clinical prognosis and tumor microenvironment (TME) immune features need to be deeply investigated in glioblastoma (GBM).

Methods

The TCGA-GBM dataset was obtained from the Cancer Genome Atlas, and the validation dataset was downloaded from Gene Expression Omnibus. Firstly, differentially expressed TGF-β genes (DEGs) were screened between GBM and normal samples. Then, univariate and multivariate Cox analyses were used to identify prognostic genes and develop the TGF-β risk model. Subsequently, the roles of TGF-β risk score in predicting clinical prognosis and immune characteristics were investigated.

Results

The TGF-β risk score signature with an independent prognostic value was successfully developed. The TGF-β risk score was positively correlated with the infiltration levels of tumor-infiltrating immune cells, and the activities of anticancer immunity steps. In addition, the TGF-β risk score was positively related to the expression of immune checkpoints. Besides, the high score indicated higher sensitivity to immune checkpoint inhibitors.

Conclusions

We first developed and validated a TGF-β risk signature that could predict the clinical prognosis and TME immune features for GBM. In addition, the TGF-β signature could guide a more personalized therapeutic approach for GBM.

Applying single cell multi-omic analyses to understand treatment resistance in pediatric high grade glioma

Despite improvements in cancer patient outcomes seen in the past decade, tumor resistance to therapy remains a major impediment to achieving durable clinical responses. Intratumoral heterogeneity related to genetic, epigenetic, transcriptomic, proteomic, and metabolic differences between individual cancer cells has emerged as a driver of therapeutic resistance. This cell to cell heterogeneity can be assessed using single cell profiling technologies that enable the identification of tumor cell clones that exhibit similar defining features like specific mutations or patterns of DNA methylation. Single cell profiling of tumors before and after treatment can generate new insights into the cancer cell characteristics that confer therapeutic resistance by identifying intrinsically resistant sub-populations that survive treatment and by describing new cellular features that emerge post-treatment due to tumor cell evolution. Integrative, single cell analytical approaches have already proven advantageous in studies characterizing treatment-resistant clones in cancers where pre- and post-treatment patient samples are readily available, such as leukemia. In contrast, little is known about other cancer subtypes like pediatric high grade glioma, a class of heterogeneous, malignant brain tumors in children that rapidly develop resistance to multiple therapeutic modalities, including chemotherapy, immunotherapy, and radiation. Leveraging single cell multi-omic technologies to analyze naïve and therapy-resistant glioma may lead to the discovery of novel strategies to overcome treatment resistance in brain tumors with dismal clinical outcomes. In this review, we explore the potential for single cell multi-omic analyses to reveal mechanisms of glioma resistance to therapy and discuss opportunities to apply these approaches to improve long-term therapeutic response in pediatric high grade glioma and other brain tumors with limited treatment options.

Partial erosion on under-methylated regions and chromatin reprogramming contribute to oncogene activation in IDH mutant gliomas

BACKGROUND

IDH1/2 hotspot mutations are well known to drive oncogenic mutations in gliomas and are well-defined in the WHO 2021 classification of central nervous system tumors. Specifically, IDH mutations lead to aberrant hypermethylation of under-methylated regions (UMRs) in normal tissues through the disruption of TET enzymes. However, the chromatin reprogramming and transcriptional changes induced by IDH-related hypermethylation in gliomas remain unclear.

RESULTS

Here, we have developed a precise computational framework based on Hidden Markov Model to identify altered methylation states of UMRs at single-base resolution. By applying this framework to whole-genome bisulfite sequencing data from 75 normal brain tissues and 15 IDH mutant glioma tissues, we identified two distinct types of hypermethylated UMRs in IDH mutant gliomas. We named them partially hypermethylated UMRs (phUMRs) and fully hypermethylated UMRs (fhUMRs), respectively. We found that the phUMRs and fhUMRs exhibit distinct genomic features and chromatin states. Genes related to fhUMRs were more likely to be repressed in IDH mutant gliomas. In contrast, genes related to phUMRs were prone to be up-regulated in IDH mutant gliomas. Such activation of phUMR genes is associated with the accumulation of active H3K4me3 and the loss of H3K27me3, as well as H3K36me3 accumulation in gene bodies to maintain gene expression stability. In summary, partial erosion on UMRs was accompanied by locus-specific changes in key chromatin marks, which may contribute to oncogene activation.

CONCLUSIONS

Our study provides a computational strategy for precise decoding of methylation encroachment patterns in IDH mutant gliomas, revealing potential mechanistic insights into chromatin reprogramming that contribute to oncogenesis.

circMMD reduction following tumor treating fields inhibits glioblastoma progression through FUBP1/FIR/DVL1 and miR-15b-5p/FZD6 signaling

BACKGROUND

Tumor treating fields (TTF) is the latest treatment for GBM. Circular RNA (circRNA) has been demonstrated to play critical roles in tumorigenesis. However, the molecular mechanism of TTF remained largely unknown and the role of circRNA in TTF was not reported. The aim of this study was to elucidate the role and mechanism of circMMD in TTF treatment of GBM.

METHODS

Divergent primer was designed to verify the existence of circMMD in GBM cells. The prognostic role of circMMD was explored in glioma specimens. The knockdown and overexpressed plasmids were used to evaluate the effect of circMMD on GBM cell proliferation and TTF efficacy. RNA pull-down and RNA immunoprecipitation were performed to identify binding proteins of circMMD. Subcutaneous and intracranial tumor models were established to validate findings in vivo.

RESULTS

The expression of circMMD was elevated in GBM and its high expression indicated poor prognoses. TTF intervention could reduce circMMD synthesis, which suppressed GBM proliferation and increased TTF-mediated apoptosis. The reduction of circMMD promoted the interaction between FUBP1 and FIR, which decreased DVL1 transcription. Meanwhile, decreased circMMD would promote the activity of miR-15b-5p to degrade FZD6. Finally, the diminished expression of DVL1 and FZD6 expression suppressed the activation of Wnt/β-catenin pathway.

CONCLUSIONS

Our study revealed a novel mechanism of TTF that TTF-mediated reduction of circMMD could inhibit Wnt/β-catenin pathway to suppress GBM proliferation.

Epigenetically regulated lncRNAs dissect the intratumoural heterogeneity and facilitate immune evasion of glioblastomas

Background: Glioblastomas are the most common and malignant central nervous system (CNS) tumors that occupied a highly heterogeneous tumor microenvironment (TIME). Long noncoding RNAs (lncRNAs), whose expression can be modified by DNA methylation, are emerging as critical regulators in the immune system. However, knowledge about the epigenetic changes in lncRNAs and their contribution to the immune heterogeneity of glioma is still lacking. Methods: In this study, we integrated paired methylome and transcriptome datasets of glioblastomas and identified 2 robust immune subtypes based on lncRNA methylation features. The immune characteristics of glioma subtypes were compared. Furthermore, immune-related lncRNAs were identified and their relationships with immune evasion were evaluated. Results: Glioma immunophenotypes exhibited distinct immune-related characteristics as well as clinical and epigenetic features. 149 epigenetically regulated (ER) lncRNAs were recognized that possessed inverse variation in epigenetic and transcriptional levels between glioma subtypes. Immune-related lncRNAs were further identified through the investigation of their correlation with immune cell infiltrations and immune-related pathways. In particular, the 'Hot' glioma subtype with higher immunoactivity while a worse survival outcome was found to character immune evasion features. We finally prioritized candidate ER lncRNAs associated with immune evasion markers and response to glioma immunotherapy. Among them, CD109-AS1 and LINC02447 were validated as novel immunoevasive biomarkers for glioma through in vitro experiments. Conclusion: In summary, our study systematically reveals the crosstalk among DNA methylation, lncRNA, and immune regulation in glioblastomas, and will facilitate the development of epigenetic immunotherapy approaches.

Harnessing epithelial-mesenchymal plasticity to boost cancer immunotherapy

Immune checkpoint blockade (ICB) therapy is a powerful option for cancer treatment. Despite demonstrable progress, most patients fail to respond or achieve durable responses due to primary or acquired ICB resistance. Recently, tumor epithelial-to-mesenchymal plasticity (EMP) was identified as a critical determinant in regulating immune escape and immunotherapy resistance in cancer. In this review, we summarize the emerging role of tumor EMP in ICB resistance and the tumor-intrinsic or extrinsic mechanisms by which tumors exploit EMP to achieve immunosuppression and immune escape. We discuss strategies to modulate tumor EMP to alleviate immune resistance and to enhance the efficiency of ICB therapy. Our discussion provides new prospects to enhance the ICB response for therapeutic gain in cancer patients.

Machine learning-based identification of SOX10 as an immune regulator of macrophage in gliomas

Gliomas, originating from the glial cells, are the most lethal type of primary tumors in the central nervous system. Standard treatments like surgery have not significantly improved the prognosis of glioblastoma patients. Recently, immune therapy has become a novel and effective option. As a conserved group of transcriptional regulators, the Sry-type HMG box (SOX) family has been proved to have a correlation with numerous diseases. Based on the large-scale machine learning, we found that the SOX family, with significant immune characteristics and genomic profiles, can be divided into two distinct clusters in gliomas, among which SOX10 was identified as an excellent immune regulator of macrophage in gliomas. The high expression of SOX10 is related to a shorter OS in LGG, HGG, and pan-cancer groups but benefited from the immunotherapy. It turned out in single-cell sequencing that SOX10 is high in neurons, M1 macrophages, and neural stem cells. Also, macrophages are found to be elevated in the SOX10 high-expression group. SOX10 has a positive correlation with macrophage cytokine production and negative regulation of macrophages' chemotaxis and migration. In conclusion, our study demonstrates the outstanding cluster ability of the SOX family, indicating that SOX10 is an immune regulator of macrophage in gliomas, which can be an effective target for glioma immunotherapy.

Injury programs shape glioblastoma

Glioblastoma is the most common and aggressive primary brain cancer in adults and is almost universally fatal due to its stark therapeutic resistance. During the past decade, although survival has not substantially improved, major advances have been made in our understanding of the underlying biology. It has become clear that these devastating tumors recapitulate features of neurodevelopmental hierarchies which are influenced by the microenvironment. Emerging evidence also highlights a prominent role for injury responses in steering cellular phenotypes and contributing to tumor heterogeneity. This review highlights how the interplay between injury and neurodevelopmental programs impacts on tumor growth, invasion, and treatment resistance, and discusses potential therapeutic considerations in view of these findings.

An update on the molecular biology of glioblastoma, with clinical implications and progress in its treatment

Glioblastoma multiforme (GBM) is the most aggressive and common malignant primary brain tumor. Patients with GBM often have poor prognoses, with a median survival of ∼15 months. Enhanced understanding of the molecular biology of central nervous system tumors has led to modifications in their classifications, the most recent of which classified these tumors into new categories and made some changes in their nomenclature and grading system. This review aims to give a panoramic view of the last 3 years' findings in glioblastoma characterization, its heterogeneity, and current advances in its treatment. Several molecular parameters have been used to achieve an accurate and personalized characterization of glioblastoma in patients, including epigenetic, genetic, transcriptomic and metabolic features, as well as age- and sex-related patterns and the involvement of several noncoding RNAs in glioblastoma progression. Astrocyte-like neural stem cells and outer radial glial-like cells from the subventricular zone have been proposed as agents involved in GBM of IDH-wildtype origin, but this remains controversial. Glioblastoma metabolism is characterized by upregulation of the PI3K/Akt/mTOR signaling pathway, promotion of the glycolytic flux, maintenance of lipid storage, and other features. This metabolism also contributes to glioblastoma's resistance to conventional therapies. Tumor heterogeneity, a hallmark of GBM, has been shown to affect the genetic expression, modulation of metabolic pathways, and immune system evasion. GBM's aggressive invasion potential is modulated by cell-to-cell crosstalk within the tumor microenvironment and altered expressions of specific genes, such as ANXA2, GBP2, FN1, PHIP, and GLUT3. Nevertheless, the rising number of active clinical trials illustrates the efforts to identify new targets and drugs to treat this malignancy. Immunotherapy is still relevant for research purposes, given the amount of ongoing clinical trials based on this strategy to treat GBM, and neoantigen and nucleic acid-based vaccines are gaining importance due to their antitumoral activity by inducing the immune response. Furthermore, there are clinical trials focused on the PI3K/Akt/mTOR axis, angiogenesis, and tumor heterogeneity for developing molecular-targeted therapies against GBM. Other strategies, such as nanodelivery and computational models, may improve the drug pharmacokinetics and the prognosis of patients with GBM.

Brain cancer stem cells: resilience through adaptive plasticity and hierarchical heterogeneity

Malignant brain tumours are complex ecosystems containing neoplastic and stromal components that generate adaptive and evolutionarily driven aberrant tissues in the central nervous system. Brain cancers are cultivated by a dynamic population of stem-like cells that enforce intratumoural heterogeneity and respond to intrinsic microenvironment or therapeutically guided insults through proliferation, plasticity and restructuring of neoplastic and stromal components. Far from a rigid hierarchy, heterogeneous neoplastic populations transition between cellular states with differential self-renewal capacities, endowing them with powerful resilience. Here we review the biological machinery used by brain tumour stem cells to commandeer tissues in the intracranial space, evade immune responses and resist chemoradiotherapy. Through recent advances in single-cell sequencing, improved models to investigate the role of the tumour microenvironment and a deeper understanding of the fundamental role of the immune system in cancer biology, we are now better equipped to explore mechanisms by which these processes can be exploited for therapeutic benefit.

Methylation associated with long- or short-term survival in glioblastoma patients from the Nordic phase 3 trial

Patients with glioblastoma (GBM) have a poor outcome, but even among patients receiving the same therapies and with good prognostic factors, one can find those with exceptionally short and long survival. From the Nordic trial, which randomized GBM patients of 60 years or older between two radiotherapy arms (60 Gy or 34 Gy) or temozolomide (TMZ), we selected 59 with good prognostic factors. These selected GBM patients were equally distributed according to treatment and MGMT promoter methylation status but had long or short survival. Methylation profiling with the Illumina Infinium Methylation EPIC BeadChip arrays was performed and utilized for methylation-based CNS tumor classification, and pathway enrichment analysis of differentially methylated CpG sites (DMCs), as well as calculation of epigenetic age acceleration with three different algorithms, to compare the long and short survival groups. Samples identified by the classifier as non-GBM IDH wildtype were excluded. DMCs between long- and short-term survivors were found in patients with methylated MGMT promoter treated with TMZ (123,510), those with unmethylated MGMT treated with 60Gy radiotherapy (4,086), and with methylated MGMT promoter treated with 34Gy radiotherapy (39,649). Long-term survivors with methylated MGMT promoter treated with TMZ exhibited hypermethylation of the Wnt signaling and the platelet activation, signaling, and aggregation pathways. The joint analysis of radiotherapy arms revealed 319 DMCs between long- and short-term survivors with unmethylated MGMT and none for samples with methylated MGMT promoter. An analysis comparing epigenetic age acceleration between patients with long- and short-term survival across all treatment arms showed a decreased epigenetic age acceleration for the latter. We identified DMCs for both TMZ and RT-treated patients and epigenetic age acceleration as a potential prognostic marker, but further systematic analysis of larger patient cohorts is necessary for confirmation of their prognostic and/or predictive properties.

The overexpression of GPX8 is correlated with poor prognosis in GBM patients

Glutathione peroxidase 8 (GPX8), located in the endoplasmic reticulum, is associated with poor prognosis in several cancers. However, the expression and functions of GPX8 in cancers remain unclear. The purpose of this study was to explore the expression and functions of GPX8 in glioblastoma (GBM). We obtained expression data of GPX8 by accessing the TCGA, CGGA, GEPIA, and TIMER2.0 databases and validated them using western blot and immunohistochemistry. The Kaplan–Meier overall survival curve and Cox regression model were used to evaluate the prognostic value of GPX8 in glioma patients. Gene ontology (GO) and function enrichment analysis were used to investigate the potential function of GPX8 in GBM. Correlation analysis was used to clarify the role of GPX8 in proneural–mesenchymal transition (PMT). We studied the correlation between GPX8 expression and GBM immune infiltration by accessing cBioPortal and TIMER2.0 databases. Here, we demonstrated that GPX8 was significantly upregulated in GBM, and was associated with IDH-wildtype and mesenchymal subtype with poor prognosis. Survival analysis results indicated that GPX8 is an independent prognostic factor for overall survival (OS) in all WHO-grade glioma patients. Through the functional studies, we found that high expression of GPX8 correlated with mesenchymal signature and negatively correlated with proneural signature, indicating that GPX8 might promote PMT in GBM. Finally, based on correlation analysis, we found that the expression of GPX8 was associated with immune infiltration and the IL1/MYD88/IRAK/NF-κB pathway in GBM. Our results show that GPX8 is a key factor affecting the prognosis of GBM patients, and its targeting has the potential to provide a novel therapeutic approach.

Morphogenetic and Imaging Characteristics in Giant Cell Glioblastoma

Giant cell glioblastoma is a rare tumor entity of IDH-wildtype glioblastoma. It is usually found in the pediatric population. We describe a particular case of a female patient diagnosed histopathologically with giant cell glioblastoma, who had two recurrences in different lobes of the same cerebral hemisphere, despite positive prognostic factors and appropriate treatment. We performed an immunohistochemical characterization of giant cell glioblastoma as well as an analysis of its aggressiveness using the cytogenetic markers TP53, CDKN2A, and TP73 using the FISH technique. The clinical picture was inconsistant, the suspicion being completely different initially. Paraclinical examination and imaging initially suggested a metastasis to the insular lobe. After surgery, histopathological and immunohistochemical examinations were the basis for the diagnosis. Despite the prognostic factors known so far in the literature, the aggressiveness denoted by multiple relapses and morphogenetic tests particularizes the case and improves the literature by bringing new information about this rare neoplasm of the central nervous system.

JUN activation modulates chromatin accessibility to drive TNFα ‐induced mesenchymal transition in glioblastoma

Chromatin dynamics as well as genetic evolution underlies the adaptability of tumour cells to environmental cues. Three subtypes of tumour cells have been identified in glioblastoma, one of the commonest malignant brain tumours in adults. During tumour progression or under therapeutic pressure, the non‐mesenchymal subtypes may progress to the mesenchymal subtype, leading to unfavourable prognosis. However, the molecular mechanisms for this transition remain poorly understood. Here taking a TNFα‐induced cellular model, we profile the chromatin accessibility dynamics during mesenchymal transition. Moreover, we identify the JUN family as one of the key driving transcription factors for the gained chromatin accessibility. Accordingly, inhibition of JUN phosphorylation and therefore its transcription activity successfully impedes TNFα‐induced chromatin remodelling and mesenchymal transition. In line with these findings based on experimental models, JUN activity is positively correlated with mesenchymal features in clinical glioblastoma specimens. Together, this study unveils a deregulated transcription regulatory network in glioblastoma progression and hopefully provides a rationale for anti‐glioblastoma therapy.

Novel cancer gene discovery using a forward genetic screen in RCAS-PDGFB-driven gliomas

BACKGROUND

Malignant gliomas, the most common malignant brain tumors in adults, represent a heterogeneous group of diseases with poor prognosis. Retroviruses can cause permanent genetic alterations that modify genes close to the viral integration site.

METHODS

Here we describe the use of a high-throughput pipeline coupled to the commonly used tissue-specific retroviral RCAS-TVA mouse tumor model system. Utilizing next-generation sequencing, we show that retroviral integration sites can be reproducibly detected in malignant stem cell lines generated from RCAS-PDGFB-driven glioma biopsies.

RESULTS

A large fraction of common integration sites contained genes that have been dysregulated or misexpressed in glioma. Others overlapped with loci identified in previous glioma-related forward genetic screens, but several novel putative cancer-causing genes were also found. Integrating retroviral tagging and clinical data, Ppfibp1 was highlighted as a frequently tagged novel glioma-causing gene. Retroviral integrations into the locus resulted in Ppfibp1 upregulation, and Ppfibp1-tagged cells generated tumors with shorter latency on orthotopic transplantation. In human gliomas, increased PPFIBP1 expression was significantly linked to poor prognosis and PDGF treatment resistance.

CONCLUSIONS

Altogether, the current study has demonstrated a novel approach to tagging glioma genes via forward genetics, validating previous results, and identifying PPFIBP1 as a putative oncogene in gliomagenesis.