A 6-gene signature identifies four molecular subgroups of neuroblastoma

Background  

There are currently three postulated genomic subtypes of the childhood tumour neuroblastoma (NB); Type 1, Type 2A, and Type 2B. The most aggressive forms of NB are characterized by amplification of the oncogene MYCN (MNA) and low expression of the favourable marker NTRK1. Recently, mutations or high expression of the familial predisposition gene Anaplastic Lymphoma Kinase (ALK) was associated to unfavourable biology of sporadic NB. Also, various other genes have been linked to NB pathogenesis.     

Delineation of the frequency and boundary of chromosomal copy number variations in paediatric neuroblastoma

Neuroblastoma, the most common solid tumour in early childhood, is characterized by very frequent chromosomal copy number variations (CNVs). While chromosome 2p amplification, 17q gain, 1p and 11q deletion in human neuroblastoma tissues are well-known, the exact frequencies and boundaries of the chromosomal CNVs have not been delineated. We analysed the publicly available single nucleotide polymorphism (SNP) array data which were originally generated by the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative, defined the frequencies and boundaries of chromosomes 2p11.2 – 2p25.3 amplification, 17q11.1-17q25.3 gain, 1p13.3-1p36.33 deletion and 11q13.3-11q25 deletion in neuroblastoma tissues, and identified chromosome 7q14.1 (Chr7:38254795-38346971) and chromosome 14q11.2 (Chr14:21637401-22024617) deletion in blood and bone marrow samples from neuroblastoma patients, but not in tumour tissues. Kaplan Meier analysis showed that double deletion of Chr7q14.1 and Chr14q11.2 correlated with poor prognosis in MYCN gene amplified neuroblastoma patients. In conclusion, the oncogenes amplified or gained and tumour suppressor genes deleted within the boundaries of chromosomal CNVs in tumour tissues should be studied for their roles in tumourigenesis and as therapeutic targets. Focal deletions of Chr7q14.1 and Chr14q11.2 together in blood and bone marrow samples from neuroblastoma patients can be used as a marker for poorer prognosis and more aggressive therapies.     

Investigation of major genetic alterations in neuroblastoma

Neuroblastoma (NB) is the most common extracranial solid tumor in childhood. This malignancy shows a wide spectrum of clinical outcome and its prognosis is conditioned by manifold biological and genetic factors. We investigated the tumor genetic profile and clinical data of 29 patients with NB by multiplex ligation-dependent probe amplification (MLPA) to assess therapeutic risk. In 18 of these tumors, MYCN status was assessed by fluorescence in situ hybridization (FISH). Copy number variation was also determined for confirming MLPA findings in two 6p loci. We found 2p, 7q and 17q gains, and 1p and 11q losses as the most frequent chromosome alterations in this cohort. FISH confirmed all cases of MYCN amplification detected by MLPA. In view of unexpected 6p imbalance, copy number variation of two 6p loci was assessed for validating MLPA findings. Based on clinical data and genetic profiles, patients were stratified in pretreatment risk groups according to international consensus. MLPA proved to be effective for detecting multiple genetic alterations in all chromosome regions as requested by the International Neuroblastoma Risk Group (INRG) for therapeutic stratification. Moreover, this technique proved to be cost effective, reliable, only requiring standard PCR equipment, and attractive for routine analysis. However, the observed 6p imbalances made PKHD1 and DCDC2 inadequate for control loci. This must be considered when designing commercial MLPA kits for NB. Finally, four patients showed a normal MLPA profile, suggesting that NB might have a more complex genetic pattern than the one assessed by presently available MLPA kits.     

UBE3A-mediated regulation of imprinted genes and epigenome-wide marks in human neurons

The dysregulation of genes in neurodevelopmental disorders that lead to social and cognitive phenotypes is a complex, multilayered process involving both genetics and epigenetics. Parent-of-origin effects of deletion and duplication of the 15q11-q13 locus leading to Angelman, Prader-Willi, and Dup15q syndromes are due to imprinted genes, including UBE3A, which is maternally expressed exclusively in neurons. UBE3A encodes a ubiquitin E3 ligase protein with multiple downstream targets, including RING1B, which in turn monoubiquitinates histone variant H2A.Z. To understand the impact of neuronal UBE3A levels on epigenome-wide marks of DNA methylation, histone variant H2A.Z positioning, active H3K4me3 promoter marks, and gene expression, we took a multi-layered genomics approach. We performed an siRNA knockdown of UBE3A in two human neuroblastoma cell lines, including parental SH-SY5Y and the SH(15M) model of Dup15q. Genes differentially methylated across cells with differing UBE3A levels were enriched for functions in gene regulation, DNA binding, and brain morphology. Importantly, we found that altering UBE3A levels had a profound epigenetic effect on the methylation levels of up to half of known imprinted genes. Genes with differential H2A.Z peaks in SH(15M) compared to SH-SY5Y were enriched for ubiquitin and protease functions and associated with autism, hypoactivity, and energy expenditure. Together, these results support a genome-wide epigenetic consequence of altered UBE3A levels in neurons and suggest that UBE3A regulates an imprinted gene network involving DNA methylation patterning and H2A.Z deposition.     

High-resolution array copy number analyses for detection of deletion,gain, amplification and copy-neutral LOH in primary neuroblastoma tumors: Four cases of homozygous deletions of the CDKN2A gene

Background

Neuroblastoma is a very heterogeneous pediatric tumor of the sympathetic nervous system showing clinically significant patterns of genetic alterations. Favorable tumors usually have near-triploid karyotypes with few structural rearrangements. Aggressive stage 4 tumors often have near-diploid or near-tetraploid karyotypes and structural rearrangements. Whole genome approaches for analysis of genome-wide copy number have been used to analyze chromosomal abnormalities in tumor samples. We have used array-based copy number analysis using oligonucleotide single nucleotide polymorphisms (SNP) arrays to analyze the chromosomal structure of a large number of neuroblastoma tumors of different clinical and biological subsets.

Results

Ninety-two neuroblastoma tumors were analyzed with 50 K and/or 250 K SNP arrays from Affymetrix, using CNAG3.0 software. Thirty percent of the tumors harbored 1p deletion, 22% deletion of 11q, 26% had MYCN amplification and 45% 17q gain. Most of the tumors with 1p deletion were found among those with MYCN amplification. Loss of 11q was most commonly seen in tumors without MYCN amplification. In the case of MYCN amplification, two types were identified. One type displayed simple continuous amplicons; the other type harbored more complex rearrangements. MYCN was the only common gene in all cases with amplification. Complex amplification on chromosome 12 was detected in two tumors and three different overlapping regions of amplification were identified. Two regions with homozygous deletions, four cases with CDKN2A deletions in 9p and one case with deletion on 3p (the gene RBMS3) were also detected in the tumors.

Conclusion

SNP arrays provide useful tools for high-resolution characterization of significant chromosomal rearrangements in neuroblastoma tumors. The mapping arrays from Affymetrix provide both copy number and allele-specific information at a resolution of 10–12 kb. Chromosome 9p, especially the gene CDKN2A, is subject to homozygous (four cases) and heterozygous deletions (five cases) in neuroblastoma tumors.

     

Identification of novel prognostic markers in relapsing localized resectable neuroblastoma

Patients with localized resectable neuroblastoma (NB) generally have an excellent prognosis and can be treated by surgery alone, but approximately 10% of them develop local recurrences or metastatic progression. The known predictive risk factors are important for the identification of localized resectable NB patients at risk of relapse and/or progression, who may benefit from early and aggressive treatment. These factors, however, identify only a subset of patients at risk, and the search for novel prognostic markers is warranted. This review focuses on the recent advances in the identification of new prognostic markers. Recently we addressed the search of novel genetic prognostic markers in a selected cohort of patients with stroma-poor localized resectable NB who underwent disease relapse or progression (group 1) or complete remission (group 2). High-resolution array-comparative genomic hybridization (CGH) DNA copy-number analysis technology was used. Chromosome 1p36.22p36.32 loss and 1q22qter gain, detected almost exclusively in group 1 patients, were significantly associated with poor event-free survival (EFS). Increasing evidence points to anaplastic lymphoma kinase (ALK) as a fundamental oncogene associated with NB. The immunohistochemical analysis of sporadic NB localized resectable primary tumors (stage 1-2) showed a correlation between aberrant ALK level of expression and tumor progression and clinical outcome. Moreover, other factors that might influence the clinical behavior of these tumors will be reviewed.     

DNA Hypomethylation Affects Cancer-Related Biological Functions and Genes Relevant in Neuroblastoma Pathogenesis

Neuroblastoma (NB) pathogenesis has been reported to be closely associated with numerous genetic alterations. However, underlying DNA methylation patterns have not been extensively studied in this developmental malignancy. Here, we generated microarray-based DNA methylation profiles of primary neuroblastic tumors. Stringent supervised differential methylation analyses allowed us to identify epigenetic changes characteristic for NB tumors as well as for clinical and biological subtypes of NB. We observed that gene-specific loss of DNA methylation is more prevalent than promoter hypermethylation. Remarkably, such hypomethylation affected cancer-related biological functions and genes relevant to NB pathogenesis such as CCND1, SPRR3, BTC, EGF and FGF6. In particular, differential methylation in CCND1 affected mostly an evolutionary conserved functionally relevant 3′ untranslated region, suggesting that hypomethylation outside promoter regions may play a role in NB pathogenesis. Hypermethylation targeted genes involved in cell development and proliferation such as RASSF1A, POU2F2 or HOXD3, among others. The results derived from this study provide new candidate epigenetic biomarkers associated with NB as well as insights into the molecular pathogenesis of this tumor, which involves a marked gene-specific hypomethylation.     

Genome‐wide promoter methylation analysis identifies epigenetic silencing of MAPK13 in primary cutaneous melanoma

The involvement of epigenetic alterations in the pathogenesis of melanoma is increasingly recognized. Here, we performed genome‐wide DNA methylation analysis of primary cutaneous melanoma and benign melanocytic nevus interrogating 14 495 genes using BeadChip technology. This genome‐wide view of promoter methylation in primary cutaneous melanoma revealed an array of recurrent DNA methylation alterations with potential diagnostic applications. Among 106 frequently hypermethylated genes, there were many novel methylation targets and tumor suppressor genes. Highly recurrent methylation of the HOXA9, MAPK13, CDH11, PLEKHG6, PPP1R3C, and CLDN11 genes was established. Promoter methylation of MAPK13, encoding p38δ, was present in 67% of primary and 85% of metastatic melanomas. Restoration of MAPK13 expression in melanoma cells exhibiting epigenetic silencing of this gene reduced proliferation, indicative of tumor suppressive functions. This study demonstrates that DNA methylation alterations are widespread in melanoma and suggests that epigenetic silencing of MAPK13 contributes to melanoma progression.     

Losses at chromosome 4q are associated with poor survival in operable ductal pancreatic adenocarcinoma

Here we tested the prognostic impact of genomic alterations in operable localized pancreatic ductal adenocarcinoma (PDAC). Fifty-two formalin-fixed and paraffin-embedded primary PDAC were laser micro-dissected and were investigated by comparative genomic hybridization after whole genome amplification using an adapter-linker PCR. Chromosomal gains and losses were correlated to clinico-pathological parameters and clinical follow-up data. The most frequent aberration was loss on chromosome 17p (65%) while the most frequent gains were detected at 2q (41%) and 8q (41%), respectively. The concomitant occurrence of losses at 9p and 17p was found to be statistically significant. Higher rates of chromosomal losses were associated with a more advanced primary tumor stage and losses at 9p and 18q were significantly associated with presence of lymphatic metastasis (chi-square: p = 0.03, p = 0.05, respectively). Deletions on chromosome 4 were of prognostic significance for overall survival and tumor recurrence (Cox-multivariate analysis: p = 0.026 and p = 0.021, respectively). In conclusion our data suggest the common alterations at chromosome 8q, 9p, 17p and 18q as well as the prognostic relevant deletions on chromosome 4q as relevant for PDAC progression. Our comprehensive data from 52 PDAC should provide a basis for future studies with a higher resolution to discover the relevant genes located within the chromosomal aberrations identified.     

snoRNPs Regulate Telomerase Activity in Neuroblastoma and Are Associated with Poor Prognosis

Amplification of the MYCN oncogene is strongly associated with poor prognosis in neuroblastoma (NB). In addition to MYCN amplification, many studies have focused on identifying patients with a poor prognosis based on gene expression profiling. The majority of prognostic signatures today are comprised of large gene lists limiting their clinical application. In addition, although of prognostic significance,most of these signatures fail to identify cellular processes that can explain their relation to prognosis. Here, we determined prognostically predictive genes in a data set containing 251 NBs. Gene Ontology analysis was performed on significant genes with a positive hazard ratio to search for cellular processes associated with poor prognosis. An enrichment in ribonucleoproteins (RNPs) was found. Genes involved in the stabilization and formation of the central small nucleolar RNP (snoRNP) complex were scrutinized using a backward conditional Cox regression resulting in an snoRNP signature consisting of three genes: DKC1, NHP2, and GAR1. The snoRNP signature significantly and independently predicted prognosis when compared to the established clinical risk factors. Association of snoRNP protein expression and prognosis was confirmed using tissue microarrays. Knockdown of snoRNP expression in NB cell lines resulted in reduced telomerase activity and an increase in anaphase bridge frequency. In addition, in patient material, expression of the snoRNP complex was significantly associated with telomerase activity, occurrence of segmental aberrations, and expression-based measurements of chromosomal instability. Together, these results underscore the prognostic value of snoRNP complex expression in NB and suggest a role for snoRNPs in telomere maintenance and genomic stability.     

Genetic predisposition to familial neuroblastoma: identification of two novel genomic regions at 2p and 12p

OBJECTIVES: The rarity of familial neuroblastoma (NB) has allowed only a few linkage studies, most of which did not show any evidence of linkage to regions involved in somatic alterations or to genes implicated in other neurocristopathies seldom associated with NB. We screened a highly informative family with recurrent NB by genome-wide linkage analysis aimed at identifying chromosomal regions for NB predisposing genes. METHODS: A genome-wide screen was performed using 382 microsatellite markers. Multipoint model-based linkage analysis was carried out under a dominant mode of inheritance for the disease using the 'affected only' approach. RESULTS: Our analysis identified two haplotypes co-segregating with the disease on chromosomes 2p and 12p, and yielded maximum lod-score values of 3.01 (p < 0.0001) for markers on both intervals. CONCLUSIONS: Evidence of linkage was reported at 16p in North American families, whereas our studies excluded this interval and indicated other loci for disease predisposition, thus confirming the remarkable genetic heterogeneity of NB. These results suggest an oligogenic inheritance in NB involving more loci in genetic determination of the disease.     

Molecular characterization of acquired tolerance of tumor cells to picropodophyllin (PPP)

Background

Picropodophyllin (PPP) is a promising novel anti-neoplastic agent that efficiently kills tumor cells in vitro and causes tumor regression and increased survival in vivo. We have previously reported that PPP treatment induced moderate tolerance in two out of 10 cell lines only, and here report the acquired genomic and expression alterations associated with PPP selection over 1.5 years of treatment.

Methodology/Principal Findings

Copy number alterations monitored using metaphase and array-based comparative genomic hybridization analyses revealed largely overlapping alterations in parental and maximally tolerant cells. Gain/ amplification of the MYC and PVT1 loci in 8q24.21 were verified on the chromosome level. Abnormalities observed in connection to PPP treatment included regular gains and losses, as well as homozygous losses in 10q24.1-q24.2 and 12p12.3-p13.2 in one of the lines and amplification at 5q11.2 in the other. Abnormalities observed in both tolerant derivatives include amplification/gain of 5q11.2, gain of 11q12.1-q14.3 and gain of 13q33.3-qter. Using Nexus software analysis we combined the array-CGH data with data from gene expression profilings and identified genes that were altered in both inputs. A subset of genes identified as downregulated (ALDH1A3, ANXA1, TLR4 and RAB5A) or upregulated (COX6A1, NFIX, ME1, MAPK and TAP2) were validated by siRNA in the tolerant or parental cells to alter sensitivity to PPP and confirmed to alter sensitivity to PPP in further cell lines.

Conclusions

Long-term PPP selection lead to altered gene expression in PPP tolerant cells with increase as well as decrease of genes involved in cell death such as PTEN and BCL2. In addition, acquired genomic copy number alterations were observed that were often reflected by altered mRNA expression levels for genes in the same regions.     

Frequent epigenetics inactivation of KEAP1 gene in non-small cell lung cancer

The KEAP1/Nrf2 pathway is a master regulator of several redox-sensitive genes implicated in resistance of tumor cells against chemotherapeutic drugs. Recent data suggest that epigenetic mechanisms may play a pivotal role in the regulation of KEAP1 expression. We performed a comprehensive genetic and epigenetic analysis of the KEAP1 gene in 47 non-small cell lung cancer tissues and normal specimens. Promoter methylation analysis was performed using a quantitative methylation specific PCR assay in real time. Methylation at the KEAP1 promoter region was detected in 22 out of the 47 NSCLCs (47%) and in none of the normal tissues analyzed. Somatic mutations were detected in 7 out of the 47 tumors (15%) and loss of heterozygosity (LOH) in 10 out of the 47 (21%) of the cases. Overall, we found at least one molecular alteration in 57% of the cases. Approximately one third of the tumors had two alterations and this feature was associated with higher risk of disease progression in univariate COX regression analysis (HR = 3.62; 95% CI 1.24–10.65, p = 0.02). This result was confirmed by Kaplan-Meier analysis, which demonstrated an association between worst outcome and KEAP1 double alterations (p = 0.01, Log rank test). Our results further suggest that deregulation of the NRF2/KEAP1 system could play a pivotal role in the cancerogenesis of NSCLC. In addition identifying patients with KEAP1 genetic and epigenetic abnormalities may contribute to disease progression prediction and response to therapy in lung cancer patients.     

Alterations of <Emphasis Type="Italic">ATM</Emphasis> and <Emphasis Type="Italic">CADM1</Emphasis> in chromosomal 11q22.3–23.2 region are associated with the development of invasive cervical carcinoma

To understand the importance of chr11q22.3–23.2 region in the development of cervical cancer, we have studied the genetic and epigenetic alterations of the candidate genes ATM, PPP2R1B, SDHD and CADM1 in cervical intraepithelial neoplasia (CIN) and cervical carcinoma (CACX) samples. Our study revealed low expression and high alterations (methylation/deletion) (55–59%) of ATM and CADM1 genes along with poor patient outcome. The alterations of ATM and CADM1 are associated with the progression of tumor from CIN to Stage I/II, thus implying their role in early invasiveness. The two genes, PPP2R1B and SDHD, lying in between ATM and CADM1, have low frequency of alterations, and majority of the alterations are in CACX samples, indicating that their alterations might be associated with disease progression. Expressions (mRNA/protein) of the genes showed concordance with their molecular alterations. Significant co-alteration of ATM and CADM1 points to their synergic action for the development of CACX. Mutation is, however, a rare phenomenon for inactivation of ATM. Association between the alteration of ATM and CHEK1 and poor survival of the patients having co-alterations of ATM and CHEK1 points to the DNA damage response pathway disruption in development of CACX. Thus, our data suggest that inactivation of ATM–CHEK1-associated DNA damage response pathway and CADM1-associated signaling network might have an important role in the development of CACX.     

Apoptosis in peripheral neuroblastic tumors. Immunohistochemical expression of bcl-2 and p53 is related to DNA fragmentation

We examined 111 cases of neuroblastoma (NB), searching for how NB relates to apoptotic control and other prognostic factors. Immunohistochemistry using avidin-biotin-peroxidase was carried out for bcl-2 and p53 proteins. Apoptosis was analyzed by in situ detection of chromosomal breakdown. DNA ladders were detected by electrophoresis and amplification of MYCN was carried out by PCR and Southern blot. Statistical analyses were performed with Pearson's chi2 and Kruskal-Wallis tests and Cox's regression. We found expression of bcl-2 protein mainly in cases of neuroblastoma without differentiation and in stages 3 and 4. Expression of p53 protein showed a correlation with bcl-2 and the apoptotic phenomenon; apoptosis was found mainly in favorable cases. Multivariate analysis showed bcl-2 protein expression to be the most independent risk factor. The study of apoptosis could be important for the design of therapies to treat neuroblastoma.