Cell-Cell Communication Inference (Ligand-Target)
Detect interactions between ligands and target cell types
1 datasets · 14 methods · 2 control methods · 2 metrics
The growing availability of single-cell data has sparked an increased interest in the inference of cell-cell communication (CCC), with an ever-growing number of computational tools developed for this purpose.
Different tools propose distinct preprocessing steps with diverse scoring functions, that are challenging to compare and evaluate. Furthermore, each tool typically comes with its own set of prior knowledge. To harmonize these, Dimitrov et al, 2022 recently developed the LIANA framework, which was used as a foundation for this task.
The challenges in evaluating the tools are further exacerbated by the lack of a gold standard to benchmark the performance of CCC methods. In an attempt to address this, Dimitrov et al use alternative data modalities, including the spatial proximity of cell types and downstream cytokine activities, to generate an inferred ground truth. However, these modalities are only approximations of biological reality and come with their own assumptions and limitations. In time, the inclusion of more datasets with known ground truth interactions will become available, from which the limitations and advantages of the different CCC methods will be better understood.
This subtask evaluates the methods’ ability to predict interactions, the corresponding of cytokines of which, are inferred to be active in the target cell types. This subtask focuses on the prediction of interactions from steady-state, or single-context, single-cell data.
Summary
Metrics
- Precision-recall AUC (Davis and Goadrich 2006): Area under the precision-recall curve for the binary classification task predicting interactions.
- Odds Ratio (Bland 2000): The odds ratio represents the ratio of true and false positives within a set of prioritized interactions (top ranked hits) versus the same ratio for the remainder of the interactions. Thus, in this scenario odds ratios quantify the strength of association between the ability of methods to prioritize interactions and those interactions assigned to the positive class.
Results
Results table of the scores per method, dataset and metric (after scaling). Use the filters to make a custom subselection of methods and datasets. The “Overall mean” dataset is the mean value across all datasets.
Details
Methods
- CellPhoneDB (max) (Efremova et al. 2020): CellPhoneDBv2 calculates a mean of ligand-receptor expression as a measure of interaction magnitude, along with a permutation-based p-value as a measure of specificity. Here, we use the former to prioritize interactions, subsequent to filtering according to p-value less than 0.05. Links: Docs.
- CellPhoneDB (sum) (Efremova et al. 2020): CellPhoneDBv2 calculates a mean of ligand-receptor expression as a measure of interaction magnitude, along with a permutation-based p-value as a measure of specificity. Here, we use the former to prioritize interactions, subsequent to filtering according to p-value less than 0.05. Links: Docs.
- Connectome (max) (Raredon et al. 2022): Connectome uses the product of ligand-receptor expression as a measure of magnitude, and the average of the z-transformed expression of ligand and receptor as a measure of specificity. Links: Docs.
- Connectome (sum) (Raredon et al. 2022): Connectome uses the product of ligand-receptor expression as a measure of magnitude, and the average of the z-transformed expression of ligand and receptor as a measure of specificity. Links: Docs.
- Log2FC (max) (Dimitrov et al. 2022): logFC (implemented in LIANA and inspired by iTALK) combines both expression and magnitude, and represents the average of one-versus-the-rest log2-fold change of ligand and receptor expression per cell type. Links: Docs.
- Log2FC (sum) (Dimitrov et al. 2022): logFC (implemented in LIANA and inspired by iTALK) combines both expression and magnitude, and represents the average of one-versus-the-rest log2-fold change of ligand and receptor expression per cell type. Links: Docs.
- Magnitude Rank Aggregate (max) (Dimitrov et al. 2022): RobustRankAggregate generates a consensus rank of all methods implemented in LIANA providing either specificity or magnitude scores. Links: Docs.
- Magnitude Rank Aggregate (sum) (Dimitrov et al. 2022): RobustRankAggregate generates a consensus rank of all methods implemented in LIANA providing either specificity or magnitude scores. Links: Docs.
- NATMI (max) (Hou et al. 2020): NATMI uses the product of ligand-receptor expression as a measure of magnitude. As a measure of specificity, NATMI proposes specificity.edge = \frac{l}{l_s} \cdot \frac{r}{r_s}; where l and r represent the average expression of ligand and receptor per cell type, and l_s and r_s represent the sums of the average ligand and receptor expression across all cell types. We use its specificity measure, as recommended by the authors for single-context predictions. Links: Docs.
- NATMI (sum) (Hou et al. 2020): NATMI uses the product of ligand-receptor expression as a measure of magnitude. As a measure of specificity, NATMI proposes specificity.edge = \frac{l}{l_s} \cdot \frac{r}{r_s}; where l and r represent the average expression of ligand and receptor per cell type, and l_s and r_s represent the sums of the average ligand and receptor expression across all cell types. We use its specificity measure, as recommended by the authors for single-context predictions. Links: Docs.
- Random Events (Open Problems for Single Cell Analysis Consortium 2022): Random generation of cell-cell communication events by random selection of ligand, receptor, source, target, and score. Links: Docs.
- SingleCellSignalR (max) (Cabello-Aguilar et al. 2020): SingleCellSignalR provides a magnitude score as LRscore = \frac{\sqrt{lr}}{\mu+\sqrt{lr}}; where l and r are the average ligand and receptor expression per cell type, and \mu is the mean of the expression matrix. Links: Docs.
- SingleCellSignalR (sum) (Cabello-Aguilar et al. 2020): SingleCellSignalR provides a magnitude score as LRscore = \frac{\sqrt{lr}}{\mu+\sqrt{lr}}; where l and r are the average ligand and receptor expression per cell type, and \mu is the mean of the expression matrix. Links: Docs.
- Specificity Rank Aggregate (max) (Dimitrov et al. 2022): RobustRankAggregate generates a consensus rank of all methods implemented in LIANA providing either specificity or magnitude scores. Links: Docs.
- Specificity Rank Aggregate (sum) (Dimitrov et al. 2022): RobustRankAggregate generates a consensus rank of all methods implemented in LIANA providing either specificity or magnitude scores. Links: Docs.
- True Events (Open Problems for Single Cell Analysis Consortium 2022): Perfect prediction of cell-cell communication events from target data. Links: Docs.
Baseline methods
- Random Events: Random generation of cell-cell communication events by random selection of ligand, receptor, source, target, and score.
- True Events: Perfect prediction of cell-cell communication events from target data.
Datasets
- Triple negative breast cancer atlas (Wu et al. 2021): Human breast cancer atlas (Wu et al., 2021), with cytokine activities, inferred using a multivariate linear model with cytokine-focused signatures, as assumed true cell-cell communication (Dimitrov et al., 2022). 42512 cells x 28078 features with 29 cell types from 10 patients.
Quality control results
✓ All checks succeeded!