KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma.

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

Piperacetazine Directly Binds to the PAX3::FOXO1 Fusion Protein and Inhibits Its Transcriptional Activity.

The tumor-specific chromosomal translocation product, PAX3::FOXO1, is an aberrant fusion protein that plays a key role for oncogenesis in the alveolar subtype of rhabdomyosarcoma (RMS). PAX3::FOXO1 represents a validated molecular target for alveolar RMS and successful inhibition of its oncogenic activity is likely to have significant clinical applications. Even though several PAX3::FOXO1 function-based screening studies have been successfully completed, a directly binding small-molecule inhibitor of PAX3::FOXO1 has not been reported. Therefore, we screened small-molecule libraries to identify compounds that were capable of directly binding to PAX3::FOXO1 protein using surface plasmon resonance technology. Compounds that directly bound to PAX3::FOXO1 were further evaluated in secondary transcriptional activation assays. We discovered that piperacetazine can directly bind to PAX3::FOXO1 protein and inhibit fusion protein-derived transcription in multiple alveolar RMS cell lines. Piperacetazine inhibited anchorage-independent growth of fusion-positive alveolar RMS cells but not embryonal RMS cells. On the basis of our findings, piperacetazine is a molecular scaffold upon which derivatives could be developed as specific inhibitors of PAX3::FOXO1. These novel inhibitors could potentially be evaluated in future clinical trials for recurrent or metastatic alveolar RMS as novel targeted therapy options. SIGNIFICANCE: RMS is a malignant soft-tissue tumor mainly affecting the pediatric population. A subgroup of RMS with worse prognosis harbors a unique chromosomal translocation creating an oncogenic fusion protein, PAX3::FOXO1. We identified piperacetazine as a direct inhibitor of PAX3::FOXO1, which may provide a scaffold for designing RMS-specific targeted therapy.

Preclinical development of a chimeric antigen receptor T cell therapy targeting FGFR4 in rhabdomyosarcoma.

Pediatric patients with relapsed or refractory rhabdomyosarcoma (RMS) have dismal cure rates, and effective therapy is urgently needed. The oncogenic receptor tyrosine kinase fibroblast growth factor receptor 4 (FGFR4) is highly expressed in RMS and lowly expressed in healthy tissues. Here, we describe a second-generation FGFR4-targeting chimeric antigen receptor (CAR), based on an anti-human FGFR4-specific murine monoclonal antibody 3A11, as an adoptive T cell treatment for RMS. The 3A11 CAR T cells induced robust cytokine production and cytotoxicity against RMS cell lines in vitro. In contrast, a panel of healthy human primary cells failed to activate 3A11 CAR T cells, confirming the selectivity of 3A11 CAR T cells against tumors with high FGFR4 expression. Finally, we demonstrate that 3A11 CAR T cells are persistent in vivo and can effectively eliminate RMS tumors in two metastatic and two orthotopic models. Therefore, our study credentials CAR T cell therapy targeting FGFR4 to treat patients with RMS.

A stem cell epigenome is associated with primary nonresponse to CD19 CAR T cells in pediatric acute lymphoblastic leukemia.

CD19 chimeric antigen receptor T-cell therapy (CD19-CAR) has changed the treatment landscape and outcomes for patients with pre-B-cell acute lymphoblastic leukemia (B-ALL). Unfortunately, primary nonresponse (PNR), sustained CD19+ disease, and concurrent expansion of CD19-CAR occur in 20% of the patients and is associated with adverse outcomes. Although some failures may be attributable to CD19 loss, mechanisms of CD19-independent, leukemia-intrinsic resistance to CD19-CAR remain poorly understood. We hypothesize that PNR leukemias are distinct compared with primary sensitive (PS) leukemias and that these differences are present before treatment. We used a multiomic approach to investigate this in 14 patients (7 with PNR and 7 with PS) enrolled in the PLAT-02 trial at Seattle Children's Hospital. Long-read PacBio sequencing helped identify 1 PNR in which 47% of CD19 transcripts had exon 2 skipping, but other samples lacked CD19 transcript abnormalities. Epigenetic profiling discovered DNA hypermethylation at genes targeted by polycomb repressive complex 2 (PRC2) in embryonic stem cells. Similarly, assays of transposase-accessible chromatin-sequencing revealed reduced accessibility at these PRC2 target genes, with a gain in accessibility of regions characteristic of hematopoietic stem cells and multilineage progenitors in PNR. Single-cell RNA sequencing and cytometry by time of flight analyses identified leukemic subpopulations expressing multilineage markers and decreased antigen presentation in PNR. We thus describe the association of a stem cell epigenome with primary resistance to CD19-CAR therapy. Future trials incorporating these biomarkers, with the addition of multispecific CAR T cells targeting against leukemic stem cell or myeloid antigens, and/or combined epigenetic therapy to disrupt this distinct stem cell epigenome may improve outcomes of patients with B-ALL.

An optimized bicistronic chimeric antigen receptor against GPC2 or CD276 overcomes heterogeneous expression in neuroblastoma.

Chimeric antigen receptor (CAR) T cell therapies targeting single antigens have performed poorly in clinical trials for solid tumors due to heterogenous expression of tumor-associated antigens (TAAs), limited T cell persistence, and T cell exhaustion. Here, we aimed to identify optimal CARs against glypican 2 (GPC2) or CD276 (B7-H3), which were highly but heterogeneously expressed in neuroblastoma (NB), a lethal extracranial solid tumor of childhood. First, we examined CAR T cell expansion in the presence of targets by digital droplet PCR. Next, using pooled competitive optimization of CAR by cellular indexing of transcriptomes and epitopes by sequencing (CITE-Seq), termed P-COCC, we simultaneously analyzed protein and transcriptome expression of CAR T cells to identify high-activity CARs. Finally, we performed cytotoxicity assays to identify the most effective CAR against each target and combined the CARs into a bicistronic "OR" CAR (BiCisCAR). BiCisCAR T cells effectively eliminated tumor cells expressing GPC2 or CD276. Furthermore, the BiCisCAR T cells demonstrated prolonged persistence and resistance to exhaustion when compared with CARs targeting a single antigen. This study illustrated that targeting multiple TAAs with BiCisCAR may overcome heterogenous expression of target antigens in solid tumors and identified a potent, clinically relevant CAR against NB. Moreover, our multimodal approach integrating competitive expansion, P-COCC, and cytotoxicity assays is an effective strategy to identify potent CARs among a pool of candidates.

Dentithecamides A-H, Diacylated Zoanthoxanthin Derivatives with PAX3-FOXO1 Inhibitory Activity from the Hydroid Dentitheca habereri.

Chemical investigation of the marine hydroid Dentitheca habereri led to the identification of eight new diacylated zoanthoxanthin alkaloids, named dentithecamides A-H (1-8), along with three previously reported analogues, zoamides B-D (9-11). The structures of compounds 1-11 were elucidated by spectroscopic and spectrometric analyses, including IR, HRESIMS, and NMR experiments, and by comparison with literature data. Compounds 1-11 are the first zoanthoxanthin alkaloids to be reported from a hydroid. Dentithecamides A (1) and B (2) along with zoamides B-D (9-11), which all share a conformationally mobile cycloheptadiene core, inhibited PAX3-FOXO1 regulated transcriptional activity and thus provided a structural framework for the potential development of more potent PAX3-FOXO1 inhibitors.

Tumor Mutation Burden, Expressed Neoantigens and the Immune Microenvironment in Diffuse Gliomas.

BACKGROUND: A consistent correlation between tumor mutation burden (TMB) and tumor immune microenvironment has not been observed in gliomas as in other cancers. METHODS: Driver germline and somatic mutations, TMB, neoantigen, and immune cell signatures were analyzed using whole exome sequencing (WES) and transcriptome sequencing of tumor and WES of matched germline DNA in a cohort of 66 glioma samples (44 IDH-mutant and 22 IDH-wildtype). RESULTS: Fourteen samples revealed a hypermutator phenotype (HMP). Eight pathogenic (P) or likely pathogenic (LP) germline variants were detected in 9 (19%) patients. Six of these 8 genes were DNA damage repair genes. P/LP germline variants were found in 22% of IDH-mutant gliomas and 12.5% of IDH-wildtype gliomas (p = 0.7). TMB was correlated with expressed neoantigen but showed an inverse correlation with immune score (R = -0.46, p = 0.03) in IDH-wildtype tumors and no correlation in IDH-mutant tumors. The Antigen Processing and Presentation (APP) score correlated with immune score and was surprisingly higher in NHMP versus HMP samples in IDH-wildtype gliomas, but higher in HMP versus NHMP in IDH-mutant gliomas. CONCLUSION: TMB was inversely correlated with immune score in IDH-wildtype gliomas and showed no correlation in IDH-mutant tumors. APP was correlated with immune score and may be further investigated as a biomarker for response to immunotherapy in gliomas. Studies of germline variants in a larger glioma cohort are warranted.

BAF complexes drive proliferation and block myogenic differentiation in fusion-positive rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of skeletal muscle lineage. The aggressive alveolar subtype is characterized by t(2;13) or t(1;13) translocations encoding for PAX3- or PAX7-FOXO1 chimeric transcription factors, respectively, and are referred to as fusion positive RMS (FP-RMS). The fusion gene alters the myogenic program and maintains the proliferative state while blocking terminal differentiation. Here, we investigated the contributions of chromatin regulatory complexes to FP-RMS tumor maintenance. We define the mSWI/SNF functional repertoire in FP-RMS. We find that SMARCA4 (encoding BRG1) is overexpressed in this malignancy compared to skeletal muscle and is essential for cell proliferation. Proteomic studies suggest proximity between PAX3-FOXO1 and BAF complexes, which is further supported by genome-wide binding profiles revealing enhancer colocalization of BAF with core regulatory transcription factors. Further, mSWI/SNF complexes localize to sites of de novo histone acetylation. Phenotypically, interference with mSWI/SNF complex function induces transcriptional activation of the skeletal muscle differentiation program associated with MYCN enhancer invasion at myogenic target genes, which is recapitulated by BRG1 targeting compounds. We conclude that inhibition of BRG1 overcomes the differentiation blockade of FP-RMS cells and may provide a therapeutic strategy for this lethal childhood tumor.

Immuno-transcriptomic profiling of extracranial pediatric solid malignancies.

We perform an immunogenomics analysis utilizing whole-transcriptome sequencing of 657 pediatric extracranial solid cancer samples representing 14 diagnoses, and additionally utilize transcriptomes of 131 pediatric cancer cell lines and 147 normal tissue samples for comparison. We describe patterns of infiltrating immune cells, T cell receptor (TCR) clonal expansion, and translationally relevant immune checkpoints. We find that tumor-infiltrating lymphocytes and TCR counts vary widely across cancer types and within each diagnosis, and notably are significantly predictive of survival in osteosarcoma patients. We identify potential cancer-specific immunotherapeutic targets for adoptive cell therapies including cell-surface proteins, tumor germline antigens, and lineage-specific transcription factors. Using an orthogonal immunopeptidomics approach, we find several potential immunotherapeutic targets in osteosarcoma and Ewing sarcoma and validated PRAME as a bona fide multi-pediatric cancer target. Importantly, this work provides a critical framework for immune targeting of extracranial solid tumors using parallel immuno-transcriptomic and -peptidomic approaches.

Genomic Classification and Clinical Outcome in Rhabdomyosarcoma: A Report From an International Consortium.

PURPOSE: Rhabdomyosarcoma is the most common soft tissue sarcoma of childhood. Despite aggressive therapy, the 5-year survival rate for patients with metastatic or recurrent disease remains poor, and beyond PAX-FOXO1 fusion status, no genomic markers are available for risk stratification. We present an international consortium study designed to determine the incidence of driver mutations and their association with clinical outcome. PATIENTS AND METHODS: Tumor samples collected from patients enrolled on Children's Oncology Group trials (1998-2017) and UK patients enrolled on malignant mesenchymal tumor and RMS2005 (1995-2016) trials were subjected to custom-capture sequencing. Mutations, indels, gene deletions, and amplifications were identified, and survival analysis was performed. RESULTS: DNA from 641 patients was suitable for analyses. A median of one mutation was found per tumor. In FOXO1 fusion-negative cases, mutation of any RAS pathway member was found in > 50% of cases, and 21% had no putative driver mutation identified. BCOR (15%), NF1 (15%), and TP53 (13%) mutations were found at a higher incidence than previously reported and TP53 mutations were associated with worse outcomes in both fusion-negative and FOXO1 fusion-positive cases. Interestingly, mutations in RAS isoforms predominated in infants < 1 year (64% of cases). Mutation of MYOD1 was associated with histologic patterns beyond those previously described, older age, head and neck primary site, and a dismal survival. Finally, we provide a searchable companion database (ClinOmics), containing all genomic variants, and clinical annotation including survival data. CONCLUSION: This is the largest genomic characterization of clinically annotated rhabdomyosarcoma tumors to date and provides prognostic genetic features that refine risk stratification and will be incorporated into prospective trials.

Outcome-Related Signatures Identified by Whole Transcriptome Sequencing of Resectable Stage III/IV Melanoma Evaluated after Starting Hu14.18-IL2.

PURPOSE: We analyzed whole transcriptome sequencing in tumors from 23 patients with stage III or IV melanoma from a pilot trial of the anti-GD2 immunocytokine, hu14.18-IL2, to identify predictive immune and/or tumor biomarkers in patients with melanoma at high risk for recurrence. EXPERIMENTAL DESIGN: Patients were randomly assigned to receive the first of three monthly courses of hu14.18-IL2 immunotherapy either before (Group A) or after (Group B) complete surgical resection of all known diseases. Tumors were evaluated by histology and whole transcriptome sequencing. RESULTS: Tumor-infiltrating lymphocyte (TIL) levels directly associated with relapse-free survival (RFS) and overall survival (OS) in resected tumors from Group A, where early responses to the immunotherapy agent could be assessed. TIL levels directly associated with a previously reported immune signature, which associated with RFS and OS, particularly in Group A tumors. In Group A tumors, there were decreased cell-cycling gene RNA transcripts, but increased RNA transcripts for repair and growth genes. We found that outcome (RFS and OS) was directly associated with several immune signatures and immune-related RNA transcripts and inversely associated with several tumor growth-associated transcripts, particularly in Group A tumors. Most of these associations were not seen in Group B tumors. CONCLUSIONS: We interpret these data to signify that both immunologic and tumoral cell processes, as measured by RNA-sequencing analyses detected shortly after initiation of hu14.18-IL2 therapy, are associated with long-term survival and could potentially be used as prognostic biomarkers in tumor resection specimens obtained after initiating neoadjuvant immunotherapy.

Case Report: Single-Cell Transcriptomic Analysis of an Anaplastic Oligodendroglioma Post Immunotherapy.

Glioma is the most common primary malignant brain tumor with a poor prognosis. Immune checkpoint inhibitors have been of great interest in investigation of glioma treatments. Here, we report single-cell transcriptomic analyses of two tumor areas from an oligodendroglioma taken from a patient who had multiple tumor recurrences, following several chemotherapies and radiation treatments. The patient subsequently received nivolumab and was considered have disease progression based on conventional diagnostic imaging after two cycles of treatment. He underwent a debulking surgical resection and pathological diagnosis was recurrent disease. During the surgery, tumor tissues were also collected from the enhancing and non-enhancing areas for a scRNAseq analysis to investigate the tumor microenvironment of these radiographically divergent areas. The scRNAseq analysis reveals a plethora of immune cells, suggesting that the increased mass observed on MRI may be partially a result of immune cell infiltration. The patient continued to receive immunotherapy after a short course of palliative radiation and remained free of disease progression for at least 12 months after the last surgery, suggesting a sustained response to immunotherapy. The scRNAseq analysis indicated that the radiological progression was in large part due to immune cell infiltrate and continued immunotherapy led to a positive clinical outcome in a patient who would have otherwise been admitted to hospice care with halting of immunotherapy. Our study demonstrates the potential of scRNAseq analyses in understanding the tumor microenvironment, which may assist the clinical decision-making process for challenging glioma cases following immunotherapy.

MEK inhibition induces MYOG and remodels super-enhancers in RAS-driven rhabdomyosarcoma.

The RAS isoforms are frequently mutated in many types of human cancers, including PAX3/PAX7 fusion-negative rhabdomyosarcoma. Pediatric RMS arises from skeletal muscle progenitor cells that have failed to differentiate normally. The role of mutant RAS in this differentiation blockade is incompletely understood. We demonstrate that oncogenic RAS, acting through the RAF-MEK [mitogen-activated protein kinase (MAPK) kinase]-ERK (extracellular signal-regulated kinase) MAPK effector pathway, inhibits myogenic differentiation in rhabdomyosarcoma by repressing the expression of the prodifferentiation myogenic transcription factor, MYOG. This repression is mediated by ERK2-dependent promoter-proximal stalling of RNA polymerase II at the MYOG locus. Small-molecule screening with a library of mechanistically defined inhibitors showed that RAS-driven RMS is vulnerable to MEK inhibition. MEK inhibition with trametinib leads to the loss of ERK2 at the MYOG promoter and releases the transcriptional stalling of MYOG expression. MYOG subsequently opens chromatin and establishes super-enhancers at genes required for late myogenic differentiation. Furthermore, trametinib, in combination with an inhibitor of IGF1R, potently decreases rhabdomyosarcoma cell viability and slows tumor growth in xenograft models. Therefore, this combination represents a potential therapeutic for RAS-mutated rhabdomyosarcoma.

Clinically Relevant Cytotoxic Immune Cell Signatures and Clonal Expansion of T-Cell Receptors in High-Risk MYCN-Not-Amplified Human Neuroblastoma.

Purpose: High-risk neuroblastoma is an aggressive disease. DNA sequencing studies have revealed a paucity of actionable genomic alterations and a low mutation burden, posing challenges to develop effective novel therapies. We used RNA sequencing (RNA-seq) to investigate the biology of this disease, including a focus on tumor-infiltrating lymphocytes (TIL).Experimental Design: We performed deep RNA-seq on pretreatment diagnostic tumors from 129 high-risk and 21 low- or intermediate-risk patients with neuroblastomas. We used single-sample gene set enrichment analysis to detect gene expression signatures of TILs in tumors and examined their association with clinical and molecular parameters, including patient outcome. The expression profiles of 190 additional pretreatment diagnostic neuroblastomas, a neuroblastoma tissue microarray, and T-cell receptor (TCR) sequencing were used to validate our findings.Results: We found that MYCN-not-amplified (MYCN-NA) tumors had significantly higher cytotoxic TIL signatures compared with MYCN-amplified (MYCN-A) tumors. A reported MYCN activation signature was significantly associated with poor outcome for high-risk patients with MYCN-NA tumors; however, a subgroup of these patients who had elevated activated natural killer (NK) cells, CD8+ T cells, and cytolytic signatures showed improved outcome and expansion of infiltrating TCR clones. Furthermore, we observed upregulation of immune exhaustion marker genes, indicating an immune-suppressive microenvironment in these neuroblastomas.Conclusions: This study provides evidence that RNA signatures of cytotoxic TIL are associated with the presence of activated NK/T cells and improved outcomes in high-risk neuroblastoma patients harboring MYCN-NA tumors. Our findings suggest that these high-risk patients with MYCN-NA neuroblastoma may benefit from additional immunotherapies incorporated into the current therapeutic strategies. Clin Cancer Res; 24(22); 5673-84. ©2018 AACR.

Frequent inactivating germline mutations in DNA repair genes in patients with Ewing sarcoma.

PURPOSE: Ewing sarcoma is a small round blue cell tumor that is highly malignant and predominantly affects the adolescent and young adult population. It has long been suspected that a genetic predisposition exists for this cancer, but the germ-line genetic underpinnings of this disease have not been well established. METHODS: We performed germline variant analysis of whole-genome or whole-exome sequencing of samples from 175 patients affected by Ewing sarcoma. RESULTS: We discovered pathogenic or likely pathogenic germline mutations in 13.1% of our cohort. Pathogenic mutations were highly enriched for genes involved with DNA damage repair and for genes associated with cancer predisposition syndromes. CONCLUSION: Our findings reported here have important clinical implications for patients and families affected by Ewing sarcoma. Genetic counseling should be considered for patients and families affected by this disease to take advantage of existing risk management strategies. Our study also highlights the importance of germline sequencing for patients enrolled in precision-medicine protocols.Genet Med advance online publication 26 January 2017.

CD19 CAR immune pressure induces B-precursor acute lymphoblastic leukaemia lineage switch exposing inherent leukaemic plasticity.

Adoptive immunotherapy using chimeric antigen receptor (CAR) expressing T cells targeting the CD19 B lineage receptor has demonstrated marked success in relapsed pre-B-cell acute lymphoblastic leukaemia (ALL). Persisting CAR-T cells generate sustained pressure against CD19 that may drive unique mechanisms of resistance. Pre-B ALL originates from a committed pre-B cell or an earlier progenitor, with potential to reprogram into other hematopoietic lineages. Here we report changes in lineage markers including myeloid conversion in patients following CD19 CAR therapy. Using murine ALL models we study the long-term effects of CD19 CAR-T cells and demonstrate partial or complete lineage switch as a consistent mechanism of CAR resistance depending on the underlying genetic oncogenic driver. Deletion of Pax5 or Ebf1 recapitulates lineage reprogramming occurring during CD19 CAR pressure. Our findings establish lineage switch as a mechanism of CAR resistance exposing inherent plasticity in genetic subtypes of pre-B-cell ALL.

MultiDimensional ClinOmics for Precision Therapy of Children and Adolescent Young Adults with Relapsed and Refractory Cancer: A Report from the Center for Cancer Research.

PURPOSE: We undertook a multidimensional clinical genomics study of children and adolescent young adults with relapsed and refractory cancers to determine the feasibility of genome-guided precision therapy. EXPERIMENTAL DESIGN: Patients with non-central nervous system solid tumors underwent a combination of whole exome sequencing (WES), whole transcriptome sequencing (WTS), and high-density single-nucleotide polymorphism array analysis of the tumor, with WES of matched germline DNA. Clinically actionable alterations were identified as a reportable germline mutation, a diagnosis change, or a somatic event (including a single nucleotide variant, an indel, an amplification, a deletion, or a fusion gene), which could be targeted with drugs in existing clinical trials or with FDA-approved drugs. RESULTS: Fifty-nine patients in 20 diagnostic categories were enrolled from 2010 to 2014. Ages ranged from 7 months to 25 years old. Seventy-three percent of the patients had prior chemotherapy, and the tumors from these patients with relapsed or refractory cancers had a higher mutational burden than that reported in the literature. Thirty patients (51% of total) had clinically actionable mutations, of which 24 (41%) had a mutation that was currently targetable in a clinical trial setting, 4 patients (7%) had a change in diagnosis, and 7 patients (12%) had a reportable germline mutation. CONCLUSIONS: We found a remarkably high number of clinically actionable mutations in 51% of the patients, and 12% with significant germline mutations. We demonstrated the clinical feasibility of next-generation sequencing in a diverse population of relapsed and refractory pediatric solid tumors. Clin Cancer Res; 22(15); 3810-20. ©2016 AACR.

MYCN controls an alternative RNA splicing program in high-risk metastatic neuroblastoma.

The molecular mechanisms underlying the aggressive behavior of MYCN driven neuroblastoma (NBL) is under intense investigation; however, little is known about the impact of this family of transcription factors on the splicing program. Here we used high-throughput RNA sequencing to systematically study the expression of RNA isoforms in stage 4 MYCN-amplified NBL, an aggressive subtype of metastatic NBL. We show that MYCN-amplified NBL tumors display a distinct gene splicing pattern affecting multiple cancer hallmark functions. Six splicing factors displayed unique differential expression patterns in MYCN-amplified tumors and cell lines, and the binding motifs for some of these splicing factors are significantly enriched in differentially-spliced genes. Direct binding of MYCN to promoter regions of the splicing factors PTBP1 and HNRNPA1 detected by ChIP-seq demonstrates that MYCN controls the splicing pattern by direct regulation of the expression of these key splicing factors. Furthermore, high expression of PTBP1 and HNRNPA1 was significantly associated with poor overall survival of stage4 NBL patients (p ≤ 0.05). Knocking down PTBP1, HNRNPA1 and their downstream target PKM2, an isoform of pro-tumor-growth, result in repressed growth of NBL cells. Therefore, our study reveals a novel role of MYCN in controlling global splicing program through regulation of splicing factors in addition to its well-known role in the transcription program. These findings suggest a therapeutically potential to target the key splicing factors or gene isoforms in high-risk NBL with MYCN-amplification.