Cell lineage and communication network inference via optimization for single-cell transcriptomics

Single-cell RNA sequencing technologies and applications: A brief overview

Single-cell RNA sequencing (scRNA-seq) technology has become the state-of-the-art approach for unravelling the heterogeneity and complexity of RNA transcripts within individual cells, as well as revealing the composition of different cell types and functions within highly organized tissues/organs/organisms. Since its first discovery in 2009, studies based on scRNA-seq provide massive information across different fields making exciting new discoveries in better understanding the composition and interaction of cells within humans, model animals and plants. In this review, we provide a concise overview about the scRNA-seq technology, experimental and computational procedures for transforming the biological and molecular processes into computational and statistical data. We also provide an explanation of the key technological steps in implementing the technology. We highlight a few examples on how scRNA-seq can provide unique information for better understanding health and diseases. One important application of the scRNA-seq technology is to build a better and high-resolution catalogue of cells in all living organism, commonly known as atlas, which is key resource to better understand and provide a solution in treating diseases. While great promises have been demonstrated with the technology in all areas, we further highlight a few remaining challenges to be overcome and its great potentials in transforming current protocols in disease diagnosis and treatment.

Computational exploration of cellular communication in skin from emerging single-cell and spatial transcriptomic data

Tissue development and homeostasis require coordinated cell-cell communication. Recent advances in single-cell sequencing technologies have emerged as a revolutionary method to reveal cellular heterogeneity with unprecedented resolution. This offers a great opportunity to explore cell-cell communication in tissues systematically and comprehensively, and to further identify signaling mechanisms driving cell fate decisions and shaping tissue phenotypes. Using gene expression information from single-cell transcriptomics, several computational tools have been developed for inferring cell-cell communication, greatly facilitating analysis and interpretation. However, in single-cell transcriptomics, spatial information of cells is inherently lost. Given that most cell signaling events occur within a limited distance in tissues, incorporating spatial information into cell-cell communication analysis is critical for understanding tissue organization and function. Spatial transcriptomics provides spatial location of cell subsets along with their gene expression, leading to new directions for leveraging spatial information to develop computational approaches for cell-cell communication inference and analysis. These computational approaches have been successfully applied to uncover previously unrecognized mechanisms of intercellular communication within various contexts and across organ systems, including the skin, a formidable model to study mechanisms of cell-cell communication due to the complex interactions between the different cell populations that comprise it. Here, we review emergent cell-cell communication inference tools using single-cell transcriptomics and spatial transcriptomics, and highlight the biological insights gained by applying these computational tools to exploring cellular communication in skin development, homeostasis, disease and aging, as well as discuss future potential research avenues.

Cell–Cell Communication Alterations via Intercellular Signaling Pathways in Substantia Nigra of Parkinson’s Disease

Parkinson’s disease (PD) is a neurodegenerative movement disorder characterized with dopaminergic neuron (DaN) loss within the substantia nigra (SN). Despite bulk studies focusing on intracellular mechanisms of PD inside DaNs, few studies have explored the pathogeneses outside DaNs, or between DaNs and other cells. Here, we set out to probe the implication of intercellular communication involving DaNs in the pathogeneses of PD at a systemic level with bioinformatics methods. We harvested three online published single-cell/single-nucleus transcriptomic sequencing (sc/snRNA-seq) datasets of human SN (GSE126838, GSE140231, and GSE157783) from the Gene Expression Omnibus (GEO) database, and integrated them with one of the latest integration algorithms called Harmony. We then applied CellChat, the latest cell–cell communication analytic algorithm, to our integrated dataset. We first found that the overall communication quantity was decreased while the overall communication strength was enhanced in PD sample compared with control sample. We then focused on the intercellular communication where DaNs are involved, and found that the communications between DaNs and other cell types via certain signaling pathways were selectively altered in PD, including some growth factors, neurotrophic factors, chemokines, etc. pathways. Our bioinformatics analysis showed that the alteration in intercellular communications involving DaNs might be a previously underestimated aspect of PD pathogeneses with novel translational potential.

The use of machine learning to discover regulatory networks controlling biological systems

Biological systems are composed of a vast web of multiscale molecular interactors and interactions. High-throughput technologies, both bulk and single cell, now allow for investigation of the properties and quantities of these interactors. Computational algorithms and machine learning methods then provide the tools to derive meaningful insights from the resulting data sets. One such approach is graphical network modeling, which provides a computational framework to explicitly model the molecular interactions within and between the cells comprising biological systems. These graphical networks aim to describe a putative chain of cause and effect between interacting molecules. This feature allows for determination of key molecules in a biological process, accelerated generation of mechanistic hypotheses, and simulation of experimental outcomes. We review the computational concepts and applications of graphical network models across molecular scales for both intracellular and intercellular regulatory biology, examples of successful applications, and the future directions needed to overcome current limitations.

Network-based integrative analysis of single-cell transcriptomic and epigenomic data for cell types

Advances in single-cell biotechnologies simultaneously generate the transcriptomic and epigenomic profiles at cell levels, providing an opportunity for investigating cell fates. Although great efforts have been devoted to either of them, the integrative analysis of single-cell multi-omics data is really limited because of the heterogeneity, noises and sparsity of single-cell profiles. In this study, a network-based integrative clustering algorithm (aka NIC) is present for the identification of cell types by fusing the parallel single-cell transcriptomic (scRNA-seq) and epigenomic profiles (scATAC-seq or DNA methylation). To avoid heterogeneity of multi-omics data, NIC automatically learns the cell–cell similarity graphs, which transforms the fusion of multi-omics data into the analysis of multiple networks. Then, NIC employs joint non-negative matrix factorization to learn the shared features of cells by exploiting the structure of learned cell–cell similarity networks, providing a better way to characterize the features of cells. The graph learning and integrative analysis procedures are jointly formulated as an optimization problem, and then the update rules are derived. Thirteen single-cell multi-omics datasets from various tissues and organisms are adopted to validate the performance of NIC, and the experimental results demonstrate that the proposed algorithm significantly outperforms the state-of-the-art methods in terms of various measurements. The proposed algorithm provides an effective strategy for the integrative analysis of single-cell multi-omics data (The software is coded using Matlab, and is freely available for academic https://github.com/xkmaxidian/NIC ).

Spatial omics: Navigating to the golden era of cancer research

The idea that tumour microenvironment (TME) is organised in a spatial manner will not surprise many cancer biologists; however, systematically capturing spatial architecture of TME is still not possible until recent decade. The past five years have witnessed a boom in the research of high-throughput spatial techniques and algorithms to delineate TME at an unprecedented level. Here, we review the technological progress of spatial omics and how advanced computation methods boost multi-modal spatial data analysis. Then, we discussed the potential clinical translations of spatial omics research in precision oncology, and proposed a transfer of spatial ecological principles to cancer biology in spatial data interpretation. So far, spatial omics is placing us in the golden age of spatial cancer research. Further development and application of spatial omics may lead to a comprehensive decoding of the TME ecosystem and bring the current spatiotemporal molecular medical research into an entirely new paradigm.

Single-cell transcriptomic analysis of zebrafish cranial neural crest reveals spatiotemporal regulation of lineage decisions during development

Neural crest (NC) cells migrate throughout vertebrate embryos to give rise to a huge variety of cell types, but when and where lineages emerge and their regulation remain unclear. We have performed single-cell RNA sequencing (RNA-seq) of cranial NC cells from the first pharyngeal arch in zebrafish over several stages during migration. Computational analysis combining pseudotime and real-time data reveals that these NC cells first adopt a transitional state, becoming specified mid-migration, with the first lineage decisions being skeletal and pigment, followed by neural and glial progenitors. In addition, by computationally integrating these data with RNA-seq data from a transgenic Wnt reporter line, we identify gene cohorts with similar temporal responses to Wnts during migration and show that one, Atp6ap2, is required for melanocyte differentiation. Together, our results show that cranial NC cell lineages arise progressively and uncover a series of spatially restricted cell interactions likely to regulate such cell-fate decisions.

Tatarakis et al. provide a single-cell transcriptomic timeline of cranial neural crest (NC) development in zebrafish and address long-standing questions surrounding the integration of NC cell migration and lineage specification. They find that lineages are specified mid-migration. These fate decisions correspond to shifts in Wnt signaling, and lineages rapidly segregate.

Single-cell RNA sequencing reveals the epithelial cell heterogeneity and invasive subpopulation in human bladder cancer

Bladder cancer represents a highly heterogeneous disease characterized by distinct histological, molecular and clinical phenotypes, and a detailed analysis of tumor cell invasion and crosstalks within bladder tumor cells has not been determined. Here, we applied droplet-based single-cell RNA sequencing (scRNA-seq) to acquire transcriptional profiles of 36 619 single cells isolated from seven patients. Single cell transcriptional profiles matched well with the pathological basal/luminal subtypes. Notably, in T1 tumors diagnosed as luminal subtype, basal cells displayed characteristics of epithelial-mesenchymal transition (EMT) and mainly located at the tumor-stromal interface as well as micrometastases in the lamina propria. In one T3 tumor, muscle-invasive tumor showed significantly higher expression of cancer stem cell markers SOX9 and SOX2 than the primary tumor. We additionally analyzed communications between tumor cells and demonstrated its relevance to basal/luminal phenotypes. Overall, our single-cell study provides a deeper insight into the tumor cell heterogeneity associated with bladder cancer progression.

Intracellular and Intercellular Gene Regulatory Network Inference From Time-Course Individual RNA-Seq

Gene regulatory network (GRN) inference is an effective approach to understand the molecular mechanisms underlying biological events. Generally, GRN inference mainly targets intracellular regulatory relationships such as transcription factors and their associated targets. In multicellular organisms, there are both intracellular and intercellular regulatory mechanisms. Thus, we hypothesize that GRNs inferred from time-course individual (whole embryo) RNA-Seq during development can reveal intercellular regulatory relationships (signaling pathways) underlying the development. Here, we conducted time-course bulk RNA-Seq of individual mouse embryos during early development, followed by pseudo-time analysis and GRN inference. The results demonstrated that GRN inference from RNA-Seq with pseudo-time can be applied for individual bulk RNA-Seq similar to scRNA-Seq. Validation using an experimental-source-based database showed that our approach could significantly infer GRN for all transcription factors in the database. Furthermore, the inferred ligand-related and receptor-related downstream genes were significantly overlapped. Thus, the inferred GRN based on whole organism could include intercellular regulatory relationships, which cannot be inferred from scRNA-Seq based only on gene expression data. Overall, inferring GRN from time-course bulk RNA-Seq is an effective approach to understand the regulatory relationships underlying biological events in multicellular organisms.

Integrating single-cell genomics pipelines to discover mechanisms of stem cell differentiation

Pluripotent stem cells underpin a growing sector that leverages their differentiation potential for research, industry, and clinical applications. This review evaluates the landscape of methods in single-cell transcriptomics that are enabling accelerated discovery in stem cell science. We focus on strategies for scaling stem cell differentiation through multiplexed single-cell analyses, for evaluating molecular regulation of cell differentiation using new analysis algorithms, and methods for integration and projection analysis to classify and benchmark stem cell derivatives against in vivo cell types. By discussing the available methods, comparing their strengths, and illustrating strategies for developing integrated analysis pipelines, we provide user considerations to inform their implementation and interpretation.

Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics

Single-cell RNA sequencing (scRNA-seq) identifies cell subpopulations within tissue but does not capture their spatial distribution nor reveal local networks of intercellular communication acting in situ. A suite of recently developed techniques that localize RNA within tissue, including multiplexed in situ hybridization and in situ sequencing (here defined as high-plex RNA imaging) and spatial barcoding, can help address this issue. However, no method currently provides as complete a scope of the transcriptome as does scRNA-seq, underscoring the need for approaches to integrate single-cell and spatial data. Here, we review efforts to integrate scRNA-seq with spatial transcriptomics, including emerging integrative computational methods, and propose ways to effectively combine current methodologies.

CellCall: integrating paired ligand-receptor and transcription factor activities for cell-cell communication

With the dramatic development of single-cell RNA sequencing (scRNA-seq) technologies, the systematic decoding of cell-cell communication has received great research interest. To date, several in-silico methods have been developed, but most of them lack the ability to predict the communication pathways connecting the insides and outsides of cells. Here, we developed CellCall, a toolkit to infer inter- and intracellular communication pathways by integrating paired ligand-receptor and transcription factor (TF) activity. Moreover, CellCall uses an embedded pathway activity analysis method to identify the significantly activated pathways involved in intercellular crosstalk between certain cell types. Additionally, CellCall offers a rich suite of visualization options (Circos plot, Sankey plot, bubble plot, ridge plot, etc.) to present the analysis results. Case studies on scRNA-seq datasets of human testicular cells and the tumor immune microenvironment demonstrated the reliable and unique functionality of CellCall in intercellular communication analysis and internal TF activity exploration, which were further validated experimentally. Comparative analysis of CellCall and other tools indicated that CellCall was more accurate and offered more functions. In summary, CellCall provides a sophisticated and practical tool allowing researchers to decipher intercellular communication and related internal regulatory signals based on scRNA-seq data. CellCall is freely available at https://github.com/ShellyCoder/cellcall.

jSRC: a flexible and accurate joint learning algorithm for clustering of single-cell RNA-sequencing data

Single-cell RNA-sequencing (scRNA-seq) explores the transcriptome of genes at cell level, which sheds light on revealing the heterogeneity and dynamics of cell populations. Advances in biotechnologies make it possible to generate scRNA-seq profiles for large-scale cells, requiring effective and efficient clustering algorithms to identify cell types and informative genes. Although great efforts have been devoted to clustering of scRNA-seq, the accuracy, scalability and interpretability of available algorithms are not desirable. In this study, we solve these problems by developing a joint learning algorithm [a.k.a. joints sparse representation and clustering (jSRC)], where the dimension reduction (DR) and clustering are integrated. Specifically, DR is employed for the scalability and joint learning improves accuracy. To increase the interpretability of patterns, we assume that cells within the same type have similar expression patterns, where the sparse representation is imposed on features. We transform clustering of scRNA-seq into an optimization problem and then derive the update rules to optimize the objective of jSRC. Fifteen scRNA-seq datasets from various tissues and organisms are adopted to validate the performance of jSRC, where the number of single cells varies from 49 to 110 824. The experimental results demonstrate that jSRC significantly outperforms 12 state-of-the-art methods in terms of various measurements (on average 20.29% by improvement) with fewer running time. Furthermore, jSRC is efficient and robust across different scRNA-seq datasets from various tissues. Finally, jSRC also accurately identifies dynamic cell types associated with progression of COVID-19. The proposed model and methods provide an effective strategy to analyze scRNA-seq data (the software is coded using MATLAB and is free for academic purposes; https://github.com/xkmaxidian/jSRC).

Applications and analytical tools of cell communication based on ligand-receptor interactions at single cell level

BACKGROUND

Cellular communication is an essential feature of multicellular organisms. Binding of ligands to their homologous receptors, which activate specific cell signaling pathways, is a basic type of cellular communication and intimately linked to many degeneration processes leading to diseases.

MAIN BODY

This study reviewed the history of ligand-receptor and presents the databases which store ligand-receptor pairs. The recently applications and research tools of ligand-receptor interactions for cell communication at single cell level by using single cell RNA sequencing have been sorted out.

CONCLUSION

The summary of the advantages and disadvantages of analysis tools will greatly help researchers analyze cell communication at the single cell level. Learning cell communication based on ligand-receptor interactions by single cell RNA sequencing gives way to developing new target drugs and personalizing treatment.

The landscape of cell-cell communication through single-cell transcriptomics

Cell-cell communication is a fundamental process that shapes biological tissue. Historically, studies of cell-cell communication have been feasible for one or two cell types and a few genes. With the emergence of single-cell transcriptomics, we are now able to examine the genetic profiles of individual cells at unprecedented scale and depth. The availability of such data presents an exciting opportunity to construct a more comprehensive description of cell-cell communication. This review discusses the recent explosion of methods that have been developed to infer cell-cell communication from non-spatial and spatial single-cell transcriptomics, two promising technologies which have complementary strengths and limitations. We propose several avenues to propel this rapidly expanding field forward in meaningful ways.

An in vitro vascularized micro-tumor model of human colorectal cancer recapitulates in vivo responses to standard-of-care therapy

Around 95% of anti-cancer drugs that show promise during preclinical study fail to gain FDA-approval for clinical use. This failure of the preclinical pipeline highlights the need for improved, physiologically-relevant in vitro models that can better serve as reliable drug-screening and disease modeling tools. The vascularized micro-tumor (VMT) is a novel three-dimensional model system (tumor-on-a-chip) that recapitulates the complex human tumor microenvironment, including perfused vasculature, within a transparent microfluidic device, allowing real-time study of drug responses and tumor-stromal interactions. Here we have validated this microphysiological system (MPS) platform for the study of colorectal cancer (CRC), the second leading cause of cancer-related deaths, by showing that gene expression, tumor heterogeneity, and treatment responses in the VMT more closely model CRC tumor clinicopathology than current standard drug screening modalities, including 2-dimensional monolayer culture and 3-dimensional spheroids.

DTFLOW: Inference and Visualization of Single-cell Pseudotime Trajectory Using Diffusion Propagation

One of the major challenges in single-cell data analysis is the determination of cellular developmental trajectories using single-cell data. Although substantial studies have been conducted in recent years, more effective methods are still strongly needed to infer the developmental processes accurately. This work devises a new method, named DTFLOW, for determining the pseudo-temporal trajectories with multiple branches. DTFLOW consists of two major steps: a new method called Bhattacharyya kernel feature decomposition (BKFD) to reduce the data dimensions, and a novel approach named Reverse Searching on k-nearest neighbor graph (RSKG) to identify the multi-branching processes of cellular differentiation. In BKFD, we first establish a stationary distribution for each cell to represent the transition of cellular developmental states based on the random walk with restart algorithm, and then propose a new distance metric for calculating pseudotime of single cells by introducing the Bhattacharyya kernel matrix. The effectiveness of DTFLOW is rigorously examined by using four single-cell datasets. We compare the efficiency of DTFLOW with the published state-of-the-art methods. Simulation results suggest that DTFLOW has superior accuracy and strong robustness properties for constructing pseudotime trajectories. The Python source code of DTFLOW can be freely accessed at https://github.com/statway/DTFLOW.

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

Single-cell technology enables study of signal transduction in a complex tissue at unprecedented resolution. We describe CytoTalk for de novo construction of cell type-specific signaling networks using single-cell transcriptomic data. Using an integrated intracellular and intercellular gene network as the input, CytoTalk identifies candidate pathways using the prize-collecting Steiner forest algorithm. Using high-throughput spatial transcriptomic data and single-cell RNA sequencing data with receptor gene perturbation, we demonstrate that CytoTalk has substantial improvement over existing algorithms. To better understand plasticity of signaling networks across tissues and developmental stages, we perform a comparative analysis of signaling networks between macrophages and endothelial cells across human adult and fetal tissues. Our analysis reveals an overall increased plasticity of signaling networks across adult tissues and specific network nodes that contribute to increased plasticity. CytoTalk enables de novo construction of signal transduction pathways and facilitates comparative analysis of these pathways across tissues and conditions.

DEEPsc: A Deep Learning-Based Map Connecting Single-Cell Transcriptomics and Spatial Imaging Data

Single-cell RNA sequencing (scRNA-seq) data provides unprecedented information on cell fate decisions; however, the spatial arrangement of cells is often lost. Several recent computational methods have been developed to impute spatial information onto a scRNA-seq dataset through analyzing known spatial expression patterns of a small subset of genes known as a reference atlas. However, there is a lack of comprehensive analysis of the accuracy, precision, and robustness of the mappings, along with the generalizability of these methods, which are often designed for specific systems. We present a system-adaptive deep learning-based method (DEEPsc) to impute spatial information onto a scRNA-seq dataset from a given spatial reference atlas. By introducing a comprehensive set of metrics that evaluate the spatial mapping methods, we compare DEEPsc with four existing methods on four biological systems. We find that while DEEPsc has comparable accuracy to other methods, an improved balance between precision and robustness is achieved. DEEPsc provides a data-adaptive tool to connect scRNA-seq datasets and spatial imaging datasets to analyze cell fate decisions. Our implementation with a uniform API can serve as a portal with access to all the methods investigated in this work for spatial exploration of cell fate decisions in scRNA-seq data. All methods evaluated in this work are implemented as an open-source software with a uniform interface.

Inferring microenvironmental regulation of gene expression from single-cell RNA sequencing data using scMLnet with an application to COVID-19

Inferring how gene expression in a cell is influenced by cellular microenvironment is of great importance yet challenging. In this study, we present a single-cell RNA-sequencing data based multilayer network method (scMLnet) that models not only functional intercellular communications but also intracellular gene regulatory networks (https://github.com/SunXQlab/scMLnet). scMLnet was applied to a scRNA-seq dataset of COVID-19 patients to decipher the microenvironmental regulation of expression of SARS-CoV-2 receptor ACE2 that has been reported to be correlated with inflammatory cytokines and COVID-19 severity. The predicted elevation of ACE2 by extracellular cytokines EGF, IFN-γ or TNF-α were experimentally validated in human lung cells and the related signaling pathway were verified to be significantly activated during SARS-COV-2 infection. Our study provided a new approach to uncover inter-/intra-cellular signaling mechanisms of gene expression and revealed microenvironmental regulators of ACE2 expression, which may facilitate designing anti-cytokine therapies or targeted therapies for controlling COVID-19 infection. In addition, we summarized and compared different methods of scRNA-seq based inter-/intra-cellular signaling network inference for facilitating new methodology development and applications.

Dissection of intercellular communication using the transcriptome-based framework ICELLNET

Cell-to-cell communication can be inferred from ligand–receptor expression in cell transcriptomic datasets. However, important challenges remain: global integration of cell-to-cell communication; biological interpretation; and application to individual cell population transcriptomic profiles. We develop ICELLNET, a transcriptomic-based framework integrating: 1) an original expert-curated database of ligand–receptor interactions accounting for multiple subunits expression; 2) quantification of communication scores; 3) the possibility to connect a cell population of interest with 31 reference human cell types; and 4) three visualization modes to facilitate biological interpretation. We apply ICELLNET to three datasets generated through RNA-seq, single-cell RNA-seq, and microarray. ICELLNET reveals autocrine IL-10 control of human dendritic cell communication with up to 12 cell types. Four of them (T cells, keratinocytes, neutrophils, pDC) are further tested and experimentally validated. In summary, ICELLNET is a global, versatile, biologically validated, and easy-to-use framework to dissect cell communication from individual or multiple cell-based transcriptomic profiles.

Bulk and single-cell transcriptomic data can be a source of novel insights into how cells interact with each other. Here the authors develop ICELLNET, a global, biologically validated, and easy-to-use framework to dissect cell communication from individual or multiple cell-based transcriptomic profiles.

Inference and analysis of cell-cell communication using CellChat

Understanding global communications among cells requires accurate representation of cell-cell signaling links and effective systems-level analyses of those links. We construct a database of interactions among ligands, receptors and their cofactors that accurately represent known heteromeric molecular complexes. We then develop CellChat, a tool that is able to quantitatively infer and analyze intercellular communication networks from single-cell RNA-sequencing (scRNA-seq) data. CellChat predicts major signaling inputs and outputs for cells and how those cells and signals coordinate for functions using network analysis and pattern recognition approaches. Through manifold learning and quantitative contrasts, CellChat classifies signaling pathways and delineates conserved and context-specific pathways across different datasets. Applying CellChat to mouse and human skin datasets shows its ability to extract complex signaling patterns. Our versatile and easy-to-use toolkit CellChat and a web-based Explorer ( http://www.cellchat.org/ ) will help discover novel intercellular communications and build cell-cell communication atlases in diverse tissues.

Deciphering cell-cell interactions and communication from gene expression

Cell-cell interactions orchestrate organismal development, homeostasis and single-cell functions. When cells do not properly interact or improperly decode molecular messages, disease ensues. Thus, the identification and quantification of intercellular signalling pathways has become a common analysis performed across diverse disciplines. The expansion of protein-protein interaction databases and recent advances in RNA sequencing technologies have enabled routine analyses of intercellular signalling from gene expression measurements of bulk and single-cell data sets. In particular, ligand-receptor pairs can be used to infer intercellular communication from the coordinated expression of their cognate genes. In this Review, we highlight discoveries enabled by analyses of cell-cell interactions from transcriptomic data and review the methods and tools used in this context.

Understanding cell-cell communication and signaling in the colorectal cancer microenvironment

Carcinomas are complex heterocellular systems containing epithelial cancer cells, stromal fibroblasts, and multiple immune cell-types. Cell-cell communication between these tumor microenvironments (TME) and cells drives cancer progression and influences response to existing therapies. In order to provide better treatments for patients, we must understand how various cell-types collaborate within the TME to drive cancer and consider the multiple signals present between and within different cancer types. To investigate how tissues function, we need a model to measure both how signals are transferred between cells and how that information is processed within cells. The interplay of collaboration between different cell-types requires cell-cell communication. This article aims to review the current in vitro and in vivo mono-cellular and multi-cellular cultures models of colorectal cancer (CRC), and to explore how they can be used for single-cell multi-omics approaches for isolating multiple types of molecules from a single-cell required for cell-cell communication to distinguish cancer cells from normal cells. Integrating the existing single-cell signaling measurements and models, and through understanding the cell identity and how different cell types communicate, will help predict drug sensitivities in tumor cells and between- and within-patients responses.

Inference of Intercellular Communications and Multilayer Gene-Regulations of Epithelial-Mesenchymal Transition From Single-Cell Transcriptomic Data

Epithelial-to-mesenchymal transition (EMT) plays an important role in many biological processes during development and cancer. The advent of single-cell transcriptome sequencing techniques allows the dissection of dynamical details underlying EMT with unprecedented resolution. Despite several single-cell data analysis on EMT, how cell communicates and regulates dynamics along the EMT trajectory remains elusive. Using single-cell transcriptomic datasets, here we infer the cell-cell communications and the multilayer gene-gene regulation networks to analyze and visualize the complex cellular crosstalk and the underlying gene regulatory dynamics along EMT. Combining with trajectory analysis, our approach reveals the existence of multiple intermediate cell states (ICSs) with hybrid epithelial and mesenchymal features. Analyses on the time-series datasets from cancer cell lines with different inducing factors show that the induced EMTs are context-specific: the EMT induced by transforming growth factor B1 (TGFB1) is synchronous, whereas the EMTs induced by epidermal growth factor and tumor necrosis factor are asynchronous, and the responses of TGF-β pathway in terms of gene expression regulations are heterogeneous under different treatments or among various cell states. Meanwhile, network topology analysis suggests that the ICSs during EMT serve as the signaling in cellular communication under different conditions. Interestingly, our analysis of a mouse skin squamous cell carcinoma dataset also suggests regardless of the significant discrepancy in concrete genes between in vitro and in vivo EMT systems, the ICSs play dominant role in the TGF-β signaling crosstalk. Overall, our approach reveals the multiscale mechanisms coupling cell-cell communications and gene-gene regulations responsible for complex cell-state transitions.

Oligonucleotide Probe for Transcriptome in Vivo Analysis (TIVA) of Single Neurons with Minimal Background

Messenger RNA (mRNA) isolated from single cells can generate powerful biological insights, including the discovery of new cell types with unique functions as well as markers potentially predicting a cell's response to various therapeutic agents. We previously introduced an oligonucleotide-based technique for site-selective, photoinduced biotinylation and capture of mRNA within a living cell called transcriptome in vivo analysis (TIVA). Successful application of the TIVA technique hinges upon its oligonucleotide probe remaining completely inert (or "caged") to mRNA unless photoactivated. To improve the reliability of TIVA probe caging in diverse and challenging biological conditions, we applied a rational design process involving iterative modifications to the oligonucleotide construct. In this work, we discuss these design motivations and present an optimized probe with minimal background binding to mRNA prior to photolysis. We assess its caging performance through multiple in vitro assays including FRET analysis, native gel comigration, and pull down with model mRNA transcripts. Finally, we demonstrate that this improved probe can also isolate mRNA from single living neurons in brain tissue slices with excellent caging control.

Divergent Resistance Mechanisms to Immunotherapy Explain Responses in Different Skin Cancers

The advent of immune checkpoint therapy for metastatic skin cancer has greatly improved patient survival. However, most skin cancer patients are refractory to checkpoint therapy, and furthermore, the intra-immune cell signaling driving response to checkpoint therapy remains uncharacterized. When comparing the immune transcriptome in the tumor microenvironment of melanoma and basal cell carcinoma (BCC), we found that the presence of memory B cells and macrophages negatively correlate in both cancers when stratifying patients by their response, with memory B cells more present in responders. Moreover, inhibitory immune signaling mostly decreases in melanoma responders and increases in BCC responders. We further explored the relationships between macrophages, B cells and response to checkpoint therapy by developing a stochastic differential equation model which qualitatively agrees with the data analysis. Our model predicts BCC to be more refractory to checkpoint therapy than melanoma and predicts the best qualitative ratio of memory B cells and macrophages for successful treatment.

Predicting cell-to-cell communication networks using NATMI

Development of high throughput single-cell sequencing technologies has made it cost-effective to profile thousands of cells from diverse samples containing multiple cell types. To study how these different cell types work together, here we develop NATMI (Network Analysis Toolkit for Multicellular Interactions). NATMI uses connectomeDB2020 (a database of 2293 manually curated ligand-receptor pairs with literature support) to predict and visualise cell-to-cell communication networks from single-cell (or bulk) expression data. Using multiple published single-cell datasets we demonstrate how NATMI can be used to identify (i) the cell-type pairs that are communicating the most (or most specifically) within a network, (ii) the most active (or specific) ligand-receptor pairs active within a network, (iii) putative highly-communicating cellular communities and (iv) differences in intercellular communication when profiling given cell types under different conditions. Furthermore, analysis of the Tabula Muris (organism-wide) atlas confirms our previous prediction that autocrine signalling is a major feature of cell-to-cell communication networks, while also revealing that hundreds of ligands and their cognate receptors are co-expressed in individual cells suggesting a substantial potential for self-signalling.

Single cell expression data allows for inferring cell-cell communication between cells expressing ligands and those expressing their cognate receptors. Here the authors present an updated and curated database of ligand-receptor pairs and a Python-based toolkit to construct and analyse communication networks from single cell and bulk expression data.

Single cell transcriptomics of human epidermis identifies basal stem cell transition states

How stem cells give rise to epidermis is unclear despite the crucial role the epidermis plays in barrier and appendage formation. Here we use single cell-RNA sequencing to interrogate basal stem cell heterogeneity of human interfollicular epidermis and find four spatially distinct stem cell populations at the top and bottom of rete ridges and transitional positions between the basal and suprabasal epidermal layers. Cell-cell communication modeling suggests that basal cell populations serve as crucial signaling hubs to maintain epidermal communication. Combining pseudotime, RNA velocity, and cellular entropy analyses point to a hierarchical differentiation lineage supporting multi-stem cell interfollicular epidermal homeostasis models and suggest that transitional basal stem cells are stable states essential for proper stratification. Finally, alterations in differentially expressed transitional basal stem cell genes result in severe thinning of human skin equivalents, validating their essential role in epidermal homeostasis and reinforcing the critical nature of basal stem cell heterogeneity.

Visualizing Single-Cell RNA-seq Data with Semisupervised Principal Component Analysis

Single-cell RNA-seq (scRNA-seq) is a powerful tool for analyzing heterogeneous and functionally diverse cell population. Visualizing scRNA-seq data can help us effectively extract meaningful biological information and identify novel cell subtypes. Currently, the most popular methods for scRNA-seq visualization are principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE). While PCA is an unsupervised dimension reduction technique, t-SNE incorporates cluster information into pairwise probability, and then maximizes the Kullback-Leibler divergence. Uniform Manifold Approximation and Projection (UMAP) is another recently developed visualization method similar to t-SNE. However, one limitation with UMAP and t-SNE is that they can only capture the local structure of the data, the global structure of the data is not faithfully preserved. In this manuscript, we propose a semisupervised principal component analysis (ssPCA) approach for scRNA-seq visualization. The proposed approach incorporates cluster-labels into dimension reduction and discovers principal components that maximize both data variance and cluster dependence. ssPCA must have cluster-labels as its input. Therefore, it is most useful for visualizing clusters from a scRNA-seq clustering software. Our experiments with simulation and real scRNA-seq data demonstrate that ssPCA is able to preserve both local and global structures of the data, and uncover the transition and progressions in the data, if they exist. In addition, ssPCA is convex and has a global optimal solution. It is also robust and computationally efficient, making it viable for scRNA-seq cluster visualization.

Random Parametric Perturbations of Gene Regulatory Circuit Uncover State Transitions in Cell Cycle

Many biological processes involve precise cellular state transitions controlled by complex gene regulation. Here, we use budding yeast cell cycle as a model system and explore how a gene regulatory circuit encodes essential information of state transitions. We present a generalized random circuit perturbation method for circuits containing heterogeneous regulation types and its usage to analyze both steady and oscillatory states from an ensemble of circuit models with random kinetic parameters. The stable steady states form robust clusters with a circular structure that are associated with cell cycle phases. This circular structure in the clusters is consistent with single-cell RNA sequencing data. The oscillatory states specify the irreversible state transitions along cell cycle progression. Furthermore, we identify possible mechanisms to understand the irreversible state transitions from the steady states. We expect this approach to be robust and generally applicable to unbiasedly predict dynamical transitions of a gene regulatory circuit.

Toward Systems Biomarkers of Response to Immune Checkpoint Blockers

Immunotherapy with checkpoint blockers (ICBs), aimed at unleashing the immune response toward tumor cells, has shown a great improvement in overall patient survival compared to standard therapy, but only in a subset of patients. While a number of recent studies have significantly improved our understanding of mechanisms playing an important role in the tumor microenvironment (TME), we still have an incomplete view of how the TME works as a whole. This hampers our ability to effectively predict the large heterogeneity of patients' response to ICBs. Systems approaches could overcome this limitation by adopting a holistic perspective to analyze the complexity of tumors. In this Mini Review, we focus on how an integrative view of the increasingly available multi-omics experimental data and computational approaches enables the definition of new systems-based predictive biomarkers. In particular, we will focus on three facets of the TME toward the definition of new systems biomarkers. First, we will review how different types of immune cells influence the efficacy of ICBs, not only in terms of their quantification, but also considering their localization and functional state. Second, we will focus on how different cells in the TME interact, analyzing how inter- and intra-cellular networks play an important role in shaping the immune response and are responsible for resistance to immunotherapy. Finally, we will describe the potential of looking at these networks as dynamic systems and how mathematical models can be used to study the rewiring of the complex interactions taking place in the TME.

Uncovering an Organ's Molecular Architecture at Single-Cell Resolution by Spatially Resolved Transcriptomics

Revealing fine-scale cellular heterogeneity among spatial context and the functional and structural foundations of tissue architecture is fundamental within biological research and pharmacology. Unlike traditional approaches involving single molecules or bulk omics, cutting-edge, spatially resolved transcriptomics techniques offer near-single-cell or even subcellular resolution within tissues. Massive information across higher dimensions along with position-coordinating labels can better map the whole 3D transcriptional landscape of tissues. In this review, we focus on developments and strategies in spatially resolved transcriptomics, compare the cell and gene throughput and spatial resolution in detail for existing methods, and highlight the enormous potential in biomedical research.

New avenues for systematically inferring cell-cell communication: through single-cell transcriptomics data

For multicellular organisms, cell-cell communication is essential to numerous biological processes. Drawing upon the latest development of single-cell RNA-sequencing (scRNA-seq), high-resolution transcriptomic data have deepened our understanding of cellular phenotype heterogeneity and composition of complex tissues, which enables systematic cell-cell communication studies at a single-cell level. We first summarize a common workflow of cell-cell communication study using scRNA-seq data, which often includes data preparation, construction of communication networks, and result validation. Two common strategies taken to uncover cell-cell communications are reviewed, e.g., physically vicinal structure-based and ligand-receptor interaction-based one. To conclude, challenges and current applications of cell-cell communication studies at a single-cell resolution are discussed in details and future perspectives are proposed.

scRCMF: Identification of Cell Subpopulations and Transition States From Single-Cell Transcriptomes

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the "fate decisions" of transition cells within the transition. The present study provides new insights into transition events during differentiation and development.

Inferring spatial and signaling relationships between cells from single cell transcriptomic data

Single-cell RNA sequencing (scRNA-seq) provides details for individual cells; however, crucial spatial information is often lost. We present SpaOTsc, a method relying on structured optimal transport to recover spatial properties of scRNA-seq data by utilizing spatial measurements of a relatively small number of genes. A spatial metric for individual cells in scRNA-seq data is first established based on a map connecting it with the spatial measurements. The cell-cell communications are then obtained by "optimally transporting" signal senders to target signal receivers in space. Using partial information decomposition, we next compute the intercellular gene-gene information flow to estimate the spatial regulations between genes across cells. Four datasets are employed for cross-validation of spatial gene expression prediction and comparison to known cell-cell communications. SpaOTsc has broader applications, both in integrating non-spatial single-cell measurements with spatial data, and directly in spatial single-cell transcriptomics data to reconstruct spatial cellular dynamics in tissues.

Environmental Carcinogenesis at the Single-Cell Level

Elucidating the mechanisms behind how exposure to environmental chemicals can lead to cancer is not easy due to the complex natures of these compounds and the challenges to establish biologically relevant experimental models to study them. Environmental chemicals often present selective mechanisms of action on different cell types and can be involved in the modulation of targeted cells and their microenvironment, including immune cells. Currently, the limitations of traditional epidemiologic correlation analyses, in vitro cell-based assays, and animal models are that they are unable to comprehensively examine cellular heterogeneity and the tissue-selective influences. To this end, we propose utilizing single-cell RNA-sequencing (scRNA-seq) to more effectively capture the subtle and complex effects of environmental chemicals and how their exposure could lead to cancer. scRNA-seq's capabilities for studying gene expression level data at a significantly higher resolution relative to bulk RNA-sequencing (RNA-seq) enable studies to evaluate how environmental chemicals regulate gene transcription on different cell types as well as how these compounds impact signaling pathways and interactions between cells in the tissue microenvironment. These studies will be valuable for evaluating environmental chemicals' carcinogenic properties at the individual cell level.See all articles in this CEBP Focus section, "Environmental Carcinogenesis: Pathways to Prevention."

A powerful and flexible weighted distance-based method incorporating interactions between DNA methylation and environmental factors on health outcomes

MOTIVATION

Deoxyribonucleic acid (DNA) methylation plays a crucial role in human health. Studies have demonstrated associations between DNA methylation and environmental factors with evidence also supporting the idea that DNA methylation may modify the risk of environmental factors on health outcomes. However, due to high dimensionality and low study power, current studies usually focus on finding differential methylation on health outcomes at CpG level or gene level combining multiple CpGs and/or finding environmental effects on health outcomes but ignoring their interactions on health outcomes. Here we introduce the idea of a pseudo-data matrix constructed with cross-product terms between CpGs and environmental factors that are able to capture their interactions. We then develop a powerful and flexible weighted distance-based method with the pseudo-data matrix where association strength was used as weights on CpGs, environmental factors and their interactions to up-weight signals and down-weight noises in distance calculations.

RESULTS

We compared the power of this novel approach and several comparison methods in simulated datasets and the Mothers and Newborns birth cohort of the Columbia Center for Children's Environmental Health to determine whether prenatal polycyclic aromatic hydrocarbons interacts with DNA methylation in association with Attention Deficit Hyperactivity Disorder and Mental Development Index at age 3.

AVAILABILITY AND IMPLEMENTATION

An R code for the proposed method Dw-M-E-int together with a tutorial and a sample dataset is available for downloading from http://www.columbia.edu/∼sw2206/softwares.htm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-cell reconstruction of differentiation trajectory reveals a critical role of ETS1 in human cardiac lineage commitment

BACKGROUND

Cardiac differentiation from human pluripotent stem cells provides a unique opportunity to study human heart development in vitro and offers a potential cell source for cardiac regeneration. Compared to the large body of studies investigating cardiac maturation and cardiomyocyte subtype-specific induction, molecular events underlying cardiac lineage commitment from pluripotent stem cells at early stage remain poorly characterized.

RESULTS

In order to uncover key molecular events and regulators controlling cardiac lineage commitment from a pluripotent state during differentiation, we performed single-cell RNA-Seq sequencing and obtained high-quality data for 6879 cells collected from 6 stages during cardiac differentiation from human embryonic stem cells and identified multiple cell subpopulations with distinct molecular features. Through constructing developmental trajectory of cardiac differentiation and putative ligand-receptor interactions, we revealed crosstalk between cardiac progenitor cells and endoderm cells, which could potentially provide a cellular microenvironment supporting cardiac lineage commitment at day 5. In addition, computational analyses of single-cell RNA-Seq data unveiled ETS1 (ETS Proto-Oncogene 1) activation as an important downstream event induced by crosstalk between cardiac progenitor cells and endoderm cells. Consistent with the findings from single-cell analysis, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq) against ETS1 revealed genomic occupancy of ETS1 at cardiac structural genes at day 9 and day 14, whereas ETS1 depletion dramatically compromised cardiac differentiation.

CONCLUSION

Together, our study not only characterized the molecular features of different cell types and identified ETS1 as a crucial factor induced by cell-cell crosstalk contributing to cardiac lineage commitment from a pluripotent state, but may also have important implications for understanding human heart development at early embryonic stage, as well as directed manipulation of cardiac differentiation in regenerative medicine.

scGEAToolbox: a Matlab toolbox for single-cell RNA sequencing data analysis

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) technology has revolutionized the way research is done in biomedical sciences. It provides an unprecedented level of resolution across individual cells for studying cell heterogeneity and gene expression variability. Analyzing scRNA-seq data is challenging though, due to the sparsity and high dimensionality of the data.

RESULTS

I developed scGEAToolbox-a Matlab toolbox for scRNA-seq data analysis. It contains a comprehensive set of functions for data normalization, feature selection, batch correction, imputation, cell clustering, trajectory/pseudotime analysis, and network construction, which can be combined and integrated to building custom workflow. While most of the functions are implemented in native Matlab, wrapper functions are provided to allow users to call the "third-party" tools developed in Matlab or other languages. Furthermore, scGEAToolbox is equipped with sophisticated graphical user interfaces (GUIs) generated with App Designer, making it an easy-to-use application for quick data processing.

AVAILABILITY

https://github.com/jamesjcai/scGEAToolbox.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.