Intratumoral lymphatic endothelial cell infiltration reflecting lymphangiogenesis is counterbalanced by immune responses and better cancer biology in the breast cancer tumor microenvironment

Lymphangiogenesis, the generation of new lymphatic vessels from existing ones, results from the dynamic interactions of lymphatic endothelial cells and the tumor microenvironment (TME). It is well known that lymphangiogenesis occurs during the initial stage of metastasis in various types of malignant tumors. However, it is currently not used as a biomarker partially because gold standard method to quantify it is labor and cost intensive. We hypothesized that the quantity of intratumoral lymphatic endothelial cells (iLECs) in the TME is an indicator of lymphangiogenesis and a predictor of metastatic potential and overall survival in breast cancer. We analyzed a total of 4145 breast cancer patients from the Cancer Genome Atlas (TCGA) and GSE96058 by quantifying their iLECs using the xCell algorithm, and correlated these scores with patient survival, tumor grade, and cancer stage. We also assessed various pro- and anti-cancer gene sets for each tumor to characterize tumor behavior and aggressiveness. As we expected, high-iLEC breast cancer demonstrated enriched lymphoangiogenesis and angiogenesis gene sets and was associated with increased expressions of related genes. Also enriched were inflammatory response and immune response-related gene sets; IL2/STAT5 pathway, IL6/JAK/STAT3 pathway, TNF alpha pathway, allograft rejection, and complement as well as cancer stemness related gene sets like Notch signaling, Hedgehog signaling, epithelial mesenchymal transition, and Wnt beta-catenin signaling. Tumors with high-iLEC showed higher proportions of stromal cells and fewer anti-cancer immune cells. On the other hand, iLEC score did not correlate with patient survival or lymph node metastasis. Surprisingly, breast cancers with fewer iLECs demonstrated enriched E2F Targets, G2M Checkpoint, MYC Targets v1, and MTORC1 signaling which are cancer cell proliferation-related gene sets and exhibited an abundance of pro-cancer immune cells. The amount of iLEC correlated inversely with Ki67 expression and histological grade, which is in agreement that low-iLEC breast cancer was associated with enhanced cancer cell proliferation. In conclusion, while iLECs can be used as a surrogate for lymphangiogenesis in breast cancer, low-iLEC tumors also exhibit features which correspond to aggressive tumor biology, which may explain why the amount of iLECs was not associated with patient survival in our cohorts.

Automated summary

Lymphangiogenesis plays a crucial role in the occurrence and development of breast cancer. High-iLEC breast cancer is associated with enriched lymphangiogenesis, angiogenesis, inflammatory response, and immune response-related gene sets, along with more stromal cells and fewer anti-cancer immune cells. However, the quantity of iLECs is not related to patient survival or lymph node metastasis. Low-iLEC breast cancer, on the other hand, is associated with cancer cell proliferation and an abundance of pro-cancer immune cells.