ULMnet: inferring physical cell-cell communication networks from scRNAseq data using univariate linear models

TL;DR

The paper introduces ULMnet, a computational method that uses scRNAseq data to infer physical cell-cell communication networks by identifying multiplets, which represent physically attached cells.

Contribution

The novel contribution is the development of ULMnet, a method using univariate linear models to infer physical cell-cell interactions from scRNAseq data.

Findings

ULMnet achieved ~56% sensitivity and ~99% precision in predicting FACS-sorted doublets with known constituents.

Applying ULMnet to partially dissociated tissues revealed physical networks matching microanatomical structures.

The method validated biologically plausible interactions using spatial transcriptomics data.

Abstract

Cells in tissues interact by direct physical contact or over short and long distances via secreted mediators. Cell-cell communication inference has now become routine in downstream scRNAseq analysis but this mostly fails to capture physical cell-cell interactions due to tissue dissociation. Multiplets (mostly doublets) in scRNAseq may represent undissociated physically attached cells that become sequenced together. Hence, identifying multiplets may serve as a good starting point to harness scRNAseq data for physical cell-cell interaction inference.

In this study, we develop a computational method, called ULMnet, which utilizes univariate linear models to identify multiplets in scRNAseq datasets, predict their cellular compositions, and infer physical cell-cell interaction networks. ULMnet was applied to multiple datasets ranging from cell pairs, partially dissociated tissues to conventional scRNAseq datasets of healthy and cancer tissues.

Indeed, our method showed good sensitivity (~56%) with excellent precision (~99%) in predicting FACS sorted doublet cell pairs with known constituents (ground truth), recording comparable or superior performance over two other existing methods. Also, applying it to scRNAseq data of partially dissociated tissues containing real multiplets unraveled physical networks which recapitulated the microanatomical structures of the tested tissues. This further underscores the accuracy in our predictions to capture biologically meaningful interactions. Finally, we tested our method on classical scRNAseq datasets and obtained biologically reasonable results. For example, when tested on classical cancer scRNAseq datasets, we recovered important interactions which followed biologically plausible cell interactions, validated by cell-cell colocalization in matched spatial transcriptomics datasets. This reassured the accuracy of our method in depicting physical interactions only between cells that were truly in close proximity in tissues.

Figures