Cited by 483Single-cell and spatial transcriptomicsGenomics and Chromatin DynamicsEpigenetics and DNA Methylation

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

Spatial epigenome–transcriptome co-profiling of mammalian tissues

Open access

Published version of a preprint

Abstract

Emerging spatial technologies, including spatial transcriptomics and spatial epigenomics, are becoming powerful tools for profiling of cellular states in the tissue context 1–5 . However, current methods capture only one layer of omics information at a time, precluding the possibility of examining the mechanistic relationship across the central dogma of molecular biology. Here, we present two technologies for spatially resolved, genome-wide, joint profiling of the epigenome and transcriptome by cosequencing chromatin accessibility and gene expression, or histone modifications (H3K27me3, H3K27ac or H3K4me3) and gene expression on the same tissue section at near-single-cell resolution. These were applied to embryonic and juvenile mouse brain, as well as adult human brain, to map how epigenetic mechanisms control transcriptional phenotype and cell dynamics in tissue. Although highly concordant tissue features were identified by either spatial epigenome or spatial transcriptome we also observed distinct patterns, suggesting their differential roles in defining cell states. Linking epigenome to transcriptome pixel by pixel allows the uncovering of new insights in spatial epigenetic priming, differentiation and gene regulation within the tissue architecture. These technologies are of great interest in life science and biomedical research.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper presents spatial ATAC–RNA-seq and spatial CUT&Tag–RNA-seq for genome-wide joint profiling of chromatin accessibility or histone modifications and mRNA on the same tissue section at near-single-cell resolution. These methods were applied to embryonic and juvenile mouse brain and adult human hippocampus, revealing spatial epigenetic–transcriptional relationships including epigenetic priming and region-specific gene regulation.

Key findings

  • In P22 mouse brain with 100 × 100 barcodes, spatial ATAC–RNA-seq yielded a median of 14,284 unique fragments per pixel, with 19% enriched at transcription start sites and 26% in peaks; RNA detected 22,914 genes, averaging 1,073 genes and 2,358 UMIs per pixel.
  • Spatial CUT&Tag–RNA-seq on P22 mouse brain gave median unique fragments per pixel of 10,644 for H3K27me3, 10,002 for H3K27ac and 2,507 for H3K4me3; paired RNA averaged 2,011, 1,513 and 1,329 genes per pixel, respectively.
  • Joint clustering of E13 mouse embryo ATAC and RNA identified a new neuronal cluster J10 not readily resolved by single modalities; in P22 mouse brain, 21,417 significant peak-to-gene linkages were detected by ATAC–RNA-seq and 19,468 from H3K27ac CUT&Tag–RNA-seq.
  • Region-specific regulation was observed: H3K27me3 anticorrelated with RNA in corpus callosum (for example, Mal, Mag and Car2 low CSS/high RNA; Grin2b and Syt1 high CSS/low RNA), while H3K27ac and H3K4me3 correlated with RNA; potential bivalency/poising was noted for genes such as Ptprz1 and Fnbp1. In human hippocampus, PROX1 was highly expressed in the granule cell layer but showed modest chromatin accessibility.

Methods

  • Frozen tissue sections were fixed, subjected to Tn5 transposition for ATAC or antibody-guided pA-Tn5 CUT&Tag, then in-tissue reverse transcription; sequential microfluidic barcode A and B ligation generated 50 × 50 or 100 × 100 pixel grids (2,500 or 10,000 pixels), with separate gDNA and cDNA libraries sequenced on Illumina NovaSeq 6000 (paired-end, 150 bp).
  • Samples included E13 mouse embryo, P21/22 mouse brain and adult human hippocampus from a 31-year-old male; pixel sizes were 50 μm for embryo and human tissue and 20 μm for juvenile mouse brain. Analyses used Signac, Seurat and ArchR, with label transfer to scRNA-seq, scATAC-seq and scCUT&Tag atlases, gene activity score (GAS), chromatin silencing score (CSS), peak-to-gene links, pseudotime and GO enrichment.
  • Reproducibility was high: Pearson correlations were r = 0.98 for ATAC and r = 0.98 for RNA in spatial ATAC–RNA-seq, and r = 0.96 for CUT&Tag (H3K27ac) and r = 0.89 for RNA in spatial CUT&Tag–RNA-seq.

Limitations

  • The paper states that low RNA expression despite open chromatin may be due in part to sequencing depth and RNA detection sensitivity, and interprets some such cases as lineage priming.
  • The method is near-single-cell, not single-cell: 50-μm pixels contained an average of 25 cells for E10 mouse embryo and 1–9 cells for human hippocampus, while 20-μm pixels contained 1–3 cells in juvenile mouse brain; when more than one cell type was present in a pixel, the major cell type was assigned.
  • Evident from the text: the human analysis used one 31-year-old male hippocampus sample, limiting generalizability, and the epigenome–transcriptome relationships are largely correlative rather than causally tested.

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

Integrity signals

Registry

Crossref · 10.1038/s41586-023-05795-1

Cites

65 of 65 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Design and evaluation of spatial epigenome–transcriptome cosequencing with E13 mouse embryo. a, Schematic workflow. b, Comparison of number of unique fragments and fraction of reads in peaks (FRiP) in spatial ATAC–RNA-seq and spatial CUT&Tag–RNA-seq. c, Gene and UMI count distribution in spatial ATAC–RNA-seq and spatial CUT&Tag–RNA-seq. Number of pixels in E13, 2,187; in human brain, 2,500; in mouse brain (ATAC), 9,215; in mouse brain (H3K27me3), 9,752; in mouse brain (H3K27ac), 9,370; in mouse brain (H3K4me3), 9,548. Box plots show the median (centre line), the first and third quartiles (box limits) and 1.5× interquartile range (whiskers). d, Spatial distribution and UMAP of all clusters for ATAC, RNA and joint clustering of ATAC and RNA data. Overlay of clusters with the tissue image shows that spatial clusters precisely match anatomic regions. Pixel size, 50 µm; scale bars, 1 mm. e, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial ATAC–RNA-seq. f, Pseudotime analysis from radial glia to postmitotic premature neurons visualized at the spatial level. g, Heatmaps delineating gene expression and GAS for marker genes. h, Dynamic changes in GAS and gene expression across pseudotime.
Fig. 1 Design and evaluation of spatial epigenome–transcriptome cosequencing with E13 mouse embryo. a, Schematic workflow. b, Comparison of number of unique fragments and fraction of reads in peaks (FRiP) in spatial ATAC–RNA-seq and spatial CUT&Tag–RNA-seq. c, Gene and UMI count distribution in spatial ATAC–RNA-seq and spatial CUT&Tag–RNA-seq. Number of pixels in E13, 2,187; in human brain, 2,500; in mouse brain (ATAC), 9,215; in mouse brain (H3K27me3), 9,752; in mouse brain (H3K27ac), 9,370; in mouse brain (H3K4me3), 9,548. Box plots show the median (centre line), the first and third quartiles (box limits) and 1.5× interquartile range (whiskers). d, Spatial distribution and UMAP of all clusters for ATAC, RNA and joint clustering of ATAC and RNA data. Overlay of clusters with the tissue image shows that spatial clusters precisely match anatomic regions. Pixel size, 50 µm; scale bars, 1 mm. e, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial ATAC–RNA-seq. f, Pseudotime analysis from radial glia to postmitotic premature neurons visualized at the spatial level. g, Heatmaps delineating gene expression and GAS for marker genes. h, Dynamic changes in GAS and gene expression across pseudotime.CC BY · PubMed Central
Fig. 2 Spatial chromatin accessibility and transcriptome co-profiling of P22 mouse brain. a, Design of microfluidic chips for 100 × 100 barcodes with 20-μm channel size. b, Spatial distribution and UMAP of all clusters for ATAC and RNA in spatial ATAC–RNA-seq of mouse brain. Pixel size, 20 µm; scale bars, 1 mm. c, Integration of ATAC data and scATAC-seq data30 from mouse brain. d, Integration of RNA data and scRNA-seq data32 from mouse brain. e, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial ATAC–RNA-seq. A list of abbreviation definitions can be found in Supplementary Table 1.
Fig. 2 Spatial chromatin accessibility and transcriptome co-profiling of P22 mouse brain. a, Design of microfluidic chips for 100 × 100 barcodes with 20-μm channel size. b, Spatial distribution and UMAP of all clusters for ATAC and RNA in spatial ATAC–RNA-seq of mouse brain. Pixel size, 20 µm; scale bars, 1 mm. c, Integration of ATAC data and scATAC-seq data30 from mouse brain. d, Integration of RNA data and scRNA-seq data32 from mouse brain. e, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial ATAC–RNA-seq. A list of abbreviation definitions can be found in Supplementary Table 1.CC BY · PubMed Central
Fig. 3 Spatial histone modification and transcriptome co-profiling of P22 mouse brain. a–c, Spatial distribution and UMAP of all clusters for H3K27me3 and RNA (a), H3K27ac and RNA (b) and H3K4me3 and RNA (c) in mouse brain. Pixel size, 20 µm; scale bars, 1 mm. d, Integration of H3K27me3 data with scCUT&Tag (H3K27me3) data37 from mouse brain. e, Integration of H3K27ac data with scCUT&Tag (H3K27ac) data37 from mouse brain. f, Integration of RNA data in spatial CUT&Tag (H3K27me3)–RNA-seq, spatial CUT&Tag (H3K27ac)–RNA-seq and spatial CUT&Tag (H3K4me3)–RNA-seq with scRNA-seq data32 from mouse brain.
Fig. 3 Spatial histone modification and transcriptome co-profiling of P22 mouse brain. a–c, Spatial distribution and UMAP of all clusters for H3K27me3 and RNA (a), H3K27ac and RNA (b) and H3K4me3 and RNA (c) in mouse brain. Pixel size, 20 µm; scale bars, 1 mm. d, Integration of H3K27me3 data with scCUT&Tag (H3K27me3) data37 from mouse brain. e, Integration of H3K27ac data with scCUT&Tag (H3K27ac) data37 from mouse brain. f, Integration of RNA data in spatial CUT&Tag (H3K27me3)–RNA-seq, spatial CUT&Tag (H3K27ac)–RNA-seq and spatial CUT&Tag (H3K4me3)–RNA-seq with scRNA-seq data32 from mouse brain.CC BY · PubMed Central
Fig. 4 Region-specific epigenetic regulation of gene expression. a–c, Correlation of H3K27me3 CSS and RNA gene expression (a), H3K27ac GAS and RNA gene expression (b) and H3K4me3 GAS and RNA gene expression (c) in corpus callosum. d, Upset plot of H3K27me3 CSS and RNA gene expression in striatum and deeper and superficial cortical layers; –, low CSS or gene expression; +, high CSS or gene expression. e, Venn diagrams showing high (+) or low (–) CSS/GAS for different histone modifications in corpus callosum with common RNA marker genes. Source Data
Fig. 4 Region-specific epigenetic regulation of gene expression. a–c, Correlation of H3K27me3 CSS and RNA gene expression (a), H3K27ac GAS and RNA gene expression (b) and H3K4me3 GAS and RNA gene expression (c) in corpus callosum. d, Upset plot of H3K27me3 CSS and RNA gene expression in striatum and deeper and superficial cortical layers; –, low CSS or gene expression; +, high CSS or gene expression. e, Venn diagrams showing high (+) or low (–) CSS/GAS for different histone modifications in corpus callosum with common RNA marker genes. Source DataCC BY · PubMed Central
Fig. 5 Spatial chromatin accessibility and transcriptome co-profiling of human hippocampus. a, Brightfield image, spatial distribution and UMAP of all clusters based on ATAC and RNA in the human hippocampus. ML, molecular layer; Pyr, pyramidal neurons. Pixel size, 50 μm, scale bars, 1 mm. b, Integration of our ATAC data with scATAC-seq data42 from human hippocampus. c, Integration of our RNA data with snRNA-seq data from human brain43. d, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA. Oligo, oligodendrocytes; astro, astrocytes; DG, dentate gyrus.
Fig. 5 Spatial chromatin accessibility and transcriptome co-profiling of human hippocampus. a, Brightfield image, spatial distribution and UMAP of all clusters based on ATAC and RNA in the human hippocampus. ML, molecular layer; Pyr, pyramidal neurons. Pixel size, 50 μm, scale bars, 1 mm. b, Integration of our ATAC data with scATAC-seq data42 from human hippocampus. c, Integration of our RNA data with snRNA-seq data from human brain43. d, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA. Oligo, oligodendrocytes; astro, astrocytes; DG, dentate gyrus.CC BY · PubMed Central
Extended Data Fig. 1 Workflow of spatial-ATAC-RNA-seq and spatial-CUT&Tag-RNA-seq. a, Schematic workflow of spatial-ATAC-RNA-seq and spatial-CUT&Tag-RNA-seq. b–c, Chemistry workflow of ATAC (b) and RNA (c) in spatial-ATAC-RNA-seq. d–e, Chemistry workflow of CUT&Tag (d) and RNA (e) in spatial-CUT&Tag-RNA-seq.
Extended Data Fig. 1 Workflow of spatial-ATAC-RNA-seq and spatial-CUT&Tag-RNA-seq. a, Schematic workflow of spatial-ATAC-RNA-seq and spatial-CUT&Tag-RNA-seq. b–c, Chemistry workflow of ATAC (b) and RNA (c) in spatial-ATAC-RNA-seq. d–e, Chemistry workflow of CUT&Tag (d) and RNA (e) in spatial-CUT&Tag-RNA-seq.CC BY · PubMed Central
Extended Data Fig. 2 Further analysis of spatial-ATAC-RNA-seq for E13 mouse embryo. a, H&E image from an adjacent tissue section of E13 mouse embryo. b, Spatial mapping of GAS and gene expression for selected marker genes in spatial-ATAC-RNA-seq. c, Integration of scRNA-seq data12 from E13.5 mouse embryos with ATAC and RNA data in spatial ATAC-RNA-seq. MOCA, Mouse Organogenesis Cell Atlas. d, Spatial mapping of cell types identified by label transfer from scRNA-seq12 to ATAC (top) and RNA (bottom). e, Genome track visualization of marker genes with peak-to-gene links for distal regulatory elements and peak co-accessibility. f, The expression level and the percentage of pixels in all clusters (marker genes for each cluster) for RNA data in spatial-ATAC-RNA-seq. g, Dot plot showing the identification of positive TF regulators. h, Spatial mapping of deviation scores for selected TF motifs. i, Annotation of marker peaks in different clusters. j, GREAT enrichment analysis of marker peaks in different clusters (Binomial and hypergeometric tests).
Extended Data Fig. 2 Further analysis of spatial-ATAC-RNA-seq for E13 mouse embryo. a, H&E image from an adjacent tissue section of E13 mouse embryo. b, Spatial mapping of GAS and gene expression for selected marker genes in spatial-ATAC-RNA-seq. c, Integration of scRNA-seq data12 from E13.5 mouse embryos with ATAC and RNA data in spatial ATAC-RNA-seq. MOCA, Mouse Organogenesis Cell Atlas. d, Spatial mapping of cell types identified by label transfer from scRNA-seq12 to ATAC (top) and RNA (bottom). e, Genome track visualization of marker genes with peak-to-gene links for distal regulatory elements and peak co-accessibility. f, The expression level and the percentage of pixels in all clusters (marker genes for each cluster) for RNA data in spatial-ATAC-RNA-seq. g, Dot plot showing the identification of positive TF regulators. h, Spatial mapping of deviation scores for selected TF motifs. i, Annotation of marker peaks in different clusters. j, GREAT enrichment analysis of marker peaks in different clusters (Binomial and hypergeometric tests).CC BY · PubMed Central
Extended Data Fig. 3 Further pseudotime analysis of radial glia and postmitotic premature neurons in spatial-ATAC-RNA-seq. a, GO enrichment analysis for genes from Fig. 1g. b,c, Pseudotime analysis from radial glia to postmitotic premature neurons with GAS (b) and gene expression (c). d, Monocle2 analyses showing different states in (b). e, Heatmap of different states along the pseudotime trajectory. f, GO analysis of genes in red box of (e) (One-sided version of Fisher’s exact test, p-value was adjusted for multiple comparisons by Benjamini & Hochberg method).
Extended Data Fig. 3 Further pseudotime analysis of radial glia and postmitotic premature neurons in spatial-ATAC-RNA-seq. a, GO enrichment analysis for genes from Fig. 1g. b,c, Pseudotime analysis from radial glia to postmitotic premature neurons with GAS (b) and gene expression (c). d, Monocle2 analyses showing different states in (b). e, Heatmap of different states along the pseudotime trajectory. f, GO analysis of genes in red box of (e) (One-sided version of Fisher’s exact test, p-value was adjusted for multiple comparisons by Benjamini & Hochberg method).CC BY · PubMed Central
Extended Data Fig. 4 Further analysis of spatial-ATAC-RNA-seq for P22 mouse brain. a, Spatial mapping of gene activity scores and gene expression for selected marker genes in spatial-ATAC-RNA-seq. b, Genome track visualization of marker genes with peak-to-gene links for distal regulatory elements and peak co-accessibility. c, Spatial mapping of cell types identified by label transfer from scATAC-seq30 to ATAC data. IT: intratelencephalic. PT: pyramidal tract. NP: near-projecting. CT: corticothalamic. L: layer. d, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA data.
Extended Data Fig. 4 Further analysis of spatial-ATAC-RNA-seq for P22 mouse brain. a, Spatial mapping of gene activity scores and gene expression for selected marker genes in spatial-ATAC-RNA-seq. b, Genome track visualization of marker genes with peak-to-gene links for distal regulatory elements and peak co-accessibility. c, Spatial mapping of cell types identified by label transfer from scATAC-seq30 to ATAC data. IT: intratelencephalic. PT: pyramidal tract. NP: near-projecting. CT: corticothalamic. L: layer. d, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA data.CC BY · PubMed Central
Extended Data Fig. 5 Further analysis of P22 mouse brain in spatial-ATAC-RNA-seq. a, Candidate TF regulators of Sox2, Pax6, and Mobp. Highlighted points are TFs with abs(regulation score) ≥ 1 (−log10 scale)33, with all other TFs shown in gray (Z test). b, Spatial mapping of deviation scores for selected TF motifs. c, Heatmaps of peak-to-gene links in spatial-ATAC-RNA-seq. d, The number of significantly correlated peaks for each gene.
Extended Data Fig. 5 Further analysis of P22 mouse brain in spatial-ATAC-RNA-seq. a, Candidate TF regulators of Sox2, Pax6, and Mobp. Highlighted points are TFs with abs(regulation score) ≥ 1 (−log10 scale)33, with all other TFs shown in gray (Z test). b, Spatial mapping of deviation scores for selected TF motifs. c, Heatmaps of peak-to-gene links in spatial-ATAC-RNA-seq. d, The number of significantly correlated peaks for each gene.CC BY · PubMed Central
Extended Data Fig. 6 Spatial chromatin accessibility and transcriptome co-sequencing of P21mouse brain. a—c, Spatial distribution and UMAP of all the clusters for ATAC (a), RNA (b), and joint clustering of ATAC and RNA (c) in spatial-ATAC-RNA-seq for the mouse brain. Pixel size, 20 µm. Scale bar, 1 mm. d, Nissl-stained image from an adjacent tissue section of P21 mouse brain. Scale bar, 1 mm. e, Integration of ATAC data in spatial-ATAC-RNA-seq with scATAC-seq data30 from mouse brain. f, Integration of RNA data in spatial-ATAC-RNA-seq with scRNA-seq data32 from mouse brain. g, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial-ATAC-RNA-seq.
Extended Data Fig. 6 Spatial chromatin accessibility and transcriptome co-sequencing of P21mouse brain. a—c, Spatial distribution and UMAP of all the clusters for ATAC (a), RNA (b), and joint clustering of ATAC and RNA (c) in spatial-ATAC-RNA-seq for the mouse brain. Pixel size, 20 µm. Scale bar, 1 mm. d, Nissl-stained image from an adjacent tissue section of P21 mouse brain. Scale bar, 1 mm. e, Integration of ATAC data in spatial-ATAC-RNA-seq with scATAC-seq data30 from mouse brain. f, Integration of RNA data in spatial-ATAC-RNA-seq with scRNA-seq data32 from mouse brain. g, Spatial mapping of GAS and gene expression for selected marker genes in different clusters for ATAC and RNA in spatial-ATAC-RNA-seq.CC BY · PubMed Central
Extended Data Fig. 7 Further analysis for spatial-CUT&Tag-RNA-seq with P22 mouse brain. a, Integration of CUT&Tag (H3K4me3) data in spatial-CUT&Tag-RNA-seq with scCUT&Tag (H3K4me3) data36 from mouse brain. b, Integration of RNA data in spatial- CUT&Tag(H3K27me3)-RNA-seq, spatial-CUT&Tag(H3K27ac)-RNA-seq, and spatial-CUT&Tag(H3K4me3)-RNA-seq with scRNA-seq data32 from mouse brain. c-e, Integration of RNA data in spatial-CUT&Tag(H3K27me3)-RNA-seq (c), RNA data in data in spatial-CUT&Tag(H3K27ac)-RNA-seq (d), and RNA data in spatial-CUT&Tag(H3K4me3)-RNA-seq (e) with scRNA-seq data32 from mouse brain. f, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA data in spatial-CUT&Tag(H3K27me3)-RNA-seq. g, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA (top) and from scRNA-seq32 to CUT&Tag (H3K27ac, bottom) data in spatial-CUT&Tag(H3K27ac)-RNA-seq. h, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA (top) and from scRNA-seq32 to CUT&Tag (H3K4me3, bottom) data in spatial-CUT&Tag(H3K4me3)-RNA-seq.
Extended Data Fig. 7 Further analysis for spatial-CUT&Tag-RNA-seq with P22 mouse brain. a, Integration of CUT&Tag (H3K4me3) data in spatial-CUT&Tag-RNA-seq with scCUT&Tag (H3K4me3) data36 from mouse brain. b, Integration of RNA data in spatial- CUT&Tag(H3K27me3)-RNA-seq, spatial-CUT&Tag(H3K27ac)-RNA-seq, and spatial-CUT&Tag(H3K4me3)-RNA-seq with scRNA-seq data32 from mouse brain. c-e, Integration of RNA data in spatial-CUT&Tag(H3K27me3)-RNA-seq (c), RNA data in data in spatial-CUT&Tag(H3K27ac)-RNA-seq (d), and RNA data in spatial-CUT&Tag(H3K4me3)-RNA-seq (e) with scRNA-seq data32 from mouse brain. f, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA data in spatial-CUT&Tag(H3K27me3)-RNA-seq. g, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA (top) and from scRNA-seq32 to CUT&Tag (H3K27ac, bottom) data in spatial-CUT&Tag(H3K27ac)-RNA-seq. h, Spatial mapping of cell types identified by label transfer from scRNA-seq32 to RNA (top) and from scRNA-seq32 to CUT&Tag (H3K4me3, bottom) data in spatial-CUT&Tag(H3K4me3)-RNA-seq.CC BY · PubMed Central
Extended Data Fig. 8 Spatial epigenome and transcriptome co-sequencing and integrative analysis of P21 mouse brain. a-c, Spatial distribution and UMAP of all the clusters for CUT&Tag (H3K27ac) (a), RNA (b), and joint clustering of CUT&Tag (H3K27ac) and RNA (c) in spatial-CUT&Tag-RNA-seq for the mouse brain. Pixel size, 20 µm. Scale bar, 1 mm. d, Nissl-stained image from an adjacent tissue section of P21 mouse brain. Scale bar, 1 mm. e, Integration of CUT&Tag (H3K27ac) data in spatial-CUT&Tag-RNA-seq with scCUT&Tag (H3K27ac) data37 from mouse brain. f, Integration of RNA data in spatial-CUT&Tag-RNA-seq with scRNA-seq data32 from mouse brain.
Extended Data Fig. 8 Spatial epigenome and transcriptome co-sequencing and integrative analysis of P21 mouse brain. a-c, Spatial distribution and UMAP of all the clusters for CUT&Tag (H3K27ac) (a), RNA (b), and joint clustering of CUT&Tag (H3K27ac) and RNA (c) in spatial-CUT&Tag-RNA-seq for the mouse brain. Pixel size, 20 µm. Scale bar, 1 mm. d, Nissl-stained image from an adjacent tissue section of P21 mouse brain. Scale bar, 1 mm. e, Integration of CUT&Tag (H3K27ac) data in spatial-CUT&Tag-RNA-seq with scCUT&Tag (H3K27ac) data37 from mouse brain. f, Integration of RNA data in spatial-CUT&Tag-RNA-seq with scRNA-seq data32 from mouse brain.CC BY · PubMed Central
Extended Data Fig. 9 Spatial mapping of CSS, GAS, and gene expression of selected genes for P22 mouse brain. a-g, Spatial mapping of CSS, GAS, and gene expression of Mal (a), Mag (b), Car2 (c), Grin2b (d), Syt1 (e), Gpr88 (f), and Ptprz1 (g) from ATAC and RNA in spatial-ATAC-RNA-seq, and CUT&Tag (H3K27me3, H3K27ac, or H3K4me3) and RNA in spatial-CUT&Tag-RNA-seq. h, Spatial mapping of CSS and gene expression of Nav3, Sncb, Ablim2, and Gng7 in spatial-CUT&Tag(H3K27me3)-RNA-seq.
Extended Data Fig. 9 Spatial mapping of CSS, GAS, and gene expression of selected genes for P22 mouse brain. a-g, Spatial mapping of CSS, GAS, and gene expression of Mal (a), Mag (b), Car2 (c), Grin2b (d), Syt1 (e), Gpr88 (f), and Ptprz1 (g) from ATAC and RNA in spatial-ATAC-RNA-seq, and CUT&Tag (H3K27me3, H3K27ac, or H3K4me3) and RNA in spatial-CUT&Tag-RNA-seq. h, Spatial mapping of CSS and gene expression of Nav3, Sncb, Ablim2, and Gng7 in spatial-CUT&Tag(H3K27me3)-RNA-seq.CC BY · PubMed Central
Extended Data Fig. 10 Further analysis of human hippocampus in spatial-ATAC-RNA-seq. a, Spatial mapping of cell types identified by label transfer from scRNA-seq43 to RNA data in spatial-ATAC-RNA-seq for human hippocampus. b, Dot plot showing the identification of positive TF regulators. c, Spatial mapping of deviation scores for selected TF motifs in spatial-ATAC-RNA-seq.
Extended Data Fig. 10 Further analysis of human hippocampus in spatial-ATAC-RNA-seq. a, Spatial mapping of cell types identified by label transfer from scRNA-seq43 to RNA data in spatial-ATAC-RNA-seq for human hippocampus. b, Dot plot showing the identification of positive TF regulators. c, Spatial mapping of deviation scores for selected TF motifs in spatial-ATAC-RNA-seq.CC BY · PubMed Central

Verification

Claims the paper makes

70 Own result · 25 Not checked yet · 22 Own interpretation · 13 Own method · 4 No basis given · 1 Partly supported · accuracy not yet measured

Verify all 45 references — free account

Show all 79 claims

  • Probably not a substantial claim. 13
  • It may state two separate claims as one. 46
56 of 56
  • The paper relies on this without support.

    The selected tissues (E13 mouse embryo, P21/22 mouse brain, adult human hippocampus) are representative models for studying spatial epigenome–transcriptome relationships in mammalian tissues.

    “We applied these techniques to embryonic and juvenile mouse brain, as well as adult human brain hippocampus, to dissect epigenetic and transcriptional states and their role in the regulation of cell types dynamics in tissue.”
  • The paper relies on this without support.

    As of the time of the study, current technologies had not been able to perform unbiased genome-wide comapping of epigenome and transcriptome on the same tissue section at the cellular level.

    “As of today, current technologies have not been able to perform unbiased genome-wide comapping of epigenome and transcriptome on the same tissue section at the cellular level.”
  • The paper relies on this without support.

    To comprehensively understand the mechanism of gene regulation, different layers of omics information need to be profiled simultaneously.

    “However, to comprehensively understand the mechanism of gene regulation, different layers of omics information need to be profiled simultaneously.”
3 of 3

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 58 references

Cited by

Sorted by citations

Items 201–300 of 482

Metadata

Authors
Zhang, D., Deng, Y., Kukanja, P., Agirre, E., Bartosovic, M., Dong, M., Ma, C., Ma, S., Su, G., Bao, S., Liu, Y., Xiao, Y., Rosoklija, G. B., Dwork, A. J., Mann, J. J., Leong, K. W., Boldrini, M., Wang, L., Haeussler, M., Raphael, B. J., Kluger, Y., Castelo-Branco, G., & Fan, R.
Affiliations
  • Yale University (US)
  • Yale Cancer Center (US)
  • Karolinska Institutet (SE)
  • Princeton University (US)
  • Broad Institute (US)
  • Columbia University (US)
  • Macedonian Academy of Sciences and Arts (MK)
  • New York State Psychiatric Institute (US)
  • Columbia University Irving Medical Center (US)
  • University of California, Santa Cruz (US)
Corresponding author
  • Deng, Yanxiang
  • Fan, Rong
Venue
Nature
Published
2023-03-15
Volume
616
Issue
7955
Pages
113-122
ISSN
  • 0028-0836 (ISSN-L)
  • 1476-4687
OpenAlex
W4324344255
Cited by
483
References
58
Citations by year
2026: 172 · 2025: 163 · 2024: 87 · 2023: 59
FWCI
42.73
Type
journal article
Language
en
Open access
Locations
  • Nature (publishedVersion, cc-by)
  • Research Square (acceptedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
Funders
  • Yale University
  • Li Ka Shing Foundation
  • Society for Immunotherapy of Cancer
  • Göran Gustafssons Stiftelser
  • European Commission (HORIZON2020, 794689)
  • Cancerfonden (190394 Pj, 2019-)
  • Svenska Sällskapet för Medicinsk Forskning (JUB2019)
  • Karolinska Institutet
  • Knut och Alice Wallenbergs Stiftelse (2019–0107)
  • Vetenskapsrådet (2019-01360)
  • National Institutes of Health (1u01ca294514-01, 1rf1mh128876-01)
  • National Institute on Aging (U54AG076043, U54AG079759)
  • National Institute of Mental Health (RF1 MH128876, RM1MH132648)
  • National Cancer Institute (U01 CA294514, UH3CA257393, UG3CA257393, U24CA248453, R01CA245313, U54CA274509)
  • National Institute of Diabetes and Digestive and Kidney Diseases (U54DK106857)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
Keywords
spatial epigenomics, spatial transcriptomics, chromatin accessibility, histone modifications, H3K27me3, H3K27ac, H3K4me3, single-cell resolution, mouse brain, human brain, cell dynamics, gene expression regulation, epigenetic priming, tissue structure
MeSH terms
Epigenome, Aging, Animals, Brain, Chromatin, Gene Expression Regulation, Histones, Humans, Mammals, Organ Specificity, Gene Expression Profiling, Epigenesis, Genetic, Mice, Epigenomics, 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-07, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-07 (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