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.

Homozygous mutation of STXBP5L explains an autosomal recessive infantile-onset neurodegenerative disorder.

We report siblings of consanguineous parents with an infantile-onset neurodegenerative disorder manifesting a predominant sensorimotor axonal neuropathy, optic atrophy and cognitive deficit. We used homozygosity mapping to identify an ∼12-Mbp interval identical by descent (IBD) between the affected individuals on chromosome 3q13.13-21.1 with an LOD score of 2.31. We combined family-based whole-exome and whole-genome sequencing of parents and affected siblings and, after filtering of likely non-pathogenic variants, identified a unique missense variant in syntaxin-binding protein 5-like (STXBP5L c.3127G>A, p.Val1043Ile [CCDS43137.1]) in the IBD interval. Considering other modes of inheritance, we also found compound heterozygous variants in FMNL3 (c.114G>C, p.Phe38Leu and c.1372T>G, p.Ile458Leu [CCDS44874.1]) located on chromosome 12. STXBP5L (or Tomosyn-2) is expressed in the central and peripheral nervous system and is known to inhibit neurotransmitter release through inhibition of the formation of the SNARE complexes between synaptic vesicles and the plasma membrane. FMNL3 is expressed more widely and is a formin family protein that is involved in the regulation of cell morphology and cytoskeletal organization. The STXBP5L p.Val1043Ile variant enhanced inhibition of exocytosis in comparison with wild-type (WT) STXBP5L. Furthermore, WT STXBP5L, but not variant STXBP5L, promoted axonal outgrowth in manipulated mouse primary hippocampal neurons. However, the FMNL3 p.Phe38Leu and p.Ile458Leu variants showed minimal effects in these cells. Collectively, our clinical, genetic and molecular data suggest that the IBD variant in STXBP5L is the likely cause of the disorder.

Detection of clinically relevant copy-number variants by exome sequencing in a large cohort of genetic disorders.

PURPOSE: Copy-number variation is a common source of genomic variation and an important genetic cause of disease. Microarray-based analysis of copy-number variants (CNVs) has become a first-tier diagnostic test for patients with neurodevelopmental disorders, with a diagnostic yield of 10-20%. However, for most other genetic disorders, the role of CNVs is less clear and most diagnostic genetic studies are generally limited to the study of single-nucleotide variants (SNVs) and other small variants. With the introduction of exome and genome sequencing, it is now possible to detect both SNVs and CNVs using an exome- or genome-wide approach with a single test. METHODS: We performed exome-based read-depth CNV screening on data from 2,603 patients affected by a range of genetic disorders for which exome sequencing was performed in a diagnostic setting. RESULTS: In total, 123 clinically relevant CNVs ranging in size from 727 bp to 15.3 Mb were detected, which resulted in 51 conclusive diagnoses and an overall increase in diagnostic yield of ~2% (ranging from 0 to -5.8% per disorder). CONCLUSIONS: This study shows that CNVs play an important role in a broad range of genetic disorders and that detection via exome-based CNV profiling results in an increase in the diagnostic yield without additional testing, bringing us closer to single-test genomics.Genet Med advance online publication 27 October 2016.

Cantú syndrome is caused by mutations in ABCC9.

Cantú syndrome is a rare disorder characterized by congenital hypertrichosis, neonatal macrosomia, a distinct osteochondrodysplasia, and cardiomegaly. Using an exome-sequencing approach applied to one proband-parent trio and three unrelated single cases, we identified heterozygous mutations in ABCC9 in all probands. With the inclusion of the remaining cohort of ten individuals with Cantú syndrome, a total of eleven mutations in ABCC9 were found. The de novo occurrence in all six simplex cases in our cohort substantiates the presence of a dominant disease mechanism. All mutations were missense, and several mutations affect Arg1154. This mutation hot spot lies within the second type 1 transmembrane region of this ATP-binding cassette transporter protein, which may suggest an activating mutation. ABCC9 encodes the sulfonylurea receptor (SUR) that forms ATP-sensitive potassium channels (K(ATP) channels) originally shown in cardiac, skeletal, and smooth muscle. Previously, loss-of-function mutations in this gene have been associated with idiopathic dilated cardiomyopathy type 10 (CMD10). These findings identify the genetic basis of Cantú syndrome and suggest that this is a new member of the potassium channelopathies.

A post-hoc comparison of the utility of sanger sequencing and exome sequencing for the diagnosis of heterogeneous diseases.

The advent of massive parallel sequencing is rapidly changing the strategies employed for the genetic diagnosis and research of rare diseases that involve a large number of genes. So far it is not clear whether these approaches perform significantly better than conventional single gene testing as requested by clinicians. The current yield of this traditional diagnostic approach depends on a complex of factors that include gene-specific phenotype traits, and the relative frequency of the involvement of specific genes. To gauge the impact of the paradigm shift that is occurring in molecular diagnostics, we assessed traditional Sanger-based sequencing (in 2011) and exome sequencing followed by targeted bioinformatics analysis (in 2012) for five different conditions that are highly heterogeneous, and for which our center provides molecular diagnosis. We find that exome sequencing has a much higher diagnostic yield than Sanger sequencing for deafness, blindness, mitochondrial disease, and movement disorders. For microsatellite-stable colorectal cancer, this was low under both strategies. Even if all genes that could have been ordered by physicians had been tested, the larger number of genes captured by the exome would still have led to a clearly superior diagnostic yield at a fraction of the cost.

Next-generation sequencing of a 40 Mb linkage interval reveals TSPAN12 mutations in patients with familial exudative vitreoretinopathy.

Familial exudative vitreoretinopathy (FEVR) is a genetically heterogeneous retinal disorder characterized by abnormal vascularisation of the peripheral retina, often accompanied by retinal detachment. To date, mutations in three genes (FZD4, LRP5, and NDP) have been shown to be causative for FEVR. In two large Dutch pedigrees segregating autosomal-dominant FEVR, genome-wide SNP analysis identified an FEVR locus of approximately 40 Mb on chromosome 7. Microsatellite marker analysis suggested similar at risk haplotypes in patients of both families. To identify the causative gene, we applied next-generation sequencing in the proband of one of the families, by analyzing all exons and intron-exon boundaries of 338 genes, in addition to microRNAs, noncoding RNAs, and other highly conserved genomic regions in the 40 Mb linkage interval. After detailed bioinformatic analysis of the sequence data, prioritization of all detected sequence variants led to three candidates to be considered as the causative genetic defect in this family. One of these variants was an alanine-to-proline substitution in the transmembrane 4 superfamily member 12 protein, encoded by TSPAN12. This protein has very recently been implicated in regulating the development of retinal vasculature, together with the proteins encoded by FZD4, LRP5, and NDP. Sequence analysis of TSPAN12 revealed two mutations segregating in five of 11 FEVR families, indicating that mutations in TSPAN12 are a relatively frequent cause of FEVR. Furthermore, we demonstrate the power of targeted next-generation sequencing technology to identify disease genes in linkage intervals.

Targeted next-generation sequencing of a 12.5 Mb homozygous region reveals ANO10 mutations in patients with autosomal-recessive cerebellar ataxia.

Autosomal-recessive cerebellar ataxias comprise a clinically and genetically heterogeneous group of neurodegenerative disorders. In contrast to their dominant counterparts, unraveling the molecular background of these ataxias has proven to be more complicated and the currently known mutations provide incomplete coverage for genotyping of patients. By combining SNP array-based linkage analysis and targeted resequencing of relevant sequences in the linkage interval with the use of next-generation sequencing technology, we identified a mutation in a gene and have shown its association with autosomal-recessive cerebellar ataxia. In a Dutch consanguineous family with three affected siblings a homozygous 12.5 Mb region on chromosome 3 was targeted by array-based sequence capture. Prioritization of all detected sequence variants led to four candidate genes, one of which contained a variant with a high base pair conservation score (phyloP score: 5.26). This variant was a leucine-to-arginine substitution in the DUF 590 domain of a 16K transmembrane protein, a putative calcium-activated chloride channel encoded by anoctamin 10 (ANO10). The analysis of ANO10 by Sanger sequencing revealed three additional mutations: a homozygous mutation (c.1150_1151del [p.Leu384fs]) in a Serbian family and a compound-heterozygous splice-site mutation (c.1476+1G>T) and a frameshift mutation (c.1604del [p.Leu535X]) in a French family. This illustrates the power of using initial homozygosity mapping with next-generation sequencing technology to identify genes involved in autosomal-recessive diseases. Moreover, identifying a putative calcium-dependent chloride channel involved in cerebellar ataxia adds another pathway to the list of pathophysiological mechanisms that may cause cerebellar ataxia.

Next-generation genetic testing for retinitis pigmentosa.

Molecular diagnostics for patients with retinitis pigmentosa (RP) has been hampered by extreme genetic and clinical heterogeneity, with 52 causative genes known to date. Here, we developed a comprehensive next-generation sequencing (NGS) approach for the clinical molecular diagnostics of RP. All known inherited retinal disease genes (n = 111) were captured and simultaneously analyzed using NGS in 100 RP patients without a molecular diagnosis. A systematic data analysis pipeline was developed and validated to prioritize and predict the pathogenicity of all genetic variants identified in each patient, which enabled us to reduce the number of potential pathogenic variants from approximately 1,200 to zero to nine per patient. Subsequent segregation analysis and in silico predictions of pathogenicity resulted in a molecular diagnosis in 36 RP patients, comprising 27 recessive, six dominant, and three X-linked cases. Intriguingly, De novo mutations were present in at least three out of 28 isolated cases with causative mutations. This study demonstrates the enormous potential and clinical utility of NGS in molecular diagnosis of genetically heterogeneous diseases such as RP. De novo dominant mutations appear to play a significant role in patients with isolated RP, having major implications for genetic counselling.

Reanalysis of exome negative patients with rare disease: a pragmatic workflow for diagnostic applications.

BACKGROUND: Approximately two third of patients with a rare genetic disease remain undiagnosed after exome sequencing (ES). As part of our post-test counseling procedures, patients without a conclusive diagnosis are advised to recontact their referring clinician to discuss new diagnostic opportunities in due time. We performed a systematic study of genetically undiagnosed patients 5 years after their initial negative ES report to determine the efficiency of diverse reanalysis strategies. METHODS: We revisited a cohort of 150 pediatric neurology patients originally enrolled at Radboud University Medical Center, of whom 103 initially remained genetically undiagnosed. We monitored uptake of physician-initiated routine clinical and/or genetic re-evaluation (ad hoc re-evaluation) and performed systematic reanalysis, including ES-based resequencing, of all genetically undiagnosed patients (systematic re-evaluation). RESULTS: Ad hoc re-evaluation was initiated for 45 of 103 patients and yielded 18 diagnoses (including 1 non-genetic). Subsequent systematic re-evaluation identified another 14 diagnoses, increasing the diagnostic yield in our cohort from 31% (47/150) to 53% (79/150). New genetic diagnoses were established by reclassification of previously identified variants (10%, 3/31), reanalysis with enhanced bioinformatic pipelines (19%, 6/31), improved coverage after resequencing (29%, 9/31), and new disease-gene associations (42%, 13/31). Crucially, our systematic study also showed that 11 of the 14 further conclusive genetic diagnoses were made in patients without a genetic diagnosis that did not recontact their referring clinician. CONCLUSIONS: We find that upon re-evaluation of undiagnosed patients, both reanalysis of existing ES data as well as resequencing strategies are needed to identify additional genetic diagnoses. Importantly, not all patients are routinely re-evaluated in clinical care, prolonging their diagnostic trajectory, unless systematic reanalysis is facilitated. We have translated our observations into considerations for systematic and ad hoc reanalysis in routine genetic care.

Accurate distinction of pathogenic from benign CNVs in mental retardation.

Copy number variants (CNVs) have recently been recognized as a common form of genomic variation in humans. Hundreds of CNVs can be detected in any individual genome using genomic microarrays or whole genome sequencing technology, but their phenotypic consequences are still poorly understood. Rare CNVs have been reported as a frequent cause of neurological disorders such as mental retardation (MR), schizophrenia and autism, prompting widespread implementation of CNV screening in diagnostics. In previous studies we have shown that, in contrast to benign CNVs, MR-associated CNVs are significantly enriched in genes whose mouse orthologues, when disrupted, result in a nervous system phenotype. In this study we developed and validated a novel computational method for differentiating between benign and MR-associated CNVs using structural and functional genomic features to annotate each CNV. In total 13 genomic features were included in the final version of a Naïve Bayesian Tree classifier, with LINE density and mouse knock-out phenotypes contributing most to the classifier's accuracy. After demonstrating that our method (called GECCO) perfectly classifies CNVs causing known MR-associated syndromes, we show that it achieves high accuracy (94%) and negative predictive value (99%) on a blinded test set of more than 1,200 CNVs from a large cohort of individuals with MR. These results indicate that this classification method will be of value for objectively prioritizing CNVs in clinical research and diagnostics.

Massively parallel sequencing of ataxia genes after array-based enrichment.

Massively parallel sequencing has tremendous diagnostic potential but requires enriched templates for sequencing. Here we report the validation of an array-based sequence capture method in genetically heterogeneous disorders. The model disorder selected was AR ataxia, using five subjects with known mutations and two unaffected controls. The genomic sequences of seven disease genes, together with two control loci were targeted on a 2-Mb sequence-capture array. After enrichment, the patients' DNA samples were analyzed using one-quarter Roche GS FLX Titanium sequencing run, resulting in an average of 65 Mb of sequence data per patient. This was sufficient for an average 25-fold coverage/base in all targeted regions. Enrichment showed high specificity; on average, 80% of uniquely mapped reads were on target. Importantly, this approach enabled automated detection of deletions and hetero- and homozygous point mutations for 6/7 mutant alleles, and greater than 99% accuracy for known SNP variants. Our results also clearly show reduced coverage for sequences in repeat-rich regions, which significantly impacts the reliable detection of genomic variants. Based on these findings we recommend a minimal coverage of 15-fold for diagnostic implementation of this technology. We conclude that massive parallel sequencing of enriched samples enables personalized diagnosis of heterogeneous genetic disorders and qualifies for rapid diagnostic implementation.

Targeted next generation sequencing reveals a novel intragenic deletion of the TPO gene in a family with intellectual disability.

BACKGROUNDS AND AIMS: Next generation sequencing (NGS) approaches have revolutionized the identification of mutations underlying genetic disorders. This technology is particularly useful for the identification of mutations in known and new genes for conditions with extensive genetic heterogeneity. In the present study we investigated a consanguineous Pakistani family with intellectual disability (ID). METHODS: Genotyping was carried out using 250k and 6k SNP microarrays in order to perform homozygosity mapping and copy number variation (CNV) analysis. Targeted NGS was performed to identify the genetic defect in this family. qPCR was performed to validate and confirm the NGS result. RESULTS: Homozygosity mapping positioned the causative defect on chromosome 2p25.3-p25.2. Subsequent targeted NGS revealed an intragenic deletion of five exons of the gene TPO. CONCLUSIONS: NGS is a powerful method to uncover submicroscopic structural variations. This result demonstrates that an unbiased screening approach such as NGS can help to identify even unexpected disease-causing mutations.

Molecular diagnosis of Usher syndrome: application of two different next generation sequencing-based procedures.

Usher syndrome (USH) is a clinically and genetically heterogeneous disorder characterized by visual and hearing impairments. Clinically, it is subdivided into three subclasses with nine genes identified so far. In the present study, we investigated whether the currently available Next Generation Sequencing (NGS) technologies are already suitable for molecular diagnostics of USH. We analyzed a total of 12 patients, most of which were negative for previously described mutations in known USH genes upon primer extension-based microarray genotyping. We enriched the NGS template either by whole exome capture or by Long-PCR of the known USH genes. The main NGS sequencing platforms were used: SOLiD for whole exome sequencing, Illumina (Genome Analyzer II) and Roche 454 (GS FLX) for the Long-PCR sequencing. Long-PCR targeting was more efficient with up to 94% of USH gene regions displaying an overall coverage higher than 25×, whereas whole exome sequencing yielded a similar coverage for only 50% of those regions. Overall this integrated analysis led to the identification of 11 novel sequence variations in USH genes (2 homozygous and 9 heterozygous) out of 18 detected. However, at least two cases were not genetically solved. Our result highlights the current limitations in the diagnostic use of NGS for USH patients. The limit for whole exome sequencing is linked to the need of a strong coverage and to the correct interpretation of sequence variations with a non obvious, pathogenic role, whereas the targeted approach suffers from the high genetic heterogeneity of USH that may be also caused by the presence of additional causative genes yet to be identified.

A de novo paradigm for mental retardation.

The per-generation mutation rate in humans is high. De novo mutations may compensate for allele loss due to severely reduced fecundity in common neurodevelopmental and psychiatric diseases, explaining a major paradox in evolutionary genetic theory. Here we used a family based exome sequencing approach to test this de novo mutation hypothesis in ten individuals with unexplained mental retardation. We identified and validated unique non-synonymous de novo mutations in nine genes. Six of these, identified in six different individuals, are likely to be pathogenic based on gene function, evolutionary conservation and mutation impact. Our findings provide strong experimental support for a de novo paradigm for mental retardation. Together with de novo copy number variation, de novo point mutations of large effect could explain the majority of all mental retardation cases in the population.

De novo mutations of SETBP1 cause Schinzel-Giedion syndrome.

Schinzel-Giedion syndrome is characterized by severe mental retardation, distinctive facial features and multiple congenital malformations; most affected individuals die before the age of ten. We sequenced the exomes of four affected individuals (cases) and found heterozygous de novo variants in SETBP1 in all four. We also identified SETBP1 mutations in eight additional cases using Sanger sequencing. All mutations clustered to a highly conserved 11-bp exonic region, suggesting a dominant-negative or gain-of-function effect.