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.

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 Adamantinomatous Craniopharyngioma (ACP) is a rare pituitary tumor of the sellar region, with a bimodal presentation pattern in children and middle-aged adults. ACPs are primarily driven by a mutation in CTNNB1 that prevents beta-catenin from degrading properly, resulting in accumulation in the nuclei of epithelial whorl-like structures sometimes referred to as cluster cells. Clinical management is challenging due to the location of the tumor, such that complete surgical resection is often not possible so radiation therapy, which can lead to high morbidity, is also used. Greater insight into these tumors can help develop less damaging treatment strategies. In this analysis we have employed single-cell RNA sequencing technology to examine seven pediatric human ACPs (28594 cells). Using consensus non-negative matrix factorization, and traditional marker genes, we have partitioned our dataset into several cell types including cluster cells, palisading epithelial cells, microglia, T-cells, B-cells, plasma cells, and fibroblasts. Though cluster cells themselves do not actively proliferate, we believe the complex inter-cellular signaling is responsible for tumor formation. We propose a model in which malignant fibroblasts form from the existing epithelium and remodel the extracellular matrix (ECM) creating a stiffened tumor microenvironment. To support this claim we show activation markers suggesting myofibroblast differentiation, the expression of EMT signatures indicating a transformation from epithelial to fibroblast, and RNA velocity analysis that demonstrates lineage tracing. We believe that the remodeling of the ECM is a wound healing response influenced by TGF-beta, SPP1, and other signaling pathways. Our analysis of ACPs reveals the role of the tumor microenvironment in the formation of these tumors.

Abstract BACKGROUND Pediatric brain cancers are the leading cause of cancer-related deaths in children. We (and others) have recently found that close to 10% of all pediatric gliomas, encompassing low-grade gliomas (LGGs) and high-grade gliomas (HGGs), harbor recurrent driver alterations in FGFR proteins, most frequently FGFR1, the gene encoding Fibroblast Growth Factor Receptor 1. In pediatric gliomas, FGFR1 alterations present as either structural variants (SVs), including kinase duplications and fusion proteins, or kinase-activating single nucleotide variants (SNVs). Recurrent FGFR1 alterations represent a promising therapeutic target for precision medicine approaches, however a challenge has been the lack of clinical trials for FGFR inhibitors. A major goal of this project is to identify preclinical FGFR inhibitor candidates for FGFR1-altered pediatric gliomas. METHODS We have generated isogenic mouse and human neural stem cell (NSC) models driven by FGFR1 and BRAF alterations commonly seen in pediatric gliomas to test the oncogenic properties and therapeutic capacities of these alterations. RESULTS We found that NSC lines grow independent of growth factors and form tumors in mice. RNA-sequencing of these lines revealed the FGFR1-altered lines are enriched in neuronal gene programs. The FGFR1-altered models are exquisitely sensitive to the four FDA-approved panFGFR inhibitors Infigratinib, Erdafitinib, Pemigatinib, and Futibatinib. However, the FGFR1 SNV + PTPN11 co-occurring SNV line is the least sensitive to FGFR inhibition, which suggests combination therapies may be needed. Combination drug studies with FGFR inhibition in our FGFR1-altered lines reveal the best synergy with mTOR inhibition. Finally, we present examples of case studies of patients with FGFR1-altered tumors treated with targeted inhibitors. CONCLUSIONS Taken together, our studies have established models that have generated key insights into the biology of FGFR1-altered gliomas and how to best target them therapeutically.

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 p53 pathway reactivation through pharmacologic inhibition of negative regulators was identified as a promising therapeutic strategy for TP53 wild-type (WT) high-grade gliomas, including Diffuse Midline Gliomas (DMGs). Despite widespread theoretical promise, the therapeutic efficacy of MDM2 inhibitors and PPM1D inhibitors as single agents have each been limited by issues, including sub-lethal effects on cancer cells and innate or emerging resistance leading to tumor progression. We carried out a genome-scale CRISPR activating (CRISPRa) screen to identify genes that, when overexpressed, confer resistance to p53 reactivation therapies. NAD(P)H quinone dehydrogenase 1 (NQO1) emerged as a top hit and was confirmed through secondary CRISPRa screens and low-throughput validations as driving a resistance phenotype in TP53 WT DMGs. Our subsequent work has aimed to elucidate the mechanisms through which NQO1 mediates the effects of p53 reactivation therapies, and assess the feasibility and impact of exploiting NQO1 in combination therapies towards a cure.

Alterations in Fibroblast growth factor receptor (FGFR)- family proteins frequently occur as oncogenes in many cancers, including a subset of pediatric gliomas. Here, we performed a genomic analysis of 11,635 gliomas across ages and found that 4.5% of all gliomas harbor FGFR alterations including structural variants (SV) and single nucleotide variants (SNV), with an incidence of almost 10% in pediatric gliomas. FGFR family members are differentially enriched by age, tumor grade, and histological subtype, with FGFR1-alterations associated with glioneuronal histologies and pediatric low-grade gliomas. Across development, we find FGFR1 expression in both neuronal and glial precursors, while FGFR3 expression is largely restricted to astrocytic lineages. Leveraging novel isogenic model systems, we confirm FGFR1 alterations to be sufficient to activate MAPK and mTOR signaling, drive gliomagenesis, activate neuronal transcriptional programs and exhibit sensitivity to MAPK pathway inhibitors, including pan-FGFR inhibitors. Models driven by FGFR1 SVs exhibited different patterns of sensitivity compared to those driven by SNVs. Finally, we performed a retrospective analysis of clinical responses in children diagnosed with FGFR-driven gliomas and found that targeted MAPK or FGFR-inhibition with currently available inhibitors is largely associated with stability of disease. This study provides key insights into the biology of FGFR1-altered gliomas, therapeutic strategies to target them and associated challenges that still need to be overcome. ### Competing Interest Statement PB serves on paid advisory boards for DayOne Biopharmaceuticals and has served on a paid advisory board for QED Therapeutics. Her lab has received grant funding from Novartis Institute of Biomedical Research. SHR has employment at Labcorp Oncology. CMVT is on advisory boards for Alexion, Bayer, and Novartis. OW is advisory board member or has received research grants from Janssen, Day One Biopharmaceuticals, and Novartis. K.L.L. disclosures: Equity: Travera Inc.; Consulting- Travera Inc., BMS, Servier, LEK, Integragen, Blaze Bioscience; Research.

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.

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.

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.

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.

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.

Abstract Pediatric low-grade gliomas (pLGGs) have excellent survival, however, with current standard of care, most patients suffer lifelong severe sequalae. pLGGs are almost exclusively driven by single activating mutations in the MAPK pathway. Clinical trials with small molecule inhibitors in BRAF-altered pLGGs are showing promising results in early clinical trials, and similar efforts are now underway for FGFR1-altered tumors, however the underlying biology and treatment response has not been thoroughly explored in a pre-clinical setting. To explore the genetic landscape of FGFR altered gliomas we assembled a cohort of 87 patients with FGFR1-4 altered gliomas across Dana-Farber Cancer Institute, Boston Children’s Hospital and Brigham and Women’s Hospital. Within this cohort we observed that pLGGs harboring FGFR1 kinase hotspot mutations (FGFR1-N546K or -K656E) frequently harbored a second alteration associated with activation of the MAPK or mTOR pathways, most commonly in the phosphatase PTPN11, NF1 or within the FGFR1 gene itself. Additionally, we observed two previously described structural variants of FGFR1, an FGFR1 internal kinase tandem duplication (FGFR-ITD) and a fusion with TACC1 (FGFR1:TACC1). The relative impact of the different FGFR1 alterations on oncogenicity, therapeutic response and resistance has not been previously explored. To address this, we have established mouse neural stem cell models overexpressing the structural variants and hot spot mutant FGFR1 alone or in combination with a second alteration. Immunoblotting revealed that the addition of a second alteration attenuated phosphorylation of ERK, AKT and S6 and influenced cell proliferation both in normal growth conditions and in absence of growth factor. Treatment with inhibitors of FGFR (Infigratinib) and MEK (Trametinib) revealed variable sensitivity both targeted therapies, suggesting that treatment of FGFR1 driven pLGG might require tailoring to the specific FGFR1 alteration.

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.

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