A single-cell transcriptomic map of the developing Atoh1 lineage identifies neural fate decisions and neuronal diversity in the hindbrain.

Proneural transcription factors establish molecular cascades to orchestrate neuronal diversity. One such transcription factor, Atonal homolog 1 (Atoh1), gives rise to cerebellar excitatory neurons and over 30 distinct nuclei in the brainstem critical for hearing, breathing, and balance. Although Atoh1 lineage neurons have been qualitatively described, the transcriptional programs that drive their fate decisions and the full extent of their diversity remain unknown. Here, we analyzed single-cell RNA sequencing and ATOH1 DNA binding in Atoh1 lineage neurons of the developing mouse hindbrain. This high-resolution dataset identified markers for specific brainstem nuclei and demonstrated that transcriptionally heterogeneous progenitors require ATOH1 for proper migration. Moreover, we identified a sizable population of proliferating unipolar brush cell progenitors in the mouse Atoh1 lineage, previously described in humans as the origin of one medulloblastoma subtype. Collectively, our data provide insights into the developing mouse hindbrain and markers for functional assessment of understudied neuronal populations.

Early rhombic lip Protogenin+ve stem cells in a human-specific neurovascular niche initiate and maintain group 3 medulloblastoma.

We identify a population of Protogenin-positive (PRTG+ve) MYChigh NESTINlow stem cells in the four-week-old human embryonic hindbrain that subsequently localizes to the ventricular zone of the rhombic lip (RLVZ). Oncogenic transformation of early Prtg+ve rhombic lip stem cells initiates group 3 medulloblastoma (Gr3-MB)-like tumors. PRTG+ve stem cells grow adjacent to a human-specific interposed vascular plexus in the RLVZ, a phenotype that is recapitulated in Gr3-MB but not in other types of medulloblastoma. Co-culture of Gr3-MB with endothelial cells promotes tumor stem cell growth, with the endothelial cells adopting an immature phenotype. Targeting the PRTGhigh compartment of Gr3-MB in vivo using either the diphtheria toxin system or chimeric antigen receptor T cells constitutes effective therapy. Human Gr3-MBs likely arise from early embryonic RLVZ PRTG+ve stem cells inhabiting a specific perivascular niche. Targeting the PRTGhigh compartment and/or the perivascular niche represents an approach to treat children with Gr3-MB.

A group 3 medulloblastoma stem cell program is maintained by OTX2-mediated alternative splicing.

OTX2 is a transcription factor and known driver in medulloblastoma (MB), where it is amplified in a subset of tumours and overexpressed in most cases of group 3 and group 4 MB. Here we demonstrate a noncanonical role for OTX2 in group 3 MB alternative splicing. OTX2 associates with the large assembly of splicing regulators complex through protein-protein interactions and regulates a stem cell splicing program. OTX2 can directly or indirectly bind RNA and this may be partially independent of its DNA regulatory functions. OTX2 controls a pro-tumorigenic splicing program that is mirrored in human cerebellar rhombic lip origins. Among the OTX2-regulated differentially spliced genes, PPHLN1 is expressed in the most primitive rhombic lip stem cells, and targeting PPHLN1 splicing reduces tumour growth and enhances survival in vivo. These findings identify OTX2-mediated alternative splicing as a major determinant of cell fate decisions that drive group 3 MB progression.

A therapeutically targetable positive feedback loop between lnc-HLX-2-7, HLX, and MYC that promotes group 3 medulloblastoma.

Recent studies suggest that long non-coding RNAs (lncRNAs) contribute to medulloblastoma (MB) formation and progression. We have identified an lncRNA, lnc-HLX-2-7, as a potential therapeutic target in group 3 (G3) MBs. lnc-HLX-2-7 RNA specifically accumulates in the promoter region of HLX, a sense-overlapping gene of lnc-HLX-2-7, which activates HLX expression by recruiting multiple factors, including enhancer elements. RNA sequencing and chromatin immunoprecipitation reveal that HLX binds to and activates the promoters of several oncogenes, including TBX2, LIN9, HOXM1, and MYC. Intravenous treatment with cerium-oxide-nanoparticle-coated antisense oligonucleotides targeting lnc-HLX-2-7 (CNP-lnc-HLX-2-7) inhibits tumor growth by 40%-50% in an intracranial MB xenograft mouse model. Combining CNP-lnc-HLX-2-7 with standard-of-care cisplatin further inhibits tumor growth and significantly prolongs mouse survival compared with CNP-lnc-HLX-2-7 monotherapy. Thus, the lnc-HLX-2-7-HLX-MYC axis is important for regulating G3 MB progression, providing a strong rationale for using lnc-HLX-2-7 as a therapeutic target for G3 MBs.

Lymphatic endothelial-like cells promote glioblastoma stem cell growth through cytokine-driven cholesterol metabolism.

Glioblastoma is the most lethal primary brain tumor with glioblastoma stem cells (GSCs) atop a cellular hierarchy. GSCs often reside in a perivascular niche, where they receive maintenance cues from endothelial cells, but the role of heterogeneous endothelial cell populations remains unresolved. Here, we show that lymphatic endothelial-like cells (LECs), while previously unrecognized in brain parenchyma, are present in glioblastomas and promote growth of CCR7-positive GSCs through CCL21 secretion. Disruption of CCL21-CCR7 paracrine communication between LECs and GSCs inhibited GSC proliferation and growth. LEC-derived CCL21 induced KAT5-mediated acetylation of HMGCS1 on K273 in GSCs to enhance HMGCS1 protein stability. HMGCS1 promoted cholesterol synthesis in GSCs, favorable for tumor growth. Expression of the CCL21-CCR7 axis correlated with KAT5 expression and HMGCS1K273 acetylation in glioblastoma specimens, informing patient outcome. Collectively, glioblastomas contain previously unrecognized LECs that promote the molecular crosstalk between endothelial and tumor cells, offering potentially alternative therapeutic strategies.

miRNA-211 maintains metabolic homeostasis in medulloblastoma through its target gene long-chain acyl-CoA synthetase 4.

The prognosis of childhood medulloblastoma (MB) is often poor, and it usually requires aggressive therapy that adversely affects quality of life. microRNA-211 (miR-211) was previously identified as an important regulator of cells that descend from neural cells. Since medulloblastomas primarily affect cells with similar ontogeny, we investigated the role and mechanism of miR-211 in MB. Here we showed that miR-211 expression was highly downregulated in cell lines, PDXs, and clinical samples of different MB subgroups (SHH, Group 3, and Group 4) compared to normal cerebellum. miR-211 gene was ectopically expressed in transgenic cells from MB subgroups, and they were subjected to molecular and phenotypic investigations. Monoclonal cells stably expressing miR-211 were injected into the mouse cerebellum. miR-211 forced expression acts as a tumor suppressor in MB both in vitro and in vivo, attenuating growth, promoting apoptosis, and inhibiting invasion. In support of emerging regulatory roles of metabolism in various forms of cancer, we identified the acyl-CoA synthetase long-chain family member (ACSL4) as a direct miR-211 target. Furthermore, lipid nanoparticle-coated, dendrimer-coated, and cerium oxide-coated miR-211 nanoparticles were applied to deliver synthetic miR-211 into MB cell lines and cellular responses were assayed. Synthesizing nanoparticle-miR-211 conjugates can suppress MB cell viability and invasion in vitro. Our findings reveal miR-211 as a tumor suppressor and a potential therapeutic agent in MB. This proof-of-concept paves the way for further pre-clinical and clinical development.

Hedgehog target genes regulate lipid metabolism to drive basal cell carcinoma and medulloblastoma.

Hedgehog (Hh) signaling is essential for development, homeostasis, and regeneration1. Misactivation of the Hh pathway underlies medulloblastoma, the most common malignant brain tumor in children, and basal cell carcinoma (BCC), the most common cancer in the United States2. Primary cilia regulate Hh signal transduction3, but target genes that drive cell fate decisions in response to ciliary ligands or oncogenic Hh signaling are incompletely understood. Here we define the Hh gene expression program using RNA sequencing of cultured cells treated with ciliary ligands, BCCs from humans, and Hh-associated medulloblastomas from humans and mice (Fig. 1a). To validate our results, we integrate lipidomic mass spectrometry and bacterial metabolite labeling of free sterols with genetic and pharmacologic approaches in cells and mice. Our results reveal novel Hh target genes such as the oxysterol synthase Hsd11ß1 and the adipokine Retnla that regulate lipid metabolism to drive cell fate decisions in response to Hh pathway activation. These data provide insights into cellular mechanisms underlying ciliary and oncogenic Hh signaling and elucidate targets to treat Hh-associated cancers.

The LIN28B-let-7-PBK pathway is essential for group 3 medulloblastoma tumor growth and survival.

Children with Group 3 medulloblastoma (G3 MB) have a very poor prognosis, and many do not survive beyond 5 years after diagnosis. A factor that may contribute to this is the lack of available targeted therapy. Expression of protein lin-28 homolog B (LIN28B), a regulator of developmental timing, is upregulated in several cancers, including G3 MB, and is associated with worse survival in this disease. Here, we investigate the role of the LIN28B pathway in G3 MB and demonstrate that the LIN28B-lethal-7 (let-7; a microRNA that is a tumor suppressor)-lymphokine-activated killer T-cell-originated protein kinase (PBK; also known as PDZ-binding kinase) axis promotes G3 MB proliferation. LIN28B knockdown in G3-MB-patient-derived cell lines leads to a significant reduction in cell viability and proliferation in vitro and in prolonged survival of mice with orthotopic tumors. The LIN28 inhibitor N-methyl-N-[3-(3-methyl-1,2,4-triazolo[4,3-b]pyridazin-6-yl)phenyl]acetamide (1632) significantly reduces G3 MB cell growth and demonstrates efficacy in reducing tumor growth in mouse xenograft models. Inhibiting PBK using HI-TOPK-032 also results in a significant reduction in G3 MB cell viability and proliferation. Together, these results highlight a critical role for the LIN28B-let-7-PBK pathway in G3 MB and provide preliminary preclinical results for drugs targeting this pathway.

A neurodevelopmental epigenetic programme mediated by SMARCD3-DAB1-Reelin signalling is hijacked to promote medulloblastoma metastasis.

How abnormal neurodevelopment relates to the tumour aggressiveness of medulloblastoma (MB), the most common type of embryonal tumour, remains elusive. Here we uncover a neurodevelopmental epigenomic programme that is hijacked to induce MB metastatic dissemination. Unsupervised analyses of integrated publicly available datasets with our newly generated data reveal that SMARCD3 (also known as BAF60C) regulates Disabled 1 (DAB1)-mediated Reelin signalling in Purkinje cell migration and MB metastasis by orchestrating cis-regulatory elements at the DAB1 locus. We further identify that a core set of transcription factors, enhancer of zeste homologue 2 (EZH2) and nuclear factor I X (NFIX), coordinates with the cis-regulatory elements at the SMARCD3 locus to form a chromatin hub to control SMARCD3 expression in the developing cerebellum and in metastatic MB. Increased SMARCD3 expression activates Reelin-DAB1-mediated Src kinase signalling, which results in a MB response to Src inhibition. These data deepen our understanding of how neurodevelopmental programming influences disease progression and provide a potential therapeutic option for patients with MB.

Targeting Microglial Metabolic Rewiring Synergizes with Immune-Checkpoint Blockade Therapy for Glioblastoma.

Glioblastoma (GBM) constitutes the most lethal primary brain tumor for which immunotherapy has provided limited benefit. The unique brain immune landscape is reflected in a complex tumor immune microenvironment (TIME) in GBM. Here, single-cell sequencing of the GBM TIME revealed that microglia were under severe oxidative stress, which induced nuclear receptor subfamily 4 group A member 2 (NR4A2)-dependent transcriptional activity in microglia. Heterozygous Nr4a2 (Nr4a2+/-) or CX3CR1+ myeloid cell-specific Nr4a2 (Nr4a2fl/flCx3cr1Cre) genetic targeting reshaped microglia plasticity in vivo by reducing alternatively activated microglia and enhancing antigen presentation capacity for CD8+ T cells in GBM. In microglia, NR4A2 activated squalene monooxygenase (SQLE) to dysregulate cholesterol homeostasis. Pharmacologic NR4A2 inhibition attenuated the protumorigenic TIME, and targeting the NR4A2 or SQLE enhanced the therapeutic efficacy of immune-checkpoint blockade in vivo. Collectively, oxidative stress promotes tumor growth through NR4A2-SQLE activity in microglia, informing novel immune therapy paradigms in brain cancer. SIGNIFICANCE: Metabolic reprogramming of microglia in GBM informs synergistic vulnerabilities for immune-checkpoint blockade therapy in this immunologically cold brain tumor. This article is highlighted in the In This Issue feature, p. 799.

Failure of human rhombic lip differentiation underlies medulloblastoma formation.

Medulloblastoma (MB) comprises a group of heterogeneous paediatric embryonal neoplasms of the hindbrain with strong links to early development of the hindbrain1-4. Mutations that activate Sonic hedgehog signalling lead to Sonic hedgehog MB in the upper rhombic lip (RL) granule cell lineage5-8. By contrast, mutations that activate WNT signalling lead to WNT MB in the lower RL9,10. However, little is known about the more commonly occurring group 4 (G4) MB, which is thought to arise in the unipolar brush cell lineage3,4. Here we demonstrate that somatic mutations that cause G4 MB converge on the core binding factor alpha (CBFA) complex and mutually exclusive alterations that affect CBFA2T2, CBFA2T3, PRDM6, UTX and OTX2. CBFA2T2 is expressed early in the progenitor cells of the cerebellar RL subventricular zone in Homo sapiens, and G4 MB transcriptionally resembles these progenitors but are stalled in developmental time. Knockdown of OTX2 in model systems relieves this differentiation blockade, which allows MB cells to spontaneously proceed along normal developmental differentiation trajectories. The specific nature of the split human RL, which is destined to generate most of the neurons in the human brain, and its high level of susceptible EOMES+KI67+ unipolar brush cell progenitor cells probably predisposes our species to the development of G4 MB.

Combined MEK and JAK/STAT3 pathway inhibition effectively decreases SHH medulloblastoma tumor progression.

Medulloblastoma (MB) is the most common primary malignant pediatric brain cancer. We recently identified novel roles for the MEK/MAPK pathway in regulating human Sonic Hedgehog (SHH) MB tumorigenesis. The MEK inhibitor, selumetinib, decreased SHH MB growth while extending survival in mouse models. However, the treated mice ultimately succumbed to disease progression. Here, we perform RNA sequencing on selumetinib-treated orthotopic xenografts to identify molecular pathways that compensate for MEK inhibition specifically in vivo. Notably, the JAK/STAT3 pathway exhibits increased activation in selumetinib-treated tumors. The combination of selumetinib and the JAK/STAT3 pathway inhibitor, pacritinib, further reduces growth in two xenograft models and also enhances survival. Multiplex spatial profiling of proteins in drug-treated xenografts reveals shifted molecular dependencies and compensatory changes following combination drug treatment. Our study warrants further investigation into MEK and JAK/STAT3 inhibition as a novel combinatory therapeutic strategy for SHH MB.

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions.

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion.

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

An OTX2-PAX3 signaling axis regulates Group 3 medulloblastoma cell fate.

OTX2 is a potent oncogene that promotes tumor growth in Group 3 medulloblastoma. However, the mechanisms by which OTX2 represses neural differentiation are not well characterized. Here, we perform extensive multiomic analyses to identify an OTX2 regulatory network that controls Group 3 medulloblastoma cell fate. OTX2 silencing modulates the repressive chromatin landscape, decreases levels of PRC2 complex genes and increases the expression of neurodevelopmental transcription factors including PAX3 and PAX6. Expression of PAX3 and PAX6 is significantly lower in Group 3 medulloblastoma patients and is correlated with reduced survival, yet only PAX3 inhibits self-renewal in vitro and increases survival in vivo. Single cell RNA sequencing of Group 3 medulloblastoma tumorspheres demonstrates expression of an undifferentiated progenitor program observed in primary tumors and characterized by translation/elongation factor genes. Identification of mTORC1 signaling as a downstream effector of OTX2-PAX3 reveals roles for protein synthesis pathways in regulating Group 3 medulloblastoma pathogenesis.

Metabolic Regulation of the Epigenome Drives Lethal Infantile Ependymoma.

Posterior fossa A (PFA) ependymomas are lethal malignancies of the hindbrain in infants and toddlers. Lacking highly recurrent somatic mutations, PFA ependymomas are proposed to be epigenetically driven tumors for which model systems are lacking. Here we demonstrate that PFA ependymomas are maintained under hypoxia, associated with restricted availability of specific metabolites to diminish histone methylation, and increase histone demethylation and acetylation at histone 3 lysine 27 (H3K27). PFA ependymomas initiate from a cell lineage in the first trimester of human development that resides in restricted oxygen. Unlike other ependymomas, transient exposure of PFA cells to ambient oxygen induces irreversible cellular toxicity. PFA tumors exhibit a low basal level of H3K27me3, and, paradoxically, inhibition of H3K27 methylation specifically disrupts PFA tumor growth. Targeting metabolism and/or the epigenome presents a unique opportunity for rational therapy for infants with PFA ependymoma.

Transcriptional Programs Define Intratumoral Heterogeneity of Ewing Sarcoma at Single-Cell Resolution.

EWSR1-FLI1, the chimeric oncogene specific for Ewing sarcoma (EwS), induces a cascade of signaling events leading to cell transformation. However, it remains elusive how genetically homogeneous EwS cells can drive the heterogeneity of transcriptional programs. Here, we combine independent component analysis of single-cell RNA sequencing data from diverse cell types and model systems with time-resolved mapping of EWSR1-FLI1 binding sites and of open chromatin regions to characterize dynamic cellular processes associated with EWSR1-FLI1 activity. We thus define an exquisitely specific and direct enhancer-driven EWSR1-FLI1 program. In EwS tumors, cell proliferation and strong oxidative phosphorylation metabolism are associated with a well-defined range of EWSR1-FLI1 activity. In contrast, a subpopulation of cells from below and above the intermediary EWSR1-FLI1 activity is characterized by increased hypoxia. Overall, our study reveals sources of intratumoral heterogeneity within EwS tumors.

ZRANB2 and SYF2-mediated splicing programs converging on ECT2 are involved in breast cancer cell resistance to doxorubicin.

Besides analyses of specific alternative splicing (AS) variants, little is known about AS regulatory pathways and programs involved in anticancer drug resistance. Doxorubicin is widely used in breast cancer chemotherapy. Here, we identified 1723 AS events and 41 splicing factors regulated in a breast cancer cell model of acquired resistance to doxorubicin. An RNAi screen on splicing factors identified the little studied ZRANB2 and SYF2, whose depletion partially reversed doxorubicin resistance. By RNAi and RNA-seq in resistant cells, we found that the AS programs controlled by ZRANB2 and SYF2 were enriched in resistance-associated AS events, and converged on the ECT2 splice variant including exon 5 (ECT2-Ex5+). Both ZRANB2 and SYF2 were found associated with ECT2 pre-messenger RNA, and ECT2-Ex5+ isoform depletion reduced doxorubicin resistance. Following doxorubicin treatment, resistant cells accumulated in S phase, which partially depended on ZRANB2, SYF2 and the ECT2-Ex5+ isoform. Finally, doxorubicin combination with an oligonucleotide inhibiting ECT2-Ex5 inclusion reduced doxorubicin-resistant tumor growth in mouse xenografts, and high ECT2-Ex5 inclusion levels were associated with bad prognosis in breast cancer treated with chemotherapy. Altogether, our data identify AS programs controlled by ZRANB2 and SYF2 and converging on ECT2, that participate to breast cancer cell resistance to doxorubicin.