ZFHX3 variants cause partial epilepsy of childhood and spasms with favorable outcomes
=====================================================================================

* Ming-Feng He
* Li-Hong Liu
* Sheng Luo
* Juan Wang
* Jia-Jun Guo
* Peng-Yu Wang
* Qiong-Xiang Zhai
* Su-Li He
* Dong-Fang Zou
* Xiao-Rong Liu
* Bing-Mei Li
* Hai-Yan Ma
* Jing-Da Qiao
* Peng Zhou
* Na He
* Yong-Hong Yi
* Wei-Ping Liao
* China Epilepsy Gene 1.0 Project

## Abstract

**Background** The *ZFHX3* gene is highly expressed in the developing brain and plays vital roles in embryonic development, cell proliferation, neuronal differentiation, and neuronal death. The association between the *ZFHX3* gene and human disease was not defined.

**Methods** Whole-exome sequencing was performed in a cohort of 378 patients with partial (focal) epilepsy of childhood. A *Drosophila Zfh2* knockdown model was used to validate the association between *ZFHX3* and epilepsy. The expression profile of the *ZFHX3* ortholog was assessed by RT-qPCR and the data from the Brainspan database.

**Results** Eight pairs of compound heterozygous variants in *ZFHX3* were identified in eight unrelated cases with partial epilepsy. The identified *ZFXH3* variants had low or no frequencies in the populations of gnomAD. The *ZFHX3* gene presented significantly higher excesses of variants in the case cohort, including a higher excess of biallelic variants than the expected number of East Asian populations, higher aggregated frequencies than that in gnomAD, and a higher frequency of compound heterozygous variants than that in the asymptomatic parent controls. In *Zfh2* knockdown flies, the incidence and duration of seizure-like behavior were significantly higher than those of the controls. The *Zfh2* knockdown flies exhibited more firing in excitatory neurons than the wild-type line in electrophysiological recordings. All patients with *ZFHX3* compound heterozygous variants presented focal seizures and focal discharges on EEGs. One patient experienced frequent nonconvulsive status epilepticus, and two patients evolved from early spasms; the three cases also had neurodevelopmental abnormalities. However, all patients achieved seizure-free. In *Drosophila*, the expression of *Zfh2* is high in larvae and decreased in pupae and early adults. In mice, *Zfhx3* is predominantly expressed in fetuses and decreases dramatically after birth. In humans, the data from the Brainspan database showed that *ZFHX3* is highly expressed in the embryonic period and decreased after birth.

**Conclusion** This study suggested that *ZFHX3* is potentially a novel causative gene of partial epilepsy of childhood and infantile spasms. The correlation between the outcome and gene expression stage provided insight into the underlying mechanism of the natural course of illness, potentially being helpful in the management of the patients.

Keywords
*   *ZFHX3* gene
*   partial epilepsy
*   trio-based WES
*   neurodevelopmental disorder
*   gene expression stage

## Introduction

Epilepsy is a common neurological disorder with age-dependent seizures. In children, the prevalence of epilepsy ranges from 0.5% to 1%,1 and approximately 68% of childhood patients present focal seizures.2 Epilepsy affects approximately 0.7% of infants,3 and approximately 10% of infantile patients with epilepsy present spasms.4 Genetic factors are believed to be the main cause of epilepsy, accounting for approximately eighty percent of cases.5 Previously, a number of causative genes were identified in children with partial (focal) epilepsy, such as *DEPDC5*6, *GRIN2A*7, *UNC13B*8, *HCFC1*9, *LAMA5*10, *BCOR*11, *CELSR1*12, *BRWD3*13, *FAT1*14, and *ELP4*15. The genes associated with infantile spasms include *KCNQ2*16, *KCNT1*17, *GRIN2B*18, *STXBP1*19, *TBC1D24*20, *SCN1A*21, *CDKL5*22, and *ARX*23. However, the etiology in most patients with childhood partial epilepsy and/or infantile spasms remains to be elucidated.

In this study, we conducted trio-based whole-exome sequencing (WES) in a cohort of 378 unrelated children with partial epilepsy. Compound heterozygous variants of *ZFHX3* were identified in eight unrelated patients with partial epilepsy, including one with frequent nonconvulsive status epilepticus and two who evolved from early spasms; the three patients with severe seizures also had neurodevelopmental abnormalities. However, all patients became seizure-free. A *Drosophila* model with *Zfh2* knockdown was established to investigate the association between *ZFHX3* and epilepsy. To explore the underlying mechanism of favorable outcomes, the genetic-dependent stages of *ZFHX3* orthologs were investigated, in which the expression of flies and mice was determined by RT-qPCR and that of humans was analyzed using data from the Brainspan database. This study suggests that the *ZFHX3* gene is potentially a novel causative gene of partial epilepsy of childhood and infantile spasms with/without developmental abnormalities. The development-dependent expression pattern of *ZFHX3* explains the natural course of the illness.

## Materials and methods

### Patients

Patients were identified in a cohort of 378 unrelated children with partial epilepsy without acquired causes who were recruited from five hospitals from 2019 to 2022, including the Second Affiliated Hospital of Guangzhou Medical University, Guangdong General Hospital, Shenzhen Children’s Hospital, Shantou Chaonan Minsheng Hospital, and the Affiliated Brain Hospital of Nanjing Medical University. The clinical information of affected individuals was collected, including sex, seizure onset age, seizure types and frequency, response to antiepileptic drugs (AEDs), family history, and general and neurological examinations. Long-term video-electroencephalography (EEG) was performed to monitor epileptic discharges. Brain magnetic resonance imaging (MRI) scans were performed to detect structural abnormalities. Epileptic seizures and epilepsy syndromes were diagnosed according to the criteria of the Commission on Classification and Terminology of the ILAE (1989, 2001, 2010, 2017, and 2022).

### Trio-based whole-exome sequencing and genetic analysis

Genomic DNA was extracted from blood samples of the probands and their parents (trios) and other available family members by using the Flexi Gene DNA kit (Qiagen, Hilden, Germany). Trios-based WES was performed with NextSeq500 sequencing instruments (Illumina, San Diego, CA, USA). Detailed sequencing methods were described in our previous studies.8 A case-by-case analytical approach was adopted to identify candidate causative variants in each trio.8 Primarily, the rare variants were prioritized with a minor allele frequency (MAF) < 0.005 in the gnomAD database ([gnomad.broadinstitute.org](http://gnomad.broadinstitute.org)). Then, potentially pathogenic variants were retained, including frameshift, nonsense, canonical splice site, initiation codon, in-frame variants, missense, and synonymous variants predicted to impact splicing. Last and importantly, possibly disease-causing variants in each individual were screened under five models: 1) epilepsy-associated gene; 2) de novo variant dominant; 3) autosomal recessive inheritance, including homozygous and compound heterozygous variants; 4) X-linked; and 5) cosegregation analysis if available. To identify novel potential epilepsy genes, genes with recurrently identified de novo, biallelic, hemizygous, and cosegregated variants were selected for further studies. *ZFHX3* appeared as a candidate gene with recurrent compound heterozygous variants in this cohort. All variants in *ZFHX3* in this study were validated by Sanger sequencing and annotated based on the transcript NM_006885.

### Analysis of the burden of variants

Three specific statistical methods were used to analyze the association between *ZFHX3* and epilepsy, including recessive burden analysis,24 aggregate frequency of variants,25 and frequency of compound heterozygous variants.

For the burden of recessive variants, the *P* value was calculated as [1-cumulative binomial probability [(n-1, N, R)] according to a previous report,24 where n is the observed biallelic variant number for *ZFHX3*, N is the number of trios (378 in this cohort), and R is the rate of *ZFHX3* variants by chance in populations. Considering that all patients were Han Chinese, the cutoff was set according to the MAF in the ExAC-East Asian population.

For aggregate frequency analysis, frequencies of identified variants between the case cohort and the controls were compared, including general and East Asian populations in the gnomAD database, general populations in the ExAC database, and the 33,444 persons without known neuropsychiatric conditions in the Epi25 WES Browser ([https://epi25.broadinstitute.org/](https://epi25.broadinstitute.org/)).26

For the control of compound heterozygous variants, we established a cohort of 1942 asymptomatic parents from trios, in whom the compound heterozygous variants were identified by detecting one of the paired variants in the child, based on the fact that one of the paired variants in a parent would transmit to the child. The frequency of identified compound heterozygous *ZFHX3* variants in the case cohort was compared with that in the controls, including the 1942 asymptomatic parent controls and the variant co-occurrence data in gnomAD27.

### Molecular structural analysis

Protein modeling was performed to predict the effects of missense variants on the molecular structure by using the I-TASSER tool ([https://zhanggroup.org/I-TASSER/](https://zhanggroup.org/I-TASSER/)). Protein structures were visualized and analyzed by using the PyMOL Molecular Graphics System (Version 2.5; Schrödinger, LLC; New York, USA). Protein stability changes of missense variants were predicted by using the I-Mutant Suite ([https://folding.biofold.org/cgi-bin/i-mutant2.0.cgi](https://folding.biofold.org/cgi-bin/i-mutant2.0.cgi)).

### *Drosophila* experiments

Flies were reared at 25°C and 60-70% humidity with a standard cornmeal diet under a 12:12 h light and dark cycle. The *UAS-Zfh2-RNAi* line (CG1449, TH01656.N) was purchased from the TsingHua Fly Center. The *tub-Gal4* driver line and *Canton-s* line were donated by Prof. Liu Ji-Yong (Guangzhou Medical University, Guangzhou, China). *Zfh2* knockdown flies were generated by crossing the *UAS-Zfh2-RNAi* line with the *tub-Gal4* driver line. *Canton-s* was the wild-type line. The knockdown efficiency was detected by using reverse transcription quantitative PCR (RT-qPCR).

To evaluate the role of Z*fh2* deficiency in development, the body length of fly larvae was assessed, as in our previous study.28 The larval stages of flies can be divided into three stages, i.e., the first (24 h), second (48 h), and third (72 h) instar larvae. The length of the larvae was measured using ImageJ NIH software ([https://imagej.nih.gov/ij/](https://imagej.nih.gov/ij/)).

To evaluate seizure-like behavior, a bang sensitivity assay was conducted on flies 3-5 days after eclosion.29 Briefly, flies were anesthetized using CO2 and transferred to clean vials for 0.5-1 h before testing. Three to seven flies were placed in each vial and stimulated mechanically using a vortex mixer (Lab Dancer, IKA, Germany) at maximum speed for 10 seconds.29 The duration and percentage of seizure-like behavior were recorded.

To determine the impact of *Zfh2* knockdown on neuronal excitability, the electrophysiological activity of projection neurons was recorded. The fly brains were prepared for electrophysiological recording as described previously.30 The standard external and internal solutions and patch pipette were prepared as described previously.31 Spontaneous activities were acquired using the 700B amplifier, 1440B Digital-Analog converter, and pClamp 10.5 software (Molecular Devices, San Jose, CA, USA). A cell with an access resistance of <30 MΩ sample was used for analysis. Spontaneous EPSP (sEPSP) data > 1 mV were analyzed by Mini Analysis software.

### Assessment of the *ZFHX3* ortholog expression profile

The mRNA expression levels of *ZFHX3* orthologs in different developmental stages in flies and mice were determined by RT-qPCR. For flies, the whole mRNA was extracted in five developmental stages, including third instar larvae, pupae, early adult (day 1), middle adult (day 5), and later adult (day 10). For mice, the mRNA of the frontal cortex was extracted in eight developmental stages, including the fetus, neonate (1 day), infant (1 week), toddle-period (2 weeks), preschool (4 weeks), juvenile-adult (10 weeks), middle-age (15 weeks), and old-age (32 weeks). The sequences of primers used in this study are listed in **Table S1**. Total RNA was extracted by using the HiPure Universal RNA Mini Kit (Magen Biotechnology, Guangzhou, China). Reverse transcription was performed with the HiScript III RT SuperMix for qPCR (+gDNA wiper) (Vazyme, Nanjing, China) kit, and then qPCR was performed by the Taq Pro Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) with the LightCycler 480 System (Roche) tool.

Human RNA-seq data at different developmental stages (from 8 postconceptional weeks to 40 years) for multiple brain areas were obtained from the Brainspan database ([http://www.brainspan.org/](http://www.brainspan.org/)). RNA expression was normalized to reads per kilobase million (RPKM). The temporal expression curve was fitted by third-order polynomial least squares to interpret the expression pattern of *ZFHX3* by GraphPad Prism 9.

### Statistical analysis

R statistical software (v4.0.2) and GraphPad Prism 9 were used for statistical analysis. All quantitative data are presented as the mean ± standard error of the mean (SEM). The burden of recessive variants was analyzed according to a previous report.24 The number and aggregated frequencies of *ZFHX3* variants in this cohort and controls were compared by a two-tailed Fisher’s exact test. Student’s t test was used to compare two independent samples, and the Mann-Whitney test was used to assess nonparametric data. The *P* value < 0.05 was considered statistically significant.

## Results

### Identification of *ZFHX3* variants

Compound heterozygous variants in *ZFHX3* were identified in eight unrelated cases with partial epilepsy, including c.314C>T/p.Pro105Leu & c.2282G>C/p.Gly761Ala, c.2419G>A/p.Glu807Lys & c.10853C>A/p.Pro3618Gln, c.2671T>C/p.Phe891Leu & c.9583\_9584insT/p.Pro3195LeufsTer44, c.2686G>A/p.Ala896Thr & c.10439C>T/p.Ala3480Val, c.4125\_4127del/p.Ala1376del & c.10853C>A/p.Pro3618Gln, c.5152A>C/p.Met1718Leu & c.10510G>C/p.Val3504Leu, c.6161C>T/p.Ala2054Val & c.10445G>T/p.Ser3482Ile, and c.10445G>T/p.Ser3482Ile & c.10853C>A/p.Pro3618Gln (**Table 1 and Figure 1**). The eight pairs of compound heterozygous variants consisted of eleven missense variants, one frameshift truncation, and one in-frame deletion. The missense variant p.Ser3482Ile was recurrently identified in two cases (Cases 7 & 8), and the variant p.Pro3618Gln was recurrently identified in three cases (Cases 2, 5, & 8). All of the *ZFHX3* compound heterozygous variants were inherited from their asymptomatic parents, consistent with the Mendelian autosomal recessive inheritance pattern.

![Figure 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/16/2023.07.16.23292551/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/F1)

Figure 1 Genetic data and representative EEG recordings of the cases with *ZFHX3* variants.
(A) Pedigrees of eight cases with compound heterozygous *ZFHX3* variants and their corresponding phenotypes. PE, partial epilepsy; FS, febrile seizures; NCSE, nonconvulsive status epilepticus; NDD, neurodevelopmental delay. (B) DNA sequencing chromatogram of *ZFHX3* variants. Arrows indicate the position of the variants. (C) EEG of case 2 showed right frontal-centro-temporal sharp-slow waves. (D) EEG of case 5 showed right temporal spike-slow waves. (E) EEG of case 6 showed right frontal spike-slow waves (left) and generalized poly-spike-slow waves (right). (F) EEG of case 8 showed bilateral occipital-temporal sharp-slow waves.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/T1)

Table 1. Clinical features of the individuals with *ZFHX3* variants

The identified variants presented no or low frequencies in the gnomAD database (MAF < 0.005) (**Table S2**). Six variants were not present in the normal control of the Epi25 WES Browser, while the other seven variants presented extremely low frequencies (MAF < 0.0005). None of the variants presented homozygous states in gnomAD.

When the burden of recessive variants was analyzed,24 the number of recessive *ZFHX3* variants identified in this cohort was significantly greater than the expected number by chance in the East Asian population (MAF < 0.005, *P* = 1.60×10-4). The aggregate frequency of variants identified in the case cohort was significantly higher than that in the controls, including the controls of the gnomAD-all population (*P* = 2.30×10-13), the controls of the gnomAD-East Asian population (*P* = 4.36×10−2), the general population of the ExAC database (*P* = 2.34×10-13), and the normal controls of the Epi25 WES Browser (*P* = 3.56×10-18) (**Table S2**). The eight pairs of compound heterozygous *ZFHX3* variants did not co-occur in gnomAD, and the frequency of compound heterozygous *ZFHX3* variants in the case cohort was significantly higher than that in the controls, including the control cohort of 1942 asymptomatic parents (8/378 vs 5/1942, *P* = 2.80×10-4) and the variant co-occurrence in gnomAD (8/378 vs 2/125748, *P* = 2.70×10-19).

None of the eight patients had pathogenic or likely pathogenic variants in the genes known to be associated with epileptic phenotypes.32

### Clinical features of the cases with *ZFHX3* variants

The detailed clinical characteristics of the eight unrelated cases with *ZFHX3* variants are summarized in **Table 1**. The onset age of seizures ranged from 5 months to 13 years old, with a median onset age of 5 years. The patients were all diagnosed with partial epilepsy, including three cases with neurodevelopmental abnormalities, of whom two presented frequent spasms (Cases 4 & 6) at the early stage of the illness and one experienced nonconvulsive status epilepticus (Case 7). The other five cases exhibited only infrequent (yearly) focal seizures or focal-origin generalized tonic-clonic seizures (Cases 1-3, 5, 8). Focal, multifocal, and/or diffuse epileptic discharges were recorded in their EEG (**Figure 1C-F**). The brain MRI was normal in all cases. A seizure-free status was achieved in all eight patients.

### Molecular effect and functional prediction of *ZFHX3* variants

The *ZFHX3* protein contains 23 zinc finger motifs and 4 homeodomains. Two variants, p.Glu807Lys and p.Ala1376del, were located in the zinc finger motifs, and other variants were scattered between zinc fingers and/or homeodomains (**Figure 2A**). The variant p.Pro3618Gln, located in the last part of the C-terminus, was repetitively identified in three of five mild cases (Cases 2, 5, & 8). Similarly, the frameshift variant (p.Pro3195LeufsTer44), which truncates the C-terminus, was also identified in a mild case (Case 3). The remaining mild case (Case 1) had two variants located in the N-terminus.

![Figure 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/16/2023.07.16.23292551/F2.medium.gif)

[Figure 2](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/F2)

Figure 2 Schematic illustration of *ZFHX3* variants.
(A) Schematic illustration of zinc finger homeobox protein 3 and locations of *ZFHX3* variants identified in this study. Two variants of the same height were a pair of biallelic variants. The red color represents cases with severe partial epilepsy with neurodevelopmental abnormalities. The blue color represents cases with mild partial epilepsy. (B) Hydrogen bond changes (red) and free energy stability changes (△△G, Kcal/mol) of the variants from the present study.

The damaging effect of missense variants was analyzed by protein modeling. Two missense variants, p.Glu807Lys and p.Ala896Thr, were predicted to alter hydrogen bonding with surrounding residues. Six variants were predicted to decrease protein stability with a LLG value less than −0.5 kcal/mol (**Figure 2B**). All the variants were predicted to be damaging by at least one *in silico* tool (**Table S3**).

### Knockdown of *Zfh2* in *Drosophila* led to increased susceptibility to seizures

To validate the association between *ZFHX3* and epilepsy, a *Drosophila* model of *Zfh2* knockdown was established, with a knockdown efficiency of approximately 66% (**Figure 3A**). The larval development of *Zfh2* knockdown flies was initially investigated. The larval body length of *Zfh2* knockdown flies was similar to that of *UAS-Zfh2-RNAi* control flies (**Figure 3B**).

![Figure 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/16/2023.07.16.23292551/F3.medium.gif)

[Figure 3](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/F3)

Figure 3 Functional studies on *Zfh2* knockdown flies
(A) Relative *Zfh2* mRNA expression levels of the knockdown flies and the *UAS-Zfh2-RNAi* controls. (B) The body length of *Zfh2* knockdown and *UAS-Zfh2-RNAi* control larvae. (C) Behavior in the Bang-sensitive test; the three phases observed in *Zfh2* knockdown flies were seizure, paralysis, and recovery. (D) Seizure-like behaviors occurred at a higher rate in *Zfh2* knockdown flies (*tub-Gal4 > UAS-Zfh2-RNAi*) than in the *UAS-Zfh2-RNAi* line. (E) The recovery time from seizure-like behaviors of *Zfh2* knockdown flies was longer than that of the controls.

A bang sensitivity assay was performed in *Zfh2* knockdown flies and *UAS-Zfh2-RNAi* control flies to assess susceptibility to seizures. The three phases of seizure activity in the Bang-sensitivity test were observed, including seizure, paralysis, and recovery (**Figure 3C**). The percentage of seizure-like behaviors in the *Zfh2* knockdown flies was significantly higher than that in the *UAS-Zfh2-RNAi* control [18.66 ± 2.09% (*n* = 15) vs 4.06 ± 0.87% (*n* = 9); *P* < 0.0001], the *Canton-s* control [18.66 ± 2.09% (*n* = 15) vs 6.88 ± 0.89% (*n* = 7); *P* = 0.0013], and the *tub-Gal4* control [18.66 ± 2.09% (*n* = 15) vs 3.56 ± 1.03% (*n* = 7); *P* < 0.0001] (**Figure 3D**). The duration of seizure-like behavior in *Zfh2* knockdown flies was also longer than that of the controls (**Figure 3E**).

The effect of *Zfh2* deficiency on the electrophysiological activity of projection neurons, which are important excitatory neurons in the central nervous system of *Drosophila*,33 was examined (**Figure 4A, B**). The frequency of sEPSPs in *Zfh2* knockdown flies was significantly higher than that in wild-type flies [0.6312 ± 0.06 Hz (*n* = 7) vs 0.09907 ± 0.01 Hz (*n* = 7); *P* < 0.0001] (**Figure 4C)**. There was no significant difference in the amplitude of sEPSPs between the *Zfh2* knockdown flies and WT flies [1.405 ± 0.2039 mV (*n* = 7) vs 1.228 ± 0.1486 mV (*n* = 7), *P* = 0.0888] (**Figure 4D**).

![Figure 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/16/2023.07.16.23292551/F4.medium.gif)

[Figure 4](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/F4)

Figure 4 *Zfh2* knockdown induces increased neural excitability in projection neurons.
(A) Whole-cell recording in the fly brain. (B) Representative traces of projection neuron sEPSPs in *Canton-s* wild-type (WT) flies and *Zfh2* knockdown flies. (C) The frequency of sEPSPs was significantly higher in *Zfh2* knockdown flies than in wild-type flies. (D) There was no difference in the amplitude of sEPSPs between *Zfh2* knockdown flies and wild-type flies.

### The temporal expression stage of *ZFHX3* orthologs

Infantile spasms are generally a severe form of epilepsy with a poor prognosis in the majority of cases.34 In this study, the two patients with early spasms presented favorable outcomes. Recent studies have indicated that the gene-dependent (expression) stage (GDS) is associated with the natural course and outcomes of illness.14, 35–37 We thus analyzed the temporal expression pattern of *ZFHX3* orthologs. In *Drosophila*, the expression of *Zfh2* was high in larvae, decreased in pupae and early adults, and increased in later adults (**Figure 5A**). In mice, the expression level of *Zfhx3* was high in the fetus, decreased dramatically after birth, and slightly increased in later adults (32 weeks) (**Figure 5B**). We further analyzed the temporal expression pattern of *ZFHX3* in the human brain by using data from BrainSpan. *ZFHX3* was also highly expressed in the embryonic stage, decreased dramatically in childhood with a nadir at approximately 10 years old, and slightly increased at approximately 30 years old (**Figure 5C**).

![Figure 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/16/2023.07.16.23292551/F5.medium.gif)

[Figure 5](http://medrxiv.org/content/early/2023/11/16/2023.07.16.23292551/F5)

Figure 5 The temporal expression profile of *ZFHX3* orthologs.
(A) Schematic illustration of the developmental stages of *Drosophila melanogaster* (left) and the temporal expression of *Zfh2* in flies (right). The mRNA levels of *Zfh2* were examined by RT-qPCR. The *Zfh2* mRNA levels were normalized to the mRNA levels of the third instar larvae. (B) Schematic illustration of the developmental stages of mice (left) and the temporal expression of *Zfhx3* in the frontal cortex of mice (right). The *Zfhx3* mRNA levels were normalized to the mRNA levels of fetal mice. (C) Temporal expression pattern of *ZFHX3* in the human brain. RPKM, reads per kilobase per million mapped reads.

## Discussion

The zinc finger homeobox 3 gene (*ZFHX3,* OMIM *104155), also known as AT motif-binding transcription factor 1 gene (*ATBF1*), encodes a transcription factor with 4 homeodomains and 23 zinc finger motifs. In this study, we identified compound heterozygous *ZFHX3* variants in eight unrelated patients with partial epilepsy, including two with early spasms. The *ZFHX3* gene presented significantly higher excesses of variants in the case cohort, including a higher excess of biallelic variants than the expected number of East Asian populations, a higher aggregated frequency than that in gnomAD, a higher frequency of compound heterozygous variants than that in the asymptomatic parent controls, and a higher frequency of compound heterozygous variants than the variant co-occurrence in gnomAD. Knockdown of *Zfh2* in flies caused increased susceptibility to seizures and abnormal firing of neurons. All patients achieved seizure-free, although two patients evolved from early spasms. The analysis of the temporal expression profile indicated that *ZFHX3* orthologs were highly expressed in the embryonic stage and decreased dramatically in childhood, which is correlated with the favorable outcomes of the patients. This study suggested that *ZFHX3* is potentially a novel causative gene of childhood partial epilepsy and infantile spasms with/without neurodevelopmental abnormalities.

The *ZFHX3* gene is highly conserved with homologs in flies, zebrafish, mice, and humans. It is ubiquitously expressed, especially in the developing brain. *ZFHX3* plays vital roles in multiple biological processes, including embryonic development, cell proliferation, neuronal differentiation, and neuronal death.38–40 In mice, homozygous knockout of *Zfhx3* led to prenatal lethality with complete penetrance, while heterozygous knockout resulted in growth restriction and postnatal lethality with incomplete penetrance.41 In zebrafish, knockdown of *Zfxh3* caused significantly increased pentylenetetrazol-induced seizures.42 In flies, previous studies showed that the majority of knockout/knockdown *Zfh2* lines exhibited preadult lethality ([http://flybase.org/reports/FBgn0004607#phenotypes](http://flybase.org/reports/FBgn0004607#phenotypes)). This study identifies biallelic *ZFHX3* variants in eight unrelated patients with partial epilepsy with/without neurodevelopmental abnormalities, and the knockdown of *Zfh2* in *Drosophila* led to increased susceptibility to seizures and abnormal firing of neurons. The recurrent epileptic and/or neurodevelopmental phenotypes across species supported the association between *ZFHX3* and epilepsy and neurodevelopmental abnormalities.

Data from variants in the human genome reveal that *ZFHX3* is intolerant to LOF variants, including constraint indexes such as the probability of being LOF intolerant (pLI),43 the LOF observed/expected upper bound fraction (LOEUF),43 and the probability of haploinsufficiency (pHaplo).44 The *ZFHX3* gene has a pLI of 0.997 (> 0.9), LOEUF of 0.147 (< 0.35), and pHaplo of 1 (≥ 0.86), all indicating that *ZFHX3* is intolerant to heterozygous LOF variants. Considering the recurrent lethal and neurodevelopmental/epileptic phenotypes in *Zfhx3* knockout/knockdown animals, it is possible that severe *ZFHX3* deficiency, such as biallelic null variants, will cause early lethality, while moderate *ZFHX3* deficiency, such as heterozygous null variants or biallelic missense variants with moderate damaging effects, is potentially associated with epilepsy and/or neurodevelopmental abnormalities. In this study, most of the variants were located outside the functional domains and predicted to be without hydrogen bond changes with surrounding residues, indicating mild damaging effects, which potentially explain the pathogenesis of compound heterozygous variants. Three of the patients with a mild phenotype had an identical variant p.Pro3618Gln, which was located in the last part of the C-terminus and presented a low MAF (0.000301) (instead of absent) in the general population, indicating a potentially mild damaging effect. One patient with a mild phenotype had a frameshift variant with truncation at the C-terminus, which was inherited from the asymptomatic father and presented with low MAF (instead of absent), thus potentially being of mild damage. These findings suggested a potential genotype-phenotype correlation.

The ZFHX3 protein, a transcription factor containing 23 C2H2-type zinc fingers and 4 homeodomains, belongs to the C2H2-type zinc finger protein family, which plays an essential role in neurodevelopment.45 Previously, *ZFHX3* de novo variants were occasionally detected in patients with neurodevelopmental disorders, including four variants in autism spectrum disorder,46–48 three variants in developmental disorder,49 and one variant in development and epileptic encephalopathy (**Table S4**),50 suggesting a potential role in neurodevelopmental disorders. In the present study, neurodevelopment abnormalities were observed in three patients. These data indicated a possible association between *ZFHX3* and neurodevelopmental disorders. Further studies are needed to determine the whole phenotypic spectrum of *ZFHX3*.

Infantile spasms are generally a severe form of epilepsy with poor outcomes.4 Clinically, approximately 33% of patients with spasms present favorable outcomes,34 but the underlying mechanism is unknown. Our recent study showed that GDS is associated with the natural course and outcomes of diseases.14, 35–37 In this study, the two patients with early spasms presented favorable outcomes. Studies on the temporal expression stage showed that *ZFHX3* orthologs were highly expressed in the embryonic stage and decreased dramatically in childhood, which is potentially one of the explanations for the favorable outcomes. Considering that the expression of *ZFHX3* was slightly increased at approximately 30 years old, long-term follow-up is needed to observe the future evolution of the disease.

In conclusion, this study suggested that *ZFHX3* is potentially a novel causative gene of partial epilepsy and infantile spasms with/without neurodevelopmental abnormalities. The correlation between the outcome and gene expression stage provided insight into the underlying mechanism of the natural course of illness, potentially being helpful in the management of the patients.

## Ethics declaration

This study was approved by the Ethics Committee of The Second Affiliated Hospital of Guangzhou Medical University (2020-h5-49), and written informed consent was obtained from all patients or their parents.

## Supplemental information

Supplemental Table S1: Primer sequences for RT-qPCR.

Supplemental Table S2: Analysis of the aggregate frequency of *ZFHX3* variants identified in this study.

Supplemental Table S3: Genetic features of the individuals with *ZFHX3* variants.

Supplemental Table S4: Previously reported *ZFHX3* variants and associated phenotypes.

Supplemental reference (for review only). *PCDH15* in self-limited partial epilepsy and the mechanism of distinct phenotypes (Revised in Annuals of Neurology).

## Supporting information

Table S1 [[supplements/292551_file07.docx]](pending:yes)

Table S2 [[supplements/292551_file08.docx]](pending:yes)

Table S3 [[supplements/292551_file09.docx]](pending:yes)

Table S4 [[supplements/292551_file10.docx]](pending:yes)

## Data Availability

Raw data were generated at Institute of Neuroscience, The Second Affiliated Hospital of Guangzhou Medical University. Derived data supporting the findings of this study are available from the corresponding author on request.

## Author contributions

W.P. L., Y.H. Y., N. H., and J.D. Q. designed this study; data collection was performed by M.F. H., L.H. L., S. L., J. W., Q.X. Z., S.L. H., D.F. Z., X.R. L., B.M. L., H.Y. M., and P. Z.; formal analysis was performed by M.F. H. and S. L.; experiments were performed by M.F. H. and J.J. G.; grant recipients were W.P. L., Y.H. Y., N. H., and X.R. L.; visualization was performed by M.F. H. and P.Y. W.; the first draft of the manuscript was written by M.F. H., L.H. L., and S. L.; writing-review and editing were performed by N. H., Y.H. Y., and W.P. L.

## Declaration of interests

The authors declare no conflicts of interest.

## Acknowledgments

We thank the patients and their families for participating in this study. We thank Prof. Liu Ji-Yong for donating the *tub-Gal4* driver line and *Canton-s* line flies and Schrödinger (New York, NY, USA) for supporting open-source PyMOL software. This study was supported by grants from the National Natural Science Foundation of China (Grant Nos. 82271505, 81870903, and 81971216); Science and Technology Project of Guangzhou (Grant No. 202201020106), Multicenter Clinical Research Fund Project of the Second Affiliated Hospital of Guangzhou Medical University (Grant Nos 010G271099, 2020-LCYJ-DZX-03, and 2021-LCYJ-DZX-02), and Scientific Research Project of Guangzhou Education Bureau (Grant No. 202235395). The funders had no role in the study design, data collection and analysis, or the decision to publish or in manuscript preparation.

## Footnotes

*   minor revision in main text and figure

*   Received July 16, 2023.
*   Revision received November 16, 2023.
*   Accepted November 16, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/)

## References

1.  1.Aaberg KM, Gunnes N, Bakken IJ, et al. Incidence and Prevalence of Childhood Epilepsy: A Nationwide Cohort Study. Pediatrics. May 2017;139(5)doi:10.1542/peds.2016-3908
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1542/peds.2016-3908&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28557750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

2.  2.Wirrell EC, Grossardt BR, Wong-Kisiel LC, Nickels KC. Incidence and classification of new-onset epilepsy and epilepsy syndromes in children in Olmsted County, Minnesota from 1980 to 2004: a population-based study. Epilepsy Res. Jun 2011;95(1-2):110–8. doi:10.1016/j.eplepsyres.2011.03.009
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eplepsyres.2011.03.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21482075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292434300012&link_type=ISI) 

3.  3.Eltze CM, Chong WK, Cox T, et al. A population-based study of newly diagnosed epilepsy in infants. Epilepsia. Mar 2013;54(3):437–45. doi:10.1111/epi.12046
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.12046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23252366&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000315709800008&link_type=ISI) 

4.  4.Zuberi SM, Wirrell E, Yozawitz E, et al. ILAE classification and definition of epilepsy syndromes with onset in neonates and infants: Position statement by the ILAE Task Force on Nosology and Definitions. Epilepsia. Jun 2022;63(6):1349–1397. doi:10.1111/epi.17239
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17239&link_type=DOI) 

5.  5.Dunn P, Albury CL, Maksemous N, et al. Next Generation Sequencing Methods for Diagnosis of Epilepsy Syndromes. Front Genet. 2018;9:20. doi:10.3389/fgene.2018.00020
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2018.00020&link_type=DOI) 

6.  6.Liu L, Chen ZR, Xu HQ, et al. DEPDC5 Variants Associated Malformations of Cortical Development and Focal Epilepsy With Febrile Seizure Plus/Febrile Seizures: The Role of Molecular Sub-Regional Effect. Front Neurosci. 2020;14:821. doi:10.3389/fnins.2020.00821
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2020.00821&link_type=DOI) 

7.  7.Lemke JR, Lal D, Reinthaler EM, et al. Mutations in GRIN2A cause idiopathic focal epilepsy with rolandic spikes. Nat Genet. Sep 2013;45(9):1067–72. doi:10.1038/ng.2728
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2728&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23933819&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

8.  8.Wang J, Qiao JD, Liu XR, et al. UNC13B variants associated with partial epilepsy with favourable outcome. Brain. Nov 29 2021;144(10):3050–3060. doi:10.1093/brain/awab164
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awab164&link_type=DOI) 

9.  9.He N, Guan BZ, Wang J, et al. HCFC1 variants in the proteolysis domain are associated with X-linked idiopathic partial epilepsy: Exploring the underlying mechanism. Clin Transl Med. Jun 2023;13(6):e1289. doi:10.1002/ctm2.1289
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ctm2.1289&link_type=DOI) 

10. 10.Luo S, Liu ZG, Wang J, et al. Recessive LAMA5 Variants Associated With Partial Epilepsy and Spasms in Infancy. Front Mol Neurosci. 2022;15:825390. doi:10.3389/fnmol.2022.825390
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnmol.2022.825390&link_type=DOI) 

11. 11.Li X, Bian WJ, Liu XR, et al. BCOR variants are associated with X-linked recessive partial epilepsy. Epilepsy Res. Nov 2022;187:107036. doi:10.1016/j.eplepsyres.2022.107036
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eplepsyres.2022.107036&link_type=DOI) 

12. 12.Chen Z, Luo S, Liu ZG, et al. CELSR1 variants are associated with partial epilepsy of childhood. Am J Med Genet B Neuropsychiatr Genet. Oct 2022;189(7-8):247–256. doi:10.1002/ajmg.b.32916
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.32916&link_type=DOI) 

13. 13.Tian MQ, Liu XR, Lin SM, et al. Variants in BRWD3 associated with X-linked partial epilepsy without intellectual disability. CNS Neurosci Ther. Feb 2023;29(2):727–735. doi:10.1111/cns.14057
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cns.14057&link_type=DOI) 

14. 14.Zou DF, Li XY, Lu XG, et al. Association of FAT1 with focal epilepsy and correlation between seizure relapse and gene expression stage. Seizure. Mar 11 2023;doi:10.1016/j.seizure.2023.03.003
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.seizure.2023.03.003&link_type=DOI) 

15. 15.Strug LJ, Clarke T, Chiang T, et al. Centrotemporal sharp wave EEG trait in rolandic epilepsy maps to Elongator Protein Complex 4 (ELP4). Eur J Hum Genet. Sep 2009;17(9):1171–81. doi:10.1038/ejhg.2008.267
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejhg.2008.267&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19172991&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

16. 16.Kato M, Yamagata T, Kubota M, et al. Clinical spectrum of early onset epileptic encephalopathies caused by KCNQ2 mutation. Epilepsia. Jul 2013;54(7):1282–7. doi:10.1111/epi.12200
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.12200&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23621294&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321110100017&link_type=ISI) 

17. 17.Ohba C, Kato M, Takahashi N, et al. De novo KCNT1 mutations in early-onset epileptic encephalopathy. Epilepsia. Sep 2015;56(9):e121–8. doi:10.1111/epi.13072
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.13072&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26140313&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

18. 18.Lemke JR, Hendrickx R, Geider K, et al. GRIN2B mutations in West syndrome and intellectual disability with focal epilepsy. Ann Neurol. Jan 2014;75(1):147–54. doi:10.1002/ana.24073
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.24073&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24272827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

19. 19.Saitsu H, Kato M, Mizuguchi T, et al. De novo mutations in the gene encoding STXBP1 (MUNC18-1) cause early infantile epileptic encephalopathy. Nat Genet. Jun 2008;40(6):782–8. doi:10.1038/ng.150
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.150&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18469812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256212000021&link_type=ISI) 

20. 20.Campeau PM, Kasperaviciute D, Lu JT, et al. The genetic basis of DOORS syndrome: an exome-sequencing study. Lancet Neurol. Jan 2014;13(1):44–58. doi:10.1016/S1474-4422(13)70265-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(13)70265-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24291220&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000329000900017&link_type=ISI) 

21. 21.Sadleir LG, Mountier EI, Gill D, et al. Not all SCN1A epileptic encephalopathies are Dravet syndrome: Early profound Thr226Met phenotype. Neurology. Sep 5 2017;89(10):1035–1042. doi:10.1212/WNL.0000000000004331
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000004331&link_type=DOI) 

22. 22.Archer HL, Evans J, Edwards S, et al. CDKL5 mutations cause infantile spasms, early onset seizures, and severe mental retardation in female patients. J Med Genet. Sep 2006;43(9):729–34. doi:10.1136/jmg.2006.041467
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiam1lZGdlbmV0IjtzOjU6InJlc2lkIjtzOjg6IjQzLzkvNzI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTEvMTYvMjAyMy4wNy4xNi4yMzI5MjU1MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

23. 23.Moey C, Topper S, Karn M, et al. Reinitiation of mRNA translation in a patient with X-linked infantile spasms with a protein-truncating variant in ARX. Eur J Hum Genet. May 2016;24(5):681–9. doi:10.1038/ejhg.2015.176
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejhg.2015.176&link_type=DOI) 

24. 24.Martin HC, Jones WD, McIntyre R, et al. Quantifying the contribution of recessive coding variation to developmental disorders. Science. Dec 7 2018;362(6419):1161–1164. doi:10.1126/science.aar6731
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjIvNjQxOS8xMTYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTEvMTYvMjAyMy4wNy4xNi4yMzI5MjU1MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

25. 25.Strande NT, Riggs ER, Buchanan AH, et al. Evaluating the Clinical Validity of Gene-Disease Associations: An Evidence-Based Framework Developed by the Clinical Genome Resource. Am J Hum Genet. Jun 1 2017;100(6):895–906. doi:10.1016/j.ajhg.2017.04.015
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2017.04.015&link_type=DOI) 

26. 26.Epi25 Collaborative. Electronic address sbuea, Epi C. Ultra-Rare Genetic Variation in the Epilepsies: A Whole-Exome Sequencing Study of 17,606 Individuals. Am J Hum Genet. Aug 1 2019;105(2):267–282. doi:10.1016/j.ajhg.2019.05.020
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2019.05.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31327507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

27. 27.Gudmundsson S, Singer-Berk M, Watts NA, et al. Variant interpretation using population databases: Lessons from gnomAD. Hum Mutat. Aug 2022;43(8):1012–1030. doi:10.1002/humu.24309
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/humu.24309&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34859531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

28. 28.Liu CQ, Qu XC, He MF, et al. Efficient strategies based on behavioral and electrophysiological methods for epilepsy-related gene screening in the Drosophila model. Front Mol Neurosci. 2023;16:1121877. doi:10.3389/fnmol.2023.1121877
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnmol.2023.1121877&link_type=DOI) 

29. 29.Parker L, Howlett IC, Rusan ZM, Tanouye MA. Seizure and epilepsy: studies of seizure disorders in Drosophila. Int Rev Neurobiol. 2011;99:1–21. doi:10.1016/B978-0-12-387003-2.00001-X
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/B978-0-12-387003-2.00001-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21906534&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

30. 30.Gu H, O’Dowd DK. Cholinergic synaptic transmission in adult Drosophila Kenyon cells in situ. J Neurosci. Jan 4 2006;26(1):265–72. doi:10.1523/JNEUROSCI.4109-05.2006
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1523/JNEUROSCI.4109-05.2006&link_type=DOI) 

31. 31.Qiao J, Yang S, Geng H, Yung WH, Ke Y. Input-timing-dependent plasticity at incoming synapses of the mushroom body facilitates olfactory learning in Drosophila. Curr Biol. Nov 21 2022;32(22):4869–4880 e4. doi:10.1016/j.cub.2022.09.054
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2022.09.054&link_type=DOI) 

32. 32.Wang J, Lin ZJ, Liu L, et al. Epilepsy-associated genes. Seizure. Jan 2017;44:11–20. doi:10.1016/j.seizure.2016.11.030
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.seizure.2016.11.030&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

33. 33.Ran D, Cai S, Wu H, Gu H. Di (2-ethylhexyl) phthalate modulates cholinergic mini-presynaptic transmission of projection neurons in Drosophila antennal lobe. Food Chem Toxicol. Sep 2012;50(9):3291–7. doi:10.1016/j.fct.2012.03.070
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.fct.2012.03.070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22490667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

34. 34.Chourasia N, Yuskaitis CJ, Libenson MH, et al. Infantile spasms: Assessing the diagnostic yield of an institutional guideline and the impact of etiology on long-term treatment response. Epilepsia. May 2022;63(5):1164–1176. doi:10.1111/epi.17209
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17209&link_type=DOI) 

35. 35.He N, Li B, Lin ZJ, Zhou P, Su T, Liao WP. Common genetic epilepsies, pathogenicity of genes/variants, and genetic dependence. Seizure. Jul 2023;109:38–39. doi:10.1016/j.seizure.2023.05.008
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.seizure.2023.05.008&link_type=DOI) 

36. 36.Li B, Bian W-J, Zhou P, et al. Genetic Dependence and Genetic Diseases. bioRxiv. 2023-08-05 2023;doi:10.1101/2023.08.02.551736
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMy4wOC4wMi41NTE3MzZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzExLzE2LzIwMjMuMDcuMTYuMjMyOTI1NTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

37. 37.Lin SM, Su T, Qiao JD, et al. PCDH15 in self-limited partial epilepsy and the mechanism of distinct phenotypes (Revised in Annuals of Neurology, attached to supplementary reference). 2023;
    
    

38. 38.Jung CG, Kim HJ, Kawaguchi M, et al. Homeotic factor ATBF1 induces the cell cycle arrest associated with neuronal differentiation. Development. Dec 2005;132(23):5137–45. doi:10.1242/dev.02098
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMToiMTMyLzIzLzUxMzciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMS8xNi8yMDIzLjA3LjE2LjIzMjkyNTUxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

39. 39.Ishii Y, Kawaguchi M, Takagawa K, et al. ATBF1-A protein, but not ATBF1-B, is preferentially expressed in developing rat brain. J Comp Neurol. Oct 6 2003;465(1):57–71. doi:10.1002/cne.10807
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cne.10807&link_type=DOI) 

40. 40.Jung CG, Uhm KO, Miura Y, et al. Beta-amyloid increases the expression level of ATBF1 responsible for death in cultured cortical neurons. Molecular neurodegeneration. Jul 5 2011;6:47. doi:10.1186/1750-1326-6-47
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1750-1326-6-47&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21729327&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

41. 41.Sun X, Fu X, Li J, et al. Heterozygous deletion of Atbf1 by the Cre-loxP system in mice causes preweaning mortality. Genesis. Nov 2012;50(11):819–27. doi:10.1002/dvg.22041
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/dvg.22041&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22644989&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

42. 42.Fuller TD, Westfall TA, Das T, Dawson DV, Slusarski DC. High-throughput behavioral assay to investigate seizure sensitivity in zebrafish implicates ZFHX3 in epilepsy. J Neurogenet. Mar-Jun 2018;32(2):92–105. doi:10.1080/01677063.2018.1445247
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/01677063.2018.1445247&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29718741&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

43. 43.Karczewski KJ, Francioli LC, Tiao G, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. May 2020;581(7809):434–443. doi:10.1038/s41586-020-2308-7
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

44. 44.Collins RL, Glessner JT, Porcu E, et al. A cross-disorder dosage sensitivity map of the human genome. Cell. Aug 4 2022;185(16):3041–3055 e25. doi:10.1016/j.cell.2022.06.036
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2022.06.036&link_type=DOI) 

45. 45.Al-Naama N, Mackeh R, Kino T. C(2)H(2)-Type Zinc Finger Proteins in Brain Development, Neurodevelopmental, and Other Neuropsychiatric Disorders: Systematic Literature-Based Analysis. Front Neurol. 2020;11:32. doi:10.3389/fneur.2020.00032
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fneur.2020.00032&link_type=DOI) 

46. 46.Wang T, Guo H, Xiong B, et al. De novo genic mutations among a Chinese autism spectrum disorder cohort. Nat Commun. Nov 8 2016;7:13316. doi:10.1038/ncomms13316
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms13316&link_type=DOI) 

47. 47.Hashimoto R, Nakazawa T, Tsurusaki Y, et al. Whole-exome sequencing and neurite outgrowth analysis in autism spectrum disorder. J Hum Genet. Mar 2016;61(3):199–206. doi:10.1038/jhg.2015.141
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/jhg.2015.141&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26582266&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

48. 48.Kosmicki JA, Samocha KE, Howrigan DP, et al. Refining the role of de novo protein-truncating variants in neurodevelopmental disorders by using population reference samples. Nat Genet. Apr 2017;49(4):504–510. doi:10.1038/ng.3789
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3789&link_type=DOI) 

49. 49.Deciphering Developmental Disorders S. Prevalence and architecture of de novo mutations in developmental disorders. Nature. Feb 23 2017;542(7642):433–438. doi:10.1038/nature21062
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature21062&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28135719&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 

50. 50.Epi KC, Epilepsy Phenome/Genome P, Allen AS, et al. De novo mutations in epileptic encephalopathies. Nature. Sep 12 2013;501(7466):217–21. doi:10.1038/nature12439
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12439&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23934111&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F16%2F2023.07.16.23292551.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000324244900040&link_type=ISI)