Advances in big data and omics: Paving the way for discovery in childhood epilepsies

The insights gained from big data and omics approaches have transformed the field of childhood genetic epilepsy. With an increasing number of individuals receiving genetic testing for seizures, we are provided with an opportunity to identify clinically relevant subgroups and extract meaningful observations from this large-scale clinical data. However, the volume of data from electronic medical records and omics (e.g., genomics, transcriptomics) is so vast that standardized methods, such as the Human Phenotype Ontology, are necessary for reliable and comprehensive characterization. Here, we explore the integration of clinical and omics data, highlighting how these approaches pave the way for discovery in childhood epilepsies.

Expanding the clinical phenotype and variant spectrum associated with RFX7

RFX7 encodes a transcription factor that is ubiquitously expressed and important for neural development. Haploinsufficiency of RFX7 is associated with intellectual disability, developmental delay, and diverse malformations of brain structures. Currently, there are only 16 clinically described individuals who have variants in RFX7. A recognizable pattern of malformation associated with mutation in RFX7 has not yet been uncovered. Here we describe the phenotypic presentation of two additional individuals who have novel de novo variants in RFX7. One of the individuals we describe is from an under-represented Afro-Caribbean population.

Clinical signatures of genetic epilepsies precede diagnosis in electronic medical records of 32,000 individuals

PURPOSE

An early genetic diagnosis can guide the time-sensitive treatment of individuals with genetic epilepsies. However, most genetic diagnoses occur long after disease onset. We aimed to identify early clinical features suggestive of genetic diagnoses in individuals with epilepsy through large-scale analysis of full-text electronic medical records.

METHODS

We extracted 89 million time-stamped standardized clinical annotations using Natural Language Processing from 4,572,783 clinical notes from 32,112 individuals with childhood epilepsy, including 1925 individuals with known or presumed genetic epilepsies. We applied these features to train random forest models to predict SCN1A-related disorders and any genetic diagnosis.

RESULTS

We identified 47,774 age-dependent associations of clinical features with genetic etiologies a median of 3.6 years before molecular diagnosis. Across all 710 genetic etiologies identified in our cohort, neurodevelopmental differences between 6 to 9 months increased the likelihood of a later molecular diagnosis 5-fold (P < .0001, 95% CI = 3.55-7.42). A later diagnosis of SCN1A-related disorders (area under the curve [AUC] = 0.91) or an overall positive genetic diagnosis (AUC = 0.82) could be reliably predicted using random forest models.

CONCLUSION

Clinical features predictive of genetic epilepsies precede molecular diagnoses by up to several years in conditions with known precision treatments. An earlier diagnosis facilitated by automated electronic medical records analysis has the potential for earlier targeted therapeutic strategies in the genetic epilepsies.

A Longitudinal Exploration of CACNA1A -related Hemiplegic Migraine in Children

Introduction

Since the initial description of CACNA1A- related hemiplegic migraine (HM), the phenotypic spectrum has expanded from mild episodes in neurotypical individuals to potentially life-threatening events frequently seen in individuals with developmental and epileptic encephalopathies. However, the overall longitudinal course throughout childhood remains unknown.

Methods

We analyzed HM and seizure history in individuals with CACNA1A -related HM, delineating frequency and severity of events in monthly increments through a standardized approach. Combining these data with medication prescription information, we assessed the response of HM to different agents.

Results

Our cohort involved 15 individuals between 3 and 29 years (163 patient years) and included 11 unique and two recurrent variants (p.R1349Q and p.V1393M; both n= 2). The age of first confirmed HM ranged from 14 months to 13 years (average 3 years). 25% of all HM events were severe (lasting >3 days) and 73% of individuals had at least 1 severe occurrence. Spacing of HM events ranged from 1 month to 14 years and changes in HM severity over time of showed increases or decreases of >2 severity levels in 12/122 events. Eight individuals had epilepsy, but severity of epilepsy did not correlate with frequency and severity of HM events. While levetiracetam ( n= 6) and acetazolamide ( n= 5) were the most frequently used medications, they did not show efficacy in HM prevention or HM severity reduction. However, verapamil ( n= 3) showed efficacy in preventing HM episodes (OR 2.68, CI 1.39-5.67).

Significance

The longitudinal course of CACNA1A -related HM lacks recognizable patterns for timing and severity of HM events or correlation with seizure patterns. Our data underscores the unpredictability of CACNA1A -related HM, highlighting the need for close surveillance for reoccurring HM events even in individuals with symptom-free periods.

Key points

24% of hemiplegic migraines (HM) in CACNA1A- related disorders are severe, involving cerebral edema and greater than 4 days to recover Timing and severity of HM are unpredictable, with large changes in severity between events, and age of onset ranging from 1-13 yearsEpilepsy occurred in 53% of individuals, with neither the timing nor severity of seizures correlated with HM.

A Comparison of Ketamine and Midazolam as First-Line Anesthetic Infusions for Pediatric Status Epilepticus

BACKGROUND

Pediatric refractory status epilepticus (RSE) often requires management with anesthetic infusions, but few data compare first-line anesthetics. This study aimed to compare the efficacy and adverse effects of midazolam and ketamine infusions as first-line anesthetics for pediatric RSE.

METHODS

Retrospective single-center study of consecutive study participants treated with ketamine or midazolam as the first-line anesthetic infusions for RSE at a quaternary care children's hospital from December 1, 2017, until September 15, 2021.

RESULTS

We identified 117 study participants (28 neonates), including 79 (68%) who received midazolam and 38 (32%) who received ketamine as the first-line anesthetic infusions. Seizures terminated more often in study participants administered ketamine (61%, 23/38) than midazolam (28%, 22/79; odds ratio [OR] 3.97, 95% confidence interval [CI] 1.76-8.98; P < 0.01). Adverse effects occurred more often in study participants administered midazolam (24%, 20/79) than ketamine (3%, 1/38; OR 12.54, 95% CI 1.61-97.43; P = 0.016). Study participants administered ketamine were younger, ketamine was used more often for children with acute symptomatic seizures, and midazolam was used more often for children with epilepsy. Multivariable logistic regression of seizure termination by first-line anesthetic infusion (ketamine or midazolam) including age at SE onset, SE etiology category, and individual seizure duration at anesthetic infusion initiation indicated seizures were more likely to terminate following ketamine than midazolam (OR 4.00, 95% CI 1.69-9.49; P = 0.002) and adverse effects were more likely following midazolam than ketamine (OR 13.41, 95% CI 1.61-111.04; P = 0.016). Survival to discharge was higher among study participants who received midazolam (82%, 65/79) than ketamine (55%, 21/38; P = 0.002), although treating clinicians did not attribute any deaths to ketamine or midazolam.

CONCLUSIONS

Among children and neonates with RSE, ketamine was more often followed by seizure termination and less often associated with adverse effects than midazolam when administered as the first-line anesthetic infusion. Further prospective data are needed to compare first-line anesthetics for RSE.

Molecular and cellular context influences SCN8A variant function

Pathogenic variants in SCN8A, which encodes the voltage-gated sodium (NaV) channel NaV1.6, associate with neurodevelopmental disorders, including developmental and epileptic encephalopathy. Previous approaches to determine SCN8A variant function may be confounded by use of a neonatally expressed, alternatively spliced isoform of NaV1.6 (NaV1.6N) and engineered mutations rendering the channel tetrodotoxin (TTX) resistant. We investigated the impact of SCN8A alternative splicing on variant function by comparing the functional attributes of 15 variants expressed in 2 developmentally regulated splice isoforms (NaV1.6N, NaV1.6A). We employed automated patch clamp recording to enhance throughput, and developed a neuronal cell line (ND7/LoNav) with low levels of endogenous NaV current to obviate the need for TTX-resistance mutations. Expression of NaV1.6N or NaV1.6A in ND7/LoNav cells generated NaV currents with small, but significant, differences in voltage dependence of activation and inactivation. TTX-resistant versions of both isoforms exhibited significant functional differences compared with the corresponding WT channels. We demonstrated that many of the 15 disease-associated variants studied exhibited isoform-dependent functional effects, and that many of the studied SCN8A variants exhibited functional properties that were not easily classified as either gain- or loss-of-function. Our work illustrates the value of considering molecular and cellular context when investigating SCN8A variants.

The clinical and genetic spectrum of paediatric speech and language disorders in 52,143 individuals

Speech and language disorders are known to have a substantial genetic contribution. Although frequently examined as components of other conditions, research on the genetic basis of linguistic differences as separate phenotypic subgroups has been limited so far. Here, we performed an in-depth characterization of speech and language disorders in 52,143 individuals, reconstructing clinical histories using a large-scale data mining approach of the Electronic Medical Records (EMR) from an entire large paediatric healthcare network. The reported frequency of these disorders was the highest between 2 and 5 years old and spanned a spectrum of twenty-six broad speech and language diagnoses. We used Natural Language Processing to assess to which degree clinical diagnosis in full-text notes were reflected in ICD-10 diagnosis codes. We found that aphasia and speech apraxia could be easily retrieved through ICD-10 diagnosis codes, while stuttering as a speech phenotype was only coded in 12% of individuals through appropriate ICD-10 codes. We found significant comorbidity of speech and language disorders in neurodevelopmental conditions (30.31%) and to a lesser degree with epilepsies (6.07%) and movement disorders (2.05%). The most common genetic disorders retrievable in our EMR analysis were STXBP1 (n=21), PTEN (n=20), and CACNA1A (n=18). When assessing associations of genetic diagnoses with specific linguistic phenotypes, we observed associations of STXBP1 and aphasia (P=8.57 × 10-7, CI=18.62-130.39) and MYO7A with speech and language development delay due to hearing loss (P=1.24 × 10-5, CI=17.46-Inf). Finally, in a sub-cohort of 726 individuals with whole exome sequencing data, we identified an enrichment of rare variants in synaptic protein and neuronal receptor pathways and associations of UQCRC1 with expressive aphasia and WASHC4 with abnormality of speech or vocalization. In summary, our study outlines the landscape of paediatric speech and language disorders, confirming the phenotypic complexity of linguistic traits and novel genotype-phenotype associations. Subgroups of paediatric speech and language disorders differ significantly with respect to the composition of monogenic aetiologies.

De novo GRIN variants in M3 helix associated with neurological disorders control channel gating of NMDA receptor

N-methyl-D-aspartate receptors (NMDARs) are members of the glutamate receptor family and participate in excitatory postsynaptic transmission throughout the central nervous system. Genetic variants in GRIN genes encoding NMDAR subunits are associated with a spectrum of neurological disorders. The M3 transmembrane helices of the NMDAR couple directly to the agonist-binding domains and form a helical bundle crossing in the closed receptors that occludes the pore. The M3 functions as a transduction element whose conformational change couples ligand binding to opening of an ion conducting pore. In this study, we report the functional consequences of 48 de novo missense variants in GRIN1, GRIN2A, and GRIN2B that alter residues in the M3 transmembrane helix. These de novo variants were identified in children with neurological and neuropsychiatric disorders including epilepsy, developmental delay, intellectual disability, hypotonia and attention deficit hyperactivity disorder. All 48 variants in M3 for which comprehensive testing was completed produce a gain-of-function (28/48) compared to loss-of-function (9/48); 11 variants had an indeterminant phenotype. This supports the idea that a key structural feature of the M3 gate exists to stabilize the closed state so that agonist binding can drive channel opening. Given that most M3 variants enhance channel gating, we assessed the potency of FDA-approved NMDAR channel blockers on these variant receptors. These data provide new insight into the structure-function relationship of the NMDAR gate, and suggest that variants within the M3 transmembrane helix produce a gain-of-function.

De novo variants in DENND5B cause a neurodevelopmental disorder

The Rab family of guanosine triphosphatases (GTPases) includes key regulators of intracellular transport and membrane trafficking targeting specific steps in exocytic, endocytic, and recycling pathways. DENND5B (Rab6-interacting Protein 1B-like protein, R6IP1B) is the longest isoform of DENND5, an evolutionarily conserved DENN domain-containing guanine nucleotide exchange factor (GEF) that is highly expressed in the brain. Through exome sequencing and international matchmaking platforms, we identified five de novo variants in DENND5B in a cohort of five unrelated individuals with neurodevelopmental phenotypes featuring cognitive impairment, dysmorphism, abnormal behavior, variable epilepsy, white matter abnormalities, and cortical gyration defects. We used biochemical assays and confocal microscopy to assess the impact of DENND5B variants on protein accumulation and distribution. Then, exploiting fluorescent lipid cargoes coupled to high-content imaging and analysis in living cells, we investigated whether DENND5B variants affected the dynamics of vesicle-mediated intracellular transport of specific cargoes. We further generated an in silico model to investigate the consequences of DENND5B variants on the DENND5B-RAB39A interaction. Biochemical analysis showed decreased protein levels of DENND5B mutants in various cell types. Functional investigation of DENND5B variants revealed defective intracellular vesicle trafficking, with significant impairment of lipid uptake and distribution. Although none of the variants affected the DENND5B-RAB39A interface, all were predicted to disrupt protein folding. Overall, our findings indicate that DENND5B variants perturb intracellular membrane trafficking pathways and cause a complex neurodevelopmental syndrome with variable epilepsy and white matter involvement.

Early life seizures and epileptic spasms in STXBP1-related disorders

OBJECTIVE

Individuals with disease-causing variants in STXBP1 frequently have epilepsy onset in the first year of life with a variety of seizure types, including epileptic spasms. However, the impact of early onset seizures and antiseizure medication (ASM) on the risk of developing epileptic spasms and impact on their trajectory are poorly understood, limiting informed and anticipatory treatment, as well as trial design.

METHODS

We retrospectively reconstructed seizure and medication histories in weekly intervals for individuals with STXBP1 developmental and epileptic encephalopathy (DEE) with epilepsy onset in the first year of life and quantitatively analyzed longitudinal seizure histories and medication response.

RESULTS

We included 61 individuals with early onset seizures, 29 of whom had epileptic spasms. Individuals with neonatal seizures were likely to have continued seizures after the neonatal period (25/26). The risk of developing epileptic spasms was not increased in individuals with neonatal seizures or early infantile seizures (21/41 vs. 8/16, odds ratio [OR] = 1, 95% confidence interval [CI] = .3-3.9, p = 1). We did not find any ASM associated with the development of epileptic spasms following prior seizures. Individuals with prior seizures (n = 16/21, 76%) had a higher risk of developing refractory epileptic spasms (n = 5/8, 63%, OR = 1.9, 95% CI = .2-14.6, p = .6). Individuals with refractory epileptic spasms had a later onset of epileptic spasms (n = 20, median = 20 weeks) compared to individuals with nonrefractory epileptic spasms (n = 8, median = 13 weeks, p = .08).

SIGNIFICANCE

We provide a comprehensive assessment of early onset seizures in STXBP1-DEE and show that the risk of epileptic spasms is not increased following a prior history of early life seizures, nor by certain ASMs. Our study provides baseline information for targeted treatment and prognostication in early life seizures in STXBP1-DEE.

Biallelic CRELD1 variants cause a multisystem syndrome, including neurodevelopmental phenotypes, cardiac dysrhythmias, and frequent infections

PURPOSE

We sought to delineate a multisystem disorder caused by recessive cysteine-rich with epidermal growth factor-like domains 1 (CRELD1) gene variants.

METHODS

The impact of CRELD1 variants was characterized through an international collaboration utilizing next-generation DNA sequencing, gene knockdown, and protein overexpression in Xenopus tropicalis, and in vitro analysis of patient immune cells.

RESULTS

Biallelic variants in CRELD1 were found in 18 participants from 14 families. Affected individuals displayed an array of phenotypes involving developmental delay, early-onset epilepsy, and hypotonia, with about half demonstrating cardiac arrhythmias and some experiencing recurrent infections. Most harbored a frameshift in trans with a missense allele, with 1 recurrent variant, p.(Cys192Tyr), identified in 10 families. X tropicalis tadpoles with creld1 knockdown displayed developmental defects along with increased susceptibility to induced seizures compared with controls. Additionally, human CRELD1 harboring missense variants from affected individuals had reduced protein function, indicated by a diminished ability to induce craniofacial defects when overexpressed in X tropicalis. Finally, baseline analyses of peripheral blood mononuclear cells showed similar proportions of immune cell subtypes in patients compared with healthy donors.

CONCLUSION

This patient cohort, combined with experimental data, provide evidence of a multisystem clinical syndrome mediated by recessive variants in CRELD1.

Molecular and Phenotypic Characterization of the RORB-Related Disorder

BACKGROUND AND OBJECTIVES

Heterozygous variants in RAR-related orphan receptor B (RORB) have recently been associated with susceptibility to idiopathic generalized epilepsy. However, few reports have been published so far describing pathogenic variants of this gene in patients with epilepsy and intellectual disability (ID). In this study, we aimed to delineate the epilepsy phenotype associated with RORB pathogenic variants and to provide arguments in favor of the pathogenicity of variants.

METHODS

Through an international collaboration, we analyzed seizure characteristics, EEG data, and genotypes of a cohort of patients with heterozygous variants in RORB. To gain insight into disease mechanisms, we performed ex vivo cortical electroporation in mouse embryos of 5 selected variants, 2 truncating and 3 missense, and evaluated on expression and quantified changes in axonal morphology.

RESULTS

We identified 35 patients (17 male, median age 10 years, range 2.5-23 years) carrying 32 different heterozygous variants in RORB, including 28 single-nucleotide variants or small insertions/deletions (12 missense, 12 frameshift or nonsense, 2 splice-site variants, and 2 in-frame deletions), and 4 microdeletions; de novo in 18 patients and inherited in 10. Seizures were reported in 31/35 (89%) patients, with a median age at onset of 3 years (range 4 months-12 years). Absence seizures occurred in 25 patients with epilepsy (81%). Nineteen patients experienced a single seizure type: absences, myoclonic absences, or absences with eyelid myoclonia and focal seizures. Nine patients had absence seizures combined with other generalized seizure types. One patient had presented with absences associated with photosensitive occipital seizures. Three other patients had generalized tonic-clonic seizures without absences. ID of variable degree was observed in 85% of the patients. Expression studies in cultured neurons showed shorter axons for the 5 tested variants, both truncating and missense variants, supporting an impaired protein function.

DISCUSSION

In most patients, the phenotype of the RORB-related disorder associates absence seizures with mild-to-moderate ID. In silico and in vitro evaluation of the variants in our cohort, including axonal morphogenetic experiments in cultured neurons, supports their pathogenicity, showing a hypomorphic effect.

A structurally precise mechanism links an epilepsy-associated KCNC2 potassium channel mutation to interneuron dysfunction

De novo heterozygous variants in KCNC2 encoding the voltage-gated potassium (K+) channel subunit Kv3.2 are a recently described cause of developmental and epileptic encephalopathy (DEE). A de novo variant in KCNC2 c.374G > A (p.Cys125Tyr) was identified via exome sequencing in a patient with DEE. Relative to wild-type Kv3.2, Kv3.2-p.Cys125Tyr induces K+ currents exhibiting a large hyperpolarizing shift in the voltage dependence of activation, accelerated activation, and delayed deactivation consistent with a relative stabilization of the open conformation, along with increased current density. Leveraging the cryogenic electron microscopy (cryo-EM) structure of Kv3.1, molecular dynamic simulations suggest that a strong π-π stacking interaction between the variant Tyr125 and Tyr156 in the α-6 helix of the T1 domain promotes a relative stabilization of the open conformation of the channel, which underlies the observed gain of function. A multicompartment computational model of a Kv3-expressing parvalbumin-positive cerebral cortex fast-spiking γ-aminobutyric acidergic (GABAergic) interneuron (PV-IN) demonstrates how the Kv3.2-Cys125Tyr variant impairs neuronal excitability and dysregulates inhibition in cerebral cortex circuits to explain the resulting epilepsy.

The Human Phenotype Ontology in 2024: phenotypes around the world

The Human Phenotype Ontology (HPO) is a widely used resource that comprehensively organizes and defines the phenotypic features of human disease, enabling computational inference and supporting genomic and phenotypic analyses through semantic similarity and machine learning algorithms. The HPO has widespread applications in clinical diagnostics and translational research, including genomic diagnostics, gene-disease discovery, and cohort analytics. In recent years, groups around the world have developed translations of the HPO from English to other languages, and the HPO browser has been internationalized, allowing users to view HPO term labels and in many cases synonyms and definitions in ten languages in addition to English. Since our last report, a total of 2239 new HPO terms and 49235 new HPO annotations were developed, many in collaboration with external groups in the fields of psychiatry, arthrogryposis, immunology and cardiology. The Medical Action Ontology (MAxO) is a new effort to model treatments and other measures taken for clinical management. Finally, the HPO consortium is contributing to efforts to integrate the HPO and the GA4GH Phenopacket Schema into electronic health records (EHRs) with the goal of more standardized and computable integration of rare disease data in EHRs.

Spliceosome malfunction causes neurodevelopmental disorders with overlapping features

Pre-mRNA splicing is a highly coordinated process. While its dysregulation has been linked to neurological deficits, our understanding of the underlying molecular and cellular mechanisms remains limited. We implicated pathogenic variants in U2AF2 and PRPF19, encoding spliceosome subunits in neurodevelopmental disorders (NDDs), by identifying 46 unrelated individuals with 23 de novo U2AF2 missense variants (including 7 recurrent variants in 30 individuals) and 6 individuals with de novo PRPF19 variants. Eight U2AF2 variants dysregulated splicing of a model substrate. Neuritogenesis was reduced in human neurons differentiated from human pluripotent stem cells carrying two U2AF2 hyper-recurrent variants. Neural loss of function (LoF) of the Drosophila orthologs U2af50 and Prp19 led to lethality, abnormal mushroom body (MB) patterning, and social deficits, which were differentially rescued by wild-type and mutant U2AF2 or PRPF19. Transcriptome profiling revealed splicing substrates or effectors (including Rbfox1, a third splicing factor), which rescued MB defects in U2af50-deficient flies. Upon reanalysis of negative clinical exomes followed by data sharing, we further identified 6 patients with NDD who carried RBFOX1 missense variants which, by in vitro testing, showed LoF. Our study implicates 3 splicing factors as NDD-causative genes and establishes a genetic network with hierarchy underlying human brain development and function.

STXBP1: fast-forward to a brighter future - a patient organization perspective

Syntaxin-binding protein 1 related disorder (STXBP1-RD) is a rare neurologic disorder associated with global neurodevelopmental delay, intellectual disability, early-onset epilepsy, motor abnormalities, and autism. The underlying pathophysiology stems from a de novo mutation in the STXBP1 gene, which codes for the STXBP1 protein. The STXBP1 protein is involved in synaptic vesicle fusion and neurotransmitter release. Pathogenic variants in the STXBP1 gene generally result in haploinsufficiency, an impairment in neurotransmitter release, and subsequent dysfunction in neuronal communication. The STXBP1 Foundation was founded in 2017 to support families of children with STXBP1-RD and accelerate the development of effective therapies and, ultimately, a cure for the disorder. The Foundation initially supported research aimed at better understanding the complex phenotypic presentation of the disease as well as the development of animal and cellular models usable by the research community to more fully characterize STXBP1 function and disease pathogenicity. In 2023, the Foundation embarked on its STXBP1 Fast Forward Strategic Plan, which includes a prospective natural history study and substantive biomarker work to drive forward the development of new precision therapies for STXBP1-RD.

Genetic variants in DDX53 contribute to Autism Spectrum Disorder associated with the Xp22.11 locus

Autism Spectrum Disorder (ASD) exhibits an ~4:1 male-to-female sex bias and is characterized by early-onset impairment of social/communication skills, restricted interests, and stereotyped behaviors. Disruption of the Xp22.11 locus has been associated with ASD in males. This locus includes the three-exon PTCHD1 gene, an adjacent multi-isoform long noncoding RNA (lncRNA) named PTCHD1-AS (spanning ~1Mb), and a poorly characterized single-exon RNA helicase named DDX53 that is intronic to PTCHD1-AS. While the relationship between PTCHD1/PTCHD1-AS and ASD is being studied, the role of DDX53 has not been examined, in part because there is no apparent functional murine orthologue. Through clinical testing, here, we identified 6 males and 1 female with ASD from 6 unrelated families carrying rare, predicted-damaging or loss-of-function variants in DDX53. Then, we examined databases, including the Autism Speaks MSSNG and Simons Foundation Autism Research Initiative, as well as population controls. We identified 24 additional individuals with ASD harboring rare, damaging DDX53 variations, including the same variants detected in two families from the original clinical analysis. In this extended cohort of 31 participants with ASD (28 male, 3 female), we identified 25 mostly maternally-inherited variations in DDX53, including 18 missense changes, 2 truncating variants, 2 in-frame variants, 2 deletions in the 3' UTR and 1 copy number deletion. Our findings in humans support a direct link between DDX53 and ASD, which will be important in clinical genetic testing. These same autism-related findings, coupled with the observation that a functional orthologous gene is not found in mouse, may also influence the design and interpretation of murine-modelling of ASD.

Stuttering associated with a pathogenic variant in the chaperone protein cyclophilin 40

Stuttering is a common speech disorder that interrupts speech fluency and tends to cluster in families. Typically, stuttering is characterized by speech sounds, words or syllables which may be repeated or prolonged and speech that may be further interrupted by hesitations or 'blocks'. Rare variants in a small number of genes encoding lysosomal pathway proteins have been linked to stuttering. We studied a large four-generation family in which persistent stuttering was inherited in an autosomal dominant manner with disruption of the cortico-basal-ganglia-thalamo-cortical network found on imaging. Exome sequencing of three affected family members revealed the PPID c.808C>T (p.Pro270Ser) variant that segregated with stuttering in the family. We generated a Ppid p.Pro270Ser knock-in mouse model and performed ex vivo imaging to assess for brain changes. Diffusion-weighted MRI in the mouse revealed significant microstructural changes in the left corticospinal tract, as previously implicated in stuttering. Quantitative susceptibility mapping also detected changes in cortico-striatal-thalamo-cortical loop tissue composition, consistent with findings in affected family members. This is the first report to implicate a chaperone protein in the pathogenesis of stuttering. The humanized Ppid murine model recapitulates network findings observed in affected family members.

Delineating clinical and developmental outcomes in STXBP1-related disorders

STXBP1-related disorders are among the most common genetic epilepsies and neurodevelopmental disorders. However, the longitudinal epilepsy course and developmental end points, have not yet been described in detail, which is a critical prerequisite for clinical trial readiness. Here, we assessed 1281 cumulative patient-years of seizure and developmental histories in 162 individuals with STXBP1-related disorders and established a natural history framework. STXBP1-related disorders are characterized by a dynamic pattern of seizures in the first year of life and high variability in neurodevelopmental trajectories in early childhood. Epilepsy onset differed across seizure types, with 90% cumulative onset for infantile spasms by 6 months and focal-onset seizures by 27 months of life. Epilepsy histories diverged between variant subgroups in the first 2 years of life, when individuals with protein-truncating variants and deletions in STXBP1 (n = 39) were more likely to have infantile spasms between 5 and 6 months followed by seizure remission, while individuals with missense variants (n = 30) had an increased risk for focal seizures and ongoing seizures after the first year. Developmental outcomes were mapped using milestone acquisition data in addition to standardized assessments including the Gross Motor Function Measure-66 Item Set and the Grasping and Visual-Motor Integration subsets of the Peabody Developmental Motor Scales. Quantification of end points revealed high variability during the first 5 years of life, with emerging stratification between clinical subgroups. An earlier epilepsy onset was associated with lower developmental abilities, most prominently when assessing gross motor development and expressive communication. We found that individuals with neonatal seizures or early infantile seizures followed by seizure offset by 12 months of life had more predictable seizure trajectories in early to late childhood compared to individuals with more severe seizure presentations, including individuals with refractory epilepsy throughout the first year. Characterization of anti-seizure medication response revealed age-dependent response over time, with phenobarbital, levetiracetam, topiramate and adrenocorticotropic hormone effective in reducing seizures in the first year of life, while clobazam and the ketogenic diet were effective in long-term seizure management. Virtual clinical trials using seizure frequency as the primary outcome resulted in wide range of trial success probabilities across the age span, with the highest probability in early childhood between 1 year and 3.5 years. In summary, we delineated epilepsy and developmental trajectories in STXBP1-related disorders using standardized measures, providing a foundation to interpret future therapeutic strategies and inform rational trial design.

Molecular and Cellular Context Influences SCN8A Variant Function

Pathogenic variants in SCN8A , which encodes the voltage-gated sodium (Na V ) channel Na V 1.6, are associated with neurodevelopmental disorders including epileptic encephalopathy. Previous approaches to determine SCN8A variant function may be confounded by the use of a neonatal-expressed alternatively spliced isoform of Na V 1.6 (Na V 1.6N), and engineered mutations to render the channel tetrodotoxin (TTX) resistant. In this study, we investigated the impact of SCN8A alternative splicing on variant function by comparing the functional attributes of 15 variants expressed in two developmentally regulated splice isoforms (Na V 1.6N, Na V 1.6A). We employed automated patch clamp recording to enhance throughput, and developed a novel neuronal cell line (ND7/LoNav) with low levels of endogenous Na V current to obviate the need for TTX-resistance mutations. Expression of Na V 1.6N or Na V 1.6A in ND7/LoNav cells generated Na V currents that differed significantly in voltage-dependence of activation and inactivation. TTX-resistant versions of both isoforms exhibited significant functional differences compared to the corresponding wild-type (WT) channels. We demonstrated that many of the 15 disease-associated variants studied exhibited isoform-dependent functional effects, and that many of the studied SCN8A variants exhibited functional properties that were not easily classified as either gain- or loss-of-function. Our work illustrates the value of considering molecular and cellular context when investigating SCN8A variants.

A quality improvement initiative to improve folic acid supplementation counseling for adolescent females with epilepsy

OBJECTIVE

We designed a quality improvement (QI) project to improve rates of documented folic acid supplementation counseling for adolescent females with epilepsy, consistent with a quality measure from the American Academy of Neurology and American Epilepsy Society. Our SMART aim was to increase the percentage of visits at which folic acid counseling was addressed from our baseline rate of 23% to 50% by July 1, 2020.

METHODS

This initiative was conducted in female patients ≥12 years old with epilepsy who were prescribed daily antiseizure medication and were seen by the 13 providers in our Neurology QI Program. Using provider interviews, we undertook a root cause analysis of low counseling rates and identified the following main factors: insufficient time during clinic visit to counsel, lack of provider knowledge, and forgetting to counsel. Countermeasures were designed to address these main root causes and were implemented through iterative plan-do-study-act (PDSA) cycles. Interventions included provider education and features within the electronic health record, which were introduced sequentially, culminating in the creation of a best practice advisory (BPA). We performed biweekly chart reviews of visits for applicable patients to establish baseline performance rate and track progress over time. We used a statistical process control p-chart to analyze the outcome measure of documented counseling. As a balancing measure, clinicians were surveyed using the Technology Adoption Model survey to assess acceptance of the BPA.

RESULTS

From September 2019 to August 2022, the QI team improved rates of documented folic acid counseling from 23% to 73% through several PDSA cycles. This level of performance has been sustained over time. The most successful and sustainable intervention was the BPA. Provider acceptance of the BPA was overall positive.

SIGNIFICANCE

We successfully used QI methodology to improve and sustain our rates of documented folic acid supplementation counseling for adolescent females with epilepsy.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Loss of Grin2a causes a transient delay in the electrophysiological maturation of hippocampal parvalbumin interneurons

N-methyl-D-aspartate receptors (NMDARs) are ligand-gated ionotropic glutamate receptors that mediate a calcium-permeable component to fast excitatory neurotransmission. NMDARs are heterotetrameric assemblies of two obligate GluN1 subunits (GRIN1) and two GluN2 subunits (GRIN2A-GRIN2D). Sequencing data shows that 43% (297/679) of all currently known NMDAR disease-associated genetic variants are within the GRIN2A gene, which encodes the GluN2A subunit. Here, we show that unlike missense GRIN2A variants, individuals affected with disease-associated null GRIN2A variants demonstrate a transient period of seizure susceptibility that begins during infancy and diminishes near adolescence. We show increased circuit excitability and CA1 pyramidal cell output in juvenile mice of both Grin2a+/- and Grin2a-/- mice. These alterations in somatic spiking are not due to global upregulation of most Grin genes (including Grin2b). Deeper evaluation of the developing CA1 circuit led us to uncover age- and Grin2a gene dosing-dependent transient delays in the electrophysiological maturation programs of parvalbumin (PV) interneurons. We report that Grin2a+/+ mice reach PV cell electrophysiological maturation between the neonatal and juvenile neurodevelopmental timepoints, with Grin2a+/- mice not reaching PV cell electrophysiological maturation until preadolescence, and Grin2a-/- mice not reaching PV cell electrophysiological maturation until adulthood. Overall, these data may represent a molecular mechanism describing the transient nature of seizure susceptibility in disease-associated null GRIN2A patients.

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

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

Natural History and Developmental Trajectories of Individuals With Disease-Causing Variants in STXBP1

BACKGROUND AND OBJECTIVES

Pathogenic variants in STXBP1 are among the major genetic causes of neurodevelopmental disorders. Despite the increasing number of individuals diagnosed without a history of epilepsy, little is known about the natural history and developmental trajectories in this subgroup and endpoints for future therapeutic studies are limited to seizure control.

METHODS

We performed a cross-sectional retrospective study using standardized questionnaires for clinicians and caregivers of individuals with STXBP1-related disorders capturing medical histories, genetic findings, and developmental outcomes. Motor and language function were assessed using Gross Motor Function Classification System (GMFCS) scores and a speech impairment score and were compared within and across clinically defined subgroups.

RESULTS

We collected data of 71 individuals with STXBP1-related disorders, including 44 previously unreported individuals. Median age at inclusion was 5.3 years (interquartile range 3.5-9.3) with the oldest individual aged 43.8 years. Epilepsy was absent in 18/71 (25%) of individuals. The range of developmental outcomes was broad, including 2 individuals presenting with close to age-appropriate motor development. Twenty-nine of 61 individuals (48%) were able to walk unassisted, and 24/69 (35%) were able to speak single words. Individuals without epilepsy presented with a similar onset and spectrum of phenotypic features but had lower GMFCS scores (median 3 vs 4, p < 0.01) than individuals with epilepsy. Individuals with epileptic spasms were less likely to walk unassisted than individuals with other seizure types (6% vs 58%, p < 0.01). Individuals with early epilepsy onset had higher speech impairment scores (p = 0.02) than individuals with later epilepsy onset.

DISCUSSION

We expand the spectrum of STXBP1-related disorders and provide clinical features and developmental trajectories in individuals with and without a history of epilepsy. Individuals with epilepsy, in particular epileptic spasms, and neonatal or early-onset presented with less favorable motor and language functional outcomes compared with individuals without epilepsy. These findings identify children at risk for severe disease and can serve as comparator for future interventional studies in STXBP1-related disorders.