We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na⁺ channel Na_v1.6, with the aim of describing clinical phenotypes related to functional effects. Six different clinical subgroups were identified: Group 1, benign familial infantile epilepsy (n = 15, normal cognition, treatable seizures); Group 2, intermediate epilepsy (n = 33, mild intellectual disability, partially pharmaco-responsive); Group 3, developmental and epileptic encephalopathy (n = 177, severe intellectual disability, majority pharmaco-resistant); Group 4, generalized epilepsy (n = 20, mild to moderate intellectual disability, frequently with absence seizures); Group 5, unclassifiable epilepsy (n = 127); and Group 6, neurodevelopmental disorder without epilepsy (n = 20, mild to moderate intellectual disability). Those in Groups 1–3 presented with focal or multifocal seizures (median age of onset: 4 months) and focal epileptiform discharges, whereas the onset of seizures in patients with generalized epilepsy was later (median: 42 months) with generalized epileptiform discharges. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin-insensitive human Na_v1.6 channels and whole-cell patch-clamping. Two variants causing developmental and epileptic encephalopathy showed a strong gain-of-function (hyperpolarizing shift of steady-state activation, strongly increased neuronal firing rate) and one variant causing benign familial infantile epilepsy or intermediate epilepsy showed a mild gain-of-function (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (reduced current amplitudes, depolarizing shift of steady-state activation, reduced neuronal firing). Functional effects were known for 170 individuals. All 136 individuals carrying a functionally tested gain-of-function variant had either focal (n = 97, Groups 1–3) or unclassifiable (n = 39) epilepsy, whereas 34 individuals with a loss-of-function variant had either generalized (n = 14), no (n = 11) or unclassifiable (n = 6) epilepsy; only three had developmental and epileptic encephalopathy. Computational modelling in the gain-of-function group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. Gain-of-function variant carriers responded significantly better to sodium channel blockers than to other anti-seizure medications, and the same applied for all individuals in 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 loss-of-function variant carriers and the extent of the electrophysiological dysfunction of the gain-of-function variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that sodium channel blockers present a treatment option in SCN8A-related focal epilepsy with onset in the first year of life.

Major Facilitator Superfamily Domain containing 2a (MFSD2A) is an essential endothelial lipid transporter at the blood-brain barrier. Biallelic variants affecting function in MFSD2A cause autosomal recessive primary microcephaly 15 (MCPH15, OMIM# 616486). We sought to expand our knowledge of the phenotypic spectrum of MCPH15 and demonstrate the underlying mechanism of inactivation of the MFSD2A transporter. We carried out detailed analysis of the clinical and neuroradiological features of a series of 27 MCPH15 cases, including eight new individuals from seven unrelated families. Genetic investigation was performed through exome sequencing (ES). Structural insights on the human Mfsd2a model and in-vitro biochemical assays were used to investigate the functional impact of the identified variants. All patients had primary microcephaly and severe developmental delay. Brain MRI showed variable degrees of white matter reduction, ventricular enlargement, callosal hypodysgenesis, and pontine and vermian hypoplasia. ES led to the identification of six novel biallelic MFSD2A variants (NG_053084.1, NM_032793.5: c.556+1G>A, c.748G>T; p.(Val250Phe), c.750_753del; p.(Cys251SerfsTer3), c.977G>A; p.(Arg326His), c.1386_1435del; p.(Gln462HisfsTer17), and c.1478C>T; p.(Pro493Leu)) and two recurrent variants (NM_032793.5: c.593C>T; p.(Thr198Met) and c.476C>T; p.(Thr159Met)). All these variants and the previously reported NM_032793.5: c.490C>A; p.(Pro164Thr) resulted in either reduced MFSD2A expression and/or transport activity. Our study further delineates the phenotypic spectrum of MCPH15, refining its clinical and neuroradiological characterization and supporting that MFSD2A deficiency causes early prenatal brain developmental disruption. We also show that poor MFSD2A expression despite normal transporter activity is a relevant pathomechanism in MCPH15.

POGZ-related disorders (also known as White-Sutton syndrome) encompass a wide range of neurocognitive abnormalities and other accompanying anomalies. Disease severity varies widely among POGZ patients and studies investigating genotype-phenotype association are scarce. Therefore, our aim was to collect data on previously unreported POGZ patients and perform a large-scale phenotype-genotype comparison from published data. Overall, 117 POGZ patients’ genotype and phenotype data were included in the analysis, including 12 novel patients. A severity scoring system was developed for the comparison. Mild and severe phenotypes were compared with the types and location of the variants and the predicted presence or absence of nonsense-mediated RNA decay (NMD). Missense variants were more often associated with mild phenotypes (p = 0.0421) and truncating variants predicted to escape NMD presented with more severe phenotypes (p < 0.0001). Within this group, variants in the prolin-rich region of the POGZ protein were associated with the most severe phenotypes (p = 0.0004). Our study suggests that gain-of-function or dominant negative effect through escaping NMD and the location of the variants in the prolin-rich domain of the protein may play an important role in the severity of manifestations of POGZ–associated neurodevelopmental disorders.

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are needed to understand the disease mechanism and improve therapeutic decision-making. We aggregated genetic and clinical data from 172 individuals with likely pathogenic/pathogenic (lp/p) SLC6A1 variants and functional data for 184 variants (14.1% lp/p). Clinical and functional data were available for a subset of 126 individuals. We explored the potential associations of variant positions on the GAT1 3D structure with variant pathogenicity, altered molecular function and phenotype severity using bioinformatic approaches. The GAT1 transmembrane domains 1, 6 and extracellular loop 4 (EL4) were enriched for patient over population variants. Across functionally tested missense variants (n = 156), the spatial proximity from the ligand was associated with loss-of-function in the GAT1 transporter activity. For variants with complete loss of in vitro GABA uptake, we found a 4.6-fold enrichment in patients having severe disease versus non-severe disease (P = 2.9 × 10⁻³, 95% confidence interval: 1.5–15.3). In summary, we delineated associations between the 3D structure and variant pathogenicity, variant function and phenotype in SLC6A1-related disorders. This knowledge supports biology-informed variant interpretation and research on GAT1 function. All our data can be interactively explored in the SLC6A1 portal (https://slc6a1-portal.broadinstitute.org/).

The most important factor that complicates the work of dysmorphologists is the significant phenotypic variability of the human face. Next-Generation Phenotyping (NGP) tools that assist clinicians with recognizing characteristic syndromic patterns are particularly challenged when confronted with patients from populations different from their training data. To that end, we systematically analyzed the impact of genetic ancestry on facial dysmorphism. For that purpose, we established the GestaltMatcher Database (GMDB) as a reference dataset for medical images of patients with rare genetic disorders from around the world. We collected 10,980 frontal facial images – more than a quarter previously unpublished - from 8,346 patients, representing 581 rare disorders. Although the predominant ancestry is still European (67%), data from underrepresented populations have been increased considerably via global collaborations (19% Asian and 7% African). This includes previously unpublished reports for more than 40% of the African patients. The NGP analysis on this diverse dataset revealed characteristic performance differences depending on the composition of training and test sets corresponding to genetic relatedness. For clinical use of NGP, incorporating non-European patients resulted in a profound enhancement of GestaltMatcher performance. The top-5 accuracy rate increased by +11.29%. Importantly, this improvement in delineating the correct disorder from a facial portrait was achieved without decreasing the performance on European patients. By design, GMDB complies with the FAIR principles by rendering the curated medical data findable, accessible, interoperable, and reusable. This means GMDB can also serve as data for training and benchmarking. In summary, our study on facial dysmorphism on a global sample revealed a considerable cross ancestral phenotypic variability confounding NGP that should be counteracted by international efforts for increasing data diversity. GMDB will serve as a vital reference database for clinicians and a transparent training set for advancing NGP technology.

GNE myopathy (GNEM) is a rare hereditary disease, but at the same time, it is the most common distal myopathy in several countries due to a founder effect of some pathogenic variants in the GNE gene. We collected the largest cohort of patients with GNEM from Russia and analyzed their mutational spectrum and clinical data. In our cohort, 10 novel variants were found, including 2 frameshift variants and 2 large deletions. One novel missense variant c.169_170delGCinsTT (p.(Ala57Phe)) was detected in 4 families in a homozygous state and in 3 unrelated patients in a compound heterozygous state. It was the second most frequent variant in our cohort. All families with this novel frequent variant were non-consanguineous and originated from the 3 neighboring areas in the European part of Russia. The clinical picture of the patients carrying this novel variant was typical, but the severity of clinical manifestation differed significantly. In our study, we reported two atypical cases expanding the phenotypic spectrum of GNEM. One female patient had severe quadriceps atrophy, hand joint contractures, keloid scars, and non-classical pattern on leg muscle magnetic resonance imaging, which was more similar to atypical collagenopathy rather than GNEM. Another patient initially had been observed with spinal muscular atrophy due to asymmetric atrophy of hand muscles and results of electromyography. The peculiar pattern of muscle involvement on magnetic resonance imaging consisted of pronounced changes in the posterior thigh muscle group with relatively spared muscles of the lower legs, apart from the soleus muscles. Different variants in the GNE gene were found in both atypical cases. Thus, our data expand the mutational and clinical spectrum of GNEM.