Single-Cell RNA Sequencing of Childhood Ependymoma Reveals Neoplastic Cell Subpopulations That Impact Molecular Classification and Etiology

Tumor Dormancy and Reactivation: The Role of Heat Shock Proteins

Tumors are a heterogeneous group of cell masses originating in various organs or tissues. The cellular composition of the tumor cell mass interacts in an intricate manner, influenced by humoral, genetic, molecular, and tumor microenvironment cues that dictate tumor growth or suppression. As a result, tumors undergo a period of a dormant state before their clinically discernible stage, which surpasses the clinical dormancy threshold. Moreover, as a genetically imprinted strategy, early-seeder cells, a distinct population of tumor cells, break off to dock nearby or extravasate into blood vessels to secondary tissues, where they form disseminated solitary dormant tumor cells with reversible capacity. Among the various mechanisms underlying the dormant tumor mass and dormant tumor cell formation, heat shock proteins (HSPs) might play one of the most important roles in how the dormancy program plays out. It is known that numerous aberrant cellular processes, such as malignant transformation, cancer cell stemness, tumor invasion, metastasis, angiogenesis, and signaling pathway maintenance, are influenced by the HSPs. An accumulating body of knowledge suggests that HSPs may be involved in the angiogenic switch, immune editing, and extracellular matrix (ECM) remodeling cascades, crucial genetically imprinted strategies important to the tumor dormancy initiation and dormancy maintenance program. In this review, we highlight the biological events that orchestrate the dormancy state and the body of work that has been conducted on the dynamics of HSPs in a tumor mass, as well as tumor cell dormancy and reactivation. Additionally, we propose a conceptual framework that could possibly underlie dormant tumor reactivation in metastatic relapse.

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.

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.

Multi-pronged analysis of pediatric low-grade glioma reveals a unique tumor microenvironment associated with BRAF alterations

Pediatric low-grade gliomas (pLGG) comprise 35% of all brain tumors. Despite favorable survival, patients experience significant morbidity from disease and treatments. A deeper understanding of pLGG biology is essential to identify novel, more effective, and less toxic therapies. We utilized single cell RNA sequencing (scRNA-seq), spatial transcriptomics, and cytokine analyses to characterize and understand tumor and immune cell heterogeneity across pLGG. scRNA-seq revealed tumor and immune cells within the tumor microenvironment (TME). Tumor cell subsets revealed a developmental hierarchy with progenitor and mature cell populations. Immune cells included myeloid and lymphocytic cells. There was a significant difference between the prevalence of two major myeloid subclusters between pilocytic astrocytoma (PA) and ganglioglioma (GG). Bulk and single-cell cytokine analyses evaluated the immune cell signaling cascade with distinct immune phenotypes among tumor samples. KIAA1549-BRAF tumors appeared more immunogenic, secreting higher levels of immune cell activators and chemokines, compared to BRAF V600E tumors. Spatial transcriptomics revealed the differential gene expression of these chemokines and their location within the TME. A multi-pronged analysis of pLGG demonstrated the complexity of the pLGG TME and differences between genetic drivers that may influence their response to immunotherapy. Further investigation of immune cell infiltration and tumor-immune interactions is warranted.

Key points

There is a developmental hierarchy in neoplastic population comprising of both progenitor-like and mature cell types in both PA and GG.A more immunogenic, immune activating myeloid population is present in PA compared to GG. Functional analysis and spatial transcriptomics show higher levels of immune mobilizing chemokines in KIAA1549-BRAF fusion PA tumor samples compared to BRAF V600E GG samples.

Importance of the Study

While scRNA seq provides information on cellular heterogeneity within the tumor microenvironment (TME), it does not provide a complete picture of how these cells are interacting or where they are located. To expand on this, we used a three-pronged approach to better understand the biology of pediatric low-grade glioma (pLGG). By analyzing scRNA-seq, secreted cytokines and spatial orientation of cells within the TME, we strove to gain a more complete picture of the complex interplay between tumor and immune cells within pLGG. Our data revealed a complex heterogeneity in tumor and immune populations and identified an interesting difference in the immune phenotype among different subtypes.

APOBEC3C is a novel target for the immune treatment of lower-grade gliomas

BACKGROUND

Apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) type 3C (A3C) has been identified as a cancer molecular biomarker in the past decade. However, the practical role of A3C in lower-grade gliomas (LGGs) in improving the clinical outcome remains unclear. This study aims to discuss the function of A3C in immunotherapy in LGGs.

METHODS

The RNA-Sequencing (RNA-seq) and corresponding clinical data were extracted from UCSC Xena and the results were verified in the Chinese Glioma Genome Atlas (CGGA). Weighted gene co-expression network analysis (WGCNA) was used for screening A3C-related genes. Comprehensive bioinformation analyses were performed and multiple levels of expression, survival rate, and biological functions were assessed to explore the functions of A3C.

RESULTS

A3C expression was significantly higher in LGGs than in normal tissues but lower than in glioblastoma (GBM), indicating its role as an independent prognosis predictor for LGGs. Twenty-eight A3C-related genes were found with WGCNA for unsupervised clustering analysis and three modification patterns with different outcomes and immune cell infiltration were identified. A3C and the A3C score were also correlated with immune cell infiltration and the expression of immune checkpoints. In addition, the A3C score was correlated with increased sensitivity to chemotherapy. Single-cell RNA (scRNA) analysis indicated that A3C most probably expresses on immune cells, such as T cells, B cells and macrophage.

CONCLUSIONS

A3C is an immune-related prognostic biomarker in LGGs. Developing drugs to block A3C could enhance the efficiency of immunotherapy and improve disease survival.Abbreviation: A3C: Apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) type 3C; LGGs: lower-grade gliomas; CGGA: Chinese Glioma Genome Atlas; WGCNA: Weighted gene co-expression network analysis; scRNA: Single-cell RNA; HGG: higher-grade glioma; OS: overall survival; TME: tumor microenvironment; KM: Kaplan-Meier; PFI: progression-free interval; IDH: isocitrate dehydrogenase; ROC: receiver operating characteristic; GS: gene significance; MM: module membership; TIMER: Tumor IMmune Estimation Resource; GSVA: gene set variation analysis; ssGSEA: single-sample gene-set enrichment analysis; PCA: principal component analysis; AUC: area under ROC curve; HAVCR2: hepatitis A virus cellular receptor 2; PDCD1: programmed cell death 1; PDCD1LG2: PDCD1 ligand 2; PTPRC: protein tyrosine phosphatase receptor type C; ACC: Adrenocortical carcinoma; BLCA: Bladder Urothelial Carcinoma;BRCA: Breast invasive carcinoma; CESC: Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOLCholangiocarcinoma; COADColon adenocarcinoma; DLBC: Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA: Esophageal carcinoma; GBM: Glioblastoma multiforme; HNSC: Head and Neck squamous cell carcinoma; KICH: Kidney Chromophobe; KIRC: Kidney renal clear cell carcinoma; KIRP: Kidney renal papillary cell carcinoma; LAML: Acute Myeloid Leukemia; LGG: Brain Lower Grade Glioma; LIHC: Liver hepatocellular carcinoma; LUAD: Lung adenocarcinoma; LUSC: Lung squamous cell carcinoma; MESO: Mesothelioma; OV: Ovarian serous cystadenocarcinoma; PAAD: Pancreatic adenocarcinoma; PCPG: Pheochromocytoma and Paraganglioma; PRAD: Prostate adenocarcinoma; READ: Rectum adenocarcinoma; SARC: Sarcoma; SKCM: Skin Cutaneous Melanoma; STAD: Stomach adenocarcinoma; TGCT: Testicular Germ Cell Tumors; THCA: Thyroid carcinoma; THYM: Thymoma; UCEC: Uterine Corpus Endometrial Carcinoma; UCS: Uterine Carcinosarcoma; UVM: Uveal Melanoma.

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.

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.

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."

Integration of single-nuclei RNA-sequencing, spatial transcriptomics and histochemistry defines the complex microenvironment of NF1-associated plexiform neurofibromas

Plexiform neurofibroma (PN) is a leading cause of morbidity in children with the genetic condition Neurofibromatosis Type 1 (NF1), often disfiguring or threatening vital structures. During formation of PN, a complex tumor microenvironment (TME) develops, with recruitment of neoplastic and non-neoplastic cell types being critical for growth and progression. Due to the cohesive cellularity of PN, single-cell RNA-sequencing is difficult and may result in a loss of detection of critical cellular subpopulations. To bypass this barrier, we performed single-nuclei RNA-sequencing (snRNA-seq) on 8 frozen PN samples, and integrated this with spatial transcriptomics (ST) in 4 PN samples and immunohistochemistry to provide morphological context to transcriptomic data. SnRNA-seq analysis definitively charted the heterogeneous cellular subpopulations in the PN TME, with the predominant fraction being fibroblast subtypes. PN showed a remarkable amount of inter-sample homogeneity regarding cellular subpopulation proportions despite being resected from a variety of anatomical locations. ST analysis identified distinct cellular subpopulations which were annotated using snRNA-seq data and correlated with histological features. Schwann cell/fibroblast interactions were identified by receptor/ligand interaction analysis demonstrating a high probability of Neurexin 1/Neuroligin 1 (NRXN1/NLGN1) receptor-ligand cross-talk predicted between fibroblasts and non-myelinated Schwann cells (NM-SC) and subtypes, respectively. We observed aberrant expression of NRXN1 and NLGN1 in our PN snRNA-seq data compared to a normal mouse sciatic nerve single-cell RNA-seq dataset. This pathway has never been described in PN and may indicate a clear and direct communication pathway between putative NM-SC cells of origin and surrounding fibroblasts, potentially driving disease progression. SnRNA-seq integrated with spatial transcriptomics advances our understanding of the complex cellular heterogeneity of PN TME and identify potential novel communication pathways that may drive disease progression, a finding that could provide translational therapy options for patients with these devastating tumors of childhood and early adulthood.

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.

Hyperglycemia increases SCO-spondin and Wnt5a secretion into the cerebrospinal fluid to regulate ependymal cell beating and glucose sensing

Hyperglycemia increases glucose concentrations in the cerebrospinal fluid (CSF), activating glucose-sensing mechanisms and feeding behavior in the hypothalamus. Here, we discuss how hyperglycemia temporarily modifies ependymal cell ciliary beating to increase hypothalamic glucose sensing. A high level of glucose in the rat CSF stimulates glucose transporter 2 (GLUT2)-positive subcommissural organ (SCO) cells to release SCO-spondin into the dorsal third ventricle. Genetic inactivation of mice GLUT2 decreases hyperglycemia-induced SCO-spondin secretion. In addition, SCO cells secrete Wnt5a-positive vesicles; thus, Wnt5a and SCO-spondin are found at the apex of dorsal ependymal cilia to regulate ciliary beating. Frizzled-2 and ROR2 receptors, as well as specific proteoglycans, such as glypican/testican (essential for the interaction of Wnt5a with its receptors) and Cx43 coupling, were also analyzed in ependymal cells. Finally, we propose that the SCO-spondin/Wnt5a/Frizzled-2/Cx43 axis in ependymal cells regulates ciliary beating, a cyclic and adaptive signaling mechanism to control glucose sensing.

Identification of HDAC4 as a potential therapeutic target and prognostic biomarker for ZFTA-fused ependymomas

Members of the HDAC family are predictive biomarkers and regulate the tumorigenesis in several cancers. However, the role of these genes in the biology of intracranial ependymomas (EPNs) remains unexplored. Here, an analysis of eighteen HDACs genes in an EPN transcriptomic dataset, revealed significantly higher levels of HDAC4 in supratentorial ZFTA fusion (ST-ZFTA) compared with ST-YAP1 fusion and posterior fossa EPNs, while HDAC7 and SIRT2 were downregulated in ST-ZFTA. HDAC4 was also overexpressed in ST-ZFTA as measured by single-cell RNA-Seq, quantitative real time-polymerase chain reaction, and immunohistochemistry. Survival analyses showed a significantly worse outcome for EPNs with higher HDAC4 and SIRT1 mRNA levels. Ontology enrichment analysis showed an HDAC4-high signature consistent with viral processes while collagen-containing extracellular matrix and cell-cell junction were enriched in those with an HDAC4-low signature. Immune gene analysis demonstrated a correlation between HDAC4 expression and low levels of NK resting cells. Several small molecules compounds targeting HDAC4 and ABCG2, were predicted by in silico analysis to be effective against HDAC4-high ZFTA. Our results provide novel insights into the biology of the HDAC family in intracranial ependymomas and reveal HDAC4 as a prognostic marker and potential therapeutic target in ST-ZFTA.

Characterizing the tumor immune microenvironment of ependymomas using targeted gene expression profiles and RNA sequencing

BACKGROUND

Defining the tumor immune microenvironment (TIME) of patients using transcriptome analysis is gaining more popularity. Here, we examined and discussed the pros and cons of using RNA sequencing for fresh frozen samples and targeted gene expression immune profiles (NanoString) for formalin-fixed, paraffin-embedded (FFPE) samples to characterize the TIME of ependymoma samples.

RESULTS

Our results showed a stable expression of the 40 housekeeping genes throughout all samples. The Pearson correlation of the endogenous genes was high. To define the TIME, we first checked the expression of the PTPRC gene, known as CD45, and found it was above the detection limit in all samples by both techniques. T cells were identified consistently using the two types of data. In addition, both techniques showed that the immune landscape was heterogeneous in the 6 ependymoma samples used for this study.

CONCLUSIONS

The low-abundant genes were detected in higher quantities using the NanoString technique, even when FFPE samples were used. RNA sequencing is better suited for biomarker discovery, fusion gene detection, and getting a broader overview of the TIME. The technique that was used to measure the samples had a considerable effect on the type of immune cells that were identified. The limited number of tumor-infiltrating immune cells compared to the high density of tumor cells in ependymoma can limit the sensitivity of RNA expression techniques regarding the identification of the infiltrating immune cells.

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.

Glial progenitor heterogeneity and key regulators revealed by single-cell RNA sequencing provide insight to regeneration in spinal cord injury

Recent studies have revealed the heterogeneous nature of astrocytes; however, how diverse constituents of astrocyte-lineage cells are regulated in adult spinal cord after injury and contribute to regeneration remains elusive. We perform single-cell RNA sequencing of GFAP-expressing cells from sub-chronic spinal cord injury models and identify and compare with the subpopulations in acute-stage data. We find subpopulations with distinct functional enrichment and their identities defined by subpopulation-specific transcription factors and regulons. Immunohistochemistry, RNAscope experiments, and quantification by stereology verify the molecular signature, location, and morphology of potential resident neural progenitors or neural stem cells in the adult spinal cord before and after injury and uncover the populations of the intermediate cells enriched in neuronal genes that could potentially transition into other subpopulations. This study has expanded the knowledge of the heterogeneity and cell state transition of glial progenitors in adult spinal cord before and after injury.

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.

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.

Common molecular features of H3K27M DMGs and PFA ependymomas map to hindbrain developmental pathways

Globally decreased histone 3, lysine 27 tri-methylation (H3K27me3) is a hallmark of H3K27-altered diffuse midline gliomas (DMGs) and group-A posterior fossa ependymomas (PFAs). H3K27-altered DMGs are largely characterized by lysine-to-methionine mutations in histone 3 at position 27 (H3K27M). Most PFAs overexpress EZH inhibitory protein (EZHIP), which possesses a region of similarity to the mutant H3K27M. Both H3K27M and EZHIP inhibit the function of the polycomb repressive complex 2 (PRC2) responsible for H3K27me3 deposition. These tumors often arise in neighboring regions of the brainstem and posterior fossa. In rare cases PFAs harbor H3K27M mutations, and DMGs overexpress EZHIP. These findings together raise the possibility that certain cell populations in the developing hindbrain/posterior fossa are especially sensitive to modulation of H3K27me3 states. We identified shared molecular features by comparing genomic, bulk transcriptomic, chromatin-based profiles, and single-cell RNA-sequencing (scRNA-seq) data from the two tumor classes. Our approach demonstrated that 1q gain, a key biomarker in PFAs, is prognostic in H3.1K27M, but not H3.3K27M gliomas. Conversely, Activin A Receptor Type 1 (ACVR1), which is associated with mutations in H3.1K27M gliomas, is overexpressed in a subset of PFAs with poor outcome. Despite diffuse H3K27me3 reduction, previous work shows that both tumors maintain genomic H3K27me3 deposition at select sites. We demonstrate heterogeneity in shared patterns of residual H3K27me3 for both tumors that largely segregated with inferred anatomic tumor origins and progenitor populations of tumor cells. In contrast, analysis of genes linked to H3K27 acetylation (H3K27ac)-marked enhancers showed higher expression in astrocytic-like tumor cells. Finally, common H3K27me3-marked genes mapped closely to expression patterns in the human developing hindbrain. Overall, our data demonstrate developmentally relevant molecular similarities between PFAs and H3K27M DMGs and support the overall hypothesis that deregulated mechanisms of hindbrain development are central to the biology of both tumors.

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.

Single-cell technologies uncover intra-tumor heterogeneity in childhood cancers

Childhood cancer is the second leading cause of death in children aged 1 to 14. Although survival rates have vastly improved over the past 40 years, cancer resistance and relapse remain a significant challenge. Advances in single-cell technologies enable dissection of tumors to unprecedented resolution. This facilitates unraveling the heterogeneity of childhood cancers to identify cell subtypes that are prone to treatment resistance. The rapid accumulation of single-cell data from different modalities necessitates the development of novel computational approaches for processing, visualizing, and analyzing single-cell data. Here, we review single-cell approaches utilized or under development in the context of childhood cancers. We review computational methods for analyzing single-cell data and discuss best practices for their application. Finally, we review the impact of several studies of childhood tumors analyzed with these approaches and future directions to implement single-cell studies into translational cancer research in pediatric oncology.

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.

Is IIIG9 a New Protein with Exclusive Ciliary Function? Analysis of Its Potential Role in Cancer and Other Pathologies

The identification of new proteins that regulate the function of one of the main cellular phosphatases, protein phosphatase 1 (PP1), is essential to find possible pharmacological targets to alter phosphatase function in various cellular processes, including the initiation and development of multiple diseases. IIIG9 is a regulatory subunit of PP1 initially identified in highly polarized ciliated cells. In addition to its ciliary location in ependymal cells, we recently showed that IIIG9 has extraciliary functions that regulate the integrity of adherens junctions. In this review, we perform a detailed analysis of the expression, localization, and function of IIIG9 in adult and developing normal brains. In addition, we provide a 3D model of IIIG9 protein structure for the first time, verifying that the classic structural and conformational characteristics of the PP1 regulatory subunits are maintained. Our review is especially focused on finding evidence linking IIIG9 dysfunction with the course of some pathologies, such as ciliopathies, drug dependence, diseases based on neurological development, and the development of specific high-malignancy and -frequency brain tumors in the pediatric population. Finally, we propose that IIIG9 is a relevant regulator of PP1 function in physiological and pathological processes in the CNS.

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.

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.

Single-Cell RNA Sequencing Unravels Upregulation of Immune Cell Crosstalk in Relapsed Pediatric Ependymoma

Ependymoma (EPN) is a malignant glial tumor occurring throughout central nervous system, which commonly presents in children. Although recent studies have characterized EPN samples at both the bulk and single-cell level, intratumoral heterogeneity across subclones remains a confounding factor that impedes understanding of EPN biology. In this study, we generated a high-resolution single-cell dataset of pediatric ependymoma with a particular focus on the comparison of subclone differences within tumors and showed upregulation of cilium-associated genes in more highly differentiated subclone populations. As a proxy to traditional pseudotime analysis, we applied a novel trajectory scoring method to reveal cellular compositions associated with poor survival outcomes across primary and relapsed patients. Furthermore, we identified putative cell–cell communication features between relapsed and primary samples and showed upregulation of pathways associated with immune cell crosstalk. Our results revealed both inter- and intratumoral heterogeneity in EPN and provided a framework for studying transcriptomic signatures of individual subclones at single-cell resolution.

Single cell transcriptomics reveals unique metabolic profiles of ependymoma subgroups

OBJECTIVE

Ependymomas are biologically diverse tumors with five major genomic subgroups. However, intratumor heterogeneity continues to be poorly understood. The present study characterized the metabolic landscapes of ependymoma subgroups at the single-cell level.

METHODS

Expression profiles from 11,200 ependymoma single cells derived from the five major subgroups and 7,200 ependymoma-derived non-neoplastic cells were computationally analyzed using a robust workflow to elucidate relative differences in metabolic pathway activities.

RESULTS

Dimensionality reduction using metabolic expression profiles exhibited clustering corresponding to each tumor subgroup, but non-neoplastic cells exhibited no discernable differences between subgroups. From the 80 metabolic pathways examined, over 75 pathways had significantly different activity scores between ependymoma subgroups. Further analysis of metabolic heterogeneity suggests that mitochondrial oxidative phosphorylation accounts for considerable metabolic variation within tumor subgroups and non-neoplastic cells of the same cell type. Drug metabolism pathways, specifically those involving cytochromes P450, were also found to be major contributors to heterogeneity.

CONCLUSIONS

Ependymoma subgroups display distinct metabolic differences as evaluated through gene expression profiles with certain pathways contributing greatly to intra-subgroup variation. These results may account for variation in tumor metabolism, treatment response, and potential targeting approaches that disrupt metabolic signalling.

The biology of ependymomas and emerging novel therapies

Ependymomas are rare central nervous system tumours that can arise in the brain's supratentorial region or posterior fossa, or in the spinal cord. In 1924, Percival Bailey published the first comprehensive study of ependymomas. Since then, and especially over the past 10 years, our understanding of ependymomas has grown exponentially. In this Review, we discuss the evolution in knowledge regarding ependymoma subgroups and the resultant clinical implications. We also discuss key oncogenic and tumour suppressor signalling pathways that regulate tumour growth, the role of epigenetic dysregulation in the biology of ependymomas, and the various biological features of ependymoma tumorigenesis, including cell immortalization, stem cell-like properties, the tumour microenvironment and metastasis. We further review the limitations of current therapies such as relapse, radiation-induced cognitive deficits and chemotherapy resistance. Finally, we highlight next-generation therapies that are actively being explored, including tyrosine kinase inhibitors, telomerase inhibitors, anti-angiogenesis agents and immunotherapy.

Neoplastic and immune single-cell transcriptomics define subgroup-specific intra-tumoral heterogeneity of childhood medulloblastoma

BACKGROUND

Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. We aimed to determine the influence of neoplastic and immune cell diversity on MB biology in patient samples and animal models.

METHODS

To better characterize cellular heterogeneity in MB we used single-cell RNA sequencing, immunohistochemistry, and deconvolution of transcriptomic data to profile neoplastic and immune populations in patient samples and animal models across childhood MB subgroups.

RESULTS

Neoplastic cells cluster primarily according to individual sample of origin which is influenced by chromosomal copy number variance. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes, including photoreceptor and glutamatergic neuron-like cells in molecular subgroups GP3 and GP4, and a specific nodule-associated neuronally differentiated subpopulation in the sonic hedgehog subgroup. We definitively chart the spectrum of MB immune cell infiltrates, which include subpopulations that recapitulate developmentally related neuron-pruning and antigen-presenting myeloid cells. MB cellular diversity matching human samples is mirrored in subgroup-specific mouse models of MB.

CONCLUSIONS

These findings provide a clearer understanding of the diverse neoplastic and immune cell subpopulations that constitute the MB microenvironment.

Cell-of-Origin and Genetic, Epigenetic, and Microenvironmental Factors Contribute to the Intra-Tumoral Heterogeneity of Pediatric Intracranial Ependymoma

Intra-tumoral heterogeneity (ITH) is a complex multifaceted phenomenon that posits major challenges for the clinical management of cancer patients. Genetic, epigenetic, and microenvironmental factors are concurrent drivers of diversity among the distinct populations of cancer cells. ITH may also be installed by cancer stem cells (CSCs), that foster unidirectional hierarchy of cellular phenotypes or, alternatively, shift dynamically between distinct cellular states. Ependymoma (EPN), a molecularly heterogeneous group of tumors, shows a specific spatiotemporal distribution that suggests a link between ependymomagenesis and alterations of the biological processes involved in embryonic brain development. In children, EPN most often arises intra-cranially and is associated with an adverse outcome. Emerging evidence shows that EPN displays large intra-patient heterogeneity. In this review, after touching on EPN inter-tumoral heterogeneity, we focus on the sources of ITH in pediatric intra-cranial EPN in the framework of the CSC paradigm. We also examine how single-cell technology has shed new light on the complexity and developmental origins of EPN and the potential impact that this understanding may have on the therapeutic strategies against this deadly pediatric malignancy.

Interrogation of the microenvironmental landscape in spinal ependymomas reveals dual functions of tumor-associated macrophages

Spinal ependymomas are the most common spinal cord tumors in adults, but their intratumoral cellular heterogeneity has been less studied, and how spinal microglia are involved in tumor progression is still unknown. Here, our single-cell RNA-sequencing analyses of three spinal ependymoma subtypes dissect the microenvironmental landscape of spinal ependymomas and reveal tumor-associated macrophage (TAM) subsets with distinct functional phenotypes. CCL2+ TAMs are related to the immune response and exhibit a high capacity for apoptosis, while CD44+ TAMs are associated with tumor angiogenesis. By combining these results with those of single-cell ATAC-sequencing data analysis, we reveal that TEAD1 and EGR3 play roles in regulating the functional diversity of TAMs. We further identify diverse characteristics of both malignant cells and TAMs that might underlie the different malignant degrees of each subtype. Finally, assessment of cell-cell interactions reveal that stromal cells act as extracellular factors that mediate TAM diversity. Overall, our results reveal dual functions of TAMs in tumor progression, providing valuable insights for TAM-targeting immunotherapy.

Targeting integrated epigenetic and metabolic pathways in lethal childhood PFA ependymomas

Childhood posterior fossa group A ependymomas (PFAs) have limited treatment options and bear dismal prognoses compared to group B ependymomas (PFBs). PFAs overexpress the oncohistone-like protein EZHIP (enhancer of Zeste homologs inhibitory protein), causing global reduction of repressive histone H3 lysine 27 trimethylation (H3K27me3), similar to the oncohistone H3K27M. Integrated metabolic analyses in patient-derived cells and tumors, single-cell RNA sequencing of tumors, and noninvasive metabolic imaging in patients demonstrated enhanced glycolysis and tricarboxylic acid (TCA) cycle metabolism in PFAs. Furthermore, high glycolytic gene expression in PFAs was associated with a poor outcome. PFAs demonstrated high EZHIP expression associated with poor prognosis and elevated activating mark histone H3 lysine 27 acetylation (H3K27ac). Genomic H3K27ac was enriched in PFAs at key glycolytic and TCA cycle–related genes including hexokinase-2 and pyruvate dehydrogenase. Similarly, mouse neuronal stem cells (NSCs) expressing wild-type EZHIP (EZHIP-WT) versus catalytically attenuated EZHIP-M406K demonstrated H3K27ac enrichment at hexokinase-2 and pyruvate dehydrogenase, accompanied by enhanced glycolysis and TCA cycle metabolism. AMPKα-2, a key component of the metabolic regulator AMP-activated protein kinase (AMPK), also showed H3K27ac enrichment in PFAs and EZHIP-WT NSCs. The AMPK activator metformin lowered EZHIP protein concentrations, increased H3K27me3, suppressed TCA cycle metabolism, and showed therapeutic efficacy in vitro and in vivo in patient-derived PFA xenografts in mice. Our data indicate that PFAs and EZHIP-WT–expressing NSCs are characterized by enhanced glycolysis and TCA cycle metabolism. Repurposing the antidiabetic drug metformin lowered pathogenic EZHIP, increased H3K27me3, and suppressed tumor growth, suggesting that targeting integrated metabolic/epigenetic pathways is a potential therapeutic strategy for treating childhood ependymomas.

Targeting fibroblast growth factor receptors to combat aggressive ependymoma

Ependymomas (EPN) are central nervous system tumors comprising both aggressive and more benign molecular subtypes. However, therapy of the high-risk subtypes posterior fossa group A (PF-A) and supratentorial RELA-fusion positive (ST-RELA) is limited to gross total resection and radiotherapy, as effective systemic treatment concepts are still lacking. We have recently described fibroblast growth factor receptors 1 and 3 (FGFR1/FGFR3) as oncogenic drivers of EPN. However, the underlying molecular mechanisms and their potential as therapeutic targets have not yet been investigated in detail. Making use of transcriptomic data across 467 EPN tissues, we found that FGFR1 and FGFR3 were both widely expressed across all molecular groups. FGFR3 mRNA levels were enriched in ST-RELA showing the highest expression among EPN as well as other brain tumors. We further identified high expression levels of fibroblast growth factor 1 and 2 (FGF1, FGF2) across all EPN subtypes while FGF9 was elevated in ST-EPN. Interrogation of our EPN single-cell RNA-sequencing data revealed that FGFR3 was further enriched in cycling and progenitor-like cell populations. Corroboratively, we found FGFR3 to be predominantly expressed in radial glia cells in both mouse embryonal and human brain datasets. Moreover, we detected alternative splicing of the FGFR1/3-IIIc variant, which is known to enhance ligand affinity and FGFR signaling. Dominant-negative interruption of FGFR1/3 activation in PF-A and ST-RELA cell models demonstrated inhibition of key oncogenic pathways leading to reduced cell growth and stem cell characteristics. To explore the feasibility of therapeutically targeting FGFR, we tested a panel of FGFR inhibitors in 12 patient-derived EPN cell models revealing sensitivity in the low-micromolar to nano-molar range. Finally, we gain the first clinical evidence for the activity of the FGFR inhibitor nintedanib in the treatment of a patient with recurrent ST-RELA. Together, these preclinical and clinical data suggest FGFR inhibition as a novel and feasible approach to combat aggressive EPN.

C11orf95-RELA fusion drives aberrant gene expression through the unique epigenetic regulation for ependymoma formation

Recurrent C11orf95-RELA fusions (RELA^FUS) are the hallmark of supratentorial ependymomas. The presence of RELA as the fusion partner indicates a close association of aberrant NF-κB activity with tumorigenesis. However, the oncogenic role of the C11orf95 has not been determined. Here, we performed ChIP-seq analyses to explore genomic regions bound by RELA^FUS and H3K27ac proteins in human 293T and mouse ependymoma cells. We then utilized published RNA-Seq data from human and mouse RELA^FUS tumors and identified target genes that were directly regulated by RELA^FUS in these tumors. Subsequent transcription factor motif analyses of RELA^FUS target genes detected a unique GC-rich motif recognized by the C11orf95 moiety, that is present in approximately half of RELA^FUS target genes. Luciferase assays confirmed that a promoter carrying this motif is sufficient to drive RELA^FUS-dependent gene expression. Further, the RELA^FUS target genes were found to be overlapped with Rela target genes primarily via non-canonical NF-κB binding sites. Using a series of truncation and substitution mutants of RELA^FUS, we also show that the activation domain in the RELA^FUS moiety is necessary for the regulation of gene expression of these RELA^FUS target genes. Lastly, we performed an anti-cancer drug screening with mouse ependymoma cells and identified potential anti-ependymoma drugs that are related to the oncogenic mechanism of RELA^FUS. These findings suggested that RELA^FUS might induce ependymoma formation through oncogenic pathways orchestrated by both C11orf95 and RELA target genes. Thus, our study unveils a complex gene function of RELA^FUS as an oncogenic transcription factor in RELA^FUS positive ependymomas.

The BET Inhibitor OTX015 Exhibits In Vitro and In Vivo Antitumor Activity in Pediatric Ependymoma Stem Cell Models

Childhood ependymomas are heterogenous chemoresistant neoplasms arising from aberrant stem-like cells. Epigenome deregulation plays a pivotal role in ependymoma pathogenesis, suggesting that epigenetic modifiers hold therapeutic promise against this disease. Bromodomain and extraterminal domain (BET) proteins are epigenome readers of acetylated signals in histones and coactivators for oncogenic and stemness-related transcriptional networks, including MYC/MYCN (Proto-Oncogene, BHLH Transcritpion Factor)-regulated genes. We explored BET inhibition as an anticancer strategy in a panel of pediatric patient-derived ependymoma stem cell models by OTX015-mediated suppression of BET/acetylated histone binding. We found that ependymoma tissues and lines express BET proteins and their targets MYC and MYCN. In vitro, OTX015 reduced cell proliferation by inducing G0/G1-phase accumulation and apoptosis at clinically tolerable doses. Mechanistically, inhibitory p21 and p27 increased in a p53-independent manner, whereas the proliferative driver, phospho-signal transducer and activator of transcription 3 (STAT3), decreased. Upregulation of apoptosis-related proteins and survivin downregulation were correlated with cell line drug sensitivity. Minor alterations of MYC/MYCN expression were reported. In vivo, OTX015 significantly improved survival in 2/3 orthotopic ependymoma models. BET proteins represent promising targets for pharmaceutical intervention with OTX015 against ependymoma. The identification of predictive determinants of sensitivity may help identify ependymoma molecular subsets more likely to benefit from BET inhibitor therapies.

Decreased content of ascorbic acid (vitamin C) in the brain of knockout mouse models of Na+,K+-ATPase-related neurologic disorders

Na⁺,K⁺-ATPase is a crucial protein responsible for maintaining the electrochemical gradients across the cell membrane. The Na⁺,K⁺-ATPase is comprised of catalytic α, β, and γ subunits. In adult brains, the α3 subunit, encoded by ATP1A3, is predominantly expressed in neurons, whereas the α2 subunit, encoded by ATP1A2, is expressed in glial cells. In foetal brains, the α2 is expressed in neurons as well. Mutations in α subunits cause a variety of neurologic disorders. Notably, the onset of symptoms in ATP1A2- and ATP1A3-related neurologic disorders is usually triggered by physiological or psychological stressors. To gain insight into the distinct roles of the α2 and α3 subunits in the developing foetal brain, whose developmental dysfunction may be a predisposing factor of neurologic disorders, we compared the phenotypes of mouse foetuses with double homozygous knockout of Atp1a2 and Atp1a3 (α2α3-dKO) to those with single knockout. The brain haemorrhage phenotype of α2α3-dKO was similar to that of homozygous knockout of the gene encoding ascorbic acid (ASC or vitamin C) transporter, SVCT2. The α2α3-dKO brain showed significantly decreased level of ASC compared with the wild-type (WT) and single knockout. We found that the ASC content in the basal ganglia and cerebellum was significantly lower in the adult Atp1a3 heterozygous knockout mouse (α3-HT) than in the WT. Interestingly, we observed a significant decrease in the ASC level in the basal ganglia and cerebellum of α3-HT in the peripartum period, during which mice are under physiological stress. These observations indicate that the α2 and α3 subunits independently contribute to the ASC level in the foetal brain and that the α3 subunit contributes to ASC transport in the adult basal ganglia and cerebellum. We propose that decreases in ASC levels may affect neural network development and are linked to the pathophysiology of ATP1A2- and ATP1A3-related neurologic disorders.