Roadmap for the next generation of Children's Oncology Group rhabdomyosarcoma trials.

Clinical trials conducted by the Intergroup Rhabdomyosarcoma (RMS) Study Group and the Children's Oncology Group have been pivotal to establishing current standards for diagnosis and therapy for RMS. Recent advancements in understanding the biology and clinical behavior of RMS have led to more nuanced approaches to diagnosis, risk stratification, and treatment. The complexities introduced by these advancements, coupled with the rarity of RMS, pose challenges to conducting large-scale phase 3 clinical trials to evaluate new treatment strategies for RMS. Given these challenges, systematic planning of future clinical trials in RMS is paramount to address pertinent questions regarding the therapeutic efficacy of drugs, biomarkers of response, treatment-related toxicity, and patient quality of life. Herein, the authors outline the proposed strategic approach of the Children's Oncology Group Soft Tissue Sarcoma Committee to the next generation of RMS clinical trials, focusing on five themes: improved novel agent identification and preclinical to clinical translation, more efficient trial development and implementation, expanded opportunities for knowledge generation during trials, therapeutic toxicity reduction and quality of life, and patient engagement.

Bromodomain 4 inhibition leads to MYCN downregulation in Wilms tumor.

BACKGROUND: Wilms tumor is the most common childhood kidney cancer. Two distinct histological subtypes of Wilms tumor have been described: tumors lacking anaplasia (the favorable subtype) and tumors displaying anaplastic features (the unfavorable subtype). Children with favorable disease generally have a very good prognosis, whereas those with anaplasia are oftentimes refractory to standard treatments and suffer poor outcomes, leading to an unmet clinical need. MYCN dysregulation has been associated with a number of pediatric cancers including Wilms tumor. PROCEDURES: In this context, we undertook a functional genomics approach to uncover novel therapeutic strategies for those patients with anaplastic Wilms tumor. Genomic analysis and in vitro experimentation demonstrate that cell growth can be reduced by modulating MYCN overexpression via bromodomain 4 (BRD4) inhibition in both anaplastic and nonanaplastic Wilms tumor models. RESULTS: We observed a time-dependent reduction of MYCN and MYCC protein levels upon BRD4 inhibition in Wilms tumor cell lines, which led to cell death and proliferation suppression. BRD4 inhibition significantly reduced tumor volumes in Wilms tumor patient-derived xenograft (PDX) mouse models. CONCLUSIONS: We suggest that AZD5153, a novel dual-BRD4 inhibitor, can reduce MYCN levels in both anaplastic and nonanaplastic Wilms tumor cell lines, reduces tumor volume in Wilms tumor PDXs, and should be further explored for its therapeutic potential.

Recurrent EML4-NTRK3 fusions in infantile fibrosarcoma and congenital mesoblastic nephroma suggest a revised testing strategy.

Infantile fibrosarcoma and congenital mesoblastic nephroma are tumors of infancy traditionally associated with the ETV6-NTRK3 gene fusion. However, a number of case reports have identified variant fusions in these tumors. In order to assess the frequency of variant NTRK3 fusions, and in particular whether the recently identified EML4-NTRK3 fusion is recurrent, 63 archival cases of infantile fibrosarcoma, congenital mesoblastic nephroma, mammary analog secretory carcinoma and secretory breast carcinoma (tumor types that are known to carry recurrent ETV6-NTRK3 fusions) were tested with NTRK3 break-apart FISH, EML4-NTRK3 dual fusion FISH, and targeted RNA sequencing. The EML4-NTRK3 fusion was identified in two cases of infantile fibrosarcoma (one of which was previously described), and in one case of congenital mesoblastic nephroma, demonstrating that the EML4-NTRK3 fusion is a recurrent genetic event in these related tumors. The growing spectrum of gene fusions associated with infantile fibrosarcoma and congenital mesoblastic nephroma along with the recent availability of targeted therapies directed toward inhibition of NTRK signaling argue for alternate testing strategies beyond ETV6 break-apart FISH. The use of either NTRK3 FISH or next-generation sequencing will expand the number of cases in which an oncogenic fusion is identified and facilitate optimal diagnosis and treatment for patients.

Dynamic network curvature analysis of gene expression reveals novel potential therapeutic targets in sarcoma.

Network properties account for the complex relationship between genes, making it easier to identify complex patterns in their interactions. In this work, we leveraged these network properties for dual purposes. First, we clustered pediatric sarcoma tumors using network information flow as a similarity metric, computed by the Wasserstein distance. We demonstrate that this approach yields the best concordance with histological subtypes, validated against three state-of-the-art methods. Second, to identify molecular targets that would be missed by more conventional methods of analysis, we applied a novel unsupervised method to cluster gene interactomes represented as networks in pediatric sarcoma. RNA-Seq data were mapped to protein-level interactomes to construct weighted networks that were then subjected to a non-Euclidean, multi-scale geometric approach centered on a discrete notion of curvature. This provides a measure of the functional association among genes in the context of their connectivity. In confirmation of the validity of this method, hierarchical clustering revealed the characteristic EWSR1-FLI1 fusion in Ewing sarcoma. Furthermore, assessing the effects of in silico edge perturbations and simulated gene knockouts as quantified by changes in curvature, we found non-trivial gene associations not previously identified.

Overcoming Clinical Resistance to EZH2 Inhibition Using Rational Epigenetic Combination Therapy.

Epigenetic dependencies have become evident in many cancers. On the basis of antagonism between BAF/SWI-SNF and PRC2 in SMARCB1-deficient sarcomas, we recently completed the clinical trial of the EZH2 inhibitor tazemetostat. However, the principles of tumor response to epigenetic therapy in general, and tazemetostat in particular, remain unknown. Using functional genomics and diverse experimental models, we define molecular mechanisms of tazemetostat resistance in SMARCB1-deficient tumors. We found distinct acquired mutations that converge on the RB1/E2F axis and decouple EZH2-dependent differentiation and cell-cycle control. This allows tumor cells to escape tazemetostat-induced G1 arrest, suggests a general mechanism for effective therapy, and provides prospective biomarkers for therapy stratification, including PRICKLE1. On the basis of this, we develop a combination strategy to circumvent tazemetostat resistance using bypass targeting of AURKB. This offers a paradigm for rational epigenetic combination therapy suitable for translation to clinical trials for epithelioid sarcomas, rhabdoid tumors, and other epigenetically dysregulated cancers. SIGNIFICANCE: Genomic studies of patient epithelioid sarcomas and rhabdoid tumors identify mutations converging on a common pathway for response to EZH2 inhibition. Resistance mutations decouple drug-induced differentiation from cell-cycle control. We identify an epigenetic combination strategy to overcome resistance and improve durability of response, supporting its investigation in clinical trials. See related commentary by Paolini and Souroullas, p. 903. This article is featured in Selected Articles from This Issue, p. 897.

Identification and Pharmacological Targeting of Treatment-Resistant, Stem-like Breast Cancer Cells for Combination Therapy.

Tumors frequently harbor isogenic yet epigenetically distinct subpopulations of multi-potent cells with high tumor-initiating potential-often called Cancer Stem-Like Cells (CSLCs). These can display preferential resistance to standard-of-care chemotherapy. Single-cell analyses can help elucidate Master Regulator (MR) proteins responsible for governing the transcriptional state of these cells, thus revealing complementary dependencies that may be leveraged via combination therapy. Interrogation of single-cell RNA sequencing profiles from seven metastatic breast cancer patients, using perturbational profiles of clinically relevant drugs, identified drugs predicted to invert the activity of MR proteins governing the transcriptional state of chemoresistant CSLCs, which were then validated by CROP-seq assays. The top drug, the anthelmintic albendazole, depleted this subpopulation in vivo without noticeable cytotoxicity. Moreover, sequential cycles of albendazole and paclitaxel-a commonly used chemotherapeutic -displayed significant synergy in a patient-derived xenograft (PDX) from a TNBC patient, suggesting that network-based approaches can help develop mechanism-based combinatorial therapies targeting complementary subpopulations.

Epigenetic targeting of PGBD5-dependent DNA damage in SMARCB1-deficient sarcomas.

Despite the potential of targeted epigenetic therapies, most cancers do not respond to current epigenetic drugs. The Polycomb repressive complex EZH2 inhibitor tazemetostat was recently approved for the treatment of SMARCB1-deficient epithelioid sarcomas, based on the functional antagonism between PRC2 and loss of SMARCB1. Through the analysis of tazemetostat-treated patient tumors, we recently defined key principles of their response and resistance to EZH2 epigenetic therapy. Here, using transcriptomic inference from SMARCB1-deficient tumor cells, we nominate the DNA damage repair kinase ATR as a target for rational combination EZH2 epigenetic therapy. We show that EZH2 inhibition promotes DNA damage in epithelioid and rhabdoid tumor cells, at least in part via its induction of the transposase-derived PGBD5. We leverage this collateral synthetic lethal dependency to target PGBD5-dependent DNA damage by inhibition of ATR but not CHK1 using elimusertib. Consequently, combined EZH2 and ATR inhibition improves therapeutic responses in diverse patient-derived epithelioid and rhabdoid tumors in vivo. This advances a combination epigenetic therapy based on EZH2-PGBD5 synthetic lethal dependency suitable for immediate translation to clinical trials for patients.

Pediatric Leukemia Roadmaps Are a Guide for Positive Metastatic Bone Sarcoma Trials.

ALL cures require many MRD therapies. This strategy should drive experiments and trials in metastatic bone sarcomas.

Chemotherapy use in elderly patients with soft tissue sarcoma: a population-based study.

Adjuvant chemotherapy for soft tissue sarcoma (STS) remains controversial while improvement in survival has never been conclusively demonstrated for metastatic STS. We identified individuals in SEER-Medicare with resected or metastatic STS, 1991-2007. Of 2,382 patients with resected STS, 106 (4.5%) received chemotherapy. High tumor grade, larger tumor size, and malignant fibrous histiocytoma subtype were associated with chemotherapy receipt. Of 365 patients with metastatic STS, 118 (32.4%) received chemotherapy. Younger age, fewer comorbidities, and being married were associated with chemotherapy receipt. Consistent with clinical trials, we found minimal use of chemotherapy. Clinical factors were associated with chemotherapy receipt in nonmetastatic STS.

A transgenic, mesodermal specific, Dkk1 mouse model recapitulates a spectrum of human congenital limb reduction defects.

Congenital limb reduction defects occurring in isolation of other developmental abnormalities continue to be an important medical problem in which little progress has been made. Herein we generated transgenic mice expressing Dkk1 in an appendicular mesodermal pattern. Prx1-Dkk1 mice recapitulate a full spectrum of human congenital limb reduction defects, without other developmental issues, and have normal life-spans. Importantly, a close examination of the inheritance pattern suggests that there is a significant degree of incomplete penetrance as progeny of phenotypically positive or phenotypically negative, but genotypically positive Prx1-Dkk1 mice, consistently give rise to both phenotypically positive mice and phenotypically normal-appearing mice. Thus, this heterogeneous phenotype is reproducible with each generation regardless of the phenotype of the parents. We further go on to identify that mesenchymal stem cells from Prx1-Dkk1 mice have limited proliferative ability, but normal differentiation potential, which may explain the mechanism for the limb reduction defects observed. We believe Prx1-Dkk1 mice may prove useful in the future to study the mechanisms underlying the development of congenital limb reduction defects.

Overcoming clinical resistance to EZH2 inhibition using rational epigenetic combination therapy.

Essential epigenetic dependencies have become evident in many cancers. Based on the functional antagonism between BAF/SWI/SNF and PRC2 in SMARCB1-deficient sarcomas, we and colleagues recently completed the clinical trial of the EZH2 inhibitor tazemetostat. However, the principles of tumor response to epigenetic therapy in general, and tazemetostat in particular, remain unknown. Using functional genomics of patient tumors and diverse experimental models, we sought to define molecular mechanisms of tazemetostat resistance in SMARCB1-deficient sarcomas and rhabdoid tumors. We found distinct classes of acquired mutations that converge on the RB1/E2F axis and decouple EZH2-dependent differentiation and cell cycle control. This allows tumor cells to escape tazemetostat-induced G1 arrest despite EZH2 inhibition, and suggests a general mechanism for effective EZH2 therapy. This also enables us to develop combination strategies to circumvent tazemetostat resistance using cell cycle bypass targeting via AURKB, and synthetic lethal targeting of PGBD5-dependent DNA damage repair via ATR. This reveals prospective biomarkers for therapy stratification, including PRICKLE1 associated with tazemetostat resistance. In all, this work offers a paradigm for rational epigenetic combination therapy suitable for immediate translation to clinical trials for epithelioid sarcomas, rhabdoid tumors, and other epigenetically dysregulated cancers.

Clonal evolution during metastatic spread in high-risk neuroblastoma.

Patients with high-risk neuroblastoma generally present with widely metastatic disease and often relapse despite intensive therapy. As most studies to date focused on diagnosis-relapse pairs, our understanding of the genetic and clonal dynamics of metastatic spread and disease progression remain limited. Here, using genomic profiling of 470 sequential and spatially separated samples from 283 patients, we characterize subtype-specific genetic evolutionary trajectories from diagnosis through progression and end-stage metastatic disease. Clonal tracing timed disease initiation to embryogenesis. Continuous acquisition of structural variants at disease-defining loci (MYCN, TERT, MDM2-CDK4) followed by convergent evolution of mutations targeting shared pathways emerged as the predominant feature of progression. At diagnosis metastatic clones were already established at distant sites where they could stay dormant, only to cause relapses years later and spread via metastasis-to-metastasis and polyclonal seeding after therapy.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma.

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma.

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease. SIGNIFICANCE: The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

A Transcriptome-Based Precision Oncology Platform for Patient-Therapy Alignment in a Diverse Set of Treatment-Resistant Malignancies.

Predicting in vivo response to antineoplastics remains an elusive challenge. We performed a first-of-kind evaluation of two transcriptome-based precision cancer medicine methodologies to predict tumor sensitivity to a comprehensive repertoire of clinically relevant oncology drugs, whose mechanism of action we experimentally assessed in cognate cell lines. We enrolled patients with histologically distinct, poor-prognosis malignancies who had progressed on multiple therapies, and developed low-passage, patient-derived xenograft models that were used to validate 35 patient-specific drug predictions. Both OncoTarget, which identifies high-affinity inhibitors of individual master regulator (MR) proteins, and OncoTreat, which identifies drugs that invert the transcriptional activity of hyperconnected MR modules, produced highly significant 30-day disease control rates (68% and 91%, respectively). Moreover, of 18 OncoTreat-predicted drugs, 15 induced the predicted MR-module activity inversion in vivo. Predicted drugs significantly outperformed antineoplastic drugs selected as unpredicted controls, suggesting these methods may substantively complement existing precision cancer medicine approaches, as also illustrated by a case study. SIGNIFICANCE: Complementary precision cancer medicine paradigms are needed to broaden the clinical benefit realized through genetic profiling and immunotherapy. In this first-in-class application, we introduce two transcriptome-based tumor-agnostic systems biology tools to predict drug response in vivo. OncoTarget and OncoTreat are scalable for the design of basket and umbrella clinical trials. This article is highlighted in the In This Issue feature, p. 1275.

A Phase I study of gemtuzumab ozogamicin (GO) in combination with busulfan and cyclophosphamide (Bu/Cy) and allogeneic stem cell transplantation in children with poor-risk CD33+ AML: a new targeted immunochemotherapy myeloablative conditioning (MAC) regimen.

Children with high-risk acute myelogenous leukemia (AML) (induction failure [IF], refractory relapse [RR], third complete remission [CR3]) have dismal outcomes. Over 80% of AML patients express CD33, a target of gemtuzumab ozogamicin (GO). GO is an active drug in childhood AML but has not been studied in a myeloablative conditioning regimen. We sought to determine the safety of GO in combination with busulfan/cyclophosphamide (Bu/Cy) conditioning before allogeneic hematopoietic stem cell transplantation (alloSCT). GO was administered on day -14 at doses of 3.0, 4.5, 6.0, and 7.5 mg/m(2), busulfan on days -7, -6, -5, -4 (12.8-16.0 mg/kg), and cyclophosphamide on days -3 and -2 (60 mg/kg/day). GVHD prophylaxis consisted of tacrolimus and mycophenolate mofetil. We enrolled 12 patients: 8 IF, 3 RR, 1 CR3; median age: 3 years (1-17); median follow-up: 1379 days (939-2305). Nine received umbilical cord blood (UCB), 2 matched unrelated donors (MUDs) and 1 HLA-matched sibling donor: 3 patients each at GO doses of 3.0, 4.5, 6.0, or 7.5 mg/m(2). No dose-limiting toxicities secondary to GO were observed. Day 100 treatment-related mortality (TRM) was 0%. Myeloid and platelet engraftment was observed in 92% and 75% of patients at median day 22 (12-40) and 42 (21-164), respectively. Median day +30 donor chimerism was 99% (85%-100%). The probability of grade II-IV acute graft-versus-host disease (aGVHD) was 42% and chronic GVHD (cGVHD) was 28%. One-year overall survival (OS) and event-free survival (EFS) was 50% (95% confidence interval [CI], 20.8-73.6). GO combined with Bu/Cy regimen followed by alloSCT is well tolerated in children with poor-risk AML. GO at 7.5 mg/m(2) in combination with Bu/Cy is currently being tested in a phase II study.

Validation of a non-oncogene encoded vulnerability to exportin 1 inhibition in pediatric renal tumors.

BACKGROUND: Malignant rhabdoid tumors (MRTs) and Wilms' tumors (WTs) are rare and aggressive renal tumors of infants and young children comprising ∼5% of all pediatric cancers. MRTs are among the most genomically stable cancers, and although WTs are genomically heterogeneous, both generally lack therapeutically targetable genetic mutations. METHODS: Comparative protein activity analysis of MRTs (n = 68) and WTs (n = 132) across TCGA and TARGET cohorts, using metaVIPER, revealed elevated exportin 1 (XPO1) inferred activity. In vitro studies were performed on a panel of MRT and WT cell lines to evaluate effects on proliferation and cell-cycle progression following treatment with the selective XPO1 inhibitor selinexor. In vivo anti-tumor activity was assessed in patient-derived xenograft (PDX) models of MRTs and WTs. FINDINGS: metaVIPER analysis identified markedly aberrant activation of XPO1 in MRTs and WTs compared with other tumor types. All MRT and most WT cell lines demonstrated baseline, aberrant XPO1 activity with in vitro sensitivity to selinexor via cell-cycle arrest and induction of apoptosis. In vivo, XPO1 inhibitors significantly abrogated tumor growth in PDX models, inducing effective disease control with sustained treatment. Corroborating human relevance, we present a case report of a child with multiply relapsed WTs with prolonged disease control on selinexor. CONCLUSIONS: We report on a novel systems-biology-based comparative framework to identify non-genetically encoded vulnerabilities in genomically quiescent pediatric cancers. These results have provided preclinical rationale for investigation of XPO1 inhibitors in an upcoming investigator-initiated clinical trial of selinexor in children with MRTs and WTs and offer opportunities for exploration of inferred XPO1 activity as a potential predictive biomarker for response. FUNDING: This work was funded by CureSearch for Children's Cancer, Alan B. Slifka Foundation, NIH (U01 CA217858, S10 OD012351, and S10 OD021764), Michael's Miracle Cure, Hyundai Hope on Wheels, Cannonball Kids Cancer, Conquer Cancer the ASCO Foundation, Cycle for Survival, Paulie Strong Foundation, and the Grayson Fund.

Single-cell analysis and functional characterization uncover the stem cell hierarchies and developmental origins of rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a common childhood cancer that shares features with developing skeletal muscle. Yet, the conservation of cellular hierarchy with human muscle development and the identification of molecularly defined tumor-propagating cells has not been reported. Using single-cell RNA-sequencing, DNA-barcode cell fate mapping and functional stem cell assays, we uncovered shared tumor cell hierarchies in RMS and human muscle development. We also identified common developmental stages at which tumor cells become arrested. Fusion-negative RMS cells resemble early myogenic cells found in embryonic and fetal development, while fusion-positive RMS cells express a highly specific gene program found in muscle cells transiting from embryonic to fetal development at 7-7.75 weeks of age. Fusion-positive RMS cells also have neural pathway-enriched states, suggesting less-rigid adherence to muscle-lineage hierarchies. Finally, we identified a molecularly defined tumor-propagating subpopulation in fusion-negative RMS that shares remarkable similarity to bi-potent, muscle mesenchyme progenitors that can make both muscle and osteogenic cells.

Therapeutic targeting of ATR in alveolar rhabdomyosarcoma.

Despite advances in multi-modal treatment approaches, clinical outcomes of patients suffering from PAX3-FOXO1 fusion oncogene-expressing alveolar rhabdomyosarcoma (ARMS) remain dismal. Here we show that PAX3-FOXO1-expressing ARMS cells are sensitive to pharmacological ataxia telangiectasia and Rad3 related protein (ATR) inhibition. Expression of PAX3-FOXO1 in muscle progenitor cells is not only sufficient to increase sensitivity to ATR inhibition, but PAX3-FOXO1-expressing rhabdomyosarcoma cells also exhibit increased sensitivity to structurally diverse inhibitors of ATR. Mechanistically, ATR inhibition leads to replication stress exacerbation, decreased BRCA1 phosphorylation and reduced homologous recombination-mediated DNA repair pathway activity. Consequently, ATR inhibitor treatment increases sensitivity of ARMS cells to PARP1 inhibition in vitro, and combined treatment with ATR and PARP1 inhibitors induces complete regression of primary patient-derived ARMS xenografts in vivo. Lastly, a genome-wide CRISPR activation screen (CRISPRa) in combination with transcriptional analyses of ATR inhibitor resistant ARMS cells identifies the RAS-MAPK pathway and its targets, the FOS gene family, as inducers of resistance to ATR inhibition. Our findings provide a rationale for upcoming biomarker-driven clinical trials of ATR inhibitors in patients suffering from ARMS.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition.

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.