De novo reconstruction of cell interaction landscapes from single-cell spatial transcriptome data with DeepLinc

Spatially Informed Graph Structure Learning Extracts Insights from Spatial Transcriptomics

Embeddings derived from cell graphs hold significant potential for exploring spatial transcriptomics (ST) datasets. Nevertheless, existing methodologies rely on a graph structure defined by spatial proximity, which inadequately represents the diversity inherent in cell-cell interactions (CCIs). This study introduces STAGUE, an innovative framework that concurrently learns a cell graph structure and a low-dimensional embedding from ST data. STAGUE employs graph structure learning to parameterize and refine a cell graph adjacency matrix, enabling the generation of learnable graph views for effective contrastive learning. The derived embeddings and cell graph improve spatial clustering accuracy and facilitate the discovery of novel CCIs. Experimental benchmarks across 86 real and simulated ST datasets show that STAGUE outperforms 15 comparison methods in clustering performance. Additionally, STAGUE delineates the heterogeneity in human breast cancer tissues, revealing the activation of epithelial-to-mesenchymal transition and PI3K/AKT signaling in specific sub-regions. Furthermore, STAGUE identifies CCIs with greater alignment to established biological knowledge than those ascertained by existing graph autoencoder-based methods. STAGUE also reveals the regulatory genes that participate in these CCIs, including those enriched in neuropeptide signaling and receptor tyrosine kinase signaling pathways, thereby providing insights into the underlying biological processes.

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.

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.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Simulating multiple variability in spatially resolved transcriptomics with scCube

A pressing challenge in spatially resolved transcriptomics (SRT) is to benchmark the computational methods. A widely-used approach involves utilizing simulated data. However, biases exist in terms of the currently available simulated SRT data, which seriously affects the accuracy of method evaluation and validation. Herein, we present scCube ( https://github.com/ZJUFanLab/scCube ), a Python package for independent, reproducible, and technology-diverse simulation of SRT data. scCube not only enables the preservation of spatial expression patterns of genes in reference-based simulations, but also generates simulated data with different spatial variability (covering the spatial pattern type, the resolution, the spot arrangement, the targeted gene type, and the tissue slice dimension, etc.) in reference-free simulations. We comprehensively benchmark scCube with existing single-cell or SRT simulators, and demonstrate the utility of scCube in benchmarking spot deconvolution, gene imputation, and resolution enhancement methods in detail through three applications.

The diversification of methods for studying cell-cell interactions and communication

No cell lives in a vacuum, and the molecular interactions between cells define most phenotypes. Transcriptomics provides rich information to infer cell-cell interactions and communication, thus accelerating the discovery of the roles of cells within their communities. Such research relies heavily on algorithms that infer which cells are interacting and the ligands and receptors involved. Specific pressures on different research niches are driving the evolution of next-generation computational tools, enabling new conceptual opportunities and technological advances. More sophisticated algorithms now account for the heterogeneity and spatial organization of cells, multiple ligand types and intracellular signalling events, and enable the use of larger and more complex datasets, including single-cell and spatial transcriptomics. Similarly, new high-throughput experimental methods are increasing the number and resolution of interactions that can be analysed simultaneously. Here, we explore recent progress in cell-cell interaction research and highlight the diversification of the next generation of tools, which have yielded a rich ecosystem of tools for different applications and are enabling invaluable discoveries.

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.

TENET: Triple-enhancement based graph neural network for cell-cell interaction network reconstruction from spatial transcriptomics

Cellular communication relies on the intricate interplay of signaling molecules, forming the Cell-cell Interaction network (CCI) that coordinates tissue behavior. Researchers have shown the capability of shallow neural networks in reconstructing CCI, given molecules' abundance in the Spatial Transcriptomics (ST) data. When encountering situations such as sparse connections in CCI and excessive noise, the susceptibility of shallow networks to these factors significantly impacts the accuracy of CCI reconstruction, resulting in subpar results. To reconstruct a more comprehensive and accurate CCI, we propose a novel method named Triple-Enhancement based Graph Neural Network (TENET). In TENET, three progressive enhancement mechanisms build upon each other, creating a cumulative effect. This approach can ensure the ability to capture valuable features in limited data and amplify the noise signal to facilitate the denoising effect. Additionally, the whole architecture guides the decoding reconstruction phase with integrated knowledge, which leverages the accumulated insights from each stage of enhancement to ensure a refined and comprehensive CCI reconstruction. The presented TENET has been implemented and tested on both real and synthetic ST datasets. Averagely, the CCI reconstruction using TENET achieves a 9.61% improvement in Average Precision (AP) and a 7.32% improvement in Area Under the Receiver Operating Characteristic (AUROC) compared to the existing state-of-the-art (SOTA) method. The source code and data are available at https://github.com/Yujian-Lee/TENET.

Single-Cell Informatics for Tumor Microenvironment and Immunotherapy

Cancer comprises malignant cells surrounded by the tumor microenvironment (TME), a dynamic ecosystem composed of heterogeneous cell populations that exert unique influences on tumor development. The immune community within the TME plays a substantial role in tumorigenesis and tumor evolution. The innate and adaptive immune cells "talk" to the tumor through ligand-receptor interactions and signaling molecules, forming a complex communication network to influence the cellular and molecular basis of cancer. Such intricate intratumoral immune composition and interactions foster the application of immunotherapies, which empower the immune system against cancer to elicit durable long-term responses in cancer patients. Single-cell technologies have allowed for the dissection and characterization of the TME to an unprecedented level, while recent advancements in bioinformatics tools have expanded the horizon and depth of high-dimensional single-cell data analysis. This review will unravel the intertwined networks between malignancy and immunity, explore the utilization of computational tools for a deeper understanding of tumor-immune communications, and discuss the application of these approaches to aid in diagnosis or treatment decision making in the clinical setting, as well as the current challenges faced by the researchers with their potential future improvements.

Revealing cell-cell communication pathways with their spatially coupled gene programs

Inference of cell-cell communication (CCC) provides valuable information in understanding the mechanisms of many important life processes. With the rise of spatial transcriptomics in recent years, many methods have emerged to predict CCCs using spatial information of cells. However, most existing methods only describe CCCs based on ligand-receptor interactions, but lack the exploration of their upstream/downstream pathways. In this paper, we proposed a new method to infer CCCs, called Intercellular Gene Association Network (IGAN). Specifically, it is for the first time that we can estimate the gene associations/network between two specific single spatially adjacent cells. By using the IGAN method, we can not only infer CCCs in an accurate manner, but also explore the upstream/downstream pathways of ligands/receptors from the network perspective, which are actually exhibited as a new panoramic cell-interaction-pathway graph, and thus provide extensive information for the regulatory mechanisms behind CCCs. In addition, IGAN can measure the CCC activity at single cell/spot resolution, and help to discover the CCC spatial heterogeneity. Interestingly, we found that CCC patterns from IGAN are highly consistent with the spatial microenvironment patterns for each cell type, which further indicated the accuracy of our method. Analyses on several public datasets validated the advantages of IGAN.

Navigating the landscapes of spatial transcriptomics: How computational methods guide the way

Spatially resolved transcriptomics has been dramatically transforming biological and medical research in various fields. It enables transcriptome profiling at single-cell, multi-cellular, or sub-cellular resolution, while retaining the information of geometric localizations of cells in complex tissues. The coupling of cell spatial information and its molecular characteristics generates a novel multi-modal high-throughput data source, which poses new challenges for the development of analytical methods for data-mining. Spatial transcriptomic data are often highly complex, noisy, and biased, presenting a series of difficulties, many unresolved, for data analysis and generation of biological insights. In addition, to keep pace with the ever-evolving spatial transcriptomic experimental technologies, the existing analytical theories and tools need to be updated and reformed accordingly. In this review, we provide an overview and discussion of the current computational approaches for mining of spatial transcriptomics data. Future directions and perspectives of methodology design are proposed to stimulate further discussions and advances in new analytical models and algorithms. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Evolution and Genomics > Computational Analyses of RNA RNA Export and Localization > RNA Localization.

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.

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.

The promising application of cell-cell interaction analysis in cancer from single-cell and spatial transcriptomics

Cell-cell interactions instruct cell fate and function. These interactions are hijacked to promote cancer development. Single-cell transcriptomics and spatial transcriptomics have become powerful new tools for researchers to profile the transcriptional landscape of cancer at unparalleled genetic depth. In this review, we discuss the rapidly growing array of computational tools to infer cell-cell interactions from non-spatial single-cell RNA-sequencing and the limited but growing number of methods for spatial transcriptomics data. Downstream analyses of these computational tools and applications to cancer studies are highlighted. We finish by suggesting several directions for further extensions that anticipate the increasing availability of multi-omics cancer data.

Adult tissue-specific stem cell interaction: novel technologies and research advances

Adult tissue-specific stem cells play a dominant role in tissue homeostasis and regeneration. Various in vivo markers of adult tissue-specific stem cells have been increasingly reported by lineage tracing in genetic mouse models, indicating that marked cells differentiation is crucial during homeostasis and regeneration. How adult tissue-specific stem cells with indicated markers contact the adjacent lineage with indicated markers is of significance to be studied. Novel methods bring future findings. Recent advances in lineage tracing, synthetic receptor systems, proximity labeling, and transcriptomics have enabled easier and more accurate cell behavior visualization and qualitative and quantitative analysis of cell-cell interactions than ever before. These technological innovations have prompted researchers to re-evaluate previous experimental results, providing increasingly compelling experimental results for understanding the mechanisms of cell-cell interactions. This review aimed to describe the recent methodological advances of dual enzyme lineage tracing system, the synthetic receptor system, proximity labeling, single-cell RNA sequencing and spatial transcriptomics in the study of adult tissue-specific stem cells interactions. An enhanced understanding of the mechanisms of adult tissue-specific stem cells interaction is important for tissue regeneration and maintenance of homeostasis in organisms.

Single-cell analysis of human prepuce reveals dynamic changes in gene regulation and cellular communications

BACKGROUND

The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells.

RESULTS

We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults.

CONCLUSIONS

Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

CLARIFY: cell-cell interaction and gene regulatory network refinement from spatially resolved transcriptomics

MOTIVATION

Gene regulatory networks (GRNs) in a cell provide the tight feedback needed to synchronize cell actions. However, genes in a cell also take input from, and provide signals to other neighboring cells. These cell-cell interactions (CCIs) and the GRNs deeply influence each other. Many computational methods have been developed for GRN inference in cells. More recently, methods were proposed to infer CCIs using single cell gene expression data with or without cell spatial location information. However, in reality, the two processes do not exist in isolation and are subject to spatial constraints. Despite this rationale, no methods currently exist to infer GRNs and CCIs using the same model.

RESULTS

We propose CLARIFY, a tool that takes GRNs as input, uses them and spatially resolved gene expression data to infer CCIs, while simultaneously outputting refined cell-specific GRNs. CLARIFY uses a novel multi-level graph autoencoder, which mimics cellular networks at a higher level and cell-specific GRNs at a deeper level. We applied CLARIFY to two real spatial transcriptomic datasets, one using seqFISH and the other using MERFISH, and also tested on simulated datasets from scMultiSim. We compared the quality of predicted GRNs and CCIs with state-of-the-art baseline methods that inferred either only GRNs or only CCIs. The results show that CLARIFY consistently outperforms the baseline in terms of commonly used evaluation metrics. Our results point to the importance of co-inference of CCIs and GRNs and to the use of layered graph neural networks as an inference tool for biological networks.

AVAILABILITY AND IMPLEMENTATION

The source code and data is available at https://github.com/MihirBafna/CLARIFY.

Deciphering plant cell-cell communications using single-cell omics data

Plants have various cell types that respond to different environmental factors, and cell-cell communication is the fundamental process that controls these plant responses. The emergence of single-cell techniques provides opportunities to explore features unique to each cell type and construct a comprehensive cell-cell communication (CCC) network. Although the most current successes of CCC inference were achieved in animal research, computational methods can also be directly applied to plants. This review describes the current major models for cell-cell communication inference and summarizes the computational tools based on single-cell omics datasets. In addition, we discuss the limitations of plant cell-cell communication research and propose new directions to expand the field in meaningful ways.

Machine Learning for Uncovering Biological Insights in Spatial Transcriptomics Data

Development and homeostasis in multicellular systems both require exquisite control over spatial molecular pattern formation and maintenance. Advances in spatially-resolved and high-throughput molecular imaging methods such as multiplexed immunofluorescence and spatial transcriptomics (ST) provide exciting new opportunities to augment our fundamental understanding of these processes in health and disease. The large and complex datasets resulting from these techniques, particularly ST, have led to rapid development of innovative machine learning (ML) tools primarily based on deep learning techniques. These ML tools are now increasingly featured in integrated experimental and computational workflows to disentangle signals from noise in complex biological systems. However, it can be difficult to understand and balance the different implicit assumptions and methodologies of a rapidly expanding toolbox of analytical tools in ST. To address this, we summarize major ST analysis goals that ML can help address and current analysis trends. We also describe four major data science concepts and related heuristics that can help guide practitioners in their choices of the right tools for the right biological questions.

Screening cell-cell communication in spatial transcriptomics via collective optimal transport

Spatial transcriptomic technologies and spatially annotated single-cell RNA sequencing datasets provide unprecedented opportunities to dissect cell-cell communication (CCC). However, incorporation of the spatial information and complex biochemical processes required in the reconstruction of CCC remains a major challenge. Here, we present COMMOT (COMMunication analysis by Optimal Transport) to infer CCC in spatial transcriptomics, which accounts for the competition between different ligand and receptor species as well as spatial distances between cells. A collective optimal transport method is developed to handle complex molecular interactions and spatial constraints. Furthermore, we introduce downstream analysis tools to infer spatial signaling directionality and genes regulated by signaling using machine learning models. We apply COMMOT to simulation data and eight spatial datasets acquired with five different technologies to show its effectiveness and robustness in identifying spatial CCC in data with varying spatial resolutions and gene coverages. Finally, COMMOT identifies new CCCs during skin morphogenesis in a case study of human epidermal development.