Didier Surdez

Systematic identification of genomic markers of drug sensitivity in cancer cells

Clinical responses to anticancer therapies are often restricted to a subset of patients. In some cases, mutated cancer genes are potent biomarkers for responses to targeted agents. Here, to uncover new biomarkers of sensitivity and resistance to cancer therapeutics, we screened a panel of several hundred cancer cell lines—which represent much of the tissue-type and genetic diversity of human cancers—with 130 drugs under clinical and preclinical investigation. In aggregate, we found that mutated cancer genes were associated with cellular response to most currently available cancer drugs. Classic oncogene addiction paradigms were modified by additional tissue-specific or expression biomarkers, and some frequently mutated genes were associated with sensitivity to a broad range of therapeutic agents. Unexpected relationships were revealed, including the marked sensitivity of Ewing’s sarcoma cells harbouring the EWS (also known as EWSR1)-FLI1 gene translocation to poly(ADP-ribose) polymerase (PARP) inhibitors. By linking drug activity to the functional complexity of cancer genomes, systematic pharmacogenomic profiling in cancer cell lines provides a powerful biomarker discovery platform to guide rational cancer therapeutic strategies.

Genomic landscape of Ewing sarcoma defines an aggressive subtype with co-association of STAG2 and TP53 mutations

UNLABELLED 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.

De novo motif identification improves the accuracy of predicting transcription factor binding sites in ChIP-Seq data analysis.

Dramatic progress in the development of next-generation sequencing technologies has enabled accurate genome-wide characterization of the binding sites of DNA-associated proteins. This technique, baptized as ChIP-Seq, uses a combination of chromatin immunoprecipitation and massively parallel DNA sequencing. Other published tools that predict binding sites from ChIP-Seq data use only positional information of mapped reads. In contrast, our algorithm MICSA (Motif Identification for ChIP-Seq Analysis) combines this source of positional information with information on motif occurrences to better predict binding sites of transcription factors (TFs). We proved the greater accuracy of MICSA with respect to several other tools by running them on datasets for the TFs NRSF, GABP, STAT1 and CTCF. We also applied MICSA on a dataset for the oncogenic TF EWS-FLI1. We discovered >2000 binding sites and two functionally different binding motifs. We observed that EWS-FLI1 can activate gene transcription when (i) its binding site is located in close proximity to the gene transcription start site (up to ∼150 kb), and (ii) it contains a microsatellite sequence. Furthermore, we observed that sites without microsatellites can also induce regulation of gene expression—positively as often as negatively—and at much larger distances (up to ∼1 Mb).

Chimeric EWSR1-FLI1 regulates the Ewing sarcoma susceptibility gene EGR2 via a GGAA microsatellite

Deciphering the ways in which somatic mutations and germline susceptibility variants cooperate to promote cancer is challenging. Ewing sarcoma is characterized by fusions between EWSR1 and members of the ETS gene family, usually EWSR1-FLI1, leading to the generation of oncogenic transcription factors that bind DNA at GGAA motifs. A recent genome-wide association study identified susceptibility variants near EGR2. Here we found that EGR2 knockdown inhibited proliferation, clonogenicity and spheroidal growth in vitro and induced regression of Ewing sarcoma xenografts. Targeted germline deep sequencing of the EGR2 locus in affected subjects and controls identified 291 Ewing-associated SNPs. At rs79965208, the A risk allele connected adjacent GGAA repeats by converting an interspaced GGAT motif into a GGAA motif, thereby increasing the number of consecutive GGAA motifs and thus the EWSR1-FLI1–dependent enhancer activity of this sequence, with epigenetic characteristics of an active regulatory element. EWSR1-FLI1 preferentially bound to the A risk allele, which increased global and allele-specific EGR2 expression. Collectively, our findings establish cooperation between a dominant oncogene and a susceptibility variant that regulates a major driver of Ewing sarcomagenesis.

Targeting the EWSR1-FLI1 Oncogene-Induced Protein Kinase PKC-β Abolishes Ewing Sarcoma Growth

Ewing sarcoma is a rare but aggressive disease most common in young adults. This cancer is driven by a unique chimeric fusion oncogene but targeted strategies have been elusive. Here we report the identification of the protein kinase PKC-ß (PRKCB) as a disease-specific druggable target for treatment of Ewing sarcoma. We found that transcriptional activation of PRKCB was directly regulated by the chimeric fusion oncogene EWSR1-FLI1 that drives this cancer. PRKCB phosphorylated histone H3T6 to permit global maintenance of H3K4 trimethylation at a variety of gene promoters. PRKCB loss induced apoptosis in vitro and prevented tumor growth in vivo. Gene expression profiling revealed a strong overlap between genes modulated by EWSR1-FLI1 and PRKCB in regulating crucial signaling pathways. Taken together, our findings offer a preclinical proof-of-concept for PRKCB as a promising therapeutic target in Ewing sarcoma.

DNA methylation heterogeneity defines a disease spectrum in Ewing sarcoma

Developmental tumors in children and young adults carry few genetic alterations, yet they have diverse clinical presentation. Focusing on Ewing sarcoma, we sought to establish the prevalence and characteristics of epigenetic heterogeneity in genetically homogeneous cancers. We performed genome-scale DNA methylation sequencing for a large cohort of Ewing sarcoma tumors and analyzed epigenetic heterogeneity on three levels: between cancers, between tumors, and within tumors. We observed consistent DNA hypomethylation at enhancers regulated by the disease-defining EWS-FLI1 fusion protein, thus establishing epigenomic enhancer reprogramming as a ubiquitous and characteristic feature of Ewing sarcoma. DNA methylation differences between tumors identified a continuous disease spectrum underlying Ewing sarcoma, which reflected the strength of an EWS-FLI1 regulatory signature and a continuum between mesenchymal and stem cell signatures. There was substantial epigenetic heterogeneity within tumors, particularly in patients with metastatic disease. In summary, our study provides a comprehensive assessment of epigenetic heterogeneity in Ewing sarcoma and thereby highlights the importance of considering nongenetic aspects of tumor heterogeneity in the context of cancer biology and personalized medicine.

Systems biology of Ewing sarcoma: a network model of EWS-FLI1 effect on proliferation and apoptosis

Ewing sarcoma is the second most frequent pediatric bone tumor. In most of the patients, a chromosomal translocation leads to the expression of the EWS-FLI1 chimeric transcription factor that is the major oncogene in this pathology. Relative genetic simplicity of Ewing sarcoma makes it particularly attractive for studying cancer in a systemic manner. Silencing EWS-FLI1 induces cell cycle alteration and ultimately leads to apoptosis, but the exact molecular mechanisms underlying this phenotype are unclear. In this study, a network linking EWS-FLI1 to cell cycle and apoptosis phenotypes was constructed through an original method of network reconstruction. Transcriptome time-series after EWS-FLI1 silencing were used to identify core modulated genes by an original scoring method based on fitting expression profile dynamics curves. Literature data mining was then used to connect these modulated genes into a network. The validity of a subpart of this network was assessed by siRNA/RT-QPCR experiments on four additional Ewing cell lines and confirmed most of the links. Based on the network and the transcriptome data, CUL1 was identified as a new potential target of EWS-FLI1. Altogether, using an original methodology of data integration, we provide the first version of EWS-FLI1 network model of cell cycle and apoptosis regulation.

Heterogeneity of neuroblastoma cell identity defined by transcriptional circuitries

Neuroblastoma is a tumor of the peripheral sympathetic nervous system, derived from multipotent neural crest cells (NCCs). To define core regulatory circuitries (CRCs) controlling the gene expression program of neuroblastoma, we established and analyzed the neuroblastoma super-enhancer landscape. We discovered three types of identity in neuroblastoma cell lines: a sympathetic noradrenergic identity, defined by a CRC module including the PHOX2B, HAND2 and GATA3 transcription factors (TFs); an NCC-like identity, driven by a CRC module containing AP-1 TFs; and a mixed type, further deconvoluted at the single-cell level. Treatment of the mixed type with chemotherapeutic agents resulted in enrichment of NCC-like cells. The noradrenergic module was validated by ChIP-seq. Functional studies demonstrated dependency of neuroblastoma with noradrenergic identity on PHOX2B, evocative of lineage addiction. Most neuroblastoma primary tumors express TFs from the noradrenergic and NCC-like modules. Our data demonstrate a previously unknown aspect of tumor heterogeneity relevant for neuroblastoma treatment strategies.

The occurrence of intracranial rhabdoid tumours in mice depends on temporal control of Smarcb1 inactivation

Rhabdoid tumours (RTs) are highly aggressive tumours of infancy, frequently localized in the central nervous system (CNS) where they are termed atypical teratoid/rhabdoid tumours (AT/RTs) and characterized by bi-allelic inactivation of the SMARCB1 tumour suppressor gene. In this study, by temporal control of tamoxifen injection in Smarcb1flox/flox;Rosa26-CreERT2 mice, we explore the phenotypes associated with Smarcb1 inactivation at different developmental stages. Injection before E6, at birth or at 2 months of age recapitulates previously described phenotypes including embryonic lethality, hepatic toxicity or development of T-cell lymphomas, respectively. Injection between E6 and E10 leads to high penetrance tumours, mainly intra-cranial, with short delays (median: 3 months). These tumours demonstrate anatomical, morphological and gene expression profiles consistent with those of human AT/RTs. Moreover, intra- and inter-species comparisons of tumours reveal that human and mouse RTs can be split into different entities that may underline the variety of RT cells of origin.

Combined experience of six independent laboratories attempting to create an Ewing sarcoma mouse model

Ewing sarcoma (ES) involves a tumor-specific chromosomal translocation that produces the EWS-FLI1 protein, which is required for the growth of ES cells both in vitro and in vivo. However, an EWS-FLI1-driven transgenic mouse model is not currently available. Here, we present data from six independent laboratories seeking an alternative approach to express EWS-FLI1 in different murine tissues. We used the Runx2, Col1a2.3, Col1a3.6, Prx1, CAG, Nse, NEFL, Dermo1, P0, Sox9 and Osterix promoters to target EWS-FLI1 or Cre expression. Additional approaches included the induction of an endogenous chromosomal translocation, in utero knock-in, and the injection of Cre-expressing adenovirus to induce EWS-FLI1 expression locally in multiple lineages. Most models resulted in embryonic lethality or developmental defects. EWS-FLI1-induced apoptosis, promoter leakiness, the lack of potential cofactors, and the difficulty of expressing EWS-FLI1 in specific sites were considered the primary reasons for the failed attempts to create a transgenic mouse model of ES.