- Letter to the Editor
- Open Access
Integrative gene network and functional analyses identify a prognostically relevant key regulator of metastasis in Ewing sarcoma
Molecular Cancer volume 21, Article number: 1 (2022)
The identification of cancer stemness genes is crucial to understand the underlying biology of therapy resistance, relapse, and metastasis. In pediatric tumors – despite being largely composed of undifferentiated stem cell-like cells – mutations generally do not involve canonical stemness/metastasis-associated genes . Ewing sarcoma (EwS), the second most common bone cancer in children and adolescents , is a highly aggressive malignancy associated with dismal outcome in metastatic disease , wherefore deciphering mechanisms of metastasis is imperative. EwS is characterized by a remarkably ‘silent’ genome with a single driver mutation generating an oncogenic fusion transcription factor (EWSR1-ETS) that promotes stem cell features [4, 5]. Thus, EwS constitutes an ideal model to study how perturbation of a transcriptional network by a dominant oncogene can mediate metastasis.
To identify highly relevant EWSR1-ETS target genes involved in the stemness/metastasis axis we implemented a systems biology approach to analyze published transcriptome profiles of 5 EwS cell lines with/without shRNA-mediated knockdown of EWSR1-FLI1 or -ERG (< 20%) for 96 h  as described in Fig. 1a. This analysis yielded 348 differentially expressed genes (DEGs) being up- or downregulated (|log2 FC| ≥ 1.5) after knockdown of EWSR1-FLI1 or -ERG across all cell lines (Supplementary Table 1). We next filtered these DEGs for genes annotated with the significantly enriched gene ontology term ‘Regulation of Cell Differentiation’ using PantherDB (P = 4.93 × 10− 12, false-discovery rate = 1.97 × 10− 9). The resulting 76 DEGs (Supplementary Table 2) were subjected to network analysis for pathway, physical, and genetic interactions using Cytoscape and GeneMania (Fig. 1b, c). To identify clinically relevant key hubs within this network we focused on transcription factors (n = 11) that were highly interconnected (> 10 interactions) with a significant (P < 0.05, Mantel-Haenszel test) association with overall survival in a cohort of 166 EwS patients with matched gene expression and clinical data . After adjustment for multiple testing (Supplementary Table 3), transcription factor 7 like 1 (TCF7L1, alias TCF3) emerged as the most promising candidate for functional follow-up, whose low expression was associated with poor patients’ overall survival (nominal P = 0.0057; P = 0.0228 Bonferroni-adjusted) (Fig. 1d). Next, we explored TCF7L1 protein levels in association with clinical outcome in tissue micro-arrays (TMAs) comprising 114 EwS tumors and found that low TCF7L1 protein expression was significantly (P = 0.035) associated with worse overall survival (Fig. 1e), supporting that low TCF7L1 expression promotes an aggressive phenotype. Notably, TCF7L1 appeared to be downregulated by EWSR1-ETS since microarray data of each respective fusion oncogene knockdown induced a ~3-fold upregulation of TCF7L1 (Supplementary Table 4), which was confirmed in 5 EwS cell lines by qRT-PCR (average ~ 6-fold upregulation) (Supplementary Fig. 1b) and validated in A673 cells xenografts with/without conditional knockdown of EWSR1-FLI1 (A673/TR/shEF1 cells) at mRNA and protein levels (Supplementary Fig. 1c).
Prior reports suggested that TCF7L1 may have either oncogenic or tumor suppressor properties depending on the cellular context. Indeed, TCF7L1 expression has been linked to promotion of cell proliferation, tumor growth and sphere formation in breast cancer , colorectal cancer [9, 10], acute lymphoblastic leukemia  and skin squamous cell carcinoma , while ectopic expression of TCF7L1 inhibits self-renewal of liver cancer stem cells . To explore its role in EwS, we first analyzed the TCF7L1 expression pattern across 18 cancer entities using microarray data from the Cancer Cell Line Encyclopedia (CCLE) and an own study that compiled well-curated microarray data from the same cancer entities . Surprisingly, TCF7L1 was very highly, but variably, expressed in EwS cell lines and primary tumors (Fig. 1f).
Weighted Gene Correlation Network Analysis (WGCNA) based on enriched gene-sets in TCF7L1-correlated genes in 166 EwS tumors showed that EwS tumors with low TCF7L1 expression were enriched in embryonic pathways (Fig. 1g). Although TCF7L1 is generally highly expressed in EwS, these data suggested that suppression of its transcription by EWSR1-FLI1 is associated with embryonic processes that may be linked with poor patient outcome. Thus, we generated two EwS cell lines (SK-N-MC and TC-71, with lowest TCF7L1 expression across the 5 EwS cell lines tested (Supplementary Fig. 1a) with a doxycycline- (DOX)-inducible re-expression of TCF7L1 (Fig. 1h). Consistent with the hypothesis that TCF7L1 may be an EWSR1-ETS-repressed downstream transcription factor (with a potentially indirect regulation, possibly partially mediated by TCF7L1 promoter hypermethylation as depicted in Supplementary Fig. 1d–g), transcriptome profiling of two EwS cell lines after either knockdown of EWSR1-ETS or re-expression of TCF7L1 showed a highly significant (P = 2.57 × 10− 165 or P = 5.34 × 10− 272) overlap of concordant DEGs (Supplementary Fig. 1e). Congruently, conditional re-expression of TCF7L1 significantly reduced cell proliferation (Fig. 1i), clonogenic growth in two-dimensional (2D) cultures (Supplementary Fig. 2b), anchorage-independent and three-dimensional (3D) growth (Fig. 1j), and local tumor growth in a pre-clinical subcutaneous xenotransplantation mouse model (Fig. 1k). The xenografts exhibiting specific re-expression of TCF7L1 (DOX(+) group) (Supplementary Fig. 2e,f) showed a reduced mitotic index (Fig. 1l), and a significant decrease of the proliferation marker Ki67 (Fig. 1m). Similar in vitro and in vivo experiments performed with both cell lines ‘empty’ controls exhibited no significant differences (Supplementary Fig. 2a, c-d, g-i).
Since stemness features, such as elevated clonogenic capacity and anchorage-independent growth, are essential for circulating tumor cells to colonize and develop into clinically apparent metastases in distant organs, we reasoned that TCF7L1 may be linked to the metastatic process in EwS. In support of this hypothesis, transcriptome profiling, subsequent gene-set enrichment and WGCNA of TCF7L1-reexpressing SK-N-MC and TC-71 EwS cells uncovered that low TCF7L1 levels led to overrepresentation of gene-sets involved in cellular migration (Fig. 2a). Accordingly, re-expression of TCF7L1 resulted in significantly reduced transwell migration (Fig. 2b), invasion, and single-cell 3D migration in fibrin gel using advanced microfluidic chambers (Fig. 2c). Next, we employed an orthotopic spontaneous in vivo metastasis model by injecting TCF7L1 re-expressing cells into the proximal tibiae of immunocompromised mice. Remarkably, while there was no significant difference in local tumor growth in the limited space of the tibial plateau (Supplementary Fig. 3a,b), we observed a strong inhibition of macrometastatic spread to liver, lungs, and kidneys upon re-expression of TCF7L1, and significantly reduced micrometastatic burden in the same organs (Fig. 2d,e, Supplementary Fig. 3c–e). In agreement with these in vivo findings, we observed that TCF7L1 was significantly (P = 0.016) lower expressed in EwS metastases (n = 7) compared to primary tumors (n = 118) in situ as defined by RNA-sequencing of tumor tissue (Fig. 2f). Strikingly, this result became even more significant when focusing on gene expression data of 4 matched pairs of metastases and primaries (P = 0.0078) (Fig. 2g).
To further investigate the mechanism of action of TCF7L1, we generated two EwS cell lines with conditional re-expression of different TCF7L1 deletion mutants for its major domains (ΔHMG: DNA binding domain; ΔCTNNB: β-catenin binding domain) and observed that only the ΔHMG mutants exhibited a normal clonogenic growth and migration capacity in vitro (Fig. 2h,i), and tumor growth in vivo (Fig. 2j). Moreover, further analysis of our WCGNA on TCF7L1-regulated signatures displayed in Fig. 2a highlighted several downregulated ‘driver target genes’ (aka leading-edge genes), including ANXA1, LMO7, SLC9A9, and TMEM71 as potential key mediators of TCF7L1 inhibition of the EwS migratory phenotype. These genes were selected for validation as they ranged among the top 10 downregulated genes by TCF7L1 and as their upregulation has been previously implicated in cancer progression [15,16,17,18]. As shown in Supplementary Fig. 3f, these genes are strongly downregulated upon re-expression of wildtype TCF7L1 in both EwS cell lines tested, but remain unaltered in the DNA-binding-domain deletion mutant (ΔHMG), further suggesting that DNA-binding of TCF7L1 is required to suppress the tumorigenic and migratory phenotype of EwS cells.
Previous studies proposed a binary model of EWSR1-FLI1-high promoting proliferation and EWSR1-FLI1-low promoting metastasis  following a traditional epithelial-mesenchymal-transition concept. However, our findings support an integrative concept of the ‘migratory stem cell’ : EwS cells may reside in a ‘metastable’ state where all EwS cells of a tumor may be equipped with both proliferative and migratory capacities . Upon different intrinsic/extrinsic cues (of which EWSR1-FLI1 may not be the only important component) these cells could underlie a certain phenotypic accentuation.
In this context, we demonstrated that TCF7L1 is a critical mediator of metastasis in EwS, which may be utilized as a prognostic biomarker. Since we could detect a strong inverse correlation of TCF7L1 levels with patients’ overall survival in both the mRNA- as well as in the TMA-cohort, we believe that the most straightforward possibility to translate TCF7L1 into a routine clinical setting would be the IHC detection and semiquantitative evaluation of TCF7L1 in primary biopsies. This technique is readily available, inexpensive, and the conditions of IHC staining of TCF7L1 in EwS tumors have been established in the current study, which would offer the possibility of further evaluating the prognostic value of TCF7L1 in ongoing and prospective clinical trials. In addition, it would be highly interesting to explore whether TCF7L1 may serve as a prognostic biomarker in other (EWSR1-rearranged) mesenchymal neoplasms, which is subject to future studies.
In summary, our study exemplifies the power of systems biology to decipher gene regulatory networks and to identify key players in the metastatic process, which may be highly relevant for, and translatable to, other oligomutated (pediatric) cancers.
Availability of data and materials
Original microarray data that support the findings of this study were deposited at the National Center for Biotechnology Information (NCBI) GEO and are accessible through the accession number GSE165929. Custom code is available from the corresponding author upon reasonable request.
Cancer Cell Line Encyclopedia
Differentially expressed genes
Epithelial to mesenchymal transition
E26 transformation specific
Ewing sarcoma breakpoint region 1
ETS transcription factor ERG
False discovery rate
FUS, EWSR1, TAF15
Formalin-fixed and paraffin-embedded
Friend leukemia virus integration 1
Gene set enrichment analysis
High mobility group
Immune reactive score
3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide
Normalized enrichment score
T-cell factor 3
Transcription factor 7 like 1
Open reading frame
Weighted gene co-expression network analysis
Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499(7457):214–8.
Grünewald TGP, Cidre-Aranaz F, Surdez D, Tomazou EM, de Álava E, Kovar H, et al. Ewing sarcoma. Nat Rev Dis Primer. 2018 05;4(1):5.
Thiel U, Wawer A, Wolf P, Badoglio M, Santucci A, Klingebiel T, et al. No improvement of survival with reduced- versus high-intensity conditioning for allogeneic stem cell transplants in Ewing tumor patients. Ann Oncol Off J Eur Soc Med Oncol. 2011;22(7):1614–21.
Tirode F, Laud-Duval K, Prieur A, Delorme B, Charbord P, Delattre O. Mesenchymal stem cell features of Ewing tumors. Cancer Cell. 2007;11(5):421–9.
Riggi N, Suvà M-L, De Vito C, Provero P, Stehle J-C, Baumer K, et al. EWS-FLI-1 modulates miRNA145 and SOX2 expression to initiate mesenchymal stem cell reprogramming toward Ewing sarcoma cancer stem cells. Genes Dev. 2010;24(9):916–32.
Orth MF, Surdez D, Marchetto A, Grossetête S, Gerke JS, Zaidi S, et al. Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation [Internet]. 2021 Jun [cited 2021 Aug 16] p. 2021.06.08.447518. Available from: https://www.biorxiv.org/content/10.1101/2021.06.08.447518v1
Musa J, Cidre-Aranaz F, Aynaud M-M, Orth MF, Knott MML, Mirabeau O, et al. Cooperation of cancer drivers with regulatory germline variants shapes clinical outcomes. Nat Commun. 2019 11;10(1):4128.
Slyper M, Shahar A, Bar-Ziv A, Granit RZ, Hamburger T, Maly B, et al. Control of breast cancer growth and initiation by the stem cell-associated transcription factor TCF3. Cancer Res. 2012;72(21):5613–24.
Murphy M, Chatterjee SS, Jain S, Katari M, DasGupta R. TCF7L1 Modulates Colorectal Cancer Growth by Inhibiting Expression of the Tumor-Suppressor Gene EPHB3. Sci Rep. 2016 23;6:28299.
Eshelman MA, Shah M, Raup-Konsavage WM, Rennoll SA, Yochum GS. TCF7L1 recruits CtBP and HDAC1 to repress DICKKOPF4 gene expression in human colorectal cancer cells. Biochem Biophys Res Commun. 2017;487(3):716–22.
Ma H, Mallampati S, Lu Y, Sun B, Wang E, Leng X, et al. The Sox4/Tcf7l1 axis promotes progression of BCR-ABL-positive acute lymphoblastic leukemia. Haematologica. 2014;99(10):1591–8.
Ku AT, Shaver TM, Rao AS, Howard JM, Rodriguez CN, Miao Q, et al. TCF7L1 promotes skin tumorigenesis independently of β-catenin through induction of LCN2. eLife. 2017;03:6.
Shan J, Shen J, Wu M, Zhou H, Feng J, Yao C, et al. Tcf7l1 acts as a suppressor for the self-renewal of liver Cancer stem cells and is regulated by IGF/MEK/ERK signaling independent of β-catenin. Stem Cells Dayt Ohio. 2019;37(11):1389–400.
Baldauf MC, Orth MF, Dallmayer M, Marchetto A, Gerke JS, Rubio RA, et al. Robust diagnosis of Ewing sarcoma by immunohistochemical detection of super-enhancer-driven EWSR1-ETS targets. Oncotarget. 2018;9(2):1587–601.
de Graauw M, van Miltenburg MH, Schmidt MK, Pont C, Lalai R, Kartopawiro J, et al. Annexin A1 regulates TGF-beta signaling and promotes metastasis formation of basal-like breast cancer cells. Proc Natl Acad Sci U S A. 2010;107(14):6340–5.
Liu X, Yuan H, Zhou J, Wang Q, Qi X, Bernal C, et al. LMO7 as an unrecognized factor promoting pancreatic Cancer progression and metastasis. Front Cell Dev Biol. 2021;9:647387.
Ueda M, Iguchi T, Masuda T, Komatsu H, Nambara S, Sakimura S, et al. Up-regulation of SLC9A9 promotes Cancer progression and is involved in poor prognosis in colorectal Cancer. Anticancer Res. 2017 May;37(5):2255–63.
Wang K-Y, Huang R-Y, Tong X-Z, Zhang K-N, Liu Y-W, Zeng F, et al. Molecular and clinical characterization of TMEM71 expression at the transcriptional level in glioma. CNS Neurosci Ther. 2019;25(9):965–75.
Franzetti G-A, Laud-Duval K, van der Ent W, Brisac A, Irondelle M, Aubert S, et al. Cell-to-cell heterogeneity of EWSR1-FLI1 activity determines proliferation/migration choices in Ewing sarcoma cells. Oncogene. 2017;36(25):3505–14.
Brabletz T, Jung A, Spaderna S, Hlubek F, Kirchner T. Opinion: migrating cancer stem cells - an integrated concept of malignant tumour progression. Nat Rev Cancer. 2005;5(9):744–9.
Sannino G, Marchetto A, Kirchner T, Grünewald TGP. Epithelial-to-mesenchymal and mesenchymal-to-epithelial transition in mesenchymal tumors: a paradox in sarcomas? Cancer Res. 2017;77(17):4556–61.
We thank Dr. Paul Northcott and Dr. Brian Gudenas for providing R scripts for weighted gene correlation network analyses and Carlos Rodriguez-Martin for assistance with script adaptation. We thank Susanne Jabar and Dr. Andreas Ranft for statistical advice. We thank Anja Heier, Andrea Sendelhofert, Lea Hofmann and Laura Dörner for expert technical assistance. We thank Dr. Martin Sill for contributing to methylation analyses.
This work was mainly supported by the Wilhelm Sander-Foundation (2016.167.1), the Barbara & Hubertus Trettner Foundation, and the Dr. Rolf M. Schwiete Foundation. In addition, the laboratory of T.G.P.G. was supported by the Matthias-Lackas Foundation, the Dr. Leopold and Carmen Ellinger Foundation, the German Cancer Aid (DKH-70112257 and DKH-70114111), the Gert und Susanna Mayer the Boehringer Ingelheim Foundation, the SMARCB1 association, and the Barbara und Wilfried Mohr Foundation. J.L. was supported by a scholarship of the Chinese Scholarship Council (CSC), and M.K. and C.M.F. by scholarships from the German Cancer Aid. M.M. was supported by the Swiss National Science Foundation (SNF 310030_179167). U.T. was supported by the Wilhelm Sander-Foundation (2021.085.1), the Dr. Sepp and Hanne Sturm Foundation and by Zukunft Gesundheit e.V.
Ethics approval and consent to participate
Human tissue samples were retrieved from the tissue archives of the Institute of Pathology of the LMU Munich (Germany), the Bone Tumor Reference Center at the Institute of Pathology of the University Hospital of Basel, or the Gerhard-Domagk Institute of Pathology of the University of Münster (Germany) upon approval of the institutional review board. All patients provided informed consent. Tissue-microarrays (TMAs) were stained and analyzed with approval of the ethics committee of the LMU Munich (approval no. 550–16 UE). Animal experiments were approved by the government of Upper Bavaria and Northbaden and performed in accordance with ARRIVE guidelines, recommendations of the European Community (86/609/EEC), and United Kingdom Coordinating Committee on Cancer Research (UKCCCR) guidelines for the welfare and use of animals in cancer research.
Consent for publication
T.G.P.G. serves as honorary consultant for Boehringer-Ingelheim International GmbH. All other authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cidre-Aranaz, F., Li, J., Hölting, T.L.B. et al. Integrative gene network and functional analyses identify a prognostically relevant key regulator of metastasis in Ewing sarcoma. Mol Cancer 21, 1 (2022). https://doi.org/10.1186/s12943-021-01470-z