Elucidating brain cell type specific gene expression patterns is critical towards a better understanding of how cell-cell communications may influence brain functions and dysfunctions. We set out to ...compare and contrast five human and murine cell type-specific transcriptome-wide RNA expression data sets that were generated within the past several years. We defined three measures of brain cell type-relative expression including specificity, enrichment, and absolute expression and identified corresponding consensus brain cell "signatures," which were well conserved across data sets. We validated that the relative expression of top cell type markers are associated with proxies for cell type proportions in bulk RNA expression data from postmortem human brain samples. We further validated novel marker genes using an orthogonal ATAC-seq dataset. We performed multiscale coexpression network analysis of the single cell data sets and identified robust cell-specific gene modules. To facilitate the use of the cell type-specific genes for cell type proportion estimation and deconvolution from bulk brain gene expression data, we developed an R package, BRETIGEA. In summary, we identified a set of novel brain cell consensus signatures and robust networks from the integration of multiple datasets and therefore transcend limitations related to technical issues characteristic of each individual study.
Genome-wide association studies (GWAS) have identified more than 40 loci associated with Alzheimer's disease (AD), but the causal variants, regulatory elements, genes and pathways remain largely ...unknown, impeding a mechanistic understanding of AD pathogenesis. Previously, we showed that AD risk alleles are enriched in myeloid-specific epigenomic annotations. Here, we show that they are specifically enriched in active enhancers of monocytes, macrophages and microglia. We integrated AD GWAS with myeloid epigenomic and transcriptomic datasets using analytical approaches to link myeloid enhancer activity to target gene expression regulation and AD risk modification. We identify AD risk enhancers and nominate candidate causal genes among their likely targets (including AP4E1, AP4M1, APBB3, BIN1, MS4A4A, MS4A6A, PILRA, RABEP1, SPI1, TP53INP1, and ZYX) in twenty loci. Fine-mapping of these enhancers nominates candidate functional variants that likely modify AD risk by regulating gene expression in myeloid cells. In the MS4A locus we identified a single candidate functional variant and validated it in human induced pluripotent stem cell (hiPSC)-derived microglia and brain. Taken together, this study integrates AD GWAS with multiple myeloid genomic datasets to investigate the mechanisms of AD risk alleles and nominates candidate functional variants, regulatory elements and genes that likely modulate disease susceptibility.
•Angiogenesis guides development and maintenance cerebral vasculature and assures adequate regional blood flow and normal brain function.•Multifaceted research implicates cerebral microvasculature as ...a potential contributor to the pathophysiology of schizophrenia.•Atypical angiogenesis gene expression signature offers support for an anomalous blood brain barrier function in schizophrenia.•The neurobehavioral consequences of cerebral angiogenesis defects in animal models relevant to schizophrenia are discussed.
Schizophrenia (SZ) is a severe mental disorder with unknown etiology and elusive neuropathological and neurobiological features have been a focus of many theoretical hypotheses and empirical studies. Current genetic and neurobiology information relevant to SZ implicates neuronal developmental and synaptic plasticity abnormalities, and neurotransmitter, microglial and oligodendrocytes dysfunction. Several recent theories have highlighted the neurovascular unit as a potential contributor to the pathophysiology of SZ. We explored the biological plausibility of a link between SZ and the neurovascular system by examining insights gained from genetic, neuroimaging and postmortem studies, which include gene expression and neuropathology analyses. We also reviewed information from animal models of cerebral angiogenesis in order to understand better the complex interplay between angiogenic and neurotrophic factors in development, vascular endothelium/blood brain barrier remodeling and maintenance, all of which contribute to sustaining adequate regional blood flow and safeguarding normal brain function. Microvascular and hemodynamic alterations in SZ highlight the importance of further research and reveal the neurovascular unit as a potential therapeutic target in SZ.
Abstract
Motivation
Most genetic variants implicated in complex diseases by genome-wide association studies (GWAS) are non-coding, making it challenging to understand the causative genes involved in ...disease. Integrating external information such as quantitative trait locus (QTL) mapping of molecular traits (e.g. expression, methylation) is a powerful approach to identify the subset of GWAS signals explained by regulatory effects. In particular, expression QTLs (eQTLs) help pinpoint the responsible gene among the GWAS regions that harbor many genes, while methylation QTLs (mQTLs) help identify the epigenetic mechanisms that impact gene expression which in turn affect disease risk. In this work, we propose multiple-trait-coloc (moloc), a Bayesian statistical framework that integrates GWAS summary data with multiple molecular QTL data to identify regulatory effects at GWAS risk loci.
Results
We applied moloc to schizophrenia (SCZ) and eQTL/mQTL data derived from human brain tissue and identified 52 candidate genes that influence SCZ through methylation. Our method can be applied to any GWAS and relevant functional data to help prioritize disease associated genes.
Availability and implementation: moloc is available for download as an R package (https://github.com/clagiamba/moloc). We also developed a web site to visualize the biological findings (icahn.mssm.edu/moloc). The browser allows searches by gene, methylation probe and scenario of interest.
Supplementary information
Supplementary data are available at Bioinformatics online.
SARS-CoV-2 infects less than 1% of cells in the human body, yet it can cause severe damage in a variety of organs. Thus, deciphering the non-cell-autonomous effects of SARS-CoV-2 infection is ...imperative for understanding the cellular and molecular disruption it elicits. Neurological and cognitive defects are among the least understood symptoms of COVID-19 patients, with olfactory dysfunction being their most common sensory deficit. Here, we show that both in humans and hamsters, SARS-CoV-2 infection causes widespread downregulation of olfactory receptors (ORs) and of their signaling components. This non-cell-autonomous effect is preceded by a dramatic reorganization of the neuronal nuclear architecture, which results in dissipation of genomic compartments harboring OR genes. Our data provide a potential mechanism by which SARS-CoV-2 infection alters the cellular morphology and the transcriptome of cells it cannot infect, offering insight to its systemic effects in olfaction and beyond.
Display omitted
•Downregulation of odor detection pathways may explain COVID-19-induced anosmia•SARS-CoV-2-mediated disruption of nuclear architecture may impair odor detection•SARS-CoV-2-mediated nuclear reorganization is non-cell autonomous
SARS-CoV-2 induces non-cell-autonomous effects in olfactory epithelium that disrupts nuclear architecture and downregulates olfactory receptor expression in olfactory sensory neurons.
The chromatin landscape of human brain cells encompasses key information to understanding brain function. Here we use ATAC-seq to profile the chromatin structure in four distinct populations of cells ...(glutamatergic neurons, GABAergic neurons, oligodendrocytes, and microglia/astrocytes) from three different brain regions (anterior cingulate cortex, dorsolateral prefrontal cortex, and primary visual cortex) in human postmortem brain samples. We find that chromatin accessibility varies greatly by cell type and, more moderately, by brain region, with glutamatergic neurons showing the largest regional variability. Transcription factor footprinting implicates cell-specific transcriptional regulators and infers cell-specific regulation of protein-coding genes, long intergenic noncoding RNAs and microRNAs. In vivo transgenic mouse experiments validate the cell type specificity of several of these human-derived regulatory sequences. We find that open chromatin regions in glutamatergic neurons are enriched for neuropsychiatric risk variants, particularly those associated with schizophrenia. Integration of cell-specific chromatin data with a bulk tissue study of schizophrenia brains increases statistical power and confirms that glutamatergic neurons are most affected. These findings illustrate the utility of studying the cell-type-specific epigenome in complex tissues like the human brain, and the potential of such approaches to better understand the genetic basis of human brain function.
Multiple lines of evidence point to glutamatergic signaling in the postsynaptic density (PSD) as a pathophysiologic mechanism in schizophrenia. Integral to PSD glutamatergic signaling is reciprocal ...interplay between GluN and mGluR5 signaling. We examined agonist-induced mGluR5 signaling in the postmortem dorsolateral prefrontal cortex (DLPFC) derived from 17 patients and age-matched and sex-matched controls. The patient group showed a striking reduction in mGluR5 signaling, manifested by decreases in Gq/11 coupling and association with PI3K and Homer compared to controls (p < 0.01 for all). This was accompanied by increases in serine and tyrosine phosphorylation of mGluR5, which can decrease mGluR5 activity via desensitization (p < 0.01). In addition, we find altered protein-protein interaction (PPI) of mGluR5 with RGS4, norbin, Preso 1 and tamalin, which can also attenuate mGluR5 activity. We previously reported molecular underpinnings of GluN hypofunction (decreased GluN2 phosphorylation) and here we show those of reduced mGluR5 signaling in schizophrenia. We find that reduced GluN2 phosphorylation can be precipitated by attenuated mGluR5 activity and that increased mGluR5 phosphorylation can result from decreased GluN function, suggesting a reciprocal interplay between the two pathways in schizophrenia. Interestingly, the patient group showed decreased mGluR5-GluN association (p < 0.01), a mechanistic basis for the reciprocal facilitation. In sum, we present the first direct evidence for mGluR5 hypoactivity, propose a reciprocal interplay between GluN and mGluR5 pathways as integral to glutamatergic dysregulation and suggest protein-protein interactions in mGluR5-GluN complexes as potential targets for intervention in schizophrenia.
The recent implication of 108 genomic loci in schizophrenia marked a great advancement in our understanding of the disease. Against the background of its polygenic nature there is a necessity to ...identify how schizophrenia risk genes interplay. As regulators of gene expression, microRNAs (miRNAs) have repeatedly been implicated in schizophrenia etiology. It is therefore of interest to establish their role in the regulation of schizophrenia risk genes in disease-relevant biological processes.
To examine the role of miRNAs in schizophrenia in the context of disease-associated genetic variation.
The basis of this study was summary statistics from the largest schizophrenia genome-wide association study meta-analysis to date (83 550 individuals in a meta-analysis of 52 genome-wide association studies) completed in 2014 along with publicly available data for predicted miRNA targets. We examined whether schizophrenia risk genes were more likely to be regulated by miRNA. Further, we used gene set analyses to identify miRNAs that are regulators of schizophrenia risk genes.
Results from association tests for miRNA targetomes and related analyses.
In line with previous studies, we found that similar to other complex traits, schizophrenia risk genes were more likely to be regulated by miRNAs (P < 2 × 10-16). Further, the gene set analyses revealed several miRNAs regulating schizophrenia risk genes, with the strongest enrichment for targets of miR-9-5p (P = .0056 for enrichment among the top 1% most-associated single-nucleotide polymorphisms, corrected for multiple testing). It is further of note that MIR9-2 is located in a genomic region showing strong evidence for association with schizophrenia (P = 7.1 × 10-8). The second and third strongest gene set signals were seen for the targets of miR-485-5p and miR-137, respectively.
This study provides evidence for a role of miR-9-5p in the etiology of schizophrenia. Its implication is of particular interest as the functions of this neurodevelopmental miRNA tie in with established disease biology: it has a regulatory loop with the fragile X mental retardation homologue FXR1 and regulates dopamine D2 receptor density.
Most common genetic risk variants associated with neuropsychiatric disease are noncoding and are thought to exert their effects by disrupting the function of
regulatory elements (CREs), including ...promoters and enhancers. Within each cell, chromatin is arranged in specific patterns to expose the repertoire of CREs required for optimal spatiotemporal regulation of gene expression. To further understand the complex mechanisms that modulate transcription in the brain, we used frozen postmortem samples to generate the largest human brain and cell-type-specific open chromatin data set to date. Using the Assay for Transposase Accessible Chromatin followed by sequencing (ATAC-seq), we created maps of chromatin accessibility in two cell types (neurons and non-neurons) across 14 distinct brain regions of five individuals. Chromatin structure varies markedly by cell type, with neuronal chromatin displaying higher regional variability than that of non-neurons. Among our findings is an open chromatin region (OCR) specific to neurons of the striatum. When placed in the mouse, a human sequence derived from this OCR recapitulates the cell type and regional expression pattern predicted by our ATAC-seq experiments. Furthermore, differentially accessible chromatin overlaps with the genetic architecture of neuropsychiatric traits and identifies differences in molecular pathways and biological functions. By leveraging transcription factor binding analysis, we identify protein-coding and long noncoding RNAs (lncRNAs) with cell-type and brain region specificity. Our data provide a valuable resource to the research community and we provide this human brain chromatin accessibility atlas as an online database "Brain Open Chromatin Atlas (BOCA)" to facilitate interpretation.