Intrinsic transcriptional heterogeneity in neuroblastoma guides mechanistic and therapeutic insights

Mutually exclusive teams-like patterns of gene regulation characterize phenotypic heterogeneity along the noradrenergic-mesenchymal axis in neuroblastoma

Neuroblastoma is the most frequent extracranial pediatric tumor and leads to 15% of all cancer-related deaths in children. Tumor relapse and therapy resistance in neuroblastoma are driven by phenotypic plasticity and heterogeneity between noradrenergic (NOR) and mesenchymal (MES) cell states. Despite the importance of this phenotypic plasticity, the dynamics and molecular patterns associated with these bidirectional cell-state transitions remain relatively poorly understood. Here, we analyze multiple RNA-seq datasets at both bulk and single-cell resolution, to understand the association between NOR- and MES-specific factors. We observed that NOR-specific and MES-specific expression patterns are largely mutually exclusive, exhibiting a "teams-like" behavior among the genes involved, reminiscent of our earlier observations in lung cancer and melanoma. This antagonism between NOR and MES phenotypes was also associated with metabolic reprogramming and with immunotherapy targets PD-L1 and GD2 as well as with experimental perturbations driving the NOR-MES and/or MES-NOR transition. Further, these "teams-like" patterns were seen only among the NOR- and MES-specific genes, but not in housekeeping genes, possibly highlighting a hallmark of network topology enabling cancer cell plasticity.

Novel miRNA-inducing drugs enable differentiation of retinoic acid-resistant neuroblastoma cells

Tumor cell heterogeneity in neuroblastoma, a pediatric cancer arising from neural crest-derived progenitor cells, poses a significant clinical challenge. In particular, unlike adrenergic (ADRN) neuroblastoma cells, mesenchymal (MES) cells are resistant to chemotherapy and retinoid therapy and thereby significantly contribute to relapses and treatment failures. Previous research suggested that overexpression or activation of miR-124, a neurogenic microRNA with tumor suppressor activity, can induce the differentiation of retinoic acid-resistant neuroblastoma cells. Leveraging our established screen for miRNA modulatory small molecules, we validated PP121, a dual inhibitor of tyrosine and phosphoinositide kinases, as a robust inducer of miR-124. A combination of PP121 and miR-132-inducing bufalin synergistically arrests proliferation, induces differentiation, and prolongs the survival of differentiated MES SK-N-AS cells for 8 weeks. RNA- seq and deconvolution analyses revealed a collapse of the ADRN core regulatory circuitry (CRC) and the emergence of novel CRCs associated with chromaffin cells and Schwann cell precursors. Using a similar protocol, we differentiated and maintained other MES neuroblastoma, as well as glioblastoma cells, over 16 weeks. In conclusion, our novel protocol suggests a promising treatment for therapy-resistant cancers of the nervous system. Moreover, these long-lived, differentiated cells provide valuable models for studying mechanisms underlying differentiation, maturation, and senescence.

Group 3 medulloblastoma transcriptional networks collapse under domain specific EP300/CBP inhibition

Chemical discovery efforts commonly target individual protein domains. Many proteins, including the EP300/CBP histone acetyltransferases (HATs), contain several targetable domains. EP300/CBP are critical gene-regulatory targets in cancer, with existing high potency inhibitors of either the catalytic HAT domain or protein-binding bromodomain (BRD). A domain-specific inhibitory approach to multidomain-containing proteins may identify exceptional-responding tumor types, thereby expanding a therapeutic index. Here, we discover that targeting EP300/CBP using the domain-specific inhibitors, A485 (HAT) or CCS1477 (BRD) have different effects in select tumor types. Group 3 medulloblastoma (G3MB) cells are especially sensitive to BRD, compared with HAT inhibition. Structurally, these effects are mediated by the difluorophenyl group in the catalytic core of CCS1477. Mechanistically, bromodomain inhibition causes rapid disruption of genetic dependency networks that are required for G3MB growth. These studies provide a domain-specific structural foundation for drug discovery efforts targeting EP300/CBP and identify a selective role for the EP300/CBP bromodomain in maintaining genetic dependency networks in G3MB.

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.

PRMT5 orchestrates EGFR and AKT networks to activate NFκB and promote EMT

Neuroblastoma remains a formidable challenge in pediatric oncology, representing 15% of cancer-related mortalities in children. Despite advancements in combinatorial and targeted treatments improving survival rates, nearly 50% of patients with high-risk neuroblastoma will ultimately succumb to their disease. Dysregulation of the epithelial-mesenchymal transition (EMT) is a key mechanism of tumor cell dissemination, resulting in metastasis and poor outcomes in many cancers. Our prior work identified PRMT5 as a key regulator of EMT via methylation of AKT at arginine 15, enhancing the expression of EMT-driving transcription factors and facilitating metastasis. Here, we identify that PRMT5 directly regulates the transcription of the epidermal growth factor receptor (EGFR). PRMT5, through independent modulation of the EGFR and AKT pathways, orchestrates the activation of NFκB, resulting in the upregulation of the pro-EMT transcription factors ZEB1, SNAIL, and TWIST1. Notably, EGFR and AKT form a compensatory feedback loop, reinforcing the expression of these EMT transcription factors. Small molecule inhibition of PRMT5 methyltransferase activity disrupts EGFR/AKT signaling, suppresses EMT transcription factor expression and ablates tumor growth in vivo . Our findings underscore the pivotal role of PRMT5 in the control of the EMT program in high-risk neuroblastoma.

Functionalized GD2 Electrochemical Immunosensor to Diagnose Minimum Residual Disease of Bone Marrow in Neuroblastoma Effectively

Neuroblastoma (NB) is known as the "king of childhood tumors" due to its highly metastatic, recurrence-prone, and difficult-to-treat characteristics. International Neuroblastoma Risk Grading Group (INRG) has recommended GD2, a disialoganglioside expressed on neuroectodermal tumor cells, as the target for detecting minimal residual disease in bone marrow metastases of high-risk neuroblastoma in children. Therefore, accurately identifying GD2-positive cells is crucial for diagnosing children with high-risk NB. Here, we designed a graphene/AuNP/GD2 Ab-functionalized electrochemical biosensor for GD2 detection. A three-electrode system was processed using a screen-printed technique with a working electrode of indium tin oxide, a counter electrode of carbon, and a reference electrode of silver/silver chloride. Graphene/AuNPs were modified on the indium tin oxide electrode using chronoamperometric scans, and then, the GD2 antibody was modified on the biosensor by electrostatic adsorption to achieve sensitive and specific detection of GD2-positive cells in bone marrow fluid. The results showed that a graphene/AuNP/GD2 Ab-functionalized electrochemical biosensor achieved GD2-positive cell detection in the range of 102 cells/mL~105 cells/mL by differential pulse voltammetry. Bone marrow fluid samples from 12 children with high-risk NB were retained for testing on our biosensor and showed 100% compliance with the clinical application of the gold-standard immunocytochemical staining technique for detecting GD2-positive cells qualitatively. The GD2-based electrochemical assay can accurately detect children with high-risk NB, providing a rapidly quantitative basis for clinical diagnosis and treatment.

T-cell inflammation is prognostic of survival in patients with high-risk neuroblastoma enriched for an adrenergic signature

Purpose

T-cell inflammation (TCI) has been shown to be a prognostic marker in neuroblastoma, a tumor comprised of cells that can exist in two epigenetic states, adrenergic (ADRN) and mesenchymal (MES). We hypothesized that elucidating unique and overlapping aspects of these biologic features could serve as novel biomarkers.

Patients and Methods

We detected lineage-specific, single-stranded super-enhancers defining ADRN and MES specific genes. Publicly available neuroblastoma RNA-seq data from GSE49711 (Cohort 1) and TARGET (Cohort 2) were assigned MES, ADRN, and TCI scores. Tumors were characterized as MES (top 33%) or ADRN (bottom 33%), and TCI (top 67% TCI score) or non-inflamed (bottom 33% TCI score). Overall survival (OS) was assessed using the Kaplan-Meier method, and differences were assessed by the log-rank test.

Results

We identified 159 MES genes and 373 ADRN genes. TCI scores were correlated with MES scores (R=0.56, p<0.001 and R=0.38, p<0.001) and anticorrelated with MYCN -amplification (R=-0.29, p<0.001 and -0.18, p=0.03) in both cohorts. Among Cohort 1 patients with high-risk, ADRN tumors (n=59), those with TCI tumors (n=22) had superior OS to those with non-inflammed tumors (n=37) (p=0.01), though this comparison did not reach significance in Cohort 2. TCI status was not associated with survival in patients with high-risk MES tumors in either cohort.

Conclusions

High inflammation scores were correlated with improved survival in some high-risk patients with, ADRN but not MES neuroblastoma. These findings have implications for approaches to treating high-risk neuroblastoma.

Genetic predisposition to neuroblastoma results from a regulatory polymorphism that promotes the adrenergic cell state

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

Modulation of DNA-protein interactions by proximal genetic elements as uncovered by interpretable deep learning

Transcription factors (TF) recognize specific motifs in the genome that are typically 6-12 bp long to regulate various aspects of the cellular machinery. Presence of binding motifs and favorable genome accessibility are key drivers for a consistent TF-DNA interaction. Although these pre-requisites may occur thousands of times in the genome, there seems to be a high degree of selectivity for the sites that are actually bound. Here, we present a deep-learning framework that identifies and characterizes the upstream and downstream genetic elements to the binding motif, for their role in enforcing the mentioned selectivity. The proposed framework is based on an interpretable recurrent neural network architecture that enables for the relative analysis of sequence context features. We apply the framework to model twenty-six transcription factors and score the TF-DNA binding at a base-pair resolution. We find significant differences in activations of DNA context features for bound and unbound sequences. In addition to standardized evaluation protocols, we offer outstanding interpretability that enables us to identify and annotate DNA sequence with possible elements that modulate TF-DNA binding. Also, differences in data processing have a huge influence on the overall model performance. Overall, the proposed framework allows for novel insights on the non-coding genetic elements and their role in facilitating a stable TF-DNA interaction.

SOX11 regulates SWI/SNF complex components as member of the adrenergic neuroblastoma core regulatory circuitry

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

Genetic Predisposition to Neuroblastoma Results from a Regulatory Polymorphism that Promotes the Adrenergic Cell State

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like "mesenchymal" cell state and a more differentiated sympathetic "adrenergic" cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional co-factor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism G/T that affects a G ATA motif in the first intron of LMO1. Here we show that wild-type zebrafish with the G ATA genotype develop adrenergic neuroblastoma, while knock-in of the protective T ATA allele at this locus reduces the penetrance of MYCN-driven tumors, which are restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrates that T ATA/ T ATA tumors also exhibit a mesenchymal cell state and are low risk at diagnosis. Thus, conversion of the regulatory G ATA to a T ATA allele in the first intron of LMO1 reduces the neuroblastoma initiation rate by preventing formation of the adrenergic cell state, a mechanism that is conserved over 400 million years of evolution separating zebrafish and humans.

Chaperonin containing TCP-1 (CCT/TRiC) is a novel therapeutic and diagnostic target for neuroblastoma

Chaperonin containing TCP1 (CCT/TRiC) is a multi-subunit protein folding complex that enables the cancer phenotype to emerge from the mutational landscape that drives oncogenesis. We and others linked increased expression of CCT subunits to advanced tumor stage and invasiveness that inversely correlates with cancer patient outcomes. In this study, we examined the expression of the second CCT subunit, CCT2, using genomic databases of adult and pediatric tumors and normal tissues, and found that it was highly expressed in pediatric cancers, showing a significant difference compared to normal tissues. Histologic staining confirmed that CCT subunits are highly expressed in tumor tissues, which was exemplified in neuroblastoma. Using two neuroblastoma cells, MYCN-amplified, IMR-32 cells, and non-amplified, SK-N-AS cells, we assessed baseline levels for CCT subunits and found expressions comparable to the highly invasive triple-negative breast cancer (TNBC) cell line, MDA-MB-231. Exogenous expression of CCT2 in both SK-N-AS and IMR-32 cells resulted in morphological changes, such as larger cell size and increased adherence, with significant increases in the CCT substrates, actin, and tubulin, as well as increased migration. Depletion of CCT2 reversed these effects and reduced cell viability. We evaluated CCT as a therapeutic target in IMR-32 cells by testing a novel peptide CCT inhibitor, CT20p. Treatment with CT20p induced cell death in these neuroblastoma cells. The use of CCT2 as a biological indicator for detection of neuroblastoma cells shed in blood was examined by spiking IMR-32 cells into human blood and using an anti-CCT2 antibody for the identification of spiked cancer cells with the CellSearch system. Results showed that using CCT2 for the detection of neuroblastoma cells in blood was more effective than the conventional approach of using epithelial markers like cytokeratins. CCT2 plays an essential role in promoting the invasive capacity of neuroblastoma cells and thus offers the potential to act as a molecular target in the development of novel therapeutics and diagnostics for pediatric cancers.