Hippocampus RNA Sequencing of Pentylenetetrazole-Kindled Rats and Upon Treatment of Novel Chemical Q808

The expression of genes altered in epilepsy remains incomplete, particularly in the hippocampus, which exhibits exquisite vulnerability to epilepsy. Q808 is an innovation chemical compound that has potent anti-convulsant effect. Exploring its mechanism can not only explore the pathogenesis of epilepsy but also provide a theoretical basis for its clinical application. The present study aimed to use RNA sequencing (RNA-seq) to reveal the gene transcriptomic profile of chronic pentylenetetrazole (PTZ)-kindled seizure rats and the difference of the PTZ model rat before and after treatment with Q808. Quantitative real-time PCR (qRT-PCR) was performed to validate the RNA-seq results. The protein level was estimated with Western blot. Hippocampal transcriptomic analysis showed that 289 differentially expressed genes (DEGs) were confirmed in the PTZ-kindled seizure group compared with the vehicle control. Gene cluster analysis identified most of the DEGs linked to neuronal apoptosis, neurogenesis, neuronal projections, and neurotransmitter regulation. After analysis across the three groups, 23 hub genes and 21 pathways were identified, and qRT-PCR analysis confirmed that most of the mRNA levels of hub genes were consistent with the RNA-seq results. Q808 treatment increased the level of ACE, a GABA-related protein. Our analysis showed the comprehensive compendium of genes and pathways differentially expressed for PTZ-kindled seizure rats and upon Q808 treatment in PTZ-kindled seizure, which may provide a theoretical basis to explore the mechanism and unique efficacy of Q808 and the pathophysiology of epilepsy in the future.

Introduction

Epilepsy is a relatively common neurological disorder that affects people of all ages. Approximately 2–3% of the population worldwide have epilepsy (Gavvala and Schuele, 2016). The typical symptoms of epilepsy include recurrent unprovoked seizures and psychiatric symptoms, such as fear and anxiety, which considerably diminish the quality of life of a person (Thijs et al., 2019). The pathophysiology of epilepsy has various mechanisms, including neurotransmitter imbalance, channelopathies, neural migration, neuronal loss, brain injuries, degenerative disorders, morphological abnormalities, cortical and/or hippocampal malformations, and genetic background (Brodie et al., 2016; Pfisterer et al., 2020). Several anti-epileptic drugs are available in the market, but about one-third of the epileptic patients fail to achieve seizure control, and they develop refractory epilepsy (Löscher et al., 2013). In addition, traditional anti-epileptic drugs show various and serious adverse reactions (Stefan and Feuerstein, 2007). Therefore, identifying a new anti-epileptic agent is necessary.

Q808, 6-(4-chlorophenoxy)-tetrazolo[5,1-a]phthalazine (Figure 1), is an innovative chemical component with an international patent, which showed a broad-spectrum anti-epileptic activity in various experimental seizure models including a maximal electroshock seizure model and chemical-induced seizure model, such as induced by PTZ, 3-mercaptopropionic acid (3-MP), thiosemicarbazide (THIO), and isoniazid (ISO) (Sun et al., 2010). In our previous study, we found that Q808 can specifically increase the level of GABA, an inhibitory neurotransmitter, and enhance the frequency of spontaneous inhibitory postsynaptic currents in the hippocampus (Li et al., 2017). However, knowledge about the molecular mechanisms of Q808 is still limited.

FIGURE 1

FIGURE 1. Chemical structure of Q808.

The hippocampus is a heterogeneous structure with distinct regional connectivity along its longitudinal axis, and it plays a significant role in epileptogenesis and seizure maintenance (Strange et al., 2014; Huberfeld et al., 2015). The chemical kindling of seizures triggered by PTZ has been used to reflect the pathogenesis of human epilepsy, which is considered an appropriate model for drug-resistant epilepsy (Löscher, 2017). In the present study, hippocampal damage, such as neuronal apoptosis, has been observed after PTZ-induced seizure (Taskıran et al., 2020). This damage in the hippocampus can not only promote the occurrence and development of epilepsy but also cause impairment of cognitive function. Sections from PTZ-treated seizure mice showed morphological changes such as vacuolation and several degenerative areas with darkly stained nuclei in the hippocampus (El-Missiry et al., 2020). Apart from the change of neuro, some neurotransmitters also significantly change in the hippocampus of the PTZ model. In our previous study, inhibitory neurotransmitters such as GABA decreased evidently in the hippocampi of the PTZ-kindled model (Li et al., 2017). These changes in hippocampus could be associated with the genetic change. Considerable evidence has shown that the mRNA expression levels of the subunit of the GABAA receptor and apoptotic genes, namely, Bax and Caspase-3, were significantly increased (Ahmadirad et al., 2014; Singh et al., 2019). Thus, characterizing the gene expression profiles of the hippocampus is crucial to elucidate the molecular mechanisms of epileptogenesis and the anti-epileptic effect of Q808.

RNA-seq is a novel and sophisticated approach to comprehensively understand the molecular mechanisms of epilepsy in recent years (Chen et al., 2020). In the previous study, transcriptomic analyses have been performed with several seizure models, such as human intractable partial epilepsy (Guelfi et al., 2019), human temporal lobe epilepsy (Conte et al., 2020), drug refractory epilepsy patients (Dixit et al., 2016), and mouse and rat kainic acid models (Dong et al., 2020; O'Leary et al., 2020). An intraperitoneal continuous low-dose injection of PTZ into an animal could induce a severe tonic–clonic seizure, and this kindling model is widely applied to investigate the pathophysiology of epilepsy and epilepsy-related genes (Shimada and Yamagata, 2018). However, no studies have reported alterations in the hippocampal transcriptome of the PTZ-kindled seizure model before and after treatment with Q808.

In the current study, the chronic PTZ-kindled seizure model and seizure model treated with Q808 for 4 weeks were used, and the anticonvulsant effect of Q808 on the PTZ-kindling model was confirmed. In addition, RNA-seq was used to investigate the alternations in mRNA expression profile in the rat hippocampus. Furthermore, differential gene expression and level of a GABA-related protein were validated. The results presented in the study may benefit identifying novel mechanisms for epilepsy.

Materials and Methods

Chemicals

The Q808 compound was a kind gift from the Academy of Chinese Medical Sciences of Jilin Province, and the synthesis conformation of Q808 has been verified (Zhang et al., 2009). The other drugs and reagents, such as PTZ, Tween-80, and CMC-Na, were purchased from Sigma-Aldrich Chemical Company (St. Louis, United States). ACE, β-tubulin, and Goat Anti-Rabbit IgG (H + L) HRP antibody were purchased from Affinity Biosciences (Affinity, United States).

Experimental Animals

Experimental procedures were carried out in accordance with the National Institutes of Health guidelines and approval by the Animal Care and Use Committee of the College of Basic Medical Sciences in Jilin University (No. 2020-24). Thirty specific pathogen–free Wistar rats (180–220 g) were purchased from Changchun Yisi Experimental Animal Technology Co., Ltd. (Changchun, China; certificate number SYXK 2018-0001). The rats were housed three per cage under controlled temperature (22–25°C) and humidity (50–60%) with a 12°h light/dark cycle. Rats had free access to water and food. One week later, the animals were randomly divided into the vehicle control group (Con.), PTZ + vehicle group (PTZ), and PTZ + Q808 group (Q808). Rats in the PTZ group and Q808 group received intraperitoneal (i.p.) injection of sub-convulsive dose (35 mg/kg) of PTZ once every 48 h for the first 28 days (14 injections) and 60 mg/kg of PTZ on the last day of treatment to achieve a fully kindled state (Mahmoudi et al., 2020). The Con. group received equal amounts of normal saline. For the next 28 days, the Q808 group received Q808. The Con. group and PTZ group received the solvents of Q808 at the same time and same route as the Q808 group.

Behavioral Assessment of Seizure

The PTZ-kindled seizures were scored according to a modified Racine scale (Racine, 1972): category 0, no evidence of convulsive activity; category 1, whisker trembling and mouth and facial twitch; category 2: head nodding; category 3: forelimb clonus; category 4: forelimb rearing or jumping; and category 5: falling and tonic–clonic seizures. Rats with at least three consecutive seizures that scaled greater than or equal to Racine category 3 within the first 28 days were considered successfully kindled (Löscher, 2017) (success rate, 70%). The seizure latency to category 2 was also recorded.

Drug Treatment

Q808 drug solutions were freshly prepared with Tween-80 and 0.5% CMC-Na, and 30 mg/kg was administered to the rats via oral gavage daily for 28 days. The dose, route, and timing of administration of Q808 were based on pharmacokinetic considerations and preliminary experiments (Zhang et al., 2009; Sun et al., 2010; Li et al., 2017). The body weights of the experimental animals were measured every other day. On the next day of the last injection of Q808, animals in PTZ and Q808 groups were injected with a sub-convulsive dose of PTZ. Behavioral assessment of seizure was done as described above.

Hippocampal Tissue Preparation

The experimental animals were euthanized by decapitation under light ether anesthesia after the last seizure evaluation. Brains were removed and washed in cold phosphate-buffered saline, meninges were removed, and the hippocampus was bilaterally dissected on an ice platform. Individual samples were collected and stored in liquid nitrogen.

RNA Extraction and Library Preparation

Total RNA was extracted from the hippocampal tissue sample (about 100 mg) using the TRIzol kit (Invitrogen, United States) according to the manufacturer’s instructions. RNA purity was checked, and quantification was evaluated using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, United States). Samples with an absorbance (260/280) ratio between 1.8 and 2.2 were considered acceptable. An Agilent 2100 Bioanalyzer (Agilent Technologies, United States) was used to assess the RNA integrity. Then, the sequencing library was generated by using the MGIEasy Universal Library Convert Kit (MGI, China), VAHTS^® Universal V6 RNA-seq Library Prep Kit for Illumina (Vazyme, China), and VAHTS^® mRNA Capture Beads (Vazyme, China) according to the manufacturer’s protocols.

RNA Sequencing and Data Analysis

The libraries were sequenced on a BGI T7 platform. Raw reads were generated first and processed using the Trimmomatic software to remove reads containing adapter and ploy-N or low-quality reads. After this step, clean reads were obtained. By using Cufflinks (Trapnell et al., 2010), fragments per kilobase of transcript per million fragments mapped (FPKM) of each gene was calculated. HTSeq-count (Anders et al., 2015) was used to obtain the read counts of each gene. Differentially expressed gene (DEG) analysis was performed using the DESeq (2012) R package. The threshold for significantly differential expression in each sample was set at p value <0.05, and log2 of the fold-changes was filtered at > 1.5 or <0.67. In addition, hierarchical cluster analysis of DEGs across the three groups was performed to explore transcripts’ expression pattern.

Enrichment Analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses are often used to perform the functional annotation of DEGs. In detail, GO analysis was usually used to evaluate the biological process, cellular component, and molecular function of DEGs (Ashburner et al., 2000). The annotations were obtained from Gene Ontology (www.geneontology.org/) and NCBI (www.ncbi.nlm.nih.gov/). For identifying the pathways which DEGs are involved in, KEGG pathway analysis was used according to the KEGG database (www.genome.jp/kegg/) (Draghici et al., 2007). The two above-mentioned enrichment analyses were performed using R based on the hypergeometric distribution. A value of false discovery rate (FDR) less than 0.05 was considered significant.

qRT-PCR Validation

RNA utilized for qRT-PCR was extracted with TRIzol as described above. First-strand cDNA synthesis was performed using the Prime Script^TM RT Master Mix (Takara Biological Engineering Company, China). Real-time PCR amplifications were performed in the QuantStudio Real-Time PCR System with QuantStudio design and analysis software in a 20-μl reaction volume. Each reaction tube contained 2 μl cDNA, 10 μl SYBR Green Master Mix (Roche, Switzerland), 1.25 μmol.L⁻¹ forward and reverse primers, and RNase-free water. The cycling parameters were as follows: an initial hold stage of 95 °C for 10 min followed by 60 cycles of 95°C for 15 s and 59°C for 30 s. The target mRNA levels were normalized to the β-actin internal control gene. The fold-change in mRNA expression was analyzed according to the 2^−△△CT method. Primers for selected genes were designed using Primer-BLAST (Primer3 Input version 0.4.0 and BLAST) and are listed in Table 1.

TABLE 1

TABLE 1. Sequences of the primers.

Western Blot

Total protein was extracted from the hippocampal tissues with radio immunoprecipitation assay lysis buffer at 4°C for 30 min, and the total protein concentration was quantified by a BCA Protein Assay Kit (Solarbio, China). Equal amounts of protein were loaded on 8% sodium dodecyl sulfate–polyacrylamide gel and then transferred to polyvinylidene fluoride membranes with ice cooling. After blocking with 5% non-fat milk, the membranes were incubated with primary antibodies, rabbit Anti-ACE1 (1:1,000, Affinity), overnight at 4°C. After washing, membranes were soaked in the secondary antibodies, Anti-rabbit IgG HRP-linked (1:5,000, Affinity), for 1 h at room temperature. β-Tubulin was used as an internal control. Images were visualized by enhanced chemiluminescence using Tanon-4200SF (Biotanon) and analyzed using ImageJ software.

Statistical Analysis

The GraphPad Prism 7 software was used for data analysis. Data were expressed as mean ± SD. Statistical analysis was carried out by the Mann–Whitney U test and unpaired t-test to compare the difference in seizure score (a discrete variable) and seizure latency (a continuous variable), respectively, between two groups. One-way ANOVA followed by Tukey’s test was used to compare differences in hub gene and protein expression among three groups, and the experiment was performed with at least three biological replicates. An adjusted p-value < 0.05 was considered to be statistically significant.

Results

Racine Score and Seizure Latency

A sub-convulsant dose of PTZ was intraperitoneally injected once every 2 days for the first 28 days of the seizure model. Then, Q808 was given via oral gavage daily for the following 28 days. Behavioral changes were evaluated at the end of the drug treatment. PTZ-kindled rats have a significant reduction (p < 0.01) in seizure score after administration of Q808 (Figure 2A). The latency to seizure was significantly increased (p < 0.05) in the Q808-treated group compared with the PTZ model group (Figure 2B). All results indicated that Q808 showed an effective anti-convulsant activity against chronic PTZ-kindled seizure rats.

FIGURE 2

FIGURE 2. (A) Seizure score and (B) latency to seizures of PTZ-kindled rats before and after treatment with Q808. Data are shown as mean ± SD. Statistical analysis was carried out by the Mann–Whitney U test and unpaired t-test, respectively. n = 5, *p < 0.05, and **p < 0.01.

DEGs in the Hippocampus and Hierarchical Clustering

Total RNA was isolated from the hippocampus of rat and subjected to transcriptomic sequencing. The raw sequencing data can be found at GEO: GSE189785. The datasets of clean data are displayed in Supplementary Table S1 and showed that the medial clean data were 6.31 Gb per sample. In addition, the sequence alignment analysis indicated that the result of RNA sequencing was reliable because the ratio of clean reads was greater than 92.34% of Q30 and beyond 47.77% of the GC content. When compared to the Con. group, 289 genes were differentially expressed (82 down-regulated, 207 up-regulated) in the PTZ group (Figure 3A; Supplementary Table S2). Compared to the PTZ group, there were 152 genes differentially expressed in the Q808 group. Among them, 39 genes were down-regulated and 113 genes were up-regulated (Figure 3B; Supplementary Table S3). 23 out of 157 genes identified in the Con. and PTZ comparison were also identified in the Q808 and PTZ comparison and defined as hub genes (Table 2). The heatmap (Figure 3C) illustrates these alternate expression levels for all hub genes across the three groups. In the PTZ group vs. Con. group, the expression of 12 genes displayed an initial increase, but there was a decrease in the Q808 group vs. PTZ group; however, the expression of 11 genes exhibited a reduction in the PTZ group vs. Con. group, but there was a subsequent increase in the Q808 group vs. PTZ group. Then, the hierarchical clustering structure (Figure 3C) was determined based on the relative expression values of hub genes. The analysis showed that the Con. and PTZ groups had significantly different gene expression profiles, whereas the Q808 treatment group and Con. group had similar gene expression profiles. Therefore, after Q808 treatment, the gene expression levels of PTZ-kindled seizure rats could reverse to healthy controls.

FIGURE 3

FIGURE 3. Volcano plot and heatmaps of hub genes in the hippocampus. (A, B) Volcano plot showing up-regulated (red dots) and down-regulated (green dots) DEGs and normal expressed genes (gray dots). (C) Heatmap showing the expression level of 23 hub genes across the three groups. Blue and red colors represent low and high levels of expression, respectively. Corresponding hierarchical clustering analysis of these hub genes is shown in black bars.

TABLE 2

TABLE 2. Information of 23 hub genes identified by RNA-seq analysis across the three groups.

GO and KEGG Enrichment Analysis

GO analysis and KEGG enrichment analysis were conducted to explore the biological functions of the DEGs. In the GO analysis, the number of DEGs annotated to GO functions was identified, and the enrichment of GO functional significance was covered. Between the Con. and PTZ groups, 2080 GO terms were significantly enriched (Supplementary Table S4). After treatment with Q808, 1070 GO terms were significantly enriched compared with the PTZ model group without Q808 treatment (Supplementary Table S5). For both groups, the GO term of the most genes was cellular process, cell, and binding for the biological process, cellular component, and molecular function, respectively (Figures 4A,B). The detailed information of the gene functions is shown in Supplementary Table S6, which suggests that most of the genes are correlated with neuro-apoptosis, neurogenesis, neuronal projections, and neurotransmitter regulation, such as neurotransmitter biosynthetic process, secretion, transporter activity, and receptor activity.

FIGURE 4

FIGURE 4. (A) GO classification of DEGs between PTZ and Con. groups. (B) GO classification of DEGs between Q808 and PTZ groups. The GO term of the most genes was cellular process, cell, and binding for the biological process, cellular component, and molecular function, respectively.

In KEGG pathway analysis between the PTZ and Con. groups, Q808 and PTZ groups revealed 35 and 33 pathways, respectively (Supplementary Table S7; Supplementary Figures S1, S2). After comprehensive analysis across the three groups, 21 significantly enriched signaling pathways were identified (Table 3), five of which, including glycine, serine, and threonine metabolism; ether lipid metabolism; renin–angiotensin system; renin secretion; and malaria pathways, were associated with hub genes and were regarded as hub pathways. The hub genes on these hub pathways were Ace, Alas2, Enpp2, Gypc, and Hbb.

TABLE 3

TABLE 3. Twenty-one hub pathways across the three groups.

Validation of Candidate DEGs

qRT-PCR was conducted to verify the expression profiles obtained from RNA-seq and DEG analyses (Figure 5). We found that the expression levels of Ace, Col18a1, Enpp2, and Man2b2 down-regulated significantly in the PTZ model compared with the Con. group. The mRNA expression levels of Alsa2, Doc2a, Fh12, Gypc, Hbb, Mex3a, Mvd, Slc25a22, and Sox 11 were significantly up-regulated in the PTZ group compared with the Con. group. After Q808 treatment, the mRNA level of Ace, Ccdc88b, Col18a1, Enpp2, Man2b2, and Shox2 increased by > 5-fold; the gene expression of Alas2, Faap100, Itpripl1, Lmo4, and Wdr54 increased by one- to threefold; the expression of Doc2a, Fbl, Fhl2, Gypc, Mex3a, Mvd, Ovol2, Slc25a22, and Sox11 remarkably decreased. Based on the RNA-seq results between the PTZ and Con. groups, the expression levels of Lmo4 and Ovol2 were significantly increased. However, when analyzing the qPCR results, no significant changes were observed. Furthermore, we found opposite changes in the qRT-PCR results compared with the RNA-seq results of the gene expression of Faap100 between the PTZ and Con. groups and Fbl, Lmo4, and Wdr54 between the PTZ and Q808 groups. The expression of LOC102556346 and LOC103692570 could not be determined because of the limited information of primers.

FIGURE 5

FIGURE 5. Altered gene expression in the hippocampus. Differentially expressed mRNAs which were identified in the RNA-seq data were further validated by qRT-PCR according to the fold-changes of qRT-PCR (2^−ΔΔCt). *p < 0.05, **p < 0.01, and ***p < 0.001. Data are expressed as mean ± SD and analyzed using one-way ANOVA followed by Tukey’s test (n = 5).

Expression Level of GABA-Related Protein

According to the RNA-seq results and previous literature, one gene, namely, Ace, related to GABA release was identified, and the protein level was analyzed by Western blot (Figure 6). When compared with the Con. group, the expression levels of ACE in the PTZ group were significantly decreased (p < 0.05). Following Q808 treatment, the expression of ACE in the hippocampus significantly increased compared with that of the PTZ group (p < 0.05). The trend of ACE protein expression change was similar to the level of Ace mRNA expression.

FIGURE 6

FIGURE 6. The protein levels of ACE in the hippocampus of rats were determined by Western blot analysis. (A) Representative Western blot bands of ACE. (B) Quantitative analysis of ACE. β-Tubulin was used as the internal control. *p < 0.05. Data are expressed as mean ± SD and analyzed using one-way ANOVA followed by Tukey’s test (n = 3).

Discussion

This study showed the RNA-seq transcriptomic profile of the chronic PTZ-kindled seizure model in rats, which identified key transcriptomic signatures in PTZ-kindled seizure. When comparing Con. with the PTZ group, RNA-seq analysis of the hippocampal tissues identified 289 DEGs and revealed their functional enrichment. Experimental modeling and clinical studies have shown that seizures can cause neuronal death in the brain, and such brain injury may contribute to epileptogenesis (Henshall and Simon, 2005). After analysis of DEGs and their GO term, we identified some DEGs, such as Bad, Ccl2, Kcnip3, Hmox1, Rras2, Grm8, Serpinf1, and Fbxw7, which were associated with the cellular process such as regulation of neuronal apoptosis. Some of these genes have been reported, while some are not. Based on previous studies, Bad (Bcl2-associated agonist of cell death), an essential component of mitochondrial-dependent apoptosis, which increased sixfold in the PTZ group, has been reported in epilepsy by regulating cell death (Lu et al., 2021). The Grm8 (glutamate metabotropic receptor 8) agonist has anti-convulsant activities, and it is recognized as a potential target for novel AEDs (Chapman et al., 1999). Considering that the mechanisms for neuro-apoptosis are not fully clarified and its important role in epileptogenesis, other genes identified in our study related to apoptosis may provide a new insight.

Apart from neuro-apoptosis, increased neurogenesis and neuronal projections were observed in seizures, and reducing epileptic neurogenesis and excessive projections could serve as a therapeutic approach to limit the occurrence of seizures and progression of epilepsy (Mo et al., 2019). Based on our RNA-seq results, some of DEGs with GO terms including neurogenesis and neuronal projections have been reported to be associated with epilepsy. For example, increased Ace (Almeida et al., 2012), EphA10 (Xia et al., 2013), Limk1 (Jiang et al., 2020), and Sgk1 (Armas-Capote et al., 19912020) were associated with the seizure susceptibility. These genes were also significantly altered in our RNA-seq analysis, which could verify the accuracy of our sequencing results to a certain extent. We also detected a number of genes that were related to neurotransmitters between the PTZ and Con. groups. Imbalance of excitatory and inhibitory neurotransmitters may lead to epileptic seizures (Lascano et al., 2016). Expressional alterations of receptors and ion channels activated by neurotransmitters can lead to pathogenesis of epilepsy (Akyuz et al., 2021). Based on our results, the expression levels of neurotransmitter signaling Slc17a8, Stx1a, Th, Shisa8, Htr1f, and Grm8 were significantly changed in the PTZ group. Examining the function of these genes in epilepsy could be a crucial goal to explain the disruption of the excitatory and inhibitory balance in the brain. In addition, some DEGs identified in our research have not been reported in epilepsy so far, and further research on these genes may improve the identification of novel therapeutic targets for treating PTZ-kindled seizure.

Our pathway analysis of genes between the PTZ and Con. groups showed that 59% of DEGs were linked to signal transduction pathways. During epilepsy, signal transduction mediated by several neurotransmitters, including serotonin, glutamate, and GABA, and inflammation of the central nervous system (CNS) could affect seizure severity and lead to some comorbidities such as mood disorders (Rocha et al., 2014). The comprehensive understanding of the underlying mechanisms is still in its infancy. The DEGs clustered in these pathways may provide comprehensive understanding of the underlying mechanisms.

Furthermore, we analyzed the transcriptomic profile across the three groups, namely, Con., PTZ, and Q808. Q808 is an innovative anti-convulsant chemical, which is currently under preclinical study. To date, no studies have characterized the gene expression response to Q808. After RNA-seq analysis across the three groups, 23 hub genes and 21 pathways were identified. Hierarchical clustering analysis of hub genes showed that changes in gene expression caused by PTZ can be reversed after Q808 treatment. qRT-PCR validation of hub genes generally supports the results obtained from the RNA-seq experiments. Through several literature searches, we found 10 hub genes were related to epilepsy (Table 4). Thus, it has a great possibility that Q808 exerts anti-convulsant activity by affecting these genes. Further analysis indicated that one gene, Ace, might be involved in modulating GABA levels in the hippocampus, and it was mediated by Q808. ACE converts angiotensin (Ang) I into Ang II and subsequently into Ang-(1–7) by ACE2 (Ramos, 19792021). Ang-(1–7)-induced GABA release in the striatum has been verified by Stragier’s team (Stragier et al., 2005). However, the effect of Ace on the GABA release in the hippocampus has not been reported. ACE is also one of the key components of the renin–angiotensin system (RAS) (Krasniqi and Daci, 2019), and the RAS is one of the hub pathways identified in this research. The RAS can modulate different neurotransmitter systems, thereby controlling the excitability of the brain (Albrecht et al., 1997). Considering that Q808 can increase the level of GABA (Li et al., 2017), ACE or RAS might be one of the targets of Q808. Our future research will focus on the effect of ACE or RAS on GABAergic interneurons, especially GABA release, in the epilepsy model, and examine its molecular mechanism. In addition, evaluating the long-time safety of AEDs is of great importance (Gaitatzis and Sander, 2013). Thus, a future research study to investigate the impact of long-term treatment with Q808 on the genes and pathways in healthy animals could be an attractive study and is indispensable for Q808 clinical use.

TABLE 4

TABLE 4. Hub genes related to epilepsy.

In summary, our study is unique as it utilizes RNA‐seq to assess the gene expression in the hippocampus among normal control rats, chronic PTZ-kindled seizure rats, and PTZ model rats after treatment with Q808. The transcriptomic changes influenced by Q808 in the hippocampus are evident. 23 hub genes and 21 hub pathways across the three groups are identified, which confirmed prior studies and provided additional potential mechanisms of epileptogenesis and therapeutic targets for Q808.

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 at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE189785.

Ethics Statement

The animal study was reviewed and approved by the Animal Care and Use Committee of the College of Basic Medical Sciences in Jilin University.

Author Contributions

XL, XC, and WL conceptualized the research idea. XL, XC, WL, QW, D-WZ, DZ, S-WZ, and ZW were involved in data investigation, XL wrote the original draft and was involved in funding acquisition. XC and WL supervised the study and reviewed and edited the paper. QW, D-WZ, DZ, and ZW were responsible for data collection and statistical analysis. S-WZ provided analysis support. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by the Education Department of Jilin Province (Grant number: JJKH20211224KJ).

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/fphar.2022.820508/full#supplementary-material

References

Ackley, B. D., Crew, J. R., Elamaa, H., Pihlajaniemi, T., Kuo, C. J., and Kramer, J. M. (2001). The NC1/endostatin Domain of Caenorhabditis elegans Type XVIII Collagen Affects Cell Migration and Axon Guidance. J. Cel Biol 152 (6), 1219–1232. doi:10.1083/jcb.152.6.1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmadirad, N., Shojaei, A., Javan, M., Pourgholami, M. H., and Mirnajafi-Zadeh, J. (2014). Effect of Minocycline on Pentylenetetrazol-Induced Chemical Kindled Seizures in Mice. Neurol. Sci. 35 (4), 571–576. doi:10.1007/s10072-013-1552-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Akyuz, E., Polat, A. K., Eroglu, E., Kullu, I., Angelopoulou, E., and Paudel, Y. N. (2021). Revisiting the Role of Neurotransmitters in Epilepsy: An Updated Review. Life Sci. 265, 118826. doi:10.1016/j.lfs.2020.118826

PubMed Abstract | CrossRef Full Text | Google Scholar

Albrecht, D., Broser, M., and Krüger, H. (1997). Excitatory Action of Angiotensins II and IV on Hippocampal Neuronal Activity in Urethane Anesthetized Rats. Regul. Pept. 70 (2-3), 105–109. doi:10.1016/s0167-0115(97)00015-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Almeida, S. S., Naffah-Mazzacoratti, M. G., Guimarães, P. B., Wasinski, F., Pereira, F. E., Canzian, M., et al. (2012). Carbamazepine Inhibits Angiotensin I-Converting Enzyme, Linking it to the Pathogenesis of Temporal Lobe Epilepsy. Transl Psychiatry 2 (3), e93. doi:10.1038/tp.2012.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq--a Python Framework to Work with High-Throughput Sequencing Data. Bioinformatics 31 (2), 166–169. doi:10.1093/bioinformatics/btu638

PubMed Abstract | CrossRef Full Text | Google Scholar

Armas-Capote, N., Maglio, L. E., Pérez-Atencio, L., Martin-Batista, E., Reboreda, A., Barios, J. A., et al. (19912020). SGK1.1 Reduces Kainic Acid-Induced Seizure Severity and Leads to Rapid Termination of Seizures. Cereb. Cortex 30 (5), 3184–3197. doi:10.1093/cercor/bhz302

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: Tool for the Unification of Biology. The Gene Ontology Consortium. Nat. Genet. 25 (1), 25–29. doi:10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodie, M. J., Besag, F., Ettinger, A. B., Mula, M., Gobbi, G., Comai, S., et al. (2016). Epilepsy, Antiepileptic Drugs, and Aggression: An Evidence-Based Review. Pharmacol. Rev. 68 (3), 563–602. doi:10.1124/pr.115.012021

PubMed Abstract | CrossRef Full Text | Google Scholar

Chapman, A. G., Nanan, K., Yip, P., and Meldrum, B. S. (1999). Anticonvulsant Activity of a Metabotropic Glutamate Receptor 8 Preferential Agonist, (R,S)-4-phosphonophenylglycine. Eur. J. Pharmacol. 383 (1), 23–27. doi:10.1016/s0014-2999(99)00615-9

CrossRef Full Text | Google Scholar

Chen, Q. L., Xia, L., Zhong, S. P., Wang, Q., Ding, J., and Wang, X. (2020). Bioinformatic Analysis Identifies Key Transcriptome Signatures in Temporal Lobe Epilepsy. CNS Neurosci. Ther. 26 (12), 1266–1277. doi:10.1111/cns.13470

PubMed Abstract | CrossRef Full Text | Google Scholar

Conte, G., Parras, A., Alves, M., Ollà, I., De Diego-Garcia, L., Beamer, E., et al. (2020). High Concordance between Hippocampal Transcriptome of the Mouse Intra-amygdala Kainic Acid Model and Human Temporal Lobe Epilepsy. Epilepsia 61 (12), 2795–2810. doi:10.1111/epi.16714

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixit, A. B., Banerjee, J., Srivastava, A., Tripathi, M., Sarkar, C., Kakkar, A., et al. (2016). RNA-seq Analysis of Hippocampal Tissues Reveals Novel Candidate Genes for Drug Refractory Epilepsy in Patients with MTLE-HS. Genomics 107 (5), 178–188. doi:10.1016/j.ygeno.2016.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, X., Hao, X., Xu, P., Fan, M., Wang, X., Huang, X., et al. (2020). RNA Sequencing Analysis of Cortex and hippocampus in a Kainic Acid Rat Model of Temporal Lobe Epilepsy to Identify Mechanisms and Therapeutic Targets Related to Inflammation, Immunity and Cognition. Int. Immunopharmacol 87, 106825. doi:10.1016/j.intimp.2020.106825

PubMed Abstract | CrossRef Full Text | Google Scholar

Draghici, S., Khatri, P., Tarca, A. L., Amin, K., Done, A., Voichita, C., et al. (2007). A Systems Biology Approach for Pathway Level Analysis. Genome Res. 17 (10), 1537–1545. doi:10.1101/gr.6202607

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Missiry, M. A., Othman, A. I., Amer, M. A., Sedki, M., Ali, S. M., and El-Sherbiny, I. M. (2020). Nanoformulated Ellagic Acid Ameliorates Pentylenetetrazol-Induced Experimental Epileptic Seizures by Modulating Oxidative Stress, Inflammatory Cytokines and Apoptosis in the Brains of Male Mice. Metab. Brain Dis. 35 (2), 385–399. doi:10.1007/s11011-019-00502-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Faigle, R., and Song, H. (2013). Signaling Mechanisms Regulating Adult Neural Stem Cells and Neurogenesis. Biochim. Biophys. Acta 1830 (2), 2435–2448. doi:10.1016/j.bbagen.2012.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaitatzis, A., and Sander, J. W. (2013). The Long-Term Safety of Antiepileptic Drugs. CNS drugs 27 (6), 435–455. doi:10.1007/s40263-013-0063-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavvala, J. R., and Schuele, S. U. (2016). New-Onset Seizure in Adults and Adolescents: A Review. Jama 316 (24), 2657–2668. doi:10.1001/jama.2016.18625

PubMed Abstract | CrossRef Full Text | Google Scholar

Guelfi, S., Botia, J. A., Thom, M., Ramasamy, A., Perona, M., Stanyer, L., et al. (2019). Transcriptomic and Genetic Analyses Reveal Potential Causal Drivers for Intractable Partial Epilepsy. Brain 142 (6), 1616–1630. doi:10.1093/brain/awz074

PubMed Abstract | CrossRef Full Text | Google Scholar

Henshall, D. C., and Simon, R. P. (2005). Epilepsy and Apoptosis Pathways. J. Cereb. Blood Flow Metab. 25 (12), 1557–1572. doi:10.1038/sj.jcbfm.9600149

CrossRef Full Text | Google Scholar

Hu, X., Tang, J., Lan, X., and Mi, X. (2019). Increased Expression of DOC2A in Human and Rat Temporal Lobe Epilepsy. Epilepsy Res. 151, 78–84. doi:10.1016/j.eplepsyres.2019.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Huberfeld, G., Blauwblomme, T., and Miles, R. (2015). Hippocampus and Epilepsy: Findings from Human Tissues. Rev. Neurol. (Paris) 171 (3), 236–251. doi:10.1016/j.neurol.2015.01.563

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamali, S., Bartolomei, F., Robaglia-Schlupp, A., Massacrier, A., Peragut, J. C., Régis, J., et al. (2006). Large-scale Expression Study of Human Mesial Temporal Lobe Epilepsy: Evidence for Dysregulation of the Neurotransmission and Complement Systems in the Entorhinal Cortex. Brain 129 (Pt 3), 625–641. doi:10.1093/brain/awl001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Q., Tang, G., Fu, J., Yang, J., Xu, T., Tan, C. H., et al. (2020). Lim Kinase1 Regulates Seizure Activity via Modulating Actin Dynamics. Neurosci. Lett. 729, 134936. doi:10.1016/j.neulet.2020.134936

PubMed Abstract | CrossRef Full Text | Google Scholar

Krasniqi, S., and Daci, A. (2019). Role of the Angiotensin Pathway and its Target Therapy in Epilepsy Management. Int. J. Mol. Sci. 20 (3). doi:10.3390/ijms20030726

PubMed Abstract | CrossRef Full Text | Google Scholar

Lascano, A. M., Korff, C. M., and Picard, F. (2016). Seizures and Epilepsies Due to Channelopathies and Neurotransmitter Receptor Dysfunction: A Parallel between Genetic and Immune Aspects. Mol. Syndromol 7 (4), 197–209. doi:10.1159/000447707

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Zhang, H. J., Wang, Q., Zhang, D. W., Wu, D., Li, W., et al. (2017). Current Study of the Mechanism of Action of the Potential Anti-epileptic Agent Q808. Molecules 22 (7). doi:10.3390/molecules22071134

PubMed Abstract | CrossRef Full Text | Google Scholar

Löscher, W. (2017). Animal Models of Seizures and Epilepsy: Past, Present, and Future Role for the Discovery of Antiseizure Drugs. Neurochem. Res. 42 (7), 1873–1888. doi:10.1007/s11064-017-2222-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Löscher, W., Klitgaard, H., Twyman, R. E., and Schmidt, D. (2013). New Avenues for Anti-epileptic Drug Discovery and Development. Nat. Rev. Drug Discov. 12 (10), 757–776. doi:10.1038/nrd4126

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Z., Tang, M., Zhang, M., Li, Y., Shi, F., Zhan, F., et al. (2021). Expression and Functional Analysis of the BCL2-Associated Agonist of Cell Death (BAD) Gene in Grass Carp (Ctenopharyngodon Idella) during Bacterial Infection. Dev. Comp. Immunol. 123, 104160. doi:10.1016/j.dci.2021.104160

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackay, D. R., Hu, M., Li, B., Rhéaume, C., and Dai, X. (2006). The Mouse Ovol2 Gene Is Required for Cranial Neural Tube Development. Dev. Biol. 291 (1), 38–52. doi:10.1016/j.ydbio.2005.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahmoudi, T., Lorigooini, Z., Rafieian-Kopaei, M., Arabi, M., Rabiei, Z., Bijad, E., et al. (2020). Effect of Curcuma Zedoaria Hydro-Alcoholic Extract on Learning, Memory Deficits and Oxidative Damage of Brain Tissue Following Seizures Induced by Pentylenetetrazole in Rat. Behav. Brain Funct. 16 (1), 7. doi:10.1186/s12993-020-00169-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Mo, M., Eyo, U. B., Xie, M., Peng, J., Bosco, D. B., Umpierre, A. D., et al. (2019). Microglial P2Y12 Receptor Regulates Seizure-Induced Neurogenesis and Immature Neuronal Projections. J. Neurosci. 39 (47), 9453–9464. doi:10.1523/JNEUROSCI.0487-19.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Molinari, F., Raas-Rothschild, A., Rio, M., Fiermonte, G., Encha-Razavi, F., Palmieri, L., et al. (2005). Impaired Mitochondrial Glutamate Transport in Autosomal Recessive Neonatal Myoclonic Epilepsy. Am. J. Hum. Genet. 76 (2), 334–339. doi:10.1086/427564

CrossRef Full Text | Google Scholar

O'Leary, H., Vanderlinden, L., Southard, L., Castano, A., Saba, L. M., and Benke, T. A. (2020). Transcriptome Analysis of Rat Dorsal Hippocampal CA1 after an Early Life Seizure Induced by Kainic Acid. Epilepsy Res. 161, 106283. doi:10.1016/j.eplepsyres.2020.106283

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfisterer, U., Petukhov, V., Demharter, S., Meichsner, J., Thompson, J. J., Batiuk, M. Y., et al. (2020). Identification of Epilepsy-Associated Neuronal Subtypes and Gene Expression Underlying Epileptogenesis. Nat. Commun. 11 (1), 5038. doi:10.1038/s41467-020-18752-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, Z., Zhou, X., Gomez-Smith, M., Pandey, N. R., Lee, K. F., Lagace, D. C., et al. (2012). LIM Domain Only 4 (LMO4) Regulates Calcium-Induced Calcium Release and Synaptic Plasticity in the hippocampus. J. Neurosci. 32 (12), 4271–4283. doi:10.1523/JNEUROSCI.6271-11.2012

CrossRef Full Text | Google Scholar

Racine, R. J. (1972). Modification of Seizure Activity by Electrical Stimulation. II. Motor Seizure. Electroencephalogr Clin. Neurophysiol. 32 (3), 281–294. doi:10.1016/0013-4694(72)90177-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramos, A. J. (19792021). Brain Angiotensin System: a New Promise in the Management of Epilepsy? Clin. Sci. (London, Engl. 135 (6), 725–730.

Google Scholar

Rocha, L., Alonso-Vanegas, M., Orozco-Suárez, S., Alcántara-González, D., Cruzblanca, H., and Castro, E. (2014). Do certain Signal Transduction Mechanisms Explain the Comorbidity of Epilepsy and Mood Disorders? Epilepsy Behav. 38, 25–31. doi:10.1016/j.yebeh.2014.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimada, T., and Yamagata, K. (2018). Pentylenetetrazole-Induced Kindling Mouse Model. J. Vis. Exp. 136. doi:10.3791/56573

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, N., Saha, L., Kumari, P., Singh, J., Bhatia, A., Banerjee, D., et al. (2019). Effect of Dimethyl Fumarate on Neuroinflammation and Apoptosis in Pentylenetetrazol Kindling Model in Rats. Brain Res. Bull. 144, 233–245. doi:10.1016/j.brainresbull.2018.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Stefan, H., and Feuerstein, T. J. (2007). Novel Anticonvulsant Drugs. Pharmacol. Ther. 113 (1), 165–183. doi:10.1016/j.pharmthera.2006.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Stragier, B., Hristova, I., Sarre, S., Ebinger, G., and Michotte, Y. (2005). In Vivo characterization of the Angiotensin-(1-7)-Induced Dopamine and Gamma-Aminobutyric Acid Release in the Striatum of the Rat. Eur. J. Neurosci. 22 (3), 658–664. doi:10.1111/j.1460-9568.2005.04188.x

CrossRef Full Text | Google Scholar

Strange, B. A., Witter, M. P., Lein, E. S., and Moser, E. I. (2014). Functional Organization of the Hippocampal Longitudinal axis. Nat. Rev. Neurosci. 15 (10), 655–669. doi:10.1038/nrn3785

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X. Y., Wei, C. X., Deng, X. Q., Sun, Z. G., and Quan, Z. S. (2010). Evaluation of the Anticonvulsant Activity of 6-(4-Chlorophenyoxy)-Tetrazolo[5,1-A]phthalazine in Various Experimental Seizure Models in Mice. Pharmacol. Rep. 62 (2), 273–277. doi:10.1016/s1734-1140(10)70266-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Taskıran, A. S., Ozdemir, E., Gumus, E., and Ergul, M. (2020). The Effects of salmon Calcitonin on Epileptic Seizures, Epileptogenesis, and Postseizure Hippocampal Neuronal Damage in Pentylenetetrazole-Induced Epilepsy Model in Rats. Epilepsy Behav. 113, 107501. doi:10.1016/j.yebeh.2020.107501

PubMed Abstract | CrossRef Full Text | Google Scholar

Thijs, R. D., Surges, R., O'Brien, T. J., and Sander, J. W. (2019). Epilepsy in Adults. Lancet 393 (10172), 689–701. doi:10.1016/S0140-6736(18)32596-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation. Nat. Biotechnol. 28 (5), 511–515. doi:10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, Y., Luo, C., Dai, S., and Yao, D. (2013). Increased EphA/ephrinA Expression in hippocampus of Pilocarpine Treated Mouse. Epilepsy Res. 105 (1-2), 20–29. doi:10.1016/j.eplepsyres.2013.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, D., Febbo, I. G., Maroteaux, M. J., Wang, H., Song, Y., Han, X., et al. (19912021). The Transcription Factor Shox2 Shapes Neuron Firing Properties and Suppresses Seizures by Regulation of Key Ion Channels in Thalamocortical Neurons. Cereb. Cortex (N Y, N.Y. 31 (7), 3194–3212.

Google Scholar

Zhang, L., Guan, L. P., Sun, X. Y., Wei, C. X., Chai, K. Y., and Quan, Z. S. (2009). Synthesis and Anticonvulsant Activity of 6-Alkoxy-[1,2,4]triazolo[3,4-A]phthalazines. Chem. Biol. Drug Des. 73 (3), 313–319. doi:10.1111/j.1747-0285.2009.00776.x

PubMed Abstract | CrossRef Full Text | Google Scholar