Comprehensive identification of essential pathways and transcription factors related to epilepsy by gene set enrichment analysis on microarray datasets

Hippocampal transcriptomic analyses reveal the potential antiapoptotic mechanism of a novel anticonvulsant agent Q808 on pentylenetetrazol-induced epilepsy in rats

Brain apoptosis is one of the main causes of epileptogenesis. The antiapoptotic effect and potential mechanism of Q808, an innovative anticonvulsant chemical, have never been reported. In this study, the seizure stage and latency to reach stage 2 of pentylenetetrazol (PTZ) seizure rat model treated with Q808 were investigated. The morphological change and neuronal apoptosis in the hippocampus were detected by hematoxylin and eosin (HE) and terminal deoxynucleotidyl transferase-mediated dUTP nick end-labeling (TUNEL) staining, respectively. The hippocampal transcriptomic changes were observed using RNA sequencing (RNA-seq). The expression levels of hub genes were verified by quantitative reverse-transcription PCR (qRT-PCR). Results revealed that Q808 could allay the seizure score and prolong the stage 2 latency in seizure rats. The morphological changes of neurons and the number of apoptotic cells in the DG area were diminished by Q808 treatment. RNA-seq analysis revealed eight hub genes, including Map2k3, Nfs1, Chchd4, Hdac6, Siglec5, Slc35d3, Entpd1, and LOC103690108, and nine hub pathways among the control, PTZ, and Q808 groups. Hub gene Nfs1 was involved in the hub pathway sulfur relay system, and Map2k3 was involved in the eight remaining hub pathways, including Amyotrophic lateral sclerosis, Cellular senescence, Fc epsilon RI signaling pathway, GnRH signaling pathway, Influenza A, Rap1 signaling pathway, TNF signaling pathway, and Toll-like receptor signaling pathway. qRT-PCR confirmed that the mRNA levels of these hub genes were consistent with the RNA-seq results. Our findings might contribute to further studies exploring the new apoptosis mechanism and actions of Q808.

Potentials of miR-9-5p in promoting epileptic seizure and improving survival of glioma patients

Background

Epilepsy affects over 70 million people worldwide; however, the underlying mechanisms remain unclear. MicroRNAs (miRNAs) have essential functions in epilepsy. miRNA-9, a brain-specific/enriched miRNA, plays a role in various nervous system diseases and tumors, but whether miRNA-9 is involved in epilepsy and glioma-associated epilepsy remains unknown. Therefore, we aimed to explore the potential role of miR-9-5p in seizures and its effect on the survival of glioma patients, in order to provide new targets for the treatment of epilepsy and glioma.

Methods

The YM500v2 database was used to validate the expression of hsa-miR-9-5p in tissues. Moreover, qRT-PCR was performed to investigate the expression of miR-9-5p in temporal lobe epilepsy patients and rats with lithium-pilocarpine-induced seizures. Recombinant adeno-associated virus containing miR-9-5p was constructed to overexpress miR-9-5p in vivo. The effects of miR-9-5p on the behavior and electroencephalographic activities of the lithium-pilocarpine rat model of epilepsy were tested. Bioinformatics analysis was used to predict the targets of miR-9-5p and explore its potential role in epilepsy and glioma-associated epilepsy.

Results

The expression of miR-9-5p increased at 6 h and 7 days after lithium-pilocarpine-induced seizures in rats. Overexpression of miR-9-5p significantly shortened the latency of seizures and increased seizure intensity at 10 min and 20 min after administration of pilocarpine (P < 0.05). Predicted targets of miR-9-5p were abundant and enriched in the brain, and affected various pathways related to epilepsy and tumor. Survival analysis revealed that overexpression of miR-9-5p significantly improved the survival of patients from with low-grade gliomas and glioblastomas. The involvement of miR-9-5p in the glioma-associated epileptic seizures and the improvement of glioma survival may be related to multiple pathways, including the Rho GTPases and hub genes included SH3PXD2B, ARF6, and ANK2.

Conclusions

miR-9-5p may play a key role in promoting epileptic seizures and improving glioma survival, probably through multiple pathways, including GTPases of the Rho family and hub genes including SH3PXD2B, ARF6 and ANK2. Understanding the roles of miR-9-5p in epilepsy and glioma and the underlying mechanisms may provide a theoretical basis for the diagnosis and treatment of patients with epilepsy and glioma.

Whole-exome sequencing with targeted analysis and epilepsy after acute symptomatic neonatal seizures

BACKGROUND

The contribution of pathogenic gene variants with development of epilepsy after acute symptomatic neonatal seizures is not known.

METHODS

Case-control study of 20 trios in children with a history of acute symptomatic neonatal seizures: 10 with and 10 without post-neonatal epilepsy. We performed whole-exome sequencing (WES) and identified pathogenic de novo, transmitted, and non-transmitted variants from established and candidate epilepsy association genes and correlated prevalence of these variants with epilepsy outcomes. We performed a sensitivity analysis with genes associated with coronary artery disease (CAD). We analyzed variants throughout the exome to evaluate for differential enrichment of functional properties using exploratory KEGG searches.

RESULTS

Querying 200 established and candidate epilepsy genes, pathogenic variants were identified in 5 children with post-neonatal epilepsy yet in only 1 child without subsequent epilepsy. There was no difference in the number of trios with non-transmitted pathogenic variants in epilepsy or CAD genes. An exploratory KEGG analysis demonstrated a relative enrichment in cell death pathways in children without subsequent epilepsy.

CONCLUSIONS

In this pilot study, children with epilepsy after acute symptomatic neonatal seizures had a higher prevalence of coding variants with a targeted epilepsy gene sequencing analysis compared to those patients without subsequent epilepsy.

IMPACT

We performed whole-exome sequencing (WES) in 20 trios, including 10 children with epilepsy and 10 without epilepsy, both after acute symptomatic neonatal seizures. Children with post-neonatal epilepsy had a higher burden of pathogenic variants in epilepsy-associated genes compared to those without post-neonatal epilepsy. Future studies evaluating this association may lead to a better understanding of the risk of epilepsy after acute symptomatic neonatal seizures and elucidate molecular pathways that are dysregulated after brain injury and implicated in epileptogenesis.

Molecular tools for the characterization of seizure susceptibility in genetic rodent models of epilepsy

Epilepsy is a chronic neurological disorder characterized by abnormal neuronal activity that arises from imbalances between excitatory and inhibitory synapses, which are highly correlated to functional and structural changes in specific brain regions. The difference between the normal and the epileptic brain may harbor genetic alterations, gene expression changes, and/or protein alterations in the epileptogenic nucleus. It is becoming increasingly clear that such differences contribute to the development of distinct epilepsy phenotypes. The current major challenges in epilepsy research include understanding the disease progression and clarifying epilepsy classifications by searching for novel molecular biomarkers. Thus, the application of molecular techniques to carry out comprehensive studies at deoxyribonucleic acid, messenger ribonucleic acid, and protein levels is of utmost importance to elucidate molecular dysregulations in the epileptic brain. The present review focused on the great diversity of technical approaches available and new research methodology, which are already being used to study molecular alterations underlying epilepsy. We have grouped the different techniques according to each step in the flow of information from DNA to RNA to proteins, and illustrated with specific examples in animal models of epilepsy, some of which are our own. Separately and collectively, the genomic and proteomic techniques, each with its own strengths and limitations, provide valuable information on molecular mechanisms underlying seizure susceptibility and regulation of neuronal excitability. Determining the molecular differences between genetic rodent models of epilepsy and their wild-type counterparts might be a key in determining mechanisms of seizure susceptibility and epileptogenesis as well as the discovery and development of novel antiepileptic agents. This article is part of the Special Issue "NEWroscience 2018".

Differentially expressed circRNA and functional pathways in the hippocampus of epileptic mice based on next-generation sequencing

Epilepsy is a clinical syndrome caused by the highly synchronized abnormal discharge of brain neurons. It has the characteristics of paroxysmal, transient, repetitive, and stereotyped. Circular RNAs (circRNAs) are a recently discovered type of noncoding RNA with diverse cellular functions related to their excellent stability; additionally, some circRNAs can bind and regulate microRNAs (miRNAs). The present study was designed to screen the differentially expressed circRNA in an acute seizure model of epilepsy in mice, analyze the related miRNA and mRNA, and study their participating functions and enrichment pathways. In order to obtain the differential expression of circRNA in epilepsy and infer their function, we used next-generation sequencing and found significantly different transcripts. CIRI (circRNA identifier) software was used to predict circRNA from the hippocampus cDNA, EdgeR was applied for the differential circRNA analysis between samples, and Cytoscape 3.7.2 software was used to draw the network diagram. A total of 10,388 differentially expressed circRNAs were identified, of which 34 were upregulated and 66 were downregulated. Among them, mm9_circ_008777 and mm9_circ_004424 were the key upregulated genes, and their expression in the epilepsy group was verified using Quantitative real-time PCR (QPCR). The analysis indicated that the extracted gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were closely related to several epilepsy-associated processes. This study determined that mm9_circ_008777 and mm9_circ_004424 are potential biomarkers of epilepsy, which play important roles in epilepsy-related pathways. These results could help improve the understanding of the biological mechanisms of circRNAs and epilepsy treatments.

Functional Genomics of Epileptogenesis in Animal Models and Humans

It has been estimated that epilepsies are among the top five neurological diseases with the highest burden of disease. In recent years, genome-wide expression studies (GWES) have been carried out in experimental models of epilepsy and in samples from human patients. In this study, I carried out meta-analyses and analyses of convergence for available GWES for epileptogenesis in humans and in mouse, rat, zebrafish and fruit fly models. Multiple lines of evidence (such as genome-wide association data and known druggable genes) were integrated to prioritize top candidate genes for epileptogenesis and a functional enrichment analysis was carried out. Several top candidate genes, which are supported by multiple lines of genomic evidence, such as GRIN1, KCNAB1 and STX1B, were identified. Druggable genes of potential interest (such as GABRA2, GRIK1, KCNAB1 and STX4) were also identified. An enrichment of genes regulated by the MEF2 and SOX5 transcription factors and the miR-106b-5p and miR-101-3p miRNAs was found. The current work is the first meta-analysis and convergent analysis of GWES for epileptogenesis in humans and in multiple animal models, integrating results from several genomic studies. Novel candidate genes and pathways for epileptogenesis were identified in this analysis.

A systems level analysis of epileptogenesis-associated proteome alterations

Despite intense research efforts, the knowledge about the mechanisms of epileptogenesis and epilepsy is still considered incomplete and limited. However, an in-depth understanding of molecular pathophysiological processes is crucial for the rational selection of innovative biomarkers and target candidates. Here, we subjected proteomic data from different phases of a chronic rat epileptogenesis model to a comprehensive systems level analysis. Weighted Gene Co-expression Network analysis identified several modules of interconnected protein groups reflecting distinct molecular aspects of epileptogenesis in the hippocampus and the parahippocampal cortex. Characterization of these modules did not only further validate the data but also revealed regulation of molecular processes not described previously in the context of epilepsy development. The data sets also provide valuable information about temporal patterns, which should be taken into account for development of preventive strategies in particular when it comes to multi-targeting network pharmacology approaches. In addition, principal component analysis suggests candidate biomarkers, which might inform the design of novel molecular imaging approaches aiming to predict epileptogenesis during different phases or confirm epilepsy manifestation. Further studies are necessary to distinguish between molecular alterations, which correlate with epileptogenesis versus those reflecting a mere consequence of the status epilepticus.

Identification of novel gene and pathway targets for human epilepsy treatment

BACKGROUND

The aim of this study was to explore epilepsy-related mechanism so as to figure out the possible targets for epilepsy treatment.

METHODS

The gene expression profile dataset GES32534 was downloaded from Gene Expression Omnibus database. We identified the differentially expressed genes (DEGs) by Affy package. Then the DEGs were used to perform gene ontology (GO) and pathway enrichment analyses. Furthermore, a protein-protein interaction (PPI) network was constructed with the DEGs followed by co-expression modules construction and analysis.

RESULTS

Total 420 DEGs were screened out, including 214 up-regulated and 206 down-regulated genes. Functional enrichment analysis revealed that down-regulated genes were mainly involved in the process of immunity regulation and biological repairing process while up-regulated genes were closely related to transporter activity. PPI network analysis showed the top ten genes with high degrees were all down-regulated, among which FN1 had the highest degree. The up-regulated and down-regulated DEGs in the PPI network generated two obvious sub-co-expression modules, respectively. In up-co-expression module, SCN3B (sodium channel, voltage gated, type III beta subunit) was enriched in GO:0006814 ~ sodium ion transport. In down-co-expression module, C1QB (complement C1s), C1S (complement component 1, S subcomponent) and CFI (complement factor I) were enriched in GO:0006955 ~ immune response.

CONCLUSION

The immune response and complement system play a major role in the pathogenesis of epilepsy. Additionally, C1QB, C1S, CFI, SCN3B and FN1 may be potential therapeutic targets for epilepsy.