The mechanisms underlying genetic susceptibility at loci discovered by genome-wide association study (GWAS) approaches in human cancer remain largely undefined. In this study we characterized the high-risk neuroblastoma association at the BRCA1-related locus, BARD1, showing that disease-associated variations correlate with increased expression of the oncogenically activated isoform, BARD1β. In neuroblastoma cells, silencing of BARD1β showed genotype-specific cytotoxic effects, including decreased substrate-adherent, anchorage-independent, and foci growth. In established murine fibroblasts, overexpression of BARD1β was sufficient for neoplastic transformation. BARD1β stabilized the Aurora family of kinases in neuroblastoma cells, suggesting both a mechanism for the observed effect and a potential therapeutic strategy. Together, our findings identify BARD1β as an oncogenic driver of high-risk neuroblastoma tumorigenesis, and more generally, they illustrate how robust GWAS signals offer genomic landmarks to identify molecular mechanisms involved in both tumor initiation and malignant progression. The interaction of BARD1β with the Aurora family of kinases lends strong support to the ongoing work to develop Aurora kinase inhibitors for clinically aggressive neuroblastoma.

Sporadic clear cell Renal Cell Carcinoma (ccRCC), the most common type of adult kidney cancer, is often associated with genomic copy number aberrations on chromosomes 3p and 5q. Aberrations on chromosome 3p are associated with inactivation of the tumor suppressor gene von-Hippel Lindau (VHL), which activates the hypoxia inducible factors HIF1α and HIF2α. In contrast, ccRCC genes on chromosome 5q remain to be defined. In this study, we performed an integrated analysis of high-density copy number and gene expression data for 54 sporadic ccRCC tumors that identified the secreted glycoprotein STC2 (stanniocalcin 2) and the proteoglycan VCAN (versican) as potential 5q oncogenes in ccRCC. In functional assays, STC2 and VCAN each promoted tumorigenesis by inhibiting cell death. Using the same approach, we also investigated the two VHL-deficient subtypes of ccRCC, which express both HIF1α and HIF2α (H1H2) or only HIF2α (H2). This analysis revealed a distinct pattern of genomic aberrations in each group, with the H1H2 group displaying, on average, a more aberrant genome than the H2 group. Together our findings provide a significant advance in understanding ccRCC by offering a molecular definition of two subtypes with distinct characteristics as well as two potential chromosome 5q oncogenes, the overexpression of which is sufficient to promote tumorigenesis by limiting cell death.

Neuroblastoma is a childhood cancer of the sympathetic nervous system that accounts for approximately 10% of all paediatric oncology deaths1,2. 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.

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.

We conducted a SNP-based genome-wide association study (GWAS) focused on the high-risk subset of neuroblastoma1. As our previous unbiased GWAS showed strong association of common 6p22 SNP alleles with aggressive neuroblastoma2, we now restricted our analysis to 397 high-risk cases compared to 2,043 controls. We detected new significant association of six SNPs at 2q35 within the BARD1 gene locus (Pallelic = 2.35×10−9 − 2.25×10−8). Each SNP association was confirmed in a second series of 189 high-risk cases and 1,178 controls (Pallelic = 7.90×10−7 − 2.77×10−4). The two most significant SNPs (rs6435862, rs3768716) were also tested 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). Significant association was also found with known BARD1 nsSNPs. These data show that common variation in BARD1 contributes to the etiology of the aggressive and most clinically relevant subset of human 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 cancer, we performed a genome-wide association study (GWAS) of CNVs in the childhood cancer neuroblastoma, a disease where SNP variations are known to influence susceptibility1,2. 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. We identified a common CNV at 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets (Pcombined = 2.97 × 10−17; OR = 2.49, 95% CI: 2.02 to 3.05). This CNV was validated by quantitative PCR, fluorescent in situ hybridization, and analysis of matched tumor specimens, and was shown to be heritable in an independent set of 713 cancer-free trios. We identified a novel transcript within the CNV which showed high sequence similarity to several “Neuroblastoma breakpoint family” (NBPF) genes3,4 and represents a new member of this gene family (NBPFX). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and 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 novel NBPF gene in early tumorigenesis of this childhood cancer.

Background

Neuroblastoma is a malignancy of the developing sympathetic nervous system that most commonly affects young children and is often lethal. The etiology of this embryonal cancer is not known.

Methods

We performed a genome-wide association study by first genotyping 1,032 neuroblastoma patients and 2,043 controls of European descent using the Illumina HumanHap550 BeadChip. Three independent groups of neuroblastoma cases (N=720) and controls (N=2128) were then genotyped to replicate significant associations.

Results

We observed highly significant association between neuroblastoma and the common minor alleles of three single nucleotide polymorphisms (SNPs) within a 94.2 kilobase (Kb) linkage disequilibrium block at chromosome band 6p22 containing the predicted genes FLJ22536 and FLJ44180 (P-value range = 1.71×10-9-7.01×10-10; allelic odds ratio range 1.39-1.40). Homozygosity for the at-risk G allele of the most significantly associated SNP, rs6939340, resulted in an increased likelihood of developing neuroblastoma of 1.97 (95% CI 1.58-2.44). Subsequent genotyping of these 6p22 SNPs in the three independent case series confirmed our observation of association (P=9.33×10-15 at rs6939340 for joint analysis). Furthermore, neuroblastoma patients homozygous for the risk alleles at 6p22 were more likely to develop metastatic (Stage 4) disease (P=0.02), show amplification of the MYCN oncogene in the tumor cells (P=0.006), and to have disease relapse (P=0.01).

Conclusion

Common genetic variation at chromosome band 6p22 is associated with susceptibility to neuroblastoma.

Background

The components of therapy required for patients with INSS Stage 3 neuroblastoma and high risk features remain controversial.

Procedure

A retrospective cohort design was used to determine if intensive chemoradiotherapy with purged autologous bone marrow rescue (ABMT) and/or 13-cis-retinoic acid (13-cis-RA) improved outcome for patients with high-risk neuroblastoma that was not metastatic to distant sites. We identified 72 patients with INSS Stage 3 neuroblastoma enrolled between 1991 and 1996 on the Phase III CCG 3891 randomized trial. Patients were analyzed on an intent-to-treat basis using a log-rank test.

Results

The 5-year event-free survival (EFS) and overall survival (OS) rates for patients with Stage 3 neuroblastoma were 55 +/- 6% and 59% +/- 6%, respectively (n=72). Patients randomized to ABMT (n=20) had 5-year EFS of 65% +/- 11% and OS of 65% +/- 11% compared to 41% +/- 11 (p=0.21) and 46% +/- 11% (p=0.23) for patients randomized to CC (n=23), respectively. Patients randomized to 13-cis-RA (n=23) had 5-year EFS of 70% +/- 10% and OS of 78% +/- 9% compared to 63% +/- 12% (p=0.67) and 67% +/- 12% (p=0.55) for those receiving no further therapy (n=16), respectively. Patients randomized to both ABMT and 13-cis-RA (n=6) had a 5-year EFS of 80% ± 11% and OS of 100%.

Conclusion

Patients with high-risk Stage 3 neuroblastoma have an overall poor prognosis despite aggressive chemoradiotherapy. Further studies are warranted to determine if myeloablative consolidation followed by 13-cis-RA maintenance therapy statistically significantly improves outcome.

SUMMARY

Survival rates for the childhood cancer neuroblastoma have not substantively improved despite dramatic escalation in chemotherapy intensity. Like most human cancers, this embryonal malignancy can be inherited, but the genetic etiology of familial and sporadically occurring neuroblastoma was largely unknown. Here we show that germline mutations in the anaplastic lymphoma kinase gene (ALK) explain the majority of hereditary neuroblastomas, and that activating mutations can also be somatically acquired. We first identified a significant linkage signal at the short arm of chromosome 2 (maximum nonparametric LOD=4.23 at rs1344063) using a whole-genome scan in neuroblastoma pedigrees. Resequencing of regional candidate genes identified three separate missense mutations in the tyrosine kinase domain of ALK (G1128A, R1192P and R1275Q) that segregated with the disease in eight separate families. Examination of 491 sporadically occurring human neuroblastoma samples showed that the ALK locus was gained in 22.8%, and highly amplified in an additional 3.3%, and that these aberrations were highly associated with death from disease (P=0.0003). Resequencing of 194 high-risk neuroblastoma samples showed somatically acquired mutations within the tyrosine kinase domain in 12.4%. Nine of the ten mutations map to critical regions of the kinase domain and were predicted to be oncogenic drivers with high probability. Mutations resulted in constitutive phosphorylation consistent with activation, and targeted knockdown of ALK mRNA resulted in profound growth inhibition of 4 of 4 cell lines harboring mutant or amplified ALK, as well as 2 of 6 wild type for ALK. Our results demonstrate that heritable mutations of ALK are the major cause of familial neuroblastoma, and that germline or acquired activation of this cell surface kinase is a tractable therapeutic target for this lethal pediatric malignancy.

Background

Neuroblastomas are characterized by chromosomal alterations with biological and clinical significance. We analyzed paired blood and primary tumor samples from 22 children with high-risk neuroblastoma for loss of heterozygosity (LOH) and DNA copy number change using the Affymetrix 10K single nucleotide polymorphism (SNP) array.

Findings

Multiple areas of LOH and copy number gain were seen. The most commonly observed area of LOH was on chromosome arm 11q (15/22 samples; 68%). Chromosome 11q LOH was highly associated with occurrence of chromosome 3p LOH: 9 of the 15 samples with 11q LOH had concomitant 3p LOH (P = 0.016). Chromosome 1p LOH was seen in one-third of cases. LOH events on chromosomes 11q and 1p were generally accompanied by copy number loss, indicating hemizygous deletion within these regions. The one exception was on chromosome 11p, where LOH in all four cases was accompanied by normal copy number or diploidy, implying uniparental disomy. Gain of copy number was most frequently observed on chromosome arm 17q (21/22 samples; 95%) and was associated with allelic imbalance in six samples. Amplification of MYCN was also noted, and also amplification of a second gene, ALK, in a single case.

Conclusions

This analysis demonstrates the power of SNP arrays for high-resolution determination of LOH and DNA copy number change in neuroblastoma, a tumor in which specific allelic changes drive clinical outcome and selection of therapy.

Neuroblastoma is a malignant neoplasm of the developing sympathetic nervous system that is notable for its phenotypic diversity. High-risk patients typically have widely disseminated disease at diagnosis and a poor survival probability, but low-risk patients frequently have localized tumors that are almost always cured with little or no chemotherapy. Our genome-wide association study (GWAS) has identified common variants within FLJ22536, BARD1, and LMO1 as significantly associated with neuroblastoma and more robustly associated with high-risk disease. Here we show that a GWAS focused on low-risk cases identified SNPs within DUSP12 at 1q23.3 (P = 2.07×10−6), DDX4 and IL31RA both at 5q11.2 (P = 2.94×10−6 and 6.54×10−7 respectively), and HSD17B12 at 11p11.2 (P = 4.20×10−7) as being associated with the less aggressive form of the disease. These data demonstrate the importance of robust phenotypic data in GWAS analyses and identify additional susceptibility variants for neuroblastoma.

Author Summary

Neuroblastoma is the most common solid tumor outside the central nervous system and is accountable for 10% of the mortality rate of all children's cancers. It has distinctive clinical behaviors and is categorized into different risk groups: high-risk, intermediate-risk, and low-risk. Genome-wide association studies have reported a number of genetic variations predisposing to high-risk neuroblastoma. This study focuses on the low-risk neuroblastoma group and identifies four novel genes (DUSP12, DDX4, IL31RA, and HSD17B12) at three distinct genomic positions that harbor disease-causing variants. This study also reports several gene sets that are enriched in overall neuroblastoma as well as in both high-risk and low-risk groups. Also of importance is that this study adopts a new computational method that identifies genes, instead of only one single nucleotide polymorphism, as disease-causing variants. Shown to have superior power of detection genome-wide association signals for neuroblastoma, the methodology presented in this study has great potential applications in case-control association studies in other diseases.