Inter-cellular communication with stromal cells is vital for cancer cells. Molecules involved in the communication are potential drug targets. To identify them systematically, we applied a systems level analysis that combined reverse network engineering with causal effect estimation. Using only observational transcriptome profiles we searched for paracrine factors sending messages from activated hepatic stellate cells (HSC) to hepatocellular carcinoma (HCC) cells. We condensed these messages to predict ten proteins that, acting in concert, cause the majority of the gene expression changes observed in HCC cells. Among the 10 paracrine factors were both known and unknown cancer promoting stromal factors, the former including Placental Growth Factor (PGF) and Periostin (POSTN), while Pregnancy-Associated Plasma Protein A (PAPPA) was among the latter. Further support for the predicted effect of PAPPA on HCC cells came from both in vitro studies that showed PAPPA to contribute to the activation of NFB signaling, and clinical data, which linked higher expression levels of PAPPA to advanced stage HCC. In summary, this study demonstrates the potential of causal modeling in combination with a condensation step borrowed from gene set analysis [Model-based Gene Set Analysis (MGSA)] in the identification of stromal signaling molecules influencing the cancer phenotype.