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.

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.

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.

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.

Panel-based NGS Reveals Novel Pathogenic Mutations in Autosomal Recessive Retinitis Pigmentosa.

Retinitis pigmentosa (RP) is a group of inherited progressive retinal dystrophies (RD) characterized by photoreceptor degeneration. RP is highly heterogeneous both clinically and genetically, which complicates the identification of causative genes and mutations. Targeted next-generation sequencing (NGS) has been demonstrated to be an effective strategy for the detection of mutations in RP. In our study, an in-house gene panel comprising 75 known RP genes was used to analyze a cohort of 47 unrelated Spanish families pre-classified as autosomal recessive or isolated RP. Disease-causing mutations were found in 27 out of 47 cases achieving a mutation detection rate of 57.4%. In total, 33 pathogenic mutations were identified, 20 of which were novel mutations (60.6%). Furthermore, not only single nucleotide variations but also copy-number variations, including three large deletions in the USH2A and EYS genes, were identified. Finally seven out of 27 families, displaying mutations in the ABCA4, RP1, RP2 and USH2A genes, could be genetically or clinically reclassified. These results demonstrate the potential of our panel-based NGS strategy in RP diagnosis.

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.

Allelic Mutations of KITLG, Encoding KIT Ligand, Cause Asymmetric and Unilateral Hearing Loss and Waardenburg Syndrome Type 2.

Linkage analysis combined with whole-exome sequencing in a large family with congenital and stable non-syndromic unilateral and asymmetric hearing loss (NS-UHL/AHL) revealed a heterozygous truncating mutation, c.286_303delinsT (p.Ser96Ter), in KITLG. This mutation co-segregated with NS-UHL/AHL as a dominant trait with reduced penetrance. By screening a panel of probands with NS-UHL/AHL, we found an additional mutation, c.200_202del (p.His67_Cys68delinsArg). In vitro studies revealed that the p.His67_Cys68delinsArg transmembrane isoform of KITLG is not detectable at the cell membrane, supporting pathogenicity. KITLG encodes a ligand for the KIT receptor. Also, KITLG-KIT signaling and MITF are suggested to mutually interact in melanocyte development. Because mutations in MITF are causative of Waardenburg syndrome type 2 (WS2), we screened KITLG in suspected WS2-affected probands. A heterozygous missense mutation, c.310C>G (p.Leu104Val), that segregated with WS2 was identified in a small family. In vitro studies revealed that the p.Leu104Val transmembrane isoform of KITLG is located at the cell membrane, as is wild-type KITLG. However, in culture media of transfected cells, the p.Leu104Val soluble isoform of KITLG was reduced, and no soluble p.His67_Cys68delinsArg and p.Ser96Ter KITLG could be detected. These data suggest that mutations in KITLG associated with NS-UHL/AHL have a loss-of-function effect. We speculate that the mechanism of the mutation underlying WS2 and leading to membrane incorporation and reduced secretion of KITLG occurs via a dominant-negative or gain-of-function effect. Our study unveils different phenotypes associated with KITLG, previously associated with pigmentation abnormalities, and will thereby improve the genetic counseling given to individuals with KITLG variants.

Comparison of Exome and Genome Sequencing Technologies for the Complete Capture of Protein-Coding Regions.

For next-generation sequencing technologies, sufficient base-pair coverage is the foremost requirement for the reliable detection of genomic variants. We investigated whether whole-genome sequencing (WGS) platforms offer improved coverage of coding regions compared with whole-exome sequencing (WES) platforms, and compared single-base coverage for a large set of exome and genome samples. We find that WES platforms have improved considerably in the last years, but at comparable sequencing depth, WGS outperforms WES in terms of covered coding regions. At higher sequencing depth (95x-160x), WES successfully captures 95% of the coding regions with a minimal coverage of 20x, compared with 98% for WGS at 87-fold coverage. Three different assessments of sequence coverage bias showed consistent biases for WES but not for WGS. We found no clear differences for the technologies concerning their ability to achieve complete coverage of 2,759 clinically relevant genes. We show that WES performs comparable to WGS in terms of covered bases if sequenced at two to three times higher coverage. This does, however, go at the cost of substantially more sequencing biases in WES approaches. Our findings will guide laboratories to make an informed decision on which sequencing platform and coverage to choose.