- Open Access
High levels of 5-hydroxymethylcytosine (5hmC) is an adverse predictor of biochemical recurrence after prostatectomy in ERG-negative prostate cancer
Clinical Epigeneticsvolume 7, Article number: 111 (2015)
Prostate cancer (PC) can be stratified into distinct molecular subtypes based on TMPRSS2-ERG gene fusion status, but its potential prognostic value remains controversial. Likewise, routine clinicopathological features cannot clearly distinguish aggressive from indolent tumors at the time of diagnosis; thus, new prognostic biomarkers are urgently needed. The DNA methylation variant 5-hydroxymethylcytosine (5hmC, an oxidized derivative of 5-methylcytosine) has recently emerged as a new diagnostic and/or prognostic biomarker candidate for several human malignancies. However, this remains to be systematically investigated for PC. In this study, we determined 5hmC levels in 311 PC (stratified by ERG status) and 228 adjacent non-malignant (NM) prostate tissue specimens by immunohistochemical analysis of a tissue microarray, representing a large radical prostatectomy (RP) cohort with long clinical follow-up. We investigated possible correlations between 5hmC and routine clinicopathological variables and assessed the prognostic potential of 5hmC by Kaplan-Meier and uni- and multivariate Cox regression analyses in ERG+ (n = 178) vs. ERG− (n = 133) PCs using biochemical recurrence (BCR) as endpoint.
We observed a borderline significant (p = 0.06) reduction in 5hmC levels in PC compared to NM tissue samples, which was explained by a highly significant (p < 0.001) loss of 5hmC in ERG− PCs. ERG status was not predictive of BCR in this cohort (p = 0.73), and no significant association was found between BCR and 5hmC levels in ERG+ PCs (p = 0.98). In contrast, high 5hmC immunoreactivity was a significant adverse predictor of BCR after RP in ERG− PCs, independent of Gleason score, pathological tumor stage, surgical margin status, and pre-operative prostate-specific antigen (PSA) level (hazard ratio (HR) (95 % confidence interval (CI)): 1.62 (1.15–2.28), p = 0.006).
This is the first study to demonstrate a prognostic potential for 5hmC in PC. Our findings highlight the importance of ERG stratification in PC biomarker studies and suggest that epigenetic mechanisms involving 5hmC are important for the development and/or progression of ERG− PC.
In 2012, more than 1 million men worldwide were diagnosed with prostate cancer (PC), accounting for 15 % of all male cancer diagnoses that year . While most PCs are indolent and rarely progress into clinically significant disease, a subset of PC patients develop highly aggressive metastatic disease with lethal outcome . At the time of diagnosis, currently available routine prognostic tools (mainly Gleason score (GS), serum prostate-specific antigen (PSA), and clinical tumor stage) are unable to accurately predict the outcome of PC. Accordingly, although widespread use of PSA testing has facilitated early detection of PC, it has also led to over diagnosis and overtreatment of many indolent tumors . Thus, novel and improved prognostic biomarkers for PC are urgently needed.
Approximately 50 % of all PCs carry the TMPRSS2-ERG gene fusion  that places the proto-oncogene ERG under androgen regulation, resulting in overexpression of this ETS family transcription factor. However, there is conflicting evidence as to whether the TMPRSS2-ERG fusion and/or the level of ERG expression have prognostic implications . Nevertheless, results from several studies indicate that distinct molecular mechanisms are at play in ERG-positive (ERG+) vs. ERG-negative (ERG−) PCs, including significant differences in the epigenome [6–14].
Methylation of the cytosine 5′ carbon (5-methylcytosine, 5mC) in CpG dinucleotides is the most extensively investigated epigenetic modification of the human genome. Whereas PC is characterized by a low frequency of somatic mutations, DNA methylation changes are considered one of its hallmarks , and several studies have shown that specific DNA methylation changes have promising potential as diagnostic and/or prognostic markers for PC [16–19]. However, the mechanisms of 5mC removal remained obscure until the recent discovery of the ten-eleven translocated (TET) proteins that convert 5mC to 5-hydroxymethylcytosine (5hmC) in an active oxidative demethylation process [20, 21]. 5hmC can be further oxidized by TETs to 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC), which are removed by base excision repair and replaced by an unmodified cytosine to complete the demethylation process [22, 23]. Active DNA demethylation is recognized as an important mechanism of plasticity in epigenetic regulation. However, accumulating evidence indicates that 5hmC is not only simply an intermediate in the process of active demethylation but also may function as a distinct epigenetic mark . Indeed, it has been demonstrated that 5mC and 5hmC interact with unique sets of chromatin binding proteins , suggesting that these epigenetic marks hold distinct roles in chromatin regulation and/or organization.
Because standard methods for DNA methylation analysis, such as genomic bisulfite sequencing, cannot distinguish between 5mC and 5hmC, the contribution of 5hmC to epigenetic regulation has been overlooked . Recent reports indicate that 5hmC levels are relatively high in terminally differentiated cells and lower in stem/progenitor cells [27–29]. Consistent with this, several studies have reported significantly reduced levels of 5hmC in cancer compared to non-malignant (NM) tissue samples [27, 29–39], and reduced levels of 5hmC have been suggested as a potential diagnostic marker for malignant transformation in several cancer types, including PC [27, 31, 32, 34, 37, 39]. However, for the latter, this is currently based only on two small-scale studies, reporting reduced 5hmC immunoreactivity in 30 PC vs. 10 NM  and 11 PC vs. 11 NM prostate tissue samples , respectively. To the best of our knowledge, there are no previous reports of a prognostic potential of 5hmC in PC, but low 5hmC levels have been associated with poor outcome in melanoma , myelodysplastic syndromes , gastric cancer , and glioblastoma , while high 5hmC levels have been associated with poor outcome in acute myeloid leukemia (AML) .
In the present study, we investigated the level of 5hmC in 311 malignant and 228 NM prostate tissue samples from a large radical prostatectomy (RP) patient cohort with 80 months median follow-up. We used a commercially available polyclonal anti-5hmC antibody with validated performance and specificity for 5hmC, as demonstrated in several published studies [21, 27–29, 37, 38, 41–44]. Our results indicate that 5hmC levels are significantly reduced in ERG− but not in ERG+ PCs, as compared to NM prostate tissue samples. Furthermore, we found that 5hmC had significant prognostic potential in ERG− but not in ERG+ PCs.
5hmC levels are significantly reduced in ERG− PC tissue samples
To systematically investigate 5hmC levels in NM and PC tissue samples, we analyzed a large RP tissue microarray by immunohistochemistry (clinical data in Table 1). Nuclear 5hmC staining intensity in prostate epithelial cells was evaluated and given a numerical grade (0, no staining; 1, moderate staining; 2, strong staining; Fig. 1a–d). A 5hmC score was then calculated as the mean grade of at least two malignant or 2 NM cores, respectively (5hmC score <1, weak; =1, moderate; >1, strong). Initially, we compared 5hmC scores for 311 patients for whom we could evaluate 5hmC staining in at least 2 malignant cores and 228 patients for whom 5hmC staining was assessable in 2 NM cores (Table 1). The vast majority (96 %) of NM samples displayed strong or intermediate 5hmC staining (Fig. 2). In the full patient set (n = 311), 5hmC levels were slightly reduced in PC compared to NM tissue samples (Fig. 1a), but the difference was only borderline significant (Fig. 2; p = 0.06, chi2 test). This finding is in agreement with results from two previous immunohistochemistry (IHC) studies that used relatively small PC patient sample sets [27, 37].
Next, we investigated 5hmC levels in PCs stratified by ERG status, as determined by an IHC-based method that has >95 % sensitivity and specificity for detecting ERG rearrangements in PC [45, 46]. For ERG− PCs (n = 133, Table 1), malignant cores showed significantly less strong (42 vs. 63 %), more moderate (46 vs. 33 %), and more weak (12 vs. 4 %) 5hmC staining compared to NM cores (Fig. 2), a difference that was highly statistically significant (p < 0.001, chi2 test; Fig. 2). In contrast, ERG+ PCs (n = 178, Table 1) displayed an almost identical score distribution as seen in NM cores (strong, 62 vs. 63 %; moderate, 34 vs. 33 %; weak, both 4 %; p = 0.98, chi2 test; Fig. 2). No difference was observed in 5hmC scores in NM tissue samples from patients with ERG+ (n = 106) vs. ERG− (n = 122) PC (p = 0.18, chi2 test; data not shown). Together, these results indicate that loss of 5hmC in PC occurs preferentially in ERG− cases.
Correlation between 5hmC levels and clinicopathological parameters in PC
Next, we investigated possible correlations between 5hmC score in malignant cores and routine clinicopathological parameters pre-operative PSA, GS, pathological tumor stage (pT), surgical margin (SM), and biochemical recurrence (BCR) status. In the full patient set, strong 5hmC staining (score >1) was significantly associated with BCR (p = 0.018, chi2 test), whereas no significant correlations were seen with any of the other parameters (p > 0.39, chi2 test; Additional file 1: Figure S1). Notably, by subgroup analysis, we found that strong 5hmC staining (score >1) was significantly associated with BCR in the ERG− (p = 0.006, chi2 test; Additional file 2: Figure S2) but not in the ERG+ subset (p = 0.45; Additional file 3: Figure S3). Strong 5hmC staining was also significantly associated with advanced pT stage (pT3-4) in ERG− (p = 0.001, chi2 test; Additional file 2: Figure S2), but not in ERG+ PCs, where instead a borderline significant trend towards reduced 5hmC staining in pT3-4 stage tumors was seen (p = 0.06, chi2 test; Additional file 3: Figure S3). There was no significant association between 5hmC score and pre-operative PSA, GS, or SM status in either the ERG+ or ERG− subset (p > 0.12; Additional file 2: Figure S2, Additional file 3: Figure S3). In summary, these findings suggest that high 5hmC levels may be associated with tumor progression in ERG− PC.
Prognostic value of 5hmC levels in PC
To assess the potential prognostic value of 5hmC levels in PC, we investigated whether 5hmC score in malignant cores was associated with time to BCR after RP. In the full patient set (n = 311), a high 5hmC score (analyzed as a continuous variable) was significantly associated with shorter time to BCR in univariate Cox regression analysis (hazard ratio (HR) (95 % confidence interval (CI)): 1.53 (1.09–2.15), p = 0.013, Table 2). Likewise, high GS, advanced pT stage, high pre-operative PSA, and positive SM status predicted early BCR in univariate analysis (p < 0.001, Table 2), whereas ERG status had no prognostic value in this patient set (p = 0.73, Table 2). Notably, 5hmC score remained significant also in multivariate Cox regression analysis (HR (95 % CI): 1.62 (1.15–2.28), p = 0.006, Table 2) together with all routine clinicopathological parameters (p ≤ 0.003, Table 2). We used Harrell’s C-index to estimate predictive accuracies of the multivariate models, which remained at 0.75 whether 5hmC score was included or not (Table 2). When 5hmC score was analyzed as a dichotomized variable (≤1 vs. >1) in the full patient set, it was only borderline significant in uni- and multivariate cox regression analysis (p = 0.059 and p = 0.066, respectively; Additional file 4: Table S1) and in Kaplan-Meier analysis (p = 0.058; Fig. 3a)
Recurrence-free survival analyses were also performed after patient stratification by ERG status. Pre-operative PSA, GS, pT stage, and SM status were significant predictors of time to BCR in univariate analysis in both patient subgroups (p < 0.05, Table 3), indicating that this is a representative RP cohort. Furthermore, a high 5hmC score (analyzed as a continuous variable) was a significant adverse predictor of BCR in the ERG− subgroup in univariate (HR (95 % CI), 1.97 (1.19–3.26), p = 0.008) as well as multivariate Cox regression analysis (HR (95 % CI), 2.02 (1.16–3.51), p = 0.013), whereas no significant association was found for ERG+ PCs (Table 3). For the ERG− PC subgroup, addition of 5hmC score to a multivariate model containing clinicopathological factors only (pre-operative PSA and SM status) improved Harrell’s C-index from 0.68 to 0.73, suggesting moderately improved predictive accuracy. Pathological T stage and GS were excluded from the final model because they failed in multivariate analysis (Table 3). Similar results were obtained when 5hmC score was analyzed as a dichotomized variable (≤1 vs. >1) in uni- and multivariate Cox regression analysis (Additional file 5: Table S2). A high 5hmC score (>1) was also significantly associated with shorter time to BCR in Kaplan-Meier analysis in ERG− PCs (p = 0.003, log-rank test; Fig. 3b), whereas no difference was observed between the high/low 5hmC score subgroups for ERG+ PCs (p = 0.95, log-rank test; Fig. 3c).
Finally, we note that despite the significant association between high 5hmC score and advanced pT stage in ERG− PCs (Additional file 2: Figure S2), both 5hmC score and pT stage remained significant predictors of time to BCR in a bivariate Cox regression model (5hmC score; HR (95 % CI), 1.67 (1.01–2.81), p = 0.044; pT stage, HR (95 % CI), 1.96 (1.14–3.35), p = 0.014). The C-index for this model was 0.66, i.e., higher than for either parameter alone (Table 3). We also performed Kaplan-Meier analysis separately for pT2 and pT3-4 stage ERG− PCs, respectively. A high 5hmC score (≥1) was significantly associated with BCR in the pT2 stage subgroup (p = 0.029, log-rank test, Additional file 6: Figure S4a) and a similar trend was observed for the pT3-4 subgroup, although statistical significance was not reached in this smaller patient subset (p = 0.35, log-rank test, Additional file 6: Figure S4b). In conclusion, our results indicate that a high 5hmC level is a significant adverse predictor of time to BCR after RP in patients with ERG− PC, independent of routine clinicopathological variables.
The present study is the first large-scale investigation of 5hmC levels in PC, as well as the first study to explore and demonstrate a prognostic potential for 5hmC in this common malignancy. By immunohistochemical analysis of a large RP cohort, we found that 5hmC levels were significantly reduced in ERG− but not in ERG+ PC as compared to NM prostate tissue samples. Furthermore, we identified 5hmC as a significant predictor of biochemical recurrence in ERG− PC, independent of routine clinicopathological factors. In contrast, 5hmC had no prognostic value in ERG+ PCs. ERG status was also not associated with BCR in this RP patient cohort, consistent with previous reports .
Our findings confirm and expand on results from two previous studies [27, 37] that reported loss of 5hmC immunoreactivity in PC, but based on analysis of only 30 PC vs. 10 NM and 11 PC vs. 11 NM tissue samples, respectively. The prognostic potential of 5hmC was not explored in these studies, nor was the association between 5hmC and ERG status. The small-scale study design also prohibited an evaluation of 5hmC levels in relation to routine clinicopathological parameters in these earlier reports [27, 37]. Correlations between reduced 5hmC levels and mutations in 5hmC regulating enzymes, such as TETs, have been identified in various other human malignancies [29, 30, 34, 35], but such mutations are rare in PC [47, 48]. Thus, further studies are needed to investigate how 5hmC levels are regulated in normal as well as in ERG− and ERG+ PC cells.
Our results highlight the importance of ERG stratification in PC biomarker studies. Although not understood in detail, the molecular pathways that operate in ERG+ tumors are relatively better described than the mechanisms involved in development/progression of ERG− PCs . It is becoming increasingly clear that ERG+ and ERG− PCs represent distinct molecular subtypes and that subtype-specific markers may be required for prediction of tumor aggressiveness . Thus, similar to our findings for 5hmC, previous studies have shown that the prognostic potential of certain PC biomarker candidates is dependent on ERG status, e.g., SPINK1 and SPOP as well as some DNA methylation marker candidates [6–12].
It is intriguing that while we observed a general loss of 5hmC in ERG− PCs, our results indicated that retaining a high global level of 5hmC is associated with ERG− PC progression. Interestingly, hypoxia has previously been associated with poor prognosis after RP  and, together with reactive oxygen species (ROS), has been shown to induce TET1 expression, leading to activation of hypoxia-inducible genes through conversion of 5mC to 5hmC in both a global and site-specific manner [50–52]. While we did not find significant differential expression of TET1 or common hypoxia response genes (HIF1A, HIF2A, CA9, PGK1, GPI, VEGFA, BNIP, ENO1) [50–52] in ERG+ vs. ERG− PCs in publicly available datasets [53–55] (data not shown), further studies are warranted to investigate whether high levels of hypoxia may exist in a subset of ERG− PC and thus potentially could be linked to elevated 5hmC levels and poor prognosis.
In line with our finding that high levels of 5hmC were associated with early BCR in ERG− PC, high levels of 5hmC have also been identified as an independent predictor of poor prognosis in AML . Furthermore, the same study reported highly variable 5hmC levels between different AML samples , consistent with our observations for both ERG− and ERG+ PCs. Likewise, 5hmC levels have been shown to vary between different types of brain cancer, with oligodendroglial tumors displaying overall high 5hmC levels, while reduced 5hmC levels were found in adult glioblastoma and anaplastic astrocytoma and associated with poor prognosis . Loss of 5hmC has also been reported to predict adverse outcome in melanoma , myelodysplastic syndromes , and gastric cancer . Together, results from these previous reports and our current results for 5hmC in PC indicate that potential prognostic implications of 5hmC are cancer (sub)type-specific, which is also in accordance with the highly tissue-specific distribution of 5hmC in normal cells [33, 43].
There are some limitations to the present study. First, the prognostic value of 5hmC in ERG− but not in ERG+ PC needs further validation in large independent patient cohorts. Moreover, we used BCR as endpoint, which is only a surrogate for tumor aggressiveness. Thus, the potential prognostic value of 5hmC should also be assessed in relation to more clinically relevant endpoints, such as metastatic progression or PC-specific and overall mortality. Due to the slowly progressing nature of PC, analysis of such endpoints would require at least 15 years of follow-up .
Also limiting our study is the lack of 5mC data for our sample set; thus, we have not investigated the correlation between 5mC and 5hmC levels. However, previous studies have reported no clear correlation between these two epigenetic marks [27, 43]. Furthermore, we have not investigated the expression levels of enzymes involved in 5hmC regulation, such as TETs or isocitrate dehydrogenases (IDHs) , which should be investigated in future studies.
Another potential limitation of our study is the use of ERG immunoreactivity as a proxy for ERG gene fusion status. However, previous studies have demonstrated that ERG immunoreactivity is highly concordant (>95 %) with ERG rearrangement status in PC tissue samples [45, 46]. Moreover, ERG overexpression is considered a PC “driver” event whether caused by translocation, copy number alteration, or other mechanisms . Furthermore, while ERG status may differ between distinct cancer foci in multifocal PC , we are confident of the match between 5hmC score and ERG status in the present study, as these were assessed on consecutive sections of the exact same cores. Finally, although this is the first large-scale evaluation of 5hmC in PC, immunohistochemical analysis is semi-quantitative and only allows assessment of global 5hmC levels. Our results warrant further studies of the genomic distribution of 5hmC in ERG− vs. ERG+ PC as well as in NM prostate tissue samples. A number of recently developed NGS-based protocols allow quantitative genome wide or whole-genome profiling of the 5hmC methylome , which could be used to increase our understanding of the epigenetic mechanisms involved in PC development and progression.
This is the first large-scale study of 5hmC in PC as well as the first study to demonstrate a prognostic potential for 5hmC in prostate adenocarcinoma. In conclusion, we found that the global level of 5hmC was significantly reduced in ERG− but not in ERG+ PC, as compared to NM prostate tissue samples. Furthermore, we found that a subgroup of ERG− tumors retained a high 5hmC level that was associated with early BCR after RP, whereas 5hmC had no significant prognostic value in ERG+ PC in our patient set. Our results highlight the importance of ERG stratification in PC biomarker studies and suggest that epigenetic mechanisms involving 5hmC are important in ERG− PC tumorigenesis and development. Future studies, using large independent PC patient cohorts, are needed to confirm our findings. Finally, in order to better understand the role of epigenetic changes in PC development and progression, it will be important to generate genome wide maps of the 5-hydroxymethylome in malignant and NM tissue samples, as well as in specimens from aggressive and non-aggressive PC, while taking ERG status into account. Such future studies could aid in the identification of important genes and pathways involved in PC development and progression, providing a deeper understanding of epigenomic reprogramming in PC.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files.
Radical prostatectomy TMA
A tissue microarray (TMA) was generated using formalin-fixed paraffin-embedded tissue samples from 552 radical prostatectomies of histological verified clinically localized PC, surgically removed with curative intent between 1998 and 2009 at Department of Urology, Aarhus University Hospital, Denmark. In all cases, Gleason scores were reassigned according to International Union Against Cancer and WHO/International Society of Urological Pathology criteria . For each patient, a trained pathologist with extensive experience in prostate histopathology (SH) identified a representative malignant tissue area in a hematoxylin and eosin-stained section from the original prostatectomy specimens. Likewise, a representative area of adjacent NM prostate tissue was selected in parallel for a subset of the patients (n = 301, chosen at random). For each patient, three cores (1 mm diameter) were punched from the selected malignant tissue area and two cores were punched from the NM tissue area. The cores were mounted in a total of 16 individual TMA blocks using the TMA master (3DHISTECH, Hungary).
The most recent clinical follow-up of time to PSA recurrence after RP was conducted in May 2015 for all 552 patients on the TMA. At this time, 6 patients had withdrawn consent, while another 101 patients were excluded due to either pre- or post-operative endocrine or radiation treatment, short (<3 months) follow-up, or BCR within 3 months from RP (Additional file 7: Figure S5). Clinical data for all patients included on the TMA, except for the six who withdrew consent, is listed in Table 1.
This study was approved by the regional ethical committee and the Danish Data Protection Agency. All patients provided written informed consent.
5hmC staining and evaluation
TMA tissue sections (4 μm) were deparaffinized (Tissue-Tek Tissue-Clear Xylene Substitute, Sakura) and re-hydrated according to standard protocols and antigen retrieval was carried out in two subsequent steps, as described previously . Briefly, heat induced epitope retrieval with citrate buffer (pH 6.0) was performed, followed by incubation in 3.5 N HCl for 15 min at room temperature. The subsequent procedure was performed on the Autostainer Link48 (DAKO, Denmark) at room temperature: rabbit polyclonal 5-hmC specific antibody (Active Motif, Carlsbad, CA, cat. no 39769) was applied at a 1:1000 dilution and incubated for 60 min. Horse radish peroxidase (HRP) conjugated rabbit secondary antibody (Envision) was applied for 30 min. TMA sections were subsequently counterstained with hematoxylin. As a negative control, the primary antibody was omitted. This primary antibody and IHC staining protocol has been shown to be highly 5hmC specific in multiple previously published reports [21, 27–29, 37, 38, 41–44], including several antigen competition experiments, thus proving the sensitivity and specificity of this antibody for 5hmC.
For each core on the TMA, 5hmC immunoreactivity was scored individually by two independent observers (SH and SHS) using the Pannoramic Viewer software (3DHISTECH, Hungary). In cases of disagreement, cores were reassessed to reach a consensus score. Kappa statistics showed very high inter-observer agreement (Kappa index = 0.95; p < 0.0001). Nuclear 5hmC staining intensity in malignant or NM prostate epithelial cells was evaluated and given a numerical grade (0, no staining; 1, moderate staining; 2, strong staining). Due to progressive TMA slicing (sections used for other projects), some cores had changed status from malignant to NM from the time of TMA construction. Thus, for 70 patients, all malignant cores were lost during TMA processing and 64 patients had less than two malignant cores that could be evaluated for 5hmC staining. Patients for whom we could not evaluate 5hmC staining intensity in at least two NM or at least two PC cores were excluded from further analysis. Then, for each patient, a 5hmC score was calculated as the mean grade of at least two malignant cores or two NM cores, respectively (score <1, weak; 1, moderate; >1, strong). In total, a 5hmC score was calculated in malignant tissue from 311 patients (Additional file 7: Figure S5) and a 5hmC score was calculated in NM tissue from 228 patients (Additional file 8: Figure S6).
TMA tissue sections (2.5 μm) were deparaffinized and epitopes were retrieved using TEG buffer, as previously described . TMA sections were stained with rabbit monoclonal ERG antibody (2805–1, Epitomics) in a 1:150 dilution in TEG. Secondary staining was performed using the EnVision + System (HRP Labelled Polymer Anti-Rabbit K4003, DakoCytomation). For all cores, ERG immunoreactivity (0,no staining; 1, moderate; 2, strong) was evaluated and scored by two independent observers (SH and ASL) as described above. Positive staining was used as a proxy for ERG fusion status as previously described [45, 46].
All statistical analyses were performed using STATA v. 11.2 (StataCorp, College Station TX, USA). In all cases, p < 0.05 was considered significant. Associations between 5hmC score and clinicopathological variables were assessed by the chi2 test. Biochemical recurrence (BCR), defined as PSA >0.2 ng/ml on two consecutive measurements, was used as the clinical endpoint in univariate and multivariate Cox regression analyses as well as in Kaplan-Meier analyses. Patients without BCR were censored at their last normal PSA measurement. Statistical significance in Kaplan-Meier analysis was evaluated using two-sided log-rank tests. Predictive accuracy was estimated using Harrell’s C-index.
acute myeloid leukemia
pathological tumor stage
International Agency for Research on Cancer, W.H.O. Globocan 2012. http://globocan.iarc.fr/Pages/fact_sheets_cancer.aspx. Accessed on May 14, 2015.
Albertsen PC, Hanley JA, Fine J. 20-year outcomes following conservative management of clinically localized prostate cancer. JAMA. 2005;293:2095–101.
Telesca D, Etzioni R, Gulati R. Estimating lead time and overdiagnosis associated with PSA screening from prostate cancer incidence trends. Biometrics. 2008;64:10–9.
Tomlins SA, Rhodes DR, Perner S, Dhanasekaran SM, Mehra R, Sun XW, et al. Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer. Science. 2005;310:644–8.
Bostrom, P.J.; Bjartell, A.S.; Catto, J.W.; Eggener, S.E.; Lilja, H.; Loeb, S.; Schalken, J.; Schlomm, T.; Cooperberg, M.R. Genomic predictors of outcome in prostate cancer. European urology 2015: doi: 10.1016/j.eururo.2015.04.008
Barbieri CE, Baca SC, Lawrence MS, Demichelis F, Blattner M, Theurillat JP, et al. Exome sequencing identifies recurrent SPOP, FOXA1 and MED12 mutations in prostate cancer. Nat Genet. 2012;44:685–9.
Brase JC, Johannes M, Mannsperger H, Falth M, Metzger J, Kacprzyk LA, et al. TMPRSS2-ERG-specific transcriptional modulation is associated with prostate cancer biomarkers and TGF-beta signaling. BMC Cancer. 2011;11:507.
Gasi Tandefelt D, Boormans JL, van der Korput HA, Jenster GW, Trapman J. A 36-gene signature predicts clinical progression in a subgroup of ERG-positive prostate cancers. Eur Urol. 2013;64:941–50.
Karnes RJ, Cheville JC, Ida CM, Sebo TJ, Nair AA, Tang H, et al. The ability of biomarkers to predict systemic progression in men with high-risk prostate cancer treated surgically is dependent on ERG status. Cancer Res. 2010;70:8994–9002.
Kron K, Liu L, Trudel D, Pethe V, Trachtenberg J, Fleshner N, et al. Correlation of ERG expression and DNA methylation biomarkers with adverse clinicopathologic features of prostate cancer. Clin Cancer Res. 2012;18:2896–904.
Tomlins SA, Rhodes DR, Yu J, Varambally S, Mehra R, Perner S, et al. The role of SPINK1 in ETS rearrangement-negative prostate cancers. Cancer Cell. 2008;13:519–28.
Vinarskaja A, Schulz WA, Ingenwerth M, Hader C, Arsov C. Association of PITX2 mRNA down-regulation in prostate cancer with promoter hypermethylation and poor prognosis. Urol Oncol. 2013;31:622–7.
Kron K, Trudel D, Pethe V, Briollais L, Fleshner N, van der Kwast T, et al. Altered DNA methylation landscapes of polycomb-repressed loci are associated with prostate cancer progression and ERG oncogene expression in prostate cancer. Clin Cancer Res. 2013;19:3450–61.
Tomlins SA, Alshalalfa M, Davicioni E, Erho N, Yousefi K, Zhao S, et al. Characterization of 1577 primary prostate cancers reveals novel biological and clinicopathologic insights into molecular subtypes. Eur Urol. 2015;68(4):555–67.
Dobosy JR, Roberts JL, Fu VX, Jarrard DF. The expanding role of epigenetics in the development, diagnosis and treatment of prostate cancer and benign prostatic hyperplasia. J Urol. 2007;177:822–31.
Chiam K, Ricciardelli C, Bianco-Miotto T. Epigenetic biomarkers in prostate cancer: current and future uses. Cancer Lett. 2014;342:248–56.
Haldrup C, Mundbjerg K, Vestergaard EM, Lamy P, Wild P, Schulz WA, et al. DNA methylation signatures for prediction of biochemical recurrence after radical prostatectomy of clinically localized prostate cancer. J Clin Oncol. 2013;31:3250–8.
Kristensen H, Haldrup C, Strand S, Mundbjerg K, Mortensen MM, Thorsen K, 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–81.
Strand SH, Orntoft TF, Sorensen KD. Prognostic DNA methylation markers for prostate cancer. Int J Mol Sci. 2014;15:16544–76.
Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, et al. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009;324:930–5.
Ito S, D'Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y. Role of Tet proteins in 5mc to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature. 2010;466:1129–33.
Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, et al. Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science. 2011;333:1300–3.
He YF, Li BZ, Li Z, Liu P, Wang Y, Tang Q, et al. Tet-mediated formation of 5-carboxylcytosine and its excision by TDG in mammalian DNA. Science. 2011;333:1303–7.
Ye C, Li L. 5-hydroxymethylcytosine: a new insight into epigenetics in cancer. Cancer Biol Ther. 2014;15:10–5.
Iurlaro M, Ficz G, Oxley D, Raiber EA, Bachman M, Booth MJ, et al. A screen for hydroxymethylcytosine and formylcytosine binding proteins suggests functions in transcription and chromatin regulation. Genome Biol. 2013;14:R119.
Huang Y, Pastor WA, Shen Y, Tahiliani M, Liu DR, Rao A. The behaviour of 5-hydroxymethylcytosine in bisulfite sequencing. PLoS One. 2010;5:e8888.
Haffner MC, Chaux A, Meeker AK, Esopi DM, Gerber J, Pellakuru LG, et al. Global 5-hydroxymethylcytosine content is significantly reduced in tissue stem/progenitor cell compartments and in human cancers. Oncotarget. 2011;2:627–37.
Haffner MC, Pellakuru LG, Ghosh S, Lotan TL, Nelson WG, De Marzo AM, et al. Tight correlation of 5-hydroxymethylcytosine and polycomb marks in health and disease. Cell Cycle. 2013;12:1835–41.
Orr BA, Haffner MC, Nelson WG, Yegnasubramanian S, Eberhart CG. Decreased 5-hydroxymethylcytosine is associated with neural progenitor phenotype in normal brain and shorter survival in malignant glioma. PLoS One. 2012;7:e41036.
Ko M, Huang Y, Jankowska AM, Pape UJ, Tahiliani M, Bandukwala HS, et al. Impaired hydroxylation of 5-methylcytosine in myeloid cancers with mutant TET2. Nature. 2010;468:839–43.
Kraus TF, Globisch D, Wagner M, Eigenbrod S, Widmann D, Munzel M, et al. Low values of 5-hydroxymethylcytosine (5hmc), the "sixth base," are associated with anaplasia in human brain tumors. Int J Cancer. 2012;131:1577–90.
Kudo Y, Tateishi K, Yamamoto K, Yamamoto S, Asaoka Y, Ijichi H, et al. Loss of 5-hydroxymethylcytosine is accompanied with malignant cellular transformation. Cancer Sci. 2012;103:670–6.
Li W, Liu M. Distribution of 5-hydroxymethylcytosine in different human tissues. J Nucleic Acids. 2011;2011:870726.
Lian CG, Xu Y, Ceol C, Wu F, Larson A, Dresser K, et al. Loss of 5-hydroxymethylcytosine is an epigenetic hallmark of melanoma. Cell. 2012;150:1135–46.
Liu X, Zhang G, Yi Y, Xiao L, Pei M, Liu S, et al. Decreased 5-hydroxymethylcytosine levels are associated with TET2 mutation and unfavorable overall survival in myelodysplastic syndromes. Leuk Lymphoma. 2013;54:2466–73.
Uribe-Lewis S, Stark R, Carroll T, Dunning MJ, Bachman M, Ito Y, et al. 5-hydroxymethylcytosine marks promoters in colon that resist DNA hypermethylation in cancer. Genome Biol. 2015;16:69.
Yang H, Liu Y, Bai F, Zhang JY, Ma SH, Liu J, et al. Tumor development is associated with decrease of TET gene expression and 5-methylcytosine hydroxylation. Oncogene. 2013;32:663–9.
Yang Q, Wu K, Ji M, Jin W, He N, Shi B, et al. Decreased 5-hydroxymethylcytosine (5-hmC) is an independent poor prognostic factor in gastric cancer patients. J Biomed Nanotechnol. 2013;9:1607–16.
Jin SG, Jiang Y, Qiu R, Rauch TA, Wang Y, Schackert G, et al. 5-Hydroxymethylcytosine is strongly depleted in human cancers but its levels do not correlate with IDH1 mutations. Cancer Res. 2011;71:7360–5.
Kroeze LI, Aslanyan MG, van Rooij A, Koorenhof-Scheele TN, Massop M, Carell T, et al. Characterization of acute myeloid leukemia based on levels of global hydroxymethylation. Blood. 2014;124:1110–8.
Globisch D, Munzel M, Muller M, Michalakis S, Wagner M, Koch S, et al. Tissue distribution of 5-hydroxymethylcytosine and search for active demethylation intermediates. PLoS One. 2010;5:e15367.
Matsuda I, Imai Y, Hirota S. Distinct global DNA methylation status in B-cell lymphomas: immunohistochemical study of 5-methylcytosine and 5-hydroxymethylcytosine. J Clin Exp Hematop. 2014;54:67–73.
Nestor CE, Ottaviano R, Reddington J, Sproul D, Reinhardt D, Dunican D, et al. Tissue type is a major modifier of the 5-hydroxymethylcytosine content of human genes. Genome Res. 2012;22:467–77.
Nettersheim D, Heukamp LC, Fronhoffs F, Grewe MJ, Haas N, Waha A, et al. Analysis of TET expression/activity and 5mc oxidation during normal and malignant germ cell development. PLoS One. 2013;8:e82881.
Park K, Tomlins SA, Mudaliar KM, Chiu YL, Esgueva R, Mehra R, et al. Antibody-based detection of erg rearrangement-positive prostate cancer. Neoplasia. 2010;12:590–8.
Braun M, Goltz D, Shaikhibrahim Z, Vogel W, Bohm D, Scheble V, et al. ERG protein expression and genomic rearrangement status in primary and metastatic prostate cancer—a comparative study of two monoclonal antibodies. Prostate Cancer Prostatic Dis. 2012;15:165–9.
Kang MR, Kim MS, Oh JE, Kim YR, Song SY, Seo SI, et al. Mutational analysis of IDH1 codon 132 in glioblastomas and other common cancers. Int J Cancer. 2009;125:353–5.
Kroeze LI, van der Reijden BA, Jansen JH. 5-hydroxymethylcytosine: an epigenetic mark frequently deregulated in cancer. Biochim Biophys Acta. 1855;2015:144–54.
Vergis R, Corbishley CM, Norman AR, Bartlett J, Jhavar S, Borre M, et al. Intrinsic markers of tumour hypoxia and angiogenesis in localised prostate cancer and outcome of radical treatment: a retrospective analysis of two randomised radiotherapy trials and one surgical cohort study. Lancet Oncol. 2008;9:342–51.
Kang KA, Piao MJ, Kim KC, Kang HK, Chang WY, Park IC, et al. Epigenetic modification of Nrf2 in 5-fluorouracil-resistant colon cancer cells: involvement of TET-dependent DNA demethylation. Cell Death Dis. 2014;5:e1183.
Mariani CJ, Vasanthakumar A, Madzo J, Yesilkanal A, Bhagat T, Yu Y, et al. TET1-mediated hydroxymethylation facilitates hypoxic gene induction in neuroblastoma. Cell Rep. 2014;7:1343–52.
Tsai YP, Chen HF, Chen SY, Cheng WC, Wang HW, Shen ZJ, et al. TET1 regulates hypoxia-induced epithelial-mesenchymal transition by acting as a co-activator. Genome Biol. 2014;15:513.
Zhu Y, Qiu P, Ji Y. TCGA-assembler: open-source software for retrieving and processing TCGA data. Nat Methods. 2014;11:599–600.
Zhu Y, Xu Y, Helseth Jr DL, Gulukota K, Yang S, Pesce LL, et al. Zodiac: a comprehensive depiction of genetic interactions in cancer by integrating TCGA data. J Natl Cancer Inst. 2015;107.
Long Q, Xu J, Osunkoya AO, Sannigrahi S, Johnson BA, Zhou W, et al. Global transcriptome analysis of formalin-fixed prostate cancer specimens identifies biomarkers of disease recurrence. Cancer Res. 2014;74:3228–37.
Adamo P, Ladomery MR. The oncogene ERG: a key factor in prostate cancer. Oncogene. 2015.
Burdova A, Bouchal J, Tavandzis S, Kolar Z. TMPRSS2-ERG gene fusion in prostate cancer. Biomed Pap Med Fac Univ Palacky Olomouc Czechoslovakia. 2014;158:502–10.
Yardimci H, Zhang Y. Charting oxidized methylcytosines at base resolution. Nat Struct Mol Biol. 2015;22:656–61.
Epstein JI, Allsbrook Jr WC, Amin MB, Egevad LL, Committee IG. The 2005 International Society of Urological Pathology (ISUP) Consensus Conference on Gleason Grading of Prostatic Carcinoma. Am J Surg Pathol. 2005;29:1228–42.
Sorensen KD, Abildgaard MO, Haldrup C, Ulhoi BP, Kristensen H, Strand S, et al. Prognostic significance of aberrantly silenced ANPEP expression in prostate cancer. Br J Cancer. 2013;108:420–8.
This work was supported by The Danish Strategic Research Council, Innovation Fund Denmark, and Aarhus University (fellowship for SHS). The Danish Cancer Biobank (DCB) is acknowledged for biological material. The authors wish to thank Pamela Celis, Kristina Mortensen, Kristina Lystlund Lauridsen, Nadia Knudsen, and Susanne Skou for excellent technical assistance.
The authors declare no competing interests.
KDS and SHS conceptualized and designed the study. ASL and SH constructed the tissue microarray. TMA scoring was done by SHS, ASL, and SH. SHS performed all statistical analyses under supervision from KDS. TMS, CH, MB, and TFØ contributed to the acquisition and interpretation of data. The manuscript was written by SHS and KDS. All authors read, commented, and approved the final version of the manuscript.
Distribution of 5hmC scores by clinicopathological parameters and BCR in the full PC patient set. Full PC patient set, n = 311. P values: chi2 test. (TIFF 145 kb)
Distribution of 5hmC scores by clinicopathological parameters and BCR status in ERG− PC. ERG− subset: n = 133. P values: chi2 test. (TIFF 146 kb)
Distribution of 5hmC scores by clinicopathological parameters and BCR status in ERG+ PC. ERG+ subset: n = 178. P values: chi2 test. (TIFF 143 kb)
Uni- and multivariate Cox regression analysis of BCR after RP in the full PC patient set. 5hmC analyzed as a dichotomized variable. (DOCX 33 kb)
Uni- and multivariate Cox regression analysis of BCR after RP in ERG−/ERG+ PC patients. 5hmC analyzed as a dichotomized variable. a Global multivariate model including all parameters. b Final multivariate model including significant variables only. c Harrell’s C-index for final model including 5hmC. d Harrell’s C-index for final model excluding 5hmC. NA: not applicable. Significant p values (p < 0.05) are highlighted in bold. (DOCX 36 kb)
Kaplan-Meier analysis: 5hmC score and time to BCR in ERG− PC stratified by pT stage. A: High 5hmC score (>1) was a significant adverse predictor of time to BCR in pT2 stage ERG− PCs (p = 0.029, log-rank test). B: The same trend was seen in pT3-4 stage ERG− PCs, but this was not statistically significant (p = 0.35, log-rank test). (TIFF 183 kb)
Flow chart illustrating the sample inclusion/exclusion process (malignant cores). A total of 552 PC patients were represented by malignant cores on the TMA. The final sample set used for all subsequent analyses consisted of 5hmC scores from 311 PC patients for whom at least 2 malignant cores could be evaluated for 5hmC staining intensity (grade). (TIFF 68 kb)
Flow chart illustrating the sample inclusion/exclusion process (NM cores). A total of 301 PC patients were represented by NM cores on the TMA. The final sample set used for all subsequent analyses consisted of 5hmC scores from 228 PC patients for whom at least 2 NM cores could be evaluated for 5hmC staining intensity (grade). (TIFF 58 kb)