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

Voltage-gated calcium channels in genetic epilepsies

Voltage-gated calcium channels (VGCC) are abundant in the central nervous system and serve a broad spectrum of functions, either directly in cellular excitability or indirectly to regulate Ca2+ homeostasis. Ca2+ ions act as one of the main connections in excitation-transcription coupling, muscle contraction and excitation-exocytosis coupling, including synaptic transmission. In recent years, many genes encoding VGCCs main α or additional auxiliary subunits have been associated with epilepsy. This review sums up the current state of knowledge on disease mechanisms and provides guidance on disease-specific therapies where applicable.

The role of candidate pharmacogenetic variants in determining valproic acid efficacy, toxicity and concentrations in patients with epilepsy

Background

Antiseizure medications (ASM) exhibit considerable interindividual variability in terms of efficacy and adverse events. Genetic variation is thought to contribute to these differences in clinical outcomes. Specifically, the response to valproic acid (VPA), a widely used ASM, is influenced by multiple pharmacogenetic factors. However, and in contrast to other ASMs such as phenytoin and carbamazepine, there is a paucity of data on the association between VPA and various gene variants. The aim of this study was hence to evaluate the influence of candidate pharmacogenetic variants on VPA efficacy, toxicity and serum concentrations in a homogeneous cohort of patients newly diagnosed with genetic generalized epilepsies (GGE).

Methods

In this prospective cohort study, demographic, clinical and treatment outcomes of GGE patients were retrieved from their medical records. Whole exome sequencing was performed in collaboration with Epi25. Gene variants associated with VPA efficacy, metabolism and toxicities were retrieved from PharmGKB. An analysis was then conducted to explore potential associations between these gene variants and VPA clinical outcomes.

Results

Of the 166 patients included, 60 (36.1%) experienced treatment failure while 106 (63.9%) achieved treatment success. After adjusting for VPA maintenance dose, carriers of the rs3892097 (CYP2D6) variant were 2.5 times more likely to experience treatment failure compared to wildtype (p = 0.026). The rs1057910 variant (CYP2C9*3) was associated with increased serum VPA concentrations (p = 0.034). Moreover, the rs1137101 variant (LEPR gene, a metabolism regulator) was significantly associated with a higher risk of weight gain (regression coefficient of 3.430 [0.674; 6.186], p = 0.015) and a higher frequency of hair loss (OR = 3.394 [1.157; 9.956], p = 0.026), while the rs4480 variant (SOD2 gene, encoding for a mitochondrial scavenging enzyme) was correlated with a lower frequency of hair loss (OR = 0.276 [0.089; 0.858], p = 0.026).

Conclusion

These findings highlight the role of genetic factors in VPA treatment and underscore the potential for developing therapeutic strategies to enhance patient outcomes and minimize adverse effects.

Exome sequencing of 20,979 individuals with epilepsy reveals shared and distinct ultra-rare genetic risk across disorder subtypes

Identifying genetic risk factors for highly heterogeneous disorders such as epilepsy remains challenging. Here we present, to our knowledge, the largest whole-exome sequencing study of epilepsy to date, with more than 54,000 human exomes, comprising 20,979 deeply phenotyped patients from multiple genetic ancestry groups with diverse epilepsy subtypes and 33,444 controls, to investigate rare variants that confer disease risk. These analyses implicate seven individual genes, three gene sets and four copy number variants at exome-wide significance. Genes encoding ion channels show strong association with multiple epilepsy subtypes, including epileptic encephalopathies and generalized and focal epilepsies, whereas most other gene discoveries are subtype specific, highlighting distinct genetic contributions to different epilepsies. Combining results from rare single-nucleotide/short insertion and deletion variants, copy number variants and common variants, we offer an expanded view of the genetic architecture of epilepsy, with growing evidence of convergence among different genetic risk loci on the same genes. Top candidate genes are enriched for roles in synaptic transmission and neuronal excitability, particularly postnatally and in the neocortex. We also identify shared rare variant risk between epilepsy and other neurodevelopmental disorders. Our data can be accessed via an interactive browser, hopefully facilitating diagnostic efforts and accelerating the development of follow-up studies.

Exome sequencing of 20,979 individuals with epilepsy reveals shared and distinct ultra-rare genetic risk across disorder subtypes

Identifying genetic risk factors for highly heterogeneous disorders like epilepsy remains challenging. Here, we present the largest whole-exome sequencing study of epilepsy to date, with >54,000 human exomes, comprising 20,979 deeply phenotyped patients from multiple genetic ancestry groups with diverse epilepsy subtypes and 33,444 controls, to investigate rare variants that confer disease risk. These analyses implicate seven individual genes, three gene sets, and four copy number variants at exome-wide significance. Genes encoding ion channels show strong association with multiple epilepsy subtypes, including epileptic encephalopathies, generalized and focal epilepsies, while most other gene discoveries are subtype-specific, highlighting distinct genetic contributions to different epilepsies. Combining results from rare single nucleotide/short indel-, copy number-, and common variants, we offer an expanded view of the genetic architecture of epilepsy, with growing evidence of convergence among different genetic risk loci on the same genes. Top candidate genes are enriched for roles in synaptic transmission and neuronal excitability, particularly postnatally and in the neocortex. We also identify shared rare variant risk between epilepsy and other neurodevelopmental disorders. Our data can be accessed via an interactive browser, hopefully facilitating diagnostic efforts and accelerating the development of follow-up studies.

Expansion of the neurodevelopmental phenotype of individuals with EEF1A2 variants and genotype-phenotype study

Translation elongation factor eEF1A2 constitutes the alpha subunit of the elongation factor-1 complex, responsible for the enzymatic binding of aminoacyl-tRNA to the ribosome. Since 2012, 21 pathogenic missense variants affecting EEF1A2 have been described in 42 individuals with a severe neurodevelopmental phenotype including epileptic encephalopathy and moderate to profound intellectual disability (ID), with neurological regression in some patients. Through international collaborative call, we collected 26 patients with EEF1A2 variants and compared them to the literature. Our cohort shows a significantly milder phenotype. 83% of the patients are walking (vs. 29% in the literature), and 84% of the patients have language skills (vs. 15%). Three of our patients do not have ID. Epilepsy is present in 63% (vs. 93%). Neurological examination shows a less severe phenotype with significantly less hypotonia (58% vs. 96%), and pyramidal signs (24% vs. 68%). Cognitive regression was noted in 4% (vs. 56% in the literature). Among individuals over 10 years, 56% disclosed neurocognitive regression, with a mean age of onset at 2 years. We describe 8 novel missense variants of EEF1A2. Modeling of the different amino-acid sites shows that the variants associated with a severe phenotype, and the majority of those associated with a moderate phenotype, cluster within the switch II region of the protein and thus may affect GTP exchange. In contrast, variants associated with milder phenotypes may impact secondary functions such as actin binding. We report the largest cohort of individuals with EEF1A2 variants thus far, allowing us to expand the phenotype spectrum and reveal genotype-phenotype correlations.

Haploinsufficiency underlies the neurodevelopmental consequences of SLC6A1 variants

Heterozygous variants in SLC6A1, encoding the GAT-1 GABA transporter, are associated with seizures, developmental delay, and autism. The majority of affected individuals carry missense variants, many of which are recurrent germline de novo mutations, raising the possibility of gain-of-function or dominant-negative effects. To understand the functional consequences, we performed an in vitro GABA uptake assay for 213 unique variants, including 24 control variants. De novo variants consistently resulted in a decrease in GABA uptake, in keeping with haploinsufficiency underlying all neurodevelopmental phenotypes. Where present, ClinVar pathogenicity reports correlated well with GABA uptake data; the functional data can inform future reports for the remaining 72% of unscored variants. Surface localization was assessed for 86 variants; two-thirds of loss-of-function missense variants prevented GAT-1 from being present on the membrane while GAT-1 was on the surface but with reduced activity for the remaining third. Surprisingly, recurrent de novo missense variants showed moderate loss-of-function effects that reduced GABA uptake with no evidence for dominant-negative or gain-of-function effects. Using linear regression across multiple missense severity scores to extrapolate the functional data to all potential SLC6A1 missense variants, we observe an abundance of GAT-1 residues that are sensitive to substitution. The extent of this missense vulnerability accounts for the clinically observed missense enrichment; overlap with hypermutable CpG sites accounts for the recurrent missense variants. Strategies to increase the expression of the wild-type SLC6A1 allele are likely to be beneficial across neurodevelopmental disorders, though the developmental stage and extent of required rescue remain unknown.

Allelic heterogeneity and abnormal vesicle recycling in PLAA-related neurodevelopmental disorders

The human PLAA gene encodes Phospholipase-A2-Activating-Protein (PLAA) involved in trafficking of membrane proteins. Through its PUL domain (PLAP, Ufd3p, and Lub1p), PLAA interacts with p97/VCP modulating synaptic vesicles recycling. Although few families carrying biallelic PLAA variants were reported with progressive neurodegeneration, consequences of monoallelic PLAA variants have not been elucidated. Using exome or genome sequencing we identified PLAA de-novo missense variants, affecting conserved residues within the PUL domain, in children affected with neurodevelopmental disorders (NDDs), including psychomotor regression, intellectual disability (ID) and autism spectrum disorders (ASDs). Computational and in-vitro studies of the identified variants revealed abnormal chain arrangements at C-terminal and reduced PLAA-p97/VCP interaction, respectively. These findings expand both allelic and phenotypic heterogeneity associated to PLAA-related neurological disorders, highlighting perturbed vesicle recycling as a potential disease mechanism in NDDs due to genetic defects of PLAA.

Gene-gene interaction network analysis indicates CNTN2 is a candidate gene for idiopathic generalized epilepsy

Twin and family studies have established the genetic contribution to idiopathic generalized epilepsy (IGE). The genetic architecture of IGE is generally complex and heterogeneous, and the majority of the genetic burden in IGE remains unsolved. We hypothesize that gene-gene interactions contribute to the complex inheritance of IGE. CNTN2 (OMIM* 615,400) variants have been identified in cases with familial adult myoclonic epilepsy and other epilepsies. To explore the gene-gene interaction network in IGE, we took the CNTN2 gene as an example and investigated its co-occurrent genetic variants in IGE cases. We performed whole-exome sequencing in 114 unrelated IGE cases and 296 healthy controls. Variants were qualified with sequencing quality, minor allele frequency, in silico prediction, genetic phenotype, and recurrent case numbers. The STRING_TOP25 gene interaction network analysis was introduced with the bait gene CNTN2 (denoted as A). The gene-gene interaction pair mode was presumed to be A + c, A + d, A + e, with a leading gene A, or A + B + f, A + B + g, A + B + h, with a double-gene A + B, or other combinations. We compared the number of gene interaction pairs between the case and control groups. We identified three pairs in the case group, CNTN2 + PTPN18, CNTN2 + CNTN1 + ANK2 + ANK3 + SNTG2, and CNTN2 + PTPRZ1, while we did not discover any pairs in the control group. The number of gene interaction pairs in the case group was much more than in the control group (p = 0.021). Taking together the genetic bioinformatics, reported epilepsy cases, and statistical evidence in the study, we supposed CNTN2 as a candidate pathogenic gene for IGE. The gene interaction network analysis might help screen candidate genes for IGE or other complex genetic disorders.

Clinical and molecular heterogeneity of VPS13D-related neurodevelopmental and movement disorders

BACKGROUND

The VPS13 family of proteins has been implicated in lipid transport and trafficking between endoplasmic reticulum and organelles, to maintain homeostasis of subcellular membranes. Recently, pathogenic variants in each human VPS13S gene, have been linked to distinct human neurodevelopmental or neurodegenerative disorders. Within the VPS13 family of genes, VPS13D is known to be implicated in mitochondria homeostasis and function.

METHODS

We investigated a Pakistani sibship affected with neurodevelopmental impairment and severe hyperkinetic (choreoathetoid) movements. Whole exome sequencing (WES) and Sanger sequencing were performed to identify potential candidate variants segregating in the family. We described clinical phenotypes and natural history of the disease during a 3-year clinical follow-up and summarized literature data related to previously identified patients with VPS13D-related neurological disorders.

RESULTS

We identified by WES an homozygous non-synonymous variant in VPS13D (c.5723 T > C; p.Ile1908Thr) as the potential underlying cause of the disease in our family. Two young siblings developed an early-onset neurological impairment characterized by global developmental delay, with impaired speech and motor milestones, associated to hyperkinetic movement disorders as well as progressive and non-progressive neurological abnormalities.

CONCLUSION

In this study we delineated the heterogeneity of VPS13D-related clinical phenotypes and described a novel VPS13D homozygous variant associated with severe neurological impairment. Further studies will be pivotal to understand the exact VPS13D function and its impact on mitochondria homeostasis, brain development and regulation of movements, to further clarify genotype-phenotype correlations and provide crucial prognostic information and potential therapeutic implications.

Delayed Graft Function in Pediatric Kidney Transplant: Risk Factors and Outcomes

OBJECTIVES

We aimed to identify risk factors and outcomes of delayed graft function in pediatric kidney transplant.

MATERIALS AND METHODS

This retrospective study included all kidney transplant recipients ≤19 years old followed up in our department for a period of 34 years, from January 1989 to December 2022.

RESULTS

We included 113 kidney transplant recipients. Delayed graft function occurred in 17 cases (15%). Posttransplant red blood cell transfusion was strongly associated with delayed graft function (adjusted odds ratio = 23.91; 95% CI, 2.889-197.915). Use of allografts with multiple arteries and cold ischemia time >20 hours were risk factors for delayed graft function (adjusted odds ratio = 52.51 and 49.4; 95% CI, 2.576-1070.407 and 1.833-1334.204, respectively). Sex-matched transplants and living donors were protective factors for delayed graft function (adjusted odds ratio = 0.043 and 0.027; 95% CI, 0.005-0.344 and 0.003-0.247, respectively). Total HLA mismatches <3 played a protective role for delayed graft function (adjusted odds ratio = 0.114; 95% CI, 0.020-0.662), whereas transplant within compatible but different blood types increased the risk of delayed graft function (adjusted odds ratio = 20.54; 95% CI, 1.960- 215.263). No significant correlation was shown between delayed graft function and allograft survival (P = .190). Our study suggested delayed graft function as a key factor in allograft rejection-free survival (adjusted odds ratio = 3.832; 95% CI, 1.186-12.377). Delayed graft function was a negative factor for early graft function; patients with delayed graft function had a lower estimated glomerular filtration rate at discharge (P = .024) and at 3 (P = .034), 6 (P = .019), and 12 months (P = .011) posttransplant.

CONCLUSIONS

Delayed graft function is a major determinant of early graft function and allograft rejection-free survival. Further research is required to establish proper preventive measures.

Neuroimaging in PRUNE1 syndrome: a mini-review of the literature

Prune exopolyphosphatase 1 (PRUNE1) is a short-chain phosphatase that is part of the aspartic acid-histidine-histidine (DHH) family of proteins. PRUNE1 is highly expressed in the central nervous system and is crucially involved in neurodevelopment, cytoskeletal rearrangement, cell migration, and proliferation. Recently, biallelic PRUNE1 variants have been identified in patients with neurodevelopmental disorders, hypotonia, microcephaly, variable cerebral anomalies, and other features. PRUNE1 hypomorphic mutations mainly affect the DHH1 domain, leading to an impactful decrease in enzymatic activity with a loss-of-function mechanism. In this review, we explored both the clinical and radiological spectrum related to PRUNE1 pathogenic variants described to date. Specifically, we focused on neuroradiological findings that, together with clinical phenotypes and genetic data, allow us to best characterize affected children with diagnostic and potential prognostic implications.

Neuroimaging features of WOREE syndrome: a mini-review of the literature

The WWOX gene encodes a 414-amino-acid protein composed of two N-terminal WW domains and a C-terminal short-chain dehydrogenase/reductase (SDR) domain. WWOX protein is highly conserved among species and mainly expressed in the cerebellum, cerebral cortex, brain stem, thyroid, hypophysis, and reproductive organs. It plays a crucial role in the biology of the central nervous system, and it is involved in neuronal development, migration, and proliferation. Biallelic pathogenic variants in WWOX have been associated with an early infantile epileptic encephalopathy known as WOREE syndrome. Both missense and null variants have been described in affected patients, leading to a reduction in protein function and stability. The most severe WOREE phenotypes have been related to biallelic null/null variants, associated with the complete loss of function of the protein. All affected patients showed brain anomalies on magnetic resonance imaging (MRI), suggesting the pivotal role of WWOX protein in brain homeostasis and developmental processes. We provided a literature review, exploring both the clinical and radiological spectrum related to WWOX pathogenic variants, described to date. We focused on neuroradiological findings to better delineate the WOREE phenotype with diagnostic and prognostic implications.

Unraveling the role of non-coding rare variants in epilepsy

The discovery of new variants has leveled off in recent years in epilepsy studies, despite the use of very large cohorts. Consequently, most of the heritability is still unexplained. Rare non-coding variants have been largely ignored in studies on epilepsy, although non-coding single nucleotide variants can have a significant impact on gene expression. We had access to whole genome sequencing (WGS) from 247 epilepsy patients and 377 controls. To assess the functional impact of non-coding variants, ExPecto, a deep learning algorithm was used to predict expression change in brain tissues. We compared the burden of rare non-coding deleterious variants between cases and controls. Rare non-coding highly deleterious variants were significantly enriched in Genetic Generalized Epilepsy (GGE), but not in Non-Acquired Focal Epilepsy (NAFE) or all epilepsy cases when compared with controls. In this study we showed that rare non-coding deleterious variants are associated with epilepsy, specifically with GGE. Larger WGS epilepsy cohort will be needed to investigate those effects at a greater resolution. Nevertheless, we demonstrated the importance of studying non-coding regions in epilepsy, a disease where new discoveries are scarce.

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.

Technological and computational approaches to detect somatic mosaicism in epilepsy

Lesional epilepsy is a common and severe disease commonly associated with malformations of cortical development, including focal cortical dysplasia and hemimegalencephaly. Recent advances in sequencing and variant calling technologies have identified several genetic causes, including both short/single nucleotide and structural somatic variation. In this review, we aim to provide a comprehensive overview of the methodological advancements in this field while highlighting the unresolved technological and computational challenges that persist, including ultra-low variant allele fractions in bulk tissue, low availability of paired control samples, spatial variability of mutational burden within the lesion, and the issue of false-positive calls and validation procedures. Information from genetic testing in focal epilepsy may be integrated into clinical care to inform histopathological diagnosis, postoperative prognosis, and candidate precision therapies.

Estimating clinical risk in gene regions from population sequencing cohort data

While pathogenic variants can significantly increase disease risk, it is still challenging to estimate the clinical impact of rare missense variants more generally. Even in genes such as BRCA2 or PALB2, large cohort studies find no significant association between breast cancer and rare missense variants collectively. Here, we introduce REGatta, a method to estimate clinical risk from variants in smaller segments of individual genes. We first define these regions by using the density of pathogenic diagnostic reports and then calculate the relative risk in each region by using over 200,000 exome sequences in the UK Biobank. We apply this method in 13 genes with established roles across several monogenic disorders. In genes with no significant difference at the gene level, this approach significantly separates disease risk for individuals with rare missense variants at higher or lower risk (BRCA2 regional model OR = 1.46 [1.12, 1.79], p = 0.0036 vs. BRCA2 gene model OR = 0.96 [0.85, 1.07] p = 0.4171). We find high concordance between these regional risk estimates and high-throughput functional assays of variant impact. We compare our method with existing methods and the use of protein domains (Pfam) as regions and find REGatta better identifies individuals at elevated or reduced risk. These regions provide useful priors and are potentially useful for improving risk assessment for genes associated with monogenic diseases.

Analysis of common genetic variation across targets of microRNAs dysregulated both in ASD and epilepsy reveals negative correlation

Genetic overlap involving rare disrupting mutations may contribute to high comorbidity rates between autism spectrum disorders and epilepsy. Despite their polygenic nature, genome-wide association studies have not reported a significant contribution of common genetic variation to comorbidity between both conditions. Analysis of common genetic variation affecting specific shared pathways such as miRNA dysregulation could help to elucidate the polygenic mechanisms underlying comorbidity between autism spectrum disorders and epilepsy. We evaluated here the role of common predisposing variation to autism spectrum disorders and epilepsy across target genes of 14 miRNAs selected through bibliographic research as being dysregulated in both disorders. We considered 4,581 target genes from various in silico sources. We described negative genetic correlation between autism spectrum disorders and epilepsy across variants located within target genes of the 14 miRNAs selected (p = 0.0228). Moreover, polygenic transmission disequilibrium test on an independent cohort of autism spectrum disorders trios (N = 233) revealed an under-transmission of autism spectrum disorders predisposing alleles within miRNAs’ target genes across autism spectrum disorders trios without comorbid epilepsy, thus reinforcing the negative relationship at the common genetic variation between both traits. Our study provides evidence of a negative relationship between autism spectrum disorders and epilepsy at the common genetic variation level that becomes more evident when focusing on the miRNA regulatory networks, which contrasts with observed clinical comorbidity and results from rare variation studies. Our findings may help to conceptualize the genetic heterogeneity and the comorbidity with epilepsy in autism spectrum disorders.

Novel GATOR1 variants in focal epilepsy

INTRODUCTION

Variants in GATOR1 genes are well established in focal epilepsy syndromes. A strong association of GATOR1 variants with drug-resistant epilepsy as well as an increased risk of sudden unexplained death in epilepsy warrants developing strategies to facilitate the identification of patients who could potentially benefit from genetic testing and precision medicine. We aimed to determine the yield of GATOR1 gene sequencing in patients with focal epilepsy typically referred for genetic testing, establish novel GATOR1 variants and determine clinical, electroencephalographic, and radiological characteristics of variant carriers.

PATIENTS AND METHODS

Ninety-six patients with clinical suspicion of genetic focal epilepsy with previous comprehensive diagnostic epilepsy evaluation in The Neurology Clinic, University Clinical Center of Serbia, were included in the study. Sequencing was performed using a custom gene panel encompassing DEPDC5, NPRL2, and NPRL3. Variants of interest (VOI) were classified according to criteria proposed by the American College of Medical Genetics and the Association for Molecular Pathology.

RESULTS

Four previously unreported VOI in 4/96 (4.2%) patients were found in our cohort. Three likely pathogenic variants were determined in 3/96 (3.1%) patients, one frameshift variant in DEPDC5 in a patient with nonlesional frontal lobe epilepsy, one splicogenic DEPDC5 variant in a patient with nonlesional posterior quadrant epilepsy, and one frameshift variant in NPRL2 in a patient with temporal lobe epilepsy associated with hippocampal sclerosis. Only one VOI, a missense variant in NPRL3, found in 1/96 (1.1%) patients, was classified as a variant of unknown significance.

CONCLUSION

GATOR1 gene sequencing was diagnostic in 3.1% of our cohort and revealed three novel likely pathogenic variants, including a previously unreported association of temporal lobe epilepsy with hippocampal sclerosis with an NPRL2 variant. Further research is essential for a better understanding of the clinical scope of GATOR1 gene-associated epilepsy.

Estimating clinical risk in gene regions from population sequencing cohort data

While pathogenic variants significantly increase disease risk in many genes, it is still challenging to estimate the clinical impact of rare missense variants more generally. Even in genes such as BRCA2 or PALB2 , large cohort studies find no significant association between breast cancer and rare germline missense variants collectively. Here we introduce REGatta, a method to improve the estimation of clinical risk in gene segments. We define gene regions using the density of pathogenic diagnostic reports, and then calculate the relative risk in each of these regions using 109,581 exome sequences from women in the UK Biobank. We apply this method in seven established breast cancer genes, and identify regions in each gene with statistically significant differences in breast cancer incidence for rare missense carriers. Even in genes with no significant difference at the gene level, this approach significantly separates rare missense variant carriers at higher or lower risk ( BRCA2 regional model OR=1.46 [1.12, 1.79], p=0.0036 vs. BRCA2 gene model OR=0.96 [0.85,1.07] p=0.4171). We find high concordance between these regional risk estimates and high-throughput functional assays of variant impact. We compare with existing methods and the use of protein domains (Pfam) as regions, and find REGatta better identifies individuals at elevated or reduced risk. These regions provide useful priors which can potentially be used to improve risk assessment and clinical management.

Risk Variants in the Exomes of Children With Critical Illness

IMPORTANCE

Diagnostic genetic testing can lead to changes in management in the pediatric intensive care unit. Genetic risk in children with critical illness but nondiagnostic exome sequencing (ES) has not been explored.

OBJECTIVE

To assess the association between loss-of-function (LOF) variants and pediatric critical illness.

DESIGN, SETTING, AND PARTICIPANTS

This genetic association study examined ES first screened for causative variants among 267 children at the Morgan Stanley Children's Hospital of NewYork-Presbyterian, of whom 22 were otherwise healthy with viral respiratory failure; 18 deceased children with bronchiolitis from the Office of the Chief Medical Examiner of New York City, of whom 14 were previously healthy; and 9990 controls from the Institute for Genomic Medicine at Columbia University Irving Medical Center. The ES data were generated between January 1, 2015, and December 31, 2020, and analyzed between January 1, 2017, and September 2, 2022.

EXPOSURE

Critical illness.

MAIN OUTCOMES AND MEASURES

Odds ratios and P values for genes and gene-sets enriched for rare LOF variants and the loss-of-function observed/expected upper bound fraction (LOEUF) score at which cases have a significant enrichment.

RESULTS

This study included 285 children with critical illness (median [range] age, 4.1 [0-18.9] years; 148 [52%] male) and 9990 controls. A total of 228 children (80%) did not receive a genetic diagnosis. After quality control (QC), 231 children harbored excess rare LOF variants in genes with a LOEUF score of 0.680 or less (intolerant genes) (P = 1.0 × 10-5). After QC, 176 children without a diagnosis harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 1.8; 95% CI, 1.3-2.5). After QC, 25 children with viral respiratory failure harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 2.8; 95% CI, 1.1-6.6). A total of 114 undiagnosed children were enriched for de novo LOF variants in genes without a known disease association (observed, 14; expected, 6.8; enrichment, 2.05).

CONCLUSIONS AND RELEVANCE

In this genetic association study, excess LOF variants were observed among critically ill children despite nondiagnostic ES. Variants lay in genes without a known disease association, suggesting future investigation may connect phenotypes to causative genes.

Sporadic hypothalamic hamartoma is a ciliopathy with somatic and bi-allelic contributions

Hypothalamic hamartoma with gelastic seizures is a well-established cause of drug-resistant epilepsy in early life. The development of novel surgical techniques has permitted the genomic interrogation of hypothalamic hamartoma tissue. This has revealed causative mosaic variants within GLI3, OFD1 and other key regulators of the sonic-hedgehog pathway in a minority of cases. Sonic-hedgehog signalling proteins localize to the cellular organelle primary cilia. We therefore explored the hypothesis that cilia gene variants may underlie hitherto unsolved cases of sporadic hypothalamic hamartoma. We performed high-depth exome sequencing and chromosomal microarray on surgically resected hypothalamic hamartoma tissue and paired leukocyte-derived DNA from 27 patients. We searched for both germline and somatic variants under both dominant and bi-allelic genetic models. In hamartoma-derived DNA of seven patients we identified bi-allelic (one germline, one somatic) variants within one of four cilia genes-DYNC2I1, DYNC2H1, IFT140 or SMO. In eight patients, we identified single somatic variants in the previously established hypothalamic hamartoma disease genes GLI3 or OFD1. Overall, we established a plausible molecular cause for 15/27 (56%) patients. Here, we expand the genetic architecture beyond single variants within dominant disease genes that cause sporadic hypothalamic hamartoma to bi-allelic (one germline/one somatic) variants, implicate three novel cilia genes and reconceptualize the disorder as a ciliopathy.

Whole genome sequences discriminate hereditary hemorrhagic telangiectasia phenotypes by non-HHT deleterious DNA variation

The abnormal vascular structures of hereditary hemorrhagic telangiectasia (HHT) often cause severe anemia due to recurrent hemorrhage, but HHT causal genes do not predict the severity of hematological complications. We tested for chance inheritance and clinical associations of rare deleterious variants in which loss-of-function causes bleeding or hemolytic disorders in the general population. In double-blinded analyses, all 104 patients with HHT from a single reference center recruited to the 100 000 Genomes Project were categorized on new MALO (more/as-expected/less/opposite) sub-phenotype severity scales, and whole genome sequencing data were tested for high impact variants in 75 HHT-independent genes encoding coagulation factors, or platelet, hemoglobin, erythrocyte enzyme, and erythrocyte membrane constituents. Rare variants (all gnomAD allele frequencies <0.003) were identified in 56 (75%) of these 75 HHT-unrelated genes. Deleteriousness assignments by Combined Annotation Dependent Depletion (CADD) scores >15 were supported by gene-level mutation significance cutoff scores. CADD >15 variants were identified in 38/104 (36.5%) patients with HHT, found for 1 in 10 patients within platelet genes; 1 in 8 within coagulation genes; and 1 in 4 within erythrocyte hemolytic genes. In blinded analyses, patients with greater hemorrhagic severity that had been attributed solely to HHT vessels had more CADD-deleterious variants in platelet (Spearman ρ = 0.25; P = .008) and coagulation (Spearman ρ = 0.21; P = .024) genes. However, the HHT cohort had 60% fewer deleterious variants in platelet and coagulation genes than expected (Mann-Whitney test P = .021). In conclusion, patients with HHT commonly have rare variants in genes of relevance to their phenotype, offering new therapeutic targets and opportunities for informed, personalized medicine strategies.

Recurrent de novo single point variant on the gene encoding Na+/K+ pump results in epilepsy

The etiology of epilepsy remains undefined in two-thirds of patients. Here, we identified a de novo variant of ATP1A2 (c.2426 T > G, p.Leu809Arg), which encodes the α2 subunit of Na⁺/K⁺-ATPase^, from a family with idiopathic epilepsy. This variant caused epilepsy with hemiplegic migraine in the study patients. We generated the point variant mouse model Atp1a2^L809R, which recapitulated the epilepsy observed in the study patients. In Atp1a2^L809R/WT mice, convulsions were observed and cognitive and memory function was impaired. This variant affected the potassium binding function of the protein, disabling its ion transport ability, thereby increasing the frequency of nerve impulses. Valproate (VPA) and Carbamazepine (CBZ) have limited therapeutic efficacy in ameliorating the epileptic syndromes of Atp1a2^L809R/WT mice. Our work revealed that ATP1A2^L809R variants cause a predisposition to epilepsy. Moreover, we provide a point variant mouse model for epilepsy research and drug screening.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

SIGNIFICANCE

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

A Phenotypic-Driven Approach for the Diagnosis of WOREE Syndrome

BACKGROUND

WOREE syndrome is a rare neurodevelopmental disorder featuring drug-resistant epilepsy and global developmental delay. The disease, caused by biallelic pathogenic variants in the WWOX gene, usually leads to severe disability or death within the first years of life. Clinicians have become more confident with the phenotypic picture of WOREE syndrome, allowing earlier clinical diagnosis. We report a boy with a peculiar clinic-radiological pattern supporting the diagnosis of WOREE syndrome.

METHODS

DNA was extracted from blood samples of the proband and his parents and subjected to Exome Sequencing (ES). Agarose gel electrophoresis, real-time quantitative PCR (Q-PCR), and array-CGH 180K were also performed.

RESULTS

ES detected a pathogenic stop variant (c.790C > T, p.Arg264*) in one allele of WWOX in the proband and his unaffected mother. A 180K array-CGH analysis revealed a 84,828-bp (g.chr16:78,360,803-78,445,630) deletion encompassing exon 6. The Q-PCR product showed that the proband and his father harbored the same deleted fragment, fusing exons 5 and 7 of WWOX.

CONCLUSIONS

Genetic testing remains crucial in establishing the definitive diagnosis of WOREE syndrome and allows prenatal interventions/parental counseling. However, our findings suggest that targeted Next Generation Sequencing-based testing may occasionally show technical pitfalls, prompting further genetic investigation in selected cases with high clinical suspicion.

Monogenic developmental and epileptic encephalopathies of infancy and childhood, a population cohort from Norway

INTRODUCTION

Developmental and epileptic encephalopathies (DEE) is a group of epilepsies where the epileptic activity, seizures and the underlying neurobiology contributes to cognitive and behavioral impairments. Uncovering the causes of DEE is important in order to develop guidelines for treatment and follow-up. The aim of the present study was to describe the clinical picture and to identify genetic causes in a patient cohort with DEE without known etiology, from a Norwegian regional hospital.

METHODS

Systematic searches of medical records were performed at Drammen Hospital, Vestre Viken Health Trust, to identify patients with epilepsy in the period 1999-2018. Medical records were reviewed to identify patients with DEE of unknown cause. In 2018, patients were also recruited consecutively from treating physicians. All patients underwent thorough clinical evaluation and updated genetic diagnostic analyses.

RESULTS

Fifty-five of 2,225 patients with epilepsy had DEE of unknown etiology. Disease-causing genetic variants were found in 15/33 (45%) included patients. Three had potentially treatable metabolic disorders (SLC2A1, COQ4 and SLC6A8). Developmental comorbidity was higher in the group with a genetic diagnosis, compared to those who remained undiagnosed. Five novel variants in known genes were found, and the patient phenotypes are described.

CONCLUSION

The results from this study illustrate the importance of performing updated genetic investigations and/or analyses in patients with DEE of unknown etiology. A genetic cause was identified in 45% of the patients, and three of these patients had potentially treatable conditions where available targeted therapy may improve patient outcome.

High-impact rare genetic variants in severe schizophrenia

Extreme phenotype sequencing has led to the identification of high-impact rare genetic variants for many complex disorders but has not been applied to studies of severe schizophrenia. We sequenced 112 individuals with severe, extremely treatment-resistant schizophrenia, 218 individuals with typical schizophrenia, and 4,929 controls. We compared the burden of rare, damaging missense and loss-of-function variants between severe, extremely treatment-resistant schizophrenia, typical schizophrenia, and controls across mutation intolerant genes. Individuals with severe, extremely treatment-resistant schizophrenia had a high burden of rare loss-of-function (odds ratio, 1.91; 95% CI, 1.39 to 2.63; P = 7.8 × 10-5) and damaging missense variants in intolerant genes (odds ratio, 2.90; 95% CI, 2.02 to 4.15; P = 3.2 × 10-9). A total of 48.2% of individuals with severe, extremely treatment-resistant schizophrenia carried at least one rare, damaging missense or loss-of-function variant in intolerant genes compared to 29.8% of typical schizophrenia individuals (odds ratio, 2.18; 95% CI, 1.33 to 3.60; P = 1.6 × 10-3) and 25.4% of controls (odds ratio, 2.74; 95% CI, 1.85 to 4.06; P = 2.9 × 10-7). Restricting to genes previously associated with schizophrenia risk strengthened the enrichment with 8.9% of individuals with severe, extremely treatment-resistant schizophrenia carrying a damaging missense or loss-of-function variant compared to 2.3% of typical schizophrenia (odds ratio, 5.48; 95% CI, 1.52 to 19.74; P = 0.02) and 1.6% of controls (odds ratio, 5.82; 95% CI, 3.00 to 11.28; P = 2.6 × 10-8). These results demonstrate the power of extreme phenotype case selection in psychiatric genetics and an approach to augment schizophrenia gene discovery efforts.

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

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

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

Background

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

Methods

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

Findings

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

Interpretation

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

Funding

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