Corrigendum: Transcriptomic profiling of 39 commonly-used neuroblastoma cell lines.

This corrects the article DOI: 10.1038/sdata.2017.33.

Evaluation of Genetic Predisposition for MYCN-Amplified Neuroblastoma.

To investigate genetic predispositions for MYCN-amplified neuroblastoma, we performed a meta-analysis of three genome-wide association studies totaling 615 MYCN-amplified high-risk neuroblastoma cases and 1869 MYCN-nonamplified non-high-risk neuroblastoma cases as controls using a fixed-effects model with inverse variance weighting. All statistical tests were two-sided. We identified a novel locus at 3p21.31 indexed by the single nucleotide polymorphism (SNP) rs80059929 (odds ratio [OR] = 2.95, 95% confidence interval [CI] = 2.17 to 4.02, Pmeta = 6.47 × 10-12) associated with MYCN-amplified neuroblastoma, which was replicated in 127 MYCN-amplified cases and 254 non-high-risk controls (OR = 2.30, 95% CI = 1.12 to 4.69, Preplication = .02). To confirm this signal is exclusive to MYCN-amplified tumors, we performed a second meta-analysis comparing 728 MYCN-nonamplified high-risk patients to identical controls. rs80059929 was not statistically significant in MYCN-nonamplified high-risk patients (OR = 1.24, 95% CI = 0.90 to 1.71, Pmeta = .19). SNP rs80059929 is within intron 16 in the KIF15 gene. Additionally, the previously reported LMO1 neuroblastoma risk locus was statistically significant only in patients with MYCN-nonamplified high-risk tumors (OR = 0.63, 95% CI = 0.53 to 0.75, Pmeta = 1.51 × 10-8; Pmeta = .95). Our results indicate that common genetic variation predisposes to different neuroblastoma genotypes, including the likelihood of somatic MYCN-amplification.

Common variants in MMP20 at 11q22.2 predispose to 11q deletion and neuroblastoma risk.

MYCN amplification and 11q deletion are two inversely correlated prognostic factors of poor outcome in neuroblastoma. Here we identify common variants at 11q22.2 within MMP20 that associate with neuroblastoma cases harboring 11q deletion (rs10895322), using GWAS in 113 European-American cases and 5109 ancestry-matched controls. The association is replicated in 44 independent cases and 1902 controls. Our study yields novel insights into the genetic underpinnings of neuroblastoma, demonstrating that the inherited common variants reported contribute to the origin of intra-tumor genetic heterogeneity in neuroblastoma.Chromosomal abnormalities such as 11q deletion are associated with poor prognosis in neuroblastoma. Here, the authors perform a genome-wide association study and identify an association between a variant within a Matrix metalloproteinase (MMP) gene member, MMP20, and 11q-deletion subtype neuroblastoma.

Common variants upstream of MLF1 at 3q25 and within CPZ at 4p16 associated with neuroblastoma.

Neuroblastoma is a cancer of the developing sympathetic nervous system that most commonly presents in young children and accounts for approximately 12% of pediatric oncology deaths. Here, we report on a genome-wide association study (GWAS) in a discovery cohort or 2,101 cases and 4,202 controls of European ancestry. We identify two new association signals at 3q25 and 4p16 that replicated robustly in multiple independent cohorts comprising 1,163 cases and 4,396 controls (3q25: rs6441201 combined P = 1.2x10-11, Odds Ratio 1.23, 95% CI:1.16-1.31; 4p16: rs3796727 combined P = 1.26x10-12, Odds Ratio 1.30, 95% CI: 1.21-1.40). The 4p16 signal maps within the carboxypeptidase Z (CPZ) gene. The 3q25 signal resides within the arginine/serine-rich coiled-coil 1 (RSRC1) gene and upstream of the myeloid leukemia factor 1 (MLF1) gene. Increased expression of MLF1 was observed in neuroblastoma cells homozygous for the rs6441201 risk allele (P = 0.02), and significant growth inhibition was observed upon depletion of MLF1 (P < 0.0001) in neuroblastoma cells. Taken together, we show that common DNA variants within CPZ at 4p16 and upstream of MLF1 at 3q25 influence neuroblastoma susceptibility and MLF1 likely plays an important role in neuroblastoma tumorigenesis.

Transcriptomic profiling of 39 commonly-used neuroblastoma cell lines.

Neuroblastoma cell lines are an important and cost-effective model used to study oncogenic drivers of the disease. While many of these cell lines have been previously characterized with SNP, methylation, and/or mRNA expression microarrays, there has not been an effort to comprehensively sequence these cell lines. Here, we present raw whole transcriptome data generated by RNA sequencing of 39 commonly-used neuroblastoma cell lines. These data can be used to perform differential expression analysis based on a genetic aberration or phenotype in neuroblastoma (e.g., MYCN amplification status, ALK mutation status, chromosome arm 1p, 11q and/or 17q status, sensitivity to pharmacologic perturbation). Additionally, we designed this experiment to enable structural variant and/or long-noncoding RNA analysis across these cell lines. Finally, as more DNase/ATAC and histone/transcription factor ChIP sequencing is performed in these cell lines, our RNA-Seq data will be an important complement to inform transcriptional targets as well as regulatory (enhancer or repressor) elements in neuroblastoma.

Genetic predisposition to neuroblastoma mediated by a LMO1 super-enhancer polymorphism.

Neuroblastoma is a paediatric malignancy that typically arises in early childhood, and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumours with excellent outcomes to widely metastatic disease in which long-term survival is approximately 40% despite intensive therapy. A previous genome-wide association study identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumour cells. Here we investigate 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. We show that SNP rs2168101 G>T is the most highly associated variant (combined P = 7.47 × 10(-29), odds ratio 0.65, 95% confidence interval 0.60-0.70), and resides 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 tumour 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 tumours, and ablates GATA3 binding (P < 0.0001). We demonstrate allelic imbalance favouring the G-containing strand in tumours heterozygous for this SNP, as demonstrated both by RNA sequencing (P < 0.0001) and reporter assays (P = 0.002). These findings indicate 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 tumour cells.

CASC15-S Is a Tumor Suppressor lncRNA at the 6p22 Neuroblastoma Susceptibility Locus.

Chromosome 6p22 was identified recently as a neuroblastoma susceptibility locus, but its mechanistic contributions to tumorigenesis are as yet undefined. Here we report that the most highly significant single-nucleotide polymorphism (SNP) associations reside within CASC15, a long noncoding RNA that we define as a tumor suppressor at 6p22. Low-level expression of a short CASC15 isoform (CASC15-S) associated highly with advanced neuroblastoma and poor patient survival. In human neuroblastoma cells, attenuating CASC15-S increased cellular growth and migratory capacity. Gene expression analysis revealed downregulation of neuroblastoma-specific markers in cells with attenuated CASC15-S, with concomitant increases in cell adhesion and extracellular matrix transcripts. Altogether, our results point to CASC15-S as a mediator of neural growth and differentiation, which impacts neuroblastoma initiation and progression.

Common genetic variants in NEFL influence gene expression and neuroblastoma risk.

The genetic etiology of sporadic neuroblastoma is still largely obscure. In a genome-wide association study, we identified single-nucleotide polymorphisms (SNP) associated with neuroblastoma at the CASC15, BARD1, LMO1, DUSP12, HSD17B12, HACE1, and LIN28B gene loci, but these explain only a small fraction of neuroblastoma heritability. Other neuroblastoma susceptibility genes are likely hidden among signals discarded by the multiple testing corrections. In this study, we evaluated eight additional genes selected as candidates for further study based on proven involvement in neuroblastoma differentiation. SNPs at these candidate genes were tested for association with disease susceptibility in 2,101 cases and 4,202 controls, with the associations found replicated in an independent cohort of 459 cases and 809 controls. Replicated associations were further studied for cis-effect using gene expression, transient overexpression, silencing, and cellular differentiation assays. The neurofilament gene NEFL harbored three SNPs associated with neuroblastoma (rs11994014: Pcombined = 0.0050; OR, 0.88; rs2979704: Pcombined = 0.0072; OR, 0.87; rs1059111: Pcombined = 0.0049; OR, 0.86). The protective allele of rs1059111 correlated with increased NEFL expression. Biologic investigations showed that ectopic overexpression of NEFL inhibited cell growth specifically in neuroblastoma cells carrying the protective allele. NEFL overexpression also enhanced differentiation and impaired the proliferation and anchorage-independent growth of cells with protective allele and basal NEFL expression, while impairing invasiveness and proliferation of cells homozygous for the risk genotype. Clinically, high levels of NEFL expression in primary neuroblastoma specimens were associated with better overall survival (P = 0.03; HR, 0.68). Our results show that common variants of NEFL influence neuroblastoma susceptibility and they establish that NEFL expression influences disease initiation and progression.

Rare variants in TP53 and susceptibility to neuroblastoma.

TP53 is the most frequently mutated gene in human malignancies; however, de novo somatic mutations in childhood embryonal cancers such as neuroblastoma are rare. We report on the analysis of three independent case-control cohorts comprising 10290 individuals and demonstrate that rs78378222 and rs35850753, rare germline variants in linkage disequilibrium that map to the 3' untranslated region (UTR) of TP53 and 5' UTR of the Δ133 isoform of TP53, respectively, are robustly associated with neuroblastoma (rs35850753: odds ratio [OR] = 2.7, 95% confidence interval [CI] = 2.0 to 3.6, P combined = 3.43×10(-12); rs78378222: OR = 2.3, 95% CI = 1.8 to 2.9, P combined = 2.03×10(-11)). All statistical tests were two-sided. These findings add neuroblastoma to the complex repertoire of human cancers influenced by the rs78378222 hypomorphic allele, which impairs proper termination and polyadenylation of TP53 transcripts. Future studies using whole-genome sequencing data are likely to reveal additional rare variants with large effect sizes contributing to neuroblastoma tumorigenesis.

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.

Common variation at 6q16 within HACE1 and LIN28B influences susceptibility to neuroblastoma.

Neuroblastoma is a cancer of the sympathetic nervous system that accounts for approximately 10% of all pediatric oncology deaths. Here, we report 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% confidence interval (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 that were 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, these data show that common variants in HACE1 and LIN28B influence neuroblastoma susceptibility and indicate that both genes likely have a role in disease progression.

Common variation at BARD1 results in the expression of an oncogenic isoform that influences neuroblastoma susceptibility and oncogenicity.

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-adherence, anchorage-independence, 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.

Replication of neuroblastoma SNP association at the BARD1 locus in African-Americans.

BACKGROUND: Neuroblastoma is an often fatal pediatric cancer more frequent in European-American than African-American children. African-American children, however, are at higher risk for the more severe form of neuroblastoma and have worse overall survival than European-American children. Genome-wide association studies (GWAS) have identified several single-nucleotide polymorphisms (SNP) associated to neuroblastoma in children of European descent. Knowledge of their association to neuroblastoma in African-American children is still lacking. METHODS: We genotyped and imputed SNPs located in three gene regions reported to be associated to neuroblastoma in children of European descent, and tested them for association in 390 African-American patients with neuroblastoma compared with 2,500 healthy, ethnically matched controls. RESULTS: SNPs in the BARD1 gene region show a similar pattern of association to neuroblastoma in African-American and European-American children. The more restricted extent of linkage disequilibrium in the African-American population suggests a smaller candidate region for the putative causal variants than previously reported. Limited association was observed at the other two gene regions tested, including LMO1 in 11p15 and FLJ22536 in 6p22. CONCLUSIONS: Common BARD1 SNPs affect risk of neuroblastoma in African-Americans. The role of other SNPs associated to neuroblastoma in children of European descent could not be confirmed, possibly due to different patterns of linkage disequilibrium or limited statistical power to detect association to variants with small effect on disease risk. Extension of GWAS to populations of African descent is important to confirm their results and validity beyond the European populations and can help to refine the location of the putative causal variants.

Phenotype restricted genome-wide association study using a gene-centric approach identifies three low-risk neuroblastoma susceptibility Loci.

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.

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.

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 approximately 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.

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 (P(allelic) = 2.35 x 10(-9)-2.25 x 10(-8)). We confirmed each SNP association in a second series of 189 high-risk cases and 1,178 controls (P(allelic) = 7.90 x 10(-7)-2.77 x 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 x 10(-18) and 2.74 x 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.

Cardiac-restricted overexpression of the A(2A)-adenosine receptor in FVB mice transiently increases contractile performance and rescues the heart failure phenotype in mice overexpressing the A(1)-adenosine receptor.

In the heart, adenosine binds to pharmacologically distinct G-protein-coupled receptors (A(1)-R, A(2A)-R, and A(3)-R). While the role of A(1)- and A(3)-Rs in the heart has been clarified, the effect of genetically manipulating the A(2A)-R has not been defined. Thus, we created mice overexpressing a cardiac-restricted A(2A)-R transgene. Mice with both low (Lo) and high (Hi) levels of A(2A)-R overexpression demonstrated an increase in cardiac contractility at 12 weeks. These changes were associated with a significantly higher systolic but not diastolic [Ca(2+)]i, higher maximal contraction amplitudes, and a significantly enhanced sarcoplasmic reticulum Ca(2+) uptake activity. At 20 weeks, the effects of A(2A)-R overexpression on cardiac contractility diminished. The positive effects elicited by A(2A)-R overexpression differ from the heart failure phenotype we observed with A(1)-R overexpression. Interestingly, coexpression of A(2A)-R TG(Hi), but not A(2A)-R TGLo, enhanced survival, prevented the development of left ventricular dysfunction and heart failure, and improved Ca(2+) handling in mice overexpressing the A(1)-R. These results suggest that adenosine-mediated signaling in the heart requires a balance between A(1)- and A(2A)-Rs--a finding that may have important implications for the ongoing clinical evaluation of adenosine receptor subtype-specific agonists and antagonists for the treatment of cardiovascular diseases.

Gene expression profiling during the transition to failure in TNF-alpha over-expressing mice demonstrates the development of autoimmune myocarditis.

Transgenic mice with cardiac-specific over-expression of tumor necrosis factor-alpha (TNF1.6) progress to dilated heart failure. A significant inflammatory response precedes functional deterioration, and may contribute to cardiac damage in this model. To evaluate the underlying molecular mechanisms, we assessed the gene expression in six groups of mouse hearts defined by age, gender, and phenotype (n = 3/group) using Affymetrix microarray analysis. Phenotype was defined as compensated (in young TNF1.6) or decompensated (in older TNF1.6) via echocardiogram. Of the >1000 transcripts altered in the compensated hearts (fold change > 2, P < 0.05 vs. wild-type (WT)), 102 were identified as immune response genes, 20 of which function in antigen presentation and processing. When comparing the compensated and decompensated hearts, >50 genes were differentially regulated, including seven immunoglobulin genes. Real-time reverse transcriptase-polymerase chain reaction and cDNA microarray confirmed the Affymetrix data. Mac3+ macrophages, CD4+ T and CD45/B220+ B-cells were identified in both compensated and decompensated hearts. However, a large amount of IgG was found deposited in areas devoid of B-lymphocytes in the myocardium of decompensated TNF1.6 mice; no such accumulation was seen in the compensated or age-matched controls. Furthermore, nuclei density analyses showed a two-fold increase in the myocardium of both compensated and decompensated TNF1.6 mice (vs. WT). This study suggests that TNF-alpha over-expression activates not only the inflammatory response, but also humoral immune responses within the transgenic hearts. The autoimmune response occurs concomitantly with cardiac decompensation and may participate in triggering the transition to failure in TNF1.6 mice.