Evaluating the pathogenic potential of genes with de novo variants in epileptic encephalopathies

MED12 variants associated with X-linked recessive partial epilepsy without intellectual disability

OBJECTIVES

The MED12 gene encodes mediator complex subunit 12, which is a component of the mediator complex involved in the transcriptional regulation of nearly all RNA polymerase II-dependent genes. MED12 variants have previously been associated with developmental disorders with or without nonspecific intellectual disability. This study aims to explore the association between MED12 variants and epilepsy.

MATERIALS AND METHODS

Trios-based whole-exome sequencing was performed in a cohort of 349 unrelated cases with partial (focal) epilepsy without acquired causes. The genotype-phenotype correlations of MED12 variants were analyzed.

RESULTS

Five hemizygous missense MED12 variants, including c.958A>G/p.Ile320Val, c.1757G>A/p.Ser586Asn, c.2138C>T/p.Pro713Leu, c.3379T>C/p.Ser1127Pro, and c.4219A>C/p.Met1407Leu were identified in five unrelated males with partial epilepsy. All patients showed infrequent focal seizures and achieved seizure free without developmental abnormalities or intellectual disability. All the hemizygous variants were inherited from asymptomatic mothers (consistent with the X-linked recessive inheritance pattern) and were absent in the general population. The two variants with damaging hydrogen bonds were associated with early-onset seizures. Further genotype-phenotype analysis revealed that congenital anomaly disorder (Hardikar syndrome) was associated with (de novo) destructive variants in an X-linked dominant inheritance pattern, whereas epilepsy was associated with missense variants in an X-linked recessive inheritance pattern. Phenotypic features of intellectual disability appeared as the intermediate phenotype in terms of both genotype and inheritance. Epilepsy-related variants were located at the MED12-LCEWAV domain and the regions between MED12-LCEWAV and MED12-POL.

CONCLUSION

MED12 is a potentially causative gene for X-linked recessive partial epilepsy without developmental or intellectual abnormalities. The genotype-phenotype correlation of MED12 variants explains the phenotypic variations and can help the genetic diagnosis.

NEXMIF variants are associated with epilepsy with or without intellectual disability

OBJECTIVES

Variants in NEXMIF had been reported associated with intellectual disability (ID) without epilepsy or developmental epileptic encephalopathy (DEE). It is unkown whether NEXMIF variants are associated with epilepsy without ID. This study aims to explore the phenotypic spectrum of NEXMIF and the genotype-phenotype correlations.

MATERIALS AND METHODS

Trio-based whole-exome sequencing was performed in patients with epilepsy. Previously reported NEXMIF variants were systematically reviewed to analyze the genotype-phenotype correlations.

RESULTS

Six variants were identified in seven unrelated cases with epilepsy, including two de novo null variants and four hemizygous missense variants. The two de novo variants were absent in all populations of gnomAD and four hemizygous missense variants were absent in male controls of gnomAD. The two patients with de novo null variants exhibited severe developmental epileptic encephalopathy. While, the patients with hemizygous missense variants had mild focal epilepsy with favorable outcome. Analysis of previously reported cases revealed that males with missense variants presented significantly higher percentage of normal intellectual development and later onset age of seizure than those with null variants, indicating a genotype-phenotype correlation.

CONCLUSION

This study suggested that NEXMIF variants were potentially associated with pure epilepsy with or without intellectual disability. The spectrum of epileptic phenotypes ranged from the mild epilepsy to severe developmental epileptic encephalopathy, where the epileptic phenotypes variability are potentially associated with patients' gender and variant type.

Biallelic structural variations within FGF12 detected by long-read sequencing in epilepsy

We discovered biallelic intragenic structural variations (SVs) in FGF12 by applying long-read whole genome sequencing to an exome-negative patient with developmental and epileptic encephalopathy (DEE). We also found another DEE patient carrying a biallelic (homozygous) single-nucleotide variant (SNV) in FGF12 that was detected by exome sequencing. FGF12 heterozygous recurrent missense variants with gain-of-function or heterozygous entire duplication of FGF12 are known causes of epilepsy, but biallelic SNVs/SVs have never been described. FGF12 encodes intracellular proteins interacting with the C-terminal domain of the alpha subunit of voltage-gated sodium channels 1.2, 1.5, and 1.6, promoting excitability by delaying fast inactivation of the channels. To validate the molecular pathomechanisms of these biallelic FGF12 SVs/SNV, highly sensitive gene expression analyses using lymphoblastoid cells from the patient with biallelic SVs, structural considerations, and Drosophila in vivo functional analysis of the SNV were performed, confirming loss-of-function. Our study highlights the importance of small SVs in Mendelian disorders, which may be overlooked by exome sequencing but can be detected efficiently by long-read whole genome sequencing, providing new insights into the pathomechanisms of human diseases.

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.

Detecting low-variant allele frequency mosaic pathogenic variants of NF1, TSC2, and AKT3 genes from blood in patients with neurodevelopmental disorders

Growing evidence indicates that early, and late postzygotic mosaicism can cause neurodevelopmental disorders (NDD), but detection of low variant allele frequency (VAF) mosaic variants from blood remains a challenge. We reviewed data of 2,162 patients with NDDs who underwent conventional genetic tests and performed a deep sequencing using specifically designed mosaic NGS panel in the patients with negative genetic test results. Forty-four patents with neurocutaneous syndrome, malformation of cortical development or nonlesional epileptic encephalopathies were included. In total, mosaic variants were detected from blood in 1.2% (25/2,162) of the patients. Using conventional NGS panels, 22 mosaic variants (VAF 8.8-29.8%) were identified in 18 different genes including TSC2, DCX, SLC2A1, PCDH19, DNM1, STXBP1, SCN2A, SCN1A, PURA, POGZ, PAFAH1B1, NF1, KIF21A, KCNQ2, GABRA1, EEF1A2, CDKL5, and ARID1B. Using a specifically designed mosaicism NGS panel, three mosaic variants of the NF1, TSC2, and AKT3 genes were identified (VAF 2.0-11.2%). Mosaic variants were found frequently in the patients who had neurocutaneous syndrome (2/7, 28.6%) whereas only one or no mosaic variant was detected for patients who had malformations of cortical development (1/20, 5%) or nonlesional epileptic encephalopathies (0%, 0/17). In summary, mosaic variants contribute to spectrum of NDDs can be detected from blood via the conventional NGS and specifically designed mosaicism NGS panels, and detection of mosaic variants using blood will increase diagnostic yield.

Insights from Spatial Measures of Intolerance to Identifying Pathogenic Variants in Developmental and Epileptic Encephalopathies

Developmental and epileptic encephalopathies (DEEs) are a group of epilepsies with early onset and severe symptoms that sometimes lead to death. Although previous work successfully discovered several genes implicated in disease outcomes, it remains challenging to identify causative mutations within these genes from the background variation present in all individuals due to disease heterogeneity. Nevertheless, our ability to detect possible pathogenic variants has continued to improve as in silico predictors of deleteriousness have advanced. We investigate their use in prioritising likely pathogenic variants in epileptic encephalopathy patients’ whole exome sequences. We showed that the inclusion of structure-based predictors of intolerance improved upon previous attempts to demonstrate enrichment within epilepsy genes.

Emerging evidence of genotype–phenotype associations of developmental and epileptic encephalopathy due to KCNC2 mutation: Identification of novel R405G

Developmental and epileptic encephalopathies (DEEs) have high genetic heterogeneity, and DEE due to the potassium voltage-gated channel subfamily C member 2 (KCNC2) variant remains poorly understood, given the scarcity of related case studies. We report on two unrelated Chinese patients, an 11-year-old boy and a 5-year-old girl, diagnosed with global developmental delay (GDD), intellectual disability (ID), and focal impaired awareness seizure characterized by generalized spike and wave complexes on electroencephalogram (EEG) in the absence of significant brain lesions. Whole-exome sequencing (WES) and electrophysiological analysis were performed to detect genetic variants and evaluate functional changes of the mutant KCNC2, respectively. Importantly, we identified a novel gain-of-function KCNC2 variant, R405G, in both patients. Previously reported variants, V471L, R351K, T437A, and T437N, and novel R405G were found in multiple unrelated patients with DEE, showing consistent genotype–phenotype associations. These findings emphasize that the KCNC2 gene is causative for DEE and facilitates treatment and prognosis in patients with DEE due to KCNC2 mutations.

Favorable Response to "Memantine" in a Child with GRIN2B Epileptic Encephalopathy

GRIN2B is a gene encoding GluN2B subunit under the family of N-methyl D-aspartate (NMDA) receptors, which is responsible for neurogenesis and cognitive processes. The role of NMDA receptor antagonists like memantine is being explored for therapies in drug-resistant epilepsies. Here, we present a case of a 20-month-old boy who presented with refractory epileptic spasms. Upon failure of multiple antiepileptic drugs, he was started on oral memantine. There was a significant reduction in average seizure episodes by ∼80%. The use of memantine along with antiepileptic drug polytherapy has proved to be beneficial in our case. Our experience with memantine and favorable outcome opens up the scope of more research into the use of NMDA receptor antagonist as a drug option for refractory epilepsies with proven genetic mutation and hence improves the overall neurodevelopmental outcome and survival chance.

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.

Detection of Deregulated miRNAs in Childhood Epileptic Encephalopathies

The term "epileptic encephalopathy" is used to describe a possible relationship between epilepsy and developmental delay. The pathogenesis of developmental encephalopathies, independent of epilepsy, can be defined by genetic control mechanisms. The aim of this study was to investigate the use of miRNAs as serum biomarkers for the determination and discrimination of epileptic encephalopathies. Whole blood samples obtained from 54 individuals in 2 groups designated as epileptic encephalopathy patients' group (n = 24) and healthy controls (n = 30) were included in this study. The expression levels of 10 miRNAs were determined using qRT-PCR. After the determination of expression levels, the correlation of upregulated miRNA levels and Ki67 index was calculated using Pearson correlation test. The comparison of epileptic encephalopathy patients' group with healthy controls revealed the upregulation of one miRNAs (hsa-miR-324-5p) and downregulation of three miRNAs (hsa-miR-146a-5p, hsa-miR-138-5p, hsa-miR-187-3p). It has been determined that miRNAs with altered expression are an important factor in the formation of epileptic seizures and seizure-induced neuronal death. The fact that processes that play a key role in epiloptogenesis are under the control of miRNAs causes miRNAs to become meta-controllers of gene expression in the brain. We thought that further studies are needed to prove that especially hsa-miR-146a-5p, hsa-miR-138-5p, and hsa-miR-187-3p can be used as epileptic encephalopathy biomarkers. The detection of disease-specific miRNAs could contribute to the development of precision treatments.

Critical Role of E1623 Residue in S3-S4 Loop of Nav1.1 Channel and Correlation Between Nature of Substitution and Functional Alteration

Objective: An overwhelming majority of the genetic variants associated with genetic disorders are missense. The association between the nature of substitution and the functional alteration, which is critical in determining the pathogenicity of variants, remains largely unknown. With a novel missense variant (E1623A) identified from two epileptic cases, which occurs in the extracellular S3-S4 loop of Na_v1.1, we studied functional changes of all latent mutations at residue E1623, aiming to understand the relationship between substitution nature and functional alteration.

Methods: Six latent mutants with amino acid substitutions at E1623 were generated, followed by measurements of their electrophysiological alterations. Different computational analyses were used to parameterize the residue alterations.

Results: Structural modeling indicated that the E1623 was located in the peripheral region far from the central pore, and contributed to the tight turn of the S3-S4 loop. The E1623 residue exhibited low functional tolerance to the substitutions with the most remarkable loss-of-function found in E1623A, including reduced current density, less steady-state availability of activation and inactivation, and slower recovery from fast inactivation. Correlation analysis between electrophysiological parameters and the parameterized physicochemical properties of different residues suggested that hydrophilicity of side-chain at E1623 might be a crucial contributor for voltage-dependent kinetics. However, none of the established algorithms on the physicochemical variations of residues could well predict changes in the channel conductance property indicated by peak current density.

Significance: The results established the important role of the extracellular S3-S4 loop in Na_v1.1 channel gating and proposed a possible effect of local conformational loop flexibility on channel conductance and kinetics. Site-specific knowledge of protein will be a fundamental task for future bioinformatics.

A recurrent de novo variant supports KCNC2 involvement in the pathogenesis of developmental and epileptic encephalopathy

Developmental and epileptic encephalopathies (DEE) are a heterogenous group of conditions characterized by the co-occurrence of epilepsy and intellectual/developmental disability. Despite several known DEE-related genes, including these encoding ion channels, still many cases remain without molecular diagnosis. Here, we present a 2-year-old girl with severe DEE in whom whole exome sequencing revealed de novo p.(Val471Leu) variant in the KCNC2 encoding Kv3.2, a voltage-gated potassium channel. To the best of our knowledge, this is the third DEE case due to KCNC2 mutation. Our clinical and molecular findings, particularly the recurrence of p.(Val471Leu) in patient with similar clinical phenotype, further support KCNC2 as a novel DEE-associated gene.

DNA Methylation Signature of Epileptic Encephalopathy-Related Pathogenic Genes Encoding Ion Channels in Temporal Lobe Epilepsy

Epilepsy is characterized by highly abnormal synchronous discharge of brain neurons, and ion channels are fundamental in the generation and modulation of neural excitability. Considering that abnormal methylation can either activate or repress genes, this study was designed to explore the DNA methylation signature of pathogenic genes encoding ion channels in temporal lobe epilepsy (TLE). In total, 38 TLE patients and 38 healthy controls were enrolled in the study, and genomic DNA and total protein of the lymphocytes were extracted from peripheral blood samples to assess methylation and protein levels. The DNA methylation levels of all 12 genes examined were significantly lower in the TLE group than in the control group. After false-positive correction, 83.3% (10/12) of these genes, namely, gamma-aminobutyric acid type A receptor subunit beta1 (GABRB1), gamma-aminobutyric acid type A receptor subunit beta2 (GABRB2), gamma-aminobutyric acid type A receptor subunit beta1 (GABRB3), glutamate ionotropic receptor NMDA type subunit 1 (GRIN1), glutamate ionotropic receptor NMDA type subunit 2A (GRIN2A), glutamate ionotropic receptor NMDA type subunit 2B (GRIN2B), hyperpolarization activated cyclic nucleotide gated potassium channel 1 (HCN1), potassium voltage-gated channel subfamily A member 2 (KCNA2), potassium voltage-gated channel subfamily B member 1 (KCNB1), and potassium sodium-activated channel subfamily T member 1 (KCNT1), were still differentially expressed. Among these ion channels, HCN1 and KCNA2 were selected to evaluate the effects of DNA methylation, and the levels of these proteins were inversely upregulated in the TLE group compared to the control group. As the genes identified as having differential methylation levels are involved in both excitatory and inhibitory ion channels, this study observed by binary logistic regression that hypermethylated GARAB1 was an independent risk factor for TLE, indicating that the overwhelming effect of ion channels on TLE is probably inhibitory from the perspective of DNA methylation. All these findings support the involvement of DNA methylation in TLE pathologies, but the mechanisms need to be further investigated.

Loss-of-function variants in DNM1 cause a specific form of developmental and epileptic encephalopathy only in biallelic state

Background

Developmental and epileptic encephalopathies (DEEs) represent a group of severe neurological disorders characterised by an onset of refractory seizures during infancy or early childhood accompanied by psychomotor developmental delay or regression. DEEs are genetically heterogeneous with, to date, more than 80 different genetic subtypes including DEE31 caused by heterozygous missense variants in DNM1.

Methods

We performed a detailed clinical characterisation of two unrelated patients with DEE and used whole-exome sequencing to identify causative variants in these individuals. The identified variants were tested for cosegregation in the respective families.

Results

We excluded pathogenic variants in known, DEE-associated genes. We identified homozygous nonsense variants, c.97C>T; p.(Gln33*) in family 1 and c.850C>T; p.(Gln284*) in family 2, in the DNM1 gene, indicating that biallelic, loss-of-function pathogenic variants in DNM1 cause DEE.

Conclusion

Our finding that homozygous, loss-of-function variants in DNM1 cause DEE expands the spectrum of pathogenic variants in DNM1. All parents who were heterozygous carriers of the identified loss-of-function variants were healthy and did not show any clinical symptoms, indicating that the type of mutation in DNM1 determines the pattern of inheritance.

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.

UNC13B variants associated with partial epilepsy with favourable outcome

The unc-13 homolog B (UNC13B) gene encodes a presynaptic protein, mammalian uncoordinated 13-2 (Munc13-2), that is highly expressed in the brain-predominantly in the cerebral cortex-and plays an essential role in synaptic vesicle priming and fusion, potentially affecting neuronal excitability. However, the functional significance of UNC13B mutation in human disease is not known. In this study we screened for novel genetic variants in a cohort of 446 unrelated cases (families) with partial epilepsy without acquired causes by trio-based whole-exome sequencing. UNC13B variants were identified in 12 individuals affected by partial epilepsy and/or febrile seizures from eight unrelated families. The eight probands all had focal seizures and focal discharges in EEG recordings, including two patients who experienced frequent daily seizures and one who showed abnormalities in the hippocampus by brain MRI; however, all of the patients showed favorable outcome without intellectual or developmental abnormalities. The identified UNC13B variants included one nonsense variant, two variants at or around a splice site, one compound heterozygous missense variant, and four missense variants that cosegregated in the families. The frequency of UNC13B variants identified in the present study was significantly higher than that in a control cohort of Han Chinese and controls of the East Asian and all populations in the Genome Aggregation Database. Computational modeling, including hydrogen bond and docking analyses, suggested that the variants lead to functional impairment. In Drosophila, seizure rate and duration were increased by Unc13b knockdown compared to wild-type flies, but these effects were less pronounced than in sodium voltage-gated channel alpha subunit 1 (Scn1a) knockdown Drosophila. Electrophysiologic recordings showed that excitatory neurons in Unc13b-deficient flies exhibited increased excitability. These results suggest that UNC13B is potentially associated with epilepsy. The frequent daily seizures and hippocampal abnormalities but ultimately favorable outcome under antiepileptic therapy in our patients indicate that partial epilepsy caused by UNC13B variant is a clinically manageable condition.

Genetic Neonatal-Onset Epilepsies and Developmental/Epileptic Encephalopathies with Movement Disorders: A Systematic Review

Despite expanding next generation sequencing technologies and increasing clinical interest into complex neurologic phenotypes associating epilepsies and developmental/epileptic encephalopathies (DE/EE) with movement disorders (MD), these monogenic conditions have been less extensively investigated in the neonatal period compared to infancy. We reviewed the medical literature in the study period 2000–2020 to report on monogenic conditions characterized by neonatal onset epilepsy and/or DE/EE and development of an MD, and described their electroclinical, genetic and neuroimaging spectra. In accordance with a PRISMA statement, we created a data collection sheet and a protocol specifying inclusion and exclusion criteria. A total of 28 different genes (from 49 papers) leading to neonatal-onset DE/EE with multiple seizure types, mainly featuring tonic and myoclonic, but also focal motor seizures and a hyperkinetic MD in 89% of conditions, with neonatal onset in 22%, were identified. Neonatal seizure semiology, or MD age of onset, were not always available. The rate of hypokinetic MD was low, and was described from the neonatal period only, with WW domain containing oxidoreductase (WWOX) pathogenic variants. The outcome is characterized by high rates of associated neurodevelopmental disorders and microcephaly. Brain MRI findings are either normal or nonspecific in most conditions, but serial imaging can be necessary in order to detect progressive abnormalities. We found high genetic heterogeneity and low numbers of described patients. Neurological phenotypes are complex, reflecting the involvement of genes necessary for early brain development. Future studies should focus on accurate neonatal epileptic phenotyping, and detailed description of semiology and time-course, of the associated MD, especially for the rarest conditions.

Genetic Effects of the Schizophrenia-Related Gene DTNBP1 in Temporal Lobe Epilepsy

Recent studies have reported patients who concurrently exhibit conditions of epilepsy and schizophrenia, indicating certain shared pathologies between them. This study aimed to investigate the genetic effects of the schizophrenia-related gene DTNBP1 in temporal lobe epilepsy (TLE). A total of 496 TLE patients and 528 healthy individuals were successfully genotyped for six DTNBP1 polymorphisms (rs760665, rs1011313, rs2619528, rs2619522, rs909706, and rs2619538), including 335 TLE patients and 325 healthy controls in cohort 1, and 161 TLE patients and 203 healthy controls in cohort 2. The frequency of the TT genotype at rs909706 T > C was lower in TLE patients than in normal controls in the initial cohort (cohort 1), which was confirmed in an independent cohort (cohort 2). However, the intronic T allele failed to be in linkage disequilibrium (LD) with any functional variations nearby; thus, together with the CCAC and TCAT haplotypes (rs1011313-rs2619528-rs2619522-rs909706) observed in the study, this allele acts only as a protective factor against susceptibility to TLE. Meanwhile, a novo mutant allele rs2619538 T > A was exclusively observed in TLE patients, and a dual-luciferase assay revealed that the mutant allele was increased by approximately 22% in the DTNBP2 promoter compared with the wild-type allele. Together with the trend of increasing DTNBP1 expression in epilepsy patients and animal models in this study, these are the first findings to demonstrate the genetic association of DTNBP1 with TLE. Homozygous mutation of rs2619538 T > A likely promotes DTNBP1 expression and facilitates subsequent processes in epilepsy pathologies. Thus, the role of DTNBP1 in TLE deserves further exploration in the future.

What are the epileptic encephalopathies?

PURPOSE OF REVIEW

To review the evolution of the concept of epileptic encephalopathy during the course of past years and analyze how the current definition might impact on both clinical practice and research.

RECENT FINDINGS

Developmental delay in children with epilepsy could be the expression of the cause, consequence of intense epileptiform activity (seizures and EEG abnormalities), or because of the combination of both factors. Therefore, the current International League Against Epilepsy classification identified three electroclinical entities that are those of developmental encephalopathy, epileptic encephalopathy, and developmental and epileptic encephalopathy (DEE). Many biological pathways could be involved in the pathogenesis of DEEs. DNA repair, transcriptional regulation, axon myelination, metabolite and ion transport, and peroxisomal function could all be involved in DEE. Also, epilepsy and epileptiform discharges might impact on cognition via several mechanisms, although they are not fully understood.

SUMMARY

The correct and early identification of cause in DEE might increase the chances of a targeted treatment regimen. Interfering with neurobiological processes of the disease will be the most successful way in order to improve both the cognitive disturbances and epilepsy that are the key features of DEE.

Leveraging correlations between variants in polygenic risk scores to detect heterogeneity in GWAS cohorts

Evidence from both GWAS and clinical observation has suggested that certain psychiatric, metabolic, and autoimmune diseases are heterogeneous, comprising multiple subtypes with distinct genomic etiologies and Polygenic Risk Scores (PRS). However, the presence of subtypes within many phenotypes is frequently unknown. We present CLiP (Correlated Liability Predictors), a method to detect heterogeneity in single GWAS cohorts. CLiP calculates a weighted sum of correlations between SNPs contributing to a PRS on the case/control liability scale. We demonstrate mathematically and through simulation that among i.i.d. homogeneous cases generated by a liability threshold model, significant anti-correlations are expected between otherwise independent predictors due to ascertainment on the hidden liability score. In the presence of heterogeneity from distinct etiologies, confounding by covariates, or mislabeling, these correlation patterns are altered predictably. We further extend our method to two additional association study designs: CLiP-X for quantitative predictors in applications such as transcriptome-wide association, and CLiP-Y for quantitative phenotypes, where there is no clear distinction between cases and controls. Through simulations, we demonstrate that CLiP and its extensions reliably distinguish between homogeneous and heterogeneous cohorts when the PRS explains as low as 3% of variance on the liability scale and cohorts comprise 50, 000 − 100, 000 samples, an increasingly practical size for modern GWAS. We apply CLiP to heterogeneity detection in schizophrenia cohorts totaling > 50, 000 cases and controls collected by the Psychiatric Genomics Consortium. We observe significant heterogeneity in mega-analysis of the combined PGC data (p-value 8.54 × 0⁻⁴), as well as in individual cohorts meta-analyzed using Fisher’s method (p-value 0.03), based on significantly associated variants. We also apply CLiP-Y to detect heterogeneity in neuroticism in over 10, 000 individuals from the UK Biobank and detect heterogeneity with a p-value of 1.68 × 10⁻⁹. Scores were not significantly reduced when partitioning by known subclusters (“Depression” and “Worry”), suggesting that these factors are not the primary source of observed heterogeneity.

Several traits, such as bipolar disease, are known to be heterogeneous and comprise distinct subtypes with unique genomic associations. For other traits such as schizophrenia, heterogeneity may be suspected, but specific subtypes are less well characterized. Furthermore, conventional mixture model-based methods to detect subtypes in genome-wide association data struggle with the high polygenicity of complex traits. We propose CLiP (Correlated Liability Predictors), a method that does not identify subtype-specific effects, but is very well-powered to detect heterogeneity of any kind within the very weak signals of GWAS. CLiP serves as a method to flag particular phenotypes for potential further study into the genomic factors driving heterogeneity, as well as a means to evaluate the transferability of polygenic risk scores. We also develop extensions of CLiP applicable to scoring heterogeneity in quantitative phenotypes and quantitative predictors such as gene expression. We apply CLiP to scoring heterogeneity in schizophrenia cohorts from the Psychiatric Genomics Consortium and neuroticism in individuals in the UK Biobank and find significant heterogeneity in both phenotypes, warranting further study.

DEPDC5 Variants Associated Malformations of Cortical Development and Focal Epilepsy With Febrile Seizure Plus/Febrile Seizures: The Role of Molecular Sub-Regional Effect

To explore the phenotype spectrum of DEPDC5 variants and the possible mechanisms underlying phenotypical variation, we performed targeted next-generation sequencing in 305 patients with focal epilepsies and 91 patients with generalized epilepsies. Protein modeling was performed to predict the effects of missense mutations. All previously reported epilepsy-related DEPDC5 variants were reviewed. The genotype–phenotype correlations with molecular sub-regional implications were analyzed. We identified a homozygous DEPDC5 mutation (p.Pro1031His) in a case with focal cortical dysplasia and eight heterozygous mutations in 11 families with mild focal epilepsies, including 13 patients in eight families with focal epilepsy with febrile seizures plus/febrile seizures (FEFS + /FS). The mutations included one termination codon mutation (p.Ser1601_Ter1604del_ext133), three truncating mutations (p.Val151Serfs^∗27, p.Arg239^∗, and p.Arg838^∗), and four missense mutations (p.Tyr7Cys, p.Tyr836Cys, p.Pro1031His, and p.Gly1545Ser) that were predicted to affect hydrogen bonds and protein stability. Analysis on epilepsy-related DEPDC5 variants revealed that malformations of cortical development (MCDs) had a tendency of higher frequency of null mutations than those without MCD. MCD-associated heterozygous missense mutations were clustered in structural axis for binding arrangement (SABA) domain and close to the binding sites to NPRL2/NPRL3 complex, whereas those associated with FEFS + /FS were a distance away from the binding sites. Evidence from four aspects and one possible evidence from sub-regional implication suggested MCD and FEFS + /FS as phenotypes of DEPDC5 variants. This study suggested that the phenotypes of DEPDC5 variants vary from mild FEFS + /FS to severe MCD. Heterozygous DEPDC5 mutations are generally less pathogenic and commonly associated with mild phenotypes. Bi-allelic mutations and second hit of somatic mutations, together with the genotype–phenotype correlation and sub-regional implication of DEPDC5 variants, explain severe phenotypes.

EEF1A2 mutations in epileptic encephalopathy/intellectual disability: Understanding the potential mechanism of phenotypic variation

EEF1A2 encodes protein elongation factor 1-alpha 2, which is involved in Guanosine triphosphate (GTP)-dependent binding of aminoacyl-transfer RNA (tRNA) to the A-site of ribosomes during protein biosynthesis and is highly expressed in the central nervous system. De novo mutations in EEF1A2 have been identified in patients with extensive neurological deficits, including intractable epilepsy, globe developmental delay, and severe intellectual disability. However, the mechanism underlying phenotype variation is unknown. Using next-generation sequencing, we identified a novel and a recurrent de novo mutation, c.294C>A; p.(Phe98Leu) and c.208G>A; p.(Gly70Ser), in patients with Lennox-Gastaut syndrome. The further systematic analysis revealed that all EEF1A2 mutations were associated with epilepsy and intellectual disability, suggesting its critical role in neurodevelopment. Missense mutations with severe molecular alteration in the t-RNA binding sites or GTP hydrolysis domain were associated with early-onset severe epilepsy, indicating that the clinical expression was potentially determined by the location of mutations and alteration of molecular effects. This study highlights the potential genotype-phenotype relationship in EEF1A2 and facilitates the evaluation of the pathogenicity of EEF1A2 mutations in clinical practice.

A 2020 View on the Genetics of Developmental and Epileptic Encephalopathies

Developmental and epileptic encephalopathies (DEEs) can be primarily attributed to genetic causes. The genetic landscape of DEEs has been largely shaped by the rise of high-throughput sequencing, which led to the discovery of new DEE-associated genes and helped identify de novo pathogenic variants. We discuss briefly the contribution of de novo variants to DEE and also focus on alternative inheritance models that contribute to DEE. First, autosomal recessive inheritance in outbred populations may have a larger contribution than previously appreciated, accounting for up to 13% of DEEs. A small subset of genes that typically harbor de novo variants have been associated with recessive inheritance, and often these individuals have more severe clinical presentations. Additionally, pathogenic variants in X-linked genes have been identified in both affected males and females, possibly due to a lack of X-chromosome inactivation skewing. Collectively, exome sequencing has resulted in a molecular diagnosis for many individuals with DEE, but this still leaves many cases unsolved. Multiple factors contribute to the missing etiology, including nonexonic variants, mosaicism, epigenetics, and oligogenic inheritance. Here, we focus on the first 2 factors. We discuss the promises and challenges of genome sequencing, which allows for a more comprehensive analysis of the genome, including interpretation of structural and noncoding variants and also yields a high number of de novo variants for interpretation. We also consider the contribution of genetic mosaicism, both what it means for a molecular diagnosis in mosaic individuals and the important implications for genetic counseling.

Optimization of in silico tools for predicting genetic variants: individualizing for genes with molecular sub-regional stratification

Genes are unique in functional role and differ in their sensitivities to genetic defects, but with difficulties in pathogenicity prediction. This study attempted to improve the performance of existing in silico algorithms and find a common solution based on individualization strategy. We initiated the individualization with the epilepsy-related SCN1A variants by sub-regional stratification. SCN1A missense variants related to epilepsy were retrieved from mutation databases, and benign missense variants were collected from ExAC database. Predictions were performed by using 10 traditional tools with stepwise optimizations. Model predictive ability was evaluated using the five-fold cross-validations on variants of SCN1A, SCN2A, and KCNQ2. Additional validation was performed in SCN1A variants of damage-confirmed/familial epilepsy. The performance of commonly used predictors was less satisfactory for SCN1A with accuracy less than 80% and varied dramatically by functional domains of Nav1.1. Multistep individualized optimizations, including cutoff resetting, domain-based stratification, and combination of predicting algorithms, significantly increased predictive performance. Similar improvements were obtained for variants in SCN2A and KCNQ2. The predictive performance of the recently developed ensemble tools, such as Mendelian clinically applicable pathogenicity, combined annotation-dependent depletion and Eigen, was also improved dramatically by application of the strategy with molecular sub-regional stratification. The prediction scores of SCN1A variants showed linear correlations with the degree of functional defects and the severity of clinical phenotypes. This study highlights the need of individualized optimization with molecular sub-regional stratification for each gene in practice.

Synaptic clustering differences due to different GABRB3 mutations cause variable epilepsy syndromes

GABRB3 is highly expressed early in the developing brain, and its encoded β3 subunit is critical for GABAA receptor assembly and trafficking as well as stem cell differentiation in embryonic brain. To date, over 400 mutations or variants have been identified in GABRB3. Mutations in GABRB3 have been increasingly recognized as a major cause for severe paediatric epilepsy syndromes such as Lennox-Gastaut syndrome, Dravet syndrome and infantile spasms with intellectual disability as well as relatively mild epilepsy syndromes such as childhood absence epilepsy. There is no plausible molecular pathology for disease phenotypic heterogeneity. Here we used a very high-throughput flow cytometry assay to evaluate the impact of multiple human mutations in GABRB3 on receptor trafficking. In this study we found that surface expression of mutant β3 subunits is variable. However, it was consistent that surface expression of partnering γ2 subunits was lower when co-expressed with mutant than with wild-type subunits. Because γ2 subunits are critical for synaptic GABAA receptor clustering, this provides an important clue for understanding the pathophysiology of GABRB3 mutations. To validate our findings further, we obtained an in-depth comparison of two novel mutations [GABRB3 (N328D) and GABRB3 (E357K)] associated with epilepsy with different severities of epilepsy phenotype. GABRB3 (N328D) is associated with the relatively severe Lennox-Gastaut syndrome, and GABRB3 (E357K) is associated with the relatively mild juvenile absence epilepsy syndrome. With functional characterizations in both heterologous cells and rodent cortical neurons by patch-clamp recordings, confocal microscopy and immunoblotting, we found that both the GABRB3 (N328D) and GABRB3 (E357K) mutations reduced total subunit expression in neurons but not in HEK293T cells. Both mutant subunits, however, were reduced on the cell surface and in synapses, but the Lennox-Gastaut syndrome mutant β3 (N328D) subunit was more reduced than the juvenile absence epilepsy mutant β3 (E357K) subunit. Interestingly, both mutant β3 subunits impaired postsynaptic clustering of wild-type GABAA receptor γ2 subunits and prevented γ2 subunits from incorporating into GABAA receptors at synapses, although by different cellular mechanisms. Importantly, wild-type γ2 subunits were reduced and less clustered at inhibitory synapses in Gabrb3+/- knockout mice. This suggests that impaired receptor localization to synapses is a common pathophysiological mechanism for GABRB3 mutations, although the extent of impairment may be different among mutant subunits. The study thus identifies the novel mechanism of impaired targeting of receptors containing mutant β3 subunits and provides critical insights into understanding how GABRB3 mutations produce severe epilepsy syndromes and epilepsy phenotypic heterogeneity.

De Novo Germline Mutations in SEMA5A Associated With Infantile Spasms

Infantile spasm (IS) is an early-onset epileptic encephalopathy that usually presents with hypsarrhythmia on an electroencephalogram with developmental impairment or regression. In this study, whole-exome sequencing was performed to detect potential pathogenic de novo mutations, and finally we identified a novel damaging de novo mutation in SEMA5A and a compound heterozygous mutation in CLTCL1 in three sporadic trios with IS. The expression profiling of SEMA5A in the human brain showed that it was mainly highly expressed in the cerebral cortex, during the early brain development stage (8 to 9 post-conception weeks and 0 to 5 months after birth). In addition, we identified a close protein-protein interaction network between SEMA5A and candidate genes associated with epilepsy, autism spectrum disorder (ASD) or intellectual disability. Gene enrichment and function analysis demonstrated that genes interacting with SEMA5A were significantly enriched in several brain regions across early fetal development, including the cortex, cerebellum, striatum and thalamus (q < 0.05), and were involved in axonal, neuronal and synapse-associated processes. Furthermore, SEMA5A and its interacting genes were associated with ASD, epilepsy syndrome and developmental disorders of mental health. Our results provide insightful information indicating that SEMA5A may contribute to the development of the brain and is associated with IS. However, further genetic studies are still needed to evaluate the role of SEMA5A in IS to definitively establish the role of SEMA5A in this disorder.

Novel WWOX deleterious variants cause early infantile epileptic encephalopathy, severe developmental delay and dysmorphism among Yemenite Jews

The human WW Domain Containing Oxidoreductase (WWOX) gene was originally described as a tumor suppressor gene. However, recent reports have demonstrated its cardinal role in the pathogenesis of central nervous systems disorders such as epileptic encephalopathy, intellectual disability, and spinocerebellar ataxia. We report on six patients from three unrelated families of full or partial Yemenite Jewish ancestry exhibiting early infantile epileptic encephalopathy and profound developmental delay. Importantly, four patients demonstrated facial dysmorphism. Exome sequencing revealed that four of the patients were homozygous for a novel WWOX c.517-2A > G splice-site variant and two were compound heterozygous for this variant and a novel c.689A > C, p.Gln230Pro missense variant. Complementary DNA sequencing demonstrated that the WWOX c.517-2A > G splice-site variant causes skipping of exon six. A carrier rate of 1:177 was found among Yemenite Jews. We provide the first detailed description of patients harboring a splice-site variant in the WWOX gene and propose that the clinical synopsis of WWOX related epileptic encephalopathy should be broadened to include facial dysmorphism. The increased frequency of the c.517-2A > G splice-site variant among Yemenite Jews coupled with the severity of the phenotype makes it a candidate for inclusion in expanded preconception screening programs.

The landscape of early infantile epileptic encephalopathy in a consanguineous population

PURPOSE

Epileptic encephalopathies (EE), are a group of age-related disorders characterized by intractable seizures and electroencephalogram (EEG) abnormalities that may result in cognitive and motor delay. Early infantile epileptic encephalopathies (EIEE) manifest in the first year of life. EIEE are highly heterogeneous genetically but a genetic etiology is only identified in half of the cases, typically in the form of de novo dominant mutations.

METHOD

This is a descriptive retrospective study of a consecutive series of patients diagnosed with EIEE from the participating hospitals. A chart review was performed for all patients. The diagnosis of epileptic encephalopathy was confirmed by molecular investigations in commercial labs. In silico study was done for all novel mutations. A systematic search was done for all the types of EIEE and their correlated genes in the literature using the Online Mendelian Inheritance In Man and PubMed databases.

RESULTS

In this case series, we report 72 molecularly characterized EIEE from a highly consanguineous population, and review their clinical course. We identified 50 variants, 26 of which are novel, causing 26 different types of EIEE. Unlike outbred populations, autosomal recessive EIEE accounted for half the cases. The phenotypes ranged from self-limiting and drug-responsive to severe refractory seizures or even death.

CONCLUSIONS

We reported the largest EIEE case series in the region with confirmed molecular testing and detailed clinical phenotyping. The number autosomal recessive predominance could be explained by the society's high consanguinity. We reviewed all the EIEE registered causative genes in the literature and proposed a functional classification.

Identification of a novel homozygous variant confirms ITPA as a developmental and epileptic encephalopathy gene

ITPA related epileptic encephalopathy (epileptic encephalopathy, early infantile, 35) is a rare inborn error of metabolism. All reported individuals with this condition, including the present case manifest global developmental delay, seizures, progressive postnatal microcephaly, hypotonia, thin corpus callosum, cerebral atrophy, delayed myelination and white matter changes in posterior limb of internal capsule. Cataract and dilated cardiomyopathy are other characteristic findings. Currently, a single publication describes this condition in four families. Three truncating and two missense variants in ITPA have been identified in these families. We hereby report another family with ITPA related disorder and review the genotype and phenotype of the reported subjects.

Mutations in ACTL6B, coding for a subunit of the neuron-specific chromatin remodeling complex nBAF, cause early onset severe developmental and epileptic encephalopathy with brain hypomyelination and cerebellar atrophy

Developmental and epileptic encephalopathies (DEEs) are genetically heterogenous conditions, often characterized by early onset, EEG interictal epileptiform abnormalities, polymorphous and drug-resistant seizures, and neurodevelopmental impairments. In this study, we investigated the genetic defects in two siblings who presented with severe DEE, microcephaly, spastic tetraplegia, diffuse brain hypomyelination, cerebellar atrophy, short stature, and kyphoscoliosis. Whole exome next-generation sequencing (WES) identified in both siblings a homozygous non-sense variant in the ACTL6B gene (NM_016188:c.820C>T;p.Gln274*) coding for a subunit of the neuron-specific chromatin remodeling complex nBAF. To further support these findings, a targeted ACTL6B sequencing assay was performed on a cohort of 85 unrelated DEE individuals, leading to the identification of a homozygous missense variant (NM_016188:c.1045G>A;p.Gly349Ser) in a patient. This variant did not segregate in the unaffected siblings in this family and was classified as deleterious by several prediction softwares. Interestingly, in both families, homozygous patients shared a rather homogeneous phenotype. Very few patients with ACTL6B gene variants have been sporadically reported in WES cohort studies of patients with neurodevelopmental disorders and/or congenital brain malformations. However, the limited number of patients with incomplete clinical information yet reported in the literature did not allow to establish a strong gene-disease association. Here, we provide additional genetic and clinical data on three new cases that support the pathogenic role of ACTL6B gene mutation in a syndromic form of DEE.