Integrative genomics identifies distinct molecular classes of neuroblastoma and shows that multiple genes are targeted by regional alterations in DNA copy number.

Neuroblastoma is remarkable for its clinical heterogeneity and is characterized by genomic alterations that are strongly correlated with tumor behavior. The specific genes that influence neuroblastoma biology and are targeted by genomic alterations remain largely unknown. We quantified mRNA expression in a highly annotated series of 101 prospectively collected diagnostic neuroblastoma primary tumors using an oligonucleotide-based microarray. Genomic copy number status at the prognostically relevant loci 1p36, 2p24 (MYCN), 11q23, and 17q23 was determined by PCR and was aberrant in 26, 20, 40, and 38 cases, respectively. In addition, 72 diagnostic neuroblastoma primary tumors assayed in a different laboratory were used as an independent validation set. Unsupervised hierarchical clustering showed that gene expression was highly correlated with genomic alterations and clinical markers of tumor behavior. The vast majority of samples with MYCN amplification and 1p36 loss of heterozygosity (LOH) clustered together on a terminal node of the sample dendrogram, whereas the majority of samples with 11q deletion clustered separately and both of these were largely distinct from the copy number neutral group of tumors. Genes involved in neurodevelopment were broadly overrepresented in the more benign tumors, whereas genes involved in RNA processing and cellular proliferation were highly represented in the most malignant cases. By combining transcriptomic and genomic data, we showed that LOH at 1p and 11q was associated with significantly decreased expression of 122 (61%) and 88 (27%) of the genes mapping to 1p35-36 and all of 11q, respectively, suggesting that multiple genes may be targeted by LOH events. A total of 71 of the 1p35-36 genes were also differentially expressed in the independent validation data set, providing a prioritized list of candidate neuroblastoma suppressor genes. Taken together, these data are consistent with the hypotheses that the neuroblastoma transcriptome is a sensitive marker of underlying tumor biology and that chromosomal deletion events in this cancer likely target multiple genes through alteration in mRNA dosage. Lead positional candidates for neuroblastoma suppressor genes can be inferred from these data, but the potential multiplicity of transcripts involved has significant implications for ongoing gene discovery strategies.