# Rare variants in the GABAA receptor subunit epsilon identified in patients with a wide spectrum of epileptic phenotypes

Fenja Markus, Chloé Angelini, Aurelien Trimouille, Gabrielle Rudolf, Gaetan

Lesca, Cyril Goizet, Eulalie Lasseaux, Benoit Arveiler, Marjon van

Slegtenhorst, Alice S. Brooks, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Fenja Markus, Chloé Angelini, Aurelien Trimouille, Gabrielle Rudolf, Gaetan Lesca, et al. Rare variants in the GABAA receptor subunit epsilon identified in patients with a wide spectrum of epileptic phenotypes. Molecular Genetics and Genomic Medicine, 2020, 8 (9), 10.1002/mgg3.1388. hal-03700156

HAL Id: hal-03700156

https://hal.science/hal-03700156

Submitted on 9 Feb 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# Rare variants in the GABA<sub>A</sub> receptor subunit ε identified in patients with a wide spectrum of epileptic phenotypes

Fenja Markus<sup>1,2</sup> | Chloé Angelini<sup>3</sup> | Aurelien Trimouille<sup>3 ID</sup> | Gabrielle Rudolf<sup>4,5</sup> |

Gaetan Lesca<sup>6</sup> | Cyril Goizet<sup>3,4</sup> | Eulalie Lasseaux<sup>3</sup> | Benoit Arveiler<sup>3</sup> |

Marjon van Slegtenhorst<sup>7</sup> | Alice S. Brooks<sup>7</sup> | Rami Abou Jamra<sup>8</sup> |

Georg-Christoph Korenke<sup>9</sup> | John Neidhardt<sup>2,10</sup> | Marta Owczarek-Lipska<sup>1,2 ID</sup>

## Correspondence

Marta Owczarek-Lipska, Junior Research Group, Genetics of childhood brain malformations, Faculty VI-School of Medicine and Health Sciences, University of Oldenburg, Ammerländer Heerstr. 114-118, 26129 Oldenburg, Germany.

Email: marta.owczarek-lipska@uni-oldenburg.de

## Abstract

**Background:** Epilepsy belongs to a group of chronic and highly heterogeneous brain disorders. Many types of epilepsy and epileptic syndromes are caused by genetic factors.

The neural amino acid y-aminobutyric acid (GABA) is a major inhibitory neurotransmitter in the mammalian central nervous system. It regulates activity of channel pores by binding to transmembrane GABA-receptors (GABRs). The GABRs are heteropeptamers assembled from different receptor subunits (α1-6, β1-3, γ1-3, δ, ε, θ, π, and ρ1-3). Several epileptic disorders are caused by mutations in genes encoding single GABRs.

**Methods:** We applied trio- and single-whole exome sequencing to search for genetic sequence variants associated with a wide range of epileptic phenotypes accompanied by intellectual disability and/or global developmental delay in the investigated patients.

**Results:** We identified four hemizygous sequence variants in the GABA<sub>A</sub> receptor subunit ε gene (GABRE), including one nonsense (NM_004961.3: c.399C>A, p.Tyr133*),

<footnote n="1"><sup>1</sup>Junior Research Group, Genetics of Childhood Brain Malformations, Faculty VI-School of Medicine and Health Sciences, University of Oldenburg, Oldenburg, Germany</footnote>

<footnote n="2"><sup>2</sup>Human Genetics, Faculty VI-School of Medicine and Health Sciences, University of Oldenburg, Oldenburg, Germany</footnote>

<footnote n="3"><sup>3</sup>Service de Génétique médicale, CHU de Bordeaux, Bordeaux, France</footnote>

<footnote n="4"><sup>4</sup>CNRS U7104, INSERM U1258, Institut de Génétique et de Biologie Moléculaire et Cellulaire, Illkirch, France</footnote>

<footnote n="5"><sup>5</sup>Service de Neurologie, Centre de Références des Maladies Neurogénétique Rares, Hôpitaux Universitaires de Strasbourg, Strasbourg, France</footnote>

<footnote n="6"><sup>6</sup>Genetics department, Lyon University Hospital and University of Lyon, Lyon, France</footnote>

<footnote n="7"><sup>7</sup>Department of Clinical Genetics, Erasmus MC University Medical Center, Rotterdam, The Netherlands</footnote>

<footnote n="8"><sup>8</sup>Institute of Human Genetics, University Medical Center Leipzig, Leipzig, Germany</footnote>

<footnote n="9"><sup>9</sup>Department of Neuropediatrics, University Children's Hospital, Oldenburg, Germany</footnote>

<footnote n="10"><sup>10</sup>Research Center Neurosensory Science, University of Oldenburg, Oldenburg, Germany</footnote>
two missense variants (NM_004961.3: c.664G>A, p.Glu222Lys; NM_004961.3: c.1045G>A, p.Val349Ile), and one variant affecting the translation initiation codon (NM_004961.3: c.1A>G, p.Met1?) in four unrelated families.

**Conclusion:** Our clinical and molecular genetic findings suggest that *GABRE* is a likely candidate gene for epilepsy. Nevertheless, functional studies are necessary to better understand pathogenicity of the *GABRE*-mutations and their associations with epileptic phenotypes.

### KEYWORDS

exome sequencing, GABA receptor type A subunit epsilon, GABRE, genetic epilepsy, novel sequence variants

# 1 | INTRODUCTION

The International League Against Epilepsy (ILAE) defined epilepsy as a prevalent and highly heterogeneous disorder of the brain (Fisher et al., 2014). This disorder leads to seizures, defining sudden erroneous and unprovoked electrical conduction between neurons (Guerreiro, 2016). The ILAE announced a three-level classification system for epileptic seizure types, epilepsy types, and epileptic syndromes that support a better clinical diagnosis and treatment of patients affected with epilepsy (Chang, Leung, Ho, & Yung, 2017). Among epileptic seizure types three major groups can be distinguished: (a) focal onset seizures derived from neuronal networks which are exclusively restricted to an unilateral brain hemisphere, (b) generalized onset seizures which originate from bilateral neuronal networks, and (c) idiopathic epilepsies without any evident etiology, frequently showing genetic predispositions (Berg et al., 2018; Chang et al., 2017; Guerreiro, 2016; Hirose, Mitsudome, Okada, & Kaneko, 2005). Additionally, epileptic disorders were classified into genetic, structural, metabolic, and unknown origin (Berg et al., 2010; Korff & Scheffer, 2013).

Genetic factors play an important role in many types of epilepsy and epileptic syndromes, often in these with unknown origin (Nieh & Sherr, 2014). Despite a fast development in bio-medical technologies and advances in next generation sequencing (NGS), resulting in a substantial progress in the discovery of genes and/or multiple-loci associated with epileptic disorders, many of the genotype-phenotype correlations remain to be elucidated (Symonds, Zuberi, & Johnson, 2017). It has been previously reported that certain forms of epilepsy can be caused by pathogenic variants in several distinct gamma-aminobutyric acid receptor (GABR) subunit genes (Butler et al., 2018; Hernandez et al., 2016b; Ishii et al., 2017; Shen et al., 2017). GABA is the major inhibitory neurotransmitter in the mammalian central nervous system (CNS), involved in numerous developmental processes in the brain (Macdonald

& Olsen, 1994; Watanabe, Maemura, Kanbara, Tamayama, & Hayasaki, 2002). It is involved in neural migration, facilitates extension of neurites, and plays an important role in a synaptic maintenance (Behar et al., 1994; Behar, Schaffner, Scott, O'Connell, & Barker, 1998; Watanabe et al., 2002). About 20% of neurons in the CNS are GABAergic (Watanabe et al., 2002). GABRs are pentameric, ligand-gated chloride channels, composed of a combination of different subunits (α1-6, β1-3, γ1-3, δ, ε, θ, π, and ρ1-3; Hernandez et al., 2016b; Macdonald & Olsen, 1994; Watanabe et al., 2002). Each GABR subunit is composed of a large extracellular N-terminal domain, four hydrophobic transmembrane domains, and a large intracellular domain, which is the most heterogeneous part of a single receptor (Moss & Smart, 2001). Ionotropic GABRs are classified based on their pharmacology into two types of receptors: GABA<sub>A</sub> and GABA<sub>C</sub> (Moss & Smart, 2001). GABA<sub>A</sub> receptors are important targets of drugs, including benzodiazepines, barbiturates, and other neuroactive medications (Sieghart, 2000). Pathogenic variants in genes encoding for several GABA<sub>A</sub> receptor subunits are known to cause epilepsies of different severity, from mild febrile seizures to severe epileptic encephalopathies (Carvill et al., 2014; Cossette et al., 2002; Hernandez et al., 2016b; Kodera et al., 2016; Moller et al., 2017; Shen et al., 2017). The majority of monogenic pathogenic variants are associated with *GABRA1*, *GABRB2/3*, and *GABRG2* genes, which form the abundant α1β2/3γ2 subtype of GABA<sub>A</sub> receptors in the CNS (Farrant & Nusser, 2005). Unlike many GABA<sub>A</sub> receptors, those containing the ε subunit have restricted distributions, form functional channels with distinct pharmaco-biochemical activities and lead to spontaneous channel opening with slowed recovery (Neelands, Fisher, Bianchi, & Macdonald, 1999; Wagner, Goldschen-Ohm, Hales, & Jones, 2005). The assembly of the ε subunit into GABA<sub>A</sub> receptors and its complete biological function is poorly understood (Bollan, Baur, Hales, Sigel, & Connolly, 2008). So far, only a few variants of unknown
Burrows-Wheeler Aligner (http://bio-bwa.sourceforge.net/). Variants were detected with Genome Analysis Toolkit (http://www.broadinstitute.org/gatk/). Subsequently, sequence variants were filtered with Cartagenia software package (Agilent Technologies), considering reads' quality (read depth ≥10), MAF ≤ 1.5% (dbSNP, ESP6500, the 1,000 Genome project, GoNL, or the ExAc database), and location (within an exon and 10 bp areas encompassing introns). Variants were further selected in regard to an autosomal recessive, X-linked or de novo dominant inheritance.

Trio-WES analyses of family 4 using DNA samples from P4 and healthy parents were performed on a HiSeq 1000 Sequencing System (Illumina). The sequencing was done with a 2 × 150 bp high output sequencing kit after enrichment with SeqCap EZ MedExome kit (Roche), according to manufacturer's specifications. Sequence quality was assessed with FastQC 0.11.5 (Babraham Bioinformatics). After sequencing, reads were mapped to the hg19/GRCh37, using BWA-MEM (http://bio-bwa.sourceforge.net), sorted and indexed (Samtools, http://samtools.sourceforge.net). Duplicates were flagged (Sambamba, https://www.open-bio.org/wiki/Sambamba) and coverage was calculated (Picard-Tools, https://broadinsti tute.github.io/picard/). Variant calling was performed with GATK 3.7 Haplotype Caller (https://gatk.broadinstitute.org/ hc/en-us) and variants were annotated with SnpEff (http:// snpeff.sourceforge.net/), dbNSFP (http://sites.google.com/ site/jpopgen/dbNSFP), GnomAD, ClinVar (https://www. ncbi.nlm.nih.gov/clinvar/), Human Gene Mutation Database (HGMD), Variome Great Middle East (http://igm.ucsd.edu/ gme/data-browser.php) Mutation Tester (http://www.mutat

iontaster.org/), Combined Annotation Dependent Depletion (CADD, https://cadd.gs.washington.edu/), as well as with an in-house-database.

## 2.4 | Primers, PCR amplification and Sanger sequencing

Sanger sequencing was performed in family 1. Primers (fwd_5'-TAAGAGAGGCAAGGTCCCAC, rev_5'-CTTG-GCAAAGCACTTACTTCTT) located in intronic regions of the human *GABRE* gene (NM_004961.3) and encompassing exon 6 were designed using Primer Input3 (http://prime r3.ut.ee/). They were additionally verified for common SNPs using SNPCheck (https://secure.ngrl.org.uk/SNPCheck/). In total, 10 ng of gDNA samples from P1 and from unaffected relatives (including parents, sister, and grandfather from maternal side) were used for PCR amplification with HotFirePol DNA Polymerase (Solis BioDyne), according to manufacturer's recommendations (Figure 1a; Figure S1a). Afterward, the amplicons were purified using enzymatic ExoI-SAP method (New England Biolabs) and bilaterally sequenced using the BigDye® Terminator v3.1 Cycle Sequencing Kit on the ABI Prism 3130xl Genetic Analyzer (Applied Biosystem). Sanger sequencing data were analyzed with SeqScape software (Applied Biosystem) and SnapGene software (GSL Biotech).

Results of high throughput sequencing from P2, P3, and P4 were not additionally verified by Sanger sequencing as the *GABRE* gene, especially the nucleotide positions of the

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** Analysis of identified hemizygous variants in *GABRE* in unrelated patients. (A) Family P1 includes healthy parents (I.I, I.II), one female healthy child (II.I), and the GE-affected child (P1, II.II). The father (I.I) showed the hemizygous reference sequence, while the mother (I.II) and the sister (II.I) were carriers of the missense variant NM_004961.3: c.664G>A. (B) Family P2 consists of healthy parents (father I.I, mother I.II), one male healthy child (II.II) and one affected child (P2, II.I, NM_004961.3: c.1A>G). The other members of this family are healthy. (C) Family P3 includes the healthy parents (father I.I, mother I.II) and two children, the healthy daughter (II.I) and P3 (II.II), who showed a C-to-A substitution, leading to a premature termination codon (NM_004961.3: c.399C>A). The other members of family P3 were healthy. (D) Family P4 consists of healthy parents (father I.II, mother I.II) and one affected child (P4, II.I, NM_004961.3: c.1045G>A). No relevant familial history was noted. All variants were detected using WES and reads were obtained from both directions, forward and reverse. Symbols: squares = males; circles = females; black solid squares with arrows = patients affected with epilepsy; solid white icons = healthy individuals; icons with dot = known carrier</figcaption></figure>
identified pathogenic sequence variants were very well cov-
ered in the trio-WES data.

**2.5 | Multiple species protein alignments and in silico predictions of GABRE protein domains**

Multiple sequence alignments were performed with the ClustalW2 tool (http://www.ebi.ac.uk/Tools/msa/clust alw2/). Amino acid sequences of GABRE proteins were compared between human (GenBank: NP_004952.2), mouse (GenBank: NP_059065.2), Norway rat (GenBank: NP_075579.1), rabbit (GenBank: XP_002721372.1), hedgehog tenrec (GenBank: XP_004715681.1), Sunda flying lemur (GenBank: XP_008589892.1), dog (GenBank: XP_549340.3), Californian Sea Lion (GenBank: XP_027465429.1), Chinese rufous horseshoe bat (GenBank: XP_019566807.1), and lesser ferret (GenBank: XP_004776709.1).

Publicly available tools to predict protein domains, such as NCBI Conserved Domain Search (https://www.ncbi.nlm.nih. gov/Structure/cdd/wrpsb.cgi), PROSITE (http://prosite.ex- pasy.org), and SMART (http://smart.embl-heidelberg.de) were used to in silico estimate locations of the mutated amino acids in the GABRE protein, found in our patients.

Additionally, MetaDome data (https://github.com/cmbi/
metadome) were used to predict pathogenicity and improve
interpretation of the newly identified GABRE-variants.
MetaDome aligns information of population variations from
the ExAc and pathogenic variants from the HGMD with Pfam
protein domains. Generally, MetaDome visualizes meta-
domain information and gene-wide profiles of genetic tol-
erance based on 56'319 human transcripts, 71'419 protein
domains, 12'164'292 genetic variants from GnomAD, and
34'076 pathogenic variants from ClinVar (Wiel et al., 2019).

3 | RESULTS

3.1 | Clinical features of epilepsy-affected patients

P1 is a 5 ½ years old boy and the only affected child of non-
consanguineous parents (Figure S1a). Pregnancy and birth at
term were normal (birth weight 3,905 g, body lengths 53 cm,
head circumference 39.0 cm). Apgar values were normal with
9-10-10. Myoclonic jerks were observed at the age of 4 weeks.
The first electroencephalography (EEG) was performed at the
age of 8 weeks and showed no epileptic discharges. At the
age of 3 months P1 was transferred to University Children's
Hospital, Klinikum Oldenburg, Oldenburg because of in-
creased infantile spasms. Because EEG showed hypsarrhyth-
mia the West syndrome was diagnosed in P1. Extensive
diagnostic investigations, including cerebral magnetic

resonance imaging (MRI), metabolic examinations with cer-
ebrospinal fluid examination with neurotransmitter analy-
sis, genetic array comparative genomic hybridization (array
CGH), and panel diagnostic for epileptic encephalopathies (79
genes) were performed. However, no abnormalities were de-
tected. An extended antiepileptic therapy over the next years
with 11 different antiepileptic drugs and a ketogenic diet were
applied. However, these treatments did not give any long-
lasting positive therapeutic effects and the patient was finally
diagnosed with pharmaco-resistant epilepsy of the West syn-
drome type (Table 1, P1). Recently, (at the age of 5 ½ years)
P1 was treated with valproic acid, sulthiame, and cannabidiol.
Renew EEG examinations showed a too slow basic rhythm,
bioccipital delta waves, multifocal sharp- and sharp-slow-
waves predominantly from temporo-parietal. Furthermore, P1
showed a strong abnormal psychomotor development. At the
age of 5½ years P1 was unable to fix or follow objects with his
eyes, he was unable to grasp, to turn around or to sit freely. He
was not able to spell out and understand any words. The head
circumference declined over the time and the patient devel-
oped a microcephaly (48 cm, first percentile, -2.9 SD). MRI
examinations were additionally performed at the age of 4 and
7 months, and 3 and 4 years. MRI findings revealed a progres-
sive global cerebral atrophy.

P2 is 21 years old and the only affected child in a non-
consanguineous family. There is a familial history of intel-
lectual deficiency in the maternal branch, with multiple male
affected with intellectual disability (Figure S1b). P2 was
born at term after an uneventful pregnancy. He presented
global developmental delay, was able to sit alone at the age
of 10 months and walked at the age of 30 months. Mild static
and kinetic ataxia were observed since the age of 3 years, and
remained stable over years. He presented partial complex epi-
leptic seizures with temporal abnormalities on EEG at the age
of 6 years (Table 1). He was effectively treated with valproic
acid, levetiracetam, and topiramate, showing normal EEGs
at the age of 10 years. He also had a bicuspid aortic valve
and small initial aorta dilatation, treated with bisoprolol.
Dilatation remained stable over years. Hyperprolinemia was
noted during an amino acid chromatography testing, which
was treated with riboflavin. He followed adapted courses
during childhood and was finally referred to a center for dis-
abled persons.

P3 is a 21-year-old man and one of two children of unre-
lated parents (Figure S1c). He was born at a gestational age
of 42 weeks with a birth weight of 3,850 g. P3 pronounced
his first words at the age of 1.5 years, in a context of bi-
lingualism (Dutch and Russian). Initial concerns appeared
around the age of 2 years, with hyperactive behavior at the
age of 3 years. He also manifested speech delay, and autistic
behaviors. Intelligence quotient (IQ) testing revealed an IQ
score of 55. At the age of 18 years he had a first tonic-clonic
seizure, without clearly provoking factors. He was treated
<table><caption>TABLE 1 Clinical and genetic features of the investigated patients</caption><thead><tr><th>Patient ID</th><th>Gender/age of onset</th><th>Likely pathogenic variant</th><th>Diagnosis</th><th>Development</th><th>Seizure type</th><th>Medication</th><th>EEG</th><th>MRI</th></tr></thead><tbody><tr><td>P1</td><td>M/2 m</td><td>c.664G&gt;A<br>p.Glu222Lys</td><td>Severe pharmacoresistant epileptic encephalopathy, West syndrome</td><td>Abnormal psychomotor development, severe global developmental delay and intellectual disability, microcephaly</td><td>Initial myoclonic seizures, 3 m: infantile spasms</td><td>Pyridoxalphosphate, pulsatile dexamethasone, prednisolone, vigabatrin, topiramate, valproic acid, levetiracetam, ketogenic diet, oxcarbazepin (aggravation), lacosamide actual: valproic acid, sulthiame, cannabidiol</td><td>Hypsarrhythmia, slow basic rhythm, bioccipital delta waves, multifocal sharp- and sharp-slow-waves</td><td>Abnormal: progressive global cerebral atrophy</td></tr><tr><td>P2</td><td>M/6 y</td><td>c.1A&gt;G<br>p.Met1?</td><td>Focal epilepsy</td><td>Global developmental delay, mild static and kinetic ataxia, intellectual disability</td><td>Focal seizures</td><td>Valproic acid, levetiracetam, topiramate, sulthiame</td><td>Temporal abnormalities</td><td>Normal</td></tr><tr><td>P3</td><td>M/18 y</td><td>c.399C&gt;A<br>p.Tyr133*</td><td>Generalized epilepsy</td><td>Autistic behaviors, dysmorphic signs, intellectual disability</td><td>Tonic-clonic seizure</td><td>Valproic acid, levetiracetam, carbamazepine, risperidone, methylphenidate</td><td>Epileptiform abnormalities fronto-polar on both sides and left occipital with low index</td><td>NORMAL</td></tr><tr><td>P4</td><td>M/3 y</td><td>c.1045G&gt;A<br>p.Val349Ile</td><td>Generalized epilepsy</td><td>Global developmental delay, intellectual disability</td><td>Focal to bilateral seizures</td><td>Ospolot, urbanyl, sulthiame, clobazam</td><td>Spikes in contralateral hemisphere, bilateral wave discharges</td><td>normal</td></tr><tr><td>a₂ (Wang, Du, et al., 2017)</td><td>M/6 m</td><td>c.1355G&gt;T<br>p.Arg452Leu</td><td>Infantile spasms</td><td>Severe global developmental delay and intellectual disability</td><td>Spasms in cluster, infantile spasms</td><td>Valproic acid, topiramate, clonazepam</td><td>Hypsarrhythmia, multiple spike, and slow wave complex</td><td>Myelin development delay, brain dysplasia</td></tr></tbody></table>

Note: Patient a₂ described by Wang, Du, et al. (2017). Sequence variants were annotated using reference sequence NM_004961.3/NP_004952.2.

Abbreviations: EEG, electroencephalography; M, male; m, month; y, year.
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Identification of sequence variations in the GABRE gene in 4 GE-affected patients. (A and B) X-chromosomal localization of the human *GABRE* gene and its distribution of exons. (C) We found variants in exon 1, 4, 6, and 8, located in different protein domains. (D) Multiple species sequence alignments of GABRE proteins. The mutated GABRE protein sequences of the patients described herein were shown on top of the alignments. Consensus symbols in the alignments represent: aa, amino acid; SP, signal peptide; LCR, low complexity region; LBD, Pfam Neurotransmitter-gated ion-channel ligand binding domain; memb, Pfam Neurotransmitter-gated ion-channel transmembrane region; (x), single; fully conserved residues; (:), strongly similar residues; (.), weakly similar residues; (-), lack of identity between residues</figcaption></figure>

with valproic acid, levetiracetam, and carbamazepine. A valproic acid medication alone was insufficient to control his epilepsy (Table 1). To improve attention abilities and better control behavioral disturbances the patient was additionally treated at the age of 8 years with risperidone and methylphenidate. He did not show any visual or auditive impairment. Physical examination identified dysmorphic signs, similar to Marfanoid habitus, with narrow and high palate, long and slender fingers with slight camptodactyly, mild pectus exca-vatum, and stretch marks on the chest, shoulder, and lower back.

P4 was 10 years old at the time of the study. He was born after an uneventful pregnancy (Figure S1d). He presented global developmental delay during childhood and manifested behavior troubles. He suffered from focal to bilateral seizures since the age of 3 years (Table 1). EEG-measurements showed highlighted spikes in contralateral hemisphere as well as bilateral wave discharges. He was treated with sulthiame and clobazam. Brain MRI examination did not detect any abnormalities. He also presented hypothyroidism.

## 3.2 | Sequencing analyses

The DNA library from P1 was multiplexed with two different DNA libraries, which were not subjected to the project. The 3-plex sequencing run resulted in 31.8 Gb raw data and showed a Q30-score of 87%, cluster density of 336 K/mm², and cluster passing filter of 74.9%. Size of sequencing data obtained from P1 was estimated to 6.70 Gb, generating 66.86 M of clusters, 133.72 M of reads, and 92.5% of 10× targeted coverage. Applying stringent selection settings, the initial number of genetic variants was reduced to 90. We verified these variants for plausible types of inheritance, including autosomal-recessive and X-linked. Only one genetic variant passed the filtering criteria of X-linked mode of inheritance. The remaining 89 variants were excluded either due to an incorrect inheritance pattern or misguided disease-associations. We also considered heterozygous variants for possibly being de novo variants. None of these heterozygous variants was found in a promising candidate gene, to further verify its potential de novo appearance in familiar co-segregation analyses
with Sanger sequencing (Table S2). According to the hg19,
this genetic variant was located on X-chromosome, at posi-
tion 151'128'431. It was a missense variant (NM_004961.3:
c.664G>A, NP_004952.2: p.Glu222Lys) located in exon 6 of
the *GABRE* gene (Figures 1 and 2). Out of a total of 28 aligned
reads, 22 were sequenced from a reverse direction, whereas
6 reads were shown from a forward direction. All sequenced
reads contained the G-to-A exchange at the nucleotide posi-
tion 664 (Figure 1a). The allele frequency of the c.664G>A
(NM_004961.3) substitution, was estimated to 0.00004409
ExAc (Table S1). In silico predictions of the glutamic acid
into lysine substitution at the amino acid position 222 in the
GABRE protein suggested a likely pathogenic variant and
showed M-CAP = 0.387 and SIFT = 0.78 (VARVIS Version
1.12; Limbus Medical Technologies GmbH). Sanger sequenc-
ing was additionally performed in P1 and the available family

members. The analyses confirmed the co-segregation of the
c.664G>A (NM_004961.3) variant in the *GABRE* gene in the
family of P1 (Figure 1a; Figure S1a).

In family 2, a trio-WES was performed. The sequenced
data for P2 resulted in 9.69 Gb raw data and showed Q30-
score of 97%, and 95% of bases were covered at 25×. After
variant filtering, only one variant was identified as a potential
candidate for GE (Table S2). This variant was located on the
X-chromosome at position 151'974'625. The sequence variant
was detected in exon 1 of the *GABRE* gene and was thought
to affect the translation initiation codon (NM_004961.3:
c.1A>G, NP_004952.2: p.Met1?). Out of 63 aligned reads,
42 were sequenced from a reverse direction, whereas 21 reads
were shown from a forward direction. The allele frequency
of this sequence variant was estimated to 0.00001550 (ExAc,
Table S1). The single A-to-G substitution at nucleotide 1 of

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** Schematic presentation of pathogenicity predictions and localization of variants in the GABRE protein. (A) Landscape of the localization of the novel and previously identified variants in the *GABRE* gene and their predicted tolerance in the specific regions of the GABRE protein performed with MetaDome (NM_004961.3, NP_004952.2). Red areas represented highly intolerant regions, whereas blue regions showed highly tolerant areas. (B) The newly identified GABRE-variants are presented as red arrows (P2, NP_004952.2: p.Met1?; P3, NP_004952.2: p.Tyr133*; P1, NP_004952.2: p.Glu222Lys, intolerant; P4, NP_004952.2: p.Val349Ile, slightly intolerant). The previous GABRE-variants reported by Hernandez et al., 2016 (a₁: p.Arg452Gly, intolerant; b: p.Arg472His, tolerant, and c: p.Ser484Leu, neutral) and by Wang, Du, et al., 2017 (a₂: p.Arg452Leu, intolerant) are shown as white arrows. TM, transmembrane region (TM1: aa 254-274; TM2: aa 281-301; TM3: aa 344-364; TM4 aa 486-506); aa, amino acid; N, N-terminus; C, C-terminus; LBD, Pfam Neurotransmitter-gated ion-channel ligand binding domain; memb, Pfam Neurotransmitter-gated ion-channel transmembrane region.</figcaption></figure>
the *GABRE* gene in P2 was predicted to likely cause an alternative initiation codon. This substitution was found in all sequencing reads. The translation initiation codon ATG (coding for methionine) was changed to GTG (coding for valine). Consequently, this variant was predicted to result in an abnormal length of the GABRE protein (NP_004952.2: p.Met1?). This variant co-segregated with the GE phenotype within the family (Figure 1b).

A gene-panel approach for intellectual disabilities that was initially applied in P3 did not reveal any obvious causative variants. Therefore, trio-WES was considered to further search for the molecular cause of GE in P3. Around 95% of sequencing reads from each WES showed a coverage of 10×. Simultaneous analyses of WES data for P3 and his healthy parent revealed a nonsense variant in the *GABRE* gene (NM_004961.3: c.399C>A, NP_004952.2: p.Tyr133*). P3 was hemizygous for the mutated allele, whereas his mother was a carrier, and the father showed the reference allele (Figure 1c).

Variant filtering was applied in P4 and his healthy parents. A hemizygous missense variant, located in exon 8 of the *GABRE* gene was only found in P4, whereas the mother was carrier for this variant. As expected, the unaffected father showed reference allele (Figure 1d). The identified *GABRE*-variant in P4 was predicted to lead to an exchange of valine into isoleucine (NM_004961.3: c.1045G>A, NP_004952.2: p.Val349Ile) and was not found in the ExAc/GnomAD database. SIFT score and Mutation Tester showed scores of 0.92 and 0.947, respectively (Table S1). The score of the CADD was estimated to 6.512. These predictions support the pathogenicity of the variant in the *GABRE* gene detected in P4. In total, 20 x coverage of each WES was obtained for about 93% of sequencing reads.

## 3.3 | Bioinformatic predictions of the identified variants in the GABRE protein

To evaluate the evolutionary conservation of the GABRE protein, multiple protein alignments of different species were performed and showed a full conservation at amino acid 1, 133, 222, and 349 in regard to the human reference protein (NP_004952.2; P2, P3, P1, P4; Figure 2).

In silico predictions of the GABRE protein domains revealed that the genetic variants identified in P1 (NP_004952.2: p.Glu222Lys) and P3 (NP_004952.2: p.Tyr133*) were both located in the neurotransmitter-gated ion-channel ligand binding domain (PF02931; Figures 2c and 3b). This protein domain plays an important role in ligand binding, which in turn leads to changes in protein conformation and channel opening (Richter et al., 2012). The inflow of chloride into the cell is thought to result an inhibitory postsynaptic potential. The variant identified in P2 (NP_004952.2: p.Met1?) was located at the signal peptide site of the GABRE protein, which

determines the localization and the transport pathway of the protein within the cell (Blobel, 1980). The missense variant found in P4 (NP_004952.2: p.Val349Ile) was mapped to the neurotransmitter-gated ion-channel transmembrane region domain (PF02932). This protein domain is essential for ion channel formation (Figures 2c and 3b; Chen et al., 2018).

Tolerance analyses of the missense variants described herein (NM_004961.3: P1 c.664G>A, NM_004961.3: P4 c.1045G>A) revealed that the mutated amino acids were located in functional regions of the GABRE protein with intolerant and slightly intolerant impacts (Figure 3a). Moreover, we also tested the impact of previously reported missense variants (NP_004952.2: p.Arg452Leu, NP_004952.2: p.Arg452Gly, and NP_004952.2: p.Arg472His) in GABRE (Hernandez et al., 2016a; Wang, Du, et al., 2017). The two variants located at arginine 452 were predicted to be intolerant, while substitutions at arginine 472 and at serine 484 were predicted to be highly tolerant or neutral, respectively (Figure 3b). These computational analyses of variants allow us to better understand the possible pathogenicity of the GABRE-variants.

# 4 | DISCUSSION

Hereditary epilepsies belong to a prevalent and highly heterogeneous group of brain disorders (Fisher et al., 2014). They lead to sudden, erroneous and unprovoked electrical currents in neuronal networks and are frequently caused by pathogenic variants in many different genes (Guerreiro, 2016; Wang, Lin, et al., 2017).

In this study, we identified four families with four novel sequence variants in the *GABRE* gene, including one nonsense (P3: NM_004961.3: c.399C>A, NP_004952.2: p.Tyr133*), two missense variants (P1: NM_004961.3: c.664G>A, NP_004952.2: p.Glu222Lys; P4: NM_004961.3: c.1045G>A, NP_004952.2: p.Val349Ile), and one affecting a translational initiation codon (P2: NM_004961.3: c.1A>G, NP_004952.2: p.Met1?) in unrelated patients with GEs.

The GABAergic signaling is a major biological pathway and its abnormal regulation in neurons is known to be associated with numerous CNS disorders (Kang, 2017; Kumar, Sharma, Kumar, & Deshmukh, 2013). Pathogenic variants in GABA~A~ receptor subunit genes, including *GABRA*, *GABRB*, and *GABRG*, impair GABAergic signaling via postsynaptic mechanisms that lead to GEs (Butler et al., 2018; Kang, 2017). Electrophysiological and biochemical studies suggested that the majority of pentamer receptors are composed of two α subunits, two β subunits, and a single γ subunit (Y. Chang, Wang, Barot, & Weiss, 1996; Tretter, Ehya, Fuchs, & Sieghart, 1997). This assembly of GABA~A~ subunits is frequently found in a human brain (Bollan et al., 2008). The ε subunit seems to be evolutionary related to the γ subunit, showing an identity of amino acid restudies above 40% (Davies, Hanna, Hales, &
Kirkness, 1997). It has been speculated that the ε subunit is capable of replacing the γ subunit in the GABR pentamer, in order to form functional channels with distinct pharmacological and biophysical properties (Davies, Kirkness, & Hales, 2001; Saxena & Macdonald, 1994). This notion is supported by the finding that a replacement of α1 at the first position or β2 at the fourth position in the subunits' assembly also resulted in functional channel pores (Bollan et al., 2008). Previous studies have further shown that pathogenic variants in α, β, and γ subunits significantly impair structures and functions of GABA<sub>A</sub> receptors and lead to GEs (Hernandez et al., 2016b; Janve, Hernandez, Verdier, Hu, & Macdonald, 2016; Wang, Du, et al., 2017). The reported deleterious sequence variants have been mainly identified in N-terminal regions and trans-membrane protein domains, emphasizing an important biological function of these protein parts (Hernandez et al., 2016b). Interestingly, the *GABRE* sequence variants reported herein were also predicted to alter these protein regions supporting their involvement in GE.

Interactions of the ε subunit with α, β, and γ subunits are poorly understood and raise the question of their potential to form specific or unique neurotransmitter binding sites. To date, the impact of sequence variants located within the neurotransmitter-gated ion-channel ligand binding domain of ε subunits on a loss or gain of binding sites remains unclear. Unlike other GABR subtypes, the receptors containing a single ε subunit were described to act agonist-independent and to spontaneously lead to slower deactivation and recovery of the receptors (Davies et al., 2001; Maksay, Thompson, & Wafford, 2003; Neelands et al., 1999; Ranna, Sinkkonen, Moykkynen, Uusi-Oukari, & Korpi, 2006; Wagner et al., 2005). Pathogenic variants in a single α, β, and γ subunit of GABA<sub>A</sub> receptors showed negative effects on opening and closing of channel pores changing their structure and/or conductive functions, and therefore, are thought to result in epileptic phenotypes (Matsumoto & Ajmonemarsan, 1964; Prince & Wilder, 1967). Since pathogenic variants in the α1 subunit lead to GEs (Burgess et al., 2019; Hernandez et al., 2016b, 2019; Krenn et al., 2019), it is likely that an exchange of α1 with a mutated ε subunit will cause epileptic phenotypes.

It has been shown that various pathogenic variants in the same GE-associated gene may result either in different forms of neurodevelopmental disorders or cause a wide spectrum of epilepsies (Johannesen et al., 2016). In the patients described herein, we observed that sequence variants in the *GABRE* gene may cause a broad range of epileptic phenotypes, from an early-onset and severe pharmaco-resistant epileptic encephalopathy to late-onset and mild forms of generalized and focal epilepsy (Table 1, Figure 3). Currently, only handful patients with *GABRE*-variants were described worldwide and only very little details were published about the clinical symptoms of these patients. So far, only one patients with a missense variant (NM_004961.3: c.1355G>A, NP_004952.2:

p.Arg452Leu) in the *GABRE* gene was described to manifest infantile spasms, hypsarrhythmia, myelin development abnormalities, and brain dysplasia (Wang, Du, et al., 2017). We also found similar clinical symptoms in P1 affected with the West syndrome. Interestingly, most of the patients presented herein showed moderate forms of epilepsy with developmental delay and intellectual disability. These are novel clinical features that were unreported in patients having sequence variants in the *GABRE* gene (Hernandez et al., 2016b; Wang, Du, et al., 2017). Interestingly, a wide spectrum of epileptic phenotypes, ranging from early-onset severe forms, to late-onset mild forms, and to cases with intellectual disability and developmental delay with other neuropsychiatric symptoms were already reported in patients with mutations in other GABA<sub>A</sub> receptor subunit genes (Maljevic et al., 2019; Oyrer et al., 2018). However, due to a limited number of cases with sequence variants in the *GABRE* gene it is currently difficult to establish a clear phenotype-genotype correlation. Moreover, modifier variants located in uncovered coding and noncoding regions of known epilepsy-related genes, as well as variants in new genes, which cannot be excluded, may also impact on the phenotypic discrepancies observed in the investigated patients.

The available clinical data also suggested that the patients may respond differently to treatments with antiepileptic drugs. Epileptic seizures in P3 and P4 were either under control or even fully silenced in P2 due to applied therapies. However, seizures were not treatable or diminished by medications in P1, who showed pharmaco-resistant epilepsy. It has been reported that many GE-affected patients, similar to P1, did not respond to any drug therapies (Pierzchala, 2010). Consequently, epilepsies and epileptic syndromes remain life-threatening neurological disorders. An identification of new pathogenic variants in epilepsy-associated genes, and a discovery of novel candidate genes associated with epilepsies and neurodevelopmental disorders are essential for a deeper understanding of the genetic heterogeneity underlying these clinical symptoms. Additionally, in vitro and in vivo functional studies will be essential to enable the understanding of the molecular functions and pathogenicity of identified sequence variants, and thus, will contribute to a development of treatment options.

In this study, we identified novel sequence variants in the *GABRE* gene in four unrelated patients affected with GEs. The previously described genomic studies on large cohorts of patients affected with different forms of GE independently described four single cases with sequence variants in the *GABRE* gene (Hernandez et al., 2016b; Wang, Du, et al., 2017). However, the interpretation of these genetic findings was difficult due to the lack of comparable cases. Private pathogenic variants or sequence variants of unknown significance could not be excluded. Our genetic and clinical findings observed in P1 were similar to the previously reported case with
GABRE-variant (Wang, Du, et al., 2017) emphasizing an association of sequence variants in the GABRE gene with GE. We also collected three patients with GABRE-variants showing milder epileptic phenotypes as well as developmental delay and intellectual disability. In summary, this is the first report that showed the wide variety of epileptic phenotypes in very rare patients with sequence variants in the GABRE gene.

We postulate that the *GABRE* gene is a likely candidate gene for epilepsies of different forms and severities. Our studies lay the grounds for a better understanding of novel molecular genetic factors related to these clinical features. They may also complement the current knowledge of pharmac biochemical activities of GABA<sub>A</sub> receptors and may be useful to improve current therapies or develop novel antiepileptic drug treatments in the future.

5 | CONSENT FOR PUBLICATION

The consent for publication is available from the corresponding author on reasonable request.

ACKNOWLEDGMENTS

We thank Mrs. Kerstin Stegmann from Junior Research Group Genetics of Childhood Brain Malformations at the University Oldenburg for technical support.

CONFLICT OF INTEREST

The authors declared that they have no conflict of interest.

AUTHOR CONTRIBUTIONS

Fenja Markus: performed molecular research investigations and analyses (P1), mainly wrote the manuscript. Chloé Angelini: collected clinical data (P2), performed exome analyses (P2), and wrote parts of the manuscript. Aurelien Trimouille: shared clinical and genetic data (P2), performed clinical investigations and exome sequencing and data analyses (P2). Gabrielle Rudolf: performed molecular research investigations and analyses (P4). Gaetan Lesca: performed molecular research investigations and analyses (P4). Cyril Goizet: performed clinical investigations, shared clinical data (P2). Eulalie Lasseau: contributed to exome sequencing and analyses (P2). Benoit Arveiler: contributed to exome sequencing and analyses (P2). Marjon van Slegtenhorst: performed exome sequencing and data analyses (P3). Alice S. Brooks: performed clinical investigations (P3), shared clinical and genetic data (P3). Rami Abou Jamra: supported interpretations of clinical and genetic data. Georg-Christoph Korenke: performed clinical investigations (P1), shared clinical data. John Neidhardt: shared clinical and genetic data (P1), contributions to the manuscript. Marta Owczarek-Lipska: performed exome sequencing and analyses (P1), coordination of data exchange with cooperation partners and coordination of this study.

DATA AVAILABILITY STATEMENT

The data sets used and/or analyzed during this study are available from the corresponding author on reasonable request.

ORCID

Aurelien Trimouille ID https://orcid.

org/0000-0002-3457-5684

Marta Owczarek-Lipska ID https://orcid.

org/0000-0001-6667-4452

REFERENCES

Behar, T. N., Schaffner, A. E., Colton, C. A., Somogyi, R., Olah, Z., Lehel, C., & Barker, J. L. (1994). GABA-induced chemokinesis and NGF-induced chemotaxis of embryonic spinal cord neurons. *Journal of Neuroscience*, *14*(1), 29–38. https://doi.org/10.1523/JNEUROSCI.14-01-00029.1994

Behar, T. N., Schaffner, A. E., Scott, C. A., O'Connell, C., & Barker, J. L. (1998). Differential response of cortical plate and ventricular zone cells to GABA as a migration stimulus. *Journal of Neuroscience*, *18*(16), 6378–6387. https://doi.org/10.1523/JNEUROSCI.18-16-06378.1998

Berg, A. T., Berkovic, S. F., Brodie, M. J., Buchhalter, J., Cross, J. H., van Emde Boas, W., ... Scheffer, I. E. (2010). Revised terminology and concepts for organization of seizures and epilepsies: Report of the ILAE Commission on Classification and Terminology, 2005–2009. *Epilepsia*, *51*(4), 676–685. https://doi.org/10.1111/j.1528-1167.2010.02522.x

Berg, A. T., Chakravorty, S., Koh, S., Grinspan, Z. M., Shellhaas, R. A., Saneto, R. P., ... Hegde, M. (2018). Why West? Comparisons of clinical, genetic and molecular features of infants with and without spasms. *PLoS One*, *13*(3), e0193599. https://doi.org/10.1371/journal.pone.0193599

Blobel, G. (1980). Intracellular protein topogenesis. *Proceedings of the National Academy of Sciences of the United States of America*, *77*(3), 1496–1500. https://doi.org/10.1073/pnas.77.3.1496

Bollan, K. A., Baur, R., Hales, T. G., Sigel, E., & Connolly, C. N. (2008). The promiscuous role of the epsilon subunit in GABAA receptor biogenesis. *Molecular and Cellular Neurosciences*, *37*(3), 610–621. https://doi.org/10.1016/j.mcn.2007.12.011

Burgess, R., Wang, S., McTague, A., Boysen, K. E., Yang, X., Zeng, Q., ... Scheffer, I. E. (2019). Genetic landscape of epilepsy of infancy with migrating focal seizures. *Annals of Neurology*, *86*(6), 821–831. https://doi.org/10.1002/ana.25619

Butler, K. M., Moody, O. A., Schuler, E., Coryell, J., Alexander, J. J., Jenkins, A., & Escayg, A. (2018). De novo variants in GABRA2 and GABRA5 alter receptor function and contribute to early-onset epilepsy. *Brain*, *141*(8), 2392–2405. https://doi.org/10.1093/brain/awy171

Carvill, G. L., Weckhuysen, S., McMahon, J. M., Hartmann, C., Moller, R. S., Hjalgrim, H., ... Mefford, H. C. (2014). GABRA1 and STXBP1: Novel genetic causes of Dravet syndrome. *Neurology*, *82*(14), 1245–1253. https://doi.org/10.1212/WNL.0000000000000291

Chang, R. S., Leung, C. Y. W., Ho, C. C. A., & Yung, A. (2017). Classifications of seizures and epilepsies, where are we? – A brief historical review and update. *Journal of the Formosan Medical Association*, *116*(10), 736–741. https://doi.org/10.1016/j.jfma.2017.06.001

Chang, Y., Wang, R., Barot, S., & Weiss, D. S. (1996). Stoichiometry of a recombinant GABAA receptor. *Journal of Neuroscience*, *16*(17), 5415–5424.
Chen, Q., Wells, M. M., Arjunan, P., Tillman, T. S., Cohen, A. E., Xu, Y., & Tang, P. (2018). Structural basis of neurosteroid anesthetic action on GABAA receptors. *Nature Communications*, 9(1), 1–10. https://doi.org/10.1038/s41467-018-06361-4

Cossette, P., Liu, L., Brisebois, K., Dong, H., Lortie, A., Vanasse, M., ... Rouleau, G. A. (2002). Mutation of GABRA1 in an autosomal dominant form of juvenile myoclonic epilepsy. *Nature Genetics*, 31(2), 184–189. https://doi.org/10.1038/ng885

Davies, P. A., Hanna, M. C., Hales, T. G., & Kirkness, E. F. (1997). Insensitivity to anaesthetic agents conferred by a class of GABA(A) receptor subunit. *Nature*, 385(6619), 820–823. https://doi.org/10.1038/385820a0

Davies, P. A., Kirkness, E. F., & Hales, T. G. (2001). Evidence for the formation of functionally distinct alphabetagammaepsilonon GABA(A) receptors. *Journal of Physiology*, 537(Pt 1), 101–113. https://doi.org/10.1111/j.1469-7793.2001.0101k.x

Farrant, M., & Nusser, Z. (2005). Variations on an inhibitory theme: Phasic and tonic activation of GABA(A) receptors. *Nature Reviews Neuroscience*, 6(3), 215–229. https://doi.org/10.1038/nrn1625

Fisher, R. S., Acevedo, C., Arzimanoglou, A., Bogacz, A., Cross, J. H., Elger, C. E., ... Wiebe, S. (2014). ILAE official report: A practical clinical definition of epilepsy. *Epilepsia*, 55(4), 475–482. https://doi.org/10.1111/epi.12550

Gnirke, A., Melnikov, A., Maguire, J., Rogov, P., LeProust, E. M., Brockman, W., ... Nusbaum, C. (2009). Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing. *Nature Biotechnology*, 27(2), 182–189. https://doi.org/10.1038/nbt.1523

Guerreiro, C. A. (2016). Epilepsy: Is there hope? *Indian Journal of Medical Research*, 144(5), 657–660. https://doi.org/10.4103/ijmr. IJMR_1051_16

Hernandez, C. C., Klassen, T. L., Jackson, L. G., Gurba, K., Hu, N., Noebels, J. L., & Macdonald, R. L. (2016a). Correction: deleterious rare variants reveal risk for loss of GABAA receptor function in patients with genetic epilepsy and in the general population. *PLoS One*, 11(11), e0167264. https://doi.org/10.1371/journal.pone.0167264

Hernandez, C. C., Klassen, T. L., Jackson, L. G., Gurba, K., Hu, N., Noebels, J. L., & Macdonald, R. L. (2016b). Deleterious rare variants reveal risk for loss of GABAA receptor function in patients with genetic epilepsy and in the general population. *PLoS One*, 11(9), e0162883. https://doi.org/10.1371/journal.pone.0162883

Hernandez, C. C., XiangWei, W., Hu, N., Shen, D., Shen, W., Lagrange, A. H., ... Macdonald, R. L. (2019). Altered inhibitory synapses in de novo GABRA5 and GABRA1 mutations associated with early onset epileptic encephalopathies. *Brain*, 142(7), 1938–1954. https://doi.org/10.1093/brain/awz123

Hirose, S., Mitsudome, A., Okada, M., & Kaneko, S. (2005). Genetics of idiopathic epilepsies. *Epilepsia*, 46(s1), 38–43. https://doi.org/10.1111/j.0013-9580.2005.461011.x

Ishii, A., Kang, J. Q., Schornak, C. C., Hernandez, C. C., Shen, W., Watkins, J. C., ... Hirose, S. (2017). A de novo missense mutation of GABRB2 causes early myoclonic encephalopathy. *Journal of Medical Genetics*, 54(3), 202–211. https://doi.org/10.1136/jmedgenet-2016-104083

Janve, V. S., Hernandez, C. C., Verdier, K. M., Hu, N., & Macdonald, R. L. (2016). Epileptic encephalopathy de novo GABRB mutations impair gamma-aminobutyric acid type A receptor function. *Annals of Neurology*, 79(5), 806–825. https://doi.org/10.1002/ana.24631

Johannesen, K., Marini, C., Pfeffer, S., Moller, R. S., Dorn, T., Niturad, C. E., ... Maljevic, S. (2016). Phenotypic spectrum of GABRA1: From generalized epilepsies to severe epileptic encephalopathies.

*Neurology*, 87(11), 1140–1151. https://doi.org/10.1212/WNL.00000 00000003087

Kang, J. Q. (2017). Defects at the crossroads of GABAergic signaling in generalized genetic epilepsies. *Epilepsy Research*, 137, 9–18. https:// doi.org/10.1016/j.epilepsyres.2017.08.013

Kodera, H., Ohba, C., Kato, M., Maeda, T., Araki, K., Tajima, D., ... Matsumoto, N. (2016). De novo GABRA1 mutations in Ohtahara and West syndromes. *Epilepsia*, 57(4), 566–573. https://doi.org/10.1111/ epi.13344

Korff, C. M., & Scheffer, I. E. (2013). Epilepsy classification: A cycle of evolution and revolution. *Current Opinion in Neurology*, 26(2), 163–167. https://doi.org/10.1097/WCO.0b013e32835ee58e

Krenn, M., Ernst, M., Tomschik, M., Treven, M., Wagner, M., Westphal, D. S., ... Aull-Watschinger, S. (2019). Phenotypic variability of GABRA1-related epilepsy in monozygotic twins. *Annals of Clinical and Translational Neurology*, 6(11), 2317–2322. https://doi.org/10.1002/acn3.50895

Kumar, K., Sharma, S., Kumar, P., & Deshmukh, R. (2013). Therapeutic potential of GABA(B) receptor ligands in drug addiction, anxiety, depression and other CNS disorders. *Pharmacology, Biochemistry and Behavior*, 110, 174–184. https://doi.org/10.1016/j.pbb.2013.07.003

Macdonald, R. L., & Olsen, R. W. (1994). GABAA receptor channels. *Annual Review of Neuroscience*, 17, 569–602. https://doi.org/10.1146/ annurev.ne.17.030194.003033

Maksay, G., Thompson, S. A., & Wafford, K. A. (2003). The pharmacology of spontaneously open alpha(1)beta(3)epsilonon GABA(A) receptor-ionophores. *Neuropharmacology*, 44(8), 994–1002. https://doi.org/10.1016/S0028-3908(03)00116-3

Maljevic, S., Moller, R. S., Reid, C. A., Perez-Palma, E., Lal, D., May, P., & Lerche, H. (2019). Spectrum of GABAA receptor variants in epilepsy. *Current Opinion in Neurology*, 32(2), 183–190. https://doi.org/10.1097/WCO.0000000000000657

Matsumoto, H., & Ajmonemarsan, C. (1964). Cellular mechanisms in experimental epileptic seizures. *Science*, 144(3615), 193–194. https:// doi.org/10.1126/science.144.3615.193

Moller, R. S., Wuttke, T. V., Helbig, I., Marini, C., Johannesen, K. M., Brilstra, E. H., ... Maljevic, S. (2017). Mutations in GABRB3: From febrile seizures to epileptic encephalopathies. *Neurology*, 88(5), 483–492. https://doi.org/10.1212/WNL.0000000000003565

Moss, S. J., & Smart, T. G. (2001). Constructing inhibitory synapses. *Nature Reviews Neuroscience*, 2(4), 240–250. https://doi.org/10.1038/35067500

Neelands, T. R., Fisher, J. L., Bianchi, M., & Macdonald, R. L. (1999). Spontaneous and gamma-aminobutyric acid (GABA)-activated GABA(A) receptor channels formed by epsilon subunit-containing isoforms. *Molecular Pharmacology*, 55(1), 168–178.

Nieh, S. E., & Sherr, E. H. (2014). Epileptic encephalopathies: New genes and new pathways. *Neurotherapeutics: the Journal of the American Society for Experimental NeuroTherapeutics*, 11(4), 796–806. https:// doi.org/10.1007/s13311-014-0301-2

Oyrer, J., Maljevic, S., Scheffer, I. E., Berkovic, S. F., Petrou, S., & Reid, C. A. (2018). Ion channels in genetic epilepsy: from genes and mechanisms to disease-targeted therapies. *Pharmacological Reviews*, 70(1), 142–173. https://doi.org/10.1124/pr.117.014456

Pierzchala, K. (2010). Pharmacoresistant epilepsy – epidemiology and current studies. *Neurologia I Neurochirurgia Polska*, 44(3), 285–290. https://doi.org/10.1016/s0028-3843(14)60043-8

Prince, D. A., & Wilder, B. J. (1967). Control mechanisms in cortical epileptogenic foci. "Surround" inhibition. *Archives of Neurology*, 16(2), 194–202. https://doi.org/10.1001/archneur.1967.00470200082007
Ranna, M., Sinkkonen, S. T., Moykkynen, T., Uusi-Oukari, M., & Korpi, E. R. (2006). Impact of epsilon and theta subunits on pharmacological properties of alpha3beta1 GABAA receptors expressed in Xenopus oocytes. *BMC Pharmacology, 6*, 1. https://doi.org/10.1186/1471-2210-6-1

Richter, L., de Graaf, C., Sieghart, W., Varagic, Z., Mörzinger, M., de Esch, I. J. P., ... Ernst, M. (2012). Diazepam-bound GABAA receptor models identify new benzodiazepine binding-site ligands. *Nature Chemical Biology, 8*(5), 455-464. https://doi.org/10.1038/nchembio.917

Saxena, N. C., & Macdonald, R. L. (1994). Assembly of GABAA receptor subunits: Role of the delta subunit. *Journal of Neuroscience, 14*(11 Pt 2), 7077-7086. https://doi.org/10.1523/JNEUROSCI.14-11-07077.1994

Shen, D., Hernandez, C. C., Shen, W., Hu, N., Poduri, A., Shiedley, B., ... Macdonald, R. L. (2017). De novo GABRG2 mutations associated with epileptic encephalopathies. *Brain, 140*(1), 49-67. https://doi.org/10.1093/brain/aww272

Sieghart, W. (2000). Unraveling the function of GABA(A) receptor subtypes. *Trends in Pharmacological Sciences, 21*(11), 411-413. https://doi.org/10.1016/S0165-6147(00)01564-9

Symonds, J. D., Zuberi, S. M., & Johnson, M. R. (2017). Advances in epilepsy gene discovery and implications for epilepsy diagnosis and treatment. *Current Opinion in Neurology, 30*(2), 193-199. https://doi.org/10.1097/WCO.0000000000000433

Tretter, V., Ehya, N., Fuchs, K., & Sieghart, W. (1997). Stoichiometry and assembly of a recombinant GABAA receptor subtype. *Journal of Neuroscience, 17*(8), 2728-2737.

Wagner, D. A., Goldschen-Ohm, M. P., Hales, T. G., & Jones, M. V. (2005). Kinetics and spontaneous open probability conferred by the epsilon subunit of the GABAA receptor. *Journal of Neuroscience, 25*(45), 10462-10468. https://doi.org/10.1523/JNEUROSCI.1658-05.2005

Wang, J., Lin, Z.-J., Liu, L., Xu, H.-Q., Shi, Y.-W., Yi, Y.-H., ... Liao, W.-P. (2017). Epilepsy-associated genes. *Seizure, 44*, 11-20. https://doi.org/10.1016/j.seizure.2016.11.030

Wang, Y. M., Du, X. N., Bin, R., Yu, S. S., Xia, Z. Z., Zheng, G., ... Wang, Y. (2017). Genetic variants identified from epilepsy of unknown etiology in Chinese children by targeted exome sequencing. *Scientific Reports, 7*(1), 1-11. https://doi.org/10.1038/srep40319

Watanabe, M., Maemura, K., Kanbara, K., Tamayama, T., & Hayasaki, H. (2002). GABA and GABA receptors in the central nervous system and other organs. *International Review of Cytology, 213*, 1-47. https://doi.org/10.1016/s0074-7696(02)13011-7

Wiel, L., Baakman, C., Gilissen, D., Veltman, J. A., Vriend, G., & Gilissen, C. (2019). MetaDome: Pathogenicity analysis of genetic variants through aggregation of homologous human protein domains. *Human Mutation, 40*(8), 1030-1038. https://doi.org/10.1002/humu.23798

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the Supporting Information section.

**How to cite this article:** Markus F, Angelini C, Trimouille A, et al. Rare variants in the GABA<sub>A</sub> receptor subunit ε identified in patients with a wide spectrum of epileptic phenotypes. *Mol Genet Genomic Med.* 2020;8:e1388.

https://doi.org/10.1002/mgg3.1388