Ewing sarcoma is a challenging cancer entity, which, besides the characteristic presence of a fusion gene, is driven by multiple alternative splicing events. So far, splice variants in Ewing sarcoma cells were mainly analyzed for EWSR1‑FLI1. The present study provided a comprehensive alternative splicing study on CADO‑ES1, an Ewing model cell line for an EWSR1‑ERG fusion gene. Based on a well‑-characterized RNA‑sequencing dataset with extensive control mechanisms across all levels of analysis, the differential spliced genes in Ewing cancer stem cells were ATP13A3 and EPB41, while the main population was defined by ACADVL, NOP58 and TSPAN3. All alternatively spliced genes were further characterized by their Gene Ontology (GO) terms and by their membership in known protein complexes. These results confirm and extend previous studies towards a systematic whole‑transcriptome analysis. A highlight is the striking segregation of GO terms associated with five basic splice events. This mechanistic insight, together with a coherent integration of all observations with prior knowledge, indicates that EWSR1‑ERG is truly a close twin to EWSR1‑FLI1, but still exhibits certain individuality. Thus, the present study provided a measure of variability in Ewing sarcoma, whose understanding is essential both for clinical procedures and basic mechanistic insight.
Keywords: CADO‑ES1; EWSR1‑ERG; Ewing sarcoma; alternative splicing; cancer stem cell; rMATS.