Lymphovascular invasion in breast cancer is associated with gene expression signatures of cell proliferation but not lymphangiogenesis or immune response.

Abstract

While the prognostic relevance of lymphovascular invasion (LVI) in breast cancer is well known, its molecular biology is poorly understood. We hypothesized that pathologically determined LVI reflects molecular features of tumors and can be discerned from their genomic and transcriptomic profiles. LVI status and Nottingham histological scores of primary breast tumors of The Cancer Genome Atlas (TCGA) project were assessed from pathology reports; other clinical and molecular data were obtained from TCGA data portals and publications. Two independent datasets (GSE5460 and GSE7849) were combined and used for validation. LVI status was determinable for 639 and 196 cases of the TCGA and validation cohorts, among whom LVI incidence was 37.8% and 37.2%, respectively. LVI was associated with high tumor Ki67 expression, advanced pathologic stage, and high Nottingham scores. LVI-positive cases had worse overall and progression-free survival regardless of cancer subtype. Surprisingly, in both cohorts, LVI was not associated with lymphangiogenesis or lymphatic vessel density as estimated from tumor expression of lymphatic endothelium-associated genes. LVI-positive tumors had higher genome copy number aberrations, aneuploidy, and homologous recombination defects, but not single-nucleotide variations or intra-tumor genome heterogeneity. Tumor immune cell composition and cytolytic activity was not associated with LVI status. On the other hand, expression of cell proliferation-related genes was significantly increased in LVI-positive tumors. Our study suggests that breast cancer with LVI is a highly proliferative cancer, and it does not correlate with gene expression markers for lymphangiogenesis or immune response.