Mutations in more than 70 genes cause hereditary anemias (HA), a highly heterogeneous group of rare/low frequency disorders in which we included: hyporegenerative anemias, as congenital dyserythropoietic anemia (CDA) and Diamond-Blackfan anemia; hemolytic anemias due to erythrocyte membrane defects, as hereditary spherocytosis and stomatocytosis; hemolytic anemias due to enzymatic defects. The study describes the diagnostic workflow for HA, based on the development of two consecutive versions of a targeted-NGS panel, including 34 and 71 genes, respectively. Seventy-four probands from 62 unrelated families were investigated. Our study includes the most comprehensive gene set for these anemias and the largest cohort of patients described so far. We obtained an overall diagnostic yield of 64.9%. Despite 54.2% of cases showed conclusive diagnosis fitting well to the clinical suspicion, the multi-gene analysis modified the original clinical diagnosis in 45.8% of patients (nonmatched phenotype-genotype). Of note, 81.8% of nonmatched patients were clinically suspected to suffer from CDA. Particularly, 45.5% of the probands originally classified as CDA exhibited a conclusive diagnosis of chronic anemia due to enzymatic defects, mainly due to mutations in PKLR gene. Interestingly, we also identified a syndromic CDA patient with mild anemia and epilepsy, showing a homozygous mutation in CAD gene, recently associated to early infantile epileptic encephalopathy-50 and CDA-like anemia. Finally, we described a patient showing marked iron overload due to the coinheritance of PIEZO1 and SEC23B mutations, demonstrating that the multi-gene approach is valuable not only for achieving a correct and definitive diagnosis, but also for guiding treatment.

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.

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.