Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

Unravelling the genetic basis of Schizophrenia

Neuronal development is a highly regulated mechanism that is central to organismal function in animals. In humans, disruptions to this process can lead to a range of neurodevelopmental phenotypes, including Schizophrenia (SCZ). SCZ has a significant genetic component, whereby an individual with an SCZ affected family member is eight times more likely to develop the disease than someone with no family history of SCZ. By examining a combination of genomic, transcriptomic and epigenomic datasets, large-scale 'omics' studies aim to delineate the relationship between genetic variation and abnormal cellular activity in the SCZ brain. Herein, we provide a brief overview of some of the key omics methods currently being used in SCZ research, including RNA-seq, the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and high-throughput chromosome conformation capture (3C) approaches (e.g., Hi-C), as well as single-cell/nuclei iterations of these methods. We also discuss how these techniques are being employed to further our understanding of the genetic basis of SCZ, and to identify associated molecular pathways, biomarkers, and candidate drug targets.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Activity-Dependent Transcriptional Program in NGN2+ Neurons Enriched for Genetic Risk for Brain-Related Disorders

BACKGROUND

Converging evidence from large-scale genetic and postmortem studies highlights the role of aberrant neurotransmission and genetic regulation in brain-related disorders. However, identifying neuronal activity-regulated transcriptional programs in the human brain and understanding how changes contribute to disease remain challenging.

METHODS

To better understand how the activity-dependent regulome contributes to risk for brain-related disorders, we profiled the transcriptomic and epigenomic changes following neuronal depolarization in human induced pluripotent stem cell-derived glutamatergic neurons (NGN2) from 6 patients with schizophrenia and 5 control participants.

RESULTS

Multiomic data integration associated global patterns of chromatin accessibility with gene expression and identified enhancer-promoter interactions in glutamatergic neurons. Within 1 hour of potassium chloride-induced depolarization, independent of diagnosis, glutamatergic neurons displayed substantial activity-dependent changes in the expression of genes regulating synaptic function. Depolarization-induced changes in the regulome revealed significant heritability enrichment for schizophrenia and Parkinson's disease, adding to mounting evidence that sequence variation within activation-dependent regulatory elements contributes to the genetic risk for brain-related disorders. Gene coexpression network analysis elucidated interactions among activity-dependent and disease-associated genes and pointed to a key driver (NAV3) that interacted with multiple genes involved in axon guidance.

CONCLUSIONS

Overall, we demonstrated that deciphering the activity-dependent regulome in glutamatergic neurons reveals novel targets for advanced diagnosis and therapy.

Convergence of the dysregulated regulome in schizophrenia with polygenic risk and evolutionarily constrained enhancers

Enhancers play an essential role in the etiology of schizophrenia; however, the dysregulation of enhancer activity and its impact on the regulome in schizophrenia remains understudied. To address this gap in our knowledge, we assessed enhancer and gene expression in 1,382 brain samples comprising cases with schizophrenia and unaffected controls. Dysregulation of enhancer expression was concordant with changes in gene expression, and was more closely associated with schizophrenia polygenic risk, suggesting that enhancer dysregulation is proximal to the genetic etiology of the disease. Modeling the shared variance of cis-coordinated genes and enhancers revealed a gene regulatory program that was highly associated with genetic vulnerability to schizophrenia. By integrating coordinated factors with evolutionary constraints, we found that enhancers acquired during human evolution are more likely to regulate genes that are implicated in neuropsychiatric disorders and, thus, hold potential as therapeutic targets. Our analysis provides a systematic view of regulome dysregulation in schizophrenia and highlights its convergence with schizophrenia polygenic risk and human-gained enhancers.

Transcriptional correlates of cocaine-associated learning in striatal ARC ensembles

Learned associations between the rewarding effects of drugs and the context in which they are experienced underlie context-induced relapse. Previous work demonstrates the importance of sparse neuronal populations - called neuronal ensembles - in associative learning and cocaine seeking, but it remains unknown whether the encoding vs. retrieval of cocaine-associated memories involves similar or distinct mechanisms of ensemble activation and reactivation in nucleus accumbens (NAc). We use ArcCreER T2 mice to establish that mostly distinct NAc ensembles are recruited by initial vs. repeated exposures to cocaine, which are then differentially reactivated and exert distinct effects during cocaine-related memory retrieval. Single-nuclei RNA-sequencing of these ensembles demonstrates predominant recruitment of D1 medium spiny neurons and identifies transcriptional properties that are selective to cocaine-recruited NAc neurons and could explain distinct excitability features. These findings fundamentally advance our understanding of how cocaine drives pathological memory formation during repeated exposures.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

DeepGAMI: deep biologically guided auxiliary learning for multimodal integration and imputation to improve genotype-phenotype prediction

BACKGROUND

Genotypes are strongly associated with disease phenotypes, particularly in brain disorders. However, the molecular and cellular mechanisms behind this association remain elusive. With emerging multimodal data for these mechanisms, machine learning methods can be applied for phenotype prediction at different scales, but due to the black-box nature of machine learning, integrating these modalities and interpreting biological mechanisms can be challenging. Additionally, the partial availability of these multimodal data presents a challenge in developing these predictive models.

METHOD

To address these challenges, we developed DeepGAMI, an interpretable neural network model to improve genotype-phenotype prediction from multimodal data. DeepGAMI leverages functional genomic information, such as eQTLs and gene regulation, to guide neural network connections. Additionally, it includes an auxiliary learning layer for cross-modal imputation allowing the imputation of latent features of missing modalities and thus predicting phenotypes from a single modality. Finally, DeepGAMI uses integrated gradient to prioritize multimodal features for various phenotypes.

RESULTS

We applied DeepGAMI to several multimodal datasets including genotype and bulk and cell-type gene expression data in brain diseases, and gene expression and electrophysiology data of mouse neuronal cells. Using cross-validation and independent validation, DeepGAMI outperformed existing methods for classifying disease types, and cellular and clinical phenotypes, even using single modalities (e.g., AUC score of 0.79 for Schizophrenia and 0.73 for cognitive impairment in Alzheimer's disease).

CONCLUSION

We demonstrated that DeepGAMI improves phenotype prediction and prioritizes phenotypic features and networks in multiple multimodal datasets in complex brains and brain diseases. Also, it prioritized disease-associated variants, genes, and regulatory networks linked to different phenotypes, providing novel insights into the interpretation of gene regulatory mechanisms. DeepGAMI is open-source and available for general use.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

Multi-omic atlas of the parahippocampal gyrus in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia worldwide, with a projection of 151 million cases by 2050. Previous genetic studies have identified three main genes associated with early-onset familial Alzheimer's disease, however this subtype accounts for less than 5% of total cases. Next-generation sequencing has been well established and holds great promise to assist in the development of novel therapeutics as well as biomarkers to prevent or slow the progression of this devastating disease. Here we present a public resource of functional genomic data from the parahippocampal gyrus of 201 postmortem control, mild cognitively impaired (MCI) and AD individuals from the Mount Sinai brain bank, of which whole-genome sequencing (WGS), and bulk RNA sequencing (RNA-seq) were previously published. The genomic data include bulk proteomics and DNA methylation, as well as cell-type-specific RNA-seq and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) data. We have performed extensive preprocessing and quality control, allowing the research community to access and utilize this public resource available on the Synapse platform at https://doi.org/10.7303/syn51180043.2 .

The stability of the myelinating oligodendrocyte transcriptome is regulated by the nuclear lamina

Oligodendrocytes are specialized cells that insulate and support axons with their myelin membrane, allowing proper brain function. Here, we identify lamin A/C (LMNA/C) as essential for transcriptional and functional stability of myelinating oligodendrocytes. We show that LMNA/C levels increase with differentiation of progenitors and that loss of Lmna in differentiated oligodendrocytes profoundly alters their chromatin accessibility and transcriptional signature. Lmna deletion in myelinating glia is compatible with normal developmental myelination. However, altered chromatin accessibility is detected in fully differentiated oligodendrocytes together with increased expression of progenitor genes and decreased levels of lipid-related transcription factors and inner mitochondrial membrane transcripts. These changes are accompanied by altered brain metabolism, lower levels of myelin-related lipids, and altered mitochondrial structure in oligodendrocytes, thereby resulting in myelin thinning and the development of a progressively worsening motor phenotype. Overall, our data identify LMNA/C as essential for maintaining the transcriptional and functional stability of myelinating oligodendrocytes.

Genome-wide methylomic regulation of multiscale gene networks in Alzheimer's disease

INTRODUCTION

Recent studies revealed the association of abnormal methylomic changes with Alzheimer's disease (AD) but there is a lack of systematic study of the impact of methylomic alterations over the molecular networks underlying AD.

METHODS

We profiled genome-wide methylomic variations in the parahippocampal gyrus from 201 post mortem control, mild cognitive impaired, and AD brains.

RESULTS

We identified 270 distinct differentially methylated regions (DMRs) associated with AD. We quantified the impact of these DMRs on each gene and each protein as well as gene and protein co-expression networks. DNA methylation had a profound impact on both AD-associated gene/protein modules and their key regulators. We further integrated the matched multi-omics data to show the impact of DNA methylation on chromatin accessibility, which further modulates gene and protein expression.

DISCUSSION

The quantified impact of DNA methylation on gene and protein networks underlying AD identified potential upstream epigenetic regulators of AD.

HIGHLIGHTS

A cohort of DNA methylation data in the parahippocampal gyrus was developed from 201 post mortem control, mild cognitive impaired, and Alzheimer's disease (AD) brains. Two hundred seventy distinct differentially methylated regions (DMRs) were found to be associated with AD compared to normal control. A metric was developed to quantify methylation impact on each gene and each protein. DNA methylation was found to have a profound impact on not only the AD-associated gene modules but also key regulators of the gene and protein networks. Key findings were validated in an independent multi-omics cohort in AD. The impact of DNA methylation on chromatin accessibility was also investigated by integrating the matched methylomic, epigenomic, transcriptomic, and proteomic data.

Dopaminoceptive D1 and D2 neurons in ventral hippocampus arbitrate approach and avoidance in anxiety

The hippocampus 1-7, as well as dopamine circuits 8-11, coordinate decision-making in anxiety-eliciting situations. Yet, little is known about how dopamine modulates hippocampal representations of emotionally-salient stimuli to inform appropriate resolution of approach versus avoidance conflicts. We here study dopaminoceptive neurons in mouse ventral hippocampus (vHipp), molecularly distinguished by their expression of dopamine D1 or D2 receptors. We show that these neurons are transcriptionally distinct and topographically organized across vHipp subfields and cell types. In the ventral subiculum where they are enriched, both D1 and D2 neurons are recruited during anxiogenic exploration, yet with distinct profiles related to investigation and behavioral selection. In turn, they mediate opposite approach/avoidance responses, and are differentially modulated by dopaminergic transmission in that region. Together, these results suggest that vHipp dopamine dynamics gate exploratory behaviors under contextual uncertainty, implicating dopaminoception in the complex computation engaged in vHipp to govern emotional states.

Interleukin-3 coordinates glial-peripheral immune crosstalk to incite multiple sclerosis

Glial cells and central nervous system (CNS)-infiltrating leukocytes contribute to multiple sclerosis (MS). However, the networks that govern crosstalk among these ontologically distinct populations remain unclear. Here, we show that, in mice and humans, CNS-resident astrocytes and infiltrating CD44hiCD4+ T cells generated interleukin-3 (IL-3), while microglia and recruited myeloid cells expressed interleukin-3 receptor-ɑ (IL-3Rɑ). Astrocytic and T cell IL-3 elicited an immune migratory and chemotactic program by IL-3Rɑ+ myeloid cells that enhanced CNS immune cell infiltration, exacerbating MS and its preclinical model. Multiregional snRNA-seq of human CNS tissue revealed the appearance of IL3RA-expressing myeloid cells with chemotactic programming in MS plaques. IL3RA expression by plaque myeloid cells and IL-3 amount in the cerebrospinal fluid predicted myeloid and T cell abundance in the CNS and correlated with MS severity. Our findings establish IL-3:IL-3RA as a glial-peripheral immune network that prompts immune cell recruitment to the CNS and worsens MS.

ZBTB7A regulates MDD-specific chromatin signatures and astrocyte-mediated stress vulnerability in orbitofrontal cortex

Hyperexcitability in the orbitofrontal cortex (OFC) is a key clinical feature of anhedonic domains of Major Depressive Disorder (MDD). However, the cellular and molecular substrates underlying this dysfunction remain unknown. Here, cell-population-specific chromatin accessibility profiling in human OFC unexpectedly mapped genetic risk for MDD exclusively to non-neuronal cells, and transcriptomic analyses revealed significant glial dysregulation in this region. Characterization of MDD-specific cis-regulatory elements identified ZBTB7A - a transcriptional regulator of astrocyte reactivity - as an important mediator of MDD-specific chromatin accessibility and gene expression. Genetic manipulations in mouse OFC demonstrated that astrocytic Zbtb7a is both necessary and sufficient to promote behavioral deficits, cell-type-specific transcriptional and chromatin profiles, and OFC neuronal hyperexcitability induced by chronic stress - a major risk factor for MDD. These data thus highlight a critical role for OFC astrocytes in stress vulnerability and pinpoint ZBTB7A as a key dysregulated factor in MDD that mediates maladaptive astrocytic functions driving OFC hyperexcitability.

Efficient differential expression analysis of large-scale single cell transcriptomics data using dreamlet

Advances in single-cell and -nucleus transcriptomics have enabled generation of increasingly large-scale datasets from hundreds of subjects and millions of cells. These studies promise to give unprecedented insight into the cell type specific biology of human disease. Yet performing differential expression analyses across subjects remains difficult due to challenges in statistical modeling of these complex studies and scaling analyses to large datasets. Our open-source R package dreamlet (DiseaseNeurogenomics.github.io/dreamlet) uses a pseudobulk approach based on precision-weighted linear mixed models to identify genes differentially expressed with traits across subjects for each cell cluster. Designed for data from large cohorts, dreamlet is substantially faster and uses less memory than existing workflows, while supporting complex statistical models and controlling the false positive rate. We demonstrate computational and statistical performance on published datasets, and a novel dataset of 1.4M single nuclei from postmortem brains of 150 Alzheimer's disease cases and 149 controls.

Efficient differential expression analysis of large-scale single cell transcriptomics data using dreamlet

Advances in single-cell and -nucleus transcriptomics have enabled generation of increasingly large-scale datasets from hundreds of subjects and millions of cells. These studies promise to give unprecedented insight into the cell type specific biology of human disease. Yet performing differential expression analyses across subjects remains difficult due to challenges in statistical modeling of these complex studies and scaling analyses to large datasets. Our open-source R package dreamlet (DiseaseNeurogenomics.github.io/dreamlet) uses a pseudobulk approach based on precision-weighted linear mixed models to identify genes differentially expressed with traits across subjects for each cell cluster. Designed for data from large cohorts, dreamlet is substantially faster and uses less memory than existing workflows, while supporting complex statistical models and controlling the false positive rate. We demonstrate computational and statistical performance on published datasets, and a novel dataset of 1.4M single nuclei from postmortem brains of 150 Alzheimer's disease cases and 149 controls.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods have been developed to estimate cell-type proportions and gene expressions in bulk-tissue samples. However, the performance of these methods and their biological applications has not been evaluated, particularly on human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk-tissue RNAseq, single-cell/nuclei (sc/sn) RNAseq, and immunohistochemistry. A total of 1,130,767 nuclei/cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expression. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk-tissue or single-cell eQTLs alone. Differential gene expression associated with multiple phenotypes were also examined using the deconvoluted data. Our findings, which were replicated in bulk-tissue RNAseq and sc/snRNAseq data, provided new insights into the biological applications of deconvoluted data.

Genetic regulation of cell-type specific chromatin accessibility shapes the etiology of brain diseases

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are major risk factors of human disease. We measured chromatin accessibility in sorted neurons and glia from 1,932 samples of human postmortem brain and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTL). Only 10.4% of caQTL are shared between neurons and glia, supporting the cell type specificity of genetic regulation of the brain regulome. Incorporating allele specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms underlying disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs, identifying the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides novel insights into brain disease etiology.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states.

Impact of schizophrenia GWAS loci converge onto distinct pathways in cortical interneurons vs glutamatergic neurons during development

Remarkable advances have been made in schizophrenia (SCZ) GWAS, but gleaning biological insight from these loci is challenging. Genetic influences on gene expression (e.g., eQTLs) are cell type-specific, but most studies that attempt to clarify GWAS loci's influence on gene expression have employed tissues with mixed cell compositions that can obscure cell-specific effects. Furthermore, enriched SCZ heritability in the fetal brain underscores the need to study the impact of SCZ risk loci in specific developing neurons. MGE-derived cortical interneurons (cINs) are consistently affected in SCZ brains and show enriched SCZ heritability in human fetal brains. We identified SCZ GWAS risk genes that are dysregulated in iPSC-derived homogeneous populations of developing SCZ cINs. These SCZ GWAS loci differential expression (DE) genes converge on the PKC pathway. Their disruption results in PKC hyperactivity in developing cINs, leading to arborization deficits. We show that the fine-mapped GWAS locus in the ATP2A2 gene of the PKC pathway harbors enhancer marks by ATACseq and ChIPseq, and regulates ATP2A2 expression. We also generated developing glutamatergic neurons (GNs), another population with enriched SCZ heritability, and confirmed their functionality after transplantation into the mouse brain. Then, we identified SCZ GWAS risk genes that are dysregulated in developing SCZ GNs. GN-specific SCZ GWAS loci DE genes converge on the ion transporter pathway, distinct from those for cINs. Disruption of the pathway gene CACNA1D resulted in deficits of Ca2+ currents in developing GNs, suggesting compromised neuronal function by GWAS loci pathway deficits during development. This study allows us to identify cell type-specific and developmental stage-specific mechanisms of SCZ risk gene function, and may aid in identifying mechanism-based novel therapeutic targets.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis.

The Neuroepigenome: Implications of Chemical and Physical Modifications of Genomic DNA in Schizophrenia

Schizophrenia is a chronic mental illness with a substantial genetic component. To unfold the complex etiology of schizophrenia, it is important to understand the interplay between genetic and nongenetic factors. Genetic factors involve variation in the DNA sequences of protein-coding genes, which directly contribute to phenotypic traits, and variation in noncoding sequences, which comprise 98% of the genome and contain DNA elements known to play a role in regulating gene expression. The epigenome refers to the chemical modifications on both DNA and the structural proteins that package DNA into the nucleus, which together regulate gene expression in specific cell types, conditions, and developmental stages. The dynamic nature of the epigenome makes it an ideal tool to investigate the relationship between inherited genetic mutations associated with schizophrenia and altered gene regulation throughout the course of brain development. In this review, we focus on the current understanding of the role of epigenetic marks and their three-dimensional nuclear organization in the developmental trajectory of distinct brain cell types to decipher the complex gene regulatory mechanisms that are disrupted in schizophrenia.

A translational genomics approach identifies IL10RB as the top candidate gene target for COVID-19 susceptibility

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes that reduce COVID-19 host susceptibility is a critical next step. Using a translational genomics approach that integrates COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX), and perturbagen signatures, we identified IL10RB as the top candidate gene target for COVID-19 host susceptibility. In a series of validation steps, we show that predicted GReX upregulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes and that in vitro IL10RB overexpression is associated with increased viral load and activation of disease-relevant molecular pathways.

Genetics of the human microglia regulome refines Alzheimer's disease risk loci

Microglia are brain myeloid cells that play a critical role in neuroimmunity and the etiology of Alzheimer's disease (AD), yet our understanding of how the genetic regulatory landscape controls microglial function and contributes to AD is limited. Here, we performed transcriptome and chromatin accessibility profiling in primary human microglia from 150 donors to identify genetically driven variation and cell-specific enhancer-promoter (E-P) interactions. Integrative fine-mapping analysis identified putative regulatory mechanisms for 21 AD risk loci, of which 18 were refined to a single gene, including 3 new candidate risk genes (KCNN4, FIBP and LRRC25). Transcription factor regulatory networks captured AD risk variation and identified SPI1 as a key putative regulator of microglia expression and AD risk. This comprehensive resource capturing variation in the human microglia regulome provides insights into the etiology of neurodegenerative disease.

Common variants contribute to intrinsic human brain functional networks

The human brain forms functional networks of correlated activity, which have been linked with both cognitive and clinical outcomes. However, the genetic variants affecting brain function are largely unknown. Here, we used resting-state functional magnetic resonance images from 47,276 individuals to discover and validate common genetic variants influencing intrinsic brain activity. We identified 45 new genetic regions associated with brain functional signatures (P < 2.8 × 10^-11), including associations to the central executive, default mode, and salience networks involved in the triple-network model of psychopathology. A number of brain activity-associated loci colocalized with brain disorders (e.g., the APOE ε4 locus with Alzheimer's disease). Variation in brain function was genetically correlated with brain disorders, such as major depressive disorder and schizophrenia. Together, our study provides a step forward in understanding the genetic architecture of brain functional networks and their genetic links to brain-related complex traits and disorders.

A complete temporal transcription factor series in the fly visual system

The brain consists of thousands of neuronal types that are generated by stem cells producing different neuronal types as they age. In Drosophila, this temporal patterning is driven by the successive expression of temporal transcription factors (tTFs)^1-6. Here we used single-cell mRNA sequencing to identify the complete series of tTFs that specify most Drosophila optic lobe neurons. We verify that tTFs regulate the progression of the series by activating the next tTF(s) and repressing the previous one(s), and also identify more complex mechanisms of regulation. Moreover, we establish the temporal window of origin and birth order of each neuronal type in the medulla and provide evidence that these tTFs are sufficient to explain the generation of all of the neuronal diversity in this brain region. Finally, we describe the first steps of neuronal differentiation and show that these steps are conserved in humans. We find that terminal differentiation genes, such as neurotransmitter-related genes, are present as transcripts, but not as proteins, in immature larval neurons. This comprehensive analysis of a temporal series of tTFs in the optic lobe offers mechanistic insights into how tTF series are regulated, and how they can lead to the generation of a complete set of neurons.

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 nucleosome to megabase-ranging domains encompassing multiple transcriptional units including heritability loci for psychiatric traits, remains largely unexplored in the human brain. Here, we construct promoter and enhancer enriched nucleosomal histone modification landscapes for adult prefrontal cortex (PFC) from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from (n=739) 388 controls and 351 subjects diagnosed with schizophrenia (SCZ) or bipolar disorder (BD). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 10⁴-10⁶ bp chromosomal organization, firmly integrated into Hi-C topologically associating domain (TAD) stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyperacetylated 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 kilo- to megabase-scaling chromosomal domains.

Non-cell-autonomous disruption of nuclear architecture as a potential cause of COVID-19-induced anosmia

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.

A bidirectional competitive interaction between circHomer1 and Homer1b within the orbitofrontal cortex regulates reversal learning

Although circular RNAs (circRNAs) are enriched in the brain, their relevance for brain function and psychiatric disorders is poorly understood. Here, we show that circHomer1 is inversely associated with relative HOMER1B mRNA isoform levels in both the orbitofrontal cortex (OFC) and stem-cell-derived neuronal cultures of subjects with psychiatric disorders. We further demonstrate that in vivo circHomer1 knockdown (KD) within the OFC can inhibit the synaptic expression of Homer1b mRNA. Furthermore, we show that circHomer1 directly binds to Homer1b mRNA and that Homer1b-specific KD increases synaptic circHomer1 levels and improves OFC-mediated behavioral flexibility. Importantly, double circHomer1 and Homer1b in vivo co-KD results in a complete rescue in circHomer1-associated alterations in both chance reversal learning and synaptic gene expression. Lastly, we uncover an RNA-binding protein that can directly bind to circHomer1 and promote its biogenesis. Taken together, our data provide mechanistic insights into the importance of circRNAs in brain function and disease.

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.

Single-nucleus transcriptome analysis of human brain immune response in patients with severe COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has been associated with neurological and neuropsychiatric illness in many individuals. We sought to further our understanding of the relationship between brain tropism, neuro-inflammation, and host immune response in acute COVID-19 cases.

METHODS

Three brain regions (dorsolateral prefrontal cortex, medulla oblongata, and choroid plexus) from 5 patients with severe COVID-19 and 4 controls were examined. The presence of the virus was assessed by western blot against viral spike protein, as well as viral transcriptome analysis covering > 99% of SARS-CoV-2 genome and all potential serotypes. Droplet-based single-nucleus RNA sequencing (snRNA-seq) was performed in the same samples to examine the impact of COVID-19 on transcription in individual cells of the brain.

RESULTS

Quantification of viral spike S1 protein and viral transcripts did not detect SARS-CoV-2 in the postmortem brain tissue. However, analysis of 68,557 single-nucleus transcriptomes from three distinct regions of the brain identified an increased proportion of stromal cells, monocytes, and macrophages in the choroid plexus of COVID-19 patients. Furthermore, differential gene expression, pseudo-temporal trajectory, and gene regulatory network analyses revealed transcriptional changes in the cortical microglia associated with a range of biological processes, including cellular activation, mobility, and phagocytosis.

CONCLUSIONS

Despite the absence of detectable SARS-CoV-2 in the brain at the time of death, the findings suggest significant and persistent neuroinflammation in patients with acute COVID-19.

Common genetic variation influencing human white matter microstructure

Brain regions communicate with each other through tracts of myelinated axons, commonly referred to as white matter. We identified common genetic variants influencing white matter microstructure using diffusion magnetic resonance imaging of 43,802 individuals. Genome-wide association analysis identified 109 associated loci, 30 of which were detected by tract-specific functional principal components analysis. A number of loci colocalized with brain diseases, such as glioma and stroke. Genetic correlations were observed between white matter microstructure and 57 complex traits and diseases. Common variants associated with white matter microstructure altered the function of regulatory elements in glial cells, particularly oligodendrocytes. This large-scale tract-specific study advances the understanding of the genetic architecture of white matter and its genetic links to a wide spectrum of clinical outcomes.

IL10RB as a key regulator of COVID-19 host susceptibility and severity

BACKGROUND

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes and readily available compounds that reduce COVID-19 host susceptibility is a critical next step.

METHODS

We integrate COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX) and perturbargen signatures to identify candidate genes and compounds that reverse the predicted gene expression dysregulation associated with COVID-19 susceptibility. The top candidate gene is validated by testing both its GReX and observed blood transcriptome association with COVID-19 severity, as well as by in vitro perturbation to quantify effects on viral load and molecular pathway dysregulation. We validate the in silico drug repositioning analysis by examining whether the top candidate compounds decrease COVID-19 incidence based on epidemiological evidence.

RESULTS

We identify IL10RB as the top key regulator of COVID-19 host susceptibility. Predicted GReX up-regulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes. In vitro IL10RB overexpression is associated with increased viral load and activation of immune-related molecular pathways. Azathioprine and retinol are prioritized as candidate compounds to reduce the likelihood of testing positive for COVID-19.

CONCLUSIONS

We establish an integrative data-driven approach for gene target prioritization. We identify and validate IL10RB as a suitable molecular target for modulation of COVID-19 host susceptibility. Finally, we provide evidence for a few readily available medications that would warrant further investigation as drug repositioning candidates.

Integration of Alzheimer's disease genetics and myeloid genomics identifies disease risk regulatory elements and genes

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.

Disruption of nuclear architecture as a cause of COVID-19 induced anosmia

Olfaction relies on a coordinated partnership between odorant flow and neuronal communication. Disruption in our ability to detect odors, or anosmia, has emerged as a hallmark symptom of infection with SARS-CoV-2, yet the mechanism behind this abrupt sensory deficit remains elusive. Here, using molecular evaluation of human olfactory epithelium (OE) from subjects succumbing to COVID-19 and a hamster model of SARS-CoV-2 infection, we discovered widespread downregulation of olfactory receptors (ORs) as well as key components of their signaling pathway. OR downregulation likely represents a non-cell autonomous effect, since SARS-CoV-2 detection in OSNs is extremely rare both in human and hamster OEs. A likely explanation for the reduction of OR transcription is the striking reorganization of nuclear architecture observed in the OSN lineage, which disrupts multi-chromosomal compartments regulating OR expression in humans and hamsters. Our experiments uncover a novel molecular mechanism by which a virus with a very selective tropism can elicit persistent transcriptional changes in cells that evade it, contributing to the severity of COVID-19.

Dissecting Hes-centred transcriptional networks in neural stem cell maintenance and tumorigenesis in Drosophila

Neural stem cells divide during embryogenesis and juvenile life to generate the entire complement of neurons and glia in the nervous system of vertebrates and invertebrates. Studies of the mechanisms controlling the fine balance between neural stem cells and more differentiated progenitors have shown that, in every asymmetric cell division, progenitors send a Delta-Notch signal to their sibling stem cells. Here, we show that excessive activation of Notch or overexpression of its direct targets of the Hes family causes stem-cell hyperplasias in the Drosophila larval central nervous system, which can progress to malignant tumours after allografting to adult hosts. We combined transcriptomic data from these hyperplasias with chromatin occupancy data for Dpn, a Hes transcription factor, to identify genes regulated by Hes factors in this process. We show that the Notch/Hes axis represses a cohort of transcription factor genes. These are excluded from the stem cells and promote early differentiation steps, most likely by preventing the reversion of immature progenitors to a stem-cell fate. We describe the impact of two of these 'anti-stemness' factors, Zfh1 and Gcm, on Notch/Hes-triggered tumorigenesis.

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.

Assessment of somatic single-nucleotide variation in brain tissue of cases with schizophrenia

The genetic architecture of schizophrenia (SCZ) includes numerous risk loci across a range of frequencies and sizes, including common and rare single-nucleotide variants and insertions/deletions (indels), as well as rare copy number variants (CNVs). Despite the clear heritability of the disease, monozygotic twins are discordant for SCZ at a significant rate. Somatic variants-genetic changes that arise after fertilization rather than through germline inheritance-are widespread in the human brain and known to contribute to risk for both rare and common neuropsychiatric conditions. The contribution of somatic variants in the brain to risk of SCZ remains to be determined. In this study, we surveyed somatic single-nucleotide variants (sSNVs) in the brains of controls and individuals with SCZ (n = 10 and n = 9, respectively). From each individual, whole-exome sequencing (WES) was performed on DNA from neuronal and non-neuronal nuclei isolated by fluorescence activated nuclear sorting (FANS) from frozen postmortem prefrontal cortex (PFC) samples, as well as DNA extracted from temporal muscle as a reference. We identified an increased burden of sSNVs in cases compared to controls (SCZ rate = 2.78, control rate = 0.70; P = 0.0092, linear mixed effects model), that included a higher rate of non-synonymous and loss-of-function variants (SCZ rate = 1.33, control rate = 0.50; P = 0.047, linear mixed effects model). Our findings suggest sSNVs in the brain may constitute an additional component of the complex genetic architecture of SCZ. This perspective argues for the need to further investigate somatic variation in the brain as an explanation of the discordance in monozygotic twins and a potential guide to the identification of novel therapeutic targets.

The Mount Sinai cohort of large-scale genomic, transcriptomic and proteomic data in Alzheimer's disease

Alzheimer's disease (AD) affects half the US population over the age of 85 and is universally fatal following an average course of 10 years of progressive cognitive disability. Genetic and genome-wide association studies (GWAS) have identified about 33 risk factor genes for common, late-onset AD (LOAD), but these risk loci fail to account for the majority of affected cases and can neither provide clinically meaningful prediction of development of AD nor offer actionable mechanisms. This cohort study generated large-scale matched multi-Omics data in AD and control brains for exploring novel molecular underpinnings of AD. Specifically, we generated whole genome sequencing, whole exome sequencing, transcriptome sequencing and proteome profiling data from multiple regions of 364 postmortem control, mild cognitive impaired (MCI) and AD brains with rich clinical and pathophysiological data. All the data went through rigorous quality control. Both the raw and processed data are publicly available through the Synapse software platform.

Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia

Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.

An atlas of chromatin accessibility in the adult human brain

Most common genetic risk variants associated with neuropsychiatric disease are noncoding and are thought to exert their effects by disrupting the function of cis 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.

Open chromatin profiling of human postmortem brain infers functional roles for non-coding schizophrenia loci

Open chromatin provides access to DNA-binding proteins for the correct spatiotemporal regulation of gene expression. Mapping chromatin accessibility has been widely used to identify the location of cis regulatory elements (CREs) including promoters and enhancers. CREs show tissue- and cell-type specificity and disease-associated variants are often enriched for CREs in the tissues and cells that pertain to a given disease. To better understand the role of CREs in neuropsychiatric disorders we applied the Assay for Transposase Accessible Chromatin followed by sequencing (ATAC-seq) to neuronal and non-neuronal nuclei isolated from frozen postmortem human brain by fluorescence-activated nuclear sorting (FANS). Most of the identified open chromatin regions (OCRs) are differentially accessible between neurons and non-neurons, and show enrichment with known cell type markers, promoters and enhancers. Relative to those of non-neurons, neuronal OCRs are more evolutionarily conserved and are enriched in distal regulatory elements. Transcription factor (TF) footprinting analysis identifies differences in the regulome between neuronal and non-neuronal cells and ascribes putative functional roles to a number of non-coding schizophrenia (SCZ) risk variants. Among the identified variants is a Single Nucleotide Polymorphism (SNP) proximal to the gene encoding SNX19. In vitro experiments reveal that this SNP leads to an increase in transcriptional activity. As elevated expression of SNX19 has been associated with SCZ, our data provide evidence that the identified SNP contributes to disease. These results represent the first analysis of OCRs and TF-binding sites in distinct populations of postmortem human brain cells and further our understanding of the regulome and the impact of neuropsychiatric disease-associated genetic risk variants.

Striatal H3K27 Acetylation Linked to Glutamatergic Gene Dysregulation in Human Heroin Abusers Holds Promise as Therapeutic Target

BACKGROUND

Opiate abuse and overdose reached epidemic levels in the United States. However, despite significant advances in animal and in vitro models, little knowledge has been directly accrued regarding the neurobiology of the opiate-addicted human brain.

METHODS

We used postmortem human brain specimens from a homogeneous European Caucasian population of heroin users for transcriptional and epigenetic profiling, as well as direct assessment of chromatin accessibility in the striatum, a brain region central to reward and emotion. A rat heroin self-administration model was used to obtain translational molecular and behavioral insights.

RESULTS

Our transcriptome approach revealed marked impairments related to glutamatergic neurotransmission and chromatin remodeling in the human striatum. A series of biochemical experiments tracked the specific location of the epigenetic disturbances to hyperacetylation of lysine 27 of histone H3, showing dynamic correlations with heroin use history and acute opiate toxicology. Targeted investigation of GRIA1, a glutamatergic gene implicated in drug-seeking behavior, verified the increased enrichment of lysine-27 acetylated histone H3 at discrete loci, accompanied by enhanced chromatin accessibility at hyperacetylated regions in the gene body. Analogous epigenetic impairments were detected in the striatum of heroin self-administering rats. Using this translational model, we showed that bromodomain inhibitor JQ1, which blocks the functional readout of acetylated lysines, reduced heroin self-administration and cue-induced drug-seeking behavior.

CONCLUSIONS

Overall, our data suggest that heroin-related histone H3 hyperacetylation contributes to glutamatergic transcriptional changes that underlie addiction behavior and identify JQ1 as a promising candidate for targeted clinical interventions in heroin use disorder.

Understanding the genetic liability to schizophrenia through the neuroepigenome

The Psychiatric Genomics Consortium-Schizophrenia Workgroup (PGC-SCZ) recently identified 108 loci associated with increased risk for schizophrenia (SCZ). The vast majority of these variants reside within non-coding sequences of the genome and are predicted to exert their effects by affecting the mechanism of action of cis regulatory elements (CREs), such as promoters and enhancers. Although a number of large-scale collaborative efforts (e.g. ENCODE) have achieved a comprehensive mapping of CREs in human cell lines or tissue homogenates, it is becoming increasingly evident that many risk-associated variants are enriched for expression Quantitative Trait Loci (eQTLs) and CREs in specific tissues or cells. As such, data derived from previous research endeavors may not capture fully cell-type and/or region specific changes associated with brain diseases. Coupling recent technological advances in genomics with cell-type specific methodologies, we are presented with an unprecedented opportunity to better understand the genetics of normal brain development and function and, in turn, the molecular basis of neuropsychiatric disorders. In this review, we will outline ongoing efforts towards this goal and will discuss approaches with the potential to shed light on the mechanism(s) of action of cell-type specific cis regulatory elements and their putative roles in disease, with particular emphasis on understanding the manner in which the epigenome and CREs influence the etiology of SCZ.

Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases

Genome-wide association studies (GWAS) have identified hundreds of cardiometabolic disease (CMD) risk loci. However, they contribute little to genetic variance, and most downstream gene-regulatory mechanisms are unknown. We genotyped and RNA-sequenced vascular and metabolic tissues from 600 coronary artery disease patients in the Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task study (STARNET). Gene expression traits associated with CMD risk single-nucleotide polymorphism (SNPs) identified by GWAS were more extensively found in STARNET than in tissue- and disease-unspecific gene-tissue expression studies, indicating sharing of downstream cis-/trans-gene regulation across tissues and CMDs. In contrast, the regulatory effects of other GWAS risk SNPs were tissue-specific; abdominal fat emerged as an important gene-regulatory site for blood lipids, such as for the low-density lipoprotein cholesterol and coronary artery disease risk gene PCSK9 STARNET provides insights into gene-regulatory mechanisms for CMD risk loci, facilitating their translation into opportunities for diagnosis, therapy, and prevention.