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.

Genomic Diagnosis of Rare Pediatric Disease in the United Kingdom and Ireland.

BACKGROUND: Pediatric disorders include a range of highly penetrant, genetically heterogeneous conditions amenable to genomewide diagnostic approaches. Finding a molecular diagnosis is challenging but can have profound lifelong benefits. METHODS: We conducted a large-scale sequencing study involving more than 13,500 families with probands with severe, probably monogenic, difficult-to-diagnose developmental disorders from 24 regional genetics services in the United Kingdom and Ireland. Standardized phenotypic data were collected, and exome sequencing and microarray analyses were performed to investigate novel genetic causes. We developed an iterative variant analysis pipeline and reported candidate variants to clinical teams for validation and diagnostic interpretation to inform communication with families. Multiple regression analyses were performed to evaluate factors affecting the probability of diagnosis. RESULTS: A total of 13,449 probands were included in the analyses. On average, we reported 1.0 candidate variant per parent-offspring trio and 2.5 variants per singleton proband. Using clinical and computational approaches to variant classification, we made a diagnosis in approximately 41% of probands (5502 of 13,449). Of 3599 probands in trios who received a diagnosis by clinical assertion, approximately 76% had a pathogenic de novo variant. Another 22% of probands (2997 of 13,449) had variants of uncertain significance in genes that were strongly linked to monogenic developmental disorders. Recruitment in a parent-offspring trio had the largest effect on the probability of diagnosis (odds ratio, 4.70; 95% confidence interval [CI], 4.16 to 5.31). Probands were less likely to receive a diagnosis if they were born extremely prematurely (i.e., 22 to 27 weeks' gestation; odds ratio, 0.39; 95% CI, 0.22 to 0.68), had in utero exposure to antiepileptic medications (odds ratio, 0.44; 95% CI, 0.29 to 0.67), had mothers with diabetes (odds ratio, 0.52; 95% CI, 0.41 to 0.67), or were of African ancestry (odds ratio, 0.51; 95% CI, 0.31 to 0.78). CONCLUSIONS: Among probands with severe, probably monogenic, difficult-to-diagnose developmental disorders, multimodal analysis of genomewide data had good diagnostic power, even after previous attempts at diagnosis. (Funded by the Health Innovation Challenge Fund and Wellcome Sanger Institute.).

Differences in the number of de novo mutations between individuals are due to small family-specific effects and stochasticity.

The number of de novo mutations (DNMs) in the human germline is correlated with parental age at conception, but this explains only part of the observed variation. We investigated whether there is a family-specific contribution to the number of DNMs in offspring. The analysis of DNMs in 111 dizygotic twin pairs did not identify a substantial family-specific contribution. This result was corroborated by comparing DNMs of 1669 siblings to those of age-matched unrelated offspring following correction for parental age. In addition, by modeling DNM data from 1714 multi-offspring families, we estimated that the family-specific contribution explains ∼5.2% of the variation in DNM number. Furthermore, we found no substantial difference between the observed number of DNMs and those predicted by a stochastic Poisson process. We conclude that there is a small family-specific contribution to DNM number and that stochasticity explains a large proportion of variation in DNM counts.

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing.

BACKGROUND: Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques. RESULTS: We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection. CONCLUSION: We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques.

Missense variants in ANKRD11 cause KBG syndrome by impairment of stability or transcriptional activity of the encoded protein.

PURPOSE: Although haploinsufficiency of ANKRD11 is among the most common genetic causes of neurodevelopmental disorders, the role of rare ANKRD11 missense variation remains unclear. We characterized clinical, molecular, and functional spectra of ANKRD11 missense variants. METHODS: We collected clinical information of individuals with ANKRD11 missense variants and evaluated phenotypic fit to KBG syndrome. We assessed pathogenicity of variants through in silico analyses and cell-based experiments. RESULTS: We identified 20 unique, mostly de novo, ANKRD11 missense variants in 29 individuals, presenting with syndromic neurodevelopmental disorders similar to KBG syndrome caused by ANKRD11 protein truncating variants or 16q24.3 microdeletions. Missense variants significantly clustered in repression domain 2 at the ANKRD11 C-terminus. Of the 10 functionally studied missense variants, 6 reduced ANKRD11 stability. One variant caused decreased proteasome degradation and loss of ANKRD11 transcriptional activity. CONCLUSION: Our study indicates that pathogenic heterozygous ANKRD11 missense variants cause the clinically recognizable KBG syndrome. Disrupted transrepression capacity and reduced protein stability each independently lead to ANKRD11 loss-of-function, consistent with haploinsufficiency. This highlights the diagnostic relevance of ANKRD11 missense variants, but also poses diagnostic challenges because the KBG-associated phenotype may be mild and inherited pathogenic ANKRD11 (missense) variants are increasingly observed, warranting stringent variant classification and careful phenotyping.

Application of metabolite set enrichment analysis on untargeted metabolomics data prioritises relevant pathways and detects novel biomarkers for inherited metabolic disorders.

Untargeted metabolomics (UM) allows for the simultaneous measurement of hundreds of metabolites in a single analytical run. The sheer amount of data generated in UM hampers its use in patient diagnostics because manual interpretation of all features is not feasible. Here, we describe the application of a pathway-based metabolite set enrichment analysis method to prioritise relevant biological pathways in UM data. We validate our method on a set of 55 patients with a diagnosed inherited metabolic disorder (IMD) and show that it complements feature-based prioritisation of biomarkers by placing the features in a biological context. In addition, we find that by taking enriched pathways shared across different IMDs, we can identify common drugs and compounds that could otherwise obscure genuine disease biomarkers in an enrichment method. Finally, we demonstrate the potential of this method to identify novel candidate biomarkers for known IMDs. Our results show the added value of pathway-based interpretation of UM data in IMD diagnostics context.

Genomic instability in non-breast/ovarian malignancies of individuals with germline pathogenic variants in BRCA1/2.

BACKGROUND: Individuals with germline pathogenic variants (gPVs) in BRCA1 or BRCA2 (BRCA1/2) are at a high risk of breast- and ovarian carcinomas (BOCs) with BRCA1/2-deficiency and homologous recombination deficiency (HRD) that can be detected by analysis of genome-wide genomic instability features such as large-scale state transitions, telomeric allelic imbalances and genomic loss-of-heterozygosity. Malignancies with HRD are more sensitive to platinum-based therapies and PARP inhibitors. Here, we aim to investigate the fraction of non-BOC malignancies that have BRCA1/2-deficiency and genomic instability features. METHODS: The full tumor history of a large historical clinic-based consecutive cohort of 2,965 individuals with gPVs in BRCA1/2 was retrieved via the Dutch nationwide pathology databank (Palga). In total, 169 non-BOC malignancies were collected and analyzed with targeted next-generation sequencing and shallow whole-genome sequencing to determine somatic second hit alterations and genomic instabilities indicative of HRD, respectively. RESULTS: BRCA1/2-deficiency was detected in 27% (21/79) and 23% (21/90) of 20 different types of non-BOC malignancies of individuals with gPVs in BRCA1 and BRCA2, respectively. These malignancies had a higher genomic instability score than BRCA1- or BRCA2-proficient malignancies (P < .001 and P < .001, respectively). CONCLUSIONS: BRCA1/2-deficiency and genomic instability features were found in 27% and 23% of a broad spectrum of non-BOC malignancies in individuals with gPVs in BRCA1 and BRCA2, respectively. Evaluation of the effectivity of PARP-inhibitors in these individuals should be focused on tumors with confirmed absence of a wild type allele.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation.

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.