Search

Pediatric low-grade gliomas (pLGGs) are the most frequent brain tumors in children and comprise a heterogeneous group of tumors with different locations, histologic subtypes, ages at presentation, and clinical behavior. Tumors frequently respond to treatment with chemotherapy or surgical removal, but they can regrow after a period of quiescence, requiring further therapy. Thus, a deeper understanding of the molecular processes involved in these tumors is required to develop therapeutic strategies that are effective against their disease mechanisms. To better understand the cellular behaviors of this heterogenous group of tumors, we have employed single-cell and single-nuclei RNA sequencing technologies to analyze a large-scale dataset (>250,000 cells) of pLGGs. Analysis of this data identified a heterogenous population of cell types and cell states, detecting mature and progenitor-like astrocytes and oligodendrocytes, as well as cells exhibiting senescence or cycling programs. Moreover, we identify a significant immune infiltrate, comprised primarily of microglia. In addition to heterogeneity within pLGG tumors, heterogeneity between LGG subtypes represents another layer that stratifies pLGG biology. We performed a compositional analysis of the cell types present in these tumors and compared transcription signatures and gene expression programs across shared cellular populations of histologically and genetically distinct pLGGs. Finally, we optimized our integration and batch correction analyses by using external 293T cells as spike in controls during our single-cell and single-nuclei data generation steps to determine the most suitable method for batch-effect removal. Our analysis of human pLGGs at the single-cell and single-nuclei resolution provides critical insight into the heterogenous biological activities that constitute these tumors.

INTRODUCTION: Adamantinomatous craniopharyngiomas (ACPs) are rare brain tumors that primarily occur in children and impact long-term morbidity and mortality. The canonical driver mutation for ACP growth occurs in CTNNB1 and leads to constitutive activation of the Wnt/β-catenin signaling pathway. In this study, we outline the genomic, transcriptomic, and structural variant (SV) landscape in a cohort of 41 ACP samples. METHODS: We performed whole-genome sequencing (WGS) and RNA-sequencing of 41 ACP samples. Matched normal samples were also characterized by WGS. Mutect2 was used to detect single nucleotide variants (SNVs) and indels, and copy number data was generated using the GATK pipeline. SvABA was used to perform SV analysis and to identify significantly recurrent breakpoints and juxtapositions. DESeq2 was used to perform differential gene expression analysis based on clinical and molecular annotation data. RESULTS: 29/41 (70%) of the ACP samples harbored missense mutations in exon 3 of CTNNB1, all of which have previously been reported in ACP tumors. SV analysis identified a median of 11.5 events per tumor. Overall, 9.7% of events were interchromosomal. Of the remainder, the majority (78.6%) were deletions. No SVs occurred within CTNNB1. A positive correlation (r = 0.533) was observed between the frequency of SVs and SNVs within samples. Analysis of significantly recurring breakpoints (SRBs) did not identify recurrent breakpoint events. Differential gene expression analysis comparing samples with and without CTNNB1 variants identified 2,143 differentially expressed genes with q-value < 0.05. CONCLUSION: This study identifies activating mutations in exon 3 of CTNNB1 in a large cohort of ACP samples. We also integrate SV and transcriptomic data to comprehensively investigate ACP tumor genomes and identify putative novel tumorigenic mechanisms that advance our understanding of ACP biology.

Pediatric low-grade gliomas (pLGGs) have excellent survival, however, with current standard of care, most patients suffer lifelong severe sequalae. pLGGs are almost exclusively driven by single activating mutations in the MAPK pathway. Clinical trials with small molecule inhibitors in BRAF-altered pLGGs are showing promising results in early clinical trials, and similar efforts are now underway for FGFR1-altered tumors, however the underlying biology and treatment response has not been thoroughly explored in a pre-clinical setting. To explore the genetic landscape of FGFR altered gliomas we assembled a cohort of 87 patients with FGFR1-4 altered gliomas across Dana-Farber Cancer Institute, Boston Children’s Hospital and Brigham and Women’s Hospital. Within this cohort we observed that pLGGs harboring FGFR1 kinase hotspot mutations (FGFR1-N546K or -K656E) frequently harbored a second alteration associated with activation of the MAPK or mTOR pathways, most commonly in the phosphatase PTPN11, NF1 or within the FGFR1 gene itself. Additionally, we observed two previously described structural variants of FGFR1, an FGFR1 internal kinase tandem duplication (FGFR-ITD) and a fusion with TACC1 (FGFR1:TACC1). The relative impact of the different FGFR1 alterations on oncogenicity, therapeutic response and resistance has not been previously explored. To address this, we have established mouse neural stem cell models overexpressing the structural variants and hot spot mutant FGFR1 alone or in combination with a second alteration. Immunoblotting revealed that the addition of a second alteration attenuated phosphorylation of ERK, AKT and S6 and influenced cell proliferation both in normal growth conditions and in absence of growth factor. Treatment with inhibitors of FGFR (Infigratinib) and MEK (Trametinib) revealed variable sensitivity both targeted therapies, suggesting that treatment of FGFR1 driven pLGG might require tailoring to the specific FGFR1 alteration.