High-resolution detection and mapping of genomic DNA alterations in neuroblastoma

17q Gain in Neuroblastoma: A Review of Clinical and Biological Implications

Neuroblastoma (NB) is the most frequent extracranial solid childhood tumor. Despite advances in the understanding and treatment of this disease, the prognosis in cases of high-risk NB is still poor. 17q gain has been shown to be the most frequent genomic alteration in NB. However, the significance of this remains unclear because of its high frequency and association with other genetic modifications, particularly segmental chromosomal aberrations, 1p and 11q deletions, and MYCN amplification, all of which are also associated with a poor clinical prognosis. This work reviewed the evidence on the clinical and biological significance of 17q gain. It strongly supports the significance of 17q gain in the development of NB and its importance as a clinically relevant marker. However, it is crucial to distinguish between whole and partial chromosome 17q gains. The most important breakpoints appear to be at 17q12 and 17q21. The former distinguishes between whole and partial chromosome 17q gain; the latter is a site of IGF2BP1 and NME1 genes that appear to be the main oncogenes responsible for the functional effects of 17q gain.

KDM6B promotes activation of the oncogenic CDK4/6-pRB-E2F pathway by maintaining enhancer activity in MYCN-amplified neuroblastoma

The H3K27me2/me3 histone demethylase KDM6B is essential to neuroblastoma cell survival. However, the mechanism of KDM6B action remains poorly defined. We demonstrate that inhibition of KDM6B activity 1) reduces the chromatin accessibility of E2F target genes and MYCN, 2) selectively leads to an increase of H3K27me3 but a decrease of the enhancer mark H3K4me1 at the CTCF and BORIS binding sites, which may, consequently, disrupt the long-range chromatin interaction of MYCN and E2F target genes, and 3) phenocopies the transcriptome induced by the specific CDK4/6 inhibitor palbociclib. Overexpression of CDK4/6 or Rb1 knockout confers neuroblastoma cell resistance to both palbociclib and the KDM6 inhibitor GSK-J4. These data indicate that KDM6B promotes an oncogenic CDK4/6-pRB-E2F pathway in neuroblastoma cells via H3K27me3-dependent enhancer-promoter interactions, providing a rationale to target KDM6B for high-risk neuroblastoma.

The histone demethylase KDM6B is reported to be essential for neuroblastoma cell survival. Here the authors show that KDM6B regulates CDK4/6-pRB-E2F pathway through H3K27me3-dependent enhancer-promoter interactions in neuroblastoma.

GAS7 Deficiency Promotes Metastasis in MYCN-Driven Neuroblastoma

One of the greatest barriers to curative treatment of neuroblastoma is its frequent metastatic outgrowth prior to diagnosis, especially in cases driven by amplification of the MYCN oncogene. However, only a limited number of regulatory proteins that contribute to this complex MYCN-mediated process have been elucidated. Here we show that the growth arrest-specific 7 (GAS7) gene, located at chromosome band 17p13.1, is preferentially deleted in high-risk MYCN-driven neuroblastoma. GAS7 expression was also suppressed in MYCN-amplified neuroblastoma lacking 17p deletion. GAS7 deficiency led to accelerated metastasis in both zebrafish and mammalian models of neuroblastoma with overexpression or amplification of MYCN. Analysis of expression profiles and the ultrastructure of zebrafish neuroblastoma tumors with MYCN overexpression identified that GAS7 deficiency led to (i) downregulation of genes involved in cell-cell interaction, (ii) loss of contact among tumor cells as critical determinants of accelerated metastasis, and (iii) increased levels of MYCN protein. These results provide the first genetic evidence that GAS7 depletion is a critical early step in the cascade of events culminating in neuroblastoma metastasis in the context of MYCN overexpression. SIGNIFICANCE: Heterozygous deletion or MYCN-mediated repression of GAS7 in neuroblastoma releases an important brake on tumor cell dispersion and migration to distant sites, providing a novel mechanism underlying tumor metastasis in MYCN-driven neuroblastoma.See related commentary by Menard, p. 2815.

Non-canonical roles of canonical telomere binding proteins in cancers

Reactivation of telomerase is a major hallmark observed in 90% of all cancers. Yet paradoxically, enhanced telomerase activity does not correlate with telomere length and cancers often possess short telomeres; suggestive of supplementary non-canonical roles that telomerase might play in the development of cancer. Moreover, studies have shown that aberrant expression of shelterin proteins coupled with their release from shortening telomeres can further promote cancer by mechanisms independent of their telomeric role. While targeting telomerase activity appears to be an attractive therapeutic option, this approach has failed in clinical trials due to undesirable cytotoxic effects on stem cells. To circumvent this concern, an alternative strategy could be to target the molecules involved in the non-canonical functions of telomeric proteins. In this review, we will focus on emerging evidence that has demonstrated the non-canonical roles of telomeric proteins and their impact on tumorigenesis. Furthermore, we aim to address current knowledge gaps in telomeric protein functions and propose future research approaches that can be undertaken to achieve this.

Large 1p36 Deletions Affecting Arid1a Locus Facilitate Mycn-Driven Oncogenesis in Neuroblastoma

Loss of heterozygosity (LOH) at 1p36 occurs in multiple cancers, including neuroblastoma (NBL). MYCN amplification and 1p36 deletions tightly correlate with markers of tumor aggressiveness in NBL. Although distal 1p36 losses associate with single-copy MYCN tumors, larger deletions correlate with MYCN amplification, indicating two tumor suppressor regions in 1p36, only one of which facilitates MYCN oncogenesis. To better define this region, we genome-edited the syntenic 1p36 locus in primary mouse neural crest cells (NCCs), a putative NBL cell of origin. In in vitro cell transformation assays, we show that Chd5 loss confers most of the MYCN-independent tumor suppressor effects of 1p36 LOH. In contrast, MYCN-driven tumorigenesis selects for NCCs with Arid1a deletions from a pool of NCCs with randomly sized 1p36 deletions, establishing Arid1a as the MYCN-associated tumor suppressor. Our findings reveal that Arid1a loss collaborates with oncogenic MYCN and better define the tumor suppressor functions of 1p36 LOH in NBL.

Somatic structural variation targets neurodevelopmental genes and identifies SHANK2 as a tumor suppressor in neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that accounts for 12% of childhood cancer deaths. Like many childhood cancers, neuroblastoma shows a relative paucity of somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) compared to adult cancers. Here, we assessed the contribution of somatic structural variation (SV) in neuroblastoma using a combination of whole-genome sequencing (WGS) of tumor-normal pairs (n = 135) and single-nucleotide polymorphism (SNP) genotyping of primary tumors (n = 914). Our study design allowed for orthogonal validation and replication across platforms. SV frequency, type, and localization varied significantly among high-risk tumors. MYCN nonamplified high-risk tumors harbored an increased SV burden overall, including a significant excess of tandem duplication events across the genome. Genes disrupted by SV breakpoints were enriched in neuronal lineages and associated with phenotypes such as autism spectrum disorder (ASD). The postsynaptic adapter protein-coding gene, SHANK2, located on Chromosome 11q13, was disrupted by SVs in 14% of MYCN nonamplified high-risk tumors based on WGS and 10% in the SNP array cohort. Expression of SHANK2 was low across human-derived neuroblastoma cell lines and high-risk neuroblastoma tumors. Forced expression of SHANK2 in neuroblastoma cells resulted in significant growth inhibition (P = 2.6 × 10-2 to 3.4 × 10-5) and accelerated neuronal differentiation following treatment with all-trans retinoic acid (P = 3.1 × 10-13 to 2.4 × 10-30). These data further define the complex landscape of somatic structural variation in neuroblastoma and suggest that events leading to deregulation of neurodevelopmental processes, such as inactivation of SHANK2, are key mediators of tumorigenesis in this childhood cancer.

Central neurocytoma: SNP array analyses, subtel FISH, and review of the literature

The central neurocytoma (CN) is a rare brain tumor with a frequency of 0.1-0.5% of all brain tumors. According to the World Health Organization classification, it is a benign grade II tumor with good prognosis. However, some CN occur as histologically "atypical" variant, combined with increasing proliferation and poor clinical outcome. Detailed genetic knowledge could be helpful to characterize a potential atypical behavior in CN. Only few publications on genetics of CN exist in the literature. Therefore, we performed cytogenetic analysis of an intraventricular neurocytoma WHO grade II in a 39-year-old male patient by use of genome-wide high-density single nucleotide polymorphism array (SNP array) and subtelomere FISH. Applying these techniques, we could detect known chromosomal aberrations and identified six not previously described chromosomal aberrations, gains of 1p36.33-p36.31, 2q37.1-q37.3, 6q27, 12p13.33-p13.31, 20q13.31-q13.33, and loss of 19p13.3-p12. Our case report contributes to the genetic knowledge about CN and to increased understanding of "typical" and "atypical" variants.

Clinical Importance of Myc Family Oncogene Aberrations in Epithelial Ovarian Cancer

Background

The Myc oncogene family has been implicated in many human malignancies and is often associated with particularly aggressive disease, suggesting Myc as an attractive prognostic marker and therapeutic target. However, for epithelial ovarian cancer (EOC), there is little consensus on the incidence and clinical relevance of Myc aberrations. Here we comprehensively investigated alterations in gene copy number, expression, and activity for Myc and evaluated their clinical significance in EOC.

Methods

To address inconsistencies in the literature regarding the definition of copy number variations, we developed a novel approach using quantitative polymerase chain reaction (qPCR) coupled with a statistical algorithm to estimate objective thresholds for detecting Myc gain/amplification in large cohorts of serous (n = 150) and endometrioid (n = 80) EOC. MYC, MYCN, and MYCL1 mRNA expression and Myc activity score for each case were examined by qPCR. Kaplan–Meier and Cox-regression analyses were conducted to assess clinical significance of Myc aberrations.

Results

Using a large panel of cancer cell lines (n = 34), we validated the statistical algorithm for determining clear thresholds for Myc gain/amplification. MYC was the most predominantly amplified of the Myc oncogene family members, and high MYC mRNA expression levels were associated with amplification in EOC. However, there was no association between prognosis and increased copy number or gene expression of MYC/MYCN/MYCL1 or with a pan-Myc transcriptional activity score, in EOC, although MYC amplification was associated with late stage and high grade in endometrioid EOC.

Conclusion

A systematic and comprehensive analysis of Myc genes, transcripts, and activity levels using qPCR revealed that although such aberrations commonly occur in EOC, overall they have limited impact on outcome, suggesting that the biological relevance of Myc oncogene family members is limited to certain subsets of this disease.

Clinical evaluation of integrated panel testing by next-generation sequencing for somatic mutations in neuroblastomas with MYCN unamplification

Neuroblastomas (NBs) exhibit heterogeneity and show clinically significant prognosis classified by genetic alterations. Among prognostic genes or genome factors, MYCN amplification (MNA) is the most established genomic marker of poor prognosis in patients with NB. However, the prognostic classification of more than 60% of patients without MNA has yet to be clarified. In this study, the application of target next-generation sequencing (NGS) was extended on the basis of a comprehensive panel of regions where copy number variations (CNVs) or point mutations occurred to improve the prognostic evaluation of these patients and obtain the sequence of 33 patients without MNA. A mean coverage depth of 887× was determined in the target regions in all of the samples, and the mapped read percentage was more than 99%. Somatic mutations in patients without MNA could be precisely defined on the basis of these findings, and 17 unique somatic aberrations, including 14 genes, were identified in 11 patients. Among these variations, most were CNVs with a number of 13. The 3-year event-free survival (EFS) of CNV(−) patients was 60.0% compared with the EFS (16.7%) of CNV(+) patients (P = 0.015, HR = 0.1344, 95%, CI = 0.027 to 0.678). CNVs were also associated with unfavorable histological characteristics (P = 0.003) and likely to occur in stage 4 (P = 0.041). These results might further indicate the role of CNVs in NB chemotherapy resistance (P = 0.059) and show CNVs as a therapeutic target. In multivariate analysis, the presence of CNVs was a clinically negative prognostic marker that impaired the outcome of patients without MNA and associated with poor prognosis in this tumor subset. Comprehensive genetic/genomic profiling instead of focusing on single genetic marker should be performed through in-depth NGS that could reveal prognostic information, improve NB target therapy, and provide a basis for investigations on NB pathogenesis.

Molecular mechanisms and therapeutic targets in neuroblastoma

Neuroblastoma is the most common extracranical tumor of childhood and the most deadly tumor of infancy. It is characterized by early age onset and high frequencies of metastatic disease but also the capacity to spontaneously regress. Despite intensive therapy, the survival for patients with high-risk neuroblastoma and those with recurrent or relapsed disease is low. Hence, there is an urgent need to develop new therapies for these patient groups. The molecular pathogenesis based on high-throughput omics technologies of neuroblastoma is beginning to be resolved which have given the opportunity to develop personalized therapies for high-risk patients. Here we discuss the potential of developing targeted therapies against aberrantly expressed molecules detected in sub-populations of neuroblastoma patients and how these selected targets can be drugged in order to overcome treatment resistance, improve survival and quality of life for these patients and also the possibilities to transfer preclinical research into clinical testing.

Role of microRNAs in epigenetic silencing of the CHD5 tumor suppressor gene in neuroblastomas

Neuroblastoma (NB), a tumor of the sympathetic nervous system, is the most common extracranial solid tumor of childhood. We and others have identified distinct patterns of genomic change that underlie diverse clinical behaviors, from spontaneous regression to relentless progression. We first identified CHD5 as a tumor suppressor gene that is frequently deleted in NBs. Mutation of the remaining CHD5 allele is rare in these tumors, yet expression is very low or absent, so expression is likely regulated by epigenetic mechanisms. In order to understand the potential role of miRNA regulation of CHD5 protein expression in NBs, we examined all miRNAs that are predicted to target the 3′-UTR using miRanda, TargetScan and other algorithms. We identified 18 miRNAs that were predicted by 2 or more programs: miR-204, -211, -216b, -17, -19ab, -20ab, -93, -106ab, -130ab, -301ab, -454, -519d, -3666. We then performed transient transfections in two NB cell lines, NLF (MYCN amplified) and SY5Y (MYCN non-amplified), with the reporter plasmid and miRNA mimic, as well as appropriate controls. We found seven miRNAs that significantly downregulated CHD5 expression in NB: miR-211, 17, -93, -20b, -106b, -204, and -3666. Interestingly, MYCN upregulates several of the candidates we identified: miR-17, -93, -106b & -20b. This suggests that miRNAs driven by MYCN and other genes represent a potential epigenetic mechanism to regulate CHD5 expression.

Identification of Aberrant Chromosomal Regions in Human Breast Cancer Using Gene Expression Data and Related Gene Information

BACKGROUND

Chromosomal instability is a hallmark of cancer. Chromosomal imbalances, like amplifications and deletions, influence the transcriptional activity of genes. These imbalances affect not only the expression of genes in the aberrant chromosomal regions, but also that of related genes, and may be relevant to the cancer status.

MATERIAL AND METHODS

Here, we used the 7 publicly available microarray studies in breast cancer tissues and propose a general and unsupervised method by using the gene expression data and related gene information to systematically identify aberrant chromosomal regions. This method aimed to identify the chromosomal regions where the genes and their related genes both show consistent changes in the expression levels. Such patterns have been reported to be associated with the chromosomal aberrations and may be used in cancer diagnosis.

RESULTS

We compared 488 tumor and 222 normal samples from 7 microarray-based human breast cancer studies and detected the amplifications of 8q11.21, 14q32.11, 4q21.23, 18q11.2, Xq28, and the deletions of 3p24.1, 10q23.2 (BSCG1), 20p11.21, 9q21.13, and 1q41, which may be involved in the novel mechanisms of tumorigenesis. In addition, several known pathogenic genes, transcription factors (TFs), and microRNAs (miRNAs) associated with breast cancer were found.

CONCLUSIONS

This approach can be applied to other microarray studies, which provide a new and useful method for exploring chromosome structural variations in different types of diseases.

The tumour suppressor CHD5 forms a NuRD-type chromatin remodelling complex

Eukaryotic gene expression is developmentally regulated, in part by chromatin remodelling, and its dysregulation has been linked to cancer. CHD5 (chromodomain helicase DNA-binding protein 5) is a tumour suppressor gene (TSG) that maps to a region of consistent deletion on 1p36.31 in neuroblastomas (NBs) and other tumour types. CHD5 encodes a protein with chromatin remodelling, helicase and DNA-binding motifs that is preferentially expressed in neural and testicular tissues. CHD5 is highly homologous to CHD3 and CHD4, which are the core subunits of nucleosome remodelling and deacetylation (NuRD) complexes. To determine if CHD5 forms a similar complex, we performed studies on nuclear extracts from NBLS, SY5Y (both with endogenous CHD5 expression), NLF (CHD5 null) and NLF cells stably transfected with CHD5 cDNA (wild-type and V5-histidine-tagged). Immunoprecipitation (IP) was performed with either CHD5 antibody or antibody to V5/histidine-tagged protein. We identified NuRD components both by GST-FOG1 (Friend Of GATA1) pull-down and by IP. We also performed MS/MS analysis to confirm the presence of CHD5 or other protein components of the NuRD complex, as well as to identify other novel proteins. CHD5 was clearly associated with all canonical NuRD components, including metastasis-associated protein (MTA)1/2, GATA zinc finger domain containing 2A (GATAD2A), histone deacetylase (HDAC)1/2, retinoblastoma-binding protein (RBBP)4/7 and methyl DNA-binding domain protein (MBD)2/3, as determined by Western blotting and MS/MS. Our data suggest CHD5 forms a NuRD complex similar to CHD4. However, CHD5-NuRD may also have unique protein associations that confer functional specificity and may contribute to normal development and to tumour suppression in NB and other cancers.

Clinical characteristics and outcome of patients with neuroblastoma presenting genomic amplification of loci other than MYCN

Background

Somatically acquired genomic alterations with MYCN amplification (MNA) are key features of neuroblastoma (NB), the most common extra-cranial malignant tumour of childhood. Little is known about the frequency, clinical characteristics and outcome of NBs harbouring genomic amplification(s) distinct from MYCN.

Methods

Genomic profiles of 1100 NBs from French centres studied by array-CGH were re-examined specifically to identify regional amplifications. Patients were included if amplifications distinct from the MYCN locus were seen. A subset of NBs treated at Institut Curie and harbouring MNA as determined by array-CGH without other amplification was also studied. Clinical and histology data were retrospectively collected.

Results

In total, 56 patients were included and categorised into 3 groups. Group 1 (n = 8) presented regional amplification(s) without MNA. Locus 12q13-14 was a recurrent amplified region (4/8 cases). This group was heterogeneous in terms of INSS stages, primary localisations and histology, with atypical clinical features. Group 2 (n = 26) had MNA as well as other regional amplifications. These patients shared clinical features of those of a group of NBs MYCN amplified (Group 3, n = 22). Overall survival for group 1 was better than that of groups 2 and 3 (5 year OS: 87.5%±11% vs 34.9%±7%, log-rank p<0.05).

Conclusion

NBs harbouring regional amplification(s) without MNA are rare and seem to show atypical features in clinical presentation and genomic profile. Further high resolution genetic explorations are justified in this heterogeneous group, especially when considering these alterations as predictive markers for targeted therapy.

Therapeutic targets for neuroblastomas

INTRODUCTION

Neuroblastoma (NB) is the most common and deadly solid tumor in children. Despite recent improvements, the long-term outlook for high-risk NB is still < 50%. Further, there is considerable short- and long-term toxicity. More effective, less toxic therapy is needed, and the development of targeted therapies offers great promise.

AREAS COVERED

Relevant literature was reviewed to identify current and future therapeutic targets that are critical to malignant transformation and progression of NB. The potential or actual NB therapeutic targets are classified into four categories: i) genes activated by amplification, mutation, translocation or autocrine overexpression; ii) genes inactivated by deletion, mutation or epigenetic silencing; iii) membrane-associated genes expressed on most NBs but few other tissues; or iv) common target genes relevant to NB as well as other tumors.

EXPERT OPINION

Therapeutic approaches have been developed to some of these targets, but many remain untargeted at the present time. It is unlikely that single targeted agents will be sufficient for long-term cure, at least for high-risk NBs. The challenge will be how to integrate targeted agents with each other and with conventional therapy to enhance their efficacy, while simultaneously reducing systemic toxicity.

Preclinical evaluation of the PARP inhibitor, olaparib, in combination with cytotoxic chemotherapy in pediatric solid tumors

BACKGROUND

Poly(ADP-ribose) polymerase (PARP) signals DNA damage and facilitates DNA repair. PARP inhibitors are being evaluated in cancers with defective DNA repair mechanisms or in combination with cytotoxic therapy or radiation. We evaluated the PARP inhibitor, olaparib, in combination with chemotherapy using in vitro and in vivo pediatric solid tumor models.

PROCEDURE

The IC50 of olaparib alone and in combination with cytotoxic agents was determined in 10 pediatric solid tumor cell lines. Synergy was assessed using the combination index of Chou-Talalay. Olaparib alone and in combination with topotecan/cyclophosphamide was evaluated in xenograft models of Ewing sarcoma (RD-ES) and neuroblastoma (NGP). PAR activity was evaluated in cell lines and tumor lysates.

RESULTS

Olaparib induced growth inhibition, median (range) IC50 = 3.6 (1-33.8) µM, and inhibited PAR activity in pediatric solid tumor cell lines. The addition of olaparib to DNA damaging agents resulted in additive to synergistic interactions. In RD-ES and NGP xenografts, olaparib inhibited PAR activity by 88-100% as a single agent and 100% when administered with cyclophosphamide/topotecan. Although the addition of olaparib did not antagonize the activity of cyclophosphamide/topotecan, clear evidence of synergy could not be demonstrated.

CONCLUSIONS

In pediatric solid tumor cell lines, clinically achievable concentrations of single agent olaparib caused growth inhibition. Although the in vitro data demonstrated synergistic efficacy of olaparib when added to the camptothecins and alkylating agents, synergy was not discernible in vivo. Clinical trials of PARP inhibitors in combination DNA damaging agents are necessary to establish the role of PARP inhibitors in childhood cancer.

Recurrent pre-existing and acquired DNA copy number alterations, including focal TERT gains, in neuroblastoma central nervous system metastases

Stage 4 neuroblastomas have a high rate of local and metastatic relapse and associated disease mortality. The central nervous system (CNS) is currently one of the most common isolated relapse sites, yet the genomic alterations that contribute to these metastases are unknown. This study sought to identify recurrent DNA copy number alterations (CNAs) and target genes relating to neuroblastoma CNS metastases by studying 19 pre-CNS primary tumors and 27 CNS metastases, including 12 matched pairs. SNP microarray analyses revealed that MYCN amplified (MYCNA) tumors had recurrent CNAs different from non-MYCNA cohorts. Several CNAs known to be prevalent among primary neuroblastomas occurred more frequently in CNS metastases, including 4p-, 7q+, 12q+, and 19q- in non-MYCNA metastases, and 9p- and 14q- irrespective of MYCNA status. In addition, novel CNS metastases-related CNAs included 18q22.1 gains in non-MYCNA pre-CNS primaries and 5p15.33 gains and 15q26.1→tel losses in non-MYCNA CNS metastases. Based on minimal common regions, gene expression, and biological properties, TERT (5p), NR2F2 (15q), ALDH1A3 (15q), CDKN2A (9p), and possibly CDH7 and CDH19 (18q) were candidate genes associated with the CNS metastatic process. Notably, the 5p15 minimal common region contained only TERT, and non-MYCNA CNS metastases with focal 5p15 gains had increased TERT expression, similar to MYCNA tumors. These findings suggest that a specific genomic lesion (18q22.1 gain) predisposes to CNS metastases and that distinct lesions are recurrently acquired during metastatic progression. Among the acquired lesions, increased TERT copy number and expression appears likely to function in lieu of MYCNA to promote CNS metastasis.

Dual CDK4/CDK6 inhibition induces cell-cycle arrest and senescence in neuroblastoma

PURPOSE

Neuroblastoma is a pediatric cancer that continues to exact significant morbidity and mortality. Recently, a number of cell-cycle proteins, particularly those within the Cyclin D/CDK4/CDK6/RB network, have been shown to exert oncogenic roles in neuroblastoma, suggesting that their therapeutic exploitation might improve patient outcomes.

EXPERIMENTAL PROCEDURES

We evaluated the effect of dual CDK4/CDK6 inhibition on neuroblastoma viability using LEE011 (Novartis Oncology), a highly specific CDK4/6 inhibitor.

RESULTS

Treatment with LEE011 significantly reduced proliferation in 12 of 17 human neuroblastoma-derived cell lines by inducing cytostasis at nanomolar concentrations (mean IC50 = 307 ± 68 nmol/L in sensitive lines). LEE011 caused cell-cycle arrest and cellular senescence that was attributed to dose-dependent decreases in phosphorylated RB and FOXM1, respectively. In addition, responsiveness of neuroblastoma xenografts to LEE011 translated to the in vivo setting in that there was a direct correlation of in vitro IC50 values with degree of subcutaneous xenograft growth delay. Although our data indicate that neuroblastomas sensitive to LEE011 were more likely to contain genomic amplification of MYCN (P = 0.01), the identification of additional clinically accessible biomarkers is of high importance.

CONCLUSIONS

Taken together, our data show that LEE011 is active in a large subset of neuroblastoma cell line and xenograft models, and supports the clinical development of this CDK4/6 inhibitor as a therapy for patients with this disease. Clin Cancer Res; 19(22); 6173-82. ©2013 AACR.

The stress protein BAG3 stabilizes Mcl-1 protein and promotes survival of cancer cells and resistance to antagonist ABT-737

Members of the Bcl-2 family of proteins are important inhibitors of apoptosis in human cancer and are targets for novel anticancer agents such as the Bcl-2 antagonists, ABT-263 (Navitoclax), and its analog ABT-737. Unlike Bcl-2, Mcl-1 is not antagonized by ABT-263 or ABT-737 and is considered to be a major factor in resistance. Also, Mcl-1 exhibits differential regulation when compared with other Bcl-2 family members and is a target for anticancer drug discovery. Here, we demonstrate that BAG3, an Hsp70 co-chaperone, protects Mcl-1 from proteasomal degradation, thereby promoting its antiapoptotic activity. Using neuroblastoma cell lines, with a defined Bcl-2 family dependence, we found that BAG3 expression correlated with Mcl-1 dependence and ABT-737 resistance. RNA silencing of BAG3 led to a marked reduction in Mcl-1 protein levels and overcame ABT-737 resistance in Mcl-1-dependent cells. In ABT-737-resistant cells, Mcl-1 co-immunoprecipitated with BAG3, and loss of Mcl-1 after BAG3 silencing was prevented by proteasome inhibition. BAG3 and Mcl-1 were co-expressed in a panel of diverse cancer cell lines resistant to ABT-737. Silencing BAG3 reduced Mcl-1 protein levels and overcame ABT-737 resistance in several of the cell lines, including triple-negative breast cancer (MDA-MB231) and androgen receptor-negative prostate cancer (PC3) cells. These studies identify BAG3-mediated Mcl-1 stabilization as a potential target for cancer drug discovery.

COPS: a sensitive and accurate tool for detecting somatic Copy Number Alterations using short-read sequence data from paired samples

Copy Number Alterations (CNAs) such as deletions and duplications; compose a larger percentage of genetic variations than single nucleotide polymorphisms or other structural variations in cancer genomes that undergo major chromosomal re-arrangements. It is, therefore, imperative to identify cancer-specific somatic copy number alterations (SCNAs), with respect to matched normal tissue, in order to understand their association with the disease. We have devised an accurate, sensitive, and easy-to-use tool, COPS, COpy number using Paired Samples, for detecting SCNAs. We rigorously tested the performance of COPS using short sequence simulated reads at various sizes and coverage of SCNAs, read depths, read lengths and also with real tumor:normal paired samples. We found COPS to perform better in comparison to other known SCNA detection tools for all evaluated parameters, namely, sensitivity (detection of true positives), specificity (detection of false positives) and size accuracy. COPS performed well for sequencing reads of all lengths when used with most upstream read alignment tools. Additionally, by incorporating a downstream boundary segmentation detection tool, the accuracy of SCNA boundaries was further improved. Here, we report an accurate, sensitive and easy to use tool in detecting cancer-specific SCNAs using short-read sequence data. In addition to cancer, COPS can be used for any disease as long as sequence reads from both disease and normal samples from the same individual are available. An added boundary segmentation detection module makes COPS detected SCNA boundaries more specific for the samples studied. COPS is available at ftp://115.119.160.213 with username “cops” and password “cops”.

Notch pathway activation induces neuroblastoma tumor cell growth arrest

BACKGROUND

Notch pathway signaling has critical roles in differentiation, proliferation, and survival, and has oncogenic or tumor suppressor effects in a variety of malignancies. The goal of this study was to evaluate the effects of Notch activation on human neuroblastoma cells.

PROCEDURE

Quantitative RT-PCR, immunoblots, and immunohistochemistry were used to determine the expression of Notch receptors (Notch1-4), cleaved Notch1 (ICN1), and downstream targets (HES1-5) in human neuroblastoma cell lines and patient tumor samples. Notch pathway signaling was induced using intracellular Notch (ICN1-3) and HES gene constructs or direct culture on Notch ligands. Quantitative methylation-specific PCR was used to quantify methylation of the HES gene promoters, and the effects of treatment with decitabine were measured.

RESULTS

Neuroblastoma cells express varying levels of Notch receptors and low levels of HES genes at baseline. However, no endogenous activation of the Notch pathway was detected in neuroblastoma cell lines or patient tumor samples. Expression of activated Notch intracellular domains and HES gene products led to growth arrest. The HES2 and HES5 gene promoters were found to be heavily methylated in most neuroblastoma lines, and HES gene expression could be induced through treatment with decitabine.

CONCLUSIONS

We report that neuroblastoma cell lines express multiple Notch receptors, which are inactive at baseline. Activation of the Notch pathway via ligand binding consistently resulted in growth arrest. HES gene expression appears to be regulated epigenetically and could be induced with decitabine. These findings support a tumor suppressor role for Notch signaling in neuroblastoma.

Progress in treatment and risk stratification of neuroblastoma: impact on future clinical and basic research

Close international collaboration between pediatric oncologists has led to marked improvements in the cure of patients, seen as a long-term overall survival rate of about 80%. Despite this progress, neuroblastoma remains a challenging disease for both clinicians and researchers. Major clinical problems include lack of acceptable cure rates in high-risk neuroblastoma and potential overtreatment of subsets of patients at low and intermediate risk of the disease. Many years of intensive international cooperation have recently led to a promising joint effort to further improve risk classification for treatment stratification, the new International Neuroblastoma Risk Group Classification System. This approach will facilitate comparison of the results of clinical trials performed by different international collaborative groups. This, in turn, should accelerate refinement of risk stratification and thereby aid selection of appropriate therapies for individual patients. To be able to identify new therapeutic modalities, it will be necessary to elucidate the pathogenesis of the different subtypes of neuroblastoma. Basic and translational research have provided new tools for molecular characterization of blood and tumor samples including high-throughput technologies for analysis of DNA, mRNAs, microRNAs and other non-coding RNAs, as well as proteins and epigenetic markers. Most of these studies are array-based in design. In neuroblastoma research they aim to refine risk group stratification through incorporation of molecular tumor fingerprints and also to enable personalized treatment modalities by describing the underlying pathogenesis and aberrant signaling pathways in individual tumors. To make optimal use of these new technologies for the benefit of the patient, it is crucial to have a systematic and detailed documentation of both clinical and molecular data from diagnosis through treatment to follow-up. Close collaboration between clinicians and basic scientists will provide access to combined clinical and molecular data sets and will create more efficient steps in response to the remaining treatment challenges. This review describes the current efforts and trends in neuroblastoma research from a clinical perspective in order to highlight the urgent clinical problems we must address together with basic researchers.

Isolation of disseminated neuroblastoma cells from bone marrow aspirates for pretreatment risk assessment by array comparative genomic hybridization

In neuroblastoma, tumor biopsies are used for prognostic evaluation and risk assessment by molecular genetic analyses such as fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (array CGH). Analysis of primary tumors by array CGH can be hampered by the lack of sufficient tumor cells due to small biopsy size or availability of invaded bone marrow only. Given the importance of accurate assessment of genetic alterations in the diagnostic work-up of patients with neuroblastoma, we evaluated the possibility to analyze bone marrow metastases in patients with disseminated disease. Disseminated neuroblastoma cells were isolated from bone marrow aspirates by using either laser capture microdissection (LCM) or magnetic activated cell sorting (MACS). The array CGH profiles of these isolated metastases were compared to array CGH profiles and/or FISH data of the corresponding primary tumor. Here, we show that the major recurrent DNA copy number alterations detected in primary neuroblastoma tumors (i.e., 1p, 3p and 11q deletion, 17q gain and MYCN amplification) can be detected, with high sensitivity and specificity, in the disseminated neuroblastoma cells isolated from the bone marrow aspirates, using an array platform with high coverage for these regions. Moreover, we demonstrate that for archived material, for example, for retrospective studies, LCM is the method of choice, while for fresh bone marrow aspirates, acquired at the time of diagnosis, MACS is superior.

Global genomic and RNA profiles for novel risk stratification of neuroblastoma

Neuroblastoma is one of the most common solid tumors in children. Its clinical behavior ranges widely from spontaneous regression to life-threatening aggressive growth. The molecular etiology of neuroblastoma is still enigmatic and the overall cure rate of advanced disease is still very poor. Recent microarray-based technology provided us with important information such as comprehensive genomic alterations and gene expression profiles to help us understand the molecular characteristics of each tumor in detail. Several retrospective studies have revealed that these signatures are strongly correlated with patient prognoses and led to the construction of new risk stratification systems, some of which are considered for evaluation in upcoming clinical studies in a prospective way. Large-scale analyses using a variety of genetic tools also discovered a major familial neuroblastoma predisposition gene ALK, as well as new candidate susceptibility genes at 6q22 and 2q35 for sporadic neuroblastoma. Of note, ALK is mutated in 6-9% of sporadic cases, and is either amplified or constitutively activated through mutations mainly within the kinase domain, promoting the possibility of new therapeutic strategies using ALK inhibitors. Additional candidates for outcome predictors such as the methylation phenotype of tumor DNA and expression profiles of microRNA have also been proposed. Such variety of information will help us understand the heterogeneity of neuroblastoma biology and further, the combined use of these signatures will be beneficial in predicting prognosis with high accuracy, as well as choosing a suitable therapy for the individual patient.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication

BACKGROUND

Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma.

METHODS

A single-gene resolution aCGH profiling was integrated with microarray-based gene expression profiling data to distinguish genetic copy number alterations that were strongly associated with transcriptional changes in two neuroblastoma cell lines. FISH analysis using a hotspot tumor tissue microarray of 37 paraffin-embedded neuroblastoma samples and in silico data mining for gene expression information obtained from previously published studies including up to 445 healthy nervous system samples and 123 neuroblastoma samples were used to evaluate the clinical significance and transcriptional consequences of the detected alterations and to identify subsequently activated gene(s).

RESULTS

In addition to the anticipated high-level amplification and subsequent overexpression of MYCN, MEIS1, CDK4 and MDM2 oncogenes, the aCGH analysis revealed numerous other genetic alterations, including microamplifications at 2p and 12q24.11. Most interestingly, we identified and investigated the clinical relevance of a previously poorly characterized amplicon at 12q24.31. FISH analysis showed low-level gain of 12q24.31 in 14 of 33 (42%) neuroblastomas. Patients with the low-level gain had an intermediate prognosis in comparison to patients with MYCN amplification (poor prognosis) and to those with no MYCN amplification or 12q24.31 gain (good prognosis) (P = 0.001). Using the in silico data mining approach, we identified elevated expression of five genes located at the 12q24.31 amplicon in neuroblastoma (DIABLO, ZCCHC8, RSRC2, KNTC1 and MPHOSPH9). Among these, DIABLO showed the strongest activation suggesting a putative role in neuroblastoma progression.

CONCLUSIONS

The presented systematic and rapid framework, which integrates aCGH, gene expression and tissue data to obtain novel targets and biomarkers for cancer, identified a low-level gain of the 12q24.31 as a potential new biomarker for neuroblastoma progression. Furthermore, results of in silico data mining suggest a new neuroblastoma target gene, DIABLO, within this region, whose functional and therapeutic role remains to be elucidated in follow-up studies.

Molecular pathogenesis of peripheral neuroblastic tumors

Neuroblastoma (NB) is an embryonal cancer of the sympathetic nervous system observed in early childhood, characterized by a broad spectrum of clinical behaviors, ranging from spontaneous regression to fatal outcome despite aggressive therapies. NB accounts for 8-10% of pediatric cancers and 15% of the deaths attributable to malignant conditions in children. Interestingly, NB may occur in various contexts, being mostly sporadic but also familial or syndromic. This review focuses on recent advances in the identification of the genes and mechanisms implicated in NB pathogenesis. Although the extensive characterization of the genomic aberrations recurrently observed in sporadic NBs provides important insights into the understanding of the clinical heterogeneity of this neoplasm, analysis of familial and syndromic cases also unravels essential clues on the genetic bases of NB. Recently, the ALK gene emerged as an important NB gene, being implicated both in sporadic and familial cases. The identification of gene expression signatures associated with patient's outcome points out the potential of using gene expression profiling to improve clinical management of patients suffering from NB. Finally, based on recent observations integrating genomic analyses, biological data and clinical information, we discuss possible evolution/progression schemes in NB.

Genomic alterations in primary breast cancers compared with their sentinel and more distal lymph node metastases: an aCGH study

Metastatic potential of breast cancer may be associated with specific genomic alterations and the earliest metastases are likely to be found in the sentinel lymph nodes (SLN). Using array comparative genomic hybridization (aCGH), we compared the genomes of primary breast invasive duct carcinomas (IDCs), their sentinel and more distal lymph node metastases, and IDCs without nodal metastasis. Thirty-three samples from 22 patients with IDC were subjected to aCGH: 8 IDC samples from patients without lymph node metastasis, 11 IDCs associated with SLN metastases out of which 7 had paired samples of metastases, and 14 samples of lymph node metastases out of which 8 were sentinel-distal pairs from 4 patients. aCGH data were analyzed by correlation of genomic profiles, cluster analysis, segmentation, and peak identification. Quantitative real-time PCR was used for data validation. We observed high genomic similarity between primary tumors and their nodal metastases as well as between metastases to the sentinel and distal lymph nodes. Several recurrent alterations were detected preferentially in IDC associated with SLN metastases compared to IDCs without metastasis. Amplification within the 17q24.1-24.2(59.96-62.76 Mb) region was associated with presence of sentinel or distal lymph node metastases; larger tumor size and higher histological grade. In our samples, there were genomic events associated with metastatic progression, which could be detected in both primary tumors and LN metastases. Gain on 17q24.1-24.2 is a candidate region for further testing as a predictor of nodal metastasis.

Mcl1 downregulation sensitizes neuroblastoma to cytotoxic chemotherapy and small molecule Bcl2-family antagonists

Neuroblastoma (NB) is a common, highly lethal pediatric cancer, with treatment failures largely attributable to the emergence of chemoresistance. The pro-survival Bcl2 homology (BH) proteins critically regulate apoptosis, and may represent important therapeutic targets for restoring drug sensitivity in NB. We used a human NB tumor tissue microarray to survey the expression of pro-survival BH proteins Mcl1 and Bcl2, and correlated expression to clinical prognostic factors and survival. Primary NB tumors heterogeneously expressed Mcl1 or Bcl2, with high expression correlating to high-risk phenotype. Co-expression is infrequent (11%), but correlates to reduced survival. Using RNA interference, we investigated the functional relevance of Mcl1 and Bcl2 in high-risk NB cell lines (SK-N-AS, IMR-5, NLF). Mcl1 knockdown induced apoptosis in all NB cell lines, while Bcl2 knockdown inhibited only NLF, suggesting functional heterogeneity. Finally, we determined the relevance of Mcl1 in resistance to conventional chemotherapy (etoposide, doxorubicin) and small molecule Bcl2-family antagonists (ABT-737 and AT-101). Mcl1 silencing augmented sensitivity to chemotherapeutics 2- to 300-fold, while Bcl2 silencing did not, even in Bcl2-sensitive NLF cells. Resistance to ABT-737, which targets Bcl2/-w/-x, was overcome by Mcl1 knockdown. AT-101, which also neutralizes Mcl1, had single-agent cytotoxicity, further augmented by Mcl1 knockdown. In conclusion, Mcl1 appears a predominant pro-survival protein contributing to chemoresistance in NB, and Mcl1 inactivation may represent a novel therapeutic strategy. Optimization of compounds with higher Mcl1 affinity, or combination with additional Mcl1 antagonists, may enhance the clinical utility of this approach.

Allelotyping, microsatellite instability, and BRAF mutation analyses in common and atypical melanocytic nevi and primary cutaneous melanomas

Loss of heterozygosity (LOH) in several chromosomal regions is found in melanoma, and it has been partially studied in nevi. BRAF mutations are found in melanoma and nevi and in colorectal cancer are linked to mismatch repair deficiency. We studied early genetic events involved in melanomagenesis through analysis of allelic loss, microsatellite instability (MSI), and BRAF mutations. DNA extracted from microdissected cells of 22 common nevi, 23 atypical nevi, and 25 primary cutaneous melanomas were examined for LOH and MSI by polymerase chain reaction-based analysis of 24 microsatellite markers and BRAF mutation. Allelic loss index was higher in atypical nevi (0.20) and melanomas (0.27) than common nevi (0.07). LOH was frequent at 9p21, 17q21, 6q23, and 5q35 in melanoma. LOH at any of this loci occurred in 27% of common nevi, 57% of atypical nevi, and 68% of melanomas. BRAF mutations were not related to MSI presence and MSI index was not related with BRAF mutational status. Similar genetic alterations in atypical nevi and melanomas support the concept of atypical nevus as melanoma precursor. Novel deletion loci at 5q35 and 17q21 (BRCA1) in atypical nevi and melanomas were identified. Mismatch repair deficiency is not a crucial event for BRAF mutation in melanocytic tumors.

Comparison of primary neuroblastoma tumors and derivative early-passage cell lines using genome-wide single nucleotide polymorphism array analysis

Stromal contamination is one of the major confounding factors in the analysis of solid tumor samples by single nucleotide polymorphism (SNP) arrays. As we propose to use genome-wide SNP microarray analysis as a diagnostic platform for neuroblastoma, the sensitivity, specificity, and accuracy of these studies must be optimized. To investigate the effects of stromal contamination, we derived early-passage cell lines from nine primary tumors and compared their genomic signature with that of the primary tumors using 100K SNP arrays. The average concordance between tumor and cell line for raw loss of heterozygosity (LOH) calls was 96% (range, 91-99%) and for raw copy number alterations, 71% (range, 43-87%). In general, there were a larger number of LOH events identified in the cell lines compared with the matched tumor samples (mean increase, 3.2% +/- 1.9%). We have developed an algorithm that shows that the presence of stroma contributes to under-reporting of LOH and copy number loss. Notable findings in this sample set were uniparental disomy of chromosome arms 11p, 1q, 14q, and 15q and a novel area of amplification on chromosome band 11p15. Our analysis shows that LOH was identified significantly more often in derived cell lines compared with the original tumor samples. Although these may in part be due to clonal selection during adaptation to tissue culture, our study indicates that stromal contamination may be a major contributing factor in underestimation of LOH and copy number loss events.

ODC1 is a critical determinant of MYCN oncogenesis and a therapeutic target in neuroblastoma

Neuroblastoma is a frequently lethal childhood tumor in which MYC gene deregulation, commonly as MYCN amplification, portends poor outcome. Identifying the requisite biopathways downstream of MYC may provide therapeutic opportunities. We used transcriptome analyses to show that MYCN-amplified neuroblastomas have coordinately deregulated myriad polyamine enzymes (including ODC1, SRM, SMS, AMD1, OAZ2, and SMOX) to enhance polyamine biosynthesis. High-risk tumors without MYCN amplification also overexpress ODC1, the rate-limiting enzyme in polyamine biosynthesis, when compared with lower-risk tumors, suggesting that this pathway may be pivotal. Indeed, elevated ODC1 (independent of MYCN amplification) was associated with reduced survival in a large independent neuroblastoma cohort. As polyamines are essential for cell survival and linked to cancer progression, we studied polyamine antagonism to test for metabolic dependence on this pathway in neuroblastoma. The Odc inhibitor alpha-difluoromethylornithine (DFMO) inhibited neuroblast proliferation in vitro and suppressed oncogenesis in vivo. DFMO treatment of neuroblastoma-prone genetically engineered mice (TH-MYCN) extended tumor latency and survival in homozygous mice and prevented oncogenesis in hemizygous mice. In the latter, transient Odc ablation permanently prevented tumor onset consistent with a time-limited window for embryonal tumor initiation. Importantly, we show that DFMO augments antitumor efficacy of conventional cytotoxics in vivo. This work implicates polyamine biosynthesis as an arbiter of MYCN oncogenesis and shows initial efficacy for polyamine depletion strategies in neuroblastoma, a strategy that may have utility for this and other MYC-driven embryonal tumors.

Characterization of amplicons in neuroblastoma: high-resolution mapping using DNA microarrays, relationship with outcome, and identification of overexpressed genes

Somatically acquired chromosomal imbalances are a key feature of neuroblastoma, a heterogeneous pediatric solid tumor. Among these alterations, genomic amplification targeting the MYCN oncogene and observed in about 25-30% of the cases, strongly correlates with advanced stage and poor outcome. In this work, we have used BAC and SNP arrays as well as gene expression arrays to characterize amplifications in neuroblastoma. Eighty-eight distinct BACs defining high-level amplification events were identified in 65 samples, including 43 tumors and 22 cell lines. Although the highest recurrence was observed on chromosome 2, clones on chromosomes 8, 12, 16, and 17 also revealed genomic amplification in several samples. A detailed analysis of the 2p22-2p25 MYCN containing region indicated highly complex patterns in a number of cases. Coamplifications involving MYCN and other regions were explored by FISH in three cell lines. High-resolution arrays then allowed us to further refine the mapping of 25 amplicons in 19 samples, either reducing the size of a single continuous amplicon or increasing the complexity by highlighting multiple noncontiguous regions of amplification. Combined analysis of gene expression profiling and array-CGH data indicated that 12 to 25% of the genes that are targeted by genomic amplification are actually over-expressed in tumor cells, several of them having already been implicated in cancer. Finally, our results suggest that the presence of amplicons localized outside of chromosome 2, in addition to MYCN amplification, may be linked to a particularly severe outcome in neuroblastoma patients.

Genomic differences between pure ductal carcinoma in situ of the breast and that associated with invasive disease: a calibrated aCGH study

PURPOSE

In the quest for new targets, genomes of ductal carcinoma in situ (DCIS) and infiltrating duct carcinoma (IDC) have been compared previously; however, genomic alterations associated with cancer progression were difficult to identify. We hypothesized that significant events can be detected by comparing lesions with a broader range of behavior: from pure DCIS to IDC associated with lymph node metastasis.

EXPERIMENTAL DESIGN

Array comparative genomic hybridization, calibrated by self-self hybridization tests, was used to study 6 cases of pure DCIS and 17 cases of DCIS paired with IDC where 8 tumors had spread to the local lymph nodes.

RESULTS

Pure DCIS exhibited a marginally higher degree of genomic complexity than DCIS and IDC components of invasive tumors. The latter two showed similarity between tumors and between components of the same tumor with several regions detected preferentially compared with pure DCIS. IDC associated with lymph node metastases showed similarity of genomic profiles as a group. Gain on 17q22-24.2 was associated with higher histologic grade, large IDC size, lymphatic/vascular invasion, and lymph node metastasis (P < 0.05).

CONCLUSIONS

Our findings suggest that DCIS and IDC are associated with specific genomic events. DCIS associated with IDC is genomically similar to the invasive component and therefore may represent either a clone with high invasive potential or invasive cancer spreading through the ducts. Specifically, gain on 17q22-24.2 is a candidate region for further testing as a predictor of invasion when detected in DCIS and predictor of nodal metastasis when detected in DCIS or IDC.

A functional screen identifies miR-34a as a candidate neuroblastoma tumor suppressor gene

MicroRNAs are small noncoding RNAs that have critical roles in regulating a number of cellular functions through transcriptional silencing. They have been implicated as oncogenes and tumor suppressor genes (oncomirs) in several human neoplasms. We used an integrated genomics and functional screening strategy to identify potential oncomirs in the pediatric neoplasm neuroblastoma. We first identified microRNAs that map within chromosomal regions that we and others have defined as frequently deleted (1p36, 3p22, and 11q23-24) or gained (17q23) in high-risk neuroblastoma. We then transiently transfected microRNA precursor mimics or inhibitors into a panel of six neuroblastoma cell lines that we characterized for these genomic aberrations. The majority of transfections showed no phenotypic effect, but the miR-34a (1p36) and miR-34c (11q23) mimics showed dramatic growth inhibition in cell lines with 1p36 hemizygous deletion. In contrast, there was no growth inhibition by these mimics in cell lines without 1p36 deletions. Quantitative reverse transcription-PCR showed a perfect correlation of absent miR-34a expression in cell lines with a 1p36 aberration and phenotypic effect after mimetic add-back. Expression of miR-34a was also decreased in primary tumors (n = 54) with 1p36 deletion (P = 0.009), but no mutations were discovered in resequencing of the miR-34a locus in 30 neuroblastoma cell lines. Flow cytometric time series analyses showed that the likely mechanism of miR-34a growth inhibition is through cell cycle arrest followed by apoptosis. BCL2 and MYCN were identified as miR-34a targets and likely mediators of the tumor suppressor phenotypic effect. These data support miR-34a as a tumor suppressor gene in human neuroblastoma.

The MYCN oncogene is a direct target of miR-34a

Loss of 1p36 heterozygosity commonly occurs with MYCN amplification in neuroblastoma tumors, and both are associated with an aggressive phenotype. Database searches identified five microRNAs that map to the commonly deleted region of 1p36 and we hypothesized that the loss of one or more of these microRNAs contributes to the malignant phenotype of MYCN-amplified tumors. By bioinformatic analysis, we identified that three out of the five microRNAs target MYCN and of these miR-34a caused the most significant suppression of cell growth through increased apoptosis and decreased DNA synthesis in neuroblastoma cell lines with MYCN amplification. Quantitative RT-PCR showed that neuroblastoma tumors with 1p36 loss expressed lower level of miR-34a than those with normal copies of 1p36. Furthermore, we demonstrated that MYCN is a direct target of miR-34a. Finally, using a series of mRNA expression profiling experiments, we identified other potential direct targets of miR-34a, and pathway analysis demonstrated that miR-34a suppresses cell-cycle genes and induces several neural-related genes. This study demonstrates one important regulatory role of miR-34a in cell growth and MYCN suppression in neuroblastoma.

KIF1Bbeta functions as a haploinsufficient tumor suppressor gene mapped to chromosome 1p36.2 by inducing apoptotic cell death

Deletion of the distal region of chromosome 1 frequently occurs in a variety of human cancers, including aggressive neuroblastoma. Previously, we have identified a 500-kb homozygously deleted region at chromosome 1p36.2 harboring at least six genes in a neuroblastoma-derived cell line NB1/C201. Among them, only KIF1Bbeta, a member of the kinesin superfamily proteins, induced apoptotic cell death. These results prompted us to address whether KIF1Bbeta could be a tumor suppressor gene mapped to chromosome 1p36 in neuroblastoma. Hemizygous deletion of KIF1Bbeta in primary neuroblastomas was significantly correlated with advanced stages (p = 0.0013) and MYCN amplification (p < 0.001), whereas the mutation rate of the KIF1Bbeta gene was infrequent. Although KIF1Bbeta allelic loss was significantly associated with a decrease in KIF1Bbeta mRNA levels, its promoter region was not hypermethylated. Additionally, expression of KIF1Bbeta was markedly down-regulated in advanced stages of tumors (p < 0.001). Enforced expression of KIF1Bbeta resulted in an induction of apoptotic cell death in association with an increase in the number of cells entered into the G2/M phase of the cell cycle, whereas its knockdown by either short interfering RNA or by a genetic suppressor element led to an accelerated cell proliferation or enhanced tumor formation in nude mice, respectively. Furthermore, we demonstrated that the rod region unique to KIF1Bbeta is critical for the induction of apoptotic cell death in a p53-independent manner. Thus, KIF1Bbeta may act as a haploinsufficient tumor suppressor, and its allelic loss may be involved in the pathogenesis of neuroblastoma and other cancers.

Increased copy number of the TERT and TERC telomerase subunit genes in cancer cells

Telomerase is a ribonucleoprotein enzyme complex that adds telomeric repeats to the ends of chromosomes. The core telomerase components are the telomerase reverse transcriptase (TERT) catalytic subunit, and the telomerase RNA (TR) template subunit. In most cancers, telomerase is expressed at levels that are substantially higher than in normal cells. A known consequence of telomerase up-regulation which is considered to play a critical role in oncogenesis is maintenance of telomere length, and thus evasion by cancer cells of the normal limits on proliferation that are associated with the steady decrease in telomere length that accompanies proliferation of normal cells. It has also been suggested that telomerase up-regulation confers other advantages on cancer cells independent of its enzymatic activity. The mechanisms responsible for up-regulation of telomerase in cancer are incompletely understood. Here we review evidence suggesting that this frequently results from increased copy number of the genes encoding telomerase components. The TERT gene is located at human chromosome band 5p15.33, and the telomerase RNA component (TERC) gene that encodes TR is at 3q26.3. Chromosomal gains and gene amplifications involving chromosome arms 5p and 3q are among the most frequent in human tumors. Increased TERT and TERC gene dosage has been detected frequently in a variety of human cancers, and clonal evolution of cells with increased TERT or TERC copy number has been observed, suggesting a growth advantage in cells with increased TERT or TERC gene dosage.

Neuroblastoma cells isolated from bone marrow metastases contain a naturally enriched tumor-initiating cell

Neuroblastoma is a heterogeneous pediatric tumor thought to arise from the embryonic neural crest. Identification of the cell responsible for propagating neuroblastomas is essential to understanding this often recurrent, rapidly progressing disease. We have isolated and characterized putative tumor-initiating cells from 16 tumors and bone marrow metastases from patients in all neuroblastoma risk groups. Dissociated cells from tumors or bone marrow grew as spheres in conditions used to culture neural crest stem cells, were capable of self-renewal, and exhibited chromosomal aberrations typical of neuroblastoma. Primary spheres from all tumor risk groups differentiated under neurogenic conditions to form neurons. Tumor spheres from low-risk tumors frequently formed large neuronal networks, whereas those from high-risk tumors rarely did. As few as 10 passaged tumor sphere cells from aggressive neuroblastoma injected orthotopically into severe combined immunodeficient/Beige mice formed large neuroblastoma tumors that metastasized to liver, spleen, contralateral adrenal and kidney, and lung. Furthermore, highly tumorigenic tumor spheres were isolated from the bone marrow of patients in clinical remission, suggesting that this population of cells may predict clinical behavior and serve as a biomarker for minimal residual disease in high-risk patients. Our data indicate that high-risk neuroblastoma contains a cell with cancer stem cell properties that is enriched in tumor-initiating capacity. These cells may serve as a model system to identify the molecular determinants of neuroblastoma and to develop new therapeutic strategies for this tumor.

Molecular biology of neuroblastoma

Neuroblastoma is one of the most frequently occurring solid tumours in children, especially in the first year of life, when it accounts for 50% of all tumours. It is the second most common cause of death in children, only preceded by accidents. The most peculiar characteristic of neuroblastoma is its clinical heterogeneity. Approximately half of the cases are classified as high risk, with overall survival rates around 40% despite intensive multimodal therapy. Nevertheless, other subsets of neuroblastomas will undergo spontaneous regression and others will show very slow progression. Despite many advances in the past three decades, neuroblastoma has remained an enigmatic challenge to clinical and basic scientists. Elucidation of the exact molecular pathways of neuroblastoma will enable researchers and clinicians to stratify the disease and adapt therapy to the risk of relapse or progression. This review focuses on recent advances in our understanding of the biology of this complex paediatric tumour. Neuroblastoma is already one of the first examples for the use of tumoral genetic markers as a tool for defining tumour behaviour and to aid clinical staging.

Genes proximal and distal to MYCN are highly expressed in human neuroblastoma as visualized by comparative expressed sequence hybridization

MYCN amplification is associated with poor prognosis in neuroblastoma disease. To improve our understanding of the influence of the MYCN amplicon and its corresponding expression, we investigated the 2p expression pattern of MYCN amplified (n = 13) and nonamplified (n = 4) cell lines and corresponding primary tumors (n = 3) using the comparative expressed sequence hybridization technique. All but one MYCN amplified cell line displayed overexpression at 2p. Expression peaks were observed frequently at 2pter and less frequently at 2p24 (MYCN locus), 2p23.3-23.2, and/or 2p23.1. Importantly, cell lines and two corresponding primary tumors displayed expression peaks at similar loci. No significant 2p24 expression level was observed for those cell lines displaying a low amplification rate (n = 3) by comparative genomic hybridization. Only the cell lines with an enhanced peak at 2p23.2-23.3 displayed coamplification of the ALK gene (2p23.2), reported to be associated with unfavorable prognosis. Finally, two of four cell lines without MYCN amplification, both derived from patients with poor outcome, also showed an expression peak at 2p23.2. These data indicate that, besides MYCN, other genes proximal and distal to MYCN are highly expressed in neuroblastoma. The prognostic significance of expression peaks at 2p23.2-23.3, independent of MYCN and ALK status, remains to be investigated.

ArrayCGH-based classification of neuroblastoma into genomic subgroups

High-resolution array comparative genomic hybridization (arrayCGH) profiling was performed on 75 primary tumors and 29 cell lines to gain further insight into the genetic heterogeneity of neuroblastoma and to refine genomic subclassification. Using a novel data-mining strategy, three major and two minor genomic subclasses were delineated. Eighty-three percent of tumors could be assigned to the three major genomic subclasses, corresponding to the three known clinically and biologically relevant subsets in neuroblastoma. The remaining subclasses represented (1) tumors with no/few copy number alterations or an atypical pattern of aberrations and (2) tumors with 11q13 amplification. Inspection of individual arrayCGH profiles showed that recurrent genomic imbalances were not exclusively associated with a specific subclass. Of particular notice were tumors with numerical imbalances typically observed in subtype 1 neuroblastoma, in association with genomic features of subtype 2A or 2B. A search for prognostically relevant genomic alterations disclosed 1q gain as a predictive marker for therapy failure within the group of subtype 2A and 2B tumors. In cell lines, a high incidence of 6q loss was observed, with a 3.87-5.32 Mb region of common loss within 6q25.1-6q25.2. Our study clearly illustrates the importance of genomic profiling in relation to tumor behavior in neuroblastoma. We propose that genome-wide assessment of copy number alterations should ideally be included in the genetic workup of neuroblastoma. Further multicentric studies on large tumor series are warranted in order to improve therapeutic stratification in conjunction with other features such as age at diagnosis, tumor stage, and gene expression signatures.

Neuroblastomas have distinct genomic DNA profiles that predict clinical phenotype and regional gene expression

Neuroblastoma is a heterogeneous neoplasm that has served as a paradigm for the clinical utility of somatically acquired genomic aberrations. DNA copy number alterations (CNA) are currently used to predict prognosis, including MYCN amplification and deletions at chromosome bands 1p36 and 11q23. We predicted that genome-wide assessment of DNA aberrations in neuroblastoma tumors would provide a more precise estimation of clinical phenotype, and could be used to predict outcome. We measured CNAs in a representative set of 82 diagnostic tumors on a customized high-resolution BAC array-based CGH platform supplemented with additional clones across 1p36, 2p24, 3p21-22, 11q14-24, and 16p12-13, and integrated these data with RNA expression data. We used an unbiased statistical method to define a set of minimal common regions (MCRs) of aberration. Unsupervised hierarchical clustering identified four distinct genomic subclasses. First, a subset of tumors with a clinically benign phenotype showed predominantly whole chromosome gains and losses. Second, tumors with MYCN amplification had a unique genomic signature of 1p deletion and 17q gain, but few other rearrangements. Third, tumors with an aggressive clinical phenotype without MYCN amplification, showed multiple structural rearrangements. Most notable were deletions of 3p, 4p, and 11q and gain of 1q, 2p, 12q, and 17q. Lastly, there was a subset of tumors with an aggressive clinical phenotype and no detectable DNA CNAs. The genomic subsets were highly correlated with patient outcome, and individual MCRs remained prognostic in a multivariable model. DNA signature patterns embed important prognostic information in diagnostic neuroblastoma samples, and can identify candidate cancer-related genes.

Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays

Genomic aberrations recurrent in a particular cancer type can be important prognostic markers for tumor progression. Typically in early tumorigenesis, cells incur a breakdown of the DNA replication machinery that results in an accumulation of genomic aberrations in the form of duplications, deletions, translocations, and other genomic alterations. Microarray methods allow for finer mapping of these aberrations than has previously been possible; however, data processing and analysis methods have not taken full advantage of this higher resolution. Attention has primarily been given to analysis on the single sample level, where multiple adjacent probes are necessarily used as replicates for the local region containing their target sequences. However, regions of concordant aberration can be short enough to be detected by only one, or very few, array elements. We describe a method called Multiple Sample Analysis for assessing the significance of concordant genomic aberrations across multiple experiments that does not require a-priori definition of aberration calls for each sample. If there are multiple samples, representing a class, then by exploiting the replication across samples our method can detect concordant aberrations at much higher resolution than can be derived from current single sample approaches. Additionally, this method provides a meaningful approach to addressing population-based questions such as determining important regions for a cancer subtype of interest or determining regions of copy number variation in a population. Multiple Sample Analysis also provides single sample aberration calls in the locations of significant concordance, producing high resolution calls per sample, in concordant regions. The approach is demonstrated on a dataset representing a challenging but important resource: breast tumors that have been formalin-fixed, paraffin-embedded, archived, and subsequently UV-laser capture microdissected and hybridized to two-channel BAC arrays using an amplification protocol. We demonstrate the accurate detection on simulated data, and on real datasets involving known regions of aberration within subtypes of breast cancer at a resolution consistent with that of the array. Similarly, we apply our method to previously published datasets, including a 250K SNP array, and verify known results as well as detect novel regions of concordant aberration. The algorithm has been fully implemented and tested and is freely available as a Java application at http://www.cbil.upenn.edu/MSA.

Cancer is a genetic disease caused by genomic mutations that confer an increased ability to proliferate and survive in a specific environment. It is now known that many regions of genomic DNA are deleted or amplified in specific cancer types. These aberrations are believed to occur randomly in the genome. If these aberrations overlap more than would be expected by chance across individual occurrences of the cancer this suggests a selective pressure on this aberration. These conserved aberrations likely represent regions that are important for the development, progression, and survival of a specific cancer type in its environment. We present a method for identifying these conserved aberrations within a class of samples. The applications for this method include accurate high resolution mapping of aberrations characteristic of cancer subtypes as well as other genetic diseases and determination of conserved copy number variations in the population. With the use of high resolution microarray methods we have profiled different tumor types. We have been able to create high resolution profiles of conserved aberrations in specific cancer types. These conserved aberrations are prime targets for cancer therapies and many of these regions have already been used to develop effective cancer therapeutics.

Telomerase reverse transcriptase haploinsufficiency and telomere length in individuals with 5p- syndrome

Telomerase, which maintains the ends of chromosomes, consists of two core components, the telomerase reverse transcriptase (TERT) and the telomerase RNA (TERC). Haploinsufficiency for TERC or TERT leads to progressive telomere shortening and autosomal dominant dyskeratosis congenita (DC). The clinical manifestations of autosomal dominant DC are thought to occur when telomeres become critically short, but the rate of telomere shortening in this condition is unknown. Here, we investigated the consequences of de novo TERT gene deletions in a large cohort of individuals with 5p– syndrome. The study group included 41 individuals in which the chromosome deletion resulted in loss of one copy of the TERT gene at 5p15.33. Telomere length in peripheral blood cells from these individuals, although within the normal range, was on average shorter than in normal controls. The shortening was more significant in older individuals suggesting an accelerated age-dependent shortening. In contrast, individuals with autosomal dominant DC due to an inherited TERC gene deletion had very short telomeres, and the telomeres were equally short regardless of the age. Although some individuals with 5p– syndrome showed clinical features that were reminiscent of autosomal dominant DC, these features did not correlate with telomere length, suggesting that these were not caused by critically short telomeres. We conclude that a TERT gene deletion leads to slightly shorter telomeres within one generation. However, our results suggest that several generations of TERT haploinsufficiency are needed to produce the very short telomeres seen in patients with DC.

Identification and characterization of DNA imbalances in neuroblastoma by high-resolution oligonucleotide array comparative genomic hybridization

Neuroblastoma (NB) is a pediatric tumor characterized by high genetic heterogeneity. Although the prognostic significance of some genomic abnormalities (i.e., MYCN amplification, 1p loss, and 17q gain) is recognized, genes that are involved in chromosome rearrangements remain largely unknown. Considerable progress has been made over the last years in characterizing DNA abnormalities by metaphase comparative genomic hybridization (mCGH) and array CGH (aCGH). Here we report a pilot study of 5 localized and 12 disseminated NB by 44,000 and of 4 out of 17 cases by 244,000 oligonucleotide aCGH. Localized tumors were predominantly characterized by losses of whole chromosomes 3, 4, 10, and 16, and gains of 6, 7, 8, 13, 17, 18, and 20, whereas disseminated tumors showed several structural aberrations including 17q gain and 3p and 11q losses. Characterization of chromosome 2p in MYCN-amplified NB revealed several structural rearrangements with regions of gain interspersed among sites of amplification, indicating that the MYCN amplicon may encompass several genes. The high-resolution zooming in chromosomal aberrant regions detected several micro-deletions and micro-amplifications in the NB genome. Our results indicate that the increased sensitivity of the aCGH also allows the identification of DNA aberrations in challenging samples (i.e., NB showing tissue and genetic heterogeneity).

Telomere instability and cancer

Telomeres are required to preserve genome integrity, chromosome stability, nuclear architecture and chromosome pairing during meiosis. Given that telomerase activity is limiting or absent in most somatic tissues, shortening of telomeres during development and aging is the rule. In vitro, telomere length operates as a mechanism to prevent uncontrolled cell growth and therefore defines the proliferation potential of a cell. In vitro, in somatic cells that have lost proliferation control, shortening of telomeres becomes the main source of genome instability leading to genetic or epigenetic changes that may allow cells to become immortal and to acquire tumor phenotypes. In vivo, mice models have indisputably shown both the protective and the promoting role of very short telomeres in cancer development. In humans, although telomere shortening and other types of telomere dysfunction probably contribute to the genome instability often detected in tumors, the specific contributions of such instability to the development of cancer remain largely undetermined.

A novel technique for measuring variations in DNA copy-number: competitive genomic polymerase chain reaction

Background

Changes in genomic copy number occur in many human diseases including cancer. Characterization of these changes is important for both basic understanding and diagnosis of these diseases. Microarrays have recently become the standard technique and are commercially available. However, it is useful to have an affordable technique to complement them.

Results

We describe a novel polymerase chain reaction (PCR)-based technique, termed competitive genomic PCR (CGP). The main characteristic of CGP is that different adaptors are added to the sample and control genomic DNAs after appropriate restriction enzyme digestion. These adaptor-supplemented DNAs are subjected to competitive PCR using an adaptor-primer and a locus-specific primer. The amplified products are then separated according to size differences between the adaptors. CGP eliminates the tedious steps inherent in quantitative PCR and achieves moderate throughput. Assays with different X chromosome numbers showed that it can provide accurate quantification. High-resolution analysis of neuroblastoma cell lines around the MYCN locus revealed novel junctions for amplification, which were not detected by a commercial array.

Conclusion

CGP is a moderate throughput technique for analyzing changes in genomic copy numbers. Because CGP can measure any genomic locus using PCR primers, it is especially useful for detailed analysis of a genomic region of interest.