Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis.

Tumor-educated blood platelets (TEPs) are implicated as central players in the systemic and local responses to tumor growth, thereby altering their RNA profile. We determined the diagnostic potential of TEPs by mRNA sequencing of 283 platelet samples. We distinguished 228 patients with localized and metastasized tumors from 55 healthy individuals with 96% accuracy. Across six different tumor types, the location of the primary tumor was correctly identified with 71% accuracy. Also, MET or HER2-positive, and mutant KRAS, EGFR, or PIK3CA tumors were accurately distinguished using surrogate TEP mRNA profiles. Our results indicate that blood platelets provide a valuable platform for pan-cancer, multiclass cancer, and companion diagnostics, possibly enabling clinical advances in blood-based "liquid biopsies".

In vivo functional investigation of oncogenes using somatic gene transfer has been successfully exploited to validate their role in tumorigenesis. For tumour suppressor genes this has proven more challenging due to technical aspects. To provide a flexible and effective method for investigating somatic loss-of-function alterations and their influence on tumorigenesis, we have established CRISPR/Cas9-mediated somatic gene disruption, allowing for in vivo targeting of TSGs. Here we demonstrate the utility of this approach by deleting single (Ptch1) or multiple genes (Trp53, Pten, Nf1) in the mouse brain, resulting in the development of medulloblastoma and glioblastoma, respectively. Using whole-genome sequencing (WGS) we characterized the medulloblastoma-driving Ptch1 deletions in detail and show that no off-targets were detected in these tumours. This method provides a fast and convenient system for validating the emerging wealth of novel candidate tumour suppressor genes and the generation of faithful animal models of human cancer.

High-risk neuroblastoma (NB) is responsible for a disproportionate number of childhood deaths due to cancer. One indicator of high-risk NB is amplification of the neural MYC (MYCN) oncogene, which is currently therapeutically intractable. Identification of anaplastic lymphoma kinase (ALK) as an NB oncogene raised the possibility of using ALK tyrosine kinase inhibitors (TKIs) in treatment of patients with activating ALK mutations. 8-10% of primary NB patients are ALK-positive, a figure that increases in the relapsed population. ALK is activated by the ALKAL2 ligand located on chromosome 2p, along with ALK and MYCN, in the "2p-gain" region associated with NB. Dysregulation of ALK ligand in NB has not been addressed, although one of the first oncogenes described was v-sis that shares > 90% homology with PDGF. Therefore, we tested whether ALKAL2 ligand could potentiate NB progression in the absence of ALK mutation. We show that ALKAL2 overexpression in mice drives ALK TKI-sensitive NB in the absence of ALK mutation, suggesting that additional NB patients, such as those exhibiting 2p-gain, may benefit from ALK TKI-based therapeutic intervention.

More than 60% of supratentorial ependymomas harbor a ZFTA-RELA (ZRfus) gene fusion (formerly C11orf95-RELA). To study the biology of ZRfus, we developed an autochthonous mouse tumor model using in utero electroporation (IUE) of the embryonic mouse brain. Integrative epigenomic and transcriptomic mapping was performed on IUE-driven ZRfus tumors by CUT&RUN, chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin sequencing, and RNA sequencing and compared with human ZRfus-driven ependymoma. In addition to direct canonical NFκB pathway activation, ZRfus dictates a neoplastic transcriptional program and binds to thousands of unique sites across the genome that are enriched with PLAGL family transcription factor (TF) motifs. ZRfus activates gene expression programs through recruitment of transcriptional coactivators (Brd4, Ep300, Cbp, Pol2) that are amenable to pharmacologic inhibition. Downstream ZRfus target genes converge on developmental programs marked by PLAGL TF proteins, and activate neoplastic programs enriched in Mapk, focal adhesion, and gene imprinting networks. SIGNIFICANCE: Ependymomas are aggressive brain tumors. Although drivers of supratentorial ependymoma (ZFTA- and YAP1-associated gene fusions) have been discovered, their functions remain unclear. Our study investigates the biology of ZFTA-RELA-driven ependymoma, specifically mechanisms of transcriptional deregulation and direct downstream gene networks that may be leveraged for potential therapeutic testing.This article is highlighted in the In This Issue feature, p. 2113.

Loss of nuclear SMARCB1 (INI1/hSNF5/BAF47) protein expression due to biallelic mutations of the SMARCB1 tumor suppressor gene is a hallmark of atypical teratoid/rhabdoid tumors (ATRT), but the presence of cytoplasmic SMARCB1 protein in these tumors has not yet been described. In a series of 102 primary ATRT, distinct cytoplasmic SMARCB1 staining on immunohistochemistry was encountered in 19 cases (19%) and was highly over-represented in cases showing pathogenic sequence variants leading to truncation or mutation of the C-terminal part of SMARCB1 (15/19 vs. 4/83; Chi-square: 56.04, p = 1.0E-10) and, related to this, in tumors of the molecular subgroup ATRT-TYR (16/36 vs. 3/66; Chi-square: 24.47, p = 7.6E-7). Previous reports have indicated that while SMARCB1 lacks a bona fide nuclear localization signal, it harbors a masked nuclear export signal (NES) and that truncation of the C-terminal region results in unmasking of this NES leading to cytoplasmic localization. To determine if cytoplasmic localization found in ATRT is due to unmasking of NES, we generated GFP fusions of one of the SMARCB1 truncating mutations (p.Q318X) found in the tumors along with a p.L266A mutation, which was shown to disrupt the interaction of SMARCB1-NES with exportin-1. We found that while the GFP-SMARCB1(Q318X) mutant localized to the cytoplasm, the double mutant GFP-SMARCB1(Q318X;L266A) localized to the nucleus, confirming NES requirement for cytoplasmic localization. Furthermore, cytoplasmic SMARCB1(Q318X) was unable to cause senescence as determined by morphological observations and by senescence-associated β-galactosidase assay, while nuclear SMARCB1(Q318X;L266A) mutant regained this function. Selinexor, a selective exportin-1 inhibitor, was effective in inhibiting the nuclear export of SMARCB1(Q318X) and caused rapid cell death in rhabdoid tumor cells. In conclusion, inhibition of nuclear export restores nuclear localization and residual tumor suppressor function of truncated SMARCB1. Therapies aimed at preventing nuclear export of mutant SMARCB1 protein may represent a promising targeted therapy in ATRT harboring truncating C-terminal SMARCB1 mutations.

Pediatric high-grade glioma (pHGG) is a major contributor to cancer-related death in children. In vitro and in vivo disease models reflecting the intimate connection between developmental context and pathogenesis of pHGG are essential to advance understanding and identify therapeutic vulnerabilities. Here we report establishment of 21 patient-derived pHGG orthotopic xenograft (PDOX) models and eight matched cell lines from diverse groups of pHGG. These models recapitulate histopathology, DNA methylation signatures, mutations and gene expression patterns of the patient tumors from which they were derived, and include rare subgroups not well-represented by existing models. We deploy 16 new and existing cell lines for high-throughput screening (HTS). In vitro HTS results predict variable in vivo response to PI3K/mTOR and MEK pathway inhibitors. These unique new models and an online interactive data portal for exploration of associated detailed molecular characterization and HTS chemical sensitivity data provide a rich resource for pediatric brain tumor research.

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

The epithelial signaling pathways involved in damage and regeneration, and neoplastic transformation are known to be similar. We noted upregulation of argininosuccinate synthetase (ASS1) in hyperproliferative intestinal epithelium. Since ASS1 leads to de novo synthesis of arginine, an important amino acid for the growth of intestinal epithelial cells, its upregulation can contribute to epithelial proliferation necessary to be sustained during oncogenic transformation and regeneration. Here we investigated the function of ASS1 in the gut epithelium during tissue regeneration and tumorigenesis, using intestinal epithelial conditional Ass1 knockout mice and organoids, and tissue specimens from colorectal cancer patients. We demonstrate that ASS1 is strongly expressed in the regenerating and Apc-mutated intestinal epithelium. Furthermore, we observe an arrest in amino acid flux of the urea cycle, which leads to an accumulation of intracellular arginine. However, loss of epithelial Ass1 does not lead to a reduction in proliferation or increase in apoptosis in vivo, also in mice fed an arginine-free diet. Epithelial loss of Ass1 seems to be compensated by altered arginine metabolism in other cell types and the liver.

Ten-Eleven Translocation (TET) proteins have recently come to light as important epigenetic regulators conserved in multicellular organisms. TET knockdown studies in rodents have highlighted the critical role of these proteins for proper brain development and function. Mutations in mammalian mTET proteins and mTET2 specifically are frequent and deregulated in leukaemia and glioma respectively. Accordingly, we examined the role of mTET2 in tumorigenesis in larval haemocytes and adult heads in Drosophila melanogaster. Our findings showed that expression of mutant and wild type mTET2 resulted in general phenotypic defects in adult flies and accumulation of abdominal melanotic masses. Notably, flies with mTET2-R43G mutation at the N-terminus of mTET2 exhibited locomotor and circadian behavioural deficits, as well as reduced lifespan. Flies with mTET2-R1261C mutation in the catalytic domain, a common mutation in acute myeloid leukaemia (AML), displayed alterations affecting haemocyte haemostasis. Using transcriptomic approach, we identified upregulated immune genes in fly heads that were not exclusive to TET2 mutants but also found in wild type mTET2 flies. Furthermore, inhibiting expression of genes that were found to be deregulated in mTET2 mutants, such as those involved in immune pathways, autophagy, and transcriptional regulation, led to a rescue in fly survival, behaviour, and hemocyte number. This study identifies the transcriptomic profile of wild type mTET2 versus mTET2 mutants (catalytic versus non-catalytic) with indications of TET2 role in normal central nervous system (CNS) function and innate immunity.

Pediatric low-grade gliomas (pLGG) show heterogeneous responses to MAPK inhibitors (MAPKi) in clinical trials. Thus, more complex stratification biomarkers are needed to identify patients likely to benefit from MAPKi therapy. Here, we identify MAPK-related genes enriched in MAPKi-sensitive cell lines using the GDSC dataset and apply them to calculate class-specific MAPKi sensitivity scores (MSSs) via single-sample gene set enrichment analysis. The MSSs discriminate MAPKi-sensitive and non-sensitive cells in the GDSC dataset and significantly correlate with response to MAPKi in an independent PDX dataset. The MSSs discern gliomas with varying MAPK alterations and are higher in pLGG compared to other pediatric CNS tumors. Heterogenous MSSs within pLGGs with the same MAPK alteration identify proportions of potentially sensitive patients. The MEKi MSS predicts treatment response in a small set of pLGG patients treated with trametinib. High MSSs correlate with a higher immune cell infiltration, with high expression in the microglia compartment in single-cell RNA sequencing data, while low MSSs correlate with low immune infiltration and increased neuronal score. The MSSs represent predictive tools for the stratification of pLGG patients and should be prospectively validated in clinical trials. Our data supports a role for microglia in the response to MAPKi.

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization.

Giant cell lesions of the jaw (GCLJ) are debilitating tumors of unknown origin with limited available therapies. Here, we analyze 58 sporadic samples using next generation or targeted sequencing and report somatic, heterozygous, gain-of-function mutations in KRAS, FGFR1, and p.M713V/I-TRPV4 in 72% (42/58) of GCLJ. TRPV4 p.M713V/I mutations are exclusive to central GCLJ and occur at a critical position adjacent to the cation permeable pore of the channel. Expression of TRPV4 mutants in HEK293 cells leads to increased cell death, as well as increased constitutive and stimulated channel activity, both of which can be prevented using TRPV4 antagonists. Furthermore, these mutations induce sustained activation of ERK1/2, indicating that their effects converge with that of KRAS and FGFR1 mutations on the activation of the MAPK pathway in GCLJ. Our data extend the spectrum of TRPV4 channelopathies and provide rationale for the use of TRPV4 and RAS/MAPK antagonists at the bedside in GCLJ.

Biallelic inactivation of SMARCB1, encoding a member of the SWI/SNF chromatin remodeling complex, is the hallmark genetic aberration of atypical teratoid rhabdoid tumors (ATRT). Here, we report how loss of SMARCB1 affects the epigenome in these tumors. Using chromatin immunoprecipitation sequencing (ChIP-seq) on primary tumors for a series of active and repressive histone marks, we identified the chromatin states differentially represented in ATRTs compared with other brain tumors and non-neoplastic brain. Re-expression of SMARCB1 in ATRT cell lines enabled confirmation of our genome-wide findings for the chromatin states. Additional generation of ChIP-seq data for SWI/SNF and Polycomb group proteins and the transcriptional repressor protein REST determined differential dependencies of SWI/SNF and Polycomb complexes in regulation of diverse gene sets in ATRTs.

Ependymal tumors across age groups are currently classified and graded solely by histopathology. It is, however, commonly accepted that this classification scheme has limited clinical utility based on its lack of reproducibility in predicting patients' outcome. We aimed at establishing a uniform molecular classification using DNA methylation profiling. Nine molecular subgroups were identified in a large cohort of 500 tumors, 3 in each anatomical compartment of the CNS, spine, posterior fossa, supratentorial. Two supratentorial subgroups are characterized by prototypic fusion genes involving RELA and YAP1, respectively. Regarding clinical associations, the molecular classification proposed herein outperforms the current histopathological classification and thus might serve as a basis for the next World Health Organization classification of CNS tumors.

Neuroblastoma, a childhood cancer with highly heterogeneous biology and clinical behavior, is characterized by genomic aberrations including amplification of MYCN. Hemizygous deletion of chromosome 11q is a well-established, independent marker of poor prognosis. While 11q22-q23 is the most frequently deleted region, the neuroblastoma tumor suppressor in this region remains to be identified. Chromosome bands 11q22-q23 contain ATM, a cell cycle checkpoint kinase and tumor suppressor playing a pivotal role in the DNA damage response. Here, we report that haploinsufficiency of ATM in neuroblastoma correlates with lower ATM expression, event-free survival, and overall survival. ATM loss occurs in high stage neuroblastoma without MYCN amplification. In SK-N-SH, CLB-Ga and GI-ME-N human neuroblastoma cells, stable ATM silencing promotes neuroblastoma progression in soft agar assays, and in subcutaneous xenografts in nude mice. This effect is dependent on the extent of ATM silencing and does not appear to involve MYCN. Our findings identify ATM as a potential haploinsufficient neuroblastoma tumor suppressor, whose inactivation mirrors the increased aggressiveness associated with 11q deletion in neuroblastoma.

Neuroblastoma (NB) is an aggressive childhood malignancy in children up to 5 years of age. High-stage tumors frequently relapse even after aggressive multimodal treatment, and then show therapy resistance, typically resulting in patient death. New molecular-targeted compounds that effectively suppress tumor growth and prevent relapse with more efficacy are urgently needed. We and others previously showed that polyamines (PA) like spermidine and spermine are essential for NB tumorigenesis and that DFMO, an inhibitor of the key PA synthesis gene product ODC, is effective both in vitro and in vivo, securing its evaluation in NB clinical trials. To find additional compounds interfering with PA biosynthesis, we tested sulfasalazine (SSZ), an FDA-approved salicylate-based anti-inflammatory and immune-modulatory drug, recently identified to inhibit sepiapterin reductase (SPR). We earlier presented evidence for a physical interaction between ODC and SPR and we showed that RNAi-mediated knockdown of SPR expression significantly reduced native ODC enzyme activity and impeded NB cell proliferation.

We recently reported that atypical teratoid rhabdoid tumors (ATRTs) comprise at least two transcriptional subtypes with different clinical outcomes; however, the mechanisms underlying therapeutic heterogeneity remained unclear. In this study, we analyzed 191 primary ATRTs and 10 ATRT cell lines to define the genomic and epigenomic landscape of ATRTs and identify subgroup-specific therapeutic targets. We found ATRTs segregated into three epigenetic subgroups with distinct genomic profiles, SMARCB1 genotypes, and chromatin landscape that correlated with differential cellular responses to a panel of signaling and epigenetic inhibitors. Significantly, we discovered that differential methylation of a PDGFRB-associated enhancer confers specific sensitivity of group 2 ATRT cells to dasatinib and nilotinib, and suggest that these are promising therapies for this highly lethal ATRT subtype.

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs, including 19 angiocentric gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in angiocentric gliomas. In vitro and in vivo functional studies show that MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression and hemizygous loss of the tumor suppressor QKI. To our knowledge, this represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

Retinoic acid (RA) induces differentiation of neuroblastoma cells in vitro and is used with variable success to treat aggressive forms of this disease. This variability in clinical response to RA is enigmatic, as no mutations in components of the RA signaling cascade have been found. Using a large-scale RNAi genetic screen, we identify crosstalk between the tumor suppressor NF1 and retinoic acid-induced differentiation in neuroblastoma. Loss of NF1 activates RAS-MEK signaling, which in turn represses ZNF423, a critical transcriptional coactivator of the retinoic acid receptors. Neuroblastomas with low levels of both NF1 and ZNF423 have extremely poor outcome. We find NF1 mutations in neuroblastoma cell lines and in primary tumors. Inhibition of MEK signaling downstream of NF1 restores responsiveness to RA, suggesting a therapeutic strategy to overcome RA resistance in NF1-deficient neuroblastomas.