Lymphovascular invasion and histologic grade are associated with specific genomic profiles in invasive carcinomas of the breast

Lymphovascular invasion (LVI) and histologic grade are clinical parameters of high prognostic value in breast cancer and indicate the level of tumor aggressiveness. Many studies have focused on the association of breast cancer subtypes with gene expression and chromosomal profiles, but considerably less genomic information is available regarding traditional prognostic factors such as histologic grade and LVI. We studied by array-CGH a group of 57 invasive ductal carcinomas of the breast to outline the DNA copy number aberration (CNA) profile linked to high histologic grades and LVI. Selected CNAs were validated using real-time quantitative PCR (qPCR). Furthermore, gene expression analysis was performed in a subset of 32 of these tumors, and findings were integrated with array-CGH data. Our findings indicated an accumulation of genomic alterations in high-grade breast tumors compared to low-grade samples. Grade III tumors showed higher number of CNAs and larger aberrations than low-grade tumors and displayed a wide range of chromosomal aberrations, which were mainly 5p, 8q, 10p, 17q12, and 19 gains, and 3p, 4, 5q proximal, 9p, 11p, 18q, and 21 losses. The presence of LVI, a well-established prognostic marker, was not significantly associated with increased genomic instability in comparison to breast tumors negative for LVI, considering the total number of chromosomal alterations. However, a slightly increase in the frequency of specific alterations could be detected in LVI-positive group, such as gains at 5p, 16p, 17q12, and 19, and losses at 8p, 11q, 18q, and 21. Three newly reported small-scale rearrangements were detected in high-risk tumors (LVI-positive grade III) harboring putative breast cancer genes (amplicons at 4q13.3 and 11p11.2, and a deletion at 12p12.3). Furthermore, gene expression analysis uncovered networks highlighting S100A8, MMP1, and MED1 as promising candidate genes involved in high-grade and LVI-positive tumors. In summary, a group of genomic regions could be associated with high-risk tumors, and expression analysis pinpointed candidate genes deserving further investigation. The data has shed some light on the molecular players involved in two highly relevant prognostic factors and may further add to the understanding of the mechanisms of breast cancer aggressiveness. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s13277-014-2786-z) contains supplementary material, which is available to authorized users.