Genotype-phenotype associations in 1018 individuals with SCN1A-related epilepsies

OBJECTIVE

SCN1A variants are associated with epilepsy syndromes ranging from mild genetic epilepsy with febrile seizures plus (GEFS+) to severe Dravet syndrome (DS). Many variants are de novo, making early phenotype prediction difficult, and genotype-phenotype associations remain poorly understood.

METHODS

We assessed data from a retrospective cohort of 1018 individuals with SCN1A-related epilepsies. We explored relationships between variant characteristics (position, in silico prediction scores: Combined Annotation Dependent Depletion (CADD), Rare Exome Variant Ensemble Learner (REVEL), SCN1A genetic score), seizure characteristics, and epilepsy phenotype.

RESULTS

DS had earlier seizure onset than other GEFS+ phenotypes (5.3 vs. 12.0 months, p < .001). In silico variant scores were higher in DS versus GEFS+ (p < .001). Patients with missense variants in functionally important regions (conserved N-terminus, S4-S6) exhibited earlier seizure onset (6.0 vs. 7.0 months, p = .003) and were more likely to have DS (280/340); those with missense variants in nonconserved regions had later onset (10.0 vs. 7.0 months, p = .036) and were more likely to have GEFS+ (15/29, χ2 = 19.16, p < .001). A minority of protein-truncating variants were associated with GEFS+ (10/393) and more likely to be located in the proximal first and last exon coding regions than elsewhere in the gene (9.7% vs. 1.0%, p < .001). Carriers of the same missense variant exhibited less variability in age at seizure onset compared with carriers of different missense variants for both DS (1.9 vs. 2.9 months, p = .001) and GEFS+ (8.0 vs. 11.0 months, p = .043). Status epilepticus as presenting seizure type is a highly specific (95.2%) but nonsensitive (32.7%) feature of DS.

SIGNIFICANCE

Understanding genotype-phenotype associations in SCN1A-related epilepsies is critical for early diagnosis and management. We demonstrate an earlier disease onset in patients with missense variants in important functional regions, the occurrence of GEFS+ truncating variants, and the value of in silico prediction scores. Status epilepticus as initial seizure type is a highly specific, but not sensitive, early feature of DS.

Polygenic burden and its association with baseline cognitive function and postoperative cognitive outcome in temporal lobe epilepsy

OBJECTIVE

Demographic and disease factors are associated with cognitive deficits and postoperative cognitive declines in adults with pharmacoresistant temporal lobe epilepsy (TLE), but the role of genetic factors in cognition in TLE is not well understood. Polygenic scores (PGS) for neurological and neuropsychiatric disorders and IQ have been associated with cognition in patient and healthy populations. In this exploratory study, we examined the relationship between PGS for Alzheimer's disease (AD), depression, and IQ and cognitive outcomes in adults with TLE.

METHODS

202 adults with pharmacoresistant TLE had genotyping and completed neuropsychological evaluations as part of a presurgical work-up. A subset (n = 116) underwent temporal lobe resection and returned for postoperative cognitive testing. Logistic regression was used to determine if PGS for AD, depression, and IQ predicted baseline domain-specific cognitive function and cognitive phenotypes as well as postoperative language and memory decline.

RESULTS

No significant findings survived correction for multiple comparisons. Prior to correction, higher PGS for AD and depression (i.e., increased genetic risk for the disorder), but lower PGS for IQ (i.e., decreased genetic likelihood of high IQ) appeared possibly associated with baseline cognitive impairment in TLE. In comparison, higher PGS for AD and IQ appeared as possible risk factors for cognitive decline following temporal lobectomy, while the possible relationship between PGS for depression and post-operative cognitive outcome was mixed.

SIGNIFICANCE

We did not observe any relationships of large effect between PGS and cognitive function or postsurgical outcome; however, results highlight several promising trends in the data that warrant future investigation in larger samples better powered to detect small genetic effects.

Healthcare utilization and clinical characteristics of genetic epilepsy in electronic health records

Understanding the clinical characteristics and medical treatment of individuals affected by genetic epilepsies is instrumental in guiding selection for genetic testing, defining the phenotype range of these rare disorders, optimizing patient care pathways and pinpointing unaddressed medical need by quantifying healthcare resource utilization. To date, a matched longitudinal cohort study encompassing the entire spectrum of clinical characteristics and medical treatment from childhood through adolescence has not been performed. We identified individuals with genetic and non-genetic epilepsies and onset at ages 0-5 years by linkage across the Cleveland Clinic Health System. We used natural language processing to extract medical terms and procedures from longitudinal electronic health records and tested for cross-sectional and temporal associations with genetic epilepsy. We implemented a two-stage design: in the discovery cohort, individuals were stratified as being 'likely genetic' or 'non-genetic' by a natural language processing algorithm, and controls did not receive genetic testing. The validation cohort consisted of cases with genetic epilepsy confirmed by manual chart review and an independent set of controls who received negative genetic testing. The discovery and validation cohorts consisted of 503 and 344 individuals with genetic epilepsy and matched controls, respectively. The median age at the first encounter was 0.1 years and 7.9 years at the last encounter, and the mean duration of follow-up was 8.2 years. We extracted 188,295 Unified Medical Language System annotations for statistical analysis across 9659 encounters. Individuals with genetic epilepsy received an earlier epilepsy diagnosis and had more frequent and complex encounters with the healthcare system. Notably, the highest enrichment of encounters compared with the non-genetic groups was found during the transition from paediatric to adult care. Our computational approach could validate established comorbidities of genetic epilepsies, such as behavioural abnormality and intellectual disability. We also revealed novel associations for genitourinary abnormalities (odds ratio 1.91, 95% confidence interval: 1.66-2.20, P = 6.16 × 10-19) linked to a spectrum of underrecognized epilepsy-associated genetic disorders. This case-control study leveraged real-world data to identify novel features associated with the likelihood of a genetic aetiology and quantified the healthcare utilization of genetic epilepsies compared with matched controls. Our results strongly recommend early genetic testing to stratify individuals into specialized care paths, thus improving the clinical management of people with genetic epilepsies.

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/).

Deep histopathology genotype-phenotype analysis of focal cortical dysplasia type II differentiates between the GATOR1-altered autophagocytic subtype IIa and MTOR-altered migration deficient subtype IIb

Focal cortical dysplasia type II (FCDII) is the most common cause of drug-resistant focal epilepsy in children. Herein, we performed a deep histopathology-based genotype-phenotype analysis to further elucidate the clinico-pathological and genetic presentation of FCDIIa compared to FCDIIb. Seventeen individuals with histopathologically confirmed diagnosis of FCD ILAE Type II and a pathogenic variant detected in brain derived DNA whole-exome sequencing or mTOR gene panel sequencing were included in this study. Clinical data were directly available from each contributing centre. Histopathological analyses were performed from formalin-fixed, paraffin-embedded tissue samples using haematoxylin-eosin and immunohistochemistry for NF-SMI32, NeuN, pS6, p62, and vimentin. Ten individuals carried loss-of-function variants in the GATOR1 complex encoding genes DEPDC5 (n = 7) and NPRL3 (n = 3), or gain-of-function variants in MTOR (n = 7). Whereas individuals with GATOR1 variants only presented with FCDIIa, i.e., lack of balloon cells, individuals with MTOR variants presented with both histopathology subtypes, FCDIIa and FCDIIb. Interestingly, 50% of GATOR1-positive cases showed a unique and predominantly vacuolizing phenotype with p62 immunofluorescent aggregates in autophagosomes. All cases with GATOR1 alterations had neurosurgery in the frontal lobe and the majority was confined to the cortical ribbon not affecting the white matter. This pattern was reflected by subtle or negative MRI findings in seven individuals with GATOR1 variants. Nonetheless, all individuals were seizure-free after surgery except four individuals carrying a DEPDC5 variant. We describe a yet underrecognized genotype-phenotype correlation of GATOR1 variants with FCDIIa in the frontal lobe. These lesions were histopathologically characterized by abnormally vacuolizing cells suggestive of an autophagy-altered phenotype. In contrast, individuals with FCDIIb and brain somatic MTOR variants showed larger lesions on MRI including the white matter, suggesting compromised neural cell migration.

Prophecy or empiricism? Clinical value of predicting versus determining genetic variant functions

The recent explosion of epilepsy genetic testing has created challenges for interpretation of gene variants. Assessments of the functional consequences of genetic variants either by predictive or experimental strategies can contribute to estimating pathogenicity, but there is no consensus on which approach is best. The Special Interest Group on Epilepsy Genetics hosted a session during the Annual American Epilepsy Society Meeting in December 2022 to discuss this topic. The session featured a debate of the relative advantages and limitations of predicting (prophecy) versus experimentally determining (empiricism) variant function using ion channel gene variants as examples. This commentary summarizes these discussions.

Classification of missense variants in the N-methyl-d-aspartate receptor GRIN gene family as gain- or loss-of-function

Advances in sequencing technology have generated a large amount of genetic data from patients with neurological conditions. These data have provided diagnosis of many rare diseases, including a number of pathogenic de novo missense variants in GRIN genes encoding N-methyl-d-aspartate receptors (NMDARs). To understand the ramifications for neurons and brain circuits affected by rare patient variants, functional analysis of the variant receptor is necessary in model systems. For NMDARs, this functional analysis needs to assess multiple properties in order to understand how variants could impact receptor function in neurons. One can then use these data to determine whether the overall actions will increase or decrease NMDAR-mediated charge transfer. Here, we describe an analytical and comprehensive framework by which to categorize GRIN variants as either gain-of-function (GoF) or loss-of-function (LoF) and apply this approach to GRIN2B variants identified in patients and the general population. This framework draws on results from six different assays that assess the impact of the variant on NMDAR sensitivity to agonists and endogenous modulators, trafficking to the plasma membrane, response time course and channel open probability. We propose to integrate data from multiple in vitro assays to arrive at a variant classification, and suggest threshold levels that guide confidence. The data supporting GoF and LoF determination are essential to assessing pathogenicity and patient stratification for clinical trials as personalized pharmacological and genetic agents that can enhance or reduce receptor function are advanced. This approach to functional variant classification can generalize to other disorders associated with missense variants.

GWAS meta-analysis of over 29,000 people with epilepsy identifies 26 risk loci and subtype-specific genetic architecture

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here we report a multi-ancestry genome-wide association study including 29,944 cases, stratified into three broad categories and seven subtypes of epilepsy, and 52,538 controls. We identify 26 genome-wide significant loci, 19 of which are specific to genetic generalized epilepsy (GGE). We implicate 29 likely causal genes underlying these 26 loci. SNP-based heritability analyses show that common variants explain between 39.6% and 90% of genetic risk for GGE and its subtypes. Subtype analysis revealed markedly different genetic architectures between focal and generalized epilepsies. Gene-set analyses of GGE signals implicate synaptic processes in both excitatory and inhibitory neurons in the brain. Prioritized candidate genes overlap with monogenic epilepsy genes and with targets of current antiseizure medications. Finally, we leverage our results to identify alternate drugs with predicted efficacy if repurposed for epilepsy treatment.

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.

Epilepsies of presumed genetic etiology show enrichment of rare variants that occur in the general population

Previous studies suggested that severe epilepsies, e.g., developmental and epileptic encephalopathies (DEEs), are mainly caused by ultra-rare de novo genetic variants. For milder disease, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 control individuals. Here, we separately analyzed three different groups of epilepsies: severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in control individuals with an allele count ≥ 1 and a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD ≥ 20), and to have an odds ratio in individuals with epilepsy ≥ 2. We identified genes enriched with QRVs primarily in NAFE (n = 72), followed by GGE (n = 32) and DEE (n = 21). This suggests that rare variants may play a more important role for causality of NAFE than for DEE. Moreover, we found that genes harboring QRVs, e.g., HSGP2, FLNA, or TNC, encode proteins that are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE that occur also in the general population, while in DEE and GGE, the contribution of such variants appears more limited.

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.

Technological and computational approaches to detect somatic mosaicism in epilepsy

Lesional epilepsy is a common and severe disease commonly associated with malformations of cortical development, including focal cortical dysplasia and hemimegalencephaly. Recent advances in sequencing and variant calling technologies have identified several genetic causes, including both short/single nucleotide and structural somatic variation. In this review, we aim to provide a comprehensive overview of the methodological advancements in this field while highlighting the unresolved technological and computational challenges that persist, including ultra-low variant allele fractions in bulk tissue, low availability of paired control samples, spatial variability of mutational burden within the lesion, and the issue of false-positive calls and validation procedures. Information from genetic testing in focal epilepsy may be integrated into clinical care to inform histopathological diagnosis, postoperative prognosis, and candidate precision therapies.

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.

Ganglioglioma with adverse clinical outcome and atypical histopathological features were defined by alterations in PTPN11/KRAS/NF1 and other RAS-/MAP-Kinase pathway genes

Exome-wide sequencing studies recently described PTPN11 as a novel brain somatic epilepsy gene. In contrast, germline mutations of PTPN11 are known to cause Noonan syndrome, a multisystem disorder characterized by abnormal facial features, developmental delay, and sporadically, also brain tumors. Herein, we performed a deep phenotype-genotype analysis of a comprehensive series of ganglioglioma (GG) with brain somatic alterations of the PTPN11/KRAS/NF1 genes compared to GG with common MAP-Kinase signaling pathway alterations, i.e., BRAFV600E. Seventy-two GG were submitted to whole exome sequencing and genotyping and 84 low grade epilepsy associated tumors (LEAT) to DNA-methylation analysis. In 28 tumours, both analyses were available from the same sample. Clinical data were retrieved from hospital files including disease onset, age at surgery, brain localization, and seizure outcome. A comprehensive histopathology staining panel was available in all cases. We identified eight GG with PTPN11 alterations, copy number variant (CNV) gains of chromosome 12, and the commonality of additional CNV gains in NF1, KRAS, FGFR4 and RHEB, as well as BRAFV600E alterations. Histopathology revealed an atypical glio-neuronal phenotype with subarachnoidal tumor spread and large, pleomorphic, and multinuclear cellular features. Only three out of eight patients with GG and PTPN11/KRAS/NF1 alterations were free of disabling-seizures 2 years after surgery (38% had Engel I). This was remarkably different from our series of GG with only BRAFV600E mutations (85% had Engel I). Unsupervised cluster analysis of DNA methylation arrays separated these tumours from well-established LEAT categories. Our data point to a subgroup of GG with cellular atypia in glial and neuronal cell components, adverse postsurgical outcome, and genetically characterized by complex alterations in PTPN11 and other RAS-/MAP-Kinase and/or mTOR signaling pathways. These findings need prospective validation in clinical practice as they argue for an adaptation of the WHO grading system in developmental, glio-neuronal tumors associated with early onset focal epilepsy.

The role of rare genetic variants enrichment in epilepsies of presumed genetic etiology

Previous studies suggested that severe epilepsies e.g., developmental and epileptic encephalopathies (DEE) are mainly caused by ultra-rare de novo genetic variants. For milder phenotypes, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 controls. Here, we separately analyzed three different groups of epilepsies : severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in controls at a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD≥20), and to have an odds ratio in epilepsy cases ≥2. We identified genes enriched with QRVs in DEE (n=21), NAFE (n=72), and GGE (n=32) - the number of enriched genes are found greatest in NAFE and least in DEE. This suggests that rare variants may play a more important role for causality of NAFE than in DEE. Moreover, we found that QRV-carrying genes e.g., HSGP2, FLNA or TNC are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE, while in DEE and GGE, the contribution of such variants appears more limited.

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.

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.

The phenotypic spectrum associated with loss-of-function variants in monogenic epilepsy genes in the general population

Variants in monogenic epilepsy genes can cause phenotypes of varying severity. For example, pathogenic variants in the SCN1A gene can cause the severe, sporadic, and drug-resistant Dravet syndrome or the milder familiar GEFS + syndrome. We hypothesized that coding variants in epilepsy-associated genes could lead to other disease-related phenotypes in the general population. We selected 127 established monogenic epilepsy genes and explored rare loss-of-function (LoF) variant associations with 3700 phenotypes across 281,850 individuals from the UK Biobank with whole-exome sequencing data. For 5.5% of epilepsy genes, we found significant associations of LoF variants with non-epilepsy phenotypes, mostly related to mental health. These findings suggest that LoF variants in epilepsy genes are associated with neurological or psychiatric phenotypes in the general population. The evidence provided may warrant further research and genetic screening of patients with atypical presentation and inform clinical care of comorbid disorders in individuals with monogenic epilepsy forms.

Clinical Features, Neuropathology, and Surgical Outcome in Patients With Refractory Epilepsy and Brain Somatic Variants in the SLC35A2 Gene

BACKGROUND AND OBJECTIVES

The SLC35A2 gene, located at chromosome Xp11.23, encodes for a uridine diphosphate-galactose transporter. We describe clinical, genetic, neuroimaging, EEG, and histopathologic findings and assess possible predictors of postoperative seizure and cognitive outcome in 47 patients with refractory epilepsy and brain somatic SLC35A2 gene variants.

METHODS

This is a retrospective multicenter study where we performed a descriptive analysis and classical hypothesis testing. We included the variables of interest significantly associated with the outcomes in the generalized linear models.

RESULTS

Two main phenotypes were associated with brain somatic SLC35A2 variants: (1) early epileptic encephalopathy (EE, 39 patients) with epileptic spasms as the predominant seizure type and moderate to severe intellectual disability and (2) drug-resistant focal epilepsy (DR-FE, 8 patients) associated with normal/borderline cognitive function and specific neuropsychological deficits. Brain MRI was abnormal in all patients with EE and in 50% of those with DR-FE. Histopathology review identified mild malformation of cortical development with oligodendroglial hyperplasia in epilepsy in 44/47 patients and was inconclusive in 3. The 47 patients harbored 42 distinct mosaic SLC35A2 variants, including 14 (33.3%) missense, 13 (30.9%) frameshift, 10 (23.8%) nonsense, 4 (9.5%) in-frame deletions/duplications, and 1 (2.4%) splicing variant. Variant allele frequencies (VAFs) ranged from 1.4% to 52.6% (mean VAF: 17.3 ± 13.5). At last follow-up (35.5 ± 21.5 months), 30 patients (63.8%) were in Engel Class I, of which 26 (55.3%) were in Class IA. Cognitive performances remained unchanged in most patients after surgery. Regression analyses showed that the probability of achieving both Engel Class IA and Class I outcomes, adjusted by age at seizure onset, was lower when the duration of epilepsy increased and higher when postoperative EEG was normal or improved. Lower brain VAF was associated with improved postoperative cognitive outcome in the analysis of associations, but this finding was not confirmed in regression analyses.

DISCUSSION

Brain somatic SLC35A2 gene variants are associated with 2 main clinical phenotypes, EE and DR-FE, and a histopathologic diagnosis of MOGHE. Additional studies will be needed to delineate any possible correlation between specific genetic variants, mutational load in the epileptogenic tissue, and surgical outcomes.

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).

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.

The genomic landscape across 474 surgically accessible epileptogenic human brain lesions

Understanding the exact molecular mechanisms involved in the etiology of epileptogenic pathologies with or without tumor activity is essential for improving treatment of drug-resistant focal epilepsy. Here, we characterize the landscape of somatic genetic variants in resected brain specimens from 474 individuals with drug-resistant focal epilepsy using deep whole-exome sequencing (>350×) and whole-genome genotyping. Across the exome, we observe a greater number of somatic single-nucleotide variants (SNV) in low-grade epilepsy-associated tumors (LEAT; 7.92 ± 5.65 SNV) than in brain tissue from malformations of cortical development (MCD; 6.11 ± 4 SNV) or hippocampal sclerosis (HS; 5.1 ± 3.04 SNV). Tumor tissues also had the largest number of likely pathogenic variant carrying cells. LEAT had the highest proportion of samples with one or more somatic copy number variants (CNV; 24.7%), followed by MCD (5.4%) and HS (4.1%). Recurring somatic whole chromosome duplications affecting Chromosome 7 (16.8%), chromosome 5 (10.9%), and chromosome 20 (9.9%) were observed among LEAT. For germline variant-associated MCD genes such as TSC2, DEPDC5, and PTEN, germline SNV were frequently identified within large loss of heterozygosity regions, supporting the recently proposed 'second hit' disease mechanism in these genes. We detect somatic variants in twelve established lesional epilepsy genes and demonstrate exome-wide statistical support for three of these in the etiology of LEAT (e.g., BRAF) and MCD (e.g., SLC35A2 and MTOR). We also identify novel significant associations for PTPN11 with LEAT and NRAS Q61 mutated protein with a complex MCD characterized by polymicrogyria and nodular heterotopia. The variants identified in NRAS are known from cancer studies to lead to hyperactivation of NRAS, which can be targeted pharmacologically. We identify large recurrent 1q21-q44 duplication including AKT3 in association with focal cortical dysplasia type 2a with hyaline astrocytic inclusions, another rare and possibly under-recognized brain lesion. The clinical genetic analyses showed that the numbers of somatic SNV across the exome and the fraction of affected cells were positively correlated with the age at seizure onset and surgery in individuals with LEAT. In summary, our comprehensive genetic screen sheds light on the genome-scale landscape of genetic variants in epileptic brain lesions, informs the design of gene panels for clinical diagnostic screening, and guides future directions for clinical implementation of epilepsy surgery genetics.

Genetic and molecular features of seizure-freedom following surgical resections for focal epilepsy: A pilot study

Objective

Seizure outcomes after brain surgery for drug-resistant epilepsy (DRE) are very heterogeneous and difficult to predict with models utilizing the current clinical, imaging, and electrophysiological variables. In this pilot study, we investigated whether genetic and molecular biomarkers (e.g., genomic, transcriptomic) can provide additional insight into differential response to surgery.

Methods

Post-operative seizure-outcomes were collected at last follow-up (>6 months) for 201 adult patients with DRE who underwent surgery between 2004 and 2020. Resected tissue was sent for miRNA sequencing (n = 132) and mRNA sequencing (n = 135). Following the selection of 10 genes (SCN1A, NBEA, PTEN, GABRA1, LGL1, DEPDC5, IL1A, ABCB1, C3, CALHM1), we investigated SNPs in those 10 genes from previously acquired exome sequencing data (n = 106). Logistic regression was performed to test for associations between individual features (mRNAs, miRNAs, and SNPs) and post-operative seizure-outcome with an exploratory FDR P < 0.25 as the threshold for significance. Post-operative time-to-seizure analyses were performed for each SNP using a Cox proportional hazards model.

Results

The majority of patients (83%) had temporal lobe epilepsy. Mean age at surgery was 38.3 years, and 56% were female. Three SNPs (rs10276036, rs11975994, rs1128503) in multi-drug resistance gene, ABCB1, were associated with post-operative seizure outcomes. Patients with alternate alleles in ABCB1 were more likely to be seizure-free at last follow-up (52–56% reduction in seizure recurrence; FDR P = 0.24). All three SNPs were in linkage disequilibrium and highly correlated with each other. Median post-operative time-to-seizure was 63 months for patients with 2 alternate alleles, 24–33 months with 1 alternate allele, and 10–11 months with 0 alternate alleles. These SNPs improved outcome prediction beyond MRI and sex alone. No independent miRNAs or mRNAs were significantly associated with seizure-outcome (P > 0.05). However, pathway analysis identified “cancer drug resistance by drug efflux” (mir-154 and mir-379) as enriched (P = 0.02), supporting the role of drug response genes in post-operative seizure recurrence.

Significance

ABCB1 may have a role in epileptogenesis and surgery outcomes independent of its drug efflux activity necessitating further investigation. SNPs in ABCB1 may serve as independent predictors of post-operative outcome.

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

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 were identified: Group 1, benign familial infantile epilepsy (n = 15, normal cognition, treatable seizures); Group 2, intermediate epilepsy (n = 33, mild intellectual disability, partially pharmaco-responsive); Group 3, developmental and epileptic encephalopathy (n = 177, severe intellectual disability, majority pharmaco-resistant); Group 4, generalized epilepsy (n = 20, mild to moderate intellectual disability, frequently with absence seizures); Group 5, unclassifiable epilepsy (n = 127); and Group 6, neurodevelopmental disorder without epilepsy (n = 20, mild to moderate intellectual disability). Those in Groups 1-3 presented with focal or multifocal seizures (median age of onset: 4 months) and focal epileptiform discharges, whereas the onset of seizures in patients with generalized epilepsy 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 developmental and epileptic encephalopathy showed a strong gain-of-function (hyperpolarizing shift of steady-state activation, strongly increased neuronal firing rate) and one variant causing benign familial infantile epilepsy or intermediate epilepsy showed a mild gain-of-function (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (reduced current amplitudes, depolarizing shift of steady-state activation, reduced neuronal firing). Functional effects were known for 170 individuals. All 136 individuals carrying a functionally tested gain-of-function variant had either focal (n = 97, Groups 1-3) or unclassifiable (n = 39) epilepsy, whereas 34 individuals with a loss-of-function variant had either generalized (n = 14), no (n = 11) or unclassifiable (n = 6) epilepsy; only three had developmental and epileptic encephalopathy. Computational modelling in the gain-of-function group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. Gain-of-function variant carriers responded significantly better to sodium channel blockers than to other anti-seizure medications, and the same applied for all individuals in 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 loss-of-function variant carriers and the extent of the electrophysiological dysfunction of the gain-of-function variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that sodium channel blockers present a treatment option in SCN8A-related focal epilepsy with onset in the first year of life.

A cross-disorder dosage sensitivity map of the human genome

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

The ILAE consensus classification of focal cortical dysplasia: An update proposed by an ad hoc task force of the ILAE diagnostic methods commission

Ongoing challenges in diagnosing focal cortical dysplasia (FCD) mandate continuous research and consensus agreement to improve disease definition and classification. An International League Against Epilepsy (ILAE) Task Force (TF) reviewed the FCD classification of 2011 to identify existing gaps and provide a timely update. The following methodology was applied to achieve this goal: a survey of published literature indexed with ((Focal Cortical Dysplasia) AND (epilepsy)) between 01/01/2012 and 06/30/2021 (n = 1349) in PubMed identified the knowledge gained since 2012 and new developments in the field. An online survey consulted the ILAE community about the current use of the FCD classification scheme with 367 people answering. The TF performed an iterative clinico-pathological and genetic agreement study to objectively measure the diagnostic gap in blood/brain samples from 22 patients suspicious for FCD and submitted to epilepsy surgery. The literature confirmed new molecular-genetic characterizations involving the mechanistic Target Of Rapamycin (mTOR) pathway in FCD type II (FCDII), and SLC35A2 in mild malformations of cortical development (mMCDs) with oligodendroglial hyperplasia (MOGHE). The electro-clinical-imaging phenotypes and surgical outcomes were better defined and validated for FCDII. Little new information was acquired on clinical, histopathological, or genetic characteristics of FCD type I (FCDI) and FCD type III (FCDIII). The survey identified mMCDs, FCDI, and genetic characterization as fields for improvement in an updated classification. Our iterative clinico-pathological and genetic agreement study confirmed the importance of immunohistochemical staining, neuroimaging, and genetic tests to improve the diagnostic yield. The TF proposes to include mMCDs, MOGHE, and "no definite FCD on histopathology" as new categories in the updated FCD classification. The histopathological classification can be further augmented by advanced neuroimaging and genetic studies to comprehensively diagnose FCD subtypes; these different levels should then be integrated into a multi-layered diagnostic scheme. This update may help to foster multidisciplinary efforts toward a better understanding of FCD and the development of novel targeted treatment options.

The role of common genetic variation in presumed monogenic epilepsies

BACKGROUND

The developmental and epileptic encephalopathies (DEEs) are the most severe group of epilepsies which co-present with developmental delay and intellectual disability (ID). DEEs usually occur in people without a family history of epilepsy and have emerged as primarily monogenic, with damaging rare mutations found in 50% of patients. Little is known about the genetic architecture of patients with DEEs in whom no pathogenic variant is identified. Polygenic risk scoring (PRS) is a method that measures a person's common genetic burden for a trait or condition. Here, we used PRS to test whether genetic burden for epilepsy is relevant in individuals with DEEs, and other forms of epilepsy with ID.

METHODS

Genetic data on 2,759 cases with DEEs, or epilepsy with ID presumed to have a monogenic basis, and 447,760 population-matched controls were analysed. We compared PRS for 'all epilepsy', 'focal epilepsy', and 'genetic generalised epilepsy' (GGE) between cases and controls. We performed pairwise comparisons between cases stratified for identifiable rare deleterious genetic variants and controls.

FINDINGS

Cases of presumed monogenic severe epilepsy had an increased PRS for 'all epilepsy' (p<0.0001), 'focal epilepsy' (p<0.0001), and 'GGE' (p=0.0002) relative to controls, which explain between 0.08% and 3.3% of phenotypic variance. PRS was increased in cases both with and without an identified deleterious variant of major effect, and there was no significant difference in PRS between the two groups.

INTERPRETATION

We provide evidence that common genetic variation contributes to the aetiology of DEEs and other forms of epilepsy with ID, even when there is a known pathogenic variant of major effect. These results provide insight into the genetic underpinnings of the severe epilepsies and warrant a shift in our understanding of the aetiology of the DEEs as complex, rather than monogenic, disorders.

FUNDING

Science foundation Ireland, Human Genome Research Institute; National Heart, Lung, and Blood Institute; German Research Foundation.

Analysing an allelic series of rare missense variants of CACNA1I in a Swedish schizophrenia cohort

CACNA1I is implicated in the susceptibility to schizophrenia by large-scale genetic association studies of single nucleotide polymorphisms. However, the channelopathy of CACNA1I in schizophrenia is unknown. CACNA1I encodes CaV3.3, a neuronal voltage-gated calcium channel that underlies a subtype of T-type current that is important for neuronal excitability in the thalamic reticular nucleus and other regions of the brain. Here, we present an extensive functional characterization of 57 naturally occurring rare and common missense variants of CACNA1I derived from a Swedish schizophrenia cohort of more than 10 000 individuals. Our analysis of this allelic series of coding CACNA1I variants revealed that reduced CaV3.3 channel current density was the dominant phenotype associated with rare CACNA1I coding alleles derived from control subjects, whereas rare CACNA1I alleles from schizophrenia patients encoded CaV3.3 channels with altered responses to voltages. CACNA1I variants associated with altered current density primarily impact the ionic channel pore and those associated with altered responses to voltage impact the voltage-sensing domain. CaV3.3 variants associated with altered voltage dependence of the CaV3.3 channel and those associated with peak current density deficits were significantly segregated across affected and unaffected groups (Fisher's exact test, P = 0.034). Our results, together with recent data from the SCHEMA (Schizophrenia Exome Sequencing Meta-Analysis) cohort, suggest that reduced CaV3.3 function may protect against schizophrenia risk in rare cases. We subsequently modelled the effect of the biophysical properties of CaV3.3 channel variants on thalamic reticular nucleus excitability and found that compared with common variants, ultrarare CaV3.3-coding variants derived from control subjects significantly decreased thalamic reticular nucleus excitability (P = 0.011). When all rare variants were analysed, there was a non-significant trend between variants that reduced thalamic reticular nucleus excitability and variants that either had no effect or increased thalamic reticular nucleus excitability across disease status. Taken together, the results of our functional analysis of an allelic series of >50 CACNA1I variants in a schizophrenia cohort reveal that loss of function of CaV3.3 is a molecular phenotype associated with reduced disease risk burden, and our approach may serve as a template strategy for channelopathies in polygenic disorders.

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.

DDGun: an untrained predictor of protein stability changes upon amino acid variants

Estimating the functional effect of single amino acid variants in proteins is fundamental for predicting the change in the thermodynamic stability, measured as the difference in the Gibbs free energy of unfolding, between the wild-type and the variant protein (ΔΔG). Here, we present the web-server of the DDGun method, which was previously developed for the ΔΔG prediction upon amino acid variants. DDGun is an untrained method based on basic features derived from evolutionary information. It is antisymmetric, as it predicts opposite ΔΔG values for direct (A → B) and reverse (B → A) single and multiple site variants. DDGun is available in two versions, one based on only sequence information and the other one based on sequence and structure information. Despite being untrained, DDGun reaches prediction performances comparable to those of trained methods. Here we make DDGun available as a web server. For the web server version, we updated the protein sequence database used for the computation of the evolutionary features, and we compiled two new data sets of protein variants to do a blind test of its performances. On these blind data sets of single and multiple site variants, DDGun confirms its prediction performance, reaching an average correlation coefficient between experimental and predicted ΔΔG of 0.45 and 0.49 for the sequence-based and structure-based versions, respectively. Besides being used for the prediction of ΔΔG, we suggest that DDGun should be adopted as a benchmark method to assess the predictive capabilities of newly developed methods. Releasing DDGun as a web-server, stand-alone program and docker image will facilitate the necessary process of method comparison to improve ΔΔG prediction.

Development and Validation of a Prediction Model for Early Diagnosis of SCN1A-Related Epilepsies

BACKGROUND AND OBJECTIVES

Pathogenic variants in the neuronal sodium channel α1 subunit gene (SCN1A) are the most frequent monogenic cause of epilepsy. Phenotypes comprise a wide clinical spectrum, including severe childhood epilepsy; Dravet syndrome, characterized by drug-resistant seizures, intellectual disability, and high mortality; and the milder genetic epilepsy with febrile seizures plus (GEFS+), characterized by normal cognition. Early recognition of a child's risk for developing Dravet syndrome vs GEFS+ is key for implementing disease-modifying therapies when available before cognitive impairment emerges. Our objective was to develop and validate a prediction model using clinical and genetic biomarkers for early diagnosis of SCN1A-related epilepsies.

METHODS

We performed a retrospective multicenter cohort study comprising data from patients with SCN1A-positive Dravet syndrome and patients with GEFS+ consecutively referred for genetic testing (March 2001-June 2020) including age at seizure onset and a newly developed SCN1A genetic score. A training cohort was used to develop multiple prediction models that were validated using 2 independent blinded cohorts. Primary outcome was the discriminative accuracy of the model predicting Dravet syndrome vs other GEFS+ phenotypes.

RESULTS

A total of 1,018 participants were included. The frequency of Dravet syndrome was 616/743 (83%) in the training cohort, 147/203 (72%) in validation cohort 1, and 60/72 (83%) in validation cohort 2. A high SCN1A genetic score (133.4 [SD 78.5] vs 52.0 [SD 57.5]; p < 0.001) and young age at onset (6.0 [SD 3.0] vs 14.8 [SD 11.8] months; p < 0.001) were each associated with Dravet syndrome vs GEFS+. A combined SCN1A genetic score and seizure onset model separated Dravet syndrome from GEFS+ more effectively (area under the curve [AUC] 0.89 [95% CI 0.86-0.92]) and outperformed all other models (AUC 0.79-0.85; p < 0.001). Model performance was replicated in both validation cohorts 1 (AUC 0.94 [95% CI 0.91-0.97]) and 2 (AUC 0.92 [95% CI 0.82-1.00]).

DISCUSSION

The prediction model allows objective estimation at disease onset whether a child will develop Dravet syndrome vs GEFS+, assisting clinicians with prognostic counseling and decisions on early institution of precision therapies (http://scn1a-prediction-model.broadinstitute.org/).

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that a combined SCN1A genetic score and seizure onset model distinguishes Dravet syndrome from other GEFS+ phenotypes.

Structural mapping of GABRB3 variants reveals genotype-phenotype correlations

PURPOSE

Pathogenic variants in GABRB3 have been associated with a spectrum of phenotypes from severe developmental disorders and epileptic encephalopathies to milder epilepsy syndromes and mild intellectual disability (ID). In this study, we analyzed a large cohort of individuals with GABRB3 variants to deepen the phenotypic understanding and investigate genotype-phenotype correlations.

METHODS

Through an international collaboration, we analyzed electro-clinical data of unpublished individuals with variants in GABRB3, and we reviewed previously published cases. All missense variants were mapped onto the 3-dimensional structure of the GABRB3 subunit, and clinical phenotypes associated with the different key structural domains were investigated.

RESULTS

We characterized 71 individuals with GABRB3 variants, including 22 novel subjects, expressing a wide spectrum of phenotypes. Interestingly, phenotypes correlated with structural locations of the variants. Generalized epilepsy, with a median age at onset of 12 months, and mild-to-moderate ID were associated with variants in the extracellular domain. Focal epilepsy with earlier onset (median: age 4 months) and severe ID were associated with variants in both the pore-lining helical transmembrane domain and the extracellular domain.

CONCLUSION

These genotype-phenotype correlations will aid the genetic counseling and treatment of individuals affected by GABRB3-related disorders. Future studies may reveal whether functional differences underlie the phenotypic differences.

Identification and quantification of oligogenic loss-of-function disorders

PURPOSE

Monogenic disorders can present clinically heterogeneous symptoms. We hypothesized that in patients with a monogenic disorder caused by a large deletion, frequently additional loss-of-function (LOF)-intolerant genes are affected, potentially contributing to the phenotype.

METHODS

We investigated the LOF-intolerant gene distribution across the genome and its association with benign population and pathogenic classified deletions from individuals with presumably monogenic disorders. For people with presumably monogenic epilepsy, we compared Human Phenotype Ontology terms in people with large and small deletions.

RESULTS

We identified LOF-intolerant gene dense regions that were enriched for ClinVar and depleted for population copy number variants. Analysis of data from >143,000 individuals with a suspected monogenic disorder showed that 2.5% of haploinsufficiency disorder-associated deletions can affect at least 1 other LOF-intolerant gene. Focusing on epilepsy, we observed that 13.1% of pathogenic and likely pathogenic ClinVar deletions <3 megabase pair, covering the diagnostically most relevant genes, affected at least 1 additional LOF-intolerant gene. Those patients have potentially more complex phenotypes with increasing deletion size.

CONCLUSION

We could systematically show that large deletions frequently affected admditional LOF-intolerant genes in addition to the established disease gene. Further research is needed to understand how additional potential disease-relevant genes influence monogenic disorders to improve clinical care and the efficacy of targeted therapies.

Common molecular mechanisms of SLC6A1 variant-mediated neurodevelopmental disorders in astrocytes and neurons

Solute carrier family 6 member 1 (SLC6A1) is abundantly expressed in the developing brain even before the CNS is formed. Its encoded GABA transporter 1 (GAT-1) is responsible for the reuptake of GABA into presynaptic neurons and glia, thereby modulating neurotransmission. GAT-1 is expressed globally in the brain, in both astrocytes and neurons. The GABA uptake function of GAT-1 in neurons cannot be compensated for by other GABA transporters, while the function in glia can be partially replaced by GABA transporter 3. Recently, many variants in SLC6A1 have been associated with a spectrum of epilepsy syndromes and neurodevelopmental disorders, including myoclonic atonic epilepsy, childhood absence epilepsy, autism, and intellectual disability, but the pathomechanisms associated with these phenotypes remain unclear. The presence of GAT-1 in both neurons and astrocytes further obscures the role of abnormal GAT-1 in the heterogeneous disease phenotype manifestations. Here we examine the impact on transporter trafficking and function of 22 SLC6A1 variants identified in patients with a broad spectrum of phenotypes. We also evaluate changes in protein expression and subcellular localization of the variant GAT-1 in various cell types, including neurons and astrocytes derived from human patient induced pluripotent stem cells. We found that a partial or complete loss-of-function represents a common disease mechanism, although the extent of GABA uptake reduction is variable. The reduced GABA uptake appears to be due to reduced cell surface expression of the variant transporter caused by variant protein misfolding, endoplasmic reticulum retention, and subsequent degradation. Although the extent of reduction of the total protein, surface protein, and the GABA uptake level of the variant transporters is variable, the loss of GABA uptake function and endoplasmic reticulum retention is consistent across induced pluripotent stem cell-derived cell types, including astrocytes and neurons, for the surveyed variants. Interestingly, we did not find a clear correlation of GABA uptake function and the disease phenotypes, such as myoclonic atonic epilepsy versus developmental delay, in this study. Together, our study suggests that impaired transporter protein trafficking and surface expression are the major disease-associated mechanisms associated with pathogenic SLC6A1 variants. Our results resemble findings from pathogenic variants in other genes affecting the GABA pathway, such as GABAA receptors. This study provides critical insight into therapeutic developments for SLC6A1 variant-mediated disorders and implicates that boosting transporter function by either genetic or pharmacological approaches would be beneficial.

Application of single cell genomics to focal epilepsies: A call to action

Focal epilepsies are the largest epilepsy subtype and associated with significant morbidity. Somatic variation is a newly recognized genetic mechanism underlying a subset of focal epilepsies, but little is known about the processes through which somatic mosaicism causes seizures, the cell types carrying the pathogenic variants, or their developmental origin. Meanwhile, the inception of single cell biology has completely revolutionized the study of neurological diseases and has the potential to answer some of these key questions. Focusing on single cell genomics, transcriptomics, and epigenomics in focal epilepsy research, circumvents the averaging artifact associated with studying bulk brain tissue and offers the kind of granularity that is needed for investigating the consequences of somatic mosaicism. Here we have provided a brief overview of some of the most developed single cell techniques and the major considerations around applying them to focal epilepsy research.

SimText: A text mining framework for interactive analysis and visualization of similarities among biomedical entities

Summary

Literature exploration in PubMed on a large number of biomedical entities (e.g. genes, diseases or experiments) can be time-consuming and challenging, especially when assessing associations between entities. Here, we describe SimText, a user-friendly toolset that provides customizable and systematic workflows for the analysis of similarities among a set of entities based on text. SimText can be used for (i) text collection from PubMed and extraction of words with different text mining approaches, and (ii) interactive analysis and visualization of data using unsupervised learning techniques in an interactive app.

Availability and implementation

We developed SimText as an open-source R software and integrated it into Galaxy (https://usegalaxy.eu), an online data analysis platform with supporting self-learning training material available at https://training.galaxyproject.org. A command-line version of the toolset is available for download from GitHub (https://github.com/dlal-group/simtext) or as Docker image (https://hub.docker.com/r/dlalgroup/simtext/tags.).

Supplementary information

Supplementary data are available at Bioinformatics online.

Shared genetic basis between genetic generalized epilepsy and background electroencephalographic oscillations

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.

Toward a better definition of focal cortical dysplasia: An iterative histopathological and genetic agreement trial

OBJECTIVE

Focal cortical dysplasia (FCD) is a major cause of difficult-to-treat epilepsy in children and young adults, and the diagnosis is currently based on microscopic review of surgical brain tissue using the International League Against Epilepsy classification scheme of 2011. We developed an iterative histopathological agreement trial with genetic testing to identify areas of diagnostic challenges in this widely used classification scheme.

METHODS

Four web-based digital pathology trials were completed by 20 neuropathologists from 15 countries using a consecutive series of 196 surgical tissue blocks obtained from 22 epilepsy patients at a single center. Five independent genetic laboratories performed screening or validation sequencing of FCD-relevant genes in paired brain and blood samples from the same 22 epilepsy patients.

RESULTS

Histopathology agreement based solely on hematoxylin and eosin stainings was low in Round 1, and gradually increased by adding a panel of immunostainings in Round 2 and the Delphi consensus method in Round 3. Interobserver agreement was good in Round 4 (kappa = .65), when the results of genetic tests were disclosed, namely, MTOR, AKT3, and SLC35A2 brain somatic mutations in five cases and germline mutations in DEPDC5 and NPRL3 in two cases.

SIGNIFICANCE

The diagnoses of FCD 1 and 3 subtypes remained most challenging and were often difficult to differentiate from a normal homotypic or heterotypic cortical architecture. Immunohistochemistry was helpful, however, to confirm the diagnosis of FCD or no lesion. We observed a genotype-phenotype association for brain somatic mutations in SLC35A2 in two cases with mild malformation of cortical development with oligodendroglial hyperplasia in epilepsy. Our results suggest that the current FCD classification should recognize a panel of immunohistochemical stainings for a better histopathological workup and definition of FCD subtypes. We also propose adding the level of genetic findings to obtain a comprehensive, reliable, and integrative genotype-phenotype diagnosis in the near future.

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.

Genome-Wide Analysis of Copy Number Variation in Latin American Parkinson's Disease Patients

BACKGROUND

Parkinson's disease is the second most common neurodegenerative disorder and affects people from all ethnic backgrounds, yet little is known about the genetics of Parkinson's disease in non-European populations. In addition, the overall identification of copy number variants at a genome-wide level has been understudied in Parkinson's patients. The objective of this study was to understand the genome-wide burden of copy number variants in Latinos and its association with Parkinson's disease.

METHODS

We used genome-wide genotyping data from 747 Parkinson's disease patients and 632 controls from the Latin American Research Consortium on the Genetics of Parkinson's disease.

RESULTS

Genome-wide copy number burden analysis showed that patients were significantly enriched for copy number variants overlapping known Parkinson's disease genes compared with controls (odds ratio, 3.97; 95%CI, 1.69-10.5; P = 0.018). PRKN showed the strongest copy number burden, with 20 copy number variant carriers. These patients presented an earlier age of disease onset compared with patients with other copy number variants (median age at onset, 31 vs 57 years, respectively; P = 7.46 × 10-7 ).

CONCLUSIONS

We found that although overall genome-wide copy number variant burden was not significantly different, Parkinson's disease patients were significantly enriched with copy number variants affecting known Parkinson's disease genes. We also identified that of 250 patients with early-onset disease, 5.6% carried a copy number variant on PRKN in our cohort. Our study is the first to analyze genome-wide copy number variant association in Latino Parkinson's disease patients and provides insights about this complex disease in this understudied population. © 2020 International Parkinson and Movement Disorder Society.

Epilepsy subtype-specific copy number burden observed in a genome-wide study of 17 458 subjects

Cytogenic testing is routinely applied in most neurological centres for severe paediatric epilepsies. However, which characteristics of copy number variants (CNVs) confer most epilepsy risk and which epilepsy subtypes carry the most CNV burden, have not been explored on a genome-wide scale. Here, we present the largest CNV investigation in epilepsy to date with 10 712 European epilepsy cases and 6746 ancestry-matched controls. Patients with genetic generalized epilepsy, lesional focal epilepsy, non-acquired focal epilepsy, and developmental and epileptic encephalopathy were included. All samples were processed with the same technology and analysis pipeline. All investigated epilepsy types, including lesional focal epilepsy patients, showed an increase in CNV burden in at least one tested category compared to controls. However, we observed striking differences in CNV burden across epilepsy types and investigated CNV categories. Genetic generalized epilepsy patients have the highest CNV burden in all categories tested, followed by developmental and epileptic encephalopathy patients. Both epilepsy types also show association for deletions covering genes intolerant for truncating variants. Genome-wide CNV breakpoint association showed not only significant loci for genetic generalized and developmental and epileptic encephalopathy patients but also for lesional focal epilepsy patients. With a 34-fold risk for developing genetic generalized epilepsy, we show for the first time that the established epilepsy-associated 15q13.3 deletion represents the strongest risk CNV for genetic generalized epilepsy across the whole genome. Using the human interactome, we examined the largest connected component of the genes overlapped by CNVs in the four epilepsy types. We observed that genetic generalized epilepsy and non-acquired focal epilepsy formed disease modules. In summary, we show that in all common epilepsy types, 1.5-3% of patients carry epilepsy-associated CNVs. The characteristics of risk CNVs vary tremendously across and within epilepsy types. Thus, we advocate genome-wide genomic testing to identify all disease-associated types of CNVs.