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

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames (ORFs). To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a stepwise approach using multiple CRISPR-Cas9 screens to elucidate non-canonical ORFs and putative microproteins implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream ORFs (uORFs) exhibited selective functionality independent of main coding sequences. A microprotein encoded by one of these ORFs, ASNSD1-uORF or ASDURF, was upregulated, associated with MYC-family oncogenes, and promoted medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future studies seeking to define new cancer targets.

Abstract Background: Pediatric brain tumors are associated with the highest rates of morbidity and mortality of all childhood cancers. Diffuse midline gliomas (DMGs) are universally fatal pediatric tumors with a median survival of less than one year. DMGs are thought to arise from stalled developmental programs in the context of histone mutations (H3 K27M) which are thought to be a critical initiating event. However, these histone mutations alone are insufficient for gliomagenesis and cooperate with other alterations, including mutations in the TP53 pathway (TP53, PPM1D), growth factor receptor pathway activation (PDGFRA, EGFR, PIK3CA), and MYC amplifications (MYC, MYC-N) to induce tumor formation. There is a high degree of genetic heterogeneity amongst DMGs, and other alterations are also critical to DMG pathogenesis. This work is focused on Forkhead Box R2 (FOXR2), a forkhead family TF that we have found to be aberrantly expressed in many cancers, including DMGs. Expression of FOXR2 is sufficient to enhance glioma formation, and we have shown that regulation of FOXR2 occurs through a previously unrecognized epigenetic mechanism. Expression of FOXR2 is regulated by novel promoters, and we have demonstrated that these promoters are necessary for FOXR2 expression in several cancer lineages using CRISPR interference technology, including DMGs, melanoma, neuroblastoma, and non-small cell lung cancers. Methods: We applied an integrative approach using transcriptomics, epigenetics, proteomics, in vitro cancer models, and in vivo mouse models to systematically evaluate how FOXR2 mediates gliomagenesis. Results: We have found that FOXR2 is activated across multiple cancer lineages, including DMGs. FOXR2 is required for proliferation in vitro, and it is sufficient to enhance gliomagenesis using an in utero electroporation in vivo mouse model. FOXR2 enhances MYC protein stability, even in the presence of cycloheximide. However, FOXR2 also exerts its oncogenic effects through MYC independent functions. FOXR2 is highly enriched at E26-transformation (ETS) motifs and specifically activates ETS transcriptional gene sets. To determine protein-protein interactions of FOXR2, we have performed immunoprecipitation and mass spectrometry to identify FOXR2 protein interactors. Moreover, we have performed both proteomic and phospho-proteomic analyses of FOXR2-expressing human neural stem cells. These proteomic studies have allowed us to identify proteins and phospho-sites that are highly enriched in FOXR2-expressing cells, and that could be potential therapeutic targets in FOXR2-expressing cancers. Conclusion: Taken together, this study elucidates how FOXR2 protein interactors mediate oncogenesis in FOXR2-expressing diffuse midline gliomas. These studies proposed here have broad applicability across a range of cancers as FOXR2 is aberrantly expressed in 8% of all cancers. Citation Format: Jessica W. Tsai, Marissa Coppola, Phonepasong Arounleut, Jared B. Collins, Patrick Hart, Hasmik Keshishian, David T. Jones, Pratiti Bandopadhayay, Timothy N. Phoenix. Delineating mediators of oncogenesis in FOXR2-expressing cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2854.

Abstract Since the KIAA1549::BRAF fusion was discovered as the most common driver of pediatric low-grade glioma (pLGG), it has been hypothesized that the fusion induces oncogenicity through BRAF activation as a result of KIAA1549 replacing the BRAF negative regulatory N-terminus. This led to the rapid translation of MAPK pathway inhibitors into the clinical setting. Despite most tumors exhibiting promising initial responses, some tumors are not sensitive and about half of responsive tumors grow back after treatment cessation. Therefore, strategies that result in sustained tumor responses are desperately needed. We have recently performed genome-scale CRISPR/Cas9 screens across isogenic neural stem cell models transduced to express pLGG-associated oncogenes (KIAA1549::BRAF, BRAFV600E, and multiple FGFR1 and MYB family alterations) to generate a dependency map of genetic vulnerabilities associated with expression of these oncogenes. We also included normal neural stem cells to allow identification of genetic dependencies specifically induced by the expression of each oncogene. Surprisingly, through these efforts, we have discovered KIAA1549::BRAF expressing cells to harbor striking and specific dependency on multiple members of an enzymatic complex that exerts its activity outside of the MAPK signaling axis*. Interestingly, this enzymatic complex has been described to modify only a few substrates, including KIAA1549, suggesting specificity. We have now validated this dependency across other isogenic models of KIAA1549::BRAF-expressing cells. Finally, our KIAA1549::BRAF-expressing models exhibit preferential sensitivity to a tool compound that targets our novel enzyme, suggesting novel therapeutic potential for KIAA1549::BRAF. These findings highlight a MAPK pathway-independent avenue for therapeutically targeting the most frequent genomic alterations in pediatric brain tumors. Consequently, we also propose that the fusion partner KIAA1549 is instrumental for the aberrant BRAF signaling driving pLGGs. (* we are currently in the process of working with our IP offices to ensure that we can disclose the names of the genes and proteins at the ISPNO meeting)

Abstract t was originally hypothesized that loss of BRAF N-terminal negative regulatory domains, which are deleted in the KIAA1549:BRAF rearrangement, results in constitutive activation of the BRAF kinase. This model suggests that BRAF fusion partners are dispensable for BRAF oncogenic signaling. Paradoxically, our data suggest expression of truncated BRAF is insufficient for transformation and identify specific domains in KIAA1549 that are necessary for transformation. These domains are critical for regulating subcellular localization of the fusion protein and highlight an aberrant pattern of cellular localization that is not observed with wildtype BRAF. Altered subcellular localization results in proteolytic cleavage of KIAA1549:BRAF, and we identify the protease responsible for this cleavage event. This protease is therapeutically tractable and clinically relevant inhibitors have been developed, presenting a MAPK-independent mechanism to target the fusion. We also highlight an unexpected role for rare fusion partners, including FAM131B, in BRAF activation. In total, these data suggest that BRAF fusion partners are not indispensable for transformation as was previously thought, presenting unexplored opportunities to therapeutically target pLGG tumors with BRAF rearrangements. We are currently working on the IP measures required to disclose the specific enzyme(s) at ISPNO.

Abstract High-risk medulloblastoma is one of the most recalcitrant pediatric cancers, and children with MYC-amplified disease frequently succumb to relapsed disease. Extensive analyses of the coding genome in this disease have characterized additional somatic events in some subsets of patients, though most tumors lack targetable mutations and do not yield insights regarding their aggressive behavior. At the same time, medulloblastoma is known to exhibit extensive rewiring of translational control in MYC-driven tumors, consistent with recent genetic evidence that the impact of this transcription factor on control of mRNA translation may be the most critical aspect of its function during tumorigenesis. Therefore, to propose previously unknown mechanisms for this disease, we have investigated the functional impact of translation of non-canonical open reading frames (ORFs) across medulloblastoma model systems. We demonstrate that these ORFs are commonly translated in medulloblastoma model systems and patient tumors, correlating with disease subtype. Using genome-wide CRISPR/Cas9 screens, we found that ORFs are frequently essential for cell survival in medulloblastoma and describe widespread reliance on upstream open reading frames (uORFs) in particular. From these, we identify a uORF in the ASNSD1 gene that is selectively upregulated and required for maintenance of cell survival by coordinating the function of the prefoldin-like complex, a poorly understood complex implicated in post-translational control. Together, our findings provide a blueprint for oncogenic uORFs as critical disease mediators both in medulloblastoma and, by extension, human cancers more broadly.

Supplementary Table from FOXR2 Is an Epigenetically Regulated Pan-Cancer Oncogene That Activates ETS Transcriptional Circuits

Pediatric low-grade gliomas (pLGGs) are the most common solid tumors in children and are associated with devastating lifelong morbidities and mortality. Recent genomic profiling efforts have revealed that these tumors are largely driven by single-driver events that activate MAPK signaling. These insights have led to early clinical trials evaluating the role of MAPK inhibitors for these children with promising initial results. However, pLGGs are not cured by MAPK inhibitors, and tumors often rapidly rebound upon cessation of treatment. Therefore, continuous dosing of these inhibitors is required throughout a child’s development–often with significant and potentially permanent toxicities. In addition, 30-40% of tumors do not respond to MAPK inhibition, indicating primary resistance. We hypothesized that pediatrlc low-grade gliomas may also harbor additional dependencies beyond the MAPK pathway that may represent potential therapeutic targets. To address this, we generated isogenic neural stem cell models that expressed pLGG relevant oncogenes (BRAF family members, FGFR family members and MYB transcription factors, and subjected them to genome-scale CRISPR-cas9 screens. Compared to vector controls, pLGG models exhibited genetic dependencies on pathways that included regulators of cell-cycle and differentiation, in addition to genes involved in regulating the mTOR and MAPK pathways. Cell lines expressing MYB family transcription factors also harbored dependencies on genes involved in the DREAM complex. Together, these genes and pathways represent potential targets for combination treatments with MAPK pathway inhibitors for the treatment of pediatric low-grade gliomas.

Neuroblastoma (NB) is a highly aggressive extracranial solid tumor in children. Due to its heterogeneity, NB remains a therapeutic challenge. Several oncogenic factors, including the Hippo effectors YAP/TAZ, are associated with NB tumorigenesis. Verteporfin (VPF) is an FDA-approved drug shown to directly inhibit YAP/TAZ activity. Our study aimed to investigate VPF’s potential as a therapeutic agent in NB. We show that VPF selectively and efficiently impairs the viability of YAP/TAZ-expressing NB GI-ME-N and SK-N-AS cells, but not of non-malignant fibroblasts. To investigate whether VPF-mediated NB cell killing is YAP-dependent, we tested VPF potency in CRISPR-mediated YAP/TAZ knock-out GI-ME-N cells, and BE(2)-M17 NB cells (a MYCN-amplified, predominantly YAP-negative NB subtype). Our data shows that VPF-mediated NB cell killing is not dependent on YAP expression. Moreover, we determined that the formation of higher molecular weight (HMW) complexes is an early and shared VPF-induced cytotoxic mechanism in both YAP-positive and YAP-negative NB models. The accumulation of HMW complexes, involving STAT3, GM130 and COX IV proteins, impaired cell homeostasis and triggered cell stress and cell death mechanisms. Altogether, our study shows significant in vitro and in vivo VPF-induced suppression of NB growth, making VPF a potential therapeutic candidate against NB.

Abstract BACKGROUND Diffuse midline gliomas (DMGs) are a universally fatal brain tumor of childhood. While histone mutations are a critical tumor initiating event, they are insufficient to drive gliomagenesis. Histone mutations co-occur with somatic alterations in other pathways including TP53, MAPK, and MYC signaling. However, the mechanisms through which these pathways are activated have not been fully elucidated. METHODS We applied an integrative approach using transcriptomics, epigenetics, proteomics, in vitro cancer models, and in vivo mouse models to systematically evaluate how FOXR2 mediates gliomagenesis. RESULTS We have recently found that a subset of DMGs aberrantly express FOXR2, a forkhead transcription factor. FOXR2 is both sufficient to enhance tumor formation, and necessary for FOXR2-expressing DMGs. While FOXR2 indeed enhances MYC protein stability, FOXR2 exerts oncogenesis through MYC-independent functions and specifically hijacks E26-transformation specific (ETS) transcriptional circuits and FOXR2 DNA-binding is highly enriched at ETS motifs. We have performed proteomic and phospho-proteomic analysis of FOXR2-expressing human neural stem cells to identify proteins and phospho-sites that are highly enriched in FOXR2-expressing cells. CONCLUSION Taken together, this study elucidates how FOXR2 interacts with ETS transcription factors to mediate oncogenesis in FOXR2-expressing diffuse midline gliomas.

Spheroid culture systems have allowed in vitro propagation of cells unable to grow in canonical cell culturing conditions, and may capture cellular contexts that model tumor growth better than current model systems. The insights gleaned from genome-wide clustered regularly interspaced short palindromic repeat (CRISPR) screening of thousands of cancer cell lines grown in conventional culture conditions illustrate the value of such CRISPR pooled screens. It is clear that similar genome-wide CRISPR screens of three-dimensional spheroid cultures will be important for future biological discovery. Here, we present a protocol for genome-wide CRISPR screening of three-dimensional neurospheres. While many in-depth protocols and discussions have been published for more typical cell lines, few detailed protocols are currently available in the literature for genome-wide screening in spheroidal cell lines. For those who want to screen such cell lines, and particularly neurospheres, we provide a step-by-step description of assay development tests to be performed before screening, as well as for the screen itself. We highlight considerations of variables that make these screens distinct from, or similar to, typical nonspheroid cell lines throughout. Finally, we illustrate typical outcomes of neurosphere genome-wide screens, and how neurosphere screens typically produce slightly more heterogeneous signal distributions than more canonical cancer cell lines. Completion of this entire protocol will take 8–12 weeks from the initial assay development tests to deconvolution of the sequencing data.

Background: Diffuse midline gliomas (DMGs) are a universally fatal brain tumor of childhood. While histone mutations are a critical tumor initiating event, they are insufficient to drive gliomagenesis. Histone mutations co-occur with somatic alterations in other pathways including TP53, MAPK, and MYC signaling. However, the mechanisms through which these pathways are activated have not been fully elucidated. Methods: We applied an integrative approach using transcriptomics, epigenetics, proteomics, in vitro cancer models, and in vivo mouse models to systematically evaluate how FOXR2 mediates gliomagenesis. Results: We have recently found that a subset of DMGs aberrantly express FOXR2, a forkhead transcription factor. FOXR2 is both sufficient to enhance tumor formation, and necessary for FOXR2-expressing DMGs. While FOXR2 indeed enhances MYC protein stability, FOXR2 exerts oncogenesis through MYC-independent functions and specifically hijacks E26-transformation specific (ETS) transcriptional circuits and FOXR2 DNA-binding is highly enriched at ETS motifs. We have performed proteomic and phospho-proteomic analysis of FOXR2-expressing human neural stem cells to identify proteins and phospho-sites that are highly enriched in FOXR2-expressing cells. Conclusion: Taken together, this study elucidates how FOXR2 interacts with ETS transcription factors to mediate oncogenesis, and further highlights a role for FOXR2 in activating ETS and MAPK signaling. Citation Format: Jessica W. Tsai, Paloma Cejas, Marissa Coppola, Dayle K. Wang, Smruti Patel, David W. Wu, Phonepasong Arounleut, Xin Wei, Ningxuan Zhou, Sudeepa Syamala, Frank P. Dubois, Kristine Pelton, Jayne Vogelzang, Cecilia Sousa, Audrey Baguette, Xiaolong Chen, Alexandra L. Condurat, Sarah E. Dixon-Clarke, Annarah Charles, Kevin N. Zhou, Sophie D. Lu, Elizabeth M. Gonzalez, Madison S. Chacon, Jeromy J. Digiacomo, Rushil Kumbhani, Dana Novikov, Maria Tsoli, David S. Ziegler, Uta Dirksen, Natalie Jager, Gnana Prakash Balasubramanian, Christof M. Kramm, Michaela Nathrath, Stefan Bielack, Suzanne J. Baker, Jinghui Zhang, James M. McFarland, Gad Getz, Francois Aguet, Nada Jabado, Olaf Witt, Stefan M. Pfister, Keith L. Ligon, Volker Hovestadt, Claudia Kleinman, Henry Long, David T. Jones, Pratiti Bandopadhayay, Timothy N. Phoenix. Dissecting mechanisms underlying FOXR2-mediated gliomagenesis in diffuse midline gliomas. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3562.

BACKGROUND: Genomic and transcriptomic studies have elucidated new insights into the landscape of diffuse intrinsic pontine glioma (DIPG). However, the role of long non-coding RNAs (lncRNAs) has not been explored at depth in these tumors, and there have not been studies focused on how lncRNAs interact with the K27M histone mutation. In a recent analysis of nearly 200 DIPGs and pediatric high-grade gliomas (pHGG), we previously detected a novel, recurring structural variant in the lncRNA CCDC26. This rearrangement occurs in nearly 10% of all DIPGs, and we have furthermore identified alterations in more than 100 lncRNAs in DIPG. METHODS: To identify lncRNAs required for proliferation of patient-derived DIPG cancer cells, we designed two custom genome-scale lncRNA libraries. We generated a genome-scale lncRNA CRISPR-Cas9 knockout pooled library, consisting of 45,766 single guide RNAs (sgRNAs). Additionally, we generated a genome-wide CRISPR interference pooled library consisting of 45,608 sgRNAs, targeting lncRNA transcription start sites (TSS). RESULTS: We utilized in vitro histone-mutant pHGG models as well as edited clones of these models with the K27M mutant corrected in order to compare lncRNA dependencies in these two contexts. We have successfully performed genome-scale CRISPR-Cas9 knockout and CRISPR interference screens targeting lncRNAs in these cell lines, revealing lncRNA dependencies. Candidate dependencies in our CRISPR-Cas9 knockout screen include LOC100507412, LOC105379524, and LINC02193. CONCLUSION: Genome-wide lncRNA CRISPR knock-out and CRISPR interference screens are a novel approach for the unbiased identification of lncRNAs that are required for pediatric high-grade glioma proliferation. Further validation of specific lncRNAs is required, and these lncRNA dependencies represent potential novel therapeutic targets.

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

The role of PPM1D mutations in de novo gliomagenesis has not been systematically explored. Here we analyze whole genome sequences of 170 pediatric high-grade gliomas and find that truncating mutations in PPM1D that increase the stability of its phosphatase are clonal driver events in 11% of Diffuse Midline Gliomas (DMGs) and are enriched in primary pontine tumors. Through the development of DMG mouse models, we show that PPM1D mutations potentiate gliomagenesis and that PPM1D phosphatase activity is required for in vivo oncogenesis. Finally, we apply integrative phosphoproteomic and functional genomics assays and find that oncogenic effects of PPM1D truncation converge on regulators of cell cycle, DNA damage response, and p53 pathways, revealing therapeutic vulnerabilities including MDM2 inhibition.

Abstract AIM: Pediatric low-grade gliomas (pLGGs) are a heterogenous group of tumors, diverse in their localization, histology, mutational landscape, clinical behavior, and treatment response. Genomic alterations impacting the MYB family of transcription factors were identified in two distinct pLGG subtypes: Angiocentric Gliomas (AG) and Diffuse Astrocytomas (DA). The molecular profiles and therapeutic vulnerabilities associated with these genomic alterations remain unexplored. In this study we highlight the use of genome-wide CRISPR/Cas9 knock-out screens for an unbiased identification of translatable therapeutic targets. METHODOLOGY: Given the lack of patient-derived pLGG cell lines, we engineered in vitro pLGG mouse and human neural stem cell (NSC) models to harbor pLGG-relevant genomic alterations. We performed single cell RNA sequencing to investigate the transcriptional profiles driven by these mutations and to dissect the central regulatory networks enabling tumorigenesis. Specific genetic dependencies associated with MYB/MYBL1 mutations were screened using the Brie genome-wide mouse CRISPR lentiviral knock-out pooled library, consisting of 78,637 single guide RNAs (sgRNAs) targeting 19,674 mouse genes. RESULTS: We have successfully generated in vitro NSC-based pLGG models crucial to deepening our knowledge on pLGG biology and the identification of translatable therapeutic targets. Genome-scale CRISPR/Cas9 knock-out screens in isogenic NSCs models, expressing distinct MYB/MYBL1 alterations or a control transgene, revealed several differential genetic dependencies. Among the top identified dependencies are regulators of cell-stress response, cell-cycle progression, and modulators of the ubiquitin-proteasome degradation pathway. CONCLUSION: Genome-wide CRISPR knock-out screens are a powerful tool for the unbiased identification of mutation-specific genetic dependencies that can be explored as candidates for precision medicine approaches.

Background: Understanding how aberrant transcription factors (TFs) hijack normal development to induce oncogenesis is a critical question in oncology. Forkhead box (FOX) proteins are a superfamily of transcriptional regulators characterized by a forkhead DNA-binding domain. Within this family, Forkhead Box R2 (FOXR2) expression has been associated with a subset of cancers including CNS and peripheral neuroblastoma. While FOXR2 has been shown to stabilize MYC isoforms, the mechanistic details through which it enhances tumor formation and the true extent of its role as an oncogene across all cancers have not been systematically evaluated. Methods: We applied an integrative approach using transcriptomics, epigenetics, in vitro cancer models, and in vivo mouse models to systematically evaluate the mechanisms by which FOXR2 is activated across human cancers. Results: We performed a pan-cancer analysis of FOXR2 activation across over 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. We identified genetic and epigenetic mechanisms that induce its expression, including hypomethylation of a novel promoter in the vast majority (78%) of FOXR2-expressing cases. We demonstrate that FOXR2 expression is both sufficient and necessary for transformation across multiple lineages, using both in vitro and in vivo models. Conclusion: Taken together, this study demonstrates the role of FOXR2 as a potent oncogene across human cancers, and highlights a novel mechanism by which its expression is activated. Citation Format: Jessica W. Tsai, Paloma Cejas, Dayle K. Wang, Smruti Patel, David W. Wu, Phonepasong Arounleut, Xin Wei, Ningxuan Zhou, Sudeepa Syamala, Frank P. Dubois, Kristine Pelton, Jayne Vogelzang, Cecilia Sousa, Audrey Baguette, Xiaolong Chen, Alexandra L. Condurat, Sarah E. Dixon-Clarke, Kevin N. Zhou, Sophie D. Lu, Elizabeth M. Gonzalez, Madison S. Chacon, Jeromy J. Digiacomo, Rushil Kumbhani, Dana Novikov, J'Ya Hunter, Maria Tsoli, David S. Ziegler, Uta Dirksen, Natalie Jager, Gnana Prakash Balasubramanian, Christof M. Kramm, Michaela Nathrath, Stefan Bielack, Suzanne J. Baker, Jinghui Zhang, James M. McFarland, Gad Getz, Francois Aguet, Nada Jabado, Olaf Witt, Stefan M. Pfister, Keith L. Ligon, Claudia Kleinman, Henry Long, David T. Jones, Pratiti Bandopadhayay, Timothy N. Phoenix. FOXR2 is an oncogenic driver across adult and pediatric cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5730.

Abstract 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. Significance: This work identifies a novel promoter that drives aberrant FOXR2 expression and delineates FOXR2 as a pan-cancer oncogene that specifically activates ETS transcriptional circuits across human cancers. See related commentary by Liu and Northcott, p. 2977

Abstract Genomic drivers of pediatric low-grade gliomas (pLGGs) converge on alterations that activate the MAPK pathway. However, expression of individual driver oncogenes fails to induce tumor formation with high penetrance and, paradoxically, expression of these oncogenes suppresses growth in vitro. This, combined with the non-monotonic tumor growth rate in patients, suggests that there are “hidden drivers” beyond a single driver oncogene that are necessary to support tumor growth. The goal of this project is to leverage high-throughput functional genomics strategies to identify these hidden drivers of pLGG. Our preliminary data indicates that genes which modulate differentiation are required for the survival of LGG cells, suggesting that these genes may be hidden drivers of LGG tumor growth. Additionally, we hypothesize that secreted factors in the tumor microenvironment regulate pLGG tumor growth, potentially by modulating differentiation. In total, genes which cooperate with pLGG driver oncogenes to promote tumor growth may represent a new class of therapeutic targets and may explain the complex patterns of tumor growth that are observed in patients.