Meningiomas and schwannomas: Molecular subgroup classification found by expression arrays

Microarray gene expression profiling is a high-throughput system used to identify differentially expressed genes and regulation patterns, and to discover new tumor markers. As the molecular pathogenesis of meningiomas and schwannomas, characterized by NF2 gene alterations, remains unclear and suitable molecular targets need to be identified, we used low density cDNA microarrays to establish expression patterns of 96 cancer-related genes on 23 schwannomas, 42 meningiomas and 3 normal cerebral meninges. We also performed a mutational analysis of the NF2 gene (PCR, dHPLC, Sequencing and MLPA), a search for 22q LOH and an analysis of gene silencing by promoter hypermethylation (MS-MLPA). Results showed a high frequency of NF2 gene mutations (40%), increased 22q LOH as aggressiveness increased, frequent losses and gains by MLPA in benign meningiomas, and gene expression silencing by hypermethylation. Array analysis showed decreased expression of 7 genes in meningiomas. Unsupervised analyses identified 2 molecular subgroups for both meningiomas and schwannomas showing 38 and 20 differentially expressed genes, respectively, and 19 genes differentially expressed between the two tumor types. These findings provide a molecular subgroup classification for meningiomas and schwannomas with possible implications for clinical practice.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

MDM2 and CDK4 immunohistochemistry is a valuable tool in the differential diagnosis of low-grade osteosarcomas and other primary fibro-osseous lesions of the bone.

Low-grade osteosarcoma is a rare malignancy that may be subdivided into two main subgroups on the basis of location in relation to the bone cortex, that is, parosteal osteosarcoma and low-grade central osteosarcoma. Their histological appearance is quite similar and characterized by spindle cell stroma with low-to-moderate cellularity and well-differentiated anastomosing bone trabeculae. Low-grade osteosarcomas have a simple genetic profile with supernumerary ring chromosomes comprising amplification of chromosome 12q13-15, including the cyclin-dependent kinase 4 (CDK4) and murine double-minute type 2 (MDM2) gene region. Low-grade osteosarcoma can be confused with fibrous and fibro-osseous lesions such as fibromatosis and fibrous dysplasia on radiological and histological findings. We investigated MDM2-CDK4 immunohistochemical expression in a series of 72 low-grade osteosarcomas and 107 fibrous or fibro-osseous lesions of the bone or paraosseous soft tissue. The MDM2-CDK4 amplification status of low-grade osteosarcoma was also evaluated by comparative genomic hybridization array in 18 cases, and the MDM2 amplification status was evaluated by fluorescence in situ hybridization or quantitative real-time polymerase chain reaction in 31 cases of benign fibrous and fibro-osseous lesions. MDM2-CDK4 immunostaining and MDM2 amplification by fluorescence in situ hybridization or quantitative real-time polymerase chain reaction were investigated in a control group of 23 cases of primary high-grade bone sarcoma, including 20 conventional high-grade osteosarcomas, two pleomorphic spindle cell sarcomas/malignant fibrous histiocytomas and one leiomyosarcoma. The results showed that MDM2 and/or CDK4 immunoreactivity was present in 89% of low-grade osteosarcoma specimens. All benign fibrous and fibro-osseous lesions and the tumors of the control group were negative for MDM2 and CDK4. These results were consistent with the MDM2 and CDK4 amplification results. In conclusion, immunohistochemical expression of MDM2 and CDK4 is specific and provides sensitive markers for the diagnosis of low-grade osteosarcomas, helping to differentiate them from benign fibrous and fibro-osseous lesions, particularly in cases with atypical radio-clinical presentation and/or limited biopsy samples.

Genetic testing in patients with obesity.

The obesity epidemic is associated with the recent availability of highly palatable and inexpensive caloric food as well as important changes in lifestyle. Genetic factors, however, play a key role in regulating energy balance and numerous twin studies have estimated the BMI heritability between 40 and 70%. While common variants, identified through genome-wide association studies (GWAS) point toward new pathways, their effect size are too low to be of any use in the clinic. This review therefore concentrates on genes and genomic regions associated with very high risks of human obesity. Although there are no consensus guidelines, we review how the knowledge on these "causal factors" can be translated into the clinic for diagnostic purposes. We propose genetic workups guided by clinical manifestations in patients with severe early-onset obesity. While etiological diagnoses are unequivocal in a minority of patients, new genomic tools such as Comparative Genomic Hybridization (CGH) array, have allowed the identification of novel "causal" loci and next-generation sequencing brings the promise of accelerated pace for discoveries relevant to clinical practice.

Copy number variation characteristics in subpopulations of patients with autism spectrum disorders.

Autism spectrum disorders (ASDs) are a heterogeneous group of disorders with a complex genetic etiology. We used high-resolution whole genome array-based comparative genomic hybridization (array-CGH) to screen 223 ASD patients for gene dose alterations associated with susceptibility for autism. Clinically significant copy number variations (CNVs) were identified in 18 individuals (8%), of which 9 cases (4%) had de novo aberrations. In addition, 20 individuals (9%) were shown to have CNVs of unclear clinical relevance. Among these, 13 cases carried rare but inherited CNVs that may increase the risk for developing ASDs, while parental samples were unavailable in the remaining seven cases. Classification of all patients into different phenotypic and inheritance pattern groups indicated the presence of different CNV patterns in different patient groups. Clinically relevant CNVs were more common in syndromic cases compared to non-syndromic cases. Rare inherited CNVs were present in a higher proportion of ASD cases having first- or second-degree relatives with an ASD-related neuropsychiatric phenotype in comparison with cases without reported heredity (P = 0.0096). We conclude that rare CNVs, encompassing potential candidate regions for ASDs, increase the susceptibility for the development of ASDs and related neuropsychiatric disorders giving us further insight into the complex genetics underlying ASDs.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Genome-wide arrays in routine diagnostics of hematological malignancies.

Over the last three decades, cytogenetic analysis of malignancies has become an integral part of disease evaluation and prediction of prognosis or responsiveness to therapy. In most diagnostic laboratories, conventional karyotyping, in conjunction with targeted fluorescence in situ hybridization analysis, is routinely performed to detect recurrent aberrations with prognostic implications. However, the genetic complexity of cancer cells requires a sensitive genome-wide analysis, enabling the detection of small genomic changes in a mixed cell population, as well as of regions of homozygosity. The advent of comprehensive high-resolution genomic tools, such as molecular karyotyping using comparative genomic hybridization or single-nucleotide polymorphism microarrays, has overcome many of the limitations of traditional cytogenetic techniques and has been used to study complex genomic lesions in, for example, leukemia. The clinical impact of the genomic copy-number and copy-neutral alterations identified by microarray technologies is growing rapidly and genome-wide array analysis is evolving into a diagnostic tool, to better identify high-risk patients and predict patients' outcomes from their genomic profiles. Here, we review the added clinical value of an array-based genome-wide screen in leukemia, and discuss the technical challenges and an interpretation workflow in applying arrays in the acquired cytogenetic diagnostic setting.

Strong smooth muscle differentiation is dependent on myocardin gene amplification in most human retroperitoneal leiomyosarcomas.

Myocardin (MYOCD), a serum response factor (SRF) transcriptional cofactor, is essential for cardiac and smooth muscle development and differentiation. We show here by array-based comparative genomic hybridization, fluorescence in situ hybridization, and expression analysis approaches that MYOCD gene is highly amplified and overexpressed in human retroperitoneal leiomyosarcomas (LMS), a very aggressive well-differentiated tumor. MYOCD inactivation by shRNA in a human LMS cell line with MYOCD locus amplification leads to a dramatic decrease of smooth muscle differentiation and strongly reduces cell migration. Moreover, forced MYOCD expression in three undifferentiated sarcoma cell lines and in one liposarcoma cell line confers a strong smooth muscle differentiation phenotype and increased migration abilities. Collectively, these results show that human retroperitoneal LMS differentiation is dependent on MYOCD amplification/overexpression, suggesting that in these well-differentiated LMS, differentiation could be a consequence of an acquired genomic alteration. In this hypothesis, these tumors would not necessarily derive from cells initially committed to smooth muscle differentiation. These data also provide new insights on the cellular origin of these sarcomas and on the complex connections between oncogenesis and differentiation in mesenchymal tumors.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Comparative genomic and phylogeographic analysis of Mycobacterium leprae.

Reductive evolution and massive pseudogene formation have shaped the 3.31-Mb genome of Mycobacterium leprae, an unculturable obligate pathogen that causes leprosy in humans. The complete genome sequence of M. leprae strain Br4923 from Brazil was obtained by conventional methods (6x coverage), and Illumina resequencing technology was used to obtain the sequences of strains Thai53 (38x coverage) and NHDP63 (46x coverage) from Thailand and the United States, respectively. Whole-genome comparisons with the previously sequenced TN strain from India revealed that the four strains share 99.995% sequence identity and differ only in 215 polymorphic sites, mainly SNPs, and by 5 pseudogenes. Sixteen interrelated SNP subtypes were defined by genotyping both extant and extinct strains of M. leprae from around the world. The 16 SNP subtypes showed a strong geographical association that reflects the migration patterns of early humans and trade routes, with the Silk Road linking Europe to China having contributed to the spread of leprosy.

Molecular features of hepatosplenic T-cell lymphoma unravels potential novel therapeutic targets.

The pathogenesis of hepatosplenic T-cell lymphoma (HSTL), a rare entity mostly derived from γδ T cells and usually with a fatal outcome, remains largely unknown. In this study, HSTL samples (7γδ and 2αβ) and the DERL2 HSTL cell line were subjected to combined gene-expression profiling and array-based comparative genomic hybridization. Compared with other T-cell lymphomas, HSTL had a distinct molecular signature irrespective of TCR cell lineage. Compared with peripheral T-cell lymphoma, not otherwise specified and normal γδ T cells, HSTL overexpressed genes encoding NK-cell-associated molecules, oncogenes (FOS and VAV3), the sphingosine-1-phosphatase receptor 5 involved in cell trafficking, and the tyrosine kinase SYK, whereas the tumor-suppressor gene AIM1 (absent in melanoma 1) was among the most down-expressed. We found highly methylated CpG islands of AIM1 in DERL2 cells, and decitabine treatment induced a significant increase in AIM1 transcripts. Syk was present in HSTL cells and DERL2 cells contained phosphorylated Syk and were sensitive to a Syk inhibitor in vitro. Genomic profiles confirmed recurrent isochromosome 7q (n = 6/9) without alterations at the SYK and AIM1 loci. Our results identify a distinct molecular signature for HSTL and highlight oncogenic pathways that offer rationale for exploring new therapeutic options such as Syk inhibitors and demethylating agents.

Beta-catenin Status in P?aediatric Medulloblastomas: correlation of immunohistochemical expression with mutational status, genetic profiles, and clinical characteristics.

Medulloblastoma is the most frequent malignant paediatric brain tumour. The activation of the Wnt/beta-catenin pathway occurs in 10-15% of medulloblastomas and has been recently described as a marker for favourable patient outcome. We report a series of 72 paediatric medulloblastomas evaluated for beta-catenin protein expression, CTNNB1 mutations, and comparative genomic hybridization. Gene expression profiles were also available in a subset of 40 cases. Immunostaining of beta-catenin showed extensive nuclear staining (>50% of the tumour cells) in six cases and focal nuclear staining (<10% of cells) in three cases. The other cases either exhibited a signal strictly limited to the cytoplasm (58 cases) or were negative (five cases). CTNNB1 mutations were detected in all beta-catenin extensively nucleopositive cases. The expression profiles of these cases documented strong activation of the Wnt/beta-catenin pathway. Remarkably, five out of these six tumours showed a complete loss of chromosome 6. In contrast, cases with focal nuclear beta-catenin staining, as well as tumours with negative or cytoplasmic staining, never demonstrated CTNNB1 mutation, Wnt/beta-catenin pathway activation or chromosome 6 loss. Patients with extensive nuclear staining were significantly older at diagnosis and were in continuous complete remission after a mean follow-up of 75.7 months (range 27.5-121.2 months) from diagnosis. All three patients with focal nuclear staining of beta-catenin died within 36 months from diagnosis. Altogether, these data confirm and extend previous observations that CTNNB1-mutated tumours represent a distinct molecular subgroup of medulloblastomas with favourable outcome, indicating that therapy de-escalation should be considered. International consensus on the definition criteria of this distinct medulloblastoma subgroup should be achieved.

A common cognitive, psychiatric, and dysmorphic phenotype in carriers of NRXN1 deletion

Deletions in the 2p16.3 region that includes the neurexin (NRXN1) gene are associated with intellectual disability and various psychiatric disorders, in particular, autism and schizophrenia. We present three unrelated patients, two adults and one child, in whom we identified an intragenic 2p16.3 deletion within the NRXN1 gene using an oligonucleotide comparative genomic hybridization array. The three patients presented dual diagnosis that consisted of mild intellectual disability and autism and bipolar disorder. Also, they all shared a dysmorphic phenotype characterized by a long face, deep set eyes, and prominent premaxilla. Genetic analysis of family members showed two inherited deletions. A comprehensive neuropsychological examination of the 2p16.3 deletion carriers revealed the same phenotype, characterized by anxiety disorder, borderline intelligence, and dysexecutive syndrome. The cognitive pattern of dysexecutive syndrome with poor working memory and reduced attention switching, mental flexibility, and verbal fluency was the same than those of the adult probands. We suggest that in addition to intellectual disability and psychiatric disease, NRXN1 deletion is a risk factor for a characteristic cognitive and dysmorphic profile. The new cognitive phenotype found in the 2p16.3 deletion carriers suggests that 2p16.3 deletions might have a wide variable expressivity instead of incomplete penetrance

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Molecular Characteristics of Neuroblastoma with Special Reference to Novel Prognostic Factors and Diagnostic Applications

Molecular Characteristics of Neuroblastoma with Special Reference to Novel Prognostic Factors and Diagnostic Applications Department of Medical Biochemistry and Genetics Annales Universitatis Turkuensis, Medica-Odontologica, 2009, Turku, Finland Painosalama Oy, Turku, Finland 2009 Background: Neuroblastoma, which is the most common and extensively studied childhood solid cancer, shows a great clinical and biological heterogeneity. Most of the neuroblastoma patients older than one year have poor prognosis despite intensive therapies. The hallmark of neuroblastoma, biological heterogeneity, has hindered the discovery of prognostic tumour markers. At present, few molecular markers, such as MYCN oncogene status, have been adopted into clinical practice. Aims: The aim of the study was to improve the current prognostic methodology of neuroblastoma, especially by taking cognizance of the biological heterogeneity of neuroblastoma. Furthermore, unravelling novel molecular characteristics which associate with neuroblastoma tumour progression and cell differentiation was an additional objective. Results: A new strictly defined selection of neuroblastoma tumour spots of highest proliferation activity, hotspots, appeared to be representative and reliable in an analysis of MYCN amplification status using a chromogenic in situ hybridization technique (CISH). Based on the hotspot tumour tissue microarray immunohistochemistry and high-resolution oligo-array-based comparative genomic hybridization, which was integrated with gene expression and in silico analysis of existing transcriptomics, a polysialylated neural cell adhesion molecule (NCAM) and poorly characterized amplicon at 12q24.31 were discovered to associate with outcome. In addition, we found that a previously considered new neuroblastoma treatment target, the mutated c-kit receptor, was not mutated in neuroblastoma samples. Conclusions: Our studies indicate polysialylated NCAM and 12q24.31 amplicon to be new molecular markers with important value in prognostic evaluation of neuroblastoma. Moreover, the presented hotspot tumour tissue microarray method together with the CISH technique of the MYCN oncogene copy number is directly applicable to clinical use. Key words: neuroblastoma, polysialic acid, neural cell adhesion molecule, MYCN, c-kit, chromogenic in situ hybridization, hotspot