MRI-localized biopsies reveal subtype-specific differences in molecular and cellular composition at the margins of glioblastoma

Glioblastomas (GBMs) diffusely infiltrate the brain, making complete removal by surgical resection impossible. The mixture of neoplastic and nonneoplastic cells that remain after surgery form the biological context for adjuvant therapeutic intervention and recurrence. We performed RNA-sequencing (RNA-seq) and histological analysis on radiographically guided biopsies taken from different regions of GBM and showed that the tissue contained within the contrast-enhancing (CE) core of tumors have different cellular and molecular compositions compared with tissue from the nonenhancing (NE) margins of tumors. Comparisons with the The Cancer Genome Atlas dataset showed that the samples from CE regions resembled the proneural, classical, or mesenchymal subtypes of GBM, whereas the samples from the NE regions predominantly resembled the neural subtype. Computational deconvolution of the RNA-seq data revealed that contributions from nonneoplastic brain cells significantly influence the expression pattern in the NE samples. Gene ontology analysis showed that the cell type-specific expression patterns were functionally distinct and highly enriched in genes associated with the corresponding cell phenotypes. Comparing the RNA-seq data from the GBM samples to that of nonneoplastic brain revealed that the differentially expressed genes are distributed across multiple cell types. Notably, the patterns of cell type-specific alterations varied between the different GBM subtypes: the NE regions of proneural tumors were enriched in oligodendrocyte progenitor genes, whereas the NE regions of mesenchymal GBM were enriched in astrocytic and microglial genes. These subtype-specific patterns provide new insights into molecular and cellular composition of the infiltrative margins of GBM.Glioma cells diffusely infiltrate the brain and intermingle with neural cells in the surrounding brain tissue, resulting in a complex mixture that includes variable proportions of glioma cells, neurons, and various lineages of reactive or recruited glia. At the infiltrative margins of glioblastoma (GBM), the nonneoplastic brain cells can far outnumber the glioma cells and, therefore, will have a significant effect on the molecular features of the tissue. Expression profiling and whole genome sequencing from hundreds of GBM specimens by The Cancer Genome Atlas (TCGA) has revealed a broad spectrum of genetic alterations and discrete expression signatures or subtypes that stratify the majority of patients (1, 2). These studies analyzed tumor samples that were removed during surgery, but were not radiographically localized and, therefore, do not address the question of how the molecular signature may vary across different regions of a tumor. Recent studies have sampled multiple regions within a GBM and shown that more than one molecular subtype can coexist within a single tumor (3). However, the effect of varying cellular composition on GBM subtype, particularly the contribution of nonneoplastic cells, has not been addressed.GBM typically appears as a contrast-enhancing mass, which represents the highly cellular core of the tumor with vascular proliferation and blood–brain barrier breakdown. This contrast-enhancing (CE) region is typically surrounded by a diffuse, nonenhancing (NE) region of abnormal T2/FLAIR signal, which represents edematous brain tissue with varying numbers of infiltrating glioma cells. The primary treatment of GBM is surgical resection, during which the surgeon removes as much of the CE mass as possible. Thus, molecular and genetic profiling of GBM, including the TCGA effort, has predominantly used samples from the CE regions of tumor. However, it is the NE regions of glioma that are left behind after surgery, which neurooncologists must treat and which inevitably give rise to recurrence. Thus, there is immense prognostic and therapeutic significance to understanding the cellular and molecular features of the NE regions of tumor, yet often these areas are not resected and, therefore, have not been directly studied.There are two major obstacles to this goal. The first is the surgical challenge of radiographically localized sampling of the NE tumor margins. The second is the issue of the complex cellular composition that characterizes these regions of diffuse infiltration. In this study, we have addressed both challenges, and associated distinct molecular and cellular features of the NE regions of GBM with the molecular subtype, as defined by the resected CE regions of the tumor.