MYC-family protein overexpression and prominent nucleolar formation represent prognostic indicators and potential therapeutic targets for aggressive high-MKI neuroblastomas: a report from the children’s oncology group

INTRODUCTION

Peripheral neuroblastic tumors (pNTs including neuroblastoma, ganglioneuroblastoma, and ganglioneuroma) represent one of the most common extra-cranial neoplasms in children [1]. They are composed of biologically distinct groups of tumors whose genomic/molecular characteristics can directly be associated with a wide range of clinical behaviors in individual patients from a benign or non-aggressive course largely due to spontaneous regression or tumor maturation to an extremely aggressive course often resulting in fatal outcome. International efforts have focused on development of the most appropriate protocols for managing pNT patients by stratifying them according to the defined risk grouping systems based on the combination of various prognostic factors [2]. For example, the Children’s Oncology Group (COG) Neuroblastoma Classification Biology Study, ANBL00B1, has been utilizing Age at Diagnosis [3], Clinical Staging (International Neuroblastoma Staging System: INSS) [4], Histopathology (International Neuroblastoma Pathology Classification: INPC) [5–8], MYCN oncogene amplification status [9, 10], DNA index [11], and Chromosomal aberrations [12] for assigning Low, Intermediate, or High Risk protocol to pNT patients from North America, Australia, and New Zealand [13].

The INPC, one of the prognostic factors critically contributing to the COG risk grouping, was established in 1999 [5, 6] and modified in 2003 [8]. It includes four morphologic indicators: Schwannian stromal development (distinguishing Neuroblastoma category, Ganglioneuroblastoma, Intermixed category, and Ganglioneuroma category), Grade of neuroblastic differentiation [distinguishing Undifferentiated (NB-U), Poorly differentiated (NB-PD), and Differentiating (NB-D) subtypes in Neuroblastoma category], Mitosis-Karyorrhexis index [MKI: distinguishing Low (< 100/5,000 cells), Intermediate (100–200/5,000 cells) and High (> 200/5,000 cells) class in Neuroblastoma category], and grossly visible neuroblastomatous nodule formation co-existing with Ganglioneuroblastoma, Intermixed or Ganglioneuroma (identifying Ganglioneuroblastoma, Nodular category). Among these morphologic indicators, there is a reproducible relationship between MYCN amplification and increased MKI classes [14, 15]: Our recent study confirmed that 68% of High-MKI, 15% of Intermediate-MKI, and only 3% of Low-MKI neuroblastomas had MYCN oncogene amplification [16].

MYCN oncogene amplification has been recognized as the strongest indicator of aggressive tumor behavior for pNT patients. Recently we reported that MYCN protein overexpression after DNA amplification could be the critical step in bringing the aggressive phenotype to pNTs [17]. Our latest study further showed that MYCN protein overexpression was detected in ~20% of patients with neuroblastoma of undifferentiated (NB-U) and poorly differentiated (NB-PD) subtypes [18]. We also discovered that MYC (aka C-myc) protein overexpression was observed in another ~10% of patients in the same cohort. Because these two subsets of neuroblastomas showed almost identical survival curves, we have combined them together and initially termed them MYC-driven neuroblastoma [18]. However, to avoid the confusion among physicians and researchers in the field, we have renamed this highly aggressive subset of neuroblastomas as “MYC family protein-driven neuroblastomas” (or MYC family-driven neuroblastomas in short). While investigating these MYC family-driven neuroblastomas, we recognized that their tumor cells were often associated with the presence of prominent nucleoli.

MYC family genes are among the most frequently activated genes in human cancer [19] and their activation leads to the elevated expression of their protein products. MYC family proteins include MYC, MYCN and MYCL, and they are known to play pivotal roles in development, differentiation, cell growth, cell death, stem cell self-renewal and probably more. MYC family proteins are basic-helix-loop-helix-leucine zipper transcription factors that activate the expression of MYC responsive genes. To activate their target genes, MYC family proteins first bind to their obligatory partner protein called MAX and form heterodimers. These active forms of MYC family proteins then bind to the specific DNA sequence, called E-Boxes. Notably, there are numerous canonical and non-canonical E-box sequences with varying affinities in the genome, and thus cells with high levels of MYC family protein expression tend to occupy more E-boxes than those with low MYC family protein expression. Upon the specific DNA binding, MYC family proteins recruit the scaffold protein, TRRAP, to the sites [19]. This molecular complex further recruits histone acetyl transferases such as GCN5. Therefore, high MYC family protein expression subsequently results in a genome-wide increase in histone acetylation and ultimately high contents of euchromatin [20].

MYC family proteins are known to activate many genes involved in rRNA synthesis, protein translation and mitochondriogenesis, glucose and glutamine metabolism, and lipid synthesis [21]. Therefore cells with high levels of MYC family protein expression tend to show hypertorophy of the cell body [22, 23], nuclei [20] and nucleoli [24]. These cellular changes brought by the augmented MYC family protein expression ultimately contribute to malignant transformation of the cells.

In this study, we focused on a subset of unfavorable neuroblastomas with a High-MKI and analyzed the pattern of MYC and MYCN protein expression compared to the distribution of conventional prognostic factors based on their MYCN oncogene status (amplified vs. non-amplified). We further examined a possible relationship between prominent nucleolar formation and augmented expression of MYC family proteins. Then we performed survival analyses to test prognostic significance of MYC family protein expression, prominent nucleolar formation, and other conventional prognostic factors. We also analyzed their prognostic relationships in this cohort of High-MKI neuroblastoma patients. Lastly, we tested the feasibility of small molecule inhibitors that specifically target hyperactive rRNA synthesis and augmented protein translation suppressing growth and elevated expression of MYC family proteins in unfavorable neuroblastoma cells.

RESULTS

In this study, we analyzed a total of 241 High-MKI cases including 120 MYCN amplified (NB-U 28, NB-PD 92) and 121 non-MYCN amplified (NB-U 4, NB-PD 117) tumors (Figure 1A and 1B), all of which had available H&E slides and unstained slides in the file of COG Neuroblastoma Pathology Reference Laboratory (see details in Materials and Methods).

Pathology review

For those 241 cases, H&E slides were re-reviewed by RN-T, RM, LW, and HS, according to the International Neuroblastoma Pathology Classification (INPC). The median number of slides reviewed was 3.80 (range 1 to 26) per case. Diagnostic criteria of INPC were previously described [5, 6]. High-MKI tumors were defined when the number of mitotic and karyorrhectic cells per 5,000 denominator cells was more than 200 [5]. Definition of prominent nucleolar formation is described in our previous reports [17, 25]. In the prominent nucleoli (+) tumors, 10% or more of neuroblastic cells had one or a few large prominent nucleoli (Figure 1C). In the prominent nucleoli (−) tumors, a majority of neuroblastic cells had conventional salt-and-pepper pattern of nuclei (none or less than 10% of neuroblastic cells with prominent nucleoli, Figure 1D).

The results of immunostaining were classified as follows (see Figure 1): negative (−); focal or sporadic, and weak nuclear staining (+/−); and diffuse and strong nuclear staining with typically heterogeneous intensity (+). However, in our previous study, the survival effect for negative (−) and focal or sporadic, and weak (±) staining were similar, therefore we combined both groups together for further analyses. Accordingly each tumor was determined to have either MYCN protein overexpression [MYCN protein (+)] or MYCN protein negative/weak expression [MYCN protein (−)/(+/−)] (Figures 1E, 1F, and 1G) and either MYC protein overexpression [MYC protein (+)] or MYC protein negative/weak expression [MYC protein (−)/(+/−)] (Figures 1H, 1I, and 1J). Appropriate positive and negative controls were stained along with those tumors.

Clinicopathological characteristics of High MKI Neuroblastomas based on MYCN oncogene status

Among MYCN amplified High-MKI cases (n = 120), 101 (84.2%) tumors were MYCN protein (+), one (0.8%) tumor was MYC protein (+), 2 tumors (1.7%) were both proteins (+), and 16 tumors (13.3%) were both proteins (−)/(+/−). In contrast, non-MYCN amplified High-MKI cases (n = 121) included 7 (5.8%) MYCN protein (+) tumors, 36 (29.8%) MYC protein (+) tumors, 3 (2.5%) both proteins (+) tumors, and 75 (62.0%) both proteins (−)/(+/−) tumors (see Figure 2A).

Controlling for MYCN amplification (see Table 1A, 1B), there were no significant differences in the distributions of conventional prognostic factors (Age at diagnosis, DNA Ploidy pattern, INSS clinical stage, and Grade of neuroblastic differentiation) between MYCN protein (+) vs. both proteins (−)/(+/−), MYC protein (+) vs. both proteins (−)/(+/−), and either protein (+) vs. both proteins (−)/(+/−).

There were, however, significant differences in the distribution of prominent nucleolar formations between all the three comparison groups within the non-MYCN amplified tumors. In each case, tumors with both proteins (−)/(+/−) were less likely to have prominent nucleolar formations. Within the MYCN amplified tumors, there were significant differences in distribution of prominent nucleolar formations between MYCN protein (+) vs. both (−)/(+/−), and either protein (+) vs. both (−)/(+/−), but not between MYC protein (+) vs. both (−)/(+/−) [Note: there were only 3 tumors with MYCN amplification and MYC protein (+)]. In both cases, MYC family-driven neuroblastomas with MYCN protein (+) or either protein (+) were more likely to have prominent nucleoli.

When the MYCN amplified and non-MYCN amplified High-MKI cases were combined (n = 241), 161 tumors were prominent nucleolar formation (+) and 139 (86.3%) of them were MYC family-driven neuroblastomas with MYCN protein (+), MYC protein (+), or both (+), as shown in Figure 2B. Eighty tumors were prominent nucleolar formation (−) and the majority of them [69 (86.2%)] did not have MYCN protein (+) and MYC protein (+). In summary, there was a significant relationship between MYCN/MYC protein overexpression and prominent nucleolar formation in this group of tumors (p < 0.0001).

Survival analysis

The median follow-up time for patients without event was 2.7 years. The median follow-up time for patients who did not die was 2.8 years. The 4-year event-free survival (EFS) and overall survival (OS) estimates for all 241 patients were 56.2+/−5.5% and 64.5+/−5.4%, respectively. As shown in Table 2A (all patients) and Figure 3, univariate analysis showed significant differences in survival between MYCN protein (+) and MYCN protein (−)/(+/−) patients, with more favorable outcomes for MYCN protein (−)/(+/−) patients. Significant differences in survival were also detected between prominent nuclear formation (+) and (−), with more favorable outcomes for prominent nucleolar formation (−) patients.

Survival of patients with MYC family-driven neuroblastomas where the tumors had MYCN protein (+) and/or MYC protein (+) was significantly worse than that of patients with non-MYC family-driven neuroblastomas where both proteins were (−)/(+/−). Four-year EFS for MYC family-driven High-MKI neuroblastomas was 50.2+/−6.2% and that for non-MYC family-driven High-MKI neuroblastomas was 66.3+/−11.1% (p = 0.0131; Table 2A). Four-year OS for MYC-family High-MKI neuroblastomas was 56.1+/−6.1% and that for non-MYC family-driven protein driven High-MKI neuroblastomas was 80.1+/−9.5% (p = 0.0064).

No other survival comparisons yielded significant results, including MYCN amplified vs. non-MYCN amplified patients and controlling for MYCN status (Table 2B). Some statistical tests were hampered by small sample sizes or, for the outcome comparisons, a lack of events. Such groups include: non-MYCN amplified and MYCN protein (+), MYCN amplified and MYCN protein (−)/ (+/−), and MYCN amplified and MYC protein (+).

Prognostic relationship among prominent nucleolar formation, MYC family protein expression and other conventional prognostic indicators

A full backward-selected Cox model in terms of EFS and OS including prominent nucleolar formation was fit, allowing any covariate to drop out of the model. The final Cox model indicated that prominent nucleolar formation [(+) vs. (−)] and INSS stage (stage 4 vs. non-stage 4) were both predictive of EFS and OS (Table 3A). A full backward-selected Cox model in terms of EFS and OS including MYCN protein was fit, allowing any covariate to drop out of the model. The final Cox model indicated that MYCN protein [(+) vs. (−)/(+/−)] and INSS stage (stage 4 vs. non-stage 4) were predictive of both EFS and OS (Table 3B). A full backward-selected Cox model in terms of EFS and OS, including MYC protein, was fit allowing any covariate to drop out of the model. The final Cox model indicated that only INSS stage (stage 4 vs. non-stage 4) was predictive of survivals for EFS and OS (Table 3C).

Finally, a full backward-selected Cox model for EFS and OS including MYC family protein expression, i.e. both proteins (−)/(+/−) vs. either protein (+), was fit, allowing any covariate to drop out of the model. The final Cox model indicated that INSS stage and MYC family protein expression were both predictive of survival for EFS and OS (Table 3D). The order of removal was MYCN amplification, Age at diagnosis, DNA ploidy, and Grade of neuroblastic differentiation for EFS; Age at diagnosis, MYCN amplification, Grade of neuroblastic differentiation, and DNA ploidy for OS, with the least statistically significant term dropping out at each step. While controlling for INSS stage, the model predicted that patients who expressed MYC family proteins were at 1.870 times greater risk of an event than patients who expressed neither protein (p = 0.0087) and at 2.245 times greater risk of death than patients who expressed neither protein (p = 0.0045). In essence, MYC family protein expression retained its prognostic significance in the presence of standard neuroblastoma risk factors in the Cox models for both EFS and OS in High MKI neuroblastomas, and it had a better predictability of adverse outcome than MYCN protein expression alone.

Inhibition of overdriven protein synthesis by inhibitors of RNA Pol I and an aminoacyl-tRNA synthetase results in MYC/MYCN down regulation in MYC family-driven neuroblastoma cells

Our results show that prominent nucleolar formation (nucleolar hypertrophy), indicative of overdriven rRNA gene transcription as well as protein translation, could provide critical prognostic and risk-stratification information for neuroblastoma. At the same time, this observation could open up an opportunity for innovative therapy for aggressive neuroblastomas. To assess this possibility, we examined the effect of two potent and specific small molecule inhibitors of RNA Pol I activity and an aminoacyl tRNA synthetase (CX-5461 [26] and Halofuginone [27], respectively) on growth and MYC family protein expression in neuroblastoma cells.

As shown in Figure 4A, these small molecule inhibitors in fact inhibited growth of neuroblastoma cell lines with various genetic profiles at low to submicromolar concentrations in vitro in forty-eight hours. Moreover, the inhibitors destabilized MYC and MYCN proteins in neuroblastoma cells (Figure 4B). Interestingly, Halofuginone showed a rapid effect on the stability of MYC family proteins at 500 nM-1 μM (< 3 hours). In contrast, CX-5461 at 1 mM was less effective at the 3-hour time. At the concentration of 100 nM-250 nM for Halofuginone and 1 μM for CX-5461 for up to forty-eight hours of the drug treatments, MYC and MYCN expression were both down-regulated, though SKNAS was less responsive to CX-5461 than other cell lines. It should be mentioned that both compounds are among the most effective small molecules that can suppress the expression of MYC family proteins. There appeared to be no apparent correlation between the drug response and the genetic profiles of cell lines used in the experiment.

Inhibition of glutamyl-prolyl-tRNA synthetase by Halofuginone suppresses growth of orthotopic Kelly neuroblastoma xenografts in mice and MYCN expression

Next, we examined the preclinical efficacy of Halofuginone in orthotopic Kelly human neuroblastoma xenografts in mice using a silk film as the drug delivery device. Since Halofuginone has a better theoretical drug retention and release profile from the silk film than CX-5461 (see below), we tested Halofuginone as a single drug in this drug delivery model in this study. As shown in Figure 4C, Halofuginone showed a significant growth suppressive effect on the xenografts, as indicated by the time periods to reach certain tumor sizes being significantly longer in the Halofuginone treated group than those in the control group. Moreover, histological examination of the xenografts showed a marked suppression of MYCN protein expression in the Halofuginone treated xenograft at day 17 of the drug treatment. Since Halofuginone (20 μg) was released from the silk film over the initial 24-hour period (Figure 5), the down regulation of MYCN in response to Halofuginone appeared to be sustained for a long time period.

DISCUSSION

This report examined MYCN and MYC protein expression and prominent nucleolar formation in High-MKI neuroblastomas and assessed the contribution of these factors to survival of this unfavorable cohort of patients in relation to the conventional prognostic factors, including MYCN amplification. As expected, the vast majority (103/120, 85.8%) of MYCN amplified High-MKI neuroblastomas overly expressed MYCN protein, including 2 tumors also having simultaneous MYCN and MYC protein expression. In contrast, more than 60% of non-MYCN amplified High-MKI neuroblastomas were MYCN protein (−)/(+/−) and MYC protein (−)/(+/−). It is interesting to note that MYC protein overexpression was associated with around 30% of High MKI neuroblastomas without MYCN amplification. In our previous study, MYC protein overexpression was shown to be associated more with rather Low- or Intermediate-MKI than High-MKI in neuroblastomas [18]. However, this study demonstrated that MYC protein overexpression could also represent a considerable proportion of non-MYCN amplified-High-MKI tumors. It is also noted that there were some non-MYCN amplified-High-MKI tumors (N = 10, 8.3%) that overly expressed MYCN protein with or without simultaneous MYC protein expression. This observation suggests that certain mechanisms exist to allow overexpression or stabilization of MYCN protein without gene amplification, such as ones described in other reports [28, 29].

Tests comparing the distributions of conventional prognostic factors, such as age at diagnosis, DNA ploidy, INSS clinical stage, and grade of neuroblastic differentiation, showed no significant difference between MYCN amplified and non-MYCN amplified groups of High-MKI neuroblastomas. However, there were significant differences in the presence of prominent nucleolar formation in both MYCN amplified and non-amplified groups, and prominent nucleolar formation was more likely to be present in MYC family-driven neuroblastomas overly expressing MYCN and/or MYC protein. These results confirmed our previous findings using a different cohort of patients [18, 30]. The observation is also consistent with the facts that ribosomal RNA gene repeats organize into nucleoli [31] and that their transcription is regulated directly and indirectly by MYC family proteins [24, 32]. It is of interest to note that in this series of High-MKI neuroblastomas, there were a total of 22 tumors with prominent nucleolar formation which did not express excess amount of MYCN or MYC protein. Further investigations will be needed to assess the involvement of other proteins, such as MYCL [33–35], in association with prominent nucleolar formation.

According to the INPC, High MKI is a sign of poor prognosis regardless of age at diagnosis for neuroblastoma [5, 6, 16]. Notably however, MYCN amplification status was not predictive of survival in this cohort. The factors shown to be predictive of outcome in this cohort of High-MKI neuroblastomas were whether or not tumors expressed augmented levels of MYC family proteins or if they exhibited prominent nucleolar formation. It is also noted that patients with MYC family-driven High-MKI neuroblastomas had a significantly lower EFS and OS than patients with non-MYC family-driven High-MKI neuroblastomas. When we analyzed prognostic relationship among the conventional prognostic factors, MYC family protein expression and prominent nucleolar formation, INSS stage, MYCN protein expression, the combined MYC family protein expression and prominent nucleolar formation were the only factors that independently predicted clinical outcome of High-MKI neuroblastoma patients.

Nuclear morphology of neuroblastic cells is often described as a “salt-and-pepper” appearance. However, as shown in this study as well as in our previous studies, there are clinically aggressive MYC family-driven neuroblastomas whose tumor cells are often characterized by the presence of one or a few prominent nucleoli as a sign of nucleolar hypertrophy. Morphological changes of the nucleus and its substructures associated with tumor aggressiveness often reflect genetic, epigenetic and biological changes. In order to exhibit nucleolar hypertrophy, indicative of hyper-activated ribosomal DNA transcription, cells need to open chromatin at ribosomal RNA gene repeats. The ribosomal RNA gene repeats then have to be transcriptionally activated by the complex of necessary transcription factors and RNA Polymerase I (Pol I). MYC-family proteins, after forming the heterodimer with the partner MAX protein, play pivotal roles in the process of ribosomal DNA transcription: the MYC/MAX complex binds E-boxes near ribosomal RNA genes and recruits histone acetyl transferases, allowing open chromatin configuration [36]. The MYC/MAX complex can also up-regulate the expression of upstream binding transcription factor (UBTF) and Pol I-specific transcription initiation factor RRN3, which are obligatory transcription factors in ribosomal RNA gene expression [37]. Therefore, the emergence of prominent nucleoli in tumor cells most likely suggests elevated expression of MYC family proteins or their functional equivalents.

It has been suggested that selective inhibition of rRNA gene transcription in the nucleolus by small molecule inhibitors could be an effective therapeutic strategy for malignant tumors [37, 38]. On the other hand, augmented protein translation, which would link to hyperactive rRNA synthesis, could be suppressed by specific inhibition of various steps involved in translation. In this study, we examined the feasibility of such strategies in neuroblastoma. CX-5461 is a small molecule that disrupts the binding of the SL1 transcription factor to the rDNA promoter and prevents the initiation of rRNA synthesis by RNA Pol I [26]. Consequently, ribosomal assembly will be halted, leading to the accumulation of unassembled ribosomal proteins [37, 38]. Free ribosomal proteins will then promote cancer-specific activation of p53 [37, 38]. Interestingly, a similar set of free ribosomal proteins could also down-regulate MYC protein by distinct mechanisms [39, 40]. Accordingly, our study showed that CX-5461 down-regulated MYC and MYCN proteins and suppressed growth of neuroblastoma cells.

Considering the inhibitory activity on RNA Pol I transcription by CX-5461, one may be concerned with severe side effects of the drug. However, the preclinical evaluation of CX- 5461 by Drygin et al. showed that CX-5461 was well-tolerated at all tested schedules as judged by the absence of significant changes in animal body weights [26]. Thus, it is expected that side effects of CX-5461 would be low, if any. In addition, it has been reported that in a Phase I study (ACTRN12613001061729), CX-5461 is well tolerated with low-grade manageable adverse events to date. Inhibition of Pol I transcription in peripheral blood mononuclear cells indicated on-target drug activity in the initial dose cohorts. This trial is continuing in dose escalation [41]. Another Phase I/II Study of CX-5461 is currently on going (https://clinicaltrials.gov/ct2/show/NCT02719977?term=CX-5461#x0026;rank=1) and these trials will answer this question in the near future.

Halofuginone specifically inhibits glutamyl-prolyl tRNA synthetase (EPRS). This results in the accumulation of uncharged prolyl tRNAs, which then leads to suppression of protein translation via the induction of the amino acid starvation response [42]. In addition, inhibition of EPRS may also liberate one of the cofactors from the human multi-tRNA-synthetase complex (i.e., AIMP2). Freed AIMP2 could then migrate into the nucleus and inhibit FBP (FUBP1), which otherwise functions as the transcriptional activator of MYC. Consequently, Halofuginone could induce down regulation of MYC family proteins through both global translational suppression and transcriptional inhibition. Our study, for the first time, showed that Halofuginone effectively down-regulated MYC and MYCN proteins in MYC family driven neuroblastoma cells.

As about half of “prominent nucleolar formation (+)” or “MYC family protein-driven High-MKI” tumors fail the current super intensive therapy, inhibition of ribosomal RNA gene transcription and aminoacyl-tRNA synthetases by small molecule inhibitors [26, 27] could be a potential therapeutic approach for this subset of unfavorable neuroblastomas as our results suggest (Figure 4). To this end, both CX-5461 and Halofuginone are currently examined for their efficacy on various human diseases in clinical trials (ACTRN12613001061729; https://clinicaltrials.gov/ct2/show/NCT02719977; https://clinicaltrials.gov/ct2/show/NCT02525302?term=halofuginone&rank=3).

Finally, results of this study along with recently accumulated data clearly prompt us to develop a new scheme of pathology classification as well as risk stratification for patients with neuroblastoma by incorporating the information on the expression of MYC family proteins in conjunction with nuclear morphology of the tumor cells. Moreover, MYC-family protein overexpression and prominent nucleolar formation should be considered for potential therapeutic targets of aggressive unfavorable neuroblastomas.

MATERIALS AND METHODS

Patient cohort

Between January 1, 2009 and June 18, 2014 a total of 3,280 pNTs were reviewed at the COG Neuroblastoma Pathology Reference Laboratory, Department of Pathology and Laboratory Medicine, Children’s Hospital Los Angeles, Los Angeles, California. Of these, 2,605 tumors were in the Neuroblastoma (Schwannian stroma-poor) category and included the undifferentiated subtype (NB-U, 85 cases), the poorly differentiated subtype (NB-PD, 2,321 cases), and the differentiating subtype (NB-D, 199 cases). In order to avoid potential confusion in identifying/detecting unique prominent nucleoli, those NB-D cases were excluded from the study, as, by definition, differentiating neuroblasts in this subtype had one prominent nucleolus as a sign of nuclear maturation. As for the NB-U and NB-PD cases, MKI classes were determined in 2,366 tumors: 1,234 were Low MKI (NB-U 5, NB-PD 1045), 660 were Intermediate MKI (NB-U 19, NB-PB 628), and 472 were High MKI (NB-U 56, NB-PD 415). Among the High-MKI cases, MYCN status was determined by FISH test in 448 tumors, and results were reported as “Amplified (309 tumors; NB-U 51, NB-PD 258)” and “Non-amplified (139 tumors; NB-U 4, NB-PD 135)”. Informed consent approved by the institutional review board was obtained for all patients at the time of enrollment in the COG biology or therapeutic study.

Immunohistochemistry

Unstained sections from the study cases were heated for 30 minutes in Bond^TM Epitope Retrieval Solution 2 (No. AR9640; Leica Biosystems Newcastle Ltd, Benton Ln, Newcastle Upon Tyne, UK) using Leica BOND-MAX^TM (Leica Microsystems Inc, Bannockburn, IL, USA), and incubated with either anti-MYCN mouse monoclonal antibody (NCM II 100 [43] at a dilution of 1:200) or anti-human MYC rabbit monoclonal antibody (clone Y69; No. #1472–1; Epitomics, Cambridge, MA, USA at a dilution of 1:200) in Bond^TM Primary Antibody Diluent (No. AR9352; Vision BioSystems Inc., Norwell, MA, USA). The counter staining with hematoxylin was performed for MYC protein staining slides, but no counter staining was performed for MYCN protein staining slides.

Statistical analysis

First, we compared MYCN amplified and non-MYCN amplified High-MKI tumors based on their protein expression pattern. Then, controlling for MYCN amplification, the distribution of prognostic factors, such as presence of prominent nucleolar formation [PN (+) vs. PN (−)], Age at diagnosis (< 18 months vs. ≥ 18 months), INSS stage (non-stage 4 vs. stage 4), DNA ploidy pattern (hyperdiploid vs. diploid), and Grade of neuroblastic differentiation (NB-U vs NB-PD), were compared with chi-squared, or if small samples, Fisher’s exact tests, between the following groups: MYCN protein (+) vs. both proteins (−)/(+/−), MYC protein (+) vs. both proteins (−)/(+/−), and either protein (+) vs. both proteins (−)/(+/−).

Event-free survival (EFS) was calculated as the number of days from diagnosis to the occurrence of an event or, if no event, the date of last follow-up. An event is defined as death, disease relapse or progression, or secondary malignancy. Overall survival (OS) was calculated as the number of days from diagnosis to death or, if the patient did not die, the date of last follow-up. Kaplan-Meier EFS and OS estimates with standard errors per the methods of Peto et al. [44, 45] were computed and log-rank tests were performed to compare EFS and OS between the groups by MYCN oncogene status, prominent nucleolar formation, and MYCN/MYC protein expression status. Also analyzed were EFS and OS by MYCN/MYC protein expression status controlling for MYCN oncogene amplification status. Survival plots and life tables were produced. Survival probabilities are presented as 4-year estimates ± standard error, written as percentages. P -values less than 0.05 were considered statistically significant.

To determine the independent prognostic strength for survival of the protein expression factors in the presence of the prognostic factors; age at diagnosis, INSS stage, MYCN amplification status, DNA ploidy pattern, grade of neuroblastic differentiation, and prominent nucleolar formation, Cox proportional hazards (PH) models with the Efron method of handling tied event times were fit. Backward selection was used to determine the most parsimonious model.

Neuroblastoma cell growth assay and Western blot assay in vitro

MYCN amplified cell lines, IMR5 and Kelly, and non-MYCN amplified cell lines, SY5Y and SKNAS, were used in this study. IMR5 and SY5Y have normal TP53 gene, whereas Kelly and SKNAS cells harbor mutations in TP53 gene. IMR5 and Kelly cells express high levels of MYCN protein, whereas SY5Y and SKNAS cells express high levels of MYC protein. IMR5, a clone of IMR32, was a gift of Roger H. Kennett (University of Pennsylvania, PA. USA) who derived this clone. Kelly and SY5Y cell lines were obtained from Sigma (St. Louis, MO, USA). SKNAS cell line was a gift of Susan L. Cohn (University of Chicago, IL, USA). Growth of these neuroblastoma cell lines in response to CX-5461 and Halofuginone was assessed by MTS assay as described in our previous study [46]. Neuroblastoma cells were treated with the small molecule inhibitors at the concentrations and durations indicated in Figure 4. Cells were then subjected to Western blot assay using anti-pan-MYC antibody, NCM II 143 [43] and anti-β-Actin antibody (Santa Cruz) as described in our previous studies [47, 48].

Mouse orthotopic neuroblastoma model

All mouse procedures were performed in accordance with University of Illinois’ recommendations for the care and use of animals and were maintained and handled under protocols approved by the Institutional Animal Care and Use Committee. All procedures were performed with female NCr nude mice (Taconic, Hudson, NY, USA) at 7 weeks of age. Procedures and ultrasound measurements (see below) were performed under general anesthesia using isoflurane inhalation. Orthotopic tumors within the adrenal gland were created as described in our previous study [49]. Briefly, a transverse incision was made on the left flank to locate the left adrenal gland, and 2 mL of phosphate buffered saline (PBS) containing 10⁶ Kelly cells were injected into the adrenal gland using a 30G needle. The fascia and skin were closed in separate layers. Tumor formation was monitored by non-invasive ultrasound measurements, and the animals euthanized when the tumor volume exceeded 1,000 mm³.

Monitoring tumor growth by high frequency ultrasound

After securing the mouse in a prone position, a VisualSonics Vevo 2100 sonographic probe (Toronto, Ontario, Canada) was applied to the left flank to locate the left adrenal gland and the tumor. Serial cross-sectional images (0.076 mm between images) were taken. The tumor volume was measured using the 3-D reconstruction tool (Vevo Software v1.6.0, Toronto, Ontario, Canada).

Silk fibroin extraction

Silk fibroin from Bomby Mori silkworm cocoons from Tajima Shoji Co (Yokohama, Japan) was extracted as previously described [50]. Briefly, 5 g of cocoons were cut into approximately 1 cm × 1 cm pieces and boiled in 0.02 M Na₂CO₃ for 30 minutes to extract the silk fibroin. After drying overnight, the extracted silk fibroin was dissolved in 9.3 M LiBr at 60°C for 3 hrs. The dissolved silk fibroin was dialyzed against deionized water for 2 days. The aqueous silk fibroin solutions was stored at 4°C for future use.

Silk film fabrication

Silk films were fabricated as previously described [49, 51, 52]. Briefly, 192 μL of 4% (w/v) silk fibroin was pipetted onto 11 cm × 17 cm polydimethylsiloxane (Sylgard^®184, Dow Corning, Auburn, MI) and allowed to dry overnight. Seven millimeter diameter disks were cut from the dried silk films and autoclaved to induce the β-sheet transition and render the materials insoluble.

Halofuginone loading and drug release studies

Halofuginone-loaded silk films were formulated via an adsorption from solution method [51, 52]. Silk films (1.2 ± 0.04 mg) were loaded with Halofuginone via adsorption presumably through electrostatic interaction between the positively charge Halofuginone (basic pKa = 9.28; ChemAxon, Cambridge, MA) and negatively charge silk fibroin (isoelectric points ranging from 4.0 – 5.2 [53]). Halofuginone was dissolved in dimethyl sulfoxide at 50 mg/mL. Films were submerged in 1 mL of 31 μg/mL, 16 μg/mL and 8.5 μg/mL Halofuginone solution (diluted in sterile deionozed water) for 2 days. Halofuginone concentration was determined via UV/visible spectroscopy at 314 nm using a standard curve method (ranging from 1.6 μg/mL to 200 μg/mL). The Halofuginone loading was calculated as follows:

For drug release studies, Halofuginone-loaded silk films were placed into 0.5 mL phosphate buffered saline (pH7.4) at 37°C in a 1.5 mL Eppendorf tube. The absorbance of the PBS solution was read at each time point and solution replaced into the Eppendorf tube for study continuation. The PBS was not changed during the experimental time course to maintain the concentration above the minimal detectable concentration (~ 2 μg/mL). For drug release studies, silk films were loaded with 25 ± 1.3 μg, 13 ± 0.4 μg and 6 ± 0.7 μg Halofuginone (Figure 5A). Halofuginone exhibited sustained release over the first 24 hours (Figure 5B).

Abbreviations

PN: prominent nucleolar or prominent nucleolus, rRAN: ribosomal RNA.

Author contributions

Naohiko Ikegaki, Contribution: Concept development, experimental design for in vitro analysis, preparation of manuscript; Hiroyuki Shimada, Contribution: Concept development, pathology review, experimental design for pathological analysis, preparation of manuscript.

CONFLICTS OF INTEREST

The authors have no conflicts of interest to disclose in relation to manuscript.

FUNDING

This work was in part supported by NIH grants; U10CA180899, U10CA98543, PO1CA081403, NS094218, and R01CA127571.

REFERENCES

1. Gurney JG, Ross JA, Wall DA, Bleyer WA, Severson RK, Robison LL. Infant cancer in the U.S.: histology-specific incidence and trends, 1973 to 1992. J Pediatr Hematol Oncol. 1997; 19:428–32.

2. Cohn SL, Pearson AD, London WB, Monclair T, Ambros PF, Brodeur GM, Faldum A, Hero B, Iehara T, Machin D, Mosseri V, Simon T, Garaventa A, et al. The International Neuroblastoma Risk Group (INRG) classification system: an INRG Task Force report. J Clin Oncol. 2009; 27:289–97.

3. London WB, Castleberry RP, Matthay KK, Look AT, Seeger RC, Shimada H, Thorner P, Brodeur G, Maris JM, Reynolds CP, Cohn SL. Evidence for an age cutoff greater than 365 days for neuroblastoma risk group stratification in the Children’s Oncology Group. J Clin Oncol. 2005; 23:6459–65.

4. Brodeur GM, Pritchard J, Berthold F, Carlsen NL, Castel V, Castelberry RP, De Bernardi B, Evans AE, Favrot M, Hedborg F. Revisions of the international criteria for neuroblastoma diagnosis, staging, and response to treatment. Journal of Clinical Oncology. 1993; 11:1466–77.

5. Shimada H, Ambros IM, Dehner LP, Hata J, Joshi VV, Roald B. Terminology and morphologic criteria of neuroblastic tumors: recommendations by the International Neuroblastoma Pathology Committee. Cancer. 1999; 86:349–63.

6. Shimada H, Ambros IM, Dehner LP, Hata J, Joshi VV, Roald B, Stram DO, Gerbing RB, Lukens JN, Matthay KK, Castleberry RP. The International Neuroblastoma Pathology Classification (the Shimada system). Cancer. 1999; 86:364–72.

7. Shimada H, Umehara S, Monobe Y, Hachitanda Y, Nakagawa A, Goto S, Gerbing RB, Stram DO, Lukens JN, Matthay KK. International neuroblastoma pathology classification for prognostic evaluation of patients with peripheral neuroblastic tumors: a report from the Children’s Cancer Group. Cancer. 2001; 92:2451–61.

8. Peuchmaur M, d’Amore ES, Joshi VV, Hata J, Roald B, Dehner LP, Gerbing RB, Stram DO, Lukens JN, Matthay KK, Shimada H. Revision of the International Neuroblastoma Pathology Classification: confirmation of favorable and unfavorable prognostic subsets in ganglioneuroblastoma, nodular. Cancer. 2003; 98:2274–81. https://doi.org/10.1002/cncr.11773.

9. Brodeur GM, Seeger RC, Schwab M, Varmus HE, Bishop JM. Amplification of N-myc in untreated human neuroblastomas correlates with advanced disease stage. Science. 1984; 224:1121–4.

10. Seeger RC, Brodeur GM, Sather H, Dalton A, Siegel SE, Wong KY, Hammond D. Association of multiple copies of the N-myc oncogene with rapid progression of neuroblastomas. New England Journal of Medicine. 1985; 313:1111–6.

11. Look AT, Hayes FA, Shuster JJ, Douglass EC, Castleberry RP, Bowman LC, Smith EI, Brodeur GM. Clinical relevance of tumor cell ploidy and N-myc gene amplification in childhood neuroblastoma: a Pediatric Oncology Group study. J Clin Oncol. 1991; 9:581–91.

12. Attiyeh EF, London WB, Mosse YP, Wang Q, Winter C, Khazi D, McGrady PW, Seeger RC, Look AT, Shimada H, Brodeur GM, Cohn SL, Matthay KK, et al. Chromosome 1p and 11q deletions and outcome in neuroblastoma. N Engl J Med. 2005; 353:2243–53.

13. Weinstein JL, Katzenstein HM, Cohn SL. Advances in the diagnosis and treatment of neuroblastoma. Oncologist. 2003; 8:278–92.

14. Shimada H, Stram DO, Chatten J, Joshi VV, Hachitanda Y, Brodeur GM, Lukens JN, Matthay KK, Seeger RC. Identification of subsets of neuroblastomas by combined histopathologic and N-myc analysis. J Natl Cancer Inst. 1995; 87:1470–6.

15. Goto S, Umehara S, Gerbing RB, Stram DO, Brodeur GM, Seeger RC, Lukens JN, Matthay KK, Shimada H. Histopathology (International Neuroblastoma Pathology Classification) and MYCN status in patients with peripheral neuroblastic tumors: a report from the Children’s Cancer Group. Cancer. 2001; 92:2699–708.

16. Teshiba R, Kawano S, Wang LL, He L, Naranjo A, London WB, Seeger RC, Gastier-Foster JM, Look AT, Hogarty MD, Cohn SL, Maris JM, Park JR, et al. Age-dependent prognostic effect by Mitosis-Karyorrhexis Index in neuroblastoma: a report from the Children’s Oncology Group. Pediatr Dev Pathol. 2014; 17:441–9. https://doi.org/10.2350/14-06-1505-OA.1.

17. Suganuma R, Wang LL, Sano H, Naranjo A, London WB, Seeger RC, Hogarty MD, Gastier-Foster JM, Look AT, Park JR, Maris JM, Cohn SL, Amann G, et al. Peripheral neuroblastic tumors with genotype-phenotype discordance: a report from the Children’s Oncology Group and the International Neuroblastoma Pathology Committee. Pediatr Blood Cancer. 2013; 60:363–70. https://doi.org/10.1002/pbc.24238.

18. Wang LL, Teshiba R, Ikegaki N, Tang XX, Naranjo A, London WB, Hogarty MD, Gastier-Foster JM, Look AT, Park JR, Maris JM, Cohn SL, Seeger RC, et al. Augmented expression of MYC and/or MYCN protein defines highly aggressive MYC-driven neuroblastoma: a Children’s Oncology Group study. Br J Cancer. 2015; 113:57–63. https://doi.org/10.1038/bjc.2015.188.

19. Stine ZE, Walton ZE, Altman BJ, Hsieh AL, Dang CV. MYC, Metabolism, and Cancer. Cancer Discov. 2015; 5:1024–39. https://doi.org/10.1158/2159-8290.CD-15-0507.

20. Knoepfler PS, Zhang XY, Cheng PF, Gafken PR, McMahon SB, Eisenman RN. Myc influences global chromatin structure. Embo J. 2006; 25:2723–34. https://doi.org/10.1038/sj.emboj.7601152.

21. Dang CV. MYC, metabolism, cell growth, and tumorigenesis. Cold Spring Harb Perspect Med. 2013; 3. https://doi.org/10.1101/cshperspect.a014217.

22. Kim S, Li Q, Dang CV, Lee LA. Induction of ribosomal genes and hepatocyte hypertrophy by adenovirus-mediated expression of c-Myc in vivo. Proc Natl Acad Sci U S A. 2000; 97:11198–202. https://doi.org/10.1073/pnas.200372597.

23. Dang CV. MYC on the path to cancer. Cell. 2012; 149:22–35. https://doi.org/10.1016/j.cell.2012.03.003.

24. van Riggelen J, Yetil A, Felsher DW. MYC as a regulator of ribosome biogenesis and protein synthesis. Nat Rev Cancer. 2010; 10:301–9. https://doi.org/10.1038/nrc2819.

25. Ambros IM, Hata J, Joshi VV, Roald B, Dehner LP, Tuchler H, Potschger U, Shimada H. Morphologic features of neuroblastoma (Schwannian stroma-poor tumors) in clinically favorable and unfavorable groups. Cancer. 2002; 94:1574–83. https://doi.org/10.1002/cncr.10359.

26. Drygin D, Lin A, Bliesath J, Ho CB, O’Brien SE, Proffitt C, Omori M, Haddach M, Schwaebe MK, Siddiqui-Jain A, Streiner N, Quin JE, Sanij E, et al. Targeting RNA polymerase I with an oral small molecule CX-5461 inhibits ribosomal RNA synthesis and solid tumor growth. Cancer Res. 2011; 71:1418–30. https://doi.org/10.1158/0008-5472.CAN-10-1728.

27. Keller TL, Zocco D, Sundrud MS, Hendrick M, Edenius M, Yum J, Kim YJ, Lee HK, Cortese JF, Wirth DF, Dignam JD, Rao A, Yeo CY, et al. Halofuginone and other febrifugine derivatives inhibit prolyl-tRNA synthetase. Nat Chem Biol. 2012; 8:311–7. https://doi.org/10.1038/nchembio.790.

28. Otto T, Horn S, Brockmann M, Eilers U, Schuttrumpf L, Popov N, Kenney AM, Schulte JH, Beijersbergen R, Christiansen H, Berwanger B, Eilers M. Stabilization of N-Myc is a critical function of Aurora A in human neuroblastoma. Cancer Cell. 2009; 15:67–78.

29. Molenaar JJ, Domingo-Fernandez R, Ebus ME, Lindner S, Koster J, Drabek K, Mestdagh P, van Sluis P, Valentijn LJ, van Nes J, Broekmans M, Haneveld F, Volckmann R, et al. LIN28B induces neuroblastoma and enhances MYCN levels via let-7 suppression. Nat Genet. 2012; 44:1199–206. https://doi.org/10.1038/ng.2436.

30. Wang LL, Suganuma R, Ikegaki N, Tang X, Naranjo A, McGrady P, London WB, Hogarty MD, Gastier-Foster JM, Look AT, Park JR, Maris JM, Cohn SL, et al. Neuroblastoma of undifferentiated subtype, prognostic significance of prominent nucleolar formation, and MYC/MYCN protein expression: a report from the Children’s Oncology Group. Cancer. 2013; 119:3718–26. https://doi.org/10.1002/cncr.28251.

31. Birnstiel ML, Wallace H, Sirlin JL, Fischberg M. Localization of the ribosomal DNA complements in the nucleolar organizer region of Xenopus laevis. Natl Cancer Inst Monogr. 1966; 23:431–47.

32. Poortinga G, Quinn LM, Hannan RD. Targeting RNA polymerase I to treat MYC-driven cancer. Oncogene. 2015; 34:403–12. https://doi.org/10.1038/onc.2014.13.

33. Jinbo T, Iwamura Y, Kaneko M, Sawaguchi S. Coamplification of the L-myc and N-myc oncogenes in a neuroblastoma cell line. Jpn J Cancer Res. 1989; 80:299–301.

34. Hogarty MH, Brodeur GM. Oncogene amplification. In: Vogelstein B and Kinzler K, eds. Genetic Basis of Human Cancer. (New York: McGraw Hill).

35. Diolaiti D, McFerrin L, Carroll PA, Eisenman RN. Functional interactions among members of the MAX and MLX transcriptional network during oncogenesis. Biochim Biophys Acta. 2015; 1849:484–500. https://doi.org/10.1016/j.bbagrm.2014.05.016.

36. Bywater MJ, Pearson RB, McArthur GA, Hannan RD. Dysregulation of the basal RNA polymerase transcription apparatus in cancer. Nat Rev Cancer. 2013; 13:299–314. https://doi.org/10.1038/nrc3496.

37. Hein N, Hannan KM, George AJ, Sanij E, Hannan RD. The nucleolus: an emerging target for cancer therapy. Trends Mol Med. 2013; 19:643–54. https://doi.org/10.1016/j.molmed.2013.07.005.

38. Lindstrom MS, Deisenroth C, Zhang Y. Putting a finger on growth surveillance: insight into MDM2 zinc finger-ribosomal protein interactions. Cell Cycle. 2007; 6:434–7. https://doi.org/10.4161/cc.6.4.3861.

39. Dai MS, Arnold H, Sun XX, Sears R, Lu H. Inhibition of c-Myc activity by ribosomal protein L11. EMBO J. 2007; 26:3332–45. https://doi.org/10.1038/sj.emboj.7601776.

40. Liao JM, Zhou X, Gatignol A, Lu H. Ribosomal proteins L5 and L11 co-operatively inactivate c-Myc via RNA-induced silencing complex. Oncogene. 2014; 33:4916–23. https://doi.org/10.1038/onc.2013.430.

41. Harrison SJ, Khot A, Brajanovski N, Cameron D, Hein N, McArthur GA, Lim JKC, O’Brien S, Ryckman DM, Yu GI, Link E, Donohoe C, Snowden A, et al. A phase 1, open-label, dose escalation, safety, PK and PD study of a first in class Pol1 inhibitor (CX-5461) in patients with advanced hematologic malignancies (HM). Journal of Clinical Oncology. 2015; 33:e22212-e. https://doi.org/10.1200/jco.2015.33.15_suppl.e22212.

42. Sundrud MS, Koralov SB, Feuerer M, Calado DP, Kozhaya AE, Rhule-Smith A, Lefebvre RE, Unutmaz D, Mazitschek R, Waldner H, Whitman M, Keller T, Rao A. Halofuginone inhibits TH17 cell differentiation by activating the amino acid starvation response. Science. 2009; 324:1334–8. https://doi.org/10.1126/science.1172638.

43. Ikegaki N, Bukovsky J, Kennett RH. Identification and characterization of the NMYC gene product in human neuroblastoma cells by monoclonal antibodies with defined specificities. Proc Natl Acad Sci U S A. 1986; 83:5929–33.

44. Kaplan EL, Meier P. Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958; 53:457–81.

45. Peto R, Pike MC, Armitage P, Breslow NE, Cox DR, Howard SV, Mantel N, McPherson K, Peto J, Smith PG. Design and analysis of randomized clinical trials requiring prolonged observation of each patient. II. analysis and examples. Br J Cancer. 1977; 35:1–39.

46. Ikegaki N, Hicks SL, Regan PL, Jacobs J, Jumbo AS, Leonhardt P, Rappaport EF, Tang XX. S(+)-ibuprofen destabilizes MYC/MYCN and AKT, increases p53 expression, and induces unfolded protein response and favorable phenotype in neuroblastoma cell lines. Int J Oncol. 2014; 44:35–43. https://doi.org/10.3892/ijo.2013.2148.

47. Torres J, Regan PL, Edo R, Leonhardt P, Jeng EI, Rappaport EF, Ikegaki N, Tang XX. Biological effects of induced MYCN hyper-expression in MYCN-amplified neuroblastomas. Int J Oncol. 2010; 37:983–91.

48. Regan PL, Jacobs J, Wang G, Torres J, Edo R, Friedmann J, Tang XX. Hsp90 inhibition increases p53 expression and destabilizes MYCN and MYC in neuroblastoma. Int J Oncol. 2011; 38:105–12.

49. Chiu B, Coburn J, Pilichowska M, Holcroft C, Seib FP, Charest A, Kaplan DL. Surgery combined with controlled-release doxorubicin silk films as a treatment strategy in an orthotopic neuroblastoma mouse model. Br J Cancer. 2014; 111:708–15. https://doi.org/10.1038/bjc.2014.324.

50. Rockwood DN, Preda RC, Yucel T, Wang X, Lovett ML, Kaplan DL. Materials fabrication from Bombyx mori silk fibroin. Nat Protoc. 2011; 6:1612–31. https://doi.org/10.1038/nprot.2011.379.

51. Coburn JM, Na E, Kaplan DL. Modulation of vincristine and doxorubicin binding and release from silk films. J Control Release. 2015; 220:229–38. https://doi.org/10.1016/j.jconrel.2015.10.035.

52. Seib FP, Coburn J, Konrad I, Klebanov N, Jones GT, Blackwood B, Charest A, Kaplan DL, Chiu B. Focal therapy of neuroblastoma using silk films to deliver kinase and chemotherapeutic agents in vivo. Acta Biomater. 2015; 20:32–8. https://doi.org/10.1016/j.actbio.2015.04.003.

53. Song ZT, Speakman PT. Isoelectric focusing of silk fibroin polypeptides. Biochemical Society Transactions. 1985; 13:387. https://doi.org/10.1042/bst0130387.