Despite gains in survival, outcomes for patients with metastatic or recurrent rhabdomyosarcoma (RMS) remain dismal. In a collaboration between the National Cancer Institute, Children's Oncology Group, and Broad Institute, we performed whole-genome, whole-exome and transcriptome sequencing to characterize the landscape of somatic alterations in 147 tumor/normal pairs. Two genotypes are evident in RMS tumors; those characterized by the PAX3 or PAX7 fusion and those that lack these fusions but harbor mutations in key signaling pathways. The overall burden of somatic mutations in RMS is relatively low, especially in tumors that harbor a PAX3/7 gene fusion. In addition to previously reported mutations of NRAS, KRAS, HRAS, FGFR4, PIK3CA, CTNNB1, we found novel recurrent mutations in FBXW7, and BCOR providing potential new avenues for therapeutic intervention. Furthermore, alteration of the receptor tyrosine kinase/RAS/PIK3CA axis affects 93% of cases providing a framework for genomics directed therapies that might improve outcomes for RMS patients.

The Ewing sarcoma family of tumors (EFT) is a group of highly malignant small round blue cell tumors occurring in children and young adults. We report here the largest genomic survey to date of 101 EFT (65 tumors and 36 cell lines). Using a combination of whole genome sequencing and targeted sequencing approaches, we discover that EFT has a very low mutational burden (0.15 mutations/Mb) but frequent deleterious mutations in the cohesin complex subunit STAG2 (21.5% tumors, 44.4% cell lines), homozygous deletion of CDKN2A (13.8% and 50%) and mutations of TP53 (6.2% and 71.9%). We additionally note an increased prevalence of the BRCA2 K3326X polymorphism in EFT patient samples (7.3%) compared to population data (OR 7.1, p = 0.006). Using whole transcriptome sequencing, we find that 11% of tumors pathologically diagnosed as EFT lack a typical EWSR1 fusion oncogene and that these tumors do not have a characteristic Ewing sarcoma gene expression signature. We identify samples harboring novel fusion genes including FUS-NCATc2 and CIC-FOXO4 that may represent distinct small round blue cell tumor variants. In an independent EFT tissue microarray cohort, we show that STAG2 loss as detected by immunohistochemistry may be associated with more advanced disease (p = 0.15) and a modest decrease in overall survival (p = 0.10). These results significantly advance our understanding of the genomic and molecular underpinnings of Ewing sarcoma and provide a foundation towards further efforts to improve diagnosis, prognosis, and precision therapeutics testing.

BackgroundTumor mutational burden (TMB), defined as the number of somatic mutations per megabase of interrogated genomic sequence, demonstrates predictive biomarker potential for the identification of patients with cancer most likely to respond to immune checkpoint inhibitors. TMB is optimally calculated by whole exome sequencing (WES), but next-generation sequencing targeted panels provide TMB estimates in a time-effective and cost-effective manner. However, differences in panel size and gene coverage, in addition to the underlying bioinformatics pipelines, are known drivers of variability in TMB estimates across laboratories. By directly comparing panel-based TMB estimates from participating laboratories, this study aims to characterize the theoretical variability of panel-based TMB estimates, and provides guidelines on TMB reporting, analytic validation requirements and reference standard alignment in order to maintain consistency of TMB estimation across platforms.MethodsEleven laboratories used WES data from The Cancer Genome Atlas Multi-Center Mutation calling in Multiple Cancers (MC3) samples and calculated TMB from the subset of the exome restricted to the genes covered by their targeted panel using their own bioinformatics pipeline (panel TMB). A reference TMB value was calculated from the entire exome using a uniform bioinformatics pipeline all members agreed on (WES TMB). Linear regression analyses were performed to investigate the relationship between WES and panel TMB for all 32 cancer types combined and separately. Variability in panel TMB values at various WES TMB values was also quantified using 95% prediction limits.ResultsStudy results demonstrated that variability within and between panel TMB values increases as the WES TMB values increase. For each panel, prediction limits based on linear regression analyses that modeled panel TMB as a function of WES TMB were calculated and found to approximately capture the intended 95% of observed panel TMB values. Certain cancer types, such as uterine, bladder and colon cancers exhibited greater variability in panel TMB values, compared with lung and head and neck cancers.ConclusionsIncreasing uptake of TMB as a predictive biomarker in the clinic creates an urgent need to bring stakeholders together to agree on the harmonization of key aspects of panel-based TMB estimation, such as the standardization of TMB reporting, standardization of analytical validation studies and the alignment of panel-based TMB values with a reference standard. These harmonization efforts should improve consistency and reliability of panel TMB estimates and aid in clinical decision-making.

PURPOSE To assess the efficacy of pembrolizumab in patients with advanced relapsed or refractory mycosis fungoides (MF) or Sézary syndrome (SS). PATIENTS AND METHODS CITN-10 is a single-arm, multicenter phase II trial of 24 patients with advanced MF or SS. Patients were treated with pembrolizumab 2 mg/kg every 3 weeks for up to 24 months. The primary end point was overall response rate by consensus global response criteria. RESULTS Patients had advanced-stage disease (23 of 24 with stage IIB to IV MF/SS) and were heavily pretreated with a median of four prior systemic therapies. The overall response rate was 38% with two complete responses and seven partial responses. Of the nine responding patients, six had 90% or more improvement in skin disease by modified Severity Weighted Assessment Tool, and eight had ongoing responses at last follow-up. The median duration of response was not reached, with a median response follow-up time of 58 weeks. Immune-related adverse events led to treatment discontinuation in four patients. A transient worsening of erythroderma and pruritus occurred in 53% of patients with SS. This cutaneous flare reaction did not result in treatment discontinuation for any patient. The flare reaction correlated with high PD-1 expression on Sézary cells but did not associate with subsequent clinical responses or lack of response. Treatment responses did not correlate with expression of PD-L1, total mutation burden, or an interferon-γ gene expression signature. CONCLUSION Pembrolizumab demonstrated significant antitumor activity with durable responses and a favorable safety profile in patients with advanced MF/SS.

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.

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). 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. We found that -not-amplified (-NA) tumors had significantly higher cytotoxic TIL signatures compared with -amplified (-A) tumors. A reported MYCN activation signature was significantly associated with poor outcome for high-risk patients with -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. 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 -NA tumors. Our findings suggest that these high-risk patients with-NA neuroblastoma may benefit from additional immunotherapies incorporated into the current therapeutic strategies. .

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 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 promoter and releases the transcriptional stalling of 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.

Purpose We undertook a multi-dimensional 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 Food and Drug Administration approved drugs. Results Fifty-nine patients in 20 diagnostic categories were enrolled from 2010 to 2014. Ages ranged from 7-months-old 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.