Distances from ligands as main predictive features for pathogenicity and functional effect of variants in NMDA receptors.

Genetic variants in genes GRIN1 , GRIN2A , GRIN2B , and GRIN2D , which encode subunits of the N-methyl-D-aspartate receptor (NMDAR), have been associated with severe and heterogeneous neurologic diseases. Missense variants in these genes can result in gain or loss of the NMDAR function, requiring opposite therapeutic treatments. Computational methods that predict pathogenicity and molecular functional effects are therefore crucial for accurate diagnosis and therapeutic applications. We assembled missense variants: 201 from patients, 631 from general population, and 159 characterized by electrophysiological readouts showing whether they can enhance or reduce the receptor function. This includes new functional data from 47 variants reported here, for the first time. We found that pathogenic/benign variants and variants that increase/decrease the channel function were distributed unevenly on the protein structure, with spatial proximity to ligands bound to the agonist and antagonist binding sites being key predictive features. Leveraging distances from ligands, we developed two independent machine learning-based predictors for NMDAR missense variants: a pathogenicity predictor which outperforms currently available predictors (AUC=0.945, MCC=0.726), and the first binary predictor of molecular function (increase or decrease) (AUC=0.809, MCC=0.523). Using these, we reclassified variants of uncertain significance in the ClinVar database and refined a previous genome-informed epidemiological model to estimate the birth incidence of molecular mechanism-defined GRIN disorders. Our findings demonstrate that distance from ligands is an important feature in NMDARs that can enhance variant pathogenicity prediction and enable functional prediction. Further studies with larger numbers of phenotypically and functionally characterized variants will enhance the potential clinical utility of this method.

Severe communication delays are independent of seizure burden and persist despite contemporary treatments in SCN1A+ Dravet syndrome: Insights from the ENVISION natural history study.

OBJECTIVE: Dravet syndrome (DS) is a developmental and epileptic encephalopathy characterized by high seizure burden, treatment-resistant epilepsy, and developmental stagnation. Family members rate communication deficits among the most impactful disease manifestations. We evaluated seizure burden and language/communication development in children with DS. METHODS: ENVISION was a prospective, observational study evaluating children with DS associated with SCN1A pathogenic variants (SCN1A+ DS) enrolled at age ≤5 years. Seizure burden and antiseizure medications were assessed every 3 months and communication and language every 6 months with the Bayley Scales of Infant and Toddler Development 3rd edition and the parent-reported Vineland Adaptive Behavior Scales 3rd edition. We report data from the first year of observation, including analyses stratified by age at Baseline: 0:6-2:0 years:months (Y:M; youngest), 2:1-3:6 Y:M (middle), and 3:7-5:0 Y:M (oldest). RESULTS: Between December 2020 and March 2023, 58 children with DS enrolled at 16 sites internationally. Median follow-up was 17.5 months (range = .0-24.0), with 54 of 58 (93.1%) followed for at least 6 months and 51 of 58 (87.9%) for 12 months. Monthly countable seizure frequency (MCSF) increased with age (median [minimum-maximum] = 1.0 in the youngest [1.0-70.0] and middle [1.0-242.0] age groups and 4.5 [.0-2647.0] in the oldest age group), and remained high, despite use of currently approved antiseizure medications. Language/communication delays were observed early, and developmental stagnation occurred after age 2 years with both instruments. In predictive modeling, chronologic age was the only significant covariate of seizure frequency (effect size = .52, p = .024). MCSF, number of antiseizure medications, age at first seizure, and convulsive status epilepticus were not predictors of language/communication raw scores. SIGNIFICANCE: In infants and young children with SCN1A+ DS, language/communication delay and stagnation were independent of seizure burden. Our findings emphasize that the optimal therapeutic window to prevent language/communication delay is before 3 years of age.

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

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 versus non-severe disease (P = 2.9 × 10-3, 95% confidence interval: 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/).

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

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.

Scanning mutagenesis of the voltage-gated sodium channel NaV1.2 using base editing.

It is challenging to apply traditional mutational scanning to voltage-gated sodium channels (NaVs) and functionally annotate the large number of coding variants in these genes. Using a cytosine base editor and a pooled viability assay, we screen a library of 368 guide RNAs (gRNAs) tiling NaV1.2 to identify more than 100 gRNAs that change NaV1.2 function. We sequence base edits made by a subset of these gRNAs to confirm specific variants that drive changes in channel function. Electrophysiological characterization of these channel variants validates the screen results and provides functional mechanisms of channel perturbation. Most of the changes caused by these gRNAs are classifiable as loss of function along with two missense mutations that lead to gain of function in NaV1.2 channels. This two-tiered strategy to functionally characterize ion channel protein variants at scale identifies a large set of loss-of-function mutations in NaV1.2.

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

MOTIVATION: 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. RESULTS: 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. AVAILABILITY AND IMPLEMENTATION: 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.

Delineation of functionally essential protein regions for 242 neurodevelopmental genes.

Neurodevelopmental disorders (NDDs), including severe paediatric 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 enrol 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 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 NDDs 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.

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

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. As a consequence, 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 3049 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 centre 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 receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1422 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 versus 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 readouts 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.

The gain of function SCN1A disorder spectrum: novel epilepsy phenotypes and therapeutic implications.

Brain voltage-gated sodium channel NaV1.1 (SCN1A) loss-of-function variants cause the severe epilepsy Dravet syndrome, as well as milder phenotypes associated with genetic epilepsy with febrile seizures plus. Gain of function SCN1A variants are associated with familial hemiplegic migraine type 3. Novel SCN1A-related phenotypes have been described including early infantile developmental and epileptic encephalopathy with movement disorder, and more recently neonatal presentations with arthrogryposis. Here we describe the clinical, genetic and functional evaluation of affected individuals. Thirty-five patients were ascertained via an international collaborative network using a structured clinical questionnaire and from the literature. We performed whole-cell voltage-clamp electrophysiological recordings comparing sodium channels containing wild-type versus variant NaV1.1 subunits. Findings were related to Dravet syndrome and familial hemiplegic migraine type 3 variants. We identified three distinct clinical presentations differing by age at onset and presence of arthrogryposis and/or movement disorder. The most severely affected infants (n = 13) presented with congenital arthrogryposis, neonatal onset epilepsy in the first 3 days of life, tonic seizures and apnoeas, accompanied by a significant movement disorder and profound intellectual disability. Twenty-one patients presented later, between 2 weeks and 3 months of age, with a severe early infantile developmental and epileptic encephalopathy and a movement disorder. One patient presented after 3 months with developmental and epileptic encephalopathy only. Associated SCN1A variants cluster in regions of channel inactivation associated with gain of function, different to Dravet syndrome variants (odds ratio = 17.8; confidence interval = 5.4-69.3; P = 1.3 × 10-7). Functional studies of both epilepsy and familial hemiplegic migraine type 3 variants reveal alterations of gating properties in keeping with neuronal hyperexcitability. While epilepsy variants result in a moderate increase in action current amplitude consistent with mild gain of function, familial hemiplegic migraine type 3 variants induce a larger effect on gating properties, in particular the increase of persistent current, resulting in a large increase of action current amplitude, consistent with stronger gain of function. Clinically, 13 out of 16 (81%) gain of function variants were associated with a reduction in seizures in response to sodium channel blocker treatment (carbamazepine, oxcarbazepine, phenytoin, lamotrigine or lacosamide) without evidence of symptom exacerbation. Our study expands the spectrum of gain of function SCN1A-related epilepsy phenotypes, defines key clinical features, provides novel insights into the underlying disease mechanisms between SCN1A-related epilepsy and familial hemiplegic migraine type 3, and identifies sodium channel blockers as potentially efficacious therapies. Gain of function disease should be considered in early onset epilepsies with a pathogenic SCN1A variant and non-Dravet syndrome phenotype.

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

BACKGROUND AND OBJECTIVES: 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 analyze 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, focal epilepsy, and myoclonic-atonic epilepsy. Of the 18 variants, 10 were de novo and 8 were classified as modifying variants. Eight drug-responsive patients became seizure-free using valproic acid as monotherapy or in combination, including severe DEE cases. Functional analysis of 4 variants demonstrated gain of function in 3 severely affected DEE cases and loss of function in 1 case with a milder phenotype (GGE) as the underlying pathomechanisms. DISCUSSION: 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.

Gene variant effects across sodium channelopathies predict function and guide precision therapy.

Pathogenic variants in the voltage-gated sodium channel gene family lead to early onset epilepsies, neurodevelopmental disorders, skeletal muscle channelopathies, peripheral neuropathies and cardiac arrhythmias. Disease-associated variants have diverse functional effects ranging from complete loss-of-function to marked gain-of-function. Therapeutic strategy is likely to depend on functional effect. Experimental studies offer important insights into channel function but are resource intensive and only performed in a minority of cases. Given the evolutionarily conserved nature of the sodium channel genes, we investigated whether similarities in biophysical properties between different voltage-gated sodium channels can predict function and inform precision treatment across sodium channelopathies. We performed a systematic literature search identifying functionally assessed variants in any of the nine voltage-gated sodium channel genes until 28 April 2021. We included missense variants that had been electrophysiologically characterized in mammalian cells in whole-cell patch-clamp recordings. We performed an alignment of linear protein sequences of all sodium channel genes and correlated variants by their overall functional effect on biophysical properties. Of 951 identified records, 437 sodium channel-variants met our inclusion criteria and were reviewed for functional properties. Of these, 141 variants were epilepsy-associated (SCN1/2/3/8A), 79 had a neuromuscular phenotype (SCN4/9/10/11A), 149 were associated with a cardiac phenotype (SCN5/10A) and 68 (16%) were considered benign. We detected 38 missense variant pairs with an identical disease-associated variant in a different sodium channel gene. Thirty-five out of 38 of those pairs resulted in similar functional consequences, indicating up to 92% biophysical agreement between corresponding sodium channel variants (odds ratio = 11.3; 95% confidence interval = 2.8 to 66.9; P < 0.001). Pathogenic missense variants were clustered in specific functional domains, whereas population variants were significantly more frequent across non-conserved domains (odds ratio = 18.6; 95% confidence interval = 10.9-34.4; P < 0.001). Pore-loop regions were frequently associated with loss-of-function variants, whereas inactivation sites were associated with gain-of-function (odds ratio = 42.1, 95% confidence interval = 14.5-122.4; P < 0.001), whilst variants occurring in voltage-sensing regions comprised a range of gain- and loss-of-function effects. Our findings suggest that biophysical characterisation of variants in one SCN-gene can predict channel function across different SCN-genes where experimental data are not available. The collected data represent the first gain- versus loss-of-function topological map of SCN proteins indicating shared patterns of biophysical effects aiding variant analysis and guiding precision therapy. We integrated our findings into a free online webtool to facilitate functional sodium channel gene variant interpretation (http://SCN-viewer.broadinstitute.org).

ATP1A3-Related Disorders: An Ever-Expanding Clinical Spectrum.

The Na+/K+ ATPases are Sodium-Potassium exchanging pumps, with a heteromeric α-ß-γ protein complex. The α3 isoform is required as a rescue pump, after repeated action potentials, with a distribution predominantly in neurons of the central nervous system. This isoform is encoded by the ATP1A3 gene. Pathogenic variants in this gene have been implicated in several phenotypes in the last decades. Carriers of pathogenic variants in this gene manifest neurological and non-neurological features in many combinations, usually with an acute onset and paroxysmal episodes triggered by fever or other factors. The first three syndromes described were: (1) rapid-onset dystonia parkinsonism; (2) alternating hemiplegia of childhood; and, (3) cerebellar ataxia, pes cavus, optic atrophy, and sensorineural hearing loss (CAPOS syndrome). Since their original description, an expanding number of cases presenting with atypical and overlapping features have been reported. Because of this, ATP1A3-disorders are now beginning to be viewed as a phenotypic continuum representing discrete expressions along a broadly heterogeneous clinical spectrum.

The spectrum of intermediate SCN8A-related epilepsy.

OBJECTIVE: Pathogenic variants in SCN8A have been associated with a wide spectrum of epilepsy phenotypes, ranging from benign familial infantile seizures (BFIS) to epileptic encephalopathies with variable severity. Furthermore, a few patients with intellectual disability (ID) or movement disorders without epilepsy have been reported. The vast majority of the published SCN8A patients suffer from severe developmental and epileptic encephalopathy (DEE). In this study, we aimed to provide further insight on the spectrum of milder SCN8A-related epilepsies. METHODS: A cohort of 1095 patients were screened using a next generation sequencing panel. Further patients were ascertained from a network of epilepsy genetics clinics. Patients with severe DEE and BFIS were excluded from the study. RESULTS: We found 36 probands who presented with an SCN8A-related epilepsy and normal intellect (33%) or mild (61%) to moderate ID (6%). All patients presented with epilepsy between age 1.5 months and 7 years (mean = 13.6 months), and 58% of these became seizure-free, two-thirds on monotherapy. Neurological disturbances included ataxia (28%) and hypotonia (19%) as the most prominent features. Interictal electroencephalogram was normal in 41%. Several recurrent variants were observed, including Ile763Val, Val891Met, Gly1475Arg, Gly1483Lys, Phe1588Leu, Arg1617Gln, Ala1650Val/Thr, Arg1872Gln, and Asn1877Ser. SIGNIFICANCE: With this study, we explore the electroclinical features of an intermediate SCN8A-related epilepsy with mild cognitive impairment, which is for the majority a treatable epilepsy.