A systems level, functional genomics analysis of chronic epilepsy

High-content screening identifies a small molecule that restores AP-4-dependent protein trafficking in neuronal models of AP-4-associated hereditary spastic paraplegia

Unbiased phenotypic screens in patient-relevant disease models offer the potential to detect therapeutic targets for rare diseases. In this study, we developed a high-throughput screening assay to identify molecules that correct aberrant protein trafficking in adapter protein complex 4 (AP-4) deficiency, a rare but prototypical form of childhood-onset hereditary spastic paraplegia characterized by mislocalization of the autophagy protein ATG9A. Using high-content microscopy and an automated image analysis pipeline, we screened a diversity library of 28,864 small molecules and identified a lead compound, BCH-HSP-C01, that restored ATG9A pathology in multiple disease models, including patient-derived fibroblasts and induced pluripotent stem cell-derived neurons. We used multiparametric orthogonal strategies and integrated transcriptomic and proteomic approaches to delineate potential mechanisms of action of BCH-HSP-C01. Our results define molecular regulators of intracellular ATG9A trafficking and characterize a lead compound for the treatment of AP-4 deficiency, providing important proof-of-concept data for future studies.

High-Content Small Molecule Screen Identifies a Novel Compound That Restores AP-4-Dependent Protein Trafficking in Neuronal Models of AP-4-Associated Hereditary Spastic Paraplegia

Unbiased phenotypic screens in patient-relevant disease models offer the potential to detect novel therapeutic targets for rare diseases. In this study, we developed a high-throughput screening assay to identify molecules that correct aberrant protein trafficking in adaptor protein complex 4 (AP-4) deficiency, a rare but prototypical form of childhood-onset hereditary spastic paraplegia, characterized by mislocalization of the autophagy protein ATG9A. Using high-content microscopy and an automated image analysis pipeline, we screened a diversity library of 28,864 small molecules and identified a lead compound, C-01, that restored ATG9A pathology in multiple disease models, including patient-derived fibroblasts and induced pluripotent stem cell-derived neurons. We used multiparametric orthogonal strategies and integrated transcriptomic and proteomic approaches to delineate putative molecular targets of C-01 and potential mechanisms of action. Our results define molecular regulators of intracellular ATG9A trafficking and characterize a lead compound for the treatment of AP-4 deficiency, providing important proof-of-concept data for future Investigational New Drug (IND)-enabling studies.

Multi-omics Integration and Epilepsy: Towards a Better Understanding of Biological Mechanisms

The epilepsies are a group of complex neurological disorders characterised by recurrent seizures. Approximately 30% of patients fail to respond to anti-seizure medications, despite the recent introduction of many new drugs. The molecular processes underlying epilepsy development are not well understood and this knowledge gap impedes efforts to identify effective targets and develop novel therapies against epilepsy. Omics studies allow a comprehensive characterisation of a class of molecules. Omics-based biomarkers have led to clinically validated diagnostic and prognostic tests for personalised oncology, and more recently for non-cancer diseases. We believe that, in epilepsy, the full potential of multi-omics research is yet to be realised and we envisage that this review will serve as a guide to researchers planning to undertake omics-based mechanistic studies.

Nrf2 is expressed more extensively in neurons than in astrocytes following an acute epileptic seizure in rats

The modulation of the nuclear factor erythroid 2-like 2 (Nrf2) activity has been reported to be implicated in the pathology of various neurological disorders, including epilepsy. Previous studies have demonstrated that Nrf2 is activated in the post-status epilepticus rat model; however, the spatiotemporal as well as cell type-specific expression of Nrf2 following brief epileptic seizures remains unclear. Here, we evaluated how an acute epileptic seizure affected the expression of Nrf2 and its downstream genes in the rats' cortex and the hippocampus up to 1 week following the induced seizure. We found that after a pentylenetetrazol-induced seizure, Nrf2 significantly increased at 24 h at the mRNA level and 3 h at the protein level in the cortex. In the hippocampus, the Nrf2 mRNA level peaked at 3 h after the seizure, and no significant changes were observed in the protein level. Interestingly, the mRNA level of Nrf2 downstream genes peaked at 3-6 h after seizure in both the cortex and the hippocampus. A significant increase in the expression of Nrf2 was observed in the neuronal population of CA1 and CA3 regions of the hippocampus, as well as in the cortex. Moreover, we observed no change in the co-localization of Nrf2 with astrocytes neither in the cortex nor in CA1 and CA3. Our results revealed that following a brief acute epileptic seizure, the expression of Nrf2 and its downstream genes is transiently increased and peaked at early timepoints after the seizure predominantly in the hippocampus, and this expression is restricted to the neuronal population.

De novo mutations within metabolism networks of amino acid/protein/energy in Chinese autistic children with intellectual disability

BACKGROUND

Autism spectrum disorder (ASD) is often accompanied by intellectual disability (ID). Despite extensive studies, however, the genetic basis for this comorbidity is still not clear. In this study, we tried to develop an analyzing pipeline for de novo mutations and possible pathways related to ID phenotype in ASD. Whole-exome sequencing (WES) was performed to screen de novo mutations and candidate genes in 79 ASD children together with their parents (trios). The de novo altering genes and relative pathways which were associated with ID phenotype were analyzed. The connection nodes (genes) of above pathways were selected, and the diagnostic value of these selected genes for ID phenotype in the study population was also evaluated.

RESULTS

We identified 89 de novo mutant genes, of which 34 genes were previously reported to be associated with ASD, including double hits in the EGF repeats of NOTCH1 gene (p.V999M and p.S1027L). Interestingly, of these 34 genes, 22 may directly affect intelligence quotient (IQ). Further analyses revealed that these IQ-related genes were enriched in protein synthesis, energy metabolism, and amino acid metabolism, and at least 9 genes (CACNA1A, ALG9, PALM2, MGAT4A, PCK2, PLEKHA1, PSME3, ADI1, and TLE3) were involved in all these three pathways. Seven patients who harbored these gene mutations showed a high prevalence of a low IQ score (< 70), a non-verbal language, and an early diagnostic age (< 4 years). Furthermore, our panel of these 9 genes reached a 10.2% diagnostic rate (5/49) in early diagnostic patients with a low IQ score and also reached a 10% diagnostic yield in those with both a low IQ score and non-verbal language (4/40).

CONCLUSION

We found some new genetic disposition for ASD accompanied with intellectual disability in this study. Our results may be helpful for etiologic research and early diagnoses of intellectual disability in ASD. Larger population studies and further mechanism studies are warranted.

Connectivity Mapping Using a Novel sv2a Loss-of-Function Zebrafish Epilepsy Model as a Powerful Strategy for Anti-epileptic Drug Discovery

Synaptic vesicle glycoprotein 2A (SV2A) regulates action potential-dependent neurotransmitter release and is commonly known as the primary binding site of an approved anti-epileptic drug, levetiracetam. Although several rodent knockout models have demonstrated the importance of SV2A for functional neurotransmission, its precise physiological function and role in epilepsy pathophysiology remains to be elucidated. Here, we present a novel sv2a knockout model in zebrafish, a vertebrate with complementary advantages to rodents. We demonstrated that 6 days post fertilization homozygous sv2a^–/– mutant zebrafish larvae, but not sv2a^+/– and sv2a^+/+ larvae, displayed locomotor hyperactivity and spontaneous epileptiform discharges, however, no major brain malformations could be observed. A partial rescue of this epileptiform brain activity could be observed after treatment with two commonly used anti-epileptic drugs, valproic acid and, surprisingly, levetiracetam. This observation indicated that additional targets, besides Sv2a, maybe are involved in the protective effects of levetiracetam against epileptic seizures. Furthermore, a transcriptome analysis provided insights into the neuropathological processes underlying the observed epileptic phenotype. While gene expression profiling revealed only one differentially expressed gene (DEG) between wildtype and sv2a^+/– larvae, there were 4386 and 3535 DEGs between wildtype and sv2a^–/–, and sv2a^+/– and sv2a^–/– larvae, respectively. Pathway and gene ontology (GO) enrichment analysis between wildtype and sv2a^–/– larvae revealed several pathways and GO terms enriched amongst up- and down-regulated genes, including MAPK signaling, synaptic vesicle cycle, and extracellular matrix organization, all known to be involved in epileptogenesis and epilepsy. Importantly, we used the Connectivity map database to identify compounds with opposing gene signatures compared to the one observed in sv2a^–/– larvae, to finally rescue the epileptic phenotype. Two out of three selected compounds rescued electrographic discharges in sv2a^–/– larvae, while negative controls did not. Taken together, our results demonstrate that sv2a deficiency leads to increased seizure vulnerability and provide valuable insight into the functional importance of sv2a in the brain in general. Furthermore, we provided evidence that the concept of connectivity mapping represents an attractive and powerful approach in the discovery of novel compounds against epilepsy.

Antiseizure medication discovery: Recent and future paradigm shifts

Despite the ever-increasing number of available options for the treatment of epilepsies and the remarkable advances on the understanding of their pathophysiology, the proportion of refractory patients has remained approximately unmodified during the last 100 years. How efficient are we translating positive outcomes from basic research to clinical trials and/or the clinical scenario? It is possible that fresh thinking and exploration of new paradigms is required to arrive at truly novel therapeutic solutions, as seemingly proven by recently approved first-in-class antiseizure medications and drug candidates undergoing late clinical trials. Here, the author discusses some approximations in line with the network pharmacology philosophy, which may result in highly innovative (and, hopefully, safer and/or more efficacious) medications for the control of seizures, as embodied with some recent examples in the field, namely tailored multi-target agents and low-affinity ligands.

Hippocampal CA3 transcriptional modules associated with granule cell alterations and cognitive impairment in refractory mesial temporal lobe epilepsy patients

In about a third of the patients with epilepsy the seizures are not drug-controlled. The current limitation of the antiepileptic drug therapy derives from an insufficient understanding of epilepsy pathophysiology. In order to overcome this situation, it is necessary to consider epilepsy as a disturbed network of interactions, instead of just looking for changes in single molecular components. Here, we studied CA3 transcriptional signatures and dentate gyrus histopathologic alterations in hippocampal explants surgically obtained from 57 RMTLE patients submitted to corticoamygdalohippocampectomy. By adopting a systems biology approach, integrating clinical, histopathological, and transcriptomic data (weighted gene co-expression network analysis), we were able to identify transcriptional modules highly correlated with age of disease onset, cognitive dysfunctions, and granule cell alterations. The enrichment analysis of transcriptional modules and the functional characterization of the highly connected genes in each trait-correlated module allowed us to unveil the modules’ main biological functions, paving the way for further investigations on their roles in RMTLE pathophysiology. Moreover, we found 15 genes with high gene significance values which have the potential to become novel biomarkers and/or therapeutic targets in RMTLE.

Gene Expression Differences Between Young Adults Based on Trauma History and Post-traumatic Stress Disorder

Background: The purpose of this study was to identify gene expression differences associated with post-traumatic stress disorder (PTSD) and trauma exposure (TE) in a three-group study design comprised of those with and without trauma exposure and PTSD. Methods: We conducted gene expression and gene network analyses in a sample (n = 45) composed of female subjects of European Ancestry (EA) with PTSD, TE without PTSD, and controls. Results: We identified 283 genes differentially expressed between PTSD-TE groups. In an independent sample of Veterans (n = 78) a small minority of these genes were also differentially expressed. We identified 7 gene network modules significantly associated with PTSD and TE (Bonferroni corrected p ≤ 0.05), which at a false discovery rate (FDR) of q ≤ 0.2, were significantly enriched for biological pathways involved in focal adhesion, neuroactive ligand receptor interaction, and immune related processes among others. Conclusions: This study uses gene network analyses to identify significant gene modules associated with PTSD, TE, and controls. On an individual gene level, we identified a large number of differentially expressed genes between PTSD-TE groups, a minority of which were also differentially expressed in the independent sample. We also demonstrate a lack of network module preservation between PTSD and TE, suggesting that the molecular signature of PTSD and trauma are likely independent of each other. Our results provide a basis for the identification of likely disease pathways and biomarkers involved in the etiology of PTSD.

miR-135a-5p inhibitor protects glial cells against apoptosis via targeting SIRT1 in epilepsy

Epilepsy is a common neurological disease that can induce severe physiological brain damage, including nerve cell apoptosis. MicroRNAs (miRs) have been widely investigated in epilepsy therapy. miR-135a-5p expression levels in children with temporal lobe epilepsy were found to be significantly increased. However, whether miR-135a-5p participates in epilepsy-induced cell apoptosis is not completely understood. In the present study, an in vitro model of epilepsy in BV2 microglia cells was induced using 6-µm kainic acid (KA). Reverse-transcription quantitative PCR was performed to analyze miR-135a-5p and sirtuin 1 (SIRT1) mRNA expression levels. Western blotting was performed to measure SIRT1 protein expression levels. BV2 cell proliferation and apoptosis were assessed by performing MTT assays and flow cytometry, respectively. A BCA protein assay kit was used to detect caspase-3 and caspase-9 activities. TargetScan and dual luciferase reporter assays were performed to investigate the interaction between miR-135a-5p and the 3'-untranslated region (UTR) of SIRT1. miR-135a-5p expression was significantly increased in the KA-induced in vitro model of epilepsy in BV2 microglia. miR-135a-5p inhibitor effectively promoted BV2 microglia proliferation and inhibited microglia apoptosis, whereas small interfering RNA targeting SIRT1 significantly repressed BV2 microglia proliferation and induced microglia apoptosis. In addition, the results demonstrated that the 3'-UTR of SIRT1 mRNA was targeted by miR-135a-5p, and SIRT1 knockdown attenuated miR-135a-5p inhibitor-mediated effects on epilepsy. In summary, the results of the present study identified the role of miR-135a-5p inhibitor pretreatment in protecting nerve cells against epilepsy-induced apoptosis and provided a novel strategy for the treatment of neural damage in seizures.

Differential Glial Activation in Early Epileptogenesis-Insights From Cell-Specific Analysis of DNA Methylation and Gene Expression in the Contralateral Hippocampus

Background and Aims: Morphological changes in mesial temporal lobe epilepsy with hippocampal sclerosis (mTLE-HS) are well-characterized. Yet, it remains elusive whether these are a consequence of seizures or originate from a hitherto unknown underlying pathology. We recently published data on changes in gene expression and DNA methylation in the ipsilateral hippocampus (ILH) using the intracortical kainate mouse model of mTLE-HS. In order to explore the effects of epileptic activity alone and also to further disentangle what triggers morphological alterations, we investigated glial and neuronal changes in gene expression and DNA methylation in the contralateral hippocampus (CLH). Methods: The intracortical kainic acid mouse model of mTLE-HS was used to elicit status epilepticus. Hippocampi contralateral to the injection site from eight kainate-injected and eight sham mice were extracted and shock frozen at 24 h post-injection. Glial and neuronal nuclei were sorted by flow cytometry. Alterations in gene expression and DNA methylation were assessed using reduced representation bisulfite sequencing and RNA sequencing. The R package edgeR was used for statistical analysis. Results: The CLH featured substantial, mostly cell-specific changes in both gene expression and DNA methylation in glia and neurons. While changes in gene expression overlapped to a great degree between CLH and ILH, alterations in DNA methylation did not. In the CLH, we found a significantly lower number of glial genes up- and downregulated compared to previous results from the ILH. Furthermore, several genes and pathways potentially involved in anti-epileptogenic effects were upregulated in the CLH. By comparing gene expression data from the CLH to previous results from the ILH (featuring hippocampal sclerosis), we derive potential upstream targets for epileptogenesis, including glial Cox2 and Cxcl10. Conclusion: Despite the absence of morphological changes, the CLH displays substantial changes in gene expression and DNA methylation. We find that gene expression changes related to potential anti-epileptogenic effects seem to dominate compared to the pro-epileptogenic effects in the CLH and speculate whether this imbalance contributes to prevent morphological alterations like neuronal death and reactive gliosis.

In vivo CRISPRa decreases seizures and rescues cognitive deficits in a rodent model of epilepsy

Epilepsy is a major health burden, calling for new mechanistic insights and therapies. CRISPR-mediated gene editing shows promise to cure genetic pathologies, although hitherto it has mostly been applied ex vivo. Its translational potential for treating non-genetic pathologies is still unexplored. Furthermore, neurological diseases represent an important challenge for the application of CRISPR, because of the need in many cases to manipulate gene function of neurons in situ. A variant of CRISPR, CRISPRa, offers the possibility to modulate the expression of endogenous genes by directly targeting their promoters. We asked if this strategy can effectively treat acquired focal epilepsy, focusing on ion channels because their manipulation is known be effective in changing network hyperactivity and hypersynchronziation. We applied a doxycycline-inducible CRISPRa technology to increase the expression of the potassium channel gene Kcna1 (encoding K_v1.1) in mouse hippocampal excitatory neurons. CRISPRa-mediated K_v1.1 upregulation led to a substantial decrease in neuronal excitability. Continuous video-EEG telemetry showed that AAV9-mediated delivery of CRISPRa, upon doxycycline administration, decreased spontaneous generalized tonic-clonic seizures in a model of temporal lobe epilepsy, and rescued cognitive impairment and transcriptomic alterations associated with chronic epilepsy. The focal treatment minimizes concerns about off-target effects in other organs and brain areas. This study provides the proof-of-principle for a translational CRISPR-based approach to treat neurological diseases characterized by abnormal circuit excitability.

Inferior Colliculus Transcriptome After Status Epilepticus in the Genetically Audiogenic Seizure-Prone Hamster GASH/Sal

The Genetic Audiogenic Seizure Hamster from Salamanca (GASH/Sal), an animal model of reflex epilepsy, exhibits generalized tonic–clonic seizures in response to loud sound with the epileptogenic focus localized in the inferior colliculus (IC). Ictal events in seizure-prone strains cause gene deregulation in the epileptogenic focus, which can provide insights into the epileptogenic mechanisms. Thus, the present study aimed to determine the expression profile of key genes in the IC of the GASH/Sal after the status epilepticus. For such purpose, we used RNA-Seq to perform a comparative study between the IC transcriptome of GASH/Sal and that of control hamsters both subjected to loud sound stimulation. After filtering for normalization and gene selection, a total of 36 genes were declared differentially expressed from the RNA-seq analysis in the IC. A set of differentially expressed genes were validated by RT-qPCR showing significant differentially expression between GASH/Sal hamsters and Syrian control hamsters. The confirmed differentially expressed genes were classified on ontological categories associated with epileptogenic events similar to those produced by generalized tonic seizures in humans. Subsequently, based on the result of metabolomics, we found the interleukin-4 and 13-signaling, and nucleoside transport as presumably altered routes in the GASH/Sal model. This research suggests that seizures in GASH/Sal hamsters are generated by multiple molecular substrates, which activate biological processes, molecular processes, cellular components and metabolic pathways associated with epileptogenic events similar to those produced by tonic seizures in humans. Therefore, our study supports the use of the GASH/Sal as a valuable animal model for epilepsy research, toward establishing correlations with human epilepsy and searching new biomarkers of epileptogenesis.

Aberrant expression of PAR bZIP transcription factors is associated with epileptogenesis, focus on hepatic leukemia factor

Epilepsy is a widespread neurological disease characterized by abnormal neuronal activity resulting in recurrent seizures. There is mounting evidence that a circadian system disruption, involving clock genes and their downstream transcriptional regulators, is associated with epilepsy. In this study, we characterized the hippocampal expression of clock genes and PAR bZIP transcription factors (TFs) in a mouse model of temporal lobe epilepsy induced by intrahippocampal injection of kainic acid (KA). The expression of PAR bZIP TFs was significantly altered following KA injection as well as in other rodent models of acquired epilepsy. Although the PAR bZIP TFs are regulated by proinflammatory cytokines in peripheral tissues, we discovered that the regulation of their expression is inflammation-independent in hippocampal tissue and rather mediated by clock genes and hyperexcitability. Furthermore, we report that hepatic leukemia factor (Hlf), a member of PAR bZIP TFs family, is invariably downregulated in animal models of acquired epilepsy, regulates neuronal activity in vitro and its overexpression in dentate gyrus neurons in vivo leads to altered expression of genes associated with seizures and epilepsy. Overall, our study provides further evidence of PAR bZIP TFs involvement in epileptogenesis and points to Hlf as the key player.

Mitochondrial Regulation of the Hippocampal Firing Rate Set Point and Seizure Susceptibility

Maintaining average activity within a set-point range constitutes a fundamental property of central neural circuits. However, whether and how activity set points are regulated remains unknown. Integrating genome-scale metabolic modeling and experimental study of neuronal homeostasis, we identified mitochondrial dihydroorotate dehydrogenase (DHODH) as a regulator of activity set points in hippocampal networks. The DHODH inhibitor teriflunomide stably suppressed mean firing rates via synaptic and intrinsic excitability mechanisms by modulating mitochondrial Ca²⁺ buffering and spare respiratory capacity. Bi-directional activity perturbations under DHODH blockade triggered firing rate compensation, while stabilizing firing to the lower level, indicating a change in the firing rate set point. In vivo, teriflunomide decreased CA3-CA1 synaptic transmission and CA1 mean firing rate and attenuated susceptibility to seizures, even in the intractable Dravet syndrome epilepsy model. Our results uncover mitochondria as a key regulator of activity set points, demonstrate the differential regulation of set points and compensatory mechanisms, and propose a new strategy to treat epilepsy.

Big data sharing and analysis to advance research in post-traumatic epilepsy

We describe the infrastructure and functionality for a centralized preclinical and clinical data repository and analytic platform to support importing heterogeneous multi-modal data, automatically and manually linking data across modalities and sites, and searching content. We have developed and applied innovative image and electrophysiology processing methods to identify candidate biomarkers from MRI, EEG, and multi-modal data. Based on heterogeneous biomarkers, we present novel analytic tools designed to study epileptogenesis in animal model and human with the goal of tracking the probability of developing epilepsy over time.

Comparison of gene co-networks analysis provide a systems view of rice (Oryza sativa L.) response to Tilletia horrida infection

The biotrophic soil-borne fungus Tilletia horrida causes rice kernel smut, an important disease affecting the production of rice male sterile lines in most hybrid rice growing regions of the world. There are no successful ways of controlling this disease and there has been little study of mechanisms of resistance to T. horrida. Based on transcriptional data of different infection time points, we found 23, 782 and 23, 718 differentially expressed genes (fragments per kilobase of transcript sequence per million, FPKM >1) in Jiangcheng 3A (resistant to T. horrida) and 9311A (susceptible to T. horrida), respectively. In order to illuminate the differential responses of the two rice male sterile lines to T. horrida, we identified gene co-expression modules using the method of weighted gene co-expression network analysis (WGCNA) and compared the different biological functions of gene co-expression networks in key modules at different infection time points. The results indicated that gene co-expression networks in the two rice genotypes were different and that genes contained in some modules of the two groups may play important roles in resistance to T. horrida, such as DTH8 and OsHop/Sti1a. Furthermore, these results provide a global view of the responses of two different phenotypes to T. horrida, and assist our understanding of the regulation of expression changes after T. horrida infection.

The Antiepileptic Drug Levetiracetam Protects Against Quinolinic Acid-Induced Toxicity in the Rat Striatum

Levetiracetam (LVT) is a relatively novel antiepileptic drug (AED) known to act through binding with the synaptic vesicular 2A (SV2A) protein, thus modulating the presynaptic neurotransmitter release. The tryptophan metabolite quinolinic acid (QUIN) acts as an excitotoxin when its brain concentrations reach toxic levels under pathological conditions. Since increased neuronal excitability induced by QUIN recruits degenerative events in the brain, and novel AED is also expected to exert neuroprotective effects in their pharmacological profiles, in this work the effect of LVT (54 mg/kg, i.p., administered for seven consecutive days) was tested as a pretreatment against the toxicity evoked by the bilateral intrastriatal injection of QUIN (60 nmol/μl) to adult rats. QUIN increased the striatal levels of peroxidized lipids and carbonylated proteins as indexes of oxidative damage 24 h after its infusion. In addition, in synaptosomal fractions isolated from QUIN-lesioned rats 24 h after the toxin infusion, γ-aminobutyric acid (GABA) release was decreased, whereas glutamate (Glu) release was increased. QUIN also decreased motor activity and augmented the rate of cell damage at 7 days post-lesion. All these alterations were significantly prevented by pretreatment of rats with LVT. The results of this study show a neuroprotective role and antioxidant action of LVT against the brain damage induced by excitotoxic events.

Transcriptional profile of hippocampal dentate granule cells in four rat epilepsy models

Global expression profiling of neurologic or psychiatric disorders has been confounded by variability among laboratories, animal models, tissues sampled, and experimental platforms, with the result being that few genes demonstrate consistent expression changes. We attempted to minimize these confounds by pooling dentate granule cell transcriptional profiles from 164 rats in seven laboratories, using three status epilepticus (SE) epilepsy models (pilocarpine, kainate, self-sustained SE), plus amygdala kindling. In each epilepsy model, RNA was harvested from laser-captured dentate granule cells from six rats at four time points early in the process of developing epilepsy, and data were collected from two independent laboratories in each rodent model except SSSE. Hierarchical clustering of differentially-expressed transcripts in the three SE models revealed complete separation between controls and SE rats isolated 1 day after SE. However, concordance of gene expression changes in the SE models was only 26–38% between laboratories, and 4.5% among models, validating the consortium approach. Transcripts with unusually highly variable control expression across laboratories provide a ‘red herring’ list for low-powered studies.

Blood DNA methylation pattern is altered in mesial temporal lobe epilepsy

Mesial temporal lobe epilepsy (MTLE) is a common epileptic disorder; little is known whether it is associated with peripheral epigenetic changes. Here we compared blood whole genomic DNA methylation pattern in MTLE patients (n = 30) relative to controls (n = 30) with the Human Methylation 450 K BeadChip assay, and explored genes and pathways that were differentially methylated using bioinformatics profiling. The MTLE and control groups showed significantly different (P < 1.03e-07) DNA methylation at 216 sites, with 164 sites involved hyper- and 52 sites hypo- methylation. Two hyper- and 32 hypo-methylated sites were associated with promoters, while 87 hyper- and 43 hypo-methylated sites corresponded to coding regions. The differentially methylated genes were largely related to pathways predicted to participate in anion binding, oxidoreductant activity, growth regulation, skeletal development and drug metabolism, with the most distinct ones included SLC34A2, CLCN6, CLCA4, CYP3A43, CYP3A4 and CYP2C9. Among the MTLE patients, panels of genes also appeared to be differentially methylated relative to disease duration, resistance to anti-epileptics and MRI alterations of hippocampal sclerosis. The peripheral epigenetic changes observed in MTLE could be involved in certain disease-related modulations and warrant further translational investigations.

mTOR and MAPK: from localized translation control to epilepsy

Background

Epilepsy is one of the most common neurological diseases characterized by excessive hyperexcitability of neurons. Molecular mechanisms of epilepsy are diverse and not really understood. All in common is the misregulation of proteins that determine excitability such as potassium and sodium channels as well as GABA receptors; which are all known as biomarkers for epilepsy. Two recently identified key pathways involve the kinases mechanistic target of rapamycin (mTOR) and mitogen-activated protein kinases (MAPK). Interestingly, mRNAs coding for those biomarkers are found to be localized at or near synapses indicating a local misregulation of synthesis and activity.

Results

Research in the last decade indicates that RNA-binding proteins (RBPs) responsible for mRNA localization, stability and translation mediate local expression control. Among others, they are affected by mTOR and MAPK to guide expression of epileptic factors. These results suggest that mTOR/MAPK act on RBPs to regulate the fate of mRNAs, indicating a misregulation of protein expression at synapses in epilepsy.

Conclusion

We propose that mTOR and MAPK regulate RBPs, thereby guiding the local expression of their target-mRNAs encoding for markers of epilepsy. Thus, misregulated mTOR/MAPK-RBP interplay may result in excessive local synthesis of ion channels and receptors thereby leading to hyperexcitability. Continuous stimulation of synapses further activates mTOR/MAPK pathway reinforcing their effect on RBP-mediated expression control establishing the basis for epilepsy. Here, we highlight findings showing the tight interplay between mTOR as well as MAPK with RBPs to control expression for epileptic biomarkers.

miRNA-187-3p-Mediated Regulation of the KCNK10/TREK-2 Potassium Channel in a Rat Epilepsy Model

Regulatory RNAs play a key role in the regulation of protein expression patterns in neurological diseases. Here we studied the regulation of miRNAs in a chronic rat model of temporal lobe epilepsy. The analysis was focused on a putative link with pharmacoresponsiveness as well as the functional implications of the regulation of a selected miRNA. The findings did not reveal a difference in hippocampal miRNA expression between phenobarbital responders and nonresponders. However, when comparing rats following status epilepticus with control rats we identified 13 differentially expressed miRNAs with miRNA-187-3p being most strongly regulated. mRNAs encoding KCNK10/TREK-2 as well as DYRK2 were confirmed as targets of miRNA-187-3p. Expression of the potassium channel protein KCNK10/TREK-2 negatively correlated with hippocampal miRNA-187-3p expression and proved to be upregulated in the chronic phase of the epilepsy model. In conclusion, our data do not suggest a relevant impact of miRNA expression patterns on pharmacoresponsiveness. However, we confirmed regulation of miRNA-187-3p and demonstrated that it impacts the expression of the two-pore domain potassium channel protein KCNK10/TREK-2. Considering evidence from brain ischemia models, KCNK10/TREK-2 upregulation might serve a protective function with a beneficial impact on astrocytic potassium and glutamate homeostasis.

Synaptic Vesicle Glycoprotein 2A Ligands in the Treatment of Epilepsy and Beyond

The synaptic vesicle glycoprotein SV2A belongs to the major facilitator superfamily (MFS) of transporters and is an integral constituent of synaptic vesicle membranes. SV2A has been demonstrated to be involved in vesicle trafficking and exocytosis, processes crucial for neurotransmission. The anti-seizure drug levetiracetam was the first ligand to target SV2A and displays a broad spectrum of anti-seizure activity in various preclinical models. Several lines of preclinical and clinical evidence, including genetics and protein expression changes, support an important role of SV2A in epilepsy pathophysiology. While the functional consequences of SV2A ligand binding are not fully elucidated, studies suggest that subsequent SV2A conformational changes may contribute to seizure protection. Conversely, the recently discovered negative SV2A modulators, such as UCB0255, counteract the anti-seizure effect of levetiracetam and display procognitive properties in preclinical models. More broadly, dysfunction of SV2A may also be involved in Alzheimer's disease and other types of cognitive impairment, suggesting potential novel therapies for levetiracetam and its congeners. Furthermore, emerging data indicate that there may be important roles for two other SV2 isoforms (SV2B and SV2C) in the pathogenesis of epilepsy, as well as other neurodegenerative diseases. Utilization of recently developed SV2A positron emission tomography ligands will strengthen and reinforce the pharmacological evidence that SV2A is a druggable target, and will provide a better understanding of its role in epilepsy and other neurological diseases, aiding in further defining the full therapeutic potential of SV2A modulation.

Integrative analysis of gene expression associated with epilepsy in human epilepsy and animal models

Epilepsy is a severe neuropsychiatric disorder, the cause of which remains to be elucidated. Genome‑wide association studies, DNA microarrays and proteomes have been widely applied to identify the candidate genes involved in epileptogenesis, and integrative analyses are often capable of extracting more detailed biological information from the data. In the present study, a total number of 1,065 genes in different animal models were collected to construct an epilepsy candidate gene database. Further analyses suggested that the response to organic substances, the intracellular signaling cascade and neurological system processes were significantly enriched biological processes, and the mitogen-activated protein kinase pathway was identified as a putative epileptogenic signaling pathway. In addition, the five key genes, growth factor receptor bound 2, amyloid β (A4) precursor protein, transforming growth factor‑β, vascular endothelial growth factor and cyclin‑dependent kinase inhibitor 1, were identified as being critical as central nodes in the protein networks. Reverse transcription‑quantitative polymerase chain reaction analysis revealed that these genes were all upregulated at the mRNA level in the epileptic loci compared with the resection margin of tissue samples from the same patients diagnosed with epilepsy. The data mining performed in the present study thus was shown to be a useful tool, which may contribute to obtaining further information on epileptic disorders and delineating the molecular mechanism of the associated genes.

Neurological Diseases from a Systems Medicine Point of View

The difficulty to understand, diagnose, and treat neurological disorders stems from the great complexity of the central nervous system on different levels of physiological granularity. The individual components, their interactions, and dynamics involved in brain development and function can be represented as molecular, cellular, or functional networks, where diseases are perturbations of networks. These networks can become a useful research tool in investigating neurological disorders if they are properly tailored to reflect corresponding mechanisms. Here, we review approaches to construct networks specific for neurological disorders describing disease-related pathology on different scales: the molecular, cellular, and brain level. We also briefly discuss cross-scale network analysis as a necessary integrator of these scales.

van der Vaart A, Web BT, Bacanu SA, Kalsi G, Kendler KS, Miles MF, Dick D, Riley BP, Dumur C, Vladimirov VI. Integrating mRNA and miRNA Weighted Gene Co-Expression Networks with eQTLs in the Nucleus Accumbens of Subjects with Alcohol Dependence

Alcohol consumption is known to lead to gene expression changes in the brain. After performing weighted gene co-expression network analyses (WGCNA) on genome-wide mRNA and microRNA (miRNA) expression in Nucleus Accumbens (NAc) of subjects with alcohol dependence (AD; N = 18) and of matched controls (N = 18), six mRNA and three miRNA modules significantly correlated with AD were identified (Bonferoni-adj. p≤ 0.05). Cell-type-specific transcriptome analyses revealed two of the mRNA modules to be enriched for neuronal specific marker genes and downregulated in AD, whereas the remaining four mRNA modules were enriched for astrocyte and microglial specific marker genes and upregulated in AD. Gene set enrichment analysis demonstrated that neuronal specific modules were enriched for genes involved in oxidative phosphorylation, mitochondrial dysfunction and MAPK signaling. Glial-specific modules were predominantly enriched for genes involved in processes related to immune functions, i.e. cytokine signaling (all adj. p≤ 0.05). In mRNA and miRNA modules, 461 and 25 candidate hub genes were identified, respectively. In contrast to the expected biological functions of miRNAs, correlation analyses between mRNA and miRNA hub genes revealed a higher number of positive than negative correlations (χ2 test p≤ 0.0001). Integration of hub gene expression with genome-wide genotypic data resulted in 591 mRNA cis-eQTLs and 62 miRNA cis-eQTLs. mRNA cis-eQTLs were significantly enriched for AD diagnosis and AD symptom counts (adj. p = 0.014 and p = 0.024, respectively) in AD GWAS signals in a large, independent genetic sample from the Collaborative Study on Genetics of Alcohol (COGA). In conclusion, our study identified putative gene network hubs coordinating mRNA and miRNA co-expression changes in the NAc of AD subjects, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

Epileptogenesis

Epileptogenesis is a chronic process that can be triggered by genetic or acquired factors, and that can continue long after epilepsy diagnosis. In 2015, epileptogenesis is not a treatment indication, and there are no therapies available in clinic to treat individuals at risk of epileptogenesis. However, thanks to active research, a large number of animal models have become available for search of molecular mechanisms of epileptogenesis. The first glimpses of treatment targets and biomarkers that could be developed to become useful in clinic are in sight. However, the heterogeneity of the epilepsy condition, and the dynamics of molecular changes over the course of epileptogenesis remain as challenges to overcome.

A new model to study neurodegeneration in ataxia oculomotor apraxia type 2

Ataxia oculomotor apraxia type 2 (AOA2) is a rare autosomal recessive cerebellar ataxia. Recent evidence suggests that the protein defective in this syndrome, senataxin (SETX), functions in RNA processing to protect the integrity of the genome. To date, only patient-derived lymphoblastoid cells, fibroblasts and SETX knockdown cells were available to investigate AOA2. Recent disruption of the Setx gene in mice did not lead to neurobehavioral defects or neurodegeneration, making it difficult to study the etiology of AOA2. To develop a more relevant neuronal model to study neurodegeneration in AOA2, we derived neural progenitors from a patient with AOA2 and a control by induced pluripotent stem cell (iPSC) reprogramming of fibroblasts. AOA2 iPSC and neural progenitors exhibit increased levels of oxidative damage, DNA double-strand breaks, increased DNA damage-induced cell death and R-loop accumulation. Genome-wide expression and weighted gene co-expression network analysis in these neural progenitors identified both previously reported and novel affected genes and cellular pathways associated with senataxin dysfunction and the pathophysiology of AOA2, providing further insight into the role of senataxin in regulating gene expression on a genome-wide scale. These data show that iPSCs can be generated from patients with the autosomal recessive ataxia, AOA2, differentiated into neurons, and that both cell types recapitulate the AOA2 cellular phenotype. This represents a novel and appropriate model system to investigate neurodegeneration in this syndrome.

Search for common targets of lithium and valproic acid identifies novel epigenetic effects of lithium on the rat leptin receptor gene

Epigenetics may have an important role in mood stabilizer action. Valproic acid (VPA) is a histone deacetylase inhibitor, and lithium (Li) may have downstream epigenetic actions. To identify genes commonly affected by both mood stabilizers and to assess potential epigenetic mechanisms that may be involved in their mechanism of action, we administered Li (N = 12), VPA (N = 12), and normal chow (N = 12) to Brown Norway rats for 30 days. Genomic DNA and mRNA were extracted from the hippocampus. We used the mRNA to perform gene expression analysis on Affymetrix microarray chips, and for genes commonly regulated by both Li and VPA, we validated expression levels using quantitative real-time PCR. To identify potential mechanisms underlying expression changes, genomic DNA was bisulfite treated for pyrosequencing of key CpG island 'shores' and promoter regions, and chromatin was prepared from both hippocampal tissue and a hippocampal-derived cell line to assess modifications of histones. For most genes, we found little evidence of DNA methylation changes in response to the medications. However, we detected histone H3 methylation and acetylation in the leptin receptor gene, Lepr, following treatment with both drugs. VPA-mediated effects on histones are well established, whereas the Li effects constitute a novel mechanism of transcriptional derepression for this drug. These data support several shared transcriptional targets of Li and VPA, and provide evidence suggesting leptin signaling as an epigenetic target of two mood stabilizers. Additional work could help clarify whether leptin signaling in the brain has a role in the therapeutic action of Li and VPA in bipolar disorder.

Which insights have we gained from the kindling and post-status epilepticus models?

Experimental animal epilepsy research got a big boost since the discovery that daily mild and short (seconds) tetanic stimulations in selected brain regions led to seizures with increasing duration and severity. This model that was developed by Goddard (1967) became known as the kindling model for epileptogenesis and has become a widely used model for temporal lobe epilepsy with complex partial seizures. During the late ninety-eighties the number of publications related to electrical kindling reached its maximum. However, since the kindling procedure is rather labor intensive and animals only develop spontaneous seizures (epilepsy) after hundreds of stimulations, research has shifted toward models in which the animals exhibit spontaneous seizures after a relatively short latent period. This led to post-status epilepticus (SE) models in which animals experience SE after injection of pharmacological compounds (e.g. kainate or pilocarpine) or via electrical stimulation of (limbic) brain regions. These post-SE models are the most widely used models in epilepsy research today. However, not all aspects of mesial temporal lobe epilepsy (MTLE) are reproduced and the widespread brain damage is often a caricature of the situation in the patient. Therefore, there is a need for models that can better replicate the disease. Kindling, although already a classic model, can still offer valid clues in this context. In this paper, we review different aspects of the kindling model with emphasis on experiments in the rat. Next, we review characteristic properties of the post-SE models and compare the neuropathological, electrophysiological and molecular differences between kindling and post-SE epilepsy models. Finally, we shortly discuss the advantages and disadvantages of these models.

Molecular alterations in areas generating fast ripples in an animal model of temporal lobe epilepsy

The molecular basis of epileptogenesis is poorly characterized. Studies in humans and animal models have identified an electrophysiological signature that precedes the onset of epilepsy, which has been termed fast ripples (FRs) based on its frequency. Multiple lines of evidence implicate regions generating FRs in epileptogenesis, and FRs appear to demarcate the seizure onset zone, suggesting a role in ictogenesis as well. We performed gene expression analysis comparing areas of the dentate gyrus that generate FRs to those that do not generate FRs in a well-characterized rat model of epilepsy. We identified a small cohort of genes that are differentially expressed in FR versus non-FR brain tissue and used quantitative PCR to validate some of those that modulate neuronal excitability. Gene expression network analysis demonstrated conservation of gene co-expression between non-FR and FR samples, but examination of gene connectivity revealed changes that were most pronounced in the cm-40 module, which contains several genes associated with synaptic function and the differentially expressed genes Kcna4, Kcnv1, and Npy1r that are down-regulated in FRs. We then demonstrate that the genes within the cm-40 module are regulated by seizure activity and enriched for the targets of the RNA binding protein Elavl4. Our data suggest that seizure activity induces co-expression of genes associated with synaptic transmission and that this pattern is attenuated in areas displaying FRs, implicating the failure of this mechanism in the generation of FRs.

Systems biology, complexity, and the impact on antiepileptic drug discovery

The number of available anticonvulsant drugs increased in the period spanning over more than a century, amounting to the current panoply of nearly two dozen so-called antiepileptic drugs (AEDs). However, none of them actually prevents/reduces the post-brain insult development of epilepsy in man, and in no less than a third of patients with epilepsy, the seizures are not drug-controlled. Plausibly, the enduring limitation of AEDs' efficacy derives from the insufficient understanding of epileptic pathology. This review pinpoints the unbalanced reductionism of the analytic approaches that overlook the intrinsic complexity of epilepsy and of the drug resistance in epilepsy as the core conceptual flaw hampering the discovery of truly antiepileptogenic drugs. A rising awareness of the complexity of epileptic pathology is, however, brought about by the emergence of nonreductionist systems biology (SB) that considers the networks of interactions underlying the normal organismic functions and of SB-based systems (network) pharmacology that aims to restore pathological networks. By now, the systems pharmacology approaches of AED discovery are fairly meager, but their forthcoming development is both a necessity and a realistic prospect, explored in this review.

PPAR agonists regulate brain gene expression: relationship to their effects on ethanol consumption

Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors that act as ligand-activated transcription factors. Although prescribed for dyslipidemia and type-II diabetes, PPAR agonists also possess anti-addictive characteristics. PPAR agonists decrease ethanol consumption and reduce withdrawal severity and susceptibility to stress-induced relapse in rodents. However, the cellular and molecular mechanisms facilitating these properties have yet to be investigated. We tested three PPAR agonists in a continuous access two-bottle choice (2BC) drinking paradigm and found that tesaglitazar (PPARα/γ; 1.5 mg/kg) and fenofibrate (PPARα; 150 mg/kg) decreased ethanol consumption in male C57BL/6J mice while bezafibrate (PPARα/γ/β; 75 mg/kg) did not. We hypothesized that changes in brain gene expression following fenofibrate and tesaglitazar treatment lead to reduced ethanol drinking. We studied unbiased genomic profiles in areas of the brain known to be important for ethanol dependence, the prefrontal cortex (PFC) and amygdala, and also profiled gene expression in liver. Genomic profiles from the non-effective bezafibrate treatment were used to filter out genes not associated with ethanol consumption. Because PPAR agonists are anti-inflammatory, they would be expected to target microglia and astrocytes. Surprisingly, PPAR agonists produced a strong neuronal signature in mouse brain, and fenofibrate and tesaglitazar (but not bezafibrate) targeted a subset of GABAergic interneurons in the amygdala. Weighted gene co-expression network analysis (WGCNA) revealed co-expression of treatment-significant genes. Functional annotation of these gene networks suggested that PPAR agonists might act via neuropeptide and dopaminergic signaling pathways in the amygdala. Our results reveal gene targets through which PPAR agonists can affect alcohol consumption behavior.

The potential of antiseizure drugs and agents that act on novel molecular targets as antiepileptogenic treatments

A major goal of contemporary epilepsy research is the identification of therapies to prevent the development of recurrent seizures in individuals at risk, including those with brain injuries, infections, or neoplasms; status epilepticus; cortical dysplasias; or genetic epilepsy susceptibility. In this review we consider the evidence largely from preclinical models for the antiepileptogenic activity of a diverse range of potential therapies, including some marketed antiseizure drugs, as well as agents that act by immune and inflammatory mechanisms; reduction of oxidative stress; activation of the mammalian target of rapamycin or peroxisome proliferator-activated receptors γ pathways; effects on factors related to thrombolysis, hematopoesis, and angiogenesis; inhibition of 3-hydroxy-3-methylglutaryl-coenzyme A reducatase; brain-derived neurotrophic factor signaling; and blockade of α2 adrenergic and cannabinoid receptors. Antiepileptogenesis refers to a therapy of which the beneficial action is to reduce seizure frequency or severity outlasting the treatment period. To date, clinical trials have failed to demonstrate that antiseizure drugs have such disease-modifying activity. However, studies in animal models with levetiracetam and ethosuximide are encouraging, and clinical trials with these agents are warranted. Other promising strategies are inhibition of interleukin 1β signaling by drugs such as VX-765; modulation of sphingosine 1-phosphate signaling by drugs such as fingolimod; activation of the mammalian target of rapamycin by drugs such as rapamycin; the hormone erythropoietin; and, paradoxically, drugs such as the α2 adrenergic receptor antagonist atipamezole and the CB1 cannabinoid antagonist SR141716A (rimonabant) with proexcitatory activity. These approaches could lead to a new paradigm in epilepsy drug therapy where treatment for a limited period prevents the occurrence of spontaneous seizures, thus avoiding lifelong commitment to symptomatic treatment.

Protection genes in nucleus accumbens shell affect vulnerability to nicotine self-administration across isogenic strains of adolescent rat

Classical genetic studies show the heritability of cigarette smoking is 0.4–0.6, and that multiple genes confer susceptibility and resistance to smoking. Despite recent advances in identifying genes associated with smoking behaviors, the major source of this heritability and its impact on susceptibility and resistance are largely unknown. Operant self-administration (SA) of intravenous nicotine is an established model for smoking behavior. We recently confirmed that genetic factors exert strong control over nicotine intake in isogenic rat strains. Because the processing of afferent dopaminergic signals by nucleus accumbens shell (AcbS) is critical for acquisition and maintenance of motivated behaviors reinforced by nicotine, we hypothesized that differential basal gene expression in AcbS accounts for much of the strain-to-strain variation in nicotine SA. We therefore sequenced the transcriptome of AcbS samples obtained by laser capture microdissection from 10 isogenic adolescent rat strains and compared all RNA transcript levels with behavior. Weighted gene co-expression network analysis, a systems biology method, found 12 modules (i.e., unique sets of genes that covary across all samples) that correlated (p<0.05) with amount of self-administered nicotine; 9 of 12 correlated negatively, implying a protective role. PCR confirmed selected genes from these modules. Chilibot, a literature mining tool, identified 15 genes within 1 module that were nominally associated with cigarette smoking, thereby providing strong support for the analytical approach. This is the first report demonstrating that nicotine intake by adolescent rodents is associated with the expression of specific genes in AcbS of the mesolimbic system, which controls motivated behaviors. These findings provide new insights into genetic mechanisms that predispose or protect against tobacco addiction.

Evaluation of seizure foci and genes in the Lgi1(L385R/+) mutant rat

Mutations in the leucine-rich, glioma inactivated 1 (LGI1) gene have been identified in patients with autosomal dominant lateral temporal lobe epilepsy (ADLTE). We previously reported that Lgi1 mutant rats, carrying a missense mutation (L385R) generated by gene-driven N-ethyl-N-nitrosourea (ENU) mutagenesis, showed generalized tonic-clonic seizures (GTCS) in response to acoustic stimuli. In the present study, we assessed clinically relevant features of Lgi1 heterozygous mutant rats (Lgi1(L385R/+)) as an animal model of ADLTE. First, to explore the focus of the audiogenic seizures, we performed electroencephalography (EEG) and brain Fos immunohistochemistry in Lgi1(L385R/+) and wild type rats. EEG showed unique seizure patterns (e.g., bilateral rhythmic spikes) in Lgi1(L385R/+) rats with GTCS. An elevated level of Fos expression indicated greater neural excitability to acoustic stimuli in Lgi1(L385R/+) rats, especially in the temporal lobe, thalamus and subthalamic nucleus. Finally, microarray analysis revealed a number of differentially expressed genes that may be involved in epilepsy. These results suggest that Lgi1(L385R/+) rats are useful as an animal model of human ADLTE.

Are alterations in transmitter receptor and ion channel expression responsible for epilepsies?

Neuronal voltage-gated ion channels and ligand-gated synaptic receptors play a critical role in maintaining the delicate balance between neuronal excitation and inhibition within neuronal networks in the brain. Changes in expression of voltage-gated ion channels, in particular sodium, hyperpolarization-activated cyclic nucleotide-gated (HCN) and calcium channels, and ligand-gated synaptic receptors, in particular GABA and glutamate receptors, have been reported in many types of both genetic and acquired epilepsies, in animal models and in humans. In this chapter we review these and discuss the potential pathogenic role they may play in the epilepsies.

Characterization of protein species and weighted protein co-expression network regulation of Escherichia coli in response to serum killing using a 2-DE based proteomics approach

Posttranslational modifications, providing covalent alterations to extend their functions, show protein species on 2-DE gels, but our knowledge on protein species is still limited. In the present study, characteristics of protein species are determined in Escherichia coli using 2-DE based proteomics. In the E. coli proteome, 691 unique proteins (representing 1096 protein spots) accounting for 15.37% of gene-coding proteins of the bacterium are identified. Out of them, 191 have 596 protein species. Proteins with higher abundance, a higher proportion of Glu, Gly, Lys, and higher pI are more likely to have protein species. Further investigation on bacterial serum resistance indicates that more proteins with protein species are found in the bacterium in response to serum stress. A weighted protein co-expression network shows that protein species are related to topological connection as a result of protein regulation. The node protein IleS is demonstrated to contribute to serum resistance using a gene-deleted mutant. These results have revealed general characteristic features of bacterial species, and also provided novel insights into the biological significance of bacterial protein species, particularly the role in serum resistance.

Complex network analysis of CA3 transcriptome reveals pathogenic and compensatory pathways in refractory temporal lobe epilepsy

We previously described - studying transcriptional signatures of hippocampal CA3 explants - that febrile (FS) and afebrile (NFS) forms of refractory mesial temporal lobe epilepsy constitute two distinct genomic phenotypes. That network analysis was based on a limited number (hundreds) of differentially expressed genes (DE networks) among a large set of valid transcripts (close to two tens of thousands). Here we developed a methodology for complex network visualization (3D) and analysis that allows the categorization of network nodes according to distinct hierarchical levels of gene-gene connections (node degree) and of interconnection between node neighbors (concentric node degree). Hubs are highly connected nodes, VIPs have low node degree but connect only with hubs, and high-hubs have VIP status and high overall number of connections. Studying the whole set of CA3 valid transcripts we: i) obtained complete transcriptional networks (CO) for FS and NFS phenotypic groups; ii) examined how CO and DE networks are related; iii) characterized genomic and molecular mechanisms underlying FS and NFS phenotypes, identifying potential novel targets for therapeutic interventions. We found that: i) DE hubs and VIPs are evenly distributed inside the CO networks; ii) most DE hubs and VIPs are related to synaptic transmission and neuronal excitability whereas most CO hubs, VIPs and high hubs are related to neuronal differentiation, homeostasis and neuroprotection, indicating compensatory mechanisms. Complex network visualization and analysis is a useful tool for systems biology approaches to multifactorial diseases. Network centrality observed for hubs, VIPs and high hubs of CO networks, is consistent with the network disease model, where a group of nodes whose perturbation leads to a disease phenotype occupies a central position in the network. Conceivably, the chance for exerting therapeutic effects through the modulation of particular genes will be higher if these genes are highly interconnected in transcriptional networks.

RMaNI: Regulatory Module Network Inference framework

Background

Cell survival and development are orchestrated by complex interlocking programs of gene activation and repression. Understanding how this gene regulatory network (GRN) functions in normal states, and is altered in cancers subtypes, offers fundamental insight into oncogenesis and disease progression, and holds great promise for guiding clinical decisions. Inferring a GRN from empirical microarray gene expression data is a challenging task in cancer systems biology. In recent years, module-based approaches for GRN inference have been proposed to address this challenge. Despite the demonstrated success of module-based approaches in uncovering biologically meaningful regulatory interactions, their application remains limited a single condition, without supporting the comparison of multiple disease subtypes/conditions. Also, their use remains unnecessarily restricted to computational biologists, as accurate inference of modules and their regulators requires integration of diverse tools and heterogeneous data sources, which in turn requires scripting skills, data infrastructure and powerful computational facilities. New analytical frameworks are required to make module-based GRN inference approach more generally useful to the research community.

Results

We present the RMaNI (Regulatory Module Network Inference) framework, which supports cancer subtype-specific or condition specific GRN inference and differential network analysis. It combines both transcriptomic as well as genomic data sources, and integrates heterogeneous knowledge resources and a set of complementary bioinformatic methods for automated inference of modules, their condition specific regulators and facilitates downstream network analyses and data visualization. To demonstrate its utility, we applied RMaNI to a hepatocellular microarray data containing normal and three disease conditions. We demonstrate that how RMaNI can be employed to understand the genetic architecture underlying three disease conditions. RMaNI is freely available at http://inspect.braembl.org.au/bi/inspect/rmani

Conclusion

RMaNI makes available a workflow with comprehensive set of tools that would otherwise be challenging for non-expert users to install and apply. The framework presented in this paper is flexible and can be easily extended to analyse any dataset with multiple disease conditions.

Gene expression profiling of a hypoxic seizure model of epilepsy suggests a role for mTOR and Wnt signaling in epileptogenesis

Microarray profiling was used to investigate gene expression in the hypoxic seizure model of acquired epilepsy in the rat, with the aim of characterizing functional pathways which are persistently activated or repressed during epileptogenesis. Hippocampal and cortical tissues were transcriptionally profiled over a one week period following an initial series of seizures induced by mild hypoxia at post-natal day 10 (P10), and the gene expression data was then analyzed with a focus on gene set enrichment analysis, an approach which emphasizes regulation of entire pathways rather than of individual genes. Animals were subjected to one of three conditions: a control with no hypoxia, hypoxic seizures, and hypoxic seizures followed by treatment with the AMPAR antagonist NBQX, a compound currently proposed to be a modulator of epileptogenesis. While temporal gene expression in the control samples was found to be consistent with known processes of neuronal maturation in the rat for the given time window, the hypoxic seizure response was found to be enriched for components of the PI3K/mTOR and Wnt signaling pathways, alongside gene sets representative of glutamatergic, synaptic and axonal processes, perhaps regulated as a downstream consequence of activation of these pathways. Wnt signaling components were also found enriched in the more specifically epileptogenic NBQX-responsive gene set. While activation of the mTOR pathway is consistent with its known role in epileptogenesis and strengthens the case for mTOR or PI3K pathway inhibitors as potential anti-epileptogenic drugs, investigation of the role of Wnt signaling and the effect of appropriate inhibitors might offer a parallel avenue of research toward anti-epileptogenic treatment of epilepsy.

The pathology of magnetic-resonance-imaging-negative epilepsy

Patients with magnetic-resonance-imaging (MRI)-negative (or 'nonlesional') pharmacoresistant focal epilepsy are the most challenging group undergoing presurgical evaluation. Few large-scale studies have systematically reviewed the pathological substrates underlying MRI-negative epilepsies. In the current study, histopathological specimens were retrospectively reviewed from MRI-negative epilepsy patients (n=95, mean age=30 years, 50% female subjects). Focal cortical dysplasia cases were classified according to the International League Against Epilepsy (ILAE) and Palmini et al classifications. The most common pathologies found in this MRI-negative cohort included: focal cortical dysplasia (n=43, 45%), gliosis (n=21, 22%), hamartia+gliosis (n=12, 13%), and hippocampal sclerosis (n=9, 9%). The majority of focal cortical dysplasia were ILAE type I (n=37) or Palmini type I (n=39). Seven patients had no identifiable pathological abnormalities. The existence of positive pathology was not significantly associated with age or temporal/extratemporal resection. Follow-up data post surgery was available in 90 patients; 63 (70%) and 57 (63%) attained seizure freedom at 6 and 12 months, respectively. The finding of positive pathology was significantly associated with seizure-free outcome at 6 months (P=0.035), but not at 12 months. In subgroup analysis, the focal cortical dysplasia group was not significantly correlated with seizure-free outcome, as compared with the negative-pathology groups at either 6 or 12 months. Of note, the finding of hippocampal sclerosis had a significant positive correlation with seizure-free outcome when compared with the negative-pathology group (P=0.009 and 0.004 for 6- and 12-month outcome, respectively). Absence of a significant histopathology in the resected surgical specimen did not preclude seizure freedom. In conclusion, our study highlights the heterogeneity of epileptic pathologies in MRI-negative epilepsies, with focal cortical dysplasia being the most common finding. The existence of positive pathology in surgical specimen may be a good indication for short-term good seizure outcome. There is a small subset of cases in which no pathological abnormalities are identified.

Emerging complexity of the HuD/ELAVl4 gene; implications for neuronal development, function, and dysfunction

Precise control of messenger RNA (mRNA) processing and abundance are increasingly being recognized as critical for proper spatiotemporal gene expression, particularly in neurons. These regulatory events are governed by a large number of trans-acting factors found in neurons, most notably RNA-binding proteins (RBPs) and micro-RNAs (miRs), which bind to specific cis-acting elements or structures within mRNAs. Through this binding mechanism, trans-acting factors, particularly RBPs, control all aspects of mRNA metabolism, ranging from altering the transcription rate to mediating mRNA degradation. In this context the best-characterized neuronal RBP, the Hu/ELAVl family member HuD, is emerging as a key component in multiple regulatory processes--including pre-mRNA processing, mRNA stability, and translation--governing the fate of a substantial amount of neuronal mRNAs. Through its ability to regulate mRNA metabolism of diverse groups of functionally similar genes, HuD plays important roles in neuronal development and function. Furthermore, compelling evidence indicates supplementary roles for HuD in neuronal plasticity, in particular, recovery from axonal injury, learning and memory, and multiple neurological diseases. The purpose of this review is to provide a detailed overview of the current knowledge surrounding the expression and roles of HuD in the nervous system. Additionally, we outline the present understanding of the molecular mechanisms presiding over the localization, abundance, and function of HuD in neurons.

Gene coexpression networks reveal key drivers of phenotypic divergence in lake whitefish

A functional understanding of processes involved in adaptive divergence is one of the awaiting opportunities afforded by high-throughput transcriptomic technologies. Functional analysis of coexpressed genes has succeeded in the biomedical field in identifying key drivers of disease pathways. However, in ecology and evolutionary biology, functional interpretation of transcriptomic data is still limited. Here, we used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify modules of coexpressed genes in muscle and brain tissue of a lake whitefish backcross progeny. Modules were connected to gradients of known adaptive traits involved in the ecological speciation process between benthic and limnetic ecotypes. Key drivers, that is, hub genes of functional modules related to reproduction, growth, and behavior were identified, and module preservation was assessed in natural populations. Using this approach, we identified modules of coexpressed genes involved in phenotypic divergence and their key drivers, and further identified a module part specifically rewired in the backcross progeny. Functional analysis of transcriptomic data can significantly contribute to the understanding of the mechanisms underlying ecological speciation. Our findings point to bone morphogenetic protein and calcium signaling as common pathways involved in coordinated evolution of trophic behavior, trophic morphology (gill rakers), and reproduction. Results also point to pathways implicating hemoglobins and constitutive stress response (HSP70) governing growth in lake whitefish.

Systems biology of complex symptom profiles: capturing interactivity across behavior, brain and immune regulation

As our thinking about the basic principles of biology and medicine continue to evolve, the importance of context and regulatory interaction is becoming increasingly obvious. Biochemical and physiological components do not exist in isolation but instead are part of a tightly integrated network of interacting elements that ensure robustness and support the emergence of complex behavior. This integration permeates all levels of biology from gene regulation, to immune cell signaling, to coordinated patterns of neuronal activity and the resulting psychosocial interaction. Systems biology is an emerging branch of science that sits as a translational catalyst at the interface of the life and computational sciences. While there is no universally accepted definition of systems biology, we attempt to provide an overview of some the basic unifying concepts and current efforts in the field as they apply to illnesses where brain and subsequent behavior are a chief component, for example autism, schizophrenia, depression, and others. Methods in this field currently constitute a broad mosaic that stretches across multiple scales of biology and physiological compartments. While this work by no means constitutes an exhaustive list of all these methods, this work highlights the principal sub-disciplines presently driving the field as well as future directions of progress.