The CXCL12/CXCR7 signaling axis, isoforms, circadian rhythms, and tumor cellular composition dictate gradients in tissue

Chemokine CXCL12 gradients drive chemotaxis in a CXCR4-dependent mechanism and have been implicated in cancer metastasis. While CXCL12 gradients are typically studied in organized, defined environments, the tumor microenvironment is disorganized. In vivo, CXCL12 gradients depend on many factors: the number and arrangement of cells secreting and degrading CXCL12, isoform-dependent binding to the extracellular matrix, diffusion, and circadian fluctuations. We developed a computational model of the tumor microenvironment to simulate CXCL12 gradient dynamics in disorganized tissue. There are four major findings from the model. First, CXCL12-β and -γ form higher magnitude (steeper) gradients compared to CXCL12-α. Second, endothelial CXCR7+ cells regulate CXCL12 gradient direction by controlling concentrations near but not far from the vasculature. Third, the magnitude and direction of CXCL12 gradients are dependent on the local composition of secreting and scavenging cells within the tumor. We theorize that “micro-regions” of cellular heterogeneity within the tumor are responsible for forming strong gradients directed into the blood. Fourth, CXCL12 circadian fluctuations influence gradient magnitude but not direction. Our simulations provide predictions for future experiments in animal models. Understanding the generation of CXCL12 gradients is crucial to inhibiting cancer metastasis.