Aggressive Phenotype of Cells Disseminated via Hematogenous and Lymphatic Route in Breast Cancer Patients

Intratumoral heterogeneity of breast cancer remains a major challenge in successful treatment. Failure of cancer therapies can also be accredited to inability to systemically eradicate cancer stem cells (CSCs). Recent evidence points to the role of epithelial-mesenchymal transition (EMT) in expanding the pool of tumor cells with CSCs features. Thus, we assessed expression level as well as heterogeneity of CSCs markers in primary tumors (PT), lymph node metastasis (LNM), and circulating tumor cells (CTCs)–enriched blood fractions in order to correlate them with signs of EMT activation as well as clinicopathological data of breast cancer patients. Level of CSCs markers (ALDH1, CD44, CD133, OCT-4, NANOG) and EMT markers was quantified in PT (N=107), LNM (N=56), and CTCs-enriched blood fractions (N=85). Heterogeneity of CSCs markers expression within each PT and LNM was assessed by calculating Gini Index. Percentage of ALDH1-positive cells was elevated in PT in comparison to LNM (P = .005). However, heterogeneity of the four CSCs markers: ALDH1 (P = .019), CD133 (P = .009), OCT-4 (P = .027), and CD44 (P < .001) was decreased in LNM. Samples classified as mesenchymal (post-EMT) showed elevated expression of CSCs markers (OCT-4 and CD44 in PT; OCT-4 in LNM; ALDH1, OCT-4, NANOG, CD44 in CTCs). Patients with mesenchymal-like CTCs had worse prognosis than patients with epithelial-like or no CTCs (P = .0025). CSCs markers are enriched in PT, LNM, and CTCs with mesenchymal features, but their heterogeneity is decreased in metastatic lymph nodes. Mesenchymal CTCs phenotype correlates with poor prognosis of the patients.     

G Protein Coupled Receptor Kinase 3 Regulates Breast Cancer Migration,Invasion, and Metastasis

Triple negative breast cancer (TNBC) is a heterogeneous disease that has a poor prognosis and limited treatment options. Chemokine receptor interactions are important modulators of breast cancer metastasis; however, it is now recognized that quantitative surface expression of one important chemokine receptor, CXCR4, may not directly correlate with metastasis and that its functional activity in breast cancer may better inform tumor pathogenicity. G protein coupled receptor kinase 3 (GRK3) is a negative regulator of CXCR4 activity, and we show that GRK expression correlates with tumorigenicity, molecular subtype, and metastatic potential in human tumor microarray analysis. Using established human breast cancer cell lines and an immunocompetent in vivo mouse model, we further demonstrate that alterations in GRK3 expression levels in tumor cells directly affect migration and invasion in vitro and the establishment of distant metastasis in vivo. The effects of GRK3 modulation appear to be specific to chemokine-mediated migration behaviors without influencing tumor cell proliferation or survival. These data demonstrate that GRK3 dysregulation may play an important part in TNBC metastasis.     

Genetic variants associated with breast cancer risk for Ashkenazi Jewish women with strong family histories but no identifiable BRCA1/2 mutation

The ability to establish genetic risk models is critical for early identification and optimal treatment of breast cancer. For such a model to gain clinical utility, more variants must be identified beyond those discovered in previous genome-wide association studies (GWAS). This is especially true for women at high risk because of family history, but without BRCA1/2 mutations. This study incorporates three datasets in a GWAS analysis of women with Ashkenazi Jewish (AJ) homogeneous ancestry. Two independent discovery cohorts comprised 239 and 238 AJ women with invasive breast cancer or preinvasive ductal carcinoma in situ and strong family histories of breast cancer, but lacking the three BRCA1/2 founder mutations, along with 294 and 230 AJ controls, respectively. An independent, third cohort of 203 AJ cases with familial breast cancer history and 263 healthy controls of AJ women was used for validation. A total of 19 SNPs were identified as associated with familial breast cancer risk in AJ women. Among these SNPs, 13 were identified from a panel of 109 discovery SNPs, including an FGFR2 haplotype. In addition, six previously identified breast cancer GWAS SNPs were confirmed in this population. Seven of the 19 markers were significant in a multivariate predictive model of familial breast cancer in AJ women, three novel SNPs [rs17663555(5q13.2), rs566164(6q21), and rs11075884(16q22.2)], the FGFR2 haplotype, and three previously published SNPs [rs13387042(2q35), rs2046210(ESR1), and rs3112612(TOX3)], yielding moderate predictive power with an area under the curve (AUC) of the ROC (receiver-operator characteristic curve) of 0.74. Population-specific genetic variants in addition to variants shared with populations of European ancestry may improve breast cancer risk prediction among AJ women from high-risk families without founder BRCA1/2 mutations.     

Significance of Lymph Node Metastasis in Cancer Dissemination of Head and Neck Cancer

Lymph node metastasis (LNM) in many solid cancers is a well-known prognostic factor; however, it has been debated whether regional LNM simply reflects tumor aggressiveness or is a source for further tumor dissemination. Similarly, the metastatic process in head and neck cancer (HNC) has not been fully evaluated. Thus, we aimed to investigate the relative significance of LNM in metastatic cascade of HNC using functional imaging of HNC patients and molecular imaging in in vivo models. First, we analyzed ¹⁸Fluorodeoxyglucose positron emission tomography (PET) parameters of 117 patients with oral cancer. The primary tumor and nodal PET parameters were measured separately, and survival analyses were conducted on the basis of clinical and PET variables to identify significant prognostic factors. In multivariate analyses, we found that only the metastatic node PET values were significant. Next, we compared the relative frequency of lung metastasis in primary ear tumors versus lymph node (LN) tumors, and we tested the rate of lung metastasis in another animal model, in which each animal had both primary and LN tumors that were expressing different colors. As a result, LN tumors showed higher frequencies of lung metastasis compared to orthotopic primary tumors. In color-matched comparisons, the relative contribution to lung metastasis was higher in LN tumors than in primary tumors, although both primary and LN tumors caused lung metastases. In summary, tumors growing in the LN microenvironment spread to systemic sites more commonly than primary tumors in HNC, suggesting that the adequate management of LNM can reduce further systemic metastasis.     

Peritumoral lymphangiogenesis induced by vascular endothelial growth factor C and D promotes lymph node metastasis in breast cancer patients

ABSTRACT: BACKGROUND: Mounting clinical and experimental data suggest that the migration of tumor cells into lymph nodes is greatly facilitated by lymphangiogenesis. Vascular endothelial growth factor (VEGF)-C and D have been identified as lymphangiogenic growth factors and play an important role in tumor lymphangiogenesis. The purpose of this study was to investigate the location of lymphangiogenesis driven by tumor-derived VEGF-C/D in breast cancer, and to determine the role of intratumoral and peritumoral lymphatic vessel density (LVD) in lymphangiogenesis in breast cancer. METHODS: The expression levels of VEGF-C/D were determined by immunohistochemistry, and intratumoral LVD and peritumoral LVD were assessed using immunohistochemistry and the D2-40 antibody in 73 patients with primary breast cancer. The associations of intratumoral LVD and peritumoral LVD with VEGF-C/D expression, clinicopathological features and prognosis were assessed. RESULTS: VEGF-C and D expression were significantly higher in breast cancer than benign disease (P < 0.01). VEGF-C (P < 0.001) and VEGF-D (P = 0.005) expression were significantly associated with peritumoral LVD, but not intratumoral LVD. Intratumoral LVD was associated with tumor size (P = 0.01). Peritumoral LVD was significantly associated with lymph node metastasis (LNM; P = 0.005), lymphatic vessel invasion (LVI; P = 0.017) and late tumor,node,metastasis(TNM) stage (P = 0.011). Moreover, peritumoral LVD was an independent risk factor for axillary lymph node metastasis, overall survival and disease-free survival in multivariate analysis. CONCLUSIONS: This study suggests that tumor-derived VEGF-C/D induce peritumoral lymphangiogenesis, which may be one mechanism that leads to lymphatic invasion and metastatic spread. Peritumoral LVD has potential as an independent prognostic factor in breast cancer patients.     

Breast cancer susceptibility loci in association with age at menarche,age at natural menopause and the reproductive lifespan

Background: Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with breast cancer risk. Some of these loci have unknown functional significance and may mediate the effects of hormonal exposures on breast cancer risk. We examined relationships between breast cancer susceptibility variants and menstrual/reproductive factors using data from two population-based studies. Methods: The first analysis was based on a sample of 1328 women age 20–74 who participated as controls in a case–control study of breast cancer conducted in three U.S. states. We evaluated the associations between age at menarche, age at natural menopause and the reproductive lifespan with 13 previously identified breast cancer variants. Associations were also examined with a genetic score created as the sum of at-risk alleles across the 13 variants. For validation, significant results were evaluated in a second dataset comprised 1353 women age 43–86 recruited as part of a cohort study in Wisconsin. Results: Neither the genetic score nor any of the 13 variants considered individually were associated with age at menarche or reproductive lifespan. Two SNPs were associated with age at natural menopause; every increase in the minor allele (A) of rs17468277 (CASP8) was associated with a 1.12 year decrease in menopause age (p = 0.02). The minor allele (G) of rs10941679 (5p12) was associated with a 1.01 year increase in age at natural menopause (p = 0.01). The results were not replicated in the validation cohort (B = −0.61, p = 0.14 and B = −0.01, p = .0.98, respectively). Conclusions: The evaluated variants and reproductive experiences may work through separate pathways to influence breast cancer risk.     

Shift work,circadian gene variants and risk of breast cancer

Circadian (clock) genes have been linked with several functions relevant to cancer, and epidemiologic research has suggested relationships with breast cancer risk for variants in NPAS2, CLOCK, CRY2 and TIMELESS. Increased breast cancer risk has also been observed among shift workers, suggesting potential interactions in relationships of circadian genes with breast cancer. Relationships with breast cancer of 100 SNPs in 14 clock-related genes, as well as potential interactions with shift work history, were investigated in a case–control study (1042 cases, 1051 controls). Odds ratios in an additive genetic model for European-ancestry participants (645 cases, 806 controls) were calculated, using a two-step correction for multiple testing: within each gene through permutation testing (10,000 permutations), and correcting for the false discovery rate across genes. Interactions of genotypes with ethnicity and shift work (<2 years vs ≥2 years) were evaluated individually. Following permutation analysis, two SNPs (rs3816360 in ARNTL and rs11113179 in CRY1) displayed significant associations with breast cancer and one SNP (rs3027188 in PER1) was marginally significant; however, none were significant following adjustment for the false discovery rate. No significant interaction with shift work history was detected. If shift work causes circadian disruption, this was not reflected in associations between clock gene variants and breast cancer risk in this study. Larger studies are needed to assess interactions with longer durations (>30 years) of shift work that have been associated with breast cancer.     

Type II Transmembrane Serine Protease Gene Variants Associate with Breast Cancer

Type II transmembrane serine proteases (TTSPs) are related to tumor growth, invasion, and metastasis in cancer. Genetic variants in these genes may alter their function, leading to cancer onset and progression, and affect patient outcome. Here, 464 breast cancer cases and 370 controls were genotyped for 82 single-nucleotide polymorphisms covering eight genes. Association of the genotypes was estimated against breast cancer risk, breast cancer–specific survival, and survival in different treatment groups, and clinicopathological variables. SNPs in TMPRSS3 (rs3814903 and rs11203200), TMPRSS7 (rs1844925), and HGF (rs5745752) associated significantly with breast cancer risk (P _trend = 0.008–0.042). SNPs in TMPRSS1 (rs12151195 and rs12461158), TMPRSS2 (rs2276205), TMPRSS3 (rs3814903), and TMPRSS7 (rs2399403) associated with prognosis (P = 0.004–0.046). When estimating the combined effect of the variants, the risk of breast cancer was higher with 4–5 alleles present compared to 0–2 alleles (P = 0.0001; OR, 2.34; 95% CI, 1.39–3.94). Women with 6–8 survival-associating alleles had a 3.3 times higher risk of dying of breast cancer compared to women with 1–3 alleles (P = 0.001; HR, 3.30; 95% CI, 1.58–6.88). The results demonstrate the combined effect of variants in TTSPs and their related genes in breast cancer risk and patient outcome. Functional analysis of these variants will lead to further understanding of this gene family, which may improve individualized risk estimation and development of new strategies for treatment of breast cancer.     

Association of glutamate carboxypeptidase II (GCPII) haplotypes with breast and prostate cancer risk

In view of the pivotal role of glutamate carboxypeptidase II (GCPII) in carcinogenesis, its expression as prostate specific membrane antigen (PSMA) and folate hydrolase (FOLH1) may be influenced by its haplotypes, contributing to the etiology of prostate and breast cancer. To test this hypothesis, breast and prostate cancer cases and controls were subjected to whole gene screening of GCPII and correlated with plasma folate levels and PSMA expression. The impact of variants on a 3-dimensional structure of GCPII was explored by in silico studies. Six novel variations i.e. V108A, P160S, Y176H, D191V, G206R and G245S; and two known variations i.e. R190W and H475Y were identified in GCPII. All-wild haplotype and a haplotype harbouring D191V showed association with breast cancer risk while haplotypes harbouring V108A and P160S reduced the risk. Haplotypes with V108A and G245S variants showed increased risk for prostate cancer due to high PSMA expression while P160S conferred protection against prostate cancer. In silico studies suggests that P160S and R190W variants result in relaxed substrate binding facilitating either rapid catalysis or exchange of substrates and products in the active site which was substantiated by high plasma folate levels associated with these variants. On the contrary, D191V was associated with very low plasma folate levels despite having a high PSMA expression. This is the first comprehensive study examining variations in GCPII in relation to breast and prostate cancer risk. Changes in the plasma folate levels and changes in PSMA expression are associated with breast and prostate cancer risk respectively.     

Duffy antigen receptor genetic variant and the association with Interleukin 8 levels

The aim of this study is to identify loci associated with circulating levels of Interleukin 8 (IL8). We investigated the associations of 121,445 single nucleotide polymorphisms (SNPs) from the Illumina 200 K CardioMetabochip with IL8 levels in 1077 controls from the Stockholm Heart Epidemiology Program (SHEEP) study, using linear regression under an additive model of inheritance.Five SNPs (rs12075A/G, rs13179413C/T, rs6907989T/A, rs9352745A/C, rs1779553T/C) reached the pre-defined threshold of genome-wide significance (p < 1.0 × 10⁻⁵) and were tested for in silico replication in three independent populations, derived from the PIVUS, MDC-CC and SCARF studies. IL8 was measured in serum (SHEEP, PIVUS) and plasma (MDC-CC, SCARF). The strongest association was found with the SNP rs12075 A/G, Asp42Gly (p = 1.6 × 10⁻⁶), mapping to the Duffy antigen receptor for chemokines (DARC) gene on chromosome 1. The minor allele G was associated with 15.6% and 10.4% reduction in serum IL8 per copy of the allele in SHEEP and PIVUS studies respectively. No association was observed between rs12075 and plasma IL8.Conclusion: rs12075 was associated with serum levels but not with plasma levels of IL8. It is likely that serum IL8 represents the combination of levels of circulating plasma IL8 and additional chemokine liberated from the erythrocyte DARC reservoir due to clotting. These findings highlight the importance of understanding IL8 as a biomarker in cardiometabolic diseases.     

Variant allele of <Emphasis Type="Italic">CHEK2</Emphasis> is associated with a decreased risk of esophageal cancer lymph node metastasis in a Chinese population

Growing evidence suggests that the checkpoint kinase 2 (CHEK2) signaling pathway occupies a central position in the signaling networks of DNA-damage signaling. Many functional and molecular epidemiological studies have evaluated the association between genetic variants of CHEK2 and various cancers. To evaluate the relationship between CHEK2 functional genetic variants and esophageal cancer risk and the risk of lymph node metastasis among a Chinese population. We genotyped CHEK2 rs738722, rs2236141 and rs2236142 single nucleotide polymorphisms (SNPs) using the matrix assisted laser desorption/ionization time-of-flight mass spectrometry assay in a case–controlled study, including 380 esophageal cancer cases and 380 healthy controls in a Chinese population. We found that none of the three polymorphisms achieved significant difference in their distributions between esophageal cancer cases and controls. Multiple logistic regression analyses revealed that esophageal cancer risk was not associated significantly with the variant genotypes of the three CHEK2 polymorphisms as compared with their wild-type genotypes. However, we found that functional variant rs738722 and rs2236142 in CHEK2 might contribute to susceptibility to lymph node metastasis. Our data did not support a significant association between CHEK2 SNPs and the risk of esophageal cancer. Functional variant CHEK2 rs738722 and rs2236142 might contribute to lymph node metastasis susceptibility. The CT allele of SNP rs738722 and the GC allele of SNP rs2236142 might be a protective factor of the risk for lymph node metastasis of esophageal cancer.     

Innate Immunity Pathways and Breast Cancer Risk in African American and European-American Women in the Women’s Circle of Health Study (WCHS)

African American (AA) women are more likely than European American (EA) women to be diagnosed with early, aggressive breast cancer. Possible differences in innate immune pathways (e.g., inflammatory responses) have received little attention as potential mechanisms underlying this disparity. We evaluated distributions of selected genetic variants in innate immune pathways in AA and EA women, and examined their associations with breast cancer risk within the Women’s Circle of Health Study (WCHS). In stage I of the study (864 AA and 650 EA women) we found that genotype frequencies for 35 of 42 tested SNPs (18 candidate genes) differed between AAs and EAs (corroborated by ancestry informative markers). Among premenopausal AA women, comparing variant allele carriers to non-carriers, reduced breast cancer risk was associated with CXCL5-rs425535 (OR=0.61, P=0.02), while among EA women, there were associations with TNFA-rs1799724 (OR =2.31, P =0.002) and CRP-rs1205 (OR=0.54, P=0.01). For postmenopausal women, IL1B-rs1143627 (OR=1.80, P=0.02) and IL1B-rs16944 (OR=1.85, P =0.02) were associated with risk among EA women, with significant associations for TNFA-rs1799724 limited to estrogen receptor (ER) positive cancers (OR=2.0, P =0.001). However, none of the SNPs retained significance after Bonferroni adjustment for multiple testing at the level of P0.0012 (0.05/42) except for TNFA-rs1799724 in ER positive cancers. In a stage II validation (1,365 AA and 1,307 EA women), we extended evaluations for four SNPs (CCL2-rs4586, CRP-rs1205, CXCL5-rs425535, and IL1RN-rs4251961), which yielded similar results. In summary, distributions of variants in genes involved in innate immune pathways were found to differ between AA and EA populations, and showed differential associations with breast cancer according to menopausal or ER status. These results suggest that immune adaptations suited to ancestral environments may differentially influence breast cancer risk among EA and AA women.     

MicroRNA Related Polymorphisms and Breast Cancer Risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88–0.96), rs1052532 (OR 0.97; 95% CI: 0.95–0.99), rs10719 (OR 0.97; 95% CI: 0.94–0.99), rs4687554 (OR 0.97; 95% CI: 0.95–0.99, and rs3134615 (OR 1.03; 95% CI: 1.01–1.05) located in the 3′ UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.     

Frondoside A Suppressive Effects on Lung Cancer Survival,Tumor Growth,Angiogenesis, Invasion,and Metastasis

A major challenge for oncologists and pharmacologists is to develop less toxic drugs that will improve the survival of lung cancer patients. Frondoside A is a triterpenoid glycoside isolated from the sea cucumber, Cucumaria frondosa and was shown to be a highly safe compound. We investigated the impact of Frondoside A on survival, migration and invasion in vitro, and on tumor growth, metastasis and angiogenesis in vivo alone and in combination with cisplatin. Frondoside A caused concentration-dependent reduction in viability of LNM35, A549, NCI-H460-Luc2, MDA-MB-435, MCF-7, and HepG2 over 24 hours through a caspase 3/7-dependent cell death pathway. The IC50 concentrations (producing half-maximal inhibition) at 24 h were between 1.7 and 2.5 µM of Frondoside A. In addition, Frondoside A induced a time- and concentration-dependent inhibition of cell migration, invasion and angiogenesis in vitro. Frondoside A (0.01 and 1 mg/kg/day i.p. for 25 days) significantly decreased the growth, the angiogenesis and lymph node metastasis of LNM35 tumor xenografts in athymic mice, without obvious toxic side-effects. Frondoside A (0.1–0.5 µM) also significantly prevented basal and bFGF induced angiogenesis in the CAM angiogenesis assay. Moreover, Frondoside A enhanced the inhibition of lung tumor growth induced by the chemotherapeutic agent cisplatin. These findings identify Frondoside A as a promising novel therapeutic agent for lung cancer.     

Genetic Variants in Hormone-Related Genes and Risk of Breast Cancer

Sex hormones play a key role in the development of breast cancer. Certain polymorphic variants (SNPs and repeat polymorphisms) in hormone-related genes are associated with sex hormone levels. However, the relationship observed between these genetic variants and breast cancer risk has been inconsistent. We conducted a case-control study nested within two prospective cohorts to assess the relationship between specific genetic variants in hormone-related genes and breast cancer risk. In total, 1164 cases and 2111 individually-matched controls were included in the study. We did not observe an association between potential functional genetic polymorphisms in the estrogen pathway, SHBG rs6259, ESR1 rs2234693, CYP19 rs10046 and rs4775936, and UGT1A1 rs8175347, or the progesterone pathway, PGR rs1042838, with the risk of breast cancer. Our results suggest that these genetic variants do not have a strong effect on breast cancer risk.     

RhoT1 and Smad4 Are Correlated with Lymph Node Metastasis and Overall Survival in Pancreatic Cancer

Cancer cell invasion and metastasis are the most important adverse prognostic factors for pancreatic cancer. Identification of biomarkers associated with outcome of pancreatic cancer may provide new approaches and targets for anticancer therapy. The aim of this study is to examine the relationship between the expression of RhoT1, Smad4 and p16 and metastasis and survival in patients with pancreatic cancer. The analysis showed that the high cytoplasmic expression levels of RhoT1, Smad4 and p16 in pancreatic cancer tissues had significantly negative correlation with lymph node metastasis (LNM) (P = 0.017, P = 0.032, P = 0.042, respectively). However, no significant association was observed between perineural invasion (PNI) and the expression of above three proteins (all P>0.05). Additionally, the survival analysis showed that the low expression levels of RhoT1 and Smad4 were significantly associated with worse survival (P = 0.034, P = 0.047, respectively). In conclusion, these results indicated that the low-expression levels of RhoT1 and Smad4 were significantly associated with LNM and shorter survival. RhoT1 may be considered as a potential novel marker for predicting the outcome in patients with pancreatic cancer.     

Tumor Invasiveness,Not Lymphangiogenesis,Is Correlated with Lymph Node Metastasis and Unfavorable Prognosis in Young Breast Cancer Patients (≤35 Years)

The morbidity rate of breast cancer is on the rise, and the age of onset appears to be trending toward a young age. Breast cancer in young women (BCYW) has a number of distinctive features that differ from breast cancer in middle-aged or elderly women (BCMEW). Lymphatic metastasis plays an important role in the spread of BCYW; however, the mechanisms of lymph node metastasis (LNM) in BCYW are not clear. This study aimed to investigate the mechanism of lymphatic metastasis in BCYW and to evaluate the relationships between lymphangiogenesis, the expression of matrix metalloproteinase 9 (MMP-9) and vascular endothelial growth factor C (VEGF-C) expression, clinicopathological characteristics, and prognosis. Using immunohistochemistry, MMP-9, VEGF-C and the level of lymphatic microvessel density (LMVD) were analyzed in 106 cases of breast invasive ductal carcinoma and 20 cases of breast proliferative lesions. Compared with BCMEW, BCYW had higher MMP-9 expression, higher LNM, and more adverse prognoses. In BCYW, high MMP-9 expression was positively correlated with LNM and impaired survival time. However, in BCMEW, MMP-9 expression was not correlated with LNM or survival time. In addition, high VEGF-C expression was positively correlated with a high level of LMVD in both BCYW and BCMEW. Nevertheless, a high level of LMVD was not correlated with LNM or survival time in the two groups. More importantly, univariate and multivariate survival analysis showed that MMP-9 expression and LNM were independent prognostic factors in BCYW. Our present study indicates that lymphangiogenesis induced by VEGF-C is augmented in breast cancer; however, a higher level of lymphangiogenesis has no significant impact on LNM or survival time. We suggest that tumor invasiveness, rather than lymphangiogenesis, plays an important role in LNM among BCYW. Moreover, MMP-9 and LNM were independent prognostic factors for BCYW.