Cyclooxygenase-2 Expression in Bladder Cancer and Patient Prognosis: Results from a Large Clinical Cohort and Meta-Analysis

Abstract

Aberrant overexpression of cyclooxygenase-2 (COX2) is observed in urothelial carcinoma of the bladder (UCB). Studies evaluating COX2 as a prognostic marker in UCB report contradictory results. We determined the prognostic potential of COX2 expression in UCB and quantitatively summarize the results with those of the literature through a meta-analysis. Newly diagnosed UCB patients recruited between 1998–2001 in 18 Spanish hospitals were prospectively included in the study and followed-up (median, 70.7 months). Diagnostic slides were reviewed and uniformly classified by expert pathologists. Clinical data was retrieved from hospital charts. Tissue microarrays containing non-muscle invasive (n = 557) and muscle invasive (n = 216) tumours were analyzed by immunohistochemistry using quantitative image analysis. Expression was evaluated in Cox regression models to assess the risk of recurrence, progression and disease-specific mortality. Meta-hazard ratios were estimated using our results and those from 11 additional evaluable studies. COX2 expression was observed in 38% (211/557) of non-muscle invasive and 63% (137/216) of muscle invasive tumors. Expression was associated with advanced pathological stage and grade (p<0.0001). In the univariable analyses, COX2 expression - as a categorical variable - was not associated with any of the outcomes analyzed. As a continuous variable, a weak association with recurrence in non-muscle invasive tumors was observed (p-value = 0.048). In the multivariable analyses, COX2 expression did not independently predict any of the considered outcomes. The meta-analysis confirmed these results. We did not find evidence that COX2 expression is an independent prognostic marker of recurrence, progression or survival in patients with UCB.

About the Authors:

Maciej J. Czachorowski

Affiliation: Genetic and Molecular Epidemiology Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain

André F. S. Amaral

Affiliation: Genetic and Molecular Epidemiology Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain

Santiago Montes-Moreno

Affiliation: Servicio de Anatomía Patológica, Hospital Universitario Marqués de Valdecilla, Santander, Spain

Josep Lloreta

Affiliations Institut Municipal d’Investigació Mèdica – Hospital del Mar, Barcelona, Spain, Departament de Patologia, Hospital del Mar – Parc de Salut Mar, Barcelona, Spain, Departament de Ciències Experimentals i de la Salut, Universitat Pompeu Fabra, Barcelona, Spain

Alfredo Carrato

Affiliations Hospital General Universitario de Elche, Elche, Spain, Hospital Ramon y Cajal, Madrid, Spain

Adonina Tardón

Affiliation: Universidad de Oviedo, Oviedo, Spain

Manuel M. Morente

Affiliation: Tumor Bank Unit, Spanish National Cancer Research Centre (CNIO), Madrid, Spain

Manolis Kogevinas

Affiliation: Centre de Recerca en Epidemiologia Ambiental (CREAL), Barcelona, Spain

Francisco X. Real

Affiliations Departament de Ciències Experimentals i de la Salut, Universitat Pompeu Fabra, Barcelona, Spain, Epithelial Carcinogenesis Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain

Núria Malats

* E-mail: [email protected]

Affiliation: Genetic and Molecular Epidemiology Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain

for the SBC/EPICURO investigators

¶Membership of the Spanish Bladder Cancer (SBC)/EPICURO Study investigators is provided in the Acknowledgments.

Introduction

Urothelial carcinoma of the bladder (UCB) is the most common bladder cancer type in developed nations [1]. UCB predominantly manifests (70–80% of patients) as a non-muscle invasive tumor (NMIBC: pTa-pT1) characterized by an overall good prognosis following transurethral resection in patients with low-grade tumors (pTaG1/2), and intravesical chemotherapy and/or Bacillus Calmette Guerin (BCG) instillation in patients with high-grade tumors (pTaG3 or pT1G2/3) [2]. Approximately 70% of NMIBC patients suffer a recurrence following treatment and a further 15% progress, developing new tumors exhibiting muscle invasion (MIBC: pT2-pT4); the risk of progression being higher among patients with high-grade tumors [2]. Due to a high rate of recurrence and the need for close follow-up over a patient’s lifetime, UCB remains one of the most expensive tumors to treat on a per patient basis [3]. A lower proportion (20–30%) of UCB patients are diagnosed with muscle invasive tumors (MIBC; pT2-pT4) characterized by poor prognosis: 50% of these patients die from their cancer [2]. Genomic profiling and gene expression analyses indicate a strong correlation between these pathologic classifications and the underlying molecular architecture of UCB [4].

Growing evidence indicates that chronic inflammation may increase the risk of UCB [5]. Studies investigating the prolonged use of cyclooxygenase-2 (COX2) inhibiting non-steroidal anti-inflammatory drugs (NSAIDs) have reported a decrease in UCB risk [6], [7]. COX2 is a prostaglandin endoperoxide synthetase that catalyzes the production of prostanoids upon induction by proinflammatory cytokines, growth factors, tumor promoters and other external stimuli [8]. COX2 activation mediates cellular processes also implicated in carcinogenesis such as angiogenesis, cell survival/proliferation and apoptosis [9]. Moreover, studies have shown that bladder tissue from patients with cystitis or UCB exhibits elevated COX2 levels in contrast to benign bladder tissue [10], [11].

While numerous groups have investigated the prognostic potential of COX2 expression in UCB [12]–[31], there is no clear consensus on its utility. The objective of this study was to assess whether COX2 protein expression in UCB cells is associated with prognosis using a large and standardized cohort of newly diagnosed bladder cancer patients. A meta-analysis was also done to summarize these results together with those from other studies published on the topic.

Materials and Methods

Study Population

A total of 773 newly diagnosed UCB cases aged 22–80 years (mean ± SD = 66±10 yrs) with a median follow-up of 70.7 months (range 0.7–117.7 months) and available tumor tissue were used in the current analysis. All cases were recruited between 1998 and 2001 from 18 hospitals in five regions of Spain as part of the Spanish Bladder Cancer (SBC)/EPIdemiology of Cancer of the UROthelium (EPICURO) study, a hospital-based case-control study described previously [32]. A pathologist review panel uniformly classified the T stage and grade (G) of each tumor biopsy according to the criteria of the TNM classification and the WHO-ISUP [33], using the three grade redefinition provided by the WHO [34], [35]. All bladder tumor samples used in the study were collected prior to the administration of any intravesical or systemic therapy. Clinical information related to diagnostic procedures, tumor characteristics and treatment was collected from medical records, and a computerized questionnaire was used for the collection of sociodemographic data. NMIBCs were removed by transurethral resection and patients received intravesical chemo- or immunotherapy (i.e. BCG) as appropriate. The majority of patients presenting with MIBCs were treated by radical cystectomy; in cases where surgery was not possible, radiotherapy or systemic chemotherapy were administered. Follow-up information was collected annually from hospital records and through direct telephone interviews by trained monitors using structured questionnaires. Among NMIBCs, recurrence was defined as the appearance of a new NMIBC following a previous negative follow-up cystoscopy, and progression, as the development of a MIBC. In patients initially presenting with MIBCs, any tumor reappearance after treatment was considered progression, regardless of whether the tumor relapse was local or distal. Tumour-specific survival was assessed only for patients with MIBCs. Informed written consent was obtained from study participants in accordance with the Ethics Committees of each participating hospital.

Immunohistochemistry

Tissue blocks of formalin-fixed, paraffin-embedded primary bladder tumors were used to construct tissue microarrays (TMA) containing tumor cores of 0.6-mm in diameter represented in duplicate and selected from the most representative regions of the tumor on which T and G were based. After deparaffinisation and heat-induced antigen retrieval, all slides were stained simultaneously at the Histology and Immunohistochemistry Core Unit of the CNIO using the PT LINK system as per manufacturer’s instructions (Dako Inc., Glostrup, Denmark). Briefly, tissue sections were incubated with anti-COX2 rabbit monoclonal antibody (ThermoFisher Scientific, Fremont, CA, USA; #RM-9121-R7; pre-diluted, ready-to-use) at room temperature, followed by visualization using the EnVision Flex Visualization system (Dako Inc., Glostrup, Denmark) and exposure to diaminobenzedine. Tissues were then counterstained with haematoxylin, dehydrated and mounted. A section of colon tissue was used as a positive control.

Evaluation of COX2 Immunostaining

COX2 expression was quantified using the Ariol SL-50 (version 3.1.2, Applied Imaging Corp., San Jose, CA, USA) high-throughput slide imaging scanner. All cores were imaged and processed using a light microscope and the accompanying TMA Multistain Imaging software. The program was trained by a pathologist (SM) to maximize the inclusion of positively stained tumor epithelium while minimizing stromal material, as described previously [36]. COX2 expression score was calculated as the product of the mean intensity of staining (by defining the background and saturation limits of the antibody and imaging sensor, respectively) and the proportion of cellular antibody-positive area divided by total cellular area. Values from replicate cores were averaged to provide a final expression score for each patient. Furthermore, one randomly selected TMA (representing 10% of all cores) was analyzed by direct visual microscopic inspection by an independent pathologist (MMM) to enable comparison with the automated scoring approach. The pathologist-derived score was calculated as the product of COX2 staining intensity (1 = weak, 2 = intermediate, 3 = strong) and a quartile of the percentage of epithelial tumor cells stained (0–4; with 0 representing 0% staining), providing a final categorical score in the range of 0–12. There was a high and significant correlation between the machine and pathologist derived scores (Spearman rho = 0.85; 95%CI = 0.79–0.90; p-value<0.00001). COX2 expression was analyzed as both a continuous variable and categorical variables partitioned at the median and extreme tertiles. Additionally, expression was examined as a categorical variable dichotomized at a threshold (0.340 arbitrary units [au]) above which COX2 expression was considered to be positive. This expression threshold was derived by comparing the pathologist’s (MMM) binary assignment of positive expression (i.e. score of 0 vs. score ≥1, as described above) to the machine-derived continuous score using receiver operating characteristic (ROC) curve analysis (area under the curve = 0.95; 86% sensitivity and 92% specificity) [37].

Meta-analysis

The meta-analysis included COX2 expression results from our own series (using the ROC-derived categorical expression variable) and relevant studies published before 1 January 2012 identified by searching PubMed and ISI Web of Knowledge. The search string used was: (cox2 OR cox-2 OR cyclooxygenase-2 OR "cyclooxygenase 2" OR ptgs2) AND (prognos* OR survival OR mortality OR recurrence OR relapse OR progression) AND ("bladder cancer"). Studies were considered eligible if: (i) they reported the effect measure (as HRs, survival curves or log-rank p-values) of COX2 protein expression on recurrence, progression or disease-specific survival; (ii) COX2 was assessed in primary tumors exhibiting homogeneity in tumor histology (≥75% UCB), and subphenotype (≥75% NMIBC or MIBC); (iii) they were written in English or Spanish (Table S1). Reviews, abstracts, non-clinical studies, and duplicate publications were excluded. HRs and 95%CIs were directly extracted from the publications whenever available. For those reporting only the log-rank p-value or the Kaplan–Meier survival curves, the HRs and 95%CIs were independently calculated by two of the co-authors (MJC, AFSA) using the spreadsheet prepared by Sydes and Tierney with any discrepancies resolved by discussion [38]. In a few indicated cases, authors were directly contacted for clarification or provision of data not shown in the published manuscripts (Table S1). The level of heterogeneity among studies was calculated by means of the I² statistic [39], and publication bias was assessed by analyzing funnel plots and Egger’s asymmetry test [40].

Statistical Analysis

Associations between demographic and clinico-pathological parameters and COX2 expression were assessed using Fisher’s exact test. In NMIBCs, expression was also assessed distinctly in low-grade/risk (pTaG1/G2) and high-grade/risk (pTa/pT2G3) tumors, based on our previous evidence suggesting differential prognostic, genetic and molecular profiles between these subgroups [41], [42]. Recurrence-free, progression-free, and overall disease-specific survival curves were generated using the Kaplan-Meier method, with statistical significance assessed using the log-rank test. Time to each endpoint was calculated from date of primary treatment to the date of event, date of last follow-up, or date of patient’s death. Individuals who did not present any event until the end of the study, those lost to follow-up, or those who died from other causes were censored either at the time of last medical visit or at death. Time to recurrence and progression were defined by applying the “mid-time” between the date of the previous disease-free visit and that when a new event was diagnosed. Survival time was measured as the time from initial treatment to death resulting from cancer. Univariable and multivariable Cox-proportional hazards analysis was used to calculate hazard ratios (HR) and 95% confidence intervals (CI). Schoenfeld residual analysis did not suggest any departure from the proportional hazards assumption in multivariable models.

All statistical analyses were done using STATA (version 10.1 SE, StataCorp, College Station, TX, USA). Statistical tests were two-sided and p-values less than 0.05 were considered significant. The REMARK [43] guidelines for prognostic studies as well as the PRISMA [44] guidelines for systematic reviews and meta-analyses were adhered to in the preparation of the manuscript.

Results

Patients and COX2 Expression in Bladder Cancer TMAs

COX2 expression was assessed in 557 patients with NMIBCs and 216 individuals with MIBCs. Median COX2 expression was 0.121 au (range 0–42.590; interquartile range 1.382) in NMIBCs, and 0.760 au (0–30.806; 3.600) in MIBCs (p-value = 4×10⁻¹²). Representative COX2 immunostaining patterns in UCBs are shown in Figure S1. Of patients with NMIBCs, 41% (230/557) were treated only by transurethral resection, with the remainder (56%) receiving endovesical BCG immunotherapy and/or chemotherapy following transurethral resection, or other treatment (3%; Table 1). Nearly half (46%) of patients with MIBCs were treated by cystectomy, with the remainder receiving systemic chemotherapy, radiotherapy, superficial or other treatment, or some combination thereof (Table 2).

[Figure omitted. See PDF.]

Table 1. Distribution of characteristics of patients with NMIBCs by COX2 expression.

https://doi.org/10.1371/journal.pone.0045025.t001

[Figure omitted. See PDF.]

Table 2. Distribution of characteristics of patients with MIBCs by COX2 expression.

https://doi.org/10.1371/journal.pone.0045025.t002

COX2 Expression and Clinicopathological Features

Two-hundred eleven (38%) NMIBCs and 137 (58%) MIBCs expressed COX2 (Tables 1 and 2, respectively), with positive expression defined as a score equal to or greater than the ROC-derived threshold of 0.340 au. Patient and tumor characteristics in the analyzed sample did not differ significantly from the initial SBC/EPICURO study population with the exception of geographic region and tumor size in NMIBC patients (data not shown). The distribution of COX2 positivity was assessed according to established bladder cancer prognosticators including tumor invasion and grade, tumor multiplicity, tumor size and treatment, among others. Demographic factors like age, gender and region were not associated with COX2 expression, nor was the type of primary treatment received by patients (Tables 1 and 2). In NMIBCs, COX2 expression was significantly associated only with T and G; being more prominent in low-grade/risk pTaG1/2 tumors than in high-grade/risk pTa/pT1G3 tumors (p-value<0.0001; Table 1). Further assessment of COX2 distribution in relevant molecular subtypes of UCB [4], revealed a greater proportion of pTaG2 than pTaG1 tumors positively expressing COX2 in low-grade NMIBCs (p<0.0001, subtype 1; Figure S2). COX2 expression did not differ among high-grade/risk NMIBCs (p = 0.075), but a greater proportion of MIBCs positively expressed COX2 than did all high-grade/risk NMIBCs combined (p<0.0001, subtype 2; Figure S2). Only tumor multiplicity was associated with positive COX2 expression in MIBC patients (p-value = 0.008; Table 2).

COX2 Expression and Prognosis in Bladder Cancer Patients

We analyzed the association of COX2 expression with tumor recurrence and progression in patients with NMIBCs and with progression and disease-specific survival in patients with MIBCs (Table 3; Figure 1). When considered as a continuous variable in the univariable analysis, COX2 expression was marginally associated with an increased risk of recurrence in NMIBCs (HR = 1.02, 95%CI = 1.00–1.04, p-value = 0.048; Table 3). However, this association disappeared upon multivariable analysis when adjusting for region, gender, tumor stage and grade, multiplicity, tumor size, and treatment. Moreover, COX2 expression was not significantly associated with recurrence in NMIBCs when considered as a categorical variable, neither in the univariable nor multivariable analyses (Figure 1A; Figure S3AC; Table 3). Lastly, no significant association between COX2 expression and progression or survival was observed in patients with NMIBCs or MIBCs, regardless of whether expression was considered as a continuous or categorical variable in non-adjusted or adjusted analyses (Figure 1B–D; Figure S3B, 3D–H; Table 3).

[Figure omitted. See PDF.]

Table 3. Analysis of COX2 expression in NMIBCs and MIBCs; univariable and multivariable analyses.

https://doi.org/10.1371/journal.pone.0045025.t003

[Figure omitted. See PDF.]

Figure 1. Kaplan-Meier survival curves corresponding to failures in superficial (A, B) and invasive (C, D) tumors for specified prognostic endpoints.

Dashed curves: patients with tumors positive for COX2 protein staining; solid curves: patients with tumors negative for COX2 protein staining. Significance values from two-sided logrank test.

https://doi.org/10.1371/journal.pone.0045025.g001

Meta-analysis of COX2 Expression and Bladder Cancer Prognosis

Twenty publications on COX2 expression and bladder cancer prognosis were identified through the literature review (Table S1) [12]–[31]. Three of them lacked prognostic data, two overlapped with other larger studies and four included patient cohorts that did not meet the eligibility criteria outlined earlier, leaving 11 evaluable publications [12]–[14], [19], [21]–[25], [28], [29] plus the current study for the meta-analysis (Figure S4). Studies were classified by the tumor subtype(s) they reported on (i.e. NMIBC or MIBC), and whether adjustment for covariates was considered for each prognostic endpoint examined (i.e. univariable or multivariable; Figures 2 and 3). Of the four meta-analyses conducted with univariable data, only the metaHR of the association between COX2 expression and recurrence in NMIBCs showed marginal significance (metaHR = 1.35, 95%CI = 1.00–1.83; Figure 2). This result was not affected by study heterogeneity (I² p-value = 0.13) but exhibited significant publication bias, as evidenced by Egger’s test (p-value = 0.019). The remaining meta-analyses considering univariable data suggested increased, albeit non-significant, risks of tumor progression in patients with NMIBCs (metaHR = 2.07, 95%CI = 0.76–5.64) and MIBCs (metaHR = 1.45, 95%CI = 0.77–2.74), and death in patients with MIBCs (metaHR = 1.13, 95%CI = 0.8–1.59; Figure 2). Notably, the summary effect for progression in NMIBCs and that observed for survival in MIBCs were both significantly affected by study heterogeneity (I² p-values: 0.006 and 0.004, respectively), with the former also significantly influenced by publication bias (Egger’s test p-value = 0.001).

[Figure omitted. See PDF.]

Figure 2. Forest plots from selected univariable studies indicating the risk of reaching the indicated prognostic endpoints in non-muscle invasive (NMIBC; two upper panels), and muscle invasive (MIBC; two lower panels) UCBs in the presence of urothelial COX2 expression.

https://doi.org/10.1371/journal.pone.0045025.g002

Due to a paucity of published prognostic studies performing multivariable analysis on patients with NMIBCs, we could only address the multivariable meta-association with progression and survival in patients with MIBCs (Figure 3). A small, non-significant increased summary risk of progression (metaHR = 1.12, 95%CI = 0.53–2.35; Figure 3) was observed in COX2 expressing MIBCs that was unaffected by study heterogeneity (I² p-value = 0.139). Similarly, a null summary effect was observed for survival (metaHR = 0.97, 95%CI = 0.69–1.36; Figure 3). This effect was influenced neither by study heterogeneity (I² p-value = 0.114) nor by publication bias (Egger’s test p-value = 0.108.

[Figure omitted. See PDF.]

Figure 3. Forest plots from selected multivariable studies indicating the risk of reaching the indicated prognostic endpoints in non-muscle invasive (NMIBC; two upper panels), and muscle invasive (MIBC; two lower panels) UCBs in the presence of urothelial COX2 expression.

https://doi.org/10.1371/journal.pone.0045025.g003

Discussion

Despite many published studies, contradictory findings prevail on COX2 expression as an independent prognostic marker in patients with UCB. The current study suggests that COX2 expression is not an independent marker associated with recurrence, progression or survival in patients with UCB.

Using the largest cohort of patients with NMIBCs evaluated for COX2 expression to date, we observed that 38% of these tumors expressed the protein. Other groups have reported frequencies ranging from 53–88%; however, these studies used different COX2 antibodies and expression evaluation techniques and had smaller sample sizes [16], [28], [29], [45], [46]. In accordance with reported results [11], [18], [45] we observed significantly higher COX2 expression in MIBCs (58%) than in NMIBCs. This frequency is similar to that observed in other large, histologically homogeneous studies [21], [28], while groups using heterogeneous cohorts of squamous and transitional cell carcinomas report frequencies different from our own [12], [29]. Collectively, these findings reiterate the importance of homogeneity, or stratification, in tumor marker studies.

The association between COX2 and clinico-pathological characteristics remains a contentious issue in the literature. The majority of studies report an association between COX2 overexpression and advanced tumor invasion and grade, but use heterogeneous populations of NMIBCs and MIBCs in their assessments [21], [25], [26], [28], [47]. Given the known disparity in COX2 expression between NMIBCs and MIBCs, an association of this type would be expected in a mixed tumor population. After pooling NMIBCs and MIBCs in our study we also observe a strong significant association between COX2 overexpression and advanced tumor invasion (p>0.0001) and grade (p>0.0001). Notably, several groups report no association between COX2 expression and T and G [14], [29], [45]; especially those working strictly with homogeneous cohorts of MIBCs [12], [23]. Similarly, in our study, COX2 expression did not differ significantly among pT2, pT3 and pT4 tumors (p = 0.896). Interestingly, we observed lower COX2 positivity in pT1 and high-grade/risk NMIBCs, than in pTa and low-grade/risk NMIBCs tumors. This result may seem counterintuitive if grade progression is considered a linear trait and COX2 expression is deemed to increase linearly with T and G. However, there is strong evidence indicating that UCB exists as two molecularly distinct subtypes, with high-grade/risk NMIBCs having a molecular signature more similar to MIBCs than to low-grade/risk NMIBCs [4], [42]. In this respect, we observed that COX2 positivity increased significantly with increasing T and G within each molecular tumor subtype (Figure S2). Shirahama et al. [26] reported a COX2 distribution similar to ours, observing 8% positivity in pT1 tumors and 50% in MIBCs when using whole section staining and a 5% expression threshold. Collectively, these results reiterate the disparity in COX2 expression between NMIBCs and MIBCs first reported by Komhoff et al. [47], and highlight the importance of considering expression within the proper molecular context.

To minimize the effects resulting from selecting an arbitrary expression threshold, we investigated COX2 protein expression as a continuous variable and three categorical variables. Only when considered as a continuous variable in the univariable analysis was COX2 expression found to be associated with a slight increase in the risk of recurrence. The meta-analysis, consisting of five other univariable studies, reiterated this association and showed a 35% increased risk of recurrence in patients with COX2 expressing NMIBCs. However, both effect estimates exhibit only marginal significance, suggesting that the observed associations may be due to chance. Moreover, the association observed in the univariable analysis did not hold after adjustment for conventional prognostic factors of recurrence in the multivariable analysis. Lastly, the summary effect observed in the meta-analysis may have been skewed by two small studies which selected only high risk NMIBCs (T1G3 [19] and Cis [24]). When a sensitivity analysis was performed removing these two studies from the meta-analysis, the association between recurrence and COX2 expression was no longer maintained (metaHR = 1.14, 95%CI = 0.94–1.38). The observed disparity between effect estimates of progression in the present study and the meta-analysis could also be attributed to the inclusion of these two studies. Upon their exclusion, the summary HR showed no association with progression (metaHR = 0.98, 95%CI = 0.47–2.03). These results do not support a role for COX2 expression in NMIBCs as an independent prognostic marker of recurrence or progression.

Several groups have investigated the ability of COX2 expression to predict outcome in patients with MIBCs. Despite wide inter-study variation in methodology, antibodies used, sample size, and adjustment parameters in the case of multivariable analyses, the majority of these studies did not identify any significant association between COX2 expression and progression or survival, consistent with our findings [13], [14], [23], [28], [29]. Shariat and Margulis and their colleagues observed a negative association between high COX2 expression and tumor progression and mortality [21], [25]. However, both studies relied on heterogeneous sample populations which included a small proportion of patients with NMIBCs; potentially accounting for the observed associations given the disparity in COX2 expression between superficial and advanced bladder tumors [45]. In another study, Wulfing et al. reported that high COX2 expression was an independent predictor of poor overall survival in a subgroup of 62 patients with MIBC treated with cisplatin-based chemotherapy [29]. We did not identify any meaningful interaction between COX2 expression and treatment (data not shown), and were unable to replicate their findings in a smaller subset of 39 patients treated with cisplatin (HR = 1.47, 95%CI = 0.48–4.51, p-value = 0.497). Aziz et al. reported a 36% survival advantage associated with increased COX2 levels in a cohort of 266 patients with MIBCs (221 with UCB) that was independent of lymph node status and neo/adjuvant chemotherapy [12]. While we also observed improved survival among patients with COX2 overexpressing MIBCs, this association did not reach significance, consistent with other univariable [23], [28] and multivariable [26] analyses.

Our study had a large sample size, included only incident cases and relied on extensive and accurately acquired follow-up information spanning ten years. Additionally, we used automated scoring of immunostained TMAs, a strategy providing a reproducible assessment of expression that correlated highly with the independent evaluation of a subset of samples by an independent pathologist. COX2 staining was done in one laboratory to avoid heterogeneity in immunohistochemical staining and scoring, and evaluated as a continuous variable in the prognostic analyses to avoid potential bias related to selection of an expression threshold. Moreover, the sample population provides an accurate representation of bladder cancer in the general population as no inclusion criteria were applied in the recruitment process which included a good mix of referral centers and county hospitals. Lastly, the recommendations of the REMARK and PRISMA studies were followed in all of the reported analyses.

Despite these considerations and attempts to accurately quantify COX2 expression only in epithelial cells, the pathologist-trained automated imaging system may have incorporated some immunostained stromal material found on the tissue core, thereby increasing type I error. To reduce potential error we averaged the expression scores from duplicate cores and also explored a method investigated by Henriksen et al. [48] in which the higher score was used (data not shown). Both methods produced similar material associations between COX2 expression and clinico-pathological parameters or HRs. Moreover, adjusted analyses for progression in NMIBCs should be interpreted cautiously given the low number of events in relation to covariates. Also, different patient management practices across recruitment hospitals could increase sample heterogeneity, necessitating the inclusion of both recruitment area and treatment regimen in our multivariable analyses.

The results presented herein focus on COX2 expression levels measured in tumor epithelial cells – only one aspect of the complex interplay between the tumor and the host immune/inflammatory response [49]. The prognostic potential of COX2 (if any) may only be revealed when considered together with other tumoral markers. When investigating several potential prognostic parameters in UCB, Hilmy et al. concluded that systemic factors of the inflammatory response such as levels of C-reactive protein were superior to tumor-based factors such as grade, COX2 expression or T-lymphocytic infiltration [18]. Moreover, in models of cervical cancer, Ferrandina et al. observed that while COX2 expression was mutually exclusive in the tumor and stromal inflammatory cells, high expression in both cell types could be used as an independent marker of poor survival [50]. Future studies investigating the prognostic value of COX2 expression in UCB should take into consideration the multi-factorial and multi-dimensional context of the inflammatory response during carcinogenesis.

The current study is the largest to investigate COX2 expression as an independent marker of outcome in a prospective cohort of UCB patients. These findings, supported by a meta-analysis that included our own data and that from other relevant studies, do not support COX2 tumor cell expression being an independent prognosticator of UCB.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Immunohistochemical staining of COX2 in primary UCBs on TMAs. Expression was scored as a product of the percentage of epithelial area stained and the staining intensity using automated imaging analysis. A score of <0.340 au was considered negative for COX2 expression, while a score of ≥0.340 was considered positive. Representative sections of a pTaG1 UCB lacking COX2 expression (A and D) and a pT2G3 UCB expressing COX2 (B and E) are shown. Normal colon tissue was used as a positive control (C and F). Upper panels show sections under 100x magnification (A-C); lower panels show sections under 200x magnification (D–F).

https://doi.org/10.1371/journal.pone.0045025.s001

(PDF)

Figure S2.

Distribution of positive COX2 expression in urothelial carcionomas of the bladder classified by their molecular and pathological stage-grade subtypes. Positive COX2 expression assessed as described in Figure S1. Statistical significance assessed using Fisher’s exact test with a 0.05 significance level. pT1G2 tumors excluded due to low sample size in the current study (n = 11), and a reported tendency to overlap both molecular subtypes.

https://doi.org/10.1371/journal.pone.0045025.s002

(PDF)

Figure S3.

Kaplan-Meier survival curves corresponding to failures in superficial (A, B, C, D) and invasive (E, F, G, H) tumors for specified prognostic endpoints and quantiles of COX2 expression. Dashed curves: patients with tumors expressing COX2 at lower specified quantiles; solid curves: patients with tumors expressing COX2 at upper specified quantiles. Significance values from two-sided logrank test.

https://doi.org/10.1371/journal.pone.0045025.s003

(PDF)

Figure S4.

Flow diagram of study selection and inclusion in meta-analysis.

https://doi.org/10.1371/journal.pone.0045025.s004

(PDF)

Table S1.

Main characteristics of eligible studies used in meta-analysis.

https://doi.org/10.1371/journal.pone.0045025.s005

(PDF)

Acknowledgments

The authors acknowledge the corresponding authors of the studies indicated in Table S1 for providing additional data used in the meta-analysis; members of MJC’s thesis committee for scientific and methodological comments to the work; the Histology and Immunohistochemistry Core Unit at the CNIO for help with immunohistochemical staining of TMAs; and the patients, technicians, field workers and coordinators who made the study possible.

Spanish Bladder Cancer/EPICURO Study investigators

Institut Municipal d’Investigació Mèdica, Universitat Pompeu Fabra, Barcelona – Coordinating Center (M. Kogevinas, N. Malats, F.X. Real, M. Sala, G. Castaño, M. Torà, D. Puente, C. Villanueva, C. Murta-Nascimento, J. Fortuny, E. López, S. Hernández, R. Jaramillo, G. Vellalta, L. Palencia, F. Fermández, A. Amorós, A. Alfaro, G. Carretero); Hospital del Mar, Universitat Autònoma de Barcelona, Barcelona (J. Lloreta, S. Serrano, L. Ferrer, A. Gelabert, J. Carles, O. Bielsa, K. Villadiego), Hospital Germans Trias i Pujol, Badalona, Barcelona (L. Cecchini, J.M. Saladié, L. Ibarz); Hospital de Sant Boi, Sant Boi de Llobregat, Barcelona (M. Céspedes); Consorci Hospitalari Parc Taulí, Sabadell (C. Serra, D. García, J. Pujadas, R. Hernando, A. Cabezuelo, C. Abad, A. Prera, J. Prat); Centre Hospitalari i Cardiològic, Manresa, Barcelona (M. Domènech, J. Badal, J. Malet); Hospital Universitario de Canarias, La Laguna, Tenerife (R. García-Closas, J. Rodríguez de Vera, A.I. Martín); Hospital Universitario Nuestra Señora de la Candelaria, Tenerife (J. Taño, F. Cáceres); Hospital General Universitario de Elche, Universidad Miguel Hernández, Elche, Alicante (A. Carrato, F. García-López, M. Ull, A. Teruel, E. Andrada, A. Bustos, A. Castillejo, J.L. Soto); Universidad de Oviedo, Oviedo, Asturias (A. Tardón); Hospital San Agustín, Avilés, Asturias (J.L. Guate, J.M. Lanzas, J. Velasco); Hospital Central Covadonga, Oviedo, Asturias (J.M. Fernández, J.J. Rodríguez, A. Herrero), Hospital Central General, Oviedo, Asturias (R. Abascal, C. Manzano, T. Miralles); Hospital de Cabueñes, Gijón, Asturias (M. Rivas, M. Arguelles); Hospital de Jove, Gijón, Asturias (M. Díaz, J. Sánchez, O. González); Hospital de Cruz Roja, Gijón, Asturias (A. Mateos, V. Frade); Hospital Alvarez-Buylla (Mieres, Asturias): P. Muntañola, C. Pravia; Hospital Jarrio, Coaña, Asturias (A.M. Huescar, F. Huergo); Hospital Carmen y Severo Ochoa, Cangas, Asturias (J. Mosquera).

Author Contributions

Conceived and designed the experiments: NM FXR MK JL AC AT. Performed the experiments: MJC NM FXR MK JL AC AT SMM MMM. Analyzed the data: MJC AFSA. Contributed reagents/materials/analysis tools: NM FXR MK AC AT. Wrote the paper: MJC AFSA NM. Revised the manuscript: SMM MMM.

Citation: Czachorowski MJ, Amaral AFS, Montes-Moreno S, Lloreta J, Carrato A, Tardón A, et al. (2012) Cyclooxygenase-2 Expression in Bladder Cancer and Patient Prognosis: Results from a Large Clinical Cohort and Meta-Analysis. PLoS ONE 7(9): e45025. https://doi.org/10.1371/journal.pone.0045025

References

1. Jemal A, Bray F (2011) Center MM, Ferlay J, Ward E, et al (2011) Global cancer statistics. CA Cancer J Clin 61: 69–90.

2. Wu XR (2005) Urothelial tumorigenesis: a tale of divergent pathways. Nat Rev Cancer 5: 713–725.

3. Botteman MF, Pashos CL, Redaelli A, Laskin B, Hauser R (2003) The health economics of bladder cancer: a comprehensive review of the published literature. Pharmacoeconomics 21: 1315–1330.

4. Lindgren D, Frigyesi A, Gudjonsson S, Sjodahl G, Hallden C, et al. (2010) Combined gene expression and genomic profiling define two intrinsic molecular subtypes of urothelial carcinoma and gene signatures for molecular grading and outcome. Cancer Res 70: 3463–3472.

5. Michaud DS (2007) Chronic inflammation and bladder cancer. Urol Oncol 25: 260–268.

6. Fortuny J, Kogevinas M, Garcia-Closas M, Real FX, Tardon A, et al. (2006) Use of analgesics and nonsteroidal anti-inflammatory drugs, genetic predisposition, and bladder cancer risk in Spain. Cancer Epidemiol Biomarkers Prev 15: 1696–1702.

7. Daugherty SE, Pfeiffer RM, Sigurdson AJ, Hayes RB, Leitzmann M, et al. (2011) Nonsteroidal antiinflammatory drugs and bladder cancer: a pooled analysis. Am J Epidemiol 173: 721–730.

8. Harris RE (2007) Cyclooxygenase-2 (cox-2) and the inflammogenesis of cancer. Subcell Biochem 42: 93–126.

9. Greenhough A, Smartt HJ, Moore AE, Roberts HR, Williams AC, et al. (2009) The COX-2/PGE2 pathway: key roles in the hallmarks of cancer and adaptation to the tumour microenvironment. Carcinogenesis 30: 377–386.

10. Wheeler MA, Hausladen DA, Yoon JH, Weiss RM (2002) Prostaglandin E2 production and cyclooxygenase-2 induction in human urinary tract infections and bladder cancer. J Urol 168: 1568–1573.

11. Shirahama T (2000) Cyclooxygenase-2 expression is up-regulated in transitional cell carcinoma and its preneoplastic lesions in the human urinary bladder. Clin Cancer Res 6: 2424–2430.

12. Aziz A, Lessard A, Moore K, Hovington H, Latulippe E, et al. (2010) Improved cancer specific-survival in patients with carcinoma invading bladder muscle expressing cyclo-oxygenase-2. BJU Int 108: 531–537.

13. Bamias A, Kyriakou F, Chorti M, Kavantzas N, Noni A, et al. (2008) Microvessel density (MVD) and cyclooxygenase-2 (COX-2)/beta-catenin interaction are associated with relapse in patients with transitional carcinoma receiving adjuvant chemotherapy with paclitaxel/carboplatin: a hellenic cooperative oncology group (HECOG) study. Anticancer Res 28: 2479–2486.

14. Diamantopoulou K, Lazaris A, Mylona E, Zervas A, Stravodimos K, et al. (2005) Cyclooxygenase-2 protein expression in relation to apoptotic potential and its prognostic significance in bladder urothelial carcinoma. Anticancer Res 25: 4543–4549.

15. Eltze E, Wulfing C, Von Struensee D, Piechota H, Buerger H, et al. (2005) Cox-2 and Her2/neu co-expression in invasive bladder cancer. Int J Oncol 26: 1525–1531.

16. Friedrich MG, Toma MI, Petri S, Huland H (2003) Cyclooxygenase-2 promotes angiogenesis in pTa/T1 urothelial bladder carcinoma but does not predict recurrence. BJU Int 92: 389–392.

17. Gudjonsson S, Bendahl PO, Chebil G, Hoglund M, Lindgren D, et al. (2011) Can tissue microarray-based analysis of protein expression predict recurrence of stage Ta bladder cancer? Scand J Urol Nephrol 45: 270–277.

18. Hilmy M, Campbell R, Bartlett JM, McNicol AM, Underwood MA, et al. (2006) The relationship between the systemic inflammatory response, tumour proliferative activity, T-lymphocytic infiltration and COX-2 expression and survival in patients with transitional cell carcinoma of the urinary bladder. Br J Cancer 95: 1234–1238.

19. Kim SI, Kwon SM, Kim YS, Hong SJ (2002) Association of cyclooxygenase-2 expression with prognosis of stage T1 grade 3 bladder cancer. Urology 60: 816–821.

20. Liedberg F, Anderson H, Chebil G, Gudjonsson S, Hoglund M, et al. (2008) Tissue microarray based analysis of prognostic markers in invasive bladder cancer: much effort to no avail? Urol Oncol 26: 17–24.

21. Margulis V, Shariat SF, Ashfaq R, Thompson M, Sagalowsky AI, et al. (2007) Expression of cyclooxygenase-2 in normal urothelium, and superficial and advanced transitional cell carcinoma of bladder. J Urol 177: 1163–1168.

22. Mokos I, Jakic-Razumovic J, Marekovic Z, Pasini J (2006) Association of cyclooxygenase-2 immunoreactivity with tumor recurrence and disease progression in superficial urothelial bladder cancer. Tumori 92: 124–129.

23. Naruse K, Yamada Y, Nakamura K, Aoki S, Taki T, et al. (2010) Potential of molecular targeted therapy of HER-2 and Cox-2 for invasive transitional cell carcinoma of the urinary bladder. Oncol Rep 23: 1577–1583.

24. Shariat SF, Kim JH, Ayala GE, Kho K, Wheeler TM, et al. (2003) Cyclooxygenase-2 is highly expressed in carcinoma in situ and T1 transitional cell carcinoma of the bladder. J Urol 169: 938–942.

25. Shariat SF, Matsumoto K, Kim J, Ayala GE, Zhou JH, et al. (2003) Correlation of cyclooxygenase-2 expression with molecular markers, pathological features and clinical outcome of transitional cell carcinoma of the bladder. J Urol 170: 985–989.

26. Shirahama T, Arima J, Akiba S, Sakakura C (2001) Relation between cyclooxygenase-2 expression and tumor invasiveness and patient survival in transitional cell carcinoma of the urinary bladder. Cancer 92: 188–193.

27. Tiguert R, Lessard A, So A, Fradet Y (2002) Prognostic markers in muscle invasive bladder cancer. World J Urol 20: 190–195.

28. Wild PJ, Kunz-Schughart LA, Stoehr R, Burger M, Blaszyk H, et al. (2005) High-throughput tissue microarray analysis of COX2 expression in urinary bladder cancer. Int J Oncol 27: 385–391.

29. Wulfing C, Eltze E, von Struensee D, Wulfing P, Hertle L, et al. (2004) Cyclooxygenase-2 expression in bladder cancer: correlation with poor outcome after chemotherapy. Eur Urol 45: 46–52.

30. Yoshimura R, Sano H, Mitsuhashi M, Kohno M, Chargui J, et al. (2001) Expression of cyclooxygenase-2 in patients with bladder carcinoma. J Urol 165: 1468–1472.

31. Youssef RF, Shariat SF, Kapur P, Kabbani W, Mosbah A, et al. (2011) Prognostic Value of Cyclooxygenase-2 Expression in Squamous Cell Carcinoma of the Bladder. J Urol 185: 6.

32. Garcia-Closas M, Malats N, Silverman D, Dosemeci M, Kogevinas M, et al. (2005) NAT2 slow acetylation, GSTM1 null genotype, and risk of bladder cancer: results from the Spanish Bladder Cancer Study and meta-analyses. Lancet 366: 649–659.

33. Epstein JI, Amin MB, Reuter VR, Mostofi FK (1998) The World Health Organization/International Society of Urological Pathology consensus classification of urothelial (transitional cell) neoplasms of the urinary bladder. Bladder Consensus Conference Committee. Am J Surg Pathol 22: 1435–1448.

34. Mostofi F, Davis C, Sesterhen I (1999) Histological typing of urinary bladder tumours. World Health Organization international classification of histological tumours. Berlin: Springer Verlag.

35. Eble JN, Sauter G, Epstein JI, Sesterhen I (2004) Pathology and genetics of tumours of the urinary system and male genital organs. WHO classification of tumours. Lyon, France IARC Press.

36. Wahlin BE, Aggarwal M, Montes-Moreno S, Gonzalez LF, Roncador G, et al. (2010) A unifying microenvironment model in follicular lymphoma: outcome is predicted by programmed death-1–positive, regulatory, cytotoxic, and helper T cells and macrophages. Clin Cancer Res 16: 637–650.

37. Metz CE (1978) Basic principles of ROC analysis. Semin Nucl Med 8: 283–298.

38. Tierney JF, Stewart LA, Ghersi D, Burdett S, Sydes MR (2007) Practical methods for incorporating summary time-to-event data into meta-analysis. Trials 8: 16.

39. Higgins JP, Thompson SG (2002) Quantifying heterogeneity in a meta-analysis. Stat Med 21: 1539–1558.

40. Egger M, Davey Smith G, Schneider M, Minder C (1997) Bias in meta-analysis detected by a simple, graphical test. BMJ 315: 629–634.

41. Hernandez S, Lopez-Knowles E, Lloreta J, Kogevinas M, Amoros A, et al. (2006) Prospective study of FGFR3 mutations as a prognostic factor in nonmuscle invasive urothelial bladder carcinomas. J Clin Oncol 24: 3664–3671.

42. Lopez-Knowles E, Hernandez S, Kogevinas M, Lloreta J, Amoros A, et al. (2006) The p53 pathway and outcome among patients with T1G3 bladder tumors. Clin Cancer Res 12: 6029–6036.

43. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, et al. (2005) REporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer 93: 387–391.

44. Moher D, Liberati A, Tetzlaff J, Altman DG, The PG (2009) Preferred reporting items for systematic reviews and meta-analyses: the PRISMA Statement. Open Med 3: e123–e130.

45. Ristimaki A, Nieminen O, Saukkonen K, Hotakainen K, Nordling S, et al. (2001) Expression of cyclooxygenase-2 in human transitional cell carcinoma of the urinary bladder. Am J Pathol 158: 849–853.

46. Mohammed SI, Knapp DW, Bostwick DG, Foster RS, Khan KN, et al. (1999) Expression of cyclooxygenase-2 (COX-2) in human invasive transitional cell carcinoma (TCC) of the urinary bladder. Cancer Res 59: 5647–5650.

47. Komhoff M, Guan Y, Shappell HW, Davis L, Jack G, et al. (2000) Enhanced expression of cyclooxygenase-2 in high grade human transitional cell bladder carcinomas. Am J Pathol 157: 29–35.

48. Henriksen KL, Rasmussen BB, Lykkesfeldt AE, Moller S, Ejlertsen B, et al. (2007) Semi-quantitative scoring of potentially predictive markers for endocrine treatment of breast cancer: a comparison between whole sections and tissue microarrays. J Clin Pathol 60: 397–404.

49. Mantovani A, Allavena P, Sica A, Balkwill F (2008) Cancer-related inflammation. Nature 454: 436–444.

50. Ferrandina G, Lauriola L, Zannoni GF, Distefano MG, Legge F, et al. (2002) Expression of cyclooxygenase-2 (COX-2) in tumour and stroma compartments in cervical cancer: clinical implications. Br J Cancer 87: 1145–1152.