- Open Access
Epigenetic silencing of MEIS2 in prostate cancer recurrence
Clinical Epigenetics volume 11, Article number: 147 (2019)
Current diagnostic and prognostic tools for prostate cancer (PC) are suboptimal, resulting in overdiagnosis and overtreatment of clinically insignificant tumors. Thus, to improve the management of PC, novel biomarkers are urgently needed.
In this study, we integrated genome-wide methylome (Illumina 450K DNA methylation array (450K)) and RNA sequencing (RNAseq) data performed in a discovery set of 27 PC and 15 adjacent normal (AN) prostate tissue samples to identify candidate driver genes involved in PC development and/or progression. We found significant enrichment for homeobox genes among the most aberrantly methylated and transcriptionally dysregulated genes in PC. Specifically, homeobox gene MEIS2 (Myeloid Ecotropic viral Insertion Site 2) was significantly hypermethylated (p < 0.0001, Mann-Whitney test) and transcriptionally downregulated (p < 0.0001, Mann-Whitney test) in PC compared to non-malignant prostate tissue in our discovery sample set, which was also confirmed in an independent validation set including > 500 PC and AN tissue samples in total (TCGA cohort analyzed by 450K and RNAseq). Furthermore, in three independent radical prostatectomy (RP) cohorts (n > 700 patients in total), low MEIS2 transcriptional expression was significantly associated with poor biochemical recurrence (BCR) free survival (p = 0.0084, 0.0001, and 0.0191, respectively; log-rank test). Next, we analyzed another RP cohort consisting of > 200 PC, AN, and benign prostatic hyperplasia (BPH) samples by quantitative methylation-specific PCR (qMSP) and found that MEIS2 was significantly hypermethylated (p < 0.0001, Mann-Whitney test) in PC compared to non-malignant prostate tissue samples (AN and BPH) with an AUC > 0.84. Moreover, in this cohort, aberrant MEIS2 hypermethylation was significantly associated with post-operative BCR (p = 0.0068, log-rank test), which was subsequently confirmed (p = 0.0067; log-rank test) in the independent TCGA validation cohort (497 RP patients; 450K data).
To the best of our knowledge, this is the first study to investigate, demonstrate, and independently validate a prognostic biomarker potential for MEIS2 at the transcriptional expression level and at the DNA methylation level in PC.
Prostate cancer (PC) is the most common non-cutaneous cancer among men in the Western world . Early detection of PC is critical, as localized PC is curable by radiation therapy or radical prostatectomy (RP), whereas metastatic disease is lethal with only palliative treatments available. Furthermore, accurate risk stratification at diagnosis is important for treatment decisions, as some PCs remain latent throughout the lifetime of the patients while others may progress to aggressive metastatic disease. However, the currently available diagnostic and prognostic tools for PC are suboptimal, and novel biomarkers are urgently needed to reduce overdiagnosis and overtreatment of clinically insignificant PCs .
A molecular hallmark for PC is aberrant DNA methylation of CpG island-containing gene promoters . Promoter hypermethylation is associated with transcriptional repression of, e.g., tumor suppressor genes, whereas promoter hypomethylation is linked with activation of, e.g., oncogenes . Thus, aberrant DNA methylation may play an important role in driving PC oncogenesis and/or progression when affecting the corresponding transcript expression. Several epigenetic candidate diagnostic markers for PC have been discovered through comparison of DNA methylation alterations in PC and non-malignant prostate tissue samples [5,6,7,8,9,10, 12]. Some of these methylation marker candidates have also shown prognostic potential for prediction of time to biochemical recurrence (BCR) [5, 6, 9, 11, 12]. In the present study, we investigated the diagnostic and prognostic biomarker potential of MEIS2 (Myeloid Ecotropic viral Insertion Site 2) in relation to PC. MEIS2 was selected as candidate gene because we found it to be significantly hypermethylated and downregulated in PC compared to non-malignant prostate tissue samples analyzed by genome-wide methylome and transcriptome profiling (see below).
MEIS2 is a homeobox gene and part of the TALE (three amino acid loop extensions) family of proteins. TALE proteins are a subtype of homeobox proteins that bind to HOX proteins and specify their transcriptional activity . MEIS2 can bind to the domain of HOXB13 that includes the G84E mutation, which in turn has been associated with elevated PC risk and aggressiveness, although it remains unclear how this mutation may affect the HOXB13-MEIS2 interaction and possibly promote the initiation and/or progression of PC [13,14,15]. Moreover, previous PC studies have shown that MEIS2 transcriptional and protein expression levels decrease gradually from non-malignant prostate to primary PC and to metastatic PC tissue samples, suggesting that MEIS2 plays a tumor suppressive role and may be involved in PC development and/or progression [16, 17]. Similarly, downregulated MEIS2 protein expression has been associated with poor overall survival in a small cohort of 83 PC patients . Another study has suggested MEIS2 as an important component of a signaling circuit with IκBα/NF-κB (p65), miR-196b-3p, and PPP3CC (protein phosphatase 3 catalytic subunit gamma) that is involved in progression to castration resistant PC (CRPC) . Together, these previous studies indicate a prognostic potential of MEIS2 in PC, but until now, only protein expression of MEIS2 has been investigated for its prognostic potential in two cohorts [16, 17].
In the present study, in an effort to identify novel candidate biomarkers for PC, we combined DNA methylation (Illumina 450K DNA methylation array (450K)) and matched RNA expression (RNA sequencing (RNAseq)) data from a set of PC and adjacent normal (AN) tissue samples from 29 patients treated by RP. Among the most differentially methylated and differentially expressed genes, we found a significant overrepresentation of homeobox and homeodomain-containing genes, including MEIS2 that was hypermethylated and downregulated in PC, indicating epigenetic silencing. Furthermore, low transcriptional expression and DNA hypermethylation of MEIS2 was significantly associated with BCR after RP in multiple large independent RP cohorts including > 700 PC patients in total. This is the first study to investigate, demonstrate, and independently validate a prognostic biomarker potential for MEIS2 transcriptional expression and DNA methylation in PC.
Materials and methods
Patient samples used for RNAseq and 450K methylation profiling (discovery)
For the discovery set, radical prostatectomy (RP) tissue specimens from 29 PC patients (Additional file 1) treated at the Department of Urology between May 2003 and October 2012 (Aarhus University Hospital, Denmark) were obtained immediately after surgery and stored at − 80 °C (fresh frozen in TissueTek). Hematoxylin and eosin (HE) stained prostate tissue sections were evaluated by a trained pathologist and adjacent normal (AN) and PC (tumor) areas were marked for laser capture microdissection (LCM, VeritasTM 704 (Arcturus)). For each sample, after LCM of the top 15–25 6-μm sections, total RNA (> 200 bp) was extracted using the RNeasy micro Kit (Qiagen) according to the manufacturer’s instructions. RNA concentration and RNA quality was assessed using the Agilent RNA 6000 Pico Chip on an Agilent 2100 Bioanalyzer (RIN ≥ 6). Similarly, after LCM of the next 15–25 6-μm sections, genomic DNA was extracted using the Puregene system (Qiagen) according to the manual provided by the manufacturer. DNA concentrations were assessed using the Quant-iT PicoGreen dsDNA Assay Kit (Life Technologies).
RNA sequencing (RNAseq) and Illumina 450K DNA methylation array (450K)
For RNAseq, directional indexed libraries were generated from 10–500 ng total RNA (> 200 bp) from 42 samples (29 PC samples and 13 AN samples) using the Scriptseq™ Complete Gold Kit Version II (Illumina). RNAseq libraries were sequenced on the Illumina HiSeq2000 (15–25 million reads/sample, 2 × 150 bp), and reads were mapped to the human genome (hg19) using the Tuxedo Suite . Counts were calculated using HTSeq . RNAseq data was analyzed in R version 3.1.2 with the EdgeR package version 3.8.5  using counts as input. A total of seven RNA samples (6 PC and 1 AN) had a poor library profile and were removed prior to the final data analysis. Moreover, one tumor sample had a low PC cell content and was also removed prior to analysis.
Genomic DNA from 43 samples (28 tumor samples and 15 AN samples) was bisulfite converted using the EpiTect Bisulfite Kit (Qiagen) and applied to the Illumina 450K DNA methylation array (450K) by service provider Aros Applied Biotechnology A/S (Aarhus, Denmark). Raw 450K DNA methylation array data was analyzed in R version 3.1.2, using the Chip Analysis Methylation Pipeline (ChAMP) package version 1.4.0 . Throughout the analyses, the argument filterXY was set to FALSE to ensure that probes from X and Y chromosomes were not removed. DNA methylation was reported as β-values (range 0-1; 0, unmethylated; 1, completely methylated). One tumor sample (also analyzed by RNAseq) had a low PC cell content and was removed prior to data analysis.
For candidate biomarker discovery, 450K DNA methylation and RNAseq datasets were merged by ENSEMBL gene name annotations, and a combined p value was calculated using Fisher’s method . Among the genes with a significant Benjamini-Hochberg (BH)-adjusted Fishers p value, The Functional Annotation Clustering tool from the Database for Annotation, Visualization and Integrated Discovery (DAVID) [24, 25] was used to identify enriched annotation terms (GO, KEGG, etc.). Furthermore, Spearman’s correlations between DNA methylation and RNA expression levels were calculated for each CpG site. DAVID analyses were also performed on the top 3000 differentially methylated genes (BH adj. p < 0.05) and for the top 2314 significantly differentially expressed genes (BH adj. p < 0.05), respectively.
Patient samples used for quantitative methylation specific PCR analyses (validation)
For validation, we used quantitative methylation-specific PCR (qMSP) to analyze a RP cohort of 264 patients treated for histologically verified clinically localized PC at Department of Urology, Aarhus University Hospital, Denmark, from 1999 to 2013. In all cases, a trained pathologist evaluated formalin-fixed paraffin-embedded (FFPE) archived prostatectomy specimens, and 1.5-mm punch biopsies were taken from representative regions with cancer (n = 254) or AN tissue (n = 37) and used for extraction of genomic DNA (see below) [5, 6, 8, 12]. As non-PC controls, FFPE transurethral resection of the prostate (TURP) tissue samples from benign prostatic hyperplasia (BPH, n = 9) patients were also included. A total of 82 samples were excluded due to insufficient DNA quality (see below). The final analysis included 195 PC, 17 AN, and 6 BPH samples (Table 1).
Quantitative methylation specific PCR (qMSP)
DNA was extracted from FFPE punch biopsies of RP specimens (PC and AN) and TURP specimens (BPH) with the gDNA Eliminator columns from the RNeasy plus micro kit (Qiagen) and bisulfite converted using the EZ-96 DNA Methylation-Gold Kit™ (Zymo Research), as previously described in detail [5, 6, 11]. For qMSP assay design, Primer3 [26, 27] and Beacon DesignerTM (Premier Biosoft) were used. Primer and probe sequences are given in Additional file 2. MEIS2 assay 1 targeted an intronic region of the MEIS2 gene (intron between exon 4 and 5, Fig. 3a), which overlapped three significant differentially methylated probes on the 450K array (cg06933370, cg23677243, and cg26708220, Fig. 3a). MEIS2 assay 2 targeted the promoter region of MEIS2 and overlapped one probe from the 450K array (cg25381383) (Fig. 3a). All qMSP reactions were run in triplicates (10 μL) with 5 ng bisulfite-converted DNA, 6 pmol of each primer, 2 pmol probe, and 5 μL Taqman universal Mastermix no UNG (Applied Biosystems). As controls, standard curves on serially diluted methylated DNA, bisulfite-converted CpGenome Universal Methylated DNA (Millipore), and two negative controls (H2O and whole-genome amplified (WGA) DNA) were included on each plate. For quality/quantity control, aluC4 and MYOD1 assays were used . AluC4 was used for normalization. Reactions were run in 384-well plates on the ViiA7 Real-Time PCR system (Applied Biosystems): 2 min at 50 °C, 10 min at 95 °C, and 40 cycles of 15 sec at 95 °C and 1 min at 56 °C. Quantities for MEIS2 assay 1 and 2, MYOD1, and aluC4 were estimated from the standard curves using QuantStudio™ Real-Time PCR Software (Applied Biosystems). Outliers (more than 2 ct values lower/higher than the ct value of the other replicates) and samples with MYOD1 ct > 38.0 in ≥ 2 of 3 replicate reactions were removed. For MEIS2 assay 2, replicates exceeding ct 38 were set to 0 (WGA cutoff). Samples were considered negative for methylation, if ≥ 2 methylation-specific reactions did not amplify.
Public cohorts (external validation)
RNAseq and clinical data for formalin-fixed paraffin-embedded (FFPE) PC tissue from 106 RP patients from Long et al.  was downloaded from GEO (GSE54460).
From The Cancer Genome Atlas (TCGA, http://cancergenome.nih.gov/), we downloaded RNA sequencing (RNAseq), 450K DNA methylation, and clinical data for 497 RP (PC) and 52 matched AN fresh frozen tissue samples . 450K data was peak corrected  and RNAseq data was mapped to hg19 and processed as previously described . DNA methylation was reported as β-values and RNAseq gene expression as counts per million (CPM). For external validation of our qMSP data, the average methylation level of the 450K probes cg06933370, cg23677243, and cg26708220 was calculated for each patient to mimic qMSP assay 1. qMSP assay 2 covered one probe from the 450K array, cg25381383, which was used for external validation of assay 2.
Normalized microarray (Affymetrix Human Exon 1.0 ST array) RNA expression and clinical data for fresh frozen PC tissue from 126 RP samples were downloaded from GEO (GSE21034) .
All statistical analyses were conducted in STATA version 13.1. To investigate the diagnostic and prognostic potential of MEIS2 DNA methylation and transcriptional expression, we used Mann-Whitney tests, ROC curve analyses, uni- and multivariate Cox regression analyses, Kaplan-Meier, and log-rank tests. In time-to-event analysis, PSA recurrence (cutoff ≥ 0.2 ng/mL) was used as endpoint. In all datasets, dichotomization of patients into high- and low-risk groups based on either MEIS2 expression or MEIS2 methylation levels was made by ROC curve analyses of BCR status at 36 months of follow-up. For the qMSP assays, the cutoffs for dichotomization were 0.109/0.162, respectively. For the MEIS2 methylation model, patients were included in the low methylation group if both qMSP assay 1 and qMSP assay 2 showed low methylation.
Methylome and transcriptome profiling of PC tissue for biomarker discovery
To identify novel epigenetically and transcriptionally deregulated biomarker candidates, we performed RNAseq and 450K DNA methylation analysis on microdissected PC and AN tissue samples from RP specimens. After quality control (see “Materials and methods” section), the final 450K dataset consisted of 27 PC and 15 AN samples and the final RNAseq dataset of 22 PC and 12 AN samples (Fig. 1). A total of 22 PC and 12 AN samples were analyzed by both methods.
The 450K analysis identified 119,519 CpG sites with a significant Benjamini-Hochberg (BH)-adjusted p value when comparing methylation in PC and AN samples (BH adj. p < 0.05, t statistics in ChAMP-package), corresponding to differential methylation at approx. 25% of all CpG sites analyzed. In a multi-dimensional scaling (MDS) plot of the 1000 most variable CpG sites (Additional file 3), PC and AN samples clustered separately, demonstrating that DNA methylation levels are substantially altered in PC tissue samples, consistent with previous reports [3, 5]. Moreover, by RNAseq, we found that transcriptional expression of 2314 genes was significantly deregulated in PC vs. AN samples after adjusting for multiple testing (BH adj. p < 0.05, exactTest in EdgeR package (negative binomial test)). The corresponding MDS plot of the 150 most variable transcripts clearly separated PC and AN samples (Additional file 3). As a technical validation of RNAseq results, RT-qPCR of MEIS2 on 7 prostate cell lines (2 benign and 5 malignant) and primary prostate epithelial cells was performed and showed that RT-qPCR results were highly comparable to RNAseq results (data not shown; Spearman’s rho 1.00, p < 0.001).
Next, to identify epigenetically deregulated gene expression, 450K and RNAseq datasets were merged by gene name and combined p values (Fisher’s p values) were calculated from the p values obtained by comparing PC and AN samples for each data type. A total of 1125 genes had a combined BH-adjusted p value < 0.05, indicating significant deregulation in DNA methylation and/or RNA expression levels in PC tissue samples. Functional annotation analysis using DAVID [24, 25] revealed significant enrichment for homeobox and homeodomain-containing genes (enrichment score of homeobox cluster, 2.79; the 4th most enriched cluster; Fig. 2) in our list of significantly epigenetically deregulated genes in PC. The homeobox cluster included genes such as HOXC6 and DLX1 which are included in the SelectMDx® urine-based test for detection of high-risk PC , as well as EN2, GLI3, and MEIS2 previously investigated for their biomarker potential and/or function in PC [18, 33, 34]. Functional annotation analyses on the top 3000 differentially methylated and 2314 differentially expressed genes, respectively, overall yielded similar results compared to the 1125 epigenetically deregulated genes (Additional file 4). Moreover, the homeobox clusters were the 6th (methylation) and 51st (expression) most enriched clusters, respectively, suggesting that transcriptional expression is not altered in all differentially methylated homeobox genes.
DNA methylation and RNA expression of MEIS2
Next, to assess correlations between DNA methylation and transcriptional expression, Spearman’s rho coefficients were calculated for each individual CpG site and its corresponding gene (using the 34 PC and AN samples analyzed by both methods). This analysis revealed significant differential methylation as well as downregulation of MEIS2, a homeobox gene previously investigated for its tumor suppressive role in PC initiation and progression [16,17,18]. More specifically, in the 450K data, MEIS2 was significantly hypermethylated in PC samples both in the promoter region and from exon 2–5 (cg06933370, BH adj. p < 0.0001; cg25181383, BH adj. p < 0.0001; Mann-Whitney test; Fig. 3a, b) and the transcriptional expression of MEIS2 was significantly downregulated in PC compared to AN samples in the RNAseq dataset (BH adj. p < 0.0001, Mann-Whitney test; Fig. 3c), consistent with epigenetic silencing through promoter hypermethylation. Moreover, we found a significant inverse correlation between MEIS2 transcriptional expression and DNA methylation both in the promoter and intronic region (cg25381383 (promoter region)/cg06933370 (intronic region), rho = − 0.5233/− 0.8197, BH adj. p = 0.0036/p < 0.0001; Spearman correlation; Fig. 3d).
This finding was subsequently confirmed in a large independent RP cohort from TCGA including 497 PC and 52 AN samples with both 450K and RNAseq data available . MEIS2 was significantly hypermethylated (cg06933370 and cg25381383, BH adj. p < 0.0001, Mann-Whitney test; Fig. 3e) and downregulated (p < 0.0001, Mann-Whitney test; Fig. 3f) in PC samples compared to AN samples in the TCGA cohorts. Likewise, DNA methylation of cg06933370 and cg25381383 was significantly inversely correlated with MEIS2 transcriptional expression also in the large TCGA cohort (rho = − 0.6739 and − 0.4864, BH adj. p < 0.0001; Spearman correlation; Fig. 3g).
Together, these results strongly indicate that aberrant DNA hypermethylation is associated with epigenetic silencing of MEIS2 transcriptional expression in PC.
Prognostic potential of MEIS2 RNA expression
Next, to test the possible prognostic potential of MEIS2 at the transcriptional level, we used data from three publicly available PC patient cohorts: RNAseq data from Long et al. (n = 106 RP patients) , microarray expression data from Taylor et al. (n = 126 RP patients) , and RNAseq data from TCGA (n = 495 RP patients) . Low MEIS2 expression was generally associated with unfavorable clinicopathological parameters (high pathological Gleason score, advanced pathological T stage, and/or positive surgical margins), although this was only statistically significant for Gleason score in the Taylor and TCGA cohort, and for pathological T-stage in the TCGA cohort (Additional file 5). Furthermore, in all three RP cohorts, low MEIS2 RNA expression was significantly associated with poor BCR-free survival in Kaplan-Meier (log-rank test, p = 0.0084, 0.0001, and 0.0191, respectively; Fig. 4a–c) and univariate cox regression analyses (p = 0.010, 0.000, and 0.022), respectively (Additional files 6, 7 and 8). After adjustment for routine clinicopathological parameters, MEIS2 expression remains a significant predictor of BCR only in the Long cohort (p = 0.005, HR = 0.39 (0.20–0.76); Additional file 6). Similar results were obtained when analyzing MEIS2 transcriptional expression as a continuous variable (Additional files 6, 7 and 8).
Together, these results indicate that low MEIS2 transcriptional expression is associated with more aggressive PC, thereby expanding on a previous report that linked low MEIS2 protein levels with poor overall survival in PC .
Prognostic potential of MEIS2 DNA methylation
Next, to examine the prognostic potential of MEIS2 DNA methylation for prediction of post-operative BCR, we designed two qMSP assays each covering a significantly hypermethylated region of the MEIS2 gene (intragenic and promoter, respectively) identified above by 450K analysis (Fig. 3a, b). Thus, qMSP assay 1 covered three probes from the 450K array within the MEIS2 gene (cg06933370, cg23677243, and cg26708220) and assay 2 covered one probe from the 450K array in the promoter region of MEIS2 (cg25381383). qMSP analyses were performed on an independent set of 195 PC, 17 AN, and 6 BPH samples (Table 1). We combined the AN and BPH samples into one non-malignant sample group, as there were no significant difference in MEIS2 methylation between these sample types. For both assays, MEIS2 was significantly hypermethylated in PC tissue samples (p < 0.0001, Mann-Whitney test; Fig. 5a, b left) and showed promising diagnostic potential with AUCs of 0.841 and 0.917, respectively (Fig. 5a, b right). For both assays, high MEIS2 DNA methylation was generally associated with adverse clinicopathological factors (high pathological Gleason score, advanced pathological T stage, positive surgical margins, and/or high CAPRA-S score (score for prediction of post-operative BCR)), although this was only statistically significant for pathological T-stage and CAPRA-S score (Additional file 9). Consistent with this, patients with low MEIS2 methylation had significantly lower CAPRA-S score (p = 0.0066 and 0.0419, respectively; Mann-Whitney test; Fig. 6). Combining the two qMSP assays into one model (low methylation in both assays vs. high methylation in at least one assay), low methylation was still significantly associated with low CAPRA-S score (p = 0.0004; Mann-Whitney test. Fig. 6).
Furthermore, for each assay, high MEIS2 methylation was significantly associated with short BCR-free survival in Kaplan-Meier analysis (qMSP assay 1/2, p = 0.0248/p = 0.0497; log-rank test; Fig. 7a, b). Similar results were obtained for the individual qMSP assays by univariate cox regression analysis but was only borderline significant for assay 2 (assay 1, p = 0.026, HR = 1.57 (1.06–2.34), C-index = 0.564; assay 2, p = 0.051, HR = 1.49 (1.00-2.22), C-index = 0.546; Additional files 10 and 11). However, the prognostic power was improved by combining the two assays, as low methylation for both assays (compared to high methylation for at least one of the assays) was associated with significantly better post-operative BCR-free survival in both Kaplan-Meier (p = 0.0068, log-rank test; Fig. 7c) and univariate cox regression analysis (p = 0.008, HR = 1.79 (1.17–2.76), C-index = 0.573; Table 2). However, the combined MEIS2 methylation model did not remain significant after adjustment for routine clinicopathological variables (p = 0.681, HR = 1.10 (0.69–1.75), c-index = 0.730, multivariate cox regression; Table 2).
Next, for external validation, we used the PC patient cohort from TCGA (450K data; assay 1, average methylation of cg06933370, cg23677243, and cg26708220; assay 2, methylation of cg25381383). In this cohort, high MEIS2 DNA methylation was generally associated with adverse clinicopathological parameters although only significant for probes corresponding to assay 1 regarding pathological T-stage and for probes corresponding to assay 2 regarding both Gleason score and pathological T-stage (Additional file 5). Moreover, low MEIS2 DNA methylation was associated with significantly better post-operative BCR-free survival in uni- and multivariate analysis ( “Assay 1”/”Assay 2”, p = 0.024/0.017, HR = 1.98 (1.09–3.59)/2.09 (1.14–3.82), C-index = 0.712/0.717, multivariate cox regression. Additional files 12 and 13).
Likewise, low methylation of the combined model was associated with significantly better BCR-free survival also in the TCGA cohort by both Kaplan-Meier (p = 0.0067, Log-rank test; Fig. 7d) and univariate cox regression analysis (p = 0.009, HR = 2.48 (1.26–4.88); Table 3). Furthermore, after adjusting for pathological Gleason score and T-stage, the combined MEIS2 methylation model remained a significant independent predictor of BCR (p = 0.042, HR = 2.02 (1.02–3.99); Table 3). In summary, this study is the first to demonstrate a significant association between aberrant MEIS2 hypermethylation and adverse clinical outcome of PC.
In this study, we performed an integrative analysis of DNA methylation and RNA expression data from PC, as this approach may identify novel candidate driver genes involved in PC development or progression. Accordingly, we conducted DNA methylation profiling (450K array) and RNAseq on PC and AN samples from 29 radical prostatectomy patients and integrated the data for biomarker discovery. Among the top candidate genes with significantly altered DNA methylation and/or RNA expression levels in PC, we found significant enrichment for homeobox genes, incl. MEIS2. RNA expression and DNA methylation of MEIS2 were inversely correlated in our discovery cohort, which was confirmed in a large independent RP patient cohort from TCGA (495 PC, 36 AN), suggesting epigenetic silencing. Furthermore, low transcriptional expression and DNA hypermethylation of MEIS2 was associated with post-operative BCR in multiple independent RP patient cohorts, including more than 700 PC patients in total. To the best of our knowledge, this is the first study to demonstrate significant prognostic value of MEIS2 epigenetic silencing in PC.
The present study is the first to examine the prognostic potential of MEIS2 transcriptional expression in PC. We used three public RNA expression dataset from RP cohorts (> 700 patients with clinical follow-up) and found that low MEIS2 expression was associated with post-operative BCR and generally correlated with adverse clinicopathological parameters. After adjustment for routine clinicopathological variables, MEIS2 expression remained a significant predictor of BCR in multivariate analysis in the Long cohort, but not in the TCGA or Taylor cohort, possible due to differences in the exact composition of the cohorts and/or in the methodologies used for expression profiling. Specifically, the Long RNAseq data was derived from FFPE tissue samples, whereas Taylor et al. used microarrays and TCGA used fresh frozen tissue samples for RNAseq analyses. Furthermore, the fraction of patients with BCR varied from 12% in TCGA to 25% in the Taylor cohort and 52% in the Long cohort, possibly affecting the statistical power. Nevertheless, our results showed that low MEIS2 transcriptional expression was associated with significantly shorter BCR-free survival in all three cohorts, which corroborates and expands on previous findings of a significant association between MEIS2 protein expression and short overall survival in PC .
We found a significant inverse correlation between MEIS2 DNA methylation and RNA expression, suggesting that MEIS2 is epigenetically silenced in PC. The observed silencing of MEIS2 could also be affected by other factors, e.g., transcriptional regulators expressed during PC development/progression, but further studies are needed to evaluate this. While our finding of frequent aberrant hypermethylation of MEIS2 in PC is consistent with results from two previous studies exploring the methylome in PC [36, 37], our demonstration of diagnostic and prognostic potential of MEIS2 DNA methylation has not been described before. Here, we found that MEIS2 hypermethylation had diagnostic potential for PC with a high AUC (AUC assay 1/2, 0.814/0.917), comparable to AUCs reported for previously published candidate methylation markers for PC (AUCs 0.794–0.980) [5, 6, 8, 9, 11, 38]. Some of these previously published DNA methylation markers could predict BCR after RP independent of routine clinicopathological parameters [5, 6]; however, this was only the case for MEIS2 in the 450K data from TCGA and not in the qMSP cohort. In addition, we found that high MEIS2 methylation levels were significantly associated with short BCR-free survival and high CAPRA-S score. MEIS2 methylation also predicted time to BCR independently of clinicopathological parameters in our external validation cohort from TCGA. However, the prognostic potential of MEIS2 DNA methylation warrants further validation in a large cohort with long clinical follow-up as the TCGA cohort used for independent validation had short clinical follow-up and few events (BCR in 12% of patients). In this study, we made a prognostic model based on MEIS2 methylation. Conceivably, adding RNA expression data to this model may potentially improve predictive accuracy, but we found no significant evidence for this in the TCGA cohort (data not shown). However, further studies on large cohorts with overlapping DNA methylation and RNA expression data and long follow-up are needed to investigate this.
MEIS2 RNA and protein expression has previously been reported as downregulated in primary PC and further downregulated in metastatic PC tissue, as compared to non-malignant prostate tissue samples [16, 17]. Furthermore, in a small cohort of 83 Gleason 6 patients, a low protein level of MEIS2 was significantly associated with short overall survival . Also, downregulation of MEIS2 has been reported to play an important functional role in progression to castration resistant PC . This corresponds with our findings of low MEIS2 transcriptional expression being associated with more aggressive PC, defined by short BCR-free survival. Although epigenetic silencing of MEIS2 has also been described in lung and hepatocellular carcinoma cell lines as well as colorectal cancer [39,40,41], the function of MEIS2 in other cancers remains unclear and may be highly disease-specific as MEIS2 is downregulated in some cancers and upregulated in others . In neuroblastoma, leukemia, and multiple myeloma, MEIS2 serves as an oncogene [42,43,44,45] and, similarly seems to play an important role in tumor cell migration and invasion in bladder and colorectal cancer [46, 47]. MEIS2 expression may also be involved in chemotherapy sensitivity, although current results are conflicting. Thus, MEIS2 knockdown increases responsiveness to chemotherapy in multiple myeloma, whereas MEIS2 is downregulated in colorectal cancer patients resistant towards oxaliplatin-based chemotherapy [41, 45]. Contradictory results also exist regarding prognosis. In ovarian cancer, high MEIS2 protein expression has been associated with improved prognosis whereas high RNA expression of MEIS2 has been associated with worse overall survival in colorectal cancer [47, 48]. Thus, the functional role of MEIS2 seems to be highly disease-specific and further studies are needed to clarify this. In the present study, we examined the prognostic potential of MEIS2 DNA methylation and transcriptional expression in five RP cohorts. Risk stratification of PC patients after RP is important to identify patients who would benefit from adjuvant therapy and to avoid treating patients with low risk of BCR. In this setting, a prognostic molecular marker as MEIS2 could improve the current management of RP patients. Additionally, it is also of clinical relevance to test the prognostic value of MEIS2 in diagnostic biopsies, as this could give important clues to whether MEIS2 DNA methylation and/or RNA expression can improve the accuracy of PC prognosis at the time of diagnosis. Currently, only pre-operative clinicopathological parameters are available at diagnosis which are upgraded and/or upstaged after RP in more than 50% of PC tumors . Thus, improved prediction of aggressiveness at diagnosis by addition of a molecular marker could improve treatment decisions and reduce overtreatment of indolent PCs.
A potential limitation to our study is that we could not discriminate between Gleason scores 3+4 and 4+3, as this information was not available for the public PC patient sets. Moreover, we have used BCR as endpoint in our survival analyses. BCR is only a surrogate marker for aggressiveness and more clinically relevant endpoints as metastatic disease or PC-specific mortality should be used in future studies instead. However, this would require > 15 years of clinical follow-up due to the generally slow natural history of PC progression .
In conclusion, we here show that the homeobox gene MEIS2 is epigenetically silenced in PC. To the best of our knowledge, this is the first study to investigate, demonstrate, and independently validate a prognostic biomarker potential for MEIS2 at the transcriptional expression level and at the DNA methylation level in PC.
Availability of data and materials
The data generated in this study is available from the corresponding author upon reasonable request.
The RNAseq, 450K, and clinical data from TCGA are available at the Cancer Genome Atlas (http://cancergenome.nih.gov/) . From Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21034), the RNAseq and clinical data from Long et al. (GSE54460)  as well as the microarray and clinical data from Taylor et al. (GSE21034)  are available.
Illumina 450K DNA methylation array
Area under the curve
Benign prostatic hyperplasia
Counts per million
Castration resistant prostate cancer
Database for Annotation, Visualization and Integrated Discovery
Hematoxylin and eosin
Quantitative methylation specific PCR
Receiver operating characteristics
Three amino acid loop extensions
The Cancer Genome Atlas
Transurethral resection of the prostate
Ferlay J, Steliarova-Foucher E, Lortet-Tieulent J, Rosso S, Coebergh JW, Comber H, et al. Cancer incidence and mortality patterns in Europe: estimates for 40 countries in 2012. Eur J Cancer. 2013;49(6):1374–403.
Loeb S, Bjurlin MA, Nicholson J, Tammela TL, Penson DF, Carter HB, et al. Overdiagnosis and overtreatment of prostate cancer. Eur Urol. 2014;65(6):1046–55.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.
Sproul D, Kitchen RR, Nestor CE, Dixon JM, Sims AH, Harrison DJ, et al. Tissue of origin determines cancer-associated CpG island promoter hypermethylation patterns. Genome Biol. 2012;13(10):R84.
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(26):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(8):2169–81.
Strand SH, Orntoft TF, Sorensen KD. Prognostic DNA methylation markers for prostate cancer. Int J Mol Sci. 2014;15(9):16544–76.
Norgaard M, Haldrup C, Storebjerg TM, Vestergaard EM, Wild PJ, Hoyer S, et al. Comprehensive evaluation of TFF3 promoter hypomethylation and molecular biomarker potential for prostate cancer diagnosis and prognosis. Int J Mol Sci. 2017;18(9).
Haldrup C, Pedersen AL, Ogaard N, Strand SH, Hoyer S, Borre M, et al. Biomarker potential of ST6GALNAC3 and ZNF660 promoter hypermethylation in prostate cancer tissue and liquid biopsies. Mol Oncol. 2018;12(4):545–60.
Banez LL, Sun L, van Leenders GJ, Wheeler TM, Bangma CH, Freedland SJ, et al. Multicenter clinical validation of PITX2 methylation as a prostate specific antigen recurrence predictor in patients with post-radical prostatectomy prostate cancer. J Urol. 2010;184(1):149–56.
Haldrup C, Lynnerup AS, Storebjerg TM, Vang S, Wild P, Visakorpi T, et al. Large-scale evaluation of SLC18A2 in prostate cancer reveals diagnostic and prognostic biomarker potential at three molecular levels. Mol Oncol. 2016;10(6):825–37.
Strand SH, Switnicki M, Moller M, Haldrup C, Storebjerg TM, Hedegaard J, et al. RHCG and TCAF1 promoter hypermethylation predicts biochemical recurrence in prostate cancer patients treated by radical prostatectomy. Oncotarget. 2017;8(4):5774–88.
Brechka H, Bhanvadia RR, VanOpstall C, Vander Griend DJ. HOXB13 mutations and binding partners in prostate development and cancer: function, clinical significance, and future directions. Genes Dis. 2017;4(2):75–87.
Xu J, Lange EM, Lu L, Zheng SL, Wang Z, Thibodeau SN, et al. HOXB13 is a susceptibility gene for prostate cancer: results from the International Consortium for Prostate Cancer Genetics (ICPCG). Hum Genet. 2013;132(1):5–14.
Storebjerg TM, Hoyer S, Kirkegaard P, Bro F, Orntoft TF, Borre M, et al. Prevalence of the HOXB13 G84E mutation in Danish men undergoing radical prostatectomy and its correlations with prostate cancer risk and aggressiveness. BJU Int. 2016;118(4):646–53.
Chen JL, Li J, Kiriluk KJ, Rosen AM, Paner GP, Antic T, et al. Deregulation of a Hox protein regulatory network spanning prostate cancer initiation and progression. Clin Cancer Res. 2012;18(16):4291–302.
Bhanvadia RR, VanOpstall C, Brechka H, Barashi NS, Gillard M, McAuley EM, et al. MEIS1 and MEIS2 expression and prostate cancer progression: a role for HOXB13 binding partners in metastatic disease. Clin Cancer Res. 2018;24(15):3668–80.
Jeong JH, Park SJ, Dickinson SI, Luo JL. A constitutive intrinsic inflammatory signaling circuit composed of miR-196b, Meis2, PPP3CC, and p65 drives prostate cancer castration resistance. Mol Cell. 2017;65(1):154–67.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics (Oxford, England). 2015;31(2):166–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England). 2010;26(1):139–40.
Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics (Oxford, England). 2014;30(3):428–30.
Fisher RA. Statistical Methods For Research Workers. Edinburgh: Oliver and Boyd; 1925.
Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Huang d W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.
Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics (Oxford, England). 2007;23(10):1289–91.
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(12):3228–37.
TCGA. The molecular taxonomy of primary prostate cancer. Cell. 2015;163(4):1011–25.
Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. Evaluation of the Infinium Methylation 450K technology. Epigenomics. 2011;3(6):771–84.
Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18(1):11–22.
Van Neste L, Hendriks RJ, Dijkstra S, Trooskens G, Cornel EB, Jannink SA, et al. Detection of high-grade prostate cancer using a urinary molecular biomarker-based risk score. Eur Urol. 2016;70(5):740–8.
Gonnissen A, Isebaert S, Perneel C, McKee CM, Verrill C, Bryant RJ, et al. Tissue microarray analysis indicates hedgehog signaling as a potential prognostic factor in intermediate-risk prostate cancer. BMC Cancer. 2017;17(1):634.
McGrath SE, Michael A, Morgan R, Pandha H. EN2 in Prostate Cancer. Adv Clin Chem. 2015;71:47–76.
Taylor BS. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18(1 SRC - GoogleScholar):11–22.
Long MD, Smiraglia DJ, Campbell MJ. The genomic impact of DNA CpG methylation on gene expression; relationships in prostate cancer. Biomolecules. 2017;7(1).
Kim JH, Dhanasekaran SM, Prensner JR, Cao X, Robinson D, Kalyana-Sundaram S, et al. Deep sequencing reveals distinct patterns of DNA methylation in prostate cancer. Genome Res. 2011;21(7):1028–41.
Bjerre MT, Strand SH, Norgaard M, Kristensen H, Rasmussen AK, Mortensen MM, et al. Aberrant DOCK2, GRASP, HIF3A and PKFP hypermethylation has potential as a prognostic biomarker for prostate cancer. Int J Mol Sci. 2019;20(5).
Kiehl S, Zimmermann T, Savai R, Pullamsetti SS, Seeger W, Bartkuhn M, et al. Epigenetic silencing of downstream genes mediated by tandem orientation in lung cancer. Sci Rep. 2017;7(1):3896.
Sun N, Zhang J, Zhang C, Shi Y, Zhao B, Jiao A, et al. Using Illumina Infinium HumanMethylation 450K BeadChip to explore genomewide DNA methylation profiles in a human hepatocellular carcinoma cell line. Mol Med Rep. 2018;18(5):4446–56.
Wang X, Ghareeb WM, Zhang Y, Yu Q, Lu X, Huang Y, et al. Hypermethylated and downregulated MEIS2 are involved in stemness properties and oxaliplatin-based chemotherapy resistance of colorectal cancer. J Cell Physiol. 2019;234(10):18180–91.
Gross A, Schulz C, Kolb J, Koster J, Wehner S, Czaplinski S, et al. Tumorigenic and antiproliferative properties of the TALE-transcription factors MEIS2D and MEIS2A in neuroblastoma. Cancer Res. 2018;78(8):1935–47.
Zha Y, Xia Y, Ding J, Choi JH, Yang L, Dong Z, et al. MEIS2 is essential for neuroblastoma cell survival and proliferation by transcriptional control of M-phase progression. Cell Death Dis. 2014;5:e1417.
Lai CK, Norddahl GL, Maetzig T, Rosten P, Lohr T, Sanchez Milde L, et al. Meis2 as a critical player in MN1-induced leukemia. Blood Cancer J. 2017;7(9):e613.
Abruzzese MP, Bilotta MT, Fionda C, Zingoni A, Soriani A, Petrucci MT, et al. The homeobox transcription factor MEIS2 is a regulator of cancer cell survival and IMiDs activity in multiple myeloma: modulation by Bromodomain and extra-terminal (BET) protein inhibitors. Cell Death Dis. 2019;10(4):324.
Xie R, Chen X, Chen Z, Huang M, Dong W, Gu P, et al. Polypyrimidine tract binding protein 1 promotes lymphatic metastasis and proliferation of bladder cancer via alternative splicing of MEIS2 and PKM. Cancer Lett. 2019;449:31–44.
Wan Z, Chai R, Yuan H, Chen B, Dong Q, Zheng B, et al. MEIS2 promotes cell migration and invasion in colorectal cancer. Oncol Rep. 2019;42(1):213–23.
Crijns AP, de Graeff P, Geerts D, Ten Hoor KA, Hollema H, van der Sluis T, et al. MEIS and PBX homeobox proteins in ovarian cancer. Eur J Cancer. 2007;43(17):2495–505.
Koksal IT, Ozcan F, Kadioglu TC, Esen T, Kilicaslan I, Tunc M. Discrepancy between Gleason scores of biopsy and radical prostatectomy specimens. Eur Urol. 2000;37(6):670–4.
Albertsen PC, Hanley JA, Fine J. 20-year outcomes following conservative management of clinically localized prostate cancer. JAMA. 2005;293(17):2095–101.
The authors thank Maria Engtoft Skjøtt and Birgitte Trolle for excellent technical assistance throughout the project and Jacob Christian Fredsøe for downloading and managing the TCGA data. Additionally, the Danish Cancer Biobank is acknowledged for biological material.
This work was supported by grants from The Danish Cancer Society, The Velux Foundation, Innovation Fund Denmark, Fabrikant Vilhelm Petersen og Hustrus Legat, and The Central Denmark Region Health Research Fund.
Ethics approval and consent to participate
The study was approved by the Central Denmark Region Committees on Health Research Ethics and by the Danish Data Protection Agency. Written informed consent was obtained from all patients.
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.
Additional file 4: Figure S2. Functional clustering of differentially methylated and expressed genes (discovery cohort) using DAVID. For the 450K data, the 3000 genes with the lowest BH-adjusted p-value (PC vs. AN) were used for the analysis. Top 10 enriched clusters are shown. For the RNAseq data, the 2314 genes with a significant BH-adjusted p-value (PC vs. AN) were used as input. Top 10 enriched clusters as well as cluster number 51 (homeobox-cluster) are shown. Barplots show enrichment scores of the clusters. A general term describing the genes/categories within each cluster is given.
Additional file 5: Figure S3. Association between RNA expression and DNA methylation of MEIS2 in public cohorts. Clinicopathological variables examined: pathological Gleason score, pathological T-stage, and surgical margin status. Public cohorts examined: Long et al. (RNAseq), Taylor (microarray), TCGA (RNAseq and 450K). P-values were calculated using Mann-Whitney tests. Path., pathological. Surg., surgical. P, p-value.
Additional file 9: Figure S4. Association between DNA methylation of MEIS2 assay 1 and 2 and clinicopathological variables in the qMSP cohort. Clinicopathological variables examined: pathological Gleason score, pathological T-stage, surgical margin status, and CAPRA-S score (low: 0-2, intermediate: 3-5, high: ≥6). P-values were calculated using Mann-Whitney tests. Path., pathological. Surg., surgical. P, p-value.
Additional file 12: Table S8. Uni- and multivariate cox regression of the average β-value of cg06933370, cg23677243, and cg26708220 (to mimic qMSP assay 1) in the TCGA 450K cohort (n=389 patients). BCR was used as end-point. Meth: Methylation. Path.: Pathologic. HR: Hazard ration. CI: Confidence interval.
About this article
Cite this article
Nørgaard, M., Haldrup, C., Bjerre, M.T. et al. Epigenetic silencing of MEIS2 in prostate cancer recurrence. Clin Epigenet 11, 147 (2019) doi:10.1186/s13148-019-0742-x
- Prostate cancer
- DNA methylation
- RNA expression
- Epigenetic silencing