- Open Access
A novel DNA methylation signature is associated with androgen receptor activity and patient prognosis in bone metastatic prostate cancer
Clinical Epigenetics volume 13, Article number: 133 (2021)
Patients with metastatic prostate cancer (PC) are treated with androgen deprivation therapy (ADT) that initially reduces metastasis growth, but after some time lethal castration-resistant PC (CRPC) develops. A better understanding of the tumor biology in bone metastases is needed to guide further treatment developments. Subgroups of PC bone metastases based on transcriptome profiling have been previously identified by our research team, and specifically, heterogeneities related to androgen receptor (AR) activity have been described. Epigenetic alterations during PC progression remain elusive and this study aims to explore promoter gene methylation signatures in relation to gene expression and tumor AR activity.
Materials and methods
Genome-wide promoter-associated CpG methylation signatures of a total of 94 tumor samples, including paired non-malignant and malignant primary tumor areas originating from radical prostatectomy samples (n = 12), and bone metastasis samples of separate patients with hormone-naive (n = 14), short-term castrated (n = 4) or CRPC (n = 52) disease were analyzed using the Infinium Methylation EPIC arrays, along with gene expression analysis by Illumina Bead Chip arrays (n = 90). AR activity was defined from expression levels of genes associated with canonical AR activity.
Integrated epigenome and transcriptome analysis identified pronounced hypermethylation in malignant compared to non-malignant areas of localized prostate tumors. Metastases showed an overall hypomethylation in relation to primary PC, including CpGs in the AR promoter accompanied with induction of AR mRNA levels. We identified a Methylation Classifier for Androgen receptor activity (MCA) signature, which separated metastases into two clusters (MCA positive/negative) related to tumor characteristics and patient prognosis. The MCA positive metastases showed low methylation levels of genes associated with canonical AR signaling and patients had a more favorable prognosis after ADT. In contrast, MCA negative patients had low AR activity associated with hypermethylation of AR-associated genes, and a worse prognosis after ADT.
A promoter methylation signature classifies PC bone metastases into two groups and predicts tumor AR activity and patient prognosis after ADT. The explanation for the methylation diversities observed during PC progression and their biological and clinical relevance need further exploration.
Prostate cancer (PC) is a very common malignancy and a major cause of cancer mortality in men worldwide (https://gco.iarc.fr). Most patients with lethal PC develop bone metastases, which are treated with androgen deprivation therapy (ADT). Initially, ADT reduces metastasis growth, but eventually castration-resistant PC (CRPC) develops. Several treatment strategies for metastatic CRPC exist, with the majority aiming at inhibiting androgen receptor (AR) signaling . Patients with CRPC show diverse responses to AR-inhibiting drugs, ranging from a strong response to complete resistance, underlining the need for complementary/alternative treatments. To guide future therapeutic developments, the tumor biology of metastatic PC must be understood in more detail.
The genetics behind CRPC have been comprehensively explored [2,3,4,5], indicating AR amplification as the most common event (seen in about 50% of the cases). Also, structural AR rearrangements and AR mutations are frequently observed together with mutations in co-factors regulating the transcriptional activity of the AR. The prognostic and therapy-predicting value of those genetic defects in relation to AR targeting therapies is still unclear [5, 6].
By exploring the gene expression in clinical samples of PC bone metastases, we have identified transcriptomic as well as proteomic profiles of suggested clinical relevance [7,8,9,10,11,12]. Specifically, we have identified three molecular subtypes of bone metastases, termed MetA–C, based on the diverged expression of genes related to AR activity, cell cycle activity, and stroma response . The MetA subtype shows high AR activity in comparison to MetB and C, while MetB shows higher cell cycle activity than the other subtypes. Accordingly, MetA patients have the best prognosis after ADT, while MetB patients have by far the worst prognosis. Interestingly, the MetA–C subtypes could be intrinsic, as selected features of MetB (low AR activity and high proliferation) were traceable back to the corresponding primary tumors . In line with this, PC subtypes with characteristics similar to the MetA–C subtypes have been identified by transcriptomic analysis of primary prostate tumors [13, 14].
The mechanistic explanation for the development of different metastasis subtypes is unknown. The subtypes could relate to diverse basal/luminal cellular origin, genetic defects as described above, but also to epigenetic defects affecting differentiation and clonal expansion, possibly influenced by the metastatic microenvironment and/or by therapy. Alterations in DNA methylation have been described to contribute to PC development and progression into metastatic disease [15,16,17]. The aim of the current study was to explore the promoter DNA methylation pattern in relation to gene expression during PC disease progression. Specific interest was put into analyzing promoter methylation levels in relation to canonical AR activity in metastasis samples and to patient prognosis.
Material and methods
Samples of bone metastases were obtained from a series of fresh-frozen biopsies collected from 70 patients with PC operated for metastatic spinal cord compression at Umeå University Hospital (2003–2013). The patient series have been previously described [7, 9, 12] and clinical characteristics of patients included in the current study are summarized in Table 1. The study also included 12 separate patients with localized PC who were treated with radical prostatectomy at Umeå University Hospital, between 2008 and 2009; mean age for these men was 60 years (range 48–68 years) and mean serum level of prostate specific antigen (PSA) was 11 ng/ml (range 3.5–26 ng/ml). Clinical local stage was T2 (n = 4) or T3 (n = 8) and Gleason score (GS) was 7 (n = 10) or 8 (n = 2).
Tissue sample preparation
Bone metastasis samples were instantly fresh-frozen in liquid nitrogen. Fresh radical prostatectomy specimens were received at the pathology department immediately after surgery and cut into 0.5 cm thick slices. From these slices, 20 samples were taken using a 0.5 cm skin punch and frozen in liquid nitrogen within 30 min after surgery. The prostate slices were formalin-fixed, embedded in paraffin, cut in 5 µm sections, whole-mounted, and stained with hematoxylin–eosin. Tissue sample composition (non-malignant or malignant) of the frozen pieces was determined from their location in the whole-mount sections.
Representative areas of fresh-frozen bone metastasis (M) samples, non-malignant (N) and malignant tumor (T) areas, and prostatectomy tissue samples were cryo-sectioned into extraction tubes. The fraction of epithelial cells in the samples was determined by examination of parallel hematoxylin–eosin-stained sections.
DNA methylation profiling
Genomic DNA was isolated by the AllPrep DNA/RNA/Protein method (Qiagen, Hilden, Germany) and DNA quality and quantity were determined by spectrophotometry (NanoDrop, Thermo Scientific, Wilmington, DE, USA) and the Qubit dsDNA BR assay kit on a Qubit 3.0 Fluorometer (Invitrogen, Carlsbad, CA). DNA (300 ng) was bisulfite converted using the EZ DNA Methylation Kit (Zymo Research, Irvine, USA) and thereafter applied to the Infinium Methylation EPIC arrays (lllumina, San Diego, CA), and operated according to the manufacturer’s instructions. Array analysis including pre-processing and normalization was performed as previously described , with some modifications. The quality of each array was evaluated with built-in controls and the matching identities of the N and T paired samples were confirmed by using the 59 built-in single nucleotide polymorphisms (SNPs) (Additional file 1: Fig. S1).
The fluorescence intensities were extracted using the Methylation Module (1.9.0) in the Genome Studio software (V2011.1), whereas pre-processing and downstream analysis was done using R (v3.4.1). An overview of the pre-processing steps is shown in Additional file 2: Fig. S2. Data was normalized using the BIMQ method to compensate for the two different bead types used in the array . Cross-reactive CpG probes that aligned to multiple loci in the genome or were located in methylation quantitative trait loci (meQTLs) [20, 21] or located less than 5 bp from a known single nucleotide polymorphism in the European population  were excluded. CpG probes with detection P value > 0.05 in any sample were also excluded. The methylation level (β-value) of each CpG site ranging from 0 (no methylation) to 1 (complete methylation) was used as the measure for methylation level in down-stream analyses.
Methylation levels (β-values) were extracted for promoter associated CpGs located in the TSS1500 (the region that covers − 200 to − 1500 nucleotides upstream of Transcription Start Site (TSS)), TSS200 (from TSS to − 200 nucleotides upstream of TSS), and 5′UTR regions, which showed an overall SD > 0.05. Differentially methylated CpG sites (DM-CpGs) were defined as a mean delta-β-value > 0.3 or < − 0.3 between compared groups. Heatmaps were produced in R (v3.4.1) without scaling and using default settings for clustering.
Principal component analysis (PCA), based on centralized β-values were used for unsupervised multivariate projection.
Copy number variation analysis
The raw signal intensity data from the HumanMethylationEPIC arrays was imported to R by the minfi package  and CNV analysis was performed using the conumee package . The 12 non-malignant samples were used as reference samples and the analysis included the AR gene as a detailed region. The CNV status of AR was determined through manual inspection. Human genome GRCh37 (NCBI)/hg19 was used for assigning all chromosome positions.
Whole genome expression analysis
The majority (96%) of the samples applied to DNA methylation analysis had been previously applied to whole genome expression array analysis using the human HT12 Illumina Beadchip technique (Illumina, San Diego, CA) [7, 9]. In 48 cases, total RNA was extracted in parallel with the DNA extractions using the AllPrep DNA/RNA/Protein method. In the other cases, RNA had been extracted from parallel tissue sections using the AllPrep or the Trizol method (Invitrogen, Carlsbad, CA). Bead chip data from two separate gene expression studies (GEO Datasets GSE29650 and GSE101607) were combined for all probes with average signals above two-times the mean background level in at least one sample per study array, leaving 15,232 gene transcripts for subsequent analysis. The gene expression arrays were individually normalized using the quantile method and data was centered by the mean for each probe. Gene transcripts were matched to corresponding promoter CpG sites by Refseq accession and Entrez gene identification numbers.
Assessment of AR activity and proliferation in metastasis samples
The AR activity of a metastasis sample was defined by its relative expression levels of genes predefined to be associated with canonical AR activity; AR, FOXA1, HOXB13, KLK2, KLK3, NKX3-1, STEAP2, and TMPRSS2, as previously described . The AR activity of each metastasis sample was represented by its score on the first PCA score vector (t), capturing the largest variation in the data as a linear combination of the selected transcript levels. Similarly, tumor cell proliferation in metastasis samples was represented by the t scores obtained from PCA of genes predefined to be associated with cell cycle proliferation ; RAD54L, CDC20, CENPF, CDKN3, PBK, TOP2A, LCMT2, ASF1B, KIF20A, CDCA8, NUSAP1, PRC1, PLK1, CDCA3, CEP55, CDC2, KIF11, BUB1B, TK1, ASPM, PTTG1, ORC6L, FOXM1, RAD51, CENPM, CDAN1, KIAA0101, MCM10.
Classification of metastasis samples into metastasis subtypes
The metastasis samples were classified into molecular metastasis subtypes MetA (n = 46), MetB (n = 12), and MetC (n = 8) by PCA and unsupervised clustering based on their transcription profiles, as previously described . Four samples were not possibly to classify in relation to MetA–C, due to lack of transcriptomic data.
Functional enrichment analysis
Functional enrichment analysis was performed by the MetaCore software (GeneGo, Thomson Reuters, New York, NY). In comparisons between different groups (T and N or M and T), analysis was based on genes showing (i) DM-CpGs (mean delta-β-value > 0.3 or < − 0.3 in promoter regions), (ii) inverse correlation between the methylation levels (β values) of CpG sites (mCpGs) and the corresponding transcript levels, and iii) significantly different transcript levels (P < 0.05). In the analysis of promoter methylation levels in relation to AR activity, analysis was based on genes showing (i) positive or negative correlations between mCpGs (β values) and AR activity score (corr > 0.4; − < 0.4), (ii) a considerable difference between samples (stdev > 0.15), and (iii) inverse correlation (> 0.4) between mCpGs (β values) and the corresponding transcript levels. Sets of genes associated with a functional pathway were determined as significantly enriched based on P values representing the probability for a process to arise by chance, considering the numbers of enriched gene products in the data versus the number of genes in the process. P values were adjusted by considering the rank of the process, given the total number of processes in the MetaCore ontology. Enriched pathways were expressed as pathway maps or as process networks, which were created by MetaCore on the basis of both pathway maps and gene ontology processes.
The Wilcoxon rank sum test was used to compare independent groups and the Wilcoxon signed rank sum test for paired data. Bivariate correlations were analyzed according to Pearson. Survival analysis was performed by Kaplan–Meier analysis and by multivariate Cox regression analysis, with death of PC as event and death by other causes as censored events. The false discovery rate (FDR) were controlled according to the Benjamini–Hochberg procedure. To evaluate the consistency of the clusters a silhouette analysis was performed using the cluster package  in R.
DNA promoter methylation associated with prostate cancer development and progression
The methylation levels (β-values) of CpG sites (mCpGs) were analyzed in pairs of non-malignant (N) and malignant primary tumor (T) areas originating from radical prostatectomy samples of 12 PC patients. In parallel, the mCpG profiles were analyzed in 70 bone metastasis (M) samples of separate patients with hormone-naive (HN, n = 14), short-term castrated (ST, n = 4), or castrate-resistant disease (CRPC, n = 52) (Table 1). To facilitate correlation studies between methylation levels of gene promoter regions and corresponding transcript levels, the analysis was focused on CpG sites in promoter regions (defined as CpGs located within the TSS1500, TSS200, 5´UTR regions), leaving 121,944 sites for down-stream analysis after filtration (Additional file 2: Fig. S2).
Principal component analysis based on the 121,944 promoter-associated mCpGs showed a clear separation between the N, T, and M samples, except for two T and two M samples (Fig. 1A). By histological re-examination of tissue sections, the miss-classified samples could not be explained by any obvious morphological characteristics.
The epithelial cell fraction was estimated by microscopy examination of tissue sections. The N samples showed a significantly lower epithelial cell content compared to the T and M samples; median (25th; 75th percentiles) epithelial fraction in N, T, and M samples were 50 (30–50), 80 (70–80), and 70 (50–80)%, respectively (P < 0.001, Fig. 1B). The M samples did not cluster related to epithelial cell fraction (Fig. 1B) or previous treatment (Fig. 1A). Based on a difference in β-value of at least ± 0.3, a total of 15,358 DM-CpGs were found for T, N, and/or M samples (Fig. 1C).
Differential methylation analysis between T and N samples resulted in 4360 DM-CpGs, most of which (86%) corresponding to CpG hypermethylation of the T samples (Fig. 1C). Those DM-CpGs (representing 1759 unique genes) clearly separated T and N samples in a hierarchical cluster analysis (Fig. 1D). The pronounced hypermethylation seen in the tumor samples was associated with significantly reduced gene expression levels in 25% (n = 442) of the DM genes (Additional file 4: Table S1). Functional enrichment analysis indicated altered processes of potential importance for cellular transformation and cancer development, including e.g. regulation of proliferation, angiogenesis and cell adhesion (Additional file 5: Table S2). In contrast, genes showing hypomethylation in combination with significantly induced expression in primary PC (n = 67, 4% of the DM genes) (Additional file 4: Table S1) showed no significant functional enrichment, according to the MetaCore software (Additional file 5: Table S2).
To further focus on DNA methylation alterations associated with progression from a primary tumor in the prostate to metastasis in the bone, differently methylated genes between these stages were identified. Among the 1183 DM-CpGs between T and M, the majority (88%) corresponded to hypomethylation of the M samples (Fig. 1C, E). Notably, only 70 of the hypomethylated genes also showed significantly increased transcript levels and for these no consensus regarding enriched ontologies could be found, while functional enrichment analysis of the few hypermethylated genes with parallel reduction in gene expression (n = 23) suggested altered apoptosis response in metastases compared to primary tumors (Additional file 5: Table S2 and Additional file 6: Table S3).
The fact that most of the DM-CpGs seen between M and T were also DM between M and N, suggests that the overlapping CpGs were likely organ-related rather than associated with bone metastasis (Fig. 1C). Further analysis was therefore focused on the 85 DM-CpGs (80 unique genes) that uniquely diverged between M and T and possibly related specifically to disease progression (Figs. 1C, 2A). Of the 85 DM-CpGs between M and T samples, 94% showed a general hypomethylation in the metastases (Figs. 1C, 2A). Hierarchical cluster analysis identified three main clusters with diversity in the degree of hypermethylation among the M samples (Fig. 2A). These clusters did not correlate to treatment status, nor to epithelial cell fraction (data not shown). Next, an integrated transcriptome analysis that aimed to identifying hyper- or hypomethylated genes with differential gene expression of potential relevance for bone metastasis was performed. Corresponding transcript levels were available for 55 of the 80 unique differently methylated genes. A set of genes showing either hypermethylation in combination with significant gene expression downregulation (SLC8B1) or hypomethylation in combination with gene overexpression (HECW2, AR, GMNN, TSPAN18, and MEFV) in metastases (Fig. 2B) were identified. The hypomethylated ubiquitin ligase HECW2, the DNA replication inhibitor GMNN (geminin), and the tetraspanin TSPAN18 showed clear overexpression in metastases, and may deserve further studies, while fold changes for the innate immunity regulator MEFV and the sodium/calcium exchanger SLC8B1 were less convincing. Of special notice was the hypomethylation of the AR gene in the M samples. Although gene amplification and overexpression of the AR in CRPC have been thoroughly described in the literature, less is known about its epigenetic regulation, and the AR was therefore chosen for further exploration.
Androgen receptor gene amplification, promoter DNA hypomethylation, and mRNA overexpression in metastases
The AR mRNA overexpression observed in the M samples (Fig. 2B) corresponded to two transcript variants (coding for the full-length AR and its shorter AR45 variant, lacking the N-terminal domain). Accordingly, hypomethylation was observed at several AR promoter-associated CpGs in M compared to N and T samples, with the most pronounced CpGs being located close to the TSS200 and 5′-UTR regions (Fig. 3A). In order to examine if AR amplification, in addition to AR promoter hypomethylation, could be associated with the high AR mRNA levels observed, the AR copy numbers of the M samples were assessed and found to be amplified in 64% of the cases. Metastases with AR amplification showed the highest AR mRNA levels, but also M samples with a normal AR copy number had significantly higher expression levels than N and T samples (Fig. 3B). As seen in Fig. 3A, all but one of the untreated M samples had a normal AR copy number, while AR promoter hypomethylation could be observed already at that stage. Interestingly, there was no difference in AR activity in-between M samples with and without AR amplification, as determined by their score on a vector constructed from mRNA levels of genes involved in or regulated through canonical AR signaling (Fig. 3B). This indicates that AR activity within M samples is regulated by other means than by AR copy number, AR promoter methylation, and AR expression levels only.
DNA promoter methylation associated with tumor androgen receptor activity and patient prognosis
To identify mCpGs related to the AR activity within PC metastases, and thus their potential AR dependence, we extracted data for all promoter CpG sites showing positive or negative correlation (correlation ≥ 0.4, StDev > 0.15) to the sample AR activity score. This Methylation Classifier for Androgen receptor activity (MCA) signature consisted of 2970 unique CpGs and separated the M samples into two main clusters named MCA negative and MCA positive (Fig. 4A). The frequency of HN, ST or CRPC samples did not significantly differ between MCA negative and positive metastasis, and neither did the frequency of AR amplified samples or the median epithelial cell fraction (Fig. 4A). Interestingly, the samples of the molecular metastasis subtype MetA were enriched in MCA positive metastasis, while MetB and MetC samples were enriched in MCA negative metastasis (P < 0.001, Chi-square test).
MCA positive metastasis showed hypomethylation and overexpression of genes associated with recurrent gene fusions in PC and AR activation and downstream signaling in PC (Fig. 4A, Table 2), indicating high AR dependence in those cases. In contrast, MCA negative showed hypermethylation and low expression of many genes related to AR activity and instead hypomethylation and high expression of genes associated with epithelial to mesenchymal transition (EMT), cytoskeletal remodeling, and immune response (Fig. 4A, Table 2). Accordingly, MCA negative metastases had significantly lower AR activity scores and a tendency of higher proliferation scores than MCA positive metastases (Fig. 4B, C). In addition, MCA negative patients showed significantly lower serum PSA levels and a worse prognosis after ADT than MCA positive patients (Fig. 4D, E). The median cancer-specific survival after ADT was 30 months for MCA negative patients compared to 60 months for MCA positive patients (P = 0.001, Fig. 4E). Analysis of cluster consistency showed an average silhouette width of 0.23 in the MCA positive and negative clusters (Additional file 3: Fig. S3).
In Cox regression analysis, the MCA negative signature, proliferation score, and the MetB subtype were significantly associated with poor survival after ADT (Table 3). Multivariate Cox analysis indicated that the MCA signature provided prognostic information that was independent from the proliferation and subtype status, with an increased risk of 2.7 for patients being MCA negative (P = 0.0017) (Table 3). No significant associations were observed between patient prognosis and AR activity score, AR amplification status, serum PSA, or age (Table 3).
This study integrates genome-wide promoter methylation data with gene expression data and describes consistent changes occurring during PC disease progression from non-malignant prostate epithelium to primary PC and further to bone metastatic disease. In metastatic samples, the study also specifically explores if promoter methylation levels are related to the sample AR activity, defined from transcript levels of genes recently described to differentiate clinically relevant molecular subtypes of PC bone metastases based on AR activity [9, 10, 12].
By pairwise comparisons of 12 tumor areas and adjacent non-malignant tissue isolated from radical prostatectomies, a pronounced hypermethylation was observed in primary PC. This is in line with a wealth of data, summarized in , suggesting not only that hypermethylation is an early event in PC tumorigenesis but also that it possibly could be used for diagnostic purpose. We observed hypermethylation of GSTP1, but also of other genes well known to be hypermethylated in PC, including AOX1, APC, BARHL2, CCDC8, CDKN2A, CYP27A1, EFS, GRASP, HOXA3, HOXC11, HOXD3, KIT, NXK2-1, NXK2-5, PHOX2A, POU3F3, PTGS2, RARB, RHCG, SIX6, TBX15, TMEM106A, WNT2, and ZNF154  (Additional file 4: Table S1). Interestingly, methylation changes in APC, HOXD3, and PTGS2 have previously been suggested to provide prognostic information [27,28,29]. The low number of primary tumor samples analyzed in the current study did, however, not allow evaluation of DM-CpGs in relation to prognosis. Instead, analysis was focused on identifying general effects of methylation on changes in gene expression during PC development. Functional pathway analysis using Gene Set Enrichment Analysis indicated hypermethylation and downregulation of genes associated with muscle contraction, such as ACTA, CNN, DES, TAGLN, and TPN, normally expressed in smooth muscle cells , but also of genes involved in diverse developmental (e.g. ANGPT1, STAT3, STAT5, TEK) and cytoskeletal (e.g. ANXA2, ANXA6, CD44, SVIL) processes, as well as in regulation of proliferation (e.g. FGF2, FGFR1, PIK3R1, PRKCA) and cell adhesion (e.g. CDH5, CDH23, PCDH18, STMN2). While reduced cell adhesion is an important hallmark for carcinomas, many of the other functional effects suggested from promoter hypermethylation in PC may be biased not only by the more abundant stroma in the T samples, but also by its different cellular content. While the stroma in the normal prostate is mainly composed of smooth muscle cells, prostate tumors show a progressive loss of smooth muscle cells in the favor of cancer-associated fibroblasts . Hypomethylation in primary PC was uncommon and not associated with specific functional processes, at least based on the literature summarized by MetaCore ontology software.
Most of the identified DM-CpGs observed between primary PC and non-malignant prostate tissue were present also in the metastatic tissue. In addition, a unique set of DM-CpGs were observed when going from primary PC to metastases. The large majority of these were hypomethylated and included sites in the AR promoter region. A detailed exploration showed a trend of demethylation at several CpGs in the AR promoter, including two alternative transcription start sites that correlated well with induced levels of both AR transcript 1 and 2 in the metastatic samples. High AR mRNA levels in CRPC have been previously associated with AR amplification [32,33,34,35,36], while methylation as a regulator of AR expression is less well described and studies have also shown inconsistent results [16, 37,38,39,40]. Our results are in line with recent findings by Zaho and co-workers demonstrating AR promoter hypomethylation in CRPC , but also add important knowledge by showing hypomethylation in previously un-treated metastasis samples suggesting that AR demethylation may be an earlier cause to AR induction in PC metastases than AR amplification, as AR amplification is generally not observed prior to ADT. Importantly, neither the large variance in AR promoter methylation levels observed among the metastatic samples nor their AR amplification status showed any clear relationship to the tumor sample AR activity.
In an attempt to study universal promoter methylation in relation to AR activity within individual metastasis samples, each sample was given an AR activity score related to its expression levels of genes involved in or regulated by canonical AR signaling, and mCpGs positively or negatively correlated to this score were taken further for functional exploration. A Methylation Classifier for Androgen receptor activity (MCA) signature was identified and used to classify the bone metastasis samples. Taken together, our results suggested that low promoter methylation levels of certain AR-regulated genes (i.e. SLC45A3, STEAP2, ELK4 and TMPRSS2) may contribute to the AR-driven tumor phenotype seen in the majority of CRPC patients with metastatic disease [2,3,4,5,6,7,8,9,10,11,12], while high promoter methylation of the same genes may contribute to the development of less AR-dependent CRPC, commonly referred to as AR-indifferent, small-cell or neuroendocrine-like CRPC . It is of importance to note, though, that histological evaluation did not confirm any enrichment of a neuroendocrine-like phenotype in samples of the MCA negative compared to the positive cluster. This was concluded from reviewing the results from immunohistochemical analysis of the neuroendocrine marker chromogranin A, previously assessed in 60 of the 66 MCA-classified samples  (data not shown). Neither could the MCA class be predicted from immunohistochemical analysis of the AR protein level  (n = 63, data not shown).
Samples of the metastasis subtype MetA previously described by us as AR-driven  primarily showed hypomethylation and high expression of canonically AR-regulated genes, while samples of the less AR-driven subtypes MetB and MetC were enriched among samples with an MCA negative signature (showing hypermethylation and low expression of many AR targeting genes). Indirectly, our result thus implies promoter methylation as a possible mechanism behind the low AR activity defining metastasis subtypes MetB and MetC. Accordingly, patients with MCA positive signature had higher serum PSA levels and a better prognosis than the smaller group of MCA negative patients, probably due to more differentiated, less proliferative tumor cells being more AR dependent and thereby showing a better response to AR targeting therapies, similar to what previously have been described for MetA in comparison to MetB–C . Notably, the MCA signature provided strong information related to patient prognosis after ADT, with MCA positive patients showing 2-times longer median survival than MCA negative patients, resulting in a 5-year cancer-specific survival of 50 and 20% respectively. Despite the obvious clinical relevance of separating metastasis samples into MCA positive or negative cases, the clusters showed variable consistencies according to the silhouette analysis. This might originate from the procedure of selecting all mCpG sites for analysis that correlated to the tumor AR activity on a continuous scale, i.e. the classifier was not originally built to dichotomize samples into two groups.
As MCA negative metastases with low AR activity scores were specified by hypermethylation of AR-regulated genes, it is tempting to speculate that they could be sensitized for AR-directed therapies by treatment with demethylating agents, such as have been shown possible by inhibitors to DNA methyl transferase 1 in different experimental systems for CRPC [42, 43]. In a recent paper, Xiao and co-workers furthermore identified the histone lysine-N-methyltransferase EZH2 as a possible epigenetic regulator of tumor response to AR-targeting therapy and demonstrated possibilities with its inhibition in the treatment of CRPC .
While the methylation levels of a relatively small set of genes were directly associated with the AR activity observed in PC bone metastases, the methylation level of a much larger gene set seemed involved in regulating diverse processes such as EMT, cytoskeletal remodeling, and immune responses, all standing in an inverse correlation to tumor AR activity.
In conclusion, this study describes general patterns of gene promoter methylation during PC disease progression, and specifically identifies a novel promoter MCA signature that is associated with canonical AR activity in PC bone metastases and to patient prognosis after ADT. Based on this methylation signature, patients with PC bone metastases could be stratified into one of two patient groups; MCA positive patients being likely to respond to ADT and possibly also other AR targeting therapies or MCA negative patients being less responsive to AR targeting therapies, but possibly suitable for treatment with AR inhibitors in combination with epigenetic modulators (currently in clinical trials, https://clinicaltrials.gov). The next step will be to verify the suggested therapy-predictive value of the MCA signature. This could be retrospectively done by analyzing methylation signatures in circulating tumor DNA collected from patients prior to treatment for metastatic PC, similarly to what is described in a recent paper by Wu et al. . For clinical implementation, however, a stable non-cluster-based classifier need to be developed. The overall cause to the two methylation types observed among PC bone metastases also remains to be identified, as it could give important implications for future developments of novel treatment strategies.
Availability of data and materials
The dataset generated (GSE174613) and analyzed (GSE29650 and GSE101607) during the current study are available at https://www.ncbi.nlm.nih.gov/geo/.
Copy number variation
Castration-resistant prostate cancer
Epithelial to mesenchymal transition
Methylation classifier for androgen receptor activity
Methylation levels (β values) of CpG sites
Metastasis subtype A–C
Methylation quantitative trait loci
Non-malignant prostate tumor tissue
Principal component analysis
Prostate specific antigen
Single nucleotide polymorphism
Malignant prostate tumor tissue
Transcription start site
Cornford P, van den Bergh RCN, Briers E, Van den Broeck T, Cumberbatch MG, De Santis M, et al. EAU-EANM-ESTRO-ESUR-SIOG guidelines on prostate cancer. Part II-2020 update: treatment of relapsing and metastatic prostate cancer. Eur Urol. 2020;20:30773–9.
Robinson D, Van Allen EM, Wu YM, Schultz N, Lonigro RJ, Mosquera JM, et al. Integrative clinical genomics of advanced prostate cancer. Cell. 2015;161:1215–28.
Quigley DA, Dang HX, Zhao SG, Lloyd P, Aggarwal R, Alumkal JJ, et al. Genomic hallmarks and structural variation in metastatic prostate cancer. Cell. 2018;175:889.
De Laere B, Oeyen S, Mayrhofer M, Whitington T, van Dam PJ, Van Oyen P, et al. TP53 outperforms other androgen receptor biomarkers to predict abiraterone or enzalutamide outcome in metastatic castration-resistant prostate cancer. Clin Cancer Res. 2019;25:1766–73.
Abida W, Cyrta J, Heller G, Prandi D, Armenia J, Coleman I, et al. Genomic correlates of clinical outcome in advanced prostate cancer. Proc Natl Acad Sci U S A. 2019;116:11428.
De Laere B, van Dam PJ, Whitington T, Mayrhofer M, Diaz EH, Van den Eynden G, et al. Comprehensive profiling of the androgen receptor in liquid biopsies from castration-resistant prostate cancer reveals novel intra-AR structural variation and splice variant expression patterns. Eur Urol. 2017;72:192–200.
Hörnberg E, Ylitalo EB, Crnalic S, Antti H, Stattin P, Widmark A, Bergh A, Wikström P. Expression of androgen receptor splice variants in prostate cancer bone metastases is associated with castration-resistance and short survival. PLoS ONE. 2011;6:e19059.
Jernberg E, Thysell E, Bovinder Ylitalo E, Rudolfsson S, Crnalic S, Widmark A, Bergh A, Wikström P. Characterization of prostate cancer bone metastases according to expression levels of steroidogenic enzymes and androgen receptor splice variants. PLoS ONE. 2013;8:e77407.
Ylitalo EB, Thysell E, Jernberg E, Lundholm M, Crnalic S, Egevad L, et al. Subgroups of castration-resistant prostate cancer bone metastases defined through an inverse relationship between androgen receptor activity and immune response. Eur Urol. 2017;71:776–87.
Nordstrand A, Bovinder Ylitalo E, Thysell E, Jernberg E, Crnalic S, Widmark A, et al. Bone cell activity in clinical prostate cancer bone metastasis and its inverse relation to tumor cell androgen receptor activity. Int J Mol Sci. 2018;19:1223.
Iglesias-Gato D, Thysell E, Tyanova S, Crnalic S, Santos A, Lima TS, et al. The proteome of prostate cancer bone metastasis reveals heterogeneity with prognostic implications. Clin Cancer Res. 2018;24:5433–44.
Thysell E, Vidman L, Ylitalo EB, Jernberg E, Crnalic S, Iglesias-Gato D, et al. Gene expression profiles define molecular subtypes of prostate cancer bone metastases with different outcomes and morphology traceable back to the primary tumor. Mol Oncol. 2019;13:1763.
You S, Knudsen BS, Erho N, Alshalalfa M, Takhar M, Al-Deen Ashab H, et al. Integrated classification of prostate cancer reveals a novel luminal subtype with poor outcome. Cancer Res. 2016;76:4948.
Zhao SG, Chang SL, Erho N, Yu M, Lehrer J, Alshalalfa M, et al. Associations of luminal and basal subtyping of prostate cancer with prognosis and response to androgen deprivation therapy. JAMA Oncol. 2017;3:1663.
Massie CE, Mills IG, Lynch AG. The importance of DNA methylation in prostate cancer development. J Steroid Biochem Mol Biol. 2017;166:1–15.
Zhao SG, Chen WS, Li H, Foye A, Zhang M, Sjöström M, et al. The DNA methylation landscape of advanced prostate cancer. Nat Genet. 2020;52:778–89.
Pomerantz MM, Qiu X, Zhu Y, Takeda DY, Pan W, Baca SC, et al. Prostate cancer reactivates developmental epigenomic programs during metastatic progression. Nat Genet. 2020;52:790–9.
Borssén M, Haider Z, Landfors M, Norén-Nyström U, Schmiegelow K, Åsberg AE, et al. DNA methylation adds prognostic value to minimal residual disease status in pediatric T-cell acute lymphoblastic leukemia. Pediatr Blood Cancer. 2016;63:1185–92.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.
Gaunt TR, Shihab HA, Hemani G, Min JL, Woodward G, Lyttleton O, et al. Systematic identification of genetic influences on methylation across the human life course. Genome Biol. 2016;17:61.
McClay JL, Shabalin AA, Dozmorov MG, Adkins DE, Kumar G, Nerella S, Clark SL, Bergen SE, Swedish Schizophrenia Consortium, Hultman CM, Magnusson PK, Sullivan PF, Aberg KA, van den Oord EJ. High density methylation QTL analysis in human blood via next-generation sequencing of the methylated genomic DNA fraction. Genome Biol. 2015;16:291.
Zhou W, Laird PW, Shen H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 2017;45(4):22.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.
Hovestadt V, Zapatka M. conumee: enhanced copy-number variation analysis using Illumina DNA methylation arrays. R package version 1.9.0, http://bioconductor.org/packages/conumee/.
Cuzick J, Swanson GP, Fisher G, Brothman AR, Berney DM, Reid JE, et al. Prognostic value of an RNA expression signature derived from cell cycle proliferation genes in patients with prostate cancer: a retrospective study. Lancet Oncol. 2011;12:245.
Maechler M, Rousseeuw P, Struyf A, Hubert M, Hornik K. cluster: cluster analysis basics and extensions. R package version 2.1.0. 2019 (2021).
Yegnasubramanian S, Kowalski J, Gonzalgo ML, Zahurak M, Piantadosi S, Walsh PC, et al. Hypermethylation of CpG islands in primary and metastatic human prostate cancer. Cancer Res. 2004;64:1975–86.
Henrique R, Ribeiro FR, Fonseca D, Hoque MO, Carvalho AL, Costa VL, et al. High promoter methylation levels of APC predict poor prognosis in sextant biopsies from prostate cancer patients. Clin Cancer Res. 2007;13:6122–9.
Stott-Miller M, Zhao S, Wright JL, Kolb S, Bibikova M, Klotzle B, et al. Validation study of genes with hypermethylated promoter regions associated with prostate cancer recurrence. Cancer Epidemiol Biomark Prev. 2014;23:1331–9.
Lee MY, Park C, Berent RM, Park PJ, Fuchs R, Syn H, et al. Smooth muscle cell genome browser: enabling the identification of novel serum response factor target genes. PLoS ONE. 2015;10(8):e0133751.
Cunha GR, Hayward SW, Wang YZ, Ricke WA. Role of the stromal microenvironment in carcinogenesis of the prostate. Int J Cancer. 2003;107:1–10.
Visakorpi T, Hyytinen E, Koivisto P, Tanner M, Keinanen R, Palmberg C, et al. In vivo amplification of the androgen receptor gene and progression of human prostate cancer. Nat Genet. 1995;9:401–6.
Palmberg C, Koivisto P, Hyytinen E, Isola J, Visakorpi T, Kallioniemi OP, Tammela T. Androgen receptor gene amplification in a recurrent prostate cancer after monotherapy with the nonsteroidal potent antiandrogen Casodex (bicalutamide) with a subsequent favorable response to maximal androgen blockade. Eur Urol. 1997;31:216–9.
Bubendorf L, Kononen J, Koivisto P, Schraml P, Moch H, Gasser TC, et al. Survey of gene amplifications during prostate cancer progression by high-throughout fluorescence in situ hybridization on tissue microarrays. Cancer Res. 1999;59:803–6.
Kaltz-Wittmer C, Klenk U, Glaessgen A, Aust DE, Diebold J, Lohrs U, Baretton GB. FISH analysis of gene aberrations (MYC, CCND1, ERBB2, RB, and AR) in advanced prostatic carcinomas before and after androgen deprivation therapy. Lab Investig J Tech Methods Pathol. 2000;80:1455–64.
Djusberg E, Jernberg E, Thysell E, Golovleva I, Lundberg P, Crnalic S, et al. High levels of the AR-V7 splice variant and co-amplification of the golgi protein coding YIPF6 in AR amplified prostate cancer bone metastases. Prostate. 2017;77:625–38.
Kinoshita H, Shi Y, Sandefur C, Meisner LF, Chang C, Choon A, et al. Methylation of the androgen receptor minimal promoter silences transcription in human prostate cancer. Cancer Res. 2000;60:3623–30.
Nakayama T, Watanabe M, Suzuki H, Toyota M, Sekita N, Hirokawa Y, et al. Epigenetic regulation of androgen receptor gene expression in human prostate cancers. Lab Investig. 2000;80:1789–96.
Friedlander TW, Roy R, Tomlins SA, Ngo VT, Kobayashi Y, Azameera A, et al. Common structural and epigenetic changes in the genome of castration-resistant prostate cancer. Cancer Res. 2012;72:616–25.
Jarrard DF, Kinoshita H, Shi Y, Sandefur C, Hoff D, Meisner LF, et al. Methylation of the androgen receptor promoter CpG island is associated with loss of androgen receptor expression in prostate cancer cells. Cancer Res. 1998;58:5310–4.
Beltran H, Prandi D, Mosquera JM, Benelli M, Puca L, Cyrta J, et al. Divergent clonal evolution of castration-resistant neuroendocrine prostate cancer. Nat Med. 2016;22:298–305.
Gravina GL, Festuccia C, Millimaggi D, Dolo V, Tombolini V, de Vito M, Vicentini C, Bologna M. Chronic azacitidine treatment results in differentiating effects, sensitizes against bicalutamide in androgen-independent prostate cancer cells. Prostate. 2008;68:793–801.
Gravina GL, Marampon F, Di Staso M, Bonfili P, Vitturini A, Jannini EA, et al. 5-Azacitidine restores and amplifies the bicalutamide response on preclinical models of androgen receptor expressing or deficient prostate tumors. Prostate. 2010;70:1166–78.
Xiao L, Tien JC, Vo J, Tan M, Parolia A, Zhang Y, et al. Epigenetic reprogramming with antisense oligonucleotides enhances the effectiveness of androgen receptor inhibition in castration-resistant prostate cancer. Cancer Res. 2018;78:5731–40.
Wu A, Cremaschi P, Wetterskog D, Conteduca V, Franceschini GM, Kleftogiannis D, et al. Genome-wide plasma DNA methylation features of metastatic prostate cancer. J Clin Investig. 2020;130:1991–2000.
Pernilla Andersson and Susanne Gidlund are acknowledged for superior technical assistance.
Open access funding provided by Umea University. The study was supported by grants from the Cancerfonden, Kempestiftelserna, Swedish Research Council, the Swedish Cancer Society, Swedish Foundation for Strategic Research, Cancer Research Foundation in Northern Sweden, the Kempe research foundation, The Swedish Prostate Cancer Federation, the Medical Faculty of Umeå University, and Uppsala-Umeå Comprehensive Cancer Consortium. Financial support was also provided through a regional agreement between Umeå University and Västerbotten County Council on cooperation in the field of Medicine, Odontology and Health.
Ethics approval and consent to participate
Patients gave their informed consent and the study was conducted in accordance with the Declaration of Helsinki. The study was approved by the local ethic review board of Umeå University (Dnr 03-158, Dnr 04-26M, Dnr 03-482, Dnr 2013-57-31M, Dnr 2017-197-32M).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Analysis of 59 built-in SNP on the HumanMethylation EPIC array to confirm identity of multiple samples taken from the same individual (N and T tissue).
Schematic flowchart of the pre-processing steps of the HumanMethylation EPIC arrays.
Silhouette analysis showing the cluster consistencies of the MCA positive and negative clusters in Figure 4A.
Table S1. Differetially methylated CpGs in promoter regions of genes in primary prostate tumor tissue (T) compared to adjacent non-malignant prostate tissue (N), their correlations to corresponding gene expression levels, and the fold change in gene expression levels between sample groups.
Table S2. Top 10 enriched process networks in localized prostate tumor tissue (T) compared to adjacent non-malignant prostate tissue (N) or in prostate cancer bone metastases (M) compared to localized tumors.
Table S3. Differetially methylated CpGs in promoter regions of genes in bone metastases (M) compared to localized prostate tumors (T), their correlations to corresponding gene expression levels, and the fold change in gene expression levels between sample groups.
Table S4. CpG sites of the AR gene shown in Figure 3A.
About this article
Cite this article
Ylitalo, E.B., Thysell, E., Landfors, M. et al. A novel DNA methylation signature is associated with androgen receptor activity and patient prognosis in bone metastatic prostate cancer. Clin Epigenet 13, 133 (2021). https://doi.org/10.1186/s13148-021-01119-0
- Prostate cancer
- DNA methylation
- Gene expression
- Androgen receptor