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.

SCN1A variants from bench to bedside-improved clinical prediction from functional characterization.

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 whole-cell 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.

Variant Score Ranker-a web application for intuitive missense variant prioritization.

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.

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.

Guideline-based and bioinformatic reassessment of lesion-associated gene and variant pathogenicity in focal human epilepsies.

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.

Differential excitatory vs inhibitory SCN expression at single cell level regulates brain sodium channel function in neurodevelopmental disorders.

The four voltage-gated sodium channels SCN1/2/3/8A have been associated with heterogeneous types of developmental disorders, each presenting with disease specific temporal and cell type specific gene expression. Using single-cell RNA sequencing transcriptomic data from humans and mice, we observe that SCN1A is predominantly expressed in inhibitory neurons. In contrast, SCN2/3/8A are profoundly expressed in excitatory neurons with SCN2/3A starting prenatally, followed by SCN1/8A neonatally. In contrast to previous observations from low resolution RNA screens, we observe that all four genes are expressed in both excitatory and inhibitory neurons, however, exhibit differential expression strength. These findings provide molecular evidence, at single-cell resolution, to support the hypothesis that the excitatory/inhibitory (E/I) neuronal expression ratios of sodium channels are important regulatory mechanisms underlying brain homeostasis and neurological diseases. Modulating the E/I expression balance within cell types of sodium channels could serve as a potential strategy to develop targeted treatment for NaV-associated neuronal developmental disorders.

Meta-Analysis of Genome-Wide Association Studies and Network Analysis-Based Integration with Gene Expression Data Identify New Suggestive Loci and Unravel a Wnt-Centric Network Associated with Dupuytren's Disease.

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.