Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: Top panel shows an iterative optimization scheme to estimate geometric transformation ( φ ) and latent feature distribution ( π ) by minimizing the normed difference (error) between the geometric and feature-transformed CCFv3 section to target MERFISH. Middle panel illustrates the application of estimated geometric transformation ( φ ) to deform the CCFv3 atlas to MERFISH coordinates (left); the application of latent feature distribution ( π ) to generate gene distributions on initial CCFv3 geometry (middle); and the application of inverse geometric transformation ( φ −1 ) to deform MERFISH genes to CCFv3 coordinates (right). Gene with the highest probability of expression at each location is shown as a MERFISH feature. Bottom panel illustrates the results of mapping CCFv3 sections to corresponding MERFISH sections. a , f 10 μm atlas sections at Z = 385 and Z = 485 out of 1320 visually chosen to match MERFISH architecture ( e , j ) rendered as meshes at 100 μm. e , j MERFISH sections rendered as meshes at 50 μm, with mRNA density depicted as a feature. b , g Geometric mappings ( φ ) of CCFv3 sections to MERFISH coordinates with the approximate determinant of the Jacobian showing areas of contraction (blue) and expansion (red). c , h Estimated mRNA density per atlas region ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${w}_{i}^{{\prime} }$$\end{document} w i ′ in ), as given by π shown in CCFv3 coordinates. d , i Estimated mRNA density shown per atlas region following geometric deformation to MERFISH coordinates.

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques: Transformation Assay, Expressing

Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: a – c Relative expression on target section of MERFISH transcriptomics data for three genes ( Gfap ( a ), Trp53i11 ( b ), Wipf3 ( c )) out of a set of twenty (shown at the right), with demonstrated spatial variability according to a computed mutual information score. e – g Predicted expression for the same three genes ( Gfap ( e ), Trp53i11 ( f ), Wipf3 ( g )) in each region of the CCFv3 section Z = 485 out of 1320, as part of the latent distribution over genes estimated in tandem with a geometric transformation to align the CCFv3 section to MERFISH section. d Gene with the highest probability in MERFISH target section. h Predicted gene with the highest probability in estimated latent distribution for the CCFv3 section.

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques: Expressing, Transformation Assay

Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: a Top row depicts MERFISH target rendered as a mesh with cell density (left) and cell type with the highest probability (right) used to summarize cell type distributions. Bottom shows predicted cell density (left) and cell type with the highest probability (right) for the latent feature distribution estimated for each CCFv3 region in the native CCFv3 coordinates. b , c Top row depicts the probability of expression for each gene out of a subset of 6 selected from a total measured set of ~500 as those with high spatial variance. Bottom row depicts the estimated probability of expression for each gene for the latent feature distribution estimated for each CCFv3 region in the native CCFv3 coordinates.

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques: Expressing

Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: a Variance in estimated cell subtype probabilities per CCFv3 region across three replicates summed overall CCFv3 regions. b Spatial variance of cell subtype probabilities, estimated empirically from three pulled-back MERFISH sections, per CCFv3 region and summed across all regions. c – e Probability of excitatory neuron subtype 2 (star in a ) for each of the three mice in CCFv3 coordinates. Yellow arrow highlights the area of the dentate gyrus with differences in excitatory neuron subtype 2 probabilities. f – h Probability for astrocyte subtypes 1,2, and 3 (stars in b ) in empirical distribution computed from all three mice in CCFv3 coordinates (most likely cell type shown in Fig. h). Yellow arrow highlights area of CA1 with differences in astrocyte probability medially to laterally in subtypes 1 and 2 but not 3. Orange arrow highlights differences medially to laterally and left and right in areas of the pons in astrocyte probability for subtypes 1 and 2 but not 3. Black lines indicate boundaries between CCFv3 regions.

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques:

Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: a , d Original CCFv3 and DevCCF ontologies at location Z = 680 out of 1320. b CCFv3 geometry with predicted DevCCF atlas ontology. Delineations of original CCFv3 partitions are outlined in gray. e DevCCF atlas geometry with predicted CCFv3 ontology. Delineations of original DevCCF partitions are outlined in gray. c , f Entropy of predicted ontologies, with higher entropy values (light) indicating less 1:1 correspondence between ontologies. g Predicted cell type with highest probability per simplex in CCFv3 atlas following mapping to MERFISH target (shown in Fig. with xIV-LDDMM. j Predicted cell type with a highest probability per simplex in DevCCF atlas following mapping to same MERFISH target with xIV-LDDMM. h , k Estimated geometric transformation, φ 1 , in each setting applied to each atlas, with areas of expansion (red) and contraction (blue) as measured by the determinant of the Jacobian. White arrow highlights differences in ontologies in amygdala and striatum designation leading to different geometric transformations. i, l Entropy of estimated cell type distribution per simplex in atlas. Circled area of the hippocampus highlights differences in atlas ontologies leading to differences in the estimated entropy of cell type distributions.

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques: Transformation Assay

Journal: Nature Communications

Article Title: Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

doi: 10.1038/s41467-024-47883-4

Figure Lengend Snippet: a Original DAPI-stained image, digitized at 2.5 μm resolution for tissue section measured with MERFISH technology (Fig. . c Image-varifold particle representation (black points) of DAPI-stained image overlaying the corresponding CCFv3 section in their respective initial coordinate spaces. Thresholded foreground pixels from ( a ) converted to particle image-varifold representation over a feature space of ~30 binned grayscale values. d Alignment of CCFv3 section to DAPI particles following diffeomorphic transformation to the DAPI coordinate space. b Alignment of CCFv3 section and DAPI particles in image format, with the deformed CCFv3 section image generated by resampling the deformed CCFv3 particles onto a regular 2.5 μm grid. White arrows highlight areas of alignment in the area of the substantia inominata (SI) and layer 1 of the cortex whereas red arrows highlight areas of questionable alignment in the area of the olfactory tubercle (OT).

Article Snippet: We demonstrate the efficacy of xIV-LDDMM for mapping tissue-scale atlases to cellular-scale data in mapping CCFv3 section Z = 675 to a section of cell-segmented MERFISH transcriptional data (courtesy of the JEFworks Lab, Johns Hopkins University) (Fig. ).

Techniques: Staining, Transformation Assay, Generated

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: Images of nuclei (DAPI), total polyA mRNA, and two MERFISH bits were obtained using epifluorescence (Epi) and spinning-disk confocal microscopy respectively. Both epifluorescence and confocal images were taken with 1 s exposure time. The MERFISH bit-1 and bit-2 images were high-pass filtered to remove cellular background.

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: Confocal Microscopy

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) A single-bit high-pass-filtered MERFISH confocal image of 242 genes in a brain tissue section taken with an exposure time of 0.1 s (left) and a magnified view of a single cell marked by the white box in the left image (right). ( b ) The correlation between the copy number of individual genes detected per field of view (FOV) using 0.1 s exposure time and those obtained using 1 s exposure time. The median ratio of the copy number and the Pearson correlation coefficient r are shown. The copy number per gene detected using 0.1 s exposure time is 24% of that detected using 1 s exposure time. ( c ) The same image as in ( a ) but after enhancement of signal-to-noise ratio (SNR) by a DL algorithm. ( d ) The same as ( b ) but after DL was used to enhance the SNR of the 0.1 s images. The copy number per gene detected with 0.1 s exposure time after DL-based enhancement is 89% of that detected using 1 s exposure time. Figure 1—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques:

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) Number of RNA molecules detected per field of view (FOV) at a single z-plane at the tissue depths of 10 µm and 90 µm in the first bit of the 242-gene MERFISH measurements in a 100-µm-thick section of the mouse cortex. ( b ) Logarithmic distribution of integrated photon counts of individual RNA molecules at the tissue depths of 10 µm and 90 µm identified in ( a ). In each boxplot, the midline represents the median value, the box represents the interquartile range (IQR), the lower whisker represents the smaller of the minimum data point or 1.5× IQR below the 25th percentile, and the upper whisker represents the greater of the maximum data point or 1.5× IQR above the 75th percentile. Molecule number (n, from left to right): 7565, 7914. ( c ) Number of RNA molecules detected per FOV at a single z-plane at tissue depths of 10 µm and 190 µm of in the first bit of the 156-gene MERFISH measurements in a 200-µm-thick section of the mouse hypothalamus. ( d ) Logarithmic distribution of integrated photon counts of individual RNA molecules at the tissue depths of 10 µm and 190 µm identified in ( c ). In each boxplot, the midline represents the median value, the box represents the IQR, the lower whisker represents the smaller of the minimum data point or 1.5× IQR below the 25th percentile, and the upper whisker represents the greater of the maximum data point or 1.5× IQR above the 75th percentile. Molecule number (n, from left to right): 9611, 7238. Figure 1—figure supplement 2—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: Whisker Assay

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) Example high-pass-filtered bit-1 images of a 242-gene MERFISH measurement in a 100-µm-thick section of mouse cortex stained with different concentrations of encoding probes. The concentration values refer to the concentration of each individual encoding probe. ( b ) Distribution of integrated photon counts of individual RNA molecules identified at different encoding probe concentrations. In each boxplot, the midline represents the median value, the box represents the interquartile range (IQR), the lower whisker represents the smaller of the minimum data point or 1.5× IQR below the 25th percentile, and the upper whisker represents the greater of the maximum data point or 1.5× IQR above the 75th percentile. Molecule number (n, from left to right): 31606, 29190, 83644. The signals from individual RNA molecules increased with the encoding probe concentration and reached saturation at ~1.0 nM per probe. We thus used 1 nM encoding probe concentrations for staining thick-tissue samples. ( c ) A 100-µm-thick mouse brain slice was stained with the 242-gene MERFISH encoding probes, followed by sequential hybridization with readout probes corresponding to the first, second, third, and fourth bit of the barcodes, each bit using a different readout probe concentration. High-pass-filtered bit-1, bit-2, bit-3, and bit-4 MERFISH images (each with a different concentration of readout probes) are shown. ( d ) Distribution of integrated photon counts of individual RNA molecules identified at different readout probe concentrations. In each boxplot, the midline represents the median value, the box represents the IQR, the lower whisker represents the smaller of the minimum data point or 1.5× IQR below the 25th percentile, and the upper whisker represents the greater of the maximum data point or 1.5× IQR above the 75th percentile. Molecule number (n, from left to right): 1939, 3815, 3869, 5248. The signal increased with readout probe concentration, but the background also increased when the probe concentration reached beyond 5 nM. We thus used 5 nM readout probe concentration for thick-tissue imaging. In addition to the probe concentrations, we also optimized readout probe incubation time. ( e ) The number of RNA molecules per field of view per z-plane and the normalized intensity of individual molecules at different tissue depths. The encoding probe concentration was 1 nM per encoding probe, the readout probe concentration was 5 nM, and the readout probe incubation time was 25 min for these measurements. Figure 2—figure supplement 1—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: Staining, Concentration Assay, Whisker Assay, Slice Preparation, Hybridization, Imaging, Incubation

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) Total copy number of decoded RNAs detected per field of view (FOV) per z-plane at different tissue depths in a 242-gene MERFISH measurement of the 100-µm-thick section. ( b ) Pearson correlation coefficients of RNA copy number of individual genes per FOV per z-plane detected at different tissue depths by MERFISH with the FPKM values measured by bulk RNA-seq. ( c ) Correlation of RNA copy number of individual genes per FOV per z-plane detected in the entire 100-µm-thick section by MERFISH with the FPKM values obtained by bulk RNA-seq. The Pearson correlation coefficient r is shown. ( d ) Example images of gel-embedded beads acquired in two rounds of imaging. Buffer exchanges were performed between imaging rounds, mimicking the MERFISH protocol. Because the gel expanded to a different extent in different rounds, the positions of beads changed from round to round in x, y, and z directions. Circles mark beads identified in both imaging rounds. Because the gel size changed, the x and y positions of the beads changed, and the brightness of these beads also changed due to the shift in their z positions. Arrows highlight beads detected in one of the imaging rounds, but not the other, due to the gel-size change, which move these beads out of focus. Figure 2—figure supplement 2—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: RNA Sequencing Assay, Imaging

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) Quantification of gel expansion factor in various buffers used in the MERFISH protocol. The initial gel size was the same as the coverslip, and the expansion factor after buffer exchange was determined as the ratio between the gel size after buffer exchange and the coverslip size. ( b ) In each round of MERFISH imaging, the sample is incubated with the readout probes in a wash buffer (containing either 10% ethylene carbonate [EC] or 10% formamide) for a duration of 15 min. Subsequently, the sample was rinsed with the wash buffer (without readout probes) to remove any excessive readout probes, followed by a treatment with imaging buffer containing glucose-oxidase-based or protocatechuate 3,4-dioxygenase rPCO-based oxygen scavenger system to reduce photobleaching. After the imaging process, the sample is treated with tris(2-carboxyethyl) phosphine buffer to cleave off the fluorescent dye linked to the readout probe through a disulfide bond, and finally washed with a solution of 2× saline-sodium citrate (SSC). All buffers, including wash, imaging, and cleavage buffers, contained 2× SSC. Gel-expansion factor in these buffers used in the MERFISH protocol was quantified and shown here. Reagents marked by * were selected for final use in the 3D thick-tissue MERFISH experiment. The dashed line highlights the expansion factor in the 2× SSC buffer alone. ( c ) XZ projection images of fiducial beads embedded in a gel undergoing buffer exchange for the indicated time period. Wash buffer containing 15% EC in 2× SSC causes noticeable gel distortion, which was recovered after 15 min in 2× SSC buffer without EC. Figure 2—figure supplement 3—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: Buffer Exchange, Imaging, Incubation, Saline

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) 3D images of DAPI and total polyA mRNA from a single field of view (FOV) in a 100-µm-thick mouse brain tissue slice (top), alongside a single z-plane at tissue depth of 50 µm marked by the yellow box in the top image (bottom). ( b ) Maximum-projection images of 10 consecutive 1 µm z-planes of individual MERFISH bits of the region marked in yellow box in the bottom panel of ( a ). ( c ) RNA molecules identified in the same region as in ( b ) with RNA molecules color coded by their genetic identities. ( d ) The RNA copy number for individual genes per unit area (100 2 µm 2 ) per z-plane detected in the 100 µm MERFISH measurements of mouse cortex versus the FPKM from bulk RNA-seq. The Pearson correlation coefficient r is shown. ( e ) The Pearson correlation between RNA copy number for individual genes per z-plane at different tissue depths detected by MERFISH and the FPKM values of individual genes from bulk RNA-seq. ( f ) Number of detected RNA molecules per FOV at different tissue depths. ( g ) The RNA copy number for individual genes per unit area (100 2 µm 2 ) per z-plane detected in the 100 µm MERFISH measurements of mouse brain sections in this work versus that detected by 10-µm-thick-tissue MERFISH measurements using an epifluorescence setup . The Pearson correlation coefficient r is shown. ( h ) The RNA copy number of individual genes per cell detected in the 100-µm-thick-tissue section versus that in individual 10-µm-thick z-ranges of the same sample. The 100-µm-thick section was evenly divided into ten 10 µm z-ranges to determine the latter. Data in this figure were generated from a 100-μm-thick section of the cortical region collected from an adult mouse. Figure 2—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: RNA Sequencing Assay, Generated

Journal: eLife

Article Title: Three-dimensional single-cell transcriptome imaging of thick tissues

doi: 10.7554/eLife.90029

Figure Lengend Snippet: ( a ) UMAP visualization of subclasses of cells identified in a 100-μm-thick section of the mouse cortex. Cells are color coded by subclass identities. IT: intratelencephalic projection neurons; ET: extratelencephalic projection neurons; CT: cortical-thalamic projection neurons; NP: near projection neurons; OPC: oligodendrocyte progenitor cells; Oligo: oligodendrocytes; Micro: microglia; Astro: astrocytes; VLMC: vascular leptomeningeal cells; Endo: endothelial cells; Peri: pericytes; SMC: smooth muscle cells; PVM: perivascular macrophages. ( b ) 3D spatial maps of the identified subclasses of excitatory neurons (left), inhibitory neurons (middle), and non-neuronal cells (right) within the 100 μm mouse cortex section. ( c ) UMAP visualization of major cell types identified in a 200-μm-thick section of the mouse anterior hypothalamus. Cells are color coded by cell type identities. ( d ) 3D spatial maps of the excitatory neuronal (left), inhibitory neuronal (middle), and non-neuronal (right) cell clusters identified in the 200-μm-thick mouse hypothalamus section. Cells are color coded by cell cluster identities in the top panels and two example excitatory neuronal (left), inhibitory neuronal (middle), and non-neuronal (right) cell clusters are shown in the bottom panels. ( e ) Boxplots showing the distributions of the nearest-neighbor distances from cells in individual inhibitory neuronal subclasses to cells in the same subclass (‘to self’) or other subclasses (‘to other’) in the mouse cortex obtained from the thick-tissue 3D MERFISH data. Cell numbers (n, from left to right): 44, 274, 125, 1161, 136, 797, 26, 330. *FDR <0.01 was determined with the Wilcoxon rank-sum one-sided test and adjusted to FDR by the Benjamini and Hochberg procedure. Only inhibitory neuronal clusters with at least 20 ‘self-self’ interacting pairs were examined and plotted. In each boxplot, the midline represents the median value, the box represents the interquartile range (IQR), the lower whisker represents the smaller of the minimum data point or 1.5× IQR below the 25th percentile, and the upper whisker represents the greater of the maximum data point or 1.5× IQR above the 75th percentile. ( f ) Distributions of the nearest-neighbor distances among all interneurons derived from the thick-tissue 3D MERFISH data of the mouse cortex. The distribution is fitted with a bimodal distribution (blue curve) with the two individual Gaussian peaks shown in red and green. ( g, h ) Same as ( e, f ) but for inhibitory neurons in the mouse anterior hypothalamus obtained from the thick-tissue 3D MERFISH data. Cell numbers in g (n, from left to right): 1404, 1222, 198, 567, 845, 1315, 866, 1263, 1149, 873, 455, 1517, 387, 1490, 409, 1321, 277, 1151, 773, 1871, 379, 830, 174, 1013, 598, 449, 811, 158, 194, 728, 103, 771, 119, 667, 217, 580, 98, 704, 1095, 1339, 186, 362, 82, 438, 41, 302. Data in this figure were generated from a 200-μm-thick section of the hypothalamic region collected from an adult mouse. Figure 3—source data 1. This source data file contains source data for .

Article Snippet: MERFISH readout probes , conjugated to either Cy5, Cy3B, or Alexa488 dye molecules through a disulfide linkage, were purchased from Bio-Synthesis, Inc.

Techniques: Whisker Assay, Derivative Assay, Generated

Journal: bioRxiv

Article Title: Data-driven fine-grained region discovery in the mouse brain with transformers

doi: 10.1101/2024.05.05.592608

Figure Lengend Snippet: Overall training and architectural scheme for CellTransformer. ( a. ) During training, a single cell is drawn (we denote this the reference cell, boxed in red). We extract the reference cell’s spatial neighbors and partition the group into a masked reference cell and its observed spatial neighbors. ( b. ) Initially, the model encoder receives information about each cell and projects those features to d- dimensional latent variable space. Features interact across cells (tokens) through the self-attention mechanism, which is repeated n times. These per-cell representations and an extra token are then aggregated into a single vector representation, which we refer to as the neighborhood representation. This representation is concatenated to a mask token which is cell type specific and chosen to represent the type of the reference cell. A shallow transformer decoder (dotted lines) further refines these representations and then a linear projection is used to output parameters of a negative binomial distribution modeling of the MERFISH probe counts for the reference cell.

Article Snippet: We downloaded the log-transformed MERFISH probe counts and metadata from the Allen Institute public beta ( https://alleninstitute.github.io/abc_atlas_access/intro.html ) access for ABC-MWB for both the Allen Institute for Brain Science and Zhuang lab mice.

Techniques: Plasmid Preparation

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , Left, schematic of experiment. Right, scRNA-seq identifies a diverse range of distinct cardiac cells that create the developing human heart as displayed by uniform manifold approximation and projection (UMAP) of ~143,000 cells. b , Schematic shows how 238 cardiac-cell-specific genes were spatially identified using MERFISH. Pseudo-coloured dots mark the location of individual molecules of ten specific RNA transcripts. c , Approximately 250,000 MERFISH-identified cardiac cells were clustered into specific cell populations as shown by UMAP and coloured accordingly in d . d , Identified MERFISH cells were spatially mapped across a frontal section of a 13 p.c.w. heart (left) and shown according to major cell classes (right). e , Joint embedding between MERFISH and age-matched scRNA-seq datasets enabled cell label transfer and MERFISH gene imputation. f , Co-occurrence heatmap shows the correspondence of cell annotations of MERFISH cells to those transferred from the 13 p.c.w. scRNA-seq dataset. g , Gene imputation performance was validated spatially by comparing normalized gene expression profiles of marker genes measured by MERFISH with the corresponding imputed gene expression profiles. Epi, epicardial; MV, mitral valve; P–RBC, platelet–red blood cell; TV, tricuspid valve. Scale bar, 250 µm ( g ). Illustration in a was created using BioRender ( https://www.biorender.com ).

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Gene Expression, Marker

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , Heatmap shows specific marker genes, as identified by NS-Forest2 classifier, for the 75 distinct cells across the developing heart. The distribution of these cells is shown according to age and region on bar graph. b , Cardiac single cells identified using the top 3,000 variable genes (left) and the 238 MERFISH genes (right) were visualized by UMAP which show that the transcriptional differences between the cell compartments (grey dashed lines) and classes (colored in a ) are preserved with a limited set of genes. aCM, atrial cardiomyocyte; BEC, blood endothelial cell; Epi, epicardial; Fibro, fibroblast; IVS, interventricular septum; LA, left atrium; LEC, lymphatic endothelial cell; LV, left ventricle; ncCM, non-chambered cardiomyocyte; p.c.w., post conception weeks; P-RBC, platelet-red blood cell; RA, right atrium; RV, right ventricle; SMC, smooth muscle cell; vCM, ventricular cardiomyocyte; WBC, white blood cell.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Marker

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , MERFISH cell boundaries were defined using CellPose with DAPI and polyA staining as input images. b , Pearson correlation of the counts of the 238 MERFISH target genes reveals strong correlation among the three replicate MERFISH experiments (Pearson correlation coefficient ( r ) > 0.95). c , Pearson correlation of the transcript counts of the 238 target genes shows that the 13 p.c.w. stage displays the highest average correlation (0.67 Pearson correlation) between the MERFISH and scRNA-seq datasets. d , MERFISH imaging was validated spatially by comparing normalized gene expression profiles of marker genes measured by single molecule FISH (smFISH) imaging with those detected by MERFISH imaging. e , Marker gene analysis identified each distinct MERFISH cell. f , Heatmap reveals that cell classes identified in the scRNA-seq dataset are detected in the MERFISH dataset, with the exception of P-RBCs. g , Table shows cellular composition similarities between the scRNA-seq and MERFISH datasets. aCM, atrial cardiomyocyte; aFibro, atrial fibroblast; adFibro, adventitial fibroblast; aEndocardial, atrial endocardial; AVC, atrioventricular canal; BEC, blood endothelial cell; CM, cardiomyocyte; EPDC, epicardial-derived cell; IFT, inflow tract; LA, left atrium; LEC, lymphatic endothelial cell; LV, left ventricle; ncCM, non-chambered cardiomyocyte; p.c.w., post conception weeks; P-RBC, platelet-red blood cell; RA, right atrium; RV, right ventricle; SMC, smooth muscle cell; vCM, ventricular cardiomyocyte; vCM-LV/RV-AV, muscular valve leaflet vCM; vEndocardial, ventricular endocardial; VEC, valve endocardial cell; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 50 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Staining, Imaging, Gene Expression, Marker, Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: The spatial mapping of each identified MERFISH cell is displayed accordingly: a , cardiomyocyte related cells, b , epicardial, EPDC, and vascular support related cells, c , endothelial related cells, d , neuronal cells, and e , blood related cells. aCM, atrial cardiomyocyte; aFibro, atrial fibroblast; adFibro, adventitial fibroblast; aEndocardial, atrial endocardial; AVC, atrioventricular canal; BEC, blood endothelial cell; EPDC, epicardial-derived cell; IFT, inflow tract; LA, left atrium; LEC, lymphatic endothelial cell; LV, left ventricle; ncCM, non-chambered cardiomyocyte; RA, right atrium; RV, right ventricle; vCM, ventricular cardiomyocyte; vCM-LV/RV-AV, muscular valve leaflet vCM; vEndocardial, ventricular endocardial; VEC, valve endocardial cell; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , Spatial mapping of identified MERFISH cells on two additional 13 p.c.w. frontal heart section replicates reveals the reproducibility of each distinct MERFISH cell and their spatial distributions. aCM, atrial cardiomyocyte; aFibro, atrial fibroblast; adFibro, adventitial fibroblast; aEndocardial, atrial endocardial; AVC, atrioventricular canal; BEC, blood endothelial cell; CM, cardiomyocyte; EPDC, epicardial-derived cell; IFT, inflow tract; LA, left atrium; LEC, lymphatic endothelial cell; LV, left ventricle; ncCM, non-chambered cardiomyocyte; p.c.w., post conception weeks; RA, right atrium; RV, right ventricle; vCM, ventricular cardiomyocyte; vCM-LV/RV-AV, muscular valve leaflet vCM; vEndocardial, ventricular endocardial; VEC, valve endocardial cell; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , Interrogation of the cell content around each individual cell identified cell zones or neighbourhoods, which formed defined CCs. b , Spatial mapping of CCs onto 13 p.c.w. hearts revealed their correspondence to distinct anatomical cardiac structures. c , The spatial location of each CC is displayed along with examples of their cellular composition and distribution (insets). d , Heatmap shows the composition of identified MERFISH cells within each defined CC. e , f , Analysis of the number of unique cell populations within each zone reveals the cellular complexity of each CC and cardiac region as displayed quantitatively ( e , violin plot) and spatially ( f , spatial complexity map). For e , the centre white dot represents the median, the bold black line represents the interquartile range, and the edges define minima and maxima of the distribution. Boxed areas in the spatial complexity map show regions of low (i) and high (ii) complexity. Insets (middle show the respective cellular composition, and magnified insets (right) show distinct identified cells). Mus. valve leaf., muscular valve leaflet. Scale bar, 250 µm ( b , c , f ).

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques:

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , MERFISH cells that constitute the ventricles (left, orange) were clustered as displayed using UMAP (right). b , Identified ventricular cells were spatially mapped in 13 p.c.w. ventricles. c , The spatial distributions of specific ventricular cells are shown for the left ventricular wall from the region outlined in the MERFISH spatial map in b . d , The ventricular wall depth distribution of ventricular cells is shown as a measured distance from the epicardial/outer surface of the ventricle for the imaged region in b . e , LV vCMs segregated into distinct vCM subpopulations. f , The molecular relationship of distinct vCMs is displayed in a connectivity map in which weighted edges between nodes represent their connectivity based on gene expression similarity. g , Heatmap shows the normalized expression of differentially expressed genes for vCMs as ordered by increasing ventricular wall depth. The coloured bar at the bottom indicates the specific vCMs as denoted in b . h , Scatter plot reveals the relationship between ventricular wall depth and pseudotime for individual vCMs in the left ventricle. i , MERFISH images of outlined regions in c ((i) and (ii)) show that specific combinations of gene markers, as shown in green and red, spatially identified specific vCMs. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Gene Expression, Expressing

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , Gene marker analysis defined MERFISH cells clustered from only the ventricles. b , MERFISH images reveal that spatial expression of genes related to specific vCMs correlate with ventricular wall depth. c , UMAP shows pseudotime of these vCMs within the left ventricular wall. d , Gene expression of specific markers for each distinct vCM is plotted along the pseudotime axis. Colored lines indicate each gene examined (see legend above plots). BEC, blood endothelial cell; EPDC, epicardial-derived cell; IVS, interventricular septum; LEC, lymphatic endothelial cell; LV, left ventricle; RV, right ventricle; vCM, ventricular cardiomyocyte; vCM-LV/RV-AV, muscular valve leaflet vCM; VEC, valve endocardial cell; vEndocardial, ventricular endocardial; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Marker, Expressing, Gene Expression, Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: The spatial mapping of each identified ventricular MERFISH cell is displayed accordingly: a , cardiomyocyte related cells, b , vascular support related cells, c , neuronal cells, d , epicardial, EPDC, and fibroblast-related cells, and e , WBC related cells. BEC, blood endothelial cell; EPDC, epicardial-derived cell; IVS, interventricular septum; LEC, lymphatic endothelial cell; LV, left ventricle; Prolif., proliferating; RV, right ventricle; vCM, ventricular cardiomyocyte; vCM-LV/RV-AV, muscular valve leaflet vCM; VEC, valve endocardial cell; vEndocardial, ventricular endocardial; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , MERFISH cells composing 15 post conception weeks (p.c.w.) ventricles were clustered as displayed on UMAP (left), and the identified ventricular cells were spatially mapped onto the ventricles and labeled in legend (right). b , Heatmap of transcriptional correlation between the MERFISH ventricular subpopulations shows that the 15 p.c.w. MERFISH dataset contained all cardiac cells previously identified by the 13 p.c.w. MERFISH dataset, except for the vCM-LV-Hybrid and vCM-RV-Hybrid cardiac cell subpopulations. c , The spatial distribution of specific ventricular cardiomyocytes is shown for the left ventricular wall from region outlined in MERFISH spatial map ( a ). d , Bar graph shows the relative cell composition of 13 p.c.w. and 15 p.c.w. ventricles. e , Bar graph of hybrid vCMs identified at specific scRNA-seq developmental stages reveals the proportion of hybrid vCMs to total vCMs in the LV from 9–15 p.c.w. f , Violin plots show the comparison of normalized ventricular wall depths of distinct ventricular cells within the apical/free walls at 13 p.c.w. and 15 p.c.w. The center dashed line represents the median, the other two dashed lines represent the interquartile range, and the edges define minima and maxima of the distribution. aFibro, atrial fibroblast; BEC, blood endothelial cell; EPDC, epicardial-derived cell; Fibro, fibroblast; IVS, interventricular septum; LEC, lymphatic endothelial cell; LV, left ventricle; Prolif., proliferating; RV, right ventricle; vCM, ventricular cardiomyocyte; vCM-AV, muscular valve leaflet vCM; vCM-LV/RV-AV, muscular valve leaflet vCM; VEC, valve endocardial cell; vEndocardial, ventricular endocardial; vFibro, ventricular fibroblast; VIC, valve interstitial cell; VSMC, vascular smooth muscle cell; WBC, white blood cell. Scale bar, 250 µm.

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques: Labeling, Comparison, Derivative Assay

Journal: Nature

Article Title: Spatially organized cellular communities form the developing human heart

doi: 10.1038/s41586-024-07171-z

Figure Lengend Snippet: a , MERFISH-identified ventricular cells assembled into nine more refined CCs within the ventricle. b , Heatmap shows the composition of distinct ventricular cells within each ventricle CC. c , MERFISH image of the outlined area in a reveals CC layers and their cell composition. Violin plot shows the ventricular wall depth distributions for distinct ventricular cells within these layers. The centre white dot represents the median, the bold black line represents the interquartile range, and the edges define minima and maxima of the distribution. Dashed lines indicate boundaries for CC layers. d , Chord diagrams reveal the strength of cell–cell signalling interactions received by specific vCMs in the inner-LV, intermediate-LV and outer-LV CCs. The size of the node represents the number of cells for a distinct ventricular cell, and the width of the edge represents the interaction strength between pairs of specific ventricular cells. e , The Venn diagram shows the number of specific and shared CCIs received by vCMs within the inner-LV, intermediate-LV and outer-LV communities. f , Dot plot shows specific signalling interactions between distinct ventricular cells within the intermediate-LV CC. g , Left, spatial map of cells participating in interactions between SEMA3C , SEMA3D , SEMA 6A or SEMA6B with PLXNA2 or PLXN4 for the intermediate-LV CC. Right, normalized ventricular wall depth distribution of these cells is shown in the histogram. h , High-resolution spatial cell map of the intermediate-LV CC shows how cells involved in interactions with SEMA3C , SEMA3D , SEMA 6A or SEMA6B with PLXNA2 or PLXN4 signalling may be spatially distributed to mediate attracting or repelling interactions. Arrows and arrowheads point to SEMA3C + SEMA3D + compact vFibro cells and SEMA6A + SEMA6B + BECs, respectively. Fibro/Epi, fibroblast and epicardial; His/mus. valve leaf., bundle of His and the muscular valve leaflet; Int., Intermediate; Out., Outer. Scale bars, 50 µm ( g , h ); 250 µm ( a ).

Article Snippet: The pipeline for processing the MERFISH dataset, including cell segmentation and assigning barcodes to cells, are available from GitHub ( https://github.com/epigen-UCSD/merfish_tools ).

Techniques:

Journal: bioRxiv

Article Title: Nicheformer: a foundation model for single-cell and spatial omics

doi: 10.1101/2024.04.15.589472

Figure Lengend Snippet: A) Nicheformer is pretrained on the SpatialCorpus-110M, a large data collection of over 110 million cells measured with dissociated and image-based spatial transcriptomics technologies. The SpatialCorpus-110M collection comprises single-cell data from Homo Sapiens and Mus Musculus across 17 distinct organs, 18 cell lines, and additional single-cell data from other anatomical systems and junctions. Shown is an exemplary UMAP visualization of a random 1% subset of the entire pretraining dataset (n=1,108,759 cells) of the non-integrated log1p-transformed normalized SpatialCorpus-110M colored by modality. B) Nicheformer includes a novel set of downstream tasks, ranging from spatial cell type, niche and region label prediction to neighborhood cell density and neighborhood composition prediction. We test our approach on large-scale, high-quality spatial transcriptomics data from the brain (mouse - MERFISH), liver (CosMx - human), lung (CosMx - human, Xenium - human), and colon (Xenium - human). Visualized are example slices of the respective datasets colored by niche labels (brain, liver, and lung) and cell density (lung and colon). C) The SpatialCorpus-110M is harmonized and mapped to orthologous gene names, as well as human and mouse-specific genes, to create the input for Nicheformer pretraining. We harmonized metadata information across all datasets, capturing species, modality, and assay. D) Each cell’s gene expression profile and metadata are fed into a gene rank tokenizer to obtain a tokenized representation for each cell. The tokenized cells serve as input for the Nicheformer transformer block to predict masked tokens. Finally, the Nicheformer embedding is generated by aggregating the gene tokens (Methods). E) The pretrained Nicheformer embedding is visualized as UMAP colored by modality. The UMAP shows a random 5% subsample of the entire Nicheformer embedding (n=4,903,086).

Article Snippet: For spatial transcriptomics, we curated image-based spatial datasets, specifically MERFISH (Vizgen MERSCOPE), 10x Genomics Xenium, Nanostring CosMx , and In Situ Sequencing (ISS) data ( , Supp.

Techniques: Transformation Assay, Gene Expression, Blocking Assay, Generated

Journal: bioRxiv

Article Title: Nicheformer: a foundation model for single-cell and spatial omics

doi: 10.1101/2024.04.15.589472

Figure Lengend Snippet: A) Single-cells resolved in space on an example slice (n=114,396 cells) of the MERFISH mouse brain dataset with niche label superimposed. B) Test-set F1-macro of niche and brain region label prediction of the fine-tuned Nicheformer model, the linear probing model, and a linear probing baseline computed based on embeddings generated with scVI and PCA, respectively. C) UMAP of dissociated scRNAseq dataset with original author cell type label superimposed. D) Nicheformer can transfer spatial niche and region labels onto dissociated single-cell data. E) Nicheformer accurately classifies cells from the dissociated motor cortex to relevant cell types (n=9 out of 33 distinct ones in the classifier) trained on the whole mouse brain MERFISH dataset. F-G) Nicheformer correctly projects dissociated single-cells to niche (F) and region (G) labels to provide spatially dependent labels. F) Nicheformer misclassified parts of L2/3 IT neurons as residing in the subpallium GABAergic niche (highlighted in the red box). Additionally, the deep cortical excitatory neurons L6b, L6 CT, L6 IT, and L6 IT Car3 (highlighted in the red box) should be classified as pallium glutamatergic niche instead of subpallium GABAergic by Nicheformer. G) Most of the non-neuronal cells (84.7 % of all non-neuronal cells, n=133) were misclassified as not belonging to the isocortex or the adjacent brain regions (highlighted in the red box). H) Cell type abundances in the scRNA-seq dataset measuring the primary motor cortex in the mouse. I-K) Classification uncertainty of label transfer of the dissociated scRNA-seq dataset to the MERFISH mouse brain data for cell type label (I), niche label (J), and region label (K) with a value of 0 representing a high uncertainty and 1 being a lower uncertainty, i.e., high certainty. K) Observed high uncertainty for parts of the Glut and GABA neurons for the region prediction of the isocortex, CTXsp and OLF, which are neighboring brain regions.

Article Snippet: For spatial transcriptomics, we curated image-based spatial datasets, specifically MERFISH (Vizgen MERSCOPE), 10x Genomics Xenium, Nanostring CosMx , and In Situ Sequencing (ISS) data ( , Supp.

Techniques: Generated

Journal: bioRxiv

Article Title: Nicheformer: a foundation model for single-cell and spatial omics

doi: 10.1101/2024.04.15.589472

Figure Lengend Snippet: A) We define the neighborhood of a cell as its local neighborhood given a radius and an index cell. The neighborhood cell density is then defined by the number of cells in the neighborhood, and the neighborhood compositions are the proportions of neighboring cell types. B) Neighborhoods are computed at multiple resolutions resulting in different neighborhood size distributions. Each barplot shows the distribution of the number of neighbors across the brain, liver, and lung datasets. We extract neighborhoods with the mean number of neighbors 10, 20, 50, and 100 for each dataset. C) The fine-tuned and linear probing Nicheformer models outperform zero-shot models trained on scVI and PCA embeddings in terms of mean absolute error across all neighborhood sizes and all three organs, the brain, liver, and lung. D) Left: Fine-tuned Nicheformer performance on the MERFISH mouse brain data grouped by index cell type. Shown are the absolute error values between predicted and observed neighborhood composition vectors for held-out test cells. For each box in (D), the centerline defines the median, the height of the box is given by the interquartile range (IQR), the whiskers are given by 1.5 × IQR, and outliers are given as points beyond the minimum or maximum whisker. Center: Index cell type abundances in the entire MERFISH mouse brain dataset. Right: UMAPs of MERFISH mouse brain Nicheformer embedding with the selected index cell type as color superimposed. E) UMAP of the Nicheformer embedding of all immune cells in the MERFISH mouse brain dataset with region label as color superimposed.

Article Snippet: For spatial transcriptomics, we curated image-based spatial datasets, specifically MERFISH (Vizgen MERSCOPE), 10x Genomics Xenium, Nanostring CosMx , and In Situ Sequencing (ISS) data ( , Supp.

Techniques: Whisker Assay

Journal: Nature Communications

Article Title: Aging in mice alters regionally enriched striatal astrocytes

doi: 10.1038/s41467-025-63429-8

Figure Lengend Snippet: a XY directional map of all cell types identified within the striatal section is shown. Medium spiny neurons, cortical neurons, astrocytes, oligodendrocytes, endothelial cells, and microglia are displayed in different colors - each dot represents a cell segmented by MERLIN. All cells within the XY coordinates of the striatum were subset, and here we display the UMAP of these subset striatal cell clusters split by major cell types and UMAP split by age (young = blue; aged = red). b Striatal astrocytes were subset from all striatal cells and clustered separately using the Louvian algorithm. The first UMAP shows striatal astrocyte subtypes from single cells, the second UMAP shows striatal astrocyte subtypes from MERFISH, and the third UMAP shows striatal astrocyte subtypes from integrated single-cell and MERFISH datasets. c – f Top 4 astrocyte subtypes (by abundance) are shown. The astrocyte subtype expression probability is quantified along the dorsal-ventral axis in 500 μm segments in young and aged mice. We divided the striatum into five 500 μm sections, starting at the base of the corpus callosum and moving ventrally. We quantified the density of each astrocyte subtype within each subregion. This astrocyte expression probability quantification was calculated by the number of astrocytes within a subcluster (A1–7), within each 500 μm subregion (0–5) ( X A1…A7 within Y 0…5 ) divided by the total number of astrocytes within that 500 μm section (Σ total ) normalized to the total number of astrocytes within each respective subcluster ( σ A1…A7 ) ([( X A1,..A7 within Y 0…5 /Σ total )/ σ A1…A7 ]). The regional change is quantified by subtracting the young astrocyte expression probability from the aged expression probability. These quantifications were statistically analyzed using a two-way repeated measures (for subregion) ANOVA. Asterisks (*) indicate significant differences ( p value < 0.05) across sub-regions, and hashtags (#) indicate significant differences across ages. Individual representative astrocyte maps for young and aged striatal sections are displayed to the right of the astrocyte density quantification, with astrocyte subtypes demarcated in their respective colors. In the graphs shown in ( c – f ), the corpus callosum is abbreviated as CC on the y -axis.

Article Snippet: MERFISH imaging was performed on an automated Vizgen Alpha Instrument using imaging buffers, hybridization buffers, and parameter files provided by Vizgen.

Techniques: Expressing

Journal: Nature Communications

Article Title: Aging in mice alters regionally enriched striatal astrocytes

doi: 10.1038/s41467-025-63429-8

Figure Lengend Snippet: a Top aging astrocyte markers were assessed using MAST differential expression, and the top 25 up and 25 downregulated transcripts are shown in the heatmap (see also Supplementary Data ). Dorsal enriched genes are marked with an asterisk (*) and ventral enriched genes are marked with a hashtag (#). b Representative image of MERFISH RNA counts shows age increases in Gfap in the dorsal striatum. c Representative images of the dorsal striatum are shown for young and aged mice. GFAP coverage per 500 μm 2 in the dorsal striatum; across young and aged mice ( n = 4–5 mice; 2-way RM ANOVA with Bonferroni post hoc (* p = 0.00032 (aged dorsal compared ventral) and * p = 0.00029 (aged dorsal compared to young dorsal), data are presented as mean values ± SEM). d Representative image of MERFISH RNA counts shows S100b expression in the dorsal striatum. e Quantification of S100β + cells per 500 μm 2 , in the dorsal, medial, and ventral striatum; across young and aged mice ( n = 4–5 mice; 2-way RM ANOVA with Bonferroni post hoc * p = 0.036, data are presented as mean values ± SEM). f IPA analysis was performed on age-induced DEGs in striatal astrocytes, and each black bar indicates the number of genes per pathway, circle size indicates the −log( p value) (right-tailed Fisher’s Exact Test), and circle color indicates the activation score number. g Aging gene score is calculated by the average expression levels of the top 20 differentially expressed genes in age on a single-cell level, subtracted by the aggregated expression of control feature sets. Each dot represents an astrocyte, and the color of the dot is the relative change in the aging score. h Heatmaps of the top 10 shared genes between our aging and A1 astrocyte subtype genes (Log2FC > 0.01) with mouse aging (Log2FC > 0.01), human striatal astrocyte aging, human striatal astrocyte Huntington’s Disease, and human Parkinson’s Disease using the DEG (see data availability and Supplementary Data ). Venn diagrams visualize the total overlap across each gene set. i UpSet plot of the overlap of murine aging striatal data and the past studies.

Article Snippet: MERFISH imaging was performed on an automated Vizgen Alpha Instrument using imaging buffers, hybridization buffers, and parameter files provided by Vizgen.

Techniques: Quantitative Proteomics, Expressing, Activation Assay, Control