SLC6A1 variant pathogenicity, molecular function, and phenotype: a genetic and clinical analysis

in Brain: a Journal of Neurology (2023)

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are ... [more ▼]

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are needed to understand the disease mechanism and improve therapeutic decision-making. We aggregated genetic and clinical data from 172 individuals with likely pathogenic/pathogenic (lp/p) SLC6A1 variants and functional data for 184 variants (14.1% lp/p). Clinical and functional data were available for a subset of 126 individuals. We explored the potential associations of variant positions on the GAT1 3D structure with variant pathogenicity, altered molecular function, and phenotype severity using bioinformatic approaches. The GAT1 transmembrane domains 1, 6, and extracellular loop 4 (EL4) were enriched for patient over population variants. Across functionally tested missense variants (n = 156), the spatial proximity from the ligand was associated with loss-of-function in the GAT1 transporter activity. For variants with complete loss of in vitro GABA uptake, we found a 4.6-fold enrichment in patients having severe disease vs. non-severe disease (P = 2.9e-3, 95% CI: 1.5 - 15.3). In summary, we delineated associations between the 3D structure and variant pathogenicity, variant function, and phenotype in SLC6A1-related disorders. This knowledge supports biology-informed variant interpretation and research on GAT1 function. All our data can be interactively explored in the SLC6A1 Portal (https://slc6a1-portal.broadinstitute.org/). [less ▲]

Genome-wide identification and phenotypic characterization of seizure-associated copy number variations in 741,075 individuals.

in Nature Communications (2023), 14(1), 4392

Copy number variants (CNV) are established risk factors for neurodevelopmental disorders with seizures or epilepsy. With the hypothesis that seizure disorders share genetic risk factors, we pooled CNV ... [more ▼]

Copy number variants (CNV) are established risk factors for neurodevelopmental disorders with seizures or epilepsy. With the hypothesis that seizure disorders share genetic risk factors, we pooled CNV data from 10,590 individuals with seizure disorders, 16,109 individuals with clinically validated epilepsy, and 492,324 population controls and identified 25 genome-wide significant loci, 22 of which are novel for seizure disorders, such as deletions at 1p36.33, 1q44, 2p21-p16.3, 3q29, 8p23.3-p23.2, 9p24.3, 10q26.3, 15q11.2, 15q12-q13.1, 16p12.2, 17q21.31, duplications at 2q13, 9q34.3, 16p13.3, 17q12, 19p13.3, 20q13.33, and reciprocal CNVs at 16p11.2, and 22q11.21. Using genetic data from additional 248,751 individuals with 23 neuropsychiatric phenotypes, we explored the pleiotropy of these 25 loci. Finally, in a subset of individuals with epilepsy and detailed clinical data available, we performed phenome-wide association analyses between individual CNVs and clinical annotations categorized through the Human Phenotype Ontology (HPO). For six CNVs, we identified 19 significant associations with specific HPO terms and generated, for all CNVs, phenotype signatures across 17 clinical categories relevant for epileptologists. This is the most comprehensive investigation of CNVs in epilepsy and related seizure disorders, with potential implications for clinical practice. [less ▲]

CNV-ClinViewer: Enhancing the clinical interpretation of large copy-number variants online

in Bioinformatics (2023)

Pathogenic copy number variants (CNVs) can cause a heterogeneous spectrum of rare and severe disorders. However, most CNVs are benign and are part of natural variation in human genomes. CNV pathogenicity ... [more ▼]

Pathogenic copy number variants (CNVs) can cause a heterogeneous spectrum of rare and severe disorders. However, most CNVs are benign and are part of natural variation in human genomes. CNV pathogenicity classification, genotype-phenotype analyses, and therapeutic target identification are challenging and time-consuming tasks that require the integration and analysis of information from multiple scattered sources by experts.Here, we introduce the CNV-ClinViewer, an open-source web-application for clinical evaluation and visual exploration of CNVs. The application enables real-time interactive exploration of large CNV datasets in a user-friendly designed interface and facilitates semi-automated clinical CNV interpretation following the ACMG guidelines by integrating the ClassifCNV tool. In combination with clinical judgment the application enables clinicians and researchers to formulate novel hypotheses and guide their decision-making process. Subsequently, the CNV-ClinViewer enhances for clinical investigators patient care and for basic scientists translational genomic research.The web-application is freely available at https://cnv-ClinViewer.broadinstitute.org and the open-source code can be found at https://github.com/LalResearchGroup/CNV-clinviewer.Supplementary data are available at Bioinformatics online. [less ▲]

Delineation of functionally essential protein regions for 242 neurodevelopmental genes

in Brain: a Journal of Neurology (2023), 146(2), 519-533

Neurodevelopmental disorders (NDDs), including severe pediatric epilepsy, autism, and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a ... [more ▼]

Neurodevelopmental disorders (NDDs), including severe pediatric epilepsy, autism, and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a pathogenic variant. For many of them, genetic therapies will be tested in this or the coming years in clinical trials. In contrast to first-generation symptomatic treatments, the new disease-modifying precision medicines require a genetic test-informed diagnosis before a patient can be enrolled in a clinical trial. However, even in 2022, most identified genetic variants in NDD genes are ‘Variants of Uncertain Significance’. To safely enroll patients in precision medicine clinical trials, it is important to increase our knowledge about which regions in NDD-associated proteins can ‘tolerate’ missense variants and which ones are ‘essential’ and will cause a NDD when mutated. In addition, knowledge about functionally indispensable regions in the three-dimensional (3D) structure context of proteins can also provide insights into the molecular mechanisms of disease variants. We developed a novel consensus approach that overlays evolutionary, and population based genomic scores to identify 3D essential sites (Essential3D) on protein structures. After extensive benchmarking of AlphaFold predicted and experimentally solved protein structures, we generated the currently largest expert curated protein structure set for 242 NDDs and identified 14,377 Essential3D sites across 189 gene disorders associated proteins. We demonstrate that the consensus annotation of Essential3D sites improves prioritization of disease mutations over single annotations. The identified Essential3D sites were enriched for functional features such as intermembrane regions or active sites and discovered key inter-molecule interactions in protein complexes that were otherwise not annotated. Using the currently largest autism, developmental disorders, and epilepsies exome sequencing studies including >360,000 NDD patients and population controls, we found that missense variants at Essential3D sites are 8-fold enriched in patients. In summary, we developed a comprehensive protein structure set for 242 neurodevelopmental disorders and identified 14,377 Essential3D sites in these. All data are available at https://es-ndd.broadinstitute.org for interactive visual inspection to enhance variant interpretation and development of mechanistic hypotheses for 242 NDDs genes. The provided resources will enhance clinical variant interpretation and in silico drug target development for NDD-associated genes and encoded proteins. [less ▲]

Genotype-phenotype correlations in SCN8A-related disorders reveal prognostic and therapeutic implications

in Brain (2022)

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of ... [more ▼]

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of describing clinical phenotypes related to functional effects. Six different clinical subgroups could be identified: 1) Benign familial infantile epilepsy (BFIE) (n = 15, normal cognition, treatable seizures), 2) intermediate epilepsy (n = 33, mild ID, partially pharmaco-responsive), 3) developmental and epileptic encephalopathy (DEE, n = 177, severe ID, majority pharmaco-resistant), 4) generalized epilepsy (n = 20, mild to moderate ID, frequently with absence seizures), 5) unclassifiable epilepsy (n = 127), and 6) neurodevelopmental disorder without epilepsy (n = 20, mild to moderate ID). Groups 1–3 presented with focal or multifocal seizures (median age of onset: four months) and focal epileptiform discharges, whereas the onset of seizures in group 4 was later (median: 42 months) with generalized epileptiform discharges. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin-insensitive human NaV1.6 channels and whole-cell patch-clamping. Two variants causing DEE showed a strong gain-of-function (GOF, hyperpolarising shift of steady-state activation, strongly increased neuronal firing rate), and one variant causing BFIE or intermediate epilepsy showed a mild GOF (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (LOF, reduced current amplitudes, depolarising shift of steady-state activation, reduced neuronal firing). Including previous studies, functional effects were known for 170 individuals. All 136 individuals carrying a functionally tested GOF variant had either focal (97, groups 1–3), or unclassifiable epilepsy (39), whereas 34 with a LOF variant had either generalized (14), no (11) or unclassifiable (6) epilepsy; only three had DEE. Computational modeling in the GOF group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. GOF variant carriers responded significantly better to sodium channel blockers (SCBs) than to other anti-seizure medications, and the same applied for all individuals of groups 1–3.In conclusion, our data reveal clear genotype-phenotype correlations between age at seizure onset, type of epilepsy and gain- or loss-of-function effects of SCN8A variants. Generalized epilepsy with absence seizures is the main epilepsy phenotype of LOF variant carriers and the extent of the electrophysiological dysfunction of the GOF variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that SCBs present a treatment option in SCN8A-related focal epilepsy with onset in the first year of life. [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 ▲]

Spectrum of Phenotypic, Genetic, and Functional Characteristics in Epilepsy Patients With KCNC2 Pathogenic Variants 10.1212/WNL.0000000000200660

in Neurology (2022)

Background: KCNC2 encodes Kv3.2, a member of the Shaw-related (Kv3) voltage-gated potassium channel subfamily, which is important for sustained high-frequency firing and optimized energy efficiency of ... [more ▼]

Background: KCNC2 encodes Kv3.2, a member of the Shaw-related (Kv3) voltage-gated potassium channel subfamily, which is important for sustained high-frequency firing and optimized energy efficiency of action potentials in the brain. The objective of this study was to analyse the clinical phenotype, genetic background, and biophysical function of disease-associated Kv3.2 variants.Methods: Individuals with KCNC2 variants detected by exome sequencing were selected for clinical, further genetic, and functional analysis. Cases were referred through clinical and research collaborations. Selected de novo variants were examined electrophysiologically in Xenopus laevis oocytes.Results: We identified novel KCNC2 variants in 18 patients with various forms of epilepsy including genetic generalized epilepsy (GGE), developmental and epileptic encephalopathy (DEE) including early-onset absence epilepsy (EOAE), focal epilepsy (FE), and myoclonic-atonic epilepsy (MAE). 10/18 variants were de novo and 8/18 variants were classified as modifying variants. 8 drug responsive cases became seizure-free using valproic acid as monotherapy or in combination including severe DEE cases. Functional analysis of four variants demonstrated gain-of-function in three severely affected DEE cases and loss-of-function in one case with a milder phenotype (GGE) as the underlying pathomechanisms.Conclusion: These findings implicate KCNC2 as a novel causative gene for epilepsy and emphasize the critical role of KV3.2 in the regulation of brain excitability. [less ▲]

Shared genetic basis between genetic generalized epilepsy and background electroencephalographic oscillations

in Epilepsia (2021), 62(7), 1518-1527

Abstract Objective Paroxysmal epileptiform abnormalities on electroencephalography (EEG) are the hallmark of epilepsies, but it is uncertain to what extent epilepsy and background EEG oscillations share ... [more ▼]

Abstract Objective Paroxysmal epileptiform abnormalities on electroencephalography (EEG) are the hallmark of epilepsies, but it is uncertain to what extent epilepsy and background EEG oscillations share neurobiological underpinnings. Here, we aimed to assess the genetic correlation between epilepsy and background EEG oscillations. Methods Confounding factors, including the heterogeneous etiology of epilepsies and medication effects, hamper studies on background brain activity in people with epilepsy. To overcome this limitation, we compared genetic data from a genome-wide association study (GWAS) on epilepsy (n = 12 803 people with epilepsy and 24 218 controls) with that from a GWAS on background EEG (n = 8425 subjects without epilepsy), in which background EEG oscillation power was quantified in four different frequency bands: alpha, beta, delta, and theta. We replicated our findings in an independent epilepsy replication dataset (n = 4851 people with epilepsy and 20 428 controls). To assess the genetic overlap between these phenotypes, we performed genetic correlation analyses using linkage disequilibrium score regression, polygenic risk scores, and Mendelian randomization analyses. Results Our analyses show strong genetic correlations of genetic generalized epilepsy (GGE) with background EEG oscillations, primarily in the beta frequency band. Furthermore, we show that subjects with higher beta and theta polygenic risk scores have a significantly higher risk of having generalized epilepsy. Mendelian randomization analyses suggest a causal effect of GGE genetic liability on beta oscillations. Significance Our results point to shared biological mechanisms underlying background EEG oscillations and the susceptibility for GGE, opening avenues to investigate the clinical utility of background EEG oscillations in the diagnostic workup of epilepsy. [less ▲]

Genotype-phenotype correlations in SCN8A-related disorders reveal prognostic and therapeutic implications

E-print/Working paper (2021)

We report detailed functional analyses and genotype-phenotype correlations in 433 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6. Five different ... [more ▼]

We report detailed functional analyses and genotype-phenotype correlations in 433 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6. Five different clinical subgroups could be identified: 1) Benign familial infantile epilepsy (BFIE) (n=17, normal cognition, treatable seizures), 2) intermediate epilepsy (n=36, mild ID, partially pharmacoresponsive), 3) developmental and epileptic encephalopathy (DEE, n=191, severe ID, majority pharmacoresistant), 4) generalized epilepsy (n=21, mild to moderate ID, frequently with absence seizures), and 5) affected individuals without epilepsy (n=25, mild to moderate ID). Groups 1-3 presented with early-onset (median: four months) focal or multifocal seizures and epileptic discharges, whereas the onset of seizures in group 4 was later (median: 39 months) with generalized epileptic discharges. The epilepsy was not classifiable in 143 individuals. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin insensitive human NaV1.6 channels and whole-cell patch clamping. Two variants causing DEE showed a strong gain-of-function (GOF, hyperpolarising shift of steady-state activation, strongly increased neuronal firing rate), and one variant causing BFIE or intermediate epilepsy showed a mild GOF (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (LOF, reduced current amplitudes, depolarising shift of steady-state activation, reduced neuronal firing). Including previous studies, functional effects were known for 165 individuals. All 133 individuals carrying GOF variants had either focal (76, groups 1-3), or unclassifiable epilepsy (37), whereas 32 with LOF variants had either generalized (14), no (11) or unclassifiable (5) epilepsy; only two had DEE. Computational modeling in the GOF group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. GOF variant carriers responded significantly better to sodium channel blockers (SCBs) than to other anti-seizure medications, and the same applied for all individuals of groups 1-3.In conclusion, our data reveal clear genotype-phenotype correlations between age at seizure onset, type of epilepsy and gain- or loss-of-function effects of SCN8A variants. Generalized epilepsy with absence seizures is the main epilepsy phenotype of LOF variant carriers and the extent of the electrophysiological dysfunction of the GOF variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that SCBs present a therapeutic treatment option in early onset SCN8A-related focal epilepsy. [less ▲]

Predicting functional effects of missense variants in voltage-gated sodium and calcium channels

in Science Translational Medicine (2020), 12(556), 6848

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases ... [more ▼]

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes. [less ▲]