A clustering approach to improve our understanding of the genetic and phenotypic complexity of chronic kidney disease

Chronic kidney disease (CKD) is a complex disorder that causes a gradual loss of kidney function, affecting approximately 9.1% of the world's population. Here, we use a soft-clustering algorithm to deconstruct its genetic heterogeneity. First, we selected 322 CKD-associated independent genetic variants from published genome-wide association studies (GWAS) and added association results for 229 traits from the GWAS catalog. We then applied nonnegative matrix factorization (NMF) to discover overlapping clusters of related traits and variants. We computed cluster-specific polygenic scores and validated each cluster with a phenome-wide association study (PheWAS) on the BioMe biobank (n = 31,701). NMF identified nine clusters that reflect different aspects of CKD, with the top-weighted traits signifying areas such as kidney function, type 2 diabetes (T2D), and body weight. For most clusters, the top-weighted traits were confirmed in the PheWAS analysis. Results were found to be more significant in the cross-ancestry analysis, although significant ancestry-specific associations were also identified. While all alleles were associated with a decreased kidney function, associations with CKD-related diseases (e.g., T2D) were found only for a smaller subset of variants and differed across genetic ancestry groups. Our findings leverage genetics to gain insights into the underlying biology of CKD and investigate population-specific associations.

Cases of trisomy 21 and trisomy 18 among historic and prehistoric individuals discovered from ancient DNA

Aneuploidies, and in particular, trisomies represent the most common genetic aberrations observed in human genetics today. To explore the presence of trisomies in historic and prehistoric populations we screen nearly 10,000 ancient human individuals for the presence of three copies of any of the target autosomes. We find clear genetic evidence for six cases of trisomy 21 (Down syndrome) and one case of trisomy 18 (Edwards syndrome), and all cases are present in infant or perinatal burials. We perform comparative osteological examinations of the skeletal remains and find overlapping skeletal markers, many of which are consistent with these syndromes. Interestingly, three cases of trisomy 21, and the case of trisomy 18 were detected in two contemporaneous sites in early Iron Age Spain (800-400 BCE), potentially suggesting a higher frequency of burials of trisomy carriers in those societies. Notably, the care with which the burials were conducted, and the items found with these individuals indicate that ancient societies likely acknowledged these individuals with trisomy 18 and 21 as members of their communities, from the perspective of burial practice.

Polygenic risk scores as a marker for epilepsy risk across lifetime and after unspecified seizure events

A diagnosis of epilepsy has significant consequences for an individual but is often challenging in clinical practice. Novel biomarkers are thus greatly needed. Here, we investigated how common genetic factors (epilepsy polygenic risk scores, [PRSs]) influence epilepsy risk in detailed longitudinal electronic health records (EHRs) of > 360k Finns spanning up to 50 years of individuals' lifetimes. Individuals with a high genetic generalized epilepsy PRS (PRSGGE) in FinnGen had an increased risk for genetic generalized epilepsy (GGE) (hazard ratio [HR] 1.55 per PRSGGE standard deviation [SD]) across their lifetime and after unspecified seizure events. Effect sizes of epilepsy PRSs were comparable to effect sizes in clinically curated data supporting our EHR-derived epilepsy diagnoses. Within 10 years after an unspecified seizure, the GGE rate was 37% when PRSGGE > 2 SD compared to 5.6% when PRSGGE < -2 SD. The effect of PRSGGE was even larger on GGE subtypes of idiopathic generalized epilepsy (IGE) (HR 2.1 per SD PRSGGE). We further report significantly larger effects of PRSGGE on epilepsy in females and in younger age groups. Analogously, we found significant but more modest focal epilepsy PRS burden associated with non-acquired focal epilepsy (NAFE). We found PRSGGE specifically associated with GGE in comparison with >2000 independent diseases while PRSNAFE was also associated with other diseases than NAFE such as back pain. Here, we show that epilepsy specific PRSs have good discriminative ability after a first seizure event i.e. in circumstances where the prior probability of epilepsy is high outlining a potential to serve as biomarkers for an epilepsy diagnosis.

GWAS meta-analysis of over 29,000 people with epilepsy identifies 26 risk loci and subtype-specific genetic architecture

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here we report a multi-ancestry genome-wide association study including 29,944 cases, stratified into three broad categories and seven subtypes of epilepsy, and 52,538 controls. We identify 26 genome-wide significant loci, 19 of which are specific to genetic generalized epilepsy (GGE). We implicate 29 likely causal genes underlying these 26 loci. SNP-based heritability analyses show that common variants explain between 39.6% and 90% of genetic risk for GGE and its subtypes. Subtype analysis revealed markedly different genetic architectures between focal and generalized epilepsies. Gene-set analyses of GGE signals implicate synaptic processes in both excitatory and inhibitory neurons in the brain. Prioritized candidate genes overlap with monogenic epilepsy genes and with targets of current antiseizure medications. Finally, we leverage our results to identify alternate drugs with predicted efficacy if repurposed for epilepsy treatment.

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.

Implementation and evaluation of personal genetic testing as part of genomics analysis courses in German universities

PURPOSE

Due to the increasing application of genome analysis and interpretation in medical disciplines, professionals require adequate education. Here, we present the implementation of personal genotyping as an educational tool in two genomics courses targeting Digital Health students at the Hasso Plattner Institute (HPI) and medical students at the Technical University of Munich (TUM).

METHODS

We compared and evaluated the courses and the students' perceptions on the course setup using questionnaires.

RESULTS

During the course, students changed their attitudes towards genotyping (HPI: 79% [15 of 19], TUM: 47% [25 of 53]). Predominantly, students became more critical of personal genotyping (HPI: 73% [11 of 15], TUM: 72% [18 of 25]) and most students stated that genetic analyses should not be allowed without genetic counseling (HPI: 79% [15 of 19], TUM: 70% [37 of 53]). Students found the personal genotyping component useful (HPI: 89% [17 of 19], TUM: 92% [49 of 53]) and recommended its inclusion in future courses (HPI: 95% [18 of 19], TUM: 98% [52 of 53]).

CONCLUSION

Students perceived the personal genotyping component as valuable in the described genomics courses. The implementation described here can serve as an example for future courses in Europe.

The role of rare genetic variants enrichment in epilepsies of presumed genetic etiology

Previous studies suggested that severe epilepsies e.g., developmental and epileptic encephalopathies (DEE) are mainly caused by ultra-rare de novo genetic variants. For milder phenotypes, 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 controls. 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 controls at a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD≥20), and to have an odds ratio in epilepsy cases ≥2. We identified genes enriched with QRVs in DEE (n=21), NAFE (n=72), and GGE (n=32) - the number of enriched genes are found greatest in NAFE and least in DEE. This suggests that rare variants may play a more important role for causality of NAFE than in DEE. Moreover, we found that QRV-carrying genes e.g., HSGP2, FLNA or TNC are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE, while in DEE and GGE, the contribution of such variants appears more limited.

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 explored^1,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.

Mono- and biallelic variant effects on disease at biobank scale

Identifying causal factors for Mendelian and common diseases is an ongoing challenge in medical genetics¹. Population bottleneck events, such as those that occurred in the history of the Finnish population, enrich some homozygous variants to higher frequencies, which facilitates the identification of variants that cause diseases with recessive inheritance^2,3. Here we examine the homozygous and heterozygous effects of 44,370 coding variants on 2,444 disease phenotypes using data from the nationwide electronic health records of 176,899 Finnish individuals. We find associations for homozygous genotypes across a broad spectrum of phenotypes, including known associations with retinal dystrophy and novel associations with adult-onset cataract and female infertility. Of the recessive disease associations that we identify, 13 out of 20 would have been missed by the additive model that is typically used in genome-wide association studies. We use these results to find many known Mendelian variants whose inheritance cannot be adequately described by a conventional definition of dominant or recessive. In particular, we find variants that are known to cause diseases with recessive inheritance with significant heterozygous phenotypic effects. Similarly, we find presumed benign variants with disease effects. Our results show how biobanks, particularly in founder populations, can broaden our understanding of complex dosage effects of Mendelian variants on disease.

Polygenic risk scores in epilepsy

An epilepsy diagnosis has large consequences for an individual but is often difficult to make in clinical practice. Novel biomarkers are thus greatly needed. Here, we give an overview of how thousands of common genetic factors that increase the risk for epilepsy can be summarized as epilepsy polygenic risk scores (PRS). We discuss the current state of research on how epilepsy PRS can serve as a biomarker for the risk for epilepsy. The high heritability of common forms of epilepsy, particularly genetic generalized epilepsy, indicates a promising potential for epilepsy PRS in diagnosis and risk prediction. Small sample sizes and low ancestral diversity of current epilepsy genome-wide association studies show, however, a need for larger and more diverse studies before epilepsy PRS could be properly implemented in the clinic.

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

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel Nav1.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 Nav1.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.

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.

The origin and legacy of the Etruscans through a 2000-year archeogenomic time transect

The origin, development, and legacy of the enigmatic Etruscan civilization from the central region of the Italian peninsula known as Etruria have been debated for centuries. Here we report a genomic time transect of 82 individuals spanning almost two millennia (800 BCE to 1000 CE) across Etruria and southern Italy. During the Iron Age, we detect a component of Indo-European–associated steppe ancestry and the lack of recent Anatolian-related admixture among the putative non–Indo-European–speaking Etruscans. Despite comprising diverse individuals of central European, northern African, and Near Eastern ancestry, the local gene pool is largely maintained across the first millennium BCE. This drastically changes during the Roman Imperial period where we report an abrupt population-wide shift to ~50% admixture with eastern Mediterranean ancestry. Last, we identify northern European components appearing in central Italy during the Early Middle Ages, which thus formed the genetic landscape of present-day Italian populations.

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-Lepr^{L536Hfs*6-1NKB}), which is located at chromosome 4, exon 11 within the CRH2-leptin-binding site. Compared with C57BL/6N mice, Lepr^L536Hfs*6 develop early onset obesity and their body weight exceeds that of Lepr^db/db mice at an age of 30 weeks. Similar to Lepr^db/db mice, the Lepr^L536Hfs*6 model is characterized by hyperphagia, obesity, lower energy expenditure and activity, hyperglycemia, and hyperinsulinemia compared with C57BL/6N mice. Crossing Lepr^db/wt with Lepr^L536Hfs*6/wt mice results in compound heterozygous Lepr^L536Hfs*6/db mice, which develop even higher body weight and fat mass than both homozygous Lepr^db/db and Lepr^L536Hfs*6 mice. Compound heterozygous Lepr deficiency affecting functionally different regions of the Lepr causes more severe obesity than the parental homozygous mutations.

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 ((AD)SHE) to epilepsy of infancy with migrating focal seizures (EIMFS) and include developmental and epileptic encephalopathies (DEE). 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 unpreviously published 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) DEE other than EIMFS (non-EIMFS DEE) (37 individuals, 17 unpublished); iii) (AD)SHE (53 patients, 14 unpublished); iv) other phenotypes (6 individuals, 2 unpublished). In our cohort of 66 new cases, the most common phenotypic features were: a) in EIMFS, heterogeneity of seizure types, including epileptic spasms, epilepsy improvement over time, no epilepsy-related deaths; b) in non-EIMFS DEE, possible onset with West syndrome, occurrence of atypical absences, possible evolution to DEE with SHE features; one case of sudden unexplained death in epilepsy (SUDEP); c) in (AD)SHE, 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 about 50% of the patients, SUDEP in one individual; d) 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 (AD)SHE-associated mutations to be clustered around the RCK2 domain in the C-terminus, distal to the NADP domain. Mutations associated with EIMFS/non-EIMFS DEE 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 DEEs as well as in focal epilepsies, namely (AD)SHE.

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.

Predicting functional effects of missense variants in voltage-gated sodium and calcium channels

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes.

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.

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.

Neolithic and medieval virus genomes reveal complex evolution of hepatitis B

The hepatitis B virus (HBV) is one of the most widespread human pathogens known today, yet its origin and evolutionary history are still unclear and controversial. Here, we report the analysis of three ancient HBV genomes recovered from human skeletons found at three different archaeological sites in Germany. We reconstructed two Neolithic and one medieval HBV genome by de novo assembly from shotgun DNA sequencing data. Additionally, we observed HBV-specific peptides using paleo-proteomics. Our results demonstrated that HBV has circulated in the European population for at least 7000 years. The Neolithic HBV genomes show a high genomic similarity to each other. In a phylogenetic network, they do not group with any human-associated HBV genome and are most closely related to those infecting African non-human primates. The ancient viruses appear to represent distinct lineages that have no close relatives today and possibly went extinct. Our results reveal the great potential of ancient DNA from human skeletons in order to study the long-time evolution of blood borne viruses.

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.

Delineating the GRIN1 phenotypic spectrum: A distinct genetic NMDA receptor encephalopathy

Objective:

To determine the phenotypic spectrum caused by mutations in GRIN1 encoding the NMDA receptor subunit GluN1 and to investigate their underlying functional pathophysiology.

Methods:

We collected molecular and clinical data from several diagnostic and research cohorts. Functional consequences of GRIN1 mutations were investigated in Xenopus laevis oocytes.

Results:

We identified heterozygous de novo GRIN1 mutations in 14 individuals and reviewed the phenotypes of all 9 previously reported patients. These 23 individuals presented with a distinct phenotype of profound developmental delay, severe intellectual disability with absent speech, muscular hypotonia, hyperkinetic movement disorder, oculogyric crises, cortical blindness, generalized cerebral atrophy, and epilepsy. Mutations cluster within transmembrane segments and result in loss of channel function of varying severity with a dominant-negative effect. In addition, we describe 2 homozygous GRIN1 mutations (1 missense, 1 truncation), each segregating with severe neurodevelopmental phenotypes in consanguineous families.

Conclusions:

De novo GRIN1 mutations are associated with severe intellectual disability with cortical visual impairment as well as oculomotor and movement disorders being discriminating phenotypic features. Loss of NMDA receptor function appears to be the underlying disease mechanism. The identification of both heterozygous and homozygous mutations blurs the borders of dominant and recessive inheritance of GRIN1-associated disorders.