Abstract

p16, encoded by the CDKN2A gene, is a tumor suppressor that has been widely studied in cancer research. However, the relationship of p16 with prognostic and clinicopathological parameters in patients with bladder cancer remains unclear. Data inclusion criteria were articles reporting on the relationship between p16 expression and the prognosis or clinicopathology in patients with bladder cancer. Meta-analyses were performed with Stata software. Hazard ratios (HRs) or odds ratios (ORs) and 95% confidence intervals (95% CI) were calculated to evaluate the relative risks. The source of heterogeneity was analyzed by subgroup analysis. A total of 37 studies with 2246 cases were included and analyzed. The results identified an important link between downregulated p16 expression and poor prognosis in patients with bladder cancer in terms of recurrence-free survival (RFS), overall survival (OS), progression-free survival (PFS), and some clinicopathological parameters including clinical staging, pathological degree, and lymph node metastasis. Subgroup analysis also showed that low p16 expression could function as a warning sign for RFS and PFS in patients with early-stage (Ta–T1) bladder cancer. In conclusion, p16 might play an essential role in the deterioration of bladder cancer and could serve as a biomarker for the prediction for patients’ progression and prognosis.

1. Introduction

Bladder cancer is the most frequent malignancy of the urinary tract and the ninth most common cancer worldwide [1]. About 95% of bladder cancers are histologically transitional cell carcinoma, with rare cases of squamous cell carcinoma and adenocarcinoma. However, the pathogenesis of bladder cancer is still unclear, and its occurrence and development appear to be affected by multiple genes [2]. Serrano et al. first cloned the cDNA of the gene encoding the tumor suppressor protein p16 (CDNK2A) in 1993; since then it has been widely studied in the field of cancer research [3]. Previous studies have reported ubiquitous downregulation of p16 gene expression in bladder cancer, as a result of various alterations including complete deletion, point mutation, or promoter methylation [46]. Furthermore, p16 could compete with cyclin D1 for binding to Cyclin Dependent Kinase (CDK) 4/6, thus blocking the phosphorylation of retinoblastoma (Rb) protein and inhibiting release of the transcription factor E2F, preventing cell conversion from G1 phase to S phase, and eventually suppressing cell proliferation. These results suggest that abnormal expression of the p16 gene in cells might be associated with tumorigenesis [6, 7].

Numerous studies to date have explored the clinicopathological and prognostic significance of p16 in patients with bladder cancer. However, as a result of differences in sample sizes, accuracies of the statistical data, study populations, and interventions, the results remain inconclusive, and evidence-based confirmation by large-scale clinical trials is still lacking. We therefore conducted an in-depth systematic review and meta-analysis to investigate the correlation between abnormal expression of p16 and clinicopathological features, as well as prognosis in patients with bladder cancer. The specific mechanisms are shown in Figure 1.

2. Materials and Methods

2.1. Literature Search

The terms and combinations including “Cyclin Dependent Kinase Inhibitor p16,” “CDKN2A Protein,” “p16INK4A Protein,” “MTS1 Protein,” “Cyclin Dependent Kinase Inhibitor 2A,” “Multiple Tumor Suppressor 1,” “Cdk4 Associated Protein p16,” “TP16,” and “urinary bladder neoplasms,” “bladder tumors,” “bladder cancers,” “bladder carcinomas,” and “,” “,” “follow-up,” “mortality,” “predict,” “course,” “outcome,” and “clinicopathological” were used to search the following domestic and international databases: PubMed, Wiley Online Library, Embase, Cochrane Central Register of Controlled Trials, Science Direct, EBSCO, Google Scholar, Ovid, LILACS, China National Knowledge Infrastructure (CNKI), China Biology Medicine disc (CBMdisc), CQVIP, and Wan Fang, with unified retrieval rules such as Boolean logic. The obtained search results were then analyzed, evaluated, reviewed, and manually screened to determine their relevance.

2.2. Inclusion and Exclusion Criteria

Inclusion criteria were as follows: (1) patients diagnosed with bladder cancer; (2) immunohistochemical (IHC) detection of p16 expression levels in the tissues; (3) relationships between abnormal expression of p16 and prognostic indicators such as recurrence-free survival (RFS), progression-free survival (PFS), and overall survival (OS) or associations between p16 and clinicopathological features that were evaluated; (4) hazard ratio (HR), odds ratio (OR), relative risk (RR), and 95% confidence intervals (CI) that could be obtained directly from the full article or indirectly calculated with relevant software based on the data provided in the graphics and tables; (5) only the newest studies or the ones with higher quality were retained if the data were repeated in different studies; and (6) studies in English or Chinese.

Exclusion criteria were as follows: (1) cell or animal studies, case reports, letters, reviews, and meta-analyses; (2) articles with similar content or using the same data or those with small sample sizes () and those with no directly or indirectly extractable HR, OR, and 95% CI data; and (3) articles that could not be understood because of language barriers.

2.3. Data Extraction

Two independent investigators (Xiaoning Gan and Rongquan He) reviewed the articles that met the criteria and extracted data on author, year of publication, nationality, sample size, patient age, detection method of p16, antibody source and dilution, clinical stage, pathological degree, other costudied prognosis-associated genes, cut-off value, outcome, and extraction method of the study subjects. Discrepancies between the two independent investigators in terms of data extraction were resolved by discussion among all the authors.

2.4. Statistical Analysis

Effects of p16 on the related prognostic indexes were detected by merging the HRs and 95% CI of the included literatures, which were evaluated through the Forest plot and related parameters after the merging. The HRs and 95% CI values mainly came from direct extraction of the original text or survival curve through extraction and calculation by software.

The relationships between p16 and the clinicopathological parameters were derived from the binary variable data extracted from the original articles. ORs and 95% CI values came from the binary variable data calculated by Stata software. The data were then combined, and their statistical significance was evaluated by Forest plot and related parameters, to clarify the relationship between p16 low-expression and clinicopathological parameters.

Heterogeneity was measured by statistics as follows: no heterogeneity: ; low heterogeneity: ; moderate heterogeneity: ; high heterogeneity: . If and , a fixed-effect model would be used in combination with HRs, ORs, and 95% CI; if and , then a random-effect model would be selected. Heterogeneity analysis was performed to assess the accuracy of the data, and subgroup and sensitivity analyses were carried out based on professional knowledge.

Publication bias was detected by Begg’s funnel plot and Egger’s test with Stata software. A two-sided value < 0.05 was considered to indicate statistical significance. Statistical analyses were carried out with StataSE 12.0, Engauge, Photoshop CS5, and Microsoft Office 2007.

3. Results

3.1. Eligible Studies

A total of 364 articles were identified from the databases, including 190 English and 174 Chinese articles, 222 of which were excluded because of discrepancies between the study theme and their abstracts. The full text of the remaining 142 articles was then reviewed for their fit with the current study, after which a further 105 articles were excluded because they met one or more of the exclusion criteria, such as the cell or animal studies, reviews, and letters and studies with identical data and no extractable HR, OR, and 95% CI data from the full text or language barrier. The remaining 37 articles [4, 5, 842] with 2246 cases were included in our study and consisted of 21 English [4, 5, 819, 2125, 41, 42] and 16 Chinese [20, 2640] articles. The screening process was demonstrated in Figure 2.

The basic features of the included studies were presented in Table 1. Among the 37 articles, 26 studies [4, 5, 821, 2834, 39, 41, 42] investigated the relationship between low expression of p16 and prognostic parameters in bladder cancer patients (RFS, OS, PFS, and DSS/CSS), and 30 studies [4, 1013, 1519, 2140] assessed the association between p16 and clinicopathological factors in patients with bladder cancer.

3.2. Relationship between Downregulated p16 Expression and RFS in Patients with Bladder Cancer

A total of 17 studies with 1032 subjects were included in the final analysis of RFS [4, 5, 8, 10, 1315, 19, 21, 2834, 39]. Low expression of p16 was related to poor RFS in patients with bladder cancer (HR = 1.63, 95% CI = 1.36~1.94, and ), with low observed heterogeneity (, ) (Figure 3(a)).

Cumulative meta-analysis based on year of publication and sample size demonstrated that the results tended to stabilize with increasing sample size, but there was no obvious relationship between the results and year of publication.

Based on sensitivity analysis, the study by Yang et al. [13] was initially excluded because of a large difference in HR compared with the overall average, which was attributed to the selection of a different calculation method in the original article. Binary variable data were extracted and the HR and 95% CI were therefore recalculated with Stata software.

Subgroup analysis based on geographic region showed that low expression of p16 was associated with RFS in patients with bladder cancer both in Asia (HR = 1.44, 95% CI = 1.15~1.81, and ) and in Europe (HR = 1.90, 95% CI = 1.13~3.19, and ). The results of American studies (HR = 1.58, 95% CI = 0.77~3.25, and ) need to be confirmed with larger sample sizes. The heterogeneity of Asian studies (, ) was lower than the overall heterogeneity, while that of Europe (, ) was higher, calculated with the random-effect model.

Subgroup analysis based on clinical stage suggested that the effect of p16 on RFS was associated with clinical stage (Tis-T1 group: HR = 1.96, 95% CI = 1.23~3.14, and ; , ; Tis–T4 group: HR = 1.41, 95% CI = 1.12~1.77, and ; , ).

Subgroup analysis based on histopathological grade showed that heterogeneity decreased from G1-G2 (HR = 4.12, 95% CI = 2.48~6.83, and ; , ), G1–G3 (HR = 1.44, 95% CI = 1.18~1.75, and ; , ), and G2-G3 (HR = 1.37, 95% CI = 0.78~2.42, and ; , ), indicating that the effects of p16 on RFS in patients with bladder cancer were closely associated with pathological grade.

Subgroup analysis also showed an effect of cut-off value on the influence of p16 on RFS (cut-off value ≤ 10%: HR = 1.83, 95% CI = 1.34~2.51, and ; , ; cut-off > 10%: HR = 1.34, 95% CI = 0.86~2.09, and ; , ).

In addition, subgroup analysis of early-stage data from 430 subjects from eight studies also demonstrated that low expression of p16 significantly affected RFS in patients with early-stage (Ta–T1) bladder cancer (HR = 1.96, 95% CI = 1.23~3.14, and ; , ).

3.3. Relationship between the Low Expression of p16 and OS in Patients with Bladder Cancer

A total of 425 subjects in eight studies were included in the final analysis of OS [9, 11, 14, 16, 19, 20, 30, 34], which showed that low expression of p16 was associated with decreased OS in patients with bladder cancer (HR = 1.70, 95% CI = 1.16~2.50, and ), with no significant observed heterogeneity (, ) (Figure 3(b)).

Cumulative meta-analysis and sensitivity analysis indicated relatively low overall heterogeneity and no study with high sensitivity.

Subgroup analysis based on geographic area showed a subtle distinction between p16 expression and OS in patients with bladder cancer in Asia (HR = 1.61, 95% CI = 0.97~2.66, and ; , ) and Europe (HR = 2.54, 95% CI = 1.05~6.15, and ; , ).

Subgroup analysis was also performed based on clinicopathological stages. However, limitations of sample size led to the impossibility of determining if the effects of p16 expression on OS were associated with these parameters in patients with bladder cancer (Ta–T1 group: HR = 1.57, 95% CI = 0.32~7.75, and ; , ; Ta–T4 group: HR = 1.59, 95% CI = 0.98~2.60, and ; , ; T2–T4 group: HR = 1.96, 95% CI = 0.99~3.88, and ; , ; low-grade group: HR = 1.41, 95% CI = 0.18~10.90, and ; , ; G1–G3 group: HR = 1.82, 95% CI = 1.16~2.84, and ; , ; high-grade group: HR = 1.41, 95% CI = 0.62–3.19, and ; , ).

Subgroup analysis based on cut-off value indicated that the effects of p16 on OS in patients with bladder cancer were associated with cut-off value (cut-off value ≤ 10%: HR = 1.83, 95% CI = 1.17~2.86, and ; , ; cut-off value > 10%: HR = 1.40, 95% CI = 0.66~2.96, and ; , ).

3.4. Relationship between Low Expression of p16 and PFS in Patients with Bladder Cancer

A total of 470 subjects in seven studies were included in the ultimate analysis of PFS [5, 9, 10, 12, 19, 41, 42]. The results showed a correlation between low expression of p16 and poor PFS in patients with bladder cancer (HR = 2.18, 95% CI = 1.37~3.48, and ), with low heterogeneity detected (, ).

Cumulative meta-analysis revealed no obvious characteristics because of the limited range of publication dates and the sample sizes.

Sensitivity analysis identified two studies [41, 42] as having the highest heterogeneities. Further investigation revealed that this heterogeneity was caused by different methods of measuring p16 (fluorescence in situ hybridization) and studying the influence of hemizygous or homozygous deletion of p16 on patient prognosis. These two studies were finally excluded because of their incompatible study objectives, leaving a total of 347 subjects from six studies in the final analysis of PFS. The results showed that low expression of p16 was correlated with poor PFS in patients with bladder cancer, and the heterogeneity was eliminated (HR = 1.84, 95% CI = 1.13~3.01, and ; , ) (Figure 3(c)).

Despite a reduced sample size, subgroup analysis of the 347 subjects from five studies [5, 9, 10, 12, 19] demonstrated that the effects of p16 expression on PFS were affected by clinical stage (Ta–T1 group: HR = 2.09, 95% CI = 1.21~3.63, and ; , ; T2–T4 group: HR = 1.14, 95% CI = 0.39~3.31, and ; , ) and geographical location (Europe: HR = 2.09, 95% CI = 1.21~3.63, and ; , ; America: HR = 1.14, 95% CI = 0.39~3.31, and ).

Subgroup analysis based on cut-off value demonstrated some relationship between cut-off value and the influence of p16 expression on PFS (cut-off value ≤ 10%: HR = 2.61, 95% CI = 1.42~4.77, and ; , ; cut-off value > 10%: HR = 0.95, 95% CI = 0.41~2.18, and ; , ).

The results from 297 subjects with early-stage (Ta–T1) bladder cancer from four studies [5, 10, 12, 19] suggested that low expression of p16 was also significantly associated with poor PFS in early-stage bladder cancer (HR = 2.09, 95% CI = 1.21~3.63, and ; , ).

3.5. Relationship between Low Expression of p16 and DSS/CSS in Patients with Bladder Cancer

A total of 187 subjects from three studies were included in the DSS/CSS analysis [9, 17, 18]; limitation of the sample size caused the impossibility of demonstrating an association between low expression of p16 and DSS/CSS (HR = 1.52, 95% CI = 0.85~2.71, and ; , ).

3.6. Relationship between Low Expression of p16 and Clinicopathological Parameters in Patients with Bladder Cancer

The relationship between low expression of p16 and clinicopathological parameters [4, 1013, 1519, 2140] was further explored by analysis of 30 studies including 1785 subjects. The results of statistical analyses were as follows: T2–T4/Ta–T1: OR = 3.13, 95% CI = 2.42~4.06, and ; , ; T1/Ta: OR = 1.55, 95% CI = 0.87~2.76, and ; , ; G3/G1-2 [43]: OR = 3.33, 95% CI = 2.51~4.42, and ; , ; and H/L [44]: OR = 1.20, 95% CI = 0.69~2.33, and ; , ; because of the high heterogeneity, a random-effects model was therefore applied. Meanwhile, these results demonstrated significant differences in the effects of low p16 expression in patients with bladder cancer between the two WHO clinical pathological grading methods in 1973 and 2004.

Analysis of the results for lymph node metastasis showed OR = 2.20, 95% CI = 1.26~3.83, and ; , . The small sample size caused the impossibility of demonstrating any significant influence of pathological parameters such as muscle invasion, tumor number (multiple/single), and tumor size on the effect of p16 expression (Table 2).

3.7. Retrospective Review

Three studies [6, 24, 25] were retrospectively reviewed because of differences between their prognosis results and the data required by the meta-analysis. As shown in Table 3, low expression of p16 was associated with poor prognosis in patients with bladder cancer. However, some of the values were <0.05 because of the small sample sizes.

3.8. Publication Bias

Publication bias was detected by Begg’s funnel plot and Egger’s test (Figure 4). The points representing studies were symmetrically arranged in a funnel shape in the funnel plot, and the values calculated from Egger’s test with higher detection effectiveness were >0.05, indicating no publication bias. The only exception was for RFS; the funnel plot was asymmetrical and with a few points outside the funnel. Publication bias was also detected by Egger’s test (G1–G3 group: ; Asia group: ), indicating potential publication bias in terms of RFS.

4. Discussion

p16, also known as tumor suppressor gene I (multiple tumor suppressor, MTS I), is located in 9p21 and is composed of two introns and three exons [45]. It is a key gene in cell cycle regulation, with its expression product being involved in the negative regulation of cell proliferation. Studies have shown that downregulation of p16 gene expression resulted in the loss of its inhibitory effects on CDK4/CDK6, which in turn may lead to malignant, abnormal cell proliferation and accelerated tumor development [7, 46, 47]. Elucidation of the relationship between low expression of p16 and prognosis and clinicopathology in patients with bladder cancer is therefore important for its early diagnosis, treatment, and prognosis.

Pan et al. performed a meta-analysis of the prognostic significance of abnormal p16 and p21 expression in bladder cancer in 2006 [48]. However, the current study analyzed a larger sample size; Pan et al.’s study included 12 articles with 975 cases, compared with 37 articles and 2246 cases in our study, leading to more accurate and reliable results. Secondly, Pan et al.’s study involved a number of mixed factors with no clear listing of each prognostic index or subgroup discussion. In contrast, the current study included subgroup analyses for the different indicators including RFS, PFS, OS, and DSS/CSS, allowing more thorough insights into the relationships between p16 expression and the prognostic and clinicopathological parameters in bladder cancer patients. Thirdly, Pan et al. found no association between p16 expression and prognosis in early Ta–T1 stage (stage I) bladder cancer, possibly because of the omission of the study by Krüger et al. [5], which explored the significance of p16 as an independent tumor predictive factor for the development of T1 bladder cancer, and demonstrated the important clinical value of low p16 expression in the early diagnosis and prognosis of patients with early-stage bladder cancer.

The current study systematically analyzed the relationships between p16 expression and prognostic index and clinicopathological parameters in patients with bladder cancer and showed that low expression of p16 was closely correlated with poor prognosis (Figure 5). However, the included studies varied in terms of study subjects, design, sample size, interventions, outcomes, time of study, and publication date. We used cumulative meta-analysis, sensitivity analysis, and subgroup analysis to explore the effects of the main variables in the included studies. Overall, the results confirmed that the relationship between low expression of p16 and prognosis in patients with bladder cancer was affected by clinicopathological stage, geographic origin of the study subjects, detection method, and cut-off values. Based on these findings, we further analyzed the relationships between p16 expression and clinicopathological parameters and demonstrated associations between low expression of p16 and clinical stage and lymph node metastasis, implying that the p16 gene tended to exert its regulatory effects during the early stage of bladder carcinogenesis. Low expression of p16 was also correlated with poor PFS and RFS in early-stage (Ta–T1) bladder cancer. These results thus confirmed an important role for p16 in the occurrence and development of bladder cancer. Meanwhile, through Phase I and II clinical trials, studies have revealed that CDK4/6 is an attractive target in p16 related pathway for anticancer therapy [4951]. Furthermore, previous study also suggested that p16 functional peptide, as a molecular targeting agent, showed effective reactions for the treatment of renal cell carcinoma [52]. The results of these researches and our current meta-analysis had the effect of mutual authentication. Therefore, a better understanding of the mechanism underlying the development and progression of bladder cancer may play a significant role in prevention, target therapy, and prognosis, particularly if more sensitive and specific correlative biomarkers can be discovered and verified.

The current study had some limitations. First, tumors are the result of both environmental and genetic factors, and p16 may thus be only one of several factors involved in the whole process of bladder carcinogenesis. Secondly, heterogeneity may result from differences in intervention measures (surgery, radiotherapy, chemotherapy, or combination), immunohistochemical techniques (different antibodies, evaluation standards, etc.), and the HR extraction methods used in the included studies. Finally, the exclusion of articles because of language barriers and of studies that were not published because of a lack of sufficient data may have led to potential publication bias.

In conclusion, the results of the current study provide evidence for a relationship between p16 expression and prognosis and clinicopathological features in patients with bladder cancer. The results of this meta-analysis will help to inform about the development of clinical guidelines promoting best medical care for patients with bladder cancer. Further studies are required to investigate the combined influence of genetic and environmental factors on the development and progression of bladder cancer.

Abbreviations

RFS:Recurrence-free survival
OS:Overall survival
PFS:Progression-free survival
DSS:Disease-specific survival
CSS:Cancer-specific survival
HR:Hazard ratio
OR:Odds ratio
CI:Confidence interval.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Xiaoning Gan, Xiaomiao Lin, Rongquan He, Xinggu Lin, Hanlin Wang, Liyan Yan, Hong Zhou, and Hui Qin performed the literature search, data extraction, and statistical analysis and drafted the paper. Gang Chen supervised the literature search, data extraction, and analysis and reviewed the paper. Xiaoning Gan, Xiaomiao Lin, Rongquan He, Xinggu Lin, Hanlin Wang, Liyan Yan, Hong Zhou, Hui Qin, and Gang Chen read and approved the final paper.