Glioblastoma is one of the most common and lethal primary neoplasms of the brain. Patient survival has not improved significantly over the past three decades and the patient median survival is just over one year. Tumor heterogeneity is thought to be a major determinant of therapeutic failure and a major reason for poor overall survival. This work aims to comprehensively define intra- and inter-tumor heterogeneity by mapping the genomic and mutational landscape of multiple areas of three primary IDH wild-type (IDH-WT) glioblastomas. Using whole exome sequencing, we explored how copy number variation, chromosomal and single loci amplifications/deletions, and mutational burden are spatially distributed across nine different tumor regions. The results show that all tumors exhibit a different signature despite the same diagnosis. Above all, a high inter-tumor heterogeneity emerges. The evolutionary dynamics of all identified mutations within each region underline the questionable value of a single biopsy and thus the therapeutic approach for the patient. Multiregional collection and subsequent sequencing are essential to try to address the clinical challenge of precision medicine. Especially in glioblastoma, this approach could provide powerful support to pathologists and oncologists in evaluating the diagnosis and defining the best treatment option.
Keywords: clonal evolution; glioblastoma; multiregional sequencing; spatial heterogeneity; temporal heterogeneity; tumor phylogeny; tumor progression.