Individual patient risk stratification of high-risk neuroblastomas using a two-gene score suited for clinical use

Intra-Tumour Genetic Heterogeneity and Prognosis in High-Risk Neuroblastoma

Spatial ITH is defined by genomic and biological variations within a tumour acquired by tumour cell evolution under diverse microenvironments, and its role in NB patient prognosis is understudied. In this work, we applied pangenomic techniques to detect chromosomal aberrations in at least two different areas of each tumour and/or in simultaneously obtained solid and liquid biopsies, detecting ITH in the genomic profile of almost 40% of HR-NB. ITH was better detected when comparing one or more tumour pieces and liquid biopsy (50%) than between different tumour pieces (21%). Interestingly, we found that patients with ITH analysed by pangenomic techniques had a significantly better survival rate that those with non-heterogeneous tumours, especially in cases without MYCN amplification. Moreover, all patients in the studied cohort with high ITH (defined as 50% or more genomic aberration differences between areas of a tumour or simultaneously obtained samples) survived after 48 months. These results clearly support analysing at least two solid tumour areas (separately or mixed) and liquid samples to provide more accurate genomic diagnosis, prognosis and therapy options in HR-NB.

Hypoxia Predicts Poor Prognosis in Neuroblastoma Patients and Associates with Biological Mechanisms Involved in Telomerase Activation and Tumor Microenvironment Reprogramming

The biological and clinical heterogeneity of neuroblastoma (NB) demands novel biomarkers and therapeutic targets in order to drive the most appropriate treatment for each patient. Hypoxia is a condition of low-oxygen tension occurring in poorly vascularized tumor tissues. In this study, we aimed to assess the role of hypoxia in the pathogenesis of NB and at developing a new clinically relevant hypoxia-based predictor of outcome. We analyzed the gene expression profiles of 1882 untreated NB primary tumors collected at diagnosis and belonging to four existing data sets. Analyses took advantage of machine learning methods. We identified NB-hop, a seven-gene hypoxia biomarker, as a predictor of NB patient prognosis, which is able to discriminate between two populations of patients with unfavorable or favorable outcome on a molecular basis. NB-hop retained its prognostic value in a multivariate model adjusted for established risk factors and was able to additionally stratify clinically relevant groups of patients. Tumors with an unfavorable NB-hop expression showed a significant association with telomerase activation and a hypoxic, immunosuppressive, poorly differentiated, and apoptosis-resistant tumor microenvironment. NB-hop defines a new population of NB patients with hypoxic tumors and unfavorable prognosis and it represents a critical factor for the stratification and treatment of NB patients.

The challenge of defining "ultra-high-risk" neuroblastoma

Given the biological and clinical heterogeneity of neuroblastoma, risk stratification is vital to determining appropriate treatment. Historically, most patients with high-risk neuroblastoma (HR-NBL) have been treated uniformly without further stratification. Attempts have been made to identify factors that can be used to risk stratify these patients and to characterize an "ultra-high-risk" (UHR) subpopulation with particularly poor outcome. However, among published data, there is a lack of consensus in the definition of the UHR population and heterogeneity in the endpoints and statistical methods used. This review summarizes our current understanding of stratification of HR-NBL and discusses the complex issues in defining UHR neuroblastoma.

Risk stratification of high-risk metastatic neuroblastoma: A report from the HR-NBL-1/SIOPEN study

BACKGROUND

Risk stratification is crucial to treatment decision-making in neuroblastoma. This study aimed to explore factors present at diagnosis affecting outcome in patients aged ≥18 months with metastatic neuroblastoma and to develop a simple risk score for prognostication.

PROCEDURE

Data were derived from the European high-risk neuroblastoma 1 (HR-NBL1)/International Society for Paediatric Oncology European Neuroblastoma (SIOPEN) trial with analysis restricted to patients aged ≥18 months with metastatic disease and treated prior to the introduction of immunotherapy. Primary endpoint was 5-year event-free survival (EFS). Prognostic factors assessed were sex, age, tumour MYCN amplification (MNA) status, serum lactate dehydrogenase (LDH)/ferritin, primary tumour and metastatic sites. Factors significant in univariate analysis were incorporated into a multi-variable model and an additive scoring system developed based on estimated log-cumulative hazard ratios.

RESULTS

The cohort included 1053 patients with median follow-up 5.5 years and EFS 27 ± 1%. In univariate analyses, age; serum LDH and ferritin; involvement of bone marrow, bone, liver or lung; and >1 metastatic system/compartment were associated with worse EFS. Tumour MNA was not associated with worse EFS. A multi-variable model and risk score incorporating age (>5 years, 2 points), serum LDH (>1250 U/L, 1 point) and number of metastatic systems (>1, 2 points) were developed. EFS was significantly correlated with risk score: EFS 52 ± 9% for score = 0 versus 6 ± 3% for score = 5 (P < 0.0001).

CONCLUSIONS

A simple score can identify an "ultra-high risk" (UHR) cohort (score = 5) comprising 8% of patients with 5-year EFS <10%. These patients appear not to benefit from induction therapy and could potentially be directed earlier to alternative experimental therapies in future trials.

Immunohistochemical analysis of PDK1, PHD3 and HIF-1α expression defines the hypoxic status of neuroblastoma tumors

Neuroblastoma (NB) is the most common solid tumor during infancy and the first cause of death among the preschool age diseases. The availability of several NB genomic profiles improves the prognostic ability, but the outcome prediction for this pathology remains imperfect. We previously produced a novel prognostic gene signature based on the response of NB cells to hypoxia, a condition of tumor microenvironment strictly connected with cancer aggressiveness. Here we attempted to further define the expression of hypoxia-modulated specific genes, looking at their protein level in NB specimens, considering in particular the hypoxia inducible factor-1α (HIF-1α), the mitochondrial pyruvate dehydrogenase kinase 1 (PDK1), and the HIF-prolyl hydroxylase domain 3 (PHD3). The evaluation of expression was performed by Western blot and immunocytochemistry on NB cell lines and by immunohistochemistry on tumor specimens. Stimulation of both HIF-1α and PDK1 and inhibition of PHD3 expression were observed in NB cell lines cultured under prolonged hypoxic conditions as well as in most of the tumors with poor outcome. Our results indicate that the immunohistochemistry analysis of the protein expression of PDK1, PHD3, and HIF-1α defines the hypoxic status of NB tumors and can be used as a simple and relevant tool to stratify high-risk patients.

Neuroblastoma patient-derived xenograft cells cultured in stem-cell promoting medium retain tumorigenic and metastatic capacities but differentiate in serum

Cultured cancer cells serve as important models for preclinical testing of anti-cancer compounds. However, the optimal conditions for retaining original tumor features during in vitro culturing of cancer cells have not been investigated in detail. Here we show that serum-free conditions are critical for maintaining an immature phenotype of neuroblastoma cells isolated from orthotopic patient-derived xenografts (PDXs). PDX cells could be grown either as spheres or adherent on laminin in serum-free conditions with retained patient-specific genomic aberrations as well as tumorigenic and metastatic capabilities. However, addition of serum led to morphological changes, neuronal differentiation and reduced cell proliferation. The epidermal growth factor (EGF) and basic fibroblast growth factor (bFGF) were central for PDX cell proliferation and MYCN expression, and also hindered the serum-induced differentiation. Although serum induced a robust expression of neurotrophin receptors, stimulation with their cognate ligands did not induce further sympathetic differentiation, which likely reflects a block in PDX cell differentiation capacity coupled to their tumor genotype. Finally, PDX cells cultured as spheres or adherent on laminin responded similarly to various cytotoxic drugs, suggesting that both conditions are suitable in vitro screening models for neuroblastoma-targeting compounds.

Artificial neural network classifier predicts neuroblastoma patients' outcome

Background

More than fifty percent of neuroblastoma (NB) patients with adverse prognosis do not benefit from treatment making the identification of new potential targets mandatory. Hypoxia is a condition of low oxygen tension, occurring in poorly vascularized tissues, which activates specific genes and contributes to the acquisition of the tumor aggressive phenotype. We defined a gene expression signature (NB-hypo), which measures the hypoxic status of the neuroblastoma tumor. We aimed at developing a classifier predicting neuroblastoma patients’ outcome based on the assessment of the adverse effects of tumor hypoxia on the progression of the disease.

Methods

Multi-layer perceptron (MLP) was trained on the expression values of the 62 probe sets constituting NB-hypo signature to develop a predictive model for neuroblastoma patients’ outcome. We utilized the expression data of 100 tumors in a leave-one-out analysis to select and construct the classifier and the expression data of the remaining 82 tumors to test the classifier performance in an external dataset. We utilized the Gene set enrichment analysis (GSEA) to evaluate the enrichment of hypoxia related gene sets in patients predicted with “Poor” or “Good” outcome.

Results

We utilized the expression of the 62 probe sets of the NB-Hypo signature in 182 neuroblastoma tumors to develop a MLP classifier predicting patients’ outcome (NB-hypo classifier). We trained and validated the classifier in a leave-one-out cross-validation analysis on 100 tumor gene expression profiles. We externally tested the resulting NB-hypo classifier on an independent 82 tumors’ set. The NB-hypo classifier predicted the patients’ outcome with the remarkable accuracy of 87 %. NB-hypo classifier prediction resulted in 2 % classification error when applied to clinically defined low-intermediate risk neuroblastoma patients. The prediction was 100 % accurate in assessing the death of five low/intermediated risk patients. GSEA of tumor gene expression profile demonstrated the hypoxic status of the tumor in patients with poor prognosis.

Conclusions

We developed a robust classifier predicting neuroblastoma patients’ outcome with a very low error rate and we provided independent evidence that the poor outcome patients had hypoxic tumors, supporting the potential of using hypoxia as target for neuroblastoma treatment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1194-3) contains supplementary material, which is available to authorized users.