ZFHX3 variants cause childhood partial epilepsy and infantile spasms with favourable outcomes

BACKGROUND

The ZFHX3 gene plays vital roles in embryonic development, cell proliferation, neuronal differentiation and neuronal death. This study aims to explore the relationship between ZFHX3 variants and epilepsy.

METHODS

Whole-exome sequencing was performed in a cohort of 378 patients with partial (focal) epilepsy. A Drosophila Zfh2 knockdown model was used to validate the association between ZFHX3 and epilepsy.

RESULTS

Compound heterozygous ZFHX3 variants were identified in eight unrelated cases. The burden of ZFHX3 variants was significantly higher in the case cohort, shown by multiple/specific statistical analyses. In Zfh2 knockdown flies, the incidence and duration of seizure-like behaviour were significantly greater than those in the controls. The Zfh2 knockdown flies exhibited more firing in excitatory neurons. All patients presented partial seizures. The five patients with variants in the C-terminus/N-terminus presented mild partial epilepsy. The other three patients included one who experienced frequent non-convulsive status epilepticus and two who had early spasms. These three patients had also neurodevelopmental abnormalities and were diagnosed as developmental epileptic encephalopathy (DEE), but achieved seizure-free after antiepileptic-drug treatment without adrenocorticotropic-hormone/steroids. The analyses of temporal expression (genetic dependent stages) indicated that ZFHX3 orthologous were highly expressed in the embryonic stage and decreased dramatically after birth.

CONCLUSION

ZFHX3 is a novel causative gene of childhood partial epilepsy and DEE. The patients of infantile spasms achieved seizure-free after treatment without adrenocorticotropic-hormone/steroids implies a significance of genetic diagnosis in precise treatment. The genetic dependent stage provided an insight into the underlying mechanism of the evolutional course of illness.

Reprint of: Recessive APC2 missense variants associated with epilepsies without neurodevelopmental disorders

OBJECTIVES

The APC2 gene, encoding adenomatous polyposis coli protein-2, is involved in cytoskeletal regulation in neurons responding to endogenous extracellular signals and plays an important role in brain development. Previously, the APC2 variants have been reported to be associated with cortical dysplasia and intellectual disability. This study aims to explore the association between APC2 variants and epilepsy.

METHODS

Whole-exome sequencing (WES) was performed in cases (trios) with epilepsies of unknown causes. The damaging effects of variants were predicted by protein modeling and in silico tools. Previously reported APC2 variants were reviewed to analyze the genotype-phenotype correlations.

RESULTS

Four pairs of compound heterozygous missense variants were identified in four unrelated patients with epilepsy without brain malformation/intellectual disability. All variants presented no or low allele frequencies in the controls. The missense variants were predicted to be damaging by silico tools, and affect hydrogen bonding with surrounding amino acids or decreased protein stability. Patients with variants that resulted in significant changes in protein stability exhibited more severe and intractable epilepsy, whereas patients with variants that had minor effect on protein stability exhibited relatively mild phenotypes. The previously reported APC2 variants in patients with complex cortical dysplasia with other brain malformations-10 (CDCBM10; MIM: 618677) were all truncating variants; in contrast, the variants identified in epilepsy in this study were all missense variants, suggesting a potential genotype-phenotype correlation.

SIGNIFICANCE

This study suggests that APC2 is potentially associated with epilepsy without brain malformation/intellectual disability. The genotype-phenotype correlation helps to understand the underlying mechanisms of phenotypic heterogeneity.

Epilepsy-associated genes: an update

PURPOSE

To provide an updated list of epilepsy-associated genes based on clinical-genetic evidence.

METHODS

Epilepsy-associated genes were systematically searched and cross-checked from the OMIM, HGMD, and PubMed databases up to July 2023. To facilitate the reference for the epilepsy-associated genes that are potentially common in clinical practice, the epilepsy-associated genes were ranked by the mutation number in the HGMD database and by case number in the China Epilepsy Gene 1.0 project, which targeted common epilepsy.

RESULTS

Based on the OMIM database, 1506 genes were identified to be associated with epilepsy and were classified into three categories according to their potential association with epilepsy or other abnormal phenotypes, including 168 epilepsy genes that were associated with epilepsies as pure or core symptoms, 364 genes that were associated with neurodevelopmental disorders as the main symptom and epilepsy, and 974 epilepsy-related genes that were associated with gross physical/systemic abnormalities accompanied by epilepsy/seizures. Among the epilepsy genes, 115 genes (68.5%) were associated with epileptic encephalopathy. After cross-checking with the HGMD and PubMed databases, an additional 1440 genes were listed as potential epilepsy-associated genes, of which 278 genes have been repeatedly identified variants in patients with epilepsy. The top 100 frequently reported/identified epilepsy-associated genes from the HGMD database and the China Epilepsy Gene 1.0 project were listed, among which 40 genes were identical in both sources.

SIGNIFICANCE

Recognition of epilepsy-associated genes will facilitate genetic screening strategies and be helpful for precise molecular diagnosis and treatment of epilepsy in clinical practice.

De novo GABRA1 variants in childhood epilepsies and the molecular subregional effects

Background

The GABRA1 gene, encoding the GABRAR subunit α1, plays vital roles in inhibitory neurons. Previously, the GABRA1 gene has been identified to be associated with developmental and epileptic encephalopathy (DEE) and idiopathic generalized epilepsy (IGE). This study aims to explore the phenotypic spectrum of GABRA1 and molecular subregional effect analysis.

Methods

Trios-based whole-exome sequencing was performed in patients with epilepsy. Previously reported GABRA1 mutations were systematically reviewed to analyze the molecular subregional effects.

Results

De novo GABRA1 mutations were identified in six unrelated patients with heterogeneous epilepsy, including three missense mutations (p.His83Asn, p.Val207Phe, and p.Arg214Cys) and one frameshift mutation (p.Thr453Hisfs*47). The two missense mutations, p.His83Asn and p.Val207Phe, were predicted to decrease the protein stability but no hydrogen bond alteration, with which the two patients also presented with mild genetic epilepsy with febrile seizures plus and achieved seizure-free status by monotherapy. The missense variant p.Arg214Cys was predicted to decrease protein stability and destroy hydrogen bonds with surrounding residues, which was recurrently identified in three cases with severe DEE. The frameshift variant p.Thr453Hisfs*47 was located in the last fifth residue of the C-terminus and caused an extension of 47 amino acids, with which the patients presented with moderated epilepsy with generalized tonic-clonic seizures alone (GTCA) but achieved seizure-free status by four drugs. The four variants were not presented in gnomAD and were evaluated as "pathogenic/likely pathogenic" according to ACMG criteria. Analysis of all reported cases indicated that patients with mutations in the N-terminal extracellular region presented a significantly higher percentage of FS and DEE, and the patients with variants in the transmembrane region presented earlier seizure onset ages.

Significance

This study suggested that GABRA1 variants were potentially associated with a spectrum of epilepsies, including EFS+, DEE, and GTCA. Phenotypic severity may be associated with the damaging effect of variants. The molecular subregional effects help in understanding the underlying mechanism of phenotypic variation.

Recessive APC2 missense variants associated with epilepsies without neurodevelopmental disorders

OBJECTIVES

The APC2 gene, encoding adenomatous polyposis coli protein-2, is involved in cytoskeletal regulation in neurons responding to endogenous extracellular signals and plays an important role in brain development. Previously, the APC2 variants have been reported to be associated with cortical dysplasia and intellectual disability. This study aims to explore the association between APC2 variants and epilepsy.

METHODS

Whole-exome sequencing (WES) was performed in cases (trios) with epilepsies of unknown causes. The damaging effects of variants were predicted by protein modeling and in silico tools. Previously reported APC2 variants were reviewed to analyze the genotype-phenotype correlations.

RESULTS

Four pairs of compound heterozygous missense variants were identified in four unrelated patients with epilepsy without brain malformation/intellectual disability. All variants presented no or low allele frequencies in the controls. The missense variants were predicted to be damaging by silico tools, and affect hydrogen bonding with surrounding amino acids or decreased protein stability. Patients with variants that resulted in significant changes in protein stability exhibited more severe and intractable epilepsy, whereas patients with variants that had minor effect on protein stability exhibited relatively mild phenotypes. The previously reported APC2 variants in patients with complex cortical dysplasia with other brain malformations-10 (CDCBM10; MIM: 618677) were all truncating variants; in contrast, the variants identified in epilepsy in this study were all missense variants, suggesting a potential genotype-phenotype correlation.

SIGNIFICANCE

This study suggests that APC2 is potentially associated with epilepsy without brain malformation/intellectual disability. The genotype-phenotype correlation helps to understand the underlying mechanisms of phenotypic heterogeneity.

HCFC1 variants in the proteolysis domain are associated with X-linked idiopathic partial epilepsy: Exploring the underlying mechanism

BACKGROUND

HCFC1 encodes transcriptional co-regulator HCF-1, which undergoes an unusual proteolytic maturation at a centrally located proteolysis domain. HCFC1 variants were associated with X-linked cobalamin metabolism disorders and mental retardation-3. This study aimed to explore the role of HCFC1 variants in common epilepsy and the mechanism underlying phenotype heterogeneity.

METHODS

Whole-exome sequencing was performed in a cohort of 313 patients with idiopathic partial (focal) epilepsy. Functional studies determined the effects of the variants on the proteolytic maturation of HCF-1, cell proliferation and MMACHC expression. The role of HCFC1 variants in partial epilepsy was validated in another cohort from multiple centers.

RESULTS

We identified seven hemizygous HCFC1 variants in 11 cases and confirmed the finding in the validation cohort with additional 13 cases and six more hemizygous variants. All patients showed partial epilepsies with favorable outcome. None of them had cobalamin disorders. Functional studies demonstrated that the variants in the proteolysis domain impaired the maturation by disrupting the cleavage process with loss of inhibition of cell growth but did not affect MMACHC expression that was associated with cobalamin disorder. The degree of functional impairment was correlated with the severity of phenotype. Further analysis demonstrated that variants within the proteolysis domain were associated with common and mild partial epilepsy, whereas those in the kelch domain were associated with cobalamin disorder featured by severe and even fatal epileptic encephalopathy, and those in the basic and acidic domains were associated with mainly intellectual disability.

CONCLUSION

HCFC1 is potentially a candidate gene for common partial epilepsy with distinct underlying mechanism of proteolysis dysfunction. The HCF-1 domains played distinct functional roles and were associated with different clinical phenotypes, suggesting a sub-molecular effect. The distinct difference between cobalamin disorders and idiopathic partial epilepsy in phenotype and pathogenic mechanism, implied a clinical significance in early diagnosis and management.

CELSR1 variants are associated with partial epilepsy of childhood

CELSR1 gene, encoding cadherin EGF LAG seven-pass G-type receptor 1, is mainly expressed in neural stem cells during the embryonic period. It plays an important role in neurodevelopment. However, the relationship between CELSR1 and disease of the central nervous system has not been defined. In this study, we performed trios-based whole-exome sequencing in a cohort of 356 unrelated cases with partial epilepsy without acquired causes and identified CELSR1 variants in six unrelated cases. The variants included one de novo heterozygous nonsense variant, one de novo heterozygous missense variant, and four compound heterozygous missense variants that had one variant was located in the extracellular region and the other in the cytoplasm. The patients with biallelic variants presented severe epileptic phenotypes, whereas those with heterozygous variants were associated with a mild epileptic phenotype of benign epilepsy with centrotemporal spikes (BECTS). These variants had no or low allele frequency in the gnomAD database. The frequencies of the CELSR1 variants in this cohort were significantly higher than those in the control populations. The evidence from ClinGen Clinical-Validity Framework suggested a strong association between CELSR1 variants and epilepsy. These findings provide evidence that CELSR1 is potentially a candidate pathogenic gene of partial epilepsy of childhood.

PRRT2 gene mutations associated with infantile convulsions induced by sucking and the genotype-phenotype correlation

Introduction

PRRT2 is a major causative gene for self-limited familial neonatal-infantile epilepsy, paroxysmal kinesigenic dyskinesia, and paroxysmal kinesigenic dyskinesia with infantile convulsions. Voluntary movement trigger is prominent in adolescence and adulthood, but the triggers are unknown in infants.

Methods

A gene panel designed for targeted next-generation sequencing (NGS) was used to screen genetic abnormalities in a cohort of 45 cases with infantile convulsions. The copy number variation was detected by a computational method based on the normalized depth of coverage and validated by a quantitative real-time polymerase chain reaction (RT-qPCR) method. The genotype-phenotype correlation of the PRRT2 mutation gene was analyzed.

Results

A de novo heterozygous PRRT2 deletion was identified in a child who had infantile convulsions induced by vigorous sucking. Seizures happened during the change of feeding behavior from breast to formula, which led to hungry and vigorous sucking. Ictal electroencephalograms recorded seizures with focal origination, which provided direct evidence of epileptic seizures in infants with PRRT2 mutations. Seizures stopped soon after the feeding behavior was changed by reducing feeding interval time and extending feeding duration. Data reanalysis on our previously reported cases with PRRT2 mutations showed that six of 18 (33.3%) patients had infantile convulsions or infantile non-convulsion seizures during feeding. The mutations included two truncating mutations (c.579dupA/p.Glu194Argfs^*6, and c.649dupC/p.Arg217Profs^*8) that were identified in each of the three affected individuals.

Conclusions

This study suggests that feeding, especially vigorous sucking, is potentially a trigger and highlights the significance of feeding behavior in preventing seizures in infants with PRRT2 mutations. Identification of PRRT2 haploinsufficiency mutations in the patients with infantile convulsions induced by sucking suggested a potential genotype-phenotype correlation.

Recessive LAMA5 Variants Associated With Partial Epilepsy and Spasms in Infancy

Objective

The LAMA5 gene encodes the laminin subunit α5, the most abundant laminin α subunit in the human brain. It forms heterotrimers with the subunit β1/β2 and γ1/γ3 and regulates neurodevelopmental processes. Genes encoding subunits of the laminin heterotrimers containing subunit α5 have been reported to be associated with human diseases. Among LAMAs encoding the laminin α subunit, LAMA1-4 have also been reported to be associated with human disease. In this study, we investigated the association between LAMA5 and epilepsy.

Methods

Trios-based whole-exome sequencing was performed in a cohort of 118 infants suffering from focal seizures with or without spasms. Protein modeling was used to assess the damaging effects of variations. The LAMAs expression was analyzed with data from the GTEX and VarCards databases.

Results

Six pairs of compound heterozygous missense variants in LAMA5 were identified in six unrelated patients. All affected individuals suffered from focal seizures with mild developmental delay, and three patients presented also spasms. These variants had no or low allele frequencies in controls and presented statistically higher frequency in the case cohort than in controls. The recessive burden analysis showed that recessive LAMA5 variants identified in this cohort were significantly more than the expected number in the East Asian population. Protein modeling showed that at least one variant in each pair of biallelic variants affected hydrogen bonds with surrounding amino acids. Among the biallelic variants in cases with only focal seizures, two variants of each pair were located in different structural domains or domains/links, whereas in the cases with spasms, the biallelic variants were constituted by two variants in the identical functional domains or both with hydrogen bond changes.

Conclusion

Recessive LAMA5 variants were potentially associated with infant epilepsy. The establishment of the association between LAMA5 and epilepsy will facilitate the genetic diagnosis and management in patients with infant epilepsy.

SHROOM4 Variants Are Associated With X-Linked Epilepsy With Features of Generalized Seizures or Generalized Discharges

Objective

SHROOM4 gene encodes an actin-binding proteins, which plays an important role in cytoskeletal architecture, synaptogenesis, and maintaining gamma-aminobutyric acid receptors-mediated inhibition. SHROOM4 mutations were reported in patients with the Stocco dos Santos type of X-linked syndromic intellectual developmental disorder (SDSX; OMIM# 300434). In this study, we investigated the association between SHROOM4 and epilepsy.

Methods

Trios-based whole-exome sequencing was performed in a cohort of 320 cases with idiopathic generalized epilepsy or idiopathic partial epilepsy. Protein modeling was used to assess the damaging effects of variations.

Results

Six hemizygous missense SHROOM4 variants, including c.13C > A/p. Pro5Thr, c.3236C > T/p.Glu1079Ala, c.3581C > T/p.Ser1194Leu, c.4288C > T/p.Arg1430Cys, c.4303G > A/p.Val1435Met, c.4331C > T/p.Pro1444Leu, were identified in six cases with idiopathic epilepsy without intellectual disability. All patients presented with features of generalized seizures or generalized discharges. These hemizygous variants had no or extremely low allele frequencies in controls and showed statistically higher frequency in the case cohort than controls. All variants were predicted to alter hydrogen bond with surrounding amino acids or decreased protein stability. The SHROOM4 variants reported in patients with SDSX were mostly destructive or duplicative variants; in contrast, the SHROOM4 variants were all missense variants, suggesting a potential genotype-phenotype correlation. The two missense variants associated with SDSX were located in the middle of SHROOM4 protein, whereas variants associated with idiopathic epilepsy were located around the N-terminal PDZ domain and the C-terminal ASD2 domain.

Significance

SHROOM4 was potentially a candidate pathogenic gene of idiopathic epilepsy without intellectual disability. The genotype-phenotype correlation and sub-regional effect helps understanding the mechanism underlying phenotypic variation.

CACNA1A Mutations Associated With Epilepsies and Their Molecular Sub-Regional Implications

PURPOSE

Previously, mutations in the voltage-gated calcium channel subunit alpha1 A (CACNA1A) gene have been reported to be associated with paroxysmal disorders, typically as episodic ataxia type 2. To determine the relationship between CACNA1A and epilepsies and the role of molecular sub-regional on the phenotypic heterogeneity.

METHODS

Trio-based whole-exome sequencing was performed in 318 cases with partial epilepsy and 150 cases with generalized epilepsy. We then reviewed all previously reported CACNA1A mutations and analyzed the genotype-phenotype correlations with molecular sub-regional implications.

RESULTS

We identified 12 CACNA1A mutations in ten unrelated cases of epilepsy, including four de novo null mutations (c.2963_2964insG/p.Gly989Argfs*78, c.3089 + 1G > A, c.4755 + 1G > T, and c.6340-1G > A), four de novo missense mutations (c.203G > T/p.Arg68Leu, c.3965G > A/p.Gly1322Glu, c.5032C > T/p.Arg1678Cys, and c.5393C > T/p.Ser1798Leu), and two pairs of compound heterozygous missense mutations (c.4891A > G/p.Ile1631Val& c.5978C > T/p.Pro1993Leu and c.3233C > T/p.Ser1078Leu&c.6061G > A/p.Glu2021Lys). The eight de novo mutations were evaluated as pathogenic or likely pathogenic mutations according to the criteria of American College of Medical Genetics and Genomics (ACMG). The frequencies of the compound heterozygous CACNA1A mutations identified in this cohort were significantly higher than that in the controls of East Asian and all populations (P = 7.30 × 10^-4, P = 2.53 × 10^-4). All of the ten cases were ultimately seizure-free after antiepileptic treatment, although frequent epileptic seizures were observed in four cases. Further analysis revealed that episodic ataxia type 2 (EA2) had a tendency of higher frequency of null mutations than epilepsies. The missense mutations in severe epileptic phenotypes were more frequently located in the pore region than those in milder epileptic phenotypes (P = 1.67 × 10^-4); de novo mutations in the epilepsy with intellectual disability (ID) had a higher percentage than those in the epilepsy without ID (P = 1.92 × 10^-3).

CONCLUSION

This study suggested that CACNA1A mutations were potentially associated with pure epilepsy and the spectrum of epileptic phenotypes potentially ranged from the mild form of epilepsies such as absence epilepsy or partial epilepsy, to the severe form of developmental epileptic encephalopathy. The clinical phenotypes variability is potentially associated with the molecular sub-regional of the mutations.

CELSR3 variants are associated with febrile seizures and epilepsy with antecedent febrile seizures

Aims

To identify novel pathogenic gene of febrile seizures (FS)/epilepsy with antecedent FS (EFS+).

Methods

The trio‐based whole‐exome sequencing was performed in a cohort of 462 cases with FS/EFS+. Silico programs, sequence alignment, and protein modeling were used to predict the damaging of variants. Statistical testing was performed to analyze gene‐based burden of variants.

Results

Five heterozygous missense variants in CELSR3 were detected in five cases (families) with eight individuals (five females, three males) affected. Two variants were de novo, and three were identified in families with more than one individual affected. All the variants were predicted to be damaging in silico tools. Protein modeling showed that the variants resulted in disappearance of multiple hydrogen bonds and one disulfide bond, which potentially caused functional impairments of protein. The frequency of CELSR3 variants identified in this study was significantly higher than that in controls. All affected individuals were diagnosed with FS/EFS+, including six patients with FS and two patients with EFS+. All cases presented favorable outcomes without neurodevelopmental disorders.

Conclusions

CELSR3 variants are potentially associated with FS/EFS+.

[Correlation study between changes in intestinal microflora structure and immune indexes in newly treated patients with pulmonary tuberculosis]

To explore the correlation between the changes of the intestinal flora of newly treated pulmonary tuberculosis patients and the immune indicators of the body, and to provide a reference for the prevention and treatment of pulmonary tuberculosis. A single-center and case-control study was adopted. From October 2020 to April 2021, 43 patients with newly diagnosed tuberculosis in the Department of Tuberculosis, Affiliated Changsha Central Hospital,University of South China were selected as the control group. 43 cases of newly treated pulmonary tuberculosis (PTB), 43 healthy control (HC) during the same period, collected fresh feces and whole blood of subjects, and used Illumina Hiseq high-throughput sequencing technology to analyze 16S of all microorganisms in feces The V4 region of rRNA was amplified and sequenced, and the structure of the intestinal flora was analyzed by QIIME software. Use flow cytometry to determine the subject's immune indicators (CD3⁺, CD4⁺, CD8⁺, CD4⁺CD25⁺CD127^-Treg, CD14⁺CD16⁺, CD14⁺CD16^-), and analyze the changes in intestinal flora and immune function in newly treated pulmonary tuberculosis patients Inherent connection. The χ² test, t test, and Wilcox rank sum test were used to analyze the differences in age, gender, α diversity, and relative abundance of the two groups of people. Compared with the HC group, the alpha diversity of the intestinal flora in the PTB group decreased (shannon index: t=3.906, P=0.000 2; simpson index: Z=553, P=0.004 7; chao1 index: t=5.395, P=0.000 0). β diversity analysis showed that there were significant differences in the structure of the intestinal flora between the two groups (P=0.000). Species difference analysis showed that at the phylum level, the relative abundance of Firmicutes in the PTB group was significantly lower than that in the HC group (Z=486.0, P=0.000 5). At the genus level, there are 15 different bacterial genera between the two groups. In the PTB group, bifidobacterium, enterococcus, lactobacillus, anaerostipes, the relative abundance of the above 5 genera of veillonella is higher than that of the HC group (P<0.05); Butyricimonas, clostridium, and broutella (blautia), coprococcus, dorea, lachnospira, roseburia, faecalibacterium, ruminococcus, the relative abundance of 10 bacterial genera including dialister was lower than that of the HC group (P<0.05). Comparison of immune indexes between groups showed that CD14⁺CD16⁺monocytes (%) in the PTB group were higher than those in the HC group (t=2.456, P=0.001 6<0.05), while CD14⁺CD16^-monocytes (%) were lower than HC (t=-4.368, P=0.000<0.05), while the differences in CD3⁺, CD4⁺, CD8⁺, CD4⁺/CD8⁺and Treg (CD4⁺CD25⁺CD127^-) were not statistically significant (P>0.05). Spearman correlation analysis showed that Firmicutes in the PTB group was negatively correlated with CD4⁺/CD8⁺, CD14⁺CD16⁺(r=-0.218, P=0.048; r=-0.245, P=0.025), and positively correlated with CD14⁺CD16^-Correlation (r=0.250, P=0.022); At the genus level, Faecalis is positively correlated with CD4⁺/CD8⁺and CD4⁺(r=0.250, P=0.023; r=0.258, P=0.019); Rosella and CD3⁺, CD8⁺and CD14⁺CD16^-are positively correlated (r=0.27, P=0.024; r=0.219, P=0.046; r=0.027, P=0.039), and negatively correlated with CD14⁺CD16⁺(r=-0.280, P= 0.01). Changes in the structure of the intestinal flora of newly treated pulmonary tuberculosis patients may be one of the influencing factors of the immune function of the body. Targeted optimization of the structure of the intestinal flora and improvement of the body's immunity may be used as an effective auxiliary treatment for pulmonary tuberculosis.

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.

CHD4 variants are associated with childhood idiopathic epilepsy with sinus arrhythmia

Aims

CHD4 gene, encoding chromodomain helicase DNA‐binding protein 4, is a vital gene for fetal development. In this study, we aimed to explore the association between CHD4 variants and idiopathic epilepsy.

Methods

Trios‐based whole‐exome sequencing was performed in a cohort of 482 patients with childhood idiopathic epilepsy. The Clinical Validity Framework of ClinGen and an evaluating method from five clinical‐genetic aspects were used to determine the association between CHD4 variants and epilepsy.

Results

Four novel heterozygous missense mutations in CHD4, including two de novo mutations (c.1597A>G/p.K533E and c.4936G>A/p.E1646K) and two inherited mutations with co‐segregation (c.856C>G/p.P286A and c.4977C>G/p.D1659E), were identified in four unrelated families with eight individuals affected. Seven affected individuals had sinus arrhythmia. From the molecular sub‐regional point of view, the missense mutations located in the central regions from SNF2‐like region to DUF1087 domain were associated with multisystem developmental disorders, while idiopathic epilepsy‐related mutations were outside this region. Strong evidence from ClinGen Clinical Validity Framework and evidences from four of the five clinical‐genetic aspects suggested an association between CHD4 variants and epilepsy.

Conclusions

CHD4 was potentially a candidate pathogenic gene of childhood idiopathic epilepsy with arrhythmia. The molecular sub‐regional effect of CHD4 mutations helped explaining the mechanisms underlying phenotypic variations.

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.

RYR2 Mutations Are Associated With Benign Epilepsy of Childhood With Centrotemporal Spikes With or Without Arrhythmia

RYR2 encodes ryanodine receptor 2 protein (RYR-2) that is mainly located on endoplasmic reticulum membrane and regulates intracellular calcium concentration. The RYR-2 protein is ubiquitously distributed and highly expressed in the heart and brain. Previous studies have identified the RYR2 mutations in the etiology of arrhythmogenic right ventricular dysplasia 2 and catecholaminergic polymorphic ventricular tachycardia. However, the relationship between RYR2 gene and epilepsy is not determined. In this study, we screened for novel genetic variants in a group of 292 cases (families) with benign epilepsy of childhood with centrotemporal spikes (BECTS) by trio-based whole-exome sequencing. RYR2 mutations were identified in five cases with BECTS, including one heterozygous frameshift mutation (c.14361dup/p.Arg4790Pro fs^∗6), two heterozygous missense mutations (c.2353G > A/p.Asp785Asn and c.8574G > A/p.Met2858Ile), and two pairs of compound heterozygous mutations (c.4652A > G/p.Asn1551Ser and c.11693T > C/p.Ile3898Thr, c.7469T > C/p.Val2490Ala and c.12770G > A/p.Arg4257Gln, respectively). Asp785Asn was a de novo missense mutation. All the missense mutations were suggested to be damaging by at least three web-based prediction tools. These mutations do not present or at low minor allele frequency in gnomAD database and present statistically higher frequency in the cohort of BECTS than in the control populations of gnomAD. Asp785Asn, Asn1551Ser, and Ile3898Thr were predicted to affect hydrogen bonds with surrounding amino acids. Three affected individuals had arrhythmia (sinus arrhythmia and occasional atrial premature). The two probands with compound heterozygous missense mutations presented mild cardiac structural abnormalities. Strong evidence from ClinGen Clinical Validity Framework suggested an association between RYR2 variants and epilepsy. This study suggests that RYR2 gene is potentially a candidate pathogenic gene of BECTS. More attention should be paid to epilepsy patients with RYR2 mutations, which were associated with arrhythmia and sudden unexpected death in previous reports.

HLA Risk Alleles in Aromatic Antiepileptic Drug-Induced Maculopapular Exanthema

To characterize human leukocyte antigen (HLA) loci as risk factors in aromatic antiepileptic drug-induced maculopapular exanthema (AED-MPE). A case-control study was performed to investigate HLA loci involved in AED-MPE in a southern Han Chinese population. Between January 2007 and June 2019, 267 patients with carbamazepine (CBZ), oxcarbazepine (OXC), or lamotrigine (LTG) associated MPE and 387 matched drug-tolerant controls from six centers were enrolled. HLA-A/B/C/DRB1 genotypes were determined using sequence-based typing. Potential risk alleles were validated by meta-analysis using data from different populations and in silico analysis of protein-drug interactions. HLA-DRB1*04:06 was significantly associated with OXC-MPE (p = 0.002, p _c = 0.04). HLA-B*38:02 was associated with CBZ-MPE (p = 0.03). When pooled, HLA-A*24:02, HLA-A*30:01, and HLA-B*35:01 additionally revealed significant association with AED-MPE. Logistic regression analysis showed a multiplicative interaction between HLA-A*24:02 and HLA-B*38:02 in CBZ-MPE. Meta-analysis of data from different populations revealed that HLA-24*:02 and HLA-A*30:01 were associated with AED-MPE (p = 0.02 and p = 0.04, respectively). In silico analysis of protein-drug interaction demonstrated that HLA-A*24:02 and HLA-A*30:01 had higher affinities with the three aromatic AEDs than the risk-free HLA-A allele. HLA-DRB1*04:06 showed relatively specific high affinity with S-monohydroxy derivative of OXC. HLA-DRB1*04:06 is a specific risk allele for OXC-induced MPE in the Southern Han Chinese. HLA-A*24:02, possibly HLA-A*30:01, are common risk factors for AED-MPE. The multiplicative risk potential between HLA-A*24:02 and HLA-B*38:02 suggests that patients with two risk alleles are at greater risk than those with one risk allele. Inclusion of these HLA alleles in pre-treatment screening would help estimating the risk of AED-MPE.

Heterozygous PGM3 Variants Are Associated With Idiopathic Focal Epilepsy With Incomplete Penetrance

Introduction

Idiopathic focal epilepsy (IFE) is a group of self-limited epilepsies. The etiology for the majority of the patients with IFE remains elusive. We thus screened disease-causing variants in the patients with IFE.

Methods

Whole-exome sequencing was performed in a cohort of 323 patients with IFE. Protein modeling was performed to predict the effects of missense variants. The genotype-phenotype correlation of the newly defined causative gene was analyzed.

Results

Four novel heterozygous variants in PGM3, including two de novo variants, were identified in four unrelated individuals with IFE. The variants included one truncating variant (c.1432C > T/p.Q478X) and three missense variants (c.478C > T/p.P160S, c.1239C > G/p.N413K, and c.1659T > A/p.N553K), which had no allele frequency in the gnomAD database. The missense variants were predicted to be damaging and affect hydrogen bonds with surrounding amino acids. Mutations Q478X, P160S, and N413K were associated with benign childhood epilepsy with centrotemporal electroencephalograph (EEG) spikes. P160S and N413K were located in the inner side of the enzyme active center. Mutation N553K was associated with benign occipital epilepsy with incomplete penetrance, located in the C-terminal of Domain 4. Further analysis demonstrated that previously reported biallelic PGM3 mutations were associated with severe immunodeficiency and/or congenital disorder of glycosylation, commonly accompanied by neurodevelopmental abnormalities, while monoallelic mutations were associated with milder symptoms like IFE.

Conclusion

The genetic and molecular evidence from the present study implies that the PGM3 variants identified in IFE patients lead to defects of the PGM3 gene, suggesting that the PGM3 gene is potentially associated with epilepsy. The genotype-phenotype relationship of PGM3 mutations suggested a quantitative correlation between genetic impairment and phenotypic severity, which helps explain the mild symptoms and incomplete penetrance in individuals with IFE.

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.

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.

HLA-B*13:01 as a Risk Allele for Antiepileptic Drugs-Induced Cutaneous Adverse Reactions: Higher Risk for Cross-Reactivity?

Antiepileptic drugs frequently cause cutaneous adverse reactions (cADRs). Numerous studies have reported associations between human leukocyte antigen (HLA) alleles and cADRs caused by single antiepileptic drug in Southern Han Chinese people. However, the relationship between the HLA allele and cADRs sequentially induced by two or more antiepileptic drugs (AEDs-induced cross-reactivity) is unclear. To explore the associations between HLA alleles and AEDs-induced cross-reactivity, we prospectively recruited patients with AEDs-induced cross-reactivity from 2009 to 2017 and performed high-resolution genotyping to detect the HLA-A, B, C, and DRB1 alleles in patients for comparison with normal controls. To verify the important genotype, we compared its presence in patients with cross-reactivity to enlarged normal controls, and its presence in patients with carbamazepine (CBZ)-induced maculopapular exanthema (MPE) to CBZ-tolerant controls. Further, the important allele was replicated by meta-analysis. Twenty-three patients with AED-induced cross-reactivity and 500 healthy individuals were enrolled from Southern China. All patients had a mild rash without mucosal or systemic involvement. The HLA-B^*13:01 allele was present in 34.78% (8/23) of patients, 14.60% (73/500) of healthy individuals, and 14.5% (763/5,270) healthy individuals, revealing a significant association (8/23 vs. 73/500; P = 0.02; OR: 3.12; 95% CI: 1.28–7.62; 8/23 vs. 763/5,270; P = 0.014; OR: 3.15; 95% CI: 1.33–7.46). HLA-B^*13:01 was presented numerically higher in CBZ-induced MPE than that in CBZ-tolerant individuals without statistical significance (33/145, 22.76%, vs. 28/179, 15.64%; P = 0.103). Meta-analysis revealed an association between HLA-B^*13:01 and cADRs induced by single AEDs or/and non-AEDs in Chinese and Thai populations (P = 0.000). This study suggests that HLA-B^*13:01 is potentially associated with AED-cADRs in general, possibly with stronger effect in cross-reactivity. Screening for HLA-B^*13:01 prior to starting AEDs therapy may help to avoid cADRs. However, this association requires further analysis in a multi-center study with a larger sample size.

[Changes of Rheb gene and protein expression in preeclampsia-like mouse model treated with pravastatin]

Objective: To explore whether pravastatin (Pra) inhibits mammalian target of rapamycin (mTOR) signal pathway by regulating Ras homolog enriched in brain (Rheb) protein through the comparison of gene and protein expression changes of Rheb in liver and placenta in preeclampsia (PE)-like mouse model treated with Pra. Methods: C57BL/6J pregnant mice were randomly divided into two groups. The PE group was established by injecting N-nitro-L-arginine methyl ester (L-NAME) daily at gestational 7-18 days, saline was injected as contol group (Con); then giving mice Pra (PE+Pra, Con+Pra group, n=8) or normal saline (PE+N, Con+N group, n=8) every day from the 8th gestational day of pregnancy. The maternal liver and placenta tissues were collected on the 18th day of pregnancy. Western blot, real-time quantitative PCR and immunohistochemistry were used to compare the levels of Rheb protein and mRNA expression in the liver and placenta. Results: (1)The results of western blot: there were no significant differences in Rheb protein expression between PE+N group (liver: 0.706±0.123; placenta: 0.866±0.128) and Con+N group (liver: 0.732±0.123; placenta: 0.909±0.097) , and the differences between PE+Pra group (liver: 0.669±0.134; placenta: 0.940±0.221) and PE+N group were not significant either in liver or in placenta (all P>0.05). (2) The results of real-time quantitative PCR: when PE+N group (liver: 1.026±0.480; placenta: 1.102±0.361) compared with Con+N group (liver: 1.058±0.389; placenta: 1.067±0.400) , PE+Pra group (liver: 0.735±0.356; placenta: 0.822±0.304) compared with PE+N group, there were no significant differences either in liver or in placenta (all P>0.05). (3) The results of immunohistochemistry: Rheb protein expression did not change significantly in maternal liver and placenta, there were no significant differences in protein expression levels between PE+N group and Con+N group, and between PE+Pra group and PE+N group (all P>0.05). Conclusion: The inhibition of Pra on mTOR signaling pathway in some PE-like model may be independent of the expression of Rheb gene and protein.

Familial paroxysmal kinesigenic dyskinesia is associated with mutations in the KCNA1 gene

Paroxysmal kinesigenic dyskinesia (PKD) is a heterogeneous movement disorder characterized by recurrent dyskinesia attacks triggered by sudden movement. PRRT2 has been identified as the first causative gene of PKD. However, it is only responsible for approximately half of affected individuals, indicating that other loci are most likely involved in the etiology of this disorder. To explore the underlying causative gene of PRRT2-negative PKD, we used a combination strategy including linkage analysis, whole-exome sequencing and copy number variations analysis to detect the genetic variants within a family with PKD. We identified a linkage locus on chromosome 12 (12p13.32-12p12.3) and detected a novel heterozygous mutation c.956 T>G (p.319 L>R) in the potassium voltage-gated channel subfamily A member 1, KCNA1. Whole-exome sequencing in another 58 Chinese patients with PKD who lacked mutations in PRRT2 revealed another novel mutation in the KCNA1 gene [c.765 C>A (p.255 N>K)] within another family. Biochemical analysis revealed that the L319R mutant accelerated protein degradation via the proteasome pathway and disrupted membrane expression of the Kv1.1 channel. Electrophysiological examinations in transfected HEK293 cells showed that both the L319R and N255K mutants resulted in reduced potassium currents and respective altered gating properties, with a dominant negative effect on the Kv1.1 wild-type channel. Our study suggests that these mutations in KCNA1 cause the Kv1.1 channel dysfunction, which leads to familial PKD. The current study further extended the genotypic spectrum of this disorder, indicating that Kv1.1 channel dysfunction maybe one of the underlying defects in PKD.

Homozygous missense TPP1 mutation associated with mild late infantile neuronal ceroid lipofuscinosis and the genotype-phenotype correlation

PURPOSE

TPP1 mutations have been identified in patients with variable phenotypes such as late infantile neuronal ceroid lipofuscinosis (LINCL), juvenile neuronal ceroid lipofuscinosis (JNCL), and spinocerebellar ataxia 7. However, the mechanism underlying phenotype variation is unknown. We screened TPP1 mutations in patients with epilepsies and analyzed the genotype-phenotype correlation to explain the phenotypic variations.

METHODS

We performed targeted next-generation sequencing in a cohort of 330 patients with epilepsies. All previously reported TPP1 mutations were systematically retrieved from the PubMed and NCL Mutation Database.

RESULTS

The homozygous missense TPP1 mutation c.646 G > A/ p.Val216Met was identified in a family with two affected siblings. The proband presented with seizures from three years of age, while no ataxia, cognitive regression, or visual abnormalities were observed. Further analysis of all reported TPP1 mutations revealed that the LINCL group had a significantly higher frequency of truncating and invariant splice-site mutations than the JNCL group. In contrast, the JNCL group had a higher frequency of variant splice-site mutations than LINCL. There was a significant correlation between phenotype severity and the frequency of destructive mutation.

CONCLUSION

This study suggested that the phenotype of mainly epilepsy can be included in the phenotypic spectrum of TPP1 mutations, which are candidate targets for genetic screening in patients with epilepsy. With the development of therapy techniques, early genetic diagnosis may enable the improvement of etiology-targeted treatments. The relationship between phenotype severity and the genotype of TPP1 mutations may help explain the phenotypic variations.

Development-related aberrations in Kv1.1 α-subunit exert disruptive effects on bioelectrical activities of neurons in a mouse model of fragile X syndrome

Kv1.1, a Shaker homologue potassium channel, plays a critical role in homeostatic regulation of neuronal excitability. Aberrations in the functional properties of Kv1.1 have been implicated in several neurological disorders featured by neuronal hyperexcitability. Fragile X syndrome (FXS), the most common form of inherited mental retardation, is characterized by hyperexcitability in neural network and intrinsic membrane properties. The Kv1.1 channel provides an intriguing mechanistic candidate for FXS. We investigated the development-related expression pattern of the Kv1.1 α-subunit by using a Fmr1 knockout (KO) mouse model of FXS. Markedly decreased protein expression of Kv1.1 was found in neonatal and adult stages when compared to age-matched wild-type (WT) mice. Immunohistochemical investigations supported the delayed development-related increases in Kv1.1 expression, especially in CA3 pyramidal neurons. By applying a Kv1.1-specific blocker, dendrotoxin-κ (DTX-κ), we isolated the Kv1.1-mediated currents in the CA3 pyramidal neurons. The isolated DTX-κ-sensitive current of neurons from KO mice exhibited decreased amplitude, lower threshold of activation, and faster recovery from inactivation. The equivalent reduction in potassium current in the WT neurons following application of the appropriate amount of DTX-κ reproduced the enhanced firing abilities of KO neurons, suggesting the Kv1.1 channel as a critical contributor to the hyperexcitability of KO neurons. The role of Kv1.1 in controlling neuronal discharges was further supported by the parallel developmental trajectories of Kv1.1 expression, current amplitude, and discharge impacts, with a significant correlation between the amplitude of Kv1.1-mediated currents and Kv1.1-blocking-induced firing enhancement. These data suggest that the expression of the Kv1.1 α-subunit has a profound pathological relevance to hyperexcitability in FXS, as well as implications for normal development, maintenance, and control of neuronal activities.

[Regulation of pravastatin on long-chain fatty acid oxidative enzyme in pre-eclampsia-like mouse model]

Objective: To investigate the modulation of long-chain 3-hydroxyacyl-CoA dehydrogenase (LCHAD) expression by pravastatin in pre-eclampsia-like mouse model. Methods: C57BL/6J mice were randomly injected with N-nitro-L-arginine methyl ester (L-NAME) as pre-eclampsia-like model group (PE) or saline as normal pregnancy control group (Con) respectively, from gestational the 7th to 18th day. For each group, pravastatin (PE+Pra, Con+Pra group) or saline (PE+N, Con+N Group) was given from the 8th to 18th day of gestation, respectively. Liver and placenta of pregnant mice were collected on gestational day 18. The LCHAD protein expression and mRNA levels of liver and placenta were detected through western blot, immunohistochemistry and real-time quantitative PCR. Results: (1) The average arterial pressure of pregnant mice increased gradually from the 8th to 18th day in PE+N group, but decreased in PE+Pra group from gestational 10th day, 24 hour urinary protein levels in PE+N group [(1 494 ± 201) μg] were significantly higher than that in Con+N group [(935±128) μg, P<0.01], and also higher than that in PE+Pra group [(981±116) μg, P<0.01].(2) The results of western blot: the expression of LCHAD was significantly lower in PE+N group (liver: 0.64±0.11, placenta: 0.48±0.06) than that in Con+N group (liver: 1.06±0.10, placenta: 0.60±0.10), and lower than that in PE+Pra group (liver: 0.99±0.04, placenta: 0.60±0.08; all P<0.01).(3)The results of real-time quantitative PCR: the levels of LCHAD mRNA in liver and placenta in PE+N group (liver: 0.621±0.128, placenta: 0.646±0.129) were significantly decreased compared with Con+N group (liver: 1.007±0.130, placenta: 1.004±0.103; all P<0.01), but there was no significant difference between PE+Pra group (liver: 0.693±0.678, placenta: 0.662±0.183; P>0.05). (4) LCHAD protein was expressed widely and evenly in liver. The expression in placental cytotrophoblast and syncytial trophoblast cells located in outer layer of villous in labyrinth layer was the most. The expression of LCHAD was significantly lower in PE+N group (liver: 0.062±0.016, placenta: 0.147±0.018) than that in Con+N group (liver: 0.126±0.013, placenta: 0.183±0.024), and lower than that in PE+Pra group (liver: 0.111±0.017, placenta: 0.174±0.027; all P<0.05). Conclusion: Pravastatin could upregulate the LCHAD protein expression of liver and placenta in the pre-eclampsia-like mouse, which may be a mechanism to improve the clinical manifestations of pre-eclampsia.

ARHGEF9 mutations in epileptic encephalopathy/intellectual disability: toward understanding the mechanism underlying phenotypic variation

ARHGEF9 resides on Xq11.1 and encodes collybistin, which is crucial in gephyrin clustering and GABA_A receptor localization. ARHGEF9 mutations have been identified in patients with heterogeneous phenotypes, including epilepsy of variable severity and intellectual disability. However, the mechanism underlying phenotype variation is unknown. Using next-generation sequencing, we identified a novel mutation, c.868C > T/p.R290C, which co-segregated with epileptic encephalopathy, and validated its association with epileptic encephalopathy. Further analysis revealed that all ARHGEF9 mutations were associated with intellectual disability, suggesting its critical role in psychomotor development. Three missense mutations in the PH domain were not associated with epilepsy, suggesting that the co-occurrence of epilepsy depends on the affected functional domains. Missense mutations with severe molecular alteration in the DH domain, or located in the DH-gephyrin binding region, or adjacent to the SH3-NL2 binding site were associated with severe epilepsy, implying that the clinical severity was potentially determined by alteration of molecular structure and location of mutations. Male patients with ARHGEF9 mutations presented more severe phenotypes than female patients, which suggests a gene-dose effect and supports the pathogenic role of ARHGEF9 mutations. This study highlights the role of molecular alteration in phenotype expression and facilitates evaluation of the pathogenicity of ARHGEF9 mutations in clinical practice.

[Efficacy and safety of cyclophosphamide as a sequential immunotherapy drug for anti-N-methyl-D-aspartate receptor encephalitis in children]

OBJECTIVE

To evaluate the efficacy and safety of cyclophosphamide as a second-line drug in the treatment of children with anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis.

METHODS

Six children with anti-NMDAR encephalitis, who showed poor response to steroids and intravenous immunoglobulin, were given cyclophosphamide as a second-line immunotherapy. Follow-up was performed to evaluate the efficacy and safety of cyclophosphamide.

RESULTS

After first-line immunotherapy for 1-4 weeks, the six patients had reduced psychiatric symptoms, seizures, and involuntary movements; three patients had an improved level of consciousness and were able to make simple conversations. However, all the patients still showed slow response, as well as cortical dysfunction symptoms such as aphasia, alexia, agraphia, acalculia, apraxia, and movement disorders. The six patients continued to receive cyclophosphamide as a sequential therapy. They were able to answer simple questions 7 days after treatment. Three school-aged patients were able to make simple calculation, had greatly improved reading and writing ability, and almost recovered self-care ability 2-3 weeks later. The cognitive function of the six patients was almost restored to the level before the onset of disease, and their living ability returned to normal 2-3 months later. During the treatment period, there were no adverse reactions or abnormal results of routine blood test and liver and kidney function tests.

CONCLUSIONS

Children with anti-NMDAR encephalitis should be given appropriate immunotherapy as soon as possible. Cyclophosphamide as a sequential therapy has good efficacy and safety.

HLA-A*24:02 as a common risk factor for antiepileptic drug-induced cutaneous adverse reactions

Objective:

To investigate the involvement of human leukocyte antigen (HLA) loci in aromatic antiepileptic drug–induced cutaneous adverse reactions.

Methods:

A case-control study was performed to detect HLA loci involved in aromatic antiepileptic drug–induced Stevens-Johnson syndrome in a southern Han Chinese population. Between January 1, 2006, and December 31, 2015, 91 cases of Stevens-Johnson syndrome induced by aromatic antiepileptic drugs and 322 matched drug-tolerant controls were enrolled from 8 centers. Important genotypes were replicated in cases with maculopapular eruption and in the meta-analyses of data from other populations. Sequence-based typing determined the HLA-A, HLA-B, HLA-C, and HLA-DRB1 genotypes.

Results:

HLA-B*15:02 was confirmed as strongly associated with carbamazepine-induced Stevens-Johnson syndrome (p = 5.63 × 10⁻¹⁵). In addition, HLA-A*24:02 was associated significantly with Stevens-Johnson syndrome induced by the aromatic antiepileptic drugs as a group (p = 1.02 × 10⁻⁵) and by individual drugs (carbamazepine p = 0.015, lamotrigine p = 0.005, phenytoin p = 0.027). Logistic regression analysis revealed a multiplicative interaction between HLA-B*15:02 and HLA-A*24:02. Positivity for HLA-A*24:02 and/or HLA-B*15:02 showed a sensitivity of 72.5% and a specificity of 69.0%. The presence of HLA-A*24:02 in cases with maculopapular exanthema was also significantly higher than in controls (p = 0.023). Meta-analysis of data from Japan, Korea, Malaysia, Mexico, Norway, and China revealed a similar association.

Conclusions:

HLA-A*24:02 is a common genetic risk factor for cutaneous adverse reactions induced by aromatic antiepileptic drugs in the southern Han Chinese and possibly other ethnic populations. Pretreatment screening is recommended for people in southern China.

A novel role of fragile X mental retardation protein in pre-mRNA alternative splicing through RNA-binding protein 14

Fragile X mental retardation protein (FMRP), an important RNA-binding protein responsible for fragile X syndrome, is involved in posttranscriptional control of gene expression that links with brain development and synaptic functions. Here, we reveal a novel role of FMRP in pre-mRNA alternative splicing, a general event of posttranscriptional regulation. Using co-immunoprecipitation and immunofluorescence assays, we identified that FMRP interacts with an alternative-splicing-associated protein RNA-binding protein 14 (RBM14) in a RNA-dependent fashion, and the two proteins partially colocalize in the nuclei of hippocampal neurons. We show that the relative skipping/inclusion ratio of the micro-exon L in the Protrudin gene and exon 10 in the Tau gene decreased in the hippocampus of Fmr1 knockout (KO) mice. Knockdown of either FMRP or RBM14 alters the relative skipping/inclusion ratio of Protrudin and Tau in cultured Neuro-2a cells, similar to that in the Fmr1 KO mice. Furthermore, overexpression of FMRP leads to an opposite pattern of the splicing, which can be offset by RBM14 knockdown. RNA immunoprecipitation assays indicate that FMRP promotes RBM14's binding to the mRNA targets. In addition, overexpression of the long form of Protrudin or the short form of Tau promotes protrusion growth of the retinoic acid-treated, neuronal-differentiated Neuro-2a cells. Together, these data suggest a novel function of FMRP in the regulation of pre-mRNA alternative splicing through RBM14 that may be associated with normal brain function and FMRP-related neurological disorders.

Epilepsy-associated genes

Development in genetic technology has led to the identification of an increasing number of genes associated with epilepsy. These discoveries will both provide the basis for including genetic tests in clinical practice and improve diagnosis and treatment of epilepsy. By searching through several databases (OMIM, HGMD, and EpilepsyGene) and recent publications on PubMed, we found 977 genes that are associated with epilepsy. We classified these genes into 4 categories according to the manifestation of epilepsy in phenotypes. We found 84 genes that are considered as epilepsy genes: genes that cause epilepsies or syndromes with epilepsy as the core symptom. 73 genes were listed as neurodevelopment-associated genes: genes associated with both brain-development malformations and epilepsy. Several genes (536) were epilepsy-related: genes associated with both physical or other systemic abnormalities and epilepsy or seizures. We found 284 additional genes putatively associated with epilepsy; this requires further verification. These integrated data will provide new insights useful for both including genetic tests in the clinical practice and evaluating the results of genetic tests. We also summarized the epilepsy-associated genes according to their function, with the goal to better characterize the association between genes and epilepsies and to further understand the mechanisms underlying epilepsy.

GAPDH-mediated posttranscriptional regulations of sodium channel Scn1a and Scn3a genes under seizure and ketogenic diet conditions

Abnormal expressions of sodium channel SCN1A and SCN3A genes alter neural excitability that are believed to contribute to the pathogenesis of epilepsy, a long-term risk of recurrent seizures. Ketogenic diet (KD), a high-fat and low-carbohydrate treatment for difficult-to-control (refractory) epilepsy in children, has been suggested to reverse gene expression patterns. Here, we reveal a novel role of GAPDH on the posttranscriptional regulation of mouse Scn1a and Scn3a expressions under seizure and KD conditions. We show that GAPDH binds to a conserved region in the 3' UTRs of human and mouse SCN1A and SCN3A genes, which decreases and increases genes' expressions by affecting mRNA stability through SCN1A 3' UTR and SCN3A 3' UTR, respectively. In seizure mice, the upregulation and phosphorylation of GAPDH enhance its binding to the 3' UTR, which lead to downregulation of Scn1a and upregulation of Scn3a. Furthermore, administration of KD generates β-hydroxybutyric acid which rescues the abnormal expressions of Scn1a and Scn3a by weakening the GAPDH's binding to the element. Taken together, these data suggest that GAPDH-mediated expression regulation of sodium channel genes may be associated with epilepsy and the anticonvulsant action of KD.

Altered intrinsic properties and bursting activities of neurons in layer IV of somatosensory cortex from Fmr-1 knockout mice

Neuroadaptations and alterations in neuronal excitability are critical in brain maturation and many neurological diseases. Fragile X syndrome (FXS) is a pervasive neurodevelopmental disorder characterized by extensive synaptic and circuit dysfunction. It is still unclear about the alterations in intrinsic excitability of individual neurons and their link to hyperexcitable circuitry. In this study, whole cell patch-clamp recordings were employed to characterize the membrane and firing properties of layer IV cells in slices of the somatosensory cortex of Fmr-1 knockout (KO) mice. These cells generally exhibited a regular spiking (RS) pattern, while there were significant increases in the number of cells that adopted intrinsic bursting (IB) compared with age-matched wild type (WT) cells. The cells subgrouped according to their firing patterns and maturation differed significantly in membrane and discharge properties between KO and WT. The changes in the intrinsic properties were consistent with highly facilitated discharges in KO cells induced by current injection. Spontaneous activities of RS neurons driven by local network were also increased in the KO cells, especially in neonate groups. Under an epileptiform condition mimicked by omission of Mg(2+) in extracellular solution, these RS neurons from KO mice were more likely to switch to burst discharges. Analysis on bursts revealed that the KO cells tended to form burst discharges and even severe events manifested as seizure-like ictal discharges. These results suggest that alterations in intrinsic properties in individual neurons are involved in the abnormal excitability of cortical circuitry and possibly account for the pathogenesis of epilepsy in FXS.

Epigenetic Downregulation of Scn3a Expression by Valproate: a Possible Role in Its Anticonvulsant Activity

Upregulation of sodium channel SCN3A expression in epileptic tissues is known to contribute to enhancing neuronal excitability and the development of epilepsy. Therefore, certain strategies to reduce SCN3A expression may be helpful for seizure control. Here, we reveal a novel role of valproate (VPA) in the epigenetic downregulation of Scn3a expression. We found that VPA, instead of carbamazepine (CBZ) and lamotrigine (LTG), could significantly downregulate Scn3a expression in mouse Neuro-2a cells. Luciferase assays and CpG methylation analyses showed that VPA induced the methylation at the -39C site in Scn3a promoter which decreased the promoter activity. We further showed that VPA downregulated the expression of methyl-CpG-binding domain protein 2 (MBD2) at the posttranscriptional level and knockdown of MBD2 increased Scn3a expression. In addition, we found that VPA induced the expression of fat mass and obesity-associated (FTO) protein and FTO knockdown abolished the repressive effects of VPA on MBD2 and Na_v1.3 expressions. Furthermore, VPA, instead of other two anticonvulsant drugs, induced the expressions of Scn3a and Mbd2 and reduced Fto expression in the hippocampus of VPA-treated seizure mice. Taken together, this study suggests an epigenetic pathway for the VPA-induced downregulation of Scn3a expression, which provides a possible role of this pathway in the anticonvulsant action of VPA.

Involvement of FMRP in Primary MicroRNA Processing via Enhancing Drosha Translation

Fragile X mental retardation protein (FMRP), associated with fragile X syndrome, is known as an RNA-binding protein to regulate gene expression at post-transcriptional level in the brain. FMRP is also involved in microRNA (miRNA) biogenesis during the process of precursor miRNA (pre-miRNA) into mature miRNA. However, there is no description of the effect of FMRP on primary miRNA (pri-miRNA) processing. Here, we uncover a novel role of FMRP in pri-miRNA processing via controlling Drosha translation. We show that the expression of DROSHA protein, instead of its messenger RNA (mRNA) transcripts, is downregulated in both the hippocampus of Fmr1-knockout mice and the FMRP-knockdown Neuro-2a cells. Overexpression or knockdown FMRP does not alter Drosha mRNA stability. Immunoprecipitation and polysome analyses demonstrate that FMRP binds to the Drosha mRNA and enhances its translation. Additionally, we show that loss of FMRP in Fmr1-deficient mice results in the accumulation of three in six analyzed pri-miRNAs and the reduction of the corresponding pre-miRNAs and mature miRNAs. Thus, our data suggest that FMRP is involved in pri-miRNA processing via enhancing DROSHA expression that may play an important role in fragile X syndrome.

A novel variant in the 3' UTR of human SCN1A gene from a patient with Dravet syndrome decreases mRNA stability mediated by GAPDH's binding

Mutations in the SCN1A gene-encoding voltage-gated sodium channel α-I subunit (Nav1.1) cause various spectrum of epilepsies including Dravet syndrome (DS), a severe and intractable form. A large number of SCN1A mutations identified from the DS patients lead to the loss of function or truncation of Nav1.1 that result in a haploinsufficiency effects, indicating that the exact expression level of SCN1A should be essential to maintain normal brain function. In this study, we have identified five variants c.*1025T>C, c.*1031A>T, c.*1739C>T, c.*1794C>T and c.*1961C>T in the SCN1A 3' UTR in the patients with DS. The c.*1025T>C, c.*1031A>T and c.*1794C>T are conserved among different species. Of all the five variants, only c.*1794C>T is a novel variant and alters the predicted secondary structure of the 3' UTR. We also show that glyceraldehyde-3-phosphate dehydrogenase (GAPDH) only binds to the 3' UTR sequence containing the mutation allele 1794U but not the wild-type allele 1794C, indicating that the mutation allele forms a new GAPDH-binding site. Functional analyses show that the variant negatively regulates the reporter gene expression by affecting the mRNA stability that is mediated by GAPDH's binding, and this phenomenon could be reversed by shRNA-induced GAPDH knockdown. These findings suggest that GAPDH and the 3'-UTR variant are involved in regulating SCN1A expression at post-transcriptional level, which may provide an important clue for further investigating on the relationship between 3'-UTR variants and SCN1A-related diseases.

Transcription of the human sodium channel SCN1A gene is repressed by a scaffolding protein RACK1

Voltage-gated sodium channel α subunit type I (Nav1.1, encoded by SCN1A gene) plays a critical role in the initiation of action potential in the central nervous system. Downregulated expression of SCN1A is believed to be associated with epilepsy. Here, we found that the SCN1A promoter (P1c), located at the 5' untranslated exon 1c, drove the reporter gene expression in human NT2 cells, and a region between nt +53 and +62 downstream of the P1c promoter repressed the promoter activity. Further analyses showed that a scaffolding protein RACK1 (receptor for activated C kinase 1) was involved in binding to this silencer. Knockdown of RACK1 expression in NT2 cells deprived the repressive role of the silencer on the P1c promoter and increased SCN1A transcription, suggesting the potential involvement of RACK1 in negatively regulating SCN1A transcription via interaction with the silencer. Furthermore, we demonstrated that the binding of the protein complex including RACK1 to the SCN1A promoter motif was decreased in neuron-like differentiation of the NT2 cells induced by retinoic acid and resulted in the upregulation of SCN1A transcription. Taken together, this study reports a novel role of RACK1 in regulating SCN1A expression that participates in retinoic acid-induced neuronal differentiation of NT2 cells.

A conserved region in the 3' untranslated region of the human LIMK1 gene is critical for proper expression of LIMK1 at the post-transcriptional level

LIM kinase 1 (LIMK1), a cytosolic serine/threonine kinase, regulates actin filament dynamics and reorganization and is involved in neuronal development and brain function. Abnormal expression of LIMK1 is associated with several neurological disorders. In this study, we performed a conservation analysis using Vector NTI (8.0) software. The dualluciferase reporter assay and real-time quantitative RT-PCR were used to assess the protein and mRNA levels of the reporter gene, respectively. We found that a region ranging from nt +884 to +966 in the human LIMK1 3' untranslated region (UTR) was highly conserved in the mouse Limk1 3' UTR and formed a structure containing several loops and stems. Luciferase assay showed that the relative luciferase activity of the mutated construct with the conserved region deleted, pGL4-hLIMK1-3U-M, in SH-SY5Y and HEK-293 cells was only ~60% of that of the wild-type construct pGL4-hLIMK1-3U, indicating that the conserved region is critical for the reporter gene expression. Real-time quantitative RT-PCR analysis demonstrated that the relative Luc2 mRNA levels in SH-SY5Y and HEK293 cells transfected with pGL4-hLIMK1-3U-M decreased to ~50% of that in cells transfected with pGL4-hLIMK1-3U, suggesting an important role of the conserved region in maintaining Luc2 mRNA stability. Our study suggests that the conserved region in the LIMK1 3' UTR is involved in regulating LIMK1 expression at the post-transcriptional level, which may help reveal the mechanism underlying the regulation of LIMK1 expression in the central nervous system and explore the relationship between the 3'-UTR mutant and neurological disorders.

Three-dimensional characterization of resorption cavity size and location in human vertebral trabecular bone

The number and size of resorption cavities in cancellous bone are believed to influence rates of bone loss, local tissue stress and strain and potentially whole bone strength. Traditional two-dimensional approaches to measuring resorption cavities in cancellous bone report the percent of the bone surface covered by cavities or osteoclasts, but cannot measure cavity number or size. Here we use three-dimensional imaging (voxel size 0.7×0.7×5.0 μm) to characterize resorption cavity location, number and size in human vertebral cancellous bone from nine elderly donors (7 male, 2 female, ages 47-80 years). Cavities were 30.10 ± 8.56 μm in maximum depth, 80.60 ± 22.23∗10(3) μm(2) in surface area and 614.16 ± 311.93∗10(3) μm(3) in volume (mean ± SD). The average number of cavities per unit tissue volume (N.Cv/TV) was 1.25 ± 0.77 mm(-3). The ratio of maximum cavity depth to local trabecular thickness was 30.46 ± 7.03% and maximum cavity depth was greater on thicker trabeculae (p<0.05, r(2)=0.14). Half of the resorption cavities were located entirely on nodes (the intersection of two or more trabeculae) within the trabecular structure. Cavities that were not entirely on nodes were predominately on plate-like trabeculae oriented in the cranial-caudal (longitudinal) direction. Cavities on plate-like trabeculae were larger in maximum cavity depth, cavity surface area and cavity volume than cavities on rod-like trabeculae (p<0.05). We conclude from these findings that cavity size and location are related to local trabecular microarchitecture.

Neuroprotection of lamotrigine on hypoxic-ischemic brain damage in neonatal rats: Relations to administration time and doses

Lamotrigine (LTG), an antiepileptic drug, has been shown to be able to improve cerebral ischemic damage by limiting the presynaptic release of glutamate. The present study investigated further the neuroprotective effect of LTG on hypoxic-ischemic brain damage (HIBD) in neonatal rats and its relations to administration time and doses. The HIBD model was produced in 7-days old SD rats by left common carotid artery ligation followed by 2 h hypoxic exposure (8% oxygen). LTG was administered intraperitoneally with the doses of 5, 10, 20, and 40 mg/kg 3 h after operation and the dose of 20 mg/kg 1 h before and 3 h, 6 h after operation. Blood and brain were sampled 24 h after operation. Nissl staining, terminal deoxynucleotidyl transferase biotin-dUTP nick end labeling (TUNEL), and neuron-specific enolase (NSE) immunohistochemical staining were used for morphological studies. Water content in left cortex and NSE concentration in serum were determined. LTG significantly reduced water content in the cerebral cortex, as well as the number of TUNEL staining neurons in the dentate gyrus and cortex in hypoxic-ischemia (HI) model. Furthermore, LTG significantly decreased the NSE level in serum and increased the number of NSE staining neurons in the cortex. These effects, except that on water content, were dose-dependent and were more remarkable in the pre-treated group than in the post-treated groups. These results demonstrate that LTG may have a neuroprotective effect on acute HIBD in neonates. The effect is more prominent when administrated with higher doses and before HI.