ABSTRACT
Therapy resistance represents an unmet challenge in the treatment of medulloblastoma. Accordingly, the identification of targets that mark drug-resistant cell populations, or drive the proliferation of resistant cells, may improve treatment strategies. To address this, we undertook a targeted approach focused on the multi-functional transcription factor YB-1. Genetic knockdown of YB-1 in Group 3 medulloblastoma cell lines diminished cell invasion in 3D in vitro assays and increased sensitivity to standard-of-care chemotherapeutic vincristine and anti-cancer agents panobinostat and JQ1. For vincristine, this occurred in part by YB-1-mediated transcriptional regulation of multi-drug resistance gene ABCB1, as determined by chromatin immunoprecipitation. Whole transcriptome sequencing of YB-1 knockdown cells identified a role for YB-1 in the regulation of tumourigenic processes, including lipid metabolism, cell death and survival and MYC and mTOR pathways. Stable cisplatin- and vincristine-tolerant Group 3 and SHH cell lines were generated to identify additional mechanisms driving resistance to standard-of-care medulloblastoma therapy. Next-generation sequencing revealed a vastly different transcriptomic landscape following chronic drug exposure, including a drug-tolerant seven-gene expression signature, common to all sequenced drug-tolerant cell lines, representing therapeutically targetable genes implicated in the acquisition of drug tolerance. Our findings provide significant insight into mechanisms and genes underlying therapy resistance in medulloblastoma.
ABSTRACT
The most common malignant brain tumour in children, medulloblastoma (MB), is subdivided into four clinically relevant molecular subgroups, although targeted therapy options informed by understanding of different cellular features are lacking. Here, by comparing the most aggressive subgroup (Group 3) with the intermediate (SHH) subgroup, we identify crucial differences in tumour heterogeneity, including unique metabolism-driven subpopulations in Group 3 and matrix-producing subpopulations in SHH. To analyse tumour heterogeneity, we profiled individual tumour nodules at the cellular level in 3D MB hydrogel models, which recapitulate subgroup specific phenotypes, by single cell RNA sequencing (scRNAseq) and 3D OrbiTrap Secondary Ion Mass Spectrometry (3D OrbiSIMS) imaging. In addition to identifying known metabolites characteristic of MB, we observed intra- and internodular heterogeneity and identified subgroup-specific tumour subpopulations. We showed that extracellular matrix factors and adhesion pathways defined unique SHH subpopulations, and made up a distinct shell-like structure of sulphur-containing species, comprising a combination of small leucine-rich proteoglycans (SLRPs) including the collagen organiser lumican. In contrast, the Group 3 tumour model was characterized by multiple subpopulations with greatly enhanced oxidative phosphorylation and tricarboxylic acid (TCA) cycle activity. Extensive TCA cycle metabolite measurements revealed very high levels of succinate and fumarate with malate levels almost undetectable particularly in Group 3 tumour models. In patients, high fumarate levels (NMR spectroscopy) alongside activated stress response pathways and high Nuclear Factor Erythroid 2-Related Factor 2 (NRF2; gene expression analyses) were associated with poorer survival. Based on these findings we predicted and confirmed that NRF2 inhibition increased sensitivity to vincristine in a long-term 3D drug treatment assay of Group 3 MB. Thus, by combining scRNAseq and 3D OrbiSIMS in a relevant model system we were able to define MB subgroup heterogeneity at the single cell level and elucidate new druggable biomarkers for aggressive Group 3 and low-risk SHH MB.
