Cell Type-Specific Intralocus Interactions Reveal Oligodendrocyte Mechanisms in MS.

Multiple sclerosis (MS) is an autoimmune disease characterized by attack on oligodendrocytes within the central nervous system (CNS). Despite widespread use of immunomodulatory therapies, patients may still face progressive disability because of failure of myelin regeneration and loss of neurons, suggesting additional cellular pathologies. Here, we describe a general approach for identifying specific cell types in which a disease allele exerts a pathogenic effect. Applying this approach to MS risk loci, we pinpoint likely pathogenic cell types for 70%. In addition to T cell loci, we unexpectedly identified myeloid- and CNS-specific risk loci, including two sites that dysregulate transcriptional pause release in oligodendrocytes. Functional studies demonstrated inhibition of transcriptional elongation is a dominant pathway blocking oligodendrocyte maturation. Furthermore, pause release factors are frequently dysregulated in MS brain tissue. These data implicate cell-intrinsic aberrations outside of the immune system and suggest new avenues for therapeutic development. VIDEO ABSTRACT.

Evolution of regulatory signatures in primate cortical neurons at cell-type resolution.

The human cerebral cortex contains many cell types that likely underwent independent functional changes during evolution. However, cell-type-specific regulatory landscapes in the cortex remain largely unexplored. Here we report epigenomic and transcriptomic analyses of the two main cortical neuronal subtypes, glutamatergic projection neurons and GABAergic interneurons, in human, chimpanzee, and rhesus macaque. Using genome-wide profiling of the H3K27ac histone modification, we identify neuron-subtype-specific regulatory elements that previously went undetected in bulk brain tissue samples. Human-specific regulatory changes are uncovered in multiple genes, including those associated with language, autism spectrum disorder, and drug addiction. We observe preferential evolutionary divergence in neuron subtype-specific regulatory elements and show that a substantial fraction of pan-neuronal regulatory elements undergoes subtype-specific evolutionary changes. This study sheds light on the interplay between regulatory evolution and cell-type-dependent gene-expression programs, and provides a resource for further exploration of human brain evolution and function.

Molecular windows into the human brain for psychiatric disorders.

Delineating the pathophysiology of psychiatric disorders has been extremely challenging but technological advances in recent decades have facilitated a deeper interrogation of molecular processes in the human brain. Initial candidate gene expression studies of the postmortem brain have evolved into genome wide profiling of the transcriptome and the epigenome, a critical regulator of gene expression. Here, we review the potential and challenges of direct molecular characterization of the postmortem human brain, and provide a brief overview of recent transcriptional and epigenetic studies with respect to neuropsychiatric disorders. Such information can now be leveraged and integrated with the growing number of genome-wide association databases to provide a functional context of trait-associated genetic variants linked to psychiatric illnesses and related phenotypes. While it is clear that the field is still developing and challenges remain to be surmounted, these recent advances nevertheless hold tremendous promise for delineating the neurobiological underpinnings of mental diseases and accelerating the development of novel medication strategies.

Substantial DNA methylation differences between two major neuronal subtypes in human brain.

The brain is built from a large number of cell types which have been historically classified using location, morphology and molecular markers. Recent research suggests an important role of epigenetics in shaping and maintaining cell identity in the brain. To elucidate the role of DNA methylation in neuronal differentiation, we developed a new protocol for separation of nuclei from the two major populations of human prefrontal cortex neurons--GABAergic interneurons and glutamatergic (GLU) projection neurons. Major differences between the neuronal subtypes were revealed in CpG, non-CpG and hydroxymethylation (hCpG). A dramatically greater number of undermethylated CpG sites in GLU versus GABA neurons were identified. These differences did not directly translate into differences in gene expression and did not stem from the differences in hCpG methylation, as more hCpG methylation was detected in GLU versus GABA neurons. Notably, a comparable number of undermethylated non-CpG sites were identified in GLU and GABA neurons, and non-CpG methylation was a better predictor of subtype-specific gene expression compared to CpG methylation. Regions that are differentially methylated in GABA and GLU neurons were significantly enriched for schizophrenia risk loci. Collectively, our findings suggest that functional differences between neuronal subtypes are linked to their epigenetic specification.

A unique gene expression signature associated with serotonin 2C receptor RNA editing in the prefrontal cortex and altered in suicide.

Editing of the pre-mRNA for the serotonin receptor 2C (5-HT2CR) by site-specific adenosine deamination (A-to-I pre-mRNA editing) substantially increases the functional plasticity of this key neurotransmitter receptor and is thought to contribute to homeostatic mechanisms in neurons. 5-HT2CR mRNA editing generates up to 24 different receptor isoforms. The extent of editing correlates with 5-HT2CR functional activity: more highly edited isoforms exhibit the least function. Altered 5-HT2CR editing has been reported in postmortem brains of suicide victims. We report a comparative analysis of the connections among 5-HT2CR editing, genome-wide gene expression and DNA methylation in suicide victims, individuals with major depressive disorder and non-psychiatric controls. The results confirm previous findings of an overrepresentation of highly edited mRNA variants (which encode hypoactive 5-HT2CR receptors) in the brains of suicide victims. A large set of genes for which the expression level is associated with editing was detected. This signature set of editing-associated genes is significantly enriched for genes that are involved in synaptic transmission, genes that are preferentially expressed in neurons, and genes whose expression is correlated with the level of DNA methylation. Notably, we report that the link between 5-HT2CR editing and gene expression is disrupted in suicide victims. The results suggest that the postulated homeostatic function of 5-HT2CR editing is dysregulated in individuals who committed suicide.

Differences in DNA methylation between human neuronal and glial cells are concentrated in enhancers and non-CpG sites.

We applied Illumina Human Methylation450K array to perform a genomic-scale single-site resolution DNA methylation analysis in neuronal and nonneuronal (primarily glial) nuclei separated from the orbitofrontal cortex of postmortem human brain. The findings were validated using enhanced reduced representation bisulfite sequencing. We identified thousands of sites differentially methylated (DM) between neuronal and nonneuronal cells. The DM sites were depleted within CpG-island-containing promoters but enriched in predicted enhancers. Classification of the DM sites into those undermethylated in neurons (neuronal type) and those undermethylated in nonneuronal cells (glial type), combined with findings of others that methylation within control elements typically negatively correlates with gene expression, yielded large sets of predicted neuron-specific and non-neuron-specific genes. These sets of predicted genes were in excellent agreement with the available direct measurements of gene expression in human and mouse. We also found a distinct set of DNA methylation patterns that were unique for neuronal cells. In particular, neuronal-type differential methylation was overrepresented in CpG island shores, enriched within gene bodies but not in intergenic regions, and preferentially harbored binding motifs for a distinct set of transcription factors, including neuron-specific activity-dependent factors. Finally, non-CpG methylation was substantially more prevalent in neurons than in nonneuronal cells.

Dependencies among editing sites in serotonin 2C receptor mRNA.

The serotonin 2C receptor (5-HT(2C)R)-a key regulator of diverse neurological processes-exhibits functional variability derived from editing of its pre-mRNA by site-specific adenosine deamination (A-to-I pre-mRNA editing) in five distinct sites. Here we describe a statistical technique that was developed for analysis of the dependencies among the editing states of the five sites. The statistical significance of the observed correlations was estimated by comparing editing patterns in multiple individuals. For both human and rat 5-HT(2C)R, the editing states of the physically proximal sites A and B were found to be strongly dependent. In contrast, the editing states of sites C and D, which are also physically close, seem not to be directly dependent but instead are linked through the dependencies on sites A and B, respectively. We observed pronounced differences between the editing patterns in humans and rats: in humans site A is the key determinant of the editing state of the other sites, whereas in rats this role belongs to site B. The structure of the dependencies among the editing sites is notably simpler in rats than it is in humans implying more complex regulation of 5-HT(2C)R editing and, by inference, function in the human brain. Thus, exhaustive statistical analysis of the 5-HT(2C)R editing patterns indicates that the editing state of sites A and B is the primary determinant of the editing states of the other three sites, and hence the overall editing pattern. Taken together, these findings allow us to propose a mechanistic model of concerted action of ADAR1 and ADAR2 in 5-HT(2C)R editing. Statistical approach developed here can be applied to other cases of interdependencies among modification sites in RNA and proteins.

Divergent single cell transcriptome and epigenome alterations in ALS and FTD patients with C9orf72 mutation.

A repeat expansion in the C9orf72 (C9) gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we investigate single nucleus transcriptomics (snRNA-seq) and epigenomics (snATAC-seq) in postmortem motor and frontal cortices from C9-ALS, C9-FTD, and control donors. C9-ALS donors present pervasive alterations of gene expression with concordant changes in chromatin accessibility and histone modifications. The greatest alterations occur in upper and deep layer excitatory neurons, as well as in astrocytes. In neurons, the changes imply an increase in proteostasis, metabolism, and protein expression pathways, alongside a decrease in neuronal function. In astrocytes, the alterations suggest activation and structural remodeling. Conversely, C9-FTD donors have fewer high-quality neuronal nuclei in the frontal cortex and numerous gene expression changes in glial cells. These findings highlight a context-dependent molecular disruption in C9-ALS and C9-FTD, indicating unique effects across cell types, brain regions, and diseases.

Ionotropic glutamate receptor mRNA editing in the prefrontal cortex: no alterations in schizophrenia or bipolar disorder.

BACKGROUND: Dysfunction of glutamate neurotransmission has been implicated in the pathology of schizophrenia and bipolar disorder, and one mechanism by which glutamate signalling can be altered is through RNA editing of ionotropic glutamate receptors (iGluRs). The objectives of the present study were to evaluate the editing status of iGluRs in the human prefrontal cortex, determine whether iGluR editing is associated with psychiatric disease or suicide and evaluate a potential association between editing and alternative splicing in the α-amino-3-hydroxy-5-methylisoxazole-4-propionate (AMPA) iGluR subunits' pre-mRNA. METHODS: We studied specimens derived from patients with antemortem diagnoses of bipolar disorder (n = 31) or schizophrenia (n = 34) who died by suicide or other causes, and from psychiatrically healthy controls (n = 34) who died from causes other than suicide. The RNA editing at all 8 editing sites within AMPA (GluA2-4 subunits) and kainate (GluK1-2 subunits) iGluRs was analyzed using a novel real-time quantitative polymerase chain reaction assay. RESULTS: No differences in editing were detected among schizophrenia, bipolar or control groups or between suicide completers and patients who died from causes other than suicide. The editing efficiency was significantly higher in the flop than in the flip splicoforms of GluA3-4 AMPA subunits (all p < 0.001). LIMITATIONS: The study is limited by the near absence of specimens from medicationnaive psychiatric patients and considerable variation in medication regimens among individuals, both of which introduce considerable uncertainty into the analysis of potential medication effects. CONCLUSION: We found that iGluR RNA editing status was not associated with bipolar disorder, schizophrenia or suicide. Differences in editing between flip and flop splicoforms suggest that glutamate sensitivity of receptors containing GluA3 and/or GluA4 flop subunits is moderated as a result of increased editing.

Cellular and genetic drivers of RNA editing variation in the human brain.

Posttranscriptional adenosine-to-inosine modifications amplify the functionality of RNA molecules in the brain, yet the cellular and genetic regulation of RNA editing is poorly described. We quantify base-specific RNA editing across three major cell populations from the human prefrontal cortex: glutamatergic neurons, medial ganglionic eminence-derived GABAergic neurons, and oligodendrocytes. We identify more selective editing and hyper-editing in neurons relative to oligodendrocytes. RNA editing patterns are highly cell type-specific, with 189,229 cell type-associated sites. The cellular specificity for thousands of sites is confirmed by single nucleus RNA-sequencing. Importantly, cell type-associated sites are enriched in GTEx RNA-sequencing data, edited ~twentyfold higher than all other sites, and variation in RNA editing is largely explained by neuronal proportions in bulk brain tissue. Finally, we uncover 661,791 cis-editing quantitative trait loci across thirteen brain regions, including hundreds with cell type-associated features. These data reveal an expansive repertoire of highly regulated RNA editing sites across human brain cell types and provide a resolved atlas linking cell types to editing variation and genetic regulatory effects.

Sex Differences in the Human Brain Transcriptome of Cases With Schizophrenia.

BACKGROUND: While schizophrenia differs between males and females in the age of onset, symptomatology, and disease course, the molecular mechanisms underlying these differences remain uncharacterized. METHODS: To address questions about the sex-specific effects of schizophrenia, we performed a large-scale transcriptome analysis of RNA sequencing data from 437 controls and 341 cases from two distinct cohorts from the CommonMind Consortium. RESULTS: Analysis across the cohorts identified a reproducible gene expression signature of schizophrenia that was highly concordant with previous work. Differential expression across sex was reproducible across cohorts and identified X- and Y-linked genes, as well as those involved in dosage compensation. Intriguingly, the sex expression signature was also enriched for genes involved in neurexin family protein binding and synaptic organization. Differential expression analysis testing a sex-by-diagnosis interaction effect did not identify any genome-wide signature after multiple testing corrections. Gene coexpression network analysis was performed to reduce dimensionality from thousands of genes to dozens of modules and elucidate interactions among genes. We found enrichment of coexpression modules for sex-by-diagnosis differential expression signatures, which were highly reproducible across the two cohorts and involved a number of diverse pathways, including neural nucleus development, neuron projection morphogenesis, and regulation of neural precursor cell proliferation. CONCLUSIONS: Overall, our results indicate that the effect size of sex differences in schizophrenia gene expression signatures is small and underscore the challenge of identifying robust sex-by-diagnosis signatures, which will require future analyses in larger cohorts.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains.

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

Neuronal and glial 3D chromatin architecture informs the cellular etiology of brain disorders.

Cellular heterogeneity in the human brain obscures the identification of robust cellular regulatory networks, which is necessary to understand the function of non-coding elements and the impact of non-coding genetic variation. Here we integrate genome-wide chromosome conformation data from purified neurons and glia with transcriptomic and enhancer profiles, to characterize the gene regulatory landscape of two major cell classes in the human brain. We then leverage cell-type-specific regulatory landscapes to gain insight into the cellular etiology of several brain disorders. We find that Alzheimer's disease (AD)-associated epigenetic dysregulation is linked to neurons and oligodendrocytes, whereas genetic risk factors for AD highlighted microglia, suggesting that different cell types may contribute to disease risk, via different mechanisms. Moreover, integration of glutamatergic and GABAergic regulatory maps with genetic risk factors for schizophrenia (SCZ) and bipolar disorder (BD) identifies shared (parvalbumin-expressing interneurons) and distinct cellular etiologies (upper layer neurons for BD, and deeper layer projection neurons for SCZ). Collectively, these findings shed new light on cell-type-specific gene regulatory networks in brain disorders.

DNA methylation differences in cortical grey and white matter in schizophrenia.

Aim: Identify grey- and white-matter-specific DNA-methylation differences between schizophrenia (SCZ) patients and controls in postmortem brain cortical tissue. Materials & methods: Grey and white matter were separated from postmortem brain tissue of the superior temporal and medial frontal gyrus from SCZ (n = 10) and control (n = 11) cases. Genome-wide DNA-methylation analysis was performed using the Infinium EPIC Methylation Array (Illumina, CA, USA). Results: Four differentially methylated regions associated with SCZ status and tissue type (grey vs white matter) were identified within or near KLF9, SFXN1, SPRED2 and ALS2CL genes. Gene-expression analysis showed differential expression of KLF9 and SFXN1 in SCZ. Conclusion: Our data show distinct differences in DNA methylation between grey and white matter that are unique to SCZ, providing new leads to unravel the pathogenesis of SCZ.

Lay abstract This study investigated the way gene activity is regulated in brain cells of patients with schizophrenia (SCZ; a severe mental illness characterized by psychosis) compared with unaffected controls. The study focuses on the differences between parts of the brain with many cell bodies (grey matter) in contrast to those parts with mainly conducting fibers (white matter). For that purpose, grey and white matter were separated from brain tissue of ten individuals with SCZ and 11 without. All brains were obtained after the patients died and donated their brains to science. Array technology was used to analyze 800,000 sections of the DNA at once. The study identified regions on four genes that can turn the genes on and off differently in schizophrenic patients compared with controls, these genes were also turned on or off depending on their location either in grey or white matter. Two of these genes showed different activation in schizophrenic patients compared with controls. Overall this study identified distinct differences between grey and white matter that are unique to SCZ, providing new leads to unravel the biology of SCZ.

Common schizophrenia risk variants are enriched in open chromatin regions of human glutamatergic 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.

TaqMan-based, real-time quantitative polymerase chain reaction method for RNA editing analysis.

Abnormal adenosine to inosine (A-to-I) messenger RNA (mRNA) editing has been linked to several disease states afflicting the central nervous system. Here we report an assay to determine RNA editing frequencies at specific sites that is based on quantitative polymerase chain reaction (qPCR) with TaqMan probes. The assay was tested by measuring the frequency of the A-to-I mRNA editing at the Q/R site of the human kainate receptor subunit GluR5 and was compared with two established methods of assessing RNA editing: sequencing of individual clones and restriction analysis. The qPCR assay displayed high sensitivity and reproducibility, demonstrated exceptional discrimination between edited and unedited transcript variants, and proved to have several advantages over the other editing methods. Due to the fact that TaqMan-based qPCR technology can be easily adapted to different editing targets, the increased capabilities afforded by this new technique should facilitate various RNA editing studies that aim to elucidate the role of this process in normal physiology and in disease.

Schizophrenia and sex associated differences in the expression of neuronal and oligodendrocyte-specific genes in individual thalamic nuclei.

Considerable evidence based on the study of postmortem brain tissue suggests deficits in both neuronal and myelin systems in schizophrenia (SZ). To date, the majority of the biochemical and molecular biological studies have focused on the cerebral cortex. Most information traveling to or from the cortex is relayed or synaptically gated through the thalamus, and numerous studies suggest structural and functional abnormalities in interconnected regions of the thalamus and cortex in SZ. The present study extends our gene expression studies of neuronal and myelin systems to the thalamus. Quantitative PCR was employed to assess the expression of 10 genes in 5 divisions of the thalamus which were precisely harvested using Laser Capture Microdissection. The divisions studied were present on coronal sections at the level of the centromedian nucleus (CMN) taken from 14 schizophrenic and 16 normal control postmortem brains. The genes examined were specific for oligodendrocytes (MAG, CNP, MBP), neurons (ENO2), glutamatergic neurons (VGlut1, VGlut2, PV, CB) or GABAergic neurons (GAD65, GAD67). Expression levels for each of these markers were quantitated and compared between diagnoses, between sexes, and across nuclei. CB was much more highly expressed in the CMN in SZs compared to NCs. No other diagnosis related differences in gene expression were observed. The expression levels of CNP and MAG, but not MBP, were highly correlated with one another and both, but not MBP, were much more highly expressed in females than in males in all thalamic divisions examined. All markers were differentially expressed across nuclei.

Altered serotonin 2C receptor RNA splicing in suicide: association with editing.

In an earlier study, we showed increases in serotonin 2C receptor (5-HT2CR) pre-mRNA editing in prefrontal cortex that were specific to suicide victims irrespective of associated psychiatric diagnoses. Here we demonstrate that the ratio between the two 5-HT2CR splice variants is increased in people who committed suicide, but does not vary among the diagnostic groups. This provides further evidence for suicide-specific neurobiology and suggests that, as it was previously shown in vitro, 5-HT2CR editing modulates its splicing in human brain. The association analysis indicates, however, that the efficiency of 5-HT2CR editing is an imperfect predictor of the splicing outcome, and that splice site selection is only partially controlled by the level of editing in vivo.

Ionotropic glutamate receptor mRNA expression in the human thalamus: absence of change in schizophrenia.

Abnormalities in glutamate neurotransmission are thought to be among the major contributing factors to the pathophysiology of schizophrenia. Although schizophrenia has been regarded mostly as a disorder of higher cortical function, the cortex and thalamus work as a functional unit. Existing data regarding alterations of glutamate receptor subunit expression in the thalamus in schizophrenia remain equivocal. This postmortem study examined mRNA expression of ionotropic glutamate receptor (iGluR) subunits and PSD95 in 5 precisely defined and dissected thalamic subdivisions (medial and lateral sectors of the mediodorsal nucleus; and the ventral lateral posterior, ventral posterior, and centromedian nuclei) of persons with schizophrenia and matched controls using quantitative PCR with normalization to multiple endogenous controls. Among 15 genes examined (NR1 and NR2A-D subunits of the NMDA receptor; GluR1-4 subunits of the AMPA receptor; GluR5-7 and KA1-2 subunits of the kainate receptor; PSD95), all but two (GluR4 and KA1) were expressed at quantifiable levels. Differences in iGluR gene expression were seen between different thalamic nuclei but not between diagnostic groups. The relative abundance of transcripts was: NR1>>NR2A>NR2B>NR2D>NR2C for NMDA, GluR2>GluR1>GluR3 for AMPA, and KA2>GluR5>GluR7>GluR6 for kainate receptors. The expression of PSD95 correlated with the expression of NR1, NR2A, NR2B, NR2D and GluR6 in all nuclei. These results provide detailed and quantitative information on iGluR subunit expression in multiple nuclei of the human thalamus but suggest that alterations in their expression are not a prominent feature of schizophrenia.

Brain Cell Type Specific Gene Expression and Co-expression Network Architectures.

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.