Oncotarget

Research Papers:

A cancer specific hypermethylation signature of the TERT promoter predicts biochemical relapse in prostate cancer: a retrospective cohort study

PDF |  HTML  |  Supplementary Files  |  How to cite

Oncotarget. 2016; 7:57726-57736. https://doi.org/10.18632/oncotarget.10639

Metrics: PDF 2674 views  |   HTML 4697 views  |   ?  

Pedro Castelo-Branco, Ricardo Leão, Tatiana Lipman, Brittany Campbell, Donghyun Lee, Aryeh Price, Cindy Zhang, Abolfazl Heidari, Derek Stephens, Stefan Boerno, Hugo Coelho, Ana Gomes, Celia Domingos, Joana D. Apolonio, Georg Schäfer, Robert G. Bristow, Michal R. Schweiger, Robert Hamilton, Alexandre Zlotta, Arnaldo Figueiredo, Helmut Klocker, Holger Sültmann and Uri Tabori _

Abstract

Pedro Castelo-Branco1,2,3,*, Ricardo Leão1,4,5,*, Tatiana Lipman1, Brittany Campbell1, Donghyun Lee1, Aryeh Price1, Cindy Zhang1, Abolfazl Heidari1, Derek Stephens1, Stefan Boerno6, Hugo Coelho5, Ana Gomes5, Celia Domingos2,3, Joana D. Apolonio2,3, Georg Schäfer11, Robert G. Bristow7, Michal R. Schweiger8,9, Robert Hamilton4, Alexandre Zlotta4,10, Arnaldo Figueiredo5, Helmut Klocker11, Holger Sültmann12, Uri Tabori1

1Arthur and Sonia Labatt Brain Tumor Research Center, The Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada

2Regenerative Medicine Program, Department of Biomedical Sciences and Medicine, University of Algarve, Faro, Portugal

3Centre for Biomedical Research (CBMR), University of Algarve, Faro, Portugal

4Division of Urology, Department of Surgical Oncology Princess Margaret Cancer Center, University of Toronto, Toronto, Ontario, Canada

5Serviço de Urologia e Transplantação Renal, Centro Hospitalar Universitário Coimbra EPE, Faculty of Medicine, University of Coimbra, Coimbra, Portugal

6Sequencing Core Facility, Max Planck Institute for Molecular Genetics, Berlin, Germany

7Department of Radiation Oncology, Princess Margaret Cancer Center, Toronto, Ontario, Canada

8Department of Vertebrate Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany

9Cologne Center for Genomics, Cologne University, Cologne, Germany

10Division of Urology, Department of Surgery, Mount Sinai Hospital, Toronto, Ontario, Canada

11Department of Urology, Medical University of Innsbruck, Innsbruck, Austria

12Cancer Genome Research, German Cancer Research Center (DKFZ) and German Consortium for Translational Cancer Research (DKTK), Heidelberg, Germany

*These authors contributed equally to this work

Correspondence to:

Uri Tabori, email: [email protected]

Pedro Castelo-Branco, email: [email protected]

Keywords: TERT, prostate cancer, biomarker, diagnostic, Gleason score

Received: March 31, 2016     Accepted: June 30, 2016     Published: July 16, 2016

ABSTRACT

The identification of new biomarkers to differentiate between indolent and aggressive prostate tumors is an important unmet need. We examined the role of THOR (TERT Hypermethylated Oncological Region) as a diagnostic and prognostic biomarker in prostate cancer (PCa).

We analyzed THOR in common cancers using genome-wide methylation arrays. Methylation status of the whole TERT gene in benign and malignant prostate samples was determined by MeDIP-Seq. The prognostic role of THOR in PCa was assessed by pyrosequencing on discovery and validation cohorts from patients who underwent radical prostatectomy with long-term follow-up data.

Most cancers (n = 3056) including PCa (n = 300) exhibited hypermethylation of THOR. THOR was the only region within the TERT gene that is differentially methylated between normal and malignant prostate tissue (p < 0.0001). Also, THOR was significantly hypermethylated in PCa when compared to paired benign tissues (n = 164, p < 0.0001). THOR hypermethylation correlated with Gleason scores and was associated with tumor invasiveness (p = 0.0147). Five years biochemical progression free survival (BPFS) for PCa patients in the discovery cohort was 87% (95% CI 73–100) and 65% (95% CI 52–78) for THOR non-hypermethylated and hypermethylated cancers respectively (p = 0.01). Similar differences in BPFS were noted in the validation cohort (p = 0.03). Importantly, THOR was able to predict outcome in the challenging (Gleason 6 and 7 (3 + 4)) PCa (p = 0.007). For this group, THOR was an independent risk factor for BPFS with a hazard-ratio of 3.685 (p = 0.0247). Finally, THOR hypermethylation more than doubled the risk of recurrence across all PSA levels (OR 2.5, p = 0.02).


INTRODUCTION

Prostate cancer (PCa) is the most frequently diagnosed cancer and the second most common cause of cancer-related mortality among men [1]. Although one sixth of men will be diagnosed with PCa during their lifetime, only one in thirty six will die from this disease [2].

Prostate cancer is a heterogeneous disease with risk that varies according to host and tumor characteristics that have not been fully elucidated [3, 4]. In view of this heterogeneous behaviour, the clinical challenge resides in maximizing patient survival without overtreatment of indolent tumors. Despite intense research, there is a lack of validated biomarkers that can help determine the natural history of PCa. This is especially true in low risk tumors (Gleason 6) where many patients will not experience tumor progression or death from disease and some intermediate risk (Gleason 7) where tumour behavior is particularly heterogeneous [2, 4, 5].

Some patients with Gleason 6 cancers will undergo active surveillance, while the majority of patients with Gleason 7 will undergo prostatectomy or radiation therapy with curative intent. Identifying which tumors will progress to advanced disease represents a major challenge in these PCa subgroups.

Several biomarkers and algorithms which utilize multiple parameters were recently described to predict PCa behavior [611]. However, many of these expression and methylation signatures are too complex to provide a clinically simple and robust tool which will be useful for the treating physician. An established oncogenic process is common in most recurrent cancers and is easy to detect without complex tools and would therefore be an attractive rational target to correlate PCa progression and patient outcome.

Cancer cells achieve limitless self-renewal capacity through the activation of telomere maintenance mainly through activation of telomerase which provides immortalization for 90% of cancers [12]. The catalytic subunit of the telomerase complex is termed Telomerase Reverse Transcriptase (TERT) and its expression has been observed in most malignant cancers including PCa [13]. Telomere shortening and telomere length has been shown to act as a predictor of disease progression in PCa [1416]. Moreover, telomerase activation stabilizes shorter telomeres and is a putative early marker for prostate carcinogenesis [1722]. Analysis of telomerase activity and TERT expression require high quality RNA and cell extracts that are challenging, especially when paraffin embedded tissues are considered. Therefore, a DNA based assay that correlates with telomerase activity would be extremely useful as a diagnostic and prognostic tool in cancer.

Epigenetic gene regulation through DNA methylation has been associated with diagnosis and prognosis in multiple cancers including brain and prostate cancer [911, 2326]. We recently identified a specific area in the TERT promoter, termed THOR (TERT Hypermethylated Oncological Region), which is hypermethylated only in cancers expressing TERT and non-hypermethylated in normal tissues and low-grade pediatric tumors, which do not express TERT. THOR predicted outcome and tumor progression in several subgroups of pediatric cancers [27].

Therefore, we postulated that THOR is hypermethylated in most cancers. In prostate cancer THOR demonstrated the diagnostic ability to differentiate cancer from normal prostate tissue and to discriminate indolent from aggressive PCa. Furthermore, THOR might be used as a marker to predict patient outcome in addition to other currently used markers.

RESULTS

THOR is highly methylated in common telomerase expressing cancers

In order to interrogate if THOR hypermethylation is observed in common adult cancers, we analyzed the CG site within THOR (CG11625005) in 11 cancer types from the Cancer Genome Atlas (n = 3056, Supplementary Table S1). THOR was hypermethylated in all cancers, which rely on telomerase activation for their telomere maintenance (Figure 1A). The beta value for prostate adenocarcinoma was 0.7, indicating a high degree of methylation at THOR, albeit with a large variance. Lower and more heterogeneous methylation statuses were observed in glioma and sarcoma. Interestingly, both cancers utilize both telomerase and the alternative lengthening of telomeres mechanism which do not exhibit THOR hypermethylation [2829]. Furthermore, the more indolent thyroid cancer exhibited the lowest levels of THOR hypermethylation.

A region in the hTERT promoter (THOR) is specifically hypermethylated in malignant prostate tissues.

Figure 1: A region in the hTERT promoter (THOR) is specifically hypermethylated in malignant prostate tissues. (A) Illumina Infinium 450 k array data obtained from the The Genome Cancer Atlas shows high THOR methylation status (cg11625005) in multiple tumors. (B) Methylation values of 51 tumours (dotted red lines) and their average (thick red line) as well as the methylation values of 53 normal prostate samples and their average (blue) are shown for 500 bp wide regions in the proximity of the TERT promoter through MeDIP-seq analysis. The methylation differences are most significant in region chr5: 1295501–1296000 (Benjamini Hochberg corrected Mann-Whitney p-value: 7.547178e-12) which matches the THOR. (C) Pyrosequencing analysis reveals that levels of THOR methylation are significantly higher in malignant prostate tissue when compared to its corresponding normals. (D) Levels of THOR methylation are significantly higher between any subgroup of Gleason scores and normal tissue and increase with Gleason scores, with statistical differences between Gleason 6 and Gleason ≥ 8 (p = 0.0416).

THOR distinguishes benign from malignant prostate tissue

In order to test whether the hypermethylation observed in the CG site within THOR in PCa is unique to this specific area of the TERT promoter we analyzed the methylation status of the whole TERT gene using MeDIP-seq on 51 PCa and 53 normal prostate tissues [30]. While methylation between normal prostate and PCa tissues is similar throughout the gene (Supplementary Figure S1), a significant difference in methylation between cancerous and normal tissue was observed only in the promoter region of TERT (p = 7.5 × 10−12), directly matching THOR (Figure 1B). To further explore THOR as a candidate cancer biomarker in PCa, we used pyrosequencing on multiple PCa samples from our discovery cohort. In 164 prostatectomies where matched benign and malignant tissues were available, THOR was significantly hypermethylated in the PCa component (p < 0.0001, Figure 1C).

THOR Hypermethylation correlates with other bio-pathological risk factors

To further elaborate on THOR’s association with known PCa risk factors, we initially compared THOR methylation with increasing Gleason scores. A significant difference between benign tissue and Gleason 6 PCa was observed (p < 0.0001). Analysis of THOR (as a continous variable) and Gleason revealed that THOR Hypermethylation is positively associated with Gleason score (X2 = 9.60; p = 0.0082). Although THOR methylation exhibited significant difference between low grade (Gleason 6) PCa and high grade (Gleason score ≥ 8) (p = 0.0416, Figure 1D) this was not observed between Gleason 6–7 tumors. We then determined the association between THOR and different established risk criteria. High-risk tumors (Gleason score ≥ 8 or PSA ≥ 20 ng/mL) had significantly higher THOR methylation than low risk PCa (Gleason score 6 and PSA < 10 ng/mL) in both the discovery cohort (p = 0.0436, Supplementary Figure S2A) and the validation cohort (p = 0.0397, Supplementary Figure S2B). Interestingly, higher THOR methylation was associated with locally advanced disease in both the discovery and validation cohorts (p = 0.0440 and p = 0.0147 respectively, Supplementary Figure S2C and S2D). In contrast, we didn’t observe any association between THOR methylation and PSA levels, age, prostate volume or TMPRSS2-ERG gene fusion [31] (Supplementary Table S2). Margin status was also evaluated as a potential prognostic factor for biochemical relapse. However, this parameter was not significantly associated with time to biochemical relapse in both cohorts (Discovery cohort p = 0.2; Validation cohort p = 0.9).

THOR as a novel risk stratification marker for low and intermediate Gleason PCa

To assess the role of THOR as a prognostic marker in PCa, we chose a threshold of 20% methylation as previously done (AUC of 0.799, p < 0.0001) [27]. The median follow-up for patients in the discovery cohort was 9,8+/−3,8 years and for patients in the validation cohort was 7,9+/−2,5 years. A total of 139 patients were analyzed for biochemical relapse and THOR methylation. Five-year biochemical progression free survival (BPFS) for the discovery cohort was 65% (95% CI 52–78) and 87% (95% CI 73–100) for THOR hypermethylated and non-hypermethylated PCa respectively (p = 0.015, Figure 2A). Similar superior BPFS was observed for non-hypermethylated PCas in the validation cohort (p = 0.0306, Figure 2B). Combining the 2 cohorts revealed similar superior BPFS for non-hypermethylated PCa (p = 0.01, Supplementary Figure S3A). To determine the ability of THOR to predict BPFS in the clinically low and intermediate PCa, we analyzed patients with lower Gleason scores. THOR non-hypermethylation was associated with improved BPFS in both Gleason 6 and 7 PCa (p = 0.016 and p = 0.02 respectively, Supplementary Figure S3B and S3C). Since Gleason 7 tumors are known to have diverse clinical outcomes and have shown heterogeneous THOR values, we divided our survival analysis to Gleason 3 + 4 and 4 + 3 subgroups. Strikingly, THOR status predicted BPFS in Gleason 7 (3 + 4) tumors but failed to do so in Gleason 7 (4 + 3) PCa (Supplementary Figure S3D, S3E). Furthermore, most Gleason 8 PCa exhibited THOR hypermethylation and THOR values did not predict survival for higher Gleason tumors (Supplementary Figure S3F). Combined, THOR methylation defines a novel risk group for PCa where in lower Gleason 6 and 7 (3 + 4) THOR predicts outcome (p = 0.007) while this signature is lost in higher Gleason scores (Figure 2C and 2D). Univariate and multivariate analysis revealed that for this risk group (Gleason 6 and Gleason 7 (3 + 4) THOR Hypermethylation has the highest risk for recurrence with HR of 6.224 (p = 0.0012; C-índex 0.91, Table 1) and is an independent risk factor for BPFS (HR: 3.684; Ci: 1.81–11.5; p = 0.0247).

Levels of THOR methylation stratify prostate cancer patients.

Figure 2: Levels of THOR methylation stratify prostate cancer patients. Biochemical recurrence reveals that patients from all Gleason scores (6, 7 and ≥ 8) with low levels of THOR methylation show significantly better biochemical progression free survival in both the discovery (A) and validation cohorts (B). Patients from both cohorts with Gleason 6 and intermediate risk (Gleason 7 (3 + 4)) and low THOR status present significantly better progression free survival when compared to patients with high levels of THOR methylation (C). Considering higher risk groups (Gleason 7 (4 + 3) and Gleason ≥ 8 patients)) the levels of THOR methylation did not show differences in terms of biochemical progression free survival (D).

Table 1: Univariate and multivariate analysis of time to biochemical recurrence among patients with Gleason 6 and Gleason 7 (3 + 4)

 

UNIVARIATE ANALYSIS

MULTIVARIATE ANALYSIS

HR

95% CI

Chi Square

P

C índex

HR

95% CI

Chi Square

P

Both Cohorts

Age

1.047

0.972 to 1.172

1.4585

0.222

0.62

1.036

0.960 to 1.17

0.083

0.3623

pT (localized vs advanced)

2.748

1.176 to 6.469

5.4423

0.0197

0.88

2.454

0.967 to 6.230

3.685

0.0589

PSA Groups (< 10 vs ≥ 10)

3.288

1.307 to 8.271

6.3995

0.0114

0.92

2.549

0.931 to 6.979

3.314

0.0687

THOR

6.224

1.314 to 9.751

6.2240

0.00126

0.91

3.685

1.181 to 11.501

5.046

0.0247

Combining PSA and THOR analysis significantly increases patients’ outcome prediction

Finally, we tested the ability of THOR to add information to the commonly used PSA as a predictor of outcome. For each PSA value, THOR hypermethylation more than doubled the risk of BPFS (OR 2.5, p = 0.02, Ci: 1.15–5.6; Figure 3). These findings were highly consistent between cohorts (Supplementary Figure S4A and S4B). All PCa with very high PSA (> 25 ng / mL) were hypermethylated and > 90% experienced tumor recurrence.

Estimated probability for biochemical relapse.

Figure 3: Estimated probability for biochemical relapse. Analysis of both cohorts together reveals that patients with high levels of THOR methylation have a significant increase in the probability of recurrence for the same values of PSA when compared to patients with low levels of THOR methylation.

DISCUSSION

In this study we extend our previous observations, which suggest a clinical role for the methylation signature of the TERT promoter in cancer. Specifically, we show that THOR hypermethylation can serve as a robust and simple tool for predicting tumor behavior in low grade PCa, where management is controversial. These findings can add important information for the clinical management of these patients.

Global DNA methylation status is highly variable during early embryogenesis but is quite stable in most tissues throughout life [3234]. Aberrant methylation of the cancer genome has major implications on gene expression and has been recently reported to refine tumor subgroups and be associated with patient survival [35, 36]. Since the assays performed are based on whole-genome arrays, changes observed are usually global and determination of methylation status is restricted to the whole gene or its promoter. Analysis of methylation regions within genes are rarely performed [37] and has not been developed as a tool for clinical diagnostics. Specifically, data on methylation status of regions within a promoter of an oncogene and its effect on gene expression in cancer are limited. Moreover, other alterations which result in telomere maintenance in cancer such as TERT promoter mutations and the presence of alternative lengthening of telomeres are rarely observed in PCa. [3839]. Our observations suggest that mapping of TERT methylation (Supplementary Figure S1) can add valuable information to the mechanisms of oncogene activation during carcinogenesis. We have previously shown that THOR methylation is a dynamic process during gliomagenesis [27]. Further studies are required to understand the causes and consequences of hypermethylation of this specific region on TERT activation and PCa progression.

As in most advanced carcinomas, aggressive PCa express telomerase, and interestingly, for tumors with high Gleason scores (≥ 8) THOR lost its prognostic value. Although this is mainly due to the vast majority of these cancers having high THOR methylation and the small amount of tumors with low methylation values to show statistical difference, this further suggests that in contrast to early stages PCa, where some tumors might lack self-renewal capacity, advanced stage tumors will maintain their telomeres by either THOR hypermethylation or other pathways to activate telomerase.

Nevertheless, currently approximately 50% of newly diagnosed PCa patients are found to have a low risk prostate cancer [45] and a significant proportion of those cancers may never become life threatening. Furthermore, the risk of biochemical recurrence after definitive surgery is highly variable and is usually poorly understood [4648].

For this challenging patient population, THOR analysis adds a new dimension for the decision making process:

First, among patients with Gleason 6 PCa, the risk of tumor recurrence is extremely low if THOR is non-hypermethylated (no progressions in 5 years, Supplementary Figure S3B) [49, 50]. This candidate biomarker may be highly valuable for patients undergoing active surveillance protocols, however, further confirmation from biopsy tissues will be required.

Second, the most relevant observation in our study is that THOR stratifies Gleason 7 tumors into 2 risk groups where the 3 + 4 group could be stratified with the Gleason 6 tumors while the 4 + 3 PCa are comparable with higher Gleason scores (Supplementary Figure S3) [51]. If confirmed, it will be possible in the future to lump the low (Gleason 6) and intermediate risk (Gleason 7 (3 + 4)) PCa to a single risk group for a more conservative management based on THOR status [52]. Indeed, for this combined group, THOR showed to be able to stratify patients into particularly low risk of recurrence (only 5% of patients with THOR non-hypermethylated PCa experienced biochemical recurrence within the first 5 years (Figure 2C)) while THOR hypermethylation increased the risk of recurrence by more than 6 fold (Table 1). This can be of clinical value for selecting patients to undergo adjuvant treatment after surgery and eventually applied to prostate biopsies establishing new criteria for active surveillance protocols.

Third, THOR status adds valuable information for each PSA level for these patients. THOR hypermethylated cancers have > 50% recurrence even with a PSA value of less than 10 ng/mL (Figure 3) while non-hypermethylated cancers carry a very low risk of recurrence in both cohorts.

Using telomerase as a biomarker has several advantages over some of potential biomarkers which have recently described for risk stratification of PCa [53, 54, 7]. Biomarkers which require RNA and use genomic information are complex, highly dependent on tissue quality and can only be performed in specialized laboratories. THOR is a DNA based marker and has proven in this and other studies to provide robust results in DNA from paraffin embedded samples and in degraded DNA. This assay could therefore be performed in most laboratories worldwide. Furthermore, telomerase represents a rational and attractive oncogene to pursue. Indeed, unlike other biomarkers, patients possessing this biomarker could have a potentially drug-targetable marker in that tumors with hypermethylated THOR could be treated using both telomerase inhibitors and demethylating enzymes.

Our studies have limitations related to retrospective cohorts. The use of biochemical relapse as an endpoint is suboptimal compared with prostate cancer-specific mortality or time to metastasis. Also, we acknowledge that our results are based on surgical specimens. However, analysis of these surgical specimens could identify within the low-intermediate risk group a subgroup of patients to whom early adjuvant treatment might be beneficial.

Furthermore, our data unveils the potential predictive value of THOR methylation when applied to tissue from prostate biopsies. THOR could stratify patients with low risk disease in both cohorts independently as the difference in outcome between cohorts did not change the general role of THOR as a candidate prognostic biomarker in PCa.

In the era of precision medicine where the aim is to determine patient-specific outcome, THOR hypermethylation represents a potential candidate biomarker for cancer diagnostics in biopsies. It can also determine aggressiveness of tumors with similar histological grade and be used as a companion biomarker for therapies using telomerase inhibitors or demethylating enzymes.

Finally, combining THOR with other biomarkers such as PSA can identify the patients with extremely low risk of recurrence where surveillance protocols could be applied. Further studies are required to verify if THOR methylation can change the current treatment paradigms in low or intermediate risk PCa.

MATERIALS AND METHODS

Patients

Tissue samples and patient data were obtained upon consent according to the Research Ethics Boards at the participating institutions. For the MeDIP-Seq analysis of TERT, tumors selected were staged pT2–pT4 and had Gleason scores ranging from 6 to 9 (Heidelberg, Germany) [30]. Two additional cohorts were definided to retrospectively study the clinical implications of THOR in PCa. These included a discovery cohort of patients from Austria (n = 164) and a separate validation cohort from Portugal (n = 103), submitted to radical prostatectomy and from which FFPE tissue was available. Patients were selected based on the availability of FFPE tissue, follow-up time and available clinical information. Gleason score classification was considered for patient selection in validation cohort (mimicking discovery cohort). Clinical outcomes were blinded at patient selection. Patients lost for follow-up or with missing values were not considered for outcome analysis. Clinical outcomes, PSA values and prostate biopsy results were not used as selection criteria.

Demographic analysis and clinical characteristics for the discovery and validation cohorts are described in Table 2. From each surgical specimen we analyzed the malignant and benign tissue separately. Pathological evaluation was performed by experienced uro-pathologists in both centres. Dominant lesions (higher Gleason score) were selected for this study, by macrodissection (where at least 70% of malignant cells were found). Benign tissue was isolated from the same surgical specimen where tumours were isolated.

Table 2: Demographic and clinical characteristics of discovery and validation cohorts

Biochemical Progression free survival (BPFS) was defined as the time from surgery to biochemical recurrence (PSA > 0, 2 ng/mL in two different measurements after an undetectable post-surgery PSA measurement). All patients submitted to any adjuvant treatment that could alter the normal evolution of the disease were excluded from outcome analysis.

Open access data

Dataset from The Cancer Genome Atlas was extracted for the Illumina Infinium 450 k array beta values. Mean and standard deviation were calculated from beta values for each tumour type. The single probe located in the THOR region (cg11625005) was used for this analysis.

MeDIP-Seq analysis

For in-depth analysis of the TERT gene methylation status MeDIP-Seq was used [30]. Library preparation 2.5 μg of genomic DNA from 51 prostate cancers or 53 normal prostate tissues were fragmented to 100 to 200 bp using the Covaris S2 system and end repaired with End Repair mix (Enzymatics) followed by a purification step (Qiagen DNA Purification Kit) and ligation of barcoded SOLiD sequencing adapters as previously described [30].

Analysis of THOR methylation

Quantitative sodium bisulfite pyrosequencing was performed for THOR as previously described [27]. In brief, targeted assays were designed using the PyroMark Assay Design Software 1.0 (Qiagen). Forward ATGATGTGGAGGTTTTGGGAATAG, reverse CCCAACCTAAAAACAACCCTAAAT and sequencing GGAGGTTTTGGGAATAG primers were used for PCR and pyrosequencing. The assay target region was 36 bp in length comprising 5 CpG sites. In our assay < 5% of the samples failed pyrosequencing analysis. Calculation of the % of THOR methylation was done as a mean value of these sites as previously described [27]. For clinical correlative studies we used the cut-off of 20% methylation with an AUC of 0.799 (p < 0.0001).

Statistical analysis

MeDIP-Seq statistical analysis was conducted using R (version 2.9.2). Read counts in 500 bp non-overlapping consecutive bins were normalised to sample wise read counts (reads per million). To assess the difference in THOR Hypermethylation between normal and malignant tissue a two-tailed Student´s t test was used. To test the association of THOR with Gleason Scores, localized prostate cancer disease and locally advanced disease Mann-Whitney U test was used. P Pearson test assessed relation between methylation and age and methylation and prostate volume. For the prognostic model we initially dichotomized into high and low-methylation for THOR groups by receiver operating characteristic (ROC) analysis. Biochemical progression free survival (BPFS) was determined by Kaplan-Meier Survival curves on 139 patients for whom all clinical information was available and were not excluded by our criteria. Patients not having experienced PSA recurrence were censored at their last PSA measurement. The proportional hazards assumption was verified by the log-negative-log survival distribution function for all variables. Univariate and multivariate Cox Proportional Hazards (CPH) regression analyses and log-rank tests were conducted. Both analyses (Univariate and Multivariate) were done for time for biochemical recurrence. To evaluate the prognostic strength of THOR, the C-index was used. To estimate the additional role of THOR methylation in predicting BPFS for patients with different values of PSA a logistic regression was defined. All statistical analyses were performed using SAS V9.3 (SAS, Cary, NC).

CONFLICTS OF INTEREST

The authors have no conflicts of interest to declare.

GRANT SUPPORT

Operating grants from The Canadian Institute of Health Research; Canadian Cancer Society Research Institute, MaRS Innovation, FCT Research Center Grant UID/BIM/04773/2013 CBMR 1334, FCT Doctoral Grant SFRH/BD/102232/2014 and Volkswagenstiftung (Lichtenberg Program to M.R.S).

REFERENCES

1. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015; 136:E359–386.

2. Wilt TJ, Brawer MK, Jones KM, Barry MJ, Aronson WJ, Fox S, Gingrich JR, Wei JT, Gilhooly P, Grob BM, Nsouli I, Iyer P, Cartagena R, et al. Radical prostatectomy versus observation for localized prostate cancer. N Engl J Med. 2012; 367:203–213.

3. Tosoian JJ, Mamawala M, Epstein JI, Landis P, Wolf S, Trock BJ, Carter HB. Intermediate and Longer-Term Outcomes From a Prospective Active-Surveillance Program for Favorable-Risk Prostate Cancer. J Clin Oncol. 2015.

4. Boutros PC, Fraser M, Harding NJ, de Borja R, Trudel D, Lalonde E, Meng A, Hennings-Yeomans PH, McPherson A, Sabelnykova VY, Zia A, Fox NS, Livingstone J, et al. Spatial genomic heterogeneity within localized, multifocal prostate cancer. Nat Genet. 2015; 47:736–745.

5. Sakr WA, Tefilli MV, Grignon DJ, Banerjee M, Dey J, Gheiler EL, Tiguert R, Powell IJ, Wood DP. Gleason score 7 prostate cancer: a heterogeneous entity? Correlation with pathologic parameters and disease-free survival. Urology. 2000; 56:730–734.

6. Hieronymus H, Schultz N, Gopalan A, Carver BS, Chang MT, Xiao Y, Heguy A, Huberman K, Bernstein M, Assel M, Murali R, Vickers A, Scardino PT, et al. Copy number alteration burden predicts prostate cancer relapse. Proc Natl Acad Sci U S A. 2014; 111:11139–11144.

7. Klein EA, Cooperberg MR, Magi-Galluzzi C, Simko JP, Falzarano SM, Maddala T, Chan JM, Li J, Cowan JE, Tsiatis AC, Cherbavaz DB, Pelham RJ, Tenggara-Hunter I, et al. A 17-gene assay to predict prostate cancer aggressiveness in the context of Gleason grade heterogeneity, tumor multifocality, and biopsy undersampling. Eur Urol. 2014; 66:550–560.

8. Lalonde E, Ishkanian AS, Sykes J, Fraser M, Ross-Adams H, Erho N, Dunning MJ, Halim S, Lamb AD, Moon NC, Zafarana G, Warren AY, Meng X, et al. Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study. Lancet Oncol. 2014; 15:1521–1532.

9. Litovkin K, Joniau S, Lerut E, Laenen A, Gevaert O, Spahn M, Kneitz B, Isebaert S, Haustermans K, Beullens M, Van Eynde A, Bollen M. Methylation of PITX2, HOXD3, RASSF1 and TDRD1 predicts biochemical recurrence in high-risk prostate cancer. J Cancer Res Clin Oncol. 2014; 140:1849–1861.

10. Kristensen H, Haldrup C, Strand S, Mundbjerg K, Mortensen MM, Thorsen K, Ostenfeld MS, Wild PJ, Arsov C, Goering W, Visakorpi T, Egevad L, Lindberg J, et al. Hypermethylation of the GABRE~miR-452~miR-224 promoter in prostate cancer predicts biochemical recurrence after radical prostatectomy. Clin Cancer Res. 2014; 20:2169–2181.

11. Haldrup C, Mundbjerg K, Vestergaard EM, Lamy P, Wild P, Schulz WA, Arsov C, Visakorpi T, Borre M, Hoyer S, Orntoft TF, Sorensen KD. DNA methylation signatures for prediction of biochemical recurrence after radical prostatectomy of clinically localized prostate cancer. J Clin Oncol. 2013; 31:3250–3258.

12. Newbold RF. The significance of telomerase activation and cellular immortalization in human cancer. Mutagenesis. 2002; 17:539–550.

13. Bodnar AG, Ouellette M, Frolkis M, Holt SE, Chiu CP, Morin GB, Harley CB, Shay JW, Lichtsteiner S, Wright WE. Extension of life-span by introduction of telomerase into normal human cells. Science. 1998; 279:349–352.

14. Heaphy CM, Yoon GS, Peskoe SB, Joshu CE, Lee TK, Giovannucci E, Mucci LA, Kenfield SA, Stampfer MJ, Hicks JL, De Marzo AM, Platz EA, Meeker AK. Prostate cancer cell telomere length variability and stromal cell telomere length as prognostic markers for metastasis and death. Cancer Discov. 2013; 3:1130–1141.

15. Meeker AK, Hicks JL, Platz EA, March GE, Bennett CJ, Delannoy MJ, De Marzo AM. Telomere shortening is an early somatic DNA alteration in human prostate tumorigenesis. Cancer Res. 2002; 62:6405–6409.

16. Hurwitz LM, Heaphy CM, Joshu CE, Isaacs WB, Konishi Y, De Marzo AM, Isaacs SD, Wiley KE, Platz EA, Meeker AK. Telomere length as a risk factor for hereditary prostate cancer. Prostate. 2014; 74:359–364.

17. Scates DK, Muir GH, Venitt S, Carmichael PL. Detection of telomerase activity in human prostate: a diagnostic marker for prostatic cancer? Br J Urol. 1997; 80:263–268.

18. Meeker AK. Telomeres and telomerase in prostatic intraepithelial neoplasia and prostate cancer biology. Urol Oncol. 2006; 24:122–130.

19. Wang Z, Ramin SA, Tsai C, Lui P, Herbert PJ, Kyeyune-Nyombi E, Ruckle HC, Beltz RE, Sands JF. Detection of telomerase activity in prostatic fluid specimens. Urol Oncol. 2000; 6:4–9.

20. Koeneman KS, Pan CX, Jin JK, Pyle JM 3rd, Flanigan RC, Shankey TV, Diaz MO. Telomerase activity, telomere length, and DNA ploidy in prostatic intraepithelial neoplasia (PIN). J Urol. 1998; 160:1533–1539.

21. Kallakury BV, Brien TP, Lowry CV, Muraca PJ, Fisher HA, Kaufman RP Jr, Ross JS. Telomerase activity in human benign prostate tissue and prostatic adenocarcinomas. Diagn Mol Pathol. 1997; 6:192–198.

22. Zhang W, Kapusta LR, Slingerland JM, Klotz LH. Telomerase activity in prostate cancer, prostatic intraepithelial neoplasia, and benign prostatic epithelium. Cancer Res. 1998; 58:619–621.

23. Liu L, Kron KJ, Pethe VV, Demetrashvili N, Nesbitt ME, Trachtenberg J, Ozcelik H, Fleshner NE, Briollais L, van der Kwast TH, Bapat B. Association of tissue promoter methylation levels of APC, TGFbeta2, HOXD3 and RASSF1A with prostate cancer progression. International Int J Cancer 2011; 129:2454–2462.

24. Cottrell S, Jung K, Kristiansen G, Eltze E, Semjonow A, Ittmann M, Hartmann A, Stamey T, Haefliger C, Weiss G. Discovery and validation of 3 novel DNA methylation markers of prostate cancer prognosis. J Urol. 2007; 177:1753–1758.

25. Olkhov-Mitsel E, Van der Kwast T, Kron KJ, Ozcelik H, Briollais L, Massey C, Recker F, Kwiatkowski M, Fleshner NE, Diamandis EP, Zlotta AR, Bapat B. Quantitative DNA methylation analysis of genes coding for kallikrein-related peptidases 6 and 10 as biomarkers for prostate cancer. Epigenetics. 2012; 7:1037–1045.

26. Kron K, Liu L, Trudel D, Pethe V, Trachtenberg J, Fleshner N, Bapat B, van der Kwast T. Correlation of ERG expression and DNA methylation biomarkers with adverse clinicopathologic features of prostate cancer. Clin Cancer Res. 2012; 18:2896–2904.

27. Castelo-Branco P, Choufani S, Mack S, Gallagher D, Zhang C, Lipman T, Zhukova N, Walker EJ, Martin D, Merino D, Wasserman JD, Elizabeth C, Alon N, et al. Methylation of the TERT promoter and risk stratification of childhood brain tumours: an integrative genomic and molecular study. Lancet Oncol. 2013; 14:534–542.

28. Dorris K, Sobo M, Onar-Thomas A, Panditharatna E, Stevenson CB, Gardner SL, Dewire MD, Pierson CR, Olshefski R, Rempel SA, Goldman S, Miles L, Fouladi M, et al. Prognostic significance of telomere maintenance mechanisms in pediatric high-grade gliomas. J Neurooncol. 2014; 117:67–76.

29. Liau JY, Tsai JH, Jeng YM, Lee JC, Hsu HH, Yang CY. Leiomyosarcoma With Alternative Lengthening of Telomeres Is Associated With Aggressive Histologic Features, Loss of ATRX Expression, and Poor Clinical Outcome. Am J Surg Pathol. 2015; 39:236–244.

30. Borno ST, Fischer A, Kerick M, Falth M, Laible M, Brase JC, Kuner R, Dahl A, Grimm C, Sayanjali B, Isau M, Rohr C, Wunderlich A, et al. Genome-wide DNA methylation events in TMPRSS2-ERG fusion-negative prostate cancers implicate an EZH2-dependent mechanism with miR-26a hypermethylation. Cancer Discov. 2012; 2:1024–1035.

31. FitzGerald LM, Agalliu I, Johnson K, Miller MA, Kwon EM, Hurtado-Coll A, Fazli L, Rajput AB, Gleave ME, Cox ME, Ostrander EA, Stanford JL, Huntsman DG. Association of TMPRSS2-ERG gene fusion with clinical characteristics and outcomes: results from a population-based study of prostate cancer. BMC Cancer. 2008; 8:230.

32. Brunner AL, Johnson DS, Kim SW, Valouev A, Reddy TE, Neff NF, Anton E, Medina C, Nguyen L, Chiao E, Oyolu CB, Schroth GP, Absher DM, et al. Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver. Genome Res. 2009; 19:1044–1056.

33. Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, Eggan K, Meissner A. DNA methylation dynamics of the human preimplantation embryo. Nature. 2014; 511:611–615.

34. Turecki G. The molecular bases of the suicidal brain. Nat Rev Neurosci. 2014; 15:802–816.

35. Dubuc AM, Remke M, Korshunov A, Northcott PA, Zhan SH, Mendez-Lago M, Kool M, Jones DT, Unterberger A, Morrissy AS, Shih D, Peacock J, Ramaswamy V, et al. Aberrant patterns of H3K4 and H3K27 histone lysine methylation occur across subgroups in medulloblastoma. Acta Neuropathol. 2013; 125:373–384.

36. White-Al Habeeb NM, Ho LT, Olkhov-Mitsel E, Kron K, Pethe V, Lehman M, Jovanovic L, Fleshner N, van der Kwast T, Nelson CC, Bapat B. Integrated analysis of epigenomic and genomic changes by DNA methylation dependent mechanisms provides potential novel biomarkers for prostate cancer. Oncotarget. 2014; 5:7858–7869. doi: 10.18632/oncotarget.2313.

37. Jones DT, Jager N, Kool M, Zichner T, Hutter B, Sultan M, Cho YJ, Pugh TJ, Hovestadt V, Stutz AM, Rausch T, Warnatz HJ, Ryzhova M, et al. Dissecting the genomic complexity underlying medulloblastoma. Nature. 2012; 488:100–105.

38. Stoehr R, Taubert H, Zinnall U, Giedl J, Gaisa NT, Burger M, Ruemmele P, Hurst CD, Knowles MA, Wullich B, Hartmann A. Frequency of TERT Promoter Mutations in Prostate Cancer. Pathobiology. 2015; 82:53–57.

39. Wu S, Huang P, Li C, Huang Y, Li X, Wang Y, Chen C, Lv Z, Tang A, Sun X, Lu J, Li W, Zhou J, et al. Telomerase reverse transcriptase gene promoter mutations help discern the origin of urogenital tumors: a genomic and molecular study. Eur Urol. 2014; 65:274–277.

40. Bougel S, Lhermitte B, Gallagher G, de Flaugergues JC, Janzer RC, Benhattar J. Methylation of the hTERT promoter: a novel cancer biomarker for leptomeningeal metastasis detection in cerebrospinal fluids. Clin Cancer Res. 2013; 19:2216–2223.

41. Miyamoto DT, Sequist LV, Lee RJ. Circulating tumour cells-monitoring treatment response in prostate cancer. Nat Rev Clin Oncol. 2014; 11:401–412.

42. Ghossein RA, Scher HI, Gerald WL, Kelly WK, Curley T, Amsterdam A, Zhang ZF, Rosai J. Detection of circulating tumor cells in patients with localized and metastatic prostatic carcinoma: clinical implications. J Clin Oncol. 1995; 13:1195–1200.

43. Roupret M, Hupertan V, Catto JW, Yates DR, Rehman I, Proctor LM, Phillips J, Meuth M, Cussenot O, Hamdy FC. Promoter hypermethylation in circulating blood cells identifies prostate cancer progression. Int J Cancer. 2008; 122:952–956.

44. Friedlander TW, Ngo VT, Dong H, Premasekharan G, Weinberg V, Doty S, Zhao Q, Gilbert EG, Ryan CJ, Chen WT, Paris PL. Detection and characterization of invasive circulating tumor cells derived from men with metastatic castration-resistant prostate cancer. Int J Cancer. 2014; 134:2284–2293.

45. Klotz L, Emberton M. Management of low risk prostate cancer-active surveillance and focal therapy. Nat Rev Clin Oncol. 2014; 11:324–334.

46. Pound CR, Partin AW, Eisenberger MA, Chan DW, Pearson JD, Walsh PC. Natural history of progression after PSA elevation following radical prostatectomy. JAMA. 1999; 281:1591–1597.

47. Brockman JA, Alanee S, Vickers AJ, Scardino PT, Wood DP, Kibel AS, Lin DW, Bianco FJ Jr, Rabah DM, Klein EA, Ciezki JP, Gao T, et al. Nomogram Predicting Prostate Cancer-specific Mortality for Men with Biochemical Recurrence After Radical Prostatectomy. Eur Urol. 2015; 67:1160–1167.

48. Simmons MN, Stephenson AJ, Klein EA. Natural history of biochemical recurrence after radical prostatectomy: risk assessment for secondary therapy. Eur Urol. 2007; 51:1175–1184.

49. Haffner MC, Mosbruger T, Esopi DM, Fedor H, Heaphy CM, Walker DA, Adejola N, Gurel M, Hicks J, Meeker AK, Halushka MK, Simons JW, Isaacs WB, et al. Tracking the clonal origin of lethal prostate cancer. J Clin Invest. 2013; 123:4918–4922.

50. Carter HB, Partin AW, Walsh PC, Trock BJ, Veltri RW, Nelson WG, Coffey DS, Singer EA, Epstein JI. Gleason score 6 adenocarcinoma: should it be labeled as cancer? J Clin Oncol. 2012; 30:4294–4296.

51. Pierorazio PM, Walsh PC, Partin AW, Epstein JI. Prognostic Gleason grade grouping: data based on the modified Gleason scoring system. BJU Int. 2013; 111:753–760.

52. Ploussard G, Isbarn H, Briganti A, Sooriakumaran P, Surcel CI, Salomon L, Freschi M, Mirvald C, van der Poel HG, Jenkins A, Ost P, van Oort IM, Yossepowitch O, et al. Members of the Prostate Cancer Working Group of the Young Academic Urologists Working Party of the European Association of U. Can we expand active surveillance criteria to include biopsy Gleason 3 + 4 prostate cancer? A multi-institutional study of 2,323 patients. Urol Oncol. 2015; 33:71 e71–79.

53. Knezevic D, Goddard AD, Natraj N, Cherbavaz DB, Clark-Langone KM, Snable J, Watson D, Falzarano SM, Magi-Galluzzi C, Klein EA, Quale C. Analytical validation of the Oncotype DX prostate cancer assay - a clinical RT-PCR assay optimized for prostate needle biopsies. BMC Genomics. 2013; 14:690.

54. Cuzick J, Berney DM, Fisher G, Mesher D, Moller H, Reid JE, Perry M, Park J, Younus A, Gutin A, Foster CS, Scardino P, Lanchbury JS, et al. Prognostic value of a cell cycle progression signature for prostate cancer death in a conservatively managed needle biopsy cohort. Br J Cancer. 2012; 106:1095–1099.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 10639