Cell type-specific inference of differential expression in spatial transcriptomics

STANCE: a unified statistical model to detect cell-type-specific spatially variable genes in spatial transcriptomics

One of the major challenges in spatial transcriptomics is to detect spatially variable genes (SVGs), whose expression patterns are non-random across tissue locations. Many SVGs correlate with cell type compositions, introducing the concept of cell type-specific SVGs (ctSVGs). Existing ctSVG detection methods treat cell type-specific spatial effects as fixed effects, leading to tissue spatial rotation-dependent results. Moreover, SVGs may exhibit random spatial patterns within cell types, meaning an SVG is not always a ctSVG, and vice versa, further complicating detection. We propose STANCE, a unified statistical model for both SVGs and ctSVGs detection under a linear mixed-effect model framework that integrates gene expression, spatial location, and cell type composition information. STANCE ensures tissue rotation-invariant results, with a two-stage approach: initial SVG/ctSVG detection followed by ctSVG-specific testing. We demonstrate its performance through extensive simulations and analyses of public datasets. Downstream analyses reveal STANCE's potential in spatial transcriptomics analysis.

STMiner: Gene-centric spatial transcriptomics for deciphering tumor tissues

Analyzing spatial transcriptomics data from tumor tissues poses several challenges beyond those of healthy samples, including unclear boundaries between different regions, uneven cell densities, and relatively higher cellular heterogeneity. Collectively, these bias the background against which spatially variable genes are identified, which can result in misidentification of spatial structures and hinder potential insight into complex pathologies. To overcome this problem, STMiner leverages 2D Gaussian mixture models and optimal transport theory to directly characterize the spatial distribution of genes rather than the capture locations of the cells expressing them (spots). By effectively mitigating the impacts of both background bias and data sparsity, STMiner reveals key gene sets and spatial structures overlooked by spot-based analytic tools, facilitating novel biological discoveries. The core concept of directly analyzing overall gene expression patterns also allows for a broader application beyond spatial transcriptomics, positioning STMiner for continuous expansion as spatial omics technologies evolve.

Precise gene expression deconvolution in spatial transcriptomics with STged

Spatially resolved transcriptomics (SRT) has transformed tissue biology by linking gene expression profiles with spatial information. However, sequencing-based SRT methods aggregate signals from multiple cell types within capture locations ("spots"), masking cell-type-specific gene expression patterns. Traditional cell-type deconvolution methods estimate cell compositions within spots but fail to resolve cell-type-specific gene expression, limiting their ability to uncover critical biological processes such as cellular interactions and microenvironmental dynamics. Here, we present STged (spatial transcriptomic gene expression deconvolution), a novel computational framework that goes beyond traditional deconvolution by reconstructing cell-type-specific gene expression profiles from mixed spots. STged integrates graph-based spatial correlations and reference-derived gene signatures using a non-negative least-squares regression framework, achieving precise and biologically meaningful deconvolution. Comprehensive simulations show that STged consistently outperforms existing methods in accuracy and robustness. Applications to human pancreatic ductal adenocarcinoma and human squamous cell carcinoma datasets reveal its capacity to identify microenvironment-specific highly variable genes, reconstruct spatial cell-cell communication networks, and resolve tissue architecture at near-single-cell resolution. In mouse kidney tissues, STged uncovers dynamic spatial gene expression patterns and distinct gene programs, advancing our understanding of tissue heterogeneity and cellular dynamics.

Mapping the topography of spatial gene expression with interpretable deep learning

Spatially resolved transcriptomics technologies provide high-throughput measurements of gene expression in a tissue slice, but the sparsity of these data complicates analysis of spatial gene expression patterns. We address this issue by deriving a topographic map of a tissue slice-analogous to a map of elevation in a landscape-using a quantity called the isodepth. Contours of constant isodepths enclose domains with distinct cell type composition, while gradients of the isodepth indicate spatial directions of maximum change in expression. We develop GASTON (gradient analysis of spatial transcriptomics organization with neural networks), an unsupervised and interpretable deep learning algorithm that simultaneously learns the isodepth, spatial gradients and piecewise linear expression functions that model both continuous gradients and discontinuous variation in gene expression. We show that GASTON accurately identifies spatial domains and marker genes across several tissues, gradients of neuronal differentiation and firing in the brain, and gradients of metabolism and immune activity in the tumor microenvironment.

Single-cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation

The R47H missense mutation of the TREM2 gene is a known risk factor for development of Alzheimer's Disease. In this study, we analyze the impact of the Trem2R47H mutation on specific cell types in multiple cortical and subcortical brain regions in the context of wild-type and 5xFAD mouse background. We profile 19 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics, a technique that enables subcellular profiling of spatial gene expression. Spatial transcriptomics and neuropathology data are analyzed using our custom pipeline to identify plaque and Trem2R47H-induced transcriptomic dysregulation. We initially analyze cell type-specific transcriptomic alterations induced by plaque proximity. Next, we analyze spatial distributions of disease associated microglia and astrocytes, and how they vary between 5xFAD and Trem2R47H; 5xFAD mouse models. Finally, we analyze the impact of the Trem2R47H mutation on neuronal transcriptomes. The Trem2R47H mutation induces consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice. Overall, our MERFISH spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

A rare PRIMER cell state in plant immunity

Plants lack specialized and mobile immune cells. Consequently, any cell type that encounters pathogens must mount immune responses and communicate with surrounding cells for successful defence. However, the diversity, spatial organization and function of cellular immune states in pathogen-infected plants are poorly understood1. Here we infect Arabidopsis thaliana leaves with bacterial pathogens that trigger or supress immune responses and integrate time-resolved single-cell transcriptomic, epigenomic and spatial transcriptomic data to identify cell states. We describe cell-state-specific gene-regulatory logic that involves transcription factors, putative cis-regulatory elements and target genes associated with disease and immunity. We show that a rare cell population emerges at the nexus of immune-active hotspots, which we designate as primary immune responder (PRIMER) cells. PRIMER cells have non-canonical immune signatures, exemplified by the expression and genome accessibility of a previously uncharacterized transcription factor, GT-3A, which contributes to plant immunity against bacterial pathogens. PRIMER cells are surrounded by another cell state (bystander) that activates genes for long-distance cell-to-cell immune signalling. Together, our findings suggest that interactions between these cell states propagate immune responses across the leaf. Our molecularly defined single-cell spatiotemporal atlas provides functional and regulatory insights into immune cell states in plants.

Categorization of 34 computational methods to detect spatially variable genes from spatially resolved transcriptomics data

In the analysis of spatially resolved transcriptomics data, detecting spatially variable genes (SVGs) is crucial. Numerous computational methods exist, but varying SVG definitions and methodologies lead to incomparable results. We review 34 state-of-the-art methods, classifying SVGs into three categories: overall, cell-type-specific, and spatial-domain-marker SVGs. Our review explains the intuitions underlying these methods, summarizes their applications, and categorizes the hypothesis tests they use in the trade-off between generality and specificity for SVG detection. We discuss challenges in SVG detection and propose future directions for improvement. Our review offers insights for method developers and users, advocating for category-specific benchmarking.

Statistical identification of cell type-specific spatially variable genes in spatial transcriptomics

An essential task in spatial transcriptomics is identifying spatially variable genes (SVGs). Here, we present Celina, a statistical method for systematically detecting cell type-specific SVGs (ct-SVGs)-a subset of SVGs exhibiting distinct spatial expression patterns within specific cell types. Celina utilizes a spatially varying coefficient model to accurately capture each gene's spatial expression pattern in relation to the distribution of cell types across tissue locations, ensuring effective type I error control and high power. Celina proves powerful compared to existing methods in single-cell resolution spatial transcriptomics and stands as the only effective solution for spot-resolution spatial transcriptomics. Applied to five real datasets, Celina uncovers ct-SVGs associated with tumor progression and patient survival in lung cancer, identifies metagenes with unique spatial patterns linked to cell proliferation and immune response in kidney cancer, and detects genes preferentially expressed near amyloid-β plaques in an Alzheimer's model.

Elucidating cardiomyocyte heterogeneity and maturation dynamics through integrated single-cell and spatial transcriptomics

The intricate development and functionality of the mammalian heart are influenced by the heterogeneous nature of cardiomyocytes (CMs). In this study, single-cell and spatial transcriptomics were utilized to analyze cells from neonatal mouse hearts, resulting in a comprehensive atlas delineating the transcriptional profiles of distinct CM subsets. A continuum of maturation states was elucidated, emphasizing a progressive developmental trajectory rather than discrete stages. This approach enabled the mapping of these states across various cardiac regions, illuminating the spatial organization of CM development and the influence of the cellular microenvironment. Notably, a subset of transitional CMs was identified, characterized by a transcriptional signature marking a pivotal maturation phase, presenting a promising target for therapeutic strategies aimed at enhancing cardiac regeneration. This atlas not only elucidates fundamental aspects of cardiac development but also serves as a valuable resource for advancing research into cardiac physiology and pathology, with significant implications for regenerative medicine.

SPathDB: a comprehensive database of spatial pathway activity atlas

Spatial transcriptomics sequencing technology deepens our understanding of the diversity of cell behaviors, fates and states within complex tissue, which is often determined by the fine-tuning of regulatory network functional activities. Therefore, characterizing the functional activity within tissue space is helpful for revealing the functional features that drive spatial heterogeneity, and understanding complex biological processes. Here, we describe a database, SPathDB (http://bio-bigdata.hrbmu.edu.cn/SPathDB/), which aims to dissect the pathway-mediated multidimensional spatial heterogeneity in the context of functional activity. We manually curated spatial transcriptomics datasets and biological pathways from public data resources. SPathDB consists of 1689 868 spatial spots of 695 slices from 84 spatial transcriptome datasets of human and mouse, which involves 36 tissues, and also diseases such as cancer, and provides interactive analysis and visualization of the functional activities of 114 998 pathways across these spatial spots. SPathDB provides five flexible interfaces to retrieve and analyze pathways with highly variable functional activity across spatial spots, the distribution of pathway functional activities along pseudo-space axis, pathway-mediated spatial intercellular communications and the associations between spatial pathway functional activity and the occurrence of cell types. SPathDB will serve as a foundational resource for identifying functional features and elucidating underlying mechanisms of spatial heterogeneity.

Spatiotemporal modeling of molecular holograms

Quantifying spatiotemporal dynamics during embryogenesis is crucial for understanding congenital diseases. We developed Spateo (https://github.com/aristoteleo/spateo-release), a 3D spatiotemporal modeling framework, and applied it to a 3D mouse embryogenesis atlas at E9.5 and E11.5, capturing eight million cells. Spateo enables scalable, partial, non-rigid alignment, multi-slice refinement, and mesh correction to create molecular holograms of whole embryos. It introduces digitization methods to uncover multi-level biology from subcellular to whole organ, identifying expression gradients along orthogonal axes of emergent 3D structures, e.g., secondary organizers such as midbrain-hindbrain boundary (MHB). Spateo further jointly models intercellular and intracellular interaction to dissect signaling landscapes in 3D structures, including the zona limitans intrathalamica (ZLI). Lastly, Spateo introduces "morphometric vector fields" of cell migration and integrates spatial differential geometry to unveil molecular programs underlying asymmetrical murine heart organogenesis and others, bridging macroscopic changes with molecular dynamics. Thus, Spateo enables the study of organ ecology at a molecular level in 3D space over time.

Identifying cell-type-specific spatially variable genes with ctSVG

Spatially variable genes (SVGs) reveal the molecular and functional heterogeneity of cells across different spatial regions of a tissue. We found that sample-wide SVGs, identified by previous methods across the whole sample, largely overlap with cell-type marker genes derived from single-cell gene expression, leaving the spatial location information largely underutilized. We developed ctSVG, a computational method specifically tailored for Visium HD spatial transcriptomics at single-cell resolution. ctSVG accurately assigns Visium squares to cells and identifies cell-type-specific SVGs. We show that cell-type-specific SVGs identified by ctSVG include many new genes that do not overlap with sample-wide SVGs or cell-type marker genes, and that these genes reveal important biological functions in real spatial datasets.

NiCo identifies extrinsic drivers of cell state modulation by niche covariation analysis

Cell states are modulated by intrinsic driving forces such as gene expression noise and extrinsic signals from the tissue microenvironment. The distinction between intrinsic and extrinsic cell state determinants is essential for understanding the regulation of cell fate in tissues during development, homeostasis and disease. The rapidly growing availability of single-cell resolution spatial transcriptomics makes it possible to meet this challenge. However, available computational methods to infer topological tissue domains, spatially variable genes, or ligand-receptor interactions are limited in their capacity to capture cell state changes driven by crosstalk between individual cell types within the same niche. We present NiCo, a computational framework for integrating single-cell resolution spatial transcriptomics with matched single-cell RNA-sequencing reference data to infer the influence of the spatial niche on the cell state. By applying NiCo to mouse embryogenesis, adult small intestine and liver data, we demonstrate the ability to predict novel niche interactions that govern cell state variation underlying tissue development and homeostasis. In particular, NiCo predicts a feedback mechanism between Kupffer cells and neighboring stellate cells dampening stellate cell activation in the normal liver. NiCo provides a powerful tool to elucidate tissue architecture and to identify drivers of cellular states in local niches.

Computational modeling for deciphering tissue microenvironment heterogeneity from spatially resolved transcriptomics

Spatial transcriptomics techniques, while measuring gene expression, retain spatial location information, aiding in situ studies of organismal tissue architecture and the progression of pathological processes. These techniques generate vast amounts of omics data, necessitating the development of computational methods to reveal the underlying tissue microenvironment heterogeneity. The main directions in spatial transcriptomics data analysis are spatial domain detection and spatial deconvolution, which can identify spatial functional regions and parse the distribution of cell types in spatial transcriptomics data by integrating single-cell transcriptomics data. In these two research directions, many computational methods have been successively proposed. This article will categorize them into three types: machine learning-based methods, probabilistic models-based methods, and deep learning-based methods. It will list and discuss the representative algorithms of each type along with their advantages and disadvantages and describe the datasets and evaluation metrics used to assess these computational methods, facilitating researchers in selecting suitable computational methods according to their research needs. Finally, combining the latest technological developments and the advantages and disadvantages of current algorithms, this article will look forward to the future directions of computational method development.

STGAT: Graph attention networks for deconvolving spatial transcriptomics data

BACKGROUND AND OBJECTIVE

Spatially resolved gene expression profiles are crucial for understanding tissue structure and function. However, the lack of single-cell resolution in these profiles demands their integration with single-cell RNA sequencing data for accurate dataset deconvolution. We propose STGAT, an innovative deconvolution method that leverages graph attention networks to enhance spatial transcriptomic (ST) data analysis.

METHODS

STGAT generates pseudo-ST data that more comprehensively represents the cell-type composition within real-ST data by using three different sampling probabilities. A comprehensive combined graph is then constructed to capture the complex relationships both across pseudo- and real-ST data and within each dataset. Moreover, integrating a graph attention network further enables STGAT to dynamically assign the weights to the connections between spots, significantly enhancing the accuracy of cell-type composition predictions.

RESULTS

Extensive comparative experiments on simulated and real-world datasets, demonstrate the superior performance of STGAT for cell-type deconvolution. The method outperforms six established methods and is robust across various biological contexts.

CONCLUSION

STGAT exhibits more precise results in cell-type composition inference that are more consistent with the known knowledge, suggesting its potential utility in improving the resolution and accuracy of spatial transcriptomics data analysis.

Sensing and guiding cell-state transitions by using genetically encoded endoribonuclease-mediated microRNA sensors

Precisely sensing and guiding cell-state transitions via the conditional genetic activation of appropriate differentiation factors is challenging. Here we show that desired cell-state transitions can be guided via genetically encoded sensors, whereby endogenous cell-state-specific miRNAs regulate the translation of a constitutively transcribed endoribonuclease, which, in turn, controls the translation of a gene of interest. We used this approach to monitor several cell-state transitions, to enrich specific cell types and to automatically guide the multistep differentiation of human induced pluripotent stem cells towards a haematopoietic lineage via endothelial cells as an intermediate state. Such conditional activation of gene expression is durable and resistant to epigenetic silencing and could facilitate the monitoring of cell-state transitions in physiological and pathological conditions and eventually the 'rewiring' of cell-state transitions for applications in organoid-based disease modelling, cellular therapies and regenerative medicine.

Cell-specific priors rescue differential gene expression in spatial spot-based technologies

Spatial transcriptomics (ST), a breakthrough technology, captures the complex structure and state of tissues through the spatial profiling of gene expression. A variety of ST technologies have now emerged, most prominently spot-based platforms such as Visium. Despite the widespread use of ST and its distinct data characteristics, the vast majority of studies continue to analyze ST data using algorithms originally designed for older technologies such as single-cell (SC) and bulk RNA-seq-particularly when identifying differentially expressed genes (DEGs). However, it remains unclear whether these algorithms are still valid or appropriate for ST data. Therefore, here, we sought to characterize the performance of these methods by constructing an in silico simulator of ST data with a controllable and known DEG ground truth. Surprisingly, our findings reveal little variation in the performance of classic DEG algorithms-all of which fail to accurately recapture known DEGs to significant levels. We further demonstrate that cellular heterogeneity within spots is a primary cause of this poor performance and propose a simple gene-selection scheme, based on prior knowledge of cell-type specificity, to overcome this. Notably, our approach outperforms existing data-driven methods designed specifically for ST data and offers improved DEG recovery and reliability rates. In summary, our work details a conceptual framework that can be used upstream, agnostically, of any DEG algorithm to improve the accuracy of ST analysis and any downstream findings.

Identifying spatially variable genes by projecting to morphologically relevant curves

Spatial transcriptomics enables high-resolution gene expression measurements while preserving the two-dimensional spatial organization of the sample. A common objective in spatial transcriptomics data analysis is to identify spatially variable genes within predefined cell types or regions within the tissue. However, these regions are often implicitly one-dimensional, making standard two-dimensional coordinate-based methods less effective as they overlook the underlying tissue organization. Here we introduce a methodology grounded in spectral graph theory to elucidate a one-dimensional curve that effectively approximates the spatial coordinates of the examined sample. This curve is then used to establish a new coordinate system that better reflects tissue morphology. We then develop a generalized additive model (GAM) to detect genes with variable expression in the new morphologically relevant coordinate system. Our approach directly models gene counts, thereby eliminating the need for normalization or transformations to satisfy normality assumptions. We demonstrate improved performance relative to current methods based on hypothesis testing, while also accurately estimating gene expression patterns and precisely identifying spatial loci where deviations from constant expression are observed. We validate our approach through extensive simulations and by analyzing experimental data from multiple platforms such as Slide-seq and MERFISH.

A spatiotemporal molecular atlas of mouse spinal cord injury identifies a distinct astrocyte subpopulation and therapeutic potential of IGFBP2

Spinal cord injury (SCI) triggers a cascade of intricate molecular and cellular changes that determine the outcome. In this study, we resolve the spatiotemporal organization of the injured mouse spinal cord and quantitatively assess in situ cell-cell communication following SCI. By analyzing existing single-cell RNA sequencing datasets alongside our spatial data, we delineate a subpopulation of Igfbp2-expressing astrocytes that migrate from the white matter (WM) to gray matter (GM) and become reactive upon SCI, termed Astro-GMii. Further, Igfbp2 upregulation promotes astrocyte migration, proliferation, and reactivity, and the secreted IGFBP2 protein fosters neurite outgrowth. Finally, we show that IGFBP2 significantly reduces neuronal loss and remarkably improves the functional recovery in a mouse model of SCI in vivo. Together, this study not only provides a comprehensive molecular atlas of SCI but also exemplifies how this rich resource can be applied to endow cells and genes with functional insight and therapeutic potential.

Recovering single-cell expression profiles from spatial transcriptomics with scResolve

Many popular spatial transcriptomics techniques lack single-cell resolution. Instead, these methods measure the collective gene expression for each location from a mixture of cells, potentially containing multiple cell types. Here, we developed scResolve, a method for recovering single-cell expression profiles from spatial transcriptomics measurements at multi-cellular resolution. scResolve accurately restores expression profiles of individual cells at their locations, which is unattainable with cell type deconvolution. Applications of scResolve on human breast cancer data and human lung disease data demonstrate that scResolve enables cell-type-specific differential gene expression analysis between different tissue contexts and accurate identification of rare cell populations. The spatially resolved cellular-level expression profiles obtained through scResolve facilitate more flexible and precise spatial analysis that complements raw multi-cellular level analysis.

Categorization of 33 computational methods to detect spatially variable genes from spatially resolved transcriptomics data

In the analysis of spatially resolved transcriptomics data, detecting spatially variable genes (SVGs) is crucial. Numerous computational methods exist, but varying SVG definitions and methodologies lead to incomparable results. We review 33 state-of-the-art methods, categorizing SVGs into three types: overall, cell-type-specific, and spatial-domain-marker SVGs. Our review explains the intuitions underlying these methods, summarizes their applications, and categorizes the hypothesis tests they use in the trade-off between generality and specificity for SVG detection. We discuss challenges in SVG detection and propose future directions for improvement. Our review offers insights for method developers and users, advocating for category-specific benchmarking.

Modulating endothelial cell dynamics in fat grafting: the impact of DLL4 siRNA via adipose stem cell extracellular vesicles

In the context of improving the efficacy of autologous fat grafts (AFGs) in reconstructive surgery, this study delineates the novel use of adipose-derived mesenchymal stem cells (ADSCs) and their extracellular vesicles (EVs) as vehicles for delivering delta-like ligand 4 (DLL4) siRNA. The aim was to inhibit DLL4, a gene identified through transcriptome analysis as a critical player in the vascular endothelial cells of AFG tissues, thereby negatively affecting endothelial cell functions and graft survival through the Notch signaling pathway. By engineering ADSC EVs to carry DLL4 siRNA (ADSC EVs-siDLL4), the research demonstrated a marked improvement in endothelial cell proliferation, migration, and lumen formation, and enhanced angiogenesis in vivo, leading to a significant increase in the survival rate of AFGs. This approach presents a significant advancement in the field of tissue engineering and regenerative medicine, offering a potential method to overcome the limitations of current fat grafting techniques.NEW & NOTEWORTHY This study introduces a groundbreaking method for enhancing autologous fat graft survival using adipose-derived stem cell extracellular vesicles (ADSC EVs) to deliver DLL4 siRNA. By targeting the delta-like ligand 4 (DLL4) gene, crucial in endothelial cell dynamics, this innovative approach significantly improves endothelial cell functions and angiogenesis, marking a substantial advancement in tissue engineering and regenerative medicine.

Unveiling contact-mediated cellular crosstalk

Cell-cell interactions orchestrate complex functions in multicellular organisms, forming a regulatory network for diverse biological processes. Their disruption leads to disease states. Recent advancements - including single-cell sequencing and spatial transcriptomics, coupled with powerful bioengineering and molecular tools - have revolutionized our understanding of how cells respond to each other. Notably, spatial transcriptomics allows us to analyze gene expression changes based on cell proximity, offering a unique window into the impact of cell-cell contact. Additionally, computational approaches are being developed to decipher how cell contact governs the symphony of cellular responses. This review explores these cutting-edge approaches, providing valuable insights into deciphering the intricate cellular changes influenced by cell-cell communication.

Craters on the melanoma surface facilitate tumor-immune interactions and demonstrate pathologic response to checkpoint blockade in humans

Immunotherapy leads to cancer eradication despite the tumor's immunosuppressive environment. Here, we used extended long-term in-vivo imaging and high-resolution spatial transcriptomics of endogenous melanoma in zebrafish, and multiplex imaging of human melanoma, to identify domains that facilitate immune response during immunotherapy. We identified crater-shaped pockets at the margins of zebrafish and human melanoma, rich with beta-2 microglobulin (B2M) and antigen recognition molecules. The craters harbor the highest density of CD8+ T cells in the tumor. In zebrafish, CD8+ T cells formed prolonged interactions with melanoma cells within craters, characteristic of antigen recognition. Following immunostimulatory treatment, the craters enlarged and became the major site of activated CD8+ T cell accumulation and tumor killing that was B2M dependent. In humans, craters predicted immune response to ICB therapy, showing response better than high T cell infiltration. This marks craters as potential new diagnostic tool for immunotherapy success and targets to enhance ICB response.

Scalable imaging-free spatial genomics through computational reconstruction

Tissue organization arises from the coordinated molecular programs of cells. Spatial genomics maps cells and their molecular programs within the spatial context of tissues. However, current methods measure spatial information through imaging or direct registration, which often require specialized equipment and are limited in scale. Here, we developed an imaging-free spatial transcriptomics method that uses molecular diffusion patterns to computationally reconstruct spatial data. To do so, we utilize a simple experimental protocol on two dimensional barcode arrays to establish an interaction network between barcodes via molecular diffusion. Sequencing these interactions generates a high dimensional matrix of interactions between different spatial barcodes. Then, we perform dimensionality reduction to regenerate a two-dimensional manifold, which represents the spatial locations of the barcode arrays. Surprisingly, we found that the UMAP algorithm, with minimal modifications can faithfully successfully reconstruct the arrays. We demonstrated that this method is compatible with capture array based spatial transcriptomics/genomics methods, Slide-seq and Slide-tags, with high fidelity. We systematically explore the fidelity of the reconstruction through comparisons with experimentally derived ground truth data, and demonstrate that reconstruction generates high quality spatial genomics data. We also scaled this technique to reconstruct high-resolution spatial information over areas up to 1.2 centimeters. This computational reconstruction method effectively converts spatial genomics measurements to molecular biology, enabling spatial transcriptomics with high accessibility, and scalability.

Answering open questions in biology using spatial genomics and structured methods

Genomics methods have uncovered patterns in a range of biological systems, but obscure important aspects of cell behavior: the shapes, relative locations, movement, and interactions of cells in space. Spatial technologies that collect genomic or epigenomic data while preserving spatial information have begun to overcome these limitations. These new data promise a deeper understanding of the factors that affect cellular behavior, and in particular the ability to directly test existing theories about cell state and variation in the context of morphology, location, motility, and signaling that could not be tested before. Rapid advancements in resolution, ease-of-use, and scale of spatial genomics technologies to address these questions also require an updated toolkit of statistical methods with which to interrogate these data. We present a framework to respond to this new avenue of research: four open biological questions that can now be answered using spatial genomics data paired with methods for analysis. We outline spatial data modalities for each open question that may yield specific insights, discuss how conflicting theories may be tested by comparing the data to conceptual models of biological behavior, and highlight statistical and machine learning-based tools that may prove particularly helpful to recover biological understanding.

Detection of allele-specific expression in spatial transcriptomics with spASE

Spatial transcriptomics technologies permit the study of the spatial distribution of RNA at near-single-cell resolution genome-wide. However, the feasibility of studying spatial allele-specific expression (ASE) from these data remains uncharacterized. Here, we introduce spASE, a computational framework for detecting and estimating spatial ASE. To tackle the challenges presented by cell type mixtures and a low signal to noise ratio, we implement a hierarchical model involving additive mixtures of spatial smoothing splines. We apply our method to allele-resolved Visium and Slide-seq from the mouse cerebellum and hippocampus and report new insight into the landscape of spatial and cell type-specific ASE therein.

Deciphering brain cellular and behavioral mechanisms: Insights from single-cell and spatial RNA sequencing

The brain is a complex computing system composed of a multitude of interacting neurons. The computational outputs of this system determine the behavior and perception of every individual. Each brain cell expresses thousands of genes that dictate the cell's function and physiological properties. Therefore, deciphering the molecular expression of each cell is of great significance for understanding its characteristics and role in brain function. Additionally, the positional information of each cell can provide crucial insights into their involvement in local brain circuits. In this review, we briefly overview the principles of single-cell RNA sequencing and spatial transcriptomics, the potential issues and challenges in their data processing, and their applications in brain research. We further outline several promising directions in neuroscience that could be integrated with single-cell RNA sequencing, including neurodevelopment, the identification of novel brain microstructures, cognition and behavior, neuronal cell positioning, molecules and cells related to advanced brain functions, sleep-wake cycles/circadian rhythms, and computational modeling of brain function. We believe that the deep integration of these directions with single-cell and spatial RNA sequencing can contribute significantly to understanding the roles of individual cells or cell types in these specific functions, thereby making important contributions to addressing critical questions in those fields. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Development RNA in Disease and Development > RNA in Disease.

Integrating Single-Cell and Spatial Transcriptomics to Uncover and Elucidate GP73-Mediated Pro-Angiogenic Regulatory Networks in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) was characterized as being hypervascular. In the present study, we generated a single-cell spatial transcriptomic landscape of the vasculogenic etiology of HCC and illustrated overexpressed Golgi phosphoprotein 73 (GP73) HCC cells exerting cellular communication with vascular endothelial cells with high pro-angiogenesis potential via multiple receptor-ligand interactions in the process of tumor vascular development. Specifically, we uncovered an interactive GP73-mediated regulatory network coordinated with c-Myc, lactate, Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/STAT3) pathway, and endoplasmic reticulum stress (ERS) signals in HCC cells and elucidated its pro-angiogenic roles in vitro and in vivo. Mechanistically, we found that GP73, the pivotal hub gene, was activated by histone lactylation and c-Myc, which stimulated the phosphorylation of downstream STAT3 by directly binding STAT3 and simultaneously enhancing glucose-regulated protein 78 (GRP78)-induced ERS. STAT3 potentiates GP73-mediated pro-angiogenic functions. Clinically, serum GP73 levels were positively correlated with HCC response to anti-angiogenic regimens and were essential for a prognostic nomogram showing good predictive performance for determining 6-month and 1-year survival in patients with HCC treated with anti-angiogenic therapy. Taken together, the aforementioned data characterized the pro-angiogenic roles and mechanisms of a GP73-mediated network and proved that GP73 is a crucial tumor angiogenesis niche gene with favorable anti-angiogenic potential in the treatment of HCC.

Unravelling infiltrating T-cell heterogeneity in kidney renal clear cell carcinoma: Integrative single-cell and spatial transcriptomic profiling

Kidney renal clear cell carcinoma (KIRC) pathogenesis intricately involves immune system dynamics, particularly the role of T cells within the tumour microenvironment. Through a multifaceted approach encompassing single-cell RNA sequencing, spatial transcriptome analysis and bulk transcriptome profiling, we systematically explored the contribution of infiltrating T cells to KIRC heterogeneity. Employing high-density weighted gene co-expression network analysis (hdWGCNA), module scoring and machine learning, we identified a distinct signature of infiltrating T cell-associated genes (ITSGs). Spatial transcriptomic data were analysed using robust cell type decomposition (RCTD) to uncover spatial interactions. Further analyses included enrichment assessments, immune infiltration evaluations and drug susceptibility predictions. Experimental validation involved PCR experiments, CCK-8 assays, plate cloning assays, wound-healing assays and Transwell assays. Six subpopulations of infiltrating and proliferating T cells were identified in KIRC, with notable dynamics observed in mid- to late-stage disease progression. Spatial analysis revealed significant correlations between T cells and epithelial cells across varying distances within the tumour microenvironment. The ITSG-based prognostic model demonstrated robust predictive capabilities, implicating these genes in immune modulation and metabolic pathways and offering prognostic insights into drug sensitivity for 12 KIRC treatment agents. Experimental validation underscored the functional relevance of PPIB in KIRC cell proliferation, invasion and migration. Our study comprehensively characterizes infiltrating T-cell heterogeneity in KIRC using single-cell RNA sequencing and spatial transcriptome data. The stable prognostic model based on ITSGs unveils infiltrating T cells' prognostic potential, shedding light on the immune microenvironment and offering avenues for personalized treatment and immunotherapy.

A Review of the Application of Spatial Transcriptomics in Neuroscience

Since spatial transcriptomics can locate and distinguish the gene expression of functional genes in special regions and tissue, it is important for us to investigate the brain development, the development mechanism of brain diseases, and the relationship between brain structure and function in Neuroscience (or Brain science). While previous studies have introduced the crucial spatial transcriptomic techniques and data analysis methods, there are few studies to comprehensively overview the key methods, data resources, and technological applications of spatial transcriptomics in Neuroscience. For these reasons, we first investigate several common spatial transcriptomic data analysis approaches and data resources. Second, we introduce the applications of the spatial transcriptomic data analysis approaches in Neuroscience. Third, we summarize the integrating spatial transcriptomics with other technologies in Neuroscience. Finally, we discuss the challenges and future research directions of spatial transcriptomics in Neuroscience.

Spatial Multi-Omics in Alzheimer's Disease: A Multi-Dimensional Approach to Understanding Pathology and Progression

Alzheimer's Disease (AD) presents a complex neuropathological landscape characterized by hallmark amyloid plaques and neurofibrillary tangles, leading to progressive cognitive decline. Despite extensive research, the molecular intricacies contributing to AD pathogenesis are inadequately understood. While single-cell omics technology holds great promise for application in AD, particularly in deciphering the understanding of different cell types and analyzing rare cell types and transcriptomic expression changes, it is unable to provide spatial distribution information, which is crucial for understanding the pathological processes of AD. In contrast, spatial multi-omics research emerges as a promising and comprehensive approach to analyzing tissue cells, potentially better suited for addressing these issues in AD. This article focuses on the latest advancements in spatial multi-omics technology and compares various techniques. Additionally, we provide an overview of current spatial omics-based research results in AD. These technologies play a crucial role in facilitating new discoveries and advancing translational AD research in the future. Despite challenges such as balancing resolution, increasing throughput, and data analysis, the application of spatial multi-omics holds immense potential in revolutionizing our understanding of human disease processes and identifying new biomarkers and therapeutic targets, thereby potentially contributing to the advancement of AD research.

spVC for the detection and interpretation of spatial gene expression variation

Spatially resolved transcriptomics technologies have opened new avenues for understanding gene expression heterogeneity in spatial contexts. However, existing methods for identifying spatially variable genes often focus solely on statistical significance, limiting their ability to capture continuous expression patterns and integrate spot-level covariates. To address these challenges, we introduce spVC, a statistical method based on a generalized Poisson model. spVC seamlessly integrates constant and spatially varying effects of covariates, facilitating comprehensive exploration of gene expression variability and enhancing interpretability. Simulation and real data applications confirm spVC's accuracy in these tasks, highlighting its versatility in spatial transcriptomics analysis.

Benchmarking spatial clustering methods with spatially resolved transcriptomics data

Spatial clustering, which shares an analogy with single-cell clustering, has expanded the scope of tissue physiology studies from cell-centroid to structure-centroid with spatially resolved transcriptomics (SRT) data. Computational methods have undergone remarkable development in recent years, but a comprehensive benchmark study is still lacking. Here we present a benchmark study of 13 computational methods on 34 SRT data (7 datasets). The performance was evaluated on the basis of accuracy, spatial continuity, marker genes detection, scalability, and robustness. We found existing methods were complementary in terms of their performance and functionality, and we provide guidance for selecting appropriate methods for given scenarios. On testing additional 22 challenging datasets, we identified challenges in identifying noncontinuous spatial domains and limitations of existing methods, highlighting their inadequacies in handling recent large-scale tasks. Furthermore, with 145 simulated data, we examined the robustness of these methods against four different factors, and assessed the impact of pre- and postprocessing approaches. Our study offers a comprehensive evaluation of existing spatial clustering methods with SRT data, paving the way for future advancements in this rapidly evolving field.

Regional analysis to delineate intrasample heterogeneity with RegionalST

MOTIVATION

Spatial transcriptomics has greatly contributed to our understanding of spatial and intra-sample heterogeneity, which could be crucial for deciphering the molecular basis of human diseases. Intra-tumor heterogeneity, e.g. may be associated with cancer treatment responses. However, the lack of computational tools for exploiting cross-regional information and the limited spatial resolution of current technologies present major obstacles to elucidating tissue heterogeneity.

RESULTS

To address these challenges, we introduce RegionalST, an efficient computational method that enables users to quantify cell type mixture and interactions, identify sub-regions of interest, and perform cross-region cell type-specific differential analysis for the first time. Our simulations and real data applications demonstrate that RegionalST is an efficient tool for visualizing and analyzing diverse spatial transcriptomics data, thereby enabling accurate and flexible exploration of tissue heterogeneity. Overall, RegionalST provides a one-stop destination for researchers seeking to delve deeper into the intricacies of spatial transcriptomics data.

AVAILABILITY AND IMPLEMENTATION

The implementation of our method is available as an open-source R/Bioconductor package with a user-friendly manual available at https://bioconductor.org/packages/release/bioc/html/RegionalST.html.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

Explore & actuate: the future of personalized medicine in oncology through emerging technologies

PURPOSE OF REVIEW

The future of medicine is aimed to equip the physician with tools to assess the individual health of the patient for the uniqueness of the disease that separates it from the rest. The integration of omics technologies into clinical practice, reviewed here, would open new avenues for addressing the spatial and temporal heterogeneity of cancer. The rising cancer burden patiently awaits the advent of such an approach to personalized medicine for routine clinical settings.

RECENT FINDINGS

To weigh the translational potential, multiple technologies were categorized based on the extractable information from the different types of samples used, to the various omic-levels of molecular information that each technology has been able to advance over the last 2 years. This review uses a multifaceted classification that helps to assess translational potential in a meaningful way toward clinical adaptation.

SUMMARY

The importance of distinguishing technologies based on the flow of information from exploration to actuation puts forth a framework that allows the clinicians to better adapt a chosen technology or use them in combination to enhance their goals toward personalized medicine.

Challenges and perspectives in computational deconvolution of genomics data

Deciphering cell-type heterogeneity is crucial for systematically understanding tissue homeostasis and its dysregulation in diseases. Computational deconvolution is an efficient approach for estimating cell-type abundances from a variety of omics data. Despite substantial methodological progress in computational deconvolution in recent years, challenges are still outstanding. Here we enlist four important challenges related to computational deconvolution: the quality of the reference data, generation of ground truth data, limitations of computational methodologies, and benchmarking design and implementation. Finally, we make recommendations on reference data generation, new directions of computational methodologies, and strategies to promote rigorous benchmarking.

Streamlining spatial omics data analysis with Pysodb

Advances in spatial omics technologies have improved the understanding of cellular organization in tissues, leading to the generation of complex and heterogeneous data and prompting the development of specialized tools for managing, loading and visualizing spatial omics data. The Spatial Omics Database (SODB) was established to offer a unified format for data storage and interactive visualization modules. Here we detail the use of Pysodb, a Python-based tool designed to enable the efficient exploration and loading of spatial datasets from SODB within a Python environment. We present seven case studies using Pysodb, detailing the interaction with various computational methods, ensuring reproducibility of experimental data and facilitating the integration of new data and alternative applications in SODB. The approach offers a reference for method developers by outlining label and metadata availability in representative spatial data that can be loaded by Pysodb. The tool is supplemented by a website ( https://protocols-pysodb.readthedocs.io/ ) with detailed information for benchmarking analysis, and allows method developers to focus on computational models by facilitating data processing. This protocol is designed for researchers with limited experience in computational biology. Depending on the dataset complexity, the protocol typically requires ~12 h to complete.

Niche-DE: niche-differential gene expression analysis in spatial transcriptomics data identifies context-dependent cell-cell interactions

Existing methods for analysis of spatial transcriptomic data focus on delineating the global gene expression variations of cell types across the tissue, rather than local gene expression changes driven by cell-cell interactions. We propose a new statistical procedure called niche-differential expression (niche-DE) analysis that identifies cell-type-specific niche-associated genes, which are differentially expressed within a specific cell type in the context of specific spatial niches. We further develop niche-LR, a method to reveal ligand-receptor signaling mechanisms that underlie niche-differential gene expression patterns. Niche-DE and niche-LR are applicable to low-resolution spot-based spatial transcriptomics data and data that is single-cell or subcellular in resolution.

Myeloid cell-derived creatine in the hypoxic niche promotes glioblastoma growth

Glioblastoma (GBM) is a malignancy dominated by the infiltration of tumor-associated myeloid cells (TAMCs). Examination of TAMC metabolic phenotypes in mouse models and patients with GBM identified the de novo creatine metabolic pathway as a hallmark of TAMCs. Multi-omics analyses revealed that TAMCs surround the hypoxic peri-necrotic regions of GBM and express the creatine metabolic enzyme glycine amidinotransferase (GATM). Conversely, GBM cells located within these same regions are uniquely specific in expressing the creatine transporter (SLC6A8). We hypothesized that TAMCs provide creatine to tumors, promoting GBM progression. Isotopic tracing demonstrated that TAMC-secreted creatine is taken up by tumor cells. Creatine supplementation protected tumors from hypoxia-induced stress, which was abrogated with genetic ablation or pharmacologic inhibition of SLC6A8. Lastly, inhibition of creatine transport using the clinically relevant compound, RGX-202-01, blunted tumor growth and enhanced radiation therapy in vivo. This work highlights that myeloid-to-tumor transfer of creatine promotes tumor growth in the hypoxic niche.

Mapping the transcriptome: Realizing the full potential of spatial data analysis

RNA sequencing in situ allows for whole-transcriptome characterization at high resolution, while retaining spatial information. These data present an analytical challenge for bioinformatics-how to leverage spatial information effectively? Properties of data with a spatial dimension require special handling, which necessitate a different set of statistical and inferential considerations when compared to non-spatial data. The geographical sciences primarily use spatial data and have developed methods to analye them. Here we discuss the challenges associated with spatial analysis and examine how we can take advantage of practice from the geographical sciences to realize the full potential of spatial information in transcriptomic datasets.

scResolve: Recovering single cell expression profiles from multi-cellular spatial transcriptomics

Many popular spatial transcriptomics techniques lack single-cell resolution. Instead, these methods measure the collective gene expression for each location from a mixture of cells, potentially containing multiple cell types. Here, we developed scResolve, a method for recovering single-cell expression profiles from spatial transcriptomics measurements at multi-cellular resolution. scResolve accurately restores expression profiles of individual cells at their locations, which is unattainable from cell type deconvolution. Applications of scResolve on human breast cancer data and human lung disease data demonstrate that scResolve enables cell type-specific differential gene expression analysis between different tissue contexts and accurate identification of rare cell populations. The spatially resolved cellular-level expression profiles obtained through scResolve facilitate more flexible and precise spatial analysis that complements raw multi-cellular level analysis.

Reconstructing Spatial Transcriptomics at the Single-cell Resolution with BayesDeep

Spatially resolved transcriptomics (SRT) techniques have revolutionized the characterization of molecular profiles while preserving spatial and morphological context. However, most next-generation sequencing-based SRT techniques are limited to measuring gene expression in a confined array of spots, capturing only a fraction of the spatial domain. Typically, these spots encompass gene expression from a few to hundreds of cells, underscoring a critical need for more detailed, single-cell resolution SRT data to enhance our understanding of biological functions within the tissue context. Addressing this challenge, we introduce BayesDeep, a novel Bayesian hierarchical model that leverages cellular morphological data from histology images, commonly paired with SRT data, to reconstruct SRT data at the single-cell resolution. BayesDeep effectively model count data from SRT studies via a negative binomial regression model. This model incorporates explanatory variables such as cell types and nuclei-shape information for each cell extracted from the paired histology image. A feature selection scheme is integrated to examine the association between the morphological and molecular profiles, thereby improving the model robustness. We applied BayesDeep to two real SRT datasets, successfully demonstrating its capability to reconstruct SRT data at the single-cell resolution. This advancement not only yields new biological insights but also significantly enhances various downstream analyses, such as pseudotime and cell-cell communication.

Single cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation

INTRODUCTION

The R47H missense mutation of the TREM2 gene is a strong risk factor for development of Alzheimer's Disease. We investigate cell-type-specific spatial transcriptomic changes induced by the Trem2R47H mutation to determine the impacts of this mutation on transcriptional dysregulation.

METHODS

We profiled 15 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics. Single-cell spatial transcriptomics and neuropathology data were analyzed using our custom pipeline to identify plaque and Trem2R47H induced transcriptomic dysregulation.

RESULTS

The Trem2R47H mutation induced consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice.

CONCLUSION

Spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

Charting the tumor microenvironment with spatial profiling technologies

In recent years technologies that can achieve readouts at cellular resolution such as single-cell RNA sequencing (scRNA-seq) have provided a comprehensive characterization of the cellular proportions and phenotypes that populate the tumor microenvironment (TME). However, because of the sample dissociation steps required by these protocols, they fail to capture information related to the intricate spatial context in which cells operate as well as their dense networks of interactions. Spatial profiling technologies have recently emerged as a valuable way to investigate the physical organization of cells crowding the TME in intact tissues. In this review we first discuss how spatial profiling technologies have propelled TME characterization, and then explore their potential to improve both diagnosis and prognosis for cancer patients in the clinic.

Integrating spatial and single-cell transcriptomics data using deep generative models with SpatialScope

The rapid emergence of spatial transcriptomics (ST) technologies is revolutionizing our understanding of tissue spatial architecture and biology. Although current ST methods, whether based on next-generation sequencing (seq-based approaches) or fluorescence in situ hybridization (image-based approaches), offer valuable insights, they face limitations either in cellular resolution or transcriptome-wide profiling. To address these limitations, we present SpatialScope, a unified approach integrating scRNA-seq reference data and ST data using deep generative models. With innovation in model and algorithm designs, SpatialScope not only enhances seq-based ST data to achieve single-cell resolution, but also accurately infers transcriptome-wide expression levels for image-based ST data. We demonstrate SpatialScope's utility through simulation studies and real data analysis from both seq-based and image-based ST approaches. SpatialScope provides spatial characterization of tissue structures at transcriptome-wide single-cell resolution, facilitating downstream analysis, including detecting cellular communication through ligand-receptor interactions, localizing cellular subtypes, and identifying spatially differentially expressed genes.

STAMarker: determining spatial domain-specific variable genes with saliency maps in deep learning

Spatial transcriptomics characterizes gene expression profiles while retaining the information of the spatial context, providing an unprecedented opportunity to understand cellular systems. One of the essential tasks in such data analysis is to determine spatially variable genes (SVGs), which demonstrate spatial expression patterns. Existing methods only consider genes individually and fail to model the inter-dependence of genes. To this end, we present an analytic tool STAMarker for robustly determining spatial domain-specific SVGs with saliency maps in deep learning. STAMarker is a three-stage ensemble framework consisting of graph-attention autoencoders, multilayer perceptron (MLP) classifiers, and saliency map computation by the backpropagated gradient. We illustrate the effectiveness of STAMarker and compare it with serveral commonly used competing methods on various spatial transcriptomic data generated by different platforms. STAMarker considers all genes at once and is more robust when the dataset is very sparse. STAMarker could identify spatial domain-specific SVGs for characterizing spatial domains and enable in-depth analysis of the region of interest in the tissue section.

CellNeighborEX: deciphering neighbor-dependent gene expression from spatial transcriptomics data

Cells have evolved their communication methods to sense their microenvironments and send biological signals. In addition to communication using ligands and receptors, cells use diverse channels including gap junctions to communicate with their immediate neighbors. Current approaches, however, cannot effectively capture the influence of various microenvironments. Here, we propose a novel approach to investigate cell neighbor-dependent gene expression (CellNeighborEX) in spatial transcriptomics (ST) data. To categorize cells based on their microenvironment, CellNeighborEX uses direct cell location or the mixture of transcriptome from multiple cells depending on ST technologies. For each cell type, CellNeighborEX identifies diverse gene sets associated with partnering cell types, providing further insight. We found that cells express different genes depending on their neighboring cell types in various tissues including mouse embryos, brain, and liver cancer. Those genes are associated with critical biological processes such as development or metastases. We further validated that gene expression is induced by neighboring partners via spatial visualization. The neighbor-dependent gene expression suggests new potential genes involved in cell-cell interactions beyond what ligand-receptor co-expression can discover.

Ataxia-telangiectasia mutated ( Atm ) disruption sensitizes spatially-directed H3.3K27M/TP53 diffuse midline gliomas to radiation therapy

Diffuse midline gliomas (DMGs) are lethal brain tumors characterized by p53-inactivating mutations and oncohistone H3.3K27M mutations that rewire the cellular response to genotoxic stress, which presents therapeutic opportunities. We used RCAS/tv-a retroviruses and Cre recombinase to inactivate p53 and induce K27M in the native H3f3a allele in a lineage- and spatially-directed manner, yielding primary mouse DMGs. Genetic or pharmacologic disruption of the DNA damage response kinase Ataxia-telangiectasia mutated (ATM) enhanced the efficacy of focal brain irradiation, extending mouse survival. This finding suggests that targeting ATM will enhance the efficacy of radiation therapy for p53-mutant DMG but not p53-wildtype DMG. We used spatial in situ transcriptomics and an allelic series of primary murine DMG models with different p53 mutations to identify transactivation-independent p53 activity as a key mediator of such radiosensitivity. These studies deeply profile a genetically faithful and versatile model of a lethal brain tumor to identify resistance mechanisms for a therapeutic strategy currently in clinical trials.