Skip to main content

Identification of ferroptosis-related genes in acute phase of temporal lobe epilepsy based on bioinformatic analysis



Epilepsy is a prevalent neurological disorder, and while its precise mechanism remains elusive, a connection to ferroptosis has been established. This study investigates the potential clinical diagnostic significance of ferroptosis-related genes (FRGs) during the acute phase of temporal lobe epilepsy.


To identify differentially expressed genes (DEGs), we accessed data from the GEO database and performed an intersection analysis with the FerrDB database to pinpoint FRGs. A protein-protein interaction (PPI) network was constructed. To assess the diagnostic utility of the discovered feature genes for the disease, ROC curve analysis was conducted. Subsequently, qRT-PCR was employed to validate the expression levels of these feature genes.


This study identified a total of 25 FRGs. PPI network analysis revealed six feature genes: IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN. ROC curve analysis demonstrated that the combination of these six feature genes exhibited the highest diagnostic potential. qRT-PCR validation confirmed the expression of these feature genes.


We have identified six feature genes (IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN) strongly associated with ferroptosis in epilepsy, suggesting their potential as biomarkers for the diagnosis of temporal lobe epilepsy.

Peer Review reports


Epilepsy, as a prevalent neurological disorder, currently affects approximately 50 million people worldwide [1]. Approximately 30% of epileptic patients cannot manage their symptoms with medication, and most of them are not suitable candidates for epilepsy surgery [2, 3]. Present antiepileptic drugs can only prevent seizures but do not impact the underlying pathological processes. Temporal lobe epilepsy (TLE), as a distinct subtype of epilepsy, is characterized by various etiological factors and complex clinical features, making diagnosis and treatment challenging [4]. Therefore, the identification of early diagnostic biomarkers in the acute phase of TLE could greatly assist in clinical diagnosis.

Ferroptosis, a novel form of regulatory cell death induced by excessive iron ion accumulation, leads to the accumulation of lipid reactive oxygen species and subsequent cell demise [5]. Recent research suggests that iron metabolism abnormalities are widespread in epilepsy [6]. Furthermore, studies indicate that the imbalance and accumulation of free iron can lead to oxidative stress-dependent damage to central nervous system neurons, resulting in epilepsy [7]. Emerging research also suggests that inhibiting ferroptosis can mitigate neuronal damage induced by epilepsy [8]. Therefore, we postulate that ferroptosis may provide valuable new insights into the pathogenesis of epilepsy.

Bioinformatics analysis has become a powerful tool in disease research, aiding in the identification of potential biomarkers and therapeutic targets. In recent years, some studies have used bioinformatics approaches to gain deeper insights into various diseases. Zhao et al. employed bioinformatics methods to explore kainic-induced TLE based on integrated bioinformatics analysis [9], while other research has also explored potential key genes of epilepsy. However, there has been limited exploration of diagnostic biomarkers related to ferroptosis in epilepsy through bioinformatics analysis.

The focus of this study is to identify diagnostic biomarkers associated with ferroptosis in epilepsy and elucidate the molecular mechanisms driving epileptogenic progression. These investigations enhance our understanding of the role of ferroptosis in the pathogenesis of epilepsy, laying a new theoretical foundation for future epilepsy treatments.

Materials and methods

Microarray data Processing

We obtained microarray expression profile datasets related to acute-phase TLE from the NCBI GEO database ( Initially, we downloaded the GSE88992 dataset, which comprises hippocampal tissues from 8 pilocarpine-induced acute-phase TLE mice and 9 control hippocampal tissues. Additionally, we obtained the GSE49030, GSE79129, and GSE143272 datasets for validation purposes. GSE49030 includes 12 acute-phase TLE mouse hippocampal tissues and 6 control hippocampal tissues. GSE79129 contains 3 acute-phase TLE mouse dentate gyrus tissues and 3 control tissues. GSE143272 comprises peripheral blood samples from 34 newly diagnosed epilepsy patients without drug treatment and 51 healthy individuals. The microarray data were normalized using the normalize between arrays function of the limma package (version 3.52.4) in R software (version 4.2.1) [10].

All of the mentioned data are freely available for download from the GEO database. This study does not require ethical approval or patient informed consent.

Data Processing and Selection of differentially expressed genes

We utilized the limma package in the R software to analyze the expression profile of the GSE88992 dataset. In order to narrow down the range of differentially expressed genes (DEGs), we applied the criteria of |log2FC| > 1 and adjusted P-value < 0.05. Additionally, we employed the ggplot2 package to create a volcano plot, enhancing the visualization of the DEGs.

Acquisition of Ferroptosis-related genes

The FerrDb database ( supplied us with a dataset comprising 259 genes associated with ferroptosis [11]. These ferroptosis-related genes were obtained from the database based on factors related to ferroptosis, including drivers, inhibitors, and markers.

To identify Ferroptosis-Related Genes (FRGs), we performed an intersection analysis between the differentially expressed genes and the ferroptosis-related genes. The results were visualized using a Venn diagram.

GO analysis and KEGG Pathway Enrichment analysis of FRGs

In order to obtain deeper insights into the biological functions and pathways of the FRGs, we conducted enrichment analysis and data visualization using R software and the ggplot2 package. Specifically, we carried out GO enrichment analysis to uncover the potential roles of FRGs in Biological Process (BP), Cellular Component (CC), and Molecular Function (MF). Additionally, KEGG pathway enrichment analysis was employed to identify the pathways linked to these differentially expressed genes [12,13,14]. The filtering criteria were established as follows: P-Value < 0.05.

Construction of PPI Network and Hub Gene Selection

All the filtered FRGs were uploaded to the STRING 11.5 database ( The STRING database is a valuable resource that encompasses both known and predicted direct physical and indirect functional interactions among proteins/genes. We utilized this database to analyze the Protein-Protein Interaction (PPI) network among the set of FRGs. In our analysis, we considered interaction scores greater than 0.7 as representing medium-confidence interactions. For the sake of clarity, we concealed disconnected nodes in the PPI network.

Furthermore, we made use of the Molecular Complex Detection (MCDOE) plugin within the Cytoscape software (version 3.9.1) to pinpoint the central node genes in the PPI network and calculate significant PPI network modules. The MCODE plugin is a widely used plugin in Cytoscape. It is designed to identify key genes within large networks based on relationships between edges and nodes, facilitating subsequent analyses.

Acquisition of biomarkers for Biomedical diagnosis

We assessed the previously identified central node genes using the Area Under the Curve (AUC) analysis. Consequently, we computed AUC values separately for the six genes and the logistic regression model to evaluate the accuracy of the diagnostic model. Finally, we calculated the AUC for the ROC curves and generated the curves using the pROC package (version 1.18.4) and ggplot2 package (version 3.4.3), both available in R.

We conducted this analysis using three separate independent validation datasets: GSE49030, GSE79129, and GSE143272, to validate the accuracy of Hub genes as diagnostic molecules.

Extraction of Hippocampus from epileptic mice

Adult male C57BL/6 mice (8–10 weeks old and weighing 20-25 g) were obtained from the Experimental Animal Center of Wenzhou Medical University. The mice were housed under specific laboratory conditions with a 12-hour light/12-hour dark cycle and a temperature of 24 ± 2 °C. They had free access to food and water. After intraperitoneal injection of pentobarbital sodium (50 mg/kg), the mice were fixed on a stereotaxic instrument. Following the coordinates previously described in the literature for the hippocampal region (AP-2.0 mm, ML-1.5 mm, V-2.0 mm), KA was slowly injected into the hippocampus while the needle was held in place for 5 min to prevent drug reflux [15]. The control group was injected with an equal volume of physiological saline. 24 h after KA injection, under anesthesia with intraperitoneal injection of pentobarbital sodium, the hippocampi were obtained from both control and KA-treated mice by intraperitoneal injection. All experimental procedures were reviewed and approved by the Ethics Committee of Wenzhou Medical University, and they were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Animals.

RNA extraction and qRT-PCR

Samples from Hippocampus were collected for qRT-PCR analysis. Total mRNA was extracted using Trizol reagent (Takara Bio, Tokyo, Japan), and cDNA was synthesized from 1 µg RNA using PrimeScript RT kit (Takara Bio, Tokyo, Japan). qRT-PCR analysis was conducted on an Opticon 2 Real-Time PCR Detection System (Bio-Rad) using SYBR®Green PCR Master Mix (Takara Bio, Tokyo, Japan) and corresponding primers.

PCR amplifcation was performed with a program of 95 °C for 5 min, followed by 40 cycles of 95 °C for 1 min, 62 °C for 30 s, and 72 °C for 30 s. The mRNA expression level of the target gene was normalized to that of the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase level (GAPDH), and calculated by the 2−ΔΔCTmethod. Results are expressed as the mean ± SEM of replicate samples from three independent experiments.

Statistical analysis

All data were analyzed using Prism 8.0.1 software (GraphPad, San Diego, CA, USA) in a blinded manner. In vivo experiments were presented as mean ± SD, while in vitro experiments were presented as mean ± SEM. Statistical analyses were performed using one-way ANOVA for multiple comparisons or Student’s t-test (and nonparametric tests). A significance level of p < 0.05 was considered statistically significant.


Differential Gene Selection

Using R software, we performed online standardization processing of the GSE88992 gene chip dataset. The resulting Differentially Expressed Genes (DEGs) were visually represented using volcano plots, as depicted in Fig. 1. Employing the filtering criteria of adjusted p-value < 0.05 and |log2FC| > 1, we identified a total of 456 DEGs in GSE88992. Among these, 366 genes were upregulated, and 90 genes were downregulated. The list of differentially expressed genes obtained from the differential analysis can be found in Supplementary Tables S1-4 for specific details.

Fig. 1
figure 1

Ferroptosis-Related DEGs of epilepsy in GSE88992. (A) Volcano plot of all DEGs. (B) Venn diagram of DEGs and ferroptosis-related genes. (C) Bubble Chart of GO enrichment analysis. (D) Bubble Chart of KEGG enrichment analysis

Identification of ferroptosis-related genes (FRGs)

We acquired a gene set from the FerrDb database, which included 259 genes. Cross-referencing this gene set with the DEGs from GSE88992 allowed us to identify 25 FRGs. These FRGs were visually represented in a Venn diagram, as illustrated in Fig. 1. Additionally, we constructed a heatmap to illustrate the differential expression patterns of these 25 Ferroptosis-Related Genes (FRGs) across the four datasets. The details can be observed in Fig. 2. The list of the 25 intersected genes can be found in Supplementary Table S5 for specific details. The heatmap was generated using the pheatmap package (version 1.0.12) in R software.

Fig. 2
figure 2

The heat map of the dataset demonstrated distinguished features between epilepsy samples and control samples. Each column represents a sample, while each row represents a gene. The expression levels of genes, ranging from high to low, are depicted by colors from red to blue, respectively. At the top of the heatmap, blue corresponds to the control group, and red represents the epilepsy group

Enrichment Analysis

We performed enrichment analysis of GO/KEGG pathways using the clusterProfiler package (version 4.8.3) in R software. The results of GO enrichment and KEGG enrichment were visualized in bubble plots. Significance was determined with an adjusted p-value < 0.05. According to the results of GO analysis, the differentially expressed FRGs were predominantly involved in cellular responses to external stimuli, responses to nutrient levels, and cellular responses to chemical stress. Based on the results of KEGG analysis, significant enrichment was observed in pathways such as the lipid and TNF signaling pathway, IL-17 signaling pathway, and HIF-1 signaling pathway among the differentially expressed FRGs.

Construction of protein-protein Interaction (PPI) Network

The STRING database was used to construct the PPI network. Interactions with a minimum interaction score > 0.7 between differentially expressed individuals were visualized in the network. The PPI network of 25 nodes and 38 edges was further analyzed using Cytoscape 3.9.1 software. We calculated the Maximum Clique Centrality (MCC) scores using the CytoHubba plugin and selected the six genes as hub genes, including IL6, PTGS2, HMOX1, NFE2L2, TLR4 and JUN. The PPI network and the key PPI network modules as shown in Fig. 3.

Fig. 3
figure 3

(A) PPI network of ferroptosis-related DEGs. (B) PPI network diagram obtained by MCODE plugin in cytoscape

Validation of feature genes

In the GSE88992 dataset, we validated the diagnostic efficacy of the feature genes by plotting ROC curves and comparing the AUC values to assess their diagnostic value. Additionally, we developed a novel diagnostic model based on the six feature genes using logistic regression. The logistic regression model exhibited an AUC value of 1.0000, indicating its high effectiveness in predicting the occurrence of epilepsy.

Furthermore, we utilized three independent datasets, GSE49030, GSE79129, and GSE143272, to validate the accuracy of the diagnostic model. We calculated AUC scores for both the feature genes and the model in the same manner. In the GSE49030 and GSE79129 datasets, the model’s AUC values were both 1.0000. In GSE143272, the model achieved an AUC value of 0.8397. Specific diagnostic performance of individual genes can be seen in the Fig. 4.

Fig. 4
figure 4

The training set of GSE88992’s feature genes (A) and model (E) ROC curve data. The validation set of GSE49030’s feature genes (B) and model (F) ROC curves. The validation set of GSE79129’s feature genes (C) and model (G) ROC curves. The validation set of GSE143272’s feature genes (D) and model (H) ROC curves. Different colored lines in the ROC curves represent distinct genes

In summary, the diagnostic model we constructed can effectively distinguish epilepsy patients from normal individuals, facilitating the diagnosis and screening of TLE during the acute phase.


To validate the expression levels of the aforementioned feature genes in the hippocampal tissue of mice with epilepsy induced by KA after 1 day, qRT-PCR was performed. Compared to the hippocampal tissue of normal mice, the expression levels of IL6, PTGS2, HMOX1, NFE2L2, TLR4 and JUN were significantly increased in the hippocampal tissue of mice with epilepsy induced by KA after 1 day, as shown in Fig. 5.

Fig. 5
figure 5

(A) mRNA expression of IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN in GSE88992. (B) mRNA expression of IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN in GSE49030. (C) mRNA expression of IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN in GSE79129. (D) mRNA expression of PTGS2, HMOX1, NFE2L2, TLR4, and JUN in GSE143272. (E) mRNA level of IL6, PTGS2, HMOX1, NFE2L2, TLR4, and JUN between sham group and kainic- induced epilepsy mouse group. The bar chart presented as mean ± standard deviation with *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 for all statistics


In our study, we aimed to explore novel biomarkers for epilepsy from the perspective of ferroptosis. In this research, we identified six characteristic genes associated with ferroptosis in epilepsy: NFE2L2, PTGS2, IL6, JUN, HMOX1, and TLR4. These characteristic genes have enriched our understanding of the relationship between ferroptosis and the pathogenesis of epilepsy. Furthermore, we found that a diagnostic model composed of these six characteristic genes exhibited excellent diagnostic performance in both animal models hippocampal tissue and human peripheral blood samples.

Although the exact mechanisms underlying epilepsy are not fully understood, some studies suggest a close association between epilepsy pathogenesis and ferroptosis. Inhibition of ferroptosis has been shown to alleviate neuronal damage and achieve neuroprotection in epilepsy [16]. Furthermore, these six genes also play a crucial role in the development of epilepsy.

NFE2L2, as the gene encoding Nrf2, is an important transcription factor that plays a crucial role in cellular antioxidant stress response and cell protection mechanisms [17]. Its main function is to regulate cellular response to oxidative stress. When cells are exposed to oxidative stress or other toxic stimuli, NFE2L2 is activated and translocates into the cell nucleus, where it binds to antioxidant response elements and initiates the transcription of a series of antioxidant stress genes, including antioxidant enzymes, detoxifying enzymes, and other cell-protective proteins [18, 19]. The expression of these genes helps cells eliminate free radicals, reduce oxidative damage, and promote cellular repair [20]. Additionally, abnormal expression or dysfunction of NFE2L2 has been implicated in the development of various diseases, including cancer, neurodegenerative diseases (such as Alzheimer’s disease and Parkinson’s disease) [21], cardiovascular diseases [22], and inflammatory disorders [23]. Furthermore, NFE2L2 also plays a role in ferroptosis. SLC7A11, as a cystine/glutamate exchange transporter, is involved in the transport of cystine and glutamate, regulating cellular glutathione (GSH) levels [24]. The expression of SLC7A11 is regulated by Nrf2, and its expression increases when Nrf2 activity is enhanced. SLC7A11 protects cells from oxidative stress damage by promoting cystine synthesis and regulating intracellular glutamate concentration. Moreover, NFE2L2/Nrf2 deficiency leads to ferritin accumulation in autophagosomes, elevated unstable intracellular iron pools, and increased sensitivity to ferroptosis [25]. Natural substances such as quercetin have been shown to alleviate neuronal death caused by seizures and preserve cognitive function by inhibiting the Nrf2-mediated ferroptosis pathway [26]. HMOX1, also known as Heme Oxygenase-1, is an enzyme with a primary function of catalyzing the release of heme molecules from hemoglobin. It then breaks down heme into carbon monoxide, iron ions, and biliverdin. In recent years, it has been extensively studied as a downstream target of Nrf2 [27, 28]. There is also growing research indicating the positive role of Nrf2/HO-1 in epilepsy protection [28,29,30]. However, there is relatively limited research regarding its connection with ferroptosis in epilepsy.

A crucial enzyme in the prostaglandin production pathway, PTGS2, referred to as cyclooxygenase-2, catalyzes the transformation of arachidonic acid into prostaglandins [31]. Under normal physiological conditions, PTGS2 expression is typically low, but it is significantly increased during inflammatory reactions, tumorigenesis, and other disease processes [32, 33]. Additionally, in certain situations, increased expression and activity of PTGS2 have been associated with the accumulation of intracellular iron ions and the occurrence of ferroptosis. Specifically, PTGS2 can induce the release and accumulation of intracellular free iron ions by producing a series of prostaglandins and other pro-inflammatory cytokines, such as PGE2. The accumulation of free iron ions further triggers oxidative stress and cell damage [34], ultimately leading to ferroptosis. Research has shown that glycyrrhetinic acid protects against neuronal damage in TLE mice by inhibiting ferroptosis through the miR-194-5p/PTGS2 axis [35].

IL-6 is a protective cytokine, and some studies have found that in certain epilepsy patients, seizure activity is related to an inflammatory response within the body. Previous reports have shown an increase in IL-6 concentration in the cerebrospinal fluid of patients with tonic-clonic seizures [36], and there are complex interactions between the immune system and the nervous system. Inflammatory factors such as IL-6 may affect the development, survival, and function of neurons through various pathways. These interactions may have an impact on the onset and progression of epilepsy [37].

Toll-like Receptor 4 (TLR4) is a crucial immune receptor and a member of the Toll-like receptor family. It plays a vital role in the human immune system, primarily in recognizing and responding to foreign pathogens and other molecules associated with infection and inflammation. Existing research has shown a significant increase in serum levels of TLR4 in cases of drug-resistant epilepsy, indicating its potential as a novel epilepsy biomarker [38]. Additionally, in animal experiments, targeting TLR4-related pathways has shown promise in influencing the onset and progression of epilepsy [39, 40].

JUN is a transcription factor and a member of the Activator Protein-1 (AP-1) protein family. JUN proteins play a crucial regulatory role within cells and are involved in various biological processes, including cell proliferation, differentiation, apoptosis, inflammation, and stress responses. The JUN protein family comprises multiple members, with the most important ones being c-Jun, JunB, and JunD [41].Studies have shown that after kainic acid-induced seizures, rat hippocampal neurons strongly express immediate early gene products, such as c-FOS and c-Jun [42]. Additionally, inhibiting c-Jun N-terminal kinase (JNK) has demonstrated significant anti-epileptic effects [43].

Our qRT-PCR results indicate an upregulation of six genes during the acute phase of TLE. Since some of the genes we discovered can play a positive role in neuroinflammation or ferroptosis, we hypothesize that the upregulation of these genes represents a compensatory response by the body to counter the harm caused by TLE. For example, the upregulation of the Nrf2/HO-1 pathway can have antioxidative effects [44], and the upregulation of TLR4 during the acute phase of TLE is a key activator of the inflammatory pathway [45]. As for the upregulation of the PTGS2 gene, we speculate that it is related to ferroptosis, as PTGS2 is a crucial factor driving ferroptosis. Additionally, the activation of the JNK signaling pathway forms the foundation for ferroptosis [46]. Current research has also found that IL-6 can induce ferroptosis [47]. Therefore, we propose that these genes promote ferroptosis during the acute phase of TLE.

However, there are several limitations to this study. Certainly, your statement stands: the results obtained from our analysis still necessitate validation through additional experimental research to confirm our conclusions. Additionally, this study utilized brain tissue samples from animal models and peripheral blood samples from humans, lacking validation on human cerebrospinal fluid or brain tissue samples. Further research and optimization of the constructed ferroptosis-related models will require additional data from human brain tissue or cerebrospinal fluid samples. Thirdly, the absence of IL6 data in GSE143272 has resulted in uncertainty regarding the applicability of the ferroptosis-related model we constructed in human peripheral blood samples. Further investigation is needed to address this limitation and assess the model’s suitability in a broader context. Lastly, while we employed kainic-induced epileptic mouse samples to verify the expression of candidate genes, it is preferable to confirm these findings through clinical samples in future studies.

In conclusion, although this study has certain strengths, further research and validation are required to determine the accuracy and reliability of ferroptosis-related biomarkers in the diagnosis of TLE.


Based on the present results, ferroptosis-related hub genes NFE2L2, PTGS2, IL6, JUN, HMOX1, and TLR4 have good diagnostic value for epilepsy and may be potential early biological diagnostic targets. However, the experimental results of this study still require larger-scale basic experiments and clinical research in the future to further confirm.

Data Availability

Data of the present study are publicly available and may also be available from the corresponding authors upon request.


  1. Collaborators GBDE. Global, regional, and national burden of Epilepsy, 1990–2016: a systematic analysis for the global burden of Disease Study 2016. Lancet Neurol. 2019;18:357–75.

    Article  Google Scholar 

  2. Chong D, Jones NC, Schittenhelm RB, Anderson A, Casillas-Espinosa PM. Multi-omics integration and Epilepsy: towards a better understanding of biological mechanisms. Prog Neurobiol. 2023;227:102480.

    Article  CAS  PubMed  Google Scholar 

  3. Rugg-Gunn F, Miserocchi A, McEvoy A. Epilepsy Surgery. Pract Neurol. 2020;20:4–14.

    PubMed  Google Scholar 

  4. Martinez B, Peplow PV. MicroRNAs as potential biomarkers in temporal lobe Epilepsy and mesial temporal lobe Epilepsy. Neural Regen Res. 2023;18:716–26.

    Article  CAS  PubMed  Google Scholar 

  5. Ni S, Yuan Y, Song S, Li X. A double-edged sword with a therapeutic target: iron and ferroptosis in immune regulation. Nutr Rev. 2023;81:587–96.

    Article  PubMed  Google Scholar 

  6. Roggenhofer E, Toumpouli E, Seeck M, Wiest R, Lutti A, Kherif F, Novy J, Rossetti AO, Draganski B. Clinical phenotype modulates brain’s myelin and iron content in temporal lobe Epilepsy. Brain Struct Funct. 2022;227:901–11.

    Article  CAS  PubMed  Google Scholar 

  7. Chen S, Chen Y, Zhang Y, Kuang X, Liu Y, Guo M, Ma L, Zhang D, Li Q. Iron Metabolism and Ferroptosis in Epilepsy. Front Neurosci. 2020;14:601193.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Jin Y, Ren L, Jing X, Wang H. Targeting ferroptosis as novel therapeutic approaches for Epilepsy. Front Pharmacol. 2023;14:1185071.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Wang L, Duan C, Wang R, Chen L, Wang Y. Inflammation-related genes and immune infiltration landscape identified in kainite-induced temporal lobe Epilepsy based on integrated bioinformatics analysis. Front Neurosci. 2022;16:996368.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W. Smyth, limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Zhou N, Bao J. FerrDb: a manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database (Oxford) 2020 (2020).

  12. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–92.

    Article  CAS  PubMed  Google Scholar 

  13. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kalozoumi G, Kel-Margoulis O, Vafiadaki E, Greenberg D, Bernard H, Soreq H, Depaulis A, Sanoudou D. Glial responses during epileptogenesis in Mus musculus point to potential therapeutic targets. PLoS ONE. 2018;13:e0201742.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Cai Y, Yang Z. Ferroptosis and its role in Epilepsy. Front Cell Neurosci. 2021;15:696889.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hallis SP, Kim JM, Kwak MK. Emerging role of NRF2 Signaling in Cancer Stem Cell phenotype. Mol Cells. 2023;46:153–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ngo V, Duennwald ML. Nrf2 and oxidative stress: a General Overview of mechanisms and implications in Human Disease. Antioxid (Basel) 11 (2022).

  19. Thiruvengadam M, Venkidasamy B, Subramanian U, Samynathan R, Ali Shariati M, Rebezov M, Girish S, Thangavel S, Dhanapal AR, Fedoseeva N, Lee J, Chung IM. Bioactive compounds in oxidative stress-mediated Diseases: targeting the NRF2/ARE signaling pathway and epigenetic regulation. Antioxid (Basel) 10 (2021).

  20. Song X, Long D. Nrf2 and ferroptosis: a New Research Direction for neurodegenerative Diseases. Front Neurosci. 2020;14:267.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Uruno A, Yamamoto M. The KEAP1-NRF2 system and neurodegenerative Diseases. Antioxid Redox Signal. 2023;38:974–88.

    Article  CAS  PubMed  Google Scholar 

  22. Jiang YQ, Yang XY, Duan DQ, Zhang YY, Li NS, Tang LJ, Peng J, Luo XJ. Inhibition of MALT1 reduces ferroptosis in rat hearts following ischemia/reperfusion via enhancing the Nrf2/SLC7A11 pathway. Eur J Pharmacol. 2023;950:175774.

    Article  CAS  PubMed  Google Scholar 

  23. Zhao L, Tao X, Wan C, Dong D, Wang C, Xi Q, Liu Y, Song T. Astaxanthin alleviates inflammatory pain by regulating the p38 mitogen-activated protein kinase and nuclear factor-erythroid factor 2-related factor/heme oxygenase-1 pathways in mice. Food Funct. 2021;12:12381–94.

    Article  CAS  PubMed  Google Scholar 

  24. Liu J, Kang R, Tang D. Signaling pathways and defense mechanisms of ferroptosis. FEBS J. 2022;289:7038–50.

    Article  CAS  PubMed  Google Scholar 

  25. Anandhan A, Dodson M, Shakya A, Chen J, Liu P, Wei Y, Tan H, Wang Q, Jiang Z, Yang K, Garcia JG, Chambers SK, Chapman E, Ooi A, Yang-Hartwich Y, Stockwell BR, Zhang DD. NRF2 controls iron homeostasis and ferroptosis through HERC2 and VAMP8. Sci Adv. 2023;9:eade9585.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Xie R, Zhao W, Lowe S, Bentley R, Hu G, Mei H, Jiang X, Sun C, Wu Y, Yueying L. Quercetin alleviates kainic acid-induced seizure by inhibiting the Nrf2-mediated ferroptosis pathway. Free Radic Biol Med. 2022;191:212–26.

    Article  CAS  PubMed  Google Scholar 

  27. Yang M, Chen Y, Huang X, Shen F, Meng Y. ETS1 ameliorates Hyperoxia-Induced Bronchopulmonary Dysplasia in mice by activating Nrf2/HO-1 mediated ferroptosis. Lung. 2023;201:425–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. An Y, Luo Q, Han D, Guan L. Abietic acid inhibits acetaminophen-induced liver injury by alleviating inflammation and ferroptosis through regulating Nrf2/HO-1 axis. Int Immunopharmacol. 2023;118:110029.

    Article  CAS  PubMed  Google Scholar 

  29. Li D, Bai X, Jiang Y, Cheng Y. Butyrate alleviates PTZ-induced mitochondrial dysfunction, oxidative stress and neuron apoptosis in mice via Keap1/Nrf2/HO-1 pathway. Brain Res Bull. 2021;168:25–35.

    Article  CAS  PubMed  Google Scholar 

  30. Dai H, Wang P, Mao H, Mao X, Tan S, Chen Z. Dynorphin activation of kappa opioid receptor protects against Epilepsy and seizure-induced brain injury via PI3K/Akt/Nrf2/HO-1 pathway. Cell Cycle. 2019;18:226–37.

    Article  CAS  PubMed  Google Scholar 

  31. Ferrer MD, Busquets-Cortes C, Capo X, Tejada S, Tur JA, Pons A, Sureda A. Cyclooxygenase-2 inhibitors as a therapeutic target in Inflammatory Diseases. Curr Med Chem. 2019;26:3225–41.

    Article  CAS  PubMed  Google Scholar 

  32. Hirano T, Hirayama D, Wagatsuma K, Yamakawa T, Yokoyama Y, Nakase H. Immunological mechanisms in inflammation-Associated Colon carcinogenesis. Int J Mol Sci 21 (2020).

  33. Hashemi Goradel N, Najafi M, Salehi E, Farhood B, Mortezaee K. Cyclooxygenase-2 in cancer: a review. J Cell Physiol. 2019;234:5683–99.

    Article  CAS  PubMed  Google Scholar 

  34. Chen L, Zhao Q, Du X, Chen XI, Jiao Q, Jiang H. Effects of oxidative stress caused by iron overload on arachidonic acid metabolites in MES23.5 cells. J Biosci 47 (2022).

  35. Yi TT, Zhang LM, Huang XN. Glycyrrhizic acid protects against temporal lobe Epilepsy in young rats by regulating neuronal ferroptosis through the miR-194-5p/PTGS2 axis. Kaohsiung J Med Sci. 2023;39:154–65.

    Article  CAS  PubMed  Google Scholar 

  36. Peltola J, Palmio J, Korhonen L, Suhonen J, Miettinen A, Hurme M, Lindholm D, Keranen T. Interleukin-6 and interleukin-1 receptor antagonist in cerebrospinal fluid from patients with recent tonic-clonic seizures. Epilepsy Res. 2000;41:205–11.

    Article  CAS  PubMed  Google Scholar 

  37. Kummer KK, Zeidler M, Kalpachidou T, Kress M. Role of IL-6 in the regulation of neuronal development, survival and function. Cytokine. 2021;144:155582.

    Article  CAS  PubMed  Google Scholar 

  38. Kamasak T, Dilber B, Yaman SO, Durgut BD, Kurt T, Coban E, Arslan EA, Sahin S, Karahan SC, Cansu A. HMGB-1, TLR4, IL-1R1, TNF-alpha, and IL-1beta: novel Epilepsy markers? Epileptic Disord. 2020;22:183–93.

    Article  PubMed  Google Scholar 

  39. Yu L, Yang J, Yu W, Cao J, Li X. Rhein attenuates PTZ–induced Epilepsy and exerts neuroprotective activity via inhibition of the TLR4-NFkappaB signaling pathway. Neurosci Lett. 2021;758:136002.

    Article  CAS  PubMed  Google Scholar 

  40. Zhou Q, Wang Q, He B, Kong H, Luo H, Wang X, Wang W. MicroRNA 322-5p reduced neuronal inflammation via the TLR4/TRAF6/NF-kappaB axis in a rat Epilepsy model. Open Med (Wars). 2022;17:907–14.

    Article  CAS  PubMed  Google Scholar 

  41. Meng Q, Xia Y. c-Jun, at the crossroad of the signaling network. Protein Cell. 2011;2:889–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Lee MC, Rho JL, Kim MK, Woo YJ, Kim JH, Nam SC, Suh JJ, Chung WK, Moon JD, Kim HI. c-JUN expression and apoptotic cell death in kainate-induced temporal lobe Epilepsy. J Korean Med Sci. 2001;16:649–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Tai TY, Warner LN, Jones TD, Jung S, Concepcion FA, Skyrud DW, Fender J, Liu Y, Williams AD, Neumaier JF, D’Ambrosio R, Poolos NP. Antiepileptic action of c-Jun N-terminal kinase (JNK) inhibition in an animal model of temporal lobe Epilepsy. Neuroscience. 2017;349:35–47.

    Article  CAS  PubMed  Google Scholar 

  44. Zhang Q, Liu J, Duan H, Li R, Peng W, Wu C. Activation of Nrf2/HO-1 signaling: an important molecular mechanism of herbal medicine in the treatment of Atherosclerosis via the protection of vascular endothelial cells from oxidative stress. J Adv Res. 2021;34:43–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhu K, Zhu X, Sun S, Yang W, Liu S, Tang Z, Zhang R, Li J, Shen T, Hei M. Inhibition of TLR4 prevents hippocampal hypoxic-ischemic injury by regulating ferroptosis in neonatal rats. Exp Neurol. 2021;345:113828.

    Article  CAS  PubMed  Google Scholar 

  46. Liang Y, Liu Z, Qu L, Wang Y, Zhou Y, Liang L, Guo Y, Tang L. Inhibition of the IRE1/JNK pathway in renal tubular epithelial cells attenuates ferroptosis in acute kidney injury. Front Pharmacol. 2022;13:927641.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li M, Jin S, Zhang Z, Ma H, Yang X. Interleukin-6 facilitates Tumor progression by inducing ferroptosis resistance in head and neck squamous cell carcinoma. Cancer Lett. 2022;527:28–40.

    Article  CAS  PubMed  Google Scholar 

Download references


We extend our heartfelt gratitude to the GEO database for generously providing their platforms and to all the contributors who have uploaded their invaluable datasets.



Author information

Authors and Affiliations



S.C. designed the study. S.C.and X.J. analyzed the data and prepared the original manuscript. S.C. conducted animal experiments. T.H. and M.Z. revised and edited the paper. H.X. supervised this work. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Huiqin Xu.

Ethics declarations

Ethics approval and consent to participate

The experiments were conducted in accordance with the Regulations of the Chinese Council on Animal Care, approved by the Institutional Animal Care and Use Committee of Wenzhou Medical University, and reported following the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, S., Jin, X., He, T. et al. Identification of ferroptosis-related genes in acute phase of temporal lobe epilepsy based on bioinformatic analysis. BMC Genomics 24, 675 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: