Rhabdomyosarcoma (RMS) is a pediatric cancer with features of skeletal muscle; patients with unresectable or metastatic RMS fare poorly due to high rates of disease recurrence. Here, we use single-cell and single-nucleus RNA sequencing to show that RMS tumors recapitulate the spectrum of embryonal myogenesis. Using matched patient samples from a clinical trial and orthotopic patient-derived xenografts (O-PDXs), we show that chemotherapy eliminates the most proliferative component with features of myoblasts within embryonal RMS; after treatment, the immature population with features of paraxial mesoderm expands to reconstitute the developmental hierarchy of the original tumor. We discovered that this paraxial mesoderm population is dependent on EGFR signaling and is sensitive to EGFR inhibitors. Taken together, these data serve as a proof of concept that targeting each developmental state in embryonal RMS is an effective strategy for improving outcomes by preventing disease recurrence.

Retinoblastoma is a childhood cancer of the developing retina that initiates with biallelic inactivation of the RB1 gene. Children with germline mutations in RB1 have a high likelihood of developing retinoblastoma and other malignancies later in life. Genetically engineered mouse models of retinoblastoma share some similarities with human retinoblastoma but there are differences in their cellular differentiation. To develop a laboratory model of human retinoblastoma formation, we make induced pluripotent stem cells (iPSCs) from 15 participants with germline RB1 mutations. Each of the stem cell lines is validated, characterized and then differentiated into retina using a 3-dimensional organoid culture system. After 45 days in culture, the retinal organoids are dissociated and injected into the vitreous of eyes of immunocompromised mice to support retinoblastoma tumor growth. Retinoblastomas formed from retinal organoids made from patient-derived iPSCs have molecular, cellular and genomic features indistinguishable from human retinoblastomas. This model of human cancer based on patient-derived iPSCs with germline cancer predisposing mutations provides valuable insights into the cellular origins of this debilitating childhood disease as well as the mechanism of tumorigenesis following RB1 gene inactivation.

Malignant peripheral nerve sheath tumor (MPNST), a highly aggressive Schwann cell (SC)-derived soft tissue sarcoma, arises from benign neurofibroma (NF); however, the identity, heterogeneity and origins of tumor populations remain elusive. Nestin+ cells have been implicated as tumor stem cells in MPNST; unexpectedly, single-cell profiling of human NF and MPNST and their animal models reveal a broad range of nestin-expressing SC lineage cells and dynamic acquisition of discrete cancer states during malignant transformation. We uncover a nestin-negative mesenchymal neural crest-like subpopulation as a previously unknown malignant stem-like state common to murine and human MPNSTs, which correlates with clinical severity. Integrative multiomics profiling further identifies unique regulatory networks and druggable targets against the malignant subpopulations in MPNST. Targeting key epithelial-mesenchymal transition and stemness regulators including ZEB1 and ALDH1A1 impedes MPNST growth. Together, our studies reveal the underlying principles of tumor cell-state evolution and their regulatory circuitries during NF-to-MPNST transformation, highlighting a hitherto unrecognized mesenchymal stem-like subpopulation in MPNST disease progression.

ProteinPaint BAM track (ppBAM) is designed to assist variant review for cancer research and clinical genomics. With performant server-side computing and rendering, ppBAM supports on-the-fly variant genotyping of thousands of reads using Smith-Waterman alignment. To better visualize support for complex variants, reads are realigned against the mutated reference sequence using ClustalO. ppBAM also supports the BAM slicing API of the NCI Genomic Data Commons (GDC) portal, letting researchers conveniently examine genomic details of vast amounts of cancer sequencing data and reinterpret variant calls.

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.