Crosstalk between tumor cells and other cells within the tumor microenvironment (TME) plays a crucial role in tumor progression, metastases, and therapy resistance. We present iTALK, a computational approach to characterize and illustrate intercellular communication signals in the multicellular tumor ecosystem using single-cell RNA sequencing data. iTALK can in principle be used to dissect the complexity, diversity, and dynamics of cell-cell communication from a wide range of cellular processes.The TME has emerged as a key modulator of tumor progression, immune evasion, and emergence of the anti-tumor therapy resistance mechanisms 1, 2 . The TME includes a diversity of cell types such as tumor cells, a heterogeneous group of immune cells, and the nonimmune stromal components. Tumor cells orchestrate and interact dynamically with these non-tumor components, and the crosstalk between them is thought to provide key signals that can direct and promote tumor cell growth and migration. Through this intercellular communication, tumor cells can elicit profound phenotypic changes in other TME cells such as tumor-associated fibroblasts, macrophages and T cells, and reprogram the TME, in order to escape from immune surveillance to facilitate survival. Therefore, a better understanding of the cell-cell communication signals may help identify novel modulating therapeutic strategies for better patient advantage. However, this has been hampered by the lack of bioinformatics tools for efficient data analysis and visualization.Here, we present iTALK (identifying and illustrating alterations in intercellular signaling network; https://github.com/Coolgenome/iTALK), an open source R package designed to profile and visualize the ligand-receptor mediated intercellular cross-talk signals from singlecell RNA sequencing data (scRNA-seq) ( Fig. 1 and Online Methods). We demonstrated that iTALK can be successfully applied to scRNA-seq data to capture highly abundant ligandreceptor gene (or transcript) pairs, identify gains or losses of cellular interactions by comparative analysis, and track the dynamic changes of intercellular communication signals in longitudinal samples. Notably, functional annotation of ligand-receptor genes is automatically added with our curated iTALK ligand-receptor database, and the output can be visualized in different formats with our efficient data visualization tool, which is implemented as part of iTALK. This approach can be applied to data sets ranging from hundreds to hundreds of thousands of cells and is not limited by sequencing platforms. It is also noteworthy that, in addition to studying the TME, iTALK can also be applied to a wide range of biomedical research fields that involve cell-cell communication.