Search

Supplementary Table 2 from 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. (Cancer Res 2006; 66(12): 6050-62)

Supplementary Table 3 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

Supplementary Table 1 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

PDF file - 189K

3016 Background: NCI-MATCH is a precision medicine trial that assigns treatment to refractory cancer patients by tumor mutation profile rather than by histology. After screening fresh tumor biopsies from nearly 6000 patients many treatment arms did not meet accrual due to the low prevalence of the eligible variants. NCI MATCH developed an approach to identify patients for the remaining arms utilizing a network of academic and commercial CLIA-certified labs that perform NGS assays as routine care at MATCH participating sites. Methods: Candidate labs were recruited through a notice in the Federal Register and posted on the NCI and ECOG ACRIN web sites. Twenty-seven labs (17 academic/10 commercial) submitted applications. After acceptance each lab analyzed a common set of 10 DNAs extracted from 8 cell lines and 2 clinical samples for concordance with the central NCI-MATCH NGS assay. Results: For the 17 labs with concordance results, a median of 8 (range 2 – 58) copy number variants (CNVs) were evaluated by the NGS assay of each DL, with the number evaluated depending on each lab’s clinical assay panel content. CNV concordance between central and DL assays, as measured by positive percent agreement (PPA), averaged 98.7% (range 87.5% - 100%) with the central assay as referent and 94.1% (range 77.8% – 100%) with the DL assay as referent. For single nucleotide variants (SNVs) and Insertion/deletions (Indels) combined, a median of 19 variants (range 11 – 26) were evaluated by each DL for concordance. PPA between central and DL assays averaged 98.0% (range 87.5% – 100%) and 98.6% (range 90.0% – 100%) with central and DL assay as referents, respectively. Strong correlations were observed between central and DL assays for both CNVs (median r = 0.93; 0.33 – 1.00) and SNV/Indels (median r = 0.98; 0.67 – 0.99). Conclusions: Our results suggest that different NGS assay platforms using diverse strategies for target enrichment and data analysis may still achieve high concordance if pre-analytical variables are minimized and the common genomic regions interrogated by each assay are well-understood. The designated lab network allows for a wider search for rare variants in tumors and provides a model for conducting future clinical trials. Clinical trial information: NCT02465060.

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase aberrantly expressed in neuroblastoma, a devastating pediatric cancer of the sympathetic nervous system. Germline and somatically acquired ALK aberrations induce increased autophosphorylation, constitutive ALK activation and increased downstream signaling. Thus, ALK is a tractable therapeutic target in neuroblastoma, likely to be susceptible to both small-molecule tyrosine kinase inhibitors and therapeutic antibodies-as has been shown for other receptor tyrosine kinases in malignancies such as breast and lung cancer. Small-molecule inhibitors of ALK are currently being studied in the clinic, but common ALK mutations in neuroblastoma appear to show de novo insensitivity, arguing that complementary therapeutic approaches must be developed. We therefore hypothesized that antibody targeting of ALK may be a relevant strategy for the majority of neuroblastoma patients likely to have ALK-positive tumors. We show here that an antagonistic ALK antibody inhibits cell growth and induces in vitro antibody-dependent cellular cytotoxicity of human neuroblastoma-derived cell lines. Cytotoxicity was induced in cell lines harboring either wild type or mutated forms of ALK. Treatment of neuroblastoma cells with the dual Met/ALK inhibitor crizotinib sensitized cells to antibody-induced growth inhibition by promoting cell surface accumulation of ALK and thus increasing the accessibility of antigen for antibody binding. These data support the concept of ALK-targeted immunotherapy as a highly promising therapeutic strategy for neuroblastomas with mutated or wild-type ALK.

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 IC 50 values of 564 nM and 548 nM, respectively. In contrast, the control lines had high micromolar IC 50 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.

Copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) are two important potential sources of phenotypic variation in humans. Until now, only SNPs have been associated with cancer, but the increasing recognition that germline DNA dosage is a critical component of human diversity raises the possibility that CNVs might also influence susceptibility to this cancer. Diskin et al. now report that a common CNV at chromosome 1q21.1 is associated with the childhood cancer neuroblastoma, and that a transcript within this CNV, the previously unknown neuroblastoma breakpoint family gene NBPF23, is involved in the early stages of tumorigenesis. Copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) are two important potential sources of phenotypic variation in humans; however, only SNPs have been associated with cancer. Here, a CNV at 1q21.1 is shown to be associated with neuroblastoma, and a transcript within this CNV, NBPF23, is implicated in early tumorigenesis of the disease. 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 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. 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) genes3,4 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.

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.

Background NCAM is a member of the immunoglobulin superfamily of cell adhesion molecules. While highly expressed on neuroblastoma cells, the relative contribution of the three major NCAM isoforms (120, 140, and 180 kDa) to neuroblastoma biology has not been investigated. Methods NCAM protein expression was measured in a neuroblastic tumor tissue microarray (N = 185) by immunohistochemistry. Relative expression of NCAM mRNA isoforms was measured in a panel of 24 human neuroblastomas and compared to fetal and adult human brain using real-time quantitative PCR and Western blot analysis. Associations with clinical and tumor biological co-variates were performed. Results NCAM protein was detected on all neuroblastic tumors and was highly expressed in all but 7/167 cases. The mRNA species predicted to encode the 120 kDa protein species was the most abundant isoform in adult brain, ganglioneuromas and ganglioneuroblastomas (P = 0.0007), but the mRNA predicted to encode the 180 kDa species was predominant in neuroblastomas (P = 0.043). Microdissected ganglion and neuroblast cells from human primary tumors confirmed these findings. Conclusion Ganglioneuromas and ganglioneuroblastomas express the adhesive 120 kDa NCAM isoform, while neuroblastomas preferentially express the 180 kDa isoform classically involved in cell motility. These data suggest a mechanism for the enhanced metastatic potential of undifferentiated neuroblastomas. Pediatr Blood Cancer 2008;51:10–16. © 2008 Wiley-Liss, Inc.

Neuroblastoma is a heterogeneous neoplasm that has served as a paradigm for the clinical utility of somatically acquired genomic aberrations. DNA copy number alterations (CNA) are currently used to predict prognosis, including MYCN amplification and deletions at chromosome bands 1p36 and 11q23. We predicted that genome-wide assessment of DNA aberrations in neuroblastoma tumors would provide a more precise estimation of clinical phenotype, and could be used to predict outcome. We measured CNAs in a representative set of 82 diagnostic tumors on a customized high-resolution BAC array-based CGH platform supplemented with additional clones across 1p36, 2p24, 3p21-22, 11q14-24, and 16p12-13, and integrated these data with RNA expression data. We used an unbiased statistical method to define a set of minimal common regions (MCRs) of aberration. Unsupervised hierarchical clustering identified four distinct genomic subclasses. First, a subset of tumors with a clinically benign phenotype showed predominantly whole chromosome gains and losses. Second, tumors with MYCN amplification had a unique genomic signature of 1p deletion and 17q gain, but few other rearrangements. Third, tumors with an aggressive clinical phenotype without MYCN amplification, showed multiple structural rearrangements. Most notable were deletions of 3p, 4p, and 11q and gain of 1q, 2p, 12q, and 17q. Lastly, there was a subset of tumors with an aggressive clinical phenotype and no detectable DNA CNAs. The genomic subsets were highly correlated with patient outcome, and individual MCRs remained prognostic in a multivariable model. DNA signature patterns embed important prognostic information in diagnostic neuroblastoma samples, and can identify candidate cancer-related genes.

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

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. (Cancer Res 2006; 66(12): 6050-62)

Neuroblastoma is a childhood cancer with considerable morbidity and mortality. Tumor-derived biomarkers may improve risk stratification.We screened 915 samples of neuroblastoma for loss of heterozygosity (LOH) at chromosome bands 1p36 and 11q23. Additional analyses identified a subgroup of cases of 11q23 LOH with unbalanced 11q LOH (unb11q LOH; defined as loss of 11q with retention of 11p). The associations of LOH with relapse and survival were determined.LOH at 1p36 was identified in 209 of 898 tumors (23 percent) and LOH at 11q23 in 307 of 913 (34 percent). Unb11q LOH was found in 151 of 307 tumors with 11q23 LOH (17 percent of the total cohort). There was a strong association of 1p36 LOH, 11q23 LOH, and unb11q LOH with most high-risk disease features (P0.001). LOH at 1p36 was associated with amplification of the MYCN oncogene (P0.001), but 11q23 LOH and unb11q LOH were not (P0.001 and P=0.002, respectively). Cases with unb11q LOH were associated with three-year event-free and overall survival rates (+/-SE) of 50+/-5 percent and 66+/-5 percent, respectively, as compared with 74+/-2 percent and 83+/-2 percent among cases without unb11q LOH (P0.001 for both comparisons). In a multivariate model, unb11q LOH was independently associated with decreased event-free survival (P=0.009) in the entire cohort, and both 1p36 LOH and unb11q LOH were independently associated with decreased progression-free survival in the subgroup of patients with features of low-risk and intermediate-risk disease (P=0.002 and P=0.02, respectively).Unb11q LOH and 1p36 LOH are independently associated with a worse outcome in patients with neuroblastoma.

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.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.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).A common genetic variation at chromosome band 6p22 is associated with susceptibility to neuroblastoma.

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. Cancer Res; 72(8); 2068–78. ©2012 AACR.

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.

A genome-wide association study (GWAS) has shown that single nucleotide variants within the LMO1 locus are associated with inherited susceptibility to neuroblastoma, a childhood cancer of the sympathetic nervous system. LMO1 encodes a transcriptional regulator previously linked to cancers. Acquired structural variation in the same locus is common in patients with neuroblastoma, suggesting that loci identified through GWAS approaches might also be prone to somatic alterations that influence tumour progression. Such studies could help to identify potential therapy targets and/or biomarkers of cancer aggressiveness. Here, single nucleotide variants within the LMO1 locus are shown to be associated with inherited susceptibility to neuroblastoma, a childhood cancer of the sympathetic nervous system. Acquired structural variation in the same locus was also frequently found in neuroblastoma patients, leading to the suggestion that loci identified through genome-wide association studies might be also prone to somatic alterations and therefore identify potential therapy targets and/or biomarkers of tumour aggressiveness. 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

PHOX2B is a homeodomain-containing protein that is involved in the development of the peripheral nervous system and is the major disease gene for the rare congenital breathing disorder congenital central hypoventilation syndrome (CCHS). Germline PHOX2B alterations were also recently discovered in neuroblastoma cases with CCHS and/or Hirschsprung disease, but a comprehensive survey for mutational frequency and functional consequence has not been performed. We therefore studied a large panel of hereditary neuroblastomas to understand the frequency and functional effects of PHOX2B mutations. Three of 47 individuals with presumed genetic predisposition to neuroblastoma showed a germline PHOX2B mutation (6.4%). Mutations were also discovered in 2 of 30 human neuroblastoma-derived cell lines, but none of 86 primary tumors from patients with sporadically occurring neuroblastoma. The vast majority of primary tumors showed abundant PHOX2B mRNA expression relative to the remainder of the transcriptome. Consistent with its role as an important neurodevelopmental gene, forced overexpression of wild-type PHOX2B in neuroblastoma cell lines suppressed cell proliferation and synergized with all-trans retinoic acid to promote differentiation. Patient-derived mutant PHOX2B constructs retained the ability to suppress cellular proliferation, but were not able to promote differentiation or activate expression of a known PHOX2B target gene in vitro. These findings show that PHOX2B alterations are a rare cause of hereditary neuroblastoma, but disruption of this neurodevelopmental pathway can interfere with transcription-dependent terminal differentiation. These data also suggest that the genetics of neuroblastoma initiation are complex, and highlight genes involved in normal noradrenergic development as candidate predisposition genes.

We used array-based comparative genomic hybridization (aCGH) to measure genomic copy number alterations (CNAs) in 42 neuroblastoma cell lines with known 1p36.3, 2p24 (MYCN), 11q23, and 17q23 allelic status. All cell lines showed CNAs, with an average of 22.0% of the genome of each sample showing evidence of gain (11.6%) or loss (10.4%). MYCN amplification was detected in 81% of cell lines, but other regions with high-level genomic amplification were observed only rarely. Gain of 17q material was present in 75% of the samples, and four discrete genomic regions at 17q23.2–17q25.3 were defined. Novel regions of gain were identified, including a 2.6-Mb subtelomeric region at 5p that includes the telomerase reverse transcriptase gene (TERT), which was found in 45% of the cell lines. Hemizygous deletions were noted at 1p36.23– 1p36.32 and 11q23.3–11q25 in 60% and 36%, respectively, of the samples, with other frequent (>25%) regions of deletion localized to 1p32.1, 3p21.31–3p22.1, 5q35.2–5q35.3, 7q31.2, 7q34, 9q22.3–9q24.1, 10q26.11–10q26.12, 16q23.1–16q24.3, 18q21.32–18q23, and 20p11.21–20p11.23. A smallest region of overlap (SRO) for CNAs was mapped across all experiments and in each case was consistent with or refined the published data. A single cell line showed a homozygous deletion at 3p22.3, which was verified, and this location was refined by FISH and PCR. There was outstanding concordance of aCGH with PCR-based CNA detection methods. Several potential cooperating loci were identified, including deletion of 11q23– 25, which was highly associated with both regional gain and loss at multiple chromosomal loci but was inversely correlated with the deletion of 1p36. Taking all of this together indicates that aCGH can accurately measure CNAs in the neuroblastoma genome and facilitate gene discovery efforts by high-throughput refinement of candidate loci. Supplementary material for this article can be found on the Genes, Chromosomes and Cancer website at http://www.interscience.wiley.com/jpages/10452257/suppmat/index.html. V V C 2005 Wiley-Liss, Inc.

Substantial genomic and functional evidence from primary tumors and cell lines indicates that a consistent region of distal chromosome 1p is deleted in a sizable proportion of human neuroblastomas, suggesting that this region contains one or more tumor suppressor genes. To determine systematically and precisely the location and extent of 1p deletion in neuroblastomas, we performed allelic loss studies of 737 primary neuroblastomas and genotype analysis of 46 neuroblastoma cell lines. Together, the results defined a single region within 1p36.3 that was consistently deleted in 25% of tumors and 87% of cell lines. Two neuroblastoma patients had constitutional deletions of distal 1p36 that overlapped the tumor-defined region. The tumor- and constitutionally-derived deletions together defined a smallest region of consistent deletion (SRD) between D1S2795 and D1S253. The 1p36.3 SRD was deleted in all but one of the 184 tumors with 1p deletion. Physical mapping and DNA sequencing determined that the SRD minimally spans an estimated 729 kb. Genomic content and sequence analysis of the SRD identified 15 characterized, nine uncharacterized, and six predicted genes in the region. The RNA expression profiles of 21 of the genes were investigated in a variety of normal tissues. The SHREW1 and KCNAB2 genes both had tissue-restricted expression patterns, including expression in the nervous system. In addition, a novel gene (CHD5) with strong homology to proteins involved in chromatin remodeling was expressed mainly in neural tissues. Together, these results suggest that one or more genes involved in neuroblastoma tumorigenesis or tumor progression are likely contained within this region.

To the Editor: We read with interest the study by Trochet and colleagues (2004), published in the April 2004 issue of The American Journal of Human Genetics, that described germline mutations of the paired-like homeobox 2B gene (PHOX2B [MIM 603851]) in neuroblastoma (MIM 256700). We have also considered PHOX2B as a candidate gene for predisposition to neuroblastoma, and we now report on a germline PHOX2B mutation in a pedigree with neuroblastoma. However, we also show that there is no evidence for mutation of this gene in eight other pedigrees with neuroblastoma screened to date. We think these data establish PHOX2B as the first bona fide gene that can predispose to neuroblastoma when mutated in the germline, and the findings further emphasize the complex genetics of this important pediatric malignancy. We previously demonstrated linkage of hereditary neuroblastoma to 16p12-13 by use of a genomewide screening strategy (Maris et al. 2002). Positional cloning of a putative 16p12-13 hereditary neuroblastoma-predisposition gene (HNB1) is ongoing, but the critical genomic region for this gene remains large. We had previously considered and excluded other genes known to be mutated in Hirschsprung disease (MIM 142623) and/or in congenital central hypoventilation syndrome (CCHS [MIM 209880]) as candidates for HNB1, because these disorders can occur coincident with both sporadic and hereditary neuroblastoma (Maris et al. 2002). Because of the recent reports that the vast majority of patients with CCHS harbor PHOX2B mutations, including two patients also affected with neuroblastoma (Amiel et al. 2003; Weese-Mayer et al. 2003), we initiated a screen for germline mutations in this gene in our series of pedigrees with neuroblastoma. Oligonucleotide primer pairs flanking the coding regions of exons 1, 2, and 3 of PHOX2B were designed by use of the program Primer 3.0; these primer pairs were used for PCR amplification and bidirectional sequencing of purified PCR products (primer sequences available on request). We screened germline DNA from the proband and an unaffected family member for each of the seven families that showed cosegregation of a 16p haplotype with disease, as well as for two pedigrees that consisted of cousins with neuroblastoma with no cosegregation of 16p marker haplotypes (see Maris et al. [2002] for details of pedigrees). We also sequenced 109 control DNA samples from the Coriell SNP500 Cancer Panel (Coriell Cell Repositories). All sequence aberrations were confirmed by repeat sequencing after cloning of purified PCR products (TOPO TA Cloning Kit [Invitrogen]), and DNA samples from the remaining available members of the pedigree were also screened for the variant. The Children’s Hospital of Philadelphia institutional review board approved this work. A heterozygous single-base deletion (676delG) was discovered in a complex pedigree with neuroblastoma (fig. 1) (see dbSNP Home Page). This family has seven members in three generations affected with neuroblastoma, and two of these individuals were also shown to have Hirschsprung disease. The proband was affected with neuroblastoma, Hirschsprung disease, and neurofibromatosis type 1 (MIM 162200). The putative nonsense mutation 676delG segregated with neuroblastoma through all three generations, and the frameshift was predicted to produce a slightly truncated protein that would no longer code for the second polyalanine tract. This family had previously been shown to cosegregate a 16p12-13 haplotype with neuroblastoma, and the proband was also shown to have an inactivating mutation in NF1 (3775delT) that was not present in either of her parents (Maris et al. 2002). Tumor material was available only for patient 1-001, and the tumor exon 3 sequence remained heterozygous for the 676delG mutation. In addition, loss-of-heterozygosity studies using microsatellite markers (D4S2912, D4S1587, D4S405, D4S2971, and D4S428) that are closely linked to the PHOX2B locus showed no evidence for allelic deletion. The only other sequence variant discovered in the remaining eight pedigrees was a putative SNP (C552T) in pedigree 12 that is not predicted to affect the resultant protein sequence (S184S) (see dbSNP Home Page). This sequence variant was present in the proband but was not detected in the patient’s affected father (no maternal DNA sample was available). It is important to note that neither sequence variant was identified in the bidirectional sequencing of 218 alleles from the control sample set. This strongly suggests that the 676delG sequence variant that segregates with the disease phenotype is a true mutation. The C552T sequence variant, which does not segregate with the disease, is more likely a very rare polymorphism, but we cannot formally exclude the possibility that there might be a functional effect of this presumably neutral polymorphism. Figure 1 Germline PHOX2B mutation in a pedigree segregating neuroblastoma and Hirschsprung disease. A, Family 1 pedigree structure. DNA samples from this family with neuroblastoma were available only for patients with a PHOX2B result shown (wt = wild type; 676delG ... Accumulated data strongly implicate PHOX2B as an essential regulator of normal autonomic nervous system development (Pattyn et al. 1999; Brunet and Pattyn 2002). The discovery of polyalanine-expansion mutations in the majority of patients with CCHS clearly defines a role for this gene in human disease (Amiel et al. 2003; Weese-Mayer et al. 2003), and there appears to be a correlation between the severity of the respiratory symptoms and the length of polyalanine expansion (Weese-Mayer et al. 2003; Matera et al. 2004). Neuroblastoma represents perhaps the most aberrant phenotype that results from abnormal adrenergic tissue development. The rare but well-described synchronous appearance of neuroblastoma with other disorders of the autonomic nervous system has suggested a common genetic etiology often referred to as a “neurocristopathy” (Gaisie et al. 1979; Nemecek et al. 2003). Although other genes implicated in Hirschsprung disease and/or CCHS have not been excluded as hereditary neuroblastoma-predisposition genes (Maris et al. 2002; Perri et al. 2002), our data further establish PHOX2B as an important gene involved in the initiation of neuroblastoma tumorigenesis. However, the fact that the majority of pedigrees studied here do not show PHOX2B mutations clearly implicates locus heterogeneity for hereditary predisposition to neuroblastoma. Assuming that our inferences of linkage to 16p are correct, and in light of the observation of two germline mutations in the proband of the family presented here, we suggest that an oligogenic mechanism for neuroblastoma initiation should be considered, as has been shown for other diseases of neural crest–derived tissues (Gabriel et al. 2002). It is not yet clear if the PHOX2B mutations discovered in patients with hereditary or sporadic neuroblastoma result in gain or loss of protein function. The hypothesis that PHOX2B functions as a tumor suppressor is supported by the potential predicted consequence of the five mutations described, to date, in patients with neuroblastoma. Weese-Mayer and colleagues discovered a nonsense mutation that predicts a significantly truncated protein that would miss most of exon 3, including all of the 20-alanine repeat motif (Weese-Mayer et al. 2003). The frameshift mutation described here is similar to that reported by Amiel and colleagues (2003) in a patient who also had CCHS, Hirschsprung disease, and neuroblastoma, and, in both cases, the changes in reading frame are predicted to abolish the polyalanine tract. Trochet and colleagues (2004) discovered two missense mutations, both of which map to a conserved portion of the homeodomain and thus may interfere with DNA binding. On the other hand, 4p12 is not a known site of frequent allelic deletion in neuroblastoma (Maris and Matthay 1999), and, to date, biallelic inactivation of PHOX2B has not been demonstrated, although far too few cases have been examined to assert this with confidence. Taken together, these data suggest that PHOX2B mutations are involved in the initiation of neuroblastoma tumorigenesis, especially in patients with associated disorders of the autonomic nervous system. Our data also indicate that germline mutational events in this gene are not involved in the majority of hereditary neuroblastoma cases and that alternative genetic events may predispose to tumorigenesis. Examination of additional patients will facilitate the definition of PHOX2B mutation frequency in the genetic and (apparently) sporadic forms of neuroblastoma and will help to clarify the role of PHOX2B mutations in tumor initiation and progression.

Regional genomic alterations resulting from single-copy allelic loss or gain have been well characterized in many human cancers and are often of prognostic relevance. Unbalanced gain of 17q material is common in malignant human neuroblastomas and typically results from unbalanced translocations. Unbalanced 17q gain may be an independent predictor of disease outcome, but technical difficulties with quantifying such gain using fluorescent in situ hybridization gives this method limited clinical applicability. We now describe a duplex genomic DNA-based quantitative polymerase chain reaction assay to determine the presence or absence of unbalanced gain of chromosome 17q in primary neuroblastoma specimens. The technique was first refined and validated in a panel of nine human neuroblastoma-derived cell lines by direct comparison with dual-color fluorescent in situ hybridization. Prospective blinded comparison of quantitative polymerase chain reaction and fluorescence in situ hybridization in 40 human neuroblastoma primary tumor samples showed a sensitivity of 96% and 100% specificity for detecting unbalanced 17q gain when a relative 17q copy number ratio of 1.3 was used to define unbalanced gain. Tumors with ratios >1.3 were highly associated with malignant tumor phenotypic features such as metastatic disease (P

Genome-wide association studies have proven to be a powerful tool to identify susceptibility loci in cancer, however the molecular basis of a majority of these association signals remain undefined. We have previously reported nine common single nucleotide polymorphism (SNP) alleles at chromosome 2q35, within and 5′ upstream of the BRCA1-associated RING domain (BARD1) gene, to be significantly associated with high-risk neuroblastoma (Capasso, M et al., Nature Genetics 2009). Here, we focus on understanding the biological underpinnings of this genetic association by investigating how common variation at this locus correlates with changes in BARD1 expression in both the germline and in neuroblastoma samples, and further how these BARD1 alterations are involved in driving neuroblastoma tumorigenesis. We first showed that the neuroblastoma associated SNP variation at 2q35 is in linkage disequilibrium (LD) with multiple genomic areas dense with regulatory domains, including many putative binding sites for the neuronal splicing regulator, FOX2. We next demonstrated that in peripheral blood mononuclear cells, these SNPs correlate both with a trend towards increased overall BARD1 expression and with differential BARD1 exon expression. Specifically, a majority of the neuroblastoma associated SNPs were significantly associated with decreased expression of the RING domain forming BARD1 exons 2 and 3 (P = 0.0008-0.034). We have identified 16 unique BARD1 isoforms in the neuroblastoma tissue of origin, sympathetic fetal ganglia, with the β-BARD1 isoform specifically splicing out exons 2 and 3. β-BARD1 was found to be ubiquitously expressed in a panel of fetal sympathetic ganglia, neuroblastoma cell lines, and high-risk primary tumors, with its expression being significantly increased in cell lines harboring a homozygous BARD1 risk allele genotype at the most significantly associated 2q35 SNP, rs6435862 (N = 22; P = 0.04). Transient transfection of a specific siRNA targeting the alternatively spliced boundary of β-BARD1 in a panel of neuroblastoma cell lines (N = 12) showed that neuroblastoma cells harboring at least one risk allele at SNP rs6435862 and the highest β-BARD1 expression showed a significant inhibition of cell proliferation upon β-BARD1 knockdown (% growth inhibition = 45-81%, P = 0.0004-0.034), while six of seven cell lines with a homozygous major allele genotype at rs6435862 and lower β-BARD1 expression did not show significant growth inhibition. These data are consistent with β-BARD1 acting as an oncogene in neuroblastoma and we hypothesize that common genomic variation drives formation of this isoform, underlying neuroblastoma predisposition on 2q35. Ongoing work will focus on further defining the role of common variation in BARD1's RNA splicing patterns and identifying the molecular basis of β-BARD1's oncogenicity in neuroblastoma. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 3867.

6 Background: Chromosome arm 11q loss of heterozygosity (LOH) is a frequent occurrence in neuroblastoma. Previous studies have shown an association between 11q LOH and high-risk disease. However, these studies lacked the statistical power to conclusively analyze the prognostic impact of 11q LOH in multivariable analyses after adjusting for the currently used prognostic factors. Methods: We screened 917 primary neuroblastomas for chromosome 11 LOH with a panel of 4–8 microsatellite markers. Further mapping with up to 52 markers was performed on selected samples. The only inclusion criteria were availability of a constitutional DNA sample and outcome data. The sample set was representative of the overall neuroblastoma population. The median follow-up was 2.7 years. Results: LOH at chromosome arm 11q and unbalanced 11q LOH (loss of 11q material with retention or gain of 11p material) were present in 302 (33%) and 148 (16%) cases, respectively. The pattern of LOH for each tumor sample was consistent with the ...

The vitamin D receptor has been found in several human organs not involved in calcium metabolism and in several malignant neoplasms found in humans. The role of the receptor in these tissues is unclear. There is, however, a relationship between the presence and quantity of the vitamin D receptor in a malignant cell line and the antineoplastic effect of vitamin D on that cell line. We found that Y-79 retinoblastoma cells have receptors specific for calcitriol (1,25-dihydroxycholecalciferol). Scatchard analysis of the receptor data shows a quantity of 56,000 receptors per retinoblastoma cell. These receptors have a dissociation constant of 1.18 nmol/L. Retinoblastoma cells treated with 10(-9) mol/L of calcitriol for nine days had 15% less cell growth than the control cells. Further studies of the effect of vitamin D on retinoblastoma may warrant its inclusion in chemotherapeutic protocols for the treatment of this childhood affliction.