Abstract
Song Zhang1, Dong Zhang1, Shanhong Yi2, Mingfu Gong1, Caibao Lu3, Yuanqing Cai1, Xuefeng Tang4, Liguang Zou1
1Department of Radiology, Xinqiao Hospital, Third Military Medical University, Chongqing 400037, China
2Department of Urology, Xinqiao Hospital, Third Military Medical University, Chongqing 400037, China
3Department of Nephrology, Xinqiao Hospital, Third Military Medical University, Chongqing 400037, China
4Department of Pathology, Xinqiao Hospital, Third Military Medical University, Chongqing 400037, China
Correspondence to:
Liguang Zou, email: [email protected], [email protected]
Keywords: lymphatic microvessel density, lymphovascular invasion, lymph node metastasis, breast cancer, meta-analysis
Received: August 04, 2016 Accepted: November 22, 2016 Published: December 01, 2016
ABSTRACT
Lymph node status is one of the key parameters used for determining the stage of breast cancer progression. The relationship of lymphatic vessel density (LVD), lymphovascular invasion (LVI), and lymph node metastasis (LNM) has not been clearly demonstrated yet. Databases of PubMed, Embase, and Web of Science were searched from inception up to 25 May 2016. Spearman correlation coefficient (r) and 95% confidence interval (CI) were used to determine the relationship within each group. Based on pre-established inclusion criteria, 28 studies involving 2920 breast cancer patients were included in this study. The r values of LVD-LVI, LVD-LNM, and LVI-LNM were 0.45 (95% CI: 0.31 to 0.57), 0.32 (95% CI: 0.23 to 0.40), and 0.24 (95% CI: 0.19 to 0.28), respectively. Compared with intratumoral LVD, peritumoral LVD showed more robust correlation with LVI (r = 0.53, 95% CI: 0.27 to 0.72) and LNM (r = 0.33, 95% CI: 0.18 to 0.46). The patients in LNM positive group presented with higher LVI detection rate of 45.85%, while in LNM negative group with detection rate of 23.85%. The results describe a triangle relationship between LVD, LVI, and LNM in breast cancer. Both LVD and LVI are indicated to be valuable predictors of LNM occurrence. Compared with intratumoral lymphatic vessels, peritumoral lymphatics might be the main disseminate route for breast tumor cells.
INTRODUCTION
Cancer metastasis is the leading cause of mortality in patients diagnosed with breast cancer and other malignant tumors [1]. Lymph node status is commonly used to identify a patient’s prognosis, tumor stage, and treatment modality [2]. Patients without lymph node metastasis (LNM) have a favorable prognosis, while with more than six positive axillary lymph nodes have a higher risk of distant metastasis [3]. The progress of lymphatic metastasis is thought to involve the proliferation of lymphatic vessels (lymphangiogenesis), lymphovascular invasion, and LNM step by step [1]. However, the mechanism leading to tumor cells spread via lymphatic vessels (lymphovascular invasion, LVI) to the regional and distant lymph nodes has not been clearly demonstrated [4].
Lymphatic vessel was thought to play a passive role in tumor metastasis, due to the absence of reliable molecular markers to distinguish lymphatic vessels and the lack of identified growth factors for the lymphatic system. During the last two decades, substantial progress within the field has rapidly lead to the recognition of the lymphatic system as an active player involved in lots of malignant tumors [5]. Lymphatic vessels not only provide an entrance for tumor cells to penetrate in [6], but also make several key contributions to tumor metastasis, such as the provision of a niche for cancer stem cells and the modulation of antitumor immune responses [4]. It is known that tumor angiogenesis and its indicator blood vessel denstiy are closely associated with the clinicopathological outcomes in breast cancer [7]. However, the clinical role of lymphangiogenesis and its indicator lymph vessel density (LVD) is needed to be further investigated [8, 9].
Despite lymphatic metastasis is thought to facilitate tumor cell dissemination in breast cancer [10], the possible relationship of LVD, LVI, and LNM remains ambiguous. Questions can be listed as follows: (1) do lymphangiogenesis facilitate the occurrence of LVI and LNM in breast cancer? (2) regarding LVI detection rate, is there any difference between LNM negative group and LNM positive group? (3) whether lymphangiogenesis is previous to LVI, or it is promoted by the chemokines secreted by tumor cells penetrated into lymphatic vessels? The answers to these questions are important to understand whether tumor-induced lymphangiogenesis is a potential target for the inhibition of distant metastasis, as well as whether high LVD and LVI presence are valuable factors to predict LNM occurrence in breast cancer. With the accumulating evidence, we conducted a systematic review and meta-analysis to estimate the relationship of LVD, LVI, and LNM in breast cancer to provide insights of the above issues.
RESULTS
Study selection process
The literature search result is shown in the flowchart of Figure 1. We initially searched 1134 potential studies from the databases. After removing the duplicated and irrelevant publications, 79 full-text publications were left over to assess the eligibility. Fifty-one papers were excluded due to not fulfill the inclusion criteria, inappropriate publication types, or insufficient data. Finally, 28 articles were included in the analysis [11–38].
Figure 1: Flowchart of study selection.
Characteristics of the included studies
The results of all included 28 studies were exhibited in Tables 1–3. The sample size of each study ranged from 29 to 374 patients, and the publication year of them ranged from 2000 to 2016. A total of 2920 breast cancer patients were adopted in this study. All patients were performed surgical treatments, and IHC staining with D2–40, podoplanin, LYVE-1, and VEGFR-3 antibodies. LVD was determined by counting the number of lymphatic vessels per area under a microscope. LVI was defined as the presence of tumor cells in lymphatic vessels. Lymph nodes, either sentinel lymph nodes or non-sentinel lymph nodes, were taken into account to determine the occurrence of LNM.
Table 1: Main characteristics and results of the studies evaluating LVD and LVI
a mean (range); b median (range); c mean ± SD; BC, breast cancer; H, high LVD; L, low LVD; NG, not given.
Table 2: Main characteristics and results of the studies evaluating LVD and LNM
a mean (range); b median (range); c mean ± SD; BC, breast cancer; H, high LVD; L, low LVD; NG, not given.
Table 3: Main characteristics and results of the studies evaluating LVI and LNM
Author, Year, Country | No. of patients | Age | Tumor type | Antibody dilution | LVI of LNM-(No.) | LVI of LNM+(No.) | r (95% CI) |
|---|---|---|---|---|---|---|---|
Kanngurn, 2013, Thailand [23] | 122 | 52 (29–86)b | invasive primary BC | D2-40 (1:200) | N (68), P (11) | N (21), P (18) | 0.35 (0.17, 0.51) |
Widodo, 2013, Indonesia [36] | 48 | 53.0 (34–75)a | breast cancer | D2-40 (1:75) | N (6), P (12) | N (7), P (23) | 0.11 (−0.17, 0.38) |
Kandemir, 2012, Turkey [22] | 69 | 54.8 (39–85)a | ductal invasive BC | D2-40 (1:50) | 0.06 ± 0.05 (15) | 0.19 ± 0.21 (28) | 0.32 (0.03, 0.57) |
Lee, 2010, Korea [25] | 46 | 47.9 ± 2.5c | microinvasive ductal BC | D2-40 (1:130) | N (33), P (4) | N (6), P (3) | 0.25 (−0.06, 0.51) |
Britto, 2009, Brazil [14] | 92 | 55 (32–77)b | BC | D2-40 (1:50) | N (44), P (17) | N (21), P (10) | 0.05 (−0.16, 0.25) |
Braun, 2008, Germany [13] | 254 | 57 (28–85)b | primary invasive BC | D2-40 (1:50) | N (114), P (20) | N (49), P (44) | 0.35 (0.22, 0.47) |
El-Gohary, 2008, USA [18] | 48 | 64 (27–89)a | invasive BC | D2-40 (1:50) | NG (24) | NG (24) | 0.52 (0.27, 0.70) |
Marinho, 2008, Brazil [26] | 123 | 52 (27–88)b | invasive BC | D2-40 (1:100) | N (32), P (9) | N (56), P (26) | 0.10 (−0.08, 0.27) |
Ito, 2007, Japan [21] | 69 | 52.1 (27–80)a | invasive BC | D2-40 (1:200) | N (37), P (7) | N (16), P (9) | 0.23 (−0.01, 0.45) |
Tezuka, 2007, Japan [32] | 132 | 55.9 (31–84)b | invasive BC | D2-40 (NG) | N (42), P (21) | N (35), P (34) | 0.16 (−0.01, 0.32) |
van den Eynden, 2006, Belgium [33] | 95 | 60.5 (33.5–86.1)a | invasive BC | D2-40 (1:100) | N (19), P (33) | N (9), P (33) | 0.16 (−0.04, 0.36) |
Kato, 2005, UK [24] | 67 | 49 (30–86)b | primary BC | LYVE-1 (1:600) | N (31), P (12) | N (8), P (12) | 0.31 (0.06, 0.52) |
Schoppmann, 2004, Austria [31] | 374 | 57.6 (median age) | invasive BC | podoplanin (1:200) | N (171), P (41) | N (98), P (64) | 0.22 (0.12, 0.32) |
a mean (range); b median (range); c mean ± SD; BC, breast cancer; N, LVI negative; P, LVI positive; NG, not given.
Data analysis
A total of 1221 patients involved in 12 studies (Table 1), which provided sufficient data to determine the relationship between LVD and LVI. According to the detection area of lymphatic vessels, the included 12 studies were divided into three subgroups as follows: total LVD (without distinguishing peritumoral and intratumoral) (n = 7) [11, 17, 24, 25, 27, 29, 31], peritumoral LVD (n = 5) [16, 18, 22, 36, 38], and intratumoral LVD (n = 3) [18, 22, 38]. The main outcomes are summarized in Figure 2. The pooled Fisher’s Z values of total LVD, peritumoral LVD, intratumoral LVD, overall LVD were 0.38 (95% CI: 0.19 to 0.57, I2 = 83.2%, P = 0.000), 0.59 (95% CI: 0.28 to 0.90, I2 = 79.2%, P = 0.001), 0.84 (95% CI: −0.06 to 1.74, I2 = 95%, P = 0.000), 0.48 (95% CI: 0.32 to 0.65, I2 = 85.6%, P = 0.000), respectively. The pooled Fisher’s Z values were converted back to r values. Both total LVD (r = 0.36, 95% CI: 0.19 to 0.52) and overall LVD (r = 0.45, 95% CI: 0.31 to 0.57) were moderately correlated with LVI, while peritumoral LVD (r = 0.53, 95% CI: 0.27 to 0.72) and intratumoral LVD (r = 0.69, 95% CI: −0.06 to 0.94) showed strong correlation with LVI (Figure 2, Table 4). However, only three sets of intratumoral LVD data were acquired. Significant evidence of heterogeneity was noted among these studies (P = 0.05, I2 = 95%).
Table 4: Pooled Z values, pooled r values, Begg’s and Egger’s results of each correlation group
Relation group | Pooled Z value(95% CI) | Pooled r value(95% CI) | Begg’s test | Egger’s test |
|---|---|---|---|---|
total LVD-LVI (n = 7) | 0.38 (0.19, 0.57) | 0.36 (0.19, 0.52) | 0.356 | 0.678 |
peritumoral LVD-LVI (n= 5) | 0.59 (0.28, 0.90) | 0.53 (0.27, 0.72) | 0.068 | 0.005 |
intratumoral LVD-LVI (n= 3) | 0.84 (−0.06, 1.74) | 0.69 (−0.06, 0.94) | 0.296 | 0.225 |
overall LVD-LVI (n= 15) | 0.48 (0.32, 0.65) | 0.45 (0.31, 0.57) | 0.700 | 0.000 |
total LVD-LNM (n = 15) | 0.33 (0.19, 0.48) | 0.32 (0.19, 0.45) | 0.023 | 0.038 |
peritumoral LVD-LNM (n = 8) | 0.34 (0.18, 0.50) | 0.33 (0.18, 0.46) | 0.711 | 0.321 |
intratumoral LVD-LNM (n = 5) | 0.30 (0.04, 0.55) | 0.29 (0.04, 0.50) | 0.806 | 0.283 |
overall LVD-LNM (n = 28) | 0.33 (0.23, 0.42) | 0.32 (0.23, 0.40) | 0.047 | 0.005 |
LVI-LNM (n= 13) | 0.24 (0.19, 0.29) | 0.24 (0.19, 0.28) | 0.428 | 0.736 |
Figure 2: Forest plot of the Fisher’s Z values for the correlation between LVD and LVI in breast cancer.
Twenty-two studies involving 2125 patients were included in the assessment of the correlation between LVD and LNM (Table 2). According to the detection area of lymphatic vessels, the studies were also divided into three subgroups of total LVD (n = 15) [11, 12, 14, 15, 17, 19, 20, 24, 25, 27–31, 37], peritumoral LVD (n = 8) [16, 18, 22, 27, 34–36, 38], and intratumoral LVD (n = 5) [18, 22, 27, 34, 38]. The pooled Fisher’s Z values were show in Figure 3, corresponding to the r values of total LVD (r = 0.32, 95% CI: 0.19 to 0.45), peritumoral LVD (r = 0.33, 95% CI: 0.18 to 0.46), intratumoral LVD (r = 0.29, 95% CI: 0.04 to 0.50), and overall LVD (r = 0.32, 95% CI: 0.23 to 0.40), respectively (Figure 3, Table 4).
Figure 3: Forest plot of the Fisher’s Z values for the correlation between LVD and LNM in breast cancer.
Among the included 13 studies, which provided sufficient data to evaluate the relationship between LVI and LNM, eleven of them [13, 14, 21, 23–26, 31–33, 36] described the detailed number of LVI occurrence within LNM negative group and LNM positive group. In LNM negative group (n = 784), 187 (23.85%) patients presented with LVI; while in LNM positive group (n = 602), 276 (45.85%) patients presented with LVI. To evaluate the correlation of LVI and LNM, all aquired data were used to obtain the r values, and then were transformed to the Fisher’s Z values. The pooled Fisher’s Z value was 0.24 (95% CI: 0.19 to 0.29, I2 = 34.8%, P = 0.104, Figure 4), and its corresponding r value was 0.24 (95% CI: 0.19 to 0.28).
Figure 4: Forest plot of the Fisher’s Z values for the correlation between LVI and LNM in breast cancer.
Sensitivity analysis and publication bias
In order to assess the stability of the results, sensitivity analyses were independently performed in the groups of LVD-LVI, LVD-LNM, and LVI-LNM. By removing individual studies in turn, sensitivity analyses demonstrated no disproportionate pooled estimates, indicating a statistically robust result of the analysis (Figure S1–S3 in Supplementary file 1). Begg’s tests and the funnel plots of the Z value against the standard error of Z value showed substantial asymmetry (Figure S4–S6 in Supplementary file 1). The results of Begg’s and Egger’s tests are displayed in Table 4.
DISCUSSION
The current meta-analysis included 28 studies with an overall population of 2920 breast cancer patients. Our study reveals the triangle relationship of LVD, LVI, and LNM in breast cancer (Figure 5). Peritumoral LVD shows the most robust correlation with LVI and LNM, while intratumoral LVD and total LVD presents with a relatively weak correlation (Figure 5). Patients in LNM positive group shows higher LVI detection rate than that of LNM negative group. The results demonstrate that both LVD and LVI are valuable predictors of the LNM occurrence in breast cancer. However, the evidence of heterogeneity was observed across the studies, which needed to be further investigated.
Figure 5: Correlation models of the meta-analysis study. (A) Relationships of LVD-LVI and LVD-LNM in the subgroups of general LVD, intratumoral LVD, and peritumoral LVD; (B) Triangle relationship model of overall LVD, LVI and LNM in breast cancer.
Heterogeneity is a potential problem to interpret the meta-analysis results. Significant heterogeneities between-studies were presented. The detailed LVD values and LVI detection rate differentiate notably among these studies. The variation might be caused by patient sources, staining techniques, antibody categories and antibody dilutions. In addition, different counting methods of lymphatic vessel density, by using different hotspots (three [19], four [34], and five [11]), magnification field (100× [27], 200× [11], 400× [19]), and measuring unit (vessels/mm2 [11], vessels/field [31]), are also accounted for the variation of results. Furthermore, the cutoff value to divide LVD and LVI as low and high is a crucial factor that cannot be ignored. Because the asset values of LVD and LVI is not normal distribution, most of the included studies chose the median value as the cutoff value, a few of them took the mean or actual value as the cutoff value. Therefore, studies with more standardized and stricter design are required for the assessment of lymphatic vessel density.
Over decades, lymphatic vessels have been described as a passive participant in metastasis and regarded as only a transportation channel of tumor cells. It is still uncertain whether a high LVD is a necessary condition for metastasis [39]. Zhang et al. demonstrated that tumor invasion, but not lymphangiogenesis, was correlated with LNM and unfavorable prognosis in young breast cancer patients [37]. Other studies found that the LVD in LVI/LNM negative group even higher than that of LVI/LNM positive group in primary breast cancer [17, 24, 38]. Nonetheless, most of the included studies showed a positive correlation between a high LVD and the presence of LVI/LNM. The meta-analysis results showed positive correlations between LVD and LVI (r = 0.45, 95% CI: 0.31 to 0.57), LVD and LNM (r = 0.32, 95% CI: 0.23 to 0.40). It indicates that lymphangiogenesis may facilitate the interaction between tumor cells and lymphatic vessels, thereby increasing the probability of tumor cells invasion and distant metastasis. In addition, lymphatic vessel may provide a safe route for cancer cell dissemination, due to the discontinuous structure of the lymphatic basement membrane, a miniature shear stress and a high concentration of hyaluronic acid [39].
The presence of intratumoral lymphatic vessels is a hotly debated issue in malignant solid tumors, particularly in breast cancer [12, 40, 41]. Previous studies reported that solid tumors did not have intratumoral lymphatic vessels [42], because of the increasing interstitial pressure induced by the proliferating tumor cells [43]. With the application of specific lymphatic vessel markers, many studies have demonstrated the presence of intratumoral lymphatic vessels [27, 44]. Moreover, intratumoral lymphatic vessels are indicated to be functional, because tumor cells can be found within the vessels [45]. Another interesting issue is to what extent intratumoral lymphatic vessels and peritumoral lymphatic vessels participate in tumor cells dissemination. Our study shows that high peritumoral LVD strongly correlates with LVI (r = 0.53, 95% CI: 0.27 to 0.72), while high intratumoral LVD moderately correlates with LVI (r = 0.33, 95% CI: 0.18 to 0.46). The result suggests that peritumoral lymphatic vessels play a more important role on metastatic dissemination in breast cancer.
The detection of LVI was mainly assessed by H&E staining method, due to the deficiency of specific markers of lymphatic vessels [46, 47]. One major challenge of this method is to distinguish LVI from the retraction artifacts caused by tissue handling and fixation. A previous meta-analysis study has summarized the data on the presence of LVI in breast cancer [48]. The results show that the detection rate of LVI widely ranges from 10% to 49% of H&E staining method, while the range is narrower (from 21% to 42%) by using IHC staining method [48]. It indicates that IHC staining method should be more reliable to identify the presence of LVI. Therefore, 13 studies, by using IHC instead of H&E, were included to investigate the relationship between LVI and LNM. Eleven of them reported the detailed LVI presence in LNM negative group (overall detection rate of 23.85%) and positive group (overall detection rate of 45.85%). The presence of LVI shows a weak correlation with the occurrence of LNM in breast cancer (r = 0.24, 95% CI: 0.19 to 0.28).
The current meta-analysis study has some strengths. It is the first study to systematically discuss the triangle relationship of LVD, LVI, and LNM in breast cancer. The included 28 studies and 2290 participants significantly enhanced the statistical power and provided more reliable results. However, some limitations should not be ignored. First, all included studies were observational studies, the sample sizes were relatively small, and several studies were excluded due to lack of sufficient data to determine the correlation coefficients. Thus, selection bias and recall bias are inevitable. Second, most of the included studies investigated lymphatic vessels without distinguishing intratumoral and peritumoral lymphatic vessels, which would confound the final results. Finally, as described above, significant between-studies heterogeneities were presented. Therefore, the standardization of LVD and LVI counting method needs to be established for the future study.
MATERIALS AND METHODS
Literature search
A systematic search of the PubMed, Embase and Web of Science databases was performed to identify all relevant articles published up to 25 May 2016. The following Medical Subject Heading (MeSH) terms or keywords were used: “breast cancer OR breast carcinoma OR breast neoplasms” AND “lymphatic vessel density OR lymphatic microvessel density OR LVD OR LMVD OR lymphangiogenesis OR lymphovascular invasion OR lymphatic vessel invasion OR lymphatic invasion OR LVI OR lymph node metastasis OR LNM”. All abstracts that indicated the correlation between LVD, LVI, and LNM in breast cancer, either prospective or retrospective, were chosen for further consideration.
Inclusion criteria
All studies were required to meet the following criteria: (1) dealt with the patients with primary breast cancer only; (2) published as a full-text research paper, rather than reviews, case reports, meeting abstracts, or animal researches; (3) clearly described the methods and procedures of tissue handling and pathological staining; (4) the determination of LVI presence was assessed by immunohistochemical (IHC) staining instead of hematoxylin and eosin (H&E) staining. When two or more articles reported duplicating data, only the study with the most recent data or the largest dataset was included. Two independent authors followed the inclusion criteria to review the publications. In case of dispute, a third author assessed the study to obtain a consensus.
Data extraction and quality assessment
Data retrieved from the articles included the first author’s name, publication year, country, number of patients, age, type of breast cancer involved, antibody and its dilution, detail data or patient number of LVD/LVI used to determine the correlation between each group. Items intended for extraction were discussed by two authors.
The quality assessment of including studies was based on the criteria of the Newcastle-Ottawa Quality Assessment scale (NOS) (Table S1 in Supplementary file 2) [49]. The scale used a star system to evaluate the study quality, including the aspects of selection, comparability, and exposure. The greater number of stars represented higher-quality studies, which performed better and strictly controlled for potential confounders.
Statistical analysis
The results showed with detailed data of means and standard deviations, or presented in two by two frequency tables, were used to obtain the Spearman correlation coefficients (r) [50]. Because some variables in the original studies were log-transformed before analysis, Spearman’s correlation coefficient, instead of Pearson’s correlation coefficient, was applied in this study [51]. Before the combination of results, a Fisher’s Z transformation was used to convert Spearman’s correlation coefficients into an approximately normal distribution (Figure S7 in Supplementary file 1). The individual Fisher’s Z values with their corresponding 95% confidence interval (95% CI) were pooled to obtain an overall estimate using STATA 12.0 software. The pooled Z value was finally transformed back to r value by an inverse Fisher’s Z transformation (Figure S7 in Supplementary file 1). Correlation coefficients are not precise but were generally classified as weak, moderate, strong, and excellent. In our study, we assumed < 0.25 to be weak or no correlation, 0.25–0.50 to be moderate, 0.50–0.75 to be strong, > 0.75 to be excellent or perfect [52].
Homogeneity test was performed with Q statistic and the I2 statistic. In this study, P < 0.05 or I2 value > 50% were considered to be statistically significant. A random-effects model or, in the absence of heterogeneity, a fixed-effects model was utilized to combine the Z values. If heterogeneity was noted, a sensitivity analysis was conducted to investigate the influence of each study on the overall estimate by omitting each study in turn. Publication bias was detected by Begg’s test and Egger’s test.
CONCLUSIONS
The study describes a triangle relationship of LVD, LVI, and LNM in primary breast cancer. Both LVD and LVI are valuable predict factors of LNM, while LVD moderately correlates with LNM and LVI weakly correlates with LNM. Compared with intratumoral LVD, peritumoral LVD shows a more robust correlation with LVI and LNM, which indicates that peritumoral lymphatic vessel is the main disseminate route for breast tumor cells. It suggests that the patients detected with high LVD or LVI presence, especially with peritumoral high LVD, should take more active treatment to prevent the aggravation and metastasis of primary breast cancer. However, further studies with larger sample sizes should be performed to validate our results.
CONFLICTS OF INTEREST
The authors have no conflicts of interest to declare.
GRANT SUPPORT
This study was supported by the National Natural Science Foundation of China (NSFC) (Nos 81401466, 81501521).
REFERENCES
1. Ran S, Volk L, Hall K, Flister MJ. Lymphangiogenesis and lymphatic metastasis in breast cancer. Pathophysiology. 2010; 17:229–251.
2. Andree C, Schmidt VJ, Munder BI, Seidenstucker K, Behrendt P, Witzel C, Horch RE, Andrews BT, Richrath P. Detecting of breast cancer metastasis by means of regional lymph node sampling during autologous breast reconstruction—a screening of 519 consecutive patients. Med Sci Monit. 2012; 18:Cr605–610.
3. Tseng HS, Chen LS, Kuo SJ, Chen ST, Wang YF, Chen DR. Tumor characteristics of breast cancer in predicting axillary lymph node metastasis. Med Sci Monit. 2014; 20:1155–1161.
4. Karaman S, Detmar M. Mechanisms of lymphatic metastasis. J Clin Invest. 2014; 124:922–928.
5. Tammela T, Alitalo K. Lymphangiogenesis: Molecular mechanisms and future promise. Cell. 2010; 140:460–476.
6. Bando H, Weich HA, Horiguchi S, Funata N, Ogawa T, Toi M. The association between vascular endothelial growth factor-C, its corresponding receptor, VEGFR-3, and prognosis in primary breast cancer: a study with 193 cases. Oncol Rep. 2006; 15:653–659.
7. Kato T, Kameoka S, Kimura T, Nishikawa T, Kobayashi M. The combination of angiogenesis and blood vessel invasion as a prognostic indicator in primary breast cancer. Br J Cancer. 2003; 88:1900–1908.
8. Agarwal B, Saxena R, Morimiya A, Mehrotra S, Badve S. Lymphangiogenesis does not occur in breast cancer. Am J Surg Pathol. 2005; 29:1449–1455.
9. Shayan R, Achen MG, Stacker SA. Lymphatic vessels in cancer metastasis: bridging the gaps. Carcinogenesis. 2006; 27:1729–1738.
10. Cote RJ, Peterson HF, Chaiwun B, Gelber RD, Goldhirsch A, Castiglione-Gertsch M, Gusterson B, Neville AM. Role of immunohistochemical detection of lymph-node metastases in management of breast cancer. Lancet. 1999; 354:896–900.
11. Abe N, Ohtake T, Saito K, Kumamoto K, Sugino T, Takenoshita S. Clinicopathological significance of lymphangiogenesis detected byimmunohistochemistry using D2-40 monoclonalantibody in breast cancer. Fukushima J Med Sci. 2016; 62:1–7.
12. Bono P, Wasenius VM, Heikkila P, Lundin J, Jackson DG, Joensuu H. High LYVE-1-positive lymphatic vessel numbers are associated with poor outcome in breast cancer. Clin Cancer Res. 2004; 10:7144–7149.
13. Braun M, Flucke U, Debald M, Walgenbach-Bruenagel G, Walgenbach KJ, Holler T, Polcher M, Wolfgarten M, Sauerwald A, Keyver-Paik M, Kuhr M, Buttner R, Kuhn W. Detection of lymphovascular invasion in early breast cancer by D2-40 (podoplanin): a clinically useful predictor for axillary lymph node metastases. Breast Cancer Res Treat. 2008; 112:503–511.
14. Britto AV, Schenka AA, Moraes-Schenka NG, Alvarenga M, Shinzato JY, Vassallo J, Ward LS. Immunostaining with D2-40 improves evaluation of lymphovascular invasion, but may not predict sentinel lymph node status in early breast cancer. BMC Cancer. 2009; 9:109.
15. Choi WWL, Lewis MM, Lawson D, Yin-Goen QQ, Birdsong GG, Cotsonis GA, Cohen C, Young AN. Angiogenic and lymphangiogenic microvessel density in breast carcinoma: correlation with clinicopathologic parameters and VEGF-family gene expression. Mod Pathol. 2005; 18:143–152.
16. Ding M, Fu X, Tan H, Wang R, Chen Z, Ding S. The effect of vascular endothelial growth factor C expression in tumor-associated macrophages on lymphangiogenesis and lymphatic metastasis in breast cancer. Mol Med Rep. 2012; 6:1023–1029.
17. El-Gendi S, Abdel-Hadi M. Lymphatic vessel density as prognostic factor in breast carcinoma: relation to clinicopathologic parameters. J Egypt Natl Canc Inst. 2009; 21:139–149.
18. El-Gohary YM, Metwally G, Saad RS, Robinson MJ, Mesko T, Poppiti RJ. Prognostic significance of intratumoral and peritumoral lymphatic density and blood vessel density in invasive breast carcinomas. Am J Clin Pathol. 2008; 129:578–586.
19. Gu Y, Qi X, Guo S. Lymphangiogenesis induced by VEGF-C and VEGF-D promotes metastasis and a poor outcome in breast carcinoma: a retrospective study of 61 cases. Clin Exp Metastas. 2008; 25:717–725.
20. Guo X, Chen L, Lang R, Fan Y, Zhang X, Fu L. Invasive micropapillary carcinoma of the breast: association of pathologic features with lymph node metastasis. Am J Clin Pathol. 2006; 126:740–746.
21. Ito M, Moriya T, Ishida T, Usami S, Kasajima A, Sasano H, Ohuchi N. Significance of pathological evaluation for lymphatic vessel invasion in invasive breast cancer. Breast Cancer-Tokyo. 2007; 14:381–387.
22. Kandemir NO, Barut F, Bektas S, Ozdamar SO. Can lymphatic vascular density be used in determining metastatic spreading potential of tumor in invasive ductal carcinomas? Pathol Oncol Res. 2012; 18:253–262.
23. Kanngurn S, Thongsuksai P, Chewatanakornkul S. Chalkley microvessel but not lymphatic vessel density correlates with axillary lymph node metastasis in primary breast cancers. Asian Pac J Cancer Prev. 2013; 14:583–587.
24. Kato T, Prevo R, Steers G, Roberts H, Leek RD, Kimura T, Kameoka S, Nishikawa T, Kobayashi M, Jackson DG, Harris AL, Gatter KC, Pezzella F. A quantitative analysis of lymphatic vessels in human breast cancer, based on LYVE-1 immunoreactivity. Brit J Cancer. 2005; 93:1168–1174.
25. Lee K, Park DJ, Choe G, Kim HH, Kim WH, Lee HS. Increased Intratumoral Lymphatic Vessel Density Correlates with Lymph Node Metastasis in Early Gastric Carcinoma. Ann Surg Oncol. 2010; 17:73–80.
26. Marinho VF, Metze K, Sanches FS, Rocha GF, Gobbi H. Lymph vascular invasion in invasive mammary carcinomas identified by the endothelial lymphatic marker D2-40 is associated with other indicators of poor prognosis. BMC Cancer. 2008; 8:64.
27. Mohammed RA, Ellis IO, Elsheikh S, Paish EC, Martin SG. Lymphatic and angiogenic characteristics in breast cancer: morphometric analysis and prognostic implications. Breast Cancer Res Treat. 2009; 113:261–273.
28. Mylona E, Nomikos A, Alexandrou P, Giannopoulou I, Keramopoulos A, Nakopoulou L. Lymphatic and blood vessel morphometry in invasive breast carcinomas: Relation with proliferation and VEGF-C and -D proteins expression. Histol Histopathol. 2007; 22:825–835.
29. Nakamura Y, Yasuoka H, Tsujimoto M, Imabun S, Nakahara M, Nakao K, Nakamura M, Mori I, Kakudo K. Lymph vessel density correlates with nodal status, VEGF-C expression, and prognosis in breast cancer. Breast Cancer Res Treat. 2005; 91:125–132.
30. Nathanson SD, Zarbo RJ, Wachna DL, Spence CA, Andrzejewski TA, Abrams J. Microvessels that predict axillary lymph node metastases in patients with breast cancer. Arch Surg. 2000; 135:586–593.
31. Schoppmann SF, Bayer G, Aumayr K, Taucher S, Geleff S, Rudas M, Kubista E, Hausmaninger H, Samonigg H, Gnant M, Jakesz R, Horvat R. Prognostic value of lymphangiogenesis and lymphovascular invasion in invasive breast cancer. Ann Surg. 2004; 240:306–312.
32. Tezuka K, Onoda N, Takashima T, Takagaki K, Ishikawa T, Wakasa T, Wakasa K, Hirakawa K. Prognostic significance of lymphovascular invasion diagnosed by lymphatic endothelium immunostaining in breast cancer patients. Oncol Rep. 2007; 17:997–1003.
33. Van den Eynden GG, Van der Auwera I, Van Laere SJ, Colpaert CG, van Dam P, Dirix LY, Vermeulen PB, Van Marck EA. Distinguishing blood and lymph vessel invasion in breast cancer: a prospective immunohistochemical study. Br J Cancer. 2006; 94:1643–1649.
34. van der Schaft DW, Pauwels P, Hulsmans S, Zimmermann M, van de Poll-Franse LV, Griffioen AW. Absence of lymphangiogenesis in ductal breast cancer at the primary tumor site. Cancer letters. 2007; 254:128–136.
35. van Iterson V, Leidenius M, von Smitten K, Bono P, Heikkila P. VEGF-D in association with VEGFR-3 promotes nodal metastasis in human invasive lobular breast cancer. Am J Clin Pathol. 2007; 128:759–766.
36. Widodo I, Ferronika P, Harijadi A, Triningsih FE, Utoro T, Soeripto. Clinicopathological significance of lymphangiogenesis and tumor lymphovascular invasion in indonesian breast cancers. Asian Pac J Cancer Prev. 2013; 14:997–1001.
37. Zhang Z-Q, Han Y-Z, Nian Q, Chen G, Cui S-Q, Wang X-Y. Tumor Invasiveness, Not Lymphangiogenesis, Is Correlated with Lymph Node Metastasis and Unfavorable Prognosis in Young Breast Cancer Patients (<= 35 Years). PLoS One. 2015; 10:e0144376.
38. Zhao YC, Ni XJ, Li Y, Dai M, Yuan ZX, Zhu YY, Luo CY. Peritumoral lymphangiogenesis induced by vascular endothelial growth factor C and D promotes lymph node metastasis in breast cancer patients. World J Surg Oncol. 2012; 10:165.
39. Paduch R. The role of lymphangiogenesis and angiogenesis in tumor metastasis. Cell Oncol. 2016; 10.1007/s13402–016–0281–9.
40. Williams CS, Leek RD, Robson AM, Banerji S, Prevo R, Harris AL, Jackson DG. Absence of lymphangiogenesis and intratumoural lymph vessels in human metastatic breast cancer. J Pathol. 2003; 200:195–206.
41. Vleugel MM, Bos R, van der Groep P, Greijer AE, Shvarts A, Stel HV, van der Wall E, van Diest PJ. Lack of lymphangiogenesis during breast carcinogenesis. J Clin Pathol. 2004; 57:746–751.
42. Padera TP, Kadambi A, di Tomaso E, Carreira CM, Brown EB, Boucher Y, Choi NC, Mathisen D, Wain J, Mark EJ. Lymphatic metastasis in the absence of functional intratumor lymphatics. Science. 2002; 296:1883–1886.
43. Clarijs R, Ruiter DJ, de Waal RM. Lymphangiogenesis in malignant tumours: does it occur? J Pathol. 2001; 193: 143–146.
44. Tezuka K, Onoda N, Takashima T, Ishikawa T, Wakasa T, Wakasa K, Hirakawa K. Clinical significance of intra-tumoral sinusoidal structures showing lympho-endothelial immunoreactivity in breast cancer. Oncol Rep. 2008; 20:25–32.
45. Pathak AP, Bhujwalla ZM, Pepper MS. Visualizing function in the tumor-associated lymphatic system. Lymphat Res Biol. 2004; 2:165–172.
46. Nathanson SD, Kwon D, Kapke A, Alford SH, Chitale D. The role of lymph node metastasis in the systemic dissemination of breast cancer. Ann Surg Oncol. 2009; 16:3396–3405.
47. Jinno H, Sakata M, Asaga S, Wada M, Shimada T, Kitagawa Y, Suzuki T, Nakahara T, Kitamura N, Kubo A, Mukai M, Ikeda T, Kitajima M. Predictors to assess non-sentinel lymph node status in breast cancer patients with sentinel lymph node metastasis. Breast J. 2008; 14:551–555.
48. Gujam FJ, Going JJ, Edwards J, Mohammed ZM, McMillan DC. The role of lymphatic and blood vessel invasion in predicting survival and methods of detection in patients with primary operable breast cancer. Crit Rev Oncol Hematol. 2014; 89:231–241.
49. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol. 2010; 25:603–605.
50. Lin J, Liu J, Davies ML, Chen W. Serum Vitamin D Level and Rheumatoid Arthritis Disease Activity: Review and Meta-Analysis. PloS one. 2016; 11:e0146351.
51. Chen L, Liu M, Bao J, Xia Y, Zhang J, Zhang L, Huang X, Wang J. The correlation between apparent diffusion coefficient and tumor cellularity in patients: a meta-analysis. PLoS One. 2013; 8:e79008.
52. Blackman AJ, Smith MV, Flanigan DC, Matava MJ, Wright RW, Brophy RH. Correlation between magnetic resonance imaging and clinical outcomes after cartilage repair surgery in the knee a systematic review and meta-analysis. Am J Sport Med. 2013; 41:1426–1434.
