Rhabdomyosarcoma (RMS) is a group of pediatric cancers with features of developing skeletal muscle. The cellular hierarchy and mechanisms leading to developmental arrest remain elusive. Here, we combined single-cell RNA sequencing, mass cytometry, and high-content imaging to resolve intratumoral heterogeneity of patient-derived primary RMS cultures. We show that the aggressive alveolar RMS (aRMS) subtype contains plastic muscle stem-like cells and cycling progenitors that drive tumor growth, and a subpopulation of differentiated cells that lost its proliferative potential and correlates with better outcomes. While chemotherapy eliminates cycling progenitors, it enriches aRMS for muscle stem-like cells. We screened for drugs hijacking aRMS toward clinically favorable subpopulations and identified a combination of RAF and MEK inhibitors that potently induces myogenic differentiation and inhibits tumor growth. Overall, our work provides insights into the developmental states underlying aRMS aggressiveness, chemoresistance, and progression and identifies the RAS pathway as a promising therapeutic target.

The majority of current cancer therapies are aimed at reducing tumour growth, but there is lack of viable pharmacological options to reduce the formation of metastasis. This is a paradox, since more than 90% of cancer deaths are attributable to metastatic progression. Integrin alpha9 (ITGA9) has been previously described as playing an essential role in metastasis; however, little is known about the mechanism that links this protein to this process, being one of the less studied integrins. We have now deciphered the importance of ITGA9 in metastasis and provide evidence demonstrating its essentiality for metastatic dissemination in rhabdomyosarcoma and neuroblastoma. However, the most translational advance of this study is to reveal, for the first time, the possibility of reducing metastasis by pharmacological inhibition of ITGA9 with a synthetic peptide simulating a key interaction domain of ADAM proteins, in experimental metastasis models, not only in childhood cancers but also in a breast cancer model.

Rhabdomyosarcoma (RMS) is the most common type of soft tissue sarcoma in children and can be divided into two main subtypes: embryonal (eRMS) and alveolar (aRMS). Among the cellular heterogeneity of tumors, the existence of a small fraction of cells called cancer stem cells (CSC), thought to be responsible for the onset and propagation of cancer, has been demonstrated in some neoplasia. Although the existence of CSC has been reported for eRMS, their existence in aRMS, the most malignant subtype, has not been demonstrated to date. Given the lack of suitable markers to identify this subpopulation in aRMS, we used cancer stem cell-enriched supracellular structures (spheres and holoclones) to study this subpopulation. This strategy allowed us to demonstrate the capacity of both aRMS and eRMS cells to form these structures and retain self-renewal capacity. Furthermore, cells contained in spheres and holoclones showed significant Hedgehog pathway induction, the inhibition of which (pharmacologic or genetic) impairs the formation of both holoclones and spheres. Our findings point to a crucial role of this pathway in the maintenance of these structures and suggest that Hedgehog pathway targeting in CSC may have great potential in preventing local relapses and metastases.

The identification of novel therapeutic targets for specific cancer molecular subtypes is crucial for the development of precision oncology. In the last few years, CRISPR/Cas9 screens have accelerated the discovery and validation of new targets associated with different tumor types, mutations, and fusions. However, there are still many cancer vulnerabilities associated with specific molecular features that remain to be explored. Here, we used data from CRISPR/Cas9 screens in 954 cancer cell lines to identify gene dependencies associated with 16 common cancer genomic amplifications. We found that high-copy-number genomic amplifications generate multiple collateral dependencies within the amplified region in most cases. Further, to prioritize candidate targets for each chromosomal region amplified, we integrated gene dependency parameters with both druggability data and subcellular location. Finally, analysis of the relationship between gene expression and gene dependency led to the identification of genes, the expression of which may constitute predictive biomarkers of dependency. In conclusion, our study provides a set of druggable targets specific for each amplification, opening the possibility to specifically target amplified tumors on this basis.