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

in EBioMedicine (2021), 72

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 ... [more ▼]

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. [less ▲]

Machine learning-based identification and characterization of 15 novel pathogenic SUOX missense mutations

in Molecular Genetics and Metabolism (2021)

Isolated sulfite oxidase deficiency (ISOD) is a rare hereditary metabolic disease caused by absence of functional sulfite oxidase (SO) due to mutations of the SUOX gene. SO oxidizes toxic sulfite and ... [more ▼]

Isolated sulfite oxidase deficiency (ISOD) is a rare hereditary metabolic disease caused by absence of functional sulfite oxidase (SO) due to mutations of the SUOX gene. SO oxidizes toxic sulfite and sulfite accumulation is associated with neurological disorders, progressive brain atrophy and early death. Similarities of these neurological symptoms to abundant diseases like neonatal encephalopathy underlines the raising need to increase the awareness for ISOD. Here we report an interdisciplinary approach utilizing exome/genome data derived from gnomAD database as well as published variants to predict the pathogenic outcome of 303 naturally occurring SO missense variants and combining these with activity determination. We identified 15 novel ISOD-causing SO variants and generated a databank of pathogenic SO missense variants to support future diagnosis of ISOD patients. We found six inactive variants (W101G, H118Y, E197K, R217W, S427W, D512Y, Q518R) and seven (D110H, P119S, G121E, G130R, Y140C, R269H, Q396P, R459Q) with severe reduction in activity. Based on the Hardy-Weinberg-equilibrium and the combination of our results with published SO missense and protein truncating variants, we calculated the first comprehensive incidence rate for ISOD of 1 in 1,377,341 births and provide a pathogenicity score to 303 naturally occurring SO missense variants. [less ▲]

Assessing the role of polygenic background on the penetrance of monogenic forms in Parkinson\textquoterights disease. 2021.06.06.21253270

E-print/Working paper (2021)

Background: Several rare and common variants are associated with Parkinson's disease. However, there is still an incomplete penetrance in the carriers of rare variants associated with Parkinson's disease ... [more ▼]

Background: Several rare and common variants are associated with Parkinson's disease. However, there is still an incomplete penetrance in the carriers of rare variants associated with Parkinson's disease. To address this issue, we investigated whether a PRS calculated from significant GWAS SNPs affects the penetrance of Parkinson's disease among carriers of rare monogenic variants in known Parkinson's disease genes and those with a family history. Methods: We calculated the PRS based on common variants and selected the carriers of rare monogenic variants by using the exome data from UK Biobank. Individuals were divided into three risk categories based on PRS: low (<10%), intermediate (10%-90%), and high (>90%) risk groups. We then compared how PRS affects Parkinson\textquoterights disease risk among carriers of rare monogenic variants and those with family-history. Results: We observed a two-fold higher odds ratio for a carrier of a monogenic variant that had a high PRS (OR 4.07,95\% CI, 1.72-8.08) compared to carriers with a low PRS (OR 1.91, 95\% CI, 0.31-6.05). In the same line, carriers with a first-degree family history and with \>90\% PRS have even a higher risk of developing PD (OR 23.53, 95\%CI 5.39-71.54) compared to those with \<90\% PRS (OR 9.54, 95\% CI 3.32-21.65). Conclusions: Our results show that PRS, carrier status, and family history contribute independently and additively to the Parkinson's disease risk. [less ▲]

Mantis: flexible and consensus-driven genome annotation

in GigaScience (2021), 10(6),

The rapid development of the (meta-)omics fields has produced an unprecedented amount of high-resolution and high-fidelity data. Through the use of these datasets we can infer the role of previously ... [more ▼]

The rapid development of the (meta-)omics fields has produced an unprecedented amount of high-resolution and high-fidelity data. Through the use of these datasets we can infer the role of previously functionally unannotated proteins from single organisms and consortia. In this context, protein function annotation can be described as the identification of regions of interest (i.e., domains) in protein sequences and the assignment of biological functions. Despite the existence of numerous tools, challenges remain in terms of speed, flexibility, and reproducibility. In the big data era, it is also increasingly important to cease limiting our findings to a single reference, coalescing knowledge from different data sources, and thus overcoming some limitations in overly relying on computationally generated data from single sources.We implemented a protein annotation tool, Mantis, which uses database identifiers intersection and text mining to integrate knowledge from multiple reference data sources into a single consensus-driven output. Mantis is flexible, allowing for the customization of reference data and execution parameters, and is reproducible across different research goals and user environments. We implemented a depth-first search algorithm for domain-specific annotation, which significantly improved annotation performance compared to sequence-wide annotation. The parallelized implementation of Mantis results in short runtimes while also outputting high coverage and high-quality protein function annotations.Mantis is a protein function annotation tool that produces high-quality consensus-driven protein annotations. It is easy to set up, customize, and use, scaling from single genomes to large metagenomes. Mantis is available under the MIT license at https://github.com/PedroMTQ/mantis. [less ▲]

Shared genetic basis between genetic generalized epilepsy and background electroencephalographic oscillations

in Epilepsia (2021), 62(7), 1518-1527

Abstract Objective Paroxysmal epileptiform abnormalities on electroencephalography (EEG) are the hallmark of epilepsies, but it is uncertain to what extent epilepsy and background EEG oscillations share ... [more ▼]

Abstract Objective Paroxysmal epileptiform abnormalities on electroencephalography (EEG) are the hallmark of epilepsies, but it is uncertain to what extent epilepsy and background EEG oscillations share neurobiological underpinnings. Here, we aimed to assess the genetic correlation between epilepsy and background EEG oscillations. Methods Confounding factors, including the heterogeneous etiology of epilepsies and medication effects, hamper studies on background brain activity in people with epilepsy. To overcome this limitation, we compared genetic data from a genome-wide association study (GWAS) on epilepsy (n = 12 803 people with epilepsy and 24 218 controls) with that from a GWAS on background EEG (n = 8425 subjects without epilepsy), in which background EEG oscillation power was quantified in four different frequency bands: alpha, beta, delta, and theta. We replicated our findings in an independent epilepsy replication dataset (n = 4851 people with epilepsy and 20 428 controls). To assess the genetic overlap between these phenotypes, we performed genetic correlation analyses using linkage disequilibrium score regression, polygenic risk scores, and Mendelian randomization analyses. Results Our analyses show strong genetic correlations of genetic generalized epilepsy (GGE) with background EEG oscillations, primarily in the beta frequency band. Furthermore, we show that subjects with higher beta and theta polygenic risk scores have a significantly higher risk of having generalized epilepsy. Mendelian randomization analyses suggest a causal effect of GGE genetic liability on beta oscillations. Significance Our results point to shared biological mechanisms underlying background EEG oscillations and the susceptibility for GGE, opening avenues to investigate the clinical utility of background EEG oscillations in the diagnostic workup of epilepsy. [less ▲]

Evaluation of the Molecular Pathogenesis of Adrenocortical Tumors by Whole-Genome Sequencing

in Journal of the Endocrine Society (2021, May 03), 5(Issue Supplement_1), 68

Pathogenesis of autonomous steroid secretion and adrenocortical tumorigenesis remains partially obscure. Our aim was to identify novel genetic alterations in adrenocortical adenomas (ACA) without somatic ... [more ▼]

Pathogenesis of autonomous steroid secretion and adrenocortical tumorigenesis remains partially obscure. Our aim was to identify novel genetic alterations in adrenocortical adenomas (ACA) without somatic mutations in known driver genes. Whole-genome sequencing was performed on 26 ACA/blood-derived DNA pairs without driver mutations in PRKACA, GNAS and CTNNB1 genes at previous WES (ENSAT study JCEM 2016). These included 12 cortisol-producing adenomas with Cushing syndrome (CS-CPAs), 7 with mild autonomous cortisol secretion (MACS-CPAs), and 7 endocrine-inactive ACAs (EIAs). Seven adrenocortical carcinomas (ACC) were added to the cohort. We developed a bioinformatics pipeline for a comprehensive genome analysis and to reveal differences in variant distribution. Strelka, VarScan2 and ANNOVAR software and an in-house confidence score were used for variant calling and functional annotation. Combined Annotation-Dependent-Depletion (CADD) values were used to prioritize pathogenic variants. Additional focus relied on variants in pathogenically known pathways (Wnt/β-catenin, cAMP/PKA pathway). NovoBreak algorithm was applied to discover structural variations. Two hypermutated CS-CPA samples were excluded from further analysis. Using different filters, we detected variants in driver genes not observed at WES (one p.S45P in CTNNB1 and one p.R206L in PRKACA in two different CS-CPAs). In total, we report 179,830 variations (179,598 SNVs; 232 indels) throughout all samples, being more abundant in ACC (88,954) compared to ACA (CS-CPAs: 31,821; MACS-CPAs: 35,008; EIAs: 29,963). Most alterations were in intergenic (>50%), followed by intronic and ncRNA intronic regions. A total of 32 predicted pathogenic variants were found in both coding (CADD values ≥ 15) and non-coding (CADD values ≥ 5) regions. We found 3,301 possibly damaging and recurrent variants (intergenic mutations removed) (CS-CPAs: 1,463; MACS-CPAs: 1,549; EIAs: 1,268; ACC: 1,660), mostly accumulated in intronic regions. Some of these were detected in members of the Wnt/β-catenin (CS-CPAs: 6; MACS-CPAs: 2; EIA: 1) and cAMP/PKA (CS-CPAs: 6; MACS-CPAs: 7; EIA: 4) pathways (e.g. ADCY1, ADCY2, GNA13, PDE11A). We also found a slightly higher number of structural variations in EIA (3,620) and ACC (3,486) compared to CS-CPAs (977) and MACS-CPAs (2,119). In conclusion, still unrevealed genetic alterations, especially in intronic regions, may accompany early adrenal tumorigenesis and/or autonomous cortisol secretion. [less ▲]

Loss of function variants in the KCNQ5 gene are associated with genetic generalized epilepsies

E-print/Working paper (2021)

Objective: De novo missense variants in KCNQ5, encoding the voltage gated K+ channel KV7.5, have been described as a cause of developmental and epileptic encephalopathy (DEE) or intellectual disability ... [more ▼]

Objective: De novo missense variants in KCNQ5, encoding the voltage gated K+ channel KV7.5, have been described as a cause of developmental and epileptic encephalopathy (DEE) or intellectual disability (ID). We set out to identify disease-related KCNQ5 variants in genetic generalized epilepsy (GGE) and their underlying mechanisms. Methods: 1292 families with GGE were studied by next-generation sequencing. Whole-cell patch-clamp recordings, biotinylation and phospholipid overlay assays were performed in mammalian cells combined with docking and homology modeling. Results: We identified three deleterious heterozygous missense variants, one truncation and one splice site alteration in five independent families with GGE with predominant absence seizures, two variants were also associated with mild to moderate ID. All three missense variants displayed a strongly decreased current density indicating a loss-of-function (LOF). When mutant channels were co-expressed with wild-type (WT) KV7.5 or KV7.5 and KV7.3 channels, three variants also revealed a significant dominant-negative effect on WT channels. Other gating parameters were unchanged. Biotinylation assays indicated a normal surface expression of the variants. The p.Arg359Cys variant altered PI(4,5)P2-interaction, presumably in the non-conducting preopen-closed state. Interpretation: Our study indicates that specific deleterious KCNQ5 variants are associated with GGE, partially combined with mild to moderate ID. The disease mechanism is a LOF partially with dominant-negative effects through functional, rather than trafficking deficits. LOF of KV7.5 channels will reduce the M-current, likely resulting in increased excitability of KV7.5-expressing neurons. Further studies on a network level are necessary to understand which circuits are affected and how the variants induce generalized seizures. [less ▲]

Burden of rare variants in synaptic genes in patients with severe tinnitus: An exome based extreme phenotype study

in EBioMedicine (2021), 66(103309),

Background: tinnitus is a heterogeneous condition associated with audio logical and/or mental disorders. Chronic, severe tinnitus is reported in 1% of the population and it shows a relevant heritability ... [more ▼]

Background: tinnitus is a heterogeneous condition associated with audio logical and/or mental disorders. Chronic, severe tinnitus is reported in 1% of the population and it shows a relevant heritability, according to twins, adoptees and familial aggregation studies. The genetic contribution to severe tinnitus is unknown since large genomic studies include individuals with self-reported tinnitus and large heterogeneity in the phenotype. The aim of this study was to identify genes for severe tinnitus in patients with extreme phenotype. Methods: for this extreme phenotype study, we used three different cohorts with European ancestry (Spanish with Meniere disease (MD), Swedish tinnitus and European generalized epilepsy). In addition, four independent control datasets were also used for comparisons. Whole-exome sequencing was performed for the MD and epilepsy cohorts and whole-genome sequencing was carried out in Swedish with tinnitus. Findings: we found an enrichment of rare missense variants in 24 synaptic genes in a Spanish cohort, the most significant being PRUNE2, AKAP9, SORBS1, ITGAX, ANK2, KIF20B and TSC2 (p < 2E−04), when they were compared with reference datasets. This burden was replicated for ANK2 gene in a Swedish cohort with 97 tinnitus individuals, and in a subset of 34 Swedish patients with severe tinnitus for ANK2, AKAP9 and TSC2 genes (p < 2E−02). However, these associations were not significant in a third cohort of 701 generalized epilepsy individuals without tinnitus. Gene ontology (GO) and gene-set enrichment analyses revealed several pathways and biological processes involved in severe tinnitus, including membrane trafficking and cytoskeletal protein binding in neurons. Interpretation: a burden of rare variants in ANK2, AKAP9 and TSC2 is associated with severe tinnitus. ANK2, encodes a cytoskeleton scaffolding protein that coordinates the assembly of several proteins, drives axonal branching and influences connectivity in neurons. [less ▲]

Genotype-phenotype correlations in SCN8A-related disorders reveal prognostic and therapeutic implications

E-print/Working paper (2021)

We report detailed functional analyses and genotype-phenotype correlations in 433 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6. Five different ... [more ▼]

We report detailed functional analyses and genotype-phenotype correlations in 433 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6. Five different clinical subgroups could be identified: 1) Benign familial infantile epilepsy (BFIE) (n=17, normal cognition, treatable seizures), 2) intermediate epilepsy (n=36, mild ID, partially pharmacoresponsive), 3) developmental and epileptic encephalopathy (DEE, n=191, severe ID, majority pharmacoresistant), 4) generalized epilepsy (n=21, mild to moderate ID, frequently with absence seizures), and 5) affected individuals without epilepsy (n=25, mild to moderate ID). Groups 1-3 presented with early-onset (median: four months) focal or multifocal seizures and epileptic discharges, whereas the onset of seizures in group 4 was later (median: 39 months) with generalized epileptic discharges. The epilepsy was not classifiable in 143 individuals. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin insensitive human NaV1.6 channels and whole-cell patch clamping. Two variants causing DEE showed a strong gain-of-function (GOF, hyperpolarising shift of steady-state activation, strongly increased neuronal firing rate), and one variant causing BFIE or intermediate epilepsy showed a mild GOF (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (LOF, reduced current amplitudes, depolarising shift of steady-state activation, reduced neuronal firing). Including previous studies, functional effects were known for 165 individuals. All 133 individuals carrying GOF variants had either focal (76, groups 1-3), or unclassifiable epilepsy (37), whereas 32 with LOF variants had either generalized (14), no (11) or unclassifiable (5) epilepsy; only two had DEE. Computational modeling in the GOF group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. GOF variant carriers responded significantly better to sodium channel blockers (SCBs) than to other anti-seizure medications, and the same applied for all individuals of groups 1-3.In conclusion, our data reveal clear genotype-phenotype correlations between age at seizure onset, type of epilepsy and gain- or loss-of-function effects of SCN8A variants. Generalized epilepsy with absence seizures is the main epilepsy phenotype of LOF variant carriers and the extent of the electrophysiological dysfunction of the GOF variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that SCBs present a therapeutic treatment option in early onset SCN8A-related focal epilepsy. [less ▲]

Statistical and Machine Learning Techniques in Human Microbiome Studies: Contemporary Challenges and Solutions

in Frontiers in Microbiology (2021), 12

The human microbiome has emerged as a central research topic in human biology and biomedicine. Current microbiome studies generate high-throughput omics data across different body sites, populations, and ... [more ▼]

The human microbiome has emerged as a central research topic in human biology and biomedicine. Current microbiome studies generate high-throughput omics data across different body sites, populations, and life stages. Many of the challenges in microbiome research are similar to other high-throughput studies, the quantitative analyses need to address the heterogeneity of data, specific statistical properties, and the remarkable variation in microbiome composition across individuals and body sites. This has led to a broad spectrum of statistical and machine learning challenges that range from study design, data processing, and standardization to analysis, modeling, cross-study comparison, prediction, data science ecosystems, and reproducible reporting. Nevertheless, although many statistics and machine learning approaches and tools have been developed, new techniques are needed to deal with emerging applications and the vast heterogeneity of microbiome data. We review and discuss emerging applications of statistical and machine learning techniques in human microbiome studies and introduce the COST Action CA18131 “ML4Microbiome” that brings together microbiome researchers and machine learning experts to address current challenges such as standardization of analysis pipelines for reproducibility of data analysis results, benchmarking, improvement, or development of existing and new tools and ontologies. [less ▲]