ORIGINAL RESEARCH article

Front. Pediatr., 27 February 2025

Sec. Pediatric Neurology

Volume 13 - 2025 | https://doi.org/10.3389/fped.2025.1471965

Evaluation of the etiology of epilepsy and/or developmental delay in children via next-generation sequencing: a single-center experience

  • 1. Department of Pediatrics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

  • 2. Division of Pediatric Genetics, Department of Pediatrics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

  • 3. Acibadem Mehmet Ali Aydinlar University Rare Diseases and Orphan Drugs Application and Research Center (ACURARE), Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

  • 4. Department of Transitional Medicine, Health Sciences Institute, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

  • 5. Acibadem Labgen Genetic Diagnosis Center, Istanbul, Türkiye

  • 6. Department of Genome Studies, Health Sciences Institute, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

  • 7. Division of Pediatric Neurology, Department of Pediatrics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye

Abstract

Background:

We aimed to understand the genetic etiology in children presenting with epilepsy and/or developmental delay by using next-generation sequencing (NGS).

Materials and methods:

We included children presenting to our pediatric neurology clinic with a diagnosis of epilepsy and/or developmental delay between January 2019 and December 2021. We evaluated the patients using the NGS equipment in our genetic laboratory.

Results:

In total, 90 patients were included in the study. Twenty (34.4%) out of 58 patients who had undergone whole-exome sequencing (WES) had pathogenic or likely pathogenic (P/LP) variants and 11 (18.9%) had variants of unknown significance (VUS). Five (41.6%) out of 12 patients who had undergone whole-genome sequencing had P/LP variants and 5 (41.6%) had VUS. Eleven (55%) out of 20 patients who had undergone WES and chromosomal microarray had P/LP variants and 2 (10%) had VUS. Twenty-six novel variants were described. Twelve patients (13.3%) were diagnosed using a known specific treatment.

Conclusion:

NGS aids in precisely diagnosing children with epilepsy and/or developmental delay. Furthermore, it provides a correct prognosis, specific treatment methods, and a multidisciplinary approach.

Introduction

The incidence of epilepsy in children ranges from 0.5 to 8 per 1,000 individuals (). Its estimated lifetime prevalence is 1%, and the highest incidence is observed during the first year of life (). Various developmental disorders, such as developmental delay (DD), intellectual disability (ID), autism spectrum disorder (ASD), and attention deficit and hyperactivity disorder (ADHD) are more common in children with epilepsy ().

In developed countries, the absence of expected psychomotor development is observed in approximately 1% of children (). DD and ID may occur alone or in conjunction with congenital malformations or neurological findings such as epilepsy, hypotonia, ataxia, electroencephalogram (EEG) abnormalities, and behavioral problems such as ASD and ADHD. DD and ID are phenotypically and genetically heterogeneous, and a specific diagnosis cannot be reached in most cases. Therefore, genetic tests play a crucial role in determining the cause of abnormality in this patient group ().

Next-generation sequencing (NGS) technologies enable the reliable identification of candidate genes responsible for epilepsy and DD in a single, cost-effective, and time-efficient process (). There are compelling reasons to make a specific diagnosis in a child presenting with DD. Diagnosis permits families to understand the long-term prognosis, access special education services, and learn about potential treatment options. Moreover, it is crucial for counseling parents regarding recurrence risk, prenatal diagnosis, and preimplantation genetic diagnosis. Furthermore, a diagnosis can help parents access research studies and connect with other families facing similar challenges ().

Chromosomal microarray (CMA) has long been utilized as a first-line genetic test in children with DD and ID. With the advancement of NGS techniques, whole-exome sequencing (WES) and whole-genome sequencing (WGS) have become pivotal diagnostic tools for children with DD and ID. The American College of Medical Genetics and Genomics (ACMG) guidelines strongly recommend considering WES/WGS as a first- or second-line test for patients with DD and ID ().

NGS has unveiled new candidate genes implicated in the pathogenesis of epilepsy. Ion channels play a crucial role in maintaining proper neuronal excitability; these include voltage-gated sodium, potassium, calcium, and chloride ion channels, and ligand-gated acetylcholine receptor-mediated ion channels. Mutations in genes encoding components of these ion channels lead to ion channel dysfunction, known as channelopathy, and form the basis for the development of epilepsy syndromes. Genetic changes causing channelopathies may contribute to the pathogenesis of epilepsy via the gain or loss of ion channel function ().

Early genetic testing and diagnosis are foundational for accurate treatment of epilepsy. However, they should be complemented with functional tests to explore the characteristics of the pathophysiological mechanisms and the functional impact underlying the detected variant. Specifically in ion channel disorders, understanding gain-of-function mutations and loss-of-function mutations is critical when making therapeutic decisions (, ).

Experience with most of the above mentioned treatments is limited to a relatively small number of patients. With a better understanding of the pathogenesis and cellular electrophysiology of genetic epilepsies, more treatment options are likely to emerge ().

Materials and methods

This retrospective cohort study included 90 patients aged 0–18 years, who were observed at the Acibadem Altunizade Hospital Pediatric Neurology Outpatient Clinic between January 2019 and December 2021 and had been diagnosed with epilepsy and/or DD. These patients had undetermined etiology despite undergoing biochemical and metabolic tests and magnetic resonance imaging (MRI). The genetic analyses were performed as part of the diagnostic process when patients first presented for evaluation. Informed consent for genetic testing was obtained by the accredited diagnostic laboratory prior to the tests. The WES and CMA results were processed at the Acibadem Labgen Genetic Diseases Evaluation Center, and the WGS results were obtained from the Centogene Laboratory (Germany).

Retrospective evaluation of the WES, CMA, and WGS results was conducted in collaboration with the Pediatric Genetics Department of Acibadem University. Patient medical records were retrospectively reviewed; these encompassed their medical history; background; family history; neurological symptoms; motor and mental developmental history; neurological examination findings; liver and kidney function test results; metabolic screening test, EEG, and MRI results; antiepileptic treatments; and genetic diagnoses.

Ethical approval for the retrospective study was obtained from the Scientific Research Ethics Committee of Acibadem University School of Medicine (Decision No: 2021/25/26).

Whole-exome sequencing

Genomic DNA isolation was performed using peripheral blood leukocytes. Enrichment was performed using the Twist® Human Core Exome kits. Prepared libraries were sequenced using the Novaseq 6000 (Illumina) instrument. Using the Burrows–Wheeler Aligner (BWA-MEM; Sentieon), sequenced data were aligned with genome reference to GRCh37 (hg19) (). Duplicate marking, indel realignment, baseline quality recalibration, and variant calling were performed using the Genome Analysis Toolkit (GATK) algorithms (Sentieon) (). Variant annotation, filtering, and interpretation were performed with the VarSome Clinical platform (Spatter). Variants in the coding and splicing regions and known pathogenic variants in the noncoding regions were included in the analysis. Variants with minor allele frequencies (below 1%) in the gnomAD database were evaluated. Filtering and prioritization of potentially disease-causing variants were performed based on the inheritance pattern, information in the ClinVar database, Human Phenotype Ontology (HPO) terms associated with the patient's findings, and in silico pathogenicity estimates. Variants were classified according to the ACMG guidelines () and ClinGen recommendations (), pathogenic (P), and likely pathogenic (LP), and clinically relevant variants of uncertain significance (VUS) associated with the patient's clinical findings were reported.

WGS

Genomic DNA was fragmented enzymatically, and libraries were generated using polymerase chain reaction-mediated addition of Illumina-compatible adapters. Libraries were sequenced with paired ends on an Illumina platform to achieve an average depth of coverage of ∼30×. Alignment of the GRCh37/hg19 genome and human mitochondrial DNA (NC_012920) to the revised Cambridge Reference Sequence (rCRS), variant calling, annotation, and extensive variant filtering were applied through an inhouse bioinformatics pipeline. The copy number variant (CNV) recall was based on Illumina's DRAGEN pipeline. All variants in the GnomAD database with minor allele frequency of <1% and disease-causing variants reported in the Human Gene Mutation Database, ClinVar, and CentoMD were evaluated. During the evaluation, the entire gene region was examined for candidate variants associated with the phenotype in exons and in ±20 intronic regions. All potential patterns were considered for the inheritance pattern. Furthermore, family history and clinical information were used to assess the pathogenicity and probability of the disease-causing effect of the detected variants. Variants were divided into five classes (pathogenic, possibly pathogenic, of uncertain clinical significance, possibly benign, and benign) according to the ACMG guidelines. All potential variants associated with the patient's phenotype were reported. CNVs of uncertain clinical significance were not reported. Variants with heteroplasmy levels of up to 15% were reported as mitochondrial variants. Centogene has established specific quality criteria and validation processes for variants detected by NGS. Variants with poor sequence quality and/or indeterminate zygosity were confirmed by orthogonal methods. As a result, a specificity of >99.9% was guaranteed for all reported variants.

Chromosomal microarray

Technical work was conducted using the Affymetrix Cytoscan 750 K Microarray Platform. Findings that fell within the resolution thresholds (at least 25 markers for decreases of 100 kb and greater and increases of 200 kb and greater in known genomic areas) were analyzed using filters in the “Chromosome Analysis Suite (CHAS)” analysis program. Findings outside the reporting standards; absence of heterozygosity regions; additional mosaic findings; findings below the resolution threshold overlapping with the Online Mendelian Inheritance in Man (OMIM) genes, above-threshold findings without the OMIM genes, the OMIM genes not detected but findings in loci associated with the phenotype; findings that cannot be directly associated with the reason for referral; and results other than the OMIM gene point mutations, small deletions or duplications, low rate mosaicism, and balanced chromosomal changes such as reciprocal and the Robertsonian translocations were not evaluated. Whether the findings were familial or not was determined only through the studies requested by the parents. When required, the segregation analysis of the CNVs detected in the mother and father via array analysis was investigated.

GRCh37/hg19 was used as the reference genome in the analysis. The ClinGen, Decipher, UNIQUE, OrphaNet, DGV, ClinVar, OMIM®, dbSNP, and HPO databases associated with the patient's clinical findings were used to assess the variants. Detected variants were assessed using the patient's clinical findings.

Detected variants were classified according to the ACMG-ClinGen Variant Classification Guidelines. Variants classified as pathogenic, possibly pathogenic, and of uncertain clinical significance associated with the patient's clinical findings were reported. Reported variants were named according to the ISCN 2020 nomenclature ().

Statistical analysis

The demographic and clinical characteristics of the patients evaluated in the study were analyzed via descriptive statistical analyses (number, percentage, mean, standard deviation, etc.). The comparisons of age, clinical findings, perinatal history, demographic characteristics, MRI, and EEG findings according to the genetic analysis results of the patients evaluated in the study were conducted using chi-square analysis. The significance level for all analyses was set at p < 0.05. The conformity of the data to normal distribution was checked with the values of kurtosis and skewness (±1.5). IBM SPSS 22.0 was used to conduct the analyses.

Results

In total, 90 patients were included in the study. Of the 90 patients, 53 (58.9%) were male, and 37 (41.1%) were female. The age of the patients at the time they were seen at our clinic was 0.04–16 years. The mean age of the patients was 3.99 ± 2.71 years, and the median age was 2.71 years. Forty-seven patients were from Turkey; the remaining 43 patients were international (Table 1). EEG abnormalities were detected in a total of 58 (64.4%) patients. EEG was not available from four patients (P68, P77, P79, P81) (Table 2).

Table 1

Patient characteristicsN%
SexFemale3741.1
Male5358.9
CountryTurkey4752.2
Abroad4347.8
Age at admission (years)<23640
2–52932.2
≥52527.8
Age at seizure onset (years)No seizures2831.1
<24853.3
≥21415.6

Demographic characteristics of all patients.

Table 2

EEG findingsEpilepsy onlyEpilepsy+DDDD
N%N%N%
FED555.62852.9310.8
G222.291700.0
FED+G111.159.400.0
SWA00.023.813.5
BS00.023.800.0
Normal111.1713.22071.4

Electroencephalography findings.

BS, burst suppression; DD, developmental delay; EEG, electroencephalography; FED, focal epileptic discharge; G, generalized epileptic discharge; SWA, slow wave activity.

Cranial MRI was abnormal in 44 patients (48.8%). The MRI results were abnormal in 27 (50%) of 54 patients with P/LP variants and clinically relevant VUS variants. MRI was not available in two patients (P68, P82) (Table 3).

Table 3

MRI findingsN*%
Corpus callosum atrophy1516.6
Ventriculomegaly1314.4
White matter changes1213.3
Cerebral atrophy1011.1
Cerebellar atrophy910
Basal ganglia and/or thalamus lesions66.6
Hypoxic ischemic brain injury55.5
Hydrocephalus33.3
Cortical dysplasia22.2
Cerebral venous sinus thrombosis11.1
Intracranial hemorrhage11.1
Hypothalamic hamartoma11.1
Holoprosencephaly11.1
Optic nerve atrophy11.1

Magnetic resonance imaging findings of all patients.

MRI, magnetic resonance imaging.

*Some patients had more than one finding.

In 15 (16.6%) patients, parental consanguinity/origin from same village was reported. There was first-degree cousin marriage between the parents in seven patients (7.8%), and second-degree cousin marriage between the parents in five patients (5.6%). In total, 24 patients (26.7%) had a family history of similar diseases. Thirty-one patients (34.4%) had abnormal perinatal history, such as asphyxia, infection, and hospitalization in the neonatal intensive care unit (Supplementary Table S1).

The patients were classified into three main groups; epilepsy only (n = 9; 10%), DD and epilepsy (n = 53; 58.9%), and only DD (n = 28; 31.1%) (Figure 1).

Figure 1

Molecular findings

Among all patients (n = 90), WES was studied in 58 patients, WGS was studied in 12 patients, and WES and CMA were studied in 20 patients. Trio-WES was performed in four, and trio WGS was performed in two families. Among the 58 patients who had undergone WES, P/LP variants were detected in 20 (34.4%), and clinically relevant VUS was detected in 11 (18.9%). Among the 12 patients who had undergone WGS, P/LP variants were detected in 5 (41.6%), and VUS was detected in 5 (41.6%). Among the 20 patients who had undergone WES and CMA, P/LP variants were detected in 11 (55%), and VUS was detected in 2 (10%) (Figure 2). A molecular diagnosis enabling a potential treatment option was established 12 patients (13.3%), (Table 4).

Figure 2

Table 4

PSGenetic testGeneZygosityIACMGDT
P2MWES, CMAPRRT2 (NM_001256442.1.)HeterozygousADPBFISOXC
P3MWESPRRT2 (NM_145239.3)HeterozygousADPBFISPB
P4MTrio-WGSKCNQ2HeterozygousAD, de novoLPDEEPB, CBZ
P5MWESSCN2A (NM_001040143.2)HeterozygousAD, de novoLPBFISLEV,CBZ
P10MWESFOLR1 (NM_016729.3)HomozygousARVUSNCFTDCZP, TPM, folinic acid
P11MTrio-WESS CN1A (NM_001353958.2)HeterozygousAD, de novoLPDravet SyndromeVPA, TPM
P41MWES CMASCN1A (NM_001202435.1)HeterozygousADLPDravet SyndromeLEV, TPM, PB, STP, FFA
P27FWESSLC2A1 (NM_6516.4)HeterozygousAD, de novoVUSGLUT-1 deficiencyKD, LEV, CLB
P25MWES CMAPNPO (NM_018129.3)HomozygousARPPNPODP5P, VPA
P35MWESALDH7A1 (NM_001201377.2)HomozygousARPEPDTPM, PB, folic acid, L-carnitine, P5P
SLC22A5 (NM_003060.4)HomozygousARLPCDSP
P33FWESKCNQ2 (NM_172107.4)HeterozygousADPDEE 7VGB, CBZ

Molecular diagnosis enabling potential treatment in epilepsy only and epilepsy+DD groups.

BFIS, benign familial infantile epilepsy; CBZ, carbamazepine; CDSP, primary carnitine deficiency; CLB, clobazam; CZP, clonazepam; D, diagnosis; DEE, developmental and epileptic encephalopathy; EPD, pyridoxine-dependent epilepsy; F, female; FFA, fenfluramine; I, inheritance; KD, ketogenic diet; LEV, levetiracetam; M, male; NCTFD, cerebral folate deficiency; OXC, oxcarbazepine; P, patient; PB, phenobarbital; PNPOD, pyridoxine 5ʹ-phosphate oxidasedeficiency; P5P, pyridoxal 5ʹ-phosphate; S, sex; STP, stiripentol; T, treatment; TPM, topiramate; VGB, vigabatrin; VPA, valproic acid.

In our cohort, several genes were frequently implicated in cases of genetic epilepsy and developmental delay, including PRRT2 (n:2), KCNQ2 (n:2), SCN1A (n:2), PLA2G6 (n:2), and KDM6A (n:2). Among molecularly diagnosed cases SCN1A, SCN2A, KCNQ2, SCN8A, and ATP1A2 were identified as channelopathy-related genes, representing 7 patients. While channelopathies were less frequent compared to other conditions in the cohort, their distinct clinical characteristics warrant further attention.

Developmental delay was present in 71.4% of patients with channelopathies, slightly lower than the 93.6% observed in patients with other conditions. This lower prevalence may reflect the primary impact of channelopathies on seizure activity and excitability rather than broader neurodevelopmental processes. Similarly, dysmorphic features were observed in 14.2% of patients with channelopathies, significantly less frequent compared to 36.1% in other genetic conditions, highlighting the subtler phenotypic manifestations of channelopathies.

Abnormal MRI findings were detected in 28.5% of patients with channelopathies, which was markedly lower than the 53.1% observed in patients with other genetic conditions. This suggests that structural brain abnormalities are less commonly associated with channelopathies, reinforcing their functional nature. In contrast, EEG abnormalities were present in 71.4% of patients with channelopathies, slightly higher than the 61.7% seen in other genetic conditions. This reflects the prominent role of electrical disturbances in the pathophysiology of channelopathies.

Table 5 provides a detailed summary of the molecular findings in the diagnosed cohort, further illustrating the differences in clinical and genetic characteristics between channelopathies and other neurodevelopmental conditions.

Table 5

PGMolecular testGene (Transcript)NucleotideProteinZygosityInheritanceVariantACMGNovel/Previously reported variantDiagnosis (#OMIM number)
Epilepsy only
P2MWES CMAPRRT2 (NM_001256442.1.)c.629dupp.Ala211Serfs*14HetADfsPClinvar 31171BFIS, 2 (#605751)
P3MWESPRRT2 (NM_145239.3)c.694dupp.Arg217Profs*8HetADfsPClinvar 65758BFIS, 2 (#605751)
P4MTrio-WGSKCNQ2, EEF1A220q13.33 (62074225_62129854)x1HetAD, de novoCNVLPNovelDEE 7 (#613720), BFNS1 (#121200), DEE 33(#616409)
P5MWESSCN2A (NM_001040143.2)c.379G>Cp.Val127LeuHetAD, de novomisLPNovelDEE, 11(#613721), BFIS, 3(#607745)
P6FWESNBEA (NM_015678.5)c.1370G>Tp.Cys457PheHetADmisVUSNovelNEDEGE (#619157)
Epilepsy and developmental delay
P10MWESFOLR1 (NM_016729.3)c.322G>Ap.Glu108LysHomARmisVUSClinvar 450195NCFTD (#613068)
P11MTrio-WESSCN1A (NM_001353958.2)c.5341T>Ap.Tyr1781AsnHetAD, de novomisLPC Depienne et al. ()Dravet syndrome (#607208)
P12MWESKNL1 (NM_170589.4)c.6700C>Tp.His2234TyrHomARmisVUSNovelMCPH 4 (#604321)
P13FWESATP1A2 (NM_000702.4)c.2429G>Ap.Gly810AspHetADmisVUSNovelDEE 98 (#619605)
P14FWGS22q11 del.22q11 (18877539-21681527)x1HetADCNVPDecipher 472470 Akgun-Dogan et al. ()22q11 deletion syndrome (#611867)
P15MWESASH1l (NM_018489.3)c.7189C>Tp.Arg2397*HetADnsLPHGMD CM1513932IDD, AD 52 (#617796)
P16MTrio-WESPCK1 (NM_002591.4)c.524G>Ap.Arg175GlnHomARmisVUSClinvar 1368065PCKDC (#261680)
P17MWES CMAABCC8 (NM_000352.6)c.4432G>Ap.Gly1478ArgHetAD, paternalmisLPClinvar 434045HHF 1 (#256450)
P18MWES CMAKAT6A (NM_006766.4)c.5828T>Cp.Met1943ThrHetADmisVUSNovelMRD 32 (#616268)
MT-ATP8 (NC_012920.1)c.98A>G*p.Tyr33CysHomoplasmicmtmisVUSNovelMC5DM2 (#516070)
P19FWGS2q37del.2q37.3 (237531981_242782258)x1HetADCNVPDecipher 324169 Akgun-Dogan et al. ()2q37 microdeletion syndrome (#600430)
3q29 dup.3q29.32q29 (176456169_197851444)x3HetADCNVPS. Goobie et al. () 2008 Akgun-Dogan et al. ()3q29 microduplication syndrome (#611936)
P20FWESGNB1 (NM_002074.5)c.233A>Gp.Lys78ArgHetADmisPClinvar 224714IDD, AD 42 (#616973)
P21MWESUBA5 (NM_024818.6)c.863T>Cp.Met288ThrHomARmisVUSClinvar 1429996DEE 44(#617132)
P22FWESPACS2 (NM_001100913.3)c.625G>Ap.Glu209LysHetADmisPClinvar 495141DEE 66 (#618067)
P23MWESPPT1 (NM_000310.4)c.538dupCp.Leu180Profs*9HomARfsPClinvar 56201NCLs 1(#256730)
P24MWESADAM22 (NM_021723.5)c.1312C>Ap.Pro438ThrHomARmisVUSMarieke M. van der Knoop et al. ()DEE 61 (#617933)
P25MWES CMAPNPO (NM_018129.3)c.674G>Ap.Arg225HisHomARmisPClinvar 223153PNPOD (#610090)
P27FWESSLC2A1 (NM_6516.4)c.275+1delp.?HetAD, de novoSVUSClinvar 1048756GLUT-1 deficiency syndrome (#606777)
P28MWESPLA2G6 (NM_003560.4)22q13.1 (38484755_38511688)x0HomARCNVPNeil V Morgan et al. ()INAD 1 (#26600)
P29FWES CMAKDM6A (NM_001291415.2)c.4180C>Tp.Leu1394PheHetXLDmisVUSNovelKabuki syndrome (#300867)
P30MWGSSPTBN4 (NM_020971.2)c.455T>Gp.Leu152ArgC/HARmisVUSAkgun-Dogan et al. ()NEDHND (#617519)
SPTBN4 (NM_020971.2)c.7149delp.Pro2384Argfs*60C/HARfsVUSClinvar 1709124 Akgun-Doganet al. ()NEDHND (#617519)
P31FWES CMAWWOX (NM_016373.4)c.716T>Gp.Leu239ArgHomARmisLPClinvar 871669DEE 28 (#616211)
P32FWGSMACF1 (NM_012090.5)c.3512T>Cp.Leu1171ProHetADmisVUSNovelLIS 9 (#618325)
P33FWESKCNQ2 (NM_172107.4)c.793G>Ap.Ala265ThrHetADmisPMilh et al. ()DEE 7 (#613720), BFNS1 (#121200)
P34FTrio-WESGFAP (NM_002055.5)c.1217G>Tp.Arg406MetHetAD, de novomisLPNovelAlexander Disease(#203450)
P35MWESALDH7A1 (NM_001201377.2)c.988C>Tp.Arg330*HomARnsPClinvar 1457545EPD (#266100)
SLC22A5 (NM_003060.4)c.1088T>Cp.Leu363ProHomARmisLPClinvar 25408CDSP (#212140)
P36FWESPLA2G6 (NM_003560.4)c.753dupp.Asn252Glnfs*130HomARfsLPSallevelt et al. ()INAD 1 (#256600)
P37FWGSNAGA (NM_000262.2)c.973G>Ap.Glu325LysC/HAR, paternalmisPClinvar 18162 Akgun-Dogan et al. ()NAGA deficiency type 3 (#609241)
NAGA (NM_000262.2)c.103C>Tp.Arg35CysC/HAR, de novomisVUSClinvar 429497 Akgun-Dogan et al. ()NAGA deficiency type 3 (#609241)
P39MTrio-WGSLSS (NM_001001438.2)c.1970T>Gp.Met657ArgC/HARmisVUSNovelAPMR 4 (#618840)
LSS (NM_001001438.2)c.2123G>Tp.Arg708MetC/HARmisVUSNovelAPMR 4 (#618840)
TUBB (NM_001293212.1)c.334T>Gp.Phe112ValHetAD,de novomisVUSNovelCDCBM 6 (#615771)
P40MWESCUX1 (NM_001913.5)c.1733_1734dupGCp.Leu579Alafs*18HetADfsLPNovelGDDI (#618330)
P41MWES CMASCN1A (NM_001202435.1)c.1070 delp. Asn357llefs*22HetADfsLPNovelDravet syndrome (#607208)
P44FWESSETD1A (NM_014712.3)c.2642G>Tp.Gly881ValHetAD,maternalmisVUSNovelEPEDD (#618832)
Developmental delay
P63MWESBCL11A (NM_022893.4)c.1486G>Tp.Glu496TerHetAD,de novonsPNovelDias-Logan Syndrome (#617101)
SCN8A (NM_014191.4)c.1743_1744delTGp.Asn581LysfsTer5HetAD,de novofsPNovelBFIS 5(#617080), CIAT (#614306), DEE 13(#614558)
P64MWGSTHOC2 (NM_001081550.1)c.521A>Gp.Lys174ArgHemizygousXLR, maternalmisVUSNovelXLID 12 (#300957)
P65MWES CMA3q13.31 mikrodel.3q12.1q13.32 (98920626_118556822)x1HetADCNVPMolin et al. ()3q13.31 deletion syndrome (#615433)
P66MWESBCL11B (NM_138576.4)c.2438_2459delp.Val813Alafs*24HetAD,de novofsLPGoos et al. ()IDDSFTA (#606558)
P67FWGSERCC6 (NM_000124.3)c.543+1G>Tp.?HomARSLPClinvar 974894 Akgun-Dogan et al. ()Cockayne syndrome (#133540)
P68MTrio-WESECEL1 (NM_004826.4)c.1507-9G>Ap.?C/HARlSLPClinvar 210907Arthrogryposis, distal type (#615065)
ECEL1 (NM_004826.4)c.1630C>Tp.Arg544CysC/HARmisVUSClinvar 816801Arthrogryposis, distal type (#615065)
P69FWES CMAKDM6AXp11.3 (44733823_44837306)x1HetXLDCNVPClinvar 564836Kabuki syndrome (#300867)
P70MWES CMATCF 20 (NM_005650.4)c.5452G>Tp.Glu1818*HetAD, de novonsPNovelDDVIBA (#618430)
P71FWESDHCR7 (NM_001360.3)c.452G>Ap.Trp151*C/HARnsPClinvar 21273Smith-Lemli-Opitz Sendromu (#270400)
DHCR7 (NM_001360.3)c.1295A>Gp.Tyr432CysC/HARmisPClinvar 1076651Smith-Lemli-Opitz Syndrome (#270400)
P72FWESGCDH (NM_000159.4)c.383G>Ap.Arg128GlnHetARmisLPClinvar 189063Glutaric acidemia type 1(#231670)
GCDH (NM_000159.4)c.1282A>Gp.lle428ValHetARmisVUSNovelGlutaric acidemia type 1(#231670)
P73MWESDARS1 (NM_001349.4)c.606C>Gp. lle202MetHomARmisVUSNovelHBSL (#615281)
ERLIN2 (NM_007171.2)c.78A>Gp.Ile26MetHetADmisVUSNovelSPG 18 (#611225)
P74FWES CMAZEB2 (NM_014795.4)c.1091C>Gp.Ser364*HetADnsLPNovelMowat-Wilson Syndrome (#235730)
P75MWES CMASOX512p12.1 (23805762_23967495)x1 delHetADCNVPLamb et al. ()Lamb-Shaffer Syndrome (#616803)
P76MWESPANK2 (NM_153638.3)c.676G>Ap.Val226IleHomARmisLPAkcakaya et al. ()NBIA 1 (#234200)
P77FWESPIK3R2 (NM_005027.4)c.1117G>Ap.Gly373ArgMozaicAD,de novomisPClinvar 39808MPPH 1 (#603387)
P78MWESSLC1A4 (NM_003038.5)c.911G>Ap.Gly304AspHomARmisVUSNovelSPATCCM (#616657)
P79MWGSLAMB1 (NM_002291.2)c.4756A>T,p.Thr1586SerC/HARmisVUSNovelLIS 5 (#615191)
LAMB1 (NM_002291.2)c.4859T>Cp.Ile1620ThrC/HARmisVUSClinvar 999947 Akgun-Dogan et al. ()LIS 5 (#615191)
P80FWES CMAKIF1A2q37.3 (240058543_242782258)x1 delHetADCNVPDecipher 411683, 277875Nescav Syndrome (#614255)

Molecular findings in the diagnosed cohort.

Significant variants detected in the patients in our study are listed along with their characteristics. Variants reported in the literature are given with Decipher, HGMD, and ClinVar numbers. Variants that had not been reported previously in the literature are called new variants. These genetic variants’ classification detected according to by June 2024.

ACMG, American College of Medical Genetics and Genomics; APMR4, alopecia-intellectual disability syndrome 4; BFIS, benign familial infantile seizures; BFNS, benign familial neonatal seizures; CDCBM 6, cortical dysplasia, complex, with other brain malformations 6; CDSP, carnitine deficiency, systemic primary; CIAT, cognitive impairment with or without cerebellar ataxia; CNV, copy number variant; DDVIBA, developmental delay with variable intellectual impairment and behavioral abnormalities; C/H, compound heterozygous; DEE, developmental and epileptic encephalopathy; EPD, epilepsy, pyridoxine-dependent; EPEDD, epilepsy, early-onset, with or without developmental delay; F, female; fs, frameshift variant, G, gender; GDDI, global developmental delay with or without impaired intellectual development; GLUT1DS1, GLUT-1 deficiency syndrome; HBSL, hypomyelination with brainstem and spinal cord involvement and leg spasticity; Het, heterozygous; Hom, homozygous; HGMD, human gene mutation database; HHF, hyperinsulinemic hypoglycemia; IDD, AD, intellectual developmental disorder, autosomal dominant; IDDSELD: epilepsi ve konuşma geriliği ile birlikte zihinsel yetersizlik; IDDSFTA, intellectual developmental disorder with dysmorphic facies, speech delay, and T-cell abnormalities; INAD, infantile neuroaxonal dystrophy; LIS 9, lissencephaly 9 with complex brainstem malformation; M, male; MCPH 4, primary microcephaly 4; MC5DM2, mitochondrial complex 5 (ATP synthase) deficiency, mitochondrial type 2; mis, missense variant; MPPH1, megalencephaly-polymicrogyria-polydactyly-hydrocephalus syndrome 1; MRD 32, mental retardation, autosomal dominant; mt, mitochondrial; NBIA1, neurodegeneration with brain iron accumulation 1; NCFTD, neurodegeneration due to cerebral folate transport deficiency; NCLs, neuronal ceroid lipofuscinosis; NEDEGE, neurodevelopmental disorder with or without early-onset generalized epilepsy; NEDHND, neurodevelopmental disorder with hypotonia, neuropathy, and deafness; ns, nonsense variant; OMIM, online Mendelian inheritance in man; PCKDC, phosphoenolpyruvate carboxykinase deficiency, cytosolic; PNPOD, pyridoxamine 5′-phosphate oxidase deficiency; SPATCCM, spastic tetraplegia, thin corpus callosum, and progressive microcephaly; SPG 18, spastic paraplegia 18; S, splice-site variant; XLID 12, intellectual developmental disorder, X-linked 12.

Dysmorphic findings were significantly more frequent in the molecularly diagnosed group (P/LP/VUS variants), whereas other features, such as MRI abnormalities, neurobehavioral findings, and ocular findings, did not show a statistically significant difference (Table 6).

Table 6

Clinical findingsGenetic resultsP
P/LP/VUS variants
n%n%
MRI findings2749.11751.40.828
Spasticity1629.61130.60.925
Hypotonia2037.0719.40.074
Dystonia47.412.80.348
Ataxia00.025.60.080
ADHD814.838.30.358
Autism11.912.80.720
Microcephaly814.8513.90.903
Macrocephaly47.400.00.095
Dysmorphism1833.325.60.002
Strabismus47.425.60.730
Nystagmus00.012.80.218
Optic atrophy23.712.80.811
Coloboma11.900.00.412
Cataract11.912.80.770
Visual impairment11.900.00.412
Microphthalmia00.012.80.218
Swallowing difficulty713.0411.10.793
Hearing loss23.712.80.811
Status epilepticus23.712.80.811
Cleft palate lip11.912.80.770
Hypothyroidism35.625.60.999
Hydrocephalus23.712.80.811
Sialorrhea47.412.80.348

Comparison of clinical characteristics according to patients' genetic results.

ADHD, Attention deficit and hyperactivity disorder; MRI, Magnetic resonance imaging.

Discussion

Several epilepsy phenotypes exhibit chromosomal imbalances or alterations in genes encoding ion channels and transcription factors. Genetic testing is crucial for elucidating the underlying cause in patients with epilepsy and DD (). Despite advanced diagnostic technologies nearly half of the affected individuals remain undiagnosed.

In the literature, the rate of diagnosis of DD and epilepsy by utilizing WES is reported to be between 25% and 55% (). Hiraide et al., have reported that WES, including CNV analysis, resulted in a diagnostic rate of 53.5% (). In our study, the percentage of LP/P variants with WES was 34.4%. Clinically relevant VUS were observed in an additional 18.9%, providing a similar diagnostic rate (53.3%) despite absence of CNV analysis. Similarly, the diagnostic rate of WGS in DD and epilepsy cohorts is reported to be 21%–50% (). Among the 12 patients who had undergone WGS, 5 (41.6%) had P/LP variants, and 5 (41.6%) had clinically relevant VUS. WGS was performed in patients without prior testing. Two out 5 P/LP WGS results showed CNV as the underlying cause which would have been missed by WES alone. WGS shortened the diagnostic duration in these patients as first-tier genetic test.

The findings of our study indicate that the most frequently observed genes are consistent with those reported in the literature (). Furthermore, although a direct comparison is not feasible due to the significant differences in the numbers, channelopathy-related genes stood out prominently in our study. Similarly, in the literature, channelopathy-related genes have also been highlighted as a significant group, in line with our data. Compared to other genetic etiologies underlying epileptic encephalopathies, channelopathies are less frequently associated with dysmorphic features and abnormal MRI findings, while EEG abnormalities are more prominently observed. This highlights the critical importance of genetic analysis in establishing an accurate diagnosis. Furthermore, genetic findings play a pivotal role in informing and tailoring treatment strategies (, ).

In this study, a molecular genetic diagnosis with available treatment was established in 13.3% of the patients. Only one patient with an actionable variant was left untreated due to absence of clinical seizures. Mahler et al. reported a treatment change rate of 6%. The higher rate of specific treatment changes in our study may be attributed to our patients having a more heterogeneous phenotype, as well as the emergence of new treatment options involving genes related to ion channels (). By re-evaluating the previously negative WES data of patient P27, a de novo VUS variant in the SLC2A1 gene was identified; this was followed by segregation analysis, which led to a diagnosis of GLUT-1 deficiency (). Seizures were controlled using a ketogenic diet, demonstrating the importance of data reanalysis in treatment modifications.

Trio WES analysis was conducted in 4 families due to economic burden. Segregation analysis was performed following initial WES analysis. 20.7% of the diagnosed individuals showed de-novo monoallelic variants. Brunet et al. reported a much higher frequency of de- novo variants (40.3%) following trio WES in their study. The lower percentage of de-novo monoallelic variants in our cohort reflects the higher rate of consanguinity leading to a higher incidence of biallelic disorders (22.2%) in comparison to 16% reported by Brunet et al. ().

Samati et al. reported abnormal MRI findings in one-third of individuals with epilepsy in a recent review. The most frequently observed abnormalities were encephalomalacia related to chronic infarcts, cerebral atrophy, disorders of neuronal migration, periventricular leukomalacia. Accompanying EEG finding was mainly focal discharges (78%) (). In our study, corpus callosum atrophy (16.6%), ventriculomegaly (14.4%), and white matter changes (13.3%) were present in 48.8% of our patients and focal discharges were most common (40%) as well. Our cohort mainly consists of unexplained patients without evident focal lesions. The relatively lower percentage of focal discharges may be related to the lower number of focal lesions on MRI in our patient group who underwent genetic testing.

The main limitation of our study is the retrospective nature and heterogeneity of the genetic tests. Some patients had only WES, some had WES and CMA, and some had WGS. Therefore, the number of patients having each test was limited. However, our study adds novel genetic variants that have not been described previously in the literature. All patients received pre-test and post-test genetic counseling from two pediatric geneticists who also contributed to deep-phenotyping and reverse phenotyping of the diagnosed cohort. Patients with negative genetic results will be recalled for reanalysis yearly. Genetic counseling is provided by clinical geneticists in Turkey. Families with a diagnosis were counseled regarding recurrence risks and future reproductive options. This is the largest single center study with a multidisciplinary approach from our country.

In conclusión, NGS helps to determine the underlying molecular etiology of neurogenetic conditions. Moreover, an increasing percentage of patients receive personalized treatment after molecular diagnosis. WES and WGS increase the diagnosis rate and shorten the diagnostic odyssey for families when applied as first-tier genetic tests. Determining the genetic diagnosis of patients with epilepsy and/or DD using NGS is crucial in terms of an accurate prognosis, potential treatment methods, and a multidisciplinary approach involving molecular and clinical genetic expertise.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Scientific Research Ethics Committee of Acibadem University School of Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin.

Author contributions

HK: Data curation, Formal Analysis, Investigation, Writing – original draft, Methodology, Project administration. OA-D: Conceptualization, Data curation, Writing – review & editing, Supervision, Validation. AY: Methodology, Writing – review & editing, Data curation. YA: Data curation, Writing – review & editing, Supervision, Conceptualization, Validation. UI: Data curation, Investigation, Project administration, Supervision, Writing – review & editing, Conceptualization.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The 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.

Supplementary material

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

References

Summary

Keywords

epilepsy, developmental delay, next-generation sequencing, pediatric neurology, pediatric geneticist

Citation

Kava H, Akgun-Dogan O, Yesilyurt A, Alanay Y and Isik U (2025) Evaluation of the etiology of epilepsy and/or developmental delay in children via next-generation sequencing: a single-center experience. Front. Pediatr. 13:1471965. doi: 10.3389/fped.2025.1471965

Received

28 July 2024

Accepted

28 January 2025

Published

27 February 2025

Volume

13 - 2025

Edited by

Piero Pavone, University of Catania, Italy

Reviewed by

Raffaele Falsaperla, Policlinico San Marco, Italy

Eulàlia Turon-Viñas, Universitat Autònoma de Barcelona, Spain

Updates

Copyright

*Correspondence: Yasemin Alanay

† These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics