Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1.

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

Integrative genomic analyses identify neuroblastoma risk genes involved in neuronal differentiation.

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes.

From the identification of actionable molecular targets to the generation of faithful neuroblastoma patient-derived preclinical models.

BACKGROUND: Neuroblastoma (NB) represents the most frequent and aggressive form of extracranial solid tumor of infants. Although the overall survival of patients with NB has improved in the last years, more than 50% of high-risk patients still undergo a relapse. Thus, in the era of precision/personalized medicine, the need for high-risk NB patient-specific therapies is urgent. METHODS: Within the PeRsonalizEd Medicine (PREME) program, patient-derived NB tumors and bone marrow (BM)-infiltrating NB cells, derived from either iliac crests or tumor bone lesions, underwent to histological and to flow cytometry immunophenotyping, respectively. BM samples containing a NB cells infiltration from 1 to 50 percent, underwent to a subsequent NB cells enrichment using immune-magnetic manipulation. Then, NB samples were used for the identification of actionable targets and for the generation of 3D/tumor-spheres and Patient-Derived Xenografts (PDX) and Cell PDX (CPDX) preclinical models. RESULTS: Eighty-four percent of NB-patients showed potentially therapeutically targetable somatic alterations (including point mutations, copy number variations and mRNA over-expression). Sixty-six percent of samples showed alterations, graded as "very high priority", that are validated to be directly targetable by an approved drug or an investigational agent. A molecular targeted therapy was applied for four patients, while a genetic counseling was suggested to two patients having one pathogenic germline variant in known cancer predisposition genes. Out of eleven samples implanted in mice, five gave rise to (C)PDX, all preserved in a local PDX Bio-bank. Interestingly, comparing all molecular alterations and histological and immunophenotypic features among the original patient's tumors and PDX/CPDX up to second generation, a high grade of similarity was observed. Notably, also 3D models conserved immunophenotypic features and molecular alterations of the original tumors. CONCLUSIONS: PREME confirms the possibility of identifying targetable genomic alterations in NB, indeed, a molecular targeted therapy was applied to four NB patients. PREME paves the way to the creation of clinically relevant repositories of faithful patient-derived (C)PDX and 3D models, on which testing precision, NB standard-of-care and experimental medicines.

FGFR1 is a potential therapeutic target in neuroblastoma.

BACKGROUND: FGFR1 regulates cell-cell adhesion and extracellular matrix architecture and acts as oncogene in several cancers. Potential cancer driver mutations of FGFR1 occur in neuroblastoma (NB), a neural crest-derived pediatric tumor arising in sympathetic nervous system, but so far they have not been studied experimentally. We investigated the driver-oncogene role of FGFR1 and the implication of N546K mutation in therapy-resistance in NB cells. METHODS: Public datasets were used to predict the correlation of FGFR1 expression with NB clinical outcomes. Whole genome sequencing data of 19 paired diagnostic and relapse NB samples were used to find somatic mutations. In NB cell lines, silencing by short hairpin RNA and transient overexpression of FGFR1 were performed to evaluate the effect of the identified mutation by cell growth, invasion and cologenicity assays. HEK293, SHSY5Y and SKNBE2 were selected to investigate subcellular wild-type and mutated protein localization. FGFR1 inhibitor (AZD4547), alone or in combination with PI3K inhibitor (GDC0941), was used to rescue malignant phenotypes induced by overexpression of FGFR1 wild-type and mutated protein. RESULTS: High FGFR1 expression correlated with low relapse-free survival in two independent NB gene expression datasets. In addition, we found the somatic mutation N546K, the most recurrent point mutation of FGFR1 in all cancers and already reported in NB, in one out of 19 matched primary and recurrent tumors. Loss of FGFR1 function attenuated invasion and cologenicity in NB cells, whereas FGFR1 overexpression enhanced oncogenicity. The overexpression of FGFR1N546K protein showed a higher nuclear localization compared to wild-type protein and increased cellular invasion and cologenicity. Moreover, N546K mutation caused the failure in response to treatment with FGFR1 inhibitor by activation of ERK, STAT3 and AKT pathways. The combination of FGFR1 and PI3K pathway inhibitors was effective in reducing the invasive and colonigenic ability of cells overexpressing FGFR1 mutated protein. CONCLUSIONS: FGFR1 is an actionable driver oncogene in NB and a promising therapy may consist in targeting FGFR1 mutations in patients with therapy-resistant NB.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene.

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.

Association of PARP1 polymorphisms with response to chemotherapy in patients with high-risk neuroblastoma.

The genetic aetiology and the molecular mechanisms that characterize high-risk neuroblastoma are still little understood. The majority of high-risk neuroblastoma patients do not take advantage of current induction therapy. So far, one of the main reasons liable for cancer therapeutic failure is the acquisition of resistance to cytotoxic anticancer drugs, because of the DNA repair system of tumour cells. PARP1 is one of the main DNA damage sensors involved in the DNA repair system and genomic stability. We observed that high PARP1 mRNA level is associated with unfavourable prognosis in 3 public gene expression NB patients' datasets and in 20 neuroblastomas analysed by qRT-PCR. Among 4983 SNPs in PARP1, we selected two potential functional SNPs. We investigated the association of rs907187, in PARP1 promoter, and rs2048426 in non-coding region with response chemotherapy in 121 Italian patients with high-risk NB. Results showed that minor G allele of rs907187 associated with induction response of patients (P = .02) and with decrease PARP1 mRNA levels in NB cell line (P = .003). Furthermore, rs907187 was predicted to alter the binding site of E2F1 transcription factor. Specifically, allele G had low binding affinity with E2F1 whose expression positively correlates with PARP1 expression and associated with poor prognosis of patients with NB. By contrast, we did not find genetic association for the SNP rs2048426. These data reveal rs907187 as a novel potential risk variant associated with the failure of induction therapy for high-risk NB.

HIF-1 transcription activity: HIF1A driven response in normoxia and in hypoxia.

BACKGROUND: HIF1A (Hypoxia-Inducible-Factor 1A) expression in solid tumors is relevant to establish resistance to therapeutic approaches. The use of compounds direct against hypoxia signaling and HIF1A does not show clinical efficiency because of changeable oxygen concentrations in solid tumor areas. The identification of HIF1A targets expressed in both normoxia and hypoxia and of HIF1A/hypoxia signatures might meliorate the prognostic stratification and therapeutic successes in patients with high-risk solid tumors. METHODS: In this study, we conducted a combined analysis of RNA expression and DNA methylation of neuroblastoma cells silenced or unsilenced for HIF1A expression, grown in normoxia and hypoxia conditions. RESULTS: The analysis of pathways highlights HIF-1 (heterodimeric transcription factor 1) activity in normoxia in metabolic process and HIF-1 activity in hypoxia in neuronal differentiation process. HIF1A driven transcriptional response in hypoxia depends on epigenetic control at DNA methylation status of gene regulatory regions. Furthermore, low oxygen levels generate HIF1A-dependent or HIF1A-independent signatures, able to stratify patients according to risk categories. CONCLUSIONS: These findings may help to understand the molecular mechanisms by which low oxygen levels reshape gene signatures and provide new direction for hypoxia targeting in solid tumor.

Fine mapping of 2q35 high-risk neuroblastoma locus reveals independent functional risk variants and suggests full-length BARD1 as tumor-suppressor.

A previous genome-wide association study (GWAS) identified common variation at the BARD1 locus as being highly associated with susceptibility to high-risk neuroblastoma, but the mechanisms underlying this association have been not extensively investigated. Here, we performed a fine mapping analysis of BARD1 locus (2q35) using GWAS data from 556 high-risk neuroblastoma patients and 2,575 controls of European-American ancestry, and identified two independent genome-wide neuroblastoma-associated loci. Functional single-nucleotide polymorphism (SNP) prioritization identified two causative variants that independently contributed to neuroblastoma risk, and each replicated robustly in multiple independent cohorts comprising 445 high-risk cases and 3,170 controls (rs17489363: combined p = 1.07 × 10-31 , OR:1.79, 95% CI:1.62-1.98 and rs1048108: combined p = 7.27 × 10-14 , OR:0.65, 95% CI:0.58-0.73). Particularly, the T risk allele of rs17489363 in the canonical promoter region of full-length BARD1 altered binding site of the transcription factor HSF1 and correlated with low expression of full-length BARD1 mRNA and protein. Low-level expression of full-length BARD1 associated with advanced neuroblastoma. In human neuroblastoma cells, attenuating full-length BARD1 increased proliferation and invasion capacity. In conclusion, we have identified two potentially causative SNPs at the BARD1 locus associated with predisposition to high-risk neuroblastoma, and have shown that full-length BARD1 may act as tumor suppressor.

Kinome expression profiling of human neuroblastoma tumors identifies potential drug targets for ultra high-risk patients.

Neuroblastoma (NBL) accounts for >7% of malignancies in patients younger than 15 years. Low- and intermediate-risk patients exhibit excellent or good prognosis after treatment, whereas for high-risk (HR) patients, the estimated 5-year survival rates is still <40%. The ability to stratify HR patients that will not respond to standard treatment strategies is critical for informed treatment decisions. In this study, we have generated a specific kinome gene signature, named Kinome-27, which is able to identify a subset of HR-NBL tumors, named ultra-HR NBL, with highly aggressive clinical behavior that not adequately respond to standard treatments. We have demonstrated that NBL cell lines expressing the same kinome signature of ultra-HR tumors (ultra-HR-like cell lines) may be selectively targeted by the use of two drugs [suberoylanilide hydroxamic acid (SAHA) and Radicicol], and that the synergic combination of these drugs is able to block the ultra-HR-like cells in G2/M phase of cell cycle. The use of our signature in clinical practice will allow identifying patients with negative outcome, which would benefit from new and more personalized treatments. Preclinical in vivo studies are needed to consolidate the SAHA and Radicicol treatment in ultra-HR NBL patients.

Proteomic Alterations in Response to Hypoxia Inducible Factor 2α in Normoxic Neuroblastoma Cells.

Hypoxia inducible factor (HIF)-2α protein expression in solid tumors promotes stem-like phenotype in cancer stem cells and increases tumorigenic potential in nonstem cancer cells. Recently, we have shown that HIF-1/2α gene expression is correlated to neuroblastoma (NB) poor survival and to undifferentiated tumor state; HIF-2α protein was demonstrated to enhance aggressive features of the disease. In this study, we used proteomic experiments on NB cells to investigate HIF-2α downstream-regulated proteins or pathways with the aim of providing novel therapeutic targets or bad prognosis markers. We verified that pathways mostly altered by HIF-2α perturbation are involved in tumor progression. In particular, HIF-2α induces alteration of central metabolism and splicing control pathways. Simultaneously, WNT, RAS/MAPK, and PI3K/AKT activity or expression are affected and may impact the sensitivity and the intensity of HIF-2α-regulated pathways. Furthermore, genes coding the identified HIF-2α-related markers built a signature able to stratify NB patients with unfavorable outcome. Taken together, our findings underline the relevance of dissecting the downstream effects of a poor survival marker in developing targeted therapy and improving patient stratification. Future prospective studies are needed to translate the use of these data into the clinical practice.

CFDP1 is a neuroblastoma susceptibility gene that regulates transcription factors of the noradrenergic cell identity.

Pleiotropic genetic factors (e.g., DNA polymorphisms) may be involved in the initiation of neuroblastoma (NB) and coronary artery disease (CAD) given their common origin from defects in neural crest development. To discover novel NB susceptibility genes, we conducted a three-stage survey including a meta-analysis of NB and CAD genome-wide association data, prioritization of NB causal variants, and validation in an independent cohort of affected individuals-control subjects. The lead SNP, rs13337397 at the 16q23.1 locus, associated with both diseases in the meta-analysis and with NB in the validation study. All the SNPs in linkage disequilibrium with rs13337397 were annotated using the H3K27ac epigenetic marker of neural crest cells (NCC) and NB cell lines. Indeed, we identified the functional SNP rs13337017, mapping within an enhancer of NCCs and NB cell lines and showing long-range interactions with CFDP1 by Hi-C analysis. Luciferase assays indicated that the risk allele of rs13337017 increased CFDP1 expression in NB cell lines. Of note, CFDP1 high expression associated with unfavorable prognostic markers in an analysis including 498 NB transcriptomes. Moreover, depletion of CFDP1 markedly decreased viability and migration and increased apoptotic rates in NB cell lines. Finally, transcriptome and qPCR analyses revealed that the depletion of CFDP1 may affect noradrenergic neuron differentiation by downregulating master regulators of sympathetic noradrenergic identity, including PHOX2B, HAND2, and GATA3. Our data strongly suggest that CFDP1 acts as oncogene in NB. In addition, we provide evidence that genetic predisposition to NB can be mediated by the alteration of noradrenergic lineage-specific gene expression.

Single-cell transcriptomics of neuroblastoma identifies chemoresistance-associated genes and pathways.

High-Risk neuroblastoma (NB) survival rate is still <50%, despite treatments being more and more aggressive. The biggest hurdle liable to cancer therapy failure is the drug resistance by tumor cells that is likely due to the intra-tumor heterogeneity (ITH). To investigate the link between ITH and therapy resistance in NB, we performed a single cell RNA sequencing (scRNAseq) of etoposide and cisplatin resistant NB and their parental cells. Our analysis showed a clear separation of resistant and parental cells for both conditions by identifying 8 distinct tumor clusters in etoposide-resistant/parental and 7 in cisplatin-resistant/parental cells. We discovered that drug resistance can affect NB cell identities; highlighting the bi-directional ability of adrenergic-to-mesenchymal transition of NB cells. The biological processes driving the identified resistant cell subpopulations revealed genes such as (BARD1, BRCA1, PARP1, HISTH1 axis, members of RPL family), suggesting a potential drug resistance due to the acquisition of DNA repair mechanisms and to the modification of the drug targets. Deconvolution analysis of bulk RNAseq data from 498 tumors with cell subpopulation signatures showed that the transcriptional heterogeneity of our cellular models reflected the ITH of NB tumors and allowed the identification of clusters associated with worse/better survival. Our study demonstrates the distinct cell populations characterized by genes involved in different biological processes can have a role in NB drug treatment failure. These findings evidence the importance of ITH in NB drug resistance studies and the chance that scRNA-seq analysis offers in the identification of genes and pathways liable for drug resistance.

Functional characterization of full-length BARD1 strengthens its role as a tumor suppressor in neuroblastoma.

BARD1 is associated with the development of high-risk neuroblastoma patients. Particularly, the expression of full length (FL) isoform, FL BARD1, correlates to high-risk neuroblastoma development and its inhibition is sufficient to induce neuroblastoma cells towards a worst phenotype. Here we have investigated the mechanisms of FL BARD1 in neuroblastoma cell lines depleted for FL BARD1 expression. We have shown that FL BARD1 expression protects the cells from spontaneous DNA damage and from damage accumulated after irradiation. We demonstrated a role for FL BARD1 as tumor suppressor to prevent unscheduled mitotic entry of DNA damaged cells and to lead to death cells that have bypassed cell cycle checkpoints. FL BARD1-depleted cells that have survived to checkpoints acquire features of aggressiveness. Overall, our results show that FL BARD1 may defend cells against cancer and prevent malignant transformation of cells.

Transcription Factors Involved in Tumorigenesis Are Over-Represented in Mutated Active DNA-Binding Sites in Neuroblastoma.

The contribution of coding mutations to oncogenesis has been largely clarified, whereas little is known about somatic mutations in noncoding DNA and their role in driving tumors remains controversial. Here, we used an alternative approach to interpret the functional significance of noncoding somatic mutations in promoting tumorigenesis. Noncoding somatic mutations of 151 neuroblastomas were integrated with ENCODE data to locate somatic mutations in regulatory elements specifically active in neuroblastoma cells, nonspecifically active in neuroblastoma cells, and nonactive. Within these types of elements, transcription factors (TF) were identified whose binding sites were enriched or depleted in mutations. For these TFs, a gene expression signature was built to assess their implication in neuroblastoma. DNA- and RNA-sequencing data were integrated to assess the effects of those mutations on mRNA levels. The pathogenicity of mutations was significantly higher in transcription factor binding site (TFBS) of regulatory elements specifically active in neuroblastoma cells, as compared with the others. Within these elements, there were 18 over-represented TFs involved mainly in cell-cycle phase transitions and 15 under-represented TFs primarily regulating cell differentiation. A gene expression signature based on over-represented TFs correlated with poor survival and unfavorable prognostic markers. Moreover, recurrent mutations in TFBS of over-represented TFs such as EZH2 affected MCF2L and ADP-ribosylhydrolase like 1 expression, among the others. We propose a novel approach to study the involvement of regulatory variants in neuroblastoma that could be extended to other cancers and provide further evidence that alterations of gene expression may have relevant effects in neuroblastoma development. SIGNIFICANCE: These findings propose a novel approach to study regulatory variants in neuroblastoma and suggest that noncoding somatic mutations have relevant implications in neuroblastoma development.

Exploring Shared Susceptibility between Two Neural Crest Cells Originating Conditions: Neuroblastoma and Congenital Heart Disease.

In the past years, genome wide association studies (GWAS) have provided evidence that inter-individual susceptibility to diverse pathological conditions can reveal a common genetic architecture. Through the analysis of congenital heart disease (CHD) and neuroblastoma (NB) GWAS data, we aimed to dissect the genetic susceptibility shared between these conditions, which are known to arise from neural crest cell (NCC) migration or development abnormalities, via identification and functional characterization of common regions of association. Two loci (2q35 and 3q25.32) harbor single nucleotide polymorphisms (SNPs) that are associated at a p-value < 10-3 with conotruncal malformations and ventricular septal defect respectively, as well as with NB. In addition, the lead SNP in 4p16.2 for atrial septal defect and the lead SNP in 3q25.32 for tetralogy of Fallot are less than 250 Kb distant from the lead SNPs for NB at the same genomic regions. Some of these shared susceptibility loci regulate the expression of relevant genes involved in NCC formation and developmental processes (such as BARD1,MSX1, and SHOX2) and are enriched in several epigenetic markers from NB and fetal heart cell lines. Although the clinical correlation between NB and CHD is unclear, our exploration of a possible common genetic basis between NB and a subset of cardiac malformations can help shed light on their shared embryological origin and pathogenetic mechanisms.

A High-Content Screening of Anticancer Compounds Suggests the Multiple Tyrosine Kinase Inhibitor Ponatinib for Repurposing in Neuroblastoma Therapy.

Novel druggable targets have been discovered in neuroblastoma (NB), paving the way for more effective treatments. However, children with high-risk NB still show high mortality rates prompting for a search of novel therapeutic options. Here, we aimed at repurposing FDA-approved drugs for NB treatment by performing a high-content screening of a 349 anticancer compounds library. In the primary screening, we employed three NB cell lines, grown as three-dimensional (3D) multicellular spheroids, which were treated with 10 µmol/L of the library compounds for 72 hours. The viability of 3D spheroids was evaluated using a high-content imaging approach, resulting in a primary hit list of 193 compounds. We selected 60 FDA-approved molecules and prioritized drugs with multi-target activity, discarding those already in use for NB treatment or enrolled in NB clinical trials. Hence, 20 drugs were further tested for their efficacy in inhibiting NB cell viability, both in two-dimensional and 3D models. Dose-response curves were then supplemented with the data on side effects, therapeutic index, and molecular targets, suggesting two multiple tyrosine kinase inhibitors, ponatinib and axitinib, as promising candidates for repositioning in NB. Indeed, both drugs showed induction of cell-cycle block and apoptosis, as well as inhibition of colony formation. However, only ponatinib consistently affected migration and inhibited invasion of NB cells. Finally, ponatinib also proved effective inhibition of tumor growth in orthotopic NB mice, providing the rationale for its repurposing in NB therapy. Mol Cancer Ther; 17(7); 1405-15. ©2018 AACR.

Inhibition of hypoxia inducible factors combined with all-trans retinoic acid treatment enhances glial transdifferentiation of neuroblastoma cells.

Neuroblastoma (NBL) is a heterogeneous tumor characterized by a wide range of clinical manifestations. A high tumor cell differentiation grade correlates to a favorable stage and positive outcome. Expression of the hypoxia inducible factors HIF1-α (HIF1A gene) and HIF2-α (EPAS1 gene) and/or hypoxia-regulated pathways has been shown to promote the undifferentiated phenotype of NBL cells. Our hypothesis is that HIF1A and EPAS1 expression represent one of the mechanisms responsible for the lack of responsiveness of NBL to differentiation therapy. Clinically, high levels of HIF1A and EPAS1 expression were associated with inferior survival in two NBL microarray datasets, and patient subgroups with lower expression of HIF1A and EPAS1 showed significant enrichment of pathways related to neuronal differentiation. In NBL cell lines, the combination of all-trans retinoic acid (ATRA) with HIF1A or EPAS1 silencing led to an acquired glial-cell phenotype and enhanced expression of glial-cell differentiation markers. Furthermore, HIF1A or EPAS1 silencing might promote cell senescence independent of ATRA treatment. Taken together, our data suggest that HIF inhibition coupled with ATRA treatment promotes differentiation into a more benign phenotype and cell senescence in vitro. These findings open the way for additional lines of attack in the treatment of NBL minimal residue disease.