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.

Summary Neuroblastoma is a pediatric malignancy that typically arises in early childhood and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumors with excellent outcomes to widely metastatic disease where long-term survival is approximately 40% despite intensive therapy1. A previous genome-wide association study (GWAS) identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumor cells2. Here we sought to discover the causal DNA variant at this locus and the mechanism by which it leads to neuroblastoma tumorigenesis. We first imputed all possible genotypes across the LMO1 locus and then mapped highly associated single nucleotide polymorphism (SNPs) to areas of chromatin accessibility, evolutionary conservation, and transcription factor binding sites. SNP rs2168101 G>T was the most highly associated variant (combined P=7.47×10-29, Odds Ratio 0.65, 95% CI: 0.60-0.70) and resided in a super-enhancer defined by extensive acetylation of histone H3 lysine 27 within the first intron of LMO1. The ancestral G-allele that is associated with tumor formation resides in a conserved GATA transcription factor binding motif. We show that the newly evolved protective TATA allele is associated with decreased total LMO1 expression (P=0.028) in neuroblastoma primary tumors and ablates GATA3 binding (P<0.0001). We demonstrate allelic imbalance favoring the G-containing strand in tumors heterozygous for this SNP as demonstrated both by RNA sequencing (P<0.0001) and reporter assays (P=0.002). These findings show that a recently evolved polymorphism within a super-enhancer element in the first intron of LMO1 influences neuroblastoma susceptibility through differential GATA transcription factor binding and direct modulation of LMO1 expression in cis, and this leads to an oncogenic dependency in tumor cells.

Neuroblastoma is a cancer of the sympathetic nervous system that accounts for approximately 10% of all pediatric oncology deaths1. Here we report on a genome-wide association study of 2,817 neuroblastoma cases and 7,473 controls. We identified two new associations at 6q16, the first within HACE1 (rs4336470; combined P = 2.7 × 10−11, odds ratio 1.26, 95% CI: 1.18–1.35) and the second within LIN28B (rs17065417; combined P = 1.2 × 10−8, odds ratio 1.38, 95% CI: 1.23–1.54). Expression of LIN28B and let-7 miRNA correlated with rs17065417 genotype in neuroblastoma cell lines, and we observed significant growth inhibition upon depletion of LIN28B specifically in neuroblastoma cells homozygous for the risk allele. Low HACE1 and high LIN28B expression in diagnostic primary neuroblastomas were associated with worse overall survival (P = 0.008 and 0.014, respectively). Taken together, we show that common variants in HACE1 and LIN28B influence neuroblastoma susceptibility and that both genes likely play a role in disease progression.

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.