Unified rhombic lip origins of group 3 and group 4 medulloblastoma.

Medulloblastoma, a malignant childhood cerebellar tumour, segregates molecularly into biologically distinct subgroups, suggesting that a personalized approach to therapy would be beneficial1. Mouse modelling and cross-species genomics have provided increasing evidence of discrete, subgroup-specific developmental origins2. However, the anatomical and cellular complexity of developing human tissues3-particularly within the rhombic lip germinal zone, which produces all glutamatergic neuronal lineages before internalization into the cerebellar nodulus-makes it difficult to validate previous inferences that were derived from studies in mice. Here we use multi-omics to resolve the origins of medulloblastoma subgroups in the developing human cerebellum. Molecular signatures encoded within a human rhombic-lip-derived lineage trajectory aligned with photoreceptor and unipolar brush cell expression profiles that are maintained in group 3 and group 4 medulloblastoma, suggesting a convergent basis. A systematic diagnostic-imaging review of a prospective institutional cohort localized the putative anatomical origins of group 3 and group 4 tumours to the nodulus. Our results connect the molecular and phenotypic features of clinically challenging medulloblastoma subgroups to their unified beginnings in the rhombic lip in the early stages of human development.

Loss of miR-144/451 alleviates ß-thalassemia by stimulating ULK1-mediated autophagy of free α-globin.

Most cells can eliminate unstable or misfolded proteins through quality control mechanisms. In the inherited red blood cell disorder ß-thalassemia, mutations in the ß-globin gene (HBB) lead to a reduction in the corresponding protein and the accumulation of cytotoxic free α-globin, which causes maturation arrest and apoptosis of erythroid precursors and reductions in the lifespan of circulating red blood cells. We showed previously that excess α-globin is eliminated by Unc-51-like autophagy activating kinase 1 (ULK1)-dependent autophagy and that stimulating this pathway by systemic mammalian target of rapamycin complex 1 (mTORC1) inhibition alleviates ß-thalassemia pathologies. We show here that disrupting the bicistronic microRNA gene miR-144/451 alleviates ß-thalassemia by reducing mTORC1 activity and stimulating ULK1-mediated autophagy of free α-globin through 2 mechanisms. Loss of miR-451 upregulated its target messenger RNA, Cab39, which encodes a cofactor for LKB1, a serine-threonine kinase that phosphorylates and activates the central metabolic sensor adenosine monophosphate-activated protein kinase (AMPK). The resultant enhancement of LKB1 activity stimulated AMPK and its downstream effects, including repression of mTORC1 and direct activation of ULK1. In addition, loss of miR-144/451 inhibited the expression of erythroblast transferrin receptor 1, causing intracellular iron restriction, which has been shown to inhibit mTORC1, reduce free α-globin precipitates, and improve hematological indices in ß-thalassemia. The beneficial effects of miR-144/451 loss in ß-thalassemia were inhibited by the disruption of Cab39 or Ulk1 genes. Together, our findings link the severity of ß-thalassemia to a highly expressed erythroid microRNA locus and a fundamental, metabolically regulated protein quality control pathway that is amenable to therapeutic manipulation.

AAV-mediated gene therapy for sialidosis.

Sialidosis (mucolipidosis I) is a glycoprotein storage disease, clinically characterized by a spectrum of systemic and neurological phenotypes. The primary cause of the disease is deficiency of the lysosomal sialidase NEU1, resulting in accumulation of sialylated glycoproteins/oligosaccharides in tissues and body fluids. Neu1-/- mice recapitulate the severe, early-onset forms of the disease, affecting visceral organs, muscles, and the nervous system, with widespread lysosomal vacuolization evident in most cell types. Sialidosis is considered an orphan disorder with no therapy currently available. Here, we assessed the therapeutic potential of AAV-mediated gene therapy for the treatment of sialidosis. Neu1-/- mice were co-injected with two scAAV2/8 vectors, expressing human NEU1 and its chaperone PPCA. Treated mice were phenotypically indistinguishable from their WT controls. NEU1 activity was restored to different extent in most tissues, including the brain, heart, muscle, and visceral organs. This resulted in diminished/absent lysosomal vacuolization in multiple cell types and reversal of sialyl-oligosacchariduria. Lastly, normalization of lysosomal exocytosis in the cerebrospinal fluids and serum of treated mice, coupled to diminished neuroinflammation, were measures of therapeutic efficacy. These findings point to AAV-mediated gene therapy as a suitable treatment for sialidosis and possibly other diseases, associated with low NEU1 expression.

Single cell transcriptomic profiling identifies tumor-acquired and therapy-resistant cell states in pediatric rhabdomyosarcoma.

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.

DRD1 suppresses cell proliferation and reduces EGFR activation and PD-L1 expression in NSCLC.

Dopamine (DA) acts in various key neurological and physiological processes as both a neurotransmitter and circulating hormone. Over the past several decades, the DA signaling network has been shown to regulate the progression of several types of solid tumors, and considerable evidence has shown it is a druggable pathway in the cancer cell context. However, the specific activity and effect of these pathway components appears to be tissue-type and cell-context-dependent. In the present study, expression and methylation of dopamine receptor D1 (DRD1) were measured using RNA sequencing (RNAseq) and reverse transcription polymerase chain reaction (RT-PCR) in non-small cell lung cancer (NSCLC) samples, and validated using publicly available datasets, including The Cancer Genome Atlas (TCGA). In vitro and in vivo functional experiments were performed for cell proliferation and tumor growth, respectively. Mechanistic analyses of the transcriptome and kinome in DRD1-modulated cells informed further experiments, which characterized the effects on the epidermal growth factor receptor (EGFR) pathway and programmed cell death 1 ligand 1 (PD-L1) proteins. Through these experiments, we identified the DRD1 gene as a negative regulator of disease progression in NSCLC. We show that DRD1, as well as other DA pathway components, are expressed in normal human lung tissue, and that loss of DRD1 expression through promoter hypermethylation is a common feature in NSCLC patients and is associated with worse survival. At the cellular level, DRD1 affects proliferation by inhibiting the activation of EGFR and mitogen-activated protein kinase 1/2 (ERK1/2). Interestingly, we also found that DRD1 regulates the expression of PD-L1 in lung cancer cells. Taken together, these results suggest that DRD1 methylation may constitute a biomarker of poor prognosis in NSCLC patients while other components of this pathway could be targeted to improve response to EGFR- and PD-L1-targeted therapies.

Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression.

Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, here we analyze 1532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n = 148) or the alternatively spliced isoform (n = 9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal ( https://cseminer.stjude.org/ ), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

Diet-induced obesity affects influenza disease severity and transmission dynamics in ferrets.

Obesity, and the associated metabolic syndrome, is a risk factor for increased disease severity with a variety of infectious agents, including influenza virus. Yet, the mechanisms are only partially understood. As the number of people, particularly children, living with obesity continues to rise, it is critical to understand the role of host status on disease pathogenesis. In these studies, we use a diet-induced obese ferret model and tools to demonstrate that, like humans, obesity resulted in notable changes to the lung microenvironment, leading to increased clinical disease and viral spread to the lower respiratory tract. The decreased antiviral responses also resulted in obese animals shedding higher infectious virus for a longer period, making them more likely to transmit to contacts. These data suggest that the obese ferret model may be crucial to understanding obesity's impact on influenza disease severity and community transmission and a key tool for therapeutic and intervention development for this high-risk population.

Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression.

Immunotherapy with CAR T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons (CSE) present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify CSE targets, we analyzed 1,532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We found 2,933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n=148) or the alternatively spliced (AS) isoform (n=9) level. Expression of selected AS targets, including the EDB domain of FN1 (EDB), and gene targets, such as COL11A1, were validated in pediatric PDX tumors. We generated CAR T cells specific to EDB or COL11A1 and demonstrated that COL11A1-CAR T-cells have potent antitumor activity. The full target list, explorable via an interactive web portal (https://cseminer.stjude.org/), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

Conditional c-MYC activation in catecholaminergic cells drives distinct neuroendocrine tumors: neuroblastoma vs somatostatinoma.

The MYC proto-oncogenes (c-MYC, MYCN , MYCL ) are among the most deregulated oncogenic drivers in human malignancies including high-risk neuroblastoma, 50% of which are MYCN -amplified. Genetically engineered mouse models (GEMMs) based on the MYCN transgene have greatly expanded the understanding of neuroblastoma biology and are powerful tools for testing new therapies. However, a lack of c-MYC-driven GEMMs has hampered the ability to better understand mechanisms of neuroblastoma oncogenesis and therapy development given that c-MYC is also an important driver of many high-risk neuroblastomas. In this study, we report two transgenic murine neuroendocrine models driven by conditional c-MYC induction in tyrosine hydroxylase (Th) and dopamine ß-hydroxylase (Dbh)-expressing cells. c-MYC induction in Th-expressing cells leads to a preponderance of Pdx1 + somatostatinomas, a type of pancreatic neuroendocrine tumor (PNET), resembling human somatostatinoma with highly expressed gene signatures of δ cells and potassium channels. In contrast, c-MYC induction in Dbh-expressing cells leads to onset of neuroblastomas, showing a better transforming capacity than MYCN in a comparable C57BL/6 genetic background. The c-MYC murine neuroblastoma tumors recapitulate the pathologic and genetic features of human neuroblastoma, express GD2, and respond to anti-GD2 immunotherapy. This model also responds to DFMO, an FDA-approved inhibitor targeting ODC1, which is a known MYC transcriptional target. Thus, establishing c-MYC-overexpressing GEMMs resulted in different but related tumor types depending on the targeted cell and provide useful tools for testing immunotherapies and targeted therapies for these diseases.

Peptide-scFv antigen recognition domains effectively confer CAR T cell multiantigen specificity.

The emergence of immune escape is a significant roadblock to developing effective chimeric antigen receptor (CAR) T cell therapies against hematological malignancies, including acute myeloid leukemia (AML). Here, we demonstrate feasibility of targeting two antigens simultaneously by combining a GRP78-specific peptide antigen recognition domain with a CD123-specific scFv to generate a peptide-scFv bispecific antigen recognition domain (78.123). To achieve this, we test linkers with varying length and flexibility and perform immunophenotypic and functional characterization. We demonstrate that bispecific CAR T cells successfully recognize and kill tumor cells that express GRP78, CD123, or both antigens and have improved antitumor activity compared to their monospecific counterparts when both antigens are expressed. Protein structure prediction suggests that linker length and compactness influence the functionality of the generated bispecific CARs. Thus, we present a bispecific CAR design strategy to prevent immune escape in AML that can be extended to other peptide-scFv combinations.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

BACKGROUND: Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. RESULTS: Here, we generate single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We develop an unsupervised machine learning approach ("automatic consensus nonnegative matrix factorization" (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirm a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly, however, this weak-mesenchymal-like program is maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 h, suggesting an uncharacterized therapy-escape mechanism. CONCLUSIONS: Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Bone morphogenetic protein (BMP) signaling determines neuroblastoma cell fate and sensitivity to retinoic acid.

Retinoic acid (RA) is a standard-of-care neuroblastoma drug thought to be effective by inducing differentiation. Curiously, RA has little effect on primary human tumors during upfront treatment but can eliminate neuroblastoma cells from the bone marrow during post-chemo consolidation therapy-a discrepancy that has never been explained. To investigate this, we treated a large cohort of neuroblastoma cell lines with RA and observed that the most RA-sensitive cells predominantly undergo apoptosis or senescence, rather than differentiation. We conducted genome-wide CRISPR knockout screens under RA treatment, which identified BMP signaling as controlling the apoptosis/senescence vs differentiation cell fate decision and determining RA's overall potency. We then discovered that BMP signaling activity is markedly higher in neuroblastoma patient samples at bone marrow metastatic sites, providing a plausible explanation for RA's ability to clear neuroblastoma cells specifically from the bone marrow, seemingly mimicking interactions between BMP and RA during normal development.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

CD47 expression is critical for CAR T-cell survival in vivo.

BACKGROUND: CD47 is an attractive immunotherapeutic target because it is highly expressed on multiple solid tumors. However, CD47 is also expressed on T cells. Limited studies have evaluated CD47-chimeric antigen receptor (CAR) T cells, and the role of CD47 in CAR T-cell function remains largely unknown. METHODS: Here, we describe the development of CD47-CAR T cells derived from a high affinity signal regulatory protein α variant CV1, which binds CD47. CV1-CAR T cells were generated from human peripheral blood mononuclear cells and evaluated in vitro and in vivo. The role of CD47 in CAR T-cell function was examined by knocking out CD47 in T cells followed by downstream functional analyses. RESULTS: While CV1-CAR T cells are specific and exhibit potent activity in vitro they lacked antitumor activity in xenograft models. Mechanistic studies revealed CV1-CAR T cells downregulate CD47 to overcome fratricide, but CD47 loss resulted in their failure to expand and persist in vivo. This effect was not limited to CV1-CAR T cells, since CD47 knockout CAR T cells targeting another solid tumor antigen exhibited the same in vivo fate. Further, CD47 knockout T cells were sensitive to macrophage-mediated phagocytosis. CONCLUSIONS: These findings highlight that CD47 expression is critical for CAR T-cell survival in vivo and is a 'sine qua non' for successful adoptive T-cell therapy.

GRP78-CAR T cell effector function against solid and brain tumors is controlled by GRP78 expression on T cells.

Lack of targetable antigens is a key limitation for developing successful T cell-based immunotherapies. Members of the unfolded protein response (UPR) represent ideal immunotherapy targets because the UPR regulates the ability of cancer cells to resist cell death, sustain proliferation, and metastasize. Glucose-regulated protein 78 (GRP78) is a key UPR regulator that is overexpressed and translocated to the cell surface of a wide variety of cancers in response to elevated endoplasmic reticulum (ER) stress. We show that GRP78 is highly expressed on the cell surface of multiple solid and brain tumors, making cell surface GRP78 a promising chimeric antigen receptor (CAR) T cell target. We demonstrate that GRP78-CAR T cells can recognize and kill GRP78+ brain and solid tumors in vitro and in vivo. Additionally, our findings demonstrate that GRP78 is upregulated on CAR T cells upon T cell activation; however, this expression is tumor-cell-line specific and results in heterogeneous GRP78-CAR T cell therapeutic response.

A murine model of hnRNPH2-related neurodevelopmental disorder reveals a mechanism for genetic compensation by Hnrnph1.

Mutations in HNRNPH2 cause an X-linked neurodevelopmental disorder with features that include developmental delay, motor function deficits, and seizures. More than 90% of patients with hnRNPH2 have a missense mutation within or adjacent to the nuclear localization signal (NLS) of hnRNPH2. Here, we report that hnRNPH2 NLS mutations caused reduced interaction with the nuclear transport receptor Kapß2 and resulted in modest cytoplasmic accumulation of hnRNPH2. We generated 2 knockin mouse models with human-equivalent mutations in Hnrnph2 as well as Hnrnph2-KO mice. Knockin mice recapitulated clinical features of the human disorder, including reduced survival in male mice, impaired motor and cognitive functions, and increased susceptibility to audiogenic seizures. In contrast, 2 independent lines of Hnrnph2-KO mice showed no detectable phenotypes. Notably, KO mice had upregulated expression of Hnrnph1, a paralog of Hnrnph2, whereas knockin mice failed to upregulate Hnrnph1. Thus, genetic compensation by Hnrnph1 may counteract the loss of hnRNPH2. These findings suggest that HNRNPH2-related disorder may be driven by a toxic gain of function or a complex loss of HNRNPH2 function with impaired compensation by HNRNPH1. The knockin mice described here are an important resource for preclinical studies to assess the therapeutic benefit of gene replacement or knockdown of mutant hnRNPH2.

Modular chimeric cytokine receptors with leucine zippers enhance the antitumour activity of CAR T cells via JAK/STAT signalling.

The limited availability of cytokines in solid tumours hinders maintenance of the antitumour activity of chimeric antigen receptor (CAR) T cells. Cytokine receptor signalling pathways in CAR T cells can be activated by transgenic expression or injection of cytokines in the tumour, or by engineering the activation of cognate cytokine receptors. However, these strategies are constrained by toxicity arising from the activation of bystander cells, by the suboptimal biodistribution of the cytokines and by downregulation of the cognate receptor. Here we show that replacement of the extracellular domains of heterodimeric cytokine receptors in T cells with two leucine zipper motifs provides optimal Janus kinase/signal transducer and activator of transcription signalling. Such chimeric cytokine receptors, which can be generated for common γ-chain receptors, interleukin-10 and -12 receptors, enabled T cells to survive cytokine starvation without induction of autonomous cell growth, and augmented the effector function of CAR T cells in vitro in the setting of chronic antigen exposure and in human tumour xenografts in mice. As a modular design, leucine zippers can be used to generate constitutively active cytokine receptors in effector immune cells.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma.

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

Haploinsufficiency of the lysosomal sialidase NEU1 results in a model of pleomorphic rhabdomyosarcoma in mice.

Rhabdomyosarcoma, the most common pediatric sarcoma, has no effective treatment for the pleomorphic subtype. Still, what triggers transformation into this aggressive phenotype remains poorly understood. Here we used Ptch1+/-/ETV7TG/+/- mice with enhanced incidence of rhabdomyosarcoma to generate a model of pleomorphic rhabdomyosarcoma driven by haploinsufficiency of the lysosomal sialidase neuraminidase 1. These tumors share mostly features of embryonal and some of alveolar rhabdomyosarcoma. Mechanistically, we show that the transforming pathway is increased lysosomal exocytosis downstream of reduced neuraminidase 1, exemplified by the redistribution of the lysosomal associated membrane protein 1 at the plasma membrane of tumor and stromal cells. Here we exploit this unique feature for single cell analysis and define heterogeneous populations of exocytic, only partially differentiated cells that force tumors to pleomorphism and promote a fibrotic microenvironment. These data together with the identification of an adipogenic signature shared by human rhabdomyosarcoma, and likely fueling the tumor's metabolism, make this model of pleomorphic rhabdomyosarcoma ideal for diagnostic and therapeutic studies.