Sharon J. Diskin

The genetic landscape of high-risk neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 affected individuals (cases) using a combination of whole-exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per Mb (0.48 nonsilent) and notably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, and an additional 7.1% had focal deletions), MYCN (1.7%, causing a recurrent p.Pro44Leu alteration) and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1 and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies that rely on frequently altered oncogenic drivers.

Copy number variation at 1q21.1 associated with neuroblastoma

Common copy number variations (CNVs) represent a significant source of genetic diversity, yet their influence on phenotypic variability, including disease susceptibility, remains poorly understood. To address this problem in human cancer, we performed a genome-wide association study of CNVs in the childhood cancer neuroblastoma, a disease in which single nucleotide polymorphism variations are known to influence susceptibility. We first genotyped 846 Caucasian neuroblastoma patients and 803 healthy Caucasian controls at ∼550,000 single nucleotide polymorphisms, and performed a CNV-based test for association. We then replicated significant observations in two independent sample sets comprised of a total of 595 cases and 3,357 controls. Here we describe the identification of a common CNV at chromosome 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets. This CNV was validated by quantitative polymerase chain reaction, fluorescent in situ hybridization and analysis of matched tumour specimens, and was shown to be heritable in an independent set of 713 cancer-free parent–offspring trios. We identified a previously unknown transcript within the CNV that showed high sequence similarity to several neuroblastoma breakpoint family (NBPF) genes and represents a new member of this gene family (NBPF23). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and the expression level was strictly correlated with CNV state in neuroblastoma cells. These data demonstrate that inherited copy number variation at 1q21.1 is associated with neuroblastoma and implicate a previously unknown neuroblastoma breakpoint family gene in early tumorigenesis of this childhood cancer.

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. (Mol Cancer Res 2008;6(5):735–42)

Adjustment of genomic waves in signal intensities from whole-genome SNP genotyping platforms

Whole-genome microarrays with large-insert clones designed to determine DNA copy number often show variation in hybridization intensity that is related to the genomic position of the clones. We found these ‘genomic waves’ to be present in Illumina and Affymetrix SNP genotyping arrays, confirming that they are not platform-specific. The causes of genomic waves are not well-understood, and they may prevent accurate inference of copy number variations (CNVs). By measuring DNA concentration for 1444 samples and by genotyping the same sample multiple times with varying DNA quantity, we demonstrated that DNA quantity correlates with the magnitude of waves. We further showed that wavy signal patterns correlate best with GC content, among multiple genomic features considered. To measure the magnitude of waves, we proposed a GC-wave factor (GCWF) measure, which is a reliable predictor of DNA quantity (correlation coefficient = 0.994 based on samples with serial dilution). Finally, we developed a computational approach by fitting regression models with GC content included as a predictor variable, and we show that this approach improves the accuracy of CNV detection. With the wide application of whole-genome SNP genotyping techniques, our wave adjustment method will be important for taking full advantage of genotyped samples for CNV analysis.

Chromosome 6p22 Locus Associated with Clinically Aggressive Neuroblastoma

BACKGROUND Neuroblastoma is a malignant condition of the developing sympathetic nervous system that most commonly affects young children and is often lethal. Its cause is not known. METHODS We performed a genomewide association study by first genotyping blood DNA samples from 1032 patients with neuroblastoma and 2043 control subjects of European descent using the Illumina HumanHap550 BeadChip. Samples from three independent groups of patients with neuroblastoma (a total of 720 patients) and 2128 control subjects were then genotyped to replicate significant associations. RESULTS We observed a significant association between neuroblastoma and the common minor alleles of three consecutive single-nucleotide polymorphisms (SNPs) at chromosome band 6p22 and containing the predicted genes FLJ22536 and FLJ44180 (P=1.71x10(-9) to 7.01x10(-10); allelic odds ratio, 1.39 to 1.40). Homozygosity for the at-risk G allele of the most significantly associated SNP, rs6939340, resulted in an increased likelihood of the development of neuroblastoma (odds ratio, 1.97; 95% confidence interval, 1.58 to 2.45). Subsequent genotyping of the three 6p22 SNPs in three independent case series confirmed our observation of an association (P=9.33x10(-15) at rs6939340 for joint analysis). Patients with neuroblastoma who were homozygous for the risk alleles at 6p22 were more likely to have metastatic (stage 4) disease (P=0.02), amplification of the MYCN oncogene in the tumor cells (P=0.006), and disease relapse (P=0.01). CONCLUSIONS A common genetic variation at chromosome band 6p22 is associated with susceptibility to neuroblastoma.

Common variations in BARD1 influence susceptibility to high-risk neuroblastoma

We conducted a SNP-based genome-wide association study (GWAS) focused on the high-risk subset of neuroblastoma. As our previous unbiased GWAS showed strong association of common 6p22 SNP alleles with aggressive neuroblastoma, we restricted our analysis here to 397 high-risk cases compared to 2,043 controls. We detected new significant association of six SNPs at 2q35 within the BARD1 locus (Pallelic = 2.35 × 10−9–2.25 × 10−8). We confirmed each SNP association in a second series of 189 high-risk cases and 1,178 controls (Pallelic = 7.90 × 10−7–2.77 × 10−4). We also tested the two most significant SNPs (rs6435862, rs3768716) in two additional independent high-risk neuroblastoma case series, yielding combined allelic odds ratios of 1.68 each (P = 8.65 × 10−18 and 2.74 × 10−16, respectively). We also found significant association with known BARD1 nonsynonymous SNPs. These data show that common variation in BARD1 contributes to the etiology of the aggressive and most clinically relevant subset of human neuroblastoma.

Integrative genomics identifies LMO1 as a neuroblastoma oncogene

Neuroblastoma is a childhood cancer of the sympathetic nervous system that accounts for approximately 10% of all paediatric oncology deaths. To identify genetic risk factors for neuroblastoma, we performed a genome-wide association study (GWAS) on 2,251 patients and 6,097 control subjects of European ancestry from four case series. Here we report a significant association within LIM domain only 1 (LMO1) at 11p15.4 (rs110419, combined P = 5.2 × 10−16, odds ratio of risk allele = 1.34 (95% confidence interval 1.25–1.44)). The signal was enriched in the subset of patients with the most aggressive form of the disease. LMO1 encodes a cysteine-rich transcriptional regulator, and its paralogues (LMO2, LMO3 and LMO4) have each been previously implicated in cancer. In parallel, we analysed genome-wide DNA copy number alterations in 701 primary tumours. We found that the LMO1 locus was aberrant in 12.4% through a duplication event, and that this event was associated with more advanced disease (P < 0.0001) and survival (P = 0.041). The germline single nucleotide polymorphism (SNP) risk alleles and somatic copy number gains were associated with increased LMO1 expression in neuroblastoma cell lines and primary tumours, consistent with a gain-of-function role in tumorigenesis. Short hairpin RNA (shRNA)-mediated depletion of LMO1 inhibited growth of neuroblastoma cells with high LMO1 expression, whereas forced expression of LMO1 in neuroblastoma cells with low LMO1 expression enhanced proliferation. These data show that common polymorphisms at the LMO1 locus are strongly associated with susceptibility to developing neuroblastoma, but also may influence the likelihood of further somatic alterations at this locus, leading to malignant progression.

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.

RNAi screen of the protein kinome identifies checkpoint kinase 1 (CHK1) as a therapeutic target in neuroblastoma

Neuroblastoma is a childhood cancer that is often fatal despite intense multimodality therapy. In an effort to identify therapeutic targets for this disease, we performed a comprehensive loss-of-function screen of the protein kinome. Thirty kinases showed significant cellular cytotoxicity when depleted, with loss of the cell cycle checkpoint kinase 1 (CHK1/CHEK1) being the most potent. CHK1 mRNA expression was higher in MYC–Neuroblastoma-related (MYCN)–amplified (P < 0.0001) and high-risk (P = 0.03) tumors. Western blotting revealed that CHK1 was constitutively phosphorylated at the ataxia telangiectasia response kinase target site Ser345 and the autophosphorylation site Ser296 in neuroblastoma cell lines. This pattern was also seen in six of eight high-risk primary tumors but not in control nonneuroblastoma cell lines or in seven of eight low-risk primary tumors. Neuroblastoma cells were sensitive to the two CHK1 inhibitors SB21807 and TCS2312, with median IC50 values of 564 nM and 548 nM, respectively. In contrast, the control lines had high micromolar IC50 values, indicating a strong correlation between CHK1 phosphorylation and CHK1 inhibitor sensitivity (P = 0.0004). Furthermore, cell cycle analysis revealed that CHK1 inhibition in neuroblastoma cells caused apoptosis during S-phase, consistent with its role in replication fork progression. CHK1 inhibitor sensitivity correlated with total MYC(N) protein levels, and inducing MYCN in retinal pigmented epithelial cells resulted in CHK1 phosphorylation, which caused growth inhibition when inhibited. These data show the power of a functional RNAi screen to identify tractable therapeutical targets in neuroblastoma and support CHK1 inhibition strategies in this disease.

The Plasmodium genome database

Designing and mining a eukaryotic genomics resource.