Skip to main content


  • Research
  • Open Access

Expression of fibroblast growth factor receptor family members is associated with prognosis in early stage cervical cancer patients

Contributed equally
Journal of Translational Medicine201614:124

  • Received: 10 November 2015
  • Accepted: 25 April 2016
  • Published:



The oncogenic role of the fibroblast growth factor receptor (FGFR) has been recognized in a number of different cancer types. However, the prognostic significance of FGFRs has not been elucidated yet in cervical cancer. In the present study, we investigate the expression of FGFRs and their prognostic value in cervical cancer patients.


FGFR1, FGFR2, FGFR3, and FGFR4 expression was determined by immunohistochemistry in conjunction with quantitative digital image analysis of 336 formalin-fixed, paraffin-embedded cervical cancer tissues and 61 normal cervical tissues, as well as NCI60 cell microarray. Subsequently, the association between clinicopathological characteristics and patient survival was assessed.


FGFRs proteins were differentially expressed in the NCI60 cell line panel and showed considerable correlation between protein and mRNA expression. The expression of FGFR1, FGFR2, and FGFR4 were higher in cancer tissues than in normal tissues, whereas the expression of FGFR3 was higher in normal tissues. FGFR1 was highly expressed in adeno-/adenosquamous carcinoma (P = 0.020), while FGFR2, FGFR3, and FGFR4 expression were more prominent in squamous cell carcinoma (P < 0.001, P < 0.001, and P = 0.020, respectively). FGFR2 expression was significantly higher in small sized tumors (P = 0.020). Additionally, high FGFR2 and FGFR4 were correlated with negative lymph node metastasis (P = 0.048 and P = 0.040, respectively). FGFR1, FGFR2, and FGFR3 were highly expressed in tumors without parametrial involvement (P = 0.030, P = 0.005, and P = 0.010, respectively). In survival analysis, high expressions of FGFR2, FGFR3, and FGFR4 was associated with longer disease-free survival (P = 0.006, P = 0.035, P = 0.001, respectively) and overall survival (P = 0.003, P = 0.002, P = 0.003, respectively). Notably, the co-expression of all three FGFRs was significantly associated with favorable disease-free survival (P < 0.001) and overall survival (P < 0.001), compared to the negative expressions of the three FGFRs. The prognostic significance persisted in the cox regression analysis.


The frequent expression of members of the FGFR family in cervical cancer suggests they may have prognostic and therapeutic relevance.


  • Fibroblast growth factor receptor
  • Image analysis
  • Immunohistochemistry
  • Prognosis
  • Survival analysis
  • Uterine cervical neoplasm


Cervical cancer is the third most common type of cancer among women worldwide, and is the most prevalent female malignancy in many developing countries [1, 2]. Although preventive vaccination and screening are good in prevention, invasive cancer continues to occur, even among women who have access to cancer screening, and the prognosis remains poor in patients with a bulky tumor or adenocarcinoma histology [35]. Clinical factors such as International Federation of Gynecology and Obstetrics (FIGO) stage, lymph node metastasis, and tumor size may serve as prognostic markers, however they are insufficient in accurately predicting recurrence and survival. Thus, biomarkers, including molecular markers, are needed. Patient care would be considerably improved if tumor behavior could be reliably prognosticated at the time of initial diagnosis.

The mammalian fibroblast growth factor (FGF) family comprises 18 ligands that exert their actions through four transmembrane tyrosine kinase receptors; fibroblast growth factor receptor (FGFR) 1, FGFR2, FGFR3, and FGFR4. The FGF/FGFR system is responsible for multiple cellular functions, including proliferation, differentiation, survival, and motility [69]. Considering these functions, it is no surprise they are susceptible to aberration in cancer cells [10]. Several types of genetic alterations have been recognized in cancers, including gene amplifications, activating mutations, chromosomal translocations, and aberrant splicing at the posttranscriptional level [10]. For instance, FGFR1 mutations and amplifications have been implicated in prostate [11], breast [12] and lung cancers [13], whereas aberrant FGFR2 expression has been observed in endometrial [14], breast [15], and gastric cancer [16]. Furthermore, activating mutations and overexpression of FGFR3 have been reported in bladder cancers [17] and multiple myeloma [18]. We have identified mutation in FGFR4, which acts as an oncogene, and suggested therapeutic targeting of FGFR4 in rhabdomyosarcomas [19]. In addition, we recently examined a prevalent occupational exposure susceptibility variants associated with increased bladder cancer risk, and observed an additive interaction for rs798766 (TMEM129-TACC3-FGFR3) with the interaction more apparent in patients with tumors positive for FGFR3 expression [20]. In cervical cancer, recent studies have shown the possible involvement of aberrant FGFR signaling with human papillomavirus (HPV) 16 E5 expression [21]. Notably, FGFR2 expression has been reported to be associated with cell growth and progression of cervical dysplasia [22, 23].

Recently, alterations in the FGFR signaling pathway have gained the spotlight due to their high incidence rate and therapeutic potential in malignancies [24]. Nevertheless, protein expression of the FGFR system and its prognostic significance in cervical cancer has not yet been elucidated. Thus, the aim of this study is to investigate the clinical significance of FGFR1, FGFR2, FGFR3 and FGFR4 expression in a well-defined cohort of cervical cancers, using immunohistochemistry and quantitative digital image analysis.


Patients and tumor samples

In the present study, we retrieved a total of 336 early stage cervical cancer patients treated in the Department of Gynecologic Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine between 2002 and 2009. Patients with rare histology such as sarcoma, malignant melanoma, and neuroendocrine carcinoma, and patients with limited availability of tissue block specimens were excluded from the tissue microarray (TMA) construction. 61 non-matched, non-adjacent normal epithelial tissues were used for the control. Tissue samples were collected from patients who had signed an informed consent form, which was approved by the Institutional Review Board at the Samsung Medical Center, Seoul, Korea (2009-09-002-002 and 2015-07-122). Some of the paraffin blocks were provided by the Korea Gynecologic Cancer Bank through Bio & Medical Technology Development Program of the Ministry of Education, Science and Technology, Korea (NRF-2012M3A9B8021800). This study was additionally approved by the Office of Human Subjects Research at the National Institute of Health.

All patients were treated primarily with radical hysterectomy with or without pelvic/para-aortic lymph node dissection. Patients received adjuvant radiotherapy with or without concurrent chemotherapy if they had risk factors. Following treatment, patients had follow-up examinations every 3 months for the first 2 years, every 6 months for the next 3 years, and once annually every year thereafter. Disease-free survival (DFS) was assessed from the date of surgery to the date of recurrence or the date of the last follow-up visit. Overall survival (OS) was measured from the date of surgery to the time of death, or for the living patients, to the date of last contact.

Tissue microarray construction and immunohistochemistry

Tissue microarrays (TMAs) were constructed from tissue blocks used for routine pathological evaluation. In each case, areas with the most representative histology were selected, and cylindrical tissue samples (0.6 mm) were cored from the donor block and extruded into the recipient array. We included three cores from separate tissue blocks of each patient to reduce sampling bias.

NCI60 cell microarray was constructed as previously described [25]. In brief, NCI60 cells were collected and then fixed in 70 % ethanol. After fixation, cells were pelleted and added 3 % low-melt agarose (SeaPlaque, Cambrex, Rockland, ME). The cells were immediately evenly dispersed in the warm by vortexing, and placed on ice for 2 min to form agar plugs. The agarose plug was transferred into a histology cassette (Tissue-TekTM, Sakura Finetek, Torrance, CA) and then processed in Tissue-Tek VIP processor (Sakura) and subsequently paraffin-embedded to form a donor block. Finally, the NCI60 cell microarray was constructed using 1.0 mm needle.

Immunohistochemical staining of FGFR1, FGFR2, FGFR3 and FGFR4 was performed on 4-μm sections of the TMA and NCI60 cell microarray, as described previously [26]. In brief, following deparaffinization and dehydration, heat-induced antigen retrieval was performed for 20 min in an antigen retrieval buffer of pH 9.0 (for FGFR3 and FGFR4) or pH 6.0 (for FGFR1 and FGFR2) (Dako, Carpinteria, CA) using a steam pressure cooker (Pascal, Dako). The endogenous peroxidase activity was blocked with 3 % H2O2 for 10 min. For FGFR1 and FGFR2, sections were incubated with protein block (Dako) for further 15 min. The sections were subsequently incubated with primary antibodies, a detailed list of which, along with dilutions, is provided in Additional file 1: Table S1. The antigen–antibody reaction was detected with Dako EnVision+Dual Link System-HRP (Dako) and DAB+ (3, 3′-Diaminobenzidine; Dako). In the negative control sections, the primary antibody was omitted. Positive controls included umbilical cord, stomach carcinoma, skin tissue, and prostate cancer for FGFR1, FGFR2, FGFR3, and FGFR4 antibodies, respectively. Tissue sections were lightly counterstained with hematoxylin and examined by light microscopy.

Quantitative evaluation of immunostaining

The evaluation of immunostaining was carried out using computer-assisted image analyzing software version (Visiopharm, Hoersholm, Denmark,). Immunohistochemically stained slides were scanned using a NanoZoomer 2.0 HT (Hamamatsu Photonics, Hamamatsu City, Japan) at 20× objective magnification (0.5-μm resolution), and captured digital images were then imported into the Visiopharm software. Each core was imported separately using the TMA workflow of the program.

Digital images consist of pixels, each of which is defined by a position and a value of intensity [27]. After training the system by digitally “painting” examples of the nucleus in the image, tumor nuclei were defined (Additional file 2: Figure S1A–B). Cytoplasm is further defined by outlining the defined nucleus (Additional file 2: Figure S1C). The mean intensity of DAB for each defined image is used for quantification of the expressions (Additional file 2: Figure S1D). The intensity of staining was categorized as 0, 1+, 2+, and 3+ according to the distribution pattern across cores. The final score was calculated by multiplying the intensity and percentage of staining resulting in a final histoscore of 0–300 [28].

Microarray gene expression profiling

To examine the prognostic significance of mRNA expression, microarray data were analyzed as described previously [29]. In brief, RNA was extracted from formalin-fixed, paraffin-embedded tissue sections with removal of nontumor elements. RNA was extracted using the High Pure RNA Paraffin kit (Roche Diagnostic, Mannheim, Germany) and subsequently, the whole genome cDNA mediated annealing selection and ligation (WG-DASL) assay was performed following the manufacturer’s instructions (Illumina, San Diego, CA). With the exclusion of inadequate samples, a total of 300 patient samples were evaluable and the expression data was deposited in the Gene Expression Omnibus (GEO) database ( GSE44001).

In-silico analysis for TCGA cervix

To evaluate the prognostic significance of FGFR mRNA expression, data from The Cancer Genome Atlas (TCGA) Research Network were also analyzed ( Pan-cancer normalized form of RNA-seq data of cervical cancers, which had been obtained using a Illumina HiSeq (Illumina, San Diego, CA, USA), were downloaded (version: 2015-02-24). For survival analysis, the mRNA expression value was dichotomized according to quartile values (lower than 25 percentile vs. higher than 75 percentile).

Statistical analysis

Statistical analysis was performed using the R software version 3.1.2. The expression level of the proteins according to the clinicopathological characteristics were analyzed using a Student’s t test. Analysis of the Spearman rho coefficient was used to assess associations and correlations between parameters. For survival analysis, expression values were dichotomized (positive vs. negative) with the cut-off values showing the most discriminative power in the univariate cox model for disease-free survival (Additional file 2: Figure S1E) (R package: survMisc). Survival distributions were estimated using the Kaplan–Meier method, and the relationship between survival and each parameter was analyzed with the log-rank test. A Cox proportional hazards model was created to identify independent predictors of survival. Statistical significance was considered to be present at values of P < 0.05.


Clinicopathological characteristics of patients

Clinicopathological characteristics of 336 patients are summarized in Table 1. Patients with IB2 or IIB who were treated primarily with radical surgery were included in the present study. One hundred and sixty patients (47.6 %) were treated with adjuvant radiation, with or without concurrent chemotherapy following radical surgery.
Table 1

Clinicopathological characteristics of the 336 cervical cancer patients





 Mean ± SD

48.9 ± 11.2


FIGO stage







Cell type







Tumor size

 ≤4 cm



 >4 cm



SCC Ag level (ng/ml)

 Median (range)

1.2 (0.1–65.1)









Depth of invasion

 <50 %



 >50 %



LN metastasis







PM involvement







Resection margin







Primary treatment

 OP only



 OP + RT









FIGO International Federation of Gynecology and Obstetrics; SCC squamous cell carcinoma; Ag antigen; AD adenocarcinoma; ASC adenosquamous cell carcinoma; LVSI lymphovascular space invasion; LN lymph node; PM parametrium; OP operation; RT radiotherapy; CCRT concurrent chemoradiotherapy

FGFRs expression in the NCI60 cell lines

The NCI60 cell lines are widely used cancer biology and drug discovery. In addition, the NCI60 cell microarray is a useful platform for antibody validation because immunohistochemistry (IHC) data of NCI60 cell microarray are used to predict antibody titer for IHC on multi-tumor TMA [25]. Thus, we performed IHC on NCI60 cell microarray to assist in assessment of antibody specificity.

A total of 58 cell lines were analyzed by hierarchical clustering with the continuous histoscore. As shown in Fig. 1, three categories were defined. Category 1 (n = 29) consists relatively of low FGFR1, FGFR2, and FGFR4 expression. In contrast, category 3 (n = 20) consists exclusively of high FGFR1, FGFR2, and FGFR4 expression. Furthermore, category 2 (n = 9) consists exclusively of low FGFR3 expression. FGFR2 and FGFR4 protein expression was correlated with mRNA expression (r = 0.308; P = 0.019 and r = 0.413; P = 0.004, respectively) (Fig. 1b). Although FGFR1 and FGFR3 did not meet statistical significance there is a positive correlation trend. These data suggest that FGFRs proteins were differentially expressed in the NCI60 cell lines and antibodies of FGFR isoforms had specificity.
Fig. 1
Fig. 1

FGFRs expression in the NCI60 cell lines. a Hierarchical clustering analysis for immunohistochemical expression of FGFR1, FGFR2, FGFR3, and FGFR4. Three categories were defined. Category 1 (n = 29) consists relatively of low FGFR1, FGFR2, and FGFR4 expression. In contrast, category 3 (n = 20) consists exclusively of high FGFR1, FGFR2, and FGFR4 expression. Category 2 (n = 9) consists exclusively of low FGFR3 expression. b Correlation between FGFRs mRNA and protein expression. mRNA expression level was measured by microarray gene expression profiling, whereas the protein expression was assessed by immunohistochemistry. FGFR2 and FGFR4 protein expression was correlated with mRNA expression (P = 0.019 and P = 0.004, respectively)

FGFR expression and its associations with clinicopathological features

FGFR1, FGFR2, FGFR3 and FGFR4 expression is mainly observed in the cytoplasm [3032], and representative examples of positive and negative staining are shown in Fig. 2. Determination of the cut-offs for expression of FGFRs was also referenced to overall survival. The P value of the predicted log-rank test of Kaplan–Meier results was plotted against the normalized staining values (Additional file 2: Figure S2) and supports the determination of cut-offs for FGFR2, FGFR3, and FGFR4 based on the histogram. FGFR1 expression lacked as clear a relationship between expression and outcome by this analysis. Among the 336 tumors investigated, the number of tumors exhibiting high FGFR expression was 88 (26.2 %, histoscore >122) for FGFR1, 167 (49.7 %, histoscore >58) for FGFR2, 211 (62.8 %, histoscore >57) for FGFR3 and 241 (71.7 %, histoscore >79) for FGFR4. When compared with normal tissues, the expression of FGFR1, FGFR2, and FGFR4 were higher and the expression of FGFR3 was lower in cancer tissues (Table 2). The expression of FGFR was cell type associated. FGFR1 was more highly expressed in adeno-/adenosquamous carcinoma, while FGFR2, FGFR3, and FGFR4 expression was more prominent in squamous cell carcinoma (P = 0.020, P < 0.001, P < 0.001, and P = 0.020, respectively) (Table 2). These results suggest that each FGFR potentially has a different role according to cell type in cervical cancers.
Fig. 2
Fig. 2

FGFR1, FGFR2, FGFR3, and FGFR4 expression in formalin-fixed, paraffin-embedded cervical cancer tissues. Representative immunohistochemical images of FGFR1 negative (a) and positive (b), FGFR2 negative (c) and positive (d), FGFR3 negative (e) and positive (f), FGFR4 negative (g) and positive (h) expression. The Scale bar represents 50 μm

Table 2

Correlation between FGFR expression and clinicopathological characteristics of cervical cancer






Mean histoscore [95 % CI]

Mean histoscore [95 % CI]

Mean histoscore [95 % CI]

Mean histoscore [95 % CI]

Normal vs. cancer


35 [27–42]

32 [23–40]

106 [89–124]

68 [57–80]


84 [77–91]

70 [64–77]

82 [76–87]

129 [122–136]

 P value







82 [75–90]

72 [64–79]

83 [76–89]

130 [121–138]


89 [70–107]

66 [48–83]

79 [65–93]

117 [97–136]

 P value





Cell type


78 [71–86]

77 [69–85]

89 [82–96]

133 [124–141]


99 [83–115]

51 [40–62]

59 [48–71]

113 [98–128]

 P value





Tumor size

 ≤4 cm

85 [77–93]

74 [667–83]

83 [76–90]

131 [122–139]

 >4 cm

79 [64–91]

59 [47–70]

79 [68–89]

120 [105–134]

 P value







84 [75–93]

73 [64–82]

84 [76–93]

137 [127–147]


82 [71–93]

67 [57–77]

78 [70–87]

114 [104–126]

 P value





Depth of invasion

 <50 %

88 [75–101]

73 [60–86]

85 [74–97]

132 [118–146]

 >50 %

81 [73–90]

70 [62–78]

80 [73–87]

126 [117–135]

 P value





LN metastasis


83 [75–91]

74 [66–82]

84 [77–91]

132 [124–141]


84 [70–98]

60 [48–72]

75 [65–83]

114 [98–129]

 P value





PM involvement


85 [78–93]

73 [66–80]

83 [77–90]

129 [121–137]


65 [48–82]

49 [34–64]

64 [50–78]

115 [91–139]

 P value





Resection margin


83 [76–90]

71 [64–78]

82 [76–88]

128 [121–136]


88 [50–127]

70 [42–99]

80 [52–108]

118 [75–161]

 P value





Primary treatment

 OP only

84 [74–94]

81 [70–91]

91 [81–100]

136 [125–147]

 OP + RT

81 [66–96]

67 [53–81]

73 [62–85]

127 [111–143]


83 [70–96]

54 [45–64]

71 [62–81]

115 [100–129]


86 [2–171]

96 [17–175]

96 [20–173]

119 [38–199]

 P value





FIGO International Federation of Gynecology and Obstetrics; SCC squamous cell carcinoma; Ag antigen; AD adenocarcinoma; ASC adenosquamous cell carcinoma; LVSI lymphovascular space invasion; LN lymph node; PM parametrium; OP operation; RT radiotherapy; CCRT concurrent chemoradiotherapy

In addition, the high expression of FGFR1, FGFR2, FGFR3 were negatively correlated with the parametrial involvement (P = 0.030, P = 0.005, and P = 0.010, respectively). Furthermore, FGFR2 expression was down-regulated in large-sized (P = 0.020) and lymph node metastatic (P = 0.048) tumors. FGFR4 expression was also down-regulated in lymphovascular space invasive (P = 0.004) and lymph node metastatic (P = 0.040) tumors. These results indicate that FGFR expression is associated with less aggressive phenotypes in cervical cancers.

In order to find the clustering of samples according to FGFR expression, a total of 336 cervical cancer cases were analyzed by hierarchical clustering with the continuous histoscore. As shown in Additional file 2: Figure S3, two categories were defined. Category 1 (n = 116) consists exclusively of high FGFR2 and FGFR3 expression. In contrast, category 2 (n = 220) consists exclusively of low FGFR1, FGFR2, and FGFR3 expression. Furthermore, there are significant differences for cell type, tumor size, lymph node metastasis and primary treatment between category 1 and category 2 (Additional file 1: Table S2). Notably, category 2 associated with advanced tumor phenotypes and had poor prognosis. Subsequently, we examined the correlation among FGFR1-4 expression. FGFR1 expression was positively correlated with FGFR2, FGFR3, or FGFR4 expression (r = 0.445, r = 0.366 and r = 0.373, respectively; all P < 0.001). Furthermore, FGFR2 expression was positively correlated with FGFR3, or FGFR4 expression (r = 0.576, and r = 0.413, respectively; all P < 0.001). There was also a significant correlation between FGFR3 and FGFR4 expression (R = 0.381, P < 0.001) (Fig. 3).
Fig. 3
Fig. 3

Relationship among FGFR1, FGFR2, FGFR3, and FGFR4 expression in cervical cancer tissues. There is a significant positive correlation between each FGFR expression (P < 0.001)

High expression of FGFR2, FGFR3, and FGFR4 predict longer survival

With a median follow-up period of 66 months (range 1–143), five-year disease-free survival and overall survival rates for the whole group were 87 % (95 % CI 83–91) and 96 % (95 % CI 93–98), respectively. A high expression of FGFR2, FGFR3, and FGFR4 was significantly associated with favorable disease-free survival (P = 0.006, P = 0.035, and P = 0.001, respectively) and overall survival (P = 0.003, P = 0.002, and P = 0.003, respectively) (Additional file 2: Figures S3A, S4A). The 5 year disease-free survival was 92, 90, and 91 % for patients with a high expression of FGFR2, FGFR3, and FGFR4, respectively, compared with 81, 82, and 76 % for patients with a low expression. Similarly, the 5 year overall survival was 99, 97, and 97 % for patients with positive FGFR2, FGFR3, and FGFR4 expression respectively, compared with 93, 93, and 91 % for patients with a negative expression (Fig. 4). Furthermore, patients with a combination of high expression of FGFR2 and FGFR3, FGFR2 and FGFR4, or FGFR3 and FGFR4 had significantly more favorable disease-free survival (P = 0.004, P < 0.001, P < 0.001, respectively) and overall survival (all P < 0.001) (Fig. 4; Additional file 2: Figures S3B, S4B). Additionally, the coexpression of all three FGFRs was significantly associated with favorable disease-free survival and overall survival (both P < 0.001) compared with the single expression of FGFR.
Fig. 4
Fig. 4

Kaplan–Meier survival curves for overall survival according to FGFR1, FGFR2, FGFR3, and FGFR4. a Cervical cancer patients with high FGFR2, FGFR3, and FGFR4 expression had longer overall survival (P = 0.003, P = 0.002, and P = 0.003, respectively) than those with low expression. b The combination of FGFR2, FGFR3, and FGFR4 was found to enhance prognostic accuracy for cervical cancer. The patients with FGFR2+/FGFR3+/FGFR4+ expression had significantly longer overall survival (P < 0.001) than those with FGFR2-/FGFR3-/FGFR4- expression. P values were obtained from log-rank tests

Using the Cox proportional hazards model, the expressions of FGFR2 and FGFR4 remained an independent prognostic factor for disease-free survival [hazard ratio = 0.51 (95 % CI 0.27–0.97), P = 0.040; hazard ratio = 0.51 (95 % CI 0.28–0.91), P = 0.020, respectively) (Table 3). With respect to overall survival, the expressions of FGFR2, FGFR3, and FGFR4 remained an independent factor [hazard ratio = 0.24 (95 % CI 0.07–0.83), P = 0.020; hazard ratio = 0.29 (95 % CI 0.11–0.78), P = 0.010; hazard ratio = 0.33 (95 % CI 0.13–0.83), P = 0.020, respectively). All the combination markers with FGFR2, FGFR3, and FGFR4 showed independent factors for disease-free survival and overall survival (Additional file 1: Table S3; Table 3). Notably, the prognostic value of FGFR expression was not influenced by cell type.
Table 3

Multivariate analysis of the association between prognostic variables and survival in cervical cancer patients


Disease-free survival

Overall survival

DFS HR [95 % CI], P value

OS HR [95 % CI], P value

FIGO stage (IB2/IIB vs. IB1/IIA)

1.72 [0.86–3.47], 0.090

2.00 [0.70–5.73], 0.190

Cell type (AD vs. SCC)

3.80 [2.10–6.86], <0.001

5.89 [2.37–14.65], <0.001

LN metastasis

4.15 [2.20–7.83], <0.001

2.89 [1.12–7.49], 0.030

Tumor size (>4 cm)

0.96 [0.49–1.88], 0.910

0.90 [0.32–2.54], 0.840

PM involvement

1.31 [0.58–2.94], 0.520

1.72 [0.52–5.69], 0.380


0.52 [0.25–1.11], 0.090

0.48 [0.16–1.50], 0.210


0.51 [0.27–0.97], 0.040

0.24 [0.07–0.83], 0.020


0.56 [0.30–1.03], 0.060

0.29 [0.11–0.78], 0.010


0.51 [0.28–0.91], 0.020

0.33 [0.13–0.83], 0.020


0.37 [0.17–0.82], 0.010

0.13 [0.03–0.60], 0.010


0.30 [0.13–0.68], 0.004

0.17 [0.05–0.63], 0.008


0.33 [0.15–0.69], 0.004

0.11 [0.03–0.44], 0.002


0.24 [1.10–0.60], 0.002

0.08 [0.02–0.40], 0.002

DFS disease-free survival; HR hazard ratio; OS overall survival; CI confidence interval; FIGO International Federation of Gynecology and Obstetrics; AD adenocarcinoma; SCC squamous cell carcinoma; LN lymph node; PM parametrium

Based on the clinical significance of FGFRs protein expressions, the clinical implications of FGFR1, FGFR2, FGFR3, and FGFR4 mRNA expressional levels were assessed. Notably, there was a similar finding between transcription and translation expression levels in cervical cancer. FGFR2 had a significantly favorable prognostic significance at the mRNA level from the GSE44001 dataset (Additional file 2: Figure S5). The RNA-seq dataset of TCGA cervix also showed a favorable prognostic significance of FGFR2 and FGFR3 (Additional file 2: Figure S6).


In the present study, we investigated the prognostic significance of FGFR1, FGFR2, FGFR3, and FGFR4 expression in a large cohort of cervical cancer patients. Here, we identified that the FGFRs were differentially expressed according to cell types of cervical cancer. Furthermore, we demonstrate for the first time that elevated expression of FGFR2, FGFR3, or FGFR4 predict favorable survival in cervical cancer patients. These results suggest that FGFR protein analysis can be a prerequisite in the diagnostic procedure of cervical cancer, and may guide the patients’ therapy. In addition, quantitative digital image analysis used in this study shows potential for more objective results with possible better prognostication. We tested methods of combining the FGFRs expression to predict outcome, as has been demonstrated in the combination of different targets within pathways [33, 34], however were unable to demonstrate a relationship of co-expression that is predictive of outcome.

Published prognostic significances of FGFRs are controversial as there are studies reporting higher recurrence rates with positive FGFR expression, and others finding significantly lower recurrence rates in FGFR expressing tumors [7, 10, 17, 3537]. FGFR2 has been found to be highly expressed in colorectal cancer and correlated with tumor growth, invasion, and angiogenesis, and stronger FGFR2 expression has been observed in the invasive front of colorectal cancer cells [35]. In addition, there was positive correlation between advanced tumor stage and a high expression of FGFR2 in rectal cancers [36]. On the obverse, down-regulation of FGFR2 in bladder cancers was associated with an adverse prognosis [7]. And, in FGFR2 amplified breast cancer cell lines, constitutive signaling appeared to confer a survival advantage over non amplified cell lines [38]. FGFR3 have also been discussed as a potential prognostic markers in several cancers, however data remain controversial [3941]. Mutations of the FGFR3 gene are one of the most frequent genetic alterations in bladder cancer and have been shown to be associated with tumors with a favorable prognosis [4244]. With respect to protein expression, correlation of FGFR3 overexpression with low-grade and low-stage bladder cancers has been reported [6, 7, 37]. Conversely, there are studies reporting an association between the overexpression of the FGFR3 protein and poor survival in breast cancer [45]. There are few studies examining the prognostic significance of FGFR4. In breast cancer, wild type FGFR4 has been proposed to be an important tumor suppressor via the regulation of genes controlling invasion like matrix metalloproteinase 1, suggesting loss of wild type FGFR4 would adversely influence disease progression [46]. Our previous report showed over expression of FGFR4 in rhabdomyosarcoma, and furthermore, activating mutation of FGFR4 can promote metastasis in human rhabdomyosarcomas [19].

These contradictory results may be attributed in part to the cell type specificity. For instance, IIIb and IIIc isoforms of FGFR1 and FGFR2 are expressed in epithelial and mesenchymal cells, respectively [47]. In addition, the cell type specificity can change when FGFRs are associated with cancer. Secondly, the prognostic significance can be different according to mechanism and nature of what is measured. Significance of protein expression can differ from other aberration such as mutation, amplification, or translocations frequently seen in the FGFR pathway. Our present results support the idea that FGFR protein expression is a strong indicator concerning prognosis and survival.

HPV is a well-known etiologic agent in cervical cancer, and persistent infection leads to a genomic instability and local immune suppression, which can lead to both the accumulation of genomic alterations in the host cell, as well as to the integration of the viral genome into the host genome [48]. Recently, genome-wide studies have described the genomic and epigenomic alterations of HPV-associated cancers, which highlighted multiple potential biomarkers and therapeutic targets [4953]. In addition to recurrent integrations in RAD51B, NR4A2, and TP63, additional genomic alterations were found in receptor tyrosine kinases, primarily FGFR2 and FGFR3, in HPV-positive head and neck squamous cell carcinoma [49]. FGFR2 and FGFR3 mutations have been identified among 17.6 % of HPV-positive tumors, and both mutations have been described in several cancer types, and are sensitive to FGFR inhibitors [54, 55]. Although there is the heterogeneity of HPV-related tumors at different anatomical sites, alteration of FGFR2 and FGFR3 could also be worthwhile to study in cervical cancers.

In the present study, we used quantitative digital image analysis for immunohistochemistry scoring. The advantage of this analysis is, once the analysis protocol has been defined for a given application, large volumes of image data can be processed with minimal user-interaction allowing a highly standardized output on a continuous scale. The prognostic significance of FGFR expression seen in this study, support the usefulness of the image analysis. As quantitative digital image analysis becomes more established, the acknowledgement of results on a continuous scale is likely to become widely accepted and beneficial for improved data mining.


In conclusion, the present study investigated the immunohistochemical expression of FGFR1, FGFR2, FGFR3, and FGFR4 in a large number of cervical cancer patients. FGFR2, FGFR3, and FGFR4 expression were elevated in squamous cell carcinoma, and a high expression of FGFR1, FGFR2, FGFR3, and FGFR4 were associated with a less aggressive phenotype. Furthermore, a high expression of FGFR2, FGFR3, and FGFR4 showed longer disease-free survival and overall survival. Cox regression analysis confirmed that FGFR2, FGFR3, and FGFR4 expression was an important prognostic indicator in cervical cancer. This information could have clinical value in identifying cervical cancer patients who are at low risk of progression.




fibroblast growth factor receptor


International Federation of Gynecology and Obstetrics


fibroblast growth factor


tissue microarray


gene expression omnibus


The Cancer Genome Atlas


squamous cell carcinoma






adenosquamous cell carcinoma


lymphovascular space invasion


lymph node








concurrent chemoradiotherapy


confidence interval


Authors’ contributions

CHC, J-YC, B-GK, and SMH conceived the study and devised the experimental design. B-GK and SMH designed and build the tissue microarrays. CHC and J-YC performed experiments. CHC, J-YC, B-GK and SMH performed data analysis for experiments or clinical records. CHC and J-YC drafted the final version of the manuscript and figure legends. B-GK and SMH revised the figures, added critical content to the discussion and were responsible for revising all portions of the submitted manuscript. All authors read and approved the final manuscript.


This study was supported in part by a Grant from the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Republic of Korea (2013R1A1A2013629) and the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.

Competing interests

The authors declare that they have no competing interests.

Open AccessThis 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.

Authors’ Affiliations

Experimental Pathology Laboratory, Laboratory of Pathology, National Cancer Institute, Center for Cancer Research, National Institutes of Health, MSC 1500, Bethesda, MD 20892, USA
Department of Obstetrics and Gynecology, Gangnam Severance Hospital, Yonsei University College of Medicine, Seoul, 135-720, Korea
Department of Obstetrics and Gynecology, Samsung Medical Center, Sungkyunkwan University School of Medicine, 50 Irwon-dong, Gangnam-gu, Seoul, 135-710, Republic of Korea


  1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.View ArticlePubMedGoogle Scholar
  2. Suh DH, Kim JW, Aziz MF, Devi UK, Ngan HY, Nam JH, et al. Asian society of gynecologic oncology workshop 2010. J Gynecol Oncol. 2010;21:137–50.View ArticlePubMedPubMed CentralGoogle Scholar
  3. Thomas GM. Improved treatment for cervical cancer–concurrent chemotherapy and radiotherapy. N Engl J Med. 1999;340:1198–200.View ArticlePubMedGoogle Scholar
  4. Gien LT, Beauchemin MC, Thomas G. Adenocarcinoma: a unique cervical cancer. Gynecol Oncol. 2010;116:140–6.View ArticlePubMedGoogle Scholar
  5. Lee YY, Choi CH, Kim TJ, Lee JW, Kim BG, Lee JH, et al. A comparison of pure adenocarcinoma and squamous cell carcinoma of the cervix after radical hysterectomy in stage IB-IIA. Gynecol Oncol. 2011;120:439–43.View ArticlePubMedGoogle Scholar
  6. Gomez-Roman JJ, Saenz P, Molina M, Cuevas Gonzalez J, Escuredo K, Santa Cruz S, et al. Fibroblast growth factor receptor 3 is overexpressed in urinary tract carcinomas and modulates the neoplastic cell growth. Clin Cancer Res. 2005;11:459–65.PubMedGoogle Scholar
  7. Tomlinson DC, Lamont FR, Shnyder SD, Knowles MA. Fibroblast growth factor receptor 1 promotes proliferation and survival via activation of the mitogen-activated protein kinase pathway in bladder cancer. Cancer Res. 2009;69:4613–20.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Dubrulle J, Pourquie O. fgf8 mRNA decay establishes a gradient that couples axial elongation to patterning in the vertebrate embryo. Nature. 2004;427:419–22.View ArticlePubMedGoogle Scholar
  9. Feldman B, Poueymirou W, Papaioannou VE, DeChiara TM, Goldfarb M. Requirement of FGF-4 for postimplantation mouse development. Science. 1995;267:246–9.View ArticlePubMedGoogle Scholar
  10. Turner N, Grose R. Fibroblast growth factor signalling: from development to cancer. Nat Rev Cancer. 2010;10:116–29.View ArticlePubMedGoogle Scholar
  11. Giri D, Ropiquet F, Ittmann M. Alterations in expression of basic fibroblast growth factor (FGF) 2 and its receptor FGFR-1 in human prostate cancer. Clin Cancer Res. 1999;5:1063–71.PubMedGoogle Scholar
  12. Theillet C, Adelaide J, Louason G, Bonnet-Dorion F, Jacquemier J, Adnane J, et al. FGFRI and PLAT genes and DNA amplification at 8p12 in breast and ovarian cancers. Genes Chromosomes Cancer. 1993;7:219–26.View ArticlePubMedGoogle Scholar
  13. Weiss J, Sos ML, Seidel D, Peifer M, Zander T, Heuckmann JM, et al. Frequent and focal FGFR1 amplification associates with therapeutically tractable FGFR1 dependency in squamous cell lung cancer. Sci Transl Med. 2010;2:62ra93.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Byron SA, Gartside M, Powell MA, Wellens CL, Gao F, Mutch DG, et al. FGFR2 point mutations in 466 endometrioid endometrial tumors: relationship with MSI, KRAS, PIK3CA, CTNNB1 mutations and clinicopathological features. PLoS One. 2012;7:e30801.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Antoniou AC, Spurdle AB, Sinilnikova OM, Healey S, Pooley KA, Schmutzler RK, et al. Common breast cancer-predisposition alleles are associated with breast cancer risk in BRCA1 and BRCA2 mutation carriers. Am J Hum Genet. 2008;82:937–48.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Jang JH, Shin KH, Park JG. Mutations in fibroblast growth factor receptor 2 and fibroblast growth factor receptor 3 genes associated with human gastric and colorectal cancers. Cancer Res. 2001;61:3541–3.PubMedGoogle Scholar
  17. van Rhijn BW, Montironi R, Zwarthoff EC, Jobsis AC, van der Kwast TH. Frequent FGFR3 mutations in urothelial papilloma. J Pathol. 2002;198:245–51.View ArticlePubMedGoogle Scholar
  18. Chesi M, Nardini E, Brents LA, Schrock E, Ried T, Kuehl WM, et al. Frequent translocation t(4;14)(p16.3;q32.3) in multiple myeloma is associated with increased expression and activating mutations of fibroblast growth factor receptor 3. Nat Genet. 1997;16:260–4.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Taylor JGt, Cheuk AT, Tsang PS, Chung JY, Song YK, Desai K, et al. Identification of FGFR4-activating mutations in human rhabdomyosarcomas that promote metastasis in xenotransplanted models. J Clin Invest. 2009;119:3395–407.PubMedGoogle Scholar
  20. Figueroa JD, Koutros S, Colt JS, Kogevinas M, Garcia-Closas M, Real FX, et al. Modification of occupational exposures on bladder cancer risk by common genetic polymorphisms. J Natl Cancer Inst. 2015;107:djv223.View ArticlePubMedGoogle Scholar
  21. Ranieri D, Belleudi F, Magenta A, Torrisi MR. HPV16 E5 expression induces switching from FGFR2b to FGFR2c and epithelial-mesenchymal transition. Int J Cancer. 2015;137(1):61–72.View ArticlePubMedGoogle Scholar
  22. Kawase R, Ishiwata T, Matsuda Y, Onda M, Kudo M, Takeshita T, et al. Expression of fibroblast growth factor receptor 2 IIIc in human uterine cervical intraepithelial neoplasia and cervical cancer. Int J Oncol. 2010;36:331–40.PubMedGoogle Scholar
  23. Kurban G, Ishiwata T, Kudo M, Yokoyama M, Sugisaki Y, Naito Z. Expression of keratinocyte growth factor receptor (KGFR/FGFR2 IIIb) in human uterine cervical cancer. Oncol Rep. 2004;11:987–91.PubMedGoogle Scholar
  24. Brooks AN, Kilgour E, Smith PD. Molecular pathways: fibroblast growth factor signaling: a new therapeutic opportunity in cancer. Clin Cancer Res. 2012;18:1855–62.View ArticlePubMedGoogle Scholar
  25. Williams R, Chung JY, Ylaya K, Whiteley G, Hewitt SM. Characterizations and validations of novel antibodies toward translational research. Proteomics Clin Appl. 2010;4:618–25.View ArticlePubMedGoogle Scholar
  26. Choi CH, Chung JY, Park HS, Jun M, Lee YY, Kim BG, et al. Pancreatic adenocarcinoma up-regulated factor expression is associated with disease-specific survival in cervical cancer patients. Hum Pathol. 2015;46:884–93.View ArticlePubMedGoogle Scholar
  27. Hamilton PW, Bankhead P, Wang Y, Hutchinson R, Kieran D, McArt DG, et al. Digital pathology and image analysis in tissue biomarker research. Methods. 2014;70:59–73.View ArticlePubMedGoogle Scholar
  28. Ishibashi H, Suzuki T, Suzuki S, Moriya T, Kaneko C, Takizawa T, et al. Sex steroid hormone receptors in human thymoma. J Clin Endocrinol Metab. 2003;88:2309–17.View ArticlePubMedGoogle Scholar
  29. Lee YY, Kim TJ, Kim JY, Choi CH, Do IG, Song SY, et al. Genetic profiling to predict recurrence of early cervical cancer. Gynecol Oncol. 2013;131:650–4.View ArticlePubMedGoogle Scholar
  30. Fukui S, Nawashiro H, Otani N, Ooigawa H, Nomura N, Yano A, et al. Nuclear accumulation of basic fibroblast growth factor in human astrocytic tumors. Cancer. 2003;97:3061–7.View ArticlePubMedGoogle Scholar
  31. So F, Daley TD, Jackson L, Wysocki GP. Immunohistochemical localization of fibroblast growth factors FGF-1 and FGF-2, and receptors FGFR2 and FGFR3 in the epithelium of human odontogenic cysts and tumors. J Oral Pathol Med. 2001;30:428–33.View ArticlePubMedGoogle Scholar
  32. Qian ZR, Sano T, Asa SL, Yamada S, Horiguchi H, Tashiro T, et al. Cytoplasmic expression of fibroblast growth factor receptor-4 in human pituitary adenomas: relation to tumor type, size, proliferation, and invasiveness. J Clin Endocrinol Metab. 2004;89:1904–11.View ArticlePubMedGoogle Scholar
  33. Kitano H, Chung JY, Ylaya K, Conway C, Takikita M, Fukuoka J, et al. Profiling of phospho-AKT, phospho-mTOR, phospho-MAPK and EGFR in non-small cell lung cancer. J Histochem Cytochem. 2014;62:335–46.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Chung JY, Hong SM, Choi BY, Cho H, Yu E, Hewitt SM. The expression of phospho-AKT, phospho-mTOR, and PTEN in extrahepatic cholangiocarcinoma. Clin Cancer Res. 2009;15:660–7.View ArticlePubMedGoogle Scholar
  35. Matsuda Y, Ishiwata T, Yamahatsu K, Kawahara K, Hagio M, Peng WX, et al. Overexpressed fibroblast growth factor receptor 2 in the invasive front of colorectal cancer: a potential therapeutic target in colorectal cancer. Cancer Lett. 2011;309:209–19.View ArticlePubMedGoogle Scholar
  36. Li CF, He HL, Wang JY, Huang HY, Wu TF, Hsing CH, et al. Fibroblast growth factor receptor 2 overexpression is predictive of poor prognosis in rectal cancer patients receiving neoadjuvant chemoradiotherapy. J Clin Pathol. 2014;67:1056–61.View ArticlePubMedGoogle Scholar
  37. Lindgren D, Liedberg F, Andersson A, Chebil G, Gudjonsson S, Borg A, et al. Molecular characterization of early-stage bladder carcinomas by expression profiles, FGFR3 mutation status, and loss of 9q. Oncogene. 2006;25:2685–96.View ArticlePubMedGoogle Scholar
  38. Turner N, Lambros MB, Horlings HM, Pearson A, Sharpe R, Natrajan R, et al. Integrative molecular profiling of triple negative breast cancers identifies amplicon drivers and potential therapeutic targets. Oncogene. 2010;29:2013–23.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Hernandez S, Lopez-Knowles E, Lloreta J, Kogevinas M, Amoros A, Tardon A, et al. Prospective study of FGFR3 mutations as a prognostic factor in nonmuscle invasive urothelial bladder carcinomas. J Clin Oncol. 2006;24:3664–71.View ArticlePubMedGoogle Scholar
  40. Lamy A, Gobet F, Laurent M, Blanchard F, Varin C, Moulin C, et al. Molecular profiling of bladder tumors based on the detection of FGFR3 and TP53 mutations. J Urol. 2006;176:2686–9.View ArticlePubMedGoogle Scholar
  41. Zieger K, Dyrskjot L, Wiuf C, Jensen JL, Andersen CL, Jensen KM, et al. Role of activating fibroblast growth factor receptor 3 mutations in the development of bladder tumors. Clin Cancer Res. 2005;11:7709–19.View ArticlePubMedGoogle Scholar
  42. Billerey C, Chopin D, Aubriot-Lorton MH, Ricol D, Diez Gil, de Medina S, Van Rhijn B, et al. Frequent FGFR3 mutations in papillary non-invasive bladder (pTa) tumors. Am J Pathol. 2001;158:1955–9.View ArticlePubMedPubMed CentralGoogle Scholar
  43. van Oers JM, Wild PJ, Burger M, Denzinger S, Stoehr R, Rosskopf E, et al. FGFR3 mutations and a normal CK20 staining pattern define low-grade noninvasive urothelial bladder tumours. Eur Urol. 2007;52:760–8.View ArticlePubMedGoogle Scholar
  44. van Rhijn BW, Lurkin I, Radvanyi F, Kirkels WJ, van der Kwast TH, Zwarthoff EC. The fibroblast growth factor receptor 3 (FGFR3) mutation is a strong indicator of superficial bladder cancer with low recurrence rate. Cancer Res. 2001;61:1265–8.PubMedGoogle Scholar
  45. Kuroso K, Imai Y, Kobayashi M, Yanagimoto K, Suzuki T, Kojima M, et al. Immunohistochemical detection of fibroblast growth factor receptor 3 in human breast cancer: correlation with clinicopathological/molecular parameters and prognosis. Pathobiology. 2010;77:231–40.View ArticlePubMedGoogle Scholar
  46. Stadler CR, Knyazev P, Bange J, Ullrich A. FGFR4 GLY388 isotype suppresses motility of MDA-MB-231 breast cancer cells by EDG-2 gene repression. Cell Signal. 2006;18:783–94.View ArticlePubMedGoogle Scholar
  47. Yan G, Fukabori Y, McBride G, Nikolaropolous S, McKeehan WL. Exon switching and activation of stromal and embryonic fibroblast growth factor (FGF)-FGF receptor genes in prostate epithelial cells accompany stromal independence and malignancy. Mol Cell Biol. 1993;13:4513–22.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Rusan M, Li YY, Hammerman PS. Genomic landscape of human papillomavirus-associated cancers. Clin Cancer Res. 2015;21:2009–19.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature. 2015;517:576–82.View ArticleGoogle Scholar
  50. Seiwert TY, Zuo Z, Keck MK, Khattri A, Pedamallu CS, Stricker T, et al. Integrative and comparative genomic analysis of HPV-positive and HPV-negative head and neck squamous cell carcinomas. Clin Cancer Res. 2015;21:632–41.View ArticlePubMedPubMed CentralGoogle Scholar
  51. Ojesina AI, Lichtenstein L, Freeman SS, Pedamallu CS, Imaz-Rosshandler I, Pugh TJ, et al. Landscape of genomic alterations in cervical carcinomas. Nature. 2014;506:371–5.View ArticlePubMedPubMed CentralGoogle Scholar
  52. Parfenov M, Pedamallu CS, Gehlenborg N, Freeman SS, Danilova L, Bristow CA, et al. Characterization of HPV and host genome interactions in primary head and neck cancers. Proc Natl Acad Sci USA. 2014;111:15544–9.View ArticlePubMedPubMed CentralGoogle Scholar
  53. Akagi K, Li J, Broutian TR, Padilla-Nash H, Xiao W, Jiang B, et al. Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability. Genome Res. 2014;24:185–99.View ArticlePubMedPubMed CentralGoogle Scholar
  54. Cappellen D, De Oliveira C, Ricol D, de Medina S, Bourdin J, Sastre-Garau X, et al. Frequent activating mutations of FGFR3 in human bladder and cervix carcinomas. Nat Genet. 1999;23:18–20.PubMedGoogle Scholar
  55. Liao RG, Jung J, Tchaicha J, Wilkerson MD, Sivachenko A, Beauchamp EM, et al. Inhibitor-sensitive FGFR2 and FGFR3 mutations in lung squamous cell carcinoma. Cancer Res. 2013;73:5195–205.View ArticlePubMedPubMed CentralGoogle Scholar


© Choi et al. 2016