A biology-driven approach identifies the hypoxia gene signature as a predictor of the outcome of neuroblastoma patients

Mutually exclusive teams-like patterns of gene regulation characterize phenotypic heterogeneity along the noradrenergic-mesenchymal axis in neuroblastoma

Neuroblastoma is the most frequent extracranial pediatric tumor and leads to 15% of all cancer-related deaths in children. Tumor relapse and therapy resistance in neuroblastoma are driven by phenotypic plasticity and heterogeneity between noradrenergic (NOR) and mesenchymal (MES) cell states. Despite the importance of this phenotypic plasticity, the dynamics and molecular patterns associated with these bidirectional cell-state transitions remain relatively poorly understood. Here, we analyze multiple RNA-seq datasets at both bulk and single-cell resolution, to understand the association between NOR- and MES-specific factors. We observed that NOR-specific and MES-specific expression patterns are largely mutually exclusive, exhibiting a "teams-like" behavior among the genes involved, reminiscent of our earlier observations in lung cancer and melanoma. This antagonism between NOR and MES phenotypes was also associated with metabolic reprogramming and with immunotherapy targets PD-L1 and GD2 as well as with experimental perturbations driving the NOR-MES and/or MES-NOR transition. Further, these "teams-like" patterns were seen only among the NOR- and MES-specific genes, but not in housekeeping genes, possibly highlighting a hallmark of network topology enabling cancer cell plasticity.

Targeting long noncoding RNAs in neuroblastoma: Progress and prospects

Neuroblastoma (NB) is the third most prevalent tumor that mostly influences infants and young children. Although different treatments have been developed for the treatment of NB, high-risk patients have been reported to have low survival rates. Currently, long noncoding RNAs (lncRNAs) have shown an attractive potential in cancer research and a party of investigations have been performed to understand mechanisms underlying tumor development through lncRNA dysregulation. Researchers have just newly initiated to exhibit the involvement of lncRNAs in NB pathogenesis. In this review article, we tried to clarify the point we stand with respect to the involvement of lncRNAs in NB. Moreover, implications for the pathologic roles of lncRNAs in the development of NB have been discussed. It seems that some of these lncRNAs have promising potential to be applied as biomarkers for NB prognosis and treatment.

Development and validation of a 21-gene prognostic signature in neuroblastoma

Survival outcomes for patients with neuroblastoma vary markedly and reliable prognostic markers and risk stratification tools are lacking. We sought to identify and validate a transcriptomic signature capable of predicting risk of mortality in patients with neuroblastoma. The TARGET NBL dataset (n = 243) was used to develop the model and two independent cohorts, E-MTAB-179 (n = 478) and GSE85047 (n = 240) were used as validation sets. EFS was the primary outcome and OS was the secondary outcome of interest for all analysis. We identified a 21-gene signature capable of stratifying neuroblastoma patients into high and low risk groups in the E-MTAB-179 (HR 5.87 [3.83-9.01], p < 0.0001, 5 year AUC 0.827) and GSE85047 (HR 3.74 [2.36-5.92], p < 0.0001, 5 year AUC 0.815) validation cohorts. Moreover, the signature remained independent of known clinicopathological variables, and remained prognostic within clinically important subgroups. Further, the signature was effectively incorporated into a risk model with clinicopathological variables to improve prognostic performance across validation cohorts (Pooled Validation HR 6.93 [4.89-9.83], p < 0.0001, 5 year AUC 0.839). Similar prognostic utility was also demonstrated with OS. The identified signature is a robust independent predictor of EFS and OS outcomes in neuroblastoma patients and can be combined with clinically utilized clinicopathological variables to improve prognostic performance.

Comparative epigenomics by machine learning approach for neuroblastoma

BACKGROUND

Neuroblastoma (NB) is the second most common pediatric solid tumor. Because the number of genetic mutations found in tumors are small, even in some patients with unfavorable NB, epigenetic variation is expected to play an important role in NB progression. DNA methylation is a major epigenetic mechanism, and its relationship with NB prognosis has been a concern. One limitation with the analysis of variation in DNA methylation is the lack of a suitable analytical model. Therefore, in this study, we performed a random forest (RF) analysis of the DNA methylome data of NB from multiple databases.

RESULTS

RF is a popular machine learning model owing to its simplicity, intuitiveness, and computational cost. RF analysis identified novel intermediate-risk patient groups with characteristic DNA methylation patterns within the low-risk group. Feature selection analysis based on probe annotation revealed that enhancer-annotated regions had strong predictive power, particularly for MYCN-amplified NBs. We developed a gene-based analytical model to identify candidate genes related to disease progression, such as PRDM8 and FAM13A-AS1. RF analysis revealed sufficient predictive power compared to other machine learning models.

CONCLUSIONS

RF is a useful tool for DNA methylome analysis in cancer epigenetic studies, and has potential to identify a novel cancer-related genes.

ViLoN-a multi-layer network approach to data integration demonstrated for patient stratification

With more and more data being collected, modern network representations exploit the complementary nature of different data sources as well as similarities across patients. We here introduce the Variation of information fused Layers of Networks algorithm (ViLoN), a novel network-based approach for the integration of multiple molecular profiles. As a key innovation, it directly incorporates prior functional knowledge (KEGG, GO). In the constructed network of patients, patients are represented by networks of pathways, comprising genes that are linked by common functions and joint regulation in the disease. Patient stratification remains a key challenge both in the clinic and for research on disease mechanisms and treatments. We thus validated ViLoN for patient stratification on multiple data type combinations (gene expression, methylation, copy number), showing substantial improvements and consistently competitive performance for all. Notably, the incorporation of prior functional knowledge was critical for good results in the smaller cohorts (rectum adenocarcinoma: 90, esophageal carcinoma: 180), where alternative methods failed.

The Role of Imaging Biomarkers to Guide Pharmacological Interventions Targeting Tumor Hypoxia

Hypoxia is a common feature of solid tumors that contributes to angiogenesis, invasiveness, metastasis, altered metabolism and genomic instability. As hypoxia is a major actor in tumor progression and resistance to radiotherapy, chemotherapy and immunotherapy, multiple approaches have emerged to target tumor hypoxia. It includes among others pharmacological interventions designed to alleviate tumor hypoxia at the time of radiation therapy, prodrugs that are selectively activated in hypoxic cells or inhibitors of molecular targets involved in hypoxic cell survival (i.e., hypoxia inducible factors HIFs, PI3K/AKT/mTOR pathway, unfolded protein response). While numerous strategies were successful in pre-clinical models, their translation in the clinical practice has been disappointing so far. This therapeutic failure often results from the absence of appropriate stratification of patients that could benefit from targeted interventions. Companion diagnostics may help at different levels of the research and development, and in matching a patient to a specific intervention targeting hypoxia. In this review, we discuss the relative merits of the existing hypoxia biomarkers, their current status and the challenges for their future validation as companion diagnostics adapted to the nature of the intervention.

Identification of a Glycosyltransferase Signature for Predicting Prognosis and Immune Microenvironment in Neuroblastoma

Neuroblastoma (NB) is one of the most common solid tumors in children. Glycosyltransferases (GTs) play a crucial role in tumor development and immune escape and have been used as prognostic biomarkers in various tumors. However, the biological functions and prognostic significance of GTs in NB remain poorly understood. The expression data from Gene Expression Omnibus (GEO) and Therapeutically Applicable Research to Generate Effective Treatments (TARGET) were collected as training and testing data. Based on a progression status, differentially expressed GTs were identified. We constructed a GTscore through support vector machine, least absolute shrinkage and selection operator, and Cox regression in NB, which included four prognostic GTs and was an independent prognostic risk factor for NB. Patients in the high GTscore group had an older age, MYCN amplification, advanced International Neuroblastoma Staging System stage, and high risk. Samples with high GTscores revealed high disialoganglioside (GD2) and neuron-specific enolase expression levels. In addition, a lack of immune cell infiltration was observed in the high GTscore group. This GTscore was also associated with the expression of chemokines (CCL2, CXCL9, and CXCL10) and immune checkpoint genes (cytotoxic T-lymphocyte–associated protein 4, granzyme H, and granzyme K). A low GTscore was also linked to an enhanced response to anti–PD-1 immunotherapy in melanoma patients, and one type of tumor was also derived from neuroectodermal cells such as NB. In conclusion, the constructed GTscore revealed the relationship between GT expression and the NB outcome, GD2 phenotype, and immune infiltration and provided novel clues for the prediction of prognosis and immunotherapy response in NB.

Revisiting the HIF switch in the tumor and its immune microenvironment

Hypoxia is a hallmark of all solid tumors and their metastases. This leads to activation of the hypoxia-inducible factor (HIF) family of transcription factors, which modulate gene expression within both tumor cells and immune cells within the tumor microenvironment, influencing tumor progression and treatment response. The best characterized HIF isoforms, HIF-1α and HIF-2α, show nonoverlapping and often antagonistic roles. With the recent availability of inhibitors that target one or both HIFs, including the first-in-class selective HIF-2α inhibitor belzutifan, the prospect of HIF-α isoform-selective targeting is now a reality. Here, we summarize current knowledge on the unique contributions of the two HIF-α isoforms to tumor progression in the context of the complex tumor immune microenvironment, highlighting important considerations for therapy.

Connectivity Map Analysis Indicates PI3K/Akt/mTOR Inhibitors as Potential Anti-Hypoxia Drugs in Neuroblastoma

Neuroblastoma (NB) is one of the deadliest pediatric cancers, accounting for 15% of deaths in childhood. Hypoxia is a condition of low oxygen tension occurring in solid tumors and has an unfavorable prognostic factor for NB. In the present study, we aimed to identify novel promising drugs for NB treatment. Connectivity Map (CMap), an online resource for drug repurposing, was used to identify connections between hypoxia-modulated genes in NB tumors and compounds. Two sets of 34 and 21 genes up- and down-regulated between hypoxic and normoxic primary NB tumors, respectively, were analyzed with CMap. The analysis reported a significant negative connectivity score across nine cell lines for 19 compounds mainly belonging to the class of PI3K/Akt/mTOR inhibitors. The gene expression profiles of NB cells cultured under hypoxic conditions and treated with the mTORC complex inhibitor PP242, referred to as the Mohlin dataset, was used to validate the CMap findings. A heat map representation of hypoxia-modulated genes in the Mohlin dataset and the gene set enrichment analysis (GSEA) showed an opposite regulation of these genes in the set of NB cells treated with the mTORC inhibitor PP242. In conclusion, our analysis identified inhibitors of the PI3K/Akt/mTOR signaling pathway as novel candidate compounds to treat NB patients with hypoxic tumors and a poor prognosis.

Dysregulation of Zinc Finger Protein 395 Contributes to the Pathogenesis of Chondrosarcoma

Introduction

The transcription factor zinc finger protein 395 (ZNF395) is involved in several cellular responses and tumorigenesis. However, the potential role and clinical significance of ZNF395 in chondrosarcoma are not well investigated. This study determines the expression profile, prognostic value and biological function of ZNF395 in human chondrosarcoma.

Methods

The mRNA and protein expressions of ZNF395 in fresh chondrosarcomas and the matched adjacent non-tumor tissues were assessed using real-time PCR and immunoblotting, respectively. The protein expression of ZNF395 in chondrosarcoma specimens was evaluated by immunohistochemistry, and the relationships among its protein level, clinicopathological parameters and prognosis were further detected. Cell viability, colony formation, migration, invasion and apoptosis assay were evaluated in chondrosarcoma cells with depletion of ZNF395.

Results

The mRNA and protein expressions of ZNF395 in chondrosarcoma tissues were significantly higher than those in the matched adjacent non-tumor tissues and benign cartilage tumors. Clinical analysis displayed that ZNF395 was expressed at higher levels in chondrosarcoma patients with higher histological grade and advanced MSTS stage. Furthermore, we demonstrated that high expression of ZNF395 correlated with a worse overall survival of chondrosarcoma patients. Multivariate Cox regression analysis indicated that ZNF395 was an independent prognostic marker in chondrosarcoma patients. Functional studies revealed that depletion of ZNF395 markedly inhibited cell viability, colony formation, migration and invasion, and promoted apoptosis in chondrosarcoma.

Conclusion

These findings suggest that dysregulation of ZNF395 contributes to chondrosarcoma development, and ZNF395 may act as a potent oncogene and serve as a independently prognostic factor, highlight the potential of ZNF395 as a novel biomarker and therapeutic target for chondrosarcoma.

Molecular Genetics in Neuroblastoma Prognosis

In recent years, much research has been carried out to identify the biological and genetic characteristics of the neuroblastoma (NB) tumor in order to precisely define the prognostic subgroups for improving treatment stratification. This review will describe the major genetic features and the recent scientific advances, focusing on their impact on diagnosis, prognosis, and therapeutic solutions in NB clinical management.

Selective Nonmethylated CpG DNA Recognition Mechanism of Cysteine Clamp Domains

Methylation of DNA at CpG sites is a major mark for epigenetic regulation, but how transcription factors are influenced by CpG methylation is not well understood. Here, we report the molecular mechanisms of how the TCF (T-cell factor) and GEF (glucose transporter 4 enhancer factor) families of proteins selectively target unmethylated DNA sequences with a C-clamp type zinc finger domain. The structure of the C-clamp domain from human GEF family protein HDBP1 (C-clamp^HDBP1) in complex with DNA was determined using NMR spectroscopy, which adopts a unique zinc finger fold and selectively binds RCCGG (R = A/G) DNA sequences with an "Arg···Trp-Lys-Lys" DNA recognition motif inserted in the major groove. The CpG base pairs are central to the binding due to multiple hydrogen bonds formed with the backbone carbonyl groups of Trp378 and Lys379, as well as the side chain ε-amino groups of Lys379 and Lys380 from C-clamp^HDBP1. Consequently, methylation of the CpG dinucleotide almost abolishes the binding. Homology modeling reveals that the C-clamp domain from human TCF1E (C-clamp^TCF1E) binds DNA through essentially the same mechanism, with a similar "Arg···Arg-Lys-Lys" DNA recognition motif. The substitution of tryptophan by arginine makes C-clamp^HDBP1 prefer RCCGC DNA sequences. The two signature DNA recognition motifs are invariant in the GEF and TCF families of proteins, respectively, from fly to human. The recognition of the CpG dinucleotide through two consecutive backbone carbonyl groups is the same as that of the CXXC type unmethylated CpG DNA binding domains, suggesting a common mechanism shared by unmethylated CpG binding proteins.

Lost in application: Measuring hypoxia for radiotherapy optimisation

The history of radiotherapy is intertwined with research on hypoxia. There is level 1a evidence that giving hypoxia-targeting treatments with radiotherapy improves locoregional control and survival without compromising late side-effects. Despite coming in and out of vogue over decades, there is now an established role for hypoxia in driving molecular alterations promoting tumour progression and metastases. While tumour genomic complexity and immune profiling offer promise, there is a stronger evidence base for personalising radiotherapy based on hypoxia status. Despite this, there is only one phase III trial targeting hypoxia modification with full transcriptomic data available. There are no biomarkers in routine use for patients undergoing radiotherapy to aid management decisions, and a roadmap is needed to ensure consistency and provide a benchmark for progression to application. Gene expression signatures address past limitations of hypoxia biomarkers and could progress biologically optimised radiotherapy. Here, we review recent developments in generating hypoxia gene expression signatures and highlight progress addressing the challenges that must be overcome to pave the way for their clinical application.

Development of a stacked, porous silk scaffold neuroblastoma model for investigating spatial differences in cell and drug responsiveness

Development of in vitro, preclinical cancer models that contain cell-driven microenvironments remains a challenge. Engineering of millimeter-scale, in vitro tumor models with spatially distinct regions that can be independently assessed to study tumor microenvironments has been limited. Here, we report the use of porous silk scaffolds to generate a high cell density neuroblastoma (NB) model that can spatially recapitulate changes resulting from cell and diffusion driven changes. Using COMSOL modeling, a scaffold holder design that facilitates stacking of thin, 200 μm silk scaffolds into a thick, bulk millimeter-scale tumor model (2, 4, 6, and 8 stacked scaffolds) and supports cell-driven oxygen gradients was developed. Cell-driven oxygen gradients were confirmed through pimonidazole staining. Post-culture, the stacked scaffolds were separated for analysis on a layer-by-layer basis. The analysis of each scaffold layer demonstrated decreasing DNA and increasing expression of hypoxia related genes (VEGF, CAIX, and GLUT1) from the exterior scaffolds to the interior scaffolds. Furthermore, the expression of hypoxia related genes at the interior of the stacks was comparable to that of a single scaffold cultured under 1% O₂ and at the exterior of the stacks was comparable to that of a single scaffold cultured under 21% O₂. The four-stack scaffold model underwent further evaluation to determine if a hypoxia activated drug, tirapazamine, induced reduced cell viability within the internal stacks (region of reduced oxygen) as compared with the external stacks. Decreased DNA content was observed in the internal stacks as compared to the external stacks when treated with tirapazamine, which suggests the internal scaffold stacks had higher levels of hypoxia than the external scaffolds. This stacked silk scaffold system presents a method for creating a single culture model capable of generating controllable cell-driven microenvironments through different stacks that can be individually assessed and used for drug screening.

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.

Prognostic significance of MYCN related genes in pediatric neuroblastoma: a study based on TARGET and GEO datasets

Background

Neuroblastoma patients with MYCN amplification are associated with poor prognosis. However, the prognostic relevance of MYCN associated genes in neuroblastoma is unclear.

Methods

The expression profiles of MYCN associated genes were identified from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and Gene Expression Omnibus (GEO) datasets. Enriched transcription factors and signaling pathways were determined using gene set enrichment analysis (GSEA). Kaplan-Meier plotter was used to identify the prognostic relevance of MYCN associated genes. Multivariate cox regression and Spearman’s correlation were used to determine the correlation coefficients of MYCN associated genes.

Results

In TARGET and GSE85047 datasets, neuroblastoma patients with MYCN amplification were associated with worse prognosis. Transcription factor MYC was positively associated with MYCN amplification in GSEA assay. We identified 13 MYC target genes which were increased in neuroblastoma patients with MYCN amplification in TARGET, GSE19274 and GSE85047 datasets. Moreover, six out of the 13 MYC target genes ARMC6, DCTPP1, EIF4G1, ELOVL6, FBL and PRMT1 were associated with adverse prognosis in TARGET and GSE85047 datasets. Transcription factor E2F1 was up-regulated by MYCN amplification and associated with the poor prognosis of neuroblastoma. Furthermore, RPS19 in ribosome signaling pathway was also associated with MYCN amplification and correlated with the poor prognosis of neuroblastoma. At last, we showed that most of MYCN target genes were correlated with each other. However, EIF4G1 was an independent prognostic marker. And the prognostic effects of the combination of MYCN amplification and EIF4G1 expression were more significant than MYCN or EIF4G1 alone.

Conclusions

MYCN target genes ARMC6, DCTPP1, EIF4G1, ELOVL6, FBL, PRMT1, E2F1 and RPS19 had significant prognostic effects in pediatric neuroblastoma. And neuroblastoma patients without MYCN amplification and low EIF4G1 expression had best prognosis.

Modelling of hypoxia gene expression for three different cancer cell lines

Gene dynamic analysis is essential in identifying target genes involved pathogenesis of various diseases, including cancer. Cancer prognosis is often influenced by hypoxia. We apply a multi-step pipeline to study dynamic gene expressions in response to hypoxia in three cancer cell lines: prostate (DU145), colon (HT29), and breast (MCF7) cancers. We identified 26 distinct temporal expression patterns for prostate cell line, and 29 patterns for colon and breast cell lines. The module-based dynamic networks have been developed for all three cell lines. Our analyses improve the existing results in multiple ways. It exploits the time-dependence nature of gene expression values in identifying the dynamically significant genes; hence, more key significant genes and transcription factors have been identified. Our gene network returns significant information regarding biologically important modules of genes. Furthermore, the network has potential in learning the regulatory path between transcription factors and the downstream genes. In addition, our findings suggest that changes in genes BMP6 and ARSJ expression might have a key role in the time-dependent response to hypoxia in breast cancer.

ARNT-dependent HIF-2 transcriptional activity is not sufficient to regulate downstream target genes in neuroblastoma

BACKGROUND

Hypoxia-inducible factor (HIF)-2α associates with poor outcome in neuroblastoma and glioblastoma, and gain-of-function mutations in the EPAS1 gene (encoding HIF-2α) have been reported in paragangliomas and pheochromocytomas. Specific targeting of a druggable hydrophobic pocket in the HIF-2α PAS-B domain with PT2385 have demonstrated promising clinical results for clear cell renal cell carcinoma (ccRCC). Here, we investigated the effect of PT2385-mediated inhibition of ARNT dependent HIF-2 activity.

METHODS

Neuroblastoma patient-derived xenograft (PDX) cells were treated with PT2385 and analyzed for HIF-2-dependent gene expression, HIF activity, HIF-2α protein localization, response to chemotherapy and orthotopic tumor growth in vivo. Two-sided student t-test was used.

RESULTS

We detected high levels of HIF-2α protein in perivascular niches in neuroblastoma PDXs in vivo and at oxygenated conditions in PDX-derived cell cultures in vitro, particularly in the cytoplasmic fraction. Nuclear HIF-2α expression was reduced following PT2385 treatment, but surprisingly, virtually no effects on tumor growth in vivo or expression of canonical HIF downstream target genes in vitro were observed. In coherence, RNA sequencing of PT2385-treated PDX cells revealed a virtually unaffected transcriptome. Treatment with PT2385 did not affect cellular response to chemotherapy. In contrast, HIF-2α protein knockdown resulted in profound downregulation of target genes.

CONCLUSIONS

The lack of effect from PT2385 treatment in combination with high cytoplasmic HIF-2α expression at normoxia suggest that HIF-2α have additional roles than acting as an ARNT dependent transcription factor. It is important to further unravel the conditions at which HIF-2α has transcriptional and non-transcriptional roles in neuroblastoma.

Transcriptome analysis defines myocardium gene signatures in children with ToF and ASD and reveals disease-specific molecular reprogramming in response to surgery with cardiopulmonary bypass

BACKGROUND

Tetralogy of Fallot (ToF) and Atrial Septal Defects (ASD) are the most common types of congenital heart diseases and a major cause of childhood morbidity and mortality. Cardiopulmonary bypass (CPB) is used during corrective cardiac surgery to support circulation and heart stabilization. However, this procedure triggers systemic inflammatory and stress response and consequent increased risk of postoperative complications. The aim of this study was to define the molecular bases of ToF and ASD pathogenesis and response to CPB and identify new potential biomarkers.

METHODS

Comparative transcriptome analysis of right atrium specimens collected from 10 ToF and 10 ASD patients was conducted before (Pre-CPB) and after (Post-CPB) corrective surgery. Total RNA isolated from each sample was individually hybridized on Affymetrix HG-U133 Plus Array Strips containing 38,500 unique human genes. Differences in the gene expression profiles and functional enrichment/network analyses were assessed using bioinformatic tools. qRT-PCR analysis was used to validate gene modulation.

RESULTS

Pre-CPB samples showed significant differential expression of a total of 72 genes, 28 of which were overexpressed in ToF and 44 in ASD. According to Gene Ontology annotation, the mostly enriched biological processes were represented by matrix organization and cell adhesion in ToF and by muscle development and contractility in ASD specimens. GSEA highlighted the specific enrichment of hypoxia gene sets in ToF samples, pointing to a role for hypoxia in disease pathogenesis. The post-CPB myocardium exhibited significant alterations in the expression profile of genes related to transcription regulation, growth/apoptosis, inflammation, adhesion/matrix organization, and oxidative stress. Among them, only 70 were common to the two disease groups, whereas 110 and 24 were unique in ToF and ASD, respectively. Multiple functional interactions among differentially expressed gene products were predicted by network analysis. Interestingly, gene expression changes in ASD samples followed a consensus hypoxia profile.

CONCLUSION

Our results provide a comprehensive view of gene reprogramming in right atrium tissues of ToF and ASD patients before and after CPB, defining specific molecular pathways underlying disease pathophysiology and myocardium response to CPB. These findings have potential translational value because they identify new candidate prognostic markers and targets for tailored cardioprotective post-surgical therapies.

Integrating tumor hypoxic stress in novel and more adaptable strategies for cancer immunotherapy

Immunotherapy is poised to become an increasingly utilized therapy in the treatment of cancer. However, several abnormalities in the tumor microenvironment (TME) that can thwart the efficacy of immunotherapies have been established. Microenvironmental hypoxia is a determining factor in shaping aggressiveness, metastatic potential and treatment resistance of solid tumors. The characterization of this phenomenon could prove beneficial for determining a patient's treatment path and for the introduction of novel targetable factors that can enhance therapeutic outcome. Indeed, the ablation of hypoxia has the potential to sensitize tumors to immunotherapy by metabolically remodeling their microenvironment. In this review, we discuss the intrinsic contributions of hypoxia to cellular plasticity, heterogeneity, stemness and genetic instability in the context of immune escape. In addition, we will shed light on how managing hypoxia can ameliorate response to immunotherapy and how integrating hypoxia gene signatures could play a role in this pursuit.

Exosomal microRNAs from Longitudinal Liquid Biopsies for the Prediction of Response to Induction Chemotherapy in High-Risk Neuroblastoma Patients: A Proof of Concept SIOPEN Study

Despite intensive treatment, 50% of children with high-risk neuroblastoma (HR-NB) succumb to their disease. Progression through current trials evaluating the efficacy of new treatments for children with HR disease usually depends on an inadequate response to induction chemotherapy, assessed using imaging modalities. In this study, we sought to identify circulating biomarkers that might be detected in a simple blood sample to predict patient response to induction chemotherapy. Since exosomes released by tumor cells can drive tumor growth and chemoresistance, we tested the hypothesis that exosomal microRNA (exo-miRNAs) in blood might predict response to induction chemotherapy. The exo-miRNAs expression profile in plasma samples collected from children treated in HR-NBL-1/SIOPEN before and after induction chemotherapy was compared to identify a three exo-miRs signature that could discriminate between poor and good responders. Exo-miRNAs expression also provided a chemoresistance index predicting the good or poor prognosis of HR-NB patients.

Decrease of Nibrin expression in chronic hypoxia is associated with hypoxia-induced chemoresistance in some brain tumour cells

Background

Solid tumours are less oxygenated than normal tissues. This is called tumour hypoxia and leads to resistance to radiotherapy and chemotherapy. The molecular mechanisms underlying such resistance have been investigated in a range of tumour types, including the adult brain tumours glioblastoma, yet little is known for paediatric brain tumours. Medulloblastoma (MB) is the most common malignant brain tumour in children. We aimed to elucidate the impact of hypoxia on the sensitivity of MB cells to chemo- and radiotherapy.

Methods

We used two MB cell line (D283-MED and MEB-Med8A) and a widely used glioblastoma cell line (U87MG) for comparison. We applied a range of molecular and cellular techniques to measure cell survival, cell cycle progression, protein expression and DNA damage combined with a transcriptomic micro-array approach in D283-MED cells, for global gene expression analysis in acute and chronic hypoxic conditions.

Results

In D283-MED and U87MG, chronic hypoxia (5 days), but not acute hypoxia (24 h) induced resistance to chemotherapy and X-ray irradiation. This acquired resistance upon chronic hypoxia was present but less pronounced in MEB-Med8A cells. Using transcriptomic analysis in D283-MED cells, we found a large transcriptional remodelling upon long term hypoxia, in particular the expression of a number of genes involved in detection and repair of double strand breaks (DSB) was altered. The levels of Nibrin (NBN) and MRE11, members of the MRN complex (MRE11/Rad50/NBN) responsible for DSB recognition, were significantly down-regulated. This was associated with a reduction of Ataxia Telangiectasia Mutated (ATM) activation by etoposide, indicating a profound dampening of the DNA damage signalling in hypoxic conditions. As a consequence, p53 activation by etoposide was reduced, and cell survival enhanced. Whilst U87MG shared the same dampened p53 activity, upon chemotherapeutic drug treatment in chronic hypoxic conditions, these cells used a different mechanism, independent of the DNA damage pathway.

Conclusion

Together our results demonstrate a new mechanism explaining hypoxia-induced resistance involving the alteration of the response to DSB in D283-MED cells, but also highlight the cell type to cell type diversity and the necessity to take into account the differing tumour genetic make-up when considering re-sensitisation therapeutic protocols.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5476-9) contains supplementary material, which is available to authorized users.

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.

Modification of Biodistribution and Brain Uptake of Copper Bis(thiosemicarbazonato) Complexes by the Incorporation of Amine and Polyamine Functional Groups

The synthesis of new bis(thiosemicarbazonato)copper(II) complexes featuring polyamine substituents via selective transamination reactions is presented. Polyamines of different lengths, with different ionizable substituent groups, were used to modify and adjust the hydrophilic/lipophilic balance of the copper complexes. The new analogues were radiolabeled with copper-64 and their lipophilicities estimated using distribution coefficients. The cell uptake of the new polyamine complexes was investigated with preliminary in vitro biological studies using a neuroblastoma cancer cell line. The in vivo biodistribution of three of the new analogues was investigated in vivo in mice using positron-emission tomography imaging, and one of the new complexes was compared to [⁶⁴Cu]Cu(atsm) in an A431 squamous cell carcinoma xenograft model. Modification of the copper complexes with various amine-containing functional groups alters the biodistribution of the complexes in mice. One complex, with a pendent ( N, N-dimethylamino)ethane functional group, displayed tumor uptake similar to that of [⁶⁴Cu]Cu(atsm) but higher brain uptake, suggesting that this compound has the potential to be of use in the diagnostic brain imaging of tumors and neurodegenerative diseases.

Nucleolin Is a Functional Binding Protein for Salinomycin in Neuroblastoma Stem Cells

The aim of this study is to illuminate a novel therapeutic approach by identifying a functional binding target of salinomycin, an emerging anticancer stem cell (CSC) agent, and to help dissect the underlying action mechanisms. By utilizing integrated strategies, we identify that nucleolin (NCL) is likely a salinomycin-binding target and a critical regulator involved in human neuroblastoma (NB) CSC activity. Salinomycin markedly suppresses NB CD34 expression and reduces CD34+ cell population in an NCL-dependent manner via disruption of the interaction of NCL with CD34 promoter. The elevated levels of NCL expression in NB tumors are associated with poor patient survival. Altogether, these results indicate that NCL is likely a novel functional salinomycin-binding target that exhibits the potential to be a prognostic marker for NB therapy.

Identification of a Prognostic Hypoxia-Associated Gene Set in IDH-Mutant Glioma

Glioma growth is often accompanied by a hypoxic microenvironment favorable for the induction and maintenance of the glioma stem cell (GSC) phenotype. Due to the paucity of cell models of Isocitrate Dehydrogenase 1 mutant (IDH1^mut) GSCs, biology under hypoxic conditions has not been sufficiently studied as compared to IDH1 wildtype (IDH1^wt) GSCs. We therefore grew well-characterized IDH1^mut (n = 4) and IDH1^wt (n = 4) GSC lines under normoxic (20%) and hypoxic (1.5%) culture conditions and harvested mRNA after 72 h. Transcriptome analyses were performed and hypoxia regulated genes were further analyzed using the expression and clinical data of the lower grade glioma cohort of The Cancer Genome Atlas (LGG TCGA) in a confirmatory approach and to test for possible survival associations. Results show that global expression changes were more pronounced in IDH1^wt than in IDH1^mut GSCs. However, when focusing on known hypoxia-regulated gene sets, enrichment analyses showed a comparable regulation in both IDH1^mut and IDH1^wt GSCs. Of 272 significantly up-regulated genes under hypoxic conditions in IDH1^mut GSCs a hypoxia-related survival score (HRS-score) of five genes (LYVE1, FAM162A, WNT6, OTP, PLOD1) was identified by the Least Absolute Shrinkage and Selection Operator (LASSO) algorithm which was able to predict survival independent of age, 1p19q co-deletion status and WHO grade (II vs. III) in the LGG TCGA cohort and in the Rembrandt dataset. Altogether, we were able to identify and validate a novel hypoxia-related survival score in IDH1^mut GSCs consisting of five hypoxia-regulated genes which was significantly associated with patient survival independent of known prognostic confounders.

Effects of hypoxic preconditioning on neuroblastoma tumour oxygenation and metabolic signature in a chick embryo model

Hypoxia episodes and areas in tumours have been associated with metastatic dissemination and poor prognosis. Given the link between tumour tissue oxygen levels and cellular metabolic activity, we hypothesised that the metabolic profile between metastatic and non-metastatic tumours would reveal potential new biomarkers and signalling cues. We have used a previously established chick embryo model for neuroblastoma growth and metastasis, where the metastatic phenotype can be controlled by neuroblastoma cell hypoxic preconditioning (3 days at 1% O₂). We measured, with fibre-optic oxygen sensors, the effects of the hypoxic preconditioning on the tumour oxygenation, within tumours formed by SK-N-AS cells on the chorioallantoic membrane (CAM) of chick embryos. We found that the difference between the metastatic and non-metastatic intratumoural oxygen levels was small (0.35% O₂), with a mean below 1.5% O₂ for most tumours. The metabolomic profiling, using NMR spectroscopy, of neuroblastoma cells cultured in normoxia or hypoxia for 3 days, and of the tumours formed by these cells showed that the effects of hypoxia in vitro did not compare with in vivo tumours. One notable difference was the high levels of the glycolytic end-products triggered by hypoxia in vitro, but not by hypoxia preconditioning in tumours, likely due to the very high basal levels of these metabolites in tumours compared with cells. In conclusion, we have identified high levels of ketones (3-hydroxybutyrate), lactate and phosphocholine in hypoxic preconditioned tumours, all known to fuel tumour growth, and we herein point to the poor relevance of in vitro metabolomic experiments for cancer research.

65 YEARS OF THE DOUBLE HELIX: Classification of endocrine tumors in the age of integrated genomics

The classification of human cancers represents one of the cornerstones of modern pathology. Over the last century, surgical pathologists established the current taxonomy of neoplasia using traditional histopathological parameters, which include tumor architecture, cytological features and cellular proliferation. This morphological classification is efficient and robust with high reproducibility and has served patients and health care providers well. The most recent decade has witnessed an explosion of genome-wide molecular genetic and epigenetic data for most cancers, including tumors of endocrine organs. The availability of this expansive multi-dimensional genomic data, collectively termed the cancer genome, has catalyzed a re-examination of the classification of endocrine tumors. Here, recent cancer genome studies of various endocrine tumors, including those of the thyroid, pituitary and adrenal glands, pancreas, small bowel, lung and skin, are presented with special emphasis on how genomic insights are impacting endocrine tumor classification.

High-risk, Expression-Based Prognostic Long Noncoding RNA Signature in Neuroblastoma

Background

Current clinical risk factors stratify patients with neuroblastoma (NB) for appropriate treatments, yet patients with similar clinical behaviors evoke variable responses. MYCN amplification is one of the established drivers of NB and, when combined with high-risk displays, worsens outcomes. Growing high-throughput transcriptomics studies suggest long noncoding RNA (lncRNA) dysregulation in cancers, including NB. However, expression-based lncRNA signatures are altered by MYCN amplification, which is associated with high-risk, and patient prognosis remains limited.

Methods

We investigated RNA-seq-based expression profiles of lncRNAs in MYCN status and risk status in a discovery cohort (n = 493) and validated them in three independent cohorts. In the discovery cohort, a prognostic association of lncRNAs was determined by univariate Cox regression and integrated into a signature using the risk score method. A novel risk score threshold selection criterion was developed to stratify patients into risk groups. Outcomes by risk group and clinical subgroup were assessed using Kaplan-Meier survival curves and multivariable Cox regression. The performance of lncRNA signatures was evaluated by receiver operating characteristic curve. All statistical tests were two-sided.

Results

In the discovery cohort, 16 lncRNAs that were differentially expressed (fold change ≥ 2 and adjusted P ≤ 0.01) integrated into a prognostic signature. A high risk score group of lncRNA signature had poor event-free survival (EFS; P < 1E-16). Notably, lncRNA signature was independent of other clinical risk factors when predicting EFS (hazard ratio = 3.21, P = 5.95E-07). The findings were confirmed in independent cohorts (P = 2.86E-02, P = 6.18E-03, P = 9.39E-03, respectively). Finally, the lncRNA signature had higher accuracy for EFS prediction (area under the curve = 0.788, 95% confidence interval = 0.746 to 0.831).

Conclusions

Here, we report the first (to our knowledge) RNA-seq 16-lncRNA prognostic signature for NB that may contribute to precise clinical stratification and EFS prediction.

CHL1 gene acts as a tumor suppressor in human neuroblastoma

Neuroblastoma is an aggressive, relapse-prone childhood tumor of the sympathetic nervous system that accounts for 15% of pediatric cancer deaths. A distal portion of human chromosome 3p is often deleted in neuroblastoma, this region may contain one or more putative tumor suppressor genes. A 2.54 Mb region at 3p26.3 encompassing the smallest region of deletion pinpointed CHL1 gene, the locus for neuronal cell adhesion molecule close homolog of L1. We found that low CHL1 expression predicted poor outcome in neuroblastoma patients. Here we have used two inducible cell models to analyze the impact of CHL1 on neuroblastoma biology. Over-expression of CHL1 induced neurite-like outgrowth and markers of neuronal differentiation in neuroblastoma cells, halted tumor progression, inhibited anchorage-independent colony formation, and suppressed the growth of human tumor xenografts. Conversely, knock-down of CHL1 induced neurite retraction and activation of Rho GTPases, enhanced cell proliferation and migration, triggered colony formation and anchorage-independent growth, accelerated growth in orthotopic xenografts mouse model. Our findings demonstrate unambiguously that CHL1 acts as a regulator of proliferation and differentiation of neuroblastoma cells through inhibition of the MAPKs and Akt pathways. CHL1 is a novel candidate tumor suppressor in neuroblastoma, and its associated pathways may represent a promising target for future therapeutic interventions.

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.

ZEB1-repressed microRNAs inhibit autocrine signaling that promotes vascular mimicry of breast cancer cells

During normal tumor growth and in response to some therapies, tumor cells experience acute or chronic deprivation of nutrients and oxygen and induce tumor vascularization. While this occurs predominately through sprouting angiogenesis, tumor cells have also been shown to directly contribute to vessel formation through vascular mimicry (VM) and/or endothelial transdifferentiation. The extrinsic and intrinsic mechanisms underlying tumor cell adoption of endothelial phenotypes, however, are not well understood. Here we show that serum withdrawal induces mesenchymal breast cancer cells to undergo VM and that knockdown of the epithelial-to-mesenchymal transition (EMT) regulator, Zinc finger E-box binding homeobox 1 (ZEB1), or overexpression of the ZEB1-repressed microRNAs (miRNAs), miR-200c, miR-183, miR-96 and miR-182 inhibits this process. We find that secreted proteins Fibronectin 1 (FN1) and serine protease inhibitor (serpin) family E member 2 (SERPINE2) are essential for VM in this system. These secreted factors are upregulated in mesenchymal cells in response to serum withdrawal, and overexpression of VM-inhibiting miRNAs abrogates this upregulation. Intriguingly, the receptors for these secreted proteins, low-density lipoprotein receptor-related protein 1 (LRP1) and Integrin beta 1 (ITGB1), are also targets of the VM-inhibiting miRNAs, suggesting that autocrine signaling stimulating VM is regulated by ZEB1-repressed miRNA clusters. Together, these data provide mechanistic insight into the regulation of VM and suggest that miRNAs repressed during EMT, in addition to suppressing migratory and stem-like properties of tumor cells, also inhibit endothelial phenotypes of breast cancer cells adopted in response to a nutrient-deficient microenvironment.

Role of hypoxia in Diffuse Large B-cell Lymphoma: Metabolic repression and selective translation of HK2 facilitates development of DLBCL

Published molecular profiling studies in patients with lymphoma suggested the influence of hypoxia inducible factor-1 alpha (HIF1α) targets in prognosis of DLBCL. Yet, the role of hypoxia in hematological malignancies remains unclear. We observed that activation of HIF1α resulted in global translation repression during hypoxic stress in DLBCL. Protein translation efficiency as measured using 35S-labeled methionine incorporation revealed a ≥50% reduction in translation upon activation of HIF1α. Importantly, translation was not completely inhibited and expression of clinically correlated hypoxia targets such as GLUT1, HK2, and CYT-C was found to be refractory to translational repression under hypoxia in DLBCL cells. Notably, hypoxic induction of these genes was not observed in normal primary B-cells. Translational repression was coupled with a decrease in mitochondrial function. Screening of primary DLBCL patient samples revealed that expression of HK2, which encodes for the enzyme hexokinase 2, was significantly correlated with DLBCL phenotype. Genetic knockdown studies demonstrated that HK2 is required for promoting growth of DLBCL under hypoxic stress. Altogether, our findings provide strong support for the direct contribution of HK2 in B-cell lymphoma development and suggest that HK2 is a key metabolic driver of the DLBCL phenotype.

Perturbation-response genes reveal signaling footprints in cancer gene expression

Aberrant cell signaling can cause cancer and other diseases and is a focal point of drug research. A common approach is to infer signaling activity of pathways from gene expression. However, mapping gene expression to pathway components disregards the effect of post-translational modifications, and downstream signatures represent very specific experimental conditions. Here we present PROGENy, a method that overcomes both limitations by leveraging a large compendium of publicly available perturbation experiments to yield a common core of Pathway RespOnsive GENes. Unlike pathway mapping methods, PROGENy can (i) recover the effect of known driver mutations, (ii) provide or improve strong markers for drug indications, and (iii) distinguish between oncogenic and tumor suppressor pathways for patient survival. Collectively, these results show that PROGENy accurately infers pathway activity from gene expression in a wide range of conditions.

c-MYC inhibition impairs hypoxia response in glioblastoma multiforme

The c-MYC oncoprotein is a DNA binding transcription factor that enhances the expression of many active genes. c-MYC transcriptional signatures vary according to the transcriptional program defined in each cell type during differentiation. Little is known on the involvement of c-MYC in regulation of gene expression programs that are induced by extracellular cues such as a changing microenvironment. Here we demonstrate that inhibition of c-MYC in glioblastoma multiforme cells blunts hypoxia-dependent glycolytic reprogramming and mitochondria fragmentation in hypoxia. This happens because c-MYC inhibition alters the cell transcriptional response to hypoxia and finely tunes the expression of a subset of Hypoxia Inducible Factor 1-regulated genes. We also show that genes whose expression in hypoxia is affected by c-MYC inhibition are able to distinguish the Proneural subtype of glioblastoma multiforme, thus potentially providing a molecular signature for this class of tumors that are the least tractable among glioblastomas.

Coexpression network analysis identifies transcriptional modules associated with genomic alterations in neuroblastoma

Neuroblastoma is a highly complex and heterogeneous cancer in children. Acquired genomic alterations including MYCN amplification, 1p deletion and 11q deletion are important risk factors and biomarkers in neuroblastoma. Here, we performed a co-expression-based gene network analysis to study the intrinsic association between specific genomic changes and transcriptome organization. We identified multiple gene coexpression modules which are recurrent in two independent datasets and associated with functional pathways including nervous system development, cell cycle, immune system process and extracellular matrix/space. Our results also indicated that modules involved in nervous system development and cell cycle are highly associated with MYCN amplification and 1p deletion, while modules responding to immune system process are associated with MYCN amplification only. In summary, this integrated analysis provides novel insights into molecular heterogeneity and pathogenesis of neuroblastoma. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

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.

Early and late effects of pharmacological ALK inhibition on the neuroblastoma transcriptome

Background

Neuroblastoma is an aggressive childhood malignancy of the sympathetic nervous system. Despite multi-modal therapy, survival of high-risk patients remains disappointingly low, underscoring the need for novel treatment strategies. The discovery of ALK activating mutations opened the way to precision treatment in a subset of these patients. Previously, we investigated the transcriptional effects of pharmacological ALK inhibition on neuroblastoma cell lines, six hours after TAE684 administration, resulting in the 77-gene ALK signature, which was shown to gradually decrease from 120 minutes after TAE684 treatment, to gain deeper insight into the molecular effects of oncogenic ALK signaling.

Aim

Here, we further dissected the transcriptional dynamic profiles of neuroblastoma cells upon TAE684 treatment in a detailed timeframe of ten minutes up to six hours after inhibition, in order to identify additional early targets for combination treatment.

Results

We observed an unexpected initial upregulation of positively regulated MYCN target genes following subsequent downregulation of overall MYCN activity. In addition, we identified adrenomedullin (ADM), previously shown to be implicated in sunitinib resistance, as the earliest response gene upon ALK inhibition.

Conclusions

We describe the early and late effects of ALK inhibitor TAE684 treatment on the neuroblastoma transcriptome. The observed unexpected upregulation of ADM warrants further investigation in relation to putative ALK resistance in neuroblastoma patients currently undergoing ALK inhibitor treatment.

A Gene Signature for Selecting Benefit from Hypoxia Modification of Radiotherapy for High-Risk Bladder Cancer Patients

Purpose: Hypoxia modification improves overall survival in muscle-invasive bladder cancer patients who undergo radiotherapy. There is evidence that hypoxic tumors benefit most from hypoxia modification. The study aimed to identify or derive a hypoxia gene signature that predicts benefit from hypoxia-modifying treatment in bladder cancer.Experimental Design: Published hypoxia signatures were tested and a new one derived by analyzing bladder cancer transcriptomic data from public databases. Tumor samples were available from the BCON phase III randomized trial of radiotherapy alone or with carbogen and nicotinamide (CON). Gene expression data were generated for 151 tumors using Affymetrix Human 1.0 Exon ST arrays and used for independent validation.Results: A 24-gene signature was derived, which was prognostic in four of six independent surgical cohorts (n = 679; meta HR, 2.32; 95% CI, 1.73-3.12; P < 0.0001). The signature was also prognostic in BCON patients receiving radiotherapy alone (n = 75; HR for local relapse-free survival, 2.37; 95% CI, 1.26-4.47; P = 0.0076). The signature predicted benefit from CON (n = 76; HR, 0.47; 95% CI, 0.26-0.86; P = 0.015). Prognostic significance (P = 0.017) and predictive significance (P = 0.058) remained after adjusting for clinicopathologic variables. A test for interaction between hypoxia status and treatment arms was significant (P = 0.0094).Conclusions: A 24-gene hypoxia signature has strong and independent prognostic and predictive value for muscle-invasive bladder cancer patients. The signature can aid identification of patients likely to benefit from the addition of carbogen and nicotinamide to radiotherapy. Clin Cancer Res; 23(16); 4761-8. ©2017 AACR.

Combined epigenetic and differentiation-based treatment inhibits neuroblastoma tumor growth and links HIF2α to tumor suppression

Neuroblastoma is a pediatric cancer characterized by variable outcomes ranging from spontaneous regression to life-threatening progression. High-risk neuroblastoma patients receive myeloablative chemotherapy with hematopoietic stem-cell transplant followed by adjuvant retinoid differentiation treatment. However, the overall survival remains low; hence, there is an urgent need for alternative therapeutic approaches. One feature of high-risk neuroblastoma is the high level of DNA methylation of putative tumor suppressors. Combining the reversibility of DNA methylation with the differentiation-promoting activity of retinoic acid (RA) could provide an alternative strategy to treat high-risk neuroblastoma. Here we show that treatment with the DNA-demethylating drug 5-Aza-deoxycytidine (AZA) restores high-risk neuroblastoma sensitivity to RA. Combined systemic distribution of AZA and RA impedes tumor growth and prolongs survival. Genome-wide analysis of treated tumors reveals that this combined treatment rapidly induces a HIF2α-associated hypoxia-like transcriptional response followed by an increase in neuronal gene expression and a decrease in cell-cycle gene expression. A small-molecule inhibitor of HIF2α activity diminishes the tumor response to AZA+RA treatment, indicating that the increase in HIF2α levels is a key component in tumor response to AZA+RA. The link between increased HIF2α levels and inhibited tumor growth is reflected in large neuroblastoma patient datasets. Therein, high levels of HIF2α, but not HIF1α, significantly correlate with expression of neuronal differentiation genes and better prognosis but negatively correlate with key features of high-risk tumors, such as MYCN amplification. Thus, contrary to previous studies, our findings indicate an unanticipated tumor-suppressive role for HIF2α in neuroblastoma.

ZNF395 Is an Activator of a Subset of IFN-Stimulated Genes

Activation of the interferon (IFN) pathway in response to infection with pathogens results in the induction of IFN-stimulated genes (ISGs) including proinflammatory cytokines, which mount the proper antiviral immune response. However, aberrant expression of these genes is pathogenic to the host. In addition to IFN-induced transcription factors non-IFN-regulated factors contribute to the transcriptional control of ISGs. Here, we show by genome wide expression analysis, siRNA-mediated suppression and Doxycycline-induced overexpression that the cellular transcription factor ZNF395 activates a subset of ISGs including the chemokines CXCL10 and CXCL11 in keratinocytes. We found that ZNF395 acts independently of IFN but enhances the IFN-induced expression of CXCL10 and CXCL11. Luciferase reporter assays revealed a requirement of intact NFκB-binding sites for ZNF395 to stimulate the CXCL10 promoter. The transcriptional activation of CXCL10 and CXCL11 by ZNF395 was abolished after inhibition of IKK by BMS-345541, which increased the stability of ZNF395. ZNF395 encodes at least two motifs that mediate the enhanced degradation of ZNF395 in response to IKK activation. Thus, IKK is required for ZNF395-mediated activation of transcription and enhances its turn-over to keep the activity of ZNF395 low. Our results support a previously unrecognized role of ZNF395 in the innate immune response and inflammation.

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.

Neuropeptide Y as a Biomarker and Therapeutic Target for Neuroblastoma

Neuroblastoma (NB) is a pediatric malignant neoplasm of sympathoadrenal origin. Challenges in its management include stratification of this heterogeneous disease and a lack of both adequate treatments for high-risk patients and noninvasive biomarkers of disease progression. Our previous studies have identified neuropeptide Y (NPY), a sympathetic neurotransmitter expressed in NB, as a potential therapeutic target for these tumors by virtue of its Y5 receptor (Y5R)-mediated chemoresistance and Y2 receptor (Y2R)-mediated proliferative and angiogenic activities. The goal of this study was to determine the clinical relevance and utility of these findings. Expression of NPY and its receptors was evaluated in corresponding samples of tumor RNA, tissues, and sera from 87 patients with neuroblastic tumors and in tumor tissues from the TH-MYCN NB mouse model. Elevated serum NPY levels correlated with an adverse clinical presentation, poor survival, metastasis, and relapse, whereas strong Y5R immunoreactivity was a marker of angioinvasive tumor cells. In NB tissues from TH-MYCN mice, high immunoreactivity of both NPY and Y5R marked angioinvasive NB cells. Y2R was uniformly expressed in undifferentiated tumor cells, which supports its previously reported role in NB cell proliferation. Our findings validate NPY as a therapeutic target for advanced NB and implicate the NPY/Y5R axis in disease dissemination. The correlation between elevated systemic NPY and NB progression identifies serum NPY as a novel NB biomarker.

Polysialic acid sustains cancer cell survival and migratory capacity in a hypoxic environment

Polysialic acid (polySia) is a unique carbohydrate polymer expressed on the surface of NCAM (neuronal cell adhesion molecule) in a number of cancers where it modulates cell-cell and cell-matrix adhesion, migration, invasion and metastasis and is strongly associated with poor clinical prognosis. We have carried out the first investigation into the effect of polySia expression on the behaviour of cancer cells in hypoxia, a key source of chemoresistance in tumours. The role of polysialylation and associated tumour cell migration and cell adhesion were studied in hypoxia, along with effects on cell survival and the potential role of HIF-1. Our findings provide the first evidence that polySia expression sustains migratory capacity and is associated with tumour cell survival in hypoxia. Initial mechanistic studies indicate a potential role for HIF-1 in sustaining polySia-mediated migratory capacity, but not cell survival. These data add to the growing body of evidence pointing to a crucial role for the polysialyltransferases (polySTs) in neuroendocrine tumour progression and provide the first evidence to suggest that polySia is associated with an aggressive phenotype in tumour hypoxia. These results have significant potential implications for polyST inhibition as an anti-metastatic therapeutic strategy and for targeting hypoxic cancer cells.

C-myc overexpression drives melanoma metastasis by promoting vasculogenic mimicry via c-myc/snail/Bax signaling

c-Myc is a well-characterized proto-oncogene that induces cellular transformation and modulates programmed cell death. While recent studies have demonstrated high expression of c-Myc protein in advanced and metastatic melanoma, the clinical and biological implications remain to be fully elucidated. In this study, we investigated the effect of c-Myc overexpression in melanoma tumorigenesis. Clinicopathological analysis demonstrated that c-Myc expression positively correlated with the formation of vasculogenic mimicry (VM) and linearly patterned programmed cell necrosis (LPPCN). Clinically, high c-Myc expression was significantly associated with distant metastasis and poor prognosis, while biologically, c-Myc overexpression led to significant increases in cell motility, invasiveness and metastasis. Moreover, c-Myc induced the formation of VM and promoted the expression of epithelial-mesenchymal transition (EMT)-associated protein Snail both in vivo and in vitro. High expression of c-Myc increased Bax expression in hypoxic conditions and induced cell apoptosis. Taken together, we conclude that c-Myc overexpression promotes the formation of VM by EMT and LPPCN in melanoma. Our improved understanding of the clinical and biological effects of c-Myc overexpression in melanoma highlights the incomplete understanding of this oncogene, and indicates that c-Myc is a potential therapeutic target of this disease.

KEY MESSAGE

High c-Myc expression is associated with tumor metastasis and poor prognosis in human melanoma. c-Myc upregulates Snail expression to promote EMT via the TGF-β/Snail/Ecadherin signal pathway. c-Myc leads to cell death by upregulating Bax expression causing a lower Bcl2/Bax ratio under severe hypoxic conditions. c-Myc promotes vasculogenic mimicry and linearly patterned programmed cell necrosis.

Transcript signatures that predict outcome and identify targetable pathways in MYCN-amplified neuroblastoma

BACKGROUND

In the pediatric cancer neuroblastoma (NB), patients are stratified into low, intermediate or high-risk subsets based in part on MYCN amplification status. While MYCN amplification in general predicts unfavorable outcome, no clinical or genomic factors have been identified that predict outcome within these cohorts of high-risk patients. In particular, it is currently not possible at diagnosis to determine which high-risk neuroblastoma patients will ultimately fail upfront therapy.

EXPERIMENTAL DESIGN

We analyzed the prognostic potential of most published gene expression signatures for NB and developed a new prognostic signature to predict outcome for patients with MYCN amplification. Network and pathway analyses identified candidate therapeutic targets for this MYCN-amplified patient subset with poor outcome.

RESULTS

Most signatures have a high capacity to predict outcome of unselected NB patients. However, the majority of published signatures, as well as most randomly generated signatures, are highly confounded by MYCN amplification, and fail to predict outcome in subpopulations of high-risk patients with MYCN-amplified NB. We identify a MYCN module signature that predicts patient outcome for those with MYCN-amplified tumors, that also predicts potential tractable therapeutic signaling pathways and targets including the DNA repair enzyme Poly [ADP-ribose] polymerase 1 (PARP1).

CONCLUSION

Many prognostic signatures for NB are confounded by MYCN amplification and fail to predict outcome for the subset of high-risk patients with MYCN amplification. We report a MYCN module signature that is associated with distinct patient outcomes, and predicts candidate therapeutic targets in DNA repair pathways, including PARP1 in MYCN-amplified NB.

Enhancing Interpretability of Gene Signatures with Prior Biological Knowledge

Biological interpretability is a key requirement for the output of microarray data analysis pipelines. The most used pipeline first identifies a gene signature from the acquired measurements and then uses gene enrichment analysis as a tool for functionally characterizing the obtained results. Recently Knowledge Driven Variable Selection (KDVS), an alternative approach which performs both steps at the same time, has been proposed. In this paper, we assess the effectiveness of KDVS against standard approaches on a Parkinson’s Disease (PD) dataset. The presented quantitative analysis is made possible by the construction of a reference list of genes and gene groups associated to PD. Our work shows that KDVS is much more effective than the standard approach in enhancing the interpretability of the obtained results.

An 18 gene expression-based score classifier predicts the clinical outcome in stage 4 neuroblastoma

Background

The prognosis of children with metastatic stage 4 neuroblastoma (NB) has remained poor in the past decade.

Patients and methods

Using microarray analyses of 342 primary tumors, we here developed and validated an easy to use gene expression-based risk score including 18 genes, which can robustly predict the outcome of stage 4 patients.

Results

This classifier was a significant predictor of overall survival in two independent validation cohorts [cohort 1 (n = 214): P = 6.3 × 10⁻⁵; cohort 2 (n = 27): P = 3.1 × 10⁻²]. The prognostic value of the risk score was validated by multivariate analysis including the established markers age and MYCN status (P = 0.027). In the pooled validation cohorts (n = 241), integration of the risk score with the age and/or MYCN status identified subgroups with significantly differing overall survival (ranging from 35 to 100 %).

Conclusion

Together, the 18-gene risk score classifier can identify patients with stage 4 NB with favorable outcome and may therefore improve risk assessment and treatment stratification of NB patients with disseminated disease.

Electronic supplementary material

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