iSCanGuide Help

Cell–Cell Communication

This page provides an overview of how to perform cell–cell communication analysis in iSCanGuide.

iSCanGuide identifies signaling interactions between cell types by determining which cell populations act as signal senders and receivers, and by identifying the associated ligand–receptor pairs. For example, in a thymus sample, the platform can detect signaling from immune cells to stromal cells through the IL1A → IL1R1 pathway and rank this interaction relative to other candidate signaling events.

For spatial transcriptomics datasets (e.g., Visium or Xenium), iSCanGuide also visualizes the spatial distribution of signaling activity within the tissue.

The analysis integrates multiple published scoring approaches and combines their rankings to generate a consensus result.

When to use this analysis

Use cell–cell communication analysis when you want to:

  • Identify which cell types are communicating with one another.

  • Obtain a ranked list of the most significant ligand–receptor interactions.

  • For spatial datasets, determine where signaling occurs within the tissue and distinguish spatially organized interactions from diffuse signaling patterns.

Before you start

This analysis requires a sample with pre-assigned cell type labels, provided either through a sample metadata column or a saved Cell Annotation.

To access the Spatial Co-expression and Inflow tabs, the sample type must be set to Spatial. Non-spatial datasets will include only the Cell Type Pairs and Interactions tabs.

Cell-Cell Communication tab in the bottom Data and Analysis panel

To open the Cell–Cell Communication panel, select the Cell–Cell Communication tab in the Data and Analysis Panel located at the bottom of the Analysis page.

Workflow

New Analysis

Cell-Cell Communication panel showing New Analysis tab

The Cell-Cell Communication panel has two tabs: Existing Results and New Analysis.

Click the New Analysis tab to configure a run.

Form fields

New Analysis form

Mandatory fields are indicated with a red asterisk.

  • Sample (required): Select the processed sample to analyze. Only pre-processed samples are available for selection, and each analysis run supports a single sample.

  • Cell type labels from (required): Specify the source of the cell type assignments.

    • Sample metadata — Cell type labels are retrieved from a column in the sample's cell metadata (the AnnData obs table). Use this option if cell type annotations were included during sample import (for example, a cell_type column).

    • Cell Annotation — Cell type labels are retrieved from an annotation saved within iSCanGuide, such as results generated by SingleR or manual annotations.

    Selecting an appropriate label source is important, as inaccurate or low-quality annotations may lead to weak or difficult-to-interpret results.

  • Cell Annotation or Metadata column (required): Select the specific annotation or metadata column containing the cell type labels.

  • Ligand–receptor resource: Choose the reference database of known ligand–receptor pairs used for scoring interactions. Default: Consensus.

Ligand-receptor resources

Resource

Notes

Consensus

Union of several public databases. Default for most analyses.

Mouse Consensus

Use for mouse samples.

CellCall

Specialty database.

CellChatDB

Use to compare to a CellChat-based publication.

CellPhoneDB

Use to compare to a CellPhoneDB-based publication.

CellTalkDB

Specialty database.

ConnectomeDB 2020

Specialty database.

EMBRACE

Specialty database.

ICELLNET

Specialty database.

LRDB

Specialty database.

Ramilowski 2015

Specialty database.

Advanced Parameters

Advanced Parameters section expanded

Click the > button next to Advanced Parameters to expand the additional settings. For most analyses, the default values are recommended.

  • Expression proportion (default: 0.10): Defines the minimum fraction of cells within a cell type that must express a gene for that gene to be considered expressed in the population. Ligands or receptors below this threshold are excluded for that cell type, and the corresponding ligand–receptor interaction is not scored.

    • Increase the threshold (e.g., 0.25) to apply stricter filtering and reduce noise.

    • Decrease the threshold (e.g., 0.05) when working with rare cell populations or sparse datasets.

  • Min cells per type (default: 5): Cell types containing fewer cells than this threshold are excluded from the analysis. This helps avoid unstable scores driven by very small populations.

    • Increase the threshold (e.g., 2050) to exclude very small clusters that may represent doublets or low-quality populations.

  • Permutations (default: 1000): Number of random permutations used to calculate p-values for non-spatial analyses. Increasing the number of permutations improves p-value precision but also increases runtime.

    The default value of 1000 is appropriate for most analyses. This parameter applies only to non-spatial datasets; spatial analyses use a different statistical approach and ignore this setting.

  • Spatial bandwidth (default: Auto-detect; spatial datasets only): Controls the size of the spatial neighbourhood (Gaussian kernel) used for Moran's I and Inflow score calculations, measured in coordinate units.

    In most cases, the recommended setting is Auto-detect, which selects an appropriate value based on the median nearest-neighbour distance within the dataset.

    • Smaller values emphasize fine-scale local signaling patterns.

    • Larger values capture broader spatial gradients.

  • Methods (default: all five): Specifies the individual scoring methods combined to generate the consensus ranking. Keeping all five methods enabled is recommended for the most robust and stable results.

Available methods:

Method

What it scores

CellPhoneDB

Mean expression of cognate L–R partners; permutation-based specificity

Connectome

Weighted product of scaled ligand and receptor expression

log2FC

Fold-change of ligand/receptor in sender vs. receiver populations

NATMI

Edge weights from expression product, specificity from contribution to global network

SingleCellSignalR

LRscore (regularised geometric mean)

Geometric Mean

Geometric mean of ligand and receptor expression

scSeqComm

scSeqComm score

CellChat

CellChat score

Submitting

Simple workflow after submitting

Click the Run Analysis button. The job runs in the background; when it finishes it appears under Existing Results. Visit Study Logs to monitor job progress.

Existing Results

Existing Results table

All completed analyses for the current study are displayed in the Existing Results tab. Each row includes the following information:

  • Sample: The analyzed sample.

  • Label source: The source of the cell type annotations.

  • Resource: The ligand–receptor reference database used for the analysis.

  • Spatial: Indicates whether the sample is spatially resolved.

  • # Cell types: Number of cell types included in the analysis.

  • # Interactions: Total number of inferred ligand–receptor interactions.

  • Created: Date and time the analysis was generated.

Click the + button on a row to expand it and view the analysis parameters.

Expanded result row showing parameters

Click View in the Actions column to open the results. Click Delete to remove the analysis.

CCC Results

The results page shows a summary at the top: Sample, Type, Resource, # Cell Types, # Interactions, # Cells, Methods. Results are organised across up to four tabs depending on the sample type.

Cell Type Pairs tab

Cell Type Pairs tab with chord diagram and table

This tab gives a macro view of cell–cell communication.

  • Chord diagram: A bipartite diagram with ribbons connecting sender and receiver cell types. Hovering over a ribbon shows the number of interactions, mean distance, and interacting status.

  • Table (right side): Every directed sender → receiver pair with:

    • # Interactions: Number of significant L–R pairs for this sender → receiver direction.

    • Mean Distance (spatial samples only): Average on-tissue distance between the two populations.

    • Interacting: Whether the two populations are co-localised on tissue (spatial samples only).

Click a row in the table to jump to the Interactions tab with that sender → receiver pair pre-selected.

Exporting

Chord diagram export panel

Hover over the chord diagram to reveal export options. Click the export button to download the chart at multiple resolutions.

Table export button

Click the export button above the table to download the results as a CSV.

Interactions tab

Interactions tab showing the ligand-receptor table

This tab shows the full ranked list of ligand–receptor pairs. Use it to find the strongest, most specific signals.

Key columns

  • Sender/Receiver: The cell types participating in the interaction, representing the signaling and receiving populations.

  • Ligand/Receptor: The ligand–receptor gene complexes involved in the interaction.

  • Magnitude/Magnitude rank: Measures the expression-based strength of the interaction. Higher magnitude values indicate stronger signaling activity, while lower magnitude ranks indicate interactions that are stronger relative to all other tested pairs.

  • Specificity/Specificity rank: Reflects the permutation-based significance of the interaction within a specific sender → receiver context. Lower specificity ranks indicate interactions that are more specific to that particular cell-type pair.

Filtering and Sorting

  • Use the column filter icons on Sender, Receiver, Ligand, or Receptor to focus on specific interactions.

  • Adjust the shared filters at the top of the table. The default thresholds are:

    • Magnitude rank ≤ 0.05

    • Specificity rank ≤ 0.05

    If no interactions are returned, consider relaxing the thresholds to ≤ 0.1.

  • Use the method dropdown on the right to display scores from a specific analysis method:

    • CellPhoneDB

    • Connectome

    • log2FC

    • NATMI

    • SingleCellSignalR

  • Sort results by:

    • Magnitude Significance — lower values indicate more statistically credible interactions.

    • Magnitude — higher values indicate stronger inferred signaling.

Click the export button to download the table as a CSV file.

Spatial Co-expression tab

Only available for spatial samples.

Spatial Co-expression tab

This tab evaluates whether significant ligand–receptor interactions are spatially co-localized within the tissue.

Columns

  • Moran's I: Measures the spatial autocorrelation of the ligand × receptor expression product. Higher values indicate stronger spatial clustering. Default filter: Moran's I ≥ 0.01.

  • Moran's p-value: Permutation-based p-value associated with Moran's I.

  • Mean Cosine: Measures the cosine similarity between ligand and receptor expression patterns across neighbouring spots or cells. Higher values indicate more similar spatial expression patterns.

  • Std Dev: Standard deviation of the cosine similarity scores.

Scatter Plot

The scatter plot on the right displays:

  • x-axis: Moran's I

  • y-axis: Mean Cosine

Interactions appearing in the upper-right region of the plot are both strongly spatially clustered and highly co-expressed.

Filtering

Enable Use interactions tab filters to apply the filtering criteria from the Interactions tab. This is recommended to restrict results to ligand–receptor pairs that are both statistically significant and spatially co-expressed.

Export Co-expression Pairs to Features

Select rows using the checkboxes, enter a name, and click Save export to push the pairs to the Features explorer in the left panel.

Each exported pair becomes a per-cell feature with value equal to the cosine similarity score for that ligand–receptor pair at each cell. Feature label format: ligand^receptor (e.g., IL1A^IL1R1).

Inflow tab

Only available for spatial samples where at least one L–R pair has Moran's I > 0.01.

Inflow tab

This tab quantifies how much signalling each receiver cell type receives from each ligand–receptor pair, accounting for spatial structure.

The table is hierarchical:

  • Parent rows (sender + L–R pair): Spatial autocorrelation stats for each inflow column (Moran's I and FDR-corrected p-value).

  • Child rows (sender + L–R pair + specific receiver): The inflow magnitude (Inflow Mean) and significance (Inflow p-value) broken down by receiver cell type.

Default filters: Moran's I ≥ 0.01 and FDR p-value ≤ 0.05.

The violin plot on the right shows the distribution of Moran's I per sender — the dashed line marks α = 0.01. Use it to compare how spatially organised each sender's output is.

The violin plot is interactive:

  • Click a sender in the violin plot to filter the view to that sender's L–R pairs and switch the chart to a bar plot ranking those pairs by Moran's I.

  • Click a bar in the bar plot to select and highlight the corresponding row in the table on the left.

Export Inflow Pairs to Features

Select rows using the checkboxes, enter a name, and click Save export to push pairs to the Features explorer. Two row types can be exported:

  • Parent row (unmasked): per-cell value is the inflow score for that sender–ligand–receptor axis across all receiver cells. Label format: Sender: ligand^receptor (e.g., Immune: IL1A^IL1R1).

  • Child row (receiver-masked): same inflow score scoped to a specific receiver cell type. Label format: Sender: ligand^receptor → ReceiverCellType (e.g., Immune: IL1A^IL1R1 → Fibroblast).

Visualising a pair on tissue

  1. From the Spatial Co-expression or Inflow tab, select the rows you want and Save export to the Features tab.

  2. Close the results panel. The exported feature (e.g., LCK^CD8A_CD8B) appears in the left panel under Cell-Cell Communication.

  3. Click the Select Label button next to the feature to add it to the plot grid.

  4. Use Side by Side blend mode to keep cluster, annotation, and L–R co-expression panels in the same coordinate system.

Quick reference

Column

What it means

# Interactions

Count of significant L–R pairs for a sender → receiver direction.

Mean Distance

Average on-tissue distance between the sender and receiver populations.

Magnitude

Score from the chosen method; higher = stronger inferred signal.

Specificity

Permutation p-value; lower = more specific to this pair.

Magnitude rank / Specificity rank

Per-method ranks normalised to [0, 1]; used by the consensus aggregator.

Moran's I

Spatial autocorrelation of ligand × receptor (or of the inflow score). Higher = more spatially clustered.

FDR p-value

Multiple-testing-corrected p-value for spatial tests.

Mean Cosine

Cosine similarity of ligand and receptor patterns in a spatial neighbourhood.

Inflow Mean / Inflow p-value

Magnitude and significance of incoming signal at the receiver.

Troubleshooting

  • Spatial Co-expression or Inflow tabs are missing — The sample is not designated as a spatial dataset. Re-import the sample with the sample type set to Spatial.

  • Interactions table is empty or contains very few interactions — Try relaxing the rank filters (for example, increasing thresholds to ≤ 0.1) or selecting an individual scoring method instead of the consensus ranking.

    If the issue persists, verify that the dataset uses HGNC gene symbols. All ligand–receptor resources expect standard gene symbols (e.g. IL1A, CD8A) rather than Ensembl IDs or other identifier formats. Check the Feature ID Column specified during sample import under Data Management. If the dataset was imported using Ensembl IDs, re-import the sample using the appropriate gene symbol column.

  • "No cell types pass the threshold" — This typically indicates that:

    • the cell type annotations are too granular, resulting in very small groups, or

    • the Min cells per type threshold is set too high.

    Lower the threshold or use a broader cell type annotation.

  • Inflow tab is missing for a spatial sample — No ligand–receptor interactions exceeded the Moran's I > 0.01 threshold required for spatial clustering. This suggests that signaling patterns in the dataset are relatively diffuse rather than spatially localized.

  • Most signaling appears in an "unassigned" category — The current annotation may be too coarse or incomplete. Re-run the analysis using a more detailed or refined cell type annotation source.

  • Spatial bandwidth warning — Switch from Auto-detect to a manually specified bandwidth value. For Visium datasets, values corresponding to approximately 1–3 spot radii are often appropriate.

  • Long runtimes on large datasets — Reduce the number of permutations or perform the analysis on a subset of relevant cell types before scaling to the full dataset.

Last modified: 21 May 2026