Background: Lymphovascular invasion (LVI), the passage of carcinoma cells through lymphatic and blood vessels, is an important early step in metastasis; however, LVI is excluded from most breast cancer (BC) clinical risk assessments. Previous studies assessed the prognostic value of LVI to estimate clinical outcomes. To gain understanding of the molecular basis of LVI, we evaluated differentially expressed genes (DEGs) between tumors with LVI versus those without LVI, stratified by the 70-gene signature (MammaPrint/MP) and 80-gene molecular subtyping signature (BluePrint/BP). Methods: The prospective, observational FLEX Study (NCT03053193) includes stage I-III BC patients who receive MP/BP testing and consent to full transcriptome and clinical data collection. Patients with LVI (n=581) and without LVI (n=600, randomly selected), enrolled from 2017 to present, were included. LVI was assessed by local pathology laboratories. Differential gene expression analysis of 44k Agilent microarray data was performed with R limma package. DEGs were compared within all samples, BP Luminal subtype, MP risk groups (Low Risk [LR]/Luminal A and High Risk [HR]/Luminal B), and by lymph node (LN) status. DEGs with FDR<0.05 were considered significant. Results: Of tumors with LVI (LVI+), 66% were MP HR; notably, 51% of tumors without LVI (LVI-) were MP HR. LVI was associated with larger T stage, LN involvement, high grade, negative ER status by IHC, and younger patient age (LVI+ vs. LVI-, p<0.05 for all comparisons). Patient ethnicity, obesity, and tumor type did not differ by LVI status; however, prevalence of type 2 diabetes trended higher in patients with LVI+ HR tumors (21%), compared with LVI- HR (15%, p=0.09) and LVI+ LR (11%, p=0.004). There were significant transcriptomic differences between LVI+ and LVI, with most DEGs evident in the Luminal B subset. DEGs in LVI+, LN-negative (LN-) tumors overlapped substantially with the overall Luminal group analysis. Functional enrichment analysis showed dysregulation of cell cycle, extracellular matrix (ECM) organization, cell adhesion, and cytokine receptor pathways. Gene sets related to insulin growth factor pathways were also enriched in LVI+ tumors. Conclusions: DEGs associated with LVI were primarily found in MP HR Luminal, LN-negative tumors; enrichment analysis suggested dysregulation of ECM organization and cell adhesion pathways, consistent with previous reports. DEGs were not associated with LVI presence in LN+ tumors, suggesting that LVI assessment may be less relevant in LN+ breast cancer. Future studies will assess clinical outcomes, as well as LVI-associated gene expression in BP Basal- and HER2-type tumors. However, the current analysis indicates few DEGs in LVI+ MP LR tumors; thus, the potential prognostic information gained from LVI-associated gene expression is likely already captured by the MP and BP signatures. Clinical trial information: NCT03053193