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.

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.

Variation in a range of mTOR-related genes associates with intracranial volume and intellectual disability

De novo mutations in specific mTOR pathway genes cause brain overgrowth in the context of intellectual disability (ID). By analyzing 101 mMTOR-related genes in a large ID patient cohort and two independent population cohorts, we show that these genes modulate brain growth in health and disease. We report the mTOR activator gene RHEB as an ID gene that is associated with megalencephaly when mutated. Functional testing of mutant RHEB in vertebrate animal models indicates pathway hyperactivation with a concomitant increase in cell and head size, aberrant neuronal migration, and induction of seizures, concordant with the human phenotype. This study reveals that tight control of brain volume is exerted through a large community of mTOR-related genes. Human brain volume can be altered, by either rare disruptive events causing hyperactivation of the pathway, or through the collective effects of common alleles.

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.

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.

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.

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.

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.