Application of Trio-Whole Exome Sequencing in Genetic Diagnosis and Therapy in Chinese Children With Epilepsy

Epilepsy is one of the most common neurological disorders in pediatric patients with other underlying neurological defects. Identifying the underlying etiology is crucial for better management of the disorder. We performed trio-whole exome sequencing in 221 pediatric patients with epilepsy. Probands were divided into seizures with developmental delay/intellectual disability (DD/ID) and seizures without DD/ID groups. Pathogenic (P) or likely pathogenic (LP) variants were identified in 71/110 (64.5%) patients in the seizures with DD/ID group and 21/111 (18.9%) patients in the seizures without DD/ID group (P < 0.001). Eighty-seven distinct P/LP single nucleotide variants (SNVs)/insertion deletions (Indels) were detected, with 55.2% (48/87) of them being novel. All aneuploidy and P/LP copy number variants (CNVs) larger than 100 Kb were identifiable by both whole-exome sequencing and copy number variation sequencing (CNVseq) in 123 of individuals (41 pedigrees). Ten of P/LP CNVs in nine patients and one aneuploidy variant in one patient (Patient #56, #47, XXY) were identified by CNVseq. Herein, we identified seven genes (NCL, SEPHS2, PA2G4, SLC35G2, MYO1C, GPR158, and POU3F1) with de novo variants but unknown pathogenicity that were not previously associated with epilepsy. Potential effective treatment options were available for 32 patients with a P/LP variant, based on the molecular diagnosis. Genetic testing may help identify the molecular etiology of early onset epilepsy and DD/ID and further aid to choose the appropriate treatment strategy for patients.

Introduction

Epilepsy is one of the most common neurological disorders with 50–100 million affected, and 2–4 million new cases diagnosed each year worldwide (Pitkänen et al., 2016). Epilepsy is a chronic disorder characterized by recurrent spontaneous seizures, and often begins in childhood. Repeated and refractory seizures cause decreased social participation, long-term cognitive impairment, and significantly lower quality of life (Nickels et al., 2016). A genetic basis for some forms of epilepsy was confirmed via gene mapping in families, and the specific mutations associated with epilepsy syndromes were identified in the 1990’s (Annegers et al., 1982; Scheffer and Berkovic, 1997; Myers and Mefford, 2015).

The genetic etiology of epilepsy may be monogenic, resulting from single-gene mutations. Mutations or variants in multiple genes are also important to cause epilepsy (Møller et al., 2015). Currently, epilepsy genetics can be broadly characterized into two categories: (i) genes and loci associated with primary epilepsy; and (ii) genes associated with neurological disorders where epilepsy may be one of the symptoms (Poduri and Lowenstein, 2011). High throughput sequencing technologies have contributed to explore novel epilepsy genes. To date, numerous pathogenic variants in several genes have been associated with epilepsy and seizures (Yang et al., 2019).

The development of next-generation sequencing have greatly increased our knowledge on the genetic changes occurring across the entire human genome, allowing for the rapid and efficient discovery of genes involved in many diseases. Whole-exome sequencing (WES) is a powerful tool for detecting variants, especially the single nucleotide variants (SNVs) and the small insertions and deletions (InDels). WES is intensively being applied to clinical practice due to its low cost, high diagnostic yields, and excellent advantages regarding the analysis of novel genes and their subsequent investigation.

Diagnostic genetic tests for these complex conditions are becoming increasingly important (Berg et al., 2019) as their clinical heterogeneity and molecular complexity pose a great challenge for their clinical diagnosis and subsequent treatment. In this study, we retrospectively analyzed the diagnostic yields of trio-WES in 221 pediatric patients with epilepsy of unclear etiology and explored novel possible pathogenic genes. We aimed to explore the P/LP variants in family and specifically focus on patients with developmental delay (DD)/intellectual disability (ID) or without DD/ID. Meanwhile, we also wanted to explore the treatment strategies based on molecular diagnosis. More importantly, analysis of novel epilepsy candidate genes was performed when no pathogenic mutations were clearly identified in the characterized genetic diseases. We identified several novel genes variations such as SCN1A, MECP2, and KCNT1, which were confined as pathogenic or likely pathogenic variants of epilepsy. Herein, our results suggest that the application of WES would benefit for defining epilepsy genetic factors and treatment strategies in the clinic.

Materials and Methods

Study Design and Sample Collection

The outline of the study design is illustrated in Figure 1. The inclusion criteria were: (1) occurrence of seizures or epilepsy before the age of 16 years-old, (2) epileptic syndromes/epileptic encephalopathy with unknown etiology, and (3) severe seizures in neonates or generalized epilepsy or intractable epilepsy in infancy with generalized tonic–clonic seizures. Patients would be excluded if the seizures were caused by non-genetic factors such as cerebral trauma, cerebral tumor, cerebral infection, cerebrovascular disorders, or diagnosed metabolic disorders. All patients underwent electroencephalogram (EEG) and magnetic resonance imaging. Epilepsy diagnoses and classifications were made by a pediatric neurologist following the criteria published by the International League Against Epilepsy. The phenotypic features of eligible patients were assessed by clinicians during the reviews of the medical records and classified according to the Human Phenotype Ontology (HPO) terms. Clinical information of the family members was obtained through face-to-face inquiries by investigators. Biological parentage was confirmed using the genomic data as described previously (Manichaikul et al., 2010).

FIGURE 1

Figure 1. Schematic summarizing the 221 probands from non-consanguineous pedigree. Two groups were considered in this study: epilepsy with DD/ID and epilepsy without DD/ID. Pathogenic or likely pathogenic variants were categorized as: De novo, incomplete penetrance, inherited from affected parents, autosomal recessive inheritance. DD, developmental delay; ID, intellectual disability; P, pathogenic; LP, likely pathogenic; VUS, variant uncertain significance; CNV, copy number variant; CNVseq, CNV sequencing; WES, whole-exome sequencing.

Patients with developmental delay (DD)/intellectual disability (ID) were diagnosed by the pediatric neurologists according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5). The developmental profile of patients <36 months of age was assessed via clinical observation, the Ages and Stages Questionnaires, Third Edition (ASD-3), and the Gesell Developmental Observation-Revised (GDO-R) assessment. The Wechsler Intelligence Scale, clinical observation, and Peabody picture vocabulary tests were used for patients aged from 3 to 10 years. ID was assessed by an IQ under 70 using the Wechsler Preschool and Primary Scale of Intelligence-Fourth Edition (WPPSI-IV) for patients between the ages of 4 and 6 years, and the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV) for patients aged over 7 years old.

Whole-Exome Sequencing (WES)

Whole-exome sequencing and bioinformatics analyses were performed following the previously proposed guidelines (MacArthur et al., 2014; Richards et al., 2015). In brief, peripheral venous blood (2∼4 mL) was collected from the patients and their family members. The genomic DNA was extracted using the Blood genome column medium extraction kit following the manufacturer’s instructions (Kangweishiji, China). Using genomic DNA, the exonic regions and flanking splice junctions of the genome were captured using the xGen Exome Research Panel v1.0 (IDT, Coralville, IA, United States). Finally, the libraries were sequenced on an Illumina NovaSeq 6000 series sequencer with the following parameters: PE150, minimum of 11.6 million reads. The sequencing was performed by the Beijing Chigene Translational Medicine Research Center Co., Ltd., Beijing, China.

Raw data were processed using the fastp tool to remove the adapters and filter out the low-quality reads. The paired-end reads were performed using a Burrows-Wheeler Aligner (BWA) against the Ensembl GRCh37/hg19 human reference genome (Li and Durbin, 2010). Both SNVs and small InDels were called using the Genomic Analysis Toolkit (GATK) software (version 4.1.7) (McKenna et al., 2010). The copy number variant (CNV) calling was based on the ExomeDepth algorithm. The total read count of the sample mapped to each exon in the same batch as described previously (Plagnol et al., 2012).

Variants were annotated using an online system, developed by Chigene which contains 35 public databases, while our in-house database contains WES data from 69015 individuals (Supplementary Table 1). Candidate SNVs/small InDels were confirmed by Sanger sequencing. A small CNV (<10 kb) would be considered if the phenotype was highly related to the candidate gene located in this CNV region; these were confirmed by quantitative polymerase chain reaction (qPCR). We classified the candidate variants according to the American College of Medical Genetics and Genomics (Richards et al., 2015) and Sequence Variant Interpretation Working Group international guidelines (SVI WG)¹.

Copy Number Variation Sequencing (CNVseq)

Copy number Variation Sequencing (CNVseq) was performed as previously described (Gao et al., 2019). Briefly, the genomic DNA was fragmented by sonication (Covaris, United States) into 200–300 bp fragments and checked using agarose gel electrophoresis. After genomic library preparation, DNA samples were subsequently sequenced on an Illumina NovaSeq 6000 series sequencer (Illumina, San Diego, CA, United States). Raw image files were processed using BclToFastq (Illumina) for the base calling and raw data generation. The reads were then mapped to the GRCh37/hg19 human reference genome using the BWA software (Li and Durbin, 2010). Variant calling for CNVs ≥100 kb was performed using an in-house pipeline, and the candidate CNVs were filtered and detected using public CNV databases (Decipher, ClinVar, OMIM, DGV, and ClinGen). The pathogenicity of CNVs was classified according to the American College of Medical Genetics and Genomics guidelines (Riggs et al., 2020).

Identification of Candidate Pathogenic de novo Variants

Variants (SNVs and InDels in coding region; canonical ± 1 or 2 splice sites) were considered to be candidate pathogenic de novo if they met the following criteria: (1) in patients with normal parental phenotype; (2) genotype call ratio > 0.3 and supporting read depth > 20; (3) minor allele frequency (MAF) < 0.0001 as reported in the Genome Aggregation Database (gnomAD); (4) Pathogenic variants were in the Ensembl canonical transcript.

Statistical Analysis

Categorical data are expressed in percentage and the comparisons between the groups were analyzed using the Pearson’s Chi-square test or a two-tailed Fisher’s exact test (for N < 40), in which a P-value smaller than 0.05 was considered to be statistically significant. We performed statistical analyses using SPSS software, version 25.0 (SPSS Inc., Chicago, IL, United States).

Results

Participant Demographics and Phenotypes

Two hundred and twenty one of unrelated patients (96 females and 125 males) and their families were recruited from our hospital between January 2016 and November 2019. Patients were from non-consanguineous families in Southeast China. The age of seizure onset ranged from 1 day after birth to 15-years old. Sixty-six of the patients had family history of seizures. Patients were divided into two groups: the seizures with DD/ID group and the seizures without DD/ID group (Figure 1). Clinical information of the patients was summarized in Supplementary Tables 2, 3.

Molecular Diagnosis Yields

We conducted WES to detect the epilepsy-associated gene variants. Pathogenic (P) or likely pathogenic (LP) variants were identified in 92 patients (92/221 = 41.6%), consisting of 87 distinct gene-level variants in eighty-two patients, 10 CNVs in nine patients, and one male patients with 47, XXY. These mutations were found in 71 patients in the group of seizures with DD/ID (71/110 = 64.5%) and 21 patients in the group of seizures without DD/ID group (21/111 = 18.9%), respectively (Figure 2A). The patients in the group of seizures with DD/ID had more P/LP mutations than those in the group of seizures without DD/ID (P < 0.001) (Figure 2B). Interestingly, the patients with seizures in DD/ID group under 1 year-old showed more P/LP variants than those of other groups. P/LP gene variants were identified in sixty patients under 1 year-old and 86.7% of the patients belong to DD/ID group. Intriguingly, 47 de novo variants and 29 novel variants were identified in 62 variants in the patients under 1-year-old (Table 1). Moreover, the number of de novo gene variants in patients with DD/ID was more than that of patients without DD/ID (Table 2). In addition, we also identified 13 variants of uncertain significance (VUS) in 4.5% (10/221) of the patients (Supplementary Table 3).

FIGURE 2

Figure 2. Molecular diagnosis yield in different ages. (A) The molecular diagnosis yield in the groups of epilepsy with DD/ID and without DD/ID during different ages of seizure onset. (B) The individuals containing gene variants in the groups of epilepsy with DD/ID and without DD/ID during different ages of seizure onset. The percentage in the diagram represents the diagnosis positive rate = (P + LP)/all individuals in this age.

TABLE 1

Table 1. Number of patients with P/LP gene variants.

TABLE 2

Table 2. Number of de novo and novel gene variants in the patient under 1 year-old.

Gene Variants

To better understand the epilepsy-associated gene variants, we categorized the gene variants as de novo, incomplete penetrance, inherited from affected parents, and autosomal recessive (AR) variants. We found 58 de novo gene variants and 48 novel variants (Table 3, Figure 3A, and Supplementary Table 2). Thirty-two genes were associated with DD/ID group. SCN1A was most frequently involved, followed by KCNQ2 and TSC2. PRRT2 was most frequently involved in the group of seizures without DD/ID (Figure 3B).

TABLE 3

Table 3. Forty-eight novel pathogenic or likely pathogenic variants.

FIGURE 3

Figure 3. Gene variants analysis. (A) The number of different variant types observed in patients. (B) Identification procedure of candidate pathogenic de novo variants. (C) Distribution of recurrent (≥2 patients) genes with pathogenic or likely pathogenic variants in the groups of epilepsy with DD/ID and without DD/ID. (D) The diagram of epilepsy pedigrees. ∘ represents female; □ represents male; 🌑 and ■ represents affected individuals; arrow represents probands; wt represents wild type.

It is interesting to explore the novel gene variants in the 129 epilepsy patients without P/LP variants. So we re-analyzed the data and selected 16 novel candidate genes from 14 individuals (Figure 3C and Supplementary Table 4); however, nine of the genes contain several de novo variants in the in-house control database, suggesting that these nine genes cannot be evaluated as the de novo genes associated with epilepsy. De novo variants in the other seven genes (NCL, SEPHS2, PA2G4, SLC35G2, MYO1C, GPR158, and POU3F1) were not found in the in-house control database. In addition, we found that GPR158 and POU3F1 are highly expressed in the nervous system according to the Genotype-Tissue Expression (GTEx) database². It suggests that these seven genes may be related to epilepsy, while the functions of them need to be further confirmed.

Moreover, we also acquired seven variants from affected parent families, including TSC1, PRRT2, TSC1, PROKR2, RYR2, GABRA1, and KCNMA1 genes (Figure 3D). Interestingly, a de novo variant, SCN2A c.668G>A was detected in patient #68 with epileptic encephalitis. Her brother also had SCN2A c.668G>A variant and showed hand clenching accompanied by slight shaking and up rolling of eyeballs. Then, we conducted ultra-deep sequencing (average deep: 20000×) and detected the SCN2A: c.668G>A variant in father’s oral formulas, urine and seminal fluid. The results showed that the mosaicism percentage of oral formulas, urine and seminal fluid were 13.14, 12.7, and 23.26%, respectively. We confirmed that SCN2A c.668G>A variant was paternal germ line mosaicism (Figure 3D and Supplementary Table 2).

CNVs by WES and CNVseq

In addition, we investigated CNVs in 123 of individuals (41 pedigrees). All aneuploidy and P/LP CNVs (>100 kb in size) were identifiable by WES and CNVseq (Figures 1, 4). Ten P/LP CNVs in nine patients and one aneuploidy variant in one patient (Patient #56, #47, XXY) were identified by CNVseq. Three CNVs were located in chromosome 16 and belong to 16p11.2 deletion syndrome. Two CNVs were duplications and eight CNVs were deletions, ranging from 411 to 12 Mb (Table 4 and Supplementary Tables 5, 6). Patient #73 carried two de novo CNVs, which manifests as developmental delay and seizures. One of the CNVs was a deletion and located in Chr4, including ZNF141, PIGG, PDE6B, and CPLX1 genes; another CNV was a duplication variant and located in Chr15, involving MEF2A, ADAMTS17, CERS3, LINS, ALDH1A3, and CHSY1 genes. It was demonstrated that CPLX1 gene functional aberration caused severe infantile myoclonic epilepsy and ID (Redler et al., 2017).

FIGURE 4

Figure 4. Locations of P/LP CNVs identified by WES and/or CNVseq. CNV, copy number variant; CNVseq, CNV sequencing; WES, whole-exome sequencing.

TABLE 4

Table 4. Pathogenic/likely pathogenic CNVs identified by WES and CNVseq in forty-one pedigrees.

Therapeutic Implications

There were specific therapeutic recommendations for 10 genes with P/LP variants in the current cohort, including ALDH7A1 (n = 1), DEPDC5 (n = 3), GRIN2A (n = 2), KCNQ2 (n = 5), SCN1A (n = 12), SCN2A (n = 5), SLC2A1 (n = 1), SLC6A1 (n = 2), TSC1 (n = 2), and TSC2 (n = 5). In this study, 32 patients were applicable drug selection based on molecular diagnosis. For example, the patient #24 was a 13 years-old female with mild ID. Prior to the genetic testing referral, she kept monthly seizures despite treatment with sodium valproate, levetiracetam, and lamotrigine. Genetic test revealed a missense variant in SLC2A1 (c.997C>T, p.R333W), which cause the GLUT1 deficiency syndrome. Then a ketogenic diet was initiated based on the genetic results. Surprisingly, the patient kept seizure-free developmental improvement (cognitive and behavioral) after treatment. She has subsequently been tapered off all of the anti-epileptic drugs. For patients with the SCN1A gene mutations, a combination of VPA and TMP improved the seizures effectively and the whole treatment process should not use sodium channel blockers. In addition, vigabatrin treatment decreased the seizure frequency and improved EEG in four of patients with TSC1 or TSC2 gene mutations. Oxcarbazepine was effective for five of patients with KCNQ2 gene mutations (Table 5).

TABLE 5

Table 5. Treatment strategies of 32 cases of patients based on WES diagnosis.

Discussion

Genetic factors were estimated to play a role in 70∼80% of epilepsy cases, especially in children and neonates (Hildebrand et al., 2013). Several studies have focused on the application of next-generation sequencing as a diagnostic tool for epilepsy (Veeramah et al., 2013; Dyment et al., 2015; Parrini et al., 2017). Recent cohort studies suggested that the diagnostic yield of WES varies from 23 to 42% in patients with epilepsy (Helbig et al., 2016; McTague et al., 2016; Costain et al., 2019; Snoeijen-Schouwenaars et al., 2019; Yang et al., 2019; Johannesen et al., 2020; Rochtus et al., 2020). In these studies, the phenotypes of patients varied widely, the inclusion and exclusion criteria of patients were also not consistent, and the pathogenic genes/pathways might be different. In the present study, the overall diagnostic yield was 41.6%. Further, we also found the diagnostic yield of the seizures with the DD/ID group to be higher than that in previous studies, especially in seizure onset under 1-year-old (∼78.4%) (Trump et al., 2016; Yang et al., 2019). It may be attributed to the following reasons. Firstly, our study analyzed SNVs, InDels, and CNVs, which can lead to a higher diagnostic yield. Secondly, the non-randomized selection/hospital-enrichment of the patients may lead to sampling bias. Pediatricians were likely to have subjective preference in the selection of patients with DD/ID for clinical genetic testing as it is easier to discover disease-related P/LP variants. Thirdly, for some patients in the seizures without DD/ID group, the patients under 4 years old might develop to ID later.

In the present study, thirteen of the patients carried P/LP variants that are inherited from unaffected parents (Supplementary Table 2), 53.8% (7/13) of them had autism spectrum disorder (ASD)/DD/ID (Figure 1). The unaffected phenotype of carriers was likely due to the incomplete penetrance, which was previously reported for the six genes: DEPDC5, SCN1A, PCDH19, PRRT2, GRIN2A, and NPRL2; 16p11.2 deletion and 16p12.1 microdeletion syndrome suggested that other modifier gene(s), as well as epigenetic or environmental factors, modulate the phenotype (Weiss et al., 2008; Girirajan et al., 2010; Dimova et al., 2012; Heron et al., 2012; Ishida et al., 2013; Lesca et al., 2013; Meng et al., 2015; Ricos et al., 2016). For example, the penetrance of DEPDC5 variants with different forms of focal epilepsy was incomplete, varying from 50 to 82% (Ishida et al., 2013; Ricos et al., 2016). In this study, the DEPDC5 gene variants were null variants (c.562+1G>T, c.2731G>T, c.484-1_c.485delGGT) carried by three patients with focal epilepsy (Patient #10, #46, and #84). Only patient #10 had epilepsy with DD. Actually the patient #10 carried two DEPDC5 variants, c.562+1G>T and c.2507A>G (p.Y836C), the latter was inherited from his asymptomatic father. We can’t confirm the DEPDC5 gene with an AR inheritance in our local database. So this bi-allelic defect may exacerbate the clinical symptoms and further studies are required to confirm the functions.

Some studies defined that therapeutic outcomes of epilepsy were mostly based on the effect of protein function, clinical observation, and literature reports (Schoonjans et al., 2017; Yang et al., 2019; Johannesen et al., 2020). Herein, we reported the choices of therapeutic intervention in 32 of patients were affected based on the genetic diagnosis and the symptoms of some patients were improved effectively (Table 5). As some types of epilepsies responded to particular antiepileptic medications, personalized therapeutic strategies will be the best choice of epilepsy therapy. The top three most frequently mutated genes were the same as reported in this cohorts, including SCN1A, KCNQ2, and TSC2 (Yang et al., 2019). In addition, we also found that PRRT2 heterogeneous variant was the most frequent mutated gene in the group of seizures without DD/ID. Moreover, we detected three cases with 16p11.2 deletion (includes PRRT2 gene). Two of the patients (#58 and #81) were with DD/ID and one patient (#78) hasn’t shown DD/ID (<4 years-old). Our results are consistent with the previous reported (Ebrahimi-Fakhari et al., 1993; Termsarasab et al., 2014). Ebrahimi-Fakhari et al. (1993) reported that the patients with PRRT2 heterogenous variants commonly exhibited epilepsy and paroxysmal movement disorders (PRRT2-associated paroxysmal movement disorders, RRT2-PxMD) without intellectual delay. While the individuals with 16p11.2 deletion, or with rare bi-allelic PRRT2 pathogenic variants exhibited DD/ID or ASD. Meanwhile, we identified seven genes with de novo variants in pathogenically uncertain patients, which included the previously reported GPR158 (OMIM: 614573) gene. GPR158 gene is related to seizures (Elmariah et al., 2014) and highly expressed in the nervous system. It may be a promising epilepsy candidate gene. The function of the other six genes (NCL, SEPHS2, PA2G4, SLC35G2, MYO1C, and POU3F1) is currently unknown in the nervous system and the gene variations were observed in only one individual. The functions of these genes need to be further confirmed. Briefly, WES could help physicians identify epilepsy-associated genes in early onset patients and further provide effective treatment in clinic and improve patients’ life quality.

In fact, this study still had several limitations. Firstly, our data revealed several de novo SNVs/InDels; however, mosaicism was not confirmed. Secondly, we didn’t identify the plausible causal mutations in more than half of the patients. It indicates that these patients may not be an aggregate of simple Mendelian disorders and therefore require further powerful tools to evaluate the disease elucidation. Thirdly, in our study, although we identified seven potential candidate genes related to epilepsy disease, there is currently not enough evidence to support their pathogenicity. Therefore, a more comprehensive testing tool and further genetic studies with larger cohorts are required to fully elucidate the underlying etiology. Meanwhile, functional tests are urgent for assessing the epilepsy-associated genes.

In conclusion, our study demonstrates that the simultaneous analysis of SNVs, InDels, and CNVs based on NGS data can provide a high diagnostic yield for epilepsy, especially for patients with DD/ID, age of seizure onset under 1-year-old. We further demonstrate the potential of genetic diagnosis impacts on choosing the optimal treatment strategy for these patients.

Data Availability Statement

The original contributions presented in the study are included in the article. Data on patients cannot be made fully accessible in accordance with local research ethics protocols. Further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the Children’s Hospital, Zhejiang University School of Medicine. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

TJ, JG, LJ, and FG designed the data collection instruments, collected the data, and carried out the initial analyses. TJ, YS, YY, and FG reviewed and revised the manuscript. LX, CZ, XS, WG, and XK conceptualized the study and coordinated and supervised data collection. All authors approved the final manuscript as submitted and agreed to be accountable for all aspects of the work.

Funding

This work is supported by the Key Research and Development Program of Zhejiang Province, Grant/Award Number: 2020C03038; the National Natural Science Foundation of China Project, Grant/Award Number: U1909209; and the Natural Science Foundation of Zhejiang Province, Grant/Award Number: LBY21H090002.

Conflict of Interest

WG, XK, YY was employed by company Beijing Chigene Translational Medical Research Center Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank the patients and their families for participating in our research. We also thank the many physicians and families for referring patients.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2021.699574/full#supplementary-material

Supplementary Table 1 | The information of the public databases used in Chigene database.

Supplementary Table 2 | Clinical and molecular features of positive patients with positive genetic diagnosis. M, male; F, female; Y, year(s); M, month(s); D, day(s). Diagnosis: BFIS, benign familial infantile seizures; CCHD, critical congenital heart disease; CPVT, catecholaminergic polymorphic ventricular tachycardia; DOOSE, Doose Syndrome; DS or DS-like, Dravet syndrome or Dravet syndrome-like; EAS, epilepsy-aphasia spectrum; EE, epileptic encephalopathy; FE, focal epilepsy; FFEVF, familial focal epilepsy with variable foci; GE, generalized epilepsy; GEFS+, genetic epilepsy with febrile seizures plus; GULT1-DS, glucose transporter type 1 deficiency syndrome; IT, immune thrombocytopenia; LGS, Lennox-Gastaut Syndrome; MMPSI, malignant migrating partial seizures of infancy; OS, Ohtahara syndrome; PKD, paroxysmal kinesigenic dyskinesia; RTT, Rett syndrome; TSC, tuberous sclerosis; WEST, West syndrome. FOS, focal onset; FS, febrile seizures; GOS, generalized onset; DD, Developmental delay: ASD, autism spectrum disorder. Treatment: ACTH, adrenocorticotropic hormone; CBZ, carbamazepine; CLB, clobazam; CZP, clonazepam; KD, ketogenic diet; LCS, lacosamide; LCT, L-carnitine; LEV, levetiracetam; LTG, lamotrigine; NZP, nitrazepam; OXC, oxcarbazepine; PA, prednisone acetate; PB, phenobarbitone; RPM, rapamycin; TPM, topiramate; VB6, vitamin B6; VGB, vigabatrin; VPA, valproate. AR, autosomal recessive; AD, autosomal dominant; XL, X-linked; XLD, X-linked dominant; XLR, X-linked recessive; P, positive; N, negative; nd, not determined.

Supplementary Table 3 | Clinical and molecular features of patients with variants of uncertain significance (VUS) identify in the cohort. M, male; F, female; Y, year(s); M, month(s); D, day(s); AR, autosomal recessive; AD, autosomal dominant; P, positive; N, negative; nd, not determined.

Supplementary Table 4 | Clinical and molecular features of patients with de novo variants identified in the negative cohort. M, male; F, female; Y, year(s); M, month(s); D, day(s); AR, autosomal recessive; AD, autosomal dominant; XLD, X-linked dominant; XLR, X-linked recessive.

Supplementary Table 5 | Sequencing data quality of whole-exome sequencing and CNVseq for the patients and their parents. CNVseq, copy number variant sequencing.

Supplementary Table 6 | Pathogenic/likely pathogenic CNVs identified by WES and CNVseq in the cohort. M, male; F, female.

Footnotes

1. ^ https://www.clinicalgenome.org/working-groups/sequence-variant-interpretation/
2. ^ https://www.gtexportal.org/home/

References

Annegers, J. F., Hauser, W. A., Anderson, V. E., and Kurland, L. T. (1982). The risks of seizure disorders among relatives of patients with childhood onset epilepsy. Neurology 32, 174–179. doi: 10.1212/wnl.32.2.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, A. T., Wusthoff, C., Shellhaas, R. A., Loddenkemper, T., Grinspan, Z. M., Saneto, R. P., et al. (2019). Immediate outcomes in early life epilepsy: A contemporary account. Epilepsy Behav. 97, 44–50. doi: 10.1016/j.yebeh.2019.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Costain, G., Cordeiro, D., Matviychuk, D., and Mercimek-Andrews, S. (2019). Clinical application of targeted next-generation sequencing panels and whole exome sequencing in childhood epilepsy. Neuroscience 418, 291–310. doi: 10.1016/j.neuroscience.2019.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Dimova, P. S., Kirov, A., Todorova, A., Todorov, T., and Mitev, V. (2012). A novel PCDH19 mutation inherited from an unaffected mother. Pediatr. Neurol. 46, 397–400. doi: 10.1016/j.pediatrneurol.2012.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dyment, D. A., Tetreault, M., Beaulieu, C. L., Hartley, T., Ferreira, P., Chardon, J. W., et al. (2015). Whole-exome sequencing broadens the phenotypic spectrum of rare pediatric epilepsy: a retrospective study. Clin. Genet. 88, 34–40. doi: 10.1111/cge.12464

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebrahimi-Fakhari, D., Moufawad El Achkar, C., and Klein, C. (1993). “PRRT2-associated paroxysmal movement disorders,” in GeneReviews(^®), © 1993-2021, University of Washington, Seattle, eds M. P. Adam, H. H. Ardinger, and R. A. Pagon (Seattle, WA: University of Washington).

Google Scholar

Elmariah, H., Garrett, M. E., Soldano, K. L., Ataga, K. I., Eckman, J. R., Telen, M. J., et al. (2014). Genes associated with survival in adult sickle cell disease. Blood 124:2719.

Google Scholar

Gao, C., Wang, X., Mei, S., Li, D., Duan, J., Zhang, P., et al. (2019). Diagnostic yields of trio-WES accompanied by CNVseq for rare neurodevelopmental disorders. Front. Genet. 10:485. doi: 10.3389/fgene.2019.00485

PubMed Abstract | CrossRef Full Text | Google Scholar

Girirajan, S., Rosenfeld, J. A., Cooper, G. M., Antonacci, F., Siswara, P., Itsara, A., et al. (2010). A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay. Nat. Genet. 42, 203–209. doi: 10.1038/ng.534

PubMed Abstract | CrossRef Full Text | Google Scholar

Helbig, K. L., Farwell Hagman, K. D., Shinde, D. N., Mroske, C., Powis, Z., Li, S., et al. (2016). Diagnostic exome sequencing provides a molecular diagnosis for a significant proportion of patients with epilepsy. Genet. Med. 18, 898–905. doi: 10.1038/gim.2015.186

PubMed Abstract | CrossRef Full Text | Google Scholar

Heron, S. E., Grinton, B. E., Kivity, S., Afawi, Z., Zuberi, S. M., Hughes, J. N., et al. (2012). PRRT2 mutations cause benign familial infantile epilepsy and infantile convulsions with choreoathetosis syndrome. Am. J. Hum. Genet. 90, 152–160. doi: 10.1016/j.ajhg.2011.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hildebrand, M. S., Dahl, H. H., Damiano, J. A., Smith, R. J., Scheffer, I. E., and Berkovic, S. F. (2013). Recent advances in the molecular genetics of epilepsy. J. Med. Genet. 50, 271–279. doi: 10.1136/jmedgenet-2012-101448

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishida, S., Picard, F., Rudolf, G., Noé, E., Achaz, G., Thomas, P., et al. (2013). Mutations of DEPDC5 cause autosomal dominant focal epilepsies. Nat. Genet. 45, 552–555. doi: 10.1038/ng.2601

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesen, K. M., Nikanorova, N., Marjanovic, D., Pavbro, A., Larsen, L. H. G., Rubboli, G., et al. (2020). Utility of genetic testing for therapeutic decision-making in adults with epilepsy. Epilepsia 61, 1234–1239. doi: 10.1111/epi.16533

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesca, G., Rudolf, G., Bruneau, N., Lozovaya, N., Labalme, A., Boutry-Kryza, N., et al. (2013). GRIN2A mutations in acquired epileptic aphasia and related childhood focal epilepsies and encephalopathies with speech and language dysfunction. Nat. Genet. 45, 1061–1066. doi: 10.1038/ng.2726

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., and Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589–595. doi: 10.1093/bioinformatics/btp698

PubMed Abstract | CrossRef Full Text | Google Scholar

MacArthur, D. G., Manolio, T. A., Dimmock, D. P., Rehm, H. L., Shendure, J., Abecasis, G. R., et al. (2014). Guidelines for investigating causality of sequence variants in human disease. Nature 508, 469–476. doi: 10.1038/nature13127

PubMed Abstract | CrossRef Full Text | Google Scholar

Manichaikul, A., Mychaleckyj, J. C., Rich, S. S., Daly, K., Sale, M., and Chen, W.-M. (2010). Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873. doi: 10.1093/bioinformatics/btq559

PubMed Abstract | CrossRef Full Text | Google Scholar

McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1101/gr.107524.110

PubMed Abstract | CrossRef Full Text | Google Scholar

McTague, A., Howell, K. B., Cross, J. H., Kurian, M. A., and Scheffer, I. E. (2016). The genetic landscape of the epileptic encephalopathies of infancy and childhood. Lancet Neurol. 15, 304–316. doi: 10.1016/s1474-4422(15)00250-1

CrossRef Full Text | Google Scholar

Meng, H., Xu, H. Q., Yu, L., Lin, G. W., He, N., Su, T., et al. (2015). The SCN1A mutation database: updating information and analysis of the relationships among genotype, functional alteration, and phenotype. Hum. Mutat. 36, 573–580. doi: 10.1002/humu.22782

PubMed Abstract | CrossRef Full Text | Google Scholar

Møller, R. S., Dahl, H. A., and Helbig, I. (2015). The contribution of next generation sequencing to epilepsy genetics. Expert Rev. Mol. Diagnost. 15, 1531–1538. doi: 10.1586/14737159.2015.1113132

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, C. T., and Mefford, H. C. (2015). Advancing epilepsy genetics in the genomic era. Genome Med. 7:91. doi: 10.1186/s13073-015-0214-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nickels, K. C., Zaccariello, M. J., Hamiwka, L. D., and Wirrell, E. C. (2016). Cognitive and neurodevelopmental comorbidities in paediatric epilepsy. Nat. Rev. Neurol. 12, 465–476. doi: 10.1038/nrneurol.2016.98

PubMed Abstract | CrossRef Full Text | Google Scholar

Parrini, E., Marini, C., Mei, D., Galuppi, A., Cellini, E., Pucatti, D., et al. (2017). Diagnostic targeted resequencing in 349 patients with drug-resistant pediatric epilepsies identifies causative mutations in 30 different genes. Hum. Mutat. 38, 216–225. doi: 10.1002/humu.23149

PubMed Abstract | CrossRef Full Text | Google Scholar

Pitkänen, A., Löscher, W., Vezzani, A., Becker, A. J., Simonato, M., Lukasiuk, K., et al. (2016). Advances in the development of biomarkers for epilepsy. Lancet Neurol. 15, 843–856. doi: 10.1016/s1474-4422(16)00112-5

CrossRef Full Text | Google Scholar

Plagnol, V., Curtis, J., Epstein, M., Mok, K. Y., Stebbings, E., Grigoriadou, S., et al. (2012). A robust model for read count data in exome sequencing experiments and implications for copy number variant calling. Bioinformatics 28, 2747–2754. doi: 10.1093/bioinformatics/bts526

PubMed Abstract | CrossRef Full Text | Google Scholar

Poduri, A., and Lowenstein, D. (2011). Epilepsy genetics–past, present, and future. Curr. Opin. Genet. Dev. 21, 325–332. doi: 10.1016/j.gde.2011.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Redler, S., Strom, T. M., Wieland, T., Cremer, K., Engels, H., Distelmaier, F., et al. (2017). Variants in CPLX1 in two families with autosomal-recessive severe infantile myoclonic epilepsy and ID. Eur. J. Hum. Genet. 25, 889–893. doi: 10.1038/ejhg.2017.52

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 17, 405–424. doi: 10.1038/gim.2015.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricos, M. G., Hodgson, B. L., Pippucci, T., Saidin, A., Ong, Y. S., Heron, S. E., et al. (2016). Mutations in the mammalian target of rapamycin pathway regulators NPRL2 and NPRL3 cause focal epilepsy. Ann. Neurol. 79, 120–131. doi: 10.1002/ana.24547

PubMed Abstract | CrossRef Full Text | Google Scholar

Riggs, E. R., Andersen, E. F., Cherry, A. M., Kantarci, S., Kearney, H., Patel, A., et al. (2020). Technical standards for the interpretation and reporting of constitutional copy-number variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet. Med. 22, 245–257. doi: 10.1038/s41436-019-0686-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Rochtus, A., Olson, H. E., Smith, L., Keith, L. G., El Achkar, C., Taylor, A., et al. (2020). Genetic diagnoses in epilepsy: the impact of dynamic exome analysis in a pediatric cohort. Epilepsia 61, 249–258. doi: 10.1111/epi.16427

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheffer, I. E., and Berkovic, S. F. (1997). Generalized epilepsy with febrile seizures plus. A genetic disorder with heterogeneous clinical phenotypes. Brain 120(Pt 3), 479–490. doi: 10.1093/brain/120.3.479

PubMed Abstract | CrossRef Full Text | Google Scholar

Schoonjans, A., Paelinck, B. P., Marchau, F., Gunning, B., Gammaitoni, A., Galer, B. S., et al. (2017). Low-dose fenfluramine significantly reduces seizure frequency in Dravet syndrome: a prospective study of a new cohort of patients. Eur. J. Neurol. 24, 309–314. doi: 10.1111/ene.13195

PubMed Abstract | CrossRef Full Text | Google Scholar

Snoeijen-Schouwenaars, F. M., van Ool, J. S., Verhoeven, J. S., van Mierlo, P., Braakman, H. M. H., Smeets, E. E., et al. (2019). Diagnostic exome sequencing in 100 consecutive patients with both epilepsy and intellectual disability. Epilepsia 60, 155–164. doi: 10.1111/epi.14618

PubMed Abstract | CrossRef Full Text | Google Scholar

Termsarasab, P., Yang, A. C., Reiner, J., Mei, H., Scott, S. A., and Frucht, S. J. (2014). Paroxysmal kinesigenic dyskinesia caused by 16p11.2 microdeletion. Tremor Other Hyperkinet. Mov. 4:274.

Google Scholar

Trump, N., McTague, A., Brittain, H., Papandreou, A., Meyer, E., Ngoh, A., et al. (2016). Improving diagnosis and broadening the phenotypes in early-onset seizure and severe developmental delay disorders through gene panel analysis. J. Med. Genet. 53, 310–317. doi: 10.1136/jmedgenet-2015-103263

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeramah, K. R., Johnstone, L., Karafet, T. M., Wolf, D., Sprissler, R., Salogiannis, J., et al. (2013). Exome sequencing reveals new causal mutations in children with epileptic encephalopathies. Epilepsia 54, 1270–1281. doi: 10.1111/epi.12201

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, L. A., Shen, Y., Korn, J. M., Arking, D. E., Miller, D. T., Fossdal, R., et al. (2008). Association between microdeletion and microduplication at 16p11.2 and autism. N. Engl. J. Med. 358, 667–675. doi: 10.1056/NEJMoa075974

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Kong, Y., Dong, X., Hu, L., Lin, Y., Chen, X., et al. (2019). Clinical and genetic spectrum of a large cohort of children with epilepsy in China. Genet. Med. 21, 564–571. doi: 10.1038/s41436-018-0091-8

PubMed Abstract | CrossRef Full Text | Google Scholar