Reconstruction of cell spatial organization from single-cell RNA sequencing data based on ligand-receptor mediated self-assembly

Simulating multiple variability in spatially resolved transcriptomics with scCube

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

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

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

Single-cell landscape reveals the epithelial cell-centric pro-inflammatory immune microenvironment in dry eye development

Dry eye disease (DED) is a prevalent chronic eye disease characterized by an aberrant inflammatory response in ocular surface mucosa. The immunological alterations underlying DED remain largely unknown. In this study, we employed single-cell transcriptome sequencing of conjunctival tissue from environment-induced DED mice to investigate multicellular ecosystem and functional changes at different DED stages. Our results revealed an epithelial subtype with fibroblastic characteristics and pro-inflammatory effects emerging in the acute phase of DED. We also found that T helper (Th)1, Th17, and regulatory T cells (Treg) were the dominant clusters of differentiation (CD)4+ T-cell types involved in regulating immune responses and identified three distinct macrophage subtypes, with the CD72+CD11c+ subtype enhancing chronic inflammation. Furthermore, bulk transcriptome analysis of video display terminal-induced DED consistently suggested the presence of the pro-inflammatory epithelial subtype in human conjunctiva. Our findings have uncovered a DED-associated pro-inflammatory microenvironment in the conjunctiva, centered around epithelial cells, involving interactions with macrophages and CD4+ T cells, which deepens our understanding of ocular surface mucosal immune responses during DED progression.

Uncover spatially informed variations for single-cell spatial transcriptomics with STew

Motivation

The recent spatial transcriptomics (ST) technologies have enabled characterization of gene expression patterns and spatial information, advancing our understanding of cell lineages within diseased tissues. Several analytical approaches have been proposed for ST data, but effectively utilizing spatial information to unveil the shared variation with gene expression remains a challenge.

Results

We introduce STew, a Spatial Transcriptomic multi-viEW representation learning method, to jointly analyze spatial information and gene expression in a scalable manner, followed by a data-driven statistical framework to measure the goodness of model fit. Through benchmarking using human dorsolateral prefrontal cortex and mouse main olfactory bulb data with true manual annotations, STew achieved superior performance in both clustering accuracy and continuity of identified spatial domains compared with other methods. STew is also robust to generate consistent results insensitive to model parameters, including sparsity constraints. We next applied STew to various ST data acquired from 10× Visium, Slide-seqV2, and 10× Xenium, encompassing single-cell and multi-cellular resolution ST technologies, which revealed spatially informed cell type clusters and biologically meaningful axes. In particular, we identified a proinflammatory fibroblast spatial niche using ST data from psoriatic skins. Moreover, STew scales almost linearly with the number of spatial locations, guaranteeing its applicability to datasets with thousands of spatial locations to capture disease-relevant niches in complex tissues.

Availability and implementation

Source code and the R software tool STew are available from github.com/fanzhanglab/STew.

Single-cell analysis of anti-BCMA CAR T cell therapy in patients with central nervous system autoimmunity

Chimeric antigen receptor (CAR) T cell immunotherapy for the treatment of neurological autoimmune diseases is promising, but CAR T cell kinetics and immune alterations after treatment are poorly understood. Here, we performed single-cell multi-omics sequencing of paired cerebrospinal fluid (CSF) and blood samples from patients with neuromyelitis optica spectrum disorder (NMOSD) treated with anti-B cell maturation antigen (BCMA) CAR T cells. Proliferating cytotoxic-like CD8+ CAR T cell clones were identified as the main effectors in autoimmunity. Anti-BCMA CAR T cells with enhanced features of chemotaxis efficiently crossed the blood-CSF barrier, eliminated plasmablasts and plasma cells in the CSF, and suppressed neuroinflammation. The CD44-expressing early memory phenotype in infusion products was potentially associated with CAR T cell persistence in autoimmunity. Moreover, CAR T cells from patients with NMOSD displayed distinctive features of suppressed cytotoxicity compared with those from hematological malignancies. Thus, we provide mechanistic insights into CAR T cell function in patients with neurological autoimmune disease.

Single-cell senescence identification reveals senescence heterogeneity, trajectory, and modulators

Cellular senescence underlies many aging-related pathologies, but its heterogeneity poses challenges for studying and targeting senescent cells. We present here a machine learning program senescent cell identification (SenCID), which accurately identifies senescent cells in both bulk and single-cell transcriptome. Trained on 602 samples from 52 senescence transcriptome datasets spanning 30 cell types, SenCID identifies six major senescence identities (SIDs). Different SIDs exhibit different senescence baselines, stemness, gene functions, and responses to senolytics. SenCID enables the reconstruction of senescent trajectories under normal aging, chronic diseases, and COVID-19. Additionally, when applied to single-cell Perturb-seq data, SenCID helps reveal a hierarchy of senescence modulators. Overall, SenCID is an essential tool for precise single-cell analysis of cellular senescence, enabling targeted interventions against senescent cells.

Single-cell transcriptomics: background, technologies, applications, and challenges

Single-cell sequencing was developed as a high-throughput tool to elucidate unusual and transient cell states that are barely visible in the bulk. This technology reveals the evolutionary status of cells and differences between populations, helps to identify unique cell subtypes and states, reveals regulatory relationships between genes, targets and molecular mechanisms in disease processes, tumor heterogeneity, the state of the immune environment, etc. However, the high cost and technical limitations of single-cell sequencing initially prevented its widespread application, but with advances in research, numerous new single-cell sequencing techniques have been discovered, lowering the cost barrier. Many single-cell sequencing platforms and bioinformatics methods have recently become commercially available, allowing researchers to make fascinating observations. They are now increasingly being used in various industries. Several protocols have been discovered in this context and each technique has unique characteristics, capabilities and challenges. This review presents the latest advancements in single-cell transcriptomics technologies. This includes single-cell transcriptomics approaches, workflows and statistical approaches to data processing, as well as the potential advances, applications, opportunities and challenges of single-cell transcriptomics technology. You will also get an overview of the entry points for spatial transcriptomics and multi-omics.

Spatial imaging of glycoRNA in single cells with ARPLA

Little is known about the biological roles of glycosylated RNAs (glycoRNAs), a recently discovered class of glycosylated molecules, because of a lack of visualization methods. We report sialic acid aptamer and RNA in situ hybridization-mediated proximity ligation assay (ARPLA) to visualize glycoRNAs in single cells with high sensitivity and selectivity. The signal output of ARPLA occurs only when dual recognition of a glycan and an RNA triggers in situ ligation, followed by rolling circle amplification of a complementary DNA, which generates a fluorescent signal by binding fluorophore-labeled oligonucleotides. Using ARPLA, we detect spatial distributions of glycoRNAs on the cell surface and their colocalization with lipid rafts as well as the intracellular trafficking of glycoRNAs through SNARE protein-mediated secretory exocytosis. Studies in breast cell lines suggest that surface glycoRNA is inversely associated with tumor malignancy and metastasis. Investigation of the relationship between glycoRNAs and monocyte-endothelial cell interactions suggests that glycoRNAs may mediate cell-cell interactions during the immune response.

A division-of-labor mode contributes to the cardioprotective potential of mesenchymal stem/stromal cells in heart failure post myocardial infarction

Background

Treatment of heart failure post myocardial infarction (post-MI HF) with mesenchymal stem/stromal cells (MSCs) holds great promise. Nevertheless, 2-dimensional (2D) GMP-grade MSCs from different labs and donor sources have different therapeutic efficacy and still in a low yield. Therefore, it is crucial to increase the production and find novel ways to assess the therapeutic efficacy of MSCs.

Materials and methods

hUC-MSCs were cultured in 3-dimensional (3D) expansion system for obtaining enough cells for clinical use, named as 3D MSCs. A post-MI HF mouse model was employed to conduct in vivo and in vitro experiments. Single-cell and bulk RNA-seq analyses were performed on 3D MSCs. A total of 125 combination algorithms were leveraged to screen for core ligand genes. Shinyapp and shinycell workflows were used for deploying web-server.

Result

3D GMP-grade MSCs can significantly and stably reduce the extent of post-MI HF. To understand the stable potential cardioprotective mechanism, scRNA-seq revealed the heterogeneity and division-of-labor mode of 3D MSCs at the cellular level. Specifically, scissor phenotypic analysis identified a reported wound-healing CD142+ MSCs subpopulation that is also associated with cardiac protection ability and CD142- MSCs that is in proliferative state, contributing to the cardioprotective function and self-renewal, respectively. Differential expression analysis was conducted on CD142+ MSCs and CD142- MSCs and the differentially expressed ligand-related model was achieved by employing 125 combination algorithms. The present study developed a machine learning predictive model based on 13 ligands. Further analysis using CellChat demonstrated that CD142+ MSCs have a stronger secretion capacity compared to CD142- MSCs and Flow cytometry sorting of the CD142+ MSCs and qRT-PCR validation confirmed the significant upregulation of these 13 ligand factors in CD142+ MSCs.

Conclusion

Clinical GMP-grade 3D MSCs could serve as a stable cardioprotective cell product. Using scissor analysis on scRNA-seq data, we have clarified the potential functional and proliferative subpopulation, which cooperatively contributed to self-renewal and functional maintenance for 3D MSCs, named as "division of labor" mode of MSCs. Moreover, a ligand model was robustly developed for predicting the secretory efficacy of MSCs. A user-friendly web-server and a predictive model were constructed and available (https://wangxc.shinyapps.io/3D_MSCs/).

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

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

Resource allocation in mammalian systems

Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.

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.

The Rise of Hypothesis-Driven Artificial Intelligence in Oncology

Cancer is a complex disease involving the deregulation of intricate cellular systems beyond genetic aberrations and, as such, requires sophisticated computational approaches and high-dimensional data for optimal interpretation. While conventional artificial intelligence (AI) models excel in many prediction tasks, they often lack interpretability and are blind to the scientific hypotheses generated by researchers to enable cancer discoveries. Here we propose that hypothesis-driven AI, a new emerging class of AI algorithm, is an innovative approach to uncovering the complex etiology of cancer from big omics data. This review exemplifies how hypothesis-driven AI is different from conventional AI by citing its application in various areas of oncology including tumor classification, patient stratification, cancer gene discovery, drug response prediction, and tumor spatial organization. Our aim is to stress the feasibility of incorporating domain knowledge and scientific hypotheses to craft the design of new AI algorithms. We showcase the power of hypothesis-driven AI in making novel cancer discoveries that can be overlooked by conventional AI methods. Since hypothesis-driven AI is still in its infancy, open questions such as how to better incorporate new knowledge and biological perspectives to ameliorate bias and improve interpretability in the design of AI algorithms still need to be addressed. In conclusion, hypothesis-driven AI holds great promise in the discovery of new mechanistic and functional insights that explain the complexity of cancer etiology and potentially chart a new roadmap to improve treatment regimens for individual patients.

Creation of a High-Throughput Microfluidic Platform for Single-Cell Transcriptome Sequencing of Cell-Cell Interactions

Cell-cell interaction is one of the major modalities for transmitting information between cells and activating the effects of functional cells. However, the construction of high-throughput analysis technologies from cell omics focusing on the impact of interactions of functional cells on targets has been relatively unexplored. Here, they propose a droplet-based microfluidic platform for cell-cell interaction sequencing (c-c-seq) and screening in vitro to address this challenge. A class of interacting cells is pre-labeled using cell molecular tags, and additional single-cell sequencing reagents are introduced to quickly form functional droplet mixes. Lastly, gene expression analysis is used to deduce the impact of the interaction, while molecular sequence tracing identifies the type of interaction. Research into the active effect between antigen-presenting cells and T cells, one of the most common cell-to-cell interactions, is crucial for the advancement of cancer therapy, particularly T cell receptor-engineered T cell therapy. As it allows for high throughput, this platform is superior to well plates as a research platform for cell-to-cell interactions. When combined with the next generation of sequencing, the platform may be able to more accurately evaluate interactions between epitopes and receptors and verify their functional relevance.

A cell transcriptomic profile provides insights into adipocytes of porcine mammary gland across development

BACKGROUND

Studying the composition and developmental mechanisms in mammary gland is crucial for healthy growth of newborns. The mammary gland is inherently heterogeneous, and its physiological function dependents on the gene expression of multiple cell types. Most studies focused on epithelial cells, disregarding the role of neighboring adipocytes.

RESULTS

Here, we constructed the largest transcriptomic dataset of porcine mammary gland cells thus far. The dataset captured 126,829 high-quality nuclei from physiological mammary glands across five developmental stages (d 90 of gestation, G90; d 0 after lactation, L0; d 20 after lactation, L20; 2 d post natural involution, PI2; 7 d post natural involution, PI7). Seven cell types were identified, including epithelial cells, adipocytes, endothelial cells, fibroblasts cells, immune cells, myoepithelial cells and precursor cells. Our data indicate that mammary glands at different developmental stages have distinct phenotypic and transcriptional signatures. During late gestation (G90), the differentiation and proliferation of adipocytes were inhibited. Meanwhile, partly epithelial cells were completely differentiated. Pseudo-time analysis showed that epithelial cells undergo three stages to achieve lactation, including cellular differentiation, hormone sensing, and metabolic activation. During lactation (L0 and L20), adipocytes area accounts for less than 0.5% of mammary glands. To maintain their own survival, the adipocyte exhibited a poorly differentiated state and a proliferative capacity. Epithelial cells initiate lactation upon hormonal stimulation. After fulfilling lactation mission, their undergo physiological death under high intensity lactation. Interestingly, the physiological dead cells seem to be actively cleared by immune cells via CCL21-ACKR4 pathway. This biological process may be an important mechanism for maintaining homeostasis of the mammary gland. During natural involution (PI2 and PI7), epithelial cell populations dedifferentiate into mesenchymal stem cells to maintain the lactation potential of mammary glands for the next lactation cycle.

CONCLUSION

The molecular mechanisms of dedifferentiation, proliferation and redifferentiation of adipocytes and epithelial cells were revealed from late pregnancy to natural involution. This cell transcriptomic profile constitutes an essential reference for future studies in the development and remodeling of the mammary gland at different stages.

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.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

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

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

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

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

Deciphering ligand-receptor-mediated intercellular communication based on ensemble deep learning and the joint scoring strategy from single-cell transcriptomic data

BACKGROUND

Cell-cell communication in a tumor microenvironment is vital to tumorigenesis, tumor progression and therapy. Intercellular communication inference helps understand molecular mechanisms of tumor growth, progression and metastasis.

METHODS

Focusing on ligand-receptor co-expressions, in this study, we developed an ensemble deep learning framework, CellComNet, to decipher ligand-receptor-mediated cell-cell communication from single-cell transcriptomic data. First, credible LRIs are captured by integrating data arrangement, feature extraction, dimension reduction, and LRI classification based on an ensemble of heterogeneous Newton boosting machine and deep neural network. Next, known and identified LRIs are screened based on single-cell RNA sequencing (scRNA-seq) data in certain tissues. Finally, cell-cell communication is inferred by incorporating scRNA-seq data, the screened LRIs, a joint scoring strategy that combines expression thresholding and expression product of ligands and receptors.

RESULTS

The proposed CellComNet framework was compared with four competing protein-protein interaction prediction models (PIPR, XGBoost, DNNXGB, and OR-RCNN) and obtained the best AUCs and AUPRs on four LRI datasets, elucidating the optimal LRI classification ability. CellComNet was further applied to analyze intercellular communication in human melanoma and head and neck squamous cell carcinoma (HNSCC) tissues. The results demonstrate that cancer-associated fibroblasts highly communicate with melanoma cells and endothelial cells strong communicate with HNSCC cells.

CONCLUSIONS

The proposed CellComNet framework efficiently identified credible LRIs and significantly improved cell-cell communication inference performance. We anticipate that CellComNet can contribute to anticancer drug design and tumor-targeted therapy.

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

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

eSPRESSO: topological clustering of single-cell transcriptomics data to reveal informative genes for spatio-temporal architectures of cells

BACKGROUND

Bioinformatics capability to analyze spatio-temporal dynamics of gene expression is essential in understanding animal development. Animal cells are spatially organized as functional tissues where cellular gene expression data contain information that governs morphogenesis during the developmental process. Although several computational tissue reconstruction methods using transcriptomics data have been proposed, those methods have been ineffective in arranging cells in their correct positions in tissues or organs unless spatial information is explicitly provided.

RESULTS

This study demonstrates stochastic self-organizing map clustering with Markov chain Monte Carlo calculations for optimizing informative genes effectively reconstruct any spatio-temporal topology of cells from their transcriptome profiles with only a coarse topological guideline. The method, eSPRESSO (enhanced SPatial REconstruction by Stochastic Self-Organizing Map), provides a powerful in silico spatio-temporal tissue reconstruction capability, as confirmed by using human embryonic heart and mouse embryo, brain, embryonic heart, and liver lobule with generally high reproducibility (average max. accuracy = 92.0%), while revealing topologically informative genes, or spatial discriminator genes. Furthermore, eSPRESSO was used for temporal analysis of human pancreatic organoids to infer rational developmental trajectories with several candidate 'temporal' discriminator genes responsible for various cell type differentiations.

CONCLUSIONS

eSPRESSO provides a novel strategy for analyzing mechanisms underlying the spatio-temporal formation of cellular organizations.

A new insight of immunosuppressive microenvironment in osteosarcoma lung metastasis

The lung is the primary organ for the metastasis of osteosarcoma. Although the application of neoadjuvant chemotherapy and surgery has remarkably improved the survival rate of patients with osteosarcoma, prognosis is still poor for those patients with metastasis. In this study, we performed further bioinformatics analysis on single-cell RNA sequencing (scRNA-seq) data published before, containing 75,317 cells from two osteosarcoma lung metastasis and five normal lung tissues. First, we classified 17 clusters, including macrophages, T cells, endothelial cells, and so on, indicating highly intratumoral heterogeneity in osteosarcoma lung metastasis. Next, we found macrophages in osteosarcoma lung metastasis did not have significant M1 or M2 polarizations. Then, we identified that T cells occupied the most abundant among all cell clusters, and found CD8+ T cells exhibited a low expression level of immune checkpoints in osteosarcoma lung metastasis. What is more, we identified C2_Malignant cells, and found CD63 might play vital roles in determining the infiltration of T cells and malignant cells in conventional-type osteosarcoma lung metastasis. Finally, we unveiled C1_Therapeutic cluster, a subcluster of malignant cells, was sensitive to oxfendazole and mevastatin, and the potential hydrogen-bond position and binding energy of oxfendazole-KIAA0907 and mevastatin-KIAA0907 were unveiled, respectively. Our results highlighted the power of scRNA-seq technique in identifying the complex tumor microenvironment of osteosarcoma lung metastasis, making it possible to devise precision therapeutic approaches.

SPIN-AI: A Deep Learning Model That Identifies Spatially Predictive Genes

Spatially resolved sequencing technologies help us dissect how cells are organized in space. Several available computational approaches focus on the identification of spatially variable genes (SVGs), genes whose expression patterns vary in space. The detection of SVGs is analogous to the identification of differentially expressed genes and permits us to understand how genes and associated molecular processes are spatially distributed within cellular niches. However, the expression activities of SVGs fail to encode all information inherent in the spatial distribution of cells. Here, we devised a deep learning model, Spatially Informed Artificial Intelligence (SPIN-AI), to identify spatially predictive genes (SPGs), whose expression can predict how cells are organized in space. We used SPIN-AI on spatial transcriptomic data from squamous cell carcinoma (SCC) as a proof of concept. Our results demonstrate that SPGs not only recapitulate the biology of SCC but also identify genes distinct from SVGs. Moreover, we found a substantial number of ribosomal genes that were SPGs but not SVGs. Since SPGs possess the capability to predict spatial cellular organization, we reason that SPGs capture more biologically relevant information for a given cellular niche than SVGs. Thus, SPIN-AI has broad applications for detecting SPGs and uncovering which biological processes play important roles in governing cellular organization.

STellaris: a web server for accurate spatial mapping of single cells based on spatial transcriptomics data

Single-cell RNA sequencing (scRNA-seq) provides insights into gene expression heterogeneities in diverse cell types underlying homeostasis, development and pathological states. However, the loss of spatial information hinders its applications in deciphering spatially related features, such as cell-cell interactions in a spatial context. Here, we present STellaris (https://spatial.rhesusbase.com), a web server aimed to rapidly assign spatial information to scRNA-seq data based on their transcriptomic similarity with public spatial transcriptomics (ST) data. STellaris is founded on 101 manually curated ST datasets comprising 823 sections across different organs, developmental stages and pathological states from humans and mice. STellaris accepts raw count matrix and cell type annotation of scRNA-seq data as the input, and maps single cells to spatial locations in the tissue architecture of properly matched ST section. Spatially resolved information for intercellular communications, such as spatial distance and ligand-receptor interactions (LRIs), are further characterized between annotated cell types. Moreover, we also expanded the application of STellaris in spatial annotation of multiple regulatory levels with single-cell multiomics data, using the transcriptome as a bridge. STellaris was applied to several case studies to showcase its utility of adding value to the ever-growing scRNA-seq data from a spatial perspective.

MLSpatial: A machine-learning method to reconstruct the spatial distribution of cells from scRNA-seq by extracting spatial features

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) technologies allow us to interrogate the state of an individual cell within its microenvironment. However, prior to sequencing, cells should be dissociated first, making it difficult to obtain their spatial information. Since the spatial distribution of cells is critical in a few circumstances such as cancer immunotherapy, we present MLSpatial, a novel computational method to learn the relationship between gene expression patterns and spatial locations of cells, and then predict cell-to-cell distance distribution based on scRNA-seq data alone.

RESULTS

We collected the drosophila embryo dataset, which contains both the fluorescence in situ hybridization (FISH) data and single cell RNA-seq (scRNA-seq) data of drosophila embryo. The FISH data provided the spatial position of 3039 cells and the expression of 84 genes for each cell. The scRNA-seq data contains the expressions of 8924 genes in 1297 high-quality cells with cell location unknown. For a comparison, we also collected the MERFISH data of 645 osteosarcoma cells with cell location and the expression status of 10,050 genes known. For each data, the cells were randomly divided into a training set and a test set, in the ratio of 7:3. The cell-to-cell distances our model extracted had a higher correspondence (i.e., correlation coefficient 0.99) with those of the real situation than those of existing methods in the FISH data of drosophila embryo. However, in the osteosarcoma data, our model captured the spatial relationship between cells, with a correlation of 0.514 to that of the real situation. We also applied the model trained using the FISH data of drosophila embryo into the single cell data of drosophila embryo, for which the real location of cells are unknown. The reconstructed pseudo drosophila embryo and the real embryo (as shown by the FISH data) had a high similarity in the spatial distribution of gene expression.

CONCLUSION

MLSpatial can accurately restore the relative position of cells from scRNA-seq data; however, the performance depends on the type of cells. The trained model might be useful in reconstructing the spatial distributions of single cells with only scRNA-seq data, provided that the scRNA-seq data and the FISH data are under similar background (i.e., the same tissue with similar disease background).

Deciphering endothelial heterogeneity in health and disease at single-cell resolution: progress and perspectives

Endothelial cells (ECs) constitute the inner lining of vascular beds in mammals and are crucial for homeostatic regulation of blood vessel physiology, but also play a key role in pathogenesis of many diseases, thereby representing realistic therapeutic targets. However, it has become evident that ECs are heterogeneous, encompassing several subtypes with distinct functions, which makes EC targeting and modulation in diseases challenging. The rise of the new single-cell era has led to an emergence of studies aimed at interrogating transcriptome diversity along the vascular tree, and has revolutionized our understanding of EC heterogeneity from both a physiological and pathophysiological context. Here, we discuss recent landmark studies aimed at teasing apart the heterogeneous nature of ECs. We cover driving (epi)genetic, transcriptomic, and metabolic forces underlying EC heterogeneity in health and disease, as well as current strategies used to combat disease-enriched EC phenotypes, and propose strategies to transcend largely descriptive heterogeneity towards prioritization and functional validation of therapeutically targetable drivers of EC diversity. Lastly, we provide an overview of the most recent advances and hurdles in single EC OMICs.

Orthogonal Light-Dependent Membrane Adhesion Induces Social Self-Sorting and Member-Specific DNA Communication in Synthetic Cell Communities

Developing orthogonal chemical communication pathways in diverse synthetic cell communities is a considerable challenge due to the increased crosstalk and interference associated with large numbers of different types of sender-receiver pairs. Herein, the authors control which sender-receiver pairs communicate in a three-membered community of synthetic cells through red and blue light illumination. Semipermeable protein-polymer-based synthetic cells (proteinosomes) with complementary membrane-attached protein adhesion communicate through single-stranded DNA oligomers and synergistically process biochemical information within a community consisting of one sender and two different receiver populations. Different pairs of red and blue light-responsive protein-protein interactions act as membrane adhesion mediators between the sender and receivers such that they self-assemble and socially self-sort into different multicellular structures under red and blue light. Consequently, distinct sender-receiver pairs come into the signaling range depending on the light illumination and are able to communicate specifically without activation of the other receiver population. Overall, this work shows how photoswitchable membrane adhesion gives rise to different self-sorting protocell patterns that mediate member-specific DNA-based communication in ternary populations of synthetic cells and provides a step towards the design of orthogonal chemical communication networks in diverse communities of synthetic cells.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

Single-cell RNA sequencing reveals intrahepatic and peripheral immune characteristics related to disease phases in HBV-infected patients

OBJECTIVE

A comprehensive immune landscape for HBV infection is pivotal to achieve HBV cure.

DESIGN

We performed single-cell RNA sequencing of 2 43 000 cells from 46 paired liver and blood samples of 23 individuals, including six immune tolerant, 5 immune active (IA), 3 acute recovery (AR), 3 chronic resolved and 6 HBV-free healthy controls (HCs). Flow cytometry and histological assays were applied in a second HBV cohort for validation.

RESULTS

Both IA and AR were characterised by high levels of intrahepatic exhausted CD8+ T (Tex) cells. In IA, Tex cells were mainly derived from liver-resident GZMK+ effector memory T cells and self-expansion. By contrast, peripheral CX3CR1+ effector T cells and GZMK+ effector memory T cells were the main source of Tex cells in AR. In IA but not AR, significant cell-cell interactions were observed between Tex cells and regulatory CD4+ T cells, as well as between Tex and FCGR3A+ macrophages. Such interactions were potentially mediated through human leukocyte antigen class I molecules together with their receptors CANX and LILRBs, respectively, contributing to the dysfunction of antiviral immune responses. By contrast, CX3CR1+GNLY+ central memory CD8+ T cells were concurrently expanded in both liver and blood of AR, providing a potential surrogate marker for viral resolution. In clinic, intrahepatic Tex cells were positively correlated with serum alanine aminotransferase levels and histological grading scores.

CONCLUSION

Our study dissects the coordinated immune responses for different HBV infection phases and provides a rich resource for fully understanding immunopathogenesis and developing effective therapeutic strategies.

Spatially resolved transcriptomics: advances and applications

Spatial transcriptomics, which is capable of both measuring all gene activity in a tissue sample and mapping where this activity occurs, is vastly improving our understanding of biological processes and disease. The field has expanded rapidly in recent years, and the development of several new technologies has resulted in spatially resolved transcriptomics (SRT) becoming highly multiplexed, high-resolution, and high-throughput. Here, we summarize and compare the major methods of SRT, including imaging-based methods, sequencing-based methods, and in situ sequencing methods. We also highlight some typical applications of SRT in neuroscience, cancer biology, developmental biology, and hematology. Finally, we discuss future possibilities for improving spatially resolved transcriptomic methods and the expected applications of such methods, especially in the adult bone marrow, anticipating that new developments will unlock the full potential of spatially resolved multi-omics in both biological research and the clinic.

Promising immunotherapeutic targets in lung cancer based on single-cell RNA sequencing

Immunotherapy has made great strides in the treatment of lung cancer, but a significant proportion of patients still do not respond to treatment. Therefore, the identification of novel targets is crucial to improving the response to immunotherapy. The tumor microenvironment (TME) is a complex niche composed of diverse pro-tumor molecules and cell populations, making the function and mechanism of a unique cell subset difficult to understand. However, the advent of single-cell RNA sequencing (scRNA-seq) technology has made it possible to identify cellular markers and understand their potential functions and mechanisms in the TME. In this review, we highlight recent advances emerging from scRNA-seq studies in lung cancer, with a particular focus on stromal cells. We elucidate the cellular developmental trajectory, phenotypic remodeling, and cell interactions during tumor progression. Our review proposes predictive biomarkers and novel targets for lung cancer immunotherapy based on cellular markers identified through scRNA-seq. The identification of novel targets could help improve the response to immunotherapy. The use of scRNA-seq technology could provide new strategies to understand the TME and develop personalized immunotherapy for lung cancer patients.

Parallel single-cell and bulk transcriptome analyses reveal key features of the gastric tumor microenvironment

BACKGROUND

The tumor microenvironment (TME) has been shown to strongly influence treatment outcome for cancer patients in various indications and to influence the overall survival. However, the cells forming the TME in gastric cancer have not been extensively characterized.

RESULTS

We combine bulk and single-cell RNA sequencing from tumors and matched normal tissue of 24 treatment-naïve GC patients to better understand which cell types and transcriptional programs are associated with malignant transformation of the stomach. Clustering 96,623 cells of non-epithelial origin reveals 81 well-defined TME cell types. We find that activated fibroblasts and endothelial cells are most prominently overrepresented in tumors. Intercellular network reconstruction and survival analysis of an independent cohort imply the importance of these cell types together with immunosuppressive myeloid cell subsets and regulatory T cells in establishing an immunosuppressive microenvironment that correlates with worsened prognosis and lack of response in anti-PD1-treated patients. In contrast, we find a subset of IFNγ activated T cells and HLA-II expressing macrophages that are linked to treatment response and increased overall survival.

CONCLUSIONS

Our gastric cancer single-cell TME compendium together with the matched bulk transcriptome data provides a unique resource for the identification of new potential biomarkers for patient stratification. This study helps further to elucidate the mechanism of gastric cancer and provides insights for therapy.

SOTIP is a versatile method for microenvironment modeling with spatial omics data

The rapidly developing spatial omics generated datasets with diverse scales and modalities. However, most existing methods focus on modeling dynamics of single cells while ignore microenvironments (MEs). Here we present SOTIP (Spatial Omics mulTIPle-task analysis), a versatile method incorporating MEs and their interrelationships into a unified graph. Based on this graph, spatial heterogeneity quantification, spatial domain identification, differential microenvironment analysis, and other downstream tasks can be performed. We validate each module's accuracy, robustness, scalability and interpretability on various spatial omics datasets. In two independent mouse cerebral cortex spatial transcriptomics datasets, we reveal a gradient spatial heterogeneity pattern strongly correlated with the cortical depth. In human triple-negative breast cancer spatial proteomics datasets, we identify molecular polarizations and MEs associated with different patient survivals. Overall, by modeling biologically explainable MEs, SOTIP outperforms state-of-art methods and provides some perspectives for spatial omics data exploration and interpretation.

Single-cell RNA sequencing reveals cellular and molecular reprograming landscape of gliomas and lung cancer brain metastases

BACKGROUND

Brain malignancies encompass gliomas and brain metastases originating from extracranial tumours including lung cancer. Approximately 50% of patients with lung adenocarcinoma (LUAD) will eventually develop brain metastases. However, the specific characteristics of gliomas and lung-to-brain metastases (LC) are largely unknown.

METHODS

We applied single-cell RNA sequencing to profile immune and nonimmune cells in 4 glioma and 10 LC samples.

RESULTS

Our analysis revealed that tumour microenvironment (TME) cells are present in heterogeneous subpopulations. LC reprogramed cells into immune suppressed state, including microglia, macrophages, endothelial cells, and CD8⁺ T cells, with unique cell proportions and gene signatures. Particularly, we identified that a subset of macrophages was associated with poor prognosis. ROS (reactive oxygen species)-producing neutrophils was found to participant in angiogenesis. Furthermore, endothelial cells participated in active communication with fibroblasts. Metastatic epithelial cells exhibited high heterogeneity in chromosomal instability (CIN) and cell population.

CONCLUSIONS

Our findings provide a comprehensive understanding of the heterogenicity of the tumor microenvironment and tumour cells and it will be crucial for successful immunotherapy development for brain metastasis of lung cancer.

PlantPhoneDB: A manually curated pan-plant database of ligand-receptor pairs infers cell-cell communication

Ligand-receptor pairs play important roles in cell-cell communication for multicellular organisms in response to environmental cues. Recently, the emergence of single-cell RNA-sequencing (scRNA-seq) provides unprecedented opportunities to investigate cellular communication based on ligand-receptor expression. However, so far, no reliable ligand-receptor interaction database is available for plant species. In this study, we developed PlantPhoneDB (https://jasonxu.shinyapps.io/PlantPhoneDB/), a pan-plant database comprising a large number of high-confidence ligand-receptor pairs manually curated from seven resources. Also, we developed a PlantPhoneDB R package, which not only provided optional four scoring approaches that calculate interaction scores of ligand-receptor pairs between cell types but also provided visualization functions to present analysis results. At the PlantPhoneDB web interface, the processed datasets and results can be searched, browsed, and downloaded. To uncover novel cell-cell communication events in plants, we applied the PlantPhoneDB R package on GSE121619 dataset to infer significant cell-cell interactions of heat-shocked root cells in Arabidopsis thaliana. As a result, the PlantPhoneDB predicted the actively communicating AT1G28290-AT2G14890 ligand-receptor pair in atrichoblast-cortex cell pair in Arabidopsis thaliana. Importantly, the downstream target genes of this ligand-receptor pair were significantly enriched in the ribosome pathway, which facilitated plants adapting to environmental changes. In conclusion, PlantPhoneDB provided researchers with integrated resources to infer cell-cell communication from scRNA-seq datasets.

Spatial transcriptome profiling of normal human liver

The comprehensive study of the spatial-cellular anatomy of the human liver is critical to addressing the cellular origins of liver disease. Here we conducted spatial transcriptomics on normal human liver tissue sections, providing detailed information of liver zonation at the transcriptional level. We present 6581 high-quality spots from normal livers of two human donors. In this dataset, cells were mainly hepatocytes, and we classified them into four sub-groups. Collectively, these data provide a reliable reference for studies on spatial heterogeneity of liver lobules.

Measurement(s)	mRNA Expression
Technology Type(s)	spatial transcriptomics
Sample Characteristic - Organism	Homo sapiens

Recent developments in application of single-cell RNA sequencing in the tumour immune microenvironment and cancer therapy

The advent of single-cell RNA sequencing (scRNA-seq) has provided insight into the tumour immune microenvironment (TIME). This review focuses on the application of scRNA-seq in investigation of the TIME. Over time, scRNA-seq methods have evolved, and components of the TIME have been deciphered with high resolution. In this review, we first introduced the principle of scRNA-seq and compared different sequencing approaches. Novel cell types in the TIME, a continuous transitional state, and mutual intercommunication among TIME components present potential targets for prognosis prediction and treatment in cancer. Thus, we concluded novel cell clusters of cancer-associated fibroblasts (CAFs), T cells, tumour-associated macrophages (TAMs) and dendritic cells (DCs) discovered after the application of scRNA-seq in TIME. We also proposed the development of TAMs and exhausted T cells, as well as the possible targets to interrupt the process. In addition, the therapeutic interventions based on cellular interactions in TIME were also summarized. For decades, quantification of the TIME components has been adopted in clinical practice to predict patient survival and response to therapy and is expected to play an important role in the precise treatment of cancer. Summarizing the current findings, we believe that advances in technology and wide application of single-cell analysis can lead to the discovery of novel perspectives on cancer therapy, which can subsequently be implemented in the clinic. Finally, we propose some future directions in the field of TIME studies that can be aided by scRNA-seq technology.

Delineating COVID-19 immunological features using single-cell RNA sequencing

Understanding the molecular mechanisms of coronavirus disease 2019 (COVID-19) pathogenesis and immune response is vital for developing therapies. Single-cell RNA sequencing has been applied to delineate the cellular heterogeneity of the host response toward COVID-19 in multiple tissues and organs. Here, we review the applications and findings from over 80 original COVID-19 single-cell RNA sequencing studies as well as many secondary analysis studies. We describe that single-cell RNA sequencing reveals multiple features of COVID-19 patients with different severity, including cell populations with proportional alteration, COVID-19-induced genes and pathways, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection in single cells, and adaptation of immune repertoire. We also collect published single-cell RNA sequencing datasets from original studies. Finally, we discuss the limitations in current studies and perspectives for future advance.

Scalable and model-free detection of spatial patterns and colocalization

The expeditious growth in spatial omics technologies enables the profiling of genome-wide molecular events at molecular and single-cell resolution, highlighting a need for fast and reliable methods to characterize spatial patterns. We developed SpaGene, a model-free method to discover spatial patterns rapidly in large-scale spatial omics studies. Analyzing simulation and a variety of spatially resolved transcriptomics data showed that SpaGene is more powerful and scalable than existing methods. Spatial expression patterns identified by SpaGene reconstruct unobserved tissue structures. SpaGene also successfully discovers ligand-receptor interactions through their colocalization.

Tumor vessel co-option: The past & the future

Tumor vessel co-option (VCO) is a non-angiogenic vascularization mechanism that is a possible cause of resistance to anti-angiogenic therapy (AAT). Multiple tumors are hypothesized to primarily rely on growth factor signaling-induced sprouting angiogenesis, which is often inhibited during AAT. During VCO however, tumors invade healthy tissues by hijacking pre-existing blood vessels of the host organ to secure their blood and nutrient supply. Although VCO has been described in the context of AAT resistance, the molecular mechanisms underlying this process and the profile and characteristics of co-opted vascular cell types (endothelial cells (ECs) and pericytes) remain poorly understood, resulting in the lack of therapeutic strategies to inhibit VCO (and to overcome AAT resistance). In the past few years, novel next-generation technologies (such as single-cell RNA sequencing) have emerged and revolutionized the way of analyzing and understanding cancer biology. While most studies utilizing single-cell RNA sequencing with focus on cancer vascularization have centered around ECs during sprouting angiogenesis, we propose that this and other novel technologies can be used in future investigations to shed light on tumor EC biology during VCO. In this review, we summarize the molecular mechanisms driving VCO known to date and introduce the models used to study this phenomenon to date. We highlight VCO studies that recently emerged using sequencing approaches and propose how these and other novel state-of-the-art methods can be used in the future to further explore ECs and other cell types in the VCO process and to identify potential vulnerabilities in tumors relying on VCO. A better understanding of VCO by using novel approaches could provide new answers to the many open questions, and thus pave the way to develop new strategies to control and target tumor vascularization.

Single-Cell Transcriptomic Analysis Reveals the Crosstalk Propensity Between the Tumor Intermediate State and the CD8+ T Exhausted State to be Associated with Clinical Benefits in Melanoma

Heterogeneous crosstalk between tumor cells and CD8+ T cells leads to substantial variation in clinical benefits from immunotherapy in melanoma. Due to spatial distribution and functional state heterogeneity, it is still unknown whether there is a crosstalk propensity between tumor cells and CD8+ T cells in melanoma, and how this crosstalk propensity affects the clinical outcome of patients. Using public single-cell transcriptome data, extensive heterogeneous functional states and ligand–receptor interactions of tumor cells and CD8+ T cells were revealed in melanoma. Furthermore, based on the association between cell–cell communication intensity and cell state activity in a single cell, we identified a crosstalk propensity between the tumor intermediate state and the CD8+ T exhausted state. This crosstalk propensity was further verified by pseudo-spatial proximity, spatial co-location, and the intra/intercellular signal transduction network. At the sample level, the tumor intermediate state and the CD8+ T exhausted state synergistically indicated better prognosis and both reduced in immunotherapy-resistant samples. The risk groups defined based on these two cell states could comprehensively reflect tumor genomic mutations and anti-tumor immunity information. The low-risk group had a higher BRAF mutation fraction as well as stronger antitumor immune response. Our findings highlighted the crosstalk propensity between the tumor intermediate state and the CD8+ T exhausted state, which may serve as a reference to guide the development of diagnostic biomarkers for risk stratification and therapeutic targets for new therapeutic strategies.

SPROUT: spectral sparsification helps restore the spatial structure at single-cell resolution

Single-cell RNA sequencing thoroughly quantifies the individual cell transcriptomes but renounces the spatial structure. Conversely, recently emerged spatial transcriptomics technologies capture the cellular spatial structure but skimp cell or gene resolutions. Ligand-receptor interactions reveal the potential of cell proximity since they are spatially constrained. Cell–cell affinity values estimated by ligand–receptor interaction can partially represent the structure of cells but falsely include the pseudo affinities between distant or indirectly interacting cells. Here, we develop a software package, SPROUT, to reconstruct the single-cell resolution spatial structure from the transcriptomics data through diminished pseudo ligand–receptor affinities. For spatial data, SPROUT first curates the representative single-cell profiles for each spatial spot from a candidate library, then reduces the pseudo affinities in the intercellular affinity matrix by partial correlation, spectral graph sparsification, and spatial coordinates refinement. SPROUT embeds the estimated interactions into a low-dimensional space with the cross-entropy objective to restore the intercellular structures, which facilitates the discovery of dominant ligand–receptor pairs between neighboring cells at single-cell resolution. SPROUT reconstructed structures achieved shape Pearson correlations ranging from 0.91 to 0.97 on the mouse hippocampus and human organ tumor microenvironment datasets. Furthermore, SPROUT can solely de novo reconstruct the structures at single-cell resolution, i.e., reaching the cell-type proximity correlations of 0.68 and 0.89 between reconstructed and immunohistochemistry-informed spatial structures on a human developing heart dataset and a tumor microenvironment dataset, respectively.

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.

Integration of computational analysis and spatial transcriptomics in single-cell study

Recent advances of single-cell transcriptomics technologies and allied computational methodologies have revolutionized molecular cell biology. Meanwhile, pioneering explorations in spatial transcriptomics have opened avenues to address fundamental biological questions in health and diseases. Here, we review the technical attributes of single-cell RNA sequencing and spatial transcriptomics, and the core concepts of computational data analysis. We further highlight the challenges in the application of data integration methodologies and the interpretation of the biological context of the findings.

Single-Cell RNA Sequencing Unravels Upregulation of Immune Cell Crosstalk in Relapsed Pediatric Ependymoma

Ependymoma (EPN) is a malignant glial tumor occurring throughout central nervous system, which commonly presents in children. Although recent studies have characterized EPN samples at both the bulk and single-cell level, intratumoral heterogeneity across subclones remains a confounding factor that impedes understanding of EPN biology. In this study, we generated a high-resolution single-cell dataset of pediatric ependymoma with a particular focus on the comparison of subclone differences within tumors and showed upregulation of cilium-associated genes in more highly differentiated subclone populations. As a proxy to traditional pseudotime analysis, we applied a novel trajectory scoring method to reveal cellular compositions associated with poor survival outcomes across primary and relapsed patients. Furthermore, we identified putative cell–cell communication features between relapsed and primary samples and showed upregulation of pathways associated with immune cell crosstalk. Our results revealed both inter- and intratumoral heterogeneity in EPN and provided a framework for studying transcriptomic signatures of individual subclones at single-cell resolution.

Discovering Immune-Mediated Mechanisms of Gastric Carcinogenesis Through Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNAseq) technology is still relatively new in the field of gastric cancer immunology but gaining significant traction. This technology now provides unprecedented insights into the intratumoral and intertumoral heterogeneities at the immunological, cellular, and molecular levels. Within the last few years, a volume of publications reported the usefulness of scRNAseq technology in identifying thus far elusive immunological mechanisms that may promote and impede gastric cancer development. These studies analyzed datasets generated from primary human gastric cancer tissues, metastatic ascites fluid from gastric cancer patients, and laboratory-generated data from in vitro and in vivo models of gastric diseases. In this review, we overview the exciting findings from scRNAseq datasets that uncovered the role of critical immune cells, including T cells, B cells, myeloid cells, mast cells, ILC2s, and other inflammatory stromal cells, like fibroblasts and endothelial cells. In addition, we also provide a synopsis of the initial scRNAseq findings on the interesting epithelial cell responses to inflammation. In summary, these new studies have implicated roles for T and B cells and subsets like NKT cells in tumor development and progression. The current studies identified diverse subsets of macrophages and mast cells in the tumor microenvironment, however, additional studies to determine their roles in promoting cancer growth are needed. Some groups specifically focus on the less prevalent ILC2 cell type that may contribute to early cancer development. ScRNAseq analysis also reveals that stromal cells, e.g., fibroblasts and endothelial cells, regulate inflammation and promote metastasis, making them key targets for future investigations. While evaluating the outcomes, we also highlight the gaps in the current findings and provide an assessment of what this technology holds for gastric cancer research in the coming years. With scRNAseq technology expanding rapidly, we stress the need for periodic review of the findings and assess the available scRNAseq analytical tools to guide future work on immunological mechanisms of gastric carcinogenesis.  

An introduction to spatial transcriptomics for biomedical research

Single-cell transcriptomics (scRNA-seq) has become essential for biomedical research over the past decade, particularly in developmental biology, cancer, immunology, and neuroscience. Most commercially available scRNA-seq protocols require cells to be recovered intact and viable from tissue. This has precluded many cell types from study and largely destroys the spatial context that could otherwise inform analyses of cell identity and function. An increasing number of commercially available platforms now facilitate spatially resolved, high-dimensional assessment of gene transcription, known as 'spatial transcriptomics'. Here, we introduce different classes of method, which either record the locations of hybridized mRNA molecules in tissue, image the positions of cells themselves prior to assessment, or employ spatial arrays of mRNA probes of pre-determined location. We review sizes of tissue area that can be assessed, their spatial resolution, and the number and types of genes that can be profiled. We discuss if tissue preservation influences choice of platform, and provide guidance on whether specific platforms may be better suited to discovery screens or hypothesis testing. Finally, we introduce bioinformatic methods for analysing spatial transcriptomic data, including pre-processing, integration with existing scRNA-seq data, and inference of cell-cell interactions. Spatial -omics methods are already improving our understanding of human tissues in research, diagnostic, and therapeutic settings. To build upon these recent advancements, we provide entry-level guidance for those seeking to employ spatial transcriptomics in their own biomedical research.

Microfluidics Facilitates the Development of Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) technology provides a powerful tool for understanding complex biosystems at the single-cell and single-molecule level. The past decade has been a golden period for the development of single-cell sequencing, with scRNA-seq undergoing a tremendous leap in sensitivity and throughput. The application of droplet- and microwell-based microfluidics in scRNA-seq has contributed greatly to improving sequencing throughput. This review introduces the history of development and important technical factors of scRNA-seq. We mainly focus on the role of microfluidics in facilitating the development of scRNA-seq technology. To end, we discuss the future directions for scRNA-seq.