Rapid exome sequencing as a first-tier test in neonates with suspected genetic disorder: results of a prospective multicenter clinical utility study in the Netherlands.

The introduction of rapid exome sequencing (rES) for critically ill neonates admitted to the neonatal intensive care unit has made it possible to impact clinical decision-making. Unbiased prospective studies to quantify the impact of rES over routine genetic testing are, however, scarce. We performed a clinical utility study to compare rES to conventional genetic diagnostic workup for critically ill neonates with suspected genetic disorders. In a multicenter prospective parallel cohort study involving five Dutch NICUs, we performed rES in parallel to routine genetic testing for 60 neonates with a suspected genetic disorder and monitored diagnostic yield and the time to diagnosis. To assess the economic impact of rES, healthcare resource use was collected for all neonates. rES detected more conclusive genetic diagnoses than routine genetic testing (20% vs. 10%, respectively), in a significantly shorter time to diagnosis (15 days (95% CI 10-20) vs. 59 days (95% CI 23-98, p < 0.001)). Moreover, rES reduced genetic diagnostic costs by 1.5% (€85 per neonate). CONCLUSION:  Our findings demonstrate the clinical utility of rES for critically ill neonates based on increased diagnostic yield, shorter time to diagnosis, and net healthcare savings. Our observations warrant the widespread implementation of rES as first-tier genetic test in critically ill neonates with disorders of suspected genetic origin. WHAT IS KNOWN: • Rapid exome sequencing (rES) enables diagnosing rare genetic disorders in a fast and reliable manner, but retrospective studies with neonates admitted to the neonatal intensive care unit (NICU) indicated that genetic disorders are likely underdiagnosed as rES is not routinely used. • Scenario modeling for implementation of rES for neonates with presumed genetic disorders indicated an expected increase in costs associated with genetic testing. WHAT IS NEW: • This unique prospective national clinical utility study of rES in a NICU setting shows that rES obtained more and faster diagnoses than conventional genetic tests. • Implementation of rES as replacement for all other genetic tests does not increase healthcare costs but in fact leads to a reduction in healthcare costs.

Phenotypic expansion of EGP5-related Vici syndrome: 15 Dutch patients carrying a founder variant.

Vici syndrome (OMIM 242840) is a very rare autosomal recessive multisystem disorder first described in 1988. In 2013, bi-allelic loss-of-function mutations in EPG5 were reported to cause Vici syndrome. Five principal diagnostic features of Vici syndrome have been proposed: agenesis of the corpus callosum, cataracts, cardiomyopathy, hypopigmentation, and combined immunodeficiency. We identified 15 patients carrying a homozygous founder missense variant in EPG5 who all exhibit a less severe clinical phenotype than classic Vici syndrome. All 15 show typical brain abnormalities on MRI. The homozygous founder variant in EPG5 they carry results in a shorter in-frame transcript and truncated, but likely still residual, EPG5 protein. We speculate that the residual EPG5 protein explains their attenuated phenotype, which is consistent with two previous observations that low expression of EPG5 can lead to an attenuated Vici syndrome phenotype. We propose renaming this condition EPG5-related neurodevelopmental disorder to emphasize the clinical variability of patients with bi-allelic mutations in EPG5.

Towards Next-Generation Sequencing (NGS)-Based Newborn Screening: A Technical Study to Prepare for the Challenges Ahead.

Newborn screening (NBS) aims to identify neonates with severe conditions for whom immediate treatment is required. Currently, a biochemistry-first approach is used to identify these disorders, which are predominantly inherited meta1bolic disorders (IMD). Next-generation sequencing (NGS) is expected to have some advantages over the current approach, for example the ability to detect IMDs that meet all screening criteria but lack an identifiable biochemical footprint. We have now designed a technical study to explore the use of NGS techniques as a first-tier approach in NBS. Here, we describe the aim and set-up of the NGS-first for the NBS (NGSf4NBS) project, which will proceed in three steps. In Step 1, we will identify IMDs eligible for NGS-first testing, based on treatability. In Step 2, we will investigate the feasibility, limitations and comparability of different technical NGS approaches and analysis workflows for NBS, eventually aiming to develop a rapid NGS-based workflow. Finally, in Step 3, we will prepare for the incorporation of this workflow into the existing Dutch NBS program and propose a protocol for referral of a child after a positive NGS test result. The results of this study will be the basis for an additional analytical route within NBS that will be further studied for its applicability within the NBS program, e.g., regarding the ethical, legal, financial and social implications.

Feasibility of predicting allele specific expression from DNA sequencing using machine learning.

Allele specific expression (ASE) concerns divergent expression quantity of alternative alleles and is measured by RNA sequencing. Multiple studies show that ASE plays a role in hereditary diseases by modulating penetrance or phenotype severity. However, genome diagnostics is based on DNA sequencing and therefore neglects gene expression regulation such as ASE. To take advantage of ASE in absence of RNA sequencing, it must be predicted using only DNA variation. We have constructed ASE models from BIOS (n = 3432) and GTEx (n = 369) that predict ASE using DNA features. These models are highly reproducible and comprise many different feature types, highlighting the complex regulation that underlies ASE. We applied the BIOS-trained model to population variants in three genes in which ASE plays a clinically relevant role: BRCA2, RET and NF1. This resulted in predicted ASE effects for 27 variants, of which 10 were known pathogenic variants. We demonstrated that ASE can be predicted from DNA features using machine learning. Future efforts may improve sensitivity and translate these models into a new type of genome diagnostic tool that prioritizes candidate pathogenic variants or regulators thereof for follow-up validation by RNA sequencing. All used code and machine learning models are available at GitHub and Zenodo.

Diagnostic yield of targeted next generation sequencing in 2002 Dutch cardiomyopathy patients.

BACKGROUND: Next-generation sequencing (NGS) is increasingly used for clinical evaluation of cardiomyopathy patients as it allows for simultaneous screening of multiple cardiomyopathy-associated genes. Adding copy number variant (CNV) analysis of NGS data is not routine yet and may contribute to the diagnostic yield. OBJECTIVES: Determine the diagnostic yield of our targeted NGS gene panel in routine clinical diagnostics of Dutch cardiomyopathy patients and explore the impact of exon CNVs on diagnostic yield. METHODS: Patients (N = 2002) referred for clinical genetic analysis underwent diagnostic testing of 55-61 genes associated with cardiomyopathies. Samples were analyzed and evaluated for single nucleotide variants (SNVs), indels and CNVs. CNVs identified in the NGS data and suspected of being pathogenic based on type, size and location were confirmed by additional molecular tests. RESULTS: A (likely) pathogenic (L)P variant was detected in 22.7% of patients, including 3 with CNVs and 25 where a variant was identified in a gene currently not associated with the patient's cardiomyopathy subtype. Only 15 out of 2002 patients (0.8%) were found to carry two (L)P variants. CONCLUSION: The yield of routine clinical diagnostics of cardiomyopathies was relatively low when compared to literature. This is likely due to the fact that our study reports the outcome of patients in daily routine diagnostics, therefore also including patients not fully fulfilling (subtype specific) cardiomyopathy criteria. This may also explain why (L)P variants were identified in genes not associated with the reported subtype. The added value of CNV analysis was shown to be limited but not negligible.

Rare functional missense variants in CACNA1H: What can we learn from Writer's cramp?

Writer's cramp (WC) is a task-specific focal dystonia that occurs selectively in the hand and arm during writing. Previous studies have shown a role for genetics in the pathology of task-specific focal dystonia. However, to date, no causal gene has been reported for task-specific focal dystonia, including WC. In this study, we investigated the genetic background of a large Dutch family with autosomal dominant‒inherited WC that was negative for mutations in known dystonia genes. Whole exome sequencing identified 4 rare variants of unknown significance that segregated in the family. One candidate gene was selected for follow-up, Calcium Voltage-Gated Channel Subunit Alpha1 H, CACNA1H, due to its links with the known dystonia gene Potassium Channel Tetramerization Domain Containing 17, KCTD17, and with paroxysmal movement disorders. Targeted resequencing of CACNA1H in 82 WC cases identified another rare, putative damaging variant in a familial WC case that did not segregate. Using structural modelling and functional studies in vitro, we show that both the segregating p.Arg481Cys variant and the non-segregating p.Glu1881Lys variant very likely cause structural changes to the Cav3.2 protein and lead to similar gains of function, as seen in an accelerated recovery from inactivation. Both mutant channels are thus available for re-activation earlier, which may lead to an increase in intracellular calcium and increased neuronal excitability. Overall, we conclude that rare functional variants in CACNA1H need to be interpreted very carefully, and additional studies are needed to prove that the p.Arg481Cys variant is the cause of WC in the large Dutch family.

CAPICE: a computational method for Consequence-Agnostic Pathogenicity Interpretation of Clinical Exome variations.

Exome sequencing is now mainstream in clinical practice. However, identification of pathogenic Mendelian variants remains time-consuming, in part, because the limited accuracy of current computational prediction methods requires manual classification by experts. Here we introduce CAPICE, a new machine-learning-based method for prioritizing pathogenic variants, including SNVs and short InDels. CAPICE outperforms the best general (CADD, GAVIN) and consequence-type-specific (REVEL, ClinPred) computational prediction methods, for both rare and ultra-rare variants. CAPICE is easily added to diagnostic pipelines as pre-computed score file or command-line software, or using online MOLGENIS web service with API. Download CAPICE for free and open-source (LGPLv3) at https://github.com/molgenis/capice .

Homozygous damaging SOD2 variant causes lethal neonatal dilated cardiomyopathy.

BACKGROUND: Idiopathic dilated cardiomyopathy (DCM) is recognised to be a heritable disorder, yet clinical genetic testing does not produce a diagnosis in >50% of paediatric patients. Identifying a genetic cause is crucial because this knowledge can affect management options, cardiac surveillance in relatives and reproductive decision-making. In this study, we sought to identify the underlying genetic defect in a patient born to consanguineous parents with rapidly progressive DCM that led to death in early infancy. METHODS AND RESULTS: Exome sequencing revealed a potentially pathogenic, homozygous missense variant, c.542G>T, p.(Gly181Val), in SOD2. This gene encodes superoxide dismutase 2 (SOD2) or manganese-superoxide dismutase, a mitochondrial matrix protein that scavenges oxygen radicals produced by oxidation-reduction and electron transport reactions occurring in mitochondria via conversion of superoxide anion (O2-·) into H2O2. Measurement of hydroethidine oxidation showed a significant increase in O2-· levels in the patient's skin fibroblasts, as compared with controls, and this was paralleled by reduced catalytic activity of SOD2 in patient fibroblasts and muscle. Lentiviral complementation experiments demonstrated that mitochondrial SOD2 activity could be completely restored on transduction with wild type SOD2. CONCLUSION: Our results provide evidence that defective SOD2 may lead to toxic increases in the levels of damaging oxygen radicals in the neonatal heart, which can result in rapidly developing heart failure and death. We propose SOD2 as a novel nuclear-encoded mitochondrial protein involved in severe human neonatal cardiomyopathy, thus expanding the wide range of genetic factors involved in paediatric cardiomyopathies.

Natural Exon Skipping Sets the Stage for Exon Skipping as Therapy for Dystrophic Epidermolysis Bullosa.

Dystrophic epidermolysis bullosa (DEB) is a devastating blistering disease affecting skin and mucous membranes. It is caused by pathogenic variants in the COL7A1 gene encoding type VII collagen, and can be inherited dominantly or recessively. Recently, promising proof-of-principle has been shown for antisense oligonucleotide (AON)-mediated exon skipping as a therapeutic approach for DEB. However, the precise phenotypic effect to be anticipated from exon skipping, and which patient groups could benefit, is not yet clear. To answer these questions, we studied new clinical and molecular data on seven patients from the Dutch EB registry and reviewed the literature on COL7A1 exon skipping variants. We found that phenotypes associated with dominant exon skipping cannot be distinguished from phenotypes caused by other dominant DEB variants. Recessive exon skipping phenotypes are generally relatively mild in the spectrum of recessive DEB. Therefore, for dominant DEB, AON-mediated exon skipping is unlikely to ameliorate the phenotype. In contrast, the overall severity of phenotypes associated with recessive natural exon skipping pivots toward the milder end of the spectrum. Consequently, we anticipate AON-mediated exon skipping for recessive DEB caused by bi-allelic null variants should lead to a clinically relevant improvement of this devastating phenotype.

Dutch genome diagnostic laboratories accelerated and improved variant interpretation and increased accuracy by sharing data.

Each year diagnostic laboratories in the Netherlands profile thousands of individuals for heritable disease using next-generation sequencing (NGS). This requires pathogenicity classification of millions of DNA variants on the standard 5-tier scale. To reduce time spent on data interpretation and increase data quality and reliability, the nine Dutch labs decided to publicly share their classifications. Variant classifications of nearly 100,000 unique variants were catalogued and compared in a centralized MOLGENIS database. Variants classified by more than one center were labeled as "consensus" when classifications agreed, and shared internationally with LOVD and ClinVar. When classifications opposed (LB/B vs. LP/P), they were labeled "conflicting", while other nonconsensus observations were labeled "no consensus". We assessed our classifications using the InterVar software to compare to ACMG 2015 guidelines, showing 99.7% overall consistency with only 0.3% discrepancies. Differences in classifications between Dutch labs or between Dutch labs and ACMG were mainly present in genes with low penetrance or for late onset disorders and highlight limitations of the current 5-tier classification system. The data sharing boosted the quality of DNA diagnostics in Dutch labs, an initiative we hope will be followed internationally. Recently, a positive match with a case from outside our consortium resulted in a more definite disease diagnosis.

No major role for rare plectin variants in arrhythmogenic right ventricular cardiomyopathy.

AIMS: Likely pathogenic/pathogenic variants in genes encoding desmosomal proteins play an important role in the pathophysiology of arrhythmogenic right ventricular cardiomyopathy (ARVC). However, for a substantial proportion of ARVC patients, the genetic substrate remains unknown. We hypothesized that plectin, a cytolinker protein encoded by the PLEC gene, could play a role in ARVC because it has been proposed to link the desmosomal protein desmoplakin to the cytoskeleton and therefore has a potential function in the desmosomal structure. METHODS: We screened PLEC in 359 ARVC patients and compared the frequency of rare coding PLEC variants (minor allele frequency [MAF] <0.001) between patients and controls. To assess the frequency of rare variants in the control population, we evaluated the rare coding variants (MAF <0.001) found in the European cohort of the Exome Aggregation Database. We further evaluated plectin localization by immunofluorescence in a subset of patients with and without a PLEC variant. RESULTS: Forty ARVC patients carried one or more rare PLEC variants (11%, 40/359). However, rare variants also seem to occur frequently in the control population (18%, 4754/26197 individuals). Nor did we find a difference in the prevalence of rare PLEC variants in ARVC patients with or without a desmosomal likely pathogenic/pathogenic variant (14% versus 8%, respectively). However, immunofluorescence analysis did show decreased plectin junctional localization in myocardial tissue from 5 ARVC patients with PLEC variants. CONCLUSIONS: Although PLEC has been hypothesized as a promising candidate gene for ARVC, our current study did not show an enrichment of rare PLEC variants in ARVC patients compared to controls and therefore does not support a major role for PLEC in this disorder. Although rare PLEC variants were associated with abnormal localization in cardiac tissue, the confluence of data does not support a role for plectin abnormalities in ARVC development.

Naturally Occurring Variants in LRP1 (Low-Density Lipoprotein Receptor-Related Protein 1) Affect HDL (High-Density Lipoprotein) Metabolism Through ABCA1 (ATP-Binding Cassette A1) and SR-B1 (Scavenger Receptor Class B Type 1) in Humans.

OBJECTIVE: Studies into the role of LRP1 (low-density lipoprotein receptor-related protein 1) in human lipid metabolism are scarce. Although it is known that a common variant in LRP1 (rs116133520) is significantly associated with HDL-C (high-density lipoprotein cholesterol), the mechanism underlying this observation is unclear. In this study, we set out to study the functional effects of 2 rare LRP1 variants identified in subjects with extremely low HDL-C levels. APPROACH AND RESULTS: In 2 subjects with HDL-C below the first percentile for age and sex and moderately elevated triglycerides, we identified 2 rare variants in LRP1: p.Val3244Ile and p.Glu3983Asp. Both variants decrease LRP1 expression and stability. We show in a series of translational experiments that these variants culminate in reduced trafficking of ABCA1 (ATP-binding cassette A1) to the cell membrane. This is accompanied by an increase in cell surface expression of SR-B1 (scavenger receptor class B type 1). Combined these effects may contribute to low HDL-C levels in our study subjects. Supporting these findings, we provide epidemiological evidence that rs116133520 is associated with apo (apolipoprotein) A1 but not with apoB levels. CONCLUSIONS: This study provides the first evidence that rare variants in LRP1 are associated with changes in human lipid metabolism. Specifically, this study shows that LRP1 may affect HDL metabolism by virtue of its effect on both ABCA1 and SR-B1.

Genetic Screening Test to Detect Translocations in Acute Leukemias by Use of Targeted Locus Amplification.

BACKGROUND: Over 500 translocations have been identified in acute leukemia. To detect them, most diagnostic laboratories use karyotyping, fluorescent in situ hybridization, and reverse transcription PCR. Targeted locus amplification (TLA), a technique using next-generation sequencing, now allows detection of the translocation partner of a specific gene, regardless of its chromosomal origin. We present a TLA multiplex assay as a potential first-tier screening test for detecting translocations in leukemia diagnostics. METHODS: The panel includes 17 genes involved in many translocations present in acute leukemias. Procedures were optimized by using a training set of cell line dilutions and 17 leukemia patient bone marrow samples and validated by using a test set of cell line dilutions and a further 19 patient bone marrow samples. Per gene, we determined if its region was involved in a translocation and, if so, the translocation partner. To balance sensitivity and specificity, we introduced a gray zone showing indeterminate translocation calls needing confirmation. We benchmarked our method against results from the 3 standard diagnostic tests. RESULTS: In patient samples passing QC, we achieved a concordance with benchmarking tests of 81% in the training set and 100% in the test set, after confirmation of 4 and nullification of 3 gray zone calls (in total). In cell line dilutions, we detected translocations in 10% aberrant cells at several genetic loci. CONCLUSIONS: Multiplex TLA shows promising results as an acute leukemia screening test. It can detect cryptic and other translocations in selected genes. Further optimization may make this assay suitable for diagnostic use.

Toward an effective exome-based genetic testing strategy in pediatric dilated cardiomyopathy.

PURPOSE: We evaluated the diagnostic yield in pediatric dilated cardiomyopathy (DCM) of combining exome sequencing (ES)-based targeted analysis and genome-wide copy-number variation (CNV) analysis. Based on our findings, we retrospectively designed an effective approach for genetic testing in pediatric DCM. METHODS: We identified 95 patients (in 85 families) with pediatric onset of DCM. We initially excluded 13 of these families because they already had a genetic diagnosis, leaving a total of 31 probands for single-nucleotide polymorphism (SNP) array and trio-ES. We used Human Phenotype Ontology (HPO)-based filtering for our data analysis. RESULTS: We reached a genetic diagnosis in 15/31 (48.4%) families. ES yielded a diagnosis in 13 probands (13/15; 86.7%), with most variants being found in genes encoding structural cardiomyocyte components. Two large deletions were identified using SNP array. If we had included the 13 excluded families, our estimated yield would have been 54%. CONCLUSION: We propose a standardized, stepwise analysis of (i) well-known cardiomyopathy genes, (ii) CNVs, (iii) all genes assigned to HPO cardiomyopathy, and (iv) if appropriate, genes assigned to other HPO terms. This diagnostic approach yields the highest increase at each subsequent step and reduces analytic effort, cost, the number of variants of unknown clinical significance, and the chance of incidental findings.

Exome sequencing and network analysis identifies shared mechanisms underlying spinocerebellar ataxia.

The autosomal dominant cerebellar ataxias, referred to as spinocerebellar ataxias in genetic nomenclature, are a rare group of progressive neurodegenerative disorders characterized by loss of balance and coordination. Despite the identification of numerous disease genes, a substantial number of cases still remain without a genetic diagnosis. Here, we report five novel spinocerebellar ataxia genes, FAT2, PLD3, KIF26B, EP300, and FAT1, identified through a combination of exome sequencing in genetically undiagnosed families and targeted resequencing of exome candidates in a cohort of singletons. We validated almost all genes genetically, assessed damaging effects of the gene variants in cell models and further consolidated a role for several of these genes in the aetiology of spinocerebellar ataxia through network analysis. Our work links spinocerebellar ataxia to alterations in synaptic transmission and transcription regulation, and identifies these as the main shared mechanisms underlying the genetically diverse spinocerebellar ataxia types.

Rapid Targeted Genomics in Critically Ill Newborns.

BACKGROUND: Rapid diagnostic whole-genome sequencing has been explored in critically ill newborns, hoping to improve their clinical care and replace time-consuming and/or invasive diagnostic testing. A previous retrospective study in a research setting showed promising results with diagnoses in 57%, but patients were highly selected for known and likely Mendelian disorders. The aim of our prospective study was to assess the speed and yield of rapid targeted genomic diagnostics for clinical application. METHODS: We included 23 critically ill children younger than 12 months in ICUs over a period of 2 years. A quick diagnosis could not be made after routine clinical evaluation and diagnostics. Targeted analysis of 3426 known disease genes was performed by using whole-genome sequencing data. We measured diagnostic yield, turnaround times, and clinical consequences. RESULTS: A genetic diagnosis was obtained in 7 patients (30%), with a median turnaround time of 12 days (ranging from 5 to 23 days). We identified compound heterozygous mutations in the EPG5 gene (Vici syndrome), the RMND1 gene (combined oxidative phosphorylation deficiency-11), and the EIF2B5 gene (vanishing white matter), and homozygous mutations in the KLHL41 gene (nemaline myopathy), the GFER gene (progressive mitochondrial myopathy), and the GLB1 gene (GM1-gangliosidosis). In addition, a 1p36.33p36.32 microdeletion was detected in a child with cardiomyopathy. CONCLUSIONS: Rapid targeted genomics combined with copy number variant detection adds important value in the neonatal and pediatric intensive care setting. It led to a fast diagnosis in 30% of critically ill children for whom the routine clinical workup was unsuccessful.

X-linked hypomyelination with spondylometaphyseal dysplasia (H-SMD) associated with mutations in AIFM1.

An X-linked condition characterized by the combination of hypomyelinating leukodystrophy and spondylometaphyseal dysplasia (H-SMD) has been observed in only four families, with linkage to Xq25-27, and recent genetic characterization in two families with a common AIFM1 mutation. In our study, 12 patients (6 families) with H-SMD were identified and underwent comprehensive assessment accompanied by whole-exome sequencing (WES). Pedigree analysis in all families was consistent with X-linked recessive inheritance. Presentation typically occurred between 12 and 36 months. In addition to the two disease-defining features of spondylometaphyseal dysplasia and hypomyelination on MRI, common clinical signs and symptoms included motor deterioration, spasticity, tremor, ataxia, dysarthria, cognitive defects, pulmonary hypertension, nystagmus, and vision loss due to retinopathy. The course of the disease was slowly progressive. All patients had maternally inherited or de novo mutations in or near exon 7 of AIFM1, within a region of 70 bp, including synonymous and intronic changes. AIFM1 mutations have previously been associated with neurologic presentations as varied as intellectual disability, hearing loss, neuropathy, and striatal necrosis, while AIFM1 mutations in this small region present with a distinct phenotype implicating bone. Analysis of cell lines derived from four patients identified significant reductions in AIFM1 mRNA and protein levels in osteoblasts. We hypothesize that AIFM1 functions in bone metabolism and myelination and is responsible for the unique phenotype in this condition.

Using the shared genetics of dystonia and ataxia to unravel their pathogenesis.

In this review we explore the similarities between spinocerebellar ataxias and dystonias, and suggest potentially shared molecular pathways using a gene co-expression network approach. The spinocerebellar ataxias are a group of neurodegenerative disorders characterized by coordination problems caused mainly by atrophy of the cerebellum. The dystonias are another group of neurological movement disorders linked to basal ganglia dysfunction, although evidence is now pointing to cerebellar involvement as well. Our gene co-expression network approach identified 99 shared genes and showed the involvement of two major pathways: synaptic transmission and neurodevelopment. These pathways overlapped in the two disorders, with a large role for GABAergic signaling in both. The overlapping pathways may provide novel targets for disease therapies. We need to prioritize variants obtained by whole exome sequencing in the genes associated with these pathways in the search for new pathogenic variants, which can than be used to help in the genetic counseling of patients and their families.