CoChat
Sign Up
Cited by 492Single-cell and spatial transcriptomicsGene expression and cancer classificationBioinformatics and Genomic Networks

The DOI matches its Crossref record; 59 of the 59 references it cites were found in a registry or an index.

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

Open access
McKenzie, Andrew T., Wang, Minghui, Hauberg, Mads E., Fullard, John F., Kozlenkov, Alexey, Keenan, Alexandra, Hurd, Yasmin L., Dracheva, Stella, Casaccia, Patrizia, Roussos, Panos, Zhang, Bin·2018·Scientific Reports
DOI 10.1038/s41598-018-27293-5PMC PMC5995803
Open this paper in CoChat — free accountOpen PDF

Abstract

Elucidating brain cell type specific gene expression patterns is critical towards a better understanding of how cell-cell communications may influence brain functions and dysfunctions. We set out to compare and contrast five human and murine cell type-specific transcriptome-wide RNA expression data sets that were generated within the past several years. We defined three measures of brain cell type-relative expression including specificity, enrichment, and absolute expression and identified corresponding consensus brain cell “signatures,” which were well conserved across data sets. We validated that the relative expression of top cell type markers are associated with proxies for cell type proportions in bulk RNA expression data from postmortem human brain samples. We further validated novel marker genes using an orthogonal ATAC-seq dataset. We performed multiscale coexpression network analysis of the single cell data sets and identified robust cell-specific gene modules. To facilitate the use of the cell type-specific genes for cell type proportion estimation and deconvolution from bulk brain gene expression data, we developed an R package, BRETIGEA. In summary, we identified a set of novel brain cell consensus signatures and robust networks from the integration of multiple datasets and therefore transcend limitations related to technical issues characteristic of each individual study.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41598-018-27293-5

Cites

59 of 59 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Explanation of the three cell type associated measures. Diagrams showing the three cell type-associated measures used in this study. (a) Cell type absolute expression simply measures the relative expression of each gene in each cell type, irrespective of the expression of that gene in other cell types. (b) Cell type enrichment measures the expression of each gene relative to the expression of that gene in all other cell types. With this measure, a gene could have relatively high expression in two cell types, and be relatively enriched in each of them compared to all other cell types. (c) Cell type specificity measures the expression of each gene relative to the highest expression of that gene in all other cell types. This measure requires that the expression of the gene is only high in one cell type; therefore, we call it “specific.”
Figure 1 Explanation of the three cell type associated measures. Diagrams showing the three cell type-associated measures used in this study. (a) Cell type absolute expression simply measures the relative expression of each gene in each cell type, irrespective of the expression of that gene in other cell types. (b) Cell type enrichment measures the expression of each gene relative to the expression of that gene in all other cell types. With this measure, a gene could have relatively high expression in two cell types, and be relatively enriched in each of them compared to all other cell types. (c) Cell type specificity measures the expression of each gene relative to the highest expression of that gene in all other cell types. This measure requires that the expression of the gene is only high in one cell type; therefore, we call it “specific.”CC BY · PubMed Central
Figure 2 Differential expression volcano plots across cell types and data sets for cell type enrichment. (a) Summary of the five data sets used in this study indicating the region in which the cell types were isolated from in both the human (left) and mouse (right) brains. (b–f): Volcano plots show the cell type enrichment differential expression calculation for each cell type in each data set (b = Darmanis et al. c = Zhang et al., d = Zhang et al., e = Zeisel et al. f = Tasic et al.). In each small multiple, the x-axis refers to the log2 fold-change of the gene in the cell type of interest compared to the reference cell set, while the y-axis shows the −log10 Benjamini-Hochberg adjusted p-value of that comparison, calculated using edgeR. Genes with adjusted p-values < 0.05 and log fold-change > 2 are colored orange. OPC = Oligodendrocyte Precursor Cell.
Figure 2 Differential expression volcano plots across cell types and data sets for cell type enrichment. (a) Summary of the five data sets used in this study indicating the region in which the cell types were isolated from in both the human (left) and mouse (right) brains. (b–f): Volcano plots show the cell type enrichment differential expression calculation for each cell type in each data set (b = Darmanis et al. c = Zhang et al., d = Zhang et al., e = Zeisel et al. f = Tasic et al.). In each small multiple, the x-axis refers to the log2 fold-change of the gene in the cell type of interest compared to the reference cell set, while the y-axis shows the −log10 Benjamini-Hochberg adjusted p-value of that comparison, calculated using edgeR. Genes with adjusted p-values < 0.05 and log fold-change > 2 are colored orange. OPC = Oligodendrocyte Precursor Cell.CC BY · PubMed Central
Figure 3 Pairwise data set comparison of cell type enrichment rankings for each cell type. Plots of the fold-enrichment of the intersections of genes ranked in the top n genes (where n = 10, 20, 50, 100, 200, 500, 1000) between pairs of data sets for the cell type enrichment measure. The data sets were merged to only include gene symbols common to both prior to calculating the fold enrichment score. A fold enrichment of 0 indicates that no genes were found in the intersection of those two sets of top genes.
Figure 3 Pairwise data set comparison of cell type enrichment rankings for each cell type. Plots of the fold-enrichment of the intersections of genes ranked in the top n genes (where n = 10, 20, 50, 100, 200, 500, 1000) between pairs of data sets for the cell type enrichment measure. The data sets were merged to only include gene symbols common to both prior to calculating the fold enrichment score. A fold enrichment of 0 indicates that no genes were found in the intersection of those two sets of top genes.CC BY · PubMed Central
Figure 4 Correlation of genewise cell type enrichment measures with PubMed text mining results within cell types. For each of the six cell types, i.e. the astrocyte (a), endothelial cell (b), microglia (c), neuron (d), oligodendrocyte (e), and oligodendrocyte precursor cell (f), the top 100 gene symbols most enriched in that cell type are plotted against the number of PubMed abstracts that mention both that gene symbol as well as the corresponding cell type. The Spearman correlation between these measures was calculated. Several gene symbols were chosen for highlighting, including gene symbols that have not been mentioned in a PubMed abstract with that cell type to date (labeled red). Note that for oligodendrocyte precursor cells (OPCs), the cell type name used in the PubMed search was “oligodendrocyte precursor.”
Figure 4 Correlation of genewise cell type enrichment measures with PubMed text mining results within cell types. For each of the six cell types, i.e. the astrocyte (a), endothelial cell (b), microglia (c), neuron (d), oligodendrocyte (e), and oligodendrocyte precursor cell (f), the top 100 gene symbols most enriched in that cell type are plotted against the number of PubMed abstracts that mention both that gene symbol as well as the corresponding cell type. The Spearman correlation between these measures was calculated. Several gene symbols were chosen for highlighting, including gene symbols that have not been mentioned in a PubMed abstract with that cell type to date (labeled red). Note that for oligodendrocyte precursor cells (OPCs), the cell type name used in the PubMed search was “oligodendrocyte precursor.”CC BY · PubMed Central
Figure 5 Intersections among the top genes for three cell type associated measures consensus rankings across all cell types. The top 100 genes ranked across both mouse and human data sets for each of the cell type measures are intersected using approximately proportional Venn diagrams for each of the astrocyte (a), endothelial cell (b), microglia (c), neuron (d), oligodendrocyte (e), and oligodendrocyte precursor cell (f) signatures. The Venn diagrams were generated using the R package Vennerable (version 3.0). The 5 genes with the top expression values in each of the cell types that intersect in all three of the top 100 gene sets are listed, with the exception of OPCs, for which all 7 of the genes with intersections between the three measures are listed.
Figure 5 Intersections among the top genes for three cell type associated measures consensus rankings across all cell types. The top 100 genes ranked across both mouse and human data sets for each of the cell type measures are intersected using approximately proportional Venn diagrams for each of the astrocyte (a), endothelial cell (b), microglia (c), neuron (d), oligodendrocyte (e), and oligodendrocyte precursor cell (f) signatures. The Venn diagrams were generated using the R package Vennerable (version 3.0). The 5 genes with the top expression values in each of the cell types that intersect in all three of the top 100 gene sets are listed, with the exception of OPCs, for which all 7 of the genes with intersections between the three measures are listed.CC BY · PubMed Central
Figure 6 Cell type relative proportion estimates from bulk RNA expression for astrocytes and microglia are associated with IHC quantifications from the Allen Brain Atlas brain bank. (a-b) Scatter plots of the relative cell type proportion estimates, generated using a singular value decomposition (SVD) of the top 50 marker genes for astrocytes (a) or microglia (b) (x-axis) with the immunohistochemistry (IHC) protein level quantifications for the astrocyte marker GFAP (a) or the microglia IBA1 (b) (y-axis) in the same donors across the four brain regions in the Allen Brain Atlas Aging, Dementia, and TBI data set. The black line is a result of a linear model fit to the data, while the grey lines represent 95% confidence intervals. (c,d) Rho values resulting from rank correlations of the RNA expression of the top 100 marker genes individually (black dots) as well as the SVD-based estimate using the cumulative astrocyte (c) or microglia (d) marker genes up to that marker gene, inclusive (red lines), with the IHC quantifications for GFAP (c) or IBA1 (d). FWM, frontal white matter; HIP, hippocampus; PCx, parietal cortex; TCx, temporal cortex.
Figure 6 Cell type relative proportion estimates from bulk RNA expression for astrocytes and microglia are associated with IHC quantifications from the Allen Brain Atlas brain bank. (a-b) Scatter plots of the relative cell type proportion estimates, generated using a singular value decomposition (SVD) of the top 50 marker genes for astrocytes (a) or microglia (b) (x-axis) with the immunohistochemistry (IHC) protein level quantifications for the astrocyte marker GFAP (a) or the microglia IBA1 (b) (y-axis) in the same donors across the four brain regions in the Allen Brain Atlas Aging, Dementia, and TBI data set. The black line is a result of a linear model fit to the data, while the grey lines represent 95% confidence intervals. (c,d) Rho values resulting from rank correlations of the RNA expression of the top 100 marker genes individually (black dots) as well as the SVD-based estimate using the cumulative astrocyte (c) or microglia (d) marker genes up to that marker gene, inclusive (red lines), with the IHC quantifications for GFAP (c) or IBA1 (d). FWM, frontal white matter; HIP, hippocampus; PCx, parietal cortex; TCx, temporal cortex.CC BY · PubMed Central
Figure 7 Relative chromatin accessibility in (a) novel and (b) known brain cell type marker genes by ATAC-seq. Relative chromatin accessibility, as measured by Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq), at the promoters of novel candidate cell marker genes in four different ATAC-seq cell populations from the adult human brain. “Marker” genes were defined as top cell type-specific genes for a particular cell type. In each panel, the boxplots show the distribution of the fraction of ATAC-seq reads regarding novel markers of a particular cell type. Astro, astrocytes.
Figure 7 Relative chromatin accessibility in (a) novel and (b) known brain cell type marker genes by ATAC-seq. Relative chromatin accessibility, as measured by Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq), at the promoters of novel candidate cell marker genes in four different ATAC-seq cell populations from the adult human brain. “Marker” genes were defined as top cell type-specific genes for a particular cell type. In each panel, the boxplots show the distribution of the fraction of ATAC-seq reads regarding novel markers of a particular cell type. Astro, astrocytes.CC BY · PubMed Central
Figure 8 Multiscale networks identified within neurons contain neuron-specific modules. (a) Multiscale network modules identified using MEGENA in the Darmanis et al. human neuron-annotated RNA-seq samples. Dots represent genes while lines represent network connections, and dots are colored by their presence in one of several modules. (b) Zoomed-in gene-level network for the genes in Module #57, which is significantly enriched in the GO term “neurotrophin TRK receptor signaling pathway”.
Figure 8 Multiscale networks identified within neurons contain neuron-specific modules. (a) Multiscale network modules identified using MEGENA in the Darmanis et al. human neuron-annotated RNA-seq samples. Dots represent genes while lines represent network connections, and dots are colored by their presence in one of several modules. (b) Zoomed-in gene-level network for the genes in Module #57, which is significantly enriched in the GO term “neurotrophin TRK receptor signaling pathway”.CC BY · PubMed Central
Figure 9 More significant module intersections are found within than between cell types across human to mouse data set comparisons. The proportion (between 0 and 0.5) of modules with significant intersections between MEGENA multiscale models is plotted for the human Darmanis et al. to mouse Tasic et al. (a) and human Darmanis et al. to mouse Zeisel et al. (b) comparisons. A module-module intersection between two cell type-associated multiscale networks was counted as a significant intersection (or overlap) if its Benjamini-Hochberg-adjusted Fisher’s Exact Test p-value was < 0.05. To generate proportions for the number of overlapping modules relative to each mouse cell type module set, the overlap matrix was normalized by column.
Figure 9 More significant module intersections are found within than between cell types across human to mouse data set comparisons. The proportion (between 0 and 0.5) of modules with significant intersections between MEGENA multiscale models is plotted for the human Darmanis et al. to mouse Tasic et al. (a) and human Darmanis et al. to mouse Zeisel et al. (b) comparisons. A module-module intersection between two cell type-associated multiscale networks was counted as a significant intersection (or overlap) if its Benjamini-Hochberg-adjusted Fisher’s Exact Test p-value was < 0.05. To generate proportions for the number of overlapping modules relative to each mouse cell type module set, the overlap matrix was normalized by column.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

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.

The self-check needs the paper’s full text.

Citations

Cites

Show all 59 references

Cited by

Sorted by citations

Related

See all

Metadata

Authors
McKenzie, A. T., Wang, M., Hauberg, M. E., Fullard, J. F., Kozlenkov, A., Keenan, A., Hurd, Y. L., Dracheva, S., Casaccia, P., Roussos, P., & Zhang, B.
Affiliations
  • Icahn School of Medicine at Mount Sinai (US)
  • Aarhus University (DK)
  • James J. Peters VA Medical Center (US)
  • City University of New York (US)
Corresponding author
Zhang, Bin
ORCID
Venue
Scientific Reports
Publisher
Springer Science and Business Media LLC
Published
2018-06-11
Volume
8
Issue
1
Pages
8868-8868
ISSN
2045-2322 (ISSN-L)
DOI
10.1038/s41598-018-27293-5
PMC
PMC5995803
OpenAlex
W2806989085
Cited by
481
References
59
Citations by year
2026: 41 · 2025: 59 · 2024: 84 · 2023: 70 · 2022: 60 · 2021: 78 · 2020: 56 · 2019: 28 · 2018: 5
FWCI
15.03
Type
journal article
Language
en
Update notices
correction (2021-09-24)
Open access
Licence
Locations
  • Scientific Reports (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Scopus (Elsevier) (submittedVersion)
Funders
  • U.S. Department of Veterans Affairs (I01 BX002876, I01BX002395, I01 BX003625)
  • Lundbeckfonden (R155-2014-1724, R248‐2017‐2003)
  • National Institutes of Health
  • National Institute on Aging (R01AG057907, F30AG052261, R01AG046170, U01AG046170, RF1AG054014, RF1AG057440, U01AG052411, R01AG050986)
  • National Institute of Mental Health (R01MH109677)
  • National Institute on Drug Abuse (P01DA008227)
  • National Institute of Neurological Disorders and Stroke (R37NS042925, R01NS052738)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
Keywords
cell type markers, co-expression network, single-cell transcriptomics, cell type deconvolution, cell type signatures, intercellular communication, postmortem brain
MeSH terms
Animals, Female, Genetic Markers, Humans, Male, Mice, Inbred C57BL, Neuroglia, Neurons, Temporal Lobe, Gene Expression Profiling, Databases, Nucleic Acid, Endothelial Cells, Mice, Gene Regulatory Networks, Single-Cell Analysis, Transcriptome, Datasets as Topic

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.

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