CytoTalk: De novo construction of signal transduction networks using single-cell transcriptomic data

CausalCCC: a web server to explore intracellular causal pathways enabling cell-cell communication

Understanding cell-cell communication (CCC) pathways from single-cell or spatial transcriptomic data is key to unraveling biological processes. Recently, multiple CCC methods have been developed but primarily focus on refining ligand-receptor (L-R) interaction scores. A critical gap for a more comprehensive picture of cellular crosstalks lies in the integration of upstream and downstream intracellular pathways in the sender and receiver cells. We report here CausalCCC, https://miic.curie.fr/causalCCC.php, an interactive web server, which addresses this need by reconstructing gene-gene interaction pathways across two or more interacting cell types from single-cell or spatial transcriptomic data. CausalCCC includes a graphical introduction and a demo dataset within the workbench page as well as a comprehensive tutorial. CausalCCC methodology integrates a robust and scalable causal network reconstruction method, multivariate information-based inductive causation, with internally computed L-R pairs using LIANA+ (including CellphoneDBv5, SingleCellSignalR, Connectome, NATMI, and Log2FC). Alternatively, user-defined L-R pairs from any CCC methods can also be uploaded. We showcase here CausalCCC on different single-cell and spatial transcriptomic datasets from three original CCC methods (NicheNet, CellChat, and Misty). CausalCCC web server offers unique interactive visualization tools dedicated to single-cell data practitioners seeking to go beyond L-R scores and explore extended CCC pathways across multiple interacting cell types.

scMINER: a mutual information-based framework for clustering and hidden driver inference from single-cell transcriptomics data

Single-cell transcriptomics data present challenges due to their inherent stochasticity and sparsity, complicating both cell clustering and cell type-specific network inference. To address these challenges, we introduce scMINER (single-cell Mutual Information-based Network Engineering Ranger), an integrative framework for unsupervised cell clustering, transcription factor and signaling protein network inference, and identification of hidden drivers from single-cell transcriptomic data. scMINER demonstrates superior accuracy in cell clustering, outperforming five state-of-the-art algorithms and excelling in distinguishing closely related cell populations. For network inference, scMINER outperforms three established methods, as validated by ATAC-seq and CROP-seq. In particular, it surpasses SCENIC in revealing key transcription factor drivers involved in T cell exhaustion and Treg tissue specification. Moreover, scMINER enables the inference of signaling protein networks and drivers with high accuracy, which presents an advantage in multimodal single cell data analysis. In addition, we establish scMINER Portal, an interactive visualization tool to facilitate exploration of scMINER results.

Longitudinal single-cell multiomic atlas of high-risk neuroblastoma reveals chemotherapy-induced tumor microenvironment rewiring

High-risk neuroblastoma, a leading cause of pediatric cancer mortality, exhibits substantial intratumoral heterogeneity, contributing to therapeutic resistance. To understand tumor microenvironment evolution during therapy, we longitudinally profiled 22 patients with high-risk neuroblastoma before and after induction chemotherapy using single-nucleus RNA and ATAC sequencing and whole-genome sequencing. This revealed profound shifts in tumor and immune cell subpopulations after therapy and identified enhancer-driven transcriptional regulators of neuroblastoma neoplastic states. Poor outcome correlated with proliferative and metabolically active neoplastic states, whereas more differentiated neuronal-like states predicted better prognosis. Proportions of mesenchymal neoplastic cells increased after therapy and a high proportion correlated with a poorer chemotherapy response. Macrophages significantly expanded towards pro-angiogenic, immunosuppressive and metabolic phenotypes. We identified paracrine signaling networks and validated the HB-EGF-ERBB4 axis between macrophage and neoplastic subsets, which promoted tumor growth through the induction of ERK signaling. These findings collectively reveal intrinsic and extrinsic regulators of therapy response in high-risk neuroblastoma.

New Insights and Implications of Cell-Cell Interactions in Developmental Biology

The dynamic and meticulously regulated networks established the foundation for embryonic development, where the intercellular interactions and signal transduction assumed a pivotal role. In recent years, high-throughput technologies such as single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) have advanced dramatically, empowering the systematic dissection of cell-to-cell regulatory networks. The emergence of comprehensive databases and analytical frameworks has further provided unprecedented insights into embryonic development and cell-cell interactions (CCIs). This paper reviewed the exponential increased CCIs works related to developmental biology from 2008 to 2023, comprehensively collected and categorized 93 analytical tools and 39 databases, and demonstrated its practical utility through illustrative case studies. In parallel, the article critically scrutinized the persistent challenges within this field, such as the intricacies of spatial localization and transmembrane state validation at single-cell resolution, and underscored the interpretative limitations inherent in current analytical frameworks. The development of CCIs' analysis tools with harmonizing multi-omics data and the construction of cross-species dynamically updated CCIs databases will be the main direction of future research. Future investigations into CCIs are poised to expeditiously drive the application and clinical translation within developmental biology, unlocking novel dimensions for exploration and progress.

Extracellular vesicle-derived miRNA-mediated cell-cell communication inference for single-cell transcriptomic data with miRTalk

MicroRNAs are released from cells in extracellular vesicles (EVs), representing an essential mode of cell-cell communication (CCC) via a regulatory effect on gene expression. Single-cell RNA-sequencing technologies have ushered in an era of elucidating CCC at single-cell resolution. Herein, we present miRTalk, a pioneering approach for inferring CCC mediated by EV-derived miRNA-target interactions (MiTIs). The benchmarking against simulated and real-world datasets demonstrates the superior performance of miRTalk, and the application to four disease scenarios reveals the in-depth MiTI-mediated CCC mechanisms. Collectively, miRTalk can infer EV-derived MiTI-mediated CCC with scRNA-seq data, providing new insights into the intercellular dynamics of biological processes.

CellPhoneDB v5: inferring cell-cell communication from single-cell multiomics data

Cell-cell communication is essential for tissue development, function and regeneration. The revolution of single-cell genomics technologies offers an unprecedented opportunity to uncover how cells communicate in vivo within their tissue niches and how disruption of these niches can lead to diseases and developmental abnormalities. CellPhoneDB is a bioinformatics toolkit designed to infer cell-cell communication by combining a curated repository of bona fide ligand-receptor interactions with methods to integrate these interactions with single-cell genomics data. Here we present a protocol for the latest version of CellPhoneDB (v5), offering several new features. First, the repository has been expanded by one-third with the addition of new interactions, including ~1,000 interactions mediated by nonpeptidic ligands such as steroidogenic hormones, neurotransmitters and small G-protein-coupled receptor (GPCR)-binding ligands. Second, we outline a new way of using the database that allows users to tailor queries to their experimental designs. Third, the update incorporates novel strategies to prioritize specific cell-cell interactions, leveraging information from other modalities such as tissue microenvironments derived from spatial transcriptomics technologies or transcription factor activities derived from a single-cell assay for transposase accessible chromatin assays. Finally, we describe the new CellPhoneDBViz module to interactively visualize and share results. Altogether, CellPhoneDB v5 enhances the precision of cell-cell communication inference, offering new insights into tissue biology in physiological microenvironments. This protocol typically takes ~15 min and requires basic knowledge of python.

Graph-Based 3-Dimensional Spatial Gene Neighborhood Networks of Single Cells in Gels and Tissues

Objective: We developed 3-dimensional spatially resolved gene neighborhood network embedding (3D-spaGNN-E) to find subcellular gene proximity relationships and identify key subcellular motifs in cell-cell communication (CCC). Impact Statement: The pipeline combines 3D imaging-based spatial transcriptomics and graph-based deep learning to identify subcellular motifs. Introduction: Advancements in imaging and experimental technology allow the study of 3D spatially resolved transcriptomics and capture better spatial context than approximating the samples as 2D. However, the third spatial dimension increases the data complexity and requires new analyses. Methods: 3D-spaGNN-E detects single transcripts in 3D cell culture samples and identifies subcellular gene proximity relationships. Then, a graph autoencoder projects the gene proximity relationships into a latent space. We then applied explainability analysis to identify subcellular CCC motifs. Results: We first applied the pipeline to mesenchymal stem cells (MSCs) cultured in hydrogel. After clustering the cells based on the RNA count, we identified cells belonging to the same cluster as homotypic and those belonging to different clusters as heterotypic. We identified changes in local gene proximity near the border between homotypic and heterotypic cells. When applying the pipeline to the MSC-peripheral blood mononuclear cell (PBMC) coculture system, we identified CD4+ and CD8+ T cells. Local gene proximity and autoencoder embedding changes can distinguish strong and weak suppression of different immune cells. Lastly, we compared astrocyte-neuron CCC in mouse hypothalamus and cortex by analyzing 3D multiplexed-error-robust fluorescence in situ hybridization (MERFISH) data and identified regional gene proximity differences. Conclusion: 3D-spaGNN-E distinguished distinct CCCs in cell culture and tissue by examining subcellular motifs.

Deep learning in single-cell and spatial transcriptomics data analysis: advances and challenges from a data science perspective

The development of single-cell and spatial transcriptomics has revolutionized our capacity to investigate cellular properties, functions, and interactions in both cellular and spatial contexts. Despite this progress, the analysis of single-cell and spatial omics data remains challenging. First, single-cell sequencing data are high-dimensional and sparse, and are often contaminated by noise and uncertainty, obscuring the underlying biological signal. Second, these data often encompass multiple modalities, including gene expression, epigenetic modifications, metabolite levels, and spatial locations. Integrating these diverse data modalities is crucial for enhancing prediction accuracy and biological interpretability. Third, while the scale of single-cell sequencing has expanded to millions of cells, high-quality annotated datasets are still limited. Fourth, the complex correlations of biological tissues make it difficult to accurately reconstruct cellular states and spatial contexts. Traditional feature engineering approaches struggle with the complexity of biological networks, while deep learning, with its ability to handle high-dimensional data and automatically identify meaningful patterns, has shown great promise in overcoming these challenges. Besides systematically reviewing the strengths and weaknesses of advanced deep learning methods, we have curated 21 datasets from nine benchmarks to evaluate the performance of 58 computational methods. Our analysis reveals that model performance can vary significantly across different benchmark datasets and evaluation metrics, providing a useful perspective for selecting the most appropriate approach based on a specific application scenario. We highlight three key areas for future development, offering valuable insights into how deep learning can be effectively applied to transcriptomic data analysis in biological, medical, and clinical settings.

Learning directed acyclic graphs for ligands and receptors based on spatially resolved transcriptomic data of ovarian cancer

To unravel the mechanism of immune activation and suppression within tumors, a critical step is to identify transcriptional signals governing cell-cell communication between tumor and immune/stromal cells in the tumor microenvironment. Central to this communication are interactions between secreted ligands and cell-surface receptors, creating a highly connected signaling network among cells. Recent advancements in in situ-omics profiling, particularly spatial transcriptomic (ST) technology, provide unique opportunities to directly characterize ligand-receptor signaling networks that power cell-cell communication. In this paper, we propose a novel statistical method, LRnetST, to characterize the ligand-receptor interaction networks between adjacent tumor and immune/stroma cells based on ST data. LRnetST utilizes a directed acyclic graph model with a novel approach to handle the zero-inflated distributions of ST data. It also leverages existing ligand-receptor regulation databases as prior information, and employs a bootstrap aggregation strategy to achieve robust network estimation. Application of LRnetST to ST data of high-grade serous ovarian tumor samples revealed both common and distinct ligand-receptor regulations across different tumors. Some of these interactions were validated through both a MERFISH dataset and a CosMx SMI dataset of independent ovarian tumor samples. These results cast light on biological processes relating to the communication between tumor and immune/stromal cells in ovarian tumors. An open-source R package of LRnetST is available on GitHub at https://github.com/jie108/LRnetST.

Unraveling cell-cell communication with NicheNet by inferring active ligands from transcriptomics data

Ligand-receptor interactions constitute a fundamental mechanism of cell-cell communication and signaling. NicheNet is a well-established computational tool that infers ligand-receptor interactions that potentially regulate gene expression changes in receiver cell populations. Whereas the original publication delves into the algorithm and validation, this paper describes a best practices workflow cultivated over four years of experience and user feedback. Starting from the input single-cell expression matrix, we describe a 'sender-agnostic' approach that considers ligands from the entire microenvironment and a 'sender-focused' approach that considers ligands only from cell populations of interest. As output, users will obtain a list of prioritized ligands and their potential target genes, along with multiple visualizations. We include further developments made in NicheNet v2, in which we have updated the data sources and implemented a downstream procedure for prioritizing cell type-specific ligand-receptor pairs. Although a standard NicheNet analysis takes <10 min to run, users often invest additional time in making decisions about the approach and parameters that best suit their biological question. This paper serves to aid in this decision-making process by describing the most appropriate workflow for common experimental designs like case-control and cell-differentiation studies. Finally, in addition to the step-by-step description of the code, we also provide wrapper functions that enable the analysis to be run in one line of code, thus tailoring the workflow to users at all levels of computational proficiency.

SpaCcLink: exploring downstream signaling regulations with graph attention network for systematic inference of spatial cell-cell communication

BACKGROUND

Cellular communication is vital for the proper functioning of multicellular organisms. A comprehensive analysis of cellular communication demands the consideration not only of the binding between ligands and receptors but also of a series of downstream signal transduction reactions within cells. Thanks to the advancements in spatial transcriptomics technology, we are now able to better decipher the process of cellular communication within the cellular microenvironment. Nevertheless, the majority of existing spatial cell-cell communication algorithms fail to take into account the downstream signals within cells.

RESULTS

In this study, we put forward SpaCcLink, a cell-cell communication analysis method that takes into account the downstream influence of individual receptors within cells and systematically investigates the spatial patterns of communication as well as downstream signal networks. Analyses conducted on real datasets derived from humans and mice have demonstrated that SpaCcLink can help in identifying more relevant ligands and receptors, thereby enabling us to systematically decode the downstream genes and signaling pathways that are influenced by cell-cell communication. Comparisons with other methods suggest that SpaCcLink can identify downstream genes that are more closely associated with biological processes and can also discover reliable ligand-receptor relationships.

CONCLUSIONS

By means of SpaCcLink, a more profound and all-encompassing comprehension of the mechanisms underlying cellular communication can be achieved, which in turn promotes and deepens our understanding of the intricate complexity within organisms.

Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data

Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.

Causal genes identification of giant cell arteritis in CD4+ Memory t cells: an integration of multi-omics and expression quantitative trait locus analysis

BACKGROUND

Giant cell arteritis (GCA) is a prevalent artery and is strongly correlated with age. The role of CD4+ Memory T cells in giant cell arteritis has not been elucidated.

METHOD

Through single-cell analysis, we focused on the CD4+ Memory T cells in giant cell arteritis. eQTL analysis and mendelian randomization analysis identified the significant genes which have a causal effect on giant cell arteritis risk. CD4+ Memory T cells were subsequently divided into gene-positive and gene-negative groups, then further single-cell analysis was conducted. Mendelian randomization of plasma proteins, blood-urine biomarkers and metabolites were also performed. Eventually, the PMA induced Jurkat cell lines were used for biological experiments to explore the specific functions of significant causal genes in CD4+ Memory T cells.

RESULTS

Similarity of CD4+ Memory T cells in GCA and old samples were explored. DDIT4 and ARHGAP15 were identified as significant risk genes via mendelian randomization. The CD4+ Memory T cells were then divided into DDIT4 ± or ARHGAP15 ± groups, and further single-cell analysis indicated the differences in aspects involving intercellular communication, functional pathways, protein activity, metabolism and drug sensitivity between positive and negative groups. In vitro experiments, including overexpression and knockdown, demonstrated that DDIT4 leading to a chronic, low-intensity inflammatory state in CD4+ Memory T cells, eventually promoting the development of GCA.

CONCLUSION

DDIT4 and ARHGAP15 have significant causal effects on giant cell arteritis risk. Specifically, DDIT4 exhibit pro-inflammatory effects on GCA via promotes chronic, low-intensity inflammatory in CD4+ Memory T cell.

CellChat for systematic analysis of cell-cell communication from single-cell transcriptomics

Recent advances in single-cell sequencing technologies offer an opportunity to explore cell-cell communication in tissues systematically and with reduced bias. A key challenge is integrating known molecular interactions and measurements into a framework to identify and analyze complex cell-cell communication networks. Previously, we developed a computational tool, named CellChat, that infers and analyzes cell-cell communication networks from single-cell transcriptomic data within an easily interpretable framework. CellChat quantifies the signaling communication probability between two cell groups using a simplified mass-action-based model, which incorporates the core interaction between ligands and receptors with multisubunit structure along with modulation by cofactors. Importantly, CellChat performs a systematic and comparative analysis of cell-cell communication using a variety of quantitative metrics and machine-learning approaches. CellChat v2 is an updated version that includes additional comparison functionalities, an expanded database of ligand-receptor pairs along with rich functional annotations, and an Interactive CellChat Explorer. Here we provide a step-by-step protocol for using CellChat v2 on single-cell transcriptomic data, including inference and analysis of cell-cell communication from one dataset and identification of altered intercellular communication, signals and cell populations from different datasets across biological conditions. The R implementation of CellChat v2 toolkit and its tutorials together with the graphic outputs are available at https://github.com/jinworks/CellChat . This protocol typically takes ~5 min depending on dataset size and requires a basic understanding of R and single-cell data analysis but no specialized bioinformatics training for its implementation.

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.

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

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

Detecting global and local hierarchical structures in cell-cell communication using CrossChat

Cell-cell communication (CCC) occurs across different biological scales, ranging from interactions between large groups of cells to interactions between individual cells, forming a hierarchical structure. Globally, CCC may exist between clusters or only subgroups of a cluster with varying size, while locally, a group of cells as sender or receiver may exhibit distinct signaling properties. Current existing methods infer CCC from single-cell RNA-seq or Spatial Transcriptomics only between predefined cell groups, neglecting the existing hierarchical structure within CCC that are determined by signaling molecules, in particular, ligands and receptors. Here, we develop CrossChat, a novel computational framework designed to infer and analyze the hierarchical cell-cell communication structures using two complementary approaches: a global hierarchical structure using a multi-resolution clustering method, and multiple local hierarchical structures using a tree detection method. This framework provides a comprehensive approach to understand the hierarchical relationships within CCC that govern complex tissue functions. By applying our method to two nonspatial scRNA-seq datasets sampled from COVID-19 patients and mouse embryonic skin, and two spatial transcriptomics datasets generated from Stereo-seq of mouse embryo and 10x Visium of mouse wounded skin, we showcase CrossChat's functionalities for analyzing both global and local hierarchical structures within cell-cell communication.

Intelligent sensing for the autonomous manipulation of microrobots toward minimally invasive cell surgery

The physiology and pathogenesis of biological cells have drawn enormous research interest. Benefiting from the rapid development of microfabrication and microelectronics, miniaturized robots with a tool size below micrometers have widely been studied for manipulating biological cells in vitro and in vivo. Traditionally, the complex physiological environment and biological fragility require human labor interference to fulfill these tasks, resulting in high risks of irreversible structural or functional damage and even clinical risk. Intelligent sensing devices and approaches have been recently integrated within robotic systems for environment visualization and interaction force control. As a consequence, microrobots can be autonomously manipulated with visual and interaction force feedback, greatly improving accuracy, efficiency, and damage regulation for minimally invasive cell surgery. This review first explores advanced tactile sensing in the aspects of sensing principles, design methodologies, and underlying physics. It also comprehensively discusses recent progress on visual sensing, where the imaging instruments and processing methods are summarized and analyzed. It then introduces autonomous micromanipulation practices utilizing visual and tactile sensing feedback and their corresponding applications in minimally invasive surgery. Finally, this work highlights and discusses the remaining challenges of current robotic micromanipulation and their future directions in clinical trials, providing valuable references about this field.

A global view of altered ligand-receptor interactions in bone marrow aging based on single-cell sequencing

Altered cell-cell communication is a hallmark of aging, but its impact on bone marrow aging remains poorly understood. Based on a common and effective pipeline and single-cell transcriptome sequencing, we detected 384,124 interactions including 2575 ligand-receptor pairs and 16 non-adherent bone marrow cell types in old and young mouse and identified a total of 5560 significantly different interactions, which were then verified by flow cytometry and quantitative real-time PCR. These differential ligand-receptor interactions exhibited enrichment for the senescence-associated secretory phenotypes. Further validation demonstrated supplementing specific extracellular ligands could modify the senescent signs of hematopoietic stem cells derived from old mouse. Our work provides an effective procedure to detect the ligand-receptor interactions based on single-cell sequencing, which contributes to understand mechanisms and provides a potential strategy for intervention of bone marrow aging.

Crosstalk of T cells within the ovarian cancer microenvironment

Ovarian cancer (OC) represents ecosystems of highly diverse tumor microenvironments (TMEs). The presence of tumor-infiltrating lymphocytes (TILs) is linked to enhanced immune responses and long-term survival. In this review we present emerging evidence suggesting that cellular crosstalk tightly regulates the distribution of TILs within the TME, underscoring the need to better understand key cellular networks that promote or impede T cell infiltration in OC. We also capture the emergent methodologies and computational techniques that enable the dissection of cell-cell crosstalk. Finally, we present innovative ex vivo TME models that can be leveraged to map and perturb cellular communications to enhance T cell infiltration and immune reactivity.

THGB: predicting ligand-receptor interactions by combining tree boosting and histogram-based gradient boosting

Ligand-receptor interaction (LRI) prediction has great significance in biological and medical research and facilitates to infer and analyze cell-to-cell communication. However, wet experiments for new LRI discovery are costly and time-consuming. Here, we propose a computational model called THGB to uncover new LRIs. THGB first extracts feature information of Ligand-Receptor (LR) pairs using iFeature. Next, it adopts a tree boosting model to obtain representative LR features. Finally, it devises the histogram-based gradient boosting model to capture high-quality LRIs. To assess the THGB performance, we compared it with three new LRI prediction models (i.e., CellEnBoost, CellGiQ, and CellComNet) and one classical protein-protein interaction inference model PIPR. The results demonstrated that THGB achieved the best overall predictions in terms of six evaluation indictors (i.e., precision, recall, accuracy, F1-score, AUC, and AUPR). To measure the effect of LR feature selection on the prediction, THGB was compared with four feature selection methods (i.e., PCA, NMF, LLE, and TSVD). The results showed that the tree boosting model was more appropriate to select representative LR features and improve LRI prediction. We also conducted ablation study and found that THGB with feature selection outperformed THGB without feature selection. We hope that THGB is a useful tool to find new LRIs and further infer cell-to-cell communication.

VGAE-CCI: variational graph autoencoder-based construction of 3D spatial cell-cell communication network

Cell-cell communication plays a critical role in maintaining normal biological functions, regulating development and differentiation, and controlling immune responses. The rapid development of single-cell RNA sequencing and spatial transcriptomics sequencing (ST-seq) technologies provides essential data support for in-depth and comprehensive analysis of cell-cell communication. However, ST-seq data often contain incomplete data and systematic biases, which may reduce the accuracy and reliability of predicting cell-cell communication. Furthermore, other methods for analyzing cell-cell communication mainly focus on individual tissue sections, neglecting cell-cell communication across multiple tissue layers, and fail to comprehensively elucidate cell-cell communication networks within three-dimensional tissues. To address the aforementioned issues, we propose VGAE-CCI, a deep learning framework based on the Variational Graph Autoencoder, capable of identifying cell-cell communication across multiple tissue layers. Additionally, this model can be applied to spatial transcriptomics data with missing or partially incomplete data and can clustered cells at single-cell resolution based on spatial encoding information within complex tissues, thereby enabling more accurate inference of cell-cell communication. Finally, we tested our method on six datasets and compared it with other state of art methods for predicting cell-cell communication. Our method outperformed other methods across multiple metrics, demonstrating its efficiency and reliability in predicting cell-cell communication.

CellMsg: graph convolutional networks for ligand-receptor-mediated cell-cell communication analysis

The role of cell-cell communications (CCCs) is increasingly recognized as being important to differentiation, invasion, metastasis, and drug resistance in tumoral tissues. Developing CCC inference methods using traditional experimental methods are time-consuming, labor-intensive, cannot handle large amounts of data. To facilitate inference of CCCs, we proposed a computational framework, called CellMsg, which involves two primary steps: identifying ligand-receptor interactions (LRIs) and measuring the strength of LRIs-mediated CCCs. Specifically, CellMsg first identifies high-confident LRIs based on multimodal features of ligands and receptors and graph convolutional networks. Then, CellMsg measures the strength of intercellular communication by combining the identified LRIs and single-cell RNA-seq data using a three-point estimation method. Performance evaluation on four benchmark LRI datasets by five-fold cross validation demonstrated that CellMsg accurately captured the relationships between ligands and receptors, resulting in the identification of high-confident LRIs. Compared with other methods of identifying LRIs, CellMsg has better prediction performance and robustness. Furthermore, the LRIs identified by CellMsg were successfully validated through molecular docking. Finally, we examined the overlap of LRIs between CellMsg and five other classical CCC databases, as well as the intercellular crosstalk among seven cell types within a human melanoma tissue. In summary, CellMsg establishes a complete, reliable, and well-organized LRI database and an effective CCC strength evaluation method for each single-cell RNA-seq data. It provides a computational tool allowing researchers to decipher intercellular communications. CellMsg is freely available at https://github.com/pengsl-lab/CellMsg.

An iron-rich subset of macrophages promotes tumor growth through a Bach1-Ednrb axis

We define a subset of macrophages in the tumor microenvironment characterized by high intracellular iron and enrichment of heme and iron metabolism genes. These iron-rich tumor-associated macrophages (iTAMs) supported angiogenesis and immunosuppression in the tumor microenvironment and were conserved between mice and humans. iTAMs comprise two additional subsets based on gene expression profile and location-perivascular (pviTAM) and stromal (stiTAM). We identified the endothelin receptor type B (Ednrb) as a specific marker of iTAMs and found myeloid-specific deletion of Ednrb to reduce tumor growth and vascular density. Further studies identified the transcription factor Bach1 as a repressor of the iTAM transcriptional program, including Ednrb expression. Heme is a known inhibitor of Bach1, and, correspondingly, heme exposure induced Ednrb and iTAM signature genes in macrophages. Thus, iTAMs are a distinct macrophage subset regulated by the transcription factor Bach1 and characterized by Ednrb-mediated immunosuppressive and angiogenic functions.

Next-generation spatial transcriptomics: unleashing the power to gear up translational oncology

The growing advances in spatial transcriptomics (ST) stand as the new frontier bringing unprecedented influences in the realm of translational oncology. This has triggered systemic experimental design, analytical scope, and depth alongside with thorough bioinformatics approaches being constantly developed in the last few years. However, harnessing the power of spatial biology and streamlining an array of ST tools to achieve designated research goals are fundamental and require real-world experiences. We present a systemic review by updating the technical scope of ST across different principal basis in a timeline manner hinting on the generally adopted ST techniques used within the community. We also review the current progress of bioinformatic tools and propose in a pipelined workflow with a toolbox available for ST data exploration. With particular interests in tumor microenvironment where ST is being broadly utilized, we summarize the up-to-date progress made via ST-based technologies by narrating studies categorized into either mechanistic elucidation or biomarker profiling (translational oncology) across multiple cancer types and their ways of deploying the research through ST. This updated review offers as a guidance with forward-looking viewpoints endorsed by many high-resolution ST tools being utilized to disentangle biological questions that may lead to clinical significance in the future.

Dimension reduction, cell clustering, and cell-cell communication inference for single-cell transcriptomics with DcjComm

Advances in single-cell transcriptomics provide an unprecedented opportunity to explore complex biological processes. However, computational methods for analyzing single-cell transcriptomics still have room for improvement especially in dimension reduction, cell clustering, and cell-cell communication inference. Herein, we propose a versatile method, named DcjComm, for comprehensive analysis of single-cell transcriptomics. DcjComm detects functional modules to explore expression patterns and performs dimension reduction and clustering to discover cellular identities by the non-negative matrix factorization-based joint learning model. DcjComm then infers cell-cell communication by integrating ligand-receptor pairs, transcription factors, and target genes. DcjComm demonstrates superior performance compared to state-of-the-art methods.

LIANA+ provides an all-in-one framework for cell-cell communication inference

The growing availability of single-cell and spatially resolved transcriptomics has led to the development of many approaches to infer cell-cell communication, each capturing only a partial view of the complex landscape of intercellular signalling. Here we present LIANA+, a scalable framework built around a rich knowledge base to decode coordinated inter- and intracellular signalling events from single- and multi-condition datasets in both single-cell and spatially resolved data. By extending and unifying established methodologies, LIANA+ provides a comprehensive set of synergistic components to study cell-cell communication via diverse molecular mediators, including those measured in multi-omics data. LIANA+ is accessible at https://github.com/saezlab/liana-py with extensive vignettes ( https://liana-py.readthedocs.io/ ) and provides an all-in-one solution to intercellular communication inference.

Deciphering cell-cell communication at single-cell resolution for spatial transcriptomics with subgraph-based graph attention network

The inference of cell-cell communication (CCC) is crucial for a better understanding of complex cellular dynamics and regulatory mechanisms in biological systems. However, accurately inferring spatial CCCs at single-cell resolution remains a significant challenge. To address this issue, we present a versatile method, called DeepTalk, to infer spatial CCC at single-cell resolution by integrating single-cell RNA sequencing (scRNA-seq) data and spatial transcriptomics (ST) data. DeepTalk utilizes graph attention network (GAT) to integrate scRNA-seq and ST data, which enables accurate cell-type identification for single-cell ST data and deconvolution for spot-based ST data. Then, DeepTalk can capture the connections among cells at multiple levels using subgraph-based GAT, and further achieve spatially resolved CCC inference at single-cell resolution. DeepTalk achieves excellent performance in discovering meaningful spatial CCCs on multiple cross-platform datasets, which demonstrates its superior ability to dissect cellular behavior within intricate biological processes.

Cell-cell communication: new insights and clinical implications

Multicellular organisms are composed of diverse cell types that must coordinate their behaviors through communication. Cell-cell communication (CCC) is essential for growth, development, differentiation, tissue and organ formation, maintenance, and physiological regulation. Cells communicate through direct contact or at a distance using ligand-receptor interactions. So cellular communication encompasses two essential processes: cell signal conduction for generation and intercellular transmission of signals, and cell signal transduction for reception and procession of signals. Deciphering intercellular communication networks is critical for understanding cell differentiation, development, and metabolism. First, we comprehensively review the historical milestones in CCC studies, followed by a detailed description of the mechanisms of signal molecule transmission and the importance of the main signaling pathways they mediate in maintaining biological functions. Then we systematically introduce a series of human diseases caused by abnormalities in cell communication and their progress in clinical applications. Finally, we summarize various methods for monitoring cell interactions, including cell imaging, proximity-based chemical labeling, mechanical force analysis, downstream analysis strategies, and single-cell technologies. These methods aim to illustrate how biological functions depend on these interactions and the complexity of their regulatory signaling pathways to regulate crucial physiological processes, including tissue homeostasis, cell development, and immune responses in diseases. In addition, this review enhances our understanding of the biological processes that occur after cell-cell binding, highlighting its application in discovering new therapeutic targets and biomarkers related to precision medicine. This collective understanding provides a foundation for developing new targeted drugs and personalized treatments.

scHyper: reconstructing cell-cell communication through hypergraph neural networks

Cell-cell communications is crucial for the regulation of cellular life and the establishment of cellular relationships. Most approaches of inferring intercellular communications from single-cell RNA sequencing (scRNA-seq) data lack a comprehensive global network view of multilayered communications. In this context, we propose scHyper, a new method that can infer intercellular communications from a global network perspective and identify the potential impact of all cells, ligand, and receptor expression on the communication score. scHyper designed a new way to represent tripartite relationships, by extracting a heterogeneous hypergraph that includes the source (ligand expression), the target (receptor expression), and the relevant ligand-receptor (L-R) pairs. scHyper is based on hypergraph representation learning, which measures the degree of match between the intrinsic attributes (static embeddings) of nodes and their observed behaviors (dynamic embeddings) in the context (hyperedges), quantifies the probability of forming hyperedges, and thus reconstructs the cell-cell communication score. Additionally, to effectively mine the key mechanisms of signal transmission, we collect a rich dataset of multisubunit complex L-R pairs and propose a nonparametric test to determine significant intercellular communications. Comparing with other tools indicates that scHyper exhibits superior performance and functionality. Experimental results on the human tumor microenvironment and immune cells demonstrate that scHyper offers reliable and unique capabilities for analyzing intercellular communication networks. Therefore, we introduced an effective strategy that can build high-order interaction patterns, surpassing the limitations of most methods that can only handle low-order interactions, thus more accurately interpreting the complexity of intercellular communications.

Crosstalk between tumor and microenvironment: Insights from spatial transcriptomics

Cancer is a dynamic disease, and clonal heterogeneity plays a fundamental role in tumor development, progression, and resistance to therapies. Single-cell and spatial multimodal technologies can provide a high-resolution molecular map of underlying genomic, epigenomic, and transcriptomic alterations involved in inter- and intra-tumor heterogeneity and interactions with the microenvironment. In this review, we provide a perspective on factors driving cancer heterogeneity, tumor evolution, and clonal states. We briefly describe spatial transcriptomic technologies and summarize recent literature that sheds light on the dynamical interactions between tumor states, cell-to-cell communication, and remodeling local microenvironment.

Predicting intercellular communication based on metabolite-related ligand-receptor interactions with MRCLinkdb

BACKGROUND

Metabolite-associated cell communications play critical roles in maintaining human biological function. However, most existing tools and resources focus only on ligand-receptor interaction pairs where both partners are proteinaceous, neglecting other non-protein molecules. To address this gap, we introduce the MRCLinkdb database and algorithm, which aggregates and organizes data related to non-protein L-R interactions in cell-cell communication, providing a valuable resource for predicting intercellular communication based on metabolite-related ligand-receptor interactions.

RESULTS

Here, we manually curated the metabolite-ligand-receptor (ML-R) interactions from the literature and known databases, ultimately collecting over 790 human and 670 mouse ML-R interactions. Additionally, we compiled information on over 1900 enzymes and 260 transporter entries associated with these metabolites. We developed Metabolite-Receptor based Cell Link Database (MRCLinkdb) to store these ML-R interactions data. Meanwhile, the platform also offers extensive information for presenting ML-R interactions, including fundamental metabolite information and the overall expression landscape of metabolite-associated gene sets (such as receptor, enzymes, and transporter proteins) based on single-cell transcriptomics sequencing (covering 35 human and 26 mouse tissues, 52 human and 44 mouse cell types) and bulk RNA-seq/microarray data (encompassing 62 human and 39 mouse tissues). Furthermore, MRCLinkdb introduces a web server dedicated to the analysis of intercellular communication based on ML-R interactions. MRCLinkdb is freely available at https://www.cellknowledge.com.cn/mrclinkdb/ .

CONCLUSIONS

In addition to supplementing ligand-receptor databases, MRCLinkdb may provide new perspectives for decoding the intercellular communication and advancing related prediction tools based on ML-R interactions.

A count-based model for delineating cell-cell interactions in spatial transcriptomics data

MOTIVATION

Cell-cell interactions (CCIs) consist of cells exchanging signals with themselves and neighboring cells by expressing ligand and receptor molecules and play a key role in cellular development, tissue homeostasis, and other critical biological functions. Since direct measurement of CCIs is challenging, multiple methods have been developed to infer CCIs by quantifying correlations between the gene expression of the ligands and receptors that mediate CCIs, originally from bulk RNA-sequencing data and more recently from single-cell or spatially resolved transcriptomics (SRT) data. SRT has a particular advantage over single-cell approaches, since ligand-receptor correlations can be computed between cells or spots that are physically close in the tissue. However, the transcript counts of individual ligands and receptors in SRT data are generally low, complicating the inference of CCIs from expression correlations.

RESULTS

We introduce Copulacci, a count-based model for inferring CCIs from SRT data. Copulacci uses a Gaussian copula to model dependencies between the expression of ligands and receptors from nearby spatial locations even when the transcript counts are low. On simulated data, Copulacci outperforms existing CCI inference methods based on the standard Spearman and Pearson correlation coefficients. Using several real SRT datasets, we show that Copulacci discovers biologically meaningful ligand-receptor interactions that are lowly expressed and undiscoverable by existing CCI inference methods.

AVAILABILITY AND IMPLEMENTATION

Copulacci is implemented in Python and available at https://github.com/raphael-group/copulacci.

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.

PathFinder: a novel graph transformer model to infer multi-cell intra- and inter-cellular signaling pathways and communications

Recently, large-scale scRNA-seq datasets have been generated to understand the complex signaling mechanisms within the microenvironment of Alzheimer's Disease (AD), which are critical for identifying novel therapeutic targets and precision medicine. However, the background signaling networks are highly complex and interactive. It remains challenging to infer the core intra- and inter-multi-cell signaling communication networks using scRNA-seq data. In this study, we introduced a novel graph transformer model, PathFinder, to infer multi-cell intra- and inter-cellular signaling pathways and communications among multi-cell types. Compared with existing models, the novel and unique design of PathFinder is based on the divide-and-conquer strategy. This model divides complex signaling networks into signaling paths, which are then scored and ranked using a novel graph transformer architecture to infer intra- and inter-cell signaling communications. We evaluated the performance of PathFinder using two scRNA-seq data cohorts. The first cohort is an APOE4 genotype-specific AD, and the second is a human cirrhosis cohort. The evaluation confirms the promising potential of using PathFinder as a general signaling network inference model.

Reading the palimpsest of cell interactions: What questions may we ask of the data?

Biological function depends on the composition and structure of the organism, the latter describing the organization of interactions between parts. While cells in multicellular organisms are capable of a remarkable degree of autonomy, most functions do require cell communication: the coordination of functions (growth, differentiation, and apoptosis), the compartmentalization of cellular processes, and the integration of cells into higher levels of structural organization. A wealth of data on putative cell interactions has become available, yet its biological interpretation depends on our expectations about the structure of interaction networks. Here, we attempt to formulate basic questions to ask when interpreting cell interaction data. We build on the understanding that cells fulfill two general functions: the integrity-maintaining and the organismal service function. We derive the expected patterns of cell interactions considering two intertwined aspects: the functional and the evolutionary. Based on these, we propose guidelines for analysis and interpretation of transcriptional cell-interactome data.

Genetic Tools for Cell Lineage Tracing and Profiling Developmental Trajectories in the Skin

The epidermis is the body's first line of protection against dehydration and pathogens, continually regenerating the outermost protective skin layers throughout life. During both embryonic development and wound healing, epidermal stem and progenitor cells must respond to external stimuli and insults to build, maintain, and repair the cutaneous barrier. Recent advances in CRISPR-based methods for cell lineage tracing have remarkably expanded the potential for experiments that track stem and progenitor cell proliferation and differentiation over the course of tissue and even organismal development. Additional tools for DNA-based recording of cellular signaling cues promise to deepen our understanding of the mechanisms driving normal skin morphogenesis and response to stressors as well as the dysregulation of cell proliferation and differentiation in skin diseases and cancer. In this review, we highlight cutting-edge methods for cell lineage tracing, including in organoids and model organisms, and explore how cutaneous biology researchers might leverage these techniques to elucidate the developmental programs that support the regenerative capacity and plasticity of the skin.

Signals that control MAIT cell function in healthy and inflamed human tissues

Mucosal-associated invariant T (MAIT) cells have a semi-invariant T-cell receptor that allows recognition of antigen in the context of the MHC class I-related (MR1) protein. Metabolic intermediates of the riboflavin synthesis pathway have been identified as MR1-restricted antigens with agonist properties. As riboflavin synthesis occurs in many bacterial species, but not human cells, it has been proposed that the main purpose of MAIT cells is antibacterial surveillance and protection. The majority of human MAIT cells secrete interferon-gamma (IFNg) upon activation, while some MAIT cells in tissues can also express IL-17. Given that MAIT cells are present in human barrier tissues colonized by a microbiome, MAIT cells must somehow be able to distinguish colonization from infection to ensure effector functions are only elicited when necessary. Importantly, MAIT cells have additional functional properties, including the potential to contribute to restoring tissue homeostasis by expression of CTLA-4 and secretion of the cytokine IL-22. A recent study provided compelling data indicating that the range of human MAIT cell functional properties is explained by plasticity rather than distinct lineages. This further underscores the necessity to better understand how different signals regulate MAIT cell function. In this review, we highlight what is known in regards to activating and inhibitory signals for MAIT cells with a specific focus on signals relevant to healthy and inflamed tissues. We consider the quantity, quality, and the temporal order of these signals on MAIT cell function and discuss the current limitations of computational tools to extrapolate which signals are received by MAIT cells in human tissues. Using lessons learned from conventional CD8 T cells, we also discuss how TCR signals may integrate with cytokine signals in MAIT cells to elicit distinct functional states.

SEnSCA: Identifying possible ligand-receptor interactions and its application in cell-cell communication inference

Multicellular organisms have dense affinity with the coordination of cellular activities, which severely depend on communication across diverse cell types. Cell-cell communication (CCC) is often mediated via ligand-receptor interactions (LRIs). Existing CCC inference methods are limited to known LRIs. To address this problem, we developed a comprehensive CCC analysis tool SEnSCA by integrating single cell RNA sequencing and proteome data. SEnSCA mainly contains potential LRI acquisition and CCC strength evaluation. For acquiring potential LRIs, it first extracts LRI features and reduces the feature dimension, subsequently constructs negative LRI samples through K-means clustering, finally acquires potential LRIs based on Stacking ensemble comprising support vector machine, 1D-convolutional neural networks and multi-head attention mechanism. During CCC strength evaluation, SEnSCA conducts LRI filtering and then infers CCC by combining the three-point estimation approach and single cell RNA sequencing data. SEnSCA computed better precision, recall, accuracy, F1 score, AUC and AUPR under most of conditions when predicting possible LRIs. To better illustrate the inferred CCC network, SEnSCA provided three visualization options: heatmap, bubble diagram and network diagram. Its application on human melanoma tissue demonstrated its reliability in CCC detection. In summary, SEnSCA offers a useful CCC inference tool and is freely available at https://github.com/plhhnu/SEnSCA.

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.

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.

Identifying potential ligand-receptor interactions based on gradient boosted neural network and interpretable boosting machine for intercellular communication analysis

Cell-cell communication is essential to many key biological processes. Intercellular communication is generally mediated by ligand-receptor interactions (LRIs). Thus, building a comprehensive and high-quality LRI resource can significantly improve intercellular communication analysis. Meantime, due to lack of a "gold standard" dataset, it remains a challenge to evaluate LRI-mediated intercellular communication results. Here, we introduce CellGiQ, a high-confident LRI prediction framework for intercellular communication analysis. Highly confident LRIs are first inferred by LRI feature extraction with BioTriangle, LRI selection using LightGBM, and LRI classification based on ensemble of gradient boosted neural network and interpretable boosting machine. Subsequently, known and identified high-confident LRIs are filtered by combining single-cell RNA sequencing (scRNA-seq) data and further applied to intercellular communication inference through a quartile scoring strategy. To validation the predictions, CellGiQ exploited several evaluation strategies: using AUC and AUPR, it surpassed six competing LRI prediction models on four LRI datasets; through Venn diagrams and molecular docking, its predicted LRIs were validated by five other popular intercellular communication inference methods; based on the overlapping LRIs, it computed high Jaccard index with six other state-of-the-art intercellular communication prediction tools within human HNSCC tissues; by comparing with classical models and literature retrieve, its inferred HNSCC-related intercellular communication results was further validated. The novelty of this study is to identify high-confident LRIs based on machine learning as well as design several LRI validation ways, providing reference for computational LRI prediction. CellGiQ provides an open-source and useful tool to decompose LRI-mediated intercellular communication at single cell resolution. CellGiQ is freely available at https://github.com/plhhnu/CellGiQ.

Role of Artificial Intelligence in Drug Discovery and Target Identification in Cancer

Drug discovery and development (DDD) is a highly complex process that necessitates precise monitoring and extensive data analysis at each stage. Furthermore, the DDD process is both timeconsuming and costly. To tackle these concerns, artificial intelligence (AI) technology can be used, which facilitates rapid and precise analysis of extensive datasets within a limited timeframe. The pathophysiology of cancer disease is complicated and requires extensive research for novel drug discovery and development. The first stage in the process of drug discovery and development involves identifying targets. Cell structure and molecular functioning are complex due to the vast number of molecules that function constantly, performing various roles. Furthermore, scientists are continually discovering novel cellular mechanisms and molecules, expanding the range of potential targets. Accurately identifying the correct target is a crucial step in the preparation of a treatment strategy. Various forms of AI, such as machine learning, neural-based learning, deep learning, and network-based learning, are currently being utilised in applications, online services, and databases. These technologies facilitate the identification and validation of targets, ultimately contributing to the success of projects. This review focuses on the different types and subcategories of AI databases utilised in the field of drug discovery and target identification for cancer.

STACCato: Supervised Tensor Analysis tool for studying Cell-cell Communication using scRNA-seq data across multiple samples and conditions

Research on cell-cell communication (CCC) is crucial for understanding biology and diseases. Many existing CCC inference tools neglect potential confounders, such as batch and demographic variables, when analyzing multi-sample, multi-condition scRNA-seq datasets. To address this significant gap, we introduce STACCato, a Supervised Tensor Analysis tool for studying Cell-cell Communication, that identifies CCC events and estimates the effects of biological conditions (e.g., disease status, tissue types) on such events, while adjusting for potential confounders. Application of STACCato to both simulated data and real scRNA-seq data of lupus and autism studies demonstrate that incorporating sample-level variables into CCC inference consistently provides more accurate estimations of disease effects and cell type activity patterns than existing methods that ignore sample-level variables. A computational tool implementing the STACCato framework is available on GitHub.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

TALKIEN: crossTALK IntEraction Network. A web-based tool for deciphering molecular communication through ligand-receptor interactions

Molecular crosstalk, the dialogue between different cell types, is attracting more attention in cancer research. On the one hand, the communication between tumor and non-tumor cells in the microenvironment or between different tumor clones has influential consequences for the progression and spread of tumors and response to treatment. On the other hand, novel techniques such as single-cell sequencing or spatial transcriptomics provide detailed information that needs to be interpreted. TALKIEN: crossTALK IntEraction Network is a simple and intuitive online R/shiny application to visualize molecular crosstalk information through the construction and analysis of a protein-protein interaction network. Taking two or more lists of genes or proteins as input, which are representative of cell lineages, TALKIEN extracts information about ligand-receptor interactions, builds a network and analyzes it using systems biology techniques such as centrality measures and component analysis, among others. Moreover, it expands the network displaying pathways downstream receptors. The application allows users to select different graphical layouts, performs functional analysis and gives information about drugs targeting receptors. In conclusion, TALKIEN allows users to detect ligand-receptor interactions generating new in silico predictions of cell-cell communication thus providing a translational rationale for future experiments. It is freely available at https://www.odap-ico.org/talkien.

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.

ESICCC as a systematic computational framework for evaluation, selection, and integration of cell-cell communication inference methods

Cell-cell communication (CCC) is critical for determining cell fates and functions in multicellular organisms. With the advent of single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), an increasing number of CCC inference methods have been developed. Nevertheless, a thorough comparison of their performances is yet to be conducted. To fill this gap, we developed a systematic benchmark framework called ESICCC to evaluate 18 ligand-receptor (LR) inference methods and five ligand/receptor-target inference methods using a total of 116 data sets, including 15 ST data sets, 15 sets of cell line perturbation data, two sets of cell type-specific expression/proteomics data, and 84 sets of sampled or unsampled scRNA-seq data. We evaluated and compared the agreement, accuracy, robustness, and usability of these methods. Regarding accuracy evaluation, RNAMagnet, CellChat, and scSeqComm emerge as the three best-performing methods for intercellular ligand-receptor inference based on scRNA-seq data, whereas stMLnet and HoloNet are the best methods for predicting ligand/receptor-target regulation using ST data. To facilitate the practical applications, we provide a decision-tree-style guideline for users to easily choose best tools for their specific research concerns in CCC inference, and develop an ensemble pipeline CCCbank that enables versatile combinations of methods and databases. Moreover, our comparative results also uncover several critical influential factors for CCC inference, such as prior interaction information, ligand-receptor scoring algorithm, intracellular signaling complexity, and spatial relationship, which may be considered in the future studies to advance the development of new methodologies.

Decoding functional cell-cell communication events by multi-view graph learning on spatial transcriptomics

Cell-cell communication events (CEs) are mediated by multiple ligand-receptor (LR) pairs. Usually only a particular subset of CEs directly works for a specific downstream response in a particular microenvironment. We name them as functional communication events (FCEs) of the target responses. Decoding FCE-target gene relations is: important for understanding the mechanisms of many biological processes, but has been intractable due to the mixing of multiple factors and the lack of direct observations. We developed a method HoloNet for decoding FCEs using spatial transcriptomic data by integrating LR pairs, cell-type spatial distribution and downstream gene expression into a deep learning model. We modeled CEs as a multi-view network, developed an attention-based graph learning method to train the model for generating target gene expression with the CE networks, and decoded the FCEs for specific downstream genes by interpreting trained models. We applied HoloNet on three Visium datasets of breast cancer and liver cancer. The results detangled the multiple factors of FCEs by revealing how LR signals and cell types affect specific biological processes, and specified FCE-induced effects in each single cell. We conducted simulation experiments and showed that HoloNet is more reliable on LR prioritization in comparison with existing methods. HoloNet is a powerful tool to illustrate cell-cell communication landscapes and reveal vital FCEs that shape cellular phenotypes. HoloNet is available as a Python package at https://github.com/lhc17/HoloNet.