No abstract

Background: Ewing sarcoma (ES) cells exhibit extreme plasticity that contributes to the cell’s survival and recurrence. Although multiple studies reveal various signaling pathways mediated by the EWSR1/FLI1 fusion, the specific transcriptional control of tumor cell resistance to doxorubicin is unknown. Understanding the molecular hubs that contribute to this behavior provides a new perspective on valuable therapeutic options against tumor cells. Methods: Single-cell RNA sequencing and LC-MS/MS-based quantitative proteomics were used. Results: A goal of this study was to identify protein hubs that would help elucidate tumor resistance which prompted ES to relapse or metastasize. Several differentially expressed genes and proteins, including adhesion, cytoskeletal, and signaling molecules, were observed between embryonic fibroblasts and control and doxorubicin-treated tumor cell lines. While several cancer-associated genes/proteins exhibited similar expression across fibroblasts and non-treated cells, upregulation of some proteins belonging to metabolic, stress response, and growth pathway activation was uniquely observed in doxorubicin-treated sarcoma cells, respectively. The novel information on differentially expressed genes/proteins provides insights into the biology of ES cells, which could help elucidate mechanisms of their recurrence. Conclusions: Collectively, our results identify a novel role of cellular proteins in contributing to tumor cell resistance and escape from doxorubicin therapy and contributing to ES progression.

The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.