Pathogenic paralogous variants can be used to apply the ACMG PS1 and PM5 variant interpretation criteria 2023.08.22.23294353

E-print/Working paper (2023)

Purpose The majority of missense variants in clinical genetic tests are classified as variants of uncertain significance. Broadening the evidence of the PS1 and PM5 criteria has the potential to increase ... [more ▼]

Purpose The majority of missense variants in clinical genetic tests are classified as variants of uncertain significance. Broadening the evidence of the PS1 and PM5 criteria has the potential to increase conclusive variant interpretation. Methods We hypothesized that incorporation of pathogenic missense variants in conserved residues across paralogous genes can increase the number of variants where ACMG PS1/PM5 criteria can be applied. We mapped over 2.5 million pathogenic and general population variants from ClinVar, HGMD, and gnomAD databases onto 9,990 genes and aligned these by gene families. Subsequently, we developed a novel framework to extend PS1/PM5 by incorporating pathogenic paralogous variants annotations (para-PS1/PM5). Results We demonstrate that para-PS1/PM5 criteria increase the number of classifiable amino acids 3.6-fold compared to PS1 and PM5. Across all gene families with at least two disease-associated genes, the calculated likelihood ratios suggest moderate evidence for pathogenicity. Moreover, for 36 genes, the extended para-PS1/PM5 criteria reach strong evidence level. Conclusion We show that single pathogenic paralogous variants incorporation at paralogous protein positions increases the applicability of the PS1 and PM5 criteria, likely leading to a reduction of variants of uncertain significance across many monogenic disorders. Future iterations of the ACMG guidelines may consider para-PS1 and para-PM5. [less ▲]

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

in Brain: a Journal of Neurology (2023), 146(3), 923-934

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without ... [more ▼]

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future. [less ▲]

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes 2022.03.23.485339

E-print/Working paper (2022)

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without ... [more ▼]

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our novel 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5 A distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and fucntional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants, and 679 electrophysiological experiments that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.Competing Interest StatementThe authors have declared no competing interest.DSSPDictionary of Protein Secondary StructuregnomADGenome aggregation DatabaseGoFGain of functionGRIN genesGRIN1, GRIN2A. GRIN2BHGMDHuman Gene Mutation DatabaseNMDA receptorN-methyl-D-aspartate receptorGABA receptorGamma-aminobutyric acid receptorLoFLoss of functionSCN genesSCN1A, SCN2A, SCN8AVCFVariant Call Format [less ▲]

Distinct gene-set burden patterns underlie common generalized and focal epilepsies

in EBioMedicine (2021), 72

Background Analyses of few gene-sets in epilepsy showed a potential to unravel key disease associations. We set out to investigate the burden of ultra-rare variants (URVs) in a comprehensive range of ... [more ▼]

Background Analyses of few gene-sets in epilepsy showed a potential to unravel key disease associations. We set out to investigate the burden of ultra-rare variants (URVs) in a comprehensive range of biologically informed gene-sets presumed to be implicated in epileptogenesis. Methods The burden of 12 URV types in 92 gene-sets was compared between cases and controls using whole exome sequencing data from individuals of European descent with developmental and epileptic encephalopathies (DEE, n = 1,003), genetic generalized epilepsy (GGE, n = 3,064), or non-acquired focal epilepsy (NAFE, n = 3,522), collected by the Epi25 Collaborative, compared to 3,962 ancestry-matched controls. Findings Missense URVs in highly constrained regions were enriched in neuron-specific and developmental genes, whereas genes not expressed in brain were not affected. GGE featured a higher burden in gene-sets derived from inhibitory vs. excitatory neurons or associated receptors, whereas the opposite was found for NAFE, and DEE featured a burden in both. Top-ranked susceptibility genes from recent genome-wide association studies (GWAS) and gene-sets derived from generalized vs. focal epilepsies revealed specific enrichment patterns of URVs in GGE vs. NAFE. Interpretation Missense URVs affecting highly constrained sites differentially impact genes expressed in inhibitory vs. excitatory pathways in generalized vs. focal epilepsies. The excess of URVs in top-ranked GWAS risk-genes suggests a convergence of rare deleterious and common risk-variants in the pathogenesis of generalized and focal epilepsies. [less ▲]

Guideline-based and bioinformatic reassessment of lesion-associated gene and variant pathogenicity in focal human epilepsies

in Epilepsia (2018)

Objective: Increasing availability of surgically resected brain tissue from patients with focal epilepsy and Focal Cortical Dysplasia (FCD) or low-grade glio-neuronal tumors has fostered large-scale ... [more ▼]

Objective: Increasing availability of surgically resected brain tissue from patients with focal epilepsy and Focal Cortical Dysplasia (FCD) or low-grade glio-neuronal tumors has fostered large-scale genetic examination. However, assessment of pathogenicity of germline and somatic variants remains difficult. Here, we present a state of the art evaluation of reported genes and variants associated with epileptic brain lesions. Methods: We critically re-evaluated the pathogenicity for all neuropathology-associated variants reported to date in PubMed and ClinVar databases including 101 neuropathology-associated missense variants encompassing 11 disease-related genes. We assessed gene variant tolerance and classified all identified missense variants according to guidelines from the American College of Medical Genetics and Genomics (ACMG). We further extended the bioinformatic variant prediction by introducing a novel gene-specific deleteriousness ranking for prediction scores. Results: Application of ACMG guidelines and in silico gene variant tolerance analysis classified only seven out of 11 genes to be likely disease-associated according to the reported a disease mechanism, while 61 (60.4%) of 101 variants of those genes were classified as of uncertain significance (VUS), 37 (36.6%) as being likely pathogenic (LP) and 3 (3%) as being pathogenic (P). Significance: We concluded that the majority of neuropathology-associated variants reported to date do not have enough evidence to be classified as pathogenic. Interpretation of lesion-associated variants remains challenging and application of current ACMG guidelines is recommended for interpretation and prediction. [less ▲]

Exome-wide analysis of mutational burden in patients with typical and atypical Rolandic Epilepsy

in European Journal of Human Genetics (2018)

Rolandic Epilepsy (RE) is the most common focal epilepsy in childhood. To date no hypothesis-free exome-wide mutational screen has been conducted for RE and Atypical RE (ARE). Here we report on whole ... [more ▼]

Rolandic Epilepsy (RE) is the most common focal epilepsy in childhood. To date no hypothesis-free exome-wide mutational screen has been conducted for RE and Atypical RE (ARE). Here we report on whole-exome sequencing of 194 unrelated patients with RE/ARE and 567 ethnically matched population controls. We identified an exome-wide significantly enriched burden for deleterious and loss-of-function variants only for the established RE/ARE gene GRIN2A. The statistical significance of the enrichment disappeared after removing ARE patients. For several disease-related gene-sets, an odds ratio > 1 was detected for loss-of-function variants. [less ▲]

Reassessment Of Lesion-Associated Gene And Variant Pathogenicity In Focal Human Epilepsies

E-print/Working paper (2017)

Purpose: Increasing availability of surgically resected brain tissue from Focal Cortical Dysplasia and low-grade epilepsy-associated tumor patients fostered large-scale genetic examination. However ... [more ▼]

Purpose: Increasing availability of surgically resected brain tissue from Focal Cortical Dysplasia and low-grade epilepsy-associated tumor patients fostered large-scale genetic examination. However, assessment of germline and somatic variant pathogenicity remains difficult. Methods: Here, we critically reevaluated the pathogenicity for all neuropathology-associated variants reported to date in the PubMed and ClinVar databases, including 12 disease-related genes and 88 neuropathology-associated missense variants. We (1) assessed evolutionary gene constraint using the pLI and missense z scores, (2) applied guidelines by the American College of Medical Genetics and Genomics (ACMG), and (3) predicted pathogenicity by using PolyPhen-2, CADD, and GERP. Results: Constraint analysis classified only seven out of 12 genes to be likely disease-associated, while 35 (40\%) of those 88 variants were classified as being variants of unknown significance (VUS) and 53 (60\%) as being likely pathogenic (LPII). Pathogenicity prediction yielded discrimination between neuropathology-associated variants (LPII and VUS) and rare variant scores obtained from individuals present in the Genome Aggregation Database (gnomAD). Conclusion: We conclude that several VUS are likely disease-associated and will be reclassified by future molecular evidence. In summary, interpretation of lesion-associated gene variants remains complex while the application of current ACMG guidelines including bioinformatic pathogenicity prediction will help improving interpretation and prediction. [less ▲]