Abstract
Cellular cross-talk in tissue microenvironments is fundamental to normal and pathological biological processes. Global assessment of cell-cell interactions (CCIs) is not yet technically feasible, but computational efforts to reconstruct these interactions have been proposed. Current computational approaches that identify CCI often make the simplifying assumption that pairwise interactions are independent of one another, which can lead to reduced accuracy. We present REMI (REgularized Microenvironment Interactome), a graph-based algorithm that predicts ligand-receptor (LR) interactions by accounting for LR dependencies on high-dimensional, small-sample size datasets. We apply REMI to reconstruct the human lung adenocarcinoma (LUAD) interactome from a bulk flow-sorted RNA sequencing dataset, then leverage single-cell transcriptomics data to increase the cell type resolution and identify LR prognostic signatures among tumor-stroma-immune subpopulations. We experimentally confirmed colocalization of CTGF:LRP6 among malignant cell subtypes as an interaction predicted to be associated with LUAD progression. Our work presents a computational approach to reconstruct interactomes and identify clinically relevant CCIs.