Ultra-Rare Genetic Variation in the Epilepsies: A Whole-Exome Sequencing Study of 17,606 Individuals

Breast cancer risks associated with missense variants in breast cancer susceptibility genes

BACKGROUND

Protein truncating variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2 are associated with increased breast cancer risk, but risks associated with missense variants in these genes are uncertain.

METHODS

We analyzed data on 59,639 breast cancer cases and 53,165 controls from studies participating in the Breast Cancer Association Consortium BRIDGES project. We sampled training (80%) and validation (20%) sets to analyze rare missense variants in ATM (1146 training variants), BRCA1 (644), BRCA2 (1425), CHEK2 (325), and PALB2 (472). We evaluated breast cancer risks according to five in silico prediction-of-deleteriousness algorithms, functional protein domain, and frequency, using logistic regression models and also mixture models in which a subset of variants was assumed to be risk-associated.

RESULTS

The most predictive in silico algorithms were Helix (BRCA1, BRCA2 and CHEK2) and CADD (ATM). Increased risks appeared restricted to functional protein domains for ATM (FAT and PIK domains) and BRCA1 (RING and BRCT domains). For ATM, BRCA1, and BRCA2, data were compatible with small subsets (approximately 7%, 2%, and 0.6%, respectively) of rare missense variants giving similar risk to those of protein truncating variants in the same gene. For CHEK2, data were more consistent with a large fraction (approximately 60%) of rare missense variants giving a lower risk (OR 1.75, 95% CI (1.47-2.08)) than CHEK2 protein truncating variants. There was little evidence for an association with risk for missense variants in PALB2. The best fitting models were well calibrated in the validation set.

CONCLUSIONS

These results will inform risk prediction models and the selection of candidate variants for functional assays and could contribute to the clinical reporting of gene panel testing for breast cancer susceptibility.

ATP6V0C Is Associated With Febrile Seizures and Epilepsy With Febrile Seizures Plus

Purpose

To identify novel genetic causes of febrile seizures (FS) and epilepsy with febrile seizures plus (EFS+).

Methods

We performed whole-exome sequencing in a cohort of 32 families, in which at least two individuals were affected by FS or EFS+. The probands, their parents, and available family members were recruited to ascertain whether the genetic variants were co-segregation. Genes with repetitively identified variants with segregations were selected for further studies to define the gene-disease association.

Results

We identified two heterozygous ATP6V0C mutations (c.64G > A/p.Ala22Thr and c.361_373del/p.Thr121Profs*7) in two unrelated families with six individuals affected by FS or EFS+. The missense mutation was located in the proteolipid c-ring that cooperated with a-subunit forming the hemichannel for proton transferring. It also affected the hydrogen bonds with surround residues and the protein stability, implying a damaging effect. The frameshift mutation resulted in a loss of function by yielding a premature termination of 28 residues at the C-terminus of the protein. The frequencies of ATP6V0C mutations identified in this cohort were significantly higher than that in the control populations. All the six affected individuals suffered from their first FS at the age of 7-8 months. The two probands later manifested afebrile seizures including myoclonic seizures that responded well to lamotrigine. They all displayed favorable outcomes without intellectual or developmental abnormalities, although afebrile seizures or frequent seizures occurred.

Conclusion

This study suggests that ATP6V0C is potentially a candidate pathogenic gene of FS and EFS+. Screening for ATP6V0C mutations would help differentiating patients with Dravet syndrome caused by SCN1A mutations, which presented similar clinical manifestation but different responses to antiepileptic treatment.

DEPDC5-related epilepsy: A comprehensive review

DEPDC5-related epilepsy, caused by pathogenic germline variants(with or without additional somatic variants in the brain) of DEPDC5 (Dishevelled, Egl-10 and Pleckstrin domain-containing protein 5) gene, is a newly discovered predominantly focal epilepsy linked to enhanced mTORC1 pathway. DEPDC5-related epilepsy includes several familial epilepsy syndromes, including familial focal epilepsy with variable foci (FFEVF) and rare sporadic nonlesional focal epilepsy. DEPDC5 has been identified as one of the more common epilepsy genes linked to infantile spasms and sudden unexpected death (SUDEP). Although intelligence usually is unaffected in DEPDC5-related epilepsy, some people have been diagnosed with intellectual disabilities, autism spectrum disorder, and other psychiatric problems. DEPDC5 variants have also been found in 20% of individuals with various brain abnormalities, challenging the traditional distinction between lesional and nonlesional epilepsies. The most exciting development of DEPDC5 variants is the possibility of precision therapeutics using mTOR inhibitors, as evidenced with phenotypic rescue in many animal models. However, more research is needed to better understand the functional impact of diverse (particularly missense or splice-region) variants, the specific involvement of DEPDC5 in epileptogenesis, and the creation and utilization of precision therapies in humans. Precision treatments for DEPDC5-related epilepsy will benefit not only a small number of people with the condition, but they will also pave the way for new therapeutic approaches in epilepsy (including acquired epilepsies in which mTORC1 activation occurs, for example, post-traumatic epilepsy) and other neurological disorders involving a dysfunctional mTOR pathway.

Association of ultra-rare coding variants with genetic generalized epilepsy: A case-control whole exome sequencing study

OBJECTIVE

We aimed to identify genes associated with genetic generalized epilepsy (GGE) by combining large cohorts enriched with individuals with a positive family history. Secondarily, we set out to compare the association of genes independently with familial and sporadic GGE.

METHODS

We performed a case-control whole exome sequencing study in unrelated individuals of European descent diagnosed with GGE (previously recruited and sequenced through multiple international collaborations) and ancestry-matched controls. The association of ultra-rare variants (URVs; in 18 834 protein-coding genes) with epilepsy was examined in 1928 individuals with GGE (vs. 8578 controls), then separately in 945 individuals with familial GGE (vs. 8626 controls), and finally in 1005 individuals with sporadic GGE (vs. 8621 controls). We additionally examined the association of URVs with familial and sporadic GGE in two gene sets important for inhibitory signaling (19 genes encoding γ-aminobutyric acid type A [GABAA ] receptors, 113 genes representing the GABAergic pathway).

RESULTS

GABRG2 was associated with GGE (p = 1.8 × 10-5 ), approaching study-wide significance in familial GGE (p = 3.0 × 10-6 ), whereas no gene approached a significant association with sporadic GGE. Deleterious URVs in the most intolerant subgenic regions in genes encoding GABAA receptors were associated with familial GGE (odds ratio [OR] = 3.9, 95% confidence interval [CI] = 1.9-7.8, false discovery rate [FDR]-adjusted p = .0024), whereas their association with sporadic GGE had marginally lower odds (OR = 3.1, 95% CI = 1.3-6.7, FDR-adjusted p = .022). URVs in GABAergic pathway genes were associated with familial GGE (OR = 1.8, 95% CI = 1.3-2.5, FDR-adjusted p = .0024) but not with sporadic GGE (OR = 1.3, 95% CI = .9-1.9, FDR-adjusted p = .19).

SIGNIFICANCE

URVs in GABRG2 are likely an important risk factor for familial GGE. The association of gene sets of GABAergic signaling with familial GGE is more prominent than with sporadic GGE.

A de novo missense variant in GABRA4 alters receptor function in an epileptic and neurodevelopmental phenotype

Variants in γ‐aminobutyric acid A (GABA_A) receptor genes cause different forms of epilepsy and neurodevelopmental disorders. To date, GABRA4, encoding the α4‐subunit, has not been associated with a monogenic condition. However, preclinical evidence points toward seizure susceptibility. Here, we report a de novo missense variant in GABRA4 (c.899C>T, p.Thr300Ile) in an individual with early‐onset drug‐resistant epilepsy and neurodevelopmental abnormalities. An electrophysiological characterization of the variant, which is located in the pore‐forming domain, shows accelerated desensitization and a lack of seizure‐protective neurosteroid function. In conclusion, our findings strongly suggest an association between de novo variation in GABRA4 and a neurodevelopmental disorder with epilepsy.

Whole-Exome Sequencing Implicates the USP34 rs777591A > G Intron Variant in Chronic Obstructive Pulmonary Disease in a Kashi Cohort

Genetic factors are important factors in chronic obstructive pulmonary disease (COPD) onset. Plenty of risk and new causative genes for COPD have been identified in patients of the Chinese Han population. In contrast, we know considerably little concerning the genetics in the Kashi COPD population (Uyghur). This study aims at clarifying the genetic maps regarding COPD susceptibility in Kashi (China). Whole-exome sequencing (WES) was used to analyze three Uyghur families with COPD in Kashi (eight patients and one healthy control). Sanger sequencing was also used to verify the WES results in 541 unrelated Uyghur COPD patients and 534 Uyghur healthy controls. WES showed 72 single nucleotide variants (SNVs), two deletions, and small insertions (InDels), 26 copy number variants (CNVs), and 34 structural variants (SVs), including g.71230620T > A (rs12449210T > A, NC_000,016.10) in the HYDIN axonemal central pair apparatus protein (HYDIN) gene and g.61190482A > G (rs777591A > G, NC_000002.12) in the ubiquitin-specific protease 34 (USP34) gene. After Sanger sequencing, we found that rs777591“AA” under different genetic models except for the dominant model (adjusted OR = 0.8559, 95%CI 0.6568–1.115, p > .05), could significantly reduce COPD risk, but rs12449210T > A was not related to COPD. In stratified analysis of smoking status, rs777591“AA” reduced COPD risk significantly among the nonsmoker group. Protein and mRNA expression of USP34 in cigarette smoke extract-treated BEAS-2b cells increased significantly compared with those in the control group. Our findings associate the USP34 rs777591“AA” genotype as a protector factor in COPD.

Genetic generalized epilepsies in adults - challenging assumptions and dogmas

Genetic generalized epilepsy (GGE) syndromes start during childhood or adolescence, and four commonly persist into adulthood, making up 15-20% of all cases of epilepsy in adults. These four GGE syndromes are childhood absence epilepsy, juvenile absence epilepsy, juvenile myoclonic epilepsy and epilepsy with generalized tonic-clonic seizures alone. However, in ~20% of patients with GGE, characteristics of more than one syndrome are present. Novel insights into the genetic aetiology, comorbidities and prognosis of the GGE syndromes have emerged and challenge traditional concepts about these conditions. Evidence has shown that the mode of inheritance in GGE is mostly polygenic. Neuropsychological and imaging studies indicate similar abnormalities in unaffected relatives of patients with GGE, supporting the concept that underlying alterations in bilateral frontothalamocortical networks are genetically determined. Contrary to popular belief, first-line anti-seizure medication often fails to provide seizure freedom in combination with good tolerability. Nevertheless, long-term follow-up studies have shown that with advancing age, many patients can discontinue their anti-seizure medication without seizure relapses. Several outcome predictors have been identified, but prognosis across the syndromes is more homogeneous than previously assumed. Overall, overlap in pathophysiology, seizure types, treatment responses and outcomes support the idea that GGEs are not separate nosological entities but represent a neurobiological continuum.

Whole-exome sequencing with targeted analysis and epilepsy after acute symptomatic neonatal seizures

BACKGROUND

The contribution of pathogenic gene variants with development of epilepsy after acute symptomatic neonatal seizures is not known.

METHODS

Case-control study of 20 trios in children with a history of acute symptomatic neonatal seizures: 10 with and 10 without post-neonatal epilepsy. We performed whole-exome sequencing (WES) and identified pathogenic de novo, transmitted, and non-transmitted variants from established and candidate epilepsy association genes and correlated prevalence of these variants with epilepsy outcomes. We performed a sensitivity analysis with genes associated with coronary artery disease (CAD). We analyzed variants throughout the exome to evaluate for differential enrichment of functional properties using exploratory KEGG searches.

RESULTS

Querying 200 established and candidate epilepsy genes, pathogenic variants were identified in 5 children with post-neonatal epilepsy yet in only 1 child without subsequent epilepsy. There was no difference in the number of trios with non-transmitted pathogenic variants in epilepsy or CAD genes. An exploratory KEGG analysis demonstrated a relative enrichment in cell death pathways in children without subsequent epilepsy.

CONCLUSIONS

In this pilot study, children with epilepsy after acute symptomatic neonatal seizures had a higher prevalence of coding variants with a targeted epilepsy gene sequencing analysis compared to those patients without subsequent epilepsy.

IMPACT

We performed whole-exome sequencing (WES) in 20 trios, including 10 children with epilepsy and 10 without epilepsy, both after acute symptomatic neonatal seizures. Children with post-neonatal epilepsy had a higher burden of pathogenic variants in epilepsy-associated genes compared to those without post-neonatal epilepsy. Future studies evaluating this association may lead to a better understanding of the risk of epilepsy after acute symptomatic neonatal seizures and elucidate molecular pathways that are dysregulated after brain injury and implicated in epileptogenesis.

OUP accepted manuscript

Epilepsy is one of the most frequent neurological diseases, with focal epilepsy accounting for the largest number of cases. The genetic alterations involved in focal epilepsy are far from being fully elucidated.

Here, we show that defective lipid signalling caused by heterozygous ultra-rare variants in PIK3C2B, encoding for the class II phosphatidylinositol 3-kinase PI3K-C2β, underlie focal epilepsy in humans. We demonstrate that patients’ variants act as loss-of-function alleles, leading to impaired synthesis of the rare signalling lipid phosphatidylinositol 3,4-bisphosphate, resulting in mTORC1 hyperactivation. In vivo, mutant Pik3c2b alleles caused dose-dependent neuronal hyperexcitability and increased seizure susceptibility, indicating haploinsufficiency as a key driver of disease. Moreover, acute mTORC1 inhibition in mutant mice prevented experimentally induced seizures, providing a potential therapeutic option for a selective group of patients with focal epilepsy.

Our findings reveal an unexpected role for class II PI3K-mediated lipid signalling in regulating mTORC1-dependent neuronal excitability in mice and humans.

Genome-Wide Association Study-Guided Exome Rare Variant Burden Analysis Identifies IL1R1 and CD3E as Potential Autoimmunity Risk Genes for Celiac Disease

Celiac disease (CeD) is a multifactorial autoimmune enteropathy characterized by the overactivation of the immune system in response to dietary gluten. The molecular etiology of CeD is still not well-understood. Therefore, this study aims to identify potential candidate genes involved in CeD pathogenesis by applying multilayered system biology approaches. Initially, we identified rare coding variants shared between the affected siblings in two rare Arab CeD families by whole-exome sequencing (WES). Then we used the STRING database to construct a protein network of rare variants and genome-wide association study (GWAS) loci to explore their molecular interactions in CeD. Furthermore, the hub genes identified based on network topology parameters were subjected to a series of computational validation analyses like pathway enrichment, gene expression, knockout mouse model, and variant pathogenicity predictions. Our findings have shown the absence of rare variants showing classical Mendelian inheritance in both families. However, interactome analysis of rare WES variants and GWAS loci has identified a total of 11 hub genes. The multidimensional computational analysis of hub genes has prioritized IL1R1 for family A and CD3E for family B as potential genes. These genes were connected to CeD pathogenesis pathways of T-cell selection, cytokine signaling, and adaptive immune response. Future multi-omics studies may uncover the roles of IL1R1 and CD3E in gluten sensitivity. The present investigation lays forth a novel approach integrating next-generation sequencing (NGS) of familial cases, GWAS, and computational analysis for solving the complex genetic architecture of CeD.

Genome sequencing of 320 Chinese children with epilepsy: a clinical and molecular study

The aim of this study is to evaluate the diagnostic value of genome sequencing in children with epilepsy, and to provide genome sequencing-based insights into the molecular genetic mechanisms of epilepsy to help establish accurate diagnoses, design appropriate treatments and assist in genetic counselling. We performed genome sequencing on 320 Chinese children with epilepsy, and interpreted single-nucleotide variants and copy number variants of all samples. The complete pedigree and clinical data of the probands were established and followed up. The clinical phenotypes, treatments, prognoses and genotypes of the patients were analysed. Age at seizure onset ranged from 1 day to 17 years, with a median of 4.3 years. Pathogenic/likely pathogenic variants were found in 117 of the 320 children (36.6%), of whom 93 (29.1%) had single-nucleotide variants, 22 (6.9%) had copy number variants and two had both single-nucleotide variants and copy number variants. Single-nucleotide variants were most frequently found in SCN1A (10/95, 10.5%), which is associated with Dravet syndrome, followed by PRRT2 (8/95, 8.4%), which is associated with benign familial infantile epilepsy, and TSC2 (7/95, 7.4%), which is associated with tuberous sclerosis. Among the copy number variants, there were three with a length <25 kilobases. The most common recurrent copy number variants were 17p13.3 deletions (5/24, 20.8%), 16p11.2 deletions (4/24, 16.7%), and 7q11.23 duplications (2/24, 8.3%), which are associated with epilepsy, developmental retardation and congenital abnormalities. Four particular 16p11.2 deletions and two 15q11.2 deletions were considered to be susceptibility factors contributing to neurodevelopmental disorders associated with epilepsy. The diagnostic yield was 75.0% in patients with seizure onset during the first postnatal month, and gradually decreased in patients with seizure onset at a later age. Forty-two patients (13.1%) were found to be specifically treatable for the underlying genetic cause identified by genome sequencing. Three of them received corresponding targeted therapies and demonstrated favourable prognoses. Genome sequencing provides complete genetic diagnosis, thus enabling individualized treatment and genetic counselling for the parents of the patients. Genome sequencing is expected to become the first choice of methods for genetic testing of patients with epilepsy.

GRIN2A Variants Associated With Idiopathic Generalized Epilepsies

Objective: The objective of this study is to explore the role of GRIN2A gene in idiopathic generalized epilepsies and the potential underlying mechanism for phenotypic variation.

Methods: Whole-exome sequencing was performed in a cohort of 88 patients with idiopathic generalized epilepsies. Electro-physiological alterations of the recombinant N-methyl-D-aspartate receptors (NMDARs) containing GluN2A mutants were examined using two-electrode voltage-clamp recordings. The alterations of protein expression were detected by immunofluorescence staining and biotinylation. Previous studies reported that epilepsy related GRIN2A missense mutations were reviewed. The correlation among phenotypes, functional alterations, and molecular locations was analyzed.

Results: Three novel heterozygous missense GRIN2A mutations (c.1770A > C/p.K590N, c.2636A > G/p.K879R, and c.3199C > T/p.R1067W) were identified in three unrelated cases. Electrophysiological analysis demonstrated R1067W significantly increased the current density of GluN1/GluN2A NMDARs. Immunofluorescence staining indicated GluN2A mutants had abundant distribution in the membrane and cytoplasm. Western blotting showed the ratios of surface and total expression of the three GluN2A-mutants were significantly increased comparing to the wild type. Further analysis on the reported missense mutations demonstrated that mutations with severe gain-of-function were associated with epileptic encephalopathy, while mutations with mild gain of function were associated with mild phenotypes, suggesting a quantitative correlation between gain-of-function and phenotypic severity. The mutations located around transmembrane domains were more frequently associated with severe phenotypes and absence seizure-related mutations were mostly located in carboxyl-terminal domain, suggesting molecular sub-regional effects.

Significance: This study revealed GRIN2A gene was potentially a candidate pathogenic gene of idiopathic generalized epilepsies. The functional quantitative correlation and the molecular sub-regional implication of mutations helped in explaining the relatively mild clinical phenotypes and incomplete penetrance associated with GRIN2A variants.

Ahmed E, Baldassari S, Achaz G, Elmugadam FA, Abdelgadir WA, Baulac S, Buratti J, Abdalla O, Gamil S, Alzubeir M, Abubaker R, Noé E, Elsayed L, Ahmed AE, Leguern E. Involvement of ADGRV1 Gene in Familial Forms of Genetic Generalized Epilepsy

Background: Genetic generalized epilepsies (GGE) including childhood absence epilepsy (CAE), juvenile absence epilepsy (JAE), juvenile myoclonic epilepsy (JME), and GGE with tonic-clonic seizures alone (GGE-TCS), are common types of epilepsy mostly determined by a polygenic mode of inheritance. Recent studies showed that susceptibility genes for GGE are numerous, and their variants rare, challenging their identification. In this study, we aimed to assess GGE genetic etiology in a Sudanese population. Methods: We performed whole-exome sequencing (WES) on DNA of 40 patients from 20 Sudanese families with GGE searching for candidate susceptibility variants, which were prioritized by CADD software and functional features of the corresponding gene. We assessed their segregation in 138 individuals and performed genotype-phenotype correlations. Results: In a family including three sibs with GGE-TCS, we identified a rare missense variant in ADGRV1 encoding an adhesion G protein-coupled receptor V1, which was already involved in the autosomal recessive Usher type C syndrome. In addition, five other ADGRV1 rare missense variants were identified in four additional families and absent from 119 Sudanese controls. In one of these families, an ADGRV1 variant was found at a homozygous state, in a female more severely affected than her heterozygous brother, suggesting a gene dosage effect. In the five families, GGE phenotype was statistically associated with ADGRV1 variants (0R = 0.9 103). Conclusion: This study highly supports, for the first time, the involvement of ADGRV1 missense variants in familial GGE and that ADGRV1 is a susceptibility gene for CAE/JAE and GGE-TCS phenotypes.

Epilepsy in the mTORopathies: opportunities for precision medicine

The mechanistic target of rapamycin signalling pathway serves as a ubiquitous regulator of cell metabolism, growth, proliferation and survival. The main cellular activity of the mechanistic target of rapamycin cascade funnels through mechanistic target of rapamycin complex 1, which is inhibited by rapamycin, a macrolide compound produced by the bacterium Streptomyces hygroscopicus. Pathogenic variants in genes encoding upstream regulators of mechanistic target of rapamycin complex 1 cause epilepsies and neurodevelopmental disorders. Tuberous sclerosis complex is a multisystem disorder caused by mutations in mechanistic target of rapamycin regulators TSC1 or TSC2, with prominent neurological manifestations including epilepsy, focal cortical dysplasia and neuropsychiatric disorders. Focal cortical dysplasia type II results from somatic brain mutations in mechanistic target of rapamycin pathway activators MTOR, AKT3, PIK3CA and RHEB and is a major cause of drug-resistant epilepsy. DEPDC5, NPRL2 and NPRL3 code for subunits of the GTPase-activating protein (GAP) activity towards Rags 1 complex (GATOR1), the principal amino acid-sensing regulator of mechanistic target of rapamycin complex 1. Germline pathogenic variants in GATOR1 genes cause non-lesional focal epilepsies and epilepsies associated with malformations of cortical development. Collectively, the mTORopathies are characterized by excessive mechanistic target of rapamycin pathway activation and drug-resistant epilepsy. In the first large-scale precision medicine trial in a genetically mediated epilepsy, everolimus (a synthetic analogue of rapamycin) was effective at reducing seizure frequency in people with tuberous sclerosis complex. Rapamycin reduced seizures in rodent models of DEPDC5-related epilepsy and focal cortical dysplasia type II. This review outlines a personalized medicine approach to the management of epilepsies in the mTORopathies. We advocate for early diagnostic sequencing of mechanistic target of rapamycin pathway genes in drug-resistant epilepsy, as identification of a pathogenic variant may point to an occult dysplasia in apparently non-lesional epilepsy or may uncover important prognostic information including, an increased risk of sudden unexpected death in epilepsy in the GATORopathies or favourable epilepsy surgery outcomes in focal cortical dysplasia type II due to somatic brain mutations. Lastly, we discuss the potential therapeutic application of mechanistic target of rapamycin inhibitors for drug-resistant seizures in GATOR1-related epilepsies and focal cortical dysplasia type II.

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

Background

Analyses of few gene-sets in epilepsy showed a potential to unravel key disease associations. We set out to investigate the burden of ultra-rare variants (URVs) in a comprehensive range of biologically informed gene-sets presumed to be implicated in epileptogenesis.

Methods

The burden of 12 URV types in 92 gene-sets was compared between cases and controls using whole exome sequencing data from individuals of European descent with developmental and epileptic encephalopathies (DEE, n = 1,003), genetic generalized epilepsy (GGE, n = 3,064), or non-acquired focal epilepsy (NAFE, n = 3,522), collected by the Epi25 Collaborative, compared to 3,962 ancestry-matched controls.

Findings

Missense URVs in highly constrained regions were enriched in neuron-specific and developmental genes, whereas genes not expressed in brain were not affected. GGE featured a higher burden in gene-sets derived from inhibitory vs. excitatory neurons or associated receptors, whereas the opposite was found for NAFE, and DEE featured a burden in both. Top-ranked susceptibility genes from recent genome-wide association studies (GWAS) and gene-sets derived from generalized vs. focal epilepsies revealed specific enrichment patterns of URVs in GGE vs. NAFE.

Interpretation

Missense URVs affecting highly constrained sites differentially impact genes expressed in inhibitory vs. excitatory pathways in generalized vs. focal epilepsies. The excess of URVs in top-ranked GWAS risk-genes suggests a convergence of rare deleterious and common risk-variants in the pathogenesis of generalized and focal epilepsies.

Funding

DFG Research Unit FOR-2715 (Germany), FNR (Luxembourg), NHGRI (US), NHLBI (US), DAAD (Germany).

Genomic analysis of "microphenotypes" in epilepsy

Large international consortia examining the genomic architecture of the epilepsies focus on large diagnostic subgroupings such as "all focal epilepsy" and "all genetic generalized epilepsy". In addition, phenotypic data are generally entered into these large discovery databases in a unidirectional manner at one point in time only. However, there are many smaller phenotypic subgroupings in epilepsy, many of which may have unique genomic risk factors. Such a subgrouping or "microphenotype" may be defined as an uncommon or rare phenotype that is well recognized by epileptologists and the epilepsy community, and which may or may not be formally recognized within the International League Against Epilepsy classification system. Here we examine the genetic structure of a number of such microphenotypes and report in particular on two interesting clinical phenotypes, Jeavons syndrome and pediatric status epilepticus. Although no single gene reached exome-wide statistical significance to be associated with any of the diagnostic categories, we observe enrichment of rare damaging variants in established epilepsy genes among Landau-Kleffner patients (GRIN2A) and pediatric status epilepticus patients (MECP2, SCN1A, SCN2A, SCN8A).

Chromosomal microarray and exome sequencing in unexplained early infantile epileptic encephalopathies in a highly consanguineous population

AIM

To identify genetic causes for early infantile epileptic encephalopathies (EIEE) in Turkish children with mostly consanguineous parents.

METHODS

In a selected EIEE group (N = 59) based on results of nongenetic and initial genetic testing with unexplained etiology, 49 patients underwent array-based comparative genomic hybridization (aCGH) and 49 patients underwent whole exome sequencing (WES) including 39 with negative aCGH results and 10 with WES-only.

RESULTS

Diagnostic yield of aCGH and WES for pathogenic or likely pathogenic variants was 14.3% and 38.8%, respectively. Including de novo variants of uncertain significance linked to compatible phenotypes, increased the diagnostic yield of WES to 61.2%. Out of 38 positive variants, 18 (47.4%) were novel and 16 (42.1%) were de novo. Twenty-one (56.8%) patients had recessive variants inherited from mostly consanguineous healthy parents (85.7%). Fourteen (37.8%) of patients with diagnostic results had positive variants in established EIEE genes. Seizures started during neonatal period in 32.4% patients. Posture or movement disorders were comorbid with EIEE in 40.5% of diagnosed patients. We identified treatable metabolic disorders in 8.1% of patients and pathogenic variants in genes which support using targeted medicine in 19% of patients.

CONCLUSIONS

Detailed electro-clinical phenotyping led to expansion of some of the known phenotypes with non-neurological and neurological findings in addition to seizures, as well as suggestion of candidate genes (SEC24B, SLC16A2 and PRICKLE2) and a copy number variant (microduplication of Xp21.1p11.4). The high ratio of recessive inheritance could be important for family counseling.

The CB2 Receptor as a Novel Therapeutic Target for Epilepsy Treatment

Epilepsy is characterized by repeated spontaneous bursts of neuronal hyperactivity and high synchronization in the central nervous system. It seriously affects the quality of life of epileptic patients, and nearly 30% of individuals are refractory to treatment of antiseizure drugs. Therefore, there is an urgent need to develop new drugs to manage and control refractory epilepsy. Cannabinoid ligands, including selective cannabinoid receptor subtype (CB₁ or CB₂ receptor) ligands and non-selective cannabinoid (synthetic and endogenous) ligands, may serve as novel candidates for this need. Cannabinoid appears to regulate seizure activity in the brain through the activation of CB₁ and CB₂ cannabinoid receptors (CB₁R and CB₂R). An abundant series of cannabinoid analogues have been tested in various animal models, including the rat pilocarpine model of acquired epilepsy, a pentylenetetrazol model of myoclonic seizures in mice, and a penicillin-induced model of epileptiform activity in the rats. The accumulating lines of evidence show that cannabinoid ligands exhibit significant benefits to control seizure activity in different epileptic models. In this review, we summarize the relationship between brain CB₂ receptors and seizures and emphasize the potential mechanisms of their therapeutic effects involving the influences of neurons, astrocytes, and microglia cells. The unique features of CB₂Rs, such as lower expression levels under physiological conditions and high inducibility under epileptic conditions, make it an important target for future research on drug-resistant epilepsy.

Recapitulation of the EEF1A2 D252H neurodevelopmental disorder-causing missense mutation in mice reveals a toxic gain of function

Heterozygous de novo mutations in EEF1A2, encoding the tissue-specific translation elongation factor eEF1A2, have been shown to cause neurodevelopmental disorders including often severe epilepsy and intellectual disability. The mutational profile is unusual; ~50 different missense mutations have been identified but no obvious loss of function mutations, though large heterozygous deletions are known to be compatible with life. A key question is whether the heterozygous missense mutations operate through haploinsufficiency or a gain of function mechanism, an important prerequisite for design of therapeutic strategies. In order both to address this question and to provide a novel model for neurodevelopmental disorders resulting from mutations in EEF1A2, we created a new mouse model of the D252H mutation. This mutation causes the eEF1A2 protein to be expressed at lower levels in brain but higher in muscle in the mice. We compared both heterozygous and homozygous D252H and null mutant mice using behavioural and motor phenotyping alongside molecular modelling and analysis of binding partners. Although the proteomic analysis pointed to a loss of function for the D252H mutant protein, the D252H homozygous mice were more severely affected than null homozygotes on the same genetic background. Mice that are heterozygous for the missense mutation show no behavioural abnormalities but do have sex-specific deficits in body mass and motor function. The phenotyping of our novel mouse lines, together with analysis of molecular modelling and interacting proteins, suggest that the D252H mutation results in a gain of function.

Paradoxical hyperexcitability from NaV1.2 sodium channel loss in neocortical pyramidal cells

Loss-of-function variants in the gene SCN2A, which encodes the sodium channel NaV1.2, are strongly associated with autism spectrum disorder and intellectual disability. An estimated 20%-30% of children with these variants also suffer from epilepsy, with altered neuronal activity originating in neocortex, a region where NaV1.2 channels are expressed predominantly in excitatory pyramidal cells. This is paradoxical, as sodium channel loss in excitatory cells would be expected to dampen neocortical activity rather than promote seizure. Here, we examined pyramidal neurons lacking NaV1.2 channels and found that they were intrinsically hyperexcitable, firing high-frequency bursts of action potentials (APs) despite decrements in AP size and speed. Compartmental modeling and dynamic-clamp recordings revealed that NaV1.2 loss prevented potassium channels from properly repolarizing neurons between APs, increasing overall excitability by allowing neurons to reach threshold for subsequent APs more rapidly. This cell-intrinsic mechanism may, therefore, account for why SCN2A loss-of-function can paradoxically promote seizure.

Recent ultra-rare inherited variants implicate new autism candidate risk genes

Autism is a highly heritable complex disorder in which de novo mutation (DNM) variation contributes significantly to risk. Using whole-genome sequencing data from 3,474 families, we investigate another source of large-effect risk variation, ultra-rare variants. We report and replicate a transmission disequilibrium of private, likely gene-disruptive (LGD) variants in probands but find that 95% of this burden resides outside of known DNM-enriched genes. This variant class more strongly affects multiplex family probands and supports a multi-hit model for autism. Candidate genes with private LGD variants preferentially transmitted to probands converge on the E3 ubiquitin-protein ligase complex, intracellular transport and Erb signaling protein networks. We estimate that these variants are approximately 2.5 generations old and significantly younger than other variants of similar type and frequency in siblings. Overall, private LGD variants are under strong purifying selection and appear to act on a distinct set of genes not yet associated with autism.

MTR3D: identifying regions within protein tertiary structures under purifying selection

The identification of disease-causal variants is non-trivial. By mapping population variation from over 448,000 exome and genome sequences to over 81,000 experimental structures and homology models of the human proteome, we have calculated both regional intolerance to missense variation (Missense Tolerance Ratio, MTR), using a sliding window of 21–41 codons, and introduce a new 3D spatial intolerance to missense variation score (3D Missense Tolerance Ratio, MTR3D), using spheres of 5–8 Å. We show that the MTR3D is less biased by regions with limited data and more accurately identifies regions under purifying selection than estimates relying on the sequence alone. Intolerant regions were highly enriched for both ClinVar pathogenic and COSMIC somatic missense variants (Mann–Whitney U test P < 2.2 × 10⁻¹⁶). Further, we combine sequence- and spatial-based scores to generate a consensus score, MTRX, which distinguishes pathogenic from benign variants more accurately than either score separately (AUC = 0.85). The MTR3D server enables easy visualisation of population variation, MTR, MTR3D and MTRX scores across the entire gene and protein structure for >17,000 human genes and >42,000 alternative alternate transcripts, including both Ensembl and RefSeq transcripts. MTR3D is freely available by user-friendly web-interface and API at http://biosig.unimelb.edu.au/mtr3d/.

Increased burden of rare variants in genes of the endosomal Toll-like receptor pathway in patients with systemic lupus erythematosus

OBJECTIVE

To compare the frequency of rare variants in genes of the pathophysiologically relevant endosomal Toll-like receptor (eTLR) pathway and any quantifiable differences in variant rarity, predicted deleteriousness, or molecular proximity in patients with systemic lupus erythematosus (SLE) and healthy controls.

PATIENTS AND METHODS

65 genes associated with the eTLR pathway were identified by literature search and pathway analysis. Using next generation sequencing techniques, these were compared in two randomised cohorts of patients with SLE (n = 114 and n = 113) with 197 healthy controls. Genetically determined ethnicity was used to normalise minor allele frequencies (MAF) for the identified genetic variants and these were then compared by their frequency: rare (MAF < 0.005), uncommon (MAF 0.005-0.02), and common (MAF >0.02). This was compared to the results for 65 randomly selected genes.

RESULTS

Patients with SLE are more likely to carry a rare nonsynonymous variant affecting proteins within the eTLR pathway than healthy controls. Furthermore, individuals with SLE are more likely to have multiple rare variants in this pathway. There were no differences in rarity, Combined Annotation Dependent Depletion (CADD) score, or molecular proximity for rare eTLR pathway variants.

CONCLUSIONS

Rare non-synonymous variants are enriched in patients with SLE in the eTLR pathway. This supports the hypothesis that SLE arises from several rare variants of relatively large effect rather than many common variants of small effect.

Assessing the role of rare genetic variants in drug-resistant, non-lesional focal epilepsy

OBJECTIVE

Resistance to antiseizure medications (ASMs) is one of the major concerns in the treatment of epilepsy. Despite the increasing number of ASMs available, the proportion of individuals with drug-resistant epilepsy remains unchanged. In this study, we aimed to investigate the role of rare genetic variants in ASM resistance.

METHODS

We performed exome sequencing of 1,128 individuals with non-familial non-acquired focal epilepsy (NAFE) (762 non-responders, 366 responders) and were provided with 1,734 healthy controls. We undertook replication in a cohort of 350 individuals with NAFE (165 non-responders, 185 responders). We performed gene-based and gene-set-based kernel association tests to investigate potential enrichment of rare variants in relation to drug response status and to risk for NAFE.

RESULTS

We found no gene or gene set that reached genome-wide significance. Yet, we identified several prospective candidate genes - among them DEPDC5, which showed a potential association with resistance to ASMs. We found some evidence for an enrichment of truncating variants in dominant familial NAFE genes in our cohort of non-familial NAFE and in association with drug-resistant NAFE.

INTERPRETATION

Our study identifies potential candidate genes for ASM resistance. Our results corroborate the role of rare variants for non-familial NAFE and imply their involvement in drug-resistant epilepsy. Future large-scale genetic research studies are needed to substantiate these findings.

Haploinsufficiency, Dominant Negative, and Gain-of-Function Mechanisms in Epilepsy: Matching Therapeutic Approach to the Pathophysiology

This review summarizes the pathogenic mechanisms that underpin the monogenic epilepsies and discusses the potential of novel precision therapeutics to treat these disorders. Pathogenic mechanisms of epilepsy include recessive (null alleles), haploinsufficiency, imprinting, gain-of-function, and dominant negative effects. Understanding which pathogenic mechanism(s) that underlie each genetic epilepsy is pivotal to design precision therapies that are most likely to be beneficial for the patient. Novel therapeutics discussed include gene therapy, gene editing, antisense oligonucleotides, and protein replacement. Discussions are illustrated and reinforced with examples from the literature.

Screening Platforms for Genetic Epilepsies-Zebrafish, iPSC-Derived Neurons, and Organoids

Recent advances in molecular and cellular engineering, such as human cell reprogramming, genome editing, and patient-specific organoids, have provided unprecedented opportunities for investigating human disorders in both animals and human-based models at an improved pace and precision. This progress will inevitably lead to the development of innovative drug-screening platforms and new patient-specific therapeutics. In this review, we discuss recent advances that have been made using zebrafish and human-induced pluripotent stem cell (iPSC)-derived neurons and organoids for modeling genetic epilepsies. We also provide our prospective on how these models can potentially be combined to build new screening platforms for antiseizure and antiepileptogenic drug discovery that harness the robustness and tractability of zebrafish models as well as the patient-specific genetics and biology of iPSC-derived neurons and organoids.

Unraveling synonymous and deep intronic variants causing aberrant splicing in two genetically undiagnosed epilepsy families

BACKGROUND

Variants identified through parent-child trio-WES yield up to 28-55% positive diagnostic rate across a variety of Mendelian disorders, there remain numerous patients who do not receive a genetic diagnosis. Studies showed that some aberrant splicing variants, which are either not readily detectable by WES or could be miss-interpreted by regular detecting pipelines, are highly relevant to human diseases.

METHODS

We retrospectively investigated the negative molecular diagnostics through trio-WES for 15 genetically undiagnosed patients whose clinical manifestations were highly suspected to be genetic disorders with well-established genotype-phenotype relationships. We scrutinized the synonymous variants from WES data and Sanger sequenced the suspected intronic region for deep intronic variants. The functional consequences of variants were analyzed by in vitro minigene experiments.

RESULTS

Here, we report two abnormal splicing events, one of which caused exon truncating due to the activation of cryptic splicing site by a synonymous variant; the other caused partial intron retention due to the generation of splicing sites by a deep intronic variant.

CONCLUSIONS

We suggest that, despite initial negative genetic test results in clinically highly suspected genetic diseases, the combination of predictive bioinformatics and functional analysis should be considered to unveil the genetic etiology of undiagnosed rare diseases.

Gabra2 is a genetic modifier of Dravet syndrome in mice

Pathogenic variants in epilepsy genes result in a spectrum of clinical severity. One source of phenotypic heterogeneity is modifier genes that affect expressivity of a primary pathogenic variant. Mouse epilepsy models also display varying degrees of clinical severity on different genetic backgrounds. Mice with heterozygous deletion of Scn1a (Scn1a^+/−) model Dravet syndrome, a severe epilepsy most often caused by SCN1A haploinsufficiency. Scn1a^+/− mice recapitulate features of Dravet syndrome, including spontaneous seizures, sudden death, and cognitive/behavioral deficits. Scn1a^+/− mice maintained on the 129S6/SvEvTac (129) strain have normal lifespan and no spontaneous seizures. In contrast, admixture with C57BL/6J (B6) results in epilepsy and premature lethality. We previously mapped Dravet Survival Modifier loci (Dsm1-Dsm5) responsible for strain-dependent differences in survival. Gabra2, encoding the GABA_A α2 subunit, was nominated as a candidate modifier at Dsm1. Direct measurement of GABA_A receptors found lower abundance of α2-containing receptors in hippocampal synapses of B6 mice relative to 129. We also identified a B6-specific single nucleotide deletion within Gabra2 that lowers mRNA and protein by nearly 50%. Repair of this deletion reestablished normal levels of Gabra2 expression. In this study, we used B6 mice with a repaired Gabra2 allele to evaluate Gabra2 as a genetic modifier of severity in Scn1a^+/− mice. Gabra2 repair restored transcript and protein expression, increased abundance of α2-containing GABA_A receptors in hippocampal synapses, and rescued epilepsy phenotypes of Scn1a^+/− mice. These findings validate Gabra2 as a genetic modifier of Dravet syndrome, and support enhancing function of α₂-containing GABA_A receptors as treatment strategy for Dravet syndrome.

Sub-genic intolerance, ClinVar, and the epilepsies: A whole-exome sequencing study of 29,165 individuals

Both mild and severe epilepsies are influenced by variants in the same genes, yet an explanation for the resulting phenotypic variation is unknown. As part of the ongoing Epi25 Collaboration, we performed a whole-exome sequencing analysis of 13,487 epilepsy-affected individuals and 15,678 control individuals. While prior Epi25 studies focused on gene-based collapsing analyses, we asked how the pattern of variation within genes differs by epilepsy type. Specifically, we compared the genetic architectures of severe developmental and epileptic encephalopathies (DEEs) and two generally less severe epilepsies, genetic generalized epilepsy and non-acquired focal epilepsy (NAFE). Our gene-based rare variant collapsing analysis used geographic ancestry-based clustering that included broader ancestries than previously possible and revealed novel associations. Using the missense intolerance ratio (MTR), we found that variants in DEE-affected individuals are in significantly more intolerant genic sub-regions than those in NAFE-affected individuals. Only previously reported pathogenic variants absent in available genomic datasets showed a significant burden in epilepsy-affected individuals compared with control individuals, and the ultra-rare pathogenic variants associated with DEE were located in more intolerant genic sub-regions than variants associated with non-DEE epilepsies. MTR filtering improved the yield of ultra-rare pathogenic variants in affected individuals compared with control individuals. Finally, analysis of variants in genes without a disease association revealed a significant burden of loss-of-function variants in the genes most intolerant to such variation, indicating additional epilepsy-risk genes yet to be discovered. Taken together, our study suggests that genic and sub-genic intolerance are critical characteristics for interpreting the effects of variation in genes that influence epilepsy.

Phenotypic homogeneity in childhood epilepsies evolves in gene-specific patterns across 3251 patient-years of clinical data

While genetic studies of epilepsies can be performed in thousands of individuals, phenotyping remains a manual, non-scalable task. A particular challenge is capturing the evolution of complex phenotypes with age. Here, we present a novel approach, applying phenotypic similarity analysis to a total of 3251 patient-years of longitudinal electronic medical record data from a previously reported cohort of 658 individuals with genetic epilepsies. After mapping clinical data to the Human Phenotype Ontology, we determined the phenotypic similarity of individuals sharing each genetic etiology within each 3-month age interval from birth up to a maximum age of 25 years. 140 of 600 (23%) of all 27 genes and 3-month age intervals with sufficient data for calculation of phenotypic similarity were significantly higher than expect by chance. 11 of 27 genetic etiologies had significant overall phenotypic similarity trajectories. These do not simply reflect strong statistical associations with single phenotypic features but appear to emerge from complex clinical constellations of features that may not be strongly associated individually. As an attempt to reconstruct the cognitive framework of syndrome recognition in clinical practice, longitudinal phenotypic similarity analysis extends the traditional phenotyping approach by utilizing data from electronic medical records at a scale that is far beyond the capabilities of manual phenotyping. Delineation of how the phenotypic homogeneity of genetic epilepsies varies with age could improve the phenotypic classification of these disorders, the accuracy of prognostic counseling, and by providing historical control data, the design and interpretation of precision clinical trials in rare diseases.

Loss-of-function variants in Kv 11.1 cardiac channels as a biomarker for SUDEP

Objective

To compare the frequency and impact on the channel function of KCNH2 variants in SUDEP patients with epilepsy controls comprising patients older than 50 years, a group with low SUDEP risk, and establish loss‐of‐function KCNH2 variants as predictive biomarkers of SUDEP risk.

Methods

We searched for KCNH2 variants with a minor allele frequency of <5%. Functional analysis in Xenopus laevis oocytes was performed for all KCNH2 variants identified.

Results

KCNH2 variants were found in 11.1% (10/90) of SUDEP individuals compared to 6.0% (20/332) of epilepsy controls (p = 0.11). Loss‐of‐function KCNH2 variants, defined as causing >20% reduction in maximal amplitude, were observed in 8.9% (8/90) SUDEP patients compared to 3.3% (11/332) epilepsy controls suggesting about threefold enrichment (nominal p = 0.04). KCNH2 variants that did not change channel function occurred at a similar frequency in SUDEP (2.2%; 2/90) and epilepsy control (2.7%; 9/332) cohorts (p > 0.99). Rare KCNH2 variants (<1% allele frequency) associated with greater loss of function and an ~11‐fold enrichment in the SUDEP cohort (nominal p = 0.03). In silico tools were unable to predict the impact of a variant on function highlighting the need for electrophysiological analysis.

Interpretation

These data show that loss‐of‐function KCNH2 variants are enriched in SUDEP patients when compared to an epilepsy population older than 50 years, suggesting that cardiac mechanisms contribute to SUDEP risk. We propose that genetic screening in combination with functional analysis can identify loss‐of‐function KCNH2 variants that could act as biomarkers of an individual’s SUDEP risk.

Modeling seizures in the Human Phenotype Ontology according to contemporary ILAE concepts makes big phenotypic data tractable

Objective:

The clinical features of epilepsy determine how it is defined, which in turn guides management. Therefore, consideration of the fundamental clinical entities that comprise an epilepsy is essential in the study of causes, trajectories, and treatment responses. The Human Phenotype Ontology (HPO) is used widely in clinical and research genetics for concise communication and modeling of clinical features, allowing extracted data to be harmonized using logical inference. We sought to redesign the HPO seizure subontology to improve its consistency with current epileptological concepts, supporting the use of large clinical data sets in high-throughput clinical and research genomics.

Methods:

We created a new HPO seizure subontology based on the 2017 International League Against Epilepsy (ILAE) Operational Classification of Seizure Types, and integrated concepts of status epilepticus, febrile, reflex, and neonatal seizures at different levels of detail. We compared the HPO seizure subontology prior to, and following, our revision, according to the information that could be inferred about the seizures of 791 individuals from three independent cohorts: 2 previously published and 150 newly recruited individuals. Each cohort’s data were provided in a different format and harmonized using the two versions of the HPO.

Results:

The new seizure subontology increased the number of descriptive concepts for seizures 5-fold. The number of seizure descriptors that could be annotated to the cohort increased by 40% and the total amount of information about individuals’ seizures increased by 38%. The most important qualitative difference was the relationship of focal to bilateral tonic-clonic seizure to generalized-onset and focal-onset seizures.

Investigation of long interspersed element-1 retrotransposons as potential risk factors for idiopathic temporal lobe epilepsy

OBJECTIVE

To determine if long interspersed element-1 (L1) retrotransposons convey risk for idiopathic temporal lobe epilepsy (TLE).

METHODS

Surgically resected temporal cortex from individuals with TLE (N = 33) and postmortem temporal cortex from individuals with no known neurological disease (N = 33) were analyzed for L1 content by Restriction Enzyme Based Enriched L1Hs sequencing (REBELseq). Expression of three KCNIP4 splice variants was assessed by droplet digital PCR (ddPCR). Protein ANalysis THrough Evolutionary Relationships (PANTHER) was used to determine ontologies and pathways for lists of genes harboring L1 insertions.

RESULTS

We identified novel L1 insertions specific to individuals with TLE, and others specific to controls. Although there were no statistically significant differences between cases and controls in the numbers of known and novel L1 insertions, PANTHER analyses of intragenic L1 insertions showed statistically significant enrichments for epilepsy-relevant gene ontologies in both cases and controls. Gene ontologies "neuron projection development" and "calcium ion transmembrane transport" were among those found only in individuals with TLE. We confirmed novel L1 insertions in several genes associated with seizures/epilepsy, including a de novo somatic L1 retrotransposition in KCNIP4 that occurred after neural crest formation in one patient. However, ddPCR results suggest this de novo L1 did not alter KCNIP4 mRNA expression.

SIGNIFICANCE

Given current data from this small cohort, we conclude that L1 elements, either rare heritable germline insertions or de novo somatic retrotranspositions, may contribute only minimally to overall genetic risk for idiopathic TLE. We suggest that further studies in additional patients and additional brain regions are warranted.

Association of SLC32A1 Missense Variants With Genetic Epilepsy With Febrile Seizures Plus

OBJECTIVE

To identify the causative gene in a large unsolved family with genetic epilepsy with febrile seizures plus (GEFS+), we sequenced the genomes of family members, and then determined the contribution of the identified gene to the pathogenicity of epilepsies by examining sequencing data from 2,772 additional patients.

METHODS

We performed whole genome sequencing of 3 members of a GEFS+ family. Subsequently, whole exome sequencing data from 1,165 patients with epilepsy from the Epi4K dataset and 1,329 Australian patients with epilepsy from the Epi25 dataset were interrogated. Targeted resequencing was performed on 278 patients with febrile seizures or GEFS+ phenotypes. Variants were validated and familial segregation examined by Sanger sequencing.

RESULTS

Eight previously unreported missense variants were identified in SLC32A1, coding for the vesicular inhibitory amino acid cotransporter VGAT. Two variants cosegregated with the phenotype in 2 large GEFS+ families containing 8 and 10 affected individuals, respectively. Six further variants were identified in smaller families with GEFS+ or idiopathic generalized epilepsy (IGE).

CONCLUSION

Missense variants in SLC32A1 cause GEFS+ and IGE. These variants are predicted to alter γ-aminobutyric acid (GABA) transport into synaptic vesicles, leading to altered neuronal inhibition. Examination of further epilepsy cohorts will determine the full genotype-phenotype spectrum associated with SLC32A1 variants.

Dravet syndrome-associated mutations in GABRA1, GABRB2 and GABRG2 define the genetic landscape of defects of GABAA receptors

Dravet syndrome is a rare, catastrophic epileptic encephalopathy that begins in the first year of life, usually with febrile or afebrile hemiclonic or generalized tonic-clonic seizures followed by status epilepticus. De novo variants in genes that mediate synaptic transmission such as SCN1A and PCDH19 are often associated with Dravet syndrome. Recently, GABAA receptor subunit genes (GABRs) encoding α1 (GABRA1), β3 (GABRB3) and γ2 (GABRG2), but not β2 (GABRB2) or β1 (GABRB1), subunits are frequently associated with Dravet syndrome or Dravet syndrome-like phenotype. We performed next generation sequencing on 870 patients with Dravet syndrome and identified nine variants in three different GABRs. Interestingly, the variants were all in genes encoding the most common GABAA receptor, the α1β2γ2 receptor. Mutations in GABRA1 (c.644T>C, p. L215P; c.640C>T, p. R214C; c.859G>A; V287I; c.641G>A, p. R214H) and GABRG2 (c.269C>G, p. T90R; c.1025C>T, p. P342L) presented as de novo cases, while in GABRB2 two variants were de novo (c.992T>C, p. F331S; c.542A>T, p. Y181F) and one was autosomal dominant and inherited from the maternal side (c.990_992del, p.330_331del). We characterized the effects of these GABR variants on GABAA receptor biogenesis and channel function. We found that defects in receptor gating were the common deficiency of GABRA1 and GABRB2 Dravet syndrome variants, while mainly trafficking defects were found with the GABRG2 (c.269C>G, p. T90R) variant. It seems that variants in α1 and β2 subunits are less tolerated than in γ2 subunits, since variant α1 and β2 subunits express well but were functionally deficient. This suggests that all of these GABR variants are all targeting GABR genes that encode the assembled α1β2γ2 receptor, and regardless of which of the three subunits are mutated, variants in genes coding for α1, β2 and γ2 receptor subunits make them candidate causative genes in the pathogenesis of Dravet syndrome.

Sodium channelopathies in neurodevelopmental disorders

The voltage-gated sodium channel α-subunit genes comprise a highly conserved gene family. Mutations of three of these genes, SCN1A, SCN2A and SCN8A, are responsible for a significant burden of neurological disease. Recent progress in identification and functional characterization of patient variants is generating new insights and novel approaches to therapy for these devastating disorders. Here we review the basic elements of sodium channel function that are used to characterize patient variants. We summarize a large body of work using global and conditional mouse mutants to characterize the in vivo roles of these channels. We provide an overview of the neurological disorders associated with mutations of the human genes and examples of the effects of patient mutations on channel function. Finally, we highlight therapeutic interventions that are emerging from new insights into mechanisms of sodium channelopathies.

A Transcriptome-Based Drug Discovery Paradigm for Neurodevelopmental Disorders

Advances in genetic discoveries have created substantial opportunities for precision medicine in neurodevelopmental disorders. Many of the genes implicated in these diseases encode proteins that regulate gene expression, such as chromatin-associated proteins, transcription factors, and RNA-binding proteins. The identification of targeted therapeutics for individuals carrying mutations in these genes remains a challenge, as the encoded proteins can theoretically regulate thousands of downstream targets in a considerable number of cell types. Here, we propose the application of a drug discovery approach originally developed for cancer called "transcriptome reversal" for these neurodevelopmental disorders. This approach attempts to identify compounds that reverse gene-expression signatures associated with disease states. ANN NEUROL 2021;89:199-211.

Investigating rare and ultrarare epilepsy syndromes with Drosophila models

One in three epilepsy cases is drug resistant, and seizures often begin in infancy, when they are life-threatening and when therapeutic options are highly limited. An important tool for prioritizing and validating genes associated with epileptic conditions, which is suitable for large-scale screening, is disease modeling in Drosophila. Approximately two-thirds of disease genes are conserved in Drosophila, and gene-specific fly models exhibit behavioral changes that are related to symptoms of epilepsy. Models are based on behavior readouts, seizure-like attacks and paralysis following stimulation, and neuronal, cell-biological readouts that are in the majority based on changes in nerve cell activity or morphology. In this review, we focus on behavioral phenotypes. Importantly, Drosophila modeling is independent of, and complementary to, other approaches that are computational and based on systems analysis. The large number of known epilepsy-associated gene variants indicates a need for efficient research strategies. We will discuss the status quo of epilepsy disease modelling in Drosophila and describe promising steps towards the development of new drugs to reduce seizure rates and alleviate other epileptic symptoms.

Multiplex gene and phenotype network to characterize shared genetic pathways of epilepsy and autism

It is well established that epilepsy and autism spectrum disorder (ASD) commonly co-occur; however, the underlying biological mechanisms of the co-occurence from their genetic susceptibility are not well understood. Our aim in this study is to characterize genetic modules of subgroups of epilepsy and autism genes that have similar phenotypic manifestations and biological functions. We first integrate a large number of expert-compiled and well-established epilepsy- and ASD-associated genes in a multiplex network, where one layer is connected through protein-protein interaction (PPI) and the other layer through gene-phenotype associations. We identify two modules in the multiplex network, which are significantly enriched in genes associated with both epilepsy and autism as well as genes highly expressed in brain tissues. We find that the first module, which represents the Gene Ontology category of ion transmembrane transport, is more epilepsy-focused, while the second module, representing synaptic signaling, is more ASD-focused. However, because of their enrichment in common genes and association with both epilepsy and ASD phenotypes, these modules point to genetic etiologies and biological processes shared between specific subtypes of epilepsy and ASD. Finally, we use our analysis to prioritize new candidate genes for epilepsy (i.e. ANK2, CACNA1E, CACNA2D3, GRIA2, DLG4) for further validation. The analytical approaches in our study can be applied to similar studies in the future to investigate the genetic connections between different human diseases.

Differential DNA Methylation Profiles in Patients with Temporal Lobe Epilepsy and Hippocampal Sclerosis ILAE Type I

Hippocampal sclerosis (HS) is one of the most prevalent pathological types of temporal lobe epilepsy (TLE), and it significantly affects patient prognoses. The methylation of DNA plays an important role in the development of epilepsy. However, few studies have focused on HS subtypes to determine DNA methylation profiles in TLE. This study aimed to determine the pathogenesis of TLE from an epigenetic perspective in patients with TLE-HS type I (TLE-HSTI) and TLE without HS (TLE-nHS) using whole-genome bisulfite sequencing (WGBS). We defined 1171 hypermethylated and 2537 hypomethylated regions and found 632 differentially methylated genes (DMG) in the promoter region that were primarily involved in the regulation of various aspects of epilepsy development. Twelve DMG overlapped with differentially expressed genes (DEG) in the promoter region, and RT-qPCR findings revealed significant overexpression of the SBNO2, CBX3, RASAL3, and TMBIM4 genes in TLE-HSTI. We present the first systematic analysis of methylation profiles of TLE-HSTI and TLE-nHS from an epigenetic perspective using WGBS. Overall, our preliminary data highlight the underlying mechanism of TLE-HSTI, providing a new perspective for guiding treatment of TLE.

Detection of Disease-Causing SNVs/Indels and CNVs in Single Test Based on Whole Exome Sequencing: A Retrospective Case Study in Epileptic Encephalopathies

Background: Epileptic encephalopathies (EEs) are a pediatric entity with highly phenotypic and genetic heterogeneity. Both single nucleotide variants (SNVs)/Indels and copy number variations (CNVs) could be the causes. Whole exome sequencing (WES) is widely applied to detect SNVs/Indels, but the bioinformatics approach for detecting CNVs is still limited and weak. In the current study, the possibility of profiling both disease-causing SNVs/Indels and CNVs in a single test based on WES in EEs was evaluated. Methods: The infants diagnosed with EEs were enrolled from a single pediatric epilepsy center between January 2018 and February 2020. Demographic and clinical data were collected. In WES data, the pathogenic SNVs were identified through an in-house pipeline, and pathogenic CNVs were identified by CNVkit. The diagnostic rate was evaluated, and the molecular findings were characterized. Results: A total of 73 infants were included; 36 (49.32%) of them were males. The median age was 7 months. Thirty-two (43.84%) infants had been diagnosed with epilepsy syndrome. The most common type of syndrome was West syndrome (22/73, 30.1%), followed by Dravet syndrome (20/77, 27.4%). Fifty-four (73.97%) had intellectual development delay. The genetic cause of EEs, pathogenic or likely pathogenic variants, were successfully discovered in 46.6% (34/73) of the infants, and 29 (39.7%) infants carried SNVs/Indels, while 5 (6.8%) carried CNVs. The majority of the disease-causing variants were inherited in de novo pattern (25, 71.4%). In addition to showing that the variants in the ion channel encoding genes accounted for the main etiology, we discovered and confirmed two new disease-causing genes, CACNA1E and WDR26. Five discovered CNVs were deletions of 2q24.3, 1p36, 15q11-q13, 16p11.2, and 17p13.3, and all were confirmed by array comparative genomic hybridization. Conclusion: The application of both SNVs/Indels and CNVs detection in a single test based on WES yielded a high diagnosis rate in EEs. WES may serve as a first-tier test with cost-effective benefit in EEs.

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.

Clinical Utility of Exome Sequencing and Reinterpreting Genetic Test Results in Children and Adults With Epilepsy

The clinical utility of genetic testing for epilepsy has been enhanced with the advancement of next-generation sequencing (NGS) technology along with the rapid updating of publicly available databases. The aim of this study was to evaluate the diagnostic yield of NGS and assess the value of reinterpreting genetic test results in children and adults with epilepsy. We performed genetic testing on 200 patients, including 82 children and 118 adults. The results were classified into three categories: positive, inconclusive, or negative. The reinterpretation of inconclusive results was conducted in April 2020. Overall, we identified disease-causing variants in 12% of the patients in the original analysis, and 14.5% at reinterpretation. The diagnostic yield for adults with epilepsy was similar to that for children (11 vs. 19.5%, p = 0.145). After reinterpretation, 9 of the 86 patients who initially had inconclusive results obtained a clinically significant change in diagnosis. Among these nine revised cases, five obtained positive diagnoses, representing a diagnosis rate of 5.8% (5/86). Manual searches for additional evidence of pathogenicity for candidate variants and updated patient clinical information were the main reasons for diagnostic reclassification. This study emphasizes the diagnostic potential of combining NGS and reinterpretation of inconclusive genetic test reports in children and adults with epilepsy.

Cognitive performance and behavior across idiopathic/genetic epilepsies in children and adolescents

We investigated the cognitive and behavioral profile of three distinct groups of epilepsies with a genetic background for intergroup differences: (1) idiopathic/genetic generalized epilepsies (IGE/GGE group); (2) idiopathic focal epilepsies (IFE group); and (3) epilepsies with proven or strongly suggested monogenic or structural/numeric chromosomal etiology (genetic epilepsies, GE group). Cognitive (total IQ and subcategories) and behavioral parameters (CBCL) were assessed at the tertiary epilepsy center of the University of Munich (Germany). We used ANOVA with post-hoc Bonferroni-correction to explore significant mean differences and Fisher’s exact test for significant proportional differences of intelligence impairment and behavioral problems. 126 (56 IGE/GGE, 26 IFE, 44 GE) patients were available. Total IQ was 89.0 ± 15.9 (95% CI 84.5–93.4) for IGE/GGE, 94.8 ± 18.1 (95% CI 87.3–102.3) for IFE and 76.4 ± 22.4 (95% CI 67.6–85.3) for GE (p = 0.001). The same trend was significant for all but one IQ subcategory. The rate of patients with an intelligence impairment (total IQ < 70) was higher for GE (40%) than for IGE/GGE (14%) and for IFE (7%) patients (p = 0.033). There were no significant differences between groups for behavior scores and behavioral problems. This study shows that the current ILAE classification of epilepsies with genetic etiology creates a heterogeneous group of patients with respect to cognitive performance but not behavior. These findings may help in further delineating epilepsies as regards cognitive performance, notwithstanding their closely related etiological classification.

The burden of rare damaging variants in hereditary atypical parkinsonism genes is increased in patients with Parkinson's disease

Increased burdens of rare coding variants in genes related to lysosomal storage disease or mitochondrial pathways were reported to be associated with idiopathic Parkinson's disease. Under a hypothesis that the burden of damaging rare coding variants is increased in causative genes for hereditary parkinsonism, we analyzed the burdens of rare coding variants with a case-control design. Two cohorts of whole-exome sequencing data and a cohort of genome-wide genotyping data of clinically validated idiopathic Parkinson's disease cases and controls, which were open to the public, were used. The sequence kernel association test-optimal was used to analyze the burden of rare variants in the hereditary parkinsonism gene set, which was constructed from the Online Mendelian Inheritance in Man database through manual curation. The hereditary parkinsonism gene set consisted of 17 genes with a locus symbol prefix for familial Parkinson's disease and 75 hereditary atypical parkinsonism genes. We detected a significant association of enriched burdens of predicted damaging rare coding variants in hereditary parkinsonism genes in all three datasets. Meta-analyses of the rare variant burden test in a subgroup of gene sets revealed an association between burdens of rare damaging variants with PD in a hereditary atypical parkinsonism gene set, but not in a subgroup gene set with a locus symbol prefix for familial Parkinson's disease. Our results highlight the roles of rare damaging variants in causative genes for hereditary atypical parkinsonian disorders. We propose that Mendelian genes associated with hereditary disorders accompanying parkinsonism are involved in Parkinson's disease-related genetic networks.

Gabra2 is a genetic modifier of Scn8a encephalopathy in the mouse

Objective

SCN8A encephalopathy is a developmental epileptic encephalopathy typically caused by de novo gain‐of‐function mutations in Na_v1.6. Severely affected individuals exhibit refractory seizures, developmental delay, cognitive disabilities, movement disorders, and elevated risk of sudden death. Patients with the identical SCN8A variant can differ in clinical course, suggesting a role for modifier genes in determining disease severity. The identification of genetic modifiers contributes to understanding disease pathogenesis and suggesting therapeutic interventions.

Methods

We generated F1 and F2 crosses between inbred mouse strains and mice carrying the human pathogenic variants SCN8A‐R1872W and SCN8A‐N1768D. Quantitative trait locus (QTL) analysis of seizure‐related phenotypes was used for chromosomal mapping of modifier loci.

Results

In an F2 cross between strain SJL/J and C57BL/6J mice carrying the patient mutation R1872W, we identified a major QTL on chromosome 5 containing the Gabra2 gene. Strain C57BL/6J carries a splice site mutation that reduces expression of Gabra2, encoding the α2 subunit of the aminobutyric acid type A receptor. The protective wild‐type allele of Gabra2 from strain SJL/J delays the age at seizure onset and extends life span of the Scn8a mutant mice. Additional Scn8a modifiers were observed in the F2 cross and in an F1 cross with strain C3HeB/FeJ.

Significance

These studies demonstrate that the SJL/J strain carries multiple modifiers with protective effects against seizures induced by gain‐of‐function mutations in Scn8a. Homozygosity for the hypomorphic variant of Gabra2 in strain C57BL/6J is associated with early seizure onset and short life span. GABRA2 is a potential therapeutic target for SCN8A encephalopathy.

[Genetics of epilepsy: successes, problems and development prospects]

The authors present a detailed review of current advances in the field of genetics of epilepsy. Separately, new views on the etiology and pathogenesis of genetic epileptic encephalopathies, focal epilepsy and idiopathic generalized epilepsies are examined. The authors emphasize the importance of genetic discoveries for the clinical practice, including the prospects in the development of patients' personalized treatment. A comparative analysis of the value of various methods of genetic research in the diagnosis of epilepsy, methods of integrating molecular genetic analyses into everyday practical medicine is presented.