Analysis of Single-Cell RNA-Seq Identifies Cell-Cell Communication Associated with Tumor Characteristics

Abstract

Tumor ecosystems are composed of multiple cell types that communicate by ligand-receptor interactions. Targeting ligand-receptor interactions (for instance, with immune checkpoint inhibitors) can provide significant benefits for patients. However, our knowledge of which interactions occur in a tumor and how these interactions affect outcome is still limited. We present an approach to characterize communication by ligand-receptor interactions across all cell types in a microenvironment using single-cell RNA sequencing. We apply this approach to identify and compare the ligand-receptor interactions present in six syngeneic mouse tumor models. To identify interactions potentially associated with outcome, we regress interactions against phenotypic measurements of tumor growth rate. In addition, we quantify ligand-receptor interactions between T cell subsets and their relation to immune infiltration using a publicly available human melanoma dataset. Overall, this approach provides a tool for studying cell-cell interactions, their variability across tumors, and their relationship to outcome.

Keywords: cancer patient samples; cell-cell communication; computational analysis; ligand-receptor interaction; single-cell RNA sequencing; syngeneic mouse models; tumor microenvironment.

Conflict of interest statement

DECLARATION OF INTERESTS

Several co-authors are employed by Merrimack Pharmaceuticals, Inc.

Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.

Figures

Figure 1.. T-SNE Visualization of Single-Cell Sequencing Data and Cell Type Classification

(A) The percentage of cells positive for a variety of immune cell markers as measured by either scRNA-seq (x axis) or flow cytometry (y axis) is consistent across tumor models.(B and C) t-distributed scholastic neighbor embedding (t-SNE) plots of cells from six syngeneic tumor models show distinct clusters predominantly determined by cell type.(B) Cells are colored by the tumor model from which the cell originated.(C) Cells are colored by the cell type label assigned at the end of the classification procedure.(D) Percentages of cell types vary across the different tumor models. To reflect the actual cell type abundances, only data from samples not enriched for CD45 are shown. See also Figures S1 and S2.

Figure 2.. Quantification of Cell-Cell Interactions Occurring in the Tumor Microenvironment

Heatmaps show selected interaction scores calculated as the product of the average ligand expression of the first cell type and average receptor expression of the second cell type. Cell type labels are written as (cell type expressing the ligand) ‒ (cell type expressing the receptor). Black dots indicate interactions that are significantly present across all tumor (one-sided Wilcoxon rank-sum test and Benjamini Hochberg false discovery rate [FDR]

Figure 3.. Interaction Scores Correlate with Relevant Characteristics of the Tumor Microenvironment

(A) Tumor volume (y axis) of treatment-naive mice measured over time (x axis) (Table S4). Dashed lines indicate the mean of a syngeneic tumor model, and shaded areas represent 1 SEM (n = 8 for Sa1N, 7 for LL2, 10 for CT26, 9 for EMT6, and 9 for MC38). Instances with no shading result from only one mouse surviving at the measured time points. Linear curves were fit to the log-normalized growth curves, and the slope of fit curves was used as a metric for tumor growth. (B) Quantified growth rates for each model. Each point represents a single mouse, and the horizontal black line indicates the median growth rate used for correlation with interaction scores.(C) Heatmap showing the Spearman correlation of interaction scores (shown in Figure 2) with tumor growth. Interactions marked with black circles indicate correlations with p

Figure 4.. Assessing Cell-Cell Interactions Occurring in Human Metastatic Melanoma

(A) Cell-cell interactions involving Tregs in human metastatic melanoma averaged across 19 tumor samples. The cell type labels are written as (cell type expressing the ligand) ‒ (cell type expressing the receptor). Black dots indicate interactions that are significantly present across all tumors (one-sided Wilcoxon rank-sum test and Benjamini Hochberg FDR