5-hydroxymethylcytosine profiling of cell-free DNA identifies bivalent genes that are prognostic of survival in high-risk neuroblastoma.

Neuroblastoma is the most common extra-cranial solid tumor in childhood and epigenetic dysregulation is a key driver of this embryonal disease. In cell-free DNA from neuroblastoma patients with high-risk disease, we found increased 5-hydroxymethylcytosine (5-hmC) deposition on Polycomb Repressive Complex 2 (PRC2) target genes, a finding previously described in the context of bivalent genes. As bivalent genes, defined as genes bearing both activating (H3K4me3) and repressive (H3K27me3) chromatin modifications, have been shown to play an important role in development and cancer, we investigated the potential role of bivalent genes in maintaining a de-differentiated state in neuroblastoma and their potential use as a biomarker. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogenous neuroblastoma cell lines. Through gene set variance analysis, we developed a clinically implementable bivalent signature. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in high-risk neuroblastoma patients. Thus, low expression of bivalent genes is a biomarker of ultra-high-risk disease and may represent a therapeutic opportunity in neuroblastoma.