Genetic regulation of cell type-specific chromatin accessibility shapes brain disease etiology.

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are risk factors for human disease. We measured chromatin accessibility in 1932 aliquots of sorted neurons and non-neurons from 616 human postmortem brains and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTLs). Only 10.4% of caQTLs are shared between neurons and non-neurons, which supports cell type-specific genetic regulation of the brain regulome. Incorporating allele-specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms that underlie disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs and identified the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides insights into disease etiology.

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.

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.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome.

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

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.

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.

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.

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.

Hypothalamic Fkbp51 is induced by fasting, and elevated hypothalamic expression promotes obese phenotypes.

To discover hypothalamic genes that might play a role in regulating energy balance, we carried out a microarray screen for genes induced by a 48-h fast in male C57Bl/6J mouse hypothalamus. One such gene was Fkbp51 (FK506 binding protein 5; Locus NP_034350). The product of this gene is of interest because it blocks glucocorticoid action, suggesting that fasting-induced elevation of this gene in the hypothalamus may reduce glucocorticoid negative feedback, leading to elevated glucocorticoid levels, thus promoting obese phenotypes. Subsequent analysis demonstrated that a 48-h fast induces Fkbp51 in ventromedial, paraventricular, and arcuate hypothalamic nuclei of mice and rats. To assess if hypothalamic Fkbp51 promotes obesity, the gene was transferred to the hypothalamus via an adeno-associated virus vector. Within 2 wk following Fkbp51 overexpression, mice on a high-fat diet exhibited elevated body weight, without hyperphagia, relative to mice receiving the control mCherry vector. Body weight remained elevated for more than 8 wk and was associated with elevated corticosterone and impaired glucose tolerance. These studies suggest that elevated hypothalamic Fkbp51 promotes obese phenotypes.

Specific preservation of biosynthetic responses to insulin in adipose tissue may contribute to hyperleptinemia in insulin-resistant obese mice.

Obesity is characterized by whole-body insulin resistance, yet the expression of many insulin-stimulated genes, including leptin, is elevated in obesity. These observations suggest that insulin resistance may depend on tissue type and gene. To address this hypothesis, we examined the regulation of immediate-early gene expression in liver and adipose tissue after injection of insulin and glucose, in lean insulin-sensitive, and in A(y)/a obese insulin-sensitive and obese insulin-resistant mice. Expression of hepatic jun-B mRNA was robustly increased after insulin injection in lean insulin-sensitive a/a mice and insulin-sensitive A(y)/a mice. In contrast, induction of hepatic jun-B and c-fos gene expression by insulin was markedly attenuated in obese insulin-resistant mice. Surprisingly, induction of adipose jun-B and c-fos gene expression by insulin was markedly enhanced in obese insulin-resistant mice. Furthermore, the expressions of jun-B and leptin were also enhanced in insulin-resistant mice after injection of glucose. Leptin mRNA was positively correlated with blood glucose levels and jun-B mRNA in lean but not insulin-resistant mice. Multiple regression analysis indicated that the correlation between leptin mRNA and jun-B mRNA was significant even after removing the effect of blood glucose, but the correlation between leptin mRNA and glucose was no longer significant after removing the effect of jun-B mRNA. These data suggest that some impairments in biosynthetic responses to insulin are manifest primarily in the liver, leading to hyperinsulinemia and stimulating the expression of some adipose insulin-stimulated genes, including leptin. These studies demonstrate the utility of immediate-early gene expression in the analysis of biosynthetic mechanisms of insulin resistance.