Cited by 477Single-cell and spatial transcriptomicsNeuroinflammation and Neurodegeneration MechanismsMolecular Biology Techniques and Applications

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

Dissection of artifactual and confounding glial signatures by single-cell sequencing of mouse and human brain

Open access

Published version of a preprint

Abstract

A key aspect of nearly all single-cell sequencing experiments is dissociation of intact tissues into single-cell suspensions. While many protocols have been optimized for optimal cell yield, they have often overlooked the effects that dissociation can have on ex vivo gene expression. Here, we demonstrate that use of enzymatic dissociation on brain tissue induces an aberrant ex vivo gene expression signature, most prominently in microglia, which is prevalent in published literature and can substantially confound downstream analyses. To address this issue, we present a rigorously validated protocol that preserves both in vivo transcriptional profiles and cell-type diversity and yield across tissue types and species. We also identify a similar signature in postmortem human brain single-nucleus RNA-sequencing datasets, and show that this signature is induced in freshly isolated human tissue by exposure to elevated temperatures ex vivo. Together, our results provide a methodological solution for preventing artifactual gene expression changes during fresh tissue digestion and a reference for future deeper analysis of the potential confounding states present in postmortem human samples.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

Enzymatic dissociation of mouse brain induces an aberrant ex vivo transcriptional signature, most prominently in microglia, that is widespread in published single-cell datasets and can confound downstream analyses. The authors provide an inhibitor-based protocol that prevents this signature and show a similar signature in human postmortem microglia/astrocytes that can be induced by delayed freezing.

Key findings

  • In sorted mouse microglia, enzymatic digestion without inhibitors produced an ex vivo “activated” microglia (exAM) cluster enriched for immediate early, stress and immune genes; ENZ-NONE samples had a significant decrease in homeostatic cluster proportion and increase in exAM proportion (*FDR ≤ 0.004; ENZ-NONE versus all other groups; 19,563 cells from n = 12 mice, n = 3 per condition).
  • Adding transcriptional/translational inhibitors eliminated the artifact: DNC-NONE versus DNC-INHIB showed no significant DEGs, and in all-CNS scRNA-seq (10,166 cells, n = 2 per group) only 18 genes were significantly upregulated and none downregulated without inhibitors.
  • Reanalysis of published datasets found the exAM signature across technologies, enzymes and brain regions; overlap was 19 of 25 genes with Van Hove et al., 16 genes with mouse muscle and 11 genes with kidney dissociation signatures.
  • Human postmortem snRNA-seq (47,505 nuclei, three donors) showed enrichment in microglia and some astrocytes; LIGER factors overlapped 79% (30 of 38) for microglia and 44% (12 of 27) for astrocytes, and a 6-h freezing delay induced the microglial factor (P < 3.2 × 10−215; 31 of 38 genes) but weaker astrocyte changes (P < 0.05; 4 of 27 genes).

Methods

  • Mouse scRNA-seq compared cold mechanical dissociation with enzymatic digestion ± an inhibitor cocktail, using FACS-sorted CD45+/CD11b+ myeloid cells and unsorted all-CNS preparations; analyses used Seurat, DESeq2 on pseudobulk “Metacells,” MAST and speckle.
  • Human snRNA-seq used postmortem brain and acutely resected peri-tumor tissue, with nuclei isolation/sorting and 10X Genomics 3′ v.3; analyses used LIGER/iNMF, gene module scoring and comparison of snap-frozen versus 6-h delayed-freezing samples.
  • The study also reanalyzed published mouse and human datasets and performed RNAscope smFISH for Ccl4/Fcrls validation.

Limitations

  • Stated: the human postmortem signature was highly variable across donors and did not correlate with PMI or most available metadata; an age correlation (Spearman r = 0.3208, *P = 0.0001, n = 48) was not robust in donors aged 75+ (r = 0.1879, P = 0.2653, n = 37), and more middle-aged samples are needed.
  • Stated: current human snRNA-seq datasets are limited by microglia proportion and lower sensitivity of earlier chemistries/nuclei isolation, and the authors caution against overinterpreting biological significance of similar signatures.
  • Evident from the text: mouse scRNA-seq used male mice at P89–90 with small group sizes (n = 2–3 per group), human PBMC mock digestion used one healthy 35-year-old male, and acute human tissue was peri-tumor rather than normal brain; the authors also note that one proliferative-cluster abundance result should be interpreted cautiously because it represented only 0.1% of cells (27 of 19,563).

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

Integrity signals

Registry

Crossref · 10.1038/s41593-022-01022-8

Cites

115 of 115 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Analysis of sorted microglia confirms profound effect of enzymatic digestion on microglial gene expression via scRNA-seq. a, Experimental design schematic for sorted mouse myeloid cells scRNA-seq experiment (Methods and Supplementary Fig. 1a–c). b, t-SNE plot for the 19,563 cells from n = 12 mice (n = 3 per group) colored and annotated by cluster (Supplementary Fig. 2). c, t-SNE plot split by experimental subgroup highlights enrichment of exAM cluster (circled) in ENZ-NONE group. d, Gene expression of several exAM cluster markers across each cluster. e, Heatmap of the mean percentage of cells in each cluster across conditions (*FDR < 0.005 for ENZ-NONE versus all other groups; #FDR < 0.05 DNC-NONE versus ENZ-INHIB; Benjamin and Hochberg correction for multiple comparisons; (Supplementary Table 3); Methods: differential abundance testing). f,g, Visualization of gene module scoring results plotted on t-SNE coordinates. f, Microglial identity score (Methods and Supplementary Table 4). g, Activation score based on consensus DEGs from ‘Metacell’ pseudobulk analysis (Methods and Supplementary Tables 4 and 5). h, Plot of microglial identity score versus activation score colored by cluster annotation from panel b. i, smFISH using RNAscope for microglial marker Fcrls (green), cluster marker Ccl4 and counterstained with DAPI; representative images from n = 2 independent experiments. Scale bar, 50 μm. Chem, Chemokine; hom, Homeostatic; mac, macrophage; mono, monocyte; prolif, proliferative.
Fig. 1 Analysis of sorted microglia confirms profound effect of enzymatic digestion on microglial gene expression via scRNA-seq. a, Experimental design schematic for sorted mouse myeloid cells scRNA-seq experiment (Methods and Supplementary Fig. 1a–c). b, t-SNE plot for the 19,563 cells from n = 12 mice (n = 3 per group) colored and annotated by cluster (Supplementary Fig. 2). c, t-SNE plot split by experimental subgroup highlights enrichment of exAM cluster (circled) in ENZ-NONE group. d, Gene expression of several exAM cluster markers across each cluster. e, Heatmap of the mean percentage of cells in each cluster across conditions (*FDR < 0.005 for ENZ-NONE versus all other groups; #FDR < 0.05 DNC-NONE versus ENZ-INHIB; Benjamin and Hochberg correction for multiple comparisons; (Supplementary Table 3); Methods: differential abundance testing). f,g, Visualization of gene module scoring results plotted on t-SNE coordinates. f, Microglial identity score (Methods and Supplementary Table 4). g, Activation score based on consensus DEGs from ‘Metacell’ pseudobulk analysis (Methods and Supplementary Tables 4 and 5). h, Plot of microglial identity score versus activation score colored by cluster annotation from panel b. i, smFISH using RNAscope for microglial marker Fcrls (green), cluster marker Ccl4 and counterstained with DAPI; representative images from n = 2 independent experiments. Scale bar, 50 μm. Chem, Chemokine; hom, Homeostatic; mac, macrophage; mono, monocyte; prolif, proliferative.CC BY · PubMed Central
Fig. 2 Enzymatic dissociation induces cell-type-specific artifactual gene expression in mice. a, Experimental design schematic for scRNA-seq of all CNS cell types (Methods). b, t-SNE plot of 10,166 cells from 4 mice (n = 2 per group), annotated and colored by cell type (Supplementary Fig. 7). c, Expression of microglial-specific genes Tmem119 and P2ry12 clearly defines two microglial clusters. d, t-SNE plot split by experimental group (±inhibitor cocktail) highlights lack of exAM cluster in samples digested with inhibitors present (bottom; circled). e, Heatmap displaying mean percentage of cells in each cluster across conditions (*FDR < 0.0005; Benjamini and Hochberg correction for multiple comparisons; Supplementary Table 18; Methods: differential abundance testing). f,g, Gene module scoring results plotted on t-SNE coordinates. f, Microglial identity score (Methods and Supplementary Table 4). g, Activation score based on DEGs from ‘Metacell’ pseudobulk analysis (Methods and Supplementary Table 4). h, Scatterplot of gene module scores from f and g colored by cluster from panel b. i, Enrichment of genes from activation score in exAM microglia cluster. CP, choroid plexus; NPC, neural progenitor cells; NSC, neural stem cells; OEC, olfactory ensheathing cells; OPC, oligodendrocyte progenitor cells; RBC, red blood cells.
Fig. 2 Enzymatic dissociation induces cell-type-specific artifactual gene expression in mice. a, Experimental design schematic for scRNA-seq of all CNS cell types (Methods). b, t-SNE plot of 10,166 cells from 4 mice (n = 2 per group), annotated and colored by cell type (Supplementary Fig. 7). c, Expression of microglial-specific genes Tmem119 and P2ry12 clearly defines two microglial clusters. d, t-SNE plot split by experimental group (±inhibitor cocktail) highlights lack of exAM cluster in samples digested with inhibitors present (bottom; circled). e, Heatmap displaying mean percentage of cells in each cluster across conditions (*FDR < 0.0005; Benjamini and Hochberg correction for multiple comparisons; Supplementary Table 18; Methods: differential abundance testing). f,g, Gene module scoring results plotted on t-SNE coordinates. f, Microglial identity score (Methods and Supplementary Table 4). g, Activation score based on DEGs from ‘Metacell’ pseudobulk analysis (Methods and Supplementary Table 4). h, Scatterplot of gene module scores from f and g colored by cluster from panel b. i, Enrichment of genes from activation score in exAM microglia cluster. CP, choroid plexus; NPC, neural progenitor cells; NSC, neural stem cells; OEC, olfactory ensheathing cells; OPC, oligodendrocyte progenitor cells; RBC, red blood cells.CC BY · PubMed Central
Fig. 3 snRNA-seq of human postmortem tissue identifies enrichment of mouse dissociation gene signatures in human microglia and astrocytes. a, Experimental design schematic for snRNA-seq of all cell types from frozen postmortem brain tissue. b, UMAP plot of 47,505 nuclei analyzed via snRNA-seq from three postmortem subjects following LIGER analysis, colored by major cell type. c, Expression of canonical marker genes delineates major cell types. d, Visualization of the gene module scoring of the mouse DEG signature on human postmortem snRNA-seq dataset. e, Gene expression of exAM signature gene FOS across clusters. f–k, Plots of mouse activation score versus sample PMI (y axis) for each of the major CNS cell classes present in the dataset: f, microglia, g, astrocytes, h, oligodendrocytes, i, oligodendrocyte precursor cells, j, excitatory neurons, k, inhibitory neurons. f-k, Percentages denote number of nuclei above enrichment threshold denoted with gray dotted line; individual points are colored by donor.
Fig. 3 snRNA-seq of human postmortem tissue identifies enrichment of mouse dissociation gene signatures in human microglia and astrocytes. a, Experimental design schematic for snRNA-seq of all cell types from frozen postmortem brain tissue. b, UMAP plot of 47,505 nuclei analyzed via snRNA-seq from three postmortem subjects following LIGER analysis, colored by major cell type. c, Expression of canonical marker genes delineates major cell types. d, Visualization of the gene module scoring of the mouse DEG signature on human postmortem snRNA-seq dataset. e, Gene expression of exAM signature gene FOS across clusters. f–k, Plots of mouse activation score versus sample PMI (y axis) for each of the major CNS cell classes present in the dataset: f, microglia, g, astrocytes, h, oligodendrocytes, i, oligodendrocyte precursor cells, j, excitatory neurons, k, inhibitory neurons. f-k, Percentages denote number of nuclei above enrichment threshold denoted with gray dotted line; individual points are colored by donor.CC BY · PubMed Central
Fig. 4 LIGER analysis independently identifies similar gene expression signatures in postmortem data that are enriched in microglia following altered sample processing. a, UMAP plot visualizing the enrichment of shared LIGER factor for 12,790 microglial nuclei from 48 samples across all postmortem datasets. b, UMAP plot visualizing the enrichment of shared LIGER factor for 23,998 astrocyte nuclei from 49 samples across all postmortem datasets. For both a and b, inset displays plot of normalized cell-specific factor loading scores across all genes in the dataset (dashed line indicates threshold cutoff for top genes for downstream analysis; Supplementary Table 22). Top loading genes, in order, are shown to the right of the inset plot. c, Experimental design schematic for experiment to analyze the effects of altered sample processing on gene expression. d–f, Visualization of gene module scoring results for score based on postmortem microglia factor, from a, in both snap-frozen (d) and 6-h delayed freezing (e) microglia nuclei on UMAP coordinates or via violin plot split by experimental group (f). g, Gene expression of top 12 loading genes in microglial factor from a split by experimental group (Supplementary Table 23 for DEG results). h–j, Visualization of gene module scoring results for score based on postmortem astrocyte factor, from b, in both snap-frozen (h) and 6-h delayed freezing (i) astrocyte nuclei on UMAP coordinates or via violin plot split by experimental group (j). k, Gene expression of top 12 loading genes in astrocyte factor from b split by experimental group (Supplementary Table 24 for DEG results). l,m, Spearman correlation of the percentage of microglia above score threshold for microglia factor score in each postmortem sample versus l, PMI and m, age of donor; graph annotations list Spearman r values and significance. n, Plot of percentage of microglia above score threshold for microglia factor score in each postmortem sample split by diagnosis; n = 48 independent samples from 5 studies/independent experiments (45 samples are from 4 previously published studies); n = 32 control/n = 16 Alzheimer’s disease. Data are presented as mean values ± s.e.m. AD, Alzheimer’s disease; NS, not significant; RT, room temperature.
Fig. 4 LIGER analysis independently identifies similar gene expression signatures in postmortem data that are enriched in microglia following altered sample processing. a, UMAP plot visualizing the enrichment of shared LIGER factor for 12,790 microglial nuclei from 48 samples across all postmortem datasets. b, UMAP plot visualizing the enrichment of shared LIGER factor for 23,998 astrocyte nuclei from 49 samples across all postmortem datasets. For both a and b, inset displays plot of normalized cell-specific factor loading scores across all genes in the dataset (dashed line indicates threshold cutoff for top genes for downstream analysis; Supplementary Table 22). Top loading genes, in order, are shown to the right of the inset plot. c, Experimental design schematic for experiment to analyze the effects of altered sample processing on gene expression. d–f, Visualization of gene module scoring results for score based on postmortem microglia factor, from a, in both snap-frozen (d) and 6-h delayed freezing (e) microglia nuclei on UMAP coordinates or via violin plot split by experimental group (f). g, Gene expression of top 12 loading genes in microglial factor from a split by experimental group (Supplementary Table 23 for DEG results). h–j, Visualization of gene module scoring results for score based on postmortem astrocyte factor, from b, in both snap-frozen (h) and 6-h delayed freezing (i) astrocyte nuclei on UMAP coordinates or via violin plot split by experimental group (j). k, Gene expression of top 12 loading genes in astrocyte factor from b split by experimental group (Supplementary Table 24 for DEG results). l,m, Spearman correlation of the percentage of microglia above score threshold for microglia factor score in each postmortem sample versus l, PMI and m, age of donor; graph annotations list Spearman r values and significance. n, Plot of percentage of microglia above score threshold for microglia factor score in each postmortem sample split by diagnosis; n = 48 independent samples from 5 studies/independent experiments (45 samples are from 4 previously published studies); n = 32 control/n = 16 Alzheimer’s disease. Data are presented as mean values ± s.e.m. AD, Alzheimer’s disease; NS, not significant; RT, room temperature.CC BY · PubMed Central

Verification

Claims the paper makes

54 Own result · 18 Own method · 15 Own interpretation · 13 Partly supported · 6 Supported · 5 No full text of the sources · 4 Supported so far · 4 No basis given · 1 Couldn’t confirm · 1 Not checked yet · accuracy not yet measured

Show all 72 claims

  • Probably not a substantial claim. 27
  • It may state two separate claims as one. 26
  • The quote only partly states this claim. 1
49 of 49
  • The paper relies on this without support.

    No prior study has examined the human brain for cell-type-specific differences in transcriptional states using snRNA-seq.

    “and to our knowledge none have done so in a cell-type-specific manner using snRNA-seq”
  • The paper relies on this without support.

    The majority of studies have utilized enzymatic dissociation to study all CNS populations at single-cell resolution.

    “Consequently, the majority of studies have utilized enzymatic dissociation to study all CNS populations at single-cell resolution.”
2 of 2

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?

Citations

Cites

Show all 108 references

Cited by

Sorted by citations

Items 101–200 of 477

Metadata

Authors
Marsh, S. E., Walker, A. J., Kamath, T., Dissing-Olesen, L., Hammond, T. R., de Soysa, T. Y., Young, A. M. H., Murphy, S., Abdulraouf, A., Nadaf, N., Dufort, C., Walker, A. C., Lucca, L. E., Kozareva, V., Vanderburg, C., Hong, S., Bulstrode, H., Hutchinson, P. J., Gaffney, D. J., Hafler, D. A., Franklin, R. J. M., Macosko, E. Z., & Stevens, B.
Affiliations
  • Broad Institute (US)
  • Wellcome/MRC Cambridge Stem Cell Institute (GB)
  • Boston Children's Hospital (US)
  • Yale University (US)
  • UK Dementia Research Institute (GB)
  • University of Cambridge (GB)
  • Wellcome Sanger Institute (GB)
Corresponding author
Stevens, Beth
Venue
Nature Neuroscience
Published
2022-03
Volume
25
Issue
3
Pages
306-316
ISSN
  • 1097-6256 (ISSN-L)
  • 1546-1726
OpenAlex
W4220661922
Cited by
431
References
108
Citations by year
2026: 102 · 2025: 71 · 2024: 103 · 2023: 83 · 2022: 61
FWCI
26.69
Type
journal article
Language
en
Open access
Locations
  • Nature Neuroscience (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (submittedVersion, cc-by-nd)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • UCL Discovery (University College London) (submittedVersion)
Funders
  • Dr. Miriam and Sheldon G. Adelson Medical Research Foundation (203151/Z/16/Z)
  • Cure Alzheimer's Fund
  • Rush University (P30AG10161, R01AG36042, R01 AG15819, R01AG30146)
  • Wellcome Trust (203151/Z/16/Z, 203151/Z/16, 203151)
  • Cancer Research UK (24423)
  • National Institute for Health and Care Research (R01AG15819, CL-2016-14-005, P30AG10161)
  • Royal College of Surgeons of England
  • Multiple Sclerosis Society
  • Academy of Medical Sciences (SGL017\1009)
  • Lundbeckfonden (R182-2014-2844, R193-2015-4103)
  • H. Lundbeck A/S
  • University of California, Irvine
  • Medical Research Council (203151/Z/16/Z)
  • National Institute on Aging (R01AG017917, R01AG030146, R01AG036836, R01AG048015, U01AG046152, R01AG036042, R01AG015819, P30AG010161, U01AG032984, U01AG061356, RC2AG036547, U01AG046161, RF1AG057473)
  • National Institute of General Medical Sciences (T32GM008313)
  • National Institute of Neurological Disorders and Stroke (T32NS007473)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Neuroscience › Neurology › Neuroinflammation and Neurodegeneration Mechanisms
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Molecular Biology Techniques and Applications
Keywords
microglia, enzymatic tissue dissociation, single-nucleus RNA sequencing, human brain
MeSH terms
Brain, Humans, Neuroglia, Sequence Analysis, RNA, Microglia, Single-Cell Analysis, Transcriptome

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