The dawn of spatial omics

Computational methods for alignment and integration of spatially resolved transcriptomics data

Most of the complex biological regulatory activities occur in three dimensions (3D). To better analyze biological processes, it is essential not only to decipher the molecular information of numerous cells but also to understand how their spatial contexts influence their behavior. With the development of spatially resolved transcriptomics (SRT) technologies, SRT datasets are being generated to simultaneously characterize gene expression and spatial arrangement information within tissues, organs or organisms. To fully leverage spatial information, the focus extends beyond individual two-dimensional (2D) slices. Two tasks known as slices alignment and data integration have been introduced to establish correlations between multiple slices, enhancing the effectiveness of downstream tasks. Currently, numerous related methods have been developed. In this review, we first elucidate the details and principles behind several representative methods. Then we report the testing results of these methods on various SRT datasets, and assess their performance in representative downstream tasks. Insights into the strengths and weaknesses of each method and the reasons behind their performance are discussed. Finally, we provide an outlook on future developments. The codes and details of experiments are now publicly available at https://github.com/YangLabHKUST/SRT_alignment_and_integration.

Single cell spatial biology over developmental time can decipher pediatric brain pathologies

Pediatric low grade brain tumors and neurodevelopmental disorders share proteins, signaling pathways, and networks. They also share germline mutations and an impaired prenatal differentiation origin. They may differ in the timing of the events and proliferation. We suggest that their pivotal distinct, albeit partially overlapping, outcomes relate to the cell states, which depend on their spatial location, and timing of gene expression during brain development. These attributes are crucial as the brain develops sequentially, and single-cell spatial organization influences cell state, thus function. Our underlying premise is that the root cause in neurodevelopmental disorders and pediatric tumors is impaired prenatal differentiation. Data related to pediatric brain tumors, neurodevelopmental disorders, brain cell (sub)types, locations, and timing of expression in the developing brain are scant. However, emerging single cell technologies, including transcriptomic, spatial biology, spatial high-resolution imaging performed over the brain developmental time, could be transformational in deciphering brain pathologies thereby pharmacology.

Decoding the spatiotemporal heterogeneity of tumor-associated macrophages

Tumor-associated macrophages (TAMs) are pivotal in cancer progression, influencing tumor growth, angiogenesis, and immune evasion. This review explores the spatial and temporal heterogeneity of TAMs within the tumor microenvironment (TME), highlighting their diverse subtypes, origins, and functions. Advanced technologies such as single-cell sequencing and spatial multi-omics have elucidated the intricate interactions between TAMs and other TME components, revealing the mechanisms behind their recruitment, polarization, and distribution. Key findings demonstrate that TAMs support tumor vascularization, promote epithelial-mesenchymal transition (EMT), and modulate extracellular matrix (ECM) remodeling, etc., thereby enhancing tumor invasiveness and metastasis. Understanding these complex dynamics offers new therapeutic targets for disrupting TAM-mediated pathways and overcoming drug resistance. This review underscores the potential of targeting TAMs to develop innovative cancer therapies, emphasizing the need for further research into their spatial characteristics and functional roles within the TME.

Open-ST: High-resolution spatial transcriptomics in 3D

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

Integrated multi-omics with machine learning to uncover the intricacies of kidney disease

The development of omics technologies has driven a profound expansion in the scale of biological data and the increased complexity in internal dimensions, prompting the utilization of machine learning (ML) as a powerful toolkit for extracting knowledge and understanding underlying biological patterns. Kidney disease represents one of the major growing global health threats with intricate pathogenic mechanisms and a lack of precise molecular pathology-based therapeutic modalities. Accordingly, there is a need for advanced high-throughput approaches to capture implicit molecular features and complement current experiments and statistics. This review aims to delineate strategies for integrating multi-omics data with appropriate ML methods, highlighting key clinical translational scenarios, including predicting disease progression risks to improve medical decision-making, comprehensively understanding disease molecular mechanisms, and practical applications of image recognition in renal digital pathology. Examining the benefits and challenges of current integration efforts is expected to shed light on the complexity of kidney disease and advance clinical practice.

Spatial landscapes of cancers: insights and opportunities

Solid tumours comprise many different cell types organized in spatially structured arrangements, with substantial intratumour and intertumour heterogeneity. Advances in spatial profiling technologies over the past decade hold promise to capture the complexity of these cellular architectures to build a holistic view of the intricate molecular mechanisms that shape the tumour ecosystem. Some of these mechanisms act at the cellular scale and are controlled by cell-autonomous programmes or communication between nearby cells, whereas other mechanisms result from coordinated efforts between large networks of cells and extracellular molecules organized into tissues and organs. In this Review we provide insights into the application of single-cell and spatial profiling tools, with a focus on spatially resolved transcriptomic tools developed to understand the cellular architecture of the tumour microenvironment and identify opportunities to use them to improve clinical management of cancers.

Spatial omics technologies for understanding molecular status associated with cancer progression

Cancer cells are generally exposed to numerous extrinsic stimulations in the tumor microenvironment. In this environment, cancer cells change their expression profiles to fight against circumstantial stresses, allowing their progression in the challenging tissue space. Technological advancements of spatial omics have had substantial influence on cancer genomics. This technical progress, especially that occurring in the spatial transcriptome, has been drastic and rapid. Here, we describe the latest spatial analytical technologies that have allowed omics feature characterization to retain their spatial and histopathological information in cancer tissues. Several spatial omics platforms have been launched, and the latest platforms finally attained single-cell level or even higher subcellular level resolution. We discuss several key papers elucidating the initial utility of the spatial analysis. In fact, spatial transcriptome analyses reveal comprehensive omics characteristics not only in cancer cells but also their surrounding cells, such as tumor infiltrating immune cells and cancer-associated fibroblasts. We also introduce several spatial omics platforms. We describe our own attempts to investigate molecular events associated with cancer progression. Furthermore, we discuss the next challenges in analyzing the multiomics status of cells, including their morphology and location. These novel technologies, in conjunction with spatial transcriptome analysis and, more importantly, with histopathology, will elucidate even novel key aspects of the intratumor heterogeneity of cancers. Such enhanced knowledge is expected to open a new path for overcoming therapeutic resistance and eventually to precisely stratify patients.

SANTO: a coarse-to-fine alignment and stitching method for spatial omics

With the flourishing of spatial omics technologies, alignment and stitching of slices becomes indispensable to decipher a holistic view of 3D molecular profile. However, existing alignment and stitching methods are unpractical to process large-scale and image-based spatial omics dataset due to extreme time consumption and unsatisfactory accuracy. Here we propose SANTO, a coarse-to-fine method targeting alignment and stitching tasks for spatial omics. SANTO firstly rapidly supplies reasonable spatial positions of two slices and identifies the overlap region. Then, SANTO refines the positions of two slices by considering spatial and omics patterns. Comprehensive experiments demonstrate the superior performance of SANTO over existing methods. Specifically, SANTO stitches cross-platform slices for breast cancer samples, enabling integration of complementary features to synergistically explore tumor microenvironment. SANTO is then applied to 3D-to-3D spatiotemporal alignment to study development of mouse embryo. Furthermore, SANTO enables cross-modality alignment of spatial transcriptomic and epigenomic data to understand complementary interactions.

Discovering cancer stem-like molecule, nuclear factor I X, using spatial transcriptome in gastric cancer

Gastric cancer (GC) is characterized by significant intratumoral heterogeneity, and stem cells are promising therapeutic targets. Despite advancements in spatial transcriptome analyses, unexplored targets for addressing cancer stemness remain unknown. This study aimed to identify Nuclear Factor IX (NFIX) as a critical regulator of cancer stemness in GC and evaluate its clinicopathological significance and function. Spatial transcriptome analysis of GC was conducted. The correlation between NFIX expression, clinicopathological factors, and prognosis was assessed using immunostaining in 127 GC cases. Functional analyses of cancer cell lines validated these findings. Spatial transcriptome analysis stratified GC tissues based on genetic profiles, identified CSC-like cells, and further refined the classification to identify and highlight the significance of NFIX, as validated by Monocle 3 and CytoTRACE analyses. Knockdown experiments in cancer cell lines have demonstrated the involvement of NFIX in cancer cell proliferation and kinase activity. This study underscores the role of spatial transcriptome analysis in refining GC tissue classification and identifying therapeutic targets, highlighting NFIX as a pivotal factor. NFIX expression is correlated with poor prognosis and drives GC progression, suggesting its potential as a novel therapeutic target for personalized GC therapies.

Mass spectrometry imaging for spatially resolved multi-omics molecular mapping

The recent upswing in the integration of spatial multi-omics for conducting multidimensional information measurements is opening a new chapter in biological research. Mapping the landscape of various biomolecules including metabolites, proteins, nucleic acids, etc., and even deciphering their functional interactions and pathways is believed to provide a more holistic and nuanced exploration of the molecular intricacies within living systems. Mass spectrometry imaging (MSI) stands as a forefront technique for spatially mapping the metabolome, lipidome, and proteome within diverse tissue and cell samples. In this review, we offer a systematic survey delineating different MSI techniques for spatially resolved multi-omics analysis, elucidating their principles, capabilities, and limitations. Particularly, we focus on the advancements in methodologies aimed at augmenting the molecular sensitivity and specificity of MSI; and depict the burgeoning integration of MSI-based spatial metabolomics, lipidomics, and proteomics, encompassing the synergy with other imaging modalities. Furthermore, we offer speculative insights into the potential trajectory of MSI technology in the future.

FUSION: A web-based application for in-depth exploration of multi-omics data with brightfield histology

Spatial -OMICS technologies facilitate the interrogation of molecular profiles in the context of the underlying histopathology and tissue microenvironment. Paired analysis of histopathology and molecular data can provide pathologists with otherwise unobtainable insights into biological mechanisms. To connect the disparate molecular and histopathologic features into a single workspace, we developed FUSION (Functional Unit State IdentificatiON in WSIs [Whole Slide Images]), a web-based tool that provides users with a broad array of visualization and analytical tools including deep learning-based algorithms for in-depth interrogation of spatial -OMICS datasets and their associated high-resolution histology images. FUSION enables end-to-end analysis of functional tissue units (FTUs), automatically aggregating underlying molecular data to provide a histopathology-based medium for analyzing healthy and altered cell states and driving new discoveries using "pathomic" features. We demonstrate FUSION using 10x Visium spatial transcriptomics (ST) data from both formalin-fixed paraffin embedded (FFPE) and frozen prepared datasets consisting of healthy and diseased tissue. Through several use-cases, we demonstrate how users can identify spatial linkages between quantitative pathomics, qualitative image characteristics, and spatial --omics.

Advances in mass spectrometry imaging for plant metabolomics-Expanding the analytical toolbox

Mass spectrometry imaging (MSI) has become increasingly popular in plant science due to its ability to characterize complex chemical, spatial, and temporal aspects of plant metabolism. Over the past decade, as the emerging and unique features of various MSI techniques have continued to support new discoveries in studies of plant metabolism closely associated with various aspects of plant function and physiology, spatial metabolomics based on MSI techniques has positioned it at the forefront of plant metabolic studies, providing the opportunity for far higher resolution than was previously available. Despite these efforts, profound challenges at the levels of spatial resolution, sensitivity, quantitative ability, chemical confidence, isomer discrimination, and spatial multi-omics integration, undoubtedly remain. In this Perspective, we provide a contemporary overview of the emergent MSI techniques widely used in the plant sciences, with particular emphasis on recent advances in methodological breakthroughs. Having established the detailed context of MSI, we outline both the golden opportunities and key challenges currently facing plant metabolomics, presenting our vision as to how the enormous potential of MSI technologies will contribute to progress in plant science in the coming years.

Multiplexed imaging to reveal tissue dendritic cell spatial localisation and function

Dendritic cells (DCs) play a pivotal role in immune surveillance, acting as sentinels that coordinate immune responses within tissues. Although differences in the identity and functional states of DC subpopulations have been identified through multiparametric flow cytometry and single-cell RNA sequencing, these methods do not provide information about the spatial context in which the cells are located. This knowledge is crucial for understanding tissue organisation and cellular cross-talk. Recent developments in multiplex imaging techniques can now offer insights into this complex spatial and functional landscape. This review provides a concise overview of these imaging methodologies, emphasising their application in identifying DCs to delineate their tissue-specific functions and aiding newcomers in navigating this field.

Points2Regions: Fast, interactive clustering of imaging-based spatial transcriptomics data

Imaging-based spatial transcriptomics techniques generate data in the form of spatial points belonging to different mRNA classes. A crucial part of analyzing the data involves the identification of regions with similar composition of mRNA classes. These biologically interesting regions can manifest at different spatial scales. For example, the composition of mRNA classes on a cellular scale corresponds to cell types, whereas compositions on a millimeter scale correspond to tissue-level structures. Traditional techniques for identifying such regions often rely on complementary data, such as pre-segmented cells, or lengthy optimization. This limits their applicability to tasks on a particular scale, restricting their capabilities in exploratory analysis. This article introduces "Points2Regions," a computational tool for identifying regions with similar mRNA compositions. The tool's novelty lies in its rapid feature extraction by rasterizing points (representing mRNAs) onto a pyramidal grid and its efficient clustering using a combination of hierarchical andk $$ k $$ -means clustering. This enables fast and efficient region discovery across multiple scales without relying on additional data, making it a valuable resource for exploratory analysis. Points2Regions has demonstrated performance similar to state-of-the-art methods on two simulated datasets, without relying on segmented cells, while being several times faster. Experiments on real-world datasets show that regions identified by Points2Regions are similar to those identified in other studies, confirming that Points2Regions can be used to extract biologically relevant regions. The tool is shared as a Python package integrated into TissUUmaps and a Napari plugin, offering interactive clustering and visualization, significantly enhancing user experience in data exploration.

Deep 3D histology powered by tissue clearing, omics and AI

To comprehensively understand tissue and organism physiology and pathophysiology, it is essential to create complete three-dimensional (3D) cellular maps. These maps require structural data, such as the 3D configuration and positioning of tissues and cells, and molecular data on the constitution of each cell, spanning from the DNA sequence to protein expression. While single-cell transcriptomics is illuminating the cellular and molecular diversity across species and tissues, the 3D spatial context of these molecular data is often overlooked. Here, I discuss emerging 3D tissue histology techniques that add the missing third spatial dimension to biomedical research. Through innovations in tissue-clearing chemistry, labeling and volumetric imaging that enhance 3D reconstructions and their synergy with molecular techniques, these technologies will provide detailed blueprints of entire organs or organisms at the cellular level. Machine learning, especially deep learning, will be essential for extracting meaningful insights from the vast data. Further development of integrated structural, molecular and computational methods will unlock the full potential of next-generation 3D histology.

SEraster: a rasterization preprocessing framework for scalable spatial omics data analysis

MOTIVATION

Spatial omics data demand computational analysis but many analysis tools have computational resource requirements that increase with the number of cells analyzed. This presents scalability challenges as researchers use spatial omics technologies to profile millions of cells.

RESULTS

To enhance the scalability of spatial omics data analysis, we developed a rasterization preprocessing framework called SEraster that aggregates cellular information into spatial pixels. We apply SEraster to both real and simulated spatial omics data prior to spatial variable gene expression analysis to demonstrate that such preprocessing can reduce computational resource requirements while maintaining high performance, including as compared to other down-sampling approaches. We further integrate SEraster with existing analysis tools to characterize cell-type spatial co-enrichment across length scales. Finally, we apply SEraster to enable analysis of a mouse pup spatial omics dataset with over a million cells to identify tissue-level and cell-type-specific spatially variable genes as well as spatially co-enriched cell types that recapitulate expected organ structures.

AVAILABILITY AND IMPLEMENTATION

SEraster is implemented as an R package on GitHub (https://github.com/JEFworks-Lab/SEraster) with additional tutorials at https://JEF.works/SEraster.

Exploring kidney allograft rejection: A proof-of-concept study using spatial transcriptomics

In this proof-of-concept study, spatial transcriptomics combined with public single-cell ribonucleic acid-sequencing data were used to explore the potential of this technology to study kidney allograft rejection. We aimed to map gene expression patterns within diverse pathologic states by examining biopsies classified across nonrejection, T cell-mediated acute rejection, interstitial fibrosis, and tubular atrophy. Our results revealed distinct immune cell signatures, including those of T and B lymphocytes, monocytes, mast cells, and plasma cells, and their spatial organization within the renal interstitium. We also mapped chemokine receptors and ligands to study immune cell migration and recruitment. Finally, our analysis demonstrated differential spatial enrichment of transcription signatures associated with kidney allograft rejection across various biopsy regions. Interstitium regions displayed higher enrichment scores for rejection-associated gene expression patterns than tubular areas, which had negative scores. This implies that these signatures are primarily driven by processes unfolding in the renal interstitium. Overall, this study highlights the value of spatial transcriptomics for revealing cellular heterogeneity and immune signatures in renal transplant biopsies and demonstrates its potential for studying the molecular and cellular mechanisms associated with rejection. However, certain limitations must be borne in mind regarding the development and future applications of this technology.

Spatial Deconvolution of Cell Types and Cell States at Scale Utilizing TACIT

Identifying cell types and states remains a time-consuming, error-prone challenge for spatial biology. While deep learning is increasingly used, it is difficult to generalize due to variability at the level of cells, neighborhoods, and niches in health and disease. To address this, we developed TACIT, an unsupervised algorithm for cell annotation using predefined signatures that operates without training data. TACIT uses unbiased thresholding to distinguish positive cells from background, focusing on relevant markers to identify ambiguous cells in multiomic assays. Using five datasets (5,000,000-cells; 51-cell types) from three niches (brain, intestine, gland), TACIT outperformed existing unsupervised methods in accuracy and scalability. Integrating TACIT-identified cell types with a novel Shiny app revealed new phenotypes in two inflammatory gland diseases. Finally, using combined spatial transcriptomics and proteomics, we discovered under- and overrepresented immune cell types and states in regions of interest, suggesting multimodality is essential for translating spatial biology to clinical applications.

Biogenic amines in the testis: sources, receptors and actions

Biogenic amines are signaling molecules with multiple roles in the central nervous system and in peripheral organs, including the gonads. A series of studies indicated that these molecules, their biosynthetic enzymes and their receptors are present in the testis and that they are involved in the regulation of male reproductive physiology and/or pathology. This mini-review aims to summarize the current knowledge in this field and to pinpoint existing research gaps. We suggest that the widespread clinical use of pharmacological agonists/antagonists of these signaling molecules, calls for new investigations in this area. They are necessary to evaluate the relevance of biogenic amines for human male fertility and infertility, as well as the potential value of at least one of them as an anti-aging compound in the testis.

Human In Vitro Models of Neuroenergetics and Neurometabolic Disturbances: Current Advances and Clinical Perspectives

Neurological conditions conquer the world; they are the leading cause of disability and the second leading cause of death worldwide, and they appear all around the world in every age group, gender, nationality, and socioeconomic class. Despite the growing evidence of an immense impact of perturbations in neuroenergetics on overall brain function, only little is known about the underlying mechanisms. Especially human insights are sparse, owing to a shortage of physiologically relevant model systems. With this perspective, we aim to explore the key steps and considerations involved in developing an advanced human in vitro model for studying neuroenergetics. We discuss biological and technological strategies to meet the requirements of a predictive model, aiming at providing a guide and inspiration for future in vitro models of neuroenergetics.

Advances in Integrated Multi-omics Analysis for Drug-Target Identification

As an essential component of modern drug discovery, the role of drug-target identification is growing increasingly prominent. Additionally, single-omics technologies have been widely utilized in the process of discovering drug targets. However, it is difficult for any single-omics level to clearly expound the causal connection between drugs and how they give rise to the emergence of complex phenotypes. With the progress of large-scale sequencing and the development of high-throughput technologies, the tendency in drug-target identification has shifted towards integrated multi-omics techniques, gradually replacing traditional single-omics techniques. Herein, this review centers on the recent advancements in the domain of integrated multi-omics techniques for target identification, highlights the common multi-omics analysis strategies, briefly summarizes the selection of multi-omics analysis tools, and explores the challenges of existing multi-omics analyses, as well as the applications of multi-omics technology in drug-target identification.

Complete spatially resolved gene expression is not necessary for identifying spatial domains

Spatially resolved transcriptomics (SRT) technologies have revolutionized the study of tissue organization. We introduce a graph convolutional network with an attention and positive emphasis mechanism, termed BINARY, relying exclusively on binarized SRT data to accurately delineate spatial domains. BINARY outperforms existing methods across various SRT data types while using significantly less input information. Our study suggests that precise gene expression quantification may not always be essential, inspiring further exploration of the broader applications of spatially resolved binarized gene expression data.

Sopa: a technology-invariant pipeline for analyses of image-based spatial omics

Spatial omics data allow in-depth analysis of tissue architectures, opening new opportunities for biological discovery. In particular, imaging techniques offer single-cell resolutions, providing essential insights into cellular organizations and dynamics. Yet, the complexity of such data presents analytical challenges and demands substantial computing resources. Moreover, the proliferation of diverse spatial omics technologies, such as Xenium, MERSCOPE, CosMX in spatial-transcriptomics, and MACSima and PhenoCycler in multiplex imaging, hinders the generality of existing tools. We introduce Sopa ( https://github.com/gustaveroussy/sopa ), a technology-invariant, memory-efficient pipeline with a unified visualizer for all image-based spatial omics. Built upon the universal SpatialData framework, Sopa optimizes tasks like segmentation, transcript/channel aggregation, annotation, and geometric/spatial analysis. Its output includes user-friendly web reports and visualizer files, as well as comprehensive data files for in-depth analysis. Overall, Sopa represents a significant step toward unifying spatial data analysis, enabling a more comprehensive understanding of cellular interactions and tissue organization in biological systems.

Advances in the joint profiling technologies of 5mC and 5hmC

DNA cytosine methylation, a crucial epigenetic modification, involves the dynamic interplay of 5-methylcytosine (5mC) and its oxidized form, 5-hydroxymethylcytosine (5hmC), generated by ten-eleven translocation (TET) DNA dioxygenases. This process is central to regulating gene expression, influencing critical biological processes such as development, disease progression, and aging. Recognizing the distinct functions of 5mC and 5hmC, researchers often employ restriction enzyme-based or chemical treatment methods for their simultaneous measurement from the same genomic sample. This enables a detailed understanding of the relationship between these modifications and their collective impact on cellular function. This review focuses on summarizing the technologies for detecting 5mC and 5hmC together but also discusses the limitations and potential future directions in this evolving field.

Spatial Deconvolution of Cell Types and Cell States at Scale Utilizing TACIT

Identifying cell types and states remains a time-consuming and error-prone challenge for spatial biology. While deep learning is increasingly used, it is difficult to generalize due to variability at the level of cells, neighborhoods, and niches in health and disease. To address this, we developed TACIT, an unsupervised algorithm for cell annotation using predefined signatures that operates without training data, using unbiased thresholding to distinguish positive cells from background, focusing on relevant markers to identify ambiguous cells in multiomic assays. Using five datasets (5,000,000-cells; 51-cell types) from three niches (brain, intestine, gland), TACIT outperformed existing unsupervised methods in accuracy and scalability. Integration of TACIT-identified cell with a novel Shiny app revealed new phenotypes in two inflammatory gland diseases. Finally, using combined spatial transcriptomics and proteomics, we discover under- and overrepresented immune cell types and states in regions of interest, suggesting multimodality is essential for translating spatial biology to clinical applications.

Unique challenges and best practices for single cell transcriptomic analysis in toxicology

The application and analysis of single-cell transcriptomics in toxicology presents unique challenges. These include identifying cell sub-populations sensitive to perturbation; interpreting dynamic shifts in cell type proportions in response to chemical exposures; and performing differential expression analysis in dose-response studies spanning multiple treatment conditions. This review examines these challenges while presenting best practices for critical single cell analysis tasks. This covers areas such as cell type identification; analysis of differential cell type abundance; differential gene expression; and cellular trajectories. Towards enhancing the use of single-cell transcriptomics in toxicology, this review aims to address key challenges in this field and offer practical analytical solutions. Overall, applying appropriate bioinformatic techniques to single-cell transcriptomic data can yield valuable insights into cellular responses to toxic exposures.

A practical guide to spatial transcriptomics

Spatial transcriptomics is revolutionizing modern biology, offering researchers an unprecedented ability to unravel intricate gene expression patterns within tissues. From pioneering techniques to newly commercialized platforms, the field of spatial transcriptomics has evolved rapidly, ushering in a new era of understanding across various disciplines, from developmental biology to disease research. This dynamic expansion is reflected in the rapidly growing number of technologies and data analysis techniques developed and introduced. However, the expanding landscape presents a considerable challenge for researchers, especially newcomers to the field, as staying informed about these advancements becomes increasingly complex. To address this challenge, we have prepared an updated review with a particular focus on technologies that have reached commercialization and are, therefore, accessible to a broad spectrum of potential new users. In this review, we present the fundamental principles of spatial transcriptomic methods, discuss the challenges in data analysis, provide insights into experimental considerations, offer information about available resources for spatial transcriptomics, and conclude with a guide for method selection and a forward-looking perspective. Our aim is to serve as a guiding resource for both experienced users and newcomers navigating the complex realm of spatial transcriptomics in this era of rapid development. We intend to equip researchers with the necessary knowledge to make informed decisions and contribute to the cutting-edge research that spatial transcriptomics offers.

CellsFromSpace: a fast, accurate, and reference-free tool to deconvolve and annotate spatially distributed omics data

Motivation

Spatial transcriptomics enables the analysis of cell crosstalk in healthy and diseased organs by capturing the transcriptomic profiles of millions of cells within their spatial contexts. However, spatial transcriptomics approaches also raise new computational challenges for the multidimensional data analysis associated with spatial coordinates.

Results

In this context, we introduce a novel analytical framework called CellsFromSpace based on independent component analysis (ICA), which allows users to analyze various commercially available technologies without relying on a single-cell reference dataset. The ICA approach deployed in CellsFromSpace decomposes spatial transcriptomics data into interpretable components associated with distinct cell types or activities. ICA also enables noise or artifact reduction and subset analysis of cell types of interest through component selection. We demonstrate the flexibility and performance of CellsFromSpace using real-world samples to demonstrate ICA's ability to successfully identify spatially distributed cells as well as rare diffuse cells, and quantitatively deconvolute datasets from the Visium, Slide-seq, MERSCOPE, and CosMX technologies. Comparative analysis with a current alternative reference-free deconvolution tool also highlights CellsFromSpace's speed, scalability and accuracy in processing complex, even multisample datasets. CellsFromSpace also offers a user-friendly graphical interface enabling non-bioinformaticians to annotate and interpret components based on spatial distribution and contributor genes, and perform full downstream analysis.

Availability and implementation

CellsFromSpace (CFS) is distributed as an R package available from github at https://github.com/gustaveroussy/CFS along with tutorials, examples, and detailed documentation.

Spatiotemporal multi-omics: exploring molecular landscapes in aging and regenerative medicine

Aging and regeneration represent complex biological phenomena that have long captivated the scientific community. To fully comprehend these processes, it is essential to investigate molecular dynamics through a lens that encompasses both spatial and temporal dimensions. Conventional omics methodologies, such as genomics and transcriptomics, have been instrumental in identifying critical molecular facets of aging and regeneration. However, these methods are somewhat limited, constrained by their spatial resolution and their lack of capacity to dynamically represent tissue alterations. The advent of emerging spatiotemporal multi-omics approaches, encompassing transcriptomics, proteomics, metabolomics, and epigenomics, furnishes comprehensive insights into these intricate molecular dynamics. These sophisticated techniques facilitate accurate delineation of molecular patterns across an array of cells, tissues, and organs, thereby offering an in-depth understanding of the fundamental mechanisms at play. This review meticulously examines the significance of spatiotemporal multi-omics in the realms of aging and regeneration research. It underscores how these methodologies augment our comprehension of molecular dynamics, cellular interactions, and signaling pathways. Initially, the review delineates the foundational principles underpinning these methods, followed by an evaluation of their recent applications within the field. The review ultimately concludes by addressing the prevailing challenges and projecting future advancements in the field. Indubitably, spatiotemporal multi-omics are instrumental in deciphering the complexities inherent in aging and regeneration, thus charting a course toward potential therapeutic innovations.

Enhancing mass spectrometry imaging accessibility using convolutional autoencoders for deriving hypoxia-associated peptides from tumors

Mass spectrometry imaging (MSI) allows to study cancer's intratumoral heterogeneity through spatially-resolved peptides, metabolites and lipids. Yet, in biomedical research MSI is rarely used for biomarker discovery. Besides its high dimensionality and multicollinearity, mass spectrometry (MS) technologies typically output mass-to-charge ratio values but not the biochemical compounds of interest. Our framework makes particularly low-abundant signals in MSI more accessible. We utilized convolutional autoencoders to aggregate features associated with tumor hypoxia, a parameter with significant spatial heterogeneity, in cancer xenograft models. We highlight that MSI captures these low-abundant signals and that autoencoders can preserve them in their latent space. The relevance of individual hyperparameters is demonstrated through ablation experiments, and the contribution from original features to latent features is unraveled. Complementing MSI with tandem MS from the same tumor model, multiple hypoxia-associated peptide candidates were derived. Compared to random forests alone, our autoencoder approach yielded more biologically relevant insights for biomarker discovery.

Identification of ligand and receptor interactions in CKD and MASH through the integration of single cell and spatial transcriptomics

BACKGROUND

Chronic Kidney Disease (CKD) and Metabolic dysfunction-associated steatohepatitis (MASH) are metabolic fibroinflammatory diseases. Combining single-cell (scRNAseq) and spatial transcriptomics (ST) could give unprecedented molecular disease understanding at single-cell resolution. A more comprehensive analysis of the cell-specific ligand-receptor (L-R) interactions could provide pivotal information about signaling pathways in CKD and MASH. To achieve this, we created an integrative analysis framework in CKD and MASH from two available human cohorts.

RESULTS

The analytical framework identified L-R pairs involved in cellular crosstalk in CKD and MASH. Interactions between cell types identified using scRNAseq data were validated by checking the spatial co-presence using the ST data and the co-expression of the communicating targets. Multiple L-R protein pairs identified are known key players in CKD and MASH, while others are novel potential targets previously observed only in animal models.

CONCLUSION

Our study highlights the importance of integrating different modalities of transcriptomic data for a better understanding of the molecular mechanisms. The combination of single-cell resolution from scRNAseq data, combined with tissue slide investigations and visualization of cell-cell interactions obtained through ST, paves the way for the identification of future potential therapeutic targets and developing effective therapies.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Crosstalk between CD8+ T cells and mesenchymal stromal cells in intestine homeostasis and immunity

The intestine represents the most complex cellular network in the whole body. It is constantly faced with multiple types of immunostimulatory agents encompassing from food antigen, gut microbiome, metabolic waste products, and dead cell debris. Within the intestine, most T cells are found in three primary compartments: the organized gut-associated lymphoid tissue, the lamina propria, and the epithelium. The well-orchestrated epithelial-immune-microbial interaction is critically important for the precise immune response. The main role of intestinal mesenchymal stromal cells is to support a structural framework within the gut wall. However, recent evidence from stromal cell studies indicates that they also possess significant immunomodulatory functions, such as maintaining intestinal tolerance via the expression of PDL1/2 and MHC-II molecules, and promoting the development of CD103+ dendritic cells, and IgA+ plasma cells, thereby enhancing intestinal homeostasis. In this review, we will summarize the current understanding of CD8+ T cells and stromal cells alongside the intestinal tract and discuss the reciprocal interactions between T subsets and mesenchymal stromal cell populations. We will focus on how the tissue residency, migration, and function of CD8+ T cells could be potentially regulated by mesenchymal stromal cell populations and explore the molecular mediators, such as TGF-β, IL-33, and MHC-II molecules that might influence these processes. Finally, we discuss the potential pathophysiological impact of such interaction in intestine hemostasis as well as diseases of inflammation, infection, and malignancies.

Approaches and challenges in identifying, quantifying, and manipulating dynamic mitochondrial genome variations

Mitochondria, the cellular powerhouses, possess their own unique genetic system, including replication, transcription, and translation. Studying these processes is crucial for comprehending mitochondrial disorders, energy production, and their related diseases. Over the past decades, various approaches have been applied in detecting and quantifying mitochondrial genome variations with also the purpose of manipulation of mitochondria or mitochondrial genome for therapeutics. Understanding the scope and limitations of above strategies is not only fundamental to the understanding of basic biology but also critical for exploring disease-related novel target(s), as well to develop innovative therapies. Here, this review provides an overview of different tools and techniques for accurate mitochondrial genome variations identification, quantification, and discuss novel strategies for the manipulation of mitochondria to develop innovative therapeutic interventions, through combining the insights gained from the study of mitochondrial genetics with ongoing single cell omics combined with advanced single molecular tools.

Unveiling aging dynamics in the hematopoietic system insights from single-cell technologies

As the demographic structure shifts towards an aging society, strategies aimed at slowing down or reversing the aging process become increasingly essential. Aging is a major predisposing factor for many chronic diseases in humans. The hematopoietic system, comprising blood cells and their associated bone marrow microenvironment, intricately participates in hematopoiesis, coagulation, immune regulation and other physiological phenomena. The aging process triggers various alterations within the hematopoietic system, serving as a spectrum of risk factors for hematopoietic disorders, including clonal hematopoiesis, immune senescence, myeloproliferative neoplasms and leukemia. The emerging single-cell technologies provide novel insights into age-related changes in the hematopoietic system. In this review, we summarize recent studies dissecting hematopoietic system aging using single-cell technologies. We discuss cellular changes occurring during aging in the hematopoietic system at the levels of the genomics, transcriptomics, epigenomics, proteomics, metabolomics and spatial multi-omics. Finally, we contemplate the future prospects of single-cell technologies, emphasizing the impact they may bring to the field of hematopoietic system aging research.

Spatial genomics: mapping human steatotic liver disease

Metabolic dysfunction-associated steatotic liver disease (MASLD, formerly known as non-alcoholic fatty liver disease) is a leading cause of chronic liver disease worldwide. MASLD can progress to metabolic dysfunction-associated steatohepatitis (MASH, formerly known as non-alcoholic steatohepatitis) with subsequent liver cirrhosis and hepatocellular carcinoma formation. The advent of current technologies such as single-cell and single-nuclei RNA sequencing have transformed our understanding of the liver in homeostasis and disease. The next frontier is contextualizing this single-cell information in its native spatial orientation. This understanding will markedly accelerate discovery science in hepatology, resulting in a further step-change in our knowledge of liver biology and pathobiology. In this Review, we discuss up-to-date knowledge of MASLD development and progression and how the burgeoning field of spatial genomics is driving exciting new developments in our understanding of human liver disease pathogenesis and therapeutic target identification.

Tertiary lymphoid structures in cancer: maturation and induction

Tertiary lymphoid structure (TLS) is an ectopic lymphocyte aggregate formed in peripheral non-lymphoid tissues, including inflamed or cancerous tissue. Tumor-associated TLS serves as a prominent center of antigen presentation and adaptive immune activation within the periphery, which has exhibited positive prognostic value in various cancers. In recent years, the concept of maturity regarding TLS has been proposed and mature TLS, characterized by well-developed germinal centers, exhibits a more potent tumor-suppressive capacity with stronger significance. Meanwhile, more and more evidence showed that TLS can be induced by therapeutic interventions during cancer treatments. Thus, the evaluation of TLS maturity and the therapeutic interventions that induce its formation are critical issues in current TLS research. In this review, we aim to provide a comprehensive summary of the existing classifications for TLS maturity and therapeutic strategies capable of inducing its formation in tumors.

The 3D Revolution in Cancer Discovery

SUMMARY

The transition from 2D to 3D spatial profiling marks a revolutionary era in cancer research, offering unprecedented potential to enhance cancer diagnosis and treatment. This commentary outlines the experimental and computational advancements and challenges in 3D spatial molecular profiling, underscoring the innovation needed in imaging tools, software, artificial intelligence, and machine learning to overcome implementation hurdles and harness the full potential of 3D analysis in the field.

Mapping cutaneous field carcinogenesis of nonmelanoma skin cancer using mesoscopic imaging of pro-inflammation cues

Nonmelanoma skin cancers remain the most widely diagnosed types of cancers globally. Thus, for optimal patient management, it has become imperative that we focus our efforts on the detection and monitoring of cutaneous field carcinogenesis. The concept of field cancerization (or field carcinogenesis), introduced by Slaughter in 1953 in the context of oral cancer, suggests that invasive cancer may emerge from a molecularly and genetically altered field affecting a substantial area of underlying tissue including the skin. A carcinogenic field alteration, present in precancerous tissue over a relatively large area, is not easily detected by routine visualization. Conventional dermoscopy and microscopy imaging are often limited in assessing the entire carcinogenic landscape. Recent efforts have suggested the use of noninvasive mesoscopic (between microscopic and macroscopic) optical imaging methods that can detect chronic inflammatory features to identify pre-cancerous and cancerous angiogenic changes in tissue microenvironments. This concise review covers major types of mesoscopic optical imaging modalities capable of assessing pro-inflammatory cues by quantifying blood haemoglobin parameters and hemodynamics. Importantly, these imaging modalities demonstrate the ability to detect angiogenesis and inflammation associated with actinically damaged skin. Representative experimental preclinical and human clinical studies using these imaging methods provide biological and clinical relevance to cutaneous field carcinogenesis in altered tissue microenvironments in the apparently normal epidermis and dermis. Overall, mesoscopic optical imaging modalities assessing chronic inflammatory hyperemia can enhance the understanding of cutaneous field carcinogenesis, offer a window of intervention and monitoring for actinic keratoses and nonmelanoma skin cancers and maximise currently available treatment options.

Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology

Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.

A revised conceptual framework for mouse vomeronasal pumping and stimulus sampling

The physiological performance of any sensory organ is determined by its anatomy and physical properties. Consequently, complex sensory structures with elaborate features have evolved to optimize stimulus detection. Understanding these structures and their physical nature forms the basis for mechanistic insights into sensory function. Despite its crucial role as a sensor for pheromones and other behaviorally instructive chemical cues, the vomeronasal organ (VNO) remains a poorly characterized mammalian sensory structure. Fundamental principles of its physico-mechanical function, including basic aspects of stimulus sampling, remain poorly explored. Here, we revisit the classical vasomotor pump hypothesis of vomeronasal stimulus uptake. Using advanced anatomical, histological, and physiological methods, we demonstrate that large parts of the lateral mouse VNO are composed of smooth muscle. Vomeronasal smooth muscle tissue comprises two subsets of fibers with distinct topography, structure, excitation-contraction coupling, and, ultimately, contractile properties. Specifically, contractions of a large population of noradrenaline-sensitive cells mediate both transverse and longitudinal lumen expansion, whereas cholinergic stimulation targets an adluminal group of smooth muscle fibers. The latter run parallel to the VNO's rostro-caudal axis and are ideally situated to mediate antagonistic longitudinal constriction of the lumen. This newly discovered arrangement implies a novel mode of function. Single-cell transcriptomics and pharmacological profiling reveal the receptor subtypes involved. Finally, 2D/3D tomography provides non-invasive insight into the intact VNO's anatomy and mechanics, enables measurement of luminal fluid volume, and allows an assessment of relative volume change upon noradrenergic stimulation. Together, we propose a revised conceptual framework for mouse vomeronasal pumping and, thus, stimulus sampling.

Digital pathology and spatial omics in steatohepatitis: Clinical applications and discovery potentials

Steatohepatitis with diverse etiologies is the most common histological manifestation in patients with liver disease. However, there are currently no specific histopathological features pathognomonic for metabolic dysfunction-associated steatotic liver disease, alcohol-associated liver disease, or metabolic dysfunction-associated steatotic liver disease with increased alcohol intake. Digitizing traditional pathology slides has created an emerging field of digital pathology, allowing for easier access, storage, sharing, and analysis of whole-slide images. Artificial intelligence (AI) algorithms have been developed for whole-slide images to enhance the accuracy and speed of the histological interpretation of steatohepatitis and are currently employed in biomarker development. Spatial biology is a novel field that enables investigators to map gene and protein expression within a specific region of interest on liver histological sections, examine disease heterogeneity within tissues, and understand the relationship between molecular changes and distinct tissue morphology. Here, we review the utility of digital pathology (using linear and nonlinear microscopy) augmented with AI analysis to improve the accuracy of histological interpretation. We will also discuss the spatial omics landscape with special emphasis on the strengths and limitations of established spatial transcriptomics and proteomics technologies and their application in steatohepatitis. We then highlight the power of multimodal integration of digital pathology augmented by machine learning (ML)algorithms with spatial biology. The review concludes with a discussion of the current gaps in knowledge, the limitations and premises of these tools and technologies, and the areas of future research.

Spatial analysis toolkits for RNA in situ sequencing

Spatial transcriptomics (ST) is featured by high-throughput gene expression profiling within their native cell and tissue context, offering a means to investigate gene regulatory networks in tissue microenvironment. In situ sequencing (ISS) is an imaging-based ST technology that simultaneously detects hundreds to thousands of genes at subcellular resolution. As a highly reproducible and robust technique, ISS has been widely adapted and undergone a series of technical iterations. As the interest in ISS-based spatial transcriptomic analysis grows, scalable and integrated data analysis workflows are needed to facilitate the applications of ISS in different research fields. This review presents the state-of-the-art bioinformatic toolkits for ISS data analysis, which covers the upstream and downstream analysis workflows, including image analysis, cell segmentation, clustering, functional enrichment, detection of spatially variable genes and cell clusters, spatial cell-cell interactions, and trajectory inference. To assist the community in choosing the right tools for their research, the application of each tool and its compatibility with ISS data are reviewed in detailed. Finally, future perspectives and challenges concerning how to integrate heterogeneous tools into a user-friendly analysis pipeline are discussed. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico.

spatialHeatmap: visualizing spatial bulk and single-cell assays in anatomical images

Visualizing spatial assay data in anatomical images is vital for understanding biological processes in cell, tissue, and organ organizations. Technologies requiring this functionality include traditional one-at-a-time assays, and bulk and single-cell omics experiments, including RNA-seq and proteomics. The spatialHeatmap software provides a series of powerful new methods for these needs, and allows users to work with adequately formatted anatomical images from public collections or custom images. It colors the spatial features (e.g. tissues) annotated in the images according to the measured or predicted abundance levels of biomolecules (e.g. mRNAs) using a color key. This core functionality of the package is called a spatial heatmap plot. Single-cell data can be co-visualized in composite plots that combine spatial heatmaps with embedding plots of high-dimensional data. The resulting spatial context information is essential for gaining insights into the tissue-level organization of single-cell data, or vice versa. Additional core functionalities include the automated identification of biomolecules with spatially selective abundance patterns and clusters of biomolecules sharing similar abundance profiles. To appeal to both non-expert and computational users, spatialHeatmap provides a graphical and a command-line interface, respectively. It is distributed as a free, open-source Bioconductor package (https://bioconductor.org/packages/spatialHeatmap) that users can install on personal computers, shared servers, or cloud systems.

Semi-automated approaches for interrogating spatial heterogeneity of tissue samples

Tissues are spatially orchestrated ecosystems composed of heterogeneous cell populations and non-cellular elements. Tissue components' interactions shape the biological processes that govern homeostasis and disease, thus comprehensive insights into tissues' composition are crucial for understanding their biology. Recently, advancements in the spatial biology field enabled the in-depth analyses of tissue architecture at single-cell resolution, while preserving the structural context. The increasing number of biomarkers analyzed, together with whole tissue imaging, generate datasets approaching several hundreds of gigabytes in size, which are rich sources of valuable knowledge but require investments in infrastructure and resources for extracting quantitative information. The analysis of multiplex whole-tissue images requires extensive training and experience in data analysis. Here, we showcase how a set of open-source tools can allow semi-automated image data extraction to study the spatial composition of tissues with a focus on tumor microenvironment (TME). With the use of Lunaphore COMET platform, we interrogated lung cancer specimens where we examined the expression of 20 biomarkers. Subsequently, the tissue composition was interrogated using an in-house optimized nuclei detection algorithm followed by a newly developed image artifact exclusion approach. Thereafter, the data was processed using several publicly available tools, highlighting the compatibility of COMET-derived data with currently available image analysis frameworks. In summary, we showcased an innovative semi-automated workflow that highlights the ease of adoption of multiplex imaging to explore TME composition at single-cell resolution using a simple slide in, data out approach. Our workflow is easily transferrable to various cohorts of specimens to provide a toolset for spatial cellular dissection of the tissue composition.

High-Resolution Spatial Transcriptomic Atlas of Mouse Soleus Muscle: Unveiling Single Cell and Subcellular Heterogeneity in Health and Denervation

Skeletal muscle is essential for both movement and metabolic processes, characterized by a complex and ordered structure. Despite its importance, a detailed spatial map of gene expression within muscle tissue has been challenging to achieve due to the limitations of existing technologies, which struggle to provide high-resolution views. In this study, we leverage the Seq-Scope technique, an innovative method that allows for the observation of the entire transcriptome at an unprecedented submicron spatial resolution. By applying this technique to the mouse soleus muscle, we analyze and compare the gene expression profiles in both healthy conditions and following denervation, a process that mimics aspects of muscle aging. Our approach reveals detailed characteristics of muscle fibers, other cell types present within the muscle, and specific subcellular structures such as the postsynaptic nuclei at neuromuscular junctions, hybrid muscle fibers, and areas of localized expression of genes responsive to muscle injury, along with their histological context. The findings of this research significantly enhance our understanding of the diversity within the muscle cell transcriptome and its variation in response to denervation, a key factor in the decline of muscle function with age. This breakthrough in spatial transcriptomics not only deepens our knowledge of muscle biology but also sets the stage for the development of new therapeutic strategies aimed at mitigating the effects of aging on muscle health, thereby offering a more comprehensive insight into the mechanisms of muscle maintenance and degeneration in the context of aging and disease.

Combined amplification-based single-molecule fluorescence in situ hybridization with immunofluorescence for simultaneous in situ detection of RNAs and proteins

We present a combined amplification-based single-molecule fluorescence in situ hybridization and immunofluorescence (asmFISH-IF) method for the detection of multiple RNAs and proteins simultaneously in cells and formaldehyde-fixed and paraffin-embedded tissue sections. We showed that performing asmFISH before immunofluorescence gives a better IF signal than the opposite. Our asmFISH-IF method could help study the interplay of RNA and protein, helping to understand their functions.

Multiplex Single-Cell Bioprinting for Engineering of Heterogeneous Tissue Constructs with Subcellular Spatial Resolution

Tissue development, function, and disease are largely driven by the spatial organization of individual cells and their cell-cell interactions. Precision engineered tissues with single-cell spatial resolution, therefore, have tremendous potential for next generation disease models, drug discovery, and regenerative therapeutics. Despite significant advancements in biofabrication approaches to improve feature resolution, strategies to fabricate tissues with the exact same organization of individual cells in their native cellular microenvironment have remained virtually non-existent to date. Here we report a method to spatially pattern single cells with up to eight cell phenotypes and subcellular spatial precision. As proof-of-concept we first demonstrate the ability to systematically assess the influence of cellular microenvironments on cell behavior by controllably altering the spatial arrangement of cell types in bioprinted precision cell-cell interaction arrays. We then demonstrate, for the first time, the ability to produce high-fidelity replicas of a patient's annotated cancer biopsy with subcellular resolution. The ability to replicate native cellular microenvironments marks a significant advancement for precision biofabricated in-vitro models, where heterogenous tissues can be engineered with single-cell spatial precision to advance our understanding of complex biological systems in a controlled and systematic manner.