Meningioma transcriptomic landscape demonstrates novel subtypes with regional associated biology and patient outcome.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype.

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Ibrutinib disrupts blood-tumor barrier integrity and prolongs survival in rodent glioma model.

In malignant glioma, cytotoxic drugs are often inhibited from accessing the tumor site due to the blood-tumor barrier (BTB). Ibrutinib, FDA-approved lymphoma agent, inhibits Bruton tyrosine kinase (BTK) and has previously been shown to independently impair aortic endothelial adhesion and increase rodent glioma model survival in combination with cytotoxic therapy. Yet additional research is required to understand ibrutinib's effect on BTB function. In this study, we detail baseline BTK expression in glioma cells and its surrounding vasculature, then measure endothelial junctional expression/function changes with varied ibrutinib doses in vitro. Rat glioma cells and rodent glioma models were treated with ibrutinib alone (1-10 µM and 25 mg/kg) and in combination with doxil (10-100 µM and 3 mg/kg) to assess additive effects on viability, drug concentrations, tumor volume, endothelial junctional expression and survival. We found that ibrutinib, in a dose-dependent manner, decreased brain endothelial cell-cell adhesion over 24 h, without affecting endothelial cell viability (p < 0.005). Expression of tight junction gene and protein expression was decreased maximally 4 h after administration, along with inhibition of efflux transporter, ABCB1, activity. We demonstrated an additive effect of ibrutinib with doxil on rat glioma cells, as seen by a significant reduction in cell viability (p < 0.001) and increased CNS doxil concentration in the brain (56 ng/mL doxil alone vs. 74.6 ng/mL combination, p < 0.05). Finally, Ibrutinib, combined with doxil, prolonged median survival in rodent glioma models (27 vs. 16 days, p < 0.0001) with brain imaging showing a - 53% versus - 75% volume change with doxil alone versus combination therapy (p < 0.05). These findings indicate ibrutinib's ability to increase brain endothelial permeability via junctional disruption and efflux inhibition, to increase BTB drug entry and prolong rodent glioma model survival. Our results motivate the need to identify other BTB modifiers, all with the intent of improving survival and reducing systemic toxicities.

A road map for the treatment of pediatric diffuse midline glioma.

Recent clinical trials for H3K27-altered diffuse midline gliomas (DMGs) have shown much promise. We present a consensus roadmap and identify three major barriers: (1) refinement of experimental models to include immune and brain-specific components; (2) collaboration among researchers, clinicians, and industry to integrate patient-derived data through sharing, transparency, and regulatory considerations; and (3) streamlining clinical efforts including biopsy, CNS-drug delivery, endpoint determination, and response monitoring. We highlight the importance of comprehensive collaboration to advance the understanding, diagnostics, and therapeutics for DMGs.

Visualizing genomic characteristics across an RNA-Seq based reference landscape of normal and neoplastic brain.

In order to better understand the relationship between normal and neoplastic brain, we combined five publicly available large-scale datasets, correcting for batch effects and applying Uniform Manifold Approximation and Projection (UMAP) to RNA-Seq data. We assembled a reference Brain-UMAP including 702 adult gliomas, 802 pediatric tumors and 1409 healthy normal brain samples, which can be utilized to investigate the wealth of information obtained from combining several publicly available datasets to study a single organ site. Normal brain regions and tumor types create distinct clusters and because the landscape is generated by RNA-Seq, comparative gene expression profiles and gene ontology patterns are readily evident. To our knowledge, this is the first meta-analysis that allows for comparison of gene expression and pathways of interest across adult gliomas, pediatric brain tumors, and normal brain regions. We provide access to this resource via the open source, interactive online tool Oncoscape, where the scientific community can readily visualize clinical metadata, gene expression patterns, gene fusions, mutations, and copy number patterns for individual genes and pathway over this reference landscape.

An RNA seq-based reference landscape of human normal and neoplastic brain.

In order to better understand the relationship between normal and neoplastic brain, we combined five publicly available large-scale datasets, correcting for batch effects and applying Uniform Manifold Approximation and Projection (UMAP) to RNA-seq data. We assembled a reference Brain-UMAP including 702 adult gliomas, 802 pediatric tumors and 1409 healthy normal brain samples, which can be utilized to investigate the wealth of information obtained from combining several publicly available datasets to study a single organ site. Normal brain regions and tumor types create distinct clusters and because the landscape is generated by RNA seq, comparative gene expression profiles and gene ontology patterns are readily evident. To our knowledge, this is the first meta-analysis that allows for comparison of gene expression and pathways of interest across adult gliomas, pediatric brain tumors, and normal brain regions. We provide access to this resource via the open source, interactive online tool Oncoscape, where the scientific community can readily visualize clinical metadata, gene expression patterns, gene fusions, mutations, and copy number patterns for individual genes and pathway over this reference landscape.

An RNA seq-based reference landscape of human normal and neoplastic brain.

In order to better understand the relationship between normal and neoplastic brain, we combined five publicly available large-scale datasets, correcting for batch effects and applying Uniform Manifold Approximation and Projection (UMAP) to RNA-seq data. We assembled a reference Brain-UMAP including 702 adult gliomas, 802 pediatric tumors and 1409 healthy normal brain samples, which can be utilized to investigate the wealth of information obtained from combining several publicly available datasets to study a single organ site. Normal brain regions and tumor types create distinct clusters and because the landscape is generated by RNA seq, comparative gene expression profiles and gene ontology patterns are readily evident. To our knowledge, this is the first meta-analysis that allows for comparison of gene expression and pathways of interest across adult gliomas, pediatric brain tumors, and normal brain regions. We provide access to this resource via the open source, interactive online tool Oncoscape, where the scientific community can readily visualize clinical metadata, gene expression patterns, gene fusions, mutations, and copy number patterns for individual genes and pathway over this reference landscape.

Expanded analysis of high-grade astrocytoma with piloid features identifies an epigenetically and clinically distinct subtype associated with neurofibromatosis type 1.

High-grade astrocytoma with piloid features (HGAP) is a recently recognized glioma type whose classification is dependent on its global epigenetic signature. HGAP is characterized by alterations in the mitogen-activated protein kinase (MAPK) pathway, often co-occurring with CDKN2A/B homozygous deletion and/or ATRX mutation. Experience with HGAP is limited and to better understand this tumor type, we evaluated an expanded cohort of patients (n = 144) with these tumors, as defined by DNA methylation array testing, with a subset additionally evaluated by next-generation sequencing (NGS). Among evaluable cases, we confirmed the high prevalence CDKN2A/B homozygous deletion, and/or ATRX mutations/loss in this tumor type, along with a subset showing NF1 alterations. Five of 93 (5.4%) cases sequenced harbored TP53 mutations and RNA fusion analysis identified a single tumor containing an NTRK2 gene fusion, neither of which have been previously reported in HGAP. Clustering analysis revealed the presence of three distinct HGAP subtypes (or groups = g) based on whole-genome DNA methylation patterns, which we provisionally designated as gNF1 (n = 18), g1 (n = 72), and g2 (n = 54) (median ages 43.5 years, 47 years, and 32 years, respectively). Subtype gNF1 is notable for enrichment with patients with Neurofibromatosis Type 1 (33.3%, p = 0.0008), confinement to the posterior fossa, hypermethylation in the NF1 enhancer region, a trend towards decreased progression-free survival (p = 0.0579), RNA processing pathway dysregulation, and elevated non-neoplastic glia and neuron cell content (p < 0.0001 and p < 0.0001, respectively). Overall, our expanded cohort broadens the genetic, epigenetic, and clinical phenotype of HGAP and provides evidence for distinct epigenetic subtypes in this tumor type.

Oncogenic role of a developmentally regulated NTRK2 splice variant.

Temporally regulated alternative splicing choices are vital for proper development, yet the wrong splice choice may be detrimental. Here, we highlight a previously unidentified role for the neurotrophin receptor splice variant TrkB.T1 in neurodevelopment, embryogenesis, transformation, and oncogenesis across multiple tumor types in humans and mice. TrkB.T1 is the predominant NTRK2 isoform across embryonic organogenesis, and forced overexpression of this embryonic pattern causes multiple solid and nonsolid tumors in mice in the context of tumor suppressor loss. TrkB.T1 also emerges as the predominant NTRK isoform expressed in a wide range of adult and pediatric tumors, including those harboring tropomyosin receptor kinase fusions. Affinity purification-mass spectrometry proteomic analysis reveals distinct interactors with known developmental and oncogenic signaling pathways such as Wnt, transforming growth factor-ß, Sonic Hedgehog, and Ras. From alterations in splicing factors to changes in gene expression, the discovery of isoform specific oncogenes with embryonic ancestry has the potential to shape the way we think about developmental systems and oncology.

Machine learning modeling of genome-wide copy number alteration signatures reliably predicts IDH mutational status in adult diffuse glioma.

Knowledge of 1p/19q-codeletion and IDH1/2 mutational status is necessary to interpret any investigational study of diffuse gliomas in the modern era. While DNA sequencing is the gold standard for determining IDH mutational status, genome-wide methylation arrays and gene expression profiling have been used for surrogate mutational determination. Previous studies by our group suggest that 1p/19q-codeletion and IDH mutational status can be predicted by genome-wide somatic copy number alteration (SCNA) data alone, however a rigorous model to accomplish this task has yet to be established. In this study, we used SCNA data from 786 adult diffuse gliomas in The Cancer Genome Atlas (TCGA) to develop a two-stage classification system that identifies 1p/19q-codeleted oligodendrogliomas and predicts the IDH mutational status of astrocytic tumors using a machine-learning model. Cross-validated results on TCGA SCNA data showed near perfect classification results. Furthermore, our astrocytic IDH mutation model validated well on four additional datasets (AUC = 0.97, AUC = 0.99, AUC = 0.95, AUC = 0.96) as did our 1p/19q-codeleted oligodendroglioma screen on the two datasets that contained oligodendrogliomas (MCC = 0.97, MCC = 0.97). We then retrained our system using data from these validation sets and applied our system to a cohort of REMBRANDT study subjects for whom SCNA data, but not IDH mutational status, is available. Overall, using genome-wide SCNAs, we successfully developed a system to robustly predict 1p/19q-codeletion and IDH mutational status in diffuse gliomas. This system can assign molecular subtype labels to tumor samples of retrospective diffuse glioma cohorts that lack 1p/19q-codeletion and IDH mutational status, such as the REMBRANDT study, recasting these datasets as validation cohorts for diffuse glioma research.

Radiogenomic modeling predicts survival-associated prognostic groups in glioblastoma.

BACKGROUND: Combined whole-exome sequencing (WES) and somatic copy number alteration (SCNA) information can separate isocitrate dehydrogenase (IDH)1/2-wildtype glioblastoma into two prognostic molecular subtypes, which cannot be distinguished by epigenetic or clinical features. The potential for radiographic features to discriminate between these molecular subtypes has yet to be established. METHODS: Radiologic features (n = 35 340) were extracted from 46 multisequence, pre-operative magnetic resonance imaging (MRI) scans of IDH1/2-wildtype glioblastoma patients from The Cancer Imaging Archive (TCIA), all of whom have corresponding WES/SCNA data. We developed a novel feature selection method that leverages the structure of extracted MRI features to mitigate the dimensionality challenge posed by the disparity between a large number of features and the limited patients in our cohort. Six traditional machine learning classifiers were trained to distinguish molecular subtypes using our feature selection method, which was compared to least absolute shrinkage and selection operator (LASSO) feature selection, recursive feature elimination, and variance thresholding. RESULTS: We were able to classify glioblastomas into two prognostic subgroups with a cross-validated area under the curve score of 0.80 (±0.03) using ridge logistic regression on the 15-dimensional principle component analysis (PCA) embedding of the features selected by our novel feature selection method. An interrogation of the selected features suggested that features describing contours in the T2 signal abnormality region on the T2-weighted fluid-attenuated inversion recovery (FLAIR) MRI sequence may best distinguish these two groups from one another. CONCLUSIONS: We successfully trained a machine learning model that allows for relevant targeted feature extraction from standard MRI to accurately predict molecularly-defined risk-stratifying IDH1/2-wildtype glioblastoma patient groups.