Clonal dynamics in pediatric B-cell precursor acute lymphoblastic leukemia with very early relapse.

INTRODUCTION: One-quarter of the relapses in children with B-cell precursor acute lymphoblastic leukemia (BCP-ALL) occur very early (within 18 months, before completion of treatment), and prognosis in these patients is worse compared to cases that relapse after treatment has ended. METHODS: In this study, we performed a genomic analysis of diagnosis-relapse pairs of 12 children who relapsed very early, followed by a deep-sequencing validation of all identified mutations. In addition, we included one case with a good initial treatment response and on-treatment relapse at the end of upfront therapy. RESULTS: We observed a dynamic clonal evolution in all cases, with relapse almost exclusively originating from a subclone at diagnosis. We identified several driver mutations that may have influenced the outgrowth of a minor clone at diagnosis to become the major clone at relapse. For example, a minimal residual disease (MRD)-based standard-risk patient with ETV6-RUNX1-positive leukemia developed a relapse from a TP53-mutated subclone after loss of the wildtype allele. Furthermore, two patients with TCF3-PBX1-positive leukemia that developed a very early relapse carried E1099K WHSC1 mutations at diagnosis, a hotspot mutation that was recurrently encountered in other very early TCF3-PBX1-positive leukemia relapses as well. In addition to alterations in known relapse drivers, we found two cases with truncating mutations in the cohesin gene RAD21. CONCLUSION: Comprehensive genomic characterization of diagnosis-relapse pairs shows that very early relapses in BCP-ALL frequently arise from minor subclones at diagnosis. A detailed understanding of the therapeutic pressure driving these events may aid the development of improved therapies.

Unravelling the Sequential Interplay of Mutational Mechanisms during Clonal Evolution in Relapsed Pediatric Acute Lymphoblastic Leukemia.

Pediatric acute lymphoblastic leukemia (ALL) is the most common pediatric malignancy and is characterized by clonal heterogeneity. Genomic mutations can increase proliferative potential of leukemic cells and cause treatment resistance. However, mechanisms driving mutagenesis and clonal diversification in ALL are not fully understood. In this proof of principle study, we performed whole genome sequencing of two cases with multiple relapses in order to investigate whether groups of mutations separated in time show distinct mutational signatures. Based on mutation allele frequencies at diagnosis and subsequent relapses, we clustered mutations into groups and performed cluster-specific mutational profile analysis and de novo signature extraction. In patient 1, who experienced two relapses, the analysis unraveled a continuous interplay of aberrant activation induced cytidine deaminase (AID)/apolipoprotein B editing complex (APOBEC) activity. The associated signatures SBS2 and SBS13 were present already at diagnosis, and although emerging mutations were lost in later relapses, the process remained active throughout disease evolution. Patient 2 had three relapses. We identified episodic mutational processes at diagnosis and first relapse leading to mutations resembling ultraviolet light-driven DNA damage, and thiopurine-associated damage at first relapse. In conclusion, our data shows that investigation of mutational processes in clusters separated in time may aid in understanding the mutational mechanisms and discovery of underlying causes.

Evidence for 28 genetic disorders discovered by combining healthcare and research data.

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

Rare genetic variants in interleukin-37 link this anti-inflammatory cytokine to the pathogenesis and treatment of gout.

OBJECTIVE: Gout is characterised by severe interleukin (IL)-1-mediated joint inflammation induced by monosodium urate crystals. Since IL-37 is a pivotal anti-inflammatory cytokine suppressing the activity of IL-1, we conducted genetic and functional studies aimed at elucidating the role of IL-37 in the pathogenesis and treatment of gout. METHODS: Variant identification was performed by DNA sequencing of all coding bases of IL37 using molecular inversion probe-based resequencing (discovery cohort: gout n=675, controls n=520) and TaqMan genotyping (validation cohort: gout n=2202, controls n=2295). Predictive modelling of the effects of rare variants on protein structure was followed by in vitro experiments evaluating the impact on protein function. Treatment with recombinant IL-37 was evaluated in vitro and in vivo in a mouse model of gout. RESULTS: We identified four rare variants in IL37 in six of the discovery gout patients; p.(A144P), p.(G174Dfs*16), p.(C181*) and p.(N182S), whereas none emerged in healthy controls (Fisher's exact p-value=0.043). All variants clustered in the functional domain of IL-37 in exon 5 (p-value=5.71×10-5). Predictive modelling and functional studies confirmed loss of anti-inflammatory functions and we substantiated the therapeutic potential of recombinant IL-37 in the treatment of gouty inflammation. Furthermore, the carrier status of p.(N182S)(rs752113534) was associated with increased risk (OR=1.81, p-value=0.031) of developing gout in hyperuricaemic individuals of Polynesian ancestry. CONCLUSION: Here, we provide genetic as well as mechanistic evidence for the role of IL-37 in the pathogenesis of gout, and highlight the therapeutic potential of recombinant IL-37 for the treatment of gouty arthritis.

Exome sequencing in routine diagnostics: a generic test for 254 patients with primary immunodeficiencies.

BACKGROUND: Diagnosis of primary immunodeficiencies (PIDs) is complex and cumbersome yet important for the clinical management of the disease. Exome sequencing may provide a genetic diagnosis in a significant number of patients in a single genetic test. METHODS: In May 2013, we implemented exome sequencing in routine diagnostics for patients suffering from PIDs. This study reports the clinical utility and diagnostic yield for a heterogeneous group of 254 consecutively referred PID patients from 249 families. For the majority of patients, the clinical diagnosis was based on clinical criteria including rare and/or unusual severe bacterial, viral, or fungal infections, sometimes accompanied by autoimmune manifestations. Functional immune defects were interpreted in the context of aberrant immune cell populations, aberrant antibody levels, or combinations of these factors. RESULTS: For 62 patients (24%), exome sequencing identified pathogenic variants in well-established PID genes. An exome-wide analysis diagnosed 10 additional patients (4%), providing diagnoses for 72 patients (28%) from 68 families altogether. The genetic diagnosis directly indicated novel treatment options for 25 patients that received a diagnosis (34%). CONCLUSION: Exome sequencing as a first-tier test for PIDs granted a diagnosis for 28% of patients. Importantly, molecularly defined diagnoses indicated altered therapeutic options in 34% of cases. In addition, exome sequencing harbors advantages over gene panels as a truly generic test for all genetic diseases, including in silico extension of existing gene lists and re-analysis of existing data.

Somatic variants in autosomal dominant genes are a rare cause of sporadic Alzheimer's disease.

INTRODUCTION: A minority of patients with sporadic early-onset Alzheimer's disease (AD) exhibit de novo germ line mutations in the autosomal dominant genes such as APP, PSEN1, or PSEN2. We hypothesized that negatively screened patients may harbor somatic variants in these genes. METHODS: We applied an ultrasensitive approach based on single-molecule molecular inversion probes followed by deep next generation sequencing of 11 genes to 100 brain and 355 blood samples from 445 sporadic patients with AD (>80% exhibited an early onset, <66 years). RESULTS: We identified and confirmed nine somatic variants (allele fractions: 0.2%-10.8%): two APP, five SORL1, one NCSTN, and one MARK4 variants by independent amplicon-based deep sequencing. DISCUSSION: Two of the SORL1 variant might have contributed to the disease, the two APP variants were interpreted as likely benign and the other variants remained of unknown significance. Somatic variants in the autosomal dominant AD genes may not be a common cause of sporadic AD, including early onset cases.

MPZL2, Encoding the Epithelial Junctional Protein Myelin Protein Zero-like 2, Is Essential for Hearing in Man and Mouse.

In a Dutch consanguineous family with recessively inherited nonsyndromic hearing impairment (HI), homozygosity mapping combined with whole-exome sequencing revealed a MPZL2 homozygous truncating variant, c.72del (p.Ile24Metfs∗22). By screening a cohort of phenotype-matched subjects and a cohort of HI subjects in whom WES had been performed previously, we identified two additional families with biallelic truncating variants of MPZL2. Affected individuals demonstrated symmetric, progressive, mild to moderate sensorineural HI. Onset of HI was in the first decade, and high-frequency hearing was more severely affected. There was no vestibular involvement. MPZL2 encodes myelin protein zero-like 2, an adhesion molecule that mediates epithelial cell-cell interactions in several (developing) tissues. Involvement of MPZL2 in hearing was confirmed by audiometric evaluation of Mpzl2-mutant mice. These displayed early-onset progressive sensorineural HI that was more pronounced in the high frequencies. Histological analysis of adult mutant mice demonstrated an altered organization of outer hair cells and supporting cells and degeneration of the organ of Corti. In addition, we observed mild degeneration of spiral ganglion neurons, and this degeneration was most pronounced at the cochlear base. Although MPZL2 is known to function in cell adhesion in several tissues, no phenotypes other than HI were found to be associated with MPZL2 defects. This indicates that MPZL2 has a unique function in the inner ear. The present study suggests that deleterious variants of Mplz2/MPZL2 affect adhesion of the inner-ear epithelium and result in loss of structural integrity of the organ of Corti and progressive degeneration of hair cells, supporting cells, and spiral ganglion neurons.

Heterozygous missense variants of LMX1A lead to nonsyndromic hearing impairment and vestibular dysfunction.

Unraveling the causes and pathomechanisms of progressive disorders is essential for the development of therapeutic strategies. Here, we identified heterozygous pathogenic missense variants of LMX1A in two families of Dutch origin with progressive nonsyndromic hearing impairment (HI), using whole exome sequencing. One variant, c.721G > C (p.Val241Leu), occurred de novo and is predicted to affect the homeodomain of LMX1A, which is essential for DNA binding. The second variant, c.290G > C (p.Cys97Ser), predicted to affect a zinc-binding residue of the second LIM domain that is involved in protein-protein interactions. Bi-allelic deleterious variants of Lmx1a are associated with a complex phenotype in mice, including deafness and vestibular defects, due to arrest of inner ear development. Although Lmx1a mouse mutants demonstrate neurological, skeletal, pigmentation and reproductive system abnormalities, no syndromic features were present in the participating subjects of either family. LMX1A has previously been suggested as a candidate gene for intellectual disability, but our data do not support this, as affected subjects displayed normal cognition. Large variability was observed in the age of onset (a)symmetry, severity and progression rate of HI. About half of the affected individuals displayed vestibular dysfunction and experienced symptoms thereof. The late-onset progressive phenotype and the absence of cochleovestibular malformations on computed tomography scans indicate that heterozygous defects of LMX1A do not result in severe developmental abnormalities in humans. We propose that a single LMX1A wild-type copy is sufficient for normal development but insufficient for maintenance of cochleovestibular function. Alternatively, minor cochleovestibular developmental abnormalities could eventually lead to the progressive phenotype seen in the families.

A Recurrent De Novo PACS2 Heterozygous Missense Variant Causes Neonatal-Onset Developmental Epileptic Encephalopathy, Facial Dysmorphism, and Cerebellar Dysgenesis.

Developmental and epileptic encephalopathies (DEEs) represent a large clinical and genetic heterogeneous group of neurodevelopmental diseases. The identification of pathogenic genetic variants in DEEs remains crucial for deciphering this complex group and for accurately caring for affected individuals (clinical diagnosis, genetic counseling, impacting medical, precision therapy, clinical trials, etc.). Whole-exome sequencing and intensive data sharing identified a recurrent de novo PACS2 heterozygous missense variant in 14 unrelated individuals. Their phenotype was characterized by epilepsy, global developmental delay with or without autism, common cerebellar dysgenesis, and facial dysmorphism. Mixed focal and generalized epilepsy occurred in the neonatal period, controlled with difficulty in the first year, but many improved in early childhood. PACS2 is an important PACS1 paralog and encodes a multifunctional sorting protein involved in nuclear gene expression and pathway traffic regulation. Both proteins harbor cargo(furin)-binding regions (FBRs) that bind cargo proteins, sorting adaptors, and cellular kinase. Compared to the defined PACS1 recurrent variant series, individuals with PACS2 variant have more consistently neonatal/early-infantile-onset epilepsy that can be challenging to control. Cerebellar abnormalities may be similar but PACS2 individuals exhibit a pattern of clear dysgenesis ranging from mild to severe. Functional studies demonstrated that the PACS2 recurrent variant reduces the ability of the predicted autoregulatory domain to modulate the interaction between the PACS2 FBR and client proteins, which may disturb cellular function. These findings support the causality of this recurrent de novo PACS2 heterozygous missense in DEEs with facial dysmorphim and cerebellar dysgenesis.

Spatial Clustering of de Novo Missense Mutations Identifies Candidate Neurodevelopmental Disorder-Associated Genes.

Haploinsufficiency (HI) is the best characterized mechanism through which dominant mutations exert their effect and cause disease. Non-haploinsufficiency (NHI) mechanisms, such as gain-of-function and dominant-negative mechanisms, are often characterized by the spatial clustering of mutations, thereby affecting only particular regions or base pairs of a gene. Variants leading to haploinsufficency might occasionally cluster as well, for example in critical domains, but such clustering is on the whole less pronounced with mutations often spread throughout the gene. Here we exploit this property and develop a method to specifically identify genes with significant spatial clustering patterns of de novo mutations in large cohorts. We apply our method to a dataset of 4,061 de novo missense mutations from published exome studies of trios with intellectual disability and developmental disorders (ID/DD) and successfully identify 15 genes with clustering mutations, including 12 genes for which mutations are known to cause neurodevelopmental disorders. For 11 out of these 12, NHI mutation mechanisms have been reported. Additionally, we identify three candidate ID/DD-associated genes of which two have an established role in neuronal processes. We further observe a higher intolerance to normal genetic variation of the identified genes compared to known genes for which mutations lead to HI. Finally, 3D modeling of these mutations on their protein structures shows that 81% of the observed mutations are unlikely to affect the overall structural integrity and that they therefore most likely act through a mechanism other than HI.

YY1 Haploinsufficiency Causes an Intellectual Disability Syndrome Featuring Transcriptional and Chromatin Dysfunction.

Yin and yang 1 (YY1) is a well-known zinc-finger transcription factor with crucial roles in normal development and malignancy. YY1 acts both as a repressor and as an activator of gene expression. We have identified 23 individuals with de novo mutations or deletions of YY1 and phenotypic features that define a syndrome of cognitive impairment, behavioral alterations, intrauterine growth restriction, feeding problems, and various congenital malformations. Our combined clinical and molecular data define "YY1 syndrome" as a haploinsufficiency syndrome. Through immunoprecipitation of YY1-bound chromatin from affected individuals' cells with antibodies recognizing both ends of the protein, we show that YY1 deletions and missense mutations lead to a global loss of YY1 binding with a preferential retention at high-occupancy sites. Finally, we uncover a widespread loss of H3K27 acetylation in particular on the YY1-bound enhancers, underscoring a crucial role for YY1 in enhancer regulation. Collectively, these results define a clinical syndrome caused by haploinsufficiency of YY1 through dysregulation of key transcriptional regulators.

Diagnostic exome sequencing in 266 Dutch patients with visual impairment.

Inherited eye disorders have a large clinical and genetic heterogeneity, which makes genetic diagnosis cumbersome. An exome-sequencing approach was developed in which data analysis was divided into two steps: the vision gene panel and exome analysis. In the vision gene panel analysis, variants in genes known to cause inherited eye disorders were assessed for pathogenicity. If no causative variants were detected and when the patient consented, the entire exome data was analyzed. A total of 266 Dutch patients with different types of inherited eye disorders, including inherited retinal dystrophies, cataract, developmental eye disorders and optic atrophy, were investigated. In the vision gene panel analysis (likely), causative variants were detected in 49% and in the exome analysis in an additional 2% of the patients. The highest detection rate of (likely) causative variants was in patients with inherited retinal dystrophies, for instance a yield of 63% in patients with retinitis pigmentosa. In patients with developmental eye defects, cataract and optic atrophy, the detection rate was 50, 33 and 17%, respectively. An exome-sequencing approach enables a genetic diagnosis in patients with different types of inherited eye disorders using one test. The exome approach has the same detection rate as targeted panel sequencing tests, but offers a number of advantages. For instance, the vision gene panel can be frequently and easily updated with additional (novel) eye disorder genes. Determination of the genetic diagnosis improved the clinical diagnosis, regarding the assessment of the inheritance pattern as well as future disease perspective.

BRCA Testing by Single-Molecule Molecular Inversion Probes.

BACKGROUND: Despite advances in next generation DNA sequencing (NGS), NGS-based single gene tests for diagnostic purposes require improvements in terms of completeness, quality, speed, and cost. Single-molecule molecular inversion probes (smMIPs) are a technology with unrealized potential in the area of clinical genetic testing. In this proof-of-concept study, we selected 2 frequently requested gene tests, those for the breast cancer genes BRCA1 and BRCA2, and developed an automated work flow based on smMIPs. METHODS: The BRCA1 and BRCA2 smMIPs were validated using 166 human genomic DNA samples with known variant status. A generic automated work flow was built to perform smMIP-based enrichment and sequencing for BRCA1, BRCA2, and the checkpoint kinase 2 (CHEK2) c.1100del variant. RESULTS: Pathogenic and benign variants were analyzed in a subset of 152 previously BRCA-genotyped samples, yielding an analytical sensitivity and specificity of 100%. Following automation, blind analysis of 65 in-house samples and 267 Norwegian samples correctly identified all true-positive variants (>3000), with no false positives. Consequent to process optimization, turnaround times were reduced by 60% to currently 10-15 days. Copy number variants were detected with an analytical sensitivity of 100% and an analytical specificity of 88%. CONCLUSIONS: smMIP-based genetic testing enables automated and reliable analysis of the coding sequences of BRCA1 and BRCA2. The use of single-molecule tags, double-tiled targeted enrichment, and capturing and sequencing in duplo, in combination with automated library preparation and data analysis, results in a robust process and reduces routine turnaround times. Furthermore, smMIP-based copy number variation analysis could make independent copy number variation tools like multiplex ligation-dependent probes amplification dispensable.

The diagnostic yield of whole-exome sequencing targeting a gene panel for hearing impairment in The Netherlands.

Hearing impairment (HI) is genetically heterogeneous which hampers genetic counseling and molecular diagnosis. Testing of several single HI-related genes is laborious and expensive. In this study, we evaluate the diagnostic utility of whole-exome sequencing (WES) targeting a panel of HI-related genes. Two hundred index patients, mostly of Dutch origin, with presumed hereditary HI underwent WES followed by targeted analysis of an HI gene panel of 120 genes. We found causative variants underlying the HI in 67 of 200 patients (33.5%). Eight of these patients have a large homozygous deletion involving STRC, OTOA or USH2A, which could only be identified by copy number variation detection. Variants of uncertain significance were found in 10 patients (5.0%). In the remaining 123 cases, no potentially causative variants were detected (61.5%). In our patient cohort, causative variants in GJB2, USH2A, MYO15A and STRC, and in MYO6 were the leading causes for autosomal recessive and dominant HI, respectively. Segregation analysis and functional analyses of variants of uncertain significance will probably further increase the diagnostic yield of WES.

Truncating de novo mutations in the Krüppel-type zinc-finger gene ZNF148 in patients with corpus callosum defects, developmental delay, short stature, and dysmorphisms.

BACKGROUND: Krüppel-type zinc finger genes (ZNF) constitute a large yet relatively poorly characterized gene family. ZNF genes encode proteins that recognize specific DNA motifs in gene promotors. They act as transcriptional co-activators or -repressors via interaction with chromatin remodeling proteins and other transcription factors. Only few ZNF genes are currently linked to human disorders and identification of ZNF gene-associated human diseases may help understand their function. Here we provide genetic, statistical, and clinical evidence to support association of ZNF148 with a new intellectual disability (ID) syndrome disorder. METHODS: Routine diagnostic exome sequencing data were obtained from 2172 patients with ID and/or multiple congenital anomalies. RESULTS: In a cohort of 2172 patient-parent trios referred for routine diagnostic whole exome sequencing for ID and/or multiple congenital anomalies (MCA) in the period 2012-2016, four patients were identified who carried de novo heterozygous nonsense or frameshift mutations in the ZNF148 gene. This was the only ZNF gene with recurrent truncating de novo mutations in this cohort. All mutations resulted in premature termination codons in the last exon of ZNF148. The number of the de novo truncating mutations in the ZNF148 gene was significantly enriched (p = 5.42 × 10-3). The newly described ZNF148-associated syndrome is characterized by underdevelopment of the corpus callosum, mild to moderate developmental delay and ID, variable microcephaly or mild macrocephaly, short stature, feeding problems, facial dysmorphisms, and cardiac and renal malformations. CONCLUSIONS: We propose ZNF148 as a gene involved in a newly described ID syndrome with a recurrent phenotype and postulate that the ZNF148 is a hitherto unrecognized but crucial transcription factor in the development of the corpus callosum. Our study illustrates the advantage of whole exome sequencing in a large cohort using a parent-offspring trio approach for identifying novel genes involved in rare human diseases.

Meta-analysis of 2,104 trios provides support for 10 new genes for intellectual disability.

To identify candidate genes for intellectual disability, we performed a meta-analysis on 2,637 de novo mutations, identified from the exomes of 2,104 patient-parent trios. Statistical analyses identified 10 new candidate ID genes: DLG4, PPM1D, RAC1, SMAD6, SON, SOX5, SYNCRIP, TCF20, TLK2 and TRIP12. In addition, we show that these genes are intolerant to nonsynonymous variation and that mutations in these genes are associated with specific clinical ID phenotypes.

Novel bioinformatic developments for exome sequencing.

With the widespread adoption of next generation sequencing technologies by the genetics community and the rapid decrease in costs per base, exome sequencing has become a standard within the repertoire of genetic experiments for both research and diagnostics. Although bioinformatics now offers standard solutions for the analysis of exome sequencing data, many challenges still remain; especially the increasing scale at which exome data are now being generated has given rise to novel challenges in how to efficiently store, analyze and interpret exome data of this magnitude. In this review we discuss some of the recent developments in bioinformatics for exome sequencing and the directions that this is taking us to. With these developments, exome sequencing is paving the way for the next big challenge, the application of whole genome sequencing.

ConBind: motif-aware cross-species alignment for the identification of functional transcription factor binding sites.

Eukaryotic gene expression is regulated by transcription factors (TFs) binding to promoter as well as distal enhancers. TFs recognize short, but specific binding sites (TFBSs) that are located within the promoter and enhancer regions. Functionally relevant TFBSs are often highly conserved during evolution leaving a strong phylogenetic signal. While multiple sequence alignment (MSA) is a potent tool to detect the phylogenetic signal, the current MSA implementations are optimized to align the maximum number of identical nucleotides. This approach might result in the omission of conserved motifs that contain interchangeable nucleotides such as the ETS motif (IUPAC code: GGAW). Here, we introduce ConBind, a novel method to enhance alignment of short motifs, even if their mutual sequence similarity is only partial. ConBind improves the identification of conserved TFBSs by improving the alignment accuracy of TFBS families within orthologous DNA sequences. Functional validation of the Gfi1b + 13 enhancer reveals that ConBind identifies additional functionally important ETS binding sites that were missed by all other tested alignment tools. In addition to the analysis of known regulatory regions, our web tool is useful for the analysis of TFBSs on so far unknown DNA regions identified through ChIP-sequencing.