Uncovering myocardial infarction genetic signatures using GWAS exploration in Saudi and European cohorts

Genome-wide association studies (GWAS) have yielded significant insights into the genetic architecture of myocardial infarction (MI), although studies in non-European populations are still lacking. Saudi Arabian cohorts offer an opportunity to discover novel genetic variants impacting disease risk due to a high rate of consanguinity. Genome-wide genotyping (GWG), imputation and GWAS followed by meta-analysis were performed based on two independent Saudi Arabian studies comprising 3950 MI patients and 2324 non-MI controls. Meta-analyses were then performed with these two Saudi MI studies and the CardioGRAMplusC4D and UK BioBank GWAS as controls. Meta-analyses of the two Saudi MI studies resulted in 17 SNPs with genome-wide significance. Meta-analyses of all 4 studies revealed 66 loci with genome-wide significance levels of p < 5 × 10-8. All of these variants, except rs2764203, have previously been reported as MI-associated loci or to have high linkage disequilibrium with known loci. One SNP association in Shisa family member 5 (SHISA5) (rs11707229) was evident at a much higher frequency in the Saudi MI populations (> 12% MAF). In conclusion, our results replicated many MI associations, whereas in Saudi-only GWAS (meta-analyses), several new loci were implicated that require future validation and functional analyses.

Identification of candidate genes for developmental colour agnosia in a single unique family

Colour agnosia is a disorder that impairs colour knowledge (naming, recognition) despite intact colour perception. Previously, we have identified the first and only-known family with hereditary developmental colour agnosia. The aim of the current study was to explore genomic regions and candidate genes that potentially cause this trait in this family. For three family members with developmental colour agnosia and three unaffected family members CGH-array analysis and exome sequencing was performed, and linkage analysis was carried out using DominantMapper, resulting in the identification of 19 cosegregating chromosomal regions. Whole exome sequencing resulted in 11 rare coding variants present in all affected family members with developmental colour agnosia and absent in unaffected members. These variants affected genes that have been implicated in neural processes and functions (CACNA2D4, DDX25, GRINA, MYO15A) or that have an indirect link to brain function, development or disease (MAML2, STAU1, TMED3, RABEPK), and a remaining group lacking brain expression or involved in non-neural traits (DEPDC7, OR1J1, OR8D4). Although this is an explorative study, the small set of candidate genes that could serve as a starting point for unravelling mechanisms of higher level cognitive functions and cortical specialization, and disorders therein such as developmental colour agnosia.

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.

Widespread genomic influences on phenotype in Dravet syndrome, a 'monogenic' condition

Dravet syndrome is an archetypal rare severe epilepsy, considered 'monogenic', typically caused by loss-of-function SCN1A variants. Despite a recognizable core phenotype, its marked phenotypic heterogeneity is incompletely explained by differences in the causal SCN1A variant or clinical factors. In 34 adults with SCN1A-related Dravet syndrome, we show additional genomic variation beyond SCN1A contributes to phenotype and its diversity, with an excess of rare variants in epilepsy-related genes as a set and examples of blended phenotypes, including one individual with an ultra-rare DEPDC5 variant and focal cortical dysplasia. The polygenic risk score for intelligence was lower, and for longevity, higher, in Dravet syndrome than in epilepsy controls. The causal, major-effect, SCN1A variant may need to act against a broadly compromised genomic background to generate the full Dravet syndrome phenotype, whilst genomic resilience may help to ameliorate the risk of premature mortality in adult Dravet syndrome survivors.

Increased prime edit rates in KCNQ2 and SCN1A via single nicking all-in-one plasmids

BACKGROUND

Prime editing (PE) is the most recent gene editing technology able to introduce targeted alterations to the genome, including single base pair changes, small insertions, and deletions. Several improvements to the PE machinery have been made in the past few years, and these have been tested in a range of model systems including immortalized cell lines, stem cells, and animal models. While double nicking RNA (dncRNA) PE systems PE3 and PE5 currently show the highest editing rates, they come with reduced accuracy as undesired indels or SNVs arise at edited loci. Here, we aimed to improve single ncRNA (sncRNA) systems PE2 and PE4max by generating novel all-in-one (pAIO) plasmids driven by an EF-1α promoter, which is especially suitable for human-induced pluripotent stem cell (hiPSC) models.

RESULTS

pAIO-EF1α-PE2 and pAIO-EF1α-PE4max were used to edit the voltage gated potassium channel gene KCNQ2 and voltage gated sodium channel gene SCN1A. Two clinically relevant mutations were corrected using pAIO-EF1α-PE2 including the homozygous truncating SCN1A R612* variant in HEK293T cells and the heterozygous gain-of-function KCNQ2 R201C variant in patient-derived hiPSC. We show that sncRNA PE yielded detectable editing rates in hiPSC ranging between 6.4% and 9.8%, which was further increased to 41% after a GFP-based fluorescence-activated cell sorting (FACS) cell sorting step. Furthermore, we show that selecting the high GFP expressing population improved editing efficiencies up to 3.2-fold compared to the low GFP expressing population, demonstrating that not only delivery but also the number of copies of the PE enzyme and/or pegRNA per cell are important for efficient editing. Edit rates were not improved when an additional silent protospacer-adjacent motif (PAM)-removing alteration was introduced in hiPSC at the target locus. Finally, there were no genome-wide off-target effects using pAIO-EF1α-PE2 and no off-target editing activity near the edit locus highlighting the accuracy of snc prime editors.

CONCLUSION

Taken together, our study shows an improved efficacy of EF-1α driven sncRNA pAIO-PE plasmids in hiPSC reaching high editing rates, especially after FACS sorting. Optimizing these sncRNA PE systems is of high value when considering future therapeutic in vivo use, where accuracy will be extremely important.

Distinct genetic basis of common epilepsies and structural magnetic resonance imaging measures

Focal and generalized epilepsies are associated with robust differences in magnetic resonance imaging (MRI) measures of subcortical structures, gray matter, and white matter. However, it is unknown whether such structural brain differences reflect the cause or consequence of epilepsy or its treatment. Analyses of common genetic variants underlying both common epilepsy risk and variability in structural brain measures can give further insights, as such inherited variants are not influenced by disease or treatment. Here, we performed genetic correlation analyses using data from the largest genome-wide association study (GWAS) on common epilepsy (n = 27 559 cases and 42 436 controls) and GWASs on MRI measures of white (n = 33 292) or gray matter (n = 51 665). We did not detect any significant genetic correlation between any type of common epilepsy and any of 280 measures of gray matter, white matter, or subcortical structures. These results suggest that there are distinct genetic bases underlying risk of common epilepsy and for structural brain measures. This would imply that the genetic basis of normal structural brain variation is unrelated to that of common epilepsy. Structural changes in epilepsy could rather be the consequence of epilepsy, its comorbidities, or its treatment, offering a cumulative record of disease.

Genome-wide copy number variant screening of Saudi schizophrenia patients reveals larger deletions in cases versus controls

Introduction

Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects.

Methods

Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs.

Results

CNVs in schizophrenia cases were twice as large on average than CNVs in controls (p = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions.

Discussion

Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not.

Using common genetic variants to find drugs for common epilepsies

Better drugs are needed for common epilepsies. Drug repurposing offers the potential of significant savings in the time and cost of developing new treatments. In order to select the best candidate drug(s) to repurpose for a disease, it is desirable to predict the relative clinical efficacy that drugs will have against the disease. Common epilepsy can be divided into different types and syndromes. Different antiseizure medications are most effective for different types and syndromes of common epilepsy. For predictions of antiepileptic efficacy to be clinically translatable, it is essential that the predictions are specific to each form of common epilepsy, and reflect the patterns of drug efficacy observed in clinical studies and practice. These requirements are not fulfilled by previously published drug predictions for epilepsy. We developed a novel method for predicting the relative efficacy of drugs against any common epilepsy, by using its Genome-Wide Association Study summary statistics and drugs' activity data. The methodological advancement in our technique is that the drug predictions for a disease are based upon drugs' effects on the function and abundance of proteins, and the magnitude and direction of those effects, relative to the importance, degree and direction of the proteins' dysregulation in the disease. We used this method to predict the relative efficacy of all drugs, licensed for any condition, against each of the major types and syndromes of common epilepsy. Our predictions are concordant with findings from real-world experience and randomized clinical trials. Our method predicts the efficacy of existing antiseizure medications against common epilepsies; in this prediction, our method outperforms the best alternative existing method: area under receiver operating characteristic curve (mean ± standard deviation) 0.83 ± 0.03 and 0.63 ± 0.04, respectively. Importantly, our method predicts which antiseizure medications are amongst the more efficacious in clinical practice, and which antiseizure medications are amongst the less efficacious in clinical practice, for each of the main syndromes of common epilepsy, and it predicts the distinct order of efficacy of individual antiseizure medications in clinical trials of different common epilepsies. We identify promising candidate drugs for each of the major syndromes of common epilepsy. We screen five promising predicted drugs in an animal model: each exerts a significant dose-dependent effect upon seizures. Our predictions are a novel resource for selecting suitable candidate drugs that could potentially be repurposed for each of the major syndromes of common epilepsy. Our method is potentially generalizable to other complex diseases.

Gene Therapies for Monogenic Autism Spectrum Disorders

Novel genome editing and transient gene therapies have been developed the past ten years, resulting in the first in-human clinical trials for monogenic disorders. Syndromic autism spectrum disorders can be caused by mutations in a single gene. Given the monogenic aspect and severity of syndromic ASD, it is an ideal candidate for gene therapies. Here, we selected 11 monogenic ASD syndromes, validated by animal models, and reviewed current gene therapies for each syndrome. Given the wide variety and novelty of some forms of gene therapy, the best possible option must be decided based on the gene and mutation.

CRISPRa-Mediated Upregulation of scn1laa During Early Development Causes Epileptiform Activity and dCas9-Associated Toxicity

Dravet syndrome (DS) is a monogenic epileptic encephalopathy caused by loss-of-function mutations in the voltage-gated sodium channel (VGSC) gene SCN1A. DS has an age of onset within the first year of life and severe disease prognosis. In the past years, it has been shown that upregulation of endogenous SCN1A can be beneficial in animal models for DS, but a complete rescue was not observed. We hypothesized that upregulation during early development that precedes onset of first symptoms might improve disease outcome. To test this hypothesis, we first evaluated the CRISPR activating method for early upregulation of voltage gated sodium channels during early development. We injected CRISPRa components, which target the proximal or distal promoter region of the VGSC gene scn1Laa in the yolk of one-cell stage zebrafish embryos. The effect of both dCas9-VPR and dCas9-VP64 was evaluated. Both CRISPRa fusions showed toxicity in the majority of embryos, with or without guide RNAs. The few embryos that survived developed normally, and dCas9-VPR induces an upregulation of scn1Laa mRNA until 24 hours after fertilization. At 5 days post fertilization, CRISPRa-injected embryos showed an epileptic phenotype, including locomotor burst movements, hyperactivity, and epileptiform activity originating from the brain. In addition to previously published scn1Laa and scn1Lab loss-of-function models, we conclude that gain of scn1Laa function can have an equally severe phenotype. Upregulation of scn1Laa in the current zebrafish model for DS, scn1Lab-KO, aggravated the disease phenotype, highlighting that early-stage upregulation using CRISPRa can lead to both toxicity and a worsening of the disease phenotype.

Assessing the role of rare genetic variants in drug-resistant, non-lesional focal epilepsy

OBJECTIVE

Resistance to antiseizure medications (ASMs) is one of the major concerns in the treatment of epilepsy. Despite the increasing number of ASMs available, the proportion of individuals with drug-resistant epilepsy remains unchanged. In this study, we aimed to investigate the role of rare genetic variants in ASM resistance.

METHODS

We performed exome sequencing of 1,128 individuals with non-familial non-acquired focal epilepsy (NAFE) (762 non-responders, 366 responders) and were provided with 1,734 healthy controls. We undertook replication in a cohort of 350 individuals with NAFE (165 non-responders, 185 responders). We performed gene-based and gene-set-based kernel association tests to investigate potential enrichment of rare variants in relation to drug response status and to risk for NAFE.

RESULTS

We found no gene or gene set that reached genome-wide significance. Yet, we identified several prospective candidate genes - among them DEPDC5, which showed a potential association with resistance to ASMs. We found some evidence for an enrichment of truncating variants in dominant familial NAFE genes in our cohort of non-familial NAFE and in association with drug-resistant NAFE.

INTERPRETATION

Our study identifies potential candidate genes for ASM resistance. Our results corroborate the role of rare variants for non-familial NAFE and imply their involvement in drug-resistant epilepsy. Future large-scale genetic research studies are needed to substantiate these findings.

Role of Common Genetic Variants for Drug-Resistance to Specific Anti-Seizure Medications

Objective: Resistance to anti-seizure medications (ASMs) presents a significant hurdle in the treatment of people with epilepsy. Genetic markers for resistance to individual ASMs could support clinicians to make better-informed choices for their patients. In this study, we aimed to elucidate whether the response to individual ASMs was associated with common genetic variation. Methods: A cohort of 3,649 individuals of European descent with epilepsy was deeply phenotyped and underwent single nucleotide polymorphism (SNP)-genotyping. We conducted genome-wide association analyses (GWASs) on responders to specific ASMs or groups of functionally related ASMs, using non-responders as controls. We performed a polygenic risk score (PRS) analyses based on risk variants for epilepsy and neuropsychiatric disorders and ASM resistance itself to delineate the polygenic burden of ASM-specific drug resistance. Results: We identified several potential regions of interest but did not detect genome-wide significant loci for ASM-specific response. We did not find polygenic risk for epilepsy, neuropsychiatric disorders, and drug-resistance associated with drug response to specific ASMs or mechanistically related groups of ASMs. Significance: This study could not ascertain the predictive value of common genetic variants for ASM responder status. The identified suggestive loci will need replication in future studies of a larger scale.

Shared genetic basis between genetic generalized epilepsy and background electroencephalographic oscillations

Objective

Paroxysmal epileptiform abnormalities on electroencephalography (EEG) are the hallmark of epilepsies, but it is uncertain to what extent epilepsy and background EEG oscillations share neurobiological underpinnings. Here, we aimed to assess the genetic correlation between epilepsy and background EEG oscillations.

Methods

Confounding factors, including the heterogeneous etiology of epilepsies and medication effects, hamper studies on background brain activity in people with epilepsy. To overcome this limitation, we compared genetic data from a genome‐wide association study (GWAS) on epilepsy (n = 12 803 people with epilepsy and 24 218 controls) with that from a GWAS on background EEG (n = 8425 subjects without epilepsy), in which background EEG oscillation power was quantified in four different frequency bands: alpha, beta, delta, and theta. We replicated our findings in an independent epilepsy replication dataset (n = 4851 people with epilepsy and 20 428 controls). To assess the genetic overlap between these phenotypes, we performed genetic correlation analyses using linkage disequilibrium score regression, polygenic risk scores, and Mendelian randomization analyses.

Results

Our analyses show strong genetic correlations of genetic generalized epilepsy (GGE) with background EEG oscillations, primarily in the beta frequency band. Furthermore, we show that subjects with higher beta and theta polygenic risk scores have a significantly higher risk of having generalized epilepsy. Mendelian randomization analyses suggest a causal effect of GGE genetic liability on beta oscillations.

Significance

Our results point to shared biological mechanisms underlying background EEG oscillations and the susceptibility for GGE, opening avenues to investigate the clinical utility of background EEG oscillations in the diagnostic workup of epilepsy.

Symptomatology of carbamazepine- and oxcarbazepine-induced hyponatremia in people with epilepsy

OBJECTIVE

To ascertain whether adverse effects experienced by people taking carbamazepine or oxcarbazepine could be attributed to carbamazepine- or oxcarbazepine-induced hyponatremia (COIH).

METHODS

We performed an observational study, collecting data between 2017 and 2019 on serum sodium levels and adverse effects retrospectively in people with epilepsy while receiving treatment with either carbamazepine (CBZ) or oxcarbazepine (OXC). We defined hyponatremia as sodium level ≤134 mEq/L and severe hyponatremia as sodium level ≤128 mEq/L. Adverse effects experienced were compared between groups of individuals with and without hyponatremia.

RESULTS

A total of 1370 people using CBZ or OXC were identified, of whom 410 had at least one episode of hyponatremia. We checked for symptoms related to the use of CBZ and OXC in 710 people (410 with and 300 without hyponatremia) and found relevant information in 688. Adverse effects occurred in 65% of people with hyponatremia compared to 21% with normal sodium levels (odds ratio [OR] 7.5, P ≤ .001) and in 83% of people with severe hyponatremia compared to 55% in those with mild hyponatremia (P ≤ .001). Significant predictors of adverse effects were the drug (OXC vs CBZ), and the number of concomitant anti-seizure medications. Dizziness (28% vs 6%), tiredness (22% vs 7%), instability (19% vs 3%), and diplopia (16% vs 4%) were reported more often in the hyponatremia group than in patients with normal levels.

SIGNIFICANCE

People with COIH had a 7-fold increased risk of developing adverse effects during treatment. Clinicians should consider ascertainment of sodium levels in patients taking CBZ and OXC and act upon findings.

GATOR1-related focal cortical dysplasia in epilepsy surgery patients and their families: A possible gradient in severity?

BACKGROUND

Variants of GATOR1-genes represent a recognised cause of focal cortical dysplasia (FCD), the most common structural aetiology in paediatric drug-resistant focal epilepsy. Reports on familial cases of GATOR1-associated FCD are limited, especially with respect to epilepsy surgery outcomes.

METHODS

We present phenotypical manifestations of four unrelated patients with drug-resistant focal epilepsy, FCD and a first-degree relative with epilepsy. All patients underwent targeted gene panel sequencing as a part of the presurgical work up. Literature search was performed to compare our findings to previously published cases.

RESULTS

The children (probands) had a more severe phenotype than their parents, including drug-resistant epilepsy and developmental delay, and they failed to achieve seizure freedom post-surgically. All patients had histopathologically confirmed FCD (types IIa, IIb, Ia). In Patient 1 and her affected father, we detected a known pathogenic NPRL2 variant. In patients 2 and 3 and their affected parents, we found novel likely pathogenic germline DEPDC5 variants. In family 4, we detected a novel variant in NPRL3. We identified 15 additional cases who underwent epilepsy surgery for GATOR1-associated FCD, with a positive family history of epilepsy in the literature; in 8/13 tested, the variant was inherited from an asymptomatic parent.

CONCLUSION

The presented cases displayed a severity gradient in phenotype with children more severely affected than the parents. Although patients with GATOR1-associated FCD are considered good surgical candidates, post-surgical seizure outcome was poor in our familial cases, suggesting that accurate identification of the epileptogenic zone may be more challenging in this subgroup of patients.

NEXMIF encephalopathy: an X-linked disorder with male and female phenotypic patterns

PURPOSE

Pathogenic variants in the X-linked gene NEXMIF (previously KIAA2022) are associated with intellectual disability (ID), autism spectrum disorder, and epilepsy. We aimed to delineate the female and male phenotypic spectrum of NEXMIF encephalopathy.

METHODS

Through an international collaboration, we analyzed the phenotypes and genotypes of 87 patients with NEXMIF encephalopathy.

RESULTS

Sixty-three females and 24 males (46 new patients) with NEXMIF encephalopathy were studied, with 30 novel variants. Phenotypic features included developmental delay/ID in 86/87 (99%), seizures in 71/86 (83%) and multiple comorbidities. Generalized seizures predominated including myoclonic seizures and absence seizures (both 46/70, 66%), absence with eyelid myoclonia (17/70, 24%), and atonic seizures (30/70, 43%). Males had more severe developmental impairment; females had epilepsy more frequently, and varied from unaffected to severely affected. All NEXMIF pathogenic variants led to a premature stop codon or were deleterious structural variants. Most arose de novo, although X-linked segregation occurred for both sexes. Somatic mosaicism occurred in two males and a family with suspected parental mosaicism.

CONCLUSION

NEXMIF encephalopathy is an X-linked, generalized developmental and epileptic encephalopathy characterized by myoclonic-atonic epilepsy overlapping with eyelid myoclonia with absence. Some patients have developmental encephalopathy without epilepsy. Males have more severe developmental impairment. NEXMIF encephalopathy arises due to loss-of-function variants.

Analysis of shared common genetic risk between amyotrophic lateral sclerosis and epilepsy

Because hyper-excitability has been shown to be a shared pathophysiological mechanism, we used the latest and largest genome-wide studies in amyotrophic lateral sclerosis (n = 36,052) and epilepsy (n = 38,349) to determine genetic overlap between these conditions. First, we showed no significant genetic correlation, also when binned on minor allele frequency. Second, we confirmed the absence of polygenic overlap using genomic risk score analysis. Finally, we did not identify pleiotropic variants in meta-analyses of the 2 diseases. Our findings indicate that amyotrophic lateral sclerosis and epilepsy do not share common genetic risk, showing that hyper-excitability in both disorders has distinct origins.

Antiepileptic Drug Teratogenicity and De Novo Genetic Variation Load

OBJECTIVE

The mechanisms by which antiepileptic drugs (AEDs) cause birth defects (BDs) are unknown. Data suggest that AED-induced BDs may result from a genome-wide increase of de novo variants in the embryo, a mechanism that we investigated.

METHODS

Whole exome sequencing data from child-parent trios were interrogated for de novo single-nucleotide variants/indels (dnSNVs/indels) and de novo copy number variants (dnCNVs). Generalized linear models were applied to assess de novo variant burdens in children exposed prenatally to AEDs (AED-exposed children) versus children without BDs not exposed prenatally to AEDs (AED-unexposed unaffected children), and AED-exposed children with BDs versus those without BDs, adjusting for confounders. Fisher exact test was used to compare categorical data.

RESULTS

Sixty-seven child-parent trios were included: 10 with AED-exposed children with BDs, 46 with AED-exposed unaffected children, and 11 with AED-unexposed unaffected children. The dnSNV/indel burden did not differ between AED-exposed children and AED-unexposed unaffected children (median dnSNV/indel number/child [range] = 3 [0-7] vs 3 [1-5], p = 0.50). Among AED-exposed children, there were no significant differences between those with BDs and those unaffected. Likely deleterious dnSNVs/indels were detected in 9 of 67 (13%) children, none of whom had BDs. The proportion of cases harboring likely deleterious dnSNVs/indels did not differ significantly between AED-unexposed and AED-exposed children. The dnCNV burden was not associated with AED exposure or birth outcome.

INTERPRETATION

Our study indicates that prenatal AED exposure does not increase the burden of de novo variants, and that this mechanism is not a major contributor to AED-induced BDs. These results can be incorporated in routine patient counseling. ANN NEUROL 2020;87:897-906.

Gene family information facilitates variant interpretation and identification of disease-associated genes in neurodevelopmental disorders

BACKGROUND

Classifying pathogenicity of missense variants represents a major challenge in clinical practice during the diagnoses of rare and genetic heterogeneous neurodevelopmental disorders (NDDs). While orthologous gene conservation is commonly employed in variant annotation, approximately 80% of known disease-associated genes belong to gene families. The use of gene family information for disease gene discovery and variant interpretation has not yet been investigated on a genome-wide scale. We empirically evaluate whether paralog-conserved or non-conserved sites in human gene families are important in NDDs.

METHODS

Gene family information was collected from Ensembl. Paralog-conserved sites were defined based on paralog sequence alignments; 10,068 NDD patients and 2078 controls were statistically evaluated for de novo variant burden in gene families.

RESULTS

We demonstrate that disease-associated missense variants are enriched at paralog-conserved sites across all disease groups and inheritance models tested. We developed a gene family de novo enrichment framework that identified 43 exome-wide enriched gene families including 98 de novo variant carrying genes in NDD patients of which 28 represent novel candidate genes for NDD which are brain expressed and under evolutionary constraint.

CONCLUSION

This study represents the first method to incorporate gene family information into a statistical framework to interpret variant data for NDDs and to discover new NDD-associated genes.

NaV1.1 and NaV1.6 selective compounds reduce the behavior phenotype and epileptiform activity in a novel zebrafish model for Dravet Syndrome

Dravet syndrome is caused by dominant loss-of-function mutations in SCN1A which cause reduced activity of Nav1.1 leading to lack of neuronal inhibition. On the other hand, gain-of-function mutations in SCN8A can lead to a severe epileptic encephalopathy subtype by over activating NaV1.6 channels. These observations suggest that Nav1.1 and Nav1.6 represent two opposing sides of the neuronal balance between inhibition and activation. Here, we hypothesize that Dravet syndrome may be treated by either enhancing Nav1.1 or reducing Nav1.6 activity. To test this hypothesis we generated and characterized a novel DS zebrafish model and tested new compounds that selectively activate or inhibit the human NaV1.1 or NaV1.6 channel respectively. We used CRISPR/Cas9 to generate two separate Scn1Lab knockout lines as an alternative to previous zebrafish models generated by random mutagenesis or morpholino oligomers. Using an optimized locomotor assay, spontaneous burst movements were detected that were unique to Scn1Lab knockouts and disappear when introducing human SCN1A mRNA. Besides the behavioral phenotype, Scn1Lab knockouts show sudden, electrical discharges in the brain that indicate epileptic seizures in zebrafish. Scn1Lab knockouts showed increased sensitivity to the GABA antagonist pentylenetetrazole and a reduction in whole organism GABA levels. Drug screenings further validated a Dravet syndrome phenotype. We tested the NaV1.1 activator AA43279 and two novel NaV1.6 inhibitors MV1369 and MV1312 in the Scn1Lab knockouts. Both type of compounds significantly reduced the number of spontaneous burst movements and seizure activity. Our results show that selective inhibition of NaV1.6 could be just as efficient as selective activation of NaV1.1 and these approaches could prove to be novel potential treatment strategies for Dravet syndrome and other (genetic) epilepsies. Compounds tested in zebrafish however, should always be further validated in other model systems for efficacy in mammals and to screen for potential side effects.

Modifier genes in SCN1A-related epilepsy syndromes

Background

SCN1A is one of the most important epilepsy‐related genes, with pathogenic variants leading to a range of phenotypes with varying disease severity. Different modifying factors have been hypothesized to influence SCN1A‐related phenotypes. We investigate the presence of rare and more common variants in epilepsy‐related genes as potential modifiers of SCN1A‐related disease severity.

Methods

87 patients with SCN1A‐related epilepsy were investigated. Whole‐exome sequencing was performed by the Beijing Genomics Institute (BGI). Functional variants in 422 genes associated with epilepsy and/or neuronal excitability were investigated. Differences in proportions of variants between the epilepsy genes and four control gene sets were calculated, and compared to the proportions of variants in the same genes in the ExAC database.

Results

Statistically significant excesses of variants in epilepsy genes were observed in the complete cohort and in the combined group of mildly and severely affected patients, particularly for variants with minor allele frequencies of <0.05. Patients with extreme phenotypes showed much greater excesses of epilepsy gene variants than patients with intermediate phenotypes.

Conclusion

Our results indicate that relatively common variants in epilepsy genes, which would not necessarily be classified as pathogenic, may play a large role in modulating SCN1A phenotypes. They may modify the phenotypes of both severely and mildly affected patients. Our results may be a first step toward meaningful testing of modifier gene variants in regular diagnostics for individual patients, to provide a better estimation of disease severity for newly diagnosed patients.

Clinical and genetic correlates of islet-autoimmune signatures in juvenile-onset type 1 diabetes

AIMS/HYPOTHESIS

Heterogeneity in individuals with type 1 diabetes has become more generally appreciated, but has not yet been extensively and systematically characterised. Here, we aimed to characterise type 1 diabetes heterogeneity by creating immunological, genetic and clinical profiles for individuals with juvenile-onset type 1 diabetes in a cross-sectional study.

METHODS

Participants were HLA-genotyped to determine HLA-DR-DQ risk, and SNP-genotyped to generate a non-HLA genetic risk score (GRS) based on 93 type 1 diabetes-associated SNP variants outside the MHC region. Islet autoimmunity was assessed as T cell proliferation upon stimulation with the beta cell antigens GAD65, islet antigen-2 (IA-2), preproinsulin (PPI) and defective ribosomal product of the insulin gene (INS-DRIP). Clinical parameters were collected retrospectively.

RESULTS

Of 80 individuals, 67 had proliferation responses to one or more islet antigens, with vast differences in the extent of proliferation. Based on the multitude and amplitude of the proliferation responses, individuals were clustered into non-, intermediate and high responders. High responders could not be characterised entirely by enrichment for the highest risk HLA-DR3-DQ2/DR4-DQ8 genotype. However, high responders did have a significantly higher non-HLA GRS. Clinically, high T cell responses to beta cell antigens did not reflect in worsened glycaemic control, increased complications, development of associated autoimmunity or younger age at disease onset. The number of beta cell antigens that an individual responded to increased with disease duration, pointing to chronic islet autoimmunity and epitope spreading.

CONCLUSIONS/INTERPRETATION

Collectively, these data provide new insights into type 1 diabetes disease heterogeneity and highlight the importance of stratifying patients on the basis of their genetic and autoimmune signatures for immunotherapy and personalised disease management.

Genomic and clinical predictors of lacosamide response in refractory epilepsies

OBJECTIVE

Clinical and genetic predictors of response to antiepileptic drugs (AEDs) are largely unknown. We examined predictors of lacosamide response in a real-world clinical setting.

METHODS

We tested the association of clinical predictors with treatment response using regression modeling in a cohort of people with refractory epilepsy. Genetic assessment for lacosamide response was conducted via genome-wide association studies and exome studies, comprising 281 candidate genes.

RESULTS

Most patients (479/483) were treated with LCM in addition to other AEDs. Our results corroborate previous findings that patients with refractory genetic generalized epilepsy (GGE) may respond to treatment with LCM. No clear clinical predictors were identified. We then compared 73 lacosamide responders, defined as those experiencing greater than 75% seizure reduction or seizure freedom, to 495 nonresponders (<25% seizure reduction). No variants reached the genome-wide significance threshold in our case-control analysis.

SIGNIFICANCE

No genetic predictor of lacosamide response was identified. Patients with refractory GGE might benefit from treatment with lacosamide.

Implications of genetic diagnostics in epilepsy surgery candidates: A single-center cohort study

OBJECTIVE

Genetic causes are increasingly identified in patients with focal epilepsy. These genetic causes may be related to the effectiveness of epilepsy surgery. We aimed to assess the use and yield of genetic testing in a large cohort of patients who were evaluated for epilepsy surgery.

METHODS

We performed a retrospective single-center consecutive cohort study of patients who were evaluated for surgery between 1990 and 2016. Within this cohort, we assessed the use of genetic testing-either before or after presurgical decision-making. We evaluated genetic results as well as the outcome of presurgical decision-making and surgery, and compared these end points for different subgroups-especially MRI-positive vs MRI-negative patients. Patients with tuberous sclerosis (TSC) and KRIT1 mutations were excluded from analysis.

RESULTS

Of the 2385 epilepsy patients who were evaluated for surgery, 1280 (54%) received surgical treatment in our center. Of the entire cohort, 325 (14%) underwent genetic testing, comprising 156 of 450 MRI-negative patients (35%) vs 169 of 1935 MRI-positive patients (9%). A genetic cause of epilepsy was found in 40 of the 325 patients (12%, 2% of the entire cohort), mainly consisting of mutations in ion channel function and synaptic transmission genes, and mTOR pathway gene mutations. Three of the seven patients with mTOR pathway gene mutations underwent surgery; two achieved complete seizure freedom. One of the 17 patients with germline mutations in ion channel function and synaptic transmission genes received resective surgery but was not rendered seizure-free; two other patients underwent invasive intracranial EEG-monitoring before being rejected.

SIGNIFICANCE

This study shows that genetic testing is increasingly applied in focal epilepsy patients who are considered for epilepsy surgery. The diagnostic yield of genetic testing is highest in next generation sequencing techniques, and the outcome of genetic testing assists selecting eligible patients for invasive intracranial monitoring and resective surgery.

Polygenic burden in focal and generalized epilepsies

Rare genetic variants can cause epilepsy, and genetic testing has been widely adopted for severe, paediatric-onset epilepsies. The phenotypic consequences of common genetic risk burden for epilepsies and their potential future clinical applications have not yet been determined. Using polygenic risk scores (PRS) from a European-ancestry genome-wide association study in generalized and focal epilepsy, we quantified common genetic burden in patients with generalized epilepsy (GE-PRS) or focal epilepsy (FE-PRS) from two independent non-Finnish European cohorts (Epi25 Consortium, n = 5705; Cleveland Clinic Epilepsy Center, n = 620; both compared to 20 435 controls). One Finnish-ancestry population isolate (Finnish-ancestry Epi25, n = 449; compared to 1559 controls), two European-ancestry biobanks (UK Biobank, n = 383 656; Vanderbilt biorepository, n = 49 494), and one Japanese-ancestry biobank (BioBank Japan, n = 168 680) were used for additional replications. Across 8386 patients with epilepsy and 622 212 population controls, we found and replicated significantly higher GE-PRS in patients with generalized epilepsy of European-ancestry compared to patients with focal epilepsy (Epi25: P = 1.64×10-15; Cleveland: P = 2.85×10-4; Finnish-ancestry Epi25: P = 1.80×10-4) or population controls (Epi25: P = 2.35×10-70; Cleveland: P = 1.43×10-7; Finnish-ancestry Epi25: P = 3.11×10-4; UK Biobank and Vanderbilt biorepository meta-analysis: P = 7.99×10-4). FE-PRS were significantly higher in patients with focal epilepsy compared to controls in the non-Finnish, non-biobank cohorts (Epi25: P = 5.74×10-19; Cleveland: P = 1.69×10-6). European ancestry-derived PRS did not predict generalized epilepsy or focal epilepsy in Japanese-ancestry individuals. Finally, we observed a significant 4.6-fold and a 4.5-fold enrichment of patients with generalized epilepsy compared to controls in the top 0.5% highest GE-PRS of the two non-Finnish European cohorts (Epi25: P = 2.60×10-15; Cleveland: P = 1.39×10-2). We conclude that common variant risk associated with epilepsy is significantly enriched in multiple cohorts of patients with epilepsy compared to controls-in particular for generalized epilepsy. As sample sizes and PRS accuracy continue to increase with further common variant discovery, PRS could complement established clinical biomarkers and augment genetic testing for patient classification, comorbidity research, and potentially targeted treatment.

GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways

Systemic sclerosis (SSc) is an autoimmune disease that shows one of the highest mortality rates among rheumatic diseases. We perform a large genome-wide association study (GWAS), and meta-analysis with previous GWASs, in 26,679 individuals and identify 27 independent genome-wide associated signals, including 13 new risk loci. The novel associations nearly double the number of genome-wide hits reported for SSc thus far. We define 95% credible sets of less than 5 likely causal variants in 12 loci. Additionally, we identify specific SSc subtype-associated signals. Functional analysis of high-priority variants shows the potential function of SSc signals, with the identification of 43 robust target genes through HiChIP. Our results point towards molecular pathways potentially involved in vasculopathy and fibrosis, two main hallmarks in SSc, and highlight the spectrum of critical cell types for the disease. This work supports a better understanding of the genetic basis of SSc and provides directions for future functional experiments.

Assessing the genetic association between vitamin B6 metabolism and genetic generalized epilepsy

Altered vitamin B6 metabolism due to pathogenic variants in the gene PNPO causes early onset epileptic encephalopathy, which can be treated with high doses of vitamin B6. We recently reported that single nucleotide polymorphisms (SNPs) that influence PNPO expression in the brain are associated with genetic generalized epilepsy (GGE). However, it is not known whether any of these GGE-associated SNPs influence vitamin B6 metabolite levels. Such an influence would suggest that vitamin B6 could play a role in GGE therapy. Here, we performed genome-wide association studies (GWAS) to assess the influence of GGE associated genetic variants on measures of vitamin B6 metabolism in blood plasma in 2232 healthy individuals. We also asked if SNPs that influence vitamin B6 were associated with GGE in 3122 affected individuals and 20,244 controls. Our GWAS of vitamin B6 metabolites reproduced a previous association and found a novel genome-wide significant locus. The SNPs in these loci were not associated with GGE. We found that 84 GGE-associated SNPs influence expression levels of PNPO in the brain as well as in blood. However, these SNPs were not associated with vitamin B6 metabolism in plasma. By leveraging polygenic risk scoring (PRS), we found suggestive evidence of higher catabolism and lower levels of the active and transport forms of vitamin B6 in GGE, although these findings require further replication.

Role of the Complement System in Chronic Central Serous Chorioretinopathy: A Genome-Wide Association Study

Importance

To date, several targeted genetic studies on chronic central serous chorioretinopathy (cCSC) have been performed; however, unbiased genome-wide studies into the genetics of cCSC have not been reported. To discover new genetic loci associated with cCSC and to better understand the causative mechanism of this disease, we performed a genome-wide association study (GWAS) on patients with cCSC.

Objective

To discover new genetic loci and pathways associated with cCSC and to predict the association of genetic variants with gene expression in patients with cCSC.

Design, Setting, and Participants

This case-control GWAS was completed in the general community, 3 referral university medical centers, and outpatient care on Europeans individuals with cCSC and population-based control participants. Genotype data was collected from May 2013 to August 2017, and data analysis occurred from August 2017 to November 2017.

Main Outcomes and Measures

Associations of single-nucleotide polymorphisms, haplotypes, genetic pathways, and predicted gene expression with cCSC.

Results

A total of 521 patients with cCSC (median age, 51 years; interquartile range [IQR], 44-59 years; 420 [80.6%] male) and 3577 European population-based control participants (median age, 52 years; IQR, 37-71 years; 1630 [45.6%] male) were included. One locus on chromosome 1 at the complement factor H (CFH) gene reached genome-wide significance and was associated with an increased risk of cCSC (rs1329428; odds ratio [OR], 1.57 [95% CI, 1.38-1.80]; P = 3.12 × 10-11). The CFH haplotypes H1 and H3 were protective for cCSC (H1: OR, 0.64 [95% CI, 0.53-0.77]; P = 2.18 × 10-6; H3: OR, 0.54 [95% CI, 0.42-0.70]; P = 2.49 × 10-6), whereas haplotypes H2, H4, H5, and the aggregate of rare CFH haplotypes conferred increased risk (H2: OR, 1.57 [95% CI, 1.30-1.89]; P = 2.18 × 10-6; H4: OR, 1.43 [95% CI, 1.13-1.80]; P = 2.49 × 10-3; H5: OR, 1.80 [95% CI, 1.36-2.39]; P = 4.61 × 10-5; rare haplotypes: OR, 1.99 [95% CI, 1.43-2.77]; P = 4.59 × 10-5). Pathway analyses showed involvement of the complement cascade and alternative open reading frame (ARF) pathway in cCSC. Using PrediXcan, we identified changes in predicted expression of complement genes CFH, complement factor H related 1 (CFHR1), complement factor related 4 (CFHR4), and membrane cofactor protein (MCP/CD46). Additionally, the potassium sodium-activated channel subfamily T member 2 (KCNT2) and tumor necrosis factor receptor superfamily member 10a (TNFRSF10A) genes were differentially expressed in patients with cCSC.

Conclusions and Relevance

In this GWAS on cCSC, we identified a locus on chromosome 1 at the CFH gene that was significantly associated with cCSC, and we report protective and risk-conferring haplotypes in this gene. Pathway analyses were enriched for complement genes, and gene expression analysis suggests a role for CFH, CFHR1, CFHR4, CD46, KCNT2, and TNFRSF10A in the disease. Taken together, these results underscore the potential importance of the complement pathway in the causative mechanisms of cCSC.

Comparative effectiveness of antiepileptic drugs in juvenile myoclonic epilepsy

OBJECTIVE

To study the effectiveness and tolerability of antiepileptic drugs (AEDs) commonly used in juvenile myoclonic epilepsy (JME).

METHODS

People with JME were identified from a large database of individuals with epilepsy, which includes detailed retrospective information on AED use. We assessed secular changes in AED use and calculated rates of response (12-month seizure freedom) and adverse drug reactions (ADRs) for the five most common AEDs. Retention was modeled with a Cox proportional hazards model. We compared valproate use between males and females.

RESULTS

We included 305 people with 688 AED trials of valproate, lamotrigine, levetiracetam, carbamazepine, and topiramate. Valproate and carbamazepine were most often prescribed as the first AED. The response rate to valproate was highest among the five AEDs (42.7%), and significantly higher than response rates for lamotrigine, carbamazepine, and topiramate; the difference to the response rate to levetiracetam (37.1%) was not significant. The rates of ADRs were highest for topiramate (45.5%) and valproate (37.5%). Commonest ADRs included weight change, lethargy, and tremor. In the Cox proportional hazards model, later start year (1.10 [1.08-1.13], P < 0.001) and female sex (1.41 [1.07-1.85], P = 0.02) were associated with shorter trial duration. Valproate was associated with the longest treatment duration; trials with carbamazepine and topiramate were significantly shorter (HR [CI]: 3.29 [2.15-5.02], P < 0.001 and 1.93 [1.31-2.86], P < 0.001). The relative frequency of valproate trials shows a decreasing trend since 2003 while there is an increasing trend for levetiracetam. Fewer females than males received valproate (76.2% vs 92.6%, P = 0.001).

SIGNIFICANCE

In people with JME, valproate is an effective AED; levetiracetam emerged as an alternative. Valproate is now contraindicated in women of childbearing potential without special precautions. With appropriate selection and safeguards in place, valproate should remain available as a therapy, including as an alternative for women of childbearing potential whose seizures are resistant to other treatments.

Influence of common SCN1A promoter variants on the severity of SCN1A-related phenotypes

Background

Pathogenic variants in SCN1A cause variable epilepsy disorders with different disease severities. We here investigate whether common variation in the promoter region of the unaffected SCN1A allele could reduce normal expression, leading to a decreased residual function of Nav1.1, and therefore to more severe clinical outcomes in patients affected by pathogenic SCN1A variants.

Methods

Five different SCN1A promoter‐haplotypes were functionally assessed in SH‐SY5Y cells using Firefly and Renilla luciferase assays. The SCN1A promoter region was analyzed in a cohort of 143 participants with SCN1A pathogenic variants. Differences in clinical features and outcomes between participants with and without common variants in the SCN1A promoter‐region of their unaffected allele were investigated.

Results

All non‐wildtype haplotypes showed a significant reduction in luciferase expression, compared to the wildtype promoter‐region (65%–80%, p = 0.039–0.0023). No statistically significant differences in clinical outcomes were observed between patients with and without common promoter variants. However, patients with a wildtype promoter‐haplotype on their unaffected SCN1A allele showed a nonsignificant trend for milder phenotypes.

Conclusion

The nonsignificant observed trends in our study warrant replication studies in larger cohorts to explore the potential modifying role of these common SCN1A promoter‐haplotypes.

Liprin alfa 2 gene expression is increased by cannabis use and associated with neuropsychological function

The relation of heavy cannabis use with decreased neuropsychological function has frequently been described but the underlying biological mechanisms are still largely unknown. This study investigates the relation of cannabis use with genome wide gene expression and subsequently examines the relations with neuropsychological function. Genome-wide gene expression in whole blood was compared between heavy cannabis users (N = 90) and cannabis naïve participants (N = 100) that were matched for psychotic like experiences. The results were validated using quantitative real-time PCR. Psychotic like experiences were assessed using the Comprehensive Assessment of Psychotic Experiences (CAPE). Neuropsychological function was estimated using four subtasks of the Wechsler Adult Intelligence Scale (WAIS). Subsequent in vitro studies in monocytes and a neuroblastoma cell line investigated expression changes in response to two major psychotropic components of cannabis; tetrahydrocannabinol (THC) and cannabidiol (CBD). mRNA expression of Protein Tyrosine Phosphatase Receptor Type F Polypeptide-Interacting-Protein Alpha-2 (PPFIA2) was significantly higher in cannabis users (LogFold Change 0.17) and confirmed by qPCR analysis. PPFIA2 expression level was negatively correlated with estimated intelligence (B=-22.9, p = 0.002) also in the 100 non-users (B=-28.5, p = 0.037). In vitro exposure of monocytes to CBD led to significant increase in PPFIA2 expression. However, exposure of monocytes to THC and neuroblastoma cells to THC or CBD did not change PPFIA2 expression. Change in PPFIA2 gene expression in response to cannabinoids is a putative mechanism by which cannabis could influence neuropsychological functions. The findings warrant further exploration of the role of PPFIA2 in cannabis induced changes of neuropsychological function, particularly in relation to CBD.

A genome-wide association study of sodium levels and drug metabolism in an epilepsy cohort treated with carbamazepine and oxcarbazepine

OBJECTIVE

To ascertain the clinical and genetic factors contributing to carbamazepine- and oxcarbazepine-induced hyponatremia (COIH), and to carbamazepine (CBZ) metabolism, in a retrospectively collected, cross-sectional cohort of people with epilepsy.

METHODS

We collected data on serum sodium levels and antiepileptic drug levels in people with epilepsy attending a tertiary epilepsy center while on treatment with CBZ or OXC. We defined hyponatremia as Na+ ≤134 mEq/L. We estimated the CBZ metabolic ratio defined as the log transformation of the ratio of metabolite CBZ-diol to unchanged drug precursor substrate as measured in serum.

RESULTS

Clinical and genetic data relating to carbamazepine and oxcarbazepine trials were collected in 1141 patients. We did not observe any genome-wide significant associations with sodium level in a linear trend or hyponatremia as a dichotomous trait. Age, sex, number of comedications, phenytoin use, phenobarbital use, and sodium valproate use were significant predictors of CBZ metabolic ratio. No genome-wide significant associations with CBZ metabolic ratio were found.

SIGNIFICANCE

Although we did not detect a genetic predictor of hyponatremia or CBZ metabolism in our cohort, our findings suggest that the determinants of CBZ metabolism are multifactorial.

Outcomes and comorbidities of SCN1A-related seizure disorders

PURPOSE

Differentiating between Dravet syndrome and non-Dravet SCN1A-related phenotypes is important for prognosis regarding epilepsy severity, cognitive development, and comorbidities. When a child is diagnosed with genetic epilepsy with febrile seizures plus (GEFS+) or febrile seizures (FS), accurate prognostic information is essential as well, but detailed information on seizure course, seizure freedom, medication use, and comorbidities is lacking for this milder patient group. In this cross-sectional study, we explore disease characteristics in milder SCN1A-related phenotypes and the nature, occurrence, and relationships of SCN1A-related comorbidities in both patients with Dravet and non-Dravet syndromes.

METHODS

A cohort of 164 Dutch participants with SCN1A-related seizures was evaluated, consisting of 116 patients with Dravet syndrome and 48 patients with either GEFS+, febrile seizures plus (FS+), or FS. Clinical data were collected from medical records, semi-structured telephone interviews, and three questionnaires: the Functional Mobility Scale (FMS), the Pediatric Quality of Life Inventory (PedsQL) Measurement Model, and the Child or Adult Behavior Checklists (CBCL/ABCL).

RESULTS

Walking disabilities and severe behavioral problems affect 71% and 43% of patients with Dravet syndrome respectively and are almost never present in patients with non-Dravet syndromes. These comorbidities are strongly correlated to lower quality-of-life (QoL) scores. Less severe comorbidities occur in patients with non-Dravet syndromes: learning problems and psychological/behavioral problems are reported for 27% and 38% respectively. The average QoL score of the non-Dravet group was comparable with that of the general population. The majority of patients with non-Dravet syndromes becomes seizure-free after 10 years of age (85%).

CONCLUSIONS

Severe behavioral problems and walking disabilities are common in patients with Dravet syndrome and should receive specific attention during clinical management. Although the epilepsy course of patients with non-Dravet syndromes is much more favorable, milder comorbidities frequently occur in this group as well. Our results may be of great value for clinical care and informing newly diagnosed patients and their parents about prognosis.

Functionally distinct ERAP1 and ERAP2 are a hallmark of HLA-A29-(Birdshot) Uveitis

Birdshot Uveitis (Birdshot) is a rare eye condition that affects HLA-A29-positive individuals and could be considered a prototypic member of the recently proposed 'MHC-I (major histocompatibility complex class I)-opathy' family. Genetic studies have pinpointed the endoplasmic reticulum aminopeptidase (ERAP1) and (ERAP2) genes as shared associations across MHC-I-opathies, which suggests ERAP dysfunction may be a root cause for MHC-I-opathies. We mapped the ERAP1 and ERAP2 haplotypes in 84 Dutch cases and 890 controls. We identified association at variant rs10044354, which mediated a marked increase in ERAP2 expression. We also identified and cloned an independently associated ERAP1 haplotype (tagged by rs2287987) present in more than half of the cases; this ERAP1 haplotype is also the primary risk and protective haplotype for other MHC-I-opathies. We show that the risk ERAP1 haplotype conferred significantly altered expression of ERAP1 isoforms in transcriptomic data (n = 360), resulting in lowered protein expression and distinct enzymatic activity. Both the association for rs10044354 (meta-analysis: odds ratio (OR) [95% CI]=2.07[1.58-2.71], P = 1.24 × 10(-7)) and rs2287987 (OR[95% CI]: =2.01[1.51-2.67], P = 1.41 × 10(-6)) replicated and showed consistent direction of effect in an independent Spanish cohort of 46 cases and 2103 controls. In both cohorts, the combined rs2287987-rs10044354 haplotype associated with Birdshot more strongly than either variant alone [meta-analysis: P=3.9 × 10(-9)]. Finally, we observed that ERAP2 protein expression is dependent on the ERAP1 background across three European populations (n = 3353). In conclusion, a functionally distinct combination of ERAP1 and ERAP2 are a hallmark of Birdshot and provide rationale for strategies designed to correct ERAP function for treatment of Birdshot and MHC-I-opathies more broadly.

Assessment of parental mosaicism in SCN1A-related epilepsy by single-molecule molecular inversion probes and next-generation sequencing

Background

Dravet syndrome is a severe genetic encephalopathy, caused by pathogenic variants in SCN1A. Low-grade parental mosaicism occurs in a substantial proportion of families (7%–13%) and has important implications for recurrence risks. However, parental mosaicism can remain undetected by methods regularly used in diagnostics. In this study, we use single-molecule molecular inversion probes (smMIP), a technique with high sensitivity for detecting low-grade mosaic variants and high cost-effectiveness, to investigate the incidence of parental mosaicism of SCN1A variants in a cohort of 90 families and assess the feasibility of this technique.

Methods

Deep sequencing of SCN1A was performed using smMIPs. False positive rates for each of the proband’s pathogenic variants were determined in 145 unrelated samples. If parents showed corresponding variant alleles at a significantly higher rate than the established noise ratio, mosaicism was confirmed by droplet digital PCR (ddPCR).

Results

Sequence coverage of at least 100× at the location of the corresponding pathogenic variant was reached for 80 parent couples. The variant ratio was significantly higher than the established noise ratio in eight parent couples, of which four (5%) were regarded as true mosaics, based on ddPCR results. The false positive rate of smMIP analysis without ddPCR was therefore 50%. Three of these variants had previously been considered de novo in the proband by Sanger sequencing.

Conclusion

smMIP technology combined withnext generation sequencing (NGS) performs better than Sanger sequencing in the detection of parental mosaicism. Because parental mosaicism has important implications for genetic counselling and recurrence risks, we stress the importance of implementing high-sensitivity NGS-based assays in standard diagnostics.

Refractory juvenile myoclonic epilepsy: a meta-analysis of prevalence and risk factors

Background and purpose

Juvenile myoclonic epilepsy (JME) is a common epilepsy syndrome for which treatment response is generally assumed to be good. We aimed to determine the prevalence and prognostic risk factors for refractoriness of JME.

Methods

We systematically searched PubMed and EMBASE and included 43 eligible studies, reporting seizure outcome after antiepileptic drug (AED) treatment in JME cohorts. We defined refractory JME as persistence of any seizure despite AED treatment and performed a random‐effects meta‐analysis to assess the prevalence of refractory JME and of seizure recurrence after AED withdrawal in individuals with well‐controlled seizures. Studies reporting potential prognostic risk factors in relation to seizure outcome were included for subsequent meta‐analysis of risk factors for refractoriness.

Results

Overall, 35% (95% confidence interval, 29–41%) of individuals (n = 3311) were refractory. There was marked heterogeneity between studies. Seizures recurred in 78% (95% confidence interval, 52–94%) of individuals who attempted to withdraw from treatment after a period of seizure freedom (n = 246). Seizure outcome by publication year suggested that prognosis did not improve over time. Meta‐analysis suggested six variables as prognostic factors for refractoriness, i.e. having three seizure types, absence seizures, psychiatric comorbidities, earlier age at seizure onset, history of childhood absence epilepsy and praxis‐induced seizures.

Conclusion

One‐third of people with JME were refractory, which is a higher prevalence than expected. Risk factors were identified and can be used to guide treatment and counselling of people with JME.

Location on the Human Genetic Linkage Map of 26 Genes Involved in Blood Coagulation

Several human genetic linkage maps have been constructed as part of the Human Genome Project. These maps show the positional order of closely linked, highly informative AC-repeat polymorphisms on each human chromosome, and are extremely useful in genetic linkage analysis of inheritable diseases. For a candidate gene approach the current linkage maps are less useful, since they consist mainly of anonymous markers rather than of specific genes. This situation also applies for inheritable disorders of blood coagulation. Numerous genes are involved in the blood coagulation cascade and its regulation, and can be considered as candidate genes for unexplained haemophilia and thrombophilia. We have selected 29 candidate genes that seem to be the ones most likely to be involved in thrombophilia. For 19 genes genotype data were already present in the CEPH database (version 7.0). We typed 7 additional genes in the CEPH reference families, i.e. the factor V, factor XII, protein C, protein S, prothrombin, thrombomodulin, and heparin cofactor II gene. The genotype data were used to integrate these 26 genes in the current genetic linkage map, and to identify closely linked AC-repeat polymorphisms. This information will benefit the investigation of inheritable disorders of blood coagulation, especially thrombophilia.

Factor V Leiden: An Additional Risk Factor for Thrombosis in Protein S Deficient Families?

We recently reported a high prevalence of the FV Leiden mutation (R506Q, responsible for Activated Protein C resistance) among symptomatic protein C deficient probands (19%), and the involvement of the FV Leiden mutation in the expression of thrombophilia in six protein C deficient families. Here, we report the results of a similar study in protein S deficient probands and families. Among 16 symptomatic protein S deficient probands the prevalence of the FV Leiden mutation was high (38%). This high prevalence is significantly different from that in the normal population, and is probably caused by the selection of probands for familial thrombosis and protein S deficiency. In 4 families, the segregation of the FV Leiden mutation and the protein S deficiency could be studied. In sibships where both abnormalities were segregating, the percentage of symptomatic individuals with both abnormalities was 80%. Three of the seven subjects with only the FV Leiden mutation, and two out of the three subjects with only protein S deficiency had developed thrombosis. These results indicate that in the families presented here the combination of the FV Leiden mutation and the protein S deficiency is associated with a high risk for thrombosis. A reliable estimate of the penetrance of the single defects is not possible, because the number of individuals with a single defect is too low.

Abnormal islet sphingolipid metabolism in type 1 diabetes

AIMS/HYPOTHESIS

Sphingolipids play important roles in beta cell physiology, by regulating proinsulin folding and insulin secretion and in controlling apoptosis, as studied in animal models and cell cultures. Here we investigate whether sphingolipid metabolism may contribute to the pathogenesis of human type 1 diabetes and whether increasing the levels of the sphingolipid sulfatide would prevent models of diabetes in NOD mice.

METHODS

We examined the amount and distribution of sulfatide in human pancreatic islets by immunohistochemistry, immunofluorescence and electron microscopy. Transcriptional analysis was used to evaluate expression of sphingolipid-related genes in isolated human islets. Genome-wide association studies (GWAS) and a T cell proliferation assay were used to identify type 1 diabetes related polymorphisms and test how these affect cellular islet autoimmunity. Finally, we treated NOD mice with fenofibrate, a known activator of sulfatide biosynthesis, to evaluate the effect on experimental autoimmune diabetes development.

RESULTS

We found reduced amounts of sulfatide, 23% of the levels in control participants, in pancreatic islets of individuals with newly diagnosed type 1 diabetes, which were associated with reduced expression of enzymes involved in sphingolipid metabolism. Next, we discovered eight gene polymorphisms (ORMDL3, SPHK2, B4GALNT1, SLC1A5, GALC, PPARD, PPARG and B4GALT1) involved in sphingolipid metabolism that contribute to the genetic predisposition to type 1 diabetes. These gene polymorphisms correlated with the degree of cellular islet autoimmunity in a cohort of individuals with type 1 diabetes. Finally, using fenofibrate, which activates sulfatide biosynthesis, we completely prevented diabetes in NOD mice and even reversed the disease in half of otherwise diabetic animals.

CONCLUSIONS/INTERPRETATION

These results indicate that islet sphingolipid metabolism is abnormal in type 1 diabetes and suggest that modulation may represent a novel therapeutic approach.

DATA AVAILABILITY

The RNA expression data is available online at https://www.dropbox.com/s/93mk5tzl5fdyo6b/Abnormal%20islet%20sphingolipid%20metabolism%20in%20type%201%20diabetes%2C%20RNA%20expression.xlsx?dl=0 . A list of SNPs identified is available at https://www.dropbox.com/s/yfojma9xanpp2ju/Abnormal%20islet%20sphingolipid%20metabolism%20in%20type%201%20diabetes%20SNP.xlsx?dl=0 .

Genetic variation in CFH predicts phenytoin-induced maculopapular exanthema in European-descent patients

OBJECTIVE

To characterize, among European and Han Chinese populations, the genetic predictors of maculopapular exanthema (MPE), a cutaneous adverse drug reaction common to antiepileptic drugs.

METHODS

We conducted a case-control genome-wide association study of autosomal genotypes, including Class I and II human leukocyte antigen (HLA) alleles, in 323 cases and 1,321 drug-tolerant controls from epilepsy cohorts of northern European and Han Chinese descent. Results from each cohort were meta-analyzed.

RESULTS

We report an association between a rare variant in the complement factor H-related 4 (CFHR4) gene and phenytoin-induced MPE in Europeans (p = 4.5 × 10-11; odds ratio [95% confidence interval] 7 [3.2-16]). This variant is in complete linkage disequilibrium with a missense variant (N1050Y) in the complement factor H (CFH) gene. In addition, our results reinforce the association between HLA-A*31:01 and carbamazepine hypersensitivity. We did not identify significant genetic associations with MPE among Han Chinese patients.

CONCLUSIONS

The identification of genetic predictors of MPE in CFHR4 and CFH, members of the complement factor H-related protein family, suggest a new link between regulation of the complement system alternative pathway and phenytoin-induced hypersensitivity in European-ancestral patients.

Neurodevelopmental Disorders Caused by De Novo Variants in KCNB1 Genotypes and Phenotypes

Importance

Knowing the range of symptoms seen in patients with a missense or loss-of-function variant in KCNB1 and how these symptoms correlate with the type of variant will help clinicians with diagnosis and prognosis when treating new patients.

Objectives

To investigate the clinical spectrum associated with KCNB1 variants and the genotype-phenotype correlations.

Design, Setting, and Participants

This study summarized the clinical and genetic information of patients with a presumed pathogenic variant in KCNB1. Patients were identified in research projects or during clinical testing. Information on patients from previously published articles was collected and authors contacted if feasible. All patients were seen at a clinic at one of the participating institutes because of presumed genetic disorder. They were tested in a clinical setting or included in a research project.

Main Outcomes and Measures

The genetic variant and its inheritance and information on the patient's symptoms and characteristics in a predefined format. All variants were identified with massive parallel sequencing and confirmed with Sanger sequencing in the patient. Absence of the variant in the parents could be confirmed with Sanger sequencing in all families except one.

Results

Of 26 patients (10 female, 15 male, 1 unknown; mean age at inclusion, 9.8 years; age range, 2-32 years) with developmental delay, 20 (77%) carried a missense variant in the ion channel domain of KCNB1, with a concentration of variants in region S5 to S6. Three variants that led to premature stops were located in the C-terminal and 3 in the ion channel domain. Twenty-one of 25 patients (84%) had seizures, with 9 patients (36%) starting with epileptic spasms between 3 and 18 months of age. All patients had developmental delay, with 17 (65%) experiencing severe developmental delay; 14 (82%) with severe delay had behavioral problems. The developmental delay was milder in 4 of 6 patients with stop variants and in a patient with a variant in the S2 transmembrane element rather than the S4 to S6 region.

Conclusions and Relevance

De novo KCNB1 missense variants in the ion channel domain and loss-of-function variants in this domain and the C-terminal likely cause neurodevelopmental disorders with or without seizures. Patients with presumed pathogenic variants in KCNB1 have a variable phenotype. However, the type and position of the variants in the protein are (imperfectly) correlated with the severity of the disorder.

Brief Report: IRF4 Newly Identified as a Common Susceptibility Locus for Systemic Sclerosis and Rheumatoid Arthritis in a Cross-Disease Meta-Analysis of Genome-Wide Association Studies

OBJECTIVE

Systemic sclerosis (SSc) and rheumatoid arthritis (RA) are autoimmune diseases that have similar clinical and immunologic characteristics. To date, several shared SSc-RA genetic loci have been identified independently. The aim of the current study was to systematically search for new common SSc-RA loci through an interdisease meta-genome-wide association (meta-GWAS) strategy.

METHODS

The study was designed as a meta-analysis combining GWAS data sets of patients with SSc and patients with RA, using a strategy that allowed identification of loci with both same-direction and opposite-direction allelic effects. The top single-nucleotide polymorphisms were followed up in independent SSc and RA case-control cohorts. This allowed an increase in the sample size to a total of 8,830 patients with SSc, 16,870 patients with RA, and 43,393 healthy controls.

RESULTS

This cross-disease meta-analysis of the GWAS data sets identified several loci with nominal association signals (P < 5 × 10(-6) ) that also showed evidence of association in the disease-specific GWAS scans. These loci included several genomic regions not previously reported as shared loci, as well as several risk factors that were previously found to be associated with both diseases. Follow-up analyses of the putatively new SSc-RA loci identified IRF4 as a shared risk factor for these 2 diseases (Pcombined  = 3.29 × 10(-12) ). Analysis of the biologic relevance of the known SSc-RA shared loci identified the type I interferon and interleukin-12 signaling pathways as the main common etiologic factors.

CONCLUSION

This study identified a novel shared locus, IRF4, for the risk of SSc and RA, and highlighted the usefulness of a cross-disease GWAS meta-analysis strategy in the identification of common risk loci.

Carbamazepine- and oxcarbazepine-induced hyponatremia in people with epilepsy

OBJECTIVE

To ascertain possible determinants of carbamazepine (CBZ)- and oxcarbazepine (OXC)-induced hyponatremia in a large cohort of people with epilepsy.

METHODS

We collected data on serum sodium levels in people with epilepsy who were attending a tertiary epilepsy center while on treatment with CBZ or OXC. We defined hyponatremia as Na+ ≤134 mEq/L and severe hyponatremia as Na+ ≤128 mEq/L.

RESULTS

We identified 1,782 people who had used CBZ (n = 1,424) or OXC (n = 358), of whom 50 were treated with both drugs. Data on sodium level measurements were available in 1,132 on CBZ and in 289 on OXC. Hyponatremia occurred in 26% of those taking CBZ and 46% of those taking OXC. This was severe in 7% in the CBZ group and 22% in the OXC group. Hyponatremia was symptomatic in 48% and led to admissions in 3%. Age over 40 years, high serum levels of CBZ and OXC, and concomitant use of other antiepileptic drugs were the main risk factors for hyponatremia in both treatment groups. Female patients on OXC were at a higher risk than male patients of hyponatremia. The risk of hyponatremia on CBZ was significantly associated with the risk of hyponatremia on OXC within a subgroup that used both drugs consecutively.

SIGNIFICANCE

Hyponatremia is a common problem in people taking CBZ or OXC. Regular ascertainment of sodium levels in those taking either drug is recommended and results should be acted on.

Identification of coagulation gene 3'UTR variants that are potentially regulated by microRNAs

MicroRNAs have been recognized as critical regulators of gene expression and might affect the risk of venous thrombosis. We aimed to identify 3' untranslated region (UTR) variants in coagulation genes that influence coagulation factor levels and venous thrombosis risk. The 3'UTR of coagulation genes were sequenced in subjects with extremely high or low plasma levels of these factors in two case-control studies. In total, 28 variants were identified. Five single nucleotide polymorphisms (SNPs) were predominantly present in one extreme level group (F2 rs1799963, F8 rs1050705 and F11 rs4253429, rs4253430 and rs1062547). Additional to F2 rs1799963, F8 rs1050705 (in men) and F11 rs4253430 were associated with an increased risk of venous thrombosis albeit confidence intervals were wide. The three F11 SNPs were in high linkage disequilibrium with functional variants rs2289252 and rs2036914. Rs1062547 and rs4253430 were associated with a significant increase of plasma FXI activity in heterozygotes and homozygotes in wild-type controls. In silico prediction revealed that these SNPs might disturb the binding sites of miR-544 and miR-513a-3p. Only miR-544 provoked a significant decrease of the luciferase activity that was not observed with a rs4253430 mutated vector. In conclusion, these results reinforce that microRNAs are candidates to play a role in haemostasis and complex disorders, such as thrombosis.

Genome-wide association study of sporadic brain arteriovenous malformations

BACKGROUND

The pathogenesis of sporadic brain arteriovenous malformations (BAVMs) remains unknown, but studies suggest a genetic component. We estimated the heritability of sporadic BAVM and performed a genome-wide association study (GWAS) to investigate association of common single nucleotide polymorphisms (SNPs) with risk of sporadic BAVM in the international, multicentre Genetics of Arteriovenous Malformation (GEN-AVM) consortium.

METHODS

The Caucasian discovery cohort included 515 BAVM cases and 1191 controls genotyped using Affymetrix genome-wide SNP arrays. Genotype data were imputed to 1000 Genomes Project data, and well-imputed SNPs (>0.01 minor allele frequency) were analysed for association with BAVM. 57 top BAVM-associated SNPs (51 SNPs with p<10(-05) or p<10(-04) in candidate pathway genes, and 6 candidate BAVM SNPs) were tested in a replication cohort including 608 BAVM cases and 744 controls.

RESULTS

The estimated heritability of BAVM was 17.6% (SE 8.9%, age and sex-adjusted p=0.015). None of the SNPs were significantly associated with BAVM in the replication cohort after correction for multiple testing. 6 SNPs had a nominal p<0.1 in the replication cohort and map to introns in EGFEM1P, SP4 and CDKAL1 or near JAG1 and BNC2. Of the 6 candidate SNPs, 2 in ACVRL1 and MMP3 had a nominal p<0.05 in the replication cohort.

CONCLUSIONS

We performed the first GWAS of sporadic BAVM in the largest BAVM cohort assembled to date. No GWAS SNPs were replicated, suggesting that common SNPs do not contribute strongly to BAVM susceptibility. However, heritability estimates suggest a modest but significant genetic contribution.