Single-cell RNA sequencing and spatial transcriptomics reveal cancer-associated fibroblasts in glioblastoma with protumoral effects.

Cancer-associated fibroblasts (CAFs) were presumed absent in glioblastoma given the lack of brain fibroblasts. Serial trypsinization of glioblastoma specimens yielded cells with CAF morphology and single-cell transcriptomic profiles based on their lack of copy number variations (CNVs) and elevated individual cell CAF probability scores derived from the expression of 9 CAF markers and absence of 5 markers from non-CAF stromal cells sharing features with CAFs. Cells without CNVs and with high CAF probability scores were identified in single-cell RNA-Seq of 12 patient glioblastomas. Pseudotime reconstruction revealed that immature CAFs evolved into subtypes, with mature CAFs expressing actin alpha 2, smooth muscle (ACTA2). Spatial transcriptomics from 16 patient glioblastomas confirmed CAF proximity to mesenchymal glioblastoma stem cells (GSCs), endothelial cells, and M2 macrophages. CAFs were chemotactically attracted to GSCs, and CAFs enriched GSCs. We created a resource of inferred crosstalk by mapping expression of receptors to their cognate ligands, identifying PDGF and TGF-ß as mediators of GSC effects on CAFs and osteopontin and HGF as mediators of CAF-induced GSC enrichment. CAFs induced M2 macrophage polarization by producing the extra domain A (EDA) fibronectin variant that binds macrophage TLR4. Supplementing GSC-derived xenografts with CAFs enhanced in vivo tumor growth. These findings are among the first to identify glioblastoma CAFs and their GSC interactions, making them an intriguing target.

The role of sex genotype in paediatric CNS tumour incidence and survival.

PURPOSE: Evidence exists, in CNS germinomas and medulloblastomas (MB), that patient sex significantly influences incidence and outcome. The role of sex genotype in other paediatric CNS tumours remains unclear. This study sought to examine the role of sex genotype in CNS tumour incidence and overall survival (OS). METHODS: Age-adjusted incidence and OS rates were collected from the Surveillance Epidemiology and End Result (SEER) registry between 2000 and 2011 for common paediatric (<=19 years) CNS tumours: pilocytic astrocytoma (PA), anaplastic astrocytoma, glioblastoma (GBM), medulloblastoma, supratentorial CNS embryonal tumour, ependymoma, and germinoma. All patients with histologically confirmed, ICD-03 coded, first tumours, were included. Kaplan-Meier and Cox regression analyses were used to calculate hazard ratios (HR). RESULTS: The total cases are as follows: males=3018 and females=2276. Highest incidence was seen in PA (n=2103). GBM displayed the worst OS, whilst PA displayed the best. Higher incidence was observed in males for all tumours, except PA. Females with ependymoma had significantly better OS compared to males, whereas males with germinomas had better OS compared to females. Females <1 year with AA had better OS than males. Increasing age significantly improved male and female survival in ependymoma and medulloblastoma. CONCLUSION: Interrogating population-based registries such as SEER minimises bias and provides credible data. Observed differences in incidence and OS between the sexes for different paediatric CNS tumours provide useful prognostic information for clinicians. Sex genotype was a significant independent prognostic factor in ependymomas and germinomas. Further investigation of possible epigenetic and hormonal differences may provide sex-specific vulnerabilities that may be exploitable for targeted therapy.

Untargeted adductomics of newborn dried blood spots identifies modifications to human serum albumin associated with childhood leukemia.

The developing fetus is exposed to chemicals, which are metabolized to electrophiles that form adducts with nucleophilic Cys34 of human serum albumin (HSA). By measuring these adducts in neonatal blood spots (NBS), we obtain information regarding fetal exposures during the last month of gestation. To discover potential risk factors for childhood leukemia resulting from in utero exposures, we used untargeted adductomics to measure HSA-Cys34 adducts in 782 archived NBS, collected from incident cases of childhood acute lymphoblastic leukemia (ALL) or acute myeloid leukemia (AML) and matched population-based controls. Among a total of 28 Cys34 modifications that were measured, we found no differences in adduct abundances between childhood leukemia cases and controls overall. However, cases of T-cell ALL had higher abundances of adducts of reactive carbonyl species and a Cys34 disulfide of homocysteine was present at lower levels in AML cases. These results suggest that oxidative stress and lipid peroxidation may be etiologic factors of T-cell ALL, and alterations in one-carbon metabolism and epigenetic changes may be predictors of AML. Future replication of the results with larger sample sizes is necessary.

Filtering procedures for untargeted LC-MS metabolomics data.

BACKGROUND: Untargeted metabolomics datasets contain large proportions of uninformative features that can impede subsequent statistical analysis such as biomarker discovery and metabolic pathway analysis. Thus, there is a need for versatile and data-adaptive methods for filtering data prior to investigating the underlying biological phenomena. Here, we propose a data-adaptive pipeline for filtering metabolomics data that are generated by liquid chromatography-mass spectrometry (LC-MS) platforms. Our data-adaptive pipeline includes novel methods for filtering features based on blank samples, proportions of missing values, and estimated intra-class correlation coefficients. RESULTS: Using metabolomics datasets that were generated in our laboratory from samples of human blood, as well as two public LC-MS datasets, we compared our data-adaptive filtering method with traditional methods that rely on non-method specific thresholds. The data-adaptive approach outperformed traditional approaches in terms of removing noisy features and retaining high quality, biologically informative ones. The R code for running the data-adaptive filtering method is provided at https://github.com/courtneyschiffman/Metabolomics-Filtering . CONCLUSIONS: Our proposed data-adaptive filtering pipeline is intuitive and effectively removes uninformative features from untargeted metabolomics datasets. It is particularly relevant for interrogation of biological phenomena in data derived from complex matrices associated with biospecimens.

Proteomic Analysis Implicates Vimentin in Glioblastoma Cell Migration.

We previously showed lithium chloride (LiCl) and other inhibitors of glycogen synthase kinase-3 (GSK-3) including 6-bromo-indirubin-3-oxime (BIO), can block glioblastoma (GBM) cell migration. To investigate the mechanisms involved we used two-dimensional difference in-gel electrophoresis (2D-DIGE) and mass spectrometry to identify proteins altered after treatment of U251 GBM cells with 20 mM LiCl. Downregulation of the intermediate filament protein vimentin was the most significant change identified. Analysis of patient tumor samples revealed that vimentin is expressed abundantly in GBM, and is prognostic especially in lower grade tumors. Additionally, siRNA-mediated vimentin knockdown impaired GBM migration. Western blotting showed that treatment with LiCl or small molecule GSK-3 inhibitors led to the rapid downregulation of detergent soluble vimentin levels across a panel of GBM-derived cells. Fluorescence reactivation after photobleaching (FRAP) microscopy studies showed a significant reduction in the ability of the vimentin cytoskeleton to recover from photo-bleaching in the presence of LiCl or BIO. Biochemical studies revealed that GSK-3 and vimentin directly interact, and analysis of vimentin revealed a GSK-3 consensus phosphorylation site. We conclude that anti-migratory compounds with the ability to inhibit GSK-3 have effects on vimentin cytoskeletal dynamics, which may play a role in their anti-invasive activity.

The functional synergism of microRNA clustering provides therapeutically relevant epigenetic interference in glioblastoma.

MicroRNA deregulation is a consistent feature of glioblastoma, yet the biological effect of each single gene is generally modest, and therapeutically negligible. Here we describe a module of microRNAs, constituted by miR-124, miR-128 and miR-137, which are co-expressed during neuronal differentiation and simultaneously lost in gliomagenesis. Each one of these miRs targets several transcriptional regulators, including the oncogenic chromatin repressors EZH2, BMI1 and LSD1, which are functionally interdependent and involved in glioblastoma recurrence after therapeutic chemoradiation. Synchronizing the expression of these three microRNAs in a gene therapy approach displays significant anticancer synergism, abrogates this epigenetic-mediated, multi-protein tumor survival mechanism and results in a 5-fold increase in survival when combined with chemotherapy in murine glioblastoma models. These transgenic microRNA clusters display intercellular propagation in vivo, via extracellular vesicles, extending their biological effect throughout the whole tumor. Our results support the rationale and feasibility of combinatorial microRNA strategies for anticancer therapies.

MiR-19b non-canonical binding is directed by HuR and confers chemosensitivity through regulation of P-glycoprotein in breast cancer.

MicroRNAs and RNA-binding proteins exert regulation on >60% of coding genes, yet interplay between them is little studied. Canonical microRNA binding occurs by base-pairing of microRNA 3'-ends to complementary "seed regions" in mRNA 3'UTRs, resulting in translational repression. Similarly, regulatory RNA-binding proteins bind to mRNAs, modifying stability or translation. We investigated post-transcriptional regulation acting on the xenobiotic pump ABCB1/P-glycoprotein, which is implicated in cancer therapy resistance. We characterised the ABCB1 UTRs in primary breast cancer cells and identified UTR sequences that responded to miR-19b despite lacking a canonical binding site. Sequences did, however, contain consensus sites for the RNA-binding protein HuR. We demonstrated that a tripartite complex of HuR, miR-19b and UTR directs repression of ABCB1/P-glycoprotein expression, with HuR essential for non-canonical miR-19b binding thereby controlling chemosensitivity of breast cancer cells. This exemplifies a new cooperative model between RNA-binding proteins and microRNAs to expand the repertoire of mRNAs that can be regulated. This study suggests a novel therapeutic target to impair P-glycoprotein mediated drug efflux, and also indicates that current microRNA binding predictions that rely on seed regions alone may be too conservative.

Disruption of the ß1L Isoform of GABP Reverses Glioblastoma Replicative Immortality in a TERT Promoter Mutation-Dependent Manner.

TERT promoter mutations reactivate telomerase, allowing for indefinite telomere maintenance and enabling cellular immortalization. These mutations specifically recruit the multimeric ETS factor GABP, which can form two functionally independent transcription factor species: a dimer or a tetramer. We show that genetic disruption of GABPß1L (ß1L), a tetramer-forming isoform of GABP that is dispensable for normal development, results in TERT silencing in a TERT promoter mutation-dependent manner. Reducing TERT expression by disrupting ß1L culminates in telomere loss and cell death exclusively in TERT promoter mutant cells. Orthotopic xenografting of ß1L-reduced, TERT promoter mutant glioblastoma cells rendered lower tumor burden and longer overall survival in mice. These results highlight the critical role of GABPß1L in enabling immortality in TERT promoter mutant glioblastoma.

Immune evasion mediated by PD-L1 on glioblastoma-derived extracellular vesicles.

Binding of programmed death ligand-1 (PD-L1) to programmed cell death protein-1 (PD1) leads to cancer immune evasion via inhibition of T cell function. One of the defining characteristics of glioblastoma, a universally fatal brain cancer, is its profound local and systemic immunosuppression. Glioblastoma has also been shown to generate extracellular vesicles (EVs), which may play an important role in tumor progression. We thus hypothesized that glioblastoma EVs may be important mediators of immunosuppression and that PD-L1 could play a role. We show that glioblastoma EVs block T cell activation and proliferation in response to T cell receptor stimulation. PD-L1 was expressed on the surface of some, but not of all, glioblastoma-derived EVs, with the potential to directly bind to PD1. An anti-PD1 receptor blocking antibody significantly reversed the EV-mediated blockade of T cell activation but only when PD-L1 was present on EVs. When glioblastoma PD-L1 was up-regulated by IFN-γ, EVs also showed some PD-L1-dependent inhibition of T cell activation. PD-L1 expression correlated with the mesenchymal transcriptome profile and was anatomically localized in the perinecrotic and pseudopalisading niche of human glioblastoma specimens. PD-L1 DNA was present in circulating EVs from glioblastoma patients where it correlated with tumor volumes of up to 60 cm3. These results suggest that PD-L1 on EVs may be another mechanism for glioblastoma to suppress antitumor immunity and support the potential of EVs as biomarkers in tumor patients.

Genomic analysis of the origins and evolution of multicentric diffuse lower-grade gliomas.

Background: Rare multicentric lower-grade gliomas (LGGs) represent a unique opportunity to study the heterogeneity among distinct tumor foci in a single patient and to infer their origins and parallel patterns of evolution. Methods: In this study, we integrate clinical features, histology, and immunohistochemistry for 4 patients with multicentric LGG, arising both synchronously and metachronously. For 3 patients we analyze the phylogeny of the lesions using exome sequencing, including one case with a total of 8 samples from the 2 lesions. Results: One patient was diagnosed with multicentric isocitrate dehydrogenase 1 (IDH1) mutated diffuse astrocytomas harboring distinct IDH1 mutations, R132H and R132C; the latter mutation has been associated with Li-Fraumeni syndrome, which was subsequently confirmed in the patient's germline DNA and shown in additional cases with The Cancer Genome Atlas data. In another patient, phylogenetic analysis of synchronously arising grade II and grade III diffuse astrocytomas demonstrated a single shared mutation, IDH1 R132H, and revealed convergent evolution via non-overlapping mutations in ATRX and TP53. In 2 cases, there was divergent evolution of IDH1-mutated and 1p/19q-codeleted oligodendroglioma and IDH1-mutated and 1p/19q-intact diffuse astrocytoma, occurring synchronously in one case and metachronously in a second. Conclusions: Each tumor in multicentric LGG cases may arise independently or may diverge very early in their development, presenting as genetically and histologically distinct tumors. Comprehensive sampling of these lesions can therefore significantly alter diagnosis and management. Additionally, somatic IDH1 R132C mutation in either multicentric or solitary LGG identifies unsuspected germline TP53 mutation, validating the limited number of published cases.

Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1.

IDH1 mutation is the earliest genetic alteration in low-grade gliomas (LGGs), but its role in tumor recurrence is unclear. Mutant IDH1 drives overproduction of the oncometabolite d-2-hydroxyglutarate (2HG) and a CpG island (CGI) hypermethylation phenotype (G-CIMP). To investigate the role of mutant IDH1 at recurrence, we performed a longitudinal analysis of 50 IDH1 mutant LGGs. We discovered six cases with copy number alterations (CNAs) at the IDH1 locus at recurrence. Deletion or amplification of IDH1 was followed by clonal expansion and recurrence at a higher grade. Successful cultures derived from IDH1 mutant, but not IDH1 wild type, gliomas systematically deleted IDH1 in vitro and in vivo, further suggestive of selection against the heterozygous mutant state as tumors progress. Tumors and cultures with IDH1 CNA had decreased 2HG, maintenance of G-CIMP, and DNA methylation reprogramming outside CGI. Thus, while IDH1 mutation initiates gliomagenesis, in some patients mutant IDH1 and 2HG are not required for later clonal expansions.

GLUT3 upregulation promotes metabolic reprogramming associated with antiangiogenic therapy resistance.

Clinical trials revealed limited response duration of glioblastomas to VEGF-neutralizing antibody bevacizumab. Thriving in the devascularized microenvironment occurring after antiangiogenic therapy requires tumor cell adaptation to decreased glucose, with 50% less glucose identified in bevacizumab-treated xenografts. Compared with bevacizumab-responsive xenograft cells, resistant cells exhibited increased glucose uptake, glycolysis, 13C NMR pyruvate to lactate conversion, and survival in low glucose. Glucose transporter 3 (GLUT3) was upregulated in bevacizumab-resistant versus sensitive xenografts and patient specimens in a HIF-1α-dependent manner. Resistant versus sensitive cell mitochondria in oxidative phosphorylation-selective conditions produced less ATP. Despite unchanged mitochondrial numbers, normoxic resistant cells had lower mitochondrial membrane potential than sensitive cells, confirming poorer mitochondrial health, but avoided the mitochondrial dysfunction of hypoxic sensitive cells. Thin-layer chromatography revealed increased triglycerides in bevacizumab-resistant versus sensitive xenografts, a change driven by mitochondrial stress. A glycogen synthase kinase-3ß inhibitor suppressing GLUT3 transcription caused greater cell death in bevacizumab-resistant than -responsive cells. Overexpressing GLUT3 in tumor cells recapitulated bevacizumab-resistant cell features: survival and proliferation in low glucose, increased glycolysis, impaired oxidative phosphorylation, and rapid in vivo proliferation only slowed by bevacizumab to that of untreated bevacizumab-responsive tumors. Targeting GLUT3 or the increased glycolysis reliance in resistant tumors could unlock the potential of antiangiogenic treatments.

A validated microRNA profile with predictive potential in glioblastoma patients treated with bevacizumab.

PURPOSE: We investigated whether microRNA expression data from glioblastoma could be used to produce a profile that defines a bevacizumab responsive group of patients. PATIENTS AND METHODS: TCGA microRNA expression data from tumors resected at first diagnosis of glioblastoma in patients treated with bevacizumab at any time during the course of their disease were randomly separated into training (n = 50) and test (n = 37) groups for model generation. MicroRNA-seq data for 51 patients whose treatment included bevacizumab in the BELOB trial were used as an independent validation cohort. RESULTS: Using penalized regression we identified 8 microRNAs as potential predictors of overall survival in the training set. We dichotomized the response score based on the most prognostic minimum of a density plot of the response scores (log-rank HR = 0.16, p = 1.2e(-5)) and validated the profile in the test cohort (one-sided log-rank HR = 0.34, p = 0.026). Analysis of the profile using all samples in the TCGA glioblastoma dataset, regardless of treatment received, (n = 473) showed that the prediction of patient benefit was not significant (HR = 0.84, p = 0.083) suggesting the profile is specific to bevacizumab. Further independent validation of our microRNA profile in RNA-seq data from patients treated with bevacizumab (alone or in combination with CCNU) at glioblastoma recurrence in the BELOB trial confirmed that our microRNA profile predicted patient benefit from bevacizumab (HR = 0.59, p = 0.043). CONCLUSION: We have identified and validated an 8-microRNA profile that predicts overall survival in patients with glioblastoma treated with bevacizumab. This may be useful for identifying patients who are likely to benefit from this agent.

Prognostic microRNAs in high-grade glioma reveal a link to oligodendrocyte precursor differentiation.

MicroRNA expression can be exploited to define tumor prognosis and stratification for precision medicine. It remains unclear whether prognostic microRNA signatures are exclusively tumor grade and/or molecular subtype-specific, or whether common signatures of aggressive clinical behavior can be identified. Here, we defined microRNAs that are associated with good and poor prognosis in grade III and IV gliomas using data from The Cancer Genome Atlas. Pathway analysis of microRNA targets that are differentially expressed in good and poor prognosis glioma identified a link to oligodendrocyte development. Notably, a microRNA expression profile that is characteristic of a specific oligodendrocyte precursor cell type (OP1) correlates with microRNA expression from 597 of these tumors and is consistently associated with poor patient outcome in grade III and IV gliomas. Our study reveals grade-independent and subtype-independent prognostic molecular signatures in high-grade glioma and provides a framework for investigating the mechanisms of brain tumor aggressiveness.

Prediction of clinical outcome in glioblastoma using a biologically relevant nine-microRNA signature.

BACKGROUND: Glioblastoma is the most aggressive primary brain tumor, and is associated with a very poor prognosis. In this study we investigated the potential of microRNA expression profiles to predict survival in this challenging disease. METHODS: MicroRNA and mRNA expression data from glioblastoma (n = 475) and grade II and III glioma (n = 178) were accessed from The Cancer Genome Atlas. LASSO regression models were used to identify a prognostic microRNA signature. Functionally relevant targets of microRNAs were determined using microRNA target prediction, experimental validation and correlation of microRNA and mRNA expression data. RESULTS: A 9-microRNA prognostic signature was identified which stratified patients into risk groups strongly associated with survival (p = 2.26e-09), significant in all glioblastoma subtypes except the non-G-CIMP proneural group. The statistical significance of the microRNA signature was higher than MGMT methylation in temozolomide treated tumors. The 9-microRNA risk score was validated in an independent dataset (p = 4.50e-02) and also stratified patients into high- and low-risk groups in lower grade glioma (p = 5.20e-03). The majority of the 9 microRNAs have been previously linked to glioblastoma biology or treatment response. Integration of the expression patterns of predicted microRNA targets revealed a number of relevant microRNA/target pairs, which were validated in cell lines. CONCLUSIONS: We have identified a novel, biologically relevant microRNA signature that stratifies high- and low-risk patients in glioblastoma. MicroRNA/mRNA interactions identified within the signature point to novel regulatory networks. This is the first study to formulate a survival risk score for glioblastoma which consists of microRNAs associated with glioblastoma biology and/or treatment response, indicating a functionally relevant signature.

MicroRNAs in cancer: biomarkers, functions and therapy.

The emergence of microRNAs has been one of the defining developments in cancer biology over the past decade, and the explosion of knowledge in this area has brought forward new diagnostic and therapeutic opportunities. The importance of microRNAs in cancer has been underlined by the identification of alterations in microRNA target binding sites and the microRNA processing machinery in tumor cells. Clinical trials utilizing microRNA profiling for patient prognosis and clinical response are now underway, and the first microRNA mimic entered the clinic for cancer therapy in 2013. In this article we review the potential applications of microRNAs for the clinical assessment of patient outcome in cancer, as well as in cancer monitoring and therapy.

microRNA-148a is a prognostic oncomiR that targets MIG6 and BIM to regulate EGFR and apoptosis in glioblastoma.

Great interest persists in useful prognostic and therapeutic targets in glioblastoma. In this study, we report the definition of miRNA (miR)-148a as a novel prognostic oncomiR in glioblastoma. miR-148a expression was elevated in human glioblastoma specimens, cell lines, and stem cells (GSC) compared with normal human brain and astrocytes. High levels were a risk indicator for glioblastoma patient survival. Functionally, miR-148a expression increased cell growth, survival, migration, and invasion in glioblastoma cells and GSCs and promoted GSC neurosphere formation. Two direct targets of miR-148a were identified, the EGF receptor (EGFR) regulator MIG6 and the apoptosis regulator BIM, which rescue experiments showed were essential to mediate the oncogenic activity of miR-148a. By inhibiting MIG6 expression, miR-148a reduced EGFR trafficking to Rab7-expressing compartments, which includes late endosomes and lysosomes. This process coincided with reduced degradation and elevated expression and activation of EGFR. Finally, inhibition of miR-148a strongly suppressed GSC and glioblastoma xenograft growth in vivo. Taken together, our findings provide a comprehensive analysis of the prognostic value and oncogenic function of miR-148a in glioblastoma, further defining it as a potential target for glioblastoma therapy.

MicroRNA-128 coordinately targets Polycomb Repressor Complexes in glioma stem cells.

BACKGROUND: The Polycomb Repressor Complex (PRC) is an epigenetic regulator of transcription whose action is mediated by 2 protein complexes, PRC1 and PRC2. PRC is oncogenic in glioblastoma, where it is involved in cancer stem cell maintenance and radioresistance. METHODS: We used a set of glioblastoma patient samples, glioma stem cells, and neural stem cells from a mouse model of glioblastoma. We characterized gene/protein expression and cellular phenotypes by quantitative PCR/Western blotting and clonogenic, cell-cycle, and DNA damage assays. We performed overexpression/knockdown studies by lentiviral infection and microRNA/small interfering RNA oligonucleotide transfection. RESULTS: We show that microRNA-128 (miR-128) directly targets mRNA of SUZ12, a key component of PRC2, in addition to BMI1, a component of PRC1 that we previously showed as a target as well. This blocks the partially redundant functions of PRC1/PRC2, thereby significantly reducing PRC activity and its associated histone modifications. MiR-128 and SUZ12/BMI1 show opposite expression in human glioblastomas versus normal brain and in glioma stemlike versus neural stem cells. Furthermore, miR-128 renders glioma stemlike cells less radioresistant by preventing the radiation-induced expression of both PRC components. Finally, miR-128 expression is significantly reduced in neural stem cells from the brain of young, presymptomatic mice in our mouse model of glioblastoma. This suggests that loss of miR-128 expression in brain is an early event in gliomagenesis. Moreover, knockdown of miR-128 expression in nonmalignant mouse and human neural stem cells led to elevated expression of PRC components and increased clonogenicity. CONCLUSIONS: MiR-128 is an important suppressor of PRC activity, and its absence is an early event in gliomagenesis.