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.

PHGDH heterogeneity potentiates cancer cell dissemination and metastasis.

Cancer metastasis requires the transient activation of cellular programs enabling dissemination and seeding in distant organs1. Genetic, transcriptional and translational heterogeneity contributes to this dynamic process2,3. Metabolic heterogeneity has also been observed4, yet its role in cancer progression is less explored. Here we find that the loss of phosphoglycerate dehydrogenase (PHGDH) potentiates metastatic dissemination. Specifically, we find that heterogeneous or low PHGDH expression in primary tumours of patients with breast cancer is associated with decreased metastasis-free survival time. In mice, circulating tumour cells and early metastatic lesions are enriched with Phgdhlow cancer cells, and silencing Phgdh in primary tumours increases metastasis formation. Mechanistically, Phgdh interacts with the glycolytic enzyme phosphofructokinase, and the loss of this interaction activates the hexosamine-sialic acid pathway, which provides precursors for protein glycosylation. As a consequence, aberrant protein glycosylation occurs, including increased sialylation of integrin αvß3, which potentiates cell migration and invasion. Inhibition of sialylation counteracts the metastatic ability of Phgdhlow cancer cells. In conclusion, although the catalytic activity of PHGDH supports cancer cell proliferation, low PHGDH protein expression non-catalytically potentiates cancer dissemination and metastasis formation. Thus, the presence of PHDGH heterogeneity in primary tumours could be considered a sign of tumour aggressiveness.

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.

EIF4EBP1 is transcriptionally upregulated by MYCN and associates with poor prognosis in neuroblastoma.

Neuroblastoma (NB) accounts for 15% of cancer-related deaths in childhood despite considerable therapeutic improvements. While several risk factors, including MYCN amplification and alterations in RAS and p53 pathway genes, have been defined in NB, the clinical outcome is very variable and difficult to predict. Since genes of the mechanistic target of rapamycin (mTOR) pathway are upregulated in MYCN-amplified NB, we aimed to define the predictive value of the mTOR substrate-encoding gene eukaryotic translation initiation factor 4E-binding protein 1 (EIF4EBP1) expression in NB patients. Using publicly available data sets, we found that EIF4EBP1 mRNA expression is positively correlated with MYCN expression and elevated in stage 4 and high-risk NB patients. In addition, high EIF4EBP1 mRNA expression is associated with reduced overall and event-free survival in the entire group of NB patients in three cohorts, as well as in stage 4 and high-risk patients. This was confirmed by monitoring the clinical value of 4EBP1 protein expression, which revealed that high levels of 4EBP1 are significantly associated with prognostically unfavorable NB histology. Finally, functional analyses revealed that EIF4EBP1 expression is transcriptionally controlled by MYCN binding to the EIF4EBP1 promoter in NB cells. Our data highlight that EIF4EBP1 is a direct transcriptional target of MYCN whose high expression is associated with poor prognosis in NB patients. Therefore, EIF4EBP1 may serve to better stratify patients with NB.

Therapeutic targeting of the PLK1-PRC1-axis triggers cell death in genomically silent childhood cancer.

Chromosomal instability (CIN) is a hallmark of cancer1. Yet, many childhood cancers, such as Ewing sarcoma (EwS), feature remarkably 'silent' genomes with minimal CIN2. Here, we show in the EwS model how uncoupling of mitosis and cytokinesis via targeting protein regulator of cytokinesis 1 (PRC1) or its activating polo-like kinase 1 (PLK1) can be employed to induce fatal genomic instability and tumor regression. We find that the EwS-specific oncogenic transcription factor EWSR1-FLI1 hijacks PRC1, which physiologically safeguards controlled cell division, through binding to a proximal enhancer-like GGAA-microsatellite, thereby promoting tumor growth and poor clinical outcome. Via integration of transcriptome-profiling and functional in vitro and in vivo experiments including CRISPR-mediated enhancer editing, we discover that high PRC1 expression creates a therapeutic vulnerability toward PLK1 inhibition that can repress even chemo-resistant EwS cells by triggering mitotic catastrophe.Collectively, our results exemplify how aberrant PRC1 activation by a dominant oncogene can confer malignancy but provide opportunities for targeted therapy, and identify PRC1 expression as an important determinant to predict the efficacy of PLK1 inhibitors being used in clinical trials.

Publisher Correction: Oncogenic hijacking of a developmental transcription factor evokes vulnerability toward oxidative stress in Ewing sarcoma.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20017-2.

mTOR Signaling and SREBP Activity Increase FADS2 Expression and Can Activate Sapienate Biosynthesis.

Cancer cells display an increased plasticity in their lipid metabolism, which includes the conversion of palmitate to sapienate via the enzyme fatty acid desaturase 2 (FADS2). We find that FADS2 expression correlates with mammalian target of rapamycin (mTOR) signaling and sterol regulatory element-binding protein 1 (SREBP-1) activity across multiple cancer types and is prognostic in some cancer types. Accordingly, activating mTOR signaling by deleting tuberous sclerosis complex 2 (Tsc2) or overexpression of SREBP-1/2 is sufficient to increase FADS2 mRNA expression and sapienate metabolism in mouse embryonic fibroblasts (MEFs) and U87 glioblastoma cells, respectively. Conversely, inhibiting mTOR signaling decreases FADS2 expression and sapienate biosynthesis in MEFs with Tsc2 deletion, HUH7 hepatocellular carcinoma cells, and orthotopic HUH7 liver xenografts. In conclusion, we show that mTOR signaling and SREBP activity are sufficient to activate sapienate metabolism by increasing FADS2 expression. Consequently, targeting mTOR signaling can reduce sapienate metabolism in vivo.

Oncogenic hijacking of a developmental transcription factor evokes vulnerability toward oxidative stress in Ewing sarcoma.

Ewing sarcoma (EwS) is an aggressive childhood cancer likely originating from mesenchymal stem cells or osteo-chondrogenic progenitors. It is characterized by fusion oncoproteins involving EWSR1 and variable members of the ETS-family of transcription factors (in 85% FLI1). EWSR1-FLI1 can induce target genes by using GGAA-microsatellites as enhancers.Here, we show that EWSR1-FLI1 hijacks the developmental transcription factor SOX6 - a physiological driver of proliferation of osteo-chondrogenic progenitors - by binding to an intronic GGAA-microsatellite, which promotes EwS growth in vitro and in vivo. Through integration of transcriptome-profiling, published drug-screening data, and functional in vitro and in vivo experiments including 3D and PDX models, we discover that constitutively high SOX6 expression promotes elevated levels of oxidative stress that create a therapeutic vulnerability toward the oxidative stress-inducing drug Elesclomol.Collectively, our results exemplify how aberrant activation of a developmental transcription factor by a dominant oncogene can promote malignancy, but provide opportunities for targeted therapy.

Pan-Cancer Analysis of Mitochondria Chaperone-Client Co-Expression Reveals Chaperone Functional Partitioning.

Metabolic reprogramming is a hallmark of cancer. Such reprogramming entails the up-regulation of the expression of specific mitochondrial proteins, thus increasing the burden on the mitochondrial protein quality control. However, very little is known about the specificity of interactions between mitochondrial chaperones and their clients, or to what extent the mitochondrial chaperone-client co-expression is coordinated. We hypothesized that a physical interaction between a chaperone and its client in mitochondria ought to be manifested in the co-expression pattern of both transcripts. Using The Cancer Genome Atlas (TCGA) gene expression data from 13 tumor entities, we constructed the mitochondrial chaperone-client co-expression network. We determined that the network is comprised of three distinct modules, each populated with unique chaperone-clients co-expression pairs belonging to distinct functional groups. Surprisingly, chaperonins HSPD1 and HSPE1, which are known to comprise a functional complex, each occupied a different module: HSPD1 co-expressed with tricarboxylic acid cycle cycle enzymes, while HSPE1 co-expressed with proteins involved in oxidative phosphorylation. Importantly, we found that the genes in each module were enriched for discrete transcription factor binding sites, suggesting the mechanism for the coordinated co-expression. We propose that our mitochondrial chaperone-client interactome can facilitate the identification of chaperones supporting specific mitochondrial pathways and bring forth a fundamental principle in metabolic adaptation.

High Specificity of BCL11B and GLG1 for EWSR1-FLI1 and EWSR1-ERG Positive Ewing Sarcoma.

Ewing sarcoma (EwS) is an aggressive cancer displaying an undifferentiated small-round-cell histomorphology that can be easily confused with a broad spectrum of differential diagnoses. Using comparative transcriptomics and immunohistochemistry (IHC), we previously identified BCL11B and GLG1 as potential specific auxiliary IHC markers for EWSR1-FLI1-positive EwS. Herein, we aimed at validating the specificity of both markers in a far larger and independent cohort of EwS (including EWSR1-ERG-positive cases) and differential diagnoses. Furthermore, we evaluated their intra-tumoral expression heterogeneity. Thus, we stained tissue microarrays from 133 molecularly confirmed EwS cases and 320 samples from morphological mimics, as well as a series of patient-derived xenograft (PDX) models for BCL11B, GLG1, and CD99, and systematically assessed the immunoreactivity and optimal cut-offs for each marker. These analyses demonstrated that high BCL11B and/or GLG1 immunoreactivity in CD99-positive cases had a specificity of 97.5% and an accuracy of 87.4% for diagnosing EwS solely by IHC, and that the markers were expressed by EWSR1-ERG-positive EwS. Only little intra-tumoral heterogeneity in immunoreactivity was observed for differential diagnoses. These results indicate that BCL11B and GLG1 may help as specific auxiliary IHC markers in diagnosing EwS in conjunction with CD99, especially if confirmatory molecular diagnostics are not available.

A comparative view on the expression patterns of PD-L1 and PD-1 in soft tissue sarcomas.

Soft tissue sarcomas (STSs) are heterogeneous cancers associated with poor prognosis due to high rates of local recurrence and metastasis. The programmed death receptor ligand 1 (PD-L1) is expressed in several cancers. PD-L1 interacts with its receptor, PD-1, on the surface of tumor-infiltrating lymphocytes (TILs), thereby attenuating anti-cancer immune response. Immune checkpoint inhibitors targeting this interaction have been established as effective anti-cancer drugs. However, studies on the PD-L1 and PD-1 expression status in STS are commonly limited by small sample size, analysis of single STS subtypes, or lack of combinatorial marker assessment. To overcome these limitations, we evaluated the expression patterns of intratumoral PD-L1, the number of TILs, their PD-1 expression, and associations with clinicopathological parameters in a large and comprehensive cohort of 225 samples comprising six STS subtypes. We found that nearly all STS subtypes showed PD-L1 expression on the tumor cells, albeit with a broad range of positivity across subtypes (50% angiosarcomas to 3% synovial sarcomas). Co-expression and correlation analyses uncovered that PD-L1 expression was associated with more PD-1-positive TILs (P < 0.001), higher tumor grading (P = 0.016), and worse patients' 5-year overall survival (P = 0.028). The results were in line with several publications on single STS subtypes, especially when comparing findings for STS with low and high mutational burden. In sum, the substantial portion of PD-L1 positivity, the co-occurrence of PD-1-positive TILs, and the association of PD-L1 with unfavorable clinical outcome provide rationales for immune checkpoint inhibition in patients with PD-L1-positive STS.

Knockout of LASP1 in CXCR4 expressing CML cells promotes cell persistence, proliferation and TKI resistance.

Chronic myeloid leukaemia (CML) is a clonal myeloproliferative stem cell disorder characterized by the constitutively active BCR-ABL tyrosine kinase. The LIM and SH3 domain protein 1 (LASP1) has recently been identified as a novel BCR-ABL substrate and is associated with proliferation, migration, tumorigenesis and chemoresistance in several cancers. Furthermore, LASP1 was shown to bind to the chemokine receptor 4 (CXCR4), thought to be involved in mechanisms of relapse. In order to identify potential LASP1-mediated pathways and related factors that may help to further eradicate minimal residual disease (MRD), the effect of LASP1 on processes involved in progression and maintenance of CML was investigated. The present data indicate that not only overexpression of CXCR4, but also knockout of LASP1 contributes to proliferation, reduced apoptosis and migration as well as increased adhesive potential of K562 CML cells. Furthermore, LASP1 depletion in K562 CML cells leads to decreased cytokine release and reduced NK cell-mediated cytotoxicity towards CML cells. Taken together, these results indicate that in CML, reduced levels of LASP1 alone and in combination with high CXCR4 expression may contribute to TKI resistance.

Focal adhesion kinase confers pro-migratory and antiapoptotic properties and is a potential therapeutic target in Ewing sarcoma.

Oncogenesis of Ewing sarcoma (EwS), the second most common malignant bone tumor of childhood and adolescence, is dependent on the expression of chimeric EWSR1-ETS fusion oncogenes, most often EWSR1-FLI1 (E/F). E/F expression leads to dysregulation of focal adhesions (FAs) enhancing the migratory capacity of EwS cells. Here, we show that, in EwS cell lines and tissue samples, focal adhesion kinase (FAK) is expressed and phosphorylated at Y397 in an E/F-dependent way involving Ezrin. Employing different EwS cell lines as in vitro models, we found that key malignant properties of E/F are mediated via substrate-independent autophosphorylation of FAK on Y397. This phosphorylation results in enhanced FA formation, Rho-dependent cell migration, and impaired caspase-3-mediated apoptosis in vitro. Conversely, treatment with the FAK inhibitor 15 (1,2,4,5-benzenetetraamine tetrahydrochloride (Y15) enhanced caspase-mediated apoptosis and EwS cell migration, independent from the respective EWSR1-ETS fusion type, mimicking an anoikis-like phenotype and paralleling the effects of FAK siRNA knockdown. Our findings were confirmed in vivo using an avian chorioallantoic membrane model and provide a first rationale for the therapeutic use of FAK inhibitors to impair metastatic dissemination of EwS.

Integrative clinical transcriptome analysis reveals TMPRSS2-ERG dependency of prognostic biomarkers in prostate adenocarcinoma.

In prostate adenocarcinoma (PCa), distinction between indolent and aggressive disease is challenging. Around 50% of PCa are characterized by TMPRSS2-ERG (T2E)-fusion oncoproteins defining two molecular subtypes (T2E-positive/negative). However, current prognostic tests do not differ between both molecular subtypes, which might affect outcome prediction. To investigate gene-signatures associated with metastasis in T2E-positive and T2E-negative PCa independently, we integrated tumor transcriptomes and clinicopathological data of two cohorts (total n = 783), and analyzed metastasis-associated gene-signatures regarding the T2E-status. Here, we show that the prognostic value of biomarkers in PCa critically depends on the T2E-status. Using gene-set enrichment analyses, we uncovered that metastatic T2E-positive and T2E-negative PCa are characterized by distinct gene-signatures. In addition, by testing genes shared by several functional gene-signatures for their association with event-free survival in a validation cohort (n = 272), we identified five genes (ASPN, BGN, COL1A1, RRM2 and TYMS)-three of which are included in commercially available prognostic tests-whose high expression was significantly associated with worse outcome exclusively in T2E-negative PCa. Among these genes, RRM2 and TYMS were validated by immunohistochemistry in another validation cohort (n = 135), and several of them proved to add prognostic information to current clinicopathological predictors, such as Gleason score, exclusively for T2E-negative patients. No prognostic biomarkers were identified exclusively for T2E-positive tumors. Collectively, our study discovers that the T2E-status, which is per se not a strong prognostic biomarker, crucially determines the prognostic value of other biomarkers. Our data suggest that the molecular subtype needs to be considered when applying prognostic biomarkers for outcome prediction in PCa.

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.

Gene expression and immunohistochemical analyses identify SOX2 as major risk factor for overall survival and relapse in Ewing sarcoma patients.

BACKGROUND: Up to 30-40% of Ewing sarcoma (EwS) patients with non-metastatic disease develop local or metastatic relapse within a time span of 2-10 years. This is in part caused by the absence of prognostic biomarkers that can identify high-risk patients and thus assign them to risk-adapted monitoring and treatment regimens. Since cancer stemness has been associated with tumour relapse and poor patient outcomes, we investigated in the current study the prognostic potential SOX2 (sex determining region Y box 2) - a major transcription factor involved in development and stemness - which was previously described to contribute to the undifferentiated phenotype of EwS. METHODS: Two independent patient cohorts, one consisting of 189 retrospectively collected EwS tumours with corresponding mRNA expression data (test-cohort) and the other consisting of 141 prospectively collected formalin-fixed and paraffin-embedded resected tumours (validation and cohort), were employed to analyse SOX2 expression levels through DNA microarrays or immunohistochemistry, respectively, and to compare them with clinical parameters and patient outcomes. Two methods were employed to test the validity of the results at both the mRNA and protein levels. FINDINGS: Both cohorts showed that only a subset of EwS patients (16-20%) expressed high SOX2 mRNA or protein levels, which significantly correlated with poor overall survival. Multivariate analyses of our validation-cohort revealed that high SOX2 expression represents a major risk-factor for poor survival (HR = 3·19; 95%CI 1·74-5·84; p < 0·01) that is independent from metastasis and other known clinical risk-factors at the time of diagnosis. Univariate analyses demonstrated that SOX2-high expression was correlated with tumour relapse (p = 0·002). The median first relapse was at 14·7 months (range: 3·5-180·7). INTERPRETATION: High SOX2 expression constitutes an independent prognostic biomarker for EwS patients with poor outcomes. This may help to identify patients with localised disease who are at high risk for tumour relapse within the first two years after diagnosis. FUNDING: The laboratory of T. G. P. Grünewald is supported by grants from the 'Verein zur Förderung von Wissenschaft und Forschung an der Medizinischen Fakultät der LMU München (WiFoMed)', by LMU Munich's Institutional Strategy LMUexcellent within the framework of the German Excellence Initiative, the 'Mehr LEBEN für krebskranke Kinder - Bettina-Bräu-Stiftung', the Walter Schulz Foundation, the Wilhelm Sander-Foundation (2016.167.1), the Friedrich-Baur foundation, the Matthias-Lackas foundation, the Barbara & Hubertus Trettner foundation, the Dr. Leopold & Carmen Ellinger foundation, the Gert & Susanna Mayer foundation, the Deutsche Forschungsgemeinschaft (DFG 391665916), and by the German Cancer Aid (DKH-111886 and DKH-70112257). J. Li was supported by a scholarship of the China Scholarship Council (CSC), J. Musa was supported by a scholarship of the Kind-Philipp foundation, and T. L. B. Hölting by a scholarship of the German Cancer Aid. M. F. Orth and M. M. L. Knott were supported by scholarships of the German National Academic Foundation. G. Sannino was supported by a scholarship from the Fritz-Thyssen Foundation (FTF-40.15.0.030MN). The work of U. Dirksen is supported by grants from the German Cancer Aid (DKH-108128, DKH-70112018, and DKH-70113419), the ERA-Net-TRANSCAN consortium (project number 01KT1310), and Euro Ewing Consortium (EEC, project number EU-FP7 602,856), both funded under the European Commission Seventh Framework Program FP7-HEALTH (http://cordis.europa.eu/), the Barbara & Hubertus Trettner foundation, and the Gert & Susanna Mayer foundation. G. Hardiman was supported by grants from the National Science Foundation (SC EPSCoR) and National Institutes of Health (U01-DA045300). The laboratory of J. Alonso was supported by Instituto de Salud Carlos III (PI12/00816; PI16CIII/00026); Asociación Pablo Ugarte (TPY-M 1149/13; TRPV 205/18), ASION (TVP 141/17), Fundación Sonrisa de Alex & Todos somos Iván (TVP 1324/15).

Evidence for an alternative fatty acid desaturation pathway increasing cancer plasticity.

Most tumours have an aberrantly activated lipid metabolism1,2 that enables them to synthesize, elongate and desaturate fatty acids to support proliferation. However, only particular subsets of cancer cells are sensitive to approaches that target fatty acid metabolism and, in particular, fatty acid desaturation3. This suggests that many cancer cells contain an unexplored plasticity in their fatty acid metabolism. Here we show that some cancer cells can exploit an alternative fatty acid desaturation pathway. We identify various cancer cell lines, mouse hepatocellular carcinomas, and primary human liver and lung carcinomas that desaturate palmitate to the unusual fatty acid sapienate to support membrane biosynthesis during proliferation. Accordingly, we found that sapienate biosynthesis enables cancer cells to bypass the known fatty acid desaturation pathway that is dependent on stearoyl-CoA desaturase. Thus, only by targeting both desaturation pathways is the in vitro and in vivo proliferation of cancer cells that synthesize sapienate impaired. Our discovery explains metabolic plasticity in fatty acid desaturation and constitutes an unexplored metabolic rewiring in cancers.