Sunitinib is an antiangiogenic therapy given as a first-line treatment for renal cell carcinoma (RCC). While treatment improves progression-free survival, most patients relapse. We hypothesized that patient relapse can stem from the development of a lymphatic network driven by the production of the main growth factor for lymphatic endothelial cells, VEGFC. In this study, we found that sunitinib can stimulate gene transcription and increase VEGFC mRNA half-life. In addition, sunitinib activated p38 MAPK, which resulted in the upregulation/activity of HuR and inactivation of tristetraprolin, two AU-rich element-binding proteins. Sunitinib stimulated a VEGFC-dependent development of lymphatic vessels in experimental tumors. This may explain our findings of increased lymph node invasion and new metastatic sites in 30% of sunitinib-treated patients and increased lymphatic vessels found in 70% of neoadjuvant treated patients. In summary, a therapy dedicated to destroying tumor blood vessels induced the development of lymphatic vessels, which may have contributed to the treatment failure..

The biology of the metastatic colonization process remains a poorly understood phenomenon. To improve our knowledge of its dynamics, we conducted a modelling study based on multi-modal data from an orthotopic murine experimental system of metastatic renal cell carcinoma. The standard theory of metastatic colonization usually assumes that secondary tumours, once established at a distant site, grow independently from each other and from the primary tumour. Using a mathematical model that translates this assumption into equations, we challenged this theory against our data that included: 1) dynamics of primary tumour cells in the kidney and metastatic cells in the lungs, retrieved by green fluorescent protein tracking, and 2) magnetic resonance images (MRI) informing on the number and size of macroscopic lesions. Critically, when calibrated on the growth of the primary tumour and total metastatic burden, the predicted theoretical size distributions were not in agreement with the MRI observations. Moreover, tumour expansion only based on proliferation was not able to explain the volume increase of the metastatic lesions. These findings strongly suggested rejection of the standard theory, demonstrating that the time development of the size distribution of metastases could not be explained by independent growth of metastatic foci. This led us to investigate the effect of spatial interactions between merging metastatic tumours on the dynamics of the global metastatic burden. We derived a mathematical model of spatial tumour growth, confronted it with experimental data of single metastatic tumour growth, and used it to provide insights on the dynamics of multiple tumours growing in close vicinity. Together, our results have implications for theories of the metastatic process and suggest that global dynamics of metastasis development is dependent on spatial interactions between metastatic lesions.