Prognostic biomarkers for oral tongue squamous cell carcinoma: a systematic review and meta-analysis.

BACKGROUND: Identifying informative prognostic biomarkers for oral tongue squamous cell carcinoma (OTSCC) is of great importance in order to better predict tumour behaviour and to guide treatment planning. Here, we summarise existing evidence regarding immunohistochemical prognostic biomarkers for OTSCC.
METHODS: A systematic search of the literature was performed using the databases of Scopus, Ovid Medline, Web of Science and Cochrane Library. All studies which had investigated the prognostic significance of immunohistochemical biomarkers in OTSCC during the period from 1985 to 2015 were retrieved. For the five most often evaluated biomarkers a random-effects meta-analysis on overall survival was performed, including those studies that provided the necessary statistical results.
RESULTS: A total of 174 studies conducted during the last three decades were found, and in these 184 biomarkers were evaluated for the prognostication of OTSCC. The five biomarkers most frequently assessed were p53, Ki-67, p16, VEGFs and cyclin D1. In the meta-analyses, the most promising results of the prognostic power for OTSCC were obtained for cyclin D1. For studies of VEGF A and C the results were equivocal, but the pooled analysis of VEGF A separately showed it to be a useful prognosticator for OTSCC. There was no sufficient evidence to support p53, Ki-67 and p16 as prognostic biomarkers for OTSCC. Limitations in the quality of the published studies (e.g., small cohorts, lack of compliance with REMARK guidelines) are widespread.
CONCLUSIONS: Numerous biomarkers have been presented as useful prognosticators for OTSCC, but the quality of the conduct and reporting of original studies is overall unsatisfactory which does not allow reliable conclusions. The value of two biomarkers (VEGF-A and cyclin D1) should be validated in a multicentre study setting following REMARK guidelines.

Oral tongue squamous cell carcinoma (OTSCC) is the most common malignancy of the oral cavity. OTSCC is increasing in incidence (Patel ; Ng ), and has an aggressive clinical behaviour with a relatively poor prognosis (Bello ). The 5-year relative survival rate was 63% in a recent report from the Netherlands (van Dijk ). During 2017, almost 16 400 new cases of tongue cancer, and 2400 deaths are projected to occur due to this cancer in the United States (Siegel ).

Predicting the outcome of OTSCC patients is important when planning treatment. In early stages (cT1–T2) of OTSCC, which are expected to have favourable prognosis, the cancer-related mortality affects about 19% of patients (Almangush ). This indicates the need for better prognostic tools. Identification of robust prognostic biomarkers that could accurately predict the behaviour of OTSCC will aid in the selection of appropriate treatment strategies.

Recent research in the field of molecular pathology has introduced thousands of tumour biomarkers, which are associated with the progression and/or prognosis of different cancers. Many of these biomarkers have been evaluated for their prognostic power in OTSCC. However, the use of any molecular biomarker for OTSCC in daily practice has not yet been approved, although several biomarkers have been presented as promising prognosticators that could provide added value upon the classical ones such as stage, tumour grade and depth of invasion.

The task of the current systematic review of literature was to retrieve original studies that have examined the prognostic value of immunohistochemical biomarkers of OTSCC and to meta-analyse the studies of the most repeatedly reported biomarkers. We also summarise the current understanding of the topic and highlight the main shortcomings of the published studies to improve future research in this field.

Materials and methods

Search strategy

A search strategy combining the terms (tongue) AND (cancer* OR squamous cell carcinoma* OR neoplas* OR tumo*) AND (prognos* OR predict* OR surviv* OR recur* OR mortal* OR metasta*) AND (immunohisto* OR protei* OR marke* OR biomark*) was developed. The search terms were entered into Scopus, Ovid Medline, Web of Science and Cochrane Library (1985–2015).

In advanced search, the following search fields were included: title, abstract, subject heading, and keyword. The Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) were utilised (Moher ).

Screening

Two independent researchers (AA & IH) examined the retrieved hits, and discarded duplicated ones. We also excluded unrelated studies through careful browsing of the title and/or abstract of each publication. Existing review articles for prognostic biomarkers of tongue cancer (Ferrari ; Bello , 2010b) were screened for papers missed in the search strategy.

Data extraction

For relevant articles, we retrieved information about the name of the first author, country, year of publication, number of patients and the immunohistochemical biomarker/s examined. For those biomarkers which were reported repeatedly, further data including the primary antibody used and its dilution, unadjusted and adjusted analyses, statistical results reported (estimated hazard ratio (HR), 95% confidence interval (CI) and P-value) were retrieved when available.

Exclusion criteria

Studies in languages other than English.

Data based on animal samples.

Studies on cancers other than SCC, or on rare histological variants of SCC.

Studies including samples from other subsites of the oral cavity, and studies that mixed samples from the oral tongue and the base of the tongue.

Studies that did not report the prognostic value of the biomarker (i.e., studies that only reported association between the biomarker and classical parameters, such as stage and grade, but did not provide results for the association between the biomarker and survival outcomes).

Quality assessment

The guidelines from Reporting Recommendations for Tumor Marker Prognostic Studies (REMARK) (McShane ; Altman ) were used to evaluate the quality of studies that were eligible for the meta-analyses of the five most often reported biomarkers. The selected guidelines taken from the REMARK criteria are summarised in Table 1.

Table 1

Evaluation criteria used to assess the quality of studies included in the meta-analysis of the five most often reported biomarkers (adapted from REMARK guidelines)

Checklist	Criteria
1- Samples	Cohort (Retrospective or Prospective) study with a well-defined study population
	Medical treatment applied to the patients was explained. Authors have explained if all patients have received the same treatment or not
2- Clinical data of the cohort	The basic clinical data such as age, gender, clinical stage and histopathologic grade was provided
3- Immunohistochemistry	Well-described staining protocol or referred to original paper
4- Prognostication	The analysed survival endpoints were defined (e.g., overall survival, disease-free survival)
5- Statistics	Cutoff point, which used to divide the cases into risk groups was well described
	Estimated effect (CI, HR) describing the relationship between the evaluated biomarker and the outcome was provided
	Adequate statistical analysis (e.g., Cox regression modelling) was performed to adjust the estimation of the effect of the biomarker for known prognostic factors
6- Classical prognostic factors	The prognostic value of the classical prognostic factors was reported
	The relationship between the evaluated biomarker/s and classical prognostic factors were reported

Statistical analysis

A meta-analysis on overall survival (OS) was performed for each of the five selected biomarkers, including only those original studies which provided an estimate of the HR and the associated 95% CI for the contrast of interest. As most studies provided only unadjusted or ‘univariate’ estimates of HR, the meta-analysis was primarily based on them. However, in case only an estimated HR was reported which was adjusted for various other prognostic factors using multiple Cox regression, such a ‘multivariate’ estimate was entered to the analysis. Depending on the biomarker the direction of HR contrast was either positive vs negative (high vs low; applied for p53, Ki-67, VEGF and cyclin D1) or negative vs positive (low vs high; applied for p16). Where an individual study reported the HR estimate in the opposite direction, the inverse HR and CI were calculated to obtain results concordant with other studies. HR and CIs were transformed into log(HR) and its standard error (s.e.). Pooled estimates of the HR were computed both by a fixed-effect and by a random-effects model based on the generic inverse variance approach (Schwartzer ). In the random-effects analysis the between-studies variance was obtained by the Sidik–Jonkmann method (Sidik and Jonkmann, 2005). The I2 statistic and tau-squared, the estimated heterogeneity variance, were used to measure the heterogeneity of HRs between studies. Analyses were performed and forest plots created using the functions ‘metagen’ and ‘forest’ contained in package ‘meta’ (Schwarzer, 2007) of the R environment (R Core Team, 2016).

Results

Search results

Our search retrieved 1817 hits, and 771 of these were relevant for our study (Figure 1). A total of 174 studies remained that evaluated a total of 184 biomarkers for OTSCC. Of these, 32 biomarkers (17.4%) were reported at least in three studies (Figure 2). Of the latter, five biomarkers including p53, Ki-67, p16, vascular endothelial growth factors (VEGFs) and cyclin D1 were reported more often than any other (Tables 2, 3, 4, 5, 6; Supplementary Tables 1–5).

Figure 1

PRISMA flowchart: studies included and excluded along the various steps.

Figure 2

Biomarkers evaluated in three studies or more for their prognostic value in OTSCC.

Table 2

Summary of studies assessing the prognostic value of p53 in OTSCC providing unadjusted or adjusted estimates of HR and their 95% CIs for one or more endpoints, the HRs contrasting positive to negative expression

(Authors, year) Country	Stage or tumour size (primary treatment)	Primary antibody	Cutoff point in percentage	No. of cases	No. of positive cases (%)	Endpoint	Unadjusted analysis	Adjusted analysis	Compliance to REMARK guidelines
(Hogmo et al, 1998) Sweden	I (S)	p53, polyclonal (1 : 700; Novocastra, Newcastle upon Tyne, UK)	Distinct brown staining, confined to the nuclei, at least in 1% of the cells	54	24 (48)	Local recurrence	NA	HR 1.05 (95% CI 0.29–3.81) P=0.936	Checklist no. 5 was not fulfilled
						Regional recurrence	P=0.03	HR 0.28 (95% CI 0.07–1.08), P=0.065
(Bova et al, 1999) Australia	T1–T4 (S only in 96 patients; and S with adjuvant RT in 52 patients)	p53, DO-7 (1 : 50; Dako, Carpinteria, CA)	10	143	89 (62.2)	DFS	HR 1.18 (95% CI 0.60–2.32), P=0.64	NA	Checklist no. 5 was not fulfilled
						OS	HR 1.46 (95% CI 0.69–3.06), P=0.32	NA
(Mineta et al, 2002) Japan	I–IV (S)	p53, DO-7 (1 : 1000; Dako, Copenhagen, Denmark)	50	109	56 (52.8)	RFS	NA	HR 1.31 (95% CI 0.77–2.23), P=0.328	Checklist no. 5 and 6 were not fulfilled
(Goto et al, 2005) Japan	I–IV (S)	p53, DO-7 (1 : 100; Dako, Glostrup, Denmark)	10	123	62 (50.4)	DFS	HRa 1.02 (95% CI 0.59–1.75)	NA	Checklist no. 5 was not fulfilled
						DSS	HRa 0.91 (95% CI 0.37–2.24)	NA
						OS	HRa 1.27 (95% CI 0.59–2.71)	NA
(Trivedi et al, 2011) India	T1–T4 (S followed by RT in 26 cases, CT in 1 case, RT with CT in 5 cases)	p53, DO-7 (1 : 50; Dako, Glostrup, Denmark)	10	43	21 (49)	RFS	HR 3.66 (95% CI 1.42–9.39), P=0.007	NA	Checklist no. 5 and 6 were not fulfilled
						OS	HR 4.25 (95% CI 1.39–12.98), P=0.001	HR 3.69 (95% CI 1.01–13.43), P=0.048
(Adduri et al, 2014) India	T1–T3 (S)	p53, DO-1 (1 μg ml−1, EMD Millipore Calbiochem, Darmstadt, Germany)	20	79	49 (62)	DSS	HR 3.35 (95% CI 1.83–6.14), P=0.0003	NA	Checklist no. 5 and 6 were not fulfilled

Abbreviations: CI=confidence interval; CT=chemotherapy; DFS=disease-free survival; DSS=disease-specific survival; HR=hazard ratio; NA=not available; OS=overall survival; RFS=relapse-free survival; RT=radiotherapy; S=surgery.

HR in Goto et al study has been calculated by contrasting negative p53 vs positive p53, and we have converted it (as 1 divided by the reported HR and 1 divided by the reported CI) to be similar to all other studies in which HR was calculated by contrasting positive p53 vs negative p53.

Table 3

Summary of studies assessing the prognostic value of Ki-67 in OTSCC providing unadjusted or adjusted estimates of HR and their 95% CIs for one or more endpoints, the HRs contrasting positive to negative expression

(Authors, year) Country	Stage or tumour size (primary treatment)	Primary antibody	Cutoff point in percentage	No. of cases	No. of positive cases (%)	Endpoint	Unadjusted analysis	Adjusted analysis	Compliance to REMARK guidelines
(Bova et al, 1999) Australia	I–IV (S in 96 patients, S+RT in 52 patients)	Ki-67, clone 7B11 (1 : 50; Zymed Laboratories, San Francisco, CA)	50	148	NA	DFS	HR 1.17 (95% CI 0.61–2.27), P=0.63	NA	Checklist no. 5 was not fulfilled
						OS	HR 1.24 (95% CI 0.62–2.47), P=0.54	NA
(Mineta et al, 2002) Japan	I–IV (Radical treatment)	MIB-1, monoclonal (1 : 50; Immunotech, Marseille, France)	10	109	52 (48)	RFS	NA	HRa 0.82 (95% CI 0.47–1.44)	Checklist no. 1 was not fulfilled
(Sakata et al, 2008) Japan	I–IV (LDR brachytherapy alone in all patients+neck dissection for 9 who had N+)	MIB-1, monoclonal (Dako, Copenhagen, Denmark)	38	68	NA	Local RFS	P=0.59	HR 0.59 (95%CI 0.0514–6.74), P=0.7	Checklist no. 5 was not fulfilled
(Wang et al, 2009) China/USA	I–IV (S)	Ki-67, (1 : 200; Abcam)	NA	42	NA	OS	HR 0.55 (95%CI 0.03–8.76), P=0.67	NA	Checklist no. 5 was not fulfilled
(Bello et al, 2011) Finland/Israel	I–IV (S in 68 patients, S+RT in 49 patients, S+CRT in 10 patients+S+CT in 1 patient)	Ki-67 polyclonal (1 : 50, Dako A/S, Denmark)	40.4	77	18 (23.3)	DSS	HR 0.89 (95% CI 0.31–2.58)	NA	Checklist no. 5 was not fulfilled
(Bockelman et al, 2011) Finland	T1–T2 (S in 40 patients, S+RT in 33 patients)	Ki-67, monoclonal (1 : 100, Dako, Cytomation, Denmark)	30	73	53 (80.3)	OS	HRb 1.31 (95% CI 0.58–2.98)	NA	Checklist no. 5 was not fulfilled
(Hwa et al, 2015) South Korea	T1–T2 (S in 17 patients, S+RT in 8 patients)	Ki-67, monoclonal (1 : 1000; Dako, Carpinteria, CA, USA)	50	25	8 (32)	DSS	HR 0.51 (95% CI 0.05–5.05), P=0.57	NA	Checklist no. 5 was not fulfilled

Abbreviations: CI=confidence interval; CT=chemotherapy; DFS=disease-free survival; DSS=disease-specific survival; HR=hazard ratio; LDR=low-dose-rate; NA=not available; OS=overall survival; RFS=relapse-free survival; RT=radiotherapy; S=surgery.

HR in Mineta et al study has been calculated by contrasting low Ki-67 vs high Ki-67, and we have converted it (as 1 divided by the reported HR and 1 divided by the reported CI) to be similar to all other studies in which HR was calculated by contrasting high Ki-67 vs low Ki-67.

HR of 1.31 (95% CI 0.58–2.98) was kindly calculated by the author for this meta-analysis.

Table 4

Summary of studies assessing the prognostic value of p16 in OTSCC providing unadjusted or adjusted estimates of HR and their 95% CIs for one or more endpoints, the HRs contrasting negative to positive expression

(Authors, year) Country	Stage or tumour size (primary treatment)	Primary antibody	Cutoff point in percentage	No. of cases	No. of positive cases (%)	Endpoint	Unadjusted analysis	Adjusted analysis	Compliance to REMARK guidelines
(Bova et al, 1999) Australia	T1–T4 (S only in 96 patients. S with adjuvant RT in 52 patients)	p16INK4A, clone ZJ11 (1 : 60; Neomarkers, Fremont, CA)	1	147	65 (45.5)	OS	HR 2.50 (95% CI 1.16–5.41), P=0.02	HR 2.32 (95% CI 1.04–5.16), P=0.04	Fulfilled items summarised in Table 1
						DFS	HR 2.73 (95% CI 1.32–5.67), P=0.007	HR 3.23 (95% CI 1.49–6.99), P=0.003
(Goto et al, 2005) Japan	I–IV (S)	p16INK4A, clone E6H4 (dilution 1 : 25, Dako, Glostrup, Denmark)	10	123	42 (34.1)	OS	HR 1.41 (95% CI 0.48–4.093), P=0.52	NA	Checklist no. 5 was not fulfilled
						DSS	HR 1.30 (95% CI 0.37–4.47), P=0.67	NA
						DFS	HR 1.84 (95%CI 0.78–4.33), P=0.15	NA
(Trivedi et al, 2011) India	T1–T4 (S followed by RT in 26 cases, CT in 1 case, RT with CT in 5 cases)	p16, Monoclonal, clone SC-166 (1 : 150; Santa Cruz Biotechnology, USA)	10	61	28 (46)	OS	HR 1.59 (95% CI 0.59–4.24), P=0.354	NA	Checklist no. 1 and 5 were not fulfilled
						RFS	HR 1.47 (95% CI 0.62–3.49), P=0.372	NA
(Harris et al, 2011) USA	40% were T1 (S in 23 cases, Radio-chemotherapy in 2 cases)	p16INK4A, monoclonal, MAB4133 (Chemicon International Company, Temecula, CA)	High= diffuse and intense staining. Low=no, weak and moderate staining	25	11 (44)	OS	HRa 3.57 (95% CI 1–12.5)	NA	Checklist no. 1 and 5 were not fulfilled
						RFS	HRa 4.35, (95% CI 1.35–14.3)	NA
(Ramshankar et al, 2014) India	I–II (S or brachytherapy)	p16, monoclonal, clone—JC8-sc-56330 (1 : 150; Santa Cruz Biotechnology Santa Cruz, CA, USA)	50	156	24 (15.4)	OS	HRa 0.32 (95% CI 0.11–0.98)	HRa 0.42 (95% CI 0.23–0.77)	Fulfilled items summarised in Table 1
						DFS	HRa 0.212 (95% CI 0.07–0.66)	HRa 0.387 (95% CI 0.22–0.69)
(Lim et al, 2014) Australia	I–IV (S, S with adjuvant RT, S with adjuvant CT, CT or RT alone)	p16, monoclonal, E6H4; (Cintec/Roche)	70	116	6 (5)	OS	HR 1.81 (95% CI: 0.44–7.47), P=0.40	NA	Checklist 5 was not fulfilled
						PFS	HR 1.14 (95% CI: 0.45–2.87), P=0.78	NA

Abbreviations: CI=confidence interval; CT=chemotherapy; DFS=disease-free survival; DSS=disease-specific survival; HR=hazard ratio; NA=not available; OS=overall survival; PFS= progression free survival; RFS=relapse-free survival; RT=radiotherapy; S=surgery.

HR in Harris et al study and Ramshankar et al study has been calculated by contrasting positive p16 vs negative p16, and we have converted it (as 1 divided by the reported HR and 1 divided by the reported CI) to be similar to all other studies in which HR was calculated by contrasting negative p16 vs positive p16.

Table 5

Summary of studies assessing the prognostic value of vascular endothelial growth factors (VEGFs) in OTSCC providing unadjusted or adjusted estimates of HRand their 95% CIs for one or more endpoints, the HRs contrasting positive to negative expression

(Authors, year) Country	Stage or tumour size (primary treatment)	Primary antibody	Cutoff point in percentage	No. of cases	No. of positive cases (%)	Endpoint	Unadjusted analysis	Adjusted analysis	Compliance to REMARK guidelines
(Mineta et al, 2002) Japan	I-IV (Radical treatment)	VEGF-A, monoclonal, M293 (1 : 1000; R&D, Abingdon, UK)	10	109	69 (63.3)	OS	P=0.0002	NA	Checklist no. 5 was not fulfilled
						RFS	P<0.0001	HRa 8.69 (95% CI 3.92–19.23)
(Shao et al, 2008) China	T1–T4 (S)	VEGF-A, monoclonal (1 : 100; ZYMED Ltd., South San Francisco, CA)	50	59	44 (74.6)	OS	HR 11.14 (95% CI 1.99–62.35), P=0.004	HR 8.74 (95% CI 2.13–53.38), P=0.045	Checklist no. 5 was not fulfilled
						Recurrence	P=0.019	NA
(Petera et al, 2015) Czech	T1–T3 (S+HDR BT in all patients, 12 had also ND)	VEGF-A, polyclonal (1 : 200; Abbiotec, San Diego, CA)	NA	30	NA	OS	HR 5.40 (95% CI 1.25–23.25), P=0.023	NA	Checklist no. 5 was not fulfilled
						DFS	HR 3.82 (95% CI 1.01–14.49), P=0.049	P=0.04
						Local recurrence	NA	NA
						Regional recurrence	HR 5.52 (95% CI 1.06–28.78), P=0.042	NA
(Imayama et al, 2015) Japan	I–IV (S)	VEGF-A (Santa Cruz Biotechnology, Inc., Texas, USA)	HSCORE ≥75 (percentage stained multiplied by intensity score)	61	30 (49.2)	DSS	VEGF-A: HR 8.24 (95% CI 1.01–67.07), P<0.05	VEGF-A: HR 3.68 (95% CI 0.19–73.31), P=0.37	All checklist items were fulfilled
		VEGF-C (Santa Cruz Biotechnology, Inc., Texas, USA)			23 (37.7)	DSS	VEGF-C: HR 3.155 (95% CI 0.75–13.27), P=0.12	NA	Checklist no. 5 was not fulfilled
(Morita et al, 2014) Japan	I–IV (S+ND)	VEGF-C, polyclonal (1 : 100; Novus Biologicals Inc, Littleton, CO, USA)	5	40	18 (45)	OS	P=0.27	HR 0.29 (95% CI 0.04–1.35), P=0.12	All checklist items were fulfilled

Abbreviations: CI=confidence interval; DFS=disease-free survival; DSS=disease-specific survival; HR=hazard ratio; HDR BT=high dose rate brachytherapy; NA=not available; ND=neck dissection; OS=overall survival; RFS=relapse-free survival; S=surgery.

HR in Mineta et al study has been calculated by contrasting VEGF negative vs VEGF positive, and we have converted it (as 1 divided by the reported HR and 1 divided by the reported CI) to be similar to all other studies in which HR was calculated by contrasting VEGF positive vs VEGF negative.

Table 6

Summary of studies assessing the prognostic value of cyclin D1 in OTSCC providing unadjusted or adjusted estimates of HR and their 95% confidence intervals for one or more endpoints, the HRs contrasting positive to negative expression

(Authors, Year) Country	Stage or tumour size (primary treatment)	Primary antibody	Cutoff point in percentage	No. of cases	No. of Positive cases (%)	Endpoint	Unadjusted analysis	Adjusted analysis	Compliance to REMARK guidelines
(Bova et al, 1999) Australia	I–IV (S in 96 patients; S with adjuvant RT in 52 patients)	NCL-cyclin D1-GM, monoclonal (1 : 50; Novocastra, Rockdale, Australia)	10	147	100 (68)	OS	HR 3.89 (95% CI 1.37–11.07), P=0.01	HR 4.2 (95% CI 1.23–14.09), P=0.02	Fulfilled items summarised in Table 1
						DFS	HR 2.50 (95% CI 1.32–5.15), P=0.006	HR 2.48 (95% CI 1.0–6.15), P=0.05
(Trivedi et al, 2011) India	T1–T4 (S followed by RT in 26 cases, CT in 1 case, RT with CT in 5 cases)	P2D11F11, monoclonal, (dilution 1 : 50; Novocastra laboratories, Newcastle, UK)	10	61	25 (41)	OS	HR 2.20 (95% CI 0.85–5.70) P=0.103	NA	Checklist no. 5 was not fulfilled
						RFS	HR 1.91 (95% CI 0.84–4.38), P=0.122	NA

Abbreviations: CI=confidence interval; CT=chemotherapy; DFS=disease-free survival; HR=hazard ratio; NA=not available; OS=overall survival; RT=radiotherapy; RFS=relapse-free survival; S=surgery.

A clear majority (86%) of the studies reported at least one biomarker as a promising prognosticator for OTSCC based on obtaining a ‘statistically significant’ result for an association of that biomarker with at least one outcome variable considered. Most of the studies on these five biomarkers were based on quite small cohorts (<100 patients), and in the analyses of the outcome they commonly mixed early and late stage cancers (Tables 2, 3, 4, 5, 6; Supplementary Tables 1–5). It was also common to exclude a biomarker from an adjusted analysis using Cox regression, when the biomarker turned out to be statistically ‘non-significant’ in an unadjusted analysis, but biomarkers ‘significant’ in an unadjusted analysis were often included in an adjusted analysis. Many authors, though, presented the results based on unadjusted analysis only.

Most of the publications of those frequently studied biomarkers did not fulfil completely REMARK guidelines (Altman ). We used the selected criteria mentioned in Table 1 to evaluate the most often studied biomarkers, on which estimated HRs and their 95% CI were reported for at least one survival endpoint in a number of studies, and more than one study analysed the OS to be included in our meta-analyses (Tables 2, 3, 4, 5, 6). We particularly noted that guideline no. 5 (related to statistics) was not fulfilled in several studies (e.g., adjusted analysis using multiple Cox regression was not conducted), and guideline no. 1 (related to the patient series) was sometimes ignored by the authors (e.g., medical treatment applied to the patients was not well-explained). In addition, guideline no. 6 (related to classical prognostic factors) was not fulfilled properly in many studies (e.g., the relationship between biomarker/s and classical prognostic factors were not reported) (Tables 2, 3, 4, 5, 6).

Results of meta-analyses

The results of the meta-analyses on OS for the five most often studied biomarkers, each being evaluated in at least two studies that reported necessary statistical data, are summarised in Figures 3 and 4. The number of eligible studies was very small ranging from two for cyclin D1 to six for p16, the remaining biomarkers having three studies each. Essential heterogeneity in the HRs across the individual studies was observed for Ki-67, VEGFs and p16, whereas the pertinent measures (I2 and tau-squared) had very low values for p53 and cyclin D1. The point estimates and the error margins of the pooled HR from the fixed effect model and from the random effects model were very similar for the two last ones, whereas for the three first ones the 95% CI of the pooled HR was substantially wider when based on the random effects model as expected. Thus, it is reasonable to focus on the results of the random effects model. On the basis of the available data there was not sufficient evidence for p53, Ki-67 and p16 to be informative prognostic biomarkers. The two available studies on cyclin D1 suggest that this biomarker could be a useful prognosticator worth further evaluation. For cyclin D1 the pooled HR estimate was 2.86 (95% CI from 1.34 to 6.08). As regards VEGFs, the results were mixed with a very wide CI for HR from the random effects model. However, two of the VEGFs studies analysed VEGF-A and one VEGF-C. A positive or high VEGF-C value was shown to be associated with an improved prognosis (Morita ), which was in sharp contrast with the two other studies analysing the expression of VEGF-A, both indicating a much worse survival for a positive VEGF-A value. When the VEGF-C study was excluded from the meta-analysis, the pooled estimate for HR was 7.34 (95% CI from 2.32 to 23.22), which provides rather strong evidence for VEGF-A being a useful prognostic biomarker.

Figure 3

Forest plots for the pooled analyses of the biomarkers that have been studied most frequently in OTSCC but did not have prognostic usefulness in OTSCC. (A) p53 studies, (B) Ki-67 studies, and (C) p16 studies.

Figure 4

Forest plots for the pooled analyses of the biomarkers that have been studied most frequently in OTSCC and have shown prognostic usefulness in OTSCC. (A) VEGFs studies and (B) cyclin D1 studies.

Discussion

Molecular biomarkers may highlight biological differences between cancers and help to prognosticate patient outcome. During the last three decades (1985–2015), more than one hundred molecular biomarkers (identified by immunohistochemistry) were introduced as prognosticators for OTSCC. Five biomarkers including p53, Ki-67, p16, VEGFs and cyclin D1 were most often reported, and their biology has been reviewed elsewhere (Oliveira and Ribeiro-Silva, 2011; Wang ). According to the findings of the present meta-analysis, the prognostic usefulness of VEGF-A and cyclin D1 is worth further evaluation. On the other hand, for p53, Ki-67 and p16 there appears to be no sufficient evidence for any prognostic value for OTSCC.

The evaluation of molecular biomarkers in different subsites of the oral cavity is common in literature. However, variations in the immunohistochemical staining results reflect variations in proteomic (and genomic) properties of SCC between different oral subsites. For example, various immunohistological biomarkers analysed in OTSCC and buccal carcinoma samples did not associate with survival in OTSCC, whereas some of them were prognostic in buccal carcinoma (Sathyan ; Trivedi ). The histological structures and the carcinogenesis are different in buccal mucosa and in oral tongue. Similarly, the base portion and the oral (mobile) portion of the tongue have differences in the etiopathogenesis of the cancer. In the base of the tongue, HPV is commonly linked with the cancer, whereas the virus is rarely founded in mobile tongue. However, still several studies combine oral and base of the tongue SCC samples, or they do not specify for which part of the tongue the analyses were done (Supplementary Tables 1, 2, 4 and 5). To avoid the effect of tumour heterogeneity, we included in the present meta-analyses studies in which the authors defined their cohort as OTSCC. Our meta-analysis and conclusions are based on those studies which included previously untreated, surgically resected, primary OTSCC.

Almost all of the most reported biomarkers (except VEGFs, which contribute to tumour angiogenesis) reflect important growth-related properties of the cancer cells. However, non-neoplastic cells of the tumour stroma, including fibroblasts, endothelial cells and inflammatory cells, seem also to have a critical role in cancer progression (Marsh ). Accordingly, biomarkers of the stromal microenvironment might even have a greater impact on prognosis than biomarkers related to tumour cells (Marsh ). However, biomarkers related to tumour microenvironment such as Activin A (Kelner ) are not yet widely studied in OTSCC. The only tumour microenvironment biomarker that was reported repeatedly in OTSCC, is the cancer-associated fibroblast identified by a α-smooth muscle actin antibody (Supplementary Table 6). Further studies should also focus on the evaluation of promising biomarkers of tumour microenvironment, such as fibronectin and tenascin-C (Sundquist ), in addition to biomarkers related to cancer cells.

Notably, 152 biomarkers (83%) have been studied only once or twice, so it is not possible to reach trustworthy conclusions on the basis of such limited evidence. About 86% of the studies claimed to have found at least one biomarker to have prognostic value. Negative or ‘non-significant’ prognostic finding as the only result was not widely published (about 14%). This might well indicate a fair amount of publication bias as noted by Soland and Brusevold (2013). However, the articles which reported new promising biomarkers have also usually evaluated the prognostic significance of previously known biomarkers (e.g., Ki-67, p53, and p16). This approach helps to validate previously published data and allows the accumulation of evidence for known biomarkers. Relationships between biomarkers and clinicopathologic manifestations, other than survival, have also been studied in several publications (Albert ; Wang ). However, even though it is important to understand such relationships, clinically the most relevant information is provided by proper survival analysis of the tested biomarker.

REMARK guidelines (Altman ) have suggested items to be reported in prognostic studies of tumour markers. Tables 2, 3, 4, 5, 6 show that many studies on p53, Ki-67, p16, VEGFs and cyclin D1 did not follow REMARK criteria accurately. This indicates shortcomings in the reporting of the biomarkers tested in these studies, and subsequently limits the possibilities to reach definitive conclusions of their usefulness.

To date, none of the most often reported biomarkers can be recommended as prognosticators valid for clinical use. This may be related to the mixing of early and late stages of OTSCC in the same analysis, the small sizes cohorts or numbers of events seen in many studies. In addition, antibodies supplied by different manufacturers, variable staining conditions, and different cut-off values were also seen. All of these factors might affect the biomarker results published.

Even though it is highly desirable to apply multiple Cox regression or similar predictive modelling (‘multivariate’ analysis) to adjust for important classical prognostic factors (like TNM stage and histologic grade) and relevant patient characteristics (like age), limited cohort size does not allow conducting a regression analysis of good predictive performance, and might affect the validity of the estimation results (Ogundimu ). Moreover, numerous studies which included small cohorts (and small numbers of events) have reported regression analyses with many variables, and such studies may run the risk of overfitting. A rule of thumb of 10 outcome events for each prognostic term involved in multiple Cox regression is widely advocated (Peduzzi ; Ogundimu ), but in many instances at least 20 events per term would be needed for reliable modelling (Ogundimu ). It seems that no consideration has been devoted to this issue in most of the studies applying Cox regression. Another common statistical shortcoming was to report the P value only and to lean on it when making judgments of the prognostic value of a biomarker depending on whether the result was ‘significant’ (when P<0.05) or ‘non-significant’ (when P⩾0.05). Such statistical malpractice meant that many studies were not eligible for inclusion in the current meta-analyses; and may have caused false inferences in some studies, thus it should be avoided in future studies. More useful results can be obtained by reporting the estimated HR, preferably an adjusted one, and its CI, as is recommended by REMARK guidelines (Altman ). When predicting mortality, it was quite common to report OS but not cancer-specific mortality. Analysis by cause of death would, however, provide more relevant prognostic information than that of OS only (Läärä ). For better statistical analysis, guidelines for improving statistical reporting (Greenland ) should be followed, the authors should involve a statistician, and journals should assign an experienced statistician as one of their reviewers.

Recently, the digital evaluation of biomarkers has been suggested to facilitate the evaluation and to avoid inter- and intra-observer variability (Bouzin ). So far only a few studies have used this technique in the prognostication of OTSCC (Hannen , 2002), and none of the studies in the meta-analysis had applied digital scorings. In the future, image analysis will most likely be more commonly applied also for OTSCC scorings.

Regarding the global distribution of the studies, the Japanese population was widely represented in studies of p53, Ki-67, VEGFs and cyclin D1 (Table 7). However, the studies included in our meta-analyses (Figures 3 and 4) were conducted in a different population (marked with an asterisk in Table 7).

Table 7

Distribution of OTSCC immunohistochemical prognostic biomarkers studies according to countries

	Number of studies (HR and CI were reported in the study for at least one survival endpoint)
Country	p53	Ki-67	p16	VEGFs	Cyclin D1
Japan	4 (2)a	4 (2)	2 (1)a	5 (3)a	5 (0)
USA	4 (0)	1 (1)	2 (1)a	–	–
India	3 (2)a	–	4 (2)b	–	2 (1)a
China	2 (0)	–	1 (0)	1 (1)a	1 (0)
Korea	1 (0)	1 (1)	–	3 (0)	–
Australia	1 (1)a	1 (1)a	2 (2)b	–	1 (1)a
Finland	1 (0)	1 (1)a	–	–	–
Iran	–	–	–	–	–
Norway	1 (0)	1 (0)	–	–	–
Sweden	1 (1)	1 (0)	1 (0)	–	–
Tunisia	1 (0)	1 (0)	–	–	–
Turkey	1 (0)	–	–	–	–
Brazil	–	2 (0)	–	–	1 (0)
Czech	–	–	1 (0)	1 (1)a	–
Spain	–	–	1 (0)	–	–
USA & China	–	1 (1)a	–	–	–
Finland & Israel	–	1 (1)	–	–	–
Total number of studies	20	15	14	10	10

The entries are total numbers of studies and in parentheses are given numbers of studies in which the hazard ratio and its 95% confidence interval were reported for at least one survival endpoint.

Only one study from this country was included in the meta-analysis of overall survival.

Two studies from India and two studies from Australia were included in the meta-analysis of p16.

The primary treatment was surgical OTSCC resection in all of the publications included in the meta-analyses. However, one of the p16 studies included cases with either surgical excision or external beam radiotherapy for patients unwilling/unfit for surgery (Ramshankar ). None of the studies in the meta-analyses mentioned any preoperative therapy, which could have led to molecular changes in the resected tumour tissue. Neck dissection, radiotherapy and/or chemotherapy are used for selected cases of OTSCC, especially those in advanced stages. Such variations in treatment modalities must have an impact on the prognosis, as well as on the prognostic analyses. Unfortunately, therefore a limitation in our meta-analyses is that the cases included were not homogenous with regard to neck dissection or other adjuvant therapies. The results of our meta-analyses must, in any event, be interpreted with due caution because of various shortcomings in the published studies and their reporting. First, for each biomarker very few original studies were available that provided the necessary statistical data. Second, essential heterogeneity was evident within and across individual studies with respect to the spectrum of patients, staining methods, and so on. Third, the numbers of patients and outcome events were mostly small implying poor statistical precision. Fourth, only crude HR estimates, unadjusted for classical prognostic factors were typically reported, leaving room for an unknown amount of bias in the assessment of the prognostic power of the pertinent marker due to uncontrolled confounding.

Conclusion

The identification of biomarkers may increase the possibility of detecting cancers with high risk for poor prognosis. Here, for the first time, we systematically reviewed the literature for immunohistochemical prognostic biomarkers for OTSCC. On the basis of our meta-analysis VEGF-A and cyclin D1 may be useful prognostic biomarkers for OTSCC. Therefore, these biomarkers, in addition to the new biomarker candidates, should be further evaluated carefully following the REMARK criteria using large, preferably multicentre, cohorts of OTSCC cases separated into early and late stages.