Next Article in Journal
Polylevolysine and Fibronectin-Loaded Nano-Hydroxyapatite/PGLA/Dextran-Based Scaffolds for Improving Bone Regeneration: A Histomorphometric in Animal Study
Previous Article in Journal
Diverse Virulence Attributes of Pantoea alfalfae sp. nov. CQ10 Responsible for Bacterial Leaf Blight in Alfalfa Revealed by Genomic Analysis
Previous Article in Special Issue
Italian Precision Medicine in Pediatric Oncology: Moving beyond Actionable Alterations
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 9
10.3390/ijms24098141
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MYCN Amplification Is Associated with Reduced Expression of Genes Encoding γ-Secretase Complex and NOTCH Signaling Components in Neuroblastoma

by
Prasoon Agarwal
1,
Aleksandra Glowacka
2,
Loay Mahmoud
3,
Wesam Bazzar
2,4,
Lars-Gunnar Larsson
2,4 and
Mohammad Alzrigat
2,4,*
1
National Bioinformatics Infrastructure Sweden (NBIS), Science for Life Laboratory, Division of Occupational and Environmental Medicine, Department of Laboratory Medicine, Lund University, 22362 Lund, Sweden
2
Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 17165 Solna, Sweden
3
Department of Cell and Molecular Biology, Karolinska Institutet, 17177 Stockholm, Sweden
4
Department of Pharmaceutical Biosciences, Biomedical Center, Uppsala University, 75124 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(9), 8141; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24098141
Submission received: 25 January 2023 / Revised: 6 April 2023 / Accepted: 26 April 2023 / Published: 2 May 2023
(This article belongs to the Special Issue Molecular Determinants of Neuroblastoma 2.0)
Browse Figures
Review Reports Versions Notes

Abstract

:
Amplification of the MYCN oncogene is found in ~20% of neuroblastoma (NB) cases and correlates with high-risk disease and poor prognosis. Despite the plethora of studies describing the role of MYCN in NB, the exact molecular mechanisms underlying MYCN’s contribution to high-risk disease are not completely understood. Herein, we implemented an integrative approach combining publicly available RNA-Seq and MYCN ChIP-Seq datasets derived from human NB cell lines to define biological processes directly regulated by MYCN in NB. Our approach revealed that MYCN-amplified NB cell lines, when compared to non-MYCN-amplified cell lines, are characterized by reduced expression of genes involved in NOTCH receptor processing, axoneme assembly, and membrane protein proteolysis. More specifically, we found genes encoding members of the γ-secretase complex, which is known for its ability to liberate several intracellular signaling molecules from membrane-bound proteins such as NOTCH receptors, to be down-regulated in MYCN-amplified NB cell lines. Analysis of MYCN ChIP-Seq data revealed an enrichment of MYCN binding at the transcription start sites of genes encoding γ-secretase complex subunits. Notably, using publicly available gene expression data from NB primary tumors, we revealed that the expression of γ-secretase subunits encoding genes and other components of the NOTCH signaling pathway was also reduced in MYCN-amplified tumors and correlated with worse overall survival in NB patients. Genetic or pharmacological depletion of MYCN in NB cell lines induced the expression of γ-secretase genes and NOTCH-target genes. Chemical inhibition of γ-secretase activity dampened the expression of NOTCH-target genes upon MYCN depletion in NB cells. In conclusion, this study defines a set of MYCN-regulated pathways that are specific to MYCN-amplified NB tumors, and it suggests a novel role for MYCN in the suppression of genes of the γ-secretase complex, with an impact on the NOTCH-target gene expression in MYCN-amplified NB.

1. Introduction

Neuroblastoma (NB) is the most common extra-cranial solid tumor form during childhood, and it is believed to arise from the neural crest [1,2,3]. It has been estimated that NB accounts for approximately 15% of cancer-related deaths in children, with the majority of cases (90%) diagnosed by the age of five years [1,2,3]. NB patients are stratified into low-, intermediate-, and high-risk groups based on age, tumor histology, clinical stage, and genetic makeup [4,5,6,7]. Low- and intermediate-risk patients have favorable outcomes with an 80–95% event-free survival rate, while high-risk patients demonstrate a <50% event-free survival rate [1,2,3]. The current treatment of high-risk patients includes traditional intensive chemotherapy followed by surgical removal, radiation, myeloablation and autologous bone marrow transplantation, and more recently, immunotherapy [8,9,10]. However, the majority of high-risk patients relapse and eventually die due to refractory disease. Thus, high-risk NB is considered a disease with unmet medical need.
Amplification of the MYCN gene is detected in ~20% of all NB cases, and in about 40% of high-risk NB cases, and it is the genetic aberration most consistently associated with high-risk disease and poor survival [11,12]. MYCN belongs to the MYC family of oncoproteins that also includes MYC (c-MYC) and MYCL. The MYC genes encode transcription factors of the basic region/helix-loop-helix/leucine zipper (bHLHZip) family and control genes involved in multiple fundamental cellular processes, such as cell cycle progression, metabolism, apoptosis, senescence, differentiation, stem cell functions, and angiogenesis [13,14,15,16,17,18]. To regulate expression of their target genes, MYC family proteins need to dimerize with another bHLHZip protein called MAX, resulting in a stable heterodimer formation (hereafter called MYC:MAX) that binds and regulates MYC target genes [19,20]. MYC has been described as a general amplifier of transcription for all active genes in a cell [21,22,23,24], while later studies suggest that this effect is indirect, and is observed subsequently as a consequence of MYC’s regulation of specific genes [25,26,27]. In addition to target gene activation, MYC has also been reported to repress the transcription of genes, often in complexes with the zinc finger protein MIZ-1 [26,28,29]. The crucial role of MYC-family oncoproteins in cancer development makes them attractive targets for therapy in cancer, including NB. Therefore, several efforts have been put forward to identify molecules that dampen MYC activity in cancer [30,31,32,33,34]. Even though some of these molecules have demonstrated promising anti-tumor activity in vitro and in vivo, no specific MYC inhibitor is yet used in clinical practice. In addition to MYCN amplification, other genetic aberrations, such as TERT rearrangements or alternative lengthening of telomeres (ALT), also define a subset of high-risk NB with a negative poor outcome [35,36,37]. Furthermore, also other genetic aberrations, such as mutations in ALK, PTPN11, ATRX, TERT, and NRAS, have been shown to correlate with advanced-stage and high-risk NB [38,39].
Although MYCN amplification is correlated with high-risk and aggressive NB, the molecular mechanisms underpinning this association are still somewhat elusive. Herein, we analyzed publicly available gene expression datasets derived from 39 commonly used human NB cell lines representing MYCN-amplified and non-MYCN-amplified NB in an effort to identify pathways that are differentially regulated in MYCN-amplified NB tumors. On the one hand, we uncovered that MYCN-amplified cell lines exhibit increased expression of genes involved in ribosome biogenesis, RNA metabolism, gene expression, and protein synthesis processes. On the other hand, the MYCN-amplified cell lines were characterized by the down-regulated expression of genes related to NOTCH receptors and amyloid protein processing, axoneme assembly, protein proteolysis, localization, and transport. These findings were supported by an in-silico analysis of publicly available gene expression data derived from primary NB tumors, demonstrating the increased expression of genes involved in ribosomal biogenesis and the low expression of genes involved in NOTCH receptor processing and signaling in MYCN-amplified tumors when compared to non-MYCN-amplified tumors. Combining gene expression data with MYCN ChIP-Seq data in NB cell lines suggested the direct role of MYCN in regulating the expression of genes involved in the differentially regulated processes. Importantly, we identified MYCN as a direct or indirect repressor of members of the γ-secretase complex, which is known to activate signaling cascades downstream of membrane-bound proteins, such as NOTCH receptors. Consequently, down-regulation of MYCN through genetic or pharmacological means resulted in up-regulation of the repressed γ-secretase genes. In keeping with this, we demonstrated the increased expression of NOTCH-target genes upon MYCN depletion, which was dampened by the chemical inhibition of the γ-secretase complex in NB cell lines. These data provide new insights into the role of MYCN in NB, identify the genes and biological processes regulated by MYCN, and further suggest MYCN as a negative modulator of genes encoding the γ-secretase complex and NOTCH-target gene expression. The role of the NOTCH signaling pathway in NB is an active field of research, and this study provides further insights into the molecular mechanisms underlying NOTCH signaling regulation in NB.

2. Results

2.1. Identification of MYCN-Regulated Gene Expression Profiles and Biological Processes in MYCN-Amplified Neuroblastoma

To identify the biological processes that are specific to MYCN-amplified NB when compared to non-MYCN-amplified NB, we analyzed gene expression data derived from 39 commonly used human NB cell lines, as described by Harneza et al. [40], and we identified 2738 differentially regulated genes (FDR ≤ 0.05) in MYCN-amplified vs. non-MYCN-amplified NB (Table S1). The gene ontology analysis of biological processes using the online Enricher analysis tool and the Molecular Signatures Database (MSigDB) uncovered that the up-regulated genes were related to ribosome biogenesis, ribosomal RNA synthesis and processing, RNA metabolism, protein translation, and gene expression (Figure 1A and Table S1). The down-regulated genes were enriched in processes related to NOTCH receptor processing, protein proteolysis and localization, peptidase activity, and axoneme assembly (Figure 1B and Table S1). To investigate whether similar correlations exist in primary tumors, we analyzed the gene expression of the SEQC study comprising RNA expression data from a 498-patient cohort with NB (r2.amc.nl/; Tumor Neuroblastoma-SEQC-498-custom-ag44kcwolf). As shown in Figure S1, the expression of genes linked to ribosomal RNA biogenesis (Figure S1A) and the NOTCH receptor processing and signaling pathway (Figure S1B) were, in general, higher and lower, respectively, in MYCN-amplified (green label) tumors when compared to non-MYCN-amplified (red label) tumors.
Next, we set out to uncover the potential direct regulation of the differentially regulated genes by MYCN and analyzed MYCN ChIP-Seq data in three MYCN-amplified human NB cell lines; KELLY, LAN-5, and COGN145 [41]. For the enrichment of the MYCN ChIP-Seq peaks, we defined ±5 kb genomic regions from known transcriptional start sites (TSSs) as a cut-off. Based on this definition, we identified 5870 annotated genes common to the three MYCN-amplified cell lines (Table S2). Overlaying these common genes with the differentially expressed genes in the MYCN-amplified NB cell lines resulted in 475 up-regulated and 497 down-regulated genes, which we define as MYCN-target genes (Figure 1C). Figure 1D,E shows the average pileup of MYCN ChIP-Seq peaks around the TSSs of up- and down-regulated genes in the KELLY cell line and indicates that the MYCN binding sites are enriched around ±1 kb of the TSSs (Figure 1D,E). To explore the transcriptional regulation of genes differentially regulated by MYCN further, we utilized the MSigDB to define the transcription factor binding motifs. This analysis revealed an enrichment of the MYC-type of E-boxes and binding motifs of other transcription factors, such as E2F, NRF1, YY1, and NFY, at the TSSs of the MYCN-bound and up-regulated genes (Figure 1F top panel), while the TSSs of the MYCN-bound and down-regulated genes were mainly enriched with binding motifs for ETS family members, such as ELK1, ETS1, ETS2, GABPB, and TEL (Figure 1F, lower panel). We also found the enrichment of the NRF2, PAX4, and MYC-type E-box binding motifs at the TSSs of the MYCN-bound and down-regulated genes (Figure 1F, lower panel).

2.2. Integrative Analysis of Gene Expression and MYCN ChIP-Seq Reveals That MYCN Directly Associates with the Suppressed ADAM17 and γ-Secretase Complex Genes

Having defined the MYCN-bound and differentially regulated genes, we sought to uncover which of these genes are potentially directly regulated by MYCN. The gene ontology analysis revealed that the 475 up-regulated and MYCN-bound genes were related to RNA synthesis and processing, ribosome biogenesis, RNA metabolism, and protein translation (Figure 2A and Table S2), while the 497 down-regulated MYCN-target genes were involved in NOTCH receptor processing, protein proteolysis and localization, endopeptidase activity, and amyloid protein metabolism (Figure 2A and Table S2). A visualization of the MYCN ChIP-Seq peaks shows that they coincide with the TSSs of the up- and down-regulated genes identified in this study, such as RPL5 and APH1A, respectively (Figure 2B, Figures S2 and S3). Intriguingly, we found the genes encoding for ADAM17 and the γ-secretase complex, which constitute the GO terms related to NOTCH receptor processing, amyloid protein processing, and membrane protein proteolysis, to be MYCN-bound and down-regulated in the MYCN-amplified cell lines (Table 1, Figure S4 and Table S2). Figure 2C,D shows the MYCN ChIP-Seq peak pileup around the TSSs of up-regulated genes involved in ribosomal biogenesis (Figure 2C) and the down-regulated ADAM17 and γ-secretase complex genes (Figure 2D) in the MYCN-amplified KELLY cell line [41], indicating that the MYCN binding sites are enriched around ±1 kb of TSSs.
Considering that MYC family proteins may affect the total RNA content per cell in certain cell types [21,22,23,24], and the fact that the publicly available RNA-Seq data from the study by Harneza et al. analyzed above did not contain spike-in control [40], we analyzed the expression levels of ADAM17, genes encoding the γ-secretase complex, as well as other components of NOTCH signaling in the publicly available gene expression dataset with spike-in RNA control [22]. This analysis revealed that ADAM17 and all γ-secretase encoding genes, with the exception of APH1B, were less expressed in the MYCN-amplified neuroblastoma cell lines KELLY, SK-N-BE(2)C, and NGP compared with the non-MYCN-amplified SK-N-AS cell line (Figure S5). Moreover, the NOTCH-target genes HES1, HES4, and HEYL (Figure S5) appeared to be less expressed in the MYCN-amplified cell lines. Further, all the NOTCH receptor genes except NOTCH4 demonstrated lower mRNA levels in the MYCN-amplified cell lines (Figure S5). In sum, these data indicate that MYCN amplification correlates with the reduced expression of genes encoding the γ-secretase complex, ADAM17, and genes that encode components of the NOTCH signaling pathway in MYCN-amplified NB cell lines.
To further document the role of MYCN in regulating the identified genes and processes, we analyzed publicly available gene expression data following the shRNA-mediated MYCN knockdown in the MYCN-amplified cell line IMR32 [42]. First, we identified 6677 differentially regulated genes (FDR ≤ 0.05), of which 3229 were up-regulated and 3448 were down-regulated following MYCN depletion (Table S3). Pathway analysis of the top 10 biological processes revealed that the up-regulated genes were enriched in genes related to axon guidance, the neurotrophin signaling pathway, endocytosis, and autophagy, as well as the NOTCH signaling pathway (Table S3). The down-regulated genes were related to DNA replication, the cell cycle, ribosome biogenesis, RNA transport, spliceosomes, and DNA repair (Table S3). When overlapping differentially regulated genes following MYCN depletion in the MYCN-amplified IMR32 cell line with our defined MYCN targets based on MYCN ChIP-Seq, we identified 1004 up-regulated genes and 1623 down-regulated genes as MYCN targets in IMR32 cells (Figure S6 and Table S3). Notably, the genes identified in this study involved in NOTCH signaling, such as APH1A, APH1B, and PSEN1, were among the MYCN-bound genes belonging to the up-regulated genes in response to MYCN depletion in the IMR32 cell line (Figure S6A and Table S3), while the MYCN-bound genes related to DNA replication, the cell cycle, RNA splicing, ribosomal biogenesis, and DNA repair were down-regulated (Figure S6B and Table S3). Moreover, the use of more stringent conditions (FDR ≤ 0.05 combined with FC ≥ 1) identified 692 up-regulated genes and 606 down-regulated genes (Table S3). A pathway analysis of these genes revealed the up-regulation of genes related to the NOTCH signaling pathway and axon guidance, while the down-regulated genes were involved in DNA replication and repair, as well as the cell cycle (Figure 2E). These data suggest that MYCN, in addition to its activation of growth-related genes, is a potential direct suppressor of genes that belong to the γ-secretase complex and other genes involved in the NOTCH pathway.

2.3. MYCN Amplification Correlates with Reduced Expression of Genes of the γ-Secretase Complex and ADAM17 in Neuroblastoma

To strengthen the analysis of our bioinformatic approach demonstrating the down-regulation of genes encoding γ-secretase and ADAM17 in the 39-human NB-cell-line dataset (Figure 3A and Figure S4), we sought to evaluate the baseline expression levels of genes encoding γ-secretase and ADAM17 in a set of human NB cell lines representative of MYCN-amplified and non-MYCN-amplified NB. First, we assessed the MYCN and c-MYC protein levels in the NB cell lines used in this study to functionally validate the gene expression levels. Western blot analysis showed that the MYCN-amplified NB cell lines IMR32, KELLY, SK-N-BE(2), and CHP212 expressed high levels of the MYCN protein, while the non-MYCN-amplified NB cell lines expressed either c-MYC (SH-SY5Y and SK-N-AS) or a low level of MYCN (SK-N-FI) relative to the other cell lines (Figure S7). Next, to validate the impact of MYCN amplification on the expression of genes encoding γ-secretase complex subunits and ADAM17, we analyzed the baseline mRNA levels of these genes in those human NB cell lines by quantitative real-time PCR (qRT-PCR). The transcript levels of genes encoding the γ-secretase complex subunits and ADAM17 were found to be significantly lower in the MYCN-amplified NB cell lines IMR32, KELLY, SK-N-BE(2), and CHP212 compared with the non-MYCN-amplified NB cell lines SK-N-FI, SH-SY5Y, and SK-N-AS (Figure 3B). To test whether the expression patterns of these genes are similar in primary human NB tumors as well, we analyzed their expression using the gene expression dataset from the Tumor Neuroblastoma—SEQC study (r2.amc.nl/; Tumor Neuroblastoma-SEQC-498-custom-ag44kcwolf). Like our gene expression analysis on human NB cell lines, we found the expression of genes encoding members of the γ-secretase complex and ADAM17 to be significantly lower in MYCN-amplified primary NB tumors when compared with non-MYCN-amplified tumors (Figure 3C). Notably, the low expressions of NCSTN, APH1B, and PSEN1 encoding for members of the γ-secretase complex correlated with poor overall survival in the analyzed cohort of 498 NB patients (Figure 3D), while other members were not indicative of any changes in overall survival.

2.4. Genetic and Pharmacological Depletion of MYCN Increases mRNA Levels of γ-Secretase Subunits and ADAM17 in Neuroblastoma

Having identified MYCN as a potential direct suppressor of the expression of ADAM17 and γ-secretase complex encoding genes, we sought to test whether MYCN genetic or pharmacological depletion in NB would increase their expression. To this end, we utilized the TET21N NB cell line, in which we could regulate the MYCN expression by doxycycline. As shown in Figure 4A, treatment of the TET21N cell line with 1 µg/mL doxycycline for 72 h completely depleted the MYCN protein levels (Figure 4A). Doxycycline-mediated MYCN depletion in the TET21N cells reduced cell proliferation and metabolic activity without affecting cell viability (Figure S8). Analysis of the mRNA levels of genes encoding the γ-secretase complex and ADAM17 by RT-qPCR revealed a significant increase in the mRNA levels following doxycycline treatment (Figure 4B). Next, we analyzed the impact of MYCN depletion using the newly described MYC:MAX inhibitor MYCMI-7 [43,44]. Utilizing the in situ proximity ligation assay (isPLA), we showed that 5 µM of MYCMI-7 inhibited endogenous MYCN:MAX interactions in the MYCN-amplified KELLY cell line 5 h post-treatment (Figure 4C,D). Treatment for 72 h with 0.4 µM MYCMI-7 reduced MYCN protein levels in the MYCN-amplified NB cell line KELLY (Figure 5E). Moreover, MYCMI-7 reduced cell proliferation and viability and induced cell death, as shown by the increased levels of cleaved-PARP protein in the KELLY cell line (Figure S8). Notably, MYCMI-7 treatment induced the expression of the γ-secretase complex encoding genes, as well as ADAM17, in the MYCN-amplified KELLY cell line (Figure 4F). To exclude effects of the MYCN status on the total RNA synthesis per cell, we measured the total RNA concentration from equal numbers of MYCN-proficient or -deficient cells. As shown in Supplementary Figure S9, the total RNA amount did not differ between the MYCN-expressing and MYCN-depleted TET21N in response to doxycycline treatment or as a consequence of the MYCMI-7 treatment of KELLY cells (Figure S9). These data indicate that MYCN depletion induces the expression of the γ-secretase complex and ADAM17 in MYCN-amplified NB cells, and they further substantiate our bioinformatic analysis suggesting that MYCN plays a role as a direct negative regulator of genes encoding the γ-secretase complex and ADAM17.

2.5. MYCN Depletion Induces Expression of NOTCH Target Genes in Neuroblastoma

Our analysis of gene expression datasets derived from human NB cell lines and primary tumors suggested that several components of the NOTCH pathway are the down-regulated in MYCN-amplified tumors. Thus, we investigated the impact of MYCN amplification on the expression of some canonical NOTCH target genes of the bHLH (basic-helix-loop-helix) family, such as HES (Hairy/enhancer of Split) and HEY (Hairy/Enhancer of Split related with YRPW motif) [45,46]. Our analysis of RNA-Seq data derived from the human NB cell lines uncovered significantly reduced expression of the NOTCH target genes HES1, HES4, HEY1, and HEYL in the MYCN-amplified cell lines (Figure 5A and Figure S5). In primary NB tumors, we also found that HES1, HES4, HEY1 and HEYL were down-regulated in MYCN-amplified tumors when compared to non-MYCN-amplified tumors (Figure 5B). Notably, the low expression of these NOTCH target genes correlated with poor overall survival in NB patients (Figure 5C). To assess the impact of MYCN depletion on the expression of NOTCH target genes, we evaluated the expression of HES1, HES4, HEY1, and HEYL in response to the genetic or pharmacological depletion of MYCN. First, we evaluated the mRNA levels of these genes in the TET21N cell line following MYCN depletion in response to doxycycline treatment. Treatment of the TET21N cells with 1 µg/mL doxycycline for 72 h led to significant increases in the mRNA levels of HES1, HES4, and HEY1 (Figure 5D). We were unable to detect any HEYL mRNA signals in the TET21N cell line. To further corroborate these findings, we analyzed the mRNA levels of these four genes in response to 72 h treatment with 0.4 µM MYCMI-7 in the KELLY cell line. Similarly, the MYCN pharmacological depletion by MYCMI-7 induced significant increases in the HES1, HES4, HEY1 and HEYL mRNA levels when compared to the DMSO control treatment (Figure 5E).
To obtain better insights into the expression of the NOTCH signaling pathway components in NB, we also analyzed the mRNA levels of all the NOTCH receptors (NOTCH1, NOTCH2, NOTCH3, and NOTCH4) in the 39-human NB-cell-line dataset. As shown in Supplementary Figure S10, we did not find significant differences in the NOTCH receptor mRNA levels between the MYCN-amplified and non-MYCN-amplified cell lines (Figure S10A). However, NOTCH receptors 1–3 were less expressed in the MYCN-amplified neuroblastoma cell lines KELLY, SK-N-BE(2)C, and NGP compared with the non-MYCN-amplified SK-N-AS cell line in another gene expression dataset with spike-in RNA control (Figure S5). All NOTCH receptors were significantly lower in expression in primary tumors with MYCN amplification when compared to non-MYCN-amplified tumors (Figure S10B). Notably, the low expressions of the NOTCH receptors were indicative of poor overall survival in NB patients (Figure S10C). In sum, these data indicate that MYCN-amplified NB cell lines and primary tumors are characterized by the low expression of NOTCH target genes, and that the inhibition of MYCN activates their expression.

2.6. Inhibition of γ-Secretase Complex Abolishes Induced Expression of NOTCH Target Genes upon MYCN Depletion in Neuroblastoma

Having shown that MYCN depletion induces the expression of γ-secretase complex and NOTCH target genes, we sought to investigate whether the chemical inhibition of the γ-secretase complex abolishes the induced expression of the NOTCH target genes in response to MYCN depletion. We utilized N-[N-(3,5-difluorophenacetyl)-l-alanyl]-S-phenylglycine t-butyl ester (DAPT), which is an inhibitor of γ-secretase that blocks the proteolytic processing and release of the NICD [47,48]. To this end, we treated the TET21N cell line with 1 µg/mL doxycycline, 5 µM DAPT, or their combination for 72 h, and then analyzed the mRNA expression levels of the NOTCH target genes HES1, HES4, and HEY1. Similar to our earlier observations, the doxycycline-mediated depletion of MYCN significantly increased the expression of HES1, HES4, and HEY1 when compared to the DMSO control treatment (Figure 6). DAPT treatment alone did not affect the mRNA levels of HES1 and HES4 but reduced HEY1 mRNA levels when compared to the DMSO control treatment (Figure 6). Notably, DAPT treatment in combination with doxycycline abrogated the increases in HES1, HES4, and HEY1 mRNA levels in response to MYCN depletion, resulting in mRNA expression levels of these genes that are comparable to the DMSO control levels (Figure 6). In conclusion, these findings suggest that the up-regulation of NOTCH target genes is mediated by the increased expression of ADAM17 and γ-secretase upon the biological depletion of MYCN.

3. Discussion

High-risk neuroblastoma (NB) is a fatal pediatric cancer that is characterized by aggressive and undifferentiated tumors, drug resistance, and poor prognosis. Therefore, better understanding the molecular mechanisms underlying high-risk NB tumor development, progression, and resistance to therapy is a necessity for improved diagnosis, management, and treatment of this deadly tumor [49,50,51]. Amplification of the MYCN oncogene is a well-established marker of poor prognosis and high-risk NB [11,12]. MYCN, as well as the other MYC family of oncoproteins, function as central hubs for almost all signaling pathways; thus, they modulate and integrate a broad range of pathways involved in fundamental biological processes, such as cell proliferation, cell death, pluripotency, differentiation, senescence, cellular plasticity, cellular energetics, and metabolism [52,53,54]. Since the discovery of MYCN amplification in NB, extensive research efforts have been undertaken to uncover the role of MYCN in NB and its contribution to high-risk tumors. These efforts have connected deregulated MYCN expression to all hallmarks of NB, such as proliferation, blocked differentiation, metastatic spread, stemness, metabolic re-wiring, the suppression of immunosurveillance mechanisms, and the emergence of drug resistance [55,56,57,58,59]. However, the direct correlation or causation between MYCN and the NB-relevant biological pathways in MYCN-amplified tumors has not been completely deciphered.
In this study, we implemented an integrative approach using publicly available gene expression datasets derived from large numbers of human NB cell lines to define the biological pathways that are specific to MYCN-amplified tumors when compared to non-MYCN-amplified tumors. This approach identified ribosomal biogenesis, ribosomal RNA synthesis and processing, as well as protein synthesis as the enriched pathways in MYCN-amplified cell lines. Supporting this, we also found a positive correlation between these pathways and MYCN amplification in the RNA expression data of primary NB tumors. This is in line with previous studies demonstrating the dependency of MYCN-amplified tumors on these pathways. Based on MYCN overexpression studies on NB cell lines, Boon and colleagues [60] demonstrated that MYCN drives the expression of a large set of genes involved in ribosome biogenesis and protein synthesis [60]. Similarly, a strong correlation between MYCN expression and ribosome biogenesis has been documented in neuroblastoma patients, with a negative impact on patient survival [61]. It has been suggested that the factors involved in ribosomal biogenesis and ribosomal RNA processing are essential factors in the pathogenesis of MYCN-amplified NB tumors, and their genetic depletion demonstrates anti-NB effects [62,63,64,65]. We demonstrate that shRNA-mediated knockdown of MYCN in the MYCN-amplified cell line IMR32 resulted in down-regulated expression levels of genes related to RNA processing and ribosomal RNA synthesis and biogenesis, many of which were defined as MYCN-target genes in this study. By integrating MYCN ChIP-Seq data, our study suggests the direct regulation of genes involved in ribosomal biogenesis and RNA synthesis, as well as protein synthesis, by MYCN in NB. Furthermore, there was an enrichment of the MYC-type of E-box motifs around the TSSs of the up-regulated genes, further supporting the potential direct regulation of these genes by MYCN, and in agreement with previous reports [26,66]. Moreover, our analysis suggests that targeting ribosomal biogenesis could be a promising therapeutic strategy in high-risk NB with MYCN amplification. Indeed, blocking ribosome RNA expression using inhibitors of RNA polymerase I has demonstrated potent anti-NB activities in vitro and in vivo [61,67], which warrants further investigation and evaluation in clinical trials.
Our bioinformatic analysis of gene expression datasets revealed down-regulated expression of genes involved in the NOTCH signaling pathway in MYCN-amplified cell lines and primary tumors. Furthermore, our MYCN ChIP-Seq analysis revealed that genes encoding ADAM17 and subunits of the γ-secretase complex are direct MYCN target genes. Thus, this study focused on evaluating the impact of MYCN on the expression of ADAM17 and the γ-secretase complex genes, as well as the expression of some NOTCH target genes in NB. It has been shown that NB cells express variable levels of NOTCH receptors, but low levels of NOTCH target genes [68], suggesting impaired NOTCH pathway activity in NB. Notably, Van Limpt, V.A., et al. [69] analyzed the expression of 21 genes involved in the NOTCH pathway in a large panel of NB cell lines, and the gene expression appeared to be the highest in the non-MYCN-amplified cell lines, such as SK-N-FI, LAN-6, SHEP, and LAN-2 [69]. Our findings demonstrating low expression of NOTCH receptors, ADAM17, γ-secretase, and NOTCH target genes in MYCN-amplified tumors, suggest that MYCN amplification contributes to the down-regulation of genes related to the NOTCH pathway in NB. This notion is supported by the demonstration that MYCN depletion induces the expression of genes encoding ADAM17, γ-secretase, and some NOTCH target genes. Further, we uncovered that the TSSs of ADAM17 and γ-secretase encoding genes are demarcated by MYCN, suggesting that MYCN is directly involved in the repression of these genes. It has been suggested that MYC/MYCN-mediated gene repression is dependent on interactions with MIZ-1 [26,29,70,71], epigenetic factors such as EZH2 histone methyltransferase and polycomb repressive complexes [71,72,73,74], histone deacetylase complexes [71,75], and DNA methyltransferase complexes [76,77]. Whether such mechanisms contribute to the reduced expression of genes encoding ADAM17, γ-secretase, and some NOTCH targets in MYCN-amplified tumors warrants further investigations. Studies in medulloblastoma have indicated that MYCN binds MIZ-1 with a lower affinity compared with MYC, with consequences for MIZ-dependent repression in different subtypes of medulloblastoma [78]. Possibly, MYCN-mediated transcriptional repression is therefore less dependent on MIZ-1. Interestingly, we found that MYCN-bound, repressed genes were enriched with DNA-binding motifs belonging to the ETS family of transcription factors, and, to a lesser extent, with the MYC-type of E-box motifs, compared with activated genes, suggesting that MYCN and ETS-family proteins may cooperate in the repression of transcription in neuroblastoma. [26,27,66]. Notably, The ETS family member ELK4 has been described as one of the factors that may drive the MYCN-mediated repression of gene transcription in NB [66]. Another ETS family member, ELK1, has been shown to negatively regulate the expression of some metalloproteases and members of the γ-secretase complex, such as PSEN1, PSEN2, and APH1A, when overexpressed in NB cells [79]. Furthermore, the γ-secretase complex is known for its wide range of substrates that act in many signal pathways [80,81]. Thus, it would be interesting to evaluate the impact of ADAM17 and the γ-secretase complex on substrates other than NOTCH receptors in MYCN-amplified tumors, and to define their biological contribution in the pathogenesis of such tumors. In addition, inhibitors targeting the γ-secretase complex have been suggested as potential anti-cancer agents [47,82,83]. The fact that our study demonstrates the low gene expression of the γ-secretase complex may argue against the therapeutic potential of γ-secretase inhibitors in MYCN-amplified NB tumors, which indeed demands further investigation.
It has been suggested that the NOTCH signaling pathway plays both oncogenic [84,85,86] and tumor suppressor [68] roles in NB. It is worth mentioning that reports describing the oncogenic functions of NOTCH signaling in NB are based on overexpression experiments, the use of chemical inhibitors, and in most cases, the utilization of non-MYCN-amplified cell lines, such as SH-SY5Y, which expresses c-MYC [84,85,86]. Cooperation between c-MYC and NOTCH during tumorigenesis has been documented in lung cancer [87], T-cell acute lymphocytic leukemia [88], and mantle cell lymphoma and chronic lymphocytic leukemia [89]. Whether NOTCH-c-MYC oncogenic cooperation occurs in non-MYCN-amplified NB with c-MYC expression demands further investigation. Recent studies investigating NB heterogeneity and cell identity based on core regulatory circuits (CRCs) and super enhancers (SEs) divided human NB cell lines and primary tumors into two main categories: undifferentiated mesenchymal neural crest-like NB (MES) and committed adrenergic cells (ADRN) [90,91]. Notably, MYCN-amplified cell lines, except for the CHP212 cell line, have demonstrated an enrichment of ADRN identity, while non-MYCN-amplified cell lines represented both identities [91]. A CRC and SE analysis of MES and ADRN in NB cell lines demonstrated an enrichment of NOTCH pathway genes in the MES but not the ADRN identity cell lines [92]. Even though MYC family members share common target genes and regulate similar pathways, and high c-MYC levels may drive high-risk NB [93], distinct transcriptional activities of MYCN and c-MYC have been suggested in NB [94]. These findings as well as our study demonstrating the reduced expression of genes involved in the NOTCH pathway in MYCN-amplified cell lines and primary tumors, point toward a possible antagonistic relationship between MYCN and the NOTCH pathway components in these NB tumors, but not in the non-MYCN-amplified NB expressing MYC. Intriguingly, high MYC expression in group 3 medulloblastoma has been shown to inversely correlate with the NOTCH signaling signature and reflects early fetal cerebellar development [95]. Therefore, the impact of NOTCH signaling on NB biology needs to be addressed carefully, as the biological outcomes of NOTCH signaling might be different depending on the subtype of NB (MYCN-amplification status), different experimental settings (physiological expression levels vs. overexpression), developmental stage and cell of origin, as well as the contribution of the tumor microenvironment (ligand activation vs. overexpression of NICD).
In summary, our study uncovers novel pathways directly regulated by MYCN in NB, and MYCN amplification correlates with the reduced expression of genes involved in the NOTCH signaling pathway compared to non-MYCN-amplified tumors. We also suggest that MYCN may function as a direct regulator dampening gene transcription of ADAM17 and the γ-secretase complex genes in NB. However, further studies are required for a deeper understanding of the interplay between MYCN and the components of the NOTCH signaling pathway in NB.

4. Materials and Methods

4.1. RNA-Seq, ChIP-Seq Data, and Gene Ontology (GO) Analysis

The RNA-Sequencing (RNA-Seq) data of the 39 neuroblastoma cell lines under the accession number GSE89413 were downloaded from the gene expression Omnibus (GEO) database. Differentially expressed genes in MYCN-amplified vs. non-MYCN-amplified cell lines were analyzed using R scripts available at https://github.com/marislab/NBL-cell-line-RNA-seq, accessed on 23 February 2020, as described previously [40]. The raw sra files were downloaded using NCBI sratools. The reads were aligned using STAR version 2.7.0a [96]. The hg19 genome was used to align the reads. The aligned BAM files were then used to make the Tag files using the HOMER v4 module makeTagDirectory [97]. The HOMER module analyzeRepeats.pl was further used to generate the raw count file, which was used as an input to the R script provided in the original manuscript [40]. Data normalization and differential analysis were performed using the tool DESeq2 v 1.38.3 [98]. An adjusted p-value ≤ 0.05 was used as a final stringent measure to define the differentially expressed genes. MYCN ChIP-Seq data for the KELLY, COGN415, and LAN5 cell lines were downloaded from GEO under the accession numbers GSE94782 and GSE138295 [41]. All raw files and input data can be found at the Sequence Read Archive under the SRA accessions SRP223942. The downloaded data were aligned using BWA version0.7.17-r1188, and the human genome assembly hg19 was used as a reference. MACS peak-caller was used to call the peaks using the default parameters [99]. The peaks were annotated using the online tool GREAT v4.04 [100]. A cutoff of ±5 kb of the transcription start site (TSS) was used to define genes as MYCN-target genes. The molecular signatures database (MSigDB) and Enrichr online gene expression analysis tool [101,102,103] were utilized to investigate the gene ontologies of the differentially regulated genes in the 39-cell-line RNA-Seq data, and MYCN-bound genes were defined. The top 10 differentially regulated biological processes with a false discovery rate (FDR) ≤ 0.05 are presented in this study.

4.2. R2: Genomics Analysis and Visualization Platform

The R2: Genomics Analysis and Visualization Platform (http://r2.amc.nl, accessed on 15 March 2022) was used to analyze the expression of the NOTCH signaling pathway genes in primary NB tumors with or without MYCN amplification in the SEQC (GSE49710) dataset, and for the analysis of the spike-in gene expression arrays of the KELLY, SK-N-BE(2)c, NGP, and SK-N-AS cell lines (GSE80149). The R2 genome browser platform was utilized to show the MYCN enrichment at the transcription start sites of selected MYCN-regulated genes relevant to the findings in this paper. Statistical analysis was performed based on the R2: Genomics Analysis and Visualization Platform recommended setting, with an FDR ≤ 0.05 as the statistical significance threshold.

4.3. Cell Lines and Cell Cultures

SK-N-FI (ATCC catalog no. CRL-2142), SK-N-AS (ATCC catalog no. CRL-2137), KELLY (DSMZ catalog no. ACC-355), and SK-N-BE(2) (ATCC catalog no. CRL-2271) were kindly provided by Professor Per Kogner at the Karolinska Institutet. The IMR32 (ATCC catalog no. CCL-127), SH-SY5Y (ATCC catalog no. CRL-2266), and CHP-212 (ATCC catalog no. CRL-2273) cell lines were kindly provided by Professor Marie Arsenian Henriksson at the Karolinska Institutet. The TET21N cell line was obtained as a gift from Werner Lutz [104]. The MYCN-amplified cell lines IMR32, KELLY, SK-N-BE(2), and CHP212, as well as the TET12N cell line with doxycycline-regulatable MYCN expression, were maintained in RPMI 1640 medium (Gibco, ThermoFisher Scientific, Waltham, MA, USA). The SK-N-FI, SK-N-AS, and SH-SY5Y non-MYCN-amplified cell lines were maintained in DMEM medium (ThermoFisher Scientific, Waltham, MA, USA). The cell culture media contained GlutaMAX and were supplemented with 10% fetal bovine serum (ThermoFisher Scientific, Waltham, MA, USA) and 1% penicillin/streptomycin (ThermoFisher Scientific, Waltham, MA, USA). The IMR32 cell line culture was supplemented with nonessential amino acids (ThermoFisher Scientific, Waltham, MA, USA). Cell lines were grown in a humidifier incubator at 37 °C in 5% CO2. Treatment of cells with inhibitors or doxycycline was performed for 72 h. In brief, ≥90% viable cells were seeded in 10 cm dishes overnight and then treated with the indicated concentrations for 72 h. DMSO was used as the control treatment. All human neuroblastoma cell lines were STR-authenticated by short tandem repeat (STR) analysis (Eurofins Genomics, Ebersberg, Germany) once acquired by our laboratory. All cell lines were routinely tested for mycoplasma using the MycoAlert Mycoplasma Detection Kit (Lonza, Basel, Switzerland).

4.4. Resazurin Assay and Assessment of Cell Proliferation and Viability

For the resazurin assay assessment of the metabolic activity, the TET12N and KELLY cell lines were seeded in 96-well plates at a 3000 cells/well cell density and incubated overnight in a humidifier incubator at 37 °C in 5% CO2. The TET21N cell line was treated with 1 µg/mL doxycycline, while the KELLY cell line was treated with 0.4 µM MYCMI-7 for 72 h. DMSO was used as the control treatment. Then, resazurin reagent (Sigma, St. Louis, MO, USA) was added at a 44 µM final concentration. Reduced resazurin was detected by measuring the absorbance at 570 nm using a 96-well plate reader. A difference of 10× between the blank and DMSO-control-treated cells was set as the threshold. The cell proliferation was evaluated by counting the cell numbers over the treatment period, and the viability was assessed at the same time using trypan blue staining.

4.5. RNA Extraction, cDNA Synthesis, and RT-qPCR

RNA was isolated using the TRIzol extraction method following the manufacturer protocol (ThermoFisher Scientific, Waltham, MA, USA). The RNA concentration was measured using a Nanodrop, and 1 µg of RNA was used for cDNA synthesis using iScript™ Reverse Transcription Supermix (BioRad, Hercules, CA, USA), following the manufacturer protocol. Real-time quantitative PCR (RT-qPCR) was performed using iTaq™ Universal SYBR® Green Supermix (BioRad, Hercules, CA, USA) in a 384-well PCR plate format. RT-qPCR reactions were performed using the C1000 Touch™ Thermal Cycler with the CFX384™ Real-time System (BioRad, Hercules, CA, USA). Analysis of the gene expression was performed using Bio-Rad CFX Maestro™ data analysis software version 1.1 (BioRad, Hercules, CA, USA) following the 2−∆∆Ct method. HPRT1 and GAPDH were used as the housekeeping genes for the normalization of the RT-qPCR analysis. The primers used in this study are listed in Supplementary Table S4.

4.6. Protein Extraction and Western Blot

Cell lines were collected and washed with 1× cold phosphate-buffered saline (PBS). Cells were lysed with radioimmunoprecipitation assay RIPA buffer (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with a protease and phosphatase inhibitor cocktail (Thermo Fisher Scientific) for 30 min on ice, sonicated (amplitude: 80; three cycles at 30 s each), and centrifuged (12,000× g for 15 min, 4 °C). The protein concentration was determined with a Pierce BCA Protein Assay (Thermo Fisher Scientific, Waltham, MA, USA). Equal amounts of protein (20 μg per sample) were separated on 4–12% Bis–Tris gradient gels (Thermo Fisher Scientific, Waltham, MA, USA) and blotted on a 0.2 µm nitrocellulose membrane (Bio-Rad, Hercules, CA, USA) with the Trans-Blot Turbo Transfer System (Bio-Rad, Hercules, CA, USA). Membranes were blocked for 1 h at room temperature, in 5% milk or bovine serum albumin (BSA) in PBS-Tween20 (PBS-T), followed by an overnight incubation with primary antibody at 4 °C. The following primary antibodies were used: anti-β-actin (A5441, Sigma-Aldrich, St. Louis, MO, USA) and anti-MYCN (sc-53993, Santa Cruz, CA, USA). Membranes were washed with PBS-T and incubated with secondary horseradish peroxidase-conjugated anti-mouse (ab97046, Abcam, Cambridge, UK) or anti-rabbit (ab97080, Abcam) antibody for 1 h at room temperature. To visualize the proteins, membranes were developed with Western Enhanced Chemiluminescence Substrates (Bio-Rad, Hercules, CA, USA). The ChemiDoc XRS+ System (Bio-Rad, Hercules, CA, USA) was used to determine the band intensity, along with Image Lab Software 5.1 for analysis.

4.7. In Situ Proximity Ligation Assay (isPLA)

The in-situ proximity ligation assay (isPLA) was performed as described previously elsewhere [37,80]. In brief, the MYCN-amplified KELLY cell line was used to detect endogenous MYCN:MAX interactions. KELLY cells were seeded in 96-well plates for 48 h, then treated with 5 µM of MYCMI-7 or DMSO for 5 h, and then washed twice with PBS and fixed with 4% paraformaldehyde for 10 min at room temperature. Cells were permeabilized with PBS with 0.05% Triton-X and incubated in a blocking buffer for 1 h at 37 °C, after which the isPLA was performed using the NaveniFlex MR kit following the manufacturer protocol (Navinci Diagnostics, Uppsala, Sweden). The following antibodies were used: mouse monoclonal anti-MYCN (sc-53993, Santa Cruz, CA, USA) and rabbit polyclonal anti-MAX (Abcam, catalog no. ab101271). Cells were stained with phalloidin and DAPI, and the isPLA signals were developed using Buffer C provided by the isPLA kit with Atto 647N. Image acquisition was performed using the ZEISS LSM 980 with an Airyscan 2 confocal microscope (Carl Zeiss Microscopy GmbH, Munich, Germany). All image stacks were acquired with comparable settings, at a resolution of 1024 × 1024 pixels and a z-stack size of 3 μm. The isPLA signals were quantified using CellProfiler 4.2.1 software. Data were analyzed using Rstudio version 4.1.0 with the following packages: tidyverse and ggplot2.

4.8. Statistical Analysis

Statistical analysis was performed using a two-tailed, unpaired Student t-test for the RT-qPCR data, and a one-way ANOVA for a comparison of the gene expression levels between the MYCN-amplified and non-MYCN-amplified cell lines and primary tumors, as well as for the isPLA data. Bonferroni correction of raw p-values was used for the Kaplan–Meier-analysis overall survival curve.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms24098141/s1.

Author Contributions

Conceptualization, M.A.; methodology, P.A., A.G., L.M. and M.A.; validation, P.A., A.G., L.M., W.B. and M.A.; formal analysis, P.A., A.G., L.M., W.B., L.-G.L. and M.A.; resources, M.A. and L.-G.L.; writing—original draft preparation, M.A.; writing—review and editing, P.A., A.G., L.M., W.B., L.-G.L. and M.A.; visualization, P.A., A.G., L.M., W.B. and M.A.; supervision, M.A.; project administration, M.A.; funding acquisition, L.-G.L. and M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Swedish Childhood Cancer Society (M.A. TJ2019-0018 and PR2019-0063; L.-G.L. PR2020-0141 and PR2022-0130), OE och Edla Johanssons vetenskapliga stiftelse (M.A. 2020-MA, 2021-MA, and 2022-MA), Lillian Sagens och Curt Ericssons Forskningsstiftelse (M.A. 2022-MA), the Swedish Cancer Society (L.-G.L. CAN 21 1817 Pj), Novo Nordisk Foundation (L.-G.L. NNF20SA0066685), KI Foundations (L.-G.L. 2020-02157), and MyCural Therapeutics (L.-G.L. MTRCA-001A2023).

Acknowledgments

The authors want to thank Marie Arsenian Henriksson, Per Kogner, and Werner Lutz for providing the cell lines. We also want to thank Urban Lendahl, Eike-Benjamin Braune, and Francesca Del Gaudio for taking the time to read and comment on the manuscript. The authors express their gratitude to all the funding agencies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maris, J.M.; Hogarty, M.D.; Bagatell, R.; Cohn, S.L. Neuroblastoma. Lancet 2007, 369, 2106–2120. [Google Scholar] [CrossRef]
  2. Maris, J.M. Recent advances in neuroblastoma. N. Engl. J. Med. 2010, 362, 2202–2211. [Google Scholar] [CrossRef]
  3. Matthay, K.K.; Maris, J.M.; Schleiermacher, G.; Nakagawara, A.; Mackall, C.L.; Diller, L.; Weiss, W.A. Neuroblastoma. Nat. Rev. Dis. Primers 2016, 2, 16078. [Google Scholar] [CrossRef] [PubMed]
  4. Van Roy, N.; De Preter, K.; Hoebeeck, J.; Van Maerken, T.; Pattyn, F.; Mestdagh, P.; Vermeulen, J.; Vandesompele, J.; Speleman, F. The emerging molecular pathogenesis of neuroblastoma: Implications for improved risk assessment and targeted therapy. Genome Med. 2009, 1, 74. [Google Scholar] [CrossRef]
  5. Newman, E.A.; Nuchtern, J.G. Recent biologic and genetic advances in neuroblastoma: Implications for diagnostic, risk stratification, and treatment strategies. Semin. Pediatr. Surg. 2016, 25, 257–264. [Google Scholar] [CrossRef] [PubMed]
  6. Cohn, S.L.; Pearson, A.D.; London, W.B.; Monclair, T.; Ambros, P.F.; Brodeur, G.M.; Faldum, A.; Hero, B.; Iehara, T.; Machin, D.; et al. The International Neuroblastoma Risk Group (INRG) classification system: An INRG Task Force report. J. Clin. Oncol. 2009, 27, 289–297. [Google Scholar] [CrossRef] [PubMed]
  7. Bosse, K.R.; Maris, J.M. Advances in the translational genomics of neuroblastoma: From improving risk stratification and revealing novel biology to identifying actionable genomic alterations. Cancer 2016, 122, 20–33. [Google Scholar] [CrossRef]
  8. Matthay, K.K.; Villablanca, J.G.; Seeger, R.C.; Stram, D.O.; Harris, R.E.; Ramsay, N.K.; Swift, P.; Shimada, H.; Black, C.T.; Brodeur, G.M.; et al. Treatment of high-risk neuroblastoma with intensive chemotherapy, radiotherapy, autologous bone marrow transplantation, and 13-cis-retinoic acid. Children’s Cancer Group. N. Engl. J. Med. 1999, 341, 1165–1173. [Google Scholar] [CrossRef]
  9. Louis, C.U.; Shohet, J.M. Neuroblastoma: Molecular pathogenesis and therapy. Annu. Rev. Med. 2015, 66, 49–63. [Google Scholar] [CrossRef]
  10. Whittle, S.B.; Smith, V.; Doherty, E.; Zhao, S.; McCarty, S.; Zage, P.E. Overview and recent advances in the treatment of neuroblastoma. Expert. Rev. Anticancer. Ther. 2017, 17, 369–386. [Google Scholar] [CrossRef] [PubMed]
  11. Brodeur, G.M.; Seeger, R.C.; Schwab, M.; Varmus, H.E.; Bishop, J.M. Amplification of N-myc in untreated human neuroblastomas correlates with advanced disease stage. Science 1984, 224, 1121–1124. [Google Scholar] [CrossRef]
  12. Seeger, R.C.; Brodeur, G.M.; Sather, H.; Dalton, A.; Siegel, S.E.; Wong, K.Y.; Hammond, D. Association of multiple copies of the N-myc oncogene with rapid progression of neuroblastomas. N. Engl. J. Med. 1985, 313, 1111–1116. [Google Scholar] [CrossRef]
  13. Gabay, M.; Li, Y.; Felsher, D.W. MYC activation is a hallmark of cancer initiation and maintenance. Cold Spring Harb. Perspect. Med. 2014, 4, a014241. [Google Scholar] [CrossRef]
  14. Stine, Z.E.; Walton, Z.E.; Altman, B.J.; Hsieh, A.L.; Dang, C.V. MYC, Metabolism, and Cancer. Cancer Discov. 2015, 5, 1024–1039. [Google Scholar] [CrossRef]
  15. Dejure, F.R.; Eilers, M. MYC and tumor metabolism: Chicken and egg. Embo J. 2017, 36, 3409–3420. [Google Scholar] [CrossRef]
  16. Rickman, D.S.; Schulte, J.H.; Eilers, M. The Expanding World of N-MYC-Driven Tumors. Cancer Discov. 2018, 8, 150–163. [Google Scholar] [CrossRef] [PubMed]
  17. Swier, L.; Dzikiewicz-Krawczyk, A.; Winkle, M.; van den Berg, A.; Kluiver, J. Intricate crosstalk between MYC and non-coding RNAs regulates hallmarks of cancer. Mol. Oncol. 2019, 13, 26–45. [Google Scholar] [CrossRef] [PubMed]
  18. Dang, C.V. MYC on the path to cancer. Cell 2012, 149, 22–35. [Google Scholar] [CrossRef] [PubMed]
  19. Blackwell, T.K.; Kretzner, L.; Blackwood, E.M.; Eisenman, R.N.; Weintraub, H. Sequence-specific DNA binding by the c-Myc protein. Science 1990, 250, 1149–1151. [Google Scholar] [CrossRef] [PubMed]
  20. Blackwood, E.M.; Eisenman, R.N. Max: A helix-loop-helix zipper protein that forms a sequence-specific DNA-binding complex with Myc. Science 1991, 251, 1211–1217. [Google Scholar] [CrossRef] [PubMed]
  21. Rahl, P.B.; Lin, C.Y.; Seila, A.C.; Flynn, R.A.; McCuine, S.; Burge, C.B.; Sharp, P.A.; Young, R.A. c-Myc regulates transcriptional pause release. Cell 2010, 141, 432–445. [Google Scholar] [CrossRef] [PubMed]
  22. Zeid, R.; Lawlor, M.A.; Poon, E.; Reyes, J.M.; Fulciniti, M.; Lopez, M.A.; Scott, T.G.; Nabet, B.; Erb, M.A.; Winter, G.E.; et al. Enhancer invasion shapes MYCN-dependent transcriptional amplification in neuroblastoma. Nat. Genet. 2018, 50, 515–523. [Google Scholar] [CrossRef] [PubMed]
  23. Lin, C.Y.; Lovén, J.; Rahl, P.B.; Paranal, R.M.; Burge, C.B.; Bradner, J.E.; Lee, T.I.; Young, R.A. Transcriptional amplification in tumor cells with elevated c-Myc. Cell 2012, 151, 56–67. [Google Scholar] [CrossRef] [PubMed]
  24. Nie, Z.; Hu, G.; Wei, G.; Cui, K.; Yamane, A.; Resch, W.; Wang, R.; Green, D.R.; Tessarollo, L.; Casellas, R.; et al. c-Myc is a universal amplifier of expressed genes in lymphocytes and embryonic stem cells. Cell 2012, 151, 68–79. [Google Scholar] [CrossRef] [PubMed]
  25. Sabò, A.; Kress, T.R.; Pelizzola, M.; de Pretis, S.; Gorski, M.M.; Tesi, A.; Morelli, M.J.; Bora, P.; Doni, M.; Verrecchia, A.; et al. Selective transcriptional regulation by Myc in cellular growth control and lymphomagenesis. Nature 2014, 511, 488–492. [Google Scholar] [CrossRef] [PubMed]
  26. Walz, S.; Lorenzin, F.; Morton, J.; Wiese, K.E.; von Eyss, B.; Herold, S.; Rycak, L.; Dumay-Odelot, H.; Karim, S.; Bartkuhn, M.; et al. Activation and repression by oncogenic MYC shape tumour-specific gene expression profiles. Nature 2014, 511, 483–487. [Google Scholar] [CrossRef]
  27. Tesi, A.; de Pretis, S.; Furlan, M.; Filipuzzi, M.; Morelli, M.J.; Andronache, A.; Doni, M.; Verrecchia, A.; Pelizzola, M.; Amati, B.; et al. An early Myc-dependent transcriptional program orchestrates cell growth during B-cell activation. EMBO Rep. 2019, 20, e47987. [Google Scholar] [CrossRef]
  28. Staller, P.; Peukert, K.; Kiermaier, A.; Seoane, J.; Lukas, J.; Karsunky, H.; Möröy, T.; Bartek, J.; Massagué, J.; Hänel, F.; et al. Repression of p15INK4b expression by Myc through association with Miz-1. Nat. Cell Biol. 2001, 3, 392–399. [Google Scholar] [CrossRef]
  29. Wanzel, M.; Herold, S.; Eilers, M. Transcriptional repression by Myc. Trends Cell Biol. 2003, 13, 146–150. [Google Scholar] [CrossRef]
  30. Soucek, L.; Whitfield, J.; Martins, C.P.; Finch, A.J.; Murphy, D.J.; Sodir, N.M.; Karnezis, A.N.; Swigart, L.B.; Nasi, S.; Evan, G.I. Modelling Myc inhibition as a cancer therapy. Nature 2008, 455, 679–683. [Google Scholar] [CrossRef]
  31. McKeown, M.R.; Bradner, J.E. Therapeutic strategies to inhibit MYC. Cold Spring Harb. Perspect. Med. 2014, 4, a014266. [Google Scholar] [CrossRef]
  32. Castell, A.; Larsson, L.G. Targeting MYC Translation in Colorectal Cancer. Cancer Discov. 2015, 5, 701–703. [Google Scholar] [CrossRef]
  33. Whitfield, J.R.; Beaulieu, M.E.; Soucek, L. Strategies to Inhibit Myc and Their Clinical Applicability. Front. Cell Dev. Biol. 2017, 5, 10. [Google Scholar] [CrossRef]
  34. Allen-Petersen, B.L.; Sears, R.C. Mission Possible: Advances in MYC Therapeutic Targeting in Cancer. BioDrugs 2019, 33, 539–553. [Google Scholar] [CrossRef] [PubMed]
  35. George, S.L.; Parmar, V.; Lorenzi, F.; Marshall, L.V.; Jamin, Y.; Poon, E.; Angelini, P.; Chesler, L. Novel therapeutic strategies targeting telomere maintenance mechanisms in high-risk neuroblastoma. J. Exp. Clin. Cancer Res. 2020, 39, 78. [Google Scholar] [CrossRef] [PubMed]
  36. Koneru, B.; Lopez, G.; Farooqi, A.; Conkrite, K.L.; Nguyen, T.H.; Macha, S.J.; Modi, A.; Rokita, J.L.; Urias, E.; Hindle, A.; et al. Telomere Maintenance Mechanisms Define Clinical Outcome in High-Risk Neuroblastoma. Cancer Res. 2020, 80, 2663–2675. [Google Scholar] [CrossRef] [PubMed]
  37. Brady, S.W.; Liu, Y.; Ma, X.; Gout, A.M.; Hagiwara, K.; Zhou, X.; Wang, J.; Macias, M.; Chen, X.; Easton, J.; et al. Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations. Nat. Commun. 2020, 11, 5183. [Google Scholar] [CrossRef]
  38. Pugh, T.J.; Morozova, O.; Attiyeh, E.F.; Asgharzadeh, S.; Wei, J.S.; Auclair, D.; Carter, S.L.; Cibulskis, K.; Hanna, M.; Kiezun, A.; et al. The genetic landscape of high-risk neuroblastoma. Nat. Genet. 2013, 45, 279–284. [Google Scholar] [CrossRef]
  39. Ackermann, S.; Cartolano, M.; Hero, B.; Welte, A.; Kahlert, Y.; Roderwieser, A.; Bartenhagen, C.; Walter, E.; Gecht, J.; Kerschke, L.; et al. A mechanistic classification of clinical phenotypes in neuroblastoma. Science 2018, 362, 1165–1170. [Google Scholar] [CrossRef]
  40. Harenza, J.L.; Diamond, M.A.; Adams, R.N.; Song, M.M.; Davidson, H.L.; Hart, L.S.; Dent, M.H.; Fortina, P.; Reynolds, C.P.; Maris, J.M. Transcriptomic profiling of 39 commonly-used neuroblastoma cell lines. Sci. Data 2017, 4, 170033. [Google Scholar] [CrossRef]
  41. Upton, K.; Modi, A.; Patel, K.; Kendsersky, N.M.; Conkrite, K.L.; Sussman, R.T.; Way, G.P.; Adams, R.N.; Sacks, G.I.; Fortina, P.; et al. Epigenomic profiling of neuroblastoma cell lines. Sci. Data 2020, 7, 116. [Google Scholar] [CrossRef]
  42. Valentijn, L.J.; Koster, J.; Haneveld, F.; Aissa, R.A.; van Sluis, P.; Broekmans, M.E.; Molenaar, J.J.; van Nes, J.; Versteeg, R. Functional MYCN signature predicts outcome of neuroblastoma irrespective of MYCN amplification. Proc. Natl. Acad. Sci. USA 2012, 109, 19190–19195. [Google Scholar] [CrossRef] [PubMed]
  43. Castell, A.; Yan, Q.; Fawkner, K.; Bazzar, W.; Zhang, F.; Wickström, M.; Alzrigat, M.; Franco, M.; Krona, C.; Cameron, D.P.; et al. MYCMI-7: A Small MYC-Binding Compound that Inhibits MYC: MAX Interaction and Tumor Growth in a MYC-Dependent Manner. Cancer Res. Commun. 2022, 2, 182–201. [Google Scholar] [CrossRef]
  44. Castell, A.; Yan, Q.; Fawkner, K.; Hydbring, P.; Zhang, F.; Verschut, V.; Franco, M.; Zakaria, S.M.; Bazzar, W.; Goodwin, J.; et al. A selective high affinity MYC-binding compound inhibits MYC:MAX interaction and MYC-dependent tumor cell proliferation. Sci. Rep. 2018, 8, 10064. [Google Scholar] [CrossRef] [PubMed]
  45. Kopan, R.; Ilagan, M.X. The canonical Notch signaling pathway: Unfolding the activation mechanism. Cell 2009, 137, 216–233. [Google Scholar] [CrossRef]
  46. Kovall, R.A.; Gebelein, B.; Sprinzak, D.; Kopan, R. The Canonical Notch Signaling Pathway: Structural and Biochemical Insights into Shape, Sugar, and Force. Dev. Cell 2017, 41, 228–241. [Google Scholar] [CrossRef]
  47. Feng, J.; Wang, J.; Liu, Q.; Li, J.; Zhang, Q.; Zhuang, Z.; Yao, X.; Liu, C.; Li, Y.; Cao, L.; et al. DAPT, a γ-Secretase Inhibitor, Suppresses Tumorigenesis, and Progression of Growth Hormone-Producing Adenomas by Targeting Notch Signaling. Front. Oncol. 2019, 9, 809. [Google Scholar] [CrossRef]
  48. Dovey, H.F.; John, V.; Anderson, J.P.; Chen, L.Z.; de Saint Andrieu, P.; Fang, L.Y.; Freedman, S.B.; Folmer, B.; Goldbach, E.; Holsztynska, E.J.; et al. Functional gamma-secretase inhibitors reduce beta-amyloid peptide levels in brain. J. Neurochem. 2001, 76, 173–181. [Google Scholar] [CrossRef]
  49. Volchenboum, S.L.; Cohn, S.L. Progress in defining and treating high-risk neuroblastoma: Lessons from the bench and bedside. J. Clin. Oncol. 2009, 27, 1003–1004. [Google Scholar] [CrossRef] [PubMed]
  50. Saarinen-Pihkala, U.M.; Jahnukainen, K.; Wikström, S.; Koivusalo, A.; Karikoski, R.; Sariola, H.; Hovi, L. Ultrahigh-risk group within the high-risk neuroblastoma category. J. Pediatr. Hematol. Oncol. 2013, 35, e254–e259. [Google Scholar] [CrossRef]
  51. Pinto, N.; Naranjo, A.; Hibbitts, E.; Kreissman, S.G.; Granger, M.M.; Irwin, M.S.; Bagatell, R.; London, W.B.; Greengard, E.G.; Park, J.R.; et al. Predictors of differential response to induction therapy in high-risk neuroblastoma: A report from the Children’s Oncology Group (COG). Eur. J. Cancer 2019, 112, 66–79. [Google Scholar] [CrossRef]
  52. Fagnocchi, L.; Zippo, A. Multiple Roles of MYC in Integrating Regulatory Networks of Pluripotent Stem Cells. Front. Cell Dev. Biol. 2017, 5, 7. [Google Scholar] [CrossRef]
  53. Meyer, N.; Penn, L.Z. Reflecting on 25 years with MYC. Nat. Rev. Cancer 2008, 8, 976–990. [Google Scholar] [CrossRef]
  54. Hann, S.R. MYC cofactors: Molecular switches controlling diverse biological outcomes. Cold Spring Harb. Perspect. Med. 2014, 4, a014399. [Google Scholar] [CrossRef]
  55. Westermark, U.K.; Wilhelm, M.; Frenzel, A.; Henriksson, M.A. The MYCN oncogene and differentiation in neuroblastoma. Semin. Cancer Biol. 2011, 21, 256–266. [Google Scholar] [CrossRef]
  56. Huang, M.; Weiss, W.A. Neuroblastoma and MYCN. Cold Spring Harb. Perspect. Med. 2013, 3, a014415. [Google Scholar] [CrossRef]
  57. Blavier, L.; Yang, R.M.; DeClerck, Y.A. The Tumor Microenvironment in Neuroblastoma: New Players, New Mechanisms of Interaction and New Perspectives. Cancers 2020, 12, 2912. [Google Scholar] [CrossRef] [PubMed]
  58. Otte, J.; Dyberg, C.; Pepich, A.; Johnsen, J.I. MYCN Function in Neuroblastoma Development. Front. Oncol. 2020, 10, 624079. [Google Scholar] [CrossRef] [PubMed]
  59. Gogolin, S.; Dreidax, D.; Becker, G.; Ehemann, V.; Schwab, M.; Westermann, F. MYCN/MYC-mediated drug resistance mechanisms in neuroblastoma. Int. J. Clin. Pharmacol. Ther. 2010, 48, 489–491. [Google Scholar] [CrossRef]
  60. Boon, K.; Caron, H.N.; van Asperen, R.; Valentijn, L.; Hermus, M.C.; van Sluis, P.; Roobeek, I.; Weis, I.; Voûte, P.A.; Schwab, M.; et al. N-myc enhances the expression of a large set of genes functioning in ribosome biogenesis and protein synthesis. Embo J. 2001, 20, 1383–1393. [Google Scholar] [CrossRef] [PubMed]
  61. Hald, Ø.H.; Olsen, L.; Gallo-Oller, G.; Elfman, L.H.M.; Løkke, C.; Kogner, P.; Sveinbjörnsson, B.; Flægstad, T.; Johnsen, J.I.; Einvik, C. Inhibitors of ribosome biogenesis repress the growth of MYCN-amplified neuroblastoma. Oncogene 2019, 38, 2800–2813. [Google Scholar] [CrossRef]
  62. Tao, T.; Sondalle, S.B.; Shi, H.; Zhu, S.; Perez-Atayde, A.R.; Peng, J.; Baserga, S.J.; Look, A.T. The pre-rRNA processing factor DEF is rate limiting for the pathogenesis of MYCN-driven neuroblastoma. Oncogene 2017, 36, 3852–3867. [Google Scholar] [CrossRef] [PubMed]
  63. O’Brien, R.; Tran, S.L.; Maritz, M.F.; Liu, B.; Kong, C.F.; Purgato, S.; Yang, C.; Murray, J.; Russell, A.J.; Flemming, C.L.; et al. MYC-Driven Neuroblastomas Are Addicted to a Telomerase-Independent Function of Dyskerin. Cancer Res. 2016, 76, 3604–3617. [Google Scholar] [CrossRef]
  64. Nakaguro, M.; Kiyonari, S.; Kishida, S.; Cao, D.; Murakami-Tonami, Y.; Ichikawa, H.; Takeuchi, I.; Nakamura, S.; Kadomatsu, K. Nucleolar protein PES1 is a marker of neuroblastoma outcome and is associated with neuroblastoma differentiation. Cancer Sci. 2015, 106, 237–243. [Google Scholar] [CrossRef] [PubMed]
  65. Zhang, X.; Liu, C.; Cao, Y.; Liu, L.; Sun, F.; Hou, L. RRS1 knockdown inhibits the proliferation of neuroblastoma cell via PI3K/Akt/NF-κB pathway. Pediatr. Res. 2022, 1–11. [Google Scholar] [CrossRef] [PubMed]
  66. Hsu, C.L.; Chang, H.Y.; Chang, J.Y.; Hsu, W.M.; Huang, H.C.; Juan, H.F. Unveiling MYCN regulatory networks in neuroblastoma via integrative analysis of heterogeneous genomics data. Oncotarget 2016, 7, 36293–36310. [Google Scholar] [CrossRef]
  67. Cortes, C.L.; Veiga, S.R.; Almacellas, E.; Hernández-Losa, J.; Ferreres, J.C.; Kozma, S.C.; Ambrosio, S.; Thomas, G.; Tauler, A. Effect of low doses of actinomycin D on neuroblastoma cell lines. Mol. Cancer 2016, 15, 1. [Google Scholar] [CrossRef]
  68. Zage, P.E.; Nolo, R.; Fang, W.; Stewart, J.; Garcia-Manero, G.; Zweidler-McKay, P.A. Notch pathway activation induces neuroblastoma tumor cell growth arrest. Pediatr. Blood Cancer 2012, 58, 682–689. [Google Scholar] [CrossRef]
  69. Van Limpt, V.A.; Chan, A.J.; Van Sluis, P.G.; Caron, H.N.; Van Noesel, C.J.; Versteeg, R. High delta-like 1 expression in a subset of neuroblastoma cell lines corresponds to a differentiated chromaffin cell type. Int. J. Cancer 2003, 105, 61–69. [Google Scholar] [CrossRef]
  70. Shostak, A.; Ruppert, B.; Ha, N.; Bruns, P.; Toprak, U.H.; Eils, R.; Schlesner, M.; Diernfellner, A.; Brunner, M. MYC/MIZ1-dependent gene repression inversely coordinates the circadian clock with cell cycle and proliferation. Nat. Commun. 2016, 7, 11807. [Google Scholar] [CrossRef]
  71. Iraci, N.; Diolaiti, D.; Papa, A.; Porro, A.; Valli, E.; Gherardi, S.; Herold, S.; Eilers, M.; Bernardoni, R.; Della Valle, G.; et al. A SP1/MIZ1/MYCN repression complex recruits HDAC1 at the TRKA and p75NTR promoters and affects neuroblastoma malignancy by inhibiting the cell response to NGF. Cancer Res. 2011, 71, 404–412. [Google Scholar] [CrossRef]
  72. Benetatos, L.; Vartholomatos, G.; Hatzimichael, E. Polycomb group proteins and MYC: The cancer connection. Cell. Mol. Life Sci. 2014, 71, 257–269. [Google Scholar] [CrossRef] [PubMed]
  73. Corvetta, D.; Chayka, O.; Gherardi, S.; D’Acunto, C.W.; Cantilena, S.; Valli, E.; Piotrowska, I.; Perini, G.; Sala, A. Physical interaction between MYCN oncogene and polycomb repressive complex 2 (PRC2) in neuroblastoma: Functional and therapeutic implications. J. Biol. Chem. 2013, 288, 8332–8341. [Google Scholar] [CrossRef]
  74. Zhao, X.; Lwin, T.; Zhang, X.; Huang, A.; Wang, J.; Marquez, V.E.; Chen-Kiang, S.; Dalton, W.S.; Sotomayor, E.; Tao, J. Disruption of the MYC-miRNA-EZH2 loop to suppress aggressive B-cell lymphoma survival and clonogenicity. Leukemia 2013, 27, 2341–2350. [Google Scholar] [CrossRef]
  75. Sun, Y.; Liu, P.Y.; Scarlett, C.J.; Malyukova, A.; Liu, B.; Marshall, G.M.; MacKenzie, K.L.; Biankin, A.V.; Liu, T. Histone deacetylase 5 blocks neuroblastoma cell differentiation by interacting with N-Myc. Oncogene 2014, 33, 2987–2994. [Google Scholar] [CrossRef]
  76. Brenner, C.; Deplus, R.; Didelot, C.; Loriot, A.; Viré, E.; De Smet, C.; Gutierrez, A.; Danovi, D.; Bernard, D.; Boon, T.; et al. Myc represses transcription through recruitment of DNA methyltransferase corepressor. Embo J. 2005, 24, 336–346. [Google Scholar] [CrossRef]
  77. Pang, Y.; Liu, J.; Li, X.; Xiao, G.; Wang, H.; Yang, G.; Li, Y.; Tang, S.C.; Qin, S.; Du, N.; et al. MYC and DNMT3A-mediated DNA methylation represses microRNA-200b in triple negative breast cancer. J. Cell Mol. Med. 2018, 22, 6262–6274. [Google Scholar] [CrossRef] [PubMed]
  78. Vo, B.T.; Wolf, E.; Kawauchi, D.; Gebhardt, A.; Rehg, J.E.; Finkelstein, D.; Walz, S.; Murphy, B.L.; Youn, Y.H.; Han, Y.G.; et al. The Interaction of Myc with Miz1 Defines Medulloblastoma Subgroup Identity. Cancer Cell. 2016, 29, 5–16. [Google Scholar] [CrossRef]
  79. Sogut, M.S.; Venugopal, C.; Kandemir, B.; Dag, U.; Mahendram, S.; Singh, S.; Gulfidan, G.; Arga, K.Y.; Yilmaz, B.; Kurnaz, I.A. ETS-Domain Transcription Factor Elk-1 Regulates Stemness Genes in Brain Tumors and CD133+ BrainTumor-Initiating Cells. J. Pers. Med. 2021, 11, 125. [Google Scholar] [CrossRef]
  80. Güner, G.; Lichtenthaler, S.F. The substrate repertoire of γ-secretase/presenilin. Semin. Cell Dev. Biol. 2020, 105, 27–42. [Google Scholar] [CrossRef] [PubMed]
  81. Hitzenberger, M.; Götz, A.; Menig, S.; Brunschweiger, B.; Zacharias, M.; Scharnagl, C. The dynamics of γ-secretase and its substrates. Semin. Cell Dev. Biol. 2020, 105, 86–101. [Google Scholar] [CrossRef]
  82. McCaw, T.R.; Inga, E.; Chen, H.; Jaskula-Sztul, R.; Dudeja, V.; Bibb, J.A.; Ren, B.; Rose, J.B. Gamma Secretase Inhibitors in Cancer: A Current Perspective on Clinical Performance. Oncologist 2021, 26, e608–e621. [Google Scholar] [CrossRef] [PubMed]
  83. López-Nieva, P.; González-Sánchez, L.; Cobos-Fernández, M.; Córdoba, R.; Santos, J.; Fernández-Piqueras, J. More Insights on the Use of γ-Secretase Inhibitors in Cancer Treatment. Oncologist 2021, 26, e298–e305. [Google Scholar] [CrossRef] [PubMed]
  84. Ferrari-Toninelli, G.; Bonini, S.A.; Uberti, D.; Buizza, L.; Bettinsoli, P.; Poliani, P.L.; Facchetti, F.; Memo, M. Targeting Notch pathway induces growth inhibition and differentiation of neuroblastoma cells. Neuro Oncol. 2010, 12, 1231–1243. [Google Scholar] [CrossRef] [PubMed]
  85. Grynfeld, A.; Påhlman, S.; Axelson, H. Induced neuroblastoma cell differentiation, associated with transient HES-1 activity and reduced HASH-1 expression, is inhibited by Notch1. Int. J. Cancer 2000, 88, 401–410. [Google Scholar] [CrossRef]
  86. Chang, H.H.; Lee, H.; Hu, M.K.; Tsao, P.N.; Juan, H.F.; Huang, M.C.; Shih, Y.Y.; Wang, B.J.; Jeng, Y.M.; Chang, C.L.; et al. Notch1 expression predicts an unfavorable prognosis and serves as a therapeutic target of patients with neuroblastoma. Clin. Cancer Res. 2010, 16, 4411–4420. [Google Scholar] [CrossRef]
  87. Allen, T.D.; Rodriguez, E.M.; Jones, K.D.; Bishop, J.M. Activated Notch1 induces lung adenomas in mice and cooperates with Myc in the generation of lung adenocarcinoma. Cancer Res. 2011, 71, 6010–6018. [Google Scholar] [CrossRef]
  88. Herranz, D.; Ambesi-Impiombato, A.; Palomero, T.; Schnell, S.A.; Belver, L.; Wendorff, A.A.; Xu, L.; Castillo-Martin, M.; Llobet-Navás, D.; Cordon-Cardo, C.; et al. A NOTCH1-driven MYC enhancer promotes T cell development, transformation and acute lymphoblastic leukemia. Nat. Med. 2014, 20, 1130–1137. [Google Scholar] [CrossRef]
  89. Ryan, R.J.H.; Petrovic, J.; Rausch, D.M.; Zhou, Y.; Lareau, C.A.; Kluk, M.J.; Christie, A.L.; Lee, W.Y.; Tarjan, D.R.; Guo, B.; et al. A B Cell Regulome Links Notch to Downstream Oncogenic Pathways in Small B Cell Lymphomas. Cell Rep. 2017, 21, 784–797. [Google Scholar] [CrossRef]
  90. van Groningen, T.; Koster, J.; Valentijn, L.J.; Zwijnenburg, D.A.; Akogul, N.; Hasselt, N.E.; Broekmans, M.; Haneveld, F.; Nowakowska, N.E.; Bras, J.; et al. Neuroblastoma is composed of two super-enhancer-associated differentiation states. Nat. Genet. 2017, 49, 1261–1266. [Google Scholar] [CrossRef]
  91. Boeva, V.; Louis-Brennetot, C.; Peltier, A.; Durand, S.; Pierre-Eugène, C.; Raynal, V.; Etchevers, H.C.; Thomas, S.; Lermine, A.; Daudigeos-Dubus, E.; et al. Heterogeneity of neuroblastoma cell identity defined by transcriptional circuitries. Nat. Genet. 2017, 49, 1408–1413. [Google Scholar] [CrossRef] [PubMed]
  92. van Groningen, T.; Akogul, N.; Westerhout, E.M.; Chan, A.; Hasselt, N.E.; Zwijnenburg, D.A.; Broekmans, M.; Stroeken, P.; Haneveld, F.; Hooijer, G.K.J.; et al. A NOTCH feed-forward loop drives reprogramming from adrenergic to mesenchymal state in neuroblastoma. Nat. Commun. 2019, 10, 1530. [Google Scholar] [CrossRef] [PubMed]
  93. Zimmerman, M.W.; Liu, Y.; He, S.; Durbin, A.D.; Abraham, B.J.; Easton, J.; Shao, Y.; Xu, B.; Zhu, S.; Zhang, X.; et al. MYC Drives a Subset of High-Risk Pediatric Neuroblastomas and Is Activated through Mechanisms Including Enhancer Hijacking and Focal Enhancer Amplification. Cancer Discov. 2018, 8, 320–335. [Google Scholar] [CrossRef] [PubMed]
  94. Westermann, F.; Muth, D.; Benner, A.; Bauer, T.; Henrich, K.O.; Oberthuer, A.; Brors, B.; Beissbarth, T.; Vandesompele, J.; Pattyn, F.; et al. Distinct transcriptional MYCN/c-MYC activities are associated with spontaneous regression or malignant progression in neuroblastomas. Genome Biol. 2008, 9, R150. [Google Scholar] [CrossRef]
  95. Williamson, D.; Schwalbe, E.C.; Hicks, D.; Aldinger, K.A.; Lindsey, J.C.; Crosier, S.; Richardson, S.; Goddard, J.; Hill, R.M.; Castle, J.; et al. Medulloblastoma group 3 and 4 tumors comprise a clinically and biologically significant expression continuum reflecting human cerebellar development. Cell Rep. 2022, 40, 111162. [Google Scholar] [CrossRef] [PubMed]
  96. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  97. Heinz, S.; Benner, C.; Spann, N.; Bertolino, E.; Lin, Y.C.; Laslo, P.; Cheng, J.X.; Murre, C.; Singh, H.; Glass, C.K. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 2010, 38, 576–589. [Google Scholar] [CrossRef]
  98. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef]
  99. Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W.; et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9, R137. [Google Scholar] [CrossRef]
  100. McLean, C.Y.; Bristor, D.; Hiller, M.; Clarke, S.L.; Schaar, B.T.; Lowe, C.B.; Wenger, A.M.; Bejerano, G. GREAT improves functional interpretation of cis-regulatory regions. Nat. Biotechnol. 2010, 28, 495–501. [Google Scholar] [CrossRef]
  101. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’ayan, A. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [PubMed]
  102. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [PubMed]
  103. Xie, Z.; Bailey, A.; Kuleshov, M.V.; Clarke, D.J.B.; Evangelista, J.E.; Jenkins, S.L.; Lachmann, A.; Wojciechowicz, M.L.; Kropiwnicki, E.; Jagodnik, K.M.; et al. Gene Set Knowledge Discovery with Enrichr. Curr. Protoc. 2021, 1, e90. [Google Scholar] [CrossRef] [PubMed]
  104. Lutz, W.; Stöhr, M.; Schürmann, J.; Wenzel, A.; Löhr, A.; Schwab, M. Conditional expression of N-myc in human neuroblastoma cells increases expression of alpha-prothymosin and ornithine decarboxylase and accelerates progression into S-phase early after mitogenic stimulation of quiescent cells. Oncogene 1996, 13, 803–812. [Google Scholar] [PubMed]
Ijms 24 08141 g001 550
Figure 1. MYCN-regulated gene ontology terms in NB. Top 10 biological processes based on analysis of up-regulated (A) and down-regulated (B) genes in human NB MYCN-amplified cell lines when compared to non-MYCN-amplified cell lines. (C) Overlap between common MYCN target genes identified in three human MYCN-amplified cell lines and differentially regulated genes in MYCN-amplified human NB cell lines. Enrichment of MYCN ChIP-Seq peaks at the TSSs of up-regulated (D) and down-regulated (E) genes in the MYCN-amplified KELLY cell line. (F) Analysis of transcription factor binding motif. (Top), MYCN-bound and up-regulated genes. (Bottom), MYCN-bound and down-regulated genes. A false discovery rate (FDR) ≤ 0.05 was used as the cut-off value to define the statistical significance.
Figure 1. MYCN-regulated gene ontology terms in NB. Top 10 biological processes based on analysis of up-regulated (A) and down-regulated (B) genes in human NB MYCN-amplified cell lines when compared to non-MYCN-amplified cell lines. (C) Overlap between common MYCN target genes identified in three human MYCN-amplified cell lines and differentially regulated genes in MYCN-amplified human NB cell lines. Enrichment of MYCN ChIP-Seq peaks at the TSSs of up-regulated (D) and down-regulated (E) genes in the MYCN-amplified KELLY cell line. (F) Analysis of transcription factor binding motif. (Top), MYCN-bound and up-regulated genes. (Bottom), MYCN-bound and down-regulated genes. A false discovery rate (FDR) ≤ 0.05 was used as the cut-off value to define the statistical significance.
Ijms 24 08141 g001
Ijms 24 08141 g002 550
Figure 2. MYCN is a potential direct regulator of genes involved in differentially regulated biological processes. (A) Top 10 gene ontologies enriched with MYCN-bound genes. (Left), up-regulated genes. (Right), down-regulated genes. (B) R2 genome browser analysis of MYCN ChIP-Seq peaks at TSSs of the RPL5 gene, representative of MYCN-bound and up-regulated genes involved in ribosomal biogenesis (Top panel), and the APH1A gene, representative of MYCN-bound and down-regulated genes involved in NOTCH signaling pathway (Bottom panel). Arrowhead indicates the direction of transcription. MYCN ChIP-Seq pileup at the TSSs of genes involved in (C) ribosomal biogenesis and (D) NOTCH signaling pathway. (E) Pathway analysis of differentially up-regulated (left) and down-regulated (right) genes in IMR32 cell line following MYCN knockdown. False discovery rate (FDR) ≤ 0.05 was used as cut-off value to define statistical significance.
Figure 2. MYCN is a potential direct regulator of genes involved in differentially regulated biological processes. (A) Top 10 gene ontologies enriched with MYCN-bound genes. (Left), up-regulated genes. (Right), down-regulated genes. (B) R2 genome browser analysis of MYCN ChIP-Seq peaks at TSSs of the RPL5 gene, representative of MYCN-bound and up-regulated genes involved in ribosomal biogenesis (Top panel), and the APH1A gene, representative of MYCN-bound and down-regulated genes involved in NOTCH signaling pathway (Bottom panel). Arrowhead indicates the direction of transcription. MYCN ChIP-Seq pileup at the TSSs of genes involved in (C) ribosomal biogenesis and (D) NOTCH signaling pathway. (E) Pathway analysis of differentially up-regulated (left) and down-regulated (right) genes in IMR32 cell line following MYCN knockdown. False discovery rate (FDR) ≤ 0.05 was used as cut-off value to define statistical significance.
Ijms 24 08141 g002
Ijms 24 08141 g003 550
Figure 3. MYCN amplification correlates with reduced expression of genes encoding γ-secretase subunits and ADAM17 in human NB cell lines and primary tumors. (A) Heatmap analysis of mRNA levels of the γ-secretase complex genes and ADAM17 in the 39-human NB-cell-line dataset. (B) RT-qPCR analysis of baseline expression levels of genes of the γ-secretase complex and ADAM17 in a set of human NB cell lines. Error bars in panel represent the standard deviation of three independent biological experiments. p-value was calculated using a two-tailed, unpaired student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; *** p-value ≤ 0.001. (C) In silico analysis of mRNA levels of γ-secretase complex genes and ADAM17 in primary NB tumors with and without MYCN amplification. Statistical analysis was performed using one-way ANOVA. (D) Kaplan–Meier-analysis overall survival curve from the Tumor Neuroblastoma-SEQC-498 study based on mRNA expressions of genes of the γ-secretase complex, APH1B, NCSTN, and PSEN1. Statistical analysis was performed using Bonferroni correction of raw p-values.
Figure 3. MYCN amplification correlates with reduced expression of genes encoding γ-secretase subunits and ADAM17 in human NB cell lines and primary tumors. (A) Heatmap analysis of mRNA levels of the γ-secretase complex genes and ADAM17 in the 39-human NB-cell-line dataset. (B) RT-qPCR analysis of baseline expression levels of genes of the γ-secretase complex and ADAM17 in a set of human NB cell lines. Error bars in panel represent the standard deviation of three independent biological experiments. p-value was calculated using a two-tailed, unpaired student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; *** p-value ≤ 0.001. (C) In silico analysis of mRNA levels of γ-secretase complex genes and ADAM17 in primary NB tumors with and without MYCN amplification. Statistical analysis was performed using one-way ANOVA. (D) Kaplan–Meier-analysis overall survival curve from the Tumor Neuroblastoma-SEQC-498 study based on mRNA expressions of genes of the γ-secretase complex, APH1B, NCSTN, and PSEN1. Statistical analysis was performed using Bonferroni correction of raw p-values.
Ijms 24 08141 g003
Ijms 24 08141 g004 550
Figure 4. MYCN depletion induces expression of genes encoding γ-secretase members and ADAM17. (A) Western blot analysis of MYCN protein levels following treatment of TET21N cell line with 1 µg/mL doxycycline (Dox) for 72 h. (B) RT-qPCR analysis of mRNA levels of genes encoding γ-secretase complex and ADAM17 in response to MYCN Dox-induced depletion in the TET21N cell line. (C,D) Detection of endogenous MYCN:MAX interactions by isPLA in the KELLY cell line. Treatment with 5 µM of MYCMI-7 for 5 h disrupts MYCN:MAX interactions. (C) Representative images of MYCN:MAX isPLA signals in the KELLY cell line. Scale bar is 20 µM. (D) Quantification of MYCN:MAX isPLA signals in the KELLY cell line. isPLA signals using MYCN or MAX antibodies alone were used as background levels. DAPI was used to stain the nucleus (Blue), Phalloidin (AF488, Green) was used to stain the cytoplasm, and isPLA signals (Atto 647N, Red) detected MYCN-MAX interactions. (E) Western blot analysis of MYCN protein levels in the KELLY cell line following 72 h treatment with 0.4 µM MYCMI-7. (F) RT-qPCR analysis of mRNA levels of genes encoding the γ-secretase complex and ADAM17 in response to MYCN pharmacological depletion in the KELLY cell line. Error bars represent the standard deviation of three independent biological experiments. For the qPCR, the p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; *** p-value ≤ 0.001. For the isPLA, the statistical significance was calculated using ANOVA in R.
Figure 4. MYCN depletion induces expression of genes encoding γ-secretase members and ADAM17. (A) Western blot analysis of MYCN protein levels following treatment of TET21N cell line with 1 µg/mL doxycycline (Dox) for 72 h. (B) RT-qPCR analysis of mRNA levels of genes encoding γ-secretase complex and ADAM17 in response to MYCN Dox-induced depletion in the TET21N cell line. (C,D) Detection of endogenous MYCN:MAX interactions by isPLA in the KELLY cell line. Treatment with 5 µM of MYCMI-7 for 5 h disrupts MYCN:MAX interactions. (C) Representative images of MYCN:MAX isPLA signals in the KELLY cell line. Scale bar is 20 µM. (D) Quantification of MYCN:MAX isPLA signals in the KELLY cell line. isPLA signals using MYCN or MAX antibodies alone were used as background levels. DAPI was used to stain the nucleus (Blue), Phalloidin (AF488, Green) was used to stain the cytoplasm, and isPLA signals (Atto 647N, Red) detected MYCN-MAX interactions. (E) Western blot analysis of MYCN protein levels in the KELLY cell line following 72 h treatment with 0.4 µM MYCMI-7. (F) RT-qPCR analysis of mRNA levels of genes encoding the γ-secretase complex and ADAM17 in response to MYCN pharmacological depletion in the KELLY cell line. Error bars represent the standard deviation of three independent biological experiments. For the qPCR, the p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; *** p-value ≤ 0.001. For the isPLA, the statistical significance was calculated using ANOVA in R.
Ijms 24 08141 g004
Ijms 24 08141 g005 550
Figure 5. MYCN depletion induces expression of NOTCH target genes. (A) In silico expression analysis of HEY1 and HEYL NOTCH target genes in 39 human NB cell lines. (B) In silico analysis of mRNA levels of NOTCH target genes HES1, HES4, HEY1, and HEYL in primary NB tumors. (C) Kaplan–Meier-analysis overall survival curve from the Tumor Neuroblastoma-SEQC-498 study based on mRNA expression of NOTCH target genes HES1, HES4, HEY1, and HEYL. (D,E) Expression of HES1, HES4, HEY1, and HEYL in response to 1 µg/mL doxycycline in (D) TET21N cell line or (E) KELLY cell line following 0.4 µM MYCMI-7 treatment. Cell lines were treated with doxycycline or MYCMI-7 for 72 h. DMSO was used as control treatment. Statistical significance was calculated using (A) unpaired Student t-test in panel, and (B) one-way ANOVA in panel. (C) Statistical analysis was performed using Bonferroni correction of raw p-values available in panel. For qPCR, p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; *** p-value ≤ 0.001.
Figure 5. MYCN depletion induces expression of NOTCH target genes. (A) In silico expression analysis of HEY1 and HEYL NOTCH target genes in 39 human NB cell lines. (B) In silico analysis of mRNA levels of NOTCH target genes HES1, HES4, HEY1, and HEYL in primary NB tumors. (C) Kaplan–Meier-analysis overall survival curve from the Tumor Neuroblastoma-SEQC-498 study based on mRNA expression of NOTCH target genes HES1, HES4, HEY1, and HEYL. (D,E) Expression of HES1, HES4, HEY1, and HEYL in response to 1 µg/mL doxycycline in (D) TET21N cell line or (E) KELLY cell line following 0.4 µM MYCMI-7 treatment. Cell lines were treated with doxycycline or MYCMI-7 for 72 h. DMSO was used as control treatment. Statistical significance was calculated using (A) unpaired Student t-test in panel, and (B) one-way ANOVA in panel. (C) Statistical analysis was performed using Bonferroni correction of raw p-values available in panel. For qPCR, p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; *** p-value ≤ 0.001.
Ijms 24 08141 g005
Ijms 24 08141 g006 550
Figure 6. Inhibition of γ-secretase complex counteracts the induced expression of NOTCH target genes upon MYCN depletion. TET21N cell line was treated with 1 µg/mL of doxycycline, 5 µM of DAPT, or a combination for 72 h. DMSO treatment was used as control. Error bars represent the standard deviation of three independent biological experiments. p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; ns: no significance.
Figure 6. Inhibition of γ-secretase complex counteracts the induced expression of NOTCH target genes upon MYCN depletion. TET21N cell line was treated with 1 µg/mL of doxycycline, 5 µM of DAPT, or a combination for 72 h. DMSO treatment was used as control. Error bars represent the standard deviation of three independent biological experiments. p-value was calculated using two-tailed, unpaired Student t-test. * p-value ≤ 0.05; ** p-value ≤ 0.01; ns: no significance.
Ijms 24 08141 g006
Table
Table 1. Description of the NOTCH receptor processing enzymes identified in this study.
Table 1. Description of the NOTCH receptor processing enzymes identified in this study.
Gene SymbolGene NameFunction
APH1AAph-1 Homolog APromotes nicastrin and presenilin association in the γ-secretase complex, and mediates γ-secretase interaction with substrates prior to cleavage
APH1BAph-1 Homolog BPresenilin-stabilizing cofactor
NCSTNNicastrinCleaves integral membrane proteins, including NOTCH receptors. It is believed to function as a stabilizing cofactor for gamma-secretase complex assembly
PSEN1Presenilin 1Protease, regulates gamma-secretase activity, cleaves the intracellular domain of NOTCH receptors leading to the production of active NOTCH intracellular domain (NICD)
PSENENPresenilin enhancerProtease, regulates gamma-secretase activity, cleaves the intracellular domain of NOTCH receptors, leading to the production of active NICD
ADAM17 Disintegrin and metalloprotease domain 17Membrane-anchored protease, mediates the initial cleavage of the extracellular domain of NOTCH receptors
All information related to the functions of these proteins were retrieved from the NCBI.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Agarwal, P.; Glowacka, A.; Mahmoud, L.; Bazzar, W.; Larsson, L.-G.; Alzrigat, M. MYCN Amplification Is Associated with Reduced Expression of Genes Encoding γ-Secretase Complex and NOTCH Signaling Components in Neuroblastoma. Int. J. Mol. Sci. 2023, 24, 8141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24098141

AMA Style

Agarwal P, Glowacka A, Mahmoud L, Bazzar W, Larsson L-G, Alzrigat M. MYCN Amplification Is Associated with Reduced Expression of Genes Encoding γ-Secretase Complex and NOTCH Signaling Components in Neuroblastoma. International Journal of Molecular Sciences. 2023; 24(9):8141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24098141

Chicago/Turabian Style

Agarwal, Prasoon, Aleksandra Glowacka, Loay Mahmoud, Wesam Bazzar, Lars-Gunnar Larsson, and Mohammad Alzrigat. 2023. "MYCN Amplification Is Associated with Reduced Expression of Genes Encoding γ-Secretase Complex and NOTCH Signaling Components in Neuroblastoma" International Journal of Molecular Sciences 24, no. 9: 8141. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24098141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop