Genome-wide analysis of differential RNA editing in epilepsy.

The recoding of genetic information through RNA editing contributes to proteomic diversity, but the extent and significance of RNA editing in disease is poorly understood. In particular, few studies have investigated the relationship between RNA editing and disease at a genome-wide level. Here, we developed a framework for the genome-wide detection of RNA sites that are differentially edited in disease. Using RNA-sequencing data from 100 hippocampi from mice with epilepsy (pilocarpine-temporal lobe epilepsy model) and 100 healthy control hippocampi, we identified 256 RNA sites (overlapping with 87 genes) that were significantly differentially edited between epileptic cases and controls. The degree of differential RNA editing in epileptic mice correlated with frequency of seizures, and the set of genes differentially RNA-edited between case and control mice were enriched for functional terms highly relevant to epilepsy, including "neuron projection" and "seizures." Genes with differential RNA editing were preferentially enriched for genes with a genetic association to epilepsy. Indeed, we found that they are significantly enriched for genes that harbor nonsynonymous de novo mutations in patients with epileptic encephalopathy and for common susceptibility variants associated with generalized epilepsy. These analyses reveal a functional convergence between genes that are differentially RNA-edited in acquired symptomatic epilepsy and those that contribute risk for genetic epilepsy. Taken together, our results suggest a potential role for RNA editing in the epileptic hippocampus in the occurrence and severity of epileptic seizures.

Describing the genetic architecture of epilepsy through heritability analysis.

Epilepsy is a disease with substantial missing heritability; despite its high genetic component, genetic association studies have had limited success detecting common variants which influence susceptibility. In this paper, we reassess the role of common variants on epilepsy using extensions of heritability analysis. Our data set consists of 1258 UK patients with epilepsy, of which 958 have focal epilepsy, and 5129 population control subjects, with genotypes recorded for over 4 million common single nucleotide polymorphisms. Firstly, we show that on the liability scale, common variants collectively explain at least 26% (standard deviation 5%) of phenotypic variation for all epilepsy and 27% (standard deviation 5%) for focal epilepsy. Secondly we provide a new method for estimating the number of causal variants for complex traits; when applied to epilepsy, our most optimistic estimate suggests that at least 400 variants influence disease susceptibility, with potentially many thousands. Thirdly, we use bivariate analysis to assess how similar the genetic architecture of focal epilepsy is to that of non-focal epilepsy; we demonstrate both significant differences (P = 0.004) and significant similarities (P = 0.01) between the two subtypes, indicating that although the clinical definition of focal epilepsy does identify a genetically distinct epilepsy subtype, there is also scope to improve the classification of epilepsy by incorporating genotypic information. Lastly, we investigate the potential value in using genetic data to diagnose epilepsy following a single epileptic seizure; we find that a prediction model explaining 10% of phenotypic variation could have clinical utility for deciding which single-seizure individuals are likely to benefit from immediate anti-epileptic drug therapy.

Purification of mouse axoplasmic proteins from dorsal root ganglia nerves for proteomics analysis.

The study of neuronal signaling ex vivo requires the identification of the proteins that are represented within the neuronal axoplasm. Here, we describe a detailed protocol to isolate the axoplasm of peripheral and central axonal branches of sciatic dorsal root ganglia neurons in mice. The axoplasm is separated by 2D gel and digestion followed by proteomics analysis with MS/MS-LC. This protocol can be applied to dissect the axoplasmic protein expression signatures before and after a sciatic nerve or a spinal cord injury. For complete details on the use and execution of this protocol, please refer to Kong et al. (2020).

Fluorescence-Activated Nuclei Negative Sorting of Neurons Combined with Single Nuclei RNA Sequencing to Study the Hippocampal Neurogenic Niche.

Adult Hippocampal Neurogenesis (AHN), which consists of a lifelong maintenance of proliferative and quiescent neural stem cells (NSCs) within the sub-granular zone (SGZ) of the dentate gyrus (DG) and their differentiation from newly born neurons into granule cells in the granule cell layer, is well validated across numerous studies. Using genetically modified animals, particularly rodents, is a valuable tool to investigate signaling pathways regulating AHN and to study the role of each cell type that compose the hippocampal neurogenic niche. To address the latter, methods combining single nuclei isolation with next generation sequencing have had a significant impact in the field of AHN to identify gene signatures for each cell population. Further refinement of these techniques is however needed to phenotypically profile rarer cell populations within the DG. Here, we present a method that utilizes Fluorescence Activated Nuclei Sorting (FANS) to exclude most neuronal populations from a single nuclei suspension isolated from freshly dissected DG, by selecting unstained nuclei for the NeuN antigen, in order to perform single nuclei RNA sequencing (snRNA-seq). This method is a potential steppingstone to further investigate intercellular regulation of the AHN and to uncover novel cellular markers and mechanisms across species.

AMPK controls the axonal regenerative ability of dorsal root ganglia sensory neurons after spinal cord injury.

Regeneration after injury occurs in axons that lie in the peripheral nervous system but fails in the central nervous system, thereby limiting functional recovery. Differences in axonal signalling in response to injury that might underpin this differential regenerative ability are poorly characterized. Combining axoplasmic proteomics from peripheral sciatic or central projecting dorsal root ganglion (DRG) axons with cell body RNA-seq, we uncover injury-dependent signalling pathways that are uniquely represented in peripheral versus central projecting sciatic DRG axons. We identify AMPK as a crucial regulator of axonal regenerative signalling that is specifically downregulated in injured peripheral, but not central, axons. We find that AMPK in DRG interacts with the 26S proteasome and its CaMKIIα-dependent regulatory subunit PSMC5 to promote AMPKα proteasomal degradation following sciatic axotomy. Conditional deletion of AMPKα1 promotes multiple regenerative signalling pathways after central axonal injury and stimulates robust axonal growth across the spinal cord injury site, suggesting inhibition of AMPK as a therapeutic strategy to enhance regeneration following spinal cord injury.

Integrated systems-genetic analyses reveal a network target for delaying glioma progression.

OBJECTIVE: To identify a convergent, multitarget proliferation characteristic for astrocytoma transformation that could be targeted for therapy discovery. METHODS: Using an integrated functional genomics approach, we prioritized networks associated with astrocytoma progression using the following criteria: differential co-expression between grade II and grade III IDH1-mutated and 1p/19q euploid astrocytomas, preferential enrichment for genetic risk to cancer, association with patient survival and sample-level genomic features. Drugs targeting the identified multitarget network characteristic for astrocytoma transformation were computationally predicted using drug transcriptional perturbation data and validated using primary human astrocytoma cells. RESULTS: A single network, M2, consisting of 177 genes, was associated with glioma progression on the basis of the above criteria. Functionally, M2 encoded physically interacting proteins regulating cell cycle processes and analysis of genome-wide gene-regulatory interactions using mutual information and DNA-protein interactions revealed the known regulators of cell cycle processes FoxM1, B-Myb, and E2F2 as key regulators of M2. These results suggest functional disruption of M2 via gene mutation or altered expression as a convergent pathway regulating astrocytoma transformation. By considering M2 as a multitarget drug target regulating astrocytoma transformation, we identified several drugs that are predicted to restore M2 expression in anaplastic astrocytoma toward its low-grade profile and of these, we validated the known antiproliferative drug resveratrol as down-regulating multiple nodes of M2 including at nanomolar concentrations achievable in human cerebrospinal fluid by oral dosing. INTERPRETATION: Our results identify M2 as a multitarget network characteristic for astrocytoma progression and encourage M2-based drug screening to identify new compounds for preventing glioma transformation.

Epigenomic signatures underpin the axonal regenerative ability of dorsal root ganglia sensory neurons.

Axonal injury results in regenerative success or failure, depending on whether the axon lies in the peripheral or the CNS, respectively. The present study addresses whether epigenetic signatures in dorsal root ganglia discriminate between regenerative and non-regenerative axonal injury. Chromatin immunoprecipitation for the histone 3 (H3) post-translational modifications H3K9ac, H3K27ac and H3K27me3; an assay for transposase-accessible chromatin; and RNA sequencing were performed in dorsal root ganglia after sciatic nerve or dorsal column axotomy. Distinct histone acetylation and chromatin accessibility signatures correlated with gene expression after peripheral, but not central, axonal injury. DNA-footprinting analyses revealed new transcriptional regulators associated with regenerative ability. Machine-learning algorithms inferred the direction of most of the gene expression changes. Neuronal conditional deletion of the chromatin remodeler CCCTC-binding factor impaired nerve regeneration, implicating chromatin organization in the regenerative competence. Altogether, the present study offers the first epigenomic map providing insight into the transcriptional response to injury and the differential regenerative ability of sensory neurons.

Author Correction: WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling.

Cardiac fibrosis is a final common pathology in inherited and acquired heart diseases that causes cardiac electrical and pump failure. Here, we use systems genetics to identify a pro-fibrotic gene network in the diseased heart and show that this network is regulated by the E3 ubiquitin ligase WWP2, specifically by the WWP2-N terminal isoform. Importantly, the WWP2-regulated pro-fibrotic gene network is conserved across different cardiac diseases characterized by fibrosis: human and murine dilated cardiomyopathy and repaired tetralogy of Fallot. Transgenic mice lacking the N-terminal region of the WWP2 protein show improved cardiac function and reduced myocardial fibrosis in response to pressure overload or myocardial infarction. In primary cardiac fibroblasts, WWP2 positively regulates the expression of pro-fibrotic markers and extracellular matrix genes. TGFß1 stimulation promotes nuclear translocation of the WWP2 isoforms containing the N-terminal region and their interaction with SMAD2. WWP2 mediates the TGFß1-induced nucleocytoplasmic shuttling and transcriptional activity of SMAD2.

A systems-level framework for drug discovery identifies Csf1R as an anti-epileptic drug target.

The identification of drug targets is highly challenging, particularly for diseases of the brain. To address this problem, we developed and experimentally validated a general computational framework for drug target discovery that combines gene regulatory information with causal reasoning ("Causal Reasoning Analytical Framework for Target discovery"-CRAFT). Using a systems genetics approach and starting from gene expression data from the target tissue, CRAFT provides a predictive framework for identifying cell membrane receptors with a direction-specified influence over disease-related gene expression profiles. As proof of concept, we applied CRAFT to epilepsy and predicted the tyrosine kinase receptor Csf1R as a potential therapeutic target. The predicted effect of Csf1R blockade in attenuating epilepsy seizures was validated in three pre-clinical models of epilepsy. These results highlight CRAFT as a systems-level framework for target discovery and suggest Csf1R blockade as a novel therapeutic strategy in epilepsy. CRAFT is applicable to disease settings other than epilepsy.

Rare and common epilepsies converge on a shared gene regulatory network providing opportunities for novel antiepileptic drug discovery.

BACKGROUND: The relationship between monogenic and polygenic forms of epilepsy is poorly understood and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy. RESULTS: We identified a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy and for common variants associated with polygenic epilepsy. The genes in the M30 network are expressed widely in the human brain under tight developmental control and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within the M30 network were preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent downregulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the downregulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons. CONCLUSIONS: Taken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy.

Systems genetics identifies a convergent gene network for cognition and neurodevelopmental disease.

Genetic determinants of cognition are poorly characterized, and their relationship to genes that confer risk for neurodevelopmental disease is unclear. Here we performed a systems-level analysis of genome-wide gene expression data to infer gene-regulatory networks conserved across species and brain regions. Two of these networks, M1 and M3, showed replicable enrichment for common genetic variants underlying healthy human cognitive abilities, including memory. Using exome sequence data from 6,871 trios, we found that M3 genes were also enriched for mutations ascertained from patients with neurodevelopmental disease generally, and intellectual disability and epileptic encephalopathy in particular. M3 consists of 150 genes whose expression is tightly developmentally regulated, but which are collectively poorly annotated for known functional pathways. These results illustrate how systems-level analyses can reveal previously unappreciated relationships between neurodevelopmental disease-associated genes in the developed human brain, and provide empirical support for a convergent gene-regulatory network influencing cognition and neurodevelopmental disease.

Systems genetics identifies Sestrin 3 as a regulator of a proconvulsant gene network in human epileptic hippocampus.

Gene-regulatory network analysis is a powerful approach to elucidate the molecular processes and pathways underlying complex disease. Here we employ systems genetics approaches to characterize the genetic regulation of pathophysiological pathways in human temporal lobe epilepsy (TLE). Using surgically acquired hippocampi from 129 TLE patients, we identify a gene-regulatory network genetically associated with epilepsy that contains a specialized, highly expressed transcriptional module encoding proconvulsive cytokines and Toll-like receptor signalling genes. RNA sequencing analysis in a mouse model of TLE using 100 epileptic and 100 control hippocampi shows the proconvulsive module is preserved across-species, specific to the epileptic hippocampus and upregulated in chronic epilepsy. In the TLE patients, we map the trans-acting genetic control of this proconvulsive module to Sestrin 3 (SESN3), and demonstrate that SESN3 positively regulates the module in macrophages, microglia and neurons. Morpholino-mediated Sesn3 knockdown in zebrafish confirms the regulation of the transcriptional module, and attenuates chemically induced behavioural seizures in vivo.