Predicting cell-to-cell communication networks using NATMI

Predictive Biomarkers for Immune Checkpoint Inhibitor Therapy in Advanced Melanomas

Biomarkers capable of predicting adverse melanoma patient responses to immune checkpoint inhibitor (ICI) therapies are an unmet need. Clinical biomarkers are largely prognostic and current response guidelines do not reflect the complex tumor-immune cell interaction dynamics attributed to ICI therapies. Validation of enhanced predictive biomarkers is dependent upon adoption of novel spatial imaging platforms capable of quantifying immune checkpoint receptor-ligand interactions within the tumor microenvironment. Assessments of these interactions at multiple points during neoadjuvant ICI regimens would inform biomarker selection based on changes in receptor-ligand interactions that best correlate with patient survival.

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.

Ligand-receptor dynamics in heterophily-aware graph neural networks for enhanced cell type prediction from single-cell RNA-seq data

Graph Neural Networks (GNNs) have emerged as powerful tools for analyzing structured data, particularly in domains where relationships and interactions between entities are key. By leveraging the inherent graph structure in datasets, GNNs excel in capturing complex dependencies and patterns that traditional neural networks might miss. This advantage is especially pronounced in the field of computational biology, where the intricate connections between biological entities play a crucial role. In this context, Our work explores the application of GNNs to single-cell RNA sequencing (scRNA-seq) data, a domain characterized by complex and heterogeneous relationships. By extracting ligand-receptor (L-R) associations from LIANA and constructing Cell-Cell association networks with varying edge homophily ratios, based on L-R information, we enhance the biological relevance and accuracy of depicting cellular communication pathways. While standard GNN models like Graph Convolutional Networks (GCN), GraphSAGE, Graph Attention Networks (GAT), and MixHop often assume homophily (similar nodes are more likely to be connected), this assumption does not always hold in biological networks. To address this, we explore advanced graph neural network methods, such asH 2 Graph Convolutional Networks and Gated Bi-Kernel GNNs (GBK-GNN), that are specifically designed to handle heterophilic data. Our study spans across six diverse datasets, enabling a thorough comparison between heterophily-aware GNNs and traditional homophily-assuming models, including Multi-Layer Perceptrons, which disregards graph structure entirely. Our findings highlight the importance of considering data-specific characteristics in GNN applications, demonstrating that heterophily-focused methods can effectively decipher the complex patterns within scRNA-seq data. By integrating multi-omics data, including gene expression profiles and L-R interactions, we pave the way for more accurate and insightful analyses in computational biology, offering a more comprehensive understanding of cellular environments and interactions.

Single-cell and spatial analyses reveal the effect of VSIG4+S100A10+TAMs on the immunosuppression of glioblastoma and anti-PD-1 immunotherapy

Therapeutic strategies aiming at the tumor immune microenvironment (TIME) hold promise for glioblastoma (GBM) treatment. However, adjuvant immunotherapies targeting checkpoint inhibitors just prove effective for a selected group of GBM patients. The extensive involvement of GBM-associated macrophages highlights their potential role in tumor behavior. In-depth exploration of the impact of macrophages on the efficacy of immunotherapy is crucial for enhancing treatment outcomes. In this study, we conducted a comprehensive analysis using bulk RNA-seq, single-cell RNA sequencing (scRNA-seq), and spatial transcriptomics to explore the heterogeneity of tumor-associated macrophages (TAMs) in GBM. Flow cytometry was employed to investigate the effects of VSIG4 on TAM phenotypes, and co-culture cellular assays were performed to evaluate its contribution to GBM malignancy. Integrating 16 patient samples, we examined the immunological significance of VSIG4+S100A10+TAMs. VSIG4 expression on macrophages is significantly upregulated and correlated with the TIME, promoting the polarization of macrophages towards M2 and facilitating GBM progression. Spatial transcriptomics and human samples multiplex immunofluorescence (mIF) confirmed the co-localization of VSIG4+S100A10+TAMs with various T cells, resulting in the inhibition of T cell immune responses and a reduction in anti-tumor immunity. Our findings demonstrate for the first time that VSIG4+S100A10+TAM is an independent prognostic indicator of poor outcome for GBM and markedly accumulates in patients exhibiting non-responsiveness to anti-PD-1 immunotherapy. Targeting this specific bifunctional subgroup can potentially open up new avenues for the immunotherapy of GBM.

scSDNE: A semi-supervised method for inferring cell-cell interactions based on graph embedding

As a fundamental characteristic of multicellular organisms, cell-cell communication is achieved through ligand-receptor (L-R) interactions, enabling the exchange of information and revealing the diversity of biological processes and cellular functions. To gain a comprehensive understanding of these complex interaction mechanisms, we constructed a manually curated L-R interaction database and developed a semi-supervised graph embedding model called scSDNE for inferring cell-cell interactions mediated by L-R interactions. scSDNE model utilizes the power of deep learning to map genes from interacting cells into a shared latent space, allowing for a nuanced representation of their relationships. Leveraging the prior information provided by database, scSDNE can infer significant L-R pairs involved in intercellular communication. Experiments on real single-cell RNA sequencing (scRNA-seq) datasets demonstrate that our method detects interactions with a high degree of reliability compared with other methods. More importantly, the model integrates gene regulation information within cells to enhance the accuracy and biological interpretability of the inferences. Our method provides a more comprehensive view of cell-cell interactions, offering new insights into complex intercellular communication.

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.

Exploring the Role of T-Cell Metabolism in Modulating Immunotherapy Efficacy for Non-Small Cell Lung Cancer Based on Clustering

BACKGROUND

Immunotherapy, especially immune checkpoint blockade (ICB) therapy, has demonstrated noteworthy advancements in the realm of non-small cell lung cancer (NSCLC). However, the efficacy of ICB therapy is limited to a small subset of patients with NSCLC, and the underlying mechanisms remain poorly understood.

STUDY DESIGN AND DISCOVERIES

In this study, we conducted a comprehensive investigation of the metabolic profiles of infiltrating T cells in NSCLC tumors and revealed the metabolic heterogeneity, which associated with the prognosis of ICB therapy, in three T-cell subtypes. After metabolic clustering, we split these metabolic clusters into two groups: Nonresponse-associated (NR) clusters that enriched with cells from nonresponders, and response-associated (R) clusters that not belonging to NR clusters. Then, we elucidated their metabolic differences and specific functions. Notably, we discovered HSPA1A was significantly downregulated in NR clusters of all three T-cell subtypes. In addition, leveraging single-cell T-cell receptor sequencing data and pseudotime series analysis, we revealed the reciprocal interconversion between R and NR metabolic clusters within the same T-cell clone. This suggests a potential metabolic reprogramming capability of T cells. Furthermore, through the analysis of intercellular communication, we identified the specific intercellular signaling in the R clusters, which might promote the activation and regulation of signal transduction pathways that affect the prognosis of ICB therapy.

CONCLUSION

In conclusion, our study offers substantial insights into the mechanisms of relationships between T-cell metabolisms and ICB therapy outcomes, shedding light on the mechanism of immunotherapy efficacy in patients with NSCLC. Such investigations will contribute to overcoming treatment resistance.

Deciphering immune predictors of immunotherapy response: A multiomics approach at the pan-cancer level

Immune checkpoint blockade (ICB) therapy has transformed cancer treatment, yet many patients fail to respond. Employing single-cell multiomics, we unveil T cell dynamics influencing ICB response across 480 pan-cancer and 27 normal tissue samples. We identify four immunotherapy response-associated T cells (IRATs) linked to responsiveness or resistance and analyze their pseudotemporal patterns, regulatory mechanisms, and T cell receptor clonal expansion profiles specific to each response. Notably, transforming growth factor β1 (TGF-β1)+ CD4+ and Temra CD8+ T cells negatively correlate with therapy response, in stark contrast to the positive response associated with CXCL13+ CD4+ and CD8+ T cells. Validation with a cohort of 23 colorectal cancer (CRC) samples confirms the significant impact of TGF-β1+ CD4+ and CXCL13+ CD4+ and CD8+ T cells on ICB efficacy. Our study highlights the effectiveness of single-cell multiomics in pinpointing immune markers predictive of immunotherapy outcomes, providing an important resource for crafting targeted immunotherapies for successful ICB treatment across cancers.

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.

Lignature: A Comprehensive Database of Ligand Signatures to Predict Cell-Cell Communication

Ligand-receptor interactions mediate intercellular communication, inducing transcriptional changes that regulate physiological and pathological processes. Ligand-induced transcriptomic signatures can be used to predict active ligands; however, the absence of a comprehensive set of ligand-response signatures has limited their practical application in predicting ligand-receptor interactions. To bridge this gap, we developed Lignature, a curated database encompassing intracellular transcriptomic signatures for 362 human ligands, significantly expanding the repertoire of ligands with available intracellular response signatures. Lignature compiles signatures from published transcriptomic datasets and established resources such as CytoSig and ImmuneDictionary, generating both gene- and pathway-based signatures for each ligand. We applied Lignature to predict active ligands driving transcriptomic changes in controlled in vitro experiments and real-world single-cell sequencing datasets. Lignature outperformed existing methods such as NicheNet, achieving higher accuracy in identifying active ligands at both the gene and pathway levels. These results establish Lignature as a robust platform for ligand signaling inference, providing a powerful tool to explore ligand-receptor interactions across diverse experimental and physiological contexts.

Expression of ENL YEATS domain tumor mutations in nephrogenic or stromal lineage impairs kidney development

Recurrent gain-of-function mutations in the histone reader protein ENL have been identified in Wilms tumor, the most prevalent pediatric kidney cancer. However, their pathological significance in kidney development and tumorigenesis in vivo remains elusive. Here, we generate mouse models mimicking ENL tumor (ENLT) mutations and show that heterozygous mutant expression in Six2+ nephrogenic or Foxd1+ stromal lineages leads to severe, lineage-specific kidney defects, both resulting in neonatal lethality. Six2-ENLT mutant kidneys display compromised cap mesenchyme, scant nephron tubules, and cystic glomeruli, indicative of premature progenitor commitment and blocked differentiation. Bulk and spatial transcriptomic analyses reveal aberrant activation of Hox and Wnt signaling genes in mutant nephrogenic cells. In contrast, Foxd1-ENLT mutant kidneys exhibit expansion in renal capsule and cap mesenchyme, with dysregulated stromal gene expression affecting stroma-epithelium crosstalk. Our findings uncover distinct pathways through which ENL mutations disrupt nephrogenesis, providing a foundation for further investigations into their role in tumorigenesis.

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.

Cellular interactions within the immune microenvironment underpins resistance to cell cycle inhibition in breast cancers

Immune evasion by cancer cells involves reshaping the tumor microenvironment (TME) via communication with non-malignant cells. However, resistance-promoting interactions during treatment remain lesser known. Here we examine the composition, communication, and phenotypes of tumor-associated cells in serial biopsies from stage II and III high-risk estrogen receptor positive (ER+ ) breast cancers of patients receiving endocrine therapy (letrozole) as single agent or in combination with ribociclib, a CDK4/6-targeting cell cycle inhibitor. Single-cell RNA sequencing analyses on longitudinally collected samples show that in tumors overcoming the growth suppressive effects of ribociclib, first cancer cells upregulate cytokines and growth factors that stimulate immune-suppressive myeloid differentiation, resulting in reduced myeloid cell- CD8 + T-cell crosstalk via IL-15/18 signaling. Subsequently, tumors growing during treatment show diminished T-cell activation and recruitment. In vitro, ribociclib does not only inhibit cancer cell growth but also T cell proliferation and activation upon co-culturing. Exogenous IL-15 improves CDK4/6 inhibitor efficacy by augmenting T-cell proliferation and cancer cell killing by T cells. In summary, response to ribociclib in stage II and III high-risk ER +  breast cancer depends on the composition, activation phenotypes and communication network of immune cells.

Mathematically mapping the network of cells in the tumor microenvironment

Cell-cell interaction (CCI) networks are key to understanding disease progression and treatment response. However, existing methods for inferring these networks often aggregate data across patients or focus on cell-type level interactions, providing a generalized overview but overlooking patient heterogeneity and local network structures. To address this, we introduce "random cell-cell interaction generator" (RaCInG), a model based on random graphs to derive personalized networks leveraging prior knowledge on ligand-receptor interactions and bulk RNA sequencing data. We applied RaCInG to 8,683 cancer patients to extract 643 network features related to the tumor microenvironment and unveiled associations with immune response and subtypes, enabling prediction and explanation of immunotherapy responses. RaCInG demonstrated robustness and showed consistencies with state-of-the-art methods. Our findings highlight RaCInG's potential to elucidate patient-specific network dynamics, offering insights into cancer biology and treatment responses. RaCInG is poised to advance our understanding of complex CCI s in cancer and other biomedical domains.

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.

Genome-scale spatial mapping of the Hodgkin lymphoma microenvironment identifies tumor cell survival factors

A key challenge in cancer research is to identify the secreted factors that contribute to tumor cell survival. Nowhere is this more evident than in Hodgkin lymphoma, where malignant Hodgkin Reed Sternberg (HRS) cells comprise only 1-5% of the tumor mass, the remainder being infiltrating immune cells that presumably are required for the survival of the HRS cells. Until now, there has been no way to characterize the complex Hodgkin lymphoma tumor microenvironment at genome scale. Here, we performed genome-wide transcriptional profiling with spatial and single-cell resolution. We show that the neighborhood surrounding HRS cells forms a distinct niche involving 31 immune and stromal cell types and is enriched in CD4+ T cells, myeloid and follicular dendritic cells, while being depleted of plasma cells. Moreover, we used machine learning to nominate ligand-receptor pairs enriched in the HRS cell niche. Specifically, we identified IL13 as a candidate survival factor. In support of this hypothesis, recombinant IL13 augmented the proliferation of HRS cells in vitro. In addition, genome-wide CRISPR/Cas9 loss-of-function studies across more than 1,000 human cancer cell lines showed that IL4R and IL13RA1, the heterodimeric partners that constitute the IL13 receptor, were uniquely required for the survival of HRS cells. Moreover, monoclonal antibodies targeting either IL4R or IL13R phenocopied the genetic loss of function studies. IL13-targeting antibodies are already FDA-approved for atopic dermatitis, suggesting that clinical trials testing such agents should be explored in patients with Hodgkin lymphoma.

MetaLigand: A database for predicting non-peptide ligand mediated cell-cell communication

Non-peptide ligands (NPLs), including lipids, amino acids, carbohydrates, and non-peptide neurotransmitters and hormones, play a critical role in ligand-receptor-mediated cell-cell communication, driving diverse physiological and pathological processes. To facilitate the study of NPL-dependent intercellular interactions, we introduce MetaLigand, an R-based and web-accessible tool designed to infer NPL production and predict NPL-receptor interactions using transcriptomic data. MetaLigand compiles data for 233 NPLs, including their biosynthetic enzymes, transporter genes, and receptor genes, through a combination of automated pipelines and manual curation from comprehensive databases. The tool integrates both de novo and salvage synthesis pathways, incorporating multiple biosynthetic steps and transport mechanisms to improve prediction accuracy. Comparisons with existing tools demonstrate MetaLigand's superior ability to account for complex biogenesis pathways and model NPL abundance across diverse tissues and cell types. Furthermore, analysis of single-nucleus RNA-seq datasets from age-related macular degeneration samples revealed that distinct retinal cell types exhibit unique NPL profiles and participate in specific NPL-mediated pathological cell-cell interactions. Finally, MetaLigand supports single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics data, enabling the visualization of predicted NPL production levels and heterogeneity at single-cell resolution.

Intercellular Tunneling Nanotubes as Natural Biophotonic Conveyors

Tunneling nanotubes (TNTs), submicrometer membranous channels that bridge and connect distant cells, play a pivotal role in intercellular communication. Organelle transfer within TNTs is crucial in regulating cell growth, signal transmission, and disease progression. However, precise control over individual organelle transport within TNTs remains elusive. In this study, we introduce an optical technique that harnesses TNTs as biophotonic conveyors for the directional transport of individual organelles between cells. By utilizing near-infrared light propagating along the TNTs, optical forces were exerted on the organelles, enabling their active transport in a predetermined direction and at a controlled velocity. As a potential application, TNT conveyors were employed to inhibit mitochondrial hijacking from immune cells to cancer cells, thereby activating immune cells and suppressing cancer cell growth. Furthermore, neural modulation was achieved by transporting mitochondria and neurotransmitter-containing vesicles between neurons via TNT conveyors and axonal conveyors, respectively. This study presents a robust and precise approach to immune activation and neural regulation through the manipulation of organelle transfer at the subcellular level.

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.

Spatiotemporal modeling of molecular holograms

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

Plasticity of Human Microglia and Brain Perivascular Macrophages in Aging and Alzheimer's Disease

The complex roles of myeloid cells, including microglia and perivascular macrophages, are central to the neurobiology of Alzheimer's disease (AD), yet they remain incompletely understood. Here, we profiled 832,505 human myeloid cells from the prefrontal cortex of 1,607 unique donors covering the human lifespan and varying degrees of AD neuropathology. We delineated 13 transcriptionally distinct myeloid subtypes organized into 6 subclasses and identified AD-associated adaptive changes in myeloid cells over aging and disease progression. The GPNMB subtype, linked to phagocytosis, increased significantly with AD burden and correlated with polygenic AD risk scores. By organizing AD-risk genes into a regulatory hierarchy, we identified and validated MITF as an upstream transcriptional activator of GPNMB, critical for maintaining phagocytosis. Through cell-to-cell interaction networks, we prioritized APOE-SORL1 and APOE-TREM2 ligand-receptor pairs, associated with AD progression. In both human and mouse models, TREM2 deficiency disrupted GPNMB expansion and reduced phagocytic function, suggesting that GPNMB's role in neuroprotection was TREM2-dependent. Our findings clarify myeloid subtypes implicated in aging and AD, advancing the mechanistic understanding of their role in AD and aiding therapeutic discovery.

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.

Identification of a novel immune-related gene signature by single-cell and bulk sequencing for the prediction of the immune landscape and prognosis of breast cancer

BACKGROUND

As a common cause of cancer-related deaths in women, BRCA (breast cancer) shows complexity and requires precise biomarkers and treatment methods. This study delves into the molecular makeup of BRCA, focusing on immune profiles, molecular subtypes, gene expression and single-cell analysis.

METHODS

XCell was used to assess immune infiltration based on TCGA (the Cancer Genome Atlas) data and the clustering analysis was made. Differentially expressed genes were examined in distinct clusters, and the WGCNA (weighted correlation network analysis) was made to establish co-expression networks. The prognostic models were developed by Cox and LASSO-Cox regression. The clustering analysis, GSEA (Gene set enrichment analysis), GSVA (gene set variation analysis) and communication analysis of the single-cell dataset GSE161529 were performed to investigate the functional relevance. Real-time polymerase chain reaction (RT-PCR) was employed for evaluating gene expression.

RESULTS

The results revealed significant differences in immune cell infiltration between two clusters (C1 and C2). C2 had poorer survival outcomes, which was associated with higher expression of immune checkpoints PD1 and PD-L1. The gene modules identified via WGCNA were correlated with the immune-based subtypes. Then, a prognostic model comprising seven genes (ACSL1, ABCB5, XG, ADH4, OPN4, NPR3, NLGN1) was used to divide patients into high- and low-risk subgroups. The high-risk group had worse prognosis and higher scores of TIDE (Tumor Immune Dysfunction and Exclusion). The single-cell analysis depicted the immune landscape. Macrophages and endothelial cells exhibited higher AUCell scores. In cellular communication analysis, notably significant ligand-receptor interactions of HLA-DRA-> CD4 and TNFSF13B-> HLA-DPB1 were observed. The proportion of endothelial cells was correlated with risk scores. Finally, RT-PCR results illustrated the expression of seven genes in BRCA specimens.

CONCLUSION

The integrative analysis provides new insights into molecular complexities of BRCA. Immune profiles and gene signatures hold potential for improving stratification of BRCA patients and guiding the development of personalized immunotherapy strategies.

Single-cell transcriptomics and tissue metabolomics uncover mechanisms underlying wooden breast disease in broilers

Accompanied by the accelerated growth rate of chickens, the quality of chicken meat has deteriorated in recent years. Wooden breast (WB) is a severe myopathy affecting meat quality, and its pathophysiology depends on gene expression and intercellular interactions of various cell types, which are not yet fully understood. We have performed a comprehensive transcriptomic and metabolomic atlas of chicken WB muscle. Our data showed a significant increase in the number of immune cells, WB muscle displayed a unique cluster of macrophages (cluster 11), distinct from the M1 and M2 macrophages. Regarding the myocytes, the most significant differences were the decrease in cell number and the intensification of fatty deposits. Satellite cells were involved in muscle repair and regeneration producing more collagen. Interestingly, the interaction network in the WB group was weaker compared to that in normal breast muscle. Additionally, we found six key differential metabolites across 22 pathways. When WB occurs, myocytes and endothelial cells undergo apoptosis, macrophages are activated and exert immune functions, satellite cells participate in muscle rebuilding and repair, and the content of metabolites undergoes significant changes. This cell transcriptome profile provides an essential reference for future studies on the development and remodeling of WB.

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.

scDCA: deciphering the dominant cell communication assembly of downstream functional events from single-cell RNA-seq data

Cell-cell communications (CCCs) involve signaling from multiple sender cells that collectively impact downstream functional processes in receiver cells. Currently, computational methods are lacking for quantifying the contribution of pairwise combinations of cell types to specific functional processes in receiver cells (e.g. target gene expression or cell states). This limitation has impeded understanding the underlying mechanisms of cancer progression and identifying potential therapeutic targets. Here, we proposed a deep learning-based method, scDCA, to decipher the dominant cell communication assembly (DCA) that have a higher impact on a particular functional event in receiver cells from single-cell RNA-seq data. Specifically, scDCA employed a multi-view graph convolution network to reconstruct the CCCs landscape at single-cell resolution, and then identified DCA by interpreting the model with the attention mechanism. Taking the samples from advanced renal cell carcinoma as a case study, the scDCA was successfully applied and validated in revealing the DCA affecting the crucial gene expression in immune cells. The scDCA was also applied and validated in revealing the DCA responsible for the variation of 14 typical functional states of malignant cells. Furthermore, the scDCA was applied and validated to explore the alteration of CCCs under clinical intervention by comparing the DCA for certain cytotoxic factors between patients with and without immunotherapy. In summary, scDCA provides a valuable and practical tool for deciphering the cell type combinations with the most dominant impact on a specific functional process of receiver cells, which is of great significance for precise cancer treatment. Our data and code are free available at a public GitHub repository: https://github.com/pengsl-lab/scDCA.git.

Lung-resident alveolar macrophages regulate the timing of breast cancer metastasis

Breast disseminated cancer cells (DCCs) can remain dormant in the lungs for extended periods, but the mechanisms limiting their expansion are not well understood. Research indicates that tissue-resident alveolar macrophages suppress breast cancer metastasis in lung alveoli by inducing dormancy. Through ligand-receptor mapping and intravital imaging, it was found that alveolar macrophages express transforming growth factor (TGF)-β2. This expression, along with persistent macrophage-cancer cell interactions via the TGF-βRIII receptor, maintains cancer cells in a dormant state. Depleting alveolar macrophages or losing the TGF-β2 receptor in cancer cells triggers metastatic awakening. Aggressive breast cancer cells are either suppressed by alveolar macrophages or evade this suppression by avoiding interaction and downregulating the TGF-β2 receptor. Restoring TGF-βRIII in aggressive cells reinstates TGF-β2-mediated macrophage growth suppression. Thus, alveolar macrophages act as a metastasis immune barrier, and downregulation of TGF-β2 signaling allows cancer cells to overcome macrophage-mediated growth suppression.

Beyond genomic studies of congenital heart defects through systematic modelling and phenotyping

Congenital heart defects (CHDs), the most common congenital anomalies, are considered to have a significant genetic component. However, despite considerable efforts to identify pathogenic genes in patients with CHDs, few gene variants have been proven as causal. The complexity of the genetic architecture underlying human CHDs likely contributes to this poor genetic discovery rate. However, several other factors are likely to contribute. For example, the level of patient phenotyping required for clinical care may be insufficient for research studies focused on mechanistic discovery. Although several hundred mouse gene knockouts have been described with CHDs, these are generally not phenotyped and described in the same way as CHDs in patients, and thus are not readily comparable. Moreover, most patients with CHDs carry variants of uncertain significance of crucial cardiac genes, further complicating comparisons between humans and mouse mutants. In spite of major advances in cardiac developmental biology over the past 25 years, these advances have not been well communicated to geneticists and cardiologists. As a consequence, the latest data from developmental biology are not always used in the design and interpretation of studies aimed at discovering the genetic causes of CHDs. In this Special Article, while considering other in vitro and in vivo models, we create a coherent framework for accurately modelling and phenotyping human CHDs in mice, thereby enhancing the translation of genetic and genomic studies into the causes of CHDs in patients.

Characterization of ligand-receptor pair in acute myeloid leukemia: a scoring model for prognosis, therapeutic response, and T cell dysfunction

Introduction

The significance of ligand-receptor (LR) pair interactions in the progression of acute myeloid leukemia (AML) has been the focus of numerous studies. However, the relationship between LR pairs and the prognosis of AML, as well as their impact on treatment outcomes, is not fully elucidated.

Methods

Leveraging data from the TCGA-LAML cohort, we mapped out the LR pair interactions and distinguished specific molecular subtypes, with each displaying distinct biological characteristics. These subtypes exhibited varying mutation landscapes, pathway characteristics, and immune infiltration levels. Further insight into the immune microenvironment among these subtypes revealed disparities in immune cell abundance.

Results

Notably, one subtype showed a higher prevalence of CD8 T cells and plasma cells, suggesting increased adaptive immune activities. Leveraging a multivariate Lasso regression, we formulated an LR pair-based scoring model, termed "LR.score," to classify patients based on prognostic risk. Our findings underscored the association between elevated LR scores and T-cell dysfunction in AML. This connection highlights the LR score's potential as both a prognostic marker and a guide for personalized therapeutic interventions. Moreover, our LR.score revealed substantial survival prediction capacities in an independent AML cohort. We highlighted CLEC11A, ICAM4, ITGA4, and AVP as notably AML-specific.

Discussion

qRT-PCR analysis on AML versus normal bone marrow samples confirmed the significant downregulation of CLEC11A, ITGA4, ICAM4, and AVP in AML, suggesting their inverse biomarker potential in AML. In summary, this study illuminates the significance of the LR pair network in predicting AML prognosis, offering avenues for more precise treatment strategies tailored to individual patient profiles.

Mapping cellular interactions from spatially resolved transcriptomics data

Cell-cell communication (CCC) is essential to how life forms and functions. However, accurate, high-throughput mapping of how expression of all genes in one cell affects expression of all genes in another cell is made possible only recently through the introduction of spatially resolved transcriptomics (SRT) technologies, especially those that achieve single-cell resolution. Nevertheless, substantial challenges remain to analyze such highly complex data properly. Here, we introduce a multiple-instance learning framework, Spacia, to detect CCCs from data generated by SRTs, by uniquely exploiting their spatial modality. We highlight Spacia's power to overcome fundamental limitations of popular analytical tools for inference of CCCs, including losing single-cell resolution, limited to ligand-receptor relationships and prior interaction databases, high false positive rates and, most importantly, the lack of consideration of the multiple-sender-to-one-receiver paradigm. We evaluated the fitness of Spacia for three commercialized single-cell resolution SRT technologies: MERSCOPE/Vizgen, CosMx/NanoString and Xenium/10x. Overall, Spacia represents a notable step in advancing quantitative theories of cellular communications.

Studying the Pulmonary Endothelium in Health and Disease: An Official American Thoracic Society Workshop Report

Lung endothelium resides at the interface between the circulation and the underlying tissue, where it senses biochemical and mechanical properties of both the blood as it flows through the vascular circuit and the vessel wall. The endothelium performs the bidirectional signaling between the blood and tissue compartments that is necessary to maintain homeostasis while physically separating both, facilitating a tightly regulated exchange of water, solutes, cells, and signals. Disruption in endothelial function contributes to vascular disease, which can manifest in discrete vascular locations along the artery-to-capillary-to-vein axis. Although our understanding of mechanisms that contribute to endothelial cell injury and repair in acute and chronic vascular disease have advanced, pathophysiological mechanisms that underlie site-specific vascular disease remain incompletely understood. In an effort to improve the translatability of mechanistic studies of the endothelium, the American Thoracic Society convened a workshop to optimize rigor, reproducibility, and translation of discovery to advance our understanding of endothelial cell function in health and disease.

Integrated multiomics revealed adenosine signaling predict immunotherapy response and regulate tumor ecosystem of melanoma

Extracellular adenosine is extensively involved in regulating the tumor microenvironment. Given the disappointing results of adenosine-targeted therapy trials, personalized treatment might be necessary, tailored to the microenvironment status of individual patients. Here, we introduce the adenosine signaling score (ADO-score) model using non-negative matrix fraction identified patient subtypes using publicly available melanoma dataset, which aimed to profile adenosine signaling-related genes and construct a model to predict prognosis. We analyzed 580 malignant melanoma samples and demonstrated its robust value for prognosis. Further investigation in immune checkpoint inhibitor dataset suggests its potential as a stratified factor of immune checkpoint inhibitor efficacy. We validated the power of the ADO-score at the protein level immunofluorescence in a melanoma cohort from Xiangya Hospital. More importantly, single-cell and spatial transcriptomic data highlighted the cell-specific expression patterns of adenosine signaling-related genes and the existence of adenosine signaling-mediated crosstalk between tumor cells and immune cells in melanoma. Our study reveals a robust connection between adenosine signaling and clinical benefits in melanoma patients and proposes a universally applicable adenosine signaling model, the ADO-score, in gene expression profiles and histological sections. This model enables us to more precisely and conveniently select patients who are likely to benefit from immunotherapy.

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.

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.

Dynamic intrauterine crosstalk promotes porcine embryo implantation during early pregnancy

Dynamic crosstalk between the embryo and mother is crucial during implantation. Here, we comprehensively profile the single-cell transcriptome of pig peri-implantation embryos and corresponding maternal endometrium, identifying 4 different lineages in embryos and 13 cell types in the endometrium. Cell-specific gene expression characterizes 4 distinct trophectoderm subpopulations, showing development from undifferentiated trophectoderm to polar and mural trophectoderm. Dynamic expression of genes in different types of endometrial cells illustrates their molecular response to embryos during implantation. Then, we developed a novel tool, ExtraCellTalk, generating an overall dynamic map of maternal-foetal crosstalk using uterine luminal proteins as bridges. Through cross-species comparisons, we identified a conserved RBP4/STRA6 pathway in which embryonic-derived RBP4 could target the STRA6 receptor on stromal cells to regulate the interaction with other endometrial cells. These results provide insight into the maternal-foetal crosstalk during embryo implantation and represent a valuable resource for further studies to improve embryo implantation.

Single-cell multiregion dissection of Alzheimer's disease

Alzheimer's disease is the leading cause of dementia worldwide, but the cellular pathways that underlie its pathological progression across brain regions remain poorly understood1-3. Here we report a single-cell transcriptomic atlas of six different brain regions in the aged human brain, covering 1.3 million cells from 283 post-mortem human brain samples across 48 individuals with and without Alzheimer's disease. We identify 76 cell types, including region-specific subtypes of astrocytes and excitatory neurons and an inhibitory interneuron population unique to the thalamus and distinct from canonical inhibitory subclasses. We identify vulnerable populations of excitatory and inhibitory neurons that are depleted in specific brain regions in Alzheimer's disease, and provide evidence that the Reelin signalling pathway is involved in modulating the vulnerability of these neurons. We develop a scalable method for discovering gene modules, which we use to identify cell-type-specific and region-specific modules that are altered in Alzheimer's disease and to annotate transcriptomic differences associated with diverse pathological variables. We identify an astrocyte program that is associated with cognitive resilience to Alzheimer's disease pathology, tying choline metabolism and polyamine biosynthesis in astrocytes to preserved cognitive function late in life. Together, our study develops a regional atlas of the ageing human brain and provides insights into cellular vulnerability, response and resilience to Alzheimer's disease pathology.

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.

Crosstalk among proximal tubular cells, macrophages, and fibroblasts in acute kidney injury: single-cell profiling from the perspective of ferroptosis

The link between ferroptosis, a form of cell death mediated by iron and acute kidney injury (AKI) is recently gaining widespread attention. However, the mechanism of the crosstalk between cells in the pathogenesis and progression of acute kidney injury remains unexplored. In our research, we performed a non-negative matrix decomposition (NMF) algorithm on acute kidney injury single-cell RNA sequencing data based specifically focusing in ferroptosis-associated genes. Through a combination with pseudo-time analysis, cell-cell interaction analysis and SCENIC analysis, we discovered that proximal tubular cells, macrophages, and fibroblasts all showed associations with ferroptosis in different pathways and at various time. This involvement influenced cellular functions, enhancing cellular communication and activating multiple transcription factors. In addition, analyzing bulk expression profiles and marker genes of newly defined ferroptosis subtypes of cells, we have identified crucial cell subtypes, including Egr1 + PTC-C1, Jun + PTC-C3, Cxcl2 + Mac-C1 and Egr1 + Fib-C1. All these subtypes which were found in AKI mice kidneys and played significantly distinct roles from those of normal mice. Moreover, we verified the differential expression of Egr1, Jun, and Cxcl2 in the IRI mouse model and acute kidney injury human samples. Finally, our research presented a novel analysis of the crosstalk of proximal tubular cells, macrophages and fibroblasts in acute kidney injury targeting ferroptosis, therefore, contributing to better understanding the acute kidney injury pathogenesis, self-repairment and acute kidney injury-chronic kidney disease (AKI-CKD) progression.

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.

Dissecting tumor microenvironment from spatially resolved transcriptomics data by heterogeneous graph learning

Spatially resolved transcriptomics (SRT) has enabled precise dissection of tumor-microenvironment (TME) by analyzing its intracellular molecular networks and intercellular cell-cell communication (CCC). However, lacking computational exploration of complicated relations between cells, genes, and histological regions, severely limits the ability to interpret the complex structure of TME. Here, we introduce stKeep, a heterogeneous graph (HG) learning method that integrates multimodality and gene-gene interactions, in unraveling TME from SRT data. stKeep leverages HG to learn both cell-modules and gene-modules by incorporating features of diverse nodes including genes, cells, and histological regions, allows for identifying finer cell-states within TME and cell-state-specific gene-gene relations, respectively. Furthermore, stKeep employs HG to infer CCC for each cell, while ensuring that learned CCC patterns are comparable across different cell-states through contrastive learning. In various cancer samples, stKeep outperforms other tools in dissecting TME such as detecting bi-potent basal populations, neoplastic myoepithelial cells, and metastatic cells distributed within the tumor or leading-edge regions. Notably, stKeep identifies key transcription factors, ligands, and receptors relevant to disease progression, which are further validated by the functional and survival analysis of independent clinical data, thereby highlighting its clinical prognostic and immunotherapy 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.

Deciphering Müller cell heterogeneity signatures in diabetic retinopathy across species: an integrative single-cell analysis

Diabetic retinopathy (DR), a leading cause of visual impairment, demands a profound comprehension of its cellular mechanisms to formulate effective therapeutic strategies. Our study presentes a comprehensive single-cell analysis elucidating the intricate landscape of Müller cells within DR, emphasizing their nuanced involvement. Utilizing scRNA-seq data from both Sprague-Dawley rat models and human patients, we delineated distinct Müller cell clusters and their corresponding gene expression profiles. These findings were further validated through differential gene expression analysis utilizing human transcriptomic data. Notably, certain Müller cell clusters displayed upregulation of the Rho gene, implying a phagocytic response to damaged photoreceptors within the DR microenvironment. This phenomenon was consistently observed across species. Additionally, the co-expression patterns of RHO and PDE6G within Müller cell clusters provided compelling evidence supporting their potential role in maintaining retinal integrity during DR. Our results offer novel insights into the cellular dynamics of DR and underscore Müller cells as promising therapeutic targets for preserving vision in retinal disorders induced by diabetes.

Bioengineered Hydrogels Recapitulate Fibroblast Heterogeneity in Cancer

Recently mapped transcriptomic landscapes reveal the extent of heterogeneity in cancer-associated fibroblasts (CAFs) beyond previously established single-gene markers. Functional analyses of individual CAF subsets within the tumor microenvironment are critical to develop more accurate CAF-targeting therapeutic strategies. However, there is a lack of robust preclinical models that reflect this heterogeneity in vitro. In this study, single-cell RNA sequencing datasets acquired from head and neck squamous cell carcinoma tissues to predict microenvironmental and cellular features governing individual CAF subsets are leveraged. Some of these features are then incorporated into a tunable hyaluronan-based hydrogel system to culture patient-derived CAFs. Control over hydrogel degradability and integrin adhesiveness enabled derivation of the predominant myofibroblastic and inflammatory CAF subsets, as shown through changes in cell morphology and transcriptomic profiles. Last, using these hydrogel-cultured CAFs, microtubule dynamics are identified, but not actomyosin contractility, as a key mediator of CAF plasticity. The recapitulation of CAF heterogeneity in vitro using defined hydrogels presents unique opportunities for advancing the understanding of CAF biology and evaluation of CAF-targeting therapeutics.

Spatial and single-cell colocalisation analysis reveals MDK-mediated immunosuppressive environment with regulatory T cells in colorectal carcinogenesis

BACKGROUND

Cell-cell interaction factors that facilitate the progression of adenoma to sporadic colorectal cancer (CRC) remain unclear, thereby hindering patient survival.

METHODS

We performed spatial transcriptomics on five early CRC cases, which included adenoma and carcinoma, and one advanced CRC. To elucidate cell-cell interactions within the tumour microenvironment (TME), we investigated the colocalisation network at single-cell resolution using a deep generative model for colocalisation analysis, combined with a single-cell transcriptome, and assessed the clinical significance in CRC patients.

FINDINGS

CRC cells colocalised with regulatory T cells (Tregs) at the adenoma-carcinoma interface. At early-stage carcinogenesis, cell-cell interaction inference between colocalised adenoma and cancer epithelial cells and Tregs based on the spatial distribution of single cells highlighted midkine (MDK) as a prominent signalling molecule sent from tumour epithelial cells to Tregs. Interaction between MDK-high CRC cells and SPP1+ macrophages and stromal cells proved to be the mechanism underlying immunosuppression in the TME. Additionally, we identified syndecan4 (SDC4) as a receptor for MDK associated with Treg colocalisation. Finally, clinical analysis using CRC datasets indicated that increased MDK/SDC4 levels correlated with poor overall survival in CRC patients.

INTERPRETATION

MDK is involved in the immune tolerance shown by Tregs to tumour growth. MDK-mediated formation of the TME could be a potential target for early diagnosis and treatment of CRC.

FUNDING

Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Science Research; OITA Cancer Research Foundation; AMED under Grant Number; Japan Science and Technology Agency (JST); Takeda Science Foundation; The Princess Takamatsu Cancer Research Fund.