Computation and visualization of cell-cell signaling topologies in single-cell systems data using Connectome

Lipid-associated macrophages reshape BAT cell identity in obesity

Obesity and type 2 diabetes cause a loss in brown adipose tissue (BAT) activity, but the molecular mechanisms that drive BAT cell remodeling remain largely unexplored. Using a multilayered approach, we comprehensively mapped a reorganization in BAT cells. We uncovered a subset of macrophages as lipid-associated macrophages (LAMs), which were massively increased in genetic and dietary model of BAT expansion. LAMs participate in this scenario by capturing extracellular vesicles carrying damaged lipids and mitochondria released from metabolically stressed brown adipocytes. CD36 scavenger receptor drove LAM phenotype, and CD36-deficient LAMs were able to increase brown fat genes in adipocytes. LAMs released transforming growth factor β1 (TGF-β1), which promoted the loss of brown adipocyte identity through aldehyde dehydrogenase 1 family member A1 (Aldh1a1) induction. These findings unfold cell dynamic changes in BAT during obesity and identify LAMs as key responders to tissue metabolic stress and drivers of loss of brown adipocyte identity.

A Thy-1-negative immunofibroblast population emerges as a key determinant of fibrotic outcomes to biomaterials

Fibrosis-associated fibroblasts have been identified across various fibrotic disorders, but not in the context of biomaterials, fibrotic encapsulation, and the foreign body response. In other fibrotic disorders, a fibroblast subpopulation defined by Thy-1 loss is strongly correlated with fibrosis yet we do not know what promotes Thy-1 loss. We have previously shown that Thy-1 is an integrin regulator enabling normal fibroblast mechanosensing, and here, leveraging nonfibrotic microporous annealed particle (MAP) hydrogels versus classical fibrotic bulk hydrogels, we demonstrate that Thy1-/- mice mount a fibrotic response to MAP gels that includes inflammatory signaling. We found that a distinct and cryptic α-smooth muscle actin-positive Thy-1- fibroblast population emerges in response to interleuklin-1β (IL-1β) and tumor necrosis factor-α (TNFα). Furthermore, IL-1β/TNFα-induced Thy-1- fibroblasts consist of two distinct subpopulations that are strongly proinflammatory. These findings illustrate the emergence of a unique proinflammatory, profibrotic fibroblast subpopulation that is central to fibrotic encapsulation of biomaterials.

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.

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.

Multimodal analyses of immune cells during bone repair identify macrophages as a therapeutic target in musculoskeletal trauma

Musculoskeletal traumatic injuries (MTI) involve soft tissue lesions adjacent to a bone fracture leading to fibrous nonunion. The impact of MTI on the inflammatory response to fracture and on the immunomodulation of skeletal stem/progenitor cells (SSPCs) remains unknown. Here, we used single cell transcriptomic analyses to describe the immune cell dynamics after bone fracture and identified distinct macrophage subsets with successive pro-inflammatory, pro-repair and anti-inflammatory profiles. Concurrently, SSPCs transition via a pro- and anti-inflammatory fibrogenic phase of differentiation prior to osteochondrogenic differentiation. In a preclinical MTI mouse model, the injury response of immune cells and SSPCs is disrupted leading to a prolonged pro-inflammatory phase and delayed resolution of inflammation. Macrophage depletion improves bone regeneration in MTI demonstrating macrophage involvement in fibrous nonunion. Finally, pharmacological inhibition of macrophages using the CSF1R inhibitor Pexidartinib ameliorates healing. These findings reveal the coordinated immune response of macrophages and skeletal stem/progenitor cells as driver of bone healing and as a primary target for the treatment of trauma-associated fibrosis.

Summary

Hachemi et al. report the immune cell atlas of bone repair revealing macrophages as pro-fibrotic regulators and a therapeutic target for musculoskeletal regeneration. Genetic depletion or pharmacological inhibition of macrophages improves bone healing in musculoskeletal trauma.

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.

Combining LIANA and Tensor-cell2cell to decipher cell-cell communication across multiple samples

In recent years, data-driven inference of cell-cell communication has helped reveal coordinated biological processes across cell types. Here, we integrate two tools, LIANA and Tensor-cell2cell, which, when combined, can deploy multiple existing methods and resources to enable the robust and flexible identification of cell-cell communication programs across multiple samples. In this work, we show how the integration of our tools facilitates the choice of method to infer cell-cell communication and subsequently perform an unsupervised deconvolution to obtain and summarize biological insights. We explain how to perform the analysis step by step in both Python and R and provide online tutorials with detailed instructions available at https://ccc-protocols.readthedocs.io/. This workflow typically takes ∼1.5 h to complete from installation to downstream visualizations on a graphics processing unit-enabled computer for a dataset of ∼63,000 cells, 10 cell types, and 12 samples.

Craniofacial chondrogenesis in organoids from human stem cell-derived neural crest cells

Knowledge of cell signaling pathways that drive human neural crest differentiation into craniofacial chondrocytes is incomplete, yet essential for using stem cells to regenerate craniomaxillofacial structures. To accelerate translational progress, we developed a differentiation protocol that generated self-organizing craniofacial cartilage organoids from human embryonic stem cell-derived neural crest stem cells. Histological staining of cartilage organoids revealed tissue architecture and staining typical of elastic cartilage. Protein and post-translational modification (PTM) mass spectrometry and snRNA-seq data showed that chondrocyte organoids expressed robust levels of cartilage extracellular matrix (ECM) components: many collagens, aggrecan, perlecan, proteoglycans, and elastic fibers. We identified two populations of chondroprogenitor cells, mesenchyme cells and nascent chondrocytes, and the growth factors involved in paracrine signaling between them. We show that ECM components secreted by chondrocytes not only create a structurally resilient matrix that defines cartilage, but also play a pivotal autocrine cell signaling role in determining chondrocyte fate.

Defining and benchmarking open problems in single-cell analysis

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

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.

Comparative analysis of cell-cell communication at single-cell resolution

Inference of cell-cell communication from single-cell RNA sequencing data is a powerful technique to uncover intercellular communication pathways, yet existing methods perform this analysis at the level of the cell type or cluster, discarding single-cell-level information. Here we present Scriabin, a flexible and scalable framework for comparative analysis of cell-cell communication at single-cell resolution that is performed without cell aggregation or downsampling. We use multiple published atlas-scale datasets, genetic perturbation screens and direct experimental validation to show that Scriabin accurately recovers expected cell-cell communication edges and identifies communication networks that can be obscured by agglomerative methods. Additionally, we use spatial transcriptomic data to show that Scriabin can uncover spatial features of interaction from dissociated data alone. Finally, we demonstrate applications to longitudinal datasets to follow communication pathways operating between timepoints. Our approach represents a broadly applicable strategy to reveal the full structure of niche-phenotype relationships in health and disease.

Computational cell-cell interaction technologies drive mechanistic and biomarker discovery in the tumor microenvironment

Complex networks of cell-cell interactions (CCIs) within the tumor microenvironment (TME) play a crucial role in cancer persistence. These communication axes represent prime targets for therapeutic intervention, but our incomplete understanding of the cellular heterogeneity and interacting partners within the TME remains a stubborn barrier to complete drug responses. This review outlines recent advances in the study of CCIs that leverage single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) technologies that can clarify TME dynamics. We anticipate that these strategies will promote discovery of CCIs critical to the tumor-immune interface and will, by extension, expand the repertoire of druggable tumor biomarkers.

Single-cell RNA sequencing to understand host-virus interactions

Single-cell RNA sequencing (scRNA-seq) has allowed for the profiling of host and virus transcripts and host-virus interactions at single-cell resolution. This review summarizes the existing scRNA-seq technologies together with their strengths and weaknesses. The applications of scRNA-seq in various virological studies are discussed in depth, which broaden the understanding of the immune atlas, host-virus interactions, and immune repertoire. scRNA-seq can be widely used for virology in the near future to better understand the pathogenic mechanisms and discover more effective therapeutic strategies.

Targeting tumor-stromal interactions in triple-negative breast cancer using a human vascularized micro-tumor model

Triple-negative breast cancer (TNBC) is highly aggressive with limited available treatments. Stromal cells in the tumor microenvironment (TME) are crucial in TNBC progression; however, understanding the molecular basis of stromal cell activation and tumor-stromal crosstalk in TNBC is limited. To investigate therapeutic targets in the TNBC stromal niche, we used an advanced human in vitro microphysiological system called the vascularized micro-tumor (VMT). Using single-cell RNA sequencing, we revealed that normal breast tissue stromal cells activate neoplastic signaling pathways in the TNBC TME. By comparing interactions in VMTs with clinical data, we identified therapeutic targets at the tumor-stromal interface with potential clinical significance. Combining treatments targeting Tie2 signaling with paclitaxel resulted in vessel normalization and increased efficacy of paclitaxel in the TNBC VMT. Dual inhibition of HER3 and Akt also showed efficacy against TNBC. These data demonstrate the potential of inducing a favorable TME as a targeted therapeutic approach in TNBC.

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.

The mycosis fungoides cutaneous microenvironment shapes dysfunctional cell trafficking, antitumor immunity, matrix interactions, and angiogenesis

Malignant T lymphocyte proliferation in mycosis fungoides (MF) is largely restricted to the skin, implying that malignant cells are dependent on their specific cutaneous tumor microenvironment (TME), including interactions with non-malignant immune and stromal cells, cytokines, and other immunomodulatory factors. To explore these interactions, we performed a comprehensive transcriptome analysis of the TME in advanced-stage MF skin tumors by single-cell RNA sequencing. Our analysis identified cell-type compositions, cellular functions, and cell-to-cell interactions in the MF TME that were distinct from those from healthy skin and benign dermatoses. While patterns of gene expression were common among patient samples, high transcriptional diversity was also observed in immune and stromal cells, with dynamic interactions and crosstalk between these cells and malignant T lymphocytes. This heterogeneity mapped to processes such as cell trafficking, matrix interactions, angiogenesis, immune functions, and metabolism that affect cancer cell growth, migration, and invasion, as well as antitumor immunity. By comprehensively characterizing the transcriptomes of immune and stromal cells within the cutaneous microenvironment of individual MF tumors, we have identified patterns of dysfunction common to all tumors that represent a resource for identifying candidates with therapeutic potential as well as patient-specific heterogeneity that has important implications for personalized disease management.

Mapping the protein-protein interactome in the tumor immune microenvironment

The cell-to-cell communication primarily occurs through cell-surface and secreted proteins, which form a sophisticated network that coordinates systemic immune function. Uncovering these protein-protein interactions (PPIs) is indispensable for understanding the molecular mechanism and elucidating immune system aberrances under diseases. Traditional biological studies typically focus on a limited number of PPI pairs due to the relative low throughput of commonly used techniques. Encouragingly, classical methods have advanced, and many new systems tailored for large-scale protein-protein screening have been developed and successfully utilized. These high-throughput PPI investigation techniques have already made considerable achievements in mapping the immune cell interactome, enriching PPI databases and analysis tools, and discovering therapeutic targets for cancer and other diseases, which will definitely bring unprecedented insight into this field.

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.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Enteric glial hub cells coordinate intestinal motility

Enteric glia are the predominant cell type in the enteric nervous system yet their identities and roles in gastrointestinal function are not well classified. Using our optimized single nucleus RNA-sequencing method, we identified distinct molecular classes of enteric glia and defined their morphological and spatial diversity. Our findings revealed a functionally specialized biosensor subtype of enteric glia that we call "hub cells." Deletion of the mechanosensory ion channel PIEZO2 from adult enteric glial hub cells, but not other subtypes of enteric glia, led to defects in intestinal motility and gastric emptying in mice. These results provide insight into the multifaceted functions of different enteric glial cell subtypes in gut health and emphasize that therapies targeting enteric glia could advance the treatment of gastrointestinal diseases.

Combining LIANA and Tensor-cell2cell to decipher cell-cell communication across multiple samples

In recent years, data-driven inference of cell-cell communication has helped reveal coordinated biological processes across cell types. While multiple cell-cell communication tools exist, results are specific to the tool of choice, due to the diverse assumptions made across computational frameworks. Moreover, tools are often limited to analyzing single samples or to performing pairwise comparisons. As experimental design complexity and sample numbers continue to increase in single-cell datasets, so does the need for generalizable methods to decipher cell-cell communication in such scenarios. Here, we integrate two tools, LIANA and Tensor-cell2cell, which combined can deploy multiple existing methods and resources, to enable the robust and flexible identification of cell-cell communication programs across multiple samples. In this protocol, we show how the integration of our tools facilitates the choice of method to infer cell-cell communication and subsequently perform an unsupervised deconvolution to obtain and summarize biological insights. We explain how to perform the analysis step-by-step in both Python and R, and we provide online tutorials with detailed instructions available at https://ccc-protocols.readthedocs.io/. This protocol typically takes ~1.5h to complete from installation to downstream visualizations on a GPU-enabled computer, for a dataset of ~63k cells, 10 cell types, and 12 samples.

Rational engineering of lung alveolar epithelium

Engineered whole lungs may one day expand therapeutic options for patients with end-stage lung disease. However, the feasibility of ex vivo lung regeneration remains limited by the inability to recapitulate mature, functional alveolar epithelium. Here, we modulate multimodal components of the alveolar epithelial type 2 cell (AEC2) niche in decellularized lung scaffolds in order to guide AEC2 behavior for epithelial regeneration. First, endothelial cells coordinate with fibroblasts, in the presence of soluble growth and maturation factors, to promote alveolar scaffold population with surfactant-secreting AEC2s. Subsequent withdrawal of Wnt and FGF agonism synergizes with tidal-magnitude mechanical strain to induce the differentiation of AEC2s to squamous type 1 AECs (AEC1s) in cultured alveoli, in situ. These results outline a rational strategy to engineer an epithelium of AEC2s and AEC1s contained within epithelial-mesenchymal-endothelial alveolar-like units, and highlight the critical interplay amongst cellular, biochemical, and mechanical niche cues within the reconstituting alveolus.

scTenifoldXct: A semi-supervised method for predicting cell-cell interactions and mapping cellular communication graphs

We present scTenifoldXct, a semi-supervised computational tool for detecting ligand-receptor (LR)-mediated cell-cell interactions and mapping cellular communication graphs. Our method is based on manifold alignment, using LR pairs as inter-data correspondences to embed ligand and receptor genes expressed in interacting cells into a unified latent space. Neural networks are employed to minimize the distance between corresponding genes while preserving the structure of gene regression networks. We apply scTenifoldXct to real datasets for testing and demonstrate that our method detects interactions with high consistency compared with other methods. More importantly, scTenifoldXct uncovers weak but biologically relevant interactions overlooked by other methods. We also demonstrate how scTenifoldXct can be used to compare different samples, such as healthy vs. diseased and wild type vs. knockout, to identify differential interactions, thereby revealing functional implications associated with changes in cellular communication status.

SARS-CoV-2 leverages airway epithelial protective mechanism for viral infection

Despite much concerted effort to better understand severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral infection, relatively little is known about the dynamics of early viral entry and infection in the airway. Here we analyzed a single-cell RNA sequencing dataset of early SARS-CoV-2 infection in a humanized in vitro model, to elucidate key mechanisms by which the virus triggers a cell-systems-level response in the bronchial epithelium. We find that SARS-CoV-2 virus preferentially enters the tissue via ciliated cell precursors, giving rise to a population of infected mature ciliated cells, which signal to basal cells, inducing further rapid differentiation. This feedforward loop of infection is mitigated by further cell-cell communication, before interferon signaling begins at three days post-infection. These findings suggest hijacking by the virus of potentially beneficial tissue repair mechanisms, possibly exacerbating the outcome. This work both elucidates the interplay between barrier tissues and viral infections and may suggest alternative therapeutic approaches targeting non-immune response mechanisms.

spaCI: deciphering spatial cellular communications through adaptive graph model

Cell-cell communications are vital for biological signalling and play important roles in complex diseases. Recent advances in single-cell spatial transcriptomics (SCST) technologies allow examining the spatial cell communication landscapes and hold the promise for disentangling the complex ligand-receptor (L-R) interactions across cells. However, due to frequent dropout events and noisy signals in SCST data, it is challenging and lack of effective and tailored methods to accurately infer cellular communications. Herein, to decipher the cell-to-cell communications from SCST profiles, we propose a novel adaptive graph model with attention mechanisms named spaCI. spaCI incorporates both spatial locations and gene expression profiles of cells to identify the active L-R signalling axis across neighbouring cells. Through benchmarking with currently available methods, spaCI shows superior performance on both simulation data and real SCST datasets. Furthermore, spaCI is able to identify the upstream transcriptional factors mediating the active L-R interactions. For biological insights, we have applied spaCI to the seqFISH+ data of mouse cortex and the NanoString CosMx Spatial Molecular Imager (SMI) data of non-small cell lung cancer samples. spaCI reveals the hidden L-R interactions from the sparse seqFISH+ data, meanwhile identifies the inconspicuous L-R interactions including THBS1-ITGB1 between fibroblast and tumours in NanoString CosMx SMI data. spaCI further reveals that SMAD3 plays an important role in regulating the crosstalk between fibroblasts and tumours, which contributes to the prognosis of lung cancer patients. Collectively, spaCI addresses the challenges in interrogating SCST data for gaining insights into the underlying cellular communications, thus facilitates the discoveries of disease mechanisms, effective biomarkers and therapeutic targets.

Single-cell transcriptomic and spatial landscapes of the developing human pancreas

Current differentiation protocols have not been successful in reproducibly generating fully functional human beta cells in vitro, partly due to incomplete understanding of human pancreas development. Here, we present detailed transcriptomic analysis of the various cell types of the developing human pancreas, including their spatial gene patterns. We integrated single-cell RNA sequencing with spatial transcriptomics at multiple developmental time points and revealed distinct temporal-spatial gene cascades. Cell trajectory inference identified endocrine progenitor populations and branch-specific genes as the progenitors differentiate toward alpha or beta cells. Spatial differentiation trajectories indicated that Schwann cells are spatially co-located with endocrine progenitors, and cell-cell connectivity analysis predicted that they may interact via L1CAM-EPHB2 signaling. Our integrated approach enabled us to identify heterogeneity and multiple lineage dynamics within the mesenchyme, showing that it contributed to the exocrine acinar cell state. Finally, we have generated an interactive web resource for investigating human pancreas development for the research community.

Comprehensive visualization of cell-cell interactions in single-cell and spatial transcriptomics with NICHES

MOTIVATION

Recent years have seen the release of several toolsets that reveal cell-cell interactions from single-cell data. However, all existing approaches leverage mean celltype gene expression values, and do not preserve the single-cell fidelity of the original data. Here, we present NICHES (Niche Interactions and Communication Heterogeneity in Extracellular Signaling), a tool to explore extracellular signaling at the truly single-cell level.

RESULTS

NICHES allows embedding of ligand-receptor signal proxies to visualize heterogeneous signaling archetypes within cell clusters, between cell clusters and across experimental conditions. When applied to spatial transcriptomic data, NICHES can be used to reflect local cellular microenvironment. NICHES can operate with any list of ligand-receptor signaling mechanisms, is compatible with existing single-cell packages, and allows rapid, flexible analysis of cell-cell signaling at single-cell resolution.

AVAILABILITY AND IMPLEMENTATION

NICHES is an open-source software implemented in R under academic free license v3.0 and it is available at http://github.com/msraredon/NICHES. Use-case vignettes are available at https://msraredon.github.io/NICHES/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Clinical Translation of Engineered Pulmonary Vascular Models

Diseases in pulmonary vasculature remain a major cause of morbidity and mortality worldwide. Numerous pre-clinical animal models were developed to understand lung vasculature during diseases and development. However, these systems are typically limited in their ability to represent human pathophysiology for the study of disease and drug mechanisms. In recent years, a growing number of studies have focused on developing in vitro experimental platforms that mimic human tissues/organs. In this chapter, we discuss the key components involved in developing engineered pulmonary vascular modeling systems and provide perspectives on ways to improve the translational potential of existing models.

Breast cancer prevention by short-term inhibition of TGFβ signaling

Cancer prevention has a profound impact on cancer-associated mortality and morbidity. We previously identified TGFβ signaling as a candidate regulator of mammary epithelial cells associated with breast cancer risk. Here, we show that short-term TGFBR inhibitor (TGFBRi) treatment of peripubertal ACI inbred and Sprague Dawley outbred rats induces lasting changes and prevents estrogen- and carcinogen-induced mammary tumors, respectively. We identify TGFBRi-responsive cell populations by single cell RNA-sequencing, including a unique epithelial subpopulation designated secretory basal cells (SBCs) with progenitor features. We detect SBCs in normal human breast tissues and find them to be associated with breast cancer risk. Interactome analysis identifies SBCs as the most interactive cell population and the main source of insulin-IGF signaling. Accordingly, inhibition of TGFBR and IGF1R decrease proliferation of organoid cultures. Our results reveal a critical role for TGFβ in regulating mammary epithelial cells relevant to breast cancer and serve as a proof-of-principle cancer prevention strategy.

CITEdb: a manually curated database of cell-cell interactions in human

MOTIVATION

The interactions among various types of cells play critical roles in cell functions and the maintenance of the entire organism. While cell-cell interactions are traditionally revealed from experimental studies, recent developments in single-cell technologies combined with data mining methods have enabled computational prediction of cell-cell interactions, which have broadened our understanding of how cells work together, and have important implications in therapeutic interventions targeting cell-cell interactions for cancers and other diseases. Despite the importance, to our knowledge, there is no database for systematic documentation of high-quality cell-cell interactions at the cell type level, which hinders the development of computational approaches to identify cell-cell interactions.

RESULTS

We develop a publicly accessible database, CITEdb (Cell-cell InTEraction database, https://citedb.cn/), which not only facilitates interactive exploration of cell-cell interactions in specific physiological contexts (e.g. a disease or an organ) but also provides a benchmark dataset to interpret and evaluate computationally derived cell-cell interactions from different tools. CITEdb contains 728 pairs of cell-cell interactions in human that are manually curated. Each interaction is equipped with structured annotations including the physiological context, the ligand-receptor pairs that mediate the interaction, etc. Our database provides a web interface to search, visualize and download cell-cell interactions. Users can search for cell-cell interactions by selecting the physiological context of interest or specific cell types involved. CITEdb is the first attempt to catalogue cell-cell interactions at the cell type level, which is beneficial to both experimental, computational and clinical studies of cell-cell interactions.

AVAILABILITY AND IMPLEMENTATION

CITEdb is freely available at https://citedb.cn/ and the R package implementing benchmark is available at https://github.com/shanny01/benchmark.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluation of cell-cell interaction methods by integrating single-cell RNA sequencing data with spatial information

Background

Cell-cell interactions are important for information exchange between different cells, which are the fundamental basis of many biological processes. Recent advances in single-cell RNA sequencing (scRNA-seq) enable the characterization of cell-cell interactions using computational methods. However, it is hard to evaluate these methods since no ground truth is provided. Spatial transcriptomics (ST) data profiles the relative position of different cells. We propose that the spatial distance suggests the interaction tendency of different cell types, thus could be used for evaluating cell-cell interaction tools.

Results

We benchmark 16 cell-cell interaction methods by integrating scRNA-seq with ST data. We characterize cell-cell interactions into short-range and long-range interactions using spatial distance distributions between ligands and receptors. Based on this classification, we define the distance enrichment score and apply an evaluation workflow to 16 cell-cell interaction tools using 15 simulated and 5 real scRNA-seq and ST datasets. We also compare the consistency of the results from single tools with the commonly identified interactions. Our results suggest that the interactions predicted by different tools are highly dynamic, and the statistical-based methods show overall better performance than network-based methods and ST-based methods.

Conclusions

Our study presents a comprehensive evaluation of cell-cell interaction tools for scRNA-seq. CellChat, CellPhoneDB, NicheNet, and ICELLNET show overall better performance than other tools in terms of consistency with spatial tendency and software scalability. We recommend using results from at least two methods to ensure the accuracy of identified interactions. We have packaged the benchmark workflow with detailed documentation at GitHub (https://github.com/wanglabtongji/CCI).

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02783-y.

CommPath: An R package for inference and analysis of pathway-mediated cell-cell communication chain from single-cell transcriptomics

Single-cell transcriptomics offers opportunities to investigate ligand-receptor (LR) interactions between heterogeneous cell populations within tissues. However, most existing tools for the inference of intercellular communication do not allow prioritization of functional LR associations that provoke certain biological responses in the receiver cells. In addition, current tools do not enable the identification of the impact on the downstream cell types of the receiver cells. We present CommPath, an open-source R package and webserver, to analyze and visualize the LR interactions and pathway-mediated intercellular communication chain with single-cell transcriptomic data. CommPath curates a comprehensive signaling pathway database to interpret the consequences of LR associations and therefore infers functional LR interactions. Furthermore, CommPath determines cell-cell communication chain by considering both the upstream and downstream cells of user-defined cell populations. Applying CommPath to human hepatocellular carcinoma dataset shows its ability to decipher complex LR interaction patterns and the associated intercellular communication chain, as well as their changes in disease versus homeostasis.

LRLoop: a method to predict feedback loops in cell-cell communication

MOTIVATION

Intercellular communication (i.e. cell-cell communication) plays an essential role in multicellular organisms coordinating various biological processes. Previous studies discovered that feedback loops between two cell types are a widespread and vital signaling motif regulating development, regeneration and cancer progression. While many computational methods have been developed to predict cell-cell communication based on gene expression datasets, these methods often predict one-directional ligand-receptor interactions from sender to receiver cells and are not suitable to identify feedback loops.

RESULTS

Here, we describe ligand-receptor loop (LRLoop), a new method for analyzing cell-cell communication based on bi-directional ligand-receptor interactions, where two pairs of ligand-receptor interactions are identified that are responsive to each other and thereby form a closed feedback loop. We first assessed LRLoop using bulk datasets and found our method significantly reduces the false positive rate seen with existing methods. Furthermore, we developed a new strategy to assess the performance of these methods in single-cell datasets. We used the between-tissue interactions as an indicator of potential false-positive prediction and found that LRLoop produced a lower fraction of between-tissue interactions than traditional methods. Finally, we applied LRLoop to the single-cell datasets obtained from retinal development. We discovered many new bi-directional ligand-receptor interactions among individual cell types that potentially control proliferation, neurogenesis and/or cell fate specification.

AVAILABILITY AND IMPLEMENTATION

An R package is available at https://github.com/Pinlyu3/LRLoop. The source code can be found at figshare (https://doi.org/10.6084/m9.figshare.20126138.v1). The datasets can be found at figshare (https://doi.org/10.6084/m9.figshare.20126021.v1).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Context-aware deconvolution of cell-cell communication with Tensor-cell2cell

Cell interactions determine phenotypes, and intercellular communication is shaped by cellular contexts such as disease state, organismal life stage, and tissue microenvironment. Single-cell technologies measure the molecules mediating cell–cell communication, and emerging computational tools can exploit these data to decipher intercellular communication. However, current methods either disregard cellular context or rely on simple pairwise comparisons between samples, thus limiting the ability to decipher complex cell–cell communication across multiple time points, levels of disease severity, or spatial contexts. Here we present Tensor-cell2cell, an unsupervised method using tensor decomposition, which deciphers context-driven intercellular communication by simultaneously accounting for multiple stages, states, or locations of the cells. To do so, Tensor-cell2cell uncovers context-driven patterns of communication associated with different phenotypic states and determined by unique combinations of cell types and ligand-receptor pairs. As such, Tensor-cell2cell robustly improves upon and extends the analytical capabilities of existing tools. We show Tensor-cell2cell can identify multiple modules associated with distinct communication processes (e.g., participating cell–cell and ligand-receptor pairs) linked to severities of Coronavirus Disease 2019 and to Autism Spectrum Disorder. Thus, we introduce an effective and easy-to-use strategy for understanding complex communication patterns across diverse conditions.

Cellular contexts such as disease state, organismal life stage and tissue microenvironment, shape intercellular communication, and ultimately affect an organism’s phenotypes. Here, the authors present Tensor-cell2cell, an unsupervised method for deciphering context-driven intercellular communication.