Cited by 103Genomics and Chromatin DynamicsGenetic Associations and EpidemiologyEpigenetics and DNA Methylation

The DOI matches its Crossref record; 88 of the 88 references it cites were found in a registry or an index; 9 of the 78 sources its claims cite were checked for support.

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

Open access
Hu, Benxia, Won, Hyejung, Mah, Won, Park, Royce B., Kassim, Bibi, Spiess, Keeley, Kozlenkov, Alexey, Crowley, Cheynna A., Pochareddy, Sirisha, The PsychENCODE Consortium, Ashley-Koch, Allison E., Crawford, Gregory E., Garrett, Melanie E., Song, Lingyun, Safi, Alexias, Johnson, Graham D., Wray, Gregory A., Reddy, Timothy E., Goes, Fernando S., Zandi, Peter, Bryois, Julien, Jaffe, Andrew E., Price, Amanda J., Ivanov, Nikolay A., Collado-Torres, Leonardo, Hyde, Thomas M., Burke, Emily E., Kleiman, Joel E., Tao, Ran, Shin, Joo Heon, Girdhar, Kiran, Jiang, Yan, Kundakovic, Marija, Brown, Leanne, Wiseman, Jennifer R., Zharovsky, Elizabeth, Jacobov, Rivka, Devillers, Olivia, Flatow, Elie, Hoffman, Gabriel E., Belmont, Judson, Del Valle, Diane, Francoeur, Nancy, Hadjimichael, Evi, Pinto, Dalila, van Bakel, Harm, Roussos, Panos, Fullard, John F., Bendl, Jaroslav, Hauberg, Mads E., Charney, Alexander W., Haroutunian, Vahram, Lipska, Barbara K., Lewis, David A., Hahn, Chang-Gyu, Mangravite, Lara M., Peters, Mette A., Chae, Yooree, Peng, Junmin, Niu, Mingming, Wang, Xusheng, Webster, Maree J., Beach, Thomas G., Chen, Chao, Jiang, Yi, Dai, Rujia, Wang, Yongjun, Xia, Yan, Shieh, Annie W., Liu, Chunyu, Grennan, Kay S., Vadukapuram, Ramu, Giase, Gina, Fitzgerald, Dominic, Cheng, Lijun, Brown, Miguel, Brown, Mimi, Brunetti, Tonya, Goodman, Thomas, Alsayed, Majd, White, Kevin P., Ray, Mohana, Polioudakis, Damon, Wamsley, Brie, Yin, Jiani, De La Torre Ubieta, Luis, Gandal, Michael J., Swarup, Vivek, Sanders, Stephan J., State, Matthew W., Werling, Donna M., An, Joon-Yong, Sheppard, Brooke, Willsey, A. Jeremy, Kefi, Amira, Mattei, Eugenio, Purcaro, Michael, Weng, Zhiping, Moore, Jill, Pratt, Henry, Huey, Jack, Borrman, Tyler, Sullivan, Patrick F., Giusti-Rodriguez, Paola, Kim, Yunjung, Szatkiewicz, Jin, Rhie, Suhn Kyong, Armoskus, Christoper, Camarena, Adrian, Farnham, Peggy J., Spitsyna, Valeria N., Witt, Heather, Schreiner, Shannon, Evgrafov, Oleg V., Knowles, James A., Gerstein, Mark, Liu, Shuang, Navarro, Fabio C. P., Warrell, Jonathan, Clarke, Declan, Emani, Prashant S., Gu, Mengting, Shi, Xu, Xu, Min, Yang, Yucheng T., Kitchen, Robert R., Gürsoy, Gamze, Zhang, Jing, Carlyle, Becky C., Nairn, Angus C., Li, Mingfeng, Skarica, Mario, Li, Zhen, Sousa, Andre M. M., Santpere, Gabriel, Choi, Jinmyung, Zhu, Ying, Gao, Tianliuyun, Miller, Daniel J., Cherskov, Adriana, Yang, Mo, Amiri, Anahita, Coppola, Gianfilippo, Mariani, Jessica, Scuderi, Soraya, Szekely, Anna, Vaccarino, Flora M., Wu, Feinan, Weissman, Sherman, Wang, Daifeng, Roychowdhury, Tanmoy, Abyzov, Alexej, Li, Yun, Dracheva, Stella, Sestan, Nenad, Akbarian, Schahram, Geschwind, Daniel H.·2021·Nature Communications

Abstract

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.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

This paper integrates Hi–C chromosome conformation from sorted human neurons (NeuN+) and glia (NeuN−) with transcriptomic and enhancer data to map cell-type-specific 3D regulatory architecture in the brain. It uses these maps to refine the cellular etiology of Alzheimer’s disease, schizophrenia, and bipolar disorder.

Key findings

  • NeuN+ and NeuN− cells show extensive cell-type-specific chromatin architecture: 287 differential FIREs (145 NeuN+, 142 NeuN−), 253 NeuN+ super-FIREs and 157 NeuN− super-FIREs with only 9 shared, and 187,674 NeuN+ versus 167,551 NeuN− promoter-based interactions.
  • AD-associated epigenetic dysregulation maps to different cell classes by direction: AD hyperacetylated H3K27ac peaks were largely active in NeuN− cells and linked to 676 genes, while AD hypoacetylated peaks were largely active in NeuN+ cells and linked to 460 genes.
  • AD genetic risk converged on microglia: AD SNP-based heritability was enriched in NeuN− enhancers (NeuN− p = 8.41 × 10−9; NeuN+ p = 0.205), and NeuN− H-MAGMA identified 181 AD risk genes, enriched in microglial co-expression modules upregulated in AD.
  • SCZ and BD showed neuronal heritability enrichment and both mapped to Glu- and GABA-specific enhancers; risk genes shared parvalbumin-expressing interneuron (In6) enrichment, while BD showed upper-layer Ex1 specificity and SCZ showed Ex2 specificity.

Methods

  • Hi–C libraries were generated from FANS-sorted NeuN+ and NeuN− nuclei from four dorsolateral prefrontal cortex (DLPFC) samples; RNA-seq was also generated from sorted nuclei, and published H3K27ac data from Glu and GABA neurons were integrated.
  • Chromatin architecture was analyzed as A/B compartments, TADs, FIREs/super-FIREs, and promoter-anchored loops; cell-type-specific enhancer–promoter relationships were built by overlapping differential H3K27ac peaks with Hi–C interactions.
  • Disease analyses used LDSC for SNP-based heritability enrichment and H-MAGMA to map GWAS risk variants to genes using NeuN−, Glu, or GABA regulatory maps; GWAS datasets included AD, SCZ, and BD.

Limitations

  • The paper states that neuronal subtype-specific findings need verification with neuronal subtype-specific Hi–C data (e.g., Glu and GABA Hi–C) and that CRISPR/Cas9 enhancer targeting could further test risk-gene enhancer roles.
  • The paper states that AD epigenetic changes may reflect either cellular composition changes or changes in regulatory landscape/cellular function, requiring future investigation.
  • Evident from the text: Hi–C was based on only four DLPFC samples, and NeuN+/NeuN− fractions are broad cell classes rather than pure cell types; Glu/GABA regulatory maps were inferred by integrating subtype H3K27ac peaks with NeuN+ Hi–C rather than subtype Hi–C.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41467-021-24243-0

Cites

88 of 88 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Differential FIREs are associated with cell-type-specific gene regulation. a The number of FIREs detected in NeuN+ and NeuN− cells. b Differential FIREs were identified in NeuN+ and NeuN− cells. c, d Differential FIREs overlap with differential H3K27ac peaks in the corresponding cell types. A neuronal gene, GRIN2B, is located in NeuN+ specific FIREs (e), while two oligodendrocytic genes, OLIG1 and OLIG2, are located in NeuN− specific FIREs (f). FIREs and significance of FIRE scores in NeuN+ and NeuN− cells are depicted in green and purple, respectively. Boxplots in the right show expression levels of GRIN2B (FDR = 3.71e−12), OLIG1 (FDR = 7.34e−26) and OLIG2 (FDR = 8.06e−23) in NeuN+ (n = 4) and NeuN− (n = 4) cells. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = 1st to 3rd quartiles (Q); minima, Q1 − 1.5 × interquartile range (IQR); maxima, Q3 + 1.5 × IQR. *FDR < 0.05 calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. Gene ontology (GO) analysis for genes assigned to differential NeuN+ (g) and NeuN− (h) FIREs. The red line denotes FDR = 0.05. i Cellular expression levels of genes assigned to differential NeuN+ and NeuN− FIREs. Center, median; box = Q1−Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. Neurons, n = 131; astrocytes (Astro), n = 62; microglia (Micro), n = 16; endothelial (Endo), n = 20; oligodendrocytes (Oligo), n = 38. Source data are provided as a Source Data file. j Genes assigned to differential NeuN+ and NeuN− FIREs are enriched in neurons and glia, respectively. Ex excitatory neurons, In inhibitory neurons.
Fig. 1 Differential FIREs are associated with cell-type-specific gene regulation. a The number of FIREs detected in NeuN+ and NeuN− cells. b Differential FIREs were identified in NeuN+ and NeuN− cells. c, d Differential FIREs overlap with differential H3K27ac peaks in the corresponding cell types. A neuronal gene, GRIN2B, is located in NeuN+ specific FIREs (e), while two oligodendrocytic genes, OLIG1 and OLIG2, are located in NeuN− specific FIREs (f). FIREs and significance of FIRE scores in NeuN+ and NeuN− cells are depicted in green and purple, respectively. Boxplots in the right show expression levels of GRIN2B (FDR = 3.71e−12), OLIG1 (FDR = 7.34e−26) and OLIG2 (FDR = 8.06e−23) in NeuN+ (n = 4) and NeuN− (n = 4) cells. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = 1st to 3rd quartiles (Q); minima, Q1 − 1.5 × interquartile range (IQR); maxima, Q3 + 1.5 × IQR. *FDR < 0.05 calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. Gene ontology (GO) analysis for genes assigned to differential NeuN+ (g) and NeuN− (h) FIREs. The red line denotes FDR = 0.05. i Cellular expression levels of genes assigned to differential NeuN+ and NeuN− FIREs. Center, median; box = Q1−Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. Neurons, n = 131; astrocytes (Astro), n = 62; microglia (Micro), n = 16; endothelial (Endo), n = 20; oligodendrocytes (Oligo), n = 38. Source data are provided as a Source Data file. j Genes assigned to differential NeuN+ and NeuN− FIREs are enriched in neurons and glia, respectively. Ex excitatory neurons, In inhibitory neurons.CC BY · PubMed Central
Fig. 2 Enhancer–promoter interactions in NeuN+ and NeuN− cells. a (Left) cell-type-specific regulatory networks were built by linking genes to NeuN+ and NeuN− specific H3K27ac peaks via Hi–C interactions in NeuN+ and NeuN− cells, respectively. (Right) The number of cell-type-specific peaks and their assigned genes in NeuN+ and NeuN− cells is described. A neuronal gene, HOMER1, is engaged with NeuN+ specific H3K27ac peaks via loops in NeuN+ cells (b), while an oligodendrocyte gene, SOX10, is engaged with NeuN− specific H3K27ac peaks via loops in NeuN− cells (c). The regions that interact with the gene promoter (gray) are highlighted in green (NeuN+) and purple (NeuN−), respectively. Boxplots in the right show expression levels of HOMER1 (FDR = 7.26e−32) and SOX10 (FDR = 1.88e−49) in NeuN+ (n = 4) and NeuN− (n = 4) cells, respectively. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. *FDR < 0.05 calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. d Genes assigned to NeuN+ specific peaks are enriched for synaptic co-expression modules, while genes assigned to NeuN− specific peaks are enriched for co-expression modules involved in transcriptional regulation and immune response during neurodevelopment. Significant enrichment (Sig.), FDR < 0.05. Fisher’s exact test was used for statistics analysis. OR, odds ratio. e Genes assigned to NeuN+ specific peaks are more highly enriched for synaptic functions such as exocytosis, intracellular signal transduction, protein cluster and structural plasticity than genes assigned to NeuN− specific peaks. Sig., FDR < 0.05. Fisher’s exact test was used for statistics analysis. f Genes assigned to NeuN+ specific peaks are highly expressed in neurons, while genes assigned to NeuN− specific peaks are highly expressed in oligodendrocytes and astrocytes. Astro astrocytes, Micro microglia, Endo Endothelial, Oligo oligodendrocytes.
Fig. 2 Enhancer–promoter interactions in NeuN+ and NeuN− cells. a (Left) cell-type-specific regulatory networks were built by linking genes to NeuN+ and NeuN− specific H3K27ac peaks via Hi–C interactions in NeuN+ and NeuN− cells, respectively. (Right) The number of cell-type-specific peaks and their assigned genes in NeuN+ and NeuN− cells is described. A neuronal gene, HOMER1, is engaged with NeuN+ specific H3K27ac peaks via loops in NeuN+ cells (b), while an oligodendrocyte gene, SOX10, is engaged with NeuN− specific H3K27ac peaks via loops in NeuN− cells (c). The regions that interact with the gene promoter (gray) are highlighted in green (NeuN+) and purple (NeuN−), respectively. Boxplots in the right show expression levels of HOMER1 (FDR = 7.26e−32) and SOX10 (FDR = 1.88e−49) in NeuN+ (n = 4) and NeuN− (n = 4) cells, respectively. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. *FDR < 0.05 calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. d Genes assigned to NeuN+ specific peaks are enriched for synaptic co-expression modules, while genes assigned to NeuN− specific peaks are enriched for co-expression modules involved in transcriptional regulation and immune response during neurodevelopment. Significant enrichment (Sig.), FDR < 0.05. Fisher’s exact test was used for statistics analysis. OR, odds ratio. e Genes assigned to NeuN+ specific peaks are more highly enriched for synaptic functions such as exocytosis, intracellular signal transduction, protein cluster and structural plasticity than genes assigned to NeuN− specific peaks. Sig., FDR < 0.05. Fisher’s exact test was used for statistics analysis. f Genes assigned to NeuN+ specific peaks are highly expressed in neurons, while genes assigned to NeuN− specific peaks are highly expressed in oligodendrocytes and astrocytes. Astro astrocytes, Micro microglia, Endo Endothelial, Oligo oligodendrocytes.CC BY · PubMed Central
Fig. 3 Cell-type-specific nature of epigenetic dysregulation in AD. a We built AD-associated gene regulatory networks by linking genes to hypoacetylated (hypo) and hyperacetylated (hyper) peaks in AD via Hi–C interactions in NeuN+ and NeuN− cells. b AD-associated hyperacetylated peaks were largely active in NeuN− cells, while AD-associated hypoacetylated peaks are largely active in NeuN+ cells in neurotypical controls. c The number of genes mapped to AD-associated hyperacetylated (top) and hypoacetylated (bottom) peaks via Hi–C interactions in NeuN+ and NeuN− cells. CACNG3 is linked to an AD-associated hypoacetylated peak (marked in yellow) in NeuN+ cells (d), while EHD1 is linked to an AD-associated hyperacetylated peak (marked in yellow) in NeuN− cells (e). CACNG3/EHD1 promoter is highlighted in gray and its interacting regions are highlighted in green and purple for NeuN+ and NeuN− cells, respectively. Boxplots in the right show expression levels of CACNG3/EHD1 in NeuN+ (n = 4) and NeuN− (n = 4) cells. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. *FDR < 0.05 (FDR = 2.81e−9 for CACNG3 and FDR = 0.004 for EHD1), calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. f NeuN+ hypoacetylated genes are highly expressed in neurons, while NeuN− hyperacetylated genes are highly expressed in glia. g NeuN− hyperacetylated genes are enriched in astrocyte-specific co-expression modules (T-M14 and T-M8) that are upregulated in AD. NeuN+ hypoacetylated genes are enriched in neuronal co-expression modules (T-M1, T-M16) that are downregulated in AD. Fisher’s exact test was used for statistical analysis. The red line denotes FDR = 0.01. Astro astrocytes, Micro microglia, Endo endothelial, Oligo oligodendrocytes.
Fig. 3 Cell-type-specific nature of epigenetic dysregulation in AD. a We built AD-associated gene regulatory networks by linking genes to hypoacetylated (hypo) and hyperacetylated (hyper) peaks in AD via Hi–C interactions in NeuN+ and NeuN− cells. b AD-associated hyperacetylated peaks were largely active in NeuN− cells, while AD-associated hypoacetylated peaks are largely active in NeuN+ cells in neurotypical controls. c The number of genes mapped to AD-associated hyperacetylated (top) and hypoacetylated (bottom) peaks via Hi–C interactions in NeuN+ and NeuN− cells. CACNG3 is linked to an AD-associated hypoacetylated peak (marked in yellow) in NeuN+ cells (d), while EHD1 is linked to an AD-associated hyperacetylated peak (marked in yellow) in NeuN− cells (e). CACNG3/EHD1 promoter is highlighted in gray and its interacting regions are highlighted in green and purple for NeuN+ and NeuN− cells, respectively. Boxplots in the right show expression levels of CACNG3/EHD1 in NeuN+ (n = 4) and NeuN− (n = 4) cells. FPKM Fragments Per Kilobase of transcript per Million mapped reads. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. *FDR < 0.05 (FDR = 2.81e−9 for CACNG3 and FDR = 0.004 for EHD1), calculated by DESeq2 (two-sided Wald test). Source data are provided as a Source Data file. f NeuN+ hypoacetylated genes are highly expressed in neurons, while NeuN− hyperacetylated genes are highly expressed in glia. g NeuN− hyperacetylated genes are enriched in astrocyte-specific co-expression modules (T-M14 and T-M8) that are upregulated in AD. NeuN+ hypoacetylated genes are enriched in neuronal co-expression modules (T-M1, T-M16) that are downregulated in AD. Fisher’s exact test was used for statistical analysis. The red line denotes FDR = 0.01. Astro astrocytes, Micro microglia, Endo endothelial, Oligo oligodendrocytes.CC BY · PubMed Central
Fig. 4 Identification and characterization of putative target genes of AD genetic risk factors by incorporating NeuN− chromatin interaction. a The SNP-based heritability enrichment of AD GWAS in differential NeuN+ and NeuN− peaks suggests glial enrichment. Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. b BIN1 promoter physically interacts with an AD GWS locus in a NeuN− specific manner. The regions that interact with BIN1 promoter (marked in gray) are highlighted in purple. c GO analysis for GWAS-guided AD risk genes identified by NeuN− H-MAGMA. The red line denotes FDR = 0.05. d AD risk genes are highly expressed in postnatal brain samples compared with prenatal samples. Pre prenatal (n = 410); Post, postnatal (n = 453). p = 6.93e−193, calculated by Wilcoxon Rank Sum test. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. Source data are provided as a Source Data file. e AD risk genes are highly expressed in microglia. f AD risk genes are significantly enriched for genes differentially regulated in AD microglia. Fisher’s exact test was used for statistics analysis. The red line denotes FDR = 0.01. g AD risk genes are enriched in a microglial co-expression module that is upregulated in AD. Fisher’s exact test was used for statistical analysis. The red line denotes FDR = 0.01. Astro astrocytes, Micro microglia, Endo endothelial, Oligo oligodendrocytes.
Fig. 4 Identification and characterization of putative target genes of AD genetic risk factors by incorporating NeuN− chromatin interaction. a The SNP-based heritability enrichment of AD GWAS in differential NeuN+ and NeuN− peaks suggests glial enrichment. Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. b BIN1 promoter physically interacts with an AD GWS locus in a NeuN− specific manner. The regions that interact with BIN1 promoter (marked in gray) are highlighted in purple. c GO analysis for GWAS-guided AD risk genes identified by NeuN− H-MAGMA. The red line denotes FDR = 0.05. d AD risk genes are highly expressed in postnatal brain samples compared with prenatal samples. Pre prenatal (n = 410); Post, postnatal (n = 453). p = 6.93e−193, calculated by Wilcoxon Rank Sum test. Center, median; box = Q1–Q3; minima, Q1 − 1.5 × IQR; maxima, Q3 + 1.5 × IQR. Source data are provided as a Source Data file. e AD risk genes are highly expressed in microglia. f AD risk genes are significantly enriched for genes differentially regulated in AD microglia. Fisher’s exact test was used for statistics analysis. The red line denotes FDR = 0.01. g AD risk genes are enriched in a microglial co-expression module that is upregulated in AD. Fisher’s exact test was used for statistical analysis. The red line denotes FDR = 0.01. Astro astrocytes, Micro microglia, Endo endothelial, Oligo oligodendrocytes.CC BY · PubMed Central
Fig. 5 Comparison of SCZ and BD risk genes. a The SNP-based heritability enrichment of SCZ and BD GWAS in differential NeuN+ and NeuN− peaks demonstrates neuronal enrichment (SCZ NeuN+, FDR = 9.33e-24; SCZ NeuN−, FDR = 2.14e-16; BD NeuN+, FDR=3.94e-15; BD NeuN−, FDR=6.09e-07). Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. b The SNP-based heritability enrichment of SCZ and BD GWAS in differential Glu and GABA peaks suggests that both Glu and GABA neurons are associated with the psychiatric disorders (SCZ GABA, FDR = 4.54e−9; SCZ Glu, FDR = 3.49e−10; BD GABA, FDR=1.10e−03; BD Glu, FDR = 3.63e−06). Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. c Neuronal subtype expression profiles of SCZ and BD risk genes detected by GABA and Glu H-MAGMA.
Fig. 5 Comparison of SCZ and BD risk genes. a The SNP-based heritability enrichment of SCZ and BD GWAS in differential NeuN+ and NeuN− peaks demonstrates neuronal enrichment (SCZ NeuN+, FDR = 9.33e-24; SCZ NeuN−, FDR = 2.14e-16; BD NeuN+, FDR=3.94e-15; BD NeuN−, FDR=6.09e-07). Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. b The SNP-based heritability enrichment of SCZ and BD GWAS in differential Glu and GABA peaks suggests that both Glu and GABA neurons are associated with the psychiatric disorders (SCZ GABA, FDR = 4.54e−9; SCZ Glu, FDR = 3.49e−10; BD GABA, FDR=1.10e−03; BD Glu, FDR = 3.63e−06). Enrichment ± standard error is depicted. The red broken line indicates heritability enrichment = 1. Source data are provided as a Source Data file. c Neuronal subtype expression profiles of SCZ and BD risk genes detected by GABA and Glu H-MAGMA.CC BY · PubMed Central

Verification

Claims the paper makes

54 Own result · 45 Not checked yet · 22 Own interpretation · 3 Partly supported · 1 Supported · accuracy not yet measured

Verify all 69 references — free account

Show all 68 claims

  • It may state two separate claims as one. 51
  • Probably not a substantial claim. 7
57 of 57
  • The paper relies on this without support.

    Changes in gene expression or histone modifications may be a consequence or compensatory in disease and not necessarily causal.

    “Changes in gene expression or histone modifications may be a consequence or compensatory in disease and not necessarily causal.”
  • The paper relies on this without support.

    How AD-associated H3K27ac changes relate to cell-type-specific vulnerability is not known.

    “Studies in AD have revealed changes in gene regulation manifested by differences in H3K27ac in bulk tissue38, but how this relates to cell-type-specific vulnerability is not known.”
  • The paper relies on this without support.

    Hi-C data from neuronal subtypes is unavailable.

    “Since Hi–C data from neuronal subtypes is yet unavailable, we constructed Glu and GABA H-MAGMA maps by integrating Glu- and GABA-enhancers with NeuN+ Hi–C data, respectively (“Methods”)”
  • The paper relies on this without support.

    Genetic risk resides at the initiation of the causal chain in disease pathogenesis.

    “Since genetic risk resides at the initiation of the causal chain in disease pathogenesis, this differential enrichment of genetic risk factors and epigenetic regulation underscores different expression signatures in AD with distinct cellular specificities and likely causal relationships.”
4 of 4

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

3 of 3

Citations

Cites

Research works printed but not linked: 1

  1. Stark, R. et al. DiffBind: differential binding analysis of ChIP-Seq peak data. R. Package Version 100, 4–3 (2011).

Cited by

Sorted by citations

See all 103

Metadata

Authors
Hu, B., Won, H., Mah, W., Park, R. B., Kassim, B., Spiess, K., Kozlenkov, A., Crowley, C. A., Pochareddy, S., The PsychENCODE Consortium, Ashley-Koch, A. E., Crawford, G. E., Garrett, M. E., Song, L., Safi, A., Johnson, G. D., Wray, G. A., Reddy, T. E., Goes, F. S., Zandi, P., Bryois, J., Jaffe, A. E., Price, A. J., Ivanov, N. A., Collado-Torres, L., Hyde, T. M., Burke, E. E., Kleiman, J. E., Tao, R., Shin, J. H., Girdhar, K., Jiang, Y., Kundakovic, M., Brown, L., Wiseman, J. R., Zharovsky, E., Jacobov, R., Devillers, O., Flatow, E., Hoffman, G. E., Belmont, J., Del Valle, D., Francoeur, N., Hadjimichael, E., Pinto, D., van Bakel, H., Roussos, P., Fullard, J. F., Bendl, J., Hauberg, M. E., Charney, A. W., Haroutunian, V., Lipska, B. K., Lewis, D. A., Hahn, C.-G., Mangravite, L. M., Peters, M. A., Chae, Y., Peng, J., Niu, M., Wang, X., Webster, M. J., Beach, T. G., Chen, C., Jiang, Y., Dai, R., Wang, Y., Xia, Y., Shieh, A. W., Liu, C., Grennan, K. S., Vadukapuram, R., Giase, G., Fitzgerald, D., Cheng, L., Brown, M., Brown, M., Brunetti, T., Goodman, T., Alsayed, M., White, K. P., Ray, M., Polioudakis, D., Wamsley, B., Yin, J., De La Torre Ubieta, L., Gandal, M. J., Swarup, V., Sanders, S. J., State, M. W., Werling, D. M., An, J.-Y., Sheppard, B., Willsey, A. J., Kefi, A., Mattei, E., Purcaro, M., Weng, Z., Moore, J., Pratt, H., Huey, J., Borrman, T., Sullivan, P. F., Giusti-Rodriguez, P., Kim, Y., Szatkiewicz, J., Rhie, S. K., Armoskus, C., Camarena, A., Farnham, P. J., Spitsyna, V. N., Witt, H., Schreiner, S., Evgrafov, O. V., Knowles, J. A., Gerstein, M., Liu, S., Navarro, F. C. P., Warrell, J., Clarke, D., Emani, P. S., Gu, M., Shi, X., Xu, M., Yang, Y. T., Kitchen, R. R., Gürsoy, G., Zhang, J., Carlyle, B. C., Nairn, A. C., Li, M., Skarica, M., Li, Z., Sousa, A. M. M., Santpere, G., Choi, J., Zhu, Y., Gao, T., Miller, D. J., Cherskov, A., Yang, M., Amiri, A., Coppola, G., Mariani, J., Scuderi, S., Szekely, A., Vaccarino, F. M., Wu, F., Weissman, S., Wang, D., Roychowdhury, T., Abyzov, A., Li, Y., Dracheva, S., Sestan, N., Akbarian, S., & Geschwind, D. H.
Affiliations
  • University of North Carolina at Chapel Hill (US)
  • Icahn School of Medicine at Mount Sinai (US)
  • James J. Peters VA Medical Center (US)
  • Yale University (US)
  • Duke University (US)
  • Johns Hopkins University (US)
  • Karolinska Institutet (SE)
  • Lieber Institute for Brain Development (US)
  • Central South University (CN)
  • National Institutes of Health (US)
  • University of Pittsburgh (US)
  • University of Pennsylvania (US)
  • Sage Bionetworks (US)
  • St. Jude Children's Research Hospital (US)
  • Stanley Medical Research Institute (US)
  • Banner Sun Health Research Institute (US)
  • SUNY Upstate Medical University (US)
  • University of Chicago (US)
  • University of California, Los Angeles (US)
  • University of California, San Francisco (US)
  • University of Illinois Chicago (US)
  • University of Massachusetts Chan Medical School (US)
  • University of Southern California (US)
  • SUNY Downstate Health Sciences University (US)
  • University of Wisconsin–Madison (US)
  • Mayo Clinic (US)
Corresponding author
  • Won, Hyejung
  • Geschwind, Daniel H.
ORCID
Venue
Nature Communications
Published
2021-06-25
Volume
12
Issue
1
Pages
3968-3968
ISSN
2041-1723 (ISSN-L)
OpenAlex
W3173670503
Cited by
103
References
86
Citations by year
2026: 6 · 2025: 21 · 2024: 24 · 2023: 22 · 2022: 24 · 2021: 5
FWCI
5.96
Type
journal article
Language
en
Open access
Locations
  • Nature Communications (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • eScholarship@UMassChan (submittedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • eScholarship (California Digital Library) (submittedVersion)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
  • Repositori digital de la UPF (Universitat Pompeu Fabra) (submittedVersion, cc-by)
Funders
  • U.S. Department of Health & Human Services | NIH | National Institute of Mental Health (R00MH113823, DP2MH122403, P50MH106438, R01MH094714, U01MH103339, R01MH110927, R01MH100027)
  • U.S. Department of Health & Human Services | NIH | National Institute on Drug Abuse (R21DA051921)
  • U.S. Department of Health & Human Services | NIH | National Institute on Aging (R01AG066871)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Associations and Epidemiology
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
Keywords
3D genome organization, chromosome conformation, neurons, glia, oligodendrocytes, microglia, Alzheimer's disease, schizophrenia, bipolar disorder, enhancers, non-coding variants, parvalbumin interneurons, glutamatergic neurons, GABAergic neurons, epigenetic dysregulation
MeSH terms
Chromatin Immunoprecipitation Sequencing, Acetylation, Alzheimer Disease, Bipolar Disorder, Chromatin, Enhancer Elements, Genetic, Gene Expression Regulation, Histones, Humans, Lysine, Neuroglia, Neurons, Promoter Regions, Genetic, Schizophrenia, Epigenesis, Genetic, Genome-Wide Association Study, GABAergic Neurons

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-10, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-10 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review