The mitotic deacetylase complex (MiDAC) is a recently identified histone deacetylase (HDAC) complex. While other HDAC complexes have been implicated in neurogenesis, the physiological role of MiDAC remains unknown. Here, we show that MiDAC constitutes an important regulator of neural differentiation. We demonstrate that MiDAC functions as a modulator of a neurodevelopmental gene expression program and binds to important regulators of neurite outgrowth. MiDAC upregulates gene expression of pro-neural genes such as those encoding the secreted ligands SLIT3 and NETRIN1 (NTN1) by a mechanism suggestive of H4K20ac removal on promoters and enhancers. Conversely, MiDAC inhibits gene expression by reducing H3K27ac on promoter-proximal and -distal elements of negative regulators of neurogenesis. Furthermore, loss of MiDAC results in neurite outgrowth defects that can be rescued by supplementation with SLIT3 and/or NTN1. These findings indicate a crucial role for MiDAC in regulating the ligands of the SLIT3 and NTN1 signaling axes to ensure the proper integrity of neurite development.

The neocortex, the center for higher brain function, first emerged in mammals and has become massively expanded and folded in humans, constituting almost half the volume of the human brain. Primary microcephaly, a developmental disorder in which the brain is smaller than normal at birth, results mainly from there being fewer neurons in the neocortex because of defects in neural progenitor cells (NPCs). Outer radial glia (oRGs), NPCs that are abundant in gyrencephalic species but rare in lissencephalic species, are thought to play key roles in the expansion and folding of the neocortex. However, how oRGs expand, whether they are necessary for neocortical folding, and whether defects in oRGs cause microcephaly remain important questions in the study of brain development, evolution, and disease. Here, we show that oRG expansion in mice, ferrets, and human cerebral organoids requires cyclin-dependent kinase 6 (CDK6), the mutation of which causes primary microcephaly via an unknown mechanism. In a mouse model in which increased Hedgehog signaling expands oRGs and intermediate progenitor cells and induces neocortical folding, CDK6 loss selectively decreased oRGs and abolished neocortical folding. Remarkably, this function of CDK6 in oRG expansion did not require its kinase activity, was not shared by the highly similar CDK4 and CDK2, and was disrupted by the mutation causing microcephaly. Therefore, our results indicate that CDK6 is conserved to promote oRG expansion, that oRGs are necessary for neocortical folding, and that defects in oRG expansion may cause primary microcephaly.

PTEN promoter hypermethylation is nearly universal and PTEN copy number loss occurs in ~25% of fusion-negative rhabdomyosarcoma (FN-RMS). Here we show Pten deletion in a mouse model of FN-RMS results in less differentiated tumors more closely resembling human embryonal RMS. PTEN loss activated the PI3K pathway but did not increase mTOR activity. In wild-type tumors, PTEN was expressed in the nucleus suggesting loss of nuclear PTEN functions could account for these phenotypes. Pten deleted tumors had increased expression of transcription factors important in neural and skeletal muscle development including Dbx1 and Pax7. Pax7 deletion completely rescued the effects of Pten loss. Strikingly, these Pten;Pax7 deleted tumors were no longer FN-RMS but displayed smooth muscle differentiation similar to leiomyosarcoma. These data highlight how Pten loss in FN-RMS is connected to a PAX7 lineage-specific transcriptional output that creates a dependency or synthetic essentiality on the transcription factor PAX7 to maintain tumor identity.

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans could cause coronavirus disease 2019 (COVID-19). Since its first discovery in Dec 2019, SARS-CoV-2 has become a global pandemic and caused 3.3 million direct/indirect deaths (2021 May). Amongst the scientific community's response to COVID-19, data sharing has emerged as an essential aspect of the combat against SARS-CoV-2. Despite the ever-growing studies about SARS-CoV-2 and COVID-19, to date, only a few databases were curated to enable access to gene expression data. Furthermore, these databases curated only a small set of data and do not provide easy access for investigators without computational skills to perform analyses. To fill this gap and advance open-access to the growing gene expression data on this deadly virus, we collected about 1,500 human bulk RNA-seq datasets from publicly available resources, developed a database and visualization tool, named CovidExpress (https://stjudecab.github.io/covidexpress). This open access database will allow research investigators to examine the gene expression in various tissues, cell lines, and their response to SARS-CoV-2 under different experimental conditions, accelerating the understanding of the etiology of this disease to inform the drug and vaccine development. Our integrative analysis of this big dataset highlights a set of commonly regulated genes in SARS-CoV-2 infected lung and Rhinovirus infected nasal tissues, including OASL that were under-studied in COVID-19 related reports. Our results also suggested a potential FURIN positive feedback loop that might explain the evolutional advantage of SARS-CoV-2.

Diffuse intrinsic pontine gliomas (DIPGs) are incurable childhood brainstem tumors with frequent histone H3 K27M mutations and recurrent alterations in PDGFRA and TP53. We generated genetically engineered inducible mice and showed that H3.3 K27M enhanced neural stem cell self-renewal while preserving regional identity. Neonatal induction of H3.3 K27M cooperated with activating platelet-derived growth factor receptor α (PDGFRα) mutant and Trp53 loss to accelerate development of diffuse brainstem gliomas that recapitulated human DIPG gene expression signatures and showed global changes in H3K27 posttranslational modifications, but relatively restricted gene expression changes. Genes upregulated in H3.3 K27M tumors were enriched for those associated with neural development where H3K27me3 loss released the poised state of apparently bivalent promoters, whereas downregulated genes were enriched for those encoding homeodomain transcription factors.

The pool of beta cell-specific CD8+ T cells in type 1 diabetes (T1D) sustains an autoreactive potential despite having access to a constant source of antigen. To investigate the long-lived nature of these cells, we established a DNA methylation-based T cell 'multipotency index' and found that beta cell-specific CD8+ T cells retained a stem-like epigenetic multipotency score. Single-cell assay for transposase-accessible chromatin using sequencing confirmed the coexistence of naive and effector-associated epigenetic programs in individual beta cell-specific CD8+ T cells. Assessment of beta cell-specific CD8+ T cell anatomical distribution and the establishment of stem-associated epigenetic programs revealed that self-reactive CD8+ T cells isolated from murine lymphoid tissue retained developmentally plastic phenotypic and epigenetic profiles relative to the same cells isolated from the pancreas. Collectively, these data provide new insight into the longevity of beta cell-specific CD8+ T cell responses and document the use of this methylation-based multipotency index for investigating human and mouse CD8+ T cell differentiation.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

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.

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that jumonji domain containing 6, arginine demethylase, and lysine hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven human neuroblastoma. JMJD6 cooperates with MYC in cellular transformation of murine neural crest cells by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a 'molecular glue' that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

Malignant peripheral nerve sheath tumors (MPNST) are aggressive sarcomas. Somatic inactivation of NF1 and cooperating tumor suppressors, including CDKN2A/B, PRC2, and p53, is found in most MPNST. Inactivation of LATS1/2 of the Hippo pathway was recently shown to cause tumors resembling MPNST histologically, although Hippo pathway mutations are rarely found in MPNST. Because existing genetically engineered mouse (GEM) models of MPNST do not recapitulate some of the key genetic features of human MPNST, we aimed to establish a GEM-MPNST model that recapitulated the human disease genetically, histologically, and molecularly.

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.

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.

The H3K27me2/me3 histone demethylase KDM6B is essential to neuroblastoma cell survival. However, the mechanism of KDM6B action remains poorly defined. We demonstrate that inhibition of KDM6B activity 1) reduces the chromatin accessibility of E2F target genes and MYCN, 2) selectively leads to an increase of H3K27me3 but a decrease of the enhancer mark H3K4me1 at the CTCF and BORIS binding sites, which may, consequently, disrupt the long-range chromatin interaction of MYCN and E2F target genes, and 3) phenocopies the transcriptome induced by the specific CDK4/6 inhibitor palbociclib. Overexpression of CDK4/6 or Rb1 knockout confers neuroblastoma cell resistance to both palbociclib and the KDM6 inhibitor GSK-J4. These data indicate that KDM6B promotes an oncogenic CDK4/6-pRB-E2F pathway in neuroblastoma cells via H3K27me3-dependent enhancer-promoter interactions, providing a rationale to target KDM6B for high-risk neuroblastoma.

Aberrant alternative pre-mRNA splicing plays a critical role in MYC-driven cancers and therefore may represent a therapeutic vulnerability. Here, we show that neuroblastoma, a MYC-driven cancer characterized by splicing dysregulation and spliceosomal dependency, requires the splicing factor RBM39 for survival. Indisulam, a “molecular glue” that selectively recruits RBM39 to the CRL4-DCAF15 E3 ubiquitin ligase for proteasomal degradation, is highly efficacious against neuroblastoma, leading to significant responses in multiple high-risk disease models, without overt toxicity. Genetic depletion or indisulam-mediated degradation of RBM39 induces significant genome-wide splicing anomalies and cell death. Mechanistically, the dependency on RBM39 and high-level expression of DCAF15 determine the exquisite sensitivity of neuroblastoma to indisulam. Our data indicate that targeting the dysregulated spliceosome by precisely inhibiting RBM39, a vulnerability in neuroblastoma, is a valid therapeutic strategy.

Histone lysine demethylases (KDMs) play critical roles in oncogenesis and therefore may be effective targets for anticancer therapy. Using a time-resolved fluorescence resonance energy transfer demethylation screen assay, in combination with multiple orthogonal validation approaches, we identified geldanamycin and its analog 17-DMAG as KDM inhibitors. In addition, we found that these Hsp90 inhibitors increase degradation of the alveolar rhabdomyosarcoma (aRMS) driver oncoprotein PAX3-FOXO1 and induce the repressive epigenetic mark H3K9me3 and H3K36me3 at genomic loci of PAX3-FOXO1 targets. We found that as monotherapy 17-DMAG significantly inhibits expression of PAX3-FOXO1 target genes and multiple oncogenic pathways, induces a muscle differentiation signature, delays tumor growth and extends survival in aRMS xenograft mouse models. The combination of 17-DMAG with conventional chemotherapy significantly enhances therapeutic efficacy, indicating that targeting KDM in combination with chemotherapy may serve as a therapeutic approach to PAX3-FOXO1-positive aRMS.

Fusion-positive rhabdomyosarcoma (FP-RMS) driven by the expression of the PAX3-FOXO1 (P3F) fusion oncoprotein is an aggressive subtype of pediatric rhabdomyosarcoma. FP-RMS histologically resembles developing muscle yet occurs throughout the body in areas devoid of skeletal muscle highlighting that FP-RMS is not derived from an exclusively myogenic cell of origin. Here we demonstrate that P3F reprograms mouse and human endothelial progenitors to FP-RMS. We show that P3F expression in aP2-Cre expressing cells reprograms endothelial progenitors to functional myogenic stem cells capable of regenerating injured muscle fibers. Further, we describe a FP-RMS mouse model driven by P3F expression and Cdkn2a loss in endothelial cells. Additionally, we show that P3F expression in TP53-null human iPSCs blocks endothelial-directed differentiation and guides cells to become myogenic cells that form FP-RMS tumors in immunocompromised mice. Together these findings demonstrate that FP-RMS can originate from aberrant development of non-myogenic cells driven by P3F.