Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Journal Sections
    • Subscriptions
    • Reviewing
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Collections
      • COVID-19 & Cancer Resource Center
      • Clinical Trials
      • Immuno-oncology
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
    • Journal Press Releases
  • COVID-19
  • Webinars
  • 10th Anniversary
  • Search More

    Advanced Search

  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Cancer Discovery
Cancer Discovery
  • Home
  • About
    • The Journal
    • AACR Journals
    • Journal Sections
    • Subscriptions
    • Reviewing
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Collections
      • COVID-19 & Cancer Resource Center
      • Clinical Trials
      • Immuno-oncology
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
    • Journal Press Releases
  • COVID-19
  • Webinars
  • 10th Anniversary
  • Search More

    Advanced Search

Research Articles

Genomic Landscape of Ewing Sarcoma Defines an Aggressive Subtype with Co-Association of STAG2 and TP53 Mutations

Franck Tirode, Didier Surdez, Xiaotu Ma, Matthew Parker, Marie Cécile Le Deley, Armita Bahrami, Zhaojie Zhang, Eve Lapouble, Sandrine Grossetête-Lalami, Michael Rusch, Stéphanie Reynaud, Thomas Rio-Frio, Erin Hedlund, Gang Wu, Xiang Chen, Gaelle Pierron, Odile Oberlin, Sakina Zaidi, Gordon Lemmon, Pankaj Gupta, Bhavin Vadodaria, John Easton, Marta Gut, Li Ding, Elaine R. Mardis, Richard K. Wilson, Sheila Shurtleff, Valérie Laurence, Jean Michon, Perrine Marec-Bérard, Ivo Gut, James Downing, Michael Dyer, Jinghui Zhang and Olivier Delattre
Franck Tirode
1INSERM U830, Laboratory of Genetics and Cancer Biology, Institut Curie, Paris, France.
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Didier Surdez
1INSERM U830, Laboratory of Genetics and Cancer Biology, Institut Curie, Paris, France.
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiaotu Ma
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthew Parker
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marie Cécile Le Deley
4Departement d'Epidémiologie et de Biostatistiques, Gustave Roussy, Villejuif, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Armita Bahrami
5Department of Pathology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhaojie Zhang
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eve Lapouble
6Unité de Génétique Somatique, Centre Hospitalier, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sandrine Grossetête-Lalami
1INSERM U830, Laboratory of Genetics and Cancer Biology, Institut Curie, Paris, France.
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Rusch
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stéphanie Reynaud
6Unité de Génétique Somatique, Centre Hospitalier, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas Rio-Frio
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Erin Hedlund
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gang Wu
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiang Chen
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gaelle Pierron
6Unité de Génétique Somatique, Centre Hospitalier, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Odile Oberlin
7Departement de Pédiatrie, Gustave Roussy, Villejuif, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sakina Zaidi
1INSERM U830, Laboratory of Genetics and Cancer Biology, Institut Curie, Paris, France.
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gordon Lemmon
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pankaj Gupta
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bhavin Vadodaria
8The Pediatric Cancer Genome Laboratory, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Easton
8The Pediatric Cancer Genome Laboratory, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marta Gut
9Centro Nacional de Análisis Genómico (CNAG), Barcelona, Spain.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Li Ding
10Department of Genetics, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
11Department of Medicine, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
12Siteman Cancer Center, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elaine R. Mardis
10Department of Genetics, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
11Department of Medicine, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
12Siteman Cancer Center, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard K. Wilson
10Department of Genetics, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
11Department of Medicine, The Genome Institute, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
12Siteman Cancer Center, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sheila Shurtleff
5Department of Pathology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Valérie Laurence
13Département d'Oncologie Medicale, Adolescents et Jeunes Adultes, Centre Hospitalier, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jean Michon
14Département d'Oncologie Pediatrique, Adolescents et Jeunes Adultes, Centre Hospitalier, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Perrine Marec-Bérard
15Institute for Paediatric Haematology and Oncology, Leon Bérard Cancer Centre, University of Lyon, Lyon, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ivo Gut
9Centro Nacional de Análisis Genómico (CNAG), Barcelona, Spain.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
James Downing
8The Pediatric Cancer Genome Laboratory, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Dyer
16Department of Developmental Neurobiology, St. Jude Children's Research Hospital, Memphis, Tennessee.
17Howard Hughes Medical Institute, Chevy Chase, Maryland.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jinghui Zhang
3Department of Computational Biology, St. Jude Children's Research Hospital, Memphis, Tennessee.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: olivier.delattre@curie.fr jinghui.zhang@stjude.org
Olivier Delattre
1INSERM U830, Laboratory of Genetics and Cancer Biology, Institut Curie, Paris, France.
2Centre de Recherche, Institut Curie, Paris, France.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: olivier.delattre@curie.fr jinghui.zhang@stjude.org
DOI: 10.1158/2159-8290.CD-14-0622 Published November 2014
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Ewing sarcoma is a primary bone tumor initiated by EWSR1–ETS gene fusions. To identify secondary genetic lesions that contribute to tumor progression, we performed whole-genome sequencing of 112 Ewing sarcoma samples and matched germline DNA. Overall, Ewing sarcoma tumors had relatively few single-nucleotide variants, indels, structural variants, and copy-number alterations. Apart from whole chromosome arm copy-number changes, the most common somatic mutations were detected in STAG2 (17%), CDKN2A (12%), TP53 (7%), EZH2, BCOR, and ZMYM3 (2.7% each). Strikingly, STAG2 mutations and CDKN2A deletions were mutually exclusive, as confirmed in Ewing sarcoma cell lines. In an expanded cohort of 299 patients with clinical data, we discovered that STAG2 and TP53 mutations are often concurrent and are associated with poor outcome. Finally, we detected subclonal STAG2 mutations in diagnostic tumors and expansion of STAG2-immunonegative cells in relapsed tumors as compared with matched diagnostic samples.

Significance: Whole-genome sequencing reveals that the somatic mutation rate in Ewing sarcoma is low. Tumors that harbor STAG2 and TP53 mutations have a particularly dismal prognosis with current treatments and require alternative therapies. Novel drugs that target epigenetic regulators may constitute viable therapeutic strategies in a subset of patients with mutations in chromatin modifiers. Cancer Discov; 4(11); 1342–53. ©2014 AACR.

This article is highlighted in the In This Issue feature, p. 1243

Introduction

Ewing sarcoma is the second most common primary malignant bone tumor in children and adolescents; the mean age at diagnosis is 15 years. Ewing sarcoma can affect any bone, but the most common primary sites are pelvis, femur, and tibia (reviewed in ref. 1). The annual incidence of Ewing sarcoma is approximately 3 per million, with a slight male bias. Histologically, Ewing sarcoma belongs to the group of small blue round-cell tumors, and the tumor cells often have abundant cytoplasmic glycogen and express CD99 on the plasma membrane (1).

Genetically, most Ewing sarcomas are characterized by a specific t(11;22)(q12;q11.2) translocation that fuses the EWSR1 gene on chromosome (chr) 22 with the FLI1 gene on chr 11 (2). In 10% to 15% of cases, EWSR1 is fused to genes encoding other members of the ETS family of transcription factors, including ERG, ETV1, E1AF, or FEV (1). In an even smaller number of cases, TAF15 and TLS/FUS, which encode the two other members of the TET family of RNA-binding proteins, may be fused to ETS family members (1). All fusions juxtapose the N-terminal domain of the TET gene family member to the DNA-binding domain of the ETS gene family member. These TET–ETS fusions are potent oncogenes that can transform NIH3T3 cells (3) by perturbing the expression of genes required for a variety of cellular processes, including cell-cycle regulation, signal transduction, and telomere maintenance (reviewed in ref. 1).

Chromosome or array-based comparative genomic hybridization (CGH) as well as SNP arrays have identified recurrent DNA copy-number alterations (CNA) in Ewing sarcoma (4–10). The most common copy-number gains occur in whole chromosomes 8 and 12 and the q (long) arm of chr 1. The long arm of chr 16 and the CDKN2A locus on chr 9p are the most common copy-number losses in Ewing sarcoma. The adverse prognosis conferred by chr 1q gain and chr 16q or CDKN2A loss has been reported, as has the negative impact of TP53 mutations (11). Finally, somatic STAG2 mutations were recently observed in a significant fraction of Ewing sarcoma cases (21%; ref. 12). Although the role of TET–ETS oncogenes in Ewing sarcoma tumorigenesis and progression has been extensively studied and the relation of copy-number changes to prognosis is emerging, relatively little is known about additional secondary genetic lesions in Ewing sarcoma beyond these chromosomal lesions.

To identify secondary genetic lesions that contribute to Ewing sarcoma tumorigenesis after formation of the TET–ETS fusion, we performed whole-genome sequencing (WGS) of 112 tumors and their matched germline DNA. The most frequent point mutations involved the STAG2 and TP53 genes, and the prognostic significance of these mutations was further demonstrated in a series of 299 cases. STAG2 mutations were significantly associated with the occurrence of structural variations and were mutually exclusive with CDKN2A deletions. In some cases, we also observed a small number of STAG2-deficient tumor cells that survived treatment and comprised the major clone in the recurrent tumors.

Results

Ewing Sarcoma Has Low Numbers of Single-Nucleotide and Structural Variants

Our discovery set for WGS comprised 112 Ewing sarcomas with matched germline DNA (Fig. 1 and Supplementary Table S1). All Ewing sarcoma tumors, with the exception of one (SJ001301) that had insufficient tumor sample for analysis, expressed EWSR1–ETS fusions: EWSR1–FLI1 in 101 cases, EWSR1–ERG in 9 cases, and EWSR1–ETV1 in one case. Tumor and germline DNA were sequenced at a median depth of 35× and 25×, respectively. Mapping, detection, and annotation of single-nucleotide variants (SNV), insertions or deletions (indels), and structural variants (SV), functional predictions, and CNAs were computed from the WGS data as previously described (13–15). Eighty percent of the tumors had >70% tumor purity, leading to a 98% power for detecting mutations present in the predominant tumor clones in this cohort (Supplementary Table S2A).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

A comprehensive profile of the genetic abnormalities in Ewing sarcoma and associated clinical information. Key clinical characteristics are indicated, including primary site, type of tissue, and metastatic status at diagnosis, follow-up, and last news. Below is the consistency of detection of gene fusions by RT-PCR and whole-genome sequencing (WGS). The numbers of structural variants (SV) and single-nucleotide variants (SNV) as well as indels are reported in grayscale. The presence of the main copy-number changes, chr 1q gain, chr 16 loss, chr 8 gain, chr 12 gain, and interstitial CDKN2A deletion is indicated. Listed last are the most significant mutations and their types. See Supplementary Table S2 for the complete lists of SNVs/indels, SVs, and CNAs. For gene mutations, “others” refers to: duplication of exon 22 leading to frameshift (STAG2), deletion of exon 2 to 11 (BCOR), and deletion of exons 1 to 6 (ZMYM3).

The median number of somatic SVs was 7 (range, 0–66) per tumor (Supplementary Table S2B). In most cases (106/112; 95%), WGS detected SVs within the previously described EWSR1 and ETS chromosome breakpoint regions (ref. 16; Fig. 1; Supplementary Table S2C; and Supplementary Fig. S1A–S1D). Five cases (SJ001303, SJ001320, IC198, IC273, and IC086) exhibited chromothripsis, including three cases with chromothripsis on chr 21 and 22 associated with EWSR1–ERG fusions (SJ001303, IC198, and IC273) and one case involving chr 22 associated with an EWSR1–FLI1 fusion (SJ001320). CNAs could be reliably analyzed from WGS data in 103 cases. Nine cases were excluded from CNA analysis due to low tumor purity or uneven sequencing coverage. The most frequent CNAs were gain of whole chr 8 (49/103; 47%), gain of whole chr 12 (22/103; 21%), gain of the long arm of chr 1 (19/103; 18%), deletion of the long arm of chr 16 (18/103; 17%), and deletion of the CDKN2A locus on the short arm of chr 9 (12/103; 12%; Fig. 1; Supplementary Fig. S2; and Supplementary Table S2D). Chr 1q gain and chr 16q loss were correlated with shorter survival (P = 2 × 10−5 and P = 0.0037, respectively, log-rank test; Fig. 2A). As chr 16q deletion and chr 1q gain were highly significantly co-associated (P = 10−8, Fisher exact test), their combination did not show an additional effect on overall survival (Fig. 2A). Chr 8 and chr 12 gains also showed a significant, although less pronounced, co-association (P = 1.63 × 10−3, Fisher exact test), but neither chr 8 nor chr 12 gains, nor their combination, were correlated with shorter survival (data not shown).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Prognostic significance of CNAs, SVs, and SNV/indels. Kaplan–Meier overall survival estimates according to (A) chr 1q gain, chr 16q loss; chr 1q gain and chr 16q loss and (B) number of SNVs/indels. Samples were stratified according to the number of genomic SNVs/indels and split into tertiles; C, a large number of SVs. The overall survival of patients whose tumors harbor an outlier number of SVs (boxplot distribution shown on the left) is compared with that of other patients. Patients with a fractured genome, low tumor purity, or death by causes other than Ewing sarcoma were excluded from the analysis.

Experimental validation by custom capture and Illumina sequencing of all WGS-predicted SNVs and indels in 19 cases showed a 95.6% verification rate. Across the entire cohort, the median number of SNVs was 319 (range, 13–1,747) per genome (Supplementary Table S2E and S2F). The background mutation rate ranged from 8.0 × 10−9 to 1.4 × 10−6 (median, 2.4 × 10−7) per base. The predominant changes were C(G) > T(A) transitions (Supplementary Fig. S3). No bias in SNV distribution was observed in cases with the highest number of SNVs. In particular, rainfall plots for 31 samples that had at least 400 somatic SNVs across the genome showed no patterns of kataegis (17). On average, there were 10 (range, 1–39) coding variants per tumor, and the ratio of missense to silent mutations was 2.4. We observed a positive correlation between age at diagnosis and the number of SNVs (r2 = 0.42; P = 2.7 × 10−5, Pearson correlation). Patients older than 20 years at diagnosis had significantly more SNVs than did younger patients (P = 0.001, Mann–Whitney U test). Survival analysis showed a negative correlation with the number (tertile-based) of SNVs/indels, i.e., a greater number of SNVs/indels was associated with shorter survival time (Fig. 2B; P = 0.04, log-rank test). Tertile-based survival did not differ significantly according to the number of SVs. However, the 8 patients whose tumors had a large number of SVs (outliers in the box plot distribution of SVs shown in Fig. 2C) had very poor outcomes (Fig. 2C; P = 0.003, log-rank test).

The Most Frequent Coding Variants Occur in STAG2, TP53, and Epigenetic Regulators

The gene most frequently carrying a somatic mutation in our cohort was STAG2 (17%, 19/112). We identified 6 nonsense mutations, 10 indels leading to frameshifts, 1 missense mutation, 1 splice-site mutation, and 1 duplication of exon 22 (Figs. 1 and 3A). As the STAG2 protein is an integral member of the cohesin complex (18) and was found to be associated with aneuploidy (12), we also investigated the relation of STAG2 mutations to the number of SVs across the discovery cohort. A significantly greater number of SVs was observed in STAG2-mutated cases (Fig. 3B; P = 0.006, Mann–Whitney U test). In contrast, STAG2 status was not associated with the number of SNVs or indels (Fig. 3C).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

STAG2 mutations and their prognostic significance in Ewing sarcomas. A, schematic of the STAG2 protein and mutations. Mutations found in tumor samples are indicated above the protein, and those observed in cell lines are indicated below. Mutation nomenclature is based on the NM_001042749 reference sequence. Exon and amino-acid numbering is indicated below the protein. The recurrent R216* mutation was observed in 7 cases. One tumor (case IC871) had two mutations (indicated in bold). SCD, stromalin conservative domain; GR, glutamine-rich region. Box plots show comparison of the number of SVs (B) and SNVs/indels (C) in wild-type (WT) and STAG2-mutated tumor samples. Samples with a fractured genome or low tumor cell content (see Fig. 1) were excluded from analysis, leaving 17 STAG2-mutated cases and 86 wild-type cases. Box represents the central 50% of data points (interquartile range). Upper and lower whiskers represent the largest and smallest observed values within 1.5 times the interquartile range from the ends of the box. Circles, individual values. P values were determined by using the Mann–Whitney U test. D, overall survival among 299 patients according to STAG2 mutation status. The number of patients in the different groups is indicated in brackets. E, overall survival of the 299 patients according to their STAG2 and/or TP53 mutation status.

TP53 was mutated in 8 cases (Fig. 1). All mutations were missense, with the exception of one nonsense mutation (p.R317* according to NM_000546), and were described in the COSMIC database. After excluding the very large genes that are recurrently mutated in most cancer genome studies (TTN, CSMD1, MACF1, and RYR2; ref. 19), the third most frequently mutated genes were EZH2, BCOR, and ZMYM3, which each presented with 3 mutations (3/112, 2.7%; Fig. 1 and Supplementary Table S2). All three EZH2 mutations were missense mutations within the SET domain (Y646F, Y646H, and A682G according to NM_004456). BCOR exhibited one missense mutation (S1083I, according to NM_017745), one indel leading to a frameshift (M1259fs), and one 116-kb intragenic deletion (Fig. 1 and Supplementary Table S2). ZMYM3 exhibited two indels (L82fs according to NM_201599) and one 17-kb intragenic deletion (Fig. 1 and Supplementary Table S2).

All other somatic gene mutations were observed in less than three cases. Mutations affecting epigenetic regulators have been found to be significantly associated with some pediatric cancers (20). In addition to the mutations in EZH2, BCOR, and ZMYM3, we identified novel somatic mutations in SETD2, MLL2, MLL3, and PRDM9 (Fig. 1 and Supplementary Table S2). Of note, two novel missense mutations were observed in EWSR1. Finally, we used the significantly mutated gene (SMG) test in the mutational significance in cancer (MuSiC) suite (21) to identify genes that are significantly enriched in somatic SNVs and indels. Only STAG2, TP53, and EZH2 were found to be significantly enriched (Supplementary Table S3).

STAG2 and CDKN2A Genetic Lesions Are Mutually Exclusive

When investigating the relationships between gene mutations, SVs, and CNAs, we found a mutually exclusive pattern of STAG2 and CDKN2A genetic alterations (Fig. 1). To confirm this mutually exclusive profile, we investigated STAG2 and CDKN2A in a panel of 19 Ewing sarcoma cell lines. STAG2 mutations and CDKN2A deletions were observed in 9 and 6 of the 19 cell lines, respectively (Table 1). The exclusive pattern of STAG2 and CDKN2A alterations shown in primary tumors (Fig. 1) was fully replicated in the cell lines (Table 1). Across the 15 cell lines that could be investigated by Western blot, all cases with STAG2 mutations but one (MHH-ES-1) expressed p16. Reciprocally, all cases with CDKN2A deletion expressed STAG2 (Supplementary Fig. S4A and S4B). When tumor and cell line results were combined, this mutually exclusive pattern of alteration was highly significant (P = 0.0079, Fisher exact test). The frequency of TP53 mutations was extremely high in the cell lines (Table 1). Altogether, all tested cell lines harbored at least one STAG2, TP53, or CDKN2A lesion.

View this table:
  • View inline
  • View popup
Table 1.

Genomic status of STAG2, CDKN2A, and TP53 in Ewing sarcoma cell lines

STAG2 and TP53 Mutations Are Co-Associated in Highly Aggressive Tumors

To determine whether STAG2 and/or TP53 mutations are associated with outcome in Ewing sarcoma, we analyzed these genes by targeted capture sequencing in an additional 199 French patients with Ewing sarcoma. Across the whole series, 30% of patients had metastatic spread at diagnosis. The presence of a metastasis was associated with a shorter overall survival time (P = 4 × 10−4, log-rank test). In total, 41 patients (13.2%) had STAG2 mutations (Fig. 3A and Supplementary Table S4) and 16 patients (5.2%) had TP53 mutations. The STAG2 mutations included 15 nonsense, 4 missense, 17 frameshift, and 4 splice-site mutations, 2 in-frame deletions, and 1 exon duplication (Figs. 1 and 3A and Supplementary Table S4). One tumor (IC871) had two distinct STAG2 mutations. Overall survival data were available for 299 patients. The presence of a STAG2 mutation was not significantly associated with dismal prognostic factors, including tumor size, response to chemotherapy, resection quality, or tumor spread. However, patients with STAG2 mutations demonstrated a significantly lower probability of survival, similar to patients with TP53-mutated tumors (Fig. 3D and E). Patients with neither STAG2 nor TP53 mutations had the highest probability of survival, and patients whose tumors carried mutations in both genes had the worst outcome (Fig. 3E). A significant decrease in overall survival of patients with either STAG2 or TP53 mutation alone was not observed. In our cohort, STAG2 and TP53 mutations were significantly co-associated (P = 2.4 × 10−4, Fisher exact test).

We also explored the CDKN2A status across these additional tumors. Expanding the CDKN2A cohort confirmed the exclusion pattern with STAG2 mutations. Indeed, we identified only 2 tumors with both STAG2 mutations and CDKN2A deletions. When compiling all our data (299 tumors and 19 cell lines), the overlap between STAG2 and CDKN2A genetic lesions was much lower than expected by chance [Fisher test: 0.0076, STAG2/CDKN2A, wild-type (WT)/WT: 221, WT/Mut: 49, Mut/WT: 46, Mut/Mut: 2]. CDKN2A status was not significantly associated with overall survival across the whole series (Supplementary Fig. S5A–S5C).

Subclonal STAG2 Mutations May Expand at Relapse

Finally, we investigated whether STAG2 mutation occurs in subclones within some tumors and whether it evolves during the course of the disease. We first took advantage of the high coverage obtained in the capture-based sequencing experiments to investigate the ratios of mutated/wild-type alleles. Seven diagnostic samples showed evidence of subclonal mutations, i.e., a mutant allele frequency <0.25 despite high tumor purity (Supplementary Table S4; example in Fig. 4A). In 21 cases, STAG2 immunostaining could be investigated in paired primary/relapse or pre-/post-therapy samples. In 18 cases, STAG2 immunostaining in paired primary/relapse or pre-/post-therapy was unaltered, of which 16 were positive and two were negative at both time points. However, in three cases, STAG2 staining at relapse revealed a reduction in STAG2-immunopositive cells (Fig. 4B). Consistent with the immunostaining result, loss-of-function STAG2 mutations were detected at relapse with high allelic fractions but were either not detected (SJEWS001303) or detected at a subclonal level at diagnosis (SJEWS014721; Supplementary Table S5).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Subclonal presence of STAG2 mutations. A, Integrative Genomics Viewer representation showing the subclonal presence of STAG2 mutations in one sample. B, evolution of STAG2 staining between diagnosis and relapse in two independent cases. Whereas only a small subset of tumor cells lacked STAG2 expression at diagnosis (see insets), the tumor cells were homogeneously negative at relapse. The few STAG2-positive stromal cells serve as an internal positive control.

Discussion

To our knowledge, the work reported here is the most comprehensive genomic analysis of Ewing sarcoma performed to date. The cases we studied met all of the criteria defining bona fide Ewing sarcoma, including clinical, pathologic, and molecular findings. The background mutation rate of Ewing sarcoma was relatively low (2.4 × 10−7), with a median of 10 coding somatic mutations per tumor. The Ewing sarcoma mutation rate is much lower than that usually observed in adult cancers and in the upper range of what is described in other pediatric solid malignancies and brain tumors, including neuroblastoma (22–24), retinoblastoma (14), rhabdomyosarcoma (25, 26), medulloblastoma (27, 28), pilocytic astrocytoma (29), pediatric glioblastoma (30), and osteosarcoma (31). We also observed a positive correlation between age at diagnosis and the number of SNVs.

The CNA most frequently detected in the present study was gain of chr 8, which was observed in close to 50% of cases, in agreement with previous CGH or SNP-array data (4, 5, 7). Loss of chr 16q and gain of chr 1q were strongly co-associated, which is fully consistent with the presence in these tumors of a derivative chr 16 resulting from an unbalanced t(1;16) translocation previously identified by cytogenetics in Ewing sarcoma (32). However, no SVs specific for this translocation were detected, consistent with the hypothesis that this t(1;16) translocation occurs within repeated elements of centromeric regions that cannot be reliably detected by WGS. In this cohort, we also report that 16q and/or 1q gains have strong negative prognostic significance.

The most frequently mutated gene in Ewing sarcoma is STAG2. STAG1 and STAG2, the human orthologs of yeast Scc3p, encode components of the cohesin multiprotein complex that plays an essential role in sister chromatin cohesion (18). STAG1 and STAG2 exist in different cohesin complexes that are essential for telomere or centromere cohesion, respectively (18). STAG2 mutations were initially observed in a diverse range of cancers, including glioblastoma, melanoma, and Ewing sarcoma (12). Subsequently, STAG2 mutations were described in a significant proportion of bladder cancers (12, 33, 34) and myeloid neoplasms (35). Although experimental systems have shown that STAG1 and STAG2 inactivation drives aneuploidy (12, 18), STAG2 mutations were not found to be associated with aneuploidy or CNAs in bladder cancer (33, 36). The case may be slightly different in Ewing sarcoma, as we observed a positive correlation between the presence of STAG2 mutation and the number of SVs. However, the interpretation of this correlation must take into account the strong co-association of STAG2 and TP53 mutations in our cohort. When cases with only one of these two mutations are considered, the positive correlation between STAG2 mutation and the number of SVs is no longer significant. The analysis of survival data must also take into account the association between STAG2 and TP53 mutations. Indeed, in our extended series of patients, the prognostic significance of STAG2 mutation appears to be strongly dependent on the coexistence of a TP53 mutation. The prognosis of cases with both STAG2 and TP53 mutations appears particularly unfavorable (Fig. 2E). Together, these data suggest that STAG2 and TP53 mutation may cooperate to increase genetic instability in a particularly aggressive subtype of Ewing sarcoma. Consistent with this hypothesis, it is noteworthy that STAG2 and TP53 mutations are much more frequent in cell lines derived mainly from aggressive cases. Finally, our results suggest that STAG2-mutated Ewing sarcoma subclones at diagnosis may evolve and become the major clone at recurrence. Further investigation of the relation of clonal expansion to tumor progression or response to therapy will be of great interest.

We observed a previously unreported, mutually exclusive pattern of STAG2 and CDKN2A mutations in Ewing sarcoma. This mutual exclusivity was observed in primary tumors and confirmed in cell lines. In addition to their role in sister chromatin cohesion, STAG2-containing cohesin complexes play an essential role in nuclear chromatin organization, particularly in the epigenetic mechanisms of insulation through direct interaction between STAG2 and CTCF, a multifunctional transcription factor that regulates chromosomal boundaries of gene expression, as recently demonstrated at the H19/Igf2 locus (37). Interestingly, CTCF has also been shown to regulate the CDKN2A locus (38), raising the possibility that STAG2 loss of function alters the epigenetic regulation of CDKN2A in CDKN2A–wild-type cases. However, as previously reported (6, 39–41), methylation is not a common mechanism for CDKN2A inactivation in Ewing sarcoma and is therefore not expected to occur in most STAG2–wild-type cases. The role of STAG2 in chromatin structure, particularly in the distribution of histone marks, and expression of the CDKN2A locus should be further investigated in depth.

Three EZH2 mutations (Y646F, Y646H, and A682G, all in the SET domain) were observed in our cohort of patients. EZH2 encodes a member of the multiprotein polycomb repressive complex 2 (PRC2), which catalyzes trimethylation of histone H3 lysine 27 (H3K27me3). Residues Y646 and to a lesser extent A682 are frequently mutated in B-cell lymphoma, and these mutations have been shown to enhance EZH2 enzymatic activity and promote malignant lymphoid transformation (42–44). Mutations of EZH2 have also been observed in a subset of acute T-cell and myeloid malignancies (15, 45). In addition to EZH2, potentially deleterious mutations in ZMYM3 and BCOR, which also encode epigenetic regulators, were reported in three cases each. In total, we observed recurrent mutations in epigenetic regulators in 17 of 112 Ewing sarcoma cases (15.2%). As described above, recent data strongly suggest that STAG2 plays a major role in epigenetic insulation and may therefore be considered an epigenetic regulator. This finding reinforces the need for studies that clarify how mutations affecting the epigenetic landscape of Ewing sarcoma may cooperate with the EWSR1–ETS fusion to promote the development of overt Ewing sarcoma.

After submission of this article, Brohl and colleagues (46) published an article describing the genomic landscape of Ewing sarcoma based mostly on exome sequencing and RNA sequencing. The observed frequency of STAG2, TP53, and CDKN2A is similar to the findings reported in this article. They also observed the association of TP53 and STAG2 mutations. However, significant correlation with clinical outcome could not be demonstrated, possibly due to the smaller size of the patient cohort. Finally, the exclusive pattern of CDKN2A and STAG2 alterations was not reported in the Brohl and colleagues dataset. The different techniques used in the two reports and the different sizes of the patient series may account for this discrepancy, which requires further investigation.

In conclusion, our comprehensive genetic analysis of Ewing sarcoma identified recurrent mutations in STAG2, TP53, and epigenetic regulators. We showed that a STAG2 mutation gains prognostic significance when associated with TP53 mutations and that a STAG2-mutated subclone may expand during the course of the disease. Finally, the mutual exclusion between STAG2 and CDKN2A loss-of-function mutations suggests that these alterations may be, at least partially, redundant.

Methods

Patients and Tumors

Our discovery cohort comprised 112 patients with Ewing sarcoma; both tumor and germline samples underwent WGS. All tumors selected for WGS were predicted to contain a large proportion of tumor cells based on pathology reports, previous CGH or SNP arrays, and/or a low Ct (cycle threshold) of the EWSR1–ETS fusion assessed by qRT-PCR. Eighteen Ewing sarcomas were obtained from the St. Jude tissue resource core facility for genome sequence analysis with St. Jude Institutional Review Board (IRB) approval for the Pediatric Cancer Genome Project. The remaining cases were those referred to Institut Curie from all over France for molecular diagnosis of Ewing sarcoma. Samples were stored in a tumor bank at the Institut Curie. The genetic study was approved by the IRB of the Institut Curie and by the Comité de Protection de Personnes Ile-de-France I (regional ethics committee; GenEwing no. IC 2009-02); specific informed consent was provided. Most patients were treated according to the EuroEwing protocol (47). An anonymization procedure was performed before compilation of clinical, histologic, and biologic information in a secure database with restricted access. All tumors included in this study were positive for the EWSR1–ETS fusion. Detailed clinicopathologic and sequencing information is provided in Fig. 1 and Supplementary Table S1.

The follow-up set comprising 199 tumor DNAs from EWSR1–ETS-positive Ewing sarcomas was distinct from the discovery set and consisted of patients treated according to the EuroEwing99 protocol.

Cell Lines

Ewing sarcoma cell lines were obtained from various sources: A673, RD-ES, SK-ES-1, and SK-NM-C from the ATCC; MHH-ES1 and TC-71 from the German Collection of Microorganisms and Cell Cultures (DSMZ); EW-1, EW-3, EW-7, EW-16, and EW-18 from the International Agency for Research on Cancer; STA-ET-1, STA-ET-3, and STA-ET-8 from Prof. Heinrich Kovar, Children's Cancer Research Institute–Childhood (Vienna, Austria); and EW-22, EW-23, MIC, ORS, and POE from the Institut Curie. Cell lines were authenticated by their TP53 genotype, which included mutations previously described.

DNA and RNA Extraction

Nucleic acids were isolated from 10 to 25 mg of snap-frozen tumor by standard proteinase K digestion and phenol or TRIzol/chloroform extraction for genomic DNA and total RNA, respectively. Germline DNA was extracted from 2 mL of whole blood using the QuickGene610L Kit (FujiFILM) according to the manufacturer's protocol. RT-PCR of tumor RNA using specific oligonucleotide primers and probe was performed as previously described (48).

WGS

WGS was performed by using the Illumina HiSeq2000 sequencing system (Illumina, Inc.). To prepare short-insert paired-end libraries, the TruSeq Sample Preparation Kit protocol (Illumina) was used with minor modifications. Briefly, 2.0 μg of genomic DNA was sheared on a Covaris E220 ultrasonicator (Covaris) and size-selected using AMPure XP beads (Agencourt; Beckman Coulter) to obtain fragments of approximately 450 bp. The fragmented DNA was end-repaired, adenylated, and ligated to Illumina-specific paired-end adaptors. Each library was sequenced in 2 × 101 bp paired-end mode on a HiSeq2000 flow-cell v3 instrument according to standard Illumina procedures, generating minimal average coverage of 35× for the tumor samples and 25× for the germline samples. Across the entire cohort, 96% of the genomic regions had ≥20× coverage.

Data are available in the European Genome-phenome Archive with the study accession numbers EGAS00001000855 (Institut Curie cohort) and EGAS00001000839 (St. Jude cohort).

Analysis of WGS Data

All samples were processed using the same analysis pipeline. Correspondence between sample and analysis numbers is indicated in Supplementary Table S2B. WGS mapping, coverage, and quality assessment; SNVs; detection of small indels; tier annotation for sequence mutations; and prediction of adverse effects of missense mutations were previously described (14, 15). SVs were analyzed by using the program CREST (13); CNAs were identified by comparing the read depth of matched tumor versus normal tissue and were analyzed by using the CONSERTING algorithm (COpy Number SEgmentation by Regression Tree In Next-Gen sequencing). The reference human genome assembly GRCh37-lite was used to map all samples. We used the program cghMCR to identify recurrent copy-number gain or loss. For this analysis, we excluded 6 cases that showed a highly fragmented copy-number variation (CNV) profile across the genome. These cases have a high number of CNAs across the genome not supported by corresponding SVs. In our experience, these readings are artifacts caused by library construction. The six cases had a total of 182,433 CNV segments, compared with the 16,354 in the remaining 108 cases. We also excluded 3 cases with likely low tumor cell content as suggested by the low number of SNVs and mutation-supporting reads (Supplementary Table S2A). We also excluded from this analysis CNVs identified in the T-cell receptor locus that are caused by physiologic rearrangements of T cells in germline samples. Thus, 16,036 CNV segments from 108 cases were used as input for this analysis. Tumor purity was estimated using loss of heterozygosity, copy-number change, and mutant allele fraction of SNVs, as previously described (25).

The background mutation rate was calculated by dividing the number of high-quality tier-3 SNVs by the total length of tier-3 regions covered at least 10× (Supplementary Table S1). The SVs detected within EWSR1 or ETS genes were consistent with the gene fusions defined by RT-PCR results in all cases. The only case (case SJ001301) that could not be investigated by RT-PCR was shown by WGS to harbor an EWSR1–FLI1 rearrangement. EWSR1 and ETS genomic rearrangements were undetectable despite positive RT-PCR results in six samples. Three of these six cases had low tumor purity, and one had uneven sequencing coverage.

Chromothripsis was analyzed using the criteria proposed by Korbel and Campbell (49). Oscillating patterns of copy-number states were manually inspected using the CIRCOS plots, and statistical tests were applied to evaluate clustering of breakpoints and randomness of DNA fragment joins. For detecting clustering of breakpoints, we applied the Bartlett goodness-of-fit test for exponential distribution to see if there was a strong departure from the null hypothesis, consistent with the chromothripsis hypothesis. For assessing randomness of DNA fragment joins, we applied the goodness-of-fit tests to evaluate if there was no significant departure from the multinomial distribution with equal probabilities, consistent with the chromothripsis hypothesis. This test was carried out for both intrachromosomal and interchromosomal breakpoints when applicable.

Validation Sequencing

For 18 tumor samples (SJEWS001301-1320), the genomic coordinates of putative alterations identified by WGS, including SNVs, SVs, and indels, were used to generate a Nimbelgen Seqcap EZ bait set for enrichment of targeted regions (Roche). The baits were hybridized to Truseq sample libraries (Illumina) prepared from amplified genomic DNA (Roche). Pooled samples were sequenced on a HiSeq 2000 by using the paired-end multiplexed 100-cycle protocol. Resulting data were converted to FASTQ files by using CASAVA 1.8.2 (Illumina) and mapped with the Burrows–Wheeler Aligner before pipeline analysis. Of the 6,659 somatic SNVs identified in the 18 cases, we were able to design a validation assay by custom capture for 6,042. Of these, 5,779 were validated as somatic mutations (overall validation rate, 95.6%).

In addition, a custom TrueSeqAmplicon (Illumina) focused on the whole coding sequence was designed for 8 genes that exhibited at least two somatic changes in the discovery cohort (STAG2, TP53, RYR2, MACF1, DIRAS1, SPTBN5, PCDH10, and CREBRF). Sequencing libraries were prepared following the manufacturer's protocol, and barcoded amplicons were multiplexed and sequenced on Illumina2500 HiSeqFast flow cells; mean target coverage was 98%, and mean number of mapped reads was 6.1 million. All mutations detected by WGS were confirmed.

CDKN2A status of cell lines was verified by PCR amplification spanning all four CDKN2A exons, as previously described (50), followed by the Sanger sequencing of the amplicons. Across the follow-up cohort, CDKN2A status was determined by real-time PCR on genomic DNA. Two sets of primers (Exon1a and Exon2) were used to detect CDKN2A: CDKN2A-ex1A_F: GGCTGGCTGGTCACCAGA, CDKN2A-ex1A_probe: FAM 5′-ATGGAGCCTTCGGCTGACTGGCT-3′BHQ1, CDKN2A-ex1A_R: CGCCCGCACCTCCTCTAC; CDKN2A-ex2_F: GGCTCTACACAAGCTTCCTTTCC, CDKN2A-ex2_probe: FAM 5′-CATGCCGGCCCCCACCCT-3′ BHQ1, CDKN2A-ex2_R: CCTGCCAGAGAGAACAGAATGG and were respectively normalized to TGFBR2 genomic levels (located on chr 3, the most stable chromosome across the Ewing sarcoma WGS cohort): TGFBR2_F: GCAAATCTGGTTGCCCTAGCAAGA, TGFBR2_probe: 5′Yakima-Yellow-CCCGTTTGCACATGAGAGGGTAAGT-3′BHQ1, TGFBR2_R: AAAGTGGGTTGGGAGTCACCTGAA. Duplex PCR (CDKN2A/TGFBR2) using TaqMan Universal PCR Master Mix (Life Technologies) was performed, and the mean of normalized CDKN2A Exon1a and Exon2 was calculated (CDKN2AEX1A-2). Ten nontumorigenic germline DNAs were used as controls (average CDKN2AEX1A-2 value set to 1). Eleven tumor samples with loss of CDKN2A and 38 CDKN2A–wild-type tumor samples from the WGS cohort were used to set CDKN2AEX1A-2 heterozygous threshold to 0.7 (no false positive for wild-type CDKN2A samples).

Statistical Analysis

Overall survival was defined as the time from diagnosis to disease-related death or last follow-up. Survival curves were analyzed according to the Kaplan–Meier method and compared using the log-rank test.

We used the SMG test in the MuSiC suite (21) to identify genes significantly enriched for somatic SNVs and indel mutations. This test assigns mutations to seven categories (AT transition, AT transversion, CG transition, CG transversion, CpG transition, CpG transversion, and indel) and then uses statistical methods based on convolution, hypergeometric distribution (Fisher test), and likelihood to combine the category-specific binomials and obtain overall P values. Genes with false discovery rate <0.1 in 2 of the 3 statistical tests were considered significantly mutated.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors' Contributions

Conception and design: F. Tirode, D. Surdez, E.R. Mardis, R.K. Wilson, J. Downing, M. Dyer, J. Zhang, O. Delattre

Development of methodology: D. Surdez, T. Rio-Frio, J. Easton, R.K. Wilson, J. Zhang

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): D. Surdez, M.C. Le Deley, A. Bahrami, E. Lapouble, S. Reynaud, T. Rio-Frio, G. Pierron, O. Oberlin, S. Zaidi, M. Gut, E.R. Mardis, R.K. Wilson, S. Shurtleff, V. Laurence, J. Michon, P. Marec-Bérard, I. Gut, M. Dyer, J. Zhang, O. Delattre

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): F. Tirode, D. Surdez, X. Ma, M. Parker, Z. Zhang, S. Grossetête-Lalami, M. Rusch, G. Wu, X. Chen, G. Lemmon, B. Vadodaria, L. Ding, R.K. Wilson, J. Downing, M. Dyer, J. Zhang, O. Delattre

Writing, review, and/or revision of the manuscript: F. Tirode, D. Surdez, M.C. Le Deley, A. Bahrami, M. Gut, E.R. Mardis, R.K. Wilson, P. Marec-Bérard, I. Gut, J. Downing, M. Dyer, J. Zhang, O. Delattre

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): D. Surdez, S. Grossetête-Lalami, M. Rusch, E. Hedlund, G. Pierron, P. Gupta, B. Vadodaria, J. Easton, J. Downing, J. Zhang, O. Delattre

Study supervision: R.K. Wilson, J. Downing, M. Dyer, J. Zhang, O. Delattre

Grant Support

This work was supported in part by a Cancer Center Support grant (CA21765) and grants to M. Dyer (EY014867, EY018599, and CA168875) from the U.S. National Institutes of Health, and by the American Lebanese Syrian Associated Charities (ALSAC). M. Dyer is a Howard Hughes Medical Institute Investigator. WGS was supported as part of the St. Jude Children's Research Hospital–Washington University Pediatric Cancer Genome Project.

The WGS of French cases was performed by the Centro Nacional de Análisis Genómico (CNAG) in Barcelona and supported by grants from the Institut National de la Santé et de la Recherche Médicale Inserm within the framework of the International Cancer Genome Consortium program. Additional sequencing was performed by the next-generation sequencing platform of the Institut Curie, supported by grants ANR-10-EQPX-03 and ANR10-INBS-09-08 from the Agence Nationale de la Recherche (investissements d'avenir) and by the Canceropôle Ile-de-France.

This work was also supported by grants from the Ligue Nationale Contre Le Cancer (Equipe labellisée), and by European PROVABES (ERA-NET TRANSCAN JTC-2011), ASSET (FP7-HEALTH-2010-259348), and EEC (HEALTH-F2-2013-602856) grants. D. Surdez is supported by the Institut Curie–SIRIC (Site de Recherche Intégrée en Cancérologie) program.

The authors also thank the following associations for their invaluable support: the Société Française des Cancers de l'Enfant, Courir pour Mathieu, Dans les pas du Géant, Olivier Chape, Les Bagouzamanon, Enfants et Santé, and les Amis de Claire.

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.

Acknowledgments

The authors thank Fabien Calvo for continuous support, Virginie Chêne, Stelly Ballet, Heather Mulder, Panduka Nagahawatte, Donald Yergeau, Yongjin Li, Michael Edmonson, Andrew Thrasher, and Carlo Lucchesi for their invaluable help, Peter Brooks for fruitful discussions, Heinrich Kovar for some of the cell lines, Alban Lermine, Nicolas Servant, Philippe Hupé, and Emmanuel Barillot for their help in processing the Next Generation Sequencing data, and UNICANCER for providing access to the clinical databases.

The authors also thank the following clinicians and pathologists for providing samples used in this work: I. Aerts, P. Anract, C. Bergeron, L. Boccon-Gibod, F. Boman, F. Bourdeaut, C. Bouvier, R. Bouvier, L. Brugières, E. Cassagnau, J. Champigneulle, C. Cordonnier, J. M. Coindre, N. Corradini, A. Coulomb-Lhermine, A. De Muret, G. De Pinieux, A.S. Defachelles, A. Deville, F. Dijoud, F. Doz, C. Dufour, K. Fernandez, N. Gaspard, L. Galmiche-Rolland, C. Glorion, A. Gomez-Brouchet, J.M. Guinebretière, H. Jouan, C. Jeanne-Pasquier, B. Kantelip, F. Labrousse, V. Laithier, F. Larousserie, G. Leverger, C. Linassier, P. Mary, G. Margueritte, E. Mascard, A. Moreau, J. Michon, C. Michot, F. Millot, Y. Musizzano, M. Munzer, B. Narciso, O. Oberlin, D. Orbach, H. Pacquement, Y. Perel, B. Petit, M. Peuchmaur, J.Y. Pierga, C. Piguet, S. Piperno-Neumann, E. Plouvier, D. Ranchere-Vince, J. Rivel, C. Rouleau, H. Rubie, H. Sartelet, G. Schleiermacher, C. Schmitt, N. Sirvent, D. Sommelet, P. Terrier, R. Tichit, J. Vannier, J. M. Vignaud, and V. Verkarre.

Footnotes

  • Note: Supplementary data for this article are available at Cancer Discovery Online (http://cancerdiscovery.aacrjournals.org/).

  • Received June 16, 2014.
  • Revision received September 4, 2014.
  • Accepted September 5, 2014.
  • ©2014 American Association for Cancer Research.

References

  1. 1.↵
    1. Lessnick SL,
    2. Ladanyi M
    . Molecular pathogenesis of Ewing sarcoma: new therapeutic and transcriptional targets. Annu Rev Pathol 2012;7:145–59.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Delattre O,
    2. Zucman J,
    3. Plougastel B,
    4. Desmaze C,
    5. Melot T,
    6. Peter M,
    7. et al.
    Gene fusion with an ETS DNA-binding domain caused by chromosome translocation in human tumours. Nature 1992;359:162–5.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. May WA,
    2. Lessnick SL,
    3. Braun BS,
    4. Klemsz M,
    5. Lewis BC,
    6. Lunsford LB,
    7. et al.
    The Ewing's sarcoma EWS/FLI-1 fusion gene encodes a more potent transcriptional activator and is a more powerful transforming gene than FLI-1. Mol Cell Biol 1993;13:7393–8.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Armengol G,
    2. Tarkkanen M,
    3. Virolainen M,
    4. Forus A,
    5. Valle J,
    6. Bohling T,
    7. et al.
    Recurrent gains of 1q, 8 and 12 in the Ewing family of tumours by comparative genomic hybridization. Br J Cancer 1997;75:1403–9.
    OpenUrlPubMed
  5. 5.↵
    1. Brisset S,
    2. Schleiermacher G,
    3. Peter M,
    4. Mairal A,
    5. Oberlin O,
    6. Delattre O,
    7. et al.
    CGH analysis of secondary genetic changes in Ewing tumors: correlation with metastatic disease in a series of 43 cases. Cancer Genet Cytogenet 2001;130:57–61.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Brownhill SC,
    2. Taylor C,
    3. Burchill SA
    . Chromosome 9p21 gene copy number and prognostic significance of p16 in ESFT. Br J Cancer 2007;96:1914–23.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Hattinger CM,
    2. Pötschger U,
    3. Tarkkanen M,
    4. Squire J,
    5. Zielenska M,
    6. Kiuru-Kuhlefelt S,
    7. et al.
    Prognostic impact of chromosomal aberrations in Ewing tumours. Br J Cancer 2002;86:1763–9.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Kovar H,
    2. Jug G,
    3. Aryee DN,
    4. Zoubek A,
    5. Ambros P,
    6. Gruber B,
    7. et al.
    Among genes involved in the RB dependent cell cycle regulatory cascade, the p16 tumor suppressor gene is frequently lost in the Ewing family of tumors. Oncogene 1997;15:2225–32.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Mackintosh C,
    2. Ordóñez JL,
    3. García-Domínguez DJ,
    4. Sevillano V,
    5. Llombart-Bosch A,
    6. Szuhai K,
    7. et al.
    1q gain and CDT2 overexpression underlie an aggressive and highly proliferative form of Ewing sarcoma. Oncogene 2012;31:1287–98.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Ozaki T,
    2. Paulussen M,
    3. Poremba C,
    4. Brinkschmidt C,
    5. Rerin J,
    6. Ahrens S,
    7. et al.
    Genetic imbalances revealed by comparative genomic hybridization in Ewing tumors. Genes Chromosomes Cancer 2001;32:164–71.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Huang H-Y,
    2. Illei PB,
    3. Zhao Z,
    4. Mazumdar M,
    5. Huvos AG,
    6. Healey JH,
    7. et al.
    Ewing sarcomas with p53 mutation or p16/p14ARF homozygous deletion: a highly lethal subset associated with poor chemoresponse. J Clin Oncol 2005;23:548–58.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Solomon DA,
    2. Kim T,
    3. Diaz-Martinez LA,
    4. Fair J,
    5. Elkahloun AG,
    6. Harris BT,
    7. et al.
    Mutational inactivation of STAG2 causes aneuploidy in human cancer. Science 2011;333:1039–43.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Wang J,
    2. Mullighan CG,
    3. Easton J,
    4. Roberts S,
    5. Heatley SL,
    6. Ma J,
    7. et al.
    CREST maps somatic structural variation in cancer genomes with base-pair resolution. Nat Methods 2011;8:652–4.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Zhang J,
    2. Benavente CA,
    3. McEvoy J,
    4. Flores-Otero J,
    5. Ding L,
    6. Chen X,
    7. et al.
    A novel retinoblastoma therapy from genomic and epigenetic analyses. Nature 2012;481:329–34.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Zhang J,
    2. Ding L,
    3. Holmfeldt L,
    4. Wu G,
    5. Heatley SL,
    6. Payne-Turner D,
    7. et al.
    The genetic basis of early T-cell precursor acute lymphoblastic leukaemia. Nature 2012;481:157–63.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Zucman J,
    2. Melot T,
    3. Desmaze C,
    4. Ghysdael J,
    5. Plougastel B,
    6. Peter M,
    7. et al.
    Combinatorial generation of variable fusion proteins in the Ewing family of tumours. EMBO J 1993;12:4481–7.
    OpenUrlPubMed
  17. 17.↵
    1. Nik-Zainal S,
    2. Alexandrov LB,
    3. Wedge DC,
    4. Van Loo P,
    5. Greenman CD,
    6. Raine K,
    7. et al.
    Mutational processes molding the genomes of 21 breast cancers. Cell 2012;149:979–93.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Remeseiro S,
    2. Cuadrado A,
    3. Losada A
    . Cohesin in development and disease. Dev Camb Engl 2013;140:3715–8.
    OpenUrl
  19. 19.↵
    1. Lawrence MS,
    2. Stojanov P,
    3. Polak P,
    4. Kryukov GV,
    5. Cibulskis K,
    6. Sivachenko A,
    7. et al.
    Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature 2013;499:214–8.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Huether R,
    2. Dong L,
    3. Chen X,
    4. Wu G,
    5. Parker M,
    6. Wei L,
    7. et al.
    The landscape of somatic mutations in epigenetic regulators across 1,000 paediatric cancer genomes. Nat Commun 2014;5:3630.
    OpenUrlPubMed
  21. 21.↵
    1. Dees ND,
    2. Zhang Q,
    3. Kandoth C,
    4. Wendl MC,
    5. Schierding W,
    6. Koboldt DC,
    7. et al.
    MuSiC: identifying mutational significance in cancer genomes. Genome Res 2012;22:1589–98.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Molenaar JJ,
    2. Koster J,
    3. Zwijnenburg DA,
    4. van Sluis P,
    5. Valentijn LJ,
    6. van der Ploeg I,
    7. et al.
    Sequencing of neuroblastoma identifies chromothripsis and defects in neuritogenesis genes. Nature 2012;483:589–93.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Pugh TJ,
    2. Morozova O,
    3. Attiyeh EF,
    4. Asgharzadeh S,
    5. Wei JS,
    6. Auclair D,
    7. et al.
    The genetic landscape of high-risk neuroblastoma. Nat Genet 2013;45:279–84.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Sausen M,
    2. Leary RJ,
    3. Jones S,
    4. Wu J,
    5. Reynolds CP,
    6. Liu X,
    7. et al.
    Integrated genomic analyses identify ARID1A and ARID1B alterations in the childhood cancer neuroblastoma. Nat Genet 2013;45:12–7.
    OpenUrlPubMed
  25. 25.↵
    1. Chen X,
    2. Stewart E,
    3. Shelat AA,
    4. Qu C,
    5. Bahrami A,
    6. Hatley M,
    7. et al.
    Targeting oxidative stress in embryonal rhabdomyosarcoma. Cancer Cell 2013;24:710–24.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Shern JF,
    2. Chen L,
    3. Chmielecki J,
    4. Wei JS,
    5. Patidar R,
    6. Rosenberg M,
    7. 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.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Pugh TJ,
    2. Weeraratne SD,
    3. Archer TC,
    4. Pomeranz Krummel DA,
    5. Auclair D,
    6. Bochicchio J,
    7. et al.
    Medulloblastoma exome sequencing uncovers subtype-specific somatic mutations. Nature 2012;488:106–10.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Robinson G,
    2. Parker M,
    3. Kranenburg TA,
    4. Lu C,
    5. Chen X,
    6. Ding L,
    7. et al.
    Novel mutations target distinct subgroups of medulloblastoma. Nature 2012;488:43–8.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Jones DTW,
    2. Hutter B,
    3. Jäger N,
    4. Korshunov A,
    5. Kool M,
    6. Warnatz H-J,
    7. et al.
    Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma. Nat Genet 2013;45:927–32.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Schwartzentruber J,
    2. Korshunov A,
    3. Liu X-Y,
    4. Jones DTW,
    5. Pfaff E,
    6. Jacob K,
    7. et al.
    Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma. Nature 2012;482:226–31.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Chen X,
    2. Bahrami A,
    3. Pappo A,
    4. Easton J,
    5. Dalton J,
    6. Hedlund E,
    7. et al.
    Recurrent somatic structural variations contribute to tumorigenesis in pediatric osteosarcoma. Cell Rep 2014;7:104–12.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Mugneret F,
    2. Lizard S,
    3. Aurias A,
    4. Turc-Carel C
    . Chromosomes in Ewing's sarcoma. II. Nonrandom additional changes, trisomy 8 and der(16)t(1;16). Cancer Genet Cytogenet 1988;32:239–45.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Balbás-Martínez C,
    2. Sagrera A,
    3. Carrillo-de-Santa-Pau E,
    4. Earl J,
    5. Márquez M,
    6. Vazquez M,
    7. et al.
    Recurrent inactivation of STAG2 in bladder cancer is not associated with aneuploidy. Nat Genet 2013;45:1464–9.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Guo G,
    2. Sun X,
    3. Chen C,
    4. Wu S,
    5. Huang P,
    6. Li Z,
    7. et al.
    Whole-genome and whole-exome sequencing of bladder cancer identifies frequent alterations in genes involved in sister chromatid cohesion and segregation. Nat Genet 2013;45:1459–63.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Kon A,
    2. Shih L-Y,
    3. Minamino M,
    4. Sanada M,
    5. Shiraishi Y,
    6. Nagata Y,
    7. et al.
    Recurrent mutations in multiple components of the cohesin complex in myeloid neoplasms. Nat Genet 2013;45:1232–7.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Taylor CF,
    2. Platt FM,
    3. Hurst CD,
    4. Thygesen HH,
    5. Knowles MA
    . Frequent inactivating mutations of STAG2 in bladder cancer are associated with low tumour grade and stage and inversely related to chromosomal copy number changes. Hum Mol Genet 2014;23:1964–74.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Xiao T,
    2. Wallace J,
    3. Felsenfeld G
    . Specific sites in the C terminus of CTCF interact with the SA2 subunit of the cohesin complex and are required for cohesin-dependent insulation activity. Mol Cell Biol 2011;31:2174–83.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Witcher M,
    2. Emerson BM
    . Epigenetic silencing of the p16(INK4a) tumor suppressor is associated with loss of CTCF binding and a chromatin boundary. Mol Cell 2009;34:271–84.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Lopez-Guerrero JA,
    2. Pellin A,
    3. Noguera R,
    4. Carda C,
    5. Llombart-Bosch A
    . Molecular analysis of the 9p21 locus and p53 genes in Ewing family tumors. Lab Invest 2001;81:803–14.
    OpenUrlPubMed
  40. 40.↵
    1. Tsuchiya T,
    2. Sekine K,
    3. Hinohara S,
    4. Namiki T,
    5. Nobori T,
    6. Kaneko Y
    . Analysis of the p16INK4, p14ARF, p15, TP53, and MDM2 genes and their prognostic implications in osteosarcoma and Ewing sarcoma. Cancer Genet Cytogenet 2000;120:91–8.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Patel N,
    2. Black J,
    3. Chen X,
    4. Marcondes AM,
    5. Grady WM,
    6. Lawlor ER,
    7. et al.
    DNA methylation and gene expression profiling of Ewing sarcoma primary tumors reveal genes that are potential targets of epigenetic inactivation. Sarcoma 2012;2012:498472.
    OpenUrlPubMed
  42. 42.↵
    1. Béguelin W,
    2. Popovic R,
    3. Teater M,
    4. Jiang Y,
    5. Bunting KL,
    6. Rosen M,
    7. et al.
    EZH2 is required for germinal center formation and somatic EZH2 mutations promote lymphoid transformation. Cancer Cell 2013;23:677–92.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Berg T,
    2. Thoene S,
    3. Yap D,
    4. Wee T,
    5. Schoeler N,
    6. Rosten P,
    7. et al.
    A transgenic mouse model demonstrating the oncogenic role of mutations in the polycomb-group gene EZH2 in lymphomagenesis. Blood 2014;123:3914–24.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Kridel R,
    2. Sehn LH,
    3. Gascoyne RD
    . Pathogenesis of follicular lymphoma. J Clin Invest 2012;122:3424–31.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Bejar R,
    2. Stevenson K,
    3. Abdel-Wahab O,
    4. Galili N,
    5. Nilsson B,
    6. Garcia-Manero G,
    7. et al.
    Clinical effect of point mutations in myelodysplastic syndromes. N Engl J Med 2011;364:2496–506.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Brohl AS,
    2. Solomon DA,
    3. Chang W,
    4. Wang J,
    5. Song Y,
    6. Sindiri S,
    7. et al.
    The genomic landscape of the Ewing sarcoma family of tumors reveals recurrent STAG2 mutation. PLoS Genet 2014;10:e1004475.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Ladenstein R,
    2. Pötschger U,
    3. Le Deley MC,
    4. Whelan J,
    5. Paulussen M,
    6. Oberlin O,
    7. et al.
    Primary disseminated multifocal Ewing sarcoma: results of the Euro-EWING 99 trial. J Clin Oncol 2010;28:3284–91.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Peter M,
    2. Gilbert E,
    3. Delattre O
    . A multiplex real-time pcr assay for the detection of gene fusions observed in solid tumors. Lab Invest 2001;81:905–12.
    OpenUrlPubMed
  49. 49.↵
    1. Korbel JO,
    2. Campbell PJ
    . Criteria for inference of chromothripsis in cancer genomes. Cell 2013;152:1226–36.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Louis-Brennetot C,
    2. Coindre J-M,
    3. Ferreira C,
    4. Pérot G,
    5. Terrier P,
    6. Aurias A
    . The CDKN2A/CDKN2B/CDK4/CCND1 pathway is pivotal in well-differentiated and dedifferentiated liposarcoma oncogenesis: an analysis of 104 tumors. Genes Chromosomes Cancer 2011;50:896–907.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Cancer Discovery: 4 (11)
November 2014
Volume 4, Issue 11
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Cancer Discovery article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Genomic Landscape of Ewing Sarcoma Defines an Aggressive Subtype with Co-Association of STAG2 and TP53 Mutations
(Your Name) has forwarded a page to you from Cancer Discovery
(Your Name) thought you would be interested in this article in Cancer Discovery.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Genomic Landscape of Ewing Sarcoma Defines an Aggressive Subtype with Co-Association of STAG2 and TP53 Mutations
Franck Tirode, Didier Surdez, Xiaotu Ma, Matthew Parker, Marie Cécile Le Deley, Armita Bahrami, Zhaojie Zhang, Eve Lapouble, Sandrine Grossetête-Lalami, Michael Rusch, Stéphanie Reynaud, Thomas Rio-Frio, Erin Hedlund, Gang Wu, Xiang Chen, Gaelle Pierron, Odile Oberlin, Sakina Zaidi, Gordon Lemmon, Pankaj Gupta, Bhavin Vadodaria, John Easton, Marta Gut, Li Ding, Elaine R. Mardis, Richard K. Wilson, Sheila Shurtleff, Valérie Laurence, Jean Michon, Perrine Marec-Bérard, Ivo Gut, James Downing, Michael Dyer, Jinghui Zhang and Olivier Delattre for the St. Jude Children's Research Hospital–Washington University Pediatric Cancer Genome Project and the International Cancer Genome Consortium
Cancer Discov November 1 2014 (4) (11) 1342-1353; DOI: 10.1158/2159-8290.CD-14-0622

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Genomic Landscape of Ewing Sarcoma Defines an Aggressive Subtype with Co-Association of STAG2 and TP53 Mutations
Franck Tirode, Didier Surdez, Xiaotu Ma, Matthew Parker, Marie Cécile Le Deley, Armita Bahrami, Zhaojie Zhang, Eve Lapouble, Sandrine Grossetête-Lalami, Michael Rusch, Stéphanie Reynaud, Thomas Rio-Frio, Erin Hedlund, Gang Wu, Xiang Chen, Gaelle Pierron, Odile Oberlin, Sakina Zaidi, Gordon Lemmon, Pankaj Gupta, Bhavin Vadodaria, John Easton, Marta Gut, Li Ding, Elaine R. Mardis, Richard K. Wilson, Sheila Shurtleff, Valérie Laurence, Jean Michon, Perrine Marec-Bérard, Ivo Gut, James Downing, Michael Dyer, Jinghui Zhang and Olivier Delattre for the St. Jude Children's Research Hospital–Washington University Pediatric Cancer Genome Project and the International Cancer Genome Consortium
Cancer Discov November 1 2014 (4) (11) 1342-1353; DOI: 10.1158/2159-8290.CD-14-0622
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Methods
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Grant Support
    • Acknowledgments
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • Poor Survival and Genetic Features of Black Patients with AML
  • Type II RAFi and MEKi to Treat MAPK-Addicted Cancers
  • ATLAS-Identified Neoantigens May Improve Tumor Immunotherapy
Show more Research Articles
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook   Twitter   LinkedIn   YouTube   RSS

Articles

  • OnlineFirst
  • Current Issue
  • Past Issues

Info For

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About Cancer Discovery

  • About the Journal
  • Editors
  • Journal Sections
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Cancer Discovery
eISSN: 2159-8290
ISSN: 2159-8274

Advertisement