Molecular analysis of archival diagnostic prostate cancer biopsies identifies genomic similarities in cases with progression post-radiotherapy, and those with de novo metastatic disease.

BACKGROUND: It is important to identify molecular features that improve prostate cancer (PCa) risk stratification before radical treatment with curative intent. Molecular analysis of historical diagnostic formalin-fixed paraffin-embedded (FFPE) prostate biopsies from cohorts with post-radiotherapy (RT) long-term clinical follow-up has been limited. Utilizing parallel sequencing modalities, we performed a proof-of-principle sequencing analysis of historical diagnostic FFPE prostate biopsies. We compared patients with (i) stable PCa (sPCa) postprimary or salvage RT, (ii) progressing PCa (pPCa) post-RT, and (iii) de novo metastatic PCa (mPCa). METHODS: A cohort of 19 patients with diagnostic prostate biopsies (n = 6 sPCa, n = 5 pPCa, n = 8 mPCa) and mean 4 years 10 months follow-up (diagnosed 2009-2016) underwent nucleic acid extraction from demarcated malignancy. Samples underwent 3'RNA sequencing (3'RNAseq) (n = 19), nanoString analysis (n = 12), and Illumina 850k methylation (n = 8) sequencing. Bioinformatic analysis was performed to coherently identify differentially expressed genes and methylated genomic regions (MGRs). RESULTS: Eighteen of 19 samples provided useable 3'RNAseq data. Principal component analysis (PCA) demonstrated similar expression profiles between pPCa and mPCa cases, versus sPCa. Coherently differentially methylated probes between these groups identified ~600 differentially MGRs. The top 50 genes with increased expression in pPCa patients were associated with reduced progression-free survival post-RT (p < 0.0001) in an external cohort. CONCLUSIONS: 3'RNAseq, nanoString and 850k-methylation analyses are each achievable from historical FFPE diagnostic pretreatment prostate biopsies, unlocking the potential to utilize large cohorts of historic clinical samples. Profiling similarities between individuals with pPCa and mPCa suggests biological similarities and historical radiological staging limitations, which warrant further investigation.

Metabolic symbiosis between oxygenated and hypoxic tumour cells: An agent-based modelling study.

Deregulated metabolism is one of the hallmarks of cancer. It is well-known that tumour cells tend to metabolize glucose via glycolysis even when oxygen is available and mitochondrial respiration is functional. However, the lower energy efficiency of aerobic glycolysis with respect to mitochondrial respiration makes this behaviour, namely the Warburg effect, counter-intuitive, although it has now been recognized as source of anabolic precursors. On the other hand, there is evidence that oxygenated tumour cells could be fuelled by exogenous lactate produced from glycolysis. We employed a multi-scale approach that integrates multi-agent modelling, diffusion-reaction, stoichiometric equations, and Boolean networks to study metabolic cooperation between hypoxic and oxygenated cells exposed to varying oxygen, nutrient, and inhibitor concentrations. The results show that the cooperation reduces the depletion of environmental glucose, resulting in an overall advantage of using aerobic glycolysis. In addition, the oxygen level was found to be decreased by symbiosis, promoting a further shift towards anaerobic glycolysis. However, the oxygenated and hypoxic populations may gradually reach quasi-equilibrium. A sensitivity analysis using Latin hypercube sampling and partial rank correlation shows that the symbiotic dynamics depends on properties of the specific cell such as the minimum glucose level needed for glycolysis. Our results suggest that strategies that block glucose transporters may be more effective to reduce tumour growth than those blocking lactate intake transporters.

New role of fat-free mass in cancer risk linked with genetic predisposition.

Cancer risk is associated with the widely debated measure body mass index (BMI). Fat mass and fat-free mass measurements from bioelectrical impedance may further clarify this association. The UK Biobank is a rare resource in which bioelectrical impedance and BMI data was collected on ~ 500,000 individuals. Using this dataset, a comprehensive analysis using regression, principal component and genome-wide genetic association, provided multiple levels of evidence that increasing whole body fat (WBFM) and fat-free mass (WBFFM) are both associated with increased post-menopausal breast cancer risk, and colorectal cancer risk in men. WBFM was inversely associated with prostate cancer. We also identified rs615029[T] and rs1485995[G] as associated in independent analyses with both PMBC (p = 1.56E-17 and 1.78E-11) and WBFFM (p = 2.88E-08 and 8.24E-12), highlighting splice variants of the intriguing long non-coding RNA CUPID1 (LINC01488) as a potential link between PMBC risk and fat-free mass.

Intra-prostatic tumour evolution, steps in metastatic spread and histogenomic associations revealed by integration of multi-region whole-genome sequencing with histopathological features.

BACKGROUND: Extension of prostate cancer beyond the primary site by local invasion or nodal metastasis is associated with poor prognosis. Despite significant research on tumour evolution in prostate cancer metastasis, the emergence and evolution of cancer clones at this early stage of expansion and spread are poorly understood. We aimed to delineate the routes of evolution and cancer spread within the prostate and to seminal vesicles and lymph nodes, linking these to histological features that are used in diagnostic risk stratification. METHODS: We performed whole-genome sequencing on 42 prostate cancer samples from the prostate, seminal vesicles and lymph nodes of five treatment-naive patients with locally advanced disease. We spatially mapped the clonal composition of cancer across the prostate and the routes of spread of cancer cells within the prostate and to seminal vesicles and lymph nodes in each individual by analysing a total of > 19,000 copy number corrected single nucleotide variants. RESULTS: In each patient, we identified sample locations corresponding to the earliest part of the malignancy. In patient 10, we mapped the spread of cancer from the apex of the prostate to the seminal vesicles and identified specific genomic changes associated with the transformation of adenocarcinoma to amphicrine morphology during this spread. Furthermore, we show that the lymph node metastases in this patient arose from specific cancer clones found at the base of the prostate and the seminal vesicles. In patient 15, we observed increased mutational burden, altered mutational signatures and histological changes associated with whole genome duplication. In all patients in whom histological heterogeneity was observed (4/5), we found that the distinct morphologies were located on separate branches of their respective evolutionary trees. CONCLUSIONS: Our results link histological transformation with specific genomic alterations and phylogenetic branching. These findings have implications for diagnosis and risk stratification, in addition to providing a rationale for further studies to characterise the genetic changes causally linked to morphological transformation. Our study demonstrates the value of integrating multi-region sequencing with histopathological data to understand tumour evolution and identify mechanisms of prostate cancer spread.

Robustness and reproducibility for AI learning in biomedical sciences: RENOIR.

Artificial intelligence (AI) techniques are increasingly applied across various domains, favoured by the growing acquisition and public availability of large, complex datasets. Despite this trend, AI publications often suffer from lack of reproducibility and poor generalisation of findings, undermining scientific value and contributing to global research waste. To address these issues and focusing on the learning aspect of the AI field, we present RENOIR (REpeated random sampliNg fOr machIne leaRning), a modular open-source platform for robust and reproducible machine learning (ML) analysis. RENOIR adopts standardised pipelines for model training and testing, introducing elements of novelty, such as the dependence of the performance of the algorithm on the sample size. Additionally, RENOIR offers automated generation of transparent and usable reports, aiming to enhance the quality and reproducibility of AI studies. To demonstrate the versatility of our tool, we applied it to benchmark datasets from health, computer science, and STEM (Science, Technology, Engineering, and Mathematics) domains. Furthermore, we showcase RENOIR's successful application in recently published studies, where it identified classifiers for SET2D and TP53 mutation status in cancer. Finally, we present a use case where RENOIR was employed to address a significant pharmacological challenge-predicting drug efficacy. RENOIR is freely available at https://github.com/alebarberis/renoir .

Antitumour effect of the mitochondrial complex III inhibitor Atovaquone in combination with anti-PD-L1 therapy in mouse cancer models.

Immune checkpoint blockade (ICB) provides effective and durable responses for several tumour types by unleashing an immune response directed against cancer cells. However, a substantial number of patients treated with ICB develop relapse or do not respond, which has been partly attributed to the immune-suppressive effect of tumour hypoxia. We have previously demonstrated that the mitochondrial complex III inhibitor atovaquone alleviates tumour hypoxia both in human xenografts and in cancer patients by decreasing oxygen consumption and consequently increasing oxygen availability in the tumour. Here, we show that atovaquone alleviates hypoxia and synergises with the ICB antibody anti-PD-L1, significantly improving the rates of tumour eradication in the syngeneic CT26 model of colorectal cancer. The synergistic effect between atovaquone and anti-PD-L1 relied on CD8+ T cells, resulted in the establishment of a tumour-specific memory immune response, and was not associated with any toxicity. We also tested atovaquone in combination with anti-PD-L1 in the LLC (lung) and MC38 (colorectal) cancer syngeneic models but, despite causing a considerable reduction in tumour hypoxia, atovaquone did not add any therapeutic benefit to ICB in these models. These results suggest that atovaquone has the potential to improve the outcomes of patients treated with ICB, but predictive biomarkers are required to identify individuals likely to benefit from this intervention.

Developing and Validating a Multivariable Prognostic-Predictive Classifier for Treatment Escalation of Oropharyngeal Squamous Cell Carcinoma: The PREDICTR-OPC Study.

PURPOSE: While there are several prognostic classifiers, to date, there are no validated predictive models that inform treatment selection for oropharyngeal squamous cell carcinoma (OPSCC).Our aim was to develop clinical and/or biomarker predictive models for patient outcome and treatment escalation for OPSCC. EXPERIMENTAL DESIGN: We retrospectively collated clinical data and samples from a consecutive cohort of OPSCC cases treated with curative intent at ten secondary care centers in United Kingdom and Poland between 1999 and 2012. We constructed tissue microarrays, which were stained and scored for 10 biomarkers. We then undertook multivariable regression of eight clinical parameters and 10 biomarkers on a development cohort of 600 patients. Models were validated on an independent, retrospectively collected, 385-patient cohort. RESULTS: A total of 985 subjects (median follow-up 5.03 years, range: 4.73-5.21 years) were included. The final biomarker classifier, comprising p16 and survivin immunohistochemistry, high-risk human papillomavirus (HPV) DNA in situ hybridization, and tumor-infiltrating lymphocytes, predicted benefit from combined surgery + adjuvant chemo/radiotherapy over primary chemoradiotherapy in the high-risk group [3-year overall survival (OS) 63.1% vs. 41.1%, respectively, HR = 0.32; 95% confidence interval (CI), 0.16-0.65; P = 0.002], but not in the low-risk group (HR = 0.4; 95% CI, 0.14-1.24; P = 0.114). On further adjustment by propensity scores, the adjusted HR in the high-risk group was 0.34, 95% CI = 0.17-0.67, P = 0.002, and in the low-risk group HR was 0.5, 95% CI = 0.1-2.38, P = 0.384. The concordance index was 0.73. CONCLUSIONS: We have developed a prognostic classifier, which also appears to demonstrate moderate predictive ability. External validation in a prospective setting is now underway to confirm this and prepare for clinical adoption.

A machine learning and directed network optimization approach to uncover TP53 regulatory patterns.

TP53, the Guardian of the Genome, is the most frequently mutated gene in human cancers and the functional characterization of its regulation is fundamental. To address this we employ two strategies: machine learning to predict the mutation status of TP53from transcriptomic data, and directed regulatory networks to reconstruct the effect of mutations on the transcipt levels of TP53 targets. Using data from established databases (Cancer Cell Line Encyclopedia, The Cancer Genome Atlas), machine learning could predict the mutation status, but not resolve different mutations. On the contrary, directed network optimization allowed to infer the TP53 regulatory profile across: (1) mutations, (2) irradiation in lung cancer, and (3) hypoxia in breast cancer, and we could observe differential regulatory profiles dictated by (1) mutation type, (2) deleterious consequences of the mutation, (3) known hotspots, (4) protein changes, (5) stress condition (irradiation/hypoxia). This is an important first step toward using regulatory networks for the characterization of the functional consequences of mutations, and could be extended to other perturbations, with implications for drug design and precision medicine.

MAGNETO: Cell type marker panel generator from single-cell transcriptomic data.

Single-cell RNA sequencing experiments produce data useful to identify different cell types, including uncharacterized and rare ones. This enables us to study the specific functional roles of these cells in different microenvironments and contexts. After identifying a (novel) cell type of interest, it is essential to build succinct marker panels, composed of a few genes referring to cell surface proteins and clusters of differentiation molecules, able to discriminate the desired cells from the other cell populations. In this work, we propose a fully-automatic framework called MAGNETO, which can help construct optimal marker panels starting from a single-cell gene expression matrix and a cell type identity for each cell. MAGNETO builds effective marker panels solving a tailored bi-objective optimization problem, where the first objective regards the identification of the genes able to isolate a specific cell type, while the second conflicting objective concerns the minimization of the total number of genes included in the panel. Our results on three public datasets show that MAGNETO can identify marker panels that identify the cell populations of interest better than state-of-the-art approaches. Finally, by fine-tuning MAGNETO, our results demonstrate that it is possible to obtain marker panels with different specificity levels.

A role for SETD2 loss in tumorigenesis through DNA methylation dysregulation.

SETD2-dependent H3 Lysine-36 trimethylation (H3K36me3) has been recently linked to the deposition of de-novo DNA methylation. SETD2 is frequently mutated in cancer, however, the functional impact of SETD2 loss and depletion on DNA methylation across cancer types and tumorigenesis is currently unknown. Here, we perform a pan-cancer analysis and show that both SETD2 mutation and reduced expression are associated with DNA methylation dysregulation across 21 out of the 24 cancer types tested. In renal cancer, these DNA methylation changes are associated with altered gene expression of oncogenes, tumour suppressors, and genes involved in neoplasm invasiveness, including TP53, FOXO1, and CDK4. This suggests a new role for SETD2 loss in tumorigenesis and cancer aggressiveness through DNA methylation dysregulation. Moreover, using a robust machine learning methodology, we develop and validate a 3-CpG methylation signature which is sufficient to predict SETD2 mutation status with high accuracy and correlates with patient prognosis.

The impact of survivorship bias in glioblastoma research.

Despite advances in the therapy of Central Nervous System (CNS) malignancies, treatment of glioblastoma (GB) poses significant challenges due to GB resistance and high recurrence rates following post-operative radio-chemotherapy. The majority of prognostic and predictive GB biomarkers are currently developed using tumour samples obtained through surgical interventions. However, the selection criteria adopted by different neurosurgeons to determine which cases are suitable for surgery make operated patients not representative of all GB cases. Particularly, geriatric and frail individuals are excluded from surgical consideration in some cancer centers. Such selection generates a survival (or selection) bias that introduces limitations, rendering the patients or data chosen for downstream analyses not representative of the entire community. In this review, we discuss the implication of survivorship bias on current and novel biomarkers for patient selection, stratification, therapy, and outcome analyses.

Homologous recombination suppresses transgenerational DNA end resection and chromosomal instability in fission yeast.

Chromosomal instability (CIN) drives cell-to-cell heterogeneity, and the development of genetic diseases, including cancer. Impaired homologous recombination (HR) has been implicated as a major driver of CIN, however, the underlying mechanism remains unclear. Using a fission yeast model system, we establish a common role for HR genes in suppressing DNA double-strand break (DSB)-induced CIN. Further, we show that an unrepaired single-ended DSB arising from failed HR repair or telomere loss is a potent driver of widespread CIN. Inherited chromosomes carrying a single-ended DSB are subject to cycles of DNA replication and extensive end-processing across successive cell divisions. These cycles are enabled by Cullin 3-mediated Chk1 loss and checkpoint adaptation. Subsequent propagation of unstable chromosomes carrying a single-ended DSB continues until transgenerational end-resection leads to fold-back inversion of single-stranded centromeric repeats and to stable chromosomal rearrangements, typically isochromosomes, or to chromosomal loss. These findings reveal a mechanism by which HR genes suppress CIN and how DNA breaks that persist through mitotic divisions propagate cell-to-cell heterogeneity in the resultant progeny.

WNT5A-ROR2 axis mediates VEGF dependence of BRAF mutant melanoma.

PURPOSE: Despite recent advances, approximately 50% of patient with metastatic melanoma eventually succumb to the disease. Patients with melanomas harboring a BRAF mutation (BRAFMut) have a worse prognosis than those with wildtype (BRAFWT) tumors. Unexpectedly, interim AVAST-M Phase III trial data reported benefit from adjuvant anti-VEGF bevacizumab only in the BRAFMut group. We sought to find mechanisms underpinning this sensitivity. METHODS: We investigated this finding in vitro and in vivo using melanoma cell lines and clones generated by BRAFV600E knock-in on a BRAFWT background. RESULTS: Compared with BRAFWT cells, isogenic BRAFV600E clones secreted more VEGF and exhibited accelerated growth rates as spheroids and xenografts, which were more vascular and proliferative. Recapitulating AVAST-M findings, bevacizumab affected only BRAFV600E xenografts, inducing significant tumor growth delay, reduced vascularity and increased necrosis. We identified 814 differentially expressed genes in isogenic BRAFV600E/BRAFWT clones. Of 61 genes concordantly deregulated in clinical melanomas ROR2 was one of the most upregulated by BRAFV600E. ROR2 was shown to be RAF-MEK regulated in BRAFV600E cells and its depletion suppressed VEGF secretion down to BRAFWT levels. The ROR2 ligand WNT5A was also overexpressed in BRAFMut melanomas, and in ROR2-overexpressing BRAFV600E cells MEK inhibition downregulated WNT5A and VEGF secretion. CONCLUSIONS: These data implicate WNT5A-ROR2 in VEGF secretion, vascularity, adverse outcomes and bevacizumab sensitivity of BRAFMut melanomas, suggesting that this axis has potential therapeutic relevance.

Machine Learning Using Gene-Sets to Infer miRNA Function.

miRNA are regulators of cell phenotype, and there is clear evidence that these small posttranscriptional modifiers of gene expression are involved in defining a cellular response across states of development and disease. Classical methods for elucidating the repressive effect of a miRNA on its targets involve controlling for the many factors influencing miRNA action, and this can be achieved in cell lines, but misses tissue and organism level context which are key to a miRNA function. Also, current technology to carry out this validation is limited in both generalizability and throughput. Methodologies with greater scalability and rapidity are required to better understand the function of these important species of RNA. To this end, there is an increasing store of RNA expression level data incorporating both miRNA and mRNA, and in this chapter, we describe how to use machine learning and gene-sets to translate the knowledge of phenotype defined by mRNA to putative roles for miRNA. We outline our approach to this process and highlight how it was done for our miRNA annotation of the hallmarks of cancer using the Cancer Genome Atlas (TCGA) dataset. The concepts we present are applicable across datasets and phenotypes, and we highlight potential pitfalls and challenges that may be faced as they are used.

Single cell RNA-sequencing: A powerful yet still challenging technology to study cellular heterogeneity.

Almost all biomedical research to date has relied upon mean measurements from cell populations, however it is well established that what it is observed at this macroscopic level can be the result of many interactions of several different single cells. Thus, the observable macroscopic 'average' cannot outright be used as representative of the 'average cell'. Rather, it is the resulting emerging behaviour of the actions and interactions of many different cells. Single-cell RNA sequencing (scRNA-Seq) enables the comparison of the transcriptomes of individual cells. This provides high-resolution maps of the dynamic cellular programmes allowing us to answer fundamental biological questions on their function and evolution. It also allows to address medical questions such as the role of rare cell populations contributing to disease progression and therapeutic resistance. Furthermore, it provides an understanding of context-specific dependencies, namely the behaviour and function that a cell has in a specific context, which can be crucial to understand some complex diseases, such as diabetes, cardiovascular disease and cancer. Here, we provide an overview of scRNA-Seq, including a comparative review of emerging technologies and computational pipelines. We discuss the current and emerging applications and focus on tumour heterogeneity a clear example of how scRNA-Seq can provide new understanding of a complex disease. Additionally, we review the limitations and highlight the need of powerful computational pipelines and reproducible protocols for the broader acceptance of this technique in basic and clinical research.

New perspectives in liquid biopsy for glioma patients.

PURPOSE OF REVIEW: Gliomas are the most common primary tumors of the central nervous system. They are characterized by a disappointing prognosis and ineffective therapy that has shown no substantial improvements in the past 20 years. The lack of progress in treating gliomas is linked with the inadequacy of suitable tumor samples to plan translational studies and support laboratory developments. To overcome the use of tumor tissue, this commentary review aims to highlight the potential for the clinical application of liquid biopsy (intended as the study of circulating biomarkers in the blood), focusing on circulating tumor cells, circulating DNA and circulating noncoding RNA. RECENT FINDINGS: Thanks to the increasing sensitivity of sequencing techniques, it is now possible to analyze circulating nucleic acids and tumor cells (liquid biopsy). SUMMARY: Although studies on the use of liquid biopsy are still at an early stage, the potential clinical applications of liquid biopsy in the study of primary brain cancer are many and have the potential to revolutionize the approach to neuro-oncology, and importantly, they offer the possibility of gathering information on the disease at any time during its history.

Obesity: a perfect storm for carcinogenesis.

Obesity-related cancers account for 40% of the cancer cases observed in the USA and obesity is overtaking smoking as the most widespread modifiable risk factor for carcinogenesis. Here, we use the hallmarks of cancer framework to delineate how obesity might influence the carcinogenic hallmarks in somatic cells. We discuss the effects of obesity on (a) sustaining proliferative signaling; (b) evading growth suppressors; (c) resisting cell death; (d) enabling replicative immortality; (e) inducing angiogenesis; (f) activating invasion and metastasis; (g) reprogramming energy metabolism; and (h) avoiding immune destruction, together with its effects on genome instability and tumour-promoting inflammation. We present the current understanding and controversies in this evolving field, and highlight some areas in need of further cross-disciplinary focus. For instance, the relative importance of the many potentially causative obesity-related factors is unclear for each type of malignancy. Even within a single tumour type, it is currently unknown whether one obesity-related factor consistently plays a predominant role, or if this varies between patients or, even in a single patient with time. Clarifying how the hallmarks are affected by obesity may lead to novel prevention and treatment strategies for the increasingly obese population.

Endothelial cell death after ionizing radiation does not impair vascular structure in mouse tumor models.

The effect of radiation therapy on tumor vasculature has long been a subject of debate. Increased oxygenation and perfusion have been documented during radiation therapy. Conversely, apoptosis of endothelial cells in irradiated tumors has been proposed as a major contributor to tumor control. To examine these contradictions, we use multiphoton microscopy in two murine tumor models: MC38, a highly vascularized, and B16F10, a moderately vascularized model, grown in transgenic mice with tdTomato-labeled endothelium before and after a single (15 Gy) or fractionated (5 × 3 Gy) dose of radiation. Unexpectedly, even these high doses lead to little structural change of the perfused vasculature. Conversely, non-perfused vessels and blind ends are substantially impaired after radiation accompanied by apoptosis and reduced proliferation of their endothelium. RNAseq analysis of tumor endothelial cells confirms the modification of gene expression in apoptotic and cell cycle regulation pathways after irradiation. Therefore, we conclude that apoptosis of tumor endothelial cells after radiation does not impair vascular structure.

The MITF regulatory network in melanoma.

Bidirectional interactions between plastic tumor cells and the microenvironment critically impact tumor evolution and metastatic dissemination by enabling cancer cells to adapt to microenvironmental stresses by switching phenotype. In melanoma, a key determinant of phenotypic identity is the microphthalmia-associated transcription factor MITF that promotes proliferation, suppresses senescence, and anticorrelates with immune infiltration and therapy resistance. What determines whether MITF can activate or repress genes associated with specific phenotypes, or how signaling regulating MITF might impact immune infiltration is poorly understood. Here, we find that MITF binding to genes associated with high MITF is via classical E/M-box motifs, but genes downregulated when MITF is high contain FOS/JUN/AP1/ATF3 sites. Significantly, the repertoire of MITF-interacting factors identified here includes JUN and ATF3 as well as many previously unidentified interactors. As high AP1 activity is a hallmark of MITFLow , invasive, slow-cycling, therapy resistant cells, the ability of MITF to repress AP1-regulated genes provides an insight into how MITF establishes and maintains a pro-proliferative phenotype. Moreover, although ß-catenin has been linked to immune exclusion, many Hallmark ß-catenin signaling genes are associated with immune infiltration. Instead, low MITF together with Notch signaling is linked to immune infiltration in both mouse and human melanoma tumors.

Liver glycogen phosphorylase is upregulated in glioblastoma and provides a metabolic vulnerability to high dose radiation.

Channelling of glucose via glycogen, known as the glycogen shunt, may play an important role in the metabolism of brain tumours, especially in hypoxic conditions. We aimed to dissect the role of glycogen degradation in glioblastoma (GBM) response to ionising radiation (IR). Knockdown of the glycogen phosphorylase liver isoform (PYGL), but not the brain isoform (PYGB), decreased clonogenic growth and survival of GBM cell lines and sensitised them to IR doses of 10-12 Gy. Two to five days after IR exposure of PYGL knockdown GBM cells, mitotic catastrophy and a giant multinucleated cell morphology with senescence-like phenotype developed. The basal levels of the lysosomal enzyme alpha-acid glucosidase (GAA), essential for autolysosomal glycogen degradation, and the lipidated forms of gamma-aminobutyric acid receptor-associated protein-like (GABARAPL1 and GABARAPL2) increased in shPYGL U87MG cells, suggesting a compensatory mechanism of glycogen degradation. In response to IR, dysregulation of autophagy was shown by accumulation of the p62 and the lipidated form of GABARAPL1 and GABARAPL2 in shPYGL U87MG cells. IR increased the mitochondrial mass and the colocalisation of mitochondria with lysosomes in shPYGL cells, thereby indicating reduced mitophagy. These changes coincided with increased phosphorylation of AMP-activated protein kinase and acetyl-CoA carboxylase 2, slower ATP generation in response to glucose loading and progressive loss of oxidative phosphorylation. The resulting metabolic deficiencies affected the availability of ATP required for mitosis, resulting in the mitotic catastrophy observed in shPYGL cells following IR. PYGL mRNA and protein levels were higher in human GBM than in normal human brain tissues and high PYGL mRNA expression in GBM correlated with poor patient survival. In conclusion, we show a major new role for glycogen metabolism in GBM cancer. Inhibition of glycogen degradation sensitises GBM cells to high-dose IR indicating that PYGL is a potential novel target for the treatment of GBMs.