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.

Variants in the SCN1A gene are associated with a wide range of disorders including genetic epilepsy with febrile seizures plus (GEFS+), familial hemiplegic migraine (FHM), and the severe childhood epilepsy Dravet syndrome (DS). Predicting disease outcomes based on variant type remains challenging. Despite thousands of SCN1A variants being reported, only a minority has been functionally assessed.We review the functional SCN1A work performed to date, critically appraise electrophysiological measurements, compare this to in silico predictions, and relate our findings to the clinical phenotype.Our results show, regardless of the underlying phenotype, that conventional in silico software correctly predicted benign from pathogenic variants in nearly 90%, however was unable to differentiate within the disease spectrum (DS vs. GEFS+ vs. FHM). In contrast, patch-clamp data from mammalian expression systems revealed functional differences among missense variants allowing discrimination between disease severities. Those presenting with milder phenotypes retained a degree of channel function measured as residual whole-cell current, whereas those without any wholecell current were often associated with DS (p = .024).These findings demonstrate that electrophysiological data from mammalian expression systems can serve as useful disease biomarker when evaluating SCN1A variants, particularly in view of new and emerging treatment options in DS. K E Y W O R D SDravet syndrome, electrophysiology, familial hemiplegic migraine, functional testing, GEFS+, patch-clamp, SCN1A *Andreas Brunklaus and Stephanie Schorge contributed equally to this study.

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.

Dupuytren´s disease, a fibromatosis of the connective tissue in the palm, is a common complex disease with a strong genetic component. Up to date nine genetic loci have been found to be associated with the disease. Six of these loci contain genes that code for Wnt signalling proteins. In spite of this striking first insight into the genetic factors in Dupuytren´s disease, much of the inherited risk in Dupuytren´s disease still needs to be discovered. The already identified loci jointly explain ~1% of the heritability in this disease. To further elucidate the genetic basis of Dupuytren´s disease, we performed a genome-wide meta-analysis combining three genome-wide association study (GWAS) data sets, comprising 1,580 cases and 4,480 controls. We corroborated all nine previously identified loci, six of these with genome-wide significance (p-value < 5x10-8). In addition, we identified 14 new suggestive loci (p-value < 10−5). Intriguingly, several of these new loci contain genes associated with Wnt signalling and therefore represent excellent candidates for replication. Next, we compared whole-transcriptome data between patient- and control-derived tissue samples and found the Wnt/β-catenin pathway to be the top deregulated pathway in patient samples. We then conducted network and pathway analyses in order to identify protein networks that are enriched for genes highlighted in the GWAS meta-analysis and expression data sets. We found further evidence that the Wnt signalling pathways in conjunction with other pathways may play a critical role in Dupuytren´s disease.

No abstract

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 2,871 gene family protein sequence alignments involving 9,990 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 65,034 missense variants from patients. With this gene family approach, we identified 398 regions enriched for patient variants spanning 33,887 amino acids in 1,058 genes. As a comparison, testing the same genes individually we identified less patient variant enriched regions involving only 2,167 amino acids and 180 genes. Next, we selected de novo variants from 6,753 patients with neurodevelopmental disorders and 1,911 unaffected siblings, and observed a 5.56-fold enrichment of patient variants in our identified regions (95% C.I. =2.76-Inf, p-value = 6.66x10 -8 ). Using an independent ClinVar variant set, we found missense variants inside the identified regions are 111-fold more likely to be classified as pathogenic in comparison to benign classification (OR = 111.48, 95% C.I = 68. 09-195.58, p-value < 2.2e -16 ). All patient variant enriched regions identified (PERs) are available online through a user-friendly platform for interactive data mining, visualization and download at http://per.broadinstitute.org. In summary, our gene family burden analysis approach identified novel patient variant enriched regions in protein sequences. This annotation can empower variant interpretation.

Summary Objective Increasing availability of surgically resected brain tissue from patients with focal epilepsy and focal cortical dysplasia or low‐grade glioneuronal tumors has fostered large‐scale genetic examination. However, assessment of pathogenicity of germ line and somatic variants remains difficult. Here, we present a state‐of‐the‐art evaluation of reported genes and variants associated with epileptic brain lesions. Methods We critically reevaluated the pathogenicity for all neuropathology‐associated variants reported to date in the PubMed and ClinVar databases, including 101 neuropathology‐associated missense variants encompassing 11 disease‐related genes. We assessed gene variant tolerance and classified all identified missense variants according to guidelines from the American College of Medical Genetics and Genomics (ACMG). We further extended the bioinformatic variant prediction by introducing a novel gene‐specific deleteriousness ranking for prediction scores. Results Application of ACMG guidelines and in silico gene variant tolerance analysis classified only seven of 11 genes to be likely disease‐associated according to the reported disease mechanism, whereas 61 (60.4%) of 101 variants of those genes were classified as of uncertain significance, 37 (36.6%) as being likely pathogenic, and 3 (3%) as being pathogenic. Significance We concluded that the majority of neuropathology‐associated variants reported to date do not have enough evidence to be classified as pathogenic. Interpretation of lesion‐associated variants remains challenging, and application of current ACMG guidelines is recommended for interpretation and prediction.

Motivation The correct classification of missense variants as benign or pathogenic remains challenging. Pathogenic variants are expected to have higher deleterious prediction scores than benign variants in the same gene. However, most of the existing variant annotation tools do not reference the score range of benign population variants on gene level. Results We present a web-application, Variant Score Ranker, which enables users to rapidly annotate variants and perform gene-specific variant score ranking on the population level. We also provide an intuitive example of how gene- and population-calibrated variant ranking scores can improve epilepsy variant prioritization. Availability and implementation http://vsranker.broadinstitute.org Supplementary information Supplementary data are available at Bioinformatics online.