Beta cell-specific CD8+ T cells maintain stem cell memory-associated epigenetic programs during type 1 diabetes.

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.

A kinase-independent function of cyclin-dependent kinase 6 promotes outer radial glia expansion and neocortical folding.

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.

Duocarmycin SA Reduces Proliferation and Increases Apoptosis in Acute Myeloid Leukemia Cells In Vitro.

Acute myeloid leukemia (AML) is a hematological malignancy that is characterized by an expansion of immature myeloid precursors. Despite therapeutic advances, the prognosis of AML patients remains poor and there is a need for the evaluation of promising therapeutic candidates to treat the disease. The objective of this study was to evaluate the efficacy of duocarmycin Stable A (DSA) in AML cells in vitro. We hypothesized that DSA would induce DNA damage in the form of DNA double-strand breaks (DSBs) and exert cytotoxic effects on AML cells within the picomolar range. Human AML cell lines Molm-14 and HL-60 were used to perform 3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide (MTT), DNA DSBs, cell cycle, 5-ethynyl-2-deoxyuridine (EdU), colony formation unit (CFU), Annexin V, RNA sequencing and other assays described in this study. Our results showed that DSA induced DNA DSBs, induced cell cycle arrest at the G2M phase, reduced proliferation and increased apoptosis in AML cells. Additionally, RNA sequencing results showed that DSA regulates genes that are associated with cellular processes such as DNA repair, G2M checkpoint and apoptosis. These results suggest that DSA is efficacious in AML cells and is therefore a promising potential therapeutic candidate that can be further evaluated for the treatment of AML.

The molecular characteristics of low-grade and high-grade areas in desmoplastic infantile astrocytoma/ganglioglioma.

AIMS: Desmoplastic infantile astrocytomas and gangliogliomas (DIA/DIGs) are rare brain tumours of infancy. A distinctive feature of their histopathology is a combination of low-grade and high-grade features. Most DIA/DIGs can be surgically resected and have a good prognosis. However, high-grade features often dominate recurrent tumours, some of which have a poor outcome. In this study, we test the hypothesis that low-grade and high-grade areas in DIA/DIGs have distinct molecular characteristics. METHODS: Tissue samples from microdissected low-grade and high-grade areas in 12 DIA/DIGs were analysed by DNA methylation profiling, whole exome sequencing, RNA sequencing and immunohistochemistry to search for potential differences at multiple molecular levels. RESULTS: Copy number variants among tumours and between the two morphologically distinct areas were infrequent. No recurrent genetic alterations were identified across the tumour series, and high-grade areas did not have additional genetic alterations to explain their distinct morphology or biological behaviour. However, high-grade areas showed relative hypomethylation in genes downstream of the transcription factors SOX9 and LEF1 and evidence of a core SOX9 transcription network alongside activation of the BMP, WNT and MAPK signalling pathways. CONCLUSIONS: This study contributes to our knowledge of molecular genetic alterations in DIA/DIGs, uncovers molecular differences between the two distinct cell populations in these tumours and suggests potential therapeutic targets among the more proliferative cell population in DIA/DIGs.

ChIPseqSpikeInFree: a ChIP-seq normalization approach to reveal global changes in histone modifications without spike-in.

MOTIVATION: The traditional reads per million normalization method is inappropriate for the evaluation of ChIP-seq data when treatments or mutations have global effects. Changes in global levels of histone modifications can be detected with exogenous reference spike-in controls. However, most ChIP-seq studies overlook the normalization that must be corrected with spike-in. A method that retrospectively renormalizes datasets without spike-in is lacking. RESULTS: ChIPseqSpikeInFree is a novel ChIP-seq normalization method to effectively determine scaling factors for samples across various conditions and treatments, which does not rely on exogenous spike-in chromatin or peak detection to reveal global changes in histone modification occupancy. Application of ChIPseqSpikeInFree on five datasets demonstrates that this in silico approach reveals a similar magnitude of global changes as the spike-in method does. AVAILABILITY AND IMPLEMENTATION: St. Jude Cloud (https://pecan.stjude.cloud/permalink/spikefree) and St. Jude Github ( https://github.com/stjude/ChIPseqSpikeInFree). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Correction to: H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

The original article can be found online.

H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

Histone H3 K27M mutation is the defining molecular feature of the devastating pediatric brain tumor, diffuse intrinsic pontine glioma (DIPG). The prevalence of histone H3 K27M mutations indicates a critical role in DIPGs, but the contribution of the mutation to disease pathogenesis remains unclear. We show that knockdown of this mutation in DIPG xenografts restores K27M-dependent loss of H3K27me3 and delays tumor growth. Comparisons of matched DIPG xenografts with and without K27M knockdown allowed identification of mutation-specific effects on the transcriptome and epigenome. The resulting transcriptional changes recapitulate expression signatures from K27M primary DIPG tumors and are strongly enriched for genes associated with nervous system development. Integrated analysis of ChIP-seq and expression data showed that genes upregulated by the mutation are overrepresented in apparently bivalent promoters. Many of these targets are associated with more immature differentiation states. Expression profiles indicate K27M knockdown decreases proliferation and increases differentiation within lineages represented in DIPG. These data suggest that K27M-mediated loss of H3K27me3 directly regulates a subset of genes by releasing poised promoters, and contributes to tumor phenotype and growth by limiting differentiation. The delayed tumor growth associated with knockdown of H3 K27M provides evidence that this highly recurrent mutation is a relevant therapeutic target.

The Histone Acetyltransferase GCN5 Expression Is Elevated and Regulated by c-Myc and E2F1 Transcription Factors in Human Colon Cancer.

The histone acetyltransferase GCN5 has been suggested to be involved in promoting cancer cell growth. But its role in human colon cancer development remains unknown. Herein we discovered that GCN5 expression is significantly upregulated in human colon adenocarcinoma tissues. We further demonstrate that GCN5 is upregulated in human colon cancer at the mRNA level. Surprisingly, two transcription factors, the oncogenic c-Myc and the proapoptotic E2F1, are responsible for GCN5 mRNA transcription. Knockdown of c-Myc inhibited colon cancer cell proliferation largely through downregulating GCN5 transcription, which can be fully rescued by the ectopic GCN5 expression. In contrast, E2F1 expression induced human colon cancer cell death, and suppression of GCN5 expression in cells with E2F1 overexpression further facilitated cell apoptosis, suggesting that GCN5 expression is induced by E2F1 as a possible negative feedback in suppressing E2F1-mediated cell apoptosis. In addition, suppression of GCN5 with its specific inhibitor CPTH2 inhibited human colon cancer cell growth. Our studies reveal that GCN5 plays a positive role in human colon cancer development, and its suppression holds a great therapeutic potential in antitumor therapy.

Synthesis of the Stable Ordered Conjugated Polymer Poly(dibromodiacetylene) from an Explosive Monomer.

Dibromobutadiyne is an extremely unstable compound that explodes at room temperature, even under inert atmosphere. This instability has limited the studies of dibromobutadiyne almost entirely to spectroscopic characterization. Here we report an approach to control the reactivity of dibromobutadiyne, via topochemical reaction in cocrystals, leading to the ordered polymer poly(dibromodiacetylene), PBDA. At low temperatures (-15 to -18 °C), dibromobutadiyne can form cocrystals with oxalamide host molecules containing either pyridyl or nitrile side groups, in which halogen bonds align the dibromobutadiyne monomers for topochemical polymerization. The cocrystals with the bis(nitrile) oxalamide host undergo complete ordered polymerization to PBDA, demonstrated by solid-state MAS-NMR, Raman, and optical absorption spectroscopy. Once formed, the polymer can be separated from the host; unlike the monomer, PBDA is stable at room temperature.

The context-dependent epigenetic and organogenesis programs determine 3D vs. 2D cellular fitness of MYC-driven cancer.

3D cellular-specific epigenetic and transcriptomic reprogramming is critical to organogenesis and tumorigenesis. Here we dissect the distinct cell fitness in 2D (normoxia vs. chronic hypoxia) vs 3D (normoxia) culture conditions. We identify over 600 shared essential genes and additional context-specific fitness genes and pathways. Knockout of the VHL-HIF1 pathway results in incompatible fitness defects under normoxia vs. 1% oxygen or 3D culture conditions. Moreover, deletion of each of the mitochondrial respiratory electron transport chain complex has distinct fitness outcomes. Notably, multicellular organogenesis signaling pathways including TGFß-SMAD specifically constrict the uncontrolled cell proliferation in 3D while inactivation of epigenetic modifiers (Bcor, Kmt2d, Mettl3 and Mettl14) has opposite outcomes in 2D vs. 3D. We further identify a 3D-dependent synthetic lethality with partial loss of Prmt5 due to a reduction of Mtap expression resulting from 3D-specific epigenetic reprogramming. Our study highlights unique epigenetic, metabolic and organogenesis signaling dependencies under different cellular settings.

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.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing associated with therapeutic response to splicing inhibitor.

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.

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma.

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

RBM39 degrader invigorates natural killer cells to eradicate neuroblastoma despite cancer cell plasticity.

The cellular plasticity of neuroblastoma is defined by a mixture of two major cell states, adrenergic (ADRN) and mesenchymal (MES), which may contribute to therapy resistance. However, how neuroblastoma cells switch cellular states during therapy remains largely unknown and how to eradicate neuroblastoma regardless of their cell states is a clinical challenge. To better understand the lineage switch of neuroblastoma in chemoresistance, we comprehensively defined the transcriptomic and epigenetic map of ADRN and MES types of neuroblastomas using human and murine models treated with indisulam, a selective RBM39 degrader. We showed that cancer cells not only undergo a bidirectional switch between ADRN and MES states, but also acquire additional cellular states, reminiscent of the developmental pliancy of neural crest cells. The lineage alterations are coupled with epigenetic reprogramming and dependency switch of lineage-specific transcription factors, epigenetic modifiers and targetable kinases. Through targeting RNA splicing, indisulam induces an inflammatory tumor microenvironment and enhances anticancer activity of natural killer cells. The combination of indisulam with anti-GD2 immunotherapy results in a durable, complete response in high-risk transgenic neuroblastoma models, providing an innovative, rational therapeutic approach to eradicate tumor cells regardless of their potential to switch cell states.

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.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication.

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.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing is associated with therapeutic response to splicing inhibitor.

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 neuroblastoma. JMJD6 cooperates with MYC in cellular transformation 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.

PAX3-FOXO1 dictates myogenic reprogramming and rhabdomyosarcoma identity in endothelial progenitors.

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.

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.

Targeting KDM4 for treating PAX3-FOXO1-driven alveolar rhabdomyosarcoma.

Chimeric transcription factors drive lineage-specific oncogenesis but are notoriously difficult to target. Alveolar rhabdomyosarcoma (RMS) is an aggressive childhood soft tissue sarcoma transformed by the pathognomonic Paired Box 3-Forkhead Box O1 (PAX3-FOXO1) fusion protein, which governs a core regulatory circuitry transcription factor network. Here, we show that the histone lysine demethylase 4B (KDM4B) is a therapeutic vulnerability for PAX3-FOXO1+ RMS. Genetic and pharmacologic inhibition of KDM4B substantially delayed tumor growth. Suppression of KDM4 proteins inhibited the expression of core oncogenic transcription factors and caused epigenetic alterations of PAX3-FOXO1-governed superenhancers. Combining KDM4 inhibition with cytotoxic chemotherapy led to tumor regression in preclinical PAX3-FOXO1+ RMS subcutaneous xenograft models. In summary, we identified a targetable mechanism required for maintenance of the PAX3-FOXO1-related transcription factor network, which may translate to a therapeutic approach for fusion-positive RMS.