FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Evaluation of polygenic scoring methods in five biobanks shows larger variation between biobanks than methods and finds benefits of ensemble learning.

Methods of estimating polygenic scores (PGSs) from genome-wide association studies are increasingly utilized. However, independent method evaluation is lacking, and method comparisons are often limited. Here, we evaluate polygenic scores derived via seven methods in five biobank studies (totaling about 1.2 million participants) across 16 diseases and quantitative traits, building on a reference-standardized framework. We conducted meta-analyses to quantify the effects of method choice, hyperparameter tuning, method ensembling, and the target biobank on PGS performance. We found that no single method consistently outperformed all others. PGS effect sizes were more variable between biobanks than between methods within biobanks when methods were well tuned. Differences between methods were largest for the two investigated autoimmune diseases, seropositive rheumatoid arthritis and type 1 diabetes. For most methods, cross-validation was more reliable for tuning hyperparameters than automatic tuning (without the use of target data). For a given target phenotype, elastic net models combining PGS across methods (ensemble PGS) tuned in the UK Biobank provided consistent, high, and cross-biobank transferable performance, increasing PGS effect sizes (ß coefficients) by a median of 5.0% relative to LDpred2 and MegaPRS (the two best-performing single methods when tuned with cross-validation). Our interactively browsable online-results and open-source workflow prspipe provide a rich resource and reference for the analysis of polygenic scoring methods across biobanks.

Epilepsies of presumed genetic etiology show enrichment of rare variants that occur in the general population.

Previous studies suggested that severe epilepsies, e.g., developmental and epileptic encephalopathies (DEEs), are mainly caused by ultra-rare de novo genetic variants. For milder disease, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 control individuals. Here, we separately analyzed three different groups of epilepsies: severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in control individuals with an allele count ≥ 1 and a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD ≥ 20), and to have an odds ratio in individuals with epilepsy ≥ 2. We identified genes enriched with QRVs primarily in NAFE (n = 72), followed by GGE (n = 32) and DEE (n = 21). This suggests that rare variants may play a more important role for causality of NAFE than for DEE. Moreover, we found that genes harboring QRVs, e.g., HSGP2, FLNA, or TNC, encode proteins that are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE that occur also in the general population, while in DEE and GGE, the contribution of such variants appears more limited.

Monoallelic intragenic POU3F2 variants lead to neurodevelopmental delay and hyperphagic obesity, confirming the gene's candidacy in 6q16.1 deletions.

While common obesity accounts for an increasing global health burden, its monogenic forms have taught us underlying mechanisms via more than 20 single-gene disorders. Among these, the most common mechanism is central nervous system dysregulation of food intake and satiety, often accompanied by neurodevelopmental delay (NDD) and autism spectrum disorder. In a family with syndromic obesity, we identified a monoallelic truncating variant in POU3F2 (alias BRN2) encoding a neural transcription factor, which has previously been suggested as a driver of obesity and NDD in individuals with the 6q16.1 deletion. In an international collaboration, we identified ultra-rare truncating and missense variants in another ten individuals sharing autism spectrum disorder, NDD, and adolescent-onset obesity. Affected individuals presented with low-to-normal birth weight and infantile feeding difficulties but developed insulin resistance and hyperphagia during childhood. Except for a variant leading to early truncation of the protein, identified variants showed adequate nuclear translocation but overall disturbed DNA-binding ability and promotor activation. In a cohort with common non-syndromic obesity, we independently observed a negative correlation of POU3F2 gene expression with BMI, suggesting a role beyond monogenic obesity. In summary, we propose deleterious intragenic variants of POU3F2 to cause transcriptional dysregulation associated with hyperphagic obesity of adolescent onset with variable NDD.

Identification of pathogenic variant enriched regions across genes and gene families.

Missense variant interpretation is challenging. Essential regions for protein function are conserved among gene-family members, and genetic variants within these regions are potentially more likely to confer risk to disease. Here, we generated 2871 gene-family protein sequence alignments involving 9990 genes and performed missense variant burden analyses to identify novel essential protein regions. We mapped 2,219,811 variants from the general population into these alignments and compared their distribution with 76,153 missense variants from patients. With this gene-family approach, we identified 465 regions enriched for patient variants spanning 41,463 amino acids in 1252 genes. As a comparison, by testing the same genes individually, we identified fewer patient variant enriched regions, involving only 2639 amino acids and 215 genes. Next, we selected de novo variants from 6753 patients with neurodevelopmental disorders and 1911 unaffected siblings and observed an 8.33-fold enrichment of patient variants in our identified regions (95% C.I. = 3.90-Inf, P-value = 2.72 × 10-11). By using the complete ClinVar variant set, we found that missense variants inside the identified regions are 106-fold more likely to be classified as pathogenic in comparison to benign classification (OR = 106.15, 95% C.I = 70.66-Inf, P-value < 2.2 × 10-16). All pathogenic variant enriched regions (PERs) identified are available online through "PER viewer," a user-friendly online platform for interactive data mining, visualization, and download. In summary, our gene-family burden analysis approach identified novel PERs in protein sequences. This annotation can empower variant interpretation.

Gene variant effects across sodium channelopathies predict function and guide precision therapy.

Pathogenic variants in the voltage-gated sodium channel gene family lead to early onset epilepsies, neurodevelopmental disorders, skeletal muscle channelopathies, peripheral neuropathies and cardiac arrhythmias. Disease-associated variants have diverse functional effects ranging from complete loss-of-function to marked gain-of-function. Therapeutic strategy is likely to depend on functional effect. Experimental studies offer important insights into channel function but are resource intensive and only performed in a minority of cases. Given the evolutionarily conserved nature of the sodium channel genes, we investigated whether similarities in biophysical properties between different voltage-gated sodium channels can predict function and inform precision treatment across sodium channelopathies. We performed a systematic literature search identifying functionally assessed variants in any of the nine voltage-gated sodium channel genes until 28 April 2021. We included missense variants that had been electrophysiologically characterized in mammalian cells in whole-cell patch-clamp recordings. We performed an alignment of linear protein sequences of all sodium channel genes and correlated variants by their overall functional effect on biophysical properties. Of 951 identified records, 437 sodium channel-variants met our inclusion criteria and were reviewed for functional properties. Of these, 141 variants were epilepsy-associated (SCN1/2/3/8A), 79 had a neuromuscular phenotype (SCN4/9/10/11A), 149 were associated with a cardiac phenotype (SCN5/10A) and 68 (16%) were considered benign. We detected 38 missense variant pairs with an identical disease-associated variant in a different sodium channel gene. Thirty-five out of 38 of those pairs resulted in similar functional consequences, indicating up to 92% biophysical agreement between corresponding sodium channel variants (odds ratio = 11.3; 95% confidence interval = 2.8 to 66.9; P < 0.001). Pathogenic missense variants were clustered in specific functional domains, whereas population variants were significantly more frequent across non-conserved domains (odds ratio = 18.6; 95% confidence interval = 10.9-34.4; P < 0.001). Pore-loop regions were frequently associated with loss-of-function variants, whereas inactivation sites were associated with gain-of-function (odds ratio = 42.1, 95% confidence interval = 14.5-122.4; P < 0.001), whilst variants occurring in voltage-sensing regions comprised a range of gain- and loss-of-function effects. Our findings suggest that biophysical characterisation of variants in one SCN-gene can predict channel function across different SCN-genes where experimental data are not available. The collected data represent the first gain- versus loss-of-function topological map of SCN proteins indicating shared patterns of biophysical effects aiding variant analysis and guiding precision therapy. We integrated our findings into a free online webtool to facilitate functional sodium channel gene variant interpretation (http://SCN-viewer.broadinstitute.org).

Analysing an allelic series of rare missense variants of CACNA1I in a Swedish schizophrenia cohort.

CACNA1I is implicated in the susceptibility to schizophrenia by large-scale genetic association studies of single nucleotide polymorphisms. However, the channelopathy of CACNA1I in schizophrenia is unknown. CACNA1I encodes CaV3.3, a neuronal voltage-gated calcium channel that underlies a subtype of T-type current that is important for neuronal excitability in the thalamic reticular nucleus and other regions of the brain. Here, we present an extensive functional characterization of 57 naturally occurring rare and common missense variants of CACNA1I derived from a Swedish schizophrenia cohort of more than 10 000 individuals. Our analysis of this allelic series of coding CACNA1I variants revealed that reduced CaV3.3 channel current density was the dominant phenotype associated with rare CACNA1I coding alleles derived from control subjects, whereas rare CACNA1I alleles from schizophrenia patients encoded CaV3.3 channels with altered responses to voltages. CACNA1I variants associated with altered current density primarily impact the ionic channel pore and those associated with altered responses to voltage impact the voltage-sensing domain. CaV3.3 variants associated with altered voltage dependence of the CaV3.3 channel and those associated with peak current density deficits were significantly segregated across affected and unaffected groups (Fisher's exact test, P = 0.034). Our results, together with recent data from the SCHEMA (Schizophrenia Exome Sequencing Meta-Analysis) cohort, suggest that reduced CaV3.3 function may protect against schizophrenia risk in rare cases. We subsequently modelled the effect of the biophysical properties of CaV3.3 channel variants on thalamic reticular nucleus excitability and found that compared with common variants, ultrarare CaV3.3-coding variants derived from control subjects significantly decreased thalamic reticular nucleus excitability (P = 0.011). When all rare variants were analysed, there was a non-significant trend between variants that reduced thalamic reticular nucleus excitability and variants that either had no effect or increased thalamic reticular nucleus excitability across disease status. Taken together, the results of our functional analysis of an allelic series of >50 CACNA1I variants in a schizophrenia cohort reveal that loss of function of CaV3.3 is a molecular phenotype associated with reduced disease risk burden, and our approach may serve as a template strategy for channelopathies in polygenic disorders.

Comprehensive characterization of amino acid positions in protein structures reveals molecular effect of missense variants.

Interpretation of the colossal number of genetic variants identified from sequencing applications is one of the major bottlenecks in clinical genetics, with the inference of the effect of amino acid-substituting missense variations on protein structure and function being especially challenging. Here we characterize the three-dimensional (3D) amino acid positions affected in pathogenic and population variants from 1,330 disease-associated genes using over 14,000 experimentally solved human protein structures. By measuring the statistical burden of variations (i.e., point mutations) from all genes on 40 3D protein features, accounting for the structural, chemical, and functional context of the variations' positions, we identify features that are generally associated with pathogenic and population missense variants. We then perform the same amino acid-level analysis individually for 24 protein functional classes, which reveals unique characteristics of the positions of the altered amino acids: We observe up to 46% divergence of the class-specific features from the general characteristics obtained by the analysis on all genes, which is consistent with the structural diversity of essential regions across different protein classes. We demonstrate that the function-specific 3D features of the variants match the readouts of mutagenesis experiments for BRCA1 and PTEN, and positively correlate with an independent set of clinically interpreted pathogenic and benign missense variants. Finally, we make our results available through a web server to foster accessibility and downstream research. Our findings represent a crucial step toward translational genetics, from highlighting the impact of mutations on protein structure to rationalizing the variants' pathogenicity in terms of the perturbed molecular mechanisms.

KCNT1-related epilepsies and epileptic encephalopathies: phenotypic and mutational spectrum.

Variants in KCNT1, encoding a sodium-gated potassium channel (subfamily T member 1), have been associated with a spectrum of epilepsies and neurodevelopmental disorders. These range from familial autosomal dominant or sporadic sleep-related hypermotor epilepsy to epilepsy of infancy with migrating focal seizures (EIMFS) and include developmental and epileptic encephalopathies. This study aims to provide a comprehensive overview of the phenotypic and genotypic spectrum of KCNT1 mutation-related epileptic disorders in 248 individuals, including 66 previously unpublished and 182 published cases, the largest cohort reported so far. Four phenotypic groups emerged from our analysis: (i) EIMFS (152 individuals, 33 previously unpublished); (ii) developmental and epileptic encephalopathies other than EIMFS (non-EIMFS developmental and epileptic encephalopathies) (37 individuals, 17 unpublished); (iii) autosomal dominant or sporadic sleep-related hypermotor epilepsy (53 patients, 14 unpublished); and (iv) other phenotypes (six individuals, two unpublished). In our cohort of 66 new cases, the most common phenotypic features were: (i) in EIMFS, heterogeneity of seizure types, including epileptic spasms, epilepsy improvement over time, no epilepsy-related deaths; (ii) in non-EIMFS developmental and epileptic encephalopathies, possible onset with West syndrome, occurrence of atypical absences, possible evolution to developmental and epileptic encephalopathies with sleep-related hypermotor epilepsy features; one case of sudden unexplained death in epilepsy; (iii) in autosomal dominant or sporadic sleep-related hypermotor epilepsy, we observed a high prevalence of drug-resistance, although seizure frequency improved with age in some individuals, appearance of cognitive regression after seizure onset in all patients, no reported severe psychiatric disorders, although behavioural/psychiatric comorbidities were reported in ∼50% of the patients, sudden unexplained death in epilepsy in one individual; and (iv) other phenotypes in individuals with mutation of KCNT1 included temporal lobe epilepsy, and epilepsy with tonic-clonic seizures and cognitive regression. Genotypic analysis of the whole cohort of 248 individuals showed only missense mutations and one inframe deletion in KCNT1. Although the KCNT1 mutations in affected individuals were seen to be distributed among the different domains of the KCNT1 protein, genotype-phenotype considerations showed many of the autosomal dominant or sporadic sleep-related hypermotor epilepsy-associated mutations to be clustered around the RCK2 domain in the C terminus, distal to the NADP domain. Mutations associated with EIMFS/non-EIMFS developmental and epileptic encephalopathies did not show a particular pattern of distribution in the KCNT1 protein. Recurrent KCNT1 mutations were seen to be associated with both severe and less severe phenotypes. Our study further defines and broadens the phenotypic and genotypic spectrums of KCNT1-related epileptic conditions and emphasizes the increasingly important role of this gene in the pathogenesis of early onset developmental and epileptic encephalopathies as well as of focal epilepsies, namely autosomal dominant or sporadic sleep-related hypermotor epilepsy.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe.

Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.

A novel compound heterozygous leptin receptor mutation causes more severe obesity than in Leprdb/db mice.

The leptin receptor (Lepr) pathway is important for food intake regulation, energy expenditure, and body weight. Mutations in leptin and the Lepr have been shown to cause early-onset severe obesity in mice and humans. In studies with C57BL/6NCrl mice, we found a mouse with extreme obesity. To identify a putative spontaneous new form of monogenic obesity, we performed backcross studies with this mouse followed by a quantitative trait locus (QTL) analysis and sequencing of the selected chromosomal QTL region. We thereby identified a novel Lepr mutation (C57BL/6N-LeprL536Hfs*6-1NKB), which is located at chromosome 4, exon 11 within the CRH2-leptin-binding site. Compared with C57BL/6N mice, LeprL536Hfs*6 develop early onset obesity and their body weight exceeds that of Leprdb/db mice at an age of 30 weeks. Similar to Leprdb/db mice, the LeprL536Hfs*6 model is characterized by hyperphagia, obesity, lower energy expenditure and activity, hyperglycemia, and hyperinsulinemia compared with C57BL/6N mice. Crossing Leprdb/wt with LeprL536Hfs*6/wt mice results in compound heterozygous LeprL536Hfs*6/db mice, which develop even higher body weight and fat mass than both homozygous Leprdb/db and LeprL536Hfs*6 mice. Compound heterozygous Lepr deficiency affecting functionally different regions of the Lepr causes more severe obesity than the parental homozygous mutations.

Adipose tissue derived bacteria are associated with inflammation in obesity and type 2 diabetes.

OBJECTIVE: Bacterial translocation to various organs including human adipose tissue (AT) due to increased intestinal permeability remains poorly understood. We hypothesised that: (1) bacterial presence is highly tissue specific and (2) related in composition and quantity to immune inflammatory and metabolic burden. DESIGN: We quantified and sequenced the bacterial 16S rRNA gene in blood and AT samples (omental, mesenteric and subcutaneous) of 75 subjects with obesity with or without type 2 diabetes (T2D) and used catalysed reporter deposition (CARD) - fluorescence in situ hybridisation (FISH) to detect bacteria in AT. RESULTS: Under stringent experimental and bioinformatic control for contaminants, bacterial DNA was detected in blood and omental, subcutaneous and mesenteric AT samples in the range of 0.1 to 5 pg/µg DNA isolate. Moreover, CARD-FISH allowed the detection of living, AT-borne bacteria. Proteobacteria and Firmicutes were the predominant phyla, and bacterial quantity was associated with immune cell infiltration, inflammatory and metabolic parameters in a tissue-specific manner. Bacterial composition differed between subjects with and without T2D and was associated with related clinical measures, including systemic and tissues-specific inflammatory markers. Finally, treatment of adipocytes with bacterial DNA in vitro stimulated the expression of TNFA and IL6. CONCLUSIONS: Our study provides contaminant aware evidence for the presence of bacteria and bacterial DNA in several ATs in obesity and T2D and suggests an important role of bacteria in initiating and sustaining local AT subclinical inflammation and therefore impacting metabolic sequelae of obesity.

Genome-wide meta-analysis identifies novel determinants of circulating serum progranulin.

Progranulin is a secreted protein with important functions in processes including immune and inflammatory response, metabolism and embryonic development. The present study aimed at identification of genetic factors determining progranulin concentrations. We conducted a genome-wide association meta-analysis for serum progranulin in three independent cohorts from Europe: Sorbs (N = 848) and KORA (N = 1628) from Germany and PPP-Botnia (N = 335) from Finland (total N = 2811). Single nucleotide polymorphisms (SNPs) associated with progranulin levels were replicated in two additional German cohorts: LIFE-Heart Study (Leipzig; N = 967) and Metabolic Syndrome Berlin Potsdam (Berlin cohort; N = 833). We measured mRNA expression of genes in peripheral blood mononuclear cells (PBMC) by micro-arrays and performed mRNA expression quantitative trait and expression-progranulin association studies to functionally substantiate identified loci. Finally, we conducted siRNA silencing experiments in vitro to validate potential candidate genes within the associated loci. Heritability of circulating progranulin levels was estimated at 31.8% and 26.1% in the Sorbs and LIFE-Heart cohort, respectively. SNPs at three loci reached study-wide significance (rs660240 in CELSR2-PSRC1-MYBPHL-SORT1, rs4747197 in CDH23-PSAP and rs5848 in GRN) explaining 19.4%/15.0% of the variance and 61%/57% of total heritability in the Sorbs/LIFE-Heart Study. The strongest evidence for association was at rs660240 (P = 5.75 × 10-50), which was also associated with mRNA expression of PSRC1 in PBMC (P = 1.51 × 10-21). Psrc1 knockdown in murine preadipocytes led to a consecutive 30% reduction in progranulin secretion. In conclusion, the present meta-GWAS combined with mRNA expression identified three loci associated with progranulin and supports the role of PSRC1 in the regulation of progranulin secretion.

Biological concepts in human sodium channel epilepsies and their relevance in clinical practice.

OBJECTIVE: Voltage-gated sodium channels (SCNs) share similar amino acid sequence, structure, and function. Genetic variants in the four human brain-expressed SCN genes SCN1A/2A/3A/8A have been associated with heterogeneous epilepsy phenotypes and neurodevelopmental disorders. To better understand the biology of seizure susceptibility in SCN-related epilepsies, our aim was to determine similarities and differences between sodium channel disorders, allowing us to develop a broader perspective on precision treatment than on an individual gene level alone. METHODS: We analyzed genotype-phenotype correlations in large SCN-patient cohorts and applied variant constraint analysis to identify severe sodium channel disease. We examined temporal patterns of human SCN expression and correlated functional data from in vitro studies with clinical phenotypes across different sodium channel disorders. RESULTS: Comparing 865 epilepsy patients (504 SCN1A, 140 SCN2A, 171 SCN8A, four SCN3A, 46 copy number variation [CNV] cases) and analysis of 114 functional studies allowed us to identify common patterns of presentation. All four epilepsy-associated SCN genes demonstrated significant constraint in both protein truncating and missense variation when compared to other SCN genes. We observed that age at seizure onset is related to SCN gene expression over time. Individuals with gain-of-function SCN2A/3A/8A missense variants or CNV duplications share similar characteristics, most frequently present with early onset epilepsy (<3 months), and demonstrate good response to sodium channel blockers (SCBs). Direct comparison of corresponding SCN variants across different SCN subtypes illustrates that the functional effects of variants in corresponding channel locations are similar; however, their clinical manifestation differs, depending on their role in different types of neurons in which they are expressed. SIGNIFICANCE: Variant function and location within one channel can serve as a surrogate for variant effects across related sodium channels. Taking a broader view on precision treatment suggests that in those patients with a suspected underlying genetic epilepsy presenting with neonatal or early onset seizures (<3 months), SCBs should be considered.

GRIN2A-related disorders: genotype and functional consequence predict phenotype.

Alterations of the N-methyl-d-aspartate receptor (NMDAR) subunit GluN2A, encoded by GRIN2A, have been associated with a spectrum of neurodevelopmental disorders with prominent speech-related features, and epilepsy. We performed a comprehensive assessment of phenotypes with a standardized questionnaire in 92 previously unreported individuals with GRIN2A-related disorders. Applying the criteria of the American College of Medical Genetics and Genomics to all published variants yielded 156 additional cases with pathogenic or likely pathogenic variants in GRIN2A, resulting in a total of 248 individuals. The phenotypic spectrum ranged from normal or near-normal development with mild epilepsy and speech delay/apraxia to severe developmental and epileptic encephalopathy, often within the epilepsy-aphasia spectrum. We found that pathogenic missense variants in transmembrane and linker domains (misTMD+Linker) were associated with severe developmental phenotypes, whereas missense variants within amino terminal or ligand-binding domains (misATD+LBD) and null variants led to less severe developmental phenotypes, which we confirmed in a discovery (P = 10-6) as well as validation cohort (P = 0.0003). Other phenotypes such as MRI abnormalities and epilepsy types were also significantly different between the two groups. Notably, this was paralleled by electrophysiology data, where misTMD+Linker predominantly led to NMDAR gain-of-function, while misATD+LBD exclusively caused NMDAR loss-of-function. With respect to null variants, we show that Grin2a+/- cortical rat neurons also had reduced NMDAR function and there was no evidence of previously postulated compensatory overexpression of GluN2B. We demonstrate that null variants and misATD+LBD of GRIN2A do not only share the same clinical spectrum (i.e. milder phenotypes), but also result in similar electrophysiological consequences (loss-of-function) opposing those of misTMD+Linker (severe phenotypes; predominantly gain-of-function). This new pathomechanistic model may ultimately help in predicting phenotype severity as well as eligibility for potential precision medicine approaches in GRIN2A-related disorders.

Targeted gene sequencing in 6994 individuals with neurodevelopmental disorder with epilepsy.

PURPOSE: We aimed to gain insight into frequencies of genetic variants in genes implicated in neurodevelopmental disorder with epilepsy (NDD+E) by investigating large cohorts of patients in a diagnostic setting. METHODS: We analyzed variants in NDD+E using epilepsy gene panel sequencing performed between 2013 and 2017 by two large diagnostic companies. We compared variant frequencies in 6994 panels with another 8588 recently published panels as well as exome-wide de novo variants in 1942 individuals with NDD+E and 10,937 controls. RESULTS: Genes with highest frequencies of ultrarare variants in NDD+E comprised SCN1A, KCNQ2, SCN2A, CDKL5, SCN8A, and STXBP1, concordant with the two other epilepsy cohorts we investigated. In only 46% of the analyzed 262 dominant and X-linked panel genes ultrarare variants in patients were reported. Among genes with contradictory evidence of association with epilepsy, CACNB4, CLCN2, EFHC1, GABRD, MAGI2, and SRPX2 showed equal frequencies in cases and controls. CONCLUSION: We show that improvement of panel design increased diagnostic yield over time, but panels still display genes with low or no diagnostic yield. With our data, we hope to improve current diagnostic NDD+E panel design and provide a resource of ultrarare variants in individuals with NDD+E to the community.

Delineating SPTAN1 associated phenotypes: from isolated epilepsy to encephalopathy with progressive brain atrophy.

De novo in-frame deletions and duplications in the SPTAN1 gene, encoding the non-erythrocyte αII spectrin, have been associated with severe West syndrome with hypomyelination and pontocerebellar atrophy. We aimed at comprehensively delineating the phenotypic spectrum associated with SPTAN1 mutations. Using different molecular genetic techniques, we identified 20 patients with a pathogenic or likely pathogenic SPTAN1 variant and reviewed their clinical, genetic and imaging data. SPTAN1 de novo alterations included seven unique missense variants and nine in-frame deletions/duplications of which 12 were novel. The recurrent three-amino acid duplication p.(Asp2303_Leu2305dup) occurred in five patients. Our patient cohort exhibited a broad spectrum of neurodevelopmental phenotypes, comprising six patients with mild to moderate intellectual disability, with or without epilepsy and behavioural disorders, and 14 patients with infantile epileptic encephalopathy, of which 13 had severe neurodevelopmental impairment and four died in early childhood. Imaging studies suggested that the severity of neurological impairment and epilepsy correlates with that of structural abnormalities as well as the mutation type and location. Out of seven patients harbouring mutations outside the α/ß spectrin heterodimerization domain, four had normal brain imaging and three exhibited moderately progressive brain and/or cerebellar atrophy. Twelve of 13 patients with mutations located within the spectrin heterodimer contact site exhibited severe and progressive brain, brainstem and cerebellar atrophy, with hypomyelination in most. We used fibroblasts from five patients to study spectrin aggregate formation by Triton-X extraction and immunocytochemistry followed by fluorescence microscopy. αII/ßII aggregates and αII spectrin in the insoluble protein fraction were observed in fibroblasts derived from patients with the mutations p.(Glu2207del), p.(Asp2303_Leu2305dup) and p.(Arg2308_Met2309dup), all falling in the nucleation site of the α/ß spectrin heterodimer region. Molecular modelling of the seven SPTAN1 amino acid changes provided preliminary evidence for structural alterations of the A-, B- and/or C-helices within each of the mutated spectrin repeats. We conclude that SPTAN1-related disorders comprise a wide spectrum of neurodevelopmental phenotypes ranging from mild to severe and progressive. Spectrin aggregate formation in fibroblasts with mutations in the α/ß heterodimerization domain seems to be associated with a severe neurodegenerative course and suggests that the amino acid stretch from Asp2303 to Met2309 in the α20 repeat is important for α/ß spectrin heterodimer formation and/or αII spectrin function.

GRIN2B encephalopathy: novel findings on phenotype, variant clustering, functional consequences and treatment aspects.

BACKGROUND: We aimed for a comprehensive delineation of genetic, functional and phenotypic aspects of GRIN2B encephalopathy and explored potential prospects of personalised medicine. METHODS: Data of 48 individuals with de novo GRIN2B variants were collected from several diagnostic and research cohorts, as well as from 43 patients from the literature. Functional consequences and response to memantine treatment were investigated in vitro and eventually translated into patient care. RESULTS: Overall, de novo variants in 86 patients were classified as pathogenic/likely pathogenic. Patients presented with neurodevelopmental disorders and a spectrum of hypotonia, movement disorder, cortical visual impairment, cerebral volume loss and epilepsy. Six patients presented with a consistent malformation of cortical development (MCD) intermediate between tubulinopathies and polymicrogyria. Missense variants cluster in transmembrane segments and ligand-binding sites. Functional consequences of variants were diverse, revealing various potential gain-of-function and loss-of-function mechanisms and a retained sensitivity to the use-dependent blocker memantine. However, an objectifiable beneficial treatment response in the respective patients still remains to be demonstrated. CONCLUSIONS: In addition to previously known features of intellectual disability, epilepsy and autism, we found evidence that GRIN2B encephalopathy is also frequently associated with movement disorder, cortical visual impairment and MCD revealing novel phenotypic consequences of channelopathies.