Search

Rhabdomyosarcoma (RMS) is a pediatric tumor that resembles undifferentiated muscle cells; yet the extent to which cell state heterogeneity is shared with human development has not been described. Using single-cell/nucleus RNA sequencing from patient tumors, patient-derived xenografts, primary in vitro cultures, and cell lines, we identify four dominant muscle-lineage cell states: progenitor, proliferative, differentiated, and ground cells. We stratify these RMS cells/nuclei along the continuum of human muscle development and show that they share expression patterns with fetal/embryonal myogenic precursors rather than postnatal satellite cells. Fusion-negative RMS (FN-RMS) have a discrete stem cell hierarchy that recapitulates fetal muscle development and contain therapy-resistant FN-RMS progenitors that share transcriptomic similarity with bipotent skeletal mesenchymal cells. Fusion-positive RMS have tumor-acquired cells states, including a neuronal cell state, that are not found in myogenic development. This work identifies previously underappreciated cell state heterogeneity including unique treatment-resistant and tumor-acquired cell states that differ across RMS subtypes., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Rhabdomyosarcoma (RMS) is a highly aggressive pediatric cancer with features of skeletal muscle differentiation. More than 80% of the high-risk patients ultimately fail to respond to chemotherapy treatment, leading to limited therapeutic options and dismal prognostic rates. The lack of response and subsequent tumor recurrence is driven in part by stem cell-like cells, the tumor subpopulation that is enriched after treatment, and characterized by expression of the AXL receptor tyrosine kinase (AXL). AXL mediates survival, migration, and therapy resistance in several cancer types; however, its function in RMS remains unclear. In this study, we investigated the role of AXL in RMS tumorigenesis, migration, and chemotherapy response, and whether targeting of AXL with small-molecule inhibitors could potentiate the efficacy of chemotherapy. We show that AXL is expressed in a heterogeneous manner in patient-derived xenografts (PDX), primary cultures and cell line models of RMS, consistent with its stem cell-state selectivity. By generating a CRISPR/Cas9 AXL knock-out and overexpressing models, we show that AXL contributes to the migratory phenotype of RMS, but not to chemotherapy resistance. Instead, pharmacologic blockade with the AXL inhibitors bemcentinib (BGB324), cabozantinib and NPS-1034 rapidly killed RMS cells in an AXL-independent manner and augmented the efficacy of the chemotherapeutics vincristine and cyclophosphamide. In vivo administration of the combination of bemcentinib and vincristine exerted strong antitumoral activity in a rapidly progressing PDX mouse model, significantly reducing tumor burden compared with single-agent treatment. Collectively, our data identify bemcentinib as a promising drug to improve chemotherapy efficacy in patients with RMS., (©2024 The Authors; Published by the American Association for Cancer Research.)

Activation of oncogenic gene expression from long-range enhancers is initiated by the assembly of DNA-binding transcription factors (TF), leading to recruitment of co-activators such as CBP/p300 to modify the local genomic context and facilitate RNA-Polymerase 2 (Pol2) binding. Yet, most TF-to-coactivator recruitment relationships remain unmapped. Here, studying the oncogenic fusion TF PAX3-FOXO1 (P3F) from alveolar rhabdomyosarcoma (aRMS), we show that a single cysteine in the activation domain (AD) of P3F is important for a small alpha helical coil that recruits CBP/p300 to chromatin. P3F driven transcription requires both this single cysteine and CBP/p300. Mutants of the cysteine reduce aRMS cell proliferation and induce cellular differentiation. Furthermore, we discover a profound dependence on CBP/p300 for clustering of Pol2 loops that connect P3F to its target genes. In the absence of CBP/p300, Pol2 long range enhancer loops collapse, Pol2 accumulates in CpG islands and fails to exit the gene body. These results reveal a potential novel axis for therapeutic interference with P3F in aRMS and clarify the molecular relationship of P3F and CBP/p300 in sustaining active Pol2 clusters essential for oncogenic transcription., (© 2023. The Author(s).)

Rhabdomyosarcoma (RMS), the most common soft-tissue sarcoma in children and adolescents, represents an aberrant form of skeletal muscle differentiation. Both skeletal muscle development, as well as regeneration of adult skeletal muscle are governed by members of the myogenic family of regulatory transcription factors (MRFs), which are deployed in a highly controlled, multi-step, bidirectional process. Many aspects of this complex process are deregulated in RMS and contribute to tumorigenesis. Interconnected loops of super-enhancers, called core regulatory circuitries (CRCs), define aberrant muscle differentiation in RMS cells. The transcriptional regulation of MRF expression/activity takes a central role in the CRCs active in skeletal muscle and RMS. In PAX3::FOXO1 fusion-positive (PF+) RMS, CRCs maintain expression of the disease-driving fusion oncogene. Recent single-cell studies have revealed hierarchically organized subsets of cells within the RMS cell pool, which recapitulate developmental myogenesis and appear to drive malignancy. There is a large interest in exploiting the causes of aberrant muscle development in RMS to allow for terminal differentiation as a therapeutic strategy, for example, by interrupting MEK/ERK signaling or by interfering with the epigenetic machinery controlling CRCs. In this review, we provide an overview of the genetic and epigenetic framework of abnormal muscle differentiation in RMS, as it provides insights into fundamental mechanisms of RMS malignancy, its remarkable phenotypic diversity and, ultimately, opportunities for therapeutic intervention.

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.

Purpose: Rhabdomyosarcomas (RMS) are rare neoplasms affecting children and young adults. Efforts to improve patient survival have been undermined by a lack of suitable disease markers. Plasma circulating tumor DNA (ctDNA) has shown promise as a potential minimally invasive biomarker and monitoring tool in other cancers; however, it remains underexplored in RMS. We aimed to determine the feasibility of identifying and quantifying ctDNA in plasma as a marker of disease burden and/or treatment response using blood samples from RMS mouse models and patients., Methods: We established mouse models of RMS and applied quantitative polymerase chain reaction (PCR) and droplet digital PCR (ddPCR) to detect ctDNA within the mouse plasma. Potential driver mutations, copy-number alterations, and DNA breakpoints associated with PAX3 / 7-FOXO1 gene fusions were identified in the RMS samples collected at diagnosis. Patient-matched plasma samples collected from 28 patients with RMS before, during, and after treatment were analyzed for the presence of ctDNA via ddPCR, panel sequencing, and/or whole-exome sequencing., Results: Human tumor-derived DNA was detectable in plasma samples from mouse models of RMS and correlated with tumor burden. In patients, ctDNA was detected in 14/18 pretreatment plasma samples with ddPCR and 7/7 cases assessed by sequencing. Levels of ctDNA at diagnosis were significantly higher in patients with unfavorable tumor sites, positive nodal status, and metastasis. In patients with serial plasma samples (n = 18), fluctuations in ctDNA levels corresponded to treatment response., Conclusion: Comprehensive ctDNA analysis combining high sensitivity and throughput can identify key molecular drivers in RMS models and patients, suggesting potential as a minimally invasive biomarker. Preclinical assessment of treatments using mouse models and further patient testing through prospective clinical trials are now warranted.