ABSTRACT
BACKGROUND: Ewing sarcoma is a malignant bone tumor mainly affecting teenagers and young adults. Its main driver mutation, the EWS-FLI1 fusion gene, has been identified more than 20 years ago, whereas its other somatic mutations have been just recently reported. METHODS: In this study, we organized the somatic mutations from 216 Ewing sarcoma cases into 216 individual protein-protein interaction networks by using interactome information. These mutation networks were then classified into five different clusters based on their structural similarities. The prognostic effect of mutation genes was evaluated according to their network features. RESULTS: The cases in cluster two exhibited remarkably high metastasis and mortality rates, and STAG2, TP53 and TTN were the three most significantly mutated genes in this cluster. Microarray data demonstrate that the expression of STAG2, TP53 and TTN are down-regulated in the EWS-FLI1-knockdown Ewing sarcoma cells. However, the mutation effect analysis shows that the somatic mutations in TTN are less damaging than those in STAG2 and TP53. The analyses of functional network modules further revealed that STAG2, TP53 and their interacting gene partners participate in the oncogenic-related biological modules such as cell cycle and regulation of transcription from RNA polymerase II promoter while TTN, TP53 and their interacting gene partners constitute the modules less relevant to oncogenesis. The results of Gene Ontology analyses demonstrated that the uniquely mutated genes associated with poor prognosis in Clusters 1, 4 and 5 were distinctively enriched in epidermal growth factor-related functions and phosphoproteins. CONCLUSIONS: Our study identified the highly lethal mutation combination cases and characterized the possible prognostic genes in Ewing sarcoma from a network perceptive.