Inhibiting phosphoglycerate dehydrogenase counteracts chemotherapeutic efficacy against MYCN-amplified neuroblastoma.

Here we sought metabolic alterations specifically associated with MYCN amplification as nodes to indirectly target the MYCN oncogene. Liquid chromatography-mass spectrometry-based proteomics identified seven proteins consistently correlated with MYCN in proteomes from 49 neuroblastoma biopsies and 13 cell lines. Among these was phosphoglycerate dehydrogenase (PHGDH), the rate-limiting enzyme in de novo serine synthesis. MYCN associated with two regions in the PHGDH promoter, supporting transcriptional PHGDH regulation by MYCN. Pulsed stable isotope-resolved metabolomics utilizing 13 C-glucose labeling demonstrated higher de novo serine synthesis in MYCN-amplified cells compared to cells with diploid MYCN. An independence of MYCN-amplified cells from exogenous serine and glycine was demonstrated by serine and glycine starvation, which attenuated nucleotide pools and proliferation only in cells with diploid MYCN but did not diminish these endpoints in MYCN-amplified cells. Proliferation was attenuated in MYCN-amplified cells by CRISPR/Cas9-mediated PHGDH knockout or treatment with PHGDH small molecule inhibitors without affecting cell viability. PHGDH inhibitors administered as single-agent therapy to NOG mice harboring patient-derived MYCN-amplified neuroblastoma xenografts slowed tumor growth. However, combining a PHGDH inhibitor with the standard-of-care chemotherapy drug, cisplatin, revealed antagonism of chemotherapy efficacy in vivo. Emergence of chemotherapy resistance was confirmed in the genetic PHGDH knockout model in vitro. Altogether, PHGDH knockout or inhibition by small molecules consistently slows proliferation, but stops short of killing the cells, which then establish resistance to classical chemotherapy. Although PHGDH inhibition with small molecules has produced encouraging results in other preclinical cancer models, this approach has limited attractiveness for patients with neuroblastoma.

Reflection of neuroblastoma intratumor heterogeneity in the new OHC-NB1 disease model.

Accurate modeling of intratumor heterogeneity presents a bottleneck against drug testing. Flexibility in a preclinical platform is also desirable to support assessment of different endpoints. We established the model system, OHC-NB1, from a bone marrow metastasis from a patient diagnosed with MYCN-amplified neuroblastoma and performed whole-exome sequencing on the source metastasis and the different models and passages during model development (monolayer cell line, 3D spheroid culture and subcutaneous xenograft tumors propagated in mice). OHC-NB1 harbors a MYCN amplification in double minutes, 1p deletion, 17q gain and diploid karyotype, which persisted in all models. A total of 80-540 single-nucleotide variants (SNVs) was detected in each sample, and comparisons between the source metastasis and models identified 34 of 80 somatic SNVs to be propagated in the models. Clonal reconstruction using the combined copy number and SNV data revealed marked clonal heterogeneity in the originating metastasis, with four clones being reflected in the model systems. The set of OHC-NB1 models represents 43% of somatic SNVs and 23% of the cellularity in the originating metastasis with varying clonal compositions, indicating that heterogeneity is partially preserved in our model system.

Minichromosome Maintenance Complex Is a Critical Node in the miR-183 Signaling Network of MYCN-Amplified Neuroblastoma Cells.

MYCN and HDAC2 jointly repress the transcription of tumor suppressive miR-183 in neuroblastoma. Enforced miR-183 expression induces neuroblastoma cell death and inhibits xenograft growth in mice. Here we aimed to focus more closely on the miR-183 signaling network using a label-free mass spectrometric approach. Analysis of neuroblastoma cells transfected with either control or miR-183 expression vectors identified 85 differentially expressed proteins. All six members of the minichromosome maintenance (MCM) complex, which is indispensable for initiation and elongation during DNA replication and transcriptionally activated by MYCN in neuroblastoma, emerged to be down-regulated by miR-183. Subsequent annotation category enrichment analysis revealed a ∼14-fold enrichment in the "MCM" protein module category, which highlighted this complex as a critical node in the miR-183 signaling network. Down-regulation was confirmed by Western blotting. MCMs 2-5 were predicted by in silico methods as direct miR-183 targets. Dual-luciferase reporter gene assays with 3'-UTR constructs of the randomly selected MCMs 3 and 5 experimentally confirmed them as direct targets of miR-183. Our results reveal the MCM complex to be a critical and directly regulated node within the miR-183 signaling network in MYCN-amplified neuroblastoma cells.

Histone deacetylase 10 promotes autophagy-mediated cell survival.

Tumor cells activate autophagy in response to chemotherapy-induced DNA damage as a survival program to cope with metabolic stress. Here, we provide in vitro and in vivo evidence that histone deacetylase (HDAC)10 promotes autophagy-mediated survival in neuroblastoma cells. We show that both knockdown and inhibition of HDAC10 effectively disrupted autophagy associated with sensitization to cytotoxic drug treatment in a panel of highly malignant V-MYC myelocytomatosis viral-related oncogene, neuroblastoma derived-amplified neuroblastoma cell lines, in contrast to nontransformed cells. HDAC10 depletion in neuroblastoma cells interrupted autophagic flux and induced accumulation of autophagosomes, lysosomes, and a prominent substrate of the autophagic degradation pathway, p62/sequestosome 1. Enforced HDAC10 expression protected neuroblastoma cells against doxorubicin treatment through interaction with heat shock protein 70 family proteins, causing their deacetylation. Conversely, heat shock protein 70/heat shock cognate 70 was acetylated in HDAC10-depleted cells. HDAC10 expression levels in high-risk neuroblastomas correlated with autophagy in gene-set analysis and predicted treatment success in patients with advanced stage 4 neuroblastomas. Our results demonstrate that HDAC10 protects cancer cells from cytotoxic agents by mediating autophagy and identify this HDAC isozyme as a druggable regulator of advanced-stage tumor cell survival. Moreover, these results propose a promising way to considerably improve treatment response in the neuroblastoma patient subgroup with the poorest outcome.

MYCN and HDAC2 cooperate to repress miR-183 signaling in neuroblastoma.

MYCN is a master regulator controlling many processes necessary for tumor cell survival. Here, we unravel a microRNA network that causes tumor suppressive effects in MYCN-amplified neuroblastoma cells. In profiling studies, histone deacetylase (HDAC) inhibitor treatment most strongly induced miR-183. Enforced miR-183 expression triggered apoptosis, and inhibited anchorage-independent colony formation in vitro and xenograft growth in mice. Furthermore, the mechanism of miR-183 induction was found to contribute to the cell death phenotype induced by HDAC inhibitors. Experiments to identify the HDAC(s) involved in miR-183 transcriptional regulation showed that HDAC2 depletion induced miR-183. HDAC2 overexpression reduced miR-183 levels and counteracted the induction caused by HDAC2 depletion or HDAC inhibitor treatment. MYCN was found to recruit HDAC2 in the same complexes to the miR-183 promoter, and HDAC2 depletion enhanced promoter-associated histone H4 pan-acetylation, suggesting epigenetic changes preceded transcriptional activation. These data reveal miR-183 tumor suppressive properties in neuroblastoma that are jointly repressed by MYCN and HDAC2, and suggest a novel way to bypass MYCN function.

Management of Busulfan-Induced Lung Injury in Pediatric Patients with High-Risk Neuroblastoma.

Background/Objectives: Integrating the cytotoxic drug busulfan into a high-dose chemotherapy regimen prior to autologous hematopoietic stem cell rescue in patients with high-risk neuroblastoma has improved the survival of children battling this deadly disease. Busulfan-induced toxicities can, however, be severe. Here, we describe the diagnosis and successful treatment of acute pulmonary injury by total-body-weight-adjusted busulfan therapy in two children with high-risk neuroblastoma. Case series: Patient 1 developed life-threatening biphasic acute respiratory failure on days +60 and +100 after busulfan therapy, requiring intubation and invasive mechanical ventilation. Despite intensive anti-inflammatory and immunomodulatory therapy, including systemic corticosteroids, topical inhalation regimens, azithromycin, nintedanib and extracorporal photopheresis, patient 1 required extended intensive care measures and non-invasive respiratory support for a total of 20 months. High-resolution computed tomography showed diffuse intra-alveolar and interstitial patterns. Patient 2 developed partial respiratory failure with insufficient oxygen saturation and dyspnea on day +52 after busulfan therapy. Symptoms were resolved after 6 months of systemic corticosteroids, topical inhalation regimens and azithromycin. High-resolution computed tomography showed atypical pneumonic changes with ground-glass opacities. While both patients fully recovered without evidence of pulmonary fibrosis, cancer therapy had to be paused and then modified until full recovery from busulfan-induced lung injury. Conclusions: Busulfan-induced lung injury requires prompt diagnosis and intervention. Symptoms and signs are nonspecific and difficult to differentiate from other causes. Therapeutic busulfan drug level monitoring and the identification of patients at risk for drug overdosing through promoter polymorphisms in the glutathione S-transferase alpha 1 gene encoding the main enzyme in busulfan metabolism are expected to reduce the risk of busulfan-induced toxicities.

Targeting KDM1A in Neuroblastoma with NCL-1 Induces a Less Aggressive Phenotype and Suppresses Angiogenesis.

Background: The KDM1A histone demethylase regulates the cellular balance between proliferation and differentiation, and is often deregulated in human cancers including the childhood tumor neuroblastoma. We previously showed that KDM1A is strongly expressed in undifferentiated neuroblastomas and correlates with poor patient prognosis, suggesting a possible clinical benefit from targeting KDM1A. Methods: Here, we tested the efficacy of NCL-1, a small molecule specifically inhibiting KDM1A in preclinical models for neuroblastoma. Results: NCL-1 mimicked the effects of siRNA-mediated KDM1A knockdown and effectively inhibited KDM1A activity in four neuroblastoma cell lines and a patient-representative cell model. KDM1A inhibition shifted the aggressive tumor cell phenotypes towards less aggressive phenotypes. The proliferation and cell viability was reduced, accompanied by the induction of markers of neuronal differentiation. Interventional NCL-1 treatment of nude mice harboring established neuroblastoma xenograft tumors reduced tumor growth and inhibited cell proliferation. Reduced vessel density and defects in blood vessel construction also resulted, and NCL-1 inhibited the growth and tube formation of HUVEC-C cells in vitro. Conclusions: Inhibiting KDM1A could attack aggressive neuroblastomas two-fold, by re-directing tumor cells toward a less aggressive, slower-growing phenotype and by preventing or reducing the vascular support of large tumors.

HDAC11 is a novel drug target in carcinomas.

Inhibition of histone deacetylase (HDAC) activity as stand-alone or combination therapy represents a promising therapeutic approach in oncology. The pan- or class I HDAC inhibitors (HDACi) currently approved or in clinical studies for oncology give rise to dose-limiting toxicities, presumably because of the inhibition of several HDACs. This could potentially be overcome by selective blockade of single HDAC family members. Here we report that HDAC11, the most recently identified zinc-dependent HDAC, is overexpressed in several carcinomas as compared to corresponding healthy tissues. HDAC11 depletion is sufficient to cause cell death and to inhibit metabolic activity in HCT-116 colon, PC-3 prostate, MCF-7 breast and SK-OV-3 ovarian cancer cell lines. The antitumoral effect induced can be mimicked by enforced expression of a catalytically impaired HDAC11 variant, suggesting that inhibition of the enzymatic activity of HDAC11 by small molecules could trigger the desired phenotypic changes. HDAC11 depletion in normal cells causes no changes in metabolic activity and viability, strongly suggesting that tumor-selective effects can be achieved. Altogether, our data show that HDAC11 plays a critical role in cancer cell survival and may represent a novel drug target in oncology.

Parallel sequencing of extrachromosomal circular DNAs and transcriptomes in single cancer cells.

Extrachromosomal DNAs (ecDNAs) are common in cancer, but many questions about their origin, structural dynamics and impact on intratumor heterogeneity are still unresolved. Here we describe single-cell extrachromosomal circular DNA and transcriptome sequencing (scEC&T-seq), a method for parallel sequencing of circular DNAs and full-length mRNA from single cells. By applying scEC&T-seq to cancer cells, we describe intercellular differences in ecDNA content while investigating their structural heterogeneity and transcriptional impact. Oncogene-containing ecDNAs were clonally present in cancer cells and drove intercellular oncogene expression differences. In contrast, other small circular DNAs were exclusive to individual cells, indicating differences in their selection and propagation. Intercellular differences in ecDNA structure pointed to circular recombination as a mechanism of ecDNA evolution. These results demonstrate scEC&T-seq as an approach to systematically characterize both small and large circular DNA in cancer cells, which will facilitate the analysis of these DNA elements in cancer and beyond.

HD-MB03 is a novel Group 3 medulloblastoma model demonstrating sensitivity to histone deacetylase inhibitor treatment.

Medulloblastomas are the most common malignant brain tumors in childhood. Emerging evidence suggests that medulloblastoma comprises at least four distinct diseases (WNT, SHH, Group 3 and 4) with different biology, clinical presentation, and outcome, with especially poor prognosis in Group 3. The tight connection of biology and clinical behavior in patients emphasizes the need for subgroup-specific preclinical models in order to develop treatments tailored to each subgroup. Herein we report on the novel cell line HD-MB03, isolated from tumor material of a patient with metastasized Group 3 medulloblastoma, and preclinical testing of different histone deacetylase inhibitors (HDACis) in this model. HD-MB03 cells grow long term in vitro and form metastatic tumors in vivo upon orthotopic transplantation. HD-MB03 cells reflect the original Group 3 medulloblastoma at the histological and molecular level, showing large cell morphology, similar expression patterns for markers Ki67, p53, and glial fibrillary acidic protein (GFAP), a gene expression profile most closely matching Group 3 medulloblastomas, and persistence of typical molecular alterations, i.e., isochromosome 17q [i(17q)] and MYC amplification. Protein expression analysis of HDACs 2, 5, 8, and 9 as well as the predictive marker HR23B showed intermediate to strong expression, suggesting sensitivity to HDACis. Indeed, treatment with HDACis Helminthosporium carbonum (HC)-toxin, vorinostat, and panobinostat revealed high sensitivity to this novel drug class, as well as a radiation-sensitizing effect with significantly increased cell death upon concomitant treatment. In summary, our data indicate that HD-MB03 is a suitable preclinical model for Group 3 medulloblastoma, and HDACis could represent a therapeutic option for this subgroup.

A Deep Dive into the Circulating ctDNA Cosmos to Vanquish Neuroblastoma.

SUMMARY: Single biopsies fail to reflect intratumor heterogeneity and tumor evolution. In this issue of Cancer Discovery, Bosse and colleagues show an important role for circulating cell-free tumor DNA sequencing to detect the genomic evolution of neuroblastoma under ALK inhibitor therapy and identify novel (sub)clonal pathogenic variants involved in disease progression under conventional therapy. See related article by Bosse et al., p. 2800 (5).

Circulating Cell-Free DNA Assessment in Biofluids from Children with Neuroblastoma Demonstrates Feasibility and Potential for Minimally Invasive Molecular Diagnostics.

Liquid biopsy strategies in pediatric patients are challenging due to low body weight. This study investigated cfDNA size distribution and concentration in blood, bone marrow, cerebrospinal fluid, and urine from 84 patients with neuroblastoma classified as low (n = 28), intermediate (n = 6), or high risk (n = 50) to provide key data for liquid biopsy biobanking strategies. The average volume of blood and bone marrow plasma provided ranged between 1 and 2 mL. Analysis of 637 DNA electropherograms obtained by Agilent TapeStation measurement revealed five different major profiles and characteristic DNA size distribution patterns for each of the biofluids. The proportion of samples containing primarily cfDNA was, at 85.5%, the highest for blood plasma. The median cfDNA concentration amounted to 6.28 ng/mL (blood plasma), 58.2 ng/mL (bone marrow plasma), 0.08 ng/mL (cerebrospinal fluid), and 0.49 ng/mL (urine) in samples. Meta-analysis of the dataset demonstrated that multiple cfDNA-based assays employing the same biofluid sample optimally require sampling volumes of 1 mL for blood and bone marrow plasma, 2 mL for cerebrospinal fluid, and as large as possible for urine samples. A favorable response to treatment was associated with a rapid decrease in blood-based cfDNA concentration in patients with high-risk neuroblastoma. Blood-based cfDNA concentration was not sufficient as a single parameter to indicate high-risk disease recurrence. We provide proof of concept that monitoring neuroblastoma-specific markers in very small blood volumes from infants is feasible.

Targeted Analysis of Cell-free Circulating Tumor DNA is Suitable for Early Relapse and Actionable Target Detection in Patients with Neuroblastoma.

PURPOSE: Treating refractory or relapsed neuroblastoma remains challenging. Monitoring body fluids for tumor-derived molecular information indicating minimal residual disease supports more frequent diagnostic surveillance and may have the power to detect resistant subclones before they give rise to relapses. If actionable targets are identified from liquid biopsies, targeted treatment options can be considered earlier. EXPERIMENTAL DESIGN: Droplet digital PCR assays assessing MYCN and ALK copy numbers and allelic frequencies of ALK p.F1174L and ALK p.R1275Q mutations were applied to longitudinally collected liquid biopsies and matched tumor tissue samples from 31 patients with high-risk neuroblastoma. Total cell-free DNA (cfDNA) levels and marker detection were compared with data from routine clinical diagnostics. RESULTS: Total cfDNA concentrations in blood plasma from patients with high-risk neuroblastoma were higher than in healthy controls and consistently correlated with neuron-specific enolase levels and lactate dehydrogenase activity but not with 123I-meta-iodobenzylguanidine scores at relapse diagnosis. Targeted cfDNA diagnostics proved superior for early relapse detection to all current diagnostics in 2 patients. Marker analysis in cfDNA indicated intratumor heterogeneity for cell clones harboring MYCN amplifications and druggable ALK alterations that were not detectable in matched tumor tissue samples in 17 patients from our cohort. Proof of concept is provided for molecular target detection in cerebrospinal fluid from patients with isolated central nervous system relapses. CONCLUSIONS: Tumor-specific alterations can be identified and monitored during disease course in liquid biopsies from pediatric patients with high-risk neuroblastoma. This approach to cfDNA surveillance warrants further prospective validation and exploitation for diagnostic purposes and to guide therapeutic decisions.

A novel human high-risk ependymoma stem cell model reveals the differentiation-inducing potential of the histone deacetylase inhibitor Vorinostat.

Incompletely resectable ependymomas are associated with poor prognosis despite intensive radio- and chemotherapy. Novel treatments have been difficult to develop due to the lack of appropriate models. Here, we report on the generation of a high-risk cytogenetic group 3 and molecular group C ependymoma model (DKFZ-EP1NS) which is based on primary ependymoma cells obtained from a patient with metastatic disease. This model displays stem cell features such as self-renewal capacity, differentiation capacity, and specific marker expression. In vivo transplantation showed high tumorigenic potential of these cells, and xenografts phenotypically recapitulated the original tumor in a niche-dependent manner. DKFZ-EP1NS cells harbor transcriptome plasticity, enabling a shift from a neural stem cell-like program towards a profile of primary ependymoma tumor upon in vivo transplantation. Serial transplantation of DKFZ-EP1NS cells from orthotopic xenografts yielded secondary tumors in half the time compared with the initial transplantation. The cells were resistant to temozolomide, vincristine, and cisplatin, but responded to histone deacetylase inhibitor (HDACi) treatment at therapeutically achievable concentrations. In vitro treatment of DKFZ-EP1NS cells with the HDACi Vorinostat induced neuronal differentiation associated with loss of stem cell-specific properties. In summary, this is the first ependymoma model of a cytogenetic group 3 and molecular subgroup C ependymoma based on a human cell line with stem cell-like properties, which we used to demonstrate the differentiation-inducing therapeutic potential of HDACi.

Spatial and temporal intratumour heterogeneity has potential consequences for single biopsy-based neuroblastoma treatment decisions.

Intratumour heterogeneity is a major cause of treatment failure in cancer. We present in-depth analyses combining transcriptomic and genomic profiling with ultra-deep targeted sequencing of multiregional biopsies in 10 patients with neuroblastoma, a devastating childhood tumour. We observe high spatial and temporal heterogeneity in somatic mutations and somatic copy-number alterations which are reflected on the transcriptomic level. Mutations in some druggable target genes including ALK and FGFR1 are heterogeneous at diagnosis and/or relapse, raising the issue whether current target prioritization and molecular risk stratification procedures in single biopsies are sufficiently reliable for therapy decisions. The genetic heterogeneity in gene mutations and chromosome aberrations observed in deep analyses from patient courses suggest clonal evolution before treatment and under treatment pressure, and support early emergence of metastatic clones and ongoing chromosomal instability during disease evolution. We report continuous clonal evolution on mutational and copy number levels in neuroblastoma, and detail its implications for therapy selection, risk stratification and therapy resistance.

P160/SRC/NCoA coactivators form complexes via specific interaction of their PAS-B domain with the CID/AD1 domain.

Transcriptional activation involves the ordered recruitment of coactivators via direct interactions between distinct binding domains and recognition motifs. The p160/SRC/NCoA coactivator family comprises three members (NCoA-1, -2 and -3), which are organized in multiprotein coactivator complexes. We had identified the PAS-B domain of NCoA-1 as an LXXLL motif binding domain. Here we show that NCoA family members are able to interact with other full-length NCoA proteins via their PAS-B domain and they specifically interact with the CBP-interaction domain (CID/AD1) of NCoA-1. Peptide competition, binding experiments and mutagenesis of LXXLL motifs point at distinct binding motif specificities of the NCoA PAS-B domains. NMR studies of different NCoA-1-PAS-B/LXXLL peptide complexes revealed similar although not identical binding sites for the CID/AD1 and STAT6 transactivation domain LXXLL motifs. In mechanistic studies, we found that overexpression of the PAS-B domain is able to disturb the binding of NCoA-1 to CBP in cells and that a CID/AD1 peptide competes with STAT6 for NCoA-1 in vitro. Moreover, the expression of an endogenous androgen receptor target gene is affected by the overexpression of the NCoA-1 or NCoA-3 PAS-B domains. Our study discloses a new, complementary mechanism for the current model of coactivator recruitment to target gene promoters.

Multiplexed Quantification of Four Neuroblastoma DNA Targets in a Single Droplet Digital PCR Reaction.

The detection and characterization of cell-free DNA (cfDNA) in peripheral blood from neuroblastoma patients may serve as a minimally invasive approach to liquid biopsy. Major challenges in the analysis of cfDNA purified from blood samples are small sample volumes and low cfDNA concentrations. Droplet digital PCR (ddPCR) is a technology suitable for analyzing low levels of cfDNA. Reported here are two quadruplexed ddPCR assay protocols that reliably quantify MYCN and ALK copy numbers in a single reaction together with the two reference genes, NAGK and AFF3, and accurately estimate ALKF1174L (exon 23 position 3522, C>A) and ALKR1275Q (exon 25 position 3824, G>A) mutant allele fractions using cfDNA as input. The separation of positive and negative droplets was optimized for detecting two targets in each ddPCR fluorescence channel by the adjustment of the probe and primer concentrations of each target molecule. The quadruplexed assays were validated using a panel of 10 neuroblastoma cell lines and paired blood plasma and primary neuroblastoma samples from nine patients. Accuracy and sensitivity thresholds in quadruplexed assays corresponded well with those from the respective duplexed assays. Presented are two robust quadruplexed ddPCR protocols applicable in the routine clinical setting and that require only minimal plasma volumes for the assessment of MYCN and ALK oncogene status.

Using droplet digital PCR to analyze MYCN and ALK copy number in plasma from patients with neuroblastoma.

The invasive nature of surgical biopsies deters sequential application, and single biopsies often fail to reflect tumor dynamics, intratumor heterogeneity and drug sensitivities likely to change during tumor evolution and treatment. Implementing molecular characterization of cell-free neuroblastoma-derived DNA isolated from blood plasma could improve disease assessment for treatment selection and monitoring of patients with high-risk neuroblastoma. We established droplet digital PCR (ddPCR) protocols for MYCN and ALK copy number status in plasma from neuroblastoma patients. Our ddPCR protocol accurately discriminated between MYCN and ALK amplification, gain and normal diploid status in a large panel of neuroblastoma cell lines, and discrepancies with reported MYCN and ALK status were detected, including a high-level MYCN amplification in NB-1, a MYCN gain in SH-SY5Y, a high-level ALK amplification in IMR-32 and ALK gains in BE(2)-C, Kelly, SH-SY5Y and LAN-6. MYCN and ALK status were also reliably determined from cell-free DNA derived from medium conditioned by the cell lines. MYCN and ALK copy numbers of subcutaneous neuroblastoma xenograft tumors were accurately determined from cell-free DNA in the mouse blood plasma. In a final validation step, we accurately distinguished MYCN and ALK copy numbers of the corresponding primary tumors using retrospectively collected blood plasma samples from 10 neuroblastoma patients. Our data justify the further development of molecular disease characterization using cell-free DNA in blood plasma from patients with neuroblastoma. This expanded molecular diagnostic palette may improve monitoring of disease progression including relapse and metastatic events as well as therapy success or failure in high-risk neuroblastoma patients.

Neuroblastoma cells depend on HDAC11 for mitotic cell cycle progression and survival.

The number of long-term survivors of high-risk neuroblastoma remains discouraging, with 10-year survival as low as 20%, despite decades of considerable international efforts to improve outcome. Major obstacles remain and include managing resistance to induction therapy, which causes tumor progression and early death in high-risk patients, and managing chemotherapy-resistant relapses, which can occur years after the initial diagnosis. Identifying and validating novel therapeutic targets is essential to improve treatment. Delineating and deciphering specific functions of single histone deacetylases in neuroblastoma may support development of targeted acetylome-modifying therapeutics for patients with molecularly defined high-risk neuroblastoma profiles. We show here that HDAC11 depletion in MYCN-driven neuroblastoma cell lines strongly induces cell death, mostly mediated by apoptotic programs. Genes necessary for mitotic cell cycle progression and cell division were most prominently enriched in at least two of three time points in whole-genome expression data combined from two cell systems, and all nine genes in these functional categories were strongly repressed, including CENPA, KIF14, KIF23 and RACGAP1. Enforced expression of one selected candidate, RACGAP1, partially rescued the induction of apoptosis caused by HDAC11 depletion. High-level expression of all nine genes in primary neuroblastomas significantly correlated with unfavorable overall and event-free survival in patients, suggesting a role in mediating the more aggressive biological and clinical phenotype of these tumors. Our study identified a group of cell cycle-promoting genes regulated by HDAC11, being both predictors of unfavorable patient outcome and essential for tumor cell viability. The data indicate a significant role of HDAC11 for mitotic cell cycle progression and survival of MYCN-amplified neuroblastoma cells, and suggests that HDAC11 could be a valuable drug target.

MYCN and HDAC5 transcriptionally repress CD9 to trigger invasion and metastasis in neuroblastoma.

The systemic and resistant nature of metastatic neuroblastoma renders it largely incurable with current multimodal treatment. Clinical progression stems mainly from the increasing burden of metastatic colonization. Therapeutically inhibiting the migration-invasion-metastasis cascade would be of great benefit, but the mechanisms driving this cycle are as yet poorly understood. In-depth transcriptome analyses and ChIP-qPCR identified the cell surface glycoprotein, CD9, as a major downstream player and direct target of the recently described GRHL1 tumor suppressor. CD9 is known to block or facilitate cancer cell motility and metastasis dependent upon entity. High-level CD9 expression in primary neuroblastomas correlated with patient survival and established markers for favorable disease. Low-level CD9 expression was an independent risk factor for adverse outcome. MYCN and HDAC5 colocalized to the CD9 promoter and repressed transcription. CD9 expression diminished with progressive tumor development in the TH-MYCN transgenic mouse model for neuroblastoma, and CD9 expression in neuroblastic tumors was far below that in ganglia from wildtype mice. Primary neuroblastomas lacking MYCN amplifications displayed differential CD9 promoter methylation in methyl-CpG-binding domain sequencing analyses, and high-level methylation was associated with advanced stage disease, supporting epigenetic regulation. Inducing CD9 expression in a SH-EP cell model inhibited migration and invasion in Boyden chamber assays. Enforced CD9 expression in neuroblastoma cells transplanted onto chicken chorioallantoic membranes strongly reduced metastasis to embryonic bone marrow. Combined treatment of neuroblastoma cells with HDAC/DNA methyltransferase inhibitors synergistically induced CD9 expression despite hypoxic, metabolic or cytotoxic stress. Our results show CD9 is a critical and indirectly druggable suppressor of the invasion-metastasis cycle in neuroblastoma.