Single-Cell RNA-Seq Reveals Cellular Hierarchies and Impaired Developmental Trajectories in Pediatric Ependymoma

TULIPs decorate the three-dimensional genome of PFA ependymoma

Posterior fossa group A (PFA) ependymoma is a lethal brain cancer diagnosed in infants and young children. The lack of driver events in the PFA linear genome led us to search its 3D genome for characteristic features. Here, we reconstructed 3D genomes from diverse childhood tumor types and uncovered a global topology in PFA that is highly reminiscent of stem and progenitor cells in a variety of human tissues. A remarkable feature exclusively present in PFA are type B ultra long-range interactions in PFAs (TULIPs), regions separated by great distances along the linear genome that interact with each other in the 3D nuclear space with surprising strength. TULIPs occur in all PFA samples and recur at predictable genomic coordinates, and their formation is induced by expression of EZHIP. The universality of TULIPs across PFA samples suggests a conservation of molecular principles that could be exploited therapeutically.

GABAergic neuronal lineage development determines clinically actionable targets in diffuse hemispheric glioma, H3G34-mutant

Diffuse hemispheric gliomas, H3G34R/V-mutant (DHG-H3G34), are lethal brain tumors lacking targeted therapies. They originate from interneuronal precursors; however, leveraging this origin for therapeutic insights remains unexplored. Here, we delineate a cellular hierarchy along the interneuron lineage development continuum, revealing that DHG-H3G34 mirror spatial patterns of progenitor streams surrounding interneuron nests, as seen during human brain development. Integrating these findings with genome-wide CRISPR-Cas9 screens identifies genes upregulated in interneuron lineage progenitors as major dependencies. Among these, CDK6 emerges as a targetable vulnerability: DHG-H3G34 tumor cells show enhanced sensitivity to CDK4/6 inhibitors and a CDK6-specific degrader, promoting a shift toward more mature interneuron-like states, reducing tumor growth, and prolonging xenograft survival. Notably, a patient with progressive DHG-H3G34 treated with a CDK4/6 inhibitor achieved 17 months of stable disease. This study underscores interneuronal progenitor-like states, organized in characteristic niches, as a distinct vulnerability in DHG-H3G34, highlighting CDK6 as a promising clinically actionable target.

Multiparametric MRI-Based Interpretable Radiomics Machine Learning Model Differentiates Medulloblastoma and Ependymoma in Children: A Two-Center Study

RATIONALE AND OBJECTIVES

Medulloblastoma (MB) and Ependymoma (EM) in children, share similarities in age group, tumor location, and clinical presentation. Distinguishing between them through clinical diagnosis is challenging. This study aims to explore the effectiveness of using radiomics and machine learning on multiparametric magnetic resonance imaging (MRI) to differentiate between MB and EM and validate its diagnostic ability with an external set.

MATERIALS AND METHODS

Axial T2 weighted image (T2WI) and contrast-enhanced T1weighted image (CE-T1WI) MRI sequences of 135 patients from two centers were collected as train/test sets. Volume of interest (VOI) was manually delineated by an experienced neuroradiologist, supervised by a senior. Feature selection analysis and the least absolute shrinkage and selection operator (LASSO) algorithm identified valuable features, and Shapley additive explanations (SHAP) evaluated their significance. Five machine-learning classifiers-extreme gradient boosting (XGBoost), Bernoulli naive Bayes (Bernoulli NB), Logistic Regression (LR), support vector machine (SVM), linear support vector machine (Linear SVC) classifiers were built based on T2WI (T2 model), CE-T1WI (T1 model), and T1 + T2WI (T1 + T2 model). A human expert diagnosis was developed and corrected by senior radiologists. External validation was performed at Sun Yat-Sen University Cancer Center.

RESULTS

31 valuable features were extracted from T2WI and CE-T1WI. XGBoost demonstrated the highest performance with an area under the curve (AUC) of 0.92 on the test set and maintained an AUC of 0.80 during external validation. For the T1 model, XGBoost achieved the highest AUC of 0.85 on the test set and the highest accuracy of 0.71 on the external validation set. In the T2 model, XGBoost achieved the highest AUC of 0.86 on the test set and the highest accuracy of 0.82 on the external validation set. The human expert diagnosis had an AUC of 0.66 on the test set and 0.69 on the external validation set. The integrated T1 + T2 model achieved an AUC of 0.92 on the test set, 0.80 on the external validation set, achieved the best performance. Overall, XGBoost consistently outperformed in different classification models.

CONCLUSION

The combination of radiomics and machine learning on multiparametric MRI effectively distinguishes between MB and EM in childhood, surpassing human expert diagnosis in training and testing sets.

Single-cell transcriptomic analysis identifies downregulated phosphodiesterase 8B as a novel oncogene in IDH-mutant glioma

Introduction

Glioma, a prevalent and deadly brain tumor, is marked by significant cellular heterogeneity and metabolic alterations. However, the comprehensive cell-of-origin and metabolic landscape in high-grade (Glioblastoma Multiforme, WHO grade IV) and low-grade (Oligoastrocytoma, WHO grade II) gliomas remains elusive.

Methods

In this study, we undertook single-cell transcriptome sequencing of these glioma grades to elucidate their cellular and metabolic distinctions. Following the identification of cell types, we compared metabolic pathway activities and gene expressions between high-grade and low-grade gliomas.

Results

Notably, astrocytes and oligodendrocyte progenitor cells (OPCs) exhibited the most substantial differences in both metabolic pathways and gene expression, indicative of their distinct origins. The comprehensive analysis identified the most altered metabolic pathways (MCPs) and genes across all cell types, which were further validated against TCGA and CGGA datasets for clinical relevance.

Discussion

Crucially, the metabolic enzyme phosphodiesterase 8B (PDE8B) was found to be exclusively expressed and progressively downregulated in astrocytes and OPCs in higher-grade gliomas. This decreased expression identifies PDE8B as a metabolism-related oncogene in IDH-mutant glioma, marking its dual role as both a protective marker for glioma grading and prognosis and as a facilitator in glioma progression.

Developmental origins shape the paediatric cancer genome

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Dissecting tumor transcriptional heterogeneity from single-cell RNA-seq data by generalized binary covariance decomposition

Profiling tumors with single-cell RNA sequencing (scRNA-seq) has the potential to identify recurrent patterns of transcription variation related to cancer progression, and produce new therapeutically relevant insights. However, the presence of strong inter-tumor heterogeneity often obscures more subtle patterns that are shared across tumors, some of which may characterize clinically relevant disease subtypes. Here we introduce a new statistical method, generalized binary covariance decomposition (GBCD), to address this problem. We show that GBCD can help decompose transcriptional heterogeneity into interpretable components - including patient-specific, dataset-specific and shared components relevant to disease subtypes - and that, in the presence of strong inter-tumor heterogeneity, it can produce more interpretable results than existing methods. Applied to data from three studies on pancreatic cancer adenocarcinoma (PDAC), GBCD produces a refined characterization of existing tumor subtypes (e.g., classical vs. basal), and identifies a new gene expression program (GEP) that is prognostic of poor survival independent of established prognostic factors such as tumor stage and subtype. The new GEP is enriched for genes involved in a variety of stress responses, and suggests a potentially important role for the integrated stress response in PDAC development and prognosis.

Identifying tumor type and cell type-specific gene expression alterations in pediatric central nervous system tumors

Central nervous system (CNS) tumors are the leading cause of pediatric cancer death, and these patients have an increased risk for developing secondary neoplasms. Due to the low prevalence of pediatric CNS tumors, major advances in targeted therapies have been lagging compared to other adult tumors. We collect single nuclei RNA-seq data from 84,700 nuclei of 35 pediatric CNS tumors and three non-tumoral pediatric brain tissues and characterize tumor heterogeneity and transcriptomic alterations. We distinguish cell subpopulations associated with specific tumor types including radial glial cells in ependymomas and oligodendrocyte precursor cells in astrocytomas. In tumors, we observe pathways important in neural stem cell-like populations, a cell type previously associated with therapy resistance. Lastly, we identify transcriptomic alterations among pediatric CNS tumor types compared to non-tumor tissues, while accounting for cell type effects on gene expression. Our results suggest potential tumor type and cell type-specific targets for pediatric CNS tumor treatment. Here we address current gaps in understanding single nuclei gene expression profiles of previously under-investigated tumor types and enhance current knowledge of gene expression profiles of single cells of various pediatric CNS tumors.

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Direct Implantation of Patient Brain Tumor Cells into Matching Locations in Mouse Brains for Patient-Derived Orthotopic Xenograft Model Development

Background: Despite multimodality therapies, the prognosis of patients with malignant brain tumors remains extremely poor. One of the major obstacles that hinders development of effective therapies is the limited availability of clinically relevant and biologically accurate (CRBA) mouse models. Methods: We have developed a freehand surgical technique that allows for rapid and safe injection of fresh human brain tumor specimens directly into the matching locations (cerebrum, cerebellum, or brainstem) in the brains of SCID mice. Results: Using this technique, we successfully developed 188 PDOX models from 408 brain tumor patient samples (both high-and low-grade) with a success rate of 72.3% in high-grade glioma, 64.2% in medulloblastoma, 50% in ATRT, 33.8% in ependymoma, and 11.6% in low-grade gliomas. Detailed characterization confirmed their replication of the histopathological and genetic abnormalities of the original patient tumors. Conclusions: The protocol is easy to follow, without a sterotactic frame, in order to generate large cohorts of tumor-bearing mice to meet the needs of biological studies and preclinical drug testing.

Development of Chromosome 1q+ Specific Treatment for Highest Risk Pediatric Posterior Fossa Ependymoma

PURPOSE

There are no effective treatment strategies for children with highest-risk posterior fossa group A ependymoma (PFA). Chromosome 1q gains (1q+) are present in approximately 25% of newly diagnosed PFA tumors, and this number doubles at recurrence. Seventy percent of children with chromosome 1q+ PFA will die because of the tumor, highlighting the urgent need to develop new therapeutic strategies for this population.

EXPERIMENTAL DESIGN

In this study, we utilize 1q+ PFA in vitro and in vivo models to test the efficacy of combination radiation and chemotherapy in a preclinical setting.

RESULTS

5-fluorouracil (5FU) enhances radiotherapy in 1q+ PFA cell lines. Specifically, 5FU increases p53 activity mediated by the extra copy of UCK2 located on chromosome 1q in 1q+ PFA. Experimental downregulation of UCK2 resulted in decreased 5FU sensitivity in 1q+ PFA cells. In in vitro studies, a combination of 5FU, retinoid tretinoin (ATRA), and radiation provided the greatest reduction in cellular proliferation and greatest increase in markers of apoptosis in 1q+ PFA cell lines compared with other treatment arms. Similarly, in vivo experiments demonstrated significant enhancement of survival in mice treated with combination radiation and 5FU and ATRA.

CONCLUSIONS

These results are the first to identify a chromosome 1q+ specific therapy approach in 1q+ PFA. Existing phase I studies have already established single-agent pediatric safety and dosages of 5FU and ATRA, allowing for expedited clinical application as phase II trials for children with high-risk PFA.

Strong OLIG2 expression in supratentorial ependymoma, ZFTA fusion-positive: A potential diagnostic pitfall

Ependymomas (EPN) are central nervous system neoplasms that exhibit an ependymal phenotype. In particular, supratentorial EPN (ST-EPN) must be differentiated from more aggressive entities such as glioblastoma, IDH-wildtype. This task is frequently addressed with the use of immunohistochemistry coupled with clinical presentation and morphological features. Here we describe the case of a young adult presenting with migraine-like symptoms and a temporoinsular-based expansile mass that was first diagnosed as a GBM, mostly based on strong and diffuse oligodendrocyte transcription factor 2 (OLIG2) expression. Molecular characterization revealed a ZFTA::RELA fusion, supporting the diagnosis of ST-EPN, ZFTA fusion-positive. OLIG2 expression is rarely reported in tumors other than GBM and oligodendrocyte-lineage committed neoplasms. The patient was treated with radiotherapy and temozolomide after surgery and was alive and well at follow-up. This report illustrates the need to assess immunostains within a broader clinical, morphological and molecular context to avoid premature exclusion of important differential diagnoses.

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

In vitro and in vivo modeling systems of supratentorial ependymomas

Ependymomas are rare brain tumors that can occur in both children and adults. Subdivided by the tumors' initial location, ependymomas develop in the central nervous system in the supratentorial or infratentorial/posterior fossa region, or the spinal cord. Supratentorial ependymomas (ST-EPNs) are predominantly characterized by common driver gene fusions such as ZFTA and YAP1 fusions. Some variants of ST-EPNs carry a high overall survival rate. In poorly responding ST-EPN variants, high levels of inter- and intratumoral heterogeneity, limited therapeutic strategies, and tumor recurrence are among the reasons for poor patient outcomes with other ST-EPN subtypes. Thus, modeling these molecular profiles is key in further studying tumorigenesis. Due to the scarcity of patient samples, the development of preclinical in vitro and in vivo models that recapitulate patient tumors is imperative when testing therapeutic approaches for this rare cancer. In this review, we will survey ST-EPN modeling systems, addressing the strengths and limitations, application for therapeutic targeting, and current literature findings.

Clinically relevant molecular hallmarks of PFA ependymomas display intratumoral heterogeneity and correlate with tumor morphology

Posterior fossa type A (PF-EPN-A, PFA) ependymoma are aggressive tumors that mainly affect children and have a poor prognosis. Histopathology shows significant intratumoral heterogeneity, ranging from loose tissue to often sharply demarcated, extremely cell-dense tumor areas. To determine molecular differences in morphologically different areas and to understand their clinical significance, we analyzed 113 PF-EPN-A samples, including 40 corresponding relapse samples. Cell-dense areas ranged from 0 to 100% of the tumor area and displayed a higher proportion of proliferating tumor cells (p < 0.01). Clinically, cell density was associated with poor progression-free and overall survival (pPFS = 0.0026, pOS < 0.01). Molecularly, tumor areas with low and high cell density showed diverging DNA methylation profiles regarding their similarity to distinct previously discovered PF-EPN-A subtypes in 9/21 cases. Prognostically relevant chromosomal changes at 1q and 6q showed spatial heterogeneity within single tumors and were significantly enriched in cell-dense tumor areas as shown by single-cell RNA (scRNA)-sequencing as well as copy number profiling and fluorescence in situ hybridization (FISH) analyses of different tumor areas. Finally, spatial transcriptomics revealed cell-dense areas of different tumors to be more similar than various different areas of the same tumor. High-density areas distinctly overexpressed genes encoding histone proteins, WNT5A, TGFB1, or IGF2. Relapsing tumors displayed a higher proportion of cell-dense areas (p = 0.036), a change in PF-EPN-A methylation subtypes (13/32 patients), and novel chromosome 1q gains and 6q losses (12/32 cases) compared to corresponding primary tumors. Our data suggest that PF-EPN-A ependymomas habor a previously unrecognized intratumoral heterogeneity with clinical implications, which has to be accounted for when selecting diagnostic material, inter alia, by histological evaluation of the proportion of cell-dense areas.

Transcriptomic and epigenetic dissection of spinal ependymoma (SP-EPN) identifies clinically relevant subtypes enriched for tumors with and without NF2 mutation

Ependymomas encompass multiple clinically relevant tumor types based on localization and molecular profiles. Tumors of the methylation class "spinal ependymoma" (SP-EPN) represent the most common intramedullary neoplasms in children and adults. However, their developmental origin is ill-defined, molecular data are scarce, and the potential heterogeneity within SP-EPN remains unexplored. The only known recurrent genetic events in SP-EPN are loss of chromosome 22q and NF2 mutations, but neither types and frequency of these alterations nor their clinical relevance have been described in a large, epigenetically defined series. Transcriptomic (n = 72), epigenetic (n = 225), genetic (n = 134), and clinical data (n = 112) were integrated for a detailed molecular overview on SP-EPN. Additionally, we mapped SP-EPN transcriptomes to developmental atlases of the developing and adult spinal cord to uncover potential developmental origins of these tumors. The integration of transcriptomic ependymoma data with single-cell atlases of the spinal cord revealed that SP-EPN display the highest similarities to mature adult ependymal cells. Unsupervised hierarchical clustering of transcriptomic data together with integrated analysis of methylation profiles identified two molecular SP-EPN subtypes. Subtype A tumors primarily carried previously known germline or sporadic NF2 mutations together with 22q loss (bi-allelic NF2 loss), resulting in decreased NF2 expression. Furthermore, they more often presented as multilocular disease and demonstrated a significantly reduced progression-free survival as compared to SP-EP subtype B. In contrast, subtype B predominantly contained samples without NF2 mutation detected in sequencing together with 22q loss (monoallelic NF2 loss). These tumors showed regular NF2 expression but more extensive global copy number alterations. Based on integrated molecular profiling of a large multi-center cohort, we identified two distinct SP-EPN subtypes with important implications for genetic counseling, patient surveillance, and drug development priorities.

Comparative Transcriptomic Analysis of Cerebellar Astrocytes across Developmental Stages and Brain Regions

Astrocytes are the most abundant glial cell type in the central nervous system, and they play a crucial role in normal brain function. While gliogenesis and glial differentiation occur during perinatal cerebellar development, the processes that occur during early postnatal development remain obscure. In this study, we conducted transcriptomic profiling of postnatal cerebellar astrocytes at postnatal days 1, 7, 14, and 28 (P1, P7, P14, and P28), identifying temporal-specific gene signatures at each specific time point. Comparing these profiles with region-specific astrocyte differentially expressed genes (DEGs) published for the cortex, hippocampus, and olfactory bulb revealed cerebellar-specific gene signature across these developmental timepoints. Moreover, we conducted a comparative analysis of cerebellar astrocyte gene signatures with gene lists from pediatric brain tumors of cerebellar origin, including ependymoma and medulloblastoma. Notably, genes downregulated at P14, such as Kif11 and HMGB2, exhibited significant enrichment across all pediatric brain tumor groups, suggesting the importance of astrocytic gene repression during cerebellar development to these tumor subtypes. Collectively, our studies describe gene expression patterns during cerebellar astrocyte development, with potential implications for pediatric tumors originating in the cerebellum.

Patient- and xenograft-derived organoids recapitulate pediatric brain tumor features and patient treatments

Brain tumors are the leading cause of cancer-related death in children. Experimental in vitro models that faithfully capture the hallmarks and tumor heterogeneity of pediatric brain cancers are limited and hard to establish. We present a protocol that enables efficient generation, expansion, and biobanking of pediatric brain cancer organoids. Utilizing our protocol, we have established patient-derived organoids (PDOs) from ependymomas, medulloblastomas, low-grade glial tumors, and patient-derived xenograft organoids (PDXOs) from medulloblastoma xenografts. PDOs and PDXOs recapitulate histological features, DNA methylation profiles, and intratumor heterogeneity of the tumors from which they were derived. We also showed that PDOs can be xenografted. Most interestingly, when subjected to the same routinely applied therapeutic regimens, PDOs respond similarly to the patients. Taken together, our study highlights the potential of PDOs and PDXOs for research and translational applications for personalized medicine.

Inferring cellular and molecular processes in single-cell data with non-negative matrix factorization using Python, R and GenePattern Notebook implementations of CoGAPS

Non-negative matrix factorization (NMF) is an unsupervised learning method well suited to high-throughput biology. However, inferring biological processes from an NMF result still requires additional post hoc statistics and annotation for interpretation of learned features. Here, we introduce a suite of computational tools that implement NMF and provide methods for accurate and clear biological interpretation and analysis. A generalized discussion of NMF covering its benefits, limitations and open questions is followed by four procedures for the Bayesian NMF algorithm Coordinated Gene Activity across Pattern Subsets (CoGAPS). Each procedure will demonstrate NMF analysis to quantify cell state transitions in a public domain single-cell RNA-sequencing dataset. The first demonstrates PyCoGAPS, our new Python implementation that enhances runtime for large datasets, and the second allows its deployment in Docker. The third procedure steps through the same single-cell NMF analysis using our R CoGAPS interface. The fourth introduces a beginner-friendly CoGAPS platform using GenePattern Notebook, aimed at users with a working conceptual knowledge of data analysis but without a basic proficiency in the R or Python programming language. We also constructed a user-facing website to serve as a central repository for information and instructional materials about CoGAPS and its application programming interfaces. The expected timing to setup the packages and conduct a test run is around 15 min, and an additional 30 min to conduct analyses on a precomputed result. The expected runtime on the user's desired dataset can vary from hours to days depending on factors such as dataset size or input parameters.

The molecular landscape of neurological disorders: insights from single-cell RNA sequencing in neurology and neurosurgery

Single-cell ribonucleic acid sequencing (scRNA-seq) has emerged as a transformative technology in neurological and neurosurgical research, revolutionising our comprehension of complex neurological disorders. In brain tumours, scRNA-seq has provided valuable insights into cancer heterogeneity, the tumour microenvironment, treatment resistance, and invasion patterns. It has also elucidated the brain tri-lineage cancer hierarchy and addressed limitations of current models. Neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis have been molecularly subtyped, dysregulated pathways have been identified, and potential therapeutic targets have been revealed using scRNA-seq. In epilepsy, scRNA-seq has explored the cellular and molecular heterogeneity underlying the condition, uncovering unique glial subpopulations and dysregulation of the immune system. ScRNA-seq has characterised distinct cellular constituents and responses to spinal cord injury in spinal cord diseases, as well as provided molecular signatures of various cell types and identified interactions involved in vascular remodelling. Furthermore, scRNA-seq has shed light on the molecular complexities of cerebrovascular diseases, such as stroke, providing insights into specific genes, cell-specific expression patterns, and potential therapeutic interventions. This review highlights the potential of scRNA-seq in guiding precision medicine approaches, identifying clinical biomarkers, and facilitating therapeutic discovery. However, challenges related to data analysis, standardisation, sample acquisition, scalability, and cost-effectiveness need to be addressed. Despite these challenges, scRNA-seq has the potential to transform clinical practice in neurological and neurosurgical research by providing personalised insights and improving patient outcomes.

Significant increase of high-risk chromosome 1q gain and 6q loss at recurrence in posterior fossa group A ependymoma: A multicenter study

BACKGROUND

Ependymoma (EPN) posterior fossa group A (PFA) has the highest rate of recurrence and the worst prognosis of all EPN molecular groups. At relapse, it is typically incurable even with re-resection and re-irradiation. The biology of recurrent PFA remains largely unknown; however, the increasing use of surgery at first recurrence has now provided access to clinical samples to facilitate a better understanding of this.

METHODS

In this large longitudinal international multicenter study, we examined matched samples of primary and recurrent disease from PFA patients to investigate the biology of recurrence.

RESULTS

DNA methylome derived copy number variants (CNVs) revealed large-scale chromosome gains and losses at recurrence in PFA. CNV changes were dominated by chromosome 1q gain and/or 6q loss, both previously identified as high-risk factors in PFA, which were present in 23% at presentation but increased to 61% at first recurrence. Multivariate survival analyses of this cohort showed that cases with 1q gain or 6q loss at first recurrence were significantly more likely to recur again. Predisposition to 1q+/6q- CNV changes at recurrence correlated with hypomethylation of heterochromatin-associated DNA at presentation. Cellular and molecular analyses revealed that 1q+/6q- PFA had significantly higher proportions of proliferative neuroepithelial undifferentiated progenitors and decreased differentiated neoplastic subpopulations.

CONCLUSIONS

This study provides clinically and preclinically actionable insights into the biology of PFA recurrence. The hypomethylation predisposition signature in PFA is a potential risk-classifier for trial stratification. We show that the cellular heterogeneity of PFAs evolves largely because of genetic evolution of neoplastic cells.

Mapping pediatric brain tumors to their origins in the developing cerebellum

BACKGROUND

Distinguishing the cellular origins of childhood brain tumors is key for understanding tumor initiation and identifying lineage-restricted, tumor-specific therapeutic targets. Previous strategies to map the cell-of-origin typically involved comparing human tumors to murine embryonal tissues, which is potentially limited due to species-specific differences. The aim of this study was to unravel the cellular origins of the 3 most common pediatric brain tumors, ependymoma, pilocytic astrocytoma, and medulloblastoma, using a developing human cerebellar atlas.

METHODS

We used a single-nucleus atlas of the normal developing human cerebellum consisting of 176 645 cells as a reference for an in-depth comparison to 4416 bulk and single-cell transcriptome tumor datasets, using gene set variation analysis, correlation, and single-cell matching techniques.

RESULTS

We find that the astroglial cerebellar lineage is potentially the origin for posterior fossa ependymomas. We propose that infratentorial pilocytic astrocytomas originate from the oligodendrocyte lineage and MHC II genes are specifically enriched in these tumors. We confirm that SHH and Group 3/4 medulloblastomas originate from the granule cell and unipolar brush cell lineages. Radiation-induced gliomas stem from cerebellar glial lineages and demonstrate distinct origins from the primary medulloblastoma. We identify tumor genes that are expressed in the cerebellar lineage of origin, and genes that are tumor specific; both gene sets represent promising therapeutic targets for future study.

CONCLUSION

Based on our results, individual cells within a tumor may resemble different cell types along a restricted developmental lineage. Therefore, we suggest that tumors can arise from multiple cellular states along the cerebellar "lineage of origin."

Multi-omic approach identifies hypoxic tumor-associated myeloid cells that drive immunobiology of high-risk pediatric ependymoma

Ependymoma (EPN) is a devastating childhood brain tumor. Single-cell analyses have illustrated the cellular heterogeneity of EPN tumors, identifying multiple neoplastic cell states including a mesenchymal-differentiated subpopulation which characterizes the PFA1 subtype. Here, we characterize the EPN immune environment, in the context of both tumor subtypes and tumor cell subpopulations using single-cell sequencing (scRNAseq, n = 27), deconvolution of bulk tumor gene expression (n = 299), spatial proteomics (n = 54), and single-cell cytokine release assays (n = 12). We identify eight distinct myeloid-derived subpopulations from which a group of cells, termed hypoxia myeloid cells, demonstrate features of myeloid-derived suppressor cells, including IL6/STAT3 pathway activation and wound healing ontologies. In PFA tumors, hypoxia myeloid cells colocalize with mesenchymal-differentiated cells in necrotic and perivascular niches and secrete IL-8, which we hypothesize amplifies the EPN immunosuppressive microenvironment. This myeloid cell-driven immunosuppression will need to be targeted for immunotherapy to be effective in this difficult-to-cure childhood brain tumor.

Transcriptional immunogenomic analysis reveals distinct immunological clusters in paediatric nervous system tumours

BACKGROUND

Cancer immunotherapies including immune checkpoint inhibitors and Chimeric Antigen Receptor (CAR) T-cell therapy have shown variable response rates in paediatric patients highlighting the need to establish robust biomarkers for patient selection. While the tumour microenvironment in adults has been widely studied to delineate determinants of immune response, the immune composition of paediatric solid tumours remains relatively uncharacterized calling for investigations to identify potential immune biomarkers.

METHODS

To inform immunotherapy approaches in paediatric cancers with embryonal origin, we performed an immunogenomic analysis of RNA-seq data from 925 treatment-naïve paediatric nervous system tumours (pedNST) spanning 12 cancer types from three publicly available data sets.

RESULTS

Within pedNST, we uncovered four broad immune clusters: Paediatric Inflamed (10%), Myeloid Predominant (30%), Immune Neutral (43%) and Immune Desert (17%). We validated these clusters using immunohistochemistry, methylation immune inference and segmentation analysis of tissue images. We report shared biology of these immune clusters within and across cancer types, and characterization of specific immune cell frequencies as well as T- and B-cell repertoires. We found no associations between immune infiltration levels and tumour mutational burden, although molecular cancer entities were enriched within specific immune clusters.

CONCLUSIONS

Given the heterogeneity of immune infiltration within pedNST, our findings suggest personalized immunogenomic profiling is needed to guide selection of immunotherapeutic strategies.

Immunotherapy Approaches in Isocitrate-Dehydrogenase-Mutant Low-Grade Glioma

Low-grade gliomas (LGGs) are slow-growing tumors in the central nervous system (CNS). Patients characteristically show the onset of seizures or neurological deficits due to the predominant LGG location in high-functional brain areas. As a molecular hallmark, LGGs display mutations in the isocitrate dehydrogenase (IDH) enzymes, resulting in an altered cellular energy metabolism and the production of the oncometabolite D-2-hydroxyglutarate. Despite the remarkable progress in improving the extent of resection and adjuvant radiotherapy and chemotherapy, LGG remains incurable, and secondary malignant transformation is often observed. Therefore, novel therapeutic approaches are urgently needed. In recent years, immunotherapeutic strategies have led to tremendous success in various cancer types, but the effect of immunotherapy against glioma has been limited due to several challenges, such as tumor heterogeneity and the immunologically "cold" tumor microenvironment. Nevertheless, recent preclinical and clinical findings from immunotherapy trials are encouraging and offer a glimmer of hope for treating IDH-mutant LGG patients. Here, we aim to review the lessons learned from trials involving vaccines, T-cell therapies, and IDH-mutant inhibitors and discuss future approaches to enhance the efficacy of immunotherapies in IDH-mutant LGG.

Cancer Stem Cells in Tumours of the Central Nervous System in Children: A Comprehensive Review

Cancer stem cells (CSCs) are a subgroup of cells found in various kinds of tumours with stem cell characteristics, such as self-renewal, induced differentiation, and tumourigenicity. The existence of CSCs is regarded as a major source of tumour recurrence, metastasis, and resistance to conventional chemotherapy and radiation treatment. Tumours of the central nervous system (CNS) are the most common solid tumours in children, which have many different types including highly malignant embryonal tumours and midline gliomas, and low-grade gliomas with favourable prognoses. Stem cells from the CNS tumours have been largely found and reported by researchers in the last decade and their roles in tumour biology have been deeply studied. However, the cross-talk of CSCs among different CNS tumour types and their clinical impacts have been rarely discussed. This article comprehensively reviews the achievements in research on CSCs in paediatric CNS tumours. Biological functions, diagnostic values, and therapeutic perspectives are reviewed in detail. Further investigations into CSCs are warranted to improve the clinical practice in treating children with CNS tumours.

Ribonucleotide reductase regulatory subunit M2 drives glioblastoma TMZ resistance through modulation of dNTP production

During therapy, adaptations driven by cellular plasticity are partly responsible for driving the inevitable recurrence of glioblastoma (GBM). To investigate plasticity-induced adaptation during standard-of-care chemotherapy temozolomide (TMZ), we performed in vivo single-cell RNA sequencing in patient-derived xenograft (PDX) tumors of GBM before, during, and after therapy. Comparing single-cell transcriptomic patterns identified distinct cellular populations present during TMZ therapy. Of interest was the increased expression of ribonucleotide reductase regulatory subunit M2 (RRM2), which we found to regulate dGTP and dCTP production vital for DNA damage response during TMZ therapy. Furthermore, multidimensional modeling of spatially resolved transcriptomic and metabolomic analysis in patients' tissues revealed strong correlations between RRM2 and dGTP. This supports our data that RRM2 regulates the demand for specific dNTPs during therapy. In addition, treatment with the RRM2 inhibitor 3-AP (Triapine) enhances the efficacy of TMZ therapy in PDX models. We present a previously unidentified understanding of chemoresistance through critical RRM2-mediated nucleotide production.

Applying single cell multi-omic analyses to understand treatment resistance in pediatric high grade glioma

Despite improvements in cancer patient outcomes seen in the past decade, tumor resistance to therapy remains a major impediment to achieving durable clinical responses. Intratumoral heterogeneity related to genetic, epigenetic, transcriptomic, proteomic, and metabolic differences between individual cancer cells has emerged as a driver of therapeutic resistance. This cell to cell heterogeneity can be assessed using single cell profiling technologies that enable the identification of tumor cell clones that exhibit similar defining features like specific mutations or patterns of DNA methylation. Single cell profiling of tumors before and after treatment can generate new insights into the cancer cell characteristics that confer therapeutic resistance by identifying intrinsically resistant sub-populations that survive treatment and by describing new cellular features that emerge post-treatment due to tumor cell evolution. Integrative, single cell analytical approaches have already proven advantageous in studies characterizing treatment-resistant clones in cancers where pre- and post-treatment patient samples are readily available, such as leukemia. In contrast, little is known about other cancer subtypes like pediatric high grade glioma, a class of heterogeneous, malignant brain tumors in children that rapidly develop resistance to multiple therapeutic modalities, including chemotherapy, immunotherapy, and radiation. Leveraging single cell multi-omic technologies to analyze naïve and therapy-resistant glioma may lead to the discovery of novel strategies to overcome treatment resistance in brain tumors with dismal clinical outcomes. In this review, we explore the potential for single cell multi-omic analyses to reveal mechanisms of glioma resistance to therapy and discuss opportunities to apply these approaches to improve long-term therapeutic response in pediatric high grade glioma and other brain tumors with limited treatment options.

Profiling spatiotemporal gene expression of the developing human spinal cord and implications for ependymoma origin

The spatiotemporal regulation of cell fate specification in the human developing spinal cord remains largely unknown. In this study, by performing integrated analysis of single-cell and spatial multi-omics data, we used 16 prenatal human samples to create a comprehensive developmental cell atlas of the spinal cord during post-conceptional weeks 5-12. This revealed how the cell fate commitment of neural progenitor cells and their spatial positioning are spatiotemporally regulated by specific gene sets. We identified unique events in human spinal cord development relative to rodents, including earlier quiescence of active neural stem cells, differential regulation of cell differentiation and distinct spatiotemporal genetic regulation of cell fate choices. In addition, by integrating our atlas with pediatric ependymomas data, we identified specific molecular signatures and lineage-specific genes of cancer stem cells during progression. Thus, we delineate spatiotemporal genetic regulation of human spinal cord development and leverage these data to gain disease insight.

Spatial transcriptomic analysis delineates epithelial and mesenchymal subpopulations and transition stages in childhood ependymoma

BACKGROUND

The diverse cellular constituents of childhood brain tumor ependymoma, recently revealed by single cell RNA-sequencing, may underly therapeutic resistance. Here we use spatial transcriptomics to further advance our understanding of the tumor microenvironment, mapping cellular subpopulations to the tumor architecture of ependymoma posterior fossa subgroup A (PFA), the commonest and most deadly childhood ependymoma variant.

METHODS

Spatial transcriptomics data from intact PFA sections was deconvoluted to resolve the histological arrangement of neoplastic and non-neoplastic cell types. Key findings were validated using immunohistochemistry, in vitro functional assays and outcome analysis in clinically-annotated PFA bulk transcriptomic data.

RESULTS

PFA are comprised of epithelial and mesenchymal histological zones containing a diversity of cellular states, each zone including co-existing and spatially distinct undifferentiated progenitor-like cells; a quiescent mesenchymal zone population, and a second highly mitotic progenitor population that is restricted to hypercellular epithelial zones and that is more abundant in progressive tumors. We show that myeloid cell interaction is the leading cause of mesenchymal transition in PFA, occurring in zones spatially distinct from hypoxia-induced mesenchymal transition, and these distinct EMT-initiating processes were replicated using in vitro models of PFA.

CONCLUSIONS

These insights demonstrate the utility of spatial transcriptomics to advance our understanding of ependymoma biology, revealing a clearer picture of the cellular constituents of PFA, their interactions and influence on tumor progression.

The Blood-Brain Barrier: Implications for Experimental Cancer Therapeutics

The blood-brain barrier is critically important for the treatment of both primary and metastatic cancers of the central nervous system (CNS). Clinical outcomes for patients with primary CNS tumors are poor and have not significantly improved in decades. As treatments for patients with extracranial solid tumors improve, the incidence of CNS metastases is on the rise due to suboptimal CNS exposure of otherwise systemically active agents. Despite state-of-the art surgical care and increasingly precise radiation therapy, clinical progress is limited by the ability to deliver an effective dose of a therapeutic agent to all cancerous cells. Given the tremendous heterogeneity of CNS cancers, both across cancer subtypes and within a single tumor, and the range of diverse therapies under investigation, a nuanced examination of CNS drug exposure is needed. With a shared goal, common vocabulary, and interdisciplinary collaboration, the field is poised for renewed progress in the treatment of CNS cancers.

Why haven't we solved intracranial pediatric ependymoma? Current questions and barriers to treatment advances

Pediatric intracranial ependymoma has seen a recent exponential expansion of biological findings, rapidly dividing the diagnosis into several subgroups, each with specific molecular and clinical characteristics. While such subdivision may complicate clinical conclusions from historical trials, this knowledge also provides an opportunity for interrogating the major clinical and biological questions preventing near-term translation into effective therapy for children with ependymoma. In this article, we briefly review some of the most critical clinical questions facing both patient management and the construct of future trials in childhood ependymoma, as well as explore some of the current barriers to efficient translation of preclinical discovery to the clinic.

Hydroxymethylation alterations in progenitor-like cell types of pediatric central nervous system tumors are associated with cell type-specific transcriptional changes

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors we utilized a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identified a preponderance differential CpG hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like HDAC4 and IGF1R, were associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric CNS tumors.

Tumor type and cell type-specific gene expression alterations in diverse pediatric central nervous system tumors identified using single nuclei RNA-seq

Central nervous system (CNS) tumors are the leading cause of pediatric cancer death, and these patients have an increased risk for developing secondary neoplasms. Due to the low prevalence of pediatric CNS tumors, major advances in targeted therapies have been lagging compared to other adult tumors. We collected single nuclei RNA-seq data from 35 pediatric CNS tumors and three non-tumoral pediatric brain tissues (84,700 nuclei) and characterized tumor heterogeneity and transcriptomic alterations. We distinguished cell subpopulations associated with specific tumor types including radial glial cells in ependymomas and oligodendrocyte precursor cells in astrocytomas. In tumors, we observed pathways important in neural stem cell-like populations, a cell type previously associated with therapy resistance. Lastly, we identified transcriptomic alterations among pediatric CNS tumor types compared to non-tumor tissues, while accounting for cell type effects on gene expression. Our results suggest potential tumor type and cell type-specific targets for pediatric CNS tumor treatment. In this study, we address current gaps in understanding single nuclei gene expression profiles of previously uninvestigated tumor types and enhance current knowledge of gene expression profiles of single cells of various pediatric CNS tumors.

Single-cell profiling of alveolar rhabdomyosarcoma reveals RAS pathway inhibitors as cell-fate hijackers with therapeutic relevance

Rhabdomyosarcoma (RMS) is a group of pediatric cancers with features of developing skeletal muscle. The cellular hierarchy and mechanisms leading to developmental arrest remain elusive. Here, we combined single-cell RNA sequencing, mass cytometry, and high-content imaging to resolve intratumoral heterogeneity of patient-derived primary RMS cultures. We show that the aggressive alveolar RMS (aRMS) subtype contains plastic muscle stem-like cells and cycling progenitors that drive tumor growth, and a subpopulation of differentiated cells that lost its proliferative potential and correlates with better outcomes. While chemotherapy eliminates cycling progenitors, it enriches aRMS for muscle stem-like cells. We screened for drugs hijacking aRMS toward clinically favorable subpopulations and identified a combination of RAF and MEK inhibitors that potently induces myogenic differentiation and inhibits tumor growth. Overall, our work provides insights into the developmental states underlying aRMS aggressiveness, chemoresistance, and progression and identifies the RAS pathway as a promising therapeutic target.

Clustering-independent estimation of cell abundances in bulk tissues using single-cell RNA-seq data

Single-cell RNA-sequencing has transformed the study of biological tissues by enabling transcriptomic characterizations of their constituent cell states. Computational methods for gene expression deconvolution use this information to infer the cell composition of related tissues profiled at the bulk level. However, current deconvolution methods are restricted to discrete cell types and have limited power to make inferences about continuous cellular processes like cell differentiation or immune cell activation. We present ConDecon, a clustering-independent method for inferring the likelihood for each cell in a single-cell dataset to be present in a bulk tissue. ConDecon represents an improvement in functionality and accuracy with respect to current deconvolution methods. Using ConDecon, we discover the implication of neurodegenerative microglial inflammatory pathways in the mesenchymal transformation of ependymoma, recapitulate spatial patterns of cell differentiation during zebrafish embryogenesis, and make temporal inferences from bulk ATAC-seq data. Overall, ConDecon significantly enhances our understanding of dynamic cellular processes within bulk tissue samples.

An ependymal cell census identifies heterogeneous and ongoing cell maturation in the adult mouse spinal cord that changes dynamically on injury

The adult spinal cord stem cell potential resides within the ependymal cell population and declines with age. Ependymal cells are, however, heterogeneous, and the biological diversity this represents and how it changes with age remain unknown. Here, we present a single-cell transcriptomic census of spinal cord ependymal cells from adult and aged mice, identifying not only all known ependymal cell subtypes but also immature as well as mature cell states. By comparing transcriptomes of spinal cord and brain ependymal cells, which lack stem cell abilities, we identify immature cells as potential spinal cord stem cells. Following spinal cord injury, these cells re-enter the cell cycle, which is accompanied by a short-lived reversal of ependymal cell maturation. We further analyze ependymal cells in the human spinal cord and identify widespread cell maturation and altered cell identities. This in-depth characterization of spinal cord ependymal cells provides insight into their biology and informs strategies for spinal cord repair.

Recurrent pediatric infratentorial ependymomas: a systematic review and meta-analysis on outcomes and molecular classification

OBJECTIVE

The aim of this study was to summarize the prognosis of recurrent infratentorial ependymomas based on treatment and molecular characterization.

METHODS

Following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, the authors searched the PubMed, Scopus, Embase, and Ovid databases for studies on recurrent infratentorial ependymomas in patients younger than 25 years of age. Exclusion criteria included case series of fewer than 5 patients and studies that did not provide time-dependent survival data.

RESULTS

The authors' database search yielded 482 unique articles, of which 18 were included in the final analysis. There were 479 recurrent infratentorial pediatric ependymomas reported; 53.4% were WHO grade II and 46.6% were WHO grade III tumors. The overall mortality for recurrent infratentorial pediatric ependymomas was 49.1% (226/460). The pooled mean survival was 30.2 months after recurrence (95% CI 22.4-38.0 months). Gross-total resection (GTR) was achieved in 243 (59.0%) patients at initial presentation. The mean survival postrecurrence for those who received initial GTR was 42.3 months (95% CI 35.7-47.6 months) versus 26.0 months (95% CI 9.6-44.6 months) for those who received subtotal resection (STR) (p = 0.032). There was no difference in the mean survival between patients who received GTR (49.3 months, 95% CI 32.3-66.3 months) versus those who received STR (41.4 months, 95% CI 11.6-71.2 months) for their recurrent tumor (p = 0.610). In the studies that included molecular classification data, there were 169 (83.3%) posterior fossa group A (PFA) tumors and 34 (16.7%) posterior fossa group B (PFB) tumors, with 28 tumors harboring a 1q gain. PFA tumors demonstrated worse mean postprogression patient survival (24.7 months, 95% CI 15.3-34.0 months) compared with PFB tumors (48.0 months, 95% CI 32.8-63.2 months) (p = 0.0073). The average postrecurrence survival for patients with 1q+ tumors was 14.7 months.

CONCLUSIONS

The overall mortality rate for recurrent infratentorial ependymomas was found to be 49.1%, with a pooled mean survival of 30.2 months in the included sample population. More than 80% of recurrent infratentorial ependymomas were of the PFA molecular subtype, and both PFA tumors and those with 1q gain demonstrated worse prognosis after recurrence.

Nuclear condensates of YAP fusion proteins alter transcription to drive ependymoma tumourigenesis

Nuclear localization of HIPPO-YAP fusion proteins has been implicated in supratentorial ependymoma development. Here, unexpectedly, we find that liquid-liquid phase separation, rather than nuclear localization, of recurrent patient-derived YAP fusions, YAP-MAMLD1 and C11ORF95-YAP, underlies ependymoma tumourigenesis from neural progenitor cells. Mutagenesis and chimaera assays demonstrate that an intrinsically disordered region promotes oligomerization of the YAP fusions into nuclear, puncta-like, membrane-less condensates. Oligomerization and nuclear condensates induced by YAP fusion with a coiled-coil domain of transcriptional activator GCN4 also promote ependymoma formation. YAP-MAMLD1 concentrates transcription factors and co-activators, including BRD4, MED1 and TEAD, in condensates while excluding transcriptional repressive PRC2, and induces long-range enhancer-promoter interactions that promote transcription and oncogenic programmes. Blocking condensate-mediated transcriptional co-activator activity inhibits tumourigenesis, indicating a critical role of liquid phase separation for YAP fusion oncogenic activity in ependymoma. YAP fusions containing the intrinsically disordered region features are common in human tumours, suggesting that nuclear condensates could be targeted to treat YAP-fusion-induced cancers.

Pioneering models of pediatric brain tumors

Among children and adolescents in the United States (0 to 19 years old), brain and other central nervous system tumors are the second most common types of cancers, surpassed in incidence only by leukemias. Despite significant progress in the diagnosis and treatment modalities, brain cancer remains the leading cause of death in the pediatric population. There is an obvious unfulfilled need to streamline the therapeutic strategies and improve survival for these patients. For that purpose, preclinical models play a pivotal role. Numerous models are currently used in pediatric brain tumor research, including genetically engineered mouse models, patient-derived xenografts and cell lines, and newer models that utilize novel technologies such as genome engineering and organoids. Furthermore, extensive studies by the Children's Brain Tumor Network (CBTN) researchers and others have revealed multiomic landscapes of variable pediatric brain tumors. Combined with such integrative data, these novel technologies have enabled numerous applicable models. Genome engineering, including CRISPR/Cas9, expanded the flexibility of modeling. Models generated through genome engineering enabled studying particular genetic alterations in clean isogenic backgrounds, facilitating the dissection of functional mechanisms of those mutations in tumor biology. Organoids have been applied to study tumor-to-tumor-microenvironment interactions and to address developmental aspects of tumorigenesis, which is essential in some pediatric brain tumors. Other modalities, such as humanized mouse models, could potentially be applied to pediatric brain tumors. In addition to current valuable models, such novel models are anticipated to expedite functional tumor biology study and establish effective therapeutics for pediatric brain tumors.

Neuronal-Activity Dependent Mechanisms of Small Cell Lung Cancer Progression

Neural activity is increasingly recognized as a critical regulator of cancer growth. In the brain, neuronal activity robustly influences glioma growth both through paracrine mechanisms and through electrochemical integration of malignant cells into neural circuitry via neuron-to-glioma synapses, while perisynaptic neurotransmitter signaling drives breast cancer brain metastasis growth. Outside of the CNS, innervation of tumors such as prostate, breast, pancreatic and gastrointestinal cancers by peripheral nerves similarly regulates cancer progression. However, the extent to which the nervous system regulates lung cancer progression, either in the lung or when metastatic to brain, is largely unexplored. Small cell lung cancer (SCLC) is a lethal high-grade neuroendocrine tumor that exhibits a strong propensity to metastasize to the brain. Here we demonstrate that, similar to glioma, metastatic SCLC cells in the brain co-opt neuronal activity-regulated mechanisms to stimulate growth and progression. Optogenetic stimulation of cortical neuronal activity drives proliferation and invasion of SCLC brain metastases. In the brain, SCLC cells exhibit electrical currents and consequent calcium transients in response to neuronal activity, and direct SCLC cell membrane depolarization is sufficient to promote the growth of SCLC tumors. In the lung, vagus nerve transection markedly inhibits primary lung tumor formation, progression and metastasis, highlighting a critical role for innervation in overall SCLC initiation and progression. Taken together, these studies illustrate that neuronal activity plays a crucial role in dictating SCLC pathogenesis in both primary and metastatic sites.

Activation of Hedgehog signaling by the oncogenic RELA fusion reveals a primary cilia-dependent vulnerability in supratentorial ependymoma

BACKGROUND

Supratentorial RELA fusion (ST-RELA) ependymomas (EPNs) are resistant tumors without an approved chemotherapeutic treatment. Unfortunately, the molecular mechanisms that lead to chemoresistance traits of ST-RELA remain elusive. The aim of this study was to assess RELA fusion-dependent signaling modules, specifically the role of the Hedgehog (Hh) pathway as a novel targetable vulnerability in ST-RELA.

METHODS

Gene expression was analyzed in EPN from patient cohorts, by microarray, RNA-seq, qRT-PCR, and scRNA-seq. Inhibitors against Smoothened (SMO) (Sonidegib) and Aurora kinase A (AURKA) (Alisertib) were evaluated. Protein expression, primary cilia formation, and drug effects were assessed by immunoblot, immunofluorescence, and immunohistochemistry.

RESULTS

Hh components were selectively overexpressed in EPNs induced by the RELA fusion. Single-cell analysis showed that the Hh signature was primarily confined to undifferentiated, stem-like cell subpopulations. Sonidegib exhibited potent growth-inhibitory effects on ST-RELA cells, suggesting a key role in active Hh signaling; importantly, the effect of Sonidegib was reversed by primary cilia loss. We, thus, tested the effect of AURKA inhibition by Alisertib, to induce cilia stabilization/reassembly. Strikingly, Alisertib rescued ciliogenesis and synergized with Sonidegib in killing ST-RELA cells. Using a xenograft model, we show that cilia loss is a mechanism for acquiring resistance to the inhibitory effect of Sonidegib. However, Alisertib fails to rescue cilia and highlights the need for other strategies to promote cilia reassembly, for treating ST-RELA tumors.

CONCLUSION

Our study reveals a crucial role for the Hh pathway in ST-RELA tumor growth, and suggests that rescue of primary cilia represents a vulnerability of the ST-RELA EPNs.

Deciphering functional tumor states at single-cell resolution

Within a tumor, cancer cells exist in different states that are associated with distinct tumor functions, including proliferation, differentiation, invasion, metastasis, and resistance to anti-cancer therapy. The identification of the gene regulatory networks underpinning each state is essential for better understanding functional tumor heterogeneity and revealing tumor vulnerabilities. Here, we review the different studies identifying tumor states by single-cell sequencing approaches and the mechanisms that promote and sustain these functional states and regulate their transitions. We also describe how different tumor states are spatially distributed and interact with the specific stromal cells that compose the tumor microenvironment. Finally, we discuss how the understanding of tumor plasticity and transition states can be used to develop new strategies to improve cancer therapy.

BAF Complex Maintains Glioma Stem Cells in Pediatric H3K27M Glioma

Diffuse midline gliomas are uniformly fatal pediatric central nervous system cancers that are refractory to standard-of-care therapeutic modalities. The primary genetic drivers are a set of recurrent amino acid substitutions in genes encoding histone H3 (H3K27M), which are currently undruggable. These H3K27M oncohistones perturb normal chromatin architecture, resulting in an aberrant epigenetic landscape. To interrogate for epigenetic dependencies, we performed a CRISPR screen and show that patient-derived H3K27M-glioma neurospheres are dependent on core components of the mammalian BAF (SWI/SNF) chromatin remodeling complex. The BAF complex maintains glioma stem cells in a cycling, oligodendrocyte precursor cell-like state, in which genetic perturbation of the BAF catalytic subunit SMARCA4 (BRG1), as well as pharmacologic suppression, opposes proliferation, promotes progression of differentiation along the astrocytic lineage, and improves overall survival of patient-derived xenograft models. In summary, we demonstrate that therapeutic inhibition of the BAF complex has translational potential for children with H3K27M gliomas.

SIGNIFICANCE

Epigenetic dysregulation is at the core of H3K27M-glioma tumorigenesis. Here, we identify the BRG1-BAF complex as a critical regulator of enhancer and transcription factor landscapes, which maintain H3K27M glioma in their progenitor state, precluding glial differentiation, and establish pharmacologic targeting of the BAF complex as a novel treatment strategy for pediatric H3K27M glioma. See related commentary by Beytagh and Weiss, p. 2730. See related article by Mo et al., p. 2906.

The landscape of tumor cell states and spatial organization in H3-K27M mutant diffuse midline glioma across age and location

Histone 3 lysine27-to-methionine (H3-K27M) mutations most frequently occur in diffuse midline gliomas (DMGs) of the childhood pons but are also increasingly recognized in adults. Their potential heterogeneity at different ages and midline locations is vastly understudied. Here, through dissecting the single-cell transcriptomic, epigenomic and spatial architectures of a comprehensive cohort of patient H3-K27M DMGs, we delineate how age and anatomical location shape glioma cell-intrinsic and -extrinsic features in light of the shared driver mutation. We show that stem-like oligodendroglial precursor-like cells, present across all clinico-anatomical groups, display varying levels of maturation dependent on location. We reveal a previously underappreciated relationship between mesenchymal cancer cell states and age, linked to age-dependent differences in the immune microenvironment. Further, we resolve the spatial organization of H3-K27M DMG cell populations and identify a mitotic oligodendroglial-lineage niche. Collectively, our study provides a powerful framework for rational modeling and therapeutic interventions.

Characterizing the biology of primary brain tumors and their microenvironment via single-cell profiling methods

Genomic and transcriptional heterogeneity is prevalent among the most common and aggressive primary brain tumors in children and adults. Over the past 20 years, advances in bioengineering, biochemistry and bioinformatics have enabled the development of an array of techniques to study tumor biology at single-cell resolution. The application of these techniques to study primary brain tumors has helped advance our understanding of their intra-tumoral heterogeneity and uncover new insights regarding their co-option of developmental programs and signaling from their microenvironment to promote tumor proliferation and invasion. These insights are currently being harnessed to develop new therapeutic approaches. Here we provide an overview of current single-cell techniques and discuss relevant biology and therapeutic insights uncovered by their application to primary brain tumors in children and adults.

Single-cell transcriptome reveals cellular hierarchies and guides p-EMT-targeted trial in skull base chordoma

Skull base chordoma (SBC) is a bone cancer with a high recurrence rate, high radioresistance rate, and poorly understood mechanism. Here, we profiled the transcriptomes of 90,691 single cells, revealed the SBC cellular hierarchies, and explored novel treatment targets. We identified a cluster of stem-like SBC cells that tended to be distributed in the inferior part of the tumor. Combining radiated UM-Chor1 RNA-seq data and in vitro validation, we further found that this stem-like cell cluster is marked by cathepsin L (CTSL), a gene involved in the packaging of telomere ends, and may be responsible for radioresistance. Moreover, signatures related to partial epithelial-mesenchymal transition (p-EMT) were found to be significant in malignant cells and were related to the invasion and poor prognosis of SBC. Furthermore, YL-13027, a p-EMT inhibitor that acts through the TGF-β signaling pathway, demonstrated remarkable potency in inhibiting the invasiveness of SBC in preclinical models and was subsequently applied in a phase I clinical trial that enrolled three SBC patients. Encouragingly, YL-13027 attenuated the growth of SBC and achieved stable disease with no serious adverse events, underscoring the clinical potential for the precision treatment of SBC with this therapy. In summary, we conducted the first single-cell RNA sequencing of SBC and identified several targets that could be translated to the treatment of SBC.

TMExplorer: A tumour microenvironment single-cell RNAseq database and search tool

MOTIVATION

The tumour microenvironment (TME) contains various cells including stromal fibroblasts, immune and malignant cells, and its composition can be elucidated using single-cell RNA sequencing (scRNA-seq). scRNA-seq datasets from several cancer types are available, yet we lack a comprehensive database to collect and present related TME data in an easily accessible format.

RESULTS

We therefore built a TME scRNA-seq database, and created the R package TMExplorer to facilitate investigation of the TME. TMExplorer provides an interface to easily access all available datasets and their metadata. The users can search for datasets using a thorough range of characteristics. The TMExplorer allows for examination of the TME using scRNA-seq in a way that is streamlined and allows for easy integration into already existing scRNA-seq analysis pipelines.

Brain cancer stem cells: resilience through adaptive plasticity and hierarchical heterogeneity

Malignant brain tumours are complex ecosystems containing neoplastic and stromal components that generate adaptive and evolutionarily driven aberrant tissues in the central nervous system. Brain cancers are cultivated by a dynamic population of stem-like cells that enforce intratumoural heterogeneity and respond to intrinsic microenvironment or therapeutically guided insults through proliferation, plasticity and restructuring of neoplastic and stromal components. Far from a rigid hierarchy, heterogeneous neoplastic populations transition between cellular states with differential self-renewal capacities, endowing them with powerful resilience. Here we review the biological machinery used by brain tumour stem cells to commandeer tissues in the intracranial space, evade immune responses and resist chemoradiotherapy. Through recent advances in single-cell sequencing, improved models to investigate the role of the tumour microenvironment and a deeper understanding of the fundamental role of the immune system in cancer biology, we are now better equipped to explore mechanisms by which these processes can be exploited for therapeutic benefit.

Pro-inflammatory cytokines mediate the epithelial-to-mesenchymal-like transition of pediatric posterior fossa ependymoma

Pediatric ependymoma is a devastating brain cancer marked by its relapsing pattern and lack of effective chemotherapies. This shortage of treatments is due to limited knowledge about ependymoma tumorigenic mechanisms. By means of single-nucleus chromatin accessibility and gene expression profiling of posterior fossa primary tumors and distal metastases, we reveal key transcription factors and enhancers associated with the differentiation of ependymoma tumor cells into tumor-derived cell lineages and their transition into a mesenchymal-like state. We identify NFκB, AP-1, and MYC as mediators of this transition, and show that the gene expression profiles of tumor cells and infiltrating microglia are consistent with abundant pro-inflammatory signaling between these populations. In line with these results, both TGF-β1 and TNF-α induce the expression of mesenchymal genes on a patient-derived cell model, and TGF-β1 leads to an invasive phenotype. Altogether, these data suggest that tumor gliosis induced by inflammatory cytokines and oxidative stress underlies the mesenchymal phenotype of posterior fossa ependymoma.