Spatial cancer systems biology resolves heterotypic interactions and identifies disruption of spatial hierarchy as a pathological driver event.

Spatially annotated single-cell datasets provide unprecedented opportunities to dissect cell-cell communication in development and disease. Heterotypic signaling includes interactions between different cell types and is well established in tissue development and spatial organization. Epithelial organization requires several different programs that are tightly regulated. Planar cell polarity (PCP) is the organization of epithelial cells along the planar axis, orthogonal to the apical-basal axis. Here, we investigate PCP factors and explore the implications of developmental regulators as malignant drivers. Utilizing cancer systems biology analysis, we derive a gene expression network for WNT-ligands (WNT) and their cognate frizzled (FZD) receptors in skin cutaneous melanoma. The profiles supported by unsupervised clustering of multiple-sequence alignments identify ligand-independent signaling and implications for metastatic progression based on the underpinning developmental spatial program. Omics studies and spatial biology connect developmental programs with oncological events and explain key spatial features of metastatic aggressiveness. Dysregulation of prominent PCP factors such as specific representatives of the WNT and FZD families in malignant melanoma recapitulates the development program of normal melanocytes but in an uncontrolled and disorganized fashion.