5-Hydroxymethylcytosine Profiles in Circulating Cell-Free DNA Associate with Disease Burden in Children with Neuroblastoma.

PURPOSE: 5-Hydroxymethylcytosine (5-hmC) is an epigenetic marker of open chromatin and active gene expression. We profiled 5-hmC with Nano-hmC-Seal technology using 10 ng of plasma-derived cell-free DNA (cfDNA) in blood samples from patients with neuroblastoma to determine its utility as a biomarker. EXPERIMENTAL DESIGN: For the Discovery cohort, 100 5-hmC profiles were generated from 34 well children and 32 patients (27 high-risk, 2 intermediate-risk, and 3 low-risk) at various time points during the course of their disease. An independent Validation cohort encompassed 5-hmC cfDNA profiles (n = 29) generated from 21 patients (20 high-risk and 1 intermediate-risk). Metastatic burden was classified as high, moderate, low, or none per Curie metaiodobenzylguanidine scores and percentage of tumor cells in bone marrow. Genes with differential 5-hmC levels between samples according to metastatic burden were identified using DESeq2. RESULTS: Hierarchical clustering using 5-hmC levels of 347 genes identified from the Discovery cohort defined four clusters of samples that were confirmed in the Validation cohort and corresponded to high, high-moderate, moderate, and low/no metastatic burden. Samples from patients with increased metastatic burden had increased 5-hmC deposition on genes in neuronal stem cell maintenance and epigenetic regulatory pathways. Further, 5-hmC cfDNA profiles generated with 1,242 neuronal pathway genes were associated with subsequent relapse in the cluster of patients with predominantly low or no metastatic burden (sensitivity 65%, specificity 75.6%). CONCLUSIONS: cfDNA 5-hmC profiles in children with neuroblastoma correlate with metastatic burden and warrants development as a biomarker of treatment response and outcome.

Maternal Embryonic Leucine Zipper Kinase (MELK), a Potential Therapeutic Target for Neuroblastoma.

Maternal embryonic leucine zipper kinase (MELK) activates pathways that mediate aggressive tumor growth and therapy resistance in many types of adult cancers. Pharmacologic and genomic inhibition of MELK impairs tumor growth and increases sensitivity to radiation and chemotherapy. On the basis of these promising preclinical studies, early-phase adult clinical trials testing the MELK inhibitor OTS167 are ongoing. To investigate whether MELK is also a therapeutic target in neuroblastoma, we analyzed MELK expression in primary tumors and cell lines, and examined the effects of OTS167 on neuroblastoma growth. In primary tumors, high levels of MELK were associated with advanced stage disease and inferior survival. Higher levels of MELK were also detected in tumorigenic versus nontumorigenic neuroblastoma cell lines, and cells with higher levels of MELK expression were more sensitive to OTS167 than low-MELK expressing cells. OTS167 suppressed the growth of neuroblastoma xenografts, and in a preclinical model of minimal residual disease, survival was prolonged with MELK inhibition. OTS167 treatment downregulated MELK and its target enhancer of zeste homolog 2 (EZH2), a component of the polycomb repressive complex 2 (PRC2) that is known to modulate the DNA damage response. We also show that OTS167 reduced the formation of collapsed replication forks induced by camptothecin or radiation. Taken together, our results indicate that MELK indirectly mediates efficient processing of replication-associated DNA lesions in neuroblastoma, and that OTS167 sensitizes cells to DNA-damaging agents by abrogating this process. Further studies evaluating the activity of combination treatment regimens with OTS167 in neuroblastoma are warranted.

5-Hydroxymethylcytosine Profiles Are Prognostic of Outcome in Neuroblastoma and Reveal Transcriptional Networks That Correlate With Tumor Phenotype.

PURPOSE: Whole-genome profiles of the epigenetic modification 5-hydroxymethylcytosine (5-hmC) are robust diagnostic biomarkers in adult patients with cancer. We investigated if 5-hmC profiles would serve as novel prognostic markers in neuroblastoma, a clinically heterogeneous pediatric cancer. Because this DNA modification facilitates active gene expression, we hypothesized that 5-hmC profiles would identify transcriptomic networks driving the clinical behavior of neuroblastoma. PATIENTS AND METHODS: Nano-hmC-Seal sequencing was performed on DNA from Discovery (n = 51), Validation (n = 38), and Children's Oncology Group (n = 20) cohorts of neuroblastoma tumors. RNA was isolated from 48 tumors for RNA sequencing. Genes with differential 5-hmC or expression between clusters were identified using DESeq2. A 5-hmC model predicting outcome in high-risk patients was established using linear discriminant analysis. RESULTS: Comparison of low- versus high-risk tumors in the Discovery cohort revealed 577 genes with differential 5-hmC. Hierarchical clustering of tumors from the Discovery and Validation cohorts using these genes identified two main clusters highly associated with established prognostic markers, clinical risk group, and outcome. Genes with increased 5-hmC and expression in the favorable cluster were enriched for pathways of neuronal differentiation and KRAS activation, whereas genes involved in inflammation and the PRC2 complex were identified in the unfavorable cluster. The linear discriminant analysis model trained on high-risk Discovery cohort tumors was prognostic of outcome when applied to high-risk tumors from the Validation and Children's Oncology Group cohorts (hazard ratio, 3.8). CONCLUSION: 5-hmC profiles may be optimal DNA-based biomarkers in neuroblastoma. Analysis of transcriptional networks regulated by these epigenomic modifications may lead to a deeper understanding of drivers of neuroblastoma phenotype.

Secreted protein acidic and rich in cysteine (SPARC) induces lipotoxicity in neuroblastoma by regulating transport of albumin complexed with fatty acids.

SPARC is a matrix protein that mediates interactions between cells and the microenvironment. In cancer, SPARC may either promote or inhibit tumor growth depending upon the tumor type. In neuroblastoma, SPARC is expressed in the stromal Schwannian cells and functions as a tumor suppressor. Here, we developed a novel in vivo model of stroma-rich neuroblastoma using non-tumorigenic SHEP cells with modulated levels of SPARC, mixed with tumorigenic KCNR cells. Tumors with stroma-derived SPARC displayed suppressed growth, inhibited angiogenesis and increased lipid accumulation. Based on the described chaperone function of SPARC, we hypothesized that SPARC binds albumin complexed with fatty acids and transports them to tumors. We show that SPARC binds albumin with Kd=18.9±2.3 uM, and enhances endothelial cell internalization and transendothelial transport of albumin in vitro. We also demonstrate that lipids induce toxicity in neuroblastoma cells and show that lipotoxicity is increased when cells are cultured in hypoxic conditions. Studies investigating the therapeutic potential of SPARC are warranted.

Integrative genomics reveals hypoxia inducible genes that are associated with a poor prognosis in neuroblastoma patients.

Neuroblastoma is notable for its broad spectrum of clinical behavior ranging from spontaneous regression to rapidly progressive disease. Hypoxia is well known to confer a more aggressive phenotype in neuroblastoma. We analyzed transcriptome data from diagnostic neuroblastoma tumors and hypoxic neuroblastoma cell lines to identify genes whose expression levels correlate with poor patient outcome and are involved in the hypoxia response. By integrating a diverse set of transcriptome datasets, including those from neuroblastoma patients and neuroblastoma derived cell lines, we identified nine genes (SLCO4A1, ENO1, HK2, PGK1, MTFP1, HILPDA, VKORC1, TPI1, and HIST1H1C) that are up-regulated in hypoxia and whose expression levels are correlated with poor patient outcome in three independent neuroblastoma cohorts. Analysis of 5-hydroxymethylcytosine and ENCODE data indicate that at least five of these nine genes have an increase in 5-hydroxymethylcytosine and a more open chromatin structure in hypoxia versus normoxia and are putative targets of hypoxia inducible factor (HIF) as they contain HIF binding sites in their regulatory regions. Four of these genes are key components of the glycolytic pathway and another three are directly involved in cellular metabolism. We experimentally validated our computational findings demonstrating that seven of the nine genes are significantly up-regulated in response to hypoxia in the four neuroblastoma cell lines tested. This compact and robustly validated group of genes, is associated with the hypoxia response in aggressive neuroblastoma and may represent a novel target for biomarker and therapeutic development.