Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Upregulation of the Mevalonate Pathway through EWSR1-FLI1/EGR2 Regulatory Axis Confers Ewing Cells Exquisite Sensitivity to Statins

Simple Summary

The objective of this project was to search for new dependencies in Ewing sarcoma, a deadly disease for which new therapeutic approaches are urgently needed. A pharmacological screening of off-patent approved drugs (FDA agency) and the investigation of downstream targets of EGR2 were performed. The two approaches showed the MVA pathway as a major dependency in Ewing sarcoma and statin, an inhibitor of this pathway, as a potential new therapeutic agent for the treatment of Ewing sarcoma.

Abstract

Ewing sarcoma (EwS) is an aggressive primary bone cancer in children and young adults characterized by oncogenic fusions between genes encoding FET-RNA-binding proteins and ETS transcription factors, the most frequent fusion being EWSR1-FLI1. We show that EGR2, an Ewing-susceptibility gene and an essential direct target of EWSR1-FLI1, directly regulates the transcription of genes encoding key enzymes of the mevalonate (MVA) pathway. Consequently, Ewing sarcoma is one of the tumors that expresses the highest levels of mevalonate pathway genes. Moreover, genome-wide screens indicate that MVA pathway genes constitute major dependencies of Ewing cells. Accordingly, the statin inhibitors of HMG-CoA-reductase, a rate-limiting enzyme of the MVA pathway, demonstrate cytotoxicity in EwS. Statins induce increased ROS and lipid peroxidation levels, as well as decreased membrane localization of prenylated proteins, such as small GTP proteins. These metabolic effects lead to an alteration in the dynamics of S-phase progression and to apoptosis. Statin-induced effects can be rescued by downstream products of the MVA pathway. Finally, we further show that statins impair tumor growth in different Ewing PDX models. Altogether, the data show that statins, which are off-patent, well-tolerated, and inexpensive compounds, should be strongly considered in the therapeutic arsenal against this deadly childhood disease.

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

Intra- and extra-cranial BCOR-ITD tumours are separate entities within the BCOR-rearranged family

BCOR-ITD tumours form an emerging family of aggressive entities with an internal tandem duplication (ITD) in the last exon of the BCOR gene. The family includes cerebral tumours, termed central nervous system BCOR-ITD (CNS BCOR-ITD), and sarcomatous types described in the kidney as clear cell sarcoma of the kidney (CCSK), in the endometrium as high-grade endometrial stromal sarcoma, and in the bone and soft tissue as undifferentiated round cell sarcoma or primitive myxoid mesenchymal tumour of infancy. Based on a series of 33 retrospective cases, including 10 CNS BCOR-ITD and 23 BCOR-ITD sarcomas, we interrogated the homogeneity of the entity regarding clinical, radiological, and histopathological findings, and molecular signatures. Whole-transcriptomic sequencing and DNA methylation profiling were used for unsupervised clustering. BCOR-ITD tumours mostly affected young children with a median age at diagnosis of 2.1 years (range 0-62.4). Median overall survival was 3.9 years and progression-free survival was 1.4 years. This dismal prognosis is shared among tumours in all locations except CCSK. Histopathological review revealed marked differences between CNS BCOR-ITD and BCOR-ITD sarcomas. These two groups were consistently segregated by unsupervised clustering of expression (n = 22) and DNA methylation (n = 21) data. Proximity between the two groups may result from common somatic changes within key pathways directly related to the novel activity of the ITD itself. Conversely, comparison of gene signatures with single-cell RNA-Seq atlases suggests that the distinction between BCOR-ITD sarcomas and CNS BCOR-ITD may result from differences in cells of origin.

Unraveling Ewing sarcoma tumorigenesis originating from patient-derived Mesenchymal Stem Cells

Ewing sarcoma (EwS) is characterized by pathognomonic translocations, most frequently fusing EWSR1 with FLI1. An estimated 30% of EwS tumors also display genetic alterations in STAG2, TP53, or CDKN2A (SPC). Numerous attempts to develop relevant EwS models from primary human cells have been unsuccessful in faithfully recapitulating the phenotypic, transcriptomic and epigenetic features of EwS. In this study, by engineering the t(11;22)(q24;q12) translocation together with a combination of SPC mutations, we generated a wide collection of immortalized cells (EWIma cells) tolerating EWSR1-FLI1 expression from primary mesenchymal stem cells (MSC) derived from an EwS patient. Within this model, SPC alterations strongly favored EwS oncogenicity. Xenograft experiments with independent EWIma cells induced tumors and metastases in mice, which displayed bona fide features of EwS. EWIma cells presented balanced but also more complex translocation profiles mimicking chromoplexy, which is frequently observed in EwS and other cancers. Collectively, these results demonstrate that bone marrow-derived MSCs are a source of origin for EwS and also provide original experimental models to investigate Ewing sarcomagenesis.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

STAG2 mutations alter CTCF-anchored loop extrusion, reduce cis-regulatory interactions and EWSR1-FLI1 activity in Ewing sarcoma

STAG2, a cohesin family gene, is among the most recurrently mutated genes in cancer. STAG2 loss of function (LOF) is associated with aggressive behavior in Ewing sarcoma, a childhood cancer driven by aberrant transcription induced by the EWSR1-FLI1 fusion oncogene. Here, using isogenic Ewing cells, we show that, while STAG2 LOF profoundly changes the transcriptome, it does not significantly impact EWSR1-FLI1, CTCF/cohesin, or acetylated H3K27 DNA binding patterns. In contrast, it strongly alters the anchored dynamic loop extrusion process at boundary CTCF sites and dramatically decreases promoter-enhancer interactions, particularly affecting the expression of genes regulated by EWSR1-FLI1 at GGAA microsatellite neo-enhancers. Down-modulation of cis-mediated EWSR1-FLI1 activity, observed in STAG2-LOF conditions, is associated with enhanced migration and invasion properties of Ewing cells previously observed in EWSR1-FLI1low cells. Our study illuminates a process whereby STAG2-LOF fine-tunes the activity of an oncogenic transcription factor through altered CTCF-anchored loop extrusion and cis-mediated enhancer mechanisms.

Abstract PR09: STAG2 mutations alter topologic organization of the genome and cis-mediated interactions

STAG2 is one of the most recurrently mutated genes in human cancer, including in Ewing sarcoma, an aggressive bone tumor driven by the chimeric EWSR1-FLI1 transcription factor. We and others have shown that STAG2 mutations occurs in approximately 20% of Ewing sarcoma, leading most frequently to complete loss of function (LOF) of its protein. We also showed that STAG2 mutations were associated with poor outcome for these patients. However, functional mechanisms associated with these mutations are still poorly understood. STAG2 encodes an integral member of the cohesin complex, a ring-shaped multiprotein structure, which is important for proper sister chromatid cohesion and release during mitosis. Cohesin, together with CTCF, is also essential to shape the architecture of the genome through its ring structure that allows for chromatin loops formation. Here, we investigated STAG2 and STAG1 functions using isogenic proficient and deficient Ewing sarcoma models. Using extensive characterization of these models by transcriptomics and epigenetics approaches, we demonstrate that STAG2 LOF has a major impact on the expression of genes associated with enhancer features, highlighting STAG2 and EWSR1-FLI1 convergent effects. Using chromatin conformation capture methods, we further show that this transcriptomic modulation results from a STAG2 LOF specific topologic alteration of the genome. This work therefore unravels a new fundamental role of STAG2 in topologic organization of the genome and in gene regulation.

This abstract is also being presented as Poster B53.

Citation Format: Didier Surdez, Sakina Zaidi, Sandrine Grossetête, Virginie Raynal, Sylvain Baulande, Véronique Hill, Olivier Delattre. STAG2 mutations alter topologic organization of the genome and cis-mediated interactions [abstract]. In: Proceedings of the AACR Special Conference on the Advances in Pediatric Cancer Research; 2019 Sep 17-20; Montreal, QC, Canada. Philadelphia (PA): AACR; Cancer Res 2020;80(14 Suppl):Abstract nr PR09.

Cooperation of cancer drivers with regulatory germline variants shapes clinical outcomes

Pediatric malignancies including Ewing sarcoma (EwS) feature a paucity of somatic alterations except for pathognomonic driver-mutations that cannot explain overt variations in clinical outcome. Here, we demonstrate in EwS how cooperation of dominant oncogenes and regulatory germline variants determine tumor growth, patient survival and drug response. Binding of the oncogenic EWSR1-FLI1 fusion transcription factor to a polymorphic enhancer-like DNA element controls expression of the transcription factor MYBL2 mediating these phenotypes. Whole-genome and RNA sequencing reveals that variability at this locus is inherited via the germline and is associated with variable inter-tumoral MYBL2 expression. High MYBL2 levels sensitize EwS cells for inhibition of its upstream activating kinase CDK2 in vitro and in vivo, suggesting MYBL2 as a putative biomarker for anti-CDK2-therapy. Collectively, we establish cooperation of somatic mutations and regulatory germline variants as a major determinant of tumor progression and highlight the importance of integrating the regulatory genome in precision medicine.

FUNGIpath: a tool to assess fungal metabolic pathways predicted by orthology

Background

More and more completely sequenced fungal genomes are becoming available and many more sequencing projects are in progress. This deluge of data should improve our knowledge of the various primary and secondary metabolisms of Fungi, including their synthesis of useful compounds such as antibiotics or toxic molecules such as mycotoxins. Functional annotation of many fungal genomes is imperfect, especially of genes encoding enzymes, so we need dedicated tools to analyze their metabolic pathways in depth.

Description

FUNGIpath is a new tool built using a two-stage approach. Groups of orthologous proteins predicted using complementary methods of detection were collected in a relational database. Each group was further mapped on to steps in the metabolic pathways published in the public databases KEGG and MetaCyc. As a result, FUNGIpath allows the primary and secondary metabolisms of the different fungal species represented in the database to be compared easily, making it possible to assess the level of specificity of various pathways at different taxonomic distances. It is freely accessible at http://www.fungipath.u-psud.fr.

Conclusions

As more and more fungal genomes are expected to be sequenced during the coming years, FUNGIpath should help progressively to reconstruct the ancestral primary and secondary metabolisms of the main branches of the fungal tree of life and to elucidate the evolution of these ancestral fungal metabolisms to various specific derived metabolisms.