Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
    • CME
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • CCR Focus Archive
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Breast Cancer
      • Clinical Trials
      • Immunotherapy: Facts and Hopes
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Clinical Cancer Research
Clinical Cancer Research
  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
    • CME
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • CCR Focus Archive
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Breast Cancer
      • Clinical Trials
      • Immunotherapy: Facts and Hopes
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

Human Cancer Biology

Chromosomal and MicroRNA Expression Patterns Reveal Biologically Distinct Subgroups of 11q− Neuroblastoma

Patrick G. Buckley, Leah Alcock, Kenneth Bryan, Isabella Bray, Johannes H. Schulte, Alexander Schramm, Angelika Eggert, Pieter Mestdagh, Katleen De Preter, Jo Vandesompele, Frank Speleman and Raymond L. Stallings
Patrick G. Buckley
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leah Alcock
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenneth Bryan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Isabella Bray
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Johannes H. Schulte
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander Schramm
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Angelika Eggert
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pieter Mestdagh
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katleen De Preter
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jo Vandesompele
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Frank Speleman
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raymond L. Stallings
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1158/1078-0432.CCR-09-3215 Published June 2010
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Purpose: The purpose of this study was to further define the biology of the 11q− neuroblastoma tumor subgroup by the integration of array-based comparative genomic hybridization with microRNA (miRNA) expression profiling data to determine if improved patient stratification is possible.

Experimental Design: A set of primary neuroblastoma (n = 160), which was broadly representative of all genetic subtypes, was analyzed by array-based comparative genomic hybridization and for the expression of 430 miRNAs. A 15-miRNA expression signature previously shown to be predictive of clinical outcome was used to analyze an independent cohort of 11q− tumors (n = 37).

Results: Loss of 4p and gain of 7q occurred at a significantly higher frequency in the 11q− tumors, further defining the genetic characteristics of this subtype. The 11q− tumors could be split into two subgroups using a miRNA expression survival signature that differed significantly in clinical outcome and the overall frequency of large-scale genomic imbalances, with the poor survival subgroup having significantly more imbalances. miRNAs from the expression signature, which were upregulated in unfavorable tumors, were predicted to target downregulated genes from a published mRNA expression classifier of clinical outcome at a higher-than-expected frequency, indicating the miRNAs might contribute to the regulation of genes within the signature.

Conclusion: We show that two distinct biological subtypes of neuroblastoma with loss of 11q occur, which differ in their miRNA expression profiles, frequency of segmental imbalances, and clinical outcome. A miRNA expression signature, combined with an analysis of segmental imbalances, provides greater prediction of event-free survival and overall survival outcomes than 11q status by itself, improving patient stratification. Clin Cancer Res; 16(11); 2971–8. ©2010 AACR.

Translational Relevance

Loss of chromosome 11q material in neuroblastoma is generally associated with poor clinical outcome. The results presented in this report, however, show that using miRNA expression profiling, the 11q− genetic subgroup can be subdivided into two distinct groups. These subgroups differ significantly in clinical outcome and with respect to the overall frequency of segmental genomic imbalances, indicating the possibility for more refined patient stratification before therapy.

Neuroblastoma, one of the most common solid tumors in children, is derived from precursor cells of the sympathetic nervous system and is responsible for approximately 15% of all childhood cancer deaths (see ref. 1 for review). These tumors are particularly noted for extensive heterogeneity in clinical behavior, ranging from spontaneous regression to aggressive clinical course and death due to disease. Tumors in infants characterized by whole chromosome gains and losses and few structural chromosome abnormalities usually have favorable clinical outcomes and, in some instances, can even spontaneously regress. In contrast, tumors from children >1.5 years of age are often metastatic at diagnosis and become refractory to treatment. These unfavorable tumor subtypes possess a broad spectrum of recurrent chromosomal gains, losses, and amplifications, many of which are highly associated with clinical outcome and form the basis for dividing neuroblastoma into genetic subtypes.

Tumors with amplification of the MYCN transcription factor represent a major genetic subtype of metastatic neuroblastoma with a particularly poor prognosis (2). Loss of chromosome 1p and gain of chromosome 17q regions are common in MYCN-amplified (MNA) tumors, often occurring by a common chromosomal mechanism involving an unbalanced t(1p;17q) (3). Tumors with hemizygous loss of a large segment on chromosome 11q represent another major genetic subtype of metastatic neuroblastoma (4, 5), which can also be associated with unfavorable clinical outcome (6, 7). Loss of chromosome 3p occurs preferentially in 11q− tumors, and interestingly, children with 3p−/11q− tumors have a higher median age at diagnosis than do those with MNA tumors (6, 8, 9). Similar to MNA tumors, the 11q− tumors usually have gain of 17q, but as a consequence of an unbalanced t(11q;17q), which simultaneously results in loss of the 11q material (10, 11). At the molecular level, MNA and 11q− tumors are clearly different, being distinguishable by both mRNA (12, 13) and microRNA (miRNA; ref. 14) genome-wide expression profiling. There are also many additional recurring large-scale chromosome imbalances found in neuroblastomas, including loss of 4p, 9p, and 14q and gain of 1q, 7q, 2p, and 11p (see ref. 15 for review). Approximately 70% to 80% of all metastatic (stage IV) neuroblastomas have neither MNA nor loss of 11q material and form a much more genetically heterogeneous group of tumors (15–17).

Although the 11q− subtype has been shown to be associated with poor clinical outcome (6, 7, 15), exactly what other genetic factors might mediate patient survival within this subgroup has remained poorly defined. Here, we have carried out an integrated analysis of DNA copy number and miRNA expressional profiling on a set of primary neuroblastomas to further define the clinical and genetic characteristics of this important tumor subtype.

Materials and Methods

Tumor samples

Neuroblastoma tumor samples were obtained retrospectively from Our Lady's Children's Hospital, Dublin, Ireland (n = 49), and the Children's Oncology Group Neuroblastoma Tumor Bank (n = 111), for a total of 160 tumors [array-based comparative genomic hybridization (aCGH) profiling was carried on the entire set of tumors, whereas miRNA profiling is available on a subset (n = 144)]. Additional 11q− tumors were also obtained retrospectively from University Children's Hospital, Essen, Germany (n = 12), and Ghent University Hospital, Ghent, Belgium (n = 11), for validation purposes. These tumors were also characterized by aCGH and by miRNA expression profiling. Patients were treated under either the European treatment protocol or the U.S. neuroblastoma treatment protocol between 1998 and 2004. Clinical details are summarized in Supplementary Table S1 and Fig. S1. Event-free survival (EFS) time was calculated from the time of enrollment on the front-line or biological study until the time of the first occurrence of relapse, progressive disease, secondary malignancy, or death, or until the time of last contact if no event occurred. Overall survival (OS) time was calculated until the time of death or until last contact. EFS and OS are presented as the estimate ± SE.

miRNA expression profiling

Tumors from all centers were profiled on the same Applied Biosystems platform at Ghent University Hospital as described by Mestdagh et al. (18) and normalized by the same methods (19).

Oligonucleotide aCGH

For tumors from Dublin and Children's Oncology Group, DNA was extracted from tumor samples using the DNeasy Kit (Qiagen). DNA was quantified using the Nanodrop spectrophotometer and stored at −80°C. Methods for array hybridization and analysis were as previously reported (20) and summarized in NimbleGen aCGH protocol v5.1, with minor modifications. Briefly, a total of 500 ng of tumor DNA and control DNA were differentially labeled with Cy3 and Cy5 (TriLink Biotechnologies), respectively. Four micrograms of labeled test and reference were hybridized to a 4-plex 72,000 probe array (Roche NimbleGen) for 18 hours in a MAUI hybridization station (BioMicro Systems) at 42°C. Posthybridization washes were done and arrays scanned using the Axon 4000B scanner and GenePix 6.0 (Molecular Devices). Fluorescent intensities were extracted and log 2 ratios calculated and qspline normalized (21) using NimbleScan Software version 2.4. Aberration breakpoints were identified using the CGH-segMNT algorithm in NimbleScan 2.4. Unaveraged log 2 ratios were written out as .gff files and used for data visualization. In general, a significance log 2 ratio threshold of <−0.25 for loss and >0.25 for gain was used to identify DNA copy number imbalances. For tumors from Essen and Ghent, a custom Agilent 44K array was applied for aCGH using standard protocols.

Bioinformatic analysis

Unsupervised computational analysis

Unsupervised partitional clustering (k-means) was used to optimize a k = 2 split in the 11q− validation data based solely on the expression within the 15-miRNA signature. Pearson's r correlation coefficient was used as a sample object similarity metric. A heatmap was subsequently produced to illustrate this split. These methods were implemented using the k-means and heatmap functions from the stats package of R version 2.10.0.

Patient survival analysis

Kaplan-Meier graphs and log-rank tests for significance were carried out using GraphPad Prism software version 5. Cox proportional hazards regression model was applied to assess the independent effect of each chromosomal aberration on survival. 17q gain was screened for independence of MNA and 11q loss for both EFS and OS. In addition to this model, a screen for independence of 1q loss and 1q gain was also done for both EFS and OS. These methods were implemented using the coxph function in the survival package in R version 2.10.0.

Statistical analysis of miRNA and mRNA signatures

Data from the mirBase miRNA targets database (Microcosm; ref. 22) were used to compile a two-dimensional miRNA/mRNA matrix containing 710 miRNAs, 34,525 transcripts, and 728,289 predicted interactions. All significant (P ≤ 0.05) predicted binding interactions connecting the 15 miRNAs to the various published gene signatures and the average prediction scores were then retrieved. The significance of these interactions was assessed by examining how often these combinations might occur by chance. For each miRNA/mRNA target combination, 100,000 submatrices of equal size were randomly generated and the numbers of connections and scores assessed. The resultant P value was thus used to evaluate the significance of relationship and possible influence the 15-miRNA signature might have on each of the published gene signatures. The above assessment was carried out using in-house generated software implemented in Java version 6 update 15.

Visualization of chromosomal variations and shortest regions of overlap

Global visualization of chromosomal variations for all 160 samples was carried out using an in-house software application implemented using Java and made use of the JFreeChart and Commons packages. Common copy number variants were excluded before visualization. The most commonly detected shortest regions of overlap (SRO) were identified and genomic positions mapped using this application.

Results

aCGH analysis

DNA copy number analysis was done on 121 primary neuroblastoma tumors using a 72k 4-plex aCGH array. This data set was then combined with a previously reported data set on 39 tumors (23) obtained with a 385k array format for a total analysis of 160 tumors. Figure 1 displays a summary of the aberrations identified across the genome in these tumors, whereas Supplementary Table S1 provides a comprehensive summary of the aCGH results and clinical details for each tumor.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Summary of DNA copy number aberrations across the genome. Regions of hemizygous loss (red bars) and gains (green bars) are plotted for each autosome with the observed frequency (percent) of each aberration detected in the tumor series plotted on the Y-axis. The frequency of MYCN and ALK amplifications is also plotted (blue).

All previously published (4, 6, 8, 9) positive and negative correlations between different genetic abnormalities were detected in our tumor cohort with statistical significance (Table 1). These included inverse correlation between MYCN amplification (MNA) and chromosome 11q loss (P < 0.0001) and positive correlations between loss of 1p and MNA (P = 0.008) and 3p loss with 11q loss (0.0006). Moreover, due to the large size of our tumor set, we were able to identify novel positive correlations between loss of 4p and gain of 7q with loss of 11q (P = 0.002 and P = 0.009, respectively, corrected for multiple comparisons), allowing further definition of the 11q− subtype (Table 1). Loss of 3p and loss 4p were not inversely related within the 11q− subgroup, indicating that these imbalances did not have a complementary effect on tumorigenesis.

View this table:
  • View inline
  • View popup
Table 1.

Association of chromosomal aberrations with 11q− and MNA tumors

Prior studies have indicated that MYCN amplification or 11q loss occurs in approximately 70% to 80% of patients presenting with stage IV disease (15–17). Tumors that lack either of these abnormalities are a much more genetically heterogeneous group that lack predictive biomarkers. The most frequent imbalances found in tumors from our data set lacking either MNA or 11q loss are displayed in Supplementary Table S2 (stage IV disease, n = 20; stages I, II, III, and IV, n = 39). The three most common abnormalities in this subgroup were 17q gain (n = 25; 40%), 16p loss (n = 20; 32%), and 19p loss (n = 19; 31%). Some imbalances, particularly 1q gain, were significantly underrepresented in non-MNA/11q− tumors (n = 1; P = 0.01, corrected for multiple comparisons). There were no discernable positive correlations between genomic imbalances in the non-MNA/11q− tumor group.

Our tumor cohort also exhibited previously published associations between clinical details, genetic abnormalities, and clinical outcome (Supplementary Figs. S1 and S2). As expected, multivariate analysis using the Cox proportional hazards regression model indicated that 1q gain (HR, 2.0; P = 0.04; OS only), 11q loss (HR, 2.1; P = 0.03; EFS only), and MNA (HR, 2.7; P = 0.006; EFS; HR, 2.7; P = 0.008; OS) were independently associated with poor patient survival, although a number of other abnormalities were significantly associated with poor survival in univariate analyses (Supplementary Fig. S2). Unlike the majority of previously published reports (24–26), 17q gain was not an independent predictor of poor survival, as assessed by the Cox proportional hazards regression model. This result, however, is consistent with that of Spitz et al. (27) who reported that 17q gain, in the absence of MNA or 11q−, was not a predictor of EFS or OS.

Shortest regions of common gain or loss of genetic material in tumors

Regions affected by DNA copy number alterations were further analyzed to identify regions commonly deleted or gained in the tumor series. Supplementary Table S3 summarizes these SROs with the approximate genomic position of each breakpoint, the observed frequency of each aberration, and comparison of SROs previously reported. Regions that were classified as copy number variant in the Database of Genomic Variants (http://projects.tcag.ca/variation/) were excluded from this analysis. Detected SROs ranged in size from 0.9 to 10.4 Mb for losses and 0.3 to 6.6 Mb for DNA copy number gains. The most frequent SRO identified was detected on 17q, encompassing a region of 5.2 Mb. This region of common gain contains 18 annotated genes including 5 members of the ATP-binding cassette, subfamily A transporter family (ABCA5, ABCA6, ABCA8, ABCA9, and ABCA10). Interestingly, overexpression of transcripts in this gene family has been previously reported in chronic lymphocytic leukemia (28), high-grade prostatic intraepithelial neoplasia (29), and breast cancer (30) and have also been shown to confer drug resistance on tumor cells (31). Other cancer-related genes located within this 17q SRO include SOX9 (32) and SSTR2 (33).

Two SROs (gains) were identified on chromosome 1q encompassing regions of 1.2 and 3 Mb, respectively. The 3-Mb SRO imbalance was further refined through integration of a previously reported common region of gain on 1q (17). This 1q SRO has been mapped to 206.4-207.7 Mb and contains the genes PLXNA2, mir-205, and LOC642587. PLXNA2 is a member of the plexin family of proteins, which function in association with semaphorin proteins during neural development (34). Although mir-205 has been previously implicated in cancer development (35–37), a screen of 430 miRNAs in 145 neuroblastoma tumors did not identify this miRNA as significant in the pathogenesis of neuroblastoma (38).

Dichotomization of the 11q− genetic subtype by miRNA expression profiling

Fischer et al. (39) recently reported that a 144-gene expression signature originally developed by Oberthuer et al. (40) could stratify 11q− tumors into two subgroups: one with favorable clinical outcome and a second group with unfavorable outcome. To determine if our cohort of 11q− tumors exhibited a similar dichotomous pattern of clinical outcome, we analyzed 11q− tumors (n = 36) using a 15-miRNA expression profiling signature that is predictive of survival, as reported by Bray et al. (38). The signature was trained to identify poor survival on a set of neuroblastomas that included all genetic subtypes of the disease, and then validated on a similar set of tumors that was also broadly representative. The set of 11q− tumors used in the current study is completely independent of the tumors used by Bray et al. (38) to derive the miRNA survival signature. The unsupervised partitional clustering algorithm, k-means, was used to divide the 11q− tumors into two subgroups on the basis of miRNA expression. The two 11q− subgroups differed significantly with respect to EFS (P = 0.04) and OS (P = 0.03; Fig. 2). For group A (n = 16), the probability for 4-year EFS was 14.5% ± 9.3 SE, whereas this increased to 54.7% ± 12.8 SE for group B. The probability of 4-year OS for group A was 25.5% ± 12.5 SE, whereas for group B the rate of survival increased to 65.0% ± 12.8 SE.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Stratification of 11q− tumors based on a miRNA survival signature. Thirty-seven 11q− tumors could split by the unsupervised k-means partitional algorithm into two groups (A and B) on the basis of a 15-miRNA survival signature previously described by Bray et al. (38). Median survival time was 800 d for group A and 1,445 d for group B. Information related to death due to disease, events, age at diagnosis, and disease stage is presented at the top of the 15-miRNA heatmap. Information on genomic imbalances present in each tumor (black) is presented below the heatmap. Red on the heatmap is associated with higher expression, and green with lower expression. The number of genomic imbalances for the unfavorable cluster (A) was significantly higher (1.5-fold) than that for the favorable cluster (B; P = 0.015). Kaplan-Meier curves for EFS and OS show a statistically significant difference for patient outcome between the two groups.

Examination of the genomic imbalances found in the favorable (group B) versus unfavorable (group A) 11q− tumor subgroups indicated that no single imbalance was overrepresented with statistical significance. Gain of 17q was present at nearly equal frequencies in the two groups and could not account for the difference in survival for the two different groups of patients with 11q− tumors, consistent with our observation that 17q gain is not independently predictive of survival. Nevertheless, the overall frequency of genomic imbalances was significantly higher in the unfavorable tumor subgroup A (mean 5.5 versus 3.8; P = 0.019; Fig. 2), indicating the possibility that an overall increase in imbalances is associated with poor survival. Interestingly, two of the patients that were identified as having favorable tumors on the basis of miRNA profiling (i.e., false negative results) had a higher number of segmental imbalances (n = 7 in each tumor) than the mean number (5.5) found in the unfavorable tumor subgroup. Thus, an analysis of segmental imbalances following miRNA expression profiling could potentially improve OS prediction. This combined approach led to a sensitivity of 82.4%, a specificity of 55.0%, a positive predictive value of 77.7%, and a negative predictive value of 75.0% on our data set. We conclude that 11q− tumors are composed of two subgroups with very divergent clinical behavior and that these subgroups can be distinguished on the basis of miRNA expression profiling and analysis of DNA copy number alterations.

There was no relationship between patient age at diagnosis and survival for the 11q− tumor data set, as only 6 patients with 11q− tumors were <1.5 years, with 3 of these cases occurring in each group. The relatively small number of patients <1.5 years of age at diagnosis, which is generally associated with better survival, was not surprising given that 11q loss is significantly correlated with older patients (6, 15). In addition, there was also no statistically significant difference in tumor stage between the unfavorable and the favorable subgroups.

Next, we determined whether any of the 15 miRNAs in our survival signature target the genes in the 144-gene signature used by Fischer et al. (39). For the downregulated miRNAs associated with poor survival, 10 potential targets with upregulated genes were identified. However, the number and significance of predicted interactions were not significantly overrepresented (P = 0.20) compared with background (see Materials and Methods). For the upregulated miRNAs associated with poor survival, 13 potential targets with downregulated genes were identified, which approached statistical significance (P = 0.06). More recently, Vermeulen et al. (41) have described a gene signature predictive of clinical outcome for neuroblastoma. A statistically significant number of interactions between downregulated miRs from our signature and upregulated genes from this signature was not evident. However, 18 putative interactions between upregulated miRNAs and downregulated genes from the gene signature were identified, which, according to our simulation studies, is significantly higher than background (P = 0.0007). Some miRNAs from our signature that are upregulated in unfavorable tumors, such as miR-572 (Supplementary Table S4), are predicted to target multiple downregulated genes from the Vermeulen signature (DDC, PRKCZ, PIK3R1, and PTPRF), further indicating that the miRNA signature regulates the gene survival signature to at least some extent. Caution is required in the interpretation of these results, however, as these are not experimentally validated interactions.

Discussion

Neuroblastoma tumors with loss of chromosome 11q material represent a major genetic subgroup of tumors with adverse clinical outcome. Loss of chromosome 3p was previously shown to occur preferentially within this tumor subgroup (6, 8, 9), and we now further show that loss of chromosome 4p and gain of chromosome 7q are significantly overrepresented in the 11q− tumor subset. Although our results further refine the shortest region of loss or gain on each of these chromosomes, the genes contributing to tumorigenesis on each of these chromosomes remain unidentified. To a lesser extent, 2p gain, 11p gain, and 14q loss are also preferentially associated with the 11q− tumor subtype.

We also show by miRNA expression profiling that 11q− tumors are actually composed of two clinical subtypes with significant differences in both EFS and OS. This finding is completely consistent with a recent report by Fischer et al. (39) showing that 11q− tumors can be split into two subgroups with distinctly different clinical outcomes on the basis of mRNA expression profiling. Moreover, we show by computational methods that miRNAs from our expression survival signature show an enrichment for predicted target sites in the 3′ untranslated region of mRNAs from two published gene survival signatures for neuroblastoma: Obertheur et al. (P = 0.06; ref. 40) and Vermeulen et al. (P < 0.0007; ref. 41). Therefore, miRNAs from the survival signature (38) could potentially influence the gene expression signatures predictive of aggressive disease pathogenesis, at least to a limited extent. It is noteworthy that Mestdagh et al. (42) have also recently shown a significant correspondence between miRNAs that are upregulated either directly or indirectly by MYCN/cMYC and genes that are significantly downregulated in unfavorable neuroblastoma (42). Their miRNA signature was based primarily on MYCN/cMYC regulated miRNAs and therefore differs significantly from the signature used by Bray et al. (38), which was trained solely on the basis of miRNAs associated with poor patient survival. This indicates that differences in miRNA profiles within the 11q− subgroups are not solely MYC driven.

Indeed, the mechanisms regulating the miRNAs in the survival signature are likely to be varied and complex. Bray et al. (38) have shown that expression levels of many miRNA loci are affected by MYCN levels and by large-scale losses and gains of genomic regions in neuroblastoma tumors. The unfavorable 11q− tumors identified by miRNA expression profiling have a significantly higher number of genomic imbalances, relative to the favorable 11q− subgroup, indicating that genomic imbalances might be contributing to the differences in miRNA expression in the two different 11q− subgroups.

Unbalanced gain of 17q, which has been reported to be predictive of poor patient survival by a number of groups (24–26), occurs at nearly equal frequencies in the favorable and unfavorable 11q− subgroups, leading us to conclude that 17q gain does not contribute to the differences in EFS and OS between these groups in a large way. This is consistent with our multivariate analysis of the entire tumor cohort indicating that 17q gain is not independently predictive of poor clinical outcome. Our result is consistent with a previous report by Spitz et al. indicating that 17q gain is not predictive of survival in tumors lacking MNA or 11q loss (27). Thus, our results do not support a role for 17q gain in contributing to an aggressive tumor phenotype, although the overall high frequency of 17q gain in neuroblastoma (including tumors that lack MNA or 11q loss) indicates that it contributes to tumorigenesis in some manner, such as early tumor initiation.

Our finding of an increased frequency of segmental imbalances in the unfavorable 11q− tumor subgroup is completely consistent with reports that the overall pattern of genomic imbalances is highly predictive of clinical outcome in neuroblastoma (43, 44). Janoueix-Lerosey et al. (43) concluded that the presence of segmental alterations is a powerful predictor of poor clinical outcome and can be used to improve patient stratification. Our demonstration that an increased frequency of segmental alterations in 11q− tumors is also correlated with lower EFS and OS provides further confirmation that the two subgroups identified by the miRNA profiling are biologically distinct.

Attiyeh et al. (7) reported following the analysis of a large neuroblastoma tumor cohort that 11q loss by itself is independently predictive of EFS and progression-free survival but is not an independent predictor of OS. Loss of 11q was also significantly associated with EFS, but not OS, in our data set. The specific application of a 15-miRNA expression signature, combined with an analysis of segmental imbalances, however, allows for the identification of poor OS for patients with 11q− tumors with a sensitivity of 82%. The signature seems to be better at identifying truly aggressive 11q− tumors than genuinely favorable tumors, as specificity was only 55%.

In conclusion, we show that two 11q− neuroblastoma tumor subgroups occur, which differ significantly in EFS and OS. These subgroups differ significantly in their miRNA expression profiles and, as shown by Fischer et al. (39), in their mRNA expression profiles. The tumor subtype with poor EFS and OS has significantly more segmental genomic imbalances than the favorable subgroup, indicating that a combination of expression profiling (miRNAs and mRNAs) with analysis of DNA copy number alterations will lead to improved prognostication of this often fatal tumor subtype.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Acknowledgments

Grant Support: Grants from Science Foundation Ireland (07/IN.1/B1776), Children's Medical and Research Foundation, the NIH (5R01CA127496), the Belgian Foundation Against Cancer, Foundation of public interest, the Fund for Scientific Research (G.0198.08), GOA (01G01910), and the VLK (Vlaamse Liga Tegen Kanker). P.G. Buckley is supported by a postdoctoral fellowship from the Irish Research Council for Science, Engineering and Technology. K. De Preter is a postdoctoral researcher with the Fund for Scientific Research-Flanders. P. Mestdagh was sponsored by BOF (01D31406).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

  • Received December 8, 2009.
  • Revision received March 31, 2010.
  • Accepted April 15, 2010.

References

  1. ↵
    1. Brodeur GM
    . Neuroblastoma: biological insights into a clinical enigma. Nat Rev Cancer 2003;3:203–16.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Brodeur GM,
    2. Seeger RC,
    3. Schwab M,
    4. Varmus HE,
    5. Bishop JM
    . Amplification of N-myc in untreated human neuroblastomas correlates with advanced disease stage. Science 1984;224:1121–4.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Van Roy N,
    2. Laureys G,
    3. Cheng NC,
    4. et al
    . 1;17 translocations and other chromosome 17 rearrangements in human primary neuroblastoma tumors and cell lines. Genes Chromosomes Cancer 1994;10:103–14.
    OpenUrlPubMed
  4. ↵
    1. Guo C,
    2. White PS,
    3. Weiss MJ,
    4. et al
    . Allelic deletion at 11q23 is common in MYCN single copy neuroblastomas. Oncogene 1999;18:4948–57.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Plantaz D,
    2. Vandesompele J,
    3. Van Roy N,
    4. et al
    . Comparative genomic hybridization (CGH) analysis of stage 4 neuroblastoma reveals high frequency of 11q deletion in tumors lacking MYCN amplification. Int J Cancer 2001;91:680–6.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Spitz R,
    2. Hero B,
    3. Ernestus K,
    4. Berthold F
    . Deletions in chromosome arms 3p and 11q are new prognostic markers in localized and 4s neuroblastoma. Clin Cancer Res 2003;9:52–8.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Attiyeh EF,
    2. London WB,
    3. Mosse YP,
    4. et al
    . Chromosome 1p and 11q deletions and outcome in neuroblastoma. N Engl J Med 2005;353:2243–53.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Vandesompele J,
    2. Van Roy N,
    3. Van Gele M,
    4. et al
    . Genetic heterogeneity of neuroblastoma studied by comparative genomic hybridization. Genes Chromosomes Cancer 1998;23:141–52.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Breen CJ,
    2. O'Meara A,
    3. McDermott M,
    4. Mullarkey M,
    5. Stallings RL
    . Coordinate deletion of chromosome 3p and 11q in neuroblastoma detected by comparative genomic hybridization. Cancer Genet Cytogenet 2000;120:44–9.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Stark B,
    2. Jeison M,
    3. Glaser-Gabay L,
    4. et al
    . der(11)t(11;17): a distinct cytogenetic pathway of advanced stage neuroblastoma (NBL)—detected by spectral karyotyping (SKY). Cancer Lett 2003;197:75–9.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Stallings RL,
    2. Carty P,
    3. McArdle L,
    4. et al
    . Molecular cytogenetic analysis of recurrent unbalanced t(11;17) in neuroblastoma. Cancer Genet Cytogenet 2004;154:44–51.
    OpenUrlCrossRefPubMed
  12. ↵
    1. McArdle L,
    2. McDermott M,
    3. Purcell R,
    4. et al
    . Oligonucleotide microarray analysis of gene expression in neuroblastoma displaying loss of chromosome 11q. Carcinogenesis 2004;25:1599–609.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Wang Q,
    2. Diskin S,
    3. Rappaport E,
    4. et al
    . Integrative genomics identifies distinct molecular classes of neuroblastoma and shows that multiple genes are targeted by regional alterations in DNA copy number. Cancer Res 2006;66:6050–62.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Chen Y,
    2. Stallings RL
    . Differential patterns of microRNA expression in neuroblastoma are correlated with prognosis, differentiation, and apoptosis. Cancer Res 2007;67:976–83.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Vandesompele J,
    2. Baudis M,
    3. De Preter K,
    4. et al
    . Unequivocal delineation of clinicogenetic subgroups and development of a new model for improved outcome prediction in neuroblastoma. J Clin Oncol 2005;23:2280–99.
    OpenUrlAbstract/FREE Full Text
    1. Michels E,
    2. Vandesompele J,
    3. De Preter K,
    4. et al
    . ArrayCGH-based classification of neuroblastoma into genomic subgroups. Genes Chromosomes Cancer 2007;46:1098–108.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Mosse YP,
    2. Diskin SJ,
    3. Wasserman N,
    4. et al
    . Neuroblastomas have distinct genomic DNA profiles that predict clinical phenotype and regional gene expression. Genes Chromosomes Cancer 2007;46:936–49.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Mestdagh P,
    2. Feys T,
    3. Bernard N,
    4. et al
    . High-throughput stem-loop RT-qPCR miRNA expression profiling using minute amounts of input RNA. Nucleic Acids Res 2008;36:e143.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Mestdagh P,
    2. Van Vlierberghe P,
    3. De Weer A,
    4. et al
    . A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol 2009;10:R64.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Selzer RR,
    2. Richmond TA,
    3. Pofahl NJ,
    4. et al
    . Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH. Genes Chromosomes Cancer 2005;44:305–19.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Workman C,
    2. Jensen LJ,
    3. Jarmer H,
    4. et al
    . A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biol 2002;3:research0048.
    OpenUrlPubMed
  21. ↵
    1. Griffiths-Jones S,
    2. Saini HK,
    3. van Dongen S,
    4. Enright AJ
    . miRBase: tools for microRNA genomics. Nucleic Acids Res 2008;36:D154–8.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Stallings RL,
    2. Nair P,
    3. Maris JM,
    4. et al
    . High-resolution analysis of chromosomal breakpoints and genomic instability identifies PTPRD as a candidate tumor suppressor gene in neuroblastoma. Cancer Res 2006;66:3673–80.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Bown N,
    2. Cotterill S,
    3. Lastowska M,
    4. et al
    . Gain of chromosome arm 17q and adverse outcome in patients with neuroblastoma. N Engl J Med 1999;340:1954–61.
    OpenUrlCrossRefPubMed
    1. Bown N,
    2. Lastowska M,
    3. Cotterill S,
    4. et al.
    ,U.K. Cancer Cytogenetics GroupU.K. Children's Cancer Study Group. 17q gain in neuroblastoma predicts adverse clinical outcome. Med Pediatr Oncol 2001;36:14–9.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Brinkschmidt C,
    2. Christiansen H,
    3. Terpe HJ,
    4. et al
    . Distal chromosome 17 gains in neuroblastomas detected by comparative genomic hybridization (CGH) are associated with a poor clinical outcome. Med Pediatr Oncol 2001;36:11–3.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Spitz R,
    2. Hero B,
    3. Ernestus K,
    4. Berthold F
    . Gain of distal chromosome arm 17q is not associated with poor prognosis in neuroblastoma. Clin Cancer Res 2003;9:4835–40.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Jelinek DF,
    2. Tschumper RC,
    3. Stolovitzky GA,
    4. et al
    . Identification of a global gene expression signature of B-chronic lymphocytic leukemia. Mol Cancer Res 2003;1:346–61.
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Hu Y,
    2. Wang M,
    3. Veverka K,
    4. Garcia FU,
    5. Stearns ME
    . The ABCA5 protein: a urine diagnostic marker for prostatic intraepithelial neoplasia. Clin Cancer Res 2007;13:929–38.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Bera TK,
    2. Iavarone C,
    3. Kumar V,
    4. Lee S,
    5. Lee B,
    6. Pastan I
    . MRP9, an unusual truncated member of the ABC transporter superfamily, is highly expressed in breast cancer. Proc Natl Acad Sci U S A 2002;99:6997–7002.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Bunting KD
    . ABC transporters as phenotypic markers and functional regulators of stem cells. Stem Cells 2002;20:11–20.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Chen Y,
    2. Takita J,
    3. Choi YL,
    4. et al
    . Oncogenic mutations of ALK kinase in neuroblastoma. Nature 2008;455:971–4.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kumar U,
    2. Grigorakis SI,
    3. Watt HL,
    4. et al
    . Somatostatin receptors in primary human breast cancer: quantitative analysis of mRNA for subtypes 1–5 and correlation with receptor protein expression and tumor pathology. Breast Cancer Res Treat 2005;92:175–86.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Negishi M,
    2. Oinuma I,
    3. Katoh H
    . Plexins: axon guidance and signal transduction. Cell Mol Life Sci 2005;62:1363–71.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Gandellini P,
    2. Folini M,
    3. Longoni N,
    4. et al
    . miR-205 exerts tumor-suppressive functions in human prostate through down-regulation of protein kinase Cϵ. Cancer Res 2009;69:2287–95.
    OpenUrlAbstract/FREE Full Text
    1. Iorio MV,
    2. Casalini P,
    3. Piovan C,
    4. et al
    . microRNA-205 regulates HER3 in human breast cancer. Cancer Res 2009;69:2195–200.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Lebanony D,
    2. Benjamin H,
    3. Gilad S,
    4. et al
    . Diagnostic assay based on hsa-miR-205 expression distinguishes squamous from nonsquamous non-small-cell lung carcinoma. J Clin Oncol 2009;27:2030–7.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Bray I,
    2. Bryan K,
    3. Prenter S,
    4. et al
    . Widespread dysregulation of miRNAs by MYCN amplification and chromosomal imbalances in neuroblastoma: association of miRNA expression with survival. PLoS One 2009;4:e7850.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Fischer M,
    2. Bauer T,
    3. Oberthur A,
    4. et al
    . Integrated genomic profiling identifies two distinct molecular subtypes with divergent outcome in neuroblastoma with loss of chromosome 11q. Oncogene 2009;29:865–75.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Oberthuer A,
    2. Berthold F,
    3. Warnat P,
    4. et al
    . Customized oligonucleotide microarray gene expression-based classification of neuroblastoma patients outperforms current clinical risk stratification. J Clin Oncol 2006;24:5070–8.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Vermeulen J,
    2. De Preter K,
    3. Naranjo A,
    4. et al
    . Predicting outcomes for children with neuroblastoma using a multigene-expression signature: a retrospective SIOPEN/COG/GPOH study. Lancet Oncol 2009;10:663–71.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Mestdagh P,
    2. Fredlund E,
    3. Pattyn F,
    4. et al
    . MYCN/c-MYC-induced microRNAs repress coding gene networks associated with poor outcome in MYCN/c-MYC-activated tumors. Oncogene 2009.
  40. ↵
    1. Janoueix-Lerosey I,
    2. Schleiermacher G,
    3. Michels E,
    4. et al
    . Overall genomic pattern is a predictor of outcome in neuroblastoma. J Clin Oncol 2009;27:1026–33.
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Schleiermacher G,
    2. Michon J,
    3. Huon I,
    4. et al
    . Chromosomal CGH identifies patients with a higher risk of relapse in neuroblastoma without MYCN amplification. Br J Cancer 2007;97:238–46.
    OpenUrlCrossRefPubMed
View Abstract
PreviousNext
Back to top
Clinical Cancer Research: 16 (11)
June 2010
Volume 16, Issue 11
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Clinical Cancer Research article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Chromosomal and MicroRNA Expression Patterns Reveal Biologically Distinct Subgroups of 11q− Neuroblastoma
(Your Name) has forwarded a page to you from Clinical Cancer Research
(Your Name) thought you would be interested in this article in Clinical Cancer Research.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Chromosomal and MicroRNA Expression Patterns Reveal Biologically Distinct Subgroups of 11q− Neuroblastoma
Patrick G. Buckley, Leah Alcock, Kenneth Bryan, Isabella Bray, Johannes H. Schulte, Alexander Schramm, Angelika Eggert, Pieter Mestdagh, Katleen De Preter, Jo Vandesompele, Frank Speleman and Raymond L. Stallings
Clin Cancer Res June 1 2010 (16) (11) 2971-2978; DOI: 10.1158/1078-0432.CCR-09-3215

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Chromosomal and MicroRNA Expression Patterns Reveal Biologically Distinct Subgroups of 11q− Neuroblastoma
Patrick G. Buckley, Leah Alcock, Kenneth Bryan, Isabella Bray, Johannes H. Schulte, Alexander Schramm, Angelika Eggert, Pieter Mestdagh, Katleen De Preter, Jo Vandesompele, Frank Speleman and Raymond L. Stallings
Clin Cancer Res June 1 2010 (16) (11) 2971-2978; DOI: 10.1158/1078-0432.CCR-09-3215
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Materials and Methods
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Acknowledgments
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • Sirt7 Promotes Colorectal Cancer Tumorigenesis
  • Contact Guidance Controls T-cell Migration in PDAC
  • MET in Papillary RCC
Show more Human Cancer Biology
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • CCR Focus Archive
  • Meeting Abstracts

Info for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About Clinical Cancer Research

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Clinical Cancer Research
eISSN: 1557-3265
ISSN: 1078-0432

Advertisement