Longitudinal analysis of individual cfDNA methylome patterns in metastatic prostate cancer
Clinical Epigenetics volume 13, Article number: 168 (2021)
Disease progression and therapeutic resistance are hallmarks of advanced stage prostate cancer (PCa), which remains a major cause of cancer-related mortality around the world. Longitudinal studies, coupled with the use of liquid biopsies, offer a potentially new and minimally invasive platform to study the dynamics of tumour progression. Our aim was to investigate the dynamics of personal DNA methylomic profiles of metastatic PCa (mPCa) patients, during disease progression and therapy administration.
Forty-eight plasma samples from 9 mPCa patients were collected, longitudinally, over 13–21 months. After circulating cell-free DNA (cfDNA) isolation, DNA methylation was profiled using the Infinium MethylationEPIC BeadChip. The top 5% most variable probes across time, within each individual, were utilised to study dynamic methylation patterns during disease progression and therapeutic response. Statistical testing was carried out to identify differentially methylated genes (DMGs) in cfDNA, which were subsequently validated in two independent mPCa (cfDNA and FFPE tissue) cohorts. Individual cfDNA global methylation patterns were temporally stable throughout the disease course. However, a proportion of CpG sites presented a dynamic temporal pattern that was consistent with clinical events, including different therapies, and were prominently associated with genes linked to immune response pathways. Additionally, study of the tumour fraction of cfDNA identified > 2000 DMGs with dynamic methylation patterns.
Longitudinal assessment of cfDNA methylation in mPCa patients unveiled dynamic patterns associated with disease progression and therapy administration, thus highlighting the potential of using liquid biopsies to study PCa evolution at a methylomic level.
Despite its 5-year survival rate of nearly 100% for localized disease, prostate cancer (PCa) is still considered a major cause of morbidity and mortality, accounting for approximately 7.6% of all cancer-related deaths worldwide . Numerous studies have highlighted the importance of epigenetic aberrations, particularly DNA methylation, in PCa [2,3,4]. However, a knowledge-gap surrounds the dynamics of epigenetic alterations during disease progression and treatment resistance. Such investigations are typically hampered by the validity of performing molecular characterisation on archival tumour tissue acquired at diagnosis, often several years before disease progression. Moreover, acquiring serial biopsies or sampling multiple lesions to study tumour dynamics remains impractical and unethical.
Longitudinal studies have long been used to decipher the temporal association of disease traits in an effort to elucidate specific biological pathways and better understand variability observed among individuals. However, there is a dearth of longitudinal epigenetic studies in cancer, as most studies have focused on associations between epigenetic modifications and aging [5, 6]. Liquid biopsies provide an opportunity to overcome challenges associated with tissue biopsies and incorporate longitudinal assessment into the management of cancer patients. Circulating tumour DNA (ctDNA) can be noninvasively evaluated in blood plasma and several studies have shown that PCa-specific methylation alterations can be studied using this approach [7,8,9], with recent evidence showing that ctDNA profiles are consistent with metastatic tissue . Additionally, since ctDNA is shed from different focal lesions, it overcomes spatial and clonal heterogeneity , thus making liquid biopsies an ideal surrogate to study the dynamics of tumour progression and real-time disease monitoring.
The aim of this study was to investigate the dynamics of individual methylomic profiles of metastatic PCa (mPCa) patients. Through longitudinal analysis of circulating cell-free DNA (cfDNA), we demonstrate the feasibility of this approach to study tumour dynamics and monitor disease-associated characteristics.
Study cohort characteristics
Longitudinal assessment of methylation dynamics in 9 mPCa patients (comprising 48 cfDNA samples) was performed over at least 4 time points for an average of 19.11 months (Table 1). The shortest follow-up time (13 months) was attributed to death from disease. During follow-up, subjects received up to four different therapies: 88.9% had luteinizing hormone-releasing hormone (LHRH) agonists/antagonists, 100% had chemotherapy (taxanes), 88.9% had androgen receptor (AR) inhibitors and 33.3% had radioisotopes.
To assess the prostate content in the cfDNA samples, we utilised the human cell-type DNA methylation atlas . Sixteen of 48 samples (33.3%), across 6 patients, had detectable prostate-derived DNA; 10 of whom had prostate content (PC) ≥ 10% (Fig. 1A). Samples with the largest prostate contribution came from two subjects, one (S019) had the highest tumour burden (as measured by prostate-specific antigen, PSA), and the second (S008) died from disease during follow-up. For the remaining samples, sources of cfDNA were largely attributed to cells of hematopoietic origin, as previously observed .
Global cfDNA methylation profiles are temporally stable
In order to understand the dynamics of cfDNA methylation during PCa progression and sequential therapies, we analysed the global methylation patterns of all samples. Principal component analysis was used to understand methylation differences between individuals, with this analysis showing little variance between samples and/or individuals (Fig. 1B). Interestingly, the biggest difference stemmed from prostate cellular contribution, as samples from S019 clearly separated from all other samples. We next investigated the degree to which global methylation changed over time within individuals. Evaluation of temporal divergence showed that most probes remained relatively stable (standard deviation—SD < 0.1), indicating a low intra-subject variation in cfDNA methylation (Fig. 1C). Again, S019 was distinct, displaying the greatest degree of temporal variation.
The dynamic patterns of the 5′-regulatory probes within the top 5% most variable probes (MVPs) (Additional file 1: Fig. S1) were analysed for congruence over time, and assessed in conjunction with the disease course for each individual (Fig. 2A and Additional file 1: Figs. S2–S9). Probes with similar patterns were aggregated into clusters. Individuals displayed as few as two probe clusters (S019) and as many as seven (S007) (Fig. 2B and Additional file 1: Figs. S2–S9). Greater numbers of probe clusters suggest that the most variably methylated genes in these patients behave differently throughout the disease course and therapeutic response. Indeed, some methylated gene clusters correlated with clinical events. For example, in subjects with high PC samples (S008 and S019), methylation dynamics mimicked PSA patterns, with methylation shifts coinciding with PSA changes observed in these patients (Fig. 2 and Additional file 1: Fig. S8). Findings also indicate that therapy administration may affect cfDNA methylation dynamics, with most subjects demonstrating a notable shift in methylation patterns following taxane administration (Fig. 2C and Additional file 1: Figs. S2–S9). This is evident, for example, by the methylation difference observed in S019 between administration of cabazitaxel at baseline (BL) and follow-up 1 (FU1) (Fig. 2C). AR inhibition may also confer temporal epigenomic shifts, as shown by methylation changes in 4/6 probe clusters in S007, after commencing enzalutamide at month 11 (Additional file 1: Fig. S7).
Overall, although the majority of the methylome remained stable over time, a proportion of 5′ regulatory CpGs had a dynamic longitudinal pattern, consistent with clinical events.
Dynamic methylation patterns during therapy administration are associated with an immune response element
To further explore the relationship between cfDNA methylation changes and therapy administration, we examined each subject’s MVPs and measured methylation changes (absolute β ≥ 0.1 difference) between time-points after therapy began. For most subjects, over 50% of the MVPs displayed a methylation change after taxane (Fig. 3A) or AR inhibitor (Fig. 3B). Notably, both drugs were administered simultaneously for S017 and S029, hindering any conclusions about methylation changes that might be associated with a specific therapy.
A detailed look at the cluster analysis also uncovered a methylation shift after taxane termination, such as the ones observed after FU1 in S019, bringing methylation patterns to a pre-therapy level (Fig. 2C; other examples are shown in Additional file 1: Figs. S2–S9). Indeed, reviewing all hyper- and hypo-MVPs showed that most were only transiently affected by taxanes, with their methylation shift largely dissipating upon therapy cessation (Additional file 1: Fig. S10).
Additionally, we investigated if the probes that changed after therapy administration were associated with any biological processes and pathways. Studying the most significantly represented processes associated with both taxane and AR therapy administration (FDR p < 0.05), we observed that a high proportion of gene ontology (GO) terms were associated with immune response (Fig. 3C, D). Other processes were also enriched, such as signaling, cellular components and neurological system. This might be expected, since most cfDNA samples had substantial contribution from cells of hematopoietic origin. Notably, S008 and S019 were distinct; their cfDNA was enriched for transcriptional regulation and multicellular processes, perhaps reflecting their higher PC and thus genes involved in PCa (Fig. 3C, D).
ctDNA temporal dynamics in mPCa patients
Irrespective of therapy administration, we also examined the tumour fraction of cfDNA methylation patterns. To do this, we used the top 5% promoter-associated MVPs to identify genes with dynamic methylation patterns among subjects with PC samples (S008, S009, S019, S029). After assessment of their methylation values across samples with no prostate content (NPC) from the whole cohort (to verify their ctDNA-specificity), a list of 3435 differentially methylated genes (DMGs) was compiled, 1991 hyper- and 1444 hypomethylated (Fig. 4A). Notably, several DMGs are already known to be involved in PCa carcinogenesis. For example, in our cohort, both APC and RASSF1 demonstrated higher methylation, and CPEB4 and EPN1 showed lower methylation, when comparing PC versus NPC samples (Fig. 4B). The majority (72%) of DMGs (1574 hyper- and 900 hypomethylated) validated in an independent mPCa cfDNA cohort, where DMGs were also compared between PC and NPC samples  (Fig. 4C). Additionally, overlap with the publicly-available COSMIC methylation dataset, showed that methylation changes in 53% and 36% of these hyper and hypomethylated genes, respectively, has been previously implicated in PCa (Fig. 4C).
Further validation was carried out in a formalin-fixed, paraffin-embedded (FFPE) tissue cohort of six distinct stages of PCa, to account for potential confounding cfDNA methylation originating from non-prostate sources. Certain hyper- and hypomethylated DMGs show similar methylation patterns across PCa stages, which might suggest that these DMGs are prostate-specific and differences can only be seen when comparing prostate with other tissues/cfDNA sources. However, the methylation values for most mPCa DMGs were similar to observations in the metastatic stage of the FFPE cohort, validating our results and suggesting an association of these epigenetic alterations and the more advanced stages of PCa (Fig. 4D). Additionally, methylation differences were observed across the different stages of PCa, emphasizing the potential role of DNA methylation in disease progression (Fig. 4D and Additional file 1: Fig. S11).
The majority of both hyper and hypomethylated DMGs (60.37% and 93.98%, respectively) were notably unique to one subject (Fig. 4A). This again demonstrates the renowned molecular heterogeneity of PCa, and the need for personalised medicine approaches such as this.
In this study, we describe the first personal methylomic longitudinal analysis of PCa using cfDNA from 9 mPCa patients. Our results reveal that distinct methylation dynamics, observed in a proportion of assessed CpG sites, were consistent with different clinical events, such as PSA progression and therapy administration. Additionally, specifically studying ctDNA methylation dynamics allowed us to identify genes previously observed as epigenetically modified in mPCa [2, 7, 14].
Longitudinal studies are defined by following individuals over time to identify disease-specific temporal patterns that relate to stage or burden of disease, or response to a given therapy. Understanding these specific patterns may result in novel diagnostic, prognostic and monitoring tools that could improve clinical management, and also help reduce the burden on the healthcare system. This is particularly relevant for advanced stage PCa, as the lack of detailed knowledge surrounding disease progression has affected clinical management. The standard use of PSA as a biomarker has been of great benefit for monitoring disease outcomes in PCa (i.e. disease progression and response to therapy) . However, it offers no insight into the molecular mechanisms driving tumour progression or resistance to therapy. Thus, there is an unmet need to investigate PCa-specific temporal patterns that might uncover the foundations of disease progression and contribute towards the establishment of a biomarker-guided and risk stratification-based clinical management of PCa.
We employed an epigenome-wide array to understand longitudinal methylation patterns and how these might be influenced by disease progression and treatment change in men with metastatic PCa. As previously reported [8, 12], the use of array platforms with liquid biopsy samples averages DNA methylation from distinctive, and potentially multiple, cell types present. We therefore employed a methylation atlas to identify the source(s) of the main cellular components in each subject’s cfDNA . One-third of the cfDNA samples in our cohort had detectable prostate content; a frequency on par with previous mPCa cfDNA studies [16,17,18]. While statistical methods have been developed that adjust methylation data for cell-type heterogeneity , most focus on whole blood cell populations or more generic cell types (i.e. epithelial cells), overlooking the necessity to analyse tumour-specific patterns. Additionally, there is a dearth of statistical methods for n-of-1 or personalised approaches, such as this. Applying these approaches in this study could affect their robustness and lead to erroneous adjustments. Thus, the longitudinal observations made in this study are largely related to the total cfDNA content.
Temporal analysis of liquid biopsy DNA methylation patterns in mPCa individuals throughout disease progression, drug treatment(s) and relapse showed that most CpG sites did not substantially change. This finding supports previous research, which demonstrated that the human methylome is relatively stable over time [5, 20, 21]. However, the proportion of the methylome that demonstrated temporal variation was notably evident on administration of two PCa therapies: taxanes and AR inhibitors. To the best of our knowledge this is the first study to interrogate how therapies might affect the PCa methylome over time. Changes in CpG methylation following therapy administration were most pronounced in genes involved with immune response, including ones that encode essential components of the JAK-STAT pathway, an important regulator of the inflammatory response . This discovery is bolstered by previous reports, which showed that chemotherapy-induced methylation changes could influence immune cells and pathways [23, 24]. We also observed that this potential therapy-induced methylation shift was transient, and was only typically present during the time of drug administration. We suggest that this ephemeral pattern may be indicative of a systemic-epigenomic response to chemotherapy, rather than an epigenetic reprogramming, as previously observed during PCa progression , which would be unlikely to disappear after treatment cessation. Conversely to systemic treatments, AR inhibitors are a targeted PCa therapy that function by antagonising AR signalling . Although AR receptors are mostly expressed in reproductive tissues, they have also been implicated in the regulation of multiple cellular processes in a variety of other organs and systems [27, 28]. In fact, recent reports suggest that AR inhibitors might have a role in modulating the immune response in PCa [29, 30]. However, the molecular mechanisms by which AR inhibitors exert this immune-effect are still unknown. Our findings suggest that an epigenetic component could be involved.
Longitudinal assessment of methylation patterns allowed us to explore the dynamics of epigenetic events during metastatic PCa. Using our cohort as a discovery set and validating in two independent cohorts, we identified 2474 DMGs, with some (i.e. AR, GSTP1, RASSF1, APC) previously linked to PCa . Two noteworthy hypomethylated genes are EPN1 and CPEB4. EPN1 is known to actively promote tumorigenesis by enhancing cell surface receptor endocytosis and upregulating tumour growth-related pathways . CPEB4, an mRNA binding protein, actively reprogrammes gene expression by acting as a translational repressor/activator and promoting an invasive phenotype . While elevated expression of both genes has been reported in several invasive cancers, including prostate [31, 33], the mechanism of their upregulation was previously unknown. Here, we show that their promoter regions are hypomethylated in cfDNA (with detectable prostate content) from mPCa patients, suggesting that a loss of promoter methylation might play a role in the upregulation of CPEB4 and EPN1 in advanced PCa.
Individualised longitudinal studies, such as this, are important for precision medicine and to enable a better understanding of how to tailor clinical decisions [34, 35]. Indeed, by analysing serial plasma samples over the course of 36 months, Chen and colleagues were able to demonstrate that personal DNA methylomes showed distinct dynamic patterns that were associated with different physiological conditions, and that those patterns were present before the onset of the disease . Our study does have some limitations. As an observational study with a small cohort size, we cannot causally attribute methylation changes to specific therapies and any attempt to do so would be an over interpretation of the data. Further studies are needed to validate our findings and provide robust translational evidence for the application of this methodology to the clinic. We prioritised promoter methylation within the top 5% MVPs in each subject, in order to study the most temporally variable probes and focus on biologically relevant regions, as promoter methylation is highly correlated with gene expression. However, this represents only ~ 1% of CpG sites represented on the EPIC array and > 0.04% of all human CpGs. A deeper dive into how the methylome (beyond the conventional promoter CpG island) evolves during disease progression and therapeutic resistance is warranted.
Overall, we show the feasibility of studying longitudinal cfDNA methylation patterns in individuals and demonstrate how these profiles might be influenced by clinical events, such as therapeutic administration. Further analysis into the identification of PCa-specific temporal patterns might improve our understanding on how to better use cfDNA-methylation changes for improving monitoring and clinical management of this malignancy.
Patients were recruited between 2015 and 2016, to a Cancer Trials Ireland sponsored translational, multicentre, longitudinal study of men with mPCa titled CTRIAL-IE (ICORG) 14-04 Irish Programme for Stratified Prostate Cancer Therapy (iPROSPECT). All patients gave informed consent for serial blood sampling.
Additionally, two independent validation cohorts were used: (1) an FFPE tissue cohort of different histological stages of PCa (n = 44, Additional file 1: Table S1 and GSE157272), (2) and a previously described cfDNA dataset (n = 181 mPCa patients) .
Sample collection and DNA methylation profiling
Peripheral blood samples were collected at baseline (BL), and subsequent follow-ups (FU), every 4 ± 1 months. Plasma was isolated, by centrifugation at 2000×g for 15 min, within 2 h of collection and stored at – 80 °C. Forty-eight plasma samples (≤ 3 mL) from 9 participants across multiple time points were used for cfDNA isolation, as previously described . cfDNA was bisulfite modified using the Zymo EZ DNA methylation kit and quality control was by real-time PCR quantification, as previously described . Bisulfite modified DNA was subjected to the Illumina Restoration kit and analysed using the Illumina Infinium Human MethylationEPIC BeadChip platform (Fig. 5A). Raw data (GSE157273), were pre-processed using minfi . After QC, normalization and probe filtering (Additional file 1), a total of 801,849 probes were retained for analysis. A similar data processing pipeline was applied to the two validation cohorts, both run on the Illumina HM450K platform.
Longitudinal methylation analysis
The top 5% MVPs within each subject over time (calculated by SD of β-values), were used to study dynamic methylation patterns. Further analysis concentrated on 5′ regulatory probes, as characterised in Illumina’s EPIC array manifesto (TSS1500, TSS200 and 5′UTR), which were used to identify clusters of probes that had similar methylation patterns over time  (Fig. 5B and Additional file 1). Additionally, the effects of different therapies, taxanes and AR inhibitors, on the methylation of cfDNA were examined using the top 5% MVPs between commencement and cessation of therapy (Fig. 5C and Additional file 1). GO was performed as previously described [39, 40]. Finally, ctDNA dynamics were investigated by analysing genes identified in subjects with PC cfDNA (S008, S009, S019 and S029). Genes were considered differentially methylated (DMGs) after comparison with NPC samples from the whole dataset (β ≥ 0.1 difference, adjusted p value < 0.05) (Additional file 1). Validation of DMGs was carried out using two independent PCa cohorts (Fig. 5D).
Unpaired t-test/Mann–Whitney U test were used to identify DMGs to study ctDNA dynamics; Kruskal–Wallis and Dunn’s multiple comparison tests were used to evaluate methylation across different histological stages of PCa carcinogenesis. Analyses were performed using Prism 6 (GraphPad) and R (v3.6.3) and deemed significant if FDR-adjusted p value < 0.05. The Benjamini–Hochberg method was used for multiple testing correction when appropriate.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GEO repository (GSE157272, GSE157273).
Circulating tumour DNA
Differentially methylated genes
Luteinizing hormone-releasing hormone
Irish Programme for Stratified Prostate Cancer Therapy
Metastatic prostate cancer
Most variable probes
Prostate specific antigen
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
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(3):616–25.
Bhasin JM, Lee BH, Matkin L, Taylor MG, Hu B, Xu Y, et al. Methylome-wide sequencing detects DNA hypermethylation distinguishing indolent from aggressive prostate cancer. Cell Rep. 2015;13(10):2135–46.
Hendriks RJ, Dijkstra S, Smit FP, Vandersmissen J, Van de Voorde H, Mulders PFA, et al. Epigenetic markers in circulating cell-free DNA as prognostic markers for survival of castration-resistant prostate cancer patients. Prostate. 2018;78(5):336–42.
Feinberg AP, Irizarry RA, Fradin D, Aryee MJ, Murakami P, Aspelund T, et al. Personalized epigenomic signatures that are stable over time and covary with body mass index. Sci Transl Med. 2010;2(49):49ra67.
Wang Y, Karlsson R, Jylhava J, Hedman AK, Almqvist C, Karlsson IK, et al. Comprehensive longitudinal study of epigenetic mutations in aging. Clin Epigenet. 2019;11(1):187.
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 Invest. 2020;130(4):1991–2000.
Silva R, Moran B, Russell NM, Fahey C, Vlajnic T, Manecksha RP, et al. Evaluating liquid biopsies for methylomic profiling of prostate cancer. Epigenetics. 2020;15(6–7):715–727.
Liu MC, Oxnard GR, Klein EA, Swanton C, Seiden MV, Liu MC, et al. Sensitive and specific multi-cancer detection and localization using methylation signatures in cell-free DNA. Ann Oncol. 2020;31(6):745–59.
Wyatt AW, Annala M, Aggarwal R, Beja K, Feng F, Youngren J, et al. Concordance of circulating tumor DNA and matched metastatic tissue biopsy in prostate cancer. J Natl Cancer Inst. 2017;109(12). https://doi.org/10.1093/jnci/djx118.
Bardelli A, Pantel K. Liquid biopsies, what we do not know (yet). Cancer Cell. 2017;31(2):172–9.
Moss J, Magenheim J, Neiman D, Zemmour H, Loyfer N, Korach A, et al. Comprehensive human cell-type methylation atlas reveals origins of circulating cell-free DNA in health and disease. Nat Commun. 2018;9(1):5068.
Gordevicius J, Krisciunas A, Groot DE, Yip SM, Susic M, Kwan A, et al. Cell-free DNA modification dynamics in abiraterone acetate-treated prostate cancer patients. Clin Cancer Res. 2018;24(14):3317–24.
Massie CE, Mills IG, Lynch AG. The importance of DNA methylation in prostate cancer development. J Steroid Biochem Mol Biol. 2017;166:1–15.
Lilja H, Ulmert D, Vickers AJ. Prostate-specific antigen and prostate cancer: prediction, detection and monitoring. Nat Rev Cancer. 2008;8(4):268–78.
Bettegowda C, Sausen M, Leary RJ, Kinde I, Wang Y, Agrawal N, et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci Transl Med. 2014;6(224):224ra24.
Vandekerkhove G, Struss WJ, Annala M, Kallio HML, Khalaf D, Warner EW, et al. Circulating tumor DNA abundance and potential utility in de novo metastatic prostate cancer. Eur Urol. 2019;75(4):667–75.
Beltran H, Romanel A, Conteduca V, Casiraghi N, Sigouros M, Franceschini GM, et al. Circulating tumor DNA profile recognizes transformation to castration-resistant neuroendocrine prostate cancer. J Clin Invest. 2020;130(4):1653–68.
Brägelmann J, Lorenzo BJ. A comparative analysis of cell-type adjustment methods for epigenome-wide association studies based on simulated and real data sets. Brief Bioinform. 2019;20(6):2055–65.
Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, Cui H, et al. Intra-individual change over time in DNA methylation with familial clustering. JAMA. 2008;299(24):2877–83.
Byun HM, Nordio F, Coull BA, Tarantini L, Hou L, Bonzini M, et al. Temporal stability of epigenetic markers: sequence characteristics and predictors of short-term DNA methylation variations. PLoS ONE. 2012;7(6):e39220.
Thomas SJ, Snowden JA, Zeidler MP, Danson SJ. The role of JAK/STAT signalling in the pathogenesis, prognosis and treatment of solid tumours. Br J Cancer. 2015;113(3):365–71.
Smith AK, Conneely KN, Pace TW, Mister D, Felger JC, Kilaru V, et al. Epigenetic changes associated with inflammation in breast cancer patients treated with chemotherapy. Brain Behav Immun. 2014;38:227–36.
Yao S, Hu Q, Kerns S, Yan L, Onitilo AA, Misleh J, et al. Impact of chemotherapy for breast cancer on leukocyte DNA methylation landscape and cognitive function: a prospective study. Clin Epigenet. 2019;11(1):45.
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(8):790–9.
Schalken J, Fitzpatrick JM. Enzalutamide: targeting the androgen signalling pathway in metastatic castration-resistant prostate cancer. BJU Int. 2016;117(2):215–25.
Rana K, Davey RA, Zajac JD. Human androgen deficiency: insights gained from androgen receptor knockout mouse models. Asian J Androl. 2014;16(2):169–77.
Chang C, Lee SO, Yeh S, Chang TM. Androgen receptor (AR) differential roles in hormone-related tumors including prostate, bladder, kidney, lung, breast and liver. Oncogene. 2014;33(25):3225–34.
Madan RA, Donahue RN, Singh H, Karzai F, Theoret MR, Heery CR, et al. Clinical and immunologic impact of short course enzalutamide without androgen deprivation therapy for biochemically recurrent prostate cancer. J Clin Oncol. 2016;34(2):214.
Graff JN, Beer TM, Alumkal JJ, Slottke RE, Redmond WL, Thomas GV, et al. A phase II single-arm study of pembrolizumab with enzalutamide in men with metastatic castration-resistant prostate cancer progressing on enzalutamide alone. J Immunother Cancer. 2020;e000642. https://doi.org/10.1136/jitc-2020-000642.
Tessneer KL, Pasula S, Cai X, Dong Y, Liu X, Yu L, et al. Endocytic adaptor protein epsin is elevated in prostate cancer and required for cancer progression. ISRN Oncol. 2013;2013:420597.
D’Ambrogio A, Nagaoka K, Richter JD. Translational control of cell growth and malignancy by the CPEBs. Nat Rev Cancer. 2013;13(4):283–90.
Xu H, Liu B. CPEB4 is a candidate biomarker for defining metastatic cancers and directing personalized therapies. Med Hypotheses. 2013;81(5):875–7.
Chen R, Xia L, Tu K, Duan M, Kukurba K, Li-Pook-Than J, et al. Longitudinal personal DNA methylome dynamics in a human with a chronic condition. Nat Med. 2018;24(12):1930–9.
Chen R, Mias GI, Li-Pook-Than J, Jiang L, Lam HY, Chen R, et al. Personal omics profiling reveals dynamic molecular and medical phenotypes. Cell. 2012;148(6):1293–307.
Campan M, Weisenberger DJ, Trinh B, Laird PW. MethyLight. Methods Mol Biol. 2009;507:325–37.
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(10):1363–9.
Mori U, Mendiburu A, Lozano J. TSdist: Distance Measures for Time Series Data. 2019; R package version 3.6.
Geeleher P, Hartnett L, Egan LJ, Golden A, Raja Ali RA, Seoighe C. Gene-set analysis is severely biased when applied to genome-wide methylation data. Bioinformatics. 2013;29(15):1851–7.
Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.
We gratefully acknowledge the support of Cancer Trials Ireland, in particular Dr. Orla Casey and Dr. Verena Murphy, and Clinical Research Development Ireland (CRDI), in particular Dr Suzanne Bracken for her work in setting up this translational study.
This work was supported by grant funding from the Irish Cancer Society and Movember (PCT13MCD and CRS17SIL), as well as Science Foundation Ireland under the Investigator Programme OPTi-PREDICT (15/IA/3104) and the Strategic Research Programme Precision Oncology Ireland (18/SPP/3522).
Ethics approval and consent to participate
Patients were recruited through an institutional review board-approved, translational, multicentre, longitudinal study. Consent to participate was collected for every patient in the study.
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.
About this article
Cite this article
Silva, R., Moran, B., Baird, AM. et al. Longitudinal analysis of individual cfDNA methylome patterns in metastatic prostate cancer. Clin Epigenet 13, 168 (2021). https://doi.org/10.1186/s13148-021-01155-w
- Metastatic prostate cancer
- DNA methylation
- Longitudinal study
- Liquid biopsy