Skip to main content

Expression profiles of VEGF-A, VEGF-D and VEGFR1 are higher in distant metastases than in matched primary high grade epithelial ovarian cancer



In many malignancies including ovarian cancer, different angiogenic factors have been related to poor prognosis. However, data on their relations to each other or importance as a prognostic factor in ovarian cancer is missing. Therefore, we investigated the expressions of VEGF-A, VEGF-C, and VEGF-D, and the receptors VEGFR1, VEGFR2, and VEGFR3 in patients with malignant epithelial ovarian neoplasms. We further compared expression levels between primary tumors and related distant omental metastases.


This study included 86 patients with malignant ovarian epithelial tumors and 16 related distant metastases. Angiogenic factor expression was evaluated using immunohistochemistry (n = 102) and qRT-PCR (n = 29).


Compared to primary high grade serous ovarian tumors, the related omental metastases showed higher expressions of VEGF-A (p = 0.022), VEGF-D (p = 0.010), and VEGFR1 (p = 0.046). In univariate survival analysis, low epithelial expression of VEGF-A in primary tumors was associated with poor prognosis (p = 0.024), and short progression-free survival was associated with high VEGF-C (p = 0.034) and low VEGFR3 (p = 0.002). The relative expressions of VEGF-D, VEGFR1, VEGFR2, and VEGFR3 mRNA determined by qRT-PCR analyses were significantly correlated with the immunohistochemically detected levels of these proteins in primary high grade serous ovarian cancer and metastases (p = 0.004, p = 0.009, p = 0.015, and p = 0.018, respectively).


The expressions of VEGF receptors and their ligands significantly differed between malignant ovarian tumors and paired distant metastases. VEGF-A, VEGF-D, and VEGFR1 protein expressions seem to be higher in distant metastases than in the primary high grade serous ovarian cancer lesions.

Peer Review reports


Among gynecological cancers, ovarian cancer has the worst prognosis and leads to over 150,000 deaths worldwide each year [1]. At the time of diagnosis, 70% of ovarian cancers are disseminated. Standard treatment strategies include cytoreductive surgery and platinum based chemotherapy. Even with optimal therapy, 30–50% of patients die within 5 years after diagnosis [2]. There are presently no curative treatment options for recurrent disease. Prognostic markers for ovarian cancer include tumor stage, residual tumor after primary surgery, the primary response to platinum-based chemotherapy, tumor grade, histological type, and patient age [3].

Improved knowledge of cancer growth and dissemination mechanisms at the molecular level has enabled targeted treatments. In ovarian cancer, targeting endothelial cells of the tumor blood vessels is already standard therapy in both recurrent and primary disease [4, 5]. Additionally, phase I–III studies are presently underway to investigate new antiangiogenic treatment options. Thus, there is an urgent need for biomarkers predicting patient prognosis and responses to specific therapies, to accelerate the development of more personalized treatments [6,7,8].

VEGF (vascular endothelial growth factor) signaling pathways play major roles in tumor angiogenesis and lymphangiogenesis [9, 10]. VEGF-A, -B, -C, and -D and PlGF signal through three tyrosine kinase receptors: VEGFR1, 2, and 3. VEGF-A binds to VEGFR1 and VEGFR2 and is the main stimulator of tumor growth and dissemination. VEGF-A/VEGFR2 signaling is considered the most important pathway for tumor angiogenesis, inducing endothelial permeability and stimulating the accumulation of malignant ascites in ovarian cancer patients [11]. VEGF-C and -D signal mainly through VEGFR3, stimulating lymphangiogenesis and contributing to the generation of angiogenic sprouts [10, 12,13,14,15].

It is essential to evaluate the behavior of angiogenic factors, and to develop new biomarkers for accurate diagnosis and prognosis. We previously reported that alterations of angiogenesis-related circulating proteins—such as elevated angiopoietin-2, high Ang-2/sVEGFR2 ratio, and low level of VEGFR2—predict poor overall survival (OS) and progression-free survival (PFS) in ovarian cancer patients [10]. Earlier studies have shown that VEGF, VEGFR2, and VEGFR3 expressions are related to poor prognosis in endometrial, colorectal, and lung adenocarcinoma [16,17,18]. In ovarian carcinoma, VEGFR1 expression is reportedly associated with shorter PFS [19]. To date, no studies have evaluated all VEGFs and their receptors at the tumor tissue level in both primary tumors and related metastases.

In our present study, we performed immunohistochemical evaluation of VEGF-A, −C, and -D and VEGFR1, R2, and R3 in 86 ovarian carcinoma primary tumors and 16 related metastatic tumors. We measured the expressions of these angiogenic factors using qRT-PCR of 29 tumors. Our aims were to determine whether the expression levels differed between primary tumors and the related metastases, and whether expression of these factors predicted the clinical course, prognosis, and survival of the epithelial ovarian carcinoma patients. To our knowledge, this is the first study to evaluate all VEGFs and their receptors, using both immunohistochemistry and qRT-PCR, in primary ovarian tumors and related metastatic tumors from the same patients, and to correlate the expression data with clinical outcomes.



This study included a total of 86 women who were diagnosed with malignant ovarian epithelial tumors at Kuopio University Hospital between 1999 and 2007. Tumor tissue samples and patient information were prospectively collected at the time of primary diagnosis, and retrospectively analyzed using immunohistochemistry and qRT-PCR. The follow-up time ended in March 2017. In a subgroup of 16 patients with high grade serous ovarian carcinoma, we also analyzed related metastases. The metastatic samples were omental metastases at the time of primary diagnosis. The patient ages ranged from 29 to 88 years (median, 58 years). Histological type and grade were evaluated according to the World Health Organization (WHO) Classification of Tumors [20]. Patients with nonepithelial neoplasms, patients treated prior to operation, unoperated patients (totally 65 patients) were excluded from this study. Epithelial ovarian carcinomas were operatively staged according to the International Federation of Gynecology and Obstetrics (FIGO) criteria. All included patients were treated with platinum-based chemotherapy. Table 1 summarizes the patient characteristics. This study was approved by the Ethical Committee of the Kuopio University Hospital [10].

Table 1 Clinicopathological data of the ovarian cancer patients


Tissue samples were embedded in paraffin and cut into 5-μm-thick sections. Next, the sections were processed for staining with hematoxylin-eosin, VEGF-A (1:250; Santa Cruz), VEGF-C (1:20; Novus), VEGF-D (1:20; Santa Cruz), VEGFR1 (1:15; Santa Cruz), VEGFR2 (1:750; Cell Signaling), and VEGFR3 (11:000, Millipore/Chemicon) [21,22,23,24,25].

During evaluation of immunostained sections, the investigators were blinded to the patients’ clinical status. Immunostaining of VEGF-A, VEGF-C, VEGF-D, VEGFR1, VEGFR2, and VEGFR3 was microscopically evaluated (Leitz Wetzlar 512,761/20; Germany). After screening the whole section, ten randomly selected microscopic fields were examined at × 200 magnification (× 20 objective lens and × 10 ocular lens; 0.74 mm2 per field) to determine the mean percentage of specific immunostained epithelial tumor cells. The percentage of positively stained cells (PP) was assigned as a numerical score: 0, negative; 1, < 10%; 2, 11–50%; 3, 51–80%; and 4, > 80% positive cells. The intensity (SI) of the immunostained areas was defined as follows: 0, negative; 1, weak; 2, moderate; 3, strong. An immunoreactive score (IRS) ranging from 0 to 12 was calculated using the following formula: IRS = PP × SI [26].

The percentage of strong epithelial staining was separately determined for the VEGFR1, 2, and 3 sections. Stromal immunostaining was evaluated by ranking the staining intensity: 0, negative; 1, weak; 2, strong. Endothelial expression was deemed positive or negative, and we also evaluated the cellular distribution of the staining (cytoplasm, nucleus, and cell membrane) in epithelial ovarian tumors. Stromal expression was positive in all samples, and thus the staining intensity was classified as weak or strong.

Quantitative real-time polymerase chain reaction (qRT-PCR) analysis

We performed qRT-PCR analysis of samples from 13 primaries and 16 metastatic tumors.

RNA was isolated using TRI-reagent (Sigma Aldrich, St. Louis, Missouri, USA), and treated with DNase (Promega, Fitchburg, Wisconsin, USA). From 5 μg total RNA, cDNA was synthesized using random hexamer primers (Promega) and RevertAid™ reverse transcriptase (Fermentas, Waltham, Massachusetts, USA). We measured the relative expressions of the mRNAs encoding VEGF-A, VEGF-C, VEGF-D, VEGFR1, VEGFR2, and VEGFR3 by qRT-PCR using specific Assays-on-Demand target mixes (Applied Biosystems, Foster City, California, USA), following the manufacturer’s protocol (StepOnePlus, Applied Biosystems). The expression levels were normalized to peptidylprolyl isomerase A [25, 27].

Statistical analysis

Statistical analyses were performed using SPSS for Windows (version 24, 1989–2016, SPSS Inc., Chicago, USA). We performed the Kruskal-Wallis test followed by the Mann-Whitney test with multiple comparisons when appropriate. For analyses of clinicopathological associations and survival analyses, PP values, IRS values, and PCR levels were each dichotomized into low and high groups using the median values as a cut-off [10]. We used the chi-squared test to analyze frequency tables, and the Wilcoxon signed-rank test to compare histology and qRT-PCR results between primary ovarian tumors and the related metastases. Univariate survival analyses were based on the Kaplan-Meier method, and survival curves were compared using the log-rank test. Variables that were significant in univariate analyses were entered in a stepwise manner into the Cox regression multivariate analysis. OS was defined as the time interval between the date of surgery and the date of death or the end of follow-up. PFS was defined as the time interval between the date of surgery and the date of identified recurrence or death. Correlations between histological parameters and qRT-PCR levels were analyzed using the Spearman’s correlation test. A p value of < 0.05 was considered significant [10].


Immunohistochemical analyses of VEGF-A, VEGF-C, VEGF-D, VEGFR1, VEGFR2, and VEGFR3 in primary high grade serous and endometrioid ovarian tumors

VEGF-A, VEGF-C, and VEGF-D were mainly detected in the cytoplasm of tumor epithelial cells. Of these VEGFs, VEGF-C showed the strongest expression in epithelial primary high grade serous and endometrioid ovarian cancer cells. VEGF-C staining intensity among epithelial cancer cells was moderate to strong in 89% of primary tumors (IRS, 8), while moderate-to-strong VEGF-A staining intensity was found in 86% of primary tumors (IRS, 6), and VEGF-D staining was at least moderate in 56% of primary tumors (IRS, 4). Of all evaluated factors, epithelial immunostaining of VEGF-D was the weakest. In primary high grade serous and endometrioid tumors, strong stromal staining was observed in 40% of VEGF-C-stained samples, 41% of VEGF-A-stained samples, and 21% of VEGF-D-stained samples. Stromal staining was most prominent in the stromal fibroblasts and endothelial staining was positive in all samples (Fig. 1).

Fig. 1
figure 1

Diagram of the immunohistochemical parameters of primary tumors and metastases. a) VEGFs and b) VEGF-receptors. IRS immunoreactive score, PP percentage of positively stained cells, SI staining intensity

VEGF receptors (R1, R2, R3) were mainly expressed in the epithelial cell membranes of primary tumors, in stromal tissue, and in vascular endothelial cells (Fig. 2). For VEGFR1, VEGFR2, and VEGFR3, respectively, the median IRS values were 6, 3, and 4; and the median PP values for epithelial cells were 75, 60, and 60%. Epithelial expression of VEGFR1 was strongest, while VEGFR2 was weakest. We detected strong stromal staining in 37% of VEGFR1-stained samples, 44% of VEGFR2-stained samples, and 65% of VEGFR3-stained samples in high grade serous and endometrioid cancer. Vascular endothelial staining was notably strong for all receptors in most samples (Fig. 2). Immunohistochemical analysis of all histologies and grades is presented in Table 2 and Fig. 1. There was no significant difference when high grade serous and endometrioid tumors were analyzed separately compared to the whole study population.

Fig. 2
figure 2

Images show VEGFR2 staining in a primary tumor samples (a) and metastasis sample (b), and VEGFR3 staining in a primary tumor samples (c) and metastasis sample (d). VEGF-receptors were mainly expressed in the epithelial cell membrane (a and c, black arrows), in stromal tissue and in vascular endothelial cells (b and d, black arrows)

Table 2 Expressions of VEGFs and their receptors in primary ovarian carcinomas and related metastases
Fig. 3
figure 3

Compared to primary tumors, the related omental metastases showed stronger expressions of VEGF-A (a and b), VEGF-D (c and d), and VEGFR1 (g and h). Primary tumors are shown on the left, and metastases in the middle column. The right column presents the mRNA levels of the analogous angiogenic factors, as determined by qRT-PCR. e and f VEGF-C expression was higher in primary tumors than in metastases

In primary ovarian tumors, VEGFR3 showed moderate but significant correlations with its ligands lymphangiogenic VEGF-C (r = 0.4, p < 0.0005) and VEGF-D (r = 0.4, p = 0.001). Similarly, VEGFR2 was significantly correlated with VEGF-D (r = 0.4, p < 0.0005). In primary tumors, VEGF-A was weakly correlated with VEGFR2 (r = 0.2, p = 0.03), but was not significantly correlated with VEGFR1. We observed only weak correlations among VEGFs, while VEGF receptors were moderately correlated with each other in primary tumors. In high grade tumors the correlations did not differ significantly from the presented ones. (data not shown).

Immunohistochemical expression of VEGF-A, VEGF-C, VEGF-D, VEGFR1, VEGFR2, and VEGFR3 in metastatic lesions of high grade serous ovarian cancer

In metastatic epithelium and stroma, the tumor tissue distribution of VEGF-A, VEGF-C, and VEGF-D expression paralleled the expression pattern in the primary cancer (Fig. 1). Of all VEGFs, VEGF-A expression was the strongest in tumor epithelium. VEGF-D was expressed in a wide area, but its staining intensity was weak. Among the growth factors, VEGF-A and VEGF-D dominated the stromal expression in metastases. (Table 2, Fig. 1). As in the primary tumors, metastases showed strong growth factor staining in vascular endothelium.

In metastases, we observed strong receptor staining in stromal tissue and vascular endothelium (Fig. 2). VEGFR1 exhibited the strongest epithelial cell expression, and VEGFR2 the weakest expression in ovarian cancer metastases (Table 2, Fig. 1). VEGF-A was significantly correlated with VEGFR1 in metastases (r = 0.502, p = 0.048).

Angiogenic factor expression in primary high grade tumors compared to the related metastases

VEGF-A and VEGF-D expressions in epithelial tumor cells were significantly stronger in omental metastases than in the related primary tumors. Compared to primary tumors, metastases showed a significantly higher median PP of VEGF-A [85% (95% CI, 63–87%) versus 70% (95% CI, 41–73%); p = 0.022], significantly higher median PP of VEGF-D [85% (95% CI, 60–90%) versus 60% (95% CI, 35–61%); p = 0.010], and significantly higher median IRS of VEGF-D [7 (95% CI, 5–9) versus 4 (95% CI, 3–6); p = 0.013]. Additionally, VEGFR1 staining intensity was significantly stronger in metastatic epithelial tumor cells (moderate-to-strong intensity in 88% of samples) than in the primary tumors (moderate-to-strong intensity in 74% of samples) (p = 0.046). In contrast to the VEGF-A and VEGF-D expressions, epithelial expression of VEGF-C was higher in the primary tumors than in the metastases (p = 0.010). Differences in other factors did not reach statistical significance (Table 2, Figs. 1, 3).

qRT-PCR levels and correlations with immunohistochemical expression

In primary high grade tumors, we found a significant correlation between VEGF-D IRS and VEGF-D mRNA (correlation coefficient, r = 0.991, p = 0.004), VEGFR1 PP and VEGFR1 mRNA (r = 0.967, p = 0.014) and VEGFR3 SI and VEGFR3 mRNA (r = − 0.655, p = 0.04).

In omental metastases, the immunohistochemical expressions of all receptors correlated to the corresponding mRNA levels. Specifically, VEGFR1 mRNA level was significantly correlated with VEGFR1 IRS (r = 0.717; p = 0.009) and VEGFR1 PP (r = 0.780; p = 0.003); VEGFR2 mRNA was significantly correlated with VEGFR2 IRS (r = 0.679; p = 0.015) and VEGFR2 PP (r = 0.810; p = 0.001); and VEGFR3 mRNA was significantly correlated with VEGFR3 IRS (r = 0.667; p = 0.018) and VEGFR3 PP with (r = 0.661; p = 0.019). The VEGFR2 mRNA level was correlated with the percentage of VEGFR2 strong epithelial staining (r = 0.726; p = 0.011). The VEGFR3 mRNA level was correlated with the percentage of VEGFR3 strong epithelial staining (r = 0.854; p = 0.007), and with VEGFR3 stromal intensity (r = 0.599, Pearson’s test p = 0.04).

The mRNA levels of VEGF-A, VEGF-C, VEGF-D, VEGFR1, −R2 and R3 did not reach statistical significance when primary tumors were compared to the related metastases.

Relation of clinicopathological data to the expression of VEGFs and VEGF receptors

Ovarian cancer recurrence was associated with low expression of VEGF-A (p = 0.005) and VEGFR1 (p = 0.01). The high histological grade was related to low expression of VEGF-C (p = 0.006), VEGF-D (p = 0.04), VEGFR1 (p < 0.001), VEGFR2 (p = 0.008) and VEGFR3 (p = 0.004) in epithelial cells of the primary ovarian cancer.

When high grade serous and endometrioid tumors were analyzed as one group, low expression of VEGF-A (p = 0.010) and high expression of VEGFR3 (p = 0.037) were associated with the recurrence of cancer. Low VEGF-C (p = 0.026) was related to no residual tumor at primary surgery in high grade tumors.

In the whole study population, VEGF expression was related to the histological subtype. Strong expressions of VEGF-C, VEGFR1, VEGFR2, and VEGFR3 were most commonly found in tumors with mucinous histology, while the weakest expressions were observed in clear cell and endometrioid tumors (p = 0.033, p < 0.001, p = 0.027, p = 0.005, respectively). Of the studied parameters, serous tumors exhibited the largest variation.

Other clinical parameters were not significantly associated with the immunohistochemical data. We obtained parallel results when analyzing associations with Pearson’s chi-squared test and as continuous parameters with the Kruskal-Wallis test (Table 3).

Table 3 Relationship between immunohistochemical results and clinicopathological data

Overall survival and progression-free survival of the ovarian cancer patients

The median follow-up time was 65 months (range, 0–198 months). At the end of follow-up, 51 patients (59%) had died. The median OS was 77 months (95% CI, 51–102 months) and the 5-year survival rate was 57% (95% CI, 46–68%). In univariate analysis, low VEGF-A expression in the primary tumors predicted poor OS (p = 0.024). Among serous primary tumors, low VEGF-A PP was associated with short OS (p = 0.005) as well as in high grade serous and endometrioid tumors (p = 0.042) (Fig. 4). Univariate survival analysis revealed that poor OS was significantly predicted by advanced stage, presence of residual primary tumor, incomplete primary response to chemotherapy, and presence of ascites. In multivariate analysis, OS remained significantly associated with the presence of residual primary tumor (p < 0.001) and incomplete response to chemotherapy (p < 0.0001) (Table 4).

Fig. 4
figure 4

Univariate analysis of immunohistochemical biomarkers as prognostic factors in ovarian cancer patients. Poor overall survival (OS) was significantly associated with low VEGF-A immunoreactive score (IRS) (a) and with low VEGF-A percentage of positively stained cells (PP) in serous carcinoma (b). Short progression-free survival (PFS) was significantly associated with high VEGF-C PP (c) and low VEGFR3 staining intensity (SI) (d)

Table 4 Overall survival and progression-free survival according to immunohistochemical staining of VEGFs and receptors and clinical characteristics

PFS analysis included a total of 67 patients, of whom 42 (62%) experienced recurrence. The median PFS was 14 months (95% CI, 9–18 months). Short PFS was significantly associated with strong stromal staining of VEGF-C in primary tumors, in both univariate analysis (p = 0.036) and multivariate analysis (p = 0.029). In high grade serous and endometrioid tumors strong stromal VEGF-C was also related to short PFS (p = 0.021). Among high grade tumors, shorter PFS was associated with high VEGF-C (p = 0.034) and low VEGFR3 SI in epithelial cells of the primary tumor (p = 0.006). In both univariate and multivariate analysis, shorter PFS was significantly predicted by advanced stage (p = 0.022 in multivariate analysis) and the presence of ascites (p = 0.045 in multivariate analysis) (Table 4).


This is the first study to include the comprehensive analysis of a large panel of VEGFs and their receptors in the same patient population as well comparing primaries to the metastatic lesions, and link the results to the clinical outcome of the ovarian cancer patients. We found that expressions of VEGF-A, VEGF-D, and VEGFR1 were higher in omental metastases than in the related primary high grade serous tumors. Furthermore, poor OS was most strongly predicted by low VEGF-A expression in the primary high grade tumors, and short PFS was predicted by high VEGF-C expression.

Compared to primary tumors, the related ovarian cancer metastases showed increased expressions of the pro-angiogenic VEGF-A, VEGF-D, and VEGFR1 proteins. Metastatic tumors seemed to exhibit accumulation of these factors, with stronger expression in stromal tissue and vessel endothelium, enabling neovascularization and lymphangiogenesis. Moreover, VEGFR2, which binds VEGF-A, showed weaker expression in metastatic epithelial cells and stronger expression in stromal tissue and vascular endothelium. This is in line with previous findings that VEGFR2 and VEGFR3 are expressed on vascular endothelium, but to a lesser extent in malignant cells, in 13 common human solid cancer types, including ovarian cancer [24].

This result supports our previous finding in a smaller population that VEGF-A is present at higher levels in omental metastasis than in the primary tumors [25]. Gadducci et al. also demonstrated increased VEGF-A immunostaining in peritoneal metastatic lesions of ovarian cancer, but without showing prognostic significance or relation to clinical outcome [28].

Compared to the primary high grade serous tumors, related metastases showed lower expression of lymphangiogenic VEGF-C, and similar VEGFR3 expression. These results may suggest that angiogenesis is more active in distantly disseminated cancer metastasis, whereas lymphangiogenic VEGF-C predominantly acts more locally via paracrine mechanisms.

The findings of our present study, indicate that VEGF-A might be the most promising angiogenic factor for clinical use as a prognostic marker. Our data showed that low IRS and PP values for VEGFs and their receptors in high grade tumor epithelial cells were associated with ovarian carcinoma recurrence, and thus with more aggressive cancer. Low VEGF-A expression also predicted a short OS in the whole cohort, as well as in analyses considering only high grade serous carcinomas. These results might be explained by changes in the tumor microenvironment, hypoxia, and metabolic reprogramming. Hypoxia-inducible factor 1 (HIF1) is the main regulator of cellular responses to hypoxia, and induces angiogenesis. This regulation mechanism might not function equally in all tumors. Notably, strong cytoplasmic p53 expression has been associated with non-angiogenic malignant lung tumors [29, 30]. VEGF distribution could also vary between tumor tissues and circulation. Interestingly, we found that the lymphangiogenic factors VEGF-C and VEGFR3 were significantly correlated with PFS.

Earlier studies have yielded contradictory and heterogeneous results regarding angiogenic factor expression and prognosis in ovarian cancer. Nishida et al. concluded that high tumor expression of VEGF-C and VEGFR2 reflects ovarian carcinoma spread and poor prognosis [31], while Yokoyama et al. found that VEGF-D expression predicts poor OS in ovarian carcinoma [32]. In contrast, Engels et al. reported shorter PFS in cases of VEGF-negative ovarian tumors [26], and Shen et al. demonstrated a better OS among patients with low-VEGF-expressing ovarian tumors [33]. Additionally, Wimberg et al. found only low rates (18–43%) of positive VEGF receptor expression in primary ovarian cancer tissue [19]. Discrepancies between studies may be due to differences in analysis and utilized antibodies, or could be related to the specific carcinoma investigated, different tumor types, or patient groups with variable treatments.

The presently observed variation between histological subtypes of ovarian cancer may indicate that the VEGF pathway plays distinct biologic roles in these different groups. Findings in lung, breast, and gastric cancers suggest that the levels of VEGF autocrine loop activity have different impacts on distinct histologic cell phenotypes, potentially influencing the clinical outcome [18]. In our study, VEGF-C and VEGF receptors were most highly expressed in mucinous tumors. However, the number of mucinous carcinomas was low, as in the general population. Tumors with serous histology exhibited the largest variation; with well-differentiated papillary tumors seeming to show stronger expression of these factors. However, the statistics did not differ significantly between the total population and the high grade serous tumors.

In our present study, expression levels were semiquantitatively scored based on staining intensity and percentage of stained tumor cells, using the immunoreactive score as described [26]. This formula has been applied in various studies with slight modifications, since proposed by Stegner in 1986 [24, 34,35,36,37]. To elucidate the nuances behind the IRS, we also examined the separate factors PP and SI. The IRS has the advantage of being a single combining value, enabling comparison with other studies using the same formula. On the other hand, as a single value, it cannot accurately indicate all contributing factors, and the score formulation can vary among different studies [34, 37].

As a limitation of this study, the Western blot analysis was not performed. Antibodies used in this study were already tested with the Western blot by the vendors. We believe that with immunohistochemical stainings, we were able to detect the expression of growth factors and their receptors reliably, and analyse the exact location of these proteins in tumor tissue. In addition, qRT-PCR was performed to confirm the results in RNA level. However, in the future directions the Western blot analyses are needed.

There are not yet any tumor angiogenic biomarkers in clinical use, but antiangiogenic therapy development has raised an urgent need for prognostic and predictive biomarkers to improve patient selection, follow-up, and treatment effect monitoring. The first antiangiogenic agent to provide clinical benefit was the humanized monoclonal VEGF antibody bevacizumab. Since then, various VEGF neutralizing antibodies, multi-target tyrosine kinase inhibitors of the VEGF signaling pathway, and soluble VEGF receptors have entered clinical use, although none of them is yet in routine clinical use in ovarian cancer [4, 5, 7]. Our patient samples are from the era when bevacizumab was not in the clinical use.

VEGF inhibition can reduce the number of tumor vessels by impeding vascular sprouting and impairing endothelial survival, leading to vessel regression and thus normalizing the remaining vessel morphology. Our present results showed stronger VEGF-A, VEGF-D, and VEGFR1 expressions in related metastases, supporting the concept of using antiangiogenic treatments in disseminated carcinoma and for treating patients with incomplete primary surgical resection [4, 5, 32].


Our results indicate that VEGF-A seems to be the most promising angiogenic marker in high grade serous ovarian cancer. VEGF-A expression was stronger in metastatic tumors, and low VEGF-A expression in primary tumors was associated with poor prognosis. In terms of lymphangiogenic factors, high VEGF-C expression and low VEGFR3 predicted shorter PFS. To improve antiangiogenic treatment efficacy, it would be beneficial to block both the angiogenic and lymphangiogenic pathways in ovarian cancer. Our present findings further suggest that it would be advantageous to examine a metastatic sample in addition to the primary tumor sample, to reveal the angiogenic status and improve the results of treatment with antiangiogenic therapies. Our data support the present and developing antiangiogenic treatment strategies, and suggest pathways for further investigation of angiogenic biomarkers.

Availability of data and materials

The datasets generated and analyzed in the current study are available from the corresponding author upon reasonable request.





Basic fibroblast growth factor


Colorectal cancer


Hypoxia-inducible factor


Immunoreactive score


Messenger ribonucleic acid




Overall survival


Platelet-derived growth factor


Progression-free survival


Percentage of positively stained cells


Quantitative real-time polymerase chain reaction


Staining intensity


Tumor growth factor β


  1. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.

    Article  PubMed  Google Scholar 

  2. du Bois A, Luck H-J, Meier W, Adams H-P, Mobus V, Costa S, Bauknecht T, Richter B, Warm M, Schroder W, Olbricht S, Nitz U, Jackisch C, Emons G, Wagner U. A randomized clinical trial of cisplatin/paclitaxel versus carboplatin/paclitaxel as first-line treatment of ovarian Cancer. J Natl Cancer Inst. 2003;95(17):1320–9.

    Article  PubMed  Google Scholar 

  3. Wimberger P, Lehmann N, Kimmig R, Burges A, Meier W, Du Bois A. Prognostic factors for complete debulking in advanced ovarian cancer and its impact on survival. An exploratory analysis of a prospectively randomized phase III study of the Arbeitsgemeinschaft Gynaekologische Onkologie ovarian Cancer study group (AGO-OVAR). Gynecol Oncol. 2007;106(1):69–74.

    Article  PubMed  Google Scholar 

  4. Burger RA, Brady MF, Bookman MA, Fleming GF, Monk BJ, Huang H, Mannel RS, Homesley HD, Fowler J, Greer BE, Boente M, Birrer MJ, Liang SX. Incorporation of bevacizumab in the primary treatment of ovarian cancer. N Engl J Med. 2011;365(26):2473–83.

    Article  CAS  PubMed  Google Scholar 

  5. Perren TJ, Swart AM, Pfisterer J, Ledermann JA, Pujade-Lauraine E, Kristensen G, Carey MS, Beale P, Cervantes A, Kurzeder C, du Bois A, Sehouli J, Kimmig R, Stähle A, Collinson F, Essapen S, Gourley C, Lortholary A, Selle F, Mirza MR, Leminen A, Plante M, Stark D, Qian W, Parmar MKB, Oza AM. A phase 3 trial of bevacizumab in ovarian cancer. N Engl J Med. 2011;365(26):2484–96.

    Article  CAS  PubMed  Google Scholar 

  6. Sawyers CL. The cancer biomarker problem. Nature. 2008;452(7187):548–52.

    Article  CAS  PubMed  Google Scholar 

  7. Monk BJ, Minion LE, Coleman RL. Anti-angiogenic agents in ovarian cancer: past, present, and future. Ann Oncol. 2016;27(suppl_1):i39.

    Google Scholar 

  8. Abdalla DR, Simoens C, Bogers J-P, Murta EFC, Marcia A. Michelin: angiogenesis markers in gynecological tumors and patents for anti- Angiogenic approach: review. Recent Patents on Anti-Cancer Drug Discovery. 2015;10(3):298–307.

    Article  CAS  PubMed  Google Scholar 

  9. Saharinen P, Eklund L, Pulkki K, Bono P, Alitalo K. VEGF and angiopoietin signaling in tumor angiogenesis and metastasis. Trends Mol Med. 2011;17(7):347–62.

    Article  CAS  PubMed  Google Scholar 

  10. Sallinen H, Heikura T, Koponen J, Kosma V-M, Heinonen S, Ylä-Herttuala S, Anttila M. Serum angiopoietin-2 and soluble VEGFR-2 levels predict malignancy of ovarian neoplasm and poor prognosis in epithelial ovarian cancer. BMC Cancer. 2014;14(1):696.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Lohela M, Bry M, Tammela T, Alitalo K. VEGFs and receptors involved in angiogenesis versus lymphangiogenesis. Curr Opin Cell Biol. 2009;21(2):154–65.

    Article  CAS  PubMed  Google Scholar 

  12. Nagy JA, Masse EM, Herzberg KT, Meyers MS, Yeo K, Yeo T, Sioussat TM, Dvorak HF. Pathogenesis of ascites tumor growth: vascular permeability factor, vascular Hyperpermeability, and ascites fluid accumulation. Cancer Res. 1995;55(2):360.

    CAS  PubMed  Google Scholar 

  13. Joukov V, Pajusola K, Kaipainen A, Chilov D, Lahtinen I, Kukk E, Saksela O, Kalkkinen N, Alitalo K. A novel vascular endothelial growth factor, VEGF-C, is a ligand for the Flt4 (VEGFR-3) and KDR (VEGFR-2) receptor tyrosine kinases. EMBO J. 1996;15(2):290–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Achen MG, Jeltsch M, Kukk E, Makinen T, Vitali A, Wilks AF, Alitalo K, Stacker SA. Vascular endothelial growth factor D (VEGF-D) is a ligand for the tyrosine kinases VEGF receptor 2 (Flk1) and VEGF receptor 3 (Flt4). Proc Natl Acad Sci U S A. 1998;95(2):548–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Wirzenius M, Alitalo K, Hellström M, Peltonen R, Suchting S, Pytowski B, Zarkada G, Laakkonen P, Duarte A, Ylä-Herttuala S, Isoniemi H, Murtomäki A, Eichmann A, Betsholtz C, Christofori G, Wallgard E, Freitas C, Waltari M, Shibuya M, Schomber T, Tammela T. Blocking VEGFR-3 suppresses angiogenic sprouting and vascular network formation. Nature. 2008;454(7204):656–60.

    Article  PubMed  Google Scholar 

  16. Wang J, Taylor A, Showeil R, Trivedi P, Horimoto Y, Bagwan I, Ewington L, Lam EW, El-Bahrawy MA. Expression profiling and significance of VEGF-A, VEGFR2, VEGFR3 and related proteins in endometrial carcinoma. Cytokine. 2014;68(2):94–100.

    Article  CAS  PubMed  Google Scholar 

  17. Martins SF, Garcia EA, Luz MAM, Pardal F, Rodrigues M, Filho AL. Clinicopathological correlation and prognostic significance of VEGF-A, VEGF-C, VEGFR-2 and VEGFR-3 expression in colorectal Cancer. Cancer Genomics Proteomics. 2013;10(2):55–68.

    CAS  PubMed  Google Scholar 

  18. Pajares MJ, Agorreta J, Larrayoz M, Vesin A, Ezponda T, Zudaire I, Torre W, Lozano MD, Brambilla E, Brambilla C, Wistuba II, Behrens C, Timsit J-F, Pio R, Field JK, Montuenga LM. Expression of tumor-derived vascular endothelial growth factor and its receptors is associated with outcome in early squamous cell carcinoma of the lung. J Clin Oncol. 2012;30(10):1129–36.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Wimberger P, Chebouti I, Kasimir-Bauer S, Lachmann R, Kuhlisch E, Kimmig R, Süleyman E, Kuhlmann JD. Explorative investigation of vascular endothelial growth factor receptor expression in primary ovarian cancer and its clinical relevance. Gynecol Oncol. 2014;133(3):467.

    Article  CAS  PubMed  Google Scholar 

  20. Hanby AM, Walker C, Tavassoli FA, Devilee P: Pathology and Genetics: Tumours of the Breast and Female Genital Organs. WHO Classification of Tumours series - volume IV. Lyon: IARC Press 2004, 6(3). Breast Cancer Res.

  21. Sahoo Subhransu S, Lombard Janine M, Yvette I, Rachel O’S, Wood Lisa G, Pravin N, Kenneth J, Tanwar Pradeep S. Adipose-derived VEGF–mTOR signaling promotes endometrial hyperplasia and Cancer: implications for obese women. Mol Cancer Res. 2018;16(2):309–21.

    Article  CAS  PubMed  Google Scholar 

  22. Dhar S, Kumar A, Zhang L, Rimando AM, Lage JM, Lewin JR, Atfi A, Zhang X, Levenson AS. Dietary pterostilbene is a novel MTA1-targeted chemopreventive and therapeutic agent in prostate cancer. Oncotarget. 2016;7(14):18469.

    Article  PubMed  PubMed Central  Google Scholar 

  23. MacDonald DD, Jensen KK, Kerns SA, Lyden D, Zhu Z, Wu Y, Riba RD, Shido K, Rafii S, Kaplan RN, Port ER, Jin DK, Ruggero D, Hicklin D, Port JL, Vincent L, Altorki N, Costa C, Shmelkov SV, Bramley AH, Zacharoulis S. VEGFR1-positive haematopoietic bone marrow progenitors initiate the pre-metastatic niche. Nature. 2005;438(7069):820–7.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Smith NR, Baker D, James NH, Ratcliffe K, Jenkins M, Ashton SE, Sproat G, Swann R, Gray N, Ryan A, Jürgensmeier JM, Womack C. Vascular endothelial growth factor receptors VEGFR-2 and VEGFR-3 are localized primarily to the vasculature in human primary solid cancers. Clin Cancer Res. 2010;16(14):3548–61.

    Article  CAS  PubMed  Google Scholar 

  25. Lindgren A, Anttila M, Rautiainen S, Arponen O, Kivelä A, Mäkinen P, Härmä K, Hämäläinen K, Kosma V, Ylä-Herttuala S, Vanninen R, Sallinen H. Primary and metastatic ovarian cancer: characterization by 3.0T diffusion-weighted MRI. Eur Radiol. 2017;27(9):4002–12.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Engels K, du Bois A, Harter P, Fisseler-Eckhoff A, Kommoss F, Stauber R, Kaufmann M, Nekljudova V, Loibl S. VEGF-A and i-NOS expression are prognostic factors in serous epithelial ovarian carcinomas after complete surgical resection. J Clin Pathol. 2009;62(5):448–54.

    Article  CAS  PubMed  Google Scholar 

  27. Heinonen SE, Kivelä AM, Huusko J, Dijkstra MH, Gurzeler E, Mäkinen PI, Leppänen P, Olkkonen VM, Eriksson U, Jauhiainen M, Ylä-Herttuala S. The effects of VEGF-A on atherosclerosis, lipoprotein profile, and lipoprotein lipase in hyperlipidaemic mouse models. Cardiovasc Res. 2013;99(4):716–23.

    Article  CAS  PubMed  Google Scholar 

  28. Gadducci A, Viacava P, Cosio S, Cecchetti D, Fanelli G, Fanucchi A, Teti G, Genazzani AR. Vascular endothelial growth factor (VEGF) expression in primary tumors and peritoneal metastases from patients with advanced ovarian carcinoma. Anticancer Res. 2003;23(3C):3001–8.

    CAS  PubMed  Google Scholar 

  29. Casillas Andrea L, Toth Rachel K, Sainz Alva G. Singh Neha, Desai Ankit a, Kraft Andrew S, Warfel noel a: hypoxia-inducible PIM kinase expression promotes resistance to antiangiogenic agents. Clin Cancer Res. 2018;24(1):169–80.

    Article  CAS  PubMed  Google Scholar 

  30. Adighibe O, Leek RD, Fernandez-Mercado M, Hu J, Snell C, Gatter KC, Harris AL, Pezzella F. Why some tumours trigger neovascularisation and others don’t: the story thus far. 癌症:英文版. 2016;35(2):69–79.

    Google Scholar 

  31. Nishida N, Yano H, Komai K, Nishida T, Kamura T, Kojiro M. Vascular endothelial growth factor C and vascular endothelial growth factor receptor 2 are related closely to the prognosis of patients with ovarian carcinoma. Cancer. 2004;101(6):1364–74.

    Article  CAS  PubMed  Google Scholar 

  32. Yokoyama Y, Charnock-Jones DS, Licence D, Yanaihara A, Hastings JM, Holland CM, Emoto M, Umemoto M, Sakamoto T, Sato S, Mizunuma H, Smith SK. Vascular endothelial growth factor-D is an independent prognostic factor in epithelial ovarian carcinoma. Br J Cancer. 2003;88(2):237–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Shen GH, Ghazizadeh M, Kawanami O, Shimizu H, Jin E, Araki T, Sugisaki Y. Prognostic significance of vascular endothelial growth factor expression in human ovarian carcinoma. Br J Cancer. 2000;83(2):196–203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mattern J, Koomägi R, Volm M. Association of vascular endothelial growth factor expression with intratumoral microvessel density and tumour cell proliferation in human epidermoid lung carcinoma. Br J Cancer. 1996;73(7):931–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kulbe H, Thompson R, Wilson JL, Robinson S, Hagemann T, Fatah R, Gould D, Ayhan A, Balkwill F. The inflammatory cytokine tumor necrosis factor- alpha generates an autocrine tumor-promoting network in epithelial ovarian cancer cells. Cancer Res. 2007;67(2):585–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Szubert S, Moszynski R, Michalak S, Nowicki M, Sajdak S, Szpurek D. The associations between serum VEGF, bFGF and endoglin levels with microvessel density and expression of proangiogenic factors in malignant and benign ovarian tumors. Microvasc Res. 2016;107:91–6.

    Article  CAS  PubMed  Google Scholar 

  37. Mehta R, Kyshtoobayeva A, Kurosaki T, Small EJ, Kim H, Stroup R, McLaren CE, Li K-T, Fruehauf JP. Independent association of angiogenesis index with outcome in prostate cancer. Clin Cancer Res. 2001;7(1):81.

    CAS  PubMed  Google Scholar 

Download references


We thank Mrs. Helena Kemiläinen for skillful technical assistance.


This study was supported by The Finnish Medical Foundation, The Ida Montin Foundation and The Emil Aaltonen Foundation in analysis, interpretation of the data and writing the manuscript and a Kuopio University Hospital VTR grant in collecting the data.

Author information

Authors and Affiliations



MS wrote the manuscript, and analyzed the tissue samples and statistics. HS and MA designed the study, collected the samples and clinical data, reviewed the data, and revised the manuscript. KH analyzed samples. AK performed qRT-PCR. SYH and VMK provided the conditions, materials, and laboratory services for the study. LKN revised the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Hanna Sallinen.

Ethics declarations

Ethics approval and consent to participate

Ethical approval of this study was obtained from the Ethics Committee of Kuopio University Hospital. All patients provided written informed consent before tissue sample collection. This study was conducted in accordance with the ethical principles of the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sopo, M., Anttila, M., Hämäläinen, K. et al. Expression profiles of VEGF-A, VEGF-D and VEGFR1 are higher in distant metastases than in matched primary high grade epithelial ovarian cancer. BMC Cancer 19, 584 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: