Glioblastoma is an aggressive form of brain cancer with well-established patterns of intra-tumoral heterogeneity implicated in treatment resistance and progression. While regional and single cell transcriptomic variations of glioblastoma have been recently resolved, downstream phenotype-level proteomic programs have yet to be assigned across glioblastoma's hallmark histomorphologic niches. Here, we leverage mass spectrometry to spatially align abundance levels of 4,794 proteins to distinct histologic patterns across 20 patients and propose diverse molecular programs operational within these regional tumor compartments. Using machine learning, we overlay concordant transcriptional information, and define two distinct proteogenomic programs, MYC- and KRAS-axis hereon, that cooperate with hypoxia to produce a tri-dimensional model of intra-tumoral heterogeneity. Moreover, we highlight differential drug sensitivities and relative chemoresistance in glioblastoma cell lines with enhanced KRAS programs. Importantly, these pharmacological differences are less pronounced in transcriptional glioblastoma subgroups suggesting that this model may provide insights for targeting heterogeneity and overcoming therapy resistance.

Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 pathogenic variants. Here, we perform whole-genome sequence (WGS) analysis of tumors from 22 patients with TP53 germline pathogenic variants. We observe somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identify near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurs earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, reveal that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising years prior to tumor diagnosis.

Radiation-induced high-grade gliomas (RIGs) are an incurable late complication of cranial radiation therapy. We performed DNA methylation profiling, RNA-seq, and DNA sequencing on 32 RIG tumors and an in vitro drug screen in two RIG cell lines. We report that based on DNA methylation, RIGs cluster primarily with the pediatric receptor tyrosine kinase I high-grade glioma subtype. Common copy-number alterations include Chromosome (Ch.) 1p loss/1q gain, and Ch. 13q and Ch. 14q loss; focal alterations include PDGFRA and CDK4 gain and CDKN2A and BCOR loss. Transcriptomically, RIGs comprise a stem-like subgroup with lesser mutation burden and Ch. 1p loss and a pro-inflammatory subgroup with greater mutation burden and depleted DNA repair gene expression. Chromothripsis in several RIG samples is associated with extrachromosomal circular DNA-mediated amplification of PDGFRA and CDK4. Drug screening suggests microtubule inhibitors/stabilizers, DNA-damaging agents, MEK inhibition, and, in the inflammatory subgroup, proteasome inhibitors, as potentially effective therapies.

Infant gliomas have paradoxical clinical behavior compared to those in children and adults: low-grade tumors have a higher mortality rate, while high-grade tumors have a better outcome. However, we have little understanding of their biology and therefore cannot explain this behavior nor what constitutes optimal clinical management. Here we report a comprehensive genetic analysis of an international cohort of clinically annotated infant gliomas, revealing 3 clinical subgroups. Group 1 tumors arise in the cerebral hemispheres and harbor alterations in the receptor tyrosine kinases ALK, ROS1, NTRK and MET. These are typically single-events and confer an intermediate outcome. Groups 2 and 3 gliomas harbor RAS/MAPK pathway mutations and arise in the hemispheres and midline, respectively. Group 2 tumors have excellent long-term survival, while group 3 tumors progress rapidly and do not respond well to chemoradiation. We conclude that infant gliomas comprise 3 subgroups, justifying the need for specialized therapeutic strategies.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.