- 1Department of Neurology, The First Affiliated Hospital of Kunming Medical University, Kunming, Yunnan, China
- 2Yunnan Provincial Clinical Research Center for Neurological Disease, Kunming, Yunnan, China
- 3Department of Emergency, The First Affiliated Hospital of Kunming Medical University, Kunming, Yunnan, China
One of the most prevalent types of epilepsy is temporal lobe epilepsy (TLE), which has unknown etiological factors and drug resistance. The detailed mechanisms underlying potassium channels in human TLE have not yet been elucidated. Hence, this study aimed to mine potassium channel genes linked to TLE using a bioinformatic approach. The results found that Four key TLE-related potassium channel genes (TERKPCGs) were identified: potassium voltage-gated channel subfamily E member (KCNA) 1, KCNA2, potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11), and KCNS1. A protein–protein interaction (PPI) network was constructed to analyze the relationship between TERKPCGs and other key module genes. The results of gene set enrichment analysis (GSEA) for a single gene indicated that the four TERKPCGs were highly linked to the cation channel, potassium channel, respiratory chain, and oxidative phosphorylation. The mRNA-TF network was established using four mRNAs and 113 predicted transcription factors. A ceRNA network containing seven miRNAs, two mRNAs, and 244 lncRNAs was constructed based on the TERKPCGs. Three common small-molecule drugs (enflurane, promethazine, and miconazole) target KCNA1, KCNA2, and KCNS1. Ten small-molecule drugs (glimepiride, diazoxide, levosimendan, and thiamylal et al.) were retrieved for KCNJ11. Compared to normal mice, the expression of KCNA1, KCNA2, KCNJ11, and KCNS1 was downregulated in the brain tissue of the epilepsy mouse model at both the transcriptional and translational levels, which was consistent with the trend of human data from the public database. The results indicated that key potassium channel genes linked to TLE were identified based on bioinformatics analysis to investigate the potential significance of potassium channel genes in the development and treatment of TLE.
Introduction
The second most prevalent neurological disorder is epilepsy, and temporal lobe epilepsy (TLE) is one of the most frequent types in adults (1). TLE comorbidity increases the occurrence of cognitive impairment and behavioral alterations, resulting in low quality of life, career, and finances (2). Although antiepileptic drug treatment, vagus nerve stimulation, and resection of the epileptogenic focus have been applied in TLE, some epilepsy cases still cannot be effectively treated with frequent recurrent seizures (3). Therefore, new strategies are required for early diagnosis and treatment of TLE.
Potassium channels in cell membranes play a critical role in regulating neuronal excitability by controlling the transport of potassium ions across cells (4–6). Several potassium channels are associated with TLE, primarily based on animal experiments (7–10). Research on the association between potassium levels and epilepsy may provide novel insights into the mechanisms and treatments of epilepsy.
However, the mechanisms underlying potassium channels in human TLE have not been elucidated. In this study, we hypothesized that key potassium channel genes are associated with epilepsy pathogenesis. The purpose of this study was to analyze the expression of potassium channels in individuals with TLE from the gene expression omnibus (GEO), and to determine the key potassium channel genes to uncover a new therapeutic target for TLE.
Methodology
Data source
The GSE134697 dataset (11), containing two normal neocortex samples from patients without epilepsy and 17 neocortex samples from patients with TLE, was obtained from GEO.1 The Molecular Signatures Database (MsigDB)2 was used to obtain the file “c5.go.mf.v7.4. symbols.gmt,” and potassium channel-related genes were extracted by referring to the following pathway: GOMF_POTASSIUM_CHANNEL_ACTIVITY, GOMF_POTASSIUM_CHANNEL_INHIBITOR_ACTIVITY and GOMF_POTASSIUM_CHANNEL_REGULATOR_ACTIVITY. Single-cell RNA sequencing (scRNA-seq) data from the GSE140393 dataset were downloaded from the GEO database.
Weighted gene co-expression network analysis
The transcriptome data of two normal neocortex samples and 17 TLE neocortex samples were analyzed using the R platform’s “WGCNA” (version 1.69) (12). The “Good Samples Genes” function was utilized to carry out clustering of the sample to detect outliers and remove them. To make the co-expression network contented with the distribution of a scale-free network, soft-thresholding power was computed using the pick Soft Threshold function and verified by the correlation between k and p(k) (13). Various modules were identified using the dynamic tree-cutting method, with at least 30 genes in each module (14). Subsequently, to merge similar modules, the threshold was set as 0.2 (15, 16). The associations between these modules the two clinical characteristics were further analyzed. The module with the highest Pearson correlation coefficient was chosen for the subsequent analysis.
Differentially expressed genes (DEGs) identification
The “limma” in the R package (version 3.44.3) with the voom and empirical Bayes (eBayes) methods was utilized to determine the differentially expressed genes (DEGs) in normal neocortex samples and TLE neocortex samples in the GSE134697 dataset (17). The log2-fold changes (log2FC), i.e., |log2FC| ≥ 0.5 and p-value ≤0.05 were deemed as cut-off criteria. A heatmap and a volcano plot were used to display the results.
Functional enrichment analysis
GO is an important bioinformatics tool for annotating genes and analyzing their products and sequences to understand fundamental biological events such as biological processes (BP), molecular functions (MF), and cellular components (CC). To perform GO analysis, “clusterProfiler” in the R package (version 3. 16.0) was used (18). p-value <0.05 was The threshold for enrichment significance was set at p < 0.05.
To explore the potential role of TLE related key potassium channel genes (TERKPCGs), Gene Set Enrichment Analysis (GSEA) was performed for a single gene in GSEA (version 2.09). The expression values for each gene were considered as phenotype files, and the correlation coefficients of every gene with other genes were ordered in the gene sets. The “c5.all.v6. 1.symbols.gmt” from MsigDB were selected as the reference gene set. |NES| > 1, NOM value of p <0.05, and FDR q-val < 0.25 were set as the thresholds for enrichment significance.
Construction and analysis of protein–protein interaction network
The Search Tool for the Retrieval of Interacting Genes (STRING) database3 was used to build a protein–protein interaction (PPI) network after hiding independent nodes (19). For the construction of the network, nodes with an interaction confidence >0.4 were used. The data were imported into the Cytoscape software (version 3.9. 1) for visual presentation.
Construction of TF-mRNA and ceRNA regulatory network
To further explore the upstream regulatory molecules of TERKPCGs expression, the upstream transcription factors (TFs) and non-coding RNAs were mined. TFs of the TERKPCGs were predicted using the online website KnockTF.4 Target miRNAs of TERKPCGs were predicted using TargetScan and miRDB in the miRWalk database5 (20), and it was determined that the miRNAs found in the two databases were the target miRNAs. The StarBase database was used to predict lncRNAs upstream of the miRNAs. Subsequently, lncRNA-miRNA-mRNA and TF-mRNA networks were constructed using the Cytoscape software.
Identification of phosphorylated sites
Protein modifications and their regulation are associated with protein function. Multiple candidate phosphorylation sites in TERKPCGs have been identified in a post-translational modification database called PhosphoSitePlus6 (21).
Small molecule drugs prediction
To further explore potential drugs that target TERKPCGs and treat TLE, the DrugBank platform,7 containing comprehensive molecular data on drugs, their mechanisms, their interactions, and their targets, was used. Small-molecule drugs targeting TERKPCGs were predicted using the DrugBank database (20).
Animals and in vivo study
Male C57 mice were fed and given water as per need and reared at a temperature of 23 C to 25 C, humidity of 45 to 55%, and a 12-h light/dark cycle in a quiet room at Kunming Medical University.
The mice were administered lithium (127 mg/kg) and atropine 18–20 h intraperitoneally (i.p.) before administration of pilocarpine (50 mg/kg, i.p.). Ninety minutes after status epilepticus (SE) establishment, mice were administered diazepam (5 mg/kg, i.p.) to relieve seizures. Behavioral seizures were evaluated according to the Racine limbic seizure categorization, and SE criteria were evaluated when the animal started exhibiting sustained behaviors or higher limbic seizures of score IV or V, with no spontaneous recovery. Similar to the experimental mice, the control mice were injected with NaCl (0.9%) rather than pilocarpine. Following acute exacerbation, the mice were administered 5% glucose and sodium chloride (5 mL, once daily, i.p.) during the incubation period. Upon entering the chronic phase (after 28 days), the mice were sacrificed by cervical dislocation and anesthetized with 1% pentobarbital sodium (50 mg/kg body mass; Wuhan Dinghui Chemical Co., Ltd., Wuhan, China). Brain tissues were collected from six mice in each group to determine the expression levels of TERKPCGs. This protocol was approved by the Committee of Ethics of Animal Experiments at Kunming Medical University (Permit Number: SCXK 2011-0015). Eight mice showed spontaneous seizures in the chronic phase and their behavior was monitored using a camera. Brain samples were collected from mice according to Racine limbic seizure categorization. Two mice died and were thus excluded from the treatment group.
Extraction of RNA and quantitative real-time polymerase chain reaction
Total RNA was extracted from six normal mouse brain tissues and six epileptic mouse brain tissues using Nuclezol LS RNA Isolation Reagent according to the manufacturer’s guidelines (ABP Biosciences Inc.). As per the given protocol, total RNA was converted into cDNA through reverse transcription employing the SureScript-First-strand-cDNA-synthesis-kit (GeneCopoeia). qPCR was then implemented using BlazeTaq™ SYBR® Green qPCR Mix 2.0 (GeneCopoeia). The thermocycling program employed for qPCR was as follows: one cycle of initial denaturation for 30 s at 95°C, followed by 40 cycles of denaturation for 10 s at 95°C, annealing for 20 s at 60°C, and extension for 30 s at 72°C. The primer sequences used are listed in Table 1. The corresponding expression levels were normalized to the endogenous control GAPDH and calculated using the 2−ΔΔCq method (22). The differences between the two groups were compared using the Student’s t-test. In the two-tailed statistical analysis, a p-value <0.05 was deemed statistically significant.
Immunohistochemistry
Hippocampal specimens from mice were polyformaldehyde-fixed, immersed in paraffin, and sliced into 4 μm sections. The sections were dewaxed, rehydrated with xylene, and washed with alcohol and PBS gradient. After heating in citric acid buffer (pH 6.0) for 15 min, antigens were retrieved. TBS/H2O2 was used to block endogenous peroxidase activity. After incubation with primary antibodies for 1 h at room temperature, the sections were incubated with secondary antibodies for 30 min at 37°C. Primary antibodies against mouse monoclonal KCNA1 (bs-8691R), KCNA2 (bs- 12182R), KCNJ11 (bs-2436R), and KCNS1 (bs- 16916R) (1:1000) were purchased from BIOSS (Beijing, China). Rabbit anti-mouse IgG H&L (HRP) (bs-0296R) incubation kit was purchased from ZSGB-Bio (Beijing, China). The antibodies were used according to the manufacturer‟s guidelines. Immunohistochemical staining was performed by exposing the sections to diaminobenzidine (DAB) as a chromogen. The sections were counterstained with hematoxylin, dehydrated in ethanol, and clarified in xylene. Observations and photographs were obtained using a microscope (Zeiss, Lab. A1). According to a previously described method (23), after staining, five high-power fields (×200) were randomly selected on each slide, and the average proportion of positive expression in each field was counted using the true color multi-functional cell image analysis management system (Image-Pro Plus; Media Cybernetics, Rockville, MD, USA).
Statistics analysis
Differences in TF expression between groups were analyzed using the Wilcoxon test. Student’s t-test was conducted for differences in the expression of key TERKPCGs in clinical mouse brain tissues and hippocampal specimens. Statistical analysis was employed with p < 0.05 as a significant difference.
Results
Identification of the key module genes for TLE
To identify the key modules highly linked to TLE, we structured WGCNA on the genes in the GSE134697 dataset. No obvious outlier samples were removed through sample clustering, and β = 6 (scale-free R2 = 0.9) was selected to develop a scale-free network (Figures 1A,B). Based on the similar expression patterns of all genes in the cohorts, average dynamic tree cutting and hierarchical clustering were conducted, in which the retrieved modules were merged to a minimum size of 30. Following eigengene calculation, 38 modules were clustered and merged (Figure 1C). Subsequently, the association between each module and two clinical characteristics (TLE and healthy controls) was assessed, suggesting that MEindianred4 (|cor| = 0.77, p = 0.0001) showed the strongest association with TLE, and the close correlation of gene significance and module significance was displayed with cor = 0.56 (p < 0.0001) (Figures 1D,E). Thus, 151 genes from this module were used for further analyses (Supplementary Table S1).
Figure 1. Gene co-expression network analysis (WGCNA) and key modules identification. (A) Sample clustering dendrogram and trait heatmap of cohorts in GSE134697 to remove outliers. (B) Analysis of the scale-free fit index (left) and the mean connectivity (right) for various soft-thresholding powers. (C) Clustering dendrogram of all genes based on similar expression patterns, with assigned module colors. (D) The correlation between different modules and clinical traits. (E) Scatterplot of Gene Significance (GS) and Module Membership (MM) in MEindianred4.
Identifying and exploration of DEGs in TLE
We further performed differential expression analysis of gene expression data from the GSE134697 dataset. Based on the established criteria, 459 DEGs (containing 144 upregulated and 315 downregulated genes) were identified between the normal and TLE samples (Figures 2A,B; Supplementary Table S2). Subsequently, 67 common genes from 151 key module genes and these DEGs were selected, namely, TLE-related DEGs (Figure 2C; Supplementary Table S3). As shown in Supplementary Table S4, in the CC ontology, the TLE-related DEGs were significantly associated with the ion channel complex, transmembrane transporter complex, and potassium channel complex. Concerning MF, these genes were related to metal ion transmembrane transporter activity, ion channel activity, and voltage-gated potassium channel activity.
Figure 2. Identification of differentially expressed genes (DEGs) in TLE. (A) The volcano map of DEGs between normal and TLE samples. (B) The heatmap of top100 DEGs. (C) Venn diagram for TLE-related DEGs.
Identification of TERKPCGs and gene set enrichment analysis
Next, we downloaded 148 potassium channel-related genes from MsigDB (Supplementary Table S5) and selected the intersection with 67 TLE-related DEGs to obtain the TLE-related key potassium channel genes, where KCNA1, KCNA2, KCNJ11, and KCNS1 were selected as TERKPCGs (Figure 3A). A PPI network was constructed to analyze the interrelationship between TERKPCGs and other key module genes (Supplementary Figure S1).
Figure 3. Identification and analysis of TLE-related key potassium channel genes (TERKPCGs). (A) Venn diagram for TERKPCGs. (B) Single-gene Gene set enrichment analysis (GSEA) of four TERKPCGs. (Top) Enrichment score (ES) line chart: the display of the ES value as it is calculated to each position when the analysis is sorted along the ranking list. The score at the highest peak (the farthest vertical distance) is the ES value of the gene set. (Middle) Lines are used to mark the position of members from the gene set in the gene sorting list, and the black line represents that the genes from the sorting gene table exist in the functional annotation gene set of the current analysis. Leading edge subset is the corresponding gene from (0,0) to the highest peak ES. (Bottom) After sorting, the rank values of all genes were distributed. The genes corresponding to the red part of the heat map were highly positively correlated to the TERKPCG, and the genes corresponding to the blue part were highly negatively correlated to the TERKPCG. The pearson correlation coefficient of each gene was displayed in a gray area map.
Furthermore, GSEA for a single gene in the expression data of the GSE134697 dataset based on Gene Ontology gene sets was conducted to evaluate the potential mechanism of TERKPCGs in TLE. The detailed enrichment results of GO terms for each TERKPCG are shown in Supplementary Table S6, and the top 10 were positively correlated with each TERKPCG, as displayed in Table 2. We noted that the functional term mostly correlating with KCNA1 was the cation channel complex, and voltage-gated cation channel activity was mostly linked to KCNA2. Meanwhile, the respiratory chain was most related to KCNJ11, and the top functional term relevant to KCNS1 was the potassium channel complex (Figure 3B).
Table 2. The top10 enrichment results of each TLE-related key potassium channel gene (TERKPCG) by Gene set enrichment analysis (GSEA).
The underlying regulatory mechanism targeting four TERKPCGs
To further investigate the underlying regulatory network targeting key TERKPCGs, the potential TFs binding sites and correlated miRNA targets were predicted using the KnockTF and miRWalk tools. As shown in Supplementary Figure S2, a TF mRNA network was established using four TERKPCGs and 113 predicted transcription factors. Interestingly, these TERKPCGs were commonly regulated by four transcription factors (TFAP4, FOXM1, POSTN, and PTEN), although their expression levels were not significantly different between the case and control groups (Supplementary Figure S3). On the other hand, six miRNAs (hsa-miR-204-5p, hsa-miR-128-3p, hsa-miR-373-3p, hsa-miR-92b-3p, hsa-miR-543, and hsa-miR-4,319) targeting KCNA1 and one miRNA, hsa-let- 7e-5p targeting KCNJ11 were predicted, and corresponding lncRNAs upstream of which were shown in Supplementary Figure S4, including the regulatory axes LINC02242-hsa-let-7e-5p-KCNJ11 and LINC01963-hsa-miR-204-5p-KCNA1. Meanwhile, the TOP4 lncRNAs common to most miRNAs are exhibited in Figure 4, for example, NEAT1, XIST, AC00728.2, and NORAD, providing new information for studying the role of TF regulation and ceRNA regulatory axes in TLE development.
Moreover, considering the influential role of protein post-translational modifications (PTMs) and potassium channels in the nervous system (24, 25), an in-depth study of PTMs is essential for elucidating the functional interactions of the key TERKPCGs in the pathogenesis of TLE. The PhosphoSitePlus online website was employed to identify the phosphorylation and ubiquitination modifications of the four TERKPCGs. Detailed modification loci for each protein are presented in Table 3. Phosphorylation modifications of serine, threonine, tyrosine, and lysine were mainly relevant to the four TERKPCGs (Figure 4).
Drugs prediction
To discover potential drugs that can function in TLE, we identified potential drugs targeting the proteins encoded by TERKPCGs, as shown in Supplementary Table S7. Considering the downregulation patterns of four TERKPCGs, seven drugs with activator actions were selected and imaged (Figure 4), where promethazine might reverse the expression of KCNA2, KCNS1, and diazoxide, levosimendan, and ritodrine were considered as inducers of KCNJ11. Among isoflurane, methoxyflurane and desflurane target KCNA1. However, further animal studies and clinical validation are required for clinical application.
The expression of TERKPCGs in the mouse model with TLE
Eventually, the expression patterns of four TERKPCGs identified above was exhibited in Figure 5A, where among the expressions of KCNA1, KCNA2, KCNJ11, and KCNS1 were decreased in TLE cohorts of GSE134697. Moreover, the expression of the four TERKPCGs at the transcriptional and translational levels was further confirmed using quantitative real-time polymerase chain reaction (qRT-PCR) and IHC analyses of clinical samples. We first established a mouse model of epilepsy and simultaneously established a control group. Brain tissues from six control mice and six epilepsy mouse models were used for the experiments. As shown in Figure 5B, KCNA1, KCNA2, KCNJ11, and KCNS1 were all downregulated at the mRNA level in the brain tissue of epileptic mice compared to control mice. IHC results indicated downregulation of the expression of the four genes at the protein level in the brain tissue of epileptic mice. This is also illustrated by the statistics on the positive rate of IHC in the group of six mice (Figure 5C). The expression profiles of KCNA1, KCNA2, KCNJ11, and KCNS1 in mice were consistent with the results of the human data from the GSE134697 dataset.
Figure 5. Expression verification of four TERKPCGs. (A) Boxplot for the expressions of four TERKPCGs in GSE134697. (B) The relative mRNA expression level of four TERKPCGs in brain tissue of control and epilepsy model mice. (C) Protein expression levels of TERKPCGs in brain tissue of control and epilepsy model mice, including the IHC images and the corresponding positive expression of IHC images. *p < 0.05, **p < 0.01, ***p < 0.001.
Discussion
Potassium channels exiting nerve cell membranes play a critical role in regulating neuronal excitability by controlling the transport of potassium ions across cells. Changes in potassium channel gene expression play an important role in the occurrence, development, treatment, and prevention of epilepsy. To date, there is a lack of studies on potassium channel genes in epilepsy. Therefore, this study aimed to probe potassium channel genes linked to TLE using a bioinformatic approach. The results indicated that key potassium channel genes linked to TLE were identified based on bioinformatics analysis to investigate the potential significance of potassium channel genes in the development and treatment of TLE.
In this study, 151 key module genes identified using WGCNA were significantly associated with TLE. A total of 459 differentially expressed genes (DEGs) were screened in TLE samples and compared with the controls. In total, 67 genes were obtained from the intersections of 151 key genes, and DEGs were considered crucial genes linked to TLE. GO functional annotations and KEGG pathway analyses of the crucial 67 genes linked to TLE were performed to uncover the mechanisms of TLE. The findings indicated that These genes mainly participated in many GO terms, including ion channel complexes, transporter complexes, transmembrane transporter complexes, and main axons. However, these crucial genes were not involved in KEGG pathway-related regulation.
Previous research has shown that three of the four key potassium channel genes (KCNA1, KCNA2, KCNJ11, and KCNS1) play vital roles in seizures and epilepsy (8–10). Four key potassium channel genes were screened based on the intersection of crucial 67 genes and genes associated with potassium channels from the MsigDB database. Based on previous studies (26–29), RT-PCR and immunohistochemistry analyses of the expression levels of four key genes in the brains of TLE mice validated the results obtained from bioinformatic analyses, showing that mRNA and protein expression of these four key genes were remarkably decreased in TLE mice compared to the control. The low expression levels of these four genes may be associated with the progression of epilepsy. An imbalance between inhibitory and excitatory neurons plays an important role in epilepsy (30, 31). Potassium channel downregulation might impair excitatory neuron repolarization and decrease the threshold for action potential (AP) firing, resulting in network hyperexcitability and increased seizure susceptibility (32).
KCNA1 encodes the alpha subunit of KV1. 1, which belongs to the Shaker subfamily (Kv1). KV1. 1 mediating the outflow of K+ regulates action potential propagation and shape, neuronal repetitive firing properties, membrane repolarization, and neurotransmitter release, which are positioned in the axons and axon terminals of the hippocampal trisynaptic circuit (33–35). Mutations in KCNA1 are associated with multiple epilepsies (36, 37). Studies have shown that KCNA1 targeted deletions result in TLE and sudden death in animal models of epilepsy (35–45). The KCNA1 gene encodes the alpha subunit of the Kv1. 1 channel (a delayed rectifier, voltage-gated potassium channel) in humans and mice, and its functional mutations are linked with partial epilepsy, episodic ataxia, sleep disorders, and sudden death in epilepsy (26, 38, 41–44, 46–49). In Kv1. 1 knockout (KO) mice, infrared telemetric actimetry was performed immediately after weaning and continued until the animals suddenly died. KCNA1-null mice with a targeted deletion of the Kv1. 1 potassium channel α subunit protein have severe spontaneous recurrent seizures (SRS) and are frequently used as epilepsy mouse models (33–39). Overexpression of Kv1. 1 can reduce the frequency of brief (1 s), and high-frequency epileptiform discharges (49). Therefore, upregulation of KCNA1 can successfully suppress seizures in a rodent model of intractable TLE (50, 51) and offers a potential new treatment option for TLE in humans. The finding in the current study that KCNA1 was downregulated in clinical case samples is consistent with previous literature, providing a theoretical basis for subsequent pathway analysis.
KCNA2, also a member of the Shaker subfamily (Kv1), encodes Kv1.2. Kv1.2 channels are widely expressed voltage-dependent potassium channels all over the central and peripheral nervous systems located within Kv1. 1 along the axons and axon terminals, as well as at presynaptic sites (52, 53). High-density clustering of Kv1. 1 and Kv1.2 is located in the juxta paranodal area neighboring the nodes of Ranvier of mammalian axons, at axonal initial segments, in the pinceau area of cerebellar basket cells, in the hippocampus, cortex, and auditory brainstem (44, 54). Kv1.2 are low-voltage activated gradually inactivating channels that open in response to small depolarizations near the resting potential (42, 55). They are involved in the initiation and shaping of action potentials, influencing the firing patterns of the action potential, controlling the excitability of neurons, and neurotransmitter release (56). Downregulation of Kv1.2 which likely increases seizure susceptibility, impairs excitatory neuron repolarization and decreases the threshold for AP firing, activating network hyperexcitability (8). Echoing the dysregulation results for Kv1.2, the mRNA and protein levels of KCNA2 were downregulated in the clinical TLE samples. In addition, Imbrici found that the E236K variant in KCNA2 exhibited voltage-dependent activation, shifting toward negative potentials and slower kinetics of deactivation and activation, causing early onset of epilepsy (32). Uysal reported that patients with mutations in KCNA2 showed early onset developmental and epileptic encephalopathies (57). KCNJ11 plays a key role in controlling K+ channel selectivity (57). More than 130 mutations have been reported in KCNJ11, which cause psychomotor developmental delay and epilepsy during the infantile period with the most severe form (58). However, the involvement of KCNJ11 and KCNS1 in epilepsy has not been reported thus far, and it was first found that the mRNA and protein levels of these two genes were downregulated in TLE samples compared to normal controls, which remains to be confirmed by further studies.
Gene set enrichment analysis (GSEA) revealed that the functions of KCNA1 and KCNA2 were enriched in signal transduction, including cation channel activity, potassium ion channel activity, and ion voltage gating. KCNJ11 functions in the respiratory chain, oxidative phosphorylation, mitochondrial intima protein complex, potassium channel complex, cation channel, and charge transfer. KCNS1 is primarily involved in the potassium channel complex and respiratory chain.
In this study, the phosphorylation modification of serine, threonine, tyrosine, and lysine was mainly relevant to the four TERKPCGs. The study found that Targeted post translational modifications of drugs have the potential to restore tetratricopeptide repeat-containing Rab8b-interacting protein phosphorylation to reduce TLE excitability (59). D-serine mitigates cell loss associated with temporal lobe epilepsy (60). Serine/threonine protein kinases are targets of rapamycin (mTOR), and mTOR activation induces the expression of voltage-gated potassium channel Kv1. 1 in TLE animal models, indirectly reflecting the correlation between TERKPCGs and serine/threonine phosphorylation (61). Transfusion of genetically modified tyrosine kinase receptor B into the hippocampus increases the susceptibility or severity of seizures (62). Controlling the expression of the silk/threonine kinase With-no-lysine kinase 3 may be a new therapeutic target for epilepsy (63). The above research on amino acid kinases provides a theoretical basis for the correlation between key TERKPCGs and amino acid phosphorylation in our study, and is expected to target key genes to further explore their case mechanisms in TLE.
To further investigate the mechanism of the transcriptional regulation of these four genes, we constructed TF mRNA and ceRNA networks. In the TF mRNA network, we found four identical transcription factors, including TFAP4, FOXM1, POSTN, and PTEN, which may regulate the transcription of the four genes. In Dravet syndrome-iPSC GABA lines, there was a remarkable upregulation of FOXM, which is a positive regulator of disrupted pathways for histone modification and cell cycle regulation (64). Glutamatergic density and strength of glutamatergic cells by PTEN. Loss of PTEN can significantly enhance the formation of excitatory synapses in neurons, making them hyperexcitable (65). POSTN promotes neural stem cell proliferation and differentiation following hypoxic–ischemic injury (66). Clarithromycin ameliorates early brain injury after subarachnoid hemorrhage via suppressing postn-related pathways in mice (67). TFAP4 focuses on various cancers, but no research has been conducted on brain injury or epilepsy (68, 69). In this study, the expression of four key TFs between the disease and control groups was not significantly different, indicating that changes in the expression of these transcription factors may not be the key factors leading to changes in the expression and mechanism of key TERKPCGs in TLE patients. Considering the small sample size of the dataset, this discovery needs to be further validated by collecting more clinical samples in the future. In the ceRNA network, possible ceRNA regulatory mechanisms of the two potassium channel genes were revealed, providing a basis for future research on specific mechanisms. Previous studies have found that microRNA-204 inhibits epileptic discharges by regulating TrkB-ERK1/2-CRB signals in cultured hippocampal neurons (66). MiR-204-5p is involved in the regulation of hippocampal sclerosis and focal cortical dysplasia, and participates in the occurrence and development of epilepsy (70). miRNA let-7e regulates the expression of caspase-3 during apoptosis in PC12 cells following anoxia/reoxygenation injury (71). TLE alters the expression levels of miRNA let-7e in the hippocampus, and deregulated miRNA let-7e may be involved in the pathogenesis of epilepsy directly or indirectly (72). Therefore, LINC02242-hsa-let-7e-5p-KCNJ11 and LINC01963-hsa-miR-204-5p-KCNA1 might be involved in the regulation of epilepsy and brain injury. However, the regulatory relationships between these factors require further exploration. In this study, the TOP4 lncRNAs common to most miRNAs were exhibited, the studies have found that The NEAT1/NORAD/XIST has-miR-204 axis regulates the development of Kawasaki disease (KD) (73). Further research has found that these three key lncRNAs together act as potential genes in interstitial lung disease (SSC ILD) and breast cancer in systemic sclerosis (74, 75). Therefore, we aimed to target the interactions between these three key genes to explore their biological significance in TLE.
The DrugBank database was used to identify the target drugs of the four key potassium channel genes to explore potential anti-epilepsy drugs. Molecular docking results showed that enflurane, promethazine, and miconazole targeted KCNA1, KCNA2, and KCNS1. Studies have found that enflurane has anticonvulsant effects in cat epilepsy models (76). There is little exacerbation of preexisting epileptic foci by enflurane, the only exception being the case of certain myoclonic-type epilepsies, such as progressive myoclonic epilepsy and photosensitive epilepsy (76). Promethazine reduces brain injury after ischemic stroke through PKC-δ/NOX/MnSOD and RIP1-RIP3 regulated activation of the NLRP3 inflammasome following ischemic stroke (77, 78). Miconazole exerts disease-modifying effects in epilepsy by suppressing neuroinflammation via the NF-κB pathway and iNOS production (79). The application of isoflurane as an anesthetic in TLE research has been reported for a long time (80). Volatile anesthetics such as isoflurane and desflurane have inhibitory effects on human Kv channels at various clinical concentrations (81). Our study suggests that their inhibitory effects on TLE may target KCNA1. Other drugs have not been further studied in TLE, and further confirmation is required.
However, this study has some limitations. First, the sample size used in this study was small, although supplemented with mouse models for testing, a larger sample of TLE patients was collected for further analysis. Second, the expression of key genes in neuronal cells has not been detected in the MsigDB database, which is inconsistent with in vitro experiments and needs to be confirmed by further research in the future. Thirdly, the binding sites among ceRNAs targeting key TERKPCGs were predicted and need to be examined using more functional experiments. Moreover, the clinical use of inducers of these TERKPCGs remains to be tested in TLE progression in the future.
Conclusion
We identified four key potassium channel genes, KCNA1, KCNA2, KCNJ11, and KCNS1, that are associated with the pathogenesis of epilepsy. Because our aim was to screen key genes related to potassium channel disorders, the four genes obtained may indirectly affect the balance between inhibitory and excitatory neurons, mainly through the relevant mechanism of potassium channel downregulation. Therapeutic strategies targeting these genes will be our future research direction.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was reviewed and approved by Committee of the Ethics of Animal Experiments of Kunming Medical University (Permit Number: SCXK 2011–0015). This study conformed to the ethical guidelines of the Science Foundation of National Natural Science of China. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
L-miZ, LC, L-meZ, and M-wL conceptualized and planned the study. LC, Y-fZ, and W-mD performed the experiments. Y-fZ and W-mD analyzed the data. L-miZ wrote the manuscript. All authors have read and approved the final version of the manuscript.
Funding
This work was supported by the Major Science and Technology Special Project of Yunnan Province under Grant (No. 202102AA100061), Nature Science Foundation of China under Grant (No. 82060252) and (No. 81960350), Yunnan Basic Research Projects under Grant (No. 2018FB115), and the Yunnan Health Training Project of High-level Talents under Grant (No. H-2018058), Yunnan Applied Basic Research Project-Union Foundation of China under Grant (No. 202201AY070001-091), and Applied Basic Research of Yunnan Neurological Disease Diagnosis and Treatment Center under Grant (No. ZX2019-03-05).
Acknowledgments
We would like to express our gratitude to Mei-xian Su, Professor of Critical Care Medicine University of the Second Affiliated Hospital of Kunming Medical University, Kunming China, for sharing their insights and expertise that greatly assisted the research and for revising and providing comments that markedly improved the manuscript.
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/fneur.2023.1175007/full#supplementary-material
Abbreviations
TLE, Temporal lobe epilepsy; GEO, Gene expression omnibus; GSEA, Gene Set Enrichment Analysis; SE, Status epilepticus; DEGs, Differentially expressed genes; PTMs, Protein post-translational modifications; qRT-PCR, Quantitative real-time polymerase chain reaction; MsigDB, Molecular signature database; STRING, Search tool for the retrieval of interacting genes; IHC, Immunohistochemistry; MF, Molecular function; CC, Cellular component; BP, Biological process.
Footnotes
1. ^https://www.ncbi.nlm.nih.gov/gds
2. ^http://www.gsea-msigdb.org/gsea/msigdb/index.jsp
4. ^http://www.licpathway.net/KnockTF/search.php
5. ^http://mirwalk.umm.uni-heidelberg.de/search_genes/
References
1. Ngugi, AK, Bottomley, C, Kleinschmidt, I, Sander, JW, and Newton, CR. Estimation of the burden of active and life-time epilepsy: a meta-analytic approach. Epilepsia. (2010) 51:883–90. doi: 10.1111/j.1528-1167.2009.02481.x
2. Vrinda, M, Arun, S, Srikumar, BN, Kutty, BM, and Shankaranarayana Rao, BS. Shankaranarayana Rao, temporal lobe epilepsy-induced neurodegeneration and cognitive deficits: implications for aging. J Chem Neuroanat. (2019) 95:146–53. doi: 10.1016/j.jchemneu.2018.02.005
3. Tramoni-Negre, E, Lambert, I, Bartolomei, F, and Felician, O. Long-term memory deficits in temporal lobe epilepsy. Rev Neurol (Paris). (2017) 173:490–7. doi: 10.1016/j.neurol.2017.06.011
4. Kuang, Q, Purhonen, P, and Hebert, H. Structure of potassium channels. Cell Mol Life Sci. (2015) 72:3677–93. doi: 10.1007/s00018-015-1948-5
5. Villa, C, and Combi, R. Potassium channels and human epileptic phenotypes: an updated overview. Front Cell Neurosci. (2016) 10:81. doi: 10.3389/fncel.2016.00081
6. Tiwari, D, Brager, DH, Rymer, JK, Bunk, AT, White, AR, Elsayed, NA, et al. MicroRNA inhibition upregulates hippocampal A-type potassium current and reduces seizure frequency in a mouse model of epilepsy. Neurobiol Dis. (2019) 130:104508. doi: 10.1016/j.nbd.2019.104508
7. Simeone, KA, Matthews, SA, Rho, JM, and Simeone, TA. Rho, loss of the Kv1. 1 potassium channel promotes pathologic sharp waves and high frequency oscillations in in vitro hippocampal slices. Epilepsia. (2016) 57:e178–82. doi: 10.1111/epi.13444
8. Liu, Y, Duan, Y, Du, D, and Chen, F. Rescuing Kv10.2 protein changes cognitive and emotional function in kainic acid-induced status epilepticus rats. Epilepsy Behav. (2020) 106:106894. doi: 10.1016/j.yebeh.2019.106894
9. Hong, Y, Deng, N, Jin, HN, Xuan, ZZ, Qian, YX, Wu, ZY, et al. Saikosaponin a modulates remodeling of Kv4.2-mediated a-type voltage-gated potassium currents in rat chronic temporal lobe epilepsy. Drug Des Devel Ther. (2018) 12:2945–58. doi: 10.2147/DDDT.S166408
10. Fan, Z, Feng, X, Fan, Z, Zhu, X, and Yin, S. Immunotherapy by targeting of VGKC complex for seizure control and prevention of cognitive impairment in a mouse model of epilepsy. Mol Med Rep. (2018) 18:169–78. doi: 10.3892/mmr.2018.9004
11. Kjær, C, Barzaghi, G, Bak, LK, Goetze, JP, Yde, CW, Woldbye, D, et al. Transcriptome analysis in patients with temporal lobe epilepsy. Brain. (2019) 142:e55. doi: 10.1093/brain/awz265
12. Langfelder, P, and Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. (2008) 9:559. doi: 10.1186/1471-2105-9-559
13. Joo, J, and Lebowitz, JL. Phys rev E behavior of susceptible-infected-susceptible epidemics on heterogeneous networks with saturation. Stat Nonlin Soft Matter Phys. (2004) 69:066105. doi: 10.1103/PhysRevE.69.066105
14. Gao, Y, Sun, B, Hu, J, Ren, H, Zhou, H, Chen, L, et al. Identification of gene modules associated with survival of diffuse large B-cell lymphoma treated with CHOP-based chemotherapy. Pharmacogenomics J. (2020) 20:705–16. doi: 10.1038/s41397-020-0161-6
15. Le, TT, Savitz, J, Suzuki, H, Misaki, M, Teague, TK, White, BC, et al. Identification and replication of RNA-Seq gene network modules associated with depression severity. Transl Psychiatry. (2018) 8. doi: 10.1038/s41398-018-0234-3
16. Jiang, Y, Leng, J, Lin, Q, and Zhou, F. Epithelial-mesenchymal transition related genes in unruptured aneurysms identified through weighted gene coexpression network analysis. Sci Rep. (2022) 12:225. doi: 10.1038/s41598-021-04390-6
17. Smyth, GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. (2004) 3:1–25. doi: 10.2202/1544-6115.1027
18. Yu, G, Wang, LG, Han, Y, and He, QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
19. Szklarczyk, D, Morris, JH, Cook, H, Kuhn, M, Wyder, S, Simonovic, M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. (2017) 45:D362–8. doi: 10.1093/nar/gkw937
20. Dweep, H, Gretz, N, and Sticht, C. miRWalk database for miRNA-target interactions. Methods Mol Biol. (2014) 1182:289–305. doi: 10.1007/978-1-4939-1062-5_25
21. Hornbeck, PV, Kornhauser, JM, Tkachev, S, Zhang, B, Skrzypek, E, Murray, B, et al. PhosphoSitePlus: a comprehensive resource for investigating the structure and function of experimentally determined post-translational modifications in man and mouse. Nucleic Acids Res. (2012) 40:D261–70. doi: 10.1093/nar/gkr1122
22. Livak, KJ, and Schmittgen, TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C (T)) method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.1262
23. Lörinc, E, Jakobsson, B, Landberg, G, and Veress, B. Ki67 and p 53 immunohistochemistry reduces interobserver variation in assessment of Barrett's oesophagus. Histopathology. (2005) 46:642–8. doi: 10.1111/j.1365-2559.2005.02139.x
24. Feng, D, Guo, YY, Wang, W, Yan, LF, Sun, T, Liu, QQ, et al. α-Subunit tyrosine phosphorylation is required for activation of the large conductance Ca2+-activated Potassium Channel in the rabbit sphincter of Oddi. Am J Pathol. (2022) 192:1725–44. doi: 10.1016/j.ajpath.2022.08.005
25. Nguyen, NH, and Brodsky, JL. The cellular pathways that maintain the quality control and transport of diverse potassium channels. Biochim Biophys Acta Gene Regul Mech. (2023) 1866:194908. doi: 10.1016/j.bbagrm.2023.194908
26. Scantlebury, MH, Chun, KC, Ma, SC, Rho, JM, and Kim, DY. Adrenocorticotropic hormone protects learning and memory function in epileptic KCNA1-null mice. Neurosci Lett. (2017) 645:14–8. doi: 10.1016/j.neulet.2017.02.069
27. Zhao, JY, Liang, L, Gu, X, Li, Z, Wu, S, Sun, L, et al. DNA methyltransferase DNMT3a contributes to neuropathic pain by repressing KCNA2 in primary afferent neurons. Nat Commun. (2017) 8:14712. doi: 10.1038/ncomms14712
28. Bowman, P, Day, J, Torrens, L, Shepherd, MH, Knight, BA, Ford, TJ, et al. Cognitive, neurological, and behavioral features in adults with KCNJ11 neonatal diabetes. Diabetes Care. (2019) 42:215–24. doi: 10.2337/dc18-1060
29. Toraman, A, Toraman, E, Özkaraca, M, and Budak, H. Increased nociceptive sensitivity is associated with periodontal inflammation and expression of chronic pain genes in gingival tissues of male rats. Chem Biol Interact. (2022) 366:110128. doi: 10.1016/j.cbi.2022.110128
30. Kong, C, Du, J, Bu, H, Huang, C, Xu, F, and Ren, H. Lnc RNA KCNA2-AS regulates spinal astrocyte activation through STAT3 to affect postherpetic neuralgia. Mol Med. (2020) 26:113. doi: 10.1186/s10020-020-00232-9
31. Akyuz, E, Polat, AK, Eroglu, E, Kullu, I, Angelopoulou, E, and Paudel, YN. Revisiting the role of neurotransmitters in epilepsy: an updated review. Life Sci. (2021) 265:118826. doi: 10.1016/j.lfs.2020.118826
32. Imbrici, P, Conte, E, Blunck, R, Stregapede, F, Liantonio, A, Tosi, M, et al. A novel KCNA2 variant in a patient with non-progressive congenital Ataxia and epilepsy: functional characterization and sensitivity to 4-Aminopyridine. Int J Mol Sci. (2021) 22:9913. doi: 10.3390/ijms22189913
33. Brew, HM, Hallows, JL, and Tempel, BL. Hyperexcitability and reduced low threshold potassium currents in auditory neurons of mice lacking the channel subunit Kv1. 1. J Physiol. (2003) 548:1–20. doi: 10.1113/jphysiol.2002.035568
34. Hsiao, CF, Kaur, G, Vong, A, Bawa, H, and Chandler, SH. Participation of Kv1 channels in control of membrane excitability and burst generation in mesencephalic V neurons. J Neurophysiol. (2009) 101:1407–18. doi: 10.1152/jn.91053.2008
35. Shu, Y, Yu, Y, Yang, J, and McCormick, DA. Selective control of cortical axonal spikes by a slowly inactivating K+ current. Proc Natl Acad Sci U S A. (2007) 104:11453–8. doi: 10.1073/pnas.0702041104
36. Eunson, LH, Rea, R, Zuberi, SM, Youroukos, S, Panayiotopoulos, CP, Liguori, R, et al. Clinical, genetic, and expression studies of mutations in the potassium channel gene KCNA1 reveal new phenotypic variability. Ann Neurol. (2000) 48:647–56. doi: 10.1002/1531-8249(200010)48:4<647::AID-ANA12>3.0.CO;2-Q
37. Zuberi, SM, Eunson, LH, Spauschus, A, De Silva, R, Tolmie, J, Wood, NW, et al. A novel mutation in the human voltate-gated potassium channel gene (Kv1. 1) associates with episodic ataxia type I and sometimes with partial epilepsy. Brain. (1999) 122:817–25. doi: 10.1093/brain/122.5.817
38. Fenoglio-Simeone, KA, Wilke, JC, Milligan, HL, Allen, CN, Rho, JM, and Maganti, RK. Ketogenic diet treatment abolishes seizure periodicity and improves diurnal rhythmicity in epileptic KCNA1-null mice. Epilepsia. (2009) 50:2027–34. doi: 10.1111/j.1528-1167.2009.02163.x
39. Wright, S, Wallace, E, Hwang, Y, and Maganti, R. Seizure phenotypes, periodicity, and sleep-wake pattern of seizures in Kcna-1 null mice. Epilepsy Behav. (2016) 55:24–9. doi: 10.1016/j.yebeh.2015.11.028
40. Iyer, SH, Matthews, SA, Simeone, TA, and Maganti, R. Kristina a Simeone accumulation of rest deficiency precedes sudden death of epileptic Kv1. 1 knockout mice, a model of sudden unexpected death in epilepsy. Epilepsia. (2018) 59:92–105. doi: 10.1111/epi.13953
41. Roundtree, HM, Simeone, TA, Johnson, C, Matthews, SA, Samson, KK, and Simeone, KA. Orexin receptor antagonism improves sleep and reduces seizures in KCNA1-null mice. Sleep. (2016) 39:357–68. doi: 10.5665/sleep.5444
42. Smart, SL, Lopantsev, V, Zhang, CL, Robbins, CA, Wang, H, Chiu, SY, et al. Deletion of the KV1. 1 Potassium Channel causes epilepsy in mice. Neuron. (1998) 20:809–19. doi: 10.1016/S0896-6273(00)81018-1
43. Glasscock, E, Yoo, JW, Chen, TT, Klassen, TL, and Noebels, JL. Kv1. 1 potassium channel deficiency reveals brain-driven cardiac dysfunction as a candidate mechanism for sudden unexplained death in epilepsy. J Neurosci. (2010) 30:5167–75. doi: 10.1523/JNEUROSCI.5591-09.2010
44. Robbins, CA, and Tempel, BL. Kv1. 1 and Kv1.2: Similar channels, different seizure models. Epilepsia. (2012) 53:134–41. doi: 10.1111/j.1528-1167.2012.03484.x
45. Persson, AS, Klement, G, Almgren, M, Sahlholm, K, Nilsson, J, Petersson, S, et al. A truncated Kv1.1 protein in the brain of the megencephaly mouse: expression and interaction. BMC Neurosci. (2005) 6:65. doi: 10.1186/1471-2202-6-65
46. Zuberi, SM, Eunson, LH, Spauschus, A, De Silva, R, Tolmie, J, Wood, NW, et al. A novel mutation in the human voltage-gated potassium channel gene (Kv1. 1) associates with episodic ataxia type 1 and sometimes with partial epilepsy. Brain. (1999) 122:817–25. doi: 10.1093/brain/122.5.817
47. Moore, BM, Jerry Jou, C, Tatalovic, M, Kaufman, ES, Kline, DD, and Kunze, DL. The Kv1. 1 null mouse, a model of sudden unexpected death in epilepsy (SUDEP). Epilepsia. (2014) 55:1808–16. doi: 10.1111/epi.12793
48. Simeone, KA, Matthews, SA, Rho, JM, and Simeone, TA. Ketogenic diet treatment increases longevity in KCNA1-null mice, a model of sudden unexpected death in epilepsy. Epilepsia. (2016) 57:e178–82. doi: 10.1111/epi.13444
49. Wykes, RC, Heeroma, JH, Mantoan, L, Zheng, K, Mac Donald, DC, Deisseroth, K, et al. Optogenetic and potassium channel gene therapy in a rodent model of focal neocortical epilepsy. Sci Transl Med. (2012) 4:161ra152. doi: 10.1126/scitranslmed.3004190
50. Snowball, A, Chabrol, E, Wykes, RC, Shekh-Ahmad, T, Cornford, JH, Lieb, A, et al. Epilepsy gene therapy using an engineered potassium channel. J Neurosci. (2019) 39:3159–69. doi: 10.1523/JNEUROSCI.1143-18.2019
51. Colasante, G, Qiu, Y, Massimino, L, Di Berardino, C, Cornford, JH, Snowball, A, et al. In vivo CRISPR a decreases seizures and rescues cognitive deficits in a rodent model of epilepsy. Brain. (2020) 143:891–905. doi: 10.1093/brain/awaa045
52. Trimmer, JS. Subcellular localization of K+ channels in mammalian brain neurons: remarkable precision in the midst of extraordinary complexity. Neuron. (2015) 85:238–56. doi: 10.1016/j.neuron.2014.12.042
53. Vacher, H, Mohapatra, DP, and Trimmer, JS. Localization and targeting of voltage-dependent ion channels in mammalian central neurons. Physiol Rev. (2008) 88:1407–47. doi: 10.1152/physrev.00002.2008
54. Willis, M, Leitner, I, Seppi, K, Trieb, M, Wietzorrek, G, Marksteiner, J, et al. Shaker-related voltage-gated potassium channels Kv1 in human hippocampus. Brain Struct Funct. (2018) 223:2663–71. doi: 10.1007/s00429-018-1653-x
55. Brew, HM, Gittelman, JX, Silverstein, RS, Hanks, TD, Demas, VP, Robinson, LC, et al. Seizures and reduced life span in mice lacking the potassium channel subunit Kv1.2, but hypoexcitability and enlarged Kv1 currents in auditory neurons. J Neurophysiol. (2007) 98:1501–25. doi: 10.1152/jn.00640.2006
56. Gong, P, Jiao, X, Zhang, Y, and Yang, Z. Generation of a human iPSC line from an epileptic encephalopathy patient with electrical status epilepticus during sleep carrying KCNA2 (p.P405L) mutation. Stem Cell Res. (2020) 49:102080. doi: 10.1016/j.scr.2020.102080
57. Uysal, B, Löffler, H, Rosa, F, Lerche, H, and Schwarz, N. Generation of an induced pluripotent stem cell (iPSC) line (HIHDNEi003-a) from a patient with developmental and epileptic encephalopathy carrying a KCNA2 (p. Thr 374Ala) mutation. Stem Cell Res. (2019) 40:101543. doi: 10.1016/j.scr.2019.101543
58. Hattersley, AT, and Ashcroft, FM. Activating mutations in Kir 6.2 and neonataldiabetes: new clinical syndromes, new scientific insights, and new therapy. Diabetes. (2005) 54:2503–13. doi: 10.2337/diabetes.54.9.2503
59. Foote, KM, Lyman, KA, Han, Y, Michailidis, IE, Heuermann, RJ, Mandikian, D, et al. Phosphorylation of the HCN channel auxiliary subunit TRIP8b is altered in an animal model of temporal lobe epilepsy and modulates channel function. J Biol Chem. (2019) 294:15743–58. doi: 10.1074/jbc.RA119.010027
60. Beesley, S, Sullenberger, T, Crotty, K, Ailani, R, D'Orio, C, Evans, K, et al. Kumar SS.D-serine mitigates cell loss associated with temporal lobe epilepsy. Nat Commun. (2020) 11:4966. doi: 10.1038/s41467-020-18757-2
61. Sosanya, NM, Brager, DH, Wolfe, S, Niere, F, and Raab-Graham, KF. Rapamycin reveals an mTOR-independent repression of Kv1.1 expression during epileptogenesis. Neurobiol Dis. (2015) 73:96–105. doi: 10.1016/j.nbd.2014.09.011
62. Lin, TW, Harward, SC, Huang, YZ, and McNamara, JO. Targeting BDNF/Trk B pathways for preventing or suppressing epilepsy. Neuropharmacology. (2020) 167:107734. doi: 10.1016/j.neuropharm.2019.107734
63. Jeong, KH, Kim, SH, Choi, YH, Cho, I, and Kim, WJ. Increased expression of WNK3 in dispersed granule cells in hippocampal sclerosis of mesial temporal lobe epilepsy patients. Epilepsy Res. (2018) 147:58–61. doi: 10.1016/j.eplepsyres.2018.09.006
64. Schuster, J, Laan, L, Klar, J, Jin, Z, Huss, M, Korol, S, et al. Transcriptomes of Dravet syndrome iPSC derived GABAergic cells reveal dysregulated pathways for chromatin remodeling and neurodevelopment. Neurobiol Dis. (2019) 132:104583. doi: 10.1016/j.nbd.2019.104583
65. Xiong, Q, Oviedo, HV, Trotman, LC, and Zador, AM. PTEN regulation of local and long range connections in mouse auditory cortex. J Neurosci. (2012) 32:1643–52. doi: 10.1523/JNEUROSCI.4480-11.2012
66. Xiang, L, Ren, Y, Li, X, Zhao, W, and Song, Y. MicroRNA-204 suppresses epileptiform discharges through regulating Trk B-ERK1/2-CREB signaling in cultured hippocampal neurons. Brain Res. (2016) 1639:99–107. doi: 10.1016/j.brainres.2016.02.045
67. Kanamaru, H, Kawakita, F, Nishikawa, H, Nakano, F, Asada, R, and Suzuki, H. Clarithromycin ameliorates early brain injury after subarachnoid hemorrhage via suppressing Periostin-related pathways in mice. Neurotherapeutics. (2021) 18:1880–90. doi: 10.1007/s13311-021-01050-5
68. Gu, Y, Jiang, J, and Liang, C. TFAP4 promotes the growth of prostate cancer cells by upregulating FOXK1. Exp Ther Med. (2021) 22:1299. doi: 10.3892/etm.2021.10734
69. Liu, JN, Kong, XS, Sun, P, Wang, R, Li, W, and Chen, QF. An integrated pan-cancer analysis of TFAP4 aberrations and the potential clinical implications for cancer immunity. J Cell Mol Med. (2021) 25:2082–97. doi: 10.1111/jcmm.16147
70. Srivastava, A, Dixit, AB, Paul, D, Tripathi, M, Sarkar, C, Chandra, PS, et al. Comparative analysis of cytokine/chemokine regulatory networks in patients with hippocampal sclerosis (HS) and focal cortical dysplasia (FCD). Sci Rep. (2017) 7:15904. doi: 10.1038/s41598-017-16041-w
71. Peng, G, Yuan, Y, He, Q, Wu, W, and Luo, BY. MicroRNA let-7e regulates the expression of caspase-3 during apoptosis of PC12 cells following anoxia/reoxygenation injury. Brain Res Bull. (2011) 86:272–6. doi: 10.1016/j.brainresbull.2011.07.017
72. Song, YJ, Tian, XB, Zhang, S, Zhang, YX, Li, X, Li, D, et al. Temporal lobe epilepsy induces differential expression of hippocampal miRNAs including let-7e and miR-23a/b. Brain Res. (2011) 1387:134–40. doi: 10.1016/j.brainres.2011.02.073
73. Xie, Y, Shi, H, and Han, B. Bioinformatic analysis of underlying mechanisms of Kawasaki disease via weighted gene correlation network analysis (WGCNA) and the least absolute shrinkage and selection operator method (LASSO) regression model. BMC Pediatr. (2023) 23:90. doi: 10.1186/s12887-023-03896-4
74. Xu, Z, and Chen, C. The downregulation of PTGS2 mediated by nc RNAs is tightly correlated with systemic sclerosis-interstitial lung disease. Front Genet. (2022) 12:795034. doi: 10.3389/fgene.2021.795034
75. Chen, Z, Feng, R, Kahlert, UD, Chen, Z, Torres-Dela Roche, LA, Soliman, A, et al. Construction of ce RNA networks associated with CD8 T cells in breast Cancer. Front Oncol. (2022) 12:883197. doi: 10.3389/fonc.2022.883197
76. Oshima, E, Urabe, N, Shingu, K, and Mori, K. Anticonvulsant actions of enflurane on epilepsy models in cats. Anesthesiology. (1985) 63:29–40. doi: 10.1097/00000542-198507000-00005
77. Jiang, Q, Wills, M, Geng, X, and Ding, Y. Chlorpromazine and promethazine reduces brain injury through RIP1-RIP3 regulated activation of NLRP3 inflammasome following ischemic stroke. Neurol Res. (2021) 43:668–76. doi: 10.1080/01616412.2021.1910904
78. Guo, S, Li, F, Wills, M, Yip, J, Wehbe, A, Peng, C, et al. Chlorpromazine and promethazine (C+P) reduce brain injury after ischemic stroke through the PKC-δ/NOX/MnSOD pathway. Mediat Inflamm. (2022) 2022:1–15. doi: 10.1155/2022/6886752
79. Gong, L, Zhu, T, Chen, C, Xia, N, Yao, Y, Ding, J, et al. Miconazole exerts disease-modifying effects during epilepsy by suppressing neuroinflammation via NF-κB pathway and iNOS production. Neurobiol Dis. (2022) 172:105823. doi: 10.1016/j.nbd.2022.105823
80. Gill, RS, Mirsattari, SM, and Leung, LS. Resting state functional network disruptions in a kainic acid model of temporal lobe epilepsy. Neuroimage Clin. (2016) 13:70–81. doi: 10.1016/j.nicl.2016.11.002
Keywords: temporal lobe epilepsy, WGCNA, potassium channel genes, target drugs, bioinformatics analyses
Citation: Zhang L-m, Chen L, Zhao Y-f, Duan W-m, Zhong L-m and Liu M-w (2023) Identification of key potassium channel genes of temporal lobe epilepsy by bioinformatics analyses and experimental verification. Front. Neurol. 14:1175007. doi: 10.3389/fneur.2023.1175007
Edited by:
Carlos De Cabo, Complejo Hospitalario Universitario de Albacete, SpainReviewed by:
Jun Hu, Fujian Medical University Union Hospital, ChinaYimin Chen, Foshan Sanshui District People's Hospital, China
Yu Yin, Kunming University of Science and Technology, China
Na He, The Second Affiliated Hospital of Guangzhou Medical University, China
Copyright © 2023 Zhang, Chen, Zhao, Duan, Zhong and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lian-mei Zhong, 13888967787@163.com; Ming-wei Liu, lmw2004210@163.com
†These authors have contributed equally to this work