Back

Modelling multicellular coordination by bridging cell-cell communication and intracellular regulation through multilayer networks

Trimbour, R.; Ramirez Flores, R. O.; Saez-Rodriguez, J.; Cantini, L.

2026-01-22 bioinformatics
10.64898/2026.01.20.700561 bioRxiv
Show abstract

In multicellular organisms, cells with various roles and locations coordinate to provide systemic and cohesive response to perturbations. These complex behaviors emerge from a complex interplay between intracellular regulation and intercellular signals that mediate cell-cell communication. While single-cell technologies opened the possibility of studying both, most methods focus solely on one of these aspects. Thus, they are only able to partially recover in vivo and multicellular behaviors. We here introduce ReCoN (REconstruction of multicellular COordination Networks from single-cell data), a framework combining intracellular gene regulation and cell-cell communication to provide insights into multicellular coordination from single-cell data. First, ReCoN infers from single-cell data a heterogeneous multilayer network containing both cell-type-specific intracellular subnetworks and ligand-receptor interactions. Through random walk with restart explorations, ReCoN then infers the response of each cell type to both intra- and extracellular perturbations, such as a gene knock-out or a cytokine, respectively. ReCoN was evaluated on predicting the in vivo response of immune cell-types to different cytokines and on recovering cardiac cell-type response in heart failure. It highlighted the role of indirect effects, where cells emit secondary messengers in response to the initial perturbation to coordinate multicellular transcriptomic responses. Additionally, ReCoN predicted distinct fibroblast states emerging in different microenvironments reconstructed from spatial data. ReCoN provides an interpretable modeling framework for multicellular systems that allows for the simulation of perturbations, including the assessment of the cellular selectivity of these treatments in vivo. Ultimately, it can help design patient-specific molecular therapies. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=107 SRC="FIGDIR/small/700561v1_ufig1.gif" ALT="Figure 1"> View larger version (34K): org.highwire.dtl.DTLVardef@7cdceforg.highwire.dtl.DTLVardef@18cfa3eorg.highwire.dtl.DTLVardef@812645org.highwire.dtl.DTLVardef@f7bce5_HPS_FORMAT_FIGEXP M_FIG C_FIG

Matching journals

The top 3 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.