Desmoplastic small round cell tumor (DSRCT) is characterized by the EWSR1–WT1 t(11;22) (p13:q12) translocation. Few additional putative drivers have been identified, and research has suffered from a lack of model systems. Next-generation sequencing (NGS) data from 68 matched tumor-normal samples, whole-genome sequencing data from 10 samples, transcriptomic and affymetrix array data, and a bank of DSRCT patient-derived xenograft (PDX) are presented. EWSR1–WT1 fusions were noted to be simple, balanced events. Recurrent mutations were uncommon, but were noted in TERT (3%), ARID1A (6%), HRAS (5%), and TP53 (3%), and recurrent loss of heterozygosity (LOH) at 11p, 11q, and 16q was identified in 18%, 22%, and 34% of samples, respectively. Comparison of tumor-normal matched versus unmatched analysis suggests overcalling of somatic mutations in prior publications of DSRCT NGS data. Alterations in fibroblast growth factor receptor 4 (FGFR4) were identified in 5 of 68 (7%) of tumor samples, whereas differential overexpression of FGFR4 was confirmed orthogonally using 2 platforms. PDX models harbored the pathognomic EWSR1–WT1 fusion and were highly representative of corresponding tumors. Our analyses confirm DSRCT as a genomically quiet cancer defined by the balanced translocation, t(11;22)(p13:q12), characterized by a paucity of secondary mutations but a significant number of copy number alterations. Against this genomically quiet background, recurrent activating alterations of FGFR4 stood out, and suggest that this receptor tyrosine kinase, also noted to be highly expressed in DSRCT, should be further investigated. Future studies of DSRCT biology and preclinical therapeutic strategies should benefit from the PDX models characterized in this study.

Implications:

These data describe the general quiescence of the desmoplastic small round cell tumor (DSRCT) genome, present the first available bank of DSRCT model systems, and nominate FGFR4 as a key receptor tyrosine kinase in DSRCT, based on high expression, recurrent amplification, and recurrent activating mutations.

This article is featured in Highlights of This Issue, p. 1097

Desmoplastic small round cell tumor (DSRCT) is a rare and aggressive small round blue cell malignancy, which occurs predominantly in adolescent and young adult males, and is characterized by abdominopelvic sarcomatosis exhibiting multilineage cellular nests of epithelial, muscular, mesenchymal, and neural differentiation admixed with desmoplastic stroma (1). The exact incidence of DSRCT is poorly defined, although because its histopathologic description at least 850 cases have been reported in the literature (2). Overall survival at 5 years from retrospective series ranges from 11% (3, 4) to 28% (5, 6). Presenting signs and symptoms are often gastrointestinal in nature, and can include abdominal pain and bloating, change in stool pattern, weight loss or weight gain secondary to ascites, and nausea. The cancer is thought to originate along the peritoneal surface and is almost universally metastatic at the time of presentation with common sites including liver, spleen, and nodal involvement above the diaphragm (6), although disease can arise in disparate sites including the testes (7–9) and the central nervous system (10–12). Several staging systems have been proposed (13), with the most recent using imaging characteristics to define intermediate (no liver involvement or ascites), high-risk (either liver involvement or ascites), and very high-risk disease (both liver involvement and ascites) (6).

Given the rarity of this disease, no randomized trials addressing treatment have been performed, and nearly all the available literature describes anecdotal or retrospective experiences. As such, no standard of care therapy has emerged, although treatment almost always includes high-dose alkylator-based chemotherapy (14), and an attempt at complete cytoreductive surgery. The inclusion of whole abdomen/pelvis radiation therapy using an intensity-modulated approach is advocated by some centers, (3, 15, 16) but not by others (5). Additional consolidative local control approaches include investigational treatment with radiolabeled antibody (17), as well as hyperthermic intraperitoneal chemotherapy (HIPEC; refs. 18–21). Recurrence despite high-intensity, multi-modality therapy is common (22–24).

DSRCT is molecularly characterized by the pathognomonic EWSR1–WT1 t(11;22) (p13:q12) translocation (25–27). Beyond this disease defining translocation, the understanding of recurrent oncogenic alterations or genomically defined subsets of DSRCT has remained limited (28). Moreover, there is a marked paucity of established and characterized preclinical DSRCT models. An early study based on hotspot mutation genotyping of 29 common oncogenes (not including FGFR4) in 24 DSRCT samples did not find any mutations (29). An additional study described that nearly one third of 135 mutated genes detected from whole exome sequencing of six DSRCT samples were related to the DNA damage response (DDR) or epithelial–mesenchymal transition (EMT) pathways, suggesting their biologic importance for this disease (30). Another series compared 35 DSRCT samples with 88 Ewing sarcoma samples using a combination of next-generation sequencing (NGS), IHC, chromogenic in situ hybridization or fluorescence in situ hybridization and detected higher expression levels of androgen receptor (AR), TUBB3, EGFR, and TOPO2A (31). The most recently published series describes 83 DSRCT samples analyzed by the FoundationOne Heme assay, and identified recurrent alterations in FGFR4, ARID1A, TP53, MSH3, and MLL3 (32).

Herein, we describe a multimodal approach to define the molecular landscape of DSRCT using hybridization capture-based NGS assay MSK-IMPACT (33) of 68 matched tumor-normal samples from 53 patients, a targeted RNA sequencing assay MSK-Solid Fusion (34) from 17 patients, and whole-genome and RNA sequencing of 10 samples from 9 patients. Furthermore, we describe MSK-IMPACT data from 24 established DSRCT patient-derived xenograft (PDX) models from 10 patients.

Patients and samples

This project was approved by the Institutional Review Board of Memorial Sloan Kettering Cancer Center (MSKCC) and was conducted in accordance with the U.S. Common Rule. Formalin-fixed paraffin-embedded DSRCT samples from a subset of patients treated at MSKCC between 2014 and 2019 were submitted for clinical sequencing using the MSK-IMPACT panel (35).

Targeted hybridization capture DNA sequencing

The MSK-IMPACT assay interrogates all exons and select introns of 468 genes to identify mutations, copy number changes, microsatellite instability status, and select structural variants in its current iteration (previous versions interrogated 341 and 410 genes), and was conducted as described previously (36). DNA from both formalin-fixed paraffin-embedded or fresh tissue biopsies as well as a matched normal blood sample was extracted and sheared. Barcoded libraries from tumor and normal samples were captured, sequenced, and subjected to a custom pipeline to identify somatic mutations. Analysis consisted of deep sequencing of all exons and selected introns of a cancer-associated gene panel. Allele-specific copy number analysis was performed with FACETS, utilizing the additional information from a tumor/matched normal design. A segment mean ≥0.2 or ≤0.2 was used to define low-level gains and losses, respectively. In an unmatched analysis of our DSRCT cohort, the matched normals were removed and the pipeline was run using a pooled normal control to find suspected somatic mutations.

Whole-genome sequencing

Tumor was sequenced to approximately 80X with a matched normal sample at 60X. Single-base substitutions were called using CaVEMan (Cancer Variants through Expectation Maximization; http://cancerit.github.io/CaVEMan/; ref. 37), Strelka2 (https://github.com/Illumina/strelka), and MuTect2 (https://github.com/broadinstitute/gatk). Copy number and cellularity information for CaVEMan were predicted with the Battenberg algorithm (37) using 1,000 Genomes (38) loci within the NGS data. Small somatic insertions and indels were identified using a modified version of Pindel (https://github.com/cancerit/cgpPindel; ref. 39), Strekla2 (https://github.com/Illumina/strelka), and MuTect2 (https://github.com/broadinstitute/gatk). Structural rearrangements were detected by BRASS (Breakpoints via Assembly; https://github.com/cancerit/BRASS), GRIDSS (https://github.com/PapenfussLab/gridss), and SvABA (https://github.com/walaj/svaba). To improve specificity, a number of custom postprocessing filters were applied. Mutational signature analysis of the substitutions was performed using the R package MutationalPatterns (https://bioconductor.org/packages/release/bioc/html/MutationalPatterns.html). Small insertion/deletions were interrogated for the presence of either short tandem repeat or microhomology at the breakpoints, as described previously (40) Data relating to this patient will be deposited in the European Genome-phenome Archive (www.ebi.ac.uk/).

Anchored multiplex RNA sequencing

Tumor specimens were analyzed using the MSK-Solid Fusion assay, a targeted RNA-based panel that utilizes the Archer Anchored Multiplex PCR (AMPTM) technology and NGS to detect genes fusions in solid tumor samples (41). Unidirectional gene-specific primers (GSP) were used to target specific exons in 62 genes known to be involved in chromosomal rearrangements. Final targeted amplicons were sequenced on an Illumina MiSeq instrument and data were analyzed using the ArcherTM Software (V4.0.10).

Affymetrix microarray analyses

Raw expression data from Filion and colleagues on 137 sarcomas samples was obtained (42). Probe set data were normalized and expression estimates were derived using robust multi-array average (RMA) as previously described (43). Specifically, FGFR4 expression in DSRCT was compared with four other fusion-positive sarcomas.

Cell viability assay

JN-DSRCT1 cells were plated in 384-well plates with a density of 500 cells/well, and were maintained in DMEM/F12 media with 10% FBS, 1% penicillin/streptomycin, and fungizone at 37°C. The cells were tested for mycoplasma as described previously (44) approximately once every 4 months (last performed October 2020). The alamar blue assay was utilized as described previously (45) with reagent added at hour 72, and results recorded at hour 96. A negative control with 1% DMSO, and a positive control with 1% killer mix in 1% DMSO were utilized. IC50 values were calculated using Prism GraphPad.

Patient demographics

The clinical characteristics of the 53 patients with DSRCT are summarized in Table 1. The median age at diagnosis was 21, with a range of 7 to 47 years. There was a strong male predominance as has been noted in prior studies with 48 (91%) of patients being male, and 5 (9%) being female.

Targeted NGS/anchored multiplex PCR

A total of 68 matched tumor/normal samples from 53 patients were analyzed by MSK-IMPACT, a hybridization capture-based NGS assay interrogating all exons and select introns of 468 genes. Almost all samples (63, 92%) represented intra-abdominal disease, whereas five samples (8%) were obtained from extra-abdominal sites. Fifty-four (79%) samples were obtained following treatment of any type, and 14 (21%) prior to any intervention. Somatic mutations as detected by MSK-IMPACT in the 68 samples from 53 patients are shown in Fig. 1A and Supplementary Table S1. Recurrent events were uncommon but included alterations in the TERT promoter (3%), ARID1A (6%), HRAS (4%), FGFR4 (7%; further described below), and TP53 (3%). Alterations in HRAS included one amplification, one deletion, and one G12V mutation known to be a hotspot. The average tumor mutational burden (TMB) was 0.87 mutations per megabase (range 0–3.9), and none of the samples showed microsatellite instability (MSI) based on analysis of the MSK-IMPACT sequencing reads using the MSIsensor algorithm (46).

Although the MSK-IMPACT hybrid capture-based NGS assay routinely detects the EWSR1–WT1 fusion of DSRCT by tiling the EWSR1 introns involved in the genomic rearrangement, a subset of these specimens (n = 17) was also analyzed via an anchored multiplex RT-PCR assay for the detection of gene fusions, the MSK-Solid Fusion assay (47). The EWSR1–WT1 fusion was orthogonally confirmed using the MSK-Solid Fusion assay in all 17 samples tested, and detected breakpoints between exons 7 of EWSR1 and exon 8 of WT1 in 12 samples, exons 9 of EWSR1 and exon 8 of WT1 in 3 samples, and exons 10 of EWSR1 and exon 8 of WT1 in 2 samples. This variability in the exon structure of the EWSR1–WT1 fusion transcript is well described (48).

Normal matched versus unmatched analysis

A recent study reported frequent mutations in multiple genes including MSH3 (14%) and MLL3 (16%) in a cohort of tumors profiled using the FoundationOne Medicine Heme assay (32). Our results using MSK-IMPACT identify a much lower overall mutation rate, with a single sample harboring an MLL3 somatic mutation (R199fs) and none with MSH3 somatic alterations. A major difference between these two studies is the use of a matched germline normal control for each sample profiled in our analysis. Therefore, we compared our MSK-IMPACT profiling results using a matched germline control versus an unmatched analysis where presumptive germline mutations were identified and removed on the basis of population databases (Fig. 1B). The mean TMB was significantly higher in the unmatched analysis (14.34 vs. 0.85; P value <0.01). Furthermore, 12 (18%) KMT2C/MLL3 and 6 (8%) MSH3 variants were identified in the unmatched analysis (whereas only one MLL3 mutation and no MSH3 mutations were found by our standard matched tumor:normal analysis). These data suggest that previous unmatched analyses may have overestimated the somatic mutation rate in DSRCT.

Copy number analysis

Copy number alterations of 65 matched tumor/normal samples were analyzed using the Fraction and Allele-Specific Copy Number Estimates and Tumor Sequencing (FACETS) algorithm as described previously (49). Three samples were excluded due to a lack of matched normal material. Recurrent loss of heterozygosity (LOH) in chromosome arms 11p, 11q, and 16q were identified in 18%, 22%, and 34% of 50 patients, respectively (Supplementary Fig. S1A). Hierarchical clustering of observed gains and losses reveals that these alterations separate into non-overlapping clusters (Supplementary Fig. S1B).

Whole genome and transcriptome sequencing

Paired DNA/RNA whole genome and transcriptome sequencing was performed on 10 DSRCT samples from 9 patients, all but one of whom were previously treated. Mutational burden in these samples was found to range from 0.32 to 3.05 mutations/Mb, with a median of 0.80 mutations/Mb. A circos plot depicting whole-genome sequencing of a representative DSRCT sample is shown in Fig. 2A–D, and the remainder are shown in Supplementary Fig. S2A. A majority of the analyzed DSRCT genomes are relatively simple, with low mutational burden, minor involvement of arm or chromosome level aneuploidies, and few structural rearrangements as exemplified in Fig. 2A–D. However, this simplicity was not universal as demonstrated in additional circos plots in Supplementary Fig. S2A.

Given recent evidence that some canonical gene fusions in Ewing sarcoma are generated from chromoplexy (50, 51), structural variant analysis was performed to determine if similar events were occurring in these DSRCT samples. This analysis revealed that the EWSR1–WT1 fusions in this subset were predominantly simple, balanced translocations as shown in the representative plot in Fig. 2E, and the remainder as shown in Supplementary Fig. S2B.

Small deletions were also observed scattered across the genome, but none were recurrent or thought to involve known tumor drivers (Supplementary Fig. S3). Although the samples overall had relatively few structural events aside from the canonical EWSR1–WT1 translocation, most samples had additional nonclustered, nonrecurrent translocations (Supplementary Fig. S3).

Clonality analysis from WGS data revealed the presence of subclones in all samples analyzed. Signature analyses comparing intratumoral subclones did not reveal significant differences between primary clones and subclones. Furthermore, no suspected secondary tumorigenic drivers of disease were noted in these subclones (data not shown).

In the 10 DSRCT samples analyzed by whole-genome sequencing, mutational signatures 1, five (associated with aging) and eight (associated with double stranded break repair) were found (52), but could not be related to presentation, treatment, or outcome.

FGFR4 analysis

Alterations in FGFR4 were identified in 5 of 67 (7%) of tumor samples analyzed by MSK-IMPACT. Amplification of FGFR4 was found in two samples, and a deep deletion in one sample. Notably, the two samples with amplification (from different patients) exhibited remarkably similar copy number profiles as shown in Fig. 3A. In addition, an exon 13 p.V550 L point mutation, previously identified as oncogenic (53), was identified in two samples (Fig. 1A).

To investigate expression levels of FGFR4 in DSRCT, whole transcriptome RNA sequencing results from a cohort of 154 pediatric malignancies, including 8 DSRCT samples, were analyzed (Fig. 3B). A clear trend toward increased transcripts per million (TPM) was noted amongst the DSRCT samples, as the 4 samples with the highest FGFR4 TPMs in the entire cohort were DSRCT. None of 8 DSRCT samples in this cohort had FGFR point mutations, amplification, or deletion. Additionally, microarray expression data from 137 tumor samples across 5 different fusion-associated sarcoma subtypes were analyzed and found to demonstrate very strong differential expression of FGFR4, with a mean of 4.9X higher expression in DSRCT when compared with all other sarcoma types in the analysis (Fig. 3C; Supplementary Table S1).

To explore the potential susceptibility of FGFR4 overexpression in DSRCT to targeted inhibition, the antiproliferative effect of a panel of FGFR4 inhibitors was assessed against the only widely available DSRCT cell line, JN-DSRCT (54). RNA sequencing revealed high expression of FGFR4 in the JN-DSRCT cell line which was similar to DSRCT patient samples, and was higher than that found in Ewing sarcoma patient samples (Supplementary Fig. S4A). Notably JN-DSRCT does not harbor FGFR4 mutation or amplification. In addition, protein expression of FGFR4 was confirmed by Western immunoblotting (Supplementary Fig. S4B). Hepatocellular carcinoma cell lines with established sensitivity and resistance (Hep 3B and H2122, respectively) and the Ewing sarcoma cell line A673 which was expected to be resistant were used as controls. The drug panel consisted of four tyrosine kinase inhibitors with varying activity against FGFR4 including ponatinib, rogaritinib, fisogatinib, and BLU9931. Only ponatinib demonstrated modest activity against the JN-DSRCT cell line with an IC50 of 0.46 μmol/L. None of the other inhibitors tested demonstrated activity against JN-DSRCT with all IC50 values >25 μmol/L (Supplementary Figs. S4C and S4D). These data suggest that FGFR4 overexpression in the absence of amplification or activating mutation is not sufficient to render DSRCT cells responsive to FGFR4 inhibition.

DSRCT PDXs

To further characterize the molecular landscape of DSRCT, as well as to produce a set of validated preclinical, in vivo models, a total of 120 freshly extracted DSRCT surgical samples from 39 patients were implanted beginning in September 2017. At the time of this publication, 44 of these samples had engrafted, 10 were pending engraftment or failure, and 66 had failed engraftment, for an overall engraftment rate of 36.6%. Tumor material from a subset of the engrafted models in their first or second passage consisting of 24 patient-derived DSRCT xenografts from 10 separate patients was analyzed using the same 468 gene MSK-IMPACT assay as described above. These data were then compared with MSK-IMPACT data from matched patient tumor samples (Fig. 4).

The pathognomonic EWSR1–WT1 fusion was observed by MSK-IMPACT in all PDX and patient samples. Two patient tumor samples harbored additional oncogenic mutations in HRAS and ARID1A, both of which were recapitulated in the corresponding PDX models. In addition, copy number alterations noted in tumor were nearly all represented in the PDX models. Several mutations were identified in PDX samples (TERT, TP53, BRAF), which were not present in matched patient tumor. MSK-IMPACT sequencing reads were reviewed and the lack of these aberrations in the original patient sample was confirmed.

The primary tumorigenic driver of DSRCT, the EWSR1–WT1 fusion, has been well described (25, 26, 48, 55), but few recent series have described any additional recurrent molecular findings (28, 30–32). As is the case for most malignancies driven by oncogenic transcription factors, knowledge of the tumorigenic fusion has yet to inform rational therapeutic strategies. Mutations in FGFR4 (7%), ARID1A (6%), and TP53 (3%) were identified in the MSK-IMPACT cohort described herein. The somatic alteration rate in this study is much lower than that described in a recent study by Chow and colleagues evaluating a cohort of DSRCT samples profiled using the FoundationOne Heme targeted NGS assay (32). Recurrent mutations in MSH3 and MLL3 were not detected in our cohort when the analysis was performed with a matched germline control but were seen at a significantly higher frequency when an unmatched analysis was performed. The use of public SNP databases for filtering of variant calls from unmatched, tumor-only sequencing is imperfect, as it does not eliminate so-called “private” SNPs that are not present in those databases (36). Therefore, we suspect that some genomic alterations reported in previous studies may represent germline SNPs, rather than tumor-specific somatic mutations.

Despite the quiescent genome suggested by this paucity of recurrent mutations, copy number alterations were noted in a striking number of DSRCT samples. This finding aligns with similar observations in other pediatric solid tumor malignancies including Ewing sarcoma (56), rhabdomyosarcoma (57), and osteosarcoma (58, 59). The clinical significance of these copy number-related subgroups is unknown, and may be elucidated through prospective studies with universal molecular profiling and standardized treatment regimens.

Whole-genome sequencing of DSRCT samples did not reveal additional information regarding secondary driver oncogenic events beyond the EWSR1–WT1 fusion. However, it did identify these driver fusions as simple and balanced translocations (Fig. 2B), distinguishing them from subsets of those found in Ewing sarcoma which have been shown to be generated from chromoplexy (50). This finding is somewhat unsurprising given the same orientation of EWSR1 and WT1 on their respective chromosomes, unlike EWSR1 and ERG, for example, which have opposite orientations on their respective chromosomes necessitating events like chromoplexy for fusion formation (51).

In this report, activating alterations in FGFR4 were noted in a subset of tumors (6%; Fig. 1A). Differential expression of FGFR4 was demonstrated by transcriptome analysis and Affymetrix array-based profiling (Fig. 3B and C; Supplementary Table S2), and did not correlate with FGFR point mutations, amplifications, or deletions. These data are aligned with a recently published series which also identified FGFR4 as a potential therapeutic target in DSRCT (60), as well as the experience in rhabdomyosarcoma in which high FGFR4 expression has not been found to be genetically encoded by FGFR4 amplification (53). The FGFRs consist of four highly conserved transmembrane receptor tyrosine kinases (FGFR1–4), and have received significant attention as anticancer targets (53, 61, 62). Specifically, FGFR4 has been of particular interest in pediatric sarcomas, including both embryonal and alveolar subtypes of rhabdomyosarcoma (53, 63). Furthermore, the EWSR1–WT1 fusion oncoprotein has been shown to contribute to FGFR4 transcriptional upregulation (64), suggesting a mechanism for the differential expression described above. However, in vitro testing of multiple FGFR inhibitors demonstrated limited activity in cell proliferation assays in the established JN-DSRCT cell line. Our data suggest that FGFR4 overexpression alone is insufficient to render sensitivity to biologic inhibition of FGFR4. Rare cases of DSRCT with FGFR4 mutations which have been established as activating in other malignancies (e.g., V550 L point mutation) may however demonstrate susceptibility to FGFR inhibition. Indeed, our finding of two patients with DSRCT with FGFR4 V550L, together with a previous report of FGFR4 V510 L in 3 patients with DSRCT (65), establish FGFR4 as the only receptor tyrosine kinase recurrently activated by somatic mutations in DSRCT, and present in 3% to 4% of cases. Additional investigation of the utility of FGFR4 inhibition in DSRCT in the setting of activating FGFR4 point mutations is warranted and would be facilitated by the generation of isogenic pair-matched cell lines. Although small molecule inhibition of the FGFR4 kinase is unlikely to provide significant antitumor effect in the majority of patients with DSRCT, we demonstrate high expression of FGFR4 across the majority of DSRCT samples in our cohort. This finding provides a rationale for investigation of FGFR4 as a therapeutic target for antibody–drug conjugates or novel immunotherapeutics including chimeric antigen receptor (CAR) T-cell or bi-specific T-cell engaging products.

Some preclinical vulnerabilities have been identified in DSRCT cells without discrete genomic data to explain them. For example, DSRCT PDX models have been shown to be exquisitely sensitive to CHK1 kinase inhibition without clear genomic rationale (66). Although one publication analyzing seven DSRCT samples suggested genomic instability due to functional pathway analysis which converged on the DDR pathway (30), a subsequent analysis in a larger cohort did not confirm this (32), nor was this suggested in the current analysis. However, the strength of the preclinical data has supported the initiation of an ongoing clinical trial evaluating CHK1 kinase inhibition in combination with irinotecan in patients with DSRCT (NCT04095221). These events highlight the gaps in knowledge regarding biologic activity, which are critically needed to dovetail with ongoing molecular discoveries.

Although therapies directed against the underlying fusion oncogene in DSRCT have proved elusive thus far, drug development against this and other similar central drivers may still be possible. Direct shRNA approaches (NCT02736565) and direct inhibition of the E26 transformation-specific (ETS) family of oncoproteins (NCT02657005) are being evaluated in Ewing sarcoma and provide guideposts around which to try to develop additional innovative clinical approaches for DSRCT. Any such effort will significantly benefit from the rich PDX bank described and molecularly characterized here, in which these approaches can be rigorously evaluated.

W.D. Tap reports personal fees from Eli Lilly, EMD Serono, Eisai, Munipharma, C4 Therapeutics, Daiichi Sankyo, Blueprint, GlaxoSmithKline, Agios, NanoCarrier, and Deciphera during the conduct of the study; and has a patent for Companion Diagnostic for CDK4 inhibitors - 14/854,329 pending to MSKCC/SKI and Companion Diagnostic for CDK4 inhibitors - 14/854,329 pending to MSKCC/SKI; scientific advisory board and stock ownership at Certis Oncology Solutions, co-founder and stock ownership at Atropos Therapeutics, and scientific advisory board at Innova Therapeutics. G. Gundem reports other support from ISABL Technologies outside the submitted work. L.H. Wexler reports personal fees from EUSA Pharma outside the submitted work. P.A. Meyers reports other support from Amgen, personal fees from Astellas, Lilly, and Salarius outside the submitted work. E. Pappaemanuil reports other support from Isabl Inc. during the conduct of the study; and is a founder and equity holder of Isabl Inc. No disclosures were reported by the other authors.

E.K. Slotkin: Conceptualization, resources, data curation, supervision, writing–original draft, writing–review and editing. A.S. Bowman: Data curation, investigation, writing–review and editing. M.F. Levine: Data curation, investigation, writing–review and editing. F. Dela Cruz: Resources. D.F. Coutinho: Data curation, investigation. G.I. Sanchez: Investigation. N. Rosales: Investigation. S. Modak: Resources, supervision. W.D. Tap: Resources. M.M. Gounder: Resources. K.A. Thornton: Resources. N. Bouvier: Project administration. D. You: Project administration. G. Gundem: Supervision. J.T. Gerstle: Resources. T.E. Heaton: Resources. M.P. LaQuaglia: Resources. L.H. Wexler: Resources. P.A. Meyers: Resources. A.L. Kung: Conceptualization, resources, data curation, formal analysis, funding acquisition, investigation, writing–review and editing. E. Papaemmanuil: Formal analysis, supervision, writing–review and editing. A. Zehir: Data curation, formal analysis, supervision, investigation, writing–review and editing. M. Ladanyi: Resources, data curation, supervision, investigation, writing–review and editing. N. Shukla: Resources, supervision, investigation, writing–review and editing.

Funding for this work was provided by the Scarlett Fund, the Paulie Strong Foundation, the Montag Family Initiative for Pediatric Cancer Research, the Grayson Fund, the Willens Family Fund, and the MSK Cancer Center Support Grant (No. P30 CA008748).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Gerald
WL
,
Miller
HK
,
Battifora
H
,
Miettinen
M
,
Silva
EG
,
Rosai
J
. 
Intra-abdominal desmoplastic small round-cell tumor. Report of 19 cases of a distinctive type of high-grade polyphenotypic malignancy affecting young individuals
.
Am J Surg Pathol
1991
;
15
:
499
513
.
2.
Mora
J
,
Modak
S
,
Cheung
N-K
,
Meyers
P
,
de Alava
E
,
Kushner
B
, et al
Desmoplastic small round cell tumor 20 years after its discovery
.
Future Oncol
2015
;
11
:
1071
81
.
3.
Lal
DR
,
Su
WT
,
Wolden
SL
,
Loh
KC
,
Modak
S
,
La Quaglia
MP
. 
Results of multimodal treatment for desmoplastic small round cell tumors
.
J Pediatr Surg
2005
;
40
:
251
5
.
4.
Forlenza
CJ
,
Kushner
BH
,
Kernan
N
,
Boulad
F
,
Magnan
H
,
Wexler
L
, et al
Myeloablative chemotherapy with autologous stem cell transplant for desmoplastic small round cell tumor
.
Sarcoma
2015
;
2015
:
269197
.
5.
Subbiah
V
,
Lamhamedi-Cherradi
SE
,
Cuglievan
B
,
Menegaz
BA
,
Camacho
P
,
Huh
W
, et al
Multimodality treatment of desmoplastic small round cell tumor: chemotherapy and complete cytoreductive surgery improve patient survival
.
Clin Cancer Res
2018
;
24
:
4865
73
.
6.
Saltsman
JA
,
Price
AP
,
Goldman
DA
,
Hammond
WJ
,
Danzer
E
,
Magnan
H
, et al
A novel image-based system for risk stratification in patients with desmoplastic small round cell tumor
.
J Pediatr Surg
2018
;
55
:
376
80
.
7.
Zhang
GM
,
Zhu
Y
,
Gan
HL
,
Ye
DW
. 
Testicular desmoplastic small round cell tumor: a case report and review of literature
.
World J Surg Oncol
2014
;
12
:
227
.
8.
Manjula
MV
,
Pawar
YS
. 
Primary desmoplastic small round cell tumor of the testis: First case in India and review of the literature
.
J Cancer Res Ther
2015
;
11
:
650
.
9.
Sedig
L
,
Geiger
J
,
Mody
R
,
Jasty-Rao
R
. 
Paratesticular desmoplastic small round cell tumors: a case report and review of the literature
.
Pediatr Blood Cancer
2017
;
64
.
10.
Neder
L
,
Scheithauer
BW
,
Turel
KE
,
Arnesen
MA
,
Ketterling
RP
,
Jin
L
, et al
Desmoplastic small round cell tumor of the central nervous system: report of two cases and review of the literature
.
Virchows Arch
2009
;
454
:
431
9
.
11.
Thondam
SK
,
Plessis
DDu
,
Cuthbertson
DJ
,
Das
KSV
,
Javadpour
M
,
MacFarlane
IA
, et al
Intracranial desmoplastic small round cell tumor presenting as a suprasellar mass
.
J Neurosurg
2015
;
122
:
773
7
.
12.
Umeda
K
,
Saida
S
,
Yamaguchi
H
,
Okamoto
S
,
Okamoto
T
,
Kato
I
, et al
Central nervous system recurrence of desmoplastic small round cell tumor following aggressive multimodal therapy: A case report
.
Oncol Lett
2016
;
11
:
856
60
.
13.
Hayes-Jordan
A
,
LaQuaglia
MP
,
Modak
S
. 
Management of desmoplastic small round cell tumor
.
Semin Pediatr Surg
2016
;
25
:
299
304
.
14.
Kushner
BH
,
LaQuaglia
MP
,
Wollner
N
,
Meyers
PA
,
Lindsley
KL
,
Ghavimi
F
, et al
Desmoplastic small round-cell tumor: prolonged progression-free survival with aggressive multimodality therapy
.
J Clin Oncol
1996
;
14
:
1526
31
.
15.
Desai
NB
,
Stein
NF
,
LaQuaglia
MP
,
Alektiar
KM
,
Kushner
BH
,
Modak
S
, et al
Reduced toxicity with intensity modulated radiation therapy (IMRT) for desmoplastic small round cell tumor (DSRCT): an update on the whole abdominopelvic radiation therapy (WAP-RT) experience
.
Int J Radiat Oncol Biol Phys
2013
;
85
:
e67
72
.
16.
Goodman
KA
,
Wolden
SL
,
La Quaglia
MP
,
Kushner
BH
. 
Whole abdominopelvic radiotherapy for desmoplastic small round-cell tumor
.
Int J Radiat Oncol Biol Phys
2002
;
54
:
170
6
.
17.
Modak
S
,
Carrasquillo
J
,
LaQuaglia
M
, et al
Intraperitoneal radioimmunotherapy for desmoplastic small round cell tumor: results of a phase I study (NCT01099644)
.
Cancer Res
2018
;
78
:
CT006
.
18.
Hayes-Jordan
AA
,
Coakley
BA
,
Green
HL
,
Xiao
L
,
Fournier
KF
,
Herzog
CE
, et al
Desmoplastic small round cell tumor treated with cytoreductive surgery and hyperthermic intraperitoneal chemotherapy: results of a Phase 2 trial
.
Ann Surg Oncol
2018
;
25
:
872
7
.
19.
Hayes-Jordan
A
,
Green
HL
,
Lin
H
,
Owusu-Agyemang
P
,
Fitzgerald
N
,
Arunkumar
R
, et al
Complete cytoreduction and HIPEC improves survival in desmoplastic small round cell tumor
.
Ann Surg Oncol
2014
;
21
:
220
4
.
20.
Hayes-Jordan
A
,
Green
H
,
Fitzgerald
N
,
Xiao
L
,
Anderson
P
. 
Novel treatment for desmoplastic small round cell tumor: hyperthermic intraperitoneal perfusion
.
J Pediatr Surg
2010
;
45
:
1000
6
.
21.
Hayes-Jordan
A
,
Anderson
P
,
Curley
S
,
Herzog
C
,
Lally
KP
,
Green
HL
, et al
Continuous hyperthermic peritoneal perfusion for desmoplastic small round cell tumor
.
J Pediatr Surg
2007
;
42
:
E29
32
.
22.
Lae
ME
,
Roche
PC
,
Jin
L
,
Lloyd
RV
,
Nascimento
AG
. 
Desmoplastic small round cell tumor: a clinicopathologic, immunohistochemical, and molecular study of 32 tumors
.
Am J Surg Pathol
2002
;
26
:
823
35
.
23.
Quaglia
MP
,
Brennan
MF
. 
The clinical approach to desmoplastic small round cell tumor
.
Surg Oncol
2000
;
9
:
77
81
.
24.
Schwarz
RE
,
Gerald
WL
,
Kushner
BH
,
Coit
DG
,
Brennan
MF
,
La Quaglia
MP
. 
Desmoplastic small round cell tumors: prognostic indicators and results of surgical management
.
Ann Surg Oncol
1998
;
5
:
416
22
.
25.
Gerald
WL
,
Rosai
J
,
Ladanyi
M
. 
Characterization of the genomic breakpoint and chimeric transcripts in the EWS-WT1 gene fusion of desmoplastic small round cell tumor
.
Proc Natl Acad Sci U S A
1995
;
92
:
1028
32
.
26.
Gerald
WL
,
Ladanyi
M
,
de Alava
E
,
Cuatrecasas
M
,
Kushner
BH
,
LaQuaglia
MP
, et al
Clinical, pathologic, and molecular spectrum of tumors associated with t(11;22)(p13;q12): desmoplastic small round-cell tumor and its variants
.
J Clin Oncol
1998
;
16
:
3028
36
.
27.
Ladanyi
M
,
Gerald
W
. 
Fusion of the EWS and WT1 genes in the desmoplastic small round cell tumor
.
Cancer Res
1994
;
54
:
2837
40
.
28.
Jiang
Y
,
Subbiah
V
,
Janku
F
,
Ludwig
JA
,
Naing
A
,
Benjamin
RS
, et al
Novel secondary somatic mutations in Ewing's sarcoma and desmoplastic small round cell tumors
.
PLoS One
2014
;
9
:
e93676
.
29.
Quayle
SN
,
Chheda
MG
,
Shukla
SA
,
Wiedemeyer
R
,
Tamayo
P
,
Dewan
RW
, et al
Integrative functional genomics identifies RINT1 as a novel GBM oncogene
.
Neuro Oncol
2012
;
14
:
1325
31
.
30.
Devecchi
A
,
De Cecco
L
,
Dugo
M
,
Penso
D
,
Dagrada
G
,
Brich
S
, et al
The genomics of desmoplastic small round cell tumor reveals the deregulation of genes related to DNA damage response, epithelial-mesenchymal transition, and immune response
.
Cancer Commun
2018
;
38
:
70
.
31.
Bulbul
A
,
Shen
JP
,
Xiu
J
,
Tamayo
P
,
Husain
H
. 
Genomic and proteomic alterations in desmoplastic small round blue-cell tumors
.
JCO Precis Oncol
2018
;
2
:
PO.17.00170.1
9
.
32.
Chow
WA
,
Yee
JK
,
Tsark
W
,
Wu
X
,
Qin
H
,
Guan
M
, et al
Recurrent secondary genomic alterations in desmoplastic small round cell tumors
.
BMC Med Genet
2020
;
21
:
101
.
33.
Cheng
DT
,
Prasad
M
,
Chekaluk
Y
,
Benayed
R
,
Sadowska
J
,
Zehir
A
, et al
Comprehensive detection of germline variants by MSK-IMPACT, a clinical diagnostic platform for solid tumor molecular oncology and concurrent cancer predisposition testing
.
BMC Med Genomics
2017
;
10
:
33
.
34.
Zhu
G
,
Benayed
R
,
Ho
C
,
Mullaney
K
,
Sukhadia
P
,
Rios
K
, et al
Diagnosis of known sarcoma fusions and novel fusion partners by targeted RNA sequencing with identification of a recurrent ACTB-FOSB fusion in pseudomyogenic hemangioendothelioma
.
Mod Pathol
2019
;
32
:
609
20
.
35.
Zehir
A
,
Benayed
R
,
Shah
RH
,
Syed
A
,
Middha
S
,
Kim
HR
, et al
Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients
.
Nat Med
2017
;
23
:
703
13
.
36.
Cheng
DT
,
Mitchell
TN
,
Zehir
A
,
Shah
RH
,
Benayed
R
,
Syed
A
, et al
Memorial Sloan kettering-integrated mutation profiling of actionable cancer targets (MSK-IMPACT): A hybridization capture-based next-generation sequencing clinical assay for solid tumor molecular oncology
.
J Mol Diagn
2015
;
17
:
251
64
.
37.
Nik-Zainal
S
,
Van Loo
P
,
Wedge
DC
,
Alexandrov
LB
,
Greenman
CD
,
Lau
KW
, et al
The life history of 21 breast cancers
.
Cell
2012
;
149
:
994
1007
.
38.
The 1000 Genomes Project Consortium
. 
An integrated map of genetic variation from 1,092 human genomes
.
Nature
2012
;
491
:
56
65
.
39.
Raine
KM
,
Hinton
J
,
Butler
AP
,
Teague
JW
,
Davies
H
,
Tarpey
P
, et al
cgpPindel: identifying somatically acquired insertion and deletion events from paired end sequencing
.
Curr Protoc Bioinformatics
2015
;
52
:
1571
12
.
40.
Nik-Zainal
S
,
Davies
H
,
Staaf
J
,
Ramakrishna
M
,
Glodzik
D
,
Zou
X
, et al
Landscape of somatic mutations in 560 breast cancer whole-genome sequences
.
Nature
2016
;
534
:
47
54
.
41.
Benayed
R
,
Offin
M
,
Mullaney
K
,
Sukhadia
P
,
Rios
K
,
Desmeules
P
, et al
High yield of RNA sequencing for targetable kinase fusions in lung adenocarcinomas with no mitogenic driver alteration detected by DNA Sequencing and low tumor mutation burden
.
Clin Cancer Res
2019
;
25
:
4712
22
.
42.
Filion
C
,
Motoi
T
,
Olshen
Ab
,
Laé
M
,
Emnett
Rj
,
Gutmann
Dh
, et al
The EWSR1/NR4A3 fusion protein of extraskeletal myxoid chondrosarcoma activates the PPARG nuclear receptor gene
.
J Pathol
2009
;
217
:
83
93
.
43.
Irizarry
RA
,
Hobbs
B
,
Collin
F
,
Beazer-Barclay
YD
,
Antonellis
KJ
,
Scherf
U
, et al
Exploration, normalization, and summaries of high density oligonucleotide array probe level data
.
Biostatistics
2003
;
4
:
249
64
.
44.
Mariotti
E
,
Mirabelli
P
,
Di Noto
R
,
Fortunato
G
,
Salvatore
F
. 
Rapid detection of mycoplasma in continuous cell lines using a selective biochemical test
.
Leuk Res
2008
;
32
:
323
6
.
45.
Kumar
P
,
Nagarajan
A
,
Uchil
PD
. 
Analysis of cell viability by the alamarblue assay
.
Cold Spring Harb Protoc
2018
;
2018
.
46.
Middha
S
,
Zhang
L
,
Nafa
K
,
Jayakumaran
G
,
Wong
D
,
Kim
HR
, et al
Reliable pan-cancer microsatellite instability assessment by using targeted next-generation sequencing data
.
JCO Precis Oncol
2017
;
2017
:
PO.17.00084
.
47.
Zheng
Z
,
Liebers
M
,
Zhelyazkova
B
,
Cao
Yi
,
Panditi
D
,
Lynch
KD
, et al
Anchored multiplex PCR for targeted next-generation sequencing
.
Nat Med
2014
;
20
:
1479
84
.
48.
Antonescu
CR
,
Gerald
WL
,
Magid
MS
,
Ladanyi
M
. 
Molecular variants of the EWS-WT1 gene fusion in desmoplastic small round cell tumor
.
Diagn Mol Pathol
1998
;
7
:
24
8
.
49.
Shen
R
,
Seshan
VE
. 
FACETS: allele-specific copy number and clonal heterogeneity analysis tool for high-throughput DNA sequencing
.
Nucleic Acids Res
2016
;
44
:
e131
.
50.
Anderson
ND
,
de Borja
R
,
Young
MD
,
Fuligni
F
,
Rosic
A
,
Roberts
ND
, et al
Rearrangement bursts generate canonical gene fusions in bone and soft tissue tumors
.
Science
2018
;
361
:
eaam8419
.
51.
Imielinski
M
,
Ladanyi
M
. 
Fusion oncogenes-genetic musical chairs
.
Science
2018
;
361
:
848
.
52.
Alexandrov
LB
,
Kim
J
,
Haradhvala
NJ
,
Huang
MiNi
,
Tian Ng
AW
,
Wu
Y
, et al
The repertoire of mutational signatures in human cancer
.
Nature
2020
;
578
:
94
101
.
53.
JGt
T
,
Cheuk
AT
,
Tsang
PS
,
Chung
J-Y
,
Song
YK
,
Desai
K
, et al
Identification of FGFR4-activating mutations in human rhabdomyosarcomas that promote metastasis in xenotransplanted models
.
J Clin Invest
2009
;
119
:
3395
407
.
54.
Nishio
J
,
Iwasaki
H
,
Ishiguro
M
,
Ohjimi
Y
,
Fujita
C
,
Yanai
F
, et al
Establishment and characterization of a novel human desmoplastic small round cell tumor cell line, JN-DSRCT-1
.
Lab Invest
2002
;
82
:
1175
82
.
55.
Gedminas
JM
,
Chasse
MH
,
McBrairty
M
,
Beddows
I
,
Kitchen-Goosen
SM
,
Grohar
PJ
. 
Desmoplastic small round cell tumor is dependent on the EWS-WT1 transcription factor
.
Oncogenesis
2020
;
9
:
41
.
56.
Lynn
M
,
Wang
Y
,
Slater
J
,
Shah
N
,
Conroy
J
,
Ennis
S
, et al
High-resolution genome-wide copy-number analyses identify localized copy-number alterations in Ewing sarcoma
.
Diagn Mol Pathol
2013
;
22
:
76
84
.
57.
Shern
JF
,
Chen
Li
,
Chmielecki
J
,
Wei
JS
,
Patidar
R
,
Rosenberg
M
, et al
Comprehensive genomic analysis of rhabdomyosarcoma reveals a landscape of alterations affecting a common genetic axis in fusion-positive and fusion-negative tumors
.
Cancer Discov
2014
;
4
:
216
31
.
58.
Sayles
LC
,
Breese
MR
,
Koehne
AL
,
Leung
SG
,
Lee
AG
,
Liu
H-Y
, et al
Genome-informed targeted therapy for osteosarcoma
.
Cancer Discov
2019
;
9
:
46
63
.
59.
Perry
JA
,
Kiezun
A
,
Tonzi
P
,
Van Allen
EM
,
Carter
SL
,
Baca
SC
, et al
Complementary genomic approaches highlight the PI3K/mTOR pathway as a common vulnerability in osteosarcoma
.
Proc Natl Acad Sci U S A
2014
;
111
:
E5564
73
.
60.
Hingorani
P
,
Dinu
V
,
Zhang
X
,
Lei
H
,
Shern
JF
,
Park
J
, et al
Transcriptome analysis of desmoplastic small round cell tumors identifies actionable therapeutic targets: a report from the Children's
Oncology Group. Sci Rep
2020
;
10
:
12318
.
61.
Babina
IS
,
Turner
NC
. 
Advances and challenges in targeting FGFR signalling in cancer
.
Nat Rev Cancer
2017
;
17
:
318
32
.
62.
Cao
L
,
Yu
Y
,
Bilke
S
,
Walker
RL
,
Mayeenuddin
LH
,
Azorsa
DO
, et al
Genome-wide identification of PAX3-FKHR binding sites in rhabdomyosarcoma reveals candidate target genes important for development and cancer
.
Cancer Res
2010
;
70
:
6497
508
.
63.
Crose
LES
,
Etheridge
KT
,
Chen
C
,
Belyea
B
,
Talbot
LJ
,
Bentley
RC
, et al
FGFR4 blockade exerts distinct antitumorigenic effects in human embryonal versus alveolar rhabdomyosarcoma
.
Clin Cancer Res
2012
;
18
:
3780
90
.
64.
Saito T
YM
,
Motoi
T
,
Iwasaki
H
,
Nagao
T
,
Ladanyi
M
,
Yao
T
. 
EWS-WT1 chimeric protein in desmoplastic small round cell tumor is a potent transactivator of FGFR4
.
Cancer Sci Ther
2012
;
4
:
335
40
.
65.
Calmette
A
. 
The protection of mankind against tuberculosis: being an address before the medico-chirurgical society of edinburgh
.
Edinb Med J
1922
;
29
:
93
104
.
66.
Lowery
CD
,
Dowless
M
,
Renschler
M
,
Blosser
W
,
VanWye
AB
,
Stephens
JR
, et al
Broad Spectrum activity of the checkpoint kinase 1 inhibitor prexasertib as a single agent or chemopotentiator across a range of preclinical pediatric tumor models
.
Clin Cancer Res
2019
;
25
:
2278
89
.