Prevalence of comorbidities in individuals with neurodevelopmental disorders from the aggregated phenomics data of 51,227 pediatric individuals.

The prevalence of comorbidities in individuals with neurodevelopmental disorders (NDDs) is not well understood, yet these are important for accurate diagnosis and prognosis in routine care and for characterizing the clinical spectrum of NDD syndromes. We thus developed PhenomAD-NDD, an aggregated database containing the comorbid phenotypic data of 51,227 individuals with NDD, all harmonized into Human Phenotype Ontology (HPO), with in total 3,054 unique HPO terms. We demonstrate that almost all congenital anomalies are more prevalent in the NDD population than in the general population, and the NDD baseline prevalence allows for an approximation of the enrichment of symptoms. For example, such analyses of 33 genetic NDDs show that 32% of enriched phenotypes are currently not reported in the clinical synopsis in the Online Mendelian Inheritance in Man (OMIM). PhenomAD-NDD is open to all via a visualization online tool and allows us to determine the enrichment of symptoms in NDD.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

BACKGROUND: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney, caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative mode of action, wherein an increased level of AFF3 resulted in pathological effects. METHODS: Evolutionary constraints suggest that other modes-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be damaging variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. RESULTS: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous Loss-of-Function (LoF) or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not rescue these phenotypes. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring + / + , KINSSHIP/KINSSHIP, LoF/ + , LoF/LoF or KINSSHIP/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the KINSSHIP/KINSSHIP or the LoF/LoF lines. While the same pathways are affected, only about one third of the differentially expressed genes are common to the homozygote datasets, indicating that AFF3 LoF and KINSSHIP variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. CONCLUSIONS: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

The detection of a strong episignature for Chung-Jansen syndrome, partially overlapping with Börjeson-Forssman-Lehmann and White-Kernohan syndromes.

Chung-Jansen syndrome is a neurodevelopmental disorder characterized by intellectual disability, behavioral problems, obesity and dysmorphic features. It is caused by pathogenic variants in the PHIP gene that encodes for the Pleckstrin homology domain-interacting protein, which is part of an epigenetic modifier protein complex. Therefore, we hypothesized that PHIP haploinsufficiency may impact genome-wide DNA methylation (DNAm). We assessed the DNAm profiles of affected individuals with pathogenic and likely pathogenic PHIP variants with Infinium Methylation EPIC arrays and report a specific and sensitive DNAm episignature biomarker for Chung-Jansen syndrome. In addition, we observed similarities between the methylation profile of Chung-Jansen syndrome and that of functionally related and clinically partially overlapping genetic disorders, White-Kernohan syndrome (caused by variants in DDB1 gene) and Börjeson-Forssman-Lehmann syndrome (caused by variants in PHF6 gene). Based on these observations we also proceeded to develop a common episignature biomarker for these disorders. These newly defined episignatures can be used as part of a multiclass episignature classifier for screening of affected individuals with rare disorders and interpretation of genetic variants of unknown clinical significance, and provide further insights into the common molecular pathophysiology of the clinically-related Chung-Jansen, Börjeson-Forssman-Lehmann and White-Kernohan syndromes.

Refining the 9q34.3 microduplication syndrome reveals mild neurodevelopmental features associated with a distinct global DNA methylation profile.

Precise regulation of gene expression is important for correct neurodevelopment. 9q34.3 deletions affecting the EHMT1 gene result in a syndromic neurodevelopmental disorder named Kleefstra syndrome. In contrast, duplications of the 9q34.3 locus encompassing EHMT1 have been suggested to cause developmental disorders, but only limited information has been available. We have identified 15 individuals from 10 unrelated families, with 9q34.3 duplications <1.5 Mb in size, encompassing EHMT1 entirely. Clinical features included mild developmental delay, mild intellectual disability or learning problems, autism spectrum disorder, and behavior problems. The individuals did not consistently display dysmorphic features, congenital anomalies, or growth abnormalities. DNA methylation analysis revealed a weak DNAm profile for the cases with 9q34.3 duplication encompassing EHMT1, which could segregate the majority of the affected cases from controls. This study shows that individuals with 9q34.3 duplications including EHMT1 gene present with mild non-syndromic neurodevelopmental disorders and DNA methylation changes different from Kleefstra syndrome.

Genome sequencing as a generic diagnostic strategy for rare disease.

BACKGROUND: To diagnose the full spectrum of hereditary and congenital diseases, genetic laboratories use many different workflows, ranging from karyotyping to exome sequencing. A single generic high-throughput workflow would greatly increase efficiency. We assessed whether genome sequencing (GS) can replace these existing workflows aimed at germline genetic diagnosis for rare disease. METHODS: We performed short-read GS (NovaSeq™6000; 150 bp paired-end reads, 37 × mean coverage) on 1000 cases with 1271 known clinically relevant variants, identified across different workflows, representative of our tertiary diagnostic centers. Variants were categorized into small variants (single nucleotide variants and indels < 50 bp), large variants (copy number variants and short tandem repeats) and other variants (structural variants and aneuploidies). Variant calling format files were queried per variant, from which workflow-specific true positive rates (TPRs) for detection were determined. A TPR of ≥ 98% was considered the threshold for transition to GS. A GS-first scenario was generated for our laboratory, using diagnostic efficacy and predicted false negative as primary outcome measures. As input, we modeled the diagnostic path for all 24,570 individuals referred in 2022, combining the clinical referral, the transition of the underlying workflow(s) to GS, and the variant type(s) to be detected. RESULTS: Overall, 95% (1206/1271) of variants were detected. Detection rates differed per variant category: small variants in 96% (826/860), large variants in 93% (341/366), and other variants in 87% (39/45). TPRs varied between workflows (79-100%), with 7/10 being replaceable by GS. Models for our laboratory indicate that a GS-first strategy would be feasible for 84.9% of clinical referrals (750/883), translating to 71% of all individuals (17,444/24,570) receiving GS as their primary test. An estimated false negative rate of 0.3% could be expected. CONCLUSIONS: GS can capture clinically relevant germline variants in a 'GS-first strategy' for the majority of clinical indications in a genetics diagnostic lab.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

Background: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney,caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative (DN) mode-of-action, wherein an increased level of AFF3 resulted in pathological effects. Methods: Evolutionary constraints suggest that other mode-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be deleterious variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. Results: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous LoF or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not complement. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring +/+, DN/DN, LoF/+, LoF/LoF or DN/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the DN/DN or the LoF/LoF lines. While the same pathways are affected, only about one-third of the differentially expressed genes are common to these homozygote datasets, indicating that AFF3 LoF and DN variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. Conclusions: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Reanalysis of whole-exome sequencing (WES) data of children with neurodevelopmental disorders in a standard patient care context.

This study aims to inform future genetic reanalysis management by evaluating the yield of whole-exome sequencing (WES) reanalysis in standard patient care in the Netherlands. Single-center data of 159 patients with a neurodevelopmental disorder (NDD), in which WES analysis and reanalysis were performed between January 1, 2014, and December 31, 2021, was retrospectively collected. Patients were included if they were under the age of 18 years at initial analysis and if this initial analysis did not result in a diagnosis. Demographic, phenotypic, and genotypic characteristics of patients were collected and analyzed. The primary outcomes of our study were (i) diagnostic yield at reanalysis, (ii) reasons for detecting a new possibly causal variant at reanalysis, (iii) unsolicited findings, and (iv) factors associated with positive result of reanalysis. In addition, we conducted a questionnaire study amongst the 7 genetic department in the Netherlands creating an overview of used techniques, yield, and organization of WES reanalysis. The single-center data show that in most cases, WES reanalysis was initiated by the clinical geneticist (65%) or treating physician (30%). The mean time between initial WES analysis and reanalysis was 3.7 years. A new (likely) pathogenic variant or VUS with a clear link to the phenotype was found in 20 initially negative cases, resulting in a diagnostic yield of 12.6%. In 75% of these patients, the diagnosis had clinical consequences, as for example, a screening plan for associated signs and symptoms could be devised. Most (32%) of the (likely) causal variants identified at WES reanalysis were discovered due to a newly described gene-disease association. In addition to the 12.6% diagnostic yield based on new diagnoses, reclassification of a variant of uncertain significance found at initial analysis led to a definite diagnosis in three patients. Diagnostic yield was higher in patients with dysmorphic features compared to patients without clear dysmorphic features (yield 27% vs. 6%; p = 0.001). CONCLUSIONS: Our results show that WES reanalysis in patients with NDD in standard patient care leads to a substantial increase in genetic diagnoses. In the majority of newly diagnosed patients, the diagnosis had clinical consequences. Knowledge about the clinical impact of WES reanalysis, clinical characteristics associated with higher yield, and the yield per year after a negative WES in larger clinical cohorts is warranted to inform guidelines for genetic reanalysis. These guidelines will be of great value for pediatricians, pediatric rehabilitation specialists, and pediatric neurologists in daily care of patients with NDD. WHAT IS KNOWN: • Whole exome sequencing can cost-effectively identify a genetic cause of intellectual disability in about 30-40% of patients. • WES reanalysis in a research setting can lead to a definitive diagnosis in 10-20% of previously exome negative cases. WHAT IS NEW: • WES reanalysis in standard patient care resulted in a diagnostic yield of 13% in previously exome negative children with NDD. • The presence of dysmorphic features is associated with an increased diagnostic yield of WES reanalysis.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

DNA methylation episignature, extension of the clinical features, and comparative epigenomic profiling of Hao-Fountain syndrome caused by variants in USP7.

PURPOSE: Hao-Fountain syndrome (HAFOUS) is a neurodevelopmental disorder caused by pathogenic variants in USP7. HAFOUS is characterized by developmental delay, intellectual disability, speech delay, behavioral abnormalities, autism spectrum disorder, seizures, hypogonadism, and mild dysmorphic features. We investigated the phenotype of 18 participants with HAFOUS and performed DNA methylation (DNAm) analysis, aiming to generate a diagnostic biomarker. Furthermore, we performed comparative analysis with known episignatures to gain more insight into the molecular pathophysiology of HAFOUS. METHODS: We assessed genomic DNAm profiles of 18 individuals with pathogenic variants and variants of uncertain significance (VUS) in USP7 to map and validate a specific episignature. The comparison between the USP7 cohort and 56 rare genetic disorders with earlier reported DNAm episignatures was performed with statistical and functional correlation. RESULTS: We mapped a sensitive and specific DNAm episignature for pathogenic variants in USP7 and utilized this to reclassify the VUS. Comparative epigenomic analysis showed evidence of HAFOUS similarity to a number of other rare genetic episignature disorders. CONCLUSION: We discovered a sensitive and specific DNAm episignature as a robust diagnostic biomarker for HAFOUS that enables VUS reclassification in USP7. We also expand the phenotypic spectrum of 9 new and 5 previously reported individuals with HAFOUS.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework.

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

Comprehensive de novo mutation discovery with HiFi long-read sequencing.

BACKGROUND: Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease. METHODS: We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing. RESULTS: We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data. CONCLUSIONS: HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

A complex structural variant near SOX3 causes X-linked split-hand/foot malformation.

Split-hand/foot malformation (SHFM) is a congenital limb defect most typically presenting with median clefts in hands and/or feet, that can occur in a syndromic context as well as in isolated form. SHFM is caused by failure to maintain normal apical ectodermal ridge function during limb development. Although several genes and contiguous gene syndromes are implicated in the monogenic etiology of isolated SHFM, the disorder remains genetically unexplained for many families and associated genetic loci. We describe a family with isolated X-linked SHFM, for which the causative variant could be detected after a diagnostic journey of 20 years. We combined well-established approaches including microarray-based copy number variant analysis and fluorescence in situ hybridization coupled with optical genome mapping and whole genome sequencing. This strategy identified a complex structural variant (SV) comprising a 165-kb gain of 15q26.3 material ([GRCh37/hg19] chr15:99795320-99960362dup) inserted in inverted position at the site of a 38-kb deletion on Xq27.1 ([GRCh37/hg19] chrX:139481061-139518989del). In silico analysis suggested that the SV disrupts the regulatory framework on the X chromosome and may lead to SOX3 misexpression. We hypothesize that SOX3 dysregulation in the developing limb disturbed the fine balance between morphogens required for maintaining AER function, resulting in SHFM in this family.

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.

The Genetics of Intellectual Disability.

Intellectual disability (ID) has a prevalence of ~2-3% in the general population, having a large societal impact. The underlying cause of ID is largely of genetic origin; however, identifying this genetic cause has in the past often led to long diagnostic Odysseys. Over the past decades, improvements in genetic diagnostic technologies and strategies have led to these causes being more and more detectable: from cytogenetic analysis in 1959, we moved in the first decade of the 21st century from genomic microarrays with a diagnostic yield of ~20% to next-generation sequencing platforms with a yield of up to 60%. In this review, we discuss these various developments, as well as their associated challenges and implications for the field of ID, which highlight the revolutionizing shift in clinical practice from a phenotype-first into genotype-first approach.

Whole genome sequencing for USH2A-associated disease reveals several pathogenic deep-intronic variants that are amenable to splice correction.

A significant number of individuals with a rare disorder such as Usher syndrome (USH) and (non-)syndromic autosomal recessive retinitis pigmentosa (arRP) remain genetically unexplained. Therefore, we assessed subjects suspected of USH2A-associated disease and no or mono-allelic USH2A variants using whole genome sequencing (WGS) followed by an improved pipeline for variant interpretation to provide a conclusive diagnosis. One hundred subjects were screened using WGS to identify causative variants in USH2A or other USH/arRP-associated genes. In addition to the existing variant interpretation pipeline, a particular focus was put on assessing splice-affecting properties of variants, both in silico and in vitro. Also structural variants were extensively addressed. For variants resulting in pseudoexon inclusion, we designed and evaluated antisense oligonucleotides (AONs) using minigene splice assays and patient-derived photoreceptor precursor cells. Biallelic variants were identified in 49 of 100 subjects, including novel splice-affecting variants and structural variants, in USH2A or arRP/USH-associated genes. Thirteen variants were shown to affect USH2A pre-mRNA splicing, including four deep-intronic USH2A variants resulting in pseudoexon inclusion, which could be corrected upon AON treatment. We have shown that WGS, combined with a thorough variant interpretation pipeline focused on assessing pre-mRNA splicing defects and structural variants, is a powerful method to provide subjects with a rare genetic condition, a (likely) conclusive genetic diagnosis. This is essential for the development of future personalized treatments and for patients to be eligible for such treatments.

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing.

PURPOSE: Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the "ClinVar low-hanging fruit" reanalysis, reasons for the failure of previous analyses, and lessons learned. METHODS: Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. RESULTS: We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). CONCLUSION: The "ClinVar low-hanging fruit" analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock.

De novo mutation hotspots in homologous protein domains identify function-altering mutations in neurodevelopmental disorders.

Variant interpretation remains a major challenge in medical genetics. We developed Meta-Domain HotSpot (MDHS) to identify mutational hotspots across homologous protein domains. We applied MDHS to a dataset of 45,221 de novo mutations (DNMs) from 31,058 individuals with neurodevelopmental disorders (NDDs) and identified three significantly enriched missense DNM hotspots in the ion transport protein domain family (PF00520). The 37 unique missense DNMs that drive enrichment affect 25 genes, 19 of which were previously associated with NDDs. 3D protein structure modeling supports the hypothesis of function-altering effects of these mutations. Hotspot genes have a unique expression pattern in tissue, and we used this pattern alongside in silico predictors and population constraint information to identify candidate NDD-associated genes. We also propose a lenient version of our method, which identifies 32 hotspot positions across 16 different protein domains. These positions are enriched for likely pathogenic variation in clinical databases and DNMs in other genetic disorders.