Informing disease modelling with brain-relevant functional genomic annotations

The past decade has seen a surge in the number of disease/trait-associated variants, largely because of the union of studies to share genetic data and the availability of electronic health records from large cohorts for research use. Variant discovery for neurological and neuropsychiatric genome-wide association studies, including schizophrenia, Parkinson's disease and Alzheimer's disease, has greatly benefitted; however, the translation of these genetic association results to interpretable biological mechanisms and models is lagging. Interpreting disease-associated variants requires knowledge of gene regulatory mechanisms and computational tools that permit integration of this knowledge with genome-wide association study results. Here, we summarize key conceptual advances in the generation of brain-relevant functional genomic annotations and amongst tools that allow integration of these annotations with association summary statistics, which together provide a new and exciting opportunity to identify disease-relevant genes, pathways and cell types in silico. We discuss the opportunities and challenges associated with these developments and conclude with our perspective on future advances in annotation generation, tool development and the union of the two

An additional k-means clustering step improves the biological features of WGCNA gene co-expression networks

Background: Weighted Gene Co-expression Network Analysis (WGCNA) is a widely used R software package for the generation of gene co-expression networks (GCN). WGCNA generates both a GCN and a derived partitioning of clusters of genes (modules). We propose k-means clustering as an additional processing step to conventional WGCNA, which we have implemented in the R package km2gcn (k-means to gene co-expression network, https://github.com/juanbot/km2gcn). Results: We assessed our method on networks created from UKBEC data (10 different human brain tissues), on networks created from GTEx data (42 human tissues, including 13 brain tissues), and on simulated networks derived from GTEx data. We observed substantially improved module properties, including: (1) few or zero misplaced genes; (2) increased counts of replicable clusters in alternate tissues (x3.1 on average); (3) improved enrichment of Gene Ontology terms (seen in 48/52 GCNs) (4) improved cell type enrichment signals (seen in 21/23 brain GCNs); and (5) more accurate partitions in simulated data according to a range of similarity indices. Conclusions: The results obtained from our investigations indicate that our k-means method, applied as an adjunct to standard WGCNA, results in better network partitions. These improved partitions enable more fruitful downstream analyses, as gene modules are more biologically meaningful

Frontotemporal dementia: insights into the biological underpinnings of disease through gene co-expression network analysis

BACKGROUND: In frontotemporal dementia (FTD) there is a critical lack in the understanding of biological and molecular mechanisms involved in disease pathogenesis. The heterogeneous genetic features associated with FTD suggest that multiple disease-mechanisms are likely to contribute to the development of this neurodegenerative condition. We here present a systems biology approach with the scope of i) shedding light on the biological processes potentially implicated in the pathogenesis of FTD and ii) identifying novel potential risk factors for FTD. We performed a gene co-expression network analysis of microarray expression data from 101 individuals without neurodegenerative diseases to explore regional-specific co-expression patterns in the frontal and temporal cortices for 12 genes (MAPT, GRN, CHMP2B, CTSC, HLA-DRA, TMEM106B, C9orf72, VCP, UBQLN2, OPTN, TARDBP and FUS) associated with FTD and we then carried out gene set enrichment and pathway analyses, and investigated known protein-protein interactors (PPIs) of FTD-genes products. RESULTS: Gene co-expression networks revealed that several FTD-genes (such as MAPT and GRN, CTSC and HLA-DRA, TMEM106B, and C9orf72, VCP, UBQLN2 and OPTN) were clustering in modules of relevance in the frontal and temporal cortices. Functional annotation and pathway analyses of such modules indicated enrichment for: i) DNA metabolism, i.e. transcription regulation, DNA protection and chromatin remodelling (MAPT and GRN modules); ii) immune and lysosomal processes (CTSC and HLA-DRA modules), and; iii) protein meta/catabolism (C9orf72, VCP, UBQLN2 and OPTN, and TMEM106B modules). PPI analysis supported the results of the functional annotation and pathway analyses. CONCLUSIONS: This work further characterizes known FTD-genes and elaborates on their biological relevance to disease: not only do we indicate likely impacted regional-specific biological processes driven by FTD-genes containing modules, but also do we suggest novel potential risk factors among the FTD-genes interactors as targets for further mechanistic characterization in hypothesis driven cell biology work

CoExp: A Web Tool for the Exploitation of Co-expression Networks

Gene co-expression networks are a powerful type of analysis to construct gene groupings based on transcriptomic profiling. Co-expression networks make it possible to discover modules of genes whose mRNA levels are highly correlated across samples. Subsequent annotation of modules often reveals biological functions and/or evidence of cellular specificity for cell types implicated in the tissue being studied. There are multiple ways to perform such analyses with weighted gene co-expression network analysis (WGCNA) amongst one of the most widely used R packages. While managing a few network models can be done manually, it is often more advantageous to study a wider set of models derived from multiple independently generated transcriptomic data sets (e.g., multiple networks built from many transcriptomic sources). However, there is no software tool available that allows this to be easily achieved. Furthermore, the visual nature of co-expression networks in combination with the coding skills required to explore networks, makes the construction of a web-based platform for their management highly desirable. Here, we present the CoExp Web application, a user-friendly online tool that allows the exploitation of the full collection of 109 co-expression networks provided by the CoExpNets suite of R packages. We describe the usage of CoExp, including its contents and the functionality available through the family of CoExpNets packages. All the tools presented, including the web front- and back-ends are available for the research community so any research group can build its own suite of networks and make them accessible through their own CoExp Web application. Therefore, this paper is of interest to both researchers wishing to annotate their genes of interest across different brain network models and specialists interested in the creation of GCNs looking for a tool to appropriately manage, use, publish, and share their networks in a consistent and productive manner

P2 nucleotide receptors on C2C12 satellite cells

In developing muscle cells environmental stimuli transmitted by purines binding to the specific receptors are crucial proliferation regulators. C2C12 myoblasts express numerous purinergic receptors representing both main classes: P2X and P2Y. Among P2Y receptors we have found the expression of P2Y1, P2Y2, P2Y4, P2Y6 and P2Y12 family members while among P2X receptors P2X4, P2X5 and P2X7 were discovered. We have been able to show that activation of those receptors is responsible for ERK class kinase activity, responsible for regulation of cell proliferation pathway. We have also demonstrated that this activity is calcium dependent suggesting Ca2+ ions as secondary messenger between receptor and kinase regulatory system. More specifically, we do suspect that in C2C12 myoblasts calcium channels of P2X receptors, particularly P2X5 play the main role in proliferation regulation. In further development of myoblasts into myotubes, when proliferation is gradually inhibited, the pattern of P2 receptors is changed. This phenomenon is followed by diminishing of the P2Y2-dependent Ca2+ signaling, while the mRNA expression of P2Y2 receptor reminds still on the high level. Moreover, P2X2 receptor mRNA, absent in myoblasts appears in myotubes. These data show that differentiation of C2C12 cell line satellite myoblasts is accompanied by changes in P2 receptors expression pattern

Heritability enrichment implicates microglia in Parkinson's disease pathogenesis

OBJECTIVE: Understanding how different parts of the immune system contribute to pathogenesis in Parkinson's disease is a burning challenge with important therapeutic implications. We studied enrichment of common variant heritability for Parkinson's disease stratified by immune and brain cell types. METHODS: We used summary statistics from the most recent meta-analysis of genome-wide association studies in Parkinson's disease and partitioned heritability using linkage disequilibrium score regression, stratified for specific cell types as defined by open chromatin regions. We also validated enrichment results using a polygenic risk score approach and intersected disease-associated variants with epigenetic data and expression quantitative loci to nominate and explore a putative microglial locus. RESULTS: We found significant enrichment of Parkinson's disease risk heritability in open chromatin regions of microglia and monocytes. Genomic annotations overlapped substantially between these two cell types, and only the enrichment signal for microglia remained significant in a joint model. We present evidence suggesting P2RY12, a key microglial gene and target for the anti-thrombotic agent clopidogrel, as the likely driver of a significant Parkinson's disease association signal on chromosome 3. INTERPRETATION: Our results provide further support for the importance of immune mechanisms in PD pathogenesis, highlight microglial dysregulation as a contributing etiological factor and nominate a targetable microglial gene candidate as a pathogenic player. Immune processes can be modulated by therapy, with potentially important clinical implications for future treatment in Parkinson's disease

Gene co-expression networks shed light into diseases of brain iron accumulation

Aberrant brain iron deposition is observed in both common and rare neurodegenerative disorders, including those categorized as Neurodegeneration with Brain Iron Accumulation (NBIA), which are characterized by focal iron accumulation in the basal ganglia. Two NBIA genes are directly involved in iron metabolism, but whether other NBIA-related genes also regulate iron homeostasis in the human brain, and whether aberrant iron deposition contributes to neurodegenerative processes remains largely unknown. This study aims to expand our understanding of these iron overload diseases and identify relationships between known NBIA genes and their main interacting partners by using a systems biology approach. We used whole-transcriptome gene expression data from human brain samples originating from 101 neuropathologically normal individuals (10 brain regions) to generate weighted gene co-expression networks and cluster the 10 known NBIA genes in an unsupervised manner. We investigated NBIA-enriched networks for relevant cell types and pathways, and whether they are disrupted by iron loading in NBIA diseased tissue and in an in vivo mouse model. We identified two basal ganglia gene co-expression modules significantly enriched for NBIA genes, which resemble neuronal and oligodendrocytic signatures. These NBIA gene networks are enriched for iron-related genes, and implicate synapse and lipid metabolism related pathways. Our data also indicates that these networks are disrupted by excessive brain iron loading. We identified multiple cell types in the origin of NBIA disorders. We also found unforeseen links between NBIA networks and iron-related processes, and demonstrate convergent pathways connecting NBIAs and phenotypically overlapping diseases. Our results are of further relevance for these diseases by providing candidates for new causative genes and possible points for therapeutic intervention

Brain iron accumulation affects myelin-related molecular systems implicated in a rare neurogenetic disease family with neuropsychiatric features

The ‘neurodegeneration with brain iron accumulation’ (NBIA) disease family entails movement or cognitive impairment, often with psychiatric features. To understand how iron loading affects the brain, we studied mice with disruption of two iron regulatory genes, hemochromatosis (Hfe) and transferrin receptor 2 (Tfr2). Inductively coupled plasma atomic emission spectroscopy demonstrated increased iron in the Hfe-/- × Tfr2mut brain (P=0.002, n =5/group), primarily localized by Perls’ staining to myelinated structures. Western immunoblotting showed increases of the iron storage protein ferritin light polypeptide and microarray and real-time reverse transcription-PCR revealed decreased transcript levels (P<0.04, n =5/group) for five other NBIA genes, phospholipase A2 group VI, fatty acid 2-hydroxylase, ceruloplasmin, chromosome 19 open reading frame 12 and ATPase type 13A2. Apart from the ferroxidase ceruloplasmin, all are involved in myelin homeostasis; 16 other myelin-related genes also showed reduced expression (P<0.05), although gross myelin structure and integrity appear unaffected (P>0.05). Overlap (P<0.0001) of differentially expressed genes in Hfe-/- × Tfr2mut brain with human gene co-expression networks suggests iron loading influences expression of NBIA-related and myelin-related genes co-expressed in normal human basal ganglia. There was overlap (P<0.0001) of genes differentially expressed in Hfe-/- × Tfr2mut brain and post-mortem NBIA basal ganglia. Hfe-/- × Tfr2mut mice were hyperactive (P<0.0112) without apparent cognitive impairment by IntelliCage testing (P>0.05). These results implicate myelin-related systems involved in NBIA neuropathogenesis in early responses to iron loading. This may contribute to behavioral symptoms in NBIA and hemochromatosis and is relevant to patients with abnormal iron status and psychiatric disorders involving myelin abnormalities or resistant to conventional treatments

Human-lineage-specific genomic elements are associated with neurodegenerative disease and APOE transcript usage.

Knowledge of genomic features specific to the human lineage may provide insights into brain-related diseases. We leverage high-depth whole genome sequencing data to generate a combined annotation identifying regions simultaneously depleted for genetic variation (constrained regions) and poorly conserved across primates. We propose that these constrained, non-conserved regions (CNCRs) have been subject to human-specific purifying selection and are enriched for brain-specific elements. We find that CNCRs are depleted from protein-coding genes but enriched within lncRNAs. We demonstrate that per-SNP heritability of a range of brain-relevant phenotypes are enriched within CNCRs. We find that genes implicated in neurological diseases have high CNCR density, including APOE, highlighting an unannotated intron-3 retention event. Using human brain RNA-sequencing data, we show the intron-3-retaining transcript to be more abundant in Alzheimer's disease with more severe tau and amyloid pathological burden. Thus, we demonstrate potential association of human-lineage-specific sequences in brain development and neurological disease

Mutations in the autoregulatory domain of β-tubulin 4a cause hereditary dystonia.

Dystonia type 4 (DYT4) was first described in a large family from Heacham in Norfolk with an autosomal dominantly inherited whispering dysphonia, generalized dystonia, and a characteristic hobby horse ataxic gait. We carried out a genetic linkage analysis in the extended DYT4 family that spanned 7 generations from England and Australia, revealing a single LOD score peak of 6.33 on chromosome 19p13.12-13. Exome sequencing in 2 cousins identified a single cosegregating mutation (p.R2G) in the β-tubulin 4a (TUBB4a) gene that was absent in a large number of controls. The mutation is highly conserved in the β-tubulin autoregulatory MREI (methionine-arginine-glutamic acid-isoleucine) domain, highly expressed in the central nervous system, and extensive in vitro work has previously demonstrated that substitutions at residue 2, specifically R2G, disrupt the autoregulatory capability of the wild-type β-tubulin peptide, affirming the role of the cytoskeleton in dystonia pathogenesis