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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Translational Strategies for Repotrectinib in Neuroblastoma.

Limited clinical data are available regarding the utility of multikinase inhibition in neuroblastoma. Repotrectinib (TPX-0005) is a multikinase inhibitor that targets ALK, TRK, JAK2/STAT, and Src/FAK, which have all been implicated in the pathogenesis of neuroblastoma. We evaluated the preclinical activity of repotrectinib monotherapy and in combination with chemotherapy as a potential therapeutic approach for relapsed/refractory neuroblastoma. In vitro sensitivity to repotrectinib, ensartinib, and cytotoxic chemotherapy was evaluated in neuroblastoma cell lines. In vivo antitumor effect of repotrectinib monotherapy, and in combination with chemotherapy, was evaluated using a genotypically diverse cohort of patient-derived xenograft (PDX) models of neuroblastoma. Repotrectinib had comparable cytotoxic activity across cell lines irrespective of ALK mutational status. Combination with chemotherapy demonstrated increased antiproliferative activity across several cell lines. Repotrectinib monotherapy had notable antitumor activity and prolonged event-free survival compared with vehicle and ensartinib in PDX models (P < 0.05). Repotrectinib plus chemotherapy was superior to chemotherapy alone in ALK-mutant and ALK wild-type PDX models. These results demonstrate that repotrectinib has antitumor activity in genotypically diverse neuroblastoma models, and that combination of a multikinase inhibitor with chemotherapy may be a promising treatment paradigm for translation to the clinic.

Prospective pan-cancer germline testing using MSK-IMPACT informs clinical translation in 751 patients with pediatric solid tumors.

The spectrum of germline predisposition in pediatric cancer continues to be realized. Here we report 751 solid tumor patients who underwent prospective matched tumor-normal DNA sequencing and downstream clinical use (clinicaltrials.gov NCT01775072). Germline pathogenic and likely pathogenic (P/LP) variants were reported. One or more P/LP variants were found in 18% (138/751) of individuals when including variants in low, moderate, and high penetrance dominant or recessive genes, or 13% (99/751) in moderate and high penetrance dominant genes. 34% of high or moderate penetrance variants were unexpected based on the patient's diagnosis and previous history. 76% of patients with positive results completed a clinical genetics visit, and 21% had at least one relative undergo cascade testing as a result of this testing. Clinical actionability additionally included screening, risk reduction in relatives, reproductive use, and use of targeted therapies. Germline testing should be considered for all children with cancer.

STAG2 loss rewires oncogenic and developmental programs to promote metastasis in Ewing sarcoma.

The core cohesin subunit STAG2 is recurrently mutated in Ewing sarcoma but its biological role is less clear. Here, we demonstrate that cohesin complexes containing STAG2 occupy enhancer and polycomb repressive complex (PRC2)-marked regulatory regions. Genetic suppression of STAG2 leads to a compensatory increase in cohesin-STAG1 complexes, but not in enhancer-rich regions, and results in reprogramming of cis-chromatin interactions. Strikingly, in STAG2 knockout cells the oncogenic genetic program driven by the fusion transcription factor EWS/FLI1 was highly perturbed, in part due to altered enhancer-promoter contacts. Moreover, loss of STAG2 also disrupted PRC2-mediated regulation of gene expression. Combined, these transcriptional changes converged to modulate EWS/FLI1, migratory, and neurodevelopmental programs. Finally, consistent with clinical observations, functional studies revealed that loss of STAG2 enhances the metastatic potential of Ewing sarcoma xenografts. Our findings demonstrate that STAG2 mutations can alter chromatin architecture and transcriptional programs to promote an aggressive cancer phenotype.