The spectrum of somatic mutation of the most aggressive forms of neuroblastoma is not completely determined. We sought to identify potential cancer drivers in clinically aggressive neuroblastoma.Whole exome sequencing was conducted on 17 germline and tumor DNA samples from high-risk patients with adverse events within 36 months from diagnosis (HR-Event3) to identify somatic mutations and deep targeted sequencing of 134 genes selected from the initial screening in additional 48 germline and tumor pairs (62.5% HR-Event3 and high-risk patients), 17 HR-Event3 tumors and 17 human-derived neuroblastoma cell lines.We revealed 22 significantly mutated genes, many of which implicated in cancer progression. Fifteen genes (68.2%) were highly expressed in neuroblastoma supporting their involvement in the disease. CHD9, a cancer driver gene, was the most significantly altered (4.0% of cases) after ALK.Other genes (PTK2, NAV3, NAV1, FZD1 and ATRX), expressed in neuroblastoma and involved in cell invasion and migration were mutated at frequency ranged from 4% to 2%.Focal adhesion and regulation of actin cytoskeleton pathways, were frequently disrupted (14.1% of cases) thus suggesting potential novel therapeutic strategies to prevent disease progression.Notably BARD1, CHEK2 and AXIN2 were enriched in rare, potentially pathogenic, germline variants.In summary, whole exome and deep targeted sequencing identified novel cancer genes of clinically aggressive neuroblastoma. Our analyses show pathway-level implications of infrequently mutated genes in leading neuroblastoma progression.

Paget disease of bone (PDB) is a skeletal disorder characterized by focal abnormalities of bone remodeling, which result in enlarged and deformed bones in one or more regions of the skeleton. In some cases, the pagetic tissue undergoes neoplastic transformation, resulting in osteosarcoma and, less frequently, in giant cell tumor of bone (GCT). We performed whole-exome sequencing in a large family with 14 PDB-affected members, four of whom developed GCT at multiple pagetic skeletal sites, and we identified the c.2810C>G (p.Pro937Arg) missense mutation in the zinc finger protein 687 gene (ZNF687). The mutation precisely co-segregated with the clinical phenotype in all affected family members. The sequencing of seven unrelated individuals with GCT associated with PDB (GCT/PDB) identified the same mutation in all individuals, unravelling a founder effect. ZNF687 is highly expressed during osteoclastogenesis and osteoblastogenesis and is dramatically upregulated in the tumor tissue of individuals with GCT/PDB. Interestingly, our preliminary findings showed that ZNF687, indicated as a target gene of the NFkB transcription factor by ChIP-seq analysis, is also upregulated in the peripheral blood of PDB-affected individuals with (n = 5) or without (n = 6) mutations in SQSTM1, encouraging additional studies to investigate its potential role as a biomarker of PDB risk.

Piwi-associated RNAs (piRNAs) are a distinct class of 24- to 30-nucleotide-long RNAs produced by a Dicer-independent mechanism, and are associated with Piwi-class Argonaute proteins. In contrast to the several hundred species of microRNAs (miRNAs) identified thus far, piRNAs consist of more than 30,000 different species in humans. Studies in flies, fish and mice implicate these piRNAs in regulating germ line development, the silencing of selfish DNA elements, and maintaining germ line DNA integrity. Most piRNAs map to unique sites in the human genome, including intergenic, intronic, and exonic sequences. However, the role of piRNAs in humans remains to be elucidated. Here, we uncover an unexpected function of the piRNA pathway in humans. We show for the first time, that the piRNA_015520, located in intron 1 of the human Melatonin receptor 1A (MTNR1A) gene, is expressed in adult human tissues (testes and brain) and in the human cell line HEK 293. Although the role of piR_015520 expression in brain tissue remains unknown, the testes-specific expression is consistent with previous findings in several species.Surprisingly, in contrast to the mechanism known for miRNA-mediated modulation of gene expression, piRNA_015520 negatively regulates MTNR1A gene expression by binding to its genomic region. This finding suggests that changes in individual piRNA levels could influence both autoregulatory gene expression and the expression of the gene in which the piRNA is located. These findings offer a new perspective for piRNAs functioning as gene regulators in humans.

BRCA1 Associated RING Domain 1 (BARD1) encodes a protein which interacts with the N-terminal region of BRCA1 in vivo and in vitro. The full length (FL) BARD1 mRNA includes 11 exons and encodes a protein comprising of six domains (N-terminal RING-finger domain, three Ankyrin repeats and two C-terminal BRCT domains) with different functions. Emerging data suggest that BARD1 can have both tumor-suppressor gene and oncogene functions in tumor initiation and progression. Indeed, whereas FL BARD1 protein acts as tumor-suppressor with and without BRCA1 interactions, aberrant splice variants of BARD1 have been detected in various cancers and have been shown to play an oncogenic role. Further evidence for a dualistic role came with the identification of BARD1 as a neuroblastoma predisposition gene in our genome wide association study which has demonstrated that single nucleotide polymorphisms in BARD1 can correlate with risk or can protect against cancer based on their association with the expression of FL and splice variants of BARD1. This review is an overview of how BARD1 functions in tumorigenesis with opposite effects in various types of cancer.

A previous genome-wide association study (GWAS) identified common variation at the BARD1 locus as being highly associated with susceptibility to high-risk neuroblastoma, but the mechanisms underlying this association have been not extensively investigated. Here, we performed a fine mapping analysis of BARD1 locus (2q35) using GWAS data from 556 high-risk neuroblastoma patients and 2,575 controls of European-American ancestry, and identified two independent genome-wide neuroblastoma-associated loci. Functional single-nucleotide polymorphism (SNP) prioritization identified two causative variants that independently contributed to neuroblastoma risk, and each replicated robustly in multiple independent cohorts comprising 445 high-risk cases and 3,170 controls (rs17489363: combined p = 1.07 × 10 , OR:1.79, 95% CI:1.62-1.98 and rs1048108: combined p = 7.27 × 10 , OR:0.65, 95% CI:0.58-0.73). Particularly, the T risk allele of rs17489363 in the canonical promoter region of full-length BARD1 altered binding site of the transcription factor HSF1 and correlated with low expression of full-length BARD1 mRNA and protein. Low-level expression of full-length BARD1 associated with advanced neuroblastoma. In human neuroblastoma cells, attenuating full-length BARD1 increased proliferation and invasion capacity. In conclusion, we have identified two potentially causative SNPs at the BARD1 locus associated with predisposition to high-risk neuroblastoma, and have shown that full-length BARD1 may act as tumor suppressor.

The genetic aetiology of sporadic neuroblastoma is still largely unknown. We have identified diverse neuroblastoma susceptibility loci by genomewide association studies (GWASs); however, additional SNPs that likely contribute to neuroblastoma susceptibility prompted this investigation for identification of additional variants that are likely hidden among signals discarded by the multiple testing corrections used in the analysis of genomewide data. There is evidence suggesting the CDKN1B, coding for the cycle inhibitor p27Kip1, is involved in neuroblastoma. We thus assess whether genetic variants of CDKN1B are associated with neuroblastoma. We imputed all possible genotypes across CDKN1B locus on a discovery case series of 2101 neuroblastoma patients and 4202 genetically matched controls of European ancestry. The most significantly associated rs34330 was analysed in an independent Italian cohort of 311 cases and 709 controls. In vitro functional analysis was carried out in HEK293T and in neuroblastoma cell line SHEP‐2, both transfected with pGL3‐CDKN1B‐CC or pGL3‐CDKN1B‐TT constructs. We identified an association of the rs34330 T allele (‐79C/T) with the neuroblastoma risk (Pcombined = 0.002; OR = 1.17). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27Kip1 promoter (P < 0.05). Three independent sets of neuroblastoma tumours carrying ‐79TT genotype showed a tendency towards lower CDKN1B mRNA levels. Our study shows that a functional variant, associated with a reduced CDKN1B gene transcription, influences neuroblastoma susceptibility.

BackgroundGlutamate is the principal excitatory neurotransmitter in the central nervous system which acts by the activation of either ionotropic (AMPA, NMDA and kainate receptors) or G-protein coupled metabotropic receptors. Glutamate is widely accepted to play a major role in the path physiology of migraine as implicated by data from animal and human studies. Genes involved in synthesis, metabolism and regulation of both glutamate and its receptors could be, therefore, considered as potential candidates for causing/predisposing to migraine when mutated.MethodsThe association of polymorphic variants of GRIA1-GRIA4 genes which encode for the four subunits (GluR1-GluR4) of the alpha-amino-3- hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor for glutamate was tested in migraineurs with and without aura (MA and MO) and healthy controls.ResultsTwo variants in the regulative regions of GRIA1 (rs2195450) and GRIA3 (rs3761555) genes resulted strongly associated with MA (P = 0.00002 and P = 0.0001, respectively), but not associated with MO, suggesting their role in cortical spreading depression. Whereas the rs548294 variant in GRIA1 gene showed association primarily with MO phenotype, supporting the hypothesis that MA and MO phenotypes could be genetically related. These variants modify binding sites for transcription factors altering the expression of GRIA1 and GRIA3 genes in different conditions.ConclusionsThis study represents the first genetic evidence of a link between glutamate receptors and migraine.

BackgroundThe prognosis of children with metastatic stage 4 neuroblastoma (NB) has remained poor in the past decade.Patients and methodsUsing microarray analyses of 342 primary tumors, we here developed and validated an easy to use gene expression-based risk score including 18 genes, which can robustly predict the outcome of stage 4 patients.ResultsThis classifier was a significant predictor of overall survival in two independent validation cohorts [cohort 1 (n = 214): P = 6.3 × 10−5; cohort 2 (n = 27): P = 3.1 × 10−2]. The prognostic value of the risk score was validated by multivariate analysis including the established markers age and MYCN status (P = 0.027). In the pooled validation cohorts (n = 241), integration of the risk score with the age and/or MYCN status identified subgroups with significantly differing overall survival (ranging from 35 to 100 %).ConclusionTogether, the 18-gene risk score classifier can identify patients with stage 4 NB with favorable outcome and may therefore improve risk assessment and treatment stratification of NB patients with disseminated disease.Electronic supplementary materialThe online version of this article (doi:10.1186/s12967-016-0896-7) contains supplementary material, which is available to authorized users.