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

Using computational approaches to enhance the interpretation of missense variants in the PAX6 gene

The PAX6 gene encodes a highly-conserved transcription factor involved in eye development. Heterozygous loss-of-function variants in PAX6 can cause a range of ophthalmic disorders including aniridia. A key molecular diagnostic challenge is that many PAX6 missense changes are presently classified as variants of uncertain significance. While computational tools can be used to assess the effect of genetic alterations, the accuracy of their predictions varies. Here, we evaluated and optimised the performance of computational prediction tools in relation to PAX6 missense variants. Through inspection of publicly available resources (including HGMD, ClinVar, LOVD and gnomAD), we identified 241 PAX6 missense variants that were used for model training and evaluation. The performance of ten commonly used computational tools was assessed and a threshold optimization approach was utilized to determine optimal cut-off values. Validation studies were subsequently undertaken using PAX6 variants from a local database. AlphaMissense, SIFT4G and REVEL emerged as the best-performing predictors; the optimized thresholds of these tools were 0.967, 0.025, and 0.772, respectively. Combining the prediction from these top-three tools resulted in lower performance compared to using AlphaMissense alone. Tailoring the use of computational tools by employing optimized thresholds specific to PAX6 can enhance algorithmic performance. Our findings have implications for PAX6 variant interpretation in clinical settings.

Structural mapping of patient-associated KCNMA1 gene variants

KCNMA1-linked channelopathy is a neurological disorder characterized by seizures, motor abnormalities, and neurodevelopmental disabilities. The disease mechanisms are predicted to result from alterations in KCNMA1-encoded BK K+ channel activity; however, only a subset of the patient-associated variants have been functionally studied. The localization of these variants within the tertiary structure or evaluation by pathogenicity algorithms has not been systematically assessed. In this study, 82 nonsynonymous patient-associated KCNMA1 variants were mapped within the BK channel protein. Fifty-three variants localized within cryoelectron microscopy-resolved structures, including 21 classified as either gain of function (GOF) or loss of function (LOF) in BK channel activity. Clusters of LOF variants were identified in the pore, the AC region (RCK1), and near the Ca2+ bowl (RCK2), overlapping with sites of pharmacological or endogenous modulation. However, no clustering was found for GOF variants. To further understand variants of uncertain significance (VUSs), assessments by multiple standard pathogenicity algorithms were compared, and new thresholds for sensitivity and specificity were established from confirmed GOF and LOF variants. An ensemble algorithm was constructed (KCNMA1 meta score (KMS)), consisting of a weighted summation of this trained dataset combined with a structural component derived from the Ca2+-bound and unbound BK channels. KMS assessment differed from the highest-performing individual algorithm (REVEL) at 10 VUS residues, and a subset were studied further by electrophysiology in HEK293 cells. M578T, E656A, and D965V (KMS+;REVEL-) were confirmed to alter BK channel properties in voltage-clamp recordings, and D800Y (KMS-;REVEL+) was assessed as benign under the test conditions. However, KMS failed to accurately assess K457E. These combined results reveal the distribution of potentially disease-causing KCNMA1 variants within BK channel functional domains and pathogenicity evaluation for VUSs, suggesting strategies for improving channel-level predictions in future studies by building on ensemble algorithms such as KMS.

ADGRL1 variants: From developmental and epileptic encephalopathy to genetic epilepsy with febrile seizures plus

AIM

To expand the phenotypic spectrum of ADGRL1 and explore the correlation between epilepsy and the ADGRL1 genotype.

METHOD

We performed whole-exome sequencing in a cohort of 115 families (including 195 males and 150 females) with familial febrile seizure or epilepsy with unexplained aetiology. The damaging effects of variants was predicted using protein modelling and multiple in silico tools. All reported patients with ADGRL1 pathogenic variants were analysed.

RESULTS

One new ADGRL1 variant (p.Pro753Leu) was identified in one family with genetic epilepsy with febrile seizures. Further analysis of 12 ADGRL1 variants in 16 patients revealed that six patients had epilepsy. Epilepsy types ranged from early-onset epileptic encephalopathy to genetic epilepsy with febrile seizures plus (GEFS+). All four variants associated with epilepsy were located in the non-helix or sheet region of ADGRL1. Three of the four epilepsy-associated variants were missense variants. Thus, all three variants located in the G-protein-coupled receptor autoproteolytic-inducing domain exhibited epilepsy.

INTERPRETATION

We found one new missense variant of ADGRL1 in one family with GEFS+. ADGRL1 may be a potential candidate or susceptibility gene for epilepsy. ADGRL1-associated epilepsy ranged from benign GEFS+ to severe epileptic encephalopathy; the genotypes and variant locations may help explain the phenotypic heterogeneity of patients with the ADGRL1 variant.

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.

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.

Investigation of FRMPD4 variants associated with X-linked epilepsy

BACKGROUND

The etiology of unexplained epilepsy in most patients remains unclear. Variants of FRMPD4 are suggested to be associated with neurodevelopmental disorders. Therefore, we screened for disease-causing FRMPD4 variants in patients with epilepsy.

METHODS

Trios-based whole-exome sequencing was conducted on a cohort of 85 patients with unexplained epilepsy, their parents, and extended family members. Additional cases with FRMPD4 variants were identified from the China Epilepsy Gene Matching Platform V.1.0. The frequency of variants was analyzed, and their subregional effects were predicted using in silico tools. The genotype-phenotype correlation of the newly defined causative genes and protein stability were analyzed using I-Mutant V.3.0 and Grantham scores.

RESULTS

Two novel missense variants of FRMPD4 were identified in two families. Using the gene matching platform, we identified three additional novel missense variants. These variants presented at low or no allele frequencies in the gnomAD database. All the variants were located outside the three FRMPD4 main domains (WW, PDZ, and FERM). In silico analyses revealed that the variants were damaging and were predicted to be the least stable. All patients eventually became seizure-free. Eight of the 21 patients with FRMPD4 variants had epilepsy, of which five (63%) had missense variants located outside the domains, two had deletions involving exon 2, and one had a frameshift variant located outside the domains. Patients with epilepsy caused by missense variants were often free of intellectual disabilities (4/5), whereas patients with epilepsy caused by truncated variants had intellectual disabilities and structural brain abnormalities (3/3).

CONCLUSIONS

The FRMPD4 gene is potentially associated with epilepsy. The genotype-phenotype correlation of FRMPD4 variants indicated that differences in variant types and locations of FRMPD4 may explain their phenotypic variation.

SCAF4 variants are associated with epilepsy with neurodevelopmental disorders

AIMS

The genetic causes of epilepsy with unknown etiology in most patients remain unknown. The aim of this study was to elucidate the phenotype of SCAF4-related epilepsy.

METHODS

Trio-based whole-exome sequencing was performed in patients with epilepsy. Silico programs and protein modeling were employed to predict the damaging of variants. Previously reported SCAF4 variants were systematically reviewed to analyze the genotype-phenotype correlations.

RESULTS

Three heterozygous variants in the SCAF4 were detected in three cases, including one missense variant and two frameshift variants. All variants were de novo. None of these variants is present in gnomAD controls. The missense variant was predicted to be damaging in silico tools. Protein modeling showed that two frameshift variants resulted in loss of domains, and the missense variant may disrupt a nearby phosphorylation site and alter the hydrogen bonds around 54C and the stability of the SCAF4 protein. Intellectual development was mildly delayed for all patients except for one with whom contact was lost. All probands experienced epilepsy as infrequent seizures, responded well to antiseizure drugs, and had a median [IQR] seizure onset age of 4 [1.75, 7.5] years. The variants in the domain-encoding exons and upstream exons exhibited a strong association with epilepsy.

CONCLUSIONS

SCAF4 is a potential causative gene of epilepsy with neurodevelopmental disorders.

DLG3 variants caused X-linked epilepsy with/without neurodevelopmental disorders and the genotype-phenotype correlation

Background

The DLG3 gene encodes disks large membrane-associated guanylate kinase scaffold protein 3, which plays essential roles in the clustering of N-methyl-D-aspartate receptors (NMDARs) at excitatory synapses. Previously, DLG3 has been identified as the causative gene of X-linked intellectual developmental disorder-90 (XLID-90; OMIM# 300850). This study aims to explore the phenotypic spectrum of DLG3 and the genotype-phenotype correlation.

Methods

Trios-based whole-exome sequencing was performed in patients with epilepsy of unknown causes. To analyze the genotype-phenotype correlations, previously reported DLG3 variants were systematically reviewed.

Results

DLG3 variants were identified in seven unrelated cases with epilepsy. These variants had no hemizygous frequencies in controls. All variants were predicted to be damaging by silico tools and alter the hydrogen bonds with surrounding residues and/or protein stability. Four cases mainly presented with generalized seizures, including generalized tonic-clonic and myoclonic seizures, and the other three cases exhibited secondary generalized tonic-clonic seizures and focal seizures. Multifocal discharges were recorded in all cases during electroencephalography monitoring, including the four cases with generalized discharges initially but multifocal discharges after drug treating. Protein-protein interaction network analysis revealed that DLG3 interacts with 52 genes with high confidence, in which the majority of disease-causing genes were associated with a wide spectrum of neurodevelopmental disorder (NDD) and epilepsy. Three patients with variants locating outside functional domains all achieved seizure-free, while the four patients with variants locating in functional domains presented poor control of seizures. Analysis of previously reported cases revealed that patients with non-null variants presented higher percentages of epilepsy than those with null variants, suggesting a genotype-phenotype correlation.

Significance

This study suggested that DLG3 variants were associated with epilepsy with/without NDD, expanding the phenotypic spectrum of DLG3. The observed genotype-phenotype correlation potentially contributes to the understanding of the underlying mechanisms driving phenotypic variation.

Structural mapping of patient-associated KCNMA1 gene variants

KCNMA1-linked channelopathy is a neurological disorder characterized by seizures, motor abnormalities, and neurodevelopmental disabilities. The disease mechanisms are predicted to result from alterations in KCNMA1-encoded BK K+ channel activity; however, only a subset of the patient-associated variants have been functionally studied. The localization of these variants within the tertiary structure or evaluation by pathogenicity algorithms has not been systematically assessed. In this study, 82 nonsynonymous patient-associated KCNMA1 variants were mapped within the BK channel protein. Fifty-three variants localized within cryo-EM resolved structures, including 21 classified as either gain-of-function (GOF) or loss-of-function (LOF) in BK channel activity. Clusters of LOF variants were identified in the pore, the AC region (RCK1), and near the Ca 2+ bowl (RCK2), overlapping with sites of pharmacological or endogenous modulation. However, no clustering was found for GOF variants. To further understand variants of uncertain significance (VUS), assessments by multiple standard pathogenicity algorithms were compared, and new thresholds for sensitivity and specificity were established from confirmed GOF and LOF variants. An ensemble algorithm was constructed (KCNMA1 Meta Score), consisting of a weighted summation of this trained dataset combined with a structural component derived from the Ca 2+ bound and unbound BK channels. KMS assessment differed from the highest performing individual algorithm (REVEL) at 10 VUS residues, and a subset were studied further by electrophysiology in HEK293 cells. M578T, E656A, and D965V (KMS+;REVEL-) were confirmed to alter BK channel properties in voltage-clamp recordings, and D800Y (KMS-;REVEL+) was assessed as benign under the test conditions. However, KMS failed to accurately assess K457E. These combined results reveal the distribution of potentially disease-causing KCNMA1 variants within BK channel functional domains and pathogenicity evaluation for VUS, suggesting strategies for improving channel-level predictions in future studies by building on ensemble algorithms such as KMS.

Functional characterization of novel NPRL3 mutations identified in three families with focal epilepsy

Focal epilepsy accounts for 60% of all forms of epilepsy, but the pathogenic mechanism is not well understood. In this study, three novel mutations in NPRL3 (nitrogen permease regulator-like 3), c.937_945del, c.1514dupC and 6,706-bp genomic DNA (gDNA) deletion, were identified in three families with focal epilepsy by linkage analysis, whole exome sequencing (WES) and Sanger sequencing. NPRL3 protein is a component of the GATOR1 complex, a major inhibitor of mTOR signaling. These mutations led to truncation of the NPRL3 protein and hampered the binding between NPRL3 and DEPDC5, which is another component of the GATOR1 complex. Consequently, the mutant proteins enhanced mTOR signaling in cultured cells, possibly due to impaired inhibition of mTORC1 by GATOR1. Knockdown of nprl3 in Drosophila resulted in epilepsy-like behavior and abnormal synaptic development. Taken together, these findings expand the genotypic spectrum of NPRL3-associated focal epilepsy and provide further insight into how NPRL3 mutations lead to epilepsy.

CD40LG-associated X-linked Hyper-IgM Syndrome (XHIGM) with pulmonary alveolar proteinosis: a case report

BACKGROUND

D40LG-associated X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis has rarely been reported, and its genotype-phenotypic correlation remains elusive.

CASE PRESENTATION

We describe a five-month-old boy with CD40LG mutation (c.516T > A, p.Tyr172Ter) X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis as the first manifestation. The patient completely recovered after immunotherapy and allogeneic hematopoietic stem cell transplantation. In addition, four previously reported patients with CD40LG mutation with pulmonary alveolar proteinosis were also analyzed. All of these patients presented with early onset of pulmonary infections and a good response to immunotherapy. The structural model of CD40LG indicated that all mutations caused the X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis to be located within the tumor necrosis factor homology domain.

CONCLUSIONS

A case was presented, and the characteristics of four cases of CD40LG-associated X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis were summarized. The variant locations may explain the phenotypic heterogeneity of patients with the CD40LG mutation.

Expanding the phenotypic spectrum of KCNK4: From syndromic neurodevelopmental disorder to rolandic epilepsy

The KCNK4 gene, predominantly distributed in neurons, plays an essential role in controlling the resting membrane potential and regulating cellular excitability. Previously, only two variants were identified to be associated with human disease, facial dysmorphism, hypertrichosis, epilepsy, intellectual/developmental delay, and gingival overgrowth (FHEIG) syndrome. In this study, we performed trio-based whole exon sequencing (WES) in a cohort of patients with epilepsy. Two de novo likely pathogenic variants were identified in two unrelated cases with heterogeneous phenotypes, including one with Rolandic epilepsy and one with the FHEIG syndrome. The two variants were predicted to be damaged by the majority of in silico algorithms. These variants showed no allele frequencies in controls and presented statistically higher frequencies in the case cohort than that in controls. The FHEIG syndrome-related variants were all located in the region with vital functions in stabilizing the conductive conformation, while the Rolandic epilepsy-related variant was distributed in the area with less impact on the conductive conformation. This study expanded the genetic and phenotypic spectrum of KCNK4. Phenotypic variations of KCNK4 are potentially associated with the molecular sub-regional effects. Carbamazepine/oxcarbazepine and valproate may be effective antiepileptic drugs for patients with KCNK4 variants.

Phenotypic expansion of KCNH1-associated disorders to include isolated epilepsy and its associations with genotypes and molecular sub-regional locations

PURPOSE

Genotype-phenotypic correlation of KCNH1 variant remains elusive. This study aimed to expand the phenotypic spectrum of KCNH1 and explore the correlations between epilepsy and molecular sub-regional locations.

METHODS

We performed whole-exome sequencing in a cohort of 98 patients with familiar febrile seizure (FS) or epilepsy with unexplained etiologies. The damaging effects of variants were predicted by protein modeling and multiple in silico tools. All reported patients with KCNH1 pathogenic variants with detailed neurological phenotypes were analyzed to evaluate the genotype-phenotype correlation.

RESULTS

Two novel KCNH1 variants were identified in three cases, including two patients with FS with inherited variant (p.Ile113Thr) and one boy with epilepsy with de novo variant (p.Arg357Trp). Variant Ile113Thr was located within the eag domain, and variant p.Arg357Trp was located in transmembrane domain 4 of KCNH1, respectively. Two patients experienced refractory status epilepticus (SE), of which one patient died of acute encephalopathy induced by SE. Further analysis of 30 variants in 51 patients demonstrated that de novo variants were associated with epileptic encephalopathy, while mosaic/somatic or germline variants cause isolated epilepsy/FS. All hotspot variants associated with epileptic encephalopathy clustered in transmembrane domain (S4 and S6), while those with isolated epilepsy/seizures or TBS/ZLS without epilepsy were scattered in the KCNH1.

CONCLUSIONS

We found two novel missense variants of KCNH1 in three individuals with isolated FS/epilepsy. Variants in the KCNH1 cause a spectrum of epileptic disorders ranging from a benign form of genetic isolated epilepsy/FS to intractable form of epileptic encephalopathy. The genotypes and variant locations help explaining the phenotypic variation of patients with KCNH1 variant.

Variant predictions in congenital adrenal hyperplasia caused by mutations in CYP21A2

CYP21A2 deficiency represents 95% of congenital adrenal hyperplasia (CAH) cases, a group of genetic disorders that affect steroid biosynthesis. The genetic and functional analysis provide critical tools to elucidate complex CAH cases. One of the most accessible tools to infer the pathogenicity of new variants is in silico prediction. Here, we analyzed the performance of in silico prediction tools to categorize missense single nucleotide variants (SNVs) of CYP21A2. SNVs of CYP21A2 characterized in vitro by functional assays were selected to assess the performance of online single and meta predictors. SNVs were tested separately or in combination with the related phenotype (severe or mild CAH form). In total, 103 SNVs of CYP21A2 (90 pathogenic and 13 neutral) were used to test the performance of 13 single-predictors and four meta-predictors. All SNVs associated with the severe phenotypes were well categorized by all tools, with an accuracy of between 0.69 (PredictSNP2) and 0.97 (CADD), and Matthews’ correlation coefficient (MCC) between 0.49 (PoredicSNP2) and 0.90 (CADD). However, SNVs related to the mild phenotype had more variation, with the accuracy between 0.47 (S3Ds&GO and MAPP) and 0.88 (CADD), and MCC between 0.18 (MAPP) and 0.71 (CADD). From our analysis, we identified four predictors of CYP21A2 variant pathogenicity with good performance, CADD, ConSurf, DANN, and PolyPhen2. These results can be used for future analysis to infer the impact of uncharacterized SNVs in CYP21A2.

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.

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.

AFF2 Is Associated With X-Linked Partial (Focal) Epilepsy With Antecedent Febrile Seizures

OBJECTIVE

AFF2 mutations were associated with X-linked intellectual developmental disorder-109 and in males with autism spectrum disorder (ASD). The relationship between AFF2 and epilepsy has not been defined.

METHOD

Trios-based whole-exome sequencing was performed in a cohort of 372 unrelated cases (families) with partial (focal) epilepsy without acquired causes.

RESULTS

Five hemizygous missense AFF2 mutations were identified in five males with partial epilepsy and antecedent febrile seizures without intellectual disability or other developmental abnormalities. The mutations did not present in the controls of general populations with an aggregate frequency significantly higher than that in the control populations. Previously, intellectual disability-associated AFF2 mutations were genomic rearrangements and CCG repeat expansion mutations mostly, whereas the mutations associated with partial epilepsy were all missense. Missense AFF2 mutations associated with epilepsy fell into the regions from N-terminal to the nuclear localization signal 1 (NLS1), while ASD-associated missense mutations fell in the regions from NLS1 to C-terminal.

CONCLUSION

AFF2 is potentially a candidate causative gene of X-link partial epilepsy with antecedent febrile seizures. The genotype-phenotype correlation and molecular sub-regional effect of AFF2 help in explaining the mechanisms underlying phenotypic variations.

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.

Biallelic ADGRV1 variants are associated with Rolandic epilepsy

OBJECTIVE

Rolandic epilepsy (RE) is among the most common focal epilepsies in childhood. For the majority of patients with RE and atypical RE (ARE), the etiology remains elusive. We thus screened patients with RE/ARE in order to detect disease-causing variants..

METHODS

A trios-based whole-exome sequencing approach was performed in a cohort of 28 patients with RE/ARE. Clinical data and EEGs were reviewed. Variants were validated by Sanger sequencing.

RESULTS

Two compound heterozygous missense variants p.Val272Ile/p.Asn3028Ser and p.Ala3657Val/p.Met4419Val of ADGRV1 were identified in two unrelated familial cases of RE/ARE. All the variants were in the calcium exchanger β domain and were suggested to be damaging by at least one web-based prediction tool. These variants are not present or are present at a very low minor allele frequency in the gnomAD database. Previously, biallelic ADGRV1 variants (p.Gly2756Arg and p.Glu4410Lys) have been observed in RE, consistent with the observation in this study and supporting the association between ADGRV1 variants and RE. Additionally, a de novo mutation, p.Asp668Asn, in GRIN2B was identified in a sporadic case of ARE, and a missense variant, p.Asn1551Ser, in RyR2 was identified in a family with RE with incomplete penetrance. These genes are all calcium homeostasis associated genes, suggesting the potential effect of calcium homeostasis in RE/ARE.

CONCLUSIONS

The results from the present study suggest that the genes ADGRV1, GRIN2B, and RyR2 are associated with RE/ARE. These data link defects in neuronal intracellular calcium homeostasis to RE/ARE pathogenesis implicating that these defects plays an important role in the development of these conditions.

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.

Bridging heterogeneous mutation data to enhance disease gene discovery

Bridging heterogeneous mutation data fills in the gap between various data categories and propels discovery of disease-related genes. It is known that genome-wide association study (GWAS) infers significant mutation associations that link genotype and phenotype. However, due to the differences of size and quality between GWAS studies, not all de facto vital variations are able to pass the multiple testing. In the meantime, mutation events widely reported in literature unveil typical functional biological process, including mutation types like gain of function and loss of function. To bring together the heterogeneous mutation data, we propose a 'Gene-Disease Association prediction by Mutation Data Bridging (GDAMDB)' pipeline with a statistic generative model. The model learns the distribution parameters of mutation associations and mutation types and recovers false-negative GWAS mutations that fail to pass significant test but represent supportive evidences of functional biological process in literature. Eventually, we applied GDAMDB in Alzheimer's disease (AD) and predicted 79 AD-associated genes. Besides, 12 of them from the original GWAS, 60 of them are supported to be AD-related by other GWAS or literature report, and rest of them are newly predicted genes. Our model is capable of enhancing the GWAS-based gene association discovery by well combining text mining results. The positive result indicates that bridging the heterogeneous mutation data is contributory for the novel disease-related gene discovery.

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.

YWHAG Mutations Cause Childhood Myoclonic Epilepsy and Febrile Seizures: Molecular Sub-regional Effect and Mechanism

YWHAG, which encodes an adapter protein 14-3-3γ, is highly expressed in the brain and regulates a diverse range of cell signaling pathways. Previously, eight YWHAG mutations have been identified in patients with epileptic encephalopathy (EE). In this study, using trios-based whole exome sequencing, we identified two novel YWHAG mutations in two unrelated families with childhood myoclonic epilepsy and/or febrile seizures (FS). The identified mutations included a heterozygous truncating mutation (c.124C>T/p.Arg42Ter) and a de novo missense mutation (c.373A>G/p.Lys125Glu). The two probands experienced daily myoclonic seizures that were recorded with ictal generalized polyspike-slow waves, but became seizure-free with simple valproate treatment. The other affected individuals presented FS. The truncating mutation was identified in the family with six individuals of mild phenotype, suggesting that YWHAG mutations of haploinsufficiency are relatively less pathogenic. Analysis on all missense mutations showed that nine mutations were located within 14-3-3γ binding groove and another mutation was located at residues critical for dimerization, indicating a molecular sub-regional effect. Mutation Arg132Cys, which was identified recurrently in five patients with EE, would have the strongest influence on binding affinity. 14-3-3γ dimers supports target proteins activity. Thus, a heterozygous missense mutation would lead to majority dimers being mutants; whereas a heterozygous truncating mutation would lead to only decreasing the number of wild-type dimer, being one of the explanations for phenotypical variation. This study suggests that YWHAG is potentially a candidate pathogenic gene of childhood myoclonic epilepsy and FS. The spectrum of epilepsy caused by YWHAG mutations potentially range from mild myoclonic epilepsy and FS to severe EE.

Targeted High-Throughput Sequencing Analysis Results of Osteogenesis Imperfecta Patients from Different Regions of Turkey

Objective: Osteogenesis imperfecta (OI) includes a group of disorders characterized by susceptibility to bone fractures with different severities. The increasing number of genes that may underlie the disorder, along with the broad phenotypic spectrum that overlaps with other skeletal diseases, provided a compelling case for the use of high-throughput sequencing (HTS) technology as an aid to OI diagnoses. The aim of this analysis was to present the data from our 5-year targeted HTS results, that includes the reporting of 9 novel and 24 known mutations, found in OI patients, from 5 different regions of Turkey. Materials and Methods: We performed a retrospective cross-sectional study, reporting the HTS results of 43 patients (23 female and 20 male; mean age: 9.5 years), directed to our center with a suspicion of OI between February 2015 and May 2020. Genetic analyses were also performed for 24 asymptomatic parents to aid the segregation analyses. We utilized an HTS panel targeting the coding regions of 57 genes associated with a reduction, increase, or abnormal development of bone mineralization. In addition, we sequenced the entire coding region of the IFITM5 gene through HTS. Results: Thirty-nine patients had at least one pathogenic/likely pathogenic variation (90.69%) in the COL1A1 (56.41%), COL1A2 (20.51%), FKBP10 (7.7%), P3H1 (5.13%), IFITM5 (5.13%), CTRAP (2.56%), or TMEM38B (2.56%) genes. Nine of the determined pathogenic/likely pathogenic variations were novel. The recurrent pathogenic mutations were c.1081C>T (p.Arg361Ter) (3/43), c.1405C>T (p.Arg469Ter) (2/43), and c.3749del (p.Gly1250AlafsTer81) in COL1A1 gene, along with c.-14C>T variation in the 5'UTR of the IFITM5 gene (2/43) and the c.890_897dup variation in the FKBP10 gene (2/43). Three out of 43 patients were carrying at least one additional variant of unknown significance, highlighting the importance of a multigene panel approach and segregation analyses. Conclusion: We suggest that a targeted HTS panel is a feasible tool for genetic diagnosis of OI in patients.

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.

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.